Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
forum:data:2023:directional_x-ray_spectroscopy_potential_bug [2023/04/08 00:10] – Charles Cardot | forum:data:2023:directional_x-ray_spectroscopy_potential_bug [2023/04/08 01:00] (current) – Charles Cardot | ||
---|---|---|---|
Line 9: | Line 9: | ||
I have been experimenting with calculating direction dependent x-ray emission using Quanty and the DFT+MLFT framework. However, I seem to be running into an issue with the mean field operators calculated from the Tight Binding (TB) object. | I have been experimenting with calculating direction dependent x-ray emission using Quanty and the DFT+MLFT framework. However, I seem to be running into an issue with the mean field operators calculated from the Tight Binding (TB) object. | ||
- | The directional dependence is encoded in the DFT step, where I distort the crystal lattice slightly (~1%) along a particular axis. When I stretch along the a-axis I get an tight binding Hamiltonian, | + | The directional dependence is encoded in the DFT step, where I distort the crystal lattice slightly (~2%) along a particular axis. |
- | a_stretched_HDFT: | + | # lattice constants; |
+ | @l@ 4.27633759 4.18633759 4.18633759 | ||
+ | # set axis angles | ||
+ | @a@ 90. 90. 90. | ||
+ | # setup Wyckoff positions | ||
+ | | ||
+ | # Now, give list of ALL !!! Wyckoff positions. | ||
+ | @1@ Ni @ 0.00000000 | ||
+ | @2@ O @ 0.00000000 | ||
+ | |||
+ | |||
+ | When I stretch along the a-axis and go through the regular Wannierization process to create a tight binding Hamiltonian, | ||
+ | |||
+ | a_stretched_HDFT: | ||
which makes sense in the context of the distorted octahedral crystal field. The //xy// and //xz// orbitals are still degenerate, while the //yz// is now lower energy. Also the //z^2-r^2// and //x^2-y^2// orbitals become mixed. The MLFT calculation within Quanty proceeds much the same way as described in the 2022 and 2019 Heidelberg tutorials for DFT + MLFT, where the mean field components Coulomb interaction is subtracted from HDFT to avoid double counting. The full code that I am running can be found here (https:// | which makes sense in the context of the distorted octahedral crystal field. The //xy// and //xz// orbitals are still degenerate, while the //yz// is now lower energy. Also the //z^2-r^2// and //x^2-y^2// orbitals become mixed. The MLFT calculation within Quanty proceeds much the same way as described in the 2022 and 2019 Heidelberg tutorials for DFT + MLFT, where the mean field components Coulomb interaction is subtracted from HDFT to avoid double counting. The full code that I am running can be found here (https:// |