Perturbation theory response

The second-order nonlinear optical processes of SHG and SFG are described correspondingly by second-order perturbation theory. In this case, two photons at the drivmg frequency or frequencies are destroyed and a photon at the SH or SF is created. This is accomplished tlnough a succession of tlnee real or virtual transitions, as shown in figure Bl.5.4. These transitions start from an occupied initial energy eigenstate g), pass tlnough intennediate states n ) and n) and return to the initial state g). A fiill calculation of the second-order response for the case of SFG yields [37]... 

The perturbation theory presented in Chapter 2 implies that orientational relaxation is slower than rotational relaxation and considers the angular displacement during a free rotation to be a small parameter. Considering J(t) as a random time-dependent perturbation, it describes the orientational relaxation as a molecular response to it. Frequent and small chaotic turns constitute the rotational diffusion which is shown to be an equivalent representation of the process. The turns may proceed via free paths or via sudden jumps from one orientation to another. The phenomenological picture of rotational diffusion is compatible with both... 

This conclusion is not unexpected. The molecule s response to weak collisions considered as a random process may be described by perturbation theory (with respect to interaction) if [53]... 

In this section we consider how to express the response of a system to noise employing a method of cumulant expansions [38], The averaging of the dynamical equation (2.19) performed by this technique is a rigorous continuation of the iteration procedure (2.20)-(2.22). It enables one to get the higher order corrections to what was found with the simplest perturbation theory. Following Zatsepin [108], let us expound the above technique for a density of the conditional probability which is the average... 

The concerns we have expressed are bound to get even more acute if the problem under study demands that we are able to adequately describe distortion effects induced in the electron distribution by external fields. The evaluation of linear (and, still more, non linear) response funetions [1] by perturbation theory then forces one to take care also of the nonoccupied portion of the complete orbital spectrum, which is entrusted with the role of representing the polarization caused by the external fields in the unperturbed electron distribution [4], ... 

The various response tensors are identified as terms in these series and are calculated using numerical derivatives of the energy. This method is easily implemented at any level of theory. Analytic derivative methods have been implemented using self-consistent-field (SCF) methods for a, ft and y, using multiconfiguration SCF (MCSCF) methods for ft and using second-order perturbation theory (MP2) for y". The response properties can also be determined in terms of sum-over-states formulation, which is derived from a perturbation theory treatment of the field operator — [iE, which in the static limit is equivalent to the results obtained by SCF finite field or analytic derivative methods. 

Although the (3-nc orbital is formally vacant in the cation, Table 3.7 shows that a small residual population (0.0346c) survives in this orbital. This occupancy can be attributed to a strong donor-acceptor interaction with the filled no orbital as depicted in Fig. 3.14. This n0 nc interaction is estimated by second-order perturbation theory (Eq. (1.24)) to stabilize the ion by 19.5 kcal mol-1, a significant delocalization that is primarily responsible for the slightly lower %p(L) value in this ion. 

In this chapter we will focus on one particular, recently developed DFT-based approach, namely on first-principles (Car-Parri-nello) molecular dynamics (CP-MD) [9] and its latest advancements into a mixed quantum mechanical/molecular mechanical (QM/MM) scheme [10-12] in combination with the calculation of various response properties [13-18] within DFT perturbation theory (DFTPT) and time-dependent DFT theory (TDDFT) [19]. 

Density-functional Perturbation Theory and the Calculation of Response Properties 21... 

There are several possible ways of deriving the equations for TDDFT. The most natural way departs from density-functional perturbation theory as outlined above. Initially it is assumed that an external perturbation is applied, which oscillates at a frequency co. The linear response of the system is then computed, which will be oscillating with the same imposed frequency co. In contrast with the standard static formulation of DFPT, there will be special frequencies cov for which the solutions of the perturbation theory equations will persist even when the external field vanishes. These particular solutions for orbitals and frequencies describe excited electronic states and energies with very good accuracy. 

The actual form of the Hamiltonian operator hp does not have to be defined at this moment. As in standard perturbation theory, it is assumed that the solution of the electronic structure problem of the combined Hamiltonian HKS +HP can be described as the solution y/(0) of HKS, corrected by a small additional linear-response wavefunction /b//(,). Only these response orbitals will explicitly depend on time - they will follow the oscillations of the external perturbation and adopt its time dependency. Thus, the following Ansatz is made for the solution of the perturbed Hamiltonian HKS +HP ... 

A Rayleigh-Schrodinger-type perturbation theory has recently been developed by Angyan [107], The consideration of external perturbations, like electric fields, permits the calculation of response functions for solvated species. 

We call the Fukui function / (r) the HOMO response. Equation 24.39 is demonstrated as follows. The PhomoW is the so-called Kohn-Sham Fukui function denoted as f (r) [32]. According to the first-order perturbation theory, one has... 

It is worth noting that screened response x/C r ) can be computed from the Kohn-Sham orbital wave functions and energies using standard first-order perturbation theory [3]... 

The exchange contribution in an ab initio perturbation theory is the only repulsive term ) around the energy minimum in most of the stable complexes and consequently we would expect no net repulsion between two closed shell molecules in semiempirical calculations. On the other hand it is known from actual calculations that intermolecular interactions are described more or less correctly by the CNDO/2 procedure. Indeed, strong repulsion is obtained between closed shell molecules. Evidently there must be another approximation which compensates for the neglect of exchange energy. With regard to the simplifications of the CNDO/2 method we find that this is in fact the case. The approximation shown in Eq. (17) is responsible for the repulsive term. 

However, until today no systematic comparison of methods based on MpUer-Plesset perturbation (MP) and Coupled Cluster theory, the SOPPA or multiconfigurational linear response theory has been presented. The present study is a first attempt to remedy this situation. Calculations of the rotational g factor of HF, H2O, NH3 and CH4 were carried out at the level of Hartree-Fock (SCF) and multiconfigurational Hartree-Fock (MCSCF) linear response theory, the SOPPA and SOPPA(CCSD) [40], MpUer-Plesset perturbation theory to second (MP2), third (MP3) and fourth order without the triples contributions (MP4SDQ) and finally coupled cluster singles and doubles theory. The same basis sets and geometries were employed in all calculations for a given molecule. The results obtained with the different methods are therefore for the first time direct comparable and consistent conclusions about the performance of the different methods can be made. 

---