Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analog quantum simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Chemist’s Guide to Solid-State Qubits: Diamond Color Centers as Embedded Vacancy-Centered Molecules

Diamond point defects or color centers are at the crux of solid-state quantum technology applications, from computing to sensing. Modeling molecular qubits has advanced by leaps and bounds mainly due to the advanced quantum mechanical (QM) tools in the arsenal of quantum chemists for modeling (isolated) molecular systems, which should be easily translatable for analogous systems in the solid state. Here, we simulate a series of diamond defects: a carbon vacancy (V C ) and second row substitutions that are adjacent to a V C or D C V C , with D = B, N, and O, to evaluate their magnetic and optical properties. We applied molecular orbital principles, an embedded cluster method, and a multiconfigurational QM theory to computationally characterize the above-mentioned defects for a multitude of spin, charge, and excited states. Our computational approach delivered high-quality QM insights that compare well with experiments. Furthermore, with this framework, we confirm the spin-triplets B C V C – and O C V C 0 to be analogous to N C V C – , with calculated bright transitions that respectively correspond to UV (3.80 eV) and blue (2.65 eV) light absorption and emission lifetimes (τ) of 2.3 and 5.6 ns, which are higher in energy and more short-lived (and thus brighter) than that for N C V C – : green (2.26 eV), τ = 9.5 ns.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2025

Lower bounds on entanglement entropy without twin copy

We discuss the possibility of estimating experimentally the von Neumann entanglement entropy S A v N of a symmetric bipartite quantum system A B by using the basic measurement counts (bitstrings) for a single copy of a prepared state. Using exact diagonalization and analog simulations performed with the publicly available QuEra facilities for chains and ladders of Rydberg atoms, we calculate the Shannon entropy S A B X associated with the bitstrings of adiabatically prepared ground states and the reduced entropies S A X and S B X obtained from the marginal probabilities in A and B . We then calculate the classical mutual information I A B X = S A X + S B X − S A B X , which is a lower bound on S A v N . We show that for a broad range of lattice spacing and detuning, I A B X is typically 20% below S A v N in regions where S A v N is large and a less close bound in regions where S A v N is low. We argue that this use of the easily available bitstrings provides a robust and efficient way to explore empirically the phase diagram of qubit-based quantum simulators and identify critical regions. Published by the American Physical Society 2025

Meurice, Yannick (ORCID:0000000209959694)↗
NASA NTRS2010

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗
DOE OSTI2024

Report for the ASCR Workshop on Basic Research Needs in Quantum Computing and Networking - 2023

Employing quantum mechanical resources in computing and networking opens the door to new computation and communication models and potential disruptive advantages over classical counterparts. However, quantifying and realizing such advantages face extensive scientific and engineering challenges. Investments by the Department of Energy (DOE) have driven progress toward addressing such challenges. Quantum algorithms have been recently developed, in some cases offering asymptotic exponential advantages in speed or accuracy, for fundamental scientific problems such as simulating physical systems, solving systems of linear equations, or solving differential equations. Empirical demonstrations on nascent quantum hardware suggest better performance than classical analogs on specialized computational tasks favorable to the quantum computing systems. However, demonstration of an end-to-end, substantial and rigorously quantifiable quantum performance advantage over classical analogs remains a grand challenge, especially for problems of practical value. The definition of requirements for quantum technologies to exhibit scalable, rigorous, and transformative performance advantages for practical applications also remains an outstanding open question, namely, what will be required to ultimately demonstrate practical quantum advantage?

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Programmable exploration of magnetic states in Lieb-kagome interpolated lattices

We investigate a hybrid modeling framework in which a quantum annealer is used to simulate magnetic interactions in molecular qubit lattices inspired by experimentally realizable systems. Using phthalocyanine assemblies as a structurally constrained prototype, we model a continuous deformation from a Lieb to a kagome lattice, revealing frustration-driven disorder and magnetic field-induced reordering in the spin structure. Here, the goal is to show how a quantum annealer can operate as a physically instantiated, programmable platform to emulate experimentally relevant lattice deformations and produce observables in a manner analogous to an experimental measurement, enabling the characterization of magnetic arrangements beyond the reach of current molecular architectures. This surrogate modeling approach offers a pathway to explore and iteratively design tunable magnetic states in synthetic materials. The synthetic design, structural characterization, and quantum simulation framework established here defines a modular and scalable paradigm for probing the limits of engineered matter across chemistry, condensed matter, and quantum information science.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Lieb-Mattis States for Robust Entangled Differential Phase Sensing

We explore a two-node, entanglement-enhanced sensor network for differential phase sensing that exploits decoherence-free subspaces to suppress common-mode noise, a primary limitation of many state-of-the-art quantum sensors. We identify a class of entangled states that, while not strictly optimal, achieve the same asymptotic sensitivity scaling as optimal states and can be prepared efficiently from initially unentangled atomic ensembles. Importantly, the preparation time decreases with increasing system size, which makes the states compatible with realistic noise processes in present-day quantum sensors that operate with large particle numbers but lack full error correction. We illustrate these ideas using two cavity-mediated preparation protocols: (i) coherent, unitary entanglement generation analogous to bosonic two-mode squeezing, yielding Heisenberg scaling, and (ii) dissipative preparation through collective emission into a shared cavity mode, providing a square-root improvement beyond the standard quantum limit. Numerical simulations show that both approaches remain effective at experimentally realistic cavity cooperativities, establishing a practical path toward scalable, quantum-enhanced differential phase sensing.

Kaubruegger, Raphael [National Institute of Standa↗
DOE OSTI2025

Ancilla-entangling Floquet kicks for accelerating quantum algorithms

Quantum simulation with adiabatic annealing can provide insight into difficult problems that are impossible to study with classical computers. However, it deteriorates when the systems scale up due to the shrinkage of the excitation gap and thus places an annealing rate bottleneck for high success probability. Here, in this study, we accelerate quantum simulation using digital multiqubit gates that entangle primary system qubits with ancillary qubits. The practical benefits originate from tuning the ancillary gauge degrees of freedom to enhance the quantum algorithm's original functionality in the system registry. For simple but nontrivial short-ranged, infinite long-ranged transverse-field Ising models, and the hydrogen molecule model after qubit encoding, we show improvement in the time to solution by one hundred percent but with higher accuracy through exact state-vector numerical simulation in a digital-analog setting. The findings are further supported by time-averaged Hamiltonian theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Studying Open Quantum Systems Relevant to Chemistry on a Trapped-Ion Quantum Simulator (Final Technical Report)

This project advances the trapped-ion quantum simulator as a versatile platform for studying open quantum system phenomena. We aim to contribute to the emerging quantum simulation toolkits and enable simulation of nanoscale energy processes. Trapped-ion platforms offer unique capabilities: their vibrational motion can be precisely manipulated, measured, and coherently coupled to auxiliary qubits. The vibrational mode can function both as a highly sensitive quantum sensor and a programmable environment bath. Using this platform, we achieved three major outcomes. First, we demonstrated using the vibrational mode as an ultrasensitive probe for testing fundamental physics, including possible nonlinear quantum mechanics effects. Second, we established that these modes can act as controllable baths in which tunable noise and loss can enhance or modify energy-transfer dynamics, providing the experimental preparation toward studying mechanisms relevant to chemical reactions and light-harvesting systems. Third, by introducing controllable nonlinear gain and loss, we showed theoretically how simulations using trapped ions can model vibrationally-assisted energy transport in a non‐Hermitian quantum system comprising a chromophore dimer weakly coupled to a vibrational mode. Exploring the non‐Hermitian dynamics of the whole system including vibrations, we found that energy transfer accompanied by absorption of phonons from a vibrational mode can be significantly enhanced near an exceptional point. This theoretical work on simulation of energy transfer processes in driven non‐Hermitian quantum systems revealed an interesting novel path to study open quantum systems dynamics under conditions of gain and loss. We then further explored the benefits of controllable gain and loss with an experimental realization of quantum analogs of nonlinear oscillators, namely, the van der Pol oscillator. Here we observed mutual synchronization mediated by collective dissipation between two oscillators. In parallel, we explored related quantum networking protocols using the same trapped-ion platform, developing fast, high-fidelity schemes for distributing entanglement. Together, these achievements show that trapped-ion vibrational modes provide a highly programmable and high-fidelity platform for investigating complex dissipative quantum behavior, while enabling new approaches to remote quantum sensing, energy science, and nonlinear quantum dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantifying the impact of precision errors on quantum approximate optimization algorithms

The quantum approximate optimization algorithm (QAOA) is a hybrid quantum-classical algorithm that seeks to achieve approximate solutions to optimization problems by iteratively alternating between intervals of controlled quantum evolution. Here, we examine the effect of analog precision errors on QAOA performance from the perspective of both algorithmic training and performance guarantees. Leveraging cumulant expansions, we recast the faulty QAOA as a control problem in which precision errors are expressed as multiplicative control noise and derive bounds on the performance of QAOA. We show using both analytical techniques and numerical simulations that fixed precision implementations of QAOA circuits are subject to an exponential degradation in performance dependent upon the number of optimal QAOA layers and magnitude of the precision error. Despite this significant reduction, we show that it is possible to mitigate precision errors in QAOA via digitization of the variational parameters at the cost of increasing circuit depth.

quantum algorithms↗
DOE OSTI2026

Analog-to-digital converter based on voltage-controlled superconducting devices

The increasing demand for cryogenic electronics in superconducting and quantum computing systems calls for ultra-energy-efficient data conversion architectures that remain functional at deep cryogenic temperatures. Here, in this work, we present the first design of a voltage-controlled superconducting flash analog-to-digital converter (ADC) based on a voltage-controlled quantum-enhanced Josephson junction field-effect transistor (JJFET). Exploiting its strong gate tunability and transistor-like behavior, the JJFET offers a scalable alternative to conventional current-controlled superconducting devices while aligning naturally with CMOS-style design methodologies. Building on our previously developed Verilog-A compact model calibrated to experimental data, we design and simulate a three-bit JJFET-based flash ADC targeted for integration within cryogenic control and readout circuitry in quantum computing. The core comparator block is realized through careful bias current selection and augmented with a three-terminal nanocryotron to precisely define reference voltages. Cascaded JJFET comparators ensure robust voltage gain, cascadability, and logic-level restoration across stages. Simulation results demonstrate accurate quantization behavior with ultra-low power dissipation, underscoring the feasibility of voltage-driven superconducting mixed-signal circuits. This work establishes a critical step toward unifying superconducting logic and data conversion, paving the way for scalable cryogenic architectures in quantum–classical co-processors, low-power artificial intelligence accelerators, and next-generation energy-constrained computing platforms.

Analog-to-digital converter↗
DOE OSTI2024

On the connection between least squares, regularization, and classical shadows

Classical shadows (CS) offer a resource-efficient means to estimate quantum observables, circumventing the need for exhaustive state tomography. Here, we clarify and explore the connection between CS techniques and least squares (LS) and regularized least squares (RLS) methods commonly used in machine learning and data analysis. By formal identification of LS and RLS ``shadows'' completely analogous to those in CS---namely, point estimators calculated from the empirical frequencies of single measurements---we show that both RLS and CS can be viewed as regularizers for the underdetermined regime, replacing the pseudoinverse with invertible alternatives. Through numerical simulations, we evaluate RLS and CS from three distinct angles: the tradeoff in bias and variance, mismatch between the expected and actual measurement distributions, and the interplay between the number of measurements and number of shots per measurement. Compared to CS, RLS attains lower variance at the expense of bias, is robust to distribution mismatch, and is more sensitive to the number of shots for a fixed number of state copies---differences that can be understood from the distinct approaches taken to regularization. Conceptually, our integration of LS, RLS, and CS under a unifying ``shadow'' umbrella aids in advancing the overall picture of CS techniques, while practically our results highlight the tradeoffs intrinsic to these measurement approaches, illuminating the circumstances under which either RLS or CS would be preferred, such as unverified randomness for the former or unbiased estimation for the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Josephson dynamics in two-dimensional ring-shaped condensates

We investigate Josephson transport in a fully closed, two-dimensional superfluid circuit formed by a ring-shaped 87 Rb Bose-Einstein condensate that contains two optical barriers acting as movable weak links. Translating these barriers at controlled speeds imposes a steady bias current, enabling direct mapping of the current-chemical-potential (𝐼−Δ⁢𝜇) characteristics. For narrow junctions (𝑤 ≈ 1µ⁢m), the circuit exhibits a pronounced dc branch that terminates at a critical current 𝐼 𝑐 = 9⁢(1) × 10 3 s −1 ; above this threshold, the system switches to an ac, resistive regime. Classical-field simulations that include the moving barriers quantitatively reproduce both the nonlinear 𝐼−Δ⁢𝜇 curve and the measured 𝐼 𝑐 , validating the underlying microscopic picture. Analysis of the ensuing phase dynamics shows that dissipation is mediated by the nucleation and traversal of vortex-antivortex pairs through the junctions, while the bulk condensate remains globally phase locked—direct evidence of the ring's topological constraint enforcing quantized circulation. These results establish a cold-atom analog of a superconducting quantum interference device in which Josephson dynamics can be resolved at the single-vortex level, providing a versatile platform for atomtronic circuit elements, nonreciprocal Josephson devices, and on-chip Sagnac interferometers for multiaxis rotation sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Interaction Effects on the Dynamical Anderson Metal-Insulator Transition Using Kicked Quantum Gases

Understanding the interplay of interaction and disorder in quantum transport poses long-standing scientific challenges for theory and experiment. While highly controlled ultracold atomic platforms combining atomic interactions with spatially disordered lattices have led to remarkable advances, the extension of such controlled studies to phenomena in high-dimensional disordered systems, such as the three-dimensional Anderson metal-insulator transition has been limited. Kicked quantum gases provide an alternate experimental platform that captures the Anderson model in momentum space and features dynamical localization as the analog of Anderson localization. Here, we utilize a momentum space lattice platform using quasiperiodically kicked ultracold atomic gases to experimentally investigate interaction effects on the three-dimensional dynamical Anderson metal-insulator transition. Here, we observe interaction-driven subdiffusion and a divergence of delocalization onset time on approaching the phase boundary. Mean-field numerical simulations show qualitative agreement with experimental observations, but with significant quantitative deviations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Third-order photon correlations extract single-nanocrystal multiexciton properties in solution

Colloidal semiconductor nanocrystals are considered promising materials for high-flux optical applications, including lasing, light-emitting diodes, biological imaging, and quantum optics. In high-flux applications, multiexcitons can significantly contribute to emission, influencing its brightness, spectral purity, and kinetics. As a result, understanding and controlling multiexciton emission in colloidal nanocrystal materials is of the utmost importance. In the past, single-nanocrystal photon correlation methods have been applied to understand biexciton and triexciton efficiencies, lifetimes, and spectra. While powerful, such methods suffer from user selection bias and require stable emission from single nanocrystals. To compensate for this shortcoming, second-order correlation methods were developed to extract sample-averaged biexciton properties from a solution of nanocrystals. Until now, however, the analogous third-order solution photon correlation methods remained unexplored. In this work, we present a pair of third-order photon correlation techniques to obtain the sample-averaged single-nanocrystal triexciton quantum yield and lifetime in a solution-phase experiment. These techniques derive from the relationship between the Poisson probability of nanocrystal photon absorption and the intrinsic probability of nanocrystal photon emission. We validate the theoretical background of these techniques by creating a numerical model to simulate the diffusion and emission of many nanocrystals in solution. Our simulations confirm that the average triexciton quantum yield and triexciton lifetime can be extracted from a solution of nanocrystals. These techniques will enable researchers to gain a better understanding of the fundamental multiexciton properties of colloidal nanocrystals.

Horowitz, Jonah R. [Massachusetts Institute of Tec↗
DOE OSTI2026

Emulating 2D Materials with Magnons

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

Ferrimagnets↗
NASA NTRS2012

Analog Exercise Hardware to Implement a High Intensity Exercise Program During Bed Rest

Background: In order to evaluate novel countermeasure protocols in a space flight analog prior to validation on the International Space Station (ISS), NASA's Human Research Program (HRP) is sponsoring a multi-investigator bedrest campaign that utilizes a combination of commercial and custom-made exercise training hardware to conduct daily resistive and aerobic exercise protocols. This paper will describe these pieces of hardware and how they are used to support current bedrest studies at NASA's Flight Analog Research Unit in Galveston, TX. Discussion: To implement candidate exercise countermeasure studies during extended bed rest studies the following analog hardware are being utilized: Stand alone Zero-Gravity Locomotion Simulator (sZLS) -- a custom built device by NASA, the sZLS allows bedrest subjects to remain supine as they run on a vertically-oriented treadmill (0-15 miles/hour). The treadmill includes a pneumatic subject loading device to provide variable body loading (0-100%) and a harness to keep the subject in contact with the motorized treadmill to provide a ground reaction force at their feet that is quantified by a Kistler Force Plate. Supine Cycle Ergometer -- a commercially available supine cycle ergometer (Lode, Groningen, Netherlands) is used for all cycle ergometer sessions. The ergometer has adjustable shoulder supports and handgrips to help stabilize the subject during exercise. Horizontal Squat Device (HSD) -- a custom built device by Quantum Fitness Corp (Stafford, TX), the HSD allows for squat exercises to be performed while lying in a supine position. The HSD can provide 0 to 600 pounds of force in selectable 5 lb increments, and allows hip translation in both the vertical and horizontal planes. Prone Leg Curl -- a commercially available prone leg curl machine (Cybex International Inc., Medway, MA) is used to complete leg curl exercises. Horizontal Leg Press -- a commercially available horizontal leg press (Quantum Fitness Corporation) is used for leg press and heel raise exercises. Minor modifications were made to the device including adding 200 lbs to the weight stack, raising the frame by 12 inches, making the footplate adjustable, and providing removable handles. Conclusion: A combination of novel and commercial exercise hardware are used to mimic the exercise hardware capabilities aboard the ISS, allowing scientific investigation of new countermeasure protocols in a space flight analog prior to flight validation

Loerch, Linda↗
DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗