aiskillstore/packmol

Packmol Skill

Build initial configurations for molecular dynamics simulations using Packmol.

What is Packmol?

Packmol creates initial configurations for MD simulations by packing molecules according to spatial constraints. It places molecules in boxes, around proteins, at interfaces, or within complex geometries (spheres, cylinders, ellipsoids) while ensuring no overlaps.

Installation

Install Packmol via pip:

pip install packmol

Verify installation:

packmol -h

For more installation options, see the Packmol website.

Quick Start

Basic Box Packing

Create a simple box of water molecules:

# water_box.inp
tolerance 2.0
filetype pdb
output water_box.pdb

structure water.pdb
  number 1000
  inside box 0. 0. 0. 40. 40. 40.
end structure

Run Packmol:

packmol < water_box.inp

Solvate a Protein

Solvate a protein with water and ions:

# solvation.inp
tolerance 2.0
filetype pdb
output solvated.pdb

structure protein.pdb
  number 1
  fixed 0. 0. 0. 0. 0. 0.
  center
end structure

structure water.pdb
  number 5000
  inside box -10. -10. -10. 50. 50. 50.
end structure

structure SOD.pdb
  number 10
  inside box -10. -10. -10. 50. 50. 50.
end structure

structure CLA.pdb
  number 10
  inside box -10. -10. -10. 50. 50. 50.
end structure

Liquid-Liquid Interface

Build a water/chloroform interface:

# interface.inp
tolerance 2.0
filetype pdb
output interface.pdb
pbc -20. -20. -30. 20. 20. 30.

structure water.pdb
  number 1000
  below plane 0. 0. 1. 0.
end structure

structure chloroform.pdb
  number 200
  above plane 0. 0. 1. 0.
end structure

Core Concepts

Input File Structure

Every Packmol input file requires:

1. tolerance: Minimum distance between atoms (Å)
2. output: Output filename
3. filetype: Format (pdb, xyz, tinker)
4. structure blocks: Define molecules to place

Structure Block Syntax

structure molecule.pdb
  number <N>                          # Number of molecules
  inside|outside <constraint>         # Spatial constraint
  [optional parameters]
end structure

Common Constraint Types

• box: inside box xmin ymin zmin xmax ymax zmax
• sphere: inside sphere xcenter ycenter zcenter radius
• cylinder: inside cylinder x1 y1 z1 dx dy dz radius length
• plane: above plane a b c d or below plane a b c d
• ellipsoid: inside ellipsoid xc yc zc xa yb zc scale

See references/constraints.md for complete constraint documentation.

Workflows

1. Basic Molecular Packing

Build boxes with multiple molecule types.

Example: Water/ethanol mixture

tolerance 2.0
output mixture.pdb
filetype pdb

structure water.pdb
  number 800
  inside box 0. 0. 0. 40. 40. 40.
end structure

structure ethanol.pdb
  number 200
  inside box 0. 0. 0. 40. 40. 40.
end structure

2. Protein Solvation

Solvate biomolecules with water and ions for neutralization.

Key parameters:

• Use fixed with center for the protein
• Add Na+/Cl- ions for neutrality and concentration
• Calculate box size based on protein + solvent shell

Automatic solvation helper:

python scripts/solvate_helper.py protein.pdb --shell 15.0 --charge +4

3. Interface Systems

Build liquid-liquid or liquid-vapor interfaces using plane constraints.

Example: Water/hexane interface

tolerance 2.0
output interface.pdb
pbc -20. -20. -30. 20. 20. 30.

structure water.pdb
  number 1000
  below plane 0. 0. 1. 0.
end structure

structure hexane.pdb
  number 200
  above plane 0. 0. 1. 0.
end structure

4. Advanced Constraints

Use spherical, cylindrical, or ellipsoidal constraints for complex geometries.

Example: Spherical vesicle

structure lipid.pdb
  number 2000
  inside sphere 0. 0. 0. 40.
  atoms 1 2 3 4
    outside sphere 0. 0. 0. 35.
  end atoms
end structure

structure water.pdb
  number 2000
  inside sphere 0. 0. 0. 35.
end structure

structure water.pdb
  number 5000
  outside sphere 0. 0. 0. 45.
end structure

Input Parameters

Required Parameters

• tolerance <distance>: Minimum intermolecular distance (Å). Default: 2.0 for all-atom
• output <filename>: Output file name
• filetype <format>: pdb, xyz, or tinker

Optional Parameters

• pbc <dimensions>: Periodic boundary conditions (e.g., pbc 30. 30. 60.)
• seed <integer>: Random seed for reproducibility
• discale <factor>: Distance scaling for optimization (default: 1.0)
• maxit <N>: Maximum iterations (default: 20)
• precision <value>: Convergence precision (default: 0.01)

See references/parameters.md for complete parameter reference.

Structure Block Options

Positioning Options

• number: Molecule count
• inside/outside: Spatial constraint
• fixed: Fix position and rotation (6 parameters: x, y, z, α, β, γ)
• center: Use center of mass for positioning

Rotation Constraints

constrain_rotation x 180. 20.  # Constrain rotation around x-axis
constrain_rotation y 180. 20.  # Constrain rotation around y-axis
constrain_rotation z 180. 20.  # Constrain rotation around z-axis

Atom Selection

Apply constraints to specific atoms within molecules:

structure molecule.pdb
  number 100
  inside box 0. 0. 0. 30. 30. 30.
  atoms 1 2 3
    inside box 0. 0. 25. 30. 30. 30.
  end atoms
end structure

Running Packmol

Basic Execution

packmol < input.inp

Output Interpretation

Success message:

------------------------------
Success!
Final objective function value: .22503E-01
Maximum violation of target distance: 0.000000
Maximum violation of the constraints: .78985E-02
------------------------------

Check that both violations are < 0.01 for a valid solution.

Validation

Check Overlaps

python scripts/check_overlaps.py output.pdb --tolerance 2.0

Verify Success

python scripts/verify_success.py input.inp output.pdb

Analyze Density

python scripts/analyze_density.py output.pdb

Validate Input

python scripts/validate_input.py input.inp

Troubleshooting

Common Issues

1. "Killed" error: System too large

   • Reduce number of molecules
   • Use restart files to build incrementally
   • See references/troubleshooting.md

2. No convergence:

   • Try discale 1.5 to scale distances
   • Reduce molecule count
   • Simplify constraints
   • Increase maxit

3. Strange geometries:

   • Add check keyword to validate constraints without packing
   • Verify constraint syntax
   • Check for conflicting constraints

4. Incorrect atom count:

   • Verify structure files are readable
   • Check for duplicate atoms in input files
   • Validate with scripts/validate_input.py

See references/troubleshooting.md for detailed solutions.

Examples

Explore example input files in the examples/ directory:

• Basic: examples/basic/ - Simple boxes and mixtures
• Solvation: examples/solvation/ - Proteins with water and ions
• Interface: examples/interface/ - Liquid-liquid interfaces
• Advanced: examples/advanced/ - Vesicles, bilayers, complex geometries

Templates

Use templates in templates/ as starting points:

• templates/basic_template.inp: Minimal template for simple packing
• templates/solvation_template.inp: Protein solvation setup
• templates/interface_template.inp: Interface systems

Helper Scripts

Use Python scripts in scripts/ for automation:

• generate_input.py: Generate inputs programmatically
• validate_input.py: Validate input syntax before running
• check_overlaps.py: Detect atomic overlaps in output
• analyze_density.py: Calculate system density
• solvate_helper.py: Automatic protein solvation setup
• verify_success.py: Verify Packmol completed successfully

Advanced Topics

Periodic Boundary Conditions

Use pbc for periodic systems:

pbc 30. 30. 60.  # or pbc xmin ymin zmin xmax ymax zmax

Restart Files

Build large systems incrementally:

structure water.pdb
  number 1000
  inside box 0. 0. 0. 40. 40. 40.
  restart_to water1.pack
end structure

Then restart:

structure water.pdb
  number 1000
  restart_from water1.pack
end structure

Atom-Specific Radii

Set different radii for multiscale models:

structure molecule.pdb
  number 100
  radius 1.5  # All atoms
end structure

structure molecule.pdb
  number 100
  atoms 1 2
    radius 1.5  # Specific atoms
  end atoms
end structure

Constraint Validation

Validate constraints without packing:

structure molecule.pdb
  number 100
  inside box 0. 0. 0. 30. 30. 30.
  check
end structure

Best Practices

1. Start simple: Test with few molecules before scaling up
2. Use appropriate tolerance: 2.0 Å for all-atom, larger for coarse-grained
3. Check constraints: Add check keyword to validate regions
4. Validate output: Use scripts to check overlaps and density
5. Reproducibility: Set seed for repeatable results
6. Large systems: Use restart files or build in stages
7. Box size: Allow 10-15 Å padding around solutes for solvation

Tips for Common Use Cases

Protein Solvation

• Add 10-15 Å solvent shell around protein
• Calculate ions for neutrality: N_ions = charge / e
• Add salt ions for desired concentration (e.g., 0.15 M NaCl)
• Use fixed with center for protein positioning

Mixed Solvents

• Calculate total number of molecules from desired molar ratios
• Use same tolerance for all components
• Test with small systems first

Membrane Systems

• Use constrain_rotation to orient lipids
• Build in stages: lipids first, then water
• Consider using specialized membrane builders for large systems

Nanotubes/Pores

• Use cylinder constraint for pore region
• Combine with outside constraint for bulk region
• May need atom selection for specific molecule orientations

Resources

• Official documentation: Packmol User Guide
• Examples: Packmol Examples
• GitHub: Packmol Repository
• Paper: Martínez et al. J Comput Chem 2009

References

For detailed information on specific topics, see:

• constraints.md - Complete constraint syntax and examples
• parameters.md - All input parameters and options
• file_formats.md - File format specifications
• troubleshooting.md - Problem-solving guide