Probing Dative and Dihydrogen Bonding in ... - ACS Publications

Jul 11, 2017 - assigned to the BN stretch and BH3 rock, respectively.7 This assignment ... crystal structure model containing nine BH3NH3 subunits rep...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCA

Probing Dative and Dihydrogen Bonding in Ammonia Borane with Electronic Structure Computations and Raman under Nitrogen Spectroscopy Katelyn M. Dreux,† Louis E. McNamara,† John T. Kelly,†,‡ Ashley M. Wright,† Nathan I. Hammer,† and Gregory S. Tschumper*,† †

Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States S Supporting Information *

ABSTRACT: Although ammonia borane is isoelectronic with ethane and they have similar structures, BH3NH3 exhibits rather atypical bonding compared to that in CH3CH3. The central bond in ammonia borane is actually a coordinate covalent or dative bond rather than the conventional covalent C−C bond in ethane where each atom donates one electron. In addition, strong intermolecular dihydrogen bonds can form between two or more ammonia borane molecules compared to the relatively weak dispersion forces between ethane molecules. As a result, ammonia borane’s physical properties are very sensitive to the environment. For example, gas-phase and solid-state ammonia borane have very different BN bond lengths and BN stretching frequencies, which led to much debate in the literature. It has been demonstrated that the use of cluster models based on experimental crystal structures led to better agreement between theory and experiment. Here, we employ a variety of cluster models to track how the interaction energies, bond lengths, and vibrational normal modes evolve with the size and structural characteristics of the clusters. The M06-2X/6-311++G(2df,2pd) level of theory was selected for this analysis on the basis of favorable comparison with CCSD(T)/aug-cc-pVTZ data for the ammonia borane monomer and dimer. Fourteen unique fully optimized molecular cluster geometries, (BH3NH3)n≤12, and nine crystal models, (BH3NH3)n≤19, were used to elucidate how the local environment impacts ammonia borane’s physical properties. Computational results for the BN stretching frequencies are also compared directly to the Raman spectrum of solid ammonia borane at 77 K using Raman under liquid nitrogen spectroscopy (RUNS). A strong linear correlation was found to exist between the BN bond length and stretching frequency, from an isolated monomer to the most distorted BH3NH3 unit in a cluster or crystal structure model. Excellent agreement was seen between the frequencies computed for the largest crystal model and the RUNS experimental spectra (typically within a few wavenumbers).

1. INTRODUCTION Since it was first synthesized in 19551 ammonia borane has generated a large amount of interest both experimentally2−20 and computationally.21−47 The coordinate covalent, often called dative, bond between boron and nitrogen is an intramolecular bond with an interaction energy of −43.91 kcal mol−1 at the estimated CCSD(T) complete basis set (CBS) limit that arises from nitrogen donating a pair of electrons into the empty porbital on boron.46 The presence of both hydridic and protonic hydrogens in BH3NH3 allows for the formation of dihydrogen bonds (BHδ−···Hδ+N) between adjacent molecules. These interactions give rise to short intermolecular H−H distances in the crystal structure ( 0.98) between the two properties regardless of the which data set is considered. Figure 5b shows the closely related correlation between changes in BN bond lengths (ΔRBN) and BN stretching frequency shifts (ΔωBN) upon complexation for all the isolated clusters as well as the optimized BH3NH3 unit of the crystal structures in Figure 3. The equations in Figure 5b have intercepts approaching 0 which merely indicates that there is no frequency shift if the bond length does not change. Figure 5c shows the increase in electronic binding energy as the number of monomers increases in the isolated clusters. The

Table 5. Low Energy Intermolecular Lattice Modes (cm−1) of Solid Ammonia Borane

9-Cs

ref 18

ref 14

ref 12

this work Raman 77 K

16C1

19C1

IR 10 K

Raman 15 K

INS 30 K

359 337 319 256 235 209 175

337 337

399 334

221 212

202

210

185

171

156

141 128

159 150 116

150

87

100

94

BH3 twist

371

397 390

352

NH3 twist

198

236 219

193

BH3 rock

162

170

NH3 rock

152

Rotation

127

Translation

111

Translation

78

170 156 147 138 136 131 115 113 90 75

158 117 116 104

G

DOI: 10.1021/acs.jpca.7b03509 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Figure 5. Plots of the correlation between various properties of ammonia borane clusters and crystal models. The table inserts for (a) and (b) contain the equations and R2 values for sets of data, including the isolated clusters (C), the crystal models (X), and the BH3NH3 monomer (M). The table inserts for (c) contain data from the isolated clusters, divided into three classes symmetry equivalent (Symm Equiv), centrally coordinated (Centrally Coord), and all isolated clusters (All). See section 3.6 for additional details.

ammonia borane and provides a more resolved spectrum, particularly in the BH and NH stretching regions. A rigorous comparison of methods and basis sets on the monomer and the C2h dimer was carried out (Tables S1 and S2 in the Supporting Information), which determined that the M06-2X/6-311+ +G(2df,2pd) level of theory provides good agreement between computational demand and an accurate description of the unusual inter- and intramolecular interactions of (BH3NH3)n. Although previous Hartree−Fock computations revealed that the BN stretching frequency in ammonia borane is exquisitely sensitive to its environment,32 the computations presented here reveal that Hartree−Fock dramatically overestimates the effect. This systematic analysis of both isolated clusters and crystal models of (BH3NH3)n with the M06-2X DFT method has demonstrated that system size (i.e., the magnitude of n) is not the main factor influencing the BN stretching frequency. Instead, the relative orientations of the BH3NH3 fragments plays a large role in determining the magnitude of the BN stretching frequency. For example, the central BH3NH3 in the

least-squares regression for different categories of isolated clusters are displayed in Figure 5c where “Symm Equiv” denotes data associated with the symmetry equivalent clusters, “Centrally Coord” denotes data associated with the centrally coordinated clusters, and “All” denotes the data for all the isolated clusters. The symmetry equivalent clusters grow in energy slightly faster than their centrally coordinated counterparts. As expected, the intercepts of the linear equations in Figure 5c approach 0 as n → 1 because the electronic binding energies is effectively 0 for a single molecule. There is no clear linear correlation between the BN frequency shifts and the number of fragments (n), the electronic binding energies or the normalized electronic binding energies of the isolated clusters. (See figures in the Supporting Information.)

4. CONCLUSIONS Raman under liquid nitrogen spectroscopy (RUNS) shows good agreement with previous low-temperature spectroscopic studies performed on the orthorhombic polymorph of H

DOI: 10.1021/acs.jpca.7b03509 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

(2) Hughes, E. W. The Crystal Structure of Ammonia-borane, H3NBH3. J. Am. Chem. Soc. 1956, 78, 502−503. (3) Lippert, E. L.; Lipscomb, W. N. The Structure of H3NBH3. J. Am. Chem. Soc. 1956, 78, 503−504. (4) Taylor, R. C.; Cluff, C. L. Vibrational Frequency Associated with the Boron-Nitrogen Dative Bond in Amine Boranes. Nature (London, U. K.) 1958, 182, 390−391. (5) Taylor, R. C. Vibrational Frequencies, Assignments, and Force Constants for some Compounds Containing Boron-Nitrogen Dative Bonds. Adv. Chem. Ser. 1964, 42, 59−70. (6) Sawodny, W.; Goubeau, J. Oscillation Spectra and Force Constants of Some X3Z←NR3-type Compounds. Z. Phys. Chem. (Muenchen, Ger.) 1965, 44, 227−241. (7) Smith, J.; Seshadri, K. S.; White, D. Infrared Spectra of Matrix Isolated Borane-Ammonia, Perdeutero Borane-Ammonia, BoraneAmmonia-D3. J. Mol. Spectrosc. 1973, 45, 327−337. (8) Suenram, R. D.; Thorne, L. R. Microwave Spectrum and Dipole Moment of Borane Monoammoniate. Chem. Phys. Lett. 1981, 78, 157−160. (9) Thorne, L. R.; Suenram, R. D.; Lovas, F. J. Microwave Spectrum, Torsional Barrier, and Structure of Borane Monoammoniate. J. Chem. Phys. 1983, 78, 167−171. (10) Carpenter, J. D.; Ault, B. S. Matrix Isolation Study of the Mechanism of the Reaction of Diborane with Ammonia: Pyrolysis of the Borane·Ammonia (H3B·NH3) Adduct. Chem. Phys. Lett. 1992, 197, 171−174. (11) Klooster, W. T.; Koetzle, T. F.; Siegbahn, P. E. M.; Richardson, T. B.; Crabtree, R. H. Study of the N-H···H-B Dihydrogen Bond Including the Crystal Structure of BH3NH3 by Neutron Diffraction. J. Am. Chem. Soc. 1999, 121, 6337−6343. (12) Allis, D. G.; Kosmowski, M. E.; Hudson, B. S. The Inelastic Neutron Scattering Spectrum of H3B:NH3 and the Reproduction of Its Solid-State Features by Periodic DFT. J. Am. Chem. Soc. 2004, 126, 7756−7757. (13) Hess, N. J.; Bowden, M. E.; Parvanov, V. M.; Mundy, C.; Kathmann, S. M.; Schenter, G. K.; Autrey, T. Spectroscopic Studies of the Phase Transition in Ammonia Borane: Raman Spectroscopy of Single Crystal NH3BH3 as a Function of Temperature from 88 to 330 K. J. Chem. Phys. 2008, 128, 034508. (14) Ziparo, C.; Colognesi, D.; Giannasi, A.; Zoppi, M. Raman Spectra of Ammonia Borane: Low Frequency Lattice Modes. J. Phys. Chem. A 2012, 116, 8827−8832. (15) Liu, A.; Song, Y. In Situ High-Pressure and Low-Temperature Study of Ammonia Borane by Raman Spectroscopy. J. Phys. Chem. C 2012, 116, 2123−2131. (16) Lin, Y.; Ma, H.; Matthews, C. W.; Kolb, B.; Sinogeikin, S.; Thonhauser, T.; Mao, W. L. Experimental and Theoretical Studies on a High Pressure Monoclinic Phase of Ammonia Borane. J. Phys. Chem. C 2012, 116, 2172−2178. (17) Kathmann, S. M.; Mundy, C. J.; Schenter, G. K.; Autrey, T.; Aeberhard, P. C.; David, B.; Jones, M. O.; Ramirez-Cuesta, T. Understanding Vibrational Anharmonicity and Phonon Dispersion in Solid Ammonia Borane. J. Phys. Chem. C 2012, 116, 5926−5931. (18) Paolone, A.; Teocoli, F.; Sanna, S.; Palumbo, O.; Autrey, T. Temperature Dependence of the Infrared Spectrum of Ammonia Borane: Librations, Rotations, and Molecular Vibrations. J. Phys. Chem. C 2013, 117, 729−734. (19) Nylén, J.; Eriksson, L.; Bensen, D.; Häussermann, U. Characterizations of a High Pressure, High Temperature Modification of Ammonia Borane (BH3NH3). J. Chem. Phys. 2013, 139, 054507. (20) Yao, Y.; Yong, X.; Tse, J. S.; Greschner, M. J. Dihydrogen Bonding in Compressed Ammonia Borane and Its Roles in Structural Stability. J. Phys. Chem. C 2014, 118, 29591−29598. (21) Umeyama, H.; Morokuma, K. Molecular Orbital Studies of Electron Donor-acceptor Complexes. 3. Energy and Charge Decomposition Analyses for Several Strong Complexes: Carbon Monoxide-Borane, Ammonia-Borane, Methylamine-Borane, Trimethylamine-Borane, and Ammonia-Boron Trifluoride. J. Am. Chem. Soc. 1976, 98, 7208−7220.

small isolated cluster 3-Cs where the fragments adopt an antiparallel orientation experiences the same +109 cm−1 shift relative to the isolated monomer as the central fragment in the largest crystal model (19-C1) where the dipole moments of the fragments are more closely aligned. When the antiparallel motif in the isolated clusters is extended, the shift can actually grow as large as +176 cm−1 in 12-C3v. The sensitivity of the BN stretching frequency to its environment can also been observed within a given centrally coordinated isolated cluster of (BH3NH3)n where the BN stretch of the central fragment is consistently at least 100 cm−1 higher than that of the peripheral BH3NH3 fragments. In contrast, the surrounding environment does not appear to affect the crystal models as strongly as the isolated clusters. For example, the BN stretching frequency associated with crystal model 9-Cs changes by less than 10 cm−1 for crystal model 16-C1 in which two 9-Cs models are effectively merged to allow two BH3NH3 subunits to relax rather than just one. The 19-C1 crystal model has a BN stretching frequency in remarkably good agreement (within 2 cm−1) with the experimental RUNS data. For the intermolecular and lattice modes of solid ammonia borane, 16-C1 presents the best agreement with experiment likely due to the optimization of two central fragments rather than just one.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpca.7b03509. Monomer and 2-C2h dimer calibration information performed in this study (Tables S1 and S2), leastsquares regression plots of several properties (Figure S1a−S1c), optimized geometries, Cartesian coordinates, harmonic vibrational frequencies, IR intensities, and Raman activities are available for all the isolated clusters and crystal models (PDF)



AUTHOR INFORMATION

Corresponding Author

*G. S. Tschumper. E-mail: [email protected]. Phone: +1 662 915 7301. Fax: +1 662 915 7300. ORCID

Nathan I. Hammer: 0000-0002-6221-2709 Gregory S. Tschumper: 0000-0002-3933-2200 Present Address ‡

Wilhelm-Ostwald-Institut für Physikalische and Theoretische Chemie Universität Leipzig, Linnéstrasse 2, 04103, Leipzig (Germany). Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported in part by the National Science Foundation (CHE-0955550, CHE-1338056, OIA-1430364, OIA-1539035, CHE-1532079). The Mississippi Center for Supercomputing Research (MCSR) is also thanked for a generous allocation of time on their computational resources.



REFERENCES

(1) Shore, S. G.; Parry, R. W. The Crystalline Compound Ammoniaborane, H3NBH3. J. Am. Chem. Soc. 1955, 77, 6084−6085. I

DOI: 10.1021/acs.jpca.7b03509 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A (22) Haaland, A. Covalent Versus Dative Bonds to Main Group Metals, a Useful Distinction. Angew. Chem., Int. Ed. Engl. 1989, 28, 992−1007. (23) Vijay, A.; Sathyanarayana, D. N. SCF and Electron Correlation Studies on Structure, Force Constants and Vibrational Spectra of Borane Monoammoniate Complex. Chem. Phys. 1995, 198, 345−352. (24) Jonas, V.; Frenking, G.; Reetz, M. T. Comparative Theoretical Study of Lewis Acid-Base Complexes of BH3, BF3, BCl3, AlCl3, and SO2. J. Am. Chem. Soc. 1994, 116, 8741−8753. (25) Richardson, T.; de Gala, S.; Crabtree, R. H.; Siegbahn, P. E. M. Unconventional Hydrogen Bonds: Intermolecular B-H···H-N Interactions. J. Am. Chem. Soc. 1995, 117, 12875−12876. (26) Popelier, P. L. A. Characterization of a Dihydrogen Bond on the Basis of the Electron Density. J. Phys. Chem. A 1998, 102, 1873−1878. (27) Merino, G.; Bakhmutov, V. I.; Vela, A. Do Cooperative ProtonHydride Interactions Explain the Gas-Solid Structural Difference of BH3NH3? J. Phys. Chem. A 2002, 106, 8491−8494. (28) Cybulski, H.; Sadlej, J. The Properties of the DihydrogenBonded Dimer (BH3NH3)2. Polish J. Chem. 2007, 81, 683−697. (29) Cramer, C. J.; Gladfelter, W. L. Ab Initio Characterization of [H3N·BH3]2, [H3N·AlH3]2, and [H3N·GaH3]2. Inorg. Chem. 1997, 36, 5358−5362. (30) Li, J.; Zhao, F.; Jing, F. B-Hδ−σ Bond as Dihydrogen Bond Acceptor: Some Theoretical Observations and Predictions. J. Chem. Phys. 2002, 116, 25−32. (31) Kar, T.; Scheiner, S. Comparison between Hydrogen and Dihydrogen Bonds among H3BNH3, H2BNH2, and NH3. J. Chem. Phys. 2003, 119, 1473−1482. (32) Dillen, J.; Verhoeven, P. The End of a 30-Year-Old Controversy? A Computational Study of the B-N Stretching Frequency of BH3-NH3 in the Solid State. J. Phys. Chem. A 2003, 107, 2570−2577. (33) Morrison, C. A.; Siddick, M. M. Density Functional Calculations: Dihydrogen Bonds in Solid BH3NH3. Angew. Chem., Int. Ed. 2004, 43, 4780−4782. (34) Meng, Y.; Zhou, Z.; Duan, C.; Wang, B.; Zhong, Q. Nonconventional Hydrogen Bonding Interaction of BH3NH3 Complexes: A Comparative Theoretical Study. J. Mol. Struct.: THEOCHEM 2005, 713, 135−144. (35) Demaison, J.; Lievin, J.; Csar, A. G.; Gutle, C. Equilibrium Structure and Torsional Barrier of BH3NH3. J. Phys. Chem. A 2008, 112, 4477−4482. (36) Mitoraj, M. P. Bonding in Ammonia Borane: An Analysis Based on the Natural Orbitals for Chemical Valence and the Extended Transition State Method (ETS-NOCV). J. Phys. Chem. A 2011, 115, 14708−14716. (37) Guerra, D.; David, J.; Restrepo, A. (H3N-BH3)4: The Ammonia Borane Tetramer. Phys. Chem. Chem. Phys. 2012, 14, 14892−14897. (38) Sams, R. L.; Xantheas, S. S.; Blake, T. A. Vapor Phase Infrared Spectroscopy and Ab Initio Fundamental Anharmonic Frequencies of Ammonia Borane. J. Phys. Chem. A 2012, 116, 3124−3136. (39) Parafiniuk, M.; Mitoraj, M. P. On the Origin of Internal Rotation in Ammonia Borane. J. Mol. Model. 2014, 20, 2−9. (40) Zapata-Escobar, A. D.; Cárcamo-Camacho, T.; Hadad, C. Z.; Restrepo, A. On the Nature of the Trimer, Tetramer, and Pentamer of Ammonia Borane. Theor. Chem. Acc. 2016, 135, 95. (41) Yuan, B.; Shin, J.-W.; Bernstein, E. R. Dynamics and Fragmentation of van der Waals and Hydrogen Bonded Cluster Cations: (NH3)n and (BH3NH3)n Ionized at 10.51 eV. J. Chem. Phys. 2016, 144, 144315. (42) Huang, F.; Jiang, J.; Wen, M.; Wang, Z.-X. Assessing the Performance of Commonly Used DFT Functionals in Studying the Chemisty of Frustrated Lewis Pairs. J. Theor. Comput. Chem. 2014, 13, 1350074−1350090. (43) Nguyen, V. S.; Matus, M. H.; Nguyen, M. T.; Dixon, D. A. Reactions of Diborane with Ammonia and Ammonia Borane: Catalytic Effects for Multiple Pathways for Hydrogen Release. J. Phys. Chem. A 2008, 112, 9946−9954.

(44) Matus, M. H.; Grant, D. J.; Nguyen, M. T.; Dixon, D. A. Fundamental Thermochemical Properties of Ammonia Borane and Dehydrogenated Derivatives (BNHn, n = 0−6). J. Phys. Chem. C 2009, 113, 16553−16560. (45) Nguyen, V. S.; Matus, M. H.; Grant, D. J.; Nguyen, M. T.; Dixon, D. A. Computational Study of the Release of N2 from Ammonia Borane Dimer (BH3NH3)2 and Its Ion Pair Isomers. J. Phys. Chem. A 2007, 111, 8844−8856. (46) Karthikeyan, S.; Sedlak, R.; Hobza, P. on the Nature of Stabilization in Weak, Medium, and Strong Changre-Transfer Complexes: CCSD(T)/CBS and SAPT Calculations. J. Phys. Chem. A 2011, 115, 9422−9428. (47) Vijayalakshmi, K. P.; Suresh, C. H. Ammonia Borane Clusters: Energetics of Dihydrogen Bonding, Cooperativity, and Role of Electrostatics. J. Phys. Chem. A 2017, 121, 2704−2714. (48) Compton, R. N.; Hammer, N. I. Raman Under Liquid Nitrogen (RUN). J. Phys.: Conf. Ser. 2014, 548, 012017. (49) Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101. (50) Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157−167. (51) Becke, A. D. Density-functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (52) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into Functional of the Electron Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (53) Møller, C.; Plesset, M. S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46, 618−622. (54) Purvis, G.; Bartlett, R. A Full Coupled-Cluster Singles and Doubles Model: The Inclusion of Disconnected Triples. J. Chem. Phys. 1982, 76, 1910−1918. (55) Dunning, T. H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. the Atoms Boron Through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (56) Frisch, M. J.; Pople, J. A.; Binkly, J. S. Self-Consistent Molecular Orbital Methods 25. Supplementary Functions for Gaussian Basis Sets. J. Chem. Phys. 1984, 80, 3265−3269. (57) Frisch, M. J.; Trucks, G.; Schlegel, H.; Scuseria, G.; Robb, M.; Cheeseman, J.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.; et al. Gaussian 09 revision D.01; Gaussian Inc.: Wallingford, CT, 2009. (58) Stanton, J. F.; Gauss, J.; Harding, M. E.; Szalay, P. G.; Auer, A. A.; Bartlett, R. J.; Benedikt, U.; Berger, C.; Bernholdt, D. E.; Bomble, Y. J.; et al. CFOUR, Coupled-Cluster Techniques for Computational Chemistry and the integral packages MOLECULE (J. Almlöf and P. R. Taylor), PROPS (P. R. Taylor), and ABACUS (T. Helgaker, H. J. Aa. Jensen, P. Jørgensen, and J. Olsen) and ECP routines by A. V. Mitin and C. van Wüllen. For the current version, see http://www.cfour.de (last accessed June 26, 2017). (59) Kesharwani, M. K.; Brauer, B.; Martin, J. M. L. Frequency and Zero-Point Vibrational Energy Scale Factors for Double-Hybrid Density Functionals (and Other Selected Methods): Can Anharmonic Force Fields Be Avoided? J. Phys. Chem. A 2015, 119, 1701−1714.

J

DOI: 10.1021/acs.jpca.7b03509 J. Phys. Chem. A XXXX, XXX, XXX−XXX