Time-Dependent Wave-Packet Approach to Reactions for Quantum Computation

arXiv:2603.26881 · nucl-th, quant-ph · Submitted 2026-03-27 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: We're looking at the paper: "Time-Dependent Wave-Packet Approach to Reactions for Quantum Computation".

Mira: The gist: This method describes a time-dependent wave-packet approach to reactions for quantum computation,

Kai: First, who's behind it and why it matters.

Title and authors: Kai: So we're looking at this paper called "Time-Dependent Wave-Packet Approach to Reactions for Quantum Computation" and the authors are Evan Rule and Ionel Stetcu from Los Alamos National Laboratory. It sounds like they're describing a way to get the scattering matrix for nuclear or chemical reactions on a finite lattice.

Mira: It’s interesting because they say you only need unitary time evolution after preparing those initial and final states as wave packets, which makes it ideal for simulations on quantum hardware.

Lev: Unitary time evolution is key, but we always have to think about the overhead of that on real hardware. What exactly are these wave packets actually built and cooled to look like?

Kai: Well, they use a first quantization mapping in Cartesian coordinates, which they say allows for efficient encoding on quantum hardware and favorable qubit scaling for describing asymptotic scattering states. It’s all about how you map the physics onto the qubits.

Mira: And they introduce this central quantity, which is a time-dependent overlap between incoming and outgoing wave packets whose Fourier transform corresponds to the scattering matrix at fixed energy, from which you can calculate elastic and inelastic cross sections for reactions involving two interacting clusters.

Lev: So they’re giving us a quantity that lets us find the cross sections after some time evolution. What assumptions are they making about the interaction?

Kai: They assume the intercluster interaction is finite-range, meaning it vanishes or is exponentially small when the relative separation between two clusters is sufficiently large. That keeps things manageable for them.

Mira: That's a big assumption because it sets a boundary on how far apart these particles can be before they interact. This paper really focuses on reactions involving two interacting clusters, which could be nuclear or chemical ones.

Lev: So if we take that finite-range interaction and the time evolution, what does that tell us about the actual computational cost?

Kai: The method is designed to require only unitary time evolution, which is great for quantum hardware because it doesn't need complex non-unitary operations to simulate the reaction dynamics.

Mira: It’s a lot of math condensed into this idea of sampling that overlap function and taking a Fourier transform classically, which they say is trivial.

Lev: Trivial classical transforms are good, but how do we actually get those overlap function samples from the quantum computer? That sounds like a major hurdle for implementation.

Kai: They suggest you sample the overlap function at sufficiently many times set by the momentum and energy resolution of the lattice, which then makes that Fourier transform trivial to perform classically.

The paper's summary: Mira: To summarize, this is about using a time-dependent wave-packet approach to reactions for quantum computation where you prepare initial and final states as wave packets and use unitary time evolution to find the scattering matrix.

Kai: Exactly, it’s setting up these states k0, alpha using internal cluster states and a relative Gaussian wave packet phi rel(k zero r), where k zero is defined by (k zero times, k 0,y, k 0,z) <ref:2603.26881#pg3>.

Lev: So they’re defining the state by the internal degrees of freedom of each cluster and a relative component described by that Gaussian wave packet phi rel(k zero r) which has spatial width sigma along the z axis and propagates with group momentum k 0,z <ref:2603.26881#pg3>.

Kai: Along the x and y directions, the relative wave function occupies momentum eigenstates with respective eigenvalues k zero times and k 0,y, which is how they encode those transverse momenta <ref:2603.26881#pg3>.

Mira: The total state is written as k0, alpha psi(one) alpha one psi(two) alpha two phi rel(k zero r), where alpha denotes the set of quantum numbers that completely specify the internal state of each cluster.

Lev: They are using a zero CM momentum state to isolate the desired relative momentum eigenstate, which means they factorize the preparation between clusters by using separate wave packets for each and projecting onto a zero CM eigenstate. That sounds like a way to simplify things computationally.

Kai: They also mention that they can reduce complexity by factoring the preparation between clusters by using separate wave packets for each cluster and then projecting onto a zero CM eigenstate to isolate the desired relative momentum eigenstate, which is page eight of this paper.

Mira: And if you sample that overlap function at enough times, you get the scattering matrix from a Fourier transform, and they extract the real and imaginary parts of that overlap function using modified Hadamard tests.

Lev: So it sounds like the main computational burden shifts to preparing those specific wave packets on the quantum computer, because once that's done, the evolution is just unitary time evolution.

Kai: Right. And they say this procedure requires propagation of the initial and final wave packets to infinite times t to plus or minus infinity, which isn't feasible.

The paper's improvements: Mira: The authors suggest a way to get around that infinite time problem, saying it suffices for time t such that... they can compute the scattering matrix.

Kai: They mention that because the wave packets are extended plane waves in x and y, the contributions to the cross section from all impact parameters are automatically included, which is page three of this paper <ref:2603.26881#pg2>.

Lev: So if we don't need infinite time, what does that mean for simulating a real reaction? Does it just mean we can get away with a finite simulation time?

Mira: It means they can compute the total and differential scattering cross sections for two to two scattering, including inelastic processes <ref:2603.26881#pg2>. They handle these by extending the formalism to a two-channel problem using a Hamiltonian involving coupling terms.

Kai: They also mention that inelastic scattering is handled by extending the formalism to a two-channel problem using a Hamiltonian involving coupling terms, which allows for calculating elastic and inelastic differential cross sections.

Lev: Extending it to two channels means you’re adding more complexity because you have those coupling terms in the Hamiltonian, right? That increases the complexity of running this on real hardware.

Mira: It does increase complexity, but it's necessary to get both elastic and inelastic differential cross sections for reactions involving composite systems.

Kai: And they also say this formalism is general because it is agnostic to the nature of the interaction, as long as it is short ranged, which means any short-ranged interaction works.

Lev: So if any short-ranged interaction works, does that mean we can apply this to anything from nuclear reactions to chemical reactions?

Mira: Yes, they state that it can be extended to simulate nuclear and chemical reactions where both target and projectile can be composite systems, which opens the way for simulating these kinds of things.

Conclusion: Kai: So we’ve talked about how this time-dependent wave-packet approach to reactions for quantum computation works, relying on unitary time evolution after preparing those specific states k0, alpha.

Mira: The main implication is that it gives us a way to obtain the scattering matrix for nuclear or chemical reactions on a finite lattice.

Lev: From my side, the real challenge is making sure that if we run this on real hardware, we can actually handle those coupling terms and get the overhead right for error correction.

Kai: So they manage to address issues like needing finite simulation time and handling impact parameters automatically, which helps make the signal distinguishable from noise.

Mira: Ultimately, this method is ideal for quantum computation because it requires only unitary time evolution, making it ideal for simulations on quantum hardware.

Lev: I just think we need to make sure the two-channel problem setup doesn't just balloon the required qubit count too much.

Kai: So that's what we have here with this paper, "Time-Dependent Wave-Packet Approach to Reactions for Quantum Computation". It gives us a solid framework for simulating these types of reactions.

Los Alamos National Laboratory

nucl-th, quant-ph

Submitted: 2026-03-27

Updated: 2026-10-19

Comments: 5+11 pages, 4 figures

Journal ref: Phys. Rev. Lett. 137, 162501 (2026)

DOI: 10.1103/kwgf-jdbn

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 74/100

The gist: The gist: This method describes a time-dependent wave-packet approach to reactions for quantum computation, which allows for obtaining the scattering matrix for nuclear or chemical reactions on a

Key concepts

Time-Dependent Overlap
This is the central quantity calculated by comparing incoming and outgoing wave packets at different times. Its Fourier transform directly corresponds to the scattering matrix at a fixed energy. By tracking this overlap through unitary time evolution, researchers can extract information about how particles scatter.
Unitary Time Evolution
The method simulates the reaction process by evolving initial and final states using the system's Hamiltonian over time. Since this evolution is unitary, it is perfectly suited for implementation on quantum hardware, making it a practical tool for quantum simulations.
Wave Packet State Preparation
The system state is prepared as a product of internal cluster states and a relative Gaussian wave packet. This involves defining the initial and final configurations using specific internal quantum numbers ($\alpha$ and $\beta$) to describe the relative motion between clusters.

Terminology

Summary

The gist: This method describes a time-dependent wave-packet approach to reactions for quantum computation, which allows for obtaining the scattering matrix for nuclear or chemical reactions on a finite lattice by utilizing unitary time evolution.

Core Methodology

The central quantity in this approach is defined as a time-dependent overlap between incoming and outgoing wave packets whose Fourier transform corresponds to the scattering matrix at fixed energy (Page 1). This technique is ideal for quantum hardware because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). The method involves preparing initial and final states as wave packets and then performing the unitary time evolution between them (Page 7).

Formalism and State Preparation

The formalism is based on the "time-dependent scattering formalism of Tannor & Weeks [22]" (Page 1). The total system Hamiltonian is given by a sum involving kinetic energy operators, internal Hamiltonians, and intercluster interaction terms (Page 1). The system state is prepared as a product of internal cluster states and a relative Gaussian wave packet: Φk0,α⟩ ≡ ψ(1)α1⟩ ψ(2)α2⟩ ϕrel(k0, r) (Page 2). To reduce complexity, the approach can factorize the preparation between clusters by using separate wave packets for each cluster and then projecting onto a zero CM eigenstate to isolate the desired relative momentum eigenstate (Page 8).

Extracting Cross Sections

The scattering matrix S is related to the dimensionless overlap function Cβ,α(k′0, k0;t) (Page 11). If this overlap function is sampled at sufficiently many times, a subsequent Fourier transform is trivial to perform classically (Page 12). The real and imaginary parts of the overlap function can then be extracted from modified Hadamard tests [23] (Page 5). The scattering amplitude fβ,α(E, θ) can be obtained by relating the S matrix to this overlap function through Fourier transform and the two-dimensional analog of Eq. (12) (Page 10).

Numerical Implementation and Results

The method is demonstrated with a series of numerical examples, for both elastic and inelastic processes, comparing against exact calculations (Page 1). For elastic scattering in two dimensions via a Gaussian potential, the total cross section λ is obtained from the forward-scattering amplitude via the optical theorem: λ = √8π/k Imf(E, 0) (Page 11). The approach shows that the resulting cross section is large enough to be distinguished from the noise because extended plane waves automatically include contributions from all impact parameters (Page 3). Inelastic scattering is handled by extending the formalism to a two-channel problem using a Hamiltonian involving coupling terms, allowing for the calculation of elastic and inelastic differential cross sections (Page 5).

General Applicability

The approach is general because it is agnostic to the nature of the interaction, as long as it is short ranged (Page 11). It can be extended to simulate nuclear and chemical reactions where both target and projectile can be composite systems (Page 1). The use of a finite lattice provides scattering information across the entire range of energies where wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The formalism also allows for the simulation of charged-particle reactions and photo-induced reactions by extending it to include Coulomb interactions (Page 11).

Summary

This paper describes a time-dependent wave-packet approach to scattering that is ideal for quantum computation in the first quantization mapping. The approach, however, is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support (Page 6). The method is ideal for simulations on quantum hardware because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1). This formalism is general—both the target and the projectile can be composite systems, opening the way for simulation of nuclear and chemical reactions (Page 1). The use of wave packets localized to a finite lattice provides scattering information across the entire range of energies where the wave packets have appreciable support, with the spatial extent of the lattice being increased to simulate lower energies (Page 6). The method is ideal for quantum computation because it requires only unitary time evolution, making this approach ideal for simulations on quantum hardware (Page 1).

Improvements for AI systems

  1. A time-dependent wave-packet approach to reactions for quantum computation allows AI systems to calculate scattering cross sections by preparing initial and final states as wave packets and performing unitary time evolution, making it ideal for simulations on quantum hardware. This enables the calculation of elastic and inelastic cross sections for reactions involving two interacting clusters.

  2. The method can be adapted to describe reactions where constituents transition between internal eigenstates, enabling the calculation of inelastic scattering amplitudes by preparing states such that the initial and final states are labeled by quantum numbers α = (α1, α2) and β = (β1, β2), respectively.

  3. AI systems can probe scattering through different angles by imparting transverse momentum to the final-state wave packets, allowing for the calculation of differential cross sections dλ/dθ at various energies.

  4. The method provides a way to obtain total and differential scattering cross sections for 2 → 2 scattering, including inelastic processes, which can be used to model nuclear or chemical reactions involving composite systems.

  5. By utilizing the overlap function Cβ,α(k'0, k0;t), AI can compute the S matrix by sampling it at sufficiently many times (set by the momentum/energy resolution of the lattice), and then performing a Fourier transform that determines the scattering matrix.

Sources

Related papers