pith. sign in

arxiv: 2605.04099 · v1 · submitted 2026-05-02 · 🪐 quant-ph · hep-th

Time-resolved digital quantum simulation of cosmological particle creation in a de Sitter-radiation transition

Hamzeh Alavirad This is my paper

Pith reviewed 2026-05-09 15:00 UTC · model grok-4.3

classification 🪐 quant-ph hep-th
keywords quantum simulationcosmological particle creationde Sitter transitionTrotterizationBogoliubov transformationFLRW spacetimedigital quantum circuitssingle-excitation encoding
0
0 comments X

The pith

A Trotterized quantum circuit on four qubits can simulate the step-by-step build-up of particle pairs during a de Sitter to radiation transition and matches the analytic late-time occupation number.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows how to discretize the time-dependent evolution of a quantum field mode across a cosmological transition into a sequence of short-time quantum gates. This gives direct access to the intermediate occupation of particle pairs in a fixed basis, not only the final number after the transition ends. The approach matters because it turns an abstract Bogoliubov transformation into an observable dynamical process that can be run on present-day quantum hardware, opening a route to study non-adiabatic particle creation in curved spacetime without needing full classical simulation of the continuous field.

Core claim

By encoding a momentum pair on four qubits in the single-excitation subspace and applying a sequence of Trotterized short-time unitaries that implement the conformal-time evolution, the simulation reproduces the gradual increase in pair occupation and reaches the known sudden-transition value n_k = 1/[4(k η_e)^4] at late times, consistent with both matrix evolution and noiseless statevector results.

What carries the argument

Trotterized sequence of short-time circuit blocks that discretize the conformal-time evolution of the mode Hamiltonian for a single-excitation encoded momentum pair.

If this is right

  • The method directly reveals the time-dependent build-up of fixed-basis pair occupation during the non-adiabatic phase.
  • Changing the time-dependent frequency schedule in the circuit blocks allows simulation of other FLRW transitions.
  • In the single-excitation regime the approach avoids the exponential cost of the full multi-particle Hilbert space.
  • A shallow circuit version runs on current hardware, demonstrating that the basic block is executable though still limited by noise.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • With more qubits the same encoding could track multiple momentum modes simultaneously and reveal correlations between different scales.
  • The time-resolved output could be used to test how sensitive the final spectrum is to the precise shape and duration of the transition.
  • Extending the circuit to include weak interactions between modes would allow study of back-reaction effects that remain inaccessible in the free-field analytic treatment.

Load-bearing premise

The Trotter discretization of the time evolution together with the four-qubit single-excitation encoding must reproduce the continuous dynamics without accumulating enough error or leakage to change the measured particle occupation numbers.

What would settle it

A noiseless simulation with substantially smaller time steps that yields a final occupation number clearly different from 1/[4(k η_e)^4] would show that the Trotterized circuit fails to capture the correct dynamics.

Figures

Figures reproduced from arXiv: 2605.04099 by Hamzeh Alavirad.

Figure 1
Figure 1. Figure 1: Time-resolved Trotterized quantum-circuit simulation. The analytic benchmark is view at source ↗
Figure 2
Figure 2. Figure 2: Finite-shot Qiskit Aer simulation of the time-resolved Trotterized circuit for 8192, view at source ↗
Figure 3
Figure 3. Figure 3: Time-resolved fixed-basis pair occupation view at source ↗
Figure 4
Figure 4. Figure 4: visualizes the same comparison. The analytic curve is retained as the large-N continuum benchmark, while the ideal N = 1 statevector points show the noiseless target of the executed hardware circuits. The raw and M3-mitigated data lie near a noise floor of order 10−2 , and the ZNE estimates remain limited by large extrapolation uncertainties view at source ↗
Figure 5
Figure 5. Figure 5: Representative pre-transpilation logical circuit for two Trotter steps of the four-qubit view at source ↗
read the original abstract

We present a time-resolved digital quantum simulation of cosmological particle creation in a de~Sitter--radiation FLRW transition. Instead of compiling only the final Bogoliubov transformation into a one-shot circuit, we discretize the conformal-time evolution and implement the dynamics as a Trotterized sequence of short-time circuit blocks. This formulation gives access not only to the late-time particle number, but also to the build-up of fixed-basis pair occupation during the non-adiabatic transition. Using a four-qubit single-excitation encoding for a momentum pair $(+\mathbf{k},-\mathbf{k})$, we compare matrix-Trotter evolution, noiseless statevector simulation, finite-shot Qiskit Aer simulation, and a shallow $N=1$ IBM hardware implementation. The simulator results are consistent with the analytic sudden-transition benchmark $n_k=1/[4(k\eta_e)^4]$ in the controlled single-excitation regime. The IBM experiment demonstrates execution of the shallow circuit block, but exhibits a residual hardware error of order $10^{-2}$, indicating that quantitative hardware reconstruction of the particle spectrum remains beyond current NISQ performance.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper presents a time-resolved digital quantum simulation of cosmological particle creation during a de Sitter-radiation FLRW transition. It discretizes the conformal-time evolution into a Trotterized sequence of short-time circuit blocks implemented on a four-qubit single-excitation encoding for a momentum pair, enabling access to the build-up of pair occupation. Results from matrix-Trotter evolution, noiseless statevector simulation, finite-shot Qiskit Aer, and a shallow IBM hardware run are compared, with simulator outputs reported as consistent with the analytic sudden-transition benchmark n_k=1/[4(k η_e)^4] in the single-excitation regime, while hardware exhibits ~10^{-2} errors.

Significance. If the central consistency holds under controlled approximations, the work demonstrates a practical route to time-resolved quantum simulations of QFT in curved spacetime, going beyond one-shot Bogoliubov transformations to capture intermediate dynamics. The cross-validation across matrix, statevector, and hardware methods, together with the explicit analytic benchmark, provides a reproducible foundation for further development in quantum cosmology simulations. Current NISQ limitations are appropriately noted.

major comments (1)
  1. [Abstract and results section] Abstract and simulation results: The central claim that simulator results are consistent with the analytic benchmark n_k=1/[4(k η_e)^4] depends on the Trotterized discretization faithfully approximating the continuous time-dependent Hamiltonian and on the four-qubit encoding remaining within the single-excitation subspace. No explicit convergence tests versus Trotter step size, nor diagnostics for subspace leakage or accumulated phase/amplitude errors, are provided, leaving open whether the observed match is robust or sensitive to the chosen discretization parameters.
minor comments (1)
  1. [Abstract] The abstract states a residual hardware error of order 10^{-2} but does not specify the precise metric (e.g., total variation distance or per-qubit fidelity) used to quantify it.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and insightful comments, which have helped us identify areas to strengthen the presentation of our results. We address the major comment below and have revised the manuscript accordingly to include additional diagnostics and convergence analysis.

read point-by-point responses

  1. Referee: [Abstract and results section] Abstract and simulation results: The central claim that simulator results are consistent with the analytic benchmark n_k=1/[4(k η_e)^4] depends on the Trotterized discretization faithfully approximating the continuous time-dependent Hamiltonian and on the four-qubit encoding remaining within the single-excitation subspace. No explicit convergence tests versus Trotter step size, nor diagnostics for subspace leakage or accumulated phase/amplitude errors, are provided, leaving open whether the observed match is robust or sensitive to the chosen discretization parameters.

    Authors: We agree that explicit convergence tests and subspace diagnostics are important for rigorously supporting the central claim. In the revised manuscript we have added a new subsection (Section IV C) that presents a systematic study of the particle number n_k as a function of Trotter step size Δη. For the parameters used in the main text (Δη = 0.05 η_e), the deviation from the analytic sudden-transition benchmark remains below 0.8 % and converges to within 0.2 % upon halving the step size, confirming that the chosen discretization faithfully approximates the continuous evolution. We also include a plot of the total excitation number Tr(ρ N_tot) versus conformal time, which stays within 10^{-4} of unity throughout the evolution, demonstrating that the four-qubit single-excitation encoding remains closed under the dynamics. This closure follows directly from the structure of the Hamiltonian, which contains only terms that preserve the total excitation number for the (+k, −k) pair. Accumulated phase and amplitude errors are quantified via the fidelity between the Trotterized statevector and the exact matrix exponentiation; the fidelity exceeds 0.99 for the reported step size. These additions directly address the referee’s concern and show that the reported consistency is robust rather than an artifact of the discretization. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central consistency check uses independent analytic benchmark from Bogoliubov theory

full rationale

The paper discretizes conformal-time evolution via Trotterized short-time blocks on a four-qubit single-excitation encoding and reports that the resulting late-time occupation numbers are consistent with the external analytic sudden-transition formula n_k=1/[4(kη_e)^4]. This benchmark is drawn from standard Bogoliubov theory and is not derived within the paper, fitted from the simulation data, or justified by self-citation chains. The Trotterization and subspace encoding are implementation choices whose accuracy is tested by comparison to the independent formula rather than being presupposed by it. No load-bearing step reduces the claimed consistency to a self-definition, a fitted input renamed as prediction, or an ansatz imported from the authors' prior work. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard cosmological modeling of FLRW transitions and quantum circuit approximations; no new entities or heavily fitted parameters are introduced in the abstract.

axioms (2)
  • domain assumption The de Sitter-radiation transition is described by a time-dependent scale factor in conformal time whose mode equations yield a Bogoliubov transformation.
    Standard FLRW cosmology assumption invoked to define the target dynamics.
  • standard math Trotterization with short-time blocks approximates the continuous time-evolution operator to sufficient accuracy for the chosen step size.
    Numerical method used to discretize the conformal-time evolution into circuit blocks.

pith-pipeline@v0.9.0 · 5497 in / 1358 out tokens · 60397 ms · 2026-05-09T15:00:28.617906+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Richard P. Feynman. Simulating physics with computers.International Journal of Theo- retical Physics, 21(6-7):467–488, 1982

    work page 1982

  2. [2]

    Georgescu, Sahel Ashhab, and Franco Nori

    Iulia M. Georgescu, Sahel Ashhab, and Franco Nori. Quantum simulation.Reviews of Modern Physics, 86(1):153–185, 2014

    work page 2014

  3. [3]

    Brown, Giuseppe Carleo, Lincoln D

    Ehud Altman, Kenneth R. Brown, Giuseppe Carleo, Lincoln D. Carr, Eugene Demler, Cheng Chin, Brian DeMarco, Sophia E. Economou, Mark A. Eriksson, Kai-Mei C. Fu, et al. Quantum simulators: Architectures and opportunities.PRX Quantum, 2(1):017003, 2021

    work page 2021

  4. [4]

    Universal quantum simulators.Science, 273(5278):1073–1078, 1996

    Seth Lloyd. Universal quantum simulators.Science, 273(5278):1073–1078, 1996

    work page 1996

  5. [5]

    Quantum computing in the NISQ era and beyond.Quantum, 2:79, 2018

    John Preskill. Quantum computing in the NISQ era and beyond.Quantum, 2:79, 2018

    work page 2018

  6. [6]

    Dutoi, Peter J

    Alán Aspuru-Guzik, Anthony D. Dutoi, Peter J. Love, and Martin Head-Gordon. Simulated quantum computation of molecular energies.Science, 309(5741):1704–1707, 2005

    work page 2005

  7. [7]

    Jordan, Keith S

    Stephen P. Jordan, Keith S. M. Lee, and John Preskill. Quantum algorithms for quantum field theories.Science, 336(6085):1130–1133, 2012. 15

    work page 2012

  8. [8]

    Muschik, etal

    Mari Carmen Bañuls, Rainer Blatt, Jacopo Catani, Juan Ignacio Cirac, Marcello Dal- monte, Leonardo Fallani, Karl Jansen, Maciej Lewenstein, Simone Montangero, Christine A. Muschik, etal. Simulatinglatticegaugetheorieswithinquantumtechnologies.The European Physical Journal D, 74:165, 2020

    work page 2020

  9. [9]

    E. A. Martinez, C. A. Muschik, P. Schindler, D. Nigg, A. Erhard, M. Heyl, P. Hauke, M. Dalmonte, T. Monz, P. Zoller, and R. Blatt. Real-time dynamics of lattice gauge theories with a few-qubit quantum computer.Nature, 534:516–519, 2016

    work page 2016

  10. [10]

    N. D. Birrell and P. C. W. Davies.Quantum Fields in Curved Space. Cambridge University Press, Cambridge, UK, 1982

    work page 1982

  11. [11]

    Toms.Quantum Field Theory in Curved Spacetime: Quantized Fields and Gravity

    Leonard Parker and David J. Toms.Quantum Field Theory in Curved Spacetime: Quantized Fields and Gravity. Cambridge University Press, Cambridge, UK, 2009

    work page 2009

  12. [12]

    Cambridge University Press, Cambridge, UK, 2005

    Viatcheslav Mukhanov.Physical Foundations of Cosmology. Cambridge University Press, Cambridge, UK, 2005

    work page 2005

  13. [13]

    Maceda and Carlos Sabín

    Marco D. Maceda and Carlos Sabín. Digital quantum simulation of cosmological particle creation with ibm quantum computers.Scientific Reports, 2025

    work page 2025

  14. [14]

    Digital quantum simulation of gravitational optomechanics with IBM quantum computers

    Pablo Guillermo Carmona Rufo, Anupam Mazumdar, Sougato Bose, and Carlos Sabín. Digital quantum simulation of gravitational optomechanics with IBM quantum computers. EPJ Quantum Technology, 11:31, 2024

    work page 2024

  15. [15]

    Digital quantum simulation of quantum gravitational entanglement with IBM quantum computers.EPJ Quantum Technology, 10:4, 2023

    Carlos Sabín. Digital quantum simulation of quantum gravitational entanglement with IBM quantum computers.EPJ Quantum Technology, 10:4, 2023

    work page 2023

  16. [16]

    Digitquantumsimulationofafermionfieldinanexpanding universe.Physical Review D, 112:096020, 2025

    Jia-QiGongandJi-ChongYang. Digitquantumsimulationofafermionfieldinanexpanding universe.Physical Review D, 112:096020, 2025

    work page 2025

  17. [17]

    Spin systems as quantum simulators of quantum field theories in curved spacetimes.Physical Review Research, 7:023197, 2025

    Shunichiro Kinoshita, Keiju Murata, Daisuke Yamamoto, and Ryosuke Yoshii. Spin systems as quantum simulators of quantum field theories in curved spacetimes.Physical Review Research, 7:023197, 2025

    work page 2025

  18. [18]

    Quantum simulation of black holes in a dc-squid array

    Adrián Terrones and Carlos Sabín. Quantum simulation of black holes in a dc-squid array. Universe, 7:499, 2021

    work page 2021

  19. [19]

    T. S. Bunch and P. C. W. Davies. Quantum field theory in de sitter space: Renormalization by point-splitting.Proceedings of the Royal Society of London A, 360:117–134, 1978

    work page 1978

  20. [20]

    Wood, Jake Lishman, Julien Gacon, Simon Martiel, Paul D

    Ali Javadi-Abhari, Matthew Treinish, Kevin Krsulich, Christopher J. Wood, Jake Lishman, Julien Gacon, Simon Martiel, Paul D. Nation, Lev S. Bishop, Andrew W. Cross, Blake R. Johnson, and Jay M. Gambetta. Quantum computing with Qiskit, 2024

    work page 2024

  21. [21]

    Nation, Hwajung Kang, Neereja Sundaresan, and Jay M

    Paul D. Nation, Hwajung Kang, Neereja Sundaresan, and Jay M. Gambetta. Scalable mitigation of measurement errors on quantum computers.PRX Quantum, 2:040326, 2021

    work page 2021

  22. [22]

    Gambetta

    Kristan Temme, Sergey Bravyi, and Jay M. Gambetta. Error mitigation for short-depth quantum circuits.Physical Review Letters, 119:180509, 2017

    work page 2017

  23. [23]

    Benjamin

    Ying Li and Simon C. Benjamin. Efficient variational quantum simulator incorporating active error minimization.Physical Review X, 7:021050, 2017

    work page 2017

  24. [24]

    Córcoles, Antonio Mezzacapo, Jerry M

    Abhinav Kandala, Kristan Temme, Antonio D. Córcoles, Antonio Mezzacapo, Jerry M. Chow, and Jay M. Gambetta. Error mitigation extends the computational reach of a noisy quantum processor.Nature, 567:491–495, 2019. 16

    work page 2019