share/hoomd/examples/dump_dcd

   1 hoomd_script::init::create_random hoomd_script::dump::dcd hoomd_script::pair::lj
   2 # Due to deficiencies in doxygen, the commands used in this example are listed explicitly here
   3 # run this script with "python -x filename" to skip the first line, or remove this header
   4
   5 # ---- dump_dcd.py ----
   6 from hoomd_script import *
   7
   8 # create 1000 random particles of name A
   9 init.create_random(N=1000, phi_p=0.01, name='A')
  10
  11 # specify Lennard-Jones interactions between particle pairs
  12 lj = pair.lj(r_cut=2.5)
  13 lj.pair_coeff.set('A', 'A', epsilon=1.0, sigma=1.0)
  14
  15 # integrate at constant temperature
  16 all = group.all()
  17 integrate.mode_standard(dt=0.005)
  18 integrate.nvt(group=all, T=1.2, tau=0.5)
  19
  20 # dump an xmle file for the structure information
  21 xml = dump.xml(filename='dump_dcd.xml', vis=True)
  22
  23 # dump a .dcd file for the trajectory
  24 dump.dcd(filename='dump_dcd.dcd', period=100)
  25
  26 # run 10,000 time steps
  27 run(10e3)
  28