Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field model”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Singlet, triplet, and pair density wave superconductivity in the doped triangular-lattice moiré system

Recent experimental progress has established the twisted bilayer transition metal dichalcogenide (TMD) as a highly tunable platform for studying many-body physics. Particularly, the homobilayer TMDs under displacement field are believed to be described by a generalized triangular-lattice Hubbard model with a spin-dependent hopping phase θ. To explore the effects of θ on the system, we perform density matrix renormalization group calculations for the relevant triangular lattice t-J model. By changing θ at small hole doping, we obtain a region of quasi-long-range superconducting order coexisting with charge and spin density wave within 0 < θ < π/3. The superconductivity is composed of a dominant spin singlet d-wave and a subdominant triplet p-wave pairing. Intriguingly, the S z = ±1 triplet pairing components feature pair-density waves. In addition, we find a region of triplet superconductivity coexisting with charge-density wave and ferromagnetism within π/3 < θ < 2π/3, which is related to the former phase at smaller θ by a combined operation of spin-flip and gauge transformation. Furthermore, our findings provide insights and directions for experimental search for exotic superconductivity in twisted TMD systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Sachdev-Ye-Kitaev models and beyond: Window into non-Fermi liquids

In this work is a review of the Sachdev-Ye-Kitaev (SYK) model of compressible quantum many-body systems without quasiparticle excitations, and its connections to various theoretical studies of non-Fermi liquids in condensed matter physics. The review is placed in the context of numerous experimental observations on correlated electron materials. Strong correlations in metals are often associated with their proximity to a Mott transition to an insulator created by the local Coulomb repulsion between the electrons. The phase diagrams of a number of models of such a local electronic correlation are explored, employing a dynamical mean-field theory in the presence of random spin exchange interactions. Numerical analyses and analytical solutions, using renormalization group methods and expansions in large spin degeneracy, lead to critical regions that display SYK physics. The models studied include the single-band Hubbard model, the t–J model, and the two-band Kondo-Heisenberg model in the presence of random spin exchange interactions. Also examined are non-Fermi liquids obtained by considering each SYK model with random four-fermion interactions to be a multiorbital atom, with the SYK atoms arranged in an infinite lattice. Connections are made to theories of sharp Fermi surfaces without any low-energy quasiparticles in the absence of spatial disorder, obtained by coupling a Fermi liquid to a gapless boson; a systematic large-N theory of such a critical Fermi surface, with SYK characteristics, is obtained by averaging over an ensemble of theories with random boson-fermion couplings. Finally, an overview of the links between the SYK model and quantum gravity is presented, and the review ends with an outlook on open questions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Divergence of Velocity Fields in Electrochemical Systems

The passage of current through electrochemical systems results in the development of concentration gradients in the electrolytic phase that can be modeled using concentrated solution theory. Application of this theory requires knowledge of three concentration-dependent transport coefficients, which are often taken to be conductivity, diffusion coefficient, and the cation transference number with respect to the solvent velocity. The governing diffusion equation for molar concentration contains two additional terms - the thermodynamic factor which is related to activity of the electrolytic species and the solvent velocity. The main advance in this paper is the derivation of an expression for the divergence of the solvent velocity. Solving this equation requires knowledge of the partial molar volume of the electrolyte. Analogous expressions are derived for the mass average and molar average velocities. These velocities occur naturally in the diffusion equation if concentration is expressed as weight fraction and mole fraction of the electrolytic species, respectively.

diffusion↗
DOE OSTI2020

Protecting Spin Coherence in a Tunable Heisenberg Model

In this work, using an ensemble of atoms in an optical cavity, we engineer a family of nonlocal Heisenberg Hamiltonians with continuously tunable anisotropy of the spin-spin couplings. Therefore, we thus gain access to a rich phase diagram, including a paramagnetic-to-ferromagnetic Ising phase transition that manifests as a diverging magnetic susceptibility at the critical point. The susceptibility displays a symmetry between Ising interactions and XY (spin-exchange) interactions of the opposite sign, which is indicative of the spatially extended atomic system behaving as a single collective spin. Images of the magnetization dynamics show that spin-exchange interactions protect the coherence of the collective spin, even against inhomogeneous fields that completely dephase the noninteracting and Ising systems. Our results underscore prospects for harnessing spin-exchange interactions to enhance the robustness of spin squeezing protocols.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Simulating spin dynamics of supersolid states in a quantum Ising magnet

Motivated by a recent experimental study on the quantum Ising magnet K 2 Co(SeO 3 ) 2 that presented spectroscopic evidence of zero-field supersolidity (Chen et al., arXiv:2402.15869), we simulate the excitation spectrum of the corresponding microscopic XXZ model for the compound, using the recently developed excitation ansatz for infinite projected entangled-pair states. Here, we map out the ground state phase diagram and compute the dynamical spin structure factors across a range of magnetic field strengths, focusing especially on the two supersolid phases found near zero and saturation fields. Our simulated excitation spectra for the zero-field supersolid “Y” phase are in excellent agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Frustrated Heisenberg J 1 – J 2 model within the stretched diamond lattice of LiYbO 2

We investigate the magnetic properties of LiYbO 2 , containing a three-dimensionally frustrated, diamondlike lattice via neutron scattering, magnetization, and heat capacity measurements. The stretched diamond network of Yb 3+ ions in LiYbO 2 enters a long-range incommensurate, helical state with an ordering wave vector k = (0.384,±0.384,0) that “locks-in” to a commensurate k = (1/3,±1/3,0) phase under the application of a magnetic field. The spiral magnetic ground state of LiYbO 2 can be understood in the framework of a Heisenberg J 1 –J 2 Hamiltonian on a stretched diamond lattice, where the propagation vector of the spiral is uniquely determined by the ratio of J 2 /|J 1 |. Here, the pure Heisenberg model, however, fails to account for the relative phasing between the Yb moments on the two sites of the bipartite lattice, and this detail as well as the presence of an intermediate, partially disordered, magnetic state below 1 K suggests interactions beyond the classical Heisenberg description of this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Paramagnon heat capacity in (Ti,Zr,Hf)NiFe x NiSn half-Heusler composites

As a measure of the temperature response of the energy of matter, the heat capacity $C_p$ is a fundamental thermodynamic property. Its dependence on magnetic field, especially at low temperatures, yields insight into the electronic, phononic, and magnetic states of condensed matter. Here, we present a set of paramagnetic and ferromagnetic (Ti, Zr, Hf)NiFe x Sn half-Heusler composites that exhibit low-field (<3 T) maxima in $C_p$ and higher-field magnetic quenching of the heat capacity at temperatures below 10 K. Using rigorous statistical analysis, we attribute the effect to the existence of paramagnons within the compounds. To explain the lowest-temperature (<4 K), low-field declines in $C_p$, we derive a magnon model up to fourth order in dispersion. While the combined paramagnon and magnon model matches the data well, the fit parameters are significantly underdetermined. Further, we provide a qualitative explanation of the secondary effect based on superconducting phases within the composites. Overall, our work highlights the insight of field-dependent heat capacity studies at fixed temperatures that cannot be as easily gleaned from the temperature-dependent heat capacity at fixed magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
DOE OSTI2023

An Interactive Spectral Analysis Tool for Chemical Identification and Quantification of Gas-Phase Species in Complex Spectra

A spectral analysis tool has been developed to interactively identify and quantify individual gas-phase species from complex infrared absorbance spectra obtained from laboratory or field data. The SpecQuant program has an intuitive graphical interface that accommodates both reference and experimental data with varying resolution and instrumental lineshape, as well as algorithms to readily align the wavenumber axis of a sample spectrum with the raster of a reference spectrum. Using a classical least squares model in conjunction with reference spectra such as those from the Pacific Northwest National Laboratory (PNNL) gas-phase infrared database or simulated spectra derived from the HITRAN line-by-line database, the mixing ratio of each identified species is determined along with its associated estimation error. After correcting the wavelength and intensity of the field data, SpecQuant displays the calculated mixing ratio versus the experimental data for each analyte along with the residual spectrum with any or all analyte fits subtracted for visual inspection of the fit and residuals. The software performance for multianalyte quantification was demonstrated using moderate resolution (0.5 cm –1 ) infrared spectra that were collected during the time-resolved infrared photolysis of methyl iodide.

Thompson, Christopher J.↗
DOE OSTI2023

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Data-driven analysis and prediction of stable phases for high-entropy alloy design

High-entropy alloys (HEAs) represent a promising class of materials with exceptional structural and functional properties. However, their design and optimization pose challenges due to the large composition-phase space coupled with the complex and diverse nature of the phase formation dynamics. In this study, a data-driven approach that utilizes machine learning (ML) techniques to predict HEA phases and their composition-dependent phases is proposed. By employing a comprehensive dataset comprising 5692 experimental records encompassing 50 elements and 11 phase categories, we compare the performance of various ML models. Our analysis identifies the most influential features for accurate phase prediction. Furthermore, the class imbalance is addressed by employing data augmentation methods, raising the number of records to 1500 in each category, and ensuring a balanced representation of phase categories. The results show that XGBoost and Random Forest consistently outperform the other models, achieving 86% accuracy in predicting all phases. Additionally, this work provides an extensive analysis of HEA phase formers, showing the contributions of elements and features to the presence of specific phases. We also examine the impact of including different phases on ML model accuracy and feature significance. Notably, the findings underscore the need for ML model selection based on specific applications and desired predictions, as feature importance varies across models and phases. This study significantly advances the understanding of HEA phase formation, enabling targeted alloy design and fostering progress in the field of materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Anisotropic domains and antiferrodistortive-transition controlled magnetization in epitaxial manganite films on vicinal SrTiO 3 substrates

We studied the microstructural evolution and magnetism of ferroelastic La 0.9 Sr 0.1 MnO 3 (LSMO) epitaxial thin films grown on SrTiO 3 (001) substrates with different miscut angles. The substrate miscut angle plays a critical role in controlling the in-plane magnetic anisotropy. The microscopic origin of such magnetic anisotropy is attributed to the formation of anisotropic stripe domains along the surface step terraces. The magnetization in the LSMO films was found to be selectively modulated by the antiferrodistortive phase transition of the SrTiO 3 substrate. This phenomenon has been qualitatively explained by a strain modified Stoner–Wohlfarth model. We conclude that the magnetization modulation by the SrTiO 3 phase transition depends on h, the ratio of applied magnetic field to the saturation field. Such modulation is only visible with h < 1. The established domain microstructure–anisotropy–magnetism correlation in manganite films can be applied to a variety of complex oxide thin films on vicinal substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraints on early dark energy from the axion weak gravity conjecture

A popular proposal for resolving the Hubble tension involves an early phase of dark energy, driven by an axion field with a periodic potential. Here, in this paper, we argue that these models are tightly constrained by the axion weak gravity conjecture: for typical parameter values, the axion decay constant must satisfy f < 0.008,M Pl , which is smaller than the axion decay constants appearing in the vast majority of early dark energy models to date. We discuss possible ways to evade or loosen this constraint, arguing that its loopholes are small and difficult to thread. This suggests that it may prove challenging to realize early dark energy models in a UV complete theory of quantum gravity.

string theory and cosmology↗