Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles simulations”

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 OSTI2025

Unconventional gate-induced superconductivity in transition-metal dichalcogenides

Superconductivity in few-layer semiconducting transition-metal dichalcogenides (TMDs) can be induced by field-effect doping through ionic-liquid gating. While several experimental observations have been collected over the years, a fully consistent theoretical picture is still missing. Here we develop a realistic framework that combines the predictive power of first-principles simulations with the versatility and insight of Bardeen-Cooper-Schrieffer gap equations to rationalize such experiments. The multivalley nature of semiconducting TMDs is taken into account, together with the doping- and momentum-dependent electron-phonon and Coulomb interactions. Consistently with experiments, we find that superconductivity occurs when the electron density is large enough that the 𝑄 valleys get occupied, as a result of a large enhancement of electron-phonon interactions. Despite being phonon driven, the superconducting state is predicted to be sensitive to Coulomb interactions, which can lead to the appearance of a relative sign difference between valleys and thus to a 𝑠 +− character. We discuss qualitatively how such scenario may account for many of the observed physical phenomena for which no microscopic explanation has been found so far, including in particular the presence of a large subgap density of states, and the sample-dependent dome-shaped dependence of 𝑇 𝑐 on accumulated electron density. Our results provide a comprehensive analysis of gate-induced superconductivity in semiconducting TMDs, and introduce an approach that will likely be valuable for other multivalley electronic systems, in which superconductivity occurs at relatively low electron density.

2-dimensional systems↗
DOE OSTI2026

Isostructural electronic transition in MoS 2 probed by solid-state high-harmonic generation spectroscopy

Studying materials under extreme pressure in diamond anvil cells (DACs) is key to discovering emergent states of matter, yet no method currently allows the direct measurement of the electronic structure in this environment. Solid-state high-harmonic generation (sHHG) offers a unique all-optical window into the electronic structure of materials. We demonstrate sHHG spectroscopy inside a DAC by probing 2H-MoS 2 , up to 30 GPa, revealing a pressure-induced crossover of the lowest direct bandgap from the K-point to the Γ-point. This transition manifests as a sharp minimum in harmonic intensity and a 30° rotation of the sHHG polarization anisotropy, despite the absence of a structural phase change. First-principles simulations attribute these features to interference between competing excitation pathways at distinct points in the Brillouin zone. Our results establish sHHG as a sensitive probe of electronic transitions at high pressure, enabling access to quantum phenomena that evade detection by conventional techniques.

Nebgen, Bailey R. [University of California, Berke↗
DOE OSTI2026

Optimizing spin qubit coherence through materials codesign

The evolution of defect-based spin qubit systems is currently transitioning from fundamental studies and proof-of-concept demonstrations into applications in the burgeoning field of quantum technology. Within this context, new challenges emerge, in particular, the need to understand and engineer the fundamental materials that form the hardware building blocks critical for the scalability and wide-scale adoption of such technologies. While earlier discussions have often focused on qubits within idealized systems, major limitations on spin coherence and optical properties arise from effects imposed by the nonideality of the surrounding host matrix. Decoherence can stem from a variety of sources, including other qubits, nuclear spins, and parasitic point- and extended defects, which interact with the qubit via magnetic and electric fields, photons, phonons, and strain. In this article, we focus on the relevant sources and mechanisms through which decoherence occurs and provide potential mitigation strategies via the synergistic integration of first-principles simulations and materials synthesis and engineering. We aim to provide a tangible link between material properties and material functions thereby enabling materials-by-design.

ab initio simulations↗
DOE OSTI2025

Effects of Minor Alloying Elements and Impurities on Long-Range Ordering in Ni-Cr-Based Alloys

There is growing concern about long-range ordered Ni2Cr phase forming during thermal aging of Ni-Cr-based alloys (e.g., 690, 625), which leads to lattice contraction, hardening, and embrittlement, possibly limiting the service life of light water reactor components. The timespan when this phase transformation occurs can vary in orders of magnitude between model and commercial alloys and is likely related to the presence of major (e.g., Fe, Mo) and minor alloying elements (e.g., Ti, Si, Nb). These minor alloying elements have been hypothesized to bind with vacancies and therefore limit diffusion. Contrarily, impurities (e.g., P, S) are believed to enhance vacancy diffusion and accelerate the Ni2Cr phase transformation because of a vacancy drag mechanism. This study implements experimental isothermal aging of model Ni-Cr binary and Ni-Cr-X ternary alloys, where X is Ti, Si, Nb, and P to help understand how each minor alloying element affects the Ni2Cr phase transformation timeline. Synchrotron-based X-ray diffraction and hardness measurements are used to quantify the evolution of the Ni2Cr phase transformation through a variety of short aging times to identify the onset of long-range order. These experimental results are compared to first principles simulations from literature to understand the mechanism by which each element impacts vacancy diffusion and the formation of the long-range ordered phase, Ni2Cr.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

First Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheets membranes are promising candidates for ammonia separation from nitrogen and hydrogen, yet questions remain on the origin of their high selectivity. Silanols, Si-OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Silanols can chemically react with ammonia, which may further contribute to the ability of zeolite nanosheet membranes to separate it from other gases. In this work, we use first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We find that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explores H-bonding and proton transfer within this system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First-Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheet membranes are promising candidates for ammonia separation from nitrogen and hydrogen gases, yet questions remain on the origin of their high selectivity. Silanols, Si–OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Acidic silanol groups can chemically react with ammonia via transfer of a proton to form ammonium ions, which may further contribute to the ability of zeolite nanosheet membranes to separate ammonia from other gases. In this work, we used first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We found that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explored H-bonding and proton transfer within this system.

ammonia↗
DOE OSTI2025

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

Hall2De Simulations with a First-Principles Electron Transport Model Based on the Electron Cyclotron Drift Instability

Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.

Mikellides, Ioannis G.↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

First Principles Based Reactive Atomistic Simulations to Understand the Effects of Molecular Hypervelocity Impact on Cassini's Ion and Neutral Mass Spectrometer

We report here on the predicted impact of species such as ice-water, CO2, CH4, and NH3, on oxidized titanium, as well as HC species on diamond surfaces. These simulations provide the dynamics of product distributions during and after a hypervelocity impact event, ionization fractions, and dissociation probabilities for the various species of interest as a function of impact velocity (energy). We are using these results to determine the relevance of the fragmentation process to Cassini INMS results, and to quantify its effects on the observed spectra.

hypervelocity collisions↗
DOE OSTI2024

Exponential Improvements in the Simulation of Lattice Gauge Theories Using Near-Optimal Techniques

We report a first-of-its-kind analysis on post-Trotter simulation of U(1), SU(2), and SU(3) lattice gauge theories including fermions in arbitrary spatial dimension. We provide explicit circuit constructions as well as T-gate counts and logical qubit counts for Hamiltonian simulation. We find a reduction of up to 25 orders of magnitude in space-time volume over Trotter methods for simulations of non-Abelian lattice gauge theories relevant to the standard model. This improvement results from our algorithm having polynomial scaling with the number of colors in the gauge theory, achieved by utilizing oracle constructions relying on the sparsity of physical operators, in contrast to the exponential scaling seen in state-of-the-art Trotter methods, which employ explicit mappings onto Pauli operators. Our work demonstrates that the use of advanced algorithmic techniques leads to dramatic reductions in the cost of simulating fundamental interactions, bringing it in step with resources required for first-principles quantum simulation of chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Age hardening response of Cu-30Ni alloys: The role of Si microalloying additions

Cu-30Ni alloys offer high strength and excellent corrosion resistance for maritime applications. While primarily a solid solution system, industrial alloys typically include microalloying additions of Fe, Mn, Nb and/or Si to enable age hardening. However, an understanding of the microstructural evolution during age hardening remains incomplete. Here, we conduct systematic ageing treatments and report the resulting hardness and microstructures. The Cu-30Ni alloy with Si aged at 650 °C for 6 h demonstrates significantly enhanced Vickers hardness, reaching ∼201 HV 1 compared to 103 HV 1 for the as-homogenised sample. Si-rich clusters and γ′ precipitates are identified, with their composition, size, and volume fraction quantified to determine their strengthening contributions. First-principles atomistic simulations elucidate the underlying formation mechanisms of these clusters and precipitates, highlighting the critical role of Si in driving their nucleation and growth. This study advances the design of high-strength, Cu alloys with the potential for corrosion resistance in demanding maritime environments.

Age hardening↗
NASA NTRS2008

Rotor Airloads Prediction Using Unstructured Meshes and Loose CFD/CSD Coupling

The FUN3D unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been modified to allow prediction of trimmed rotorcraft airloads. The trim of the rotorcraft and the aeroelastic deformation of the rotor blades are accounted for via loose coupling with the CAMRAD II rotorcraft computational structural dynamics code. The set of codes is used to analyze the HART-II Baseline, Minimum Noise and Minimum Vibration test conditions. The loose coupling approach is found to be stable and convergent for the cases considered. Comparison of the resulting airloads and structural deformations with experimentally measured data is presented. The effect of grid resolution and temporal accuracy is examined. Rotorcraft airloads prediction presents a very substantial challenge for Computational Fluid Dynamics (CFD). Not only must the unsteady nature of the flow be accurately modeled, but since most rotorcraft blades are not structurally stiff, an accurate simulation must account for the blade structural dynamics. In addition, trim of the rotorcraft to desired thrust and moment targets depends on both aerodynamic loads and structural deformation, and vice versa. Further, interaction of the fuselage with the rotor flow field can be important, so that relative motion between the blades and the fuselage must be accommodated. Thus a complete simulation requires coupled aerodynamics, structures and trim, with the ability to model geometrically complex configurations. NASA has recently initiated a Subsonic Rotary Wing (SRW) Project under the overall Fundamental Aeronautics Program. Within the context of SRW are efforts aimed at furthering the state of the art of high-fidelity rotorcraft flow simulations, using both structured and unstructured meshes. Structured-mesh solvers have an advantage in computation speed, but even though remarkably complex configurations may be accommodated using the overset grid approach, generation of complex structured-mesh systems can require months to set up. As a result, many rotorcraft simulations using structured-grid CFD neglect the fuselage. On the other hand, unstructured-mesh solvers are easily able to handle complex geometries, but suffer from slower execution speed. However, advances in both computer hardware and CFD algorithms have made previously state-of-the-art computations routine for unstructured-mesh solvers, so that rotorcraft simulations using unstructured grids are now viable. The aim of the present work is to develop a first principles rotorcraft simulation tool based on an unstructured CFD solver.

Biedron, Robert T.↗
NASA NTRS1990

First-principles calculations of shear moduli for Monte Carlo-simulated Coulomb solids

The paper presents a first-principles study of the shear modulus tensor for perfect and imperfect Coulomb solids. Allowance is made for the effects of thermal fluctuations for temperatures up to the melting conditions. The present theory treats the cases of the long-range Coulomb interaction, where volume fluctuations should be avoided in the Ewald sums.

Ogata, Shuji↗
DOE OSTI2025

Metastability and Ostwald step rule in the crystallisation of diamond and graphite from molten carbon

Abstract Experimental challenges in determining the phase diagram of carbon at temperatures and pressures near the graphite-diamond-liquid triple point are often related to the persistence of metastable crystalline or glassy phases, superheated crystals, or supercooled liquids. A deeper understanding of the crystallisation kinetics of diamond and graphite is crucial for effectively interpreting the outcomes of these experiments. Here, we reveal the microscopic mechanisms of diamond and graphite nucleation from liquid carbon through molecular simulations with first-principles machine learning potentials. Our simulations accurately reproduce the experimental phase diagram of carbon near the triple point and show that liquid carbon crystallises spontaneously upon cooling. Metastable graphite crystallises in the domain of diamond thermodynamic stability at pressures above the triple point. Furthermore, whereas diamond crystallises through a classical nucleation pathway, graphite follows a two-step process in which low-density fluctuations forego ordering. Calculations of the nucleation rates of the two competing phases confirm this result and reveal a manifestation of Ostwald’s step rule, where the strong metastability of graphite hinders the transformation to the stable diamond phase. Our results provide a key to interpreting melting and recrystallisation experiments and shed light on nucleation kinetics in polymorphic materials with deep metastable states.

Science & Technology - Other Topics↗
NASA NTRS2008

Simulating Coronal Mass Ejections in the Heliosphere

Coronal mass ejections (CMEs) and their attendant eruptive flares (EFs) are the most explosive events in the solar system. Magnetic reconnection, arguably the most universal of heliophysical processes, is widely accepted as playing crucial roles in these phenomena. In this paper we summarize some recent advances in understanding CMEs and EFs developed from first-principles numerical simulations of slowly evolving configurations that transition suddenly to violent eruptions. The study emphasizes the role of magnetic complexity in the solar sources of the most energetic events, which can be exploited by magnetic reconnection to unleash the mass, magnetic flux, and energies of the CME and the EF on the Heliosphere

DeVore, Richard C.↗
NASA NTRS2019

Hybrid RANS/LES Simulation of Jet Surface Interaction Noise

NASA’s Commercial Supersonic Technology (CST) project has formulated a technical challenge to design a quiet propulsion system for a low boom supersonic aircraft that meets Federal Aviation Authority’s airport noise regulations with sufficient margin. Several proposed configurations take advantage of shielding from the wing or other air-frame components. Development of carefully validated computational tools are necessary for critically evaluating installation concepts that are currently being proposed to meet the technical challenge. Semi-empirical models that predict the noise reduction potential of arbitrary shielding surfaces are yet to mature. Another key challenge is the systematic assessment of additional noise from the interaction between high speed jet turbulence and a surface in it’s vicinity. As a first step towards predicting noise reduction due to radical installation concepts from first principles, we simulate the noise generated by a high speed turbulent round jet near a simple planar surface. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center. Sensitivity of far-field noise predictions to grid resolution is systematically documented. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing both the jet and the shielding surface is used to predict far-field noise from the simulated flowfield. Details of the structured overset grids, numerical discretization, and turbulence model are provided. Near-field comparisons to PIV data and far-field comparisons to microphone array measurements are discussed. Excellent agreement for an initial validation study on an isolated free round jet was obtained and the findings were utilized in the jet surface interaction study. The split between shielded and reflected side of the microphone array was captured with good agreement, as well as the peak in the noise spectra due to scattering of turbulent energy into sound by the trailing edge of the surface.

Stich, Gerrit-Daniel↗