:orphan:
.. _relccsd:
Kramers unrestricted Coupled-Cluster methods
============================================
See :cite:`Visscher1996` for initial CCSD(T) implementation. The Fock-space CC
implementation is described in :cite:`Visscher2001`.
The old namelist style input used up to DIRAC10 is now no longer supported.
The input below follows the regular DIRAC style conventions.
.. index:: **RELCC
.. _**RELCC:
=========
\*\*RELCC
=========
Specification of reference determinant, type of calculation, and general
settings.
.. index:: .FOCKSP
.. _RELCC_.FOCKSP:
.FOCKSP
-------
Activate the Fock space module. This option should be used for
multireference calculations. See further :ref:`*CCFSPC`.
Because the first sector of Fock space gives the same result as a
regular CCSD calculation, the latter calculation is switched off. If you
do wish to perform also a regular CC calculation (e.g. to get the CCSD(T) energy)
you need to activate this explicitly via .ENERGY (see below).
.. index:: .ENERGY
.. _RELCC_.ENERGY:
.ENERGY
-------
Activate the energy calculation. This is the preferred option for
calculations on closed shell or simple open shell systems and need not
be specified explictly in such cases. For Fock space calculations the
keyword switches on a separate single reference calculation done prior
to the FS calculation.
*Default:*
::
Perform energy calculation..
.. index:: .GRADIENT
.. _RELCC_.GRADIENT:
.GRADIENT
---------
Calculate the effective 1-particle density matrix (currently only for a
closed shell MP2 wave function). This option can be used to calculate
molecular properties. See :cite:`vanStralen2005`.
*Default:*
::
No gradient calculation.
.. index:: .NELEC
.. _RELCC_.NELEC:
.NELEC
------
Number of active electrons. This variable determines the reference determinant to be used in the exponential expansion of the wave function.
Since the default values correspond to the information passed on by the MOLTRA code on basis of the Hartree-Fock occupations
and chosen range of active orbitals in MOLTRA, **there is usually no need to specify this variable manually**.
Notable exception is the treatment of a high-spin open shell state with the single reference CCSD(T) ansatz.
In such a case one first needs to determine the number of closed shell electrons in each irrep and
then add to these occupations the open shell electrons (see corresponding tutorial examples).
*Arguments:*
::
Integer (NELEC(I),I=1,NFSYM*2).
*Default:*
::
Active electrons in these irreps (written by **MOLTRA).
keyword(NEL\_F1)
Number of electrons in the *gerade* irreps of the Abelian symmetry group. See below for more information.
keyword(NEL\_F2)
Number of electrons in the *ungerade* irreps of the Abelian symmetry group.
.. index:: .PRINT
.. _RELCC_.PRINT:
.PRINT
------
Print level.
*Default:*
::
.PRINT
0
.. index:: .TIMING
.. _RELCC_.TIMING:
.TIMING
-------
Print detailed timing information.
*Default:*
::
Only limited timing information is printed.
.. index:: .DEBUG
.. _RELCC_.DEBUG:
.DEBUG
------
Print debug information.
*Default:*
::
Debug information is not printed.
.. index:: *CCENER
.. _*CCENER:
========
\*CCENER
========
Covers options related to energy.
.. index:: .NOMP2
.. _CCENER_.NOMP2:
.NOMP2
------
Deactivate MP2 calculation.
.. index:: .NOSD
.. _CCENER_.NOSD:
.NOSD
-----
Deactivate CCSD calculation.
.. index:: .NOSDT
.. _CCENER_.NOSDT:
.NOSDT
------
Deactivate the calculation of perturbative triples. This is potentially
useful when running into memory problems for very big calculations and will also
save some CPU time.
.. index:: .MAXIT
.. _CCENER_.MAXIT:
.MAXIT
------
Set maximum number of iterations allowed to solve the CC equations.
.. index:: .MAXDIM
.. _CCENER_.MAXDIM:
.MAXDIM
-------
Set maximum number of amplitude vectors used in the DIIS extrapolation.
.. index:: .NTOL
.. _CCENER_.NTOL:
.NTOL
-----
Specify requested convergence (10^-NTOL) in the amplitudes.
.. index:: .NOSING
.. _CCENER_.NOSING:
.NOSING
-------
Eliminate T1 amplitudes in the calculation (only interesting for test
purposes, this gives no computational speed-up).
.. index:: .NODOUB
.. _CCENER_.NODOUB:
.NODOUB
-------
Eliminate T2 amplitudes in the calculation (only interesting for test
purposes, this gives no computational speed-up).
Deactivate contribution from doubles; corresponds to a CCS calculation.
.. index:: *CCFOPR
.. _*CCFOPR:
========
\*CCFOPR
========
Calculate first-order properties (expectation values) for the MP2 wave
function.
.. index:: .NATORB
.. _CCFOPR_.NATORB:
.NATORB
-------
Calculate natural orbitals (currently only for MP2 density matrix)
.. index:: .UNRELAXED
.. _CCFOPR_.UNRELAXED:
.UNRELAXED
----------
Use unrelaxed density matrix (computationally cheaper but less accurate)
.. index:: *CCFSPC
.. _*CCFSPC:
========
\*CCFSPC
========
Perform a Fock space MRCC calculation in which a model space is
correlated and then diagonalized to give CC energies for a set of
states.
.. index:: .DOIH
.. _CCFSPC_.DOIH:
.DOIH
-----
Use the Intermediate Hamiltonian formalism in which an auxiliary space
is used to prevent the "intruder state" problem. Default: IH formalism
not used.
.. index:: .DOEA
.. _CCFSPC_.DOEA:
.DOEA
-----
Calculate electron affinities (add one electron to the reference state,
allowing occupation of the active virtual orbitals)
.. index:: .DOIE
.. _CCFSPC_.DOIE:
.DOIE
-----
Calculate ionization energies (remove one electron from the reference
state, allowing depletion of the active occupied orbitals)
.. index:: .DOEA2
.. _CCFSPC_.DOEA2:
.DOEA2
------
Calculate second electron affinities (add two electrons to the reference
state, allowing occupation of the active virtual orbitals)
.. index:: .DOIE2
.. _CCFSPC_.DOIE2:
.DOIE2
------
Calculate second ionization energies (remove two electrons from the
reference state, allowing depletion of the active occupied orbitals)
.. index:: .DOEXC
.. _CCFSPC_.DOEXC:
.DOEXC
------
Calculate excitation energies (allow excitation from the set of active
occupied orbitals to the set of active virtual orbitals)
.. index:: .NACTH
.. _CCFSPC_.NACTH:
.NACTH
------
Specification of the set of active hole orbitals (from which
ionization/excitation takes place)
.. index:: .NACTP
.. _CCFSPC_.NACTP:
.NACTP
------
Specification of the set of active particle orbitals (to which electron
attachment/excitation takes place)
.. index:: .MAXIT
.. _CCFSPC_.MAXIT:
.MAXIT
------
Maximum number of iterations allowed to solve the FSCC equations
.. index:: .MAXDIM
.. _CCFSPC_.MAXDIM:
.MAXDIM
-------
Set maximum number of amplitude vectors used in the DIIS extrapolation.
.. index:: .NTOL
.. _CCFSPC_.NTOL:
.NTOL
-----
Specify requested convergence (10^-NTOL) in the amplitudes.
.. index:: .GESTAT
.. _CCFSPC_.GESTAT:
.GESTAT
-------
Specify the state number in the last active sector to pick the energy
from (remember to account for degeneracies) for a state-specific FSCC
geometry optimization based on a numerical gradient.
.. index:: *CCIH
.. _*CCIH:
======
\*CCIH
======
Options for intermediate hamiltonian in FSCC.
.. index:: .EHMIN
.. _CCIH_.EHMIN:
.EHMIN
------
Minimum orbital energy of occupied orbitals forming the auxiliary (Pi)
space. Orbitals with energies lower than this energy are taken in the
secundary (Q) space and do not contribute to the model space.
low limit of orbital energies of active occupied
orbitals, which constitute the secondary Pi space. Could be
used in (1,0), (2,0) and (1,1) sectors. Arguments: real.
.. index:: .EHMAX
.. _CCIH_.EHMAX:
.EHMAX
------
Maximum orbital energy of occupied orbitals forming the auxiliary (Pi)
space. Orbitals with energies higher than this energy are taken in the
primary (Pm) space.
This is upper limit of one-electronic energies of
active occupied orbitals, which constitute the secondary Pi space.
Could be used in (1,0), (2,0) and (1,1) sectors. Arguments: real.
.. index:: .EPMIN
.. _CCIH_.EPMIN:
.EPMIN
------
Minimum orbital energy of virtual orbitals forming the auxiliary (Pi)
space. Orbitals with energies lower than this energy are taken in the
primary (Pm) space.
This is the low limit of orbital energies of active virtual
orbitals, which constitute the secondary Pi space.
Could be used in (0,1), (0,2) and (1,1) sectors. Arguments: real.
.. index:: .EPMAX
.. _CCIH_.EPMAX:
.EPMAX
------
Maximum orbital energy of virtual orbitals forming the auxiliary (Pi)
space. Orbitals with energies higher than this energy are taken in the
secundary (Q) space and do not contribute to the model space.
This is the upper limit of one-electronic energies of active virtual
orbitals, which constitute the secondary Pi space. Could be
used in (0,1), (0,2) and (1,1) sectors. Arguments: real.
====================================================
Other Intermediate Hamiltonian (IH) input parameters
====================================================
For *experts* only.
Following keyowrds belong to the CCIH namelist section.
.. index:: .IHSCHEME
.. _CCIH_.IHSCHEME:
.IHSCHEME
---------
Choose particular IH scheme. Arguments: Integer IHSCHEME = 1, or 2.
The IHSCHEME =1 corresponds to the extrapolated IH
(XIH) approach, described in the paper :cite:`Eliav2005`.
Main idea: proper modification of the energetic denominators,
containing problematic Pi -> Q transition.
The original denominator 1/(E_Pi - E_Q) , used during CC
iterations, is substituted by the following expression (1)
.. math::
\frac{(1-[\frac{AIH*SHIFT}{(E_{Pi} - E_{Q} + SHIFT)}]^{NIH})}{\frac{(1-AIH*SHIFT}{(E_{Pi} - E_{Q} + SHIFT))}},
here AIH, SHIFT,NIH are parameters, specially chosen for
overcoming of the intruder states problem. These parameters
could be used in the procedure of the extrapolation of
the "exact" effective Hamiltonian solutions from
corresponding IH CC energies and wave functions.
IHSCHEME =2 corresponds to the simplified IH-2 approach, described in
the paper :cite:`Landau2004`.
Here the problematic denominators :math:`1/(E_{Pi} - E_{Q})` are
substituted simply by the factor 0.
Default: IHSCHEME = 2
Next key options are used only in case of XIH (IHSCHEME = 1).
.. index:: .SHIFTH11
.. _CCIH_.SHIFTH11:
.SHIFTH11
---------
Energy shift for the one-electronic excitations in (1,0) sector. Arguments: real.
.. index:: .SHIFTH12
.. _CCIH_.SHIFTH12:
.SHIFTH12
---------
Energy shift for the two-electronic excitations in (1,0) sector. Arguments: real.
.. index:: .SHIFTH2
.. _CCIH_.SHIFTH2:
.SHIFTH2
--------
Energy shift for the two-electronic excitations in (2,0) sector. Arguments: real.
.. index:: .SHIFTP11
.. _CCIH_.SHIFTP11:
.SHIFTP11
---------
Energy shift for the one-electronic excitations in (0,1) sector. Arguments: real.
.. index:: .SHIFTP12
.. _CCIH_.SHIFTP12:
.SHIFTP12
---------
Energy shift for the two-electronic excitations in (0,1) sector. Arguments: real.
.. index:: .SHIFTP2
.. _CCIH_.SHIFTP2:
.SHIFTP2
--------
Energy shift for the two-electronic excitations in (0,2) sector. Arguments: real.
Usually we choose the approximate difference between the highest orbital energy belonging to Pi and the lowest
orbital energy belonging to the Pm space.
Works only with the old style of RELCC input.
.. index:: .SHIFTHP
.. _CCIH_.SHIFTHP:
.SHIFTHP
--------
Energy shift for the two-electronic excitations in (1,1) sector. Arguments: real
.. index:: .AIH
.. _CCIH_.AIH:
.AIH
----
Compensation factor, used in expression (1). Arguments: real positive, not greater then 1.0.
.. index:: .NIH
.. _CCIH_.NIH:
.NIH
----
Compensation power, used in expression (1). Arguments: integer.
In the case of the limit: AIH=1.0 and NIH -> "infinity" ( NIH>100, in practice)
we have so called "full compensation" method, corresponding to the
extrapolation of the effective Hamiltonian from the intermediate one.
.. index:: *CCSORT
.. _*CCSORT:
========
\*CCSORT
========
Specialist options related to the sorting of two-electron integrals
and the calculation of the reference Fock matrix.
.. index:: .NORECMP
.. _CCSORT_.NORECMP:
.NORECMP
--------
Do not recompute the Fock matrix, but assume a diagonal matrix with the orbital energies
taken from the SCF program on the dioagonal. This is usually not recommended as the latter
correspond to a restricted open shell expression and RELCCSD uses an unrestricted formalism.
For closed shell systems the two expressions are identical and this option merely suppresses
a build-in check on the accuracy of transformed integrals.
.. index:: .USEOE
.. _CCSORT_.USEOE:
.USEOE
------
Ignore recomputed Fock matrix and use orbital energies supplied by the SCF program. This option is
sometimes useful for degenerate open shell cases in which case the perturbation theory for the unrestricted formalism
is not invariant for rotations among degenerate orbitals. It should only change the outcome
of the [T], (T) and -T energy corrections.