Next Article in Journal
Pharmaceutical Potential of Synthetic and Natural Pyrrolomycins
Next Article in Special Issue
Thermal Conductivity of Epoxy Resin Composites Filled with Combustion Synthesized h-BN Particles
Previous Article in Journal
Influence of Excipients and Spray Drying on the Physical and Chemical Properties of Nutraceutical Capsules Containing Phytochemicals from Black Bean Extract
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetry Breaking of B2N(−, 0, +): An Aspect of the Electric Potential and Atomic Charges

Department of Chemistry, Science and Research Branch, Islamic Azad University, Tehran 801, Iran
*
Author to whom correspondence should be addressed.
Molecules 2015, 20(12), 21636-21657; https://doi.org/10.3390/molecules201219769
Submission received: 6 August 2015 / Revised: 16 October 2015 / Accepted: 29 October 2015 / Published: 3 December 2015
(This article belongs to the Special Issue Boron Nitride: Synthesis and Application)

Abstract

:
In this study, the three forms of B2N(−, 0, +)—radical, anion and cation—have been compared in terms of electric potential and atomic charges, ESP, rather than the well-known cut of the potential energy surface (PES). We have realized that the double minimum of the BNB radical is related to the lack of the correct permutational symmetry of the wave function and charge distribution. The symmetry breaking (SB) for B2N(0, +) exhibits energy barrier in the region of (5–150) cm−1. The SB barrier goes through a dynamic change with no centrosymmetric form which depends on the wave function or charge distribution. In spite of A ˜ 2 Σ g + exited state, the   B ˜ 2 g excited configuration contributes to the ground state (   B ˜ 2 g X ˜ 2 Σ u + ) for forming radicals. The SB did not occur for the anion form (B2N(−)) in any electrostatic potential and charges distribution. Finally, we have modified the Columbic term of the Schrödinger equation to define the parameters “αα' and ββ'” in order to investigate the SBs subject.

1. Introduction

Linear triatomic structures, particularly those of the X-Y-X type are one of the simplest molecular systems with “high” symmetry in which symmetry breaking (SB) may occur. The point group of the BNB in the linear symmetric X-Y-X species is “Dh” whereas its asymmetrically distorted form has the symmetry of the Cv group. Tri-atomic B2N(−, 0, +) molecules have been subsequently studied using a variety of calculations and spectroscopic methods. It is a deep challenge to measure the real or artifactual SB effects due to its capability to a display a pseudo second-order Jahn–Teller effect, which results in a structure with unequal BN bond lengths [1,2,3,4,5,6,7,8,9,10,11,12,13,14].
Obviously, the adiabatic system of the Born–Oppenheimer approximation breaks down when the electronic states are degenerate, so the SB can be associated with degeneracy. SB may result from first-order Jahn–Teller effects or Renner–Teller effects for degenerate electronic states, while non-degenerate states result from interactions between different states, i.e., the so-called pseudo Jahn–Teller (SOJT) effects.
The first important investigation for BNB was done by Martin et al. in 1989. Based on an UHF/6-31G* geometry optimization, they predicted that the B2N has a symmetric linear regulation in its ground state ( X ˜ 2 Σ u + ) ) with an unusually low bending frequency (73 cm−1) [1].
They noted that the full valence CASSCF wave function for the X ˜ 2 Σ u + state of linear B–N(0)–B is unstable with respect to symmetry lowering, i.e., the Cv structure (r1(BN)r2(BN)) yields a lower total energy than the Dh symmetric structure (r1(BN) = r2(BN)) [10].
In 1992, Knight et al. found two minimum structures, i.e., the linear B2N(0) and cyclic B2N [4]. The UHF/6-31G* level of theory predicts the cyclic 2B2 state to be at the global minimum, while the correlated methods predict the X ˜ 2 Σ ( u ) + state of linear B–N–B to be at the global minimum and the cyclic 2B2 state to be at the local minimum.
There are two conceivable valence bond structures for open-shell system of the BNB radical in the ground state. The first one is the localized form where the symmetry is Cv{2Σ+:1σ2, 2σ2, 3σ2, 4σ2, 5σ2, 1π4, 6σ2, 7σ1}, and the second is the resonance form where the symmetry is Dh { ( 2 Σ u + ) : 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 1 π u 4 , 4 σ g 2 , 3 σ u 1 }. Obviously describing the electronic structure of the open-shell system (B2N(0)) is much more difficult than the closed-shell system of B2N(−, +) and in many cases, the description of an open-shell system is challenging in computational quantum chemistry calculations. This is primarily due to the presence of static correlation effects (requiring a multireference-type description) [13,14].
In general, researchers agree that BNB is linear in its ground electronic state of ( X ˜ 1 Σ g + ) and ( X   ˜ 2 Σ u + ) for anionic and neutral forms, respectively [2,13].The ground state of the BN(0)B has been examined via the developed reduced multireference coupled cluster method with singles and doubles that is perturbatively corrected for triples [RMR CCSD (T)] using the correlation consistent basis sets (cc-pVDZ, cc-pVTZ and cc-pVQZ) by J. Paldus [9]. They showed that the ground state has an asymmetric structure Cv with two BN bonds of unequal length.
In an experimental study, the state of A   ˜ 2 Σ g + ( B 2 N ) was observed in photoelectron spectroscopic studies (PES) by Asmis et al., placing the zero-point level of the 6330 ± 40 cm−1 above the ground state of the X ˜ 2 Σ ( u ) + [2]. He showed that the observed signal in the 355 and 266 nm photoelectron spectra of B2N has been indicated to a photodetachment from the anion ground state ( X ˜ 1 Σ g + ) to the ground and lowest excited states of neutral B2N i.e., X ˜ 2 Σ u +   and   A ˜ 2 Σ g +   with a linear symmetry and is assigned to the X ˜ 1 Σ g +   X ˜ 2 Σ u + +   e and X ˜ 1 Σ g +   A ˜   2 Σ g + +   e transitions.
Electron spin resonance (ESR) experiments indicate the unpaired electron between the two B atoms implying symmetric geometry [4]. However, this result cannot exclude the presence of a small barrier to asymmetry in point view of the time scales involved. A further, matrix infrared study [7] indicated that the ground state of BNB is cyclic, whereas a more recent photoelectron spectroscopic study of B2N(−) anion [2] confirmed a linear symmetric geometry. Furthermore, the infrared absorptions were also observed in the cryogenic argon matrix near 6000 cm−1 of the electronic band system, due to the A   ˜ 2 Σ g + X   ˜ 2 Σ u + [7].
Recently, the compounds of BNB, especially their electronic structures have been considered [15,16,17,18,19]. In 2009, a series of multi-reference approaches based on the SA-CASSCF wave function, i.e., CASPT2, MRCI, and MRAQCC, have been employed by Boggs and coworkers to investigate the SB in the ground state X ˜ 2 Σ ( u ) + of the triatomic B2N(0) radical [20]. Their results show that B2N in its ground state has a linear non-centrosymmetric structure with two equivalent global minima of the adiabatic potential energy surface, including two oppositely directed dipole moments, respectively. They accepted that the PJT effect involving vibronic interaction with the first excited state A ˜ 2 Σ ( g ) + via the asymmetric stretching vibrations is the major reason for the double-minimum. On the other hand, the large-scale multi-reference configuration interaction calculations, CASSCF+1+2 predicted an asymmetric configuration, while the SACASSCF+1+2 predicted a symmetric Dh ground state [21].
In another study in 2010, Stanton has discussed an unusually large non-adiabatic error in the BNB molecule as well as non-adiabatic corrections to energy level [22]. He illustrated those non-adiabatic corrections to energy levels should fall out only when the affected vibrational frequency is large enough to be of comparable magnitude to the energy gap. In other words, non-adiabatic corrections should be given as much weight as issues such as high-level electron correlation, relativistic corrections, etc.
In other words, calculations using larger and larger basis sets, and more and more advanced methods of electron correlation, are doomed to approach the “wrong” limit for the vibronic levels of BNB if the Born-Oppenheimer approximation is applied, so the electrostatic potential charges which are based on wave functions of the BNB systems and are related to quantum mechanics phenomena can be used as suitable tools through the significant approach reported in this work.
In 2013, Kalemos [21] tried to approach the SB problem by using high level multi-reference variation and full configuration interaction methods. He indicated that the (SB) problem is related to the lack of the “correct” permutation symmetry of the wave-functions adopted to attack the problem and is by no means a real effect. Furthermore, he indicated the wave-function which is not invariant under all symmetry operations of the point group should be symmetrically broken (SB). SB (in classical mechanics) occurs when a stable minimum undergoes splits into two stable minima. He checked the MRCI results by FCI (9e)/[3s2p] calculations and found no trace of SB in qualitative disagreement with all previous theoretical investigations that predicted a barrier to a centro-symmetric structure either of 20 cm−1 (based on MRCI methods of Boggs) or of 100–160 cm−1 (based on CC methods) [9,22,23].
In this study, we compare three forms of B2N(−, 0, +) (anion, radical and cation) in terms of electrostatic potential charges “ESP” rather than the cut of the potential energy surface (ESP is changed using the trial wave functions). For a charged system with charge Q, the density | ψ ( x ) | 2 multiplied by the atomic charge yields the charge density | Q ( x ) | 2 . Large points for fitting of various situations have been used to calculate the atomic charges and electrostatic potential of the systems. As a result, the possibility of an asymmetric ground state may not be eliminated or it seems that an asymmetric geometry is a rather comfortable situation; the double minimum nature of BNB is related to the lack of the correct permutation symmetry of the wave functions (Scheme 1).
Scheme 1. A simplified of spontaneous symmetry breaking, (a) indicates high energy level which the ball settles in the center, and the result is symmetrical; (b) and (c) are lower in energy levels and the overall “rules” remain symmetrical, however, as the potential comes into effect, the “local symmetry” is inevitably broken since eventually the ball must roll one way (at random) and not another; (A) is a non-linear combination of two wave function and (B) is a non-linear combination of three wave function.
Scheme 1. A simplified of spontaneous symmetry breaking, (a) indicates high energy level which the ball settles in the center, and the result is symmetrical; (b) and (c) are lower in energy levels and the overall “rules” remain symmetrical, however, as the potential comes into effect, the “local symmetry” is inevitably broken since eventually the ball must roll one way (at random) and not another; (A) is a non-linear combination of two wave function and (B) is a non-linear combination of three wave function.
Molecules 20 19769 g005
The small barrier energy including the two equal minima in some parts of B2N(−, 0, +) systems have been estimated using symmetrical linear combination of wave functions (SLC-WFs). These barriers are dependent on the charge distribution, SLC-WFs and correlation effects between | α   or   | β levels which are not global minima.
There are 16, 17, and 18 electrons in B2N(−, 0, +) systems. Therefore, based on Walsh prediction [24], it appears that the linear p-block molecule holds up quite well for all B2N(−, 0, +) ions and the radical. This work has focused on a spontaneous symmetry breaking (SSB) [25,26,27] for B2N(−, 0, +) systems in view of ESP.

2. Wave Functions and Symmetry Breaking

SB is generally applied to two phenomena—the first is the failure of the electronic wave function to transform as an irreducible representation of the molecular point group and the second is the preference of the nuclear framework for lower-symmetry geometry [28]. The failure of the electronic wave function is purely artifactual in that the exact wave function necessarily obeys the symmetry properties of the molecular point group and the symmetry-contaminated wave functions was described by Lowdin as the symmetry dilemma [29]. In such cases, a broken-symmetry wave function may yield a solution with lower energy than the symmetry-adapted one.
On the other hand, in the absence of external fields, the motion of nuclei is governed by a function of nuclear coordinates which is a sum of the Coulomb repulsion between nuclei and the electronic effective potential terms, resulting from the Born Oppenheimer approximation and wave function, so the meaning of the wave function can be interpreted by the charge density distributions of electrons in a quantum system [30].
It is obvious that every wave function must be satisfactory before being submitted to the Schrödinger equation; however, it has been shown by Wigner that exact electronic wave functions satisfy the Pauli Exclusion Principle [31]. The physical root of charge density is one of the main questions in electronic wave function origin and the key to unveiling the meaning of the wave function is to find the physical root of the charge density [32,33,34,35].
Many questions have been raised from interpretation of the wave function in the theory which is dependent on SB problem such as regarding the wave function as a field similar to electromagnetic field [36], an active information field [37], and a field carrying energy and momentum [38].
In a time independent system using Born-Oppenheimer approximation, the Hamiltonian of boron and nitrogen nuclei are sum of the kinetic energy term and effective potentials (Equation (1)):
H ^ n = j = 1 3 2 2 M j Δ j + V ( R B 1 , R N , R B 2 )
The B2N(−, 0, +) forms will be stable by minimization of “V”. So, it depends on sum of Coulomb repulsions between nuclei and the electronic effective potential:
V ( R ^ B , R ^ N , R ^ B ) = j = 1 2 l = j + 1 3 q j q l 4 π ε 0 | R j R l |   + E ( R B 1 , R N ,   R B 2 )
The dynamic representation was first postulated by Wigner in 1930 with the justification that its action does not change the relative distances between nuclei and is therefore the symmetry of the dynamical problem [31]. Damljanovic modified the Wigner postulate in another way for dynamical representation, while for every molecule; there exists at least one normal mode that belongs to the totally symmetric irreducible representation of the point group of that molecule [34]. In this approach, the problem of finding minima of V becomes another case of the spontaneous symmetry breaking phenomenon and it allows finding approximate relations between bond lengths in a molecule and its vibrational frequencies. The Schrödinger equation for electrons moving in the field generated by nuclei in a molecule subjected to no external fields can be written as follows:
[ j = 1 N e 2 2 m e . ( r j a ) 2 j = 1 N e l = 1 N e q ( l ) 4 π ε 0 | r j a R j a | +   j = 1 N e 1 l = j + 1 N e e 2 4 π ε 0 | r j a r l a | ] ψ ( r a , R a ) = ψ ( r a , R a )
where r a = ( r 1 a ,   r N e a ) T and R a =   ( R 1 a ,   R N e a ) T and T is the transposition vector in the 3N-dimesional real vector space. In this equation, the first term in the left is kinetic energy of the electrons, the second term is Coulomb attraction between electrons and nuclei and the third term is Coulomb repulsion between electrons and coordinates of nuclei that are the parameters on which the eigenvalue E is dependent “E = E( R a )”. Damljanovic has discussed that any two configurations of nuclei can be obtained from each other by rotation, translation or a permutation of one configuration as a whole correspond to the same value of V or E [39].
Since Coulomb repulsion between nuclei diverges at origin, the function “V” cannot be expanded by the use of Taylor formula around that point. On the other hand, “E” is finite at the point R = 0. It is equal to the electronic energy of an atom (called united atom) having charge of the nucleus equal to sum of charges of all nuclei in the molecule under investigation. Moreover, stable configurations of relatively simple molecules are confined in space in the small volume around an origin.
We have modified our systems based on definition of “V” in order to find minima as follows:
V = α [ ( R N R B 1 ) 2 ] + β [ ( R N R B 2 ) 2 ] + γ [ ( R B 1 R B 2 ) 2 ] + 1 4 π ε 0 [ α q N q B 1 | R N R B 1 | + β q N q B 2 | R N R B 2 | + γ q B 1 q B 1 | R B 1 R B 2 | ]
where R N , R B 1 ,   R B 2 and charges qN, q B 1 , q B 2 are bound distances and charge distribution of B2N(−, 0, +) and λ 1 , λ 2 and λ 3 are defined as:
λ 1 = R N R B 1 ,   λ 2 = R N R B 2  and  λ 3 = R B 1 + R B 2 + R N  and V  = α λ 1 2 + β λ 2 2 + γ ( λ 1 λ 2 ) 2 + 1 4 π ε 0 [ α q N q B 1 | λ 1 | + β q N q B 2 | λ 2 | + γ q B 1 q B 1 | λ 1 λ 2 | ]
λ 1 = λ 2 = λ 0 , in the center of mass of BNB molecules, while:
V ( λ ) = 2   α | λ 2 | + 2 4 π ε 0 α q N q B | λ |  and in minimum : 4 α | λ 0 | = 2 4 π ε 0 α q N q B | λ 0 | q N q B = 8 π ε 0 α α | λ 0 3 |
Depending on the radical, cation and anion forms of BNB, qN × qB multiplication can be positive or negative and consequently, ⎢ λ 0 = q N ( 0 ) q B ( 0 ) 8 π ε 0 α α 3 would be either positive or negative.
In our calculations, the charges for nitrogen were changed between +0.75 < qN <+0.87 (for the anion form) and in equilibrium the charges of the nitrogen and the two borons are qN = 0.868802 and q B 1 = q B 2 = q B = 0.934401 . Therefore, the λ 0 = 1.329 and αα' = ββ' = 1.506 can be yielded in the center.
λ 1 λ 2 ,   λ 1 ( 0 ) and λ 2 ( 0 ) denote a stable configuration, and V can be expanded around this point up to the second order as Equation (10):
1 | λ + δ λ | 1 | λ | λ ( 0 ) λ | λ 3 | 1 2 ( δ λ ) 2 | λ 3 | +   3 2 ( λ ( 0 ) δ λ ) 2 | λ 5 |
{ ( 2 α + 2 β + 2 γ α q N q B 1 4 π ε 0 | λ 1 ( 0 ) 3 | γ q B 1 q B 2 4 π ε 0 | ( λ 1 ( 0 ) λ 2 ( 0 ) ) 3 | ) } ( λ 1 ( 0 ) ) +   { ( 2 γ γ q B 1 q B 2 4 π ε 0 | ( λ 1 ( 0 ) λ 2 ( 0 ) ) 3 | ) } ( λ 2 ( 0 ) ) = 0
{ ( 2 γ γ q B 1 q B 2 4 π ε 0 | ( λ 1 ( 0 ) λ 2 ( 0 ) ) 3 | ) } ( λ 1 ( 0 ) ) +   { ( 2 α + 2 β + 2 γ β q N q B 2 4 π ε 0 | λ 2 ( 0 ) 3 | γ q B 1 q B 2 4 π ε 0 | ( λ 1 ( 0 ) λ 2 ( 0 ) ) 3 | ) } ( λ 2 ( 0 ) ) = 0
so, the distance between B and N is equal to | λ B N | = q N ( n ) q B ( n ) 8 π ε 0 α α 3 = q N ( n ) q B ( n ) 8 π ε 0 β β 3 and the distance between the two borons is | λ B B | = q B 1 ( m ) q B 2 ( m ) 8 π ε 0 γ γ 3 .
In the B2N(−, 0, +) radical, cation and anion forms, the charges of atoms always localize in a definite position in space. In fact, for a charged quantum system, it has been described by the wave function. Thus, the charges distribution with a certain amount in space and different distributions between borons and nitrogen atoms are important for the understanding of real or artifactual SBs problems of radical and ion BNB forms.

3. Computational Details

Even though various methods and basis sets (both large and medium) have been employed in this study, among them, the EPR-III and EPR-II basis sets of Barone [40] show accurate results for electrostatic potential (ESP) fitting. EPR-II is a double-ζ basis set with a single set of polarization functions and an enhanced “s” part: (6, 1)/ [41] for H and (10, 5, 1)/ [41,42] for B to F. EPR-III is a triple-ζ basis set including diffuse functions, double d-polarizations, and a single set of f-polarization functions. Also in this case, the s-part is improved to better describe the nuclear region, (6, 2)/[37,38] for H and [40,41,42,43,44] for B up to F.
The active space for the CASSCF methods was composed of all valence electrons and orbitals of these atoms, i.e., 11 active electrons and 12 active orbitals for B2N(0) and 10 and 12 electrons for B2N(+)and B2N(−) respectively. In some part of our discussion, the BNB has been optimized via various levels of theory such as CASSCF (11, 12)/cc-pvqz and CASSCF (11, 12)/AUG-cc-pvqz (for radical) and CASSCF (10, 12)/cc-pvqz for cation. Approximation spin orbit coupling between two spin states has been computed during CASSCF calculations [45,46].
A Quadratic CI calculation including single and double substitutions [47] has been used to evaluate various one-electron properties including NBO, bonding analysis, atoms in molecules (AIM) [48], multi-pole moment, natural population analysis, electrostatic potentials, and electrostatic potential-derived charge using the Merz-Kollman-Singh [49], chelp [50], or ChelpG [51].
Polarizabilities and hyper-polarizabilities have been computed by CISD and QCISD and CASSCF methods and double numerical differentiation of energies have been used by with the pol=En only keyword in some cases. The AIM properties have been used to request molecular properties predicted via the theory of atoms in molecules [48]. The AIM keyword is used to compute atomic charges of atoms in molecules, covalent bonds, localized orbitals, and critical points.
Atomic charges have been calculated from electrostatic potentials using a grid-based method, ChelpG, which was developed by Breneman and Wiberg [51]. Atomic charges are fitted to reproduce the molecular electrostatic potential (MESP) at a number of points around the molecule [52,53]. The charge calculation methods based on molecular electrostatic potential (MESP) fitting (including CHELPG) are not well-suited for treating larger systems where some of the innermost atoms are located far away from the points at which the MESP is computed. In such a circumstance, variations of the innermost atomic charges will not lead to significant changes of the MESP outside of the molecule, meaning that the accurate values for the innermost atomic charges are not well-determined by MESP outside of the molecule [52,53].
The representative atomic charges for flexible molecules hence should be computed as average values over several molecular conformations. A number of alternative MESP charge schemes have been developed, such as those employing Connolly surfaces or geodesic point selection algorithms, in order to improve the rotational invariance by increasing the point selection density and reducing anisotropies in the sampled points on the MESP surface. A detailed overview of the effects of the basis set and the Hamiltonian on the charge distribution can be found in reference [54]. CHELPG charges can be computed using the popular ab initio quantum chemical packages such as Gaussian or GAMESS-US.
Indeed in our study, because of the large number of calculations in various ESP simulation situations, using expensive basis sets and methods such as MRCI was difficult and also not necessary for CHELPG and ESP calculations. Therefore, with medium methods in terms of computational cost, we have obtained accurate results for our approach. All the calculations were performed using the Gaussian program package [55] and the optimization was done along with a frequency calculation to verify that the geometry was a real minimum without any imaginary frequency.

4. Results and Discussion

The radical form of BNB is linear in its ground electronic state ( X   ˜ 2 Σ u + ) ) with an orbital occupancy of 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 1 π u 4 , 4 σ g 2 , 3 σ u 1 while the lowest electronically excited state is predicted to be A ˜ 2 Σ ( g ) + , with an orbital occupancy of 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 1 π u 4 , 4 σ g 1 , 3 σ u 2 .
Although the geometric data in Table 1 show boron-nitrogen bonds within a molecule that differ from each other by less than 0.001 angstrom, the BNB in the three radical, cation and anion forms in the ground state are linear. Moreover, the other structures have angles that differ from 180 degrees by less than a degree or in many cases less than a tenth of a degree and the reason depend on the methods used.
In addition, the   B ˜ 2 g exited state with an orbital occupancy of 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , π u 4 , 4 σ g 2 , 1 π g 1 (above the A   ˜ 2 ) is subject to the Renner-Teller effect and further exited states depend on the ( 4 g )   of triplet form. The geometric structures and electronic energies in ground and exited states are listed in Table 1.
In addition, NBO, atomic occupancies, Fock Matrix and IFCC [F (Δ), F (θ)] are listed in Table 2 and Table 3. The total energies of | α and | β spin orbitals for the ground state of strata/stratum, formed with symmetry C2V/CV (global minima of X   ˜ 2 Σ+), are E ( | α ) = 34.87083 and E ( | β ) = 34.15046 Hartree respectively, while these energies for CV/Dh (local minima) are E ( | α ) = 34.87079 and E ( | β ) = 34.15042 , respectively. Although the energy of A ˜ 2 Σ ( g ) + state is near the ground state, this excited configuration does not contribute to the ground state wave function.
In the symmetric D h geometry, the unpaired electron of BNB is delocalized, while in the asymmetric C v geometry it is localized on either one of the B atoms. Broken symmetry C v structures will be stabilized by this interaction relative to the symmetric D h geometry. Physically, the second-order Jahn-Teller interaction permits the unpaired electron to localize on a single boron atom, rather than being delocalized. The other two, which correspond to localizing the unpaired electron on one or both of the boron atoms (when the bond lengths are unequal), do not transform as an irreducible representation of D h for BNB radical via any changing of charges on N and either one of the B atoms (Figure 1a,a',b).
Figure 1. Relative energies of B2N(−, 0, +) versus B-N-B bond distance in various level of methods (a, a') cation; (b, b') radical and (c) anion.
Figure 1. Relative energies of B2N(−, 0, +) versus B-N-B bond distance in various level of methods (a, a') cation; (b, b') radical and (c) anion.
Molecules 20 19769 g001
When the molecule has D h symmetry, the real wave-function must transform as an irreducible representation of the D h point group. However, when the two B–N bonds are asymmetrically stretched, 4 σ g 2 and 3 σ u 1 become near degenerate, while 6 σ and 7 σ MOs have the same symmetry. Because of this vicinity, the singly excited state of [ c o r e ] 6 σ , 7 σ has a rather strong interaction with single and triple excitations. It is known however, that approximate electronic structure methods could suffer from an artifactual symmetry-breaking effect which would thus be confused as a real Jahn-Teller distortion. DFT methods such as B3LYP incorrectly underestimate the second-order Jahn-Teller distortion which leads the B3LYP calculations to predict a symmetric structure with too much high frequency for the anti-symmetric stretch. However, at this point, the “UHF” solution for the ground state wave-function exists [23].
Table 1. Geometric and electronic structures of B2N(−,0,+) in ground and exited states.
Table 1. Geometric and electronic structures of B2N(−,0,+) in ground and exited states.
State E e ( H a r t r e e ) | α Configuration | β Configuration r e ( B 1 N ) A 1 ( 2 , 1 , 3 , 2 , 1 )
(*Ne)Total Energy of | α * Total Energy of | β * , Virtual ** r e ( N B 2 ) A 2 ( 2 , 1 , 3 , 1 , 2 )
( H o m o L u m o )   **
X   ˜ 2Σ+ 104.0781959   a [ A ] , π 4 , σ 2 , σ 1 | α = 0.44641   a [ B ] σ 1 | β = 0.26180   a 1.3189   a A 1 = 179.9499   a
( C v ) 104.0820328   a * E ( | α ) ( 34.87083 )   a * ( 34.15046 )   a 1.3185   a A 2 = 179.9602   a
(*17e) 104.0754917   a * * ( 0.48614 )   a * * ( 0.17697 )   a
X   ˜ 2 Σ u + 104.0781959   b [ A ] π u 4 , 4 σ g 2 , 3 σ u 1 | α = 0.44641   b [ B ] , 4 σ g 1 | β = 0.26180   b 1.3187   b A 1 = 180.0   b
( D h ) * ( 34.87079 )   b * ( 34.15042 )   b 1.3187   b A 2 = 180.0   b
* * ( 0.48614 )   b * * ( 0.17697 )   b
X   ˜ 2 Σ u + 103.6396773   d [ A ] π u 4 , 4 σ g 2 , 3 σ u 1 | α =   0.42477   d [ B ] , 4 σ g 1 | β =-0.24844 d 1.317 6 d A 1 = 180.0   d
( D h ) * ( 34.74888 )   d * ( 34.06924 )   d 1.3176   d A 2 = 180.0   d
(*17e) * * ( 0.49245 )   d **   ( 0.18742 ) d
X   ˜ 2 Σ u + 104.159145   f [ A ] π u 4 , 4 σ g 2 , 3 σ u 1 | α =   0.44701   f [ B ] , 4 σ g 1 | β = 0.26341 1.3275   f A 1 = 180.0   f
( D h ) 104.1355512   n * ( 34.86501 ) * ( 34.14329 ) 1.3275 A 2 = 180.0
* * ( 0.48249 ) * * ( 0.17779 )
A   ˜ 2 Σ g + 104.1047582   k [ A ] π u 4 , 4 σ g 1 , 3 σ u 2 | α = 0.44638   u [ A ] π u 4 ,   4 σ u 1 , = 0.26134 1.3154   k A 1 = 180.0   k
(*17e) 104.0781729   K * ( 34.8744 )   u * ( 34.15408 ) 1.3154   k A 2 = 180.0   k
* * ( 0.48626 )   u * * ( 0.17666 )
  B ˜ 4 g 104.0297022   h [ A ] π u 4 , π g , 2 σ g 1 | α = 0.26899   h [ A ] , π u 2 | β 0.49646 1.3079   h A 1 = 180.0   h
(*17e) 104.0141196   a * ( 35.04484 )   h * ( 33.74801 ) 1.3079   h A 2 = 180.0   h
* ( 0.30125 )   h * * ( 0.50595 )
X   ˜ 1 Σ g + 104.1965676   a [ A ] 1 π u 4 , 4 σ g 2 , 3 σ u 2 = 0.13904   a 1.3291   a A 1 = 179.9322   a
(*18e) 104.2019114   a * E ( | α ) = 32.48647   a 1.3291 A 2 = 179.9462
104.1950169   a * * 0.36286
X   ˜ 1 Σ g + 104.1965676   b * E ( | α ) = 32.48644   b 1.3291   b A 1 = 180.0   b
(*18e) 1.3291 A 2 = 180.0
X   ˜ 1 Σ g + 104.1145486   c [ A ] π u 4 , σ g 2 , σ u 2 = 0.13454   c 1.3459   c A 1 = 179.8967   c
(*18e) 104.1162881   c 1.3459   c A 2 = 179.9181   c
104.112568   c
  A ˜ 3 Π u 104.0884685   a [C] π u 4 , 3 σ u 1 , 1 π g 1 | α = 0.03986   h [ C ] π u , 4 , 4 σ g 1 | β = 0.01559   h ) 1.3422   h A 1 = 180.0   h
(*18e) 104.1164696   h * E ( | α ) = 32.79721   h * E ( | β ) = 32.26335   h 1.3422 A 2 = 180.0
104.0809447   a * * 0.25941   h * * 0.14777   h
104.0792592   a
X   ˜ 1 Σ g + 103.7454365   a [ D ] ( 1 π u 4 ) , ( 4 σ g 2 = 0.57788   a ) 1.2938   a A 1 = 180.7911   a
(*16e) 103.8054934   a * E ( | α ) = 36.52419   a 1.2938   a A 2 = 180.6276   a
103.7545006   a * * 0.17926 a
X   ˜ 1 Σ g + 103.6023122   m [ D ] ( 1 π u 4 ) , ( 4 σ g 2 = 0.57803   g ) 1.3156   m   A 1 = 179.981 m
(*16e) 103.3012055   g * E ( | α ) = 36.51985   g 1.3156   m A 2 = 179.985 m
103.8377401   f * * 0.17944   g 1.3003   h
103.7903312   h 1.3004   h
B   ˜ 3Σg 103.7609202   a { [ A ] 4 σ g 1 , 3 σ u 1 , π u 4 , }   a { [ A ]   , π u 2 | β = 0.75222 }   a 1.2976   a A 1 = 180.0   a
(*16e) 103.7767141   h π u 2 | α = 0.74323   a * E ( | β ) ( 35.72475 ) a   1.2976   a A 2 = 180.0
103.7628505   a * E ( | α ) ( 37.33372 )   a * * ( 0.52131 )   a
103.7590557   a * * ( 0.53385 )   a
(a) QCISD/EPR-II, (a'') MP4D/EPR-III//QCISD/EPR-III, (m) CASSCF (a'') MP4SDQ/EPR-III//QCISD/EPR-III, (d) CASSCF (11, 12)/UHF, (g) CASSCF(10,12)rohf AUG-cc-pvqz, (b) QCISD/EPR-III( θ c o n s t = 180.0 ) , [   A   ] : 1 σ 2 , 2 σ 2 , 3 σ 2 , 4 σ 2 , 5 σ 2 , (c) QCISD/EPR-II, [ A ] : 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 ; (c') MP4D/EPR-II//QCISD/EPR-II, [ B   ] :   σ 1 , σ 1 , σ 1 , σ 1 , σ 1 , π 2 ; (c'') MP4SDQ/EPR-II//QCISD/EPR-II, [ B ]: 1 σ g 1 , 2 σ g 1 , 1 σ u 1 , 3 σ g 1 , 2 σ u 1 , 1 π u 2 , (f) CBS-lq, (h) QCISD(T)/EPR-III, [C]: 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 4 σ g 2 , [ C ] : 1 σ g , 1 σ u , 2 σ g , 3 σ g , 2 σ u , (n) CBS4O, (u) TD/EPR-II (k) TD/EPR-III//QCISD (T)/EPR-III, (k') TD/EPR-III//QCISD /EPR-III, [D]: 1 σ g 2 , 2 σ g , 2 1 σ u 2 , 3 σ g 2 , 2 σ u 2 . Total energy of | β and the energy of virtual orbital are shown with “*” and “**” symbols, respectively.
Table 2. NBO, electric potential, gradient of electric potential ( ) , atomic occupancies, and Fock Matrix data of B2N(−, 0, +) in ground and exited states.
Table 2. NBO, electric potential, gradient of electric potential ( ) , atomic occupancies, and Fock Matrix data of B2N(−, 0, +) in ground and exited states.
State B E p 1 N E p 2 B E p 3 = E P ( B ) E P ( N ) R ( B N ) Hybrids Coefficient& E a c c e p t o r ( j ) E D o n o r ( i ) Atomic Occupancies
(*Ne) * F o c k   M a t r i x   ( F i , j , a . u . )
X   ˜ 2 Σ u + N E p 2 = 18.3638   a ( 5.32 ) ,   ( 5.32 ) | ψ B D ( 1 ) = 0.91 S P N 1 1.01 + 0.41   S P B 2 2.46 a | ψ B D * ( 2 ) | ψ B D ( 1 ) = 2.01 | α N : 2 s 0.77 2 P x 0.84 2 P y 0.84 2 P z 0.87
(*17e) B E p 1 = B E p 3 11.3464   a | ψ B D ( 2 ) = 0.96 S P N 1 1.0 + 0.29 S P B 2 1.0 a * 0.094 | α B : 2 s 0.76 2 P x 0.79 2 P y 0.79 2 P z 0.40
| ψ B D ( 3 ) = 0.96 S P N 1 1.0 + 0.29 S P B 2 1.0 a | ψ B D * ( 3 ) | ψ B D ( 1 ) = 1.00 | β N : 2 s 0.80 2 P x 0.86 2 P y 0.86 2 P z 0.83
| ψ B D * ( 1 ) = 0.28 S P N 1 1.0 0.95   S P B 2 1.0 a * 0.089 | β B : 2 s 0.50 2 P x 0.69 2 P y 0.69 2 P z 0.16
| ψ B D * ( 2 ) = 0.42 S P N 1 1.01 0.90   S P B 3 0.95 a | ψ B D * ( 1 ) | ψ B D ( 2 ) = 0.7
| ψ B D * ( 3 ) = 0.70 S P B 2 0.97 0.70   S P B 3 0.97 a * 0.029
X   ˜ 1 Σ g + N E p 2 = 18.6013   a ( 5.25 ) ,   ( 5.25 ) | ψ B D ( 1 ) = 0.90 S P N 1 1.01 + 0.43   S P B 2 2.26 a | ψ B D * ( 3 ) | ψ B D ( 1 ) = 1.94 | α , β N : 2 s 1.59 2 P x 1.70 2 P y 1.70 2 P z 1.66
(*18e) B E p 1 = B E p 3 11.6190   a | ψ B D ( 2 ) = 0.66 S P N 1 1.0 + 0.75 S P B 2 1.0 a * 0.107 | α , β B : 2 s 0.97 2 P x 0.14 2 P y 0.14 2 P z 0.37
| ψ B D ( 3 ) = 0.95 S P N 1 1.0 + 0.32 S P B 2 1.0 a | ψ B D * ( 1 ) | ψ B D ( 1 ) = 0.73
| ψ B D * ( 1 ) = 0.44 S P N 1 1.01 0.90   S P B 2 2.24 a *0.021
| ψ B D * ( 2 ) = 0.75 S P N 1 1.0 0.67   S P B 3 1.0 a | ψ B D * ( 3 ) | ψ B D ( 3 ) = 1.94
| ψ B D * ( 3 ) = 0.32 S P N 1 1.0 0.94   S P B 3 1.0 a * 0.107
B   ˜ 3Σ+ N E p 2 = 18.10657   h ( 5.43 ) ,   ( 5.43 ) | ψ B D ( 1 ) = 0.90 S P N 1 1.0 + 0.42   S P B 0.72   h | ψ B D * ( 3 ) | ψ B D ( 1 ) = 1.98 | α N : 2 s 0.78 2 P x 0.81 2 P y 0.81 2 P z 0.87
(*16e) B E p 1 = B E p 3 11.0485   h | ψ B D ( 2 ) = 0.96 S P N 1 1.0 + 0.26 S P B 1.0 h * 0.119 | α B : 2 s 0.72 2 P x 0.92 2 P y 0.92 2 P z 0.44
| ψ B D ( 3 ) = 0.96 S P N 1 1.0 + 0.26 S P B 1.0 h | ψ B D * ( 1 ) | ψ B D ( 2 ) = 0.77 | β N : 2 s 0.77 2 P x 0.87 2 P y 0.87 2 P z 0.86
| ψ B D * ( 1 ) = 0.42 S P N 1 1.0 0.90   S P B 2 1.0 h * 0.028 | β B : 2 s 0.10 2 P x 0.62 2 P y 0.62 2 P z 0.08
| ψ B D * ( 2 ) = 0.26 S P N 1 1.01 0.96   S P B 3 0.95 h | ψ B D * ( 1 ) | ψ B D ( 2 ) = 0.7
| ψ B D * ( 3 ) = 0.42 S P N 1 1.0 0.90   S P B 3 0.72 h * 0.029
  A ˜ 3 Π u N E p 2 = 18.5665   h ( 5.19 ) ,   ( 5.19 ) | ψ B D ( 1 ) = 0.91 S P N 1 1.0 + 0.40   S P B 2 1.62 h | ψ B D * ( 3 ) | ψ B D ( 1 ) = 1.94 | α N : 2 s 0.75 2 P x 0.81 2 P y 0.60 2 P z 0.86
(*18e) B E p 1 = B E p 3 11.5955   h | ψ B D ( 2 ) = 0.97 S P N 1 1.0 + 0.26 S P B 2 1.0 h * 0.107 | α B : 2 s 0.74 2 P x 0.09 2 P y 0.70 2 P z 0.43
| ψ B D ( 3 ) = 0.91 S P N 1 1.0 + 0.40 S P B 3 1.62 h | ψ B D * ( 4 ) | ψ B D ( 1 ) = 1.01 | β N : 2 s 0.80 2 P x 0.87 2 P y 0.90 2 P z 0.86
| ψ B D ( 4 ) = 0.70 S P B 2 0.51 + 0.70 S P B 3 0.51 * 0.064 | β B : 2 s 0.46 2 P x 0.06 2 P y 0.04 2 P z 0.16
| ψ B D * ( 1 ) = 0.4 S P N 1 1.0 0.91   S P B 2 1.62 h | ψ B D * ( 1 ) | ψ B D ( 2 ) = 0.73
| ψ B D * ( 2 ) = 0.26 S P N 1 1.0 0.97   S P B 3 1.0 h * 0.021
| ψ B D * ( 3 ) = 0.40 S P N 1 1.0 0.91   S P B 3 1.62 h
| ψ B D * ( 4 ) = 0.70 S P B 2 0.51 0.70   S P B 3 051 h
The Fock Matrix ( F i , j , a . u . ) data have been shown with star symbol “*”. (a) QCISD/EPR-III, (h) QCISD(T)/EPR-III.
Table 3. Isotropic Fermi contact coupling (IFCC) of B2 N ( 0 ) in ground ( X   ˜ 2 Σ u + ) and exited (   B ˜ 4 g ) states.
Table 3. Isotropic Fermi contact coupling (IFCC) of B2 N ( 0 ) in ground ( X   ˜ 2 Σ u + ) and exited (   B ˜ 4 g ) states.
State Δ = r 2 r 1 I F C C [ f ( Δ ) ] Δ θ = θ 1 θ 1 I F C C [ f ( θ ) ]
(*Ne) θ c o n s = 180.0 N , B 1 ,   B 2 r c o n s = 1.3176 N , B 1 ,   B 2
X   ˜ 2 Σ u + (*17e) Δ   = 0.000 −29.81, 428.6, 428.6 Δ θ = 0.0 −29.81, 428.6, 428.6
Δ   = 0.010 −29.76, 386.2, 469.7 Δ θ = 2.0 −29.81, 428.71, 428.71
Δ   = 0.020 −29.11, 348.2, 508.1 Δ θ = 3.0 −29.81, 428.77, 428.77
Δ   = 0.030 −28.02, 316.1, 541.9 Δ θ = 5.0 −29.81, 428.95, 428.95
Δ   = 0.040 −26.64, 290.3, 570.2 Δ θ = 10.0 −29.83, 429.84, 429.84
Δ   = 0.050 −25.12, 270.7, 592.6 Δ θ = 20.0 −29.87, 433.44, 433.44
Δ   = 0.060 −23.54, 256.7, 609.5 Δ θ = 30.0 −29.86, 439.68, 439.68
Δ   = 0.070 −21.98, 247.6, 621.3 Δ θ = 40.0 −29.61, 448.75, 448.75
Δ   = 0.080 −20.46, 242.7, 628.7 Δ θ = 50.0 −28.85, 460.82, 460.82
Δ   = 0.082 −20.10, 242.1, 629.8 Δ θ = 60.0 −27.16, 475.86, 475.86
Δ   = 0.085 −19.72, 241.6, 630.9 Δ θ = 70.0 −24.06, 493.31, 493.31
Δ   = 0.086 4.35, 60.3, 899.9 Δ θ = 80.0 −18.57, 512.55, 512.55
Δ   = 0.087 4.46, 59.4, 901.2 Δ θ = 90.0 −10.65, 533.75, 533.75
Δ   = 0.090 4.79, 56.6, 904.7
  B ˜ 4 g Δ   = 0.000 20.26, 377.8, 377.8
Δ   = 0.010 20.26, 381.2, 377.9
Δ   = 0.020 20.26, 384.5, 377.9
Δ   = 0.030 20.28, 387.8, 377.9
Δ   = 0.040 20.30, 391.2, 377.9
Δ   = 0.050 20.30, 394.4, 377.9
Δ   = 0.060 20.38, 397.8, 377.9
Δ   = 0.070 20.44, 401.0, 377.8
Δ   = 0.080 20.50, 404.27, 377.7
It is prudent to employ a highly correlated method which can use a large number of reference determinants to recover dynamic and static correlations. In this work, the results of EPR-III basis sets are monotonous through the comparison between different situations. Although (in Table 1) the difference between the two positions of global minima and local minima for both | α and | β is 8.77 cm−1, our calculations show that the total energies for both of them are the same (i.e., 104.0781959 ). This is due to the fact that the spin orbital energies are related to the small bending angles of A1 and A2 which have an extremely low bending frequency (70 cm−1). Harmonic frequencies were determined at the QCISD/EPR-III//prop=EPR and characterized by 228.79 cm−1 ( ϑ 1 = ϑ 1 ,   bending mode π u ), 1178.64 cm−1 ( ϑ 2 , symmetric stretching “ σ g ”) and 2146.42 ( ϑ 3 , asymmetric stretching “ σ u ”). The IR and Raman intensities for ϑ 3 are 10165.0 and 0.00, respectively, while the ϑ 2 mode has intensity in Raman (51.0) but zero intensity in the IR region.
As it is shown in Table 1, the energy difference between two states of ( A   ˜ 2 Σ g + ) and ( X   ˜ 2 Σ u + ) {(k − a) and (k –  k )} are 5829.75 cm−1 and 5834.79 cm−1, respectively. Those values are close to the photoelectron spectroscopy calculation results which Asmis et al. have shown [2]. They have discussed that the signal observed in the 355 nm and 266 nm photoelectron spectra of   B 2   N has been indicated as due to a photodetachment from the X ˜ 1   Σ g +   ) to the ground and lowest excited state of neutral B2N { X ˜ 2   Σ u +   ) and (   A ˜ 2   Σ g +   ) with a linear symmetry and assigned to the X ˜ 1   Σ g +   X ˜ 2   Σ u + +   e and X ˜ 1   Σ g +   A ˜   2   Σ g + +   e transitions {the (   A ˜ 2   Σ g +   ) term energy T0 is 0.785 eV or 6331.77 cm−1}.
The 1 π u orbital is a bonding combination of all 2 p   π orbitals on all three atoms, while 4 σ g 2 and 3 σ u 1 orbitals are close lying and not strongly bonding in character. The small separation of these two orbitals accounts for the small energy required to promote the molecule from the ground state to the first exited state at 5829.75 cm−1. The difference between (k − a) and (k − k ) is 5 cm−1 which is near 8.77 cm−1 (different between C 2 V / C V and C V / D h of strata/stratum). In addition   B ˜ 2 g (exited state above the A   ˜ 2 Σ g + ) is subject to the Renner-Teller effect, leading to a complicated pattern of bending vibrational levels. Our calculation shows the existence of a larger gap between 3 σ u and 1 π g orbitals, thereby placing the transition to [ A ] π u 4 , 4 σ g 2 , 1 π g 1 with   B ˜ 2 g state much higher in energy. Analysis of the vibronic structure of the   B ˜ 2 g ( X   ˜ 2 Σ u + ) band system shows the transition to B ˜ 2 g at 19,452 cm−1. Nevertheless, the nonbonding character of 1 π g and 3 σ u orbitals implies no significant change in B–N bond lengths in this transition, as it is observed by Ding et al. [56]. Therefore, in spite of A ˜ 2 Σ ( g ) + , the “   B ˜ 2 g ” excited configuration does contribute to the ground state wave function as a subject of Renner-Teller effect.
It is notable that the electron configuration sequences for the | α and | β situations are not the same and as shown in Table 1. These sequences correspond to [ A ] and [ B ] for | α and | β , respectively, which result the sequences of 1 σ g > 1 σ u > 2 σ g for | α and 1 σ g > 2 σ g > 1 σ u for | β . They indicate the relation between Centro-symmetric and electron correlations in various levels from UHF and DFT to the CASSCF and MRCI methods for the magnitude of the barrier energies. As a result, if the molecule has a barrier at the Centro-symmetric structure, the barrier is not extraordinarily high, and if a high barrier were present, the vibrational level of the ground state would be thermally populated and transitions from this level would be sharp. Ding shows that in order to decide whether a barrier exists at the Centro-symmetric configuration, the energies of the ( ϑ 3 ) levels in the ground state would have to be measured, and the resulting set of levels used to deduce the potential function along the anti-symmetric stretching coordinate, Q3. For even values of v3, this would be possible in principle by either dispersed fluorescence or stimulated emission pumping spectroscopy. In practice, however, it appears that these levels will be difficult to reach owing to poor Franck-Condon factors. For odd values of ϑ 3 , a direct infrared absorption study provides the best method for the accurate measurement of those levels [56].
Although the observation of excitations involving uneven quanta of the anti-symmetric stretching mode, v3, indicates a breakdown of the Franck–Condon (FC) approximation, it cannot be the only results from Herzberg–Teller vibronic coupling between the ( X ˜ 2   Σ u + ) ) and ( ( A ˜ 2   Σ g + ) ) states involving the v3 mode.
In Table 2, the NBO calculation shows that the 2sN orbital is considered to be primarily core-like, forming the 3 σ g orbital, though, of course, some mixing of the 2sN orbital into the other σ g orbitals ix expected. NBO analysis of the orbital containing the unpaired electron in BNB shows that most of the spin density is located in the boron sp orbitals. The boron atomic orbitals are best described as a “sp” hybrid, directed away from the nitrogen atom however, bonding with respect the σ P orbital on the nitrogen atom. The 2 σ u orbital is a bonding combination of the 2 p σ orbital on the central nitrogen atom with 2 s p σ hybrid orbitals on the two borons. The 1 π u orbital is a bonding combination of all 2 p π orbitals on all three atoms, while 4 σ g and 3 σ u orbitals are close lying and not strongly bonding in character.
The hyperfine parameters were calculated for the linear geometry with a bond length of 1.3176 Å via a CASSCF optimization at several levels of configuration interaction and exited states (Table 3).
A i s o in the ground state for boron atoms varies from 428.6 MHz at Δ   = 0.0 to 241.6 MHz at Δ   = 0.085 for B1 and from 428.6 MHz at Δ   = 0.0 to 630.9 MHz at Δ   = 0.085 for B2 while the nitrogen A i s o varies smoothly from −29.81 MHz ( Δ   = 0.0 ) to −19.72.0 MHz ( Δ   = 0.085 ). There is a critical point for A i s o (both in B and N) between Δ   = 0.085 and Δ   = 0.086 where the data reverses to 4.35, 60.3and 899.9 for N, B1 and B2 respectively (Table 3), indicating the symmetry breaking in point charges or distances in this region. The dipolar hyperfine coupling constants exhibited negligible dependence on bound. The highest occupied orbital containing the unpaired electron is   σ u orbital with most of the electron density on the boron atoms. Symmetry constrains this orbital to have only “ σ P ” orbital contributions from nitrogen with no “s” character so that the isotropic hyperfine parameter from nitrogen is small and arises mostly from spin polarization effects.
The molecular isotropic hyperfine ( A i s o ) values for BNB can be obtained from the experimental “ A ” and “ A ”: A i s o = ( 2 A + A ) 3 = 8 π 3   g e g N β e β N < ( r ) > (10)
As a result, the experimental data are B1 = B2 = 451 and N = −14.0 [4] which are close to our calculated results. Isotropic hyperfine interaction including “g” tensor for BNB shows the large boron isotropic amount. Thus, the properties of the g tensor eliminate the possibility of a low exited 2Πg state for the radical. The averaged vibrations of B–N bond lengths, I F C C [ f ( Δ ) ] and I F C C [ f ( θ ) ] can be shown, even if the ground state is precisely linear, quasi-linear or the geometry at the potential minimum went through a symmetry breaking to form a C v structure.
In this study, we have focused on the charge distribution of boron and nitrogen atoms to exhibit the charge breaking (Table 1). As it is shown in Figure 1 and Figure 2, the symmetry breaking in point charge distribution does not only indicate a barrier energy in the regions of 100–160 cm−1 [9,22,23] or 20 cm−1 [20], but also it creates several SB through the asymmetry stretching (or interaction between asymmetry stretching and bending) with different barrier energies from high to small values (about 5 cm−1) (Figure 1).
Figure 2. Relative energies of B2N(−, 0, +) vs. B-N-B bond angle in various level of methods for anion, cation and radical forms of BNB respectively
Figure 2. Relative energies of B2N(−, 0, +) vs. B-N-B bond angle in various level of methods for anion, cation and radical forms of BNB respectively
Molecules 20 19769 g002
Therefore, the symmetry-breaking barrier has a dynamic changing with no Centro-symmetric form, and it depends on wave function or charge distribution. Furthermore, a large barrier can be estimated via fixed-node diffusion Monte Carlo methods in which D h and C v ROHF WFs (with symmetric stretching) have energies separated by a gap as big as 0.5 milli Hartree approximately [57]. In the case of C V symmetry solutions, there are two degenerate solutions which correspond to a single electron, localized on either one of the two “different” boron atoms, as can be observed from a natural bond analysis. A linear combination between these two solutions will restore the D h symmetry of the WF, but this WF would be different from the original D h solution where the single electron is delocalized over the two boron atoms [57].
As shown in Figure 3, based on Knight’s reports [4] {(in which the variation of energy with bond angle for finding a minimum around (100°) have been discussed}, a cyclic radical or anion B2N(−, 0) in our calculations has not been observed, though, for cation, there is a bulge in the curve at 90° in MP4DQ and MP4DSQ methods which indicates a cyclic B2N+.
Figure 3. Isotropic Fermi contact coupling (IFCC) of B2N(0): (a) Function of angles changing; (b) Function of distances changing. The axes for nitrogen and borons have different scales; (a') IFCC for Nin a short scale between −30 to −29.8 indicates of IFCC breaking and (a'') different IFCC of B and N in a short scale.
Figure 3. Isotropic Fermi contact coupling (IFCC) of B2N(0): (a) Function of angles changing; (b) Function of distances changing. The axes for nitrogen and borons have different scales; (a') IFCC for Nin a short scale between −30 to −29.8 indicates of IFCC breaking and (a'') different IFCC of B and N in a short scale.
Molecules 20 19769 g003
At the SCF level, the lowest energy corresponds to a bent molecule with an angle of 100°, however, for the QCISD (T), MP4DQ, MP4DSQ and HF/aug-cc-pVTZ calculations (Figure 2) the linear structure clearly has the lowest energies for radical and anion structures. Martin [7] has shown a cyclic B2N (2B2) via reactions of pulsed laser produced boron and nitrogen atoms in a condensed argon stream (at higher laser power reactions) and has discussed that the vibration 882.3 cm−1 must be considerably an-harmonic. This possibility receives substantial support from the five combination bands observed in the 3000–6500 cm−1 regions. The 882.3 cm−1 one is assigned to the anti-symmetric B-N stretching fundamental v3(b2) of cyclic B2N and the 1998.4 cm−1 combination band is the sum of v1 (a1), the symmetric B-N stretching fundamental, and v3. The difference 1998 − 882 = 1116 cm−1 can help to measure the v1.
The electronic transition energy σg → σu excitation is less than 6000 cm−1, which indicates that the higher overtones of the cyclic B2N (2B2 state) vibrations will display significant vibronic interaction effects. The failure to observe cyclic B2N in the 2B2 state by ESR [4] is most likely due to the differences in production and relaxation of the energized evaporated species. The fact that the radical BNB might be converted to cation BNB towards the cyclic B2N (due to the laser ionization effect) can be predicted from Martin’s study. In addition, Becker et al. [58] used laser ionization mass spectrometry to study the formation of B n N m + clusters ions in laser plasma which resulted in production of BNB cation whereas our calculations resulted in production of cyclic B2N+.
In a simple form, the best Lewis structure representation (Scheme 2) is the pair of resonance structures, i.e., “.B=N=B:” and “: B=N=B.” These structures have a formal +1 charge on nitrogen and a formal −1/2 charge on each boron, but the difference in electronegativity gives rise to a net-negative charge density on the nitrogen atom. Radical of B2 N ( 0 ) is stabilized by six resonance structures of the linear forms and three resonance structures of cyclic forms (with various distributions) as follows:
Scheme 2. Lewis structures representation of BNB in three forms of radical, anion and cation.
Scheme 2. Lewis structures representation of BNB in three forms of radical, anion and cation.
Molecules 20 19769 g006
Although in some structures of the six resonances, the nitrogen is negative and the two borons are positive, in our ESP calculations for B2 N ( 0 ) with different levels of theory, nitrogen is always positive and close to zero while the two boron atoms are always negative near zero and the converse appears for the B2N(+) species. The atomic charges via 4516 points for one of the ESP fitting at QCISD/EPR-III level of theory for radical shows the values of 0.076055, −0.037580 and −0.038476 for N and two boron atoms, respectively.
All in all, the nuclear hyperfine of various states are shown in Table 3 and Figure 3. The nitrogen and one of the boron atoms in   B ˜ 4 g state are independent of Δ during hyperfine coupling calculation. Although these two states   B ˜ 4 g and ground state) can interact along the anti-symmetric stretching as a reason of SB problem in BNB radical, the lack of the “correct” per-mutational symmetry of the wave-functions which arises due to the oversimplification of the wave function is a major reason for spontaneous symmetry breaking.
For the sextuplets contaminating C V symmetry 2Σ+:   1 σ 2 , 2 σ 2 , 3 σ 2 , 4 σ 2 , 5 σ 2 , 6 σ 2 , 7 σ 1 , 1 π 2 , 1 π 2 and 6Σ+: : 1 σ 2 , 2 σ 2 , 3 σ 2 , 4 σ 2 , 5 σ 2 , 6 σ 2 , 1 π 2 , 7 σ 1 , 1 π 2 the dynamic correlation would be even smaller [1]. With a cc-pVQZ basis set and high correlation, the single reference CCSD (T) energy gap between the symmetric and asymmetric configurations is 136 cm−1 [23] which is reduced to 99 cm−1 in the RMR CCSD (T) method [14]. An additional extended discussion would be reported in a subsequent publication concerning strata/stratum (S/s) configuration.
We have reinvestigated the anion form at the QCISD (T), QCISD, CCSD (T), MP4SDQ, MP4D and full space CASSCF levels of theory employing Aug-CC-PVTZ and EPR-III basis sets. Although by our calculations the cyclic  B 2   N anion structure cannot be predicted to lie above the X   ˜ 1   Σ g + ground state, it is notable that the lowest stable bent solution for B2 N ( ) should be a 3B2 state and not a single state.
Harmonic frequencies were determined at the QCISD/EPR-III//prop = EPR and characterized by 224.70 cm−1 ( ϑ 1 = ϑ 1 ,   bending mode “ π u ”) and 1203.42 cm−1 ( ϑ 2 , symmetric stretching “ σ g ”) and 1837.50 ( ϑ 3 , asymmetric stretching “ σ u ”). The IR and Raman intensities for ϑ 3 are 1494.6 and 0.00, respectively, while the ϑ 2 ” mode has intensity in Raman (39.97), but zero intensity in IR region. The valence orbital occupancy of ground state ( X   ˜ 1   Σ g + ) is: 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 1 π u 4 , 4 σ g 2 , 3 σ u 2 , while the lowest excited triplet state in the D h representation for anion form has calculated (Table 1) and is predicted to be a   A ˜ 3 Π u state, lying 2.94 eV {   A ˜ 3 Π u ( a ) X ˜ 1 Σ g + ( a ) in Table 1} above the X   ˜ 1 Σ g + state (with the calculated value of 2.7 eV anion with photoelectron spectroscopy of B2N(−)).
The energies of | α and | β valance occupied MOs for the first triplet of exited state ( 1 σ g 2 , 1 σ u 2 , 2 σ g 2 , 3 σ g 2 , 2 σ u 2 , 4 σ g 2 , 1 π u 4 , 4 σ u 1 , 1 π g 1 ) are 1 π g 1 | α = 0.03986 , 4 σ u 1 | α = −0.23701 and { 1 π u 2 | α = 0.25466 ,   0.25136 }   meanwhile, for | β are { 1 π u 2 | β = 0.28802 , 0.27730   } , and 4 σ g 2 | β = 0.01559 . The total energy for | α is more than the total energy for | β which indicates more correlation in anion form compared to the other two forms of B N ( + ,    0 ) .
However, the unrestricted QCISD (T) wave function is considerably spin-contaminated and characterized by a large T1 value. The A ˜ 3 Σ g + excited state corresponding to the promotion of an electron from the highest occupied molecular orbital (HOMO), 3 σ u , to the 4 σ u MO is predicted to be considerably above the   A ˜ 3 Π u state.
4516 points have been used for fitting the atomic charges to electrostatic potential and the charges of ESP in this fitting are: {N = 0.891292, B = −0.945425, B = −0.945867}. Nitrogen is always positive while the two boron atoms are always negative and it is appeared conversely for the B2 N ( + ) in which the nitrogen in cation is negative and two borons are positive in all ESP calculations. Both of the two highest occupied orbitals are predominantly nonbonding with the electron density localized mainly on the “terminal” boron atoms.
Isotropic Fermi contact coupling (IFCC) of B2N(0) as a function of angles and distances is shown in Table 4 and Figure 3. Although the symmetry breaking cannot be seen in the IFCC { f ( Δ r ) } , it can be seen in the IFCC { f ( Δ θ ) } for the nitrogen in range of angles between Δ θ = (50, −50) (Figure 3a'). On the other hand, the SB has not been observed for the Δ ( Δ θ ) = Δ θ B Δ θ N in range of Δ θ = (5, −5) (Figure 3a''). It seems that nitrogen in SB problem plays a major role and it depends on asymmetry bond changing and angle deformation interaction ( π u + σ u ) .
Table 4. {TQ#: Traceless Quadrupole moment (Debye-Ang)}; (§) Charges from ESP fitting and Isotropic Fermi Contact Couplings (MHz) (IFCC).
Table 4. {TQ#: Traceless Quadrupole moment (Debye-Ang)}; (§) Charges from ESP fitting and Isotropic Fermi Contact Couplings (MHz) (IFCC).
State E e ( H a r t r e e ) B δ q 1 N δ q 2 B δ q 3 ( § ) T Q x x #
(* Ne)IFCC(N,B1,B2) T Q y y
T Q z z
X   ˜ 2Σ+ 104.0781959   a δ q 1 = δ q 3 = 0.015   a 1.8265   a
( C v ) 104.0820328   a δ q 2 = 0.03   a 0.8644   a
(*17e) 104.0754917   a '' 29.8 ,   427.2 ,   429.6   a 0.9621   a
X   ˜ 2 Σ u + 104.0781959     b δ q 1 = δ q 3 = 0.015   b 0.9621   b
( D h ) δ q 2 = 0.03   b 0.9621   b
104.159145   f 29.9 ,   428.2 ,   428.2   b 1.9241   b
X   ˜ 2 Σ u + 103.6396773   d δ q 1 = δ q 3 = 0.0165   d Q x x = 16.5   d
( D h ) δ q 2 = 0.033 d Q y y = 16.5   d
(*17e) Q z z = 19.8   d
X   ˜ 2 Σ u + 104.159145   f δ q 1 = δ q 3 = 0.0045   f 0.8928   f
( D h ) 104.1355512   n δ q 2 = 0.009 0.8928
1.7855
A   ˜ 2 104.1047582   k 137.9 ,   817.0 ,   817.0   u 17.87   u
(*17e) 104.0781729   K 17.87   u
20.80   u
  B ˜ 4 g 104.0297022   h 20.26 ,   377.84 ,   377.84   h 1.0334
(*17e) 104.0141196   a -2.4251
1.3917
X   ˜ 1 104.1965676   a δ q 1 = δ q 3 = 0.934 14.8029   a
(*18e) 104.2019114   a δ q 2 = 0.868   a 7.4015   a
104.1950169   a '' −14.98, 211.91, 211.91 7.4015   a
X   ˜ 1 δ q 1 = δ q 3 = 0.934 7.4015   b
(*18e) 104.1965676   b δ q 2 = 0.868   b 7.4015   b
14.803   b
X   ˜ 1 104.1145486   c δ q 1 = δ q 3 = 0.987 15.7996   c
(*18e) 104.1162881 c δ q 2 = 0.974   c 7.8992   c
104.112568   c '' 7.9004   c
  A ˜ 3 Π u 104.0884685   a δ q 1 = δ q 3 = 0.65   h 2.444   h
(*18e) 104.1164696   h δ q 2 = 0.3 6.5747
104.0809447   a 9.0187
X   ˜ 1 103.7454365   a δ q 1 = δ q 3 = 0.805   a 8.0329   a
(*16e) 103.8054934   a δ q 2 = 0.611 4.0164
103.7545006   a 4.0164
X   ˜ 1 103.6023122   m δ q 1 = δ q 3 = 0.828   f 7.009   m
(*16e) 103.3012055   g δ q 2 = 0.656   f 3.505   m
103.8377401   f 3.505   m
B   ˜ 3Σg 103.7609202   a δ q 1 = δ q 3 = 0.763   a 3.6992   a
(*16e) 103.7767141   h δ q 2 = 0.525   a   3.6992
103.7628505   a 7.3983
(a) QCISD/EPR-III; (d) CASSCF(11,12)/UHF; (a') MP4D/EPR-III//QCISD/EPR-III; (m) CASSCF(10,12)/EPR-II; (a'') MP4SDQ/EPR-III//QCISD/EPR-III; (g) CASSCF(10,12)rohfAUG-cc-pvqz; (b) QCISD/EPR-III; (c) QCISD/EPR-II; (c') MP4D/EPR-II//QCISD/EPR-II; (c'') MP4SDQ/EPR-II//QCISD/EPR-II; (f) CBS-lq; (h) QCISD(T)/EPR-III; (n) CBS4O; (u) TD/EPR-II; (k) TD/EPR-III//QCISD(T)/EPR-III; (K') TD/EPR-III//QCISD/EPR-III.
The role of nitrogen can also be discussed regarding the quadrupole moment. In Table 1, the traceless quadrupole moments in several ground and exited states are listed and it can be seen that the tensors of T Q x x and T Q y y are equal and positive, while T Q z z is negative when the symmetry of radical is D h (in ground state). However, all three components are negative in the exited state for (   A ˜ 2 Σ g +   ) and their values are in a large scale in comparison to those in ground state while for   B ˜ 4 g they are irregular in the range of ground states. So, B ˜ 4 g X   ˜ 2 Σ+ interaction should be stronger than A ˜ 2 Σ g + X   ˜ 2 Σ + . These results are the same for the anion in ground state of   A ˜ 3 Π u and the excited state of   A ˜ 3 Π u and the suitable interaction of   A ˜ 3 Π u   A ˜ 3 Π u can be predicted (it is ˜ 1 Σ +   ˜ 3Σg for cation). In other words, the quadrupole moment for BNB is highly sensitive to angle deformation and bond distance changing, so it can be discussed for any SBs in terms of quadrupole moment.
Based on our previous work [59], we have modified the Columbic term of the Schrödinger equation for the definition of parameters α α  and  β β . The left part in the chart for the three curves of anion, cation and radical of these parameters are straight towards up and down (Figure 4) and their values are 1.293, 1.310 and 1.329 for cation, radical and anion, respectively. In the right part there is a point splitting between two curves of α and β (1.303, 1.316 and 1.339 for cation, radical and anion, respectively) indicating a proper region for charge stability distribution in the SB problem.
Figure 4. Parameters of Charge correction coefficients, vs. changing of boron and nitrogen distances for B2N(−, 0, +) for anion , cation and radical.
Figure 4. Parameters of Charge correction coefficients, vs. changing of boron and nitrogen distances for B2N(−, 0, +) for anion , cation and radical.
Molecules 20 19769 g004

5. Conclusions

In this study, we have shown that the SB problem is not a real phenomenon; it is a hidden function depending on various variables such as charge distribution, bond length, IFCC, primitive Gaussians, trial wave-function properties, frozen core electrons and most importantly, non-Born–Oppenheimer approximation approach. While the nuclear repulsion energy of B2N(−, 0, +) in Born–Oppenheimer approximation mostly depends on the variables such as B-N bond length, using large and larger basis sets and more and more electron correlation are doomed to result in wrong limit for the energy level of SB barriers or SB estimation. It is prudent to employ a highly correlated method which can use a large number of reference determinants to recover dynamic and static correlations. We have shown that the SB is generally applied to the failure of the electronic wave function in order to be transformed as an irreducible representation of the molecular point group, so the failure of the electronic wave function is purely artifactual. It is not wise to conclude that the only special level of theory on the symmetry breaking for BNB is real (which has been concluded in reference [9]) because there exist some hidden variables in the electronic wave functions which should be considered. We have found that the symmetry breaking (SB) for some hidden variables (such as charge distribution) not only exhibit an energy barrier, it also creates several SBs through the asymmetry stretching with bending mode interaction.

Author Contributions

Majid Monajjemi: The Idea and designed of the whole research, S.B.: wrote and typed the manuscript, N.M.: Some of the calculations, M.Z. & N.S. & N.A.: Designed the Tables and Figures, M.A., M.S.M. & H.A. analyzed some of the data to collaborate with Majid Monajjemi , M.S., M.S.K., G.N. checked the dictation and revised some section.

Conflicts of Interest

The authors declare no conflict of interest.

References and Notes

  1. Martin, J.M.L.; Francois, J.P.; Gijbels, R.J. Ab initio study of boron, nitrogen, and boron-nitrogen clusters. I. Isomers and thermochemistry of B3, B2N, BN2, and N3. J. Chem. Phys. 1989, 90, 6469–6985. [Google Scholar] [CrossRef]
  2. Asmis, K.R.; Taylor, T.R.; Neumark, D.M. Anion photoelectron spectroscopy of B2N. J. Chem. Phys. 1999, 111, 8838–8851. [Google Scholar] [CrossRef]
  3. Monajjemi, M.; Lee, V.S.; Khaleghian, M.; Honarparvar, B.; Mollaamin, F. Theoretical Description of Electromagnetic Nonbonded Interactions of Radical, Cationic, and Anionic NH2BHNBHNH2 Inside of the B18N18 Nanoring. J. Phys. Chem. C 2010, 114, 15315–15330. [Google Scholar] [CrossRef]
  4. Knight, L.B., Jr.; Hill, D.W.; Kirk, T.J.; Arrington, C.A. Matrix Isolation ESR and Theoretical Investigations of l1Bl4Nl1Band loBl4Nl1B: Laser Vaporization Generation. J. Phys. Chem. 1992, 96, 555–561. [Google Scholar] [CrossRef]
  5. Monajjemi, M.; Khaleghian, M. EPR Study of Electronic Structure of [CoF6]3 2and B18N18 Nano Ring Field Effects on Octahedral Complex. J. Clust. Sci. 2011, 22, 673–692. [Google Scholar] [CrossRef]
  6. Monajjemi, M. Non bonded interaction between BnNn (stator) and BNB (rotor) systems: A quantum rotation in IR region. Chem. Phys. 2013, 425, 29–45. [Google Scholar] [CrossRef]
  7. Andrews, L.; Hassanzadeh, P.; Burkholder, T.R.; Martin, J.M.L. Reactions of pulsed laser produced boron and nitrogen atoms in a condensing argon stream. J. Chem. Phys. 1993, 98, 922–931. [Google Scholar] [CrossRef]
  8. Monajjemi, M.; Boggs, J.E. A New Generation of BnNn Rings as a Supplement to Boron NitrideTubes and Cages. J. Phys. Chem. A 2013, 117, 1670–1684. [Google Scholar] [CrossRef] [PubMed]
  9. Li, X.; Paldus, J.J. Real or artifactual symmetry breaking in the BNB radical: A multireference coupled cluster viewpoint. Chem. Phys. 2007, 126. [Google Scholar] [CrossRef] [PubMed]
  10. Martin, J.M.L.; El-Yazal, J.; Francois, J.-P.; Gijbels, R. The structure and energetics of B3N2, B2N3, and BN4 Symmetry breaking effects in B3N2. Mol. Phys. 1995, 85, 527–537. [Google Scholar] [CrossRef]
  11. Jalilian, H.; Monajjemi, H. Capacitor simulation including of X-doped graphene (X = Li, Be, B) as two electrodes and (h-BN)m (m = 1–4) as the insulator. Jpn. J. Appl. Phys. 2015, 54. [Google Scholar] [CrossRef]
  12. Monajjemi, M.; Jafari Azan, M.; Mollaamin, F. Density Functional Theory Study on B30N20 Nanocage in Structural Properties and Thermochemical Outlook. Fuller. Nanotub. Carbon Nanostruct. 2013, 21, 503–515. [Google Scholar] [CrossRef]
  13. Paldus, J.; Li, X.A. Critical Assessment of Coupled Cluster Method in Quantum Chemistry. Adv. Chem. Phys. 1999, 110. [Google Scholar] [CrossRef]
  14. Paldus, J. Handbook of Molecular Physics and Quantum Chemistry, Part 3; Wilson, S., Ed.; Wiley: Chichester, UK, 2003; Volume 2, p. 272. [Google Scholar]
  15. Esmkhani, R.; Monajjemi, M. Electronic Structural Investigation of Boron Nitride Nano Cage (B30N20) in Point of Exchange and Correlation Energy. J. Comput. Theor. Nanosci. 2015, 12, 652–659. [Google Scholar] [CrossRef]
  16. Wang, Z.X.; Huang, M.B.; Schleyer, P.V.R. Theoretical Study of Boron−Ammonia Reactions. J. Phys. Chem. A 1999, 103, 6475–6484. [Google Scholar] [CrossRef]
  17. Monajjemi, M.; Hosseini, M.S.; Mollaamin, F. Theoretical Study of Boron Nitride Nanotubes with Armchair Forms. Fuller. Nanotub. Carbon Nanostruct. 2013, 21, 381–393. [Google Scholar] [CrossRef]
  18. Asmis, K.R.; Taylor, T.R.; Neumark, D.M. Anion photoelectron spectroscopy of small boron nitride clusters: Adiabatic detachment energies and vibrational frequencies of low-lying electronic states in B2N and B3N. Eur. Phys. J. D 1999, 9, 257–261. [Google Scholar] [CrossRef]
  19. Monajjemi, M.; Hosseini, M.S. Non Bonded Interaction of B16N16 Nano Ring with Copper Cations in Point of Crystal Fields. J. Comput. Theor. Nanosci. 2013, 10, 2473–2477. [Google Scholar] [CrossRef]
  20. Liu, Y.; Zou, W.; Bersuker, I.B.; Boggs, J.E. Symmetry breaking in the ground state of BNB: A high levelmultireference study. J. Chem. Phys. 2009, 130. [Google Scholar] [CrossRef] [PubMed]
  21. Kalemos, A. Symmetry breaking in a nutshell: The odyssey of a pseudo problem in molecular physics. The X ˜ 2 Σ u +   BNB case revisited. J. Chem. Phys. 2013, 138. [Google Scholar] [CrossRef] [PubMed]
  22. Stanton, J.F. An unusually large nonadiabatic error in the BNB molecule. J. Chem. Phys. 2010, 133. [Google Scholar] [CrossRef] [PubMed]
  23. Gwaltney, S.R.; Head-Gordon, M. Calculating the equilibrium structure of the BNB molecule: Real vs. artifactual symmetry breaking. Phys. Chem. Chem. Phys. 2001, 3, 4495–4500. [Google Scholar] [CrossRef]
  24. Walsh, A.D. The electronic orbitals, shapes, and spectra of polyatomic molecules. Part I. AH2 molecules. J. Chem. Soc. 1953, 2260–2266. [Google Scholar] [CrossRef]
  25. Miransky, V.A. Dynamical Symmetry Breaking in Quantum Field Theories; World Scientific Publishing Co Pte Ltd: Singapore, 1994; p. 15. [Google Scholar]
  26. The spontaneous symmetry breaking has fall out when a system has a symmetry group, but requires that one element of the group be distinct. It is not important to bring up which member is distinct, only that one is. So if this element really is distinct, with the condition that any results found in this way must be re-symmetrized by tacking the average of each of the elements of the group being the distinct one. In the SSB processing of the BNB (radical) the system in a symmetrical state (Dh), ends up to a lower symmetrical state (Cv). In SSB the equation of motions obeycertain symmetries, but the lowest energy solutions do not exhibit that symmetry.
  27. Gasser, J.; Leutwyler, H. Chiral perturbation theory to one loop. Ann. Phys. 1984, 158, 142–210. [Google Scholar] [CrossRef]
  28. Davidson, E.R.; Borden, W.T. Symmetry breaking in polyatomic molecules: real and artifactual. J. Phys. Chem. 1983, 87, 4783–4790. [Google Scholar] [CrossRef]
  29. Lowdin, P.O. Wave and Reaction Operators in the Quantum Theory of Many-Particle Systems. Rev. Mod. Phys. 1963, 35, 702–707. [Google Scholar] [CrossRef]
  30. Adler, S.L. Comments on Proposed Gravitational Modifications of Schrodinger Dynamics and their Experimental Implications. J. Phys. A 2007, 40, 755–764. [Google Scholar] [CrossRef]
  31. Wigner, E. Einige Folgerungen aus der Schrödingerschen Theorie für die Term. Gött. Nachr. 1930, 133–146. [Google Scholar]
  32. Diósi, L. A Gravitation and quantum-mechanical localization of macro-objects. Phys. Lett. 1984, 105, 199–202. [Google Scholar] [CrossRef]
  33. Penrose, R. On Gravity’s role in Quantum State Reduction. Gen. Relativ. Grav. 1996, 28, 581–600. [Google Scholar] [CrossRef]
  34. Penrose, R. Quantum computation, entanglement and state reduction. Phil. Trans. R. Soc. Lond. A 1998, 356, 1927–1939. [Google Scholar]
  35. Aharonov, Y.; Englert, B.G.; Scully, M.O. Protective measurements and Bohm trajectories. Phys. Lett. A 1999, 263, 137–146. [Google Scholar] [CrossRef]
  36. Bohm, D. A Suggested Interpretation of the Quantum Theory in Terms of “Hidden” Variables. I. Phys. Rev. 1952, 85, 166–193. [Google Scholar] [CrossRef]
  37. Bohm, D.; Hiley, B.J. 1993 the Undivided Universe: An Ontological Interpretation of Quantum Theory; Routledge: London, UK, 2003. [Google Scholar]
  38. Holland, P. The Quantum Theory of Motion: An Account of the de Broglie-Bohm Causal Interpretation of Quantum Mechanics; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  39. Damljanovic, V. Structure and dynamics of Xn-type clusters (n = 3, 4, 6) from spontaneous symmetry breaking theory. Phys. Scr. 2013. [Google Scholar] [CrossRef]
  40. Barone, V. Recent Advances in Density Functional Methods, Part I; Chong, D.P., Ed.; World Scientific Publ. Co.: Singapore, 1996. [Google Scholar]
  41. Chopra, N.G.; Luyken, R.J.; Cherrey, K.; Crespi, V.H.; Cohen, M.L.; Louie, S.G.; Zettl, A. Boron nitrid nanotubes. Science 1995, 269, 966–967. [Google Scholar] [CrossRef] [PubMed]
  42. Blasé, X.; Rubio, A.; Louie, S.G.; Cohen, M.L. Stability and band gap constancy of boron nitride nanotubes. Eur. Phys. Lett. 1994, 28, 335–340. [Google Scholar] [CrossRef]
  43. Blasé, X.; Charlier, J.C.; de Vita, A.; Car, R. Theory of composite BxCyNz nanotube heterojunction. Appl. Phys. Lett. 1997, 70, 197–199. [Google Scholar] [CrossRef]
  44. Chopra, N.G.; Benedict, L.X.; Crespi, V.H.; Cohen, M.L.; Louie, S.G.; Zettl, A. Fully Collapsed Carbon Nanotubes. Nature 1995, 377, 135–138. [Google Scholar] [CrossRef]
  45. Walker, T.E.H.; Richards, W.G. Molecular spin-orbit coupling constants. The role of core polarization. J. Chem. Phys. 1970, 52, 1311–1314. [Google Scholar]
  46. Koseki, S.; Schmidt, M.W.; Gordon, M.S. MCSCF/6-31G(d,p) calculations of one-electron spin-orbit coupling constants in diatomic molecules. J. Phys. Chem. 1992, 96, 10768–10772. [Google Scholar] [CrossRef]
  47. Pople, J.A.; Head-Gordon, M.; Raghavachari, K. Quadratic configuration interaction. A general technique for determining electron correlation energies. J. Chem. Phys. 1987, 87, 5968–5975. [Google Scholar] [CrossRef]
  48. Bader, R.F.W. Atoms in Molecule: A quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  49. Besler, B.H.; Merz, K.M.; Kollman, P.A. Atomic Charges Derived from Semiempirical Methods. J. Comp. Chem. 1990, 11, 431–439. [Google Scholar] [CrossRef]
  50. Chirlian, L.E.; Francl, M.M. Atomic charges derived from electrostatic potentials: A detailed study. J. Comp. Chem. 1987, 8, 894–905. [Google Scholar] [CrossRef]
  51. Brneman, C.M.; Wiberg, K.B. Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis. J. Comp. Chem. 1990, 11, 361–373. [Google Scholar] [CrossRef]
  52. The approach of CHELPG method is shown to be considerably less dependent upon molecular orientation than the original CHELP program.
  53. In the Charges from Electrostatic Potentials using a Grid based method (CHELPG), atomic charges are fitted to reproduce the molecular electrostatic potential (MEP) at a number of points around the molecule.
  54. Martin, F.; Zipse, H. Charge distribution in the water molecule—A comparison of methods. J. Comp. Chem. 2005, 26, 97–105. [Google Scholar] [CrossRef] [PubMed]
  55. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scusera, G.E.; Robb, M.A.; Cheeseman, J.R.; Pople, J.A. Gaussian 98; Revision A.7; Gaussian, Inc.: Pittsburgh, PA, USA, 1998. [Google Scholar]
  56. Ding, H.; Morse, M.D.; Apetrei, C.; Chacaga, L.; Aier, J.P. Resonant two-photon ionization spectroscopy of BNB. J. Chem. Phys. 2006, 125. [Google Scholar] [CrossRef] [PubMed]
  57. Al-Saidi, W.A. Ground state structure of BNB using fixed-node diffusion Monte Carlo. Chem. Phys. Lett. 2012, 543, 41–44. [Google Scholar] [CrossRef]
  58. Becker, S.; Dietze, H.-J. Clusterionenbildung von bornitrid in laserplasmen. Int. J. Mass Spectrum. Ion Process. 1986, 73, 157–166. [Google Scholar] [CrossRef]
  59. Monajjemi, M. Non-covalent attraction of B2N( 2,0) and repulsion of B2N(+)in the BnNn ring: A quantum rotatory due to an external field. Theor. Chem. Acc. 2015, 134, 77, 1–22. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Monajjemi, M.; Bagheri, S.; Moosavi, M.S.; Moradiyeh, N.; Zakeri, M.; Attarikhasraghi, N.; Saghayimarouf, N.; Niyatzadeh, G.; Shekarkhand, M.; Khalilimofrad, M.S.; et al. Symmetry Breaking of B2N(−, 0, +): An Aspect of the Electric Potential and Atomic Charges. Molecules 2015, 20, 21636-21657. https://doi.org/10.3390/molecules201219769

AMA Style

Monajjemi M, Bagheri S, Moosavi MS, Moradiyeh N, Zakeri M, Attarikhasraghi N, Saghayimarouf N, Niyatzadeh G, Shekarkhand M, Khalilimofrad MS, et al. Symmetry Breaking of B2N(−, 0, +): An Aspect of the Electric Potential and Atomic Charges. Molecules. 2015; 20(12):21636-21657. https://doi.org/10.3390/molecules201219769

Chicago/Turabian Style

Monajjemi, Majid, Samira Bagheri, Matin S. Moosavi, Nahid Moradiyeh, Mina Zakeri, Naime Attarikhasraghi, Nastaran Saghayimarouf, Ghorban Niyatzadeh, Marzie Shekarkhand, Mohammad S. Khalilimofrad, and et al. 2015. "Symmetry Breaking of B2N(−, 0, +): An Aspect of the Electric Potential and Atomic Charges" Molecules 20, no. 12: 21636-21657. https://doi.org/10.3390/molecules201219769

APA Style

Monajjemi, M., Bagheri, S., Moosavi, M. S., Moradiyeh, N., Zakeri, M., Attarikhasraghi, N., Saghayimarouf, N., Niyatzadeh, G., Shekarkhand, M., Khalilimofrad, M. S., Ahmadin, H., & Ahadi, M. (2015). Symmetry Breaking of B2N(−, 0, +): An Aspect of the Electric Potential and Atomic Charges. Molecules, 20(12), 21636-21657. https://doi.org/10.3390/molecules201219769

Article Metrics

Back to TopTop