Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice 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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Breakdown of order fractionalization in the CPT model

Here, we present an analysis of the half-filled CPT model, an analytically tractable Kondo lattice model with Yao-Lee spin-spin interactions on a 3D hyperoctagon lattice, proposed by Coleman, Panigrahi, and Tsvelik. Previous studies have established that the CPT model exhibits odd-frequency triplet superconductivity and order fractionalization. Through asymptotic analyses in the small- and large- Kondo coupling limits, we identify a quantum critical point at , marking a transition from a superconductor to a Kondo insulator. By estimating the vison gap energy to account for thermal gauge fluctuations, we determine the energy scales governing the thermal breakdown of order fractionalization. Moreover, at large the Kondo insulator undergoes orbital decoupling, leading to the formation of a decoupled Kitaev orbital liquid. These findings and analogies with the -gauged model lead us to propose a tentative phase diagram for the CPT model at half-filling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Emergent Superconductivity and Competing Charge Orders in Hole-Doped Square-Lattice t – J Model

The square-lattice Hubbard and closely related t-J models are considered as basic paradigms for understanding strong correlation effects and unconventional superconductivity (SC). Recent large-scale density matrix renormalization group simulations on the extended t-J model have identified d-wave SC on the electron-doped side (with the next-nearest-neighbor hopping t 2 > 0) but a dominant charge density wave (CDW) order on the hole-doped side (t 2 < 0), which is inconsistent with the SC of hole-doped cuprate compounds. We re-examine the ground-state phase diagram of the extended t-J model by employing the state-of-the-art density matrix renormalization group calculations with much enhanced bond dimensions, allowing more accurate determination of the ground state. On six-leg cylinders, while different CDW phases are identified on the hole-doped side for the doping range δ = 1/16 – 1/8, a SC phase emerges at a lower doping regime, with algebraically decaying pairing correlations and d-wave symmetry. On the wider eight-leg systems, the d-wave SC also emerges on the hole-doped side at the optimal 1=8 doping, demonstrating the winning of SC over CDW by increasing the system width. Furthermore, our results not only suggest a new path to SC in general t-J model through weakening the competing charge orders, but also provide a unified understanding on the SC of both hole- and electron-doped cuprate superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2024

Improved information criteria for Bayesian model averaging in lattice field theory

Bayesian model averaging is a practical method for dealing with uncertainty due to model specification. Use of this technique requires the estimation of model probability weights. Here, we revisit the derivation of estimators for these model weights. Use of the Kullback-Leibler divergence as a starting point leads naturally to a number of alternative information criteria suitable for Bayesian model weight estimation. We explore three such criteria, known to the statistics literature before, in detail: a Bayesian analog of the Akaike information criterion which we call the BAIC, the Bayesian predictive information criterion, and the posterior predictive information criterion (PPIC). We compare the use of these information criteria in numerical analysis problems common in lattice field theory calculations. We find that the PPIC has the most appealing theoretical properties and can give the best performance in terms of model-averaging uncertainty, particularly in the presence of noisy data, while the BAIC is a simple and reliable alternative.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Symmetry breaking in an extended O(2) model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term − h q cos ( q φ ) . Its coupling h q allows to smoothly interpolate between the O(2) model ( h q = 0 ) and a q -state clock model ( h q → ∞ ). In the latter case, a q -state clock model can also be defined for noninteger values of q . Thus, such a limit can also be considered as an analytic continuation of an ordinary q -state clock model to noninteger q . In previous work, we established the phase diagram for noninteger q in the infinite coupling limit ( h q → ∞ ). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger q , the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer q = 2 , 3, 4, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for q = 3 , 4 vary with the coupling, and for q = 4 the transition may turn into a Berezinskii-Kosterlitz-Thouless transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

Symmetry Breaking in an Extended-O(2) Model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we consider an extended-O(2) model that differs from the ordinary O(2) model by an explicit symmetry breaking term. Its coupling allows to smoothly interpolate between the O(2) model (zero coupling) and a $q$-state clock model (infinite coupling). In the latter case, a $q$-state clock model can also be defined for non-integer values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to non-integer $q$. The phase diagram of the extended-O(2) model in the infinite coupling limit was established in our previous work, where it was shown that for non-integer $q$, there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we investigate the model at finite values of the coupling using Monte Carlo and tensor methods. The results may be relevant for configurable Rydberg-atom arrays.

Hostetler, Leon↗
DOE OSTI2025

LSM and CPT

We study a number of 1+1d lattice models with anti-unitary symmetries that simultaneously reflect space and reverse time. Some of these symmetries are anomalous, leading to Lieb-Schultz-Mattis-type constraints, thus excluding a trivially gapped phase. Examples include a mod 8 anomaly in the Majorana chain and various mod 2 anomalies in the spin chain. In some cases, there is an exact, non-anomalous lattice symmetry that flows in the continuum to CPT. In some other cases, the CPT symmetry of the continuum theory is emergent or absent. Depending on the model, the anomaly of the lattice model is matched in the continuum in different ways. In particular, it can be mapped to an emergent anomaly of an emanant symmetry.

Anomalies in Field and String Theories↗
DOE OSTI2020

Theory of spin waves in a hexagonal antiferromagnet

We construct a field-theoretic description of spin waves in hexagonal antiferromagnets with three magnetic sublattices and coplanar 120° magnetic order. The three Goldstone modes can be separated by point-group symmetry into a singlet α 0 and a doublet β. The α 0 singlet is described by the standard theory of a free relativistic scalar field. The field theory of the β doublet is analogous to the theory of elasticity of a two-dimensional isotropic solid with distinct longitudinal and transverse “speeds of sound.” The well-known Heisenberg models on the triangular and kagome lattices with nearest-neighbor exchange turn out to be special cases with accidental degeneracy of the spin-wave velocities. The speeds of sound can be readily calculated for any lattice model. We apply this approach to the compounds of the Mn 3 X family with stacked kagome layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1985

A remark about pointed bubbles

The polymer expansion is a formal algebraic identity between a partition function and logarithm in statistical physics problems. The expansion gives a systematic method to control the free energy or to establish exponential tree-graph decay of connected correlations. Here, the convergence properties of the polymer expansion are analyzed in connection with three practical examples, including: intersecting bonds in chemical polymer chains; a connected closed hypersurface built from the (d-1)-faces of the d-dimensional unit cubes; and the set of Feynamn diagrams in the perturbation series of the Euclidean field theory partition function Z. The example of connected polymer chains is generalized to apply to other lattice models, including n-state Ising models at high temperature; short range lattice gases at high temperature; and weak coupling lattice field and gauge theories.

Garabedian, P. R.↗
DOE OSTI2024

Optically induced umklapp shift currents in striped cuprates

Motivated by recent experiments that observed low-frequency second-order optical responses in doped striped superconductors, here we investigate the nonlinear electrodynamics of systems exhibiting a charge density wave (CDW) order parameter. Due to the Bragg scattering off the CDW order, an incoming spatially homogeneous electric field in addition to zero momentum current generates umklapp currents that are modulated in space at momenta of the reciprocal CDW lattice. In particular, here we predict and microscopically evaluate the umklapp shift current, a finite momentum analog of the regular shift current which represents the second-order optical process that downconverts homogeneous AC electric field into low-frequency, zero momentum current. Specifically, we evaluate real-time response functions within mean-field theory via the Keldysh technique and use the Peierls substitution to compute observables at finite momenta in lattice models. We find that systems with certain lattice symmetries (such as inversion symmetry), where the regular shift current is disallowed, may give rise to the umklapp one. We apply our framework to investigate lattice symmetries in layered materials with helical-like stripes and show that both types of shift currents provide insight into the nature of intertwined phases of matter. Finally, we discuss the relation of our findings to recent experiments in striped superconductors. Published by the American Physical Society 2024

Materials Science↗
NASA NTRSDate not supplied

Exploring Quantum State Preparation Using Tensor Networks and Sparse Wavefunction Simulations

The variational quantum eigenvalue solver is a powerful hybrid quantum-classical approach that has been suggested as a candidate method to run on near-term quantum hardware for computing ground state electronic energies of molecular systems. However, even for small molecules, the number of variational parameters and qubits required to minimize the electronic energy is beyond the reach of current quantum computers except for small basis sets. We explore a new paradigm for state preparation where we test how much of the optimization can be approximately prepared with classical computers to reduce the number of optimization steps performed using a quantum device. By adapting a recent algorithm for the factorized form of the UCC ansatz, we can study molecular electronic structure problems with up to 64 qubits. In addition, we also test a related approach of using tensor networks to optimize quantum circuits in order to benchmark various lattice models. We present results using these approaches and discuss strategies for incorporating these ideas into variational algorithms involving near-term quantum computers. Our results help demonstrate the strength of the UCC ansatz and address pressing questions about optimal initial parameterizations and circuit construction.

quantum computing↗
DOE OSTI2020

Laser pulse driven control of charge and spin order in the two-dimensional Kondo lattice

Fast and dynamical control of quantum phases is highly desired in the application of quantum phenomena in devices. On the experimental side, ultrafast laser pulses provide an ideal platform to induce femtosecond dynamics in a variety of materials. Here we show that a laser pulse driven heavy fermion system can be tuned to dynamically evolve into a phase, which is not present in equilibrium. Using the state-of-the-art time-dependent variational Monte Carlo method, we perform numerical simulations of realistic laser pulses applied to the Kondo lattice model. By tracking spin and charge degrees of freedom we identify a dynamical phase transition, within the partially Kondo screened phase, from a charge ordered into a metallic state, while preserving the spin order of the system. Furthermore, we propose using high-harmonic generation and nonequilibrium optical conductivity to identify the dynamical phase transition in an ultrafast experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Importance learning estimator for the site-averaged turnover frequency of a disordered solid catalyst

For disordered catalysts such as atomically dispersed “single-atom” metals on amorphous silica, the active sites inherit different properties from their quenched-disordered local environments. The observed kinetics are site-averages, typically dominated by a small fraction of highly active sites. Standard sampling methods require expensive ab initio calculations at an intractable number of sites to converge on the siteaveraged kinetics. We present a new method that efficiently estimates the site-averaged turnover frequency (TOF). The new estimator uses the same importance learning algorithm [Vandervelden et al., React. Chem. Eng. 5, 77 (2020)] that we previously used to compute the siteaveraged activation energy. We demonstrate the method by computing the site-averaged TOF for a simple disordered lattice model of an amorphous catalyst. The results show that with the importance learning algorithm, the site-averaged TOF and activation energy can now be obtained concurrently with orders of magnitude reduction in required ab initio calculations.

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

Performance of the rigorous renormalization group for first-order phase transitions and topological phases

Expanding and improving the repertoire of numerical methods for studying quantum lattice models is an ongoing focus in many-body physics. While the density matrix renormalization group (DMRG) has been established as a practically useful algorithm for finding the ground state in one-dimensional systems, a provably efficient and accurate algorithm remained elusive until the introduction of the rigorous renormalization group (RRG) by Landau [Nat. Phys. 11, 566 (2015)1745-247310.1038/nphys3345]. In this paper, we study the accuracy and performance of a numerical implementation of RRG at first-order phase transitions and in symmetry-protected topological phases. Our study is motivated by the question of when RRG might provide a useful complement to the more established DMRG technique. In particular, despite its general utility, DMRG can give unreliable results near first-order phase transitions and in topological phases, since its local update procedure can fail to adequately explore (near-)degenerate manifolds. Furthermore, the rigorous theoretical underpinnings of RRG, meanwhile, suggest that it should not suffer from the same difficulties. We show this optimism is justified, and that RRG indeed determines well-ordered, accurate energies even when DMRG does not. Moreover, our performance analysis indicates that in certain circumstances seeding DMRG with states determined by coarse runs of RRG may provide an advantage over simply performing DMRG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Modeling of Wake-vortex Aircraft Encounters

There are more people passing through the world's airports today than at any other time in history. With this increase in civil transport, airports are becoming capacity limited. In order to increase capacity and thus meet the demands of the flying public, the number of runways and number of flights per runway must be increased. In response to the demand, the National Aeronautics and Space Administration (NASA), in conjunction with the Federal Aviation Administration (FAA), airport operators, and the airline industry are taking steps to increase airport capacity without jeopardizing safety. Increasing the production per runway increases the likelihood that an aircraft will encounter the trailing wake-vortex of another aircraft. The hazard of a wake-vortex encounter is that heavy load aircraft can produce high intensity wake turbulence, through the development of its wing-tip vortices. A smaller aircraft following in the wake of the heavy load aircraft will experience redistribution of its aerodynamic load. This creates a safety hazard for the smaller aircraft. Understanding this load redistribution is of great importance, particularly during landing and take-off. In this research wake-vortex effects on an encountering 10% scale model of the B737-100 aircraft are modeled using both strip theory and vortex-lattice modeling methods. The models are then compared to wind tunnel data that was taken in the 30ft x 60ft wind tunnel at NASA Langley Research Center (LaRC). Comparisons are made to determine if the models will have acceptable accuracy when parts of the geometry are removed, such as the horizontal stabilizer and the vertical tail. A sensitivity analysis was also performed to observe how accurately the models could match the experimental data if there was a 10% error in the circulation strength. It was determined that both models show accurate results when the wing, horizontal stabilizer, and vertical tail were a part of the geometry. When the horizontal stabilizer and vertical tail were removed there were difficulties modeling the sideforce coefficient and pitching moment. With the removal of only the vertical tail unacceptable errors occurred when modeling the sideforce coefficient and yawing moment. Lift could not be modeled with either the full geometry or the reduced geometry attempts.

Smith, Sonya T.↗
NASA NTRS2013

Coupled Vortex-Lattice Flight Dynamic Model with Aeroelastic Finite-Element Model of Flexible Wing Transport Aircraft with Variable Camber Continuous Trailing Edge Flap for Drag Reduction

This paper presents a coupled vortex-lattice flight dynamic model with an aeroelastic finite-element model to predict dynamic characteristics of a flexible wing transport aircraft. The aircraft model is based on NASA Generic Transport Model (GTM) with representative mass and stiffness properties to achieve a wing tip deflection about twice that of a conventional transport aircraft (10% versus 5%). This flexible wing transport aircraft is referred to as an Elastically Shaped Aircraft Concept (ESAC) which is equipped with a Variable Camber Continuous Trailing Edge Flap (VCCTEF) system for active wing shaping control for drag reduction. A vortex-lattice aerodynamic model of the ESAC is developed and is coupled with an aeroelastic finite-element model via an automated geometry modeler. This coupled model is used to compute static and dynamic aeroelastic solutions. The deflection information from the finite-element model and the vortex-lattice model is used to compute unsteady contributions to the aerodynamic force and moment coefficients. A coupled aeroelastic-longitudinal flight dynamic model is developed by coupling the finite-element model with the rigid-body flight dynamic model of the GTM.

Flexible Aicraft↗
NASA NTRS1995

Macroscopic and bulk-controlled elastic modes in an interaction of interstitial alcali metal cations within a face-centered cubic crystalline fullerine

The calculations in a majority of previous works for the fulleride (AqC-60) crystals were performed within the framework of the rigid-lattice model, neglecting the distortion relaxation of the host fullerene (C-60) crystal caused by the interstitial alkali-metal (A) cations. However, an each cation is a source of a static distortion field, and the resulting field is a superposition of such fields generated by all cations. This is a reason why the host-crystal distortions depend on the A-cations configurations, i.e. on a type of a spatial bulk distribution of interstitial cations. The given paper seeks to find a functional relation between the amplitudes of the doping-induced structure-distortion waves and of static concentration ones. A semiphenomenological model is constructed here within the scope of statistical-thermodynamic treatment and using the lattice-statistics simulation method(*). In this model the effects due to the presence of q solute A cations over available interstices (per unit cell) on the static inherent reorientation and/or displacements of the solvent molecules from the 'average-lattice' sites' as well as on the lattice parameter a of a elastically-anysotropic 'cubic' C-60 crystal are taken into account.

Tatarenko, Valentine A.↗