Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “derived lattice”

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.

547 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Boson localization and universality in YBa2Cu(3-x)M(x)O(7-delta)

We consider a two component mixture of charged fermions on neutralizing background with all sign combinations and arbitrarily small mass ratios. In the two impurity limit for the heavier component we show that the pair forms a bound state for all charge combinations. In the lowest order approximation we derive a closed form expression Veff(r) for the binding potential which has short-range repulsion followed by attraction. In the classical limit, when the mass of embedded particles is large m2 much greater than m, we can calculate from Veff(r) also the cohesive energy E and the bond length R of a metallic crystal such as lithium. The lowest order result is R = 3.1 A, E = -0.9 eV, not entirely different from the experimental result for lithium metal. The same interaction for two holes on a parabolic band with m2 greater than m gives the quantum mechanical bound state which one may interpret as a boson or local pair in the case of high-Te and heavy fermion superconductors. We also show that for compounds of the type YBa2Cu(3 - x)M(x)O(7 - delta) one can understand most of the experimental results for the superconducting and normal states with a single temperature dependent boson breaking function f(T) for each impurity content x governing the decay of bosons into pairing fermions. In the normal state f(T) turns out to be a linear, universal function, independent of the impurity content I and the oxygen content delta. We predict with universality a depression in Tc(x) with slight down bending in agreement with experiment. As a natural consequence of the model the bosons become localized slightly above Tc due to the Wigner crystallization, enhanced with lattice local field minima. The holes remain delocalized with a linearly increasing concentration in the normal state, thus explaining the rising Hall density. The boson localization temperature T(sub BL) shows up as a minimum in the Hall density R(sub ab)(exp -1). We also give explanation for very recently observed scaling of temperature dependent Hall effect in La(2 - x)Sr(x)CuO4.

Kallio, A.↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. As a result, it also paves the way for many-body computations of the phonon spectrum in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interplay between magnetism and band topology in the kagome magnets R Mn 6 Sn 6

Kagome-lattice magnets R Mn 6 Sn 6 recently emerged as a new platform to exploit the interplay between magnetism and topological electronic states. Some of the most exciting features of this family are the dramatic dependence of the easy magnetization direction on the rare-earth specie, despite other magnetic and electronic properties being essentially unchanged, and the kagome geometry of the Mn planes that in principle can generate flat bands and Dirac points; gapping of the Dirac points by spin-orbit coupling has been suggested recently to be responsible for the observed anomalous Hall response in the member TbMn 6 Sn 6 . In this paper, we address both issues with ab initio calculations. We have discovered the significant role played by higher-order crystal-field parameters and rare-earth magnetic anisotropy constants in these systems. Further, we demonstrate that the microscopic origin of rare-earth magnetic anisotropy can also be quantified and understood at various levels: ab initio, phenomenological, and analytical. In particular, using a simple and physically transparent analytical model based on perturbation theory, we are able to explain, with full quantitative agreement, the evolution of rare-earth magnetic anisotropy across the series. We analyze in detail the topological properties of Mn-dominated bands and demonstrate how they emerge from the multiorbital planar kagome model. We further show that, despite this fact, most of the topological features at the Brillouin zone corner K are strongly 3D and therefore cannot explain the observed quasi-2D anomalous Hall effect, while the most pronounced quasi-2D dispersion are too far removed from the Fermi level. By employing self-consistent calculations with ab initio many-body approaches, we demonstrate that the exchange-correlation effects beyond the density functional theory for itinerant Mn- d electrons do not significantly alter the obtained electronic and magnetic structure. Therefore, we conclude that, contrary to previous claims, the most pronounced 2D kagome-derived topological band features bear little relevance to transport in RMn 6 Sn 6 , albeit they may possibly be brought to focus by electron or hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Theory of the charge density wave in A V 3 Sb 5 kagome metals

The family of metallic kagome compounds AV 3 Sb 5 (A = K, Rb, Cs) was recently discovered to exhibit both superconductivity and charge order. The nature of the charge density wave (CDW) phase is presently unsettled, which complicates the interpretation of the superconducting ground state. In this paper, we use group theory and density functional theory (DFT) to derive and solve a phenomenological Landau model for this CDW state. The DFT results reveal three unstable phonon modes with the same in-plane momentum but different out-of-plane momenta, whose frequencies depend strongly on the electronic temperature. This is indicative of an electronically driven CDW, stabilized by features of the in-plane electronic dispersion. Motivated by the DFT analysis, we construct a Landau free-energy expansion for coupled CDW order parameters with wave vectors at the M and L points of the hexagonal Brillouin zone. We find an unusual trilinear term coupling these different order parameters, which can promote the simultaneous condensation of both CDWs even if the two modes are not nearly degenerate. We classify the different types of coupled multi-Q CDW orders, focusing on those that break the sixfold rotational symmetry and lead to a unit-cell doubling along all three crystallographic directions, as suggested by experiments. Here, we determine a region in parameter space, characterized by large nonlinear Landau coefficients, where these phases—dubbed staggered trihexagonal and staggered Star of David—are the leading instabilities of the system. Finally, we discuss the implications of our results for the kagome metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2016

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗