:orphan:
 

.. index:: *QUADRATIC RESPONSE
.. _*QUADRATIC RESPONSE:

====================
\*QUADRATIC RESPONSE
====================

This section gives directives for the calculation of quadratic response
functions :cite:`Saue2002a`.


General control statements
==========================

.. index:: .PRINT
.. _QUADRATIC_RESPONSE_.PRINT:

.PRINT
------

Print level.

*Default:*

::

    .PRINT
     0

Definition of the quadratic response function
=============================================

.. index:: .DIPLEN
.. _QUADRATIC_RESPONSE_.DIPLEN:

.DIPLEN
-------

Specification of dipole operators for A, B, and C
(see :ref:`one_electron_operators` for details).

.. index:: .A OPERATOR
.. _QUADRATIC_RESPONSE_.A OPERATOR:

.A OPERATOR
-----------

Specification of the A operator
(see :ref:`one_electron_operators` for details).

.. index:: .B OPERATOR
.. _QUADRATIC_RESPONSE_.B OPERATOR:

.B OPERATOR
-----------

Specification of the B operator
(see :ref:`one_electron_operators` for details).

.. index:: .C OPERATOR
.. _QUADRATIC_RESPONSE_.C OPERATOR:

.C OPERATOR
-----------

Specification of the C operator
(see :ref:`one_electron_operators` for details).


.. index:: .B FREQ
.. _QUADRATIC_RESPONSE_.B FREQ:

.B FREQ
-------

Specify frequencies of operator B.

*Example:* 3 different frequencies.

::

    .B FREQ
     3
     0.001
     0.002
     0.01

*Default:* Static case.

::

    .B FREQ
     1
     0.0

.. index:: .C FREQ
.. _QUADRATIC_RESPONSE_.C FREQ:

.C FREQ
-------

Specify frequencies of operator C (see :ref:`QUADRATIC_RESPONSE_.B FREQ`).

.. index:: .ALLCMB
.. _QUADRATIC_RESPONSE_.ALLCMB:

.ALLCMB
-------

Evaluate all nonzero quadratic response functions and thereby
disregarding analysis of overall permutational symmetry.

*Default:* Evaluate only unique, nonzero, response functions.

Excited state properties
========================

**This page describes unreleased functionality. The keywords may not be
available in your version of DIRAC.**

First order properties of excited states can be computed from the
quadratic response function.

.. index:: .EXCPRP
.. _QUADRATIC_RESPONSE_.EXCPRP:

.EXCPRP
-------

Give the number of "left" and "right" states in each boson symmetry.

*Example*:

::

    .EXCPRP
    5 5 5 5
    0 0 0 0

Compute the excited state expectation values \|\\langle i\|\\hat{A}\|i\\rangle|,
where i goes from 1 to 5 in each symmetry (four symmetries in this case). The
zeros can be substituted with positive integers to generate transition state moments \|\\langle i\|\\hat{A}\|j\\rangle|.

Control variational parameters
==============================

.. index:: .SKIPEE
.. _QUADRATIC_RESPONSE_.SKIPEE:

.SKIPEE
-------

Exclude all rotations between occupied positive-energy and virtual
positive-energy orbitals.

.. index:: .SKIPEP
.. _QUADRATIC_RESPONSE_.SKIPEP:

.SKIPEP
-------

Exclude all rotations between occupied positive-energy and virtual
negative-energy orbitals.

Control reduced equations
=========================

.. index:: .MAXITR
.. _QUADRATIC_RESPONSE_.MAXITR:

.MAXITR
-------

Maximum number of iterations.

*Default:*

::

    .MAXITR
     30

.. index:: .MAXRED
.. _QUADRATIC_RESPONSE_.MAXRED:

.MAXRED
-------

Maximum dimension of matrix in reduced system.

*Default:*

::

    .MAXRED
     100

.. index:: .THRESH
.. _QUADRATIC_RESPONSE_.THRESH:

.THRESH
-------

Threshold for convergence of reduced system.

*Default:*

::

    .THRESH
     1.0D-5

Control integral contributions
==============================

The user is encouraged to experiment with these options since they may
have an important effect on run time.

.. index:: .INTFLG
.. _QUADRATIC_RESPONSE_.INTFLG:

.INTFLG
-------

Specify what two-electron integrals to include
(default: :ref:`HAMILTONIAN_.INTFLG` under :ref:`**HAMILTONIAN`).

.. index:: .CNVINT
.. _QUADRATIC_RESPONSE_.CNVINT:

.CNVINT
-------

Set threshold for convergence before adding SL and SS integrals to
SCF-iterations.

*2 (real) Arguments:*

::

    .CNVINT
     CNVXQR(1) CNVXQR(2)

*Default:* Very large numbers.

.. index:: .ITRINT
.. _QUADRATIC_RESPONSE_.ITRINT:

.ITRINT
-------

Set the number of iterations before adding SL and SS integrals to
SCF-iterations.

*Default:*

::

    .ITRINT
     1 1

Control trial vectors
=====================

.. index:: .XQRNRM
.. _QUADRATIC_RESPONSE_.XQRNRM:

.XQRNRM
-------

Normalize trial vectors. Using normalized trial vectors will reduce
efficiency of screening.

*Default:* Use un-normalized vectors.

Advanced/debug flags
====================

.. index:: .NOPREC
.. _QUADRATIC_RESPONSE_.NOPREC:

.NOPREC
-------

No preconditioning of initial trial vectors.

*Default:* Preconditioning of trial vectors.

.. index:: .RESFAC
.. _QUADRATIC_RESPONSE_.RESFAC:

.RESFAC
-------

New trial vector will be generated only for variational parameter
classes whose residual has a norm that is larger than a fraction
1/RESFAC of the maximum norm.

*Default:*

::

    .RESFAC
     1000.0