Precision Measurement of the Speciﬁc Activity of 39 Ar in Atmospheric Argon with the DEAP-3600 Detector

The speciﬁc activity of the beta decay of 39 Ar in atmospheric argon is measured using the DEAP-3600 detector. DEAP-3600, located 2 km underground at SNOLAB, uses a total of (3269 ± 24) kg of liquid argon distilled from the atmosphere to search for dark matter. This detector with very low background uses pulseshape discrimination to differentiate between nuclear recoils and electron recoils and is well-suited to measure the decay of 39 Ar. With 167 live-days of data, the measured speciﬁc activity at the time of atmospheric extraction is [0.964 ± 0.001 (stat) ± 0.024 (sys)] Bq/kg atmAr which is consistent with results from other experiments. A cross-check analysis using different event selection criteria provides a consistent result.

While the production of 39 Ar in the atmosphere is in equilibrium, measurements of ice cores and tree rings by Gu et al. [14] show the 39 Ar/Ar ratio has varied by as much as 17% in the past 2500 years.Recent measurements of the specific activity of 39 Ar in atmospheric argon, S Ar39 , were realized by the WARP collaboration with a result of S Ar39 = (1.01 ± 0.02 stat ± 0.08 sys ) Bq/kg atmAr [15] and by the ArDM collaboration with S Ar39 = (0.95 ± 0.05) Bq/kg atmAr [16].
This paper describes the measurement of the 39 Ar specific activity using the DEAP-3600 detector [9], located 2 km underground in Creighton Mine at SNOLAB in Sudbury, Ontario, Canada.DEAP-3600 is a dark matter experiment with a liquid argon (LAr) target which achieves low-background levels due to both its use of low-background construction materials and implementation of pulse-shape discrimination (PSD).The PSD technique is able to differentiate between nuclear recoils and electron recoils; it achieves an expected a Deceased.b deap-papers@snolab.cac Currently at SNOLAB, Lively, Ontario, P3Y 1M3, Canada d Currently at Nikhef and the University of Amsterdam, Science Park, 1098XG Amsterdam, Netherlands leakage of electron recoil events into the dark matter search region of interest of fewer than 1 event per year of data.The large mass of atmospheric argon and the very low background levels achieved with this experiment [17,18] enable the precision specific activity measurement presented here.The specific activity is calculated by estimating the total number of 39 Ar decays N within a certain live-time T live as shown in Eq. 1.
where m LAr is the mass of LAr in the detector.
A brief description of DEAP-3600 is provided in Section 2. A dedicated estimate of the LAr mass in DEAP-3600 is presented in Section 3. The dataset and livetime calculation, as well as the event selection are described in Section 4. The measurement of N is presented in Section 5, alongside details on the S Ar39 calculation, the systematic uncertainties, and the results.Section 6 briefly describes an updated version of the measurement from Ref. [19] which appears here as a cross-check.Concluding remarks are given in Section 7.

The Detector and Data Acquisition System
The DEAP-3600 experiment operated an ultra-pure LAr target of over 3 tonnes held in a spherical acrylic vessel (AV) with 85 cm radius from November 2016 to April 2020.The atmospheric argon for the LAr target was procured from Air Liquide.
Connected to the top of the AV is an acrylic neck surrounding a liquid nitrogen (LN 2 ) filled stainless steel cooling coil which condenses the gaseous argon (GAr) contained in the top of the AV.The AV was partially filled with the GAr/LAr interface approximately 55 cm above the equator.After filling with LAr the detector was sealed and this volume of argon remained within the AV for the duration of data taking.Cooling of the LAr was achieved by continuous circulation of LN 2 within the cooling coil.
Coated on the inner surface of the AV is a layer of tetraphenyl butadiene (TPB).The TPB wavelength-shifts the 128 nm ultraviolet (UV) scintillation light from the LAr target into the visible spectrum with a peak at 420 nm [20].This light is detected by 255 photomultiplier tubes (PMTs) which point inward and are optically coupled to the AV by acrylic light guides.The PMTs are distributed in rings around the AV, with the PMTs in each ring having the same vertical position.The AV and PMTs are enclosed in a stainless steel shell which is continuously flushed with radon-scrubbed nitrogen gas.Installed on the outer surface of the shell are 48 PMTs which point outward and, combined with the water held within a cylindrical tank surrounding the shell, act as a muon veto system.This muon veto system detects Cherenkov light produced by muon interactions within the water.A schematic of the detector is shown in Figure 1.
Temperature sensors are placed around the AV at 85 locations along filler blocks which are mounted in the spaces between the PMTs.The sensors are spread around the AV and placed at distances of 0.9 m, 1.1 m, or 1.3 m from the centre of the AV.These sensors, along with the temperature and pressure within the LN 2 cooling system, are monitored and logged using a slow control system.
Within the data acquisition system (DAQ) each PMT is connected to a channel on a custom-built signal conditioning board (SCB) which achieves the high voltage decoupling and shapes the signals.The SCB outputs are transmitted to high-gain (CAEN V1720) and low-gain (CAEN V1740) waveform digitizers.These digitizers convert a continuous analogue signal to a discrete digital signal using analogueto-digital converters (ADCs).The summed input from each SCB is also passed to a digitizer and trigger module (DTM) which resolves the trigger criteria based on two rolling charge integrals: a narrow integral Q n over a 177 ns window and a wide integral Q w over a 3.1 µs window.The promptness of the signal is computed by the Q n /Q w fraction.Five trigger regions are defined based on these three variables.A prescaling factor of 100 is applied to events in the energy range of Q n ≈ [50, 565] keV ee in the low Q n /Q w region.This prescaling predominantly affects 39 Ar decays and reduces the available statistics by storing only the observed PMT waveforms for precisely 1 out of every 100 events.In a 24 hour period, roughly 2.7 × 10 6 39 Ar events remain after the prescaling.The timestamp of every event, included those which are prescaled, is recorded in the data.
The DTM makes the decision to trigger based on the summed value of Q n from all 255 PMTs and sends a trigger signal to the digitizers if this value passes a threshold of 19 PE.Each digitizer channel records PMT waveforms for 16 µs upon receiving a trigger signal, including a pre-trigger window of 2.4 µs.The data acquisition system is operated by MIDAS [21] and the data are analyzed with RAT [22], a software framework built on GEANT4 [23] and ROOT [24].The observed charge in each PMT is integrated over a window of [-28, 10000] ns relative to the event time.This charge is divided by the single photoelectron (PE) charge for each PMT measured through independent calibration [25].The resulting PE number provides the energy estimator for the data.The PSD variable F prompt distinguishes nuclear recoil events at high F prompt from 39 Ar decays and electron recoil backgrounds (ERB) at low F prompt .For this measurement it is defined as the fraction of PE detected in a time window of [-28, 150] ns around the event time and is calculated as The ERB is composed of events generated by both γ-rays emitted by trace radioactivity in detector components and β -decays which scatter on electrons in the LAr.
A more detailed description of the DEAP-3600 detector can be found in Ref. [9].

Liquid Argon Mass Estimate
The LAr mass is determined by evaluating both its density and its volume within the AV.This method previously resulted in a LAr mass of (3279 ± 96) kg [17].That result has been refined for this work.
Two inputs are required to evaluate the volume of LAr in the detector: the AV radius and the LAr height within the AV.The internal radius of the AV was measured during its construction.After correcting for the thermal contraction that occurred during cool-down using a temperature-dependent coefficient measured in Ref. [26] the AV radius is determined to be (845.6± 0.9) mm.The LAr height is measured by taking advantage of the total internal reflection of the UV light at the GAr/LAr interface.The TPB re-emits light isotropically and so the photon detection rates for each PMT depend on the area of visible TPB immersed in the LAr.The rates for every PMT ring are averaged and the distribution is fit with an analytic model of the corresponding immersed area.This method is validated by comparing the data to Monte Carlo simulations of 39 Ar decays within the LAr while varying the simulated LAr height.The best fit is found at a LAr height of (550 ± 10) mm above the equator and is stable across the dataset.The systematic uncertainty on the LAr height is the dominant source of uncertainty for the LAr mass estimate.A cross-check using the position reconstruction of detected events to evaluate the LAr height provides a consistent result.In this cross-check, a template fit in the reconstructed vertical position of 39 Ar decay events is performed by comparing simulations with different LAr height values to the data histogram.
The LAr density is a function of its temperature.This temperature is constrained by the liquid-vapor transition of the argon in the AV and by the liquid-vapor transition of the nitrogen in the cooling coil.As the pressure in both systems is constantly recorded, the average LAr temperature is known within a few degrees K, and thus the effective density can be established to 0.5% precision.
The possibility of argon bubbles is also investigated, the presence of which would reduce the total mass of LAr.Using the behavior of nitrogen as a reference [27] a limiting case is considered where all of the exterior heat entering the LAr creates bubbles.This worst-case scenario indicates that at most 6.3 kg of LAr is displaced by bubbles.
A toy Monte Carlo sampling the probability distribution functions (PDFs) of the AV radius, the LAr height, the LAr density, and bubble displacement is used to determine the central value of the LAr mass and its uncertainty.Flat PDFs are used for the constraints on the LAr density and the bubble displacement, while the AV radius and LAr height PDFs are considered Gaussian.According to this method, during the data-taking period of this measurement the DEAP-3600 AV contained m LAr = (3269 ± 24) kg of LAr.

Run Selection
The dataset is divided into discrete runs during which signals from the LAr are recorded.A single run is typically about 22 hours long, though this can vary between just a few minutes and up to about 2 days.The runs examined here are from the same 2016-2017 dataset used for the dark matter search published by the DEAP collaboration [17] with the additional restriction that runs are at least 18 hours long.This requirement is imposed to ensure sufficient statistics to fit the γ-dominated region of the ERB spectrum in each run.
The selection of runs is based on stability criteria concerning the cooling system of the AV, the charge distributions in the PMTs, and the efficiency of the trigger.A data cleaning cut is applied to each run to reject events occurring within δt cut = 32 µs of the previous event, which removes δt i ≤ 32µs of livetime per event i; the total number of events removed by this cut is N DCcut .Low-level cuts are then applied to reject events recorded from pulse injections by periodic monitoring triggers and events with inconsistent data acquisition readouts such as busy signals, for a total of N LLcut events.The events removed by these cuts, along with all remaining physics triggers N phys , are taken into account in the run-dependent livetime calculation shown in Eq. 3.
Here, T live is the livetime for a run, T run is the total time of that run, and δt int = 10 µs corresponds to the charge integration window during which the detector can record a pile-up event, while the time between δt int and δt cut is dead time.The value of N phys includes the prescaled triggers as the timestamp of each of these events is stored.Testing of the algorithm was performed using values of δt cut ranging from 20 µs up to 250 µs for a selection of data runs.For each δt cut value the livetime and specific activity of each run were calculated.We observed negligibly small variations in the measured specific activity as a function of δt cut , as expected.
An offline reduction is applied where precisely 1 out of every 100 events from outside the prescaled trigger region is kept in order to remove boundary effects and obtain a smooth spectrum.

Event Selection
In addition to the data cleaning and low-level cuts described in the previous section, event selection cuts are applied.Pileup needs to be taken into account given the high rate of 39 Ar decays and the length of the event window: approximately 5% of recorded events are expected to contain 2 or more decays.Additionally, a triggered event can follow an energy deposit which occurred during time in which DAQ was busy and unable to record (deadtime).The late scintillation light from this previous, unrecorded energy deposit can reach into the beginning of the triggered event.Since the full energy of the previous energy deposit is not visible in the digitized trace, this type of pile-up is hard to model.While this analysis endeavours to keep pile-up events and account for them in the specific activity calculation, events with this pre-trigger pileup are not suitable analysis candidates.
To select events without pre-trigger pile-up the time at which the event occurred within the trigger window must be in the range [2250, 2700] ns, and it is required that fewer than 4 pulses are recorded by the PMTs in the first 1600 ns of the event.These cuts do not remove a significant number of events, and the majority of the events removed are at very low energies.These removed events are mainly outside the range of the fits described in Section 5.1.
Electron recoil events, which are dominated by 39 Ar decays at lower energies and γ backgrounds at higher energies, are selected with the requirement 0.1 ≤ F prompt ≤ 0.5.These events, along with 39 Ar-39 Ar pileup events and 39 Ar signal events, are shown in Figure 2. A more in-depth discussion of the electron recoil events can be found in Ref. [18].This analysis also makes use of a peak-finding algorithm based which examines the PMT waveforms to count the number of "sub-events" within the trigger window.The algorithm counts pulses from each PMT in the event window to look for statistically significant increases in the pulse count and is able to identify sub-events separated by as little as 50 ns.When tested using MC the algorithm was able to correctly identify 96% of 39 Ar-39 Ar pile-up events, and only 0.1% of single 39 Ar were incorrectly identified as having multiple subevents.The number of sub-events is used to select pile-up candidates in order to perform a data-driven estimate of the double 39 Ar pile-up cut efficiency as described in Section 5.2.

Specific Activity Measurement
The specific activity of 39 Ar is measured individually for each run in the dataset.Each measurement is based on a fit to the low F prompt energy spectrum and consists of an 39 Ar β -decay spectrum (single 39 Ar), a spectrum with two 39 Ar decays occurring within the same trigger window (double 39 Ar pile-up), and a spectrum containing all non-39 Ar ERB events scaled to the activities measured in Ref. [18].The ERB and double 39 Ar pile-up input spectra are generated by simulating events within the DEAP-3600 detector using the RAT software.The single 39 Ar component is built directly from the theoretical model provided by Kostensalo et al. [28].Each of the three model components is normalized using a parameter in the fits.Energy scale PE and detector resolution effects σ (PE) in the form of a Gaussian term are applied to all three model components, parameterized as The constant energy scale parameter p 0 is fixed in the fits to a value obtained by measuring PMT baselines.The number of 39 Ar decays is split into two main components as with the number of single 39 Ar decays N single and the number of 39 Ar decays which are part of a pile-up event N pile−up .The latter number includes double 39 Ar decays, triple 39 Ar decays (three 39 Ar decays in one trigger window), pile-up of 39 Ar decays with ERB decays, and pile-up of 39 Ar decays with high F prompt events (F prompt > 0.5).

Fitting the Energy Spectrum
This analysis uses Minuit in ROOT [29] to fit the three input spectra to data.The fit performs a chi-square minimization as with the number of bins n b in the data histogram, the data content D i in bin i and the model contribution M i .The parameter P is a penalty term applied to a shape nuisance parameter a 0 , which corrects the theoretical 39 Ar input spectrum linearly in energy to fit the data; it is constructed to account for the differences ∆ a 0 = 0.01 observed between the Kostensalo et al. [28] and the Behrens & Janecke [30] 39 Ar β -shapes and is calculated as The closer a 0 is to zero, the more the spectrum fits the shape by Kostensalo et al.The fit model is adapted from the model used for the energy response fits in Ref. [17].
The fit range for this measurement is [200, 11000] PE and is chosen to avoid trigger efficiency effects at low PE and to provide a handle for the fit to scale n double and n ERB beyond the 39 Ar endpoint at high PE.This range includes the 40 K γ-emission peak at 1460 keV [18] (approximately 10,500 PE) and allows the ERB spectrum normalization to be determined in the fit.This analysis is performed with n b corresponding to a bin width of b = 20 PE which was chosen to provide sufficient statistics to define the 40 K peak.The 39 Ar β -shape nuisance parameter is an output of the fit.Details of the fit inputs and outputs are provided in Table 1, along with the other parameter values taken as input to the specific activity measurement.Figure 2 shows an example fit using this model for one data run.The parameters from each fit are examined to look for trends across the dataset and for issues such as getting stuck at the boundaries of their allowed ranges: no such issues are observed.

Calculating the Specific Activity
The number of 39 Ar single decays N single is obtained from the single 39 Ar spectrum fit result integral n as with the bin width b of the fitted data histogram, a trigger prescaling correction factor a presc = 100 and the cut efficiency ε.The main, data-driven method used to estimate this cut efficiency involves defining a loose event selection for the denominator spectrum of events present before the cuts with an F prompt < 0.7 requirement; this loose event selection  Fig. 2 The top panel shows an example fit on one run including the 39 Ar, ERB, and 39 Ar pile-up components which form the fit function.
The fit range from 200-11,000 PE is shown by the vertical dashed lines.The bottom panel shows the residual between the fit function and data normalized to the square root of the contents in the observed PE distribution.The fit is extrapolated below the lower bound to count events in the low-energy region where the trigger efficiency is not 100%.The slight mismodelling of the ERB spectrum as apparent in the residual plot does not significantly affect the final result.The reduced chi-square for this run is given by χ 2 /ndf = 687.42/531= 1.29.The run shown here is approximately 28.5 hours long.
removes the unwanted Cherenkov and nuclear recoil events.The numerator spectrum for the efficiency calculation contains those events which pass the selection cuts described in Section 4.2.First, the efficiency ε lowerPE is calculated over the range [300, 3000] PE, which is dominated by single 39 Ar events, by taking the ratio of the two spectra bin-by-bin.Then, to extract the efficiency ε for single 39 Ar events, a correction is applied to ε lowerPE to account for the presence of double 39 Ar pile-up in the sample.The cut efficiency ε double is calculated binby-bin in the PE histogram of each run over the range [300, 3000] PE with an additional event selection requirement for the numerator and denominator spectra to select events which contain exactly two sub-events.The efficiency for single 39 Ar events is calculated bin-by-bin by next solving for ε i using the bin-dependent ε lowerPE and ε double values in the following equation as where x i is the fraction of single 39 Ar events measured by the fits in bin i.The value for ε is then calculated as the average of the ε i values as any energy dependence here is negligible.The resulting data-driven value of ε is calculated run-by-run and each value is used in the calculations for its respective run.The average ε is 0.983.As a cross-check to this data-driven method, the cut efficiency values are evaluated using the Monte Carlo simulated samples of the 39 Ar decays and the double 39 Ar pile-up de-scribed earlier.While the simulations do not describe the data perfectly, this method yields clean spectra of these two event classes which can be individually analyzed.With the Monte Carlo method ε MC = 0.999.We evaluate the specific activity S Ar39 using both ε and ε MC as inputs to Eq. 8 and take the difference as a systematic uncertainty.This difference is the dominant source of systematic uncertainty for this measurement.
The number of 39 Ar decays that are part of pile-up events is split into the different components as Here, N double is the number of 39 Ar decays that are part of a double 39 Ar pile-up event, N triple is the number of 39 Ar decays that are part of a triple pile-up event, N ERB,Ar39 is the number of 39 Ar decays which pile-up with a ERB recoil, and N hFp,Ar39 is the number of 39 Ar decays which pile-up with a high F prompt process such as Cherenkov light or a nuclear recoil.N double is obtained from the double 39 Ar spectrum fit result integral n double as Here, ε double is the cut efficiency on double 39 Ar pile-up events described above, and the factor 2 corrects for 2 39 Ar decays in 1 double pile-up event.The energy-dependence of ε double over the wider range of the pile-up spectrum is taken into account by applying this correction bin-by-bin.N double is utilized to calculate the single 39 Ar rate R Ar39 from a first-order pile-up calculation as shown in Eq. 12.
R Ar39 is used to determine the remaining pile-up components which are estimated with first-order pile-up approximations as where the factor of 3 in N triple accounts for the 3 39 Ar decays in each of these pile-up events.The ERB rate R ERB = (10.5±0.6) Hz and the high F prompt rate R hFp = (270 ± 3) Hz are established from the fit output n ERB and from the rate of events observed in the high F prompt window in the dataset, respectively.The pile-up rates can be calculated by dividing the quantities in Eq. 13 by T live .
Beyond the ERB measured in Ref. [18], the dataset considered in this analysis may contain a small number of 85 Kr β -decay events.The 85 Kr beta spectrum has an endpoint energy of 687.0 keV; this is in the region dominated by the double 39 Ar pileup events.Uncertainty in the amplitude of a peak at 600.66 keV from the 226 Ra chain makes obtaining the 85 Kr from fitting the energy spectrum challenging.The 85 Kr activity is studied a posteriori by repeating the fit including the 85 Kr β -shape from Ref. [31] with a normalization parameter n Kr85 , while also varying the energy response parameters and the 39 Ar endpoint within their uncertainties.No cuts are made to remove the double 39 Ar pileup so that both the nominal fits and the fits including a 85 Kr spectrum are performed on the same data.These fit results suggest that at most 0.01 Bq/kg atmAr of 85 Kr is present in our dataset.This limit is considered as an additional source of systematic uncertainty.

Results
The specific activity of 39 Ar is evaluated for each run by combining Eqs. 8, 11 and 13 with Eq. 1.The run-by-run results are presented in Figure 3 which includes an exponential fit to the measured specific activity over time.This fit is used to determine the specific activity value at the start of the dataset.
Uncertainties due to the liquid argon mass estimate and related to the determination of cut efficiencies were discussed in Sections 3 and 5.2 respectively.Additional systematic uncertainties on the specific activity measurement are evaluated as follows.For each run, the fit is repeated with the linear energy scale parameter p 1 fixed to its central value for that run plus or minus 0.15 PE/keV.The uncertainties on the other energy scale and resolution parameters p 0 , p 3 and p 4 are likewise propagated to the measurement by repeating the fits with fixed parameters set according to their uncertainties determined in the energy response measurement described in Ref. [17].Uncertainties due to the choice of the histogram bin width (varied to 10 PE and to 40 PE) and the choice of the fit range (lower bound increased to 500 PE) are evaluated in a similar manner.Theoretical β -shape uncertainties are accounted for by repeating the fit with a 0 fixed to 0, and then fixed to the median value found over the entire dataset.The systematic uncertainty due to the ERB normalization is negligible, and so any systematics associated with the MC generation of the spectra used in the fits is similarly negligible.Optical model uncertainties within the MC do not affect the pile-up spectrum shape used in the fits.The systematic uncertainty due to the double 39 Ar pile-up spectrum shape and normalization are negligible.The impact of each source of systematic uncertainty on the result is detailed in Table 1.
The statistical uncertainty of 0.001 Bq/kg atmAr shown in Figure 3 is calculated by propagating the uncertainties on the ERB and the high F prompt background rates, and the fit uncertainties on the single 39 Ar and the double 39 Ar pile-up normalization parameters.The fit uncertainty of the 39 Ar spectrum dominates the statistical uncertainty.
A correction is applied to the measured specific activity determined from the exponential fit to account for the age of Fig. 3 The measured specific activity of 39 Ar versus run time for the entire dataset.The exponential trendline fit is shown, with the average statistical uncertainty depicted as an error band.The systematic uncertainty band is wider than the y-axis range shown here.
the argon.The correction factor is calculated as where t age = (1.0 ± 0.5) y is the average time between atmospheric extraction of the argon and the start of the data-taking period.Multiplying by η t corrects for the approximately 0.26 % drop in the activity before data were taken.Cosmogenic activation of 39 Ar during the time after the argon was extracted from the atmosphere is negligible.The specific activity of 39 Ar in atmospheric argon is measured to be S Ar39 = (0.964 ± 0.001 stat ± 0.024 sys ) Bq/kg atmAr .

Cross-Check Analysis
Here we present a cross-check to our result which is an update of an earlier analysis, the details of which are described in Ref. [19].This analysis used the Bayesian Analysis Toolkit (BAT) [32] software to fit the input spectra to the data and extract the model parameters.BAT uses Markov Chain Monte Carlo to generate posterior probability distributions of the fit parameters based on prior probability distributions and a likelihood function input by the user.This cross-check also differs from the analysis presented in previous sections by applying a different set of event selection cuts than those described in Section 4.2.A cut on the event time within the trigger window was not applied, and the peak-finding algorithm to count "sub-events" was used to remove the majority of pileup events.The data cleaning cut to remove events close in time to a previous event was not applied.Otherwise, the same criteria described in Section 4 were applied and a fit was performed on each run in the dataset.
The detector response model in the fit included a constant energy scale parameter p 0 , a linear energy scale parameter p 1 ,

Measurement
Specific activity [Bq/kg atmAr ] WARP [15] 1.01 ± 0.02 stat ± 0.08 sys ArDM [16] 0.95 ± 0.05 DEAP-3600 (this work) 0.964 ± 0.001 stat ± 0.024 sys a quadratic energy scale parameter p 2 , and a linear resolution parameter p 3 .The quadratic resolution parameter p 4 was not considered in these fits.The p 1 , p 2 , and p 3 model parameters were given flat priors in the fits and allowed to float.The p 0 parameter was fixed in the fits.
The nominal input 39 Ar spectrum used was from Behrens & Janecke [30].Each fit returned the normalization of this spectrum and was given a flat prior.In addition to the 39 Ar spectrum, the inputs to the fit were an ERB spectrum and an MC-generated spectrum of double 39 Ar pileup events which survive the cuts.Each fit returned the normalization parameters for these spectra; at the input stage these were given Gaussian priors with a mean value of 1, which corresponded to a normalization based on an assumed event rate and the known runtime.All three normalization parameters were allowed to float in the fits, and the posterior values were used to calculate the 39 Ar specific activity.The single 39 Ar events counted through the fit outputs of this cross-check method and that described in Section 5.1 do not differ significantly.An additional set of fits was performed using the Kostensalo et al. spectrum [28].The measured specific activity differed by a negligible amount between these fits and those using the nominal spectrum.
The result in Ref. [19] has been updated here to include the updated LAr mass, the new data-driven cut efficiency estimates, the revised livetime calculation, and the correction for the age of the argon.This method yields the following value for the specific activity of 39 Ar at the time of atmospheric extraction: (0.97 ± 0.001 stat ± 0.03 sys ) Bq/kg atmAr .

Conclusion
A measurement of the specific activity of 39 Ar in atmospheric argon using the LAr target of the DEAP-3600 detector has been presented.This result is the most precise measurement of the specific activity of 39 Ar in atmospheric argon to date and agrees with existing measurements which are summarized in Table 2.The high precision of this measurement is owing to a combination of factors including the low-background nature of DEAP-3600, the large number of decays observed in each data run, and the precise measurement of the LAr target mass.The statistical uncertainties here are much smaller than the systematic uncertainties due to the high statistics of the data.The dominant systematic uncertainties arise from the event selection cut efficiencies, the energy scale and energy resolution parameters, and the possible presence of 85 Kr within the LAr.
This precision measurement is an important input to the background models of experiments operating with argon as a medium.It will benefit current experiments, help to inform the design of future detectors, and support measurements in radiometric dating which use the 39 Ar/Ar ratio as an input.

Fig. 1
Fig.1A cross-section of the DEAP-3600 detector components located inside the water Cherenkov muon veto detector (not shown).

Table 1
Parameters, their values and constraints, and the resulting contributions to the uncertainty for the specific activity measurement.Negligibly small systematic uncertainties are indicated with '-'.The dominant uncertainty on S Ar39 arises from the uncertainties on event selection cut efficiency values as determined with the data-driven method (d-d) and the Monte Carlo method (MC).

Table 2
Summary of specific activity measurements of 39 Ar by different collaborations.