Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamics (MD)”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Relaxation dynamics of water in the vicinity of cellulose nanocrystals

Water molecules near cellulose nanocrystals (CNCs; produced via the sulfuric acid-catalyzed hydrolysis of wood pulp) are believed to relax slower than those in the bulk liquid, which may result in unique properties of CNC aqueous dispersions. Herein this study analyzed the polarization behavior of water molecules in CNC aqueous dispersions and other reference samples using a dielectric relaxation spectroscopy (DRS) technique in the microwave frequency range (0.2–20 GHz). As the CNC concentration increases, two slow relaxation components become prominent. The comparison with DRS data of aqueous dispersions of nanoporous silica, polyvinyl alcohol (PVA), and hairy CNC (HCNC) with amorphous chains protruding from both ends suggested that these slow relaxation modes of water near CNC surfaces cannot be attributed to direct hydrogen bonding interactions with the hydroxyl (OH) groups exposed and immobilized at the solid surface. Instead, they are similar to the water molecules interacting with OH groups attached to flexible polymer chains. Molecular dynamics (MD) simulations of the polarization behavior of water near the (110) facet of cellulose Iβ crystals confirmed that the interactions of water molecules with the cellulose crystal surface do not cause slower relaxations in the frequency range studied via the DRS. These results indicated that the CNC surface cannot be depicted with the crystallographic facets of cellulose Iβ; instead, it resembles a polymer-brush surface on which the short glucan residues or fragments of the strong acid-catalyzed hydrolysis process are swollen and extended into the aqueous phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enhanced liquid metal wetting on oxide surfaces via patterned particles

Liquid metals typically do not intrinsically wet oxide surfaces due to chemical dissimilarity, preventing the formation of well-bonded interfaces in many applications. This paper investigates how the wetting and spreading of liquid silver (Ag) droplets on oxide surfaces can be enhanced via patterned nickel (Ni) particles using multiscale modeling. First, a force field (FF) for an Ag/yttria-stabilized-zirconia system was developed by fitting the interfacial binding energy to density-functional theory (DFT) calculations. Then, molecular dynamics (MD) simulations using these FF parameters revealed that the good wetting of Ag on Ni caused Ag droplets to be attracted to the Ni particles, promoting spreading. Thus, the optimal hexagonal Ni pattern was shown to enhance the wetting area by 224% compared to a bare oxide surface. Second, an analytical model predicting the maximum possible enhancement in wetting area as a function of the intrinsic wetting angle, size of the Ag droplet, and size of the Ni particles was derived and validated with MD simulations, to generalize this method to other substrates. Here, both the analytical model and MD results suggested that the ability of Ni patterns to enhance the Ag wetting area becomes more significant as the wettability of substrate becomes worse.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhanced resistance to helium irradiations through unusual interaction between high-entropy-alloy and helium

Finding high performance plasma-facing materials (PFMs) is one of the most important and challenging tasks for realizing the commercial application of fusion reactors. Herein, we found the CrMoTaWV high entropy alloy (HEA) is highly resistant to low-energy and high-flux He plasma exposure. The nanochannel HEA film has 20 times higher initial fluence for the formation of fuzz and a remarkable 8.9 times slower fuzz growth rate than those of W. Combining the in-situ TEM observation and the Molecular dynamics (MD) simulation of the He bubble growth process, a new mechanism for the enhanced radiation resistance in HEA with the unusual interaction between HEA and He is found, where, differing from traditional metal, bubble growth in HEA leads to non-directional emission of interstitial atoms while HEA greatly suppress the growth of He bubbles. Additionally, the special nanochannel structure further rise the radiation resistance through releasing He out of the HEA film and reducing the He concentration. This new nanochannel refractory HEA material presents a promising choice as the PFMs with excellent performance and a much longer serving lifetime for future commercial fusion reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transition in helium bubble strengthening of copper from quasi-static to dynamic deformation

Damage from low-temperature irradiation and the subsequent degradation of materials performance pose significant challenges for the storage of radioactive materials and for peripheral components in some nuclear reactor designs. Fully understanding the mechanical behavior of such materials requires test data for strain rates in both the quasi-static (< 10/s) and dynamic (>> 10/s) regimes. While dynamic testing has generally been avoided in the past for neutron irradiated (contamination concerns) and ion irradiated (insufficient volume) materials, surface-sensitive Richtmyer-Meshkov instability (RMI) tests were used in the present work to overcome these limitations. Here, nanopillar compression, nanoindentation, and RMI testing data from a helium implanted surface layer (~10 µm thick) were compiled to explore the effects of helium bubbles on the materials strength of high-purity copper at strain rates of 0.001/s – 10 8 /s. While nano-mechanical testing revealed increases in yield strength and hardness with increasing helium dose from 1000 to 4000 appm He, RMI indicated no significant changes in strength as compared to unimplanted copper. Here, this discrepancy in behavior was rationalized through a combination of recent literature and follow-on molecular dynamics (MD) simulations, leading to the conclusion that the nanoscale helium bubbles acting as dispersed barriers to dislocation motion at quasi-static strain rates collapse under shock loading and cease to be effective barriers at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2024

High-pressure phase transitions in a laser directed energy deposited Fe-33Cu Alloy

Here, additively manufactured Fe-Cu alloys contain both equilibrium face-centered cubic (FCC) and metastable body-centered cubic (BCC) crystal structure Cu precipitates depending on their size. However, the stability of these nanoscale precipitates under extreme conditions such as high pressures has not been reported. This study investigates the phase transformations and microstructural stability of laser directed energy deposition (DED-LB) made Fe67Cu33 alloy (nominal composition in at.%) under high static pressure deformation using in-situ synchrotron X-ray diffraction under pressure, postmortem high-resolution scanning transmission electron microscopy (HR-STEM), and molecular dynamics (MD) simulations. In-situ XRD results reveal a reversible phase transformation from BCC to the hexagonal close-packed (HCP) structure in the Fe grains at an onset pressure of 16.4 GPa, significantly higher than reported in the literature for pure Fe. Although no high-pressure phase transition was observed in the FCC Cu grains through XRD, HR-STEM analysis uncovers a phase transition to the HCP structure in nanoscale metastable BCC Cu precipitates within the BCC Fe matrix. After decompression, the Fe matrix reverted back to the BCC structure with periodic lath martensite, while regions of the nanoscale BCC Cu precipitates retained the metastable HCP structure. MD simulations support the BCC → HCP transition in the embedded nanoscale coherent Cu precipitates, consistent with the classical Burgers mechanism. Thus, by leveraging in-situ XRD observation, postmortem high-resolution S/TEM microscopy, and MD simulations, this study offers profound insights into the distinctive phase transformations induced by high pressure in DED-LB Fe67Cu33 alloy, which is distinguished by its hierarchical microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2020

The carboxylation status of osteocalcin has important consequences for its structure and dynamics

The carboxylation status of Osteocalcin (Ocn) not only influences formation and structure in bones but also has important endocrine functions affecting energy metabolism and expenditure. In this study, the role of γ-carboxylation of the glutamate residues in the structure-dynamics-function relationship in Ocn is investigated. Three forms of Ocn, differentially carboxylated at the Glu-17, 21 and 24 residues, along with a mutated form of Ocn carrying Glu/Ala mutations, are modeled and simulated using molecular dynamics (MD) simulation in the presence of calcium ions. Characterization of the global conformational dynamics of Ocn, described in terms of the orientational variations within its 3-helical domain, highlights large structural variations in the non-carboxylated osteocalcin (nOcn). The bi-carboxylated Ocn (bOcn) and tri-carboxylated (tOcn) species, in contrast, display relatively rigid tertiary structures, with the dynamics of most regions strongly correlated. Radial distribution functions calculated for both bOcn and tOcn show long-range ordering of the calcium ion distribution around the carboxylated glutamate (γGlu) residues, likely playing an important role in promoting stability of these Ocns. Additionally, the same calcium ions are observed to coordinate with neighboring γGlu, better shielding their negative charges and in turn stabilizing these systems more than do the singly coordinating calcium ions observed in the case of nOcn. bOcn is also found to exhibit a more helical C-terminal structure, that has been shown to activate its cellular receptor GPRC6A, highlighting the allosteric role of Ocn carboxylation in modulating the stability and binding potential of the active C-terminal. Furthermore, the carboxylation status of Ocn as well and its calcium coordination appear to have a direct influence on Ocn structure and dynamics, possibly leading to the known differences in Ocn biological function.

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

Unlocking the secret of lignin-enzyme interactions: Recent advances in developing state-of-the-art analytical techniques

Bioconversion of renewable lignocellulosics to produce liquid fuels and chemicals is one of the most effective ways to solve the problem of fossil resource shortage, energy security, and environmental challenges. Among the many biorefinery pathways, hydrolysis of lignocellulosics to fermentable monosaccharides by cellulase is arguably the most critical step of lignocellulose bioconversion. In the process of enzymatic hydrolysis, the direct physical contact between enzymes and cellulose is an essential prerequisite for the hydrolysis to occur. However, lignin is considered one of the most recalcitrant factors hindering the accessibility of cellulose by binding to cellulase unproductively, which reduces the saccharification rate and yield of sugars. This results in high costs for the saccharification of carbohydrates. The various interactions between enzymes and lignin have been explored from different perspectives in literature, and a basic lignin inhibition mechanism has been proposed. However, the exact interaction between lignin and enzyme as well as the recently reported promotion of some types of lignin on enzymatic hydrolysis is still unclear at the molecular level. Multiple analytical techniques have been developed, and fully unlocking the secret of lignin-enzyme interactions would require a continuous improvement of the currently available analytical techniques. This review summarizes the current commonly used advanced research analytical techniques for investigating the interaction between lignin and enzyme, including quartz crystal microbalance with dissipation (QCM-D), surface plasmon resonance (SPR), attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy, atomic force microscopy (AFM), nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy (FLS), and molecular dynamics (MD) simulations. Interdisciplinary integration of these analytical methods is pursued to provide new insight into the interactions between lignin and enzymes. Finally, this review will serve as a resource for future research seeking to develop new methodologies for a better understanding of the basic mechanism of lignin-enzyme binding during the critical hydrolysis process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Molecular simulation and artificial intelligence for the circular economy of bioenergy and bioproducts

The concept of the circular bioeconomy is a carbon neutral, sustainable system with zero waste. One vision for such an economy is based upon lignocellulosic biomass. This lignocellulosic circular bioeconomy requires CO 2 absorption from biomass growth and the efficient deconstruction of recalcitrant biomass into solubilized and fractionated biopolymers, which are then used as precursors for the sustainable production of high-quality liquid fuels, chemical bioproducts, and bio-based materials. Here, in this study, we summarize the roles that molecular dynamics (MD) simulations and machine learning (ML) are playing in overcoming several fundamental challenges hindering the adoption of a circular bioeconomy. Specifically, we discuss the role of MD and ML/AI in overcoming lignocellulose recalcitrance by designing biomass pretreatment methods to efficiently produce solubilized cellulose/lignin/hemicellulose and of that in improving energy-intensive manufacturing of biomass-based materials and their structural and mechanical properties. Quantum mechanical methods and MD simulations, in addition to offering a mechanistic understanding of biomass deconstruction and biomaterials design, can provide meaningful structural, energetics, and physiochemical properties as inputs to train AI/ML models. The ML models can guide the experimental prioritization of materials/solvents and process parameters that significantly accelerate the development of biofuel and biomaterial components of the circular bioeconomy.

Smith, Jeremy C. [Oak Ridge National Laboratory (O↗
DOE OSTI2024

Computationally guided synthesis of carbon coated mesoporous silica materials

Mesoporous silica materials (MSMs) are well-suited for biomedical applications due to their unique features, including a large surface area and tunable pore size. To enhance their durability, the small pores in MSMs are filled with carbon precursors and then carbonized to prevent them from interacting with unreacted silicic acid. Here, in this study, we synthesized and healed MSMs using a combination of non-reactive and reactive molecular dynamics (MD) simulations. The non-reactive MD simulation revealed that the self-assembly of Pluronic® L64 polymers in water resulted in nearly 80 % hydrogen bonds between the hydrophilic sections of the micelle and water. In the bond-boosted ReaxFF MD simulations, silicic acid precursors were condensed on the micelle surface, with over 60 % of them leading to the creation of periodic mesoporous silica within the system. Condensation of silicic acid precursors at 300 K with bond-boosting and at 1500 K without it both significantly promoted the polymerization of Si(OH) 4 , with the latter doubling the rate compared to the former. Subsequently, we healed the MSM surface by carbonizing carbon precursors inside an MSM pore. Polyethylene (PE) and high-rank lignite were identified as the most suitable precursors due to their ability to form turbostratic graphene structures. High-rank lignite exhibited the highest carbon conversion ratio to 6-membered rings, closely followed by PE, in the carbon ring formation analysis. Additionally, the production of gases, such as H 2 , increased significantly for PE at both 2200 K and 2600 K, indicating the conversion of a considerable portion of carbon into graphitic or turbostratic structures. The carbonization of PE primarily led to the formation of planar (sp 2 ) structures, while sucrose yielded the least planar structures. Finally, we studied the protective blocking of unreacted silicic acid precursor by considering a PET turbostratic graphene structure in a silica mesopore formed at 2600 K. The trajectory analysis showed that the surface of the silica was effectively coated with PET tar, preventing unreacted silicic acid from interacting with the inner silica pore surface. These findings offer valuable insights into the synthesis and carbonization-based healing processes of MSMs, enhancing their potential for various biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

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

On the integration of molecular dynamics, data science, and experiments for studying solvent effects on catalysis

Computational workflows that combine molecular dynamics (MD) simulations and emerging data-centric (DC) methods can accelerate the screening and analysis of solvent systems experimentally and computationally. Here, MD simulations provide atomic positions and velocities of reactant, solvent, and catalyst materials that can be manipulated into data representations that in turn can be used by DC techniques to conduct predictive modeling, feature extraction, and experimental design. For liquid-phase catalytic applications, emerging DC techniques such as Convolutional and Graph Neural Networks (CNN/GNN), Topological Data Analysis (TDA), and Active Learning (AL) can leverage MD and experimental data to quickly predict solvent effects on reaction outcomes. For instance, in recent studies, 3D solvent environments obtained with MD have been exploited by CNNs to predict experimental reaction rates for homogeneous acid-catalyzed lignocellulosic processes. In this perspective, we discuss basic principles of DC methods and how these can be combined with MD to enable high-throughput screening of solvent selection for diverse catalysis applications.

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

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗