Search NASA⌕ Search

SEARCH · Search NASA

Results for “isomer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theoretical Study of the Reaction of the Methylidyne Radical (CH; X 2 Π) with 1-Butyne (CH 3 CH 2 CCH; X 1 A')

In this study, ab initio CCSD(T)-F12/cc-pVTZ-f12//ωB97X-D/6-311G(d,p) + ZPE(ωB97X-D/6-311G(d,p)) calculations were carried out to unravel the area of the C 5 H 7 potential energy surface accessed by the reaction of the methylidyne radical with 1-butyne. The results were utilized in Rice-Ramsperger-Kassel-Marcus calculations of the product branching ratios at the zero pressure limit. The preferable reaction mechanism has been shown to involve (nearly) instantaneous decomposition of the initial reaction adducts, whose structures are controlled by the isomeric form of the C 4 H 6 reactant. If CH adds to the triple CΞC bond in the entrance reaction channel, the reaction is predicted to predominantly form the methylenecyclopropene + methyl (CH 3 ) and cyclopropenylidene + ethyl (C 2 H 5 ) products roughly in a 2:1 ratio. CH insertion into a C-H bond in the methyl group of 1-butyne is anticipated to preferentially form ethylene + propargyl (C 3 H 3 ) by the C-C bond β-scission in the initial complex, whereas CH insertion into C-H of the CH 2 group would predominantly produce vinylacetylene + methyl (CH 3 ) also by the C-C bond β-scission in the adduct. The barrierless and highly exoergic CH + 1-butyne reaction, facile in cold molecular clouds, is not likely to lead to the carbon skeleton molecular growth but generates C 4 H 4 isomers methylenecyclopropene, vinylacetylene, and 1,2,3-butatriene and smaller C 2 and C 3 hydrocarbons such as methyl, ethyl, and propargyl radicals, ethylene, and cyclopropenylidene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Water Network Shape-Dependence of Local Interactions with the Microhydrated –NO 2 – and –CO 2 – Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH 3 NO 2 – radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH 3 NO 2 – ·(H 2 O) n=1–6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO 2 – ·(H 2 O) n=5–7 (R = CD 3 CD 2 ), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a “book” isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor–acceptor–donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient –XO 2 H–OH–(H 2 O) 2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH 3 NO 2 – ·(H 2 O) n=1–6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. Furthermore, this behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Substituent Effects on the Electronic Spectroscopy of Four-Carbon Criegee Intermediates

Atmospheric ozonolysis of biogenic and anthropogenic alkenes generates zwitterionic carbonyl oxide intermediates (R 1 R 2 C=O + O – ), known as Criegee intermediates, with different structural motifs and conformations. This study reports a systematic laboratory study of substituent effects on the electronic spectroscopy of four-carbon Criegee intermediates (CIs) with methyl–ethyl (MECI) and isopropyl (IPCI) groups, which are isomers produced in ozonolysis of asymmetric branched alkenes. The four-carbon CIs are separately generated by an alternative synthetic route, and spectroscopically characterized on the strong π* ← π transition associated with the carbonyl oxide group in a pulsed supersonic expansion with VUV photoionization at 118 nm and UV-induced depletion of the m/z 88 signal. The resultant broad and unstructured UV spectral features for MECI and IPCI are peaked at ca. 320 and 330 nm, respectively, with large absorption cross-sections of ca. 10 –17 cm 2 . Comparisons are made with the four-carbon CIs formed in isoprene ozonolysis, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide), which have the same backbone connectivity as MECI and IPCI but have extended conjugation across the vinyl and carbonyl groups. A remarkable 50 nm shift of the peak absorption to longer wavelength is observed for MVK-oxide and MACR-oxide compared to MECI and IPCI, respectively. Furthermore, vertical excitation energies computed theoretically agree well with the experimental findings, confirming that the spectral shifts are caused by the extended π conjugation in the isoprene-derived Criegee intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. In this study, we observe the excited state dynamics following excitation with 266 nm UV light to the 1 ππ* (S 2 ) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the 1 nπ* (S 1 ) and 3 ππ* (T 1 ) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576—16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Electronic Structure and Anion Photoelectron Spectroscopy of Uranium–Gold Clusters UAu n – , n = 3–7

A collaborative effort between experiment and theory towards elucidating the electronic and molecular structures of uranium-gold clusters is presented. Anion photoelectron spectra of UAu n – (n=3-7) were taken at third (355 nm) and fourth (266 nm) harmonics of a Nd:YAG laser, as well as excimer (ArF 193nm) photon energies, where the experimental adiabatic electron affinities (AEA) and vertical detachment energies (VDE) values were measured. Complementary first principles calculations were subsequently carried out to corroborate experimentally determined electron detachment energies and to determine the geometry and electronic structure for each cluster. Except for the ring-like neutral isomer of UAu 6 where one unpaired electron is spread over the Au atoms, all other neutral and anionic UAu n clusters (n = 3 - 7) were calculated to possess open-shell electrons with the unpaired electrons localized on the central U atom. The smaller clusters closely resemble the analogous UF n species, but significant deviations are seen starting with UAu 5 where a competition between U-Au and Au-Au bonding begins to become apparent. The UAu 6 system appears to mark a transition where Au-Au interactions begin to dominate, where both a ring-like and two heavily distorted octahedral structures around the central U atom are to be nearly isoenergetic. With UAu 7 , only ring-like structures are calculated. Altogether, the calculated electron detachment energies are in good agreement with the experimental values

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Stress-Alteration Enhancement of the Reactivity of Aluminum Nanoparticles in the Catalytic Decomposition of exo -Tetrahydrodicyclopentadiene (JP-10)

High energy density aluminum nanoparticles (AlNP) upon thermal annealing followed by super-quenching results in elevated stress levels in the metallic core and reduced surface energy at the core-shell interface. Isomer-selective vacuum ultraviolet based photoionization mass spectrometry coupled to a high temperature chemical microreactor reveals that these stress-altered (SA) AlNPs exhibit distinctive temperature-dependent reactivities towards catalytic decomposition of the hydrocarbon jet fuel exo-tetrahydrodicyclopentadiene (JP-10, C 10 H 16 ) compared to untreated (UN) AlNPs. SA-AlNPs show a delayed initiation of the decomposition for JP-10 by 200 K relative to the UN-AlNPs, however, the full decomposition is achieved at 100 K lower temperature. Furthermore, there is less oxygenated products that are generated from the alumina surface-induced heterogeneous oxidation process, and a larger fraction of closed and open-shell hydrocarbons. Furthermore, chemical insight bridging the reactivity order of SA-AlNPs at low and high temperatures, simultaneously, is obtained via a detailed examination of the product branching ratios obtained in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Dynamics of HCN, HNC, and HNCO Formation in the 193 nm Photodissociation of Formamide

Formamide (NH 2 CHO) is the simplest molecule containing a peptide linkage [-NH-C($=$O)-], and it plays an essential role in the study of prebiotic chemistry. Exposure to UV irradiation allows formamide to decompose and act as a prebiotic feedstock in the formation of nucleobases and other necessary starting materials. The photodissociation mechanism of gaseous formamide at 193 nm is studied using (a) chirped-pulse Fourier transform millimeter-wave spectroscopy in the 260-290 GHz spectral region in a room-temperature flow-tube reactor at 1 mu bar pressure, (b) a combination of electronic structure theory, transition state theory, and quasiclassical trajectories, and (c) the Active Thermochemical Tables. The HCN and HNC photoproducts of hydrogenated (NH 2 CHO) and deuterated (NH 2 CDO and ND 2 CHO) formamide precursors are examined to gain insight into the photodissociation mechanism. Here, the theoretical investigation has characterized the main pathway leading to each of the HCN/HNC isomers from the precursor isotopologues. The theoretical branching ratio [HNC]/[HCN] = 2.1 for nascent photofragments agrees with the experiment. The effect of the postphotolysis HNC ↔ HCN isomerization on the [HNC]/[HCN] ratio is predicted. We report the experimental branching ratio [HNCO]: ([HNC] + [HCN]) = 12 ± 3 and propose that most of HNCO originates from dissociation on the S-1 electronic state of formamide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2025

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗
DOE OSTI2020

Prediction of Structures and Atomization Energies of Coinage Metals, (M)n, n < 20: Extrapolation of Normalized Clustering Energies to Predict the Cohesive Energy

The geometries of the group 11 coinage metals (n = 2-20) were optimized to determine the lowest energy isomers for each cluster size, singlets for even numbers and doublets for odd numbers. For copper and silver, 2-D (planar) geometries were favored up to n = 6. For gold, 2D (planar) geometries were favored up to n = 13. Normalized clustering energies were plotted as a function of cluster size (n-1/3, for n = 4-20) with various DFT functionals and the CCSD(T)-F12b method and were extrapolated to predict the bulk cohesive energy. In the case of copper and silver, there is excellent agreement between the cohesive energies predicted at the CCSD(T)-F12b level of theory and the experimental values. For gold, the CCSD(T)-F12b values needed to be corrected for spin-orbit relativistic effects to obtain good agreement with experiment. Electronic properties including the HOMO-LUMO gaps for the even clusters and the spin densities for the odd clusters were calculated. The lowest gap is predicted to occur for n = 16 where the HOMO and LUMO are very similar in shape.

Persuad, Rudradatt R.↗
DOE OSTI2020

Competition between Solvation and Intramolecular Hydrogen-Bonding in Microsolvated Protonated Glycine and β-Alanine

Infrared predissociation (IRPD) spectroscopy is used to reveal and compare the microsolvation motifs of GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . The chemical structure of these amino acids differ only in the length of the carbon chain connecting the amine and carboxyl terminals, which nonetheless leads to a significant difference in the strength of the intramolecular C=O···H-N hydrogen bond in the unsolvated ions. This difference makes them useful in our studies of the competition between solvation and internal hydrogen bonding interactions. Analysis of the IRPD results reveals that the sequential addition of water molecules leads to similar effects on the intramolecular interaction in both GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . Solvation of the -NH 3 + group leads to a weakening of the C=O···H-N hydrogen bond, while solvation of the carboxyl -OH leads to a strengthening of this bond. Additionally, we have found that for β-AlaH + , the addition of a H 2 O to the second solvation shell can still influence the strength of the C=O···H-N hydrogen bonding interaction. Lastly because the C=O···H-N interaction in β-AlaH + is stronger than that in GlyH + , more solvent molecules are needed to sufficiently weaken the intramolecular hydrogen bond such that isomers without this bond begin to be energetically competitive; this occurs at n=5 for β-AlaH + and n=1 for GlyH + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Understanding the Reactivity, Selectivity, and Deactivation of Frustrated Lewis Pairs for Semihydrogenation of Acetylene

A good catalyst for semihydrogenation of alkynes must preclude both over-hydrogenation of alkene to alkane and isomerization to the other alkene isomer. In addition, it should balance the trade-off between selectivity and activity. In 2013, the Repo and Pápai groups reported a frustrated Lewis pair (FLP) (1-NMe 2 -2-B(C 6 F 5 ) 2 -C 6 H 4 ), 1, which is a metal-free catalyst and for the first time shows excellent reactivity for the hydrogenation of internal alkynes. However, it is unreactive for terminal alkynes. In this work, we have designed 13 FLPs, a–m, based on 1 by varying the Lewis base site with N and P and the Lewis acid site with B, Al, Ga, and In and replacing pentafluorophenyl with 1,3,5-trifluorophenyl, phenyl, or trifluoromethyl. We apply density functional theory to study the activity, selectivity, and deactivation of FLP 1-m for acetylene semihydrogenation. The catalytic cycle consists of three steps: (1) alkyne insertion, (2) H 2 heterolysis, and (3) intramolecular protonation. We found the activity does not change much by the modification of bulky ligands, while it decreases with the direct replacement of LA and LB sites. The overall activity depends on steps 1 and 3, which are, respectively, positively and negatively linear correlated with the charge of the Lewis acid site. Most of the FLPs in this work show comparable or better selectivity for semihydrogenation of acetylene than 1. FLP deactivation is due to the strong binding of acetylene and the elimination of electron-withdrawing bulky ligands at the preactivated catalyst rather than at activated catalysts. Furthermore, taking the selectivity and stability of FLPs into account, we predict d and k are potentially active for terminal alkynes.

02 PETROLEUM↗
DOE OSTI2023

Mechanism of Stoichiometrically Governed Titanium Oxide Brownian Tree Formation on Stepped Au(111)

Previously observed formation of substoichiometric titanium oxide dendritic structures across terraces of Au(111) is computationally studied and shown to follow the classical fractal formation mechanism of diffusion-limited aggregation (DLA). Globally optimized gas-phase oxide cluster structures are sampled in a variety of landing formations on gold surfaces and shown to favor isomers driving polymerization to Brownian tree fractal structures. Mobility of Ti 3 O 5 monomers is shown to be extremely high, with diffusion barriers of 0.21 eV or less. Through bonding stabilization, polymerization of these monomers is energetically favorable and irreversible on the 111 terrace but geometrically impossible to propagate along the step edge. Simulated scanning tunneling microscopy (STM) images show strong similarity to experiment. By contrast, observation of Ti 3 O 6 aggregating as wires along step edges is explained by the affinity of oxygen to step edges and statistical arguments for aggregation entropy at the step, in addition to low barriers for monomer diffusion and polymerization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural Evolution of Mixed-Addenda Keggin Polyoxometalate Anions with Atom-by-Atom Substitution

Polyoxometalates (POMs) are molecular metal oxides with distinctive electronic properties that make them promising materials for applications in energy, sensors, and memory devices. One of the most promising methods of tuning the stability, photochromic, redox, and electron-spin properties of POMs is through the substitution of the metal “addenda” atoms that, along with oxygen, constitute their cage-like structures. Because traditional synthesis methods typically produce a distribution of POMs, the isolation and characterization of multimetallic POMs with predetermined stoichiometry remains challenging. The presence of multiple energetically accessible isomers further complicates the experimental characterization and theoretical modeling of multimetallic POMs. Herein, we leverage the distinguishing mass-selection capabilities of ion soft landing to prepare stoichiometrically selected Keggin PMo x W 12-x O 40 3- (x = 0 – 6, 8, 10, and 12) POMs on self-assembled monolayer surfaces free of the solvent molecules and counterions that often confound characterization of complex species at interfaces. The structures of the supported POMs are characterized with atom-by-atom precision using in situ infrared (IR) reflection absorption spectroscopy complemented by detailed density functional theory calculations. Our joint experimental and theoretical results reveal an almost linear shift in the positions of the IR bands towards lower wavenumbers with an increase in the number of lighter molybdenum atoms compared to heavier W atoms in PMoxW 12-x O 40 3- . The theoretical calculations also indicate that numerous isomeric structures may be populated at the experimental conditions and, consequently, contribute to the overall IR spectra. In conclusion, our findings indicate that in addition to the number of substituted addenda atoms and the presence of multiple isomeric structures, interactions with the surface play an important role in determining the IR spectra and structure of supported bimetallic POMs.

Prabhakaran, Venkateshkumar [Pacific Northwest Nat↗
DOE OSTI2026

Transient Pulse-Response Time-of-Flight Mass Spectrometry for Complex, Deactivating Heterogeneous Catalytic Systems: Application to Ethane Dehydroaromatization

The study of complex, multistep bond-forming and -breaking reactions in heterogeneous catalytic systems often encounters challenges associated with the involvement of large numbers of intermediates among branching pathways. Kinetic information obtained from traditional steady-state measurements can be complemented with that from time-resolved methods to uncover details of the underlying chemistry. Herein, we describe an approach for tracking the complete time-resolved chemical composition (ca. 4–200 u) of a reactor effluent in response to a reactant pulse. We use a six-port rotary valve with a metered sampling loop to pulse reactants at ambient pressure into a flow reactor packed with a catalyst bed within the isothermal region of a heated furnace. The temporal evolution of effluent species is tracked using time-resolved molecular-beam time-of-flight mass spectrometry. We highlight the possibilities that this method has to offer by studying the complex bifunctional mechanism of ethane dehydroaromatization over an HZSM-5-supported platinum catalyst. We demonstrate that energy-tunable ionization sources, which facilitate isomer resolution, enable the measurement of the full mass spectral time-dependent system response. This includes the evolution of major products and mechanistically relevant reactive intermediates such as 1,3-butadiene and cyclopentadiene; these species have not previously been observed from this reaction. In additional studies, we also assess the role of platinum in the catalyst by examining temporal responses to ethane and ethylene feeds. Results show that two temporally distinct formation pathways exist for methane and benzene, and that their importance depends on both catalyst composition and reactant identity. Additionally, characteristics of catalyst deactivation are uniquely observable in the time-resolved mass spectral response, including the selective deactivation of a benzene formation pathway. The combination of time-of-flight mass spectrometry with tunable ionization enables the simultaneous observation of all effluents in a complex mixture of intermediates with isomer/isobar differentiation capabilities that can be applied to any complex reaction system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. Furthermore, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt 13 cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗
DOE OSTI2022

Manifesting Direction-Specific Complexation in [HFIP –H ·H 2 O 2 ] – : Exclusive Formation of a High-Lying Conformation

Size-selective, negative ion photoelectron spectroscopy in conjunction with quantum chemical calculations is employed to investigate the geometric and electronic structures of a protype system in catalytic olefin epoxidation research, i.e., deprotonated hexafluoroisopropanol ([HFIP -H ] - ) complexed with hydrogen peroxide (H 2 O 2 ). Spectral assignments and molecular electrostatic surface analyses unveil a surprising prevalent existence of a high-lying isomer with asymmetric dual hydrogenbonding configuration that is preferably formed driven by influential directionspecific electrostatic interactions upon H 2 O 2 approaching [HFIP -H ] - anion. Subsequent inspections of molecular orbitals, charge and spin density distributions indicate the occurrence of partial charge transfer from [HFIP -H ] - to H 2 O 2 upon hydrogen bonding interactions. Accompanying with electron detachment, a proton transfer occurs to form the neutral complex of [HFIP•HOO • ] structure. In conclusion, this work conspicuously illustrates the importance of directionality encoded in intermolecular interactions involving asymmetric and complex molecules, while the produced hydroperoxyl radical HOO • offers a possible new pathway in olefin epoxidation chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗