Photon rings of spherically symmetric black holes and robust tests of non-Kerr metrics

Maciek Wielgus∗ Black Hole Initiative at Harvard University, 20 Garden Street, Cambridge, MA 02138, USA and Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA (Dated: September 23, 2021) Under very general assumptions on the accretion flow geometry, images of a black hole illuminated by electromagnetic radiation display a sequence of photon rings (demagnified and rotated copies of the direct image) which asymptotically approach a purely theoretical critical curve – the outline of the black hole photon shell. To a distant observer, these images appear dominated by the direct emission, which forms a ring whose diameter is primarily determined by the effective radius of the emitting region. For that reason, connecting the image diameter seen by a distant observer to the properties of the underlying spacetime crucially relies on a calibration that necessarily depends on the assumed astrophysical source model. On the other hand, the diameter of the photon rings depends more on the detailed geometry of the spacetime than on the source structure. As such, a photon ring detection would allow for the spacetime metric to be probed in a less model-dependent way, enabling more robust tests of General Relativity (GR) and the Kerr hypothesis. Here we present the photon ring structure of several spherically symmetric black hole spacetimes and perform comparisons with the Schwarzschild/Kerr case. We offer our perspective on future tests of the spacetime metric with photon rings, discussing challenges and opportunities involved.

I. INTRODUCTION Additionally, the observable ring-like direct emission fea- ture is commonly conflated with the critical curve in the The observational appearance of a black hole (BH) has literature. One of the aims of this paper is to clarify that been studied extensively for many years from a purely issue. theoretical standpoint [1–4]. This appearance depends As noticed by [10], the technique of very long base- on numerous astrophysical details, such as geometry and line interferometry (VLBI) is uniquely suited to probe optical depth of the emission region [5–7]. The feature the photon ring structure of black hole images. As sharp depending exclusively on the spacetime geometry is the image domain features, photon rings die-off slowly in the critical curve [8], an outline of the photon shell on a dis- Fourier domain, dominating the interferometric signal in tant observer’s screen. However, since the null geodesics the range of very high spatial frequencies, corresponding comprising the photon shell are unstable orbits [9], the to baseline lengths of tens and hundreds billion wave- critical curve remains a purely mathematical concept. lengths. While such long baselines were never observed Unless the accretion flow itself is spherically symmetric so far and this feat will surely require significant progress [4,6] (which is typically not an astrophysically plausible and new developments in the field of space VLBI, the scenario), the image of a BH contains a sequence of dis- projects that would enable robust detection of the pho- crete photon rings [8, 10]. Each photon ring is a copy ton ring structure are already being discussed and pro- of the direct image of the emission region, corresponding posed, e.g., [7, 10, 11, 17–19], and algorithms to enable to photon trajectories lensed around the photon shell, photon ring detection in VLBI data are already being with the n-th photon ring photons executing n/2 loops implemented and tested with synthetic data [11, 20, 21]. before reaching the observer’s screen [5, 10], see Fig.1. Only recently first resolved observational images of a This sequence of photon rings quickly converges to the BH were obtained by the Event Horizon Telescope VLBI critical curve, and the diameter and shape of the n = 2 array (EHT; [22–27]), presenting the supermassive BH photon ring is already strongly dominated by the space- arXiv:2109.10840v1 [gr-qc] 22 Sep 2021 M 87∗ in the center of the M 87 galaxy with a resolution time geometry – hence it can be used for precise tests of corresponding to 3 r (Schwarzschild radii; [24–26]). the theory of gravity, with only very general assumptions S The measured diameter∼ of a bright ring feature found on the astrophysical model of the accretion flow [10, 11]. in the EHT images was used to test the null hypoth- In the case of non-Kerr spacetimes, most of the previ- esis of consistency with the Kerr metric [27, 28] under ous research has focused on the theoretical critical curve the BH mass prior derived from stellar dynamics [29]. (commonly referred to as shadow, occasionally also as To formulate this test, the scaling relationship between photon ring, which has a different meaning in this paper), the ring diameter and the BH mass assuming Kerr met- which has been characterized for a broad range of non- ric was evaluated using 103 synthetic images generated Kerr spacetimes, e.g., [12–16]. Relatively less work has with general relativistic magnetohydrodynamic simula- been dedicated to the observable photon ring structure. tions (GRMHD; [26, 27]). These synthetic images, rep- resenting variety of viable theoretical GRMHD models of BH accretion, allowed to perform diameter-mass calibra- tion and quantify the null hypothesis test. The results of ∗ [email protected] this test were later on reinterpreted as EHT constraints IFT-UAM/CSIC-21-100

Scalar Cerenkovˇ radiation from high-energy cosmic rays

Charles Dalang,1, ∗ Pierre Fleury,2, † and Lucas Lombriser1, ‡ 1Université de Genève, Département de Physique Théorique and Center for Astroparticle Physics, 24 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland 2Instituto de Física Teórica UAM-CSIC, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain As first noted by Robert Wagoner in the nineteen-seventies, if a scalar field is nonminimally coupled to the Ricci scalar and propagates at subluminal speeds, then there exists the possibility of scalar Cˇ erenkov radiation from a moving particle. The mere observation of high-energy cosmic rays could in principle rule out the existence of such scalar fields since any particle moving faster than scalar perturbations would lose energy in the form of scalar waves until it moves slower than those. We compute in detail the energy loss to scalar waves and find that it scales with the square of the ultra-violet (UV) cutoff frequency of the effective field theory (EFT) of gravity. For dark-energy-motivated EFTs, the UV cutoff can be low, in which case that energy loss could always be negligible. In contrast, if viewed as a covariant theory valid at all scales or as an EFT valid at higher energies, perhaps even all the way up to the Planck scale, as may be the case if motivated by quantum-gravity perspectives, then the energy loss to scalar waves may diverge or become dramatically large. In this case, high-energy cosmic rays of extragalactic origin stringently constrain any conformally coupled scalar fields with non-canonical kinetic terms.

I. INTRODUCTION of modified gravity and dark energy [15–19]. This constraint arises because many of these scalar-tensor interactions im- here exists, in the context of gravitational physics, an ply non-luminal propagation of GWs due to the nontrivial T ambiguity between missing matter and describing the time dependence of the background scalar field configuration. Universe with an incomplete theory of gravity. As an example, Another, perhaps even more exciting consequence of the ob- unexplained perturbations of the orbit of Uranus led the French servation of gravitational waves, is that those do not seem to astronomer A. Bouvard to predict, in 1821, the existence of a decay into scalar waves, as predicted by degenerate higher- missing planet, known as Neptune, which was observed and order scalar-tensor (DHOST) theories, thereby eliminating their identified as such for the first time more than two decades later cosmological viability [20]. by J. C. Adams, in 1843 [1]. Conversely, it was working with As a follow-up question and since those constraints have the inadequate theory of Newtonian gravity and the observation been so powerful, can we learn anything about gravitational of the perihelion advance of Mercury that led U. Le Verrier to theories from the propagation of other fields such as high- think of a missing unobserved planet, named Vulcan, which energy cosmic particles? It turns out that for metric theories this time, was never discovered [2]. The theory of gravity which possess a direct coupling of the scalar field to the Ricci was incomplete and a more fundamental theory of gravity, scalar, scalar perturbations are sourced by the trace of the general relativity, was shown to solve this missing-matter issue. energy-momentum tensor of matter fields. Their equation of With this historical perspective in mind, and the contemporary motion is reminiscent of the electromagnetic potential which problems of cosmology, such as the inferred missing dark matter is sourced by electrically charged particles that lose energy and mysterious dark energy when adopting general relativity, to electromagnetic waves if they travel faster than the speed it seems reasonable to question the currently accepted theory of light in a medium. This emitted electromagnetic signal, of gravitation. known as Cˇ erenkov radiation [21] compensates the energy loss Alternative theories of gravity of present interest include per unit distance of those particles, as predicted by I. Franck scalar-tensor theories and their generalizations, including Horn- and I. Tamm [22], leading to their Nobel Prize shared with deski theories [3] and beyond [4–7]. A subset of those are P. Cˇ erenkov in 1958. Tight constraints on the speed of spin-2 arXiv:2109.10812v1 [gr-qc] 22 Sep 2021 able to affect cosmological dynamics and evade, via so-called GWs via GW Cˇ erenkov radiation have already been inferred screening mechanisms [8–10], tight precision tests of gravity from the observation of high-energy particles [23]. An obvious which have been conducted within the Solar System [11, 12] follow-up question is then whether high-energy particles lose and within galactic environments with binary pulsars [13, 14], their energy by emitting scalar waves if they travel faster than for example. A powerful test of those theories, which circum- those [24]. After all, high-energy particles are observed all vents the problem of screening, comes from the measurement the time. Could it be that their mere observation rules out a of the speed of gravitational waves (GWs) through cosmologi- significant class of scalar-tensor theories? cal distances, which practically eliminates a large parameter In Sec. II, we explore for the first time the energy loss per unit space of gravitational modifications for effective field theories distance of a massive particle to scalar waves due to it traveling faster than the speed of scalar waves in a conformally coupled k-essence theory [25, 26]. In Sec. III, we discuss constraints on the ultra-violet (UV) cutoff of the effective field theory (EFT) ∗ [email protected] of gravity for conformally coupled k-essence models and also † pierre.fl[email protected] discuss more generic phenomenological models. Finally, we ‡ [email protected] present conclusions in Sec. IV. DESY 21-141, TTP21-030, P3H-21-062

Stellar Evolution confronts Axion Models

Luca Di Luzioa,b,c, Marco Fedeled, Maurizio Giannottie, Federico Mesciaf , Enrico Nardig

a Dipartimento di Fisica e Astronomia ‘G. Galilei’, Universit`adi Padova, Italy b Istituto Nazionale Fisica Nucleare, Sezione di Padova, Italy c DESY, Notkestraße 85, D-22607 Hamburg, Germany d Institut f¨urTheoretische Teilchenphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany e Physical Sciences, Barry University, 11300 NE 2nd Ave., Miami Shores, FL 33161, USA f Departament de F´ısica Qu`antica i Astrof´ısica, Institut de Ci`enciesdel Cosmos (ICCUB), Universitat de Barcelona, Mart´ıi Franqu`es1, E-08028 Barcelona, Spain g INFN, Laboratori Nazionali di Frascati, C.P. 13, 100044 Frascati, Italy

Axion production from astrophysical bodies is a topic in continuous develop- ment, because of theoretical progress in the estimate of stellar emission rates and, especially, because of improved stellar observations. We carry out a com- prehensive analysis of the most informative astrophysics data, revisiting the bounds on axion couplings to photons, nucleons and electrons, and reassessing the significance of various hints of anomalous stellar energy losses. We con- front the performance of various theoretical constructions in accounting for these hints, while complying with the observational limits on axion couplings. We identify the most favorable models, and the regions in the mass/couplings parameter space which are preferred by the global fit. Finally, we scruti- nize the discovery potential for such models at upcoming helioscopes, namely IAXO and its scaled versions. arXiv:2109.10368v1 [hep-ph] 21 Sep 2021

1

Contents

1 Introduction 3

2 Axion effective Lagrangian 4

3 Astrophysical observables and cooling anomalies 5 3.1 Axion production in stars ...... 6 3.2 Axion Bounds from White Dwarfs ...... 8 3.2.1 White Dwarf luminosity function ...... 9 3.2.2 White Dwarf Variables ...... 9 3.3 Axion Bounds from the tip of RGB stars in globular cluster ...... 10 3.4 Helium burning stars ...... 11 3.5 Supernovae ...... 12 3.6 Neutron stars ...... 12 3.7 Model-independent fit ...... 13

4 A representative sample of axion models 14 4.1 Universal DFSZ models ...... 15 4.1.1 DFSZ1 ...... 16 4.1.2 DFSZ2 ...... 16 4.2 Non-universal DFSZ models ...... 16 4.2.1 Non-universal 2+1 models ...... 17 4.2.2 1+1+1 models ...... 20 4.2.3 3HDM model ...... 21 4.3 Summary of axion models ...... 22

5 Axion models in the light of cooling anomalies 23 5.1 Axion fits to cooling anomalies ...... 23 5.2 Discovery potential of meV-scale axion experiments ...... 26

6 Summary and conclusions 31

A E/N factors in non-universal DFSZ models 32

2 Pole inflation from non-minimal coupling to gravity

Sotirios Karamitsos and Alessandro Strumia

Dipartimento di Fisica, Universit`adi Pisa, Pisa, Italia

Abstract

Transforming canonical scalars to the Einstein frame can give a multi- field generalization of pole inflation (namely, a scalar with a divergent kinetic term) at vanishing field-dependent Planck mass. However, to obtain an attractor, the scalar potential must obey certain non-generic conditions. These are automatically satisfied in Quantum Field The- ories with dimension-less couplings. The resulting models of pole in- flation have special inflationary predictions determined by the full RG running of couplings. Acceptable predictions for the tensor/scalar ra- tio arise for perturbative but moderately large couplings, so we explore the possible QFT runnings: to confinement, to an IR fixed point, and to a UV fixed point.

Contents

1 Introduction2 2 Overview of pole inflation3 3 Multi-field pole inflation from graviton canonicalisation5 3.1 Two-field pole inflation: hyperbolic singularity ...... 9 3.2 Two-field pole inflation: elliptic singularity ...... 9 3.3 Two-field pole inflation: linear singularity ...... 10 3.4 Point singularity ...... 10 4 Pole inflation in dimension-less quantum theories 11 arXiv:2109.10367v1 [hep-th] 21 Sep 2021 4.1 Inflation with dimension-less theories at small coupling ...... 12 4.2 Inflation with strong coupling in the UV ...... 14 4.3 Inflation with strong coupling in the IR ...... 16 4.4 Inflation with an IR fixed-point ...... 16 4.5 Inflation with an UV fixed-point ...... 19 5 Conclusions 20

1 Prepared for submission to JCAP

Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields

Anna Ijjas,a,b,1 Frans Pretorius,c Paul J. Steinhardtc and David Garfinkled

aMax Planck Institute for Gravitational Physics, Hanover, 30167, Germany bInstitute for Theory and Computation, Harvard University, Cambridge, MA, 02138, USA cDepartment of Physics, Princeton University, Princeton, NJ, 08544, USA dDepartment of Physics, Oakland University, Rochester, MI 48309, USA E-mail: [email protected]

Abstract. We present non-perturbative numerical relativity simulations of slowly contract- ing spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble den- sity fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann- Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant. arXiv:2109.09768v1 [gr-qc] 20 Sep 2021

1Corresponding Author Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra Application: Chemical Homogeneity in Open Cluster M67

Aarya A. Patil,1, 2 Jo Bovy,1, 2 Gwendolyn Eadie,1, 3 and Sebastian Jaimungal3

1David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto ON M5S 3H4, Canada 2Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St George Street, Toronto, ON M5S 3H4, Canada 3Department of Statistical Sciences, University of Toronto, 100 St George St, Toronto, ON M5S 3G3, Canada

(Received August 22, 2021; Revised September 21, 2021) Submitted to ApJ

ABSTRACT High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic Archaeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18, 933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 Functional Principal Components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe . 0.02 dex; C . 0.03 dex; O, Mg, Si, Ni . 0.04 dex; Ca . 0.05 dex; N, Al . 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

Keywords: Astrostatistics (1882); Spectroscopy (1558); Chemical abundances (224); Open star clusters (1160); Galaxy structure (622); Milky Way disk (1050); Dimensionality reduction (1943)

1. INTRODUCTION 2017), GALAH (De Silva et al. 2015), and Gaia-ESO The physical processes that govern galaxy disk for- (Gilmore et al. 2012) provide medium to high-resolution 5 mation and evolution encode themselves in the struc- (R ≈ 10, 000) spectra of ∼10 stars over a wide range ture of stellar populations. Discerning this structure of Galactocentric distances (3 < R < 15 kpc) and dis- arXiv:2109.10891v1 [astro-ph.GA] 22 Sep 2021 places key constraints on models of galactic disks and tances from the midplane (|z| < 1 kpc). Future surveys provides data for cosmological simulations to replicate. such as WEAVE (Dalton et al. 2012), MOONS (Cira- In the Milky Way, we have the ability to resolve individ- suolo et al. 2014), SDSS-V (Kollmeier et al. 2017), and ual stars and characterize disk stellar populations using 4MOST (de Jong et al. 2019) will increase this number large samples (Rix & Bovy 2013). Modern wide-field to millions. Using these data, we can estimate radial spectroscopic surveys such as APOGEE (Majewski et al. velocities, stellar parameters, chemical abundances, and ages for a high proportion of Galactic disk stars. Accurate and precise stellar abundance estimates Corresponding author: Aarya A. Patil through high-resolution spectroscopy allow us to probe [email protected] the detailed chemical and dynamical evolution of the Milky Way. We can understand this using a simplified MNRAS 000,1–13 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

MUSE observations of the ionized gas in a galaxy group at z = 0.45 with 21-cm H i absorption towards PKS 1610-771

Simon Weng,1,2,3? Elaine M. Sadler,1,2,3 Caroline Foster,1,3 Céline Péroux,4,5 Elizabeth K. Mahony,2 James R. Allison,3,6 Vanessa A. Moss,2,1 Renzhi Su,2,7 Matthew T. Whiting2 and Hyein Yoon1,3 1Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 2ATNF, CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW 1710, Australia 3ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) 4European Southern Observatory, Karl-Schwarzschildstrasse 2, D-85748 Garching bei München, Germany 5Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, F-13388 Marseille, France 6Sub-Dept. of Astrophysics, Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Rd., Oxford, OX1 3RH, UK 7Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We present results from MUSE observations following up on a 21-cm H i absorption system detected with the Australian Square Kilometre Array Pathfinder radio tele- scope at redshift 0.4503 towards the quasar PKS 1610-771. This H i absorber has a 20 −2 column density NH i = 2.7 0.1 10 [Ts/100K] cm , making it a likely Damped Lyman-α (DLA) system. We± identify× a· galaxy group· with four members (A, B, X and Y) at the same redshift as the H i absorption system, with impact parameters ranging from less than 10 kpc to almost 200 kpc from the quasar sightline. Ca ii and Na i ab- sorption is also seen in the MUSE spectrum of the background QSO, with velocities coinciding with the initial H i 21-cm detection, but tracing less dense and potentially warmer gas. This metal-line component aligns with the rotating ionized disc of galaxy B (impact parameter 18 kpc from the QSO) and appears to be co-rotating with the galaxy disc, although outflowing gas cannot be directly excluded. In contrast, the 21- cm H i absorber is blueshifted relative to the galaxies nearest the absorber and has the opposite sign to the velocity field of galaxy B. Since galaxies A and B are separated by only 17 kpc on the sky and 70 km s−1 in velocity, it is likely that the 21-cm de- tection traces extragalactic clouds of gas formed from their interaction. This system represents a first case study of the cold gas detected in galaxy groups by future large 21-cm absorption surveys, such as the First Large Absorption Survey in H i. Key words: intergalactic medium – galaxies: kinematics and dynamics – quasars: absorption lines – radio lines: ISM arXiv:2109.10875v1 [astro-ph.GA] 22 Sep 2021 1 INTRODUCTION affect gas abundance such as inflows and outflows also re- quire characterisation from observations as we pursue more The reservoirs of cold gas within galaxies serve as fuel for complete simulations of the baryon cycle. star formation, which in turn drives the morphological and kinematic properties of galaxies. It follows that the avail- The amount and distribution of neutral hydrogen (H i) ability of hydrogen and phenomena that affect this availabil- is well-constrained at low redshift (z . 0.1) from H i ity play an important role in how galaxies evolve (Péroux emission-line surveys Rosenberg & Schneider(2002); Zwaan & Howk 2020). Placing observational constraints on the et al.(2005); Hoppmann et al.(2015); Giovanelli & Haynes amount and distribution of gas across all redshifts will im- (2015); Jones et al.(2018). Beyond this redshift, 21-cm emis- prove our understanding of galaxy evolution. Processes that sion is too faint for detection without stacking techniques (e.g. Kanekar et al. 2016; Rhee et al. 2018; Chowdhury et al. 2020) as the transition rate is highly forbidden. While ? E-mail: [email protected] (KTS) H i emission line stacking can extend the search for neu-

© 2021 The Authors Prepared for submission to JCAP

Clustering of Gravitational Wave and Supernovae events: a multitracer analysis in Luminosity Distance Space

S. Libanore,a,b M.C. Artale,c D. Karagiannis,d M. Liguori,a,b N. Bartolo,a,b Y. Bouffanais,a,b M. Mapelli,a,b,e S. Matarresea,b,e,f

aDipartimento di Fisica e Astronomia G. Galilei, Università degli Studi di Padova Via Marzolo 8, 35131 Padova, Italy bINFN, Sezione di Padova, Via Marzolo 8, I−35131, Padova, Italy cInstitut fur¨ Astro- und Teilchenphysik, Universit¨at Innsbruck Technikerstrasse 25/8, 6020 Innsbruck, Austria dDepartment of Physics and Astronomy, University of the Western Cape Cape Town 7535, South Africa eINAF, Osservatorio Astronomico di Padova Vicolo dell’Osservatorio 5, I−35122, Padova, Italy f Gran Sasso Science Institute (GSSI), Viale F. Crispi 7, I−67100, L’Aquila, Italy E-mail: [email protected]

Abstract. We study the clustering of Gravitational Wave (GW) merger events and Super- novae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar veloci- ties and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodinamical N-body simulations to determine the bias fiducial model and considering future observa- tions from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a

arXiv:2109.10857v1 [astro-ph.CO] 22 Sep 2021 3 detector network configuration. We find that adding SN to the GW merger dataset con- siderably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case. MNRAS 000,1–13 (2022) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

A Differential RR Lyrae Line-of-Sight Distance Between M31 and M33

Ata Sarajedini,1★ 1Department of Physics, Florida Atlantic University, Boca Raton, Florida, USA

Accepted 2021 September 22. Received 2021 September 21; in original form 2021 September 9

ABSTRACT We present a purely differential line-of-sight distance between M31 and M33 using ab-type RR Lyrae variables observed in each galaxy by the Hubble Space Telescope Advanced Camera for Surveys in the F606W filter. Using 1501 RR Lyraes in 13 M31 fields and 181 RR Lyraes in 6 M33 fields, and placing all of these stars on a uniform photometric scale with internally consistent corrections for metal abundance and extinction, we find a relative absolute distance modulus of Δ(m-M)표 = –0.298 ± 0.016 in the sense of (m-M)표,푀31 – (m-M)표,푀33. Adopting an absolute distance modulus of (m-M)표=24.46±0.10 for M31 places M33 115 kpc beyond M31 in line-of-sight distance. Key words: Stars:horizontal branch, stars:variables:RR Lyrae, galaxies:interactions

1 INTRODUCTION

The distance of a celestial object is arguably its most important astronomical property. In the case of Local Group galaxies, there is an extensive literature on efforts to use the distances of the Large Magellanic Cloud (LMC) and M31 as anchors in the cosmic distance ladder (e.g. Riess et al. 2012; de Grijs et al. 2014; de Grijs & Bono 2014). In addition, relative distance measurements of M31 and M33 have figured prominently in the realization that M33 appears to be a satellite of M31 (i.e. the Andromeda system) and that their past dynamical interactions have influenced their detailed morphological and stellar population properties (McConnachie et al. 2009, Tepper-Garcia et al. 2020). In fact, the present-day distance between M31 and M33 is one of the most important parameters in constraining dynamical models of the system (e.g. Semczuk et al. 2018; van der Marel et al. 2012; 2019, Tepper-Garcia et al. 2020), which, in the big picture, can then shed light on the details of the ΛCDM hierarchical formation scenario (e.g. Lewis et al. 2013). A number of different and complimentary techniques have been employed to infer the distances of M31 and M33. de Grijs & Bono (2014) present a comprehensive summary of these efforts. They surveyed the literature from 1918 to 2014 and compiled a list of 168 distance measurements to M31 and 90 for M33. Based on these lists and the techniques employed in each case, de Grijs & Bono (2014) derive a “recommended" distance modulus of (m-M)표=24.46±0.10 for M31 and (m-M)표=24.67±0.07 for M33 leading to a difference of –0.21±0.12 mag in distance modulus. More recent observations of distance indicators to M33 and M31 such as those by Tanakul et al. (2017) and Tanakul & Sarajedini (2018), respectively, using RR Lyrae variables find consistent results with those from de Grijs & Bono (2014). These studies are discussed more below. In this work, we present the results of a purely differential method to measure the line-of-sight distance between M31 and M33 using ab-type RR Lyrae variables in both galaxies. This paper is organized as follows. The next section describes the nature of the observations employed herein while section 3 presents details of the photometric corrections made to the observations to ensure a uniform and internally consistent photometric scale. The results are presented in Section 4 and comparisons are made to previous work in Section 5. The conclusions arXiv:2109.10853v1 [astro-ph.GA] 22 Sep 2021 are presented in Section 6.

2 OBSERVATIONS

Details of the data-sets used in the present study are given in Table 1, which lists the publications from which the data for M31 and M33 ab-type RR Lyrae variables are taken. The sky locations of the M31 fields are shown in Figs. 1 and 2 of Tanakul & Sarajedini (2018), while the M33 fields are illustrated in Fig. 1 of Tanakul et al. (2017). We are focused exclusively on observations obtained by the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) taken in the F606W filter. Because the photometric performance of ACS has been closely tracked since it began operations on HST, we can make detailed corrections for the variation in its photometric zeropoint over time. As such, this is likely to be the most precise photometry available for these stars.

★ E-mail: [email protected]

© 2022 The Authors The NASA Multi-Messenger Astrophysics Science Support Center

Rita M. Sambrunaa,∗, Joshua E. Schliedera,∗, Daniel Kocevskib,∗, Regina Caputoa, Michelle C. Huib, Craig B. Markwardta, Brian P. Powella, Judith L. Racusina, Christopher Robertsa, Leo P. Singera, Alan P. Smalea, Tonia M. Ventersa, Colleen A. Wilson-Hodgeb

aNASA Goddard Space Flight Center, 8800 Greenbelt Rd., Greenbelt, 20771, MD, USA bNASA Marshall Space Flight Center, Martin Rd. SW, Huntsville, 35808, AL, USA

Abstract The era of multi-messenger astrophysics has arrived, leading to key new discoveries and revealing a need for coordination, collab- oration, and communication between world-wide communities using ground and space-based facilities. To fill these critical needs, NASA’s Goddard Space Flight Center and Marshall Space Flight Center are jointly proposing to establish a virtual Multi-Messenger Astrophysics Science Support Center that focuses entirely on community-directed services. In this article, we describe the baseline plan for the virtual Support Center which will position the community and NASA as an Agency to extract maximum science from multi-messenger events, leading to new breakthroughs and fostering increased coordination and collaboration. Keywords: Astroinformatics, Astronomy web services, High energy astrophysics, Gravitational waves, Neutrino astronomy, Cosmic ray astronomy

1. Introduction Center (SSC), with 100% community-directed services. Here we describe the baseline plan for the virtual MMA SSC which Multi-messenger astrophysics (MMA) has come of age will provide: thanks to the detection of gravitational wave (GW) sources from the ground with the Advanced LIGO (LIGO Scientific Collabo- • A community access portal (CAP) website to act as a one- ration et al., 2015) and Virgo (Acernese et al., 2015) observato- stop-shop for information, tools, and support. ries, and of an extra-galactic neutrino source with the IceCube Neutrino Observatory (Aartsen et al., 2017). Together with the • A real-time board for instant community collaboration, co- concurrent observations of coincident gamma-ray photons fol- ordination, and communication and dissemination of ob- lowed by photons at other electromagnetic (EM) wavelengths, serving plans. these discoveries provide new insights into the physics of the Universe. While at this time the 2020 Astrophysics Decadal • A modern alert system for NASA missions and other cur- Survey report is still to be released, it is expected that strong rent Gamma-ray Coordinates Network/Transient Astron- recommendations will be made for MMA. omy Network (GCN/TAN1) streams, built on technologies The advent of advanced ground-based observatories in a few embraced by the MMA community (e.g., Apache Kafka). years will expand the discovery horizon and drastically increase the number of sources needing prompt EM follow-up from the • A Guest Observer Facility (GOF)-like service for com- ground and in space. The needs of the MMA community will munity support in observing and proposing to NASA and increase many-fold. This includes the need for coordination, ground based facilities for MMA targets. collaboration, and communication (the 3Cs) between space and • A suite of tools for analysis and interpretation of data from ground-based facilities; the need for adequate infrastructure— NASA missions. data analysis and interpretation tools, efficient alert systems, proposer and observer support, rapid data transmission links, • A curated archive of data serving the specific needs of arXiv:2109.10841v1 [astro-ph.IM] 22 Sep 2021 etc.; and the need for common and frequent transfer of ideas the MMA community and the development of relevant between communities to anticipate future needs and provide so- analysis tools, including automation of certain functions lutions. and analysis through artificial intelligence/machine learn- To fill these critical needs, NASA’s Goddard Space Flight ing (AI/ML). Center (GSFC) and Marshall Space Flight Center (MSFC) are jointly proposing to establish a virtual MMA Science Support • A service where new and existing MMA community mem- bers can obtain consultation, tools, and expertise to im- prove MMA cross-integration and to design new missions. ∗Corresponding author Email addresses: rita.m.sambruna@nasa.gov (Rita M. Sambruna), [email protected] (Joshua E. Schlieder), [email protected] (Daniel Kocevski) 1https://gcn.gsfc.nasa.gov/

Preprint submitted to Astronomy and Computing September 23, 2021 Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

A Snowball in Hell: The Potential Steam Atmosphere of TOI-1266c

C. E. Harman ,1 Ravi Kumar Kopparapu ,2 Guðmundur Stefánsson ,3, 4 Andrea S.J. Lin ,5, 6 Suvrath Mahadevan ,5, 6 Christina Hedges ,7, 8 and Natasha E. Batalha 1

1Planetary Systems Branch, Space Science and Astrobiology Division, NASA Ames Research Center, Moffett Field, CA 94035, USA 2Planetary Environments Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 3Princeton University, Princeton, NJ 08540, USA 4Henry Norris Russell Fellow 5Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 6Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 16802, USA 7Bay Area Environmental Research Institute, Moffett Field, CA 94035, USA 8Astrophysics Branch, Space Science and Astrobiology Division, NASA Ames Research Center, Moffett Field, CA 94035, USA

(Received June 18, 2021; Revised September 13, 2021) Submitted to PSJ

ABSTRACT TOI-1266c is a recently discovered super-Venus in the radius valley orbiting an early M dwarf. However, its notional bulk density (∼2.2 g cm−3) is consistent with a large volatile fraction, suggesting that it might have volatile reservoirs that have survived billions of years at more than twice the Earth’s insolation. On the other hand, the upper mass limit paints a picture of a cool super Mercury dominated by >50% iron core (∼9.2 g cm−3) that has tiptoed up to the collisional stripping limit and into the radius gap. Here, we examine several hypothetical states for TOI-1266c using a combination of new and updated open-source atmospheric escape, radiative-convective, and photochemical models. We find that water-rich atmospheres with trace amounts of H2 and CO2 are potentially detectable (SNR >∼ 5) in less than 20 hours of JWST observing time. We also find that water vapor spectral features are not substantially impacted by the presence of high-altitude water or ice clouds due the presence of a significant amount of water above the cloud-deck, although further work with self-consistent cloud models is needed. Regardless of its mass, however, TOI-1266c represents a unique proving ground for several hypotheses related to the evolution of sub-Neptunes and Venus-like worlds, particularly those near the radius valley.

Keywords: Exoplanet atmospheres (487), Hot Neptunes (754), Extrasolar ice giants (2024), Extrasolar rocky planets (511), Exoplanets (498), Star-planet interactions (2177), Super Earths (1655), Exoplanet atmospheric composition (2021)

1. INTRODUCTION Bakos et al. 2004; Hellier et al. 2012), TRAPPIST (Je- arXiv:2109.10838v1 [astro-ph.EP] 22 Sep 2021 The list of exoplanet detections from space-based ob- hin et al. 2011; Gillon et al. 2017; Delrez et al. 2018), servatories like the Kepler Space Telescope (e.g., Borucki and the Habitable-zone Planet Finder (HPF) Spectro- et al. 2010; Twicken et al. 2016) and the Transiting graph (Mahadevan et al. 2012, 2014), is rapidly grow- Exoplanet Survey Satellite (TESS; Ricker et al. 2014; ing. These detections enhance our understanding of Barclay et al. 2018), as well as ground-based endeavors planetary occurrence rates (Batalha 2014) as well as en- such as WASP (e.g. Pollacco et al. 2006), HATNet (e.g. able robust statistical insights into planet populations (e.g., Dressing & Charbonneau 2013; Burke et al. 2015; Fulton et al. 2017; Hardegree-Ullman et al. 2019). In Corresponding author: C. E. Harman particular, the presence of a gap in the radius distribu- [email protected] tion of planets (Rogers 2015; Fulton et al. 2017; Fulton & Petigura 2018) highlights the cumulative effects of a MNRAS 000,1–11 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

SCSS-Net: Solar Corona Structures Segmentation by Deep Learning

Šimon Mackovjak1?, Martin Harman2, Viera Maslej-Krešnákovᡠ2, and Peter Butka2 1Department of Space Physics, Institute of Experimental Physics, Slovak Academy of Sciences, Košice, Slovakia 2Department of Cybernetics and Artificial Intelligence, Faculty of Electrical Engineering and Informatics, Technical University of Košice, Košice, Slovakia

Accepted 2021 September 03. Received 2021 August 11; in original form 2021 May 14

ABSTRACT Structures in the solar corona are the main drivers of space weather processes that might directly or indirectly affect the Earth. Thanks to the most recent space-based solar observatories, with capabilities to acquire high-resolution images continuously, the structures in the solar corona can be monitored over the years with a time resolution of minutes. For this purpose, we have developed a method for automatic segmentation of solar corona structures observed in EUV spectrum that is based on a deep learning approach utilizing Convolutional Neural Networks. The available input datasets have been examined together with our own dataset based on the manual annotation of the target structures. Indeed, the input dataset is the main limitation of the developed model’s performance. Our SCSS-Net model provides results for coronal holes and active regions that could be compared with other generally used methods for automatic segmentation. Even more, it provides a universal procedure to identify structures in the solar corona with the help of the transfer learning technique. The outputs of the model can be then used for further statistical studies of connections between solar activity and the influence of space weather on Earth. Key words: Sun: corona – methods: data analysis – techniques: image processing – software: development

1 INTRODUCTION rithms were introduced. The Spatial Possibilistic Clustering Algo- rithm (SPoCA) (Barra et al. 2009; Verbeeck et al. 2014) or Coro- Solar activity has been quantified using various indices such as the nal Hole Identification via Multi-thermal Emission Recognition Al- sunspot number (Clette et al. 2014), the F10.7 index (Tapping 2013), gorithm (CHIMERA) (Garton et al. 2018) have been found to be coronal index (Rybanský et al. 2005), and others (Ermolli et al. very effective and are widely used in online solar data visualization 2014) over more than six decades. These indices represent inte- tools1, 2. The SPoCA also provides entries for catalogues of coro- grated quantities for specific processes in the solar atmosphere as nal holes and active regions within the Heliophysics Events Knowl- they are measured by ground-based instruments. Thanks to the most edgebase (HEK) (Hurlburt et al. 2012) that is commonly used in the recent space missions such as e.g. SOHO (Domingo et al. 1995), SolarSoft (Freeland & Handy 1998) and SunPy (The SunPy Com- Hinode (Kosugi et al. 2007), SDO (Pesnell et al. 2012), Solar Or- munity et al. 2020) frameworks. As will be presented later, these al- biter (Müller et al. 2013), and Parker Solar Probe (Fox et al. 2016), gorithms still have limitations for the precise segmentation of struc- solar activity can be monitored in a much more detailed way. The tures in the solar corona. Due to the advances in computer vision temporal resolution in terms of seconds and spatial resolution down in recent years, approaches based on machine learning techniques to hundreds of kilometers in multiple spectral bands enable individ- are able to extend the standard methods (Aschwanden 2010). Con- ual events to be studied with particular consequences for Sun–Earth ventional machine learning techniques as Support Vector Machine relations (Müller et al. 2020; Lörincíkˇ et al. 2021a). However, the (SVM), Decision Tree, or Random Forest could improve the de- challenging task is to effectively process the huge amount of images arXiv:2109.10834v1 [astro-ph.SR] 22 Sep 2021 tection of coronal holes as they provide automated distinguishing that are automatically acquired with a very high resolution. This task from filaments in EUV solar images (Reiss et al. 2015; Delouille requires sophisticated algorithms which are not as straightforward as et al. 2018). On the assumption that the techniques based on Con- the methods used for the creation of standard solar indices’ time se- volutional Neural Networks (CNN) are the most convenient tech- ries. niques for image segmentation tasks (Lecun et al. 2015), a new The initially used automated methods for detection of solar ac- era in automated processing of solar images has begun. Illarionov tivity in the solar images are based on the predefined rules in- & Tlatov(2018) demonstrated that for the segmentation of coro- ferred from the appearance and usual characteristics of the struc- nal holes, CNN provides quantitatively comparable results as the tures in the solar atmosphere (Henney & Harvey 2005; Krista & SPoCA and CHIMERA algorithms. Jarolim et al.(2021) provided a Gallagher 2009; Pérez-Suárez et al. 2011). However, it is not possi- method called CHRONOS based on CNN with progressively grow- ble to capture all the unique solar structures using generally defined ing network approach for robust segmentation of coronal holes. rules. Therefore, methods with more advanced mathematical algo-

1 https://helioviewer.org ? E-mail: [email protected] 2 https://www.solarmonitor.org

© 2021 The Authors MNRAS 000,1–10 (xxxx) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Charting nearby stellar systems: The intensity of Galactic cosmic rays for a sample of solar-type stars

D. Rodgers-Lee1?, A. A. Vidotto1 and A. L. Mesquita1, 1 School of Physics, Trinity College Dublin, University of Dublin, College Green, Dublin 2, D02 PN40, Ireland

Accepted 2021 September 22. Received 2021 August 30; in original form 2021 June 14.

ABSTRACT Cosmic rays can penetrate planetary atmospheres driving the formation of prebiotic molecules, which are important for the origin of life. We calculate the Galactic cosmic ray fluxes in the habitable zone of five nearby, well-studied solar-type stars and at the orbits of 2 known exoplanets. We model the propagation of Galactic cosmic rays through the stellar winds using a combined 1.5D stellar wind and 1D cosmic ray transport model. We find that the habitable zone of 61 Cyg A has comparable Galactic cosmic ray fluxes to present-day Earth values. For the other four systems ( Eri,  Ind, ξ Boo B and π1 UMa), the fluxes are orders of magnitude smaller than Earth values. Thus, it is unlikely that any as-of-yet undetected Earth-like planets in their habitable zones would receive a higher radiation dose than is received on Earth.  Ind b, a Jupiter-like planet orbiting at ∼11 au, receives higher Galactic cosmic ray fluxes than Earth. We find the suppression of Galactic cosmic rays is influenced by whether diffusion or advection dominates at GeV energies and at distances where the wind has reached its’ terminal velocity. For advectively-dominated winds (∼younger systems), varying the astrospheric size influences the suppression significantly. For diffusion-dominated systems (∼older systems) the astrospheric size, and therefore knowledge of the ISM properties, are not very important. This reduces the Galactic cosmic ray flux uncertainties in the habitable zone for diffusion-dominated systems. Whether a system is advection- or diffusion-dominated can be determined from the stellar wind properties. Key words: diffusion – cosmic rays – methods: numerical – stars: solar-type – stars: winds, outflows – planetary systems

1 INTRODUCTION JWST, it is important to model Galactic cosmic ray propa- gation for different stellar systems. Galactic cosmic rays ionise and dissociate molecules as they The flux of Galactic cosmic rays reaching an exoplanet’s penetrate exoplanetary atmospheres. This can lead to the orbit depends on two factors: the properties of the stellar production of prebiotic molecules (Dartnell 2011). It has also wind and the Galactic cosmic ray spectrum in the interstel- recently been suggested that cosmic rays left an imprint on lar medium (ISM). For individual stellar systems knowledge the helicity of DNA via their secondary particles, which are of the local ISM properties, combined with Lyman-α obser- on-average spin polarised (Globus & Blandford 2020; Globus vations, can be used to constrain the star’s mass-loss rate et al. 2021). Thus, Galactic cosmic rays may have been im- (e.g. Wood 2004; Wood et al. 2014, 2021). Additionally, the portant for the origin of life on Earth and on exoplanets. A ISM properties surrounding a stellar system, combined with arXiv:2109.10831v1 [astro-ph.SR] 22 Sep 2021 number of fingerprint ions, indicative of ionisation by ener- the stellar wind properties, determine the size of the astro- getic particles (i.e. cosmic rays and stellar energetic parti- sphere. γ−ray observations indicate that the Galactic cosmic cles) in exoplanetary atmospheres, have recently been iden- ray spectrum in the ISM is similar throughout the Galaxy tified (such as H O+, Helling & Rimmer 2019; Barth et al. 3 (with the exception of the Galactic centre, Aharonian et al. 2021). Emission from fingerprint ions may be detectable using 2020). Thus, the local interstellar spectrum (LIS) of Galac- transmission spectroscopy with the James Webb Space Tele- tic cosmic rays (constrained by Voyager 1 & 2 observations scope (JWST, Gardner et al. 2006). This will place limits on outside of the heliosphere, Cummings et al. 2016; Stone et al. the chemical abundances of fingerprint ions, and thus, on the 2019) is often adopted as the spectrum for other stellar sys- flux of energetic particles in exoplanetary atmospheres. In or- tems. der to perform chemical modelling of the impact of Galactic cosmic rays on exoplanetary atmospheres (e.g., Grießmeier The LIS is suppressed in an energy-dependent way as the et al. 2015; Barth et al. 2021) that will be characterised with cosmic rays travel through a stellar wind. This was first stud- ied in the context of the solar system (e.g. Parker 1965). The suppression of Galactic cosmic rays varies spatially within the ? E-mail: [email protected] solar system and temporally with the magnetic solar cycle

© xxxx The Authors MNRAS 000,1–17 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

AGN and Star Formation at Cosmic Noon: Comparison of Data to Theoretical Models

Jonathan Florez,1★ Shardha Jogee,1 Yuchen Guo,1 Sofía A. Cora,2,3 Rainer Weinberger,4 Romeel Davé,5 Lars Hernquist,4 Mark Vogelsberger,6 Robin Ciardullo,7,8 Steven L. Finkelstein,1 Caryl Gronwall,7,8 Lalitwadee Kawinwanichakij,9 Gene C.K. Leung,1 Stephanie LaMassa,10 Casey Papovich,11 Matthew L. Stevans,1 and Isak Wold12 1Department of Astronomy, University of Texas at Austin, Austin, TX 78712, USA 2Instituto de Astrofísica de La Plata (CCT La Plata, CONICET, UNLP), Observatorio Astronómico, Paseo del Bosque, B1900FWA La Plata, Argentina 3Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Observatorio Astronómico, Paseo del Bosque, B1900FWA La Plata, Argentina 4Center for Astrophysics – Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 5Institute for Astronomy, Royal Observatory, University of Edinburgh, Edinburgh EH9 3HJ, UK 6Department of Physics, Kavli Institute for Astrophysics and Space Research, MIT, Cambridge, MA 02139, USA 7Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 8The Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 9Kavli Institute for the Physics and Mathematics of the Universe, The University of Tokyo, Kashiwa, Japan 277-8583 (Kavli IPMU, WPI) 10Space Telescope Science Institute, 3700 San Martin Dr, Baltimore, MD 21218, USA 11Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA 12NASA Goddard Space Flight Center, Greenbelt, MD 20771

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT In theoretical models of galaxy evolution, AGN and star formation (SF) activity are closely linked and AGN feedback is routinely invoked to regulate galaxy growth. In order to constrain such models, we compare the hydrodynamical simulations IllustrisTNG and SIMBA, and the semi-analytical model SAG to the empirical results on AGN and SF at cosmic noon (0.75 < 푧 < 2.25) reported in Florez et al.(2020). The empirical results are based on a large mass-complete sample drawn from 93,307 galaxies with and without high X-ray luminosity 44 −1 2 3 AGN (퐿푋 & 10 erg s ), selected from a 11.8 deg area (∼ 0.18 Gpc comoving volume at 푧 = 0.75 − 2.25). The main results of our comparisons are: (i) SAG and IllustrisTNG both qualitatively reproduce the empirical result that galaxies with high X-ray luminosity AGN have higher mean SFR, at a given stellar mass, than galaxies without such AGN. SAG, however, strongly over-produces the number density of high X-ray luminosity AGN by a factor of 10 to 100, while IllustrisTNG shows a lack of high X-ray luminosity AGN at high stellar mass 11 (푀∗ > 10 푀 ) at 푧 ∼ 2. (ii) In SIMBA, the mean SFR of galaxies with high X-ray luminosity AGN is lower than the SFR of galaxies without such AGN. Contrary to the data, many high X-ray luminosity AGN in SIMBA have quenched SF, suggesting that AGN feedback, or other feedback modes in galaxies with such AGN, might be too efficient in SIMBA. arXiv:2109.10828v1 [astro-ph.GA] 22 Sep 2021 Key words: galaxies: evolution – galaxies: star formation – galaxies: general – quasars: general

1 INTRODUCTION lution. Observations have shown that the central BH mass of a galaxy correlates with the galaxy bulge stellar mass (Magorrian Understanding the connection between galaxy properties and their et al. 1998; McLure & Dunlop 2002) and galaxy bulge velocity central massive black holes (BHs) currently remains one of the dispersion (Ferrarese & Merritt 2000; Gebhardt et al. 2000; Kor- biggest challenges to formulating theoretical models of galaxy evo- mendy & Ho 2013), suggesting a clear link between galaxy and BH growth. Furthermore, the cosmic black hole accretion rate (BHAR) density and star formation rate (SFR) density are observed to peak ★ E-mail: jfl[email protected]

© 2021 The Authors Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

Mapping the 3D Kinematical Structure of the Gas Disk of HD 169142

Haochuan Yu ,1, 2 Richard Teague ,3 Jaehan Bae ,4, 5, ∗ and Karin Oberg¨ 3

1Department of Astronomy, Beijing Normal University, Beijing 100875, People’s Republic of China 2Department of Physics, University of Oxford, Oxford OX1 3RH, UK 3Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 4Department of Astronomy, University of Florida, Gainesville, FL 32611, USA 5Earth and Planets Laboratory, Carnegie Institution for Science, 5241 Broad Branch Road NW, Washington, DC 20015, USA

ABSTRACT The disk around HD 169142 has been suggested to host multiple embedded planets due to the range of structures observed in the dust distributions. We analyze archival ALMA observations of 12CO (2 − 1), 13CO (2 − 1), and C18O (2 − 1) to search for large-scale kinematic structures associated with other embedded planets in the outer disk. At 125 au, we identify a coherent flow from the disk surface to the midplane, traced by all three CO isotopologues, and interpret it as a meridional flow, potentially driven by an embedded planet. We use changes in the rotation speed of the gas to characterize the physical structure across this region, finding that at 125 au the CO emission traces regions of increased gas pressure, despite being at a surface density minimum. Developing a simple analytical model, we demonstrate that the physical structure of the gap can have non-trivial responses to changes in the surface density, consistent both with previous thermo-chemical models, and the conditions inferred observationally. Applying this technique to a range of sources will allow us to directly confront theoretical models of gap-opening in protoplanetary disks.

Keywords: kinematics and dynamics – planet–disk interactions – protoplanetary disks

1. INTRODUCTION velocity, vrot, radial velocity, vrad and vertical velocity, With the Atacama Large (sub-) Millimetre Array vz, given by (ALMA), we have the ability to characterize the physi- cal and chemical structure of protoplanetary disks, the v = v sin(i) cos(φ)+v sin(i) sin(φ)+v cos(i)+v , formation location of exoplanets. Particularly, spatially LOS rot rad z LSR (1) and spectrally resolved molecular line emission enables where φ is the azimuthal angle in the frame of reference us to study the gas structure of the disk. of the disk, i is the inclination of the disk. If we assume The dynamical structure of the disk encodes infor- that the disk is azimuthally symmetric, then we can dis- mation about any potential planet-disk interactions, as entangle the velocity structure of the disk according to embedded planets are predicted to drive large pertur- their different dependence on φ (e.g., Teague et al. 2019; bations in the gas velocity structure during their for- Disk Dynamics Collaboration et al. 2020). mation (e.g., Perez et al. 2015; P´erez et al. 2018; Pinte As protoplanetary disks are expected to be in hydro- et al. 2019; Disk Dynamics Collaboration et al. 2020; static equilibrium, the rotational velocity, vrot, is closely arXiv:2109.10822v1 [astro-ph.EP] 22 Sep 2021 Bae et al. 2021). The line-of-sight velocity, v of the LOS linked to the physical structure of the disk (e.g., Rosen- gas, is a combination of the systemic velocity, v , and LSR feld et al. 2012), through the internal motions of the disk, including the rotational v2 GM r 1 ∂P rot = star + , (2) 3/2 ∗ NASA Hubble Fellowship Program Sagan Fellow r (r2 + z2) ρgas ∂r

where P = ngaskT is the gas pressure, ρgas = µgasngas is the gas density, µgas is the mean molecular mass of the gas. As embedded planets are predicted to open large gaps in the gas surface density (e.g., Papaloizou & Lin 1984), leading to large radial changes in the gas pressure ©2021 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works. 1 DOI: 10.1109/TTHZ.2021.3095429 Terahertz Band-Pass Filters for Wideband Superconducting On-chip Filter-bank Spectrometers Alejandro Pascual Laguna, Kenichi Karatsu, David J. Thoen, Vignesh Murugesan, Bruno T. Buijtendorp, Akira Endo, and Jochem J. A. Baselmans

Abstract—A superconducting microstrip half-wavelength reso- frequency fgap ≈ 1.1 THz. The former version of DESHIMA nator is proposed as a suitable band-pass filter for broadband [5], despite showing technological readiness for astronomy at moderate spectral resolution spectroscopy for terahertz (THz) a ground-based telescope [12], had several issues in its filters. astronomy. The proposed filter geometry has a free spectral range of an octave of bandwidth without introducing spurious The reported performance of the 49 channels sampling the resonances, reaches a high coupling efficiency in the pass-band 332–377 GHz spectrum averaged to a spectral resolution of and shows very high rejection in the stop-band to minimize hRi ≈ 300 and a coupling efficiency of 8% with respect reflections and cross-talk with other filters. A spectrally sparse to the signal entering the filter-bank. Both the low coupling prototype filter-bank in the band 300–400 GHz has been de- efficiency and the low spectral resolution were caused by veloped employing these filters as well as an equivalent circuit model to anticipate systematic errors. The fabricated chip has the radiative losses in its co-planar filters and the ohmic been characterized in terms of frequency response, reporting an (conductive and dielectric) losses in the bridges balancing the average peak coupling efficiency of 27% with an average spectral potentials of the ground planes. Moreover, the filter design resolution of 940. used was incompatible with broadband operation. In order Index Terms—Astronomy, band-pass filter, DESHIMA, filter- to allow for truly wideband and highly-efficient filter-bank bank, microwave kinetic inductance detector (MKID), on-chip, spectrometers a new band-pass filter design was required. As spectrometer, superconducting, terahertz (THz), wideband we shall see later in the manuscript, co-planar technology is not suitable for THz filter-banks because it incurs in radiation problems. For this reason we opted for microstrip technology to develop the new filters. I.INTRODUCTION In this paper we propose the design of a microstrip band- EXT generation imaging spectrometers for far-infrared pass filter suitable for mid-resolution broadband THz on-chip N astronomy [1]–[4] will benefit from the development of spectrometers. We target filters with a coupling efficiency to compact on-chip terahertz (THz) filter-banks [5]–[7] with a the detectors as high as possible, an octave of instantaneous spectral resolution matched to the typical emission linewidth bandwidth free from spurious resonances and a moderate of extra-galactic sources of R ≈ 500 [1]. This advancement spectral resolution of R = 500. The rest of the article is will allow the construction of true integral field units [8], organized as follows. In SectionII, we revisit the concept of which are crucial for the tomographic mapping of the early the superconducting on-chip filter-bank spectrometer. In Sec- Universe [9]–[11], but are yet nonexistent in the far-infrared tion III we investigate the foundations of resonators as band- regime. pass filters and construct a circuit model, largely inspired by The single-pixel on-chip THz filter-bank spectrometer [13], [14], to study large filter-banks. In SectionIV, we explain DESHIMA is based upon NbTiN, a superconducting material why co-planar waveguide (CPW) technology should not be with a critical temperature of Tc ≈ 15 K, which allows used for THz filter-banks and motivate the use of microstrips circuitry with negligible conductor loss up to the the gap instead. In SectionV, a microstrip band-pass filter geometry is proposed and used to design a filter-bank sparsely sampling This work was supported by the European Research Council Consolidator arXiv:2109.10800v1 [astro-ph.IM] 22 Sep 2021 Grant (ERC-2014-CoG MOSAIC), no. 648135, and in part by the Netherlands the spectrum 300–400 GHz, whose fabrication details and Organization for Scientific Research NWO under Vidi grant no. 639.042.423. measurements are discussed in SectionVI. The concluding (Corresponding author: Alejandro Pascual Laguna.) remarks are outlined in Section VII. A. Pascual Laguna, K. Karatsu, and J. J. A. Baselmans are with SRON Netherlands Institute for Space Research, 2333 CA Leiden, The Nether- lands, and also with the THz Sensing group, Microelectronics Department, Delft University of Technology, 2628 CD Delft, The Netherlands (e-mail: II.ON-CHIP FILTER-BANK SPECTROMETER [email protected]; [email protected]; [email protected]). V. Murugesan is with SRON Netherlands Institute for Space Research, 2333 In this paper, we specifically target a fully sampled super- CA Leiden, The Netherlands (e-mail: [email protected]). conducting on-chip THz filter-bank spectrometer with a con- B. T. Buijtendorp is with the THz Sensing group, Faculty of Electrical Engi- stant spectral resolution R = 500. Such a spectrometer couples neering, Mathematics and Computer Science, Delft University of Technology, 2628 CD Delft, The Netherlands (e-mail: [email protected]). THz radiation to the chip with a broadband antenna, which D. J. Thoen and A. Endo are with the THz Sensing group, Faculty of transforms the incoming radiation onto a guided mode. From Electrical Engineering, Mathematics and Computer Science, Delft University the antenna, a superconducting transmission line, hereinafter of Technology, 2628 CD Delft, The Netherlands, and also with the Kavli Institute of NanoScience, Delft University of Technology, 2628 CJ Delft, The called through-line, guides the wideband signal with virtually Netherlands (e-mail: [email protected]; [email protected]). no loss of power [15] to the filter-bank, where it is sorted into Astronomy & Astrophysics manuscript no. Observational_errors_arxiv ©ESO 2021 September 23, 2021

Towards a systematic treatment of observational uncertainties in forward asteroseismic modelling of gravity-mode pulsators Dominic M. Bowman and Mathias Michielsen

Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium e-mail: [email protected]

Received 6 July 2021 / accepted 22 Sept 2021

ABSTRACT

Context. In asteroseismology the pulsation mode frequencies of a star are the fundamental data that are compared to theoretical pre- dictions to determine a star’s interior physics. Recent significant advances in the numerical, theoretical and statistical asteroseismic methods applied to main sequence stars with convective cores have renewed the interest in investigating the propagation of observa- tional uncertainties within a forward asteroseismic modelling framework. Aims. We aim to quantify the impact of various choices made throughout the observational aspects of extracting pulsation mode frequencies in main sequence stars with gravity modes. Methods. We use a well-studied benchmark slowly pulsating B star, KIC 7760680, to investigate the sensitivity of forward astero- seismic modelling to various sources of observational uncertainty that affect the precision of the input pulsation mode frequencies. Results. We quantify the impact of the propagation of the observational uncertainties involved in forward asteroseismic modelling. We find that one of the largest sources of uncertainty in our benchmark star is in the manual building of period spacing patterns, such that the inclusion of a potentially ambiguous pulsation mode frequency may yield differences in model parameters of up to 10% for mass and age depending on the radial order of the mode. Conclusions. We conclude that future asteroseismic studies of main sequence stars with a convective core should quantify and include observational uncertainties introduced by the light curve extraction, iterative pre-whitening and the building of period spacing patterns, as these propagate into the final modelling results. Key words. asteroseismology – stars: early-type – stars: oscillations – stars: evolution – stars: rotation – stars: fundamental parame- ters – stars: massive

1. Introduction 1993; Townsend 2005; Walczak et al. 2015), which defines the group of slowly pulsating B (SPB) stars (Waelkens 1991). Such In recent decades, asteroseismology has become an estab- a coherent pulsation excitation mechanism is predicted to pro- lished method to accurately probe the interior physics of stars duce pulsation modes with essentially infinite lifetimes from the from quantitative comparison of their observed pulsation fre- observer’s perspective. The modes are observed as periodic vari- quencies with those predicted by theoretical models of stel- ations in a star’s surface flux and can be extracted by means of lar structure and evolution (Chaplin & Miglio 2013; Hekker Fourier analysis (Aerts et al. 2010). In reality, however, the ma- & Christensen-Dalsgaard 2017; Bowman 2020; Aerts 2021). jority of coherent pulsators are not observed to have perfectly Across the Hertzsprung–Russell (HR) diagram, stars are able to periodic pulsation modes (see e.g. Bowman et al. 2016; Van excite a variety of different types of pulsation modes, which are Beeck et al. 2021), which introduces additional uncertainties in standing waves with properties defined by stellar structure. Stel- the practical analysis of a pulsating star’s light curve. lar pulsations are typically grouped into pressure (p) modes and gravity (g) modes based on their dominant restoring forces being In the last decade significant advances in our understanding the pressure force and buoyancy, respectively (Aerts et al. 2010). of main sequence B-type stars have been made because of as- However, for rotating stars, the Coriolis force is also important, teroseismology. A veritable revolution in the amount and quality which gives rise to gravito-inertial modes (Aerts et al. 2019). of time series data assembled by the CoRoT (Auvergne et al. Stars can be approximated as spheres, so stellar pulsations are 2009), Kepler/K2 (Borucki et al. 2010; Koch et al. 2010; How- arXiv:2109.10776v1 [astro-ph.SR] 22 Sep 2021 typically described using spherical harmonics. The radial order, ell et al. 2014), BRITE (Weiss et al. 2014, 2021) and now TESS n, defines the radial geometry of a pulsation mode, and the angu- (Ricker et al. 2015) space missions have provided the necessary lar degree, `, and the azimuthal order, m, together define the sur- continuous, high-precision and long-term data for detailed aster- face geometry. Identified observed pulsation mode frequencies oseismic modelling studies (Aerts 2015; Bowman 2020). Light in terms of {n, `, m} are then quantitatively compared to a grid of curves from space missions have driven increasingly more so- predicted frequencies calculated from models via forward aster- phisticated g-mode modelling methods to be developed, which oseismic modelling (Aerts 2021). have tested the impact of rotation, interior mixing, convection In main sequence late-B stars, which are born with a convec- and opacity tables on the resultant stellar pulsation frequencies tive core (i.e. 3 . M . 8 M ), an opacity mechanism operating (Moravveji et al. 2015, 2016; Buysschaert et al. 2018; Szewczuk in the Z bump is efficient in exciting high radial order (asymp- & Daszynska-Daszkiewicz´ 2018; Wu & Li 2019; Mombarg et al. totic) g modes with periods of order days (Dziembowski et al. 2019; Szewczuk et al. 2021; Pedersen et al. 2021; Michielsen

Article number, page 1 of 21 arXiv:2109.10759v1 [astro-ph.IM] 22 Sep 2021 oee,tersaeas rsnscalne o data for challenges presents also scale their astronomers; mining however, data for resource rich exabyte a undoubt- provide and will edly petabyte surveys These Ar- the respectively. in Kilometre regimes, data scale Square collect will the (SKA) or Ru- ray C. (LSST) Vera next the Observatory The as bin such sky. telescopes the survey of of regions generation sur- large col- scale observe large data which to of veys, objects method specific targeting main from a the lection change enabled to This shift paradigm dramatically. astronomers ca- and of collecting data detectors pabilities the improved have telescopes, analysis computer in advances Technological Introduction 1 model. and provenance tool associated the the both designing challenges, in involved requirements, opportunities the and discuss to are subsequently basis is which the analysis as cases used This use astronomers. ulterior to repro- interest other and of as trust astronomi- well of the as issues to the ducibility addresses specific which analysis, needs case cal use tool a generation present therefor provenance We pipelines. automated astronomical/data-processing for an making for data results. final scientific its key its system, ultimately is and the processing products, within data confidence time establish real to this pipelines automated of complex, provenance anal- of the and use stored the be Through solu- can volumes novel ysed. data demands these that how 600 pro- on world of data tions the civil order of biggest the machine the of be ducing products will data SKA and the TB/s Pbytes/year, streams 10 data to up time of real Rubin Producing C. Vera (LSST). the or Observatory (SKA) Array Kilometre as Square such observatories the generation shift next paradigm from data a handle undergoing to is astronomy decade this In Abstract h neto fti ae st a h foundation the lay to is paper this of intention The ihe .C Johnson C. A. Michael rse Lackeos Kristen 1 nttt fDt cec (DLR) Science Data of Institute srnmclPpln rvnne s aeEvaluation Case Use A Provenance: Pipeline Astronomical 2 2 1

asRie Kl¨ocknerHans-Rainer [email protected] ik Schindler Sirko acsParadies Marcus 2 e a lnkIsiuefrRdoAstronomy Radio for Institute Planck Max 1 eeaino eecps n h urn ako a of astronomi- lack the current for the tool generation and provenance suitable telescopes, next approaching of the of generation timeliness current the tronomy, annotations, users. its manual on responsibility to the which putting according implementation provenance independent language gathers a YesWork- is datasets. flow astronomical of demands scale the to meet unable to is it however scripts, Python provenance from directly generates automatically NoWorklow indepen- (language dent). [8] YesWorkflow [10] NoWorkflow and only) include (Python languages scripting for cally to and knowledge. is provenance are specific it they domain where community, require them not this meet for to en- tool important to therefore any order of In adoption languages. courage scripting native proficiency their their to within due requirement this meet re- to a respectiveluctance have the Astronomers within systems. pipelines management workflow their re- build systems to these users of quire not Implementation is pipelines practice. astronomical standard within provenance distribu- of and recording tion [1]) Chimera Kepler and (e.g., [11], Taverna pipelines [4], scientific for generation designed provenance of systems availability bene- current the the Despite and products. fits data and processing data reproducibility the and trust establishing in instrumental is Common [5]. the (CASA) and Application Observatory; Software Astronomy Rubin data C. processing Vera for the responsible from is which the [6] stack include: LSST processing data pythonic of examples astronomical, notable this Some Python. as script such within languages ing created typically are datasets these analyse the pipelines. of complexity required the size necessary of nature the the distributed also and only but not themselves from datasets stem challenges generation. provenance These and analysis, data management, h nraigiprac fpoeac ihnas- within provenance of importance increasing The specifi- tools generation provenance of Examples discipline scientific any within recording Provenance to designed are which pipelines astronomical The 1 2 1 ai .Champion J. David at Dembska Marta 2 1 in - MNRAS 000,1–9 (2020) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

The Relative Supernovae Contribution to the Chemical Enrichment History of Abell 1837

M. K. Erdim,1,2★ C. Ezer,1 O. Ünver1,3 F. Hazar1 and M. Hudaverdi1† 1Yıldız Technical University, Faculty of Science and Art, Department of Physics, Istanbul 34220, Turkey 2Yıldız Technical University, Graduate School of Natural and Applied Sciences, Istanbul 34220, Turkey 3Beykent University, Faculty of Engineering and Architecture, Istanbul 34396, Turkey

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT In this paper, we report the relative SNe contributions on the metal budget of the ICM of Abell 1837 galaxy cluster at redshift z = 0.069. For this purpose, we analysed the hot ICM of the cluster and obtained radial metal distributions by using XMM-Newton archival data with total exposure ∼100 ks. These metal measurements consist of Mg, Si, S, Fe and Ni within 0.7 R500 radius which is divided into three concentric annuli. In order to explain the observed metal abundance pattern in terms of relative supernova contributions, we used our newly developed code SNeRatio which utilizes theoretical nucleosynthesis models. This study covers the most recent 3D SNIa and SNcc yield tables. All combinations of these theoretical yields were fitted with our measured abundance ratios and statistically acceptable ones were selected. Each of these models were found to predict a uniform SNIa percentage contribution to the total SNe from the cluster center to the outskirts and form an SNIa ratio distribution with a mean 39 ± 14%. This uniformity is consistent with the early enrichment scenario which assumes that the metal production processes begin in early phase of cluster formation, namely proto-cluster phase at epoch z ≥ 2. Key words: galaxies: clusters: intracluster medium – supernovae: general – X-rays: galaxies: clusters

1 INTRODUCTION sponsible for the production and ejection of heavier metals such as Fe and Ni. Intermediate elements namely Si and S are presumably Clusters of galaxies are remarkable probes of dynamical history produced by both types. High resolution X-ray spectroscopy allows and chemical enrichment since their evolution time is comparable us to observe these metals via their emission lines which have been to the age of the Universe. These massive objects can be portrayed synthesised by billions of SNe inside the cluster galaxies. A com- by stars (∼3%) and X-ray emitting hot dilute gas (∼12%) trapped in parison of data and theoretical models, i.e SNe yield tables, is a well deep dark matter (∼85%) potential well (for a review, see Pratt et al. known method to obtain the fraction of each SNe types contribut- (2019)). Therefore, most of the baryonic matter resides in the form ing to the enrichment. This fraction provides important information of hot gas constituting the so-called intra-cluster medium (ICM). about the chemical evolution of the Universe. Observations reveal that the ICM is enriched in terms of metal abundance. Supernovae (SNe) explosions are regarded as arXiv:2109.10738v1 [astro-ph.CO] 22 Sep 2021 main mechanisms to be responsible for creating this metal-rich en- vironment. The main proposed transport mechanisms responsible for spreading these produced elements into the inter-galactic space After the launch of ASCA, the metal abundances in the ICM are ram-pressure stripping (Gunn & Gott 1972), galactic winds were studied extensively (e.g. Mushotzky et al. 1996). The spa- (De Young 1978) and Active Galactic Nuclei (AGN) outflows (e.g. tially resolved metal measurements (e.g. Si, S, Fe) from center to Blanton et al. 2001; McNamara et al. 2001; Fabian et al. 2006; larger radii for relaxed clusters of galaxies pointed toward an in- McNamara & Nulsen 2007). Core-collapse supernovae (SNcc) ex- creasing SNcc dominance at clusters’ outskirts (Finoguenov et al. plosions account for both creation and dispersion of the 훼 elements 2000, 2001). These first results supported a scenario in which SNcc such as O, Ne and Mg; whereas Type Ia supernovae (SNIa) are re- products are mostly produced before clustering and well mixed throughout the ICM, but SNIa products increase as a result of a late contribution in the central region after cluster collapse. The ★ E-mail: [email protected] (MKE) studies with better spectral resolution also verified similar enrich- † E-mail:[email protected] (MH) ment processes, e.g. XMM-Newton observations of central regions

© 2020 The Authors AUTOCORRELATION FUNCTIONS, COSMOLOGY AND INVESTIGATING THE CMB COLD SPOT WITH EMU-ASKAP RADIO CONTINUUM SURVEY Syed Faisal ur Rahman1,2

(1)University of Karachi, Karachi, [email protected] (2)Institute of Business Management, Karachi

ABSTRACT Galaxy angular-power spectrum and autocorrelation functions (ACFs) provide information about the distribution of matter by using galaxy counts as a proxy. In this study, we are going to estimate autocorrelation angular power spectrum and angular autocorrelation function for EMU-ASKAP 5 sigma sources and then compare them with results from NVSS. We will also use SUMSS data to compare ACF results using Landy-Szalay estimator. EMU-ASKAP will provide excellent opportunity to observe universe with high sensitivity and is likely going to observe millions of high redshift sources which will help in studying the clustering of the large scale structures, constraining cosmological parameters and exploring mysteries like the existence of a cosmic cold spot or the CMB cold spot as observed by both Planck and WMAP probes. We will discuss some possible ways, the CMB cold spot puzzle can be explored further by using the galaxy clustering, integral source count and galaxy bias analysis with a highly sensitive survey like EMU-ASKAP.

Keywords: cosmology, radio astronomy, EMU-ASKAP, CMB cold spot, galaxy surveys INTRODUCTION

Galaxy angular-power spectrum provides information about the distribution of matter by using galaxy counts as a proxy [1][2]. In this study, we are going to estimate autocorrelation angular power spectrum and angular autocorrelation function (ACF) for EMU 5 sigma sources [3][4][5][6] and then compare them with results from NVSS [1] [7]. We will also use SUMSS data to compare ACF results. ANGULAR POWER SPECTRUM AND AUTOCORRELATION FUNCTION

We can write the angular power spectrum, Clgg as [1]: 푘푚푎푥 푑푘 2 2 Clgg= 4휋 ∫ Δ (푘){푊푙 (푘)} Equation (1) 푘푚푖푛 푘 푔 Here Δ2(푘) is the logarithmic matter power spectrum, which can be calculated as: 푘3 Δ2(k)= 푃(푘) 2휋2 In this, P(k) is the matter power spectrum. 푊푙푔(푘) is galaxy window function. For this analysis we adopted Limber’s approximation as discussed in [8][9]. This will give an expression for Clgg as: 푧∗ 푑푧 퐻(푧) 2 Clgg= 4휋 ∫ {푊푙 (푧)} 푃(푘) Equation (2) 0 푐 Χ2(푧) 푔

With k being approximated as [10]: 1 푙 + 푘 = 2 휒2(푧) With the window function defined as [9]: 2 푧∗ ( ) ( ) 푑푁 3Ω푚 퐻0 ′ Χ 푧′ − Χ 푧 ′ 푑푁 푊푙푔(푧) = 푏(푧) + (1 + 푧)Χ(푧) ∫ 푑푧 ( ) (훼(푧 ) − 1) ′ 푑푧 2푐 퐻(푧) 푧 Χ(푧′) 푑푧 Equation (3) arXiv:2109.10730v1 [astro-ph.SR] 22 Sep 2021 ★ cetdXX eevdYY noiia omZZZ form original in YYY; Received XXX. Accepted rdeu.Mroe,tems osrt sas eyuncertai very also (Table is quoted rate commonly loss estimates mass t the reach Moreover, to takes up. it time dredge the with comparison in much not is this oie al. et Hoai ( ( variability its of amplitude remark large stars, the AGB observed for most the wi of one nascent A), the (Mira Ceti of Mira generation the governing know mechanisms our improving the and of reviewing at aims article present The INTRODUCTION 1 † itneo 0 c 0mssann u While au. 1 spanning mas 10 pc, 100 of distance a npriua,ashrclysmercwind. symmetric spherically a particular, in by obtained been have results negative Nhung P.T. archetype atypical Ceti, Mira copne yadsatcmain(iaB, (Mira companion distant a by accompanied 12 ai aaeesaelse nTable in listed are parameters Basic pn yMr nteABi hrfr neti;i sa leas at G is by detected it wake uncertain; turbulent therefore the ( is by AGB testified as the years, on 30’000 A Mira Th up. by dredge third spent the experienced already have to supposed 8 on ucVe,NhaD,CuGa,H o,Vietnam Noi, Ha Giay, Cau Do, Nghia Viet, Quoc Hoang Center 18, Space National Vietnam Astrophysics, of Department atne l 2007 al. et Martin 1991 CO/ -al [email protected] E-mail: -al [email protected] E-mail: aedtce eheimi h iaAsetu,conflictin spectrum, A Mira the in technetium detected have ) 13 Oaudnertoo 12 of ratio abundance CO ( 2020 ,hratrrfre oa pprI,maue a measured I”, “paper as to referred hereafter ), ,wihtae t neato ihteIM but ISM; the with interaction its traces which ), ★ ..Hoai D.T. , ± 1 1 ,arte o au o star a for value low rather a 2, ntepeetatcew adopt we article present the In . en nietadassuming, and indirect being ) ai ie:stars lines: radio evtoso h msino h i(-)adC(-)molec CO(3-2) and SiO(5-4) thin optically the r the of possibly of and emission analysis morpho-kin extend the review, the we of of Ceti, servations complexity Mira the of unravelling (CSE) envelope of aim the With ABSTRACT tradfrabodSOln-it within line-width SiO broad a for and star aeegh,w ieeiec o ofieeto es a v gas dense a of confinement for evidence give we wavelengths, e words: Key stars. AGB oxygen-rich ver of the CSE of the models in developing stake for suited better are companion n ae h omo lmshvn eylwSOC abundanc SiO/CO low very a (8 as having rate clumps loss mass of current form the takes and asABsas ihams osrt ttelvlo e 10 few a of level the at rate loss mass arche a good a with not stars, is AGB Ceti mass Mira eject that grain mass remark dust a We of from observations. terms the results in interpretation but earlier wind an of nascent support the quadran south-western of the generation in of observed emission SiO the that † ∼ Kipper .Ta-n,P arua,PN ip ..No,adTT Th T.T. and Ngoc, N.B. Diep, P.N. Darriulat, P. Tuan-Anh, P. , a) n nw obe to known and mag), 8 ∼ 0t 0a away). au 80 to 70 ( 1992 tr:ABadps-G icmtla atr–sas indi stars: – matter circumstellar – post-AGB and AGB stars: atr tal. et Vanture ita cdm fSineadTechnology, and Science of Academy Vietnam , ;moreover ); ethird he time e ,the n, ledge dof nd able alex g t ± 4) 13 × O32 msin nareetwt bevtosa shorte at observations with agreement In emission. CO(3-2) 10 al 1. Table htw umrz reytemi eut ntepeetsect present the in results Table main in list the and briefly l summarize published we the of that irrele fraction largely important be an such to represent i expected they little have are to features known is related it away, binarity far is B, Mira Mir companion, of mech its wind the nascent of the of exploration formation an the governing on currently focus we article, history present past the and In features related Binarity 1.1 oaiga ihvlct n t ihyvral Vemissio UV variable generally highly is its B and Mira velocity accretion VLT. thin high a the at by by surrounded rotating Dwarf infrared White a the be to in considered and i Chandra, HST by X-ray br the in clear observed by been A has wavelengths. stars pai other two the many the and linking at HST later the resolved with been 1997 has in obtained was companion its and − usto period Pulsation asls rate Mass-loss 8 pcrltype Spectral Temperature M h rtiaesoigacersprto ewe iaA Mira between separation clear a showing image first The Quantity Distance ⊙ oebscparameters. basic Some yr ∼ 5a.W hwta h asls aei episodic is rate loss mass the that show We au. 15 − 1 uhsalrta sal sue.W argue We assumed. usually than smaller much , 2 h orsodn references. corresponding the 90t 20K 3200 to 2900 10 M5-9IIIe+DA iaCt,ayia archetype atypical Ceti, Mira 110 92 − 7 few A Value 3 d 333 ± M ± 0pc 10 o htocre lvnyasbefore years eleven occurred that ion pc 9 ⊙ yr sntrltdt h mechanism the to related not is t ope a-utceityat chemistry gas-dust complex y − ugsigSOmlcls in molecules, SiO outgassing s yei em fiswn:lower wind: its of terms in type 1 − lm within olume 7 vstAM ilmteob- millimetre ALMA evisit mtc ftecircumstellar the of ematics lrlns nldn new a including lines, ular M epeo Karovska & Templeton ⊙ ai.W vlaethe evaluate We ratio. e ye&Schöier & Ryde odu tal. et Woodruff ea Hony & Heras yr a Leeuwen van ai tal. et Haniff − 1 iul iaCt – Ceti Mira vidual: Skiff Reference n rfrbyno preferably and ∼ ( 0a rmthe from au 50 2014 ( ( 1995 ( 2007 ( ( 2005 ) 2004 2001 at Yet, vant. iterature ( ) 2009 ) .As A. a mpact: ) ) ) anism has n UV n idge disc ion ) 1 r r ai MNRAS 000,1–4 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

What does lie at the Milky Way centre? Insights from the S2 star orbit precession

C. R. Argüelles,1,2,3? M. F. Mestre,1,4E. A. Becerra-Vergara,2,4,5 V. Crespi,1 A. Krut,2,3 J. A. Rueda2,3,6,7 R. Ruffini2,3,6,7 1Fac. de Ciencias Astron. y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque, B1900FWA La Plata, Argentina 2ICRANet, Piazza della Repubblica 10, 65122 Pescara, Italy 3ICRA, Dip. di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy 4Instituto de Astrofísica de La Plata, UNLP & CONICET, Paseo del Bosque, B1900FWA La Plata, Argentina 5GIRG, Escuela de Física, Universidad Industrial de Santander, 680002 Bucaramanga, Colombia 6ICRANet-Ferrara, Dip. di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy 7INAF, Istituto de Astrofisica e Planetologia Spaziali, Via Fosso del Cavaliere 100, 00133 Rome, Italy

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT It has been recently demonstrated that both, a classical Schwarzschild black hole (BH), and a dense concentration of self-gravitating fermionic dark matter (DM) placed at the Galaxy centre, can explain the precise astrometric data (positions and radial velocities) of the S-stars orbiting Sgr A*. This result encompasses the 17 best resolved S-stars, and includes the test of general relativistic effects such as the gravitational redshift in the S2-star. In addition, the DM model features another remarkable result: the dense core of fermions is the central region of a continuous density distribution of DM whose diluted halo explains the Galactic rotation curve. In this Letter, we complement the above findings by analyzing in both models the relativistic periapsis precession of the S2-star orbit. While the Schwarzschild BH scenario predicts a unique prograde precession for S2, in the DM scenario it can be either retrograde or prograde, depending on the amount of DM mass enclosed within the S2 orbit, which in turn is a function of the DM fermion mass. We show that all the current and publicly available data of S2 can not discriminate between the two models, but upcoming S2 astrometry close to next apocentre passage could potentially establish if Sgr A* is governed by a classical BH or by a quantum DM system. Key words: Galaxy: centre – Galaxy: structure – dark matter – stars: kinematics and dynamics

1 INTRODUCTION North and Subaru Telescopes, using the Adaptative Optics technique. This group has recently confirmed the detection of The most effective method to explore the nature of the su- the S2-star relativistic redshift (Do et al. 2019), converging permassive compact object at the centre of our Galaxy, Sgr with the first group in an estimate of the central object mass A*, has been the tracing of the orbits of the S-cluster stars. of about 4 × 106M . This astrometric data acquisition has been performed by two leading groups, one of which started a significant progress in From this two observational campaigns, the inference on constraining the central object mass (Eckart & Genzel 1997;

arXiv:2109.10729v1 [astro-ph.GA] 22 Sep 2021 the nature of Sgr A* has been reached on the ground of Schödel et al. 2002; Gillessen & et al. 2009; Genzel & et al. novel theoretical understandings. A recent important re- 2010; Gillessen et al. 2017). More recently this group has in- sult has been obtained in Matsumoto et al.(2020), by re- corporated the GRAVITY instrument of the VLT, allowing to considering the flare emissions around SgrA*, emphasizing detect the gravitational redshift of the S2-star (Gravity Col- both, that their motion is not purely geodesic, and estab- laboration et al. 2018a); the detection of flares or hot spots lishing a limit on the spin of a putative Kerr BH mass of close to Sgr A* (Gravity Collaboration et al. 2018b); and |a| < 0.5. Soon after, in Becerra-Vergara et al.(2020), it was the detection of the relativistic precession of the S2-star or- introduced an alternative model to the classical BH in Sgr A* bit (Gravity Collaboration et al. 2020). In parallel a second by re-interpreting it as a high concentration of quantum self- group reached similar constraints to the mass of the central gravitating DM made of fermions of about 56 keV/c2 rest compact object in Sgr A* (Ghez & et al. 1998, 2005, 2008; mass. This alternative approach can explain the astromet- Boehle et al. 2016), mainly operating with the Keck, Gemini ric data of both the S2-star and the G2 object with similar accuracy than the Schwarzschild BH scenario, but without introducing a drag force on G2 which is needed in the BH ? E-mail: [email protected] (CRA) case to reconcile it with the G2 post-pericentre passage ve-

© 2021 The Authors Phenomenological Gravitational Phase Transition: Early and Late Modifications

Nima Khosravi1, ∗ and Marzieh Farhang1, † 1Department of Physics, Shahid Beheshti University, 1983969411, Tehran, Iran (Dated: September 23, 2021) In this work we generalize the idea of a gravitational phase transition or GPT at late times [1] to allow for a modified gravity scenario in the early universe as well. The original GPT was shown to simultaneously relax the H0 and σ8 tensions and Planck internal inconsistencies. However, the primary results from GPT predic- tions implied disagreement with the data from the baryonic acoustic oscillations (BAO). Here we investigate whether the generalized GPT scenario could address the Hubble tension in the presence of the BAO data as well. We find that, despite the vastly enlarged gravitational parameter space, the BAO data strongly prefer the ΛCDM paradigm with a low inferred Hubble constant, 68.03 ± 0.76 km/s/Mpc, in tension with the local H0 measurements from [2], H0 = 74.03 ± 1.42 km/s/Mpc. This failure of the generalized GPT scenario to host the P18-BAO-R19 trio is of significance and noticeably shrinks the space of possible gravitational solutions to the Hubble tension.

I. INTRODUCTION AND MOTIVATIONS In a recent work [1] the authors proposed a phenomeno- logical phase transition in the gravity sector by assuming that The recent cosmological tensions between early and late the early standard ΛCDM transits to a modified gravity model time observations in the context of the standard model of both at the background and perturbation levels. We verified that the model could resolve the CMB-lensing, low- and high- cosmology, the so-called ΛCDM, have been extensively ad- dressed by many cosmologists in the last few years. However, `’s, σ8 and H0 tensions simultaneously (at different levels). a truly convincing theoretical model for their explanation is The BAO datasets were not included in the main analysis and parameter estimation of [1]. However, the predictions of GPT still missing. The main tension is the well-known H -tension 0 for BAO observables were shown to disagree with some BAO which is the 4.2σ discrepancy between the local measure- ments of the∼ Hubble constant [2,3,4] and its inferred value data points. In this work we generalize the model to allow for early late from the cosmic microwave background (CMB) data in the transition between two general and gravity models at some intermediate redshift. We will study this model in the ΛCDM framework [5]1. There have also been many attempts presence of BAO datasets. These standard rulers provide us for H estimation from other independent observations, with 0 with independent observations at low redshifts, where a tran- higher H values [8,9, 10, 11, 12]. Beside the Hubble tension, 0 sition in the gravity sector was shown to relax the tensions there is the milder σ8 tension with its local measurements im- pling less clumpiness in matter distribution compared to the between CMB and local Hubble measurements. On one hand, CMB-based inferences [5]. Certain internal inconsistencies BAO data are expected to break the parameter degeneracies are also reported in Planck results, specifically in the CMB and improve parameter measurements. On the other hand, the enhanced space of gravitational degrees of freedom at early lensing amplitude and the low- and high-` discrepancy. If not due to systematics, and if considered statistically signif- times may have the potential to accommodate both BAO data and local values. This is the driving motivation of this pa- icant, theoretical models beyond ΛCDM are required to ex- H0 plain them. per: to explore the huge subset of gravity models effectively modelled by a transition at some redshift, with possible si- Many different models are proposed to address these cos- multaneous deviations from GR in early and late times, and mological tensions [see 13, for a comprehensive review]. The search for a possible sweet spot where CMB, BAO and local models could in general be categorized into late and early, H measurements can all co-exist in peace. based on the era of the modification. The motivation behind 0 the late time modifications is clear and there are many at- tempts [see, e.g., 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, II. THE GPT MODEL 25, 26, 27, 28, 29, 30, 31]. The early modifications try to

arXiv:2109.10725v1 [astro-ph.CO] 22 Sep 2021 reduce the sound horizon at the last scattering surface which The gravitational phase transition, or GPT model, intro- results in a higher value for H0 [32, 33, 34, 35]. There are duced in [1], suggests a transition in the gravitational theory, many other different ideas to approach the H0 tension, such as dark matter [36], primordial magnetic fields [37] and etc. phenomenologically modelled by a "tanh" function. It was Neither of these models are consistent with all the datasets, in assumed that the early phase of gravitation is governed by particular when data from the baryonic acoustic oscillations the standard Einstein-Hilbert model, referred to as GR, which (BAO) are included [38]. forms the gravity model in the ΛCDM scenario. In this work we relax this assumption and study the consequences. We assume our Universe is isotropic, homogeneous and flat, where the background evolution is given by ∗Electronic address: [email protected] †   Electronic address:m [email protected] H2(z) = H2 Ω (1 + z)4 + Ω (1 + z)3 + Ω (z) (1) 1 For some controversy on the significance of the Hubble tension see [6,7]. 0 r m Λ MNRAS 000,1–15 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Using strong lensing to understand the microJy radio emission in two radio quiet quasars at redshift 1.7

P. Hartley,1,2★ N. Jackson,2 S. Badole,2 J. P. McKean,3,4 D. Sluse,5 H. Vives-Arias6 1SKAO, Jodrell Bank, Lower Withington, Macclesfield, Cheshire, SK11 9FT, UK 2Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK 3ASTRON, Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD, Dwingeloo, the Netherlands 4Kapteyn Astronomical Institute, University of Groningen, P.O. Box 800, 9700AV Groningen, the Netherlands 5STAR Institute, Quartier Agora - Allée du six Août, 19c B-4000, Liège, Belgium 6Centro de Astrobiología (CSIC-INTA), Carretera de Ajalvir km. 4, Torrejón de Ardoz, E28850, Madrid, Spain

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT 23 −1 The radio quasar luminosity function exhibits an upturn around 퐿6 GHz = 10 W Hz that is well-modelled by a star-forming host galaxy population. This distribution leads some authors to cite star formation as the main radio emission mechanism in so-called radio-quiet quasars (RQQs). Understanding the origin of RQQ radio emission is crucial for our understanding of quasar feedback mechanisms – responsible for the regulation of star-formation in the host galaxy – and for understanding galaxy evolution as a whole. By observing RQQs that have been magnified by strong gravitational lensing, we have direct access to the RQQ population out to cosmic noon, where evidence for twin mini-jets has recently been found in a sub-µJy RQQ. Here we present radio observations of two lensed RQQs using the VLA at 5 GHz, the latest objects to be observed in a sample of quadruply-imaged RQQs above -30◦. In SDSS J1004+4112 we find strong evidence for AGN-related radio emission in the variability of the source. In PG 1115+080 we find tentative evidence for AGN-related emission, determined by comparing the radio luminosity with modelled dust components. If confirmed in the case of PG 1115+080, which lies on the radio–FIR correlation, the result would reinforce the need for caution when applying the correlation to rule out jet activity and when assuming no AGN heating of FIR- emitting dust when calculating star formation rates. Our programme so far has shown that two of the faintest radio sources ever imaged show strong evidence for AGN-dominated radio emission. Key words: gravitational lensing: strong – galaxies: individual: SDSS J1004+4112 – galaxies: individual: PG 1115+080 – galaxies: active – galaxies: star formation – quasars: general

1 INTRODUCTION and higher resolution (e.g. Claeskens et al. 2006; Spingola et al. 2019). In particular this allows studies with radio interferometers of Strong gravitational lens systems involve the multiple imaging of objects with intrinsically faint radio emission, whose detailed inves- background galaxies or quasars, typically by a foreground galaxy

arXiv:2109.10720v1 [astro-ph.GA] 22 Sep 2021 tigation would otherwise only become possible with next generation close to the line of sight (see Courbin et al. 2002; Meylan et al. 2006; telescopes such as the Square Kilometre Array (SKA, Orienti et al. Bartelmann 2010; Jackson 2013 for reviews). Since the deflection 2015). By using strong lenses to study the faint radio emission from of light is dependent on the distribution of matter, modelling of radio-quiet quasars (RQQs), we can obtain direct evidence to ad- strong lens systems can map out dark matter in the lensing galaxy, dress one of the oldest problems in radio astronomy: the radio life variously probing large (Treu & Koopmans 2002, 2003, 2004; Son- of quasars. nenfeld et al. 2012, 2013, 2015) and small (Mao & Schneider 1998; Metcalf & Madau 2001; Kochanek & Dalal 2004; Vegetti et al. 2012; Nierenberg et al. 2017) scales of dark matter substructure. Although the first quasars were discovered by identifying Moreover, the consequent magnification of the background object strong radio sources associated with highly redshifted optical point- allows us to study it at some combination of higher signal-to-noise sources (Schmidt 1963), quasars only rarely appear to exist in a “radio-loud” state, with more than 90% (Sandage 1965) displaying only weak radio emission. Understanding the physical mechanism ★ E-mail: [email protected] for this weak emission is crucial for our understanding of galaxy

© 2021 The Authors arXiv:2109.10711v1 [astro-ph.EP] 22 Sep 2021 marco orsodn uhr ac Fenucci Marco author: Corresponding hc per ooclaefra bevrpae nthe Fig. on in placed illustrated observer as an motion, for Earth, orbital oscillate particular to its appears which to according chosen been aei opsto fHwia od ( words (fragment), Hawaiian of Its composition a 27. is April survey name 2016 1 on Pan-STARRS Observatory Haleakala the at by telescope discovered been has that and unstable. dynamically influence located typically of are are sphere therefore quasi-satellites gravitational Earth of contrast the In orbits outside time. satellites, of period true a over to they Earth the that to such close planet, stay our with configuration co-orbital 1 491)Km’aeai ml at quasi-satellite, Earth small a is Kamo’oalewa (469219) yee sn L using 2021 Typeset 23, September version Draft Earth’s ore NASA/JPL Source: h oeo h akvk ffc nteln-emdnmc fa of dynamics long-term the in effect Yarkovsky the of role The [email protected] ml htKm’aeai h otsal at’ oobtlojc k Keywords: scales. object time co-orbital long Earth’s 0.5 on t stable least also compared most at but for perspective the when companion short-term is faster Earth Kamo’oalewa bit an that a remain eff imply Yarkovsky could region the still co-orbital that it Earth showed Nevertheless, simulations the Our from t of exit orbit. orbits assum propagating nominal effect by Yarkovsky the drift the of axis caused of semi-major Kamo’oalewa the magnitude wor of introduced the this and motion of In long-term estimation dynamics. the an long-term used in the changes in the role significant studying of a diameter play small its may b to effect characterized Due is configurations. dynamics horseshoe and short-term satellite Its target. mission space erErhatri 491)Km’aea(k 06H3 sa Ea an is HO3) 2016 (aka Kamo’oalewa (469219) asteroid Near-Earth quasi-satellites a o)and (of) A T 1. E X io lnt,atris niiul 491)Kamo’oalewa (469219) individual: asteroids: planets, minor INTRODUCTION twocolumn https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=3752445#cont lewa r bet naseictp of type specific a in objects are 1 t oscillate) (to . tl nAASTeX63 in style 1 eateto srnm,Fclyo Mathematics of Faculty Astronomy, of Department ac Fenucci Marco tdnsitg1,100Blrd,Serbia Belgrade, 11000 16, trg Studentski ka 1 Rcie ;Rvsd– cetd–) Accepted –; Revised –; (Received ,adi has it and ), (the), nvriyo Belgrade of University umte oAJ to Submitted mo’o ABSTRACT 1 and oa Novakovi Bojan h usqetepoaino h anbl comet main-belt the and of Kamo’oalewa exploration 311P/PANSTARRS to mission subsequent by con- sample-return designed it the a and been Technology in Space has sists of ZhangHe, Academy Kamo’oalewa; Chinese the of the one, surface under- the second and of the mass, composition the the estimate stand model, shape a develop h efracso h ai-cec xeietaiming experiment radio-science the of performances the 2019 ( in- asteroid mining for future and technologies new for redirection of candidate test and good exploration situ a is Kamo’oalewa bility, ofr h rtoe h erErhAtri Char- Asteroid Near-Earth the one, (NEACO; Observation proposed first and been the acterization already have far: sci- missions so from two attention of and lots entists, attracted has it reasons, these 2 u oispoiiyt h at n t ria sta- orbital its and Earth the to proximity its to Due u h eodtre a encagdrecently. changed been has target 133P/Elst-Pi second comet the main-belt but visit to was plan original The ent ,i ocp td ie omptesurface, the map to aimed study concept a is ), eidcsicigbtenquasi- between switching periodic a y nyaot3 ees h Yarkovsky the meters, 36 about only s atce ersnigteclones the representing particles est n ieetsraecompositions surface different ing c ´ ueygaiainlmodel. gravitational purely a o 1 c a as aoolw to Kamo’oalewa cause may ect tri 491)Kamo’oalewa (469219) steroid yi h uue u results Our future. the in My t oobtladapotential a and co-orbital rth ons a,ntol rma from only not far, so nown yteYrosyeet We effect. Yarkovsky the by ,w drse hsiseby issue this addressed we k, 2 oti ups,smltosof simulations purpose, this To . nhn mm 2018 Emami & Anthony eial tal. et Venigalla .For ). zarro, arXiv:2109.10685v1 [astro-ph.IM] 22 Sep 2021 ya order an by n h following the and density 2010 (D Mat´ern case the for stability evaluation numerical likelihood of efficient terms to in route advantages alternative Dampe Mat´ern an and popular offering the cesses, for representations model space in tion endvsdta lo ieiodeauto in evaluation likelihood allow that devised been oea-akye l 2017 al. et Foreman-Mackey L using 2021 Typeset 23, September version Draft where where tt-pc ersnaino a´r n apdSml Ha Simple Damped Mat´ern and of representation State-space hr sacrepnec ewe ttoayGusa Proces Gaussian stationary a between correspondence a is There asinpoessaenwwdl sdi h nlsso astronomic of analysis the in used widely now are processes Gaussian ; Saat¸ci 2012 x W O inlseta este aesaesaerpeettoswh representations state-space have densities spectral tional Keywords: damped the processes. for 5/2 representations and 3/2 space Mat´ern 1/2, state analytic calculate We S ( = ( ( asinpoess(P)aeue ieyi h nlsso astrono of analysis the in widely used are (GPs) processes Gaussian saWee rcs Bona oin ihvariance with motion) (Brownian process Wiener a is n ω 1 prtosvateKla le ( filter Kalman the via operations ) autdd neirı inis nvria dloIb´ Adolfo Universidad Ciencias, Ingenier´ıa y de Facultad hti ainlfnto of function rational a is that ) f p A ( T iersohsi ieeta qaingvnby given equation differential stochastic linear t E ) , X 3 a-lnkIsiu ¨rEtaersrsh hsk Gie Physik, f¨ur Extraterrestrische Max-Planck-Institut d .Equation ). d measurement RNAAS f t iesre nlss(1916) analysis series Time , . . . , d d p t Andr´es Jord´an, − p tl nAASTeX631 in style − 1 n eeecsteen.Lna asinsaesaemdl allo models space state Gaussian Linear therein). references and , f 1 1 ), seuvln oasaesaemdlwt h following the with model space state a to equivalent is equation L (0 = d d p t 2 f p , A ilnimIsiuefrAtohsc,Chile Astrophysics, for Institute Millennium 0 ,2 1, + , . . . , = a ω = 1 uaaEyheramendy, Susana i.e. ,       d d b p 1) 0 t d − − amn1960 Kalman p d d − x ⊺ 0 0 1 0 1 d t . . . a W INTRODUCTION f 1 S t p and = ( + ABSTRACT O ω + − A ( ) · · · a y . n b x ∝ . p 1 prtos where operations, ) = . − ( d + t d p 1 2 + ) b t ( W a 2 ⊺ ω − · · · · · · · · · · · · .I hsnt epeetaayi ierGusa state Gaussian linear analytic present we note this In ). ˜ e,A.Daoa a ors24,Pe˜nalol´en, 2640, Torres Chile las Diagonal a˜nez, Av. p x ) d L /q d + f t snahtas,878Grhn,Germany Garching, 85748 ssenbachstrasse, dW ipeHroi siltr(SO asinpro- Gaussian (DSHO) Oscillator Harmonic Simple d ,2 1, ( · · · . σ ω oea-aky22,proa communication). personal 2021, Foreman-Mackey . 0 1 0 0 . a w 2 c allow ich G)over (GP) s with ) . 2 ( n oansBuchner Johannes and + o hs rcse.Ti ersnaincnoffer can representation This processes. these for t ( ) − rcwl ai 2013 Davis & Brockwell b q . . . a d 1 mncOclao asinprocesses Gaussian Oscillator rmonic d q p ipehroi siltradthe and oscillator harmonic simple       t +1 ltm eis ffiin loihshave algorithms Efficient series. time al q and O +1 , W ia iesre.Gswt ra- with GPs series. time mical n ( n R , vlaino h likelihood. the of evaluation ) q stenme fosrain (see observations of number the is ihkernel with oyoil,adtepirimplied prior the and polynomials, state 3 equation k ; ( atkie Sarkka & Hartikainen · n oe spectral power and ) ieiodevalua- likelihood w (1) (2) Astronomy & Astrophysics manuscript no. BCmass ©ESO 2021 September 23, 2021

Letter to the Editor

The mass of β Pictoris c from β Pictoris b orbital motion S. Lacour1, 2, J. J. Wang3 ?, L. Rodet4, M. Nowak5, J. Shangguan6, H. Beust7, A.-M. Lagrange7, 1, R. Abuter2, A. Amorim8, 9, R. Asensio-Torres10, M. Benisty7, J.-P. Berger7, S. Blunt3, A. Boccaletti1, A. Bohn11, M.-L. Bolzer6, M. Bonnefoy7, H. Bonnet2, G. Bourdarot6, 7, W. Brandner10, F. Cantalloube12, P. Caselli6, B. Charnay1, G. Chauvin7, E. Choquet12, V. Christiaens13, Y. Clénet1, V. Coudé du Foresto1, A. Cridland11, R. Dembet2, J. Dexter14, P. T. de Zeeuw11, 6, A. Drescher6, G. Duvert7, A. Eckart15, 16, F. Eisenhauer6, F. Gao17, P. Garcia9, 18, R. Garcia Lopez19, 10, E. Gendron1, R. Genzel6, S. Gillessen6, J. H. Girard20, X. Haubois21, G. Heißel1, Th. Henning10, S. Hinkley22, S. Hippler10, M. Horrobin15, M. Houllé12, Z. Hubert7, L. Jocou7, J. Kammerer20, M. Keppler10, P. Kervella1, L. Kreidberg10, V. Lapeyrère1, J.-B. Le Bouquin7, P. Léna1, D. Lutz6, A.-L. Maire23, 10, A. Mérand2, P. Mollière10, J. D. Monnier24, D. Mouillet7, E. Nasedkin10, T. Ott6, G. P. P. L. Otten12, 25, C. Paladini21, T. Paumard1, K. Perraut7, G. Perrin1, O. Pfuhl2, E. Rickman26, L. Pueyo20, J. Rameau7, G. Rousset1, Z. Rustamkulov20, M. Samland27, T. Shimizu6, D. Sing20, J. Stadler6, T. Stolker11, O. Straub6, C. Straubmeier15, E. Sturm6, L. J. Tacconi6, E.F. van Dishoeck11, 6, A. Vigan12, F. Vincent1, S. D. von Fellenberg6, K. Ward-Duong28, F. Widmann6, E. Wieprecht6, E. Wiezorrek6, J. Woillez2, S. Yazici6, A. Young6, and the GRAVITY Collaboration (Affiliations can be found after the references)

September 23, 2021

ABSTRACT

Aims. We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. Methods. We combined previous astrometry of β Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of β Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, β Pictoris c, was also reobserved at a separation of 96 mas, confirming the previous orbital estimations. +4.53 Results. From the astrometry of planet b only, we can (i) detect the presence of β Pictoris c and (ii) constrain its mass to 10.04−3.10 MJup. If one +1.08 adds the astrometry of β Pictoris c, the mass is narrowed down to 9.15−1.06 MJup. The inclusion of radial velocity measurements does not affect the +0.75 orbital parameters significantly, but it does slightly decrease the mass estimate to 8.89−0.75 MJup. With a semimajor axis of 2.68 ± 0.02 au, a period of 1221 ± 15 days, and an eccentricity of 0.32 ± 0.02, the orbital parameters of β Pictoris c are now constrained as precisely as those of β Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets’ dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet. Key words. Exoplanets – Astrometry and celestial mechanics – Instrumentation: interferometers – Techniques: high angular resolution

1. Introduction 70 b (Wang et al. 2021), and HR8799 e (Brandt et al. 2021b). One of the main difficulties of a direct measurement is the long The formation and evolution of giant exoplanets is an intense orbital period of directly imaged exoplanets. Another is the fact field of research. Several formation scenarios, ranging from that young stars are often pulsating (for example, β Pictoris is gravitational instabilities (e.g., Boss 2011; Nayakshin 2017) to a δ Scuti variable), which makes accurate radial velocity (RV) a variety of core-accretion models (e.g., Alibert et al. 2005; Em- measurements difficult. senhuber et al. 2020), are still actively debated. The issue is the lack of observables with which to distinguish between the differ- An efficient technique for measuring dynamical masses is ent scenarios. One solution is to analyze the atmospheric com- observing multi-planetary systems and detecting the mutual in- position and search for formation signatures (Öberg et al. 2011; fluence of planets. The main historical example is the predic- arXiv:2109.10671v1 [astro-ph.EP] 22 Sep 2021 Mordasini et al. 2016; Madhusudhan et al. 2017). Another so- tion of Neptune from the irregularities in the orbit of Uranus lution consists in measuring the energy dissipation during for- by Le Verrier (1846). A more recent example is the first exo- mation as a function of the final exoplanetary mass (Mordasini planetary system detected around PSR B1257+12 (Wolszczan 2013; Marleau & Cumming 2014). & Frail 1992), where the mutual interactions of the planets were This paper focuses on the latter. Obtaining a dynamical mass used to confirm their masses (e.g., Konacki & Wolszczan 2003). for young directly imaged exoplanets is difficult. Only a handful Last, transit timing variations (Agol et al. 2005) has become a of these objects have published dynamical masses: β Pictoris b key technique for obtaining the mass of transiting planets (e.g., (Snellen & Brown 2018; Dupuy et al. 2019; Lagrange et al. Demory et al. 2020). In this paper we demonstrate that optical 2020; Vandal et al. 2020), β Pictoris c (Nowak et al. 2020), PDS interferometry is now able to detect and measure the mass of an exoplanet solely from the astrometry of another exoplanet at ? 51 Pegasi b Fellow. larger separation.

Article number, page 1 of 11 arXiv:2109.10655v1 [astro-ph.GA] 22 Sep 2021 [email protected] orsodn uhr uamdA Latif A. Muhammad author: Corresponding clmse fafwhnrdslrmse n omin form and ( masses halo solar per hundred one few isolation, a of masses ical uni- high-redshift the of and understanding verse. question our open to Their an central still is metals. is of (IMF) function absence mass the initial in cooling line vibrational sarsl fgaiainlclas nue yH by induced collapse gravitational of result a as tla lc oe.PpIIsasaeepce oform to expected are stars III 10 Pop in holes. first black the of chemi-stellar formation cosmic early the and the reionization, enrichment, enrichment ended cosmic cal Bang began and Big Ages the Dark after years million 200 h rtsuissgetdta o I tr aetyp- have stars III Pop that suggested studies first The yee sn L using 2021 Typeset 23, September version Draft h mrec fPplto I PpII tr 0 - 100 stars III) (Pop III Population of emergence The 5 − hrceitcmse f11 M 1-10 of masses characteristic > assi u oesrnefo .-7M of 0.1-37 fraction from the range and number, models spin our halo in format disk, Star masses the years. onto thousand star few accretion 23 a mass for to d formation up outflows, star with strong delays generates and time, stars over the increase from feedback numbers Radiative their t st that find III We and Pop simulations. minihalo investigate cosmological we radiation-hydrodynamical Here 3D with constraints. previous to observational t limitations any to of of because importance uncertain their remain of properties spite their in However, enrichment. chemical and Keywords: alo which, interactions masses. three-body by final disks their their from ejected are 10 ouainII(o I)saseddtecsi akAe n beg and Ages Dark cosmic the ended stars III) (Pop III Population 3 0M 10 nttt o srnm,Uiest fEibrh oa Ob Royal Edinburgh, of University Astronomy, for Institute 6 1 M hsc eatet olg fSine ntdAa Emirat Arab United Science, of College Department, Physics A ⊙ ⊙ T 1. 2 E n oto h tesae11 M 1-10 are others the of most and rmrilmnhlsat minihalos primordial nttt fCsooyadGaiain nvriyo Port of University Gravitation, and Cosmology of Institute X ehd:nmrcl—eryuies aais ihrdhf dar — high-redshift galaxies: — universe early — numerical methods: rtstars first INTRODUCTION twocolumn rm ta.2002 al. et Bromm uamdA Latif, A. Muhammad tl nAASTeX63 in style h it asFnto fPpIIStars III Pop of Function Mass Birth The ⊙ n o-ev M ihdN/dM with IMF top-heavy a and z ; ble al. et Abel 20 = 1 ailWhalen, Daniel umte oApJ to Submitted 2 − ⊙ ABSTRACT n fte5 tr htfr normdl wleare twelve models our in form that stars 55 the of and , ro- ⊙ 30 u iuain hssgetta o I tr have stars III Pop that suggest thus simulations Our . yrdnmclsmltosby simulations hydrodynamical ht3 ikprilsfr noemnhl n h most the M and 20 found to minihalo grows and one star in only massive form protostar particles sink massive 37 most that the from back for- the protostar. after the years 10 of survived third mation a only and merged 2002 o la o aito rmmlil tr ntedisk studies. the these in of outcome stars the multiple affect from would radiation how clear disk M not the 600 in to form grows to star allowed the was that star found one only which in ( 2014 ( pro- al. et can Hirano that ( clumps stars multiple multiple duce into of center fragments the halo at disk the accretion the that show simulations httermse ag rm1 00M found 1000 - and 10 minihalo from range per masses forming their stars that single only ered in(e also (see tion evtr,BakodHl,EibrhE93J UK 3HJ, EH9 Edinburgh Hill, Blackford servatory, 2016 2012 2 sUiest,P o 55,A-i,UAE Al-Ain, 15551, Box PO University, es ; n aehKhochfar Sadegh and efre omlgclsmlto ihfeed- with simulation cosmological a performed ) later fragments of thirds two almost that found ) ohd ta.2008 al. et Yoshida muh otmuhP13X UK 3FX, PO1 Portsmouth smouth, frstesronigpootla disk, protostellar surrounding the eforms o ae aywt aoaddpn on depend and halo with vary rates ion rfraini v rmrilhalos primordial five in formation ar gwt oiigU edak limits feedback, UV ionizing with ng ueia iuain n h lack the and simulations numerical eeouino h al Universe, early the of evolution he ∝ nerycsooia reionization cosmological early an asv tr ntehl.Stellar halo. the in stars massive ai clihr2015 Schleicher & Latif M a utpesasfr neach in form stars multiple hat omn noeo h halos. the of one in forming s ∗ ; − 1 lr ta.2011 al. et Clark hrae l 2020 al. et Sharda . 18 pt 0 ftestars the of 70% to Up . .Hwvr eetnumerical recent However, ). ⊙ iaoe al. et Hirano gs reionization, ages, k 3 ntefis y.Radiation kyr. 5 first the in ooaae al. et Hosokawa ⊙ ; .Hwvr tis it However, ). ycupmigra- clump by ai ta.2013 al. et Latif ). ⊙ ( 2015 . tc tal. et Stacy ri tal. et Greif consid- ) ( 2016 ) ; MNRAS 000,1–9 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Nonlinear particle reacceleration by multiple shocks

T. Vieu,1★ S. Gabici,1 V. Tatischeff2 1Université de Paris, CNRS, Astroparticule et Cosmologie, F-75013 Paris, France 2Université Paris-Saclay, CNRS/IN2P3, IJCLab, 91405 Orsay, France

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT When the pressure of particles accelerated at shock waves is no longer negligible compared to the kinetic pressure of the gas, the linear theory of diffusive shock acceleration breaks down. This is expected in particular when the shock sweeps up preexisting cosmic rays, or when multiple shocks reaccelerate successively the same particles. To describe these systems, one has to account for the nonlinear backreaction of the particles on the magnetohydrodynamic flow. Using an up-to-date semi-analytical model of particle reacceleration at nonlinear shocks, we show that the presence of prexisting energetic particles strongly affects the shock profile, in such a way that the reacceleration of non thermal particles or the acceleration of particles from the thermal bath becomes less efficient. We further describe the evolution of the distribution of particles after several shocks and study the properties of the asymptotic solution. We detail the case of identical shocks as well as more realistic scenarios, including the heating of the medium or superbubble environments. When the particles are efficiently confined in the acceleration region, it is generally found that the spectrum converges toward a concave solution after a few tens of shocks, with a spectral index around 3.5 at the highest energy. The postshock cosmic ray pressure reaches an asymptotic value of about 4 – 5% of the ram pressure of one shock. Most of the shock pressure is transferred to escaping particles. Key words: acceleration of particles – shock waves – cosmic rays

1 INTRODUCTION of the SN explosion. In these cases, an accurate DSA computation requires to consider nonlinear effects. The main sources of cosmic rays (CR) in our galaxy are believed to be supernova remnant (SNR) shocks. Particle acceleration at shocks Since the pioneer works of Eichler(1979); Blandford(1980); has been investigated for decades under the framework of diffusive Drury & Voelk(1981), there have been numerous attempts to model shock acceleration (DSA) (Axford et al. 1977; Krymskii 1977; Bell nonlinear particle acceleration at shocks, including Monte-Carlo 1978a; Blandford & Ostriker 1978). A standard simplifying assump- simulations (Ellison & Eichler 1984), two-fluid models (e.g. Malkov tion made in DSA computations is that the accelerated CRs are not & Voelk 1996), semi-analytical solutions (Blasi 2002; Caprioli et al. energetic enough to influence the fluid dynamics. This gives rise 2009), numerical resolutions (e.g. Kang & Ryu 2011), hydrody- to the theory of linear DSA, which has the remarkable property of namic simulations (e.g. Caprioli & Spitkovsky 2014). It is now predicting a universal power-law in momentum for the distribution understood that the backreaction of the accelerated particles cre- function of CR accelerated at strong shocks ( 푓 (푝) ∝ 푝−4). ates a precursor in the gas upstream of the shock. The compression factor at the subshock decreases while the total compression factor arXiv:2109.10639v1 [astro-ph.HE] 22 Sep 2021 The linear theory is however expected to break down in most astrophysical environments. The CR luminosity of the galaxy in- increases, in such a way that the CR spectrum tends to develop a deed suggests that a fraction of about 10% of the energy provided concave shape instead of the universal power law predicted by the by supernovae (SNe) is converted into CRs (Strong et al. 2010). This linear theory. requires high injection efficiencies at the shock, such that the CR All the aforementioned works assumed that no seeds preexist pressure cannot be neglected anymore. Besides, CR are in equipar- in the medium before the acceleration process, i.e., the shock is tition with the gas pressure in the interstellar medium: not only the isolated. Multiple shock acceleration is however expected to take shock accelerates fresh CRs, it also sweeps up preexisting non ther- place in several physical environments, including, e.g., stellar winds mal particles (hereafter called seeds). When SNRs expand in low (White 1985), accretion disks (Spruit 1988; Achterberg 1990), su- density media, they span a large volume, such that the total energy perbubbles (Bykov 2001; Parizot et al. 2004), or clusters of galaxies of the seeds may become of the order of 1051 erg, that is, the energy (Gabici & Blasi 2003; Kang 2021). Within the framework of linear DSA, it was already shown by Bell(1978b) that multiple shock acceleration efficiently produce hard spectra. Including adiabatic ★ E-mail: [email protected] decompression, Melrose & Pope(1993) demonstrated that the spec-

© 2021 The Authors arXiv:2109.10624v1 [astro-ph.SR] 22 Sep 2021 ulse ne raieCmoslcneC YS 4.0 BY-SA CC license Commons Creative under Published ora fteKra srnmclSociety Astronomical 00 Korean the of Journal ino h aa npriua,i a o ere- be are not users may ALMA Re- it principle, Maximal In particular, so-called (MRS). the In interpreta- number, Scale coverable single reliable a knowl- data. a to precise for ducible the the mandatory of is acceptance, tion which array of of the morphology edge of complex a map high in collected the result the producing data may configurations However, baselines with antenna short (ACA). ALMA, possible observa- array of of dish diversity compact case single the the with with using in coped data or, typically array tions, is the and SSP merging 2012, The by al. et therein). “ghosts” (Chandra the references as distributions” of known “invisible artefacts coverage produce or Walterbos incomplete to & known Braun particular, is 2018; plane In numer- al. of 1985). et subject (Faridani the been studies problem has ous short-spacing and known so-called well This is base- (SSP) away target. short flux the detected reso- of of from losses detriment angular significant in the high resulting lines, at favoured of pos- configurations availability often of the lution has choice particular, baselines broad In long a configurations. offering unprecedented sible antennas, an with of a interferometry number radio open of has era Large new Atacama (ALMA) the Array of operation Millimeter/submillimeter the years, recent the In orsodn author: Corresponding oTiHa,Pa unAh hmTytNug ireDarriula Pierre Nhung, Tuyet Pham Anh, Tuan Pham Hoai, Thi Do eateto srpyis ita ainlSaeCenter Space National Vietnam Astrophysics, of Department 1 : tr ehius mg rcsig interferometric processing, image techniques: – stars eyipratadmyesl eudrsiae,o vnoverlook even or recove underestimated, to impa be the The array easily of maxima observations. central observations may effective Array) and the as Compact important only (Atacama adopt very ACA use should or and one dish distribu array cases single baseline extended such double-humped the M In a with the in groups. as result s two such may a the array, number, by the single complemented of a configurations, centre by ante antenna measurable of extended simply whic number particular not for large form, situations very generates complex The (ALMA) target. Array the Millimeter/sub-millimeter r from a causing distances interferometers, radio angular of performance the of limitation Abstract: h ulse ieaue r reydsrbdadcmae ihRD au R 15 words: with Key within Hya. compared W Hya and to W described similar briefly of very are envelope presente literature, is circumstellar star published the the the of from away Properties observations the scope. of reliability the of ∼ 9 01xx 2021 99, ALMA ucVe,NhaD,CuGa,H o,Vietnam. Noi, Ha Giay, Cau Do, Nghia Viet, Quoc .1 I 1. 1. h ako hr aeie,rfre oa hr-pcn problem short-spacing as to referred baselines, short of lack The tr:ABadps-G icmtla atr–sas individual: stars: – matter circumstellar – post-AGB and AGB stars: BEVTOSOF OBSERVATIONS NTRODUCTION 29 ..Ha n .Tuan-Anh P. and Hoai D.T. i(-)ln msino G trWHaa nilsrto.Aciia disc critical A illustration. an as Hya W star AGB of emission line SiO(8-7) eevdAgs x 01 cetdx x 2021 xx, xx accepted 2021; xx, August Received ihNo,adTa h Thai Thi Tran and Ngoc, Bich H W uv YA 1 : ita cdm fSineadTcnlg,1,Hoang 18, Technology, and Science of Academy Vietnam , 04) t ansqec aswsbten1ad1.5 and 1 between was mass M sequence main Its 2014a). u vnLewn20) ti eirglrvariable semi-regular of a is period It a 2007). with Leeuwen (van Sun e utda ia(ezle ta.20)belong- is rate 2005) mass-loss Its al. M7.5e-M9ep. et ∼ class (Lebzelter spectral Mira to a ing as quoted ten icmtla neoe(S)o y,a oxygen- only of an distance Hya, a at W star of AGB rich (CSE) envelope circumstellar or nsuc vrtresprt lcsbten30 of between blocks total separate a three for over source out on carried hours were Hya which W Takigawa), A. of (PI: ADS/JAO.ALMA#2015.1.01446.S observations archival project from uses work present The fteeiso fthe illus- of observations We emission over- ALMA the archival to real. of using danger argument is the the distortions trate and imaging significant measure the look underestimate this to of easy often importance is it however, practice, l hpe fteAM ehia handbook technical exam- ALMA spe- for the (see the of 7 used for Chapter being process ple configuration imaging antenna pro- the to cific of encouraged are simulations and fact duce this of warned properly 1 MATO ISN BASELINES MISSING OF IMPACT https://almascience.nao.ac.jp/documents-and- tools/cycle7/alma-technical-handbook 10 toivs.r.n [email protected] [email protected], ⊙ − Kor ta.21b aioihe l 2017). al. et Danilovich 2014b; al. et (Khouri .2 O 2. 2. 7 M dcino h eetdflxa increasing at flux detected the of eduction ⊙ yr oehrwt omnso broader a of comments with together d https://doi.org/10.5303/JKAS.2021.00.0.1 algopo lsb nensa the at antennas closeby of group mall BSERVATIONS isn u,a n ol owith do would one as flux, missing r nsoeae nteAaaaLarge Atacama the in operated nnas − d h rsn td ssarchival uses study present The ed. inwt infiatgpbetween gap significant a with tion 1 eoeal cl htassociated that scale recoverable l r nABsa aigproperties having star AGB an or, to h isn aeie a be can baselines missing the of ct Mece ta.20;Kor tal. et Khouri 2008; al. et (Maercker h mato h S ae a takes SSP the of impact the h IS:1225-4614 pISSN: ∼ rmtesa,ntmnindin mentioned not star, the from 6 as(au ta.21) of- 2017), al. et (Samus days 361 xmlRcvrbeSae In Scale. Recoverable aximal ,Pa gcDe,Nguyen Diep, Ngoc Pham t, 29 IMAGING SiO( , SP,i well-known a is (SSP), AARDCINAND REDUCTION DATA y ai lines: radio – Hya W ν =0, http://jkas.kas.org J · ∼ 87 ieb the by line =8-7) IS:2288-890X eISSN: 104 − +14 11 ussion ct the to pc 1 .In ). ∼ 2 MNRAS 000,1–1 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Multi-resolution angular momentum measurements of 푧 ∼ 1.5 − 2 star-forming galaxies

Juan M. Espejo Salcedo1★, Karl Glazebrook1, Deanne B. Fisher1, Sarah M. Sweet3,4, Danail Obreschkow2, A. M. Swinbank5, Steven Gillman5,6,7 and Alfred L. Tiley2 1Centre for Astrophysics & Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, VIC 3122, Australia 2International Centre for Radio Astronomy Research, University of Western Australia, 7 fairway, Crawley, WA 6009, Australia 3School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 4ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) 5Centre for Extra-galactic Astronomy, Durham University, South Road, Durham DH1 3LE UK 6Cosmic Dawn Centre (DAWN), Copenhagen, Denmark 7DTU-Space, Technical University of Denmark, Elektrovej 327, DK-2800 Kgs. Lyngby, Denmark

Accepted 2021 September 20. Received 2021 September 6; in original form 2021 May 5

ABSTRACT

We present detailed stellar specific angular momentum ( 푗∗) measurements of ten star- forming galaxies at 푧 ∼ 1.5 − 2 using both high and low spatial resolution integral field spectroscopic data. We developed a code that simultaneously models the adaptive optics (AO) assisted observations from OSIRIS/SINFONI along with their natural seeing (NS) counterparts from KMOS at spatial resolutions of [0.1 − 0.4] arcsec and [0.6 − 1.0] arcsec respectively. The AO data reveals 2/10 systems to be mergers and for the remaining eight the mean uncertainties Δ¯ 푗∗ decrease from 49% (NS), and 26.5% (AO), to 16% in the combined analysis. These 푗∗ measurements agree within 20% with simple estimates ( 푗˜∗) calculated from the Hubble Space Telescope photometry and NS kinematics, however higher resolution kinematics are required to first identify these disks. We find that the choice of surface mass density model and the measurement of effective radius from photometry are the key sources of systematic effects in the measurement of 푗∗ between different analyses. Fitting the 푗∗ vs 푀∗ relations (Fall, 1983) with a fixed power-law slope of 훽 = 2/3, we find a zero-point consistent with prior NS results at 푧 ≥ 1 within ∼ 0.3 dex. Finally, we find a ∼ 0.38 dex scatter about that relation that remains high despite the AO data so we conclude it is intrinsic to galaxies at 푧 > 1. This compares to a scatter of ≤ 0.2 dex for disks at 푧 ' 0 pointing to a settling of the Fall relation with cosmic time. Key words: galaxies: disks – galaxies: kinematics and dynamics – galaxies: evolution

1 INTRODUCTION mentum 푗∗ = 퐽∗/푀∗ was lower than that of local disks (Law et al. 2009; Obreschkow et al. 2015; Teklu et al. 2015). Approximately half of all the stars in the Universe were formed dur-

arXiv:2109.10614v1 [astro-ph.GA] 22 Sep 2021 The morphological diversity of regular galaxies at low red- ing the time period when it was 3-6 Gyr old (1 < 푧 < 3), with the shift, typically classified by the Hubble sequence (Hubble 1926), peak of star formation happening at around 푧 ∼ 2 (Madau & Dick- emerged at 푧 ∼ 2 (Szomoru et al. 2011), where the transition inson 2014). At this epoch, the majority of the galaxy population from the high-redshift turbulent disks to the grand design spirals appears very different from the rapidly rotating disks that make up we observe in the local Universe begins to take place. Both ob- the bulk of the star-forming galaxies in the local Universe. At 푧 ∼ 2, servations and theoretical models suggest that galaxy-wide insta- star-forming galaxies had clumpy morphologies (Glazebrook et al. bilities play an important role in driving the evolution of galax- 1995; Conselice et al. 2000; Elmegreen et al. 2005; Livermore et al. ies during this epoch (Dekel et al. 2009; Dekel & Burkert 2014; 2015), their gas fractions ( 푓 ) were high compared to local ana- gas Tacchella et al. 2016; Genzel et al. 2008). These strong instabilities logues (Bouché et al. 2007; Tacconi et al. 2013; Decarli et al. 2019) can be produced in major merger events and gravitational interac- and for a fixed stellar mass (푀∗), their stellar specific angular mo- tions which are more frequent at 푧 > 1.5 (Ventou et al. 2017, 2019). In the absence of mergers, a galaxy-wide instability in a rotating sys- tem can be generated by a combination of high gas surface density ★ E-mail: [email protected] and low angular momentum.

© 2021 The Authors MNRAS 000,1–10 (2020) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

A GHz Peaked Spectrum radio galaxy born at one-tenth the current age of the Universe

Tao An 1,3★, Ailing Wang 1,2†, Yingkang Zhang1, J.N.H.S. Aditya1, Xiaoyu Hong1, Lang Cui3 1Shanghai Astronomical Observatory, Chinese Academy of Sciences, Nandan Road 80, Shanghai 200030, China 2College of Astronomy and Space Sciences, University of Chinese Academy of Sciences, 19A Yuquanlu, Beijing 100049, China 3 Xinjiang Astronomical Observatory, Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, 150 Science 1-Street, Urumqi 830011, China

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Studies of high redshift radio galaxies can shed light on the activity of active galactic nuclei (AGN) in massive elliptical galaxies, and on the assembly and evolution of galaxy clusters in the Universe. The vast majority of observed high redshift (푧 > 4.5) AGN are quasars, with very few galaxies. J1606+3124 is a radio galaxy at a redshift of 4.56, at an era of one-tenth of the current age of the Universe. Very long baseline interferometry (VLBI) images reveal a two-sided jet structure with edge-brightened terminal hotspots separated by about 68 parsecs. No evidence for the presence of extended emission (relic of past activity) is found, suggesting that it might be a nascent radio source. We study the radio properties of J1606+3124 including radio spectrum, variability, core brightness temperature, jet proper motion. All observations consistently indicate that it is a GHz Peaked Spectrum (GPS) source, making it the highest redshift GPS galaxy known to date. The expansion velocity of the hotspots and the turnover in the radio spectrum suggest that J1606+3124 is a young (kinematic age of ∼3000 years) and developing radio source. Its ultra-high jet power gives it a good chance to grow into a large-scale double-lobe radio galaxy. Infrared observations reveal a gas- and dust-rich host galaxy environment, which may hinder the jet growth. Key words: galaxies: high-redshift–galaxies: active – galaxies: jets – galaxies: nuclei– proper motions – instrumentation: high angular resolution

1 INTRODUCTION is a radio galaxy at 푧 = 4.56 (Malkin 2018). It is very faint in the optical band (Stanghellini et al. 1993) but bright in the infrared High redshift radio galaxies (HzRGs) are the progenitors of the (named WISEA J160608.53+312446.5, Cutri et al. 2013) and radio dominant elliptical galaxies in galaxy clusters in the local Universe, bands (e.g., Owen & Mufson 1977). J1606+3124 has been detected and clusters are the building blocks of the Universe (Miley & De in several radio surveys, including the Green Bank 91 m telescope Breuck 2008). The study of HzRGs sheds light on the assembly survey at 4.85 GHz (Gregory & Condon 1991; Becker et al. 1991), and evolution of large-scale structures in the Universe. Moreover, the Faint Images of the Radio Sky at Twenty Centimeters (FIRST, HzRGs are also among the most massive galaxies at their redshifts Becker et al. 1995) of the Very Large Array (VLA) and NRAO VLA (Overzier et al. 2009). They are therefore fascinating objects on their arXiv:2109.10599v1 [astro-ph.GA] 22 Sep 2021 Sky Survey at 1.4 GHz (NVSS, Condon et al. 1998), the VLBA own to study the activity of the active galactic nuclei (AGN) and the Calibrator Survey (VCS, Beasley et al. 2002) and the Combined accretion of supermassive black holes (SMBHs). However, the chal- Radio All-Sky Targeted Eight GHz Survey of the VLA at 8.4 GHz lenges of observing HzRGs are significant, resulting in only a few (CRATES, Healey et al. 2007). It has also been detected in the have been observed at 푧 4.5, when the Universe was about one- & ongoing surveys, such as the Very Large Array Sky Survey at 3 GHz tenth its current age (Saxena et al. 2018). Despite the observational (VLASS, Gordon et al. 2020) and the Australian Square Kilometre difficulties, an in-depth study of the few individual HzRGs that have Array Pathfinder (ASKAP) continuum survey at 0.888 GHz (the been identified and confirmed is well worthwhile and necessary due first large-area survey named RACS, McConnell et al. 2020). to the critical value of the HzRGs in astrophysical research. The source FIRST J160613.1+312623 (in short, J1606+3124) There was once a controversy in the radio classification of J1606+3124. In early observations made three decades ago, this source was identified as a gigahertz-peaked spectrum (GPS) radio ★ E-mail: [email protected] source (Stanghellini et al. 1993). However, some other subsequent † E-mail: [email protected] studies identified it as a flat spectrum source (Torniainen et al.

© 2020 The Authors MNRAS 000,1–13 (2020) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Chemical abundances in Seyfert galaxies – VIII. Argon abundance estimates

A. F. Monteiro,1,2? O. L. Dors,1 1Universidade do Vale do Paraíba. Av. Shishima Hifumi, 2911, CEP: 12244-000, São José dos Campos, SP, Brazil 2Instituto Federal do Maranhão. Av. Newton Bello s/n, CEP: 65906-335, Imperatriz, MA, Brazil

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT

For the first time, the argon abundance relative to hydrogen abundance (Ar/H) in the narrow line region of a sample of Seyfert 2 nuclei has been derived. In view of this, optical narrow emission line intensities of a sample of 64 local Seyfert 2 nuclei (z < 0.25) taken from Sloan Digital Sky Survey DR7 and measured by the MPA/JHU group were considered. We adopted the Te-method for AGNs, which is based on direct determination of the electron temperature, together with a grid of photoionization model results, built with the Cloudy code, to obtain a method for the derivation of the < < Ar/H abundance. We find that for a metallicity range of 0.2 ∼ (Z/Z ) ∼ 2.0, Seyfert 2 nuclei present Ar/H abundance ranging from ∼ 0.1 to ∼ 3 times the argon solar value, adopting log(O/H) = −3.31 and log(Ar/H) = −5.60. These range of values < < < < correspond to 8.0 ∼ (12+log(O/H) ∼ 9.0 and 5.4 ∼ (12+log(Ar/H) ∼ 6.9, respectively. The range of Ar/H and Ar/O abundance values obtained from our sample are in consonance with estimations from extrapolations of the radial abundance gradients to the central parts of the disk for four spiral galaxies. We combined our abundance results with estimates obtained from a sample of H ii galaxies, which were taken from the literature, and found that the Ar/O abundance ratio decreases slightly as the O/H abundance increases. Key words: galaxies: Seyfert – galaxies: active – galaxies: abundances –ISM: abun- dances –galaxies: evolution –galaxies: nuclei

1 INTRODUCTION drogen (O/H, e.g. McGaugh 1991; Yates et al. 2012; Kewley et al. 2019). Other elements such as the noble gases (e.g. Active Galactic Nuclei (AGNs) present strong metal Ne, Ar) present emission-lines in the optical spectrum emit- emission-lines in their optical spectra, which when combined ted by only few of their ions (e.g. Ne2+, Ar2+), which make with hydrogen recombination lines, make it possible to esti- it necessary to apply Ionization Correction Factors (ICFs) mate the abundance of heavy elements and the metallicity in proposed by Peimbert & Costero(1969) (see Stasińska 2002; the gas phase of these objects. AGNs play an essential role Delgado-Inglada et al. 2014; Dors et al. 2013, 2016) in order in chemical abundance studies of nearby objects and of the to account for the unobserved ions in the estimation of the arXiv:2109.10590v1 [astro-ph.GA] 22 Sep 2021 early stages of galaxy formation due to the aforementioned total abundance. The use of ICFs can introduce uncertain- feature and to their high luminosity. ties in order of 20% (e.g. Henry et al. 1996; Alexander & Among the heavy elements, oxygen presents strong Balick 1997; Croxall et al. 2016) in the resulting total abun- emission lines (i.e. [O ]λ3726, λ3729; [O ]λ5007) emitted ii iii dance (for a detailed discussion on ICF uncertainties see by its most abundant ions (O+, O2+) in the optical spec- Delgado-Inglada et al. 2014). Despite this drawback, noble trum of gaseous nebulae (H regions, Planetary Nebulae) ii gases can be used to derive the metallicity with some advan- and AGNs (e.g. Koski 1978; van Zee et al. 1998; Kennicutt tages over oxygen, as they are also useful elements for deter- et al. 2003; Maciel et al. 2007; Dopita et al. 2015; Flury & mining constraints in stellar nucleosynthesis studies. Noble Moran 2020; Dors et al. 2020a). Therefore, the total metal- gas atoms can not combine in molecules formation and they licity (Z) of the gas phase from emission lines emitter objects can not be trapped in dust grains due to their quantum con- is commonly traced by the oxygen abundance relative to hy- figuration, unlike the oxygen which is depleted onto dust in order of 0.1 dex (e.g. Izotov et al. 2006; Pilyugin et al. 2007). ? E-mail: [email protected] The Te-method, which is based on direct estimation of

© 2020 The Authors Astronomy & Astrophysics manuscript no. aanda ©ESO 2021 September 23, 2021

Beyond the dips of V807 Tau, a spectropolarimetric study of a dipper’s magnetosphere? K. Pouilly1, J. Bouvier1, E. Alecian1, S.H.P. Alencar2, A.-M. Cody3, J.-F. Donati4, K. Grankin5, L. Rebull6, and C.P. Folsom7, 8

1 Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France e-mail: [email protected] 2 Departamento de Fisica - ICEx - UFMG, Av. Antônio Carlos 6627, 30270-901 Belo Horizonte, MG, Brazil 3 SETI Institute, 189 N Bernardo Ave #200, Mountain View, CA 94043 USA 4 Univ. de Toulouse, CNRS, IRAP, 14 avenue Belin, 31400 Toulouse, France 5 Crimean Astrophysical Observatory, Nauchny, Crimea 298409 6 Infrared Science Archive (IRSA), IPAC, California Institute of Technology, 1200 E. California Blvd, MS 100-22, Pasadena, CA 91125 USA 7 Department of Physics & Space Science, Royal Military College of Canada, PO Box 17000 Station Forces, Kingston, ON, Canada K7K 0C6 8 Tartu Observatory, University of Tartu, Observatooriumi 1, Tõravere, 61602 Tartumaa, Estonia

Received March 22, 2021; accepted September 13, 2021

ABSTRACT

Context. The so-called dippers are pre-main-sequence objects that accrete material from their circumstellar disks through the stellar magnetosphere. Their unique type of variability allows us to probe the magnetic star-disk interaction processes in young stellar objects. Aims. We aim to characterize the magnetospheric accretion process in the young stellar object V807 Tau, one of the most stable dippers revealed by K2 in the Taurus star forming region. Methods. We performed photometric and spectropolarimetric follow-up observations of this system with CFHT/ESPaDOnS in order to investigate the variability of the system over several rotational periods. Results. We derive a 4.38 day period from the K2 dipper light curve. This period is also seen in the radial velocity variations, which we ascribe to spot modulation. The slightly redshifted narrow component of the He I 5876 Å line as well as the high velocity red wing of the Hβ and Hγ emission line profiles also vary in intensity with the same periodicity. The former traces the accretion shock at the stellar surface, and the latter is a signature of an accretion funnel flow crossing the line of sight. We derive a surface brightness map and the topology of the surface magnetic field from the modeling of Stokes I and V profiles, respectively, for photospheric lines and for the He I emission line. The latter reveals a bright spot at the stellar surface, located at a latitude of 60◦, and a maximum field strength of ∼2 kG at this location. The topology of the magnetic field at the stellar surface is dominated by a dipolar component inclined by about 40◦ onto the spin axis. Variable blueshifted absorption components seen in the Balmer line profiles suggest episodic outflows. Despite of its clear and stable dipper behavior, we derive a relatively low inclination of 40◦ to 50◦ for this system, which calls question the origin of the dips. The low inclination we infer is also consistent with the absence of deep inverse P Cygni components in the line profiles. Conclusions. We conclude that magnetospheric accretion is ongoing in V807 Tau, taking place through non-axisymmetric accretion funnel flows controlled by a strong, tilted, and mainly dipolar magnetic topology. Whether an inner disk warp resulting from this process can account for the dipper character of this source remains to be seen, given the low inclination of the system. Key words. Stars: variables: T Tauri - Stars: pre-main sequence - Accretion, accretion disk - Stars: magnetic field - Stars: individual: V807 Tau - Stars: starspots - Stars: winds

1. Introduction the inner disk, leading to accretion funnel flows controlled by the magnetic field lines (Bouvier et al. 2007b; Hartmann et al. arXiv:2109.10589v1 [astro-ph.SR] 22 Sep 2021 The study of classical T Tauri stars (CTTSs) provides impor- 2016). The magnetospheric star-disk interaction takes place on a tant constraints on the evolution of young stellar objects. These scale of ≤0.1 au, which cannot be resolved by direct imaging. low-mass stars actively accrete material from their circumstel- lar disks and are characterized by strong photometric and spec- tral variability, part of which arises from the star-disk interaction A common approach for investigating the star-disk interac- process. Their strong surface magnetic field is able to truncate tion process is to monitor the system over several rotational pe- riods, using photometry, spectroscopy, spectropolarimetry, and, ? Based on observations obtained at the Canada-France-Hawaii Tele- more recently, interferometry (e.g., Bouvier et al. 2007a; Do- scope (CFHT), which is operated by the National Research Council of nati et al. 2010, 2019; Alencar et al. 2012, 2018; Pouilly et al. Canada, the Institut National des Sciences de l’Univers of the Centre 2020; Bouvier et al. 2020a,b). Such studies have revealed strong National de la Recherche Scientifique of France, and the University of magnetic fields, on the order of 1 kG at the stellar surface, of- Hawaii. ten dominated by a dipolar component, as well as cool magnetic

Article number, page 1 of 18 Astronomy & Astrophysics manuscript no. AA_arxiv ©ESO 2021 September 23, 2021

Discs and outflows in the early phases of massive star formation: influence of magnetic fields and ambipolar diffusion B. Commerçon1 , M. González2, R. Mignon-Risse3, P. Hennebelle2 and N. Vaytet4

1 Univ Lyon, Ens de Lyon, Univ Lyon 1, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, 69007, Lyon, France e-mail: [email protected] 2 AIM, CEA, CNRS, Université Paris-Saclay, Université de Paris, 91191 Gif-sur-Yvette, France 3 Université de Paris, CNRS, AstroParticule et Cosmologie, F-75013 Paris, France 4 Data Management and Software Centre, European Spallation Source ERIC, Ole Maaløes Vej 3, 2200, Copenhagen, Denmark ...; ...

ABSTRACT

Context. Massive star formation remains one of the most challenging problems in astrophysics, as illustrated by the fundamental issues of the radiative pressure barrier and the initial fragmentation. The wide variety of physical processes involved, in particular the protostellar radiative feedback, increases the complexity of massive star formation in comparison with its low-mass counterpart. Aims. We aim to study the details of mass accretion and ejection in the vicinity of massive star forming cores using high-resolution (5 au) 3-dimensional numerical simulations. We investigate the mechanisms at the origin of outflows (radiative force versus mag- netic acceleration). We characterise the properties of the disc forming around massive protostars depending on the physics included: hydrodynamics, magnetic fields, and ambipolar diffusion. Methods. We use state-of-the-art 3-dimensional adaptive-mesh-refinement models of massive dense cores collapse which integrate the equations of (resistive) grey radiation-magnetohydrodynamics, and include sink particles evolution. For the first time, we include both protostellar radiative feedback via pre-main sequence evolutionary tracks and magnetic ambipolar diffusion. To determine the role of magnetic fields and AD play in the formation of outflows and discs, we studied 3 different cases: a purely hydrodynamical run, a magnetised simulation under the ideal approximation (perfect coupling), and a calculation with ambipolar diffusion (resistive case). In the most micro-physically complex model (resistive MHD), we also investigated the effect the initial amplitude of both magnetic field and solid body rotation have on the final properties of the massive protostellar system. We use simple criteria to identify the outflow and disc material and follow their evolution as the central star accretes mass up to 20 M in most of our models. The radiative, magnetic, and hydrodynamical properties of the outflows and discs are quantitatively measured and cross compared-between models. Results. Massive stars form in all our models, together with outflows and discs. The outflow is completely different when magnetic fields are introduced, so that magneto-centrifugal processes are the main driver of the outflow up to stellar masses of 20 M . Then, the disc properties heavily depend on the physics included. In particular, the disc formed in the ideal and resistive runs show opposite properties in terms of plasma beta, i.e. the ratio of thermal to magnetic pressures, and of magnetic fields topology. While the disc in the ideal case is dominated by the magnetic pressure and the toroidal magnetic fields, the one formed in the resistive runs is dominated by the thermal pressure and has essentially a vertical magnetic fields in the inner regions (R < 100 − 200 au). Conclusions. We find that magnetic processes dominate the early evolution of massive protostellar systems (M? < 20 M ) and shapes the accretion/ejection as well as the disc formation. Ambipolar diffusion is mainly at work at disc scales and regulates its properties. We predict magnetic fields topology within the disc and outflows, as well as disc masses and radii to be compared with observations. Last, our finding for the outflow and disc properties are reminiscent of the low-mass star formation framework, suggesting that accretion and ejection in young massive and low-mass protostars are regulated by the same physical processes at the early stages. Key words. hydrodynamics – Magnetohydrodynamics (MHD) – Radiative transfer - Methods: numerical – Stars: formation

1. Introduction rier and the initial fragmentation. Thanks to the formidable in- crease of observational capabilities, it is now possible to probe Massive stars dominate from birth to death the energy in- massive star forming regions with an unprecedented angular res- put in the interstellar medium. As such, their formation mecha- arXiv:2109.10580v1 [astro-ph.SR] 22 Sep 2021 olution corresponding to sub-50 au scales (e.g. Beuther et al. nism needs to be better understood in order to put constraints on 2019; Maud et al. 2019). In parallel, heavy numerical develop- their impact on star formation and galaxy evolution as a whole. It ments are undertaken in order to describe accurately the physics is currently admitted that massive stars form in giant molecular of star formation in numerical models (see for a recent review clouds, which exhibit turbulent motions and magnetic fields (e.g. Teyssier & Commerçon 2019). Tan et al. 2014; Motte et al. 2018). The mechanisms by which massive stars get their mass in this giant, turbulent, and magne- A number of observational studies focus on the fragmenta- tised complex remain poorly constrained though. In particular, tion of massive cores which are thought to be the progenitors two major issues make it difficult to assemble a large reservoir of massive stars (e.g. Bontemps et al. 2010; Palau et al. 2013; of mass within a single massive star: the radiative pressure bar- Zhang et al. 2015; Fontani et al. 2016; Csengeri et al. 2017; Figueira et al. 2018; Sanhueza et al. 2019; Sanna et al. 2019a). Send offprint requests to: B. Commerçon There is no clear consensus on which physical processes is regu-

Article number, page 1 of 21 MNRAS 000,1–11 (2015) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Periodic activity from fast radio burst FRB180916 explained in the frame of the orbiting asteroid model

Guillaume Voisin1,2★, Fabrice Mottez1, Philippe Zarka3, 1 LUTh, Observatoire de Paris, PSL Research University, CNRS, Université de Paris, Sorbonne Université, 5 place Jules Janssen, 92190 Meudon, France 2 Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, The University of Manchester, Manchester M19 9PL, UK 3 LESIA & USN, Observatoire de Paris, CNRS, PSL Research University, SU/UP/UO, 92195 Meudon, France

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Observation of fast radio bursts (FRBs) are rising very quickly with the advent of specialised instruments and surveys, and it has recently been shown that some of them repeat quasi-periodically. In particular, evidence of a 푃 = 16.35 day period has been reported for FRB 180916.J0158+65. We seek an explanation within the frame of our orbiting asteroid model, whereby FRBs are produced in the plasma wake of asteroids immersed in the wind of a pulsar or a magnetar. We used the data reported by the CHIME/FRB collaboration in order to infer the orbital characteristics of asteroid swarms, and performed parametric studies to explore the possible characteristics of the pulsar, its wind, and of the asteroids, under the constraint that the latter remain dynamically and thermally stable. We found a plausible configuration in which a young pulsar is orbited by a main −3 ∼ 10 푀 companion with a period 3푃 = 49d, three times longer than the apparent periodicity 푃. Asteroids responsible for FRBs are located in three dynamical swarms near the L3, L4 and L5 Lagrange points, in a 2:3 orbital resonance akin to the Hildas class of asteroids in the Solar system. In addition, asteroids could be present in the Trojan swarms at the L4 and L5 Lagrange points. Together these swarms form a carousel that explains the apparent 푃 periodicity and dispersion. We estimated that the presence of at least a few thousand asteroids, of size ∼ 20km, is necessary to produce the observed burst rate. We show how radius-to-frequency mapping in the wind and small perturbations by turbulence can suffice to explain downward-drifting sub-pulses, micro-structures, and narrow spectral occupancy. Key words: (Stars:) pulsar – Minor planets, asteroids: general – Relativistic processes – Radio continuum: general

1 INTRODUCTION perfect periodicity has not been observed so far, we believe that this theory has a few characteristics that are appealing for any theory of Fast radio bursts (FRBs) consist in short, typically a few millisec- FRBs. Alfvén wings are akin to a plasma wake left by an electrically onds, and intense flashes that have so far been observed only in conducting object immersed in a magnetised wind. They have been radio bands (Petroff et al. 2019). One of their most puzzling proper- originally theorised and observed in the context of the Jupiter-Io ties is the large electron column density that the signal has crossed, interaction (Neubauer 1980). These wings are favourable sites for the so-called dispersion measure (DM), which is encoded in their radio-emitting plasma instabilities such as the cyclotron maser insta- dynamic spectrum (e.g. Lyne & Graham-Smith(2012)). The mea- bility (Mottez & Zarka 2014). In the particular case of a pulsar the sured DM is compatible with extragalactic and even cosmological wind is ultra-relativistic with the immediate consequence that any distances. If the observed DM could in principle be caused by a radiation in the wing is highly collimated within a cone of aperture dense environment, the extra-galactic origin of 12 FRBs has been arXiv:2109.10577v1 [astro-ph.HE] 22 Sep 2021 ∼ 1/훾 where 훾 . 106 is the Lorentz factor of the wind. Thus, a confirmed by the identification of their host galaxies (see Bochenek few interesting characteristics follow: i) compared to isotropic emis- et al. 2021, and references therein), and in particular the host of FRB sion scenarios that require tremendous amounts of energy, here little 180916.J0158+65 (Marcote et al. 2020). In addition, the repetition is necessary to produce the observed radio flux, ii) no high-energy of FRB 180916.J0158+65 and at least 19 others 1 rules out, in those counterpart is required which is in agreement with most of the ob- cases, theories appealing to cataclysmic events such as mergers or servations, iii) we know that objects orbiting pulsars are common, collisions. Fast radio bursts caused by Alvén wings of planets orbit- only is it unlikely to observe a galactic FRB due to their very a nar- ing pulsars (Mottez & Zarka 2014) have been, to our knowledge, the row beam, and iv) this very narrow beam explains the short burst first theory predicting a periodic repetition of FRBs. Although such duration. Concerning ii) and iii), we note that the FRB-like event recently observed from the Galactic magnetar SGR 1935+2154 (An- dersen et al. 2020) was coincident with soft gamma-ray emissions. ★ E-mail: [email protected]; However, the lack of radio bursts during other high-energy events [email protected] (Lin et al. 2020) suggests that the two components could result from 1 Twenty repeaters were identified by CHIME/FRB by June 2021: https: two different mechanisms. //www.chime-frb.ca/repeaters

© 2015 The Authors arXiv:2109.10570v1 [astro-ph.GA] 22 Sep 2021 srnm etr,22,Vl 7 o7 p 454–473. pp. 7, No 47, Vol. 2021, Letters, Astronomy ntr,i eemndb h vlto fteGlci oeta.Ti e This o formation potential. the Galactic as the such ones, of secular evolution including the processes, by many on determined is turn, in netgto sapwru en fsuyn h hsclcondition physical the evolution. studying its of of formatio comparable course galaxy means is the earliest powerful and the a Gyr of is 13 witnesses investigation reaches are age they Basically, whose dista verse. objects Galactocentric clus oldest at globular distributed the 150 and are than center more Galactic numbers the Galaxy around Way Milky our present, At p evolving INTRODUCTION and static t both provide We in derived static. poten clusters deemed errors globular static be the of may by a c potential parameters caused of Galactic the uncertainties the case of statistical scales the the sizes time of to changing limits compared and the changes masses into changing minor average, both with on potential a of fdsatgoua lsesadtoewt ag pcne distan apocenter to up large scales a time on with that those show compensa We and maximally ch changes. param clusters is The orbital globular influence the distant components. on this of act Galactic distances to Galactocentric the shown small are of At components sizes se the of the of influence in sizes the and change studied th have a integrating we and when time masses first potentials the evolving orbital For the and backward. compare Gyr static We in time. derived in backward clusters Navarro–Frenk integrated spherical are a orbits and G The disk, the Miyamoto–Nagai Ha of flat in a potential bulge, described model algorithm Va three-component the changin axisymmetric and used with have adopted Baumgardt we potential from components, model Galactic distances pro evolving the with mean semicosmological clusters and a globular catalogue 152 construct of EDR3 motion Gaia orbital the the on Galaxy Way Milky 1 eyipratapc fivsiaigGsi td fterorbit their of study a is GCs investigating of aspect important very A Abstract -al [email protected] e-mail: td fteIflec fa vligGlci Potential Galactic Evolving an of Influence the of Study nteObtlPoete f12Goua clusters Globular 152 of Properties Orbital the on W aesuidteiflec fa vliggaiainlpotentia gravitational evolving an of influence the studied have —We ukv srnmclOsraoy usa cdm fScie of Academy Russian Observatory, Astronomical Pulkovo ihDt rmteGi D3Catalogue EDR3 Gaia the from Data with .T Bajkova T. A. ukvkes.6,S.Ptrbr,164 Russia 196140 Petersburg, St. 65, sh. Pulkovskoe 1 .A mro n .V Bobylev V. V. and Smirnov A. A. , − 1 y h riso lblrcutr ntecase the in clusters globular of orbits the Gyr 5 lx nldsaspherical a includes alaxy Wiedr atrhalo. matter dark –White pc n,hne their hence, and, epoch n tr nteopst way. opposite the in eters a PrzVlea et (Perez-Villegas bar a f csu o20kc GCs kpc. 200 to up nces these on So, data. the in aaeyacag nthe in change a parately il hs hne fit changes These tial. eudrotebiggest the undergo ce otentials. aaeeso globular of parameters nteUies over Universe the in s oteaeo h Uni- the of age the to assadszsof sizes and masses g h ta.(05.The (2015). al. et ghi be ihteorbital the with ables risfr5ad12 and 5 for orbits e moet undergo, omponents ne ntemasses the in anges nces, e o.Teorbits The for. ted es(C)orbiting (GCs) ters iiv(01.To (2021). siliev e oin from motions per lhsoy which, history, al ouindepends volution fthe of l arXiv:2109.10546v1 [astro-ph.GA] 22 Sep 2021 ★ © hsclcniin n hmclaudne nP -6 Re 2-36. M observations PN echelle in deep abundances chemical and conditions Physical ieec sqatfidb en fteAudneDiscrepanc Abundance the of T ( means higher. (ADF) by Factor being always quantified ORLs systematical is from yield difference those ORLs with and values, CELs different from estimated that abundances found been has it sets, both using abundances determine 1 epeetasetu ftepaeaynbl -6obtained 2-36 M nebula planetary the of spectrum a present We ABSTRACT ZZZ form original in YYY; Received XXX. Accepted MNRAS acaRjs&Etbn2007 Esteban & García-Rojas 2017 al. et Peimbert https://www.nebulousresearch.org/adfs/ ADF al h aeo nihetb jcafo eta iayst binary central from ejecta by enrichment of case the tably ( composition chemical different with phases inhomogengas thermal discrepancy: abundance ( the of origin the some found, been have (e.g. values 100 extreme to close and (2–10) moderate PNe for (see used be also can (ORLs) lines recombination optical o ion for tteVr ag eecp VT.46eiso ie r det are lines emission 446 (VLT). Telescope Large Very the at utpeeeto eprtr ( temperature electron multiple eie h lsia E bnacs eas rsn robus O present of also ORLs we abundances, CEL classical the Besides O ratios, (CEL) line excited collisionally .0 rmaltedffrn ierto htw td,w n t si find to we unable study, are we models that photoionization two-phased ratios that line find different the all From 0.50. CL)o ev lmns ihasizeable a with elements, heavy of (CELs) H c of study The INTRODUCTION 1 words: chemi Key and physical its since resolution, spectral diagno high line suc with density reproduce and adequately temperature to observed the able reproduce model three-phased a find we 1 oéN Espíritu N. José emet1967 Peimbert nttt eAtooí,UiesddNcoa uóoade Autónoma Nacional Universidad Astronomía, de Instituto oe esnhscmie ito bet ihmaue ADF measured with objects of list a compiled has Wesson Roger ii 01TeAuthors The 2021 -al [email protected] E-mail: ubro hsclseaishv enpooe oexplain to proposed been have scenarios physical of number A H ii ein a enmsl ae ncliinlyectdline excited collisionally on based mostly been has regions ( - ein yial xii oeaeAFvle,aon (e. 2 around values, ADF moderate exhibit typically regions - + 000 8 + ) 8 . = , ++ 1 – = = S:audne lntr eua:idvda:M2-36. M individual: nebulae: planetary – abundances ISM: eia bnacsi lntr eua Pe and (PNe) nebulae planetary in abundances hemical eoti e au o h bnac iceac atr( Factor Discrepancy Abundance the for value new a obtain we , i ta.2001 al. et Liu 20 ( ; ( - - emet&Cseo1969 Costero & Peimbert cabe l 2013 al. et McNabb (2021) + + 8 8 n eeecsteen.Weei spsil to possible is it Where therein). references and , ) ) ORL CEL 1 ★ , noi Peimbert Antonio & ): ; eaGerr ta.2012 al. et Peña-Guerrero ) 4 n est ( density and ) ). 1 ;ceitneo multiple of coexistence ); ubro bet where objects of number ++ pia eobnto ie OL) n esrn h inte the measuring and (ORLs), lines recombination optical i ta.2000 al. et Liu 1 = 4 éio po otl724 iddd éio050 México 04510, México de Ciudad 70264, Postal Apdo. México, diagnostics. ) ;wiein while ); r nPNe in ars chemical example rpit2 etme 01Cmie sn NA L MNRAS using Compiled 2021 September 23 Preprint eities ,no- ), a tutr a emr opiae hnpeiul thoug previously than complicated more be may structure cal his (1) ly g. tc.Oraayi hw hti sipratt td ihA high study to important is it that shows analysis Our stics. s: utnosyrpoueciia O critical reproduce multaneously y s oi bnacsfo Rso ev lmns rmCL and CELs From elements. heavy of ORLs from abundances ionic t a h betcno eceial ooeeu;moreover, homogeneous; chemically be cannot object the hat ce.W efr naayi fteceia opsto us composition chemical the of analysis an perform We ected. sn h lr iltadVsa cel pcrgah(UVE Spectrograph Echelle Visual and Violet Ultra the using ais hl ecnie hst eatymdl ti beto able is it model, toy a be to this consider we While ratios. h ) 4 e opoieisgtit hspolm h eeto fCE of detection the problem, this t most into clear of insight is contribution provide it a to Nonetheless, to der due scenarios. physical be may proposed cause the real The object. ery ( acsaebsdo E oOLln aiswihso strong a show which ratios line ORL temperature. to on CEL dependence on based are dances 2018 ( distributions energy agrAF itt h rsneo hmclinhomogeneitie chemical of presence the to hint ADFs larger an ( remains astrophysics discrepancy contemporary abundance in the problem of emissi origin of the behavior hence observed and the reproduced fully has narios ( stars binary period short from fields radiation ing ettv hmclcmoiino htinzdrgo,i region, defined. o be photoionized average a the of find composition to chemical important sentative is s it of reason, this calibration For and methods. formation galactic as such trophysics, oiino ohH both of position Dshsdmne ehnigo hsapproach. this of rethinking a demanded has ADFs rsneo hra nooeete ( inhomogeneities thermal of presence nw H known eete,snerto frcmiainlnsaenal t nearly ( are lines independent in recombination thermal of small ratios of since presence geneities, the in composition chemical the ord ta.2015 al. et Corradi and ti eylkl htteAFmyhv ieetoii nev- in origin different a have may ADF the that likely very is It eua bnacspa nipratrl nohrfilso a of fields other in role important an play abundances Nebular rdtoal,CL aebe sdt siaetechemical the estimate to used been have CELs Traditionally, ntecs fPe ml D aus( values ADF small PNe, of case the In o hmclyhmgnospoooie ein eg alm (e.g. regions photoionized homogeneous chemically For .Dsieteeipratavne,nn ftepooe sc proposed the of none advances, important these Despite ). = 4 ii r optduigavreyo ehd,including methods, of variety a using computed are ein) Rscudpoieabte ersnainof representation better a provide could ORLs regions), D)o hsnbl,bigADF(O being nebula, this of ADF) emet&Pibr 2013 Peimbert & Peimbert ii ; ein n N,hwvrtedsoeyo large of discovery the however PNe, and regions esne l 2018 al. et Wesson ihlse l 2012 al. et Nicholls ii n [O and emete l 2014 al. et Peimbert acaRjse l 2019 al. et García-Rojas st fteBle jump. Balmer the of nsity ;nnMxela electron non-Maxwellian ); ≤ , iii ;mawie E abun- CEL meanwhile, ); 2013 )cnb xlie ythe by explained be can 5) ierto.However, ratios. line ] ut from sults ;advrigioniz- varying and ); atsa&Ahmed & Bautista A T E ++ tl l v3.0 file style X ) emperature = rn line trong a nor- in hat hscan this f ht. FPNe DF ;while ); nlines on 6.76 repre- r s. sand Ls homo- s all ost com- open ). ing we S) of e- s- ± MNRAS 000,1–14 (2020) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Identifying Potential Exomoon Signals with Convolutional Neural Networks

Alex Teachey1,2★ & David Kipping2,3 1Academia Sinica Institute of Astronomy & Astrophysics, Taipei, Taiwan R.O.C. 2Department of Astronomy, Columbia University in the City of New York, USA 3Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010 USA

Accepted 15 September 2021. Received 12 August 2021; in original form 23 October 2020

ABSTRACT Targeted observations of possible exomoon host systems will remain difficult to obtain and time-consuming to analyze in the foreseeable future. As such, time-domain surveys such as Kepler, K2 and TESS will continue to play a critical role as the first step in identifying candidate exomoon systems, which may then be followed-up with premier ground- or space- based telescopes. In this work, we train an ensemble of convolutional neural networks (CNNs) to identify candidate exomoon signals in single-transit events observed by Kepler. Our training set consists of ∼27,000 examples of synthetic, planet-only and planet+moon single transits, injected into Kepler light curves. We achieve up to 88% classification accuracy with individual CNN architectures and 97% precision in identifying the moons in the validation set when the CNN ensemble is in total agreement. We then apply the CNN ensemble to light curves from 1880 Kepler Objects of Interest with periods > 10 days (∼57,000 individual transits), and further test the accuracy of the CNN classifier by injecting planet transits into each light curve, thus quantifying the extent to which residual stellar activity may result in false positive classifications. We find a small fraction of these transits contain moon-like signals, though we caution against strong inferences of the exomoon occurrence rate from this result. We conclude by discussing some ongoing challenges to utilizing neural networks for the exomoon search. Key words: planets and satellites: detection

1 INTRODUCTION Recent work (e.g Shallue & Vanderburg 2018; Pearson et al. 2018; Dattilo et al. 2019; Yu et al. 2019; Osborn et al. 2020) with The identification of candidate exomoon signals in time-domain the Kepler, K2, and TESS dataset have shown that convolutional photometry remains a significant challenge. With such an enormous neural networks (CNNs) can be effective in distinguishing genuine data volume to comb, and rigorous Bayesian model selection so exoplanet transits from false positives. These efforts employed two computationally expensive (e.g. Kipping et al. 2015), how does one light curve inputs – a “local” and a “global” view – and in some cases efficiently identify potentially interesting signals? Increasingly, the also incorporated additional observables as inputs, achieving up to astronomical community has turned to the use of machine learning 96% accuracy and confirming a number of new planet candidates. for approaching such problems (see a comprehensive review from Fluke & Jacobs 2020). Indeed, mentions of the term in refereed In most of the works highlighted above (the exception being Osborn et al. 2020), the networks were trained on a sample of pre-

arXiv:2109.10503v1 [astro-ph.EP] 22 Sep 2021 publication abstracts has grown exponentially in the last decade, from about 200 articles in 2010 to more than 4000 in 2020 – a growth viously confirmed systems. This is certainly a reasonable approach, rate far exceeding that of the yearly literature increase. In some since it incorporates the great diversity of systems so far discovered applications, machine learning may be able to identify patterns that and avoids the requirement of producing simulated systems. At the human beings cannot (in, for example, higher dimensional spaces same time, there are some downsides to this methodology, namely, where a system state is determined by more than three variables). In that it 1) runs the risk of inadvertently teaching the networks about other contexts, machine learning algorithms may be only as good (and therefore encoding) the inherent biases of previous detection as what the human eye can discern (for example, computer vision), / confirmation methods, 2) limits the training sample to those sys- but are nevertheless employed for their speed (e.g. autonomous tems that have already been detected, and 3) ground-truth values vehicles) and / or capacity for handling vast amounts of data (e.g. are not known with perfect precision or accuracy, adding noise to mail sorting). the training set. We note, however, that there is some evidence that CNNs can be robust against even high levels of label noise (Rolnick et al. 2017), making this last point not especially worrisome. ★ E-mail: [email protected] In this work, we apply an ensemble of CNNs to the problem

© 2020 The Authors MNRAS 000,1–24 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

The density distribution and the physical origins of density intermittency in sub- to trans-Alfv´enic supersonic turbulence

James R. Beattie 1?, Philip Mocz 2>,3 , Christoph Federrath 1,4 and Ralf S. Klessen 5,6 1Research School of Astronomy and Astrophysics, Australian† National University, Canberra,‡ ACT 2611, Australia 2Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA 3Lawrence Livermore National Laboratory, 4 Ivy Lane, 7000 East Ave, Livermore, CA 94550, USA 4Australian Research Council Centre of Excellence in All Sky Astrophysics (ASTRO3D), Canberra, ACT 2611, Australia 5Universit¨at Heidelberg, Zentrum fur¨ Astronomie, Institut fur¨ Theoretische Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany 6Universit¨at Heidelberg, Interdisziplin¨ares Zentrum fur¨ Wissenschaftliches Rechnen, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany >Einstein Fellow

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT

The probability density function (PDF) of the logarithmic density contrast, s = ln(ρ/ρ0), with gas density ρ and mean density ρ0, for hydrodynamical supersonic turbulence is well-known to have significant signatures of intermittency that monotonically increase with the turbulent Mach number, . By studying the mass- and volume-weighted s-PDF M for an ensemble of 16 sub- to trans-Alf´venic mean-field, supersonic, isothermal turbulence simulations, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the intermittency of s. Using four independent measures of the intermittency we find hydrodynamical-like intermittency in the highly magnetised plasma for 4. However, for 4, the signatures of intermittency disappear, leaving approximately M . M & lognormal s-statistics – exactly the opposite of hydrodynamical turbulence in the high- limit. To understand the M 4 intermittency we use one-dimensional (1D) pencil beams to explore the dynamics along and across the mean M . magnetic field, B 0. We discuss kinetic, density and magnetic field fluctuations from the pencil beams, and identify physical sources of intermittency as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 and create large, volume-poor under-densities. These under-densities contribute significantly to the skewness of the s-PDF. We confirm this result independently using 1D fluid shock simulations. We discuss the Gaussianisation of the 4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenisation of M & the dynamical timescales between the over- and under-dense regions in the compressible gas. Key words: MHD – turbulence – ISM: kinematics and dynamics – ISM: magnetic fields – ISM: structure

1 INTRODUCTION an approximately normal volume-weighted s-PDF1, where s ln(ρ/ρ0), and ρ is the cloud density, with ρ0 the value of The PDF of the gas density is a valuable tool for understand- ≡ the volume-weighted mean. It follows ing the nature of compressible turbulence and star forma- ( 2 ) tion in the interstellar medium of galaxies (Vazquez-Semadeni 2 1 (s s0,V ) pN,V (s; σs,V ) = exp − 2 , (1) 1994; Padoan et al. 1997; Passot & V´azquez-Semadeni 1998; q 2 arXiv:2109.10470v1 [astro-ph.GA] 22 Sep 2021 2πσ − 2σs,V Federrath et al. 2008b; Federrath et al. 2009; Brunt et al. s,V 2 2010a; Brunt et al. 2010b; Krumholz & McKee 2005; Hen- σs,V s0,V = , (2) nebelle et al. 2011; Federrath & Klessen 2012; Burkhart & − 2 Lazarian 2012; Konstandin et al. 2012b; Molina et al. 2012; 2 σs,V = f( , A, b, γ, Γ). (3) Hopkins 2013; Nolan et al. 2015; Federrath & Banerjee 2015; M M 2 Squire & Hopkins 2017; Pan et al. 2019; Mocz & Burkhart The log-density variance, σs , captures the density variations 2019; Menon et al. 2021, 2020; Khullar et al. 2021). Simple induced by different physical processes in a MC. It is a func- models of molecular clouds (MCs) in the interstellar medium tion of (i) the turbulent Mach number,

(ISM), which are supersonic and magnetised, and have not = σV /cs, (4) yet started to collapse under their own self-gravity have M where σV is the velocity dispersion on system scale L, and cs is the sound speed (Vazquez-Semadeni 1994; Padoan et al. ? E-mail: [email protected] E-mail: [email protected] † E-mail: [email protected] 1 Or equivalently, a lognormal ρ-PDF. ‡

© 2021 The Authors Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

ALMA Imaging of a Galactic Molecular Outflow in NGC 4945

Alberto D. Bolatto,1, 2, 3 Adam K. Leroy,4 Rebecca C. Levy,1 David S. Meier,5, 6 Elisabeth A. C. Mills,7 Todd A. Thompson,4 Kimberly L. Emig,8, ∗ Sylvain Veilleux,1, 3 Jürgen Ott,6 Mark Gorski,9 Fabian Walter,10, 6 Laura A. Lopez,4 and Laura Lenkić1

1Department of Astronomy, University of Maryland, College Park, MD 20742, USA 2Visiting Scholar at the Flatiron Institute, Center for Computational Astrophysics, NY 10010, USA 3Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA 4Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA 5Department of Physics, New Mexico Institute of Mining and Technology, 801 Leroy Pl., Socorro, NM, 87801, USA 6National Radio Astronomy Observatory, P. O. Box O, 1003 Lopezville Rd., Socorro, NM, 87801, USA 7Department of Physics and Astronomy, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, KS 66045, USA 8National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, USA 9Chalmers University of Technology, Department of Space, Earth, and Environment, SE-412 96 Gothenburg, Sweden 10Max-Planck-Institut für Astronomie, Königstuhl 17, D-69120 Heidelberg, Germany

ABSTRACT We present the ALMA detection of molecular outflowing gas in the central regions of NGC 4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO J = 3 − 2 at ∼ 0.300 resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the disk. The fastest and brightest of these plumes has emission reaching observed line-of-sight projected velocities of over 450 km s−1 beyond systemic, equivalent to an estimated physical outflow velocity −1 v & 600 km s for the fastest emission. Most of these plumes have corresponding emission in HCN or HCO+ J = 4 − 3. We discuss a kinematic model for the outflow emission where the molecular gas has the geometry of the ionized gas cone and shares the rotation velocity of the galaxy when ejected. We use this model to explain the velocities we observe, constrain the physical speed of the ejected material, and account for the fraction of outflowing gas that is not detected due to confusion with the −1 galaxy disk. We estimate a total molecular mass outflow rate M˙ mol ∼ 20 M yr flowing through a surface within 100 pc of the disk midplane, likely driven by a combination of the central starburst and AGN.

Keywords: galaxies: individual (NGC 4945) — galaxies: starburst — ISM: kinematics — ISM: molecules

1. INTRODUCTION and NGC 253, NGC 4945 is one of the brightest far- NGC 4945, located in the southern hemisphere, is one infrared sources outside the Local Group (Sanders et al. of the nearest massive galaxies (3.8 ± 0.3 Mpc corre- 2003). And like these two examples, it is one of the very arXiv:2109.10437v1 [astro-ph.GA] 21 Sep 2021 sponding to a scale of 18.4 pc per arcsecond; Karachent- few starbursts seen in gamma-ray emission by Fermi, sev et al. 2007). The galaxy harbors a central star- although it is unclear if the emission is due to the star burst (Bendo et al. 2016; Emig et al. 2020) as well as a formation activity (and the ensuing supernovae) or to Seyfert 2 Active Galactic Nucleus (AGN, Schurch et al. the central AGN (Ackermann et al. 2012). 2002). Similar to the nearby starburst galaxies M 82 NGC 4945 is viewed almost edge-on, with an axis ratio of nearly 5, corresponding to an inclination i > 80◦. In- spection of the extinction apparent in optical and near- Corresponding author: Alberto D. Bolatto infrared images shows that the southwest rim of the [email protected] galaxy lies closest to us. Thus the galaxy disk obscures ∗ Jansky Fellow of the National Radio Astronomy Observatory much of what happens south of the nucleus while direc- Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

Precise Masses and Orbits for Nine Radial Velocity Exoplanets

Yiting Li ,1 Timothy D. Brandt ,1 G. Mirek Brandt ,1, ∗ Trent J. Dupuy ,2 Daniel Michalik ,3 Rebecca Jensen-Clem ,4 Yunlin Zeng ,5 Jacqueline Faherty ,6 and Elena L. Mitra 7

1Department of Physics, University of California, Santa Barbara, Santa Barbara, CA 93106, USA 2Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 3European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands 4Astronomy & Astrophysics Department, University of California, Santa Cruz, CA 95064, USA 5School of Physics, Georgia Institute of Technology, 837 State Street, Atlanta, Georgia, USA 6American Museum of Natural History, New York, NY 7CUNY Hunter College, New York, NY 10065, USA

ABSTRACT Radial velocity (RV) surveys have discovered hundreds of exoplanetary systems but suffer from a fundamental degeneracy between planet mass Mp and orbital inclination i. In this paper we break this degeneracy by combining RVs with complementary absolute astrometry taken from the Gaia EDR3 version of the cross-calibrated Hipparcos-Gaia Catalog of Accelerations (HGCA). We use the Markov Chain Monte Carlo orbit code orvara to simultaneously fit literature RVs and absolute astrometry from the HGCA. We constrain the orbits, masses, and inclinations of nine single and massive RV companions orbiting nearby G and K stars. We confirm the planetary nature of six companions: +0.67 +1.0 +0.31 HD 29021 b (4.47−0.65MJup), HD 81041 b (7.24−0.37MJup), HD 87883 b (6.31−0.32MJup), HD 98649 b +2.3 +0.78 +3.6 (9.7−1.9MJup), HD 106252 b (10.00−0.73MJup), and HD 171238 b (8.8−1.3MJup). We place one com- +2.5 panion, HD 196067 b (12.5−1.8MJup) on the planet-brown dwarf boundary, and two companions in the +1.5 +2.0 low mass brown dwarf regime: HD 106515 Ab (18.9−1.4MJup), and HD 221420 b (20.6−1.6MJup). The +0.74 +3.0 brown dwarf HD 221420 b, with a semi-major axis of 9.99−0.70 AU, a period of 27.7−2.5 years, and +0.035 an eccentricity of 0.162−0.030 represents a promising target for high-contrast imaging. The RV orbits of HD 87883 b, HD 98649 b, HD 171238 b, and HD 196067 b are not fully constrained yet because of insufficient RV data. We find two possible inclinations for each of these orbits due to difficulty in separating prograde from retrograde orbits, but we expect this will change decisively with future Gaia data releases.

Keywords: —

1. INTRODUCTION The RV method measures the Doppler shift in the The radial velocity (RV) technique is one of the ear- spectrum of a star perturbed by one or multiple un- liest and most commonly used indirect methods in the seen orbiting companions (Walker 2012; Plavchan et al. detection of exoplanets in both nearby and distant so- 2015; Fischer et al. 2016; Wright 2018). The RV semi- lar systems (Hatzes 2016; Wei 2018; Wright 2018). 51 amplitude due to a companion scales with the compan- arXiv:2109.10422v1 [astro-ph.EP] 21 Sep 2021 Pegasi b, the first exoplanet around a Sun-like star, was ion’s mass, the sine of the orbital inclination (sin i), and discovered via RV in 1995 (Mayor & Queloz 1995). Since the inverse of the square root of the semi-major axis. then, RV instruments have overcome many design and This means that massive and close-in RV companions performance challenges to reach extreme m/s precisions, are usually more easily and precisely characterized ow- at once revealing the richness and diversity of the popu- ing to complete orbital coverage. Much like 51 Pegasi lation of exoplanets (e.g. Udry et al. 2007; Akeson et al. b, these are often short-period Jupiter analogs orbiting 2013; Rice 2014). solar-type stars (Wright et al. 2003). Thanks to the growing temporal baselines of RV sur- veys, detection sensitivity is now being extended beyond ∗ NSF Graduate Research Fellow the ice line (& 3 AU) (Ford 2014; Wittenmyer et al. Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

The Population of Interstellar Objects Detectable with the LSST and Accessible for In Situ Rendezvous with Various Mission Designs

Devin Hoover ,1 Darryl Z. Seligman ,2 and Matthew J. Payne 3

1Dept. of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637 2Dept. of the Geophysical Sciences, University of Chicago, Chicago, IL 60637 3Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 51, Cambridge, MA 02138, USA

ABSTRACT The recently discovered population of interstellar objects presents us with the opportunity to char- acterize material from extrasolar planetary and stellar systems up close. The forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) will provide an unprecedented increase in sen- sitivity to these objects compared to the capabilities of currently operational observational facilities. In this paper, we generate a synthetic population of interstellar objects drawn from their galactic kine- matics, and identify the distribution of impact parameters, eccentricities, hyperbolic velocities and sky locations of objects detectable with the LSST. This population is characterized by a clustering of tra- jectories in the direction of the solar apex and anti-apex, centered at orbital inclinations of ∼ 90◦. We identify the ecliptic or solar apex as the optimal sky locations to search for future interstellar objects as a function of survey limiting magnitude. Moreover, we identify the trajectories of detectable objects that will be reachable for in-situ rendezvous with a dedicated mission with the capabilities of the forthcoming Comet Interceptor or proposed BRIDGE concept. By scaling our fractional population statistics with the inferred spatial number density, we estimate that the LSST will detect of order ∼ 50 interstellar objects over the course of its ∼ 10 year observational campaign. Furthermore, we find that there should be of order ∼ 10 and ∼ 0.05 reachable targets for missions with propulsion capabilities comparable to BRIDGE and Comet Interceptor, respectively. These number estimates will be readily updateable when the number density and size frequency distribution of interstellar objects is better constrained.

Keywords: Interstellar Objects

1. INTRODUCTION on the spatial distribution, size-frequency distribution The recent discoveries of the first two interstellar ob- and levels of cometary activity of these objects. The jects (ISOs), 1I/2017 U1 (‘Oumuamua) and C/2019 future detections and characterizations of these objects 2I (Borisov), imply that a galactic-wide population of will yield insights into the physical mechanisms at work roughly 1026 similar bodies exist, with spatial number in extra-solar planetary systems and star formation re- −1 −3 gions throughout the galaxy. densities of order no ∼ 2 × 10 au (Trilling et al. 2017; Laughlin & Batygin 2017; Do et al. 2018; Moro- The first interstellar object, ‘Oumuamua, defied all ex- Mart´ın 2019, 2018; Levine et al. 2021). These number pectations for an interstellar comet. It exhibited a dis- arXiv:2109.10406v1 [astro-ph.EP] 21 Sep 2021 densities have been estimated by incorporating the de- tinct lack of dust and typical cometary volatiles (Meech tailed detection capabilities of current observational fa- et al. 2017; Jewitt et al. 2017; Trilling et al. 2018), an cilities. Forthcoming observatories, such as the Rubin extreme 6 : 6 : 1 aspect ratio (Knight et al. 2017; Drahus Observatory Legacy Survey of Space and Time (LSST), et al. 2018; Bolin et al. 2018; McNeill et al. 2018; Bel- will provide unprecedented completeness for populations ton et al. 2018; Mashchenko 2019), a reddened reflection of minor bodies in the Solar System (Jones et al. 2009; spectrum (Masiero 2017; Fitzsimmons et al. 2018; Ban- Ivezi´cet al. 2019). It has been estimated that the LSST nister et al. 2017; Ye et al. 2017), and a non-gravitational will detect between 1−10 interstellar objects every year acceleration (Micheli et al. 2018). Curiously, it also ex- (Engelhardt et al. 2017; Cook et al. 2016; Trilling et al. hibited a surprisingly low incoming velocity with respect 2017; Seligman & Laughlin 2018). It is important to to the Local Standard of Rest (LSR) in the direction note that the detection frequency depends sensitively of the solar apex (Mamajek 2017), implying a young ICRC 2021 THE ASTROPARTICLE PHYSICS CONFERENCE ONLINE ICRC 2021Berlin | Germany THE ASTROPARTICLE PHYSICS CONFERENCE th Berlin37 International| Germany Cosmic Ray Conference 12–23 July 2021

The Compton Spectrometer and Imager Project for MeV Astronomy

John A. Tomsick푎,∗ on behalf of the COSI Collaboration (a complete list of authors can be found at the end of the proceedings) 푎Space Sciences Laboratory, 7 Gauss Way, University of California, Berkeley, CA 94720-7450, USA E-mail: [email protected]

The Compton Spectrometer and Imager (COSI) is a 0.2–5 MeV Compton telescope capable of imaging, spectroscopy, and polarimetry of astrophysical sources. Such capabilities are made possible by COSI’s germanium cross-strip detectors, which provide high efficiency, high resolution spectroscopy and precise 3D positioning of photon interactions. Science goals for COSI include studies of 0.511 MeV emission from antimatter annihilation in the Galaxy, mapping radioactive elements from nucleosynthesis, determining emission mechanisms and source geometries with polarization, and detecting and localizing multimessenger sources. The instantaneous field of view (FOV) for the germanium detectors is >25% of the sky, and they are surrounded on the sides and bottom by active shields, providing background rejection as well as allowing for detection of gamma-ray bursts or other gamma-ray flares over >50% of the sky. We have completed a Phase A concept study to consider COSI as a Small Explorer (SMEX) satellite mission, and here we discuss the advances COSI-SMEX provides for astrophysics in the MeV bandpass. arXiv:2109.10403v1 [astro-ph.IM] 21 Sep 2021

37th International Cosmic Ray Conference (ICRC 2021) July 12th – 23rd, 2021 Online – Berlin, Germany

∗Presenter

© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). https://pos.sissa.it/ Quenching of star formation from a lack of inflowing gas to galaxies

Katherine E. Whitaker1,2, Christina C. Williams3, Lamiya Mowla4, Justin S. Spilker5,6, Sune Toft2,7, Desika Narayanan8,2, Alexandra Pope1, Georgios E. Magdis2,9,7, Pieter G. van Dokkum10, Mohammad Akhshik11, Rachel Bezanson12, Gabriel B. Brammer2,7, Joel Leja13,14,15, Allison Man16, Erica J. Nelson17, Johan Richard18, Camilla Pacifici19, Keren Sharon20, Francesco Valentino2,7

Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old1. Very low amounts of molecular gas appear responsible for this, at least in some cases2-7, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases8-11. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly2-7, or stacked emission with its problems of averaging8,9. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts12-14. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies.

1 Department of Astronomy, University of Massachusetts, Amherst, MA, 01003 USA 2 Cosmic Dawn Center (DAWN), Copenhagen, Denmark 3 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 4 Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St George St, Toronto, ON M5S 3H4, Canada 5 Department of Astronomy, University of Texas at Austin, 2515 Speedway, Stop C1400, Austin, TX 78712, USA 6 NASA Hubble Fellow 7 Niels Bohr Institute, University of Copenhagen, Jagtvej 128, DK2200, Copenhagen, Denmark 8 Department of Astronomy, University of Florida, 211 Bryant Space Sciences Center, Gainesville, FL 32611 USA 9 DTU-Space, Technical University of Denmark, Elektrovej 327, DK-2800 Kgs. Lyngby, Denmark 10 Astronomy Department, Yale University, 52 Hillhouse Ave, New Haven, CT 06511, USA 11 Department of Physics, University of Connecticut, Storrs, CT 06269, USA 12 Department of Physics and Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA 13 Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 14 Institute for Computational & Data Sciences, The Pennsylvania State University, University Park, PA, USA 15 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 16 Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada 17 Department of Astrophysical and Planetary Sciences, 391 UCB, University of Colorado, Boulder, CO 80309-0391, USA 18 Univ Lyon, Univ Lyon1, Ens de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F-69230, Saint-Genis-Laval, France 19 Space Telescope Science Institute, Baltimore, MD 21218, USA 20 Department of Astronomy, University of Michigan, 1085 S. University Ave, Ann Arbor, MI 48109, USA Normal, Dust-Obscured Galaxies in the Epoch of Reionization

Y. Fudamoto1,2,3, P. A. Oesch1,4, S. Schouws5, M. Stefanon5, R. Smit6, R. J. Bouwens5, R. A. A. Bowler7, R. Endsley8, V. Gonzalez9,10, H. Inami11, I. Labbe12, D. Stark8, M. Aravena13, L. Barrufet1, E. da Cunha14,15, P. Dayal16, A. Ferrara17, L. Graziani18,20,27, J. Hodge5, A. Hutter16, Y. Li21,22, I. De Looze23,24, T. Nanayakkara12, A. Pallottini17, D. Riechers25, R. Schneider18,19,26,27, G. Ucci16, P. van der Werf5, C. White8

1Department of Astronomy, University of Geneva, 51 Ch. Pegasi, 23Sterrenkundig Observatorium, Ghent University, Krijgslaan 281 1290 Versoix, Switzerland - S9, 9000 Gent, Belgium 2Research Institute for Science and Engineering, Waseda Univer- 24Dept. of Physics & Astronomy, University College Lon- sity, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan don,Gower Street, London WC1E 6BT, UK 3National Astronomical Observatory of Japan, 2-21-1, Osawa, Mi- 25Cornell University, 220 Space Sciences Building, Ithaca, NY taka, Tokyo, Japan 14853, USA 4Cosmic Dawn Center (DAWN), Niels Bohr Institute, University 26Sapienza School for Advanced Studies, Viale Regina Elena 291, of Copenhagen, Jagtvej 128, København N, DK-2200, Denmark 00161 Roma Italy 5Leiden Observatory, Leiden University, PO Box 9500, 2300 RA 27INFN, Sezione di Roma 1, P.le Aldo Moro 2, 00185 Roma, Leiden, The Netherlands 6Astrophysics Research Institute, Liverpool John Moores Univer- Over the past decades, rest-frame ultraviolet (UV) ob- sity, 146 Brownlow Hill, Liverpool L3 5RF, UK servations have provided large samples of UV luminous 7Sub-department of Astrophysics, The Denys Wilkinson Building, galaxies at redshift (z) greater than 61–3, during the so- University of Oxford, Keble Road, Oxford, OX1 3RH, UK called epoch of reionization. While a few of these UV 8Steward Observatory, University of Arizona, 933 N Cherry Ave, identified galaxies revealed significant dust reservoirs4–7, Tucson, AZ 85721, USA very heavily dust-obscured sources at these early times 9Departmento de Astronomia, Universidad de Chile, Casilla 36-D, have remained elusive. They are limited to a rare pop- Santiago 7591245, Chile ulation of extreme starburst galaxies8–12, and compan- 13, 14 10Centro de Astrofisica y Tecnologias Afines (CATA), Camino del ions of rare quasars . These studies conclude that Observatorio 1515, Las Condes, Santiago 7591245, Chile the contribution of dust-obscured galaxies to the cosmic 11Hiroshima Astrophysical Science Center, Hiroshima Univer- star formation rate density at z > 6 is sub-dominant. sity, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Recent ALMA and Spitzer observations have identi- Japan fied a more abundant, less extreme population of ob- 15, 16 12Centre for Astrophysics & Supercomputing, Swinburne Univer- scured galaxies at z = 3 − 6 . However, this popu- sity of Technology, PO Box 218, Hawthorn, VIC 3112, Australia lation has not been confirmed in the reionization epoch 13Nucleo de Astronomia, Facultad de Ingenieria y Ciencias, Uni- so far. Here, we report the discovery of two dust- versidad Diego Portales, Av. Ejercito 441, Santiago, Chile obscured star forming galaxies at z = 6.6813 ± 0.0005 and z = 7.3521 ± 0.0005. These objects are not de- 14International Centre for Radio Astronomy Research, University of Western Australia, 35 Stirling Hwy, Crawley, WA 6009, Aus- tected in existing rest-frame UV data, and were only tralia discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV- 15ARC Centre of Excellence for All Sky Astrophysics in 3 Dimen- sions (ASTRO 3D) luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation 16Kapteyn Astronomical Institute, University of Groningen, PO rates than extreme starbursts, in line with typical star- Box 800, NL-9700 AV Groningen, the Netherlands 17 forming galaxies at z ∼ 7. This population of heavily Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, dust-obscured galaxies appears to contribute 10-25 per

arXiv:2109.10378v1 [astro-ph.GA] 21 Sep 2021 Italy cent to the z > 6 cosmic star formation rate density. 18Dipartimento di Fisica, Sapienza, Universita di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy As part of the ongoing ALMA large program REBELS 19INAF/Osservatorio Astronomico di Roma, via Frascati 33, 00078 (Reionization-Era Bright Emission Line Survey), we are Monte Porzio Catone, Roma, Italy observing 40 UV-luminous primary targets at z > 6.517. 20INAF/Osservatorio Astrofisico di Arcetri, Largo E. Femi 5, I- Among these targets are REBELS-12 and REBELS-29. 50125 Firenze, Italy When inspecting the ALMA data cube of these two sources, 21Department of Astronomy & Astrophysics, The Pennsylvania we identified two additional emission line neighbors. State University, 525 Davey Lab, University Park, PA 16802, USA The primary targets of our ALMA observing program, 22Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA REBELS-12 and REBELS-29, represent some of the most UV-luminous galaxies in this redshift range, and have

1 Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX631

AMBER: A Semi-Numerical Abundance Matching Box for the Epoch of Reionization

Hy Trac,1, 2 Nianyi Chen,1 Ian Holst,3 Marcelo A. Alvarez,4, 5 and Renyue Cen6

1McWilliams Center for Cosmology, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA 2NSF AI Planning Institute for Physics of the Future, Carnegie Mellon University, Pittsburgh, PA 15213, USA 3Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA 4Berkeley Center for Cosmological Physics, Department of Physics, University of California, Berkeley, CA 94720, USA 5Lawrence Berkeley National Laboratory, One Cyclotron Road, Berkeley, CA 94720, USA 6Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA

ABSTRACT The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism, but takes the novel approach of calculating the reionization-redshift field zre(x) assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fractionx ¯i(z). The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint zmid, duration ∆z, and asymmetry Az input parameters. The reionization process is further controlled through the minimum halo mass Mmin for galaxy formation and the radiation mean free path lmfp for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR.

1. INTRODUCTION Our current understanding of when the Epoch of Cosmic dawn is a fascinating period in the first billion Reionization (EoR) occurred primarily comes from “A years of the Universe and is a frontier topic for both the- Tale of Two Optical Depths”. The Thomson optical oretical and observational explorations and discoveries. depth τ (CMB photons scattering with free electrons) The reionization of hydrogen by the first stars, galaxies, and the Gunn-Peterson optical depth τGP (Lyman-alpha and quasars drastically converts the cold and neutral photons scattering with neutral hydrogen) have long gas into a warm and highly ionized medium. Galaxies provided two major observational constraints. However, most likely provided the bulk of the ionizing photons the latest measurements tell a starkly different tale than (e.g. Bouwens et al. 2015a; Finkelstein et al. 2015), but early observations (e.g. Fan et al. 2002; Kogut et al. there may have been early contributions from Popula- 2003). Planck Collaboration et al.(2020) recently in- ferred τ = 0.054 ± 0.007 from measurements of the

arXiv:2109.10375v1 [astro-ph.CO] 21 Sep 2021 tion III stars (e.g. Wu et al. 2021) and late contributions from active galactic nuclei (AGN; e.g. Madau & Haardt CMB temperature and polarization angular power spec- 2015). On large scales, higher-density regions near radi- tra, implying a late reionization midpoint at redshift ation sources are generally reionized earlier than lower- z ≈ 7.7 ± 0.6 (e.g. Glazer et al. 2018). There are now density regions far from sources in this inhomogeneous multiple evidence of dark Lyα troughs extending down process (e.g. Barkana & Loeb 2004; Trac et al. 2008). to z ≈ 5.5 in the spectra of high-redshift quasars (e.g. Becker et al. 2015; Bosman et al. 2018; Eilers et al. 2018; Bosman et al. 2021). This suggests that reionization Corresponding author: Hy Trac could have ended at z < 6 (e.g. Keating et al. 2020; [email protected] Nasir & D’Aloisio 2020), later than previously assumed. MNRAS 000,1–11 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Polarization signatures of the head-on collision model for Type Ia supernovae: How much asymmetry is too much?

Ran Livneh,1★ Boaz Katz1 1Department of Particle Physics & Astrophysics, The Weizmann Institute of Science, Rehovot 76100, Israel

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT

In a previous paper, we showed that the asymmetric ejecta produced by (zero impact parameter) head-on collisions of carbon- oxygen white dwarfs allow these progenitor models for Type Ia supernovae (SNe Ia) to cover the observed two-dimensional (2D) distribution of Si ii line depths (Branch plot). In this paper, we study the polarization signature associated with the 2D asymmetric ejecta of the collision model and a double-detonation model using similar Tardis radiative transfer simulations along different lines of sight with a spherical photosphere, combined with a new 3D Monte Carlo polarization code. We show that the polarization 푄 can be parametrized as a product 푄 = 푄max푄x of a radial structure component 푄max which is insensitive to the model specifics and is shown to be universally around 푄max ∼ 5%, and a cancellation component 푄x which depends on the asymmetry details. The continuum polarization is found to be low for both the collision and double-detonation models with 푄 ∼ 0.5%. However, the irregular Si distribution in the 2D head-on collision model results in Si ii line polarization reaching 푄 ∼ 3% (푄x . 50%) in tension with observations (mostly . 1.2%). In contrast, we show that the double-detonation model also covers the Branch plot, and yet results in low line polarization 푄 . 0.7% (푄x ∼ 10%) consistent with previous results and most SNe Ia. These results strengthen the case for asymmetric explosions as progenitors of SNe Ia, emphasizing an additional requirement for large polarization cancellations to account for the low observed line polarizations. Key words: supernovae: general – radiative transfer – polarization

1 INTRODUCTION Deviations from linear relations in the Stokes 푄 − 푈 plane as a func- tion of wavelength have also been observed (Cikota et al. 2019), There is strong evidence that Type Ia supernovae (SNe Ia) are the evidence that at least in some cases the SN ejecta does not exhibit product of thermonuclear explosions of white dwarfs (WDs), yet the axis-symmetry. Additionally, some correlations of Si ii line polariza- nature of the progenitor systems and the mechanism that triggers the tion with Δ푚 and with the pEW and velocity of the Si ii features explosion remain long-standing open questions (see e.g. Maoz et al. 15 have been reported (Wang et al. 2007; Meng et al. 2017; Cikota et al. 2014; Livio & Mazzali 2018; Soker 2019). 2019). The geometry of a SN explosion can help discriminate between various proposed progenitor systems. This geometry cannot be di- Several studies regarding polarization in SNe Ia have been pub- rectly observed but leaves traces in the polarization signature that lished. These range from the analysis of simple synthetic models can be sensed and analyzed. Thus, modeling the continuum and line of high-velocity feature polarization (Kasen et al. 2003) to the addition of polarization modeling to radiative transfer simulators

arXiv:2109.10371v1 [astro-ph.HE] 21 Sep 2021 polarization of SN spectra can be an important tool in the quest for understanding SNe Ia progenitors. such as Artis (Bulla et al. 2015). In Bulla et al.(2016b) a sub- Comparatively little observational polarization data exists, as spec- Chandrasekhar mass double-detonation model (Fink et al. 2010) and tropolarimetric observations consume significantly more resources a Chandrasekhar mass delayed-detonation model (Seitenzahl et al. than regular spectroscopic observations. Within the existing data, 2013) were simulated and shown to produce low levels of Si ii line observed spectra show low continuum polarization of the order of polarization (. 1%) consistent with most of the observed SNe Ia 0.3% and as high as 0.61% (Cikota et al. 2019). These low polariza- (. 1.2%), but inconsistent with higher polarization observations tion values have led to estimated ejecta asymmetries of 10% − 20% such as SN 2004dt (∼ 1.6%). In Bulla et al.(2016a) the violent (e.g. Hoflich 1991; Kasen et al. 2003). More significant polarization merger of a 1.1 and 0.9 M WD binary system (Pakmor et al. 2012) of up to 2% is observed in strong spectral lines such as Silicon, was simulated, resulting in Si ii line polarizations of ∼ 0.5% − 3.2% Calcium, Sulfur, and Magnesium (Wang & Wheeler 2008), suggest- over 35 viewing angles. This was concluded to be too high to explain ing that asymmetries exist in the element distribution in the ejecta. the low polarization levels commonly observed in normal Type Ia supernovae. This last result may seem to imply that SNe Ia progenitor mod- ★ E-mail: [email protected] els displaying significant asymmetry should not be considered as

© 2021 The Authors The Astrophysical Journal Letters, 916:L24 (7pp), 2021 August 1 https://doi.org/10.3847/2041-8213/ac1357 © 2021. The American Astronomical Society. All rights reserved.

Void Galaxy Distribution: A Challenge for ΛCDM

Saeed Tavasoli Physics Department, Kharazmi University, Tehran, Iran; [email protected] Received 2021 April 21; revised 2021 July 7; accepted 2021 July 9; published 2021 August 4

Abstract We extract void catalogs from the Sloan Digital Sky Survey Data Release 16 (SDSS DR16) survey and also from the Millennium simulation. We focus our comparison on distribution of galaxies brighter than Mr < −18 inside voids and study the mean separation of void galaxies, distance from the void center, and the radial density profile. We find that mean separation of void galaxies depends on void size, as bigger voids have lower mean separation in both samples. However, void galaxies in the observation sample seem to have generally larger mean–distance than simulated ones at any given void size. In addition, observed void galaxies tend to reside closer to the void center than those in the simulation. This discrepancy is also shown in the density profile of voids. Regardless of the void size, the central densities of real void profiles are higher than the ones in the predicted simulated catalog. Unified Astronomy Thesaurus concepts: Voids (1779)

1. Introduction gas-rich (Kreckel et al. 2012), and have lower gas-phase ( The internal structure of deepest under-densities of the metallicities than those in denser regions Mouhcine et al. 2007; Nicholls et al. 2014; Kreckel et al. 2015; Douglass & Cosmic Web provides an ideal testing ground for the theories ) fi of cosmological structure formation. Void environments as the Vogeley 2017 . However, there is no signi cant difference largest observable structures were first reported by Jöeveer between the faint end slopes in the luminosity function between ( ) et al. (1978), Tarenghi et al. (1978), Tifft & Gregory (1978), void and wall galaxies Croton et al. 2005 . Observational and Tully & Fisher (1978) in the three-dimensional (3D) studies demonstrate that statistically void galaxies have stellar distribution of a spectroscopic galaxy sample. Their sizes span disks with smaller radii compared to a control sample of late- tens to hundreds of Mpc, and close to 60% of the volume of the type galaxies (van de Weygaert & Platen 2011). Rojas et al. universe is occupied by these huge structures (da Costa et al. (2005) indicated that void galaxies are forming stars at a high 1988; Geller & Huchra 1989), yet they contain only about 7% rate and are fainter compared to the ones in the cluster regions. of the observed galaxies in the Universe (Pan et al. 2012). However, the existence of low-mass galaxies in voids seems These under-dense regions that exist between the filaments, to be in conflict with the prediction of ΛCDM models (e.g., cluster, and walls of galaxies are powerful cosmological Dekel & Silk 1986; Hoffman et al. 1992; Tikhonov & laboratories and have been widely used in recent years to Klypin 2009; Tavasoli et al. 2013). This problem was termed extract cosmological information (see, e.g., Pisani et al. 2015; the “void phenomenon” by Peebles (2001). Hamaus et al. 2016; Planck Collaboration et al. 2016, and In this Letter we study galaxy distributions inside voids to references therein). They are expected to be less evolved and discern whether or not void environment has an influence on therefore better preserve the memory of the initial Universe. the geometrical position of void galaxies. We attempt to For instance, because voids are in the linear regime, this makes address this issue by calculating the distance of void galaxies them attractive to theoretical modeling and tests. Also, several from center of the host void and also by finding the minimum key results such as measurements of the velocity and density distance of void galaxies in each void separately. To achieve profiles of cosmic voids (Paz et al. 2013; Hamaus et al. 2016), this we apply the Minimal Spanning Tree (MST) algorithm on the void auto-correlation function and clustering bias (Clampitt a 3D void catalog extracted from observations and simulations. et al. 2016), the depth, abundance, and general properties of The observational measurements of galaxy position inside ( ) voids Sutter et al. 2012; Nadathur et al. 2019 , and the Alcock- voids can be compared to the numerical orsemi-analytic ( Paczynski test Sutter et al. 2012; Hamaus et al. 2016; Mao simulations to test the ΛCDM model. We study this problem ) et al. 2017 , have been obtained from these observations and by analyzing voids in the Sloan Digital Sky Survey Data simulation samples. Release 16 (SDSS DR16; Ahumada et al. 2020) and by On the other hand, the low frequency of galaxy mergers carrying out a parallel and comparative analysis on the inside these regions makes them an ideal environment for the ( ) / Millennium simulation Springel et al. 2005 . study of the environmental dependence of galaxy formation The outline of this Letter is as follows. In Section 2,we evolution with the minimal impact from mergers (see van de briefly describe the simulation and observation sample that we Weygaert & Platen 2011, for a review). Observationally, void use in this Letter and the method that we apply to find voids. In galaxies have shown different properties than galaxies in more Section 3, we introduce two parameters to investigate the average or denser regions of the Universe. They are expected to ( distribution of void galaxies in observation and simulation have lower stellar mass Croton et al. 2005; Hoyle et al. 2005; fi Moorman et al. 2015), tend to be bluer and have a later sample and compare the density pro le of void galaxies in both morphological type (Grogin & Geller 2000; Rojas et al. 2004; sample. We present our conclusions in Section 4. Patiri et al. 2006; Park et al. 2007; von Benda-Beckmann & Müller 2008; Hoyle et al. 2012; Tavasoli et al. 2015), are more

1 ICRC 2021 THE ASTROPARTICLE PHYSICS CONFERENCE ONLINE ICRC 2021Berlin | Germany THE ASTROPARTICLE PHYSICS CONFERENCE th Berlin37 International| Germany Cosmic Ray Conference 12–23 July 2021

Study of 26Al in the COSI 2016 superpressure balloon flight

Jacqueline Beechert푎,∗ on behalf of the COSI Collaboration (a complete list of authors can be found at the end of the proceedings) 푎Space Sciences Laboratory, University of California, Berkeley, 7 Gauss Way, Berkeley, CA 94720, USA E-mail: [email protected]

The Compton Spectrometer and Imager (COSI) is a balloon-borne compact Compton telescope designed to survey the 훾-ray sky from 0.2 to 5 MeV. COSI’s wide field-of-view (FOV) and excellent energy resolution from high-purity germanium detectors make it uniquely capable of probing this under-explored energy regime. In particular, it can facilitate understanding of stellar nucleosynthesis through studies of diffuse emission from the radioisotope 26Al at 1.809 MeV. In 2016, COSI was launched from Wanaka, New Zealand on a NASA superpressure balloon and flew for 46 days. The flight was a technologic and scientific success, boasting live detection and polarization studies of GRB160530A, spectral analysis of the Crab Nebula and the 511-keV positron annihilation emission at the Galactic Center, and detection of Cygnus X-1. This article details the first maximum-likelihood search for the 1.809 MeV signature of Galactic 26Al in the 2016 data. The analysis reveals a promising excess around the expected energies of an 26Al signature with 3.7휎 significance and a measured flux of (17.0 ± 4.9) × 10−4 ph cm−2 s−1. Further exploration is currently underway to solidify the measurement. arXiv:2109.10365v1 [astro-ph.HE] 21 Sep 2021

37th International Cosmic Ray Conference (ICRC 2021) July 12th – 23rd, 2021 Online – Berlin, Germany

∗Presenter

© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). https://pos.sissa.it/

26Al in the COSI 2016 SPB flight Jacqueline Beechert

1. Introduction

Aluminum-26 (26Al) is a radioactive isotope that decays to an excited state of Magnesium-26 (26Mg) with a half-life time of 0.7 Myr. The excited 26Mg de-excites to its ground state, emitting a 1.809 MeV 훾-ray. Because of its long-lived nature, 26Al is likely to decay after ejection from its stellar production sites into the interstellar medium (ISM), facilitating both studies of the stellar conditions responsible for nucleosynthesis and the hot phase of the ISM. Hence, 26Al traces the life-cycle of elements in the Milky Way from their initial synthesis to their feedback into the ISM. The 1.809 MeV line was first detected in the Milky Way by the High Energy Astronomy Observatory (HEAO-3) satellite in 1984 [1]. The Compton Telescope (COMPTEL) on board the Compton Gamma-Ray Observatory (CGRO) advanced studies of 26Al in the 1990s by mapping the diffuse emission [2], and the spectrometer SPI on board the International Gamma-ray Astrophysics Laboratory (INTEGRAL) satellite is currently measuring the signal. SPI has mapped 26Al [3] and detected the line at 1809.83 ± 0.04 keV. The full-sky flux was found to be (1.7 ± 0.1) × 10−3 ph cm−2 s−1 [4]. Kinematic measurements reveal a systematic shift of the line with Galactic longitude, implying a velocity of 26Al emission up to 200 km s−1 greater than that expected from Galactic rotation [5,6]. Ejection of stellar material carrying 26Al into merging, low-density HI supershells leading the spiral arms of the Galaxy has been proposed as an explanation of the observed velocity [7]. A complete review of current understanding of 26Al, with an emphasis on the contributions from SPI, is provided in [8]. The Compton Spectrometer and Imager (COSI) is a balloon-borne compact Compton telescope with excellent spectral resolution of 0.24% FWHM at 1.8 MeV, making it uniquely capable of new measurements required to resolve the primary sources of 26Al emission and its distribution throughout the ISM. In this proceedings, we search for a signal of 26Al in the COSI 2016 flight data and compare the results to expectations from previous studies of 26Al. Section2 provides an overview of the COSI 2016 flight, Section3 details the analysis, and Section4 discusses the results. This study establishes the scientific potency of modern Compton telescopes in nucleosynthesis studies and is key proof-of-concept for such studies on the Compton telescope spacecraft mission, COSI-SMEX, currently in competitive NASA Phase A development [9].

2. The Compton Spectrometer and Imager

COSI was launched as a science payload on a NASA superpressure balloon from Wanaka, New Zealand on May 17, 2016. Twelve high-purity cross-strip Ge detectors (each 8 × 8 × 1.5 cm3) arranged in a 2 × 2 × 3 array comprise a 972 cm3 active volume within which incident photons of energies 0.2–5 MeV are measured. The energy and three-dimensional position of each interaction are recorded, enabling reconstruction of the photon path and a localization of the incident photon direction to a circle on the sky. Detailed explanations of calibrations, event reconstruction, and imaging algorithms employed in Compton telescopes are provided in [10] and [11]. COSI has a wide FOV of ∼1휋 sr defined by the five anti-coincidence CsI shields that surround the four sides and bottom of the detector array. It is always pointed at zenith and sweeps the sky through the Earth’s rotation during flight. The shields actively veto the dominant 훾-ray and cosmic-ray background incident from below the instrument, namely that of atmospheric Earth albedo radiation.

2 Axion-like particle production in neutron star mergers: gamma-ray signal detectability

Damiano F. G. Fiorillo1, 2 and Fabio Iocco1, 2 1Dipartimento di Fisica ”Ettore Pancini”, Universit`adegli studi di Napoli “Federico II”, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy 2INFN - Sezione di Napoli, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy (Dated: September 23, 2021) Axion-like particles can in principle be produced in the dense environment of neutron star merg- ers and core-collapse supernovae, and converted in gamma–rays in the magnetic field intervening between the event and the Earth. While this process has been studied for supernovae, there is no estimate of the gamma-ray signal expected from neutron star mergers. In this work we explore this new, exotic signature of the merger, and find that for a large region of the ALP parameter space its magnitude is comparable with that of a signal of similar nature from a SN. We show detection forecasts for these events, finding that they could be observable with the e-ASTROGAM telescope, thus opening a new window into both the astrophysics of these cataclysmic events, and of new particles beyond the standard model.

I. INTRODUCTION we find that SNe are still a brighter source than NSMs in terms of ALP production: for this reason, we also Axion-like particles (ALPs) can in principle be pro- determine the projected constraints (see Fig.3) that e- duced in the hot environments of extreme astrophysical ASTROGAM could obtain from the non-observation of sources, such as magnetars [1,2], pulsars [3], neutron gamma-rays from a future extragalactic SN. stars [4], compact binaries [5], black holes [6], gamma- ray bursts [7–10], active galactic nuclei [11–14], and su- II. ALP FLUXES FROM SN AND NSM pernovae (SNe) [15–33]. Such extreme environments now also list neutron star We consider ALPs produced via nucleon-nucleon ALP mergers (NSMs). The actual observation of a NSM bremsstrahlung, a process has been thoroughly studied in through gravitational waves [34] produced in the event the literature [26, 46, 47, 51]. Since strong interactions in has prompted the exploration of a new realm of densi- dense nucleon environments are poorly known, they are ties and temperatures, as addressed in [35–43]. In such often treated at the level of one-pion exchange (OPE): unique conditions, the bremsstrahlung process of neu- with this approach, Ref. [51] determines the ALP emis- trons N + N → N + N + a could be an efficient way sivity in the hypermassive neutron star which forms for of producing ALPs, a process that has been studied in a short time at the end of a NSM. Recently, Ref. [26] has the context of core–collapse SNe [23, 26, 44–47], lead- shown that going beyond the OPE approximation can ing to constraints on the ALPs properties from the non- lead to a reduction by roughly an order of magnitude in observation from the SN1987A [15, 18, 19, 21, 48, 49], the ALP emissivity. However, for this work we adopt the and to forecasts on future SN detections [30, 31, 33]. emissivity from [51]: as discussed in Sec.IV, this proce- The possibility that ALP production could take place dure does not change the qualitative conclusions of our in NSMs has been addressed in [50, 51], where the authors work1. We use the emissivity for the full phase-space cal- study the impact of ALPs streaming through the dense culation for non-relativistic neutrons, assuming a density merger, provoking an additional cooling of the material equal to the saturation density of neutrons. For relativis- and modifying its structure. This suggests the possibil- tic neutrons and for higher densities the emissivity would ity of constraining ALP couplings by imposing the limi- be larger, so that our choice is conservative. tations arising from current observations of gravitational Furthermore, Ref. [51] only considers the ALP pro- arXiv:2109.10364v1 [astro-ph.HE] 21 Sep 2021 waves. duction due to neutron-neutron ALP bremsstrahlung. Here we ask whether there is a region of the parameter space –still allowed by observations– where ALP emission is sizable, and if so whether any observational signature 1 A possible discrepancy between that adopted here and the one in would be left by ALPs being converted into photons in Ref. [26] would not affect our conclusions on the detectability. In the magnetic field intervening between us and the merger. fact, the normalization of the ALP flux depends on the coupling We find that for a sizable region in the ALP parameter Gan, for which the main constraints come from the ALP emission space the merger structure is preserved, and a signal in in SN1987A (either because of cooling or because of ALP to gamma conversion). Therefore a discrepancy in the emissivity is soft gamma–rays of energy is potentially observable with reabsorbed in a redefinition of the coupling Gan: in other words, the future gamma-ray detector e-ASTROGAM [52] (see if the emissivity were 10 times lower, the constraints on Gan Fig.2), thus opening a new exciting window in the multi- would be a factor of 3 weaker, and therefore the signal we obtain messenger study of NSMs. Due to their longer duration, would still be allowed by the present parameter space. Chemical Abundances in Gaseous Nebulae: from the Milky Way to the early universe AAA Workshop Series nnn, 2021 M. Cardaci & G. Hägele, eds.

Contributed Paper

Self-consistent population spectral synthesis with FADO: mean stellar metallicity of galaxies in spectral synthesis methods

Ciro Pappalardo1,2 1Instituto de Astrofísica e Ciências do Espaço, Universidade de Lisboa - OAL, Tapada da Ajuda, PT1349-018 Lisboa, Portugal 2Departamento de Física, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande, PT1749-016 Lisboa, Portugal

Abstract. In this work, we investigate the reliability of spectral syn- thesis methods in the estimation of the mean stellar age and metallicity, addressing the question of which signal-to-noise ratios (S/N) are needed to determine these quantities and how these depend on the tool used to model the spectra. To address this problem we used realistic simulated spectra containing stellar and nebular emission, reproducing the evolution of a galaxy for a constant and exponentially declining star formation law. The spectra have been degraded to different S/N and analysed with three different spectral synthesis codes: Fado, Steckmap, and Starlight assuming similar fitting set-ups and the same spectral bases. For S/N ≤ 5 all tools considered show a large diversity in the results. Fado and Starlight find median differences in light-weighted mean stel- lar ages of ∼0.1 dex, while Steckmap shows a higher value of ∼0.2 dex. For S/N > 50 the median differences in Fado are ∼0.03 dex (∼ 7%), a factor 3 and 4 lower than the 0.08 dex (∼20%) and 0.11 dex (∼30%) obtained from Starlight and Steckmap, respectively. Our results indicate that phases of high specific star formation rate (sSFR) in galaxies require analysis tools that do not neglect the nebular contin- uum emission in the fitting process, since purely stellar models would have strong problems in the estimation of star formation history, even in presence of high S/N spectra. The median value of these differences are of the order of 7% (Fado), 20% (Starlight), and 30% (Steckmap) for arXiv:2109.10363v1 [astro-ph.GA] 21 Sep 2021 light-weighted quantities, and 20% (Fado), 60% (Starlight), and 20% (Steckmap) for mass-weighted quantities. That implies a severe overes- timation of the mass-to-light ratio and stellar mass, even in the presence of a mild contribution from the nebular continuum. Our work under- lines once more the importance of a self-consistent treatment of nebular emission, which is the only route towards a reliable determination of the assembly of any high-sSFR galaxy at high and low redshift.

Key words: galaxies: abundances — galaxies: ISM — techniques: spec- troscopy 1 MNRAS 000,1–21 (2015) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Combined analysis of neutron star natal kicks using proper motions and parallax measurements for radio pulsars and Be X-ray binaries

Andrei P. Igoshev,1★ Martyna Chruslinska,2 Andris Dorozsmai, 3 and Silvia Toonen 3,4 1Department of Applied Mathematics, University of Leeds, LS2 9JT Leeds, UK 2Institute of Mathematics, Astrophysics and Particle Physics, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands 3Institute of Gravitational Wave Astronomy, School of Physics and Astronomy, University of Birmingham, Birmingham, B15 2TT, United Kingdom 4Anton Pannekoek Institute for Astronomy, University of Amsterdam, 1090 GE Amsterdam, The Netherlands

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Supernova explosion and the associated neutron star natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain neutron star natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries. Our best model for natal kicks is consistent with both measurements and contains a fraction of 푤 = 0.2±0.1 weak natal kicks +25 with 휎1 = 45−15 km/s, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: 휎2 = 336 km/s. We found that currently used models for natal kicks of neutron stars produced by electron capture supernova (combination of maxwellian 휎 = 265 km/s and 휎 = 30 km/s for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and Be X-ray binaries. Key words: binaries: general – binaries: close – pulsars: general – stars: neutron – stars: massive – methods: statistical – X-rays: binaries

1 INTRODUCTION quently disrupt the binary and (2) different formation channels could lead to different velocity distributions. In studies of observed veloc- Neutron stars (NSs) often receive a significant birth velocity (natal ity distribution of radio pulsars e.g. Hobbs et al.(2005); Verbunt kick) at the moment of supernova explosion (Lyne & Lorimer 1994). et al.(2017); Igoshev(2020) young pulsar velocities were assumed On one hand, this kick is frequently strong with amplitude |푣| > to be equal to the isolated NS natal kick velocities and the influence 100 km s−1 (in some cases reaching 800 km s−1) while typical of binarity of the NS progenitors was never analysed. velocity dispersion of NS progenitors are < 20 km s−1. On the other Individual velocities of radio pulsars (mode of activity typical hand, formation of certain Galactic double neutron stars (DNSs) arXiv:2109.10362v1 [astro-ph.HE] 21 Sep 2021 for young, magnetised NSs with fast rotation, see e.g. Deller et al. systems and presence of pulsars in globular clusters (with escape 2019) are estimated based on VLBI measurements of parallaxes and velocities .50 km/s) clearly require small natal kicks (Tauris et al. proper motions. Individual velocities for high-mass X-ray binaries 2017). It is unclear whether the NS natal kick velocity distribution (HMXBs) and double neutron stars are estimated based on analysis is broad and extends from ≈ 0 – 1000 km s−1, or if NSs with some of orbital properties (see e.g. Dewi et al. 2005; Beniamini & Piran particular characteristics are formed with much weaker natal kicks 2016 and Tauris et al. 2017). (0 − 30 km s−1). On the one hand, young radio pulsars with reliably measured The problem of estimating the NS natal kick distribution has parallaxes and proper motions seem to be isolated (Igoshev & Perets two parts, which were not sufficiently addressed in the previous 2019). Also, recent analysis of all radio pulsars recorded in the research on the subject: (1) NSs receiving weak natal kick should ATNF pulsar catalogue1 (Manchester et al. 2005) found an upper be less common among isolated radio pulsars because they less fre-

★ E-mail: [email protected] 1 ATNF pulsar catalogue http://www.atnf.csiro.au/research/pulsar/psrcat

© 2015 The Authors arXiv:2109.10361v1 [astro-ph.SR] 21 Sep 2021 ★ hnpeitdb ubln ovcini h ht dwarfs’ white the in convection com turbulent dwarf by white predicted with than (eMSPs) pulsars millisecond Eccentric ABSTRACT ZZZ form original in YYY; Received XXX. Accepted ceti ilscn usr yrsnn convection resonant by pulsars millisecond Eccentric edsuiso tr hwsm vdnefrepcal orde especially for evidence some forci show tidal stars the of also studies is field which period, so spin flows giant’s convective red by locked driven resonantly are eccentricities † MNRAS nse ta.2015 al. et Knispel 2010 Bailes 1995 al. et Rappaport © a scenario overflow companion Roche-lobe differ basic dwarf dramatic same describ the the white as eccentricity, despite relation Their that suggesting period d. MSPs, mass–orbital 30 ventional same – the an 20 along 0.1 of fall – periods 0.03 orbital of eccentricities ilar similar have (eMSPs) pulsars 10 from range typically epicycli the ( between eddies convective the reached and is perturbs equipartition stochastically energy until field quadrupole fluctuating eddie measurable convective ant’s but from forcing small gravitational a by to maintained down eccentricity orbital the damp agrecnrcte aebe icvrd( discovered been have eccentricities pulsars larger millisecond field several however, years, recent In re nt.W hwta eoatcneto osseccent boosts convection resonant that show We unity. order e words: Key eas hwhwvrain nteed unvrtm arising resonance, time the turnover of eddy width the effective in variations how show also We austa hrceiecnetoa Ss ngo agreem good in MSPs, conventional characterize that values eidadtems ftehlu ht wr httegatle giant the that dwarf white the helium ( the correlates behind of mass condition the filling and rel th Roche-lobe period Combining a the radius. prescribes giant’s with core the lation and helium mass bu inert core’s shell giant’s the Nuclear between red star. the ma neutron transfers a around and to ing lobe momentum Roche angular its though spin overflows giant are and red companions a dwarf when white form with pulsars Millisecond INTRODUCTION 1 % 2 1 ia Ginzburg Sivan eateto at n lntr cec,Uiest fCa of University Science, Planetary and Earth of Department eateto srnm,Uiest fClfri,Berkel California, of University Astronomy, of Department 1Pgs Fellow. b Pegasi 51 01TeAuthors The 2021 -al [email protected] E-mail: h togtdsta h eto trrie nteconvecti the on raises star neutron the that tides strong The = 0–3 ,rmral ls oterdgatcnetv dytu eddy convective giant red the to close remarkably d, 30 – 20 ose l 1987 al. et Joss 000 ; , are l 2013 al. et Barr 1 – usr:gnrl–wiedaf iais eea convec – general binaries: – dwarfs white – general pulsars: 6 (2021) ; ca ta.2018 al. et Octau ; 1 ars&Svne1999 Savonije & Tauris ★ − ; 7 † aoij 1987 Savonije hne 1992 Phinney o10 to n ueeChiang Eugene and ; eeae l 2013 al. et Deneva − 3 eedn nteobtlperiod. orbital the on depending , Δ .Teeecnrcmillisecond eccentric These ). % ; .Tersda eccentricities residual The ). / hne ukri1994 Kulkarni & Phinney % ≈ hmine l 2008 al. et Champion ). 0 ; . aioe l 2015 al. et Camilo 4. y A97031,USA 94720-3411, CA ey, 1 , 2 iona ekly A97046,USA 94720-4767, CA Berkeley, lifornia, ihmuch with .Tegi- The s. motion c h orbit the pisto pplies ne in ences egiant ve scon- es sim- d orbital rpit2 etme 01Cmie sn NA L MNRAS using Compiled 2021 September 23 Preprint value sre- is ation aves to t gpro.Nmrclsmltoso oaigrdgat and giants red rotating of simulations Numerical period. ng rn- ss eo aechrn hntetroe iemthstetidal the matches time turnover the when coherent made mehow n ihosrain ( observations with ent s e in rgntr.Teobtlproso MP cluster eMSPs of periods orbital The progenitors. giant red e o atrswe h ovcieRsb number Rossby convective the when patterns flow red ; ; ; iiisb atrof factor a by ricities ain xii ria cetiiisodr fmagnitud of orders eccentricities orbital exhibit panions rmrdgatmtliiyvrain a erdc h obs the reproduce can variations metallicity giant red from ohppltos(see populations both idyecnrcPRJ1146 PSR eccentric mildly 1 o eccentrici formation system’s the the excite to of can lead top which and on disc gas loss circumbinary residual mass of drive burning dwarf unstable proto-white to due flashes clear 2019 ( evolution different progenitors qualitatively massive a re in transfer more mass unstable had of phase these a undergone have dwarfs; may which white ONeMg or CO ( e eccentricity have may interactions dynamical where clusters, globular 2021 a ern submission, nearing was ha envelope extended companion’s pated. the after presumably delayed, avo until be To rotation, must eccentricity. rele collapse orbital sudden the large the circularization, a scenario, induces this energy In binding pulsar. millisecond binar a close a form in dwarf white (super-Chandrasekhar) massive Tauris & Freire hs rniinfo eto trt tag ur star quark strange a to star neutron a also from transition phase a eso ihbt fteeseais hc rdc ahrs rather predict which scenarios, ( masses these pulsar of both with tension ain n steeoeacerecpin( exception clear a therefore is and panion, eccentricity w de h aafrom data the added (we eecuePRJ9302,wihhsaslrms main-sequ solar-mass a has which J1903+0327, PSR exclude We eea da aebe u ot oepaneMSPs. explain to forth put been have ideas Several epo h bevdssesi Fig. in systems observed the plot We nvrtime rnover in ta.2021 al. et Jiang ; therein). references and , h ta.2019 al. et Zhu in tal. et Jiang tion 4 = ai ege1995 Heggie & Rasio noidse l 2016 al. et Antoniadis 7 ( . 2014 4 C C eddy ( ( nuc × C 2015 nuc 10 .Sm MPmse,hwvr emt ein be to seem however, masses, eMSP Some ). ). oie tal. et Lorimer oie ta.2021 al. et Lorimer a of collapse accretion-induced that argued ) epooeta h nmlul large anomalously the that propose We . stegatsncerbrigtime-scale). burning nuclear giant’s the is u ta.2018 al. et Hui − / Antoniadis 3 rpsdasmlrmcaiminvolving mechanism similar a proposed ) % − hc sjs eo u 10 our below just is which , 60wt nobtlperiod orbital an with 6610 ) 1 / 2 1 ≈ 3 ( 2021 ( with systems exclude also We ). hmine l 2008 al. et Champion 2014 × hne ukri1994 Kulkarni & Phinney ; ; are l 2017 al. et Barr 10 none h icvr fthe of discovery the announced ) a i2021 Li & Han oFig. to ugse htthermonu- that suggested ) 3 1 vrterandom-walk the over eecueplasin pulsars exclude We . 1 ). − A % T 2 E MPthreshold eMSP = tl l v3.0 file style X ; % .A hsletter this As ). ; tvl tal. et Stovall 62 / in tal. et Jiang C magnetic . eddy ctdthe xcited n an and d 8 necom- ence yrapid by larger e around could y dissi- s dtidal id y The ty. sulting pecific s of ase erved ). sof is (see the a f ly Robust marginalization of baryonic effects for cosmological inference at the field level

Francisco Villaescusa-Navarro,1, 2, ∗ Shy Genel,2, 3 Daniel Angl´es-Alc´azar,4, 2 David N. Spergel,2, 1 Yin Li,2 Benjamin Wandelt,5, 2 Leander Thiele,6 Andrina Nicola,1 Jose Manuel Zorrilla Matilla,1 Helen Shao,1 Sultan Hassan,2, 7 Desika Narayanan,8, 9 Romeel Dave,10, 7, 11 and Mark Vogelsberger12 1Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton NJ 08544, USA 2Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, New York, NY, 10010, USA 3Columbia Astrophysics Laboratory, Columbia University, New York, NY, 10027, USA 4Department of Physics, University of Connecticut, 196 Auditorium Road, Storrs, CT, 06269, USA 5Sorbonne Universite, CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis boulevard Arago, 75014 Paris, France 6Department of Physics, Princeton University, Princeton, NJ 08544, USA 7Department of Physics & Astronomy, University of the Western Cape, Cape Town 7535, South Africa 8Department of Astronomy, University of Florida, Gainesville, FL, USA 9University of Florida Informatics Institute, 432 Newell Drive, CISE Bldg E251, Gainesville, FL, USA 10Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 11South African Astronomical Observatories, Observatory, Cape Town 7925, South Africa 12Kavli Institute for Astrophysics and Space Research, Department of Physics, MIT, Cambridge, MA 02139, USA (Dated: September 23, 2021) We train neural networks to perform likelihood-free inference from (25 h−1Mpc)2 2D maps con- taining the total mass surface density from thousands of hydrodynamic simulations of the CAMELS project. We show that the networks can extract information beyond one-point functions and power −1 spectra from all resolved scales (& 100 h kpc) while performing a robust marginalization over baryonic physics at the field level: the model can infer the value of Ωm(±4%) and σ8(±2.5%) from simulations completely different to the ones used to train it.

The wealth and quality of data from upcoming surveys training or testing the model on a particular type of sim- such as DESI [1], eRosita [2], Euclid [3], PFS [4], SKA [5], ulation. We will show that neural networks trained on Roman [6], and Rubin [7] is expected to bring cosmology total matter surface density maps from one suite of simu- to a different level. The data from these surveys contain lations can infer the correct cosmology even when tested the key to answering fundamental questions in physics on maps from another suite, and vice versa. It is not obvi- such as what the nature of dark energy or the sum of ous that this should work, because the two suites follow a the neutrino masses are. The main obstacle to achieving different set of hydrodynamic equations, use distinct nu- this is not the amount and quality of data, but two sig- merical methods to solve them, and employ different sub- nificant theoretical challenges: 1) the optimal summary grid models. We also show that rather than simply ex- statistic(s) needed to extract the maximum amount of tracting cosmological information only from large scales information from data in the mildly and fully non-linear to avoid bias, our networks successfully extract cosmo- regimes is unknown, and 2) astrophysical processes such logical information from the smallest scales resolved in as feedback from supernovae and active galactic nuclei the maps. These results have important consequences for (AGN) are expected to impact the non-linear regime in weak lensing surveys, paving the way to a method that a poorly understood manner. can extract all information from the field while marginal- Ideally, one would want to use an estimator that is able izing over baryonic effects. to extract the maximum amount of information directly We use data from the Cosmology and Astrophysics from the field while marginalizing over baryonic effects at with MachinE Learning Simulations (CAMELS) project the same time. In [8] we demonstrated with toy exam- [10] at z = 0. CAMELS contains two simulation suites, ples that neural networks can achieve both goals. In our IllustrisTNG and SIMBA, which have been run using the arXiv:2109.10360v1 [astro-ph.CO] 21 Sep 2021 companion paper [9] we have shown that neural networks AREPO [11] and GIZMO [12] codes, employing the sub- can extract cosmological information embedded in state- grid physics implemented in the IllustrisTNG [13, 14] and of-the-art hydrodynamic simulations from complex 2D SIMBA [15] simulations, respectively. For each suite, we fields such as gas temperature, electron number density, use the latin-hypercube (LH) set, which contains 1,000 and magnetic fields, while marginalizing over baryonic simulations, each with a different value of the cosmologi- effects at the same time. cal parameters Ωm and σ8, the fractional matter density In this letter we show, for the first time, that this pro- and the amplitude of fluctuations in spheres of radius 8 cedure can be robust. In other words, it does not rely on h−1Mpc, and four astrophysical parameters controlling the efficiency of supernova and AGN feedback (ASN1, ASN2, AAGN1, AAGN2) as well as a different value of the random seed used to generate the initial conditions. The ∗Electronic address: [email protected] value of the cosmological and astrophysical parameters MNRAS 000,1–15 (2021) Preprint 23 September 2021 Compiled using MNRAS LATEX style file v3.0

Formation and evolution of young massive clusters in galaxy mergers: a SMUGGLE view

Hui Li1,2★† , Mark Vogelsberger2 , Greg L. Bryan1,3 , Federico Marinacci4 , Laura V. Sales5 , Paul Torrey6 1Department of Astronomy, Columbia University, New York, NY 10027, USA 2Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 3Center for Computational Astrophysics, Flatiron Institute, 162 5th Ave, New York, NY 10003, USA 4Department of Physics & Astronomy “Augusto Righi”, University of Bologna, via Gobetti 93/2, 40129 Bologna, Italy 5Department of Physics & Astronomy, University of California, Riverside, 900 University Avenue, Riverside, CA 92521, USA 6Department of Astronomy, University of Florida, 211 Bryant Space Sciences Center, Gainesville, FL 32611, USA

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength 0.71 correlates positively with YMC mass, 휆tid ∝ 푀 , which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the center of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day. Key words: galaxies: interactions – galaxies: ISM – galaxies: structure – ISM: clouds – methods: numerical

1 INTRODUCTION Mortlock et al. 2013), trigger the formation of ultraluminous in- frared galaxies (e.g. Sanders et al. 1988; Melnick & Mirabel 1990), The interactions of dark matter halos over cosmic time are the key and enhance atomic/molecular gas contents (e.g. Braine & Combes physical process of hierarchical structure formation in the Λ-cold 1993; Ueda et al. 2014; Larson et al. 2016; Violino et al. 2018; Elli- dark matter cosmology (Somerville & Davé 2015). Mergers of halos son et al. 2018). Moreover, major mergers not only change the global arXiv:2109.10356v1 [astro-ph.GA] 21 Sep 2021 naturally lead to the mergers of their central galaxies and inevitably properties of the galaxies but also strongly affect various spatially- leave strong imprints on their evolution. resolved properties of the interstellar medium (ISM) and modes Over the last several decades, observations of nearby and high of star formation. Several high-resolution optical observations of redshift galaxies suggest that major mergers alter star formation nearby merger systems (e.g. Holtzman et al. 1992; Whitmore & activities (e.g. Ellison et al. 2008; Wong et al. 2011; Patton et al. Schweizer 1995; Holtzman et al. 1996; Schweizer & Seitzer 1998; 2013; Knapen et al. 2015; Davies et al. 2015), modify the mor- Zepf et al. 1999; Forbes & Hau 2000; Schweizer 2004; Reines et al. phology of galaxies (e.g. Conselice et al. 2003; Lotz et al. 2008; 2008; Miah et al. 2015) have discovered a large number of young 4 massive clusters (YMCs) with masses larger than 10 M and ages ★ E-mail: [email protected] younger than 100 Myr (see a recent review by Portegies Zwart et al. † NHFP Hubble Fellow 2010a). The total number of YMCs are much larger and the cluster

© 2021 The Authors Draft version September 23, 2021 Typeset using LATEX twocolumn style in AASTeX63

Rapid stellar and binary population synthesis with COMPAS

Team COMPAS,1 Jeff Riley,2, 3 Poojan Agrawal,4, 3 Jim W. Barrett,5 Kristan N. K. Boyett,6 Floor S. Broekgaarden,7 Debatri Chattopadhyay,4, 3 Sebastian M. Gaebel,8 Fabian Gittins,9 Ryosuke Hirai,2, 3 George Howitt,10, 3 Stephen Justham,11, 12, 13 Lokesh Khandelwal,11 Floris Kummer,11 Mike Y. M. Lau,2, 3 Ilya Mandel,2, 3, 5 Selma E. de Mink,13, 11, 7 Coenraad Neijssel,5, 3 Tim Riley,2, 3 Lieke van Son,7, 11, 13 Simon Stevenson,4, 3 Alejandro Vigna-Gomez,´ 14, 15 Serena Vinciguerra,11 Tom Wagg,7, 13 Reinhold Willcox,2, 3

1The public COMPAS code is a product of work by the entire COMPAS collaboration over many years; we therefore kindly request that, in recognition of this team effort, the paper is cited as “Team COMPAS: J. Riley et al.” 2School of Physics and Astronomy, Monash University, Clayton, Victoria 3800, Australia 3OzGrav, Australian Research Council Centre of Excellence for Gravitational Wave Discovery, Australia 4Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, VIC 3122, Australia 5Institute of Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Birmingham, B15 2TT, United Kingdom 6Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK 7Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., Cambridge, MA 02138, USA 8Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Callinstrasse 38, D-30167 Hannover, Germany 9Mathematical Sciences and STAG Research Centre, University of Southampton, Southampton SO17 1BJ, UK 10School of Physics, University of Melbourne, Parkville, Victoria, 3010, Australia 11Anton Pannekoek Institute of Astronomy and GRAPPA, Science Park 904, University of Amsterdam, 1098XH Amsterdam, The Netherlands 12School of Astronomy & Space Science, University of the Chinese Academy of Sciences, Beijing 100012, China 13Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany 14DARK, Niels Bohr Institute, University of Copenhagen, Jagtvej 128, 2200, Copenhagen, Denmark 15Niels Bohr International Academy, The Niels Bohr Institute, Blegdamsvej 17, 2100 Copenhagen, Denmark

(Dated: September 23, 2021)

ABSTRACT COMPAS (Compact Object Mergers: Population Astrophysics and Statistics) is a public rapid binary pop- ulation synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrised assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github arXiv:2109.10352v1 [astro-ph.IM] 21 Sep 2021 repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00.

Keywords: stars: stellar evolution, stars: binaries, black holes, gravitational waves

jeff[email protected], fl[email protected], [email protected], [email protected] Draft version September 23, 2021 Typeset using LATEX default style in AASTeX631

Competitive X-ray and Optical Cooling in the Collisionless Shocks of WR 140

A. M. T. Pollock,1 M. F. Corcoran,2, 3 I. R. Stevens,4 C. M. P. Russell,3 K. Hamaguchi,2, 5 P. M. Williams,6 A. F. J. Moffat,7 G. Weigelt,8 V. Shenavrin,9 N. D. Richardson,10 D. Espinoza,3 and S. A. Drake2

1Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, UK 2CRESST and X-ray Astrophysics Laboratory, NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 3The Catholic University of America, 620 Michigan Ave., N.E. Washington, DC 20064, USA 4School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK 5Department of Physics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 6Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 7D´epartement de physique and Centre de Recherche en Astrophysique du Qu´ebec (CRAQ), Universit´ede Montr´eal, C.P. 6128, Succ. Centre-Ville, Montr´eal,Qu´ebec, H3C 3J7, Canada 8Max-Planck-Institut f¨urRadioastronomie, Auf dem H¨ugel69, 53121 Bonn, Germany 9Sternberg Astronomical Institute, Moscow State University, Universitetskii pr. 13, Moscow, 119992 Russia 10Embry-Riddle Aeronautical University, Physics and Astronomy Department, 3700 Willow Creek Road, Prescott, AZ 86301, USA

ABSTRACT WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally well- determined orbital and stellar parameters. Bright, variable X-ray emission is generated in shocks produced by the collision of the winds of the WC7pd+O5.5fc component stars. We discuss the varia- tions in the context of the colliding-wind model using broad-band spectrometry from the RXTE, Swift, and NICER observatories obtained over 20 years and nearly 1000 observations through 3 consecutive 7.94-year orbits including 3 periastron passages. The X-ray luminosity varies as expected with the inverse of the stellar separation over most of the orbit: departures near periastron are produced when cooling shifts to excess optical emission in C III λ5696 in particular. We use X-ray absorption to esti- mate mass-loss rates for both stars and to constrain the system morphology. The absorption maximum coincides closely with inferior conjunction of the WC star and provides evidence of the ion-reflection mechanism that underlies the formation of collisionless shocks governed by magnetic fields probably generated by the Weibel instability. Comparisons with K-band emission and He I λ10830 absorption show that both are correlated after periastron with the asymmetric X-ray absorption. Dust appears within a few days of periastron suggesting formation within shocked gas near the stagnation point. X-ray flares seen in η Car have not occurred in WR 140, suggesting the absence of large-scale wind inhomogeneities. Relatively constant soft emission revealed during the X-ray minimum is probably not from recombining plasma entrained in outflowing shocked gas.

Keywords: plasmas — shock waves — stars: Wolf-Rayet — stars: individual (WR 140) — X-rays: stars arXiv:2109.10350v1 [astro-ph.HE] 20 Sep 2021 1. WR 140: A COLLIDING-WIND BINARY LABORATORY Massive stars, though extremely rare in the local Universe, are important cosmic engines which generate those chemicals critical for planetary habitability. There are many open questions about how these unusual stars evolve and die, mostly having to do with mass and angular momentum loss during their short lives (e.g. Meynet & Maeder 2003), and how such losses depends on the (evolving) underlying stellar parameters of effective gravity, temperature and chemical abundance, and on initial rotation rate and other complicating factors. Additionally, the majority of massive stars either have binary companions or were part of a multiple system at some point in their history (Sana et al. 2012), so that binary-driven mass exchange and loss can play an important evolutionary role (Vanbeveren et al. 1979; Eldridge et al. 2008), and may play a key role in forming massive gamma-ray burst progenitors by enhancing the stellar rotation rate (Langer et al. 2010). Massive binaries are the progenitors of the gravitational wave mergers seen