Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

From trivial to topological paramagnets: The case of Z 2 and Z 2 3 symmetries in two dimensions

Using quantum Monte Carlo simulations, we map out the phase diagram of Hamiltonians interpolating between trivial and nontrivial bosonic symmetry-protected topological phases, protected by Z 2 and Z 2 3 symmetries, in two dimensions. In all cases, we find that the trivial and the topological phases are separated by an intermediate phase in which the protecting symmetry is spontaneously broken. Depending on the model, we identify a variety of magnetic orders on the triangular lattice, including ferromagnetism, 3 × 3 order, and stripe orders (both commensurate and incommensurate). Critical properties are determined through a finite-size scaling analysis. Possible scenarios regarding the nature of the phase transitions are discussed.

2-dimensional systems↗
DOE OSTI2021

Evidence for deconfined U(1) gauge theory at the transition between toric code and double semion

Building on quantum Monte Carlo simulations, we study the phase diagram of a one-parameter Hamiltonian interpolating between trivial and topological Ising paramagnets in two dimensions, which are dual to the toric code and the double semion. We discover an intermediate phase with stripe order which spontaneously breaks the protecting Ising symmetry. Remarkably, we find evidence that this intervening phase is gapless due to the incommensurability of the stripe pattern and that it is dual to a $U(1)$ gauge theory exhibiting Cantor deconfinement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Composite anyons on a torus

An adiabatic approach put forward by Greiter and Wilczek interpolates between the integer quantum Hall effects of electrons and composite fermions by varying the statistical flux bound to electrons continuously from zero to an even integer number of flux quanta, such that the intermediate states represent anyons in an external magnetic field with the same “effective” integer filling factor. We consider such anyons on a torus, and construct representative wave functions for their ground as well as excited states. These wave functions involve higher Landau levels in general, but can be explicitly projected into the lowest Landau level for many parameters. We calculate the variational energy gap between the first excited state and ground state and find that it remains open as the statistical phase is varied. Lastly, we obtain from these wave functions, both analytically and numerically, various topological quantities, such as ground-state degeneracy, the Chern number, and the Hall viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Strain-tunable metamagnetic critical endpoint in Mott insulating rare-earth titanates

Rare-earth titanates are Mott insulators whose magnetic ground state—antiferromagnetic (AFM) or ferromagnetic (FM)—can be tuned by the radius of the rare-earth element. We combine phenomenology and first-principles calculations to shed light on the generic magnetic phase diagram of a chemically substituted titanate on the rare-earth site that interpolates between an AFM and a FM state. Octahedral rotations present in these perovskites cause the AFM order to acquire a small FM component—and vice-versa—removing any multicritical point from the phase diagram. However, for a wide parameter range, a first-order metamagnetic transition line terminating at a critical endpoint survives inside the magnetically ordered phase. Like the liquid-gas transition, a Widom line emerges from the endpoint, characterized by enhanced fluctuations. In contrast to metallic FMs, this metamagnetic transition involves two symmetry-equivalent and insulating canted spin states. Moreover, instead of a magnetic field, we show that uniaxial strain can be used to tune this transition to zero temperature, inducing a quantum critical endpoint.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3

In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of localization on photoluminescence and zero-field splitting of silicon color centers

The study of defect centers in silicon has been recently reinvigorated by their potential applications in optical quantum information processing. A number of silicon defect centers emit single photons in the telecommunication O-band, making them promising building blocks for quantum networks between computing nodes. The two-carbon G-center, self-interstitial W-center, and spin-1/2 T-center are the most intensively studied silicon defect centers, yet despite this, there is no consensus on the precise configurations of defect atoms in these centers, and their electronic structures remain ambiguous. Here, in this work, we employ ab initio density functional theory to characterize these defect centers, providing insight into the relaxed structures, band structures, and photoluminescence spectra, which are compared to experimental results. Motivation is provided for how these properties are intimately related to the localization of electronic states in the defect centers. In particular, we present the calculation of the zero-field splitting for the excited triplet state of the G-center defect as the structure is linearly interpolated from the A-configuration to the B-configuration, showing a sudden increase in the magnitude of the D zz component of the zero-field-splitting tensor. By performing projections onto the local orbital states of the defect, we analyze this transition in terms of the symmetry and bonding character of the G-center defect, which sheds light on its potential application as a spin-photon interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Ab initio study of electron-phonon interaction in phosphorene

The monolayer of black phosphorus, or “phosphorene,” has recently emerged as a two-dimensional semiconductor with intriguing highly anisotropic transport properties. Existing calculations of its intrinsic phonon-limited electronic transport properties so far rely on the deformation potential approximation, which is in general not directly applicable to anisotropic materials since the deformation along one specific direction can scatter electrons traveling in all directions. We perform a first-principles calculation of the electron-phonon interaction in phosphorene based on density functional perturbation theory and Wannier interpolation. Our calculation reveals that (1) the high anisotropy provides extra phase space for electron-phonon scattering, and (2) optical phonons have appreciable contributions. Both effects cannot be captured by the deformation potential calculations. Furthermore, our simulation predicts carrier mobilities ~170 cm 2 /Vs for both electrons and holes at 300 K, and a thermoelectric figure of merit $zT$ of up to 0.14 in $p$-type impurity-free phosphorene at 500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Quenched crystal-field disorder and magnetic liquid ground states in Tb 2 Sn 2 - x Ti x O 7 [Crystal field disorder in the quantum spin ice ground state of Tb2Sn2-xTixO7]

Solid-solutions of the "soft" quantum spin ice pyrochlore magnets Tb 2 B 2 O 7 with B=Ti and Sn display a novel magnetic ground state in the presence of strong B-site disorder, characterized by a low susceptibility and strong spin fluctuations to temperatures below 0.1 K. These materials have been studied using ac-susceptibility and muSR techniques to very low temperatures, and time-of-flight inelastic neutron scattering techniques to 1.5 K. Remarkably, neutron spectroscopy of the Tb 3+ crystal field levels appropriate to at high B-site mixing (0.5 < x < 1.5 in Tb 2 Sn 2-x Ti x O 7 ) reveal that the doublet ground and first excited states present as continua in energy, while transitions to singlet excited states at higher energies simply interpolate between those of the end members of the solid solution. The resulting ground state suggests an extreme version of a random-anisotropy magnet, with many local moments and anisotropies, depending on the precise local configuration of the six B sites neighboring each magnetic Tb 3+ ion.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Spin excitations used to probe the nature of exchange coupling in the magnetically ordered ground state of Pr 0.5 Ca 0.5 MnO 3

We have used time-of-flight inelastic neutron scattering to measure the spin wave spectrum of the canonical half-doped manganite Pr 0.5 Ca 0.5 MnO 3 in its magnetic and orbitally ordered phase. Comparison of the data, which cover multiple Brillouin zones and the entire energy range of the excitations, with several different models shows that only the CE-type ordered state provides an adequate description of the magnetic ground state, provided interactions beyond nearest neighbor are included. We are able to rule out a ground state in which there exist pairs of dimerized spins which interact only with their nearest neighbors. The Zener polaron ground state, which comprises strongly bound magnetic dimers, can be ruled out on the basis of gross features of the observed spin wave spectrum. A model with weaker dimerization reproduces the observed dispersion but can be ruled out on the basis of subtle discrepancies between the calculated and observed structure factors at certain positions in reciprocal space. Adding further neighbor interactions results in almost no dimerization, i.e. interpolating back to the CE model. These results are consistent with theoretical analysis of the degenerate double exchange model for half-doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

First-principles mode-by-mode analysis for electron-phonon scattering channels and mean free path spectra in GaAs

We present a first-principles framework to investigate the electron scattering channels and transport properties for polar materials by combining the exact solution of the linearized electron-phonon (e-ph) Boltzmann transport equation in its integral-differential form associated with the e-ph coupling matrices obtained from the polar Wannier interpolation scheme. No ad hoc parameter is required throughout this calculation, and GaAs, a well-studied polar material, is used as an example to demonstrate this method. In this work, the long-range and short-range contributions as well as the intravalley and intervalley transitions in the e-ph interactions (EPIs) have been quantitatively addressed. Promoted by such mode-by-mode analysis, we find that in GaAs, the piezoelectric scattering is comparable to deformation-potential scattering for electron scatterings by acoustic phonons in EPI even at room temperature, and it makes a significant contribution to mobility. Furthermore, we achieved good agreement with experimental data for the mobility, and we identified that electrons with mean free paths between 130 and 210 nm provide the dominant contribution to the electron transport at 300 K. As a result, such information provides a deeper understanding of the electron transport in GaAs, and the presented framework can be readily applied to other polar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Spin-orbit coupled systems in the atomic limit: rhenates, osmates, iridates

Motivated by RIXS experiments on a wide range of complex heavy oxides, including rhenates, osmates, and iridates, we discuss the theory of RIXS for site-localized t 2 g orbital systems with strong spin-orbit coupling. For such systems, we present exact diagonalization results for the spectrum at different electron fillings, showing that it accesses “single-particle” and “multiparticle” excitations. This leads to a simple picture for the energies and intensities of the RIXS spectra in Mott insulators such as double perovskites which feature highly localized electrons, and yields estimates of the spin-orbit coupling and Hund's coupling in correlated 5d oxides. We present new higher resolution RIXS data at the Re L 3 edge in Ba 2 YReO 6 which finds a previously unresolved peak splitting, providing further confirmation of our theoretical predictions. Using ab initio electronic structure calculations on Ba 2 MReO 6 (with M=Re, Os, Ir) we show that while the atomic limit yields a reasonable effective Hamiltonian description of the experimental observations, effects such as t 2 g –e g interactions and hybridization with oxygen are important. Our ab initio estimate for the strength of the intersite exchange coupling shows that, compared to the d 3 systems, the exchange is one or two orders of magnitude weaker in the d 2 and d 4 materials, which may partly explain the suppression of long-range magnetic order in the latter compounds. As a way to interpolate between the site-localized picture and our electronic structure band calculations, we discuss the spin-orbital levels of the MO 6 cluster. In conclusion, this suggests a possible role for intracluster excitons in Ba 2 YIrO 6 which may lead to a weak breakdown of the atomic J eff = 0 picture and to small magnetic moments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Mixing A -type and G -type B -site antiferromagnetism in A M n 1 - x F e x O 3 ( A = La , Nd )

Using a classical Heisenberg model with anisotropic nearest-neighbor superexchange and uniaxial anisotropy, I investigate the magnetic order of the transition metals on the B sites in AMn 1-x Fe x O 3 (A=La, Nd; x=[0,1]) with Monte Carlo methods. Magnetic parameters are extracted from preexisting experimental data when available and otherwise estimated with calculations. Superexchange energies are compared to density functional theory and anisotropy energies to ligand field theory, and the potential for doping-induced strain affecting these parameters is considered. I find coexistence regions of different types of magnetic order (A-type and G-type antiferromagnetism), regions with distinct magnetic transitions for the two order parameters, noncollinear ground states for intermediate values of x that decrease the amplitude of the ordered moment in a way consistent with experiment for low dopings, and a reproduction of the experimentally determined, nonmonotonic interpolation of the ordering temperatures. The detailed shape of the phase boundaries in (x,T) space is found to be diagnostic of the magnetic interactions. Finally, these results also suggest a ferrimagnetic ground state for the A 2 MnFeO 6 double perovskite with T N ≈270K for A=La and Nd.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

One-dimensional Hubbard-Holstein model with finite-range electron-phonon coupling

The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling g to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction U drives the fermions into a Mott antiferromagnet. At low g and U, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range ξ, interpolating between the Holstein and Fröhlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range ξ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction U will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. Furthermore, the phase separation region is once again present in the U ≠ 0 case, even for small values of the coupling range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Applying Bayesian parameter estimation to relativistic heavy-ion collisions: Simultaneous characterization of the initial state and quark-gluon plasma medium

We quantitatively estimate properties of the quark-gluon plasma created in ultrarelativistic heavy-ion collisions utilizing Bayesian statistics and a multiparameter model-to-data comparison. The study is performed using a recently developed parametric initial condition model, T R ENTo, which interpolates among a general class of particle production schemes, and a modern hybrid model which couples viscous hydrodynamics to a hadronic cascade.We calibrate the model to multiplicity, transverse momentum, and flow data and report constraints on the parametrized initial conditions and the temperature-dependent transport coefficients of the quark-gluon plasma. We show that initial entropy deposition is consistent with a saturation-based picture, extract a relation between the minimum value and slope of the temperature-dependent specific shear viscosity, and find a clear signal for a nonzero bulk viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗