Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
ΜΟΡΙΑΚΗ ∆ΥΝΑΜΙΚΗ ΑΝΑΛΥΣΗ<br />
156<br />
Παράµετρος Έκφραση<br />
Μήκος<br />
*<br />
r = r r0<br />
Χρόνος *<br />
t = t ( r0m D)<br />
Μάζα<br />
Ενέργεια<br />
4.3 Αρχικές συνθήκες<br />
*<br />
m m m<br />
= = 1<br />
*<br />
Ε =Ε<br />
*<br />
∆ύναµη F = F ( D r )<br />
Ταχύτητα *<br />
v = v ( D m)<br />
*<br />
Θερµοκρασία T = T ( D k )<br />
Πίνακας 4.1: Αδιάστατες εξισώσεις<br />
Το πρώτο βήµα σύµφωνα µε τη µέθοδο της Μ∆ ανάλυσης [1]-[3] για την δηµιουργία µιας<br />
προσοµοίωσης Μοριακής ∆υναµικής είναι οι αρχικές συνθήκες, η αρχική κατάσταση δηλαδή,<br />
οι οποίες περιγράφουν το µοντέλο. Η βάση για κάθε µοριακή προσοµοίωση είναι η<br />
προσοµοίωση του υλικού, το οποίο θα µελετηθεί. Πρώτα απαιτείται η περιγραφή της αρχικής<br />
θέσης κάθε ατόµου, η ατοµική δοµή δηλαδή του υλικού. Έπειτα, η αλληλεπίδραση<br />
(interaction) µεταξύ των ατόµων που περιέχονται στον όγκο προσοµοίωσης είναι επίσης µία<br />
σηµαντική παράµετρος που πρέπει να περιγραφεί. Η αλληλεπίδραση αυτή δεν είναι τίποτα<br />
άλλο παρά η δυναµική ενέργεια των ατόµων. Η δυναµική ενέργεια στην συγκεκριµένη<br />
µελέτη προσοµοιώνεται από την συνάρτηση δυναµικού του Morse (Morse Potential<br />
Function) που εξαρτάται από τις θέσεις (συντεταγµένες) των ατόµων, οι οποίες έχουν ήδη<br />
ορισθεί. Κατόπιν η κινητική ενέργεια κάθε ατόµου δίνεται από την κλασική εξίσωση της<br />
µηχανικής. Προφανώς η κινητική ενέργεια εξαρτάται από τις ταχύτητες που στο στάδιο της<br />
αρχικής κατάστασης θεωρούνται µηδενικές. Το άθροισµα της δυναµικής και κινητικής<br />
ενέργειας είναι η ολική ενέργεια του κάθε ατόµου και έτσι χρησιµοποιείται µια κατά<br />
Hamilton περιγραφή του συστήµατος, όπου οι συντεταγµένες του χώρου κατάστασης είναι οι<br />
θέσεις και οι ορµές των ατόµων του συστήµατος. Τέλος απαιτείται και ο ορισµός των<br />
συνοριακών συνθηκών, οι οποίες εξασφαλίζουν ότι το προσοµοιωµένο υλικό είναι στην<br />
κατάσταση που απαιτείται, συγκεκριµένα στην παρούσα φάση στην στερεά. Συνεπώς για την<br />
περιγραφή των αρχικών συνθηκών απαιτούνται:<br />
D<br />
0<br />
B