csipair.xmo

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

XMD (Version 2.5.29 Jan 20 2002)

Not using pthread library.
Not using asm("finit"); patch in thread routines.
#######################################################################
#                                                                     #
#  Example:  Using auxillary pair potentials with Tersoff's C-Si      #
#                                                                     #
#######################################################################
#                                                                     #
#    This file gives an example of how use auxillary pair potentials  #
#  with Tersoff's C-Si potential.  The energies and forces due to     #
#  these pair potentials are summed with the energy and forces due    #
#  to the Tersoff C-Si potential.  Pair potentials may be specified   #
#  between one or more pairs atom types - each such pair of atoms     #
#  will require its own potential.                                    #
#                                                                     #
#    The pair potential is entered as a polynomial of the form        #
#                                                                     #
#                       2         3                n+1                #
#          f(r) =   a  u   +  a  u   + ... + a    u                   #
#                    2         3              n+1                     #
#                                                                     #
#  where u = rc-r.  This polyomial form provides a simple and         #
#  compact specification while insuring that the potential            #
#  has a smooth cutoff of the r=rc.                                   #
#                                                                     #
#    The syntax of the pair potential command is                      #
#                                                                     #
#      POTENTIAL PAIR t1 t2 n rc a2 [a3 .. an+1]                      #
#                                                                     #
#    where                                                            #
#                                                                     #
#      t1,t2  are the atom types.                                     #
#      n      is the number of polymial terms.                        #
#      rc     is the pair potential cutoff in units of angtroms.      #
#      a2     is the first term, it is the coefficient of u^2.        #
#             It's units are (eunit/cm^2) where eunit is specified    #
#             by the EUNIT command.  The default value for eunit      #
#             is Kelvin.                                              #
#      a3     is the optional second term, the coefficient of u^3     #
#             It's units are (eunit/cm^3).                            #
#      an+1   is the optional n'th term, the coefficient of u^(n+1)   #
#             It's units are (eunit/cm^(n+1)).                        #
#                                                                     #
#######################################################################

#  Switch on Tersoff C-Si potential
potential set tersoff

   Tersoff's Silicon-Carbon potential
   reference:
   "Modeling solid-state chemistry:  Interatomic potentials
   for multicomponent systems",
   J. Tersoff, Phys Rev. B 39 (8), 5566-5568 (15 March 1989)


#  Set up some useful constants for this potential
calc  C=1
calc  Si=2
calc  Al=3
calc  MassC=12.01
calc  MassSi=28.086
calc  MassAl=26.982
calc  A0=4.32
calc  DTIME=0.8e-15

#  Add pair potentials
calc NCOEFF=1
calc CUTOFF=3.0
calc COEFF=5e18
potential pair C  Al NCOEFF CUTOFF COEFF
potential pair Si Al NCOEFF CUTOFF COEFF
potential pair Al Al NCOEFF CUTOFF COEFF



#  Form diamond unit cell
box 2 2 2
particle  8
dup  1  1 0 0
dup  1  0 1 0
dup  1  0 0 1

#  Scale to the unit cell length of diamond CSi
scale 4.32

#  Set particle masses
select type C
***  NUMBER SELECTED 24
mass MassC
select type Si
***  NUMBER SELECTED 24
mass MassSi
select type Al
***  NUMBER SELECTED 16
mass MassAl

#  Set simulation step size
dtime DTIME

#  Instruct XMD to save energy every CMD step in file csi.e
esave 1 csipair.e

#  Initialize velocities to 1000K
itemp 1000

#  Set clamp off (adiabatic simulation)
clamp off

# Equilibrate particles for 500 step
write cor +csipair.cor
repeat 50
   cmd 10
   write cor +csipair.cor
end
***   Current step is 10
***   Current step is 20
***   Current step is 30
***   Current step is 40
***   Current step is 50
***   Current step is 60
***   Current step is 70
***   Current step is 80
***   Current step is 90
***   Current step is 100
***   Current step is 110
***   Current step is 120
***   Current step is 130
***   Current step is 140
***   Current step is 150
***   Current step is 160
***   Current step is 170
***   Current step is 180
***   Current step is 190
***   Current step is 200
***   Current step is 210
***   Current step is 220
***   Current step is 230
***   Current step is 240
***   Current step is 250
***   Current step is 260
***   Current step is 270
***   Current step is 280
***   Current step is 290
***   Current step is 300
***   Current step is 310
***   Current step is 320
***   Current step is 330
***   Current step is 340
***   Current step is 350
***   Current step is 360
***   Current step is 370
***   Current step is 380
***   Current step is 390
***   Current step is 400
***   Current step is 410
***   Current step is 420
***   Current step is 430
***   Current step is 440
***   Current step is 450
***   Current step is 460
***   Current step is 470
***   Current step is 480
***   Current step is 490
***   Current step is 500

TIME INFORMATION
   Ending   Time: 01/20/2002   11:57:32
   Elapsed  Time: 0 min 02 sec

DYNAMICS STATISTICS
   Number of MD Steps:                 500
   Number of Neighbor Searchs:         103
   Time spent on MD steps:               2 (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