Article Resonance in Interacting Induced-Dipole Polarizing Force Fields: Application to Force-Field Derivatives

The Silberstein model of the molecular polarizability of diatomic molecules, generalized by Applequist et al. for polyatomic molecules, is analyzed. The atoms are regarded as isotropically polarizable points located at their nuclei, interacting via the fields of their induced dipoles. The use of additive values for atom polarizabilities gives poor results, in some cases leading to artificial predictions of absorption bands. The molecular polarizability of methane and its derivative are computed. The agreement with experimental mean molecular polarizabilities is within 1–5%. A hypothesis is indispensable for a suitable representation of polarizability derivative.


Introduction
A simple way to understand the origin of various optical effects in many substances is based on an atom-dipole interaction model.The atoms in a molecule are regarded as isotropic particles, which interact by way of the dipole moments induced in them by an external field.

OPEN ACCESS
In a diatomic molecule AB whose atoms have isotropic polarizabilities  A and  B , respectively, the additive model of polarizability allows calculating the isotropic molecular polarizability  mol =  A +  B .The diagonal form of  mol has only one component: In the interacting induced-dipole polarization model, the molecular polarizability of AB can be placed in a simple and explicit form.The diagonal form of anisotropic  mol has two distinct components  || and   , parallel and perpendicular, respectively, to the bond axis.Siberstein's equations for this case follow [1][2][3]: (2) (3) The mean polarizability  and anisotropy  are defined by: As an example of the manner in which the molecular polarizability depends on the atom polarizabilities, the isotropic molecule methane (CH 4 ) is considered [4].Certain features of the dependence are illustrated by the surface in Figure 1.Additivity of atom polarizabilities would require that the surface be a plane, which is approximately true in the vicinity of the origin, where interactions are small.However, the experimental polarizability of CH 4 is 2.62Å 3 , which is reached in only regions of the surface where the influence of interactions is quite marked.Notice a curve of discontinuity along which the polarizability approaches ±.The behaviour is seen in the polarizability surfaces of several molecules that have been similarly explored.Its origin for diatomic molecules can be seen in Equations ( 2) and (3), where the denominators vanish when  A  B approaches r 6 /4 or r 6 , respectively.For this case the  A and  B are inversely related along the curve of discontinuity, which in the  H - C plane follows the relation  H  C = 0.193Å 6 , which is of the form expected for diatomic molecules, though the numerical constant is not predictable from Equations (2) or (3).
The significance of a polarizability of ± is that the molecule is in a state of resonance and absorbs energy from the applied field, which occurs in spite of the fact that any absorption properties of the atoms are not introduced [4].The behaviour of the model can be understood from its close relation to the classical system of N coupled oscillators, which likewise shows resonance under conditions other than the resonance conditions of the isolated oscillators [5].The absorption properties [6][7][8] of some types of systems can be predicted from the point-dipole interaction approach used in the present report, but it seems doubtful that this could be done reliably for molecules with the model [7]; the neglect of electron exchange between atoms is probably serious.In the immediate vicinity of the resonance condition, the polarizability is bound to be in error because the model does not take into account damping effects, which would prevent the polarizability from going to infinity.The resonance conditions, wherever they appear, are regarded simply as indications that the coupling between atoms far exceeds the extent that can be treated by the model.The purpose of the present report is to present some new derivations of fundamental equations relating molecular electric multipole moments, polarizabilities and their derivatives to the electrostatic energy, and to draw some new conclusions regarding the general characteristics of the polarizability tensors and derivatives.

Molecular Dipolar Polarizability Model
The derivative of a molecular dipolar-polarizability tensor ', with respect to the normal coordinates of vibration Q, is obtained as: where (A jk )' is the derivative of A jk with respect to Q, A jk the inverse polarizability matrix A= E,  the induced dipole moment vector, E the applied electric field vector, and B = A -1 the many-body polarizability matrix.The (A jk )' can be calculated for j ≠ k, from the matrix of transformation from a Cartesian coordinate system in normal coordinates, for the considered vibration mode (A jk = T jk ; j ≠ k, where T jk is the gradient tensor of the electric field E. However, the diagonal elements (A ii )' will be calculated, since the model must take into account the dependence of the atomic polarizabilities vs. the length of the bonds between atoms, which leads us to impose an indispensable supplementary hypothesis: the use of the vibrational correction of atomic polarizabilities to average molecular vibrational stretching modes.The modified atomic polarizabilities are calculated as: where R n is the displacement of the n-th bond length from its experimental value in equilibrium, and  j º the atomic polarizability of the j unit at the rest configuration of the molecule.
By optimizing the d j /dR n relations to represent the best the experimental values of the derivative of the molecular polarizability tensor ', to model the molecular polarizability variations during bond vibrations is possible.

Molecular Dipolar Polarizability of Methane
Equation ( 8) is an indispensable hypothesis for a suitable representation of ', without which the registered variations will be of opposite sign to the experimental magnitudes (cf. Figure 2).a Atomic isotropic and mean additive molecular polarizabilities (Å 3 ).
c Mean effective atomic and mean molecular polarizabilities (Å 3 ).
d Reference bond length C-H (Å) to calculate '.
e Polarizability derivatives for vibrational modes H-C and C-H (Å 2 ).
A set of strategies devised for applications to large systems is reported.The following improvements are implemented in the model: 1.To build up the many-body polarizability matrix A, the atomic polarizability tensors given by 2. Damping functions are used in the calculation of the interaction tensor to prevent the polarizability from going to infinity in the models by Applequist and Birge.
a. Böttcher proposed [25]:  3r pq, r pq, where T (2) is the gradient tensor of the electric field E and  represents the Kronecker delta function: b. Thole proposed [26]: where v pq = r pq /s pq if r pq < s pq , otherwise v pq = 1.The term s is defined as s = 1.662( p  q ) 1/6 .The values of  p are parametrized by ab initio coupled Hartree-Fock.
3. The interactions between bonded atoms and between atoms within a distance lying in an interval defined by [r inf , r sup ] is neglected.The starting values for this interval are [0, 10 30 ] a.u. and r inf is incremented if resonance conditions are detected.
4. The following tests indicating a resonance condition are implemented: a. test whether matrix B -1 is singular, where B is the many-body polarizability matrix; b. test whether matrix B -1 is not defined positive; c. test whether matrix B is not defined positive; d. test whether some effective matrix B p is not defined positive.5. Atomic polarizabilities are improved with vibrational corrections to average molecular vibrational stretching modes.
6.The atomic polarizabilities allow selecting between optimized ab initio and experimental values.
7. The atomic classes can be either read from input or assigned by the program.
Applications are carried out for the following systems.( 1)  000), (b) a rotational isomer in which each unit is rotated by a fixed angle  with respect to the previous one, all rotations performed in the same direction (+++) (0º ≤  ≤ 10º), and (c) rotations performed in the alternate directions (+-+) (0º ≤  ≤ 10º).

Conclusions
From the present results and discussion, the following conclusions can be drawn: 1.The determined atom polarizabilities are suitable for calculating reasonably accurate molecular polarizabilities of polyatomic molecules and, probably, have greater validity in general than additive values for calculating interactions among atoms in molecules.The anisotropies of a wide variety of polyatomic molecules can be largely attributed to atom dipole interactions, in support of Silberstein's [1-3] and Applequist's [4] suggestions.The variety of molecular properties, which can be accounted for, at least approximately suggests that the atom-dipole interaction model is a reasonable facsimile of real molecules.The model is appealing because it requires relatively little information for practical calculations.It is to be expected that further research on the model will deal to refinements and extensions to other properties; e.g., the theory has been extended to include effects because of the nonlinear response of atoms to an applied field, and to consider the role of higher multipole interactions in optical rotation.Both of these developments may well have significant effects on the properties discussed in this report.Even if one is interested in only the dipole response of the system, which is expressed by appropriate elements of the A -1 matrix, that response is influenced by all orders of multipole response, by virtue of the dependence of each element of A -1 on all elements of A.
2. A hypothesis is indispensable for a suitable representation of polarizability derivative, without which the registered variations will be of opposite sign to the experimental magnitudes.
3. The point-atom model is not only a mathematical approximation for real molecules, but also a well-defined system, whose behaviour can be calculated to a high degree of accuracy on the basis of known physical laws.Where the behaviour of a real molecule is similar to that of the model, an insight is gained into the origins of the behaviour.Where discrepancies, e.g., resonances, are found, it can be concluded that the molecule is not like the model.
4. The results of the present work clearly indicate that, because of the observed differences between the molecular polarizability results of our program POLAR and our version of program PAPID, it may become necessary to recalibrate POLAR.The values calculated with POLAR are quite larger than those obtained with the reference program PAPID, and POLAR polarizability values do not increase monotonically with bond distances.The recalibration means to obtain better atomic polarizabilities to start with in Equation ( 13).It appears that the results of good-quality ab initio calculations might be suitable as primary standards for such a calibration.
Work is in progress on the recalibration of POLAR.The difference between POLAR polarizability and PAPID polarizability results is because of the different parametrization schemes used for the initial atomic polarizabilities.PAPID uses atomic polarizability values fitted to high-quality calculations, while POLAR performs a simpler individual computation for each molecule, in order to exploit the difference among different atoms in different functional groups in different molecular environments.

Experimental Section
The dipole-dipole polarizability  is calculated with the interacting induced-dipole polarization model [1][2][3], which calculates effective anisotropic point polarizability tensors by the method of Applequist et al. [4].The molecular polarizability,  ab mol , is defined as the linear response to an external electric field: where  a ind is the induced molecular dipole moment and a, b, c… denote Cartesian components [9][10][11][12][13][14][15][16][17][18][19].Considering a set of N interacting atomic polarizabilities, the atomic induced-dipole moment has a contribution also from the other atoms: where T pq,bc (2) is the gradient tensor of the electric field: where r pq is the distance between atoms p and q, and  represents the Kronecker delta function: (a,b) = 1 if a = b, and (a,b) = 0 if a ≠ b.The molecular polarizability can then be written as: where  eff p is the effective polarizability of atom p and B is the relay matrix defined as (in supermatrix notation): The difference between our program POLAR and our version of program PAPID is in the different parametrization scheme used for the initial atomic polarizabilities.PAPID uses atomic  p values fitted to high-quality calculations, while POLAR performs a simpler individual computation for each molecule, in order to exploit the difference among different atoms in different functional groups in different molecular environments.The following improvements have been implemented in the polarization model: 1.In order to build the relay matrix B, the atomic polarizability tensor  p =  p  +  p  has been used instead of the scalar polarizability  p .The -atomic polarizability tensor is calculated according to: where the z-axis is defined as the bond direction for each bond.The diagonal form of   has two distinct components noted as   || and    , parallel and perpendicular to the bond axis, respectively.

The parameter  
|| /   = 1.676 has been obtained fitting the isotropic bonding polarizabilities of Vogel [20].The bonding polarizabilities have been implemented in the database of program SIBFA [21].The scalar -atomic polarizabilities have been evaluated with a modified Hückel molecular orbital (HMO) method.A wave-like scheme is assumed for the variation of the Hückel -parameter with the torsion angle.In first approximation,  between p z orbitals twisted from coplanarity by an angle  can be evaluated as  =  o cos, where  o is taken as the  value for benzene.An improved estimate of the  function is  =  o cos 1.15  [22].The -atomic polarizability tensor is calculated as: where the xy plane is the -plane.The   || and    are parallel and perpendicular to the -plane, respectively.The parameter    /  || = 1.741 has been obtained fitting the experimental polarizabilities of aromatic molecules.
2. Damping functions have been used in the calculation of the interaction tensor to prevent the polarizability from going to infinity [23][24][25].Thole proposed: where v pq = r pq /s pq if r pq < s pq ; otherwise v pq = 1 [26].The term s is defined as s = 1.662( p  q ) 1/6 .Miller proposed: with  pq =  p +  q (sum of the van der Waals radii) [27,28].
3. The interactions between bonded atoms and between atoms within a distance lying in an interval defined by [r inf ,r sup ] have been neglected.The starting values for this interval are [0, 10 30 ] a.u.The calculation can start from a predefined value of r inf (e.g., r inf = 0) and increment r inf by one atomic unit (r = 1 a.u.), each time resonance conditions are detected, up to r sup .
4. Four tests indicating a resonance condition have been implemented: test whether matrix B -1 is singular, whether matrix B -1 is not positive definite, whether matrix B is not defined positive and whether some effective matrix B p is not defined positive.
5. Vibrational corrections are obtained by a modification of atomic polarizability increments (Equation 8).
6.The optimized experimental polarization parameters have been included in the database for a number of atomic classes [29].

Figure 1 .
Figure 1.Polarizability of CH 4 as a function of polarizabilities of H and C [4].Units are Å 3 .Drawn with subroutine THREED (H.Jespersen, Iowa State University Computation Center).The vertical coordinate was truncated at 7 and -2.

Table 1
lists the vibrational corrections for CH 4 effective atomic and mean molecular polarizabilities.

Table 1 .
Vibrational corrections for methane effective atomic and mean molecular polarizabilities.