Results for “first principles simulations” — Page 10

DOE OSTI2023

Atomic-scale polarization switching in wurtzite ferroelectrics

Ferroelectric wurtzites have the potential to revolutionize modern microelectronics because they are easily integrated with multiple mainstream semiconductor platforms. However, the electric fields required to reverse their polarization direction and unlock electronic and optical functions need substantial reduction for operational compatibility with complementary metal-oxide semiconductor (CMOS) electronics. To understand this process, we observed and quantified real-time polarization switching of a representative ferroelectric wurtzite (Al 0.94 B 0.06 N) at the atomic scale with scanning transmission electron microscopy. The analysis revealed a polarization reversal model in which puckered aluminum/boron nitride rings in the wurtzite basal planes gradually flatten and adopt a transient nonpolar geometry. Independent first-principles simulations reveal the details and energetics of the reversal process through an antipolar phase. Furthermore, this model and local mechanistic understanding are a critical initial step for property engineering efforts in this emerging material class.

36 MATERIALS SCIENCE↗

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 OSTI2023

Machine learning and LHC event generation

First-principle simulations are at the heart of the high-energy physics research program. They link the vast data output of multi-purpose detectors with fundamental theory predictions and interpretation. This review illustrates a wide range of applications of modern machine learning to event generation and simulation-based inference, including conceptional developments driven by the specific requirements of particle physics. New ideas and tools developed at the interface of particle physics and machine learning will improve the speed and precision of forward simulations, handle the complexity of collision data, and enhance inference as an inverse simulation problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Gallium Nitride Superjunction Transistor (Continued Funding Report)

The initial goal of this project is to develop the gallium vacancy assisted diffusion (GAID) method in GaN and apply it to demonstrate 2-terminal, 3-terminal, and large area Superjunction-enhanced modules. These devices will find use in high efficiency power conversion applications with the capability of operating at higher voltages, higher frequency, and less loss than comparable stat-ofthe-art devices. The approach of this project has been to utilize simulation tools to drive diffusion experiments and inform semiconductor device fabrication processes. VASP has been relied upon to simulate first principles diffusion mechanisms while SILVACO TCAD has been utilized to understand tradeoffs of the the critical dimensions in GaN devices.

42 ENGINEERING↗

DOE OSTI2023

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides (Final Report)

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides. Final Report

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

First Principles Modeling of Tritium Trapping in γ-LiAlO 2 Nanovoids

Experimental STEM and APT investigations both indicate the presence of nanovoids in γ-LiAlO 2 pellets. While those nanovoids can affect the diffusion of tritium in pellets, first principles simulations have been conducted to evaluate the potential for tritium trapping by single and double cation and anion vacancies. It was found that tritium trapping allows to recover a significant amount of the formation energy cost of cation vacancy, while the trapping of tritium by anion vacancies further increase the energy cost of forming a vacancy. The simulations also indicated the potential for tritium to be trapped in various forms by forming bonds such as Al—T, O—T or T 2 . These results raise further questions regarding the speciation and dynamics of tritium trapped in nanovoids.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

Revolutionizing Materials Design: The Intersection of Quantum Mechanics and Data Modeling

The field of materials design is currently experiencing a notable evolution, driven by the convergence of sophisticated computational methodologies based on first principles and data-driven modeling approaches. I will review our recent endeavors employing AI/ML to expedite first-principles simulations and mitigate traditional methods' temporal and spatial limitations. Central to our efforts is developing and utilizing ML interatomic potentials (MLPs) across a diverse spectrum of materials. We show that MLPs serve as invaluable tools for navigating the complexities of the simulations, such as understanding the behavior of MgO at extreme environments of ~1 terapascal and temperatures >10,000 Kelvin. Moreover, we show that MLPs can provide precise details of the intricate dynamics governing the oxidation processes of binary alloy systems due to the competition between surface segregation and reconstruction tendencies. In summation, advancements in MLPs open the door to fresh possibilities in material modeling and, ultimately, discovery.

Saidi, Wissam↗

DOE OSTI2022

A Novel Dynamic Polymer Synthesis via Chlorinated Solvent Quenched Depolymerization

Dynamic polymers with both physical interactions and dynamic covalent bonds exhibit superior performance while achieving such dry polymers in an efficient manner remains a challenge. Herein, we report a novel organic solvent quenched polymer synthesis using a natural molecule, i.e., thioctic acid (TA), which has both dynamic disulfide bond and carboxylic acid. The effects of the solvent types and concentrations along with reaction times on the proposed reaction have been thoroughly explored for polymer synthesis. Solid-state proton nuclear magnetic resonance ( 1 H NMR) and first-principles simulations are carried out to investigate the reaction mechanism. It shows that the chlorinated solvent can efficiently stabilize and mediate the depolymerization of poly(TA) which is more kinetically favorable upon lowering the temperature. Attributed to the numerous dynamic covalent disulfide bonds and noncovalent hydrogen bonding, the obtained poly(TA) shows high extensibility, self-healable and re-processable properties. It can also be employed as an efficient adhesive even on a Teflon surface and 3D printed using the fused deposition modeling (FDM) technique. This new polymer synthesis approach on using organic solvents as catalysts along with the unique reaction mechanism provides a new pathway for efficient polymer synthesis, especially for those multi-functional dynamic polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Ultrafast Martensitic Phase Transition Driven by Intense Terahertz Pulses

We report on an ultrafast nonequilibrium phase transition with a strikingly long-lived martensitic anomaly driven by above-threshold single-cycle terahertz pulses with a peak field of more than 1 MV/cm. A nonthermal, terahertz-induced depletion of low-frequency conductivity in Nb 3 Sn indicates increased gap splitting of high-energy 12 bands by removal of their degeneracies, which induces the martensitic phase above their equilibrium transition temperature. In contrast, optical pumping leads to a Γ 12 gap thermal melting. Such light-induced nonequilibrium martensitic phase exhibits a substantially enhanced critical temperature up to ~100 K, i.e., more than twice the equilibrium temperature, and can be stabilized beyond technologically relevant, nanosecond time scales. Together with first-principle simulations, we identify a compelling terahertz tuning mechanism of structural order via Γ 12 phonons to achieve the ultrafast phase transition to a metastable electronic state out of equilibrium at high temperatures far exceeding those for equilibrium states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Striped Jets in Post–Neutron Star Merger Systems

Abstract Models invoking magnetic reconnection as the particle acceleration mechanism within relativistic jets often adopt a gradual energy dissipation profile within the jet. However, such a profile has yet to be reproduced in first-principles simulations. Here we perform a suite of 3D general relativistic magnetohydrodynamic simulations of post–neutron star merger disks with an initially purely toroidal magnetic field. We explore the variations in both the microphysics (e.g., nuclear recombination, neutrino emission) and system parameters (e.g, disk mass). In all of our simulations, we find the formation of magnetically striped jets. The stripes result from the reversals in the poloidal magnetic flux polarity generated in the accretion disk. The simulations display large variations in the distributions of stripe duration, τ , and power, 〈 P Φ 〉. We find that more massive disks produce more powerful stripes, the most powerful of which reaches 〈 P Φ 〉 ∼ 10 49 erg s −1 at τ ∼ 20 ms. The power and variability that result from the magnetic reconnection of the stripes agree with those inferred in short-duration gamma-ray bursts. We find that the dissipation profile of the cumulative energy is roughly a power law in both radial distance, z , and τ , with a slope in the range of ∼1.7–3; more massive disks display larger slopes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Using Diffuse Scattering to Observe X-Ray-Driven Nonthermal Melting

We present results from the SPring-8 Angstrom Compact free electron LAser (SACLA) XFEL facility, using a high intensity ($\sim\!10^{20}\,$W/cm$^2$) X-ray pump X-ray probe scheme to observe changes in the ionic structure of silicon induced by X-ray heating of the electrons. By avoiding Laue spots in the scattering signal from a single crystalline sample, we observe a rapid rise in diffuse scattering, which we attribute to a loss of lattice order and a transition to a liquid state within 100 fs of irradiation, a timescale which agrees well with first principles simulations, but is faster than that predicted by purely inertial behavior. This method is capable of observing liquid scattering without masking or filtering of signal from the ambient solid, allowing the liquid structure to be measured throughout and beyond the phase change.

Hartley, Nicholas J.↗

DOE OSTI2023

Temperature Dependence of Band Gap Renormalization in High-T Sensor Materials via First-Principles and Experimental Corroboration

Understanding the temperature dependence of functional properties of high-T gas sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we assess the temperature dependence of band gap renormalization in metal oxides and perovskites by employing Allen-Heine-Cardona theory with first-principles simulations and corroborate with experimental observation. The calculated temperature-dependent band gap changes of these materials studied are in good agreement with in-house experimental data, proving that the theory can adequately predict renormalization on the band gap in the system of interest. The predicted and measured band gap variations are characterized using an analytical model, which can provide useful insights on the simulated zero-temperature band gaps. Based on the available data, a set of 53 metal oxides and perovskites were identified as potential high-T gas sensors. A machine learning model has been developed to predict the band-gap change by capturing the overall trend of the empirical parameters with respect to a reduced feature obtained by transforming the set of available physical features.

Park, Jongwoo↗

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↗