units metal # Angstrom, eV, ps, bar boundary p p p # Περιοδικό σύστημα. atom_style atomic # Τύπος σωματιδίων στην προσομοίωση. #----------------------------------------------------------------------- # Κρυσταλλικό σύστημα και πλεγματική σταθερά. #----------------------------------------------------------------------- lattice fcc 3.63 #----------------------------------------------------------------------- # Επανάληψη στις τρείς διευθύνσεις και τοποθέτηση ατόμων. #----------------------------------------------------------------------- region my_box block 0 5 0 5 0 5 create_box 1 my_box create_atoms 1 box displace_atoms all move 0.9 0.9 0.9 units box #----------------------------------------------------------------------- # Ορισμός δυναμικού αλληλεπίδρασης. #----------------------------------------------------------------------- pair_style eam pair_coeff * * Cu.eam #----------------------------------------------------------------------- # Ποιές θερμοδυναμικές πληροφορίες θα εμφανιστούν στην οθόνη. #----------------------------------------------------------------------- thermo_style custom step temp pe press lx ly lz thermo 1 #----------------------------------------------------------------------- # Ορισμός στατιστικού συνόλου (ισόθερμο-ισοβαρές NPT). #----------------------------------------------------------------------- fix my_npt all npt temp 300 300 0.01 iso 0 0 1 #----------------------------------------------------------------------- # Ορισμός αρχικών ταχυτήτων. #----------------------------------------------------------------------- velocity all create 600 4928459 mom yes rot yes dist gaussian #----------------------------------------------------------------------- # Χρονικό βήμα (σε ps) #----------------------------------------------------------------------- timestep 0.002 #----------------------------------------------------------------------- # Πως να γράφεται η τροχιά. #----------------------------------------------------------------------- dump my_trajectory all custom 5 t.lammpstrj id x y z type dump_modify my_trajectory sort id #----------------------------------------------------------------------- # Αρχικός υπολογισμός ενέργειας και θερμοδυναμικών ποσοτήτων. #----------------------------------------------------------------------- run 0 #----------------------------------------------------------------------- # Αρχική ισορροπία του συστήματος. #----------------------------------------------------------------------- ##run 10000 #----------------------------------------------------------------------- # Παραγωγή τροχιάς. #----------------------------------------------------------------------- ##run 10000