Modar Tutorial
Example 3: Ionic Liquid Simulation
     step7.2_merge.inp

 

#load force field

include "loadff.inc"

 

#load solute with contour ions

loadmsf file="step6_equil_npt.msf"

 

#load water box by appending

loadmsf file="step7_waterbox.msf" append=true

 

#load PBC crystal setup

include "pbc.inc"

 

#remove overlaping waters

defcls name=solute select="heavyatom && ! resn TIP*"

delete atom select=".byres. ( around solute cut 2.8 && (resn TIP* && type O*) )"

delete atom select=".byres. ( resn tip* && z >-22 && z <22 )"

 

#resort sequence

atom resortseq  start=1 select="resn tip*"

#tell center select=all

#tell size   select=all

 

tell msf_info

echo -hs "| "

echo -hs "|  Summary of MD system:"

echo -hs "|    number of segments:    $nseg"

echo -hs "|    number of residues:    $nres"

echo -hs "|    number of atoms:       $natom"

echo -hs "|    number of bonds:       $nbond"

echo -hs "|    number of angles:      $nangl"

echo -hs "|    number of dihedrals:   $ndihe"

echo -hs "|    number of impropers:   $nimpr"

echo -hs "|    number of cross terms: $ncmap"

echo -hs "|    net charge:  $netcharge"

echo -hs "|    total mass:  $totmass"

 

tell boxsize

 

echo -hs "|    PBC box: $boxtype, $boxa x $boxb x $boxc"

echo -hs "| "

 

#checkmsf

checkmsf

 

#remove fixes if have

atom fix select="none"

 

savemsf file="step7_two_phases.msf"

savepdb fmt=pdb file="step7_two_phases.pdb"

 

echo -hs "|  output files:"

echo -hs "|    step7_two_phases.msf"

echo -hs "|    step7_two_phases.pdb"

echo -hs "| "

  • Please read the comments first.

 

 

 

  Contact us

  Phone: 400-660-8656
  Email: support@beemd.org

 

       我们长期和北京市计算中心合作提供计算培训服务,承接托管计算业务,如有需求请随时联系我们。