A Package for calculating elastic tensors of tetragonal Phases by using second-order derivative with Wien2k Package
|
|
- Alexandra Jody Douglas
- 5 years ago
- Views:
Transcription
1 IR ELAST + WIEN2k A Package for calculating elastic tensors of tetragonal Phases by using second-order derivative with Wien2k Package User s guide, Tetra-elastic_13.2 (Release ) Morteza Jamal Ghods City-Tehran-Iran 1
2 MANDATORY CONDITIONS: In any publication in the scientific literature please reference the program as follows: M. Jamal, Tetra-elastic, (2012). ACKNOWLEDGMENT I gratefully appreciate B.Z. Yanchitsky for fruitful discussions, P. Blaha and S. Jalali Asadabadi for suggestions, A. H. Reshak for calculations, and Carol Phillips for editing. For suggestions or bug reports please contact the author by m_jamal57@yahoo.com 2
3 1- Introduction Tetra-elastic is a Package for finding elastic constants of tetragonal symmetries with Wien2k. This Package calculates elastic constants by second-order derivative ( E (δ) ) of Polynomial fit ( E=E(δ) ) of Energy vs. strains (δ) at zero strain (δ =0). This called energy approach[1]. 2- Background theory (energy approach) Elastic constants are defined by means of a Taylor expansion of the total energy for the system, with respect to a small strain ( ) of the lattice. If we consider the bravais lattice vectors of a tetragonal crystal structure as a matrix form the distortion of the lattice ( ) is expressed by multi plying with a symmetric ( ) distortion matrix i.e. ( ), which is written as, And in Voigt notation ( It is often convenient to change to the Voigt notation in order to reduce the number of indices. The Voigt notation replaces 1, 2, 3, (and ) 4, (and ) 5, (and ) 6 ) 3
4 We express the energy of the strained system by means of a Taylor expansion in the distortion parameters, The linear terms vanish if the strain causes no changes in the volume of the crystal. Otherwise, are related to the strain on the crystal and are elastic constants and is the volume of unstrained tetragonal system and we use it to evaluate the elastic constants. To obtain elastic constants of tetragonal structure we have used the method has been obtained by Mehl et al [2]. In this method elastic constants were calculated by applying small strains to the unstrained lattice. There are six independent elastic constants for a tetragonal symmetry, called C 11, C 12, C 13, C 33, C 44, and C 66. Since we have six independent elastic constants, we need six different strains to determine these elastic constants. The six distortions used in the tetra-elastic Package are described below. The first three distortions are written as and These three distortions change the lattice parameter in the a, b and c directions. The symmetry of the strained lattice is therefore still tetragonal, however the volume of the distortion lattice changes by using D 1 and D 3 and the energy for these distortions can be obtained as 4
5 and The second three type of distortions are orthorhombic ( D 4 and monoclinic (D 5 ) distortions and written as D 6 ) and and and the energy for these distortions can be obtained as 5
6 and respectively. 6
7 3- File structure and program flow The following table describes input and output files for each program of the Tetra-elastic Package. Program needs generates T_set_elast_lapw case.struct init.struct runcommand1 runcommand2 pwdname command_init_lapw T_command_run_lapw T_setup11m12, T_setupc33 T_setupc1112, T_setupczz T_setupc44, T_setupc66 getcalljobt makestructt init.struct pwdname runcommand1 runcommand2 init.struct.styp TETRA.job auto_init_lapw runcommand1 runcommand2.vper.styp TETRA.job number.strain StypX_Y.struct vol.optimize T_modifyjob_lapw T_calljob_lapw TETRA.job VstVene T_fitdivELC number.strain ELCorder.fit vol.optimize ELC.output.styp ELC.fit VstVene T_ana_elastc_lapw VstVene case.outputeos ELC.fit vol.optimize ELC.ps T_ana_elast_lapw ELC.output ELC-matrix case.output_elastic T_InverseELC ELC-matrix INVELC-matrix MassRho init.struct.rho vol.optimize T_ana_elastorder_lapw ELCorder.fit output-order sgroupcheck_lapw StypX_Y.struct case.struct command_intso_lapw.infso command_initu_lapw.infldau auto_initso_lapw.infso case.inso auto_initu_lapw.infldau case.inorb case.indm/c Bold font is OPTIONAL Italic bold font means it is the user s choice 7
8 3-1- Short description for input and output files case.struct init.struct pwdname runcommand1/2 Is a Wien2k standard struct file. Is a copy of the case.struct file. Contains the name of the present work directory. Contains the run commands for running. It looks similar to: run_lapw ec p in1new 2 auto_init_lapw A C-shell program which automatically runs the initialization. It looks similar to:.vper.styp TETRA.job #!/bin/csh -f set RM = not if ( $RM == 'not' ) then init_lapw -vxc 13 -ecut -6 -mix 0.2 -numk b else init_lapw -red 0 -vxc 13 -ecut -6 -mix 0.2 -numk b endif Defines the percent of changes for different strains. Defines the type of strain. A C-shell program which calculates the energy for each strain by using the Wien2k Package. It looks similar to: #!/bin/csh -f #STRAIN TYPE IS 1 #Modify this script according to your needs unalias rm set co = 1 set name set bj set file = `pwd` set file = $file:t if (-e VstVene ) then set i=`/bin/ls VstVene* wc ` echo " saving pervious VstVene to VstVene_$i[1]" cp VstVene VstVene_$i[1] rm VstVene endif # # to reuse previous scf runs (without a new scf run) set answscf=y # and use the same "savename". # When you make modifications (RKmax, k-mesh, XC-potentials) choose # answscf=no and a new savename (eg. "_pbe_rk8_1000k"). set answscf=y set savename= 8
9 if (-e cscl.clmsum && \! -z cscl.clmsum) then x dstart -super endif if (-e cscl.clmup && \! -z cscl.clmup ) then x dstart -super -up x dstart -super -dn endif foreach i ( \ ) echo "*******************************" echo $i set name=`echo "$name $i"` echo "*******************************" Styp1_-3.0 \ Styp1_-2.0 \ Styp1_-1.0 \ Styp1 0.0 \ Styp1 1.0 \ Styp1 2.0 \ Styp1 3.0 \ StypX_Y.struct VstVene number.strain case.outputeos A Wien2k struct file for each value of changes and for each strain type where X and Y denote type of strain and value of changes, respectively. The main information file, contains values of changes (strains) and energies for each type of strain, for the calculation of the elastic constant. Contains the number of strains. A Wien2k output of equation of states(eos). For finding the best values of elastic constants, find EOS and then copy the case.outputeos file in the "case" directory within the c11+c12, c11-c12, c44, c33, czz, and c66 directories. Otherwise, it sets the optimized volume from the original struct file i.e. case.struct vol.optimize Contains the optimized volume. ELCorder.fit Contains the elastic constants for different values of order of fit. ELC.output Contains the elastic constants for order of fit =2 ELC.fit Contains the data to plot a curve of energy vs value of changes (strains) for each strain type. ELC-matrix Defines the elastic constant matrix for each symmetry. case.output_elastic Contains the final elastic constant values. INVELC-matrix Defines the inverse of elastic constant matrix. 9
10 .rho.infso.infldau Contains density of mass and atomic volume. Contains information for making the case.inso file for running spin-orbit coupling. Contains information for making the case.inorb and case.indm/c files for LDA+U calculations Flow and short description for programs The Tetra-elastic Package consists of several FORTRAN and SHELL SCRIPTS which are described below. A flowchart of the program is shown in the following diagram. T_set_elast_lapw : Makes an elast-constant directory in the present work directory ( PWD ) and c11+c12, c11-c12, c33, c44, c66, and czz directories in the elast-constant directory. The T_set_elast_lapw program also copies information of the "PWD" into the c11+c12, c11- c12, c33, c44, c66, and czz directories and calls "command_init_lapw", T_command_run_lapw, T_setupc1112, T_setupc11m12, T_setupc33, T_setupc44, T_setupc66, and T_setupczz programs.. command_init_lapw : Gets information for making "auto_init_lapw". T_command_run_lapw : Gets the run commands for making TETRA.job. T_setupcX (X=1112, 11m12, 44,.) : Gets the type of strain and calls the getcalljobt program. getcalljobt : Calls makestructt program and makes the TETRA.job file. makestructt : Makes the StypX_Y.struct files where X and Y stand for the type of strain and value of changes, respectively and the vol.optimize file. T_modifyjob_lapw : Edits the job files according to the user s needs. T_calljob_lapw : Calls the TETRA.job files for running. 10
11 T_ana_elast_lapw : Calls the T_ana_elastc_lapw program for calculating elastic constants then calculates the Voigt, Reuss, and Hill bulk, shear, and the Young modulus as well as the Poisson ratio. After that it calls the T_InverseELC and MassRho programs and calculates sound velocity and Debye temperature then makes two output files in the elastconstant directory with the name case.output_elastic and the INVELC-matrix which is the Elastic compliance constants generated by inverting the elastic constant matrix. At the end it calls T_ana_elastorder_lapw program. T_ana_elastc_lapw : Calls the T_fitdivELC program with appropriate libraries for calculating C11-C12, C11+C12, C44, and. T_InverseELC : Makes the Elastic compliance constants generated by inverting the elastic constant matrix. MassRho : Finds density of mass and atomic volume. T_ana_elastorder_lapw : Checks the sensitivity of the elastic constants to the order of fit. sgroupcheck_lapw : Finds the best value of tol in the sgroup[3] program and copies case.struct_sgroup as case.struct. 11
12 T_set_elast_lapw command_init_lapw generates auto_init_lapw T_command_run_lapw generates commandrun1/2 T_setupcX X=1112, 11m12,. getcalljobt generates TETRA.job makestructt generates StypX_Y.struct sgroupcheck_lapw TETRA.job generates VstVene auto_init_lapw T_ana_elastorder_lapw generates ELCorder.fit T_ana_elast_lapw T_ana_elastc_lapw T_InverseELC MassRho ELASTIC CONSTANTS IS READY generates case.output_elastic and INVELC-matrix 12 T_fitdivELC calls Libraries generates C 11 -C 12, C 11 +C 12,. Program flow in Tetra-elastic Dash arrow means user must run
13 4 Elastic constants calculation 1. Create a struct file and validate it by running "sgroupcheck_lapw". 2. If Spin-Orbit calculations are required run "command_initso_lapw. 3. If LDA+U calculations are required run "command_initu_lapw" and then "auto_initu_lapw". 4. Run "T_set_elast_lapw.. 5. Now you must adapt the job files according to your needs (you can run "T_modifyjob_lapw" in Terminal ). It is not necessary to do step 5 if you defined the COMMAND RUN commands in step Now you must run the job files (you can run "T_calljob_lapw" ). It will take some time. 7. Run "T_ana_elast_lapw".. This package calculates elastic constants by second-order derivative ( E (δ) ) of Polynomial fit ( E=E(δ) ) of Energy vs. strains ( є) at zero strain (δ =0) so, you must use values of strain around zero and from the viewpoint of fit convergence, we usually expect to see a minimum when we plot Energy vs. strain ( this Package plots it ). It is recommended that the sensitivity of the results is checked to the order of fit. This program shows them. 4-1 Notes about elastic constants calculation After using distortions for the calculation of C11-C12, C44, and C66 the symmetry of the tetragonal compound changes and usually the number of atoms change. So when you run "command_initso_lapw" or "command_initu_lapw", in the section name of an atom, type "all <name of atom>" ( for example: all Mn). With this command, you use SO or LDA+U calculations for example for all Mn atoms. When you want to rerun job files with modifications in (RKmax, k-mesh, XC-potentials ) call the command_init_lapw and after that choose "answscf=no" in the TETRA.job files and a new "savename" (eg. "_use_pbe_rk8"). 13
14 Optionally you can specify more cases by rerunning T_setupcX (X=1112, 11m12, 44, see section 5-3 ). Specify also your old cases. The old results will then be taken automatically into account without recalculation (unless you modify job files i.e: set answscf=no ). For the calculation of the best values of elastic constants, please find EOS and then copy case.outputeos in the "case" directory within the c11+c12, c11-c12, c44, and directories. Otherwise, it sets the optimized volume from the original struct file i.e. case.struct. 4-2 One calculation To calculate C 33 or C 44 calculation of C 44 : the following steps should be performed for example for 1. Make a directory for example c Make a "case" directory in c44 directory. 3. Make a "case.struct" file in the "case" directory and name it "init.struct". Create a "pwdname" file and write in it "case." and save it. 4. Run the command_init_lapw 5. chmod +x auto_init_lapw 6. For SO calculations, run the command_initso_lapw. 7. For LDA+U calculations, run the command_initu_lapw and auto_initu_lapw To avoid step 10, you can run T_command_run_lapw for setting the COMMAND RUN commands for making TETRA.job. 8. Run the T_setupc44 program. 9. chmod +x TETRA.job file. 10. Modify the TETRA.job file. 11. Call TETRA.job 12. Call T_ana_elastc_lapw 4-3 Run with more data points Optionally you can specify more data points, for the calculation of the elastic constants, by rerunning T_setupcX (X=1112, 11m12, 44, ). Specify also your old data points. The old results will then be taken automatically into account without recalculation ( unless you modify the job files i.e: set answscf=no ). Please do the following steps for this goal for example for c Cd to the elast-constant directory. 2. Cd to the c44 directory. 3. cd the case directory. 3-1) To avoid step 6, you can run T_command_run_lapw for setting the 14
15 COMMAND RUN commands for making the TETRA.job. 4. Run the T_setupc44 program. 5. If you want to rerun the job files with modifications in (RKmax, k-mesh, XCpotentials ) call command_init_lapw and then choose "answscf=no" in TETRA.job files and a new "savename" (eg. "_use_pbe_rk8"). 6. Modify the TETRA.job file. 7. Call TETRA.job 8. Call T_ana_elastc_lapw 4-4 Elastic constants calculation for α-pt 2 Si α-pt 2 Si compound is a test case for elastic constants calculation. The α-pt 2 Si structure is described in detail in the following: alpha Pt2Si B LATTICE,NONEQUIV.ATOMS: 2139_I4/mmm MODE OF CALC=RELA unit=bohr ATOM -1: X= Y= Z= MULT= 2 ISPLIT=-2-1: X= Y= Z= Pt NPT= 781 R0= RMT= Z: 78.0 LOCAL ROT MATRIX: ATOM -2: X= Y= Z= MULT= 1 ISPLIT=-2 Si NPT= 781 R0= RMT= Z: 14.0 LOCAL ROT MATRIX: NUMBER OF SYMMETRY OPERATIONS Select Xc = PBE-GGA, R_Kmax = 7, L_max = 8, and nkpoint = 5000 In the following examples you can find the percents that were used for strains. 15
16 ####################################### # T_ana_elast_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # ####################################### ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ########################################## =============================================================== Order of fit: 2 C11+C12 is: GPa, RMS: E-03 Order of fit: 3 C11+C12 is: GPa, RMS: E-04 Order of fit: 4 C11+C12 is: GPa, RMS: E-04 Order of fit: 5 C11+C12 is: GPa, RMS: E-10 ****************************************** Polynomial fit for C11+C12 done A RMS of E-03 was achieved using a polynome of degree : 2 At volume= bohr^3 C11+C12 is: a.u or GPa ****************************************** Analyze done... Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ##########################################
17 =============================================================== Order of fit: 2 C33 is: GPa, RMS: E-04 Order of fit: 3 C33 is: GPa, RMS: E-04 Order of fit: 4 C33 is: GPa, RMS: E-04 Order of fit: 5 C33 is: GPa, RMS: E-04 Order of fit: 6 C33 is: GPa, RMS: E-10 ****************************************** Polynomial fit for C33 done A RMS of E-04 was achieved using a polynome of degree : 2 At volume= bohr^3 C33 is: a.u or GPa ****************************************** Analyze done... Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ########################################## =============================================================== Order of fit: 2 Czz is: GPa, RMS: E-03 Order of fit: 3 Czz is: GPa, RMS: E-04 Order of fit: 4 Czz is: GPa, RMS: E-11 ****************************************** Polynomial fit for Czz done A RMS of E-03 was achieved using a polynome of degree : 2 At volume= bohr^3 Czz is: a.u or GPa ****************************************** Analyze done... Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** 17
18 ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ########################################## =============================================================== Order of fit: 2 C11-C12 is: GPa, RMS: E-04 Order of fit: 3 C11-C12 is: GPa, RMS: E-04 Order of fit: 4 C11-C12 is: GPa, RMS: E+00 ****************************************** Polynomial fit for C11-C12 done A RMS of E-04 was achieved using a polynome of degree : 2 At volume= bohr^3 C11-C12 is: a.u or GPa ****************************************** Analyze done... Press <enter> to continue Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ########################################## =============================================================== Order of fit: 2 C44 is: GPa, RMS: E-04 Order of fit: 3 C44 is: GPa, RMS: E-04 Order of fit: 4 C44 is: GPa, RMS: E-10 ****************************************** Polynomial fit for C44 done A RMS of E-04 was achieved using a polynome of degree : 2 18
19 At volume= bohr^3 C44 is: a.u or GPa ****************************************** Analyze done... Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** ########################################## # T_ana_elastc_lapw analyses Elastic # # constant # # C(2012) by Morteza Jamal # # using case.outputeos # # VstVene # # which have been created by # # TETRA.job # ########################################## =============================================================== Order of fit: 2 C66 is: GPa, RMS: E-05 Order of fit: 3 C66 is: GPa, RMS: E-05 Order of fit: 4 C66 is: GPa, RMS: E-10 ****************************************** Polynomial fit for C66 done A RMS of E-05 was achieved using a polynome of degree : 2 At volume= bohr^3 C66 is: a.u or GPa ****************************************** Analyze done... Do you want a hardcopy? (y/n) *************************************** You can find data in ELC.output file. *************************************** Printing final Elastic constant At voulme= bohr^3. ======================================================================= C11+C12 = GPa C11-C12 = GPa C33 = GPa C44 = GPa Czz = C11+C12+2C33-4C13 = GPa C66 = GPa ======================================================================= Atom name = Pt Atomic Mass from Periodic table = (gr/mol) 19
20 Atomic Mass from Periodic table = *10^(-23) (gr) Atom name = Si Atomic Mass from Periodic table = (gr/mol) Atomic Mass from Periodic table = *10^(-23) (gr) Volume in unit of cm^3 = *10^(-24) (cm^3) Mass of Compound : *10^(-23) (gr) Density of Compound : (gr/cm^3) ======================================================================= C11 = GPa C12 = GPa C13 = GPa C33 = GPa C44 = GPa C66 = GPa Prediction VOIGT Bulk modulus by using elastic constant values Prediction REUSS Bulk modulus by using elastic constant values Prediction HILL Bulk modulus by using elastic constant values = (GPa) = (GPa) = (GPa) Prediction VOIGT Shear modulus by using elastic constant values = (GPa) Prediction REUSS Shear modulus by using elastic constant values = (GPa) Prediction HILL Shear modulus by using elastic constant values = (GPa) Prediction VOIGT Young modulus by using elastic constant values = (GPa) Prediction REUSS Young modulus by using elastic constant values = (GPa) Prediction HILL Young modulus by using elastic constant values = (GPa) Prediction VOIGT Poisson's coefficient by using elastic constant values =.334 Prediction REUSS Poisson's coefficient by using elastic constant values =.356 Prediction HILL Poisson's coefficient by using elastic constant values =.345 ======================================================================= By using HILL data Transverse elastic wave velocity = (m/s) Longitudinal elastic wave velocity = (m/s) The average wave velocity = (m/s) Debye Temperature = (K) ======================================================================= Melting Temperature = K or K ======================================================================= Press enter key to continue... ################################################ # T_ana_elastorder_lapw checks the sensitivity # # of elastic constants to the order of fit # # by using ELCorder.fit file # # C(2012) by Morteza Jamal # # # ################################################ CHECK THE SENSITIVITY OF YOUR RESULT TO THE ORDER OF FIT 20
21 Press enter key to continue... Order of fit for calculations were 4,6,4,5, 4, and 4 We select minimum value for ORDER OF FIT i.e. 4 Press enter key to continue... ######## ORDER OF FIT IS : 2, At volume = (bohr^3) ######## (c11-c12) = (GPa) (c11+c12) = (GPa) (c33) = (GPa) (c44) = (GPa) (czz=c11+c12+2c33-4c13) = (GPa) (c66)= (GPa) c11 = (GPa) c12 = (GPa) c13 = (GPa) c33 = (GPa) c44 = (GPa) c66 = (GPa) ======================================================================= Prediction VOIGT Bulk modulus by using elastic constant values = (GPa) Prediction REUSS Bulk modulus by using elastic constant values = (GPa) Prediction HILL Bulk modulus by using elastic constant values = (GPa) Prediction VOIGT Shear modulus by using elastic constant values = (GPa) Prediction REUSS Shear modulus by using elastic constant values = (GPa) Prediction HILL Shear modulus by using elastic constant values = (GPa) Prediction VOIGT Young modulus by using elastic constant values = (GPa) Prediction REUSS Young modulus by using elastic constant values = (GPa) Prediction HILL Young modulus by using elastic constant values = (GPa) Prediction VOIGT Poisson's coefficient by using elastic constant values =.334 Prediction REUSS Poisson's coefficient by using elastic constant values =.356 Prediction HILL Poisson's coefficient by using elastic constant values =.345 ======================================================================= ######## ORDER OF FIT IS : 3, At volume = (bohr^3) ######## (c11-c12) = (GPa) (c11+c12) = (GPa) (c33) = (GPa) (c44) = (GPa) (czz=c11+c12+2c33-4c13) = (GPa) (c66)= (GPa) c11 = (GPa) c12 = (GPa) c13 = (GPa) c33 = (GPa) c44 = (GPa) c66 = (GPa) ======================================================================= Prediction VOIGT Bulk modulus by using elastic constant values = (GPa) Prediction REUSS Bulk modulus by using elastic constant values = (GPa) Prediction HILL Bulk modulus by using elastic constant values = (GPa) 21
22 Prediction VOIGT Shear modulus by using elastic constant values = (GPa) Prediction REUSS Shear modulus by using elastic constant values = (GPa) Prediction HILL Shear modulus by using elastic constant values = (GPa) Prediction VOIGT Young modulus by using elastic constant values = (GPa) Prediction REUSS Young modulus by using elastic constant values = (GPa) Prediction HILL Young modulus by using elastic constant values = (GPa) Prediction VOIGT Poisson's coefficient by using elastic constant values =.346 Prediction REUSS Poisson's coefficient by using elastic constant values =.366 Prediction HILL Poisson's coefficient by using elastic constant values =.356 ======================================================================= ######## ORDER OF FIT IS : 4, At volume = (bohr^3) ######## (c11-c12) = (GPa) (c11+c12) = (GPa) (c33) = (GPa) (c44) = (GPa) (czz=c11+c12+2c33-4c13) = (GPa) (c66)= (GPa) c11 = (GPa) c12 = (GPa) c13 = (GPa) c33 = (GPa) c44 = (GPa) c66 = (GPa) ======================================================================= Prediction VOIGT Bulk modulus by using elastic constant values = (GPa) Prediction REUSS Bulk modulus by using elastic constant values = (GPa) Prediction HILL Bulk modulus by using elastic constant values = (GPa) Prediction VOIGT Shear modulus by using elastic constant values = (GPa) Prediction REUSS Shear modulus by using elastic constant values = (GPa) Prediction HILL Shear modulus by using elastic constant values = (GPa) Prediction VOIGT Young modulus by using elastic constant values = (GPa) Prediction REUSS Young modulus by using elastic constant values = (GPa) Prediction HILL Young modulus by using elastic constant values = (GPa) Prediction VOIGT Poisson's coefficient by using elastic constant values =.313 Prediction REUSS Poisson's coefficient by using elastic constant values =.322 Prediction HILL Poisson's coefficient by using elastic constant values =.318 ======================================================================= You can find these data in the output-order file. C11 = GPa C12 = GPa C13 = GPa C33 = GPa C44 = 71.0 GPa C66 = Gpa 22
23 α-pt 2 Si 4 (others) α-pt 2 Si (our Cal) Method FPLMTO FPLAPW Exchange LDA PBE a (Å) c (Å) C C C C C C
24 5 Installation of the Tetra-elastic package The Tetra-elastic package comes as a compressed tar file namely tetra-elastic-deriativemethod-2012-p.tar.gz or tetra-elastic.tar.gz. To install the package firstly copy the file to a directory of your choice. Now, uncompress and expand it as: tar zxvf tetra-elastic-deriative-method-2012-p.tar.gz cd tetra-elastic 24
25 Run buildtirelast_lapw This program helps you to create the "Makefile" and then compile tetra-elastic. By default, the Makefile expects the lapack_lapw and blas_lapw to be in the location../src_lib. This should be changed to the correct location by modifying the FOPT parameter as shown below. This Program helps you to define Fortran compiler, Fortran options, and Library options if you have installed WIEN2k. As you can see here this program defines Fortran compiler, Fortran options, and Library options as automatically. Otherwise you can define compiler and linker options as well as the path of mkl library depending on the your selected system. /home/mylib/mkl/lib/em64t is the path of my mkl library. To make Makefile by the lapack_lapw and blas_lapw libraries in the location../src_lib and gfortran use the following options: Fortran compiler: gfortran Fortran options: -ffree-form Library options (Lapack and BLAS): $(FOPT) L/home/physicsprogram/SRC_lib lpthread static llapack_lapw lblas_lapw the location../src_lib should be changed to the correct location by modifying the FOPT parameter as shown above. PS: To install with -ffree-form, you should compile the lapack_lapw and blas_lapw libraries with -ffree-form options. Otherwise it might was caused error. 25
26 If you view the OPTIONS file of the WIEN2k package you can use the FOPT, LDFLAGS, and R_LIBS of it for compiling. After defining the Fortran compiler, Fortran options, and Library options press Enter key. 26
27 The Environment Variable ELASTT_PATH is then defined and added to the end of the.bashrc file. Thus you will be able to call tetra-elastic s programs for any location. 27
28 If you view the.bashrc file you can see Now, logout from your Linux system and then login. 28
29 6 References [1] R. Stadler, W. Wolf, R. Podloucky, G. Kresse, J. Furthmller, J. Hafner, Phys. Rev. B 54 (1996) [2] M. J. Mehl, J. E. Osburn, D. A. Papaconstantopoulus, B. M. Klein, Phys. Rev. B 41 (1990) [3] B. Z. Yanchitsky, A. N. Timoshevskii, Determination of the space group and unit cell for a periodic solid, Comp. Phys. Comm. 139 (2001) [4] O. Beckstein, J. E. Klepeis, G. L. W. Hart, and O. Pankratov, Phys. Rev. B 64 (2001)
A Package for calculating elastic tensors of rhombohedral Phases by using second-order derivative with Wien2k Package
IR ELAST + WIEN2k A Package for calculating elastic tensors of rhombohedral Phases by using second-order derivative with Wien2k Package User s guide, Rhom-elastic_13.1 (Release 16.09.2013) Morteza Jamal
More informationA Package for calculating elastic tensors of cubic Phases by using second-order derivative with WIEN2k Package
IR E L A S T + W I E N 2k A Package for calculating elastic tensors of cubic Phases by using second-order derivative with WIEN2k Package User s guide, Cubic-elastic_13.2 (Release 27.08.2013) Morteza Jamal
More informationA Package for calculating elastic tensors of hexagonal Phases by using second-order derivative with Wien2k Package
IR E L A S T + W I E N 2k A Package for calculating elastic tensors of hexagonal Phases by using second-order derivative with Wien2k Package User s guide, Hex-elastic_13.2 (Release 27.08.2013) Morteza
More informationThis package performs a convenient 2-Dimensional structure optimization ( Volume and c/a ) for hexagonal or tetragonal spacegroups.
Package 2Doptimize -------------------------- This package performs a convenient 2-Dimensional structure optimization ( Volume and c/a ) for hexagonal or tetragonal spacegroups. If you want to use this
More informationCOMPUTATIONAL TOOL. Fig. 4.1 Opening screen of w2web
CHAPTER -4 COMPUTATIONAL TOOL Ph.D. Thesis: J. Maibam CHAPTER: 4 4.1 The WIEN2k code In this work, all the calculations presented are performed using the WIEN2k software package (Blaha et al., 2001). The
More informationColloque National sur les Techniques de Modélisation et de Simulation en Science des Matériaux, Sidi Bel-Abbès Novembre 2009
Colloque National sur les Techniques de Modélisation et de Simulation en Science des Matériaux, Sidi Bel-Abbès. 23-24 Novembre 2009 Elastic, electronic and optical properties of SiGe 2N 4 under pressure
More informationExercises: In the following you find some suggestions for exercises, which teach you various tasks one may perform with WIEN2k.
Exercises: In the following you find some suggestions for exercises, which teach you various tasks one may perform with WIEN2k. New WIEN2k users should start with the first basic exercises (1-5), covering:
More informationAnalysis of an Electric-Field Gradient (EFG): the EFG-switch in LAPW2
Analysis of an Electric-Field Gradient (EFG): the EFG-switch in LAPW2 Katrin Koch Stefaan Cottenier August 10, 2011 0.1 In brief: the EFG and how it is obtained in LAPW The EFG is a traceless symmetric
More informationComputational Material Science Part II-1: introduction. Horng-Tay Jeng ( 鄭弘泰 ) Institute of Physics, Academia Sinica
Computational Material Science Part II-1: introduction Horng-Tay Jeng ( 鄭弘泰 ) Institute of Physics, Academia Sinica Outline Introduction of Computational Material Science (CMS) Density Functional Theory
More informationGaAs -- MLWF. Special thanks to Elias Assmann (TU Graz) for the generous help in preparation of this tutorial
GaAs -- MLWF + + Special thanks to Elias Assmann (TU Graz) for the generous help in preparation of this tutorial YouTube video: https://youtu.be/r4c1yhdh3ge 1. Wien2k SCF Create a tutorial directory, e.g.
More informationChemical bonding, elasticity, and valence force field models: A case study for -Pt 2 Si and PtSi
PHYSICAL REVIEW B, VOLUME 64, 550 Chemical bonding, elasticity, and valence force field models: A case study for -Pt 2 Si and PtSi J. E. Klepeis, O. Beckstein, 2 O. Pankratov, 2 and G. L. W. Hart 3 Lawrence
More informationAs an example we consider the stacking and analysis of a 2-ply symmetric laminate. First we clear the Matlab space and close figures.
Program Lam 1 Introduction The collection of Matlab command and function files in the package lam allows the stacking and analysis of a laminate. In fact only a representative piece of a laminate is considered,
More informationEOS-FIT V6.0 R.J. ANGEL
EOS-FIT V6. R.J. AGEL Crystallography Laboratory, Dept. Geological Sciences, Virginia Tech, Blacksburg, VA46, USA http://www.geol.vt.edu/profs/rja/ ITRODUCTIO EosFit started as a program to fit equations
More informationExercises: In the following you find some suggestions for exercises, which teach you various tasks one may perform with WIEN2k.
Exercises: In the following you find some suggestions for exercises, which teach you various tasks one may perform with WIEN2k. Please note, that often calculational parameters are set to minimal cpu-time
More informationPredicting the Structure of Solids by DFT
Questions? Hudson Hall 235 or Hudson Hall 1111 Predicting the Structure of Solids by DFT Hands-On Instructions Contents 1 Cohesive Energy for Bulk Phases of Si 11 Setting up the Structures 12 Structure-Dependent
More informationSpin transport in Magnetic Tunnel Junctions
Spin transport in Magnetic Tunnel Junctions Tutorial on spin transport in Fe-MgO-Fe Version 2015.2 Spin transport in Magnetic Tunnel Junctions: Tutorial on spin transport in Fe-MgO-Fe Version 2015.2 Copyright
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics Elasticity M.P. Vaughan Overview Overview of elasticity Classical description of elasticity Speed of sound Strain Stress Young s modulus Shear modulus Poisson ratio
More informationconnect/setup of w2web environment
connect/setup of w2web environment Connect to pleiades.bc.edu using NoMachine as indicated in the instructions and with the username/pw you got during registration. open a terminal and connect to your
More information1.1 Contribution to the total energy For the LDA+U and orbital polarization methods a correction to the total energy appears E corr = E orb T r(v orb
Orbital package in WIEN P. Novak Institute of Physics, Prague, Czech Republic, novakp@fzu.cz April 2002 1 Description of SRC orb package Orbital package calculates the orbitally dependent potentials V
More informationExample: H 2 O (the car file)
Example: H 2 O (the car file) As a practical example of DFT methods we calculate the energy and electronic properties of the water molecule. In order to carry out the DFT calculation you will need a set
More informationLecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc.
Lecture 4 Implementing material models: using usermat.f Implementing User-Programmable Features (UPFs) in ANSYS 1 Lecture overview What is usermat.f used for? Stress, strain and material Jacobian matrix
More informationOPENATOM for GW calculations
OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW
More informationA two-dimensional FE truss program
A two-dimensional FE truss program 4M020: Design Tools Eindhoven University of Technology Introduction The Matlab program fem2d allows to model and analyze two-dimensional truss structures, where trusses
More informationPROMAL2012 SOFTWARE PACKAGE A USER GUIDE
PROMAL2012 SOFTWARE PACKAGE A USER GUIDE 1. This manual is only for VISTA, WINDOWS 7 and WINDOWS 8 users. The PROMAL2012 software and manual are available at http://www.eng.usf.edu/~kaw/promal2012/ 2.
More informationMustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A.
Ab Initio Studies On Phase Behavior of Barium Titanate Mustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A. 1 Physics Department,
More informationPWSCF First examples
PWSCF First examples (much more in espresso 3.1.1/examples directory!) Guido Fratesi (Università di Milano) Urbana, August 2006 A pw.x input file &CONTROL / &SYSTEM / &ELECTRONS title calculation restart_mode
More informationarxiv:cond-mat/ v2 [cond-mat.mtrl-sci] 10 Jun 2001
First principles elastic constants and electronic structure of α-pt 2 Si and PtSi arxiv:cond-mat/0008200v2 [cond-mat.mtrl-sci] 0 Jun 200 O. Beckstein, () J. E. Klepeis, (2) G. L. W. Hart, (3) and O. Pankratov
More informationTable of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures
Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state
More informationMcMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007
McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007 RANDOM POLYCRYSTALS OF GRAINS WITH CRACKS: MODEL OF ELASTIC BEHAVIOR FOR FRACTURED SYSTEMS James G. Berryman Earth Sciences Division
More informationSpin-orbit coupling in Wien2k
Spin-orbit coupling in Wienk Robert Laskowski rolask@ihpc.a-star.edu.sg Institute of High Performance Computing Singapore Dirac Hamiltonian Quantum mechanical description of electrons, consistent with
More informationLindgren CRYSTAL SYMMETRY AND ELASTIC CONSTANTS MICHAEL WANDZILAK. S.B., Massachusetts Institute of Technology (196'7)
CRYSTAL SYMMETRY AND ELASTIC CONSTANTS by MICHAEL WANDZILAK S.B., Massachusetts Institute of Technology (196'7) Submitted in partial fulfillment of the requirements for the degree of Master of Science
More informationInvestigation of mechanical stability of possible structures of PtN using first-principles computations
Investigation of mechanical stability of possible structures of PtN using first-principles computations S. K. R. Patil, S. V. Khare, * B. R. Tuttle, 3 J. K. Bording, 4 and S. Kodambaka 5 Department of
More informationEBIS-PIC 2D. 2D EBIS Simulation Code. Version 0.1. Copyright ( ) January FAR-TECH, Inc Science Center Dr., Ste 150. San Diego CA 92121
EBIS-PIC 2D 2D EBIS Simulation Code Version 0.1 Copyright ( ) January 2017 by 10350 Science Center Dr., Ste 150 San Diego CA 92121 Phone 858-455-6655 Email support@far-tech.com URL http://www.far-tech.com
More informationCHAPTER: 8. ELECTRONIC STRUCTURE AND ELASTIC PROPERTIES OF CrC AND CrN. 8.1 Introduction. Ph.D. Thesis: J. Maibam
CHAPTER -8 CHAPTER: 8 ELECTRONIC STRUCTURE AND ELASTIC PROPERTIES OF CrC AND CrN 8.1 Introduction In this chapter, we have selected CrC and CrN from group VIB transition metal carbides and nitrides for
More informationStrain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity
Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Torino, Italy, September 4-9, 2016 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it
More informationTABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA TABLE OF CONTENTS 1. INTRODUCTION TO COMPOSITE MATERIALS 1.1 Introduction... 1.2 Classification... 1.2.1
More information3rd International Conference on Machinery, Materials and Information Technology Applications (ICMMITA 2015)
3rd International Conference on Machinery, Materials and Information Technology Applications (ICMMITA 2015) Study on Energy Band-gap Calculation of CuGaS 2 Liuyang YUa, Yong Xub, Kegao LIUc* School of
More informationCHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50
CHAPTER 3 WIEN2k WIEN2k is one of the fastest and reliable simulation codes among computational methods. All the computational work presented on lanthanide intermetallic compounds has been performed by
More informationFIRST PRINCIPLES STUDY OF ELECTRONIC, ELASTIC AND OPTICAL PROPERTIES OF SnO UNDER IN-PLANE BIAXIAL STRAIN
Digest Journal of Nanomaterials and Biostructures Vol. 14, No. 1, January - March 2019, p. 127-137 FIRST PRINCIPLES STUDY OF ELECTRONIC, ELASTIC AND OPTICAL PROPERTIES OF SnO UNDER IN-PLANE BIAXIAL STRAIN
More informationMaterials that you may find helpful when working through this exercise
Detailed steps illustrating how to use VASP on the Suns in Fitz 177 For use in lab: 11/10/2009 (Original file by Dr. Rachel Getman, 11/18/2007. Editted for use by Dorrell McCalman 11/09/2009.) Note on
More informationOn Dynamic and Elastic Stability of Lanthanum Carbide
Journal of Physics: Conference Series On Dynamic and Elastic Stability of Lanthanum Carbide To cite this article: B D Sahoo et al 212 J. Phys.: Conf. Ser. 377 1287 Recent citations - Theoretical prediction
More informationThe Augmented Spherical Wave Method
Introduction Institut für Physik, Universität Augsburg Electronic Structure in a Nutshell Outline 1 Fundamentals Generations 2 Outline 1 Fundamentals Generations 2 Outline Fundamentals Generations 1 Fundamentals
More informationDynamics of the Atmosphere GEMPAK Lab 3. 3) In-class exercise about geostrophic balance in the real atmosphere.
Dynamics of the Atmosphere GEMPAK Lab 3 Goals of this lab: 1) Learn about Linux scripts. 2) Learn how to combine levels in GEMPAK functions. 3) In-class exercise about geostrophic balance in the real atmosphere.
More informationAb Initio Study of Electronic, Structural, Thermal and Mechanical Characterization of Cadmium Chalcogenides 65
Ab Initio Study of Electronic, Structural, Thermal and Mechanical Characterization of Cadmium Chalcogenides 65 Devi Prasadh P.S. 1, a, B.K. Sarkar 2, b 1 Department of Physics, Dr. Mahalingam College of
More informationELECTRONIC AND MAGNETIC PROPERTIES OF BERKELIUM MONONITRIDE BKN: A FIRST- PRINCIPLES STUDY
ELECTRONIC AND MAGNETIC PROPERTIES OF BERKELIUM MONONITRIDE BKN: A FIRST- PRINCIPLES STUDY Gitanjali Pagare Department of Physics, Sarojini Naidu Govt. Girls P. G. Auto. College, Bhopal ( India) ABSTRACT
More informationElasticity Constants of Clay Minerals Using Molecular Mechanics Simulations
Elasticity Constants of Clay Minerals Using Molecular Mechanics Simulations Jin-ming Xu, Cheng-liang Wu and Da-yong Huang Abstract The purpose of this paper is to obtain the elasticity constants (including
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationWannier functions. Macroscopic polarization (Berry phase) and related properties. Effective band structure of alloys
Wannier functions Macroscopic polarization (Berry phase) and related properties Effective band structure of alloys P.Blaha (from Oleg Rubel, McMaster Univ, Canada) Wannier functions + + Wannier90: A Tool
More informationA DFT Study on Electronic Structures and Elastic Properties of AgX (X=C, N) in Rock Salt Structure
Invertis Journal of Jameson Science Maibam, and Technology, Kh. Kabita, Vol. B. Indrajit 7, No. 2, Sharma, 2014. R.K. ; pp. Thapa 114-118 and R.K. Brojen Singh A DFT Study on Electronic Structures and
More informationDMDW: A set of tools to calculate Debye-Waller factors and other related quantities using dynamical matrices.
DMDW: A set of tools to calculate Debye-Waller factors and other related quantities using dynamical matrices. DMDW is a set of tools developed to calculate Debye-Waller (DW) factors and other related quantities
More informationTable of Contents. Table of Contents Initialize from a converged state. Introduction Examples
Table of Contents Table of Contents Initialize from a converged state Introduction Examples Water molecule example Looping over parameters Stepping up in bias Special example: anti-parallel spin in MTJs
More informationarxiv:physics/ v2 [physics.atom-ph] 31 May 2004
arxiv:physics/0405136v2 [physics.atom-ph] 31 May 2004 Pure spin angular momentum coefficients for non scalar one particle operators in jj coupling G. Gaigalas a and S. Fritzsche b a Institute of Theoretical
More informationFIRST PRINCIPLES STUDY OF AlBi
Available at: http://publications.ictp.it IC/2008/025 United Nations Educational, Scientific and Cultural Organization and International Atomic Energy Agency THE ABDUS SALAM INTERNATIONAL CENTRE FOR THEORETICAL
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationModule 2: Quantum Espresso Walkthrough
Module 2: Quantum Espresso Walkthrough Energy and Geometry Optimization of the H 2 Molecule We will be using the PWSCF code for quantum mechanical calculations of extended systems. The PWSCF program is
More informationMATERIAL ELASTIC ANISOTROPIC command
MATERIAL ELASTIC ANISOTROPIC command.. Synopsis The MATERIAL ELASTIC ANISOTROPIC command is used to specify the parameters associated with an anisotropic linear elastic material idealization. Syntax The
More informationSpacegroup P4 2 /mnm. Structure given by: spacegroup lattice parameter positions of atoms (basis)
Spacegroup P4 2 /mnm Structure given by: spacegroup lattice parameter positions of atoms (basis) Rutile TiO 2 : P4 2 /mnm (136) a=8.68, c=5.59 bohr Ti: (0,0,0) O: (0.304,0.304,0) Specify: Structure generator
More informationand strong interlayer quantum confinement
Supporting Information GeP3: A small indirect band gap 2D crystal with high carrier mobility and strong interlayer quantum confinement Yu Jing 1,3, Yandong Ma 1, Yafei Li 2, *, Thomas Heine 1,3 * 1 Wilhelm-Ostwald-Institute
More informationSupporting Information. Potential semiconducting and superconducting metastable Si 3 C. structures under pressure
Supporting Information Potential semiconducting and superconducting metastable Si 3 C structures under pressure Guoying Gao 1,3,* Xiaowei Liang, 1 Neil W. Ashcroft 2 and Roald Hoffmann 3,* 1 State Key
More informationCrystalline and Magnetic Anisotropy of the 3d Transition-Metal Monoxides
Crystalline and of the 3d Transition-Metal Monoxides Institut für Festkörpertheorie und -optik Friedrich-Schiller-Universität Max-Wien-Platz 1 07743 Jena 2012-03-23 Introduction Crystalline Anisotropy
More informationLab 3: Handout Quantum-ESPRESSO: a first principles code, part 2.
1 Lab 3: Handout Quantum-ESPRESSO: a first principles code, part 2. In this lab, we will be using Quantum-ESPRESSO as our first-principles code again. In problem 1, we will compare energy between allotropes
More informationSpeed-up of ATK compared to
What s new @ Speed-up of ATK 2008.10 compared to 2008.02 System Speed-up Memory reduction Azafulleroid (molecule, 97 atoms) 1.1 15% 6x6x6 MgO (bulk, 432 atoms, Gamma point) 3.5 38% 6x6x6 MgO (k-point sampling
More informationStatus report on the showering of Alpgen events with Herwig++ for ATLAS
Status report on the showering of Alpgen events with Herwig++ for ATLAS University of Oslo December 202 Introduction We want to set up and validate the processing of Alpgen events with Herwig++ in ATLAS,
More informationDue: since the calculation takes longer than before, we ll make it due on 02/05/2016, Friday
Homework 3 Due: since the calculation takes longer than before, we ll make it due on 02/05/2016, Friday Email to: jqian@caltech.edu Introduction In this assignment, you will be using a commercial periodic
More informationTable of Contents. Table of Contents Electronic structure of NiO with DFT+U. Introduction The electronic structure of NiO calculated with DFT
Table of Contents Table of Contents Electronic structure of NiO with DFT+U Introduction The electronic structure of NiO calculated with DFT Setting up the calculation Performing the calculation Analysing
More informationCalculating Vibrational Spectra from Molecular Dynamics
Calculating Vibrational Spectra from Molecular Dynamics A Simulating a Trajectory with Wannier Centers To calculate IR spectra from Molecular Dynamics, it is necessary to have dipole information for the
More informationResidual Stress analysis
Residual Stress analysis Informations: strains Macro elastic strain tensor (I kind) Crystal anisotropic strains (II kind) Fe Cu C Macro and micro stresses Applied macro stresses Residual Stress/Strain
More informationTwo-dimensional Phosphorus Carbide as Promising Anode Materials for Lithium-ion Batteries
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry A. This journal is The Royal Society of Chemistry 2018 Supplementary Material for Two-dimensional Phosphorus Carbide as Promising
More informationQUANTUM CHEMISTRY PROJECT 3: PARTS B AND C
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 6, 2017 QUANTUM CHEMISTRY PROJECT 3: PARTS B AND C PART B: POTENTIAL CURVE, SPECTROSCOPIC CONSTANTS, AND DISSOCIATION ENERGY OF DIATOMIC HYDROGEN (20
More informationCalculating anisotropic physical properties from texture data using the MTEX open source package
Calculating anisotropic physical properties from texture data using the MTEX open source package David Mainprice 1, Ralf Hielscher 2, Helmut Schaeben 3 February 10, 2011 1 Geosciences Montpellier UMR CNRS
More informationThe quasi-harmonic approximation (QHA)
The quasi-harmonic approximation (QHA) M. Palumbo 19/01/2017 Trieste, Italy Limitations of the harmonic approximation E tot(r I, u I )=E tot(r I )+ I,α E tot u u Iα + 1 2 E tot Iα 2 I,α u u Iα u Iα u Jβ
More informationIntroduction to Hartree-Fock calculations in Spartan
EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,
More informationSurface stress and relaxation in metals
J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S0953-8984(00)11386-4 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown
More informationFULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY
FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY C.A. Madu and B.N Onwuagba Department of Physics, Federal University of Technology Owerri, Nigeria
More informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationAvogadro with YAeHMOP Manual
Avogadro with YAeHMOP Manual Patrick Avery, Herbert Ludowieg, Jochen Autschbach, and Eva Zurek Department of Chemistry University at Buffalo Buffalo, NY 14260-3000, USA September 11, 2017 Contents 1 Installation
More informationFirst-principles thermoelasticity of bcc iron under pressure
First-principles thermoelasticity of bcc iron under pressure Xianwei Sha and R. E. Cohen Carnegie Institution of Washington, 5251 Broad Branch Road, NW, Washington, D.C. 20015, USA Received 17 May 2006;
More informationQuantum ESPRESSO. PWSCF: first steps
Quantum ESPRESSO PWSCF: first steps What can I learn in this tutorial? What can I learn in this tutorial? How to run PWscf (pw.x) in self-consistent mode for Silicon How to get the band structure of Silicon
More informationChapter: 22. Visualization: Making INPUT File and Processing of Output Results
Chapter: 22 Visualization: Making INPUT File and Processing of Output Results Keywords: visualization, input and output structure, molecular orbital, electron density. In the previous chapters, we have
More informationChapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.
Chapter - Macromechanical Analysis of a Lamina Exercise Set. The number of independent elastic constants in three dimensions are: Anisotropic Monoclinic 3 Orthotropic 9 Transversely Orthotropic 5 Isotropic.
More informationScientific Programming in C XIII. Shell programming
Scientific Programming in C XIII. Shell programming Susi Lehtola 11 December 2012 Introduction Often in scientific computing one needs to do simple tasks related to renaming of files file conversions unit
More informationN = Shear stress / Shear strain
UNIT - I 1. What is meant by factor of safety? [A/M-15] It is the ratio between ultimate stress to the working stress. Factor of safety = Ultimate stress Permissible stress 2. Define Resilience. [A/M-15]
More informationHigh Temperature High Pressure Properties of Silica From Quantum Monte Carlo
High Temperature High Pressure Properties of Silica From Quantum Monte Carlo K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation:
More informationElasticity of hcp nonmagnetic Fe under pressure
PHYSICAL REVIEW B 68, 054103 003 Elasticity of hcp nonmagnetic Fe under pressure S. L. Qiu Department of Physics, Alloy Research Center, Florida Atlantic University, Boca Raton, Florida 33431-0991, USA
More informationDFT calculation of pressure induced phase transition in silicon
DFT calculation of pressure induced phase transition in silicon Michael Scherbela scherbela@student.tugraz.at 2015-12-07 Contents 1 Introduction 2 2 Structure of the input file 2 3 Calculation of phase
More informationEquilibrium state of a metal slab and surface stress
PHYSICAL REVIEW B VOLUME 60, NUMBER 23 15 DECEMBER 1999-I Equilibrium state of a metal slab and surface stress P. M. Marcus IBM Research Division, T. J. Watson Research Center, Yorktown Heights, New York
More informationTHERMOPHYSICAL PROPERTIES OF THORIUM COMPOUNDS FROM FIRST PRINCIPLES
THERMOPHYSICAL PROPERTIES OF THORIUM COMPOUNDS FROM FIRST PRINCIPLES Vinayak Mishra a,* and Shashank Chaturvedi a a Computational Analysis Division, Bhabha Atomic Research Centre, Visakhapatnam 530012,
More informationSupporting Information
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry C. This journal is The Royal Society of Chemistry 2016 Supporting Information Metal-Free Half-Metallicity in a High Energy Phase
More informationBoris Mantisi. Contents. Context 2. How to install APoGe 3. How APoGe works 4. Input File 5. Crystal File 8. How to use APoGe 10
Boris Mantisi Contents Context 2 How to install APoGe 3 How APoGe works 4 Input File 5 Crystal File 8 How to use APoGe 10 Context Polycrystalline structure plays an important role in the macroscopic properties
More informationProject. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release
Project First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release Contents Units Model (A4, B4) o Geometry! Solid Bodies! Parts! Parts! Body Groups! Parts! Parts
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More informationDetermining elastic constants of graphene using density functional theory
UNIVERSITY OF GRAZ Determining elastic constants of graphene using density functional theory by Max Niederreiter 01510759 A thesis submitted in partial fulfillment for the degree of Bachelor Of Science
More informationLeaf Spring (Material, Contact, geometric nonlinearity)
00 Summary Summary Nonlinear Static Analysis - Unit: N, mm - Geometric model: Leaf Spring.x_t Leaf Spring (Material, Contact, geometric nonlinearity) Nonlinear Material configuration - Stress - Strain
More informationElectronic structure calculations results from LDA+U method
Electronic structure calculations results from LDA+U method Vladimir I. Anisimov Institute of Metal Physics Ekaterinburg, Russia LDA+U method applications Mott insulators Polarons and stripes in cuprates
More informationStructural and Optical Properties of ZnSe under Pressure
www.stmjournals.com Structural and Optical Properties of ZnSe under Pressure A. Asad, A. Afaq* Center of Excellence in Solid State Physics, University of the Punjab Lahore-54590, Pakistan Abstract The
More informationTHEORETICAL STUDY OF THE STRUCTURAL, ELASTIC AND EQUATION OF STATE OF CHALCOPYRITE STRUCTURE AgAlS 2
Chalcogenide Letters Vol. 14, No. 7, July 2017, p. 251-257 THEORETICAL STUDY OF THE STRUCTURAL, ELASTIC AND EQUATION OF STATE OF CHALCOPYRITE STRUCTURE AgAlS 2 H. J. HOU a,*, T. J. LI b, G. CHENG a, X.
More informationBasics of Electronic-Structure Theory
A guideline through the tutorial 23.06.2009 For each exercise in this tutorial, a separate directory should be created in your home directory in which the input files are generated and the calculation
More informationTransition states and reaction paths
Transition states and reaction paths Lab 4 Theoretical background Transition state A transition structure is the molecular configuration that separates reactants and products. In a system with a single
More informationComputational Workshop: Running R-matrix Codes for Atomic Processes
Computational Workshop: Running R-matrix Codes for Atomic Processes - Prof. Sultana N. Nahar Astronomy, Ohio State U, Columbus, Ohio, USA Emails, Web: nahar.1@osu.edu http://www.astronomy.ohio-state.edu/
More information