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,207 records · Page 67

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.↗

Simulation of adiabatic quantum computing for molecular ground states

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.

Kremenetski, Vladimir↗

Multifunctional Chiral Chemically‐Powered Micropropellers for Cargo Transport and Manipulation

Practical applications of synthetic self-propelled nano and microparticles for microrobotics, targeted drug delivery, and manipulation at the nanoscale are rapidly expanding. However, fabrication limitations often hinder progress, resulting in relatively simple shapes and limited functionality. Here, taking advantage of 3D nanoscale printing, chiral micropropellers powered by the hydrogen peroxide reduction reaction are fabricated. Due to their chirality, the propellers exhibit multifunctional behavior controlled by an applied magnetic field: spinning in place (loitering), directed migration in the prescribed direction, capture, and transport of polymer cargo particles. Design parameters of the propellers are optimized by computation modeling based on mesoscale molecular dynamics. It is predicted by computer simulations, and confirmed experimentally, that clockwise rotating propellers attract each other and counterclockwise repel. These results shed light on how chirality and shape optimization enhance the functionality of synthetic autonomous micromachines.

36 MATERIALS SCIENCE↗

Mechanical behavior in living cells consistent with the tensegrity model

Alternative models of cell mechanics depict the living cell as a simple mechanical continuum, porous filament gel, tensed cortical membrane, or tensegrity network that maintains a stabilizing prestress through incorporation of discrete structural elements that bear compression. Real-time microscopic analysis of cells containing GFP-labeled microtubules and associated mitochondria revealed that living cells behave like discrete structures composed of an interconnected network of actin microfilaments and microtubules when mechanical stresses are applied to cell surface integrin receptors. Quantitation of cell tractional forces and cellular prestress by using traction force microscopy confirmed that microtubules bear compression and are responsible for a significant portion of the cytoskeletal prestress that determines cell shape stability under conditions in which myosin light chain phosphorylation and intracellular calcium remained unchanged. Quantitative measurements of both static and dynamic mechanical behaviors in cells also were consistent with specific a priori predictions of the tensegrity model. These findings suggest that tensegrity represents a unified model of cell mechanics that may help to explain how mechanical behaviors emerge through collective interactions among different cytoskeletal filaments and extracellular adhesions in living cells.

NASA Discipline Cell Biology↗

Warm dense gas in the reflection nebula NGC 2023

A photodissociation region (PDR) with a known UV excitation source was studied by examining emission from atomic and molecular material. In agreement with models of externally excited PDRs, the photodissociation region southwest of HD 37903 contains both atomic and excited molecular material and arises close to the exciting star than does the bulk of the molecular cloud. The C(+) emission arises closest to the star but is 1.5 x 10 to the 17th cm or less from the peaks in the excited CO and fluorescent H2 emission regions. The (C-18)O J = 2 - 1 line, which traces molecular column density, has its emission peak about 4 x 10 to the 17th cm beyond the C(+) peak. Even in a source with an incident UV flux as low as about 1000 times the mean interstellar radiation field, there is a significant amount of CO 7 - 6 emission arising from regions with T(gas) greater than 85 K which is substantially greater than T(dust).

Jaffe, D. T.↗

Data Needs for Stellar Atmosphere and Spectrum Modeling

The main data need for stellar atmosphere and spectrum modeling remains atomic and molecular transition data, particularly energy levels and transition cross-sections. We emphasize that data is needed for bound-free (b - f) as well as bound-bound (b - b), and collisional as well as radiative transitions. Data is now needed for polyatomic molecules as well as atoms, ions, and diatomic molecules. In addition, data for the formation of, and extinction due to, liquid and solid phase dust grains is needed. A prioritization of species and data types is presented, and gives emphasis to Fe group elements, and elements important for the investigation of nucleosynthesis and Galactic chemical evolution, such as the -elements and n-capture elements. Special data needs for topical problems in the modeling of cool stars and brown dwarfs are described.

Short, C. I.↗

Orientation-dependent etching of silicon by fluorine molecules: A quantum chemistry computational study

Anisotropic etching is a widely used process in semiconductor manufacturing, in particular, for micro- and nanoscale texturing of silicon surfaces for black silicon production. The typical process of plasma-assisted etching uses energetic ions to remove materials in the vertical direction, creating anisotropic etch profiles. Plasmaless anisotropic etching, considered here, is a less common process that does not use ions and plasma. The anisotropy is caused by the unequal etching rates of different crystal planes; the etching process, thus, proceeds in a preferred direction. In this paper, we have performed quantum chemistry modeling of gas-surface reactions involved in the etching of silicon surfaces by molecular fluorine. The results confirm that orientation-dependent etch rates are the reason for anisotropy. The modeling of F 2 dissociative chemisorption on F-terminated silicon surfaces shows that Si–Si bond breaking is slow for the Si(111) surface, while it is fast for Si(100) and Si(110) surfaces. Both Si(100) and Si(110) surfaces incorporate a larger number of fluorine atoms resulting in Si–Si bonds having a larger amount of positive charge, which lowers the reaction barrier of F 2 dissociative chemisorption, yielding a higher etch rate for Si(100) and Si(110) surfaces compared to Si(111) surfaces. Molecular dynamics modeling of the same reactions has shown that the chosen reactive bond order potential does not accurately reproduce the lower reaction barriers for F 2 dissociative chemisorption on Si(100) and Si(100) surfaces. Furthermore, reparameterization is necessary to model the anisotropic etching process that occurs at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Evolution of Mass and Physical Properties of Atmospheric Organic Aerosol under Solar Irradiance

Organic aerosol (OA) particles constitute a substantial fraction of sub-micron particulate mass in the atmosphere and play a critical role in climate system. OA undergoes dynamic aging processes in the atmosphere, with photolytic aging induced by ultraviolet solar irradiance being an important yet poorly characterized mechanism. Knowledge gaps persist regarding the role of volatility transformations during photolytic aging on the OA mass decay kinetics and the evolution of climate-relevant properties, such as hygroscopicity, hindering the model evaluation of OA spatiotemporal distributions and atmospheric budgets. In this study, we conduct isothermal photolytic aging experiments on both laboratory-generated secondary organic aerosols and ambient-collected particles from urban Atlanta, utilizing a high-sensitivity Quartz Crystal Microbalance. Our results reveal that photolytic aging reduces 40–66% of the low-volatility OA mass with lifetimes ranging from 8 to 200 hours under solar irradiance, and 44–92% of the photolytic mass loss is through slow evaporation of semi- or intermediate-volatile products, kinetically limited by their volatility. We observe up to ±50% changes in OA hygroscopicity with the transformation of fresh OA to photo-recalcitrant low-volatility products, associated with changes in oxygen-to-carbon ratio and molecular weight. A kinetic model incorporating photolytic volatility transformation provides a cohesive explanation for the observed photolysis-induced changes in mass, volatility, and hygroscopicity. Our results can help constrain model representation of the dynamic evolutions of mass and climate-relevant properties during photolytic aging processes of the ambient OA, improving our understanding of OA atmospheric behavior and climate impact.

Bai, Bin↗

Nonadiabatic dissociation of molecular Bose-Einstein condensates: Competition between chemical reactions

Here, we provide a framework to solve generic models describing the dissociation of multiple molecular Bose Einstein condensates in a nonadiabatic regime. The competition between individual chemical reactions can lead to non-trivial dependence on critical components such as path interference and symmetries, thus, affecting the final distribution of atomic population. We find an analytical solution for an illustrative example model involving four atomic modes. When the system parameters satisfy CPT symmetry, where C is charge conjugation, P is parity, and T is time-reversal symmetry, our solution predicts a population imbalance between atomic modes that is exponentially sensitive to system parameters. However, a weakly broken symmetry alters the population in each atomic mode and can reverse the population imbalance. Our solution also demonstrates a strong quantum correlation between atomic modes that leads to the spontaneous production of atoms in a multi-mode squeezed state. Moreover, in our framework, a time-dependent non-Hermitian quantum mechanics naturally manifests which can alternatively be realized experimentally in photonic systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA Space Radiation Program Integrative Risk Model Toolkit

NASA Space Radiation Program Element scientists have been actively involved in development of an integrative risk models toolkit that includes models for acute radiation risk and organ dose projection (ARRBOD), NASA space radiation cancer risk projection (NSCR), hemocyte dose estimation (HemoDose), GCR event-based risk model code (GERMcode), and relativistic ion tracks (RITRACKS), NASA radiation track image (NASARTI), and the On-Line Tool for the Assessment of Radiation in Space (OLTARIS). This session will introduce the components of the risk toolkit with opportunity for hands on demonstrations. The brief descriptions of each tools are: ARRBOD for Organ dose projection and acute radiation risk calculation from exposure to solar particle event; NSCR for Projection of cancer risk from exposure to space radiation; HemoDose for retrospective dose estimation by using multi-type blood cell counts; GERMcode for basic physical and biophysical properties for an ion beam, and biophysical and radiobiological properties for a beam transport to the target in the NASA Space Radiation Laboratory beam line; RITRACKS for simulation of heavy ion and delta-ray track structure, radiation chemistry, DNA structure and DNA damage at the molecular scale; NASARTI for modeling of the effects of space radiation on human cells and tissue by incorporating a physical model of tracks, cell nucleus, and DNA damage foci with image segmentation for the automated count; and OLTARIS, an integrated tool set utilizing HZETRN (High Charge and Energy Transport) intended to help scientists and engineers study the effects of space radiation on shielding materials, electronics, and biological systems.

Kim, Myung-Hee Y.↗

Development of Thermal Scattering Law and Cross Sections of Hydrogen in Paraffin Wax

Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Implications of SWAS Observations for Interstellar Chemistry and Star Formation

A long standing prediction of steady state gas-phase chemical theory is that H2O and O2 are important reservoirs of elemental oxygen and major coolants of the interstellar medium. Analysis of SWAS observations has set sensitive upper limits on the abundance Of O2 and has provided H2O abundances toward a variety of star forming regions. Based on these results, we show that gaseous H2O and O2 are not dominant carriers of elemental oxygen in molecular clouds. Instead the available oxygen is presumably frozen on dust grains in the form of molecular ices, with a significant portion potentially remaining in atomic form, along with CO, in the gas phase. H2O and O2 are also not significant coolants for quiescent molecular gas. In the case of H2O, a number of known chemical processes can locally elevate its abundance in regions with enhanced temperatures, such as warm regions surrounding young stars or in hot shocked gas. Thus, water can be a locally important coolant. The new information provided by SWAS, when combined with recent results from the Infrared Space Observatory, also provide several hard observational constraints for theoretical models of the chemistry in molecular clouds and we discuss various models that satisfy these conditions.

Bergin, Edwin A.↗

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗

Negative capacitance regime in ferroelectrics demystified from nonequilibrium molecular dynamics

Negative capacitance regime in ferroelectrics is associated with the negative slope in the dependence of the polarization on electric field and is prohibited by the laws of equilibrium thermodynamics. Very recently, it developed into a controversial but very productive area with a promise for next generation low-power transistors. At least, in part, the controversies arise from the lack of direct experimental insight into the mechanism for such phenomena. Although the experimental insight is, indeed, challenging at present, we probe the effect in ferroelectric PbTiO 3 bulk and a thin film with realistic first-principles-based nonequilibrium molecular dynamics simulations which model most closely experimental conditions. Here, the following insights emerge from the simulations: (i) the ferroelectric never enters the thermodynamically forbidden region even under extreme rates of electric-field application; (ii) the negative slope naturally arises when the compensating charge is either scarce or slowed down by the external circuitry; (iii) the negative slope is well controlled by the interplay between the applied step pulse field and the amount of the screening charge. The latter can be used to design the negative capacitance regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗

Design of Efficient Molecular Electrocatalysts for Water and Carbon Dioxide Reduction Using Predictive Models of Thermodynamic Properties

Objectives: 1) To design and synthesize water soluble transition metal complexes that function as electrocatalysts for the reduction of water to hydrogen or carbon dioxide to formate 2) To experimentally measure thermochemical properties, such as hydride donor ability, and use these quantities to optimize catalyst performance in solutions of varying pH 3) To incorporate results from detailed mechanistic and kinetic studies to direct improvements in catalyst design 4) To use thermochemical data from objective 2 to generate predictive models for the properties of new complexes to apply to the development of catalysts for other reductive reactions, with an emphasis on carbon dioxide reduction to methanol Description: Wide-spread implementation of renewable but intermittent energy sources, such as solar, requires the development of efficient methods for energy storage. The high energy density of chemical bonds makes chemical fuels an ideal solution for energy storage. Hydrogen and reduced carbon compounds have been proposed as ideal candidates for chemical fuels. However, the generation of chemical fuels from electricity requires competent electrocatalysts. The proposed research focuses on developing electrocatalysts for the reduction of water to hydrogen, and carbon dioxide to formate. Both products can be used directly as an energy carrier in fuel cells, or as a reductant for more saturated chemical fuels. Hydrogen is the most common fuel used in current commercial fuel cells. Additionally, there is interest in formic acid as a liquid carrier for hydrogen because of its increased storage density and ease of dehydrogenation. Formate is also an intermediate in the sequential reduction of CO 2 to methanol, another potential chemical fuel with high energy storage density. Despite the utility of formate as a chemical fuel or precursor, there are very few examples of electrocatalysts for its production from CO 2 , and even fewer demonstrate high product selectivity. In heterogeneous catalysis, the Sabatier principle is used to generate volcano curves that describe optimal thermochemical properties for key intermediates that result in peak catalytic activity. In most cases, such as hydrogen production and oxidation, the most favorable metals (Pt, Re, Rh, and Ir) are rare and expensive. Molecular inorganic complexes provide an opportunity to use electronic and steric ligand effects to tune the critical thermodynamic parameters of abundant metals to values comparable to key surface intermediates on precious metals. This principle will be applied to the design of aqueous homogeneous catalysts for the reduction of H 2 O and CO 2 optimized to function at specific pH ranges. The critical intermediate in both of these reactions is a metal hydride. The strength of this bond, or hydricity (ΔGH - ) dictates the overall thermodynamics of the reduction of H + to H 2 and the sequential reduction of C1 substrates, such as CO 2 , CO, and H 2 CO (the latter to CH 3 OH). ΔGH - will be systematically measured for a series of first row metal complexes to form predictive models for metal and ligand electronic effects. Benefits and Outcomes: This approach addresses thermodynamic requirements essential for energy-efficient catalyst design and transition state barrier considerations for faster catalysis. The advancements made from this research will include the discovery of efficient, stable, and fast catalysts for hydrogen and formate generation from electricity. The ability to store energy efficiently is critical to the widespread use of renewable energy schemes. This provides a viable solution to produce high energy density fuels for stationary and mobile energy applications. The fundamental research conducted is also directly applicable to the design of catalysts for the generation of further reduced products such as methanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Density Amorphous Ice, the Frost on Interstellar Grains

Most water ice in the universe is in a form which does not occur naturally on Earth and of which only minimal amounts have been made in the laboratory. We have encountered this 'high-density amorphous ice' in electron diffraction experiments of low-temperature (T less than 30 K) vapor-deposited water and have subsequently modeled its structure using molecular dynamics simulations. The characteristic feature of high-density amorphous ice is the presence of 'interstitial' oxygen pair distances between 3 and 4 A. However, we find that the structure is best described as a collapsed lattice of the more familiar low-density amorphous form. These distortions are frozen in at temperatures below 38 K because, we propose, it requires the breaking of one hydrogen bond, on average, per molecule to relieve the strain and to restructure the lattice to that of low-density amorphous ice. Several features of astrophysical ice analogs studied in laboratory experiments are readily explained by the structural transition from high-density amorphous ice into low-density amorphous ice. Changes in the shape of the 3.07 gm water band, trapping efficiency of CO, CO loss, changes in the CO band structure, and the recombination of radicals induced by low-temperature UV photolysis all covary with structural changes that occur in the ice during this amorphous to amorphous transition. While the 3.07 micrometers ice band in various astronomical environments can be modeled with spectra of simple mixtures of amorphous and crystalline forms, the contribution of the high-density amorphous form nearly always dominates.

Jenniskens, P.↗

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗