-
Physics of the low momentum diffusivity regime in tokamaks and its experimental applicability
Authors:
Haomin Sun,
Justin Ball,
Stephan Brunner,
Anthony Field,
Bhavin Patel,
Alessandro Balestri,
Daniel Kennedy,
Colin Roach,
Diego Jose Cruz-Zabala,
Fernando Puentes Del Pozo,
Eleonora Viezzer,
Manuel Garcia Munoz
Abstract:
Strong $E\times B$ plasma flow shear is beneficial for reducing turbulent transport. However, traditional methods of driving flow shear do not scale well to large devices such as future fusion power plants. In this paper, we use a large number of nonlinear gyrokinetic simulations to study a novel approach to increase flow shear: decreasing the momentum diffusivity to make the plasma ``easier to pu…
▽ More
Strong $E\times B$ plasma flow shear is beneficial for reducing turbulent transport. However, traditional methods of driving flow shear do not scale well to large devices such as future fusion power plants. In this paper, we use a large number of nonlinear gyrokinetic simulations to study a novel approach to increase flow shear: decreasing the momentum diffusivity to make the plasma ``easier to push''. We first use an idealized circular geometry and find that one can obtain low momentum diffusivity at tight aspect ratio, low safety factor, high magnetic shear and low temperature gradient. This is the so-called Low Momentum Diffusivity (LMD) regime. To drive intrinsic momentum flux, we then tilt the flux surface, making it up-down asymmetric. In the LMD regime, this intrinsic momentum flux drives strong flow shear that can significantly reduce the heat flux and increase the critical temperature gradient. We also consider the actual experimental geometry of the MAST tokamak to illustrate that this strategy can be practical and create experimentally significant flow shear. Lastly, a preliminary prediction for the SMART tokamak is made.
△ Less
Submitted 22 August, 2024;
originally announced August 2024.
-
Gyrokinetic flux-driven simulations in mixed TEM/ITG regime using a delta-f PIC scheme with evolving background
Authors:
Moahan Murugappan,
Laurent Villard,
Stephan Brunner,
Giovanni Di Giannatale,
Ben Fynney McMillan,
Alberto Bottino
Abstract:
In the context of global gyrokinetic simulations of turbulence using a Particle-In-Cell framework, verifying the delta-f assumption with a fixed background distribution becomes challenging when determining quasi-steady state profiles corresponding to given sources over long time scales, where plasma profiles can evolve significantly. The advantage of low relative sampling noise afforded by the del…
▽ More
In the context of global gyrokinetic simulations of turbulence using a Particle-In-Cell framework, verifying the delta-f assumption with a fixed background distribution becomes challenging when determining quasi-steady state profiles corresponding to given sources over long time scales, where plasma profiles can evolve significantly. The advantage of low relative sampling noise afforded by the delta-f scheme is shown to be retained by considering the background as a time-evolving Maxwellian with time-dependent density and temperature profiles. Implementation of this adaptive scheme to simulate electrostatic collisionless flux-driven turbulence in tokamak plasmas show small and non-increasing sampling noise levels, which would otherwise increase indefinitely with a stationary background scheme. The adaptive scheme furthermore allows one to reach numerically converged results of quasi-steady state with much lower marker numbers.
△ Less
Submitted 5 August, 2024; v1 submitted 24 June, 2024;
originally announced June 2024.
-
On the importance of slow ions in the kinetic Bohm criterion
Authors:
Alessandro Geraldini,
Stephan Brunner
Abstract:
Between a plasma and a solid target lies a positively charged sheath of several Debye lengths $λ_{\rm D}$ in width, typically much smaller than the characteristic length scale $L$ of the main plasma. This scale separation implies that the asymptotic limit $ε= λ_{\rm D} / L \rightarrow 0$ is useful to solve for the plasma-sheath system. In this limit, the Bohm criterion must be satisfied at the she…
▽ More
Between a plasma and a solid target lies a positively charged sheath of several Debye lengths $λ_{\rm D}$ in width, typically much smaller than the characteristic length scale $L$ of the main plasma. This scale separation implies that the asymptotic limit $ε= λ_{\rm D} / L \rightarrow 0$ is useful to solve for the plasma-sheath system. In this limit, the Bohm criterion must be satisfied at the sheath entrance. A new derivation of the kinetic criterion, admitting a general ion velocity distribution, is presented. It is proven that for $ε\rightarrow 0$ the distribution of the velocity component normal to the target, $v_x$, and its first derivative must vanish for $|v_x| \rightarrow 0$ at the sheath entrance. These two conditions can be subsumed into a third integral one after it is integrated by parts twice. A subsequent interchange of the limits $ε\rightarrow 0$ and $|v_x| \rightarrow 0$ is invalid, leading to a divergence which underlies the misconception that the criterion gives undue importance to slow ions.
△ Less
Submitted 14 August, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Self-interaction of turbulent eddies in tokamaks with low magnetic shear
Authors:
Arnas Volčokas,
Justin Ball,
Stephan Brunner
Abstract:
Using local nonlinear gyrokinetic simulations, we demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns in tokamaks with weak or zero magnetic shear $\hat{s}$. We observe that this parallel eddy length scales inversely with magnetic shear and at $\hat{s}=0$ is limited by the thermal speed of electrons $v_{th,e}$. We examine the consequences of these…
▽ More
Using local nonlinear gyrokinetic simulations, we demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns in tokamaks with weak or zero magnetic shear $\hat{s}$. We observe that this parallel eddy length scales inversely with magnetic shear and at $\hat{s}=0$ is limited by the thermal speed of electrons $v_{th,e}$. We examine the consequences of these "ultra long" eddies on turbulent transport, in particular, how field line topology mediates strong parallel self-interaction. Our investigation reveals that, through this process, field line topology can strongly affect transport. It can cause transitions between different turbulent instabilities and in some cases triple the logarithmic gradient needed to drive a given amount of heat flux. We also identify a novel "eddy squeezing" effect, which reduces the perpendicular size of eddies and their ability to transport energy, thus representing a novel approach to improve confinement. Finally, we investigate the triggering mechanism of Internal Transport Barriers (ITBs) using low magnetic shear simulations, shedding light on why ITBs are often easier to trigger where the safety factor has a low-order rational value.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
System size scaling of triangularity effects on global temperature gradient-driven gyrokinetic simulations
Authors:
Giovanni Di Giannatale,
Alberto Bottino,
Stephan Brunner,
Moahan Murugappan,
Laurent Villard
Abstract:
In this work, we explore the triangularity effects on turbulent transport employing global gyrokinetic simulations performed with the ORB5 code. Numerous experiments on the Tokamak à Configuration Variable (TCV) and, more recently, on the DIII-D machine, have demonstrated superior confinement properties in L-mode of negative triangularity (NT) over positive triangularity (PT) configuration. This p…
▽ More
In this work, we explore the triangularity effects on turbulent transport employing global gyrokinetic simulations performed with the ORB5 code. Numerous experiments on the Tokamak à Configuration Variable (TCV) and, more recently, on the DIII-D machine, have demonstrated superior confinement properties in L-mode of negative triangularity (NT) over positive triangularity (PT) configuration. This presents a particularly attractive scenario, as L-mode operation eliminates or significantly mitigates the presence of hazardous edge-localized modes (ELMs). However, a full theoretical understanding of all these observations remains elusive. Specifically, questions remain about how NT improvements can extend to the core where triangularity is very small, and whether these improvements can scale to larger devices. This paper addresses these two questions. Our analysis is divided into two parts: we first demonstrate that the confinement enhancement in NT configurations arises from the interdependent edge-core dynamics, and then we present the results of a system size scan. Crucially, we find that the relative turbulent transport reduction of NT over PT appears not to be contingent on machine dimensions or fluctuation scales and is moreover robust with respect to variations in plasma profiles. This insight underscores the fundamental nature of the NT confinement advantage and paves the way for its potential application in future fusion devices, regardless of their size.
△ Less
Submitted 25 February, 2024;
originally announced February 2024.
-
The XENONnT Dark Matter Experiment
Authors:
XENON Collaboration,
E. Aprile,
J. Aalbers,
K. Abe,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
M. Balata,
L. Baudis,
A. L. Baxter,
M. Bazyk,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui
, et al. (170 additional authors not shown)
Abstract:
The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in…
▽ More
The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in cryostat). The experiment is expected to extend the sensitivity to WIMP dark matter by more than an order of magnitude compared to XENON1T, thanks to the larger active mass and the significantly reduced background, improved by novel systems such as a radon removal plant and a neutron veto. This article describes the XENONnT experiment and its sub-systems in detail and reports on the detector performance during the first science run.
△ Less
Submitted 15 February, 2024;
originally announced February 2024.
-
Sheath constraints on turbulent magnetised plasmas
Authors:
Alessandro Geraldini,
Stephan Brunner,
Felix I. Parra
Abstract:
A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale $λ_{\rm D}$. In magnetised plasmas where the magnetic field is inclined a…
▽ More
A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale $λ_{\rm D}$. In magnetised plasmas where the magnetic field is inclined at an oblique angle $α$ with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale $ρ_{\rm S} \cos α$, caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness $l_{\rm ms} \sim \max (λ_{\rm D}, ρ_{\rm S} \cos α)$ is smaller than the characteristic length scale $L$ of spatial fluctuations in the bulk plasma, such that the limit $l_{\rm ms} / L \rightarrow 0$ is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle $ α\sim δ\equiv ρ_{\rm S} / L \ll 1$, electric fields tangential to the target transport ions towards the target via ExB drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm-Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.
△ Less
Submitted 24 June, 2024; v1 submitted 14 January, 2024;
originally announced January 2024.
-
A new quasilinear model for turbulent momentum transport in tokamaks with flow shear and plasma shaping
Authors:
Haomin Sun,
Justin Ball,
Stephan Brunner,
Arnas Volčokas
Abstract:
In tokamak experiments, sufficiently strong $E\times B$ flow shear reduces turbulent transport, thereby improving the prospects for fusion power plants. It is therefore of great importance to efficiently explore parameter space to find where strong plasma flow can be achieved. To this end, we propose a new, physically motivated quasi-linear model for estimating momentum transport from turbulence i…
▽ More
In tokamak experiments, sufficiently strong $E\times B$ flow shear reduces turbulent transport, thereby improving the prospects for fusion power plants. It is therefore of great importance to efficiently explore parameter space to find where strong plasma flow can be achieved. To this end, we propose a new, physically motivated quasi-linear model for estimating momentum transport from turbulence in the presence of toroidal flow shear and plasma shaping. The method gives good estimates of momentum transport for up-down asymmetric geometries as well as low magnetic shear and tight aspect ratio. The results are benchmarked with high-fidelity nonlinear GENE simulations, demonstrating that it provides a fast and accurate estimate of momentum transport.
△ Less
Submitted 29 January, 2024; v1 submitted 13 November, 2023;
originally announced November 2023.
-
Design and performance of the field cage for the XENONnT experiment
Authors:
E. Aprile,
K. Abe,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
M. Bazyk,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai,
J. M. R. Cardoso,
D. Cichon
, et al. (139 additional authors not shown)
Abstract:
The precision in reconstructing events detected in a dual-phase time projection chamber depends on an homogeneous and well understood electric field within the liquid target. In the XENONnT TPC the field homogeneity is achieved through a double-array field cage, consisting of two nested arrays of field shaping rings connected by an easily accessible resistor chain. Rather than being connected to t…
▽ More
The precision in reconstructing events detected in a dual-phase time projection chamber depends on an homogeneous and well understood electric field within the liquid target. In the XENONnT TPC the field homogeneity is achieved through a double-array field cage, consisting of two nested arrays of field shaping rings connected by an easily accessible resistor chain. Rather than being connected to the gate electrode, the topmost field shaping ring is independently biased, adding a degree of freedom to tune the electric field during operation. Two-dimensional finite element simulations were used to optimize the field cage, as well as its operation. Simulation results were compared to ${}^{83m}\mathrm{Kr}$ calibration data. This comparison indicates an accumulation of charge on the panels of the TPC which is constant over time, as no evolution of the reconstructed position distribution of events is observed. The simulated electric field was then used to correct the charge signal for the field dependence of the charge yield. This correction resolves the inconsistent measurement of the drift electron lifetime when using different calibrations sources and different field cage tuning voltages.
△ Less
Submitted 21 September, 2023;
originally announced September 2023.
-
Search for events in XENON1T associated with Gravitational Waves
Authors:
XENON Collaboration,
E. Aprile,
K. Abe,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antoń Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
M. Bazyk,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai,
J. M. R. Cardoso
, et al. (138 additional authors not shown)
Abstract:
We perform a blind search for particle signals in the XENON1T dark matter detector that occur close in time to gravitational wave signals in the LIGO and Virgo observatories. No particle signal is observed in the nuclear recoil, electronic recoil, CE$ν$NS, and S2-only channels within $\pm$ 500 seconds of observations of the gravitational wave signals GW170104, GW170729, GW170817, GW170818, and GW1…
▽ More
We perform a blind search for particle signals in the XENON1T dark matter detector that occur close in time to gravitational wave signals in the LIGO and Virgo observatories. No particle signal is observed in the nuclear recoil, electronic recoil, CE$ν$NS, and S2-only channels within $\pm$ 500 seconds of observations of the gravitational wave signals GW170104, GW170729, GW170817, GW170818, and GW170823. We use this null result to constrain mono-energetic neutrinos and Beyond Standard Model particles emitted in the closest coalescence GW170817, a binary neutron star merger. We set new upper limits on the fluence (time-integrated flux) of coincident neutrinos down to 17 keV at 90% confidence level. Furthermore, we constrain the product of coincident fluence and cross section of Beyond Standard Model particles to be less than $10^{-29}$ cm$^2$/cm$^2$ in the [5.5-210] keV energy range at 90% confidence level.
△ Less
Submitted 27 October, 2023; v1 submitted 20 June, 2023;
originally announced June 2023.
-
Probing hydrodynamic fluctuation-induced forces with an oscillating robot
Authors:
Steven W. Tarr,
Joseph S. Brunner,
Daniel Soto,
Daniel I. Goldman
Abstract:
We study the dynamics of an oscillating, free-floating robot that generates radially expanding gravity capillary waves at a fluid surface. In open water, the device does not self-propel; near a rigid boundary, it can be attracted or repelled. Visualization of the wave field dynamics reveals that when near a boundary, a complex interference of generated and reflected waves induces a wave amplitude…
▽ More
We study the dynamics of an oscillating, free-floating robot that generates radially expanding gravity capillary waves at a fluid surface. In open water, the device does not self-propel; near a rigid boundary, it can be attracted or repelled. Visualization of the wave field dynamics reveals that when near a boundary, a complex interference of generated and reflected waves induces a wave amplitude fluctuation asymmetry. Attraction increases as wave frequency increases or robot-boundary separation decreases. Theory on confined gravity-capillary wave radiation dynamics developed by Hocking in the 1980s captures the observed parameter dependence due to these "Hocking fields." The flexibility of the robophysical system allows detailed characterization and analysis of locally generated nonequilibrium fluctuation-induced forces [M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 (1999)].
△ Less
Submitted 20 February, 2024; v1 submitted 7 May, 2023;
originally announced May 2023.
-
Searching for Heavy Dark Matter near the Planck Mass with XENON1T
Authors:
E. Aprile,
K. Abe,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
M. Bazyk,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai,
J. M. R. Cardoso,
D. Cichon
, et al. (142 additional authors not shown)
Abstract:
Multiple viable theoretical models predict heavy dark matter particles with a mass close to the Planck mass, a range relatively unexplored by current experimental measurements. We use 219.4 days of data collected with the XENON1T experiment to conduct a blind search for signals from Multiply-Interacting Massive Particles (MIMPs). Their unique track signature allows a targeted analysis with only 0.…
▽ More
Multiple viable theoretical models predict heavy dark matter particles with a mass close to the Planck mass, a range relatively unexplored by current experimental measurements. We use 219.4 days of data collected with the XENON1T experiment to conduct a blind search for signals from Multiply-Interacting Massive Particles (MIMPs). Their unique track signature allows a targeted analysis with only 0.05 expected background events from muons. Following unblinding, we observe no signal candidate events. This work places strong constraints on spin-independent interactions of dark matter particles with a mass between 1$\times$10$^{12}\,$GeV/c$^2$ and 2$\times$10$^{17}\,$GeV/c$^2$. In addition, we present the first exclusion limits on spin-dependent MIMP-neutron and MIMP-proton cross-sections for dark matter particles with masses close to the Planck scale.
△ Less
Submitted 21 April, 2023;
originally announced April 2023.
-
First Dark Matter Search with Nuclear Recoils from the XENONnT Experiment
Authors:
XENON Collaboration,
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
M. Bazyk,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai
, et al. (141 additional authors not shown)
Abstract:
We report on the first search for nuclear recoils from dark matter in the form of weakly interacting massive particles (WIMPs) with the XENONnT experiment which is based on a two-phase time projection chamber with a sensitive liquid xenon mass of $5.9$ t. During the approximately 1.1 tonne-year exposure used for this search, the intrinsic $^{85}$Kr and $^{222}$Rn concentrations in the liquid targe…
▽ More
We report on the first search for nuclear recoils from dark matter in the form of weakly interacting massive particles (WIMPs) with the XENONnT experiment which is based on a two-phase time projection chamber with a sensitive liquid xenon mass of $5.9$ t. During the approximately 1.1 tonne-year exposure used for this search, the intrinsic $^{85}$Kr and $^{222}$Rn concentrations in the liquid target were reduced to unprecedentedly low levels, giving an electronic recoil background rate of $(15.8\pm1.3)~\mathrm{events}/(\mathrm{t\cdot y \cdot keV})$ in the region of interest. A blind analysis of nuclear recoil events with energies between $3.3$ keV and $60.5$ keV finds no significant excess. This leads to a minimum upper limit on the spin-independent WIMP-nucleon cross section of $2.58\times 10^{-47}~\mathrm{cm}^2$ for a WIMP mass of $28~\mathrm{GeV}/c^2$ at $90\%$ confidence level. Limits for spin-dependent interactions are also provided. Both the limit and the sensitivity for the full range of WIMP masses analyzed here improve on previous results obtained with the XENON1T experiment for the same exposure.
△ Less
Submitted 5 August, 2023; v1 submitted 26 March, 2023;
originally announced March 2023.
-
Verification of the Fourier-enhanced 3D finite element Poisson solver of the gyrokinetic full-f code PICLS
Authors:
Annika Stier,
Alberto Bottino,
Mathias Boesl,
Martin Campos Pinto,
Thomas Hayward-Schneider,
David Coster,
Andreas Bergmann,
Moahan Murugappan,
Stephan Brunner,
Laurent Villard,
Frank Jenko
Abstract:
We introduce and derive the Fourier-enhanced 3D electrostatic field solver of the gyrokinetic full-f PIC code PICLS. The solver makes use of a Fourier representation in one periodic direction of the domain to make the solving of the system easily parallelizable and thus save run time. The presented solver is then verified using two different approaches of manufactured solutions. The test setup use…
▽ More
We introduce and derive the Fourier-enhanced 3D electrostatic field solver of the gyrokinetic full-f PIC code PICLS. The solver makes use of a Fourier representation in one periodic direction of the domain to make the solving of the system easily parallelizable and thus save run time. The presented solver is then verified using two different approaches of manufactured solutions. The test setup used for this effort is a pinch geometry with ITG-like electric potential, containing one non-periodic and two periodic directions, one of which will be discrete Fourier transformed. The results of these tests show that in all three dimensions the L2-error decreases with a constant rate close to the ideal prediction, depending on the degree of the chosen basis functions.
△ Less
Submitted 2 March, 2023;
originally announced March 2023.
-
The Triggerless Data Acquisition System of the XENONnT Experiment
Authors:
E. Aprile,
J. Aalbers,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
R. Biondi,
A. Bismark,
E. J. Brookes,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai,
J. M. R. Cardoso
, et al. (140 additional authors not shown)
Abstract:
The XENONnT detector uses the latest and largest liquid xenon-based time projection chamber (TPC) operated by the XENON Collaboration, aimed at detecting Weakly Interacting Massive Particles and conducting other rare event searches. The XENONnT data acquisition (DAQ) system constitutes an upgraded and expanded version of the XENON1T DAQ system. For its operation, it relies predominantly on commerc…
▽ More
The XENONnT detector uses the latest and largest liquid xenon-based time projection chamber (TPC) operated by the XENON Collaboration, aimed at detecting Weakly Interacting Massive Particles and conducting other rare event searches. The XENONnT data acquisition (DAQ) system constitutes an upgraded and expanded version of the XENON1T DAQ system. For its operation, it relies predominantly on commercially available hardware accompanied by open-source and custom-developed software. The three constituent subsystems of the XENONnT detector, the TPC (main detector), muon veto, and the newly introduced neutron veto, are integrated into a single DAQ, and can be operated both independently and as a unified system. In total, the DAQ digitizes the signals of 698 photomultiplier tubes (PMTs), of which 253 from the top PMT array of the TPC are digitized twice, at $\times10$ and $\times0.5$ gain. The DAQ for the most part is a triggerless system, reading out and storing every signal that exceeds the digitization thresholds. Custom-developed software is used to process the acquired data, making it available within $\mathcal{O}\left(10\text{ s}\right)$ for live data quality monitoring and online analyses. The entire system with all the three subsystems was successfully commissioned and has been operating continuously, comfortably withstanding readout rates that exceed $\sim500$ MB/s during calibration. Livetime during normal operation exceeds $99\%$ and is $\sim90\%$ during most high-rate calibrations. The combined DAQ system has collected more than 2 PB of both calibration and science data during the commissioning of XENONnT and the first science run.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Low-energy Calibration of XENON1T with an Internal $^{37}$Ar Source
Authors:
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
R. Biondi,
A. Bismark,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
T. K. Bui,
C. Cai,
C. Capelli,
J. M. R. Cardoso
, et al. (139 additional authors not shown)
Abstract:
A low-energy electronic recoil calibration of XENON1T, a dual-phase xenon time projection chamber, with an internal $^{37}$Ar source was performed. This calibration source features a 35-day half-life and provides two mono-energetic lines at 2.82 keV and 0.27 keV. The photon yield and electron yield at 2.82 keV are measured to be (32.3$\pm$0.3) photons/keV and (40.6$\pm$0.5) electrons/keV, respecti…
▽ More
A low-energy electronic recoil calibration of XENON1T, a dual-phase xenon time projection chamber, with an internal $^{37}$Ar source was performed. This calibration source features a 35-day half-life and provides two mono-energetic lines at 2.82 keV and 0.27 keV. The photon yield and electron yield at 2.82 keV are measured to be (32.3$\pm$0.3) photons/keV and (40.6$\pm$0.5) electrons/keV, respectively, in agreement with other measurements and with NEST predictions. The electron yield at 0.27 keV is also measured and it is (68.0$^{+6.3}_{-3.7}$) electrons/keV. The $^{37}$Ar calibration confirms that the detector is well-understood in the energy region close to the detection threshold, with the 2.82 keV line reconstructed at (2.83$\pm$0.02) keV, which further validates the model used to interpret the low-energy electronic recoil excess previously reported by XENON1T. The ability to efficiently remove argon with cryogenic distillation after the calibration proves that $^{37}$Ar can be considered as a regular calibration source for multi-tonne xenon detectors.
△ Less
Submitted 21 March, 2023; v1 submitted 25 November, 2022;
originally announced November 2022.
-
An approximate likelihood for nuclear recoil searches with XENON1T data
Authors:
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
B. Andrieu,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
R. Biondi,
A. Bismark,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
B. Cimmino
, et al. (129 additional authors not shown)
Abstract:
The XENON collaboration has published stringent limits on specific dark matter -nucleon recoil spectra from dark matter recoiling on the liquid xenon detector target. In this paper, we present an approximate likelihood for the XENON1T 1 tonne-year nuclear recoil search applicable to any nuclear recoil spectrum. Alongside this paper, we publish data and code to compute upper limits using the method…
▽ More
The XENON collaboration has published stringent limits on specific dark matter -nucleon recoil spectra from dark matter recoiling on the liquid xenon detector target. In this paper, we present an approximate likelihood for the XENON1T 1 tonne-year nuclear recoil search applicable to any nuclear recoil spectrum. Alongside this paper, we publish data and code to compute upper limits using the method we present. The approximate likelihood is constructed in bins of reconstructed energy, profiled along the signal expectation in each bin. This approach can be used to compute an approximate likelihood and therefore most statistical results for any nuclear recoil spectrum. Computing approximate results with this method is approximately three orders of magnitude faster than the likelihood used in the original publications of XENON1T, where limits were set for specific families of recoil spectra. Using this same method, we include toy Monte Carlo simulation-derived binwise likelihoods for the upcoming XENONnT experiment that can similarly be used to assess the sensitivity to arbitrary nuclear recoil signatures in its eventual 20 tonne-year exposure.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
Moment-Based Approach to the Flux-Tube linear Gyrokinetic Model
Authors:
B. J. Frei,
A. C. D. Hoffmann,
P. Ricci,
S. Brunner,
Z. Tecchiolli
Abstract:
This work reports on the development and numerical implementation of the linear electromagnetic gyrokinetic (GK) model in a tokamak flux-tube geometry using a moment approach based on the expansion of the perturbed distribution function on a velocity-space Hermite-Laguerre polynomials basis. A hierarchy of equations of the expansion coefficients, referred to as the gyro-moments (GM), is derived. W…
▽ More
This work reports on the development and numerical implementation of the linear electromagnetic gyrokinetic (GK) model in a tokamak flux-tube geometry using a moment approach based on the expansion of the perturbed distribution function on a velocity-space Hermite-Laguerre polynomials basis. A hierarchy of equations of the expansion coefficients, referred to as the gyro-moments (GM), is derived. We verify the numerical implementation of the GM hierarchy in the collisionless limit by performing a comparison with the continuum GK code GENE, recovering the linear properties of the ion-temperature gradient, trapped electron, kinetic ballooning, and microtearing modes, as well as the collisionless damping of zonal flows. The present investigation reveals the ability of the GM approach to describe fine velocity-space scale structures appearing near the trapped and passing boundary and kinetic effects associated with parallel and perpendicular particle drifts. In addition, the effects of collisions are studied using advanced collision operators, including the GK Coulomb collision operator. The main findings are that the number of GMs necessary for convergence decreases with plasma collisionality and is lower for pressure gradient-driven modes, such as in H-mode pedestal regions, compared to instabilities driven by trapped particles and magnetic gradient drifts often found in the core. The accuracy of approximations often used to model collisions (relative to the GK Coulomb operator) is studied, showing differences between collision operator models that increase with collisionality and electron temperature gradient in the case of TEM. The present linear analysis demonstrates that the GM approach efficiently describes the plasma dynamics for typical parameters of the tokamak boundary, ranging from the low-collisionality banana H-mode to the high-collisionality Pfirsch-Schlüter conditions.
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
Gyrokinetic simulations of turbulence and zonal flows driven by steep profile gradients using a delta-f approach with an evolving background Maxwellian
Authors:
Moahan Murugappan,
Laurent Villard,
Stephan Brunner,
Ben Fynney McMillan,
Alberto Bottino
Abstract:
Long global gyrokinetic turbulence simulations are particularly challenging in situations where the system deviates strongly from its initial state and when fluctuation level are high e.g. in strong gradient regions. For Particle-in-Cell simulations, statistical sampling noise accumulation from large marker weights due to large deviations from the control variate of a delta-f scheme make such simu…
▽ More
Long global gyrokinetic turbulence simulations are particularly challenging in situations where the system deviates strongly from its initial state and when fluctuation level are high e.g. in strong gradient regions. For Particle-in-Cell simulations, statistical sampling noise accumulation from large marker weights due to large deviations from the control variate of a delta-f scheme make such simulations often impractical. An adaptive control variate in the form of a flux-surface-averaged Maxwellian with a time-dependent temperature profile is introduced in an attempt to alleviate the former problem. Under simplified collisionless physics, this adaptive delta-f scheme is shown to reduce noise accumulation in the zonal flows and the simulated heat flux in a quasi-steady turbulent state. The method also avoids the collapse of the signal-to-noise ratio which occurs in the standard non-adaptive scheme, and therefore, allows one to reach numerically converged results even with lower marker numbers.
△ Less
Submitted 31 August, 2022;
originally announced August 2022.
-
Local gyrokinetic simulations of tokamaks with non-uniform magnetic shear
Authors:
Justin Ball,
Stephan Brunner
Abstract:
In this work, we modify the standard flux tube simulation domain to include arbitrary ion gyroradius-scale variation in the radial profile of the safety factor. To determine how to appropriately include such a modification, we add a strong ion gyroradius-scale source (inspired by electron cyclotron current drive) to the Fokker-Planck equation, then perform a multi-scale analysis that distinguishes…
▽ More
In this work, we modify the standard flux tube simulation domain to include arbitrary ion gyroradius-scale variation in the radial profile of the safety factor. To determine how to appropriately include such a modification, we add a strong ion gyroradius-scale source (inspired by electron cyclotron current drive) to the Fokker-Planck equation, then perform a multi-scale analysis that distinguishes the fast electrons driven by the source from the slow bulk thermal electrons. This allows us to systematically derive the needed changes to the gyrokinetic model. We find new terms that adjust the ion and electron parallel streaming to be along the modified field lines. These terms have been successfully implemented in a gyrokinetic code (while retaining the typical Fourier representation), which enables flux tube studies of non-monotonic safety factor profiles and the associated profile shearing. As an illustrative example, we investigate tokamaks with positive versus negative triangularity plasma shaping and find that the importance of profile shearing is not significantly affected by the change in shape.
△ Less
Submitted 13 December, 2022; v1 submitted 25 August, 2022;
originally announced August 2022.
-
Ultra Long Turbulent Eddies, Magnetic Topology, and the Triggering of Internal Transport Barriers in Tokamaks
Authors:
A. Volčokas,
J. Ball,
S. Brunner
Abstract:
Local nonlinear gyrokinetic simulations of tokamak plasmas demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear is very small. By accurately modeling different field line topologies (e.g. low-order, almost rational, or irrational safety factor), we show that the parallel self-interaction of such "ultra long" eddies can drama…
▽ More
Local nonlinear gyrokinetic simulations of tokamak plasmas demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear is very small. By accurately modeling different field line topologies (e.g. low-order, almost rational, or irrational safety factor), we show that the parallel self-interaction of such "ultra long" eddies can dramatically reduce heat transport. This reveals novel strategies to improve confinement, constitutes experimentally testable predictions, and illuminates past observations of internal transport barriers.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics
Authors:
J. Aalbers,
K. Abe,
V. Aerne,
F. Agostini,
S. Ahmed Maouloud,
D. S. Akerib,
D. Yu. Akimov,
J. Akshat,
A. K. Al Musalhi,
F. Alder,
S. K. Alsum,
L. Althueser,
C. S. Amarasinghe,
F. D. Amaro,
A. Ames,
T. J. Anderson,
B. Andrieu,
N. Angelides,
E. Angelino,
J. Angevaare,
V. C. Antochi,
D. Antón Martin,
B. Antunovic,
E. Aprile,
H. M. Araújo
, et al. (572 additional authors not shown)
Abstract:
The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut…
▽ More
The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
Application and modeling of an online distillation method to reduce krypton and argon in XENON1T
Authors:
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
A. Bernard,
R. Biondi,
A. Bismark,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
B. Cimmino
, et al. (129 additional authors not shown)
Abstract:
A novel online distillation technique was developed for the XENON1T dark matter experiment to reduce intrinsic background components more volatile than xenon, such as krypton or argon, while the detector was operating. The method is based on a continuous purification of the gaseous volume of the detector system using the XENON1T cryogenic distillation column. A krypton-in-xenon concentration of…
▽ More
A novel online distillation technique was developed for the XENON1T dark matter experiment to reduce intrinsic background components more volatile than xenon, such as krypton or argon, while the detector was operating. The method is based on a continuous purification of the gaseous volume of the detector system using the XENON1T cryogenic distillation column. A krypton-in-xenon concentration of $(360 \pm 60)$ ppq was achieved. It is the lowest concentration measured in the fiducial volume of an operating dark matter detector to date. A model was developed and fit to the data to describe the krypton evolution in the liquid and gas volumes of the detector system for several operation modes over the time span of 550 days, including the commissioning and science runs of XENON1T. The online distillation was also successfully applied to remove Ar-37 after its injection for a low energy calibration in XENON1T. This makes the usage of Ar-37 as a regular calibration source possible in the future. The online distillation can be applied to next-generation experiments to remove krypton prior to, or during, any science run. The model developed here allows further optimization of the distillation strategy for future large scale detectors.
△ Less
Submitted 14 June, 2022; v1 submitted 22 December, 2021;
originally announced December 2021.
-
Emission of Single and Few Electrons in XENON1T and Limits on Light Dark Matter
Authors:
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
A. Bernard,
R. Biondi,
A. Bismark,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
B. Cimmino
, et al. (130 additional authors not shown)
Abstract:
Delayed single- and few-electron emissions plague dual-phase time projection chambers, limiting their potential to search for light-mass dark matter. This paper examines the origins of these events in the XENON1T experiment. Characterization of the intensity of delayed electron backgrounds shows that the resulting emissions are correlated, in time and position, with high-energy events and can effe…
▽ More
Delayed single- and few-electron emissions plague dual-phase time projection chambers, limiting their potential to search for light-mass dark matter. This paper examines the origins of these events in the XENON1T experiment. Characterization of the intensity of delayed electron backgrounds shows that the resulting emissions are correlated, in time and position, with high-energy events and can effectively be vetoed. In this work we extend previous S2-only analyses down to a single electron. From this analysis, after removing the correlated backgrounds, we observe rates < 30 events/(electron*kg*day) in the region of interest spanning 1 to 5 electrons. We derive 90% confidence upper limits for dark matter-electron scattering, first direct limits on the electric dipole, magnetic dipole, and anapole interactions, and bosonic dark matter models, where we exclude new parameter space for dark photons and solar dark photons.
△ Less
Submitted 2 September, 2024; v1 submitted 22 December, 2021;
originally announced December 2021.
-
Material radiopurity control in the XENONnT experiment
Authors:
E. Aprile,
K. Abe,
F. Agostini,
S. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
E. Angelino,
J. R. Angevaare,
V. C. Antochi,
D. Antón Martin,
F. Arneodo,
L. Baudis,
A. L. Baxter,
L. Bellagamba,
R. Biondi,
A. Bismark,
A. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
B. Cimmino,
M. Clark
, et al. (128 additional authors not shown)
Abstract:
The selection of low-radioactive construction materials is of the utmost importance for rare-event searches and thus critical to the XENONnT experiment. Results of an extensive radioassay program are reported, in which material samples have been screened with gamma-ray spectroscopy, mass spectrometry, and $^{222}$Rn emanation measurements. Furthermore, the cleanliness procedures applied to remove…
▽ More
The selection of low-radioactive construction materials is of the utmost importance for rare-event searches and thus critical to the XENONnT experiment. Results of an extensive radioassay program are reported, in which material samples have been screened with gamma-ray spectroscopy, mass spectrometry, and $^{222}$Rn emanation measurements. Furthermore, the cleanliness procedures applied to remove or mitigate surface contamination of detector materials are described. Screening results, used as inputs for a XENONnT Monte Carlo simulation, predict a reduction of materials background ($\sim$17%) with respect to its predecessor XENON1T. Through radon emanation measurements, the expected $^{222}$Rn activity concentration in XENONnT is determined to be 4.2$\,(^{+0.5}_{-0.7})\,μ$Bq/kg, a factor three lower with respect to XENON1T. This radon concentration will be further suppressed by means of the novel radon distillation system.
△ Less
Submitted 26 January, 2023; v1 submitted 10 December, 2021;
originally announced December 2021.
-
Effect of collisions on non-adiabatic electron dynamics in ITG-driven microturbulence
Authors:
Ajay C. J.,
Stephan Brunner,
Justin Ball
Abstract:
Non-adiabatic electron response leads to significant changes in Ion Temperature Gradient (ITG) eigenmodes, leading in particular to fine-structures that are significantly extended along the magnetic field lines at corresponding Mode Rational Surfaces (MRSs). These eigenmodes can nonlinearly interact with themselves to drive zonal flows via the so-called self-interaction mechanism. In this paper, t…
▽ More
Non-adiabatic electron response leads to significant changes in Ion Temperature Gradient (ITG) eigenmodes, leading in particular to fine-structures that are significantly extended along the magnetic field lines at corresponding Mode Rational Surfaces (MRSs). These eigenmodes can nonlinearly interact with themselves to drive zonal flows via the so-called self-interaction mechanism. In this paper, the effect of collisions on these processes are studied. In presence of non-adiabatic electrons, the linear growth rate of ITG eigenmodes decreases with increasing collisionality. Detailed velocity space analysis of the distribution function shows that this results from collisions leading to a more adiabatic-like response of electrons away from MRSs. In linear simulations, collisions are furthermore found to broaden the radial width of the fine-structures, which translates to narrower tails of the eigenmode in extended ballooning space. The characteristic parallel scale length associated to these tails is shown to scale with the mean free path of electrons. In nonlinear turbulence simulations accounting for physically relevant values of collisionality, the fine structures located at MRSs, together with the associated drive of zonal flows via self-interaction, are shown to persist and play a significant role.
△ Less
Submitted 27 April, 2021; v1 submitted 26 April, 2021;
originally announced April 2021.
-
Nonlocal effects in negative triangularity TCV plasmas
Authors:
G. Merlo,
Z. Huang,
C. Marini,
S. Brunner,
S. Coda,
D. Hatch,
D. Jarema,
F. Jenko,
O. Sauter,
L. Villard
Abstract:
Global gradient driven GENE gyrokinetic simulations are used to investigate TCV plasmas with negative triangularity. Considering a limited L-mode plasma, corresponding to an experimental triangularity scan, numerical results are able to reproduce the actual transport level over a major fraction of the plasma minor radius for a plasma with $δ_{\rm LCFS}=-0.3$ and its equivalent with standard positi…
▽ More
Global gradient driven GENE gyrokinetic simulations are used to investigate TCV plasmas with negative triangularity. Considering a limited L-mode plasma, corresponding to an experimental triangularity scan, numerical results are able to reproduce the actual transport level over a major fraction of the plasma minor radius for a plasma with $δ_{\rm LCFS}=-0.3$ and its equivalent with standard positive triangularity $δ$. For the same heat flux, a larger electron temperature gradient is sustained by $δ<0$, in turn resulting in an improved electron energy confinement. Consistently with the experiments, a reduction of the electron density fluctuations is also seen. Local flux-tube simulations are used to gauge the magnitude of nonlocal effects. Surprisingly, very little differences are found between local and global approaches for $δ>0$, while local results yield a strong overestimation of the heat fluxes when $δ<0$. Despite the high sensitivity of the turbulence level with respect to the input parameters, global effects appear to play a crucial role in the negative triangularity plasma and must be retained to reconcile simulations and experiments. Finally, a general stabilizing effect of negative triangularity, reducing fluxes and fluctuations by a factor dependent on the actual profiles, is recovered.
△ Less
Submitted 14 January, 2021;
originally announced January 2021.
-
A non-twisting flux tube for local gyrokinetic simulations
Authors:
Justin Ball,
Stephan Brunner
Abstract:
Local gyrokinetic simulations use a field-aligned domain that twists due to the magnetic shear of the background magnetic equilibrium. However, if the magnetic shear is strong and/or the domain is long, the twist can become so extreme that it fails to properly resolve the turbulence. In this work, we derive and implement the "non-twisting flux tube," a local simulation domain that remains rectangu…
▽ More
Local gyrokinetic simulations use a field-aligned domain that twists due to the magnetic shear of the background magnetic equilibrium. However, if the magnetic shear is strong and/or the domain is long, the twist can become so extreme that it fails to properly resolve the turbulence. In this work, we derive and implement the "non-twisting flux tube," a local simulation domain that remains rectangular at all parallel locations. Convergence and runtime tests indicate that it can calculate the heat flux more efficiently than the conventional flux tube. For one test case, it was 30 times less computationally expensive and we found no case for which it was more expensive. It is most advantageous when the magnetic shear is high and the domain includes at least two regions of turbulent drive (e.g. stellarator simulations, pedestal simulations, tokamak simulations with several poloidal turns). Additionally, it more accurately models the inboard midplane when the magnetic shear is large. Lastly, we show how the non-twisting flux tube can be generalized to allow further optimization and control of the simulation domain.
△ Less
Submitted 6 May, 2021; v1 submitted 8 December, 2020;
originally announced December 2020.
-
$^{222}$Rn emanation measurements for the XENON1T experiment
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
E. Angelino,
J. R. Angevaare,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon
, et al. (118 additional authors not shown)
Abstract:
The selection of low-radioactive construction materials is of utmost importance for the success of low-energy rare event search experiments. Besides radioactive contaminants in the bulk, the emanation of radioactive radon atoms from material surfaces attains increasing relevance in the effort to further reduce the background of such experiments. In this work, we present the $^{222}$Rn emanation me…
▽ More
The selection of low-radioactive construction materials is of utmost importance for the success of low-energy rare event search experiments. Besides radioactive contaminants in the bulk, the emanation of radioactive radon atoms from material surfaces attains increasing relevance in the effort to further reduce the background of such experiments. In this work, we present the $^{222}$Rn emanation measurements performed for the XENON1T dark matter experiment. Together with the bulk impurity screening campaign, the results enabled us to select the radio-purest construction materials, targeting a $^{222}$Rn activity concentration of 10 $μ$Bq/kg in 3.2 t of xenon. The knowledge of the distribution of the $^{222}$Rn sources allowed us to selectively eliminate critical components in the course of the experiment. The predictions from the emanation measurements were compared to data of the $^{222}$Rn activity concentration in XENON1T. The final $^{222}$Rn activity concentration of (4.5 $\pm$ 0.1) $μ$Bq/kg in the target of XENON1T is the lowest ever achieved in a xenon dark matter experiment.
△ Less
Submitted 25 November, 2020; v1 submitted 29 September, 2020;
originally announced September 2020.
-
Radon daughter removal from PTFE surfaces and its application in liquid xenon detectors
Authors:
S. Bruenner,
D. Cichon,
G. Eurin,
P. Herrero Gómez,
F. Jörg,
T. Marrodán Undagoitia,
H. Simgen,
N. Rupp
Abstract:
Long-lived radon daughters are a critical background source in experiments searching for low-energy rare events. Originating from radon in ambient air, radioactive polonium, bismuth and lead isotopes plate-out on materials that are later employed in the experiment. In this paper, we examine cleaning procedures for their capability to remove radon daughters from PTFE surfaces, a material often used…
▽ More
Long-lived radon daughters are a critical background source in experiments searching for low-energy rare events. Originating from radon in ambient air, radioactive polonium, bismuth and lead isotopes plate-out on materials that are later employed in the experiment. In this paper, we examine cleaning procedures for their capability to remove radon daughters from PTFE surfaces, a material often used in liquid xenon TPCs. We found a large difference between the removal efficiency obtained for the decay chains of $^{222}$Rn and $^{220}$Rn, respectively. This indicates that the plate-out mechanism has an effect on the cleaning success. While the long-lived $^{222}$Rn daughters could be reduced by a factor of ~2, the removal of $^{220}$Rn daughters was up to 10 times more efficient depending on the treatment. Furthermore, the impact of a nitric acid based PTFE cleaning on the liquid xenon purity is investigated in a small-scale liquid xenon TPC.
△ Less
Submitted 18 September, 2020;
originally announced September 2020.
-
Projected WIMP Sensitivity of the XENONnT Dark Matter Experiment
Authors:
The XENON collaboration,
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
E. Angelino,
J. R. Angevaare,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon
, et al. (115 additional authors not shown)
Abstract:
XENONnT is a dark matter direct detection experiment, utilizing 5.9 t of instrumented liquid xenon, located at the INFN Laboratori Nazionali del Gran Sasso. In this work, we predict the experimental background and project the sensitivity of XENONnT to the detection of weakly interacting massive particles (WIMPs). The expected average differential background rate in the energy region of interest, c…
▽ More
XENONnT is a dark matter direct detection experiment, utilizing 5.9 t of instrumented liquid xenon, located at the INFN Laboratori Nazionali del Gran Sasso. In this work, we predict the experimental background and project the sensitivity of XENONnT to the detection of weakly interacting massive particles (WIMPs). The expected average differential background rate in the energy region of interest, corresponding to (1, 13) keV and (4, 50) keV for electronic and nuclear recoils, amounts to $12.3 \pm 0.6$ (keV t y)$^{-1}$ and $(2.2\pm 0.5)\times 10^{-3}$ (keV t y)$^{-1}$, respectively, in a 4 t fiducial mass. We compute unified confidence intervals using the profile construction method, in order to ensure proper coverage. With the exposure goal of 20 t$\,$y, the expected sensitivity to spin-independent WIMP-nucleon interactions reaches a cross-section of $1.4\times10^{-48}$ cm$^2$ for a 50 GeV/c$^2$ mass WIMP at 90% confidence level, more than one order of magnitude beyond the current best limit, set by XENON1T. In addition, we show that for a 50 GeV/c$^2$ WIMP with cross-sections above $2.6\times10^{-48}$ cm$^2$ ($5.0\times10^{-48}$ cm$^2$) the median XENONnT discovery significance exceeds 3$σ$ (5$σ$). The expected sensitivity to the spin-dependent WIMP coupling to neutrons (protons) reaches $2.2\times10^{-43}$ cm$^2$ ($6.0\times10^{-42}$ cm$^2$).
△ Less
Submitted 17 November, 2020; v1 submitted 17 July, 2020;
originally announced July 2020.
-
Monte Carlo Simulation Variance Reduction Techniques for Photon Transport in Liquid Xenon Detectors
Authors:
S. Bruenner,
A. P. Colijn,
M. P. Decowski,
O. V. Kesber
Abstract:
Monte Carlo simulations are a crucial tool for the analysis and prediction of various background components in liquid xenon (LXe) detectors. With improving shielding in new experiments, the simulation of external backgrounds, such as induced by gamma rays from detector materials, gets more computationally expensive. We introduce and validate an accelerated Monte Carlo simulation technique for phot…
▽ More
Monte Carlo simulations are a crucial tool for the analysis and prediction of various background components in liquid xenon (LXe) detectors. With improving shielding in new experiments, the simulation of external backgrounds, such as induced by gamma rays from detector materials, gets more computationally expensive. We introduce and validate an accelerated Monte Carlo simulation technique for photon transport in liquid xenon detectors. The method simulates photon-induced interactions within a defined geometry and energy range with high statistics while interactions outside of the region of interest are not simulated directly but are taken into account by means of probability weights. For a simulation of gamma induced backgrounds in an exemplary detector geometry we achieve a three orders of magnitude acceleration compared to a standard simulation of a current ton-scale LXe dark matter experiment.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Solar Neutrino Detection Sensitivity in DARWIN via Electron Scattering
Authors:
J. Aalbers,
F. Agostini,
S. E. M. Ahmed Maouloud,
M. Alfonsi,
L. Althueser,
F. Amaro,
J. Angevaare,
V. C. Antochi,
B. Antunovic,
E. Aprile,
L. Arazi,
F. Arneodo,
M. Balzer,
L. Baudis,
D. Baur,
M. L. Benabderrahmane,
Y. Biondi,
A. Bismark,
C. Bourgeois,
A. Breskin,
P. A. Breur,
A. Brown,
E. Brown,
S. Brünner,
G. Bruno
, et al. (141 additional authors not shown)
Abstract:
We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would ben…
▽ More
We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2θ_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$σ$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.
△ Less
Submitted 20 December, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size
Authors:
Ajay C. J.,
Stephan Brunner,
Ben McMillan,
Justin Ball,
Julien Dominski,
Gabriele Merlo
Abstract:
Self-interaction is the process by which a microturbulence eigenmode that is extended along the direction parallel to the magnetic field interacts with itself non-linearly. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics. Self-interaction is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near…
▽ More
Self-interaction is the process by which a microturbulence eigenmode that is extended along the direction parallel to the magnetic field interacts with itself non-linearly. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics. Self-interaction is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low order mode rational surfaces [Weikl et. al., Phys. Plasmas 25, 072305 (2018)]. We find however that it also plays a significant role in generating fluctuating zonal flows, which is critical to regulating transport throughout the radial extent. Unlike the usual picture of zonal flow drive where microturbulence eigenmodes coherently amplify the flow, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In flux-tube simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.
△ Less
Submitted 2 August, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
Energy resolution and linearity of XENON1T in the MeV energy range
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
E. Angelino,
J. Angevaare,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon
, et al. (113 additional authors not shown)
Abstract:
Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XEN…
▽ More
Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $σ/μ$ is as low as (0.80$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{ββ}$. We also present a new signal correction method to rectify the saturation effects of the signal readout system, resulting in more accurate position reconstruction and indirectly improving the energy resolution. The very good result achieved in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.
△ Less
Submitted 9 September, 2020; v1 submitted 8 March, 2020;
originally announced March 2020.
-
Eliminating turbulent self-interaction through the parallel boundary condition in local gyrokinetic simulations
Authors:
Justin Ball,
Stephan Brunner,
Ajay C. J
Abstract:
In this work, we highlight an issue that may reduce the accuracy of many local nonlinear gyrokinetic simulations - turbulent self-interaction through the parallel boundary condition. Given a sufficiently long parallel correlation length, individual turbulent eddies can span the full domain and "bite their own tails," thereby altering their statistical properties. Such self-interaction is only mode…
▽ More
In this work, we highlight an issue that may reduce the accuracy of many local nonlinear gyrokinetic simulations - turbulent self-interaction through the parallel boundary condition. Given a sufficiently long parallel correlation length, individual turbulent eddies can span the full domain and "bite their own tails," thereby altering their statistical properties. Such self-interaction is only modeled accurately when the simulation domain corresponds to a full flux surface, otherwise it is artificially strong. For Cyclone Base Case parameters and typical domain sizes, we find that this mechanism modifies the heat flux by roughly 40% and it can be even more important. The effect is largest when using kinetic electrons, low magnetic shear, and strong turbulence drive (i.e. steep background gradients). It is found that parallel self-interaction can be eliminated by increasing the parallel length and/or the binormal width of the simulation domain until convergence is achieved.
△ Less
Submitted 13 March, 2020; v1 submitted 19 December, 2019;
originally announced December 2019.
-
Collisional gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS
Authors:
Mathias Boesl,
Andreas Bergmann,
Alberto Bottino,
Stephan Brunner,
David Coster,
Frank Jenko
Abstract:
Applying gyrokinetic simulations for theoretical turbulence and transport studies to the plasma edge and scrape-off layer (SOL) presents significant challenges. To in particular account for steep density and temperature gradients in the SOL, the "full-f" code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code and is based on an electrostatic model with a linearized field equat…
▽ More
Applying gyrokinetic simulations for theoretical turbulence and transport studies to the plasma edge and scrape-off layer (SOL) presents significant challenges. To in particular account for steep density and temperature gradients in the SOL, the "full-f" code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code and is based on an electrostatic model with a linearized field equation and uses kinetic electrons. In previously published results we were applying PICLS to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) in the SOL without collisions. As an extension to this collision-less case and in preparation for 3D simulations, in this work a collisional model will be introduced. The implemented Lenard-Bernstein collision operator and its Langevin discretization will be shown. Conservation properties of the collision operator as well as a comparison of the collisional and non-collisional case will be discussed.
△ Less
Submitted 9 September, 2019;
originally announced September 2019.
-
Gyrokinetic Simulations on Many- and Multi-core Architectures with the Global Electromagnetic Particle-In-Cell Code ORB5
Authors:
Noé Ohana,
Claudio Gheller,
Emmanuel Lanti,
Andreas Jocksch,
Stephan Brunner,
Laurent Villard
Abstract:
Gyrokinetic codes in plasma physics need outstanding computational resources to solve increasingly complex problems, requiring the effective exploitation of cutting-edge HPC architectures. This paper focuses on the enabling of ORB5, a state-of-the-art, first-principles-based gyrokinetic code, on modern parallel hybrid multi-core, multi-GPU systems. ORB5 is a Lagrangian, Particle-In-Cell (PIC), fin…
▽ More
Gyrokinetic codes in plasma physics need outstanding computational resources to solve increasingly complex problems, requiring the effective exploitation of cutting-edge HPC architectures. This paper focuses on the enabling of ORB5, a state-of-the-art, first-principles-based gyrokinetic code, on modern parallel hybrid multi-core, multi-GPU systems. ORB5 is a Lagrangian, Particle-In-Cell (PIC), finite element, global, electromagnetic code, originally implementing distributed parallelism through MPI-based on domain decomposition and domain cloning. In order to support multi/many cores devices, the code has been completely refactored. Data structures have been re-designed to ensure efficient memory access, enhancing data locality. Multi-threading has been introduced through OpenMP on the CPU and adopting OpenACC to support GPU acceleration. MPI can still be used in combination with the two approaches. The performance results obtained using the full production ORB5 code on the Summit system at ORNL, on Piz Daint at CSCS and on the Marconi system at CINECA are presented, showing the effectiveness and performance portability of the adopted solutions: the same source code version was used to produce all results on all architectures.
△ Less
Submitted 14 February, 2020; v1 submitted 6 August, 2019;
originally announced August 2019.
-
Search for Light Dark Matter Interactions Enhanced by the Migdal effect or Bremsstrahlung in XENON1T
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
E. Angelino,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
D. Coderre
, et al. (109 additional authors not shown)
Abstract:
Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, whic…
▽ More
Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.
△ Less
Submitted 18 August, 2020; v1 submitted 30 July, 2019;
originally announced July 2019.
-
XENON1T Dark Matter Data Analysis: Signal Reconstruction, Calibration and Event Selection
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
F. Arneodo,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
D. Coderre,
A. P. Colijn,
J. Conrad
, et al. (103 additional authors not shown)
Abstract:
The XENON1T experiment at the Laboratori Nazionali del Gran Sasso is the most sensitive direct detection experiment for dark matter in the form of weakly interacting particles (WIMPs) with masses above $6\,$GeV/$c^2$ scattering off nuclei. The detector employs a dual-phase time projection chamber with 2.0 metric tons of liquid xenon in the target. A one metric $\mathrm{ton}\times\mathrm{year}$ exp…
▽ More
The XENON1T experiment at the Laboratori Nazionali del Gran Sasso is the most sensitive direct detection experiment for dark matter in the form of weakly interacting particles (WIMPs) with masses above $6\,$GeV/$c^2$ scattering off nuclei. The detector employs a dual-phase time projection chamber with 2.0 metric tons of liquid xenon in the target. A one metric $\mathrm{ton}\times\mathrm{year}$ exposure of science data was collected between October 2016 and February 2018. This article reports on the performance of the detector during this period and describes details of the data analysis that led to the most stringent exclusion limits on various WIMP-nucleon interaction models to date. In particular, signal reconstruction, event selection and calibration of the detector response to nuclear and electronic recoils in XENON1T are discussed.
△ Less
Submitted 4 November, 2019; v1 submitted 11 June, 2019;
originally announced June 2019.
-
The XENON1T Data Acquisition System
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
F. Arneodo,
D. Barge,
L. Baudis,
B. Bauermeister,
L. Bellagamba,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
L. Bütikofer,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
D. Coderre
, et al. (108 additional authors not shown)
Abstract:
The XENON1T liquid xenon time projection chamber is the most sensitive detector built to date for the measurement of direct interactions of weakly interacting massive particles with normal matter. The data acquisition system (DAQ) is constructed from commercial, open source, and custom components to digitize signals from the detector and store them for later analysis. The system achieves an extrem…
▽ More
The XENON1T liquid xenon time projection chamber is the most sensitive detector built to date for the measurement of direct interactions of weakly interacting massive particles with normal matter. The data acquisition system (DAQ) is constructed from commercial, open source, and custom components to digitize signals from the detector and store them for later analysis. The system achieves an extremely low signal threshold below a tenth of a photoelectron using a parallelized readout with the global trigger deferred to a later, software stage. The event identification is based on MongoDB database queries and has over 97% efficiency at recognizing interactions at the analysis energy threshold. A readout bandwidth over 300 MB/s is reached in calibration modes and is further expandable via parallelization. This DAQ system was successfully used during three years of operation of XENON1T.
△ Less
Submitted 7 August, 2019; v1 submitted 3 June, 2019;
originally announced June 2019.
-
A High Count-Rate and Depth-of-Interaction Resolving Single Layered One-Side Readout Pixelated Scintillator Crystal Array for PET Applications
Authors:
J. M. C. Brown,
S. E. Brunner,
D. R. Schaart
Abstract:
Organ-specific, targeted Field-of-View (FoV) Positron Emission Tomography (PET)/Magnetic Resonance Imaging (MRI) inserts are viable solutions for a number of imaging tasks where whole-body PET/MRI systems lack the necessary sensitivity and resolution. To meet the required PET detector performance of these systems, high count-rates and effective spatial resolutions on the order of a few mm, a novel…
▽ More
Organ-specific, targeted Field-of-View (FoV) Positron Emission Tomography (PET)/Magnetic Resonance Imaging (MRI) inserts are viable solutions for a number of imaging tasks where whole-body PET/MRI systems lack the necessary sensitivity and resolution. To meet the required PET detector performance of these systems, high count-rates and effective spatial resolutions on the order of a few mm, a novel two-axis patterned reflector foil pixelated scintillator crystal array design is developed and its proof-of-concept illustrated in-silico with the Monte Carlo radiation transport modelling toolkit Geant4. It is shown that the crystal surface roughness and phased open reflector cross-section patterns could be optimised to maximise either the PET radiation detector's effective spatial resolution, or count rate before event pile up. In addition, it was illustrated that these two parameters had minimal impact on the energy and time resolution of the proposed PET radiation detector design. Finally, it is shown that a PET radiation detector with balance performance could be constructed using ground crystals and phased open reflector cross-section pattern corresponding to the middle of the tested range.
△ Less
Submitted 18 September, 2019; v1 submitted 27 May, 2019;
originally announced May 2019.
-
ORB5: a global electromagnetic gyrokinetic code using the PIC approach in toroidal geometry
Authors:
E. Lanti,
N. Ohana,
N. Tronko,
T. Hayward-Schneider,
A. Bottino,
B. F. McMillan,
A. Mishchenko,
A. Scheinberg,
A. Biancalani,
P. Angelino,
S. Brunner,
J. Dominski,
P. Donnel,
C. Gheller,
R. Hatzky,
A. Jocksch,
S. Jolliet,
Z. X. Lu,
J. P. Martin Collar,
I. Novikau,
E. Sonnendrücker,
T. Vernay,
L. Villard
Abstract:
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization d…
▽ More
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid'' model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.
△ Less
Submitted 6 May, 2019;
originally announced May 2019.
-
First principles gyrokinetic analysis of electromagnetic plasma instabilities
Authors:
Natalia Tronko,
Alberto Bottino,
Cristel Chandre,
Eric Sonnendrücker,
Emmanuel Lanti,
Noé Ohana,
Stephan Brunner,
Laurent Villard
Abstract:
A two-fold analysis of electromagnetic core tokamak instabilities in the framework of the gyrokinetic theory is presented. First principle theoretical foundations of the gyrokinetic theory are used to explain and justify the numerical results obtained with the global electromagnetic particle-in-cell code Orb5 whose model is derived from the Lagrangian formalism. The energy conservation law corresp…
▽ More
A two-fold analysis of electromagnetic core tokamak instabilities in the framework of the gyrokinetic theory is presented. First principle theoretical foundations of the gyrokinetic theory are used to explain and justify the numerical results obtained with the global electromagnetic particle-in-cell code Orb5 whose model is derived from the Lagrangian formalism. The energy conservation law corresponding to the Orb5 model is derived from the Noether theorem and implemented in the code as a diagnostics for energy balance and conservation verification. An additional Noether theorem based diagnostics is implemented in order to analyse destabilising mechanisms for the electrostatic and the electromagnetic Ion Temperature Gradient (ITG) instabilities in the core region of the tokamak. The transition towards the Kinetic Ballooning Modes (KBM) at high electromagnetic $β$ is also investigated.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
XENON1T Dark Matter Data Analysis: Signal & Background Models, and Statistical Inference
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
V. C. Antochi,
F. Arneodo,
L. Baudis,
B. Bauermeister,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
D. Coderre,
A. P. Colijn,
J. Conrad,
J. P. Cussonneau
, et al. (101 additional authors not shown)
Abstract:
The XENON1T experiment searches for dark matter particles through their scattering off xenon atoms in a 2 tonne liquid xenon target. The detector is a dual-phase time projection chamber, which measures simultaneously the scintillation and ionization signals produced by interactions in target volume, to reconstruct energy and position, as well as the type of the interaction. The background rate in…
▽ More
The XENON1T experiment searches for dark matter particles through their scattering off xenon atoms in a 2 tonne liquid xenon target. The detector is a dual-phase time projection chamber, which measures simultaneously the scintillation and ionization signals produced by interactions in target volume, to reconstruct energy and position, as well as the type of the interaction. The background rate in the central volume of XENON1T detector is the lowest achieved so far with a liquid xenon-based direct detection experiment. In this work we describe the response model of the detector, the background and signal models, and the statistical inference procedures used in the dark matter searches with a 1 tonne$\times$year exposure of XENON1T data, that leaded to the best limit to date on WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c$^2$.
△ Less
Submitted 28 February, 2019; v1 submitted 28 February, 2019;
originally announced February 2019.
-
The effect of background flow shear on gyrokinetic turbulence in the cold ion limit
Authors:
Justin Ball,
Stephan Brunner,
Ben F. McMillan
Abstract:
The cold ion limit of the local gyrokinetic model is rigorously taken to produce a nonlinear system of fluid equations that includes background flow shear. No fluid closure is required. By considering a simple slab geometry with magnetic drifts, but no magnetic shear, these fluid equations reduce to the Charney-Hasegawa-Mima model in the presence of flow shear. Analytic solutions to this model are…
▽ More
The cold ion limit of the local gyrokinetic model is rigorously taken to produce a nonlinear system of fluid equations that includes background flow shear. No fluid closure is required. By considering a simple slab geometry with magnetic drifts, but no magnetic shear, these fluid equations reduce to the Charney-Hasegawa-Mima model in the presence of flow shear. Analytic solutions to this model are found to study the impact of ExB flow shear on the stability of a single Parallel Velocity Gradient (PVG) driven mode. Additionally, the model is used to investigate the effect of background ExB flow shear on the basic three-mode nonlinear coupling, which reveals differences between zonal and non-zonal modes. These analytic results agree with gyrokinetic simulations and can serve to benchmark the numerical implementation of flow shear and nonlinear coupling.
△ Less
Submitted 26 March, 2019; v1 submitted 20 December, 2018;
originally announced December 2018.
-
Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality
Authors:
R. Jorge,
P. Ricci,
S. Brunner,
S. Gamba,
V. Konovets,
N. F. Loureiro,
L. M. Perrone,
N. Teixeira
Abstract:
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength…
▽ More
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard-Bernstein and Dougherty collision operators, and the electron-ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode.
△ Less
Submitted 5 March, 2019; v1 submitted 26 November, 2018;
originally announced November 2018.
-
Pullback scheme implementation in ORB5
Authors:
Alexey Mishchenko,
Alberto Bottino,
Alessandro Biancalani,
Roman Hatzky,
Thomas Hayward-Schneider,
Noe Ohana,
Emmanuel Lanti,
Stephan Brunner,
Laurent Villard,
Matthias Borchardt,
Ralf Kleiber,
Axel Koenies
Abstract:
The pullback scheme is implemented in the global gyrokinetic particle-in-cell code ORB5 [S. Jolliet et al, Comp. Phys. Comm., 177, 409 (2007)] to mitigate the cancellation problem in electromagnetic simulations. The equations and the discretisation used by the code are described. Numerical simulations of the Toroidal Alfven Eigenmodes are performed in linear and nonlinear regimes to verify the sch…
▽ More
The pullback scheme is implemented in the global gyrokinetic particle-in-cell code ORB5 [S. Jolliet et al, Comp. Phys. Comm., 177, 409 (2007)] to mitigate the cancellation problem in electromagnetic simulations. The equations and the discretisation used by the code are described. Numerical simulations of the Toroidal Alfven Eigenmodes are performed in linear and nonlinear regimes to verify the scheme. A considerable improvement in the code efficiency is observed. For the internal kink mode, it is shown that the pullback mitigation efficiently cures a numerical instability which would make the simulation more costly otherwise.
△ Less
Submitted 13 November, 2018;
originally announced November 2018.
-
Dark Matter Search Results from a One Tonne$\times$Year Exposure of XENON1T
Authors:
E. Aprile,
J. Aalbers,
F. Agostini,
M. Alfonsi,
L. Althueser,
F. D. Amaro,
M. Anthony,
F. Arneodo,
L. Baudis,
B. Bauermeister,
M. L. Benabderrahmane,
T. Berger,
P. A. Breur,
A. Brown,
A. Brown,
E. Brown,
S. Bruenner,
G. Bruno,
R. Budnik,
C. Capelli,
J. M. R. Cardoso,
D. Cichon,
D. Coderre,
A. P. Colijn,
J. Conrad
, et al. (95 additional authors not shown)
Abstract:
We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra…
▽ More
We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.
△ Less
Submitted 13 September, 2018; v1 submitted 31 May, 2018;
originally announced May 2018.
-
Simulating background shear flow in local gyrokinetic simulations
Authors:
B. F. McMillan,
J. Ball,
S. Brunner
Abstract:
Local gyrokinetic simulations solve the gyrokinetic equations with homogeneous background gradients, typically using a doubly periodic domain in the (x,y) plane (i.e. perpendicular to the field line). Spatial Fourier representations are almost universal in local gyrokinetic codes, and the wavevector-remap method was introduced in [Hammett et. al., Bull Am Phys Soc VP1 136, (2006)] as a simple meth…
▽ More
Local gyrokinetic simulations solve the gyrokinetic equations with homogeneous background gradients, typically using a doubly periodic domain in the (x,y) plane (i.e. perpendicular to the field line). Spatial Fourier representations are almost universal in local gyrokinetic codes, and the wavevector-remap method was introduced in [Hammett et. al., Bull Am Phys Soc VP1 136, (2006)] as a simple method for expressing the local gyrokinetic equations with a background shear flow in a Fourier representation. Although extensively applied, the wavevector-remap method has not been formally shown to converge, and suffers from known unphysicality when the solutions are plotted in real space [Fox et. al. PPCF 59, 044008]. In this work, we use an analytic solution in slab geometry to demonstrate that wavevector-remap leads to incorrect smeared non-linear coupling between modes. We derive a correct, relatively simple method for solving local gyrokinetics in Fourier space with a background shear flow, and compare this to the wavevector-remap method. This allows us to show that the error in wavevector-remap can be seen as an incorrect rounding in wavenumber space in the nonlinear term. By making minor modifications to the nonlinear term, we implement the corrected wavevector-remap scheme in the GENE[25] code and compare results of the original and corrected wavevector-remap for standard nonlinear benchmark cases. Certain physical phenomena are impacted by the errors in the original remap scheme, and these numerical artefacts do not reduce as system size increases: that is, original wavevector-remap scheme does not converge to the correct result.
△ Less
Submitted 15 January, 2019; v1 submitted 10 November, 2017;
originally announced November 2017.