screw.xmo

来自「一个很好的分子动力学程序」· XMO 代码 · 共 74 行

XMD (Version 2.5.29 Nov 27 2000)

Not using pthread library.
Not using asm("finit"); patch in thread routines.
#######################################################################
#                                                                     #
#  Example:  Making a screw dislocation with the SCREW command        #
#                                                                     #
#######################################################################
#                                                                     #
#    The SCREW command displaces the *selected* atoms to form a       #
#  screw dislocation in a previously perfect lattice.                 #
#                                                                     #
#######################################################################
#                                                                     #
# The format of the SCREW command is                                  #
#                                                                     #
#   SCREW xburgers yburgers zburgers xorg yorg zorg [xref yref zref]  #
#                                                                     #
#  where   xburgers    the components of the burgers vector           #
#          yburgers                                                   #
#          zburgers                                                   #
#          xorg        a point through which the screw axis passes    #
#          yorg                                                       #
#          zorg                                                       #
#          xref        optional, a vector non-parallel with the       #
#          yref         burgers vector.  Those atoms which lie in     #
#          zref         plane containing this and the burgers vector  #
#                       will have a zero displacement.  If this       #
#                       option is not specified, then the plane of    #
#                       zero displacement will be defined by an       #
#                       arbitrary vector                              #
#                                                                     #
#######################################################################
#
#
#  Allocate storage
allocate 2000 40000
WARNING:  Allocate command is now obsolete.
#
#
#  Create an FCC lattice with unit cell length of 1 angstrom
#
box 4 4 4
particle 4
dup  3  1 0 0
dup  3  0 1 0
dup  3  0 0 1
#
#
#  Write coordinates
#
step 0
write cor screw.cor
#
#
#  Insert screw dislocation
#
screw  0 0 1  2.25 2.25 2.25  1 0 0
#
#
#  Write coordinates
step 1
write cor +screw.cor

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

ERROR STATISTICS
   Number of Fatal Errors:               0
   Number of Unknown Command Errors:     0
   Number of Misc. Warnings:             1   <=======  !!!