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

Turbulent magneto-genesis in a collisionless plasma

F. Pucci,1, 2 M. Viviani,3 F. Valentini,3 G. Lapenta,1 W. H. Matthaeus,4 and S. Servidio3

1Centre for Mathematical Plasma Astrophysics, Department of Mathematics, KU Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium 2Istituto per la Scienza e Teconologia dei Plasmi, Consiglio Nazionale delle Ricerche (ISTP-CNR), Via Amendola 122/D, 70126 Bari, Italy 3Dipartimento di Fisica, Universit`adella Calabria, I-87036 Cosenza, Italy 4Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA

ABSTRACT We investigate an efficient mechanism for generating magnetic fields in turbulent, collisionless plas- mas. By using fully kinetic, particle-in-cell simulations of an initially non-magnetized plasma, we inspect the genesis of magnetization, in a nonlinear regime. The complex motion is initiated via a Taylor-Green vortex, and the plasma locally develops strong electron temperature anisotropy, due to the strain tensor of the turbulent flow. Subsequently, in a domino effect, the anisotropy triggers a Weibel instability, localized in space. In such active wave-particle interaction regions, the magnetic field seed grows exponentially and spreads to larger scales due to the interaction with the underlying stirring motion. Such a self-feeding process might explain magneto-genesis in a variety of astrophysical plasmas, wherever turbulence is present.

Keywords: Turbulence, dynamo, plasma

1. INTRODUCTION The generation of seed magnetic fields is relevant not Magnetic fields permeate the Universe, but their ori- only in astrophysics but also in laboratory settings such gin still represents an open question Kulsrud & Zweibel as laser-plasma experiments Schoeffler et al.(2016). A (2008), despite years of studies (see Widrow et al.(2012) proposed explanation invokes effects beyond standard and references therein). The magnetic fields of most magnetohydrodynamics, such as the Biermann battery stars and galaxies are believed to be sustained and am- Biermann(1950), which relies on the different inertial plified by hydromagnetic dynamo action Parker(1970), response of electrons and ions to a pressure gradient. an essential element of which is turbulence Kraichnan & In the case of weakly-collisional media, one of the most Nagarajan(1967); Pouquet et al.(1976); Mininni et al. efficient candidates for the so-called magneto-genesis is (2003); Brandenburg(2018). In turn, turbulent motion the kinetic Weibel instability Weibel(1959). This pro- is stirred by the evolution of baryonic and dark mat- cess, based on the instability of strongly anisotropic par- ter, with stellar collisions being suggested as a possible ticle distribution functions, has been verified in a vari- mechanism for increasing local shearing motions Colgate ety of simulations and plasma experiments Huntington et al.(2001). For extragalactic plasmas, also dynamos et al.(2015); Schoeffler & Silva(2020). In the collision- in non-collisional and weakly collisional plasmas have less case, it dominates the dynamics at small scales at been investigated Rincon et al.(2016); St-Onge & Kunz ` ∼ de, where de is the electron inertial length. Although such kinetic instability is a very powerful magnetic field

arXiv:2109.07432v1 [physics.plasm-ph] 15 Sep 2021 (2018); Rincon(2019); Pusztai et al.(2020); St-Onge et al.(2020). In both these scenarios, an initial seed source, it requires an existing highly non-Maxwellian of magnetic field is assumed. It is therefore essential plasma, namely an ad-hoc, unstable velocity distribu- to investigate the origin of such seeds from which mag- tion. In the Weibel scenario, such anisotropic distribu- netic fields can emerge and grow as we can observe them tion is a given ingredient. today on large and small scales. More recently, in a sequence of inspiring works, it has been suggested that collisionless plasmas might de- velop large temperature anisotropy, and hence non- [email protected] Maxwellian distributions, not only via an existing mag- arXiv:2109.07423v1 [hep-ph] 15 Sep 2021 n tutr omto 2.Isaudnecudas edet the to due g [3]. be condensate a also gluon to a [1 could by extends abundance potential e.g., Its wavelength QCD generated, effective Broglie potential [2]. the formation de of structure whose minimum and condensate the mec around bosonic misalignment field light the axion from the arises of population cosmological whose axion where eedn nteszso h uaaculn n asso h sca s the of analytic masses V. approximate and Section coupling IV, Yukawa in the and conclude freez of we the III sizes th f with the Sections general constitutes on connection the In depending may in describe fermions it we constraints. thermal result, II, (BBN) to Section a coupling In (qu field As fermion boson model. thermal the state. a of coherent to parameters the the couples in of boson term evolution source scalar and a mass that temperature-dependent generates Suppose 5]. [4, coupling obtains ymtybekn rmahge iesoa operator dimensional higher a from breaking symmetry eciigtehmgnoschrn silto olw h evolutio the follows oscillation coherent homogeneous the describing saS ig-emo-emo prtr ntetemlbackgrou thermal the In operator. Higgs-fermion-fermion SM a is itiuino h fermion the of distribution where ∗ lcrncades [email protected] address: Electronic nti ril,w ics oe ehns o h eeaino a of generation the for mechanism novel a discuss we article, this In easm ih boson light a assume We n fatatv akmte addt salgtbsni coheren a in boson light a is candidate matter dark attractive of One f H ol eaquark a be could y steHbl aaee,adtetemleetecddin encoded effect thermal the and parameter, Hubble the is φ : ooi akmte nachrn tt rvnb hra fer thermal by driven state coherent a in matter dark Bosonic etrso hsmcaimadaayehwi ok nvariou numbers: in PACS obtained. and works be coupling it can the how solutions of analyze analytic range and where wide mechanism limits thro a this in of generated candidate features be matter can dark oscillation good coherent its background, fteei clrbsnfil neatn oiatywt a with dominantly interacting field boson scalar a is there If I HRA FETI H VLTO FABSNFIELD BOSON A OF EVOLUTION THE IN EFFECT THERMAL II. q charged-lepton , φ f olwn h aclto nafiietmeaue(4 ]adreferen and 5] ([4, temperature finite a in calculation the Following . ihmass with 1 oe nttt o dacdSuy eu 25,Korea 02455, Seoul Study, Advanced for Institute Korea m φ ¨ T 2 ( t 3 + ) = m g 24 f H φ l y rneutrino or , φ 2 htculst tnadMdl(M emo hog iyYuka tiny a through fermion (SM) Model Standard a to couples that L φ ˙ T .INTRODUCTION I. ( φff t 2 ( + ) , ugJnChun Jin Eung = and m y φ φ 2 φ hL + φ ( ν f L ¯ hsitrcini upsdt rs fe h electroweak the after arise to supposed is interaction This . φff L m SM f T 2 R i lsia qaino oingvrigtecosmological the governing motion of equation classical ) / T + hr sactffsaeand scale cutoff a is Λ where Λ φ ( = 1, f t aue ftecsooia vlto faclassical a of evolution cosmological the of eatures ¯ = ) R ∗ y do h emos h lsia oo field boson classical the fermions, the of nd g h hra ffc n eoe a becomes and effect thermal the ugh equation: n f a oo edadtetemlfrin Finally, fermion. thermal the and field boson lar φ L ltosaedrvdcnieigvroslimits various considering derived are olutions bevddr atrdniydpnigon depending density matter dark observed e oeetbsncsaei hra background thermal a in state bosonic coherent iainetpouto yacraneffective certain a by production misalignment φ -nmcaimadBgBn Nucleosynthesis Bang Big and mechanism e-in r rlpo) hnistemldistribution thermal its Then lepton). or ark ) ∂φ iutoscnieigasymptotic considering situations s ∂  aimwihdie h oeetoscillation coherent the drives which hanism lci cl ffcigtedr atrhalo matter dark the affecting scale alactic g ur ralpo ntethermal the in lepton a or quark ass edsrb h general the describe We masses. m hL f .Aohritrsigeapei ultra- a is example interesting Another ]. m tt.Acasceapei h QCD the is example classic A state. t T 2 24 φff f and T i 2 T  hL φff i T rssfo h thermal the from arises L mions SM e hri) one therein), ces = f ¯ L Φ f R + φ h.c. (1) (2) (3) wa ( t ) Prepared for submission to JCAP

Structure Formation in the Effective Field Theory of Holographic Dark Energy

Alexander Ganz, Chunshan Lin

aFaculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, 30-348 Krakow, Poland arXiv:2109.07420v1 [gr-qc] 15 Sep 2021 E-mail: [email protected], [email protected]

Abstract. We investigate the structure formation in the effective field theory of the holo- graphic dark energy. The equation of motion for the energy contrast δm of the cold dark matter is the same as the one in the general relativity up to the leading order in the small scale limit k  aH, provided the equation of state is Quintessence-like. Our effective field theory breaks down while the equation of state becomes phantom-like. We propose a solution to this problem by eliminating the scalar graviton. SciPost Physics Submission

An introduction to axions and their detection

Igor G. Irastorza

Center for Astroparticle and High Enery Physics (CAPA), Universidad de Zaragoza, 50009 Zaragoza, Spain. [email protected]

September 16, 2021

Abstract In these notes I try to introduce the reader to the topic of axions: their theoretical motivation and expected phenomenology, their role in astrophysics and as dark matter candidate, and the experimental techniques to detect them. Special emphasis is made in this last point, for which a relatively updated review of worldwide efforts and future prospects is made. The material is intended as an introduction to the topic, and it was prepared as lecture notes for Les Houches summer school 2021. Abundant references are included to direct the reader to deeper insight on the different aspects of axion physics.

Contents

1 Introduction 2

2 Introduction to axion theory and phenomenology 3 2.1 The strong CP problem 3 2.2 The Peccei-Quinn mechanism 4 2.3 Axion mixing and mass 5 2.4 Axion-photon coupling 5 2.5 Axion models 6 2.6 Axion-nucleon couplings 7 2.7 Axion-electron coupling 8 2.8 Other phenomenology 8 2.9 Axion-like particles 9 arXiv:2109.07376v1 [hep-ph] 15 Sep 2021 3 Axions in cosmology 9 3.1 Axions as cold dark matter 11 3.1.1 The axion potential and the vacuum realignment (VR) mechanism 11 3.2 Decay of topological defects (TD) 13 3.3 The domain wall problem 14 3.4 Preferred ma values for axion dark matter 15 3.5 ALP dark matter 15 3.6 Other cosmological phenomenology and constraints on axion/ALP properties 16

4 Axions in astrophysics 17 4.1 ALPs and axions in stellar evolution 18 4.1.1 ALP-photon coupling 18

1 Performance evaluation of GAGG(Ce)/LFS scintillator + MPPC array readout with ASIC

Daiki Nobashia,∗, Kazutaka Yamaokaa, Hiroyasu Tajimaa, Kazuya Itob

aInstitute for Space-Earth Environmental Research (ISEE), Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan bElectronic Information Technology Office, Technical center, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan

Abstract We constructed a gamma-ray detector by combining two types of scintil- lator array detectors with an MPPC array and evaluated the spectral per- formance by reading out the signals from the MPPC with a low-power inte- grated circuit (ASIC) manufactured by IDEAS in Norway. One of the two types of scintillators is a GAGG(Ce) (Ce-doped Gd3Al2Ga3O12) scintillator, and the other is an LFS scintillator. The scintillator array is 2.5 cm × 2.5 cm in size and is coated with BaSO4-based white paint for GAGG(Ce) and an enhanced specular reflector (ESR) for LFS except for the side optically coupled to the MPPC. The spectra derived from the array are affected by the MPPC photon saturations and light leakage from the adjacent pixels, and we carefully corrected for both effects in our data analysis. The energy resolution of 662 keV at 20 ◦C is 6.10±0.04% for the GAGG(Ce) scintillator array and 8.57±0.15% for the LFS scintillator array, this is equivalent to the typical energy resolution found in the references. The energy resolution depends on the temperature: the energy resolution improves as the temper- ature decreases. We found that the contribution of thermal noise from the MPPCs to the energy resolution is negligible within the range of –20 to 40 ◦C, and the energy resolution is mainly determined by the light yield of the crystals.

arXiv:2109.07362v1 [physics.ins-det] 15 Sep 2021 Keywords: MPPC, GAGG(Ce) scintillator array, LFS scintillator array, ASIC

∗Corresponding author Email address: [email protected] (Daiki Nobashi)

Preprint submitted to Journal of LATEX Templates September 16, 2021 Cosmology of an Axion-Like Majoron

A.J. Cuesta,(1) M.E. Gómez,(2) J.I. Illana,(3) M. Masip(3)

(1)Departamento de Física Universidad de Córdoba, E-14071 Córdoba, Spain

(2)Departamento de Ciencias Integradas Universidad de Huelva, E-21071 Huelva, Spain

(3)CAFPE and Departamento de Física Teórica y del Cosmos Universidad de Granada, E-18071 Granada, Spain

[email protected], [email protected], jillana,[email protected] September 16, 2021

Abstract

We propose a singlet majoron model that defines an inverse seesaw mechanism

in the ν sector. The majoron φ has a mass mφ ≈ 0.5 eV and a coupling to the τ lepton similar to the one to neutrinos. In the early universe it is initially in thermal equilibrium, then it decouples at T ≈ 500 GeV and contributes with just

∆Neff = 0.026 during BBN. At T = 26 keV (final stages of BBN) a primordial magnetic field induces resonant γ ↔ φ oscillations that transfer 6% of the photon

energy into majorons, implying ∆Neff = 0.55 and a 4% increase in the baryon

arXiv:2109.07336v1 [hep-ph] 15 Sep 2021 to photon ratio. At T ≈ mφ the majoron enters in thermal contact with the heaviest neutrino and it finally decays into νν¯ pairs near recombination, setting

∆Neff = 0.85. This boost in the expansion rate at later times solves the Hubble tension, while the neutrino–majoron interactions suppress the ν free streaming and make the model consistent with large scale structure observations. Its lifetime and the fact that it decays into neutrinos instead of photons lets this axion-like majoron avoid the strong bounds that affect other axion-like particles of similar mass and coupling to photons. Eur. Phys. J. C manuscript No. (will be inserted by the editor)

A measurement of the mean electronic excitation energy of liquid xenon

Laura Baudis, Patricia Sanchez-Lucas, Kevin Thiemea Department of Physics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland

Received: date / Accepted: date

Abstract Detectors using liquid xenon as target are widely dark photons [7–9] as well as for WIMP (weakly interact- deployed in rare event searches. Conclusions on the interact- ing massive particle) dark matter in the GeV to TeV mass ing particle rely on a precise reconstruction of the deposited range. In particular, TPCs operated in dual-phase mode with energy which requires calibrations of the energy scale of the a gaseous xenon (GXe) layer at the top are among the lead- detector by means of radioactive sources. However, a micro- ing technologies in the past, present and near-future hunt for scopic calibration, i.e. the translation from the number of ex- WIMPs [10–16]. citation quanta into deposited energy, also necessitates good A particle that deposits energy in LXe yields scintilla- knowledge of the energy required to produce single scintil- tion photons in the vacuum-ultraviolet (VUV) range with lation photons or ionisation electrons in liquid xenon. The a peak centred at 178nm wavelength, ionisation electrons sum of these excitation quanta is directly proportional to the and heat via atomic motion, where only the former two deposited energy in the target. The proportionality constant processes are detectable with dual-phase TPCs. Photosen- is the mean excitation energy and is commonly known as sors detect the prompt scintillation light (S1) which is com- W-value. Here we present a measurement of the W-value posed of direct scintillation of the xenon atoms and light with electronic recoil interactions in a small dual-phase from the recombination of xenon ions with electrons. In both xenon time projection chamber with a hybrid (photomulti- mechanisms, scintillation light is produced by de-excitation plier tube and silicon photomultipliers) photosensor config- of xenon dimers that form from xenon excitons and xenon uration. Our result is based on calibrations at O(1−10keV) atoms. Electrons that do not recombine are vertically drifted with internal 37Ar and 83mKr sources and single electron and extracted to the gas phase, by means of an electric drift +0.2 events. We obtain a value of W = 11.5 −0.3 (syst.)eV, with and extraction field, where they collide with xenon atoms. In negligible statistical uncertainty, which is lower than previ- this process, a secondary scintillation signal (S2) is produced ously measured at these energies. If further confirmed, our by electroluminescence which is also detected by photosen- result will be relevant for modelling the absolute response sors. Due to energy conservation, the size of the signals from of liquid xenon detectors to particle interactions. scintillation and ionisation are anti-correlated [17]. The re- combination fraction and thus the energy distribution be- tween these is drift field dependent. We collectively refer 1 Introduction to scintillation and ionisation as excitation in this work. The number of scintillation photons nγ and ionisation Liquid xenon (LXe) is widely used as both sensitive target electrons ne− in an electronic recoil interaction is linearly arXiv:2109.07151v1 [physics.ins-det] 15 Sep 2021 and radiation source in time projection chambers (TPCs) related to the deposited energy E by the constant W [18]: for present and future rare event searches. These include E = (nγ + ne− )W . (1) e.g. searches for neutrinoless double beta decay [1–3], the extremely rare decay of 124Xe via double electron cap- The W-value can be regarded as the average energy needed ture [4,5], the measurement of low-energy solar neutri- to produce a single free quantum in LXe and the above ex- nos [6], searches for solar axions, axion-like particles and pression as its defining equation. It determines the underly- ing recombination-independent energy scale in a LXe detec- ae-mail: [email protected] tor that detects both scintillation light and ionisation charge. Spatiotemporal Characterization of VIIRS Night Light

Christopher Small Lamont Doherty Earth Observatory Columbia University Palisades, NY USA

[email protected]

Abstract

The Visible Infrared Imaging Radiometer Suite (VIIRS) Day Night Band (DNB) on board the Suomi NPP satellite now provides almost a decade of daily observations of night light. The temporal frequency of sampling, without the degree of temporal averaging of annual composites, makes it necessary to consider the distinction between apparent temporal changes of night light related to the imaging process and actual changes in the underlying sources of the night light being imaged. The most common approach to night light change detection involves direct attribution of observed changes to the phenomenon of interest. Implicit in this approach is the assumption that other forms of actual and apparent change in the light source are negligible or non-existent. An alternative approach is to characterize the spatiotemporal variability prior to deductive attribution of causation so that the attribution can be made in the context of the full range of spatial and temporal variation. The primary objective of this study is to characterize night light variability over a range of spatial and temporal scales to provide a context for interpretation of night light changes observed on both subannual and interannual time scales. This analysis is based on a combination of temporal moments, spatial correlation and Empirical Orthogonal Function (EOF) analysis. A key result of this study is the pervasive heteroskedasticity of VIIRS monthly mean night light. Specifically, the monotonic decrease of temporal variability with increasing mean brightness. Anthropogenic night light is remarkably stable on subannual time scales. The overall variance partition derived from the eigenvalues of the spatiotemporal covariance matrix are 88%, 2% and 2% for spatial, seasonal and interannual variance (respectively) in the most diverse geographic region on Earth (Eurasia). Heteroskedasticity is a pervasive characteristic of monthly VIIRS composites and is present in all areas for all months of the year, suggesting that much, if not most, of observed month-to-month variability may be related to luminance of otherwise stable sources subjected to multiple aspects of the imaging process varying in time. Given the skewed distribution of all night light arising from radial peripheral dimming of bright sources subject to atmospheric scattering, even aggregate metrics using thresholds must be interpreted in light of the fact that much larger numbers of more variable low luminance pixels may statistically overwhelm smaller numbers of stable higher luminance pixels and cause apparent changes related to the imaging process to be interpreted as actual changes in the light sources.

Introduction

The Visible Infrared Imaging Radiometer Suite (VIIRS) Day Night Band (DNB) on board the Suomi NPP satellite now provides almost a decade of daily observations of night light. In comparison to the Defense Meteorological Satellite Program (DMSP) Operational Line Scanner (OLS) night light imagery, VIIRS DNB provides greater dynamic range (14 vs 6 bit), higher

1 Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX63

Local Environments of Low-Redshift Supernovae

Serena A. Cronin,1, 2, 3 Dyas Utomo,4, 1, 5, ∗ Adam K. Leroy,1, 5 Erica A. Behrens,6, 1 Jer´ emy´ Chastenet,7 Tyler Holland-Ashford,1, 5 Eric W. Koch,8 Laura A. Lopez,1, 5 Karin M. Sandstrom,9 and Thomas G. Williams10

1Department of Astronomy, The Ohio State University, 4055 McPherson Laboratory, 140 West 18th Avenue, Columbus, OH 43210, USA 2Department of Astronomy, University of Maryland, College Park, MD 20742, USA 3National Radio Astronomy Observatory, 1003 Lopezville Road, Socorro, NM 87801, USA 4National Radio Astronomy Observatory, 520 Edgemont Rd., Charlottesville, VA 22903, USA 5Center for Cosmology and AstroParticle Physics, 191 West Woodruff Avenue, Columbus, OH 43210, USA 6Department of Astronomy, University of Virginia, 530 McCormick Rd., Charlottesville, VA, 22904, USA 7Sterrenkundig Observatorium, Ghent University, Kri-jgslaan 281 − S9, 9000 Gent, Belgium 8Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., 02138 Cambridge, MA, USA 9Center for Astrophysics and Space Sciences, Department of Physics, University of California, San Diego 9500 Gilman Drive, La Jolla, CA 92093, USA 10Max Planck Institut f¨urAstronomie, K¨onigstuhl17, 69117 Heidelberg, Germany

(Received June 21, 2021; Revised August 29, 2021; Accepted September 14, 2021) Submitted to ApJ

ABSTRACT We characterize the local (2-kpc sized) environments of Type Ia, II, and Ib/c supernovae (SNe) that have recently occurred in nearby (d . 50 Mpc) galaxies. Using ultraviolet (UV, from GALEX) and infrared (IR, from WISE) maps of 359 galaxies and a sample of 472 SNe, we measure the star formation rate surface density (ΣSFR) and stellar mass surface density (Σ?) in a 2-kpc beam centered on each SN site. We show that core-collapse SNe are preferentially located along the resolved galactic star-forming main sequence, whereas Type Ia SNe are extended to lower values of ΣSFR at fixed Σ?, indicative of locations inside quiescent galaxies or quiescent regions of galaxies. We also test how well the radial distribution of each SN type matches the radial distributions of UV and IR light in each host galaxy. We find that, to first order, the distributions of all types of SNe mirror that of both near-IR light (3.4 and 4.5 µm, tracing the stellar mass distribution) and mid-IR light (12 and 22 µm, tracing emission from hot, small grains), and also resemble our best-estimate ΣSFR. All types of SNe appear more radially concentrated than the UV emission of their host galaxies. In more detail, the distributions of Type II SNe show small statistical differences from that of near-IR light. We attribute this overall structural uniformity to the fact that within any individual galaxy, ΣSFR and Σ? track one another well, with variations in ΣSFR/Σ? most visible when comparing between galaxies.

Keywords: star formation, galaxies, supernovae arXiv:2109.07453v1 [astro-ph.GA] 15 Sep 2021 1. INTRODUCTION Kochanek et al. 2017), the Zwicky Transient Facility Surveys like the Lick Observatory Supernova Search (ZTF; Bellm et al. 2019), and the Asteroid Terrestrial- (LOSS; Filippenko et al. 2001), the All-Sky Automatic impact Last Alert System (ATLAS; Tonry et al. 2018) Survey for Supernovae (ASAS-SN; Shappee et al. 2014; have led to a surge in supernova (SN) detections in re- cent years. The identification of many new supernovae (SNe) has improved our knowledge of their sub-types Corresponding author: Serena A. Cronin and locations. Nonetheless, fundamental questions re- [email protected] main about the progenitors of each SN type and the ∗ Jansky Fellow of the National Radio Astronomy Observatory environments in which they explode. Astronomy & Astrophysics manuscript no. psmap ©ESO 2021 September 16, 2021

A new method to perform data/model comparison in Fermi-LAT analysis P. Bruel

Laboratoire Leprince-Ringuet, CNRS, Ecole polytechnique, Institut Polytechnique de Paris, F-91128 Palaiseau, France e-mail: [email protected]

Received 15 June 2021; accepted 15 September 2021

ABSTRACT

Context. The analysis of Fermi Large Area Telescope (LAT) gamma-ray data in a given Region Of Interest (RoI) usually consists of performing a binned log-likelihood fit in order to determine the sky model that, after convolution with the instrument response, best accounts for the distribution of observed counts. Aims. While tools are available to perform such a fit, it is not easy to check the goodness-of-fit. The difficulty of the assessment of the data/model agreement is twofold. First of all, the observed and predicted counts are binned in three dimensions (two spatial dimensions and one energy dimension) and comparing two 3D maps is not straightforward. Secondly, gamma-ray source spectra generally decrease with energy as the inverse of the energy square. As a consequence the number of counts above several GeV generally falls into the Poisson regime, which precludes performing a simple χ2 test. Methods. We propose a method that overcomes these two obstacles by producing and comparing spatially integrated count spectra for data and model at each pixel of the analysed RoI. The comparison is performed following a log-likelihood approach that extends the χ2 test to histograms with low statistics. This method can take into account likelihood weights that are used to account for systematic uncertainties. Results. We optimize the new method so that it provides a fast and reliable tool to assess the goodness-of-fit of Fermi-LAT data and we use it to check the latest gamma-ray source catalog on 10 years of data. Key words. gamma rays: general — methods: data analysis — methods: statistical

1. Introduction source spectral parameters that give the best agreement between the observed and predicted 3D count maps. Since its launch in June 2008, the Fermi Large Area Tele- After the fit has converged, it is important to check that the scope (LAT, Atwood et al. 2009) has been continuously observ- obtained agreement is satisfactory but the user faces the difficulty ing the gamma-ray sky in the energy range between 30 MeV and 1 of comparing two 3D maps. The simplest way would be to de- 2 TeV. The Fermi-LAT data are public and have been widely fine an energy band and produce the corresponding 2D residual used by the gamma-ray community. The Fermi-LAT latest gen- (data − model) count map. But the energy band definition and eral catalog, 4FGL-DR2 (Abdollahi et al. 2020; Ballet et al. the analysis of the residual map is not straightforward because 2020), reports almost 5800 sources, among which are hundreds of the energy dependence of the Point Spread Function (PSF) of pulsars, tens of supernovae remnants and pulsar wind nebulae, of the instrument: its 68% containment angle varies from 5◦ at and thousands of blazars. 100 MeV to 0.1◦ at 30 GeV (Atwood et al. 2009, 2013) 4. As a The usual way to analyse LAT data in a given Region of consequence the binning of the residual map should depend on Interest (RoI) is to bin the data in a three dimensional map the energy range definition. Furthermore, since the spectral char- (two spatial dimensions and one energy dimension 2) and then acteristics of a potential data/model mismatch are unknown, it is search for the sky model that best predicts the number of gamma not obvious to choose a priori an energy range, leading often to rays observed in this RoI. The sky model is a list of gamma- look at several energy bands. ray sources, whose position, spatial nature (point-like or ex- The Fermitools provide a method to quantify the tended), as well as their energy spectrum are provided. The data/model agreement which is frequently used. It consists of arXiv:2109.07443v1 [astro-ph.HE] 15 Sep 2021 3 Fermi-LAT analysis package, Fermitools , is used to convolve computing a Test Statistic (TS) map: at each pixel, the presence the emission of each source with the Instrument Response Func- of an additional source is tested by computing twice the differ- tions (IRFs, Ackermann et al. 2012) in order to predict the num- ence in log-likelihood obtained with and without the source. A ber of observed counts. As a consequence, one of the main steps TS = 25 corresponds to ∼ 4σ significance (Mattox et al. 1996). of the analysis is generally a log-likelihood fit that finds the The drawbacks of this method are twofold: it is computationally intensive and, above all, it is not sensitive to negative deviations 1 https://fermi.gsfc.nasa.gov/ssc (data

Article number, page 1 of 13 Concepts and status of Chinese space gravitational wave detection projects

Yungui Gong ,1, ∗ Jun Luo,2, 3 and Bin Wang 4, 5, † 1School of Physics, Huazhong University of Science and Technology, Wuhan, China 2MOE Key Laboratory of TianQin Mission, TianQin Research Center for Gravitational Physics & School of Physics and Astronomy, Frontiers Science Center for TianQin, CNSA Research Center for Gravitational Waves, Sun Yat-sen University, Zhuhai, China 3Center for Gravitational Experiments, School of Physics, MOE Key Laboratory of Fundamental Physical Quantities Measurement & Hubei Key Laboratory of Gravitation and Quantum Physics, PGMF, Huazhong University of Science and Technology, Wuhan, China 4Shanghai Frontier Science Research Center for Gravitational Wave Detection, School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai, China 5Center for Gravitation and Cosmology, Yangzhou University, Yangzhou, China. Gravitational wave (GW) detection in space probes GW spectrum that is inaccessible from the Earth. In addition to LISA project led by European Space Agency, and the DECIGO detector proposed by the Japan Aerospace Exploration Agency, two Chinese space-based GW observatories– TianQin and Taiji–are planned to be launched in the 2030s. TianQin has a unique concept in its design with a geocentric orbit. Taiji’s design is similar to LISA, but is more ambitious with longer arm distance. Both facilities are complementary to LISA, considering that TianQin is sensitive to higher frequencies and Taiji probes similar frequencies but with higher sensitivity. In this Perspective we explain the concepts for both facilities and introduce the development milestones of TianQin and Taiji projects in testing extraordinary technologies to pave the way for future space-based GW detections. Considering that LISA, TianQin and Taiji have similar scientific goals, all are scheduled to be launched around the 2030s and will operate concurrently, we discuss possible collaborations among them to improve GW source localization and characterization.

I. INTRODUCTION vanced LIGO14,15 and Virgo16 detectors, the Kamioka Gravitational Wave Detector (KAGRA)17,18 has joined After the first direct detection of gravitational wave the network of ground-based detectors. In the near fu- (GW) in September 20151,2, there have been 50 reported ture, LIGO India will also join the network. GW detections by the ground-based Laser Interferome- Ground-based detectors are not sensitive to GWs be- ter Gravitational-Wave Observatory (LIGO) and Virgo low 1Hz because of terrestrial gravity gradient noise. A observatory3,4. These GW events were identified as space-based detector is free from such noise and can be merges of binary black holes (BBH)1,2,5–10, binary neu- made very large, thereby expanding the frequency range −4 tron stars11,12, BH-neutron star mergers or the coales- downwards to 10 Hz, where exciting GW sources are cence of a BH with a mystery compact object13. waiting to be explored. Space-based observations can For ground-based GW detectors, their high frequency probe dynamic coalescing systems, characterize how they range (10-1,000Hz) implies that even stellar-mass binary grow, pair up and merge, and estimate the space den- mergers emit GWs for only short durations. Multiple sity and environments of massive BHs (MBHs). Cover- interferometers are needed to confidently detect GWs ing different GW frequency ranges, space-based GW ob- and eliminate false alarm events. Using the time delay servations will provide a unique tool for testing gravity. between multiple detectors in different locations allows Some highly relativistic events, such as BBH coalescences 5 the sky location of the GW source to be pinpointed us- with masses below 10 solar masses (M ), last a year or ing triangulation. The large distance between the ad- longer. This long duration allows a single space-based vanced Virgo and the advanced LIGO detectors helps detector to provide directional information as it orbits in them achieve different orientations due to the Earth’s the sky during the observation. The equilateral triangle arXiv:2109.07442v1 [astro-ph.IM] 15 Sep 2021 curvature and because their interferometer arms are not design of space-based GW interferometers can provide in- parallel to each other it is expected that the polariza- dependent information on the GW polarization. Further- tion of GWs can be extracted from LIGO-Virgo net- more, the improved sensitivities of space-based detectors, work data7. A global network of ground-based GW afforded due to the increased interferometer arm length, detectors will play a key role in exploring the popula- will enable improved searches for the cosmological GW tion of BHs and neutron stars with GWs in the fre- background. Note that if the arm length is comparable quency band around 100 Hz. In addition to the ad- or even longer than the wavelength of GWs, then the sensitivity deteriorates. Proposed projects for space-based GW observato- ries include the Laser Interferometer Space Antenna ∗ [email protected] (LISA)19,20 led by European Space Agency (ESA), the † wang [email protected] Deci-hertz Interferometer Gravitational Wave Observa- arXiv:2109.07427v1 [astro-ph.GA] 15 Sep 2021 Cm 04 n eeecsteen.Mr eety loex also al. recently, C et that found More Berné observations, therein). emission mid-infrared references ploiting o and astrophysical various 2014, in then (Cami 2010), al. et Sellgren 2010; i-nrrdeiso bevtosrvae h presence the revealed C fullerene observations emission Mid-infrared Introduction 1. swl.A otefleeecto C cation fullerene the to As well. as 06,;Wle ta.21,21,21) fe eea test several After 2017). C of 201 2016, identification al. 2015, the et al. (Campbell et Walker wavelengths 2016a,b; di near-infrared of ident at carrier was observed the it as 1997), recently 1994, fied Ehrenfreund & (Foing anticipation Crie ta.21,adrfrne hri) demonstrat C therein), of references presence and 2019, al. et (Cordiner bevditrtla ulrnsC fullerenes interstellar observed C n eeecsteen.A neapeo plcto,the application, of example an As therein). C al. references et (Brieva and environments astrophysical various probe to 95.Brée l 21)eautdteC the evaluated (2017) al. Tie et & Berné (Bakes 1995). recombination it electron since and density photoionization electron to and (ISRF) field radiation stellar ueIMad rmteosre i-nrrdeiso bands emission mid-infrared observed C the of from and, ISM fuse tained etme 6 2021 16, September Astronomy ⋆ 60 − 60 nI msinsetao h trfrigrgo C348. IC region star-forming the of spectra emission IR in , :C e-mail: h bqiyo C of ubiquity The h C The 60 C Conclusions. Results. Methods. Aims. Context. e words. Key non-detecti the observations. despite and data notable laboratory be may clouds interstellar h irncbnso h ain h ouaino eta m neutral of population The cation. the of bands vibronic di the main the of characteristics published from Received ad oa o orl nabgosI msinmeasurements emission IR ambiguous on rely to not as so bands ac n hresaeo hi components. their of state charge and dance ad fC of bands oeueioae nN c ocmlt them. complete to ice Ne in isolated molecule 60 + 2 1 n h esrmn fDB trbtdt C to attributed DIBs of measurement the and 60 nvriyo il,CR,UR82 hA aoaor eP de Laboratoire - PhLAM - 8523 UMR Germany CNRS, Lille, Jena, f of D-07743 University Institute 3, Planck Helmholtzweg Physics, Max State the Solid of Group Astrophysics Laboratory ai nteIMgvsifraino h oa inter- local the on information gives ISM the in ratio :C [email protected] 60 + 60 epooea vlaino h C the of evaluation an propose We 60 & ratio. npaeayadrflcinnblefis Cm tal. et (Cami first nebulae reflection and planetary in 60 + esaedw h siltrsrntsfrasrto bands absorption for strengths oscillator the down scale We nih notecniin htdietepyisadchemi and physics the drive that conditions the into Insight < speeti h di the in present is h ai sdtrie yaayigacia pcr n li and spectra archival analyzing by determined is ratio The Astrophysics 60 60 date ntedi the in r o ruh u yosrain.W eieteoscillato the revise We observations. by out brought not are srceity–IM bnacs–IM ie n ad ISM – bands and lines ISM: – abundances ISM: – astrochemistry :C ecnld httefato fnurlmlclsi h buc the in molecules neutral of fraction the that conclude We / > 60 60 + 60 + Accepted ff ulrn ae ta trciespecies attractive an it makes fullerene stecriri osdrdt efirm be to considered is carrier the as s S.Lsl,Ilsa-rt (2019) Iglesias-Groth Lastly, ISM. use ai nDfueadTasuetItrtla Clouds Interstellar Translucent and Diffuse in Ratio aucitn.roui2108 no. manuscript alRouillé Gaël < ff 60 ff date s neselrmdu (ISM) medium interstellar use s neselrbns(DIBs) bands interstellar use n C and > 60 + fe w eae of decades two after , 70 60 n h osC ions the and , 60 :C 1 ⋆ :C ff 60 + eg .Krasnokutski A. Serge , s neselrbnsatiue oC to attributed bands interstellar use 60 + ai o h dif- the for ratio ai ndi in ratio 60 + hyob- they , n the ing (2017) relates 60 + ff bjects no h xetdvboi ad.Mr etit ilrequi will certainty More bands. vibronic expected the of on 2016, s n rnlcn neselrcod htepot elec exploits that clouds interstellar translucent and use ABSTRACT lclsi ecie ntrso pe ii eas h rel the because limit upper of terms in described is olecules lens and of 5, s, i- - . rAtooya h rerc cilrUiest ea Ins Jena, University Schiller Friedrich the at Astronomy or isefleee ti hrceie ya pe ii bec limit C upper of an bands by absorption characterized the is it minsterfullerene, C that measurem detailed DIBs published of by us to provided is population fteC the of of range the al. and et (1995) Berné Tielens by & the proposed of Bakes values 2000), by limits prediction (Herbig contr upper the in 0.0075 to with populations, ionized down and strongly indicating 1988) 2017) two last al. Campbell et & (Léger (Maier 1 0.1 from range ues xetda aeegh er40 n 00Å ontappear not do Å, 6000 and 4000 near wavelengths at expected hs fC of di those the probing for suited ter rbtosb nlgu pce Grí-enne ta.2 2019). al. al. et (García-Hernández et Krasnokutski species analogous by tributions ain a earmd,aohriseaie.Ide,teIR the Indeed, C arises. to issue attributed another bands remedy, obser- emission a emission its col be improved of the may though measurements vations of Even from evaluation spectrum. molecule emission the neutral IR in the of a chiefly density that originates umn uncertainty al. et large Berné The by 5.6. of ratio iae rmFg nBks&Tees19) hti,aC a is, that 1995), Tielens & Bakes in 5 Fig. from timated hstenurlseisi osbymr bnatta the than abundant more possibly C di is the in species cation neutral the Thus tyi neselrcod sgie rmdtriigtea the determining from gained is clouds interstellar in stry yiu e aesAoe tMlcls -90 il,Fran Lille, F-59000 Molécules, et Atomes Lasers des hysique 60 eauedt.Ifraino h ainpplto sobta is population cation the on Information data. terature fC of :C aigti iceac,w rps nte determination another propose we discrepancy, this Facing lcrncasrto ad r ycnrs pcfi,as b also specific, contrast by are bands absorption Electronic boueaudne ..,wihpitt ntvr re- very (not a to C point strictive) which [...], abundances absolute tegh eotdfrC for reported strengths r 1 60 + n vi Carpentier Yvain and , 60 :C 60 mntrulrn ouaino di of population kminsterfullerene :C n n nuprlmto prxmtl . o the for 1.3 approximately of limit upper an find and molecules : 60 − 60 60 + 60 + eaieaudne f06:.102 peetyes- (presently 0.62:0.11:0.27 of abundances relative oeaut h C the evaluate to n bopincosscin led eotdfor reported already sections cross absorption and ai ntedi the in ratio 60 ff 60 + s S,i gemn ihamdlpredicting model a with agreement in ISM, use vrC over are.A otepplto fnurlbuck- neutral of population the to As carries. 60 + 60 ai agn ewe . o6. to 0.3 between ranging ratio 60 ff ff htaems rmsn sgauges, as promising most are that n esr h pcrmo the of spectrum the measure and s S,hneatmt oexploit to attempts hence ISM, use s S.Ifraino h cation the on Information ISM. use 60 1 60 , :C 2 a culycmrs con- comprise actually may 60 + ril ubr ae1o 13 of 1 page number, Article ai.Rsligrtoval- ratio Resulting ratio. ff s n translucent and use ff rncabsorption tronic cstertogiven ratio the ects vn vibronic evant eimproved re © iueof titute S 2021 ESO bun- ined adiction ce 60 013; ause :C ents et- 60 + - Solar Physics DOI: 10.1007/•••••-•••-•••-••••-•

Radio Interferometric Observations of the Sun Using Commercial Dish TV Antennas

G. V. S. Gireesh1 · C. Kathiravan1 · Indrajit V. Barve1 · R. Ramesh1

© Springer ••••

Abstract The radio astronomy group in the Indian Institute of Astrophysics (IIA) has been carrying out routine observations of radio emission from the solar corona at low frequencies (≈40-440MHz) at the Gauribidanur observa- tory, about 100km north of Bangalore. Since IIA has been performing regular observations of the solar photosphere and chromosphere using different optical telescopes in its Kodaikanal Solar Observatory (KSO) also∗, the possibilities of obtaining two-dimensional radio images of the solar chromosphere using low-cost instrumentation to supplement the optical observations are being explored. As a part of the exercise, recently the group had developed prototype instrumentation for interferometric observations of radio emission from the solar chromosphere at high frequencies (≈11.2GHz) using two commercial dish TV antennas. The hardware set-up and initial observations are presented.

Keywords: Sun, radio emission. Sun, eclipse. Radio telescopes, interferometry.

1. Introduction

Commercial dish TV antennas are parabolic structures designed to receive radio waves from a communication satellite. The antennas and the associated front end receiver systems have improved with advances in the TV systems. They

G. V. S. Gireesh [email protected] C. Kathiravan [email protected] Indrajit V. Barve [email protected] R. Ramesh [email protected] arXiv:2109.07422v1 [astro-ph.SR] 15 Sep 2021 1 Indian Institute of Astrophysics, 2nd Block, Koramangala, Bangalore 560034 ∗see https://www.iiap.res.in/kodai.htm

SOLA: manuscript.tex; 16 September 2021; 0:43; p. 1 J. Astrophys. Astr. (0000) 000: DOI

An Improved Approach to Orbital Determination and Prediction of Near- Earth Asteroids: Computer Simulation, Modeling and Test Measure- ments

Muhammad Farae1*, Cameron Woo1, Anka Hu1

1Department of Physics, University of Colorado Boulder, Boulder CO 80309 USA

*Corresponding author. E-mail: [email protected]

Abstract. In this article, theory-based analytical methodologies of astrophysics employed in the modern era are suitably operated alongside a test research-grade telescope to image and determine the orbit of a near-earth asteroid from original observations, measurements, and calculations. Subsequently, its intrinsic orbital path has been cal- culated including the chance it would likely impact Earth in the time ahead. More so specifically, this case-study incorporates the most effective, feasible, and novel Gauss’s Method in order to maneuver the orbital plane compo- nents of a planetesimal, further elaborating and extending our probes on a selected near-earth asteroid (namely the 12538-1998 OH) through the observational data acquired over a six week period. Utilizing the CCD (Charge Cou- pled Device) snapshots captured, we simulate and calculate the orbit of our asteroid as outlined in quite detailed explanations. The uncertainties and deviations from the expected values are derived to reach a judgement whether our empirical findings are truly reliable and representative measurements by partaking a statistical analysis based systematic approach. Concluding the study by narrating what could have caused such discrepancy of findings in the first place, if any, measures are put forward that could be undertaken to improve the test-case for future investiga- tions. Following the calculation of orbital elements and their uncertainties using Monte Carlo analysis, simulations were executed with various sample celestial bodies to derive a plausible prediction regarding the fate of Asteroid 1998 OH. Finally, the astrometric and photometric data, after their precise verification, were officially submitted to the Minor Planet Center: an organization hosted by the Center for Astrophysics, Harvard and Smithsonian and funded by NASA, for keeping track of the asteroid’s potential trajectories.

Keywords. Near-Earth Asteroid—Keplerian Orbital Elements—Gauss’s Method—Astrometry—Photometry— Simulation—Two-Body Problem—Least Squares Theory.

1. Introduction As a result of the ubiquity of asteroids residing in our Solar System, tons of space debris enter the Earth’s atmosphere each year, most of which are not massive The Solar System, in its current state, is rife with plan- enough to cause fatal destruction. However, some as- etesimals forming via the collision, spallation, and it- teroids have potential to impact the earth with immense arXiv:2109.07397v1 [astro-ph.EP] 15 Sep 2021 erative fragmentation of larger planetary bodies. The energy and cataclysmic out-turns (Chodas et al. 2001). paramount constituents of these minor celestial bodies With the exponential evolution of imaging technology are namely asteroids and comets, each with a uniquely and the surge in astronomical, and specifically orbital, defined orbit (Milone & Wilson 2014). This study research over the past years, enormous breakthroughs focuses specifically on NEAs (Near-Earth Asteroids) have been witnessed in the orbital determination of which have a perihelion less than 1.3 AU and an aphe- NEAs (Raol & Sinha 1985; Farnocchia et al. 2013). In lion greater than or equal to 0.983 AU from the Sun order to make the future behaviour of asteroids in the (NASA 2007). An Earth-crossing asteroid is one which Solar System more predictable, multiple search cam- has an orbit intersecting the Earth’s around the Sun, and paigns have been established to map the night sky (el- a Mars-crosser is defined analogously. Showk 2017; Yeomans et al. 1997). Our research work

© Indian Academy of Sciences 1 Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX62

Exploring the sky localization and early warning capabilities of third generation gravitational wave detectors in three-detector network configurations.

YUFENG LI,1, 2, 3 Ik Siong Heng,3 Man Leong Chan,4 Chris Messenger,3 and XILONG FAN5

1Key Laboratory for Computational Astrophysics, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China 2School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China 3SUPA, School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom 4Department of Applied Physics, Fukuoka University, Nanakuma 8-19-1, Fukuoka 814-0180, Japan 5School of Physics and Technology, Wuhan University, Wuhan 430072, China

ABSTRACT This work characterises the sky localization and early warning performance of different networks of third generation gravitational wave detectors, consisting of different combinations of detectors with either the Einstein Telescope or Cosmic Explorer configuration in sites in North America, Europe and Australia. Using a Fisher matrix method which includes the effect of earth rotation, we estimate the sky localization uncertainty for 1.4M -1.4M binary neutron star mergers at distances 40Mpc, 200Mpc, 400Mpc, 800Mpc, 1600Mpc, to characterize each network’s performance for binary neutron star observations at a given distance. We also characterize each network’s sky localization capability for an assumed astrophysical population up to redshift ≤ 2. Furthermore, we also study the capabilities for the different networks to localize a binary neutron star merger prior to merger (early warning) and characterise the network performance for sky localization uncertainty between 1 and 30 square degrees. We find that, for example, for binary neutron star mergers at 200Mpc and a network consisting of the Einstein Telescope, Cosmic Explorer and an extra Einstein Telescope-like detector in Australia (2ET1CE), the upper limit of the size of the 90% credible region for the best localized 90% signals is 0.51deg2. For the simulated astrophysical distribution, this upper limit is 183.58deg2. If the Einstein Telescope-like detector in Australia is replaced with a Cosmic Explorer-like detector (1ET2CE), for signals at 200Mpc, the size of the 90% credible region for the best localized 90% signals is 0.36deg2, while the corresponding value for the best localized 90% sources following the delay time distribution is 113.55deg2. We note that the 1ET2CE network can detect 7.2% more of the simulated astrophysical population than the 2ET1CE network. In terms of early warning performance, we find that a network of 2ET1CE and 1ET2CE networks can both provide early warnings of the order of 1 hour prior to merger with sky localization uncertainties of 30 square degrees or less. In some cases, the 2ET1CE network is capable of estimating the sky location with an uncertainty of 5 square degrees or less on timescales of about 1 hour prior to merger. Our study concludes that the 1ET2CE network is a good compromise between binary neutron stars detection rate, sky localization and early warning capabilities. arXiv:2109.07389v1 [astro-ph.IM] 14 Sep 2021 Keywords: gravitational wave sources (677); Gravitational wave detectors (676);

1. INTRODUCTION verse (Abbott et al. 2016). During the first two ob- The first directly detected gravitational wave (GW) servation runs (O1 and O2) of Advanced LIGO, and event GW150914 by Advanced LIGO in September, Advanced Virgo (joined in O2) (Acernese et al. 2015), a 2015 opened a new window on exploring the Uni- total of ten binary black holes (BBHs) mergers and one binary neutron star (BNS) merger were identified (Ab- bott et al. 2019). The follow-up observations of the Corresponding author: Ik Siong Heng electromagnetic (EM) counterparts of the BNS event [email protected] have initiated the GW and EM multi-messenger astron- MNRAS 000,1–12 (2021) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

Cross-correlation of Planck CMB lensing with DESI galaxy groups

Zeyang Sun,1★ Ji Yao,1,2 Fuyu Dong,3 Xiaohu Yang,1,2,4 Le Zhang,5 Pengjie Zhang1,2,4† 1Department of Astronomy, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai, 200240, China 2Shanghai Key Laboratory for Particle Physics and Cosmology, Shanghai, 200240, China 3School of Physics, Korea Institute for Advanced Study (KIAS), 85 Hoegiro, Dongdaemun-gu, Seoul, 02455, Republic of Korea 4Division of Astronomy and Astrophysics, Tsung-Dao Lee Institute, Shanghai Jiao Tong University, Shanghai, 200240, China 5School of Physics and Astronomy, Sun Yat-Sen University, 2 Daxue Road, Tangjia, Zhuhai, 519082, China

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We measure the cross-correlation between galaxy groups constructed from DESI Legacy Imaging Survey DR8 (Yang et al. 2020) and Planck CMB lensing, over overlapping sky area of 16876 deg2. The detections are significant and consistent with the expected signal of the large scale structure of the universe, over group samples of various redshift, mass and richness 푁g and over various scale cuts. The overall S/N is 39 for a conservative sample with 푁g ≥ 5, and increases to 48 for the sample with 푁g ≥ 2. Adopting the Planck 2018 cosmology, we constrain the density bias of groups with 푁g ≥ 5 as 푏g = 1.31 ± 0.10, 2.22 ± 0.10, 3.52 ± 0.20 at 0.1 < 푧 ≤ 0.33, 0.33 < 푧 ≤ 0.67, 0.67 < 푧 ≤ 1 respectively. The value-added group catalog allows us to detect the dependence of bias on group mass with high significance. It also allows us to compare the measured bias with the theoretically predicted one using the estimated group mass. We find excellent agreement for the two high redshift bins. However, it is lower than the theory by ∼ 3휎 for the lowest redshift bin. Another interesting finding is the significant impact of the thermal Sunyaev Zel’dovich (tSZ). It contaminates the galaxy group-CMB lensing cross-correlation at ∼ 30% level, and must be deprojected first in CMB lensing reconstruction. Key words: large scale structure of Universe – cosmology: observations – cosmic background radiation

1 INTRODUCTION most recently a lensing map over 2100 deg2 by ACT (Darwish et al. 2021). There is a wealth of cosmological information in secondary anisotropies resulting from perturbations to Cosmic Microwave The CMB lensing signal is an integral of all deflections sourced Background (CMB) light after the time of recombination. A partic- by the large-scale structure (LSS) between the last-scattering surface ularly interesting source of secondary anisotropy is the gravitational and the observer, thus it contains rich cosmological information with lensing (Seljak 1996; Seljak & Zaldarriaga 1999; Hu 2000; Hu & its auto-correlation. Meanwhile, it is physically correlated with other Okamoto 2002; Okamoto & Hu 2003), which deflects the paths of LSS tracers such as galaxy and galaxy cluster distribution, and cosmic photons from the last-scattering surface. These deflections, on the shear. On one hand, these cross-correlations are immune to certain order of a few arc-minutes, alter the CMB primary anisotropies by systematics in the auto-correlation, such as the additive errors in re-distributing power across different angular scales and producing both CMB lensing and cosmic shear. On the other hand, they provide a non-Gaussian and anisotropic component to the primordial dis- essential information, such as the redshift information, to improve tribution of temperature anisotropies. Through such non-Gaussian the cosmological applications of CMB lensing. CMB lensing-cosmic and anisotropic features, the gravitational field can be reconstructed arXiv:2109.07387v1 [astro-ph.CO] 15 Sep 2021 shear cross-correlation has been detected by various data sets (Hand from the observed CMB maps (Seljak & Zaldarriaga 1999; Hu & et al. 2015; Liu & Hill 2015; Kirk et al. 2016; Harnois-Déraps et al. Okamoto 2002; Okamoto & Hu 2003; Lewis & Challinor 2006). The 2016; Singh et al. 2017; Harnois-Déraps et al. 2017; Omori et al. first detection of CMB lensing is by the Atacama Cosmology Tele- 2019b; Namikawa et al. 2019; Marques et al. 2020; Robertson et al. scope (ACT) (Das et al. 2011). High signal-to-noise measurements 2021; Heymans et al. 2021; DES Collaboration et al. 2021). Cross- of CMB lensing have been obtained by several collaborations, in- correlations between CMB lensing and galaxies (e.g. recent works cluding a full-sky map of lensing convergence by the Planck satellite 2 by Singh et al.(2017); Schmittfull & Seljak(2018); Bianchini & mission (Planck Collaboration et al. 2020a,b), a 2500 deg CMB Reichardt(2018); Baxter et al.(2018); Raghunathan et al.(2018); Planck lensing map from combined South Pole Telescope (SPT) and Omori et al.(2019a,b); Krolewski et al.(2020); Marques & Bernui temperature data (Omori et al. 2017), a measurement of the lensing (2020); Darwish et al.(2021); Kitanidis & White(2021); Hang et al. B-mode CMB polarization by POLARBEAR (Ade et al. 2014), and (2021); Shekhar Saraf et al.(2021); Fang et al.(2021); Krolewski et al.(2021); Alonso et al.(2021); Yan et al.(2021)), galaxy groups and clusters (Madhavacheril et al. 2015; Baxter et al. 2015; Geach & ★ E-mail: [email protected] Peacock 2017; Baxter et al. 2018; Madhavacheril et al. 2020) have † E-mail: [email protected] also been detected. Recently, the cross-correlation between CMB

© 2021 The Authors SCIENCE CHINA Physics, Mechanics & Astronomy

. Article . XXXX 2021 Vol. xx No. x: 000000 doi: 1

Spectral Analysis of New Black Hole Candidate AT2019wey Observed by NuSTAR

Ye Feng1,2, Xueshan Zhao1,2, Lijun Gou1,2, Yufeng Li1,2, James F. Steiner3, Javier A. García4,5, Yuan Wang1,2, Nan Jia1,2, Zhenxuan Liao1,2, and Huixian Li6

1Key Laboratory for Computational Astrophysics, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China; 2School of Astronomy and Space Sciences, University of Chinese Academy of Sciences, Beijing 100049, China; 3Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA; 4ahill Centre for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA; 5Dr. Karl Remeis-Observatory and Erlangen Centre for Astroparticle Physics, Bamberg, Germany; 6School of Science, Jimei University, Xiamen 361021, China;

Received Month date, Year; accepted Month date, Year

AT2019wey is a new galactic X-ray binary that was first discovered as an optical transient by the Australia Telescope Large Area Survey (ATLAS) on December 7, 2019. AT2019wey consists of a black hole candidate as well as a low-mass companion star (푀star . 0.8푀 ) and is likely to have a short orbital period (푃orb . 8 h). Although AT2019wey began activation in the X-ray band during almost the entire outburst on March 8, 2020, it did not enter the soft state during the entire outburst. In this study, we present a detailed spectral analysis of AT2019wey in the low/hard state during its X-ray outburst on the basis of Nuclear Spectroscopic Telescope Array observations. We obtain tight constraints on several of its important physical parameters by applying the State-of-art relxill relativistic reflection model family. In particular, we determine that the measured inner radius of the accretion disk is = +0.23 most likely to have extended to the innermost stable circular orbit (ISCO) radius, i.e., 푅in 1.38−0.16 푅ISCO. Hence, assuming 푅in=푅ISCO, we find the spin of AT2019wey to be 푎∗ ∼ 0.97, which is close to the extreme and an inner disk inclination angle of 푖 ∼ ◦ 22 . Additionally, according to our adopted models, AT2019wey tends to have a relatively high iron abundance of 퐴Fe ∼ 5 퐴Fe, and a high disk ionization state of log 휉 ∼ 3.4.

NuSTAR, black hole physics, X-rays: binaries–stars: individual: AT2019wey

PACS number(s): 01.55.+b, 04.70.-s, 07.85.Fv

1 Introduction

According to the no-hair theorem of black holes, spin (i.e., angular momentum) is one of the three basic physical quantities of mass, spin, and charge used to determine all properties of a black hole. This parameter is closely related to strong gravity and might be closely associated with jets occurring near the black hole [1]. Because the spin of an isolated black hole is difficult to arXiv:2109.07357v1 [astro-ph.HE] 15 Sep 2021 estimate, analysis of the accreting black hole X-ray binary (BHXRB) is thought to be a more appropriate method for studying black hole physics. In addition to the continuum-fitting method [2], the iron-line profile-fitting method, first proposed by Fabian et al.[3] and also known as the Fe K훼 fitting method, is used to measure the spin of a black hole by modeling the reflection component in a spectrum. This method has increased in maturity during the past decades and is widely applied to dozens of both stellar-mass

1*Corresponding author: [email protected], [email protected], [email protected]

© Science China Press and Springer-Verlag Berlin Heidelberg 2021 phys.scichina.com link.springer.com

A Bidirectional Reflectance Distribution Function for VisorSat

Calibrated with 10,628 Magnitudes from the MMT-9 Database

Anthony Mallama [email protected] 2021 September 15

Abstract

A BRDF for the VisorSat model of Starlink satellites is described. The parameter coefficients were determined by least squares fitting to more than 10,000 magnitudes recorded by the MMT-9 robotic observatory. The BRDF is defined in a satellite-centered coordinate system (SCCS) which corresponds to the physical shape of the spacecraft and to the direction of the Sun. The three parameters of the model in the SCCS are the elevations of the Sun and of MMT-9 along with the azimuth of MMT-9 relative to that of the Sun. The mean VisorSat magnitude at a standardized distance of 1,000 km is 6.84 and the RMS of the distribution around that mean is 1.05. After the magnitudes are adjusted with the BRDF, the RMS reduces to 0.51. The set of MMT-9 observations transformed to the SCCS is available from the author.

1

Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX62

First results on RR Lyrae stars with the TESS space telescope: untangling the connections between mode content, colors and distances L. Molnar´ ,1, 2, 3 A. Bodi´ ,1, 2, 3 A. Pal´ ,1, 3, 4 A. Bhardwaj,5 F–J. Hambsch,6, 7 J. M. Benko˝,1, 2 A. Derekas,8, 9, 1 M. Ebadi,10 M. Joyce,11, 12, 13 A. Hasanzadeh,10 K. Kolenberg,14, 15, 16 M. B. Lund,17 J. M. Nemec,18 H. Netzel,1, 2, 19 C.–C. Ngeow,20 J. Pepper,21 E. Plachy,1, 2 Z. Prudil,22 R. J. Siverd,23 M. Skarka,24, 25 R. Smolec,19 A.´ Sodor´ ,1 S. Sylla,26, 15 P. Szabo´,1, 2, 27 R. Szabo´,1, 2, 3 H. Kjeldsen,28, 29 J. Christensen-Dalsgaard,28 and G. R. Ricker4

1Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, E¨otv¨osLor´andResearch Network (ELKH), Konkoly Thege Mikl´os´ut15-17, H-1121 Budapest, Hungary 2MTA CSFK Lend¨uletNear-Field Cosmology Research Group, 1121, Budapest, Konkoly Thege Mikl´os´ut15-17, Hungary 3ELTE E¨otv¨osLor´andUniversity, Institute of Physics, 1117, P´azm´anyP´eters´et´any1/A, Budapest, Hungary 4MIT Kavli Institute for Astrophysics and Space Research, 70 Vassar St, Cambridge, MA 02109, USA 5Korea Astronomy and Space Science Institute, Daedeokdae-ro 776, Yuseong-gu, Daejeon 34055, Republic of Korea 6Vereniging Voor Sterrenkunde (VVS), Oostmeers 122 C, 8000 Brugge, Belgium 7Bundesdeutsche Arbeitsgemeinschaft f¨urVer¨anderlicheSterne e.V. (BAV), Munsterdamm 90, D-12169 Berlin, Germany 8ELTE E¨otv¨os Lor´andUniversity, Gothard Astrophysical Observatory, Szent Imre h. u. 112, 9700, Szombathely, Hungary 9MTA-ELTE Exoplanet Research Group, Szent Imre h. u. 112, 9700 Szombathely, Hungary 10Institute of Geophysics, University of Tehran, Tehran, Iran 11Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 12Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia 13ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia 14Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, B-3001 Heverlee, Belgium 15Physics Department, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 16Physics and Astronomy Department, Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussel 17Caltech/IPAC, 1200 E. California Blvd. Pasadena, CA 91125, USA 18Department of Physics & Astronomy, Camosun College, Victoria, BC, Canada 19Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warsaw, Poland 20Graduate Institute of Astronomy, National Central University, Jhongli 32001, Taiwan 21Department of Physics, Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015, USA 22Astronomisches Rechen-Institut, Zentrum f¨urAstronomie der Univ¨ersitatHeidelberg, M¨onchhofstr. 12-14, D-69120 Heidelberg, Germany 23Gemini Observatory/NSF’s NOIRLab, 670 N. A‘ohoku Place, Hilo, HI, 96720, USA 24Department of Theoretical Physics and Astrophysics, Masaryk University, Kotl´aˇrsk´a2, 61137 Brno, Czech Republic 25Astronomical Institute, Czech Academy of Sciences, Friˇcova298, 25165, Ondˇrejov,Czech Republic 26Universit`eCheikh Anta Diop Dakar-Fann, S`en`egal 27Magdalene College, University of Cambridge, CB3 0AG, Cambridge, UK 28Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 29Astronomical Observatory, Institute of Theoretical Physics and Astronomy, Vilnius University, Saul˙etekioav. 3, 10257, Vilnius, Lithuania

(Dated: Submitted: 23/06/2021; Accepted: 14/09/2021) arXiv:2109.07329v1 [astro-ph.SR] 15 Sep 2021

ABSTRACT The TESS space telescope is collecting continuous, high-precision optical photometry of stars throughout the sky, including thousands of RR Lyrae stars. In this paper, we present results for an initial sample of 118 nearby RR Lyrae stars observed in TESS Sectors 1 and 2. We use differential- image photometry to generate light curves and analyse their mode content and modulation properties. We combine accurate light curve parameters from TESS with parallax and color information from the

Corresponding author: L. Moln´ar [email protected], [email protected] Astronomy & Astrophysics manuscript no. paper_main ©ESO 2021 September 16, 2021

Euclid: Constraining ensemble photometric redshift distributions with stacked spectroscopy? M.S. Cagliari1??, B.R. Granett2???, L. Guzzo1,2,3, M. Bolzonella4, L. Pozzetti4, I. Tutusaus5,6, S. Camera7,8,9, A. Amara10, N. Auricchio4, R. Bender11,12, C. Bodendorf12, D. Bonino9, E. Branchini13,14, M. Brescia15, V. Capobianco9, C. Carbone16, J. Carretero17, F.J. Castander5,6, M. Castellano18, S. Cavuoti15,19,20, A. Cimatti21,22, R. Cledassou23,24, G. Congedo25, C.J. Conselice26, L. Conversi27,28, Y. Copin29, L. Corcione9, M. Cropper30, H. Degaudenzi31, M. Douspis32, F. Dubath31, S. Dusini33, A. Ealet29, S. Ferriol29, N. Fourmanoit34, M. Frailis35, E. Franceschi4, P. Franzetti16, B. Garilli16, C. Giocoli36,37, A. Grazian33, F. Grupp11,12, S.V.H. Haugan38, H. Hoekstra39, W. Holmes40, F. Hormuth41,42, P. Hudelot43, K. Jahnke42, S. Kermiche44, A. Kiessling40, M. Kilbinger45, T. Kitching30, M. Kümmel11, M. Kunz46, H. Kurki-Suonio47, S. Ligori9, P.B. Lilje38, I. Lloro48, E. Maiorano4, O. Mansutti35, O. Marggraf49, K. Markovic40, R. Massey50, M. Meneghetti4,51,52, E. Merlin18, G. Meylan53, M. Moresco4,21, L. Moscardini4,21,52, S.M. Niemi54, C. Padilla17, S. Paltani31, F. Pasian35, K. Pedersen55, W.J. Percival56,57,58, V. Pettorino45, S. Pires45, M. Poncet24, L. Popa59, F. Raison12, R. Rebolo60,61, J. Rhodes40, H.-W. Rix42, M. Roncarelli4,21, E. Rossetti21, R. Saglia11,12, R. Scaramella18,62, P. Schneider49, M. Scodeggio16, A. Secroun44, G. Seidel42, S. Serrano5,6, C. Sirignano63,64, G. Sirri52, D. Tavagnacco35, A.N. Taylor25, I. Tereno65,66, R. Toledo-Moreo67, E.A. Valentijn68, L. Valenziano4,52, Y. Wang69, N. Welikala25, J. Weller11,12, G. Zamorani4, J. Zoubian44, M. Baldi4,21,52, R. Farinelli70, E. Medinaceli4, S. Mei71, G. Polenta72, E. Romelli35, T. Vassallo11, A. Humphrey73 (Affiliations can be found after the references)

ABSTRACT

Context. The ESA Euclid mission will produce photometric galaxy samples over 15 000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. Aims. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment the broad-band photometric infor- mation to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. Methods. We model the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. Results. We select mock photometric galaxy samples at redshift z > 1 using the self-organizing map algorithm. Considering the idealized case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of the Euclid near-infrared spectroscopy. Key words. method: data analysis – method: statistical – galaxies: distances and redshifts – large-scale structure of Universe

1. Introduction spectral resolution. Imaging surveys are limited to photomet- ric measurements such as broadband colours to infer the red- arXiv:2109.07303v1 [astro-ph.CO] 15 Sep 2021 The next generation of photometric surveys will produce un- shifts of galaxies (Connolly et al. 1995; Bolzonella et al. 2000; precedented galaxy statistics that will fuel large-scale structure Benítez 2000). The minimum error in a photometric redshift es- studies (LSST Science Collaboration et al. 2009; Laureijs et al. timate with optical and near-infrared broadband photometry is 2011; Benitez et al. 2014; Dark Energy Survey Collaboration σz/(1+z) ∼ 0.02 due to fundamental degeneracies in the colour- et al. 2016). Compared with their spectroscopic counterparts (Le redshift parameter space (Salvato et al. 2019). Nevertheless, with Fèvre et al. 2005; Driver et al. 2011; Guzzo et al. 2014; DESI the promise of large sample sizes this precision is often accept- Collaboration et al. 2016), photometric surveys go deeper and able for large-scale structure studies based on galaxy clustering faster; however, the surveying efficiency comes at the cost of and weak lensing analyses. The redshift of individual galaxies is

? not required for these analyses, but instead precise knowledge of This paper is published on behalf of the Euclid Consortium. the redshift distribution of the sample is needed to properly inter- ?? e-mail: [email protected] ??? e-mail: [email protected] pret the statistics. A systematic error in the redshift distribution

Article number, page 1 of 15 arXiv:2109.07290v1 [astro-ph.HE] 15 Sep 2021 rdcino Mrdaindpnso h rwhmdl fB n it and BH of models growth the on depends radiation EM of production ocuin r eyitrsig eew ilts hs w growt two radiation these EM test no 13 will ta we expects [12, data. we Here one in Insight-HXMT if Therefore, interesting. analysed very BH, time. further are a finite was conclusions within onto topic pre-e BH falling This a the matters into metric. of ar for background metric it process However, the EM stationary the 10]. on of a [9, describe Generation BHs uses to merging one from metric [8]. ga If radiation BHs short EM BH. strong charged events, of predicts of GW models g ra mergers growth fast a [7]. from the a in mergers expected alike g happens BH flare be the merger stellar-mass radio also transient BBHs on from a if follow depends by that possibly accompanied merger argued be BBHs can been in process has radiation It EM environment. of obs generation counterpart general, its and wave weak gravitational the joint [3], parameters from unde cosmological [5] to [2], violations important dimension also extra is probe the band can th in (EM) found gravitational electromagnetic been in have test counterpart relativity to general from opportunities deviations provide significant GWs observations. eto I,w ettegot oeso Hb onigLG obser LIGO jointing by BH V. of section in models results growth our the summarize test briefly we III, section † 3 ∗ ‡ lcrncades [email protected] address: Electronic ainlLclJitEgneigLbrtr fNwEnergy New of Laboratory Engineering Joint National-Local [email protected] author Corresponding h eeto fgaiainlwvs(W)b h IOCollaborat LIGO the by (GWs) waves gravitational of detection The hte iaybakhls(Bs egrpoue electromagne produces merger (BBHs) holes black binary Whether h eto h ae sognsda olw.I eto I ebriefl we II, section In follows. as organised is paper the of rest The h rwho Hi eyipratrsac oi,sm aspects some topic, research important very a is BH of growth The ettegot oeso lc oeb onigLG n Insi and LIGO jointing by hole black of models growth the Test 2 ee e a fOtcEetoi nomto n Material and Information Optic-Electronic of Lab Key Hebei ASnumbers: m PACS stationary test the We find time. We observations. short Insight-HXMT very observations. and within Insight-HXMT hole and black a LIGO horizo final metric, event the time-dependent of the the horizon In outside observer. accumulated external be distant must matter falling ogJaYang Rong-Jia h rwhmdlo lc oei tl otoesa topic controversial a still is hole black of model growth The 1 olg fPyia cec n ehooy ee Universi Hebei Technology, and Science Physical of College 4 5 e aoaoyo ihpiiinCmuainadApplicat and Computation High-pricision of Laboratory Key il hoyo ee rvne ee nvriy adn 0 Baoding University, Hebei Province, Hebei of Theory Field e aoaoyo atceAtohsc,Isiueo High of Institute Astrophysics, Particle of Laboratory Key ∗ , ,2 ,4, 3, 2, 1, hns cdm fSine,Biig104,China 100049, Beijing Sciences, of Academy Chinese I RWHMDL FBAKHOLE BLACK OF MODELS GROWTH II. † hwiGuo, Zhiwei .INTRODUCTION I. observations 1, † agagZheng, Yaoguang hteeti eie,HbiUiest,Boig071002, Baoding University, Hebei Devices, Photoelectric oeso Hb onigLG bevtosand observations LIGO jointing by BH of models h li-aln atrcnfl noteevent the into fall can matter in-falling ll ,HbiUiest,Boig010,China 071002, Baoding University, Hebei s, ervations. oisbyn eea eaiiy pt o,no now, to Up relativity. general beyond eories ekfil eie naG vn,ietfigthe identifying event, GW a In regime. weak-field sadtefl hsclpoes o xml,we example, For process. physical full the rstand i us 6.I a hw htE inl can signals EM that shown was It [6]. burst dio qiaec rnil 4,adLrnzinvariance Lorentz and [4], principle equivalence itn Ht ecietegot fB,one BH, of growth the describe to BH xisting udi 1]ta esol s time-dependent use should we that [11] in gued ftebakhl,a lce ya by clocked as hole, black the of n .I hsmdl l nfligmtescnfall can matters in-falling all model, this In ]. eit con h ffc fi-aln matter in-falling of effect the account into ke ain n nih-XTdt.Fnly we Finally, data. Insight-HXMT and vations hs w rwhmdl yjoining by models growth two these i onepr ssila pnpolm In problem. open an still is counterpart tic niomn.I hsscin ewl briefly will we section, this In environment. s o 1 pn e idwfrastronomical for window new a opens [1] ion ntesainr erc l in- all metric, stationary the In . fi r tl nla.I Bsmre,the merger, BBHs In unclear. still are it of sosds rudasprasv H the BH, supermassive a around disk aseous rmmrigBs hs w contradictory two These BHs. merging from aito nBH egras eed on depends also merger BBHs in radiation eiwtetogot oeso H In BH. of models growth two the review y m-a uss n atrdobrt can bursts radio fast and bursts, mma-ray 1 ot oeso lc oe(H n its and (BH) hole black of models rowth y adn 702 China 071002, Baoding ty, dli nossetwt LIGO with inconsistent is odel n hagNnZhang Shuang-Nan and nryPhysics, Energy o fQuantum of ion 10,China 71002, ght-HXMT ‡ 5 China arXiv:2109.07235v1 [astro-ph.HE] 15 Sep 2021 iheeg bv 0GVpoaaigoe h omlgcld cosmological the over propagating GeV 10 d above process energy This with production. pair electron-positron an extragalac to the ing with interaction the through absorption the e words: Key of 21)sostebs efrac o h rpgto of model propagation particu the EBL In for the performance method. to best BP corresponding the the layer shows to (2010) hidden superior the be in to neurons found 40 is perfo RBF The the algorithms. ing (RBF) st function employs radial-basis model and EBL (BP) each op and for inputs ANN-performance as of redshift energy, timization using ANN the divided of randomly dif testing are and three estimates using depth 1 to optical 0.01 These range models. the in redshifts at sources the ∗ cosmic of propagation (ANN network the neural probe artificial an of potential the explore We Abstract rica erlnetwork neural artificial os eetmt h pia et ausfor values depth optical the estimate We tons. eed nthe on depends a teuto scaatrzdb hsclqatt cal quantity physical a by characterized is attenuation ray rpitsbitdt e srnm 6Spebr2021 September 16 Astronomy New to submitted Preprint a b mi address: Email author. Corresponding γ srpyia cecsDvso,Bah tmcRsac C Research Atomic Bhabha Division, Sciences Astrophysical rica erlNtok o omcgamma-ray cosmic for Networks Neural Artificial hsc eatet nvriyo h reSae Bloemfon State, Free the of University Department, Physics ry mte rmadsatsuc iebaa oteobserv the to blazar like source distant a from emitted -rays am-as xrglci akrudlgt photon-phot light; background extragalactic gamma-rays; .K Singh K. K. γ rypoo nry esito h ore n est fth of density and source, the of redshift energy, photon -ray [email protected] rpgto nteUniverse the in propagation a , γ ∗ rypoosi h xrglci nvre h journey The Universe. extragalactic the in photons -ray .K Dhar K. V. , India γ ryeege bv 0GVeitdfrom emitted GeV 10 above energies -ray K .Singh). K. (K. a .J Meintjes J. P. , i akrudlgt(B) lead- (EBL), light background tic led mntsfrgmaryphotons ray gamma for ominates γ eetadms rmsn EBL promising most and ferent ry nteUniverse. the in -rays mneo h N oe us- model ANN the of rmance ae ehditliec to intelligence method based ) notodt esfrtraining for sets data two into pia depth optical ia et sotu.Teop- The output. as depth tical ra h at sipddby impeded is Earth the at er en-90,SuhAfrica South 9300, - tein nr,Mma-40085, 400 Mumbai- entre, sacs h fetof effect The istances. nadback-propagation andard rpsdb ik tal. et Finke by proposed a,teRFANwith RBF-ANN the lar, ninteraction; on b hc strongly which , B pho- EBL e γ - DRAFTVERSION SEPTEMBER 16, 2021 Typeset using LATEX default style in AASTeX63

In search of infall motion in molecular clumps III: HCO+ (1-0) and H13CO+ (1-0) mapping observations toward the confirmed infall sources

YANG YANG,1, 2, 3 ZHIBO JIANG,1, 2, 4 ZHIWEI CHEN,1, 2 YIPING AO,1, 2, 3 AND SHULING YU1, 2, 3

1Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023, China 2Key Laboratory of Radio Astronomy, Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023, China 3University of Science and Technology of China, Chinese Academy of Sciences, Hefei 230026, China 4Center astronomy and Space Sciences, China Three Gorges University, Yichang 443002, China

(Received April 26, 2021; Revised August 26, 2021; Accepted August 30, 2021) Submitted to ApJ

ABSTRACT The study of infall motion helps us to understand the initial stages of star formation. In this paper, we use the IRAM 30-m telescope to make mapping observations of 24 infall sources confirmed in previous work. The lines we use to track gas infall motions are HCO+ (1-0) and H13CO+ (1-0). All 24 sources show HCO+ emissions, while 18 sources show H13CO+ emissions. The HCO+ integrated intensity maps of 17 sources show clear clumpy structures; for the H13CO+ line, 15 sources show clumpy structures. We estimated the column + 13 + + density of HCO and H CO using the RADEX radiation transfer code, and the obtained [HCO ]/[H2] and [H13CO+]/[HCO+] of these sources are about 10−11–10−7 and 10−3–1, respectively. Based on the asymmetry of the line profile of the HCO+, we distinguish these sources: 19 sources show blue asymmetric profiles, and the other sources show red profiles or symmetric peak profiles. For eight sources that have double-peaked blue line profiles and signal-to-noise ratios greater than 10, the RATRAN model is used to fit their HCO+ (1-0) lines, −1 and to estimate their infall parameters. The mean Vin of these sources are 0.3 – 1.3 km s , and the M˙ in are −3 −4 −1 about 10 – 10 M yr , which are consistent with the results of intermediate or massive star formation in −1 −3 previous studies. The Vin estimated from the Myers model are 0.1 – 1.6 km s , and the M˙ in are within 10 – −5 −1 10 M yr . In addition, some identified infall sources show other star-forming activities, such as outflows and maser emissions. Especially for those sources with a double-peaked blue asymmetric profile, most of them have both infall and outflow evidence.

Keywords: Star formation — Interstellar medium — Molecular clouds — Collapsing clouds

1. INTRODUCTION The gravitational infall in molecular clouds mainly occurs in the beginning stage of star formation (e.g. Bachiller 1996), and it is an important part of the star-forming process. Stars are believed to form through the inside-out gravitational collapse of dense molecular cloud cores (e.g. Shu et al. 1987). In the subsequent evolution, the accretion flow has been maintained to feed the forming stars (e.g. Allen et al. 2004). The study of this process helps us to understand the initial stage of the star-formation arXiv:2109.07221v1 [astro-ph.GA] 15 Sep 2021 process. There exist several studies focusing on the infall motions: Mardones et al.(1997) used the IRAM 30-m, SEST 15-m, and Haystack 37-m radio telescopes to observe some candidate protostars, and reported 15 infall candidates. Most of them are associated with clumps and have been found to have turbulent motions. Klaassen & Wilson(2008) observed seven massive star-forming regions (SFRs), and found that four of them have infall signatures based on their HCO+ emission profiles. He et al. + + (2015) selected more than 400 compact sources from the MALT90 survey, and used the HCO (1-0), HNC (1-0), and N2H (1-0) data to identify 131 infall candidates. Spitzer’s mid-infrared data were also used to limit the evolutionary stages of the candidates as a supplement, and to analyze the relationship between the infall motion and evolution stage of these candidates.

Corresponding author: Zhibo Jiang [email protected] Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX63

Multipoint ICME events during the first year of combined Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A observations

Christian Mostl¨ ,1, 2 Andreas J. Weiss,1, 2, 3 Martin A. Reiss,1, 2 Tanja Amerstorfer,1 Rachel L. Bailey,4, 1 Jurgen¨ Hinterreiter,1, 3 Maike Bauer,1, 3 David Barnes,5 Jackie A. Davies,5 Richard A. Harrison,5 Johan L. Freiherr von Forstner,6 Emma E. Davies,7 Daniel Heyner,8 Tim Horbury,9 and Stuart D. Bale10

1Space Research Institute, Austrian Academy of Sciences, Schmiedlstraße 6, 8042 Graz, Austria 2Institute of Geodesy, Graz University of Technology, Steyrergasse 30, 8010 Graz, Austria 3Institute of Physics, University of Graz, Universit¨atsplatz5, 8010 Graz, Austria 4Zentralanstalt f¨ur Meteorologie und Geodynamik, Hohe Warte 38, 1190 Vienna, Austria 5RAL Space, Rutherford Appleton Laboratory, Harwell Campus, UK 6Institute of Experimental and Applied Physics, University of Kiel, Kiel, Germany 7Space Science Center, Institute for the Study of Earth, Oceans, and Space, University of New Hampshire, Durham, NH, USA 8TU Braunschweig, Germany 9Imperial College, UK 10Physics Department and Space Sciences Laboratory, University of California, Berkeley, CA, USA

(Received 2021 September 15; Accepted -) Submitted to ApJL

ABSTRACT We report the result of the first search for multipoint in situ and imaging observations of inter- planetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter data in April 2020 to April 2021. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A, in connection with coronagraph and heliospheric imaging observations from STEREO-Ahead/SECCHI and SOHO/LASCO. We identify ICME events that could be unambigu- ously followed with the STEREO-A heliospheric imagers during their interplanetary propagation to their impact at the aforementioned spacecraft, and look for events where the same ICME is seen in situ by widely separated spacecraft. We highlight two events: (1) a small streamer blowout CME on 2020 June 23 observed with a triple lineup by Parker Solar Probe, BepiColombo and Wind, guided by imaging with STEREO-A, and (2) the first fast CME of solar cycle 25 (≈ 1600 km s−1) on 2020 Nov 29 observed in situ by Parker Solar Probe and STEREO-A. These results are useful for modeling the magnetic structure of ICMEs, the interplanetary evolution and global shape of their flux ropes and shocks, and for studying the propagation of solar energetic particles. The combined data from these missions is already turning out to be a treasure trove for space weather research and is expected to become even more valuable with an increasing number of ICME events expected during the rise and maximum of solar cycle 25. arXiv:2109.07200v1 [astro-ph.SR] 15 Sep 2021 Keywords: editorials, notices — miscellaneous — catalogs — surveys

1. INTRODUCTION tions observed in the solar corona and their signatures Two major problems in space weather research con- measured in situ by spacecraft in interplanetary space, cern (1) establishing a connection between solar erup- and (2) finding events where more than one spacecraft experiences the eruptive event at different locations in the inner heliosphere. Magnetic field observations by Corresponding author: Christian M¨ostl magnetometers are currently the only available method [email protected] to measure the magnetic field structure of the flux ropes Noname manuscript No. (will be inserted by the editor)

UV lines as a tracers for the XUV-fluxes of stars and the PLATOspec project

Eike Wolf Guenther · Petr Kabath · Leonardo Vanzi

Received: date / Accepted: date

Abstract Observations in the UV-regime are very important for exoplanet research, because many diagnostically important lines for studying stellar ac- tivity are in this regime. Studying stellar activity is not only important because of its negative effects on the determination planetary parameters, but also be- cause the XUV-radiation from the host stars affects the photochemistry and the erosion of planetary atmospheres . Unfortunately, the XUV-region is only accessible from space. However, since the XUV-radiation is correlated with

the Ca II H&K-lines, we can use these lines to study the XUV radiation in-

directly. The Ca II H&K-lines for relatively bright stars can be observed with PLATOspec, a new high-resolution echelle spectrograph in development for the ESO 1.5m telescope at La Silla. One advantage compared to instruments on larger telescopes will be that large programs can be carried out. There will be two modes for obtaining precise RV-measurements. In the future, the CUBES instrument on the VLT will be able to study the same lines to probe the XUV-radiation in much fainter targets.

Keywords extrasolar planet · stellar activity · CaII,HK-lines · PLATO · ARIEL

E.W. Guenther Th¨uringerLandessterwarte Tautenburg, 07778 Tautenburg, Germany Tel.: +49-36427-86355 Fax: +49-36427-86329 E-mail: [email protected] Petr Kab´ath Astronomical Institute of the Czech Academy of Sciences, Friˇcova 298, 25166, Ondˇrejov, Czech Republic Leonardo Vanzi Centre of Astro-Engineering, Pontificia Universidad Catolica de Chile, Av. Vicu˜naMackenna 4860, Santiago, Chile, and Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Av. Vicu˜naMackenna 4860, Santiago, Chile arXiv:2109.07168v1 [astro-ph.EP] 15 Sep 2021 MNRAS 000,1–12 (2021) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

Compensation of differential dispersion: application to multiband stellar interferometry.

Cyril Pannetier1,2?, Denis Mourard1, Frédéric Cassaing2,3, Stéphane Lagarde1, Jean-Baptiste Le Bouquin4, John Monnier5, Judit Sturmann6, and Theo Ten Brummelaar6. 1Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, Bd de l’Observatoire, CS 34229, 06304 Nice cedex 4, France 2DOTA, ONERA, Université Paris Saclay, F-92322 Châtillon, France 3Now at DTIS, ONERA, Université Paris Saclay, F-91123 Palaiseau, France 4Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France 5University of Michigan, Ann Arbor, MI 48109 USA 6CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT With the aim of pushing the limiting magnitude of interferometric instruments, the need for wide-band detection channels and for a coordinated operation of different instruments has considerably grown in the field of long-baseline interferometry. For this reason, the Center for High Angular Resolution Astronomy (CHARA), an array of six telescopes, requires a new configuration of longitudinal dispersion compensators to keep the fringe contrast above 95% simultaneously in all spectral bands, while preserving the transmission above 85 %. In this paper, we propose a new method for defining the longitudinal dispersion compensators (LDC) suited for multi-band observations. A literal approximation of the contrast loss resulting from the dispersion residues enables us to define a general criterion for fringe contrast maximisation on several bands simultaneously. The optimization of this criterion leads to a simple solution with only two LDC stages per arm and existing differential delay lines, to the glass choice and a simple linear formula for thickness control of all these media. A refined criterion can also take into account glass transmission. After presenting this criterion, we give the optimal solution (medium, configuration) and its expected performance for the planned observing modes on CHARA. Key words: atmospheric effects – instrumentation: interferometers – techniques: interferometric – methods: observational

1 INTRODUCTION This problem has been identified and solved from the very first stellar interferometer with the introduction of LDC, usually made Ground-based long-baseline optical interferometers combine the of pieces of glass of variable thickness to compensate the air thick- light from an array of telescopes distant by tens to hundreds of ness (Labeyrie 1975). Another solution, currently being implemented meters. The combination of the electromagnetic fields creates an at Magdalena Ridge Optical Interferometer (Creech-Eakman et al. interferometric pattern whose characteristics depend on the inter- 2018), is to use evacuated delay lines. ferometer configuration and the source properties. It gives access to the angular resolution of a single telescope of diameter equal to The recent years have seen important developments of the inter- the maximum distance between the telescopes in the array. Because ferometric instrumentation (Gravity Collaboration et al. 2017; Lopez they are operated on broadband sources, performance of these in- et al. 2014; Anugu et al. 2020; Monnier et al. 2018; Mourard et al. arXiv:2109.07163v1 [astro-ph.IM] 15 Sep 2021 terferometers intended to measure the spatial coherence is critically 2018). For measuring visibilities on fainter objects, the need for large affected by temporal coherence. The current generation of optical spectral channels instruments (Gravity Collaboration et al. 2017; long-baseline interferometers, schematically represented in Fig.1, Mourard et al. 2018) and for a coordinated operation of different in- uses main delay lines (MDL) to equalize the optical paths, which struments (Matter et al. 2010; Mourard et al. 2017; Lopez et al. 2014) achromatic difference is mostly introduced in vacuum outside the at- has considerably grown. Indeed, this coordination not only enables mosphere. However, if this equalization is done in air, the dispersion to track fringes in one band while integrating in all other bands but law of this medium introduces a dispersion at the interferometric also increases the available observing time of all instruments. For focus. The equalization of the optical path becomes chromatic and example, the Gravity Fringe Tracker (Lacour et al. 2019) working in not only geometric: the coherence envelop multiplying the interfer- the K band is now used to stabilize the fringes of the instrument MA- ometric fringes is chromatically shifted and the fringes are blurred TISSE (Lopez et al. 2014) in the L and M band, with the introduction within a spectral channel as it gets wider. of a correction of the chromatic group-delay with the internal delay lines of MATISSE. CHARA is located at Mount Wilson (California, USA). It has ? E-mail: [email protected] (CP) been designed by the University of Georgia and entered into service

© 2021 The Authors arXiv:2109.07156v1 [astro-ph.GA] 15 Sep 2021 ap,ec hsapoc fcoeiseto fglci sy galactic of inspection close of approach This etc. warps, ihsia r ic nl,cnrlvlct iprinan dispersion velocity central angle, pitch arm spiral with ihhl rpris ..ms,sz,ec,cudgenerate could etc., size, mass, i.e. properties, halo with eoiy oee,sm eerhr aepoie evidenc provided have researchers some However, velocity. ncnieigec atclrglx,oecudosreth observe could one galaxy, particular each considering In Introduction 1 ieo eerhlast aatcsaiglw.I sreason ori is the It understanding better laws. toward scaling guidelines galactic valuable to as leads research of of properties global line of investigation system the various be of could behavior approach dynamical the method, this Using apeo 8sia aaisadrpr orltosbetwe correlations report and galaxies spiral 48 of sample a conne could one which [ Navarro by scenario characters. observational suitable a defining is issue [ systems [ Brandt by b link proposed a is describes rule primarily relation cons notab Tully-Fisher usually is the is relation luminosity that this NIR in since scatter addition, In the that infrared. is BTFR of favour in erro of sources the of phys issue the the studying is in rules employed scaling var be the on deriving could depend evolution to the expected of indeed t trend is with parameters compared scaling be the can of results the thus, dependent; measurable redshift easily some with astronomers provide relations ewe aaysms n iewt hs fhl r on to found are der halo in of importance those critical with of size and also mass is galaxy’s halo between matter dark cold of [ xml,i sepce ntecneto h tnadmdlo model standard the of context the in expected is it example, ue r eie ae ntosml supin eae t related - assumptions simple [ two Brandt on and based derived are rules orltstebroi aswt hrceitcvelocit characteristic a with mass baryonic the correlates 1 rvdsfu necnetdsaigrltosbtenb between relations scaling interconnected four provides ] * -al [email protected] E-mail: n ftems cuaesaigrlsi eotdb ul an Tully by reported is rules scaling accurate most the of One h rgno aatcsaigrlsi osdrdt eroot be to considered is rules scaling galactic of origin The h iblt ftesaigrltosi etdadi ssho is it and tested is relations scaling the of viability The velocity R bevtoa aa ehv end ae nteproperties the on based defined, have we data, observational vlto fglci rpriswt esit hsbehav works. future This in model redshift. dynamical modified with present systems. properties the galactic of of characteristics evolution general pl the relation of Tully-Fisher baryonic description of set the that concluded eainbtensraedniyadcrua eoiyo d rela of velocity scaling circular scalin and These density The surface between conditions. relation model. stability dynamical and modified equilibrium applying a of context the eq eateto hsc,Fros nvriyo aha,P.O Mashhad, of University Ferdowsi Physics, of Department 11 cuzglci cln ue [ rules scaling galactic Schulz and , 10 12 aatcSaigRlsi oie yaia Model Dynamical Modified a in Rules Scaling Galactic V V eq eain ysreigteltrtr,oemyesl n o find easily may one literature, the surveying By relation. ] eq ]. rvd h lss t otetertclpeitos eh We predictions. theoretical the to fits closest the provide fsses yepoigteemaue flnt n veloci and length of measures these employing By systems. of 10 sacreainbtenms,vlct n ie( size and velocity mass, between correlation a as ] 14 1 rus ae nNbd n yrdnmclsmltos t simulations, hydrodynamical and N-body on based argues, ] ,wihicuebroi ul-ihrrlto,hv been have relation, Tully-Fisher baryonic include which ], etme 6 2021 16, September osi Shenavar Hossein ,ie h a eoiy fglci ytm [ systems galactic of velocity, flat the i.e. y, tenms n oainlvlct ( velocity rotational and mass etween aso h eta lc hole. black central the of mass d Abstract s aais ocektevldt ftesaigrelations scaling the of validity the check To galaxies. isk drda h etmaueo h tla as n ol conc could one mass, stellar the of measure best the as idered nta,cmae osm te esrso eghadveloci and length of measures other some to compared that, wn sms-ierlto rms-eoiy)poie nappro an provides ) mass-velocity or ( relation mass-size us aaissc slmnst,ms,vlct,sz,ec Thi etc. size, velocity, mass, luminosity, as such galaxies ue r lootie ydmninlaayi fa integ an of analysis dimensional by obtained also are rules g i n vlto fglxe.A ti xlie y[ by explained is it As galaxies. of evolution and gin si a endsusdb [ by discussed been has it as r ttsial oiae,caatr fglci systems galactic of characters motivated, statistically , tdfeetssesdslydsic etrssc sbars as such features distinct display systems different at aty ti hw htteesaigrltospeitace a predict relations scaling these that shown is it Lastly, bet ru htacrt aatcsaiglw ol eus be could laws scaling galactic accurate that argue to able o rvdsacac oeaietecsi vlto fthe of evolution cosmic the examine to chance a provides ior ysalri erifae NR oprdt lrvoe o ultraviolet to compared (NIR) near-infrared in smaller ly eotoeo ueia iuain.Mroe,teevoluti the Moreover, simulations. numerical of outcome he nspoto ayncTlyFse eain(TR whic (BTFR) relation Tully-Fisher baryonic a of support in e nmxmmds oainlvlct adhl’ akmatter dark halo’s (and velocity rotational disk maximum en aatcsaigrltos[ relations scaling galactic 1 ftemdl h rpreulbimsize equilibrium proper the model, the of roi as aisadvlct fglxe.Teescalin These galaxies. of velocity and radius mass, aryonic vn h cln ue.I hs iuain,tecorrelat the simulations, these In rules. scaling the iving ih eudrto ndti.Hwvr nte productiv another However, detail. in understood be might s osmcaim fglx omto;teeoe h genera the therefore, formation; galaxy of mechanisms ious h hnmnlg fglxe n hyicueBTFR include they and - galaxies of phenomenology the o ttepoete ftedr n ayncmte otepure the to matter baryonic and dark the of properties the ct tm ol edt nitiaecasfiaino galaxies. of classification intricate an to lead would stems c fglx omto.Hwvr oiatcalnein challenge dominant a However, formation. galaxy of ics ehgl o-ier e loDvse l [ al. et Davis also See non-linear. highly be omlg htsseai aito fglci features galactic of variation systematic that cosmology f ihr[ Fisher d in r eie yepoigtevra hoe and theorem virial the employing by derived are tions dwti h oenn yaiso h ytm.For systems. the of dynamics governing the within ed v oprdorrslswt ro ok n have and works prior with results our compared ave * y PR aa[ data SPARC ty, hrwrsmninn iia eut nvarious in results similar mentioning works ther 5 sarlto ewe uioiyadcircular and luminosity between relation a as ] o 46 aha,Iran. Mashhad, 1436, Box . M − 13 4 2 V extensively. ] sue oaayeteresults. the analyze to used is ] , uvydi hswr within work this in surveyed − 14 R .I hscs,a essential an case, this In ]. M fglxe.As,Schulz Also, galaxies. of ) − a h efsmlrnature self-similar the hat 6 R V , eq 7 [ ) , n equilibrium and 8 1 , .Aohrscaling Another ]. 9 .Teargument The ]. 15 ae on based h survey who ] priate 3 rtain hc are which ,scaling ], ral second s ty, boxes, , far- r lude ions on ed g h e ) l Astronomy & Astrophysics manuscript no. Canopus_2021 ©ESO 2021 September 16, 2021

Refined fundamental parameters of Canopus from combined near-IR interferometry and spectral energy distribution ? A. Domiciano de Souza1,??, J. Zorec2, F. Millour1, J.-B. Le Bouquin3, A. Spang1, and F. Vakili1

1 Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, France 2 Sorbonne Université CNRS, UMR 7095, Institut d’Astrophysique de Paris, F-75014 Paris, France 3 Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France

Received ...; accepted ...

ABSTRACT

Context. Canopus, the brightest and closest yellow supergiant to our Solar System, offers a unique laboratory for understanding the physics of evolved massive stars. Aims. We aim at quantitatively exploring a large space of fundamental parameters of Canopus based on the combined analysis of its spectral energy distribution (SED) and optical-IR long baseline interferometry. Methods. We use the most recent high resolution near-IR data from the VLTI focal beam combiners PIONIER (H and K bands) and AMBER (K band), together with precise spectrophotometric measures that cover the SED of Canopus, from the UV to the IR, taken from ground and space observatories. Results. The accurate and precise PIONIER data allowed us to simultaneously measure the angular diameter and the limb darkening (LD) profile using different analytical laws. We found that the power-law LD, being also in agreement with predictions from stellar atmosphere models, reproduces the interferometric data well. For this model we measured an angular diameter of 7.184 ± 0.0017 ± 0.029 mas and an LD coefficient of 0.1438 ± 0.0015, which are respectively & 5 and ∼ 15 − 25 more precise than in our previous A&A paper on Canopus from 2008. From a dedicated analysis of the interferometric data, we also provide new constraints on the putative presence of weak surface inhomogeneities. Additionally, we analyzed the SED in a innovative way by simultaneously fitting the reddening-related parameters and the stellar effective temperature and gravity. We find that a model based on two effective temperatures is much better at reproducing the whole SED, from which we derived several parameters, including a new bolometric −6 −2 −1 flux estimate: fbol = (59.22 ± 2.45) × 10 erg cm s . We were also able to estimate the stellar mass from these measurements, and it is shown to be in agreement with additional predictions from evolutionary models, from which we inferred the age of Canopus as well. Conclusions. The Canopus angular diameter and LD measured in this work with PIONIER are the most precise to date, with a direct impact on several related fundamental parameters. Moreover, thanks to our joint analysis, we were able to determine a set of fundamental parameters that simultaneously reproduces both high-precision interferometric data and a good quality SED and, at the same time, agrees with stellar evolution models. This refined set of fundamental parameters constitutes a careful balance between the different methodologies used, providing invaluable observationally based constraints to models of stellar structure and evolution, which still present difficulties in simulating stars such as Canopus in detail. Key words. Stars, individual: Canopus – Stars: massive – Stars: supergiants – Methods: observational, numerical – Techniques: high angular resolution, interferometric, spectroscopic, photometric

1. Introduction One of the crucial quantities for constraining stellar funda- mental parameters is the angular diameter, / , of the star. The Canopus (α Car, HD 45348) is the brightest star of the south- first direct measure of the angular diameter of Canopus, used ern night sky, and, as such, it has been extensively studied by many authors, was made by Hanbury Brown et al. (1974a), and observed. Its reported spectral type turns around the yel- based on the intensity interferometry technique. Several other low supergiant class, for example A9 II (SIMBAD) and F0 Ib- direct measurements – notably from optical-IR long baseline in- II (Jerzykiewicz & Molenda-Zakowicz 2000). Canopus is thus a terferometry (OLBI) – and indirect measurements exist in the lit- arXiv:2109.07153v1 [astro-ph.SR] 15 Sep 2021 rare, massive (∼ 10M ; e.g., Jerzykiewicz & Molenda-Zakowicz erature, with reported angular diameters around ∼ 6−7 milliarc- 2000; Smiljanic et al. 2006) yellow supergiant that can be well seconds (mas). Cruzalèbes et al. (2013) provide a compilation observed from Earth, providing a unique opportunity to measure of many of these angular diameter estimates. It is important to its fundamental parameters precisely, which is a key step toward note that a precise estimation of / from high angular resolution understanding in detail the physical structure and evolution of techniques is intimately related to the estimation of the center- massive stars beyond the main sequence. to-limb intensity profile – limb darkening (LD) – at the spectral domain of the observation. These LD estimations also provide important constraints on the physical structure of the stellar at- ? Based on observations performed at ESO, Chile under program mosphere, which depend on several physical parameters (e.g., IDs 60.A-9237 and 092.D-0366 for PIONIER, and 084.D-0151 (VISA- mass, temperature, gravity, chemical composition). These obser- CNRS GTO) for AMBER. ?? e-mail: [email protected] vational constraints are particularly important for yellow super-

Article number, page 1 of 14 The NEVOD-EAS air-shower array

M.B. Amelchakov1, N.S. Barbashina1, A.G. Bogdanov1, A. Chiavassa1,2,3, D.M. Gromushkin1, S.S. Khokhlov1, V.V. Kindin1, R.P. Kokoulin1, K.G. Kompaniets1, A.Yu. Konovalova1, V.V. Ovchinnikov1, N.A. Pasyuk1, A.A. Petrukhin1, I.A. Shulzhenko1, V.V. Shutenko1, I.I. Yashin1, K.O. Yurin1

1 National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), 115409 Moscow, Russia 2 Dipartimento di Fisica dell’ Università degli Studi di Torino, 10125 Torino, Italy 3 Sezione di Torino dell’ Istituto Nazionale di Fisica Nucleare – INFN, 10125 Torino, Italy

Abstract The Experimental complex NEVOD includes several different setups for studying various components of extensive air showers (EAS) in the energy range from 1010 to 1018 eV. The NEVOD-EAS array for detection of the EAS electron-photon component began its data taking in 2018. It is a distributed system of scintillation detectors installed over an area of about 104 m2. A distinctive feature of this array is its cluster organization with different-altitude layout of the detecting elements. The main goal of the NEVOD-EAS array is to obtain an estimation of the primary particle energy for events measured by various detectors of the Experimental complex NEVOD. This paper describes the design, operation principles and data processing of the NEVOD-EAS array. The criteria for the event selection and the accuracy of the EAS parameters reconstruction obtained on the simulated events are discussed. The results of the preliminary analysis of experimental data obtained during a half-year operation are presented. Keywords: cosmic rays, primary particles, extensive air showers, scintillation detectors 1. Introduction The Experimental complex (ЕС) NEVOD is located at the campus of the National Research Nuclear University MEPhI (Moscow, Russia) and is designed to study various components of cosmic rays on the Earth's surface. The basis of the EC NEVOD is a 2-kiloton Cherenkov water detector with a spatial array of quasi-spherical optical modules [1] which is surrounded by a large-area coordinate-tracking detector DECOR [2]. Since 2002, an experiment for studying muon component of inclined extensive air showers (EAS) using the method of local muon density spectra (LMDS) has been carried out at the Experimental Complex NEVOD. As a result of this experiment, an excess of muons in bundles in comparison with the calculations for all existing models of hadronic interactions has been discovered [3]. This phenomenon was named the “muon puzzle” [4, 5]. To calibrate the LMDS method, a facility for detection of EAS electromagnetic component (air-shower array) has been constructed around the building of the Cherenkov water detector. To detect EAS particles in this array, the scintillation counters which previously operated in the EAS-TOP [6] and KASCADE-Grande [7] experiments are used. During the array constructing, the main difficulty was the high building density of the MEPhI campus. In this regard, at the development stage the array structure has been organized by a cluster approach. According to this approach, the array consists of independent clusters installed on the roofs of the MEPhI laboratory buildings and on the ground surface between them.

1

MNRAS 000,1–13 (2021) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

Creating Jackknife and Bootstrap estimates of the covariance matrix for the two-point correlation function

Faizan G. Mohammad,1,2★ Will J. Percival,1,2,3 1Waterloo Center for Astrophysics, University of Waterloo, Waterloo, ON N2L 3G1, Canada 2 Department of Physics and Astronomy, University of Waterloo, Waterloo, ON N2L 3G1, Canada 3 Perimeter Institute for Theoretical Physics, 31 Caroline St. North, Waterloo, ON N2L 2Y5, Canada

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We present correction terms that allow delete-one Jackknife and Bootstrap methods to be used to recover unbiased estimates of the data covariance matrix of the two-point correlation function 휉 (r). We demonstrate the accuracy and precision of this 3 new method using a large set of 1000 QUIJOTE simulations that each cover a comoving volume of 1h−1Gpc . The corrected resampling techniques accurately recover the correct amplitude and structure of the data covariance matrix as represented by its principal components. Our corrections for the internal resampling methods are shown to be robust against the intrinsic clustering of the cosmological tracers both in real- and redshift space using two snapshots at 푧 = 0 and 푧 = 1 that mimic two samples with significantly different clustering. We also analyse two different slicing of the simulation volume into 푛sv = 64 or 125 sub-samples and show that the main impact of different 푛sv is on the structure of the covariance matrix due to the limited number of independent internal realisations that can be made given a fixed 푛sv. Key words: keyword1 – keyword2 – keyword3

1 INTRODUCTION accurate simulations of the non-linear effects are computationally expensive. Fast approximate simulations are usually employed to The two-point correlation function of the galaxy distribution contains provide thousands of copies of real data. However they can fail to a wealth of information about the gravity-driven growth of cosmo- accurately predict the non-linear effects in the growth of structure, logical structure and the expansion rate of the Universe. As such the and their accuracy on small scales is debatable (e.g. Zhao et al. 2021). two-point correlation function and the power spectrum, its Fourier- The recovered matrix is drawn from a Wishart distribution, and this space counterpart, have become an essential statistics to measure contributes to the recovered model parameter errors (Hartlap et al. and fit with models of galaxy clustering expected in galaxy redshift 2007; Dodelson & Schneider 2013; Percival et al. 2014; Sellentin & surveys (e.g. Alam et al. 2017, 2021). Heavens 2016; Percival et al. 2021). For the analytical calculations Ongoing and planned experiments such as Euclid (Laureijs et al. of the covariance matrix required for (ii), the simplifications required 2011) and the Dark Energy Spectroscopic Instrument (DESI; DESI to be able to calculate the covariance may result in a biased estimate Collaboration et al. 2016b,a) will provide an unprecedented amount of the data covariance matrix, and the inclusion of window function of data enabling a percent level precision on the estimates of the and other effects is complicated (Howlett & Percival 2017; Blake key cosmological parameters. As with any experiment, the data anal- et al. 2018; Li et al. 2019). Internal estimators on the other hand ysis requires an equally accurate and precise estimate of the data use appropriate resampling of the observed data and as such are covariance matrix, a key ingredient of model fitting. In cosmological less prone to unknown physics. However, they are found to provide arXiv:2109.07071v1 [astro-ph.CO] 15 Sep 2021 surveys this is a challenging task as we are provided with only one estimates of the covariance matrix that can overestimate the true realisation of the observed data. covariance matrix by anywhere between 25-60% (Norberg et al. Three different techniques have been proposed in the past literature 2009; Friedrich et al. 2016; Favole et al. 2020). Schemes to provide to obtain estimates of the data covariance matrix: i) use of simulated better convergence are generally complicated, such as the marked datasets that provide an ensemble of synthetic replicas of the real point bootstrap method of Loh(2008). In this paper we focus on (iii), observations that mimic the underlying probability distribution of exploring the application of the Jackknife and Bootstrap methods to the observed data (e.g. Zhao et al. 2021); ii) Analytical covariances analyses of the correlation function. that rely on a number of assumptions and theoretical models of the galaxy clustering in the Universe (e.g. Blake et al. 2018); iii) Internal Complications due to the window function make internal methods estimators that resample directly the observed data using specific difficult to apply to the power spectrum, where subsampling would schemes to generate multiple copies of the observations (e.g. Norberg change the window effects (e.g. Wilson et al. 2017). This is not the et al. 2009). Each of these techniques has its pros and cons: For (i), case for the correlation function, where the estimator removes win- dow effects, although one still have to be careful to subsample in a way that preserves volume scaling. Even so, the application is not ★ E-mail: [email protected] straightforward because of the pairwise nature of the measurement

© 2021 The Authors Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX631

Inpainting hydrodynamical maps with deep learning

Faizan G. Mohammad ,1, 2 Francisco Villaescusa-Navarro ,3, 4 Shy Genel ,4, 5 Daniel Angles-Alc´ azar´ ,6, 4 and Mark Vogelsberger 7

1Waterloo Center for Astrophysics, University of Waterloo, Waterloo, ON N2L 3G1, Canada 2Department of Physics and Astronomy, University of Waterloo, Waterloo, ON N2L 3G1, Canada 3Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton NJ 08544-0010, USA 4Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, 10010, New York, NY, USA 5Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY, 10027, USA 6Department of Physics, University of Connecticut, 196 Auditorium Road, U-3046, Storrs, CT, 06269, USA 7Kavli Institute for Astrophysics and Space Research, Department of Physics, MIT, Cambridge, MA 02139, USA

ABSTRACT From 1,000 hydrodynamic simulations of the CAMELS project, each with a different value of the cosmological and astrophysical parameters, we generate 15,000 gas temperature maps. We use a state- of-the-art deep convolutional neural network to recover missing data from those maps. We mimic the missing data by applying regular and irregular binary masks that cover either 15% or 30% of the area of each map. We quantify the reliability of our results using two summary statistics: 1) the distance between the probability density functions (pdf), estimated using the Kolmogorov-Smirnov (KS) test, and 2) the 2D power spectrum. We find an excellent agreement between the model prediction and the unmasked maps when using the power spectrum: better than 1% for k < 20 h/Mpc for any irregular mask. For regular masks, we observe a systematic offset of 5% when covering 15% of the maps while the results become unreliable when 30% of the data is missing.∼ The observed KS-test p-values favor the null hypothesis that the reconstructed and the ground-truth maps are drawn from the same underlying distribution when irregular masks are used. For regular-shaped masks on the other hand, we find a strong evidence that the two distributions do not match each other. Finally, we use the model, trained on gas temperature maps, to perform inpainting on maps from completely different fields such as gas mass, gas pressure, and electron density and also for gas temperature maps from simulations run with other codes. We find that visually, our model is able to reconstruct the missing pixels from the maps of those fields with great accuracy, although its performance using summary statistics depends strongly on the considered field.

Keywords: large-scale structure of universe – methods: numerical – methods: statistical

1. INTRODUCTION open questions like the sum of the neutrino masses, the Cosmology is in a transformative stage. Nowadays, community has spent billions of dollars in surveys like we know the value of the main cosmological parameters the Dark Energy Spectroscopic Instrument (DESI, DESI arXiv:2109.07070v1 [astro-ph.CO] 15 Sep 2021 with a relatively high precision. This has allowed us Collaboration et al. 2016), Euclid (Laureijs et al. 2011), to claim, with high confidence, the existence of a sub- Prime Focus Spectrograph (PFS, Tamura et al. 2016), stance that is responsible for the accelerated expansion extended ROentgen Survey with an Imaging Telescope of the Universe: dark energy. The nature and proper- Array (eROSITA, Merloni et al. 2012), Roman Obser- ties of dark energy remain the biggest mysteries in mod- vatory (Spergel et al. 2015), Rubin Observatory (The ern physics. In order to shed light on these and other LSST Dark Energy Science Collaboration et al. 2018), Square Kilometer Array (SKA, Square Kilometre Array Cosmology Science Working Group et al. 2018), and Si- Corresponding author: Faizan G. Mohammad mons Observatory (Ade et al. 2019) whose data may [email protected] contain the answer to all these fundamental questions. arXiv:2109.07068v1 [astro-ph.EP] 15 Sep 2021 uigteerypaeo lntr irlnigexperi- microlensing ( planetary OGLE of example, phase for early ments, the During Introduction 1. h xsec fpae-nue nmle ntelgtcurv light the in anomalies severa planet-induced was of events existence to lensing inspected the thoroughly of be could rate events individual detection dozens, annual the which ( etme 6 2021 16, September Astronomy loke l 1993 al. et Alcock ree Mróz Przemek yu-o Kim Hyoun-Woo ihe .Albrow D. Michael eevd;accepted ; Received assaei h ag f0 of range the in are masses e words. Key high late-time of role the sources. discuss a We stars. in p low-mass planets microlensing The around other Jupiter. many with is, together that systems, system, planetary solar the of planet asv hntesn nteohrhn,tepaes ihmas with planets, the hand, other the On sun. the than massive Results. Methods. Aims. ffidn lnt htwr isddet h eitoso pl of deviations the to due missed were that planets finding of rsig ffitsuc tr vrtersnn asisf caustics com resonant a the over have stars events source the faint of of curves crossings light lensing The BLG-0299Lb. 15 14 11 13 10 12 nG Shin In-Gu 6 5 4 3 9 1 2 8 7 zmnKozłowski Szymon nvriyo atruy eateto hsc n Astron and Physics of Department Canterbury, of University siuoNzoaed iiaNcer,Szoed aoi V Napoli, di Sezione Nucleare, Fisica di Nazionale Istituto eateto hsc,Uiest fWrik ibtHl R Hill Gibbet Warwick, of University Physics, of Department iiino hsc,Mteais n srnm,Californ Astronomy, and Mathematics, Physics, of Division iatmnod iia“.R aailo,Uiest iSa di Università Caianiello”, R. “E. Fisica di Dipartimento etrfrAstrophysics for Center colo pc eerh yn e nvriy ogn Kye Yongin, University, Hee Kyung Research, Space of School eateto atcePyisadAtohsc,Weizmann Astrophysics, and Physics Particle of Department oe srnm n pc cec nttt,Deo 34055, Daejon Institute, Science Space and Astronomy Korea eateto srnm,Tigu nvriy ejn 100 Beijing University, Tsinghua Astronomy, of Department eateto srnm,TeOi tt nvriy 4 .1 W. 140 University, State Ohio The Astronomy, of Department eateto hsc,Cugu ainlUiest,Cheo University, e-mail: National Chungbuk Physics, of Department srnmclOsraoy nvriyo asw l Ujazd Al. Warsaw, of University Observatory, Astronomical a lnkIsiuefrAtooy öisul1,D-69117 17, Königstuhl Astronomy, for Institute Planck Max oe nvriyo cec n ehooy 1 Gajeong-ro 217 Technology, and Science of University Korea ecnutdapoeto enetgtn h 0721 mic 2017–2019 the reinvestigating of project a conducted We & o l lntr vns h esn ouin r uniquel are solutions lensing the events, planetary all For he an-oremcoesn lnt eetdvia detected planets microlensing faint-source Three h rjc e st n he lnt nldn KMT-2017-B including planets three find to us led project The Astrophysics [email protected] 2 rvttoa irlnig–paesadstlie:dete satellites: and planets – microlensing gravitational , 3 hogoHan Cheongho 14 ,EO ( EROS ), os Shvartzvald Yossi , 3 ihłK Szyma K. Michał , en-e Kim Seung-Lee , 6 u-uChung Sun-Ju , aucitn.ms no. manuscript 2 uor ta.1993 al. et Aubourg ryzo Ulaczyk Krzysztof , | . avr mtsna 0Gre t,Cmrde A018 US 02138, MA Cambridge, St., Garden 60 Smithsonian & Harvard dlk ta.1994 al. et Udalski 45 1 . nre Udalski Andrzej , M / 10 3 eoatcutcchannel resonant-caustic M , nski ´ 7 enfrC Yee C. Jennifer , ⊙ ogJoLee Dong-Joo , 3 . , 2 7 0 nrwGould Andrew , a Skowron Jan , . 9 hc orsod oeryMt aeKdaf,adtu the thus and dwarfs, K late to M early to corresponds which 59, TeKTe Collaboration), KMTNet (The uvy,for surveys, ) TeOL Collaboration) OGLE (The ,MACHO ), 15 ais Gromadzki Mariusz ryzo .Rybicki A. Krzysztof , 2 Laigauthors) (Leading oo Kim Doeon , re ygatpaeslctdna h isenrnso hos of rings Einstein the near located planets giant by ormed check lsseslebyn h nwlnso h ot,adtu the thus and hosts, the of lines snow the beyond lie systems ll aeaysses upr htmsiegsgatpaesa planets gas-giant massive that support systems, lanetary ABSTRACT 3 rslto mgn ncaiyn eoatiaelne w lenses resonant-image clarifying in imaging -resolution 11 ntr inl rmtetpclfr fsottr anomali short-term of form typical the from signals anetary en,VaGoan al I -48 icao(A,Italy (SA), Fisciano I-84084 II, Paolo Giovanni Via lerno, aCni,I816Npl,Italy Napoli, I-80126 Cintia, ia ogekLee Yongseok , es aIsiueo ehooy aaea A915 USA 91125, CA Pasadena, Technology, of Institute ia wke4 048Wrzw,Poland Warszawa, 00-478 4, owskie 2 e nterneo 2 of range the in ses ecegZang Weicheng , m,PiaeBg40,Crscuc 00 e Zealand New 8020, Christchurch 4800, Bag Private omy, l 8,China 084, gu264 eulco Korea of Republic 28644, ngju nttt fSine eoo 60,Israel 76100, Rehovot Science, of Institute aołwPoleski Radosław , a,Cvnr,C47L UK 7AL, CV4 Coventry, oad, uen-u ajo,313 eulco Korea of Republic 34113, Daejeon, Yuseong-gu, , o hrceitcta h lntr inl eeproduce were signals planetary the that characteristic mon edleg Germany Heidelberg, eulco Korea of Republic 8 eemndwtotaydgnrc.I setmtdta th that estimated is It degeneracy. any without determined y ng 70,Rpbi fKorea of Republic 17104, onggi , t v. oubs H420 USA 43210, OH Columbus, Ave., 8th oesn aacletdb h ihcdnesreswt t with surveys high-cadence the by collected data rolensing inl o oeeet a saedetection. plan escape event, may each events for some data even for of signals of quantity number large increased the greatly with the together With year. per 3000 than ( KMTNet ftelnigeet.Wt h deto ihcdnesurve high-cadence of ( advent MOA the is, With that events. lensing the of 9 ction y-aHwang Kyu-Ha , 1 G20L,OL-07BG19L,adOGLE-2019- and OGLE-2017-BLG-1099Lb, LG-2509Lb, onHu Ryu Yoon-Hyun , 2 i ta.2016 al. et Kim 3 2 , 13 arkIwanek Patryk , uie l 2013 al. et Sumi 12 . 1 yogGnPark Byeong-Gon , agMkCha Sang-Mok , . 2 grSoszy Igor , M 3 onKlJung Kil Youn , / M J 3 A aei Bozza Valerio , ,tedtcinrt a ordt more to soared has rate detection the ), . 6 . ,aehairta h heaviest the than heavier are 2, ,OL-V( OGLE-IV ), 2 acnWrona Marcin , nski ´ 3 , 13 2 3 , ril ubr ae1o 10 of 1 page number, Article ae Pietrukowicz Paweł , 3 7 ogJnKim Dong-Jin , hn-kLee Chung-Uk , ihr .Pogge W. Richard , otsasaeless are stars host 4 ecommonplace re , dlk ta.2015 al. et Udalski 5 stars. t t eyfaint very ith es. discovered 2 and , © ythe by d eaim he S 2021 ESO host e 3 3 , , etary 2 9 ys, , ts ), The relativistic Pythagorean three-body problem The relativistic Pythagorean three-body problem Tjarda C. N. Boekholta) Rudolf Peierls Centre for Theoretical Physics, Clarendon Laboratory, Parks Road, OX1 3PU, Oxford, UK. Arend Moerman and Simon F. Portegies Zwartb) Leiden Observatory, Leiden University, PO Box 9513, 2300 RA, Leiden, The Netherlands. (Dated: 16 September 2021) We study the influence of relativity on the chaotic properties and dynamical outcomes of an unstable triple system; the Pythagorean three-body problem. To this end, we extend the Brutus N-body code to include Post-Newtonian pairwise terms up to 2.5 order, and the first order Taylor expansion to the Einstein-Infeld-Hoffmann equations of motion. The degree to which our system is relativistic depends on the scaling of the total mass (the unit size was 1 parsec). Using the Brutus method of convergence, we test for time-reversibility in the conservative regime, and demonstrate that we are able to obtain definitive solutions to the relativistic three-body problem. It is also confirmed that the minimal required numerical accuracy for a successful time-reversibility test correlates with the amplification factor of an initial perturbation, as was found previously for the Newtonian case. When we take into account dissipative effects through gravitational wave emission, we find that the duration of the resonance, and the amount of exponential growth of small perturbations depend on the mass scaling. For a unit mass ≤ 10 M , the system behavior is indistinguishable from Newton’s equations of motion, and the resonance always ends in a 7 binary and one escaping body. For a mass scaling up to 10 M , relativity gradually becomes more prominent, but the majority of the systems still dissolve in a single body and an isolated 5 7 9 binary. The first mergers start to appear for a mass of ∼ 10 M , and between 10 M and 10 M all systems end prematurely in a merger. These mergers are preceded by a gravitational 9 wave driven in-spiral. For a mass scaling ≥ 10 M , all systems result in a gravitational wave merger upon the first close encounter. Relativistic three-body encounters thus provide an efficient pathway for resolving the final parsec problem. The onset of mergers at the 7 characteristic mass scale of 10 M potentially leaves an imprint in the mass function of supermassive black holes.

Keywords: Chaos – Relativity – N-body simulations – Black holes

I. INTRODUCTION subsequent evolution, the bodies fall towards the center

1 of mass of the triple, after which the prolonged, chaotic The Newtonian three-body problem is one of the stan- interaction takes place. After a first close encounter be- dard examples for illustrating chaos. This line of research 2 tween the three bodies, they experience multiple close dates back to Poincar´e , but more recent demonstrations encounters until the least massive body escapes, leaving have been made using computers. For instance, the se- 3 the other two in a stable binary orbit. For a detailed ries of papers started by Hut and Bahcall focuses on discussion on the orbital evolution, we refer the reader to binary-single scattering experiments, in which a single Szebehely and Peters 5 .

arXiv:2109.07013v1 [astro-ph.IM] 14 Sep 2021 body approaches a binary system with a certain impact parameter and relative velocity at infinity. They mea- The Pythagorean problem exhibits exponential sensi- sure the cross section for the occurrence of a democratic tivity to small changes in the initial condition6. Over resonance, e.g. a prolonged interaction between all three the lifetime of the interaction, small perturbations grow bodies with similar pair-wise gravitational forces. An by about 9 orders of magnitude7,8, equivalent to about example of the orbital chaos during such a resonance is 20 e-folding time scales. This is one of the reasons why given in their Fig. 3. the Pythagorean three-body problem is excellently suit- A textbook example for a chaotic triple system is the able for studying chaotic N-body systems numerically; Pythagorean problem, also called Burrau’s problem4. In the time scale on which perturbations grow is short com- this case, the orbital chaos is not initiated by a binary- pared to the system’s lifetime. As a result, reprehensible8 single encounter, but by three bodies performing a cold N-body codes have considerable difficulty in reaching collapse, e.g. three bodies with masses 3, 4, and 5 are a converged solution to such chaotic N-body problems, positioned at rest in a planar, right-angled triangle at whereas for an arbitrary precise N-body code this can be positions (1, 3), (-2, -1), and (1, -1), respectively. In the achieved rather effortlessly9. Note that within this con- ,

Testing Self-Organized Criticality Across the Main Sequence using Stellar Flares from TESS

Adina D. Feinstein∗ Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA

Darryl Z. Seligman Department of the Geophysical Sciences, University of Chicago, Chicago, IL 60637, USA

Maximilian N. G¨unther Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology (MIT), Cambridge, MA 02139, USA and European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands†

Fred C. Adams Physics Department, University of Michigan, Ann Arbor, MI 48109 and Astronomy Department, University of Michigan, Ann Arbor, MI 48109 (Dated: September 16, 2021) Stars produce explosive flares, which are believed to be powered by the release of energy stored in coronal magnetic field configurations. It has been shown that solar flares exhibit energy distributions typical of self-organized critical systems. This study applies a novel flare detection technique to data obtained by NASA’s TESS mission and identifies ∼ 106 flaring events on ∼ 105 stars across spectral types. Our results suggest that magnetic reconnection events that maintain the topology of the magnetic field in a self-organized critical state are ubiquitous among stellar coronae.

Introduction — The theory of self-organized critical- of earthquake energies [16–18], ity [1] describes a class of dissipative dynamical systems dN which remain at a critical point with no intrinsic length ∼ E−α , (1) or time scale. The existence of the critical state re- dE quires a local instability, which occurs when some pa- where N is the number of earthquakes, E is the en- rameter exceeds its critical value and results in a dissi- ergy released in the earthquake (where the earthquake pative transport process where this same parameter in- magnitude m ∝ log E), and the power-law exponent creases at neighboring sites. A simple physical analogy 1.25 < α ∼ 1.5. A 3-dimensional slip-stick model of is a sand pile. As sand particles are added, the difference tectonic events produces a critical exponent of α = 1.35 in height between neighboring sites on the pile increases. [16], in reasonable agreement with the observed law. When the additional sand particles make the new height Even the solar coronal magnetic field may reside in difference exceed a critical threshold, avalanche events a self-organized critical state [19]. This hypothesis nat- occur. This system maintains a critical slope, represent- urally explains the power-law dependence of the magni- ing a dynamical attractor that is insensitive to the initial tude of solar flares, which is of the same form as Equation conditions. This critical state is maintained via nonlin- (1), where α ≈ 1.4 [19]. This characteristic exponent and ear avalanche events spanning all length scales triggered the subsequent temporal clustering are universal between by perturbations. earthquakes and solar flares [20]. arXiv:2109.07011v1 [astro-ph.SR] 14 Sep 2021 Explosive events on the Sun are believed to be driven by the energy stored in twisted coronal magnetic field While the sand pile analogy is simplistic by construc- lines [21]. Such field configurations are generated through tion, self-organized criticality naturally manifests in a va- dynamo action in the fluid-dominated interior [22], and riety of physical systems. Applications have been found through convective and coriolis driven vortical subsurface in hydrodynamical turbulence, forest fires and other per- plasma flows [23–26]. The concentration of magnetic field colation systems [2], landslides [3, 4], neuroscience [5, 6], lines can lead to the release and subsequent dissipation climate fluctuations [7], rainfall [8], accretion disks [9], of energy via the process of magnetic reconnection [27]. traffic jams [10], evolution [11], extinction [12], finan- These reconnection events could be triggered when the cial markets [13] and even Conway’s game of Life [14], angle, θ, between neighboring magnetic field vectors is to name a non-comprehensive list. The theory also ex- greater than a critical value, θc. The reconnection event plains the Gutenberg-Richter [15] law for the distribution changes the angles for neighboring field lines and can

universe

Article Search for and Identification of Young Compact Galactic Supernova Remnants Using THOR

Sujith Ranasinghe * , Denis Leahy and Jeroen Stil

Department of Physics and Astronomy, University of Calgary, Calgary, AB T2N 1N4, Canada; [email protected] (D.L.); [email protected] (J.S.) * Correspondence: [email protected]

Abstract: Young Supernova remnants (SNRs) with smaller angular sizes are likely missing from existing radio SNR catalogues, caused by observational constraints and selection effects. In order to find new compact radio SNR candidates, we searched the high angular resolution (2500) THOR radio survey of the first quadrant of the galaxy. We selected sources with non-thermal radio spectra. HI absorption spectra and channel maps were used to identify which sources are galactic and to estimate their distances. Two new compact SNRs were found: G31.299-0.493 and G18.760-0.072, of which the latter was a previously suggested SNR candidate. The distances to these SNRs are 5.0 0.3 kpc and 4.7 0.2 kpc, respectively. Based on the SN rate in the galaxy or on the statistics of ± ± known SNRs, we estimate that there are 15–20 not-yet detected compact SNRs in the galaxy and that the THOR survey area should contain three or four. Our detection of two SNRs (half the expected number) is consistent with the THOR sensitivity limit compared with the distribution of integrated flux densities of SNRs.

ꢀꢀꢀ ꢂꢀꢃꢄꢅꢆꢇ Keywords: supernova remnants; radio continuum; radio lines ꢀ ꢀ ꢂꢃꢄꢅꢆ

Citation: Ranasinghe, S.; Leahy, D.; Stil, J. Search for and Identification of Young Compact Galactic Supernova 1. Introduction Remnants Using THOR. Universe For a better understanding of the interstellar medium (ISM) and the evolution of our 2021, 7, 338. https://doi.org/ galaxy, we need not only to study the properties of individual supernova remnants (SNRs) 10.3390/universe7090338 but also to obtain an as complete as possible sample of SNRs. The current observed number of galactic SNRs is 294 [1]. There is an apparent deficit between the number of observed Academic Editor: Jacco Th. van Loon and expected galactic radio SNRs [2,3]. 00 Received: 29 July 2021 From the [1] catalogue, the SNR with the smallest angular radius has a radius 45 . Accepted: 2 September 2021 SNR G1.9+0.3 is one of the smallest galactic SNRs, and [4] finds that it is near the galactic center. With the distance to the galactic center R 0 = 8.34 0.16 kpc [5], the radius is 1.8 pc. Published: 9 September 2021 ± Many other small SNRs with known distances have radii of 2 pc. ∼ Publisher’s Note: MDPI stays neutral Assuming a mean lifetime of 60,000 yr for radio SNRs [6] and a supernova (SN) rate ∼ with regard to jurisdictional claims in of one per 40 10 yr [7], the expected number of SNRs is 1500. The reasons for the lack published maps and institutional affil- of observed SNRs± is likely the result of selection effects acting∼ against the discovery of old, iations. faint, large remnants as well as very young, small remnants due to poor sensitivity and spa- tial resolution [3,8]. The identification of faint, old, and large remnants remains a difficult task because of their low surface brightness and the complex structure of background radio emission. However, recent observations from the low-frequency GLEAM survey have Copyright: © 2021 by the authors. produced a number of SNRs and candidates [9,10]. The majority of these SNRs and SNR Licensee MDPI, Basel, Switzerland. candidates have radii ranging from 20 to 500. The SNR candidates presented by [11] ∼ ∼ This article is an open access article using the MAGPIS data have radii of a few arcmin. In this work, we searched for compact distributed under the terms and SNR candidates (radii < few arcmin). With a 2500 angular resolution [12] and sensitive up conditions of the Creative Commons to 12000 [13], the THOR survey should be ideal for identifying compact remnants. Attribution (CC BY) license (https:// For catalogued SNRs with distance estimates (Green’s Galactic SNR catalogue [1]), creativecommons.org/licenses/by/ 9 have a radius <3 pc and 20 have radius <5 pc. These correspond to angular radii 4.0/). ∼ ∼ of 3900 and 6400 at 16 kpc, respectively, or 7500 and 12300 at a more common distance of

Universe 2021, 7, 338. https://doi.org/10.3390/universe7090338 https://www.mdpi.com/journal/universe Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX631

Self-Similarity of k-Nearest Neighbor Distributions in Scale-Free Simulations

Lehman H. Garrison ,1 Tom Abel ,2, 3, 4 and Daniel J. Eisenstein5

1Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, New York, NY 10010, USA 2Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA 3Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA 4SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 5Center for Astrophysics | Harvard & Smithsonian 60 Garden Street, Cambridge, MA 02138, USA

ABSTRACT We use the k-nearest neighbor probability distribution function (kNN-PDF, Banerjee & Abel 2021) to assess convergence in a scale-free N-body simulation. Compared to our previous two-point analysis, the kNN-PDF allows us to quantify our results in the language of halos and numbers of particles, while also incorporating non-Gaussian information. We find good convergence for 32 particles and greater at densities typical of halos, while 16 particles and fewer appears unconverged. Halving the softening length extends convergence to higher densities, but not to fewer particles. Our analysis is less sensitive to voids, but we analyze a limited range of underdensities and find evidence for convergence at 16 particles and greater even in sparse voids.

1. INTRODUCTION have a well-defined continuum value—consider that tak- Cosmological N-body simulations model the evolution ing  → 0 will increase scattering encounters between of the continuous Vlasov-Poisson distribution function particles, which does not occur in the continuous Vlasov- with a set of N particles interacting under mutual grav- Poisson system. itational attraction. The discrete nature of the particles Another class of tests exists that exploits the self- imposes a small-scale cutoff to the resolution of the sim- similar nature of gravity, known as scale-free simula- ulation, although precisely the length and mass scale tions. With a power-law power spectrum of index ns of such a cutoff and the mechanism by which it oper- and an ΩM = 1 background cosmology, only one scale is ates is not so clear. In some limits, such as the weakly present in the continuum version of the problem, with perturbative regime, it is possible to show that the cor- time and space related by a power law whose index is respondence of the dynamics to the continuum solution determined by ns. In other words, a small-scale prop- quickly degrades as one approaches the mean interpar- erty of the simulation at one time must be equal to the ticle separation, ` = N 1/3/L (for box size L, Garrison same property measured on large-scales at a later time. et al. 2016). However, in the strongly non-linear regime, Where such self-similarity is not observed, the simula- resolution is often obtained at scales many times smaller tion may be said to differ from the continuum solution. than the mean particle separation. The relevant cut- Scale-free simulations have a rich history in the N- off is therefore not `, but something less obvious, and body literature (e.g. Efstathiou et al. 1988; Colombi the question becomes the scale of such a cutoff and the et al. 1996; Jain & Bertschinger 1998; Scoccimarro & mechanism by which it manifests. And to complicate Frieman 1999; Smith et al. 2003; Widrow et al. 2009; Or- arXiv:2109.06991v1 [astro-ph.CO] 14 Sep 2021 matters, additional parameters, such as softening, are ban & Weinberg 2011), alternately being used to build often introduced to regularize the particle interactions. ΛCDM predictions and as tests of N-body dynamics. Their role in setting the resolution limit of a simulation Scale-free tests do not establish the absolute convergence must also be assessed. of a simulation, but do guarantee that any remaining Traditional convergence tests operate by tuning a dis- errors must themselves be of a self-similar nature. For creteness parameter, such as N, towards its continuum many classes of error, this is sufficient, as UV cutoffs value, in this case infinity. But such tests are limited like N and , as well as IR cutoffs like L, impose a pre- in their dynamic range, as the computational expense ferred scale. A notable exception is certain classes of necessarily increases towards the continuum. Further- time stepping; for example, a fixed time step in log-a more, parameters such as the softening length  do not will yield self-similar results for any time step size, even implausibly large. Even time step schemes that are not Mon. Not. R. Astron. Soc. 000, 000–000 (0000) Printed 16 September 2021 (MN LATEX style file v2.2)

A Speckle-Imaging Search for Close and Very Faint Companions to the Nearest and Brightest Wolf-Rayet Stars

Michael M. Shara1?, Steve B. Howell2, Elise Furlan3, Crystal L. Gnilka2, Anthony F.J. Moffat 4, Nicholas J. Scott2, and David Zurek1 1Department of Astrophysics, American Museum of Natural History, Central Park West at 79th Street, New York, NY 10024, USA 2NASA Ames Research Center, Moffett Field, CA 94035, USA 3NASA Exoplanet Science Institute, Caltech/IPAC, Mail Code 100-22, 1200 E. California Blvd., Pasadena, CA 91125, USA 4D´epartement de Physique et Centre de Recherche en Astrophysique du Qu´ebec, Universit´ede Montr´eal,Montr´eal,QC H3C 3J7, Canada

Accepted Received

ABSTRACT Gravitationally bound companions to stars enable determinations of their masses, and offer clues to their formation, evolution and dynamical histories. So motivated, we have carried out a speckle imaging survey of eight of the nearest and brightest Wolf-Rayet (WR) stars to directly measure the frequency of their resolvable companions, and to search for much fainter companions than hitherto possible. We found one new, close companion to each of WR 113, WR 115 and WR 120 in the separation range ∼ 0.2” - 1.2”. Our results provide more evidence that similar-brightness, close companions to WR stars are common. More remarkably, they also demonstrate that the predicted, but much fainter and thus elusive companions to WR stars are now within reach of modern speckle cameras on 8m class telescopes by finding the first example. The new companion to WR 113 is just 1.16” distant from it, and is ∼ 8 magnitudes fainter than the WR star. The empirical probability of a chance line-of-sight of the faint companion at the position of WR 113 is < 0.5%, though we cannot yet prove or disprove if the two stars are gravitationally bound. If these three new detections are physical companions we suggest, based on their narrowband magnitudes, colours, reddenings and GAIA distances that the companions to WR113, WR 115 and WR 120 are an F-type dwarf, an early B-type dwarf, and a WNE-type WR star, respectively. Key words: binaries:close – stars:Wolf-Rayet – techniques:high angular resolution

1 INTRODUCTION AND MOTIVATION Woosley, Sukhbold, & Kasen 2021). WR stars are thus im- portant to both stellar astrophysics and supernova theory. 1.1 Wolf-Rayet stars The descendants of O stars are Wolf-Rayet (WR) stars,

arXiv:2109.06975v1 [astro-ph.SR] 14 Sep 2021 whose large initial masses (typically 20 - 100 M at birth), short lifetimes (typically ∼ 300 kyr) and strong emission lines make them excellent tracers of recent star forma- The intense and inhomogeneous stellar winds of WR −5 tion (Crowther 2007). Single Nitrogen-rich (WN type) and stars (>10 M /yr (Barlow, Smith, & Willis 1981; Carbon-rich (WC and even hotter but rare WO types) WR Schmutz, Hamann, & Wessolowski 1989; Hillier 1991; Nugis, stars likely die as type Ib (Gaskell et al. 1986) and Ic (Corsi Crowther, & Willis 1998; Hillier & Miller 1998; Vink & et al. 2012) supernovae, respectively. A large fraction of mas- de Koter 2005; Puls, Vink, & Najarro 2008; Smith 2014; sive stars are born in binaries that eventually exchange mass Sander et al. 2018; Poniatowski et al. 2021)) add a signif- (see references below), so that a major fraction of type Ib icant amount of C and O, and some N, to the interstellar and Ic SNe may be produced by the low mass, stripped cores medium (ISM), and contribute a significant fraction of the of formerly massive stars (Nomoto, Iwamoto, & Suzuki 1995; Galaxy’s ISM energy and momentum budgets. WR stars dis- Woosley, Langer, & Weaver 1995; Tauris, Langer, & Pod- play extremely strong emission lines (up to 1000 A˚ Equiva- siadlowski 2015; Prentice et al. 2019; Dessart et al. 2020; lent Width), which are probes of the rapidly expanding (∼ 3000 km/s (Underhill 1962; Abbott 1978)) and blob-filled (Moffat et al. 1988) atmospheres of these very hot, evolved ? E-mail: [email protected] stars.

© 0000 RAS Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX62

Optical Observations and Modeling of the Superluminous Supernova 2018lfe

Yao Yin,1 Sebastian Gomez,2, 1 Edo Berger,1 Griffin Hosseinzadeh,3, 1 Matt Nicholl,4 and Peter K. Blanchard5

1Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 2Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 3Steward Observatory, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA 4Birmingham Institute for Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 5Center for Interdisciplinary Exploration and Research in Astrophysics and Department of Physics and Astronomy, Northwestern University, 1800 Sherman Ave, Evanston, IL 60201, USA

ABSTRACT We present optical imaging and spectroscopy of SN 2018lfe, which we classify as a Type I superlumi- nous supernova (SLSN-I) at a redshift of z = 0.3501 with a peak absolute magnitude of Mr ≈ −22.1 mag, one of the brightest SLSNe discovered. SN 2018lfe was identified for follow-up using our FLEET machine learning pipeline. Both the light curve and the spectra of SN 2018lfe are consistent with the broad population of SLSNe. We fit the light curve with a magnetar central engine model and find an ejecta mass of Mej ≈ 3.8 M , a magnetar spin period of P ≈ 2.9 ms and a magnetic field strength of 14 B⊥ ≈ 2.8×10 G. The magnetic field strength is near the top of the distribution for SLSNe, while the spin period and ejecta mass are near the median values of the distribution for SLSNe. From late-time imaging and spectroscopy we find that the host galaxy of SN 2018lfe has an absolute magnitude of ∗ Mr ≈ −17.85, (LB ≈ 0.029 L ), and an inferred metallicity of Z ≈ 0.3 Z , and star formation rate of −1 ≈ 0.8 M yr .

Keywords: supernova: general – supernova: individual (SN 2018lfe)

1. INTRODUCTION Several possible mechanisms that can power SLSNe-I Superluminous supernovae (SLSNe) are a rare class of have been explored in the literature: large radioactive 56 core-collapse supernovae (SNe) that can exceed the lu- Ni mass (Pastorello et al. 2010), circumstellar interac- minosities of normal SNe by two orders of magnitude tion (Chevalier & Irwin 2011), and a magnetar central (Chomiuk et al. 2011; Quimby et al. 2011). Type I engine (Kasen & Bildsten 2010; Woosley 2010). The SLSNe (hereafter, SLSNe-I) are classified based on their bulk of the observational evidence points to a magne- hydrogen-free spectra, strong O II absorption lines at tar engine as the dominant mechanism. This includes early times, and blue continuum (Chomiuk et al. 2011; the diverse light curve behavior (Nicholl et al. 2017), Quimby et al. 2011; Nicholl et al. 2014a; Gal-Yam 2019; the early-time UV spectra (Mazzali et al. 2016; Nicholl Inserra et al. 2013). Despite their increasing discovery et al. 2017; Howell et al. 2013; Dessart 2019; Yan et al. rate over the past decade thanks to systematic wide-field 2017), the late-time light curve flattening (Nicholl et al. arXiv:2109.06970v1 [astro-ph.HE] 14 Sep 2021 optical surveys, SLSNe-I still remain a small population, 2018; Blanchard et al. 2021), and their nebular spec- with only about 150 spectroscopically classified events tra (Dessart et al. 2012; Nicholl et al. 2019; Jerkstrand to date (Gomez et al. 2020); moreover, not all of these 2017). The preference of SLSNe-I to low metallicity host events have well observed light curves and spectroscopic galaxies (similar to the hosts of the engine-powered long time series. gamma-ray bursts) also supports a magnetar engine en- ergy source (Lunnan et al. 2014; Leloudas et al. 2015; Ørum et al. 2020; Angus et al. 2019; Perley et al. 2016; Corresponding author: Yao Yin Chen et al. 2017; Schulze et al. 2018). [email protected] Here, we present detailed optical observations of SN 2018lfe, a transient first detected by the Zwicky Transient Facility (ZTF; Bellm et al. 2019) on 2018 Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus- like Atmosphere

Yuni Lee1,2,3,4, Chuanfei Dong5,6, and Valeriy Tenishev7

Abstract

Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO2-dominated atmosphere of a Venus-like rocky exoplanet, where dissociative

+ recombination of O2 ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D

Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. The extended hot O corona may lead to a more extended H exosphere through

1 Corresponding author [email protected]

2 University of Maryland Baltimore County, Baltimore, MD, USA 3 Center for Research and Exploration in Space Science and Technology II, Baltimore, MD, USA 4 NASA Goddard Space Flight Center, Planetary Environments Laboratory, code 699, Greenbelt, MD, USA 5 Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ, USA 6 Department of Astrophysical Sciences, Princeton University, Princeton, NJ, USA 7 Department of Climate and Space Sciences and Engineering, University of Michigan, Ann Arbor, MI, USA

1 manuscript submitted to Space Weather

OSPREI: A Coupled Approach to Modeling CME-Driven Space Weather with Automatically-Generated, User-Friendly Outputs

C. Kay1,2, M. L. Mays1, Y. M. Collado-Vega1

1Heliophysics Science Division, NASA Goddard Space Flight Center, Greenbelt, MD, USA 2Dept. of Physics, The Catholic University of America, Washington DC, USA

Key Points:

• OSPREI performs ensemble simulations of the Sun-to-Earth behavior of CMEs on time scales relevant for future forecasting

• A standardized set of automatically-generated visualizations can provide essen- tial information for forecasters in an intuitive manner

• Mimicking a forecasting approach, we apply OSPREI to the 22 April 2021 CME to illustrate its capabilities arXiv:2109.06960v1 [astro-ph.SR] 14 Sep 2021

Corresponding author: Christina Kay, [email protected]

–1–

manuscript submitted to Space Weather

Abstract Coronal Mass Ejections (CMEs) drive space weather activity at Earth and throughout the solar system. Current CME-related space weather predictions rely on information reconstructed from coronagraphs, sometimes from only a single viewpoint, to drive a sim- ple interplanetary propagation model, which only gives the arrival time or limited ad- ditional information. We present the coupling of three established models into OSPREI (Open Solar Physics Rapid Ensemble Information), a new tool that describes Sun-to- Earth CME behavior, including the location, orientation, size, shape, speed, arrival time, and internal thermal and magnetic properties, on the timescale needed for forecasts. First, ForeCAT describes the trajectory that a CME takes through the solar corona. Second, ANTEATR simulates the propagation, including expansion and deformation, of a CME in interplanetary space and determines the evolution of internal properties via conser- vation laws. Finally, FIDO produces in situ profiles for a CME’s interaction with a syn- thetic spacecraft. OSPREI includes ensemble modeling by varying each input param- eter to probe any uncertainty in their values, yielding probabilities for all outputs. Stan- dardized visualizations are automatically generated, providing easily-accessible, essen- tial information for space weather forecasting. We show OSPREI results for a CME ob- served in the corona on 2021 April 22 and at Earth on 2021 April 25. We approach this CME as a forecasting proof-of-concept, using information analogous to what would be available in real time rather than fine-tuning input parameters to achieve a best fit for a detailed scientific study. The OSPREI “prediction” shows good agreement with the arrival time and in situ properties.

Plain Language Summary Coronal mass ejections (CMEs) are large plasma structures that erupt from the Sun. These eruptions occur frequently but, depending on their exact properties (such as speed and magnetic field strength and orientation), can have negative effects on human tech- nologies both in space or on the surface of the Earth. It is important to know where a CME is going, when it will arrive, and its properties upon impact. Currently, most CME forecasts are limited to the time of arrival, and potentially the speed at that time, but have little to no additional information. We present OSPREI, which couples existing mod- els to simulate the full Sun-to-Earth evolution of a CME. OSPREI simulates the change in the path of a CME near the Sun, a range of effects during it’s transit to Earth, and what it would look like as it passes over a spacecraft. OSPREI runs very quickly so we can see how results change as the inputs change. To facilitate forecasts, we have devel- oped a set of standardized figures that present OSPREI results (including CME loca- tion, time of arrival, and internal properties) in an easily accessible format displaying the most likely values and the range of possibilities.

1 Introduction Forecasting space weather has gained increased importance in the public eye with recent emphasis from NASA’s Moon to Mars initiative and the US government’s National Space Weather Strategy and Action Plan in 2019 and the Promoting Research and Ob- servations of Space Weather to Improve the Forecasting of Tomorrow (PROSWIFT) Act in 2020. As CMEs are key drivers of space weather, both near the Earth and through- out the rest of the solar system, understanding and forecasting their behavior should be a part of any space weather predictions. Currently, most predictions focus simply on the arrival time of a CME. The severity of CME-driven space weather effects depends on the actual properties of the CME, such as its magnetic field strength and orientation, but currently these are not routinely included in the forecasts.

The first step to determining the potential impact of a CME is to determine where it exists in three-dimensional space. Many eruptions have a visible source on the solar

–2– MNRAS 000,1–5 (2020) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

Exploring the role of binarity in the origin of the bimodal rotational velocity distribution in stellar clusters

Sebastian Kamann,1? Nate Bastian,1 Christopher Usher,2 Ivan Cabrera-Ziri,3,1, Sara Saracino1 1Astrophysics Research Institute, Liverpool John Moores University, IC2 Liverpool Science Park, 146 Brownlow Hill, Liverpool L3 5RF, UK 2Department of Astronomy, Oskar Klein Centre, Stockholm University, AlbaNova University Centre, SE-106 91 Stockholm, Sweden 3 Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstraße 12-14, D-69120 Heidelberg, Germany

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Many young and intermediate age massive stellar clusters host bimodal distributions in the rotation rates of their stellar popula- tions, with a dominant peak of rapidly rotating stars and a secondary peak of slow rotators. The origin of this bimodal rotational distribution is currently debated and two main theories have been put forward in the literature. The first is that all/most stars are born as rapid rotators and that interacting binaries brake a fraction of the stars, resulting in two populations. The second is that the rotational distribution is a reflection of the early evolution of pre-main sequence stars, in particular, whether they are able to retain or lose their protoplanetary discs during the first few Myr. Here, we test the binary channel by exploiting multi-epoch VLT/MUSE observations of NGC 1850, a 100 Myr massive cluster in the LMC, to search for differences in the ∼ binary fractions of the slow and fast rotating populations. If binarity is the cause of the rotational bimodality, we would expect that the slowly rotating population should have a much larger binary fraction than the rapid rotators. However, in our data we detect similar fractions of binary stars in the slow and rapidly rotating populations (5.9 1.1% and 4.5 0.6%, respectively). ± ± Hence, we conclude that binarity is not a dominant mechanism in the formation of the observed bimodal rotational distributions. Key words: galaxies: star clusters: individual: NGC 1850 – stars: rotation – binaries: spectroscopic

1 INTRODUCTION suggest that the key parameter is the time scale on which pre-main sequence stars lose their proto-stellar discs, with shorter (longer) re- It has recently been established that resolved young massive clus- moval times resulting in fast rotators (slow rotators). These authors ters host stars with a bimodal rotational distribution (e.g., Dupree show that if a bimodality is set at early times it will persist for the et al. 2017; Kamann et al. 2018) with a dominant peak of rapid ro- subsequent 2 Gyr (for stars that are not magnetically braked, i.e. tators and a smaller peak of slow rotators (e.g., Marino et al. 2018; ∼ for M? & 1.5M ). Bastian et al. 2018; Kamann et al. 2020). This rotational bimodality The pre-main sequence origin of the bimodality can be tested manifests itself as a split main sequence (D’Antona et al. 2015) and through observations of the period distributions of stars in very an extended main sequence turnoff (eMSTO - Bastian & de Mink young clusters, across a wide range of stellar masses, and ongoing 2009) in colour magnitude diagrams of massive clusters with ages surveys (e.g., of the Galactic cluster Westerlund 2 - Sabbi et al. 2020) of 300 Myr and < 2 Gyr, respectively. . should be able to confirm or refute this hypothesis. The origin of this rotational bimodality is still debated in the lit- erature. D’Antona et al.(2015, 2017) suggest that all/most stars (at Alternatively, a possible link between binarity and the origin of arXiv:2109.06947v1 [astro-ph.SR] 14 Sep 2021 least those above 1.5M ) are born rapidly rotating, and that the the bimodal distribution can be tested through time-series observa- minor peak consists∼ of stars that were braked during their main se- tions of massive clusters over a wide range of ages (from forming today to 2 Gyr). In this scenario, tight binaries are needed to quence lifetimes. Motivated by an observed trend among early F ∼ to B stars in the Galactic field, which show lower rotational veloc- brake stars, and tight binaries are usually not destroyed via three- ities when being part of binary systems with periods of a few to body encounters, even in massive clusters (e.g. Heggie & Hut 2003; 500 days (Abt & Boonyarak 2004), D’Antona et al.(2015) pro- Ivanova et al. 2005; Fregeau et al. 2009). To illustrate this, we fol- pose∼ that the slow rotators observed in clusters were braked via tides low Ivanova et al.(2005) and define the hardness η of a binary as its in similarly tight binary stars. Such a scenario has also been advo- binding energy relative to the typical kinetic energy of a cluster star, 2 i.e. η = Gm1 m2/( m σ a). Here, m1 and m2 are the masses of the cated by Sun et al.(2019) for the specific case of NGC 2287, an open h i cluster in the Milky Way. binary components, orbiting each other with a semi-major axis a, m is the average stellar mass inside the cluster, σ its velocity dis- Alternatively, Bastian et al.(2020) argue that such a rotational bi- h i modality may be imprinted at very early times in a cluster’s life, and persion, and G the gravitational constant. If we consider a circular orbit with a period of 500 days (the upper limit for which a decrease in vsini is observed in the Galactic field, cf. Abt & Boonyarak 2004), ? E-mail: [email protected] typical masses of m1 = 2M and m2 = m , and a velocity disper- h i © 2020 The Authors Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX63

Thermal history of matrix forsterite grains from Murchison based on high-resolution tomography

Giulia Perotti ,1 Henning O. Sørensen ,2, 3 Henning Haack ,4 Anja C. Andersen ,1 Dario Ferreira Sanchez ,5 Elishevah M. M. E. van Kooten ,6 Esther H. R. Tsai ,5, 7 Kim N. Dalby ,3, 8 Mirko Holler ,5 Daniel Grolimund ,5 and Tue Hassenkam 3, 6

1Niels Bohr Institute, Centre for Star and Planet Formation, University of Copenhagen, Øster Voldgade 5–7, 1350 Copenhagen, Denmark 2Department of Physics, Technical University of Denmark, Fysikvej, 2800 Kongens Lyngby, Denmark 3Nano-Science Center, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark 4Maine Mineral and Gem Museum, 99 Main St, Bethel, Maine 04217, USA 5Paul Scherrer Institut, Forschungsstrasse 111, 5232 Villigen PSI, Switzerland 6Globe Institute, University of Copenhagen, Øster Voldgade 5–7, 1350 Copenhagen, Denmark 7Brookhaven National Laboratory, Upton, New York 11973, USA 8Haldor Topsøe A/S, Haldor Topsøes All´e1, 2800 Kongens Lyngby, Denmark

(Received –; Revised –; Accepted –) Submitted to ApJ

ABSTRACT Protoplanetary disks are dust- and gas-rich structures surrounding protostars. Depending on the distance from the protostar, this dust is thermally processed to different degrees and accreted to form bodies of varying chemical compositions. The primordial accretion processes occurring in the early protoplanetary disk such as chondrule formation and metal segregation are not well understood. One way to constrain them is to study the morphology and composition of forsteritic grains from the matrix of carbonaceous chondrites. Here, we present high-resolution ptychographic X-ray nanotomography and multimodal chemical micro-tomography (X-ray diffraction and X-ray fluorescence) to reveal the early history of forsteritic grains extracted from the matrix of the Murchison CM2.5 chondrite. The 3D electron density maps revealed, at unprecedented resolution (64 nm), spherical inclusions containing Fe-Ni, very little silica-rich glass and void caps (i.e., volumes where the electron density is consistent with conditions close to vacuum) trapped in forsterite. The presence of the voids along with the overall composition, petrological textures and shrinkage calculations is consistent with the grains experiencing one or more heating events with peak temperatures close to the melting point of forsterite (∼2100 K), and subsequently cooled and contracted, in agreement with chondrule-forming conditions.

Keywords: Meteorites – Methods: laboratory: solid state – Methods: laboratory: molecular – stars: formation – protoplanetary disks – planets and satellites: formation

1. INTRODUCTION (once molten silicate droplets that formed during tran- sient heating events in the solar nebula) and refractory arXiv:2109.06942v1 [astro-ph.EP] 14 Sep 2021 Carbonaceous chondrites contain materials that pro- vide key insights onto the physical and chemical prop- inclusions (i.e., the first formed solids; Connelly et al. erties of solar nebula dust grains that later agglom- 2012). erated into the building blocks of the planets in the The primordial accretion processes in the early pro- Solar System (e.g., Scott 2007). They consist of a toplanetary disk such as metal segregation, chondrule fine-grained matrix which harbours spherical chondrules formation and recycling are still poorly understood. By following a multimodal approach, we performed a high-resolution pilot study of µm-sized particles from Corresponding author: Giulia Perotti the Murchison CM2.5 chondrite to gain new insights [email protected] into their formation processes. In particular, we used a novel powerful method, Ptychographic X-ray Com- DRAFTVERSION 16TH SEPTEMBER, 2021 Typeset using LATEX twocolumn style in AASTeX63

Keck/OSIRIS Paβ high-contrast imaging and updated constraints on PDS 70b

TAICHI UYAMA ,1, 2, 3 CHEN XIE ,4 YUHIKO AOYAMA ,5, 6 CHARLES A.BEICHMAN ,2, 1 JUN HASHIMOTO ,7, 8 RUOBING DONG ,9 YASUHIRO HASEGAWA ,10 MASAHIRO IKOMA ,3 DIMITRI MAWET ,11, 10 MICHAEL W. MCELWAIN ,12 JEAN-BAPTISTE RUFFIO ,11 KEVIN R.WAGNER ,13, 14 , * JASON J.WANG ,11 ,† AND YIFAN ZHOU 15 ,‡

1Infrared Processing and Analysis Center, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125, USA 2NASA Exoplanet Science Institute, Pasadena, CA 91125, USA 3National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 4Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France 5Institute for Advanced Study, Tsinghua University, Beijing 100084, People’s Republic of China 6Department of Astronomy, Tsinghua University, Beijing 100084, People’s Republic of China 7Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 8Subaru Telescope, National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 9Department of Physics & Astronomy, University of Victoria, Victoria, BC, V8P 1A1, Canada 10Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA, 91109, USA 11Department of Astronomy, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125, USA 12Exoplanets and Stellar Astrophysics Laboratory, Code 667, Goddard Space Flight Center, Greenbelt, MD 20771, USA 13Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 14NASA NExSS Alien Earths Team, USA 15Department of Astronomy/McDonald Observatory, The University of Texas, 2515 Speedway, Stop C1400 Austin, TX 78712, USA

Submitted to AAS Journals

ABSTRACT We present a high-contrast imaging search for Paβ line emission from protoplanets in the PDS 70 system with Keck/OSIRIS integral field spectroscopy. We applied the high-resolution spectral differential imaging technique to the OSIRIS J-band data but did not detect the Paβ line at the level predicted using the parameters of Hashimoto et al.(2020). This lack of Pa β emission suggests the MUSE-based study may have overestimated the line width of Hα. We compared our Paβ detection limits with the previous Hα flux and Hβ limits and estimated AV to be ∼ 0.9 and 2.0 for PDS 70 b and c respectively. In particular, PDS 70 b’s AV is much smaller than implied by high-contrast near-infrared studies, which suggests the infrared-continuum photosphere and the hydrogen-emitting regions exist at different heights above the forming planet.

Keywords: Planet Formation, Exoplanet Astronomy

1. INTRODUCTION between various planet formation mechanisms. Characteri- A variety of theoretical and observational studies have zation of the physical and atmospheric parameters of proto- investigated planet formation, yet the mechanisms are still planets at specific ages helps in assessing the initial condi- arXiv:2109.06930v1 [astro-ph.EP] 14 Sep 2021 poorly understood. High-contrast imaging at infrared (IR) tions of their formation (e.g. Bonnefoy et al. 2014). Further- wavelengths can detect the thermal emission of young exo- more, addressing such problems as reconciling the evolution- planets directly and thus provide key insights to distinguish ary cooling models (hot/warm/cold-start; Spiegel & Burrows 2012) with the relevant physical processes (e.g. core accre- tion and disk instability; Pollack et al. 1996; Boss 1997) are Corresponding author: Taichi Uyama essential to improving our understanding of planet formation. [email protected] One of the ways to probing planet formation is to observe * NASA Hubble Fellowship Program – Sagan Fellow hydrogen emission originating in active mass accretion onto † 51 Pegasi b Fellow protoplanets (Aoyama et al. 2018). ‡ Harlan J. Smith McDonald Observatory Fellow Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX63

Robustness of Synthetic Observations in Producing Observed Core Properties: Predictions for the TolTEC Clouds to Cores Legacy Survey

S. K. Betti ,1 R. Gutermuth ,1 S. Offner ,2 G. Wilson,1 A. Sokol,1 and R. Pokhrel 3

1Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 2Department of Astronomy, University of Texas at Austin, Austin, TX 78712, USA 3Department of Astronomy, University of Toledo, Toledo, OH 43606, USA

ABSTRACT We use hydrodynamical simulations of star-forming gas with stellar feedback and sink particles (proxies for young stellar objects, i.e., YSOs) to produce and analyze synthetic 1.1 mm continuum observations at different distances (150 - 1000 pc) and ages (0.49 - 1.27 Myr). We characterize how the inferred core properties, including mass, size, and clustering with respect to diffuse natal gas structure, change with distance, cloud evolution, and the presence of YSOs. We find that atmospheric filtering and core segmentation treatments have distance-dependent impacts on the resulting core properties for d < 300 pc and 500 pc, respectively, which dominate over evolutionary differences. Concentrating on synthetic observations at further distances (650-1000 pc), we find a growing separation between the inferred sizes and masses of cores with and without YSOs in the simulations, which is not seen in recent observations of the Mon R2 cloud at 860 pc. We find that the synthetic cores cluster in smaller groups, and their mass densities are correlated with gas column density over a much narrower range, than the MonR2 observations. Such differences limit applicability of the evolutionary predictions we report here and motivate future efforts to adapt our synthetic observation and analysis framework to next generation simulations such as STARFORGE. These predictions and systematic characterizations will help guide analysis of cores for the upcoming TolTEC Clouds to Cores Legacy Survey on the Large Millimeter Telescope Alfonso Serrano (LMT).

Keywords: stars:formation - ISM:clouds - stars:protostars

1. INTRODUCTION 104−5 cm−3 in density, 0.03-0.1 pc in size, and T < 12 Dense cores within giant molecular clouds (GMC) K in temperature (di Francesco et al. 2007). Through are the birthplaces of stars (Bergin & Tafalla 2007; gravitational instabilities, some prestellar cores may col- di Francesco et al. 2007; Ward-Thompson et al. 2007). lapse to form protostars. These cores are part of a hierarchical structure of frag- In the millimeter to submillimeter regime, thermal ments within GMCs (Pokhrel et al. 2018) that form emission from cores is optically thin, which allows the as a result of a variety of physical processes, includ- flux density to trace the total dust mass, making them ing self-gravity (Heyer et al. 2009; Ballesteros-Paredes ideal targets of study for both ground-based and space- et al. 2011, 2012), magnetohydrodynamic turbulence based submm telescopes. Various surveys of nearby arXiv:2109.06916v1 [astro-ph.GA] 14 Sep 2021 (Mac Low & Klessen 2004; Hennebelle & Falgarone Gould Belt molecular clouds (d < 500 pc) have pro- 2012), supersonic interstellar turbulence (e.g. Pudritz & duced censuses of dense cores (e.g., K¨onyves et al. 2015; Kevlahan 2013), and ionization of molecular gas (Whit- Marsh et al. 2016; Bresnahan et al. 2018; Benedettini worth et al. 1994; Dale et al. 2009). The densest filamen- et al. 2018; Ladjelate et al. 2020; Enoch et al. 2006, tary structures within clouds have been observed to host 2007, 2008). These surveys also confirmed that the cores (Andr´eet al. 2010; Polychroni et al. 2013) of order cores predominantly are found within relatively dense ∼ 0.1 pc filamentary structures in clouds (Arzoumanian et al. 2011, 2019; Andr´eet al. 2010, 2014; Men’shchikov Corresponding author: S. K. Betti et al. 2010; Polychroni et al. 2013). These cores have [email protected] been observed in varying stages of gravitational stabil- ity, from those that are gravitationally collapsing and MNRAS 000,1–14 (2021) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

The Eddington ratio-dependent ‘changing look’ events in NGC 2992

Muryel Guolo1,2★, Daniel Ruschel-Dutra1, Dirk Grupe3, Bradley M. Peterson4,5,6, Thaisa Storchi-Bergmann7, Jaderson Schimoia8, Rodrigo Nemmen9, Andrew Robinson10 1Departamento de Física - CFM - Universidade Federal de Santa Catarina, 476, 88040-900 Florianópolis, SC, Brazil 2Department of Physics and Astronomy, Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA 3Space Science Center, Morehead State University, 235 Martindale Drive, Morehead, KY 40351, USA 4 Department of Astronomy, The Ohio State University, 140 W 18th Ave, Columbus, OH 43210, USA 5 Center for Cosmology and AstroParticle Physics, 191 Woodruff Ave., Columbus, OH 43210, USA 6 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 7 Departamento de Astronomia, Universidade Federal do Rio Grande do Sul. Av. Bento Goncalves 9500, 91501-970 Porto Alegre, RS, Brazil 8Departamento de Física - CCNE - Universidade Federal de Santa Maria, 97105-90, Santa Maria, RS, Brazil 9 Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, São Paulo, SP 05508-090, Brazil 10 School of Physics and Astronomy, Rochester Institute of Technology, 84 Lomb Memorial Drive, Rochester, NY 14623, USA

Accepted 2021 September 06. Received 2021 August 24; in original form 2020 November 03

ABSTRACT

We present an analysis of historical multi-wavelength emission of the Changing Look (CL) Active Galactic Nucleus (AGN) in NGC 2992, covering epochs ranging from 1978 to 2021, as well as new X-ray and optical spectra. The galaxy presents multiple Seyfert type transitions from type 2 to intermediate-type, losing and regaining its H훼 BEL recurrently. In X-rays, the source shows intrinsic variability with the absorption corrected luminosity varying by a factor of ∼ 40. We rule-out tidal disruption events or variable obscuration as causes of the type transitions, and show that the presence and the flux of the broad H훼 emission line are 42 −2 −1 directly correlated with the 2-10 keV X-ray luminosity (L2−10): the component disappears at L2−10 ≤ 2.6×10 erg cm s , this value translates into an Eddington ratio (휆Edd) of ∼ 1%. The 휆Edd in which the BEL transitions occur is the same as the critical value at which there should be a state transition between a radiatively inefficient accretion flow (RIAF) and a thin accretion disk, such similarity suggests that the AGN is operating at the threshold mass accretion rate between the two accretion modes. We find a correlation between the narrow Fe K훼 flux and 휆Edd, and an anti-correlation between full-width at half maximum of H훼 BEL and 휆Edd, in agreement with theoretical predictions. Two possible scenarios for type transitions are compatible with our results: either the dimming of the AGN continuum, which reduces the supply of ionising photons available to excite the gas in the Broad Line Region (BLR), or the fading of the BLR structure itself occurs as the low accretion rate is not able to sustain the required cloud flow rate in a disk-wind BLR model. Key words: galaxies: Seyfert – galaxies: nuclei – X-rays: galaxies – optical: galaxies

1 INTRODUCTION 103 km s−1) permitted emission lines and narrow forbidden lines in their optical spectra; while Type 2 AGNs lack BELs. Intermediate

arXiv:2109.06914v1 [astro-ph.GA] 14 Sep 2021 Supermassive black holes (SMBH) at the centre of galaxies can AGN types were added by Osterbrock(1981), including types 1.8 become active by the accretion of surrounding matter (see Storchi- and 1.9, classified, respectively, according to the presence of a weak, Bergmann & Schnorr-Müller 2019, for a review on the feeding mech- or absent, H훽 BEL, while retaining an H훼 BEL. anisms). Such active galactic nuclei (AGN) can be identified by fea- tures in their optical spectra - such as broad emission lines (BEL) AGN unification models attribute differences in AGN type to the and/or narrow high ionisation emission lines due to the non-stellar viewing angle toward an axisymmetric, parsec scale obscuring struc- ionising continuum - or by their nuclear X-ray emission. Emission ture, which can block our direct view of the AGN central engine (An- lines produced by the gas present in the vicinity of the SMBH are tonucci 1993; Urry & Padovani 1995) and is responsible for the lack Doppler broadened, forming the broad-line region (BLR), while the of BEL in Type 2 AGNs. However, a relatively small number of AGN narrow-line region (NLR) extends from a few parsecs to hundreds show transitions from one type to another in a few years; these objects of parsecs beyond the accretion disk. Historically, Type 1 AGNs are are known as “changing look” (CL) AGNs. CL-AGNs have changed those showing both broad (> 103 , < 104 km s−1) and narrow (< the widely accepted AGN paradigm, not only in the orientation-based unified model but also standard accretion disk models (see Lawrence 2018, for a review of the so-called Quasar Viscosity Crisis). Their ★ E-mail:[email protected] extreme and fast variability has motivated alternatives and/or modifi-

© 2021 The Authors Draft version September 16, 2021 Typeset using LATEX default style in AASTeX63

Where to find over-massive brown dwarfs: new benchmark systems for binary evolution

Dorsa Majidi,1, 2 John C. Forbes,3 and Abraham Loeb4

1Department of Astronomy, University of Wisconsin, 475 North Charter Street, Madison, WI 53706 USA 2UCLA Physics & Astronomy Department, 475 Portola Plaza, Los Angeles, CA, 90095 USA 3Flatiron Institute, 162 5th Avenue, New York, NY, 10010 4Center for Astrophysics | Harvard and Smithsonian, 60 Garden Street, Cambridge, MA, 02138

Submitted to AAS Journals

ABSTRACT Under the right conditions brown dwarfs that gain enough mass late in their lives to cross the hydrogen burning limit will not turn into low-mass stars, but rather remain essentially brown dwarf-like. While these objects, called either beige dwarfs or over-massive brown dwarfs, may exist in principle, it remains unclear exactly how they would form astrophysically. We show that accretion from AGB winds, aided by the wind Roche lobe overflow mechanism, is likely to produce a substantial population of observable overmassive brown dwarfs, though other mechanisms are still plausible. Specifically we predict that sun-like stars born with a massive brown dwarf companion on an orbit with a semi-major axis of order 10 AU will likely produce overmassive brown dwarfs, which may be found today as companions to the donor star’s remnant white dwarf. The identification and characterization of such an object would produce unique constraints on binary evolution because there is a solid upper limit on the brown dwarf’s initial mass.

1. INTRODUCTION Brown dwarfs were predicted on theoretical grounds by S. Kumar in the early 1960s (Kumar 1962a,b, 1963a,b, 2003). It took several decades before searches for brown dwarfs (e.g. Latham et al. 1989) yielded a definitive discovery (Nakajima et al. 1995). In the intervening decades, surveys have identified thousands of brown dwarf candidates, particularly in the infrared and taking advantage of proper motion measurements (Oppenheimer et al. 1995; Kirkpatrick et al. 1999; Lucas & Roche 2000; Gizis et al. 2002, 2000; Cruz et al. 2003; Close et al. 2003; Faherty et al. 2009; Kirkpatrick et al. 2011; West et al. 2011). Brown dwarfs have proved enigmatic as the “missing link” between planets and low-mass stars, with some ambiguity remaining about whether brown dwarfs tend to form like scaled-up planets (Stamatellos & Whitworth 2009), or scaled- down stars (Reipurth & Clarke 2001; Bate et al. 2002). The atmospheres of brown dwarfs have also proven incredibly challenging to model as the result of chemistry and cloud physics (Marley et al. 1996; Burrows et al. 1997; Allard et al. 1997, 2011, 2012), which in turn affects their evolution (Chabrier & Baraffe 1997, 2000; Chabrier et al. 2000; Burrows et al. 2001; Baraffe et al. 2003) and hinders the inference of basic properties of field brown dwarfs (Filippazzo et al.

arXiv:2109.06899v1 [astro-ph.SR] 14 Sep 2021 2015). Around the time of their discovery, brown dwarfs and even brown dwarfs exceeding the hydorgen burning limit (HBL) were proposed as constituents of dark matter (Salpeter 1992; Hansen 1999; Lynden-Bell & Tout 2001), but stringent limits on massive compact halo objects from microlensing (e.g. Alcock et al. 1993, 1997) have made this explanation for dark matter untenable. Nonetheless, the possibility that brown dwarfs more massive than the HBL exist remains intriguing but unrealized thus far observationally. Recently (Forbes & Loeb 2019, hereafter FL19) re-discovered this possibility as a consequence of the diagram of equilibrium hydrogen burning as a function of a low-mass star’s central

Corresponding author: Dorsa Majidi [email protected] MNRAS 000,1–19 (2020) Preprint 16 September 2021 Compiled using MNRAS LATEX style file v3.0

Stellar feedback in a clumpy galaxy at 푧 ∼ 3.4

E. Iani1,2,3★, A. Zanella4, J. Vernet3, J. Richard5, M. Gronke6†, C. M. Harrison7, F. Arrigoni-Battaia8, G. Rodighiero2, A. Burkert9,10, M. Behrendt9,10, Chian-Chou Chen11, E. Emsellem3,12, J. Fensch5, P. Hibon13, M. Hilker3, E. Le Floc’h14, V. Mainieri3, A. M. Swinbank15,16, F. Valentino17,18, E. Vanzella19 and M. A. Zwaan3 Affiliations are listed at the end of the paper

Accepted 2021 July 27. Received 2021 July 8; in original form 2021 February 24

ABSTRACT Giant star-forming regions (clumps) are widespread features of galaxies at 푧 ≈ 1 − 4. Theory predicts that they can play a crucial role in galaxy evolution if they survive to stellar feedback for > 50 Myr. Numerical simulations show that clumps’ survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps’ outflows strength and gas removal timescale are still uncertain. In this context, we study a line-emitting galaxy at redshift 푧 ' 3.4 lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes . 280 pc and representative of the low-mass end of clumps’ mass distribution (stellar 8 masses . 2 × 10 M ) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population (∼ 10 Myr). The properties of the Lyman-훼 (Ly훼) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of ∼ 200 − 300 km/s. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor (휂) for the clumps of ∼ 1.8 − 2.4 consistent with results obtained from hydro-dynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semi-closed box model), the estimates of 휂 suggest that the timescale 8 over which the outflows expel the molecular gas reservoir (' 7 × 10 M ) of the four detected low-mass clumps is . 50 Myr. Key words: galaxies: evolution – galaxies: high-redshift – galaxies: irregular – galaxies: ISM – galaxies: star formation

1 INTRODUCTION dominate their light profile, thus making the clumps’ host galaxies generally referred to as clumpy galaxies. Deep rest-frame ultraviolet (UV) and optical observations In the last decade, many efforts have been devoted to the un- (e.g. Driver et al. 1995; Glazebrook et al. 1995; van den Bergh et al. derstanding of clumps’ nature and properties. Clumps have been 1996; Driver et al. 1998; Elmegreen et al. 2007, 2009; Overzier et al. detected in rest-frame UV imaging (e.g. Guo et al. 2012, 2015, 2010; Swinbank et al. 2010; Förster Schreiber et al. 2011b; Genzel 2018; Livermore et al. 2012; Shibuya et al. 2016; Soto et al. 2017; et al. 2011; Wuyts et al. 2012; Guo et al. 2012; Conselice 2014;

arXiv:2109.06898v1 [astro-ph.GA] 14 Sep 2021 Dessauges-Zavadsky & Adamo 2018; Messa et al. 2019; Vanzella Tadaki et al. 2014; Murata et al. 2014; Guo et al. 2015; Shibuya et al. 2021) as well as in maps of Balmer (e.g. H훼,H훽, see Liver- et al. 2016; Soto et al. 2017; Fisher et al. 2017; Guo et al. 2018) more et al. 2012; Mieda et al. 2016; Fisher et al. 2017; Zanella et al. have revealed that galaxies at the cosmic noon (redshift 푧 ∼ 1 − 3) 2019; Whitmore et al. 2020) and Paschen (e.g. Pa훼, see Larson typically display higher gas fractions (Daddi et al. 2010; Tacconi et al. 2020) transitions. They also have been found to contribute to et al. 2010, 2013; Genzel et al. 2015), star formation rates (SFRs, their host galaxies optical continuum (e.g. Elmegreen et al. 2009; e.g. Genzel et al. 2006; Förster Schreiber et al. 2006; Genzel et al. Förster Schreiber et al. 2011b,a) and CO emissions (e.g. Jones et al. 2008), and velocity dispersions (Elmegreen & Elmegreen 2005; 2010; Dessauges-Zavadsky et al. 2017). Förster Schreiber et al. 2006) than local star-forming galaxies. Fur- Observations showed that clumps have sizes . 1 kpc thermore, bright concentrations of light, the so-called clumps, often (e.g. Elmegreen et al. 2007; Förster Schreiber et al. 2011b), es- 7 9 timated stellar masses (M★) of ∼ 10 − 10 M (e.g. Förster Schreiber et al. 2011a; Guo et al. 2012; Soto et al. 2017), and SFR ★ E-mail: [email protected] from 0.1 − 10 M /yr (e.g. Guo et al. 2012; Soto et al. 2017). † Hubble Fellow Evidence also suggests that clumps are starbursting, i.e. they have

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

The effect of inhomogeneous reionisation on the Lyman-훼 forest power spectrum at redshift 푧 > 4: implications for thermal parameter recovery

Margherita Molaro1★, Vid Iršič2, James S. Bolton1, Laura C. Keating3, Ewald Puchwein3, Prakash Gaikwad2, Martin G. Haehnelt2, Girish Kulkarni4 & Matteo Viel5,6,7,8 1School of Physics and Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 2Kavli Institute for Cosmology and Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA, UK 3Leibniz-Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 4Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 5SISSA - International School for Advanced Studies, Via Bonomea 265, I-34136 Trieste, Italy 6IFPU, Institute for Fundamental Physics of the Universe, Via Beirut 2, I-34151 Trieste, Italy 7INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, I-34131 Trieste, Italy 8INFN - National Institute for Nuclear Physics, Via Valerio 2, I-34127 Trieste, Italy

16 September 2021

ABSTRACT We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the Ly훼 forest transmitted flux at redshifts 4.2 ≤ 푧 ≤ 5. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, 푘 ∼ 0.1 km−1 s, by ∼ 10 per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, 푘 < 0.03 km−1 s, the power spectrum is instead enhanced by 10–50 per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly훼 forest power spectrum at 4.1 ≤ 푧 ≤ 5.4 in a form that can be easily applied within other parameter inference frameworks. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly훼 forest data at 푧 > 4, 0.006 km−1 s ≤ k ≤ 0.2 km−1 s, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of ∼ 10 per cent. However, for 5 per cent uncertainties, ∼ 1휎 shifts between the estimated and true parameters occur. Key words: cosmology: large-scale structure of Universe – methods: numerical – galaxies: intergalactic medium – QSOs: absorption lines

arXiv:2109.06897v1 [astro-ph.CO] 14 Sep 2021 1 INTRODUCTION et al. 2020; Keating et al. 2020; Nasir & D’Aloisio 2020; Qin et al. 2021). Additionally, because the Ly훼 forest closely tracks the dark The study of Ly훼 absorption features in the spectra of bright, high- matter down to scales of ∼ 100 ckpc, it is also sensitive to the redshift quasars offers a valuable insight into the physical properties suppression of the matter power spectrum on small scales (see e.g. of the intergalactic medium (IGM) in the early Universe (Fan et al. Seljak et al. 2006; Viel et al. 2008; Boyarsky et al. 2009). The power 2006; Mortlock et al. 2011; Becker et al. 2015; Eilers et al. 2017; spectrum of the Ly훼 forest transmitted flux thus provides one of the Bosman et al. 2018; Yang et al. 2020). Cumulatively referred to as tightest lower limits on the mass of a putative warm dark matter the “Ly훼 forest”, these absorption features provide constraints on thermal relic (Viel et al. 2013; Iršič et al. 2017a; Garzilli et al. the ionisation and thermal state of the IGM during and immediately 2019; Palanque-Delabrouille et al. 2020; Rogers & Peiris 2021). after the final stages of reionisation at redshifts 푧 ' 5 − 7 (Oñorbe et al. 2017b; Kulkarni et al. 2019; Walther et al. 2019; Gaikwad Existing constraints on reionisation and warm dark matter from the Ly훼 forest typically rely on high quality quasar spectra com- bined with numerical models for the distribution of matter in the ★ E-mail: [email protected] IGM (see e.g. McQuinn 2016). On the observational side, the num-

© 2021 The Authors MNRAS 000,1–16 (2021) Preprint 16 September 2021

Cosmic metal density evolution in neutral gas: insights from observations and cosmological simulations

Robert M. Yates,1★ Céline Péroux,2,3 Dylan Nelson4

1Department of Physics, University of Surrey, Stag Hill, Guildford, GU2 7XH, UK 2European Southern Observatory, Karl-Schwarzschildstrasse 2, D-85748 Garching bei München, Germany 3Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 4Universität Heidelberg, Zentrum für Astronomie, Institut für theoretische Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany

ABSTRACT We contrast the latest observations of the cosmic metal density in neutral gas (휌met,neu) with three cosmological galaxy evolution simulations: L-Galaxies 2020, TNG100, and EAGLE. We find that the fraction of total metals that are in neutral gas is < 40 per cent at 3 . 푧 . 5 in these simulations, whereas observations of damped Lyman-훼 (DLA) systems suggest & 85 per cent. In all three simulations, hot, low-density gas is also a major contributor to the cosmic metal budget, even at high redshift. By considering the evolution in cosmic SFR density (휌sfr), neutral gas density (휌HI), and mean gas-phase metallicity ([hM/Hi]neu), we determine two possible ways in which the absolute 휌met,neu observed in DLAs at high redshift can be matched by simulations: (a) the 휌sfr at 푧 & 3 is greater than inferred from current FUV observations, or (b) current high-redshift DLA metallicity samples have a higher mean host mass than the overall galaxy population. If the first is correct, TNG100 would match the ensemble data best, however there would be an outstanding tension between the currently observed 휌sfr and 휌met,neu. If the second is correct, L-Galaxies 2020 would match the ensemble data best, but would require an increase in neutral gas mass inside subhaloes above 푧 ∼ 2.5. If neither is correct, EAGLE would match the ensemble data best, although at the expense of over-estimating [hM/Hi]neu. Modulo details related to numerical resolution and Hi mass modelling in simulations, these incompatibilities highlight current tensions between key observed cosmic properties at high redshift. Key words: galaxies: formation – galaxies: evolution – galaxies: abundances – methods: numerical

1 INTRODUCTION lower than the total expected metal density in the Universe (휌met,tot) at 푧 ∼ 2.5, and constitutes & 85 per cent at 푧 & 4. These estimates are The evolution of the cosmic metal density (휌 ) and cosmic star- met higher than previously predicted from comparisons between DLAs formation rate density (휌 ) are key diagnostics for understanding sfr and Lyman-break galaxies (LBGs) (Rafelski et al. 2014). galaxy evolution. Combined, they encode the relative significance of all the main physical processes driving galaxy evolution, including Such a result is perhaps surprising, given that the presence of gas accretion, star formation, stellar evolution, metal enrichment, strong galactic winds in highly-star-forming galaxies and AGN hosts and outflows driven by supernovae (SNe) and active galactic nuclei at high redshift (e.g. Adelberger et al. 2003; Steidel et al. 2010; (AGN). Bischetti et al. 2019; Sugahara et al. 2019; Ginolfi et al. 2020; Spilker Recently, observational studies of the metal content in and around et al. 2020) suggest an efficient removal of metals from galaxies. arXiv:2109.06888v1 [astro-ph.GA] 14 Sep 2021 galaxies have suggested that nearly all the metals in the Universe at This material should enrich the circumgalactic medium (CGM) and 푧 & 2.5 reside in regions dominated by neutral gas (Péroux & Howk intracluster medium (ICM) to high levels by 푧 ∼ 2 (e.g. Mushotzky & 2020, hereafter PH20). These regions have Hi column densities of Loewenstein 1997; Tozzi et al. 2003; Balestra et al. 2007; Anderson 푁 (Hi) > 2 × 1020 cm−2 and are probed in absorption against bright et al. 2009; Baldi et al. 2012; McDonald et al. 2016), or at least be background sources as damped Lyman-훼 (DLA) systems. DLAs are transferred to a hot, ionized phase, undetected by DLAs. The recent expected to trace the cold and warm neutral medium within galaxies results on 휌met,neu seem to be in tension with this scenario. (Wolfe et al. 2005; Bird et al. 2014; Rahmati et al. 2015; Peeples A precise census of the distribution of metals in the Universe at et al. 2019), and do not suffer from the selection biases or metallicity high redshift is difficult to obtain from current observations, as the calibration issues that hamper emission-based samples. An analysis full range of dominant contributors to the total metal budget have of this DLA data by PH20 argues that the cosmic, dust-corrected, yet to be comprehensively measured (see PH20). Therefore, in this metal mass density in neutral gas (휌met,neu) is only a factor of a few work, we turn to large-volume cosmological galaxy evolution simu- lations to assess the apportionment of metals in and around galaxies back to high redshift from a theoretical perspective. This allows us ★ E-mail: [email protected] to determine under which conditions 휌met,neu dominates the univer-

© 2021 The Authors Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX62

Far-Ultraviolet Spectra of Main-Sequence O Stars at Extremely Low Metallicity

O. Grace Telford,1 John Chisholm,2 Kristen B. W. McQuinn,1 and Danielle A. Berg2

1Rutgers University, Department of Physics and Astronomy, 136 Frelinghuysen Road, Piscataway, NJ 08854, USA 2Department of Astronomy, The University of Texas at Austin, 2515 Speedway, Stop C1400, Austin, TX 78712-1205, USA

ABSTRACT Metal-poor massive stars dominate the light we observe from star-forming dwarf galaxies and may have produced the bulk of energetic photons that reionized the universe at high redshift. Yet, the rarity of observations of individual O stars below the 20% solar metallicity (Z ) of the Small Magellanic Cloud (SMC) hampers our ability to model the ionizing fluxes of metal-poor stellar populations. We present new Hubble Space Telescope far-ultraviolet (FUV) spectra of three O-dwarf stars in the galaxies Leo P (3% Z ), Sextans A (6% Z ), and WLM (14% Z ). We quantify equivalent widths of photospheric metal lines and strengths of wind-sensitive features, confirming that both correlate with metallicity. We infer the stars’ fundamental properties by modeling their FUV through near-infrared spectral energy distributions and identify stars in the SMC with similar properties to each of our targets. Comparing to the FUV spectra of the SMC analogs suggests that (1) the star in WLM has an SMC- like metallicity, and (2) the most metal-poor star in Leo P is driving a much weaker stellar wind than its SMC counterparts. We measure projected rotation speeds and find that the two most metal-poor stars have high v sin(i) ≥ 290 km s−1, and estimate just a 3−6% probability of finding two fast rotators if the metal-poor stars are drawn from the same v sin(i) distribution observed for O dwarfs in the SMC. These observations suggest that models should include the impact of rotation and weak winds on ionizing flux to accurately interpret observations of metal-poor galaxies in both the near and distant universe. 1. INTRODUCTION Robertson et al. 2015; Finkelstein et al. 2019), but their 1.1. Uncertainties in Modeling Metal-Poor Massive ionizing fluxes remain uncertain. Stellar Populations The neutral intergalactic medium is opaque to H- ionizing photons (λ ≤ 912 A)˚ during the epoch of reion- Massive stars dominate the spectral energy distribu- ization (e.g., Becker et al. 2001). Observations of high-z tions (SEDs) of star-forming galaxies. Throughout their galaxies at longer rest wavelengths must therefore be lives and during their explosive deaths as supernovae, modeled with stellar population synthesis (SPS) codes they deposit energy, momentum, and ionizing photons to infer the ionizing flux from massive stars. SPS models into the surrounding interstellar medium (ISM). This combine template stellar spectra, stellar evolution mod- feedback from massive stellar populations is widely rec- els, an initial mass function (IMF), and a star formation ognized as an important regulator of star formation history (SFH) parameterization to predict the intrinsic and galaxy evolution processes (e.g., Somerville & Dav´e SED of a stellar population (Tinsley 1980; Conroy 2013). 2015; Naab & Ostriker 2017). Bursts of star forma- From this, the ionizing flux from massive stars can be tion can efficiently create low-density channels in the

arXiv:2109.06885v1 [astro-ph.SR] 14 Sep 2021 inferred, then combined with observational constraints ISM through which ionizing photons produced by young, on ionizing photon escape fractions (e.g., Izotov et al. massive stars can escape, particularly from low-mass 2016; Steidel et al. 2018) and the galaxy luminosity func- galaxies in the early universe (e.g., Wise & Cen 2009; tion at high z (e.g., Livermore et al. 2017; Oesch et al. Trebitsch et al. 2017). Thus, massive stars in metal- 2018) to predict the total ionizing photons escaping from poor dwarf galaxies are likely drivers of the reionization the predominantly metal-poor galaxies in the early uni- of the universe at high redshift, z & 6 (Ouchi et al. 2009; verse. Our understanding of the contribution of galaxies to cosmic reionization is therefore highly sensitive to the [email protected] stellar evolution models and theoretical spectra adopted in SPS codes. MNRAS submitted,1–17 (2021) 16 September 2021

Predictions for anisotropic X-ray signatures in the circumgalactic medium: imprints of supermassive black hole driven outflows

Nhut Truong1★, Annalisa Pillepich1, Dylan Nelson2, Norbert Werner3, and Lars Hernquist4 1Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany 2Universität Heidelberg, Zentrum für Astronomie, Institut für theoretische Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany 3Department of Theoretical Physics and Astrophysics, Faculty of Science, Masaryk University, Kotlářská 2, Brno, 611 37, Czech Republic 4Institute for Theory and Computation, Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

ABSTRACT The circumgalactic medium (CGM) encodes signatures of the galaxy-formation process, including the interaction of galactic outflows driven by stellar and supermassive black hole (SMBH) feedback with the gaseous halo. Moving beyond spherically symmetric radial profiles, we study the angular dependence of CGM properties around 푧 = 0 massive galaxies in the IllustrisTNG simulations. We characterize the angular signal of density, temperature, and metallicity of the CGM as a function of galaxy stellar mass, halo mass, distance, and SMBH mass, via stacking. TNG predicts that the CGM is anisotropic in its thermodynamical 10−11.5 properties and chemical content over a large mass range, 푀∗ ∼ 10 푀 . Along the minor axis directions, gas density is diluted, whereas temperature and metallicity are enhanced. These feedback-induced anisotropies in the CGM have a magnitude 10.8 of 0.1−0.3 dex, extend out to the halo virial radius, and peak at Milky Way-like masses, 푀∗ ∼ 10 푀 . In TNG, this mass scale corresponds to the onset of efficient SMBH feedback and the production of strong outflows. By comparing the anisotropic signals predicted by TNG versus other simulations – Illustris and EAGLE – we find that each simulation produces distinct signatures and mass dependencies, implying that this phenomenon is sensitive to the underlying physical models. Finally, we explore X-ray emission as an observable of this CGM anistropy, finding that future X-ray observations, including the eROSITA all-sky survey, will be able to detect and characterize this signal, particularly in terms of an angular modulation of the X-ray hardness. Key words: galaxies: evolution – galaxies: formation – galaxies: haloes – galaxies: circumgalactic medium–galaxies: super- massive black holes — X-ray: galaxies — methods: numerical

1 INTRODUCTION Without AGN feedback, simulations produce overly-massive galax- ies inconsistent with observations (Khalatyan et al. 2008; Le Brun The circumgalactic medium (CGM) is the gaseous atmosphere that et al. 2014; Pillepich et al. 2018a). In addition, AGN feedback is extends far beyond the stellar body of a galaxy, up to hundreds of necessary for simulations to reproduce the observed thermodynami- kiloparsecs in galactocentric distance, and plays an essential role cal properties of halo gas, from galaxies to groups and cluster mass in the process of galaxy evolution (see Tumlinson et al. 2017 for a scales (McCarthy et al. 2010; Le Brun et al. 2014; Henden et al. review). The CGM is a unique laboratory where the cosmic baryon 2018; Nelson et al. 2018b; Truong et al. 2018, 2020). cycle of gas inflows and gas outflows plays out. While gas accretion Observations have inferred that most, if not all, massive galaxies provides the fuel for future star formation and the overall growth of contain a SMBH at their centers (Kormendy & Ho 2013; van den galaxies, gaseous outflows act in opposition, removing baryons from Bosch 2016). Nonetheless, both observationally and theoretically, the the dense star-forming interstellar medium. Galactic-scale outflows, exact nature of how AGN feedback interacts with the galaxy, much also referred to as galactic winds, are sourced primarily by stellar less the surrounding CGM, remains unclear. feedback in low-mass galaxies (Zhang 2018), active galactic nuclei

arXiv:2109.06884v1 [astro-ph.GA] 14 Sep 2021 Observationally, it is suggested that AGN feedback appears in two (AGN) feedback in more massive systems (Morganti 2017), or a modes: “quasar” and “radio” (Croton et al. 2006; Cattaneo et al. 2009; combination of the two at intermediate masses. Fabian 2012; McNamara & Nulsen 2012; Harrison 2017; Morganti AGN feedback, powered by gas accretion onto supermassive black 2017). The quasar (or radiative) mode is mostly found in galaxies holes (SMBHs), is an important ingredient in current theoretical hosting luminous AGN, where the SMBH has a high accretion rate models of galaxy formation and evolution (Springel 2005; Springel & 1 − 10% of Eddington. In this case, feedback energy is mainly et al. 2005; Booth & Schaye 2009; Sijacki et al. 2015; Weinberger released in the form of radiation from the accretion disk. In contrast, et al. 2017; Davé et al. 2019). It enables cosmological simulations a radio mode may instead operate for SMBHs with low accretion to reproduce various observed galaxy properties, both of the stellar rates. In this state, observations suggest that feedback is released body as well as their gaseous atmospheres. It is thought to be the in the form of mechanical energy by means of bubbles, cavities, or physical mechanism that quenches star formation in massive galax- relativistic jets, with large amounts of energy (∼ 1060 erg/s) deposited ies (Silk & Rees 1998; Di Matteo et al. 2005; Man & Belli 2018). into the surrounding environments via e.g. turbulence and/or shocks (Nulsen et al. 2009; Randall et al. 2011; Zhuravleva et al. 2014; Hlavacek-Larrondo et al. 2015). ★ E-mail: [email protected] From the perspective of galaxy evolution simulations, successful

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

Small-scale double quasars discovered in the DECam Legacy Survey

Yu-Ching Chen1,2,3★ 1Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 West Green Street, Urbana, IL 61801, USA 2National Center for Supercomputing Applications, 1205 West Clark Street, Urbana, IL 61801, USA 3Center for AstroPhysical Surveys, National Center for Supercomputing Applications, Urbana, IL, 61801, USA

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Dual quasars are precursors of binary supermassive black holes, which are important objects for gravitational wave study, galaxy evolution, and cosmology. I report four double quasars, one of which is newly discovered, with separations of 1–2 arcsec selected from the DECam Legacy Survey. The Gemini optical spectra confirm that both sources are quasars at the same redshifts. J0118-0104 and J0932+0722 are classified as dual quasars, whereas J0037+2058 could be either a dual quasar or a lensed quasar. I estimate the physical properties such as black hole masses and Eddington ratios for the four double quasars. The newly discovered system supplements the incomplete sample of dual quasars at ∼10 kpc at 푧 > 0.5. The results demonstrate that new high-quality imaging surveys with existing spectroscopic data could reveal additional small-scale double quasars. Combined with the dual quasars from the literature, I find a power-law index of 1.45±0.48 for the distribution of dual quasars as a function of separation between 5–20 kpc, steeper than that expected from dynamical friction, though it is likely due to the incomplete sample at <15 kpc. I also discuss the implications for future searches with new imaging surveys such as the Dark Energy Survey and the Legacy Survey of Space and Time. Key words: quasars: general – surveys – galaxies: active – galaxies: nuclei – galaxies: high-redshift – black hole physics

1 INTRODUCTION 2017). However, the angular diameter distance increases dramati- cally from local universe (푧 = 0) to high-redshift universe (푧 = 2), Observational evidences support that most massive galaxies contain It is difficult to identify close dual quasars systematically in op- supermassive black holes (SMBHs; Kormendy & Richstone 1995; tical surveys due to the angular resolution limits of ground-based Kormendy & Ho 2013). Thus, galaxy mergers should lead to the telescopes. Only handful dual quasars with projected separation formation of binary supermassive black holes (BSBHs; Begelman 푟 푝 .10 kpc at redshift 푧 > 0.5 are discovered (Junkkarinen et al. et al. 1980). BSBHs are important objects for gravitational wave 2001; Inada et al. 2008, 2012; More et al. 2016; Schechter et al. study, galaxy evolution and cosmology. While the formation of 2017; Anguita et al. 2018; Lemon et al. 2018, 2020; Silverman et al. BSBHs is inevitable, direct evidence has so far been elusive. There 2020; Shen et al. 2021). are only two candidates BSBHs with separation of pc scales that were directly imaged by the Very Long Baseline Array (Rodriguez

arXiv:2109.06881v1 [astro-ph.GA] 14 Sep 2021 et al. 2006; Kharb et al. 2017) and hundreds of candidates discovered New imaging surveys with excellent image quality and im- by the indirect methods (such as periodic light curves and radial proved classification pipeline enable the discovery of small- velocity shifts, Eracleous et al. 2012; Shen et al. 2013; Liu et al. separation double quasars that are not found before. In the paper, I 2014; Graham et al. 2015; Charisi et al. 2016; Runnoe et al. 2017; observe eight close (1–2 arcsec) double quasar candidates selected Wang et al. 2017; Guo et al. 2019; Liu et al. 2019; Chen et al. 2020; from DECaLS imaging survey with Gemini and report four spec- Liao et al. 2021). troscopically confirmed double quasars. Among the four confirmed In the past decade, thanks to large spectroscopic surveys such double quasars, one is newly discovered and the other three have as the Sloan Digital Sky Survey, significant progress has been made been reported (Inada et al. 2008; More et al. 2016; Hutsemékers in finding dual quasars, the progenitors of the BSBHs when they are et al. 2020; Lemon et al. 2020). The paper is organized as fol- active and still evolving in the potential of the (merged) host galaxy lows. In section 2, I describe the data and the target selection for (Liu et al. 2011; Hennawi et al. 2006, 2010; Eftekharzadeh et al. this study. I present the main results including DECaLS images and Gemini/GMOS optical spectra and classify the individual systems in section 3. The implications are discussed in section 4. I summarize ★ Contact e-mail: [email protected] the findings in section 5. All physical separations are the projected

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

Anomalous Hydrogen Recombination-Line Ratios in Ultraluminous Infrared Galaxies

Kenichi Yano,1, 2 Shunsuke Baba ,3, ∗ Takao Nakagawa ,2 Matthew Malkan,4 Naoki Isobe,2 Mai Shirahata,2 Ryosuke Doi,1, 2 and Vanshree Bhalotia 5

1Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 2Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210, Japan 3National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 4Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1547, USA 5Department of Physics and Astronomy, University of Hawai’i at M¯anoa, 2505 Correa Rd, Honolulu, HI 96822, USA

(Received Dec 25, 2020; Revised Aug 16, 2021; Accepted Sep 13, 2021) Submitted to ApJ

ABSTRACT We conducted systematic observations of the H I Brα (4.05 µm) and Brβ (2.63 µm) lines in 52 nearby (z < 0.3) ultraluminous infrared galaxies (ULIRGs) with AKARI. Among 33 ULIRGs wherein the lines are detected, three galaxies show anomalous Brβ/Brα line ratios (∼ 1.0), which are significantly higher than those for case B (0.565). Our observations also show that ULIRGs have a tendency to exhibit higher Brβ/Brα line ratios than those observed in Galactic H II regions. The high Brβ/Brα line ratios cannot be explained by a combination of dust extinction and case B since dust extinction reduces the ratio. We explore possible causes for the high Brβ/Brα line ratios and show that the observed ratios can be explained by a combination of an optically thick Brα line and an optically thin Brβ line. We simulated the H II regions in ULIRGs with the Cloudy code, and our results show that the high Brβ/Brα line ratios can be explained by high-density conditions, wherein the Brα line becomes optically thick. To achieve a column density large enough to make the Brα line optically thick within a single H II region, the gas density must be as high as n ∼ 108 cm−3. We therefore propose an ensemble of H II regions, in each of which the Brα line is optically thick, to explain the high Brβ/Brα line ratio.

Keywords: Active galaxies — Starburst galaxies — Ultraluminous infrared galaxies

1. INTRODUCTION 1996). However, the large amount of dust harbored in a Ultraluminous infrared galaxies (ULIRGs) are charac- ULIRG makes it difficult to investigate the energy source observationally, and thus it is important to determine terized by their enormous infrared luminosities LIR (8– 12 the amount of dust extinction in a ULIRG in order to 1000 µm), which exceed 10 L (Sanders et al. 1988). Large infrared luminosity is produced by thermal radi- correct the observed quantities. ation from heated dust. This indicates that powerful One of the most widely used indicators of dust extinc- arXiv:2109.06880v1 [astro-ph.GA] 14 Sep 2021 energy sources are hidden behind dust. Since the dis- tion is the ratio of hydrogen recombination lines. These covery of ULIRGs in the 1980s, it has been debated lines have been extensively studied and are widely used whether the dominant energy source in a ULIRG is due as tracers of ionized gas because hydrogen is the simplest to starburst activity and/or to an active galactic nucleus and most abundant element in the universe (e.g., Seaton (AGN) (e.g., Sanders et al. 1988; Sanders & Mirabel 1959; Johnson 1972; Hummer & Storey 1987; Storey & Hummer 1995). The ratio of the H I line fluxes from photoionized gas can be calculated numerically for the Corresponding author: Shunsuke Baba so-called the “case B” model, in which all the Lyman- [email protected] line photons are assumed to be absorbed by other hy- ∗ JSPS Fellow (PD) drogen atoms, and all other H I lines are assumed to be optically thin. The line ratios calculated for case B Draft version September 16, 2021 Typeset using LATEX twocolumn style in AASTeX631

Multiwavelength Analysis of A1240, the Double Radio Relic Merging Galaxy Cluster Embedded in a ∼80 Mpc-long Cosmic Filament

Hyejeon Cho ,1 M. James Jee ,1, 2 Rory Smith ,3, 4 Kyle Finner ,5, 1 and Wonki Lee 1

1Department of Astronomy, Yonsei University, 50 Yonsei-ro, Seodaemun-gu, Seoul 03722, Korea 2Department of Physics and Astronomy, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA 3Korea Astronomy and Space Science Institute (KASI), 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Korea 4University of Science and Technology (UST), 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Korea 5Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA

ABSTRACT We present a multiwavelength study of the double radio relic cluster A1240 at z = 0.195. Our Subaru-based weak lensing analysis detects three mass clumps forming a ∼4 Mpc filamentary struc- +0.51 14 ture elongated in the north-south orientation. The northern (M200 = 2.61−0.60 10 M ) and mid- +0.34 14 × dle (M200 = 1.09 10 M ) mass clumps separated by ∼1.3 Mpc are associated with A1240 −0.43 × and co-located with the X-ray peaks and cluster galaxy overdensities revealed by Chandra and +0.44 14 MMT/Hectospec observations, respectively. The southern mass clump (M200 = 1.78 10 M ), −0.55 × ∼1.5 Mpc to the south of the middle clump, coincides with the galaxy overdensity in A1237, the A1240 companion cluster at z = 0.194. Considering the positions, orientations, and polarization fractions of the double radio relics measured by the LOFAR study, we suggest that A1240 is a post-merger binary system in the returning phase with the time-since-collision ∼1.7 Gyr. With the SDSS DR16 data analysis, we also find that A1240 is embedded in the much larger-scale (∼80 Mpc) filamentary structure whose orientation is in remarkable agreement with the hypothesized merger axis of A1240.

Keywords: Abell clusters (9); Cosmic web (330); Dark matter (353); Intracluster medium (858); Radio continuum emission (1340); Weak gravitational lensing (1797); X-ray astronomy (1810)

1. INTRODUCTION utilize merging clusters as dark matter laboratories, ra- Radio relics, also known as cluster radio shocks, re- dio relic clusters provide powerful environments because fer to Mpc-scale diffuse synchrotron radio emission fea- the information on their merger histories is more read- tures observed in the outskirts of merging galaxy clus- ily available (e.g., Ng et al. 2015; Golovich et al. 2016; ters. Tracing the merger shocks, they provide unique Monteiro-Oliveira et al. 2017; Finner et al. 2021). information, such as viewing angle, collision velocity, Double radio relics are a subclass of radio relics that merger axis, etc., on the merging scenarios that cannot exhibit two relics on the opposite sides. The first detec- be directly accessed with other wavelength data alone tion of double relics dates back to 1997, when R¨ottgering (e.g., Enßlin et al. 1998; see van Weeren et al. 2019 et al.(1997) found two extended diffuse radio sources in arXiv:2109.06879v1 [astro-ph.CO] 14 Sep 2021 for a recent review and references therein). There is a the northwest and southeast peripheral regions of A3667 growing interest in studying radio relic clusters to under- using ATCA and MOST data. Although in principle ev- stand the formation and evolution of the merger shocks, ery collision between clusters generates shocks in pairs, the mechanism of cosmic ray particle acceleration, the observationally there are only a dozen or so double radio interaction between plasma turbulence and intracluster relic systems reported to date. The rarity indicates that magnetic field, etc. In addition, to those who desire to the merger perhaps should meet some delicate criteria in order to produce two observable radio relics. Given that our understanding of the exact mechanism Corresponding author: M. James Jee in relic creation and particle acceleration is still highly [email protected], [email protected] incomplete, the double radio relic system serves as a powerful testbed for different models because the two