======================================= | OUTPUT OF PROGRAM "STATICS" | | --------------------------- | | (Yuri Mishin, Feb.-March 1996) | | [revised in June 1998] | ======================================= For comments, suggestions and bugs contact mishin@vt.edu REMEMBER: in this program the energy is measured in eV, distances in Angstroems, forces in eV/Angstroem, planar defect energy in eV/Angstroem^2 CHEMICAL ELEMENTS: TYPE ELEMENT 1 Fe 2 Fe THE FOLLOWING POTENTIAL FILES ARE USED: pfe-1.dat palni.dat pni_m.dat ffe-1.dat for021.dat Ffe-1.dat frni.dat CUT-OFF RADIUS OF ATOMIC INTERACTION: 5.8037 INITIAL CONFIGURATION IS GENERATED BY 'CRYSTAL' USING FILE: crystal.gen PRIMITIVE TRANSLATION VECTORS: 2.87000 0.00000 0.00000 0.00000 2.87000 0.00000 0.00000 0.00000 2.87000 NUMBER OF BASIS ATOMS = 2 COORDINATES ELEMENT 0.000000 0.000000 0.000000 Fe 1.435000 1.435000 1.435000 Fe EDGE VECTORS OF THE SIMULATION BLOCK: 1.00000 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 1.00000 NUMBER OF PLANES BEFORE LATTICE REPEATS IN BLOCK DIRECTIONS: 1 1 1 RADIUS OF SEARCH = 7.0000 OUTER BLOCK IS 17 BY 17 BY 17 BUFFER THICKNESS 3 BY 3 BY 3 INNER BLOCK IS 4 BY 4 BY 4 DIMENSIONS OF THE INNER BLOCK: 11.48000 11.48000 11.48000 THE PERFECT LATTICE BLOCK CONTAINS: TOTAL NUMBER OF ATOMS: NTOT = 9826 NUMBER OF FREE ATOMS: NFREE = 128 NUMBER OF BUFFER ATOMS: 1872 NUMBER OF FIXED ATOMS: 7826 --------------------------------------------------- COHESIVE ENERGY OF THE PERFECT LATTICE: -4.280000 ELECTRON DENSITY IN THE PERFECT LATTICE: LCO ELEMENT EL. DENSITY 1 Fe 1.000000 2 Fe 1.000000 -------------------------------------------------- LATTICE DEFECTS ARE GENERATED USING FILE: crystal.def THE BLOCK CONTAINS A VACANCY AT SITE 27 ======= COMPOSITION OF THE DEFECTED BLOCK: ======= THE BLOCK CONTAINS 127 ATOMS AND 1 VACANCIES 127 ATOMS ARE Fe 0 ATOMS ARE Fe =================================================== THE LINK CELL MAP IS 6 by 6 by 6 CONSTRUCTING TABLE OF NEIGHBOURS MAXIMUM NUMBER OF NEIGHBOURS IS 112 INITIAL ENERGY OF THE DEFECTED BLOCK: -0.8553769D+04 CRACK INFORMATION IS READ FROM FILE: crystal.crack N K_I ENERGY CRACK TIP POSITION NF 1 0.0000D+00 0.5988157D+01 X =-0.1000D+11 Y = 0.0000D+00 0 *** THE MINIMIZATION ALGORITHM SAFELY CONVERGED *** Congratulations !!! TOTAL NUMBER OF INTERATONS MADE: 10 TOTAL NUMBER OF FUNCTION CALLS: 21 MAXIMUM ATOMIC DISPLACEMENT DUE TO RELAXATION: 0.08651 MAXIMUM FORCE AFTER RELAXATION: 0.000184 ---------------------------------------------------- MINIMUM ENERGY OF THE BLOCK: -0.8554012D+04 DEFECT ENERGY: 0.5988157D+01 ---------------------------------------------------- THE FINAL CONFIGURATION HAS BEEN SAVED IN FILE: fevac.xyz ... END OF OUTPUT.