Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular processes”

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 OSTI2022

Molecular-orbital-based machine learning for open-shell and multi-reference systems with kernel addition Gaussian process regression

We introduce a novel machine learning strategy, kernel addition Gaussian process regression (KA-GPR), in molecular-orbital-based machine learning (MOB-ML) to learn the total correlation energies of general electronic structure theories for closed- and open-shell systems by introducing a machine learning strategy. The learning efficiency of MOB-ML(KA-GPR) is the same as the original MOB-ML method for the smallest criegee molecule, which is a closed-shell molecule with multi-reference characters. In addition, the prediction accuracies of different small free radicals could reach the chemical accuracy of 1 kcal/mol by training on one example structure. Accurate potential energy surfaces for the H10 chain (closed-shell) and water OH bond dissociation (open-shell) could also be generated by MOB-ML(KA-GPR). To explore the breadth of chemical systems that KA-GPR can describe, we further apply MOB-ML to accurately predict the large benchmark datasets for closed- (QM9, QM7b-T, and GDB-13-T) and open-shell (QMSpin) molecules.

Chemistry↗
DOE OSTI2021

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 μs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

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

Analysis of vibronic coupling in a 4f molecular magnet with FIRMS

Vibronic coupling, the interaction between molecular vibrations and electronic states, is a fundamental effect that profoundly affects chemical processes. In the case of molecular magnetic materials, vibronic, or spin-phonon, coupling leads to magnetic relaxation, which equates to loss of magnetic memory and loss of phase coherence in molecular magnets and qubits, respectively. The study of vibronic coupling is challenging, and most experimental evidence is indirect. Here we employ far-infrared magnetospectroscopy to directly probe vibronic transitions in [Yb(trensal)] (where H 3 trensal = 2,2,2-tris(salicylideneimino)trimethylamine). We find intense signals near electronic states, which we show arise due to an “envelope effect” in the vibronic coupling Hamiltonian, which we calculate fully ab initio to simulate the spectra. We subsequently show that vibronic coupling is strongest for vibrational modes that simultaneously distort the first coordination sphere and break the C 3 symmetry of the molecule. With this knowledge, vibrational modes could be identified and engineered to shift their energy towards or away from particular electronic states to alter their impact. Hence, these findings provide new insights towards developing general guidelines for the control of vibronic coupling in molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Out of the Dark and Into the Light: A New View of Phytochrome Photobodies

Light is a critical environmental stimulus for plants, serving as an energy source via photosynthesis and a signal for developmental programming. Plants perceive light through various light-responsive proteins, termed photoreceptors. Phytochromes are red-light photoreceptors that are highly conserved across kingdoms. In the model plant Arabidopsis thaliana , phytochrome B serves as a light and thermal sensor, mediating physiological processes such as seedling germination and establishment, hypocotyl growth, chlorophyll biogenesis, and flowering. In response to red light, phytochromes convert to a biologically active form, translocating from the cytoplasm into the nucleus and further compartmentalizes into subnuclear compartments termed photobodies. PhyB photobodies regulate phytochrome-mediated signaling and physiological outputs. However, photobody function, composition, and biogenesis remain undefined since their discovery. Based on photobody cellular dynamics and the properties of internal components, photobodies have been suggested to undergo liquid-liquid phase separation, a process by which some membraneless compartments form. Here, we explore photobodies as environmental sensors, examine the role of their protein constituents, and outline the biophysical perspective that photobodies may be undergoing liquid-liquid phase separation. Understanding the molecular, cellular, and biophysical processes that shape how plants perceive light will help in engineering improved sunlight capture and fitness of important crops.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗
DOE OSTI2024

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Environmental Molecular Sciences Laboratory (EMSL) 2020 Fact Sheet

EMSL, the Environmental Molecular Sciences Laboratory, is a U.S. Department of Energy (DOE) Office of Science User Facility supported by the DOE Office of Biological and Environmental Research (BER). Research at EMSL supports BER’s mission to understand the biological, biogeochemical, and physical principles needed to predict a continuum of processes occurring at the molecular and genomics-controlled smallest scales to environmental and Earth system change at the largest scales. EMSL is in Richland, WA, at Pacific Northwest National Laboratory.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Population balance models for polymer upcycling: signatures of the mechanism in the molecular weight evolution

Chemical and catalytic upcycling processes could help realize a circular plastics economy, but current models for testing mechanistic hypotheses and designing catalysts remain primitive. This work shows how proposed catalytic mechanisms can be incorporated into population balance models to predict the time evolution of molecular weight distributions. We develop models for homogeneous and heterogeneous catalysts, including catalysts that cut at chain ends and catalysts that cut at random locations. For heterogeneous catalysts, we illustrate the effect of adsorption constants that depend on polymer chain length. Here, we discuss ongoing efforts and challenges in measuring and modeling the time evolving molecular weight distributions in polymer upcycling processes.

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

Processes in DNA damage response from a whole-cell multi-omics perspective

Technological advances have made it feasible to collect multi-condition multi-omic time courses of cellular response to perturbation, but the complexity of these datasets impedes discovery due to challenges in data management, analysis, visualization, and interpretation. Here, we report a whole-cell mechanistic analysis of HL-60 cellular response to bendamustine. We integrate both enrichment and network analysis to show the progression of DNA damage and programmed cell death over time in molecular, pathway, and process-level detail using an interactive analysis framework for multi-omics data. Our framework, Mechanism of Action Generator Involving Network analysis (MAGINE), automates network construction and enrichment analysis across multiple samples and platforms, which can be integrated into our annotated gene-set network to combine the strengths of networks and ontology-driven analysis. Taken together, our work demonstrates how multi-omics integration can be used to explore signaling processes at various resolutions and demonstrates multi-pathway involvement beyond the canonical bendamustine mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Cryo-EM Structures of Human Drosha and DGCR8 in Complex with Primary MicroRNA

Metazoan microRNAs require specificmaturation steps initiated by Microprocessor, comprising Drosha and DGCR8. Lack of structural information for the assembled complex has hindered an understanding of how Microprocessor recognizes primary microRNA transcripts (pri-miRNAs). In this paper we present a cryoelectron microscopy structure of human Microprocessor with a pri-miRNA docked in the active site, poised for cleavage. The basal junction is recognized by a four-way intramolecular junction in Drosha, triggered by the Belt and Wedge regions that clamp over the ssRNA. The belt is important for efficiency and accuracy of pri-miRNA processing. TwodsRBDs form a molecular ruler tomeasure the stem length between the twodsRNA-ssRNAjunctions. The specific organization of thedsRBDs near the apical junction is independent of Drosha core domains, as observed in a second structure in the partially docked state. Collectively, we derive a molecular model to explain how Microprocessor recognizes a pri-miRNA and accurately identifies the cleavage site.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

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

Modeling the sequential dissociative double ionization of O 2 by ultrashort intense infrared laser pulses

A density matrix approach for sequential double ionization (DM-SDI) of molecules has been developed recently and was applied to the N 2 molecule. In this article, we extended the DM-SDI model to O 2 , which is a more complicated system to model than N 2 , due to its electronic structures and spin-orbit and laser couplings in the manifold of doubly charged states. We obtained a good agreement on the kinetic energy release spectrum of O + + O + from previous experiments. Thanks to the low computational cost of the model, we explored the mechanism behind the ionization and dissociation dynamics as well as the effects of lasers on the spectrum. Furthermore, this work will pave the way to model sequential dissociative double ionization of larger molecules and to probe molecular dynamics by measuring kinetic energy release spectra from this process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Fast neutrino cooling of nuclear pasta in neutron stars: Molecular dynamics simulations

The direct Urca process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crusts of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct Urca reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct Urca reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified Urca process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Molecular Specificity and Proton Transfer Mechanisms in Aerosol Prenucleation Clusters Relevant to New Particle Formation

Atmospheric aerosol particles influence Earth’s radiative energy balance and cloud properties either directly through interacting with the solar radiation or indirectly through acting as cloud condensation nuclei (CCN), thus impacting Earth’s climate change, air quality, and human health. Regarding both the important scientific and overarching practical implications of aerosols, the past two decades have seen an extensive research effort, with emphasis on the chemical compositions and underlying mechanisms of aerosol formation. New particle formation (NPF) has been recognized as an important pathway, contributing up to 50% of atmospheric aerosols. Nowadays, the general consensus in NPF is that it proceeds via two distinct stages—the nucleation from gaseous precursors to form critical nuclei of sub-1–2 nanometer size, and the subsequent growth into large particles. However, a fundamental understanding of both the NPF process and the molecular-level characterization of the critical size aerosol clusters is still largely missing, hampering the development of reliable and predictive aerosol nucleation and climate models.

Hou, Gao-Lei↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultrafast infrared nano-imaging of far-from-equilibrium carrier and vibrational dynamics

Ultrafast infrared nano-imaging has demonstrated access to ultrafast carrier dynamics on the nanoscale in semiconductor, correlated-electron, or polaritonic materials. However, mostly limited to short-lived transient states, the contrast obtained has remained insufficient to probe important long-lived excitations, which arise from many-body interactions induced by strong perturbation among carriers, lattice phonons, or molecular vibrations. Here, we demonstrate ultrafast infrared nano-imaging based on excitation modulation and sideband detection to characterize electron and vibration dynamics with nano- to micro-second lifetimes. As an exemplary application to quantum materials, in phase-resolved ultrafast nano-imaging of the photoinduced insulator-to-metal transition in vanadium dioxide, a distinct transient nano-domain behavior is quantified. In another application to lead halide perovskites, transient vibrational nano-FTIR spatially resolves the excited-state polaron-cation coupling underlying the photovoltaic response. These examples show how heterodyne pump-probe nano-spectroscopy with low-repetition excitation extends ultrafast infrared nano-imaging to probe elementary processes in quantum and molecular materials in space and time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗