Showing 1–50 of 59 results for author: Davies, E

The importance of stochasticity in determining galaxy emissivities and UV LFs during cosmic dawn and reionization

Authors: Ivan Nikolić, Andrei Mesinger, James E. Davies, David Prelogović

Abstract: The stochastic nature of star formation and photon propagation in high-redshift galaxies can result in sizable galaxy-to-galaxy scatter in their properties. Ignoring this scatter by assuming mean quantities can bias estimates of their emissivity and corresponding observables. We construct a flexible, semi-empirical model, sampling scatter around the following mean relations: (i) the conditional ha… ▽ More The stochastic nature of star formation and photon propagation in high-redshift galaxies can result in sizable galaxy-to-galaxy scatter in their properties. Ignoring this scatter by assuming mean quantities can bias estimates of their emissivity and corresponding observables. We construct a flexible, semi-empirical model, sampling scatter around the following mean relations: (i) the conditional halo mass function (CHMF); (ii) the stellar-to-halo mass relation (SHMR); (iii) galaxy star formation main sequence (SFMS); (iv) fundamental metallicity relation (FMR); (v) conditional intrinsic luminosity; and (vi) photon escape fraction. In our fiducial model, ignoring scatter in these galaxy properties overestimates the duration of the EoR, delaying its completion by up to $Δz$ ~ 2. We quantify the relative importance of each of the above sources of scatter in determining the ionizing, soft-band X-ray and Lyman Werner (LW) emissivities as a function of scale and redshift. We find that scatter around the SFMS is important for all bands, especially at the highest redshifts where the emissivity is dominated by the faintest, most "bursty" galaxies. Ignoring this scatter would underestimate the mean emissivity and its standard deviation computed over 5 cMpc regions by factors of up to $\sim$2-10 at $5< z < 15$. Scatter around the X-ray luminosity to star formation rate relation is important for determining X-ray emissivity, accounting for roughly half of its mean and standard deviation. The importance of scatter in the ionizing escape fraction depends on its functional form, while scatter around the SHMR contributes at the level of ~10-20%. Although scatter does flatten the UV luminosity functions, shifting the bright end by 1-2 magnitudes, the level of scatter in our fiducial model is insufficient to fully explain recent estimates from JWST photometry (consistent with previous studies). △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 17 pages, 13 figures, submitted to A&A

Euclid. III. The NISP Instrument

Authors: Euclid Collaboration, K. Jahnke, W. Gillard, M. Schirmer, A. Ealet, T. Maciaszek, E. Prieto, R. Barbier, C. Bonoli, L. Corcione, S. Dusini, F. Grupp, F. Hormuth, S. Ligori, L. Martin, G. Morgante, C. Padilla, R. Toledo-Moreo, M. Trifoglio, L. Valenziano, R. Bender, F. J. Castander, B. Garilli, P. B. Lilje, H. -W. Rix , et al. (412 additional authors not shown)

Abstract: The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the proc… ▽ More The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Flux rope modeling of the 2022 Sep 5 CME observed by Parker Solar Probe and Solar Orbiter from 0.07 to 0.69 au

Authors: Emma E. Davies, Hannah T. Rüdisser, Ute V. Amerstorfer, Christian Möstl, Maike Bauer, Eva Weiler, Tanja Amerstorfer, Satabdwa Majumdar, Phillip Hess, Andreas J. Weiss, Martin A. Reiss, Lucie M. Green, David M. Long, Teresa Nieves-Chinchilla, Domenico Trotta, Timothy S. Horbury, Helen O'Brien, Edward Fauchon-Jones, Jean Morris, Christopher J. Owen, Stuart D. Bale, Justin C. Kasper

Abstract: As both Parker Solar Probe (PSP) and Solar Orbiter (SolO) reach heliocentric distances closer to the Sun, they present an exciting opportunity to study the structure of CMEs in the inner heliosphere. We present an analysis of the global flux rope structure of the 2022 September 5 CME event that impacted PSP at a heliocentric distance of only 0.07 au and SolO at 0.69 au. We compare in situ measurem… ▽ More As both Parker Solar Probe (PSP) and Solar Orbiter (SolO) reach heliocentric distances closer to the Sun, they present an exciting opportunity to study the structure of CMEs in the inner heliosphere. We present an analysis of the global flux rope structure of the 2022 September 5 CME event that impacted PSP at a heliocentric distance of only 0.07 au and SolO at 0.69 au. We compare in situ measurements at PSP and SolO to determine global and local expansion measures, finding a good agreement between magnetic field relationships with heliocentric distance, but significant differences with respect to flux rope size. We use PSP/WISPR images as input to the ELEvoHI model, providing a direct link between remote and in situ observations; we find a large discrepancy between the resulting modeled arrival times, suggesting that the underlying model assumptions may not be suitable when using data obtained close to the Sun, where the drag regime is markedly different in comparison to larger heliocentric distances. Finally, we fit the SolO/MAG and PSP/FIELDS data independently with the 3DCORE model and find that many parameters are consistent between spacecraft, however, challenges are apparent when reconstructing a global 3D structure that aligns with arrival times at PSP and Solar Orbiter, likely due to the large radial and longitudinal separations between spacecraft. From our model results, it is clear the solar wind background speed and drag regime strongly affects the modeled expansion and propagation of CMEs and needs to be taken into consideration. △ Less

Submitted 17 May, 2024; originally announced May 2024.

The ratio of [Eu/$α$] differentiates accreted/in-situ Milky Way stars across metallicities, as indicated by both field stars and globular clusters

Authors: Stephanie Monty, Vasily Belokurov, Jason L. Sanders, Terese T. Hansen, Charli M. Sakari, Madeleine McKenzie, GyuChul Myeong, Ellot Y. Davies, Anke Ardern-Arentsen, Davide Massari

Abstract: We combine stellar orbits with the abundances of the heavy, $r$-process element europium and the light, $α$-element, silicon to separate in-situ and accreted populations in the Milky Way across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in-situ, the levels of [Eu/Si] are lower.… ▽ More We combine stellar orbits with the abundances of the heavy, $r$-process element europium and the light, $α$-element, silicon to separate in-situ and accreted populations in the Milky Way across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in-situ, the levels of [Eu/Si] are lower. These systematically different levels of [Eu/Si] in the MW and the accreted halo imply that the scatter in [Eu/$α$] within a single galaxy is smaller than previously thought. At the lowest metallicities, we find that both accreted and in-situ populations trend down in [Eu/Si], consistent with enrichment via neutron star mergers. Through compiling a large dataset of abundances for 46 globular clusters (GCs), we show that differences in [Eu/Si] extend to populations of in-situ/accreted GCs. We interpret this consistency as evidence that in $r$-process elements, GCs trace the star formation history of their hosts, motivating their use as sub-Gyr timers of galactic evolution. Furthermore, fitting the trends in [Eu/Si] using a simple galactic chemical evolution model, we find that differences in [Eu/Si] between accreted and in-situ MW field stars cannot be explained through star formation efficiency alone. Finally, we show that the use of [Eu/Si] as a chemical tag between GCs and their host galaxies extends beyond the Local Group, to the halo of M31 - potentially offering the opportunity to do Galactic Archaeology in an external galaxy. △ Less

Submitted 24 May, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: 15 pages plus appendix. Submitted to MNRAS, comments welcome!

Understanding the effects of spacecraft trajectories through solar coronal mass ejection flux ropes using 3DCOREweb

Authors: Hannah Theresa Rüdisser, Andreas Jeffrey Weiss, Justin Le Louëdec, Ute V. Amerstorfer, Christian Möstl, Emma E. Davies, Helmut Lammer

Abstract: This study investigates the impact of spacecraft positioning and trajectory on in situ signatures of coronal mass ejections (CMEs). Employing the 3DCORE model, a 3D flux rope model that can generate in situ profiles for any given point in space and time, we conduct forward modeling to analyze such signatures for various latitudinal and longitudinal positions, with respect to the flux rope apex, at… ▽ More This study investigates the impact of spacecraft positioning and trajectory on in situ signatures of coronal mass ejections (CMEs). Employing the 3DCORE model, a 3D flux rope model that can generate in situ profiles for any given point in space and time, we conduct forward modeling to analyze such signatures for various latitudinal and longitudinal positions, with respect to the flux rope apex, at 0.8~au. Using this approach, we explore the appearance of the resulting in situ profiles for different flux rope types, with different handedness and inclination angles, for both high and low twist CMEs. Our findings reveal that CMEs exhibit distinct differences in signatures between apex hits and flank encounters, with the latter displaying elongated profiles with reduced rotation. However, constant, non-rotating in situ signatures are only observed for flank encounters of low twist CMEs, suggesting the existence of untwisted magnetic field lines within CME legs. Additionally, our study confirms the unambiguous appearance of different flux rope types in in situ signatures, contributing to the broader understanding and interpretation of observational data. Given the model assumptions, this may refute trajectory effects to be the cause for mismatching flux rope types as identified in solar signatures. While acknowledging limitations inherent in our model, such as the assumption of constant twist and non-deformable torus-like shape, we still draw relevant conclusions within the context of global magnetic field structures of CMEs and the potential for distinguishing flux rope types based on in situ observations. △ Less

Submitted 10 July, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: submitted to ApJ 2024 May 3, revised 2024 June 7, revised 2024 July 10, 23 pages, 21 figures, 2 tables

A Survey of Coronal Mass Ejections Measured In Situ by Parker Solar Probe During 2018-2022

Authors: Tarik M. Salman, Teresa Nieves-Chinchilla, Lan K. Jian, Noé Lugaz, Fernando Carcaboso, Emma E. Davies, Yaireska M. Collado-Vega

Abstract: We present a statistical investigation of the radial evolution of 28 interplanetary coronal mass ejections (ICMEs), measured in situ by the Parker Solar Probe (PSP) spacecraft from 2018 October to 2022 August. First, by analyzing the radial distribution of ICME classification based on magnetic hodograms, we find that coherent configurations are more likely to be observed close to the Sun. In contr… ▽ More We present a statistical investigation of the radial evolution of 28 interplanetary coronal mass ejections (ICMEs), measured in situ by the Parker Solar Probe (PSP) spacecraft from 2018 October to 2022 August. First, by analyzing the radial distribution of ICME classification based on magnetic hodograms, we find that coherent configurations are more likely to be observed close to the Sun. In contrast, more complex configurations are observed farther out. We also notice that the post-ICME magnetic field is more impacted following an ICME passage at larger heliocentric distances. Second, with a multi-linear robust regression, we derive a slower magnetic ejecta (ME) expansion rate within 1~au compared to previous statistical estimates. Then, investigating the magnetic field fluctuations within ICME sheaths, we see that these fluctuations are strongly coupled to the relative magnetic field strength gradient from the upstream solar wind to the ME. Third, we identify ME expansion as an important factor in forming sheaths. Finally, we determine the distortion parameter (DiP) which is a measure of magnetic field asymmetry in an ME. We discover lower overall asymmetries within MEs. We reveal that even for expanding MEs, the time duration over which an ME is sampled does not correlate with DiP values, indicating that the aging effect is not the sole contributor to the observed ME asymmetries. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Euclid preparation: TBD. The pre-launch Science Ground Segment simulation framework

Authors: Euclid Collaboration, S. Serrano, P. Hudelot, G. Seidel, J. E. Pollack, E. Jullo, F. Torradeflot, D. Benielli, R. Fahed, T. Auphan, J. Carretero, H. Aussel, P. Casenove, F. J. Castander, J. E. Davies, N. Fourmanoit, S. Huot, A. Kara, E. Keihanen, S. Kermiche, K. Okumura, J. Zoubian, A. Ealet, A. Boucaud, H. Bretonniere , et al. (251 additional authors not shown)

Abstract: The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previous… ▽ More The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previously as well as faster algorithms for the large-scale production of the expected Euclid data products. In this paper, we present the Euclid SGS simulation framework as applied in a large-scale end-to-end simulation exercise named Science Challenge 8. Our simulation pipeline enables the swift production of detailed image simulations for the construction and validation of the Euclid mission during its qualification phase and will serve as a reference throughout operations. Our end-to-end simulation framework starts with the production of a large cosmological N-body & mock galaxy catalogue simulation. We perform a selection of galaxies down to I_E=26 and 28 mag, respectively, for a Euclid Wide Survey spanning 165 deg^2 and a 1 deg^2 Euclid Deep Survey. We build realistic stellar density catalogues containing Milky Way-like stars down to H<26. Using the latest instrumental models for both the Euclid instruments and spacecraft as well as Euclid-like observing sequences, we emulate with high fidelity Euclid satellite imaging throughout the mission's lifetime. We present the SC8 data set consisting of overlapping visible and near-infrared Euclid Wide Survey and Euclid Deep Survey imaging and low-resolution spectroscopy along with ground-based. This extensive data set enables end-to-end testing of the entire ground segment data reduction and science analysis pipeline as well as the Euclid mission infrastructure, paving the way to future scientific and technical developments and enhancements. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 38 pages, 25 figures, A&A submitted

Hints of a disrupted binary dwarf galaxy in the Sagittarius stream

Authors: Elliot Y. Davies, Stephanie Monty, Vasily Belokurov, Adam M. Dillamore

Abstract: In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less-massive galaxy will often mostly be found in the extended stream and less-so in the central remnant. Motivated by these simulations, we use APOGEE DR17 chemical data from approximately 1100 stars… ▽ More In this work, we look for evidence of a non-unity mass ratio binary dwarf galaxy merger in the Sagittarius stream. Simulations of such a merger show that, upon merging with a host, particles from the less-massive galaxy will often mostly be found in the extended stream and less-so in the central remnant. Motivated by these simulations, we use APOGEE DR17 chemical data from approximately 1100 stars in both the Sagittarius remnant and stream to look for evidence of contamination from a second dwarf galaxy. This search is initially justified by the idea that disrupted binary dwarf galaxies provide a possible explanation of the Sagittarius bifurcation, and the location of the massive, chemically peculiar globular cluster NGC 2419 found within the stream of Sagittarius. We separate the Sagittarius data into its remnant and stream and compare the [Mg/Fe] content of the two populations. In particular, we select [Mg/Fe] to search for hints of unique star formation histories among our sample stars. Comparing the stream and remnant populations, we find regions have distinct [Mg/Fe] distributions for fixed [Fe/H], in addition to distinct chemical tracks in [Mg/Fe] -- [Fe/H] abundance space. We show that there are large regions of the tracks for which the probability of the two samples being drawn from the same distribution is very low (p < 0.05). Furthermore, we show that the two tracks can be fit with unique star formation histories using simple, one zone galactic chemical evolution models. While more work must be done to discern whether the hypothesis presented here is true, our work hints at the possibility that Sagittarius may consist of two dwarf galaxy progenitors. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 8 pages, 8 figures. Submitted to MNRAS. Comments welcome!

Properties of an interplanetary shock observed at 0.07 and 0.7 Astronomical Units by Parker Solar Probe and Solar Orbiter

Authors: D. Trotta, A. Larosa, G. Nicolaou, T. S. Horbury, L. Matteini, H. Hietala, X. Blanco-Cano, L. Franci, C. H. K. Chen, L. Zhao, G. P. Zank, C. M. S. Cohen, S. D. Bale, R. Laker, N. Fargette, F. Valentini, Y. Khotyaintsev, R. Kieokaew, N. Raouafi, E. Davies, R. Vainio, N. Dresing, E. Kilpua, T. Karlsson, C. J. Owen , et al. (1 additional authors not shown)

Abstract: The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with… ▽ More The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with the opportunity to study the shock properties at so different heliocentric distances. We characterize the shock, investigate its typical parameters and compare its small-scale features at both locations. Using the PSP observations, we investigate how magnetic switchbacks and ion cyclotron waves are processed upon shock crossing. We find that switchbacks preserve their V--B correlation while compressed upon the shock passage, and that the signature of ion cyclotron waves disappears downstream of the shock. By contrast, the SolO observations reveal a very structured shock transition, with a population of shock-accelerated protons of up to about 2 MeV, showing irregularities in the shock downstream, which we correlate with solar wind structures propagating across the shock. At SolO, we also report the presence of low-energy ($\sim$ 100 eV) electrons scattering due to upstream shocklets. This study elucidates how the local features of IP shocks and their environments can be very different as they propagate through the heliosphere. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: In review in ApJ

The Width of Magnetic Ejecta Measured Near 1 au: Lessons from STEREO-A Measurements in 2021--2022

Authors: Noé Lugaz, Bin Zhuang, Camilla Scolini, Nada Al-Haddad, Charles J. Farrugia, Réka M. Winslow, Florian Regnault, Christian Möstl, Emma E. Davies, Antoinette B. Galvin

Abstract: Coronal mass ejections (CMEs) are large-scale eruptions with a typical radial size at 1 au of 0.21 au but their angular width in interplanetary space is still mostly unknown, especially for the magnetic ejecta (ME) part of the CME. We take advantage of STEREO-A angular separation of 20$^\circ$-60$^\circ$ from the Sun-Earth line from October 2020 to August 2022, and perform a two-part study to cons… ▽ More Coronal mass ejections (CMEs) are large-scale eruptions with a typical radial size at 1 au of 0.21 au but their angular width in interplanetary space is still mostly unknown, especially for the magnetic ejecta (ME) part of the CME. We take advantage of STEREO-A angular separation of 20$^\circ$-60$^\circ$ from the Sun-Earth line from October 2020 to August 2022, and perform a two-part study to constrain the angular width of MEs in the ecliptic plane: a) we study all CMEs that are observed remotely to propagate between the Sun-STEREO-A and the Sun-Earth lines and determine how many impact one or both spacecraft in situ, and b) we investigate all in situ measurements at STEREO-A or at L1 of CMEs during the same time period to quantify how many are measured by the two spacecraft. A key finding is that, out of 21 CMEs propagating within 30$^\circ$ of either spacecraft, only four impacted both spacecraft and none provided clean magnetic cloud-like signatures at both spacecraft. Combining the two approaches, we conclude that the typical angular width of a ME at 1 au is $\sim$ 20$^\circ$-30$^\circ$, or 2-3 times less than often assumed and consistent with a 2:1 elliptical cross-section of an ellipsoidal ME. We discuss the consequences of this finding for future multi-spacecraft mission designs and for the coherence of CMEs. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 21 pages, revision submitted to ApJ

Discrepancies in the Properties of a Coronal Mass Ejection on Scales of 0.03~au as Revealed by Simultaneous Measurements at Solar Orbiter and Wind: The 2021 November 3--5 Event

Authors: F. Regnault, N. Al-Haddad, N. Lugaz, C. J. Farrugia, W. Yu, B. Zhuang, E. E. Davies

Abstract: Simultaneous in situ measurements of coronal mass ejections (CMEs), including both plasma and magnetic field, by two spacecraft in radial alignment have been extremely rare. Here, we report on one such CME measured by Solar Orbiter (SolO) and Wind on 2021 November 3--5, while the spacecraft were radially separated by a heliocentric distance of 0.13 au and angularly by only 2.2°. We focus on the ma… ▽ More Simultaneous in situ measurements of coronal mass ejections (CMEs), including both plasma and magnetic field, by two spacecraft in radial alignment have been extremely rare. Here, we report on one such CME measured by Solar Orbiter (SolO) and Wind on 2021 November 3--5, while the spacecraft were radially separated by a heliocentric distance of 0.13 au and angularly by only 2.2°. We focus on the magnetic cloud (MC) part of the CME. We find notable changes in the R and N magnetic field components and in the speed profiles inside the MC between SolO and Wind. We observe a greater speed at the spacecraft further away from the Sun without any clear compression signatures. Since spacecraft are close to each other and computing fast magnetosonic wave speed inside the MC we rule out temporal evolution as the reason on the observed differences suggesting that spatial variations over 2.2° of the MC structure are at the heart of the observed discrepancies. Moreover, using shock properties at SolO, we forecast an arrival time 2h30 too late for a shock that is just 5h31 away hours from Wind. Predicting the north-south component of the magnetic field at Wind from SolO measurements leads to a relative error of 55 %. These results show that even angular separations as low as 2.2° (or 0.03 au in arclength) between spacecraft can have a large impact on the observed CME properties, rising up the issue of the resolutions of current CME models and potentially affecting our forecasting capabilities. △ Less

Submitted 23 November, 2023; originally announced November 2023.

The effect of magnetic field line topology on ICME-related GCR modulation

Authors: Emma E. Davies, Camilla Scolini, Réka M. Winslow, Andrew P. Jordan, Christian Möstl

Abstract: The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to affect the galactic cosmic ray (GCR) flux measured in situ by spacecraft, causing temporary decreases known as Forbush decreases (Fds). In some ICMEs, the magnetic ejecta exhibits a magnetic flux rope (FR) structure; the strong magnetic field strength and closed field line geometry of such ICME FR… ▽ More The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to affect the galactic cosmic ray (GCR) flux measured in situ by spacecraft, causing temporary decreases known as Forbush decreases (Fds). In some ICMEs, the magnetic ejecta exhibits a magnetic flux rope (FR) structure; the strong magnetic field strength and closed field line geometry of such ICME FRs has been proposed to act as a shield to GCR transport. In this study, we identify four ICMEs near Earth that drove Fds with similar mean magnetic field strengths (20 - 25 nT); two ICMEs with more typical mean speeds (~400 km/s), and two fast (~750 km/s) ICMEs. Within each speed pairing, we identify an ICME that exhibited an open magnetic field line topology and compare its effect on the GCR flux to that which exhibited a mostly closed topology. We investigate the different mechanisms that contribute to the resulting ICME-related Fds and their recovery, and determine which properties, if any, play a more important role than others in driving Fds. We find that much of the GCR response to the ICME events in this study is independent of the open or closed magnetic field line topology of the flux rope, and that features such as the fluctuations in speed, magnetic field structure, and expansion within the FR may play more of a role in determining the smaller-scale structure of the Fd profile. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Investigating the Magnetic Structure of Interplanetary Coronal Mass Ejections using Simultaneous Multi-Spacecraft In situ Measurements

Authors: F. Regnault, N. Al-Haddad, N. Lugaz, C. J. Farrugia, W. Yu, E. E. Davies, A. B. Galvin, B. Zhuang

Abstract: In situ measurements from spacecraft typically provide a time series at a single location through coronal mass ejections (CMEs) and they have been one of the main methods to investigate CMEs. CME properties derived from these in situ measurements are affected by temporal changes that occur as the CME passes over the spacecraft, such as radial expansion and ageing, as well as spatial variations wit… ▽ More In situ measurements from spacecraft typically provide a time series at a single location through coronal mass ejections (CMEs) and they have been one of the main methods to investigate CMEs. CME properties derived from these in situ measurements are affected by temporal changes that occur as the CME passes over the spacecraft, such as radial expansion and ageing, as well as spatial variations within a CME. This study uses multi-spacecraft measurements of the same CME at close separations to investigate both the spatial variability (how different a CME profile is when probed by two spacecraft close to each other) and the so-called ageing effect (the effect of the time evolution on in situ properties). We compile a database of 19 events from the past four decades measured by two spacecraft with a radial separation <0.2 au and an angular separation <10°. We find that the average magnetic field strength measured by the two spacecraft differs by 18% of the typical average value, which highlights non-negligible spatial or temporal variations. For one particular event, measurements taken by the two spacecraft allow us to quantify and significantly reduce the ageing effect to estimate the asymmetry of the magnetic field strength profile. This study reveals that single-spacecraft time series near 1 au can be strongly affected by ageing and that correcting for self-similar expansion does not capture the whole ageing effect. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: Accepted for publication at ApJ. 18 pages, 6 figures, 2 tables

New Observations Needed to Advance Our Understanding of Coronal Mass Ejections

Authors: Erika Palmerio, Benjamin J. Lynch, Christina O. Lee, Lan K. Jian, Teresa Nieves-Chinchilla, Emma E. Davies, Brian E. Wood, Noé Lugaz, Réka M. Winslow, Tibor Török, Nada Al-Haddad, Florian Regnault, Meng Jin, Camilla Scolini, Fernando Carcaboso, Charles J. Farrugia, Vincent E. Ledvina, Cooper Downs, Christina Kay, Sanchita Pal, Tarik M. Salman, Robert C. Allen

Abstract: Coronal mass ejections (CMEs) are large eruptions from the Sun that propagate through the heliosphere after launch. Observational studies of these transient phenomena are usually based on 2D images of the Sun, corona, and heliosphere (remote-sensing data), as well as magnetic field, plasma, and particle samples along a 1D spacecraft trajectory (in-situ data). Given the large scales involved and th… ▽ More Coronal mass ejections (CMEs) are large eruptions from the Sun that propagate through the heliosphere after launch. Observational studies of these transient phenomena are usually based on 2D images of the Sun, corona, and heliosphere (remote-sensing data), as well as magnetic field, plasma, and particle samples along a 1D spacecraft trajectory (in-situ data). Given the large scales involved and the 3D nature of CMEs, such measurements are generally insufficient to build a comprehensive picture, especially in terms of local variations and overall geometry of the whole structure. This White Paper aims to address this issue by identifying the data sets and observational priorities that are needed to effectively advance our current understanding of the structure and evolution of CMEs, in both the remote-sensing and in-situ regimes. It also provides an outlook of possible missions and instruments that may yield significant improvements into the subject. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: White Paper submitted to the Heliophysics 2024-2033 Decadal Survey, 9 pages, 4 figures

Journal ref: Bulletin of the AAS, 55(3), 307, 2023

The eruption of a magnetic flux rope observed by \textit{Solar Orbiter} and \textit{Parker Solar Probe}

Authors: David M. Long, Lucie M. Green, Francesco Pecora, David H. Brooks, Hanna Strecker, David Orozco-Suárez, Laura A. Hayes, Emma E. Davies, Ute V. Amerstorfer, Marilena Mierla, David Lario, David Berghmans, Andrei N. Zhukov, Hannah T. Rüdisser

Abstract: Magnetic flux ropes are a key component of coronal mass ejections, forming the core of these eruptive phenomena. However, determining whether a flux rope is present prior to eruption onset and, if so, the rope's handedness and the number of turns that any helical field lines make is difficult without magnetic field modelling or in-situ detection of the flux rope. We present two distinct observatio… ▽ More Magnetic flux ropes are a key component of coronal mass ejections, forming the core of these eruptive phenomena. However, determining whether a flux rope is present prior to eruption onset and, if so, the rope's handedness and the number of turns that any helical field lines make is difficult without magnetic field modelling or in-situ detection of the flux rope. We present two distinct observations of plasma flows along a filament channel on 4 and 5 September 2022 made using the \textit{Solar Orbiter} spacecraft. Each plasma flow exhibited helical motions in a right-handed sense as the plasma moved from the source active region across the solar disk to the quiet Sun, suggesting that the magnetic configuration of the filament channel contains a flux rope with positive chirality and at least one turn. The length and velocity of the plasma flow increased from the first to the second observation, suggesting evolution of the flux rope, with the flux rope subsequently erupting within $\sim$5~hours of the second plasma flow. The erupting flux rope then passed over the \textit{Parker Solar Probe} spacecraft during its Encounter 13, enabling \textit{in-situ} diagnostics of the structure. Although complex and consistent with the flux rope erupting from underneath the heliospheric current sheet, the \textit{in-situ} measurements support the inference of a right-handed flux rope from remote-sensing observations. These observations provide a unique insight into the eruption and evolution of a magnetic flux rope near the Sun. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 19 pages, 11 figures, accepted for publication in The Astrophysical Journal

Disrupted dwarf binary merger as the possible origin of NGC 2419 and Sagittarius stream substructure

Authors: Elliot Y. Davies, Vasily Belokurov, Stephanie Monty, N. Wyn Evans

Abstract: Using $N$-body simulations, we demonstrate that satellite dwarf galaxy pairs which undergo significant mixing ($\sim 6$ Gyr) can have their respective most bound particles separated great distances upon subsequently merging with a more massive host. This may provide an explanation as to the origin of the complex globular cluster NGC 2149, which is found within the tail of the Sagittarius dwarf sph… ▽ More Using $N$-body simulations, we demonstrate that satellite dwarf galaxy pairs which undergo significant mixing ($\sim 6$ Gyr) can have their respective most bound particles separated great distances upon subsequently merging with a more massive host. This may provide an explanation as to the origin of the complex globular cluster NGC 2149, which is found within the tail of the Sagittarius dwarf spheroidal galaxy, yet separated from its central remnant by over 100 kpc. Dynamical investigations could support the chemical evidence which already points to the NGC 2419 being a nuclear star cluster. Motivated by the distinct nature of NGC 2419, we run a suite of simulations whereby an initial pre-infall merger of two satellites is followed by a post-infall merger of the remnant into a MW-like host potential. We present a striking example from our suite in this work, in which this separation is reproduced by the most bound particles of the two pre-infall satellites. Additionally, this double merger scenario can induce unusual on-sky features in the tidal debris of the post-infall merger, such as clouds, overdensities, and potentially new arms that could resemble the bifurcation observed in Sagittarius. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 5 pages, 4 figures, submitted to MNRAS. Comments welcome!

SimSpin v2.6.0 -- Constructing synthetic spectral IFU cubes for comparison with observational surveys

Authors: K. E. Harborne, A. Serene, E. J. A. Davies, C. Derkenne, S. Vaughan, A. I. Burdon, C. del P. Lagos, R. McDermid, S. O'Toole, C. Power, A. S. G. Robotham, G. Santucci, R. Tobar

Abstract: In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code,… ▽ More In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations. △ Less

Submitted 30 August, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 20 pages, 15 figures, 2 tables. Accepted for publication in PASA. 30/08/23

Stellar halo striations from assumptions of axisymmetry

Authors: Elliot Y. Davies, Adam M. Dillamore, Vasily Belokurov, N. Wyn Evans

Abstract: Motivated by the LMC's impact on the integral of motion space of the stellar halo, we run an $N$-body merger simulation to produce a population of halo-like stars. We subsequently move to a test particle simulation, in which the LMC perturbs this debris. When an axisymmetric potential is assumed for the final snapshot of the $N$-body merger remnant, a series of vertical striations in $(L_z, E)$ sp… ▽ More Motivated by the LMC's impact on the integral of motion space of the stellar halo, we run an $N$-body merger simulation to produce a population of halo-like stars. We subsequently move to a test particle simulation, in which the LMC perturbs this debris. When an axisymmetric potential is assumed for the final snapshot of the $N$-body merger remnant, a series of vertical striations in $(L_z, E)$ space form as the LMC approaches its pericentre. These result from the formation of overdensities in angular momentum owing to a relationship between the precession rate of near radial orbits and the torquing of these orbits by the LMC. This effect is heavily dependent on the shape of the inner potential. If a quadrupole component of the potential is included these striations become significantly less apparent due to the difference in precession rate between the two potentials. The absence of these features in data, and the dramatic change in orbital plane precession rate, discourages the use of an axisymmetric potential for highly eccentric orbits accreted from a massive GSE-like merger. Given the link between appearance of these striations and the shape of the potential, this effect may provide a new method of constraining the axisymmetry of the halo. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: 14 pages, 17 figures, submitted to MNRAS

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σ_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Stellar halo substructure generated by bar resonances

Authors: Adam M. Dillamore, Vasily Belokurov, N. Wyn Evans, Elliot Y. Davies

Abstract: Using data from the Gaia satellite's Radial Velocity Spectrometer Data Release 3 (RVS, DR3), we find a new and robust feature in the phase space distribution of halo stars. It is a prominent ridge at constant energy and with angular momentum $L_z>0$. We run test particle simulations of a stellar halo-like distribution of particles in a realistic Milky Way potential with a rotating bar. We observe… ▽ More Using data from the Gaia satellite's Radial Velocity Spectrometer Data Release 3 (RVS, DR3), we find a new and robust feature in the phase space distribution of halo stars. It is a prominent ridge at constant energy and with angular momentum $L_z>0$. We run test particle simulations of a stellar halo-like distribution of particles in a realistic Milky Way potential with a rotating bar. We observe similar structures generated in the simulations from the trapping of particles in resonances with the bar, particularly at the corotation resonance. Many of the orbits trapped at the resonances are halo-like, with large vertical excursions from the disc. The location of the observed structure in energy space is consistent with a bar pattern speed in the range $Ω_\mathrm{b}\approx35-40$ km s$^{-1}$ kpc$^{-1}$. Overall, the effect of the resonances is to give the inner stellar halo a mild, net spin in the direction of the bar's rotation. As the distribution of the angular momentum becomes asymmetric, a population of stars with positive mean $L_z$ and low vertical action is created. The variation of the average rotational velocity of the simulated stellar halo with radius is similar to the behaviour of metal-poor stars in data from the APOGEE survey. Though the effects of bar resonances have long been known in the Galactic disc, this is strong evidence that the bar can drive changes even in the diffuse and extended stellar halo through its resonances. △ Less

Submitted 20 July, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 13 pages, 15 figures, published in MNRAS

Accelerated phase-mixing in the stellar halo due to a rotating bar

Authors: Elliot Y. Davies, Adam M. Dillamore, Eugene Vasiliev, Vasily Belokurov

Abstract: In a galaxy merger, the stars tidally stripped from the satellite and accreted onto the host galaxy undergo phase mixing and form finely-grained structures in the phase space. However, these fragile structures may be destroyed in the subsequent galaxy evolution, in particular, by a rotating bar that appears well after the merger is completed. In this work, we investigate the survivability of phase… ▽ More In a galaxy merger, the stars tidally stripped from the satellite and accreted onto the host galaxy undergo phase mixing and form finely-grained structures in the phase space. However, these fragile structures may be destroyed in the subsequent galaxy evolution, in particular, by a rotating bar that appears well after the merger is completed. In this work, we investigate the survivability of phase-space structures in the presence of a bar. We find that a bar with amplitude and pattern speed similar to those of the Milky Way would blur and destroy a substantial amount of the substructure that consists of particles with pericentre radii comparable to the bar length. While this appears to be in tension with the recent discovery of phase-space chevrons in \textit{Gaia} DR3 data, the most prominent chevrons in our simulations can still be recovered when applying the same analysis procedure as in observations. Moreover, the smoothing effect is less pronounced in the population of stars whose angular momenta have the opposite sign to the bar pattern speed. △ Less

Submitted 3 March, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 5 pages, 7 figures, accepted by MNRAS

Characterizing ICME-related Forbush Decreases at Mercury using MESSENGER Observations: Identification of a One or Two-Step Structure

Authors: Emma E. Davies, Réka M. Winslow, David J. Lawrence

Abstract: The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to cause decreases in the galactic cosmic ray (GCR) flux measured in situ by spacecraft, known as Forbush decreases (Fds). We use measurements of the GCR count rate obtained by MESSENGER during its orbital phase around Mercury to identify such Fds related to the passage of ICMEs and characterize thei… ▽ More The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to cause decreases in the galactic cosmic ray (GCR) flux measured in situ by spacecraft, known as Forbush decreases (Fds). We use measurements of the GCR count rate obtained by MESSENGER during its orbital phase around Mercury to identify such Fds related to the passage of ICMEs and characterize their structure. Of the 42 ICMEs with corresponding high-quality GCR data, 79% are associated with a Fd. Thus a total of 33 ICME-related Fds were identified, 24 of which (73%) have a two-step structure. We use a superposed epoch analysis to build an average Fd profile at MESSENGER and find that despite the variability of individual events, a two-step structure is produced and is directly linked with the magnetic boundaries of the ICME. By using results from previous studies at Earth and Mars, we also address whether two-step Fds are more commonly observed closer to the Sun; we found that although likely, this is not conclusive when comparing to the wide range of results of previous studies conducted at Earth. Finally, we find that the percentage decrease in GCR flux of the Fd is greater at MESSENGER on average than at Earth and Mars, decreasing with increasing heliocentric distance. The relationship between the percentage decrease and maximum hourly decrease is also in agreement with previous studies, and follows trends relating to the expansion of ICMEs as they propagate through the heliosphere. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Successive interacting coronal mass ejections: How to create a perfect storm?

Authors: Gordon J. Koehn, Ravindra T. Desai, Emma E. Davies, Robert J. Forsyth, Jonathan P. Eastwood, Stefaan Poedts

Abstract: Coronal mass ejections (CMEs) are the largest type of eruptions on the Sun and the main driver of severe space weather at the Earth. In this study, we implement a force-free spheromak CME description within 3-D magnetohydrodynamic simulations to parametrically evaluate successive interacting CMEs within a representative heliosphere. We explore CME-CME interactions for a range of orientations, laun… ▽ More Coronal mass ejections (CMEs) are the largest type of eruptions on the Sun and the main driver of severe space weather at the Earth. In this study, we implement a force-free spheromak CME description within 3-D magnetohydrodynamic simulations to parametrically evaluate successive interacting CMEs within a representative heliosphere. We explore CME-CME interactions for a range of orientations, launch time variations and CME handedness and quantify their geo-effectiveness via the primary solar wind variables and empirical measures of the disturbance storm time index and subsolar magnetopause standoff distance. We show how the interaction of two moderate CMEs between the Sun and the Earth can translate into extreme conditions at the Earth and how CME-CME interactions at different radial distances can maximise different solar wind variables that induce different geophysical impacts. In particular, we demonstrate how the orientation and handedness of a given CME can have a significant impact on the conservation and loss of magnetic flux, and consequently B$_z$, due to magnetic reconnection with the interplanetary magnetic field. This study thus implicates identification of CME chirality in the solar corona as an early diagnostic for forecasting geomagnetic storms involving multiple CMEs. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 18 pages, 11 figures, 2 tables. The Astrophysical Journal, accepted 10 November 2022

Connecting Solar and Stellar Flares/CMEs: Expanding Heliophysics to Encompass Exoplanetary Space Weather

Authors: B. J. Lynch, B. E. Wood, M. Jin, T. Török, X. Sun, E. Palmerio, R. A. Osten, A. A. Vidotto, O. Cohen, J. D. Alvarado-Gómez, J. J. Drake, V. S. Airapetian, Y. Notsu, A. Veronig, K. Namekata, R. M. Winslow, L. K. Jian, A. Vourlidas, N. Lugaz, N. Al-Haddad, W. B. Manchester, C. Scolini, C. J. Farrugia, E. E. Davies, T. Nieves-Chinchilla , et al. (3 additional authors not shown)

Abstract: The aim of this white paper is to briefly summarize some of the outstanding gaps in the observations and modeling of stellar flares, CMEs, and exoplanetary space weather, and to discuss how the theoretical and computational tools and methods that have been developed in heliophysics can play a critical role in meeting these challenges. The maturity of data-inspired and data-constrained modeling of… ▽ More The aim of this white paper is to briefly summarize some of the outstanding gaps in the observations and modeling of stellar flares, CMEs, and exoplanetary space weather, and to discuss how the theoretical and computational tools and methods that have been developed in heliophysics can play a critical role in meeting these challenges. The maturity of data-inspired and data-constrained modeling of the Sun-to-Earth space weather chain provides a natural starting point for the development of new, multidisciplinary research and applications to other stars and their exoplanetary systems. Here we present recommendations for future solar CME research to further advance stellar flare and CME studies. These recommendations will require institutional and funding agency support for both fundamental research (e.g. theoretical considerations and idealized eruptive flare/CME numerical modeling) and applied research (e.g. data inspired/constrained modeling and estimating exoplanetary space weather impacts). In short, we recommend continued and expanded support for: (1.) Theoretical and numerical studies of CME initiation and low coronal evolution, including confinement of "failed" eruptions; (2.) Systematic analyses of Sun-as-a-star observations to develop and improve stellar CME detection techniques and alternatives; (3.) Improvements in data-inspired and data-constrained MHD modeling of solar CMEs and their application to stellar systems; and (4.) Encouraging comprehensive solar--stellar research collaborations and conferences through new interdisciplinary and multi-agency/division funding mechanisms. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 9 pages, 5 figures, white paper submitted to the Heliophysics 2024--2033 Decadal Survey

Journal ref: Bulletin of the American Astronomical Society, 55(3), 254, 2023

Ironing the folds: The phase space chevrons of a GSE-like merger as a dark matter subhalo detector

Authors: Elliot Y. Davies, Eugene Vasiliev, Vasily Belokurov, N. Wyn Evans, Adam M. Dillamore

Abstract: Recent work uncovered features in the phase space of the Milky Way's stellar halo which may be attributed to the last major merger. When stellar material from a satellite is accreted onto its host, it phase mixes and appears finely substructured in phase space. For a high-eccentricity merger, this substructure most clearly manifests as numerous wrapping chevrons in $(v_r, r)$ space, corresponding… ▽ More Recent work uncovered features in the phase space of the Milky Way's stellar halo which may be attributed to the last major merger. When stellar material from a satellite is accreted onto its host, it phase mixes and appears finely substructured in phase space. For a high-eccentricity merger, this substructure most clearly manifests as numerous wrapping chevrons in $(v_r, r)$ space, corresponding to stripes in $(E, θ_r)$ space. We introduce the idea of using this substructure as an alternative subhalo detector to cold stellar streams. We simulate an N-body merger akin to the GSE and assess the impact of subhaloes on these chevrons. We examine how their deformation depends on the mass, pericentre, and number of subhaloes. To quantify the impact of perturbers, we utilise the appearance of chevrons in $(E, θ_r)$ space to introduce a new quantity -- the ironing parameter. We show that: (1) a single flyby of a massive ($\sim 10^{10}$ M$_{\odot}$) subhalo with pericentre comparable to, or within, the shell's apocentre smooths out the substructure, (2) a single flyby of a low mass ($\lesssim 10^8$ M$_{\odot}$) has negligible effect, (3) multiple flybys of subhalos derived from a subhalo mass function between $10^7-10^{10}$ M$_{\odot}$ cause significant damage if deep within the potential, (4) the effects of known perturbers (e.g. Sagittarius) should be detectable and offer constraints on their initial mass. The sensitivity to the populations of subhaloes suggests that we should be able to place an upper limit on the Milky Way's subhalo mass function. △ Less

Submitted 5 December, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 18 pages, 20 figures, accepted by MNRAS

Energy wrinkles and phase-space folds of the last major merger

Authors: Vasily Belokurov, Eugene Vasiliev, Alis J. Deason, Sergey E. Koposov, Azadeh Fattahi, Adam M. Dillamore, Elliot Y. Davies, Robert J. J. Grand

Abstract: Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-s… ▽ More Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures are predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Additionally, the observed energy and angular momentum (E, L_z) distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E -- which, we surmise, may be evidence of the dwarf's rapid sinking -- and riddled with wrinkles and bumps. If these small-scale energy inhomogeneities have been seeded during or immediately after the interaction with the Milky Way, and are not due to the spatial restriction of our study, then making use of the (v_r,r) chevrons to constrain the time of the merger becomes cumbersome. Nonetheless, we demonstrate that similar phase-space and (E,L_z) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor's disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy's most important interaction, but also can be used as a novel tool to map out the Milky Way's current gravitational potential and its perturbations. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 16 pages, 14 figures, submitted to MNRAS. Comments are welcome!

The Transport of Extrusive Volcanic Deposits on Jezero Crater Through Paleofluvial Processes

Authors: Antonio Paris, Kate Morgan, Evan Davies

Abstract: Jezero, an impact crater in the Syrtis Major quadrangle of Mars, is generally thought to have amassed a large body of liquid water in its ancient past. NASA spectra of the proposed paleolake interpret the youngest surface unit as olivine-bearing minerals crystallized from magma. In early 2021, the Perseverance rover landed at the leading edge of a fan-delta deposit northwest of Jezero, an area arg… ▽ More Jezero, an impact crater in the Syrtis Major quadrangle of Mars, is generally thought to have amassed a large body of liquid water in its ancient past. NASA spectra of the proposed paleolake interpret the youngest surface unit as olivine-bearing minerals crystallized from magma. In early 2021, the Perseverance rover landed at the leading edge of a fan-delta deposit northwest of Jezero, an area argued to have experienced two distinct periods of fluvial activity. Surface imagery obtained by Perseverance reveal partially buried and unburied vesicular and non-vesicular rocks that appear volcanic in origin, emplaced sometime during the Noachian-Hesperian boundary. The absence of volcanic extrusive features along the fan-delta deposit, however, have made the origin of these ballast deposits a matter of contention among planetary scientists. To establish the origin of these basalt-like rocks, a comparison was made between analogous deposits on the Moenkopi Plateau in Arizona with similar deposits imaged by Perseverance on Jezero. The search for geologic analogs along the Moenkopi Plateau were guided by observable similarities in surface geomorphology, influenced and modified by fluvial, eolian, and past volcanic activity, primarily from the Late Pleistocene-Holocene boundary. By analyzing surface imagery taken by Perseverance and comparing it with the analogue site, we hypothesize that the exposed vesicular rocks imaged by Perseverance were likely transported into the paleolake by geomorphic interactions, specifically fluvial processes similarly to the deposits that were transported along drainage patterns we observed on the Moenkopi Plateau. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 12 pages, 12 figures

Journal ref: Journal of the Washington Academy of Sciences, Volume 108, Summer 2022

A Coronal Mass Ejection and Magnetic Ejecta Observed In Situ by STEREO-A and Wind at 55$^\circ$ Angular Separation

Authors: Noé Lugaz, Tarik M. Salman, Charles J. Farrugia, Wenyuan Yu, Bin Zhuang, Nada Al-Haddad, Camilla Scolini, Réka M. Winslow, Christian Möstl, Emma E. Davies, Antoinette B. Galvin

Abstract: We present an analysis of {\it in situ} and remote-sensing measurements of a coronal mass ejection (CME) that erupted on 2021 February 20 and impacted both the Solar TErrestrial RElations Observatory (STEREO)-A and the {\it Wind} spacecraft, which were separated longitudinally by 55$^\circ$. Measurements on 2021 February 24 at both spacecraft are consistent with the passage of a magnetic ejecta (M… ▽ More We present an analysis of {\it in situ} and remote-sensing measurements of a coronal mass ejection (CME) that erupted on 2021 February 20 and impacted both the Solar TErrestrial RElations Observatory (STEREO)-A and the {\it Wind} spacecraft, which were separated longitudinally by 55$^\circ$. Measurements on 2021 February 24 at both spacecraft are consistent with the passage of a magnetic ejecta (ME), making this one of the widest reported multi-spacecraft ME detections. The CME is associated with a low-inclined and wide filament eruption from the Sun's southern hemisphere, which propagates between STEREO-A and {\it Wind} around E34. At STEREO-A, the measurements indicate the passage of a moderately fast ($\sim 425$~km\,s$^{-1}$) shock-driving ME, occurring 2--3 days after the end of a high speed stream (HSS). At {\it Wind}, the measurements show a faster ($\sim 490$~km\,s$^{-1}$) and much shorter ME, not preceded by a shock nor a sheath, and occurring inside the back portion of the HSS. The ME orientation measured at both spacecraft is consistent with a passage close to the legs of a curved flux rope. The short duration of the ME observed at {\it Wind} and the difference in the suprathermal electron pitch-angle data between the two spacecraft are the only results that do not satisfy common expectations. We discuss the consequence of these measurements on our understanding of the CME shape and extent and the lack of clear signatures of the interaction between the CME and the HSS. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: accepted to ApJ

Multi-spacecraft observations of the structure of the sheath of an interplanetary coronal mass ejection and related energetic ion enhancement

Authors: E. K. J. Kilpua, S. W. Good, N. Dresing, R. Vainio, E. E. Davies, R. J. Forsyth, J. Gieseler, B. Lavraud, E. Asvestari, D. E. Morosan, J. Pomoell, D. J. Price, D. Heyner, T. S. Horbury, V. Angelini, H. O'Brien, V. Evans, J. Rodriguez-Pacheco, R. Gómez Herrero, G. C. Ho, R. Wimmer-Schweingruber

Abstract: Sheaths ahead of coronal mass ejections (CMEs) are large heliospheric structures that form with CME expansion and propagation. Turbulent and compressed sheaths contribute to the acceleration of particles in the corona and in interplanetary space, but the relation of their internal structures to particle energization is still relatively little studied. In particular, the role of sheaths in accelera… ▽ More Sheaths ahead of coronal mass ejections (CMEs) are large heliospheric structures that form with CME expansion and propagation. Turbulent and compressed sheaths contribute to the acceleration of particles in the corona and in interplanetary space, but the relation of their internal structures to particle energization is still relatively little studied. In particular, the role of sheaths in accelerating particles when the shock Mach number is low is a significant open problem. This work seeks to provide new insights on the internal structure of CME sheaths with regard to energetic particle enhancements. A good opportunity to achieve this aim was provided by observations of a sheath made by radially aligned spacecraft at 0.8 and $\sim$ 1 AU (Solar Orbiter, Wind, ACE and BepiColombo) on 19-21 April 2020. The sheath was preceded by a weak shock. Energetic ion enhancements occurred at different locations within the sheath structure at Solar Orbiter and L1. Magnetic fluctuation amplitudes at inertial-range scales increased in the sheath relative to the upstream wind. However, when normalised to the local mean field, fluctuation amplitudes did not increase significantly; magnetic compressibility of fluctuation also did not increase. Various substructures were embedded within the sheath at the different spacecraft, including multiple heliospheric current sheet (HCS) crossings and a small-scale flux rope. At L1, the ion flux enhancement was associated with the HCS crossings, while at Solar Orbiter, the enhancement occurred within the rope. Substructures that are swept from the upstream solar wind and compressed in the sheath can act as particularly effective acceleration sites. A possible acceleration mechanism is betatron acceleration associated with the small-scale flux rope and the warped HCS in the sheath. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: 14 pages, 12 figures; published in Astronomy & Astrophysics, Solar Orbiter First Results (Cruise Phase) special issue

Causes and Consequences of Magnetic Complexity Changes within Interplanetary Coronal Mass Ejections: a Statistical Study

Authors: Camilla Scolini, Réka M. Winslow, Noé Lugaz, Tarik M. Salman, Emma E. Davies, Antoinette B. Galvin

Abstract: We present the first statistical analysis of complexity changes affecting the magnetic structure of interplanetary coronal mass ejections (ICMEs), with the aim of answering the questions: How frequently do ICMEs undergo magnetic complexity changes during propagation? What are the causes of such changes? Do the in situ properties of ICMEs differ depending on whether they exhibit complexity changes?… ▽ More We present the first statistical analysis of complexity changes affecting the magnetic structure of interplanetary coronal mass ejections (ICMEs), with the aim of answering the questions: How frequently do ICMEs undergo magnetic complexity changes during propagation? What are the causes of such changes? Do the in situ properties of ICMEs differ depending on whether they exhibit complexity changes? We consider multi-spacecraft observations of 31 ICMEs by MESSENGER, Venus Express, ACE, and STEREO between 2008 and 2014 while radially aligned. By analyzing their magnetic properties at the inner and outer spacecraft, we identify complexity changes which manifest as fundamental alterations or significant re-orientations of the ICME. Plasma and suprathermal electron data at 1 au, and simulations of the solar wind enable us to reconstruct the propagation scenario for each event, and to identify critical factors controlling their evolution. Results show that ~65% of ICMEs change their complexity between Mercury and 1 au and that interaction with multiple large-scale solar wind structures is the driver of these changes. Furthermore, 71% of ICMEs observed at large radial (>0.4 au) but small longitudinal (<15 degrees) separations exhibit complexity changes, indicating that propagation over large distances strongly affects ICMEs. Results also suggest ICMEs may be magnetically coherent over angular scales of at least 15 degrees, supporting earlier theoretical and observational estimates. This work presents statistical evidence that magnetic complexity changes are consequences of ICME interactions with large-scale solar wind structures, rather than intrinsic to ICME evolution, and that such changes are only partly identifiable from in situ measurements at 1 au. △ Less

Submitted 24 November, 2021; originally announced November 2021.

A Catalog of Interplanetary Coronal Mass Ejections Observed by Juno between 1 and 5.4 AU

Authors: Emma E. Davies, Robert J. Forsyth, Réka M. Winslow, Christian Möstl, Noé Lugaz

Abstract: We use magnetic field measurements by the Juno spacecraft to catalog and investigate interplanetary coronal mass ejections (ICMEs) beyond 1 AU. During its cruise phase, Juno spent about 5 years in the solar wind between 2011 September and 2016 June, providing measurements of the interplanetary magnetic field (IMF) between 1 and 5.4 AU. Juno therefore presents the most recent opportunity for a stat… ▽ More We use magnetic field measurements by the Juno spacecraft to catalog and investigate interplanetary coronal mass ejections (ICMEs) beyond 1 AU. During its cruise phase, Juno spent about 5 years in the solar wind between 2011 September and 2016 June, providing measurements of the interplanetary magnetic field (IMF) between 1 and 5.4 AU. Juno therefore presents the most recent opportunity for a statistical analysis of ICME properties beyond 1 AU since the Ulysses mission (1990-2009). Our catalog includes 80 such ICME events, 32 of which contain associated flux-rope-like structures. We find that the dependency of the mean magnetic field strength of the magnetic flux ropes decreases with heliocentric distance as $r^{-1.24 \pm 0.43}$ between 1 and 5.4 AU, in good agreement with previous relationships calculated using ICME catalogs at Ulysses. We combine the Juno catalog with the HELCATS catalog to create a dataset of ICMEs covering 0.3-5.4 AU. Using a linear regression model to fit the combined dataset on a double-logarithmic plot, we find that there is a clear difference between global expansion rates for ICMEs observed at lesser heliocentric distances and those observed farther out beyond 1 AU. The cataloged ICMEs at Juno present a good basis for future multispacecraft studies of ICME evolution between the inner heliosphere, 1 AU, and beyond. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: Table 1 is available in its entirety in machine-readable format. Please contact the corresponding author if you wish to use the database

Multipoint interplanetary coronal mass ejections observed with Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A

Authors: C. Möstl, A. J. Weiss, M. A. Reiss, T. Amerstorfer, R. L. Bailey, J. Hinterreiter, M. Bauer, D. Barnes, J. A. Davies, R. A. Harrison, J. L. Freiherr von Forstner, E. E. Davies, D. Heyner, T. Horbury, S. D. Bale

Abstract: We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP),… ▽ More We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP), Wind and STEREO-A, in connection with coronagraph and heliospheric imaging observations from STEREO-A/SECCHI and SOHO/LASCO. We identify ICME events that could be unambiguously followed with the STEREO-A heliospheric imagers during their interplanetary propagation to their impact at the aforementioned spacecraft, and look for events where the same ICME is seen in situ by widely separated spacecraft. We highlight two events: (1) a small streamer blowout CME on 2020 June 23 observed with a triple lineup by PSP, BepiColombo and Wind, guided by imaging with STEREO-A, and (2) the first fast CME of solar cycle 25 ($ \approx 1600$ km s$^{-1}$) on 2020 November 29 observed in situ by PSP and STEREO-A. These results are useful for modeling the magnetic structure of ICMEs and the interplanetary evolution and global shape of their flux ropes and shocks, and for studying the propagation of solar energetic particles. The combined data from these missions are already turning out to be a treasure trove for space weather research and are expected to become even more valuable with an increasing number of ICME events expected during the rise and maximum of solar cycle 25. △ Less

Submitted 4 January, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: in press at ApJ Letters (submitted 2021 September 15, revised 2021 November 24, accepted 2021 December 9). 11 pages, 3 figures, 1 table

Solar origins of a strong stealth CME detected by Solar Orbiter

Authors: Jennifer O'Kane, Lucie M. Green, Emma E. Davies, Christian Möstl, Jürgen Hinterreiter, Johan L. Freiherr von Forstner, Andreas J. Weiss, David M. Long, Tanja Amerstorfer

Abstract: Aims.We aim to locate the origin of a stealth coronal mass ejection (CME) detected in situ by the MAG instrument on board Solar Orbiter, and make connections between the CME observed at the Sun, and the interplanetary CME (ICME) measured in situ. Methods. Remote sensing data are analysed using advanced image processing techniques to identify the source region of the stealth CME, and the global mag… ▽ More Aims.We aim to locate the origin of a stealth coronal mass ejection (CME) detected in situ by the MAG instrument on board Solar Orbiter, and make connections between the CME observed at the Sun, and the interplanetary CME (ICME) measured in situ. Methods. Remote sensing data are analysed using advanced image processing techniques to identify the source region of the stealth CME, and the global magnetic field at the time of the eruption is examined using Potential Field Source Surface (PFSS) models. The observations of the stealth CME at the Sun are compared with the magnetic field measured by the Solar Orbiter spacecraft, and plasma properties measured by the Wind spacecraft. Results. The source of the CME is found to be a quiet Sun cavity in the northern hemisphere. We find that the stealth CME has a strong magnetic field in situ, despite originating from a quiet Sun region with an extremely weak magnetic field. Conclusions. The interaction of the ICME with its surrounding environment is the likely cause of a higher magnetic field strength measured in situ. Stealth CMEs require multi-wavelength and multi-viewpoint observations in order to confidently locate the source region, however their elusive signatures still pose many problems for space weather forecasting. The findings have implications for Solar Orbiter observing sequences with instruments such as EUI that are designed to capture stealth CMEs △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 6 pages, 7 figures, accepted for publication in A&A

Multi point analysis of coronal mass ejection flux ropes using combined data from Solar Orbiter, BepiColombo and Wind

Authors: A. J. Weiss, C. Moestl, E. E. Davies, T. Amerstorfer, M. Bauer, J. Hinterreiter, M. Reiss, R. L. Bailey, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heiner, I. Richter, H-U. Auster, W. Magnes, D. Fischer, W. Baumjohann

Abstract: The recent launch of Solar Orbiter and BepiColombo opened a brief window in which these two spacecraft were positioned in a constellation that allows for the detailed sampling of any Earth-directed CMEs. Fortunately, two such events occurred with in situ detections of an ICME by Solar Orbiter on the 19th of April and the 28th of May 2020. These two events were subsequently also observed in situ by… ▽ More The recent launch of Solar Orbiter and BepiColombo opened a brief window in which these two spacecraft were positioned in a constellation that allows for the detailed sampling of any Earth-directed CMEs. Fortunately, two such events occurred with in situ detections of an ICME by Solar Orbiter on the 19th of April and the 28th of May 2020. These two events were subsequently also observed in situ by BepiColombo and Wind around a day later. We attempt to reconstruct the observed in situ magnetic field measurements for all three spacecraft simultaneously using an empirical magnetic flux rope model. This allows us to test the validity of our flux rope model on a larger and more global scale and allows for cross-validation of the analysis with different spacecraft combinations. Finally, we can also compare the results from the in situ modeling to remote observations obtained from the STEREO-A heliospheric imagers. We make use of the 3D coronal rope ejection model in order to simulate the ICME evolution. We adapt a previously developed ABC-SMC fitting algorithm for the application to multi point scenarios. We show that we are able to generally reconstruct the flux ropes signatures at three different spacecraft positions simultaneously using our model in combination with the flux rope fitting algorithm. For the well-behaved 19th of April ICME our approach works very well. The 28th of May ICME, on the other hand, shows the limitations of our approach. Unfortunately, the usage of multi-point observations for these events does not appear to solve inherent issues, such as the estimation of the magnetic field twist or flux rope aspect-ratios due to the specific constellation of the spacecraft positions. As our general approach can be used for any fast forward simulation-based model we give a blueprint for future studies using more advanced ICME models. △ Less

Submitted 17 May, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: 11 pages, 7 figures

Journal ref: A&A 656, A13 (2021)

Multi-spacecraft Study of the Solar Wind at Solar Minimum: Dependence on Latitude and Transient Outflows

Authors: R. Laker, T. S. Horbury, S. D. Bale, L. Matteini, T. Woolley, L. D. Woodham, J. E. Stawarz, E. E. Davies, J. P. Eastwood, M. J. Owens, H. O'Brien, V. Evans, V. Angelini, I. Richter, D. Heyner, C. J. Owen, P. Louarn, A. Federov

Abstract: The recent launches of Parker Solar Probe (PSP), Solar Orbiter (SO) and BepiColombo, along with several older spacecraft, have provided the opportunity to study the solar wind at multiple latitudes and distances from the Sun simultaneously. We take advantage of this unique spacecraft constellation, along with low solar activity across two solar rotations between May and July 2020, to investigate h… ▽ More The recent launches of Parker Solar Probe (PSP), Solar Orbiter (SO) and BepiColombo, along with several older spacecraft, have provided the opportunity to study the solar wind at multiple latitudes and distances from the Sun simultaneously. We take advantage of this unique spacecraft constellation, along with low solar activity across two solar rotations between May and July 2020, to investigate how the solar wind structure, including the Heliospheric Current Sheet (HCS), varies with latitude. We visualise the sector structure of the inner heliosphere by ballistically mapping the polarity and solar wind speed from several spacecraft onto the Sun's source surface. We then assess the HCS morphology and orientation with the in situ data and compare with a predicted HCS shape. We resolve ripples in the HCS on scales of a few degrees in longitude and latitude, finding that the local orientation of sector boundaries were broadly consistent with the shape of the HCS but were steepened with respect to a modelled HCS at the Sun. We investigate how several CIRs varied with latitude, finding evidence for the compression region affecting slow solar wind outside the latitude extent of the faster stream. We also identified several transient structures associated with HCS crossings, and speculate that one such transient may have disrupted the local HCS orientation up to five days after its passage. We have shown that the solar wind structure varies significantly with latitude, with this constellation providing context for solar wind measurements that would not be possible with a single spacecraft. These measurements provide an accurate representation of the solar wind within $\pm 10^{\circ}$ latitude, which could be used as a more rigorous constraint on solar wind models and space weather predictions. In the future, this range of latitudes will increase as SO's orbit becomes more inclined. △ Less

Submitted 22 June, 2021; v1 submitted 27 February, 2021; originally announced March 2021.

Comments: Accepted version

Journal ref: A&A 652, A105 (2021)

Radial Evolution of the April 2020 Stealth Coronal Mass Ejection between 0.8 and 1 AU -- A Comparison of Forbush Decreases at Solar Orbiter and Earth

Authors: Johan L. Freiherr von Forstner, Mateja Dumbović, Christian Möstl, Jingnan Guo, Athanasios Papaioannou, Robert Elftmann, Zigong Xu, Jan Christoph Terasa, Alexander Kollhoff, Robert F. Wimmer-Schweingruber, Javier Rodríguez-Pacheco, Andreas J. Weiss, Jürgen Hinterreiter, Tanja Amerstorfer, Maike Bauer, Anatoly V. Belov, Maria A. Abunina, Timothy Horbury, Emma E. Davies, Helen O'Brien, Robert C. Allen, G. Bruce Andrews, Lars Berger, Sebastian Boden, Ignacio Cernuda Cangas , et al. (18 additional authors not shown)

Abstract: Aims. We present observations of the first coronal mass ejection (CME) observed at the Solar Orbiter spacecraft on April 19, 2020, and the associated Forbush decrease (FD) measured by its High Energy Telescope (HET). This CME is a multispacecraft event also seen near Earth the next day. Methods. We highlight the capabilities of HET for observing small short-term variations of the galactic cosmic r… ▽ More Aims. We present observations of the first coronal mass ejection (CME) observed at the Solar Orbiter spacecraft on April 19, 2020, and the associated Forbush decrease (FD) measured by its High Energy Telescope (HET). This CME is a multispacecraft event also seen near Earth the next day. Methods. We highlight the capabilities of HET for observing small short-term variations of the galactic cosmic ray count rate using its single detector counters. The analytical ForbMod model is applied to the FD measurements to reproduce the Forbush decrease at both locations. Input parameters for the model are derived from both in situ and remote-sensing observations of the CME. Results. The very slow (~350 km/s) stealth CME caused a FD with an amplitude of 3 % in the low-energy cosmic ray measurements at HET and 2 % in a comparable channel of the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) on the Lunar Reconnaissance Orbiter, as well as a 1 % decrease in neutron monitor measurements. Significant differences are observed in the expansion behavior of the CME at different locations, which may be related to influence of the following high speed solar wind stream. Under certain assumptions, ForbMod is able to reproduce the observed FDs in low-energy cosmic ray measurements from HET as well as CRaTER, but with the same input parameters, the results do not agree with the FD amplitudes at higher energies measured by neutron monitors on Earth. We study these discrepancies and provide possible explanations. Conclusions. This study highlights that the novel measurements of the Solar Orbiter can be coordinated with other spacecraft to improve our understanding of space weather in the inner heliosphere. Multi-spacecraft observations combined with data-based modeling are also essential to understand the propagation and evolution of CMEs as well as their space weather impacts. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 656, A1 (2021)

In situ multi-spacecraft and remote imaging observations of the first CME detected by Solar Orbiter and BepiColombo

Authors: E. E. Davies, C. Möstl, M. J. Owens, A. J. Weiss, T. Amerstorfer, J. Hinterreiter, M. Bauer, R. L. Bailey, M. A. Reiss, R. J. Forsyth, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heyner, I. Richter, H-U. Auster, W. Magnes, W. Baumjohann, D. Fischer, D. Barnes, J. A. Davies, R. A. Harrison

Abstract: On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ i… ▽ More On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ in longitude. The source of the CME, which was launched on April 15th, was an almost entirely isolated streamer blowout. STEREO-A observed the event remotely from -75.1$^\circ$ longitude, which is an exceptionally well suited viewpoint for heliospheric imaging of an Earth directed CME. The configuration of the four spacecraft has provided an exceptionally clean link between remote imaging and in situ observations of the CME. We have used the in situ observations of the CME at Solar Orbiter, Wind, and BepiColombo, and the remote observations of the CME at STEREO-A in combination with flux rope models to determine the global shape of the CME and its evolution as it propagated through the inner heliosphere. A clear flattening of the CME cross-section has been observed by STEREO-A, and further confirmed by comparing profiles of the flux rope models to the in situ data, where the distorted flux rope cross-section qualitatively agrees most with in situ observations of the magnetic field at Solar Orbiter. Comparing in situ observations of the magnetic field between spacecraft, we find that the dependence of the maximum (mean) magnetic field strength decreases with heliocentric distance as $r^{-1.24 \pm 0.50}$ ($r^{-1.12 \pm 0.14}$), in disagreement with previous studies. Further assessment of the axial and poloidal magnetic field strength dependencies suggests that the expansion of the CME is likely neither self-similar nor cylindrically symmetric. △ Less

Submitted 23 February, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Journal ref: A&A 656, A2 (2021)

Three Dimensional Simulations of Solar Wind Preconditioning and the 23 July 2012 Interplanetary Coronal Mass Ejection

Authors: Ravindra T. Desai, Han Zhang, Emma E. Davies, Julia E. Stawarz, Joan Mico-Gomez, Pilar Iváñez-Ballesteros

Abstract: Predicting the large-scale eruptions from the solar corona and their propagation through interplanetary space remains an outstanding challenge in solar- and helio-physics research. In this article, we describe three dimensional magnetohydrodynamic simulations of the inner heliosphere leading up to and including the extreme interplanetary coronal mass ejection (ICME) of 23 July 2012, developed usin… ▽ More Predicting the large-scale eruptions from the solar corona and their propagation through interplanetary space remains an outstanding challenge in solar- and helio-physics research. In this article, we describe three dimensional magnetohydrodynamic simulations of the inner heliosphere leading up to and including the extreme interplanetary coronal mass ejection (ICME) of 23 July 2012, developed using the code PLUTO. The simulations are driven using the output of coronal models for Carrington rotations 2125 and 2126 and, given the uncertainties in the initial conditions, are able to reproduce an event of comparable magnitude to the 23 July ICME, with similar velocity and density profiles at 1 au. The launch-time of this event is then varied with regards to an initial 19 July ICME and the effects of solar wind preconditioning are found to be significant for an event of this magnitude and to decrease over a time-window consistent with the ballistic refilling of the depleted heliospheric sector. These results indicate that the 23 July ICME was mostly unaffected by events prior, but would have travelled even faster had it erupted closer in time to the 19 July event where it would have experienced even lower drag forces. We discuss this systematic study of solar wind preconditioning in the context of space weather forecasting. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 17 pages, 5 figures, 1 table. Solar Physics accepted 26 August 2020

Stacking Redshifted 21cm Images of HII Regions Around High Redshift Galaxies as a Probe of Early Reionization

Authors: James E Davies, Rupert A. Croft, Tiziana Di-Matteo, Bradley Greig, Yu Feng, Stuart Wyithe

Abstract: A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21cm line. Using the \textsc{BlueTides} simulation, we investigate the measurement of an \textit{average} ionised region towards the beginning of reionization by stacking redshifted 21cm images around optically identified bright galaxies u… ▽ More A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21cm line. Using the \textsc{BlueTides} simulation, we investigate the measurement of an \textit{average} ionised region towards the beginning of reionization by stacking redshifted 21cm images around optically identified bright galaxies using mock observations. We find that with an SKA 1000 hour observation, assuming perfect foreground subtraction, a $5σ$ detection of a stacked HII region can be made with 30 images around some of the brightest galaxies in \textsc{bluetides} (brighter than $M_{UV} < -22.75$) at $z=9$ (corresponding to a neutral fraction of 90.1 \% in our model). We present simulated relationships between the UV magnitude of galaxies, the sizes of the ionised regions they reside in, and the shape of the stacked profiles. These mock observations can also distinguish between scenarios where the IGM is in net emission or absorption of 21cm photons. Once 21cm foreground contamination is included, we find that even with up to 200 images around these rare, bright galaxies, only a tentative $> 1σ$ detection will be possible. However, partial foreground subtraction substantially improves signal-to-noise. For example, we predict that reducing the area of Fourier space dominated by foregrounds by 50 (80) percent will allow $> 3σ$ ($> 5σ$) detections of ionised regions at $z=9$. △ Less

Submitted 19 November, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 12 pages, 17 figures, Accepted for publication in MNRAS

Robustness of slow contraction to cosmic initial conditions

Authors: Anna Ijjas, William G. Cook, Frans Pretorius, Paul J. Steinhardt, Elliot Y. Davies

Abstract: We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spa… ▽ More We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results. △ Less

Submitted 8 July, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 41 pages, 18 figures; accepted for publication in JCAP

Journal ref: JCAP 08 (2020) 030

Prospective Lava Tubes at Hellas Planitia

Authors: Antonio Paris, Evan Davies, Laurence Tognetti, Carly Zahniser

Abstract: After surveying 1,500 images from MRO, this investigation has identified three candidate lava tubes in the vicinity of Hadriacus Mons as prospective sites for manned exploration. This investigation, therefore, concluded that terrestrial lava tubes can be leveraged for radiation shielding and, accordingly, that the candidate lava tubes on Mars (as well as known lava tubes on the lunar surface) can… ▽ More After surveying 1,500 images from MRO, this investigation has identified three candidate lava tubes in the vicinity of Hadriacus Mons as prospective sites for manned exploration. This investigation, therefore, concluded that terrestrial lava tubes can be leveraged for radiation shielding and, accordingly, that the candidate lava tubes on Mars (as well as known lava tubes on the lunar surface) can serve as natural radiation shelters and habitats for a prospective crewed mission to the planet. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: 16

Journal ref: Journal of the Washington Academy of Sciences, 2019

El Bahr: A Prospective Impact Crater

Authors: Antonia Paris, Shalabiea Osama M, Ahmed Mohamed, Evan Davies

Abstract: This preliminary investigation addresses the discovery of an unidentified crater located south of the Sahara Desert between Qaret Had El Bahr and Qaret El Allafa, Egypt. The unidentified crater (hereafter tentatively named El Bahr Crater) was discovered during a terrain analysis of the Sahara Desert. El Bahr Crater is located Southwest Al-Jiza Giza and is approximately 327 meters across, has a rim… ▽ More This preliminary investigation addresses the discovery of an unidentified crater located south of the Sahara Desert between Qaret Had El Bahr and Qaret El Allafa, Egypt. The unidentified crater (hereafter tentatively named El Bahr Crater) was discovered during a terrain analysis of the Sahara Desert. El Bahr Crater is located Southwest Al-Jiza Giza and is approximately 327 meters across, has a rim with a circumference of approximately 1,027 meters, and occupies a surface area of approximately 83,981 square meters. Preliminary spectral and topographic analysis reveal features characteristic of an impact crater produced by a hypervelocity event of extraterrestrial origin, including a bowl-shaped rim and a crater wall. No proximal and/or distal ejecta, however, are visible from Landsat imagery. Moreover, the geomorphic features, along with the fact that the El Bahr basalts are known to be rich in orthopyroxene while the surrounding basalts are not, imply an impact as the most plausible explanation. The El Bahr Crater is not indexed in the Earth Impact Database, and an analysis of impact structures in Africa did not identify it as either a confirmed, proposed or disproved impact crater. In collaboration with the University of Cairo, therefore, an expedition has been organized to conduct an in-situ investigation of El Bahr Crater, to ascertain if planar formations, shatter cones, and shock metamorphic and or other meteoritic properties are present. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: 7 Pages, 8 Figures, Peer-reviewed and published

Journal ref: Journal of the Washington Academy of Sciences, 2018

Silicate Melting and Vaporization during Rocky Planet Formation

Authors: Erik J. Davies, Phil J. Carter, Seth Root, Richard G. Kraus, Dylan K. Spaulding, Sarah T. Stewart, Stein B. Jacobsen$^{4}$

Abstract: Collisions that induce melting and vaporization can have a substantial effect on the thermal and geochemical evolution of planets. However, the thermodynamics of major minerals are not well known at the extreme conditions attained during planet formation. We obtained new data at the Sandia Z Machine and use published thermodynamic data for the major mineral forsterite (Mg$_2$SiO$_4$) to calculate… ▽ More Collisions that induce melting and vaporization can have a substantial effect on the thermal and geochemical evolution of planets. However, the thermodynamics of major minerals are not well known at the extreme conditions attained during planet formation. We obtained new data at the Sandia Z Machine and use published thermodynamic data for the major mineral forsterite (Mg$_2$SiO$_4$) to calculate the specific entropy in the liquid region of the principal Hugoniot. We use our calculated specific entropy of shocked forsterite, and revised entropies for shocked silica, to determine the critical impact velocities for melting or vaporization upon decompression from the shocked state to 1 bar and the triple points, which are near the pressures of the solar nebula. We also demonstrate the importance of the initial temperature on the criteria for vaporization. Applying these results to $N$-body simulations of terrestrial planet formation, we find that up to 20 to 40% of the total system mass is processed through collisions with velocities that exceed the criteria for incipient vaporization at the triple point. Vaporizing collisions between small bodies are an important component of terrestrial planet formation. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 15 pages, 8 figures, Accepted for publication in JGR: Planets

The Shock Physics of Giant Impacts: Key Requirements for the Equations of State

Authors: Sarah T. Stewart, Erik J. Davies, Megan S. Duncan, Simon J. Lock, Seth Root, Joshua P. Townsend, Richard G. Kraus, Razvan Caracas, Stein B. Jacobsen

Abstract: We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an upda… ▽ More We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model. △ Less

Submitted 16 January, 2020; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: M-ANEOS source code public release (https://github.com/isale-code/M-ANEOS); forsterite EOS public release (https://github.com/ststewart/aneos-forsterite-2019); manuscript accepted (no changes): Stewart, S., et al. (accepted). In J. Lane, T. Germann, and M. Armstrong (Eds.), 21st Biennial APS Conference on Shock Compression of Condensed Matter (SCCM19). AIP Publishing

Fuzzy Dark Matter Soliton Cores around Supermassive Black Holes

Authors: Elliot Yarnell Davies, Philip Mocz

Abstract: We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing' effect on the soliton density profile, decreasing the core… ▽ More We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a `squeezing' effect on the soliton density profile, decreasing the core radius and increasing the central density. In the limit of large black hole mass, the solution approaches one akin to the hydrogen atom, with radius inversely proportional to the black hole mass. By applying our analysis to two specific galaxies (M87 and the Milky Way) and pairing it with known observational limits on the amount of centrally concentrated dark matter, we obtain a constraint on the FDM particle mass, finding that the range $10^{-22.12}$ eV $\lesssim m \lesssim 10^{-22.06}$ eV should be forbidden (taking into account additional factors concerning the life-time of the soliton in the vicinity of a black hole). Improved observational mass measurements of the black hole and total enclosed masses will significantly extend the lower-bound on the excluded FDM mass region, while self-consistent theoretical modeling of the soliton--black hole system can extend the upper-bound. △ Less

Submitted 22 January, 2020; v1 submitted 13 August, 2019; originally announced August 2019.

Comments: 10 pages, 7 figures, accepted to MNRAS

Dark-Ages Reionisation & Galaxy Formation Simulation XVI: The Thermal Memory of Reionisation

Authors: James E. Davies, Simon J. Mutch, Yuxiang Qin, Andrei Mesinger, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how diff… ▽ More Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how different thermal statistics of the intergalactic medium can be used to constrain reionisation. Temperature is highly correlated with redshift of reionisation for a period of time after the gas is heated. However as the gas cools, thermal memory of reionisation is lost, and a power-law temperature-density relation is formed, $T = T_0(1+δ)^{1-γ}$ with $γ\approx 1.5$. Constraining our model against observations of electron optical depth and temperature at mean density, we find that reionisation likely finished at $z_{\rm{reion}} = 6.8 ^{+ 0.5} _{-0.8}$ with a soft spectral slope of $α= 2.8 ^{+ 1.2} _{-1.0}$. By restricting spectral slope to the range $[0.5,2.5]$ motivated by population II synthesis models, reionisation timing is further constrained to $z_{\rm{reion}} = 6.9 ^{+ 0.4} _{-0.5}$. We find that, in the future, the degeneracies between reionisation timing and background spectrum can be broken using the scatter in temperatures and integrated thermal history. △ Less

Submitted 23 August, 2019; v1 submitted 12 April, 2019; originally announced April 2019.

Comments: 17 pages, 17 figures, Accepted for publication in MNRAS

Hydrogen Clouds from Comets 266P Christensen and P2008 Y2 (Gibbs) are Candidates for the Source of the 1977 WOW!Signal

Authors: Antonio Paris, Evan Davies

Abstract: On 1977 August 15, the Ohio State University Radio Observatory detected a strong narrowband signal northwest of the globular star cluster M55 in the constellation Sagittarius. The frequency of the signal, which closely matched the hydrogen line (1420.40575177 MHz), peaked at approximately 2316 EDT. Since then, several investigations into the Wow signal have ruled out the source as terrestrial in o… ▽ More On 1977 August 15, the Ohio State University Radio Observatory detected a strong narrowband signal northwest of the globular star cluster M55 in the constellation Sagittarius. The frequency of the signal, which closely matched the hydrogen line (1420.40575177 MHz), peaked at approximately 2316 EDT. Since then, several investigations into the Wow signal have ruled out the source as terrestrial in origin or other objects such as satellites, planets and asteroids. From 1977 July 27 to 1977 August 15, comets 266P Christensen and P2008 Y2 (Gibbs) were transiting in the neighborhood of the Chi Sagittarii star group. Because the frequency for the Wow! signal fell close to the hydrogen line, and the hydrogen clouds of 266P/Christensen and P/2008 Y2 (Gibbs) were in the proximity of the right ascension and declination values of the Wow signal, the comet(s) and or their hydrogen clouds are strong candidates for the source of the 1977 Wow signal. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 7 pages

Report number: Volume 101, Number 4, Winter 2015

Journal ref: Journal of the Washington Academy of Sciences, Volume 101, Number 4, Winter 2015

Metallicity and Age of the Stellar Stream Around the Disk Galaxy NGC 5907

Authors: Seppo Laine, Carl J. Grillmair, Peter Capak, Richard G. Arendt, Aaron J. Romanowsky, David Martinez-Delgado, Matthew L. N. Ashby, James E. Davies, Stephen R. Majewski, Jean P. Brodie, R. Jay GaBany, Jacob A. Arnold

Abstract: Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron d… ▽ More Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream. △ Less

Submitted 12 July, 2016; originally announced July 2016.

Comments: 16 pages, 11 figures, 1 table, accepted for publication in the Astronomical Journal