Derivation of force constants based on the electric field gradient

The present work devotes to studying the electric properties: electric quadrupole moment and electric field gradient (EFG) of molecules PdH and (PdH) based on the full relativistic theory. It is the first time to explore that the force constants are essentially caused by electric field gradients, and indirectly with spectroscopic data. If EFG is positive, zero or negative, then the 2 will be positive, zero or negative. Therefore, second order force constants are adjustable to changing the intensity of EFG. k


INTRODUCTION
The force constants of diatomic molecules, in general, are derived from spectroscopic data [1], or may be determined by a least-square fitting to the data of quantum mechanic calculations [2], and the latter is time tediously consuming and complicated.If the force constants and dissociation energy are known, the calculations of diatomic potential are numerically straightforward.
It is explored that the force constants are essentially connected with the electric field gradient (EFG), and indirectly with spectroscopic data.There are some works for the calculations of the electric field gradient using non-relativistic or relativistic quantum mechanic [3][4][5].The present work has evaluated the EFG of molecules PdH and (PdH) −2 based on the full relativistic quantum mechanic.

THE FULL RELATIVISTIC THEORY
Symmetry plays an important role for quantum mechanical theory.There are three levels for quantum me-chanical method: Non-relativistic based on the single point group, relativistic based on the double group [6,7] and relativistic based on the full symmetry group or quaternion symmetry [8,9].
In consideration of the equivalence of space and time, for the relativistic theory based on the full symmetry group or quaternion symmetry, the time reversal symmetry operation is included in the symmetry group.It is well known that the product of two space symmetry operators is represented by the product of the corresponding unitary matrices in the theory of group representation, however, time reversal symmetry operation is an antiunitary operator.The inclusion both of space and time operations, the group is called the full symmetry group, and their group representation is called a co-representation, this nomenclature is introduced by Wigner [8,9].The co-representation is not unitary representation; however, it is still possible to be broken down to irreducible form.The third level, i.e., relativistic quantum mechanical theory based on the full symmetry group or quaternion symmetry is applied by the work of T. Saue and the others.
Time reversal operator or called Kramer operator [7,10,11] is defined as follows If Hamiltonian  is a real function of r, then The commutator (3) is the necessary and sufficient condition for two operators  and having the same set of eigenfunctions, that is, the state is invariant to time reversal (motion reversal). It is proved that time reversal operator is antiunitary [10,11] It is also anti-linear operator.
Kramer theorem [7,8,10,11]: It can be proved that corresponding to J value of half odd integer, i.e. fermions, there will be the new double degeneracy by time reversal; and corresponding to J value of integer, i.e. bosons, there will be no double degeneracy, where identity.
For the time-independent Dirac equation where where, 0 And Pauli spins matrices

If
, it is the Dirac equation of free electrons.ˆ0 V  Now, time reversal operator is defined using another equivalent method, i.e.Kramer theorem [8,9], in which , that is operator [8,9] is defined as follows where It is easy to know That is, A is a Hermitian and B is antihermitian.For Kramer's partners, Dirac operator ˆD h can be expressed as where, ˆD h is expressed as the quaternion, which exhibits the time reversal symmetry.
Quaternion algebra [12] was developed by Hamilton and Fresenius, however, limited applications of quarternion in quantum mechanics have been made in recent years.
Quaternion is expressed as ; ; where , , are quaternion units, i-imaginary, 0 , 1 , 2 and 3 V are real part, Pauli spins matrices in (7).Quaternion includes 3-dimension complex space and one dimension real space.It is known from (7) to (11) that quaternion algebra includes both time and space reversal symmetry, then, it is called the full symmetry group, its matrix representation is called co-representation.The ˆD h in (12) is the Dirac operator of full symmetry group.The Dirac equation can therefore be rewritten into quaternion algebra form q q q D h     (15)

THE CALCULATIONS OF ELECTRIC PROPERTIES
The calculations of energy and electric properties: electric quadrupole moment, electric field gradients of molecules PdH based on the full relativistic quantum mechanical theory with the basis dyall.v2zusing of program DIRAC10.The results are shown as Tables 1 and  2.
The charge distributions of dipole for molecule PdH is calculated as follows      Where r is the nuclear distance and is its equilibrium nuclear distance, and 2 , 3 and 4 k are the quadratic, cubic and quartic force constants, respectively.The negative derivative of potential with respect to bond distance r is the electric field force of nuclei, therefore e r k k where E is the electric field strength, q is the electric charge.
However, its quadrupole is the second-rank tensor with two dimension charge distribution (vertical plane V  ), Figure 1 in which the unit of charge distributions is 10 19 C, C-coulomb, If r = r e , then, the quadratic force constant will be The zz q and   xx yy q of PdH can be calculated as where q is negative and take its absolute value.If EFG is positive, zero or negative, then, the 2 will be positive, zero or negative.If we want the 2 to be zero or negative, it is possible to make the EFG descended or tend to negative.
Suppose there is 1 au of electric charge and at 1au of electric field gradient, the quadratic force constant 2 k For neutral molecule PdH, i.e., electron number of 47, and the = 2.059396 for Pd and = 0.184936 for H, and their Mulliken Atomic Charges: Pd = 0.037613, H = −0.037613.Then zz q zz q OPEN ACCESS q zz = +2.144107q zz = +2.144107q xx(yy) = −1.0721 where 1.97719645 is the reduced mass of PdD.All the electric field gradients (EFG), quadratic force constants 2 and Harmonic frequencies k  for mole- cules are listed in Table 2.

CONLUSIONS
From Equation (19), it is proved that the second order force constant is directly derived from electric field gradients EFG.Therefore, it is the first time to explore that the force constant is essentially caused by electric field gradients, and the second order force constant and its harmonic oscillator frequency have been derived from electric field gradients.However, the force constants are usually derived from spectroscopic data; however, it is indirect with spectroscopic data.Therefore, this is the sup-plemental to the spectral values.
Momentously, the second order force constants are adjustable to changing the intensity of EFG.If EFG is positive, zero or negative, then, the 2 will be positive, zero or negative.It is obvious to see from The linear charge distributions, ±0.73647/10 19 C, of dipole moment for molecule PdH is much less than its two-dimension charge distribution of quadrupole moment q zz = 2.144107/10 19 C and q xx(yy) = 1.0721/1019 C. The charge distribution of PdH could be comparable with that of molecules HBr or HCl; however, it does not like HF, to see Table 3 of references [13], and therefore quadrupole  , in which the state are recovered by twice of time reversal, however, the wave function is changed to the negative; If , in which the state are recovered by twice of time reversal, the wave function is unchanged.Then  and  are said to form a Kramers pair.Therefore, the operator can be expanded in terms of Kramer's partners    and    .For example, Hermitian of one-electron

Table 1 .
The energy, electric dipole and quadrupole moment of molecule PdH.

Table 2 .
The electric field gradients (EFG), k 2 and  for molecules.  The quadrupole charge distribution of PdH.