Tip
All input files can be downloaded: Files
.
pseudopotential
This keyword defines the pseudopotentials used for quantum chemistry calculations. You can define pseudopotentials in several flexible ways.
Warning
pseudopotential
only contains pseudopotential and does not contain any valence electron basis set information. You must assign basis sets in basis
. The valence basis and core pseudopotential must match. See basis.
Below is a typical example. H is represented by a def2-TZVP basis set. For Au, the valence electrons are represented by a def2-TZVP basis set, and core electrons are repesented by def2-ecp, i.e. Stuggart-Cologne pseudopotential.
1basis
2 def2-TZVP
3end
4pseudopotential
5 def2-ecp
6end
7# The Karlsruhe basis sets for Au was developed with Stuggart-Cologne pseudopotential.
8mol
9 Au 0. 0. 0. # SDD+def2-TZVP
10 H 0. 0. 1. # def2-TZVP
11end
Warning
If you write it in this WRONG way:
1# A wrong input!
2basis
3 def2-TZVP
4end
5mol
6 Au 0. 0. 0.
7 H 0. 0. 1.
8end
Then, NO pseudopotential is applied for any element!
Also, keep in mind that the valence basis and core pseudopotential must match. The following combinations are usually accepted:
Valence Basis Set |
Pseudopotential |
---|---|
def2-X |
def2-ecp |
(aug-)cc-X-pp |
cc-ecp |
lanlX |
lanl-ecp |
Using Built-in Pseudopotentials
A lot of important pseudopotentials have been provided in a folder pseudopotential
in the same path of Qbics. The files are named after their names well-known in computational chemistry community. For example, pseudopotential/def2-ecp
contains the Stuggart-Cologne pseudopotentials. All files are named in small cases.
To use them, simple write down the basis set name. It is case-insensitive. For example, to use def2-ecp:
1pseudopotential
2 def2-ecp
3end
Qbics will extract pseudopotential information from pseudopotential/def2-ecp
for all atoms that have pseudopotentials. For example, you molecule contains only C, H, N, Ce, and F. Since in pseudopotential/sdd
, there is only pseudopotential for Ce, then for C, H, N, and F, no pseudopotentials are applied.
Explicit Pseudopotential Definitions
You can also explicitly define your pseudopotentials. For example, you want to apply pseudopotentials for Rb and Sr, then their pseudopotentials can be defined in this way:
1pseudopotential
2 RB 0
3 RB-ECP 3 28
4 f POTENTIAL
5 1
6 2 3.84311400 -12.31690000
7 s-f POTENTIAL
8 3
9 2 5.03655100 89.50019800
10 2 1.97084900 0.49376100
11 2 3.84311400 12.31690000
12 p-f POTENTIAL
13 3
14 2 4.25834100 58.56897400
15 2 1.47070900 0.43179100
16 2 3.84311400 12.31690000
17 d-f POTENTIAL
18 3
19 2 3.02312700 26.22489800
20 2 0.65038300 0.96283900
21 2 3.84311400 12.31690000
22 ****
23 SR 0
24 SR-ECP 3 28
25 f POTENTIAL
26 1
27 2 4.63397500 -15.80599200
28 s-f POTENTIAL
29 3
30 2 7.40007400 135.47943000
31 2 3.60637900 17.53446300
32 2 4.63397500 15.80599200
33 p-f POTENTIAL
34 3
35 2 6.48486800 88.35970900
36 2 3.28805300 15.39437200
37 2 4.63397500 15.80599200
38 d-f POTENTIAL
39 3
40 2 4.62284100 29.88898700
41 2 2.24690400 6.65941400
42 2 4.63397500 15.80599200
43 ****
44end
The analyitcal expression of pseudopotential is:
The pseudopotential definition is of standard Gaussian94 format:
The definition of the pseudopotential for each atom ends with 4 asterisks, i.e.
****
.The definition starts with the element name like
Rb
and a0
. Currently0
has no meaning.Then, three parameters are given: pseudopotential name, maximum angular momentum \(L\), and number of core electrons.
Then, the semi-local part (\(V_l(r) (0\le l < L)\), defined by
s-f POTENTIAL
, etc.) and the local one (\(V_L(r)\), defined byf POTENTIAL
) are listed. Each pseudopotential has the form :The first line is a comment.
The contraction degree \(K\).
Then, each line defines the power \(n_{kl}\) , the exponent \(\xi_{kl}\), and the contraction coefficient \(d_{kl}\). They are 3 real numbers.
Hint
Pseudopotentials in Gaussian94 format can be obtained from:
But, remember to replace D
to E
since the former is not recognized by Qbics. Also, add ****
between elements.
Using Self-defined Pseudopotential Files
You can also put your explicit pseudopotential definitions into some files, say /home/zhang/userdef/my-own-pseudopotential
. Qbics will automatically read it if you give explicit file name including path.
1pseudopotential
2 /home/zhang/userdef/my-own-pseudopotential
3end
Theoretical Background
Pseudopotential is widely used in the fields of quantum chemistry, atomic physics, solid state physics, and computational materials science. The basic idea is to replace the influence of the nucleus and core electrons on the valence electrons with an effective potential. It can reduce the size of the basis set and also include relativistic effects on a non-relativistic level.
Input Examples
Some examples are also given in basis.
Example: Using Explicitly Defined Pseudopotentials for Ce(H2O)83+
Below is an example of calculating Ce(H2O)83+ using explicitly defined pseudopotentials for Ce and standard def2-svp for Ce, H, and O.
Note that, since the pseudopotential has been applied to describe 47 core electrons of Ce ([Kr]4d104f1), the apparant spin multiplicity is 1
, greatly simplifying the calculation.
1basis
2 def2-svp
3end
4
5pseudopotential
6 Ce 0
7 ECP47MWB 4 47
8 G-Komponente
9 1
10 2 1.000000 0.000000
11 S-G
12 2
13 2 3.522200 95.842155
14 2 1.761100 -3.775040
15 P-G
16 2
17 2 3.017700 68.092779
18 2 1.508900 -0.966756
19 D-G
20 2
21 2 2.144300 36.381848
22 2 1.072200 0.190447
23 F-G
24 1
25 2 4.278500 -40.585328
26 ****
27end
28
29grimmedisp
30 type bj
31end
32
33scf
34 charge +3
35 spin2p1 1
36end
37
38mol
39 Ce -0.0000088 0.0000037 -0.0006635
40 O 1.8174770 -1.1442674 -1.3565307
41 O 2.0940885 0.4762325 1.3551096
42 O 0.4762019 -2.0941013 1.3550928
43 O -1.1442759 -1.8174502 -1.3565612
44 O -1.8174735 1.1443092 -1.3565265
45 O -2.0941079 -0.4762528 1.3550846
46 O 1.1442883 1.8174686 -1.3565204
47 O -0.4762411 2.0940836 1.3551277
48 H -2.7384080 -1.1961061 1.2400375
49 H -2.4318208 0.0300465 2.1143243
50 H -1.6982474 1.7415082 -2.1154550
51 H -2.7821029 1.0904008 -1.2419254
52 H 0.0300785 2.4317909 2.1143563
53 H -1.1960900 2.7383918 1.2400978
54 H 2.7383924 1.1960849 1.2400788
55 H 2.4317951 -0.0300811 2.1143427
56 H 1.6982634 -1.7414504 -2.1154737
57 H 2.7821042 -1.0903611 -1.2419111
58 H -1.0903471 -2.7820826 -1.2419964
59 H -1.7414853 -1.6982065 -2.1154789
60 H 1.7415073 1.6982406 -2.1154327
61 H 1.0903589 2.7820987 -1.2419353
62 H 1.1960541 -2.7384032 1.2400470
63 H -0.0301159 -2.4318311 2.1143124
64end
65
66task
67 energy b3lyp
68end