www.nature.

com/scientificreports

OPEN Probing the Fluxional

Bonding Nature of Rapid Cope
rearrangements in Bullvalene
C10H10 and Its Analogs C8H8, C9H10,
and C8BH9
Yuan-Yuan Ma, Miao Yan, Hai-Ru Li, Yan-Bo Wu   , Xin-Xin Tian, Hai-Gang Lu & Si-Dian Li   *

Bullvalene C10H10 and its analogs semibullvalene C8H8, barbaralane C9H10, and 9-Borabarbaralane C8BH9
are prototypical fluxional molecules with rapid Cope rearrangements at finite temperatures. Detailed
bonding analyses performed in this work reveal the existence of two fluxional π-bonds (2 2c-2e π → 2
3c-2e π → 2 2c-2e π) and one fluxional σ-bond (1 2c-2e σ → 1 4c-2e σ → 1 2c-2e σ) in their ground
states and transition states, unveiling the universal π + σ double fluxional bonding nature of these
fluctuating cage-like species. The highest occupied natural bond orbitals (HONBOs) turn out to be
typical fluxional bonds dominating the dynamics of the systems. The 13C-NMR and 1H-NMR shielding
tensors and chemical shifts of the model compound C8BH9 are computationally predicted to facilitate
future experiments.

Chemical bond is the most fundamental and important concept in chemistry. Classical bonds include localized
two-center-two-electron (2c-2e) bonds and delocalized multi-center-two-electron (mc-2e, m ≥ 3) bonds. Our
group predicted the existence of fluxional σ- and π-bonds (FBs) in planar B182− and B19−, half-sandwich KB18−,
tubular Ta@B20−, Ta@B21, and Ta@B22+, and cage-like B39− in four recent papers1–4. Multicenter FBs in these fluc-
tuating boron nanoclusters form and break constantly in concerted mechanisms at room temperatures. It is these
FBs that facilitate the fluxional behaviors of these electron-deficient boron-based nanoclusters which possess
energy barriers lower than the differences of the corresponding zero-point energy corrections. However, boron
nanoclusters are known to be unstable in air and moisture and have hitherto been observed and characterized in
gas phase only. Fluxional bonds in stable systems beyond boron which fluctuate rapidly and reversibly at finite
temperatures remain to be explored in chemistry.
Prototypical fluxional molecules in organic chemistry include the norcaradiene-cycloheptatriene system, var-
ious annulenes, and homotropilidenes. Bridged homotropilidenes with degenerate valence-bond tautomerisms,
such as the cage-like bullvalene C10H10, semibullvalene C8H8, barbaralane C9H10, and 9-borabarbaralane C8BH9,
are of particular interest which exhibit reversible fluxionalities in rapid Cope rearrangements through a transition
state with a bis-homoaromatic array of orbitals5–17. A topological analysis of experimental electron densities of the
ground-state C3v bullvalene was reported in 199618. C10H10, C8H8, and C9H10 have the experimental free energy
barriers of ΔG≠ = 12.8 kcal/mol at 100 °C, 5.5 kcal/mol at −143 °C, and 7.8 kcal/mol at −77 °C in NMR meas-
urements, respectively11, while the model compound C8BH9 has the calculated ΔG≠ = 10.36 kcal/mol at 27 °C17.
Semibullvalene C8H8 has proven to have the lowest fluxional energy barrier, fastest rearrangement rate, and low-
est fluctuating temperature in the series11. Despite their differences in compositions and ground state structures,
cage-like C10H10, C8H8, C9H10, and C8BH9 have similar transition state structures in rapid Cope rearrangements
which have obvious multicenter bonding characteristics. However, the specific bonding patterns and fluxional
bonding nature which facilitate the fluctuating behaviors of these intriguing molecules still remain unknown to
date in both theory and experiments.
We aim to tackle the problem at first-principles theory level in this work. Detailed bonding analyses reveal a
universal bonding pattern with two fluxional π-bonds and one fluxional σ-bond in the ground states and transi-
tion states of the C10H10, C8H8, C9H10, and C8BH9 series, unveiling the σ + π double fluxional bonding nature of

Institute of Molecular Science, Shanxi University, Taiyuan, 030006, China. *email: [email protected]

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 1

www.nature.com/scientificreports/ www.nature.com/scientificreports

these rapidly and reversibly fluctuating species. Their highest occupied natural bond orbitals appear to be typical
fluxional bonds which dominate the fluxional behaviors of the systems in Cope rearrangements. We have also
calculated the 13C-NMR and 1H-NMR shielding tensors and chemical shifts of the model compound C8BH9 to
facilitate future NMR measurements.

Theoretical Procedure.  The ground-state (GS) and transition-state (TS) structures of the concerned spe-
cies were fully optimized at density functional theory (DFT) level of PBE019 with the basis sets of 6–311 + G(d)20.
Frequency checks were performed to make sure all the optimized structures are true GMs or TSs. All the PBE0
structural optimizations and coupled cluster CCSD(T)21–23 single-point calculations in this work were performed
using the Gaussian 09 package24. Detailed bonding analyses were performed on the concerned species using the
adaptive natural density partitioning (AdNDP)25–27 method. The AdNDP approach recovers both the localized
and delocalized bonding elements of the concerned systems and has been successfully applied to a wide range of
nanoclusters and molecules1–4,28–38. Natural bonding orbital (NBO) analyses were performed utilizing the NBO
6.0 program39. The nuclear magnetic resonance (NMR) shielding tensors are calculated using the Continuous Set
of Gauge Transformations (CSGT) method40–42 implemented in Gaussian09 program.

Results and Discussions

Structures and stabilities.  We start from the optimized structures of the GSs and TSs of concerned species
first. As shown in Fig. 1, C3v C10H10 (1), Cs C8H8 (4), Cs C9H10 (7), and Cs C8BH9 (10) as true minima of the systems
possess cage-like structures with the lowest vibrational frequencies of 227, 303, 285, and 194 cm−1 at PBE0 level
(Fig. 1), respectively. They all contain three C-C σ single bonds in the C3 triangle on the top and two C−⃛ C σ + π
double bonds (C3−⃛ C5 and C4−⃛ C6) on the two long edges of the C7 heptagon in the front, with a mirror plane
perpendicular to the paper surface. Their equivalent counterparts GSs′ C10H10 (3), Cs C8H8 (6), Cs C9H10 (9), and
Cs C8BH9 (12) with a C3 triangle at the bottom are degenerate in energy with the GSs discussed above. Obviously,
there exists no delocalized bonding interaction in the true minima GSs and GSs′ in which each carbon atom fol-
lows the octet rule. In contrast, the more open high-symmetry transition states C2v C10H10 (2), C2v C8H8 (5), C2v
C9H10 (8), and C2v C8BH9 (11) with one imaginary frequency at −386i, −336i, −351i, and −456i cm−1 at PBE0,
respectively, all feature two effective C−⃛ C−⃛ C multicenter π-bonding interactions over C1−⃛ C3−⃛ C5 and C2−⃛ C4−⃛
C6 units along the two long edges of the C8 octohedron in the front (with rc1—c3 = rc2—c4 = rc3—c5 = rc4—c6 = 1.39 Å).
They lie 12.9, 9.0, 10.2, and 14.0 kcal/mol higher in energy than their ground states at CCSD(T)//PBE0 level at
298 K, respectively. Such energy barriers appear to be much higher than that previously reported in boron nano-
clusters1–4. This can be qualitatively understood based on the fact that, due to its prototypical electron-deficiency,
boron has the strong propensity to form delocalized σ and π bonds in highly reactive boron nanoclusters with
extremely small energy barriers1–4, while the fluxional processes in 1, 4, 7, and 10 possess much higher energy
barriers because they involve the formations and breakages of C-C interactions in stable organic species. The
C1-C2 single bond with rc1-c2 = 1.53~1.59 Å on the top in the C3v or Cs GSs has been elongated to rc1–c2 = rc5–c6 
= 1.92~2.04 Å in the C2v TSs. The calculated C1--C2 and C5--C6 distances in the C2v TSs appear to be about 0.5 Å
longer than the sum of the single-bond covalent radii of two carbon atoms (rc-c = 1.50 Å)43, indicating that the
C1--C2 and C5--C6 interactions across the two long edges in C2v TSs are much weaker than a usual C-C single
bond. Such C--C distances also appear to be much longer than the C-C single bond (1.579 Å) observed between
the two inverted carbon atoms in propellane44,45. The calculated distances of rc3-c4 = 2.9~3.2 Å in C2v TSs (2, 5, 8,
11) clearly show that there exists no bonding interaction between C3-C4. These transition states with two weak
C--C interactions (C1--C2 and C5--C6) on the top and at the bottom of the C8 octahedron are at the critical
points of Cope intramolecular rearrangements, where the original C1-C2 single σ-bond in the GS is to be broken
while the C5-C6 σ-interaction in GS′ is to be formed simultaneously in the same process and vice versa. The six
carbon atoms (1–6) in the front of the C2v TSs can be divided into two equivalent groups weakly bonded together,
with two effective parallel C−⃛ C−⃛ C multicenter bonds (C1−⃛ C3−⃛ C5 and C2−⃛ C4−⃛ C6) along the two long edges of
the C8 octagon and two weak C--C interactions (C1--C2 and C5-C6) on the top and at the bottom between them.
C10H10, C8H8, C9H10, and C8BH9 possess the calculated free energy barriers of ΔG≠ = 13.32 kcal/mol at 100 °C,
5.94 kcal/mol at −143 °C, 7.86 kcal/mol at −77 °C, and 10.84 kcal/mol at 27 °C at PBE0 level, respectively, well in
line with the corresponding values previously reported for these species at finite temperatures11,17.

AdNDP bonding analyses.  The calculated AdNDP natural bond orbital energy levels of the GSs/GSs′ and
TSs of C10H10 and C8BH9 are comparatively shown in Fig. 2, with that of C8H8 and C9H10 depicted in Fig. S2.
These natural bond orbital energy levels reveal the bonding patterns of the concerned molecules clearly, exhibit
the symmetries of the concerned species perfectly, and show the relative energies of the symmetrically distributed
chemical bonds of the systems directly. The localized AdNDP natural bond orbitals have the advantage over the
delocalized canonical molecular orbitals (CMOs) in providing a pictorial representation of the relative energies
of the concerned chemical bonds and their electron density distributions in space, well in line with chemical
intuitions. As anticipated, C3v GM C10H10 (1) possesses 3 equivalent 2c-2e C-C π bonds with the occupation
numbers of ON = 1.94 along the three long edges as its highest occupied natural bond orbitals (HONBOs) and 3
equivalent 2c-2e C-C σ bonds with ON = 1.93 on the top C3 triangle as the second highest occupied natural bond
orbitals (HONBO-1), together with the remaining 9 2c-2e C-C σ bonds and 10 2c-2e C-H σ bonds to form the GS
in an overall bonding symmetry of C3v (Fig. 3a). From C3v GS to C2v TS, two π-HONBOs of the GS in the front
(2 2c-2e π bonds over C3-C5 and C4-C6) are converted into 2 3c-2e π bonds as HONBO-2 of the C2v TS over C1−⃛
C3−⃛ C5 and C2−⃛ C4−⃛ C6 with ON = 1.96 on the two long edges, one σ-HONBO-1 of the GS in the front (1 2c-2e
σ bond on C1-C2 on the top C3 triangle) is transferred into 1 4c-2e σ-bond with ON = 1.95, the HONBO of the
TS, which is evenly distributed on C1–C4 and C5–C6 with obvious bonding/antibonding characteristics, while
the remaining 1 2c-2e π bond, 11 2c-2e C-C σ bonds, and 10 C-H σ bonds remain basically unchanged. The

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 2

www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 1.  Optimized structures of the ground states (GSs/GSs′) and transition states (TSs) of (a) C10H10, (b)
C8H8, (c) C9H10, and (d) C8BH9, with the lowest vibrational frequencies νmin and relative energies ΔEa indicated
at PBE0 and CCSD(T) levels, respectively. Typical calculated C-C bond lengths are indicated in Å.

delocalized 4c-2e σ-bond is a σ + π mixture between two sets of titled pz-pz pair interactions, with the major
contribution from a head-to-head σ-overlap and minor contribution from a shoulder-by-shoulder π-overlap. An
opposite process occurs from C2v TS to the second minimum C3v GS′. Thus, as clearly shown in Fig. 3a, in a full
fluxional process C3v GS → C2v TS → C3v GS′ → C2v TS′ → C3v GS, C10H10 undergoes a π-fluctuation of 2 2c-2e π
(HONBOs) → 2 3c-2e π (HONBO-2) → 2 2c-2e π′ (HONBOs) → 2 3c-2e π′(HONBO-2) → 2 2c-2e π (HONBOs)
and a σ-fluctuation of 1 2c-2e σ (HONBO-1) → 1 4c-2e σ (HONBO) → 1 2c-2e σ′ (HONBO-1) → 1 4c-2e σ′
(HONBO) → 1 2c-2e σ (HONBO-1) simultaneously in a concerted mechanism. Such a bonding fluctuation

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 3

www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 2.  AdNDP natural bond orbital energy levels and bonding patterns of the ground states (GSs/GSs′) and
transition states (TSs) of (a) C10H10 and (b) C8BH9 at PBE0/6-311 + G (d) level, with the two fluxional π-bonds
and one fluxional σ-bond interlinked by arrowed lines from GS, TS, to GS′.

process occurs randomly in three equivalent directions perpendicular to the three equivalent C7 heptagons
around the C3 molecular axis in both C3v GS and GS′, generating 10!/3 equivalent isomers (~1.2 million) in total
for C10H10, making all the ten H atoms magnetically equivalent with one signal observed in NMR measurements
above 100 °C15.
Similarly, a π-fluctuation of 2 2c-2e π (HONBOs) → 2 3c-2e π (HONBO-1) → 2 2c-2e π′(HONBOs) → 2 3c-2e
π′ (HONBO-1) → 2 2c-2e π (HONBOs) and a σ-fluctuation of 1 2c-2e σ (HONBO-3) → 1 4c-2e σ (HONBO) → 1
2c-2e σ′ (HONBO-3) → 1 4c-2e σ′ (HONBO) → 1 2c-2e σ (HONBO-3) occur in C8BH9 in a full fluxional circle
(Fig. 3b). The two AdNDP π bonds and one AdNDP σ bond of C10H10 and C8BH9 have their origin from the two
CMO π orbitals and one CMO σ orbital in them in both the ground state and transition state, as shown in Fig. S1,
well supporting the AdNDP bonding patterns presented in Figs. 2 and 3. As shown in Fig. S3, C8H8 and C9H10
also exhibit two similar fluxional π bonds (2 2c-2e π → 2 3c-2e π → 2 2c-2e π) and one fluxional σ-bond (1 2c-2e
σ → 1 4c-2e σ → 1 2c-2e σ) in reversible Cope rearrangements. However, different from C3v C10H10 which pos-
sesses three equivalent fluctuating directions, Cs C8H8, Cs C9H10, and Cs C8BH9 can only fluctuate in one direction
in the front C7 heptagon perpendicular to the mirror plane to form two equivalent isomers (Fig. 1). In particular,
it is noticed that the HONBOs which have the highest natural bond orbital energies and highest relative reactivity
in both the GSs and TSs in Cope intramolecular rearrangements are typical fluxional bonds which dominate the
fluxional behaviors of the concerned systems (Figs. 2, 3 and S3). The C3v/Cs GSs all have π-HONBOs, while the
C2v TSs possess σ-HONBOs (Fig. 2 and S2). Such a natural bond orbital energy order renders low stability and
high reactivity to the C2v TSs.

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 4

www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 3.  π- and σ-bonding fluctuations of (a) C10H10 and (b) C8BH9 in a full circle
GS → TS → GS′ → TS′ → GS), with the two fluxional π-bonds and one fluxional σ-bond fluctuating up and
down in opposite directions indicated by red arrows. The ON values represent the calculated occupation
numbers of corresponding bonds.

The calculated electron numbers from the involved C atoms to the fluxional 3c-2e π-bonds and fluxional 4c-2e
σ-bond in Table 1 indicate that, in the C2v transition states of these fluxional molecules, the two central C atoms
(C3 and C4) on the two long edges each contribute one electron to the respective 3c-2e π bond over C1−⃛ C3−⃛ C5
or C2−⃛ C4−⃛ C6, while the four C atoms on the top (C1 and C2) and at the bottom (C5 and C6) each contribute
approximately half an electron. Meanwhile, the four C atoms on the top (C1 and C2) and at the bottom (C5 and
C6) each contribute half an electron to the delocalized 4c-2e σ-bond, forming a half σ-bond on the top over
C1--C2 and a half σ-bond at the bottom over C5--C6. The two separated half σ-bonds over C1-C2 and C5-C6 are
antibonding in nature. As shown in Figs. 3a,b and S3a,b and Video S1, the two fluxional π bonds and one flux-
ional σ-bond in each species fluctuate up and down continuously and reversibly in opposite directions to keep the
balance of σ- and π-bond densities in the fluctuating molecules.
The calculated NBO bond orders of the C2v transition states of these molecules in Fig. S4 also well support the
bonding patterns presented above, with C1--C2, C1−⃛ C3, C3−⃛ C5, and C5--C6 interactions possessing the bond
orders of 0.35, 1.46, 1.46, and 0.35 in C2v C10H10 and 0.43, 1.43, 1.43, and 0.43 in C2v C8BH9, respectively. The
simultaneous formation of both the 2 3c-2e fluxional π-bonds and 1 4c-2e fluxional σ-bond in C2v TS is a natural

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 5

www.nature.com/scientificreports/ www.nature.com/scientificreports

TSs Atoms 3c-2e π-bonds 4c-2e σ-bond

C1(C2) 0.47 0.49
C2v C10H10 C3(C4) 1.02 —
C5(C6) 0.47 0.49
C1(C2) 0.45 0.49
C2v C8BH9 C3(C4) 1.03 —
C5(C6) 0.45 0.49

Table 1.  Calculated electron numbers from specific carbon atoms (1–6) contributed to the respective fluxional
3c-2e π-bonds and fluxional 4c-2e σ-bond in the transition states C2v C10H10 (2) and C2v C8BH9 (11) at PBE0/6-
311 + G(d) level.

C-NMR δ/
13
C-NMR
13 1
H-NMR δ/ 1
H-NMR
ppm Δδ/ppm ppm Δδ/ppm
130.39 (3) 5.07 (3)
54.62 (3) 25.33 (3)
[128.5] [5.70]
130.29 (3) 4.92 (3)
C3v 54.72 (3) 25.48 (3)
[128.3] [5.62]
C10H10
(1) 32.75 (1) 1.41 (1)
152.27 (1) 28.99(1)
[31.0] [2.13]
22.28 (3) 1.38(3)
162.74 (3) 29.02(3)
[21.0] [2.07]
55.23 (2) 129.78 (2) 23.04 (1) 7.36 (1)
60.28 (2) 124.73 (2) 25.49 (2) 4.92 (2)
Cs
C8BH9 144.00 (2) 41.01 (2) 25.75 (2) 4.65 (2)
(10) 144.77 (1) 40.24 (1) 27.85 (1) 2.56 (1)
(298 K)
160.20 (1) 24.82 (1) 28.41 (2) 1.99(2)
29.19 (1) 1.22 (1)

Table 2.  Calculated absolute 13C-NMR and 1H-NMR shielding tensors δ and chemical shifts Δδ with TMS as
internal reference in C3v C10H10 (1) and Cs C8BH9 (10) at PBE0/6-311 + G(d) level. The experimental 13C and
1
H chemical shifts (∆δ) of C10H10 (1) at −59.9 °C and −59.2 °C are cited in square brackets for comparison,
respectively15.

requirement to convert GS and GS′ backward and forward in a continuously and reversibly fluctuating process
(Figs. 3 and S3).

NMR shielding constants.  NMR has proved to be a powerful tool for the determination of the energy
barriers and rate constants of molecules with fluxional bonds in rapid Cope rearrangements11. The calculated
absolute 13C- and 1H-NMR shielding tensors δ and chemical shifts ∆δ relative to tetramethylsilane (TMS) are
tabulated for C3v C10H10 (1) and Cs C8BH9 (10) in Table 2 in ppm. Our calculated 13C and 1H chemical shifts
(∆δ) of bullvalene C10H10 agree well with that measured in NMR experiments at −59.9 °C and −59.2 °C, respec-
tively15. (Table 2). The predicted 13C-NMR spectrum of the GM Cs C8BH9 at 298 K exhibits five kinds of C
atoms with the absolute magnetic shielding tensors of δ = 55.23, 60.28, 144.00, 144.77, and 160.20 ppm in the
intensity ratios of 2:2:2:1:1, respectively, while the corresponding 1H-NMR shielding tensors are calculated to
be at δ = 23.04, 25.49, 25.75, 27.85, 28.41, and 29.19 ppm in the ratios of 1:2:2:1:2:1 (with 1H-B having the lowest
1
H-NMR shielding constant). The B atom in Cs C8BH9 has the calculated 11B-NMR shielding tensor of δ = 24.41
ppm.
In summary, detailed AdNDP bonding analyses performed in this work reveal the existence of two fluxional
π-bonds and one fluxional σ-bond in bullvalene C10H10 and its analogs C8H8, C9H10, and C8BH9. These fluxional
bonds form and break constantly and reversibly in rapid Cope rearrangements at finite temperatures. Such a
universal π + σ double fluxional bonding pattern reflects both the structural characteristics and fluxional bond-
ing nature of these rapidly fluctuating species. Their HONBOs with the highest relative energies and reactivity
belong to typical fluxional bonds which dominate the dynamics of the systems. The fluxional behaviors of these
organic molecules are different in nature from that of the classical fluxional molecules like iron pentacarbonyl
(Fe(CO)5), phosphorus pentafluoride (PF5), and dimethylformamide (CH3)2NC(O)H which undergo Berry
pseudo-rotations46–51 via bond bending, swing, or stretching of localized 2c-2e σ bonds without the breaks or for-
mations of chemical bonds involved in the fluxional process. Explorations of fluxional bonds in more complicated
intramolecular rearrangements known in chemistry are currently underway.

Received: 8 May 2019; Accepted: 31 October 2019;

Published: xx xx xxxx

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 6

www.nature.com/scientificreports/ www.nature.com/scientificreports

References
1. Yan, M. et al. Fluxional Bonds in Planar B19−, Tubular Ta@B20−, and Cage-Like B39−. J. Comput. Chem. 40, 966–970 (2019).
2. Yan, M. et al. Fluxional bonds in quasi-planar B182− and half-sandwich MB18− (M = K, Rb, and Cs). J. Comput. Chem. 40, 1227–1232
(2019).
3. Zhao, X. Y., Luo, X. M., Tian, X. X., Lu, H. G. & Li, S. D. NiB10, NiB11−, NiB12, and NiB13+: Half-Sandwich Complexes with the
Universal Coordination Bonding Pattern of σ Plus π Double Delocalization. J. Clust. Sci. 30, 115–121 (2019).
4. Li, H. R. et al Structural transition in metal-centered boron clusters: from tubular molecular rotors Ta@B21 and Ta@B22+ to cage-like
endohedral metalloborospherene Ta@B22, https://doi.org/10.1007/s10876-019-01646-0. (2019).
5. Uchida, K., Mou, Z. Y., Kertesz, M. & Kubo, T. Fluxional σ-Bonds of the 2,5,8-Trimethylphenalenyl Dimer: Direct Observation of
the Sixfold σ-Bond Shift via a π-Dimer. J. Am. Chem. Soc. 138, 4665–4672 (2016).
6. Schröder, G. Synthese und Eigenschaften von Tricyclo[3.3.2.04.6]deca-2.7.9-trien (Bullvalen). Angew. Chem. 75, 722 (1963).
7. Schröder, G. Chem. Ber. 97, 3140–3149 (1964).
8. Döring, W. Zh. Vsesouzn. Khim. Obshch. 7, 30 (1962).
9. Doering, W. V. E. & Roth, W. R. A rapidly reversible degenerate cope rearrangement: Bicyclo [5.1.0] octa-2,5-diene. Tetrahedron 19,
715–737 (1963).
10. Thermische Umlagerungsreaktionen. Angew. Chem. 75, 27–35 (1963).
11. Cheng, A. K., Anet, F. A. L., Mioduski, J. & Meinwald, J. Determination of the fluxional barrier in semibullvalene by proton and
carbon-13 nuclear magnetic resonance spectroscopy. J. Am. Chem. Soc. 96, 2887–2891 (1974).
12. Gompper, R., Schwarzensteiner, M. L. & Wagner, H. U. Non-isotopic perturbation of fast equilibrating systems. Valence tautomerism
in donor-acceptor-substituted semibullvalenes. Tetrahedron Lett. 26, 611–614 (1985).
13. Zimmerman, H. E. & Grunewald, G. J. The Chemistry of Barrelene. III. A Unique Photoisomerization to Semibullvalene. J. Am.
Chem. Soc. 88, 183–184 (1966).
14. Roth, H. D. & Ablet, C. J. Radical cations of bridged bicyclo[5.1.0]octa-2,5-dienes: open-shell bis-homoaromatic systems. J. Am.
Chem. Soc. 108, 2013–2019 (1986).
15. Oth, J. F. M., Mullen, K., Gilles, J. M. & Schroder, G. Comparison of 13C- and 1H-Magnetic Resonance Spectroscopy as Techniques
for the Quantitative Investigation of Dynamic Processes. The Cope Rearrangement in Bullvalene. Helv. Chim. Acta. 57, 1415–1433
(1974).
16. Schröder, G., Oth, J. F. M. & Merényi, R. Molecules Undergoing Fast, Reversible Valence‐Bond Isomerization. (Molecules with
Fluctuating Bonds). Angew. Chem. Int. Ed. 4, 752–761 (1965).
17. Herberich, G. L., Marx, H. W., Moss, S., Schleyer, Pv. R. & Wagner, T. 9-Borabarbaralanes. Chem. Eur. L. 2, 458–461 (1996).
18. Koritsanszky, T., Buschmann, J. & Luger, P. Topological Analysis of Experimental Electron Densities. 1. The Different C-C Bonds in
Bullvalene. J. Phys. Chem. 100, 10547–10553 (1996).
19. Adamo, C. & Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys.
110, 6158–6170 (1999).
20. Krishnan, R., Binkley, J. S., Seeger, R. & Pople, J. A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave
functions. J. Chem. Phys. 72, 650–654 (1980).
21. Cizek, J. On the Use of the Cluster Expansion and the Technique of Diagrams in Calculations of Correlation Effects in Atoms and
Molecules. Adv. Chem. Phys. 14, 35–89 (1969).
22. Purvis, G. D. III & Bartlett, R. J. A full coupled-cluster singles and doubles model: The inclusion of disconnected triples. J. Chem.
Phys. 76, 1910–1918 (1982).
23. Raghavachari, K., Trucks, G. W., Pople, J. A. & Head-Gordon, M. A fifth-order perturbation comparison of electron correlation
theories. Chem. Phys. Lett. 157, 479–483 (1989).
24. Frisch, M. J. et al. GAUSSIAN 09, Revision D.01, Gaussian, Inc., Wallingford, CT (2009).
25. Zubarev, D. Y. & Boldyrev, A. I. Developing paradigms of chemical bonding: adaptive natural density partitioning. Phys. Chem.
Chem. Phys. 10, 5207–5217 (2008).
26. Zubarev, D. Y. & Boldyrev, A. I. Revealing intuitively assessable chemical bonding patterns in organic aromatic molecules via
adaptive natural density partitioning. J. Org. Chem. 73, 9251–9258 (2008).
27. Lu, T. & Chen, F. W. Multiwfn: a multifunctional wavefunction analyzer. J. Comput. Chem. 33, 580–592 (2012).
28. Wang, L. S. Photoelectron spectroscopy of size-selected boron clusters: from planar structures to borophenes and borospherenes.
Int. Rev. Phys. Chem. 35, 69–142 (2016).
29. Sergeeva, A. P. et al. Understanding boron through size-selected clusters: structure, chemical bonding, and fluxionality. Accounts
Chem. Res. 47, 1349–1358 (2014).
30. Zhai, H. J. et al. Observation of an all-boron fullerene. Nat. Chem. 6, 727–731 (2014).
31. Wang, Y. J. et al. Observation and characterization of the smallest borospherene, B28− and B28. J. Chem. Phys. 144, 064307 (2016).
32. Li, H. R. et al. Competition between quasi-planar and cage-like structures in the B29− cluster: photoelectron spectroscopy and ab
initio calculations. Phys. Chem. Chem. Phys. 18, 29147–29155 (2016).
33. Chen, Q. et al. B33– and B34–: Aromatic Planar Boron Clusters with a Hexagonal Vacancy. Eur. J. Inorg. Chem. 38–39, 4546–4551
(2017).
34. Chen, Q. et al. Planar B38− and B37− clusters with a double-hexagonal vacancy: molecular motifs for borophenes. Nanoscale 9,
4550–4557 (2017).
35. Romanescu, C., Galeev, T. R., Li, W. L., Boldyrev, A. I. & Wang, L. S. Transition-Metal-Centered Monocyclic Boron Wheel Clusters
(M©Bn): A New Class of Aromatic Borometallic Compounds. Acc. Chem. Res. 46, 350–358 (2013).
36. Moreno, D. et al. B182−: a quasi-planar bowl member of the Wankel motor family. Chem. Commun. 50, 8140–8143 (2014).
37. Jian, T. et al. Probing the Structures and Bonding of Size-Selected Boron and Doped-Boron Clusters. Chem. Soc. Rev. 48, 3550–3591
(2019).
38. VandeVondele, J. et al. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves
approach. Comput. Phys. Commun. 167, 103–128 (2005).
39. Glendening, E. D., et al, NBO 6.0, Theoretical Chemistry Institute, University of Wisconsin, Madison (2013).
40. Keith, T. A. & Bader, R. F. W. Calculation of magnetic response properties using atoms in molecules. Chem. Phys. Lett. 194, 1–8
(1992).
41. Keith, T. A. & Bader, R. F. W. Calculation of magnetic response properties using a continuous set of gauge transformations. Chem.
Phys. Lett. 210, 223–231 (1993).
42. Cheeseman, J. R., Trucks, G. W., Keith, T. A. & Frisch, M. J. A comparison of models for calculating nuclear magnetic resonance
shielding tensors. J. Chem. Phys. 104, 5497–5509 (1996).
43. Pyykko, P. & Atsumi, M. Molecular Double-Bond Covalent Radii for Elements Li-E112. Chem. Eur. J. 15, 12770–12779 (2009).
44. Messerschmidt, M. et al. Electron Density and Bonding at Inverted Carbon Atoms: An Experimental Study of a [1.1.1] Propellane
Derivative. Angew. Chem. Int. Ed. 44, 3925–3928 (2005).
45. Mata, R. A. & Suhm, M. A. Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction? Angew. Chem. Int.
Ed. 56, 11011–11018 (2017).
46. Berry, R. S. Correlation of rates of intramolecular tunneling processes, with application to some group V compounds. J. Chem. Phys.
32, 933–938 (1960).

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 7

www.nature.com/scientificreports/ www.nature.com/scientificreports

47. Gutowsky, H. S. & Holm, C. H. Rate processes and nuclear magnetic resonance spectra. II. Hindered internal rotation of amides. J.
Chem. Phys. 25, 1228–1234 (1956).
48. Cass, M. E., Hii(Mimi), K. K. & Rzepa, H. S. Mechanisms that Interchange Axial and Equatorial Atoms in Fluxional processes:
Illustration of the Berry Pseudorotation, the Turnstile and the Lever Mechanisms via animation of transition state normal
vibrational modes. J. Chem. Educ. 83, 336 (2006).
49. Kramer, G. M. CH5+ Stability and Mass Spectrometry. Science 286, 1051 (1999).
50. White, E. T., Tang, J. & Oka, T. CH5+: the infrared spectrum observed. Science 284, 135–137 (1999).
51. Marx, D. & Parrinello, M. CH5+: The cheshire cat smiles. Science 284, 59–61 (1999).

Acknowledgements
This work was supported by the National Natural Science Foundation of China (21720102006 to S.-D. Li).

Author contributions
S.-D. Li, H.-G. Lu and Y.-B. Wu designed the research. Y.-Y. Ma, M. Yan, H.-R. Li and X.-X. Tian performed
the calculations. All authors contributed to the interpretation and discussion of the data and participated in the
preparation of the manuscript.

Competing interests
The authors declare no competing interests.

Additional information
Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-53488-5.
Correspondence and requests for materials should be addressed to S.-D.L.
Reprints and permissions information is available at www.nature.com/reprints.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019

Scientific Reports | (2019) 9:17074 | https://doi.org/10.1038/s41598-019-53488-5 8