Showing 1–50 of 73 results for author: Richardson, G

Towards quantum-enhanced long-baseline optical/near-IR interferometry

Authors: Jayadev K. Rajagopal, Ryan M. Lau, Isack Padilla, Stephen T. Ridgway, Chaohan Cui, Brittany McClinton, Aqil Sajjad, Stuartt Corder, Mark Rawlings, Fredrik Rantakyro, J. Gabriel Richardson, Amit Ashok, Saikat Guha

Abstract: Microarcsecond resolutions afforded by an optical-NIR array with kilometer-baselines would enable breakthrough science. However significant technology barriers exist in transporting weakly coherent photon states over these distances: primarily photon loss and phase errors. Quantum telescopy, using entangled states to link spatially separated apertures, offers a possible solution to the loss of pho… ▽ More Microarcsecond resolutions afforded by an optical-NIR array with kilometer-baselines would enable breakthrough science. However significant technology barriers exist in transporting weakly coherent photon states over these distances: primarily photon loss and phase errors. Quantum telescopy, using entangled states to link spatially separated apertures, offers a possible solution to the loss of photons. We report on an initiative launched by NSF NOIRLab in collaboration with the Center for Quantum Networks and Arizona Quantum Initiative at the University of Arizona, Tucson, to explore these concepts further. A brief description of the quantum concepts and a possible technology roadmap towards a quantum-enhanced very long baseline optical-NIR interferometric array is presented. An on-sky demonstration of measuring spatial coherence of photons with apertures linked through the simplest Gottesman protocol over short baselines and with limited phase fluctuations is envisaged as the first step. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Proceeding of SPIE Conference "Astronomical Telescopes + Instrumentation" (June 2024)

Report number: Paper No. 13095-58

Imaging of single barium atoms in a second matrix site in solid xenon for barium tagging in a $^{136}$Xe double beta decay experiment

Authors: M. Yvaine, D. Fairbank, J. Soderstrom, C. Taylor, J. Stanley, T. Walton, C. Chambers, A. Iverson, W. Fairbank, S. Al Kharusi, A. Amy, E. Angelico, A. Anker, I. J. Arnquist, A. Atencio, J. Bane, V. Belov, E. P. Bernard, T. Bhatta, A. Bolotnikov, J. Breslin, P. A. Breur, J. P. Brodsky, E. Brown, T. Brunner , et al. (112 additional authors not shown)

Abstract: Neutrinoless double beta decay is one of the most sensitive probes for new physics beyond the Standard Model of particle physics. One of the isotopes under investigation is $^{136}$Xe, which would double beta decay into $^{136}$Ba. Detecting the single $^{136}$Ba daughter provides a sort of ultimate tool in the discrimination against backgrounds. Previous work demonstrated the ability to perform s… ▽ More Neutrinoless double beta decay is one of the most sensitive probes for new physics beyond the Standard Model of particle physics. One of the isotopes under investigation is $^{136}$Xe, which would double beta decay into $^{136}$Ba. Detecting the single $^{136}$Ba daughter provides a sort of ultimate tool in the discrimination against backgrounds. Previous work demonstrated the ability to perform single atom imaging of Ba atoms in a single-vacancy site of a solid xenon matrix. In this paper, the effort to identify signal from individual barium atoms is extended to Ba atoms in a hexa-vacancy site in the matrix and is achieved despite increased photobleaching in this site. Abrupt fluorescence turn-off of a single Ba atom is also observed. Significant recovery of fluorescence signal lost through photobleaching is demonstrated upon annealing of Ba deposits in the Xe ice. Following annealing, it is observed that Ba atoms in the hexa-vacancy site exhibit antibleaching while Ba atoms in the tetra-vacancy site exhibit bleaching. This may be evidence for a matrix site transfer upon laser excitation. Our findings offer a path of continued research toward tagging of Ba daughters in all significant sites in solid xenon. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 9 pages, 8 figures

Supernova Electron-Neutrino Interactions with Xenon in the nEXO Detector

Authors: nEXO Collaboration, S. Hedges, S. Al Kharusi, E. Angelico, J. P. Brodsky, G. Richardson, S. Wilde, A. Amy, A. Anker, I. J. Arnquist, P. Arsenault, A. Atencio, I. Badhrees, J. Bane, V. Belov, E. P. Bernard, T. Bhatta, A. Bolotnikov, J. Breslin, P. A. Breur, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Q. Cao , et al. (121 additional authors not shown)

Abstract: Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-beta decay detector (~5-tonne, 90% ${}^{136}$Xe, 10% ${}^{134}$Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the MARLEY event generator. We find good agreement between MARLEY's predictions and existing theor… ▽ More Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-beta decay detector (~5-tonne, 90% ${}^{136}$Xe, 10% ${}^{134}$Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the MARLEY event generator. We find good agreement between MARLEY's predictions and existing theoretical calculations of the inclusive cross sections at supernova neutrino energies. The interactions modeled by MARLEY were simulated within the nEXO simulation framework and were run through an example reconstruction algorithm to determine the detector's efficiency for reconstructing these events. The simulated data, incorporating the detector response, were used to study the ability of nEXO to reconstruct the incident electron-neutrino spectrum and these results were extended to a larger xenon detector of the same isotope enrichment. We estimate that nEXO will be able to observe electron-neutrino interactions with xenon from supernovae as far as 5 to 8 kpc from earth, while the ability to reconstruct incident electron-neutrino spectrum parameters from observed interactions in nEXO is limited to closer supernovae. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 17 pages, 16 figures

Report number: LLNL-JRNL-864783-DRAFT

Report from the Workshop on Xenon Detector $0νββ$ Searches: Steps Towards the Kilotonne Scale

Authors: A. Anker, A. Avasthi, M. Brodeur, T. Brunner, N. K. Byrnes, N. R. Catarineu, A. Cottle, P. Englezos, W. Fairbank, D. González Díaz, R. Guenette, S. J. Haselschwardt, S. Hedges, M. Heffner, J. D. Holt, A. Jamil, B. J. P. Jones, N. Kawada, S. Leardini, B. G. Lenardo, A. Marc, J. Masbou, K. Mistry, B. Mong, B. Monreal , et al. (13 additional authors not shown)

Abstract: These proceedings summarize the program and discussions of the ``Workshop on Xenon Detector $0νββ$ Searches: Steps Towards the Kilotonne Scale'' held on October 25-27 2023 at SLAC National Accelerator Laboratory. This workshop brought together experts from the communities of neutrinoless double-beta decay and dark matter detection, to discuss paths forward for the realization of monolithic experim… ▽ More These proceedings summarize the program and discussions of the ``Workshop on Xenon Detector $0νββ$ Searches: Steps Towards the Kilotonne Scale'' held on October 25-27 2023 at SLAC National Accelerator Laboratory. This workshop brought together experts from the communities of neutrinoless double-beta decay and dark matter detection, to discuss paths forward for the realization of monolithic experiments with xenon approaching the kilotonne scale. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 7 pages, 2 figures. Proceedings

Improved modelling of SEP event onset within the WSA-Enlil-SEPMOD framework

Authors: Erika Palmerio, Janet G. Luhmann, M. Leila Mays, Ronald M. Caplan, David Lario, Ian G. Richardson, Kathryn Whitman, Christina O. Lee, Beatriz Sánchez-Cano, Nicolas Wijsen, Yan Li, Carlota Cardoso, Marco Pinto, Daniel Heyner, Daniel Schmid, Hans-Ulrich Auster, David Fischer

Abstract: Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories, but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best-suited to capture a particular phase of an SEP event, rather… ▽ More Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories, but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best-suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain -- usually set at a heliocentric distance of 20-30 $R_{\odot}$. This means that particle acceleration at CME-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a "fixed source" of particle injection at the outer boundary of the coronal domain of the coupled WSA-Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is SEPMOD, a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD's fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere -- Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 31 pages, 8 figures, 4 tables, accepted for publication in Journal of Space Weather and Space Climate

Solar Energetic Particle-Associated Coronal Mass Ejections Observed by the Mauna Loa Solar Observatory Mk3 and Mk4 Coronameters

Authors: I. G. Richardson, O. C. St Cyr, J. T. Burkepile, H. Xie, B. J. Thompson

Abstract: We report on the first comprehensive study of the coronal mass ejections (CMEs) associated with $\sim$25 MeV solar energetic proton (SEP) events in 1980-2013 observed in the low/inner corona by the Mauna Loa Solar Observatory (MLSO) Mk3 and Mk4 coronameters. Where possible, these observations are combined with spacebased observations from the Solar Maximum Mission C/P, P78-1 SOLWIND or SOHO/LASCO… ▽ More We report on the first comprehensive study of the coronal mass ejections (CMEs) associated with $\sim$25 MeV solar energetic proton (SEP) events in 1980-2013 observed in the low/inner corona by the Mauna Loa Solar Observatory (MLSO) Mk3 and Mk4 coronameters. Where possible, these observations are combined with spacebased observations from the Solar Maximum Mission C/P, P78-1 SOLWIND or SOHO/LASCO coronagraphs. The aim of the study is to understand directly-measured (rather than inferred from proxies) CME motions in the low to middle corona and their association with SEP acceleration, and hence attempt to identify early signatures that are characteristic of SEP acceleration in ground-based CME observations that may be used to warn of impending SEP events. Although we find that SEP events are associated with CMEs that are on average faster and wider than typical CMEs observed by MLSO, a major challenge turns out to be determining reliable estimates of the CME dynamics in the low corona from the 3-minute cadence Mk3/4 observations since different analysis techniques can produce inconsistent results. This complicates the assessment of what early information on a possible SEP event is available from these low coronal observations △ Less

Submitted 24 August, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: To be published in Solar Physics

CME Propagation Through the Heliosphere: Status and Future of Observations and Model Development

Authors: M. Temmer, C. Scolini, I. G. Richardson, S. G. Heinemann, E. Paouris, A. Vourlidas, M. M. Bisi, writing teams, :, N. Al-Haddad, T. Amerstorfer, L. Barnard, D. Buresova, S. J. Hofmeister, K. Iwai, B. V. Jackson, R. Jarolim, L. K. Jian, J. A. Linker, N. Lugaz, P. K. Manoharan, M. L. Mays, W. Mishra, M. J. Owens, E. Palmerio , et al. (9 additional authors not shown)

Abstract: The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field follo… ▽ More The ISWAT clusters H1+H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. SIRs, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field followed by high-speed streams that recur at the ca. 27 day solar rotation period. Short-term reconfigurations of the lower coronal magnetic field generate flare emissions and provide the energy to accelerate enormous amounts of magnetised plasma and particles in the form of CMEs into interplanetary space. The dynamic interplay between these phenomena changes the configuration of interplanetary space on various temporal and spatial scales which in turn influences the propagation of individual structures. While considerable efforts have been made to model the solar wind, we outline the limitations arising from the rather large uncertainties in parameters inferred from observations that make reliable predictions of the structures impacting Earth difficult. Moreover, the increased complexity of interplanetary space as solar activity rises in cycle 25 is likely to pose a challenge to these models. Combining observational and modeling expertise will extend our knowledge of the relationship between these different phenomena and the underlying physical processes, leading to improved models and scientific understanding and more-reliable space-weather forecasting. The current paper summarizes the efforts and progress achieved in recent years, identifies open questions, and gives an outlook for the next 5-10 years. It acts as basis for updating the existing COSPAR roadmap by Schrijver+ (2015), as well as providing a useful and practical guide for peer-users and the next generation of space weather scientists. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Advances in Space Research

Statistical Relationship Between Long-duration High-Energy Gamma-Ray Emission and Solar Energetic Particles

Authors: Alessandro Bruno, Georgia A. de Nolfo, James M. Ryan, Ian G. Richardson, Silvia Dalla

Abstract: Large solar eruptions are often associated with long-duration gamma-ray emission extending well above 100 MeV. While this phenomenon is known to be caused by high-energy ions interacting with the solar atmosphere, the underlying dominant acceleration process remains under debate. Potential mechanisms include continuous acceleration of particles trapped within large coronal loops or acceleration at… ▽ More Large solar eruptions are often associated with long-duration gamma-ray emission extending well above 100 MeV. While this phenomenon is known to be caused by high-energy ions interacting with the solar atmosphere, the underlying dominant acceleration process remains under debate. Potential mechanisms include continuous acceleration of particles trapped within large coronal loops or acceleration at coronal mass ejection (CME)-driven shocks, with subsequent back-propagation towards the Sun. As a test of the latter scenario, previous studies have explored the relationship between the inferred particle population producing the high-energy gamma-rays, and the population of solar energetic particles (SEPs) measured in situ. However, given the significant limitations on available observations, these estimates unavoidably rely on a number of assumptions. In an effort to better constrain theories of the gamma-ray emission origin, we re-examine the calculation uncertainties and how they influence the comparison of these two proton populations. We show that, even accounting for conservative assumptions related to gamma-ray flare, SEP event and interplanetary scattering modeling, their statistical relationship is only poorly/moderately significant. However, though the level of correlation is of interest, it does not provide conclusive evidence for or against a causal connection. The main result of this investigation is that the fraction of the shock-accelerated protons required to account for the gamma-ray observations is >20-40% for six of the fourteen eruptions analyzed. Such high values argue against current CME-shock origin models, predicting a <2% back-precipitation, hence the computed numbers of high-energy SEPs appear to be greatly insufficient to sustain the measured gamma-ray emission. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: Accepted for publication in the Astrophysical Journal

The effect of the ambient solar wind medium on a CME-driven shock and the associated gradual solar energetic particle event

Authors: Nicolas Wijsen, David Lario, Beatriz Sánchez-Cano, Immanuel C. Jebaraj, Nina Dresing, Ian G. Richardson, Angels Aran, Athanasios Kouloumvakos, Zheyi Ding, Antonio Niemela, Erika Palmerio, Fernando Carcaboso, Rami Vainio, Alexandr Afanasiev, Marco Pinto, Daniel Pacheco, Stefaan Poedts, Daniel Heyner

Abstract: We present simulation results of a gradual solar energetic particle (SEP) event detected on 2021 October 9 by multiple spacecraft, including BepiColombo (Bepi) and near-Earth spacecraft such as the Advanced Composition Explorer (ACE). A peculiarity of this event is that the presence of a high speed stream (HSS) affected the low-energy ion component ($\lesssim 5$ MeV) of the gradual SEP event at bo… ▽ More We present simulation results of a gradual solar energetic particle (SEP) event detected on 2021 October 9 by multiple spacecraft, including BepiColombo (Bepi) and near-Earth spacecraft such as the Advanced Composition Explorer (ACE). A peculiarity of this event is that the presence of a high speed stream (HSS) affected the low-energy ion component ($\lesssim 5$ MeV) of the gradual SEP event at both Bepi and ACE, despite the HSS having only a modest solar wind speed increase. Using the EUHFORIA (European Heliospheric FORecasting Information Asset) magnetohydrodynamic model, we replicate the solar wind during the event and the coronal mass ejection (CME) that generated it. We then combine these results with the energetic particle transport model PARADISE (PArticle Radiation Asset Directed at Interplanetary Space Exploration). We find that the structure of the CME-driven shock was affected by the non-uniform solar wind, especially near the HSS, resulting in a shock wavefront with strong variations in its properties such as its compression ratio and obliquity. By scaling the emission of energetic particles from the shock to the solar wind compression at the shock, an excellent match between the PARADISE simulation and in-situ measurements of $\lesssim 5$ MeV ions is obtained. Our modelling shows that the intricate intensity variations observed at both ACE and Bepi were influenced by the non-uniform emission of energetic particles from the deformed shock wave and demonstrates the influence of even modest background solar wind structures on the development of SEP events. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal

On the onset delays of solar energetic electrons and protons: Evidence for a common accelerator

Authors: R. D. Strauss, N. Dresing, I. G. Richardson, J. P. van den Berg, P. J. Steyn

Abstract: The processes responsible for the acceleration of solar energetic particles (SEPs) are still not well understood, including whether SEP electrons and protons are accelerated by common or separate processes. Using a numerical particle transport model that includes both pitch-angle and perpendicular spatial diffusion, we simulate, amongst other quantities, the onset delay for MeV electrons and proto… ▽ More The processes responsible for the acceleration of solar energetic particles (SEPs) are still not well understood, including whether SEP electrons and protons are accelerated by common or separate processes. Using a numerical particle transport model that includes both pitch-angle and perpendicular spatial diffusion, we simulate, amongst other quantities, the onset delay for MeV electrons and protons and compare the results to observations of SEPs from widely-separated spacecraft. Such observations have previously been interpreted, in a simple scenario assuming no perpendicular diffusion, as evidence for different electron and proton sources. We show that, by assuming a common particle source together with perpendicular diffusion, we are able to simultaneously reproduce the onset delays for both electrons and protons. We argue that this points towards a common accelerator for these particles. Moreover, a relatively broad particle source is required in the model to correctly describe the observations. This is suggestive of diffusive shock acceleration occurring at large shock structures playing a significant role in the acceleration of these SEPs. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted to ApJ

An integrated online radioassay data storage and analytics tool for nEXO

Authors: R. H. M. Tsang, A. Piepke, S. Al Kharusi, E. Angelico, I. J. Arnquist, A. Atencio, I. Badhrees, J. Bane, V. Belov, E. P. Bernard, A. Bhat, T. Bhatta, A. Bolotnikov, P. A. Breur, J. P. Brodsky, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Q. Cao, D. Cesmecioglu, C. Chambers, E. Chambers, B. Chana, S. A. Charlebois , et al. (135 additional authors not shown)

Abstract: Large-scale low-background detectors are increasingly used in rare-event searches as experimental collaborations push for enhanced sensitivity. However, building such detectors, in practice, creates an abundance of radioassay data especially during the conceptual phase of an experiment when hundreds of materials are screened for radiopurity. A tool is needed to manage and make use of the radioassa… ▽ More Large-scale low-background detectors are increasingly used in rare-event searches as experimental collaborations push for enhanced sensitivity. However, building such detectors, in practice, creates an abundance of radioassay data especially during the conceptual phase of an experiment when hundreds of materials are screened for radiopurity. A tool is needed to manage and make use of the radioassay screening data to quantitatively assess detector design options. We have developed a Materials Database Application for the nEXO experiment to serve this purpose. This paper describes this database, explains how it functions, and discusses how it streamlines the design of the experiment. △ Less

Submitted 20 June, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Neutrinoless Double Beta Decay

Authors: C. Adams, K. Alfonso, C. Andreoiu, E. Angelico, I. J. Arnquist, J. A. A. Asaadi, F. T. Avignone, S. N. Axani, A. S. Barabash, P. S. Barbeau, L. Baudis, F. Bellini, M. Beretta, T. Bhatta, V. Biancacci, M. Biassoni, E. Bossio, P. A. Breur, J. P. Brodsky, C. Brofferio, E. Brown, R. Brugnera, T. Brunner, N. Burlac, E. Caden , et al. (207 additional authors not shown)

Abstract: This White Paper, prepared for the Fundamental Symmetries, Neutrons, and Neutrinos Town Meeting related to the 2023 Nuclear Physics Long Range Plan, makes the case for double beta decay as a critical component of the future nuclear physics program. The major experimental collaborations and many theorists have endorsed this white paper. This White Paper, prepared for the Fundamental Symmetries, Neutrons, and Neutrinos Town Meeting related to the 2023 Nuclear Physics Long Range Plan, makes the case for double beta decay as a critical component of the future nuclear physics program. The major experimental collaborations and many theorists have endorsed this white paper. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: white paper submitted for the Fundamental Symmetries, Neutrons, and Neutrinos Town Meeting in support of the US Nuclear Physics Long Range Planning Process

A non-parametric approach to the relation between the halo mass function and internal dark matter structure of haloes

Authors: T. R. G. Richardson, P. -S. Corasaniti

Abstract: Context. Galaxy cluster masses are usually defined as the mass within a spherical region enclosing a given matter overdensity (in units of the critical density). Converting masses from one overdensity definition to another can have several useful applications. Aims. In this article we present a generic non-parametric formalism that allows one to accurately map the halo mass function between two di… ▽ More Context. Galaxy cluster masses are usually defined as the mass within a spherical region enclosing a given matter overdensity (in units of the critical density). Converting masses from one overdensity definition to another can have several useful applications. Aims. In this article we present a generic non-parametric formalism that allows one to accurately map the halo mass function between two different mass overdensity definitions using the distribution of halo sparsities defined as the ratio of the two masses. We show that changing mass definitions reduces to modelling the distribution of halo sparsities. Methods. Using standard transformation rules of random variates, we derive relations between the halo mass function at different overdensities and the distribution of halo sparsities. Results. We show that these relations reproduce the N-body halo mass functions from the Uchuu simulation within the statistical errors at a few percent level. Furthermore, these relations allow the halo mass functions at different overdensities to be related to parametric descriptions of the halo density profile. In particular, we discuss the case of the concentration-mass relation of the Navarro-Frenk-White profile. Finally, we show that the use of such relations allows us to predict the distribution of sparsities of a sample of haloes of a given mass, thus opening the way to inferring cosmological constraints from individual galaxy cluster sparsity measurements. △ Less

Submitted 10 May, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: 13 pages, 9 figures, Accepted in A&A

Journal ref: A&A 674, A173 (2023)

Performance of novel VUV-sensitive Silicon Photo-Multipliers for nEXO

Authors: G. Gallina, Y. Guan, F. Retiere, G. Cao, A. Bolotnikov, I. Kotov, S. Rescia, A. K. Soma, T. Tsang, L. Darroch, T. Brunner, J. Bolster, J. R. Cohen, T. Pinto Franco, W. C. Gillis, H. Peltz Smalley, S. Thibado, A. Pocar, A. Bhat, A. Jamil, D. C. Moore, G. Adhikari, S. Al Kharusi, E. Angelico, I. J. Arnquist , et al. (140 additional authors not shown)

Abstract: Liquid xenon time projection chambers are promising detectors to search for neutrinoless double beta decay (0$νββ$), due to their response uniformity, monolithic sensitive volume, scalability to large target masses, and suitability for extremely low background operations. The nEXO collaboration has designed a tonne-scale time projection chamber that aims to search for 0$νββ$ of \ce{^{136}Xe} with… ▽ More Liquid xenon time projection chambers are promising detectors to search for neutrinoless double beta decay (0$νββ$), due to their response uniformity, monolithic sensitive volume, scalability to large target masses, and suitability for extremely low background operations. The nEXO collaboration has designed a tonne-scale time projection chamber that aims to search for 0$νββ$ of \ce{^{136}Xe} with projected half-life sensitivity of $1.35\times 10^{28}$~yr. To reach this sensitivity, the design goal for nEXO is $\leq$1\% energy resolution at the decay $Q$-value ($2458.07\pm 0.31$~keV). Reaching this resolution requires the efficient collection of both the ionization and scintillation produced in the detector. The nEXO design employs Silicon Photo-Multipliers (SiPMs) to detect the vacuum ultra-violet, 175 nm scintillation light of liquid xenon. This paper reports on the characterization of the newest vacuum ultra-violet sensitive Fondazione Bruno Kessler VUVHD3 SiPMs specifically designed for nEXO, as well as new measurements on new test samples of previously characterised Hamamatsu VUV4 Multi Pixel Photon Counters (MPPCs). Various SiPM and MPPC parameters, such as dark noise, gain, direct crosstalk, correlated avalanches and photon detection efficiency were measured as a function of the applied over voltage and wavelength at liquid xenon temperature (163~K). The results from this study are used to provide updated estimates of the achievable energy resolution at the decay $Q$-value for the nEXO design. △ Less

Submitted 25 November, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

CME Evolution in the Structured Heliosphere and Effects at Earth and Mars During Solar Minimum

Authors: Erika Palmerio, Christina O. Lee, Ian G. Richardson, Teresa Nieves-Chinchilla, Luiz F. G. Dos Santos, Jacob R. Gruesbeck, Nariaki V. Nitta, M. Leila Mays, Jasper S. Halekas, Cary Zeitlin, Shaosui Xu, Mats Holmström, Yoshifumi Futaana, Tamitha Mulligan, Benjamin J. Lynch, Janet G. Luhmann

Abstract: The activity of the Sun alternates between a solar minimum and a solar maximum, the former corresponding to a period of "quieter" status of the heliosphere. During solar minimum, it is in principle more straightforward to follow eruptive events and solar wind structures from their birth at the Sun throughout their interplanetary journey. In this paper, we report analysis of the origin, evolution,… ▽ More The activity of the Sun alternates between a solar minimum and a solar maximum, the former corresponding to a period of "quieter" status of the heliosphere. During solar minimum, it is in principle more straightforward to follow eruptive events and solar wind structures from their birth at the Sun throughout their interplanetary journey. In this paper, we report analysis of the origin, evolution, and heliospheric impact of a series of solar transient events that took place during the second half of August 2018, i.e. in the midst of the late declining phase of Solar Cycle 24. In particular, we focus on two successive coronal mass ejections (CMEs) and a following high-speed stream (HSS) on their way towards Earth and Mars. We find that the first CME impacted both planets, whilst the second caused a strong magnetic storm at Earth and went on to miss Mars, which nevertheless experienced space weather effects from the stream interacting region (SIR) preceding the HSS. Analysis of remote-sensing and in-situ data supported by heliospheric modelling suggests that CME--HSS interaction resulted in the second CME rotating and deflecting in interplanetary space, highlighting that accurately reproducing the ambient solar wind is crucial even during "simpler" solar minimum periods. Lastly, we discuss the upstream solar wind conditions and transient structures responsible for driving space weather effects at Earth and Mars. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 27 pages, 7 figures, 1 table, accepted for publication in Space Weather

The Influence of Visual Provenance Representations on Strategies in a Collaborative Hand-off Data Analysis Scenario

Authors: Jeremy E. Block, Shaghayegh Esmaeili, Eric D. Ragan, John R. Goodall, G. David Richardson

Abstract: Conducting data analysis tasks rarely occur in isolation. Especially in intelligence analysis scenarios where different experts contribute knowledge to a shared understanding, members must communicate how insights develop to establish common ground among collaborators. The use of provenance to communicate analytic sensemaking carries promise by describing the interactions and summarizing the steps… ▽ More Conducting data analysis tasks rarely occur in isolation. Especially in intelligence analysis scenarios where different experts contribute knowledge to a shared understanding, members must communicate how insights develop to establish common ground among collaborators. The use of provenance to communicate analytic sensemaking carries promise by describing the interactions and summarizing the steps taken to reach insights. Yet, no universal guidelines exist for communicating provenance in different settings. Our work focuses on the presentation of provenance information and the resulting conclusions reached and strategies used by new analysts. In an open-ended, 30-minute, textual exploration scenario, we qualitatively compare how adding different types of provenance information (specifically data coverage and interaction history) affects analysts' confidence in conclusions developed, propensity to repeat work, filtering of data, identification of relevant information, and typical investigation strategies. We see that data coverage (i.e., what was interacted with) provides provenance information without limiting individual investigation freedom. On the other hand, while interaction history (i.e., when something was interacted with) does not significantly encourage more mimicry, it does take more time to comfortably understand, as represented by less confident conclusions and less relevant information-gathering behaviors. Our results contribute empirical data towards understanding how provenance summarizations can influence analysis behaviors. △ Less

Submitted 7 August, 2022; originally announced August 2022.

Comments: to be published in IEEE Vis 2022

Journal ref: IEEE Transactions on Visualization and Computer Graphics 2022

Effect of boundary conditions on a high-performance isolation hexapod platform

Authors: Alessandro Stabile, Vladimir V. Yotov, Guglielmo S. Aglietti, Pasquale De Francesco, Guy Richardson

Abstract: Isolation of spacecraft microvibrations is essential for the successful deployment of instruments relying on high-precision pointing. Hexapod platforms represent a promising solution, but the difficulties associated with attaining desirable 3D dynamics within acceptable mass and complexity budgets have led to a minimal practical adoption. This paper addresses the influence of strut boundary condit… ▽ More Isolation of spacecraft microvibrations is essential for the successful deployment of instruments relying on high-precision pointing. Hexapod platforms represent a promising solution, but the difficulties associated with attaining desirable 3D dynamics within acceptable mass and complexity budgets have led to a minimal practical adoption. This paper addresses the influence of strut boundary conditions (BCs) on system-level mechanical disturbance suppression. Inherent limitations of the traditional all-rotational joint configuration are highlighted and shown to originate in link mass and rotational inertia. A pin-slider BC alternative is proposed and analytically proven to alleviate them in both 2D and 3D. The advantages of the new BC hold for arbitrary parallel manipulators and are demonstrated for several hexapod geometries through numerical tests. A configuration with favourable performance is suggested. Finally, a novel planar joint that allows the physical realisation of the proposed BC is described and validated. Consequently, this work enables the development of platforms for microvibration attenuation that do not require active control. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: Manuscript as accepted for publication in Mechanism and Machine Theory on 04.07.22; 23 pages, 17 figures

ACM Class: I.2.9

Journal ref: Mech. Mach. Theory, 177 (2022), 105020

Forecasting cosmological parameter constraints using multiple sparsity measurements as tracers of the mass profiles of dark matter haloes

Authors: P. S. Corasaniti, A. M. C. Le Brun, T. R. G. Richardson, Y. Rasera, S. Ettori, M. Arnaud, G. W Pratt

Abstract: The dark matter halo sparsity, i.e. the ratio between spherical halo masses enclosing two different overdensities, provides a non-parametric proxy of the halo mass distribution which has been shown to be a sensitive probe of the cosmological imprint encoded in the mass profile of haloes hosting galaxy clusters. Mass estimations at several overdensities would allow for multiple sparsity measurement… ▽ More The dark matter halo sparsity, i.e. the ratio between spherical halo masses enclosing two different overdensities, provides a non-parametric proxy of the halo mass distribution which has been shown to be a sensitive probe of the cosmological imprint encoded in the mass profile of haloes hosting galaxy clusters. Mass estimations at several overdensities would allow for multiple sparsity measurements, that can potentially retrieve the entirety of the cosmological information imprinted on the halo profile. Here, we investigate the impact of multiple sparsity measurements on the cosmological model parameter inference. For this purpose, we analyse N-body halo catalogues from the Raygal and M2Csims simulations and evaluate the correlations among six different sparsities from Spherical Overdensity halo masses at $Δ=200,500,1000$ and $2500$ (in units of the critical density). Remarkably, sparsities associated to distinct halo mass shells are not highly correlated. This is not the case for sparsities obtained using halo masses estimated from the Navarro-Frenk-White (NFW) best-fit profile, that artificially correlates different sparsities to order one. This implies that there is additional information in the mass profile beyond the NFW parametrization and that it can be exploited with multiple sparsities. In particular, from a likelihood analysis of synthetic average sparsity data, we show that cosmological parameter constraints significantly improve when increasing the number of sparsity combinations, though the constraints saturate beyond four sparsity estimates. We forecast constraints for the CHEX-MATE cluster sample and find that systematic mass bias errors mildly impact the parameter inference, though more studies are needed in this direction. △ Less

Submitted 1 August, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 16 pages, 11 figures, novel analyses including marginalisation over baryon systematics, updated to match MNRAS accepted version

Journal ref: MNRAS 516, 437-452 (2022)

CMEs and SEPs During November-December 2020: A Challenge for Real-Time Space Weather Forecasting

Authors: Erika Palmerio, Christina O. Lee, M. Leila Mays, Janet G. Luhmann, David Lario, Beatriz Sánchez-Cano, Ian G. Richardson, Rami Vainio, Michael L. Stevens, Christina M. S. Cohen, Konrad Steinvall, Christian Möstl, Andreas J. Weiss, Teresa Nieves-Chinchilla, Yan Li, Davin E. Larson, Daniel Heyner, Stuart D. Bale, Antoinette B. Galvin, Mats Holmström, Yuri V. Khotyaintsev, Milan Maksimovic, Igor G. Mitrofanov

Abstract: Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric sy… ▽ More Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric system, and multi-point measurements of solar transients can provide useful insights and validations for prediction models. In this work, we aim to analyse the whole inner heliospheric context between two eruptive flares that took place in late 2020, i.e. the M4.4 flare of November 29 and the C7.4 flare of December 7. This period is especially interesting because the STEREO-A spacecraft was located ~60° east of the Sun-Earth line, giving us the opportunity to test the capabilities of "predictions at 360°" using remote-sensing observations from the Lagrange L1 and L5 points as input. We simulate the CMEs that were ejected during our period of interest and the SEPs accelerated by their shocks using the WSA-Enlil-SEPMOD modelling chain and four sets of input parameters, forming a "mini-ensemble". We validate our results using in-situ observations at six locations, including Earth and Mars. We find that, despite some limitations arising from the models' architecture and assumptions, CMEs and shock-accelerated SEPs can be reasonably studied and forecast in real time at least out to several tens of degrees away from the eruption site using the prediction tools employed here. △ Less

Submitted 5 May, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 50 pages, 16 figures, 3 tables, accepted for publication in Space Weather (v2 contains updated CCMC links after their website update)

A Continuum of Physics-Based Lithium-Ion Battery Models Reviewed

Authors: Ferran Brosa Planella, Weilong Ai, Adam M. Boyce, Abir Ghosh, Ivan Korotkin, Smita Sahu, Valentin Sulzer, Robert Timms, Thomas G. Tranter, Maxim Zyskin, Samuel J. Cooper, Jacqueline S. Edge, Jamie M. Foster, Monica Marinescu, Billy Wu, Giles Richardson

Abstract: Physics-based electrochemical battery models derived from porous electrode theory are a very powerful tool for understanding lithium-ion batteries, as well as for improving their design and management. Different model fidelity, and thus model complexity, is needed for different applications. For example, in battery design we can afford longer computational times and the use of powerful computers,… ▽ More Physics-based electrochemical battery models derived from porous electrode theory are a very powerful tool for understanding lithium-ion batteries, as well as for improving their design and management. Different model fidelity, and thus model complexity, is needed for different applications. For example, in battery design we can afford longer computational times and the use of powerful computers, while for real-time battery control (e.g. in electric vehicles) we need to perform very fast calculations using simple devices. For this reason, simplified models that retain most of the features at a lower computational cost are widely used. Even though in the literature we often find these simplified models posed independently, leading to inconsistencies between models, they can actually be derived from more complicated models using a unified and systematic framework. In this review, we showcase this reductive framework, starting from a high-fidelity microscale model and reducing it all the way down to the Single Particle Model (SPM), deriving in the process other common models, such as the Doyle-Fuller-Newman (DFN) model. We also provide a critical discussion on the advantages and shortcomings of each of the models, which can aid model selection for a particular application. Finally, we provide an overview of possible extensions to the models, with a special focus on thermal models. Any of these extensions could be incorporated into the microscale model and the reductive framework re-applied to lead to a new generation of simplified, multi-physics models. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 50 pages, 9 figures

Interaction of coronal mass ejections and the solar wind. A force analysis

Authors: Dana-Camelia Talpeanu, Stefaan Poedts, Elke D'Huys, Marilena Mierla, Ian G. Richardson

Abstract: Aims. Our goal is to thoroughly analyse the dynamics of single and multiple solar eruptions, as well as a stealth ejecta. The data were obtained through self-consistent numerical simulations performed in a previous study. We also assess the effect of a different background solar wind on the propagation of these ejecta to Earth. Methods. We calculated all the components of the forces contributing… ▽ More Aims. Our goal is to thoroughly analyse the dynamics of single and multiple solar eruptions, as well as a stealth ejecta. The data were obtained through self-consistent numerical simulations performed in a previous study. We also assess the effect of a different background solar wind on the propagation of these ejecta to Earth. Methods. We calculated all the components of the forces contributing to the evolution of the numerically modelled consecutive coronal mass ejections (CMEs) obtained with the 2.5D magnetohydrodynamics (MHD) module of the code MPI-AMRVAC. We analysed the thermal and magnetic pressure gradients and the magnetic tension dictating the formation of several flux ropes in different locations in the aftermath of the eruptions. These three components were tracked in the equatorial plane during the propagation of the CMEs to Earth. Their interaction with other CMEs and with the background solar wind was also studied. Results. We explain the formation of the stealth ejecta and the plasma blobs (or plasmoids) occurring in the aftermath of solar eruptions. We also address the faster eruption of a CME in one case with a different background wind, even when the same triggering boundary motions were applied, and attribute this to the slightly different magnetic configuration and the large neighbouring arcade. The thermal pressure gradient revealed a shock in front of these slow eruptions, formed during their propagation to 1 AU. The double-peaked magnetic pressure gradient indicates that the triggering method affects the structure of the CMEs and that a part of the adjacent streamer is ejected along with the CME. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Journal ref: A&A 663, A32 (2022)

Development of a $^{127}$Xe calibration source for nEXO

Authors: B. G. Lenardo, C. A. Hardy, R. H. M. Tsang, J. C. Nzobadila Ondze, A. Piepke, S. Triambak, A. Jamil, G. Adhikari, S. Al Kharusi, E. Angelico, I. J. Arnquist, V. Belov, E. P. Bernard, A. Bhat, T. Bhatta, A. Bolotnikov, P. A. Breur, J. P. Brodsky, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Cao, B. Chana, S. A. Charlebois , et al. (103 additional authors not shown)

Abstract: We study a possible calibration technique for the nEXO experiment using a $^{127}$Xe electron capture source. nEXO is a next-generation search for neutrinoless double beta decay ($0νββ$) that will use a 5-tonne, monolithic liquid xenon time projection chamber (TPC). The xenon, used both as source and detection medium, will be enriched to 90% in $^{136}$Xe. To optimize the event reconstruction and… ▽ More We study a possible calibration technique for the nEXO experiment using a $^{127}$Xe electron capture source. nEXO is a next-generation search for neutrinoless double beta decay ($0νββ$) that will use a 5-tonne, monolithic liquid xenon time projection chamber (TPC). The xenon, used both as source and detection medium, will be enriched to 90% in $^{136}$Xe. To optimize the event reconstruction and energy resolution, calibrations are needed to map the position- and time-dependent detector response. The 36.3 day half-life of $^{127}$Xe and its small $Q$-value compared to that of $^{136}$Xe $0νββ$ would allow a small activity to be maintained continuously in the detector during normal operations without introducing additional backgrounds, thereby enabling in-situ calibration and monitoring of the detector response. In this work we describe a process for producing the source and preliminary experimental tests. We then use simulations to project the precision with which such a source could calibrate spatial corrections to the light and charge response of the nEXO TPC. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 24 pages, 16 figures

Timing the last major merger of galaxy clusters with large halo sparsity

Authors: T. R. G. Richardson, P. -S. Corasaniti

Abstract: Numerical simulations have shown that massive dark matter haloes, which today host galaxy clusters, assemble their mass over time alternating periods of quiescent accretion and phases of rapid growth associated with major merger episodes. Observations of such events in clusters can provide insights on the astrophysical processes that characterise the properties of the intra-cluster medium, as well… ▽ More Numerical simulations have shown that massive dark matter haloes, which today host galaxy clusters, assemble their mass over time alternating periods of quiescent accretion and phases of rapid growth associated with major merger episodes. Observations of such events in clusters can provide insights on the astrophysical processes that characterise the properties of the intra-cluster medium, as well as the gravitational processes that contribute to their assembly. It is therefore of prime interest to devise a fast and reliable way of detecting such perturbed systems. We present a novel approach to identifying and timing major mergers in clusters characterised by large values of halo sparsity. Using halo catalogues from the MultiDark-Planck2 simulation, we show that major merger events disrupt the radial mass distribution of haloes, thus leaving a distinct universal imprint on the evolution of halo sparsity over a period not exceeding two dynamical times. We exploit this feature using numerically calibrated distributions to test whether an observed galaxy cluster with given sparsity measurements has undergone a recent major merger and to eventually estimate when such an event occurred. We implement these statistical tools in a specifically developed public python library \textsc{lammas}, which we apply to the analysis of Abell 383 and Abell 2345 as test cases. Finding that, for example, Abell 2345 had a major merger about $2.1\pm 0.2$ Gyr ago. This work opens the way to detecting and timing major mergers in galaxy clusters solely through measurements of their mass at different radii. △ Less

Submitted 31 March, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 18 pages, 14 figures, submitted to MNRAS, comments welcome

Journal ref: MNRAS 513, 4951-4967 (2022)

Understanding the Origins of Problem Geomagnetic Storms Associated With "Stealth" Coronal Mass Ejections

Authors: Nariaki V. Nitta, Tamitha Mulligan, Emilia K. J. Kilpua, Benjamin J. Lynch, Marilena Mierla, Jennifer O'Kane, Paolo Pagano, Erika Palmerio, Jens Pomoell, Ian G. Richardson, Luciano Rodriguez, Alexis P. Rouillard, Suvadip Sinha, Nandita Srivastava, Dana-Camelia Talpeanu, Stephanie L. Yardley, Andrei N. Zhukov

Abstract: Geomagnetic storms are an important aspect of space weather and can result in significant impacts on space- and ground-based assets. The majority of strong storms are associated with the passage of interplanetary coronal mass ejections (ICMEs) in the near-Earth environment. In many cases, these ICMEs can be traced back unambiguously to a specific coronal mass ejection (CME) and solar activity on t… ▽ More Geomagnetic storms are an important aspect of space weather and can result in significant impacts on space- and ground-based assets. The majority of strong storms are associated with the passage of interplanetary coronal mass ejections (ICMEs) in the near-Earth environment. In many cases, these ICMEs can be traced back unambiguously to a specific coronal mass ejection (CME) and solar activity on the frontside of the Sun. Hence, predicting the arrival of ICMEs at Earth from routine observations of CMEs and solar activity currently makes a major contribution to the forecasting of geomagnetic storms. However, it is clear that some ICMEs, which may also cause enhanced geomagnetic activity, cannot be traced back to an observed CME, or, if the CME is identified, its origin may be elusive or ambiguous in coronal images. Such CMEs have been termed "stealth CMEs." In this review, we focus on these "problem" geomagnetic storms in the sense that the solar/CME precursors are enigmatic and stealthy. We start by reviewing evidence for stealth CMEs discussed in past studies. We then identify several moderate to strong geomagnetic storms (minimum Dst < -50 nT) in solar cycle 24 for which the related solar sources and/or CMEs are unclear and apparently stealthy. We discuss the solar and in situ circumstances of these events and identify several scenarios that may account for their elusive solar signatures. These range from observational limitations (e.g., a coronagraph near Earth may not detect an incoming CME if it is diffuse and not wide enough) to the possibility that there is a class of mass ejections from the Sun that have only weak or hard-to-observe coronal signatures. In particular, some of these sources are only clearly revealed by considering the evolution of coronal structures over longer time intervals than is usually considered. We also review a variety of numerical modelling approaches... △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 60 pages, 25 figures, 5 tables. Accepted for publication in Space Science Reviews

Journal ref: Space Sci Rev 217, 82 (2021)

A System-on-Chip for Closed-loop Optogenetic Sleep Modulation

Authors: Xilin Liu, Andrew G. Richardson

Abstract: Stimulation of target neuronal populations using optogenetic techniques during specific sleep stages has begun to elucidate the mechanisms and effects of sleep. To conduct closed-loop optogenetic sleep studies in untethered animals, we designed a fully integrated, low-power system-on-chip (SoC) for real-time sleep stage classification and stage-specific optical stimulation. The SoC consists of a 4… ▽ More Stimulation of target neuronal populations using optogenetic techniques during specific sleep stages has begun to elucidate the mechanisms and effects of sleep. To conduct closed-loop optogenetic sleep studies in untethered animals, we designed a fully integrated, low-power system-on-chip (SoC) for real-time sleep stage classification and stage-specific optical stimulation. The SoC consists of a 4-channel analog front-end for recording polysomnography signals, a mixed-signal machine-learning (ML) core, and a 16-channel optical stimulation back-end. A novel ML algorithm and innovative circuit design techniques improved the online classification performance while minimizing power consumption. The SoC was designed and simulated in 180 nm CMOS technology. In an evaluation using an expert labeled sleep database with 20 subjects, the SoC achieves a high sensitivity of 0.806 and a specificity of 0.947 in discriminating 5 sleep stages. Overall power consumption in continuous operation is 97 uW. △ Less

Submitted 16 July, 2021; originally announced August 2021.

Comments: 8 pages, 9 figures, 2 tables

Investigating Remote-sensing Techniques to Reveal Stealth Coronal Mass Ejections

Authors: Erika Palmerio, Nariaki V. Nitta, Tamitha Mulligan, Marilena Mierla, Jennifer O'Kane, Ian G. Richardson, Suvadip Sinha, Nandita Srivastava, Stephanie L. Yardley, Andrei N. Zhukov

Abstract: Eruptions of coronal mass ejections (CMEs) from the Sun are usually associated with a number of signatures that can be identified in solar disc imagery. However, there are cases in which a CME that is well observed in coronagraph data is missing a clear low-coronal counterpart. These events have received attention during recent years, mainly as a result of the increased availability of multi-point… ▽ More Eruptions of coronal mass ejections (CMEs) from the Sun are usually associated with a number of signatures that can be identified in solar disc imagery. However, there are cases in which a CME that is well observed in coronagraph data is missing a clear low-coronal counterpart. These events have received attention during recent years, mainly as a result of the increased availability of multi-point observations, and are now known as 'stealth CMEs'. In this work, we analyse examples of stealth CMEs featuring various levels of ambiguity. All the selected case studies produced a large-scale CME detected by coronagraphs and were observed from at least one secondary viewpoint, enabling a priori knowledge of their approximate source region. To each event, we apply several image processing and geometric techniques with the aim to evaluate whether such methods can provide additional information compared to the study of "normal" intensity images. We are able to identify at least weak eruptive signatures for all events upon careful investigation of remote-sensing data, noting that differently processed images may be needed to properly interpret and analyse elusive observations. We also find that the effectiveness of geometric techniques strongly depends on the CME propagation direction with respect to the observers and the relative spacecraft separation. Being able to observe and therefore forecast stealth CMEs is of great importance in the context of space weather, since such events are occasionally the solar counterparts of so-called 'problem geomagnetic storms'. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 26 pages, 8 figures, 1 table, accepted for publication in Frontiers in Astronomy and Space Sciences

A new ideality factor for perovskite solar cells and an analytical theory for their impedance spectroscopy response

Authors: Laurence J. Bennett, Antonio J. Riquelme, Nicola E. Courtier, Juan A. Anta, Giles Richardson

Abstract: Impedance spectroscopy (IS) is a relatively straightforward experimental technique that is commonly used to obtain information about the physical and chemical characteristics of photovoltaic devices. However, the non-standard physical behaviour of perovskite solar cells (PSC), which are heavily influenced by the motion of mobile ion vacancies, has hindered efforts to obtain a consistent theory to… ▽ More Impedance spectroscopy (IS) is a relatively straightforward experimental technique that is commonly used to obtain information about the physical and chemical characteristics of photovoltaic devices. However, the non-standard physical behaviour of perovskite solar cells (PSC), which are heavily influenced by the motion of mobile ion vacancies, has hindered efforts to obtain a consistent theory to interpret PSC impedance data. This work rectifies this omission by deriving a simple analytic model of the impedance response of a PSC from the underlying drift-diffusion model of charge carrier dynamics and ion vacancy motion. Extremely good agreement is shown between the analytic model and the much more complex drift-diffusion model in regimes (including maximum power point) where the applied voltage is close to the open circuit voltage $V_{oc}$. Both models show good qualitative agreement to experimental IS data in the literature and predict many of the observed anomalous features found in impedance measurements on PSCs, such as `the giant low frequency capacitance` and `inductive arcs' in the Nyquist plots. Where the physical properties of the PSC are already known the analytic model can be used to predict the recombination current $j_{rec}$ and the high and low frequency resistances and capacitances of the cell, $R_{HF}$, $C_{HF}$, $R_{LF}$ and $C_{LF}$. In scenarios where the physical properties of the cell are unknown the analytic model can also used to extract physical parameters from experimental PSC impedance data. {A novel physical parameter of particular significance to PSC physics is identified. This is termed the electronic ideality factor, $n_{el}$, and (as opposed to the standard ideality factor) can be used to deduce the dominant source of recombination in a PSC, independent of its ionic properties. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 64 pages, 24 figures

Heat Generation and a Conservation Law for Chemical Energy in Lithium-ion Batteries

Authors: G. Richardson, I. Korotkin

Abstract: Present theories of irreversible energy losses and heat generation within Li-ion cells are unsatisfactory because they are not compatible with energy conservation. This work aims to provide a consistent theoretical treatment of energy transport and losses in such devices. An energy conservation law is derived from the Doyle-Fuller-Newman (DFN) model of a Li-ion cell using a rigorous mathematical a… ▽ More Present theories of irreversible energy losses and heat generation within Li-ion cells are unsatisfactory because they are not compatible with energy conservation. This work aims to provide a consistent theoretical treatment of energy transport and losses in such devices. An energy conservation law is derived from the Doyle-Fuller-Newman (DFN) model of a Li-ion cell using a rigorous mathematical approach. The resulting law allows irreversible chemical energy losses to be located to seven different regions of the cell, namely: (i) the electrolyte, (ii) the anode particles, (iii) the cathode particles, (iv) the solid parts of the anode (ohmic losses), (v) the solid parts of the cathode (ohmic losses), (vi) the surfaces of the anode particles (polarisation losses), and (vii) the surfaces of the cathode particles (polarisation losses). Numerical solutions to the DFN model are used to validate the conservation law in the cases of a drive cycle and constant current discharges, and to compare the energy losses occurring in different locations. It is indicated how cell design can be improved, for a specified set of operating conditions, by comparing the magnitude of energy losses in the different regions of the cell. △ Less

Submitted 4 April, 2021; originally announced April 2021.

DandeLiion v1: An extremely fast solver for the Newman model of lithium-ion battery (dis)charge

Authors: Ivan Korotkin, Smita Sahu, Simon O'Kane, Giles Richardson, Jamie M. Foster

Abstract: DandeLiion (available at dandeliion.com) is a robust and extremely fast solver for the Doyle Fuller Newman (DFN) model, the standard electrochemical model for (dis)charge of a planar lithium-ion cell. DandeLiion conserves lithium, uses a second order spatial discretisation method (enabling accurate computations using relatively coarse discretisations) and is many times faster than its competitors.… ▽ More DandeLiion (available at dandeliion.com) is a robust and extremely fast solver for the Doyle Fuller Newman (DFN) model, the standard electrochemical model for (dis)charge of a planar lithium-ion cell. DandeLiion conserves lithium, uses a second order spatial discretisation method (enabling accurate computations using relatively coarse discretisations) and is many times faster than its competitors. The code can be used `in the cloud' and does not require installation before use. The difference in compute time between DandeLiion and its commercial counterparts is roughly a factor of 100 for the moderately-sized test case of the discharge of a single cell. Its linear scaling property means that the disparity in performance is even more pronounced for bigger systems, making it particularly suitable for applications involving multiple coupled cells. The model is characterised by a number of phenomenological parameters and functions, which may either be provided by the user or chosen from DandeLiion's library. This library contains data for the most commonly used electrolyte (LiPF6) and a number of common active material chemistries including graphite, lithium iron phosphate (LFP), nickel cobalt aluminium (NCA), and a variant of nickel cobalt manganese (NMC). △ Less

Submitted 8 February, 2021; originally announced February 2021.

Non-Halo Structures and their Effects on Gravitational Lensing

Authors: T. R. G. Richardson, J. Stücker, R. E. Angulo, O. Hahn

Abstract: Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_χ> 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be… ▽ More Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_χ> 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be correct as other structures, such as filaments and pancakes, exist and make up a significant fraction of the mass in the universe, ranging between 5$\%$ -- 50$\%$ depending on the dark matter model. Using novel fragmentation-free simulations of 1 and 3keV WDM cosmologies we study these "non-halo" structures and estimate their impact on flux-ratio observations. We find that these structures display sharp density gradients with short correlation lengths, and can contribute more to the lensing signal than all haloes up to the half-mode mass combined, thus reducing the differences expected among WDM models. We estimate that non-halo structures can be the dominant cause of line-of-sight flux-ratio anomalies in very warm, but already excluded, $m_x \sim 1 \rm{keV}$ scenarios. For colder cases $m_x \gtrsim 3 \rm{keV}$, we estimate that non-haloes can contribute about $5 - 10\%$ of the total flux-ratio signal. △ Less

Submitted 10 February, 2022; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: 14 pages, 10 figures. MNRAS submitted. Comments welcome

Empirical Model of 10-130 MeV Solar Energetic Particle Spectra at 1 AU Based on Coronal Mass Ejection Speed and Direction

Authors: Alessandro Bruno, Ian G. Richardson

Abstract: We present a new empirical model to predict solar energetic particle (SEP) event-integrated and peak intensity spectra between 10 and 130 MeV at 1 AU, based on multi-point spacecraft measurements from the Solar TErrestrial RElations Observatory (STEREO), the Geostationary Operational Environmental Satellites (GOES) and the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PA… ▽ More We present a new empirical model to predict solar energetic particle (SEP) event-integrated and peak intensity spectra between 10 and 130 MeV at 1 AU, based on multi-point spacecraft measurements from the Solar TErrestrial RElations Observatory (STEREO), the Geostationary Operational Environmental Satellites (GOES) and the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment. The analyzed data sample includes 32 SEP events occurring between 2010 and 2014, with a statistically significant proton signal at energies in excess of a few tens of MeV, unambiguously recorded at three spacecraft locations. The spatial distributions of SEP intensities are reconstructed by assuming an energy-dependent 2D Gaussian functional form, and accounting for the correlation between the intensity and the speed of the parent coronal mass ejection (CME), and the magnetic field line connection angle. The CME measurements used are from the Space Weather Database Of Notifications, Knowledge, Information (DONKI). The model performance, including its extrapolations to lower/higher energies, is tested by comparing with the spectra of 20 SEP events not used to derive the model parameters. Despite the simplicity of the model, the observed and predicted event-integrated and peak intensities at Earth and at the STEREO spacecraft for these events show remarkable agreement, both in the spectral shapes and their absolute values. △ Less

Submitted 23 February, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Journal ref: Sol Phys 296, 36 (2021)

Bayesian structural equation modeling for data from multiple cohorts

Authors: Khue-Dung Dang, Louise M. Ryan, Tugba Akkaya-Hocagil, Richard J. Cook, Gale A. Richardson, Nancy L. Day, Claire D. Coles, Heather Carmichael Olson, Sandra W. Jacobson, Joseph L. Jacobson

Abstract: While it is well known that high levels of prenatal alcohol exposure (PAE) result in significant cognitive deficits in children, the exact nature of the dose response is less well understood. In particular, there is a pressing need to identify the levels of PAE associated with an increased risk of clinically significant adverse effects. To address this issue, data have been combined from six longi… ▽ More While it is well known that high levels of prenatal alcohol exposure (PAE) result in significant cognitive deficits in children, the exact nature of the dose response is less well understood. In particular, there is a pressing need to identify the levels of PAE associated with an increased risk of clinically significant adverse effects. To address this issue, data have been combined from six longitudinal birth cohort studies in the United States that assessed the effects of PAE on cognitive outcomes measured from early school age through adolescence. Structural equation models (SEMs) are commonly used to capture the association among multiple observed outcomes in order to characterise the underlying variable of interest (in this case, cognition) and then relate it to PAE. However, it was not possible to apply classic SEM software in our context because different outcomes were measured in the six studies. In this paper we show how a Bayesian approach can be used to fit a multi-group multi-level structural model that maps cognition to a broad range of observed variables measured at multiple ages. These variables map to several different cognitive subdomains and are examined in relation to PAE after adjusting for confounding using propensity scores. The model also tests the possibility of a change point in the dose-response function. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Edge Deep Learning for Neural Implants

Authors: Xilin Liu, Andrew G. Richardson

Abstract: Implanted devices providing real-time neural activity classification and control are increasingly used to treat neurological disorders, such as epilepsy and Parkinson's disease. Classification performance is critical to identifying brain states appropriate for the therapeutic action. However, advanced algorithms that have shown promise in offline studies, in particular deep learning (DL) methods,… ▽ More Implanted devices providing real-time neural activity classification and control are increasingly used to treat neurological disorders, such as epilepsy and Parkinson's disease. Classification performance is critical to identifying brain states appropriate for the therapeutic action. However, advanced algorithms that have shown promise in offline studies, in particular deep learning (DL) methods, have not been deployed on resource-restrained neural implants. Here, we designed and optimized three embedded DL models of commonly adopted architectures and evaluated their inference performance in a case study of seizure detection. A deep neural network (DNN), a convolutional neural network (CNN), and a long short-term memory (LSTM) network were designed to classify ictal, preictal, and interictal phases from the CHB-MIT scalp EEG database. After iterative model compression and quantization, the algorithms were deployed on a general-purpose, off-the-shelf microcontroller. Inference sensitivity, false positive rate, execution time, memory size, and power consumption were quantified. For seizure event detection, the sensitivity and FPR (h-1) for the DNN, CNN, and LSTM models were 87.36%/0.169, 96.70%/0.102, and 97.61%/0.071, respectively. Predicting seizures for early warnings was also feasible. The implemented compression and quantization achieved a significant saving of power and memory with an accuracy degradation of less than 0.5%. Edge DL models achieved performance comparable to many prior implementations that had no time or computational resource limitations. Generic microcontrollers can provide the required memory and computational resources, while model designs can be migrated to ASICs for further optimization. The results suggest that edge DL inference is a feasible option for future neural implants to improve classification performance and therapeutic outcomes. △ Less

Submitted 1 June, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Comments: 25 pages, 7 figures, 3 tables

Journal ref: J Neural Eng. 2021 Apr 26;18(4)

A Fully Integrated Sensor-Brain-Machine Interface System for Restoring Somatosensation

Authors: Xilin Liu, Hongjie Zhu, Tian Qiu, Srihari Y. Sritharan, Dengteng Ge, Shu Yang, Milin Zhang, Andrew G. Richardson, Timothy H. Lucas, Nader Engheta, Jan Van der Spiegel

Abstract: Sensory feedback is critical to the performance of neural prostheses that restore movement control after neurological injury. Recent advances in direct neural control of paralyzed arms present new requirements for miniaturized, low-power sensor systems. To address this challenge, we developed a fully-integrated wireless sensor-brain-machine interface (SBMI) system for communicating key somatosenso… ▽ More Sensory feedback is critical to the performance of neural prostheses that restore movement control after neurological injury. Recent advances in direct neural control of paralyzed arms present new requirements for miniaturized, low-power sensor systems. To address this challenge, we developed a fully-integrated wireless sensor-brain-machine interface (SBMI) system for communicating key somatosensory signals, fingertip forces and limb joint angles, to the brain. The system consists of a tactile force sensor, an electrogoniometer, and a neural interface. The tactile force sensor features a novel optical waveguide on CMOS design for sensing. The electrogoniometer integrates an ultra low-power digital signal processor (DSP) for real-time joint angle measurement. The neural interface enables bidirectional neural stimulation and recording. Innovative designs of sensors and sensing interfaces, analog-to-digital converters (ADC) and ultra wide-band (UWB) wireless transceivers have been developed. The prototypes have been fabricated in 180nm standard CMOS technology and tested on the bench and in vivo. The developed system provides a novel solution for providing somatosensory feedback to next-generation neural prostheses. △ Less

Submitted 17 October, 2020; originally announced October 2020.

Comments: 12 pages, 17 figures

Journal ref: IEEE Sensors Journal, 2020

An Energy-efficient Wireless Neural Recording System with Compressed Sensing and Encryption

Authors: Xilin Liu, Andrew G. Richardson, Jan Van der Spiegel

Abstract: This paper presents a wireless neural recording system featuring energy-efficient data compression and encryption. An ultra-high efficiency is achieved by leveraging compressed sensing (CS) for simultaneous data compression and encryption. CS enables sub-Nyquist sampling of neural signals by taking advantage of its intrinsic sparsity. It simultaneously encrypts the data with the sampling matrix be… ▽ More This paper presents a wireless neural recording system featuring energy-efficient data compression and encryption. An ultra-high efficiency is achieved by leveraging compressed sensing (CS) for simultaneous data compression and encryption. CS enables sub-Nyquist sampling of neural signals by taking advantage of its intrinsic sparsity. It simultaneously encrypts the data with the sampling matrix being the cryptographic key. To share the key over an insecure wireless channel, we implement an elliptic-curve cryptography (ECC) based key exchanging protocol. The CS operation is executed in a custom-designed IC fabricated in 180nm CMOS technology. Mixed-signal circuits are designed to optimize the power efficiency of the matrix-vector multiplication (MVM) of the CS operation. The ECC algorithm is implemented in a low-power Cortex-M0 microcontroller (MCU). To be protected from timing and power analysis attacks, the implementation avoids possible data-dependent branches and also employs a randomized ECC initialization. At a compression ratio of 8x, the average correlated coefficient between the reconstructed signals and the uncompressed signals is 0.973, while the ciphertext-only attacks (CoA) achieve no better than 0.054 over 200,000 attacks. The prototype achieves a 35x power saving compared with conventional implementation in low-power MCUs. This work demonstrates a promising solution for future chronic neural recording systems with requirements in high energy efficiency and security. △ Less

Submitted 26 February, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 10 pages, 10 figures

A Hierarchical Meta-Analysis for Settings Involving Multiple Outcomes across Multiple Cohorts

Authors: Tugba Akkaya Hocagil, Louise M. Ryan, Richard J. Cook, Gale A. Richardson, Nancy L. Day, Claire D. Coles, Heather Carmichael Olson, Sandra W. Jacobson, Joseph L. Jacobson

Abstract: Evidence from animal models and epidemiological studies has linked prenatal alcohol exposure (PAE) to a broad range of long-term cognitive and behavioral deficits. However, there is virtually no information in the scientific literature regarding the levels of PAE associated with an increased risk of clinically significant adverse effects. During the period from 1975-1993, several prospective longi… ▽ More Evidence from animal models and epidemiological studies has linked prenatal alcohol exposure (PAE) to a broad range of long-term cognitive and behavioral deficits. However, there is virtually no information in the scientific literature regarding the levels of PAE associated with an increased risk of clinically significant adverse effects. During the period from 1975-1993, several prospective longitudinal cohort studies were conducted in the U.S., in which maternal reports regarding alcohol use were obtained during pregnancy and the cognitive development of the offspring was assessed from early childhood through early adulthood. The sample sizes in these cohorts did not provide sufficient power to examine effects associated with different levels and patterns of PAE. To address this critical public health issue, we have developed a hierarchical meta-analysis to synthesize information regarding the effects of PAE on cognition, integrating data on multiple endpoints from six U.S. longitudinal cohort studies. Our approach involves estimating the dose-response coefficients for each endpoint and then pooling these correlated dose-response coefficients to obtain an estimated `global' effect of exposure on cognition. In the first stage, we use individual participant data to derive estimates of the effects of PAE by fitting regression models that adjust for potential confounding variables using propensity scores. The correlation matrix characterizing the dependence between the endpoint-specific dose-response coefficients estimated within each cohort is then run, while accommodating incomplete information on some endpoints. We also compare and discuss inferences based on the proposed approach to inferences based on a full multivariate analysis △ Less

Submitted 2 September, 2020; originally announced September 2020.

Deducing the key physical properties of a perovskite solar cell from its impedance response: insights from drift-diffusion modelling

Authors: Antonio Riquelme, Laurence J. Bennett, Nicola E. Courtier, Matthew J. Wolf, Lidia Contreras-Bernal, Alison Walker, Giles Richardson, Juan A. Anta

Abstract: Interpreting the impedance response of perovskite solar cells (PSC) is significantly more challenging than for most other photovoltaics. This is for a variety of reasons, of which the most significant are the mixed ionic-electronic conduction properties of metal halide perovskites and the difficulty in fabricating stable, and reproducible, devices. Experimental studies, conducted on a variety of P… ▽ More Interpreting the impedance response of perovskite solar cells (PSC) is significantly more challenging than for most other photovoltaics. This is for a variety of reasons, of which the most significant are the mixed ionic-electronic conduction properties of metal halide perovskites and the difficulty in fabricating stable, and reproducible, devices. Experimental studies, conducted on a variety of PSCs, produce a variety of impedance spectra shapes. However, they all possess common features, the most noteworthy of which is that they have at least two signals, at high and low frequency, with different characteristic responses to temperature, illumination and electrical bias. It is shown, by a combination of experiment and drift-diffusion modelling of the ion and charge carrier transport and recombination within the cell, that these common features are well reproduced by the simulation. In addition, we show that the high frequency response contains all the key information relating to the steady-state performance of a PSC, i.e. it is a signature of the recombination mechanisms and provides a measure of charge collection efficiency. Moreover, steady-state performance is significantly affected by the distribution of mobile ionic charge within the perovskite layer. Comparison between the electrical properties of different devices should therefore be made using high frequency impedance measurements performed in the steady-state voltage regime in which the cell is expected to operate. △ Less

Submitted 16 March, 2020; originally announced March 2020.

Comments: 36 pages, 7 figures on the main text and 9 on the supporting information

Charge transport modelling of lithium ion batteries

Authors: Giles W Richardson, Jamie M Foster, Rahifa Ranom, Colin P Please, Angel M Ramos

Abstract: This paper presents the current state of mathematical modelling of the electrochemical behaviour of lithium-ion batteries as they are charged and discharged. It reviews the models developed by Newman and co-workers, both in the cases of dilute and moderately-concentrated electrolytes and indicates the modelling assumptions required for their development. Particular attention is paid to the interfa… ▽ More This paper presents the current state of mathematical modelling of the electrochemical behaviour of lithium-ion batteries as they are charged and discharged. It reviews the models developed by Newman and co-workers, both in the cases of dilute and moderately-concentrated electrolytes and indicates the modelling assumptions required for their development. Particular attention is paid to the interface conditions imposed between the electrolyte and the active electrode material; necessary conditions are derived for one of these, the Butler-Volmer relation, in order to ensure physically realistic solutions. Insight into the origin of the differences between various models found in the literature is revealed by considering formulations obtained by using different measures of the electric potential. Materials commonly used for electrodes in lithium ion batteries are considered and the various mathematical models used to describe lithium transport in them discussed. The problem of up-scaling from models of behaviour at the single electrode particle scale to the cell scale is addressed using homogenisation techniques resulting in the pseudo 2D model commonly used to describe charge transport and discharge behaviour in lithium-ion cells. Numerical solution to this model is discussed and illustrative results for a common device are computed. △ Less

Submitted 14 January, 2020; originally announced February 2020.

MSC Class: 35B27; 35Q99; 35M33

Direct comparison of high voltage breakdown measurements in liquid argon and liquid xenon

Authors: L. Tvrznikova, E. P. Bernard, S. Kravitz, K. O'Sullivan, G. Richardson, Q. Riffard, W. L. Waldron, J. Watson, D. N. McKinsey

Abstract: As noble liquid time projection chambers grow in size their high voltage requirements increase, and detailed, reproducible studies of dielectric breakdown and the onset of electroluminescence are needed to inform their design. The Xenon Breakdown Apparatus (XeBrA) is a 5-liter cryogenic chamber built to characterize the DC high voltage breakdown behavior of liquid xenon and liquid argon. Electrode… ▽ More As noble liquid time projection chambers grow in size their high voltage requirements increase, and detailed, reproducible studies of dielectric breakdown and the onset of electroluminescence are needed to inform their design. The Xenon Breakdown Apparatus (XeBrA) is a 5-liter cryogenic chamber built to characterize the DC high voltage breakdown behavior of liquid xenon and liquid argon. Electrodes with areas up to 33~cm$^2$ were tested while varying the cathode-anode separation from 1 to 6~mm with a voltage difference up to 75~kV. A power-law relationship between breakdown field and electrode area was observed. The breakdown behavior of liquid argon and liquid xenon within the same experimental apparatus was comparable. △ Less

Submitted 16 December, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Generalised single particle models for high-rate operation of graded lithium-ion electrodes: systematic derivation and validation

Authors: G. Richardson, I. Korotkin, R. Ranom. M. Castle, J. M. Foster

Abstract: A derivation of the single particle model (SPM) is made from a porous electrode theory model (or Newman model) of half-cell (dis)charge for an electrode composed of uniformly sized spherical electrode particles of a single chemistry. The derivation uses a formal asymptotic method based on the disparity between the size of the thermal voltage and that of the characteristic change in overpotential t… ▽ More A derivation of the single particle model (SPM) is made from a porous electrode theory model (or Newman model) of half-cell (dis)charge for an electrode composed of uniformly sized spherical electrode particles of a single chemistry. The derivation uses a formal asymptotic method based on the disparity between the size of the thermal voltage and that of the characteristic change in overpotential that occurs during (de)lithiation. Comparison is made between solutions to the SPM and to the porous electrode theory (PET) model for NMC, graphite and LFP. These are used to identify regimes where the SPM gives accurate predictions. For most chemistries, even at moderate (dis)charge rates, there are appreciable discrepancies between the PET model and the SPM which can be attributed to spatial non-uniformities in the electrolyte. This motivates us to calculate a correction term to the SPM. Once this has been incorporated into the model its accuracy is significantly improved. Generalised versions of the SPM, that can describe graded electrodes containing multiple electrode particle sizes (or chemistries), are also derived. The results of the generalised SPM, with the correction term, compare favourably to the full PET model where the active electrode material is either NMC or graphite. △ Less

Submitted 26 June, 2019; originally announced July 2019.

Comparing Long-Duration Gamma-Ray Flares and High-Energy Solar Energetic Particles

Authors: G. A. de Nolfo, A. Bruno, J. M. Ryan, S. Dalla, J. Giacalone, I. G. Richardson, E. R. Christian, S. J. Stochaj, G. A. Bazilevskaya, M. Boezio, M. Martucci, V. V. Mikhailov, R. Munini

Abstract: Little is known about the origin of the high-energy and sustained emission from solar Long-Duration Gamma-Ray Flares (LDGRFs), identified with the Compton Gamma Ray Observatory (CGRO), the Solar Maximum Mission (SMM), and now Fermi. Though Fermi/Large Area Space Telescope (LAT) has identified dozens of flares with LDGRF signature, the nature of this phenomenon has been a challenge to explain both… ▽ More Little is known about the origin of the high-energy and sustained emission from solar Long-Duration Gamma-Ray Flares (LDGRFs), identified with the Compton Gamma Ray Observatory (CGRO), the Solar Maximum Mission (SMM), and now Fermi. Though Fermi/Large Area Space Telescope (LAT) has identified dozens of flares with LDGRF signature, the nature of this phenomenon has been a challenge to explain both due to the extreme energies and long durations. The highest-energy emission has generally been attributed to pion production from the interaction of >300 MeV protons with the ambient matter. The extended duration suggests that particle acceleration occurs over large volumes extending high in the corona, either from stochastic acceleration within large coronal loops or from back precipitation from coronal mass ejection driven shocks. It is possible to test these models by making direct comparison between the properties of the accelerated ion population producing the gamma-ray emission derived from the Fermi/LAT observations, and the characteristics of solar energetic particles (SEPs) measured by the Payload for Matter-Antimatter Exploration and Light Nuclei Astrophysics (PAMELA) spacecraft in the energy range corresponding to the pion-related emission detected with Fermi. For fourteen of these events we compare the two populations -- SEPs in space and the interacting particles at the Sun -- and discuss the implications in terms of potential sources. Our analysis shows that the two proton numbers are poorly correlated, with their ratio spanning more than five orders of magnitude, suggesting that the back precipitation of shock-acceleration particles is unlikely the source of the LDGRF emission. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJ

Journal ref: de Nolfo et al. 2019, ApJ 879 90

Spectral Analysis of the September 2017 Solar Energetic Particle Events

Authors: A. Bruno, E. R. Christian, G. A. de Nolfo, I. G. Richardson, J. M. Ryan

Abstract: An interval of exceptional solar activity was registered in early September 2017, late in the decay phase of solar cycle 24, involving the complex Active Region 12673 as it rotated across the western hemisphere with respect to Earth. A large number of eruptions occurred between 4-10 September, including four associated with X-class flares. The X9.3 flare on 6 September and the X8.2 flare on 10 Sep… ▽ More An interval of exceptional solar activity was registered in early September 2017, late in the decay phase of solar cycle 24, involving the complex Active Region 12673 as it rotated across the western hemisphere with respect to Earth. A large number of eruptions occurred between 4-10 September, including four associated with X-class flares. The X9.3 flare on 6 September and the X8.2 flare on 10 September are currently the two largest during cycle 24. Both were accompanied by fast coronal mass ejections and gave rise to solar energetic particle (SEP) events measured by near-Earth spacecraft. In particular, the partially-occulted solar event on 10 September triggered a ground level enhancement (GLE), the second GLE of cycle 24. A further, much less energetic SEP event was recorded on 4 September. In this work we analyze observations by the Advanced Composition Explorer (ACE) and the Geostationary Operational Environmental Satellites (GOES), estimating the SEP event-integrated spectra above 300 keV and carrying out a detailed study of the spectral shape temporal evolution. Derived spectra are characterized by a low-energy break at few/tens of MeV; the 10 September event spectrum, extending up to ~1 GeV, exhibits an additional rollover at several hundred MeV. We discuss the spectral interpretation in the scenario of shock acceleration and in terms of other important external influences related to interplanetary transport and magnetic connectivity, taking advantage of multi-point observations from the Solar Terrestrial Relations Observatory (STEREO). Spectral results are also compared with those obtained for the 17 May 2012 GLE event. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Space Weather; 24 pages, 10 figures, 1 table

Journal ref: Space Weather 17(3) 419 (2019)

A Detailed Model of the Irish High Voltage Power Network for Simulating GICs

Authors: Seán P. Blake, Peter T. Gallagher, Joan Campanyà, Colin Hogg, Ciarán D. Beggan, Alan W. P. Thomson, Gemma S. Richardson, David Bell

Abstract: Constructing a power network model for geomagnetically induced current (GIC) calculations requires information on the DC resistances of elements within a network. This information is often not known, and power network models are simplified as a result, with assumptions used for network element resistances. Ireland's relatively small, isolated network presents an opportunity to model a complete pow… ▽ More Constructing a power network model for geomagnetically induced current (GIC) calculations requires information on the DC resistances of elements within a network. This information is often not known, and power network models are simplified as a result, with assumptions used for network element resistances. Ireland's relatively small, isolated network presents an opportunity to model a complete power network in detail, using as much real-world information as possible. A complete model of the Irish 400, 275, 220 and 110 kV network was made for GIC calculations, with detailed information on the number, type and DC resistances of transformers. The measured grounding resistances at a number of substations were also included in the model, which represents a considerable improvement on previous models of the Irish power network for GIC calculations. Sensitivity tests were performed to show how calculated GIC amplitudes are affected by different aspects of the model. These tests investigated: (1) How the orientation of a uniform electric field affects GICs. (2) The effect of including/omitting lower-voltage elements of the power network. (3) How the substation grounding resistances assumptions affected GIC values. It was found that changing the grounding resistance value had a considerable effect on calculated GICs at some substations, and no discernible effect at others. Finally, five recent geomagnetic storm events were simulated in the network. It was found that heavy rainfall prior to the 26-28 August 2015 geomagnetic storm event may have had a measurable impact on measured GIC amplitudes at a 400/220 kV transformer ground. △ Less

Submitted 9 November, 2018; originally announced November 2018.

Comments: 13 pages

Journal ref: Space Weather, 16 (2018)

Solar energetic particle events observed by the PAMELA mission

Authors: A. Bruno, G. A. Bazilevskaya, M. Boezio, E. R. Christian, G. A. de Nolfo, M. Martucci, M. Merge', V. V. Mikhailov, R. Munini, I. G. Richardson, J. M. Ryan, S. Stochaj, O. Adriani, G. C. Barbarino, R. Bellotti, E. A. Bogomolov, M. Bongi, V. Bonvicini, S. Bottai, F. Cafagna, D. Campana, P. Carlson, M. Casolino, G. Castellini, C. De Santis , et al. (33 additional authors not shown)

Abstract: Despite the significant progress achieved in recent years, the physical mechanisms underlying the origin of solar energetic particles (SEPs) are still a matter of debate. The complex nature of both particle acceleration and transport poses challenges to developing a universal picture of SEP events that encompasses both the low-energy (from tens of keV to a few hundreds of MeV) observations made by… ▽ More Despite the significant progress achieved in recent years, the physical mechanisms underlying the origin of solar energetic particles (SEPs) are still a matter of debate. The complex nature of both particle acceleration and transport poses challenges to developing a universal picture of SEP events that encompasses both the low-energy (from tens of keV to a few hundreds of MeV) observations made by space-based instruments and the GeV particles detected by the worldwide network of neutron monitors in ground-level enhancements (GLEs). The high-precision data collected by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment offer a unique opportunity to study the SEP fluxes between $\sim$80 MeV and a few GeV, significantly improving the characterization of the most energetic events. In particular, PAMELA can measure for the first time with good accuracy the spectral features at moderate and high energies, providing important constraints for current SEP models. In addition, the PAMELA observations allow the relationship between low and high-energy particles to be investigated, enabling a clearer view of the SEP origin. No qualitative distinction between the spectral shapes of GLE, sub-GLE and non-GLE events is observed, suggesting that GLEs are not a separate class, but are the subset of a continuous distribution of SEP events that are more intense at high energies. While the spectral forms found are to be consistent with diffusive shock acceleration theory, which predicts spectral rollovers at high energies that are attributed to particles escaping the shock region during acceleration, further work is required to explore the relative influences of acceleration and transport processes on SEP spectra. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: 26 pages, 13 figures, 2 tables

Journal ref: The Astrophysical Journal, 862:97, 2018 (17pp)

Evidence of energy and charge sign dependence of the recovery time for the December 2006 Forbush event measured by the PAMELA experiment

Authors: R. Munini, M. Boezio, A. Bruno, E. C. Christian, G. A. de Nolfo, V. Di Felice, M. Martucci, M. Merge, I. G. Richardson, J. M. Ryan, S. Stochaj, O. Adriani, G. C. Barbarino, G. A. Bazilevskaya, R. Bellotti, M. Bongi, V. Bonvicini, S. Bottai, F. Cafagna, D. Campana, P. Carlson, M. Casolino, G. Castellini, C. De Santis, A. M. Galper , et al. (33 additional authors not shown)

Abstract: New results on the short-term galactic cosmic ray (GCR) intensity variation (Forbush decrease) in December 2006 measured by the PAMELA instrument are presented. Forbush decreases are sudden suppressions of the GCR intensities which are associated with the passage of interplanetary transients such as shocks and interplanetary coronal mass ejections (ICMEs). Most of the past measurements of this phe… ▽ More New results on the short-term galactic cosmic ray (GCR) intensity variation (Forbush decrease) in December 2006 measured by the PAMELA instrument are presented. Forbush decreases are sudden suppressions of the GCR intensities which are associated with the passage of interplanetary transients such as shocks and interplanetary coronal mass ejections (ICMEs). Most of the past measurements of this phenomenon were carried out with ground-based detectors such as neutron monitors or muon telescopes. These techniques allow only the indirect detection of the overall GCR intensity over an integrated energy range. For the first time, thanks to the unique features of the PAMELA magnetic spectrometer, the Forbush decrease commencing on 2006 December 14, following a CME at the Sun on 2006 December 13 was studied in a wide rigidity range (0.4 - 20 GV) and for different species of GCRs detected directly in space. The daily averaged GCR proton intensity was used to investigate the rigidity dependence of the amplitude and the recovery time of the Forbush decrease. Additionally, for the first time, the temporal variations in the helium and electron intensities during a Forbush decrease were studied. Interestingly, the temporal evolutions of the helium and proton intensities during the Forbush decrease were found in good agreement, while the low rigidity electrons (< 2 GV) displayed a faster recovery. This difference in the electron recovery is interpreted as a charge-sign dependence introduced by drift motions experienced by the GCRs during their propagation through the heliosphere. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Journal ref: ApJ, 853, 76, 2018

A fast and robust numerical scheme for solving models of charge carrier transport and ion vacancy motion in perovskite solar cells

Authors: N. E. Courtier, G. Richardson, J. M. Foster

Abstract: Drift-diffusion models that account for the motion of both electronic and ionic charges are important tools for explaining the hysteretic behaviour and guiding the development of metal halide perovskite solar cells. Furnishing numerical solutions to such models for realistic operating conditions is challenging owing to the extreme values of some of the parameters. In particular, those characterisi… ▽ More Drift-diffusion models that account for the motion of both electronic and ionic charges are important tools for explaining the hysteretic behaviour and guiding the development of metal halide perovskite solar cells. Furnishing numerical solutions to such models for realistic operating conditions is challenging owing to the extreme values of some of the parameters. In particular, those characterising (i) the short Debye lengths (giving rise to rapid changes in the solutions across narrow layers), (ii) the relatively large potential differences across devices and (iii) the disparity in timescales between the motion of the electronic and ionic species give rise to significant stiffness. We present a finite difference scheme with an adaptive time step that is posed on a non-uniform staggered grid that provides second order accuracy in the mesh spacing. The method is able to cope with the stiffness of the system for realistic parameters values whilst providing high accuracy and maintaining modest computational costs. For example, a transient sweep of a current-voltage curve can be computed in only a few minutes on a standard desktop computer. △ Less

Submitted 17 January, 2018; originally announced January 2018.

Comments: 22 pages, 8 figures

Binder migration during drying of lithium-ion battery electrodes: modelling and comparison to experiment

Authors: Francesc Font, Bartosz Protas, Giles Richardson, Jamie M. Foster

Abstract: The drying process is a crucial step in electrode manufacture as it can affect the component distribution within the electrode. Phenomena such as binder migration can have negative effects in the form of poor cell performance (e.g. capacity fade) or mechanical failure (e.g. electrode delamination from the current collector). We present a mathematical model that tracks the evolution of the binder c… ▽ More The drying process is a crucial step in electrode manufacture as it can affect the component distribution within the electrode. Phenomena such as binder migration can have negative effects in the form of poor cell performance (e.g. capacity fade) or mechanical failure (e.g. electrode delamination from the current collector). We present a mathematical model that tracks the evolution of the binder concentration in the electrode during drying. Solutions to the model predict that low drying rates lead to a favourable homogeneous binder profile across the electrode film, whereas high drying rates result in an unfavourable accumulation of binder near the evaporation surface. These results show strong qualitative agreement with experimental observations and provide a cogent explanation for why fast drying conditions result in poorly performing electrodes. Finally, we provide some guidelines on how the drying process could be optimised to offer relatively short drying times whilst simultaneously maintaining a roughly homogeneous binder distribution. △ Less

Submitted 4 January, 2018; originally announced January 2018.

Asymptotic models for transport in large aspect ratio nanopores

Authors: B. Matejczyk, J. -F. Pietschmann, G. Richardson, M. -T. Wolfram

Abstract: Ion flow in charged nanopores is strongly influenced by the ratio of the Debye length to the pore radius. We investigate the asymptotic behaviour of solutions to the Poisson-Nernst-Planck (PNP) system in narrow pore like geometries and study the influence of the pore geometry and surface charge on ion transport. The physical properties of real pores motivate the investigation of distinguished asym… ▽ More Ion flow in charged nanopores is strongly influenced by the ratio of the Debye length to the pore radius. We investigate the asymptotic behaviour of solutions to the Poisson-Nernst-Planck (PNP) system in narrow pore like geometries and study the influence of the pore geometry and surface charge on ion transport. The physical properties of real pores motivate the investigation of distinguished asymptotic limits, in which either the Debye length and pore radius are comparable or the pore length is very much greater than its radius. This results in a Quasi-1D PNP model which can be further simplified, in the physically relevant limit of strong pore wall surface charge, to a fully one-dimensional model. Favourable comparison is made to the two-dimensional PNP equations in typical pore geometries. It is also shown that, for physically realistic parameters, the standard 1D Area Averaged PNP model for ion flow through a pore is a very poor approximation to the (real) two-dimensional solution to the PNP equations. This leads us to propose that the Quasi-1D PNP model derived here, whose computational cost is significantly less than two-dimensional solution of the PNP equations, should replace the use of the 1D Area Averaged PNP equations as a tool to investigate ion and current flows in ion pores. △ Less

Submitted 3 May, 2018; v1 submitted 2 January, 2018; originally announced January 2018.

Journal ref: Eur. J. Appl. Math 30 (2019) 557-584

Systematic derivation of a surface polarization model for planar perovskite solar cells

Authors: N. E. Courtier, J. M. Foster, S. E. J. O'Kane, A. B. Walker, G. Richardson

Abstract: Increasing evidence suggests that the presence of mobile ions in perovskite solar cells can cause a current-voltage curve hysteresis. Steady state and transient current-voltage characteristics of a planar metal halide CH$_3$NH$_3$PbI$_3$ perovskite solar cell are analysed with a drift-diffusion model that accounts for both charge transport and ion vacancy motion. The high ion vacancy density withi… ▽ More Increasing evidence suggests that the presence of mobile ions in perovskite solar cells can cause a current-voltage curve hysteresis. Steady state and transient current-voltage characteristics of a planar metal halide CH$_3$NH$_3$PbI$_3$ perovskite solar cell are analysed with a drift-diffusion model that accounts for both charge transport and ion vacancy motion. The high ion vacancy density within the perovskite layer gives rise to narrow Debye layers (typical width $\sim$2nm), adjacent to the interfaces with the transport layers, over which large drops in the electric potential occur and in which significant charge is stored. Large disparities between (I) the width of the Debye layers and that of the perovskite layer ($\sim$600nm) and (II) the ion vacancy density and the charge carrier densities motivate an asymptotic approach to solving the model, while the stiffness of the equations renders standard solution methods unreliable. We derive a simplified surface polarisation model in which the slow ion dynamic are replaced by interfacial (nonlinear) capacitances at the perovskite interfaces. Favourable comparison is made between the results of the asymptotic approach and numerical solutions for a realistic cell over a wide range of operating conditions of practical interest. △ Less

Submitted 30 August, 2017; originally announced August 2017.

Comments: 32 pages, 7 figures

Global Energetics of Solar Flares: V. Energy Closure in Flares and Coronal Mass Ejections

Authors: Markus J. Aschwanden, Amir Caspi, Christina M. S. Cohen, Gordon Holman, Ju Jing, Matthieu Kretzschmar, Eduard P. Kontar, James M. McTiernan, Richard A. Mewaldt, Aidan O'Flannagain, Ian G. Richardson, Daniel Ryan, Harry P. Warren, Yan Xu

Abstract: In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy o… ▽ More In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events. △ Less

Submitted 7 February, 2017; v1 submitted 4 January, 2017; originally announced January 2017.

Comments: 35 pages, 10 Figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Vol. 836, Issue 1, 17 (17pp); 2017 February 10