Analog Quantum Simulation

Analog quantum simulation uses the continuously acting, engineered dynamics of one quantum system to reproduce selected states, observables, or processes of a target quantum model. Instead of decomposing a target evolution into a universal gate sequence, the experiment arranges physical degrees of freedom and controls so that its effective Hamiltonian or open-system generator belongs to the desired model family.

The word analog does not mean uncontrolled, classical, or merely qualitative. Modern analog simulators can have site-resolved preparation, programmable geometry, time-dependent fields, randomized measurements, and closed-loop calibration. What distinguishes them is that a many-body interaction acts as a native physical process rather than as a sequence of logically specified elementary gates. This often gives large system sizes and direct access to dynamics, at the cost of a restricted model family and no generic error-correction layer between the target and the apparatus.

Required background. What Is Quantum Simulation? introduces model correspondence; unitary evolution and density operators supply the propagators and expectation values used below.

Target observables and laboratory predictions

An analog-simulation claim should identify more than a Hamiltonian name. A useful task specification is

$$\mathfrak A = \left( \mathcal H_{\mathrm{tar}}, \Omega_{\mathrm{state}}, \Omega_{\mathrm{obs}}, \Theta, \mathcal T, \epsilon, \delta \right).$$

Here:

• $\mathcal H_{\mathrm{tar}}$ is the target Hilbert space;
• $\Omega_{\mathrm{state}}$ is the family of target inputs to be represented;
• $\Omega_{\mathrm{obs}}$ is the family of target observables to be estimated;
• $\Theta$ is the target parameter region;
• $\mathcal T$ is the set of target times or preparation protocols;
• $\epsilon$ is the promised accuracy for the reported quantities;
• $\delta$ is the allowed statistical failure probability.

Restricting the state and observable families is not a defect. It is often the reason a specialized simulator can be useful. A device may faithfully reproduce local spin correlations following a quench without representing every state in the target Hilbert space or permitting arbitrary global measurements. The restriction must, however, be stated before examining the result.

Let the completely positive trace-preserving map $\mathcal E_{\mathrm{s}}$ encode target states into the laboratory system, let $\mathcal T_\tau^{\boldsymbol\lambda}$ be the laboratory evolution at control setting $\boldsymbol\lambda$, and let $\mathcal E_{\mathrm{o}}(O)$ be the laboratory observable associated with $O$. If target parameters are related to controls by $\boldsymbol\theta=f(\boldsymbol\lambda)$ and target time is $t=\alpha\tau$ with $\alpha>0$, the operational prediction is

$$\widetilde\mu_O(t) = \operatorname{Tr} \left[ \mathcal E_{\mathrm{o}}(O) \, \mathcal T_{t/\alpha}^{\boldsymbol\lambda} \left( \mathcal E_{\mathrm{s}}(\rho) \right) \right].$$

The ideal target value is

$$\mu_O(t) = \operatorname{Tr} \left[ O\, \Phi_t^{\boldsymbol\theta}(\rho) \right].$$

For a bounded observable family (for example, $\lVert O\rVert\leq1$), an observable-specific analog guarantee can therefore be written

$$\sup_{\substack{ \rho\in\Omega_{\mathrm{state}}\\ O\in\Omega_{\mathrm{obs}}\\ \boldsymbol\theta\in\Theta,\, t\in\mathcal T }} \left| \widetilde\mu_O(t)-\mu_O(t) \right| \le \epsilon.$$

This form makes three facts explicit. The guarantee can be local rather than global, it applies only over a declared regime, and it includes preparation and measurement as well as evolution. A statement that a device “implements the Ising model” is incomplete until these domains and interfaces are given.

The Hamiltonian Correspondence

For closed-system simulation, the central reduction often takes the form

$$P H_{\mathrm{lab}}(\boldsymbol\lambda) P = \alpha(\boldsymbol\lambda) V H_{\mathrm{tar}} \left( \boldsymbol\theta(\boldsymbol\lambda) \right) V^\dagger + \beta(\boldsymbol\lambda)P + \Delta H(\boldsymbol\lambda).$$

The ingredients have distinct meanings:

Symbol	Role
$V$	isometric encoding from the target Hilbert space into the simulator
$P=VV^\dagger$	projector onto the intended encoded subspace
$\alpha$	conversion between target and laboratory energy scales
$\beta$	constant energy shift within the encoded subspace
$\Delta H$	in-subspace mismatch: omitted terms, disorder, drift, and calibration error
$Q=I-P$	leakage sector outside the intended description

The projected equation is not enough by itself. Couplings $QH_{\mathrm{lab}}P$ can move population out of the encoded space, and virtual excursions into $Q$ can modify the effective Hamiltonian even when the final leakage population is small. The model reduction should therefore state the spectral separation, coupling scale, truncation order, and regime over which $P$ remains meaningful.

Time and energy rescaling

For fixed controls and time-independent Hamiltonians, suppose $\Delta H=0$ and the encoded subspace is invariant: $[H_{\mathrm{lab}},P]=0$. Because $PV=V$, one has

$$H_{\mathrm{lab}}V = \alpha V H_{\mathrm{tar}} + \beta V.$$

Expanding the exponential or using functional calculus gives

$$e^{-iH_{\mathrm{lab}}\tau/\hbar}V = e^{-i\beta\tau/\hbar} V e^{-iH_{\mathrm{tar}}t/\hbar}, \qquad t=\alpha\tau.$$

The $\beta$ term contributes a global phase to an encoded pure state and cancels from density-operator evolution. The scale $\alpha$ is operationally important: it converts laboratory coherence time into reachable target time, laboratory thermal energy $k_{\mathrm B}T_{\mathrm{lab}}$ into target thermal energy $k_{\mathrm B}T_{\mathrm{lab}}/\alpha$ when the prepared state is a Gibbs state of this encoded Hamiltonian, and frequency uncertainty into target-parameter uncertainty. A dimensionless temperature still requires division by a target coupling or other target energy scale.

The same mapping must hold for observables. In the ideal encoded subspace,

$$P O_{\mathrm{lab}} P = V O_{\mathrm{tar}}V^\dagger + \Delta O.$$

A correct Hamiltonian with an uncalibrated detector does not yield a correct target observable. Conversely, a detector may access a useful target observable even when it cannot reconstruct the full state.

Effective Hamiltonians are conditional

Many analog mappings are not exact microscopic identities. They are low-energy, rotating-frame, high-frequency, dispersive, or adiabatically eliminated descriptions. Schematically, if $H_0$ separates $P$ and $Q$ by a gap $\Delta_{\mathrm{sep}}$ and $W$ couples them, eliminating the $Q$ component of an eigenvector gives an energy-dependent reduced operator. Expanding its resolvent gives

$$H_{\mathrm{eff}}(E) = P(H_0+W)P - PWQ \left( QH_0Q-E \right)^{-1} QWP + \cdots.$$

The correction is typically of order $\lVert W\rVert^2/\Delta_{\mathrm{sep}}$, but coefficients and signs depend on the convention and spectrum. Schrieffer–Wolff Transformation instead constructs an energy-independent block-diagonal Hamiltonian order by order, with a dressed encoding. States and measured operators must be transformed with the same dressing. Adiabatic Elimination and Average Hamiltonian Theory explain two other common reductions.

The reduction needs a small expansion parameter, a specified subspace, and a time range over which omitted terms remain small. A virtual coupling can produce a second-order energy shift even when its excited-level population stays small. Exercise 8 makes this distinction explicit.

Native Dynamics and Restricted Programmability

A convenient control model is

$$H_{\mathrm{lab}} \left( \boldsymbol\lambda,t \right) = H_0 + \sum_{a=1}^{m} u_a(\boldsymbol\lambda,t)H_a + H_{\mathrm{parasitic}} \left( \boldsymbol\lambda,t \right).$$

The drift $H_0$ and control generators $H_a$ define a physically reachable family. Geometry, laser frequency, magnetic field, interaction dressing, coupler bias, and trap depth may tune that family. They rarely tune every coefficient independently. One optical intensity can change a desired coupling, a light shift, and a scattering rate together; one geometric move can change several pair interactions at once.

An analog device is programmable when controls select a meaningful family of models, schedules, states, or observables. The term does not imply gate-model universality or arbitrary $k$-local Hamiltonian synthesis. Reports should distinguish:

• fixed: one Hamiltonian with a small number of global settings;
• tunable: continuous control of several global parameters;
• reconfigurable: geometry or coupling graph can be changed between runs;
• locally programmable: site- or bond-dependent parameters can be set;
• time programmable: ramps, quenches, Floquet drives, or feedback are available;
• universal: an explicit construction approximates a broad computational or Hamiltonian class with controlled resources.

These capabilities are different axes of control, not automatically a nested hierarchy; accuracy and scaling must accompany each description. A nominally local control is not useful if its crosstalk grows unacceptably with system size. A universal construction based on exponentially large interaction strengths or exponentially precise controls is not a scalable simulator.

Dimensionless coordinates matter

Phase structure and dynamics usually depend on ratios such as

$$\frac{U}{J}, \qquad \frac{h}{J}, \qquad \frac{k_{\mathrm B}T}{J}, \qquad \frac{\gamma}{J/\hbar}, \qquad \frac{tJ}{\hbar},$$

not on a coupling in hertz by itself. Reporting dimensionless coordinates exposes whether a parameter scan changes only the intended target coordinate. For example, increasing lattice depth can simultaneously reduce tunnelling, change interaction matrix elements, alter heating, and modify adiabatic loading. The calibrated path through parameter space is then curved rather than a one-parameter line.

State Preparation

Analog evolution begins from a physical state, not from a symbolic ket. Common preparation routes include:

Route	Intended state	Principal difficulty
product initialization and quench	simple nonequilibrium input	SPAM bias and uncertain quench origin
adiabatic or quasi-adiabatic ramp	ground state connected to an easy state	small gaps, diabatic defects, heating
cooling and equilibration	thermal or low-energy ensemble	thermometry and equilibration evidence
dissipative preparation	steady state or dark state	generator identification and uniqueness
postselection	constrained or fixed-number sector	acceptance bias and resource accounting
local injection	excitation or response probe	pulse shape, spectral selectivity, back-action

For an interpolation $H(s)$ traversed at constant schedule speed $s=\tau/T$ over total laboratory time $T$, a standard adiabatic diagnostic compares matrix elements of $\partial_s H$ with squared gaps:

$$\max_{\substack{s\in[0,1]\\m\ne 0}} \frac{ \hbar \left| \langle m(s)| \partial_s H(s) |0(s)\rangle \right| }{ T \left| E_m(s)-E_0(s) \right|^2 } \ll 1.$$

This is a useful scale estimate, not a universal certificate. Degeneracies, many-body gap closing, nonsmooth schedules, finite temperature, and open dynamics require more careful statements. Experimental evidence may combine forward-and-reverse ramps, defect scaling, energy estimates, local thermometry, symmetry checks, and agreement with solvable sizes.

A ground-state simulator does not become a finite-temperature simulator by failing to cool. Likewise, a finite-temperature experiment is not defective when the target question explicitly concerns a Gibbs ensemble and the temperature is estimated with uncertainty.

Evolution Protocols

Analog simulation includes more than passive evolution under one static Hamiltonian.

Quenches

A sudden parameter change prepares a simple state under one Hamiltonian and evolves it under another. Quenches probe transport, correlation spreading, prethermalization, relaxation, confinement, and dynamical phase structure. The switching time must be short relative to target dynamics but slow enough to avoid unmodeled levels when those levels are outside the target.

Ramps and anneals

Slow schedules aim to follow instantaneous low-energy states. Their outcome depends on the minimum gap, matrix elements, noise spectrum, schedule shape, and initial ensemble. An observed low energy does not by itself establish ground-state fidelity, especially when many states have similar energy.

Adiabatic Quantum Computation explains the ideal closed-system path and its relation to circuits. Quantum Annealing treats driver–problem schedules, open dynamics, and endpoint samples. Applying either description to an experiment also requires the Hamiltonian and state mapping developed here.

Periodic driving

Floquet protocols use periodic controls to generate an effective stroboscopic Hamiltonian. The drive frequency, micromotion, prethermal window, and heating rate are part of the simulation specification. A Floquet Hamiltonian is not a static Hamiltonian valid at arbitrary observation times.

Engineered dissipation

Reservoirs can cool, stabilize constraints, prepare steady states, or realize an open-system target. Uncontrolled decoherence becomes a target process only when its operators, rates, correlations, and domain match the desired generator. Lindblad–GKSL Equation explains the Markovian formalism and its assumptions.

How Hamiltonian Error Reaches an Observable

Let

$$H_{\mathrm{act}} = H_{\mathrm{tar}} + \Delta H$$

after encoding and unit conversion, with finite-dimensional Hermitian operators (or suitable common-domain and bounded-perturbation assumptions), and write

$$U(t)=e^{-iH_{\mathrm{tar}}t/\hbar}, \qquad \widetilde U(t)=e^{-iH_{\mathrm{act}}t/\hbar}.$$

Duhamel’s identity gives

$$\widetilde U(t)-U(t) = -\frac{i}{\hbar} \int_0^t \widetilde U(t-s) \Delta H U(s) \, ds.$$

Unitary invariance and submultiplicativity of the operator norm imply

$$\left\| \widetilde U(t)-U(t) \right\| \le \frac{|t|}{\hbar} \left\| \Delta H \right\|.$$

For an initial density operator $\rho$ and bounded observable $O$,

$$\begin{aligned} \left| \widetilde\mu_O(t)-\mu_O(t) \right| &= \left| \operatorname{Tr} \left[ O \left( \widetilde U\rho\widetilde U^\dagger - U\rho U^\dagger \right) \right] \right| \\ &\le 2 \left\|O\right\|_\infty \left\| \widetilde U(t)-U(t) \right\| \\ &\le \frac{2|t|}{\hbar} \left\|O\right\|_\infty \left\|\Delta H\right\|. \end{aligned}$$

This bound explains why a small generator mismatch can accumulate over time. It is also often pessimistic. For an extensive many-body Hamiltonian, $\lVert\Delta H\rVert$ can grow with system size even when the observable is local. Locality, Lieb–Robinson bounds, perturbative response, measured susceptibilities, or direct short-time tests can give sharper observable-specific estimates. The Lieb–Robinson Bound supplies the relevant locality structure.

Calibration uncertainty

Suppose a reported observable is modeled as $\mu(\boldsymbol\theta)$ and calibrated parameters have covariance matrix $\Sigma_\theta$. Linear uncertainty propagation gives

$$\operatorname{Var}_{\mathrm{cal}}(\mu) \approx \nabla_{\boldsymbol\theta}\mu^{\mathsf T} \Sigma_\theta \nabla_{\boldsymbol\theta}\mu.$$

Correlations matter. If one laser intensity shifts two couplings together, adding independent parameter error bars in quadrature is wrong. Near a critical point or resonance, the linear approximation may fail because susceptibilities are large or the response is nonlinear. Parameter ensembles, bootstrap propagation, or posterior predictive distributions are then preferable.

Leakage is not an in-subspace perturbation

Define the leakage probability at time $t$ by

$$p_{\mathrm{leak}}(t) = \operatorname{Tr} \left[ Q\widetilde\rho(t) \right].$$

Postselecting the $P$ sector changes the estimand. The conditional state

$$\rho_P(t) = \frac{ P\widetilde\rho(t)P }{ 1-p_{\mathrm{leak}}(t) }$$

is defined when $p_{\mathrm{leak}}<1$ and assumes an ideal projective sector measurement. It can have accurate conditional observables even while the unconditional experiment has substantial failure probability. Both the conditional result and acceptance rate should be reported, along with whether leakage is detected without disturbing retained runs.

Error Classification

The same physical effect can be signal or error depending on the target.

Effect	Target feature when	Simulation error when
finite temperature	the target specifies a thermal ensemble	a ground-state claim silently uses a warm state
disorder	the disorder realization or ensemble is prescribed	uncontrolled inhomogeneity changes the intended clean model
dissipation	jump operators and rates define the target	decoherence is unmodeled or incorrectly correlated
long-range interaction	the target includes the calibrated tail	a nearest-neighbor model is claimed without bounding the tail
finite size and boundary	the finite geometry is the target	thermodynamic conclusions ignore size and edge effects
postselection	the conditional task is declared	discarded runs and resource cost are hidden
Floquet micromotion	observables are defined at physical times	only a stroboscopic effective model is reported as exact

This classification should be made at the level of the scientific question, not by labeling some hardware effects “good” and others “bad.” A noisy device does not automatically simulate open physics, and a long-range device does not automatically simulate a nearest-neighbor model.

Measurement and Inference

Analog simulators return detector records. The inference chain may include state-dependent fluorescence, atom or photon loss, parity projection, finite-resolution imaging, basis rotations, thresholding, and symmetry postselection.

For a discrete readout with ideal probabilities $\boldsymbol p$ and calibrated confusion matrix $M$, the observed distribution is

$$\boldsymbol q = M\boldsymbol p.$$

Inverting $M$ can remove bias under the model, but it amplifies variance when $M$ is ill-conditioned and can produce unphysical estimates at finite sample size. Regularization or constrained likelihood methods introduce their own bias. The raw counts, calibration data, estimator, and propagated uncertainty are needed to reproduce the uncertainty calculation.

The useful output is usually modest:

• local densities or spin populations;
• one- and two-point correlation functions;
• structure factors and susceptibilities;
• full counting statistics in an accessible basis;
• response functions after a calibrated perturbation;
• subsystem entropies or overlaps from randomized measurements;
• energy estimates from a known decomposition;
• phase-boundary or scaling estimates with finite-size qualifications.

State tomography is neither necessary nor generally scalable. Measuring more observables does not automatically produce stronger evidence if they are all derived from the same untested calibration model.

Worked Example: A Rydberg Array as an Ising Simulator

Consider two internal states per atom: a ground state $|g\rangle$ and a Rydberg state $|r\rangle$. Choose Pauli conventions

$$Z_i|r_i\rangle=|r_i\rangle, \qquad Z_i|g_i\rangle=-|g_i\rangle, \qquad n_i=|r_i\rangle\langle r_i| = \frac{I+Z_i}{2}.$$

In a rotating-frame and rotating-wave description, a standard laboratory Hamiltonian is

$$H_{\mathrm{Ry}} = \sum_i \left( \frac{\hbar\Omega_i}{2}X_i - \hbar\Delta_i n_i \right) + \sum_{i<j} V_{ij}n_in_j.$$

For van der Waals interactions away from resonances,

$$V_{ij} \approx \frac{C_6}{R_{ij}^6},$$

with a sign and angular dependence set by the chosen states and geometry. The platform-specific interaction calculation and blockade approximations belong on Rydberg Atoms and Rydberg Blockade.

Exact spin expansion within the two-level model

Using

$$n_in_j = \frac{1}{4} \left( I+Z_i+Z_j+Z_iZ_j \right),$$

the Hamiltonian becomes

$$H_{\mathrm{Ry}} = C I + \sum_i \left( \frac{\hbar\Omega_i}{2}X_i + h_i^z Z_i \right) + \sum_{i<j} \frac{V_{ij}}{4}Z_iZ_j,$$

where

$$C = -\frac{\hbar}{2} \sum_i\Delta_i + \frac{1}{4} \sum_{i<j}V_{ij},$$

and

$$h_i^z = -\frac{\hbar\Delta_i}{2} + \frac{1}{4} \sum_{j\ne i}V_{ij}.$$

The interaction therefore produces both Ising couplings and site-dependent longitudinal fields. For a desired zero longitudinal field, the detuning must satisfy

$$\hbar\Delta_i = \frac{1}{2} \sum_{j\ne i}V_{ij}.$$

A single global detuning cancels this field only when the interaction sum is site independent. Edges, vacancies, geometry errors, and anisotropic interactions spoil that simplification. This is a direct example of why the control-to-model map must include geometry and the realized occupancy pattern.

Long-range tails define the actual target

Suppose the intended interpretation retains only a nearest-neighbor interaction $V_1$. The omitted coherent term is

$$\Delta H_{\mathrm{tail}} = \frac{1}{4} \sum_{\substack{i<j\\|i-j|>1}} V_{ij}Z_iZ_j + \text{associated longitudinal fields}.$$

A worst-case bound is

$$\left\| \Delta H_{\mathrm{tail}} \right\| \le \frac{1}{4} \sum_{\substack{i<j\\|i-j|>1}} \left|V_{ij}\right| + \left\| \Delta H_z \right\|.$$

For an infinite one-dimensional chain with $V_r=V_1/r^6$, the sum of pair-coupling coefficients per site beyond nearest neighbors, divided by its nearest-neighbor value, is

$$\sum_{r=2}^{\infty}\frac{1}{r^6} = \zeta(6)-1 \approx 0.01734.$$

This is a coefficient ratio, not an energy ratio for an arbitrary state: the actual energy weights each pair by $\langle n_i n_j\rangle$ (or by the corresponding spin correlation). The ratio is small but not universally negligible. It can shift phase boundaries, lift degeneracies, alter constrained dynamics, or matter at long times. The scientifically clean alternatives are to include the calibrated tail in the target, demonstrate that the chosen observable is insensitive to it, or carry its effect as model discrepancy.

The blockade limit is another effective model

When nearest-neighbor excitation pairs cost an energy much larger than the drive and detuning scales, the dynamics can be projected into the no-adjacent- excitation subspace. The leading constrained Hamiltonian has the form

$$H_{\mathrm{PXP}} \approx \frac{\hbar\Omega}{2} \sum_i P_{i-1}X_iP_{i+1} - \hbar\Delta \sum_i n_i,$$

where $P_i=I-n_i$ projects site $i$ onto $|g_i\rangle$; at an open edge the missing neighbor projector is the identity. This formula also assumes longer-range interactions are negligible on the time scale of interest. They survive the nearest-neighbor blockade projection as $\sum_{|i-j|>1}V_{ij}n_in_j$. For example, $V_1/(\hbar\Omega)=100$ implies $V_2/(\hbar\Omega)=100/64\simeq1.56$ in a van der Waals chain: good nearest-neighbor blockade alone does not give the pure PXP model. Finite blockade, decay, Doppler shifts, and preparation defects introduce further corrections. “PXP simulator” therefore names an effective regime, not an exact identity of the microscopic apparatus.

Platform Families

Analog simulation is a method, not a single hardware architecture.

Platform	Native degrees of freedom and controls	Common target families	Characteristic model risk
optical lattices and quantum gases	atomic motion, internal states, lattice depth, scattering length, geometry	Bose–Hubbard, Fermi–Hubbard, spin, gauge, and continuum models	multiband terms, trap envelope, thermometry, loss, heating
neutral-atom and Rydberg arrays	two-level atoms, geometry, drive, detuning, pair interactions	Ising, XY, constrained, dimer, and optimization models	interaction tails, vacancies, decay, Doppler shifts, edge fields
trapped ions	internal spins, collective phonons, optical forces, mode shaping	long-range Ising and XY models, spin–boson dynamics	residual spin–motion entanglement, mode drift, interaction-matrix mismatch
superconducting circuits	qubits, resonators, tunable couplers, microwave drives	spin lattices, Bose–Hubbard variants, synthetic gauge and driven models	frequency disorder, crosstalk, photon loss, higher levels
polar molecules and magnetic atoms	rotational or spin states, dipolar interactions, external fields	long-range spin and lattice models	loss, tensor interactions, state-dependent trapping
photonic and polaritonic arrays	modes, hopping, drive, dissipation, nonlinearities	bosonic transport, synthetic dimensions, driven–dissipative matter	preparation and loss are inseparable from steady-state interpretation

Optical Lattices, Trapped-Ion Control, Neutral-Atom and Rydberg Qubits, and Superconducting Qubits give the platform-specific control and error physics. Lattice Models Overview and its linked model pages define the target Hamiltonians and observables.

Verification and Validation

Verification asks whether the apparatus implements the declared target task. Validation asks whether that target model answers the intended physical question. Verification of Quantum Simulation compares these protocols in depth; the essential analog workflow is given here.

Independent tests of the correspondence

1. Component calibration: characterize controls, couplings, detector response, loss, and drift independently where possible.
2. Generator checks: infer or constrain effective Hamiltonian and dissipative terms from dynamics, spectroscopy, response, or conservation laws.
3. Solvable overlap: compare with exact diagonalization, free limits, tensor networks, Monte Carlo without a sign problem, or perturbation theory where those methods are reliable.
4. Convergence and sensitivity: vary size, boundary, preparation, evolution time, interaction cutoff, and analysis assumptions.
5. Held-out prediction: calibrate on one set of states, times, or observables and predict another set without refitting.
6. Independent reproduction: compare platforms, implementations, observables, or laboratories with distinct systematic errors.
7. Frontier inference: report only the observable and regime supported by the preceding layers, with uncertainty and classical baselines.

Agreement in the calibration set is not yet prediction. A sufficiently flexible model can absorb error into fitted parameters. Held-out times, observables, initial states, or geometries test whether the inferred generator transfers beyond the data that selected it.

Hamiltonian learning

For a structured ansatz

$$H(\boldsymbol\theta) = \sum_{a=1}^{m}\theta_a H_a,$$

Hamiltonian learning estimates $\boldsymbol\theta$ from measured dynamics, steady states, response, or conserved quantities. It can expose coherent terms that component calibration misses. Its conclusion is conditional on the ansatz, identifiability, data range, SPAM model, and statistical uncertainty. A small residual within a misspecified ansatz does not prove that unmodeled operators are absent.

Generator-learning methods can propagate uncertainty in inferred Hamiltonian and dissipative parameters to observable bounds; Kraft and colleagues provide one construction. Those bounds remain conditional on the model class, measurement assumptions, and tested time range. Analog Quantum Simulation research directions discusses open problems in extending such methods.

Cross-platform comparison

Independent simulators can compare local observables, probability distributions, or state overlaps inferred from randomized measurements. Agreement becomes stronger evidence when the platforms have different microscopic errors. It still does not prove correctness if both implement the same mistaken target map or share the same analysis assumptions.

Why extrapolation is difficult

Classical comparison often stops where the analog simulator becomes most interesting. Trust in that frontier regime must be earned by overlap: calibrated generators, size and time trends, conservation laws, solvable limits, alternative observables, and reproducible data. A single qualitative match at small size does not validate arbitrary larger systems or later times.

Experimental time and sampling cost

Continuous evolution is only part of a run. Cooling, loading, rearrangement, calibration, postselection, and repeated destructive measurement can dominate wall-clock time. The important costs therefore include occupied sites or modes, independent controls, preparation and evolution time, accepted fraction, measurement repetitions, and classical fitting. Geometry design, waveform optimization, and parameter inversion are also real computational work.

Analog, Digital, and Hybrid Boundaries

Question	Analog simulation	Digital simulation
What is programmed?	a physical Hamiltonian, channel, geometry, or control schedule	a logical circuit or algorithm
How does evolution occur?	native continuous many-body dynamics	discrete compiled operations
Typical approximation knob	detuning, scale separation, protection, truncation, or control precision	step count, polynomial degree, synthesis precision, or code distance
Flexibility	restricted but often physically natural family	broad under a universal gate and access model
Typical dominant risk	target–device mismatch and uncorrected physical noise	algorithmic overhead plus gate and logical error
Verification focus	generator, preparation, observable map, and held-out behavior	circuit implementation, logical error, and algorithmic approximation

The boundary is functional rather than cosmetic. Piecewise control does not make a protocol digital if each segment uses an uncompiled many-body interaction as its primitive. Conversely, a device built from analog pulses can implement digital gates when those pulses realize a calibrated discrete instruction set.

Digital–analog protocols compose native interaction blocks with discrete operations. Variational protocols use measured quantum states inside a classical feedback loop. Hybrid Quantum Simulation explains block composition, projected dynamics, self-consistent embedding, and feedback stability. What Is Quantum Simulation? compares the broader simulation approaches.

Claims About Quantum Advantage

Analog simulation can be scientifically useful before any complexity advantage is established. It can realize clean model experiments, reveal unexpected observables, test approximations, and guide theory in regimes where classical methods remain available.

A claim of advantage requires a sharper comparison:

1. define the input family and output quantity;
2. state accuracy and confidence;
3. count preparation, calibration, rejected runs, sampling, and classical inference;
4. compare with the best applicable classical methods for the same task;
5. justify correctness where full classical calculation is unavailable;
6. distinguish a hardware-scale milestone from a useful scientific result.

An exponentially large Hilbert space, a Monte Carlo sign problem, or failure of one tensor-network calculation is not by itself a proof. Classical tractability depends on model structure, dimension, entanglement, observable, time, and requested precision. Verification of Quantum Advantage explains what different advantage claims establish.

Exercises

1. Time rescaling and the constant energy shift

Assume

$$H_{\mathrm{lab}}V = \alpha V H_{\mathrm{tar}} + \beta V.$$

Show directly from the power series that

$$e^{-iH_{\mathrm{lab}}\tau/\hbar}V = e^{-i\beta\tau/\hbar} V e^{-iH_{\mathrm{tar}}\alpha\tau/\hbar}.$$

Explain why $\beta$ is irrelevant for density-operator dynamics but can matter if the encoded sector interferes with a reference sector having a different energy shift.

Solution

The assumed intertwining relation implies by induction that

$$H_{\mathrm{lab}}^nV = V \left( \alpha H_{\mathrm{tar}}+\beta I \right)^n.$$

Substituting into the exponential series gives

$$\begin{aligned} e^{-iH_{\mathrm{lab}}\tau/\hbar}V &= V e^{-i(\alpha H_{\mathrm{tar}}+\beta I)\tau/\hbar} \\ &= e^{-i\beta\tau/\hbar} V e^{-iH_{\mathrm{tar}}\alpha\tau/\hbar}. \end{aligned}$$

For an encoded density operator, the phase multiplies the ket side and its complex conjugate multiplies the bra side, so it cancels. If a coherent reference occupies another sector with shift $\beta'$, the relative phase $(\beta-\beta')\tau/\hbar$ is observable and cannot be discarded.

2. Derive the generator-error bound

Starting from

$$F(s) = \widetilde U(t-s)U(s),$$

differentiate with respect to $s$, integrate from $0$ to $t$, and derive Duhamel’s identity and

$$\left\| \widetilde U(t)-U(t) \right\| \le \frac{|t|}{\hbar} \left\|\Delta H\right\|.$$

Solution

With $i\hbar\,\partial_tU=H_{\mathrm{tar}}U$ and $i\hbar\,\partial_t\widetilde U=(H_{\mathrm{tar}}+\Delta H)\widetilde U$, the product rule gives

$$\frac{dF}{ds} = \frac{i}{\hbar} \widetilde U(t-s) \Delta H U(s).$$

Since $F(t)=U(t)$ and $F(0)=\widetilde U(t)$,

$$U(t)-\widetilde U(t) = \frac{i}{\hbar} \int_0^t \widetilde U(t-s) \Delta H U(s)\,ds.$$

Rearranging gives the stated Duhamel identity. Both evolution operators have operator norm one, so

$$\begin{aligned} \left\| \widetilde U(t)-U(t) \right\| &\le \frac{1}{\hbar} \int_0^{|t|} \left\|\Delta H\right\|\,ds \\ &= \frac{|t|}{\hbar} \left\|\Delta H\right\|. \end{aligned}$$

3. Rydberg-to-Ising mapping

Starting from

$$H = \sum_i \left( \frac{\hbar\Omega_i}{2}X_i - \hbar\Delta_i n_i \right) + \sum_{i<j}V_{ij}n_in_j,$$

with $n_i=(I+Z_i)/2$, derive the constant $C$, longitudinal fields $h_i^z$, and Ising couplings. For which detunings do the longitudinal fields vanish?

Solution

Expanding each number operator gives

$$-\hbar\Delta_i n_i = -\frac{\hbar\Delta_i}{2}I - \frac{\hbar\Delta_i}{2}Z_i$$

and

$$V_{ij}n_in_j = \frac{V_{ij}}{4} \left( I+Z_i+Z_j+Z_iZ_j \right).$$

Collecting terms,

$$H = C I + \sum_i \left( \frac{\hbar\Omega_i}{2}X_i+h_i^zZ_i \right) + \sum_{i<j}\frac{V_{ij}}{4}Z_iZ_j,$$

with

$$C = -\frac{\hbar}{2}\sum_i\Delta_i + \frac{1}{4}\sum_{i<j}V_{ij},$$

and

$$h_i^z = -\frac{\hbar\Delta_i}{2} + \frac{1}{4}\sum_{j\ne i}V_{ij}.$$

Thus $h_i^z=0$ when

$$\hbar\Delta_i = \frac{1}{2}\sum_{j\ne i}V_{ij}.$$

A global detuning suffices only if the interaction sum is the same at every site.

4. Quantify a van der Waals tail

For an infinite one-dimensional array with spacing $a$ and $V_r=V_1/r^6$, compute the pair-coupling sum per site from all distances $r\ge 2$ relative to the nearest-neighbor coefficient. State when it also equals the ratio of interaction energies. Explain why the result does not alone establish observable-level accuracy.

Solution

Counting each pair once, the coupling sum per site is $\sum_{r\ge1}V_r$. It equals the interaction energy per site in the fully excited Rydberg product state, for which every $\langle n_i n_j\rangle=1$. Therefore

$$\frac{ \sum_{r=2}^{\infty}V_r }{ V_1 } = \sum_{r=2}^{\infty}\frac{1}{r^6} = \zeta(6)-1 \approx 0.017343.$$

This is about $1.73\%$ of the nearest-neighbor coupling sum. In a blockaded state the nearest-neighbor pair occupation may vanish while farther-neighbor pair occupations remain finite, so it need not be a $1.73\%$ correction to that state’s interaction energy. It is not a universal observable-error bound. Long-time phases can accumulate, near-degenerate states can be sensitive to weak terms, and some observables respond much more strongly than the energy density. The tail should be included in the modeled generator or propagated to the particular observable.

5. Correlated calibration error

Let

$$\mu(J,h) = \frac{h}{J},$$

and suppose the calibrated pair $(J,h)$ has standard deviations $\sigma_J,\sigma_h$ and correlation coefficient $\rho$. Use linear propagation to find $\operatorname{Var}(\mu)$.

Solution

The gradient is

$$\nabla\mu = \left( -\frac{h}{J^2}, \frac{1}{J} \right).$$

With covariance $\operatorname{Cov}(J,h)=\rho\sigma_J\sigma_h$,

$$\operatorname{Var}(\mu) \approx \frac{h^2}{J^4}\sigma_J^2 + \frac{1}{J^2}\sigma_h^2 - \frac{2h}{J^3} \rho\sigma_J\sigma_h.$$

Positive correlation can reduce uncertainty in the ratio because common-mode scale fluctuations partly cancel. Treating $J$ and $h$ as independent would miss this effect.

6. Readout inversion

A binary detector reports the wrong state with symmetric probability $e$. Write its confusion matrix, invert the relation between the observed magnetization $\langle Z\rangle_{\mathrm{obs}}$ and the ideal magnetization, and identify the variance-amplification factor.

Solution

In the ordered basis of the two outcomes,

$$M = \begin{pmatrix} 1-e & e\\ e & 1-e \end{pmatrix}.$$

Symmetric flips reverse the sign of $Z$ with probability $e$, so

$$\langle Z\rangle_{\mathrm{obs}} = (1-2e) \langle Z\rangle_{\mathrm{ideal}}.$$

For $e\ne 1/2$,

$$\langle Z\rangle_{\mathrm{ideal}} = \frac{ \langle Z\rangle_{\mathrm{obs}} }{ 1-2e }.$$

Consequently, the sampling variance of the corrected estimator is multiplied by $1/(1-2e)^2$, before including uncertainty in the calibration of $e$. Bias correction therefore does not restore the information lost by the detector.

7. Adiabatic bottleneck

Suppose the smallest relevant gap along a ramp scales as $\Delta_{\min}(L)\propto L^{-z}$ while the largest relevant matrix element of $\partial_sH$ scales as $L^\kappa$. Using the stated adiabatic diagnostic, find the scaling required of the ramp duration $T$.

Solution

The diagnostic requires

$$\frac{ \hbar L^\kappa }{ T\Delta_{\min}^2 } \ll 1.$$

Since $\Delta_{\min}^2\propto L^{-2z}$, the ratio scales as $L^{\kappa+2z}/T$. Thus $T=C L^{\kappa+2z}$ with a sufficiently large constant $C$ keeps this diagnostic small at a fixed target tolerance. Making the diagnostic vanish as $L\to\infty$ instead requires $T/L^{\kappa+2z}\to\infty$. The result is only a scaling diagnostic: optimized schedules, selection rules, degeneracies, open-system relaxation, and the desired local rather than global accuracy can change the requirement.

8. Small leakage, appreciable phase error

A retained state $|0\rangle$ couples with real energy $g$ to a discarded state $|\ell\rangle$ of energy $\Delta>0$:

$$H=\begin{pmatrix}0&g\\g&\Delta\end{pmatrix}, \qquad |g|\ll\Delta.$$

Find the lower eigenvalue and the maximum discarded population from initial $|0\rangle$. If another retained state $|1\rangle$ remains at zero energy and uncoupled, why can neglecting the discarded level still give a large relative-phase error?

Solution

Writing $\omega_E=\sqrt{\Delta^2+4g^2}$, the lower eigenvalue is

$$E_-=\frac{\Delta-\omega_E}{2} =-\frac{g^2}{\Delta}+O(g^4/\Delta^3).$$

Direct two-level evolution gives

$$p_\ell(t)=\frac{4g^2}{\Delta^2+4g^2} \sin^2\left(\frac{\omega_E t}{2\hbar}\right).$$

The maximum population is only $O(g^2/\Delta^2)$, but the predominantly retained eigenstate acquires phase $-E_-t/\hbar$. Relative to $|1\rangle$, this is order one at times $t\sim\hbar\Delta/g^2$. A superposition of the retained levels therefore reveals a large coherent error even though leakage remains small. The effective Hamiltonian must include the virtual energy shift, together with the small dressing of states and observables.

References

• E. Bairey, I. Arad, and N. H. Lindner, “Learning a local Hamiltonian from local measurements,” Physical Review Letters 122, 020504 (2019), doi:10.1103/PhysRevLett.122.020504.
• H. Bernien et al., “Probing many-body dynamics on a 51-atom quantum simulator,” Nature 551, 579–584 (2017), doi:10.1038/nature24622.
• I. Bloch, J. Dalibard, and S. Nascimbène, “Quantum simulations with ultracold quantum gases,” Nature Physics 8, 267–276 (2012), doi:10.1038/nphys2259.
• S. Bravyi, D. P. DiVincenzo, and D. Loss, “Schrieffer–Wolff transformation for quantum many-body systems,” Annals of Physics 326, 2793–2826 (2011), doi:10.1016/j.aop.2011.06.004.
• A. Browaeys and T. Lahaye, “Many-body physics with individually controlled Rydberg atoms,” Nature Physics 16, 132–142 (2020), doi:10.1038/s41567-019-0733-z.
• T. Cubitt, A. Montanaro, and S. Piddock, “Universal quantum Hamiltonians,” Proceedings of the National Academy of Sciences 115, 9497–9502 (2018), doi:10.1073/pnas.1804949115.
• A. J. Daley et al., “Practical quantum advantage in quantum simulation,” Nature 607, 667–676 (2022), doi:10.1038/s41586-022-04940-6.
• A. Elben et al., “Cross-platform verification of intermediate scale quantum devices,” Physical Review Letters 124, 010504 (2020), doi:10.1103/PhysRevLett.124.010504.
• P. Hauke, F. M. Cucchietti, L. Tagliacozzo, I. Deutsch, and M. Lewenstein, “Can one trust quantum simulators?” Reports on Progress in Physics 75, 082401 (2012), doi:10.1088/0034-4885/75/8/082401.
• M. Kliesch and I. Roth, “Theory of quantum system certification,” PRX Quantum 2, 010201 (2021), doi:10.1103/PRXQuantum.2.010201.
• T. Kraft et al., “Bounded-error quantum simulation via Hamiltonian and Lindbladian learning,” Physical Review X 16, 031037 (2026), doi:10.1103/s96t-n8tx.
• C. Monroe et al., “Programmable quantum simulations of spin systems with trapped ions,” Reviews of Modern Physics 93, 025001 (2021), doi:10.1103/RevModPhys.93.025001.