On the CSL Scalar Field Relativistic Collapse Model

Daniel Bedingham; Philip Pearle

arxiv: 1906.11510 · v1 · pith:CIPZOXFHnew · submitted 2019-06-27 · 🪐 quant-ph

On the CSL Scalar Field Relativistic Collapse Model

Daniel Bedingham , Philip Pearle This is my paper

Pith reviewed 2026-05-25 15:10 UTC · model grok-4.3

classification 🪐 quant-ph

keywords CSL collapserelativistic invariancescalar fieldcoherent statesparticle productiondensity matrixdynamical reduction

0 comments

The pith

The relativistic CSL model with scalar field collapse operator causes a superposition of two particle clumps to collapse to one clump with high probability.

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

This paper adapts the CSL dynamical collapse model to a relativistically invariant setting by using a one-dimensional massive scalar field as the collapse-generating operator. For an initial normalized superposition of two widely separated coherent-state clumps, the exact solution of the density matrix when the Hamiltonian is zero (and the short-time solution otherwise) shows that the state evolves toward an eigenstate of the field operator close to one of the initial clump profiles. The result is essentially one surviving clump rather than a persistent superposition. Over longer times the same dynamics generates particles uniformly across all momenta, which the authors note makes the model experimentally unrealistic even though the collapse behavior itself is consistent.

Core claim

When the collapse-generating operator is the one-dimensional scalar field of mass m and the initial state is a superposition of two coherent clumps whose mean densities are N times squared Gaussians separated by much more than their width, the density matrix solution favors with high probability eigenstates whose eigenvalues are close to one of the initial clump profiles, thereby reducing the superposition to essentially a single clump of particles.

What carries the argument

The collapse-generating operator given by the one-dimensional scalar field operator acting on the initial coherent-state superposition.

If this is right

The model produces sensible collapse to a single clump on short timescales.
Particle production occurs at a constant rate per unit volume in every momentum mode.
The density matrix admits an exact solution for the zero-Hamiltonian case and short times otherwise.
The construction is relativistically invariant yet ruled out by the particle-production effect.

Where Pith is reading between the lines

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

The particle-production problem may be generic to relativistic collapse models that use local field operators.
Alternative choices of collapse operator could be explored to reduce or eliminate the unwanted energy generation while preserving the clump-selection behavior.
The short-time equivalence between zero and nonzero Hamiltonian cases suggests the collapse dynamics can be analyzed independently of the unitary evolution in the early regime.

Load-bearing premise

The collapse is generated by the scalar field operator on an initial state consisting of a superposition of two widely separated coherent-state clumps.

What would settle it

Preparing a macroscopic superposition of two such clumps and checking whether the final state is localized to one clump location while particles appear at a uniform rate across all momenta.

read the original abstract

The CSL dynamical collapse structure, adapted to the relativistically invariant model where the collapse-generating operator is a one-dimensional scalar field $\hat\phi(x,t)$ (mass $m$) is discussed. A complete solution for the density matrix is given, for an initial state $|\psi,0\rangle=\frac{1}{\sqrt{2}}[|L\rangle+|R\rangle]$ when the Hamiltonian $\hat H$ is set equal to 0, and when $\hat H$ is the free field Hamiltonian. Here $|L\rangle, |R\rangle$ are coherent states which represent clumps of particles, with mean particle number density $N\chi_{i}^{2}(x)$, where $\chi_{1}(x),\chi_{1}(x) $ are gaussians of width $\sigma>>m^{-1}$ with mean positions separated by distance $>>\sigma$. It is shown that, with high probability, the solution for $\hat H=0$ (identical to the short time solution for $\hat H\neq 0$) favors collapse toward eigenstates of the scalar field whose eigenvalues are close to $\sim\chi_{i}(x)$. Thus, this collapse dynamics results in essentially one clump of particles. However, eventually particle production dominates the density matrix since, as is well known, the collapse generates energy/sec-volume of every particle momentum in equal amounts. Because of the particle production, this is not an experimentally viable physical theory but, as is emphasized by the discussion, it is a sound relativistic collapse model, with sensible collapse behavior.

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.

Desk Editor's Note private letter to a colleague

The paper gives an explicit density-matrix solution for the H=0 relativistic scalar CSL model on a two-clump initial state and shows the expected short-time collapse to one clump.

read the letter

The useful piece is the closed-form density matrix for the zero-Hamiltonian case, which is also the short-time limit when the free field Hamiltonian is restored. Starting from a normalized superposition of two distant coherent-state clumps, the off-diagonal terms damp at a rate set by the squared difference of the field eigenvalues, so the state settles into a mixture of configurations close to one clump or the other. That matches the claim in the abstract and follows directly from integrating the master equation; the separation of the Gaussians makes the damping clean and the result transparent. The paper also states plainly that particle production eventually dominates, rendering the model non-viable, and treats that fact as already established in the literature rather than re-deriving it here. That is a limitation but not a contradiction with the short-time result they actually compute. The initial-state choice (clumps far apart, width much larger than Compton length) is reasonable for isolating the collapse dynamics, though it leaves open how the behavior changes for overlapping or more general states. No internal inconsistency appears in the central calculation. This is technical work aimed at people who need the explicit solution for this particular relativistic collapse construction. It supplies a missing concrete step even though the model is already known to fail on other grounds. I would send it to referees because the derivation is self-contained and the result is falsifiable in principle; a serious editor could usefully ask for more detail on the long-time particle-production term or a check on overlapping clumps, but the paper as written deserves that review.

Referee Report

1 major / 2 minor

Summary. The paper adapts the CSL collapse model to a relativistically invariant setting with the collapse-generating operator taken to be the one-dimensional scalar field φ̂(x,t) of mass m. It supplies an explicit solution of the master equation for the density matrix starting from a normalized superposition of two spatially separated coherent-state clumps |L⟩ and |R⟩ (mean densities Nχ_i²(x), Gaussian width σ ≫ m⁻¹, separation ≫ σ) both for Ĥ = 0 and at short times for the free-field Hamiltonian. The solution is shown to damp off-diagonal elements at a rate set by the squared eigenvalue difference of the collapse operator, driving the state toward a mixture of configurations close to one or the other χ_i(x) and thereby producing essentially a single clump; particle production is noted separately as eventually dominating and rendering the model non-viable, while the collapse dynamics itself is argued to be sound.

Significance. If the derivation is correct, the work supplies a parameter-free, explicit demonstration that the relativistic CSL dynamics produces the expected localization onto one clump without circularity or fitted parameters. The direct integration yielding the density-matrix evolution and the explicit short-time equivalence to the Ĥ = 0 case are concrete strengths that can be checked against the master equation.

major comments (1)

[Derivation of the density matrix (Ĥ = 0 case)] The central claim that the evolved density matrix favors eigenstates close to ∼χ_i(x) (and therefore one clump) rests on the explicit solution obtained by direct integration; the manuscript must display the key intermediate expression for the off-diagonal damping factor (proportional to the squared difference of the collapse-operator eigenvalues on the disjoint supports) so that the rate and the resulting mixture can be verified.

minor comments (2)

The initial-state definition should state the normalization factor and the precise functional form of χ_i(x) explicitly (including the relation between N and the coherent-state amplitude) to allow immediate reproduction of the setup.
The assertion that particle production 'is well known' would be strengthened by a one-sentence recap of the energy-generation rate per unit volume and momentum, even if only by reference to the standard CSL literature.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive evaluation and recommendation of minor revision. We address the single major comment below.

read point-by-point responses

Referee: [Derivation of the density matrix (Ĥ = 0 case)] The central claim that the evolved density matrix favors eigenstates close to ∼χ_i(x) (and therefore one clump) rests on the explicit solution obtained by direct integration; the manuscript must display the key intermediate expression for the off-diagonal damping factor (proportional to the squared difference of the collapse-operator eigenvalues on the disjoint supports) so that the rate and the resulting mixture can be verified.

Authors: We agree that an explicit display of the intermediate damping factor will make the derivation easier to verify. The off-diagonal elements of the density matrix acquire a factor exp(−½∫[Δλ(t′)]² dt′), where Δλ(t′) is the difference between the eigenvalues of the collapse operator φ̂ on the two disjoint supports of |L⟩ and |R⟩. In the revised manuscript we will insert this expression immediately after the statement of the master-equation solution, together with the resulting mixture weights for the two clumps. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper derives its central result—an explicit solution for the density matrix at Ĥ=0 (and short-time equivalence for Ĥ≠0)—by direct integration of the CSL master equation starting from the normalized superposition of two coherent states. The damping of off-diagonal elements follows immediately from the squared difference of the collapse-operator eigenvalues for disjoint Gaussian supports, with no reduction to a fitted parameter, self-defined quantity, or load-bearing self-citation. The subsequent observation of particle production is presented as a known separate consequence rather than part of the collapse derivation itself. The derivation chain is therefore self-contained against the model equations and initial state.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities beyond the standard CSL framework; the scalar-field operator and coherent-state preparation are taken as given inputs of the model construction.

pith-pipeline@v0.9.0 · 5800 in / 1171 out tokens · 33352 ms · 2026-05-25T15:10:31.341346+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The density matrix in the Schrödinger picture satisfies the Lindblad evolution equation dρ̂(t)/dt = −i[Ĥ,ρ̂(t)] − λ/2 ∫ dx φ̂(x,0)[φ̂(x,0),ρ̂(t)].
IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean absolute_floor_iff_bare_distinguishability unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

It is shown that, with high probability, the solution for Ĥ=0 favors collapse toward eigenstates of the scalar field whose eigenvalues are close to ∼χi(x).

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

9 extracted references · 9 canonical work pages

[1]

(2.5) We will need to know |0⟩k and |ℓ⟩k in the |Xk⟩|pk⟩ basis

Therefore, we have to solve Eq.(1.7) for ˆρ11k(t), ˆρ12k(t), ˆρ21k(t), ˆρ22k(t), with corresponding initial conditions ˆρss′k(0) = |ℓs⟩k k⟨ℓs′|. (2.5) We will need to know |0⟩k and |ℓ⟩k in the |Xk⟩|pk⟩ basis. For |0⟩k, since ( ˆXk +i 1 2 ˆPk)|0⟩k = 0, (ˆpk − i 1 2 ˆxk)|0⟩k = 0: ⟨Xk|⟨pk|0⟩k = √ 2 πe− X 2 ke− p2 k. (2.6) For |ℓ⟩k, we apply the Campbell-Bake...

work page
[2]

Density Matrix in the Position Representation We now proceed to calculate the matrix element of (A10), using eSˆa† LˆaR|0⟩⟨0|= ∑∞ n=0Sn|n⟩⟨n|, whose matrix elements in the position representation are given by a w ell-known identity involving Hermite polynomials [5]. We also use the Campbell-Baker-Hauss dorf theorem, obtaining: ⟨X|ˆρ(t)|X ′⟩ =C(t)e(1− S)γ1...

work page
[3]

Short Time and Long Time Limits For short times, S ≈ λ 2ωt<< 1, (A12) becomes ⟨X|ˆρ(t)|X ′⟩ ≈ √ 2 πe− λt ω [X− X ′]2 e− [X− γ1]2 e− [X ′− γ2]2 . (A13) It is consistent to neglect the exponent in the last factor in (A12), e S 1− S [γ1− γ2]2 ≈ e− λt 2ω [γ1− γ2]2 ≈ 1 to accompany the approximation Sγ i << γi, so that, even in this approximation, the property...

work page
[4]

Harmonic Oscillator and Collapse Generated by ˆp The other harmonic oscillator problem in (1.6), d dt ˆρ(t) = −iω [ 1 4 ˆx2 + ˆp2, ˆρ(t)] − λ ω [ˆp, [ˆp, ˆρ(t)]] } = −iω [ˆa†ˆa,ρ (t)] − λ 4ω [ˆa† + ˆa, [ˆa† + ˆa, ˆρ(t)], (A16) subject to the initial condition ˆρ(0) = e− 1 2 γ ′2 1 eγ ′ 1ˆa† |0⟩⟨0|eγ ′ 2ˆae− 1 2 γ ′2 2 , (A17) has precisely the same soluti...

work page
[5]

Pearle, P., ‘Combining stochastic dynamical state-vec tor reduction with spontaneous localiza- tion’, Phys. Rev. A 39, 2277 (1989)

work page 1989
[6]

C., Pearle, P

Ghirardi, G. C., Pearle, P. and Rimini, A., ‘Markov proce sses in Hilbert space and continuous spontaneous localization of systems of identical particle s’, Phys. Rev. A 42, 78 (1990)

work page 1990
[7]

Pearle, P., ‘Toward a Relativistic Theory of Statevecto r Reduction’ in Sixty-Two Years of Uncertainty, ed. A. Miller (Plenum, New York), p. 193 (1990). 17

work page 1990
[8]

C., Grassi, R

Ghirardi, G. C., Grassi, R. and Pearle, P., ‘Relativisti c Dynamical Reduction Models: General Framework and Examples,’ Found. Phys. 20, 1271 (1990)

work page 1990
[9]

However, t he usual annihilation operator is a = 1√ 2 [ ˆX + i ˆP ] whereas the one used here is a = [ ˆX + i 1 2 ˆP ]

The identity (e.g., in the last section of Wikipedia’s ar ticle on Hermite polynomials) is ∞∑ n=0 Snψn(x)ψn(y) = 1 √ π[1 − S2] [ e− 1− S 1+S (x+y)2 4 + e− 1+S 1− S (x− y)2 4 ] (A19) where ψn(x) is the usual harmonic oscillator wave function. However, t he usual annihilation operator is a = 1√ 2 [ ˆX + i ˆP ] whereas the one used here is a = [ ˆX + i 1 2 ˆ...

work page

[1] [1]

(2.5) We will need to know |0⟩k and |ℓ⟩k in the |Xk⟩|pk⟩ basis

Therefore, we have to solve Eq.(1.7) for ˆρ11k(t), ˆρ12k(t), ˆρ21k(t), ˆρ22k(t), with corresponding initial conditions ˆρss′k(0) = |ℓs⟩k k⟨ℓs′|. (2.5) We will need to know |0⟩k and |ℓ⟩k in the |Xk⟩|pk⟩ basis. For |0⟩k, since ( ˆXk +i 1 2 ˆPk)|0⟩k = 0, (ˆpk − i 1 2 ˆxk)|0⟩k = 0: ⟨Xk|⟨pk|0⟩k = √ 2 πe− X 2 ke− p2 k. (2.6) For |ℓ⟩k, we apply the Campbell-Bake...

work page

[2] [2]

Density Matrix in the Position Representation We now proceed to calculate the matrix element of (A10), using eSˆa† LˆaR|0⟩⟨0|= ∑∞ n=0Sn|n⟩⟨n|, whose matrix elements in the position representation are given by a w ell-known identity involving Hermite polynomials [5]. We also use the Campbell-Baker-Hauss dorf theorem, obtaining: ⟨X|ˆρ(t)|X ′⟩ =C(t)e(1− S)γ1...

work page

[3] [3]

Short Time and Long Time Limits For short times, S ≈ λ 2ωt<< 1, (A12) becomes ⟨X|ˆρ(t)|X ′⟩ ≈ √ 2 πe− λt ω [X− X ′]2 e− [X− γ1]2 e− [X ′− γ2]2 . (A13) It is consistent to neglect the exponent in the last factor in (A12), e S 1− S [γ1− γ2]2 ≈ e− λt 2ω [γ1− γ2]2 ≈ 1 to accompany the approximation Sγ i << γi, so that, even in this approximation, the property...

work page

[4] [4]

Harmonic Oscillator and Collapse Generated by ˆp The other harmonic oscillator problem in (1.6), d dt ˆρ(t) = −iω [ 1 4 ˆx2 + ˆp2, ˆρ(t)] − λ ω [ˆp, [ˆp, ˆρ(t)]] } = −iω [ˆa†ˆa,ρ (t)] − λ 4ω [ˆa† + ˆa, [ˆa† + ˆa, ˆρ(t)], (A16) subject to the initial condition ˆρ(0) = e− 1 2 γ ′2 1 eγ ′ 1ˆa† |0⟩⟨0|eγ ′ 2ˆae− 1 2 γ ′2 2 , (A17) has precisely the same soluti...

work page

[5] [5]

Pearle, P., ‘Combining stochastic dynamical state-vec tor reduction with spontaneous localiza- tion’, Phys. Rev. A 39, 2277 (1989)

work page 1989

[6] [6]

C., Pearle, P

Ghirardi, G. C., Pearle, P. and Rimini, A., ‘Markov proce sses in Hilbert space and continuous spontaneous localization of systems of identical particle s’, Phys. Rev. A 42, 78 (1990)

work page 1990

[7] [7]

Pearle, P., ‘Toward a Relativistic Theory of Statevecto r Reduction’ in Sixty-Two Years of Uncertainty, ed. A. Miller (Plenum, New York), p. 193 (1990). 17

work page 1990

[8] [8]

C., Grassi, R

Ghirardi, G. C., Grassi, R. and Pearle, P., ‘Relativisti c Dynamical Reduction Models: General Framework and Examples,’ Found. Phys. 20, 1271 (1990)

work page 1990

[9] [9]

However, t he usual annihilation operator is a = 1√ 2 [ ˆX + i ˆP ] whereas the one used here is a = [ ˆX + i 1 2 ˆP ]

The identity (e.g., in the last section of Wikipedia’s ar ticle on Hermite polynomials) is ∞∑ n=0 Snψn(x)ψn(y) = 1 √ π[1 − S2] [ e− 1− S 1+S (x+y)2 4 + e− 1+S 1− S (x− y)2 4 ] (A19) where ψn(x) is the usual harmonic oscillator wave function. However, t he usual annihilation operator is a = 1√ 2 [ ˆX + i ˆP ] whereas the one used here is a = [ ˆX + i 1 2 ˆ...

work page