Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulation”

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

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Designing polar textures with ultrafast neuromorphic features from atomistic simulations

This review summarizes recent works, all using a specific atomistic approach, that predict and explain the occurrence of key features for neuromorphic computing in three archetypical dipolar materials, when they are subject to THz excitations. The main ideas behind such atomistic approach are provided, and illustration of model relaxor ferroelectrics, antiferroelectrics, and normal ferroelectrics are given, highlighting the important potential of polar materials as candidates for neuromorphic computing. Some peculiar emphases are made in this Review, such as the connection between neuromorphic features and percolation theory, local minima in energy path, topological transitions and/or anharmonic oscillator model, depending on the material under investigation. By considering three different and main polar material families, this work provides a complete and innovative toolbox for designing polar-based neuromorphic systems.

THz pulses↗
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Erratum: “DFTB+, a software package for efficient approximate density functional theory based atomistic simulations” [J. Chem. Phys. 152, 124101 (2020)]

The implementation of the GPU support in DFTB+, as described in Sec. III C of the original publication,1 was developed based on a previous unpublished implementation by Jacek Jakowski. In order to acknowledge his work on this first implementation, the authors of the original publication wish to include J. Jakowski as co-author. The scientific content of the original publication is not affected.

42 ENGINEERING↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2023

Mechanism of monovalent and divalent ion mobility in Nafion membrane: An atomistic simulation study

Polymer cation exchange membranes (CEMs) are widely used in water treatment processes. The fundamental factors that control the separation and selectivity of cations with different valences in CEMs are not fully understood. In this study, we use atomistic molecular dynamics simulations to investigate the underlying molecular mechanisms that control the mobility of cations with different valences in Nafion membranes. Our results indicate structural differences in binding of monovalent and divalent cations, which gives rise to differences in ion diffusion in Nafion. Monovalent cations are shown to be “territorially” bound, i.e., bound while partially hydrated, to the fixed charge groups whereas divalent cations are shown to be “site” bound, i.e., bound while fully dehydrated, to the charge groups on the polymer. This difference in binding structure gives rise to differences in transport characteristics of cations in Nafion.

Chemistry↗
DOE OSTI2025

Unraveling Exciton Trap Dynamics and Nonradiative Loss Pathways in Quantum Dots via Atomistic Simulations

Surface defects in colloidal quantum dots are a major source of nonradiative losses, yet the microscopic mechanisms underlying exciton trapping and recombination remain elusive. Here, we develop a model Hamiltonian based on atomistic electronic calculations to investigate exciton dynamics in CdSe/CdS core/shell QDs containing a single hole trap introduced by an unpassivated sulfur atom. By systematically varying the defect depth and reorganization energy, we uncover how defect-induced excitonic states mediate energy relaxation pathways. Our simulations reveal that a single localized defect can induce a rich spectrum of excitonic states, leading to multiple dynamical regimes, from slow, energetically off-resonant trapping to fast, cascaded relaxation through in-gap defect states. Crucially, we quantify how defect-induced polaron shifts and exciton-phonon couplings govern the balance between efficient radiative emission and rapid nonradiative decay. These insights clarify the microscopic origin of defect-assisted loss channels and suggest pathways for tailoring QD optoelectronic properties via surface and defect engineering.

Defects↗
DOE OSTI2023

Dislocation loop bias and void swelling in irradiated α-iron from mesoscale and atomistic simulations

Abstract Dislocation loops are ubiquitous in irradiated materials, and dislocation loop bias plays a critical role in void swelling. However, due to complicated interactions between dislocation loops and point defects, it is challenging to evaluate the bias factors of dislocation loops. Here, we determine the bias of sessile < 100 > loops in α-iron using a recently developed atomistic approach based on the lifetime of point defects. We establish a mechanistic understanding of the loop interaction based on the diffusion tendency of point defects near the loop core region. Mobile self-interstitial atoms tend to be absorbed from the edge of the loop, and a trapping region perpendicular to the habit plane of the loop exists. The dislocation loop bias is found to be substantially lower than those of straight dislocations in α-iron and should be included in swelling rate estimates. With the obtained sink strength and bias values, agreement is achieved with experimental results for both absolute values and temperature dependence.

Yu, Ziang↗
DOE OSTI2023

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute–solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Chemistry↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Structural and thermodynamic characterization of CuZr metallic glass nanoparticles: Insights from atomistic simulations

Metallic crystalline nanoparticles (NPs) have been shown to display intriguing size-dependent properties. However, the properties of metallic glass (MG) NPs remain largely unexplored. Here, using molecular dynamics simulations, we produce Cu 64 Zr 36 MG NPs ranging in size from 1 to 20 nm by cooling molten systems at a relatively slow rate of 10 9 K/s. Results indicate that NPs are coated with a Cu-rich layer with a thickness that increases with NP size. We employ various simulated rates of heating, from 10 11 to 10 13 K/s, to determine the melting points of these NPs. The results show a significant decrease in the solidus temperature for NPs smaller than 10 nm. Furthermore, as NP size decreases, the fraction of Cu full icosahedra within the NPs increases, suggesting that smaller NPs are stronger and harder. These findings provide insights for designing heterogeneous metallic nanoglass materials that leverage the size-dependent properties of MG NPs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insights into the Properties of MXenes and MXene Analogs from Atomistic Simulation

We review the properties that have been predicted for MXenes and MXene analogs from computational simulation with a focus on structural and electronic properties, energy storage, and ion transport. Methods considered range from quantum mechanical approaches to classical molecular dynamics. We conclude by reviewing current limitations and outstanding questions for simulation and MXene properties that have been little explored to-date.

Muraleedharan, Murali Gopal↗
DOE OSTI2025

Atomistic Simulations of Thermal and Chemical Expansions of PrNi x Co 1‐x O 3‐δ Accelerated by Machine Learning Potentials

The electrodes and solid-state electrolytes in protonic ceramic electrochemical cells (PCECs) experience significant lattice expansions when exposed to high steam concentrations at elevated temperatures. In this paper, phonon calculations based on a new machine learning potential (MLP) are employed to elucidate the volume expansions of the proton-conducting PrNi x Co 1-x O 3-δ (PNC) lattices, manifested under a combined influence of oxygen vacancies (V$^{\cdot\cdot}_O$ ) and proton uptake (OH$^{\cdot}_O$ ) in the bulk at varying Ni/Co occupancies. It is revealed that the Ni/Co occupancy contributes to thermal and chemical expansions differently, where thermal expansions are related to Co occupancy. In contrast, chemical expansions are more closely associated with the Ni occupancy. Both V$^{\cdot\cdot}_O$ and OH$^{\cdot}_O$ lead to higher thermal expansions when compared to the pristine PNC. The temperature increase will negatively impact the hydration-induced chemical expansions. For combined thermal and chemical expansions, it is predicted that the strategies that boost the PCEC's electrochemical performance may harm the electrode–electrolyte interfacial stability, when the Ni occupancy is high, due to severe chemical expansions. Mitigating chemical expansions of the Ni-abundant PNC will benefit the interfacial stability. Finally, the presented computational methods for phonon calculations, based on emerging machine learning interatomic potential techniques are anticipated to have a lasting impact on future PCEC development.

computational prediction↗
DOE OSTI2022

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗