pith. machine review for the scientific record. sign in

arxiv: 2604.10447 · v2 · submitted 2026-04-12 · ⚛️ physics.ins-det

Recognition: unknown

The SpinQuest Microwave System for Dynamic Nuclear Polarization

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:27 UTC · model grok-4.3

classification ⚛️ physics.ins-det
keywords dynamic nuclear polarizationmicrowave controlreinforcement learningmonte carlo simulationpolarized targetsspinquest experimentextended interaction oscillatorai automation
0
0 comments X

The pith

A Monte Carlo digital twin trains AI controllers to autonomously tune a 140 GHz microwave system and maintain optimal dynamic nuclear polarization despite radiation-induced drifts.

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

The paper describes the integration of an extended interaction oscillator microwave source with continuous-wave NMR feedback and cryogenic monitoring in the SpinQuest polarized ammonia target. It builds an automation layer around a Monte Carlo simulation that captures DNP rate equations, frequency dependence, dose drifts, beam depolarization, and measurement noise. This simulation serves as a testbed for heuristic, reinforcement learning, and unsupervised control algorithms that adjust cavity tuning and anode voltage in real time. A sympathetic reader cares because the resulting remote, adaptive operation reduces manual exposure to radiation while sustaining higher average polarization for spin-structure measurements.

Core claim

The authors claim that their Monte Carlo digital twin of the DNP process, incorporating rate-equation dynamics, frequency-dependent steady-state behavior, dose-induced frequency drift, beam-induced depolarization, and realistic NMR noise, can be used to design and validate control strategies. These strategies enable autonomous frequency tuning through motorized cavity adjustment, combined with anode-voltage modulation for power control, which together improve polarization ramp-up speed and maintain near-optimal values under evolving high-radiation conditions.

What carries the argument

The Monte Carlo digital twin of the DNP process, which models the full set of rate equations, frequency dependence, drifts, depolarization, and noise to benchmark and train feedback and reinforcement-learning controllers for microwave frequency and power.

If this is right

  • Remote operation of the microwave system through automated tuning reduces personnel exposure to radiation.
  • Simultaneous control of cavity frequency and anode voltage avoids power nonuniformities and better matches the broad Larmor distribution in irradiated targets.
  • The AI methods increase polarization ramp-up efficiency and sustain near-optimal values as radiation conditions change.
  • The overall framework supplies a template that can be scaled to other polarized-target experiments and cryogenic high-field applications.

Where Pith is reading between the lines

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

  • If the digital twin proves accurate, the same simulation-plus-RL pipeline could be adapted to other DNP setups where manual frequency adjustment is impractical.
  • Real-time data from the physical system could be fed back to update the twin, creating a continuously improving model.
  • The approach might shorten experiment downtime by predicting optimal settings before radiation damage accumulates.
  • One could test generalization by applying the controller to a different target material or field strength and checking whether polarization performance remains close to simulation.

Load-bearing premise

The Monte Carlo digital twin accurately reproduces the real DNP rate equations, frequency-dependent steady-state behavior, dose-induced frequency drift, beam-induced depolarization, and NMR noise so that control strategies tested in simulation transfer to the physical system.

What would settle it

Deploy the trained reinforcement-learning controller on the physical SpinQuest EIO system and measure whether the resulting target polarization under actual beam conditions matches the polarization levels and stability predicted by the digital twin.

Figures

Figures reproduced from arXiv: 2604.10447 by Dustin Keller, Jordan D. Roberts, Vibodha Bandara.

Figure 1
Figure 1. Figure 1: EIO microwave generator mounted on the moving table of the target [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: EIO mounted on the target lifter, with microwave output coupled to [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Motorized actuator assembly connected to the EIO for precise [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Measured EIO output frequency as a function of tuning-shaft position. [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: LabVIEW interface for microwave-frequency control and monitoring. The main strip charts display target polarization, microwave frequency, and [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: LabVIEW interface for the microwave-frequency-control DNP simulator. The Gaussian-like frequency distributions on the left, shown in red for [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Simulation of NH3 from the polarization-data generator, showing polarization ramp-up with a continuous electron beam turned on at about 70 min. The response exhibits a sudden drop from beam heating followed by gradual decay due to radiation damage. underlying polarization Pn(tk) and returns Pk = Pn(tk) + δk, where δk is drawn from a Gaussian distribution with zero mean and configurable standard deviation. … view at source ↗
Figure 9
Figure 9. Figure 9: Positive polarization achieved in NH3 during the final commissioning run using the basic heuristic automation. Overall, the automated controller performed significantly better than manual control because it continuously evaluates the polarization slope over successive sets of measurements. Performance was characterized using the standard NH3 sim￾ulation. In maintenance mode, the average error rate was appr… view at source ↗
read the original abstract

The SpinQuest experiment at Fermilab employs a dynamically polarized solid ammonia target to probe the spin structure of the proton, requiring stable, optimized microwave-driven Dynamic Nuclear Polarization (DNP) under high radiation conditions. We present the design, operation, and automation of a 140 GHz microwave system based on an extended interaction oscillator (EIO), integrated with real-time polarization feedback from a continuous-wave NMR system and cryogenic diagnostics. The system enables fine frequency control through motorized cavity tuning and is operated remotely to mitigate radiation exposure. To continuously optimize target polarization, we develop an automation framework supported by a Monte Carlo (digital twin) of the DNP process. The simulation incorporates rate-equation dynamics, frequency-dependent steady-state behavior, dose-induced frequency drift, beam-induced depolarization, and realistic NMR noise. This framework is used to design and benchmark control strategies, including a heuristic feedback algorithm, reinforcement learning (RL), and unsupervised RL approaches. These methods enable autonomous frequency tuning, improve ramp-up efficiency, and maintain near-optimal polarization under evolving conditions. We also demonstrate integration of EIO power-supply control into the feedback loop via anode voltage modulation, providing an additional degree of freedom for simultaneous control of microwave frequency and RF power. This combined control of cavity tuning and anode voltage allows the system to avoid frequency-dependent power nonuniformities and to better match broad Larmor distributions in irradiated targets. The results establish a scalable framework for AI-driven control of complex microwave systems in polarized-target experiments, with implications for future spin-physics measurements and other cryogenic, high-field applications.

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

2 major / 1 minor

Summary. The manuscript presents the design, remote operation, and automation of the 140 GHz EIO microwave system for the SpinQuest DNP target at Fermilab. It describes integration with continuous-wave NMR polarization feedback and cryogenic diagnostics, motorized cavity tuning, and an automation framework built around a Monte Carlo digital twin that incorporates rate-equation dynamics, frequency-dependent steady-state polarization, dose-induced frequency drift, beam depolarization, and NMR noise. Control strategies (heuristic feedback, RL, unsupervised RL) and anode-voltage modulation are developed and benchmarked inside this twin to achieve autonomous frequency tuning, faster polarization ramp-up, and sustained near-optimal polarization; the work concludes that the combined framework offers a scalable approach for AI-driven microwave control in polarized-target experiments.

Significance. If the digital twin is shown to transfer to hardware, the paper supplies a concrete, reproducible engineering template for closed-loop optimization of complex, high-field microwave systems under radiation, directly relevant to future spin-structure measurements and other cryogenic DNP applications.

major comments (2)
  1. [Abstract and Monte Carlo model section] Abstract and §3 (Monte Carlo model description): the central claims of improved ramp-up efficiency and near-optimal polarization maintenance are supported only by benchmarks performed inside the Monte Carlo twin. No measured polarization curves, ramp-up times, or stability metrics from the physical 140 GHz EIO system during SpinQuest beam operations are reported, nor is any quantitative comparison (e.g., RMS deviation or correlation coefficient) between simulated and experimental data provided. This leaves the transferability of the RL and heuristic policies untested and makes the scalability assertion rest on an unverified modeling assumption.
  2. [Abstract and control strategy section] Abstract and §4 (control strategy benchmarks): the statements that the methods “enable autonomous frequency tuning” and “maintain near-optimal polarization under evolving conditions” are demonstrated exclusively in simulation. Without at least one closed-loop run on the real hardware (or a clear statement that such runs are outside the present scope), the load-bearing claim that the framework improves performance in the actual irradiated target cannot be evaluated.
minor comments (1)
  1. [Introduction or Conclusions] The manuscript would benefit from an explicit statement in the introduction or conclusions clarifying whether the reported performance numbers are simulation-only or include any hardware validation data.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. We address each major comment below, acknowledging that the performance claims rest on simulation benchmarks within the Monte Carlo digital twin. We will revise the manuscript to clarify the scope of these results.

read point-by-point responses
  1. Referee: [Abstract and Monte Carlo model section] Abstract and §3 (Monte Carlo model description): the central claims of improved ramp-up efficiency and near-optimal polarization maintenance are supported only by benchmarks performed inside the Monte Carlo twin. No measured polarization curves, ramp-up times, or stability metrics from the physical 140 GHz EIO system during SpinQuest beam operations are reported, nor is any quantitative comparison (e.g., RMS deviation or correlation coefficient) between simulated and experimental data provided. This leaves the transferability of the RL and heuristic policies untested and makes the scalability assertion rest on an unverified modeling assumption.

    Authors: We agree that the reported improvements in ramp-up efficiency and polarization maintenance are demonstrated exclusively via benchmarks inside the Monte Carlo digital twin, with no experimental polarization curves, ramp-up times, stability metrics, or quantitative sim-to-experiment comparisons provided from the physical 140 GHz EIO system. The manuscript's focus is the design of the microwave system, development of the digital twin incorporating the listed physics, and benchmarking of control strategies within that simulation. We will revise the abstract and §3 to state explicitly that these metrics are simulation results, to discuss key modeling assumptions (rate equations, dose-induced drift, beam depolarization, NMR noise), and to note that hardware validation under beam conditions is planned future work. This will qualify the scalability assertion accordingly. revision: yes

  2. Referee: [Abstract and control strategy section] Abstract and §4 (control strategy benchmarks): the statements that the methods “enable autonomous frequency tuning” and “maintain near-optimal polarization under evolving conditions” are demonstrated exclusively in simulation. Without at least one closed-loop run on the real hardware (or a clear statement that such runs are outside the present scope), the load-bearing claim that the framework improves performance in the actual irradiated target cannot be evaluated.

    Authors: We agree that the demonstrations of autonomous frequency tuning and sustained near-optimal polarization are shown only through simulation benchmarks in the digital twin, with no closed-loop runs on the physical hardware reported. The present work centers on the automation framework and its evaluation inside the twin. We will revise the abstract and §4 to specify that these capabilities are simulation results and to include an explicit statement that real-hardware closed-loop implementation and testing lie outside the current scope and are planned for subsequent studies. This will ensure the claims are appropriately scoped. revision: yes

Circularity Check

0 steps flagged

No circularity; engineering integration with simulation benchmarking only

full rationale

The paper describes hardware design, remote operation, and a Monte Carlo digital twin used to benchmark control algorithms (heuristic, RL, unsupervised RL) plus anode-voltage modulation. No derivation chain is presented that claims to predict new quantities from first principles; the simulation is explicitly a design tool whose outputs are not asserted as independent experimental results. No self-citations, fitted parameters renamed as predictions, or ansatzes smuggled via prior work appear in the load-bearing steps. The central claim is an engineering framework whose validity rests on future physical transfer, not on any internal reduction to its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents identification of specific free parameters or axioms; the simulation is stated to incorporate rate-equation dynamics and realistic NMR noise but no explicit fitted constants or unstated assumptions are listed.

pith-pipeline@v0.9.0 · 5577 in / 1110 out tokens · 41333 ms · 2026-05-10T16:27:52.946468+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

38 extracted references · 14 canonical work pages

  1. [1]

    E1039 fnal proposal,

    A. Klein, D. Keller, K. Liuet al., “E1039 fnal proposal,”SEAQUEST Document 1720-v3, 2016

  2. [2]

    Dynamic nuclear polarization,

    M. Abraham, M. A. H. McCausland, and F. N. H. Robinson, “Dynamic nuclear polarization,”Phys. Rev. Lett., vol. 2, pp. 449–451, Jun 1959. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett.2.449

  3. [3]

    Principles of dynamic nuclear polar- isation,

    A. Abragam and M. Goldman, “Principles of dynamic nuclear polar- isation,”Reports on Progress in Physics, vol. 41, pp. 395–467, 1978, printed in Great Britain

  4. [4]

    A new principle in nuclear resonance. applications,

    A. Abragam and C. S. . G. sur Yvette (France), “A new principle in nuclear resonance. applications,”International Atomic Energy Agency,

  5. [5]

    Available: https://inis.iaea.org/records/503js-9bm91

    [Online]. Available: https://inis.iaea.org/records/503js-9bm91

  6. [6]

    Possible new effects in superconductive tunnelling,

    A. Abragam, M. Borghini, P. Catillon, J. Coustham, P. Roubeau, and J. Thirion, “Diffusion de protons polarises de 20 mev par une cible de protons polarises et mesure preliminaire du parametre cnn,” Physics Letters, vol. 2, no. 7, pp. 310–311, 1962. [Online]. Available: https://www.sciencedirect.com/science/article/pii/0031916362901221

  7. [7]

    Faddeev and L

    M. Borghini, “Nuclear spin relaxation and dynamic polarization versus electron spin-spin relaxation,”Physics Letters A, vol. 26, no. 6, pp. 242–244, 1968. [Online]. Available: https://www.sciencedirect.com/ science/article/pii/0375960168906269

  8. [8]

    Spin-temperature model of nuclear dynamic polarization using free radicals,

    m. Borghini, “Spin-temperature model of nuclear dynamic polarization using free radicals,”Phys. Rev. Lett., vol. 20, no. 9, pp. 419–421, Feb

  9. [9]

    Available: https://link.aps.org/doi/10.1103/PhysRevLett

    [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett. 20.419

  10. [10]

    Solid polarized targets for nuclear and particle physics experiments,

    D. Crabb and W. Meyer, “Solid polarized targets for nuclear and particle physics experiments,”Annual Review of Nuclear and Particle Science, vol. 47, no. 1, pp. 67–109, 1997

  11. [11]

    T. O. Niinikoski,The Physics of Polarized Targets. Cambridge University Press, 2020

  12. [12]

    Polarization of nuclei in metals,

    A. W. Overhauser, “Polarization of nuclei in metals,”Phys. Rev., vol. 92, pp. 411–415, Oct 1953. [Online]. Available: https: //link.aps.org/doi/10.1103/PhysRev.92.411

  13. [13]

    A polarized proton and deuteron target for precise polarization measurement,

    S. Hiramatsu, S. Isagawa, S. Ishimoto, A. Masaike, and K. Morimoto, “A polarized proton and deuteron target for precise polarization measurement,”Japanese Journal of Applied Physics, vol. 19, no. 1, p. 161, jan 1980. [Online]. Available: https://doi.org/10.1143/JJAP.19.161

  14. [14]

    Acceleration of polarized protons to 22 gev/c and the measurement of spin-spin effects inp ↑+p↑→p+p,

    F. Z. Khiariet al., “Acceleration of polarized protons to 22 gev/c and the measurement of spin-spin effects inp ↑+p↑→p+p,” Phys. Rev. D, vol. 39, pp. 45–85, Jan 1989. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevD.39.45

  15. [15]

    The millimeter-wave extended interaction oscillator,

    W. R. Day and J. A. Noland, “The millimeter-wave extended interaction oscillator,”Proceedings of the IEEE, vol. 54, no. 4, pp. 539–543, Apr

  16. [16]

    Available: https://doi.org/10.1109/PROC.1966.4771

    [Online]. Available: https://doi.org/10.1109/PROC.1966.4771

  17. [17]

    D. G. Crabb, C. B. Higley, A. D. Krisch, R. S. Raymond, T. Roser, J. A. Stewart, and G. R. Court, “Observation of a 96

  18. [18]

    The virginia/basel/slac polarized target: operation and performance during experiment e143 at slac,

    D. Crabb and D. Day, “The virginia/basel/slac polarized target: operation and performance during experiment e143 at slac,”Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 356, no. 1, pp. 9–19, 1995, proceedings of the Seventh International Workshop on Polarized Target Mate...

  19. [19]

    A study of lithium deuteride as a material for a polarized target,

    S. Bültmannet al., “A study of lithium deuteride as a material for a polarized target,”Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 425, no. 1, pp. 23–36, 1999. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0168900298013412

  20. [20]

    A polarized target for the clas detector,

    C. Keithet al., “A polarized target for the clas detector,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 501, no. 2, pp. 327–339, 2003. [Online]. Available: https: //www.sciencedirect.com/science/article/pii/S0168900203004297

  21. [21]

    Dynamically polarized target for the g2p and gep experiments at jefferson lab,

    J. Pierceet al., “Dynamically polarized target for the g2p and gep experiments at jefferson lab,”Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 738, pp. 54–60, 2014. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0168900213016999

  22. [22]

    Operation of a Longitudinally Polarized Solid Nuclear Target in CLAS12,

    P. Pandey, J. Brock, T. Kageya, C. Keith, S. Kuhn, V . Lagerquist, J. Maxwell, and X. Wei, “Operation of a Longitudinally Polarized Solid Nuclear Target in CLAS12,”PoS, vol. SPIN2023, p. 211, 2024

  23. [23]

    Measurement of the Longitudinal Spin Asymmetry of the Deuteron in the Resonance Region,

    K. L. Kovacs, “Measurement of the Longitudinal Spin Asymmetry of the Deuteron in the Resonance Region,” Ph.D. dissertation, Virginia U., 2010

  24. [24]

    Deeply virtual compton scattering on the neutron with a longitudinally polarized deuteron target,

    S. Niccolai, G. Charles, R. Dupré, M. Guidal, D. Marchand, C. Munoz Camacho, E. V outier, A. Biselli, C. Keith, H. Avakian, V . Burkert, F. Girod, L. Elouadrhiri, V . Kubarovsky, K. Park, P. Rossi, S. Stepanyan, M. Ungaro, S. Pisano, V . Lucherini, M. Mirazita, D. Sokhan, B. McKinnon, G. Murdoch, M. Battaglieri, A. Celentano, R. De Vita, E. Fanchini, M. O...

  25. [25]

    A polarized target measurement of the electric form factor of the neutron at jlab,

    N. Savvinov, “A polarized target measurement of the electric form factor of the neutron at jlab,” inConference: A polarized target measurement of the electric form factor of the neutron at Jlab. Thomas Jefferson National Accelerator Facility, Newport News, V A, 09 2004. [Online]. Available: https://www.osti.gov/biblio/883635

  26. [26]

    Sane experiment,

    H. Baghdasaryan and S. Collaboration, “Sane experiment,”AIP Conference Proceedings, vol. 1423, no. 1, pp. 214–218, 02 2012. [Online]. Available: https://doi.org/10.1063/1.3688805

  27. [27]

    A polarized target for the clas detector,

    C. D. Keith, M. Anghinolfi, M. Battaglieri, D. Branford, S. Bultmann, V . D. Burkert, S. A. Comer, D. G. Crabb, R. D. Vita, G. Dodgeet al., “A polarized target for the clas detector,”Elsevier Science, 09 2002. [Online]. Available: https://www.osti.gov/biblio/801753

  28. [28]

    JLab experiment 01-006: Resonances-prime spin structure,

    O. A. Rondon-Aramayo, “JLab experiment 01-006: Resonances-prime spin structure,”Fizika B, vol. 13, pp. 57–64, 2004

  29. [29]

    Iniewski,Radiation Effects in Semiconductors

    K. Iniewski,Radiation Effects in Semiconductors. CRC Press, 2011

  30. [30]

    Microwave ORC: HV Short Cable,

    Fermi National Accelerator Laboratory, “Microwave ORC: HV Short Cable,” Fermi National Accelerator Laboratory, Batavia, Illinois, USA, Tech. Rep. ORC-1968, 2021, Operational Readiness Clearance (ORC) document

  31. [31]

    Extended Interaction Klystron Technology at Millimeter and Sub-Millimeter Wavelengths,

    B. Steer, A. Roitman, P. Horoyski, M. Hyttinen, R. Dobbs, and D. Berry, “Extended Interaction Klystron Technology at Millimeter and Sub-Millimeter Wavelengths,” Communications & Power Industries Canada Inc., Tech. Rep. [Online]. Available: https://www.cpii.com/docs/related/40/EIK%20Technology% 20at%20MMW%20%26%20SubMMW%20Wavelengths.pdf

  32. [32]

    CW Extended In- teraction Klystron Oscillator Final Test Results: Model VKT2438P2, Serial E0480B6,

    Communications & Power Industries Canada Inc., “CW Extended In- teraction Klystron Oscillator Final Test Results: Model VKT2438P2, Serial E0480B6,” Communications & Power Industries Canada Inc., Tech. Rep., May 2016, manufacturer final test sheet, dated May 9, 2016

  33. [33]

    ——,VKT2438P6M Extended Interaction Oscillator (EIO): Instruction Manual, Communications & Power Industries Canada Inc., Georgetown, Ontario, Canada, n.d., Manufacturer’s manual

  34. [34]

    Varian Canada Inc.,VPW 2838A2 EIO Power Supply: Instruction Manual, Varian Canada Inc., Georgetown, Ontario, Canada, n.d., Man- ufacturer’s manual

  35. [35]

    Solid-effect rate for dynamic nuclear polarization of spin-1/2 and spin-1 nuclei,

    L. JiZhi, “Solid-effect rate for dynamic nuclear polarization of spin-1/2 and spin-1 nuclei,”Communications in Theoretical Physics, vol. 31, no. 4, p. 619, jun 1999. [Online]. Available: https: //doi.org/10.1088/0253-6102/31/4/619

  36. [36]

    Dynamic orientation of nuclei by forbidden transitions in paramagnetic resonance,

    C. D. Jeffries, “Dynamic orientation of nuclei by forbidden transitions in paramagnetic resonance,”Phys. Rev., vol. 117, pp. 1056–1069, Feb 1960. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRev.117.1056

  37. [37]

    An asme-compliant helium-4 evaporation refrigerator for the spinquest experiment,

    J. D. Roberts, V . Bandara, K. Nakano, and D. Keller, “An asme-compliant helium-4 evaporation refrigerator for the spinquest experiment,”Unpublished, 2025. [Online]. Available: https://arxiv.org/ abs/2511.09689

  38. [38]

    Dynamic polarization of nuclei by electron-nuclear dipolar coupling in crystals,

    O. S. Leifson and C. D. Jeffries, “Dynamic polarization of nuclei by electron-nuclear dipolar coupling in crystals,”Phys. Rev., vol. 122, pp. 1781–1795, Jun 1961. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRev.122.1781