Tutorial -- Electric Field Gradient (EFG) in Periodic Solids

Performing a parallel relativistic calculation of solid-state EFG tensor is a simple two-step process

respect --pscf --inp=my-input-file --nt=4 --scratch=/path/to/scratch/directory
respect --pefg --inp=my-input-file --nt=4 --scratch=/path/to/scratch/directory

The first call of respect with the argument --pscf performs the periodic self-consistent field (pSCF) optimization of ground-state crystal orbitals, which is a necessary prerequisite for the following periodic EFG (pEFG) calculation. The second call of respect with the argument --pefg actually performs the periodic solid-state EFG simulation.

The other arguments mandatory to respect are

  • --inp
    specifies a name of the input file;

  • --nt
    specifies the number of physical cores required for parallel execution;

  • --scratch
    specifies a path to the scratch directory.

Note that the EFG tensor is related to the Nuclear Electric Quadrupole Hyperfine Interaction tensor (Q-tensor). Q-tensor elements (in MHz) equal the EFG tensor elements (in a.u.) times 234.9647 times the nuclear quadrupole moment (NQM in barn units, 1 barn = \(10^{-28}m^{2}\) ) and divided by \(2I(2I-1)\), where \(I\) is the nuclear spin. The Nuclear Quadrupole Coupling Constant (NQCC) (in MHz) is the largest value of the principal values of the EFG (in a.u.) times 234.9647 times the nuclear quadrupole moment (in barn units).

The input file my-input-file.inp should contain the input block pscf: with some pSCF-specific keywords, and the input block pefg: with some pEFG-specific keywords. A comprehensive list of all keywords can be found for pSCF here and for pEFG-tensor here. A simple example of my-input-file.inp looks like

#periodic SCF input block
#4c Dirac-Kohn-Sham DFT (PBE) calculation of 1D LiH system 
pscf:

   geo-units:  bohr
   geometry:
        Li   0.000000   0.000000   0.000000
        H    3.925886   0.000000   0.000000

   kmesh:   31x1x1
   lattice:
        a1: 7.851772

   cscale:        0.1      #artificially reduced speed of light
   method:        mdks/pbe
   basis:         pob-tzvp

   maxiterations: 30
   convergence:   1.0e-7
   nc-model:      gauss    #nuclear charge model
   initmo:        atomic

   peri:
        threshold: 1.0e-10


#periodic EFG input block
pefg:
     active-atoms: 1    #eval. EFG tensor for the first atom (Li)

As a final note, there are several important and worth-to-remember aspects associated with the input syntax, namely

  • the input is case-insensitive
    This means that the program does not distinguish between uppercase and lowercase letters.

  • the input is insensitive to the number of blank lines and/or comment lines
    All comments begin with the number sign (#), can start anywhere on a line and continue until the end of the line.

  • the input is compliant with the dictionary syntax of the YAML markup language
    This means that each input line is represented either by a single block: statement or by a simple keyword:value pair, such as

  • block1:
        keyword1: value1 
        keyword2: value2  
        ...
        block2: 
          keyword3: value3
          keyword4: value4
          ...
        block3: 
               keyword5: value5
               keyword6: value6
               ...
    

    It is essential to remember that all members of one block: are lines beginning at the same indentation level. Whitespace indentation is used to denote the block structure; however, tab characters are never allowed as indentation. The only exception to the YAML-based input syntax is the block geometry: which utilizes a simple xyz format for the molecular geometry specification.


TIPS & TRICKS

Q: Is there a way to launch periodic SCF and EFG calculations without the need to explicitly setup the scratch path by "--scratch=/path/to/scratch/directory"?

Yes, the argument "--scratch=/path/to/scratch/directory" can be saved to the file .respectrc in your home directory. If both the file and the command line argument exist, then ReSpect takes the scratch directory setting from the command line.


Latest Posts


Useful Links

Our Contacts

Hylleraas Centre
Department of Chemistry
UiT The Arctic University of Norway
Tromsø, NO-9037 Norway
Email: info@respectprogram.eu