Time-dependent polarization functions

We start with the time dependent polarization function P(t). This quantity may be expressed in the form of an integral equation ... 

Chandra and his coworkers have developed analytical theories to predict and explain the interfacial solvation dynamics. For example, Chandra et al. [61] have developed a time-dependent density functional theory to predict polarization relaxation at the solid-liquid interface. They find that the interfacial molecules relax more slowly than does the bulk and that the rate of relaxation changes nonmonotonically with distance from the interface They attribute the changing relaxation rate to the presence of distinct solvent layers at the interface. Senapati and Chandra have applied theories of solvents at interfaces to a range of model systems [62-64]. 

One may write the time-dependent wave function in the polar form, viz.,... 

Debye obtained his result by solving a forced diffusion equation Ci.e., with torque of the applied field included) for the distribution of dipole coordinate p - pcosS, with 6 the polar angle between the dipole axis and tSe field, and the same result for the model follows very simply from equation (3) using the time dependent distribution function in the absence of the field (5). The relaxation time is given by td = 1/2D, which for a molecular sphere of volume v rotating in fluid of viscosity n becomes... 

The derivation of formulae for the frequency-dependent nonlinear susceptibilities of nonlinear optics from the time-dependent response functions can be found in a number of sources, (Bloembergen, Ward and New, Butcher and Cotter, Flytzanis ). Here it is assumed that the susceptibilities can be expressed in terms of frequency-dependent quantities that connect individual (complex) Fourier components of the polarization with simple products of the Fourier components of the field. What then has to be shown is how the quantities measured in various experiments can be reduced to these simpler parameters. 

Photoinduced birefringence illustrates what is happening in bulk, but use of time-dependent polarization modulation infrared spectroscopy can offer a detailed insight at the molecular level. Different infrared bands can be monitored in time, and the change in their orientation function allows one almost to watch the different groups move in real time. In order to do this, the process must be slowed down considerably, which is achieved by using a fraction of the pump intensity. 

The drastical enhancement of the TPA cross section in the presence of the solvent for the two-photon polymerization initiator [4-trans-[p-(N,N-Di- -butylamino)-p-stilbenyl vinyl piridinc (DBASVP) has been illustrated in a recent work by Wong et al. [112]. The DBASVP is the typical D-tt-A molecule exhibiting the positive solvatochromism (scheme 7). Hence, the lowest excited state of the DBASVP molecule has been found to be a CT state, which completely dominates the linear absorption spectrum. Wong et al. have combined the time-dependent density functional theory and the polarized continuum model (PCM) to evaluate the solvatochromic shift, TPA cross-section, and oxidation potential of the DBASVP molecule in different solutions. 

To illustrate the exchange of the phase information between the atomic transition and the multipole field, consider the electric dipole Jaynes-Cummings model (34). Assume that the field consists of two circularly polarized components in a coherent state each. The atom is supposed to be initially in the ground state. Then, the time-dependent wave function of the system has the form [53]... 

The simplest polarization propagator corresponds to choosing an HF reference and including only the h2 operator, known as the Random Phase Approximation (RPA), which is identical to Time-Dependent Hartree-Fock (TDHF), with the corresponding density functional version called Time-Dependent Density Functional Theory (TDDFT). For the static case co= 0) the resulting equations are identical to those obtained from a coupled Hartree-Fock approach (Section 10.5). When used in conjunction with coupled cluster wave functions, the approach is usually called Equation Of Motion (EOM) methods. ... 

In a later paper, Casida et used this formalism to calculate the excitation energies of some smaller molecules (N2, CO, CH2, and C2H4). In Table 12 we have collected their results for N2 and in Table 13 those for CO. Those for N2 can be compared directly with those of Table 11 obtained with an exact-exchange method. The results of both tables show that the time-dependent density-functional methods give results that are almost as accurate as those of the sophisticated correlation methods (like coupled-cluster, configuration-interaction, multiple-configuration, or polarization-propagator methods) and considerably... 

A special problem where the time-dependent density-functional theory could be useful is that of calculating the polarizability and hyperpolarizability (Section 12). It turned out that although accurate results could be achieved for smaller molecules (partly, however, requiring a careful choice of the approximate density functional), severe problems could turn up (but did not always) when considering extended systems. It might mean that the current density functionals are lacking an explicit dependence on the polarization, but further studies are needed in order to clarify this point. 

On the contrary, when the time-dependent electric field varies on a time scale faster than the relaxation time of one or more molecular degrees of freedom there is not time to reach at any moment a time-dependent polarization which is in equilibrium with the electric field. In this regime, which is called non-equilibrium polarization, the actual value of polarization will also depend values of the electric field at previous time, and the relation between the polarization of a dielectric medium and the time-dependent polarizing field is phenomenologically described in terms of the whole specuiim of the dielectric permittivity as a function of the frequency co of the oscillating electric field. 

Time dependent polarization charges, 24, 29 Time dependent quasi fl ee-energy functional, 25... 

We will later on in Chapter 7 see that the measurable molecular properties related to the response functions or polarization propagators are normally defined to be frequency dependent and not time dependent. We define therefore the Fourier transform of the time dependent polarization propagator as... 

However, Eq. (3.108) is only the definition of the Fourier transform, which then has to be applied to the expression for the time-dependent polarization propagator in Eq. (3.107) with t —t replaced by t. This leads us [see Exercise 3.8] to the spectral representation of the polarization propagator or hnear response function... 

Excitation energies and transition moments can in principle be obtained as poles and residua of polarization propagators as discussed in Section 7.4. However, only in the case that the set of operators hn in Eq. (7.77) is restricted to single excitation and de-excitation operators q i,qai is it computationally feasible to determine all excitation energies. This restricts this approach to single-excitation-based methods like the random phase approximation (RPA) discussed in Sections 10.3 and 11.1 or time-dependent density functional theory (TD-DFT). 

This approximation is better known as the time-dependent Hartree—Fock approximation (TDHF) (McLachlan and Ball, 1964) (see Section 11.1) or random phase approximation (RPA) (Rowe, 1968) and can also be derived as the linear response of an SCF wavefunction, as described in Section 11.2. Furthermore, the structure of the equations is the same as in time-dependent density functional theory (TD-DFT), although they differ in the expressions for the elements of the Hessian matrix E22. The polarization propagator in the RPA is then given as... 

The formulation of approximate response theory based on an exponential parame-trization of the time-dependent wave function, Eq. (11.36), and the Ehrenfest theorem, Eq. (11.40), can also be used to derive SOPPA and higher-order Mpller-Plesset perturbation theory polarization propagator approximations (Olsen et al., 2005). Contrary to the approach employed in Chapter 10, which is based on the superoperator formalism from Section 3.12 and that could not yet be extended to higher than linear response functions, the Ehrenfest-theorem-based approach can be used to derive expressions also for quadratic and higher-order response functions. In the following, it will briefly be shown how the SOPPA linear response equations, Eq. (10.29), can be derived with this approach. 

Thereby Snse(Q, 0 is the effective time-dependent polarization detected at the instmment. We note that the conversion of the ideal intensities to S(Q, f) under real experimental conditions need to consider polarization losses in the instmment that are measured in a reference experiment on a purely elastic scatterer. Such a measurement serves as a determination of the experimental resolution function. Since the experiments are performed in Fourier space instead of a deconvolution as required in resolution corrections merely requires a division of the sample data by the reference measurement. 

Fig. 4.6. Piezoelectric pulse diagrams can be used to obtain explicit representations of the time dependent electric fields in piezoelectric substances. The magnitudes and orientations of these electric fields are critical to development of shock-induced conduction. As an example, the diagram on the left shows the polarization and displacement relations for a location at the input electrode. The same functions for a location within the crystal is shown on the right (after Davison and Graham [79D01]).

The linear piezoeleetrie model can be used to demonstrate that the magnitude of the electric field encountered for a given polarization function is a sensitive function of the thickness of the sample. This behavior can be demonstrated by noting that the electric displacement at a given time is inversely proportional to the thickness. Thus, the thickness of the sample is an important variable for investigating effects such as conductivity that depend upon the magnitude of the electric field. Conversely, various input stress wave shapes can be used to cause various field distributions at fixed thicknesses. 

All these features were observed experimentally for solutions of 3-amino-/V-methylphthalimide, 4-amino-/V-methylphthalimide, and for nonsubstituted rhoda-mine. The results were observed for cooled, polar solutions of phthalimides, in which the orientational relaxation is delayed. Exactly the same spectral behavior was observed [50] by picosecond spectroscopy for low viscosity liquid solutions at room temperature, in which the orientational relaxation rate is much higher. All experimental data indicate that correlation functions of spectral shifts Av-l(t), which are used frequently for describing the Time Dependent Stokes Shift, are essentially the functions of excitation frequency. 

Another pertinent observation is the fact that the reaction proceeded twice as fast in -butyraldehyde (polar) as in benzene (nonpolar), even though the catalyst concentration was reduced to only one-third the comparable level. A graphic illustration of this effect is given in Fig. 9. The rate of gas uptake is plotted as a function of time for a reaction conducted in benzene and again for a second reaction conducted in butyraldehyde. The rate of reaction in the polar solvent was initially fast and decreased with time. The rate in the nonpolar benzene was initially slow, became faster as the solvent became more polar with the presence of product aldehyde, and then subsequently diminished with time. When the data were replotted as the log of unreacted olefin vs. time, the polar medium reaction showed first-order dependence on olefin concentration, whereas the nonpolar solvent reaction showed no definite order, owing to the constantly changing polarity. 

---