Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer systems organization”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An empirical investigation of organic software product lines

Abstract Software product line engineering is a best practice for managing reuse in families of software systems that is increasingly being applied to novel and emerging domains. In this work we investigate the use of software product line engineering in one of these new domains, synthetic biology. In synthetic biology living organisms are programmed to perform new functions or improve existing functions. These programs are designed and constructed using small building blocks made out of DNA. We conjecture that there are families of products that consist of common and variable DNA parts, and we can leverage product line engineering to help synthetic biologists build, evolve, and reuse DNA parts. In this paper we perform an investigation of domain engineering that leverages an open-source repository of more than 45,000 reusable DNA parts. We show the feasibility of these new types of product line models by identifying features and related artifacts in up to 93.5% of products, and that there is indeed both commonality and variability. We then construct feature models for four commonly engineered functions leading to product lines ranging from 10 to 7.5 × 10 20 products. In a case study we demonstrate how we can use the feature models to help guide new experimentation in aspects of application engineering. Finally, in an empirical study we demonstrate the effectiveness and efficiency of automated reverse engineering on both complete and incomplete sets of products. In the process of these studies, we highlight key challenges and uncovered limitations of existing SPL techniques and tools which provide a roadmap for making SPL engineering applicable to new and emerging domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A mathematical analysis of carbon fixing materials that grow, reinforce, and self-heal from atmospheric carbon dioxide

Carbon fixing materials are a new class of self-healing, self-reinforcing materials we have introduced that utilize ambient CO 2 to chemically add to an ever extending carbon backbone. This class of materials can utilize biological or non-biological photocatalysts and support a wide range of potential backbone chemistries. However, there is no analysis to date that describes their fundamental limits in terms of chemical kinetics and mass transfer. In this computational study, we employ a reaction engineering, and materials science analysis to answer basic questions about the maximum growth rate, photocatalytic requirements and limits of applicable materials. Our proposed mathematical framework envelops three main functions required for carbon fixing materials: (1) adsorption of CO 2 from air, (2) photocatalytic reduction of CO 2 into selected monomers, and (3) polymerization of CO 2 -derived monomers. First, by performing a Damköhler number analysis, we derive criteria for the cross over from kinetic control to mass transfer limited growth, setting upper limits on performance of a potential photocatalyst. Next, we analyze photocatalytic reduction of CO 2 to single carbon products, using known catalytic pathways and kinetic data. We identify formaldehyde as a C1 intermediate having unique potential for incorporation into the material backbone of carbon fixing materials. As an example, we find that a diamond-shaped reaction network graph for CO 2 reduction, passing through CO and HCOOH intermediates, accurately describes kinetic data for cobalt-promoted TiO 2 nanoparticles at room temperature and 1 atm CO 2 . Finally, as an applied case study, we consider and analyze the example of a carbon fixing poly(oxymethylene) system with embedded catalyst promoting the photocatalytic reduction of CO 2 to formaldehyde. The latter then trimerizes to a trioxane monomer which subsequently polymerizes to polyoxymethylene. This reaction engineering analysis introduces benchmarks for carbon fixing materials with respect to achievable rates of photocatalysis, adsorption, and polymerization. These results should prove valuable for the design, evaluation, and benchmarking of this emergent and new class of environmentally sustainable materials.

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

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Elucidating the Role of Electric Fields in Fe Oxidation via an Environmental Atom Probe

We quantify the effects of intensely applied electric fields on the Fe oxidation mechanism. The specimen are pristine Fe single crystals exposing a variety of surface structures identified by field ion microscopy. These crystals are simultaneously exposed to low pressures of pure oxygen gas, on the order of 10 −7 mbar, while applying intense electric fields on their surface of several tens of volts per nanometer. The local composition of the different surface structures is probed directly and in real time using an Environmental Atom Probe and successfully compared with first principles-based models. We found that rough Fe{244} and Fe{112} facets are more reactive toward oxygen than compact Fe{024} and Fe{011} facets. Results demonstrate that the influence of an electric field on the oxidation kinetics depends on the timescales that are involved as the system evolves toward equilibrium. The initial oxidation kinetics show that strong increases in electric fields facilitate the formation of an oxide. However, as one approaches equilibrium, high field values mitigate this formation. Ultimately, this study elucidates how high externally applied electric fields can be used to dynamically exploit reaction dynamics at the nanoscale towards desired products in a catalytic reaction at mild reaction conditions.

09 BIOMASS FUELS↗
DOE OSTI2022

The role of energy transfer and competing bimolecular reactions in characterizing the unimolecular dissociations of allylic radicals

Resonantly stabilized radicals (RSR's) such as allylic radicals typically have bond energies exceeding 45 kcal/mol and are therefore reasonably long-lived up to quite high temperatures. Consequently, kinetics of the unimolecular dissociations and bimolecular reactions of these long-lived RSR's assume significance in high temperature chemically reacting systems such as in flames. Our recent analyses (J. Cho et al., Proc. Comb. Inst., 2022) of the uncertainty propagated by the energy transfer parameters for 1-methylallyl (1MA) dissociation to the laminar flame speed, indicates an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Here, this work extends such analyses to other C 3 – C 5 allylic radicals: allyl, 2-methylallyl, 1,1-dimethylallyl, 1,2-dimethylallyl, and 1,3-dimethylallyl. Theoretical kinetics for the C 3 – C 4 allylic radicals were taken from prior literature studies (J.A. Miller et al., J. Phys. Chem., 2008; J. Cho et al., Proc. Comb. Inst., 2022; R.S. Tranter et al., Proc. Comb. Inst., 2017). Potential energy surfaces were newly computed for the dissociations of the larger C5 allylic radicals and the reactions of these radicals with H-atoms. These PES's were then used to calculate rate constants for the relevant unimolecular and bimolecular reactions. Energy transfer has a crucial role to play in these pressure-dependent reactions and therefore the present kinetics calculations relied on a-priori predictions of collisional energy transfer parameters to best characterize the branching between chain propagating unimolecular dissociations of these allylic radicals and chain terminating bimolecular reactions of these RSR's with H-atoms. Prompt dissociations are now established to be a universal feature of free radical dissociations. Here we have also characterized the propensity for prompt dissociations in these more-stable allylic RSR's and their role in combustion simulations.

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

Preventing Loss of Selectivity during the Oxidative Dehydrogenation of Propane over Supported Vanadium Catalysts

Supported vanadium materials are promising catalysts for the oxidative dehydrogenation of propane to propylene (ODHP), but a lack of mechanistic understanding limits the rational design of catalysts with improved propylene selectivity. Adding Ta to V/SiO 2 increases the propylene selectivity, as well as the activity, leading to superior performance compared to state-of-the-art boron-based systems. In this contribution, we utilize this surprising promotional effect of Ta to elucidate key elements of the mechanistic cycle. Through a combination of characterization techniques, computational modeling, and kinetic experiments, we show that the catalytic cycle over V/SiO 2 likely involves the formation of an isopropyl alcohol intermediate, the fate of which is in kinetic competition between subsequent dehydration to propylene or further oxidation. Furthermore, we show that the relatively facile propylene overoxidation observed for these materials occurs via the epoxidation of propylene by a proposed peroxovanadium intermediate, rather than the abstraction of propylene’s allylic C–H bond as previously assumed. Using these key mechanistic features, we rationalize the enhanced selectivity and activity of Ta promotion. In conclusion, our mechanistic framework offers avenues for future catalyst development to improve supported vanadium materials for ODHP.

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

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

California Energy Systems for the 21st Century (CES-21) Program, CRADA TC02200 (Final Report)

The following report representsthe culmination of a five-year Cybersecurity Research and Development (R&D) program performed by California’s Investor-Owned Utilities (IOUs) and Lawrence Livermore National Laboratory (LLNL) and funded by the IOUs’ electricity ratepayers as authorized by California Legislation and the California Public Utilities Commission (CPUC). The cybersecurity R&D, which was intensely technical in nature, was broken into three major workstreams: 1. The development of a modeling & simulation platform, to explore the potential effects of various threat and response scenarios at grid scale; 2. The establishment of a physical testbed with separate substation instances from each of the IOUs, to evaluate threats and responses on actual substation equipment; 3. The development of a research package consisting of several capabilities to support the industry’s evolution towards automated threat response and other next-generation cybersecurity techniques Throughout the program, there was extensive collaboration between the program team and national laboratories, federal departments, academic institutions and industry organizations. Several of the tools developed through the program have been made available to the open source community, to enable faster adoption and continued development of important cybersecurity capabilities. While this program began to develop much of the foundation for automated threat response, much work remains to be done, and this report recommends a series of next steps. While the Executive Summary below has been simplified as much as possible, the very technical nature of the subject matter would make further reduction difficult without causing meaning to be lost. Due to the sensitive nature of the research, some topics are intentionally covered at a high level to protect the research and therefore the electric grid and customers of California and the Nation. For further detail, clarification or questions we suggest discussion with the program team.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Solvent and Temperature Effects on Photoinduced Proton-Coupled Electron Transfer in the Marcus Inverted Region

Concerted proton-coupled electron transfer (PCET) in the Marcus inverted region was recently demonstrated. Understanding the requirements for such reactivity is fundamentally important and holds promise as a design principle for solar energy conversion systems. Herein, we investigate the solvent polarity and temperature dependence of photoinduced proton-coupled charge separation (CS) and charge recombination (CR) in anthracene–phenol–pyridine triads: 1 (10-(4-hydroxy-3-(4-methylpyridin-2-yl)benzyl)anthracene-9-carbonitrile) and 2 (10-(4-hydroxy-3-(4-methoxypyridin-2-yl)benzyl)anthracene-9-carbonitrile). Both the CS and CR rate constants increased with increasing polarity in acetonitrile:n-butyronitrile mixtures. The kinetics were semi-quantitatively analyzed where changes in dielectric and refractive index, and thus consequently changes in driving force (–ΔG°) and reorganization energy (λ), were accounted for. The results were further validated by fitting the temperature dependence, from 180 to 298 K, in n-butyronitrile. The analyses support previous computational work where transitions to proton vibrational excited states dominate the CR reaction with a distinct activation free energy (ΔG* CR ~ 140 meV). However, the solvent continuum model fails to accurately describe the changes in ΔG° and λ with temperature via changes in dielectric constant and refractive index. Satisfactory modeling was obtained using the results of a molecular solvent model, which predicts that λ decreases with temperature, opposite to that of the continuum model. To further assess the solvent polarity control in the inverted region, the reactions were studied in toluene. Nonpolar solvents decrease both ΔG° CR and λ, slowing CR into the nanosecond time regime for 2 in toluene at 298 K. This demonstrates how PCET in the inverted region may be controlled to potentially use proton-coupled CS states for efficient solar fuel production and photoredox catalysis.

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

A guided ion beam investigation of UO 2 + thermodynamics and f orbital participation: Reactions of U + + CO 2 , UO + + O 2 , and UO + + CO

A guided ion beam tandem mass spectrometer was employed to study the reactions of U + + CO 2 , UO + + O 2 , and the reverse of the former, UO + + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U + + CO 2 proceeds to form UO + + CO with an efficiency of 118% ± 24% as well as generating UO 2 + + C and UCO + + O. The reaction of UO + + O 2 forms UO 2 + in an exothermic, barrierless process and also results in the collision-induced dissociation of UO + to yield U + . In the UO + + CO reaction, the formation of UO 2 + in an endothermic process is the dominant reaction, but minor products of UCO + + O and U + + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D 0 (U + –O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D 0 (U + –CO) = 0.73 ± 0.13 eV, and D 0 (OU + –O) = 7.56 ± 0.12 eV. The values obtained for D 0 (U + –O) and D 0 (OU + –O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D 0 (U + –CO). Furthermore, an analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO 2 + relative to ThO 2 + and especially transition metal dioxide cations.

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

Toward the Development of a High-Voltage Mg Cathode Using a Chromium Sulfide Host

Development of Mg-ion batteries as advanced electrochemical energy storage systems relies on the design and discovery of high-voltage positive electrode (cathode) materials. To date, a variety of sulfide cathodes have been reported (e.g., Mg x Mo 6 S 8 , Mg x Ti 2 S4, etc.), but the voltages of these materials are too low to prepare a high energy density Mg cell. Theoretical computations predicted that Mg x Cr 2 S 4 operating with the high-voltage Cr 3+/4+ redox couple would serve as a suitable cathode candidate, but experimental attempts to extract Mg 2+ from the lattice have been largely unsuccessful. Here, we show that reversible electrochemical activity within a thiospinel framework (AB 2 S 4 ) relies on a redox-active transition metal present in the B site; otherwise, anionic redox activity triggers decomposition of the spinel structure. Since Cr and S states are highly coupled in MgCr 2 S 4 , the Cr 3+/4+ redox couple is inaccessible so that reversible (de)intercalation of Mg 2+ cannot occur and charging leads to dissolution of the active material. These findings point to an insufficiency in the screening criteria previously used to identify MgCr 2 S 4 as a promising Mg cathode. Thus, a computable descriptor based on the electronic structure of the discharged material is proposed to predict the prevalence of cation vs anion redox and improve future surveys of potential candidates. It is unlikely that the high-voltage Cr redox couple will be accessible to oxidation in the presence of sulfur within the restrictions of a spinel framework; however, it is possible that a suitable layered Mg x CrS 2 structure could serve as a reversible high-voltage Mg cathode.

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

Chemical-Free Lithium Separation from High-Salinity Brines Using Model-Informed and Machine Learning-Optimized Multi-Column Zwitterionic Chromatography

Direct Lithium Extraction (DLE) technologies often struggle to produce high-purity lithium salts from high-salinity brines, as current approaches require chemical-based elution, regeneration, and precipitation steps, resulting in significant environmental footprints. A novel salt fractionation approach using carboxybetaine resin, known as zwitterionic chromatography (ZIC), has demonstrated that lithium ions can be separated from divalent cations under high-salinity conditions using only water as eluent, with no regeneration required. To enable continuous and scalable deployment of this approach, we developed a chemical-free Multi-column Zwitterionic Chromatography (MZC) process and its theoretical and process models. To predict and optimize this nontraditional separation system, we introduced a novel anti-Langmuir isotherm, and the isotherm parameters were estimated through a machine learning-driven optimization based on artificial neural network ensembles with numerical feasibility assessment. Using machine learning-driven optimization, the MZC process achieved 98.0% lithium recovery, 99.5 % Li/(Li + Mg + Ca) purity, a 31.3% productivity increase, and a 33% reduction in water use compared to batch operation. The proposed MZC process enables lithium separation at $0.6-1.2 kg-1 Li, with costs dominated by resin manufacturing, while offering lower separation costs and carbon footprint compared with conventional carbonation. Overall, these findings position the MZC process as an effective polishing step within scalable and sustainable lithium production pipelines.

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

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

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

Harnessing dimethyl ether with ultra-low-grade heat for scaling-resistant brine concentration and fractional crystallization

Solvent-driven separations may enable scalable concentration of hypersaline brines, supporting a circular resource economy from the extraction of lithium and rare earth elements from spent battery and magnet leachates. This work analyses a novel solvent-driven water extraction (SDWE) system employing dimethyl ether (DME) and ultra-low-grade heat for brine concentration and fractional crystallization. SDWE exploits DME’s unique properties: (1) a low dielectric constant that promotes water solubility over charged solutes by a factor of 10 3 , and (2) a high volatility that facilitate efficient DME reconcentration with ultra-low-grade heat. The techno-economic viability of SDWE is assessed with a computational framework that encompasses a liquid–liquid separator and a solvent concentrator. We integrate the extended universal quasichemical model with the virial equation of state to predict the compositions of the complex three-phase DME-water mixture at vapor–liquid and liquid–liquid equilibrium. Subsequently, we optimize the thermodynamic and economic performance of SDWE, by controlling the interstage flash pressure, heat source temperature, and the number of concentrating stages. DME-based SDWE concentrates an input saline feed to 5.5 M and regenerates over 99 % of the DME using ultra-low-grade heat below 50 °C, with a DME/water selectivity ratio of 125. Here our calculations reveal that optimal performance is achieved at interstage flash pressures of 0.4 – 0.5 bar for heat source temperatures between 323–373 K, with improved exergetic efficiencies at lower temperatures. At a heat source temperature of 323 K and an interstage pressure of 0.489 bar, DME-driven SDWE achieves an optimal thermodynamic efficiency of 20.5 % and a projected specific cost of US$ 1.93 m -3 . These specific costs suggest that SDWE is competitive with commercialized thermal distillation technologies, while mitigating the traditional risks associated with scaling in heat and mass exchangers with hypersaline brines.

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

Influence of Terminal Carboxyl Groups on the Structure and Reactivity of Functionalized m-Carboranethiolate Self-Assembled Monolayers

The structure and function of self-assembled monolayers (SAMs) at the nanoscale are determined by the steric and electronic effects of their building blocks. Carboranethiol molecules form pristine monolayers that provide tunable two-dimensional systems to probe lateral and interfacial interactions. Additional ω-functionality, such as carboxyl groups, can be introduced to change the properties of the exposed surfaces. Here, two geometrically similar isomeric m-carborane analogues of m-mercaptobenzoic acid, 1-COOH-7-SH-1,7-C 2 B 10 H 10 and racem-1-COOH-9-SH-1,7-C 2 B 10 H 10 , are characterized and their SAMs on Au{111} are examined. The latter isomer belongs to the rare group of chiral cage molecules and becomes, to our knowledge, the first example assembled on Au{111}. Although different in symmetry, molecules of both isomers assemble into similar hexagonal surface patterns. The nearest-neighbor spacing of 8.4 ± 0.4 Å is larger than that of non-carboxylated isomers, consistent with the increased steric demands of the carboxyl groups. Computational modeling reproduced this spacing and suggests a tilt relative to the surface normal. However, tilt domains are not observed experimentally, suggesting the presence of strong lateral interactions. Analyses of the influence of the functional groups through the pseudo-aromatic m-carborane skeleton showed that the thiol group attached to either carbon or boron atoms increases the carboxyl group acidity in solution. In contrast, the acidity of the exposed carboxyl group in the SAMs decreases upon surface attachment; computational analyses suggest that the driving force of this shift is the dielectric of the environment in the monolayer as a result of confined intermolecular interactions, proximity to the Au surface, and partial desolvation.

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

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗