CRYSTAL & phono3py calculation¶
CRYSTAL program package has a robust built-in phonon calculation workflow for harmonic phonon properties. However, combining CRYSTAL with Phono3py enables convenient access to anharmonic phonon properties.
An example for CRYSTAL is found in the
To invoke the CRYSTAL interface,
--crystal option has to be always
% phono3py --crystal [options] [arguments]
When the file name of the unit cell is different from the default one
(see Default unit cell file name),
is used to specify the file name. CRYSTAL unit cell file parser used in
phono3py is the same as that in phonopy. It reads a limited number of
keywords that are documented in the phonopy web site
In this example (Si-CRYSTAL), the CRYSTAL output file is crystal.o. This is the default file name for the CRYSTAL interface, so the -c crystal.o parameter is not needed.
Create supercells with displacements (4x4x4 for 2nd order FC, 2x2x2 for 3rd order FC)
% phono3py --crystal --dim="2 2 2" --dim-fc2="4 4 4" -d
57 supercell files (
supercell-xxx.d12/.ext) for the third order force constants are created. In addition, one supercell file (
supercell_fc2-00001.d12/.ext) is created for the second order force constants.
To make valid CRYSTAL input files, there are two possible options:
Manually: modify the generated supercell-xxx.d12 and supercell_fc2-xxxxx.d12 files by replacing the line
Insert basis sets and parameters herewith the basis set and computational parameters.
Recommended option: before generating the supercells, include files named
TEMPLATE3in the current directory. These files should contain the basis sets and computational parameters for CRYSTAL (see the example). When phono3py finds these files it automatically generates complete CRYSTAL input files in the step 1.
Note that supercells with displacements must not be relaxed in the force calculations, because atomic forces induced by a small atomic displacement are what we need for phonon calculation. To get accurate forces, TOLDEE parameter should be 10 or higher. Phono3py includes this parameter and the necessary GRADCAL keyword automatically in the inputs.
Then, CRYSTAL supercell calculations are executed to obtain forces on atoms, e.g., as follows:
% runcry17 supercell-001.d12
Collect forces in
% phono3py --crystal --cf3 supercell-*o % phono3py --crystal --cf2 supercell_fc2-*o
disp_fc2.yamlare used to create
FORCES_FC2, therefore they must exist in current directory.
Calculate 3rd and 2nd order force constants in files
% phono3py --crystal --dim="2 2 2" --dim-fc2="4 4 4" --sym-fc
--sym-fcis used to symmetrize second- and third-order force constants.
Thermal conductivity calculation:
% phono3py --crystal --fc3 --fc2 --dim="2 2 2" --dim-fc2="4 4 4" --mesh="20 20 20" --br
--brinvokes the Relaxation Time Approximation. Add
--isotopefor isotope scattering. Check the effect of
--nacfor polar systems. Carefully test the convergence with respect to