input

variable lattconst equal 5
variable repetition equal 10
log Al_temps4_10.log append

# ---------- Initialize Simulation ---------------------
clear
units metal # pour avoir des unités en Angstroms et eV, entre autres
dimension 3
boundary p p p # pour définir que les limites dans les directions de chaque dimension sont périodiques
atom_style atomic
atom_modify map array

# ---------- Create Atoms ---------------------
lattice 	fcc ${lattconst} # type de maille suivi du lattice constant
region	box block 0 1 0 1 0 1 units lattice
create_box	1 box

lattice	fcc ${lattconst} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1  # un latice unit par direction
create_atoms 1 box
replicate ${repetition} ${repetition} ${repetition}  # réplique chaque cellule périodique dans toutes directions

# ---------- Define Interatomic Potential ---------------------
pair_style eam/alloy # détermine quel type de potentiel interatomique est utilisé
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.
neighbor 2.0 bin
neigh_modify delay 0 every 1 check yes

# ---------- Run Minimization ---------------------

timestep 0.001
#fix 1 all nvt temp 1000.0 1000.0 0.1
fix 2 all box/relax iso 0.0


thermo 100
thermo_style custom step temp enthalpy press cella cellb cellc
thermo_modify flush yes
thermo_modify norm yes

min_style cg
minimize 1e-25 1e-25 5000 10000

fix 1 all npt temp 1000.0 1000.0 0.1 iso 1000.0 1000.0 0.1

variable iterations equal 10000
#dump mydump all xyz 1000 C:\Bin\LAMMPS\Exo7\dump5.xyz
run ${iterations}



# définition de variables
variable natoms equal "count(all)"
variable latticeconstant equal "cella/10"

# affichage des variables définies
print "Number of atoms = ${natoms};"
print "Lattice constant = ${latticeconstant}"

print "All done!"