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.

138 records · Page 8

Quantum Quench Dynamics–Crossover Phenomena in Nonequilibrium Correlated Quantum Systems (Final Technical Report)

In this project we have developed analytical and numerical tools to study the interplay between interactions, dimensionality, fluctuations and disorder in quantum impurity systems. We have studied the transient dynamics following a quench, i.e. the sudden change in a control variable of the model Hamiltonian, which offers a unique way to examine the system’s spatial, temporal and energetic structure. Specifically, we have investigated how crossover scales, i.e. given by Kondo couplings, are affected by disorder. We have addressed a wide range of phenomena, such as the effects of static disorder on dynamical screening of magnetic impurities, and the influence of a fully time-fluctuating bath environment on the dephasing of system qubits. We have found that many of these systems exhibit fascinating counterintuitive behavior, such as e.g. the possibility of asymptotic unitary time evolution in specifically conditioned system-bath Hamiltonians. From a technical point of view, we have developed a supersymmetric formalism to analyze scaling functions of disordered quantum impurity spin chains, we generalized the strong-coupling real-space renormalization group method to probe the quench dynamics of higher-dimensional quantum impurity systems, and we devised a master equation approach that included the effects of a fully fluctuating environment. This way, we have addressed a very broad spectrum of transient phenomena in quantum impurity systems, with disorder ranging from the quenched to the annealed limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extracting Higher Central Charge from a Single Wave Function

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_ is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_ to gapping out the edge. Here, in this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν = 1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1) 2 and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ADG: automated generation and evaluation of many-body diagrams

The goal of the present paper is twofold. First, a novel expansion many-body method applicable to superfluid open-shell nuclei, the so-called Bogoliubov in-medium similarity renormalization group (BIMSRG) theory, is formulated. This generalization of standard single-reference IMSRG theory for closed-shell systems parallels the recent extensions of coupled cluster, self-consistent Green’s function or many-body perturbation theory. Within the realm of IMSRG theories, BIMSRG provides an interesting alternative to the already existing multi-reference IMSRG (MR-IMSRG) method applicable to open-shell nuclei. The algebraic equations for low-order approximations, i.e., BIMSRG(1) and BIMSRG(2), can be derived manually without much difficulty. However, such a methodology becomes already impractical and error prone for the derivation of the BIMSRG(3) equations, which are eventually needed to reach high accuracy. Based on a diagrammatic formulation of BIMSRG theory, the second objective of the present paper is thus to describe the third version (v3.0) of the code that automatically (1) generates all valid BIMSRG(n) diagrams and (2) evaluates their algebraic expressions in a matter of seconds. This is achieved in such a way that equations can easily be retrieved for both the flow equation and the Magnus expansion formulations of BIMSRG. Expanding on this work, the first future objective is to numerically implement BIMSRG(2) (eventually BIMSRG(3)) equations and perform ab initio calculations of mid-mass open-shell nuclei.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗

Higher Hall conductivity from a single wave function: Obstructions to symmetry-preserving gapped edge of (2+1)-dimensional topological order

A (2+1)D topologically ordered phase with U(1) symmetry may or may not have a symmetric gapped edge state, even if both thermal and electric Hall conductivity are vanishing. It has recently been discovered that there are “higher” versions of Hall conductivity valid for fermionic fractional quantum Hall (FQH) states that obstruct symmetry-preserving gapped edge states beyond thermal and electric Hall conductivity. In this paper, we show that one can extract higher Hall conductivity from a single wave function of an FQH state, by evaluating the expectation value of the “partial rotation” unitary, which is a combination of partial spatial rotation and a U(1) phase rotation. This result is verified numerically with the fermionic Laughlin state with 𝜈=1/3 and 1/5, as well as the non-Abelian Moore-Read state. Together with topological entanglement entropy, we prove that the expectation values of the partial rotation completely determine if a bosonic/fermionic Abelian topological order with U(1) symmetry has a symmetry-preserving gappable edge state or not. We also show that thermal and electric Hall conductivity of Abelian topological order can be extracted by partial rotations. Even in non-Abelian FQH states, partial rotation provides the Lieb-Schultz-Mattis type theorem constraining the low-energy spectrum of the bulk-boundary system. The generalization of higher Hall conductivity to the case with Lie group symmetry is also presented.

2-dimensional systems↗

A fragmentation approach to jet flavor

An intuitive definition of the partonic flavor of a jet in quantum chromodynamics is often only well-defined in the deep ultraviolet, where the strong force becomes a free theory and a jet consists of a single parton. However, measurements are performed in the infrared, where a jet consists of numerous particles and requires an algorithmic procedure to define their phase space boundaries. To connect these two regimes, we introduce a novel and simple partonic jet flavor definition in the infrared. We define the jet flavor to be the net flavor of the partons that lie exactly along the direction of the Winner-Take-All recombination scheme axis of the jet, which is safe to all orders under emissions of soft particles, but is not collinear safe. Collinear divergences can be absorbed into a perturbative fragmentation function that describes the evolution of the jet flavor from the ultraviolet to the infrared. The evolution equations are linear and a small modification to traditional DGLAP and we solve them to leading-logarithmic accuracy. The evolution equations exhibit fixed points in the deep infrared, we demonstrate quantitative agreement with parton shower simulations, and we present various infrared and collinear safe observables that are sensitive to this flavor definition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum magnetism of iron-based ladders: Blocks, spirals, and spin flux

Motivated by increasing experimental evidence of exotic magnetism in low-dimensional iron-based materials, we present a comprehensive theoretical analysis of magnetic states of the multiorbital Hubbard ladder in the orbital-selective Mott phase (OSMP). Here, the model we used is relevant for iron-based compounds of the AFe 2 X 3 family (where A=Cs, Rb, Ba, K are alkali metals and X=S, Se are chalcogenides). To reduce computational effort, and obtain almost exact numerical results in the ladder geometry, we utilize a low-energy description of the Hubbard model in the OSMP—the generalized Kondo-Heisenberg Hamiltonian. Our main result is the doping vs interaction magnetic phase diagram. We reproduce the experimental findings on the AFe 2 X 3 materials, especially the exotic block magnetism of BaFe 2 Se 3 (antiferromagnetically coupled 2×2 ferromagnetic islands of the ↑↑↓↓ form). As in recent studies of the chain geometry, we also unveil block magnetism beyond the 2×2 pattern (with block sizes varying as a function of the electron doping) and also an interaction-induced frustrated block-spiral state (a spiral order of rigidly rotating ferromagnetic islands). Moreover, we predict new phases beyond the one-dimensional system: a robust regime of phase separation close to half filling, incommensurate antiferromagnetism for weak interaction, and a quantum spin-flux phase of staggered plaquette spin currents at intermediate doping. Finally, exploiting the bonding/antibonding band occupations, we provide an intuitive physical picture giving insight into the structure of the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gluon double-spin asymmetry in the longitudinally polarized p + p collisions

We derive the first-ever small-x expression for the inclusive gluon production cross section in the central rapidity region of the longitudinally polarized proton-proton collisions. The cross section depends on the polarizations of both protons, therefore comprising the numerator of the longitudinal double-spin asymmetry ALL for the produced gluons. The cross section is calculated in the shock wave formalism and is expressed in terms of the polarized dipole scattering amplitudes on the projectile and target protons. We show that the small-x evolution corrections are included into our cross section expression if one evolves these polarized dipole amplitudes using the double-logarithmic helicity evolution derived in [1–4]. Our calculation is performed for the gluon sector only, with the quark contribution left for future work. When that work is complete, the resulting formula will be applicable to longitudinally polarized proton-proton and proton-nucleus collisions, as well as to polarized semi-inclusive deep inelastic scattering (SIDIS) on a proton or a nucleus. Our results should allow one to extend the small-x helicity phenomenology analysis of [5] to the jet/hadron production data reported for the longitudinally polarized proton-proton collisions at RHIC and to polarized SIDIS measurements at central rapidities to be performed at the EIC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Projection-to-Born-improved subtractions at NNLO

While the current frontier in fixed-order precision for collider observables is N 3 LO, important steps are necessary to consolidate NNLO cross-section predictions with improved stability and efficiency. Slicing methods have been successfully applied to obtain NNLO and N 3 LO predictions, but have shown poor performance in the presence of fiducial cuts due to large kinematical power corrections. In this paper we implement Projection-to-Born-improved q T (P2B q T ) and jettiness (P2Bτ 0 ) subtractions for a large class of color singlet processes in MCFM. This method allows for the efficient evaluation of fiducial power corrections in any non-local subtraction scheme using a Projection-to-Born subtraction. We demonstrate the significant numerical improvements of this method based on fiducial Drell-Yan and Higgs cross-sections. Moreover, with fiducial power corrections removed via this method, the leading-logarithmic power corrections that have only been calculated without fiducial cuts can be included, further improving the calculations. For di-photon production with photon isolation, we devise a novel method in combination with P2B-improved subtractions, which we name P2B γ τ 0 , and P2B γ q T for the two subtraction schemes, respectively. This method allows the inclusion of both fiducial power corrections due to kinematic cuts on the photons and a set of isolation power corrections in the fragmentation channel where a quark may enter the isolation cone. We find significant improvements in the convergence of NNLO di-photon cross-sections with photon isolation cuts, demonstrating that it is possible to achieve a stable and efficient calculation of di-photon cross-sections using slicing methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kaon physics without new physics in $$ \varepsilon _K$$

Abstract Despite the observation of significant suppressions of $$b\rightarrow s\mu ^+\mu ^-$$ b → s μ + μ - branching ratios no clear sign of New Physics (NP) has been identified in $$\Delta F=2$$ Δ F = 2 observables $$\Delta M_{d,s}$$ Δ M d , s , $$\varepsilon _K$$ ε K and the mixing induced CP asymmetries $$S_{\psi K_S}$$ S ψ K S and $$S_{\psi \phi }$$ S ψ ϕ . Assuming negligible NP contributions to these observables allows to determine CKM parameters without being involved in the tensions between inclusive and exclusive determinations of $$|V_{cb}|$$ | V cb | and $$|V_{ub}|$$ | V ub | . Furthermore this method avoids the impact of NP on the determination of these parameters present likely in global fits. Simultaneously it provides SM predictions for numerous rare K and B branching ratios that are most accurate to date. Analyzing this scenario within $$Z^\prime $$ Z ′ models we point out, following the 2009 observations of Monika Blanke and ours of 2020, that despite the absence of NP contributions to $$\varepsilon _K$$ ε K , significant NP contributions to $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ , $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ , $$K_S\rightarrow \mu ^+\mu ^-$$ K S → μ + μ - , $$K_L\rightarrow \pi ^0\ell ^+\ell ^-$$ K L → π 0 ℓ + ℓ - , $$\varepsilon '/\varepsilon $$ ε ′ / ε and $$\Delta M_K$$ Δ M K can be present. In the simplest scenario, this is guaranteed, as far as flavour changes are concerned, by a single non-vanishing imaginary left-handed $$Z^\prime $$ Z ′ coupling $$g^L_{sd}$$ g sd L . This scenario implies very stringent correlations between the Kaon observables considered by us. In particular, the identification of NP in any of these observables implies automatically NP contributions to the remaining ones under the assumption of non-vanishing flavour conserving $$Z^\prime $$ Z ′ couplings to $$q{\bar{q}}$$ q q ¯ , $$\nu {\bar{\nu }}$$ ν ν ¯ , and $$\mu ^+\mu ^-$$ μ + μ - . A characteristic feature of this scenario is a strict correlation between $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ and $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ branching ratios on a branch parallel to the Grossman-Nir bound. Moreover, $$\Delta M_K$$ Δ M K is automatically suppressed as seems to be required by the results of the RBC-UKQCD lattice QCD collaboration. Furthermore, there is no NP contribution to $$K_L\rightarrow \mu ^+\mu ^-$$ K L → μ + μ - which otherwise would bound NP effects in $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ . Of particular interest are the correlations of $$K^+\rightarrow \pi ^+\nu {\bar{\nu }}$$ K + → π + ν ν ¯ and $$K_{L}\rightarrow \pi ^0\nu {\bar{\nu }}$$ K L → π 0 ν ν ¯ branching ratios and of $$\Delta M_K$$ Δ M K with the ratio $$\varepsilon '/\varepsilon $$ ε ′ / ε . We investigate the impact of renormalization group effects in the context of the SMEFT on this simple scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗