Ab-Initio Calculations of 27 Electronic States of the BP + Ion-Molecule

The theoretical investigation of the potential energy curves, in the representation 2s+1Λ(+/−), of the 27 low-lying Doublet and Quartet electronic states of the BP+ molecular ion has been performed with the methods in quantum chemistry, the Complete Active Space Self Consistent Field (CASSCF) and the Multireference Configuration Interaction (MRCI) calculations. The harmonic vibrational frequency ωe, the inter-nuclear distance at equilibrium Re, the rotational constant Be, the electronic energy with respect to the minimum ground state energy Te, and the permanent dipole moment have also been calculated. Twenty-three new electronic states have been investigated here for the first time. The comparison between the values of the present work and those available in the literature for several electronic states shows a good agreement. These investigated data can be a conducive to further work on BP+ molecular ion in both experimental and theoretical research.


Introduction
The theoretical calculations and potential energy function of diatomic molecules have been one of several molecular physics important researches [1].Although a solid description of molecular properties of molecular potential energy function is approximated at Born-Oppenheimer, it remains important for the analysis of dissociation and stability [2].In 1972, Gingerich and Piacente [3] only studied the dissociation energy of BP molecule by means of Knudsen cell mass spectrometry; they measured the dissociation energy ( ) 0 0 D BP = 82.0& 4 kcal mol −l or, and 0 K was used as the reference temperature, but they did not give other experimental data on the BP radical.BP + ion-molecule, attracted much attention as refractory semiconductor due to its physical and chemical properties [4]- [11].Linguerri et al. [11] found that the ground state of BP + was X 4 Σ − and they determined their spectroscopic constants T e , R e , ω e , ω e χ e and B e for the first time by using MRCI calculations.The physical properties of BP + took less attention than the BP molecule [6]- [9], but recently the ground state and the lowest excited state of BP + ion-molecule acquired much focus.Multi-reference configuration interaction method was applied by Guo et al. [12] to investigate the BP + ground and the lowest 5 excited electronic states.
The present work is devoted to an accurate description of the ground and electronically 26 excited doublet and quartet electronic states.In order to emphasize the accuracy of our work we used the multi-reference configuration interaction (MRCI+Q) model expansion.The potential energy curves, the spectroscopic constants R e , T e , B e , and ω e for each of the corresponding electronic states have been investigated along with their static dipole moment.

Method of Calculations
The Hartree-Fock method relies on averages; it does not consider the instantaneous electrostatic interactions between electrons, nor does it take into account the quantum mechanical effects on electron distributions because the effect of n-1 electrons, on an electron of interest, is treated in an average way.This method ignores the electron correlation which is the difference in energy between the HF and the lowest possible energy.
A singly, doubly, triply, quadruply, …, N-tuply excited states is called a Configuration State Function (CSF).The exact ground and excited-states wavefuntions can be expressed as a linear combination of all possible Nelectron Slater determinants arising from a complete set of spin orbitals where all the C's are expansion coefficients that must converge to zero and this method is called Configuration Interaction (CI).A simple formula known as the Davidson Correction has been proposed for estimating the contribution ΔE Q of quadruply excited determinants to the correlation energy The development of efficient multiconfiguration self-consistent field method MCSCF is particularly important for excited states.One such scheme is the Complete Active-Apace Self-Consistent Field method (CASSCF).In this approach, the spinorbitals are divided into three classes: 1) set of inactive orbitals composed of the lowest energy spinorbitals which are doubly occupied in all determinants, 2) a set of virtual orbitals of very high energy spinorbitals which are unoccupied in all determinants, 3) a set of active orbitals which are energetically intermediate between the inactive and virtual orbitals.In the conventional CI methods, the Hartree-Fock SCF wave function Ψ 0 is used as a reference configuration and the configuration state functions are formed by moving electrons out of the occupied spinorbitals of Ψ 0 into unoccupied spinorbitals.In multireference configuration interaction (MRCI) a set of reference configurations is created from which excited determinants are formed for use in a CI calculation.Large scale MRCI wavefunctions can generate very accurate wavefunctions but are also computationally very intensive.

Results and Discussion
The potential energy curves of the 26 doublet and quartet electronic states, in the representation 2s+1 Λ ± of the molecule BP + , has been performed for 189 internuclear distances.These curves are given in Figures 1-4.
One can notice some avoided crossings between different doublet and quartet electronic states.The repulsion between the two atoms, when they are brought close to each other, is related to the Pauli principle: when the electronic clouds surrounding the atoms start to overlap, the energy of the system increases abruptly.Moreover, some crossings, of abscissa R c , and avoided crossings, of abscissa R ac , have been obtained between the potential energy curves of Figures 1-4.If the two curves correspond to states of different symmetry cross, they will ignore each other.Their wavefunctions Φ 1 and Φ 2 remain unperturbed and they are the adiabatic solutions of the Schrödinger equation.The crossing in this case is strictly allowed.But if these wavefunctions have the same symmetry, they will mix with each other to give two adiabatic solutions which no longer cross and the crossing    becomes avoided.The adiabatic solutions of the Schrödinger equation Ψ 1 and Ψ 2 are obtained by linear combinations of the diabatic ones where the variation method is used to solve for the coefficients.If H 12 is the matrix element of the electronic hamiltonian H of the system between Φ 1 and Φ 2 , the energy gap E ac between the two adiabatic states is twice the coupling term H 12 at R ac .The calculated values of R ac and E ac are given in Table 1.Such crossings or avoided crossings can dramatically alter the stability of the molecule.
The spectroscopic constants such as the vibrational harmonic constant ω e , the internuclear distance at equilibrium R e , the rotational constant B e , and the electronic transition energy with respect to the ground state T e have been calculated by fitting the energy values around the equilibrium position to a polynomial in terms of the internuclear distance by using the data obtained from the ab initio calculation.These values are given in Table 2 along with the rare available data in literature on this ion-molecule.The comparison of our calculated data with those in literature for the 4 electronic states X 4 Σ − , (1) 2 Σ − , (1) 2 Δ and (1) 2 Π one can find that an excellent agreement is obtained for the values of R e with the relative difference 0.5% < ΔR e /R e < 1.0% and a very good agreement for the values of B e with the relative difference 1.1% < ΔB e /B e < 4.9% with a larger relative difference for the state (1) 2 Π with the relative difference 8.2%.By comparing our calculated values of ω e with those in literature we found a very good agreement with relative difference 0.3%< ΔB e /B e < 4.0% for the 3 states X 4 Σ − , (1) 2 Σ − , (1) 2 Δ, while the agreement deteriorate for the state (1) 2 Π with the relative difference 23%.For our values of T e the relative difference is good for the (1) 2 Σ − with the relative difference 5.4% and becomes lager for the states (1) 2 Δ and (1) 2 Π with the relative differences 15% and 27% respectively.No comparison with other investigated states since, up to our knowledge, they are given here for the first time.

Permanent Dipole Moment
The permanent dipole moments are investigated for the considered electronic states of the ion-molecule BP + by taking the boron atom at the origin, and the positions of the phosphorus atom are variable along the positive internuclear z-axis.All the calculations were performed with the MOLPRO [21] program.The dipole moment operator is among the most reliably predicted physical properties, because the quantum mechanical operator is a simple sum of one-electron operators.The expectation value of this operator is sensitive to the nature of the least energetic and most chemically relevant valence electrons.The positive sign of the dipole moment corresponds to a charge transfer from the P atom towards the B atom.To obtain the best accuracy, MRCI wavefunctions were constructed using MCSCF active space.The values of the permanent dipole moments, for the investigated electronic states, are given in atomic unites (a.u.) as a function of the internuclear distance R in Figures 5-8.One can notice that: 1) some of dipole moment curves, as (1) 2 Σ − , (2) 2 Σ + , (1) 2 Δ, (1,2,3) 2 Π, (1) 4 Σ − , (3) 4 Π, tend to ± ∞ when R tends the ∞ which explains the ionic bounding of all these in these regions, 2) some other dipole moment curves, as (2) 4 Σ − , (1,2,4) 4 Π, (1) 4 Δ, (1) 4 Σ + , (3) 4 Π, tend to constant limit where the bounding is covalent, 3) the agreement between the positions of the avoided crossings of the potential energy curves of 2 electronic states and the positions of the crossing of the dipole moment curves of these 2 states.This agreement may confirm the validity and the accuracy of the present calculation.

Conclusion
In the present work, ab initio calculation for 27 low-lying electronic states, in the representation 2s+1 Λ (+/−) , of the BP + ion-molecule were computed via CASSCF/MRSDCI methods for the electronic excited The potential energy curves, the electronic energy with respect to the ground state T e , the harmonic frequency ω e , the internuclear distance R e have been calculated along with the rotational constant B e .The comparisons of the present results with the available values in the literature show an overall very good agreement.To the best of our knowledge, 23 new electronic states have been investigated for the first time through this work.These new excited electronic states may provide a reliable theoretical basis and information for the experimental spectral properties related to the electronic structure for BP + ion-molecule and the relevant sub-system [22] [23] study as well.

Figure 1 .
Figure 1.Potential energy curves of the doublet electronic 2 Σ ± and 2 Δ states of the ion-molecule BP + .

Figure 2 .
Figure 2. Potential energy curves of the doublet electronic 2 Π and 2 Φ states of the ion-molecule BP + .

Figure 3 .
Figure 3. Potential energy curves of the quartet electronic 4 Σ ± and 4 Δ states of the ion-molecule BP + .

Figure 4 .
Figure 4. Potential energy curves of the quartet electronic 4 Π and 4 Φ states of the ion-molecule BP + .

Figure 5 .
Figure 5. Dipole moment curves of the doublet electronic 2 Σ ± and 2 Δ states of the ion-molecule BP + .

Figure 6 .
Figure 6.Dipole moment curves of the doublet electronic 2 Π states of the ion-molecule BP + .

Figure 7 .
Figure 7. Dipole moment curves of the quartet electronic 4 Σ ± and 4 Δ states of the ion-molecule BP + .

Figure 8 .
Figure 8. Dipole moment curves of the quartet electronic 4 Π states of the ion-molecule BP + .

Table 1 .
Avoided crossings between the different electronic states of the ion-molecule BP + .

Table 2 .
Spectroscopic constants of the potential wells of the electronic states of the BP + molecule.