Ab Initio and Electron Propagator Theory Study of Boron Hydrides

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Photoelectron spectra and electronic structure of boron diacetate formazanates

The electronic structure and cationic states of two 1,5-diphenylformazanes and two boron diacetate (B(OAc)₂) formazanates were modeled using the outer valence Green's function (OVGF) and density functional theory (DFT) methods. Comparison of data of the OVGF and ultraviolet photoelectron spectroscopy (UPS) methods made it possible to determine an effect of functional groups and complexing agents on energies of cationic states. Addition of NO₂-group at the γ-position of the chelate cycle causes stabilization of levels the five upper occupied molecular orbitals (MO) and destabilization of the bonding orbital π₃^Ph + π₃ level. The levels of MOs π₃^Ph-π₃ and n_- are stabilized due to influence of the complexing agent B(OAc)₂, with a difference in the shift of 0.67 eV. The ionization energies (I_n) changes for the π-orbitals of benzene rings are within the error of the OVGF method. Under methylation of phenyl groups, the differences between the calculated I_n, corresponding to the π-orbitals of aromatic substituents, are in good agreement with the experimental I_n shifts at transition from benzene to toluene. According to the OVGF method, in all the studied complexes the lowest unoccupied molecular orbital (LUMO) is localized mainly on the chelate cycle and has a strong acceptor character, which should contribute the low-lying charge-transfer electronic excitations. Moreover, an application of the DFT analog of the Koopmans' theorem with the BHHLYP and B2PLYP functionals made it possible to determine qualitatively a sequence of cationic states and energy intervals between them in the spectral range up to 10 eV. The DFT/wB97x/cc-pVTZ method data on the energy gap between the highest occupied molecular orbital (HOMO) and LUMO levels correlate with the OVGF/cc-pVTZ calculation results.

11B NMR Chemical Shift Predictions via Density Functional Theory and Gauge-Including Atomic Orbital Approach: Applications to Structural Elucidations of Boron-Containing Molecules

11B nuclear magnetic resonance (NMR) spectroscopy is a useful tool for studies of boron-containing compounds in terms of structural analysis and reaction kinetics monitoring. A computational protocol, which is aimed at an accurate prediction of 11B NMR chemical shifts via linear regression, was proposed based on the density functional theory and the gauge-including atomic orbital approach. Similar to the procedure used for carbon, hydrogen, and nitrogen chemical shift predictions, a database of boron-containing molecules was first compiled. Scaling factors for the linear regression between calculated isotropic shielding constants and experimental chemical shifts were then fitted using eight different levels of theory with both the solvation model based on density and conductor-like polarizable continuum model solvent models. The best method with the two solvent models yields a root-mean-square deviation of about 3.40 and 3.37 ppm, respectively. To explore the capabilities and potential limitations of the developed protocols, classical boron-hydrogen compounds and molecules with representative boron bonding environments were chosen as test cases, and the consistency between experimental values and theoretical predictions was demonstrated.

Convergence of nuclear magnetic shieldings and one-bond 1 J(11 B1 H) indirect spin-spin coupling constants in small boron molecules

Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH₃ molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual ¹¹ B and ¹ H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced ¹¹ B shielding in BH molecule.

B2(BO)6 0/- and B 2(BS) 6 0/- doubly bridged structures containing BO or BS as ligands

The investigation on the geometrical and electronic properties of B(2)(BO)(6) (0/-) and B(2)(BS)(6) (0/-) has been performed by density functional theory (DFT) using the B3LYP and BP86 methods. The chemical bonding in B(2)A(6) (A = H, BO, and BS) series is elucidated through the recently developed adaptive natural density partitioning (AdNDP). D(2h) B(2)(BO)(6) and B(2)(BS)(6) were found to possess two bridging η (2)-BO or η (2)-BS groups, as well as four terminal BO or BS groups that are analogs of diborane B(2)H(6). D(2)h B(2)(BO)(6) (-) and B(2)(BS)(6) (-) with two bridging η (2)-BO or η (2)-BS groups which are more stable than their corresponding D(3d) structures. The binding energy of B(2)(BO)(6) and B(2)(BS)(6) with respect to B(2)(BO)(6) (D2h) → 2B(BO)(3) (D(3h)) and B(2)(BS)(6)(D(2h)) → 2B(BS)(3) (D(3h)) are estimated to be (△)E = 19.8 and 40.6 kcal mol(-1) at CCSD(T)//B3LYP level, respectively. This finding advances the boronyl chemistry and helps establish the isolobal analogy between boron-rich oxide clusters and boranes.

Computational study of the initial stage of diborane pyrolysis

The rate constants for the association of two boranes to form diborane are investigated using several methods. The most sophisticated method is the variable reaction coordinate-variational transition state theory (VRC-VTST) which has been developed to handle reactions with no enthalpic barriers. The calculated rate constant of 8.2 × 10(-11) cm(3)·molecule(-1)·s(-1) at 545 K is in good agreement with experiment. The rate constant was also computed using conventional VTST with the G4 composite method. Two variations of the multistep mechanisms for diborane pyrolysis are presented. One is initiated by the step B2H6 ⇄ 2 BH3 while the other begins with 2 B2H6 ⇄ B3H9 + BH3 as the initial elementary step. Both variations are 3/2 order in diborane and have the same activation energy (G4, 28.65 kcal/mol at 420 K). In contrast, the traditional mechanism involving a B3H9 intermediate with C3v symmetry has a higher activation energy (33.37 kcal/mol). The two variations involve a C2-symmetry penta-coordinate B3H9 structure that, while an electronic minimum, is not a stationary point on the free energy path between B2H6 + BH3 → B3H7 + H2. While the calculated activation barrier is higher than the recently determined experimental barrier, the variation in reported values is large (22.0-29.0 kcal/mol). We discuss possible sources of disagreement between experiment and theory.

Heats of Formation of Boron Hydride Anions and Dianions and Their Ammonium Salts [BnHmy-][NH4+]y with y = 1−2

Thermochemical parameters of the closo boron hydride BnHn2- dianions, with n=5-12, the B3H8- and B11H14- anions, and the B5H9 and B10H14 neutral species were predicted by high-level ab initio electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T)/complete basis set (CBS) extrapolated energies, plus additional corrections were used to predict the heats of formation of the simplest BnHmy- species in the gas phase in kcal/mol at 298 K: DeltaHf(B3H8-)=-23.1+/-1.0; DeltaHf(B5H52-)=119.4+/-1.5; DeltaHf(B6H62-)=64.1+/-1.5; and DeltaHf(B5H9)=24.1+/-1.5. The heats of formation of the larger species were evaluated by the G3 method from hydrogenation reactions (values at 298 K, in kcal/mol with estimated error bars of+/-3 kcal/mol): DeltaHf(B7H72-)=51.8; DeltaHf(B8H82-)=46.1; DeltaHf(B9H92-)=24.4; DeltaHf(B10H102-)=-12.5; DeltaHf(B11H112-)=-11.8; DeltaHf(B12H122-)=-86.3; DeltaHf(B11H14-)=-57.3; and DeltaHf(B10H14)=18.7. A linear correlation between atomization energies of the dianions and energies of the BH units was found. The heats of formation of the ammonium salts of the anions and dianions were predicted using lattice energies (UL) calculated from an empirical expression based on ionic volumes. The UL values (0 K) of the BnHn2- dianions range from 319 to 372 kcal/mol. The values of UL for the B3H8- and B11H14- anions are 113 and 135 kcal/mol, respectively. The calculated lattice energies and gas-phase heats of formation of the constituent ions were used to predict the heats of formation of the ammonium crystal salts [BnHmy-][NH4+]y. These results were used to evaluate the thermodynamics of the H2 release reactions from the ammonium hydro-borate salts.

Conformation Effects on the Electronic Structures of β-Alanine

Ten low-lying conformers of beta-alanine have been studied by the hybrid density functional B3LYP/aug-cc-pVDZ method. Energetic extrapolation calculations at the MP3 and MP4(SDQ) levels of theory and the theoretical photoelectron spectra simulated with the electron propagation theory demonstrate that there are at least three gauche conformers (G1, G2, and G3) in gas-phase experiments. The calculated ionization potentials are in good agreement with the experimental data available in the literature. Natural bond orbital and atoms-in-molecules analyses exhibit a remarkable influence on the molecular electronic structures by the strong intramolecular hydrogen bonding O-H...N in the neutral conformer G2. Remarkable internal rotations of the COOH group are found in the cationic G1+ and G3+ with respect to the neutral conformers. A distonic [NH3+-(CH2)2-COO*] radical can be formed through the spontaneous intramolecular proton transfer in G2+. A novel intramolecular hydrogen bonding, C-H...O, is found in the anti A1+ cation.

Novel diborane-analogue transition structures for borane reactions with alkyl halides

Ab initio and DFT calculations were performed to examine the mechanisms of reduction of alkyl halides and formaldehyde by borane. With alkyl halides, the optimized transition structure geometry resembled diborane, with a pair of hydrogen atoms bridging the boron and carbon atoms by three-center-two-electron bonds. A similar transition structure was found for the reduction of formaldehyde, although it was not the lowest-energy transition structure. Solvation by dimethyl ether or dimethyl sulfide disrupted this bridging with chloromethane, while both ligands dissociated from borane during the reduction of formaldehyde. The high calculated activation free energies of alkyl halide reduction are consistent with their observed lack of reactivity with borane.

Theoretical Study of Isoelectronic Molecules: B6H10, 2-CB5H9, 2,3-C2B4H8, 2,3,4-C3B3H7, and 2,3,4,5-C4B2H6

Isoelectronic molecules regarding B6H10, 2-CB5H9, 2,3-C2B4H8, 2,3,4-C3B3H7, and 2,3,4,5-C4B2H6 are studied by the density functional B3LYP/6-311G(d,p) method and the electron propagator theory in the partial third-order quasiparticale approximation, as well as the extrapolated calculation with the coupled-cluster CCSD(T) theory. The calculated ionization potentials are in good agreement with the experimental data from photoelectron spectroscopy. Valence structures are characterized with natural orbital bond (NBO) theory, exhibiting the multiple three-center two-electron bonds B-H-B, B-B-B, C-B-B, B-C-B, and C-B-C, and chemical bond rearrangements in the cations.