SEARCH · Search NASA
Results for “materials and molecular modeling”
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.
Molecular dynamics simulation of vacancy and void effects on strain-induced martensitic transformations in Fe-50 at.% Ni model concentrated solid solution alloy
Strain-induced martensitic transformations can improve the strength and ductility of some face centered cubic (FCC) metals and alloys. Irradiation-induced defects such as vacancies, dislocation loops, and voids activate martensitic transformations over a wider range of mechanical loading conditions than in pristine material. Furthermore, the mechanisms underlying irradiation-enabled martensite transformations remain unclear. In this work, we use molecular dynamics simulations to study the effect of vacancies and voids on strain-induced martensitic transformations in a model concentrated solid solution alloy Fe-50 at.% Ni. It has been found that single vacancies have no resolvable effect on the transformation because they reduce the stacking fault energy by a relatively insignificant margin and do so only if the vacancy is located on the stacking fault plane. Voids, however, activate the martensite transformation through shear strain accumulation around the void due to dislocation pileup. The larger the void, the more pronounced this effect.
U.S. Department of Energy: National Virtual Biotechnology Laboratory (Technical Report)
With funding from the CARES Act, the U.S Department of Energy (DOE) established the National Virtual Biotechnology Laboratory (NVBL) in March 2020 to address key challenges associated with the COVID-19 crisis. The NVBL brought together the broad scientific and technical expertise and resources of DOE’s 17 national laboratories to help tackle medical supply shortages, discover potential drugs to fight the virus, develop and validate COVID-19 testing methods, model disease spread and impact across the nation, and understand virus transport in buildings and the environment. National laboratory resources leveraged for this effort include a suite of world-leading user facilities broadly available to the research community, such as light and neutron sources, nanoscale science research centers, sequencing and biocharacterization facilities, and high-performance computing facilities. As part of the NVBL framework, DOE rapidly assembled five project teams to (1) identify new targets for medical therapeutics; (2) develop innovations in testing capabilities; (3) provide epidemiological and logistical support; (4) understand viral fate and transport in the environment; and (5) address supply chain bottlenecks by harnessing extensive additive manufacturing capabilities. Each research team was charged with defining high-impact projects that could be completed in a 6-month sprint while coordinating their developments with academia, other government agencies, and the private sector. Within months, NVBL teams used DOE’s high-performance computers and light and neutron sources to identify promising candidates for antibodies and antivirals that universities and drug companies are now evaluating. NVBL researchers also developed new diagnostic targets and sample collection approaches, and supported efforts by the U.S. Food and Drug Administration, Centers for Disease Control and Prevention, and U.S. Department of Defense to establish national guidelines used in administering millions of tests. Researchers used artificial intelligence and high-performance computing to produce near-real-time data analysis to forecast disease transmission, stress on public health infrastructure, and economic impact, which supported decision-makers at the local, state, and national levels. To minimize virus uptake and protect human health, NVBL teams studied how to control indoor virus movement. Researchers also produced innovations in materials and advanced manufacturing that mitigated shortages in test kits and personal protective equipment, creating nearly 1,000 new jobs. Through its NVBL framework, DOE has contributed significantly to the nation’s COVID response, demonstrating in only a few months the critical impact of its national laboratories. NVBL’s accomplishments demonstrate not only the powerful resource represented by DOE’s national laboratories working together to meet national needs, but also the effectiveness of the integrated NVBL framework for rapidly responding to emergencies with research and development solutions. Going forward, the NVBL is poised to apply the unique capabilities and expertise of the national laboratory complex to future national and international emergencies, both natural and engineered. Through this framework, DOE will continue to be an integral component of agency-wide efforts to prepare for and respond to biorisks and other crises. This technical report describes the goals, progress, and results of NVBL’s five project teams—Molecular Design for COVID-19 Therapeutics, COVID-19 Testing, Epidemiological Modeling, Viral Fate and Transport, and Materials and Manufacturing of Critical Supplies—and lists each team’s publications and research output.
N2H(+) in warm and cold clouds
A topic of importance in both cometary and instellar astronomy is the N2/NH3 abundance ratio in star forming regions. A key test for solar system formation models is the N2/NH3 abundance ratio in comets, which are thought to be made up of material from the protostar nebula, and in molecular clouds, from which stars are thought to be approx. 0.1. The N2/NH3 abundance ratio has never been measured in the interstellar medium, but is estimated from theoretical models to be in the range of 100 to 200. If the theoretical models are truly representative of dense molecular clouds then the difference in the abundance ratios of N2/NH3 between comets and the ISM would imply that significant chemical processing must have occurred in the protostar nebula that somehow altered the amount of N2 relative to NH3. A systematic survey was made of star forming regions to establish N2 abundances and N2/NH3 ratios for comparison with cometary data.
Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel
Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.
Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces
As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.
Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon
The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.
MISPR : an open-source package for high-throughput multiscale molecular simulations
Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.
Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework
Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.
Envelope structure of deeply embedded young stellar objects in the Serpens Molecular Cloud
Aperture-synthesis and single-dish (sub-) millimeter molecular-line and continuum observations reveal in great detail the envelope structure of deeply embedded young stellar objects (SMM 1 = FIRS 1, SMM 2, SMM 3, SMM 4) in the densely star-forming Serpens Molecular Cloud. SMM 1, 3, and 4 show partially resolved (>2" = 800 AU) continuum emission in the beam of the Owens Valley Millimeter Array at lambda = 3.4-1.4 mm. The continuum visibilities accurately constrain the density structure in the envelopes, which can be described by a radial power law with slope -2.0 +/- 0.5 on scales of 300 to 8000 AU. Inferred envelope masses within a radius of 8000 AU are 8.7, 3.0, and 5.3 Msolar for SMM 1, 3, and 4, respectively. A point source with 20%-30% of the total flux at 1.1 mm is required to fit the observations on long baselines, corresponding to warm envelope material within approximately 100 AU or a circumstellar disk. No continuum emission is detected interferometrically toward SMM 2, corresponding to an upper limit of 0.2 Msolar assuming Td = 24 K. The lack of any compact dust emission suggests that the SMM 2 core does not contain a central protostar. Aperture-synthesis observations of the 13CO, C18O, HCO+, H13CO+, HCN, H13CN, N2H+ 1-0, SiO 2-1, and SO 2(2)-1(1) transitions reveal compact emission toward SMM 1, 3, and 4. SMM 2 shows only a number of clumps scattered throughout the primary field of view, supporting the conclusion that this core does not contain a central star. The compact molecular emission around SMM 1, 3, and 4 traces 5"-10" (2000-4000 AU) diameter cores that correspond to the densest regions of the envelopes, as well as material directly associated with the molecular outflow. Especially prominent are the optically thick HCN and HCO+ lines that show up brightly along the walls of the outflow cavities. SO and SiO trace shocked material, where their abundances may be enhanced by 1-2 orders of magnitude over dark-cloud values. A total of 31 molecular transitions have been observed with the James Clerk Maxwell and Caltech Submillimeter telescopes in the 230, 345, 490, and 690 GHz atmospheric windows toward all four sources, containing, among others, lines of CO, HCO+, HCN, H2CO, SiO, SO, and their isotopomers. These lines show 20-30 km s-1 wide line wings, deep and narrow (1-2 km s-1) self-absorption, and 2-3 km s-1 FWHM line cores. The presence of highly excited lines like 12CO 4-3 and 6-5, 13CO 6-5, and several H2CO transitions indicates the presence of material with temperatures > or approximately 100 K. Monte Carlo calculations of the molecular excitation and line transfer show that the envelope model derived from the dust emission can successfully reproduce the observed line intensities. The depletion of CO in the cold gas is modest compared to values inferred in objects like NGC 1333 IRAS 4, suggesting that the phase of large depletions through the entire envelope is short lived and may be influenced by the local star formation density. Emission in high-excitation lines of CO and H2CO requires the presence of a small amount of approximately 100 K material, comprising less than 1% of the total envelope mass and probably associated with the outflow or the innermost region of the envelope. The derived molecular abundances in the warm (Tkin > 20 K) envelope are similar to those found toward other class 0 YSOs like IRAS 16293-2422, though some species appear enhanced toward SMM 1. Taken together, the presented observations and analysis provide the first comprehensive view of the physical and chemical structure of the envelopes of deeply embedded young stellar objects in a clustered environment on scales between 1000 and 10,000 AU.
Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal
Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.
Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations
Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.
Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials
This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.
A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes
Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.
Comparison of water nanodroplet properties on different graphite-based substrates
The molecular structure and dynamics of water differ considerably at various interfaces. We compare the interfacial water structure–property relationship on three different carbon substrates, namely, amorphous carbon, compressed expanded natural graphite, and pure graphite by utilizing atomistic molecular dynamics simulations. The effect of different substrates on the structural and dynamical properties of water can readily be observed. Here, the density distributions parallel and normal to the substrates show oblate droplet structures. The normal to the substrate water distribution shows a strong hydration layer at the interface that does not vary with substrates. However, the disparity in the structure and dynamics on three different substrates shows that the surface morphologies of the substrates are critical for determining nanoscale water properties. Furthermore, it is observed that the formation of an interfacial water layer or the hydration layer is a direct consequence of both water “confinement” at the nanoscale and “attraction” between water molecules and the carbon substrates.
High dopant activation in arsenic doped single-crystal CdTe thin films: Insights from MBE growth and rapid thermal processing
Single-crystal model systems are valuable tools to investigate fundamental material properties. In this work, we use molecular beam epitaxy to deposit in situ arsenic (As) doped single-crystal CdTe films on large area Si substrates to better understand As doping for photovoltaic applications. We found that As incorporation is highly temperature dependent: a substrate temperature difference of 50 °C can lead to several orders of magnitude difference in As concentration. Cd overpressure during in situ doping may limit out-diffusion of As but decrease As incorporation, especially at lower growth temperatures. Carrier concentrations greater than 10 16 cm −3 can be achieved with or without Cd overpressure when annealed at temperatures above 500 °C. However, unlike the low (∼1% to 5%) dopant activation commonly observed in polycrystalline CdTe, our films achieve significantly higher activation ratios—exceeding 50%, and in some cases approaching 80%. These values are consistent with or exceed prior reports in single-crystal CdTe systems. In addition to as-deposited arsenic concentrations, we also consider arsenic distribution after different rapid thermal processing temperatures. We propose a detailed definition and description of how arsenic incorporation is considered and calculated. Due to carrier concentration saturation, As incorporation also needs to be controlled to average levels of 10 17 cm −3 to achieve high activation. These findings suggest that higher annealing temperature regimes may be beneficial to polycrystalline CdTe based PV devices.
Scattering theory of frequency-entangled biphoton states facilitated by cavity polaritons
The use of quantum light to probe exciton properties in semiconductor and molecular nanostructures typically occurs in the low-intensity regime. A substantial enhancement of exciton-photon coupling can be achieved with photonic cavities, where excitons hybridize with cavity modes to form polariton states. Here, to provide a theoretical framework for interpreting emerging experimental efforts in this direction, we develop a scattering theory describing the interaction of frequency-entangled photon pairs with cavity polariton and bipolariton states under various coupling regimes. Employing the Tavis-Cummings model in combination with our scattering approach, we present a quantitative analysis of how the interaction of the entangled photon pair with the polariton or bipolariton modifies its joint spectral amplitude (JSA). Specifically, we examine the effects of the cavity-mode steady-state population, exciton-cavity coupling strength, and different forms of the input photon JSA. Our results show that the entanglement entropy of the scattered photons is highly sensitive to the interplay between the input JSA and the spectral line shapes of the polariton resonances, emphasizing the cavity filtering effects. We suggest that biphoton-scattering quantum light spectroscopy best serves as a sensitive probe of polariton and bipolariton states in the photon-vacuum cavity state. Our approach is not only robust to various regimes of cavity-exciton coupling, but also amenable to extensions beyond the Tavis-Cummings model, enabling the representation of a broad class of molecular systems and solid state quantum materials.
Advanced Physics Models for Particle-Particle Interactions (UCD final report)
The LDRD "Advanced Physics Models for Particle-Particle Interactions” collaborative research with Professor Jennifer Sinclair Curtis and PhD student Chanh Nguyen at the University of California, Davis involved meso-scale, discrete element method (DEM) modeling of micron-sized particles impacting a substrate. Using DEM, physical insight into the details of particle response upon impact can be obtained. Such details provide the type of information necessary to enhance existing continuum or macro-scale simulation tools so that modeling for material strength and particle breakage models is based on particle physics versus empirical relationships. While molecular dynamics simulations of atoms can also provide physical insight into the particle response to impact, the scale of particles that are treated is on the order of angstroms, not microns. This study gives evidence to the ability of the DEM simulation technique to provide such fundamental, detailed insight for micron-sized particles as a function of particle properties and impact velocity.