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 1,189 records · Page 66

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗

Determining critical overlap concentration of polyethylene oxide to support excipient safety assessment of opioid products

Intravenous administration of abuse-deterrent opioid products poses high safety risks, in part due to the presence of high molecular weight polymeric excipients. Previous in vivo studies in animal models have shown that the higher molecular weight (Mw) polymeric excipients like polyethylene oxide (PEO) were directly linked to such adverse responses as intravenous hemolysis and kidney damage. PEO polymers have been widely used in abuse-deterrent formulations (ADF) of opioid products, adding to concerns over the general safety of the opioid category due to the unknown safety risk from abuse via unintended routes. Here, the current study focused on the determination of the critical overlap concentration (c*) at various PEO molecular weights to aid in explaining differences in observed adverse responses from previous animal studies on the intravenous administration of PEO solutions. Adverse in vivo responses may be related to the viscoelastic regime of the polymer solution, which depends not only on Mw but also on concentration. Having a localized polymer concentration in the blood above the c*, i.e., the transition from the dilute to semi-dilute entangled viscoelastic regime, may influence the flow behavior and interactions of cells in the blood. The relationship of c* to this combination of physical, chemical, and rheological effects is a possible driving force behind adverse in vivo responses.

60 APPLIED LIFE SCIENCES↗

Low-temperature formation of pyridine and (iso)quinoline via neutral–neutral reactions

Aromatic molecules represent fundamental building blocks in prebiotic chemistry and are contemplated as vital precursors to DNA and RNA nitrogen bases. However, despite the identification of some 300 molecules in extraterrestrial environments, the pathways to pyridine (C 5 H 5 N), pyridinyl (C 5 H 4 N•), and (iso)quinoline (C 9 H 7 N) – the simplest represen-tative of mono and bicyclic aromatic molecule carrying nitrogen – are elusive. Here, we afford compelling evidence on the gas-phase formation of methylene amidogen (H 2 CN•) and cyanomethyl (H 2 CCN•) radicals via molecular beam studies and electronic structure calcula-tions. The modeling of the chemistries of Taurus Molecular Cloud (TMC-1) and Titan’s atmosphere contemplates a complex chain of reactions synthesizing pyridine, pyridinyl, and (iso)quinoline from H 2 CN• and H 2 CCN• at levels of up to 75%. Furthermore, this study affords unique entry points to precursors of DNA and RNA nitrogen bases in hydrocarbon-rich extraterrestrial environments thus changing the way we think about the origin of prebiotic molecules in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Solvent effects on triplet yields in BODIPY-based photosensitizers

Here, we employ molecular dynamics simulations and quantum rate theories to elucidate the complex condensed-phase dynamics underpinning triplet-state formation in organic photosensitizers. Using models informed by first-principles calculations complete with a molecular representation of solvents of different polarities, we elucidate the interplay of the internal and environmental interactions underlying triplet yield. We find that triplet yields depend sensitively on the dielectric stabilization of the charge transfer intermediate that facilitates a transition into the triplet manifold. Our results illustrate the importance of molecularly detailed models in understanding the excited-state internal charge-transfer dynamics of photochemically relevant organic molecules.

Coello Escalante, Leonardo [University of Californ↗

Photochemical Modeling of CH3 Abundances in the Outer Solar System

Recent measurements of methyl radicals (CH3) in the upper atmospheres of Saturn and Neptune by the Infrared Space Observatory (ISO) provide new constraints to photochemical models of hydrocarbon chemistry in the outer solar system. The derived column abundances of CH3 on Saturn above 10 mbar and Neptune above the 0.2 mbar pressure level are (2.5 - 6.0) x 10(exp 13) / sq cm and (0.7 - 2.8) x 10(exp 13) / sq cm, respectively. We use the updated Caltech/Jet Propulsion Laboratory photochemical model, which incorporates hydrocarbon photochemistry, vertical molecular and bulk atmospheric eddy diffusion, and realistic radiative transfer modeling, to study the CH3 abundances in the upper atmosphere of the giant planets and Titan. We identify the key reactions that control the concentrations of CH3 in the model, such as the three-body recombination reaction, CH3 + CH3 + M yields C2H6 + M. We evaluate and extrapolate the three-body rate constant of this reaction to the low-temperature limit (1.8 x 10(exp -16) T(sup -3.75) e(sup -300/T), T < 300 K) and compare methyl radical abundances in five atmospheres: Jupiter, Saturn, Uranus, Neptune, and Titan. The sensitivity of our models to the rate coefficients for the reactions H + CH3 + M yields CH4 + M, H + C2H3 yields C2H2 + H2, (sup 1)CH2 + H2 yields CH3 + H, and H + C2H5 yields 2CH3, the branching ratios of CH4 photolysis, vertical mixing in the five atmospheres, and Lyman alpha photon enhancement at the orbit of Neptune have all been tested. The results of our model CH3 abundances for both Saturn (5.1 x 10(exp 13) / sq cm) and Neptune (2.2 x 10(exp 13) / sq cm) show good agreement with ISO Short Wavelength Spectrometer measurements. Using the same chemical reaction set, our calculations also successfully generate vertical profiles of stable hydrocarbons consistent with Voyager and ground-based measurements in these outer solar system atmospheres. Predictions of CH3 column concentrations (for p <= 0.2 mbar) in the atmospheres of Jupiter (3.3 x 10(exp 13) /sq cm), Uranus (2.5 x 10(exp 12) / sq cm), and Titan (1.9 x 10(exp 15) / sq cm) may be checked by future observations.

Lee, Anthony Y. T.↗

Molecular beam mass spectrometer development

An analytical model, based on the kinetics theory of a drifting Maxwellian gas is used to determine the nonequilibrium molecular density distribution within a hemispherical shell open aft with its axis parallel to its velocity. The concept of a molecular shield in terrestrial orbit above 200 km is also analyzed using the kinetic theory of a drifting Maxwellian gas. Data are presented for the components of the gas density within the shield due to the free stream atmosphere, outgassing from the shield and enclosed experiments, and atmospheric gas scattered off a shield orbiter system. A description is given of a FORTRAN program for computating the three dimensional transition flow regime past the space shuttle orbiter that employs the Monte Carlo simulation method to model real flow by some thousands of simulated molecules.

Brock, F. J.↗

Molecular Hydrogen in Diffuse Interstellar Clouds of Arbitrary Three-Dimensional Geometry

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen in diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photodissociation rate of H2 in such clouds was computed with a 26 ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H2 line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 micron continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 micron continuum intensity to H I column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be the largest. Using a standard chi square analysis to assess the goodness of fit of our models, we find (at the 60 a level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be much smaller than its extent in the plane of the sky. If D is the distance to G236 + 39, and given standard assumptions about the rate of grain-catalyzed H2 formation, we find that the cloud has an extent along the line of sight that is 0.9 +/- 0.1 times its mean extent projected onto the plane of the sky and a gas density of 53 +/- 8 (100 pc/D) H nuclei/cc and is illuminated by a radiation field of 1.1 +/- 0.2 (100 pc/D) times the mean interstellar radiation field. The derived 100 micron emissivity per nucleon is 1.13 +/- 0.06 x 10(exp -20) MJy/sr sq cm.

Spaans, Marco↗

Observation-constrained kinetic modeling of isoprene SOA formation in the atmosphere

Isoprene has the largest global non-methane hydrocarbon emission, and the oxidation of isoprene plays a crucial role in the formation of secondary organic aerosol (SOA). Two primary processes are known to contribute to SOA formation from isoprene oxidation: (1) the reactive uptake of isoprene-derived epoxides on acidic or aqueous particle surfaces and (2) the absorptive gas–particle partitioning of low-volatility oxidation products. In this study, we developed a new multiphase condensed isoprene oxidation mechanism that includes these processes with key molecular intermediates and products. The new mechanism was applied to simulate isoprene gas-phase oxidation products and SOA formation from previously published chamber experiments under a variety of conditions and atmospheric observations during the Southern Oxidant and Aerosol Studies (SOAS) field campaign. Our results show that SOA formation from most of the chamber experiments is reasonably reproduced using our mechanism, except when the concentration ratios of initial nitric oxide to isoprene exceed ~ 2, the formed SOA is significantly underpredicted. The SOAS simulations also reasonably agree with the measurements regarding the diurnal pattern and concentrations of different product categories, while the total isoprene SOA remains underestimated. The molecular compositions of the modeled SOA indicate that multifunctional low-volatility products contribute to isoprene SOA more significantly than previously thought, with a median mass contribution of ~ 57 % to the total modeled isoprene SOA. However, this contribution is intricately intertwined with IEPOX-derived SOA (IEPOX: isoprene-derived epoxydiols), posing challenges for their differentiation using bulk aerosol composition analysis (e.g., the aerosol mass spectrometer with positive matrix factorization). Furthermore, the SOA from these pathways may vary greatly, mainly dependent on the volatility estimation and treatment of particle-phase processes (i.e., photolysis and hydrolysis). Our findings emphasize that the various pathways to produce these low-volatility species should be considered in models to more accurately predict isoprene SOA formation. The new condensed isoprene chemical mechanism can be further incorporated into regional-scale air quality models, such as the Community Multiscale Air Quality Modelling System (CMAQ), to assess isoprene SOA formation on a larger scale.

54 ENVIRONMENTAL SCIENCES↗

Ion Pairing, Clustering and Transport in a LiFSI-TMP Electrolyte as Functions of Salt Concentration using Molecular Dynamics Simulations

Battery capacity is highly related to ion-pairing mechanisms in electrolytes, since a cluster formation can lead to dead Li formation, reducing the number of charge carriers and leading to capacity fading. We use molecular dynamics simulations to model an electrolyte comprising trimethyl phosphate (TMP) solvent and a lithium bis(fluorosulfonyl)imide (LiFSI) salt, exploring effects of salt concentration on solvation and ion-transport. We simulate the LiFSI-TMP electrolyte for salt concentrations of 0.7, 1.43 and 3.82 molar. A statistical analysis was performed to study ion-pairing, clustering, diffusivity, conductivity, and coordination of Li-ions, providing insights into relations between molecular structures and transport properties. Molecular structure of ionic components changes as concentration increases, from a predominant solvent separated ion pair (SSIP) and contact ion pair (CIP) to aggregate salt (AGG) and ionic cluster formation. Given the formation of the ionic cluster, the diffusion mechanism followed by Li-ions changes from a hopping/exchange to a vehicular mechanism as concentration increases; this is reflected in a decrease of ionic conductivities. Ionicity was also calculated to reveal how the ionic motion changes from an uncorrelated to a correlated one as the salt concentration increases. Furthermore, we also compared our results with experimental calculations performed for similar electrolyte systems

25 ENERGY STORAGE↗

Modeling the aqueous interface of amorphous TiO 2 using deep potential molecular dynamics

Amorphous titanium dioxide (a-TiO 2 ) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO 2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO 2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO 2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO 2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO 2 , leading to an ~10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti 2 –O b H) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti–O w H) due to fast Ti–O w H 2 → Ti–O w H proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO 2 in electrochemical environments. Moreover, the procedure of generating the a-TiO 2 -interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of molecular anisotropy on the intensity and degree of polarization of light scattered from model atmospheres

Computations of the properties of sunlight scattered from models of the earth-atmosphere system are presented to show the effect of molecular anisotropy on the intensity, flux, and degree of polarization of the scattered light. The values of these parameters change significantly when the anisotropy factor is neglected in the molecular optical thickness and scattering phase matrix. However, if the Rayleigh scattering optical thickness is kept constant and the molecular anisotropy factor is included only in the Rayleigh phase matrix, the flux does not change, the intensity changes by a small amount, but the changes in the degree of polarization are still significant.

Bahethi, O. P.↗

Theoretical Modeling and Computer Simulations for the Origins and Evolution of Reproducing Molecular Systems and Complex Systems with Many Interactive Parts

Our research effort has produced nine publications in peer-reviewed journals listed at the end of this report. The work reported here are in the following areas: (1) genetic network modeling; (2) autocatalytic model of pre-biotic evolution; (3) theoretical and computational studies of strongly correlated electron systems; (4) reducing thermal oscillations in atomic force microscope; (5) transcription termination mechanism in prokaryotic cells; and (6) the low glutamine usage in thennophiles obtained by studying completely sequenced genomes. We discuss the main accomplishments of these publications.

Liang, Shoudan↗

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Development of Catalysts for Upgrading of Pyrolysis Vapor for Refinery Feedstocks and Intermediates (CRADA Final Report)

Catalytic fast pyrolysis (CFP) is a versatile technology platform to convert biomass into fungible hydrocarbon transportation fuels and chemical co-products. Key technical barriers to reaching this goal include increasing the product yields and achieving the desired fuel properties for gasoline, diesel, and jet range fuels or blendstocks that would be suitable for introduction into existing refinery unit operations. Overcoming these barriers will require durable catalysts that are effective at upgrading and stabilizing biomass pyrolysis vapors. Towards these goals, this CRADA leveraged NREL experience as a leader in biomass pyrolysis research and Johnson Matthey's (JM) experience as a leader in the production of advanced catalytic materials. The scope spanned CFP catalyst development, characterization, multi-scale reaction testing, and computational modeling. CRADA benefits to DOE, Participant, and U.S. Taxpayer: Assists laboratory in achieving programmatic scope, Uses the laboratory’s core competencies. The purpose of this CRADA was to develop and deploy catalysts for biomass CFP to help achieve cost-competitive biofuels and bio-based products. This was accomplished through a close collaboration between biomass conversion researchers at NREL and catalyst development researchers at JM. Summary of Research Results: Focus Area 1. Foundational research on catalytic conversion and deactivation: Key interactions between pyrolysis vapors and heterogeneous catalysts were probed through catalyst characterization, model compound reaction testing, and atomistic-scale computational modeling. Catalyst development focused on multifunctional materials, which include zeolites, oxides, carbides, and nitrides. Computational modeling identified reaction mechanisms and elucidated surface chemistry to test hypotheses regarding mechanisms of deoxygenation, coupling, cracking, dehydration, coke formation, hydrogen transfer, and aromatic ring reactions. This information was used to design multifunctional catalysts to increase product yields, control product selectivity, and reduce deactivation during CFP and downstream processing steps. The results served to increase fundamental understanding of key catalyst attributes and durability features for the upgrading of biomass pyrolysis vapors. Model compound experiments confirmed the importance of metal-acid bifunctionality for the deoxygenation of lignin-derived phenolic species under hydrodeoxygenation conditions. This insight led to the development of catalysts such as Pt/TiO2 and Mo2C, which were confirmed as high-performing materials during subsequent bench-scale experiments using biomass-derived pyrolysis vapors. This focus area also led to the identification of important catalyst deactivation mechanisms associated with the deposition of inorganic contaminants such as potassium. The molecular-level insight from model compound experiments and computational modeling, shown in Figure 1, informed the development of regeneration procedures that have been shown to be effective for restoration of > 90% of initial catalyst activity. This understanding has subsequently been translated to other catalyst systems, including zeolite materials that can be operated without requirements for co-fed hydrogen.

09 BIOMASS FUELS↗

Targeting tissues via dynamic human systems modeling in generative design

Drug discovery is a complex, costly process with high failure rates. A successful drug should bind to a target, be deliverable to an intended site of activity, and promote a desired pharmacological effect without causing toxicity. Typically, these factors are evaluated in series over the course of a pipeline where the number of candidates is sequentially whittled down from a very large initial pool. One promise of AI-driven discovery is the opportunity to evaluate multiple facets of drug performance in parallel. However, despite ML-driven advancements, current models for pharmacological property prediction are exclusively trained to predict molecular properties, ignoring important, dynamic biodistribution and bioactivity effects. Here, we present our progress towards incorporating quantitative systems physiology models into an AI-enabled drug design and molecular generation pipeline. Within a genetic algorithm, we include human-relevant physiologically based pharmacokinetic (PBPK) models. These PBPK models leverage properties that are predicted by a fine-tuned molecular language model. Together, these models will aid in capturing the mapping between molecules and therapeutic outcomes that is necessary to accelerate the drug discovery process.

Fox, Zach↗

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A comprehensive model of ion diffusion and charge exchange in the cold Io torus

A comprehensive analytic model of radial diffusion in the cold Io torus is developed. The model involves a generalized molecular cloud theory of SO2 and its dissociation fragments SO, O2, S, and O, which are formed at a relatively large rate by solar UV photodissociation of SO2. The key component of the new theory is SO, which can react with S(+) through a near-resonant charge exchange process that is exothermic. This provides a mechanism for the rapid depletion of singly ionized sulfur in the cold torus and can account for the large decrease in the total flux tube content inward of Io's orbit. The model is used to demonstrate quantitatively the effects of radial diffusion in a charge exchange environment that acts as a combined source and sink for ions in various charge states. A detailed quantitative explanation for the O(2+) component of the cold torus is given, and insight is derived into the workings of the so-called plasma 'ribbon'.

Barbosa, D. D.↗