Triple-Bonded Boron≡Phosphorus Molecule: Is That Possible? - ACS

Jan 4, 2018 - The effect of substitution on the potential energy surfaces of RB≡PR (R = H, F, OH, SiH3, and CH3) is studied using density functional...
0 downloads 7 Views 2MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article Cite This: ACS Omega 2018, 3, 76−85

Triple-Bonded BoronPhosphorus Molecule: Is That Possible? Jia-Syun Lu,† Ming-Chung Yang,† and Ming-Der Su*,†,‡ †

Department of Applied Chemistry, National Chiayi University, Chiayi 60004, Taiwan Department of Medicinal and Applied Chemistry, Kaohsiung Medical University, Kaohsiung 80708, Taiwan



S Supporting Information *

ABSTRACT: The effect of substitution on the potential energy surfaces of RBPR (R = H, F, OH, SiH3, and CH3) is studied using density functional theories (M06-2X/Def2-TZVP, B3PW91/Def2TZVP, and B3LYP/LANL2DZ+dp). There is significant theoretical evidence that RBPR compounds with smaller substituents are fleeting intermediates, so they would be difficult to be detected experimentally. These theoretical studies using the M06-2X/Def2-TZVP method demonstrate that only the triply bonded R′BPR′ molecules bearing sterically bulky groups (R′ = Tbt (=C6H2-2,4,6-{CH(SiMe3)2}3), SiMe(SitBu3)2, Ar* (=C6H3-2,6-(C6H2-2,4,6-i-Pr3)2), and SiiPrDis2) are significantly stabilized and can be isolated experimentally. Using the simple valence-electron bonding model and some sophisticated theories, the bonding character of R′BPR′ should be viewed as R′BI PR′. The present theoretical observations indicate that both the electronic and the steric effect of bulkier substituent ligands play a key role in making triply bonded R′BPR′ species synthetically accessible and isolable in a stable form.

I. INTRODUCTION Molecules that do not contain a metal atom but which feature a triple bond are very rare, probably because the nonmetallic elements do not have valence d orbitals. Therefore, it is quite difficult to use nonmetallic atoms to prepare or synthesize triply bonded molecules.1−8 The most famous example of a triply bonded nonmetallic compound is ethyne, HCCH. During the last two decades, heavy homonuclear acetylene analogues, RE14E14R (E14 = Si, Ge, Sn, and Pb), have been successfully isolated and structurally characterized.9−18 Heteronuclear alkyne analogues, RE14E′14R (E′14 = group 14 elements), have also been the subject of many studies in attempts to increase the understanding of the structural features and their high reactivity.19−31 However, the chemistry of the triply bonded RE13E15R (E13 = group 13 atoms and E15 = group 15 atoms) molecules, which are isoelectronic with acetylene, has not been the subject of as much study as the RE14E14R species. In fact, experimental studies of molecules that feature a E13E15 triple bond are few.32−35 Significantly little is known about the structural and electronic properties of these triply bonded RE13E15R systems, but there have been many studies of these new types of acetylene analogues in many disciplines, such as structural chemistry, catalysis, inorganic chemistry, and organometallic chemistry.1−8,32−38 Since the pioneering work of Coates and Livingstone on the preparation of phosphinoborane (Ph2BPPh2) in 1961,39 many new examples of the B−P single bond have been found.40,41 These are reviewed in the literature.42−44 In addition to that of the molecules that feature a B−P single bond, the synthesis of compounds that feature the BP double bond character has attracted the interest of inorganic and organometallic chemists.42−45 Synthesis of these molecules, which have either © 2018 American Chemical Society

the B−P single bond or the BP double bond, is of interest because they are respectively isoelectronic with ethane and ethylene, from the valence-electron viewpoint. It is therefore expected that the structural, physical, and chemical properties of BP species are all similar to those of CC. However, there is no such similarity of reactivity between CC and BP molecules.46 The polarity and the relative weakness of B−P bonds make the BP molecules much more reactive, in additions or substitutions,46 than comparable CC compounds. Although molecules that feature either a B−P single or a BP double bond have been the subject of many experimental and theoretical studies of their structures and reactivity,47−73 compounds that feature a BP triple bond remain elusive. To the authors’ best knowledge, there have been no experimental or theoretical studies on this subject. This is the first theoretical examination to explore the possibility of forming triply bonded RBPR molecules with respect to small and large substituents. For details about the methodology and computational results, see the Supporting Information. The substituent groups that are studied are H, F, OH, SiH3, CH3 and SiiPrDis2, Ar* (=C6H3-2,6-(C6H2-2,4,6-iPr 3 ) 2 ), SiMe(SitBu 3 ) 2 , and Tbt (=C 6 H 2 -2,4,6-{CH(SiMe3)2}3)74−76 (see Scheme 1). The conclusions of this theoretical study will give experimental synthetic chemists a useful insight into a new research area: the BP triply bonded species. Received: October 4, 2017 Accepted: December 18, 2017 Published: January 4, 2018 76

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega

together to produce a triply bonded RBPR compound in the singlet ground state. One is [R−B]1 + [R−P]1 → [RBPR]1 (path [1]) and the other is [R−B]3 + [R−P]3 → [RBPR]1 (path [2]). Their valence-bond orbital interactions are schematically shown in Figure 1. Figure 1 shows that the pathway that is used to yield RB PR depends on the promotion energies of the fragments. As will be shown in the next section, the present density functional theory (DFT) calculations show that the R−B component has a singlet ground state, whereas the R−P unit has a triplet ground state. Therefore, if the excitation energy from the triplet ground state to the singlet excited state for the R−P unit is smaller than that from the singlet ground state to the triplet excited state for the R−B component, then R−B and R−P follow path [1] to form a triply bonded RBPR molecule. However, if the promotion energy from the singlet ground state to the triplet excited state for the R−B fragment is smaller than that from the triplet ground state to the singlet excited state for the R−P fragment, then R−B and R−P follow path [2] to produce the RBPR species with a BP triple bond. As seen for path [1] in Figure 1, the triple bond in the RB PR compound can be viewed as two donor−acceptor π bonds and one donor−acceptor σ bond. Accordingly, the bonding property of RBPR can be regarded as RB PR. Similarly, as seen for path [2] in Figure 1, the triple bond in the RBPR compound can be described as one customary σ bond, one customary π bond, and one donor−acceptor π bond, so the bonding character of RBPR can be considered as RB PR. Using these theoretical analyses, the bonding properties of the RBPR species are interpreted in the following section.

Scheme 1

II. GENERAL CONSIDERATIONS According to the valence bond theory,77 the most convenient way to understand the chemical bonding natures of the BP triple bond is to partition the RBPR molecule into two units: R−P and R−B (see Figure 1). From the spin viewpoint, there are two reaction routes that combine the two fragments

III. RESULTS AND DISCUSSION III.I. Small Ligands on Substituted RBPR. The small groups, R (=H, F, OH, SiH3, and CH3), are used to examine the effect of substituents on the stability of triply bonded RB PR molecules. Three types of DFT calculations (M06-2X/ Def2-TZVP, B3PW91/Def2-TZVP, and B3LYP/LANL2DZ +dp) are also used to determine their reliability. Some physical properties and the important geometrical parameters of RB PR are collected in Table 1. From Table 1, the DFT computations anticipate that the BP triple bond lengths (Å) are between 1.713 and 1.777 (M06-2X/Def2-TZVP), 1.714 and 1.781 (B3PW91/Def2TZVP), and 1.725 and 1.793 (B3LYP/LANL2DZ+dp). According to the available experimental data, the bond distances (Å) for the B−P single bond and the BP double bond are 1.94−2.015 and 1.786−1.859,42−46,78 respectively. The bond lengths of the BP triple bonds that are obtained from DFT calculations are shorter than those that are acquired experimentally.42−46,78 Additionally, the DFT computations given in Table 1 indicate that the promotion energy from the singlet ground state to the triplet state for the R−B unit is calculated to be at least 22 kcal/mol. On the other hand, the excitation energy from the triplet ground state to the singlet state for the R−P fragment is estimated to be at least 15 kcal/mol. Namely, the theoretical evidence demonstrates that the formation of the triply bonded RBPR molecule with small ligands should follow path [1], as shown in Figure 1. In other words, the bonding character of a triply bonded RBPR molecule possessing small ligands can be represented as RB PR. Because the covalent radii of B and P are reported to be 82 and 106 pm,79 respectively, the overlapping populations between

Figure 1. Valence-bond bonding paths [1] and [2] for the triply bonded RBPR compound. ΔE1 = E(triplet state for R−P) − E(singlet state for R−P) and ΔE2 = E(triplet state for R−B) − E(singlet state for R−B). 77

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega Table 1. Important Geometrical Parametersa R BP (Å)

∠R−P−B (deg)

∠B−P−R (deg)

∠R−P−B−R (deg)

QBb

QPc

ΔEB for R−B (kcal/mol)d

ΔEP for R−P (kcal/mol)e

HOMO−LUMO (kcal/mol)

BE (kcal/mol)f

WBIg

F

OH

H

CH3

SiH3

1.777 (1.781) [1.793] 173.5 (173.1) [174.3] 98.82 (99.24) [96.42] 180.0 (180.0) [180.0] 0.6484 (0.6101) [0.6673] 0.3809 (0.3772) [0.3842] 81.01 (73.60) [73.97] −28.91 (−33.35) [−31.76] 217.2 (210.1) [194.7] 127.2 (120.8) [121.1] 1.868 (1.876) [1.874]

1.770 (1.775) [1.784] 177.2 (176.3) [177.5] 97.56 (98.51) [96.19] 174.1 (178.5) [180.0] 0.5520 (0.5181) [0.5992] 0.2947 (0.2875) [0.2829] 68.97 (62.12) [64.90] −17.53 (−21.29) [−20.24] 219.8 (208.8) [200.3] 123.4 (117.7) [117.4] 1.863 (1.871) [1.862]

1.718 (1.723) [1.736] 176.9 (179.6) [178.9] 72.53 (77.20) [76.40] 180.0 (180.0) [180.0] 0.1535 (0.1664) [0.1097] −0.0975 (−0.0618) [−0.0281] 28.74 (27.65) [25.39] −30.75 (−35.49) [−33.16] 136.8 (128.3) [121.7] 140.7 (145.9) [128.6] 2.231 (2.176) [2.179]

1.726 (1.728) [1.743] 172.1 (173.5) [173.7] 101.9 (104.6) [99.30] 180.0 (178.1) [179.6] 0.4269 (0.4007) [0.4689] −0.0564 (−0.0633) [−0.0386] 38.47 (32.69) [36.79] −26.47 (−30.26) [−29.21] 228.6 (219.0) [215.1] 129.7 (124.8) [121.9] 2.078 (2.047) [2.056]

1.713 (1.714) [1.725] 174.7 (173.3) [175.0] 72.21 (73.04) [76.42] 180.0 (180.0) [179.7] 0.3578 (0.3283) [0.2387] −0.0484 (−0.0561) [−0.1344] 22.33 (21.77) [22.28] −5.80 (−8.67) [−14.46] 156.2 (144.4) [144.9] 139.5 (132.9) [146.0] 2.220 (2.184) [2.189]

a The Wiberg bond index (WBI), the natural charge densities (QP and QB), the highest occupied molecular orbital (HOMO)−lowest unoccupied molecular orbital (LUMO) energy gaps, the singlet−triplet energy splitting (ΔEB and ΔEP), and the binding energies (BEs) for RBPR using the B3PW91/Def2-TZVP (in round brackets), M06-2X/Def2-TZVP, and B3LYP/LANL2DZ+dp (in square brackets) levels of theory. bThe natural charge density on the boron atom. cThe natural charge density on the phosphorus atom. dΔEB = E(triplet state for R−B) − E(singlet state for R−B). e ΔEP = E(triplet state for R−P) − E(singlet state for R−P). fBE = E(singlet state for R−B) + E(singlet state for R−P) − E(singlet state for RB PR). gThe Wiberg bond index (WBI) for the BP bond: see refs 80, 81.

not stabilize triply bonded RBPR species, sterically bulky ligands (R′), such as SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*74−76 (Scheme 1), are used to determine the possibility of forming a triply bonded R′BPR′ molecule. Recently, Liptrot and Power demonstrated that London dispersion forces, which are nonvalent interactions between the bulkier ligands, significantly influence the stability and structure of sterically congested inorganic systems.82 Therefore, the M06-2X/Def2TZVP method83 is used in this work to obtain precise results. Similarly to RBPR that are substituted with small ligands, as discussed previously, the 1,2-migration reactions, R′BPR′ → R2′BP: and R′BPR′ → :BPR2′ (Scheme 2), are used to determine the feasibility of producing triply bonded R′B PR′ compounds, using the M06-2X/Def2-TZVP method. Clearly, as seen in Table 2, the M06-2X/Def2-TZVP data demonstrates that because of a large steric effect the doubly bonded isomers (:BPR2′ and R2′BP:) are all respectively higher in energy by at least 89 and 73 kcal/mol than their corresponding triply bonded R′BPR′ molecule. There is strong theoretical evidence that sterically crowded substituents protect the weak BP triple bond.

boron and phosphorus should be small. Indeed, as seen in Table 1, the Wiberg bond index (WBI)80,81 for the BP bond is predicted to be less than 2.2. However, the WBI for the C C bond in acetylene is shown to be 2.99. The 1,2-migration reactions are used to determine the kinetic stability of the RBPR species that are substituted with small ligands at the three DFT levels of theory. The computational results are listed in Figure 2, which shows that the 1,2-shift reactions have two reaction routes: RBPR → TS1 → R2B P: and RBPR → TS2 → :BPR2. From Figure 2, the three DFT calculations all demonstrate that regardless of whether the small substituent groups are electronegative or electropositive, the triply bonded RBPR molecules are neither thermodynamically nor kinetically stable on the 1,2-shift potential energy surfaces. These theoretical observations provide strong evidence that all of the RBPR triply bonded molecules possessing small groups should be unstable and must spontaneously rearrange to other doubly bonded isomers. III.II. Large Ligands on Substituted R′BPR′. Because these theoretical conclusions show that small ligands (R) do 78

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega

Figure 2. Relative Gibbs free energy surfaces for RBPR (R = H, F, OH, SiH3, and CH3). These energies are calculated in kcal/mol and are calculated at the M06-2X/Def2-TZVP, B3PW91/Def2-TZVP, and B3LYP/LANL2DZ+dp levels of theory. For details, see the text and Table 1.

Figure 1, the bonding character of a triply bonded R′BPR′ compound featuring bulky ligands can be considered to be R′B PR′. In other words, its BP triple bond consists of one traditional σ bond, one donor−acceptor π bond, and one traditional π bond. However, for the cases of Ar*, our computational values shown in Table 2 indicate that ΔEB′ for the B−R′ fragment (35 kcal/mol) is larger than ΔEP′ for the P−R′ fragment (at least 31 kcal/mol). According to Figure 1, the triply bonded (Ar*)BP(Ar*) molecule can be considered (Ar*) (path [1]). to be (Ar*)BP From Figure 1, it is noteworthy that the lone-pair orbital of the R′−P fragment includes both s and p valence orbitals. This means that the overlap coupling between the lone-pair orbital of R′−P and the pure p orbital of R′−B is feeble. As a result, the bond order of the BP triple bond must not be large. Indeed, the M06-2X/Def2-TZVP results in Table 2 demonstrate that the WBI values for the BP triple bond are

As seen in Table 2, the triple bond lengths for these triply bonded R′BPR′ molecules are in the range of 1.74−2.02 Å, which are somewhat larger than the BP triple bond distances of the substituted RBPR compounds bearing small ligands shown in Table 1. The reason for the longer bond lengths for the R′BPR′ species should be the steric effects. Table 2 shows that the singlet−triplet energy splitting (ΔEB′) for the R′−B unit is at least 25 kcal/mol, but the modulus of ΔEP′ for the R′−P moiety is predicted to be at least 30 kcal/ mol on the basis of the M06-2X/Def2-TZVP method. Namely, it is easier for the R′−B unit to jump from the singlet ground state to the triplet state than for the R′−P monomer, which is promoted from the triplet ground state to the singlet state. Therefore, the two bulkily substituted fragments, R′−B and R′−P, follow path [2] (Figure 1) to yield a triply bonded R′BPR′ molecule. That is, [R′−B]3 + [R′−P]3 → [R′B PR′]1. According to the valence-orbital interactions shown in 79

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega

is No. 224 (HOMO) orbital, demonstrating that the former donates electrons (0.216e) to the latter mainly through the HOMO orbital. In addition, the largest contributions from a (SiMe(SitBu3)2)−B component to a (SiMe(SitBu3)2)−P component is No. 223 (HOMO-1) orbital, indicating that the latter donates electrons (0.0682e) to the former principally through the HOMO-1 orbital. After considering the electron donations from various orbitals of both (SiMe(SitBu3)2)−P and (SiMe(SitBu3)2)−B, one may obtain the net amount of electron transfer, which is calculated to be negative (−0.226), suggesting that a (SiMe(SitBu3)2)−P unit donates more electrons to a (SiMe(SitBu3)2)−B unit. This conclusion based on the CDA method agrees well with the valence-electron bonding model (Figure 3). Again, the theoretical information indicates that the bonding nature of R′BPR′ can be represented as R′B PR′. The natural bond orbital (NBO)80,81 and natural resonance theory (NRT)85−87 are used to analyze the electronic densities of (SiMe(SitBu3)2)BP(SiMe(SitBu3)2), (Tbt)BP(Tbt), (SiiPrDis2)BP(SiiPrDis2), and (Ar*)BP(Ar*), whose theoretical results are collected in Table 4. On the basis of the WBI values in Table 4, the bond orders for the BP bond are 2.39, 1.96, 2.15, and 1.97 for SiMe(SitBu3)2, Tbt, SiiPrDis2, and Ar*, respectively. This data is similar to the values that are calculated using NRT, of 2.55, 2.47, and 2.27, respectively. This theoretically computed data strongly suggests that steric hindrance shelters the intramolecular rearrangements and gives the central BP bond a greater bond order. The NBO/NRT analyses in Table 4 show that the bulkily substituted R′BPR′ compounds that are studied in this work all have an analogous electronic structure. For instance, the triple bond in (SiMe(SitBu3)2)BP(SiMe(SitBu3)2) is composed of a σ component, which is slightly polarized toward the P element (57%) and two nondegenerate π components that are heavily polarized toward the P center (π⊥, 76% and π∥, 63%).88 Both the σ and π NBO orbitals are equally filled, with approximately 1.98 (σ), 1.92 (π⊥), and 1.93 (π∥) electrons, which is a feature of localized bonds. Figure 3 schematically

Scheme 2

calculated to be 2.388 (SiMe(SitBu3)2), 2.152 (SiiPrDis2), 1.963 (Tbt), and 1.966 (Ar*), respectively. Again, the theoretical findings show that the bulky ligands can effectively protect the central weak BP triple bond and increase its bond order. The charge decomposition analysis (CDA), which was published by Dapprich and Frenking,84 has been utilized in this work to explain the interactions between two components (i.e., (SiMe(SitBu3)2)−B and (SiMe(SitBu3)2)−P) of a (SiMe(SitBu3)2)BP(SiMe(SitBu3)2) molecule. On the basis of the M06-2X/Def2-TZVP method, the computational data are given in Table 3. From Table 3, the largest contribution from a (SiMe(SitBu3)2)−P fragment to a (SiMe(SitBu3)2)−B fragment

Table 2. Bond Lengths (Å), Bond Angles (deg), Singlet−Triplet Energy Splitting (ΔEB′ and ΔEP′), Natural Charge Densities (QB′ and QP′), Binding Energies (BEs), the Wiberg Bond Index (WBI), HOMO−LUMO Energy Gaps, and Some Reaction Enthalpies for R′BPR′ at the M06-2X/Def2-TZVP Level of Theorya R′

SiMe(SitBu3)2

SiiPrDis2

Tbt

Ar*

BP (Å) ∠R′−B−P (deg) ∠B−P−R′ (deg) ∠R′−B−P−R′ (deg) QB′b QP′c ΔEB′ for R′−B (kcal/mol)d ΔEP′ for R′−P (kcal/mol)e HOMO−LUMO (kcal/mol) BE (kcal/mol)f ΔH1 (kcal/mol)g ΔH2 (kcal/mol)g WBIh

1.736 157.2 122.0 174.7 −0.2574 −0.1824 25.92 −33.10 73.76 89.54 73.75 80.53 2.388

2.021 166.0 112.5 165.5 −0.1395 −0.3922 24.86 −37.47 43.44 90.37 86.65 77.67 2.152

2.023 164.4 121.3 168.9 0.2718 0.2260 28.76 −29.74 47.10 85.42 87.89 101.7 1.963

2.021 166.6 123.3 169.5 0.3520 0.2522 34.64 −30.52 41.60 71.43 87.59 88.01 1.966

See also Scheme 2. bThe natural charge density on the boron atom. cThe natural charge density on the phosphorus atom. dΔEB′ (kcal mol−1) = E(triplet state for R′−B) − E(singlet state for R′−B). eΔEP′ (kcal mol−1) = E(triplet state for R′−P) − E(singlet state for R′−P). fBE (kcal mol−1) = E(triplet state for R′−B) + E(triplet state for R′−P) − E(singlet for R′BPR′). gSee Scheme 2. hThe Wiberg bond index (WBI) for the BP bond: see refs 80, 81. a

80

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega

Table 3. Charge Decomposition Analysis (CDA) Resultsb for the R′BPR′ (R′ = SiMe(SitBu3)2) System Based on M06-2X Orbitalsa

HOMO LUMO sumc

orbital

occupancy

A

B

A−B

M

214 215 216 217 218 219 220 221 222 223 224 225 226

2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 448.000000

0.001584 0.008839 0.002683 0.000235 −0.000147 0.003282 0.004845 −0.000176 0.001401 0.006280 −0.021614 0.000000 0.000000 0.063156

0.051729 0.005614 0.004550 0.004887 0.000075 0.000852 0.002735 0.004588 0.018757 0.068203 −0.039334 0.000000 0.000000 0.289323

−0.050145 0.003225 −0.001867 −0.004653 −0.000222 0.002430 0.002110 −0.004764 −0.017356 −0.061923 0.017720 0.000000 0.000000 −0.226167

−0.047241 −0.026216 −0.017145 −0.009183 −0.004893 −0.021957 −0.003281 −0.012462 −0.048796 −0.111756 −0.026170 0.000000 0.000000 −0.189901

a

The A term shows the number of electrons donated from the R′−B fragment to the R′−P fragment. The B term shows the number of electrons back-donated from the R′−P fragment to the R′−B fragment. The M term shows the number of electrons involved in repulsive polarization. bFor clearness, only the A, B, and M terms are listed for HOMO (No. 224) − 10 ∼ LUMO + 2. cSummation of contributions from all unoccupied and occupied orbitals.

Figure 3. Natural BP π bonding orbitals (π⊥ and π∥ for (i) and (ii), respectively) for (SiMe(SitBu3)2)BP(SiMe(SitBu3)2), compared with path [2] in Figure 1.

shows the two nondegenerate π bonding orbitals (π⊥ and π∥). These are similar to those in path [2], shown in Figure 1.

bonded RBPR species. In terms of these substituted RBPR systems, the following conclusions are drawn. (1) Calculations at three different levels of DFT are in excellent agreement. All three DFT computations demonstrate that regardless of the electronegativity of the substituents that are attached,88 triply bonded RB PR molecules, which bear the small substituents, are unstable and readily undergo intramolecular rearrangement to yield the correspondingly doubly bonded isomers. In other words, the triply bonded RBPR

IV. CONCLUSIONS This study uses theoretical computations to determine the possibility of producing triply bonded RBPR molecules using different types of substituents. The present work provides the first theoretical examination of the effect of substituents on the stability and electronic bonding structures of these triply 81

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega

Table 4. Natural Bond Orbital (NBO) and Natural Resonance Theory (NRT) Analysis for R′BPR′ Molecules Possessing Bulky Ligands (R′ = SiMe(SitBu3)2, Tbt, SiiPrDis2, and Ar*) at the M06-2X/Def2-TZVP Level of Theorya,b NBO analysis

NRT analysis

R′BPR′

WBI

occupancy

R′=SiMe(SitBu3)2

2.39

σ: 1.98

σ: 0.6538 B (sp1.20) + 0.7567 P (sp0.80)

π⊥: 1.92

π⊥: 0.4861 B (sp17.12) + 0.8739 P (sp9.30)

π∥: 1.93

π∥: 0.6062 B (sp99.99) + 0.7953 P (sp1.00)

σ: 1.93

σ: 0.6141 B (sp2.20) + 0.7892 P (sp1.20)

π⊥: 1.86

π⊥: 0.5402 B (sp2.52) + 0.8415 P (sp4.70)

π∥: 1.94

π∥: 0.5208 B (sp99.99) + 0.8537 P (sp1.00)

σ: 1.92

σ: 0.7076 B (sp0.94) + 0.7066 P (sp11.73)

π⊥: 1.86

π⊥: 0.5382 B (sp63.76) + 0.8428 P (sp68.98)

π∥: 1.77

π∥: 0.5317 B (sp99.99) + 0.7737 P (sp1.00)

σ: 1.93

σ: 0.7171 B (sp0.74) + 0.6969 P (sp12.85)

π⊥: 1.84

π⊥: 0.5842 B (sp99.99) + 0.8516 P (sp99.99)

π∥: 1.76

π∥: 0.5402 B (sp99.99) + 0.8426 P (sp1.00)

R′ = SiiPrDis2

R′ = Tbt

R′ = Ar*

2.15

1.96

1.97

hybridization

molecules possessing the small ligands are so fragile that their experimental detection is unlikely. (2) Only sterically large substituents (R′) significantly stabilize triply bonded R′BPR′ compounds. These theoretical observations demonstrate that these bulky ligands hinder the intramolecular rearrangement and obstruct the intermolecular associations. Accordingly, the present theoretical proof demonstrates that it should be possible to synthesize molecules that feature a BP triple bond possessing both bulkier ligands. (3) On the basis of the valence electron bonding model and some theoretical bonding analyses (i.e., NBO and NRT), it is concluded that the bonding nature of triply bonded R′BPR′ species that feature bulky groups is R′B PR′. Namely, the BP triple bond comprises one donor−acceptor π bond, one conventional σ bond, and one conventional π bond. Experimentalists are encouraged to confirm these theoretical predictions.





resonance weight

polarization

total/covalent/ionic

42.74% 57.26% 23.63% 76.37% 36.74% 63.26% 37.71% 62.29% 29.19% 70.81% 27.13% 72.87% 50.07% 49.93% 28.97% 71.03% 26.44% 73.56% 51.43% 48.57% 27.48% 72.52% 25.74% 74.26%

2.19/1.52/0.67

B−P: 6.07% BP: 68.87% BP: 25.06%

2.52/1.23/1.29

B−P: 10.18% BP: 73.31% BP: 16.51%

2.12/1.52/0.60

B−P: 7.29% BP: 73.23% BP: 19.48%

2.11/1.57/0.54

B−P: 2.63% BP: 83.80% BP: 13.57%

(B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P) (B) (P)

B3PW91/Def2-TZVP, and B3LYP/LANL2DZ+dp levels of theory (PDF)

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Jia-Syun Lu: 0000-0001-6776-0889 Ming-Der Su: 0000-0002-5847-4271 Author Contributions

J.-S.L. and M.-C.Y. conducted all of the theoretical computations and analyzed the results. M.-D.S. supervised the research activities and contributed to the manuscript preparation. Three authors regularly discussed the progress of the research, reviewed the manuscript, and gave approval for the final version. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors are grateful to the National Center for HighPerformance Computing of Taiwan for generous amounts of computing time, and the Ministry of Science and Technology of Taiwan for the financial support. Special thanks are also due to reviewers 1, 2, and 3 for very helpful suggestions and comments.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.7b01480. The optimized geometries and the absolute energies (in Hartrees) for all of the points on the potential energy surfaces of RBPR (R = F, OH, H, CH3, and SiH3) and R′InSbR′ (R′ = SiMe(SitBu3 ) 2, SiiPrDis 2 , Tbt (=C6H2-2,4,6-{CH(SiMe3)2}3), and Ar* (=C6H3-2,6(C 6 H 2 -2,4,6-i-Pr 3 ) 2 )) at the M06-2X/Def2-TZVP,



REFERENCES

(1) Power, P. P. π-Bonding and the Lone Pair Effect in Multiple Bonds between Heavier Main Group Elements. Chem. Rev. 1999, 99, 3463−3504. 82

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega (2) Power, P. P. Silicon, Germanium, Tin and Lead Analogues of Acetylenes. Chem. Commun. 2003, 2091−2101. (3) Lein, M.; Krapp, A.; Frenking, G. Why Do the Heavy-Atom Analogues of Acetylene E2H2 (E = Si−Pb) Exhibit Unusual Structures? J. Am. Chem. Soc. 2005, 127, 6290−6299 and references therein.. (4) Sekiguchi, A.; Ichinohe, M.; Kinjo, R. The Chemistry of Disilyne with a Genuine Si−Si Triple Bond: Synthesis, Structure, and Reactivity. Bull. Chem. Soc. Jpn. 2006, 79, 825−832. (5) Power, P. P. Bonding and Reactivity of Heavier Group 14 Element Alkyne Analogues. Organometallics 2007, 26, 4362−4372. (6) Fischer, R. C.; Power, P. P. π-Bonding and the Lone Pair Effect in Multiple Bonds Involving Heavier Main Group Elements: Developments in the New Millennium. Chem. Rev. 2010, 110, 3877−3923 and references therein.. (7) Sasamori, T.; Han, J. S.; Hironaka, K.; Takagi, N.; Nagase, S.; Tokitoh, N. Synthesis and Structure of Stable 1,2-diaryldisilyne. Pure Appl. Chem. 2010, 82, 603−612. (8) Peng, Y.; Fischer, R. C.; Merrill, W. A.; Fischer, J.; Pu, L.; Ellis, B. D.; Fettinger, J. C.; Herber, R. H.; Power, P. P. Substituent Effects in Ditetrel Alkyne Analogues: Multiple vs. Single Bonded Isomers. Chem. Sci. 2010, 1, 461−468. (9) Sekiguchi, A.; Kinjo, R.; Ichinohe, M. A Stable Compound Containing A Silicon−Silicon Triple Bond. Science 2004, 305, 1755− 1757. (10) Wiberg, N.; Vasisht, S. K.; Fischer, G.; Mayer, P. Disilynes. III [1] a Relatively Stable Disilyne RSiSiR (R = SiMe(SitBu3)2). Z. Anorg. Allg. Chem. 2004, 630, 1823−1828. (11) Sasamori, T.; Hironaka, K.; Sugiyama, T.; Takagi, N.; Nagase, S.; Hosoi, Y.; Furukawa, Y.; Tokitoh, N. Synthesis and Reactions of a Stable 1,2-diaryl-1,2-dibromodisilene: A Precursor for Substituted Disilenes and 1,2-diaryldisilyne. J. Am. Chem. Soc. 2008, 130, 13856− 13857. (12) Stender, M.; Phillips, A. D.; Wright, R. J.; Power, P. P. Synthesis and Characterization of a Digermanium Analogue of an Alkyne. Angew. Chem., Int. Ed. 2002, 41, 1785−1787. (13) Stender, M.; Phillips, A. D.; Power, P. P. Formation of [Ar*Ge{CH2C(Me)C(Me)CH2}CH2C(Me)N]2 (Ar* = C6H3-2,6Trip2; Trip = C6H2-2,4,6-i-Pr3) Via Reaction of Ar*GeGeAr* with 2,3-dimethyl-1,3-butadiene: Evidence for the Existence of a Germanium Analogue of an Alkyne. Chem. Commun. 2002, 1312− 1313. (14) Pu, L.; Phillips, A. D.; Richards, A. F.; Stender, M.; Simons, R. S.; Olmstead, M. M.; Power, P. P. Germanium and Tin Analogues of Alkynes and Their Reduction Products. J. Am. Chem. Soc. 2003, 125, 11626−11636. (15) Sugiyama, Y.; Sasamori, T.; Hosoi, Y.; Furukawa, Y.; Takagi, N.; Nagase, S.; Tokitoh, N. Synthesis and Properties of a New Kinetically Stabilized Digermyne: New Insights for a Germanium Analogue of an Alkyne. J. Am. Chem. Soc. 2006, 128, 1023−1031. (16) Spikes, G. H.; Power, P. P. Lewis Base Induced Tuning of the Ge−Ge Bond Order in a “digermyne”. Chem. Commun. 2007, 85−87. (17) Phillips, A. D.; Wright, R. J.; Olmstead, M. M.; Power, P. P. Synthesis and Characterization of 2,6-Dipp2-H3C6SnSnC6H3-2,6Dipp2 (Dipp = C6H3-2,6-Pri2): A Tin Analogue of an Alkyne. J. Am. Chem. Soc. 2002, 124, 5930−5931. (18) Pu, L.; Twamley, B.; Power, P. P. Synthesis and Characterization of 2,6-Trip2H3C6PbPbC6H3-2,6-Trip2 (Trip = C6H2-2,4,6-i-Pr3): a Stable Heavier Group 14 Element Analogue of an Alkyne. J. Am. Chem. Soc. 2000, 122, 3524−3525. (19) Bino, A.; Ardon, M.; Shirman, E. Formation of a Carbon− Carbon Triple Bond by Coupling Reactions in Aqueous Solution. Science 2005, 308, 234−235. (20) Su, P.; Wu, J.; Gu, J.; Wu, W.; Shaik, S.; Hiberty, P. C. Bonding Conundrums in the C2 Molecule: A Valence Bond Study. J. Chem. Theory Comput. 2011, 7, 121−130. (21) Ploshnik, E.; Danovich, D.; Hiberty, P. C.; Shaik, S. The Nature of the Idealized Triple Bonds Between Principal Elements and the σ Origins of Trans-Bent GeometriesA Valence Bond Study. J. Chem. Theory Comput. 2011, 7, 955−968.

(22) Seidu, I.; Seth, M.; Ziegler, T. Role Played by Isopropyl Substituents in Stabilizing the Putative Triple Bond in Ar′EEAr′ [E = Si, Ge, Sn; Ar′ = C6H3-2,6-(C6H3-2,6-Pri2)2] and Ar*PbPbAr* [Ar* = C6H3-2,6-(C6H2-2,4,6-Pri3)2]. Inorg. Chem. 2013, 52, 8378−8388. (23) Danovich, D.; Bino, A.; Shaik, S. Formation of Carbon−Carbon Triply Bonded Molecules from Two Free Carbyne Radicals via a Conical Intersection. J. Phys. Chem. Lett. 2013, 4, 58−64. (24) Karni, M.; Apeloig, Y.; Schröder, D.; Zummack, W.; Rabezzana, R.; Schwarz, H. HCSiF and HCSiCl: The First Detection of Molecules with Formal CSi Triple Bonds. Angew. Chem., Int. Ed. 1999, 38, 331−335 and related references therein.. (25) Danovich, D.; Ogliaro, F.; Karni, M.; Apeloig, Y.; Cooper, D. L.; Shaik, S. Silynes (RCSiR′) and Disilynes (RSiSiR′): Why Are Less Bonds Worth Energetically More? Angew. Chem., Int. Ed. 2001, 40, 4023−4026. (26) Gau, D.; Kato, T.; Saffon-Merceron, N.; De Cózar, A.; Cossío, F. P.; Baceiredo, A. Synthesis and Structure of a Base-Stabilized CPhosphino-Si-AminoSilyne. Angew. Chem., Int. Ed. 2010, 49, 6585− 6588. (27) Lühmann, N.; Müller, T. A Compound with a Si−C Triple Bond. Angew. Chem., Int. Ed. 2010, 49, 10042−10044. (28) Liao, H.-Y.; Su, M.-D.; Chu, S.-Y. A Stable Species with a Formal GeC Triple Bond − A Theoretical Study. Chem. Phys. Lett. 2001, 341, 122−128. (29) For CGe, see: Wu, P.-C.; Su, M.-D. Theoretical Designs for Germaacetylene (RCGeR′): A New Target For Synthesis. Dalton Trans. 2011, 40, 4253−4259. (30) For CSn, see: Wu, P.-C.; Su, M.-D. Triply Bonded Stannaacetylene (RCSnR): Theoretical Designs and Characterization. Inorg. Chem. 2011, 50, 6814. (31) For CPb, see: Wu, P.-C.; Su, M.-D. A New Target for Synthesis of Triply Bonded Plumbacetylene (RCPbR): A Theoretical Design. Organometallics 2011, 30, 3293−3301. (32) Paetzold, P. Iminoboranes. Adv. Inorg. Chem. 1987, 31, 123− 170. (33) Paetzold, P. In Born Chemistry, Proceedings of the 6th International Meeting on Born Chemistry, Reactions at the BoronNitrogen Triple Bond; Hermanek, S.World Scientific: Singapore, 1987; pp 446−475. (34) (a) Paetzold, P. New Perspectives in Boron Nitrogen Chemistry − I. Pure Appl. Chem. 1991, 63, 345. (b) Paetzold, P. Boron−Nitrogen Analogues of Cyclobutadiene, Benzene and Cyclooctatetraene: Interconversions. Phosphorus Sulfur Silicon Relat. Elem. 1994, 93, 39−50. (35) Wright, R. J.; Phillips, A. D.; Allen, T. L.; Fink, W. H.; Power, P. P. Synthesis and Characterization of the Monomeric Imides Ar′MNAr“ (M = Ga or In; Ar′ or Ar” = Terphenyl Ligands) with Two-Coordinate Gallium and Indium. J. Am. Chem. Soc. 2003, 125, 1694−1695. (36) Su, J.; Li, X.-W.; Crittendon, R. C.; Robinson, G. H. How Short is a −Ga⋮Ga− Triple Bond? Synthesis and Molecular Structure of Na2[Mes*2C6H3-Ga⋮Ga-C6H3Mes*2] (Mes* = 2,4,6-i-Pr3C6H2): The First Gallyne. J. Am. Chem. Soc. 1997, 119, 5471−5472. (37) Robinson, G. H. Gallanes, Gallenes, Cyclogallenes, and Gallynes: Organometallic Chemistry about the Gallium−Gallium Bond. Acc. Chem. Res. 1999, 32, 773−782. (38) Xie, Y.; Grev, R. S.; Gu, J.; Schaefer, H. F., III; Schleyer, P. V. R.; Su, J.; Li, X.-W.; Robinson, G. H. The Nature of the Gallium−Gallium Triple Bond. J. Am. Chem. Soc. 1998, 120, 3773−3780. (39) Coates, G. E.; Livingstone, J. G. Aminoditolylborane and the Preparation of Diarylborinic Acids. J. Chem. Soc. 1961, 4909−4911. (40) However, a fully characterized monomeric phosphinoborane was not published until 1986. See: Feng, X.; Olmstead, M. N.; Power, P. P. The Structure of A Monomeric Phosphinoborane: Synthesis and X-Ray Crystal Structure of (Diphenylphosphino)dimesitylborane. Inorg. Chem. 1986, 25, 4615−4616. (41) Arif, A. M.; Cowley, A. H.; Pakulski, M.; Power, J. M. Diphospha- and diarsa-diboretanes. Four-membered Rings Containing 83

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega Boron and Phosphorus or Arsenic. J. Chem. Soc., Chem. Commun. 1986, 889−890. (42) For reviews concerning BP and BP bondings, see: Power, P. P. Boron−Phosphorus Compounds and Multiple Bonding. Angew. Chem., Int. Ed. Engl. 1990, 29, 449−460. (43) Power, P. P.; Moezzi, A.; Pestana, D. C.; Petrie, M. A.; Shoner, S. C.; Waggoner, K. M. Multiple Bonding, π-bonding Contributions and Aromatic Character in Isoelectronic Boron−Phosphorus, Boron− Arsenic, Aluminum−Nitrogen and Zinc−sulfur Compounds. Pure Appl. Chem. 1991, 63, 859−866. (44) Paine, R. T.; Noeth, H. Recent Advances in Phosphinoborane Chemistry. Chem. Rev. 1995, 95, 343−379. (45) Pestana, D. C.; Power, P. P. Nature of the Boron−Phosphorus Bond in Monomeric Phosphinoboranes and Related Compounds. J. Am. Chem. Soc. 1991, 113, 8426−8437. (46) Shih, T.-W.; Li, M.-C.; Su, M.-D. Doubly Bonded E13P and BE15 Molecules and Their Reactions with H2, Acetonitrile, Benzophenone, and 2,3-Dimethylbutadiene. Inorg. Chem. 2015, 54, 5154−5161. (47) Nyhlén, J.; Privalov, T. “Frustration” of Orbital Interactions in Lewis Base/Lewis Acid Adducts: A Computational Study of H2 Uptake by Phosphanylboranes R2PBR′2. Eur. J. Inorg. Chem. 2009, 2009, 2759−2764. (48) Karsch, H. H.; Hanika, G.; Huber, B.; Meindl, K.; König, S.; Krüger, C.; Müller, G. Monomeric Phosphinoboranes: the Role of the B- and P-Substituents for PB Multiple Bonding. J. Chem. Soc., Chem. Commun. 1989, 373−375. (49) Kubo, K.; Kawanaka, T.; Tomioka, M.; Mizuta, T. Synthesis and Crystal Structures of P-Iron-Substituted Phosphinoborane Monomers. Organometallics 2012, 31, 2026−2034. (50) Allen, T. L.; Scheiner, A. C.; Schaefer, H. F., III Theoretical Studies of Borylphosphine, Its Conjugate Base, and the Lithium Salt of Its Conjugate Base. The Use of Orbital Kinetic Energies to Determine the Origin of the Driving Force for Changes in Molecular Geometry. Inorg. Chem. 1990, 29, 1930−1936. (51) Gilbert, T. M.; Bachrach, S. M. Computational Studies of [2+2] and [4+2] Pericyclic Reactions between Phosphinoboranes and Alkenes. Steric and Electronic Effects in Identifying a Reactive Phosphinoborane that Should Avoid Dimerization. Organometallics 2007, 26, 2672−2678. (52) Privalov, T. On the Possibility of Conversion of Alcohols to Ketones and Aldehydes by Phosphinoboranes R2PBR′R″: A Computational Study. Chem. − Eur. J. 2009, 15, 1825−1829. (53) Kenward, A. L.; Piers, W. E. Heterolytic H2 Activation by Nonmetals. Angew. Chem., Int. Ed. 2008, 47, 38−41. (54) Coolidge, M. B.; Borden, W. T. Ab Initio Calculations on Borylphosphines: Prediction of a Synergistic Substituent Effect in Diborylphosphine. J. Am. Chem. Soc. 1990, 112, 1704−1706. (55) Geier, S. J.; Gilbert, T. M.; Stephan, D. W. Activation of H2 by Phosphinoboranes R2PB(C6F5)2. J. Am. Chem. Soc. 2008, 130, 12632− 12633. (56) Paetzold, P. The P-centered Chemistry of [(Di-tertbutylphosphanyl)imino]-(2,2,6,6-tetramethylpiperidino)borane. Adv. Inorg. Chem. 1987, 31, 123−170. (57) Grant, D. J.; Dixon, D. A. Thermodynamic Properties of Molecular Borane Phosphines, Alane Amines, and Phosphine Alanes and the [BH4−][PH4+], [AlH4−][NH4+], and [AlH4−][PH4+] Salts for Chemical Hydrogen Storage Systems from ab Initio Electronic Structure Theory. J. Phys. Chem. A 2005, 109, 10138−10147. (58) Grant, D. J.; Dixon, D. A. σ- and π-Bond Strengths in Main Group 3−5 Compounds. J. Phys. Chem. A 2006, 110, 12955−12962. (59) Ditty, M. J. T.; Power, W. P. 31P Chemical Shielding as a Function of Bond Angle: Phosphorus Chemical Shielding Surface of Phosphinoborane R2PBR′2. Can. J. Chem. 1999, 77, 1951−1961. (60) Jemmis, E. D.; Subramanian, G. Torsional Barriers in H2B:XCH3 (X = NH, O, PH, and S). J. Phys. Chem. 1994, 98, 8937−8939.

(61) Allen, T. L.; Fink, W. H. Theoretical Studies of Borylphosphine and Its Conjugate Base. 2. Internal Rotation and Inversion Transition States. Inorg. Chem. 1992, 31, 1703−1705. (62) Magnusson, E. Conformational Energies and Structural Relaxation in Singly Bonded Silicon, Phosphorus and SulfurCompounds. Aust. J. Chem. 1986, 39, 735−745. (63) Geier, S. J.; Gilbert, T. M.; Stephan, D. W. Synthesis and Reactivity of the Phosphinoboranes R2PB(C6F5)2. Inorg. Chem. 2011, 50, 336−344. (64) Stephan, D. W.; Erker, G. Frustrated Lewis Pairs: Metal-free Hydrogen Activation and More. Angew. Chem., Int. Ed. 2010, 49, 46− 76. (65) Schulz, S. Group 13/15 Organometallic CompoundsSynthesis, Structure, Reactivity and Potential Applications. Adv. Organomet. Chem. 2003, 49, 225−317. (66) Pons, V.; Baker, R. T. Soluble Boron−Nitrogen High Polymers from Metal-Complex-Catalyzed Amine Borane Dehydrogenation. Angew. Chem., Int. Ed. 2008, 47, 9600−9602. (67) Clark, T. J.; Lee, K.; Manners, I. Transition-Metal-Catalyzed Dehydrocoupling: A Convenient Route to Bonds between MainGroup Elements. Chem. − Eur. J. 2006, 12, 8634−8648. (68) Jaska, C. A.; Bartole-Scott, A.; Manners, I. Metal-catalyzed Routes to Rings, Chains and Macromolecules Based on Inorganic Elements. Dalton Trans. 2003, 4015−4021. (69) Maya, L. Boron Nitride Precursors? A Perspective. Appl. Organomet. Chem. 1996, 10, 175−182. (70) Manners, I. Polymers and the Periodic Table: Recent Developments in Inorganic Polymer Science. Angew. Chem., Int. Ed. Engl. 1996, 35, 1602−1621. (71) Whittell, G. R.; Manners, I. Synthesis of Ammonia Borane Nanoparticles and the Diammoniate of Diborane by Direct Combination of Diborane and Ammonia. Angew. Chem., Int. Ed. 2011, 50, 10288−10289. (72) Hügle, T.; Hartl, M.; Lentz, D. The Route to a Feasible Hydrogen-Storage Material: MOFs versus Ammonia Borane. Chem. − Eur. J. 2011, 17, 10184−10207. (73) Price, A. N.; Nichol, G. S.; Cowley, M. J. Phosphaborenes: Accessible Reagents for the Synthesis of C−C/P−B Isosteres. Angew. Chem., Int. Ed. 2017, 56, 9953−9957. (74) Kobayashi, K.; Nagase, S. Silicon−Silicon Triple Bonds: Do Substituents Make Disilynes Synthetically Accessible? Organometallics 1997, 16, 2489−2491. (75) Kobayashi, K.; Takagi, N.; Nagase, S. Do Bulky Aryl Groups Make Stable Silicon−Silicon Triple Bonds Synthetically Accessible? Organometallics 2001, 20, 234−236. (76) Takagi, N.; Nagase, S. Substituent Effects on Germanium− Germanium and Tin−Tin Triple Bonds. Organometallics 2001, 20, 5498−5500. (77) For instance, see: Chen, J.-Y.; Su, M.-D. Diboration of the E kE Double Bond by [2]Metallocenophanes (E = N, P, As, Sb, and Bi): A Theoretical Study. Organometallics 2010, 29, 5812−5820. (78) Bartlett, R. A.; Feng, X.; Power, P. P. Synthesis and Characterization of the Phosphinidene Borate Complexes [Li(Et2O)2PRB(Mes)2] and [Li(12-crown-4)2][RPB(Mes)2].THF [R = Ph, C6H11 or Mes (Mes = 2,4,6-Me3C6H2)]: The First Structurally Characterized Boron-Phosphorus Double Bonds. J. Am. Chem. Soc. 1986, 108, 6817−6819. (79) Greenwood, N. N.; Earnshaw, A. Chemistry of the Elements; Pergamon: Oxford, England, 1984; pp 452−513. (80) (a) The Wiberg bond index, which is used to screen atom pairs for the possible bonding in the natural bonding orbital (NBO) search, are performed with the NBO program. For details, see: http://www. chem.wisc.edu/~nbo5. (b) Also, see: Wiberg, K. B. Application of the Pople-Santry-Segal CNDO Method to the Cyclopropylcarbinyl and Cyclobutyl Cation and to Bicyclobutane. Tetrahedron 1968, 24, 1083− 1096. (81) Reed, A. E.; Curtiss, L. A.; Weinhold, F. Intermolecular Interactions from a Natural Bond Orbital, Donor−Acceptor Viewpoint. Chem. Rev. 1988, 88, 899−926. 84

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85

Article

ACS Omega (82) Liptrot, D. J.; Power, P. P. London Dispersion Forces in Sterically Crowded Inorganic and Organometallic Molecules. Nat. Rev. Chem. 2017, 1, No. 0004. (83) Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157−167. (84) Dapprich, S.; Frenking, G. Investigation of Donor−Acceptor Interactions: A Charge Decomposition Analysis Using Fragment Molecular Orbitals. J. Phys. Chem. 1995, 99, 9352−9362. (85) Glendening, E. D.; Weinhold, F. Natural resonance theory: I. general formalism. J. Comput. Chem. 1998, 19, 593−609. (86) Glendening, E. D.; Weinhold, F. Natural resonance theory: II. natural bond order and valency. J. Comput. Chem. 1998, 19, 610−627. (87) Glendening, E. D.; Badenhoop, J. K.; Weinhold, F. Natural resonance theory: III. Chemical Applications. J. Comput. Chem. 1998, 19, 628−646. (88) The electronegativities of boron and phosphorus are 2.051 and 2.253, respectively. See: Mann, J. B.; Meek, T. L.; Allen, L. C. Configuration Energies of the Main Group Elements. J. Am. Chem. Soc. 2000, 122, 2780−2783 and references cited therein..

85

DOI: 10.1021/acsomega.7b01480 ACS Omega 2018, 3, 76−85