Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Data Generation for Machine Learning Interatomic Potentials and Beyond

The field of data-driven chemistry is undergoing an evolution, driven by innovations in machine learning models for predicting molecular properties and behavior. Recent strides in ML-based interatomic potentials have paved the way for accurate modeling of diverse chemical and structural properties at the atomic level. The key determinant defining MLIP reliability remains the quality of the training data. A paramount challenge lies in constructing training sets that capture specific domains in the vast chemical and structural space. This Review navigates the intricate landscape of essential components and integrity of training data that ensure the extensibility and transferability of the resulting models. We delve into the details of active learning, discussing its various facets and implementations. We outline different types of uncertainty quantification applied to atomistic data acquisition and the correlations between estimated uncertainty and true error. The role of atomistic data samplers in generating diverse and informative structures is highlighted. Furthermore, we discuss data acquisition via modified and surrogate potential energy surfaces as an innovative approach to diversify training data. The Review also provides a list of publicly available data sets that cover essential domains of chemical space.

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

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2001

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Neutral composition in the thermosphere.

Data obtained from the Ogo 6 neutral atmospheric composition (NAC) experiment for a period of more than 1 year are used to compare the average constituent composition at an altitude of 400 km with that predicted by the Jacchia 1971 and 1965 models. The comparison shows that the Jacchia 1971 model underestimates the molecular nitrogen densities at 400 km by a factor of 2. An atmosphere is constructed down to 120 km by means of the Stein and Walker technique. A fit is made with (1) the 400-km total densities from drag measurements and the composition from Ogo 6 NAC, (2) the 250-km measured molecular nitrogen densities, and (3) the 150-km total densities from drag. This fit shows that the Jacchia 1971 model overestimates the atomic oxygen content at 150 km.

Taeusch, D. R.↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Perturbative theoretical model of electronic transient circular dichroism spectroscopy of molecular aggregates

A chiral analog of transient absorption spectroscopy, transient circular dichroism (TCD) spectroscopy is an emerging time-resolved method. Both spectroscopic methods can probe the electronic transitions of a sample, and TCD is additionally sensitive to the dynamic aspects of chirality, such as those induced by molecular excitons. Here, we develop a theoretical description of TCD for electronic multi-level models in which the pump pulse is linearly polarized and probe pulse is alternately left- and right-circularly polarized. We derive effective response functions analogous to those often used to describe other four-wave mixing methods and then simulate and analyze TCD spectra for three representative multi-level electronic model systems. We elaborate on the presence and detection of the spectral signatures of electronic coherences.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1986

Multistep modeling of protein structure: application to bungarotoxin

Modelling of bungarotoxin in atomic details is presented in this article. The model-building procedure utilizes the low-resolution crystal coordinates of the c-alpha atoms of bungarotoxin, sequence homology within the neurotoxin family, as well as high-resolution x-ray diffraction data of cobratoxin and erabutoxin. Our model-building procedure involves: (a) principles of comparative modelling, (b) embedding procedures of distance geometry, and (c) use of molecular mechanics for optimizing packing. The model is not only consistent with the c-alpha coordinates of crystal structure, but also agrees with solution conformational features of the triple-stranded beta sheet as observed by NOE measurements.

NASA Discipline Exobiology↗
NASA NTRS1986

Translucent molecular clouds: Theory and observations

Few suitable stars behind molecular clouds have been identified. A limited survey was performed of interstellar lines toward highly reddened stars in the southern sky, using the ESO 1.4 m CAT telescope with a Reticon detector, and the Cerro Tololo 4 m telescope equipped with a GEC charge coupled device (CCD) detector. Because of the reduced extinction at longer wavelengths, molecules were searched for with transitions in the red part of the spectrum such as C2 and CN. For some lines-of-sight for which C2 was detected, the 4300 A line of CH was also observed. Absorption lines of interstellar C2 around 8750 A were detected in the spectra of about 1/4 of the 36 observed stars. The inferred C2 column densities range between 10 to the 13th power and 10 to the 14th power sq. cm., and are up to an order of magnitude larger than those found for diffuse clouds. The observed column densities of CH correlate very well with those of C2 over this range. In contrast, the measured column densities of CN vary by orders of magnitude between the various regions, and they do not correlate with those of C2 and CH. The observed rotational population distribution of C2 also provides information about the physical conditions in the clouds. Models of translucent molecular clouds have been constructed along the lines described by van Dishoeck and Black (1986) for diffuse clouds. The models compute accurately the fractions of atomic and molecular hydrogen as functions of depth into the clouds, as well as the excitation of H2 by ultraviolet pumping. They also incorporate a detailed treatment of the photodissociation processes of the molecules (cf. van Dishoeck 1986), which play an important role in the chemistry up to depths of about 3 mag.

Vandishoeck, E. F.↗
NASA NTRS1994

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS1967

Dipole properties of molecular nitrogen.

Optical absorption, refractivity and electron scattering used to construct model dipole spectrum of molecular nitrogen and calculate dipole properties

MOLECULAR GAS↗
DOE OSTI2025

The Birth of Minerals: From Single Step to Multiple Step Mechanisms

Mineral formation from ions in aqueous solutions begins with complex initial stages, where amorphous and liquid-like precursors play pivotal roles before crystalline growth occurs. Both classical and non-classical nucleation and growth theories, introduced in previous chapters, offer explanations, each with their own strengths and limitations, for the complex intermediate phases observed in experimental research. Analytical techniques play a critical role in detecting and characterizing precursor phases, offering valuable insights into nucleation and growth mechanisms across various temporal and spatial scales. In conclusion, molecular dynamics and modelling provide in-depth perspectives on these phases, allowing for a closer examination of their nucleation and growth mechanisms at the molecular level, and revealing the intricate processes that govern their behaviour.

amorphous↗
DOE OSTI2024

Understanding the role of negative charge in the scaffold of an artificial enzyme for CO 2 hydrogenation on catalysis

Here, we have approached the construction of an artificial enzyme by employing a robust protein scaffold, lactococcal multidrug resistance regulator, LmrR, providing a structured secondary and outer coordination spheres around a molecular rhodium complex, [Rh I (P Et2 N gly P Et2 ) 2 ] - . Previously, we demonstrated a 2–3 fold increase in activity for one Rh-LmrR construct by introducing positive charge in the secondary coordination sphere. In this study, a series of variants was made through site-directed mutagenesis where the negative charge is located in the secondary sphere or outer coordination sphere, with additional variants made with increasingly negative charge in the outer coordination sphere while keeping a positive charge in the secondary sphere. Placing a negative charge in the secondary or outer coordination sphere demonstrates decreased activity by a factor of two compared to the wild-type Rh-LmrR. Interestingly, addition of positive charge in the secondary sphere, with the negatively charged outer coordination sphere restores activity. Vibrational and NMR spectroscopy suggest minimal changes to the electronic density at the rhodium center, regardless of inclusion of a negative or positive charge in the secondary sphere, suggesting another mechanism is impacting catalytic activity, explored in the discussion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Photodissociation regions

Photodissociation regions are interstellar regions of predominantly neutral gas where the FUV radiation field plays a significant role in the chemistry and/or the heating. Photodissociation, grain attenuation of the FUV flux and grain photoelectric heating lead to copious emission of C II (158 microns), O I (63 microns), Si II (35 microns), C I (370,609 microns), H2 vibrotational and CO rotational transitions, and IR continuum. Theoretical models compared with observations diagnose such physical parameters as the density and temperature structure, the elemental abundances, and the FUV radiation field. Applications are made to Orion, M17 and galactic nuclei. Theoretical photodissociation models can explain the correlation in the C II (158 microns) and CO J = 1-0 emission and the correlation of the CO J = 1-0 luminosity with the molecular mass. Theoretical models also point to feedback mechanisms which may control the rate of star formation in galaxies and which may regulate the column density through giant molecular clouds.

Hollenbach, David J.↗
NASA NTRS2003

Nanolayered Features of Collagen-like Peptides

We have been investigating collagen-like model oligopeptides as molecular bases for complex ordered biomimetic materials. The collagen-like molecules incorporate aspects of native collagen sequence and secondary structure. Designed modifications to native primary and secondary structure have been incorporated to control the nanostructure and microstructure of the collagen-like materials produced. We find that the collagen-like molecules form a number of lyotropic rod liquid crystalline phases, which because of their strong temperature dependence in the liquid state can also be viewed as solvent intercalated thermotropic liquid crystals. The liquid crystalline phases formed by the molecules can be captured in the solid state by drying off solvent, resulting in solid nanopatterned (chemically and physically) thermally stable (to greater than 100 C) materials. Designed sequences which stabilize smectic phases have allowed a variety of nanoscale multilayered biopolymeric materials to be developed. Preliminary investigations suggest that chemical patterns running perpendicular to the smectic layer plane can be functionalized and used to localize a variety of organic, inorganic, and organometallic moieties in very simple multilayered nanocomposites. The phase behavior of collagen-like oligopeptide materials is described, emphasizing the correlation between mesophase, molecular orientation, and chemical patterning at the microscale and nanoscale. In many cases, the textures observed for smectic and hexatic phase collagens are remarkably similar to the complex (and not fully understood) helicoids observed in biological collagen-based tissues. Comparisons between biological morphologies and collagen model liquid crystalline (and solidified materials) textures may help us understand the molecular features which impart order and function to the extracellular matrix and to collagen-based mineralized tissues. Initial studies have utilized synthetic collagen-like peptides while future work will also focus on similar sequences generated via genetic engineering methods.

Valluzzi, Regina↗
DOE OSTI2025

Protein-ligand binding affinity prediction using multi-instance learning with docking structures

Recent advances in 3D structure-based deep learning approaches demonstrate improved accuracy in predicting protein-ligand binding affinity in drug discovery. These methods complement physics-based computational modeling such as molecular docking for virtual high-throughput screening. Despite recent advances and improved predictive performance, most methods in this category primarily rely on utilizing co-crystal complex structures and experimentally measured binding affinities as both input and output data for model training. Nevertheless, co-crystal complex structures are not readily available and the inaccurate predicted structures from molecular docking can degrade the accuracy of the machine learning methods. We introduce a novel structure-based inference method utilizing multiple molecular docking poses for each complex entity. Our proposed method employs multi-instance learning with an attention network to predict binding affinity from a collection of docking poses. We validate our method using multiple datasets, including PDBbind and compounds targeting the main protease of SARS-CoV-2. The results demonstrate that our method leveraging docking poses is competitive with other state-of-the-art inference models that depend on co-crystal structures. This method offers binding affinity prediction without requiring co-crystal structures, thereby increasing its applicability to protein targets lacking such data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Relationships between the Stellar, Gaseous, and Star Formation Disks in LITTLE THINGS Dwarf Irregular Galaxies: Indirect Evidence for Substantial Fractions of Dark Molecular Gas

The stellar, gaseous and young stellar disks in the LITTLE THINGS sample of nearby dwarf irregular galaxies are fitted with functions to search for correlations between the parameters. We find that the H i radial profiles are generally flatter in the center and fall faster in the outer regions than the V-band profiles, while young stars are more centrally concentrated, especially if the H i is more centrally flat. This pattern suggests that the H i is turning into molecules in the center, and the molecular clouds are forming stars and FUV. A model that assumes the molecular surface density is proportional to the total gas surface density to a power of 1.5 or 2, in analogy with the Kennicutt–Schmidt relation, reproduces the relationship between the ratio of the visible to the H i scale length and the H i Sérsic index. The molecular fraction is estimated as a function of radius for each galaxy by converting the FUV to a molecular surface density using conventional calibrations. The average molecular fraction inside 3R {sub D} is 23% ± 17%. However, the break in the stellar surface brightness profile has no unified tracer related to star formation.

74 ATOMIC AND MOLECULAR PHYSICS↗