WRT NMRCUR: Difference between revisions
No edit summary |
No edit summary |
||
| Line 5: | Line 5: | ||
---- | ---- | ||
{{not released}} | |||
In conjunction with {{TAGO|LCHIMAG|True}}, {{TAG|WRT_NMRCUR}} allows to write the current response on the fine grid {{TAG|NGXF}} x {{TAG|NGYF}} x {{TAG|NGZF}} in atomic units (hartree bohr<math>^{-2}</math>) to an external magnetic field within [[:Category:NMR|linear response NMR]]. The output is written to {{FILE|NMRCURBX}}, {{FILE|NMRCURBY}}, and/or {{FILE|NMRCURBZ}} depending on the selected direction of the perturbing <math>\mathbf{B}</math> field: | In conjunction with {{TAGO|LCHIMAG|True}}, {{TAG|WRT_NMRCUR}} allows to write the current response on the fine grid {{TAG|NGXF}} x {{TAG|NGYF}} x {{TAG|NGZF}} in atomic units (hartree bohr<math>^{-2}</math>) to an external magnetic field within [[:Category:NMR|linear response NMR]]. The output is written to {{FILE|NMRCURBX}}, {{FILE|NMRCURBY}}, and/or {{FILE|NMRCURBZ}} depending on the selected direction of the perturbing <math>\mathbf{B}</math> field: | ||
Revision as of 11:10, 7 March 2025
WRT_NMRCUR = 0 | 1 | 2 | 3 | 4
Default: WRT_NMRCUR = 0
Description: Allows to write the NMR current response in atomic units to file.
Warning: Not yet released!
This page contains information about a feature that will be available in a future release of VASP. In other words, currently you cannot use it even with the latest version of VASP. The information may change significantly until it is released.
In conjunction with LCHIMAG = True, WRT_NMRCUR allows to write the current response on the fine grid NGXF x NGYF x NGZF in atomic units (hartree bohr[math]\displaystyle{ ^{-2} }[/math]) to an external magnetic field within linear response NMR. The output is written to NMRCURBX, NMRCURBY, and/or NMRCURBZ depending on the selected direction of the perturbing [math]\displaystyle{ \mathbf{B} }[/math] field:
WRT_NMRCUR = 0: no current response written to file (default)WRT_NMRCUR = 1: [math]\displaystyle{ B_x }[/math]WRT_NMRCUR = 2: [math]\displaystyle{ B_y }[/math]WRT_NMRCUR = 3: [math]\displaystyle{ B_z }[/math]WRT_NMRCUR = 4: all three directions of [math]\displaystyle{ \mathbf{B}=(B_x,B_y,B_z)^T }[/math]
It is also written to vaspout.h5, if compiled with HDF5 support. You can find the data groups
/results/nmrcurbx Group
/results/nmrcurbx/grid Dataset {3}
/results/nmrcurbx/structure Group
/results/nmrcurbx/structure/position Group
/results/nmrcurbx/structure/position/direct_coordinates Dataset {SCALAR}
/results/nmrcurbx/structure/position/ion_sha256 Dataset {1}
/results/nmrcurbx/structure/position/ion_types Dataset {1}
/results/nmrcurbx/structure/position/lattice_vectors Dataset {3, 3}
/results/nmrcurbx/structure/position/number_ion_types Dataset {1}
/results/nmrcurbx/structure/position/position_ions Dataset {2, 3}
/results/nmrcurbx/structure/position/scale Dataset {SCALAR}
/results/nmrcurbx/structure/position/system Dataset {SCALAR}
/results/nmrcurbx/values Dataset {3, 24, 24, 24}
and use py4vasp to access these, e.g., using
import py4vasp as pv
calc = pv.Calculation.from_path(".")
py4vasp.current_density.plot("NMR(bx)")
Warning: For bulk calculations you must switch off the use of symmetry. In other words, set ISYM <= 0 if there is more than a single k point at zero (the Γ point).
|
Tip: Consider switching on current augmentation (LLRAUG = True).
|