dc_dft
Densitycorrected DFT.
This command runs a densitycorrected DFT calculation, in which an SCF is first performed at a KohnSham DFT or HartreeFock level, specified by the first dft subcommand, then followed by a singleshot energy evaluation at another DFT level, specified by the second dft subcommand, using the density obtained from the first SCF calculation.
For example, a DCDFT calculation using HF for obtaining the converged density and PBE for evaluating the energy might be run using:
dc_dft(
structure( molecule = water )
dft(
xc = HF
ao = '631G*'
)
dft(
xc = PBE
ao = '631G*'
)
)
Subcommands
Options
ansatz

SCF wavefunction ansatz.
 The type is string
 The default is default
 The value must be one of:
default
 default ansatz for given number of electronsr
 spinrestricted closedshell ansatzu
 spinunrestricted openshell ansatzro
 spinrestricted openshell ansatz. By default, the Davidson coupling scheme is used. See rohf_scheme for different choices of the coupling schemes.cu
constrained spinunrestricted openshell ansatz of Tsuchimochi and Scuseria.
This ansatz gives the same total energy and wavefunction as the restricted openshell ansatz ro, but exhibits more robust SCF convergence and gives physically meaningful orbital energies, i.e. satisfying Koopmans' theorem, and the aufbau principle.
Another difference between this ansatz and the ro ansatz is that it yields two different sets of spin orbitals, where as the ro ansatz gives only one set of orbitals for the alpha and beta spins.
However, one can obtain the ro orbitals from the cu wavefunction by setting ro_orbitals_from_cu to true. For example:
hf( structure( formula = OH bond_length = 1.0 angstrom ) ao = 'Def2SVP' ansatz = cu ro_orbitals_from_cu = true )
charge

Total molecular charge. No constraint is imposed, but clearly the value has to be such that for a given molecule the number of electrons is nonnegative, and that the number of electrons does not exceed the number that can be described with the given atomicorbital basis. Note that noninteger values are possible as well.
This option is deprecated. Please use the charge option inside structure instead.
This option is deprecated.
 The type is real
 There is no default value.
check_fock_norm

Check the norm of the Fock matrix against the numerical derivative of the energy before starting SCF. Only for debugging purposes.
 The type is bool
 The default is false
diis

Method of DIIS.
 The type is stringlowered
 The default is adiis+cdiis
 The value must be one of:
none
 DIIS is not used.cdiis
 This is the original DIIS algorithm by P. Pulay. P. Pulay, Chem. Phys. Lett. 73, 393 (1980).ediis
 Energybased DIIS which minimizes the HartreeFock energy functional. Kudin, Scuseria and Cances, J. Chem. Phys. 116, 8255 (2002).adiis
 DIIS based on the augmented RoothaanHall energy. X. Hu and W. Yang, J. Chem. Phys. 132, 054109 (2010).ediis+cdiis
 Combination of these two algorithms; see diis_switch_threshold_gradadiis+cdiis
 Combination of these two algorithms; see diis_switch_threshold_grad
diis_max

Maximum number of vectors for DIIS.
 The type is int
 The default is 10
 The value must be nonnegative
diis_shift

Diagonal shift for regularization of DIIS.
 The type is real
 The default is 0
 The value must be nonnegative
diis_switch_threshold_grad

In DIIS modes ADIIS+CDIIS and EDIIS+CDIIS the algorithm is switched to CDIIS if the gradient falls below this threshold.
 The type is real
 The default is 1e1
 The value must be nonnegative
energy_threshold

SCF convergence threshold using the absolute value of the energy difference between iterations.
 The type is real
 The default is 1e6
 The value must be nonnegative
guess

Initial guess for the SCF wavefunction.
 The type is string
 The default is SAD
 The value must be one of:
SAD
 Superposition of Atomic Densities (SAD).H0
 Using eigenvectors of the core Hamiltonian as the initial guess. This is equivalent to setting the initial density to zero.
level_shift

Turn on level shifting to improve SCF convergence. Raises the energy of virtual orbitals by level shifting value. Whether to remove the level shift at the end is controlled by option remove_level_shift.
 The type is real
 The default is 0
 The value must be nonnegative
load

Name of assigned result (containing orbitals) to start the SCF from.
 The type is string
 There is no default value.
max_iter

Maximum number of SCF iterations.
 The type is int
 The default is 128
 The value must be nonnegative
multiplicity

Spin multiplicity \(2S + 1\)  only valid for integerelectron calculations. For nonintegerelectron number, or for higher multiplcities, use the spin option.
By default entos selects the lowest available spin multiplicity.
 The type is string
 There is no default value.
 The value must be one of:
singlet
 2S + 1 = 1; S = 0doublet
 2S + 1 = 2; S = 1/2triplet
 2S + 1 = 3; S = 1quartet
 2S + 1 = 4; S = 3/2quintet
 2S + 1 = 5; S = 2sextet
 2S + 1 = 6; S = 5/2
name

Specify the name of a set of results.
This option is deprecated.
 The type is string
 There is no default value.
no_fail

If set to true, the program will continue its execution even if the SCF does not converge. This is useful for frozenorbital calculations.
 The type is bool
 The default is false
orbital_grad_threshold

SCF convergence threshold using the infinity norm of the orbital gradient.
 The type is real
 The default is 1e5
 The value must be nonnegative
orthogonalization

Method to orthogonalize the atomic orbital basis. Note that in case of a nearly linearly dependent basis (see overlap_threshold) the canonical algorithm is used irrespective of this option.
 The type is stringlowered
 The default is symmetric
 The value must be one of:
symmetric
 Symmetric orthogonalization of the AO basis. Uses the symmetric inverse square root of the overlap matrix for the orthogonalization.canonical
 Canonical orthogonalization of the AO basis.cholesky
 Orthogonalization through Cholesky decomposition of the overlap matrix.
overlap_threshold

Threshold for throwing away the smallest eigenvalues of the overlap matrix. This is a safety measure for calculations with nearlinearly dependent basis sets. Usually, the default value should be fine.
 The type is real
 The default is 1e7
 The value must be nonnegative
print_level

Print level.
 The type is int
 There is no default value.
 The value must be one of:
2
 No output1
 Minimum output0
 Output that doesn't scale with system size1
 Output that scales linearly with system size2
 (Debugging) output that scales quadratically with system size3
 (Debugging) output that scales cubically with system size
pseudo_diagonalization

Whether to use pseudodiagonalization instead of full diagonalization of the Fock matrix. Recommended for calculations on very large systems or with a very large basis.
 The type is bool
 The default is false
remove_level_shift

Remove level shift from final orbitals.
 The type is bool
 The default is true
ro_orbitals_from_cu

Whether to constuct ROHF/ROKS orbitals and orbital energies from a converged CUHF/CUKS wavefunction.
 The type is bool
 The default is true
rohf_scheme

The coupling scheme for constructing the ROHF effective Fock matrix. Note that different choices give the same SCF total energy, but they give different sets of orbitals and orbital energies. See Plakhutin, Gorelik, and Breslavskaya, J. Chem. Phys. 125, 204110 (2006) for a comprehensive discussion on different coupling schemes.
 The type is stringlowered
 The default is davidson
 The value must be one of:
mcweeny_diercksen
 Use the couping scheme of McWeeny and Diercksen. See McWeeny and Diercksen, J. Chem. Phys. 49, 4852 (1968)davidson
 Use the couping scheme of Davidson. See Davidson, Chem. Phys. Lett. 21, 565 (1973)
spin

For spinunrestricted calculations: spin is \(2S\), where S is the total spin angular momentum
Notice that the value of this option need not be an integer. See the documentation on
spin_z
for an explanation.By default entos selects the lowest available spin multiplicity.
 The type is real
 There is no default value.
 The value must be nonnegative
spin_z

For spinunrestricted calculations, \(2m_s\), the number of alphaspin electrons minus the number of betaspin electrons.
If omitted,
spin_z
defaults to be the same as the total spin angular momentum.Notice that the value of this option need not be an integer, and need not be the same as the value of
spin
. See this paper for an explanation of noninteger charge.The thinking is similar to thinking about noninteger charge. A fractionalcharge (fc) density matrix is a mixed state of the density matrices of two Slater determinants (each with integer charge)
\[ \rho_\text{fc}(N+1/2) = \rho(N)/2 + \rho(N+1)/2 . \]A fractionalspin (fs) density matrix is a mixed state of two density matrices with different numbers of alpha and beta electrons
\[ \rho_\text{fs}(S=S, m_s=0) = \rho(N_\alpha = N_\text{el}/2 + S, N_\beta = N_\text{el}/2  S )/2 + \rho(N_\alpha = N_\text{el}/2  S, N_\beta = N_\text{el}/2 + S )/2 . \] The type is real
 There is no default value.
temperature

Specify the temperature for the electrons and perform the SCF calculation in the canonicalensemble. Convergence is often poor with the default convergence settings, and using the original Pulay DIIS method is recommended (see diis option).
 The type is quantity
 The default is 0 kelvin