relax.xmo

一个很好的分子动力学程序

XMD (Version 2.5.29 Nov 27 2000)

Not using pthread library.
Not using asm("finit"); patch in thread routines.
#######################################################################
#                                                                     #
#  Example:  Relaxing a lattice using CMD and QUENCH commands         #
#                                                                     #
#######################################################################
#                                                                     #
#    This example shows a method for relaxing a lattice to a low      #
#  energy state.  This is done by first performing a molecular        #
#  dynamics run at a moderate or low temperature (in this case 200K)  #
#  followed by a QUENCH command.  The QUENCH command acts like a      #
#  normal adiabatic (constant energy, with CLAMP set to -1)           #
#  molecular dynamics run only while the potential energy drops.      #
#  When the potential energy starts to rise, all the atom velocities  #
#  are automatically set to zero.  This removes the momentum that     #
#  would have carried the lattice uphill in potential energy,         #
#  allowing the lattice to beginning falling downhill again.          #
#                                                                     #
#######################################################################

# Relax lattice

# Read potential for nial
read ../nial.txt

# Make repeating box and lattice (in units of a0)
box 6 6 6
particle 2
dup 5   1 0 0
dup 5   0 1 0
dup 5   0 0 1

# Scale up to units of angstroms (2.8712 unit cell)
scale 2.8712

# Save energies from every dynamics step in file "timestep.e'
esave 10 relax.e

# Set particle masses (in atomic mass units)
select type 1
***  NUMBER SELECTED 216
mass 58.71
select type 2
***  NUMBER SELECTED 216
mass 26.982

# Set adiabtic simulation at starting temperature of 200K
clamp -1
itemp 200

#  This may not be the optimum step size
dtime  1e-15

#  Equilibrate at 200K for 500 steps
write rcv relax.rcv
repeat 5
   cmd 100
   write rcv +relax.rcv
end
***   Current step is 100
***   Current step is 200
***   Current step is 300
***   Current step is 400
***   Current step is 500

#  Quench for 500 steps
repeat 5
   quench 100
   write rcv +relax.rcv
end
***   Current step is 600
***   Current step is 700
***   Current step is 800
***   Current step is 900
***   Current step is 1000

TIME INFORMATION
   Ending   Time: 03/19/2001   13:41:03
   Elapsed  Time: 0 min 30 sec

DYNAMICS STATISTICS
   Number of MD Steps:                1000
   Number of Neighbor Searchs:           1
   Time spent on MD steps:              29 (secs)
   Time spent on Neighbor Search:        0 (secs)

ERROR STATISTICS
   Number of Fatal Errors:               0
   Number of Unknown Command Errors:     0
   Number of Misc. Warnings:             0