Cyclotrons as Drivers for Precision Neutrino Measurements

As we enter the age of precision measurement in neutrino physics, improved flux sources are required. These must have a well-defined flavor content with energies in ranges where backgrounds are low and cross section knowledge is high. Very few sources of neutrinos can meet these requirements. However, pion/muon and isotope decay-at-rest sources qualify. The ideal drivers for decay-at-rest sources are cyclotron accelerators, which are compact and relatively inexpensive. This paper describes a scheme to produce decay-at-rest sources driven by such cyclotrons, developed within the DAEdALUS program. Examples of the value of the high precision beams for pursuing Beyond Standard Model interactions are reviewed. New results on a combined DAEdALUS--Hyper-K search for CP-violation that achieve errors on the mixing matrix parameter of 4 degrees to 12 degrees are presented.


Introduction
As we reach the 100th anniversary of the birth of Bruno Pontecorvo, neutrino physics is facing a transition. Neutrino oscillations are well established, albeit in a different form from what Pontecorvo expected [1,2]. We have a data-driven "Neutrino Standard Model," (νSM) which, despite questions about its underlying theoretical description, is remarkably predictive. Now, the neutrino community must pivot from "searches" to "precision-measurements," in which we can test the νSM. The transition requires new and better tools for these measurements and further calls for original approaches to experiments.
The mixings are described with a 3×3 matrix, commonly called the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, that connects the mass eigenstates (ν 1 , ν 2 , ν 3 ) to the flavor eigenstates (ν e , ν µ , ν τ ) : where the ranges indicate our knowledge of each of the entries [17]. Together with the mass splittings, the mixing matrix is pictorially represented in Figure 1, in which the lengths of the colored bars are proportional to the squared moduli of the matrix elements, |U αi | 2 . highly efficient. Decay-at-rest (DAR) sources satisfy these requirements and provide an opportunity for the precision measurements required for the future of neutrino physics. This paper explores cyclotrons as a relatively low cost means of producing DAR sources at or near underground laboratories. We begin by discussing the pros and cons of DAR for neutrino physics. Next, we review the history of and development of cyclotrons. We then describe the machines under development within the DAEδALUS program. The final sections of this paper provide examples of the precision science opened up by the DAR sources. We explore tests in neutrino oscillations and neutrino scattering. We also discuss the potential impact beyond particle physics.

Decay-at-rest Sources of Neutrinos
The most common source presently used in neutrino physics is the "conventional muon-neutrino beam." Such a source is produced with GeV-scale protons striking a target, resulting in pions and kaons which decay-in-flight (DIF) to produce neutrinos. The energy distribution and relevant backgrounds of a given DIF beam are dependent on the design of the beamline. For example, the characteristics of a beam are quite sensitive to the magnetic focusing and decay region geometry as well as primary/secondary hadron production and interaction physics in the target. These complications make first principles predictions of the flavor-dependent neutrino energy distributions for both signal and background from DIF beams quite difficult.
A number of techniques are available for DIF-based experiments to understand the neutrino flux. Experiments with very high interaction rates can use data to constrain the flux. For example, the MiniBooNE experiment has successfully used ν µ interactions, which come largely via pion DIF from the Booster Neutrino Beamline (BNB) at Fermi National Laboratory, to constrain the ν e backgrounds that are due to the subsequent decay of the pion's daughter muon [29]. Long baseline experiments use both a near detector and a far detector to reduce flux uncertainties. This comparison works well for charged current interactions, in which the neutrino energy can be fully reconstructed. Neutrinos which are within the acceptance of the far detector are considered in measuring the flux of the near detector. However, this approach is not possible for neutral current background events, because the neutrino energy cannot be fully reconstructed. Lastly, a wide range of event topologies are produced by conventional beams, which range in energy from hundreds of MeV (e.g. JPARC [30], BNB [29]) up to tens of GeV (e.g. NuMI [31], CNGS [32]). The cross sections for neutral current events and topologies involving pion production are not well measured and understood [33] and the pions (and their decay products) produced in the events can lead to electron-like backgrounds.
DAR neutrino sources, produced through π → µ → ν and isotope decay, offer a precision alternative to DIF beams. DAR flux sources have well defined flavor content and energy distributions, as can be seen with the flux from a π → µ → ν source shown in Figure 3. However, the neutrino energy is low compared to DIF-ranging from a few MeV to 52.8 MeV. The low energy of a DAR source is both an advantage and a disadvantage. A great advantage is that two of the best-known neutrino cross sections, each with less than 1% uncertainty, are accessible at DAR beam energies. The first is the inverse beta decay interaction (IBD,ν e + p → e + + n). This has a cross section that is well known because it is connected to neutron decay, which is measured very precisely [34]. IBD can be efficiently observed by requiring a coincidence between the prompt positron and delayed neutron capture signal [12,[18][19][20]. The coincidence signature also allows the signal to be easily distinguished from background, especially if the background is low as is the case for an underground detector. The second is neutrino-electron elastic scattering (ν + e → ν + e). This cross section is well constrained by Standard Model measurements of e + e − scattering [35]. Although this interaction lacks a coincidence signal, it is highly directional, even at DAR energies. On the other hand, the low energy neutrinos from a DAR source means that the relevant interactions have low absolute cross sections, leading to the high flux requirement. A DAR source therefore has the overall disadvantage of requiring a very intense source that can be installed at or near a detector in an underground location. Below, we will show that cyclotrons, as DAR neutrino drivers, have sufficiently high intensity and small enough size to overcome these disadvantages. DAR sources range in energy from up to a few MeV from isotope decay, where we use 8 Li decay as our example (see Figure 2), to 52.8 MeV from the π + → µ + chain (see Figure 3). The flux from isotope decay is pure in flavor, while the pion/muon decay has well-defined flavor ratios.
The pion/muon DAR beam is best produced by impinging low energy ∼ 800 MeV protons on a light target to produce a high rate of pions through the ∆-resonance. The target must be surrounded by heavy material to stop the outgoing pions before DIF. In this case, the neutrinos originate primarily from π + or µ + decay. The negatively charged pions and muons stop and capture on nuclei before they can decay to produce neutrinos [37]. The production of kaons or heavier mesons, which could produce unwanted backgrounds, is negligible if the primary proton energy is below about 1 GeV. Theν e flux can be maintained at the level ∼ 5 × 10 −4 of theν µ flux in the 20 < E ν < 52.8 MeV energy range. As a result, the source is well suited to search forν e appearance through oscillations [38], as discussed below. An isotope DAR source produces a pure electron-flavor flux through β-decay. Such a source can be produced through high-intensity, low-energy protons impinging on a beryllium target. These and subsequent interactions result in a flood of neutrons that are captured on surrounding material to produce the isotope of interest. Precision experiments are best performed using neutrinos above 3 MeV, where environmental backgrounds are low. As a result, high Q-value isotopes are favored. Below, we discuss the use of 8 Li as the decaying isotope, produced by neutron capture in a 99.99% pure 7 Li sleeve surrounding the target. This process produces a very pureν e flux of well-defined energy that can be used for scattering and neutrino disappearance experiments.
The IBD interaction has a large cross section at these MeV-scale energies, as shown in Figure 4. We note, however, that as the IBD interaction requires free protons as an interaction target, this approach is relevant for water-and scintillator-based detectors only. The IBD interaction has a distinguished history. Pontecorvo himself first suggested a search for this interaction, in a 1946 report to the National Research Council of Canada [39,40], and it was the first type of neutrino interaction observed [41]. IBD is the signal interaction now widely used in reactor experiments. The reactor IBD range is shown in Figure 4 in comparison to the DAR fluxes.

Cyclotrons as DAR Source Drivers
Cyclotrons represent ideal drivers for the DAR sources discussed above. These machines are compact and low-cost compared to most particle physics accelerators. The size, power, and cooling requirements are sufficiently modest that it is possible to install cyclotrons at underground laboratories where these sources can be paired with existing large scintillator and water detectors. Cyclotrons date back to Pontecorvo's era. However, and as in the case of neutrino physics, cyclotrons have come a long way since their origin. Modern cyclotrons are capable of producing the very high intensity flux required for modern precision neutrino measurements.
This section reviews the history of cyclotrons, with some discussion about how cyclotrons work. We consider important examples in use today and then discuss their development within the DAEδALUS project.

A brief history of cyclotrons
Unbeknownst to Ernest O. Lawrence, the cyclotron was first invented by Leo Szilard, who received a German patent for the device in 1929, but Szilard never attempted a practical realization of his idea. Lawrence's own invention stemmed from his study of a paper by Rolf Wideroe on resonant acceleration in linear structures using radio frequency (RF) voltages. Although Lawrence could not understand German, he was able to understand enough of the concept from the drawings and equations to come to his own invention.
In the Wideroe linear accelerator, a beam of ions was accelerated through a series of small gaps between hollow metal tubes, called drift tubes, connected in series to the RF-voltage generator. At any instant of time successive gaps carry a voltage of opposite sign. The voltage changes sign during the time it takes the ion to traverse the tube, allowing the ion to increase its energy. Since the nonrelativistic ions increase their velocity in passing through the gaps, each successive drift tube must be longer for the ion motion to remain in synchronism with the RF-generator. Higher and higher energy meant the length of the accelerator increased non-linearly, at least initially.
To reach energies of 1 MeV or more, Lawrence's insight was that for non-relativistic particles injected into a dipole magnetic field perpendicular to the particle velocity, the revolution frequency is independent of the particle energy. Higher energy particles travel on larger orbits, maintaining synchronous (or isochronous) motion. If the ion orbits are contained within two hollow "D"-shaped electrodes ("dees") which are connected to an RF-voltage source, one can accelerate ions to energies not possible with DC-voltage structures. Beam is injected at the center of the device and spirals outward.
Lawrence's ideas were soon realized in practice for protons by his student M. Stanley Livingston. As has been the case since the first accelerators, the Lawrence team pushed forward on two fronts: particle energy and beam intensity. The strong limitations on intensity were imposed by losses of ions on the vertical surfaces of the cyclotron's vacuum chamber. This problem was mitigated when magnetic field shims were introduced to provide a radial component to the magnetic field. The radial field increased with distance from the center of the cyclotron. The result was vertical focusing that confined the beam to the median horizontal plane. With vertical focusing, many bunches of beam particles, each with a different kinetic energy, could be accelerated in the cyclotron. The cyclotron delivered a continuous train of RF-bunches of ions.
Trying to increase the proton energy significantly beyond 10 MeV provided a different difficulty. The revolution frequency (cyclotron frequency) began to decrease due to the relativistic increase in the proton mass. Even a change as small as 1 to 2% is enough for the revolution frequency to be outside of the frequency bandwidth of the RF-generator and therefore to lose synchronism with the RF-voltage. The solution to this limitation seems simple: modulate the RF to lower values to maintain synchronism with the highest-energy beam particles. The frequency modulated cyclotron (or synchrocyclotron) allowed the Lawrence team to achieve energies as high as 340 MeV with their 184 inch cyclotron. At the University of Chicago, Enrico Fermi's team reached 450 MeV with the slightly smaller 170 inch cyclotron. In the early 1950s this machine produced copious pions via the ∆(3,3)-resonance. The age of accelerator-generated neutrinos had begun. It is noteworthy that Fermi and Szilard were also responsible for the invention of another important source of man-made neutrinos: the nuclear reactor.
Even larger machines were built in Russia. Russia's effort culminated in the giant 1 GeV proton synchrocyclotron at the Leningrad Institute for Nuclear Physics in Gatchina. This machine was built around the world's largest one-piece electromagnet with a pole diameter of 7 m and iron weight of almost 8000 tonnes. Maintaining synchronism through frequency modulation had come at a price. Only one bunch, i.e., protons of only a single energy, could be accelerated at any one time. Therefore, the Gatchina machine could only accelerate 0.2 µA, hundreds of times less than classical fixed-frequency cyclotrons.
Generating continuous trains of ion bunches at energies exceeding ∼10 MeV required yet another invention that would allow particle synchronism despite the acceleration being provided by fixed frequency RF-power. The solution published by Llewellyn Thomas in 1938 showed that conditions for resonance and focusing could be maintained if the vertical magnetic field varied with polar angle. This variation introduces an additional focusing effect and leads to scalloped rather than simple spiral orbits. If in addition to the azimuthal variation, the vertical field increases in strength at larger radii, the particle motion can be made isochronous independent of energy over the full range of operation of the cyclotron. In isochronous cyclotrons, the beam current can be increased to the mA range. Indeed, 1 mA proton beams can be produced by commercially available machines designed for radioisotope production today.
Donald Kerst introduced a further improvement of Thomas' scheme of azimuthal variation by radial sectors. Kerst suggested using spiral sectors to increase axial focusing of the beam even more through the application of the alternating gradient principle, which was by then being designed into synchrotrons. Spiral sectors are now used in almost all cyclotrons over ∼40 MeV, enormously increasing both the energies and the intensities available and thereby providing a factor of ∼1000 more intense beams for π, µ, n and neutrino production at low energies.
The energy and intensity (or current) range provided by historical and present-day cyclotrons, compared to other types of circular accelerators, is shown in Figure 5. Various types of cyclotrons are noted: FF is the Fixed Field or Classical Cyclotron; FM is the Frequency Modulation (Synchro-) Cyclotron; and AVF is the Azimuthal Varying Field Cyclotron. One can see that, at the low energies needed for DAR beams, cyclotrons are ideal drivers. Linear accelerators are also an option, but require much higher power and have much higher cost per unit energy than cyclotrons.
Research on two existing cyclotrons has provided important intellectual input for the 800 MeV machine ultimately envisioned for the DAEδALUS cyclotron program. These are the Ring Cyclotron at the Paul Scherrer Institut (PSI) in Villigen, Switzerland [43], and the Superconducting Ring Cyclotron (SRC) [44] at RIKEN, Wako, Japan, which, although designed for high-energy highly stripped heavyion beams, represents an engineering "proof-of-principle" design for a cyclotron magnet applicable for DAEδALUS. With an energy of 590 MeV and beam current of 2.4 mA, the PSI Ring Cyclotron is currently the world's most powerful accelerator in this energy range, delivering 1.4 MW of protons [43], as seen in Figure 5. The PSI complex routinely achieves 99.98% extraction efficiency, and this sets the bar for future accelerators such as those in the DAEδALUS program. The RIKEN SRC is the world's first ring cyclotron that uses superconducting magnets, and has the strongest beam bending force among the cyclotrons. The magnet design for the 800 MeV/n DAEδALUS SRC is based on RIKEN. RIKEN does not appear on Figure 5 because it is a heavy ion rather than a proton machine. As such, the current from the RIKEN machine is limited by the available shielding, and not by the machine design. RIKEN can boost the ion beam energy up to 440 MeV/nucleon for light ions and 350 MeV/nucleon for very heavy ions, such as uranium nuclei, to produce intense radioactive beams. The ring cyclotron consists of six major superconducting sector magnets with a maximum field of 3.8 T. The total stored energy is 235 MJ, and its overall dimensions are 19 m diameter, 8 m height and 8,300 tons. The magnet system assembly was completed in August 2005 and successfully reached the maximum field in November 2005. After magnetic field measurements for two months, the superconducting magnets was installed and the first beam was extracted from the SRC in December 2006.

Cyclotrons as pion/muon factories
Cyclotrons have been used to produce pions and muons for many years; what is novel about DAEδALUS is their application as drivers for DAR sources, in which the pions and muons come to rest and decay to neutrinos. In fact, two out of three of the major "meson factories" commissioned in the 1970s were cyclotron-based. These were TRIUMF, in Vancouver, BC, Canada [45], and SIN [46] (now PSI) in Villigen, Switzerland. The competing technology was LAMPF [47] (now LANSCE) at Los Alamos, which was an 800 MeV linear accelerator. These two cyclotron facilities remain at the forefront of precision pion and muon studies to this day. TRIUMF, a 500 MeV H − cyclotron, produces several hundred µA. This program has expanded to also become a world-leading laboratory for radioactive ion beams. PSI is currently the world's most powerful accelerator in this energy range with 590 MeV protons. At 1.4 MW of beam power, this cyclotron is a shining example of high-quality beams with extremely high extraction efficiency (99.98%) and low beam losses (<200 W per cyclotron vault). An example of the beautiful muon physics now being published from PSI is the precision measurement of G F from the MuLan experiment [48]. However, the PSI program with the primary beam is now evolving toward being primarily directed at production of low-energy (meV) neutrons via the spallation process for neutron scattering and diffraction studies of materials. These cyclotrons, which have been a tremendous asset to the field, inform the DAEδALUS design below.
Although neither the TRIUMF nor PSI machines have been applied to the neutrino field, DAR experiments were conducted at Los Alamos with the competing LAMPF beam [49,50]. Also, it should be noted that other low-energy synchrotrons have also hosted or are considering important neutrino deca-at-rest experiments, namely the 700 MeV ISIS machine at the Rutherford Appleton Lab in UK (KARMEN [51]) and the SNS at Oak Ridge [52].

The DAEδALUS cyclotrons
The DAEdALUS cyclotron system accelerates ions through a series of two cyclotrons. The full system is called an "accelerating module." Figure 6 shows the schematic layout of one of the H + 2 accelerating modules. The DAEδAUS injector cyclotron (DIC) captures up to 5 mA (electrical) of H + 2 and accelerates the beam to about 60 MeV/n. This beam is extracted electrostatically. This beam can then be used for a stand-alone program (IsoDAR) or for injection into the DSRC for further acceleration (DAEδALUS). This second machine consists of 8 wedge-shaped superconducting magnets and 6 RF cavities (4 of the PSI single-gap type, 2 double-gap). The stripper foil is located at the outer radius, at the trailing edge of one of the sector magnets, and the extraction channel comes out roughly along one of the valleys about 270 o away. Figure 6 also shows schematically one of the beam dumps, a graphite block with a hole shaped to correspond to the beam profile so the energy is uniformly distributed over a wide area. The graphite is surrounded by a copper, water-cooled jacket and is expected to dissipate 6 MW of beam power.
A key aspect of the design is acceleration of H + 2 . This novel choice of ion was selected by L. Calabretta [53] in response to a suggestion by Carlo Rubbia in the 1990s to use high-current, ∼ 1 GeV cyclotrons for driving thorium reactors [54]. Since most cyclotrons accelerate protons or H − , it is worth examining the pros and cons of this choice.
A drawback of H + 2 is that the higher rigidity of this ion (q/A=0.5) compared to bare protons or H − (both with |q|/A=1) requires a cyclotron of relatively large radius. However, the size of the machine is practical given the higher fields available from superconducting magnets. In fact, the RIKEN SRC is close to the field and size specifications required.
By choosing the H + 2 ion, one can use a stripping foil to cleanly extract the beam from the cyclotron. Although stripping extraction is not available to proton accelerators, it has been used extremely effectively in lower-energy cyclotrons that accelerate H − beams. The value of stripping extraction is that turn separation is no longer an issue. All ions will pass through the stripper foil, even if turn separation is not clean. Ions will be stripped from several turns, probably not more than two or three, and the protons will carry the energy associated with their turn number. The extraction channel, which will pass through the central region of the DSRC (as the protons are bent in rather than out as is the case with H − ), will need adequate momentum acceptance to transmit all the protons from the turn where they are stripped.
It is the lower binding energy (0.7 eV) of the H − ion that renders it unusable in a high-energy machine like the DSRC. This makes H − susceptible to Lorentz stripping in fields as low as 2 T and energies below about 70 MeV. The higher binding energy of the H + 2 ion (2.7 eV, at least in its ground state) renders it more stable, and able to survive to 800 MeV in the highest 6 T fields anticipated in the DSRC.
Lastly, a great advantage of H + 2 is reduction of space-charge effects. The space charge of the particle beam produces a repulsive force inside the beam, which generates detuning effects. A measure of the strength of this effect, the "generalized perveance," is defined by where q, I, m, c, γ, and β are the charge, current, rest mass, speed of light, and the relativistic parameters of the particle beam, respectively [55]. The higher the value of K, the stronger the space-charge detuning effects. According to Eq. 2, the space-charge effects for the 5 mA of H + 2 beam in the DSRC are equivalent to a 2.5 mA proton beam with the same γ. Consequently, they are similar to the space-charge effects present in the 2.4 mA proton beam currently being accelerated at PSI. Another degree of freedom to reduce space-charge effects are the choices injection energy and acceleration voltage. Given this premise, we have carried out precise beam dynamics studies, including the 3D space-charge effects (excluding the central region of the DIC), with self-consistent 3D models implemented in the code Object Oriented Parallel Accelerator Library (OPAL) [56]. The beam dynamics model is described in detail in [57,58]. For the DSRC, we have implemented a simple stripper model into OPAL in order to study the complex extraction trajectories of the stripped protons.
DAEδALUS R&D related to H + 2 acceleration has begun at a test stand at Best Cyclotrons, Inc., in Vancouver, BC, Canada. These studies employ the VIS, or Versatile Ion Source, a non-resonant ECR (Electron Cyclotron Resonance) source [59] built at the Laboratori Nazionali del Sud (LNS) in Catania, Italy, which has been shipped to Vancouver for these tests. As with any ion source, protons, H + 2 and even H + 3 will be emitted. We have begun studies of emittance, inflection into the cyclotron, capture, and acceleration of H + 2 . In summary, the DAEδALUS program has developed a plan for the production of two high power cyclotrons, one producing beam at 60 MeV/n and the other at 800 MeV/n. The former provides the injector to the latter, which then can produce pion/muon DAR neutrino fluxes. As discussed below, the injector can also be used by itself to produce isotope DAR beams.

Application of the DIC Cyclotron to Isotope Production
The remainder of this paper will describe the value of the DAEδALUS cyclotrons for basic research in neutrino physics. However, it is worth pausing to note that these cyclotrons have practical applications. This has attracted industry involvement in DAEδALUS development. Examples range from medical applications to accelerator driven systems for thorium reactors. We concentrate on the former here, and note that Ref. [60] provides a more extensive description of isotope production for medicine using the DIC cyclotron.
Even from Ponetcorvo's days, isotopes produced by cyclotrons were being investigated for suitability in medical diagnostics and therapeutics. The first direct therapeutic use of beams occurred in 1937 with trials using neutrons produced from the 60 MeV Crocker Cyclotron at Berkeley. Incidentally, this machine is still in use, primarily for proton treatments of eye tumors at UC Davis. Higher energy cyclotrons were built in the 1940's, the 800 MeV 184 inch Synchrocyclotron at Berkeley, and two similar-sized machines at Dubna and St. Petersburg. Medical treatments played prominent roles in all three, first with stereotactic microsurgery with very fine beams for pituitary ablation, then when diagnostic tools such as CT scanners became available that could carefully measure the density of material upstream of a tumor to accurately predict the stopping point of the beam, Bragg-peak therapy with alpha-particle beams was instituted at the 184 inch Synchrocyclotron.
Cyclotrons are now used extensively for therapy with proton beams, with several commercial companies marketing highly effective cyclotron-based systems for proton therapy with beam energies of around 250 MeV [61]. As beam currents needed for radiotherapy are only in the nA range, even with the losses inherent in degrading the fixed-energy beams down to match the required range in the patient, adequate beam brightness can be achieved with beams of no more than a µA of protons extracted.
Meanwhile, radioactive isotopes produced with cyclotrons of energies of 30 MeV or less have become widely used in medical diagnostics and therapy, with an ever-increasing demand as techniques are refined and results improved. Beam currents in the range of 750 µA to 2 mA are now being extracted from commercial and research isotope-producing cyclotrons; the limiting factor is often heat-dissipation in the complex targets that are needed for effective isotope production. Increased production capacity is being obtained by multiple extraction ports enabled by acceleration of H − beams that can be extracted by stripping. Sharing the total beam power between two target stations enables greater production capacity.
The DAEδALUS injector cyclotron, used for IsoDAR, will become a powerful tool for isotope production along two different directions. As a source of 60 MeV protons at 600 kW, beam powers  82 Rb, a potassium analogue This isotope is also directly used as a PET isotope for heart imaging.  [42].
are substantially higher than existing isotope machines. This could enable either significantly greater yield on a single target, should the technology be developed to use all this beam power on a target or by sharing the beam between many targets to increase the versatility of the isotope factory. As H + 2 ions are extracted from the cyclotron via a conventional septum, a narrow stripper can be placed over a portion of the beam to convert ions passing through the stripper into protons that can then be cleanly separated from the body of the beam and transferred to a production target. The remaining beam is transported to further stripping stations, each peeling of a small portion of the beam to deliver to a different target. In this way the power limits on any given target will not be exceeded, and high efficiency for use of the whole beam maintained. Examples of the isotopes which can be produced, and their applications, are shown in Table 1.
A second isotope application of the H + 2 cyclotron is that ions of the same charge-to-mass ratios can also be accelerated. Specifically, He ++ (alpha-particle) beams can be accelerated at currents limited only by the availability of such He ++ ion sources. There are many isotopes that have tremendous application potential and are limited today only by the very restricted availability of suitable highcurrent alpha beams. In fact, the first prototype cyclotron to be built for testing injection of the high-current H + 2 beams, to be built at the LNS in Catania, Italy, is being designed to be used directly following the H + 2 injection tests as a dedicated alpha-particle cyclotron for producing radiotheraputic isotopes. One example will be 211 At, which is in short supply for even long term clinical studies The DAEδALUS Superconducting Ring Cyclotron, in extending the performance of today's recordholding PSI by increasing energy from 590 to 800 MeV and a factor of five in current, becomes a member of the GeV -10 MW -class of machines. Many such machines have been designed and proposed but cost has been an impediment to their realization. To date, only one such project has progressed to the advanced R&D and construction phase: MYRRHA [62] to be sited in Mol, Belgium. These projects all fall within the ADS (Accelerator-Driven Systems) category, such as nuclear waste transmutation, driving of sub-critical thorium-based reactors, tritium production, and others.
Along with the physics possibilities previously described, the DAEδALUS cyclotrons provide new opportunities in this field by offering beams at a substantially reduced cost over the linear accelerators which until now have been viewed as the only viable technology to reach these levels of beam power in the GeV energy range. With successful development of these cyclotrons, a substantial growth in the ADS field can be anticipated, with the cost hurdle having been surpassed.

IsoDAR
IsoDAR is the first proposed step of DAEδALUS. IsoDAR represents both a novel concept of application to neutrino physics measurements and a demonstration of the 60 MeV/n injector cyclotron relevant for the larger program.
The baseline cyclotron design for IsoDAR is a 5 mA H + 2 machine that will accelerate beam to 60 MeV per nucleon. Beam would be injected at 70 keV (35 keV/amu) via a spiral inflector. An appropriate ion source already exists and its testing will be completed in Summer 2013. The current plan for IsoDAR is to locate a cyclotron accelerator underground in an experimental hall close to the KamLAND detector, in Kamioka, Japan. This is a continuous-wave source with a 90% duty cycle to allow for machine maintenance. The resulting beam will be transported a short distance up the drift at KamLAND to a target located in a room. The end of the beam dump is assumed to be 16 m from the detector center.
We continue to optimize the target for the production of 8 Li, a β-emitter and the source of the ν e for the measurement. The baseline design for the target is a cylinder of 9 Be, 20 cm long and 20 cm in diameter. This cylinder is surrounded by an additional 5 cm of D 2 O which works to both moderate neutrons and to provide target cooling. The D 2 O is then surrounded by a cylindrical sleeve of 99.99% pure 7 Li, 150 cm long and 200 cm in outer diameter. Some 8 Li is produced directly in the 9 Be target but the majority of the 8 Li is produced by the many neutrons made in the 9 Be target capturing on 7 Li. The isotopic purity of the 7 Li sleeve is needed to avoid production of tritium by neutrons on 6 Li. Further, this production cross section is several orders of magnitude larger than neutron capture on 7 Li and therefore reduces the production of 8 Li severely. The needed level of 7 Li purity for the sleeve is readily available from a number of sources as it is commonly used in the nuclear industry. A nominal running period of five years with a 90% duty cycle produces 1.29×10 23 antineutrinos from the decay of 8 Li.
When paired with a liquid scintillator detector, this isotope DAR flux opens a number of opportunities for precision neutrino measurements. This paper presents two examples. The first is a high sensitivity sterile neutrino search. The second is a search for new physics in the neutrino sector from neutrino-electron scattering. Both cases describe pairing with KamLAND, to provide specific information on rates and backgrounds. An example involving the detection of coherent neutrino scattering is also provided although such a measurement would require a dark-matter-style detector sensitive to keV-scale excitations.

Sterile neutrino searches
Searches for light sterile neutrinos with mass ∼1 eV are motivated by observed anomalies in several experiments. Intriguingly, these results come from a wide range of experiments covering neutrinos, anti-neutrinos, different flavors, and different energies. Short baseline accelerator neutrino oscillation experiments [50,63], short baseline reactor experiments [64,65], and even the radioactive source experiments, which were originally intended as calibrations for the chemical solar neutrino experiments [4,5], have all observed anomalies that can be interpreted as due to one or more sterile neutrinos. To understand these anomalies in terms of the νSM for neutrino oscillations, one or more sterile neutrinos are added to the oscillation probability calculation [66]. These extended models are referred to as "3+1", "3+2", or "3+3" neutrino models depending on the number of additional sterile neutrinos. Global fits to these data indicate that there are regions in this extended parameter space where the anomalies can be reconciled with each other as well as negative results from other experiments in and near the relevant parameter space [28]. The global fits tend to prefer models with two or three sterile neutrinos.
The exciting diversity of experiments showing these anomalous results has motivated a number of proposals to address them. Suggestions range from repeating the source experiments, to specially designed reactor antineutrino experiments, to accelerator-based ones. Many of these proposals, however, do not have sufficient sensitivity to make a definitive > 5σ statement about the existence of sterile neutrinos in all of the relevant parameter space. The experiments that are designed to make a definitive measurement are based on pion or isotope DAR sources. Notably, the full DAEδALUS complex could be used to generate a pion DAR beam for such a measurement. However, the IsoDAR concept calls for simply the DAEδALUS injector cyclotron to be used to generate an isotope DAR source. Such a complex situated next to a kiloton-scale scintillator detector such as KamLAND would enable a definitive search for sterile neutrinos by observing a deficit of antineutrinos as a function of the distance L and antineutrino energy E across the detector-the definitive signature of neutrino oscillation. This is the concept behind the IsoDAR proposal [67].
The proposed IsoDAR target is to be placed adjacent to the KamLAND detector. The antineutrinos propagate 9.5 m through a combination of rock, outer muon veto, and buffer liquid to the active scintillator volume of KamLAND. The scintillator is contained in a nylon balloon 6.5 m in radius bringing the total distance from target to detector center to 16 m. The antineutrinos are then detected via the IBD interaction. This interaction has a well known cross section with an uncertainty of 0.2% [68], and creates a distinctive coincidence signal between a prompt positron signal, E e + = Eν e − 0.78 MeV, and a delayed neutron capture giving a 2.2 MeV gamma ray within ∼200 µs.
KamLAND was designed to efficiently detect IBD. A standard analysis has a 92% efficiency for identifying IBD events [69]. In IsoDAR's nominal 5 year run, 8.2 × 10 5 IBD events are expected. The largest background comes from the 100 reactor antineutrino IBD events detected by KamLAND per year [70]. The reactor antineutrino rate is dependent on the operation of the nuclear reactors in Japan which has been significantly lower in 2012 and 2013 [71]. The sterile neutrino analysis uses an energy threshold of 3 MeV. Due to the effective background rejection efficiency provided by the IBD delayed coincidence signal, this threshold enables use of the full KamLAND fiducial volume, R<6.5 m and 897 tons, with negligible backgrounds from sources other than from the aforementioned reactor antineutrinos.
The sterile neutrino analysis makes use of neutrino oscillations's L/E signature. Therefore, the energy and vertex resolutions are essential in determining sensitivity. The KamLAND detector has a vertex reconstruction resolution of 12 cm/ E(MeV) and an energy resolution of 6.4%/ E(MeV) [69]. Example data sets for reasonable 3+1 and 3+2 sterile are shown in Figure 7 for the nominal detector parameters, summarized in Table 2. In most currently favored oscillation scenarios, the L/E signal is observable. Furthermore, separation of the various 3+N models may be possible as exemplified by Figure 7 (right).
To understand the sensitivity relative to other proposals, the IsoDAR 95% CL is compared to other electron antineutrino disappearance experiments in the two neutrino oscillation parameter space in Figure 8. In just five years of running, IsoDAR rules out the entire global 3+1 allowed region, sin 2 2θ new = 0.067 and ∆m 2 = 1 eV 2 at 20σ. This is the most definitive measurement among the proposals in the most probable parameter space of ∆m 2 between 1 − 10 eV 2 .

Precision electroweak tests of the Standard Model
In addition to the 8.2 × 10 5 IBD interactions, the IsoDAR neutrino source [72], when combined with the KamLAND detector [73], can collect the largest sample of low-energyν e -electron (ES) scatters that has been observed to date. More than 7200 ES events will be collected above a 3 MeV visible energy threshold over a 5 year run, and both the total rate and the visible energy can be measured. In the Standard Model, the ES differential cross section is given by where T ∈ 0, 2E 2 ν me+2Eν is electron recoil energy, E ν is the energy of the incomingν e , and the weak coupling constants g R and g L are given at tree level by g R = sin 2 θ W and g L = 1 2 + sin 2 θ W . Eq. 3 can also be expressed in terms of the vector and axial weak coupling constants, g V and g A , using the relations g R = 1 2 (g V − g A ) and g L = 1 2 (g V − g A ). The ES cross section can be therefore be used as a probe of the weak couplings, g V and g A , as well as sin 2 θ W , a fundamental parameter of the Standard Model as described in Ref. [78]. Although sin 2 θ W has been determined to high precision [79], there is a longstanding discrepancy [16] between the value obtained by e + e − collider experiments and the value obtained by NuTeV, a precision neutrino-quark scattering experiment [80]. Despite having lower statistics than the NuTeV, IsoDAR would measure sin 2 θ W using the purely leptonic ES interaction, which does not involve any nuclear dependence. This could therefore shed some light on the value of sin 2 θ W measured by neutrino scattering experiments.
The ES cross section is also sensitive to new physics in the neutrino sector arising from nonstandard interactions (NSIs), which are included in the theory via dimension six, four-fermion effective operators. NSIs give rise to weak coupling corrections and modify the Standard Model ES cross section given in Eq. 3 to whereg R = g R + eR ee andg L = g L + eL ee . The NSI parameters eLR eµ and eLR eτ are associated with flavor-changing-neutral currents, whereas eLR ee are called non-universal parameters. We can estimate IsoDAR's sensitivity to these parameters by fitting Eq. 4 to the measured ES cross section, assuming the Standard Model value for sin 2 θ W . In general, lepton flavor violating processes are tightly  Table 3: Total signal and background events in KamLAND with E vis between 3-12 MeV given the IsoDAR assumptions in Table 2 and the selection cuts outlined in text.
constrained so we focus only on IsoDAR's sensitivity to the two non-universal parameters eLR ee . The ES interaction used for these electroweak tests of the Standard Model is very different than the IBD interaction used for the sterile neutrino search. The IBD signal consists of a delayed coincidence of a positron and a 2.2 MeV neutron capture γ, whereas the ES signal consists of isolated events in the detector. Another difference is that at IsoDAR energies, the IBD cross section is several orders of magnitude larger than the ES cross section. In fact, if just 1% of IBD events are mis-identified as ES events, they will be the signal largest background. On the other hand, as was suggested in Ref. [81], the IBD signal can also be used to reduce the normalization uncertainty on the ES signal to about 0.7% . A final difference is that the incomingν e energy for IBD interactions in KamLAND can be inferred from the visible energy on an event-by-event basis, while the incomingν e energy for ES interactions in KamLAND cannot. Therefore, the differential ES cross section is measured in visible energy bins corresponding to the kinetic energy of the recoil electron, and the dependence on the incomingν e energy is integrated out according to the IsoDAR flux.
The backgrounds to the ES signal can be grouped into beam-related backgrounds, which are dominated by mis-identified IBD events, and non-beam backgrounds, arising from solar neutrino interactions, muon spallation, and environmental sources. We adopt a strategy similar to the one outlined in Ref. [82] to reduce the non-beam backgrounds. First, a cosmic muon veto is applied to reduce the background due to radioactive light isotopes produced in muon spallation inside the detector. This reduces the live time by 62.4%. Next, a visible energy threshold of 3 MeV is employed to reduce the background from environmental source which pile up at low energies. Finally, a fiducial radius of 5 m is used to reduce the background from external gamma rays emanating from the rock or stainless steel surrounding the detector. To reduce the beam-related backgrounds, an IBD veto is employed to reject any ES candidate that is within 2 ms of a subsequent event with visible energy > 1.8 MeV in a 6 m fiducial radius. The IBD veto is estimated to have an efficiency of 99.75%±0.02%, where the uncertainty is assumed to come from the statistical uncertainty on measuring the IBD selection efficiency with 50,000 AmBe calibration source events. Table 3 shows the expected signal and background event totals assuming a nominal 5 year IsoDAR run with a 90% duty factor. We assume that the energy spectrum of the non-beam backgrounds can be measured with 4.5 years of KamLAND data before the IsoDAR source turns on. The energy spectrum of the non-beam background, mis-identified IBD events, can be extracted from beam-on data with a dedicated delayed coincidence selection. Given these assumptions, Table 4 gives the IsoDAR sensitivity to sin 2 θ W from a combined fit to the rate and spectral "shape" of the differential ES cross section, as well as each individually. Sensitivities are also shown for the case of a 50% background reduction and for the case of a 100% background reduction.   Figure 9: IsoDAR's sensitivity to g V and g A along with allowed regions from other neutrino scattering experiments and the electroweak global best fit point taken from Ref. [16]. The IsoDAR, LSND, and TEXONO contours are all at 1σ and are all plotted in terms of g νµe V,A = g νee V,A − 1 to compare with ν µ scattering data. The ν µ e/ν µ e contour is at 90% C.L.
To compare the sensitivity of IsoDAR with that of other experiments, the fits to the ES cross section can also be done in terms of g V and g A . Fig 9 shows the IsoDAR 1σ contour in the g Vg A plane as well as contours from other experiments. IsoDAR would be the most sensitive ν e e/ν e e experiment to date and could test the consistency of ν e e/ν e e couplings with ν µ e/ν µ e couplings.
Finally, we can also estimate IsoDAR's sensitivity to the non-universal NSI parameters eL ee and eR ee , assuming the Standard Model value for sin 2 θ W = 0.238. The results are shown in Fig. 10 along with the current global allowed region [83]. In the region around eL ee and eR ee ∼ 0, the IsoDAR 90% confidence interval significantly improves the global picture.

Coherent neutrino scattering at IsoDAR
As discussed in the previous section, an intense source of neutrinos provides an immense opportunity for a number of physics measurements other than a sterile neutrino search. Along with the weak mixing angle measurement and sensitivity to non-standard neutrino interactions, such a source could allow the first detection of and subsequent high statistics sampling of coherent, neutrino-nucleus scattering events. Although the process is well predicted by the Standard Model and has a comparatively large cross section in the relevant energy region (∼10-15 MeV), neutral current coherent scattering has never been observed before as the low energy nuclear recoil signature is difficult to observe. A modest sample of a few hundred events collected with a keV-scale-sensitive dark matter style detector could improve upon existing non-standard neutrino interaction parameter sensitivities by an order of magnitude or more. A deviation from the ∼5% predicted cross section could be an indication of new physics. Furthermore, the cross section is relevant for understanding the evolution of core collapse supernovae as well as for characterizing future burst supernova neutrino events collected with terrestrial detectors.
A dark matter style detector with keV-scale sensitivity to nuclear recoil events, perhaps based on germanium crystal or single phase liquid argon/neon technology, in combination with an intense proton source such as IsoDAR could perform the physics discussed above. The technology currently exists for such a detector (with requisite passive and active shielding) to be deployed on the surface or underground. Note that the KamLAND detector is not sensitive enough for such a measurement. Figure 11 shows the expected rates in terms of neutrino energy and nuclear recoil energy for an IsoDAR source (2.58×10 22 ν e /year) in combination with a 1000 kg argon detector at a 10 m average baseline from the source with a 1 keV nuclear recoil energy threshold and 20% energy resolution. Given these assumptions, about 1200 events per year could be collected for a high statistics sampling of this event class. A first observation of the process is clearly possible with a more modest size detector as well. the light neutrino sector would be a first hint of such effects in the early universe where GUT-scale Majorana neutrinos can have CP violating decays that lead to the matter-antimatter asymmetry that we now observe. This process is called "leptogenesis." [84][85][86] The parameter δ CP is accessible through the muon-to-electron neutrino flavor oscillation probability. For oscillations in a vacuum, the probability is given by [87]: P µ→e = sin 2 θ 23 sin 2 2θ 13 sin 2 ∆ 31 ∓ sin δ sin 2θ 13 sin 2θ 23 sin 2θ 12 sin 2 ∆ 31 sin ∆ 21 + cos δ sin 2θ 13 sin 2θ 23 sin 2θ 12 sin ∆ 31 cos ∆ 31 sin ∆ 21 + cos 2 θ 23 sin 2 2θ 12 sin 2 ∆ 21 (5) where ∆ ij = ∆m 2 ij L/4E ν . In the second term, the −(+) refers to neutrino (antineutrino) operation. A critical parameter for measuring CP violation is the size of the mixing angle θ 13 , which determines the size of the first three terms in Eq. 5. Recently, several reactor neutrino disappearance experiments (Double Chooz, Daya Bay, and RENO) have made precision measurements of θ 13 giving a global average of θ 13 = 8.75 • ± 0.43 • [88]. The fact that θ 13 is now known to be fairly large makes the search for CP violation viable and a key next step in particle physics.

CP violation searches
For long baseline experiments, searches for CP violation rely on comparing neutrino and antineutrino oscillation probabilities, thus, exploiting the above change of sign in order to isolate δ CP . This type of measurement is complicated by matter effects, in which the forward scattering amplitude for neutrinos and antineutrinos differs due to the presence of electrons, rather than positrons, in matter. The matter effects result in a modification of Eq. 5 giving In this equation, a = G F Ne √ 2 and ∓ refers to neutrinos (antineutrinos). Matter effects only appear when L is large, because a ≈ (3500 km) −1 (with ρY e = 3.0 g/cm 3 ) is small. Short-baseline experiments, such as DAEδALUS and moderate baseline experiments such as T2K at L = 295 km suffer negligible matter effects. On the other hand, long baseline experiments such as NOνA and LBNE have significant matter effects [89]. The terms that are modified by the matter effects also depend on sign(∆m 2 31 ), making the corrections dependent on knowing this sign, commonly called the "mass hierarchy." For long baseline accelerator oscillation experiments, gathering sufficient antineutrino data sets is difficult due to the reduced negative pion production rate by accelerator protons and by the reduced interaction cross section of antineutrinos. The current event estimates for the LBNE experiment with a 34 kton liquid argon detector at 1300 km from the Fermilab site is shown in Figure 12 [89] for the normal mass hierarchy.
In contrast, the DAEδALUS experiment will be a search only in the antineutrino modeν µ →ν e with no matter effects, and with reduced backgrounds and systematic uncertainties, plus a unique experimental layout in which several low-cost neutrino sources are at different distances from one large detector. With an antineutrino-only beam, the oscillation probability is given by Eq. 5, and the sensitivity to CP violation comes about through the interference between ∆ 12 and ∆ 13 transitions, which have a distinctive L dependence. Specifically, DAEδALUS will search forν µ →ν e oscillations using neutrinos from three stoppedpion DAR sources, which interact in a single large 200 to 500 kton Gd-doped water Cherenkov or a large 50 kton liquid scintillator detector. The spectrum ofν µ 's that can oscillate intoν e 's is shown in Figure 4. The detection of the electron antineutrinos is done through the IBD process where the outgoing positron is once again required to have a delayed coincidence with a neutron capture on Gd for the water detector or on hydrogen for the scintillator detector. This process has a high cross section at ≈ 50 MeV, but requires either Gd doping for a water detector or a scintillator detector to detect the outgoing neutron and separate the IBD events from the preponderance of charged-current ν e events.
The accelerators will be positioned at 1.5, 8, and 20 km from the large detector as shown in Figure 13 and are all above ground to reduce the installation and running complexity. Each accelerator provides different physics data for the CP violation search. The 1.5-km accelerator allows measurement of the beam-on backgrounds and the normalization. The 8 km site is at an oscillation wavelength of about π/4 at 50 MeV and the 20 km site is at oscillation maximum for this energy. Each site will be run for 20% of the time so that the events from a given source distance can be identified by their time-stamp with respect to this running cycle. This will leave 40% of the time for beam-off running to measure the non-beam backgrounds and provide other physics data. The baseline plan is for a ten year run with 1 MW, 2 MW, and 5 MW neutrino sources at the 1.5, 8, and 20 km sites, respectively.
Combining the data from the three accelerators helps to minimize the systematic uncertainties associated with the beam and detector and leads to a highly sensitive search for CP violation. The shape of the DAR flux with energy is known to high precision and is common among the various distances; thus shape comparisons will have small uncertainties. The interaction and detector systematic errors are low since all events are detected in a single detector. The fiducial volume error on the IBD events is also small due to the extreme volume-to-surface-area ratio of the ultra-large detector. Therefore, the main errors for the measurements are related to the statistics of the data and to normalization uncertainties. The normalization uncertainties are dominated by the neutron  tagging efficiency, assumed to be 0.5%, and the antineutrino flux uncertainties that are constrained as described next.
The DAEδALUS CP violation analysis follows three steps. First, the absolute normalization of the flux from the near accelerator is measured using the >21,000 neutrino-electron scatters from that source in the detector, for which the cross section is known to 1%. The relative flux normalization between the sources is then determined using the comparative rates of charged current ν e -oxygen (or ν e -carbon) interactions in the the detector. Since this is a relative measurement, the cross section uncertainty does not come in but the high statistics is important. Once the normalizations of the accelerators are known, then the IBD data can be fit to extract the CP -violating parameter δ CP . The fit needs to include all the above systematic uncertainties along with the physics parameter uncertainties associated with, for example, the knowledge of sin 2 2θ 13 and sin 2 θ 23 , which are assumed to be known with an error of ±0.005 and ±0.01, respectively.
DAEδALUS must be paired with water or scintillator detectors that have free proton targets. The original case was developed for a 300 kt Gd doped water detector at Homestake, in coordination with LBNE [91]. Subsequently, DAEδALUS was incorporated into a programa with the 50 kt LENA detector [92] (called "DAEδALUS@LENA"). This paper introduces a new study, where DAEδALUS is paired with the Gd-doped 560 kt Hyper-K [93] ("DAEδALUS@Hyper-K"). This results in inprecedented sensitivity to CP violation when "DAEδALUS@Hyper-K" data is combined with data from Hyper-K running with the 750 kW JPARC beam. ("DAEδALUS/JPARC@Hyper-K"). In this scenario, JPARC provides a pure ν µ flux, rather than running in neutrino and antineutrino mode. This plays to the strength of the JPARC conventional beam, while DAEδALUS provides a high statisticsν µ flux with no ν µ contamination. A summary of the assumptions for the various configuration scenarios is provided in Table 5.
CP violation sensitivities have been estimated for 10 year baseline data sets for all the configurations given in Table 5 using a ∆χ 2 fit with pull parameters for each of the systematic uncertainties. For the DAEδALUS configurations, data from all three neutrino sources are included along with the neutrino-electron and ν e -oxygen (or ν e -carbon) normalization samples. As an example, Table 6 and Figure 14 presents a summary of the events by category for the DAEδALUS@Hyper-K configuration. The precision for measuring the δ CP parameter in the DAEδALUS@Hyper-K configuration is given in Table 7 for sin 2 2θ 13 = 0.10 [88], both for the total and statistical-only uncertainty. The distribution of the uncertainty as a function of δ CP is shown in Figure 15. From these estimates, it is clear that, even with the large Hyper-K detector, the measurement is dominated by statistical uncertainty. Also shown in the table are estimates of the measurement uncertainties for the proposed Hyper-K [93] and LBNE [94] experiments for ten year runs with the proposed upgraded beam intensities (0.75 MW for HyperK and an average 0.85 MW for LBNE). Depending on the true value, DAEδALUS has comparable sensitivity for measuring δ CP but has very different systematic uncertainties. Thus, DAEδALUS could provide key information that can be used in conjunction with the other experiments to reduce the global measurement uncertainty.   Table 7: DAEδALUS@Hyper-K 1σ measurement uncertainty (in degrees) on δ CP for sin 2 θ 13 = 0.10 assuming the baseline 10 year data sample with a 560 kton Gd-doped water detector. (Statistical only errors are shown in parentheses.) Also shown is an estimate for the JPARC@Hyper-K sensitivity for a 560 kton water detector run for 7.5 MWyrs (3 years ν and 7 yearsν) assuming 5% systematic errors and the LBNE experiment with a 35 kton liquid argon detector run for 8.5 MWyrs (5 years ν and 5 yearsν).
The DAEδALUS high-statistics antineutrino data can be combined with a neutrino-only long baseline measurement to provide improved sensitivity for measuring δ CP . One possibility is a ten year neutrino-only run of the JPARC@Hyper-K configuration combined with a ten year DAEδALUS@Hyper-K exposure. The complementarity of the two experiments allows for a very precise search for CP violation with uncertainties estimated to be around 5 • . For this discussion we make the same assumptions are were used for Table 6: a 560 kton Gd-doped water detector, sin 2 2θ 13 = 0.1, and θ 23 = 49 • [88].
The power of the combined run is shown in Figure 15. The top plot shows the expectation for the two experiments individually. Nominal JPARC@Hyper-K running assumes three years of running in neutrino mode. This data set would yield the uncertainty indicates by the green diamonds. This would be followed by seven years of running in antineutrino mode. This data set, alone, results in the curve indicated by the green × symbols. One clearly sees that the strength of JPARC@Hyper-K is in neutrino running, as one would expect from a conventional neutrino beam. Combining these two data sets gives the green solid curve with triangles. DAEδALUS@Hyper-K alone, with a 10 year run, results in the solid red curve. One can see that DAEδALUS@Hyper-K has a similar shape to the JPARC@Hyper-K antineutrino running, where the differences come from the additional purity of the ν e flux and the lack of a mass hierarchy effect in the antineutrino data sample from DAEδALUS. We are proposing to combine a JPARC@Hyper-K run in neutrino-mode only with the DAEδALUS@Hyper-K pure antineutrino data set. The result is shown on Figure 15, bottom, by the black curve. The individual contributions of the experiments are also shown. One can see the complementarity, where DAEδALUS@Hyper-K provides the strong reach for δ cp < 0 and JPARC@Hyper-K provides the strong reach for δ cp > 0 Finally, Figure 16 shows the cross comparison of the experimental configurations shown in Table 5. The the combination DAEδALUS/JPARC(nu-only)@Hyper-K configuration is compared to the two DAEδALUS only configurations using the LENA and Hyper-K detectors in Figure 16 (top) and to JPARC@Hyper-K and LBNE in Figure 16 (bottom). From this figure, it is clear that the DAEδALUS/JPARC(nu-only)@Hyper-K configuration has impressive sensitivity to δ CP with a significantly smaller measurement error as compared to any of the other scenarios. Figure 17 shows a comparison of the δ CP regions where an experiment can discover CP violation by excluding the δ CP = 0 • or 180 • at 3σ or 5σ. Again, the DAEδALUS/JPARC(nu-only)@Hyper-K experiment clearly has substantially better coverage.

Other physics with DAEδALUS
A number of other potential physics opportunities complement the main goal of a measurement of the neutrino CP violating phase with DAEδALUS. These experiments can be located near any of the three cyclotron locations. We provide three examples here.
A large water detector, used for a sin 2 θ W measurement using ν e -electron scattering [95], would be complementary to theν e -electron search in two ways. First, it explores differences in neutrinos versus antineutrinos that can be introduced by new physics processes. Second, with the ultra-high statistics of the DAEδALUS machines, an energy-dependent analysis, rather than a rate analysis, becomes possible.
A search for ν µ → ν e appearance at high ∆m 2 , hence short-baseline, can be performed to address the LSND [50] and MiniBooNE [63] signals, if the DAEδALUS configuration uses a large scintillator detector such as LENA [96]. In this case, the cyclotron must be located underground, within tens of meters of the detector. Like the ν e disappearance search of IsoDAR, this study searches for the oscillation wave across the detector. Thus it would be powerful confirmation of a sterile-neutrinorelated oscillation signal as the source of this high ∆m 2 anomaly. The ν e flux can also be used for oscillation studies via the disappearance channel [96].
Lastly, a discovery of coherent neutrino scattering is also possible at a DAEδALUS cyclotron. Notably, a cyclotron can provide the source of neutrinos for a coherent discovery at a deep underground detector 1.5 km away from the source [97] with only a small effect on such a detector's dark matter search exposure. Furthermore, hundreds of coherent events can be collected with a dark-matter-style detector close to such a source for non-standard neutrino interaction sensitivity. A sensitive, unique neutral-current-based sterile neutrino search using coherent events can also be accomplished [36].

Conclusion
At the 100th anniversary of Pontecorvo's birth, neutrino physics is entering a new "precision era." To achieve our goals for the next 100 years, improved flux sources are needed. Decay-at-rest sources, driven by cyclotron accelerators, offer neutrino beams of well-defined flavor content and with energies in ranges  Table 5 for the description of each configuration.  Table 5 for the description of each configuration.
where backgrounds are low and knowledge of the cross section is high. This paper describes schemes to produce isotope and pion/muon decay-at-rest sources, developed as a part of the DAEδALUS program. This paper has provided examples of the value of the high precision beams for pursuing new physics. In particular, new results on a combined DAEδALUS-Hyper-K search for CP -violation are presented. This study shows that errors on the mixing matrix parameter ranging from 4% to 12% are achievable. While this result is a centerpiece of the program, DAEδALUS, and its early phase program, IsoDAR, allow for a wide range of important measurements and searches. Many examples have been presented here, focussing primarily on searches for beyond Standard Model Physics through oscillations and non-standard interactions. This establishment of this rich new program is an great way for today's neutrino physicists-the intellectual descendants of Pontecorvo-to celebrate the anniversary of his birth.