-
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
Authors:
The Astropy Collaboration,
Adrian M. Price-Whelan,
Pey Lian Lim,
Nicholas Earl,
Nathaniel Starkman,
Larry Bradley,
David L. Shupe,
Aarya A. Patil,
Lia Corrales,
C. E. Brasseur,
Maximilian Nöthe,
Axel Donath,
Erik Tollerud,
Brett M. Morris,
Adam Ginsburg,
Eero Vaher,
Benjamin A. Weaver,
James Tocknell,
William Jamieson,
Marten H. van Kerkwijk,
Thomas P. Robitaille,
Bruce Merry,
Matteo Bachetti,
H. Moritz Günther,
Thomas L. Aldcroft
, et al. (111 additional authors not shown)
Abstract:
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as…
▽ More
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
First direct detection constraints on Planck-scale mass dark matter with multiple-scatter signatures using the DEAP-3600 detector
Authors:
P. Adhikari,
R. Ajaj,
M. Alpízar-Venegas,
D. J. Auty,
H. Benmansour,
C. E. Bina,
W. Bonivento,
M. G. Boulay,
M. Cadeddu,
B. Cai,
M. Cárdenas-Montes,
S. Cavuoti,
Y. Chen,
B. T. Cleveland,
J. M. Corning,
S. Daugherty,
P. DelGobbo,
P. Di Stefano,
L. Doria,
M. Dunford,
E. Ellingwood,
A. Erlandson,
S. S. Farahani,
N. Fatemighomi,
G. Fiorillo
, et al. (72 additional authors not shown)
Abstract:
Dark matter particles with Planck-scale mass ($\simeq10^{19}\text{GeV}/c^2$) arise in well-motivated theories and could be produced by several cosmological mechanisms. Using a blind analysis of data collected over a 813 d live time with DEAP-3600, a 3.3 t single-phase liquid argon-based dark matter experiment at SNOLAB, a search for supermassive dark matter was performed, looking for multiple-scat…
▽ More
Dark matter particles with Planck-scale mass ($\simeq10^{19}\text{GeV}/c^2$) arise in well-motivated theories and could be produced by several cosmological mechanisms. Using a blind analysis of data collected over a 813 d live time with DEAP-3600, a 3.3 t single-phase liquid argon-based dark matter experiment at SNOLAB, a search for supermassive dark matter was performed, looking for multiple-scatter signals. No candidate signal events were observed, leading to the first direct detection constraints on Planck-scale mass dark matter. Leading limits constrain dark matter masses between $8.3\times10^{6}$ and $1.2\times10^{19} \text{GeV}/c^2$, and cross sections for scattering on $^{40}$Ar between $1.0\times10^{-23}$ and $2.4\times10^{-18} \text{cm}^2$. These are used to constrain two composite dark matter models.
△ Less
Submitted 5 January, 2022; v1 submitted 20 August, 2021;
originally announced August 2021.
-
Pulseshape discrimination against low-energy Ar-39 beta decays in liquid argon with 4.5 tonne-years of DEAP-3600 data
Authors:
The DEAP Collaboration,
P. Adhikari,
R. Ajaj,
M. Alpízar-Venegas,
P. -A. Amaudruz,
D. J. Auty,
M. Batygov,
B. Beltran,
H. Benmansour,
C. E. Bina,
J. Bonatt,
W. Bonivento,
M. G. Boulay,
B. Broerman,
J. F. Bueno,
P. M. Burghardt,
A. Butcher,
M. Cadeddu,
B. Cai,
M. Cárdenas-Montes,
S. Cavuoti,
M. Chen,
Y. Chen,
B. T. Cleveland,
J. M. Corning
, et al. (104 additional authors not shown)
Abstract:
The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from $^{39}$Ar beta decays and is suppressed using pulseshape discrimination (PSD).
We use two types of PSD algorithm: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window ar…
▽ More
The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from $^{39}$Ar beta decays and is suppressed using pulseshape discrimination (PSD).
We use two types of PSD algorithm: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the charge of a single photoelectron, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulseshape and for afterpulsing in the light detectors.
The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected.
△ Less
Submitted 6 April, 2021; v1 submitted 22 March, 2021;
originally announced March 2021.
-
Constraints on dark matter-nucleon effective couplings in the presence of kinematically distinct halo substructures using the DEAP-3600 detector
Authors:
P. Adhikari,
R. Ajaj,
C. E. Bina,
W. Bonivento,
M. G. Boulay,
M. Cadeddu,
B. Cai,
M. Cárdenas-Montes,
S. Cavuoti,
Y. Chen,
B. T. Cleveland,
J. M. Corning,
S. Daugherty,
P. DelGobbo,
P. Di Stefano,
L. Doria,
M. Dunford,
A. Erlandson,
S. S. Farahani,
N. Fatemighomi,
G. Fiorillo,
D. Gallacher,
E. A. Garcés,
P. García Abia,
S. Garg
, et al. (59 additional authors not shown)
Abstract:
DEAP-3600 is a single-phase liquid argon detector aiming to directly detect Weakly Interacting Massive Particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon spin-independent, isoscalar cross section. This study reinterprets this result within a Non-Relativistic Effect…
▽ More
DEAP-3600 is a single-phase liquid argon detector aiming to directly detect Weakly Interacting Massive Particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon spin-independent, isoscalar cross section. This study reinterprets this result within a Non-Relativistic Effective Field Theory framework, and further examines how various possible substructures in the local dark matter halo may affect these constraints. Such substructures are hinted at by kinematic structures in the local stellar distribution observed by the Gaia satellite and other recent astronomical surveys. These include the Gaia Sausage (or Enceladus), as well as a number of distinct streams identified in recent studies. Limits are presented for the coupling strength of the effective contact interaction operators $\mathcal{O}_1$, $\mathcal{O}_3$, $\mathcal{O}_5$, $\mathcal{O}_8$, and $\mathcal{O}_{11}$, considering isoscalar, isovector, and xenonphobic scenarios, as well as the specific operators corresponding to millicharge, magnetic dipole, electric dipole, and anapole interactions. The effects of halo substructures on each of these operators are explored as well, showing that the $\mathcal{O}_5$ and $\mathcal{O}_8$ operators are particularly sensitive to the velocity distribution, even at dark matter masses above 100 GeV/$c^2$.
△ Less
Submitted 5 January, 2022; v1 submitted 29 May, 2020;
originally announced May 2020.
-
EinsteinPy: A Community Python Package for General Relativity
Authors:
Shreyas Bapat,
Ritwik Saha,
Bhavya Bhatt,
Shilpi Jain,
Akshita Jain,
Sofía Ortín Vela,
Priyanshu Khandelwal,
Jyotirmaya Shivottam,
Jialin Ma,
Gim Seng Ng,
Pratyush Kerhalkar,
Hrishikesh Sudam Sarode,
Rishi Sharma,
Manvi Gupta,
Divya Gupta,
Tushar Tyagi,
Tanmay Rustagi,
Varun Singh,
Saurabh Bansal,
Naman Tayal,
Abhijeet Manhas,
Raphael Reyna,
Gaurav Kumar,
Govind Dixit,
Ratin Kumar
, et al. (24 additional authors not shown)
Abstract:
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increas…
▽ More
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.
△ Less
Submitted 22 May, 2020;
originally announced May 2020.
-
Waldmeier Effect in Stellar Cycles
Authors:
Suyog Garg,
Bidya Binay Karak,
Ricky Egeland,
Willie Soon,
Sallie Baliunas
Abstract:
One of the most robust features of the solar magnetic cycle is that the stronger cycles rise faster than the weaker ones. This is popularly known as the Waldmeier Effect, which is known for more than 80 years. This fundamental feature of the solar cycle has not only practical implications, e,g., in predicting the solar cycle, but also implications in understanding the solar dynamo. Here we ask the…
▽ More
One of the most robust features of the solar magnetic cycle is that the stronger cycles rise faster than the weaker ones. This is popularly known as the Waldmeier Effect, which is known for more than 80 years. This fundamental feature of the solar cycle has not only practical implications, e,g., in predicting the solar cycle, but also implications in understanding the solar dynamo. Here we ask the question whether the Waldmeier Effect exists in other Sun-like stars. To answer this question, we analyze the Ca \II{} H \& K S-index from Mount Wilson Observatory for 21 Sun-like G--K stars. We specifically check two aspects of Waldmeier Effect, namely, WE1: the anti-correlation between the rise times and the peaks and WE2: the positive correlation between rise rates and amplitudes. We show that except HD~16160, HD~81809, HD~155886 and HD~161239, all stars considered in the analysis show WE2. While WE1 is found to be present only in some of the stars studied. Further, the WE1 correlation is weaker than the WE2. Both WE1 and WE2 exist in the solar S-index as well. Similar to the solar cycles, the magnetic cycles of many stars are asymmetric about their maxima. The existence of the Waldmeier Effect and asymmetric cycles in Sun-like stars suggests that the dynamo mechanism which operates in the Sun is also operating in other stars.
△ Less
Submitted 26 September, 2019;
originally announced September 2019.
-
Electromagnetic Backgrounds and Potassium-42 Activity in the DEAP-3600 Dark Matter Detector
Authors:
R. Ajaj,
G. R. Araujo,
M. Batygov,
B. Beltran,
C. E. Bina,
M. G. Boulay,
B. Broerman,
J. F. Bueno,
P. M. Burghardt,
A. Butcher,
M. Cárdenas-Montes,
S. Cavuoti,
M. Chen,
Y. Chen,
B. T. Cleveland,
K. Dering,
F. A. Duncan,
M. Dunford,
A. Erlandson,
N. Fatemighomi,
G. Fiorillo,
A. Flower,
R. J. Ford,
D. Gallacher,
P. García Abia
, et al. (66 additional authors not shown)
Abstract:
The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $β$ and…
▽ More
The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $β$ and $γ$-rays originating from internal and external radioactivity in the detector material. A background model of the ER interactions in DEAP-3600 was developed and is described in this work. The model is based on several components which are expected from radioisotopes in the LAr, from ex-situ material assay measurements, and from dedicated independent in-situ analyses. This prior information is used in a Bayesian fit of the ER components to a 247.2 d dataset to model the radioactivity in the surrounding detector materials. While excellent discrimination between ERs and NRs is reached with pulse shape discrimination, utilizing the large difference between fast and slow components of LAr scintillation light, detailed knowledge of the ER background and activity of detector components, sets valuable constraints on other key types of backgrounds in the detector: neutrons and alphas. In addition, the activity of $^{42}$Ar in LAr in DEAP-3600 is determined by measuring the daughter decay of $^{42}$K. This cosmogenically activated trace isotope is a relevant background at higher energies for other rare event searches using atmospheric argon e.g. DarkSide-20k, GERDA or LEGEND. The specific activity of $^{42}$Ar in the atmosphere is found to be $40.4 \pm 5.9$ $μ$Bq/kg of argon.
△ Less
Submitted 14 May, 2019;
originally announced May 2019.
-
Search for dark matter with a 231-day exposure of liquid argon using DEAP-3600 at SNOLAB
Authors:
R. Ajaj,
P. -A. Amaudruz,
G. R. Araujo,
M. Baldwin,
M. Batygov,
B. Beltran,
C. E. Bina,
J. Bonatt,
M. G. Boulay,
B. Broerman,
J. F. Bueno,
P. M. Burghardt,
A. Butcher,
B. Cai,
S. Cavuoti,
M. Chen,
Y. Chen,
B. T. Cleveland,
D. Cranshaw,
K. Dering,
J. DiGioseffo,
L. Doria,
F. A. Duncan,
M. Dunford,
A. Erlandson
, et al. (77 additional authors not shown)
Abstract:
DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment, operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3279 kg of LAr contained in a spherical acrylic vessel. This paper reports on the analysis of a 758 tonne\cdot day exposure taken over a period of 231 live-days during the first year of operation. No candidate signal events are obs…
▽ More
DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment, operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3279 kg of LAr contained in a spherical acrylic vessel. This paper reports on the analysis of a 758 tonne\cdot day exposure taken over a period of 231 live-days during the first year of operation. No candidate signal events are observed in the WIMP-search region of interest, which results in the leading limit on the WIMP-nucleon spin-independent cross section on a LAr target of $3.9\times10^{-45}$ cm$^{2}$ ($1.5\times10^{-44}$ cm$^{2}$) for a 100 GeV/c$^{2}$ (1 TeV/c$^{2}$) WIMP mass at 90\% C. L. In addition to a detailed background model, this analysis demonstrates the best pulse-shape discrimination in LAr at threshold, employs a Bayesian photoelectron-counting technique to improve the energy resolution and discrimination efficiency, and utilizes two position reconstruction algorithms based on PMT charge and photon arrival times.
△ Less
Submitted 4 July, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
Bounds on Slow Roll and the de Sitter Swampland
Authors:
Sumit K. Garg,
Chethan Krishnan
Abstract:
The recently introduced swampland criterion for de Sitter (arXiv:1806.08362) can be viewed as a (hierarchically large) bound on the smallness of the slow roll parameter $ε_V$. This leads us to consider the other slow roll parameter $η_V$ more closely, and we are lead to conjecture that the bound is not necessarily on $ε_V$, but on slow roll itself. A natural refinement of the de Sitter swampland c…
▽ More
The recently introduced swampland criterion for de Sitter (arXiv:1806.08362) can be viewed as a (hierarchically large) bound on the smallness of the slow roll parameter $ε_V$. This leads us to consider the other slow roll parameter $η_V$ more closely, and we are lead to conjecture that the bound is not necessarily on $ε_V$, but on slow roll itself. A natural refinement of the de Sitter swampland conjecture is therefore that slow roll is violated at ${\cal O}(1)$ in Planck units in any UV complete theory. A corollary is that $ε_V$ need not necessarily be ${\cal O}(1)$, if $η_V \lesssim -{\cal O}(1)$ holds. We consider various tachyonic tree level constructions of de Sitter in IIA/IIB string theory (as well as closely related models of inflation), which superficially violate arXiv:1806.08362, and show that they are consistent with this refined version of the bound. The phrasing in terms of slow roll makes it plausible why both versions of the conjecture run into trouble when the number of e-folds during inflation is high. We speculate that one way to evade the bound could be to have a large number of fields, like in $N$-flation.
△ Less
Submitted 26 July, 2018; v1 submitted 13 July, 2018;
originally announced July 2018.
-
The New Minimal Supersymmetric GUT : Spectra, RG analysis and Fermion Fits
Authors:
Charanjit S. Aulakh,
Sumit K. Garg
Abstract:
The supersymmetric SO(10) GUT based on the ${\bf{210\oplus 10\oplus 120\oplus 126\oplus {\bar {126}}}}$ Higgs system provides a minimal framework for the emergence of the R-parity exact MSSM at low energies and a viable supersymmetric seesaw explanation for the observed neutrino masses and mixing angles. We present formulae for MSSM decomposition of the superpotential invariants, tree level light…
▽ More
The supersymmetric SO(10) GUT based on the ${\bf{210\oplus 10\oplus 120\oplus 126\oplus {\bar {126}}}}$ Higgs system provides a minimal framework for the emergence of the R-parity exact MSSM at low energies and a viable supersymmetric seesaw explanation for the observed neutrino masses and mixing angles. We present formulae for MSSM decomposition of the superpotential invariants, tree level light charged fermion effective Yukawa couplings, Weinberg neutrino mass generation operator, and the $d=5,ΔB=ΔL \neq 0$ effective superpotential in terms of GUT parameters. We use them to determine fits of the 18 available fermion mass-mixing data in terms of the superpotential parameters of the NMSGUT and SUGRY(NUHM) type soft supersymmetry breaking parameters ($\{m_{\tilde f},m_{1/2},A_0,M^2_{H,\bar H}\} $) specified at the MSSM one loop unification scale $M_X^0=10^{16.33} $ GeV. Our fits are compatible with electroweak symmetry breaking and Unification constraints and yield right-handed neutrino masses in the leptogenesis relevant range : $10^8-10^{13} $ GeV. Matching the SM data requires lowering the strange and down quark Yukawas in the MSSM via large $\tanβ$ driven threshold corrections and characteristic soft Susy breaking spectra. The Susy spectra have light pure Bino LSP, heavy exotic Higgs(inos) and large $ μ,A_0,M_{H,\bar{H}}$ parameters $\sim 100$ TeV. Typically third generation sfermions are much \emph{heavier} than the first two generations. The smuon is often the lightest charged sfermion thus offering a Bino-CDM co-annihilation channel. The parameter sets obtained are used to calculate B violation rates which are found to be generically much faster($\sim 10^{-28}\, yr^{-1}$) than the current experimental limits. Improvements which may allow acceptable B violation rates are identified.
△ Less
Submitted 3 January, 2012; v1 submitted 7 July, 2008;
originally announced July 2008.
