• ### Cosmogenic neutron production at Daya Bay(1711.00588)

March 23, 2018 hep-ex, physics.ins-det
Neutrons produced by cosmic ray muons are an important background for underground experiments studying neutrino oscillations, neutrinoless double beta decay, dark matter, and other rare-event signals. A measurement of the neutron yield in the three different experimental halls of the Daya Bay Reactor Neutrino Experiment at varying depth is reported. The neutron yield in Daya Bay's liquid scintillator is measured to be $Y_n=(10.26\pm 0.86)\times 10^{-5}$, $(10.22\pm 0.87)\times 10^{-5}$, and $(17.03\pm 1.22)\times 10^{-5}~\mu^{-1}~$g$^{-1}~$cm$^2$ at depths of 250, 265, and 860 meters-water-equivalent. These results are compared to other measurements and the simulated neutron yield in Fluka and Geant4. A global fit including the Daya Bay measurements yields a power law coefficient of $0.77 \pm 0.03$ for the dependence of the neutron yield on muon energy.
• ### Seasonal Variation of the Underground Cosmic Muon Flux Observed at Daya Bay(1708.01265)

Jan. 8, 2018 hep-ex, physics.ins-det
The Daya Bay Experiment consists of eight identically designed detectors located in three underground experimental halls named as EH1, EH2, EH3, with 250, 265 and 860 meters of water equivalent vertical overburden, respectively. Cosmic muon events have been recorded over a two-year period. The underground muon rate is observed to be positively correlated with the effective atmospheric temperature and to follow a seasonal modulation pattern. The correlation coefficient $\alpha$, describing how a variation in the muon rate relates to a variation in the effective atmospheric temperature, is found to be $\alpha_{\text{EH1}} = 0.362\pm0.031$, $\alpha_{\text{EH2}} = 0.433\pm0.038$ and $\alpha_{\text{EH3}} = 0.641\pm0.057$ for each experimental hall.
• ### Phosphorus-bearing molecules in the Galactic Center(1712.07006)

Dec. 19, 2017 astro-ph.GA
Phosphorus (P) is one of the essential elements for life due to its central role in biochemical processes. Recent searches have shown that P-bearing molecules (in particular PN and PO) are present in star-forming regions, although their formation routes remain poorly understood. In this Letter, we report observations of PN and PO towards seven molecular clouds located in the Galactic Center, which are characterized by different types of chemistry. PN is detected in five out of seven sources, whose chemistry is thought to be shock-dominated. The two sources with PN non-detections correspond to clouds exposed to intense UV/X-rays/cosmic-ray radiation. PO is detected only towards the cloud G+0.693$-$0.03, with a PO/PN abundance ratio of $\sim$1.5. We conclude that P-bearing molecules likely form in shocked gas as a result of dust grain sputtering, while are destroyed by intense UV/X-ray/cosmic ray radiation.
• ### The 95zr(n, gamma)96zr cross section from the surrogate ratio method and its effect on the s-process nucleosynthesis(1709.04635)

Sept. 14, 2017 nucl-ex
The 95Zr(n,gamma)96Zr reaction cross section is crucial in the modelling of s-process nucleosynthesis in asymptotic giant branch stars because it controls the operation of the branching point at the unstable 95Zr and the subsequent production of 96Zr. We have carried out the measurement of the 94Zr(18O,16O) and 90Zr(18O,16O) reactions and obtained the gamma-decay probability ratio of 96Zr* and 92Zr* to determine the 95Zr(n,gamma)96Zr reaction cross sections with the surrogate ratio method. Our deduced maxwellian-averaged cross section of 66+-16 mb at 30 keV is close to the value recommended by Bao et al. (2000), but 30% and more than a factor of two larger than the values proposed by Toukan & Kappeler (1990) and Lugaro et al. (2014), respectively, and routinely used in s-process models. We tested the new rate in stellar models with masses between 2 and 6 Msun and metallicities 0.014 and 0.03. The largest changes - up 80% variations in 96Zr - are seen in models of mass 3-4 Msun, where the 22Ne neutron source is mildly activated. The new rate can still provide a match to data from meteoritic stardust silicon carbide grains, provided the maximum mass of the parent stars is below 4 Msun, for a metallicity of 0.03.
• ### Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay(1704.01082)

June 20, 2017 hep-ex, nucl-ex, physics.ins-det
The Daya Bay experiment has observed correlations between reactor core fuel evolution and changes in the reactor antineutrino flux and energy spectrum. Four antineutrino detectors in two experimental halls were used to identify 2.2 million inverse beta decays (IBDs) over 1230 days spanning multiple fuel cycles for each of six 2.9 GW$_{\textrm{th}}$ reactor cores at the Daya Bay and Ling Ao nuclear power plants. Using detector data spanning effective $^{239}$Pu fission fractions, $F_{239}$, from 0.25 to 0.35, Daya Bay measures an average IBD yield, $\bar{\sigma}_f$, of $(5.90 \pm 0.13) \times 10^{-43}$ cm$^2$/fission and a fuel-dependent variation in the IBD yield, $d\sigma_f/dF_{239}$, of $(-1.86 \pm 0.18) \times 10^{-43}$ cm$^2$/fission. This observation rejects the hypothesis of a constant antineutrino flux as a function of the $^{239}$Pu fission fraction at 10 standard deviations. The variation in IBD yield was found to be energy-dependent, rejecting the hypothesis of a constant antineutrino energy spectrum at 5.1 standard deviations. While measurements of the evolution in the IBD spectrum show general agreement with predictions from recent reactor models, the measured evolution in total IBD yield disagrees with recent predictions at 3.1$\sigma$. This discrepancy indicates that an overall deficit in measured flux with respect to predictions does not result from equal fractional deficits from the primary fission isotopes $^{235}$U, $^{239}$Pu, $^{238}$U, and $^{241}$Pu. Based on measured IBD yield variations, yields of $(6.17 \pm 0.17)$ and $(4.27 \pm 0.26) \times 10^{-43}$ cm$^2$/fission have been determined for the two dominant fission parent isotopes $^{235}$U and $^{239}$Pu. A 7.8% discrepancy between the observed and predicted $^{235}$U yield suggests that this isotope may be the primary contributor to the reactor antineutrino anomaly.
• ### 15N Fractionation in Infrared-Dark Cloud Cores(1705.04082)

May 11, 2017 astro-ph.GA, astro-ph.SR
Nitrogen is one of the most abundant elements in the Universe and its 14N/15N isotopic ratio has the potential to provide information about the initial environment in which our Sun formed. Recent findings suggest that the Solar System may have formed in a massive cluster since the presence of short-lived radioisotopes in meteorites can only be explained by the influence of a supernova. The aim of this project is to determine the 14N/15N ratio towards a sample of cold, massive dense cores at the initial stages in their evolution. We have observed the J=1-0 transitions of HCN, H13CN, HC15N, HN13C and H15NC toward a sample of 22 cores in 4 Infrared-Dark Clouds (IRDCs). IRDCs are believed to be the precursors of high-mass stars and star clusters. Assuming LTE and a temperature of 15K, the column densities of HCN, H13CN, HC15N, HN13C and H15NC are calculated and their 14N/15N ratio is determined for each core. The 14N/15N ratio measured in our sample of IRDC cores range between ~70 and >763 in HCN and between ~161 and ~541 in HNC. They are consistent with the terrestrial atmosphere (TA) and protosolar nebula (PSN) values, and with the ratios measured in low-mass pre-stellar cores. However, the 14N/15N ratios measured in cores C1, C3, F1, F2 and G2 do not agree with the results from similar studies toward the same massive cores using nitrogen bearing molecules with nitrile functional group (-CN) and nitrogen hydrides (-NH) although the ratio spread covers a similar range. Amongst the 4 IRDCs we measured relatively low 14N/15N ratios towards IRDC G which are comparable to those measured in small cosmomaterials and protoplanetary disks. The low average gas density of this cloud suggests that the gas density, rather than the gas temperature, may be the dominant parameter influencing the initial nitrogen isotopic composition in young PSN.
• ### Improved Measurement of the Reactor Antineutrino Flux and Spectrum at Daya Bay(1607.05378)

Jan. 9, 2017 hep-ex, nucl-ex, physics.ins-det
A new measurement of the reactor antineutrino flux and energy spectrum by the Daya Bay reactor neutrino experiment is reported. The antineutrinos were generated by six 2.9~GW$_{\mathrm{th}}$ nuclear reactors and detected by eight antineutrino detectors deployed in two near (560~m and 600~m flux-weighted baselines) and one far (1640~m flux-weighted baseline) underground experimental halls. With 621 days of data, more than 1.2 million inverse beta decay (IBD) candidates were detected. The IBD yield in the eight detectors was measured, and the ratio of measured to predicted flux was found to be $0.946\pm0.020$ ($0.992\pm0.021$) for the Huber+Mueller (ILL+Vogel) model. A 2.9~$\sigma$ deviation was found in the measured IBD positron energy spectrum compared to the predictions. In particular, an excess of events in the region of 4-6~MeV was found in the measured spectrum, with a local significance of 4.4~$\sigma$. A reactor antineutrino spectrum weighted by the IBD cross section is extracted for model-independent predictions.
• ### Limits on Active to Sterile Neutrino Oscillations from Disappearance Searches in the MINOS, Daya Bay, and Bugey-3 Experiments(1607.01177)

Oct. 17, 2016 hep-ex
Searches for a light sterile neutrino have been performed independently by the MINOS and the Daya Bay experiments using the muon (anti)neutrino and electron antineutrino disappearance channels, respectively. In this Letter, results from both experiments are combined with those from the Bugey-3 reactor neutrino experiment to constrain oscillations into light sterile neutrinos. The three experiments are sensitive to complementary regions of parameter space, enabling the combined analysis to probe regions allowed by the LSND and MiniBooNE experiments in a minimally extended four-neutrino flavor framework. Stringent limits on $\sin^2 2\theta_{\mu e}$ are set over 6 orders of magnitude in the sterile mass-squared splitting $\Delta m^2_{41}$. The sterile-neutrino mixing phase space allowed by the LSND and MiniBooNE experiments is excluded for $\Delta m^2_{41} < 0.8$ eV$^2$ at 95% CL$_s$.
• ### Measurement of electron antineutrino oscillation based on 1230 days of operation of the Daya Bay experiment(1610.04802)

Oct. 16, 2016 hep-ex, nucl-ex, physics.ins-det
A measurement of electron antineutrino oscillation by the Daya Bay Reactor Neutrino Experiment is described in detail. Six 2.9-GW$_{\rm th}$ nuclear power reactors of the Daya Bay and Ling Ao nuclear power facilities served as intense sources of $\overline{\nu}_{e}$'s. Comparison of the $\overline{\nu}_{e}$ rate and energy spectrum measured by antineutrino detectors far from the nuclear reactors ($\sim$1500-1950 m) relative to detectors near the reactors ($\sim$350-600 m) allowed a precise measurement of $\overline{\nu}_{e}$ disappearance. More than 2.5 million $\overline{\nu}_{e}$ inverse beta decay interactions were observed, based on the combination of 217 days of operation of six antineutrino detectors (Dec. 2011--Jul. 2012) with a subsequent 1013 days using the complete configuration of eight detectors (Oct. 2012--Jul. 2015). The $\overline{\nu}_{e}$ rate observed at the far detectors relative to the near detectors showed a significant deficit, $R=0.949 \pm 0.002(\mathrm{stat.}) \pm 0.002(\mathrm{syst.})$. The energy dependence of $\overline{\nu}_{e}$ disappearance showed the distinct variation predicted by neutrino oscillation. Analysis using an approximation for the three-flavor oscillation probability yielded the flavor-mixing angle $\sin^22\theta_{13}=0.0841 \pm 0.0027(\mathrm{stat.}) \pm 0.0019(\mathrm{syst.})$ and the effective neutrino mass-squared difference of $\left|{\Delta}m^2_{\mathrm{ee}}\right|=(2.50 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$. Analysis using the exact three-flavor probability found ${\Delta}m^2_{32}=(2.45 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ assuming the normal neutrino mass hierarchy and ${\Delta}m^2_{32}=(-2.56 \pm 0.06(\mathrm{stat.}) \pm 0.06(\mathrm{syst.})) \times 10^{-3}\ {\rm eV}^2$ for the inverted hierarchy.
• ### Improved Search for a Light Sterile Neutrino with the Full Configuration of the Daya Bay Experiment(1607.01174)

Oct. 11, 2016 hep-ex
This Letter reports an improved search for light sterile neutrino mixing in the electron antineutrino disappearance channel with the full configuration of the Daya Bay Reactor Neutrino Experiment. With an additional 404 days of data collected in eight antineutrino detectors, this search benefits from 3.6 times the statistics available to the previous publication, as well as from improvements in energy calibration and background reduction. A relative comparison of the rate and energy spectrum of reactor antineutrinos in the three experimental halls yields no evidence of sterile neutrino mixing in the $2\times10^{-4} \lesssim |\Delta m^{2}_{41}| \lesssim 0.3$ eV$^{2}$ mass range. The resulting limits on $\sin^{2}2\theta_{14}$ are improved by approximately a factor of 2 over previous results and constitute the most stringent constraints to date in the $|\Delta m^{2}_{41}| \lesssim 0.2$ eV$^{2}$ region.
• ### Study of the wave packet treatment of neutrino oscillation at Daya Bay(1608.01661)

Aug. 5, 2016 hep-ph, hep-ex
The disappearance of reactor $\bar{\nu}_e$ observed by the Daya Bay experiment is examined in the framework of a model in which the neutrino is described by a wave packet with a relative intrinsic momentum dispersion $\sigma_\text{rel}$. Three pairs of nuclear reactors and eight antineutrino detectors, each with good energy resolution, distributed among three experimental halls, supply a high-statistics sample of $\bar{\nu}_e$ acquired at nine different baselines. This provides a unique platform to test the effects which arise from the wave packet treatment of neutrino oscillation. The modified survival probability formula was used to fit Daya Bay data, providing the first experimental limits: $2.38 \cdot 10^{-17} < \sigma_{\rm rel} < 0.23$. Treating the dimensions of the reactor cores and detectors as constraints, the limits are improved: $10^{-14} \lesssim \sigma_{\rm rel} < 0.23$, and an upper limit of $\sigma_{\rm rel} <0.20$ is obtained. All limits correspond to a 95\% C.L. Furthermore, the effect due to the wave packet nature of neutrino oscillation is found to be insignificant for reactor antineutrinos detected by the Daya Bay experiment thus ensuring an unbiased measurement of the oscillation parameters $\sin^22\theta_{13}$ and $\Delta m^2_{32}$ within the plane wave model.
• ### The Detector System of The Daya Bay Reactor Neutrino Experiment(1508.03943)

Jan. 8, 2016 hep-ex, physics.ins-det
The Daya Bay experiment was the first to report simultaneous measurements of reactor antineutrinos at multiple baselines leading to the discovery of $\bar{\nu}_e$ oscillations over km-baselines. Subsequent data has provided the world's most precise measurement of $\rm{sin}^22\theta_{13}$ and the effective mass splitting $\Delta m_{ee}^2$. The experiment is located in Daya Bay, China where the cluster of six nuclear reactors is among the world's most prolific sources of electron antineutrinos. Multiple antineutrino detectors are deployed in three underground water pools at different distances from the reactor cores to search for deviations in the antineutrino rate and energy spectrum due to neutrino mixing. Instrumented with photomultiplier tubes (PMTs), the water pools serve as shielding against natural radioactivity from the surrounding rock and provide efficient muon tagging. Arrays of resistive plate chambers over the top of each pool provide additional muon detection. The antineutrino detectors were specifically designed for measurements of the antineutrino flux with minimal systematic uncertainty. Relative detector efficiencies between the near and far detectors are known to better than 0.2%. With the unblinding of the final two detectors' baselines and target masses, a complete description and comparison of the eight antineutrino detectors can now be presented. This paper describes the Daya Bay detector systems, consisting of eight antineutrino detectors in three instrumented water pools in three underground halls, and their operation through the first year of eight detector data-taking.
• ### JUNO Conceptual Design Report(1508.07166)

Sept. 28, 2015 hep-ex, physics.ins-det
The Jiangmen Underground Neutrino Observatory (JUNO) is proposed to determine the neutrino mass hierarchy using an underground liquid scintillator detector. It is located 53 km away from both Yangjiang and Taishan Nuclear Power Plants in Guangdong, China. The experimental hall, spanning more than 50 meters, is under a granite mountain of over 700 m overburden. Within six years of running, the detection of reactor antineutrinos can resolve the neutrino mass hierarchy at a confidence level of 3-4$\sigma$, and determine neutrino oscillation parameters $\sin^2\theta_{12}$, $\Delta m^2_{21}$, and $|\Delta m^2_{ee}|$ to an accuracy of better than 1%. The JUNO detector can be also used to study terrestrial and extra-terrestrial neutrinos and new physics beyond the Standard Model. The central detector contains 20,000 tons liquid scintillator with an acrylic sphere of 35 m in diameter. $\sim$17,000 508-mm diameter PMTs with high quantum efficiency provide $\sim$75% optical coverage. The current choice of the liquid scintillator is: linear alkyl benzene (LAB) as the solvent, plus PPO as the scintillation fluor and a wavelength-shifter (Bis-MSB). The number of detected photoelectrons per MeV is larger than 1,100 and the energy resolution is expected to be 3% at 1 MeV. The calibration system is designed to deploy multiple sources to cover the entire energy range of reactor antineutrinos, and to achieve a full-volume position coverage inside the detector. The veto system is used for muon detection, muon induced background study and reduction. It consists of a Water Cherenkov detector and a Top Tracker system. The readout system, the detector control system and the offline system insure efficient and stable data acquisition and processing.
• ### Proton widths of the low-lying 16F states from the 15N(7Li, 6Li)16N reaction(1403.6912)

March 27, 2014 nucl-ex
All the 16F levels are unbound by proton emission. To date the four low-lying 16F levels below 1 MeV have been experimentally identified with well established spin-parity values and excitation energies with an accuracy of 4 - 6 keV. However, there are still considerable discrepancies for their level widths. The present work aims to explore these level widths through an independent method. The angular distributions of the 15N(7Li, 6Li)16N reaction leading to the first four states in 16N were measured using a high-precision Q3D magnetic spectrograph. The neutron spectroscopic factors and the asymptotic normalization coefficients for these states in 16N were then derived based on distorted wave Born approximation analysis. The proton widths of the four low-lying resonant states in 16F were obtained according to charge symmetry of strong interaction.
• ### Spectroscopic factors for low-lying 16N levels and the astrophysical 15N(n,gamma)16N reaction rate(1401.0067)

Dec. 31, 2013 nucl-ex
Fluorine is a key element for nucleosynthetic studies since it is extremely sensitive to the physical conditions within stars. The astrophysical site to produce fluorine is suggested to be asymptotic giant branch (AGB) stars. In these stars the 15N(n, g)16N reaction could affect the abundance of fluorine by competing with 15N(a, g)19F. The 15N(n, g)16N reaction rate depends directly on the neutron spectroscopic factors of the low-lying states in 16N. The angular distributions of the 15N(7Li, 6Li)16N reaction populating the ground state and the first three excited states in 16N are measured using a Q3D magnetic spectrograph and are used to derive the spectroscopic factors of these states based on distorted wave Born approximation (DWBA) analysis. The spectroscopic factors of these four states are extracted to be 0.96+-0.09, 0.69+-0.09, 0.84+-0.08 and 0.65+-0.08, respectively. Based on the new spectroscopic factors we derive the 15N(n,g)16N reaction rate. The accuracy and precision of the spectroscopic factors are enhanced due to the first application of high-precision magnetic spectrograph for resolving the closely-spaced 16N levels which can not be achieved in most recent measurement. The present result demonstrates that two levels corresponding to neutron transfers to the 2s1/2 orbit in 16N are not so good single-particle levels although 15N is a closed neutron-shell nucleus. This finding is contrary to the shell model expectation. The present work also provides an independent examination to shed some light on the existing discrepancies in the spectroscopic factors and the 15N(n, g)16N rate.
• ### Coherent control at its most fundamental: CEP-dependent electron localization in photodissoziation of a H2+ molecular ion beam target(1306.1755)

Measurements and calculations of the absolute carrier-envelope phase (CEP) effects in the photodissociation of the simplest molecule, H2+, with a 4.5-fs Ti:Sapphire laser pulse at intensities up to (4 +- 2)x10^14 Watt/cm^2 are presented. Localization of the electron with respect to the two nuclei (during the dissociation process) is controlled via the CEP of the ultra-short laser pulses. In contrast to previous CEP-dependent experiments with neutral molecules, the dissociation of the molecular ions is not preceded by a photoionization process, which strongly influences the CEP dependence. Kinematically complete data is obtained by time- and position-resolved coincidence detection. The phase dependence is determined by a single-shot phase measurement correlated to the detection of the dissoziation fragments. The experimental results show quantitative agreement with ab inito 3D-TDSE calculations that include nuclear vibration and rotation.
• ### Carrier-envelope phase control over pathway interference in strong-field dissociation of H$_2^+$(1306.4010)

The dissociation of an H$_2^+$ molecular-ion beam by linearly polarized, carrier-envelope-phase-tagged 5 fs pulses at 4$\times10^{14}$W/cm$^2$ with a central wavelength of 730 nm was studied using a coincidence 3D momentum imaging technique. Carrier-envelope-phase-dependent asymmetries in the emission direction of H$^+$ fragments relative to the laser polarization were observed. These asymmetries are caused by interference of odd and even photon number pathways, where net-zero photon and 1-photon interference predominantly contributes at H$^+$+H kinetic energy releases of 0.2 -- 0.45 eV, and net-2-photon and 1-photon interference contributes at 1.65 -- 1.9 eV. These measurements of the benchmark H$_2^+$ molecule offer the distinct advantage that they can be quantitatively compared with \textit{ab initio} theory to confirm our understanding of strong-field coherent control via the carrier-envelope phase.
• ### New determination of the 13C(a, n)16O reaction rate and its influence on the s-process nucleosynthesis in AGB stars(1208.0714)

Aug. 3, 2012 nucl-ex, astro-ph.SR
We present a new measurement of the $\alpha$-spectroscopic factor ($S_\alpha$) and the asymptotic normalization coefficient (ANC) for the 6.356 MeV 1/2$^+$ subthreshold state of $^{17}$O through the $^{13}$C($^{11}$B, $^{7}$Li)$^{17}$O transfer reaction and we determine the $\alpha$-width of this state. This is believed to have a strong effect on the rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction, the main neutron source for {\it slow} neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. Based on the new width we derive the astrophysical S-factor and the stellar rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction. At a temperature of 100 MK our rate is roughly two times larger than that by \citet{cau88} and two times smaller than that recommended by the NACRE compilation. We use the new rate and different rates available in the literature as input in simulations of AGB stars to study their influence on the abundances of selected $s$-process elements and isotopic ratios. There are no changes in the final results using the different rates for the $^{13}$C($\alpha$, $n$)$^{16}$O reaction when the $^{13}$C burns completely in radiative conditions. When the $^{13}$C burns in convective conditions, as in stars of initial mass lower than $\sim$2 $M_\sun$ and in post-AGB stars, some changes are to be expected, e.g., of up to 25% for Pb in our models. These variations will have to be carefully analyzed when more accurate stellar mixing models and more precise observational constraints are available.
• ### Asymptotic normalization coefficient from the 12C(7Li,6He)13N reaction and the astrophysical 12C(p,g)13N reaction rate(0905.1530)

May 11, 2009 nucl-ex
Angular distribution of the 12C(7Li,6He)13N reaction at E(7Li) = 44.0 MeV was measured at the HI-13 tandem accelerator of Beijing, China. Asymptotic normalization coefficient (ANC) of 13N = 12C + p was derived to be 1.64 $\pm$ 0.11 fm$^{-1/2}$ through distorted wave Born approximation (DWBA) analysis. The ANC was then used to deduce the astrophysical $S(E)$ factors and reaction rates for direct capture in 12C(p,g)13N at energies of astrophysical relevance.
• ### The origin of the red luminescence in Mg-doped GaN(cond-mat/0602104)

Feb. 4, 2006 cond-mat.mtrl-sci
Optically-detected magnetic resonance (ODMR) and positron annihilation spectroscopy (PAS) experiments have been employed to study magnesium-doped GaN layers grown by metal-organic vapor phase epitaxy. As the Mg doping level is changed, the combined experiments reveal a strong correlation between the vacancy concentrations and the intensity of the red photoluminescence band at 1.8 eV. The analysis provides strong evidence that the emission is due to recombination in which electrons both from effective mass donors and from deeper donors recombine with deep centers, the deep centers being vacancy-related defects.