This script compiles with python 3 or higher version.
ENERGY Vs STRAIN APPROACH: Two scripts are provided for different deformation schemes. Please read the documentation as well as the details on the screen.
To generate input files for VASP. It uses the lagrangian strain to strain the system with user specified deformation constant. This value should be in the harmonic regime otherwise this approximation will fail. Large deformation leads to phase transitions.
NB: This code is the modification of the script provided with the exciting code. All rights belongs to the original author, if there is an error in the code or bug please contact the exciting team or me. Two approaches differes in the way deformation atrix is defined:
**Deformation 1:
I have modified the script to read CONTCAR optimized file (already minimized with IBRION=2, ISIF=3)
and then print out the deformed POSCAR files with various strain.
Just follow the recipe on the screen and it will guide you to the rest of the process.
For more information please visit exciting website (Elastic code). The literature behind can be found at the exciting website.
|ฮท[0] ฮท[5]/2 ฮท[4]/2|
ฮท = |ฮท[5]/2 ฮท[1] ฮท[3]/2|
|ฮท[4]/2 ฮท[3]/2 ฮท[2] |
D' = I + ฮท")
**Deformation 2:
Cubic_strain.py script uses another deformation matrix technique. Please note it only pertians to
cubic system. It can be extended to include other crystal deformation matrix.
0 => (ฮท, ฮท, ฮท,0,0,0) | volumetric strain [3C11+6C12=9B0]
1 => (ฮท,-ฮท,1/1-ฮท**2 - 1,0,0,0) | volume-conserving orthorhombic distortion 2[C11-C12]
2 => (0, 0,1/1-ฮท**2 - 1,0,0,2ฮท)| volume-conserving monoclinic distortion [2C44]
[1] ref: exciting
[2] ref: materialsproject
**Another interesting appraoch is: STRESS Vs STRAIN **
[3] ref: ELASTIC
[L.D. Landau, E.M. Lifszyc, Theory of elasticity, Elsevier (1986) ; ISBN: 075062633X, 9780750626330]
#============================
Deform Crystal structure for Monoclinic structure. It reads only POSCAR file.
USAGE::
Launch the code and it will prompt you to enter the structural details. Remember to rad the code and understand what it is doing and if you find some bug please let me know.
NB: No guarranty is implied before using for your work ot publication please do check with other codes. Exciting and Elastic code on my github are distributed freely.
๐ข๐ป๐ฒ ๐ณ๐ถ๐ฒ๐น๐ฑ ๐ผ๐ณ ๐๐ผ๐ฟ๐ธ ๐ถ๐ป ๐๐ต๐ถ๐ฐ๐ต ๐๐ต๐ฒ๐ฟ๐ฒ ๐ต๐ฎ๐ ๐ฏ๐ฒ๐ฒ๐ป ๐๐ผ๐ผ ๐บ๐๐ฐ๐ต ๐๐ฝ๐ฒ๐ฐ๐๐น๐ฎ๐๐ถ๐ผ๐ป ๐ถ๐ ๐ฐ๐ผ๐๐บ๐ผ๐น๐ผ๐ด๐. ๐ง๐ต๐ฒ๐ฟ๐ฒ ๐ฎ๐ฟ๐ฒ ๐๐ฒ๐ฟ๐ ๐ณ๐ฒ๐ ๐ต๐ฎ๐ฟ๐ฑ ๐ณ๐ฎ๐ฐ๐๐ ๐๐ผ ๐ด๐ผ ๐ผ๐ป, ๐ฏ๐๐ ๐๐ต๐ฒ๐ผ๐ฟ๐ฒ๐๐ถ๐ฐ๐ฎ๐น ๐๐ผ๐ฟ๐ธ๐ฒ๐ฟ๐ ๐ต๐ฎ๐๐ฒ ๐ฏ๐ฒ๐ฒ๐ป ๐ฏ๐๐๐ ๐ฐ๐ผ๐ป๐๐๐ฟ๐๐ฐ๐๐ถ๐ป๐ด ๐๐ฎ๐ฟ๐ถ๐ผ๐๐ ๐บ๐ผ๐ฑ๐ฒ๐น๐ ๐ณ๐ผ๐ฟ ๐๐ต๐ฒ ๐จ๐ป๐ถ๐๐ฒ๐ฟ๐๐ฒ, ๐ฏ๐ฎ๐๐ฒ๐ฑ ๐ผ๐ป ๐ฎ๐ป๐ ๐ฎ๐๐๐๐บ๐ฝ๐๐ถ๐ผ๐ป๐ ๐๐ต๐ฎ๐ ๐๐ต๐ฒ๐ ๐ณ๐ฎ๐ป๐ฐ๐. ๐ง๐ต๐ฒ๐๐ฒ ๐บ๐ผ๐ฑ๐ฒ๐น๐ ๐ฎ๐ฟ๐ฒ ๐ฝ๐ฟ๐ผ๐ฏ๐ฎ๐ฏ๐น๐ ๐ฎ๐น๐น ๐๐ฟ๐ผ๐ป๐ด. ๐๐ ๐ถ๐ ๐๐๐๐ฎ๐น๐น๐ ๐ฎ๐๐๐๐บ๐ฒ๐ฑ ๐๐ต๐ฎ๐ ๐๐ต๐ฒ ๐น๐ฎ๐๐ ๐ผ๐ณ ๐ป๐ฎ๐๐๐ฟ๐ฒ ๐ต๐ฎ๐๐ฒ ๐ฎ๐น๐๐ฎ๐๐ ๐ฏ๐ฒ๐ฒ๐ป ๐๐ต๐ฒ ๐๐ฎ๐บ๐ฒ ๐ฎ๐ ๐๐ต๐ฒ๐ ๐ฎ๐ฟ๐ฒ ๐ป๐ผ๐. ๐ง๐ต๐ฒ๐ฟ๐ฒ ๐ถ๐ ๐ป๐ผ ๐ท๐๐๐๐ถ๐ณ๐ถ๐ฐ๐ฎ๐๐ถ๐ผ๐ป ๐ณ๐ผ๐ฟ ๐๐ต๐ถ๐. ๐ง๐ต๐ฒ ๐น๐ฎ๐๐ ๐บ๐ฎ๐ ๐ฏ๐ฒ ๐ฐ๐ต๐ฎ๐ป๐ด๐ถ๐ป๐ด, ๐ฎ๐ป๐ฑ ๐ถ๐ป ๐ฝ๐ฎ๐ฟ๐๐ถ๐ฐ๐๐น๐ฎ๐ฟ, ๐พ๐๐ฎ๐ป๐๐ถ๐๐ถ๐ฒ๐ ๐๐ต๐ฎ๐ ๐ฎ๐ฟ๐ฒ ๐ฐ๐ผ๐ป๐๐ถ๐ฑ๐ฒ๐ฟ๐ฒ๐ฑ ๐๐ผ ๐ฏ๐ฒ ๐ฐ๐ผ๐ป๐๐๐ฎ๐ป๐๐ ๐ผ๐ณ ๐ป๐ฎ๐๐๐ฟ๐ฒ ๐บ๐ฎ๐ ๐ฏ๐ฒ ๐๐ฎ๐ฟ๐๐ถ๐ป๐ด ๐๐ถ๐๐ต ๐ฐ๐ผ๐๐บ๐ผ๐น๐ผ๐ด๐ถ๐ฐ๐ฎ๐น ๐๐ถ๐บ๐ฒ. ๐ฆ๐๐ฐ๐ต ๐๐ฎ๐ฟ๐ถ๐ฎ๐๐ถ๐ผ๐ป๐ ๐๐ผ๐๐น๐ฑ ๐ฐ๐ผ๐บ๐ฝ๐น๐ฒ๐๐ฒ๐น๐ ๐๐ฝ๐๐ฒ๐ ๐๐ต๐ฒ ๐บ๐ผ๐ฑ๐ฒ๐น ๐บ๐ฎ๐ธ๐ฒ๐ฟ๐." ๐๐ถ๐ฟ๐ฎ๐ฐ, ๐ฃ๐ฎ๐๐น. ๐ข๐ป ๐บ๐ฒ๐๐ต๐ผ๐ฑ๐ ๐ถ๐ป ๐๐ต๐ฒ๐ผ๐ฟ๐ฒ๐๐ถ๐ฐ๐ฎ๐น ๐ฝ๐ต๐๐๐ถ๐ฐ๐. (๐ง๐ฟ๐ถ๐ฒ๐๐๐ฒ. ๐๐๐ป๐ฒ ๐ญ๐ต๐ฒ๐ด .)