Skip to content

Instantly share code, notes, and snippets.

@martok
Created November 6, 2017 14:22
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save martok/9102b3e174e960e2f8a0be1fd8c7e3c4 to your computer and use it in GitHub Desktop.
Save martok/9102b3e174e960e2f8a0be1fd8c7e3c4 to your computer and use it in GitHub Desktop.
units metal
dimension 3
boundary p p p
atom_style atomic
timestep 0.001
variable min2th index 20
variable max2th index 90
lattice bcc 3.147
region box block 0 1 0 1 0 1 units lattice
create_box 1 box
create_atoms 1 region box
mass 1 95.95
pair_style none
atom_modify sort 0 0.0
compute XRD all xrd 1.541838 "Mo" 2Theta ${min2th} ${max2th} echo LP 1 c 1 1 1
fix 1 all ave/histo/weight 1 1 1 ${min2th} ${max2th} $((v_max2th-v_min2th)/0.1) c_XRD[1] c_XRD[2] mode vector file hist.xrd
run 0
unfix 1
uncompute XRD
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment