welcome: please sign in
location: Diff for "xuanyuan"
Differences between revisions 22 and 23
⇤ ← Revision 22 as of 2018-12-27 03:10:42 →
Size: 4598
Editor: bsuo
Comment:
← Revision 23 as of 2018-12-27 05:46:44 → ⇥
Size: 4485
Editor: lzd
Comment:
Deletions are marked like this. Additions are marked like this.
Line 7: Line 7:
Line 9: Line 8:

Line 16: Line 13:

Line 21: Line 16:
Ask for integral direct calculations.  Ask for integral direct calculations.
Line 23: Line 18:
Line 29: Line 23:
Line 35: Line 30:

Line 41: Line 34:
Examples:
Line 42: Line 36:
Examples:
Line 46: Line 39:
  512MW    512MW
Line 49: Line 42:
Line 52: Line 44:
  Range separation ERIs required.    Range separation ERIs required.
Line 55: Line 47:
Examples:
Line 56: Line 49:
Examples:
Line 63: Line 55:
Line 70: Line 61:
Examples:
Line 71: Line 63:
Examples:
Line 79: Line 70:
Line 82: Line 72:
Socint is a keyword to turn on soc integral calculations in post-SCF steps.
Default option for hsoc is 0 (only 1e-soc int). 		  is a keyword to turn on soc integral calculations in post-SCF steps.
Default option for hsoc is 0 (only 1e-soc int).
Line 85: Line 75:
Other options are used in soint_util/somf2e.F90 for choosing different combinations of so1e and so2e operators. Other options are used in soint_util/somf2e.F90 for choosing different combinations of so1e and mean-feild so2e (SOMF) operators.
Line 90: Line 80:

2 so-1e + somf-1c (one-center approximation)
Line 103: Line 90:
}}}
Examples:
Line 104: Line 93:
}}}

Examples:
Line 117: Line 103:
Line 126: Line 111:
Line 135: Line 119:
Examples:
Line 136: Line 121:
Examples:
Line 143: Line 127:
Line 145: Line 128:
Line 152: Line 134:
Line 157: Line 138:
Line 162: Line 142:
Line 167: Line 146:
= Depend Files =
||Filename ||Description ||Format ||
|| || || ||
Line 168: Line 150:
= Depend Files =

|| Filename || Description || Format ||
|| || || ||

Xuanyuan

Contents

  1. Xuanyuan
    1. General keywords
      1. Direct
      2. Schwarz
      3. Maxmem
      4. RS
      5. Scalar & Heff
      6. Soint & Hsoc
      7. Nuclear & Inuc
      8. Cholesky
    2. Expert keywords
      1. NoCheck
      2. NRDebug
      3. Keyword3
      4. Keyword4
  2. Depend Files
  3. Examples
  4. References

Xuanyuan is used to calculate one electron and two electron integrals. It is named after Chinese ancestor Xuanyuan Huangdi.

General keywords

Direct

Ask for integral direct calculations.

Schwarz

Used with direct, ask for Schwarz equality prescreening.

Examples: 

$xuanyuan
Direct
Schwarz
$end

Maxmem

Set maximum memory used in integral calculation. Unit can be MW and GW, i.e. Mega Words and Giga Words

Examples: 

$xuanyuan
Maxmem
  512MW
$end

RS

  • Range separation ERIs required. No default value. Suggested value: 0.33.

Examples: 

$xuanyuan
Rs
 0.33
$end

Scalar & Heff

Scalar is a keyword to turn on scalar relativistic effects using sf-X2C (Heff=3) by default.

Other options for Heff are (0, nonrelativistic; 1, sf-ZORA; 2, sf-IORA; 3/4, sf-X2C; 5, sf-X2C+so-DKH3 (spin-free); 21, sf-X2C with gradients)

Examples: 

$xuanyuan
scalar
heff
3
$end

Soint & Hsoc

  • is a keyword to turn on soc integral calculations in post-SCF steps.

Default option for hsoc is 0 (only 1e-soc int).

Other options are used in soint_util/somf2e.F90 for choosing different combinations of so1e and mean-feild so2e (SOMF) operators.

0 so-1e

1 so-1e + somf (two-electron spin-orbit interaction is included via an effective fock operator) 3 so-1e + somf-1c / no soo (turn off spin-other-orbit contributions)

4 so-1e + somf-1c / no soo + WSO_XC (use dft xc functional as soo part)

5 so-1e + somf-1c / no soo + WSO_XC(-2x: following Neese's paper scale dft part by -2 to mimic soo part)

These options plus 10 gives the operators in BP approximations. In practice, hsoc=1 is the most accurate, and hsoc=2 is preferred for large molecules.

Note if heff=5, then the one-electron part will be calculated in xuanyuan and stored in disc for so-DKH3 type one-electron spin-orbit term. The accuracy of such operator requires further tests.

Examples: 

$xuanyuan
scalar
heff
3
soint
hsoc
2
$end

Nuclear & Inuc

Inuc defines the nuclear charge distribution used in the V and pVp integrals, which can be -1 for point charge model (debug only), 0 for point charge model (default), 1 for finite nucleus model by an s-type Gaussian function, and other finite nucleus models (N.Y.I.).

For Za < 110, the nuclear charge radii are taken from Ref.[Visscher1997] (in a.u).

For Za ≥ 110, the nuclear charge radius is 0.57 + 0.836 * A1/3 (in fm), where the isotope mass number A is estimated by Za according to the relationship A(Za) = 0.004467 * Za2 + 2.163 * Za - 1.168. See Appendix A in Ref.[Andrae2000] and Ref.[Andrae2002].

NOTE: the finite nucleus model has been implemented only in scalar calculations at present, but will be used in SOC calculations soon.

Cholesky

  • The following line contains a string and a float number. Set method and threshold of ERI Cholesky decomposition. S-CD for standard CD. 1c-CD for one center Cholesky decomposition.

Examples: 

$xuanyuan
Cholesky
  S-CD 1.d-5
$end

Expert keywords

NoCheck

For Heff=21 only: check inverse variational collapse (IVC; see Ref.[Liu2007]). Stop (0; default) or not (1) in the case of IVC.

IVC may lead to numerical instability, which may be serious in geometry optimization.

NRDebug

In relativistic calculations, use a C-light of 10^8 to reproduce non-relativistic results (for debug only).

Keyword3

xxx

Keyword4

xxx

Depend Files

Filename

Description

Format

Examples

N.A.

References

  • [Andrae2000] D. Andrae, Phys. Rep. 336, 414 (2000).
  • [Andrae2002] D. Andrae, Nuclear charge density distributions in quantum chemistry, in Relativistic Electronic Structure Theory, Part 1: Fundamentals, P. Schwerdtfeger Ed., Theoretical and Computational Chemistry, Vol. 11, Elsevier, 2002.
  • [Liu2007] W. Liu and W. Kutzelnigg, J. Chem. Phys. 126, 114107 (2007).
  • [Visscher1997] L. Visscher and K. G. Dyall, At. Data and Nucl. Data Tables 67, 207 (1997).

xuanyuan (last edited 2024-01-31 09:23:43 by wzou)