Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multidimensional”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2018

Classical Optimal Control for Energy Minimization Based On Diffeomorphic Modulation under Observable-Response-Preserving Homotopy

In this work, we introduce the so-called “Classical Optimal Control Optimization” (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass $m(t;β)$ and dipole $μ(r,t;β)$ is evolved classically on the potential energy surface $V(r)$ coupled to an electric field $E(t;β)$, as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters $β$ defining $m(t;β), μ(r,t;β),$ and $E(t;β)$ are optimized by following the gradients of the energy with respect to $β,$ adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.

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

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

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

High-Temperature Dissociation of Neopentanol: Shock Tube/Photoionization Mass Spectrometry Studies

Herein, the pyrolysis mechanism of 2,2-dimethylpropan-1-ol (neopentanol) has been investigated at high temperatures (1128–1401K) and high pressures (5 and 15 bar). The experiments were performed in a miniature shock tube coupled to a time-of-flight mass spectrometer. Cations were generated by tunable vacuum ultraviolet photoionization resulting in multidimensional data sets containing mass and photoionization spectra and the time histories of species. At the elevated temperatures and pressures of this work, neopentanol was determined to dissociate primarily by the scission of a C–C bond yielding tert-butyl and hydroxymethyl radicals. These promptly form isobutene and formaldehyde by H-atom elimination. In the structurally similar molecule neopentane, roaming radical reactions have previously been found to be important under conditions close to the present work (1260–1459 K, 1.1 bar). There are two possible roaming radical reactions for neopentanol. However, no experimental evidence for these reactions was found at the elevated pressures in this study, and the dissociation of neopentanol is dominated by bond scission yielding radical products.

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

MQCT. I. Inelastic Scattering of Two Asymmetric-Top Rotors with Application to H 2 O + H 2 O

A mixed quantum/classical theory (MQCT) for the inelastic collision of two asymmetric-top rotor molecules is developed. Here, the quantum state-to-state transitions between the rotational states of molecules (internal) are treated quantum mechanically using the time-dependent Schrodinger equation, whereas their relative translational motion (responsible for scattering) is treated classically, using the average trajectory approach. Two versions of the formula for transition matrix elements are presented: a straightforward approach that uses numerical multidimensional quadrature over all the internal degrees of freedom and a more standard analytic approach that uses the expansion of the PES over the basis set of spherical harmonics. Adaptation to the case of identical molecules scattering is presented and is applied to the rotational excitation of two water molecules, H 2 O + H 2 O, using the PES from recent literature. Calculations of collisional excitation from the ground state of the system into a number of low-lying excited rotational states are carried out in a broad range of energies. Analysis of computed opacity functions shows a rather unusual scattering regime, dominated by a strong anisotropic long-range interaction (dipole–dipole). The coupled-states (CS) approximation is tested and found to agree semiquantitatively (within a factor of 2) with the fully coupled version of the method. Differential cross sections for the elastic scattering suggest a very narrow forward scattering peak.

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

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Phylloquinone (PhQ) is a molecule involved in photosynthetic electron transfer, where the local protein environment is known to tune the properties of PhQ through various noncovalent interactions. In this work, we determine how hydrogen bonding, one of these noncovalent interactions, alters the vibrational potential energy surface of PhQ. In addition, this work demonstrates how hydrogen bonding to PhQ manifests in ultrafast multidimensional vibrational spectra, which is an important step towards using ultrafast vibrational based spectroscopies to probe PhQ embedded in the binding pocket of photosynthetic reaction centers.

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

Molecular dynamics simulations of hydrogen diffusion in aluminum

In this study, hydrogen diffusion impacts the performance of solid-state hydrogen storage materials and contributes to the embrittlement of structural materials under hydrogen-containing environments. In atomistic simulations, the diffusion energy barriers are usually calculated using molecular statics simulations where a nudged elastic band method is used to constrain a path connecting the two end points of an atomic jump. This approach requires prior knowledge of the “end points”. For alloy and defective systems, the number of possible atomic jumps with respect to local atomic configurations is tremendous. Even when these jumps can be exhaustively studied, it is still unclear how they can be combined to give an overall diffusion behavior seen in experiments. Here we describe the use of molecular dynamics simulations to determine the overall diffusion energy barrier from the Arrhenius equation. This method does not require information about atomic jumps, and it has additional advantages, such as the ability to incorporate finite temperature effects and to determine the pre-exponential factor. As a test case for a generic method, we focus on hydrogen diffusion in bulk aluminum. We find that the challenge of this method is the statistical variation of the results. However, highly converged energy barriers can be achieved by an appropriate set of temperatures, output time intervals (for tracking hydrogen positions), and a long total simulation time. Our results help elucidate the inconsistencies of the experimental diffusion data published in the literature. The robust approach developed here may also open up future molecular dynamics simulations to rapidly study diffusion properties of complex material systems in multidimensional spaces involving composition and defects.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Comparison of Anion Reorientational Dynamics in MCB 9 H 10 and M 2 B 10 H 10 (M = Li, Na) via Nuclear Magnetic Resonance and Quasielastic Neutron Scattering Studies

The disordered phases of the 1-carba-closo-decaborates LiCB9H10 and NaCB9H10 exhibit the best solid-state ionic conductivities to date among all known polycrystalline competitors, likely facilitated in part by the highly orientationally mobile CB9H10- anions. We have undertaken both NMR and quasielastic neutron scattering (QENS) measurements to help characterize the monovalent anion reorientational mobilities and mechanisms associated with these two compounds and to compare their anion reorientational behaviors with those for the divalent B10H102- anions in the related Li2B10H10 and Na2B10H10 compounds. NMR data show that the transition from the low-T ordered to the high-T disordered phase for both LiCB9H10 and NaCB9H10 is accompanied by a nearly two-orders-of-magnitude increase in the reorientational jump rate of CB9H10- anions. QENS measurements of the various disordered compounds indicate anion jump correlation frequencies on the order of 1010-1011 s-1 and confirm that NaCB9H10 displays jump frequencies about 60% to 120% higher than those for LiCB9H10 and Na2B10H10 at comparable temperatures. The Q-dependent quasielastic scattering suggests similar small-angular-jump reorientational mechanisms for the different disordered anions, changing from more uniaxial in character at lower temperatures to more multidimensional at higher temperatures, although still falling short of full three-dimensional rotational diffusion below 500 K within the nanosecond neutron window.

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

Photoinduced Charge Transfer and Vibronic Coherence in CdSe Quantum Dots with Methyl Viologen Acceptors

We show herein that photoinduced charge transfer from CdSe quantum dots (QDs) to surface-bound methyl viologen (MV 2+ ) acceptors is mediated by a vibronically coherent, nonadiabatic mechanism. Broadband multidimensional electronic spectra and an analysis of coherences show that a mixed QD–MV charge-transfer (CT) state is populated on the <50 fs time scale after optical preparation of the X3 (1P 3/2 –1P e ) state, well prior to the appearance of the one-electron photoreduced ground state (MV +• ). A partial redistribution of charge from the core of the QD to the acceptor is revealed by excited-state coherences of an out-of-plane vibrational mode local to MV 2+ and of a low-frequency mode mixing a MV 2+ mode with the longitudinal optical (LO) phonon of the QD core. The ultrafast damping of these coherences indicates that excited-state wavepackets travel from the optically prepared, Franck–Condon structure through a conical intersection to reach the CT state. These results suggest that vibronically coherent processes generating CT intermediates can be exploited to improve the efficiency of QD-based solar cells and photocatalysts.

Cadmium selenide↗
DOE OSTI2018

Tracking the Progression of Anion Reorientational Behavior between α-Phase and β-Phase Alkali-Metal Silanides, MSiH 3 , by Quasielastic Neutron Scattering

Quasielastic neutron scattering (QENS) measurements over a wide range of energy resolutions were used to probe the reorientational behavior of the pyramidal SiH 3 - anions in the monoalkali silanides (MSiH 3 , where M = K, Rb, and Cs) within the low-temperature ordered β-phases, and for CsSiH 3 , the high-temperature disordered α-phase and intervening hysteretic transition region. Maximum jump frequencies of the β-phase anions near the β-α transitions range from around 10 9 s -1 for β-KSiH 3 to 10 10 s -1 and higher for β-RbSiH 3 and β-CsSiH 3 . The β-phase anions undergo uniaxial 3-fold rotational jumps around the anion quasi-C 3 symmetry axis. CsSiH 3 was the focus of further studies to map out the evolving anion dynamical behavior at temperatures above the β-phase region. As in α-KSiH 3 and α-RbSiH 3 , the highly mobile anions (with reorientational jump frequencies approaching and exceeding 10 12 s -1 ) in the disordered α-CsSiH 3 are all adequately modeled by H jumps between 24 different locations distributed radially around the anion center of gravity, although even higher anion reorientational disorder cannot be ruled out. QENS data for CsSiH 3 in the transition region between the α- and β-phases corroborated the presence of dynamically distinct intermediate ( i -) phase. The SiH 3 - anions within i -phase appear to undergo uniaxial small-angular-jump reorientations that are more akin to the lower-dimensional β-phase anion motions rather than to the multidimensional α-phase anion motions. Moreover, they possess orientational mobilities that are an order-of-magnitude lower than those for α-phase anions but also an order-of-magnitude higher than those for β-phase anions. Combined QENS and neutron powder diffraction results strongly suggest that this i -phase is associated chiefly with the more short-range-ordered, nanocrystalline portions (invisible to diffraction) that appear to dominate the CsSiH 3 .

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

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. Here, we show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

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

Visualizing Stereodynamics in Cold Collisions through Shape Resonance Wavefunctions

Shape resonances in cold collisions are often strongly affected by stereodynamics. These resonances can sometimes be assigned by single-channel quantization along the scattering coordinate. However, sensitive steric control of collision implies strong anisotropy of the interaction potential energy surface, which usually leads to coupling among multiple scattering channels in the strongly interacting region. Hence, the resonances might be the result of quantization in a multidimensional space. Using the cold rotationally inelastic collision between para-H 2 (v 1 = 1, j 1 = 2) and HF (v 2 = 0, j 2 = 0) as an example, we analyze four low-lying shape resonances via diagonalization of the full-dimensional Hamiltonian with a stabilization method. While some resonances can indeed be assigned with a single partial wave, others apparently involve more than one scattering channel. Furthermore, a new model based on these resonance wavefunctions is developed to better understand the stereodynamics through shape resonances in cold collisions.

Collisions↗
DOE OSTI2016

Optical Dedoping Mechanism for P3HT:F4TCNQ Mixtures

Doping-induced solubility control (DISC) is a recently introduced photolithographic technique for semiconducting polymers, which utilizes reversible changes in polymer solubility upon doping to allow the polymer to function as its own photoresist. Central to this process is a wavelength sensitive optical dedoping reaction, which is poorly understood but generates subdiffraction-limited topographic features and provides optical control of the polymer doping level. In this work, we examine the mechanism of optical dedoping in the semiconducting polymer poly-3-hexylthiophene (P3HT) doped by 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ), via a combination of ultrafast and steady-state spectroscopy, ab initio calculations, and multidimensional NMR. A simple photoinduced back electron transfer mechanism from reduced F4TCNQ to oxidized P3HT does not explain the observed photophysics. Instead, photoexcited F4TCNQ* reacts with THF solvent molecules to form a neutral, nondoping, and highly soluble F4TCNQ-THF complex. Hence, ionized F4TCNQ is removed from the P3HT indirectly by depletion of the neutral F4TCNQ. Because the reaction involves only the dopant and similar photoreactivity would expected for most other dopant molecules, we expect optical DISC patterning should be generalizable to a wide range of polymer:dopant systems.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Machine-Learning Enabled New Insights into the Coil-to-Globule Transition of Thermosensitive Polymers Using a Coarse-Grained Model

In this paper, we present a computational framework that integrates coarse-grained (CG) molecular dynamics (MD) simulations and a data-driven machine-learning (ML) method to gain insights into the conformations of polymers in solutions. We employ this framework to study conformational transition of a model thermosensitive polymer, poly(N-isopropylacrylamide) (PNIPAM). Here, we have developed the first of its kind, a temperature-independent CG model of PNIPAM that can accurately predict its experimental lower critical solution temperature (LCST) while retaining the tacticity in the presence of an explicit water model. The CG model was extensively validated by performing CG MD simulations with different initial conformations, varying the radius of gyration of chain, the chain length, and the angle between the adjacent monomers of the initial configuration of PNIPAM (total simulation time = 90 μs). Moreover, for the first time, we utilize the nonmetric multidimensional scaling (NMDS) method, a data-driven ML approach, to gain further insights into the mechanisms and pathways of this coil-to-globule transition by analyzing CG MD simulation trajectories. In conclusion, NMDS analysis provides entirely new insights and shows multiple metastable states of PNIPAM during its coil-to-globule transition above the LCST.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Revealing Hot and Long-Lived Metastable Spin States in the Photoinduced Switching of Solvated Metallogrid Complexes with Femtosecond Optical and X-ray Spectroscopies

An atomistic understanding of the photoinduced spin-state switching (PSS) within polynuclear systems of d 4 –d 7 transition metal ion complexes is required for their rational integration into light-driven reactions of chemical and biological interest. However, in contrast to mononuclear systems, the multidimensional dynamics of the PSS in solvated molecular arrays have not yet been elucidated due to the expected complications associated with the connectivity between the metal centers and the strong interactions with the surroundings. Herein, the PSS in a solvated triiron(II) metallogrid complex is characterized using transient optical absorption and X-ray emission spectroscopies on the femtosecond time scale. The complementary measurements reveal the photoinduced creation of energy-rich (hot) and long-lived quintet states, whose dynamics differ critically from their mononuclear congeners. We report this finding opens major prospects for developing novel schemes in solution-phase spin chemistry that are driven by the dynamic PSS process in compact oligometallic arrays.

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

A Preprocessing Tool for Enhanced Ion Mobility–Mass Spectrometry-Based Omics Workflows

The ability to improve the data quality of ion mobility–mass spectrometry (IM-MS) measurements is of great importance for enabling modular and efficient computational workflows and gaining better qualitative and quantitative insights from complex biological and environmental samples. We developed the PNNL PreProcessor, a standalone and user-friendly software housing various algorithmic implementations to generate new MS-files with enhanced signal quality and in the same instrument format. Different experimental approaches are supported for IM-MS based on Drift-Tube (DT) and Structures for Lossless Ion Manipulations (SLIM), including liquid chromatography (LC) and infusion analyses. The algorithms extend the dynamic range of the detection system, while reducing file sizes for faster and memory-efficient downstream processing. Specifically, multidimensional smoothing improves peak shapes of poorly defined low-abundance signals, and saturation repair reconstructs the intensity profile of high-abundance peaks from various analyte types. Further, other functionalities are data compression and interpolation, IM demultiplexing, noise filtering by low intensity threshold and spike removal, and exporting of acquisition metadata. Several advantages of the tool are illustrated, including an increase of 19.4% in lipid annotations and a two-times faster processing of LC-DT IM-MS data-independent acquisition spectra from a complex lipid extract of a standard human plasma sample. The software is freely available at https://omics.pnl.gov/software/pnnl-preprocessor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Photoinitiated Olefin Metathesis and Stereolithographic Printing of Polydicyclopentadiene

Recent progress in photoinitiated ring-opening metathesis polymerization (photoROMP) has enabled the lithographic production of patterned films from olefinic resins. Recently, we reported the use of a latent ruthenium catalyst (HeatMet) in combination with a photosensitizer (2-isopropylthioxanthone) to rapidly photopolymerize dicyclopentadiene (DCPD) formulations upon irradiation with UV light. While this prior work was limited in terms of catalyst and photosensitizer scope, a variety of alternative catalysts and photosensitizers are commercially available that could allow for tuning of thermomechanical properties, potlifes, activation rates, and irradiation wavelengths. Herein, 14 catalysts and 8 photosensitizers are surveyed for the photoROMP of DCPD and the structure–activity relationships of the catalysts examined. Properties relevant to stereolithography additive manufacturing (SLA AM)-potlife, irradiation dose required to gel, conversion-are characterized to develop catalyst and photosensitizer libraries to inform development of SLA AM resin systems. Two optimized catalyst/photosensitizer systems are demonstrated in the rapid SLA printing of complex, multidimensional pDCPD structures with microscale features under ambient conditions.

36 MATERIALS SCIENCE↗