Perturbation methods linear response

More recently, Filatov has developed a Unear response theory for the isomer shift and used it in conjunction with ab initio methods [71-73]. In many respects, this theory is more rigorous than the cahbration approach described earlier. Hence, it will be briefly outhned here. Filatov considered a linear response treatment in which the perturbation was taken as the radial expansion of a finite sized nucleus. The resulting equation for the isomer shift is ... 

We present and analyze the most important simplified free energy methods, emphasizing their connection to more-rigorous methods and the underlying theoretical framework. The simplified methods can all be superficially defined by their use of just one or two simulations to compare two systems, as opposed to many simulations along a complete connecting pathway. More importantly, the use of just one or two simulations implies a common approximation of a near-linear response of the system to a perturbation. Another important theme for simplified methods is the use, in many cases, of an implicit description of solvent usually a continuum dielectric model, often supplemented by a simple description of hydrophobic effects [11]. 

The above theory is usually called the generalized linear response theory because the linear optical absorption initiates from the nonstationary states prepared by the pumping process [85-87]. This method is valid when pumping pulse and probing pulse do not overlap. When they overlap, third-order or X 3 (co) should be used. In other words, Eq. (6.4) should be solved perturbatively to the third-order approximation. From Eqs. (6.19)-(6.22) we can see that in the time-resolved spectra described by x"( ), the dynamics information of the system is contained in p(Af), which can be obtained by solving the reduced Liouville equations. Application of Eq. (6.19) to stimulated emission monitoring vibrational relaxation is given in Appendix III. 

The study of behavior of many-electron systems such as atoms, molecules, and solids under the action of time-dependent (TD) external fields, which includes interaction with radiation, has been an important area of research. In the linear response regime, where one considers the external held to cause a small perturbation to the initial ground state of the system, one can obtain many important physical quantities such as polarizabilities, dielectric functions, excitation energies, photoabsorption spectra, van der Waals coefficients, etc. In many situations, for example, in the case of interaction of many-electron systems with strong laser held, however, it is necessary to go beyond linear response for investigation of the properties. Since a full theoretical description based on accurate solution of TD Schrodinger equation is not yet within the reach of computational capabilities, new methods which can efficiently handle the TD many-electron correlations need to be explored, and time-dependent density functional theory (TDDFT) is one such valuable approach. 

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. 

Static charge-density susceptibilities have been computed ab initio by Li et al (38). The frequency-dependent susceptibility x(r, r cd) can be calculated within density functional theory, using methods developed by Ando (39 Zang-will and Soven (40 Gross and Kohn (4I and van Gisbergen, Snijders, and Baerends (42). In ab initio work, x(r, r co) can be determined by use of time-dependent perturbation techniques, pseudo-state methods (43-49), quantum Monte Carlo calculations (50-52), or by explicit construction of the linear response function in coupled cluster theory (53). Then the imaginary-frequency susceptibility can be obtained by analytic continuation from the susceptibility at real frequencies, or by a direct replacement co ico, where possible (for example, in pseudo-state expressions). 

Several methods can be distinguished within the framework of the perturbative approach. Some [29-37] are based on a multipolar expansion of the operator i.e. the interaction potential of the two species, others rely on the linear response theory [38,39]. 

Banerjee and Harbola [69] have worked out a variation perturbation method within the hydrodynamic approach to the time-dependent density functional theory (TDDFT) in order to evaluate the linear and nonlinear responses of alkali metal clusters. They employed the spherical jellium background model to determine the static and degenerate four-wave mixing (DFWM) y and showed that y evolves almost linearly with the number of atoms in the cluster. 

We saw in Section III that the polarization propagator is the linear response function. The linear response of a system to an external time-independent perturbation can also be obtained from the coupled Hartree-Fock (CHF) approximation provided the unperturbed state is the Hartree-Fock state of the system. Thus, RPA and CHF are the same approximation for time-independent perturbing fields, that is for properties such as spin-spin coupling constants and static polarizabilities. That we indeed obtain exactly the same set of equations in the two methods is demonstrated by Jorgensen and Simons (1981, Chapter 5.B). Frequency-dependent response properties in the... 

The minor disturbance or perturbation method relies on equilibrium theory too and was suggested, for example, by Reilley, Hildebrand, and Ashley (1962). As known from linear chromatography and exploited above already frequently, the retention time of the response to a small pulse injected into a column filled with pure eluent can be used to obtain the initial slope of the isotherm. This approach can be expanded to cover the whole isotherm range. For the example of a singlecomponent system the procedure is as follows (Figure 6.25) the column is equilibrated with a concentration and, once the plateau is established, a small pulse is injected at a time tstarta and a pulse of a different concentration is detected at the corresponding retention time (r a. 

Fig. 1 2.4 The difference of errors of the evaluated self-replication rate constant (output error) evaluated from emulated experiments of type (a), with linearization, and type (b), without linearization, respectively. The figure shows that the variations of the input perturbation and of the experimental error have a different effect on the two types of response methods. The biggest difference occurs for large perturbations, because for large perturbations the linear approach is very inaccurate. (From [11].)...

---