Search NASA⌕ Search

SEARCH · Search NASA

Results for “alkane”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Surface characterization and methane activation on SnO x /Cu 2 O/Cu(111) inverse oxide/metal catalysts

To activate methane at low or medium temperatures is a difficult task and a pre-requisite for the conversion of this light alkane into high value chemicals. In this work, we report the preparation and characterizations of novel SnO x /Cu 2 O/Cu(111) interfaces that enable low-temperature methane activation. Scanning tunneling microscopy identified small, well-dispersed SnO x nanoclusters on the Cu 2 O/Cu(111) substrate with an average size of 8 Å, and such morphology was sustained up to 450 K in UHV annealing. Ambient pressure X-ray photoelectron spectroscopy showed that hydrocarbon species (CH x groups), the product of methane activation, were formed on SnO x /Cu 2 O/Cu(111) at a temperature as low as 300 K. An essential role of the SnO x –Cu 2 O interface was evinced by the SnO x coverage dependence. Systems with a small amount of tin oxide, 0.1–0.2 ML coverage, produced the highest concentration of adsorbed CH x groups. Calculations based on density functional theory showed a drastic reduction in the activation barrier for C–H bond cleavage when going from Cu 2 O/Cu(111) to SnO x /Cu 2 O/Cu(111). On the supported SnO x , the dissociation of methane was highly exothermic (ΔE ~ –35 kcal mol –1 ) and the calculated barrier for activation (~20 kcal mol –1 ) could be overcome at 300–500 K, target temperatures for the conversion of methane to high value chemicals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lignin monomer conversion into biolubricant base oils

Despite progress in the depolymerization of lignin, only a few studies convert the obtained monomers to value-added products. Here we introduce a strategy to synthesize branched benzene lubricant (BBL) and branched cyclic lubricant (BCL) base oils from lignin-derived monomers and aldehyde. We perform carbon–carbon coupling via Brønsted acid-catalyzed hydroxyalkylation/alkylation (HAA) then hydrodeoxygenation (HDO). Optimum HAA reaction conditions achieve up to 90% guaiacol conversion and an HAA product containing 76% BBL and 24% enal condensation product over a P-SiO 2 catalyst. Subsequent HDO of HAA products over an Ir-ReO x /SiO 2 catalyst produces a lubricant-ranged mixture of BCL (C 24 ) up to yield (82%) and small fractions of dodecyl cyclohexane and C 10 and C 15 carbons alkanes. The kinematic viscosity, viscosity index, and Noack volatility of these base oils are comparable to commercial petroleum-derived poly α-olefin Group IV and refrigerant base oils. In conclusion, this approach provides a sustainable pathway for replacing petroleum-derived base oils.

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

Polyolefin plastic waste hydroconversion to fuels, lubricants, and waxes: a comparative study

Hydroconversion technologies have surged to the forefront of deconstructing plastic waste. Recent studies have been performed over several catalysts with varying conditions and plastics that make comparisons difficult. We compile and compare data from the literature by introducing various metrics and perform a simple energy analysis. We draw mechanistic similarities to and differences from the past literature on small alkane hydroconversion and leverage the former to propose standard approaches to tune product selectivity. We exemplify the plastics materials gap and the challenges it creates. Lastly, we discuss the current limitations and suggest future work.

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

Mixed oxides as multi-functional reaction media for chemical looping catalysis

Over the past two decades, chemical looping combustion (CLC) has been extensively investigated as a promising means to produce electric power while generating a concentrated carbon dioxide stream for sequestration. We note that the chemical looping strategy can be extended well outside of combustion-based carbon capture. In fact, application of the chemical looping strategy in areas beyond combustion can result in somewhat unexpected energy and carbon dioxide savings without producing a concentrated CO 2 stream at all. Furthermore, it allows the looping-based technologies to tap into applications such as chemical production – a $\$$4 trillion per year industrial sector with high energy and carbon intensities. Here, the key resides in the design of effective oxygen carriers, also known as redox catalysts in the context of selective chemical conversion through chemical looping catalysis (CLCa). This contribution focuses on the design and applications of mixed oxides as multi-function reaction media in CLCa. Since typical mixed oxide oxygen carriers tend to be nonselective for hydrocarbon conversion, the first part of this article presents generalized design principles for surface modification of mixed oxides to improve their selectivity and catalytic activity. Applications of these redox catalysts in chemical looping – oxidative dehydrogenation (CL-ODH) of a variety of light alkanes and alkyl-benzenes are presented. This is followed with a discussion of computation assisted mixed oxide design based upon thermodynamic criteria. Finally, a few new directions for the chemical looping technologies are introduced.

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

Binding of saturated and unsaturated C 6 -hydrocarbons to the electrophilic anion [B 12 Br 11 ] – : a systematic mechanistic study

The highly reactive gaseous ion [B 12 Br 11 ] – is a metal-free closed-shell anion which spontaneously forms covalent bonds with hydrocarbon molecules, including alkanes. Herein, we systematically investigate the reaction mechanism for binding of [B 12 Br 11 ] – to the five hexane isomers yielding [B 12 Br 11 (C 6 H 14 )] – , as well as to cyclohexane and several hexene isomers (yielding [B 12 Br 11 (C 6 H 12 )] – ) using collision-induced dissociation (CID), infrared photodissociation spectroscopy (IRPD) and computational methods. CID of the different [B 12 Br 11 (C 6 H 14 )] – ions results in distinct fragmentation patterns dependent on the structure of the hexane isomer. Further, the observed fragmentation reactions provide insights into the addition mechanism of [B 12 Br 11 ] – to hexane. Based on the observed CID patterns, we identified that either B–C bond formation through heterolytic C–C or C–H bond cleavages or B–H bond formation through heterolytic C–H cleavage occur dependent on the structure of the hexane isomer. Meanwhile, we observe identical CID spectra of adducts originating from isomers of C 6 H 12 . Spectroscopic investigations of adducts of 1-hexene and cyclohexane indicate the same product structure with an open C 6 chain. Computational investigations evidenced that low lying transition states are present, which enable a ring opening reaction of cyclohexane when binding to [B 12 Br 11 ] – .

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

Effect of nanostructuring on the interaction of CO 2 with molybdenum carbide nanoparticles

Transition metal carbides are increasingly used as catalysts for the transformation of CO 2 into useful chemicals. Recently, the effect of nanostructuring of such carbides has started to gain relevance in tailoring their catalytic capabilities. Catalytic materials based on molybdenum carbide nanoparticles (MoC y ) have shown a remarkable ability to bind CO 2 at room temperature and to hydrogenate it into oxygenates or light alkanes. However, the involved chemistry is largely unknown. In the present work, a systematic computational study is presented aiming to elucidate the chemistry behind the bonding of CO 2 with a representative set of MoC y nanoparticles of increasing size, including stoichiometric and non-stoichiometric cases. Furthermore, the obtained results provide clear trends to tune the catalytic activity of these systems and to move towards more efficient CO 2 transformation processes.

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

Diffusion of Hydrophilic to Hydrophobic Forms of Nile Red in Aqueous C 12 EO 10 Gels by Variable Area Fluorescence Correlation Spectroscopy

In this work, solute diffusion within lyotropic liquid crystal gels prepared from a series of water and decaethylene glycol monododecyl ether (C 12 EO 10 ) mixtures was explored by variable area fluorescence correlation spectroscopy. Aqueous C 12 EO 10 gels were prepared in concentrations ranging from 55:45 to 70:30 wt% of surfactant and water. Small angle X-ray scattering revealed that these gels comprise hexagonal mesophases of cylindrical micelles. Micelle spacing was found to decrease with increasing C 12 EO 10 concentration. Three different Nile red (NR) dyes were employed as model solutes and were separately doped into the gels at nanomolar levels. These include a hydrophilic form of NR incorporating an anionic sulfonate group (NRSO 3 - ), a hydrophobic form incorporating a fourteen-carbon alkane tail (NRC 14 ), and commercial NR as an intermediate case. FCS data acquired from the gels revealed that NRSO 3 - diffused primarily in 3D. Its diffusion coefficient exhibited a monotonic decrease with increasing gel concentration and micelle packing density, consistent with confinement of its motions as a result of its exclusion from the micelle cores. NRC 14 exhibited the smallest diffusion coefficient, most likely due to its larger size and enhanced interactions with the micelle cores. NR yielded an intermediate diffusion coefficient and the most anomalous behavior of the three dyes, attributable to its facile partitioning between core and corona regions, and greater participation by 1D diffusion. The results of these studies afford an improved understanding of molecular mass transport through soft-matter nanomaterials like those being developed for use in drug delivery and membrane based chemical separations.

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

Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetanyl peroxy radicals

2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n-pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O 2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300–1000 K and from 0.01–100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO 2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. Here, the rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.

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

Critical fluctuations in liquid–liquid extraction organic phases controlled by extractant and diluent molecular structure

Extractant aggregation in liquid–liquid extraction organic phases impacts extraction energetics and is related to the deleterious efficiency-limiting liquid–liquid phase transition known as third phase formation. Using small angle X-ray scattering, we find that structural heterogeneities across a wide range of compositions in binary mixtures of malonamide extractants and alkane diluents are well described by Ornstein–Zernike scattering. This suggests that structure in these simplified organic phases originates from the critical point associated with the liquid–liquid phase transition. To confirm this, we measure the temperature dependence of the organic phase structure, finding critical exponents consistent with the 3D Ising model. Molecular dynamics simulations were also consistent with this mechanism for extractant aggregation. Due to the absence of water or any other polar solutes required to form reverse-micellar-like nanostructures, these fluctuations are inherent to the binary extractant/diluent mixture. Further, we also show how the molecular structure of the extractant and diluent modulate these critical concentration fluctuations by shifting the critical temperature: critical fluctuations are suppressed by increasing extractant alkyl tail lengths or decreasing diluent alkyl chain lengths. This is consistent with how extractant and diluent molecular structure are known to impact metal and acid loading capacity in many-component LLE organic phases, suggesting phase behavior of practical systems may be effectively studied in simplified organic phases. Overall, the explicit connection between molecular structure, aggregation and phase behavior demonstrated here will enable the design of more efficient separations processes.

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

Carbon-negative hydrogen: aqueous phase reforming (APR) of glycerol over NiPt bimetallic catalyst coupled with CO 2 sequestration

Herein we report the production of high-pressure (19.3 bar), carbon-negative hydrogen (H 2 ) from glycerol with a purity of 98.2 mol% H 2 , 1.8 mol% light hydrocarbons (mainly methane), and 400 ppm of CO. Aqueous phase reforming (APR) of 10 wt% glycerol solution was studied with a series of NiPt alumina bimetallic catalysts supported on alumina. The Ni 8 Pt 1 -450 catalyst had the highest hydrogen selectivity (95.6%) and the lowest alkanes selectivity (3.7%) of the tested catalysts. The hydrogen selectivity decreased in the order of Ni 8 Pt 1 -450 > Ni 8 Pt 1 -260 > Ni 1 Pt 1 -260 > Pt-260. The CO 2 was sequestered with CaO adsorbent which formed CaCO 3 . We measured the adsorption capacity of the CaO adsorbent at different temperatures. Life cycle analysis showed that the APR of glycerol coupled with CO 2 capture has net negative CO 2 equivalent greenhouse gas emissions. The CO 2 emissions are –9.9 kg CO 2 eq./kg H 2 and –50.1 kg CO 2 eq./kg H 2 when grid electricity and renewable electricity are used, respectively, and the CO 2 is allocated respectively to the mass of products produced. The cost of this H 2 (denoted as “green-emerald”) was estimated to be 2.4 USD per kg H 2 when grid electricity is used and 2.7 USD per kg H 2 when using renewable electricity. The cost of glycerol has the highest contribution of 1.71 USD per kg H 2 . As a result, participation in the carbon credit markets can further decrease the price of the produced H 2 .

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

Stabilization of reactive rare earth alkyl complexes through mechanistic studies

Rare earth tris(alkyl) complexes such as M(CH 2 SiMe 3 ) 3 (sol) n are widely used as precursors for many compounds and as homogeneous catalysts for alkene polymerization and alkane functionalization. However, the thermal instability of those most conveniently made from the commercially available lithium salt of the neosilyl anion, LiCH 2 SiMe 3 , Li(r), restricts their utility. We present a new range of synthetically useful, more kinetically stable rare earth neosilyl solvates, derived from a full kinetic study of the various possible decomposition mechanisms of 7 known and 12 new solvated rare earth neosilyl complexes M(CH 2 SiMe 3 ) 3 (sol) n M = Sc(III), Y(III), Lu(III), Sm(III), and sol = THF; TMEDA; DMPE; diglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O, G 2 ), triglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 3 O, G 3 ). Surprisingly, simply using higher-denticity donors to sterically disfavor neosilyl γ-H elimination is not effective. While Sc(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) has a half-life, t 1/2 , of 258.1 h, six times longer than for Sc(r) 3 (C 4 H 8 O) 2 (t 1/2 = 43 h), Lu(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) and Y(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) do not show the expected, analogous increased t 1/2 . This is because new decomposition pathways appear for poorly fitting donors. Finally, kinetic studies demonstrate the impact of small, and increasing amounts of LiCl on the kinetics of the reactivity of the smaller alkyls Y(r) 3 (THF) 2 and Lu(r) 3 (THF) 2 ; molecules used in hydrocarbon chemistry and catalysis for fifty years. A new route to pure Y(r) 3 (THF) 2 , which avoids the traditional use of Li(r), is presented.

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

Computational insights into hydrogen adsorption energies on medium-entropy oxides

High entropy oxides (HEOs) have emerged as promising catalysts for several important chemical transformations including alkane activation. Hydrogen adsorption energy (HAE) has been used as a key descriptor for many reactions including methane C–H activation and hydrogen evolution reactions. Hence, understanding the relationship between HAEs and the surface chemistry of HEO surfaces could lay the foundation for meaningful correlations among methane C–H activation, HAE, and the complex, local environment of HEO surfaces. Here, we used a medium-entropy oxide as a prototypical system – Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O with a rock-salt structure – to interrogate these relationships. We sampled 2000 different surfaces of its (100) plane and calculated the HAEs at randomly chosen surface O sites using density functional theory (DFT). Our analysis of the 2000 data points reveals that the HAEs at the surface O sites are significantly influenced by the local environment around the adsorption sites, particularly the nature of the metal atom directly below the surface O site where H adsorbs. After comparing several popular graph-neural-network-based machine learning models, we found that the DimeNet++ model performed best achieving satisfactory accuracy in predicting HAEs for both Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O and slightly varied compositions. Our work underscores the promise of such models and the need for further refinement to address the complexity of HEOs.

Song, Haohong [Vanderbilt Univ., Nashville, TN (Un↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580–6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Kumar, Ashutosh (ORCID:0000000175896030)↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

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

Characterization of a new ultra-high pressure shock tube facility for combustion and propulsion studies

A new shock tube facility has been designed, constructed, and characterized at the University of Central Florida. This facility is capable of withstanding pressures of up to 1000 atm, allowing for combustion diagnostics of extreme conditions, such as in rocket combustion chambers or in novel power conversion cycles. For studies with toxic gas impurities, the high initial pressures required the development of a gas delivery system to ensure the longevity of the facility and the safety of the personnel. Data acquisition and experimental propagation were implemented with remote access to ensure safety, paired with a LabVIEW- and Python-based user interface. Thus far, test pressures of 270 atm, blast pressures of 730 atm, and temperatures approaching 10 000 K have been achieved. The extreme limitations of this facility allow for emission spectroscopy to be performed during the oxidation of fuel mixtures, e.g., alkanes diluted in argon and carbon dioxide. Ignition delay times were determined and compared to simulations using chemical kinetic mechanisms. The design, experimental procedures, processes of analysis, and uncertainty determination are outlined, and typical pressure profiles are compared with a new gas dynamics solver and empirical correlations developed across multiple shock tube facilities. Preliminary reactive mixture analyses are included with further investigation of the mixtures outlined.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

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

Interpretation of autoencoder-learned collective variables using Morse–Smale complex and sublevelset persistent homology: An application on molecular trajectories

Dimensionality reduction often serves as the first step toward a minimalist understanding of physical systems as well as the accelerated simulations of them. In particular, neural network-based nonlinear dimensionality reduction methods, such as autoencoders, have shown promising outcomes in uncovering collective variables (CVs). However, the physical meaning of these CVs remains largely elusive. In this work, we constructed a framework that (1) determines the optimal number of CVs needed to capture the essential molecular motions using an ensemble of hierarchical autoencoders and (2) provides topology-based interpretations to the autoencoder-learned CVs with Morse–Smale complex and sublevelset persistent homology. Furthermore, this approach was exemplified using a series of n-alkanes and can be regarded as a general, explainable nonlinear dimensionality reduction method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗