Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical renormalization group”

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

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Plaquette valence bond state in the spin-$\frac{1}{2}$ J 1 –J 2 XY model on the square lattice

Here, we studied the ground-state phase diagram of spin-1/2 J 1 –J 2 XY model on the square lattice with first-neighbor J 1 and second-neighbor J 2 antiferromagnetic interactions using both infinite density matrix renormalization group (iDMRG) and DMRG approaches. We show that a plaquette valence bond phase is realized in an intermediate region 0.50 ≤ J 2 /J 1 < 0.54 between a Néel magnetic ordered phase at J 2 /J 1 < 0.50 and a stripy magnetic ordered phase at J 2 /J 1 ≥ 0.54 . The plaquette valence bond phase is characterized by finite dimer orders in both the horizontal and vertical directions. Contrary to the spin-1/2 J 1 –J 2 Heisenberg model, we do not find numerical evidence for a quantum spin liquid phase in the J 1 –J 2 XY model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Superconductivity in the doped quantum spin liquid on the triangular lattice

Abstract Broad interest in quantum spin liquid (QSL) phases was triggered by the notion that they can be viewed as insulating phases with preexisting electron pairs, such that upon light doping they might automatically yield high temperature superconductivity. Yet despite intense experimental and numerical efforts, definitive evidence showing that doping QSLs leads to superconductivity has been lacking. We address the problem of a lightly doped QSL through a large-scale density-matrix renormalization group study of the t - J model on finite-circumference triangular cylinders with a small but nonzero concentration of doped holes. We provide direct evidences that doping QSL can naturally give rise to d -wave superconductivity. Specifically, we find power-law superconducting correlations with a Luttinger exponent, K s c ≈ 1, which is consistent with a strongly diverging superconducting susceptibility, $${\chi }_{sc} \,\sim\, {T}^{-(2\,-\,{K}_{sc})}$$ χ s c ~ T − ( 2 − K s c ) as the temperature T → 0. The spin–spin correlations—as in the undoped QSL state—fall exponentially which suggests that the superconducting pair-pair correlations evolve smoothly from the insulating parent state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Emergent pair density wave order across a Lifshitz transition

We numerically investigate the telltale signs of pair-density-wave order (PDW) in the Kondo-Heisenberg chain by focusing on the momentum resolved spectrum in different parameter regimes. Density matrix renormalization group calculations reveal that this phase is characterized by a dispersion with two minima and four Fermi points, indicating the emergence of an effective next-nearest-neighbor hopping that arises as a second-order effect to avoid magnetic frustration. The pairs appear in the spectrum as in-gap bound states with weight concentrated in the hole pockets. The low-energy physics can be understood by means of a generalized t−J model with next-nearest-neighbor hopping. Our results offer a guide for searching for experimental signatures, and for other models that can realize PDW physics.

Yang, Luhang [University of Tennessee Knoxville]↗
DOE OSTI2023

Global quantum phase diagram and non-Abelian chiral spin liquid in a spin- 3 2 square-lattice antiferromagnet

Since strong quantum fluctuations are essential for the emergence of quantum spin liquids, there have been extensive exploration and identification of spin liquid candidates in spin-$\frac{1}{2}$ systems, while such activities are rare in higher spin systems. Here we report an example of non-Abelian chiral spin liquid emerging in a spin-$\frac{3}{2}$ Heisenberg model on a square lattice. By tuning Heisenberg exchange interaction and scalar chirality interaction, we map out a quantum phase diagram enclosing three conventional magnetic orders and a chiral spin liquid based on density-matrix renormalization group studies. The nature of the spin liquid is identified as a long-sought bosonic version of the Read-Rezayi state that supports non-Abelian Fibonacci anyonic statistics, identified by the ground state entanglement spectrum. Significantly, we establish that the non-Abelian chiral spin liquid emerges through the enlarged local degrees of freedom and enhanced quantum fluctuations near the classical phase boundaries of competing magnetic orders. Finally, our numerical discovery of an exotic quantum spin liquid in a spin-$\frac{3}{2}$ system suggests a route for discovering fractionalized quantum phases in frustrated higher spin magnetic compounds.

2-dimensional systems↗
DOE OSTI2023

Enhanced superconductivity by near-neighbor attraction in the doped extended Hubbard model

A recent experiment has unveiled an anomalously strong electron-electron attraction in the one-dimensional copper-oxide chain Ba 2–x Sr x CuO 3+δ . While the effect of the near-neighbor electron attraction V in the one-dimensional extended Hubbard chain has been examined recently, its effect in the Hubbard model beyond the one-dimensional chain remains unclear. Here, in this work, we report a density matrix renormalization group study of the extended Hubbard model on long four-leg cylinders on the square lattice. We find that the near-neighbor electron attraction V can notably enhance the long-distance superconducting correlations while simultaneously suppressing the charge-density-wave (CDW) correlations. Specifically, for a modestly strong electron attraction, the superconducting correlations become dominant over the CDW correlations with a Luttinger exponent K sc ~1 and strong divergent superconducting susceptibility. Our results provide a promising way to realize long-range superconductivity in the doped Hubbard model. The relevance of our numerical results to cuprate materials is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Enhancing detection of topological order by local error correction

The exploration of topologically-ordered states of matter is a long-standing goal at the interface of several subfields of the physical sciences. Such states feature intriguing physical properties such as long-range entanglement, emergent gauge fields and non-local correlations, and can aid in realization of scalable fault-tolerant quantum computation. However, these same features also make creation, detection, and characterization of topologically-ordered states particularly challenging. Motivated by recent experimental demonstrations, we introduce a paradigm for quantifying topological states—locally error-corrected decoration (LED)—by combining methods of error correction with ideas of renormalization-group flow. Our approach allows for efficient and robust identification of topological order, and is applicable in the presence of incoherent noise sources, making it particularly suitable for realistic experiments. We demonstrate the power of LED using numerical simulations of the toric code under a variety of perturbations. We subsequently apply it to an experimental realization, providing new insights into a quantum spin liquid created on a Rydberg-atom simulator. Finally, we extend LED to generic topological phases, including those with non-abelian order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

SMEFT ATLAS: The landscape beyond the Standard Model

The Standard Model Effective Field Theory (SMEFT) based on the unbroken gauge group SU(3) 𝐶 ⊗SU(2) 𝐿 ⊗U(1) 𝑌 and containing only particles of the Standard Model (SM) has developed in the last decade to a mature field. It is the framework to be used in the energy gap from scales sufficiently higher than the electroweak scale up to the lowest energy scale at which new particles show up. We summarize the present status of this theory with a particular emphasize on its role in the indirect search for new physics (NP). While flavour physics of both quarks and leptons is the main topic of our review, we also discuss electric dipole moments, anomalous magnetic moments (𝑔−2) 𝜇,𝑒 , 𝑍-pole observables, Higgs observables and high-𝑝𝑇 scattering processes within the SMEFT. We group the observables into ten classes and list for each class the most relevant operators and the corresponding renormalization group equations (RGEs). We exhibit the correlations between different classes implied both by the operator mixing and the SU(2) 𝐿 gauge symmetry. Our main goal is to provide an insight into the complicated operator structure of this framework which hopefully will facilitate the identification of valid ultraviolet completions behind possible anomalies observed in future data. Numerous colourful charts, and 85 tables, while representing rather complicated RG evolution from the NP scale down to the electroweak scale, beautify the involved SMEFT landscape. Over 1000 references to the literature underline the importance and the popularity of this field. We discuss both top-down and bottom-up approaches as well as their interplay. This allows us eventually to present an atlas of different landscapes beyond the SM that includes heavy gauge bosons and scalars, vector-like quarks and leptons and leptoquarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2023

Auger spectroscopy beyond the ultra-short core-hole relaxation time approximation

Abstract We present a time-dependent computational approach to study Auger electron spectroscopy (AES) beyond the ultra-short core-hole relaxation time approximation and, as a test case, we apply it to the paradigmatic example of a one-dimensional Mott insulator represented by a half-filled Hubbard chain. The Auger spectrum is usually calculated by assuming that, after the creation of a core-hole, the system thermalizes almost instantaneously. This leads to a relatively simple analytical expression that uses the ground-state with a core-hole as a reference state and ignores all the transient dynamics related to the screening of the core-hole. In this picture, the response of the system can be associated to the pair spectral function. On the other hand, in our numerical calculations, the core hole is created by a light pulse, allowing one to study the transient dynamics of the system in terms of the pulse duration and in the non-perturbative regime. Time-dependent density matrix renormalization group calculations reveal that the relaxation process involves the creation of a polarization cloud of doublon excitations that have an effect similar to photo-doping. As a consequence, there is a leak of spectral weight to higher energies into what otherwise would be the Mott gap. For longer pulses, these excited states, mostly comprised of doublons, can dominate the spectrum. By changing the duration of the light-pulse, the entire screening process can be resolved in time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Gravitational Fluctuations as an Alternative to Inflation III. Numerical Results

Power spectra play an important role in the theory of inflation, and their ability to reproduce current observational data to high accuracy is often considered a triumph of inflation, largely because of a lack of credible alternatives. In previous work we introduced an alternative picture for the cosmological power spectra based on the nonperturbative features of the quantum version of Einstein’s gravity, instead of currently popular inflation models based on scalar fields. The key ingredients in this new picture are the appearance of a nontrivial gravitational vacuum condensate (directly related to the observed cosmological constant), and a calculable renormalization group running of Newton’s G on cosmological scales. More importantly, one notes the absence of any fundamental scalar fields in this approach. Results obtained previously were largely based on a semi-analytical treatment, and thus, while generally transparent in their implementation, often suffered from the limitations of various approximations and simplifying assumptions. In this work, we extend and refine our previous calculations by laying out an updated and extended analysis, which now utilizes a set of suitably modified state-of-the-art numerical programs (ISiTGR, MGCAMB and MGCLASS) developed for observational cosmology. As a result, we are able to remove some of the approximations employed in our previous studies, leading to a number of novel and detailed physical predictions. These should help in potentially distinguishing the vacuum condensate picture of quantum gravity from that of other models such as scalar field inflation. Here, besides the matter power spectrum P m (k) , we work out, in detail, predictions for what are referred to as the TT, TE, EE, BB angular spectra, as well as their closely related lensing spectra. However, the current limited precision of observational data today (especially on large angular scales) does not allow us yet to clearly prove or disprove either set of ideas. Nevertheless, by exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with an expected future influx of increasingly accurate observational data, one can hope that the new quantum gravitational picture can be subjected to further stringent tests in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient simulation of moiré materials using the density matrix renormalization group

Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Triplet pair density wave superconductivity on the π -flux square lattice

Pair-density waves (PDWs) are superconducting states that spontaneously break translation symmetry in systems with time-reversal symmetry (TRS). Evidence for PDWs has been seen in several recent experiments, as well as in the pseudogap regime in cuprates. Theoretical understanding of PDWs has been largely restricted to phenomenological and numerical studies, while microscopic theories typically require strong coupling or fine tuning. In this work, we provide a symmetry-based mechanism under which PDWs emerge as a weak-coupling instability of a two-dimensional TRS metal. Combining mean-field and renormalization-group analyses, we identify a weak-coupling instability towards a triplet PDW realized in the π-flux square lattice model with on-site repulsion and moderate nearest-neighbor attraction when the Fermi level crosses Van Hove singularities at 1/4 and 3/4 fillings. Further, this PDW is protected by the magnetic translation symmetries characteristic of Hofstadter systems, of which the π-flux lattice is a special time-reversal-symmetric case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths U and for degenerate orbitals. At half filling and large U, the ferromagnetic Hund's coupling, J H , generates effective spin-1 moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we accurately study the system's behavior under light hole-doping. For U = 1.6 in units of the noninteracting bandwidth and for J H /U ≳ 0.275, we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity.

1-dimensional systems↗
DOE OSTI2024

A numerical study of bounds in the correlations of fractional quantum Hall states

We numerically compute the guiding center static structure factor \bar S(k) S ‾ ( k ) of various fractional quantum Hall (FQH) states to \mathcal{O}(k\ell)^6 𝒪 ( k ℓ ) 6 where k k is the wavenumber and \ell ℓ is the magnetic length. Employing density matrix renormalization group on an infinite cylinder of circumference L_y L y , we study the two-dimensional limit using L_y/\xi \gg 1 L y / ξ ≫ 1 , where \xi ξ is the correlation length. The main findings of our work are: 1) the ground states that deviate away from the ideal conformal block wavefunctions, do not saturate the Haldane bound, and 2) the coefficient of O(k\ell)^6 O ( k ℓ ) 6 term appears to be bounded above by a value predicted by field theories proposed in the literature. The first finding implies that the graviton mode is not maximally chiral for experimentally relevant FQH states.

Kumar, Prashant (ORCID:0000000258002768)↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗