---
title: Center-of-Mass Imaging in Heavy-Ion Collisions
url: https://www.emergentmind.com/topics/center-of-mass-imaging
type: topic
---

# Center-of-Mass Imaging in Heavy-Ion Collisions

Center-of-Mass Imaging is a femtoscopic source-imaging strategy for identical-pion correlations in heavy-ion collisions that reconstructs the source directly in the source center-of-mass frame (CMFS) rather than in the conventional pair center-of-mass frame (CMFP). Its motivation is that standard imaging in the CMFP is mathematically clean because the pair energy difference vanishes there, but each pion pair has its own CMFP, so the reconstructed source is assembled from many pair-dependent Lorentz frames. The CMFS formulation instead aims to recover a source image in a single fixed frame, more directly tied to the actual spatial scale and morphology of the emitter, by suppressing the residual temporal phase term through a small-\(\Delta E\) pair selection [1806.05061, 2511.01251].

## 1. Frame choice and physical meaning

In two-pion femtoscopy, the object of interest is the source function, namely the distribution of relative separations between emitted pions. Historically, Brown–Danielewicz-type imaging reconstructs this source from measured correlation functions without imposing a Gaussian ansatz. In the conventional formulation, the reconstruction is carried out in the CMFP, the frame in which the total three-momentum of a given pion pair vanishes. That choice is exact for the standard inversion formula, but it is also geometrically awkward because different pairs have different CMFP boosts; the resulting image therefore mixes many different pair-dependent Lorentz frames rather than representing a single unified source frame [1806.05061].

Center-of-Mass Imaging replaces that pair-wise frame assignment by the CMFS, the global frame associated with the emitting source itself. In that frame the source geometry is fixed, so the reconstructed image is intended to reflect the emitter more directly. Later work sharpened this point by arguing that part of the long non-Gaussian tail seen in RHIC pion imaging is methodological rather than purely physical: CMFP imaging can itself broaden the source because the reconstructed separations are evaluated after different Lorentz transformations for different pairs [2511.01251].

The frame distinction is not merely terminological. For a spherical source of radius \(10\) fm, the 2018 study found that model source functions in the CMFS are nearly identical for two different momentum distributions, whereas in the CMFP they differ visibly; it also emphasized that the CMFS source vanishes beyond roughly the source diameter, while the CMFP source can extend to much larger separations, even beyond \(25\) fm [1806.05061].

## 2. Correlation formalism and the CMFS obstruction

Neglecting final-state interactions, the two-pion correlation function is written as
\[
C(\mathbf q)=1+\int d^4x\,\cos(\mathbf q\cdot \mathbf r-\Delta E\,\Delta t)\, d(x),
\]
where \(d(x)\) is the relative space-time separation distribution, \(\mathbf q\) is the relative momentum, \(\mathbf r\) is the relative spatial separation, and \(\Delta t\) is the emission-time difference. In the CMFP one has \(\Delta E=0\), so the correlation reduces to
\[
C(\mathbf q)=1+\int \cos(\mathbf q\cdot \mathbf r)\,S(\mathbf r)\,d\mathbf r,
\qquad
S(\mathbf r)=\int d(x)\,dt.
\]
After angle averaging, the one-dimensional imaging relation used in both papers is
\[
\mathcal R(q)\equiv C(q)-1 =4\pi\int \frac{1}{q}\sin(qr)\, r\, S(r)\,dr,
\]
with inversion
\[
S(r)=\frac{1}{2\pi^2}\frac{1}{r}\int \mathcal R(q)\, q\,\sin(qr)\,dq.
\]
The 2025 paper describes this as the one-dimensional Koonin-Pratt–type relation in the simplified noninteracting, spherically symmetric case, while the 2018 paper situates it within standard Brown–Danielewicz imaging [1806.05061, 2511.01251].

The CMFS difficulty is immediate: in general \(\Delta E\neq 0\), so the phase \(\Delta E\,\Delta t\) cannot be dropped exactly. Physically, the time extent of emission contaminates the apparent spatial separation. The standard CMFP treatment avoids this because the kernel depends only on space there; the CMFS treatment does not. A plausible implication is that, in a moving pair frame, unequal emission times shift the effective separation by a velocity-times-time term, but the CMFS program avoids reconstructing that pair-dependent quantity directly and instead tries to suppress the temporal phase itself [1806.05061].

## 3. Operational definition of the method

The central approximation is to select identical pion pairs with very small energy difference
\[
\Delta E\equiv |E_1-E_2|,
\]
so that
\[
\Delta E\,\Delta t \ll 1
\]
and hence
\[
\cos(\mathbf q\cdot \mathbf r-\Delta E\,\Delta t)\approx \cos(\mathbf q\cdot \mathbf r).
\]
Under that condition, the CMFS correlation function approximately obeys the same inversion logic as in the CMFP:
\[
C(\mathbf q)\approx 1+\int \cos(\mathbf q\cdot \mathbf r)\,S(\mathbf r)\,d\mathbf r.
\]
This small-\(\Delta E\) condition is the formal basis of Center-of-Mass Imaging in the pion-source context [1806.05061].

Operationally, one works in the CMFS, selects identical pion pairs, computes \(\mathbf q\), and imposes a cut such as
\[
\Delta E<30~\text{MeV}
\qquad\text{or}\qquad
\Delta E<10~\text{MeV}.
\]
Using only those pairs, one constructs \(C(q)\) in the CMFS, angle-averages to obtain \(\mathcal R(q)\), and reconstructs the one-dimensional source image \(S(r)\) by Fourier inversion. As in earlier imaging work, the benchmark “model source” for comparison is obtained from generated pion pairs with relative momentum less than \(60\) MeV/\(c\) [1806.05061].

The 2025 refinement places the conceptual emphasis on staying in one global source frame throughout. Its practical prescription remains the same in substance—construct the two-pion correlation function in the CMFS and retain only pairs satisfying
\[
\Delta E<10~\text{MeV},
\]
so that the neglected \(\Delta E\,\Delta t\) term is small—but it interprets the gain explicitly as elimination of frame mixing and reduction of kinematic correlations generated by varying pair boosts [2511.01251].

## 4. Model tests and AMPT validation

The 2018 paper first verified that standard imaging behaves as expected in the conventional frame. Using the string-melting AMPT model for central Au+Au collisions at \(\sqrt{s_{NN}}=200\) GeV with impact parameter \(b=0\) fm and \(200\) events, it found good agreement between the directly computed source function and the reconstructed one in the CMFP [1806.05061].

It then tested CMFS imaging in a controlled toy model: a homogeneous sphere of radius \(R=10\) fm, pion momenta drawn from a Boltzmann distribution with \(T_f=158\) MeV, and a one-sided Gaussian emission-time distribution
\[
f(t)=
\begin{cases}
\sqrt{2/\pi}\,\exp\!\left(-t^2/2\tau^2\right)/\tau, & t>0,\\
0, & t<0,
\end{cases}
\qquad \tau=10~\text{fm}/c.
\]
Using \(10^7\) pion pairs, the paper showed that CMFS imaging without any \(\Delta E\) cut gives a clear discrepancy from the true source. Once the cut is applied, the reconstructed source approaches the true source; for \(\Delta E<10\) MeV, the difference is described as “almost ignored.” Replacing the Boltzmann momentum distribution by a uniform one from \(0\) to \(1\) GeV/\(c\) produced the same trend, indicating that the effectiveness of the cut is only weakly influenced by the underlying momentum spectrum [1806.05061].

The same study then examined bias from the cut itself in AMPT. Comparing the full CMFS source function with the source restricted to pairs satisfying \(\Delta E<10\) MeV, it found that the energy cut alters the source only slightly. It further compared directly computed model sources with reconstructed CMFS images for three impact parameters—\(b=0\) fm with \(200\) events, \(b=5\) fm with \(500\) events, and \(b=10\) fm with \(10000\) events—and reported generally good agreement. The only noticeable mismatch occurred around \(r\approx 30\) fm in the \(b=0\) fm case, attributed to the discreteness of the measured correlation function rather than to failure of the method. The paper also reported that using an even smaller \(\Delta E\) cut did not materially change the result, indicating that \(\Delta E<10\) MeV was already sufficient in that study [1806.05061].

## 5. Non-Gaussian tails and the 2025 reformulation

The 2025 paper reinterprets a long-standing RHIC observation: source images reconstructed from two-pion correlations often show a pronounced non-Gaussian tail. Earlier analyses had usually associated that tail with an extended halo, for example from long-lived resonance decays. The 2025 result does not reject physical non-Gaussianity, but it separates two contributions: genuine source structure and frame-dependent distortion produced by CMFP imaging itself [2511.01251].

To expose the issue, the paper used a simulated Gaussian source
\[
S(r,t)= \frac{1}{\pi^2 r_0^3 \tau}
\exp\left( -\frac{r^2}{2r_0^2} -\frac{t^2}{2\tau^2} \right),
\]
with \(r_0=5\) fm, \(\tau=5\) fm/\(c\), Boltzmann momenta corresponding to \(T_f=130\) MeV, and a pair cut
\[
q<60~\text{MeV}/c.
\]
In the CMFP, the imaged source agreed with the model source in that frame, but the result displayed a stronger tail than a Gaussian source in the source frame would suggest. In the CMFS, restricting to pairs with \(\Delta E<10\) MeV produced an imaged source that agreed very well with the model source evaluated in the CMFS and was much closer to Gaussian [2511.01251].

The paper then turned to AMPT string-melting events for Au+Au collisions at \(\sqrt{s_{NN}}=200\) GeV, \(b=0\) fm, and identical \(\pi^+\pi^+\) pairs. It compared CMFP results obtained with \(200\) simulated events to CMFS results obtained with \(10000\) events, explicitly attributing the much larger CMFS event count to the restrictive small-\(\Delta E\) selection. The reconstructed CMFS source had a much smaller long-range tail than the CMFP source. The interpretation advanced there is that the CMFP tail contains both physical and frame-artifact contributions, whereas the CMFS tail reflects only the physical non-Gaussianity of the pion source. The evidence is presented visually rather than through goodness-of-fit measures, moments, or tail exponents [2511.01251].

## 6. Scope, limitations, and status

The method is explicitly limited to one-dimensional imaging. Both papers use a noninteracting correlation formalism, so in practice the construction applies to correlation functions built from Coulomb-corrected data. The 2018 paper describes the method as preliminary, and the 2025 paper states that the small-\(\Delta E\) approximation is validated empirically rather than by a more rigorous analytical error estimate [1806.05061, 2511.01251].

Statistics are a central practical constraint. The cut
\[
\Delta E<10~\text{MeV}
\]
strongly reduces the available pair sample, which is why large event counts appear in the AMPT tests, especially in the 2025 CMFS study. The 2025 paper also notes the absence of a systematic optimization study for the \(\Delta E\) threshold, and it emphasizes that its comparisons are primarily morphological. The 2018 paper adds that resolution and discreteness of the measured \(C(q)\) can affect the inversion, particularly at large \(r\) [1806.05061, 2511.01251].

Within those limits, Center-of-Mass Imaging establishes a clear methodological alternative to pair-frame imaging. Its defining claim is that, once the \(\Delta E\,\Delta t\) phase is sufficiently suppressed, the reconstructed source in the CMFS is tied to a single physical frame, is less sensitive to the pair momentum distribution, and yields a more direct representation of source geometry than the conventional CMFP image. The later reformulation further argues that this single-frame construction substantially reduces the artificial non-Gaussian tail induced by CMFP frame mixing, while leaving the genuinely physical tail available for interpretation [1806.05061, 2511.01251].

Source: https://www.emergentmind.com/topics/center-of-mass-imaging