Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Automated detection of double nuclei galaxies using GOTHIC and the discovery of a large sample of dual AGN

We present a novel algorithm to detect double nuclei galaxies (DNG) called Gothic (Graph-bOosTed iterated HIll Climbing) – that detects whether a given image of a galaxy has two or more closely separated nuclei. Our aim is to test for the presence of dual/multiple active galactic nuclei (AGN) in galaxies that visually represent a DNG. Although galaxy mergers are common, the detection of dual AGN is rare. Their detection is very important as they help us understand the formation of supermassive black hole (SMBH) binaries, SMBH growth and AGN feedback effects in multiple nuclei systems. There is thus a need for an algorithm to do a systematic survey of existing imaging data for the discovery of DNGs and dual AGNs. We have tested Gothic on an established sample of DNGs with a 100 per cent detection rate and subsequently conducted a blind search of 1 million SDSS DR16 galaxies (with spectroscopic data available) lying in the redshift range of z = 0 to 0.75. From the list of candidate DNGs found, we have detected 159 dual AGNs, of which 2 are triple AGN systems. Our results show that dual AGNs are not common, and triple AGN even rarer. The colour (u–r) magnitude plots of the DNGs indicate that star formation is quenched as the nuclei come closer and as the AGN fraction increases. The quenching is especially prominent for dual/triple AGN galaxies that lie at the extreme end of the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analysis of the impact of broad absorption lines on quasar redshift measurements with synthetic observations

ABSTRACT Accurate quasar classifications and redshift measurements are increasingly important to precision cosmology experiments. Broad absorption line (BAL) features are present in 15–20 per cent of all quasars, and these features can introduce systematic redshift errors, and in extreme cases produce misclassifications. We quantitatively investigate the impact of BAL features on quasar classifications and redshift measurements with synthetic spectra that were designed to match observations by the Dark Energy Spectroscopic Instrument (DESI) survey. Over the course of 5 yr, DESI aims to measure spectra for 40 million galaxies and quasars, including nearly three million quasars. Our synthetic quasar spectra match the signal-to-noise ratio and redshift distributions of the first year of DESI observations, and include the same synthetic quasar spectra both with and without BAL features. We demonstrate that masking the locations of the BAL features decreases the redshift errors by about 1 per cent and reduces the number of catastrophic redshift errors by about 80 per cent. We conclude that identifying and masking BAL troughs should be a standard part of the redshift determination step for DESI and other large-scale spectroscopic surveys of quasars.

(galaxies:) quasars: absorption↗
DOE OSTI2025

Landau-level composition of bound exciton states in magnetic field

Here, we present a theory that studies the state composition of a bound exciton in magnetic field. Using a basis set made of products of free electron and hole wave functions in Landau gauge, we derive a secular equation which shows the relation between Landau levels (LLs) of the electron and hole when a bound exciton is formed. Focusing on excitons in the light cone, we establish a scattering selection rule for the interaction of an electron in LL 𝑛 e with a hole in LL 𝑛 h . We solve the resulting secular equation and identify a simple pairing law, 𝑛 e =𝑛 h +𝑙, which informs us on the construction of a bound exciton state with magnetic quantum number 𝑙, and on the interaction of the exciton magnetic moment with magnetic field. We obtain good agreement between theory results and recent measurements of the diamagnetic shifts of exciton states in WSe 2 monolayers.

2-dimensional systems↗
DOE OSTI2022

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2021

Magnetic skyrmion diode: Unidirectional skyrmion motion via symmetry breaking of potential energy barriers

In this work, we realize a magnetic skyrmion diode operated by a unidirectional skyrmion transport that flows in only one direction, which is highly significant for information processing in spintronic and nanoelectronic devices. We easily control the skyrmion transport by engineering asymmetric shapes of geometric structures. Here, we present a simple method to describe the underlying mechanism behind the unidirectional skyrmion transport by characterizing the topography of potential energy surfaces from a purely geometric perspective. Our approach enables a deeper physical insight into skyrmion transport manipulation and efficient design of skyrmion-based devices in geometric structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetoelastic coupling anisotropy in the Kitaev material α-RuCl 3

Here we show that the Kitaev material α-RuCl 3 is among the most prominent candidates to host a quantum spin-liquid state endowed with fractionalized excitations. Recent experimental and theoretical investigations have separately revealed the importance of both the magnetoelastic coupling and the magnetic anisotropy in dependence of the applied magnetic-field direction. In this combined theoretical and experimental research, we investigate the anisotropic magnetic and magnetoelastic properties for magnetic fields applied along the main crystallographic axes as well as for fields canted out of the honeycomb plane. We found that the magnetostriction anisotropy is unusually large compared to the anisotropy of the magnetization, which is related to the strong magnetoelastic $\tilde{Γ}'$-type coupling in our ab-initio derived model. We observed large, nonsymmetric magnetic anisotropy for magnetic fields canted out of the honeycomb ab plane in opposite directions, namely, toward the +c* or -c* axes, respectively. The observed directional anisotropy is explained by considering the relative orientation of the magnetic field with respect to the coaligned RuCl 6 octahedra. Magnetostriction measurements in canted fields support this nonsymmetric magnetic anisotropy; however, these experiments are affected by magnetic torque effects. Comparison of theoretical predictions with experimental findings allow us to recognize the significant contribution of torque effects in experimental setups where α-RuCl 3 is placed in canted magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

One-dimensional Holstein model revisited

Here, we analyze the global ground-state (quantum) phase diagram of the one-dimensional spinful Holstein model at half filling as a function of the strength of the electron-phonon coupling (represented by the strength of the phonon-induced attraction, U) and the phonon frequency ω 0 . In addition to reanalyzing the various asymptotic regimes, we carry out density-matrix renormalization group simulations to correct previous inferences concerning the antiadiabatic (large ω 0 ) and strong-coupling (large U) regimes. There are two distinct phases—a fully gapped commensurate charge-density-wave phase and a spin-gapped Luther-Emery phase with a gapless charge mode—separated by a phase boundary, with a shape that reflects different microscopic physics in the weak- and strong-coupling limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Variation of carrier density in semimetals via short-range correlation: A case study with nickelate NdNiO 2

Carrier density is one of the key controlling factors of material properties, particularly in controlling the essential correlations in strongly correlated materials. Typically, carrier density is externally tuned by doping or gating and remains fixed below room temperature. Strangely, the carrier density in correlated semimetals is often found to vary sensitively against weak external controls such as temperature, magnetic field, and pressure. Here, we develop a realistic simulation scheme that incorporates interatomic noncollinear magnetic correlation without a long-range order. Using the recently discovered nickelate superconductor as an example, we demonstrate a rather generic low-energy mechanism that in semimetals short-range correlation can reversely modulate the carrier density as well. Such a mutual influence between correlation and carrier density provides an extra ingredient for sensitive bifurcating behavior. This special feature of correlated semimetals explains their versatile carrier density at low energy and opens up new possibilities of functionalizing these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Spinon continuum in the Heisenberg quantum chain compound Sr 2 ⁢V 3 ⁢O 9

Magnetic excitations in the spin chain candidate Sr 2 ⁢V 3 ⁢O 9 have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of ~22 meV is observed along the chain formed by alternating magnetic V 4+ and nonmagnetic V 5+ ions, which reveals the importance of the orbital degree of freedom in determining the chain axis as identified by prior electronic structure calculations. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at T N ~ 5.3 K, while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group calculations, and effective field theories confirm Sr 2 ⁢V 3 ⁢O 9 as a host of weakly coupled S = 1/2 chains dominated by antiferromagnetic intrachain interactions of ~7.1(1) meV.

1-dimensional spin chains↗
DOE OSTI2024

Extending the Takagi-Taupin equations for x-ray nanobeam Bragg coherent diffraction

We present an approach for simulating x-ray nanobeam Bragg coherent diffraction patterns based on the Takagi-Taupin equations. Compared to conventional methods, the current approach can be universally applied to any weakly strained system including semi-infinite crystals that diffract dynamically. It addresses issues such as the curved wave front and redivergence of the focused incident beam. In conclusion, we show excellent agreement with experimental data for a strained La 0.7 ⁢Sr 0.3⁢ MnO 3 thin film on a SrTiO 3 substrate and a path to extracting physical information using automatic differentiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical Study of a Dual Representation of the Integer Quantum Hall Transition

Here, we study the critical properties of the noninteracting integer quantum Hall to insulator transition (IQHIT) in a “dual” composite-fermion (CF) representation. A key advantage of the CF representation over electron coordinates is that at criticality CF states are delocalized at all energies. The CF approach thus enables us to study the transition from a new vantage point. Using a lattice representation of CF mean-field theory, we compute the critical and multifractal exponents of the IQHIT. We obtain ν = 2.56 ± 0.02 and η = 0.51 ± 0.01, both of which are consistent with the predictions of the Chalker-Coddington network model formulated in the electron representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗