A general method is presented for the calculation of molecular properties to arbitrary order at theKohn–Sham density functional level of theory. The quasienergy and Lagrangian formalisms arecombined to derive response functions and their residues by straightforward differentiation of thequasienergy derivative Lagrangian using the elements of the density matrix in the atomic orbitalrepresentation as variational parameters. Response functions and response equations are expressedin the atomic orbital basis, allowing recent advances in the field of linear-scaling methodology to beused. Time-dependent and static perturbations are treated on an equal footing, and atomic basissets that depend on the applied frequency-dependent perturbations may be used, e.g.,frequency-dependent London atomic orbitals. The 2n+1 rule may be applied if computationallyfavorable, but alternative formulations using higher-order perturbed density matrices are alsoderived. These may be advantageous in order to minimize the number of response equations thatneeds to be solved, for instance, when one of the perturbations has many components, as is the casefor the first-order geometrical derivative of the hyperpolarizability.

A density matrix-based quasienergy formulation of Kohn-Sham density functional response theory using perturbation- and time-dependent basis sets

CORIANI, Sonia
2008-01-01

Abstract

A general method is presented for the calculation of molecular properties to arbitrary order at theKohn–Sham density functional level of theory. The quasienergy and Lagrangian formalisms arecombined to derive response functions and their residues by straightforward differentiation of thequasienergy derivative Lagrangian using the elements of the density matrix in the atomic orbitalrepresentation as variational parameters. Response functions and response equations are expressedin the atomic orbital basis, allowing recent advances in the field of linear-scaling methodology to beused. Time-dependent and static perturbations are treated on an equal footing, and atomic basissets that depend on the applied frequency-dependent perturbations may be used, e.g.,frequency-dependent London atomic orbitals. The 2n+1 rule may be applied if computationallyfavorable, but alternative formulations using higher-order perturbed density matrices are alsoderived. These may be advantageous in order to minimize the number of response equations thatneeds to be solved, for instance, when one of the perturbations has many components, as is the casefor the first-order geometrical derivative of the hyperpolarizability.
http://dx.doi.org/10.1063/1.2996351
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11368/1862914
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 97
  • ???jsp.display-item.citation.isi??? ND
social impact