SEARCH · Search NASA

Results for “Lattice Integrable Models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2017

An exactly solvable quench protocol for integrable spin models

Quantum quenches in continuum field theory across critical points are known to display different scaling behaviours in different regimes of the quench rate. We extend these results to integrable lattice models such as the transverse field Ising model on a one-dimensional chain and the Kitaev model on a two-dimensional honeycomb lattice using a nonlinear quench protocol which allows for exact analytical solutions of the dynamics. Our quench protocol starts with a finite mass gap at early times and crosses a critical point or a critical region, and we study the behaviour of one point functions of the quenched operator at the critical point or in the critical region as a function of the quench rate. For quench rates slow compared to the initial mass gap, we find the expected Kibble-Zurek scaling. In contrast, for rates fast compared to the mass gap, but slow compared to the inverse lattice spacing, we find scaling behaviour similar to smooth fast continuum quenches. For quench rates of the same order of the lattice scale, the one point function saturates as a function of the rate, approaching the results of an abrupt quench. The presence of an extended critical surface in the Kitaev model leads to a variety of scaling exponents depending on the starting point and on the time where the operator is measured. We discuss the role of the amplitude of the quench in determining the extent of the slow (Kibble-Zurek) and fast quench regimes, and the onset of the saturation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Hydrodynamics of operator spreading and quasiparticle diffusion in interacting integrable systems

In this work, we address the hydrodynamics of operator spreading in interacting integrable lattice models. In these models, operators spread through the ballistic propagation of quasiparticles, with an operator front whose velocity is locally set by the fastest quasiparticle velocity. In interacting integrable systems, this velocity depends on the density of the other quasiparticles, so equilibrium density fluctuations cause the front to follow a biased random walk, and therefore to broaden diffusively. Ballistic front propagation and diffusive front broadening are also generically present in nonintegrable systems in one dimension; thus, although the mechanisms for operator spreading are distinct in the two cases, these coarse-grained measures of the operator front do not distinguish between the two cases. We present an expression for the front-broadening rate; we explicitly derive this for a particular integrable model (the “Floquet-Fredrickson-Andersen” model), and argue on kinetic grounds that it should apply generally. In conclusion, our results elucidate the microscopic mechanism for diffusive corrections to ballistic transport in interacting integrable models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Enhanced negative energy with a massless Dirac field

Motivated by traversable wormhole constructions that require large amounts of negative energy, we explore constraints on the amount of negative energy that can be carried by a free Dirac field in a slab-shaped region between two parallel spatial planes. Specifically, we ask what is the minimum possible uniform energy density that can exist at some time, considering all possible states and all possibilities for the physics outside the slab. The vacuum state where we identify the two sides of the slab with antiperiodic boundary conditions gives one possible state with uniform negative energy, but we argue that states with more negative energy exist above 1+1 dimensions. Technically, we reduce the problem to studying a massive Dirac field on an interval in 1+1 dimensions and numerically search for states with uniform energy density in a lattice regulated model. We succeed in finding states with enhanced negative energy (relative to the antiperiodic vacuum) which also appear to have a sensible continuum limit. Our results for the mass-dependence of the minimum uniform energy density in 1+1 dimensions suggest that for a 3+1 dimensional massless Dirac fermion, it is possible to have states with arbitrarily large uniform negative energy density in an arbitrarily wide slab.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Lattice realizations of topological defects in the critical (1+1)-d three-state Potts model

Topological/perfectly-transmissive defects play a fundamental role in the analysis of the symmetries of two dimensional conformal field theories (CFTs). In the present work, spin chain regularizations for these defects are proposed and analyzed in the case of the three-state Potts CFT. In particular, lattice versions for all the primitive defects are presented, with the remaining defects obtained from the fusion of the primitive ones. The defects are obtained by introducing modified interactions around two given sites of an otherwise homogeneous spin chain with periodic boundary condition. The various primitive defects are topological on the lattice except for one, which is topological only in the scaling limit. The lattice models are analyzed using a combination of exact diagonalization and density matrix renormalization group techniques. Low-lying energy spectra for different defect Hamiltonians as well as entanglement entropy of blocks located symmetrically around the defects are computed. The latter provides a convenient way to compute the g-function which characterizes various defects. Finally, the eigenvalues of the line operators in the “crossed channel” and fusion of different defect lines are also analyzed. The results are all in agreement with expectations from conformal field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Axially symmetric e-lens based on McMillan map

In this work, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potential of this device. We classify possible regimes with stable trajectories and provide set of canonical action-angle variables, along with an evaluation of the dynamical aperture, Poincar\'e rotation numbers as functions of amplitudes, and spread in nonlinear tunes.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Axially Symmetric E-Lens Based on McMillan Map

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Recovery Act: An Integrated Experimental and Numerical Study: Developing a Reaction Transport Model that Couples Chemical Reactions of Mineral Dissolution/Precipitation with Spatial and Temporal Flow Variations.

A total of 12 publications and 23 abstracts were produced as a result of this study. In particular, the compilation of a thermodynamic database utilizing consistent, current thermodynamic data is a major step toward accurately modeling multi-phase fluid interactions with solids. Existing databases designed for aqueous fluids did not mesh well with existing solid phase databases. Addition of a second liquid phase (CO2) magnifies the inconsistencies between aqueous and solid thermodynamic databases. Overall, the combination of high temperature and pressure lab studies (task 1), using a purpose built apparatus, and solid characterization (task 2), using XRCT and more developed technologies, allowed observation of dissolution and precipitation processes under CO2 reservoir conditions. These observations were combined with results from PIV experiments on multi-phase fluids (task 3) in typical flow path geometries. The results of the tasks 1, 2, and 3 were compiled and integrated into numerical models utilizing Lattice-Boltzmann simulations (task 4) to realistically model the physical processes and were ultimately folded into TOUGH2 code for reservoir scale modeling (task 5). Compilation of the thermodynamic database assisted comparisons to PIV experiments (Task 3) and greatly improved Lattice Boltzmann (Task 4) and TOUGH2 simulations (Task 5). PIV (Task 3) and experimental apparatus (Task 1) have identified problem areas in TOUGHREACT code. Additional lab experiments and coding work has been integrated into an improved numerical modeling code.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Reconstruction and uncertainty quantification of lattice Hamiltonian model parameters from observations of microscopic degrees of freedom

The emergence of scanning probe and electron beam imaging techniques has allowed quantitative studies of atomic structure and minute details of electronic and vibrational structure on the level of individual atomic units. These microscopic descriptors, in turn, can be associated with local symmetry breaking phenomena, representing the stochastic manifestation of the underpinning generative physical model. In this work, we explore the reconstruction of exchange integrals in the Hamiltonian for a lattice model with two competing interactions from observations of microscopic degrees of freedom and establish the uncertainties and reliability of such analysis in a broad parameter-temperature space. In contrast to other approaches, we specifically specify a loss function inherent to thermodynamic systems and utilize it to estimate uncertainty in simulated realizations of different models. As an ancillary task, we develop a machine learning approach based on histogram clustering to predict phase diagrams efficiently using a reduced descriptor space. We further demonstrate that reconstruction is possible well above the phase transition and in the regions of parameter space when the macroscopic ground state of the system is poorly defined due to frustrated interactions. This suggests that this approach can be applied to the traditionally complex problems of condensed matter physics such as ferroelectric relaxors and morphotropic phase boundary systems, spin and cluster glasses, and quantum systems once the local descriptors linked to the relevant physical behaviors are known.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Phase controlled metal–insulator transition in multi-leg quasiperiodic optical lattices

A tight-binding model of a multi-leg ladder network with a continuous quasiperiodic modulation in both the site potential and the inter-arm hopping integral is considered. The model mimics optical lattices where ultra-cold fermionic or bosonic atoms are trapped in double well potentials. It is observed that, the relative phase difference between the on-site potential and the inter-arm hopping integral, which can be controlled by the tuning of the interfering laser beams trapping the cold atoms, can result in a mixed spectrum of one or more absolutely continuous subband(s) and point like spectral measures. This opens up the possibility of a re-entrant metal–insulator transition. The subtle role played by the relative phase difference mentioned above is revealed, and we corroborate it numerically by working out the multi-channel electronic transmission for finite two-, and three-leg ladder networks. The extension of the calculation beyond the two-leg case is trivial, and is discussed in the work. - Graphical abstract: ▪ - Highlights: • Phase controlled metal–insulator transition is discussed. • An analytical prescription is given to understand MI transition. • Our work provides a way of designing experiments involving laser beams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Thermalization and light cones in a model with weak integrability breaking

Here, we employ equation-of-motion techniques to study the nonequilibrium dynamics in a lattice model of weakly interacting spinless fermions. Our model provides a simple setting for analyzing the effects of weak integrability-breaking perturbations on the time evolution after a quantum quench. We establish the accuracy of the method by comparing results at short and intermediate times to time-dependent density matrix renormalization group computations. For sufficiently weak integrability-breaking interactions we always observe prethermalization plateaus , where local observables relax to nonthermal values at intermediate time scales. At later times a crossover towards thermal behavior sets in. We determine the associated time scale, which depends on the initial state, the band structure of the noninteracting theory, and the strength of the integrability-breaking perturbation. Our method allows us to analyze in some detail the spreading of correlations and in particular the structure of the associated light cones in our model. We find that the interior and exterior of the light cone are separated by an intermediate region, the temporal width of which appears to scale with a universal power law t 1/3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Quasilocal charges and the generalized Gibbs ensemble in the Lieb-Liniger model

Here, we consider the construction of a generalized Gibbs ensemble composed of complete bases of conserved charges in the repulsive Lieb-Liniger model. We will show that it is possible to construct these bases with varying locality as well as demonstrating that such constructions are always possible provided one has in hand at least one complete basis set of charges. This procedure enables the construction of bases of charges that possess well-defined, finite expectation values given an arbitrary initial state. We demonstrate the use of these charges in the context of two different quantum quenches: a quench where the strength of the interactions in a one-dimensional gas is switched suddenly from zero to some finite value and the release of a one-dimensional cold atomic gas from a confining parabolic trap. While we focus on the Lieb-Liniger model in this paper, the principle of the construction of these charges applies to all integrable models, both in continuum and lattice form.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Derivative-free optimization of a rapid-cycling synchrotron

Here, we develop and solve a constrained optimization model to identify an integrable optics rapid-cycling synchrotron lattice design that performs well in several capacities. Our model encodes the design criteria into 78 linear and nonlinear constraints, as well as a single nonsmooth objective, where the objective and some constraints are defined from the output of Synergia, an accelerator simulator. We detail the difficulties of the 23-dimensional simulation-constrained decision space and establish that the space is nonempty. We use a derivative-free manifold sampling algorithm to account for structured nondifferentiability in the objective function. Our numerical results quantify the dependence of solutions on constraint parameters and the effect of the form of objective function.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗