Exact Simulation Research Articles

Does the Correlation between 2MRS Galaxies and the CMB Indicate an Unmodeled CMB Foreground?

We revisit the claimed detection of a new cosmic microwave background (CMB) foreground based on the correlation between low-redshift Two Micron All Sky Survey Redshift Survey (2MRS) galaxies and CMB temperature maps from the Planck and Wilkinson Microwave Anisotropy Probe missions. We reproduce the reported measurements but argue that the original analysis significantly underestimated the uncertainties. We cross-correlate the 2MRS galaxy positions with simulated CMB maps and show that the correlation measured with the real data for late-type spiral galaxies at angular scales θ ≥ 0.°1 and redshift cz < 4500 km s−1 is consistent with zero at the 1.7σ level or less, depending on the exact CMB map and simulation construction. This was the sample that formed the basis for the original detection claim. For smaller angular separations the results are not robust to galaxy type or CMB cleaning method, and we are unable to draw firm conclusions. The original analysis did not propose a specific, falsifiable physical correlation mechanism, and it is impossible to rule out any contribution from an underlying physical effect. However, given our calculations, the lack of signal from expanding the redshift range, and the lack of corroboration from other galaxy surveys, we do not find the evidence for a new CMB foreground signal compelling.

Bath-induced interactions and transient dynamics in open quantum systems at strong coupling: Effective Hamiltonian approach.

Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.

Understanding Central Spin Decoherence Due to Interacting Dissipative Spin Baths.

We propose a new approach to simulate the decoherence of a central spin coupled to an interacting dissipative spin bath with cluster-correlation expansion techniques. We benchmark the approach on generic 1D and 2D spin baths and find excellent agreement with numerically exact simulations. Our calculations show a complex interplay between dissipation and coherent spin exchange, leading to increased central spin coherence in the presence of fast dissipation. Finally, we model near-surface nitrogen-vacancy centers in diamond and show that accounting for bath dissipation is crucial to understanding their decoherence. Our method can be applied to a variety of systems and provides a powerful tool to investigate spin dynamics in dissipative environments.

Fast Exact Simulation of the First Passage of a Tempered Stable Subordinator Across a Non-Increasing Function

Fast Exact Simulation of the First Passage of a Tempered Stable Subordinator Across a Non-Increasing Function

Efficient Large-Scale Many-Body Quantum Dynamics via Local-Information Time Evolution

During time evolution of many-body systems entanglement grows rapidly, limiting exact simulations to small-scale systems or small timescales. Quantum information tends, however, to flow towards larger scales without returning to local scales, such that its detailed large-scale structure does not directly affect local observables. This allows for the removal of large-scale quantum information in a way that preserves all local observables and gives access to large-scale and large-time quantum dynamics. To this end, we use the recently introduced to organize quantum information into different scales, allowing us to define and that we employ to systematically discard long-range quantum correlations in a controlled way. Our approach relies on decomposing the system into subsystems up to a maximum scale and time evolving the subsystem density matrices by solving the subsystem von Neumann equations in parallel. Importantly, the information flow needs to be preserved during the discarding of large-scale information. To achieve this without the need to make assumptions about the microscopic details of the information current, we introduce a second scale at which information is discarded, while using the state at the maximum scale to accurately obtain the information flow. The resulting algorithm, which we call local-information time evolution, is highly versatile and suitable for investigating many-body quantum dynamics in both closed and open quantum systems with diverse hydrodynamic behaviors. We present results for the energy transport in the mixed-field Ising model and the magnetization transport in the XX spin chain with onsite dephasing where we accurately determine the power-law exponent and the diffusion coefficients. Furthermore, the information lattice framework employed here promises to offer insightful results about the spatial and temporal behavior of entanglement in many-body systems. Published by the American Physical Society 2024

Simplifying the simulation of local Hamiltonian dynamics

Local Hamiltonians Hk describe nontrivial k-body interactions in quantum many-body systems. Here, we address the dynamical simulatability of a k-local Hamiltonian by a simpler one Hk′ with k′&lt;k, under the realistic constraint that both Hamiltonians act on the same Hilbert space. When it comes to exact simulation, we build upon known methods to derive examples of Hk and Hk′ that simulate the same physics. We also address the most realistic case of approximate simulation. There, we upper-bound the error up to which a Hamiltonian can simulate another one, regardless of their internal structure, and show an example suggesting that the accuracy of a (k′=2)-local Hamiltonian to simulate Hk with k&gt;2 increases with k. Finally, we propose a method to search for the k′-local Hamiltonian that simulates, with the highest possible precision, the short time dynamics of a given Hk Hamiltonian. Published by the American Physical Society 2024

Classical and quantum theory of fluctuations for many‐particle systems out of equilibrium

AbstractCorrelated classical and quantum many‐particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynamics, this is not the case for quantum systems. Alternatively, one can use many‐particle approaches such as hydrodynamics, kinetic theory, or nonequilibrium Green functions (NEGF). However, NEGF exhibit a very unfavorable cubic scaling of the CPU time with the number of time steps. An alternative is the G1–G2 scheme [N. Schlünzen et al., Phys. Rev. Lett. 124, 076601 (2020)] which allows for NEGF simulations with time linear scaling, however, at the cost of large memory consumption. The reason is the need to store the two‐particle correlation function. This problem can be overcome for a number of approximations by reformulating the kinetic equations in terms of fluctuations – an approach that was developed, for classical systems, by Yu.L. Klimontovich [JETP 33, 982 (1957)]. Here, we present an overview of his ideas and extend them to quantum systems. In particular, we demonstrate that this quantum fluctuations approach can reproduce the nonequilibrium GW approximation [E. Schroedter et al., Cond. Matt. Phys. 25, 23401 (2022)] promising high accuracy at low computational cost which arises from an effective semiclassical stochastic sampling procedure. We also demonstrate how to extend the approach to the two‐time exchange‐correlation functions and the density response properties. [E. Schroedter et al., Phys. Rev. B 108, 205109 (2023)]. The results are equivalent to the Bethe–Salpeter equation for the two‐time exchange‐correlation function when the generalized Kadanoff‐Baym ansatz with Hartree‐Fock propagators is applied [E. Schroedter and M. Bonitz, phys. stat. sol. (b) 2024, 2300564].

Hilbert Space Delocalization under Random Unitary Circuits.

The unitary dynamics of a quantum system initialized in a selected basis state yield, generically, a state that is a superposition of all the basis states. This process, associated with the quantum information scrambling and intimately tied to the resource theory of coherence, may be viewed as a gradual delocalization of the system's state in the Hilbert space. This work analyzes the Hilbert space delocalization under the dynamics of random quantum circuits, which serve as a minimal model of the chaotic dynamics of quantum many-body systems. We employ analytical methods based on the replica trick and Weingarten calculus to investigate the time evolution of the participation entropies which quantify the Hilbert space delocalization. We demonstrate that the participation entropies approach, up to a fixed accuracy, their long-time saturation value in times that scale logarithmically with the system size. Exact numerical simulations and tensor network techniques corroborate our findings.

Investigation of the Effect of Heat Transfer during Friction Stir Welding (FSW) of AZ80A Mg Alloy Plates using a Pin Tool by Conducting Finite Elements Analysis

Friction stir welding (FSW) is an innovative solid-state welding process that has attracted substantial attention due to its potential for combining problematic materials such magnesium alloys, such as AZ80A. In order to better understand the impact of heat transport during FSW of AZ80A magnesium alloy plates using a pin tool, this study used finite element analysis (FEA). The welding process's thermal features, such as temperature distribution, thermal stresses, and material flow patterns, are the major focus of this analysis. The first step of the study is to conduct a comprehensive literature evaluation to lay a firm groundwork and pinpoint knowledge gaps. The thermal conductivity, specific heat, density, and mechanical characteristics of AZ80A magnesium alloy are measured and recorded as part of the material characterisation process. To ensure an exact simulation of real-world welding circumstances, a comprehensive 3D model of the welding setup is built, including the AZ80A magnesium alloy plates and the pin tool. In order to accurately record temperature variations, a tiny mesh is used, particularly in the welding zone. By include boundary conditions that mimic the real-world welding characteristics, such as the rotation of the pin tool and the clamping or fixturing of the plates, finite element analysis is used to model the FSW procedure. To simulate the heat input produced by FSW, a heat source or heat production model is used.

Accuracy vs memory advantage in the quantum simulation of stochastic processes

Many inference scenarios rely on extracting relevant information from known data in order to make future predictions. When the underlying stochastic process satisfies certain assumptions, there is a direct mapping between its exact classical and quantum simulators, with the latter asymptotically using less memory. Here we focus on studying whether such quantum advantage persists when those assumptions are not satisfied, and the model is doomed to have imperfect accuracy. By studying the trade-off between accuracy and memory requirements, we show that quantum models can reach the same accuracy with less memory, or alternatively, better accuracy with the same memory. Finally, we discuss the implications of this result for learning tasks.

An atomic boson sampler.

A boson sampler implements a restricted model of quantum computing. It is defined by the ability to sample from the distribution resulting from the interference of identical bosons propagating according to programmable, non-interacting dynamics1. An efficient exact classical simulation of boson sampling is not believed to exist, which has motivated ground-breaking boson sampling experiments in photonics with increasingly many photons2-12. However, it is difficult to generate and reliably evolve specific numbers of photons with low loss, and thus probabilistic techniques for postselection7 or marked changes to standard boson sampling10-12 are generally used. Here, we address the above challenges by implementing boson sampling using ultracold atoms13,14 in a two-dimensional, tunnel-coupled optical lattice. This demonstration is enabled by a previously unrealized combination of tools involving high-fidelity optical cooling and imaging of atoms in a lattice, as well as programmable control of those atoms using optical tweezers. When extended to interacting systems, our work demonstrates the core abilities required to directly assemble ground and excited states in simulations of various Hubbard models15,16.

A remark on exact simulation of tempered stable Ornstein–Uhlenbeck processes

AbstractQu, Dassios, and Zhao (2021) suggested an exact simulation method for tempered stable Ornstein–Uhlenbeck processes, but their algorithms contain some errors. This short note aims to correct their algorithms and conduct some numerical experiments.

Exact Simulation of the Multifactor Ornstein–Uhlenbeck Driven Stochastic Volatility Model

Exact Simulation of the Multifactor Ornstein–Uhlenbeck Driven Stochastic Volatility Model

Numerically Exact Simulation of Photodoped Mott Insulators.

A description of long-lived photodoped states in Mott insulators is challenging, as it needs to address exponentially separated timescales. We demonstrate how properties of such states can be computed using numerically exact steady state techniques, in particular, the quantum Monte Carlo algorithm, by using a time-local ansatz for the distribution function with separate Fermi functions for the electron and hole quasiparticles. The simulations show that the Mott gap remains robust to large photodoping, and the photodoped state has hole and electron quasiparticles with strongly renormalized properties.

Quantum dynamics of a fully blockaded Rydberg atom ensemble

Classical simulation of quantum systems plays an important role in the study of many-body phenomena and in the benchmarking and verification of quantum technologies. Exact simulation is often limited to small systems because the dimension of the Hilbert space increases exponentially with the size of the system. For systems that possess a high degree of symmetry, however, classical simulation can reach much larger sizes. Here we consider an ensemble of strongly interacting atoms with permutation symmetry, enabling the computation of certain collective observables for hundreds of atoms at arbitrarily long evolution times. The system is realized by an ensemble of three-level atoms, where one of the levels corresponds to a highly excited Rydberg state. In the limit of all-to-all Rydberg blockade, the Hamiltonian is invariant under permutation of the atoms. Using techniques from representation theory, we construct a block-diagonal form of the Hamiltonian, where the size of the largest block increases only linearly with the system size. We apply this formalism to derive efficient pulse sequences to prepare arbitrary permutation-invariant quantum states. Moreover, we study the quantum dynamics following a quench, uncovering a parameter regime in which the system thermalizes slowly and exhibits pronounced revivals. Our results create opportunities for the experimental and theoretical study of large interacting and nonintegrable quantum systems. Published by the American Physical Society 2024

Exact simulation of the Hull and White stochastic volatility model

We show how to simulate exactly the asset price and the variance under the Hull and White stochastic volatility model. We derive analytical formulas for the Laplace transform of the time integral of volatility conditional on the variance level at the endpoint of the time interval and the Laplace transform of integrated variance conditional on both integrated volatility and variance. Based on these results, we simulate the model through a nested-conditional factorization approach, where Laplace transforms are inverted through the (conditional) Fourier-cosine (COS) method. Under this model, our approach can be used to generate unbiased estimates for the price of derivatives instruments. We propose some variants of the exact simulation scheme for computing unbiased estimates of option prices and sensitivities, a difficult task in the Hull and White model. These variants also allow for a significant reduction in the Monte Carlo simulation estimator's variance (around 93-98%) and the computing time (around 22%) when pricing options. The performances of the proposed algorithms are compared with various benchmarks. Numerical results demonstrate the faster convergence rate of the error in our method, which achieves an O(s−1/2) convergence rate, where s is the total computational budget, largely outperforming the benchmark.

Benchmarking highly entangled states on a 60-atom analogue quantum simulator

Quantum systems have entered a competitive regime in which classical computers must make approximations to represent highly entangled quantum states1,2. However, in this beyond-classically-exact regime, fidelity comparisons between quantum and classical systems have so far been limited to digital quantum devices2–5, and it remains unsolved how to estimate the actual entanglement content of experiments6. Here, we perform fidelity benchmarking and mixed-state entanglement estimation with a 60-atom analogue Rydberg quantum simulator, reaching a high-entanglement entropy regime in which exact classical simulation becomes impractical. Our benchmarking protocol involves extrapolation from comparisons against an approximate classical algorithm, introduced here, with varying entanglement limits. We then develop and demonstrate an estimator of the experimental mixed-state entanglement6, finding our experiment is competitive with state-of-the-art digital quantum devices performing random circuit evolution2–5. Finally, we compare the experimental fidelity against that achieved by various approximate classical algorithms, and find that only the algorithm we introduce is able to keep pace with the experiment on the classical hardware we use. Our results enable a new model for evaluating the ability of both analogue and digital quantum devices to generate entanglement in the beyond-classically-exact regime, and highlight the evolving divide between quantum and classical systems.

Composite Qdrift-product formulas for quantum and classical simulations in real and imaginary time

Recent study has shown that it can be advantageous to implement a composite channel that partitions the Hamiltonian H for a given simulation problem into subsets A and B such that H=A+B, where the terms in A are simulated with a Trotter-Suzuki channel and the B terms are randomly sampled via the Qdrift algorithm. Here we extend Qdrift and composite product formulas to imaginary time, formulating candidate classical algorithms for quantum Monte Carlo calculations. We upper bound the induced Schatten-1→1 norm on both imaginary-time Qdrift and composite channels. Another recent result demonstrated that simulations of lattice Hamiltonians containing geometrically local interactions can be improved using a Lieb-Robinson argument to decompose H into subsets that contain only terms supported on that subset of the lattice. Here, we provide a quantum algorithm by unifying this result with the composite approach into “local composite channels” and we upper bound the diamond distance. We provide exact numerical simulations of algorithmic cost by counting the number of gates of the form e−iHjt and e−Hjβ to meet a certain error tolerance ε. In doing so, we optimize the partitioning into sets A and B using gradient boosted tree models from machine learning. These numerical studies are important given that product formulas have been historically known to outperform analytic upper bounds. We show constant factor advantages for a variety of interesting Hamiltonians, the maximum of which is a ≈20-fold speedup that occurs in the simulation of Jellium. Published by the American Physical Society 2024

The dilute 3D and three site dimer Heisenberg antiferromagnet: review of some magnetic properties.

Many properties of the spin–½ kagome Heisenberg antiferromagnet remains controversial and receiving continuous attention from theorists and experimentalists due to their fascinating magnetic properties and high heterogeneity coax by the geometrical frustration. In this research we seek and compare the results obtained in our recent article [1] with that of Napel[2] where we elucidate some of the identified magnetic behaviour using the spin–½ Heisenbeg model with the application of Exact Diagonalization and classical Monte-Carlo simulation on three site dimer and 3D simple cubic lattice simulation. The results obtained shows excellent agreement in both papers reviewed with [1] establishing firmly the FM to AFM transition, having a clear transition point referred to mixed phased [2].

Exact simulation of classical heat engine cycles using single-ion phonon laser

Heat engines are essential devices in modern industry, converting heat energy into useful mechanical work via their working substances. Here we experimentally simulate the conventional heat engine by employing the vibrational mode of a single trapped ion as the working substance. In contrast to simply employing the ion in thermal motion, we consider coherently stimulating the ion’s vibrational motion as the phonon laser, which helps acquire clearer results by effectively suppressing the thermal fluctuation. As such, we demonstrate in an exact and high signal-to-noise way the standard steps of both the Otto and Carnot cycles in a single ion, and compare their maximum efficiencies by monitoring the amplitude and frequency of the vibration. Our work witnesses an interesting single-atom thermal engine using coherently controlled phonons. It would be the smallest platform for simulating or demonstrating classical thermodynamic laws and phenomena at a single ion scale via optical manipulation techniques for phonon lasers.