Showing 1–50 of 59 results for author: Mitra, P

Onset and evolution of solar flares: Application of 2D and 3D models of magnetic reconnection

Authors: Bhuwan Joshi, Prabir K. Mitra, Astrid M. Veronig, R. Bhattacharyya

Abstract: The contemporary multi-wavelength observations have revealed various important features during solar flares which, on one hand, support the two-dimensional (2D) "standard flare model" while, on other hand, also urge for the exploration of three-dimensional (3D) magnetic field topologies involved in flares. Traditionally, the formation of parallel ribbons on both side of the polarity inversion line… ▽ More The contemporary multi-wavelength observations have revealed various important features during solar flares which, on one hand, support the two-dimensional (2D) "standard flare model" while, on other hand, also urge for the exploration of three-dimensional (3D) magnetic field topologies involved in flares. Traditionally, the formation of parallel ribbons on both side of the polarity inversion line (PIL) and associated overlying loop arcades have been recognized as the most prominent features of eruptive flares which has formed the basis for the development of the standard model providing a 2D description of the flare-associated phenomena. The actual flare, however, occurs in a more complicated 3D magnetic structure. Thus, despite the general applicability, the standard model has limited or no scope in explaining some of the features which exclusively requires a 3D description. In this context, the observations of "circular ribbon flares" stand out where one of the ribbons presents an almost fully closed quasi-circular or quasi-ellipsoidal shape, evidencing the involvement of a typical fan-spine magnetic configuration. In this article, we discuss observational features vis à vis theoretical understanding of solar flares in view of 2D and 3D models of magnetic reconnection. We highlight a few complex cases of circular ribbon flares exhibiting parallel ribbons, a coronal jet, and/or an erupting magnetic flux rope. Exploring various 3D topologies also enables us to probe similarities between the circumstances that govern the onset of jets, confined flares and CME-producing eruptive flares. △ Less

Submitted 30 September, 2023; originally announced October 2023.

Comments: 21 pages, 5 figures, invited review in the 3rd BINA workshop, to be published in the Bulletin of the Liège Royal Society of Sciences

The Giant Radio Array for Neutrino Detection (GRAND) Collaboration -- Contributions to the 38th International Cosmic Ray Conference (ICRC 2023)

Authors: GRAND Collaboration, Rafael Alves Batista, Aurélien Benoit-Lévy, Teresa Bister, Mauricio Bustamante, Yiren Chen, LingMei Cheng, Simon Chiche, Jean-Marc Colley, Pablo Correa, Nicoleta Cucu Laurenciu, Zigao Dai, Beatriz de Errico, Sijbrand de Jong, João R. T. de Mello Neto, Krijn D. de Vries, Peter B. Denton, Valentin Deocoene, Kaikai Duan, Bohao Duan, Ralph Engel, Yizhong Fan, Arsène Ferrière, QuanBu Gou, Junhua Gu , et al. (74 additional authors not shown)

Abstract: The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the at… ▽ More The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground. In particular, for ultra-high-energy neutrinos, the future final phase of GRAND aims to be sensitive enough to discover them in spite of their plausibly tiny flux. Presently, three prototype GRAND radio arrays are in operation: GRANDProto300, in China, GRAND@Auger, in Argentina, and GRAND@Nancay, in France. Their goals are to field-test the design of the radio-detection units, understand the radio background to which they are exposed, and develop tools for diagnostic, data gathering, and data analysis. This list of contributions to the 38th International Cosmic Ray Conference (ICRC 2023) presents an overview of GRAND, in its present and future incarnations, and a look at the first data collected by GRANDProto13, the first phase of GRANDProto300. △ Less

Submitted 27 July, 2023; originally announced August 2023.

Comments: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

Offline Signal Identification with GRANDProto300

Authors: Pragati Mitra

Abstract: The GRANDProto300 (GP300) array is a pathfinder for the Giant Radio Array for Neutrino Detection (GRAND) project. Serving as a test bench, the GP300 array is expected to pioneer techniques of autonomous radio detection including identification and reconstruction of nearly horizontal cosmic-ray (CR) air showers, and shed light in understanding the interesting `transition region' from the galactic t… ▽ More The GRANDProto300 (GP300) array is a pathfinder for the Giant Radio Array for Neutrino Detection (GRAND) project. Serving as a test bench, the GP300 array is expected to pioneer techniques of autonomous radio detection including identification and reconstruction of nearly horizontal cosmic-ray (CR) air showers, and shed light in understanding the interesting `transition region' from the galactic to extragalactic CR sources. An offline analysis of signal identification over ambient noise is crucial at this stage, where very relaxed self-triggering thresholds of radio antennas will be used for study purposes. In this work, we show results and efficiency of signal identification with classical approaches using a wide set of simulated realistic signal templates and also validated against measured background recorded by deployed prototypes. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Search for Daily Modulation of MeV Dark Matter Signals with DAMIC-M

Authors: I. Arnquist, N. Avalos, D. Baxter, X. Bertou, N. Castello-Mor, A. E. Chavarria, J. Cuevas-Zepeda, A. Dastgheibi-Fard, C. De Dominicis, O. Deligny, J. Duarte-Campderros, E. Estrada, N. Gadola, R. Gaior, T. Hossbach, L. Iddir, B. J. Kavanagh, B. Kilminster, A. Lantero-Barreda, I. Lawson, S. Lee, A. Letessier-Selvon, P. Loaiza, A. Lopez-Virto, K. J. McGuire , et al. (15 additional authors not shown)

Abstract: Dark Matter (DM) particles with sufficiently large cross sections may scatter as they travel through Earth's bulk. The corresponding changes in the DM flux give rise to a characteristic daily modulation signal in detectors sensitive to DM-electron interactions. Here, we report results obtained from the first underground operation of the DAMIC-M prototype detector searching for such a signal from D… ▽ More Dark Matter (DM) particles with sufficiently large cross sections may scatter as they travel through Earth's bulk. The corresponding changes in the DM flux give rise to a characteristic daily modulation signal in detectors sensitive to DM-electron interactions. Here, we report results obtained from the first underground operation of the DAMIC-M prototype detector searching for such a signal from DM with MeV-scale mass. A model-independent analysis finds no modulation in the rate of 1$e^-$ events with periods in the range 1-48 h. We then use these data to place exclusion limits on DM in the mass range [0.53, 2.7] MeV/c$^2$ interacting with electrons via a dark photon mediator. Taking advantage of the time-dependent signal we improve by $\sim$2 orders of magnitude on our previous limit obtained from the total rate of 1$e^-$ events, using the same data set. This daily modulation search represents the current strongest limit on DM-electron scattering via ultralight mediators for DM masses around 1 MeV/c$^2$. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Reconstructing Air Shower Parameters with MGMR3D

Authors: P. Mitra, O. Scholten, T. N. G. Trinh, S. Buitink, J. Bhavani, A. Corstanje, M. Desmet, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, N. Karastathis, G. K. Krampah, K. Mulrey, A. Nelles, H. Pandya, S. Thoudam, K. D. de Vries, S. ter Veen

Abstract: Measuring the radio emission from cosmic ray particle cascades has proven to be a very efficient method to determine their properties such as the mass composition. Efficient modeling of the radio emission from air showers is crucial in order to extract the cosmic ray physics parameters from the measured radio emission. MGMR3D is a fast semi-analytic code that calculates the complete radio footprin… ▽ More Measuring the radio emission from cosmic ray particle cascades has proven to be a very efficient method to determine their properties such as the mass composition. Efficient modeling of the radio emission from air showers is crucial in order to extract the cosmic ray physics parameters from the measured radio emission. MGMR3D is a fast semi-analytic code that calculates the complete radio footprint, i.e.\ intensity, polarization, and pulse shapes, for a parametrized shower-current density and can be used in a chi-square optimization to fit a given radio data. It is many orders of magnitude faster than its Monte Carlo counterparts. We provide a detailed comparative study of MGMR3D to Monte Carlo simulations, where, with improved parametrizations, the shower maximum $\Xmax$ is found to have very strong agreement with a small dependency on the incoming zenith angle of the shower. Another interesting feature we observe with MGMR3D is sensitivity to the shape of the longitudinal profile in addition to $\Xmax$. This is achieved by probing the distinguishable radio footprint produced by a shower having a different longitudinal profile than usual. Furthermore, for the first time, we show the results of reconstructing shower parameters for LOFAR data using MGMR3D, and obtaining a $\Xmax$ resolution of 22 g/cm$^2$ and energy resolution of 19\%. △ Less

Submitted 9 July, 2023; originally announced July 2023.

Constraining the cosmic-ray mass composition by measuring the shower length with SKA

Authors: S. Buitink, A. Corstanje, J. Bhavani, M. Desmet, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, N. Karasthatis, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, K. Nivedita, H. Pandya, J. P. Rachen, O. Scholten, S. Thoudam, G. Trinh, S. ter Veen

Abstract: The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of… ▽ More The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Proceedings 9th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities - ARENA2022, 7-10 June 2022, Santiago de Compostela, Spain (8 pages)

A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

Authors: A. Corstanje, S. Buitink, M. Desmet, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, V. B. Jhansi, N. Karastathis, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, K. Nivedita, H. Pandya, O. Scholten, K. Terveer, S. Thoudam, G. Trinh, S. ter Veen

Abstract: Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly $60,000$ antennas envisioned to be used for air shower studies at the Square Kilometre Array… ▽ More Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly $60,000$ antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs. We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid. Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods. We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings. We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy. △ Less

Submitted 16 August, 2023; v1 submitted 23 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures. Accepted for publication in JINST (Journal of Instrumentation)

Journal ref: JINST 18 P09005, 2023

Prospects for measuring the longitudinal particle distribution of cosmic-ray air showers with SKA

Authors: A. Corstanje, S. Buitink, J. Bhavani, M. Desmet, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, N. Karasthatis, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, K. Nivedita, H. Pandya, J. P. Rachen, O. Scholten, S. Thoudam, G. Trinh, S. ter Veen

Abstract: We explore the possibilities of measuring the longitudinal profile of individual air showers beyond $X_{\rm max}$ when using very dense radio arrays such as SKA. The low-frequency part of the Square Kilometre Array, to be built in Australia, features an enormous antenna density of about $50,000$ antennas in the inner core region of radius 500 m, with a frequency band from 50 to 350 MHz. From CoREA… ▽ More We explore the possibilities of measuring the longitudinal profile of individual air showers beyond $X_{\rm max}$ when using very dense radio arrays such as SKA. The low-frequency part of the Square Kilometre Array, to be built in Australia, features an enormous antenna density of about $50,000$ antennas in the inner core region of radius 500 m, with a frequency band from 50 to 350 MHz. From CoREAS simulations, a SKA-Low antenna model plus noise contributions, and adapted LOFAR analysis scripts, we obtain a resolution in the shower maximum $X_{\rm max}$ and energy that is considerably better than at LOFAR. Already from this setup, we show that at least one additional parameter of the longitudinal profile can be measured. This would improve mass composition analysis by measuring an additional composition-dependent quantity. Moreover, it would offer an opportunity to discriminate between the different predictions of hadronic interaction models, hence contributing to hadronic physics at energy levels beyond man-made accelerators. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 8 pages, 6 figures. Presented at ARENA Conference 2022, Santiago de Compostela, Spain

First Constraints from DAMIC-M on Sub-GeV Dark-Matter Particles Interacting with Electrons

Authors: I. Arnquist, N. Avalos, D. Baxter, X. Bertou, N. Castello-Mor, A. E. Chavarria, J. Cuevas-Zepeda, J. Cortabitarte Gutierrez, J. Duarte-Campderros, A. Dastgheibi-Fard, O. Deligny, C. De Dominicis, E. Estrada, N. Gadola, R. Gaıor, T. Hossbach, L. Iddir, L. Khalil, B. Kilminster, A. Lantero-Barreda, I. Lawson, S. Lee, A. Letessier-Selvon, P. Loaiza, A. Lopez-Virto , et al. (20 additional authors not shown)

Abstract: We report constraints on sub-GeV dark matter particles interacting with electrons from the first underground operation of DAMIC-M detectors. The search is performed with an integrated exposure of 85.23 g days, and exploits the sub-electron charge resolution and low level of dark current of DAMIC-M Charge-Coupled Devices (CCDs). Dark-matter-induced ionization signals above the detector dark current… ▽ More We report constraints on sub-GeV dark matter particles interacting with electrons from the first underground operation of DAMIC-M detectors. The search is performed with an integrated exposure of 85.23 g days, and exploits the sub-electron charge resolution and low level of dark current of DAMIC-M Charge-Coupled Devices (CCDs). Dark-matter-induced ionization signals above the detector dark current are searched for in CCD pixels with charge up to 7 e-. With this data set we place limits on dark matter particles of mass between 0.53 and 1000 MeV/c2, excluding unexplored regions of parameter space in the mass ranges [1.6,1000] MeV/c2 and [1.5,15.1] MeV/c2 for ultra-light and heavy mediator interactions, respectively. △ Less

Submitted 5 February, 2023; originally announced February 2023.

Search for inelastic dark matter-nucleus scattering with the PICO-60 CF$_{3}$I and C$_{3}$F$_{8}$ bubble chambers

Authors: E. Adams, B. Ali, I. J. Arnquist, D. Baxter, E. Behnke, M. Bressler, B. Broerman, C. J. Chen, K. Clark, J. I. Collar, P. S. Cooper, C. Cripe, M. Crisler, C. E. Dahl, M. Das, S. Fallows, J. Farine, R. Filgas, A. García Viltres, G. Giroux, O. Harris, T. Hillier, E. W. Hoppe, C. M. Jackson, M. Jin , et al. (30 additional authors not shown)

Abstract: PICO bubble chambers have exceptional sensitivity to inelastic dark matter-nucleus interactions due to a combination of their extended nuclear recoil energy detection window from a few keV to $O$(100 keV) or more and the use of iodine as a heavy target. Inelastic dark matter-nucleus scattering is interesting for studying the properties of dark matter, where many theoretical scenarios have been dev… ▽ More PICO bubble chambers have exceptional sensitivity to inelastic dark matter-nucleus interactions due to a combination of their extended nuclear recoil energy detection window from a few keV to $O$(100 keV) or more and the use of iodine as a heavy target. Inelastic dark matter-nucleus scattering is interesting for studying the properties of dark matter, where many theoretical scenarios have been developed. This study reports the results of a search for dark matter inelastic scattering with the PICO-60 bubble chambers. The analysis reported here comprises physics runs from PICO-60 bubble chambers using CF$_{3}$I and C$_{3}$F$_{8}$. The CF$_{3}$I run consisted of 36.8 kg of CF$_{3}$I reaching an exposure of 3415 kg-day operating at thermodynamic thresholds between 7 and 20 keV. The C$_{3}$F$_{8}$ runs consisted of 52 kg of C$_{3}$F$_{8}$ reaching exposures of 1404 kg-day and 1167 kg-day running at thermodynamic thresholds of 2.45 keV and 3.29 keV, respectively. The analysis disfavors various scenarios, in a wide region of parameter space, that provide a feasible explanation of the signal observed by DAMA, assuming an inelastic interaction, considering that the PICO CF$_{3}$I bubble chamber used iodine as the target material. △ Less

Submitted 21 January, 2023; originally announced January 2023.

Comments: 7 pages, 3 figures

Origin of extreme solar eruptive activity from the active region NOAA 12673 and the largest flare of solar cycle 24

Authors: Bhuwan Joshi, Prabir K. Mitra

Abstract: During 2017, when the Sun was moving toward the minimum phase of solar cycle 24, an exceptionally eruptive active region (AR) NOAA 12673 emerged on the Sun during August 28-September 10. During the highest activity level, the AR turned into a delta-type sunspot region, which manifests the most complex configuration of magnetic fields from the photosphere to the coronal heights. The AR 12673 produc… ▽ More During 2017, when the Sun was moving toward the minimum phase of solar cycle 24, an exceptionally eruptive active region (AR) NOAA 12673 emerged on the Sun during August 28-September 10. During the highest activity level, the AR turned into a delta-type sunspot region, which manifests the most complex configuration of magnetic fields from the photosphere to the coronal heights. The AR 12673 produced four X-class and 27 M-class flares, along with numerous C-class flares, making it one of the most powerful ARs of solar cycle 24. Notably, it produced the largest flare of solar cycle 24, namely, the X9.3 event on 2017 September 6. In this work, we highlight the results of our comprehensive analysis involving multi-wavelength imaging and coronal magnetic field modeling to understand the evolution and eruptivity from AR 12673. We especially focus on the morphological, spectral and kinematical evolution of the two X-class flares on 6 September 2017. We explore various large- and small-scale magnetic field structures of the active region which are associated with the triggering and subsequent outbursts during the powerful solar transients. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: IAU S372 Proceedings Series (8 pages, 5 figures)

The DAMIC-M Experiment: Status and First Results

Authors: I. Arnquist, N. Avalos, P. Bailly, D. Baxter, X. Bertou, M. Bogdan, C. Bourgeois, J. Brandt, A. Cadiou, N. Castelló-Mor, A. E. Chavarria, M. Conde, N. J. Corso, J. Cortabitarte Gutiérrez, J. Cuevas-Zepeda, A. Dastgheibi-Fard, C. De Dominicis, O. Deligny, R. Desani, M. Dhellot, J-J. Dormard, J. Duarte-Campderros, E. Estrada, D. Florin, N. Gadola , et al. (47 additional authors not shown)

Abstract: The DAMIC-M (DArk Matter In CCDs at Modane) experiment employs thick, fully depleted silicon charged-coupled devices (CCDs) to search for dark matter particles with a target exposure of 1 kg-year. A novel skipper readout implemented in the CCDs provides single electron resolution through multiple non-destructive measurements of the individual pixel charge, pushing the detection threshold to the eV… ▽ More The DAMIC-M (DArk Matter In CCDs at Modane) experiment employs thick, fully depleted silicon charged-coupled devices (CCDs) to search for dark matter particles with a target exposure of 1 kg-year. A novel skipper readout implemented in the CCDs provides single electron resolution through multiple non-destructive measurements of the individual pixel charge, pushing the detection threshold to the eV-scale. DAMIC-M will advance by several orders of magnitude the exploration of the dark matter particle hypothesis, in particular of candidates pertaining to the so-called "hidden sector." A prototype, the Low Background Chamber (LBC), with 20g of low background Skipper CCDs, has been recently installed at Laboratoire Souterrain de Modane and is currently taking data. We will report the status of the DAMIC-M experiment and first results obtained with LBC commissioning data. △ Less

Submitted 25 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 10 pages, 6 figures, Submission to SciPost Physics Proceedings: 14th International Conference on Identification of Dark Matter (IDM) 2022

Precision measurement of Compton scattering in silicon with a skipper CCD for dark matter detection

Authors: D. Norcini, N. Castello-Mor, D. Baxter, N. J. Corso, J. Cuevas-Zepeda, C. De Dominicis, A. Matalon, S. Munagavalasa, S. Paul, P. Privitera, K. Ramanathan, R. Smida, R. Thomas, R. Yajur, A. E. Chavarria, K. McGuire, P. Mitra, A. Piers, M. Settimo, J. Cortabitarte Gutierrez, J. Duarte-Campderros, A. Lantero-Barreda, A. Lopez-Virto, I. Vila, R. Vilar , et al. (19 additional authors not shown)

Abstract: Experiments aiming to directly detect dark matter through particle recoils can achieve energy thresholds of $\mathcal{O}(1\,\mathrm{eV})$. In this regime, ionization signals from small-angle Compton scatters of environmental $γ$-rays constitute a significant background. Monte Carlo simulations used to build background models have not been experimentally validated at these low energies. We report a… ▽ More Experiments aiming to directly detect dark matter through particle recoils can achieve energy thresholds of $\mathcal{O}(1\,\mathrm{eV})$. In this regime, ionization signals from small-angle Compton scatters of environmental $γ$-rays constitute a significant background. Monte Carlo simulations used to build background models have not been experimentally validated at these low energies. We report a precision measurement of Compton scattering on silicon atomic shell electrons down to 23$\,$eV. A skipper charge-coupled device (CCD) with single-electron resolution, developed for the DAMIC-M experiment, was exposed to a $^{241}$Am $γ$-ray source over several months. Features associated with the silicon K, L$_{1}$, and L$_{2,3}$-shells are clearly identified, and scattering on valence electrons is detected for the first time below 100$\,$eV. We find that the relativistic impulse approximation for Compton scattering, which is implemented in Monte Carlo simulations commonly used by direct detection experiments, does not reproduce the measured spectrum below 0.5$\,$keV. The data are in better agreement with $ab$ $initio$ calculations originally developed for X-ray absorption spectroscopy. △ Less

Submitted 2 July, 2022; originally announced July 2022.

Comments: 12 pages, 10 figures

Report number: Phys. Rev. D 106, 092001

Journal ref: 3 November 2022

Determining the bubble nucleation efficiency of low-energy nuclear recoils in superheated C$_3$F$_8$ dark matter detectors

Authors: B. Ali, I. J. Arnquist, D. Baxter, E. Behnke, M. Bressler, B. Broerman, K. Clark, J. I. Collar, P. S. Cooper, C. Cripe, M. Crisler, C. E. Dahl, M. Das, D. Durnford, S. Fallows, J. Farine, R. Filgas, A. García-Viltres, F. Girard, G. Giroux, O. Harris, E. W. Hoppe, C. M. Jackson, M. Jin, C. B. Krauss , et al. (32 additional authors not shown)

Abstract: The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. The PICO Collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C$_3$F$_8$ from calibration data taken with 5 distinct… ▽ More The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. The PICO Collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C$_3$F$_8$ from calibration data taken with 5 distinct neutron spectra at various thermodynamic thresholds ranging from 2.1 keV to 3.9 keV. Instead of assuming any particular functional forms for the nuclear recoil efficiency, a generalized piecewise linear model is proposed with systematic errors included as nuisance parameters to minimize model-introduced uncertainties. A Markov-Chain Monte-Carlo (MCMC) routine is applied to sample the nuclear recoil efficiency for fluorine and carbon at 2.45 keV and 3.29 keV thermodynamic thresholds simultaneously. The nucleation efficiency for fluorine was found to be $\geq 50\, \%$ for nuclear recoils of 3.3 keV (3.7 keV) at a thermodynamic Seitz threshold of 2.45 keV (3.29 keV), and for carbon the efficiency was found to be $\geq 50\, \%$ for recoils of 10.6 keV (11.1 keV) at a threshold of 2.45 keV (3.29 keV). Simulated data sets are used to calculate a p-value for the fit, confirming that the model used is compatible with the data. The fit paradigm is also assessed for potential systematic biases, which although small, are corrected for. Additional steps are performed to calculate the expected interaction rates of WIMPs in the PICO-60 detector, a requirement for calculating WIMP exclusion limits. △ Less

Submitted 7 November, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: 17 pages, 22 figures, 5 tables

Results on photon-mediated dark matter-nucleus interactions from the PICO-60 C$_{3}$F$_{8}$ bubble chamber

Authors: B. Ali, I. J. Arnquist, D. Baxter, E. Behnke, M. Bressler, B. Broerman, C. J. Chen, K. Clark, J. I. Collar, P. S. Cooper, C. Cripe, M. Crisler, C. E. Dahl, M. Das, D. Durnford, S. Fallows, J. Farine, R. Filgas, A. García-Viltres, G. Giroux, O. Harris, T. Hillier, E. W. Hoppe, C. M. Jackson, M. Jin , et al. (30 additional authors not shown)

Abstract: Many compelling models predict dark matter coupling to the electromagnetic current through higher multipole interactions, while remaining electrically neutral. Different multipole couplings have been studied, among them anapole moment, electric and magnetic dipole moments, and millicharge. This study sets limits on the couplings for these photon-mediated interactions using non-relativistic contact… ▽ More Many compelling models predict dark matter coupling to the electromagnetic current through higher multipole interactions, while remaining electrically neutral. Different multipole couplings have been studied, among them anapole moment, electric and magnetic dipole moments, and millicharge. This study sets limits on the couplings for these photon-mediated interactions using non-relativistic contact operators in an effective field theory framework. Using data from the PICO-60 bubble chamber leading limits for dark matter masses between 2.7 GeV/c$^2$ and 24 GeV/c$^2$ are reported for the coupling of these photon-mediated dark matter-nucleus interactions. The detector was filled with 52 kg of C$_3$F$_8$ operating at thermodynamic thresholds of 2.45 keV and 3.29 keV, reaching exposures of 1404 kg-day and 1167 kg-day, respectively. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Homologous compact major blowout-eruption solar flares and their production of broad CMEs

Authors: Suraj Sahu, Bhuwan Joshi, Alphonse C. Sterling, Prabir K. Mitra, Ronald L. Moore

Abstract: We analyze the formation mechanism of three homologous broad coronal mass ejections (CMEs) resulting from a series of solar blowout-eruption flares with successively increasing intensities (M2.0, M2.6, and X1.0). The flares originated from active region NOAA 12017 during 2014 March 28-29 within an interval of approximately 24 hr. Coronal magnetic field modeling based on nonlinear-force-free-field… ▽ More We analyze the formation mechanism of three homologous broad coronal mass ejections (CMEs) resulting from a series of solar blowout-eruption flares with successively increasing intensities (M2.0, M2.6, and X1.0). The flares originated from active region NOAA 12017 during 2014 March 28-29 within an interval of approximately 24 hr. Coronal magnetic field modeling based on nonlinear-force-free-field extrapolation helps to identify low-lying closed bipolar loops within the flaring region enclosing magnetic flux ropes. We obtain a double flux rope system under closed bipolar fields for the all the events. The sequential eruption of the flux ropes led to homologous flares, each followed by a CME. Each of the three CMEs formed from the eruptions gradually attain a large angular width, after expanding from the compact eruption-source site. We find these eruptions and CMEs to be consistent with the 'magnetic-arch-blowout' scenario: each compact-flare blowout eruption was seated in one foot of a far-reaching magnetic arch, exploded up the encasing leg of the arch, and blew out the arch to make a broad CME. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: Accepted for publication in The Astrophysical Journal

Determining Atmospheric Electric Fields using MGMR3D

Authors: T. N. G. Trinh, O. Scholten, S. Buitink, K. D. de Vries, P. Mitra, T. Phong Nguyen, D. T. Si

Abstract: Cosmic-ray particles impinging on the atmosphere induce high-energy particle cascades in air, an Extensive Air Shower (EAS), emitting coherent radio emission. This emission is affected by the presence of strong electric fields during thunderstorm conditions. To reconstruct the atmospheric electric field from the measured radio footprint of the EAS we use an analytic model for the calculation of th… ▽ More Cosmic-ray particles impinging on the atmosphere induce high-energy particle cascades in air, an Extensive Air Shower (EAS), emitting coherent radio emission. This emission is affected by the presence of strong electric fields during thunderstorm conditions. To reconstruct the atmospheric electric field from the measured radio footprint of the EAS we use an analytic model for the calculation of the radio emission, MGMR3D. In this work we make an extensive comparison between the results of a microscopic model for radio emission, CoREAS, to obtain an improved parametrization for MGMR3D in the presence of atmospheric electric fields, as well as confidence intervals. The approach to extract the electric field structure is applied successfully to an event with a complicated radio footprint measured by LOFAR during thunderstorm conditions. This shows that, with the improved parametrization, MGMR3D can be used to extract the structure of the atmospheric electric field. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: The paper is accepted for publication in Physical Review D

EXCESS workshop: Descriptions of rising low-energy spectra

Authors: P. Adari, A. Aguilar-Arevalo, D. Amidei, G. Angloher, E. Armengaud, C. Augier, L. Balogh, S. Banik, D. Baxter, C. Beaufort, G. Beaulieu, V. Belov, Y. Ben Gal, G. Benato, A. Benoît, A. Bento, L. Bergé, A. Bertolini, R. Bhattacharyya, J. Billard, I. M. Bloch, A. Botti, R. Breier, G. Bres, J-. L. Bret , et al. (281 additional authors not shown)

Abstract: Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was… ▽ More Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was initiated. In its first iteration in June 2021, ten rare event search collaborations contributed to this initiative via talks and discussions. The contributing collaborations were CONNIE, CRESST, DAMIC, EDELWEISS, MINER, NEWS-G, NUCLEUS, RICOCHET, SENSEI and SuperCDMS. They presented data about their observed energy spectra and known backgrounds together with details about the respective measurements. In this paper, we summarize the presented information and give a comprehensive overview of the similarities and differences between the distinct measurements. The provided data is furthermore publicly available on the workshop's data repository together with a plotting tool for visualization. △ Less

Submitted 4 March, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 44 pages, 20 figures; Editors: A. Fuss, M. Kaznacheeva, F. Reindl, F. Wagner; updated copyright statements and funding information

Journal ref: SciPost Phys. Proc. 9, 001 (2022)

Multiwavelength Signatures of Episodic Null Point Reconnection in a Quadrupolar Magnetic Configuration and the Cause of Failed Flux Rope Eruption

Authors: Prabir K. Mitra, Bhuwan Joshi, Astrid M. Veronig, Thomas Wiegelmann

Abstract: In this paper, we present multiwavelength observations of the triggering of a failed-eruptive M-class flare from the active region NOAA 11302, and investigate the possible reasons for the associated failed eruption. Photospheric observations and Non-Linear Force Free Field extrapolated coronal magnetic field revealed that the flaring region had a complex quadrupolar configuration with a pre-existi… ▽ More In this paper, we present multiwavelength observations of the triggering of a failed-eruptive M-class flare from the active region NOAA 11302, and investigate the possible reasons for the associated failed eruption. Photospheric observations and Non-Linear Force Free Field extrapolated coronal magnetic field revealed that the flaring region had a complex quadrupolar configuration with a pre-existing coronal null point situated above the core field. Prior to the onset of the M-class flare, we observed multiple periods of small-scale flux enhancements in GOES and RHESSI soft X-ray observations, from the location of the null point. The pre-flare configuration and evolution reported here are similar to the ones presented in the breakout model but at much lower coronal heights. The core of the flaring region was characterized by the presence of two flux ropes in a double-decker configuration. During the impulsive phase of the flare, one of the two flux ropes initially started erupting but resulted in a failed eruption. Calculation of the magnetic decay index revealed a saddle-like profile where decay index initially increased to the torus unstable limits within the heights of the flux ropes but then decreased rapidly reaching to negative values, which was most likely responsible for the failed eruption of the initially torus unstable flux rope. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 10 figures, 1 table, accepted for publication in the Astrophysical Journal

A Transient Coronal Sigmoid in Active Region NOAA 11909: Build-up Phase, M-class Eruptive Flare, and Associated Fast Coronal Mass Ejection

Authors: Hema Kharayat, Bhuwan Joshi, Prabir K. Mitra, P. K. Manoharan, Christian Monstein

Abstract: In this article, we investigate the formation and disruption of a coronal sigmoid from the active region (AR) NOAA 11909 on 07 December 2013, by analyzing multi-wavelength and multi-instrument observations. Our analysis suggests that the formation of `transient' sigmoid initiated $\approx$1 hour before its eruption through a coupling between two twisted coronal loop systems. A comparison between c… ▽ More In this article, we investigate the formation and disruption of a coronal sigmoid from the active region (AR) NOAA 11909 on 07 December 2013, by analyzing multi-wavelength and multi-instrument observations. Our analysis suggests that the formation of `transient' sigmoid initiated $\approx$1 hour before its eruption through a coupling between two twisted coronal loop systems. A comparison between coronal and photospheric images suggests that the coronal sigmoid was formed over a simple $β$-type AR which also possessed dispersed magnetic field structure in the photosphere. The line-of-sight photospheric magnetograms also reveal moving magnetic features, small-scale flux cancellation events near the PIL, and overall flux cancellation during the extended pre-eruption phase which suggest the role of tether-cutting reconnection toward the build-up of the flux rope. The disruption of the sigmoid proceeded with a two-ribbon eruptive M1.2 flare (SOL2013-12-07T07:29). In radio frequencies, we observe type III and type II bursts in meter wavelengths during the impulsive phase of the flare. The successful eruption of the flux rope leads to a fast coronal mass ejection (with a linear speed of $\approx$1085 km s -1 ) in SOHO/LASCO field-of-view. During the evolution of the flare, we clearly observe typical "sigmoid-to-arcade" transformation. Prior to the onset of the impulsive phase of the flare, flux rope undergoes a slow rise ($\approx$15 km s -1 ) which subsequently transitions into a fast eruption ($\approx$110 km s -1 ). The two-phase evolution of the flux rope shows temporal associations with the soft X-ray precursor and impulsive phase emissions of the M-class flare, respectively, thus pointing toward a feedback relationship between magnetic reconnection and early CME dynamics. △ Less

Submitted 2 May, 2021; originally announced May 2021.

Comments: 28 pages, 14 figures, accepted in Solar Physics Journal

Two-stage evolution of an extended C-class eruptive flaring activity from sigmoid active region NOAA 12734: SDO and Udaipur-CALLISTO observations

Authors: Bhuwan Joshi, Prabir K. Mitra, R. Bhattacharyya, Kushagra Upadhyay, Divya Oberoi, K. Sasikumar Raja, Christian Monstein

Abstract: We present a multi-wavelength investigation of a C-class flaring activity that occurred in the active region NOAA 12734 on 8 March 2019. The investigation utilises data from AIA and HMI on board the SDO and the Udaipur-CALLISTO solar radio spectrograph of the Physical Research Laboratory. This low intensity C1.3 event is characterised by typical features of a long duration event (LDE), viz. extend… ▽ More We present a multi-wavelength investigation of a C-class flaring activity that occurred in the active region NOAA 12734 on 8 March 2019. The investigation utilises data from AIA and HMI on board the SDO and the Udaipur-CALLISTO solar radio spectrograph of the Physical Research Laboratory. This low intensity C1.3 event is characterised by typical features of a long duration event (LDE), viz. extended flare arcade, large-scale two-ribbon structures and twin coronal dimmings. The eruptive event occurred in a coronal sigmoid and displayed two distinct stages of energy release, manifested in terms of temporal and spatial evolution. The formation of twin dimming regions are consistent with the eruption of a large flux rope with footpoints lying in the western and eastern edges of the coronal sigmoid. The metric radio observations obtained from Udaipur-CALLISTO reveals a broad-band ($\approx$50-180 MHz), stationary plasma emission for $\approx$7 min during the second stage of the flaring activity that resemble a type IV radio burst. A type III decametre-hectometre radio bursts with starting frequency of $\approx$2.5 MHz precedes the stationary type IV burst observed by Udaipur-CALLISTO by $\approx$5 min. The synthesis of multi-wavelength observations and Non-Linear Force Free Field (NLFFF) coronal modelling together with magnetic decay index analysis suggests that the sigmoid flux rope underwent a zipping-like uprooting from its western to eastern footpoints in response to the overlying asymmetric magnetic field confinement. The asymmetrical eruption of the flux rope also accounts for the observed large-scale structures viz. apparent eastward shift of flare ribbons and post flare loops along the polarity inversion line (PIL), and provides an evidence for lateral progression of magnetic reconnection site as the eruption proceeds. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 30 pages, 14 figures, 1 table; accepted for publication in Solar Physics

Depth of shower maximum and mass composition of cosmic rays from 50 PeV to 2 EeV measured with the LOFAR radio telescope

Authors: A. Corstanje, S. Buitink, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, O. Scholten, S. ter Veen, S. Thoudam, G. Trinh, T. Winchen

Abstract: We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ o… ▽ More We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of $X_{\rm max}$ versus primary energy, and analyze the $X_{\rm max}$-dataset at distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of $X_{\rm max}$-distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit $23$ to $39$ $\%$ protons plus helium, depending on the choice of hadronic interaction model. The fraction of intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on $X_{\rm max}$ now lowered to 7 to $9$ $\mathrm{g/cm^2}$. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average $X_{\rm max}$. The values are in tension with those found at Pierre Auger Observatory, but agree with results from other cosmic ray observatories based in the Northern hemisphere. △ Less

Submitted 2 May, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 24 pages, 14 figures. Accepted for publication in Phys. Rev. D

Journal ref: Phys. Rev. D 103, 102006 (2021)

Evolutionary stages and triggering process of a complex eruptive flare with circular and parallel ribbons

Authors: Navin Chandra Joshi, Bhuwan Joshi, Prabir K. Mitra

Abstract: We report multi-wavelength study of a complex M-class solar eruptive flare that consists of three different sets of flare ribbons, viz. circular, parallel, and remote ribbons. Magnetic field modeling of source active region NOAA 12242 exhibits the presence of 3D null-point magnetic topology which encompasses an inner bipolar region. The event initiates with the faint signatures of the circular rib… ▽ More We report multi-wavelength study of a complex M-class solar eruptive flare that consists of three different sets of flare ribbons, viz. circular, parallel, and remote ribbons. Magnetic field modeling of source active region NOAA 12242 exhibits the presence of 3D null-point magnetic topology which encompasses an inner bipolar region. The event initiates with the faint signatures of the circular ribbon along with remote brightening right from the pre-flare phase which points toward the ongoing slow yet persistent null-point reconnection. We first detected flux cancellation and an associated brightening, which are likely signatures of tether-cutting reconnection that builds the flux rope near the polarity inversion line (PIL) of the inner bipolar region. In the next stage, with the onset of M8.7 flare, there is a substantial enhancement in the brightening of circular ribbon which essentially suggests an increase in the rate of ongoing null-point reconnection. Finally, the eruption of underlying flux rope triggers ``standard flare reconnection" beneath it producing an abrupt rise in the intensity of the parallel ribbons as well as enhancing the rate of null-point reconnection by external forcing. We show that within the the fan dome, the region with magnetic decay index n>1.5 borders the null-point QSL. Our analysis suggests that both the torus instability and the breakout model have played role toward the triggering mechanism for the eruptive flare. This event is a nice example of the dynamical evolution of a flux rope initially confined in a null-point topology, that subsequently activates and erupts with the progression of the circular -- cum -- parallel ribbon flare. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Journal ref: MNRAS, 501, 4703 (2021)

Successive occurrences of quasi-circular ribbon flares in a fan-spine-like configuration involving hyperbolic flux tube

Authors: Prabir K. Mitra, Bhuwan Joshi

Abstract: We present a comprehensive analysis of the formation and evolution of a fan-spine-like configuration that developed over a complex photospheric configuration where dispersed negative polarity regions were surrounded by positive polarity regions. This unique photospheric configuration, analogous to the geological "atoll" shape, hosted four homologous flares within its boundary. Computation of the d… ▽ More We present a comprehensive analysis of the formation and evolution of a fan-spine-like configuration that developed over a complex photospheric configuration where dispersed negative polarity regions were surrounded by positive polarity regions. This unique photospheric configuration, analogous to the geological "atoll" shape, hosted four homologous flares within its boundary. Computation of the degree of squashing factor (Q) maps clearly revealed an elongated region of high Q-values between the inner and outer spine-like lines, implying the presence of an hyperbolic flux tube (HFT). The coronal region associated with the photospheric atoll configuration was distinctly identified in the form of a diffused dome-shaped bright structure directly observed in EUV images. A filament channel resided near the boundary of the atoll region. The activation and eruption of flux ropes from the filament channel led to the onset of four eruptive homologous quasi-circular ribbon flares within an interval of $\approx$11 hours. During the interval of the four flares, we observed continuous decay and cancellation of negative polarity flux within the atoll region. Accordingly, the apparent length of the HFT gradually reduced to a null-point-like configuration before the fourth flare. Prior to each flare, we observed localised brightening beneath the filaments which, together with flux cancellation, provided support for the tether-cutting model of solar eruption. The analysis of magnetic decay index revealed favourable conditions for the eruption, once the pre-activated flux ropes attained the critical heights for torus instability. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 20 pages, 11 figures, 3 tables; accepted for publication in Monthly Notices of the Royal Astronomical Society

Results on Low-Mass Weakly Interacting Massive Particles from an 11 kg-day Target Exposure of DAMIC at SNOLAB

Authors: A. Aguilar-Arevalo, D. Amidei, D. Baxter, G. Cancelo, B. A. Cervantes Vergara, A. E. Chavarria, J. C. D'Olivo, J. Estrada, F. Favela-Perez, R. Gaior, Y. Guardincerri, E. W. Hoppe, T. W. Hossbach, B. Kilminster, I. Lawson, S. J. Lee, A. Letessier-Selvon, A. Matalon, P. Mitra, C. T. Overman, A. Piers, P. Privitera, K. Ramanathan, J. Da Rocha, Y. Sarkis , et al. (7 additional authors not shown)

Abstract: We present constraints on the existence of weakly interacting massive particles (WIMPs) from an 11 kg-day target exposure of the DAMIC experiment at the SNOLAB underground laboratory. The observed energy spectrum and spatial distribution of ionization events with electron-equivalent energies $>$200 eV$_{\rm ee}$ in the DAMIC CCDs are consistent with backgrounds from natural radioactivity. An exces… ▽ More We present constraints on the existence of weakly interacting massive particles (WIMPs) from an 11 kg-day target exposure of the DAMIC experiment at the SNOLAB underground laboratory. The observed energy spectrum and spatial distribution of ionization events with electron-equivalent energies $>$200 eV$_{\rm ee}$ in the DAMIC CCDs are consistent with backgrounds from natural radioactivity. An excess of ionization events is observed above the analysis threshold of 50 eV$_{\rm ee}$. While the origin of this low-energy excess requires further investigation, our data exclude spin-independent WIMP-nucleon scattering cross sections $σ_{χ-n}$ as low as $3\times 10^{-41}$ cm$^2$ for WIMPs with masses $m_χ$ from 7 to 10 GeV$c^{-2}$ . These results are the strongest constraints from a silicon target on the existence of WIMPs with $m_χ$$<$9 GeV$c^{-2}$ and are directly relevant to any dark matter interpretation of the excess of nuclear-recoil events observed by the CDMS silicon experiment in 2013. △ Less

Submitted 25 December, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 6 pages, 4 figures

Report number: FERMILAB-PUB-20-284-E

Journal ref: Phys. Rev. Lett. 125, 241803 (2020)

Eruptive-Impulsive Homologous M-class Flares Associated with Double-Decker Flux Rope Configuration in Mini-Sigmoid of NOAA 12673

Authors: Prabir K. Mitra, Bhuwan Joshi, Astrid M. Veronig, Ramesh Chandra, K. Dissauer, Thomas Wiegelmann

Abstract: We present a multiwavelength analysis of two homologous, short lived, impulsive flares of GOES class M1.4 and M7.3, that occurred from a very localized mini-sigmoid region within the active region NOAA 12673 on 2017 September 7. Both flares were associated with initial jet-like plasma ejection which for a brief amount of time moved toward east in a collimated manner before drastically changing dir… ▽ More We present a multiwavelength analysis of two homologous, short lived, impulsive flares of GOES class M1.4 and M7.3, that occurred from a very localized mini-sigmoid region within the active region NOAA 12673 on 2017 September 7. Both flares were associated with initial jet-like plasma ejection which for a brief amount of time moved toward east in a collimated manner before drastically changing direction toward southwest. Non-linear force-free field extrapolation reveals the presence of a compact double-decker flux rope configuration in the mini-sigmoid region prior to the flares. A set of open field lines originating near the active region which were most likely responsible for the anomalous dynamics of the erupted plasma, gave the earliest indication of an emerging coronal hole near the active region. The horizontal field distribution suggests a rapid decay of the field above the active region, implying high proneness of the flux rope system toward eruption. In view of the low coronal double-decker flux ropes and compact extreme ultra-violet (EUV) brightening beneath the filament along with associated photospheric magnetic field changes, our analysis supports the combination of initial tether-cutting reconnection and subsequent torus instability for driving the eruption. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 25 pages, 11 figures, 2 tables; accepted for publication in the Astrophysical Journal

Cosmogenic activation of silicon

Authors: R. Saldanha, R. Thomas, R. H. M. Tsang, A. E. Chavarria, R. Bunker, J. L. Burnett, S. R. Elliott, A. Matalon, P. Mitra, A. Piers, P. Privitera, K. Ramanathan, R. Smida

Abstract: The production of $^{3}$H, $^{7}$Be, and $^{22}$Na by interactions of cosmic-ray particles with silicon can produce radioactive backgrounds in detectors used to search for rare events. Through controlled irradiation of silicon CCDs and wafers with a neutron beam that mimics the cosmic-ray neutron spectrum, followed by direct counting, we determined that the production rate from cosmic-ray neutrons… ▽ More The production of $^{3}$H, $^{7}$Be, and $^{22}$Na by interactions of cosmic-ray particles with silicon can produce radioactive backgrounds in detectors used to search for rare events. Through controlled irradiation of silicon CCDs and wafers with a neutron beam that mimics the cosmic-ray neutron spectrum, followed by direct counting, we determined that the production rate from cosmic-ray neutrons at sea level is ($112 \pm 24$) atoms/(kg day) for $^{3}$H, ($8.1 \pm 1.9 $) atoms/(kg day) for $^{7}$Be, and ($43.0 \pm 7.1 $) atoms/(kg day) for $^{22}$Na. Complementing these results with the current best estimates of activation cross sections for cosmic-ray particles other than neutrons, we obtain a total sea-level cosmic-ray production rate of ($124 \pm 24$) atoms/(kg day) for $^{3}$H, ($9.4 \pm 2.0 $) atoms/(kg day) for $^{7}$Be, and ($49.6 \pm 7.3 $) atoms/(kg day) for $^{22}$Na. These measurements will help constrain background estimates and determine the maximum time that silicon-based detectors can remain unshielded during detector fabrication before cosmogenic backgrounds impact the sensitivity of next-generation rare-event searches. △ Less

Submitted 22 November, 2020; v1 submitted 20 July, 2020; originally announced July 2020.

Comments: 21 pages, 16 figures, 10 tables. Accepted version of the article published in Physical Review D

Journal ref: Phys. Rev. D 102, 102006 (2020)

The initial stage of cloud lightning imaged in high-resolution

Authors: O. Scholten, B. M. Hare, J. Dwyer, C. Sterpka, I. Kolmašová, O. Santolík, R. Lán, L. Uhlíř, S. Buitink, A. Corstanje, H. Falcke, T. Huege, J. R. Hörandel, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, A. Pel, J. P. Rachen, T. N. G. Trinh, S. ter Veen, S. Thoudam, T. Winchen

Abstract: With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the ch… ▽ More With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Submitted to Journal of geophysics research: Atmospheres

Radio emission from negative lightning leader steps reveals inner meter-scale structure

Authors: B. M. Hare, O. Scholten, J. Dwyer, U. Ebert, S. Nijdam, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, T. Huege, J. R. Hörandel, G. K. Krampah, P. Mitra, K. Mulrey, B. Neijzen, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, T. N. G. Trinh, S. ter Veen, T. Winchen

Abstract: We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ μ$s). This is evidence for streamer formation during corona flashes that occur with each leader st… ▽ More We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ μ$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 5 pages

Journal ref: Phys. Rev. Lett. 124, 105101 (2020)

Solar Radio Observation Using CALLISTO at the USO/PRL, Udaipur

Authors: Kushagra Upadhyay, Bhuwan Joshi, Prabir K. Mitra, Ramit Bhattacharyya, Divya Oberoi, Christian Monstein

Abstract: This paper presents a detailed description of various subsystems of CALLISTO solar radio spectrograph installed at the USO/PRL. In the front-end system, a log periodic dipole antenna (LPDA) is designed for the frequency range of 40-900 MHz. In this paper LPDA design, its modifications, and simulation results are presented. We also present some initial observations taken by CALLISTO at Udaipur. This paper presents a detailed description of various subsystems of CALLISTO solar radio spectrograph installed at the USO/PRL. In the front-end system, a log periodic dipole antenna (LPDA) is designed for the frequency range of 40-900 MHz. In this paper LPDA design, its modifications, and simulation results are presented. We also present some initial observations taken by CALLISTO at Udaipur. △ Less

Submitted 3 July, 2020; originally announced July 2020.

Comments: Published in IEEE. https://ieeexplore.ieee.org/document/9118669

Journal ref: 2019 IEEE MTT-S International Microwave and RF Conference (IMARC), Electronic ISBN: 978-1-7281-4040-7, Print on Demand(PoD) ISBN: 978-1-7281-4041-4

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

Authors: P. Mitra, A. Bonardi, A. Corstanje, S. Buitink, G. K Krampah, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, O. Scholten, S. ter Veen, T. N. G. Trinh, T. Winchen

Abstract: The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. How… ▽ More The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.02859

On the cosmic-ray energy scale of the LOFAR radio telescope

Authors: K. Mulrey, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, G. K. Krampah, P. Mitra, A. Nelles, H. Pandya, J. P. Rachen, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-ev… ▽ More Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-event and absolute scale uncertainties. The resulting energies reconstructed with each method are shown to be in good agreement, and because the radio-based reconstructed energy has smaller uncertainty on an event-to-event basis, LOFAR analyses will use that technique in the future. Second, we present the radiation energy of air showers measured at LOFAR and demonstrate how radiation energy can be used to compare the energy scales of different experiments. The radiation energy scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LORA particle-based energy scale and that of the Pierre Auger Observatory. They are shown to agree to within (6$\pm$20)% for a radiation energy of 1 MeV, where the uncertainty on the comparison is dominated by the antenna calibrations of each experiment. This study motivates the development of a portable radio array which will be used to cross-calibrate the energy scales of different experiments using radiation energy and the same antennas, thereby significantly reducing the uncertainty on the comparison. △ Less

Submitted 29 September, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Development of a Confined Circular-cum-parallel Ribbon Flare and Associated Pre-flare Activity

Authors: Pooja Devi, Bhuwan Joshi, Ramesh Chandra, Prabir K. Mitra, Astrid M. Veronig3, Reetika Joshi

Abstract: We study a complex GOES M1.1 circular ribbon flare and related pre-flare activity on 26 January 2015 [SOL26-01-2015] in solar active region NOAA 12268. This flare activity was observed by the AIA on board SDO and the RHESSI. The examination of photospheric magnetograms during the extended period, prior to the event, suggests the successive development of a so-called 'anemone' type magnetic configu… ▽ More We study a complex GOES M1.1 circular ribbon flare and related pre-flare activity on 26 January 2015 [SOL26-01-2015] in solar active region NOAA 12268. This flare activity was observed by the AIA on board SDO and the RHESSI. The examination of photospheric magnetograms during the extended period, prior to the event, suggests the successive development of a so-called 'anemone' type magnetic configuration. NLFFF extrapolation reveals a fan-spine magnetic configuration with the presence of a coronal null-point. We found that the pre-flare activity in the active region starts ~15 min prior to the main flare in the form of localized bright patches at two locations. A comparison of locations and spatial structures of the pre-flare activity with magnetic configuration of the corresponding region suggests onset of magnetic reconnection at the null-point along with the low-atmosphere magnetic reconnection caused by the emergence and the cancellation of the magnetic flux. The main flare of M1.1 class is characterized by the formation of a well-developed circular ribbon along with a region of remote brightening. Remarkably, a set of relatively compact parallel ribbons formed inside the periphery of the circular ribbon which developed lateral to the brightest part of the circular ribbon. During the peak phase of the flare, a coronal jet is observed at the north-east edge of the circular ribbon which suggests interchange reconnection between large-scale field lines and low-lying closed field lines. Our investigation suggests a combination of two distinct processes in which ongoing pre-flare null-point reconnection gets further intensified as the confined eruption along with jet activity proceeded from within the circular ribbon region which results to the formation of inner parallel ribbons and corresponding post-reconnection arcade. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 24 pages, 11 figures, Accepted for publication in Solar Physics

HXR emission from an activated flux rope and subsequent evolution of an eruptive long duration solar flare

Authors: Suraj Sahu, Bhuwan Joshi, Prabir K. Mitra, Astrid M. Veronig, V. Yurchyshyn

Abstract: In this paper, we present a comprehensive study of the evolutionary phases of a major M6.6 long duration event (LDE) with special emphasize on its pre-flare phase. The event occurred in NOAA 12371 on 2015 June 22. A remarkable aspect of the event was an active pre-flare phase lasting for about an hour during which a hot EUV coronal channel was in build-up stage and displayed co-spatial hard X-ray… ▽ More In this paper, we present a comprehensive study of the evolutionary phases of a major M6.6 long duration event (LDE) with special emphasize on its pre-flare phase. The event occurred in NOAA 12371 on 2015 June 22. A remarkable aspect of the event was an active pre-flare phase lasting for about an hour during which a hot EUV coronal channel was in build-up stage and displayed co-spatial hard X-ray (HXR) emission up to energies of 25 keV. As such, this is the first evidence of HXR coronal channel. The coronal magnetic field configuration based on NLFFF modeling clearly exhibited a magnetic flux rope (MFR) oriented along the polarity inversion line (PIL) and co-spatial with the coronal channel. We observed significant changes in the AR's photospheric magnetic field during an extended period of $\approx$42 hours in the form of rotation of sunspots, moving magnetic features, and flux cancellation along the PIL. Prior to the flare onset, the MFR underwent a slow rise phase ($\approx$14 km s$^{-1}$) for $\approx$12 min which we attribute to the faster build-up and activation of the MFR by tether-cutting reconnection occurring at multiple locations along the MFR itself. The sudden transition in the kinematic evolution of the MFR from the phase of slow to fast rise ($\approx$109 km s$^{-1}$ with acceleration $\approx$110 m s$^{-2}$) precisely divides the pre-flare and impulsive phase of the flare, which points toward the feedback process between the early dynamics of the eruption and the strength of the flare magnetic reconnection. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted for publication in The Astrophysical Journal

Identification of Pre-flare Processes and Their Possible Role in Driving a Large-scale Flux Rope Eruption with Complex M-class Flare in the Active Region NOAA 12371

Authors: Prabir K. Mitra, Bhuwan Joshi, Avijeet Prasad

Abstract: In this article, we study the origin of precursor flare activity and investigate its role towards triggering the eruption of a flux rope which resulted into a dual-peak M-class flare (SOL2015-06-21T02:36) in the active region NOAA 12371. The flare evolved in two distinct phases with peak flux levels of M2.1 and M2.6 at an interval of $\approx$54 min. The active region exhibited striking moving mag… ▽ More In this article, we study the origin of precursor flare activity and investigate its role towards triggering the eruption of a flux rope which resulted into a dual-peak M-class flare (SOL2015-06-21T02:36) in the active region NOAA 12371. The flare evolved in two distinct phases with peak flux levels of M2.1 and M2.6 at an interval of $\approx$54 min. The active region exhibited striking moving magnetic features (MMFs) along with sunspot rotation. Non-linear force free field (NLFFF) modelling of the active region corona reveals a magnetic flux rope along the polarity inversion line in the trailing sunspot group which is observationally manifested by the co-spatial structures of an active region filament and a hot channel identified in the 304 and 94 Å images, respectively, from the Atmospheric Imaging Assembly (AIA). The active region underwent a prolonged phase of flux enhancement followed by a relatively shorter period of flux cancellation prior to the onset of the flare which led to the build up and activation of the flux rope. Extreme ultra-violet (EUV) images reveal localised and structured pre-flare emission, from the region of MMFs, adjacent to the location of the main flare. Our analysis reveals strong, localised regions of photospheric currents of opposite polarities at the precursor location, thereby making the region susceptible to small-scale magnetic reconnection. Precursor reconnection activity from this location most likely induced a slipping reconnetion towards the northern leg of the hot channel which led to the destabilization of the flux rope. The application of magnetic virial theorem suggests that there was an overall growth of magnetic free energy in the active region during the prolonged pre-flare phase which decayed rapidly after the hot channel eruption and its successful transformation into a halo coronal mass ejection (CME). △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 31 pages, 14 figures, 1 table. This is a pre-print of an article published in Solar Physics. The final authenticated version is available online at: https://doi.org/10.1007/s11207-020-1596-2

Contributions of the LOFAR Cosmic Ray Key Science Project to the 36th International Cosmic Ray Conference (ICRC 2019)

Authors: A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Horandel, T. Huege, G. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, O. Scholten, S. ter Veen, T. N. G. Trinh, T. Winchen

Abstract: This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection.… ▽ More This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate. △ Less

Submitted 7 November, 2019; originally announced November 2019.

Comments: 42 pages, 17 figures, Collection of papers presented at the ICRC2019

Pre-flare processes, flux rope activation, large-scale eruption and associated X-class flare from the active region NOAA 11875

Authors: Prabir K. Mitra, Bhuwan Joshi

Abstract: We present a multi-wavelength analysis of the eruption of a hot coronal channel associated with an X1.0 flare (SOL2013-10-28T02:03) from the active region NOAA 11875 by combining observations from AIA/SDO, HMI/SDO, RHESSI, and HiRAS. EUV images at high coronal temperatures indicated the presence of a hot channel at the core of the active region from the early pre-flare phase evidencing the pre-exi… ▽ More We present a multi-wavelength analysis of the eruption of a hot coronal channel associated with an X1.0 flare (SOL2013-10-28T02:03) from the active region NOAA 11875 by combining observations from AIA/SDO, HMI/SDO, RHESSI, and HiRAS. EUV images at high coronal temperatures indicated the presence of a hot channel at the core of the active region from the early pre-flare phase evidencing the pre-existence of a quasi-stable magnetic flux rope. The hot channel underwent an activation phase after a localized and prolonged pre-flare event occurring adjacent to one of its footpoints. Subsequently, the flux rope continued to rise slowly for $\approx$16 min during which soft X-ray flux gradually built-up characterizing a distinct precursor phase. The flux rope transitioned from the state of slow rise to the eruptive motion with the onset of the impulsive phase of the X1.0 flare. The eruptive expansion of the hot channel is accompanied by a series of type III radio bursts in association with impulsive rise of strong hard X-ray non-thermal emissions that included explicit hard X-ray sources of energies up to $\approx$50 keV from the coronal loops and $\approx$100 keV from their footpoint locations. Our study contains evidence that pre-flare activity occurring within the spatial extent of a stable flux rope can destabilize it toward eruption. Moreover, sudden transition of the flux rope from the state of slow rise to fast acceleration precisely bifurcated the precursor and the impulsive phases of the flare which points toward a feedback relationship between early CME dynamics and the strength of the large-scale magnetic reconnection. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: Accepted for publication in The Astrophysical Journal. 20 pages, 10 figures, 1 table

Constraints on Light Dark Matter Particles Interacting with Electrons from DAMIC at SNOLAB

Authors: A. Aguilar-Arevalo, D. Amidei, D. Baxter, G. Cancelo, B. A. Cervantes Vergara, A. E. Chavarria, E. Darragh-Ford, J. R. T. de Mello Neto, J. C. D'Olivo, J. Estrada, R. Gaïor, Y. Guardincerri, T. W. Hossbach, B. Kilminster, I. Lawson, S. J. Lee, A. Letessier-Selvon, A. Matalon, V. B. B. Mello, P. Mitra, Y. S. Mobarak, J. Molina, S. Paul, A. Piers, P. Privitera , et al. (9 additional authors not shown)

Abstract: We report direct-detection constraints on light dark matter particles interacting with electrons. The results are based on a method that exploits the extremely low levels of leakage current of the DAMIC detector at SNOLAB of 2-6$\times$10$^{-22}$ A cm$^{-2}$. We evaluate the charge distribution of pixels that collect $<10~\rm{e^-}$ for contributions beyond the leakage current that may be attribute… ▽ More We report direct-detection constraints on light dark matter particles interacting with electrons. The results are based on a method that exploits the extremely low levels of leakage current of the DAMIC detector at SNOLAB of 2-6$\times$10$^{-22}$ A cm$^{-2}$. We evaluate the charge distribution of pixels that collect $<10~\rm{e^-}$ for contributions beyond the leakage current that may be attributed to dark matter interactions. Constraints are placed on so-far unexplored parameter space for dark matter masses between 0.6 and 100 MeV$c^{-2}$. We also present new constraints on hidden-photon dark matter with masses in the range $1.2$-$30$ eV$c^{-2}$. △ Less

Submitted 8 April, 2020; v1 submitted 29 July, 2019; originally announced July 2019.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 123, 181802 (2019)

Velocity independent constraints on spin-dependent DM-nucleon interactions from IceCube and PICO

Authors: IceCube Collaboration, M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, C. Alispach, K. Andeen, T. Anderson, I. Ansseau, G. Anton, C. Argüelles, J. Auffenberg, S. Axani, P. Backes, H. Bagherpour, X. Bai, A. Balagopal V., A. Barbano, S. W. Barwick, B. Bastian, V. Baum, S. Baur, R. Bay , et al. (398 additional authors not shown)

Abstract: Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from… ▽ More Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $σ_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $τ^+ τ^-$ ($b\bar{b}$) with a local density of $ρ_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $ρ_{\mathrm{DM}}$ and are independent of the DM velocity distribution. △ Less

Submitted 25 May, 2020; v1 submitted 29 July, 2019; originally announced July 2019.

Comments: 9 pages, 3 figures, 1 table. To appear in Eur.Phys.J. C

Journal ref: Eur. Phys. J. C 80 (2020) 819

Data-Driven Modeling of Electron Recoil Nucleation in PICO C$_3$F$_8$ Bubble Chambers

Authors: C. Amole, M. Ardid, I. J. Arnquist, D. M. Asner, D. Baxter, E. Behnke, M. Bressler, B. Broerman, G. Cao, C. J. Chen, S. Chen, U. Chowdhury, K. Clark, J. I. Collar, P. S. Cooper, C. B. Coutu, C. Cowles, M. Crisler, G. Crowder, N. A. Cruz-Venegas, C. E. Dahl, M. Das, S. Fallows, J. Farine, R. Filgas , et al. (54 additional authors not shown)

Abstract: The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by… ▽ More The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$. △ Less

Submitted 25 November, 2020; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: 18 pages, 12 figures

Journal ref: Phys. Rev. D 100, 082006 (2019)

Cosmic Ray Physics with the LOFAR Radio Telescope

Authors: T Winchen, A Bonardi, S Buitink, A Corstanje, H Falcke, B M Hare, J R Hörandel, P Mitra, K Mulrey, A Nelles, J P Rachen, L Rossetto, P Schellart, O Scholten, S ter Veen, S Thoudam, T N G Trinh

Abstract: The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a d… ▽ More The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Status of the Lunar Detection Mode for Cosmic Particles of LOFAR

Authors: T. Winchen, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh

Abstract: Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum fr… ▽ More Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Calibration of the LOFAR low-band antennas using the Galaxy and a model of the signal chain

Authors: K. Mulrey, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, P. Mitra, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information abo… ▽ More The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: 23 pages, 10 figures

Dark Matter Search Results from the Complete Exposure of the PICO-60 C$_3$F$_8$ Bubble Chamber

Authors: C. Amole, M. Ardid, I. J. Arnquist, D. M. Asner, D. Baxter, E. Behnke, M. Bressler, B. Broerman, G. Cao, C. J. Chen, U. Chowdhury, K. Clark, J. I. Collar, P. S. Cooper, C. B. Coutu, C. Cowles, M. Crisler, G. Crowder, N. A. Cruz-Venegas, C. E. Dahl, M. Das, S. Fallows, J. Farine, I. Felis, R. Filgas , et al. (47 additional authors not shown)

Abstract: Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total backgr… ▽ More Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. D 100, 022001 (2019)

Successive flux rope eruptions from $δ$-Sunspots region of NOAA 12673 and associated X-class eruptive flares on 2017 September 6

Authors: Prabir K. Mitra, Bhuwan Joshi, Avijeet Prasad, Astrid M. Veronig, R. Bhattacharyya

Abstract: In this paper, we present a multi-wavelength analysis of two X-class solar eruptive flares of classes X2.2 and X9.3 that occurred in the sigmoidal active region NOAA 12673 on 2017 September 6, by combining observations of Atmospheric Imaging Assembly and Helioseismic Magnetic Imager instruments on board the Solar Dynamics Observatory. On the day of the reported activity, the photospheric structure… ▽ More In this paper, we present a multi-wavelength analysis of two X-class solar eruptive flares of classes X2.2 and X9.3 that occurred in the sigmoidal active region NOAA 12673 on 2017 September 6, by combining observations of Atmospheric Imaging Assembly and Helioseismic Magnetic Imager instruments on board the Solar Dynamics Observatory. On the day of the reported activity, the photospheric structure of the active region displayed a very complex network of $δ$-sunspots that gave rise to the formation of a coronal sigmoid observed in the hot EUV channels. Both X-class flares initiated from the core of the sigmoid sequentially within an interval of $\sim$3 hours and progressed as a single "sigmoid--to--arcade" event. Differential emission measure analysis reveals strong heating of plasma at the core of the active region right from the pre-flare phase which further intensified and spatially expanded during each event. The identification of a pre-existing magnetic null by non-force-free-field modeling of the coronal magnetic fields at the location of early flare brightenings and remote faint ribbon-like structures during the pre-flare phase, which were magnetically connected with the core region, provide support for the breakout model of solar eruption. The magnetic extrapolations also reveal flux rope structures prior to both flares which are subsequently supported by the observations of the eruption of hot EUV channels. The second X-class flare diverged from the standard flare scenario in the evolution of two sets of flare ribbons, that are spatially well separated, providing firm evidence of magnetic reconnections at two coronal heights. △ Less

Submitted 31 October, 2018; originally announced October 2018.

Comments: 24 pages, 13 figures. Accepted for publication by The Astrophysical Journal (ApJ)

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017)

Authors: A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017) Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017) △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: compendium of 8 proceedings for the ICRC 2017

Observational and model analysis of a two-ribbon flare possibly induced by a neighbouring blowout jet

Authors: Bhuwan Joshi, Julia K. Thalmann, Prabir K. Mitra, Ramesh Chandra, Astrid M. Veronig

Abstract: In this paper, we present unique observations of a blowout coronal jet that possibly triggered a two-ribbon confined C1.2 flare in a bipolar solar active region NOAA 12615 on 2016 December 5. The jet activity initiates at chromospheric/transition-region heights with a small brightening that eventually grows in a larger volume with well developed standard morphological jet features, viz., base and… ▽ More In this paper, we present unique observations of a blowout coronal jet that possibly triggered a two-ribbon confined C1.2 flare in a bipolar solar active region NOAA 12615 on 2016 December 5. The jet activity initiates at chromospheric/transition-region heights with a small brightening that eventually grows in a larger volume with well developed standard morphological jet features, viz., base and spire. The spire widens up with a collimated eruption of cool and hot plasma components, observed in the 304 and 94 A channels of AIA, respectively. The speed of the plasma ejection, which forms the jet's spire, was higher for the hot component (~200 km/s) than the cooler one (~130 km/s). The NLFF model of coronal fields at pre- and post-jet phases successfully reveal opening of previously closed magnetic field lines with a rather inclined/low-lying jet structure. The peak phase of the jet emission is followed by the development of a two-ribbon flare that shows coronal loop emission in HXRs up to ~25 keV energy. The coronal magnetic fields rooted at the location of EUV flare ribbons, derived from the NLFF model, demonstrate the pre-flare phase to exhibit an "X-type" configuration while the magnetic fields at the post-flare phase are more or less parallel oriented. The comparisons of multi-wavelength measurements with the magnetic field extrapolations suggest that the jet activity likely triggered the two-ribbon flare by perturbing the field in the interior of the active region. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 18 pages, 12 figures, accepted for publication in The Astrophysical Journal (ApJ)

Thunderstorm electric fields probed by extensive air showers through their polarized radio emission

Authors: T. N. G. Trinh, O. Scholten, A. Bonardi, S. Buitink, A. Corstanje, U. Ebert, J. E. Enriquez, H. Falcke, J. R. Horandel, B. M. Hare, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, C. Rutjes, P. Schellart, S. Thoudam, S. ter Veen, T. Winchen

Abstract: We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function o… ▽ More We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed. △ Less

Submitted 20 March, 2017; v1 submitted 14 March, 2017; originally announced March 2017.

Dark Matter Search Results from the PICO-60 C$_3$F$_8$ Bubble Chamber

Authors: C. Amole, M. Ardid, I. J. Arnquist, D. M. Asner, D. Baxter, E. Behnke, P. Bhattacharjee, H. Borsodi, M. Bou-Cabo, P. Campion, G. Cao, C. J. Chen, U. Chowdhury, K. Clark, J. I. Collar, P. S. Cooper, M. Crisler, G. Crowder, C. E. Dahl, M. Das, S. Fallows, J. Farine, I. Felis, R. Filgas, F. Girard , et al. (37 additional authors not shown)

Abstract: New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less th… ▽ More New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results. △ Less

Submitted 2 August, 2017; v1 submitted 24 February, 2017; originally announced February 2017.

Comments: 6 pages, 6 figures

Journal ref: Phys. Rev. Lett. 118, 251301 (2017)

The effect of the atmospheric refractive index on the radio signal of extensive air showers

Authors: A. Corstanje, A. Bonardi, S. Buitink, H. Falcke, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, G. Trinh, T. Winchen

Abstract: For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are… ▽ More For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 13 pages, 5 figures. Accepted for publication in Astroparticle Physics