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 “
D∞h” whereas its asymmetrically distorted form has the symmetry of the
C∞v group. Tri-atomic B
2N
(−, 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 B
2N has a symmetric linear regulation in its ground state (
) with an unusually low bending frequency (73 cm
−1) [
1].
They noted that the full valence CASSCF wave function for the
state of linear B–N
(0)–B is unstable with respect to symmetry lowering,
i.e., the
C∞v structure (
r1(BN) ≠
r2(BN)) yields a lower total energy than the
D∞h symmetric structure (
r1(BN) =
r2(BN)) [
10].
In 1992, Knight
et al. found two minimum structures,
i.e., the linear B
2N
(0) and cyclic B
2N [
4]. The UHF/6-31G* level of theory predicts the cyclic
2B
2 state to be at the global minimum, while the correlated methods predict the
state of linear B–N–B to be at the global minimum and the cyclic
2B
2 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
C∞v{
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
D∞h {
}. Obviously describing the electronic structure of the open-shell system (B
2N
(0)) is much more difficult than the closed-shell system of B
2N
(−, +) 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 (
) and (
for anionic and neutral forms, respectively [
2,
13].The ground state of the B
N(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-pV
DZ, cc-pV
TZ and cc-pV
QZ) by J. Paldus [
9]. They showed that the ground state has an asymmetric structure
C∞v with two BN bonds of unequal length.
In an experimental study, the state of
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
[
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
to the ground and lowest excited states of neutral B
2N
i.e.,
and
with a linear symmetry and is assigned to the
and
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 B
2N
(−) 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
[
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
of the triatomic B
2N
(0) radical [
20]. Their results show that B
2N 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
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
D∞h 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
−)/[3
s2
p] 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 B
2N
(−, 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 |
multiplied by the atomic charge yields the charge density |
. 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.
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 B
2N
(−, 0, +) systems. Therefore, based on Walsh prediction [
24], it appears that the linear
p-block molecule holds up quite well for all B
2N
(−, 0, +) ions and the radical. This work has focused on a spontaneous symmetry breaking (SSB) [
25,
26,
27] for B
2N
(−, 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)):
The B
2N
(−, 0, +) forms will be stable by minimization of “
V”. So, it depends on sum of Coulomb repulsions between nuclei and the electronic effective potential:
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:
where
and
and
T is the transposition vector in the 3
N-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(
)”. 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:
where
and charges
qN,
are bound distances and charge distribution of B
2N
(−, 0, +) and
,
and
are defined as:
, in the center of mass of BNB molecules, while:
Depending on the radical, cation and anion forms of BNB, qN × qB multiplication can be positive or negative and consequently, ⎢ 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 . Therefore, the 1.329 and αα' = ββ' = 1.506 can be yielded in the center.
,
and
denote a stable configuration, and V can be expanded around this point up to the second order as Equation (10):
so, the distance between B and N is equal to
and the distance between the two borons is
.
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.
4. Results and Discussion
The radical form of BNB is linear in its ground electronic state () with an orbital occupancy of while the lowest electronically excited state is predicted to be , with an orbital occupancy of .
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
exited state with an orbital occupancy of
(above the
) is subject to the Renner-Teller effect and further exited states depend on the
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/
C∞V (global minima of
Σ
+), are
and
Hartree respectively, while these energies for
C∞V/
D∞h (local minima) are
and
, respectively. Although the energy of
state is near the ground state, this excited configuration does not contribute to the ground state wave function.
In the symmetric
geometry, the unpaired electron of BNB is delocalized, while in the asymmetric
geometry it is localized on either one of the B atoms. Broken symmetry
structures will be stabilized by this interaction relative to the symmetric
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
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.
When the molecule has
symmetry, the real wave-function must transform as an irreducible representation of the
point group. However, when the two B–N bonds are asymmetrically stretched,
and
become near degenerate, while
and
MOs have the same symmetry. Because of this vicinity, the singly excited state of
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 | | Configuration | Configuration | | |
---|
(*Ne) | Total Energy of | Total Energy of , Virtual ** | | |
---|
| ** | | | |
---|
2Σ+ | | | | | |
| | | | | |
(*17e) | | | | | |
| | | | | |
| | | | | |
| | | | | |
| | | =-0.24844 d | 6 d | |
| | | | | |
(*17e) | | | ** | | |
| | | | | |
| | | | | |
| | | | | |
| | | | | |
(*17e) | | | | | |
| | | | | |
| | | | | |
(*17e) | | | | | |
| | | | | |
| | | | | |
(*18e) | | | | | |
| | | | | |
| | | | | |
(*18e) | | | | | |
| | | | | |
(*18e) | | | | | |
| | | | | |
| | [C]= | = | | |
(*18e) | | | | | |
| | | | | |
| | | | | |
| | = | | | |
(*16e) | | | | | |
| | | | | |
| | = | | | |
(*16e) | | | | | |
| | | | | |
| | | | | |
3Σg | | | | | |
(*16e) | | | | | |
| | | | | |
| | | | | |
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 | | | Hybrids Coefficient& | | Atomic Occupancies |
---|
(*Ne) | | | | | |
---|
| | , | | | |
(*17e) | | | | * 0.094 | |
| | | | | |
| | | | * 0.089 | |
| | | | | |
| | | | * 0.029 | |
| | , | | | |
(*18e) | | | | * 0.107 | |
| | | | | |
| | | | *0.021 | |
| | | | | |
| | | | * 0.107 | |
3Σ+ | | , | | | |
(*16e) | | | | * 0.119 | |
| | | | | |
| | | | * 0.028 | |
| | | | | |
| | | | * 0.029 | |
| | , | | | |
(*18e) | | | | * 0.107 | |
| | | | | |
| | | | * 0.064 | |
| | | | | |
| | | | * 0.021 | |
| | | | | |
| | | | | |
Table 3.
Isotropic Fermi contact coupling (IFCC) of B2 in ground () and exited () states.
Table 3.
Isotropic Fermi contact coupling (IFCC) of B2 in ground () and exited () states.
State | | | | |
---|
(*Ne) | | , | | , |
---|
(*17e) | | −29.81, 428.6, 428.6 | | −29.81, 428.6, 428.6 |
| −29.76, 386.2, 469.7 | | −29.81, 428.71, 428.71 |
| −29.11, 348.2, 508.1 | | −29.81, 428.77, 428.77 |
| −28.02, 316.1, 541.9 | | −29.81, 428.95, 428.95 |
| −26.64, 290.3, 570.2 | | −29.83, 429.84, 429.84 |
| −25.12, 270.7, 592.6 | | −29.87, 433.44, 433.44 |
| −23.54, 256.7, 609.5 | | −29.86, 439.68, 439.68 |
| −21.98, 247.6, 621.3 | | −29.61, 448.75, 448.75 |
| −20.46, 242.7, 628.7 | | −28.85, 460.82, 460.82 |
| −20.10, 242.1, 629.8 | | −27.16, 475.86, 475.86 |
| −19.72, 241.6, 630.9 | | −24.06, 493.31, 493.31 |
| 4.35, 60.3, 899.9 | | −18.57, 512.55, 512.55 |
| 4.46, 59.4, 901.2 | | −10.65, 533.75, 533.75 |
| 4.79, 56.6, 904.7 | | |
| | 20.26, 377.8, 377.8 | | |
| 20.26, 381.2, 377.9 | | |
| 20.26, 384.5, 377.9 | | |
| 20.28, 387.8, 377.9 | | |
| 20.30, 391.2, 377.9 | | |
| 20.30, 394.4, 377.9 | | |
| 20.38, 397.8, 377.9 | | |
| 20.44, 401.0, 377.8 | | |
| 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.,
). This is due to the fact that the spin orbital energies are related to the small bending angles of A
1 and A
2 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 (
bending mode
), 1178.64 cm
−1 (
, symmetric stretching “
”) and 2146.42 (
, asymmetric stretching “
”). The IR and Raman intensities for
are 10165.0 and 0.00, respectively, while the
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 (
) and (
{(k − a) and (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
has been indicated as due to a photodetachment from the
to the ground and lowest excited state of neutral B
2N {
and
with a linear symmetry and assigned to the
and
transitions {the
term energy
T0 is 0.785 eV or 6331.77 cm
−1}.
The
orbital is a bonding combination of all
orbitals on all three atoms, while
and
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 −
) is
5 cm
−1 which is near 8.77 cm
−1 (different between
and of strata/stratum). In addition
(exited state above the
) 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
and
orbitals, thereby placing the transition to
with
state much higher in energy. Analysis of the vibronic structure of the
band system shows the transition to
at 19,452 cm
−1. Nevertheless, the nonbonding character of
and
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
, the “
” 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
and
for
and
, respectively, which result the sequences of
for
and
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
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
, 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 () and () states involving the v3 mode.
In
Table 2, the NBO calculation shows that the 2
sN orbital is considered to be primarily core-like, forming the
orbital, though, of course, some mixing of the 2
sN orbital into the other
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
orbital on the nitrogen atom. The
orbital is a bonding combination of the
orbital on the central nitrogen atom with
hybrid orbitals on the two borons. The
orbital is a bonding combination of all
orbitals on all three atoms, while
and
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).
in the ground state for boron atoms varies from 428.6 MHz at
to 241.6 MHz at
for B1 and from 428.6 MHz at
to 630.9 MHz at
for B2 while the nitrogen
varies smoothly from −29.81 MHz (
) to −19.72.0 MHz (
). There is a critical point for
(both in B and N) between
and
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
orbital with most of the electron density on the boron atoms. Symmetry constrains this orbital to have only “
” 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 values for BNB can be obtained from the experimental “” and “”: (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,
and
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
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
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
and
ROHF WFs (with symmetric stretching) have energies separated by a gap as big as 0.5 milli Hartree approximately [
57]. In the case of
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
symmetry of the WF, but this WF would be different from the original
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 B
2N
(−, 0) in our calculations has not been observed, though, for cation, there is a bulge in the curve at 90° in MP
4DQ and MP
4DSQ methods which indicates a cyclic B
2N
+.
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.
At the SCF level, the lowest energy corresponds to a bent molecule with an angle of 100°, however, for the QCISD (T), MP
4DQ, MP
4DSQ 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 B
2N (
2B
2) 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(b
2) of cyclic B
2N and the 1998.4 cm
−1 combination band is the sum of
v1 (a
1), 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 B
2N (
2B2 state) vibrations will display significant vibronic interaction effects. The failure to observe cyclic B
2N 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
clusters ions in laser plasma which resulted in production of BNB cation whereas our calculations resulted in production of cyclic B
2N
+.
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 B
2 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.
Although in some structures of the six resonances, the nitrogen is negative and the two borons are positive, in our ESP calculations for B2 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
state are independent of
during hyperfine coupling calculation. Although these two states
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
symmetry
2Σ
+:
and
6Σ
+:
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 anion structure cannot be predicted to lie above the ground state, it is notable that the lowest stable bent solution for B2 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 (
bending mode “
”) and 1203.42 cm
−1 (
, symmetric stretching “
”) and 1837.50 (
, asymmetric stretching “
”). The IR and Raman intensities for
are 1494.6 and 0.00, respectively, while the
” mode has intensity in Raman (39.97), but zero intensity in IR region. The valence orbital occupancy of ground state (
is:
, while the lowest excited triplet state in the
representation for anion form has calculated (
Table 1) and is predicted to be a
state, lying 2.94 eV {
in
Table 1} above the
state (with the calculated value of 2.7 eV anion with photoelectron spectroscopy of B
2N
(−)).
The energies of and valance occupied MOs for the first triplet of exited state are , = −0.23701 and meanwhile, for are , and . 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.
However, the unrestricted QCISD (T) wave function is considerably spin-contaminated and characterized by a large T1 value. The excited state corresponding to the promotion of an electron from the highest occupied molecular orbital (HOMO), , to the MO is predicted to be considerably above the 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 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 B
2N
(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
, it can be seen in the IFCC
for the nitrogen in range of angles between
(50, −50) (
Figure 3a'). On the other hand, the SB has not been observed for the
=
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
.
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 | | | |
---|
(* Ne) | | IFCC(N,B1,B2) | |
---|
| | | |
---|
2Σ+ | | | |
| | | |
(*17e) | | | |
| | | |
| | | |
| | | |
| | | |
| | | |
(*17e) | | | |
| | | |
| | | |
| | | |
| | | |
(*17e) | | | |
| | | |
| | | 1.0334 |
(*17e) | | | -2.4251 |
| | | 1.3917 |
| | | |
(*18e) | | | |
| | −14.98, 211.91, 211.91 | |
| | | |
(*18e) | | | |
| | | |
| | | |
(*18e) | | | |
| | | |
| | | |
(*18e) | | | |
| | | |
| | | |
(*16e) | | | |
| | | |
| | | |
(*16e) | | | |
| | | |
3Σg | | | |
(*16e) | | | |
| | | |
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
and
are equal and positive, while
is negative when the symmetry of radical is
(in ground state). However, all three components are negative in the exited state for
and their values are in a large scale in comparison to those in ground state while for
they are irregular in the range of ground states. So,
Σ+ interaction should be stronger than
. These results are the same for the anion in ground state of
and the excited state of
and the suitable interaction of
can be predicted (it is
−
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
. 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.