PROGRAM MD_2D ******************************************************************* * Programm zur Durchfuehrung von Molekulardynamiksimuationen * mit dem Verlet Algorithmus. Propagation eines Metallteilchens * mit der Masse von Al (26,981538 amu) in einem periodischen * zweidimensionalen Potential * Zeitschritt 1fs, Gitterkonstante A= 3 Angstrom ****************************************************************** COMMON /DATA/ A,PI,POT0 OPEN (10,FILE="Trajektorie_Al.dat", & FORM="FORMATTED",STATUS="UNKNOWN") OPEN (25,FILE="Trajektorie_Al_all.dat", & FORM="FORMATTED",STATUS="UNKNOWN") POT0=1. A=3. PI=ACOS(-1.) AMASS=26.981538 DEL=1. C X0 x-Koordinate des Anfangspunktes in der Einheitszelle in Einheiten von a C Y0 y-Koordinate des Anfangspunktes in der Einheitszelle in Einheiten von a C EKIN kinetische Energie am Anfang in eV C ANG Winkel der Anfangsrichtung C AMASS Masse in amu des Teilchens C DEL Zeitschritt in fs C N Anzahl der Zeitschritte Write (*,*) "Molekulardynamiksimulationen" CAG WRITE (*,*) "Gebe die Anfangskoordinaten X0 und Y0 in der" CAG WRITE (*,*) "Einheitszelle an (0