SEARCH · Engineering Papers
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.
electron density profile as a function of normalized poloidal flux: psinorm ne(10^20/m^3) dne/dpsiN Pressure profile as a function of normalized poloidal flux: psinorm ptot(KPa) dptot/dpsiN scaled to oftain beta=5% at magnetic axis Magnetic safety factor profile as a function of normalized poloidal flux: psinorm q_profile All needed to perform the NOVA calculations.
Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.
Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.
Here, we use high-resolution angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) to investigate the electronic structure of trigonal phase PtBi 2 (𝑡 −PtBi 2 ), a proposed Weyl semimetal that is expected to exhibit topological Fermi arcs. Our ARPES data elucidates the topography of these objects and confirms their Fermi arc character. The arcs are formed by surface bands that have fairly flat bottoms located very close to the chemical potential ∼ 6–8 meV, before they merge with bulk bands at higher binding energy. Comparison of the ARPES data with DFT calculations shows good agreement about their location and topography. Data acquired at low temperatures does not show any signatures of superconductivity down to 3 K in terms of expected changes in dispersion due to formation of Bogoliubov quasiparticles or superconducting gap in the form of reliable shifts of peaks in energy distribution curves.
Recent studies have focused on the relationship between charge density wave (CDW) collective electronic ground states and nontrivial topological states. YTe 3 , a nonmagnetic quasi-two-dimensional chalcogenide, has been reported to exhibit a CDW state below 334 K. Using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT), we establish that YTe 3 is a CDW-induced Kramers nodal line (KNL) metal, a recently proposed topological state of matter. Scanning tunneling microscopy and low energy electron diffraction reveal two orthogonal domains, each with a unidirectional CDW and a similar wave vector (𝐪 CDW ). When the influence of twin domains is considered, the effective band structure (EBS) computations that utilize DFT-calculated bands using a noncentrosymmetric structure determined by x-ray crystallography, show excellent agreement with ARPES. The noncentrosymmetry of YTe 3 is established by Raman spectroscopy. The Fermi surface and ARPES intensity plots show weak shadow bands displaced by 𝐪 CDW from the main bands. Furthermore, these are linked to CDW modulation, as the EBS calculation confirms. Bilayer split main and shadow bands suggest the existence of crossings, according to theory and experiment. DFT bands, including spin-orbit coupling, indicate existence of a KNL along the Σ direction from multiple crossings of bands dispersing perpendicular to it. Additionally, doubly degenerate bands are only found along the KNL at all energies, with some bands dispersing through the Fermi level.
Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.
The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.
Here, we include spin–orbit coupling (SOC) effects in linearized pair-density functional theory (L-PDFT), which is a multistate extension of multiconfiguration pair-density functional theory (MC-PDFT). Both 1-electron and 2-electron SOC integrals are computed using Breit-Pauli and Douglas–Kroll–Hess Hamiltonians in the atomic mean-field approximation. SO-L-PDFT removes the unphysical J-symmetry breaking observed in MC-PDFT. The accuracy of SO-L-PDFT is validated by calculations of zero-field splittings, fine-structure excitation energies, and low-energy excited-state spectra for a diverse group of atoms and molecules spanning the whole range of the periodic table, including atoms of groups 3, 11, and 13–17, the Ce 3+ and U 5+ ions, group 16 monohydrides, group 17 monoxides, lanthanide hexachlorides ([CeCl 6 ] 3− , [PrCl 6 ] 3− , and [NdCl 6 ] 3− ), actinyl ions ([UO 2 ] + , [NpO 2 ] 2+ ), and tricarbonatoactinyl complexes ([UO 2 (CO 3 ) 3 ] 5− , [NpO 2 (CO 3 ) 3 ] 4− ). We also compare the results to new spin–orbit-inclusive calculations by single-state and multistate multireference perturbation theory.
X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this context, a key property is given by the Rayleigh weight that describes the electronic localization around the ions. Here, we show that it is possible to extract the Rayleigh weight directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility [Döppner et al., Nature 618, 270–275 (2023)]. We demonstrate that experimental results for the Rayleigh weight open up new avenues for the interpretation of XRTS experiments by matching the measurement with ab initio simulations such as density functional theory or path integral Monte Carlo. Interestingly, this new procedure leads to significantly lower density compared to previously used chemical models.
The large van der Waals gap in transition metal dichalcogenides (TMDs) offers an avenue to tune the ground state of 2D materials through the intercalation of magnetic atoms. Here, we investigate the charge correlations in Fe 1/3 TaS 2 , Co 1/3 TaS 2 , and Fe 0.35 NbS 2 by combining angle-resolved photoemission spectroscopy (ARPES), x-ray scattering, magnetometry, and density functional theory (DFT). We find that, while short-range charge fluctuations develop in Ta-based compounds, Fe 0.35 NbS 2 exhibits long-range charge order which is strongly coupled with magnetic order and tunable by external magnetic field. Our electronic structure analysis reveals that intercalation reconstructs the Fermi surface via charge transfer and band renormalization, yet does not generate the nesting conditions compatible with the observed ordering vectors. Complementary phonon calculations further exclude a conventional electron-phonon origin of charge order. Together, these results establish magnetoelastic coupling as the dominant mechanism behind charge ordering in Fe 0.35 NbS 2 and highlight the contrasting role of Nb and Ta hosts in stabilizing correlated ground states in intercalated TMDs.
High-quality, multifunctional two-dimensional (2D) titanium oxynitide (TiNO) thin films and one-dimensional (1D) TiNO nanowires have been synthesized using a pulsed laser deposition, a simple, fast, and congruent evaporation method. First-principles calculations as a function of surface orientation and termination indicate that surface oxidation of TiNO nanowires can stabilize the (110) orientation observed experimentally. The specific capacitance value for the TiNO nanowire samples (2725 mF/cm 2 ) has been found to be nearly six times more than that of the TiNO thin film samples (400 mF/cm 2 ), which is attributed to the high packing density of TiNO nanowires over a given area. The nanowire samples have also been found to exhibit a significantly higher energy density (1.35 μWh/cm 2 ) than the TiNO thin-film samples (0.33 μWh/cm 2 ). Thus, the TiNO material system in thin-film and nanowire forms has been demonstrated to be a promising candidate for use as an electrode material in supercapacitors and other charge-storage applications.
We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.
Here, we predict two different universal scaling regimes for the quantum transmission of metallic nanodevices following the addition of a small amount of uncorrelated disorder. A nanodevice is connected to two thin semi-infinite uniform leads, and the Non-Equilibrium Green’s Function (NEGF) methodology yields the electron transmission $\mathscr{T}(E)$ as a function of the injected electron energy $E$. Ballistic nanodevices have no disorder and have $\mathscr{T}(E)=1$for all $E$ that allow electron propagation in the leads. Quantum dragon nanodevices can have extremely strong properly correlated disorder, and still have $\mathscr{T}(E)=1$for all $E$. Additional uncorrelated site disorder leads to Fano resonances in $\mathscr{T}(E)$. Averaging over the uncorrelated disorder we predict using perturbation theory two universal scaling regimes for $\mathscr{T}_{ave}(E)$. The functional form of both universal scaling regimes depend on the device length and width, energy, and variance of the uncorrelated disorder. The second scaling regime, valid for small but somewhat larger uncorrelated disorder than the first scaling regime, also has the form dependent on the density of states of the system. These two scaling regimes are demonstrated to be valid via large scale computer calculations.
Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.
We have investigated the electronic and magnetic structures of topological kagome Fe 1−𝑥 Mn 𝑥 Sn (0 ≤ 𝑥 ≤ 0.3) thin films via neutron diffraction, electronic transport measurements, and ab initio density functional theory (DFT) to understand the interplay between hole doping, magnetism, and the electronic structures. Temperature-dependent neutron diffraction measurements on parent FeSn reveal the Néel temperature to be 𝑇 N ∼ 355 K and the underlying A-type antiferromagnetic ordering is associated with a wave vector 𝒒 = (001/2). Upon Mn doping to 𝑥 = 0.15, 𝑇 N decreases slightly while the magnetic ordering vector remains the same. Resistivity measurements show metallic characteristics and in-plane anisotropy down to 10 K for all the investigated samples. The effects of hole doping are mapped in terms of electronic ground state calculations via DFT which show that the Dirac point is moved closer to the Fermi level (𝐸 F ) and the flat bands get pushed away from 𝐸 F upon hole doping. However, a comparison between hole-doped Fe 1−𝑥 Mn 𝑥 Sn and electron-doped Fe 1−𝑥 Co 𝑥 Sn indicates that the Néel temperature does not scale with the position of 𝐸 F relative to the flat band. Furthermore, our results establish the antiferromagnetic state of FeSn and Fe 1−𝑥 Mn 𝑥 Sn films at room temperature, laying the groundwork for future studies of magnetism in kagome heterostructures.
The recently synthesized monolayer fullerene network in a quasi-hexagonal phase (qHP-C 60 ) exhibits superior electron mobility and optoelectronic properties compared to molecular fullerene (C 60 ), making it highly promising for a variety of applications. However, the microscopic carrier dynamics of qHP-C 60 remain unclear, particularly in realistic environments, which are of significant importance for applications in optoelectronic devices. Unfortunately, traditional ab initio methods are prohibitive for capturing the real-time carrier dynamics of such large systems due to their high computational cost. In this work, we present the first real-time electron–nuclear dynamics study of qHP-C 60 using velocity-gauge density functional tight binding, which enables us to perform several picoseconds of excited-state electron–nuclear dynamics simulations for nanoscale systems with periodic boundary conditions. When applied to C 60 , qHP-C 60 , and their solvated counterparts, we demonstrate that water/moisture significantly increases the electron–hole recombination time in C 60 but has little impact on qHP-C 60 . Our excited-state electron–nuclear dynamics calculations show that qHP-C 60 is extremely unique and enable exploration of time-resolved dynamics for understanding excited-state processes of large systems in complex, solvated environments.
A quantum error-correcting code with a nonzero error threshold undergoes a mixed-state phase transition when the error rate reaches that threshold. We explore this phase transition for Haar-random quantum codes, in which the logical information is encoded in a random subspace of the physical Hilbert space. We focus on the spectrum of the encoded system density matrix as a function of the rate of uncorrelated, single-qudit errors. For low error rates, this spectrum consists of well-separated bands, representing errors of different weights. As the error rate increases, the bands for high-weight errors merge. The evolution of these bands with increasing error rate is well described by a simple analytic ansatz. Using this ansatz, as well as an explicit calculation, we show that the threshold for Haar-random quantum codes saturates the hashing bound, and thus coincides with that for random stabilizer codes. For error rates that exceed the hashing bound, typical errors are uncorrectable, but postselected error correction remains possible until a much higher detection threshold. Postselection can in principle be implemented by projecting onto subspaces corresponding to low-weight errors, which remain correctable past the hashing bound.
Metastable states of matter are of great interest as they offer the promise of novel functionality. They are often a natural consequence of exposure to nonequilibrium environments. In oxides, metastability can take the form of new polymorphs, chemical disorder, and even amorphization. While significant attention has been given to the impact those changes have on the atomic properties of the material, the corresponding changes in the electronic structure have received less attention. Here, using density functional theory, we consider how the electronic structure varies with potential metastable structures in two classes of lanthanide-bearing oxides—pyrochlores and interlanthanide sesquioxides. We find that the changes depend strongly on both the crystal structure and crystal chemistry of the compound with, for example, disordering and amorphization either increasing or decreasing the bandgap depending on the chemistry. For the 𝐴 2 𝐵 2 O 7 pyrochlores, we find different dependencies of the bandgap on the 𝐴 = 𝐿𝑛 cations as the 𝐵 cation is changed, which we relate to the nature of the density of states at the conduction band minimum for different 𝐵 chemistries. Our calculations are validated by electron energy loss spectroscopy measurements for two pyrochlore compounds in which amorphization does reduce the bandgap, consistent with our calculations on these two compounds. In conclusion, our results highlight the relationship between atomic and electronic structure and how radiation can be used to modify and potentially control the electronic properties of oxides.