Polarizabilities second-order

Electrostatic terms other than the simple charge interactions above are commonly included in molecular mechanics calculations. particularly dipole-dipole interactions. More recently, second-order electrostatic interactions like those describing polarizability have been added to some force fields. 

Examination of the energy expression 17.7 shows that the polarizability components are the second-order energies and a little analysis shows that (for example)... 

The second-order property, the dipole polarizability, as given by the derivative formula eq. (10.31), is... 

The poles con espond to excitation energies, and the residues (numerator at the poles) to transition moments between the reference and excited states. In the limit where cj —> 0 (i.e. where the perturbation is time independent), the propagator is identical to the second-order perturbation formula for a constant electric field (eq. (10.57)), i.e. the ((r r))Q propagator determines the static polarizability. 

These concepts play an important role in the Hard and Soft Acid and Base (HSAB) principle, which states that hard acids prefer to react with hard bases, and vice versa. By means of Koopmann s theorem (Section 3.4) the hardness is related to the HOMO-LUMO energy difference, i.e. a small gap indicates a soft molecule. From second-order perturbation theory it also follows that a small gap between occupied and unoccupied orbitals will give a large contribution to the polarizability (Section 10.6), i.e. softness is a measure of how easily the electron density can be distorted by external fields, for example those generated by another molecule. In terms of the perturbation equation (15.1), a hard-hard interaction is primarily charge controlled, while a soft-soft interaction is orbital controlled. Both FMO and HSAB theories may be considered as being limiting cases of chemical reactivity described by the Fukui ftinction. 

The second procedure, several aspects of which are reviewed in this paper, consists of directly computing the asymptotic value by employing newly-developed polymeric techniques which take advantage of the one-dimensional periodicity of these systems. Since the polarizability is either the linear response of the dipole moment to the field or the negative of the second-order term in the perturbation expansion of the energy as a power series in the field, several schemes can be proposed for its evaluation. Section 3 points out that several of these schemes are inconsistent with band theory summarized in Section 2. In Section 4, we present the main points of the polymeric polarization propagator approaches we have developed, and in Section 5, we describe some of their characteristics in applications to prototype systems. 

Remembering that the polarizabilities are expressed in units of it follows that the largest positive eigenvalue corresponds to the energetically most favorable second-order bond distortion, and the type... 

The second part of this paper concerns the choice of the atomic basis set and especially the polarization functions for the calculation of the polarizability, o , and the hyperpo-larizabiliy, 7. We propose field-induced polarization functions (6) constructed from the first- and second-order perturbed hydrogenic wavefunctions respectively for a and 7. In these polarization functions the exponent ( is determined by optimization with the maximum polarizability criterion. These functions have been successfully applied to the calculation of the polarizabilities, a and 7, for the He, Be and Ne atoms and the molecule. 

This calculation has shown the importance of the basis set and in particular the polarization functions necessary in such computations. We have studied this problem through the calculation of the static polarizability and even hyperpolarizability. The very good results of the hyperpolarizabilities obtained for various systems give proof of the ability of our approach based on suitable polarization functions derived from an hydrogenic model. Field—induced polarization functions have been constructed from the first- and second-order perturbed hydrogenic wavefunctions in which the exponent is determined by optimization with the maximum polarizability criterion. We have demonstrated the necessity of describing the wavefunction the best we can, so that the polarization functions participate solely in the calculation of polarizabilities or hyperpolarizabilities. 

A computer program for the theoretical determination of electric polarizabilities and hyperpolarizabilitieshas been implemented at the ab initio level using a computational scheme based on CHF perturbation theory [7-11]. Zero-order SCF, and first-and second-order CHF equations are solved to obtain the corresponding perturbed wavefunctions and density matrices, exploiting the entire molecular symmetry to reduce the number of matrix element which are to be stored in, and processed by, computer. Then a /j, and iap-iS tensors are evaluated. This method has been applied to evaluate the second hyperpolarizability of benzene using extended basis sets of Gaussian functions, see Sec. VI. 

Cyvin, S. J., Rauch, J. E. and Decius, J. C. (1965) Theory of hyper-Raman effects (nonlinear inelastic light scattering) selection rules and depolarization ratios for the second-order polarizability. 

The determination of the number of the SHG active complex cations from the corresponding SHG intensity and thus the surface charge density, a°, is not possible because the values of the molecular second-order nonlinear electrical polarizability, a , and molecular orientation, T), of the SHG active complex cation and its distribution at the membrane surface are not known [see Eq. (3)]. Although the formation of an SHG active monolayer seems not to be the only possible explanation, we used the following method to estimate the surface charge density from the SHG results since the square root of the SHG intensity, is proportional to the number of SHG active cation com-... 

As electric fields and potential of molecules can be generated upon distributed p, the second order energies schemes of the SIBFA approach can be directly fueled by the density fitted coefficients. To conclude, an important asset of the GEM approach is the possibility of generating a general framework to perform Periodic Boundary Conditions (PBC) simulations. Indeed, such process can be used for second generation APMM such as SIBFA since PBC methodology has been shown to be a key issue in polarizable molecular dynamics with the efficient PBC implementation [60] of the multipole based AMOEBA force field [61]. 

A novel second-order nonlinear optical medium which should offer considerable fabrication flexibility has been described. The physics of alignment of the highly nonlinearly polarizable moiety was discussed. However, observation of complex dynamical and thermal behavior indicates that an important role is played by the polymer liquid crystalline host. Additional properties of modified members of this family of lc polymers were consequently investigated. The explanations of guest alignment stabilization and thermal dependence of the alignability remain unresolved issues. 

The proportionality constants a and (> are the linear polarizability and the second-order polarizability (or first hyperpolarizability), and x(1) and x<2) are the first- and second-order susceptibility. The quadratic terms (> and x<2) are related by x(2) = (V/(P) and are responsible for second-order nonlinear optical (NLO) effects such as frequency doubling (or second-harmonic generation), frequency mixing, and the electro-optic effect (or Pockels effect). These effects are schematically illustrated in Figure 9.3. In the remainder of this chapter, we will primarily focus on the process of second-harmonic generation (SHG). 

In order to describe second-order nonlinear optical effects, it is not sufficient to treat (> and x<2) as a scalar quantity. Instead the second-order polarizability and susceptibility must be treated as a third-rank tensors 3p and Xp with 27 components and the dipole moment, polarization, and electric field as vectors. As such, the relations between the dipole moment (polarization) vector and the electric field vector can be defined as ... 

---