Advertisement
Guest User

input

a guest
Jul 23rd, 2019
214
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.90 KB | None | 0 0
  1. variable lattconst equal 5
  2. variable repetition equal 10
  3. log Al_temps4_10.log append
  4.  
  5. # ---------- Initialize Simulation ---------------------
  6. clear
  7. units metal # pour avoir des unités en Angstroms et eV, entre autres
  8. dimension 3
  9. boundary p p p # pour définir que les limites dans les directions de chaque dimension sont périodiques
  10. atom_style atomic
  11. atom_modify map array
  12.  
  13. # ---------- Create Atoms ---------------------
  14. lattice fcc ${lattconst} # type de maille suivi du lattice constant
  15. region box block 0 1 0 1 0 1 units lattice
  16. create_box 1 box
  17.  
  18. lattice fcc ${lattconst} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 # un latice unit par direction
  19. create_atoms 1 box
  20. replicate ${repetition} ${repetition} ${repetition} # réplique chaque cellule périodique dans toutes directions
  21.  
  22. # ---------- Define Interatomic Potential ---------------------
  23. pair_style eam/alloy # détermine quel type de potentiel interatomique est utilisé
  24. pair_coeff * * Al99.eam.alloy Al # détermine où se trouvent les coefficients du type de potentiel, l'extension peut servir de pointeur vers celui à utiliser.
  25. neighbor 2.0 bin
  26. neigh_modify delay 0 every 1 check yes
  27.  
  28. # ---------- Run Minimization ---------------------
  29.  
  30. timestep 0.001
  31. #fix 1 all nvt temp 1000.0 1000.0 0.1
  32. fix 2 all box/relax iso 0.0
  33.  
  34.  
  35. thermo 100
  36. thermo_style custom step temp enthalpy press cella cellb cellc
  37. thermo_modify flush yes
  38. thermo_modify norm yes
  39.  
  40. min_style cg
  41. minimize 1e-25 1e-25 5000 10000
  42.  
  43. fix 1 all npt temp 1000.0 1000.0 0.1 iso 1000.0 1000.0 0.1
  44.  
  45. variable iterations equal 10000
  46. #dump mydump all xyz 1000 C:\Bin\LAMMPS\Exo7\dump5.xyz
  47. run ${iterations}
  48.  
  49.  
  50.  
  51. # définition de variables
  52. variable natoms equal "count(all)"
  53. variable latticeconstant equal "cella/10"
  54.  
  55. # affichage des variables définies
  56. print "Number of atoms = ${natoms};"
  57. print "Lattice constant = ${latticeconstant}"
  58.  
  59. print "All done!"
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement