Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular data”

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 559 records · Page 31

Comparative transcriptomics uncovers poplar and fungal genetic determinants of ectomycorrhizal compatibility

Ectomycorrhizal symbiosis supports tree growth and is crucial for nutrient cycling and temperate and boreal ecosystems functioning. The establishment of functional ectomycorrhiza (ECM) first requires the association of compatible partners. However, host and fungal genetic determinants governing mycorrhizal compatibility are unknown. To identify such factors in poplar and its fungal associates, we mined existing and de novo tree and fungal transcriptional datasets. We identified co-expressed genes enabling ECM symbiosis at early and mature stages of the interaction. These sets of genes can be divided into general fungal-sensing and ECM-specific components. We highlight the importance of fungal modulation of plant JA-related defenses and the regulation of secretory pathways for ECM compatibility, including upregulation of key fungal small secreted proteins, the downregulation of plant secreted peroxidases, and the downregulation of plant cell wall remodeling proteins concomitantly with the upregulation of fungal glycosyl hydrolases acting on pectin. Not only gene regulation, but also its temporal scale and dynamics seem to play a crucial role for mycorrhizal compatibility. The expression profile of the host Common Symbiosis Pathway and nutrient transporters was also studied, revealing constitutive levels of expression and moderate upregulation in compatible ECM interactions. Overall, these results underscore the importance of novel biological functions during the establishment of ECM symbiosis, help us gain insights into the molecular events determining mycorrhiza compatibility, and serve as a data-rich transcriptomic resource to open new research questions in the field.

Marqués‐Gálvez, José Eduardo↗

Benchmarking the performance of uncertainty quantification methods for neural network-based interatomic potentials

Machine-learned interatomic potentials (ML-IAPs) continue to gain popularity as accurate, computationally efficient replacements for traditional, physics-based interatomic potentials and expensive ab initio methods. Uncertainty quantification (UQ) of ML-IAPs is a growing area of research as UQ is critical in many applications of IAPs, such as developing curated datasets, active learning-based data augmentation, self-improving models, and estimating the uncertainty of molecular dynamics simulations. In this paper, we construct and benchmark a series of different neural network potentials (NNPs) with varying network architectures to determine the performance of these models with respect to both the mean and uncertainty calibration error. Each NNP method is specifically designed to predict either epistemic or aleatoric uncertainty with particular focus on the differences in behavior between the epistemic and aleatoric uncertainty estimates. We benchmark these methods using multiple datasets common in the ML-IAP literature. The results show that the aleatoric uncertainty from single-shot model architectures is a competitive alternative to ensemble-based epistemic uncertainty predictions in regions of sufficient data-density. However, in regions where the representative data is sparse, aleatoric uncertainty models tend to overpredict and epistemic methods tend to underpredict the actual model error. We conclude that the type of UQ is crucial when discussing performance of probabilistic model results as different methods have different performance characteristics depending on the regime in which they are evaluated. Therefore, the type of UQ method should be carefully evaluated against both the data characteristics and requirements for the intended application.

97 MATHEMATICS AND COMPUTING↗

Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces

An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.

36 MATERIALS SCIENCE↗

Chemical Reactivity Through Adaptive Quantum Mechanics/Many-Body Representations: Theoretical Development, Software Implementation, and Applications (Final Report)

The main objective of this research project was the development and application of a new theoretical/computational framework to model chemical transformations and electronic excitations in fluid mixtures across different phases. Our theoretical/computational framework combines our data-driven many-body (DD-MB) potentials representing molecular interactions with adaptive schemes for modeling chemical reactions in solution. The combination of these two components resulted in an adaptive quantum mechanics/many-body (adQM/MB) method that largely suppresses discontinuities between QM and MM regions, which affect existing QM/MM methods, and thus provides an accurate representation of both quantum mechanical and environmental effects through a rigorous description of mutual polarization between QM and MM regions. The implementation of our DD-MB potentials and adQM/MB method in popular software enabled computer simulations of solvation phenomena, reactive processes, and electronic excitations in fluid mixtures with chemical and spectroscopic accuracy, representing a major step toward realistic computer simulations of a wide range of molecular systems relevant to the DOE mission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 4.0 (MSTDB-TP V.4.0)

A central repository of thermophysical and thermochemical properties of molten salt compositions of relevance to molten salt reactors (MSRs) is vital in supporting the broad community of MSR developers, who are at various stages of developing and deploying their reactor designs. In general, these MSR designs differ significantly from developer to developer (e.g., with respect to the hardness of the neutron spectra, level of fissile loading, target multicomponent temperatures and power levels, and moderating capabilities). Therefore, the fuel and coolant salts being considered vary greatly: they may be chlorides or fluorides, they utilize different actinides at different ratios, and the cations in the melt are selected based on perceived advantages and disadvantages. Considering the general need for thermal properties, and the vastness of the array of potential candidate salt mixtures, the Molten Salt Thermal Properties Database (MSTDB) was initiated in 2018 with the goal of providing thermophysical and thermochemical characterization of key molten salt compounds and mixtures across their temperature and compositional domains. The MSTDB is thus divided into the thermophysical arm (MSTDB-TP) and the thermochemical arm (MSTDB-TC). The MSTDB is an effort funded by the Department of Energy, Office of Nuclear Energy (DOE-NE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, and the MSR Campaign. This report provides an overview of the MSTDB-TP v4.0 in terms of the data contained within, the state of the tools used to access the data, the availability of predictive models that leverage the raw data in the database, the preliminary status of developmental efforts that are currently underway, and an account of future goals for MSTDB-TP. The primary goal for the update from MSTDB-TP v.3.1 to v4.0 was the incorporation of surface tension data into the database; this property is important for thermal hydraulics modeling and species transport in other tools that have been developed under the NEAMS program. A breakdown of the surface tension data that have been added into MSTDB-TP v4.0 is provided herein, and the manner in which the quality of the data has been assessed is also documented. For MSTDB-TP v4.0, newly published thermophysical property data—primarily from collaborative experimental efforts under the MSR Campaign—have been incorporated into the database, and the resulting expansion is documented here. Because of the size to which MSTDB-TP has grown, the raw data format has now been recast into JavaScript Object Notation (JSON) format for easier connection with the MSTDB-TP application programming interface (API). Saline; the pre-existing comma-separated value (CSV) format has been deprecated but is still maintained, accessible, and up to date. As a final effort in packaging the MSTDB-TP v4.0 update, the graphical user interface (GUI) for MSTDB has been updated to allow full accessibility to the density and viscosity predictive models, which are based on Redlich-Kister expansions of MSTDB-TP raw data. Some other major aspects of this report, in terms of preliminary and future work, include: (1) documentation of the formalism and preliminary testing of a kinetic theory model that may act as a predictive model for thermal conductivity; (2) documentation of the candidate predictive models that may be considered in the future for surface tension, making use of the surface tension data now in MSTDB-TP v4.0; (3) a preliminary account of a data collection process that will enable the filling of additional gaps within MSTDB-TP, namely with data which have been collected computationally (e.g., through ab initio molecular dynamics).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Developing predictive models for µ opioid receptor binding using machine learning and deep learning techniques

Opioids exert their analgesic effect by binding to the µ opioid receptor (MOR), which initiates a downstream signaling pathway, eventually inhibiting pain transmission in the spinal cord. However, current opioids are addictive, often leading to overdose contributing to the opioid crisis in the United States. Therefore, understanding the structure-activity relationship between MOR and its ligands is essential for predicting MOR binding of chemicals, which could assist in the development of non-addictive or less-addictive opioid analgesics. This study aimed to develop machine learning and deep learning models for predicting MOR binding activity of chemicals. Chemicals with MOR binding activity data were first curated from public databases and the literature. Molecular descriptors of the curated chemicals were calculated using software Mold2. The chemicals were then split into training and external validation datasets. Random forest, k-nearest neighbors, support vector machine, multi-layer perceptron, and long short-term memory models were developed and evaluated using 5-fold cross-validations and external validations, resulting in Matthews correlation coefficients of 0.528–0.654 and 0.408, respectively. Furthermore, prediction confidence and applicability domain analyses highlighted their importance to the models’ applicability. Our results suggest that the developed models could be useful for identifying MOR binders, potentially aiding in the development of non-addictive or less-addictive drugs targeting MOR.

Research & Experimental Medicine↗

Intergranular Spin Dependent Tunneling Dominated Magnetoresistance in Helimagnetic Manganese Phosphide Thin Films

Helical magnets are emerging as a novel class of materials for spintronics and sensor applications; however, research on their charge- and spin-transport properties in a thin film form is less explored. Herein, we report the temperature and magnetic field-dependent charge transport properties of a highly crystalline MnP nanorod thin film over a wide temperature range (2 K < T < 350 K). The MnP nanorod films of ~100 nm thickness were grown on Si substrates at 500 °C using molecular beam epitaxy. The temperature-dependent resistivity ρ(T) data exhibit a metallic behavior (dρ/dT > 0) over the entire measured temperature range. However, large negative magnetoresistance (Δρ/ρ) of up to 12% is observed below ~50 K at which the system enters a stable helical (screw) magnetic state. In this temperature regime, the Δρ(H)/ρ(0) dependence also shows a magnetic field-manipulated CONE + FAN phase coexistence. The observed magnetoresistance is dominantly governed by the intergranular spin dependent tunneling mechanism. These findings pinpoint a correlation between the transport and magnetism in this helimagnetic system.

36 MATERIALS SCIENCE↗

Validated ligand geometries for macromolecular refinement restraints and molecular-mechanics force fields

In macromolecular structure refinement, the low observation-to-parameter ratio and the lack of high-resolution data are countered by using a priori information in the form of restraints. Having accurate geometries of the chemical entities in the sample is paramount for generating accurate chemical restraints and, therefore, accurate macromolecular structures. In particular, it is desirable to have accurate restraints for known and novel ligand entities. Quantum mechanics (QM) can minimize the energy of a ligand by adjusting its geometry, and these geometries can be used to generate restraints for macromolecular refinement. This article describes a library of approximately 37 000 small molecules extracted from the Chemical Component Dictionary in the Protein Data Bank and minimized by density-functional QM. The library includes restraint files for use in crystallography or cryo-EM refinement, along with files suitable for molecular-dynamics simulation. Because the geometries are validated using the Cambridge Structural Database, the restraints library provides users with both functional restraints and minimized geometries. This work also provides procedures for generating new and accurate restraints.

Amber↗

Dramatic effect of electrode type on tunnel junction based molecular spintronic devices

A new class of molecular spintronic devices can be fabricated by chemically bonding magnetic molecular channels to the electrodes of a prefabricated tunnel junction with exposed side edges. Prior experimental studies showed that the cyanide-bridged octametallic molecular cluster, [(pzTp)FeIII(CN)3]4[NiII(L)]4¬[O3SCF3]4 [(pzTp) = tetra(pyrazol-1-yl)borate; L = 1-S(acetyl)tris(pyrazolyl)decane] molecule impact depended on the type of metallic electrodes used in the tunnel junction testbed. Experimental magnetization and transport studies showed a dramatic difference in molecule response on tunnel junctions with different combinations of metallic electrodes. Transport via paramagnetic molecular channels on a tunnel junction involving paramagnetic and ferromagnetic metal electrodes was dramatically different than the suppressed current state observed on tunnel junctions involving two ferromagnetic electrodes. We conducted theoretical studies to understand the experimental data and explore a wide range of electrode materials on tunnel junction-based molecular spintronics devices (TJMSD). Here, we report a Monte Carlo simulation study that focuses on understanding the effect of electrodes on the magnetic and physical properties of TJMSD. A 3D Heisenberg model of cross-junction-shaped TJMSD was used for the simulation study. We studied the effects of ferromagnetic, paramagnetic, and antiferromagnetic electrode materials. This study provides insights for designing and understanding futuristic molecular spintronics devices.

36 MATERIALS SCIENCE↗

19-LW-045 Full Length Final Report. Molecular Mechanisms of Bacterial Pathogenesis: Waging the Arms Race with Superbugs

As the current global pandemic makes abundantly clear, we need a better understanding of infectious disease to safeguard human health, the economy and global security. Modern omics techniques hold the promise of providing a comprehensive understanding of the molecular mechanisms of life, including causes of pathogenesis from infectious disease at the molecular level, but we there is a serious gap in annotation of gene function. For as much as half of the genes and gene products encoded in genomes the molecular and/or cellular function is unknown or only partially understood. Recent innovations in fluorescence microscopy for live cell imaging and genetic engineering make it possible to determine the temporal correlation between molecular events, such as a gene being expressed due to host-pathogen interaction, and cellular events, such as bacterial invasion of immune cells. This is turn allows us to gain new insight as to the molecular and cellular role of individual genes and will enable the discovery and validation of new molecular mechanisms essential for infectious disease. Knowing the molecular mechanisms of disease processes will provide new therapeutic targets or novel countermeasure strategies. We aimed to develop a lattice light sheet fluorescence microscope as a unique resource at LLNL for long time course live cell imaging experiments; to develop the reagents and cell lines needed to monitor molecular events during the course pathogenic bacteria infecting mammalian immune cells; and to demonstrate that we could capture molecular events during an infection. We fully commissioned the LLNL lattice light sheet microscope and conducted initial proof of principle imaging experiments on mammalian immune cells and pathogenic bacteria. It is clear from the experience gained that long time course live cell imaging has tremendous potential to help elucidate molecular mechanisms of host-pathogen interactions and to help annotate gene function, which would establish a basis for new countermeasures. It is also clear that if live cell imaging is to realize its full potential new data processing and analysis tools will need to be developed to facilitate analysis of molecular events within cells; new sample chambers and stages could facilitate studies with a wider range of cell and tissue types; and alternative molecular tagging methods need to be explored to enable more facile engineering of cells labeled with molecular specificity.

59 BASIC BIOLOGICAL SCIENCES↗

Examining a hadronic γ-ray scenario for the radiative shell and molecular clouds of the old GeV supernova remnant G298.6−0.0

Abstract Based on the 13.7 yr Fermi-LAT data, Yeung, Bamba, and Sano (2023, PASJ, 75, 384) claimed detection of two γ-ray sources (namely Src-NE and Src-NW) associated with the supernova remnant (SNR) G298.6−0.0, and interpreted it as an old GeV SNR interacting with molecular clouds (MCs). In this follow-up study, we refine the flux measurements below 2 GeV with Fermi-LAT event types of better angular reconstruction. Then, we report our γ-ray spectral fittings and cosmic-ray phenomenology in a hadronic scenario, considering both the shell and MC regions of SNR G298.6−0.0. We confirm that the γ-ray spectra of both Src-NE and Src-NW exhibit spectral breaks, at $1.50_{-0.50}^{+0.60}$ and $0.68_{-0.11}^{+0.32}\:$GeV, respectively. Src-NW has a harder broad-band photon index than Src-NE, suggesting an appreciable difference between the physical separations of their respective emission sites from SNR G298.6−0.0. The cosmic-ray spectrum responsible for Src-NE starts with a minimum energy $E_\mathrm{CR,min}=1.38_{-0.16}^{+0.47}\:$GeV, and has a proton index $\Gamma _\mathrm{CR}=2.57_{-0.21}^{+0.18}$ below the exponential cutoff energy $E_\mathrm{CR,max}=240_{-150}^{+240}\:$GeV. Accordingly, we argue that Src-NE is dominated by the SNR shell, while only a minor portion of lower-energy emission is contributed by the MCs interacting with the SNR. The cosmic-ray population for Src-NW starts at a higher energy such that the ECR, min ratio of Src-NW to Src-NE is ≳2. The high ECR, min, as well as the high cosmic-ray energy density required (∼26 eV cm−3), supports the interpretation that Src-NW is predominantly the γ-ray emission from the farther MCs being bombarded by protons that had earlier escaped from SNR G298.6−0.0. By comparing the high-energy features of G298.6−0.0 with those of analogical SNRs, especially SNR W 28 and SNR W 44, we further constrain the age of SNR G298.6−0.0 to be 10–30 kyr, under the assumption of a purely hadronic scenario.

Astronomy & Astrophysics↗

Ab initio calculations in atoms, molecules, and solids, treating spin–orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin–orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin–orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br 2 and I 2 , and solid Bi.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗

Nanoporous Tio2 Water Training Data

Data and input files used to train a Deep Potential (DP) model for the nanoporous TiO2-water interface. The DeepMD-kit code was used to train the DP. Information about data format, and how to use DeepMD-kit can be found at https://docs.deepmodeling.com/projects/deepmd/en/master/

08 HYDROGEN↗

Non-Adiabatic Excited State Molecular Dynamics Methodologies: comparison and convergence

Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock compression of vanadium at extremes: Theory and experiment

Here, the equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50 000 K, for densities between ρ = 3 and 15 g / cm 3 , with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Statistical Characterization of Temperature and Pressure Vertical Profiles for the Analysis of Laser Heterodyne Radiometry Data

The statistical analysis of historic pressure and temperature profiles from radiosonde launches for use in the fitting of molecular oxygen line shapes is presented. As the O 2 mixing ratio is nearly constant throughout the lower atmosphere, only variations in pressure and temperature profiles will affect the fit of observed O 2 features in Laser Heterodyne Radiometry (LHR) spectra. Radiosonde temperature and pressure data are extracted from the Integrated Global Radiosonde Archive (IGRA) for a given station, date, and launch time. Data may be extracted for a single launch, for the same date over several years, and/or within a window centered on a target date. The temperature and pressure profiles are further characterized by the statistical variation in coefficients of polynomial fits in altitude. The properties of the probability distributions for each coefficient are used to constrain fits of O 2 line shapes through Nelder–Mead optimization. The refined temperature and pressure profiles are then used in the retrieval of vertically resolved mixing ratios for greenhouse gases (GHGs) measured in the same instrument. In continuous collections, each vertical profile determination may be treated as a Bayesian prior to inform subsequent measurements and provide an estimate of uncertainties.

47 OTHER INSTRUMENTATION↗

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗