Case Study 2: A Day in the Life of a Nuclear Physics Experiment

The Experiment: Measuring the First Excited State of ${}^{54}\text{Ca}$

This case study follows a (composite but realistic) nuclear physics experiment from conception to publication, based on the kind of measurement routinely carried out at radioactive beam facilities. The experiment aims to measure the energy of the first $2^+$ excited state of ${}^{54}\text{Ca}$ ($Z = 20$, $N = 34$) — a nucleus that tests whether $N = 34$ is a new magic number in calcium.

The physics motivation is sharp: shell model calculations with modern interactions predict that the $N = 34$ subshell closure produces a large energy gap between the neutron $p_{3/2}$ and $p_{1/2}$ orbitals, which should manifest as a high $2^+_1$ energy in ${}^{54}\text{Ca}$ — higher than in ${}^{52}\text{Ca}$ ($N = 32$), even though ${}^{54}\text{Ca}$ has more neutrons. This would be a direct signature of a new magic number.

Phase 1: The Proposal (Month 0 – Month 6)

The Idea

Dr. Elena Vasquez, an assistant professor at a midwestern university, has been studying shell evolution in neutron-rich calcium isotopes for three years. Her theoretical collaborators have published predictions: $E(2^+_1)$ in ${}^{54}\text{Ca}$ should be approximately 2.0 MeV if $N = 34$ is magic, compared to $\sim 1.0$ MeV if it is not. The measurement would decisively test the prediction.

Elena assembles a collaboration: 28 physicists from 9 institutions, including experimentalists with expertise in gamma-ray spectroscopy (her own group and a German group), reaction theorists, and the GRETINA operations team.

The Proposal Document

The collaboration writes a 15-page proposal specifying:

• Physics case: Three pages reviewing the theoretical predictions, existing data on calcium isotopes through ${}^{52}\text{Ca}$, and why ${}^{54}\text{Ca}$ is the critical test case.
• Reaction mechanism: The $2^+$ state will be populated by proton inelastic scattering ${}^{54}\text{Ca}(p,p')$ at $\sim 70$ MeV/u using a liquid hydrogen target. The de-excitation gamma ray will be detected by GRETINA.
• Rate estimate: Using FRIB's predicted production rate of $\sim 500$ ${}^{54}\text{Ca}$ ions/second (from fragmentation of ${}^{76}\text{Ge}$), and an inelastic scattering cross section of $\sim 20$ mb, the collaboration estimates 2.5 gamma-ray events per hour in GRETINA. To achieve the needed statistics ($\sim 200$ events for a clear peak), they request 80 hours of beam on target — plus 48 hours for setup, calibration, and contingency. Total request: 5.5 days.
• Beam specifications: ${}^{54}\text{Ca}$ beam at 60–80 MeV/u, purity $> 80\%$, minimum rate 200/s.
• Detector setup: GRETINA at the target position, with the S800 spectrograph downstream for coincident particle identification.
• Simulations: GEANT4 simulations of the setup showing the expected gamma-ray spectrum (Doppler-corrected), the estimated peak-to-background ratio, and the detection efficiency ($\sim 6\%$ at 2 MeV for GRETINA at this position).

The PAC

The proposal is submitted to FRIB's Program Advisory Committee (PAC), which meets twice a year. The PAC consists of 12–15 nuclear physicists from institutions worldwide, selected for their broad expertise. They evaluate each proposal on:

1. Scientific merit and impact: Is this an important question? Will the result be published in a high-impact journal?
2. Technical feasibility: Are the rate estimates realistic? Is the analysis strategy sound?
3. Efficiency: Is the requested beam time well justified?

Elena presents the proposal in a 15-minute talk followed by 10 minutes of questions. The committee is enthusiastic about the physics but skeptical about the rate estimate — they ask for a more conservative calculation assuming only 300 ${}^{54}\text{Ca}$/s. Elena shows that 100 hours of beam time would still yield $\sim 150$ events, enough for a clear measurement but with larger uncertainties.

The PAC recommends approval with A priority (schedule as soon as possible), allocating 7 days of beam time.

Phase 2: Preparation (Month 6 – Month 12)

Detector Setup

Six months before the beam time, the preparation accelerates:

• GRETINA configuration: The 12 GRETINA detector modules are arranged in a hemisphere around the target position. Each module is tested: high-voltage bias applied, energy calibration with ${}^{152}\text{Eu}$ and ${}^{60}\text{Co}$ sources, timing calibration, segment response verification.
• Liquid hydrogen target: A 30-mm-thick target cell with thin Kapton windows, operated at 18 K. The target must be leak-tested, the cryocooler installed, and the safety interlocks verified.
• S800 spectrograph: Configured for the expected ${}^{54}\text{Ca}$ beam rigidity, with focal-plane detectors (ionization chamber for $\Delta E$, plastic scintillators for TOF) calibrated.
• Data acquisition: The digital data acquisition system is configured to read out GRETINA (12 modules $\times$ 4 crystals $\times$ 36 segments = 1,728 channels) and the S800 focal plane in coincidence. The trigger requires at least one GRETINA module in coincidence with a particle in the S800.
• GEANT4 simulations: Updated with the actual detector positions (measured by laser survey) to produce the final efficiency curves and Doppler correction parameters.

Analysis Software

Elena's graduate student, Marcus, has spent six months developing the analysis pipeline: - Raw signal processing (trapezoidal filtering, timing, pileup rejection) - Gamma-ray tracking algorithm (reconstruct Compton scattering sequences) - Doppler correction (using the beam velocity measured by the S800 and the gamma-ray emission angle from tracking) - Particle identification cuts (to ensure the gamma ray came from ${}^{54}\text{Ca}$ and not a contaminant) - Background subtraction and peak fitting

All of this is coded, tested on simulated data, and ready before the beam arrives.

Phase 3: The Beam Time (7 Days)

Day 1: Beam Development

The FRIB operations team tunes the primary beam: ${}^{76}\text{Ge}^{30+}$ accelerated to 150 MeV/u. The beam strikes the beryllium production target, and the ARIS fragment separator is tuned to transmit ${}^{54}\text{Ca}$ — adjusting dozens of magnets, slits, and the wedge degrader to optimize the secondary beam purity and rate.

By mid-afternoon, a secondary beam is delivered to the experimental area. The particle identification plot on the S800 focal plane shows a clear ${}^{54}\text{Ca}$ cluster, along with neighboring isotopes (${}^{53}\text{Ca}$, ${}^{55}\text{Sc}$, etc.). The ${}^{54}\text{Ca}$ rate is 350 ions/second — less than the proposal estimate but within the range the PAC considered.

The team fills the liquid hydrogen target (a carefully choreographed safety procedure), verifies the cryogenic stability, and takes calibration data with a ${}^{152}\text{Eu}$ source to check GRETINA's performance in the experimental configuration.

Day 2: First Physics Data

At 6:00 AM, the shift team (three physicists, working 8-hour shifts around the clock) starts recording physics data. The first Doppler-corrected gamma-ray spectrum appears on the monitor after about four hours, showing a possible peak near 2.0 MeV — but with only $\sim 15$ counts, it is too early to be sure.

Shift work: Every 8 hours, a new team takes over. The shift leader monitors: - Beam rate (plotted in real time from the S800 particle rate) - Detector live time (the fraction of time the DAQ is not busy — typically $> 90\%$) - Gamma-ray spectrum (accumulating slowly) - Cryogenic systems (liquid hydrogen temperature and pressure)

Shift logs are recorded in an electronic logbook, and a twice-daily phone meeting keeps the full collaboration informed.

Day 3: The Crisis

At 2:30 AM, the night shift reports a sudden drop in beam rate from 350 to 40 ions/second. The FRIB operators diagnose the problem: the ECR ion source has lost plasma stability, reducing the primary beam current by a factor of 10.

For three hours, the experiment collects almost no data while the source group works to restore the plasma. By dawn, the beam is back — but the lost time represents $\sim 40$ gamma-ray events that will not be recovered.

Elena sends a message to the PAC liaison requesting 12 additional hours of beam time if available. The request is noted but cannot be guaranteed — the next experiment is already scheduled.

Days 4–6: Steady Running

The beam stabilizes at 300–400 ions/second. Data accumulates steadily. The gamma-ray spectrum grows:

By Day 6, the peak at $E_\gamma = 2.04 \pm 0.02$ MeV is clearly visible, sitting above a smooth background of Compton events and target-related gamma rays. The $2^+$ energy is consistent with the theoretical prediction for a shell closure at $N = 34$.

Day 7: Final Data and Teardown

The last 12 hours are split between final physics data collection (bringing the total to $\sim 205$ events in the peak) and systematic checks: empty-target runs (to characterize target-out background), calibration source runs, and timing calibration.

At 8:00 PM, the FRIB operator announces "beam off." The experiment is over. The team drains the liquid hydrogen target, powers down the detectors, and begins the long process of packing equipment. Some collaborators catch a late flight home. Others stay to begin the analysis.

Phase 4: Analysis (Month 12 – Month 24)

Signal Extraction

Marcus spends three months processing the raw data:

1. Gamma-ray tracking: Each GRETINA event is decomposed into individual gamma-ray interaction points, and the Compton scattering sequence is reconstructed to determine the gamma-ray energy and emission direction.

2. Doppler correction: Using the measured beam velocity ($\beta = 0.38 \pm 0.01$) and the gamma-ray direction from tracking, each photon energy is corrected to the projectile rest frame: $$E_\gamma^{0} = E_\gamma^{\text{lab}} \cdot \gamma(1 - \beta\cos\theta_{\text{lab}})$$

3. Particle identification gate: Only gamma rays in coincidence with a ${}^{54}\text{Ca}$ ion identified in the S800 focal plane are included. This eliminates $> 95\%$ of the background.

4. Background subtraction: A background spectrum from neighboring isotopes (${}^{53}\text{Ca}$) is scaled and subtracted to remove residual contamination.

The final Doppler-corrected spectrum shows a peak at $E_\gamma = 2043 \pm 19$ keV with 192 $\pm$ 18 counts (after background subtraction), corresponding to the $2^+_1 \to 0^+_1$ transition.

Systematic Uncertainties

The energy uncertainty is dominated by the Doppler correction: the 1% uncertainty in $\beta$ produces a $\sim 15$ keV uncertainty in the reconstructed gamma-ray energy. The GRETINA calibration contributes $\sim 3$ keV. These are added in quadrature to give a total systematic uncertainty of $\sim 15$ keV.

The number of counts has a statistical uncertainty of $\sqrt{192} \approx 14$ ($\sim 7\%$) plus a systematic uncertainty from the background subtraction ($\sim 5\%$).

The Result

The measurement: $E(2^+_1; {}^{54}\text{Ca}) = 2043 \pm 19 (\text{stat}) \pm 15 (\text{syst})$ keV.

This is compared to: - $E(2^+_1; {}^{52}\text{Ca}) = 2563$ keV (previously measured) - Shell model prediction with the GXPF1A interaction: 2020 keV (agreeing within uncertainties) - Prediction without $N = 34$ closure: $\sim 1100$ keV (ruled out at $> 5\sigma$)

The high $2^+$ energy in ${}^{54}\text{Ca}$ — though lower than in ${}^{52}\text{Ca}$ — confirms that $N = 34$ is indeed a significant subshell closure in calcium, driven by the tensor component of the nuclear force that increases the gap between the neutron $p_{3/2}$ and $p_{1/2}$ orbitals when protons fill the $f_{7/2}$ shell.

Phase 5: Publication (Month 24 – Month 30)

Elena writes the paper — a Physical Review Letters article of 5 pages. The draft circulates among the 28 co-authors for two months of revisions. It is submitted to PRL, sent to two referees (one experimentalist, one theorist), receives positive reviews with minor revisions, and is published 4 months after submission.

The abstract reads: "We report the first measurement of the $2^+_1$ state in ${}^{54}\text{Ca}$ at $2043 \pm 24$ keV, obtained via proton inelastic scattering in inverse kinematics at FRIB. The result provides direct evidence for a subshell closure at $N = 34$ in calcium isotopes, consistent with shell model predictions including the tensor force."

The paper receives 120 citations in its first two years. The data point — one number, $E(2^+_1) = 2043$ keV — joins the nuclear data tables, adding one more constraint to our understanding of the nuclear force.

The Timeline in Summary

Discussion Questions

1. Rate estimates. The proposal predicted 500 ions/s; the experiment received 300–400. Why are rate predictions uncertain? What factors could cause the discrepancy? How should a proposal be written to account for this?

2. Human factors. A nuclear physics beam time runs 24/7 for a week or more. What challenges does shift work present for data quality and decision-making? How do collaborations mitigate these risks?

3. One number, $2.5$ years. Is the effort proportionate to the result — a single energy measurement? Discuss the value of each data point in the context of building a comprehensive understanding of nuclear structure across the chart of nuclides.

4. Reproducibility. This experiment measured $\sim 200$ gamma-ray events. In a field where beam time is scarce and experiments are often one-of-a-kind, how is reproducibility ensured? What role do systematic checks, calibrations, and independent analysis play?

5. Careers. Elena is a junior professor; Marcus is a graduate student. How does a 2.5-year experiment fit into the timescales of academic careers, promotion, and doctoral programs? What pressures does this create?

Quantitative Exercises Based on This Case Study

CS30.2.1 The proposal estimates a proton inelastic scattering cross section of $\sigma = 20$ mb. The liquid hydrogen target is 30 mm thick ($\rho = 0.071$ g/cm$^3$). Calculate the target areal density in atoms/cm$^2$ and the reaction probability per beam ion.

CS30.2.2 If GRETINA has a photopeak efficiency of 6% at 2 MeV and the gamma-ray angular distribution is approximately isotropic, how many inelastic scattering events must occur to detect 200 gamma rays in the photopeak? At 350 ions/s beam rate and the reaction probability from CS30.2.1, how many hours of beam time are required?

CS30.2.3 The Doppler-corrected gamma-ray resolution depends on the velocity uncertainty. If $\beta = 0.38 \pm 0.01$ (the uncertainty arises from the target thickness — the reaction can occur at any depth), calculate the energy uncertainty of the Doppler correction at $\theta = 90°$ for a 2.04 MeV gamma ray.

CS30.2.4 Marcus's analysis pipeline processes 1,728 detector channels per event. If the data acquisition records 5,000 events per second (including background and all trigger types) and each event is 4 kB, how much raw data is produced per day? Per 7-day beam time?

Reflection: What Makes an Experimentalist

This case study illustrates that experimental nuclear physics is not merely a technical exercise but a deeply human endeavor. The experiment described here required:

• Scientific judgment: Choosing the right physics question, the right observable, and the right technique to address it.
• Technical mastery: Designing the detector setup, writing the analysis software, and understanding the systematic uncertainties.
• Project management: Coordinating 28 collaborators across 9 institutions, managing a 7-day beam time with round-the-clock shifts, and navigating the PAC process.
• Resilience: Diagnosing a 3-hour beam failure at 2:30 AM and deciding in real time how to recover lost data.
• Patience: Spending 12 months on analysis for a result that can be stated in a single sentence.

The combination of these skills — physics insight, technical ability, teamwork, and perseverance — defines the experimental nuclear physicist. Every data point on the chart of nuclides represents this investment of human effort, and every point extends our knowledge of the physical world by one hard-won increment.