arXiv:2109.10336v1 [gr-qc] 21 Sep 2021 ntene o netnino etnster to theory Newton’s of extension topologies an non-Euclidean for effect? need relativistic the a On cosmology in backreaction Is nFR akrudepnin oee,ti supini o c a not a is solved assumption is this However, model, standard expansion. background the FLRW under an structures, th expansion of Friedmann that the dynamics implies and the relativity, Friedmann-Lemaˆıtr This general (homogeneous) of the scale. solution by (FLRW) homogeneity given affec is the not domains do than such scale of bigger Mo this size Standard under of The inhomogeneities w domains the isotropic. above that and scale assumes homogeneous a then as exists considered there be that can asserts principle cosmological The Introduction 1. Abstract. 2021 September 22 E-mail: France Lyon, F–69007, UMR5574, Lyon de Toru´n, Poland 87-100 2 5, Grudziadzka University, Copernicus 1 uni Vigneron Quentin hster sytt ejsie rmannrltvsi ii fgenera of limit non-relativistic a from justified Univ be our of to topology yet the is on non-zer theory depend is this might backreaction and Newtonian that non-relativistic locally such show be is we might for which topology, account theory non-Euclidean a a to non-Euc define on generalised we a univers where has be a calculation but im therefore to heuristic scales not result apply should small is not this on theorem this does theory that The that implication Newton’s said this show by that will generally described sense We well is the effect. in It case, relativistic the purely relativity. gravitation, a general of is in theory backreaction Newton’s case t using the stating surroun calculated 1997) be res question Ehlers, important when acceleratio & most zero main the (Buchert recent exactly of theorem The the One Buchert-Ehlers explain Energy. the could Dark is Universe. as it known not the also or factor, of whether expansion is phenomenon global the on nvLo,Esd yn nvLo1 NS eted ehrh As Recherche de Centre CNRS, Lyon1, Univ Lyon, de Ens Lyon, Univ Info and Astronomy Physics, of Faculty Astronomy, of Institute omlgclbcrato orsod oteeeto inhomoge of effect the to corresponds backreaction Cosmological uinrngalcm [email protected] [email protected], 1 , 2 eito around deviation a s e-Robertson-Walker as hti why is That laws. h yaisof dynamics the t ihteUniverse the hich neuneo the of onsequence e fCosmology of del mtc,Nicolaus rmatics, a akecinis backreaction hat l nti subject this on ult structure of neities ,maigta it that meaning o, cnro na In scenario. a ientopology. lidean relativity. l hc a not may which re However, erse. hc sstill is which e ftescale the of n expansion e trophysique u defined but necessarily le that plies igthis ding Electric Penrose process: high-energy acceleration of ionized particles by non-rotating weakly charged black hole

Arman Tursunov,∗ Bakhtinur Juraev, ZdenˇekStuchl´ık, and Martin Koloˇs Research Centre for Theoretical Physics and Astrophysics, Institute of Physics, Silesian University in Opava, Bezruˇcovon´am.13,CZ-74601 Opava, Czech Republic

In many astrophysical scenarios the charge of the black hole is often neglected due to unrealistically large values of the charge required for the Reissner-Nordstr¨omspacetime metric. However, black holes may possess small electric charge due to various selective accretion mechanisms. In this paper we investigate the effect of a small hypothetical electric charge of a Schwarzschild black hole on the ionization of a freely falling neutral particle and subsequent escape of the ionized particle from the black hole. We show that the energy of ionized particle can grow ultra-high and discuss distinguishing signatures of particle acceleration by weakly charged black holes. We also discuss a possible application of the proposed mechanism as an alternative cosmic ray acceleration scenario. In particular we show that the Galactic centre supermassive black hole is capable to act as a PeVatron of protons. The presented mechanism can serve as a simple toy model of a non-rotating compact object acting as a particle accelerator with a potential astrophysical implementations related to the cosmic ray physics and beyond.

I. WEAKLY CHARGED BLACK HOLE hole mass, spin and electric charge [5]. The later is usu- ally neglected in astrophysical scenarios justified, on the Recently, it has been pointed out that the ionization or one hand, by unrealistically large values of the charge re- decay of neutral particles in the vicinity of rotating Kerr quired for its visible effect on the spacetime metric and black hole immersed into external magnetic field can lead on the other hand by quick discharge of any charge ex- to the acceleration of ionized particles to ultra-high en- cess by an accretion of a plasma surrounding black hole. ergies, with the Lorentz γ-factors of particles that may Indeed, one can compare the gravitational radius of a exceeding 1012 near supermassive black holes in realis- black hole with the characteristic length of the charge tically plausible conditions [1, 2]. The formalism of the QG of the Reissner-Nordstr¨omblack hole, which gives 1 mentioned acceleration mechanism is based on the mag- the maximum charge of the black hole per solar mass netic Penrose process [3, 4] in its novel, ultra-efficient r Q2 G 2GM regime, in which the energy of ionized particle drives G 4 = 2 , (1) away the rotational energy of the black hole through elec- c c 1/2 30 M tromagnetic interaction. It has been claimed that the ⇒ QG = 2G M ≈ 10 Fr. (2) mechanism might be responsible for the production of M the highest-energy cosmic ray particles [1] with energy This value of the charge is unattainable in any known exceeding 1020eV when applied to realistic supermassive astrophysically relevant scenario. Moreover, a neutral- black hole candidates. The process, however, requires ization of such a hypothetical charge QG would require the rotation of the black hole and the presence of exter- an accretion of a net charge with the total mass of ac- nal magnetic field. creted charged particles of MQ = mp,eQG/e, where the In this paper we investigate whether the acceleration indices denote proton and electron, respectively, and e of ionized particles can be achieved in a more simpli- is an elementary charge. Luminosity of black hole sur- fied settings, namely, in the vicinity of a non-rotating rounded by plasma or accretion disk can be derived from arXiv:2109.10288v1 [gr-qc] 21 Sep 2021 Schwarzschild black hole with a radial test electric field. infalling matter as L = Mc˙ 2, where M˙ accretion rate By a test electric field we denote the field, whose energy- and  is the fraction of the rest mass energy radiated momentum tensor can be neglected in the description of away. On the other hand, from the balance of gravita- the gravitational field of the black hole. This implies that tional force and radiation pressure in the vicinity of a the electric field influences on the dynamics of charged black hole one can derive the Eddington luminosity for particles only, being negligible for the geodesics of neu- fully ionized hydrogen plasma surrounding a black hole tral particles. Such a simplified setup is motivated by the in the form following reasons.   4πGMmpc 38 M First of all, the no-hair theorem of black hole physics LEdd = ≈ 1.26 × 10 erg/s. (3) states that the spacetime around black holes can be fully σT M described by at most three metric parameters – black

1 Hereafter in this section we use the esu-cgs system of units, in which the electrostatic unit of charge is given by 1 Fr ≡ 1 esu = ∗ [email protected] 1 cm3/2 g1/2 s−1. In SI system of units, 1 C = 3 × 109 Fr. Life, the universe and the hidden meaning of everything

1,3 2, Zhi-Wei Wang † and Samuel L. Braunstein ∗ 1College of Physics, Jilin University, Changchun, 130012, People’s Republic of China 2Computer Science, University of York, York YO10 5GH, United Kingdom 3Department of Physics and Astronomy, University of Lethbridge, Lethbridge AB T1K 3M4, Canada †[email protected] and ∗[email protected]

It is hard to look at the universe and not wonder about the meaning, of, well, everything. A natural question is whether what we see is a sign of intelligent design. The antithesis of design would be a random universe or, assuming laws of physics, one whose fundamental physical parameters were randomly selected, but conditioned on life (ourselves) being here to observe it. In unpublished work, the British physicist Dennis Sciama argued that such a randomly selected universe would display a statistical signature. He concluded that a random universe would almost certainly have parameters only just allowing for the possibility of life. Here we consider whether this signature is definitive. We find that with plausible additional assumptions Sciama’s signature would appear to reverse: Were our universe random, it could give the false impression of being intelligently designed, with the fundamental constants appearing to be fine tuned to a strong probability for life to emerge and be maintained.

Whatever might be eventually concluded about a uni- would be the same as us being located in one of the uni- versal definition for life, we can certainly agree that the verses within the multiverse where humans are possible. universe we inhabit has so far supported the emergence, Such ‘explanations’ are said to invoke the weak anthropic evolution and continued sustenance of human beings. principle,6 yet they explain nothing and fail to provide Despite our having grown collectively more powerful than any real resolution. Are they at least predictive? most known species, within the universe we are very frag- Dennis Sciama, considered to be one of the fathers of ile and maintain a precarious hold on existence. We modern cosmology, argued that were our universe ran- are carbon based, requiring a planetary body to live on, dom it would almost certainly have a low probability for which follows a comfortable and steady orbit around a life as we know it.8 single and not too energetic star. Sciama assumed that the feature distinguishing differ- These constraints already place tight bounds on the ent potential universes was the set of specific values taken fundamental constants of the universe. To ensure that by the fundamental constants; the underlying physical a population of yellow dwarf stars like our sun exist, laws themselves being fixed. We can then envision the the fine structure constant must be tuned to within a human-compatible universes as an ‘island’ within a ‘sea’ percent or two of its current value (outside this narrow of more general possibilities. Each point on the island range almost all stars would either be blue giants, or or in the sea describes a unique universe that is de- red dwarfs).1 The cosmological constant must be between scribed by a distinct set of fundamental constants. The 120-124 orders-of-magnitude smaller than its naive quan- dimensionality of this space of points is naively given by tum field-theoretic value. The upper bound ensures that the number of fundamental constants. Thus the human- bodies can form gravitationally2,3 and the lower bound compatible island of universes corresponds to some shape ensures that nascent life will not be extinguished by prox- in a high dimensional space. The shoreline of the island imity to gamma ray bursts.4 To ensure that nuclear reac- corresponds to the boundary separating universes with a tions within stars can form carbon, but not have the pro- chance for human life to form from those where this is cess bypassed leaving only oxygen, a remarkable set of co- impossible. Thus, the shoreline itself will be made up of incidences is required among the fundamental constants universes with an exactly vanishing probability for life. so that one resonant energy level exists, yet another level Assuming continuity, as one moves inland, this probabil- just fails to be resonant.5 Many of the fundamental con- ity will increase reaching a maximum presumably some- stants therefore seem to be boxed into a narrow range of where far in from the shoreline. arXiv:2109.10241v1 [physics.hist-ph] 10 Sep 2021 6 values compatible with our existence. This probability landscape is different from the chance How can we understand our being in such a human- of randomly selecting a universe. Because the range of compatible universe? It has been suggested7 that the parameters consistent with human life is quite small, one fundamental constants may have been selected ‘ran- might expect any smooth measure for randomly select- domly’ among all possible values. If that were the case, ing universes, to be approximately uniform across the then such compatibility is merely a condition consistent island. Sciama’s argument now follows from a well- with our being here to observe it. The conditional proba- known concentration-of-measure phenomenon.9 For any bility for such a fact would be one. Equivalently, if these uniform-density shape in a high-dimensional space, the ‘random’ selections were individual universes within a weight is almost entirely concentrated within a thin layer multiverse, then our universe being human-compatible at the surface. Figuratively, a high-dimensional elephant

2

is essentially all skin. Applied to the high-dimensional is- narrow shore, with a low chance for life like us. land of human-compatible universes, a randomly selected Notwithstanding this, where on the island the universe universe will then almost certainly be found in a narrow appears to lie can depend on the island’s shape. This band on the shore, where the probability for life would can be the case whenever our knowledge of the list of be expected to be low, since it vanishes at the shoreline. fundamental constants is incomplete. In this case, we This prediction is in contrast to that of intelligent de- would consider the island and its surrounding sea to be a sign where one might expect a universe further inland lower-dimensional space than it actually is. Our view of closer to, or possibly achieving, the greatest chance for the island would be one that projects out the unknown life. constants. This may be visualized as an ‘X-ray’ of the Is this space of universes really high dimensional? In actual n-dimensional island onto a lower m-dimensional 1936 Eddington counted four fundamental constants.10 island describing the known constants. See Fig. 1. This count excludes Newton’s gravitational constant, the In the case that the island has the shape of a uniform- speed of light, Planck’s constant and the permittivity of weight n-dimensional cube (a hypercube), with indepen- the vacuum, all used to provide scales for dimensional dent bounds on each constant, the X-ray is simply lower- quantities like length, time, mass and electric charge.10 dimensional uniform-weight hypercube. See Fig. 1(a). Just a few years ago this count had grown to 26 for the Again the lower-dimensional shore contains the greatest ‘standard model’ including the cosmological constant for weight. However, for a uniform-weight hyperball shaped gravity.11 Today, if we add three neutrino masses, the island, the X-ray, integrating out many dimensions, leads count would be 29. However, our current model of the to a narrow Gaussian with the weight concentrated at the universe hardly explains everything. There remain nu- center of the island, far inland from the shoreline (see merous long-standing open questions, many cosmological Appendix). See Fig. 1(b). in nature, such as matter-antimatter asymmetry, dark Surprisingly, the result we find for a hyperball-shaped matter, dark energy and more. Thus it would be surpris- island shown in Fig. 1(b) may well be the generic result. ing if the total number of the fundamental constants in Indeed, the coincidences found among the fundamen- a complete theory of the universe were not much larger. tal constants when determining the human-compatible (a) (b) island’s shoreline,5,6,14 suggest that sharp corners or edges are unlikely, leaving a smooth and potentially convex island. Further, the projective central limit For n constants, For n constants, 푛 12,13 1 − 0.9 of the 1 − 0.9 푛 of the theorem ensures that the projection of any such high points are within points are within n-dimensional uniform-weight shape will be well approx- 5% of the boundary. 5% of the boundary. imated by a Gaussian with variance scaling as 1/n.

X-ray of the (upper) Probability to be within 5% of the boundary X-ray of the (upper) hyperball; points hypercube; points now concentrate at are again uniformly the center. distributed

FIG. 1: Random sampling from the human-compatible island = 100 = 250 can look different dependent on its shape, if one has only lim- 𝑛𝑛 ited access to the fundamental constants. The high-weight 𝑛𝑛 = 42 region (typically the shoreline) is shown in white, with the remaining low-weight contribution in gray. The island of the 𝑛𝑛 accessible fundamental constants is obtained by integrating out those constants which are unknown or unobserved; this is visualized as an ‘X-ray’ of the actual island. (a) For an island which is an n-dimensional hypercube (upper), an X-ray re- 𝑚𝑚 duces to a uniform-measure hypercube in a lower dimension FIG. 2: Probability for the m accessible constants of a random𝑛𝑛 (lower). (b) For an island which is a uniformly-distributed universe to be within 5% of the life-denying boundary, versus n-dimensional ball (upper) with many unknown constants, the fraction m/n. The calculation assumes the island has the its X-ray is well approximated by a narrow Gaussian concen- shape of a hyperball with n fundamental constants, though trated at the center of the human-compatible island (lower). only m are known about. We consider n 42, 100, 250 . ∈ { } Although it did not figure into Sciama’s original argu- In each case, unless at least 80% of the total number of con- stants are accessible (m/n 0.8), the chance of being near ment, we shall see that the shape of the island plays a the boundary is less than ≥0.35. crucial role in the possible apparent reversal of Sciama’s ' conclusion. Note that the ‘elephant skin’ result itself is However, the projective central limit theorem only tells essentially independent of this shape, which follows sim- us that the distribution is peaked far inland from the ac- ply from the scaling of the ‘hypervolume’ with dimen- cessible shoreline. What about the probability for the ac- sionality. Thus, there is no question about a randomly cessible parameters nevertheless being found on the pro- selected universe compatible with human life having a set jected shore near the boundary? In Fig. 2 we compute of fundamental constants that almost certainly lie on the the probability to be within 5% of the projected shoreline arXiv:2109.09826v1 [hep-ph] 20 Sep 2021 a ec neulbimsaeo pitn n uin hc scle the called is which fusion, and splitting of state equilibrium th density, an gluon nonline reach the may JIMWLK/BK of value the critical on a at based that research predicted equations QCD The unitarity. fast too the be of cannot gluons of increase this However, gluo splitting. of gluon number the general, In energy. high at g processes since hadronic important, are nucleon in distributions gluon of knowledge The Introduction 1 lnighdo oldr eodteLHC the beyond colliders hadron planning keywords ann:Temn am-a-ussin gamma-ray-bursts mini The Warning: eateto hsc,Es hn omlUiest,Saga 2 Shanghai University, Normal China East Physics, of Department nlsso omcgmaryseta efidta the that t find condensation we gluon spectra, gamma-ray the cosmic exceed of analyses proton-proto energies the collision increase the that greatly will condensation gluon u oteguncnesto ffc.Sc rica iigamm mini detectors. artificial the Such spac damage effect. narrow may condensation a laboratory gluon in the emitted to be due increasing will colliders condensa gamma-rays gluon hadron strong the of extremely of generation region next energy the the for to that close are LHC the at losmycneg oasal tt taciia oetmi momentum critical a at state stable a to converge may Gluons lo odnain am-a-uss H;Le;EIC LHeC; LHC; Gamma-ray-bursts; condensation; Gluon : e h,ZiiCiadJahn Ruan Jianhong and Cui Zhiyi Zhu, Wei energies Abstract 1 p − h olso nris the energies, collision the b P rs eto provided section cross n rsod ae nthe on Based hreshold. and ineet ewarn We effect. tion si ulo rw with grows nucleon in ns fteaccelerator the of e -a-ussi the in a-ray-bursts oo-hre gluons color-charged e b P u oterestriction the to due ulo.This nucleon. n − un oiaethe dominate luons 04,China 00241, b P collisions revolution ar oo Glass Color Dynamics and Observational Signatures of Shell-like Black Hole Mimickers

Ulf Danielsson,1 Luis Lehner,2 and Frans Pretorius3 1Uppsala University, Uppsala, Sweden 2Perimeter Institute for Theoretical Physics, 31 Caroline St., Waterloo, ON, N2L 2Y5, Canada 3Princeton University, NJ, USA (Dated: September 22, 2021) We undertake the task of studying the non-linear dynamics of quantum gravity motivated alterna- tives to black holes that in the classical limit appear as ultra-compact shells of matter. We develop a formalism that should be amenable to numerical solution in generic situations. For a concrete model we focus on the spherically symmetric AdS black bubble — a shell of matter at the Buchdahl radius separating a Schwarzschild exterior from an AdS interior. We construct a numerical code to study the radial dynamics of and accretion onto AdS black bubbles, with exterior matter provided by scalar fields. In doing so we develop numerical methods that could be extended to future studies beyond spherical symmetry. Regarding AdS black bubbles in particular, we find that the original prescription for the internal matter fluxes needed to stabilize the black bubble is inadequate in dynamical settings, and we propose a two parameter generalization of the flux model to fix this. To allow for more efficient surveys of parameter space, we develop a simpler numerical model adapted to spherically symmetric bubble dynamics. We identify regions of parameter space that do allow for stable black bubbles, and moreover allow control to a desired end-state after an accretion episode. Based on these results, and evolution of scalar fields on black bubble backgrounds, we speculate on some observational consequences if what are currently presumed to be black holes in the universe were actually black bubbles.

CONTENTS B. Numerical results from the Einstein-Klein-Gordon-Hydrodynamic I. Introduction 1 system 16 1. Accuracy and convergence 16 II. AdS black bubbles and General considerations 3 2. Interior energy 17 A. Formulation 4 3. Gradual release of internal energy 18 B. Scalar field and shell material 4 C. Matter equations of motion 5 VII. Discussion 19

III. Restricting to Spherical Symmetry 5 Acknowledgments 22 A. Einstein-Klein-Gordon-Hydrodynamic system 6 A. Dynamical effects on proper acceleration 22 B. Evolution equations 6 C. Constraint Equations and Initial Data 7 B. Linear perturbation analysis 22 D. Boundary conditions 8 a. Particular solution 23 E. Simplified Shell Dynamics with an External b. Impulse response 24 Source 8 F. Simplified dissipation 9 C. x(r,t) map 24

IV. AdS Black Bubble Matter 10 D. Weak-form integration stencil 24 arXiv:2109.09814v1 [gr-qc] 20 Sep 2021 A. Specific flux model 11 1. Internal energy exchange and stability 11 E. Rate of Energy Loss from AdS interior 25 2. Alternative flux model 12 References 26 V. Implementation specifics 12 A. Einstein-Klein-Gordon-Hydrodynamic system 12 I. INTRODUCTION 1. Weak form integration 13 2. Mapping between radial metric and code Since the recognition of black hole entropy and for- coordinates 13 mulation of the black hole information paradox, many 3. Initial data 14 efforts have focused on interpreting and reconciling such B. Simplified Shell model 14 puzzling aspects of black holes. These, and related at- tempts to tame the singularities inside black holes, moti- VI. Applications/dynamics 14 vated consideration of “extensions” to black holes, alter- A. Numerical results from the shell model 14 ing their structure in the vicinity of the classical horizon Prepared for submission to JHEP PREPRINT

More on the SW-QNM correspondence

Massimo Bianchi,a Dario Consoli,a Alfredo Grillo,a Jos`eFrancisco Morales.a aDipartimento di Fisica, Universit`adi Roma “Tor Vergata” & Sezione INFN Roma2, Via della ricerca scientifica 1, 00133, Roma, Italy

Abstract: We exploit the recently proposed correspondence between gravitational per- turbations and quantum Seiberg-Witten curves to compute the spectrum of quasi-normal modes of asymptotically flat Kerr Newman black holes and establish detailed gauge/gravity dictionaries for a large class of black holes, D-branes and fuzzballs in diverse dimensions.

QNM frequencies obtained from the quantum periods of SU(2) N = 2 SYM with Nf = 3 flavours are compared against numerical results, WKB (eikonal) approximation and geode- tic motion showing remarkable agreement. Starting from the master example relating

quasi-normal modes of Kerr-Newman black holes in AdS4 to SU(2) gauge theory with Nf = 4, we illustrate the procedure for some simple toy-models that allow analytic so- lutions. We also argue that the AGT version of the gauge/gravity correspondence may give precious hints as to the physical/geometric origin of the quasi-normal modes/Seiberg- Witten connection and further elucidate interesting properties (such as tidal Love numbers and grey-body factors) that can help discriminating black holes from fuzzballs. arXiv:2109.09804v1 [hep-th] 20 Sep 2021 Mergers of Maximally Charged Primordial Black Holes

Konstantinos Kritos1, ∗ and Joseph Silk2, 1, 3, † 1Department of Physics & Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA 2Institut d’Astrophysique de Paris, UMR 7095 CNRS & UPMC, Sorbonne Universite´, F-75014 Paris, France 3Beecroft Institute for Particle Astrophysics and Cosmology, Department of Physics, University of Oxford, Oxford OX1 3RH, UK (Dated: September 22, 2021) Near-extremal primordial black holes stable over cosmological timescales may constitute a signif- icant fraction of the dark matter. Due to their charge the coalescence rate of such black holes is enhanced inside clusters and the non-extremal merger remnants are prone to Hawking evaporation. We demonstrate that if these clusters of near-extremal holes contain a sufficient number of mem- bers to survive up to low redshift, the hard photons from continued evaporation begin to dominate the high energy diffuse background. We find that the diffuse photon flux can be observed for a monochromatic mass spectrum of holes lighter than about 1013g. We place upper bounds on their abundance respecting the current bounds set by gamma ray telescopes. Furthermore, the gravita- tional wave background induced at the epoch of primordial black hole formation may be detectable by future planned and proposed ground-based and space-borne gravitational wave observatories op- erating in the mHz to kHz frequency range and can be an important tool for studying light charged primordial black holes over masses in the range 1013g − 1019g.

I. INTRODUCTION mass black holes formed via some astrophysical mecha- nism such as from core collapse of massive stars are pre- The nature of dark matter (DM) is an outstanding sumed to be electrically neutral. Depending on the initial mystery in modern cosmology. One possible explanation conditions in the early universe and the density perturba- is that DM is composed of massive compact halo objects, tions that gave rise to primordial black holes, it is possi- and black holes serve as one of the candidates. These are ble that PBHs formed with some residual electric charge. assumed to form in the early universe, and are more ap- Since the total charge in the density field averaged over propriately termed primordial black holes (PBHs). Ac- some comoving volume is expected to be zero, there is cording to the standard scenario, PBHs form during the correspondingly an equal number of positively and nega- radiation-dominated era by the gravitational collapse of tively charged primordial black holes per unit comoving large over-densities in the density field shortly after the volume. More generally, depending on the power spec- Big Bang and have been hypothesized to account for the trum of curvature perturbations that is responsible for DM in the universe [1–5]. Interest in primordial black PBH formation, PBHs may have any mass within a range holes was revived after the first gravitational wave events spanning more than twenty orders of magnitude. observed by the LIGO and Virgo collaborations [6] were The abundance of PBHs is constrained by various ob- found to involve higher mass black holes than routinely servations, for instance gravitational microlensing, cos- observed in X-ray binaries. Data analysis of these grav- mic microwave background radiation and dynamical con- itational waves indicate that they correspond to the co- siderations, and it is generally concluded that PBHs do alescence of stellar-mass black hole binaries. Some or not comprise the totality of DM in certain mass inter- even all of them could be of primordial origin and might vals [15]. For light neutral PBHs with masses below explain the merger rate of the detected events [7–13]. ∼ 1017g, the most stringent constraints come from evap- Black holes are thought to be the simplest objects oration arguments [16–20]. Black holes evaporate as

arXiv:2109.09769v1 [gr-qc] 20 Sep 2021 in the Universe in the sense that they can be fully de- thermal blackbodies [21] and those with a mass smaller scribed in terms of only three parameters: mass, spin than ' 5 × 1014g are expected to have evaporated on angular momentum and electric charge, a proposition a timescale smaller than the age of the universe [22]. that is often referred to as the no-hair theorem. The If primordial black holes are cosmologically stable, then most general geometry describing such a black hole is they can avoid evaporation and such a constraint can be the Kerr-Newman metric which reduces to the Reissner- evaded. Nordstr¨ommetric for zero spin and the Kerr metric for a An extremal black hole holds the largest possible neutral rotating black hole. So far, most interest in PBHs charge, has a vanishing Hawking temperature and does has been devoted to the case of uncharged black holes. not evaporate. But extremal black holes promptly de- Nevertheless, in principle black holes are physically al- viate from the extremality condition due to Schwinger lowed to carry an electric charge [14], even though stellar discharge [23] which is a form of athermal emission by spontaneous electron-positron pair production driven by the large electric field near the event horizon. To consider ∗ [email protected] extremal black holes that are stable against Schwinger † [email protected]; [email protected]; [email protected] discharge, we assume the black holes to be charged un- Friendship in the Axiverse: Late-time direct and astrophysical signatures of early-time nonlinear axion dynamics

David Cyncynates,∗ Tudor Giurgica-Tiron,† Olivier Simon,‡ and Jedidiah O. Thompson§ Stanford Institute for Theoretical Physics Stanford University 382 Via Pueblo, Stanford, CA 94305, USA

A generic low-energy prediction of string theory is the existence of a large collection of axions, commonly known as a string axiverse. Axions also have a natural cosmological production mech- anism, vacuum misalignment, making them well-motivated dark matter (DM) candidates. Much work on axion production has considered the case of a single free axion, but in a realistic axiverse, string axions are expected to be distributed densely over many orders of magnitude in mass, and to interact with one another through their joint potential. In this paper, we show that non-linearities in this potential lead to a new type of resonant energy transfer between axions with nearby masses. This resonance generically transfers energy from axions with larger decay constants to those with smaller decay constants, and leads to a multitude of signatures. These include enhanced direct detection prospects for a resonant pair comprising even a small subcomponent of dark matter, and boosted small-scale structure if the pair is the majority of DM. Near-future iterations of experiments such as ADMX and DM Radio will be sensitive to this scenario, as will astrophysical probes of DM substructure.

CONTENTS 3. Stability of autoresonance 27 4. Chaotic parameter space 28 I. Introduction1 C. Perturbations in detail 29 II. Friendly zero-mode dynamics4 1. General formalism 29 2. Analytical approximations 29 III. Spatial fluctuations9 A. Invitation: A single axion model of D. Simulations of non-perturbative structure growth perturbation growth9 during radiation domination 31 B. Perturbation growth during autoresonance 11 1. Metric perturbations and the equations of C. Nonperturbative structures during motion 31 autoresonance 12 2. Numerical methods 32 D. Newtonian evolution and gravitational 3. Future directions 33 collapse 14 References 33 IV. Signatures 16 A. Enhanced direct detection prospects 16 I. INTRODUCTION B. Gravitational signatures of substructure 18 C. Superradiance signatures and constraints 19 Among the best-motivated extensions of the Standard V. Repulsive self-interactions 20 Model (SM) are axions, periodic pseudoscalar fields with an approximate shift symmetry that protects their mass VI. Discussion and future directions 22 from radiative corrections.1 The most well-known exam- arXiv:2109.09755v1 [hep-ph] 20 Sep 2021 ple is the QCD axion, which was originally proposed as a Acknowledgments 24 solution to the puzzling smallness of the neutron’s elec- tric dipole moment [1–4]. This is not the only type of A. Mass vs. interaction basis 24 axion that can exist though: axions can be quite generic in UV completions of the SM with compact extra dimen- B. A more detailed study of autoresonance 25 sions and nontrivial topologies, the principle example of 1. Adiabatic evolution of resonance curves 25 which is string theory [5–7]. The combined motivation 2. Expected relic density ratio 26 of the QCD axion and string theory lead to predictions of a plenitude of string axions with mass scales spanning many orders of magnitude, a possibility referred to as the String Axiverse [8]. ∗ [email protected] † [email protected] ‡ [email protected] § [email protected] 1 The term “axion-like particles” is also used in the literature. MITP-21-041

Charged B Mesogenesis

Fatemeh Elahi,1, ∗ Gilly Elor,1, † and Robert McGehee2, ‡ 1PRISMA+ Cluster of Excellence & Mainz Institute for Theoretical Physics Johannes Gutenberg University, 55099 Mainz, Germany 2Leinweber Center for Theoretical Physics, Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA We leverage the CP violation in charged B meson decays to generate the observed baryon asymme- + try and dark matter at O(10 MeV) temperatures. We realize this in two scenarios: Bc Mesogenesis + ± ± and B Mesogenesis. In the first, CP violating Bc decays to B mesons are followed by decays to dark and Standard Model baryons. In the second, CP violating B± decays to lighter charged mesons are accompanied by the latter’s decays to dark and Standard Model leptons, which then + scatter into the baryon asymmetry. Bc Mesogenesis is actively being probed at Belle and LHCb, while B+ Mesogenesis can be tested at colliders and sterile neutrino searches.

I. INTRODUCTION paradigm of low-scale baryogenesis has been proposed: Mesogenesis [30, 31]. In this framework, an out-of- How did we come to be here? Said quantitatively: equilibrium scalar decays to SM quarks which hadronize what are the origins of the measured baryon asymmetry at low temperatures. The resulting SM mesons undergo (BAU) and dark matter? The answer to this fundamental known CP-violating processes and decay into dark sec- question still eludes us after decades of effort. tor particles carrying SM baryon or lepton number. The Explanations of the BAU usually fall into one of two decays conserve baryon and lepton number and thus gen- broad categories: electroweak baryogenesis [1–11] or lep- erate an equal and opposite baryon or lepton asymmetry togenesis [12]. These ideas satisfy the three Sakharov between the dark and visible sectors. In the latter sce- conditions [13] – baryon number violation, C and CP nario, dark sector interactions then convert the lepton Violation (CPV), and departure from thermal equilib- asymmetry into a baryon asymmetry. rium – in unique ways. Electroweak baryogenesis at- What is most compelling about Mesogenesis is that tempts to explain the BAU and satisfy these criteria it revives an original, but long dead [32–35], hope of using a strongly first order electroweak phase transition electroweak baryogenesis: that the requisite CP-violating (EWPT), while leptogenesis uses out-of-equilibrium de- processes already reside in the SM. In Mesogenesis, we cays of heavy neutrinos already motivated by the seesaw look to the mesons. The BAU is directly proportional to mechanism [14–16]. CP-violating experimental observables, making Mesoge- However, each of these answers suffers significant draw- nesis testable at current experiments. Furthermore, the backs. Many models of electroweak baryogenesis require dark sector typically contains a dark matter candidate fine tuning [17]; construct extended Higgs sectors, but whose abundance will be generated along with the baryon still can’t make the EWPT strongly first order [18, 19]; asymmetry. All of these mechanisms do not violate B or fail to actually produce the observed baryon asymme- L, but rather, hide equal and opposite asymmetries in try [20]; or are outright excluded by increasingly precise the dark sector. experimental results [21]. Often, they simultaneously ne- In neutral B Mesogenesis [30], the baryon asymmetry 1 0 glect dark matter. Though the original formulations of is generated by leveraging the CPV in Bq particle/anti- electroweak baryogenesis were minimal, Nature increas- particle oscillations, while in D+ Mesogenesis [31], the ingly seems to disfavor this now less-than-simple asym- CPV of D± meson decays is utilized. Both these fla- metry generator, perhaps in favor of a similar mechanism vors of Mesogenesis necessarily occur at low (5-20 MeV) in the dark sector [23–28]. Leptogenesis, in contrast, suf- scales, and are generically testable.2 Indeed for the case arXiv:2109.09751v1 [hep-ph] 20 Sep 2021 fers not from experimental exclusion but from exclusion of neutral B Mesogenesis, some experimental searches of experiment – there is no hope of directly probing the and proposals are already in various stages at LHCb heavy states in leptogenesis models [29]. It may be the [38, 39] and Belle-I and -II [40]. true origin of the BAU, but humans may never know. But some of the simplest and most compelling Mesoge- These substantial disadvantages of decades-old ideas nesis stories have, until now, been overlooked. Charged should not be ignored. They portend eventual failure and B mesons contain a large amount of CPV in their decays sound a call to innovation. Answering this call, a new to SM final states [41]. In the present work, we introduce two scenarios of Charged B Mesogenesis, in which baryo-

∗ [email protected] † [email protected] ‡ [email protected] 2 See [36] for a detailed study of experimental implications of de- 1 Sometimes, dark matter may be explained “after the fact” (see cays into dark sector baryons and [37] for additional indirect e.g. [22]). signals. arXiv:2109.10334v1 [astro-ph.GA] 21 Sep 2021 mi:[email protected] Email: orsodn uhr .K Bhadari K. N. author: Corresponding urudn niomn ssila pnrsac topic research their open with astrophysics. an stars in still massive is of environment interaction surrounding the of study the seterve ril by ways numerous article in review formation the star (see trigger may and size, ewe h oie einadtesronigneu- surrounding H the an and gas, region tral/molecular ionized the between uhmsiesasinz h yrgnaosadcre- and atoms hydrogen H the ate ionize stars massive such clsaeo h neselrmdu IM.Teex- The (ISM). medium phys- (EUV; interstellar the ultraviolet control the treme to of them state allow ical which energetics, ful .K Bhadari, K. N. 4 yee sn L using 2021 Typeset 22, September version Draft asv Btp tr (M stars OB-type Massive eateto srnm n srpyis aaIsiueo Institute Tata Astrophysics, and Astronomy of Department ii hls hc ih aeoiiae u ogaiainlisaiiyi th in Keywords: instability gravitational to due originated have might which shells, oeua oe nS0 r enajcn ooeaohraon t around (mass another clumps one dust to respe adjacent and geometry, seen molecular shell are dense neutral S305 spaced outer in and zones re inner molecular The the component. trace redshifted spatially brighter the have and blue-skewed ne eta hl.Tevlct hne asadposition-velocity and maps [C channel compact velocity a The of shell. neutral inner ne n ue eta hls aigattlms of mass total a having shells) neutral outer and inner omdsqetal.Te[C The sequentially. formed v nO.Vsa tiscnr n a opc iglk appearance. ring-like compact [C a diffuse has and and extended centre more its with at star O9.5V an iedt,wihaeue oeaietegsdnmc n structure and dynamics gas the examine to [C used integrated are The which data, line 3 ein.Det h hra rsuedifference pressure thermal the to Due regions. nttt fApidPyiso h usa cdm fScien of Academy Russian the of Physics Applied of Institute exp erpr nosrainlsuyo h aatcH Galactic the of study observational an report We ∼ . ms km 1.3 rbn a ieaisadPRsrcueaon -yesasi S in stars O-type around structure PDR and kinematics gas Probing A T ,2 1, 5 1. E rahtaRsac nttt fOsrainlSciences, Observational of Institute Research Aryabhatta X omto tr:pemi sequence pre–main stars: – formation ut xicin–HIrgos–IM lus–IM niiulobject individual ISM: – clouds ISM: – regions HII – extinction dust, .K Dewangan, K. L. INTRODUCTION twocolumn − 2 ii 1 ninIsiueo ehooyGnhngrPlj Gandhin Palaj, Gandhinagar Technology of Institute Indian h ue hl per ob le hnteinrsel itn tha hinting shell, inner the than older be to appears shell outer The . ii hl meddi h iueotrsel n ohtesel emto seem shells the both and shell, outer diffuse the in embedded shell ] ii hν msinmpa 3.,4.]k s km 49.5] [39.4, at map emission ] lere 1998 Elmegreen 1 eincniust rwin grow to continues region hsclRsac aoaoy arnpr,Amdbd380 Ahmedabad Navrangpura, Laboratory, Research Physical ∗ ≥ tl nAASTeX62 in style 36e)poosfrom photons eV) 13.6 > M 8 ii 1 rfie r xmndtwr 35 hc r ihrdouble-peake either are which S305, toward examined are profiles ] .M Zemlyanukha, M. P. ⊙ rv power- drive ) ii .However, ). msin otn nO.Vsa tiscnr n urud the surrounds and centre its at star O8.5V an hosting emission, ] ABSTRACT udmna eerh oiBah od ubi400,In 400005, Mumbai Road, Bhabha Homi Research, Fundamental f ∼ 3 10–10 .K Ojha, K. D. e 6U’ao t. 090Nzn ogrd Russia. Novgorod, Nizhny 603950 str., Ul’yanov 46 ces − ∼ rb fknmtcadrdaieitrcino massive of interaction radiative and kinematic of probe a eg,H (e.g., gas eta nPR,btsilicue h rcr fsingly of tracers C the primarily (e.g., includes is species still gas ionized but The PDRs, limit. po- in Lyman neutral ionization the having than less elements tential those photoionize and rcro PH msin(e.g., emission (PAH) drocarbon ueln finzdcro,[C carbon, the ionized of PDRs, line in ture coolant major a 1999 1985a 1999 Hollenbach & Tielens h rniinzn ewe oie n oeua gas in molecular 1 and Figure ionized (see between zone transition the sr ftegs r eerdt sphotodissocia- as to e.g., referred (PDRs; chem- regions are and/or photon-dominated gas, or heating regions the the tion in of role istry significant a eV plays 6 (FUV; ultraviolet ii 1 6 M 565 h urudnso h H the of surroundings The einS -0/35uigte[C the using 2-305/S305 Sh region ptal rcstosellk tutrs(i.e., structures shell-like two traces spatially 3 ; molecules the dissociate can photons FUV The ). aoaPa,Niia 602 India. 263002, Nainital Peak, Manora M afa ta.2006 al. et Kaufman ⊙ ⊙ r netgtdaon h neutral the around investigated are ) h ne eta hl encompasses shell neutral inner The . sitdadbusitdcomponents blueshifted and dshifted 4 2 tvl.Teinzd eta,and neutral, ionized, The ctively. .I Zinchenko, I. I. gr325,India. 382355, agar n O,adPlcci rmtcHy- Aromatic Polycyclic and CO), and eOtp tr.Teregularly The stars. O-type he oee,teotrseli seen is shell outer the However, igascnr h presence the confirm diagrams iln olnah1985a Hollenbach & Tielens hl fcletdmaterials. collected of shell e 0,India 009, fpooiscainregions. photodissociation of -0 I region HII 2-305 h + S , ≤ .Hne h Dsspecify PDRs the Hence, ). + S -0)–stars: – 2-305) (Sh Si , 3 hν , ii b ii n arb Sharma Saurabh and hs hlsare shells these t ; 158 ] 2 + ein,weetefar the where regions, P n Fe and , olnah&Tielens & Hollenbach Tielens & Hollenbach ≤ 3 xadwith expand / 36e)radiation eV) 13.6 2 µ ii – ffr unique a offers m 158 ] 2 P 1 + / 2 or d ,molecular ), µ n struc- fine m .Being ). dia. 5 MNRAS 000,1–10 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

Detectability of Population III stellar remnants as X-ray binaries from tidal captures in the local Universe

Rabia Husain,1,2★ Boyuan Liu1 and Volker Bromm1 1Department of Astronomy, University of Texas, TX 78712, USA 2Jefferson Laboratory of Physics, Harvard University, Cambridge, MA 02138

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We assess the feasibility of detecting the compact object remnants from Population III (Pop III) stars in nearby dense star clusters, where they become luminous again as X-ray binaries (XRBs) and tidal disruption events (TDEs) via strong tidal encounters. Analytically modelling the formation of Pop III stars, coupled with a top-heavy initial mass function predicted by numerical simulations, we derive the number of (active) Pop III XRBs and TDEs in the present-day Milky Way (MW) nuclear star cluster as ∼ 0.06 − 0.3 and . 4 × 10−6, rendering any detection unlikely. The detection probability, however, can be significantly boosted when surveying all massive star clusters from the MW and neighboring galaxy clusters. Specifically, we predict ∼ 1.5 − 6.5 and ∼ 40 − 2800 active Pop III XRBs in the MW and the Virgo cluster, respectively. Our Pop III XRBs are dominated (∼ 99%) by black holes with a typical mass and 36 −1 luminosity of ∼ 45 M and ∼ 10 erg s . Deep surveys of nearby (. 30 − 300 Mpc) galaxy clusters for such Pop III XRBs are well within reach of next-generation X-ray telescopes, such as ATHENA and LYNX. Key words: early universe - dark ages – reionization – first stars - neutron stars, black holes

1 INTRODUCTION (Rydberg et al. 2013). Rather, a futuristic telescope such as the Ultimately Large Telescope (ULT), a 100m liquid-mirror telescope Heralding the end of the cosmic dark ages, the formation of the first on the Moon, would be necessary to image these exceedingly faint stars, the so-called Population III (Pop III), signaled the beginning objects (Schauer et al. 2020). of complex structure formation in the Universe (Bromm & Larson Given the challenges to detect the first stars directly at high 2004). These first stars formed out of the primordial H and He redshifts, with any telescope capable of doing so, such as the ULT, synthesized in the Big Bang (Haiman et al. 1996; Liu & Bromm 6 lying far in the future, we must turn to indirect ways to elucidate their 2018), within dark matter minihaloes with virial masses ∼ 10 M properties, and to test our respective theories. One approach aims at 푧 ∼ 20−30 (Yoshida et al. 2006). As a consequence of this lack of at their compact, neutron star (NS) and black hole (BH), remnants, metals, their properties are different from present-day stars, causing which persist throughout cosmic history, allowing us to probe them them to be more massive, typically & 10 − 20 M , implying that in the local Universe. For instance, the gravitational waves from they end their lives in energetic supernova (SN) explosions, leaving the mergers of Pop III remnants have been intensively studied (e.g. massive remnants behind (Abel et al. 2002; Bromm et al. 2002). The Kinugawa et al. 2021b; Hijikawa et al. 2021; Liu & Bromm 2021), SN explosions in turn created metal-rich environments in the early and found to be able to explain the most massive events detected by

arXiv:2109.10321v1 [astro-ph.HE] 21 Sep 2021 intergalactic medium (IGM) from which Population I/II (Pop I/II) the Laser Interferometer Gravitational-Wave Observatory (LIGO) stars were able to form at later times (Ostriker & Gnedin 1996; campaign (Abbott et al. 2020c; Nitz et al. 2021; The LIGO Scientific Heger & Woosley 2002; Karlsson et al. 2013; Jeon et al. 2015). Collaboration et al. 2021). Cosmological hydrodynamic simulations show that Pop III In addition to gravitational wave sources, Pop III compact star formation likely ended around 푧 ∼ 6 (Johnson et al. 2013; Xu object remnants1 can also become luminous in X-rays when they et al. 2016; Sarmento et al. 2018a; Liu & Bromm 2020a). This may experience strong tidal encounters with a companion star to form be within the capabilities of next-generation telescopes, depending an X-ray binary (XRB) or a tidal disruption event (TDE). In this on the intrinsic brightness of these sources. However, even frontier study, we investigate the feasibility of detecting Pop III XRBs/TDEs telescopes that are set to launch soon, such as the James Webb Space Telescope (JWST), will not be able to detect individual minihaloes at 푧 > 10 (Zackrisson et al. 2011, 2015), due to their low luminosities 1 The Pop III initial mass function is also predicted to reach the low mass regime of white dwarf progenitors. Here, we will not consider white dwarf ★ E-mail: [email protected] remnants (but see Chabrier 2004).

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

The influence of the environment on the spin evolution of low-mass stars. I. External photoevaporation of circumstellar disks

J. Roquette,1★, S. P. Matt1, A. J. Winter2, L. Amard1, and S. Stasevic1 1University of Exeter, Physics and Astrophysics dept, Stoker Road, EX44QL, Exeter, UK 2 Institut für Theoretische Astrophysik, Zentrum für Astronomie, Heidelberg University, Albert Ueberle Str. 2, 69120 Heidelberg, Germany

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Massive stars are strong sources of far-ultraviolet radiation that can be hostile to the evolution of protoplanetary disks, driving mass loss by external photoevaporation and shortening disk-dissipation timescales. Their effect may also reduce the timescale of angular momentum exchanges between the disk and host star during the early pre-main-sequence phase. To improve our understanding of the environmental influence on the rotational history of stars, we developed a model that considers the influence of the local far-ultraviolet radiation on the spin evolution of low mass stars. Our model includes an assumption of disk-locking, which fixes the rotation rate during the star-disk-interaction phase, with the duration of this phase parametrised as a function of the local far-ultraviolet radiation and stellar mass (in the range 0.1–1.3 M ). In this way, we demonstrate how the feedback from massive stars can significantly influence the spin evolution of stars and explain the mass-dependency observed in period-mass distributions of young regions like Upper Sco and NGC 2264. The high far-ultraviolet environments of high-mass stars can skew the period distribution of surrounding stars towards fast-rotation, explaining the excess of fast-rotating stars in the open cluster h Per. The proposed link between rotation and the pre-main-sequence environment opens new avenues for interpreting the rotational distributions of young stars. For example, we suggest that stellar rotation may be used as a tracer for the primordial ultraviolet irradiation for stars up to ∼1 Gyr, which offers a potential method to connect mature planetary systems to their birth environment. Key words: stars: evolution – stars: rotation – stars: low-mass – stars:solar-type – stars: pre-main-sequence

1 INTRODUCTION clusters, depending on their proximity to massive stars (Balog et al. 2007; Guarcello et al. 2010; Fang et al. 2012; Guarcello et al. 2016; Over the last couple of decades, a growing number of studies have Ansdell et al. 2017; Eisner et al. 2018). Because of its importance been exploring the influence of the environment on the evolution on constraining the timescales for planet-formation (e.g., Armitage and dissipation of protoplanetary disks (e.g., Bally et al. 1998; John- 2000; Najita & Kenyon 2014; Concha-Ramírez et al. 2019b; Ansdell stone et al. 1998; Adams et al. 2004; Adams 2010; Facchini et al. et al. 2020), the dependence of the disk-dissipation timescales with 2016; Nicholson et al. 2019; Haworth & Clarke 2019; Sellek et al. the local FUV flux has been recently theoretically investigated by a 2020; Parker et al. 2021). Of particular importance is the existence number of authors (e.g., Adams et al. 2004; Facchini et al. 2016; Ha- of massive stars in a cluster’s population (Adams et al. 2004; Win- worth et al. 2018b; Winter et al. 2018; Concha-Ramírez et al. 2019a; ter et al. 2018). These massive stars can radically influence their Winter et al. 2020a; Sellek et al. 2020; Parker et al. 2021). surrounding environment due to their strong far-ultraviolet (FUV) arXiv:2109.10296v1 [astro-ph.SR] 21 Sep 2021 emission. FUV photons (6eV< ℎ휈 < 13.6eV) can dissociate H2 External photoevaporation of disks will also shorten the star-disk- molecules and trigger the external photoevaporation of the material interaction (SDI) phase. During this phase, the magnetic interaction in protoplanetary disks around neighbouring low mass stars, short- between stars and their accretion disks leads to an exchange of mass ening the timescales of disk evolution and dissipation (Adams et al. and angular momentum between the two (e.g., Ghosh & Lamb 1979; 2004; Winter et al. 2020a; Concha-Ramírez et al. 2021). A number Koenigl 1991; Ireland et al. 2021). Stars lose the vast majority of of studies have observationally explored the indirect consequences their angular momentum during their first ∼10 Myr, and it is thought of external photoevaporation of disks, such as the study of proplyd that the SDI is the primary driver of the stars’ angular momentum objects in the neighbourhood of massive stars (Smith et al. 2003; evolution (Bodenheimer 1995; Mathieu 2004; Bouvier et al. 2014). Balog et al. 2006; Rigliaco et al. 2009; Wright et al. 2012; Mann Physical models for the magnetic SDI show that the net result of et al. 2014; Guarcello et al. 2014; Kim et al. 2016; Haworth et al. this interaction can either spin up or spin down the star, driving 2021) and variations of the disk-fraction and disk-masses of stars in the star’s rotation toward an equilibrium spin rate (Ghosh & Lamb 1979; Matt & Pudritz 2005; Zanni & Ferreira 2009, 2013; Ireland et al. 2021). In this theoretical equilibrium state, the rotation rate of ★ E-mail: [email protected] the star is dictated by a balance of torques, which is determined by

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

Iron Rain: measuring the occurrence rate and origin of small iron meteoroids at Earth

Tristan Mills1,2★, P. G. Brown1,2, M. J. Mazur1,2, D. Vida1,2, Peter S. Gural3 and Althea V. Moorhead4 1Department of Physics and Astronomy, University of Western Ontario, London, Ontario, N6A 3K7, Canada 2 Institute for Earth and Space Exploration (IESX), The University of Western Ontario, London, Ontario N6A 3K7, Canada 3Gural Software and Analysis LLC, Lovettsville, VA USA 4NASA Meteoroid Environment Office, Marshall Space Flight Center EV44, Huntsville, Alabama 35812, USA

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We report results of a four-year survey using Electron Multiplied Charged Coupled Device (EMCCD) cameras recording 34761 two-station video meteor events complete to a limiting magnitude of +6. The survey goal was to characterize probable iron meteoroids. Using only physical properties of the meteor trajectories including early peaking light curves, short lumi- nous trajectories, and high energies accumulated per area at beginning, we identified 1068 iron meteors. Our iron candidates are most abundant at slow speeds < 15 km/s, where they make up ≈20% of the mm-sized meteoroid population. They are overwhelmingly on asteroidal orbits, and have particularly low orbital eccentricities and smaller semi-major axes when compared to non-irons between 10-20 km/s. Our iron population appears to be more numerous at fainter magnitudes, comprising 15% of slow (10-15 km/s) meteors with peak brightness of +3 with the fraction rising to 25% at +6 to +7, our survey limit. The iron orbits are most consistent with an asteroidal source and are in highly evolved orbits, suggesting long collisional lifetimes (107 years). Metal-rich chondrules (nodules) found in abundance in EL Chondrites are one possible source for this population. We also propose a possible technique using R-band colours to more robustly identify fainter iron meteors with very high confidence. Key words: meteors – meteoroids

1 INTRODUCTION dominated precursor meteoroids, recovered in Antarctica is < 2% (Taylor et al. 2000). These results suggest that iron meteoroids are The composition and origin of extraterrestrial dust impacting Earth a small fraction of the impactor population at Earth. are varied. At large sizes, the many distinct classes of meteorites However, more recently, a spectral survey performed by Vo- demonstrate directly that asteroidal parent bodies have experienced jáček et al.(2015, 2019) classified 15 out of 152 meteors as iron diverse evolutionary histories. The distinct populations of interplan- at mm-sizes. Campbell-Brown(2015) reported a population of low, etary dust particles and micrometeorites (Bradley 2007; Genge et al. slow meteors at mm and smaller meteoroid sizes that appeared to 2008) also show that many different classes of solar system bodies have properties consistent with iron, and appeared to increase in deliver dust to Earth (Carrillo-Sánchez et al. 2016) and that these

arXiv:2109.10265v1 [astro-ph.EP] 21 Sep 2021 abundance with decreasing mass. This suggests that the fraction of parent populations differ significantly. iron meteoroids as a function of size and orbit type or entry speed Among all the material impacting Earth, iron mete- may vary and raises the question: what fraction of the meteoroid oroids/meteorites are among the most processed and distinct. While population closer to the mass influx peak at Earth (near 10−8 kg comprising 4% of all observed meteorite falls (Borovička et al. Love & Brownlee(1993)) are iron meteoroids? 2019), the fraction of apparently iron meteoroids producing fire- In addition to questions of the origin of such a unique compo- balls is much smaller, of order 1% (Borovička 2006). The differ- sitional meteoroid population, ‘irons’ pose a greater risk to space- ence is likely due to recovery biases, such as the higher density of craft. Single-plate penetration equations can be used to assess how a iron meteorites, their resistance to weathering and their ability to spacecraft would withstand hypervelocity impacts, and many, such be magnetically isolated. At still smaller sizes, the abundance of as the Cours-Palais equation, have a dependence on impactor den- I-type micrometeorite spherules, believed to be derived from iron- sity (Hayashida & Robinson 1991). Due to the higher density of iron, the predicted damage caused by iron meteoroids of similar ★ E-mail:[email protected] velocity and mass is greater. Understanding this population of small

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

Imprints of the jittering jets explosion mechanism in the morphology of the supernova remnant SNR 0540-69.3 Noam Soker1

1Department of Physics, Technion, Haifa, 3200003, Israel; [email protected]

ABSTRACT I identify a point-symmetric structure in recently published velocity maps of different elements in a plane along the line of sight at the center of the supernova remnant SNR 0540-69.3, and argue that jittering jets that exploded this core collapse supernova shaped this point-symmetric structure. The four pairs of two opposite clumps that compose this point symmetric structure suggest that two to four pairs of jittering jets shaped the inner ejecta in this plane. In addition, intensity images of several spectral lines reveal a faint strip (the main jet-axis) that is part of this plane of jittering jets and its similarity to morphological features in a few other SNRs and in some planetary nebulae further suggests shaping by jets. My interpretation implies that in addition to instabilities, jets also mix elements in the ejecta of core collapse supernovae. Based on the point-symmetric structure and under the assumption that jittering jets exploded this supernova, I estimate the component of the neutron star natal kick velocity on the plane of the sky to be ' 235 km s−1, and at an angle of ' 47◦ to the direction of the main jet-axis. I analyse this natal kick direction together with other 12 SNRs in the frame of the jittering jets explosion mechanism.

Keywords: ISM: supernova remnants - stars: jets - supernovae: general - supernovae: individual (SNR 0540-69.3)

1. INTROCUTION stabilities supply the stochastic angular momentum to Core collapse supernova (CCSN) remnants (SNRs) form the intermittent accretion disk that launches the have inhomogeneous structures of filaments, arcs, jittering jets (e.g., Papish et al. 2015a; Shishkin & Soker clumps, and ‘Ears’ (two opposite protrusions from the 2021). main SNR). Examples include the SNRs Vela (images In some cases instabilities alone cannot account for by, e.g., Aschenbach et al. 1995; Garc´ıa et al. 2017), the filamentary structure of the ejecta and it seems SNR G292.0 + 1.8 (e.g., Park et al. 2002, 2007), and that jets as well shape the ejecta. Consider the SNR W49B (e.g., Lopez et al. 2013; Sano et al. 2021). clumpy/filamentary structure of the ejecta of SN 1987A This holds for the inner ejecta of CCSNe that have inho- (e.g., Fransson et al. 2015, 2016; Larsson et al. 2016; mogeneous structures of filaments, clumps, and rings of Abell´an et al. 2017; Matsuura et al. 2017), now the different heavy elements such as oxygen, silicon, sul- SNR 1987A. Although there are claims that the non- fur, argon and iron. According to the delayed neutrino symmetric explosion of SN 1987A is due to instabili- explosion mechanism of CCSNe instabilities that are in- ties alone (e.g., Kjær et al. 2010), recent studies suggest herently exist in this explosion mechanism cause these that this is not the case. Abell´anet al.(2017) find that arXiv:2109.10230v1 [astro-ph.HE] 21 Sep 2021 filamentary structures of the inner ejecta (e.g., Janka none of the neutrino driven explosion models they com- et al. 2017; Wongwathanarat et al. 2017; Gabler et al. pare to fit all observations of SN 1987A, a conclusion 2021; Sandoval et al. 2021). According to the jittering that supports a similar earlier claim by Soker(2017). jets explosion mechanism of CCSNe both instabilities Soker(2017) compared the Fe structure of SN 1987A and jittering jets shape the ejecta (e.g., Papish & Soker from Larsson et al.(2016) with the numerical simula- 2014a; Papish et al. 2015a; Gilkis & Soker 2016). In tions by Wongwathanarat et al.(2015) and concluded the jittering jets explosion mechanism instabilities in- that the neutrino-driven explosion mechanism cannot herently exist, in particular the spiral standing accre- account for the properties of SN 1987A. Bear & Soker tion shock instability (for studies of this instability see, (2018b) use the observations of Abell´anet al.(2017) to e.g., Blondin & Mezzacappa 2007; Rantsiou et al. 2011; compare some morphological features of SN 1987A with Fern´andez 2015; Kazeroni et al. 2017), because these in- morphological features of other SNRs and with plane- arXiv:2109.10196v1 [astro-ph.SR] 21 Sep 2021 ytm a lola otdlprubto rectto fteros their of excitation or perturbation tidal to lead also may systems sators r spsil.Ti a o enahee o aytpso pulsat of types many s for on achieved constraints been many now δ has as This wield possible. to as order ory in stars, pulsating containing fselreouini re ogieterfieeto hs models these of refinement the guide to the order of in predictions evolution the to stellar compared of be can properties their as value rprissc sterdniis gsadrttoa profiles rotational and ages mod densities, theoretical their to as compared such be star may properties pulsating objects the these is in structure, frequencies interior their particular in stars normal of radii and masses r ieyt rs ntescnaysa.R y soeo h f Introduction the of one is pulsations. Lyn of RR propert types well-established star. multiple have exhibit secondary to the p known in these binaries Ly Th RR of eclipsing arise all perturbed). in to or likely pulsations tidally Most are excited origin. be instrumental tidally may of be alternatively which may of frequencies lower (some pulsations type a opsto n naeo y,adadsac lgtysh slightly distance a the and from Gyr, ch 1 solar expected of a age for than an models and theoretical composition with agreement cal good b a upwards find temperatures we effective published adjusting After yeplain,adteitreit rqece to frequencies intermediate the and pulsations, type RLnsoscereiec o ustos emaue3 pul to 35 frequencies measure higher We the pulsations. attribute for and evidence frequencies clear sation shows Lyn RR s h ih uv rmthe from curve metal- light possibly the one and use metallic-lined one A-stars: two ing h ytmt ihpeiin efidmse f1 of masses find We precision. prope high the to determine system to the measurements velocity radial lished 01Dec 2021 1 . ct stars Scuti novosga st obn hs w ye faayi ystudying by analysis of types two these combine to is goal obvious An nte yeo betwl utdt rbn h hsclproper physical the probing to suited well object of type Another measuremen direct of source primary our are stars binary Eclipsing RLni eahdelpigbnr iha99 ri conta orbit d 9.95 a with binary eclipsing detached a is Lyn RR 510 ET CT,GMADRDSADTIDALLY-PERTURBED AND DORADUS GAMMA SCUTI, DELTA srpyisGop el nvriy tffrsie T 5 ST5 Staffordshire, University, Keele Group, Astrophysics 17–19 ± EICSINO CISN IAIS AE VII. PAPER BINARIES. ECLIPSING OF REDISCUSSION 0 . n e inswt oa-ieoscillations solar-like with giants red and 0 M 003 10–13 , γ ⊙ oau stars Doradus n ai f2 of radii and , USTOSI RLYNCIS RR IN PULSATIONS yJh Southworth John By Gaia .Southworth J. 1,2 14,15 TESS eahdelpigbnre dB)aeof are (dEBs) binaries eclipsing Detached . D3prla.Telgtcreof curve light The parallax. EDR3 lwyplaigB-stars slowly-pulsating , . 564 ± aelt n w eso pub- of sets two and satellite 0 . 1 n 1 and 019 20–22 h iaiyo these of binarity The . 5 eetdoscillation Detected . . . 939 613 6–9 16 . l oconstrain to els rtclmodels oretical , G UK BG, ± ± cillations β γ iso stars, of ties rincluding or 0 ehipul- Cephei 3,4 e n to and ies 0 Doradus- or We poor. ulsations . . elrthe- tellar 0 and 007 . 1 R 013 0 K 200 y δ so the of ts te of rties Scuti- orter dEBs emi- or n 23–26 ew in- ⊙ e 1 - . . Draft version September 22, 2021 Typeset using LATEX twocolumn style in AASTeX62

The phantom dark matter halos of the Local Volume in the context of modified Newtonian dynamics

P.-A. Oria,1 B. Famaey,1 G. F. Thomas,2, 3 R. Ibata,1 J. Freundlich,1 L. Posti,1 M. Korsaga,1, 4 G. Monari,1 O. Muller,¨ 1 N. I. Libeskind,5, 6 and M. S. Pawlowski5

1Universit´ede Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, F-67000 Strasbourg, France 2Instituto de Astrof´ısica de Canarias, E-38205 La Laguna, Tenerife, Spain 3Universidad de La Laguna, Dpto. Astrof´ısica, E-38206 La Laguna, Tenerife, Spain 4Laboratoire de Physique et de Chimie de l’Environnement, Observatoire d’Astrophysique de l’Universit´eJoseph Ki-Zerbo (ODAUO), 03 BP 7021,Ouaga 03, Burkina Faso 5Leibniz-Institut fur Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany 6Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, IP2I Lyon / IN2P3, IMR 5822, F-69622, Villeurbanne, France

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

ABSTRACT We explore the predictions of Milgromian gravity (MOND) in the Local Universe by considering the distribution of the ‘phantom’ dark matter (PDM) that would source the MOND gravitational field in Newtonian gravity, allowing an easy comparison with the dark matter framework. For this, we specifically deal with the quasi-linear version of MOND (QUMOND). We compute the ‘stellar-to- (phantom)halo-mass relation’ (SHMR), a monotonically increasing power-law resembling the SHMR observationally deduced from spiral galaxy rotation curves in the Newtonian context. We show that the gas-to-(phantom)halo-mass relation is flat. We generate a map of the Local Volume in QUMOND, highlighting the important influence of distant galaxy clusters, in particular Virgo. This allows us to explore the scatter of the SHMR and the average density of PDM around galaxies in the Local Volume, Ωpdm ≈ 0.1, below the average cold dark matter density in a ΛCDM Universe. We provide a model of the Milky Way in its external field in the MOND context, which we compare to an observational estimate of the escape velocity curve. Finally, we highlight the peculiar features related to the external field effect in the form of negative PDM density zones in the outskirts of each galaxy, and test a new analytic formula for computing galaxy rotation curves in the presence of an external field in QUMOND. While we show that the negative PDM density zones would be difficult to detect dynamically, we quantify the weak lensing signal they could produce for lenses at z ∼ 0.3.

Keywords: –

1. INTRODUCTION discs (e.g. Posti et al. 2019a; Marasco et al. 2020), the The missing mass problem is one of the most press- planes of satellite galaxies problem (Pawlowski 2018), or ing questions in present day (astro)physics. The cur- the prevalence of cold stellar kinematics and the absence of bulges (or massive stellar halos) in most disc galaxies

arXiv:2109.10160v1 [astro-ph.GA] 21 Sep 2021 rent dominant Λ Cold Dark Matter (ΛCDM) paradigm – while reproducing with impressive accuracy a large (e.g., Peebles 2020). This justifies exploring alternative number of observations on all scales – still has some frameworks, which can range from modifications of the pending issues, both at cosmological scales (e.g., Riess dark matter properties to radical modifications of grav- et al. 2019; Bowman et al. 2018) and on small scales ity. (e.g., Bullock & Boylan-Kolchin 2017). This includes In galaxies – especially rotationally-supported ones challenges in explaining the diversity of galaxy rotation – the observed dynamics can be predicted surprisingly curve shapes (e.g. Oman et al. 2015), the suprisingly low well based on the distribution of baryons alone, through scatter of the Baryonic Tully-Fisher Relation (BTFR, Milgrom’s law (Milgrom 1983) which is at the heart of McGaugh et al. 2000; Lelli et al. 2016a; Desmond 2017a; the Modified Newtonian Dynamics (MOND) paradigm Lelli et al. 2019), the high baryon fraction in massive (Milgrom 1983; Famaey & McGaugh 2012). Milgrom’s arXiv:2109.10151v1 [astro-ph.IM] 21 Sep 2021 [email protected] orsodn uhr ila Cotton William author: Corresponding enrttso h k ihprlatcangle parallactic pat- with antenna the sky as the axis on the rotates in off mountedtern variations sources individ- alt-az of gain the direction axis dependent the of 2 time in pattern introduce Asymmetries beam antennas the antennas. by ual modified as sky fwae ore.I swrhntn htacrt cor- accurate that noting worth is It sources. imaging weaker the corrupt of differ- seriously and from can asymmetries antennas Artifacts among beam ences time. by same caused the sources differ- at bright a even see even sky antennas effective an different ent the are an- sky, as case of the extreme arrays on more patterns arrays, dissimilar heterogeneous with of tennas use and The derived source the varying image. in artifacts time non-convolutional a introduce will observing of equivalent the 1 ai nefrmti raspoueiae fthe of images produce arrays interferometric Radio yee sn L using 2021 Typeset 22, September version Draft SA a xsmutt lmnt hseffect this eliminate to mount axis 3 a ( has ASKAP ontne l 2008 al. et Johnston Keywords: tim described. run test computational realistic additional a dynam also in The high are 50% sources. a properties in polarized polarization sources partly Source the and recovered. ante of fr well m indices those 15 are spectral of and model level and 13.5 of the densities array near flux mixed The to a a U artifacts of Q, simulations of use I, level demonstrations parameters the Stokes reduce used all to in antennas data images among correct differences can technique and The patterns antenna in metries edsrb n eosrt ehiu o orcigdirection correcting for technique a demonstrate and describe We A T 1. E X srnmcltcnqe,Itreoer,Polarimetry Interferometry, techniques, Astronomical INTRODUCTION orcino ai nefrmti mgn o nen Pa Antenna for Imaging Interferometric Radio of Correction twocolumn ). 2 ot fia ai srnm bevtr,2FrS. Obser St., Fir 2 Observatory, Astronomy Radio African South Rcie ? eie ? cetd? rf etme 2 20 22, September draft ?? Accepted ??; Revised ??; (Received tl nAASTeX63 in style 1 ainlRdoAtooyObservatory Astronomy Radio National .D Cotton D. W. hrotsil,V 20,USA 22903, VA Charlottesville, 1 hsis This . umte oPASP to Submitted 2 deotRoad Edgemont 520 ABSTRACT 1 n .Mauch T. and a erA+(35mMeKTdse ls1 SKA and m 15 patterns plus dishes MeerKAT antenna ar- m heterogeneous (13.5 MeerKAT+ of the ray of effects simulations using the tests presents for images rived antenna problem. heterogeneous this useful, with complicates otherwise array further the, designs An of beam. portion significant antenna limited a of over measurements This use polarization effect. render this eliminate can not geom- effect but (observing reduce, angle to parallactic tends of etry) range a Observa- over polarizations. tions instrumental com- larger, with pattern or antenna parable, the source of dwarf regions per- can in few response polarization instrumental a this be less, inten- can or polarization total cent antennas fractional to source the beam response As area, antenna sity. polarized limited the spurious this a of from Away introduce area sky. limited the a on only over mance of pointing accurate requires antennas. models the beam for rection hspprdsrbsatcnqefrcretn de- correcting for technique a describes paper This perfor- polarization optimize can design Antenna mtebsciaigtcnqe The technique. imaging basic the om 2 eoee ntssuigunpolarized using tests in recovered nrdoitrfrmti imaging. interoferometric radio in dVadi hw ihsimulated with shown is and V nd eedn riat u oasym- to due artifacts dependent o tksIcreto sabout is correction I Stokes for e aoy ot Africa South vatory, nswt smercpatterns. asymmetric with nnas crneraitcsmltdsky simulated realistic range ic 21) tterns MNRAS 000,1–19 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

Ionizing radiation from AGNs at z > 3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

Ikuru Iwata,1,2,3,4? Marcin Sawicki,4† Akio K. Inoue,5,6 Masayuki Akiyama,7 Genoveva Micheva,8 Toshihiro Kawaguchi,9 Nobunari Kashikawa,10 Stephen Gwyn,11 Stephane Arnouts,12 Jean Coupon,13 Guillaume Desprez13 1Subaru Telescope, National Astronomical Observatory of Japan, 650 North A’ohoku Place, Hilo, Hawaii 96720, USA 2TMT Project, National Astronomical Observatory of Japan, 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 3Department of Astronomical Science, The Graduate University for Advanced Studies (Sokendai), 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 4Department of Astronomy and Physics and Institute for Computational Astrophysics, Saint Mary’s University, 923 Robie Street, Halifax, Nova Scotia B3H 3C3, Canada 5 Department of Physics, School of Advanced Science and Engineering, Faculty of Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku, Tokyo 169-8555, Japan 6 Waseda Research Institute for Science and Engineering, Faculty of Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku, Tokyo 169-8555, Japan 7 Astronomical Institute, Tohoku University, Aramaki, Aoba-ku, Sendai, Miyagi 980-8578, Japan 8 Leibniz-Institute for Astrophysics Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 9 Department of Economics, Management and Information Science, Onomichi City University, Hisayamada 1600-2, Onomichi, Hiroshima 722-8506, Japan 10 Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 11 NRC-Herzberg, 5071 West Saanich Road, Victoria, British Columbia V9E 2E7, Canada 12 Aix Marseille Universit´e, CNRS, Laboratoire d’Astrophysique de Marseille, UMR 7326, F-13388 Marseille, France 13 Astronomy Department, University of Geneva, Chemin d’Ecogia 16, CH-1290 Versoix, Switzerland

Accepted 2021 September 20. Received 2021 September 14; in original form 2021 January 26

ABSTRACT We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum; LyC) escape fraction from AGNs at z ∼ 3 − 4. For 94 AGNs with spectroscopic redshifts at 3.3 < z < 4.0, we use their U- band / i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at z > 3.3 is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find out the average flux ratio of LyC and non-ionizing UV photons escaping from the objects (fLyC/fUV) = 0.182 ± 0.043 for AGNs at 3.3 < z < 3.6, which corresponds to LyC escape fraction fesc = 0.303 ± 0.072 if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function of AGNs, we argue that UV-selected AGNs’ contribution to the LyC emissivity at the epoch is minor, although the size of their

arXiv:2109.10138v1 [astro-ph.GA] 21 Sep 2021 contribution largely depends on the shape of the UV luminosity function. Key words: galaxies: active – galaxies: evolution – galaxies:high-redshift – intergalactic medium – cosmology: obser- vations.

1 INTRODUCTION the Universe. Understanding the relative contributions from these two populations to the ionizing photon budget over cos- Active Galactic Nuclei (AGNs) and massive stars in star- mic time is deeply connected to our understanding of AGN forming galaxies are the two primary sources of hydrogen- and star-formation activity in galaxies at different epochs. ionizing radiation (Lyman Continuum; LyC hereafter) in In light of measurements of the faint-end galaxy UV lu- minosity function (UVLF) at high redshift, it has been sug- ? E-mail: [email protected] gested that faint galaxies are the population primarily re- † Canada Research Chair sponsible for reionizing the Universe (e.g., Inoue et al. 2006;

© 2021 The Authors Astronomy & Astrophysics manuscript no. DiskQuadrants5 ©ESO 2021 September 22, 2021

Quadrant polarization parameters for the scattered light of circumstellar disks. Analysis of debris disk models and observations of HR 4796A

H.M. Schmid

ETH Zurich, Institute for Particle Physics and Astrophysics, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland, e-mail: [email protected]

Received ...; accepted ...

ABSTRACT

Context. Modern imaging polarimetry provides spatially resolved observations for many circumstellar disks and quantitative results for the measured polarization which can be used for comparisons with model calculations and for systematic studies of disk samples. Aims. This paper introduces the quadrant polarization parameters Q000, Q090, Q180, Q270 for Stokes Q and U045, U135, U225, U315 for Stokes U for circumstellar disks and describes their use for the polarimetric characterization of the dust in debris disks. Methods. We define the quadrant polarization parameters Qxxx and Uxxx and illustrate their properties with measurements of the debris disk around HR 4796A from Milli et al. (2019). We calculate quadrant parameters for simple models of rotationally symmetric and optically thin debris disks and the results provide diagnostic diagrams for the determination of the scattering asymmetry of the dust. This method is tested with data for HR 4796A and compared with detailed scattering phase curve extractions in the literature. Results. The parameters Qxxx and Uxxx are ideal for a well-defined and simple characterization of the azimuthal dependence of the polarized light from circumstellar disk because they are based on the “natural” Stokes Q and U quadrant pattern produced by cir- cumstellar scattering. For optically thin and rotationally symmetric debris disks the quadrant parameters normalized to the integrated azimuthal polarization Qxxx/Qφ and Uxxx/Qφ or quadrant ratios like Q000/Q180 depend only on the disk inclination i and the polarized scattering phase function fφ(θ) of the dust, and they do not depend on the radial distribution of the scattering emissivity. Because the disk inclination i is usually well known for resolved observations, we can derive the shape of fφ(θ) for the phase angle range θ sampled by the polarization quadrants. This finding also applies to models with vertical extensions as observed for debris disks. Diagnostic diagrams are calculated for all normalized quadrant parameters and several quadrant ratios for the determination of the asymmetry parameter g of the polarized Henyey-Greenstein scattering phase function fφ(θ, g). We apply these diagrams to the measurement of HR 4796A, and find that a phase function with only one parameter does not reproduce the data well. We find a better solution with a three-parameter phase function fφ(θ, g1, g2, w), but it is also noted that the well-observed complex disk of HR 4796A cannot be described in full detail with the simple quadrant polarization parameters. Conclusions. The described quadrant polarization parameters are very useful for quantifying the azimuthal dependence of the scat- tering polarization of spatially resolved circumstellar disks illuminated by the central star. They provide a simple test of the deviations of the disk geometry from axisymmetry and can be used to constrain the scattering phase function for optically thin disks without detailed model fitting of disk images. The parameters are easy to derive from observations and model calculations and are therefore well suited to systematic studies of the dust scattering in circumstellar disks. Key words. stars: pre-main sequence – planetary systems: debris disks – planetary systems: protoplanetary disks – star: individual object: HR 4796A – Techniques: polarimetric

1. Introduction With AO systems, the observational point spread function (PSF) depends to a large extent on the atmospheric turbulence Circumstellar disks reflect the light from the central star, and and is highly variable (e.g., Cantalloube et al. 2019). For this the produced scattered intensity and polarization contain a lot of reason, the disks are often only detected in polarized light and it information about the disk geometry and the scattering dust par- is not possible to disentangle the disk intensity signal from the arXiv:2109.10099v1 [astro-ph.SR] 21 Sep 2021 ticles. The scattered light of disks is usually only a contribution strong, variable PSF of the central star (see e.g., Esposito et al. of a few percent or less to the direct light from the central star 2020). Therefore, analyses of the scattered light from the disk and therefore requires observations with sufficiently high reso- are often based on the differential polarization alone and only in lution and contrast to resolve the disk from the star. This was favorable cases can one combine this with measurements of the achieved in recent years for many circumstellar disks with adap- disk intensity. For the data analysis, the observed polarization tive optics (AO) systems at large telescopes using polarimetry, must first be corrected for instrumental polarization effects, and a powerful high-contrast technique, to disentangle the scattered this is relatively difficult for complex AO systems (e.g., Schmid and therefore polarized light of the disk from the direct and typ- 2021, and references therein). For this reason, the first generation ically unpolarized light of the central star (e.g., Apai et al. 2004; of AO systems with polarimetric mode provided useful qualita- Oppenheimer et al. 2008; Quanz et al. 2011; Hashimoto et al. tive polarimetric results but hardly any quantitative results. 2011; Muto et al. 2012).

Article number, page 1 of 21 September 22, 2021 0:45 main

Journal of Astronomical Instrumentation © World Scientific Publishing Company

Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case

J. Cumner1†, E. De Lera Acedo1,2, D.I.L. de Villiers3, D. Anstey1, C. I. Kolitsidas4, B. Gurdon1, N. Fagnoni1, P. Alexander1, G. Bernardi5,6,7, H.T.J. Bevins1, S. Carey1, J. Cavillot8, R. Chiello9, C. Craeye8, W. Croukamp3, J.A. Ely1, A. Fialkov2,10, T. Gessey-Jones1, Q. Gueuning1, W. Handley1,2, R. Hills1, A.T. Josaitis1, G. Kulkarni11, A. Magro12, R. Maiolino1,2, P. D. Meerburg13, S. Mittal11, J.R. Pritchard14, E. Puchwein15, N. Razavi-Ghods1, I.L.V. Roque1, A. Saxena13, K.H. Scheutwinkel1, E. Shen1, P.H. Sims16,17, O. Smirnov6,7, M. Spinelli18,19,20, K. Zarb-Adami9,12 1Cavendish Astrophysics, University of Cambridge, Cambridge, UK, 2Kavli Institute for Cosmology in Cambridge, University of Cambridge, Cambridge, UK 3Department of Electrical and Electronic Engineering, Stellenbosch University, Stellenbosch, South Africa 4Ericsson AB, Torshamnsgatan 21, 164 40 Kista, Sweden 5INAF-Istituto di Radio Astronomia, via Gobetti 101, 40129 Bologna, Italy 6Department of Physics and Electronics, Rhodes University, PO Box 94, Grahamstown, 6140, South Africa 7South African Radio Astronomy Observatory, Black River Park, 2 Fir Street, Observatory, Cape Town, 7925, South Africa 8Antenna Group, Universit´ecatholique de Louvain, Belgium 9Physics Department, University of Oxford, Oxford, UK 10Institute of Astronomy, University of Cambridge, Cambridge, UK 11Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 12Institute of Space Sciences and Astronomy, University of Malta, Malta 13Faculty of Science and Engineering, University of Groningen, Groningen, The Netherlands 14Department of Physics, Imperial College London, London, UK 15Leibniz Institute for Astrophysics, Potsdam, Germany 16McGill Space Institute, McGill University, Montr´eal,Canada 17Department of Physics, McGill University, Montr´eal,Canada 18INAF-Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, I-34143 Trieste, Italy 19 IFPU - Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy 20 Department of Physics and Astronomy, University of the Western Cape, Robert Sobukhwe Road, Bellville, 7535, South Africa

Received (to be inserted by publisher); Revised (to be inserted by publisher); Accepted (to be inserted by publisher);

Following the reported detection of an absorption profile associated with the 21 cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21 cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance arXiv:2109.10098v1 [astro-ph.IM] 21 Sep 2021 reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis.

Keywords: 21 cm cosmology; radio antenna; global experiment.

1. Introduction The formation of the first stars in the universe, and associated events, lead to a small but distinct global cosmological 21-cm signature (Furlanetto et al., 2006). The key quantity in deciding this is the spin tem- perature (quantifying the relative population of hyperfine levels of a neutral hydrogen atom), which in

†[email protected]

1 Astronomy & Astrophysics manuscript no. sample ©ESO 2021 September 22, 2021

The GRAVITY Young Stellar Object survey VII. Gas and dust faint inner rings in the hybrid disk of HD141569

GRAVITY Collaboration?: V. Ganci1, 2, L. Labadie1, L. Klarmann3, A. de Valon4, K. Perraut4, M. Benisty4, 5, W. Brandner3, A. Caratti o Garatti3, 6, 7, 16, C. Dougados4, F. Eupen1, R. Garcia Lopez3, 6, 7, R. Grellmann1, J. Sanchez-Bermudez3, 8, A. Wojtczak1, P. Garcia9, 10, A. Amorim9, 11, M. Bauböck12, J.-P. Berger4, P. Caselli12, Y. Clénet13, V. Coudé du Foresto13, P.T. de Zeeuw12, 14, A. Drescher12, G. Duvert4, A. Eckart1, 2, F. Eisenhauer12, M. Filho9, 10, F. Gao12, E. Gendron13, R. Genzel12, S. Gillessen12, G. Heissel13, T. Henning3, S. Hippler3, M. Horrobin1, Z. Hubert4, A. Jiménez-Rosales12, L. Jocou4, P. Kervella13, S. Lacour13, V. Lapeyrère13, J.-B. Le Bouquin4, P. Léna13, T. Ott12, T. Paumard13, G. Perrin13, O. Pfuhl15, G. Heißel13, G. Rousset13, S. Scheithauer3, J. Shangguan12, T. Shimizu12, J. Stadler12, O. Straub12, C. Straubmeier1, E. Sturm12, E. van Dishoeck12, 14, F. Vincent13, S.D. von Fellenberg12, F. Widmann12, and J. Woillez15 (Affiliations can be found after the references)

Received xx, xxxx; accepted xx, xxxx

ABSTRACT

Context. The formation and evolution of planetary systems impact the evolution of the primordial accretion disk in its dust and gas content. HD 141569 is a peculiar object in this context as it is the only known pre-main sequence star characterized by a hybrid disk. Observations with 8 m class telescopes probed the outer-disk structure showing a complex system of multiple rings and outer spirals. Furthermore, interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. Aims. The goal of this work was to explore with new high-resolution interferometric observations the geometry, properties, and dynamics of the dust and gas in the internal regions of HD 141569. Methods. We observed HD 141569 on milliarcsecond scales with GRAVITY/VLTI in the near-infrared (IR) at low (R ∼20) and high (R ∼4000) spectral resolution. We interpreted the interferometric visibilities and spectral energy distribution with geometrical models and through radiative transfer techniques using the code MCMax to constrain the dust emission. We analyzed the high spectral resolution quantities (visibilities and differential phases) to investigate the properties of the Brackett-γ (Brγ) line emitting region. Results. Thanks to the combination of three different epochs, GRAVITY resolves the inner dusty disk in the K band with squared visibilities down to V2 ∼ 0.8. A differential phase signal is also detected in the region of the Brγ line along most of the six baselines. Data modeling shows that an IR excess of about 6% is spatially resolved and that the origin of this emission is confined in a ring of material located at a radius of ∼1 au −8 from the star with a width . 0.3 au. The MCMax modeling suggests that this emission could originate from a small amount (1.4 × 10 M⊕) of quantum-heated particles, while large silicate grain models cannot reproduce at the same time the observational constraints on the properties of near-IR and mid-IR fluxes. The high spectral resolution differential phases in the Brγ line clearly show an S-shape that can be best reproduced with a gaseous disk in Keplerian rotation, confined within 0.09 au (or 12.9 R?). This is also hinted at by the double-peaked Brγ emission line shape, known from previous observations and confirmed by GRAVITY. The modeling of the continuum and gas emission shows that the inclination and position angle of these two components are consistent with a system showing relatively coplanar rings on all scales. Conclusions. With a new and unique observational dataset on HD 141569, we show that the complex disk of this source is composed of a multitude of rings on all scales. This aspect makes HD 141569 a potentially unique source to investigate planet formation and disk evolution in intermediate-mass pre-main sequence stars. Key words. Proto-planetary disk – interferometry – near infra-red – stars: individual (HD 141569)

1. Introduction up to a few hundred au (e.g., Lodato et al. 2019; Benisty et al. 2017). The inner regions (at ∼au scale) of such disks are also The formation and evolution of protoplanetary disks are di- arXiv:2109.10070v1 [astro-ph.SR] 21 Sep 2021 of prime interest since key processes like gas accretion flows, rectly linked to planet formation. The outer disk features of winds, outflows, and dust sublimation take place. All these pro- young stellar objects (YSOs) have been thoroughly studied in the cesses affect the dynamics and evolution of the first few au re- past through scattered light imaging (e.g., with SPHERE/VLT, gions where terrestrial planets may form and/or migrate over few Beuzit et al. 2019) and with ALMA in the (sub-)millimeter range million years. Constraints on these processes can be derived in- (ALMA Partnership et al. 2015). With both techniques, disk ob- directly through spectroscopic studies, but at typical distances servations have shown rings, gaps, and asymmetric structures of a few hundred parsecs only observations with milliarcsec- ? GRAVITY is developed in a collaboration by the Max Planck Insti- ond (mas) resolution, which are required to probe sub-au scales, tute for Extraterrestrial Physics, LESIA of Paris Observatory and IPAG can discriminate between competing models. Numerous inter- of Université Grenoble Alpes / CNRS, the Max Planck Institute for As- ferometric studies have been conducted in the past in the near- tronomy, the University of Cologne, the Centro de Astrofísica e Gravi- and the mid-infrared (IR), for instance with IOTA (Millan-Gabet tação and the European Southern Observatory.

Article number, page 1 of 25 SF2A 2021 A. Siebert, K. Bailli, E. Lagadec, N. Lagarde, J. Malzac , J.-B. Marquette, M. NDiaye, J. Richard, O. Venot (eds)

EXPLORING THE COSMIC DAWN WITH NENUFAR

F.G. Mertens1,2, B. Semelin1 and L.V.E. Koopmans2

Abstract. The exploration of the Cosmic Dawn, the period of the Universe during which the first stars and galaxies were formed, is one of the last frontiers of modern astronomy and cosmology. The redshifted 21-cm line emission from neutral hydrogen is a unique probe and can open this era for astrophysical and cosmological studies. The tentative detection of the 21-cm global signal by the EDGES team at z ∼ 17 underlines the need for an interferometric detection of this signal to discriminate between the numerous models trying to explain this unexpected discovery. The NenuFAR Cosmic-Dawn Key-Science Program aims to perform this detection in the redshift range z ∼ 15 − 31 with a novel SKA precursor, the NenuFAR radio telescope located at the Station de Radioastronomie de Nanay and that started operating in 2019. Due to its compactness it is particularly sensitive to the large scale of the 21-cm signal. Only 100 hours of observation are needed to reach the level of the most extreme models, while 1000 hours are needed for the more standard models. Observations have already started, accumulating to almost 500 hours on the North Celestial Pole field. In this contribution, we introduce the project, our first results and the developments in calibration and RFI mitigation specific to this new instrument.

Keywords: cosmology: dark ages, reionization, first stars; cosmology: observations; techniques: interfero- metric; methods: data analysis

1 Introduction

The observation of the 21-cm line of neutral hydrogen is the most promising method to study the Cosmic Dawn and Epoch of Reionization (EoR) in detail. Its detection can inform us about the timing and mechanisms of the formation of the first stars, as well as the impact on the physics of the interstellar medium and the intergalactic medium (IGM) of the radiation emitted by these first sources of light. Many observational programs are underway to detect it. The so-called “global” experiments, such as EDGES ∗ (Bowman et al. 2018), aim at detecting the 21-cm signal average in the sky. In addition to the global signal, spatial fluctuations in the brightness temperature of the 21-cm signal can be detected by radio-interferometric imaging, using interferometers such as NenuFAR † (Zarka et al. 2020), LOFAR ‡ (van Haarlem et al. 2013), MWA § (Tingay et al. 2013), HERA ¶ (DeBoer et al. 2017), and in a near future SKA ‖ (Koopmans et al. 2015). The power spectrum of the image cube thus obtained allows in particular to characterize the fluctuations of the spin temperature, of the hydrogen density, and of the ionization fraction, on scales ranging from a few arc minutes to a few degrees. The tentative detection reported by the EDGES team of an absorption spectrum at redshift z ∼ 14 − 21 (Bowman et al. 2018), in the range where the Cosmic Dawn is predicted, has unexpected features. This observed absorption signal (-500 mK ± 200 mK, 99 % IC) is considerably deeper than expected (Fraser et al. arXiv:2109.10055v1 [astro-ph.CO] 21 Sep 2021 2018). Standard models predict a maximum absorption of the order of -200 mK and, to explain this observation, a new process must be introduced, either by supercooling the Inter Galactic Medium (IGM) gas (Barkana 2018,

1 LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universit´e,F-75014 Paris, France 2 Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands ∗Experiment to Detect the Global Epoch of Reionization Signature, https://loco.lab.asu.edu/edges †New Extension in Nanay Upgrading LOFAR, https://nenufar.obs-nancay.fr ‡Low-Frequency Array, www.lofar.org §Murchison Widefield Array, www.mwatelescope.org ¶Hydrogen Epoch of Reionization Array, https://reionization.org ‖Square Kilometre Array, www.skatelescope.org

© Soci´et´eFrancaise d’Astronomie et d’Astrophysique (SF2A) 2021 MNRAS 000,1–10 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

Structure, stability and optical absorption spectra of small Ti푛C푥 clusters: a first-principles approach

Sergio Gámez-Valenzuela,1 Julio A. Alonso,2 Gonzalo Santoro3 and José I. Martínez3★ 1Department of Physical Chemistry, University of Malaga, Campus de Teatinos s/n, ES-29071, Malaga, Spain 2Department of Theoretical and Atomic Physics and Optics, University of Valladolid, ES-47011 Valladolid, Spain 3ESISNA Group, Institute of Materials Science of Madrid (ICMM-CSIC), ES-28049 Madrid, Spain

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Titanium-carbide molecular clusters are thought to form in the circumstellar envelopes (CSEs) of carbon-rich Asymptotic Giant Branch stars (AGBs) but, to date, their detection has remained elusive. To facilitate the astrophysical identification of those clusters in AGBs and post-AGBs environments, the molecular structures and optical absorption spectra of small Ti푛C푥 clusters, with n = 1–4 and x = 1–4, and some selected larger clusters, Ti3C8, Ti4C8, Ti6C13, Ti7C13, Ti8C12, Ti9C15, and Ti13C22, have been calculated. The density functional formalism, within the B3LYP approximation for electronic exchange and correlation, was used to find the lowest energy structures. Except the clusters having a single Ti atom, the rest exhibit three-dimensional structures. Those are formed by a Ti fragment surrounded in general by carbon dimers. The optical spectra of Ti푛C푥, computed by time-dependent density functional theory, using the corrected CAM-B3LYP functional, show absorption features in the visible and near infrared regions which may help in the identification of these clusters in space. In addition, most of the clusters have sizable electric dipole moments, allowing their detection by radioastronomical observations. Key words: astrochemistry – methods: numerical – molecular data – stars: AGB and post-AGB – stars: carbon – (stars:) circumstellar matter

1 INTRODUCTION On the other hand, TiC nanocrystals were suggested as carriers of the 21 휇m emission feature in C-rich protoplanetary nebulae (PPNe) Titanium carbide (TiC) dust is predicted to form in the inner re- (Helden et al. 2000). However, the low abundance of Ti in these gions of the circumstellar envelopes (CSEs) of carbon rich (C-rich) environments is not consistent with this carrier (Chigai et al. 2003; Asymptotic Giant Branch stars (AGBs) (Agúndez, M. et al. 2020). Zhang et al. 2009) and the possibility of a hydrocarbon-based carrier Due to its high refractory nature, TiC is considered to be the first con- has been considered instead, although coated TiC grains could cir- densate in the atmosphere of C-rich AGBs and may act as seed for cumvent the constrains imposed by the Ti abundance in post-AGB the heterogeneous nucleation of dust. The identification of TiC crys- environments (Sloan et al. 2014). tals embedded in graphite spherules extracted from the Murchinson meteorite supports this prediction (Bernatowicz et al. 1991; Berna- The solid-state phase diagram of the Ti-C alloy (Okamoto 1998) towicz et al. 1996) but, to date, the detection of molecules comprised shows a very stable compound centered at Ti55C45 with a wide of titanium and carbon that may serve as precursors for the formation composition range. This compound has a high melting temperature of TiC dust has remained elusive. of 3066◦C. These features reveal that the formation of strong Ti-C In the inner regions of the CSE of AGBs, titanium is in atomic form bonds is favorable, and this is also the case for Ti–C nanoclusters, arXiv:2109.10050v1 [astro-ph.GA] 21 Sep 2021 presenting fairly low abundances. The reaction of Ti with acetylene which exhibit highly stable structures at particular stoichiometries. (C2H2) has been suggested as a possible pathway for the formation of stoichiometric TiC molecules (Chigai et al. 1999; Gail & Sedlmayr Three decades ago, by reacting metal vapors with hydrocarbons, 2013), and chemical equilibrium calculations of AGBs atmospheres Castleman and coworkers discovered a class of highly stable cage- show that the depletion of atomic Ti starts at stellar radii in which like clusters with formula M8C12, the metallocarbohedrenes, or met- cars (Guo et al. 1992a; Wei et al. 1992; Guo et al. 1992b; Rohmer et al. most carbon is locked into C2H2. However, the predicted abundance of the diatomic molecule TiC is too low and the high stability of large 2000), where M denotes an early transition metal (mainly, Ti, V, Zr, Nb, Hf, Mo, Cr, and Fe). Later on, Yu et al.(1995) performed similar clusters such as Ti8C12, which might be formed through the reaction of smaller Ti푛C푥 clusters, suggests that it is a reasonable candidate experiments by reacting laser-vaporized Ti with CH4, and observed as molecular precursor of TiC dust. Indeed, chemical equilibrium the formation of Ti8C12 along with several other clusters including calculations show that Ti C replaces atomic Ti as the main Ti Ti푛C2푛 stoichiometries. The mass spectrometric studies of anionic 8 12 − reservoir at 2–3 stellar radii (Agúndez, M. et al. 2020). Ti푛C푥 clusters performed by Wang et al.(1997, 1998) revealed new stable anionic clusters and evidenced the importance of C2 dimers in stabilizing the structure of these clusters. These authors measured ★ E-mail: [email protected] the photoelectron spectra of the most stable anionic clusters, namely:

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

Unequal-mass mergers of dark matter haloes with rare and frequent self-interactions

Moritz S. Fischer,1★ Marcus Brüggen,1 Kai Schmidt-Hoberg,2 Klaus Dolag,3,4 Antonio Ragagnin,5,6 Andrew Robertson7 1Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 2Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, D-22607 Hamburg, Germany 3Universitäts-Sternwarte München, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, D-81679 München, Germany 4Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, D-85748 Garching, Germany 5INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, I-34143 Trieste, Italy 6IFPU - Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy 7Institute for Computational Cosmology, Department of Physics, Durham University, South Road, Durham DH1 3LE, UK

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Dark matter self-interactions have been proposed to solve problems on small length scales within the standard cold dark matter cosmology. Here we investigate the effects of dark matter self-interactions in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform 푁-body dark matter-only simulations of idealised setups to study the effects of dark matter self-interactions that are elastic and velocity-independent. We go beyond the commonly adopted assumption of large-angle (rare) dark matter scatterings, paying attention to the impact of small-angle (frequent) scatterings on astrophysical observables and related quantities. Specifically, we focus on dark matter-galaxy offsets, galaxy-galaxy distances, halo shapes, morphology and the phase-space distribution. Moreover, we compare two methods to identify peaks: one based on the gravitational potential and one based on isodensity contours. We find that the results are sensitive to the peak finding method, which poses a challenge for the analysis of merging systems in simulations and observations, especially for minor mergers. Large dark matter-galaxy offsets can occur in minor mergers, especially with frequent self-interactions. The subhalo tends to dissolve quickly for these cases. While clusters in late merger phases lead to potentially large differences between rare and frequent scatterings, we believe that these differences are non-trivial to extract from observations. We therefore study the galaxy/star populations which remain distinct even after the dark matter haloes have coalesced. We find that these collisionless tracers behave differently for rare and frequent scatterings, potentially giving a handle to learn about the micro-physics of dark matter. Key words: astroparticle physics – methods: numerical – galaxies: haloes – dark matter

1 INTRODUCTION well by 푁-body simulations assuming only CDM, there seem to be discrepancies between predictions and observations on scales of In the standard cosmological model, structures in the Universe such as galaxies. This is often referred to as the small-scale crisis of ΛCDM, galaxies and galaxy clusters are thought to form hierarchically, mean- which manifests itself in the core-cusp problem, the diversity prob- ing that small objects merge to form larger ones (e.g. Lacey & Cole lem, the too-big-to-fail-problem and the plane of satellites problem 1993). In particular, cosmological 푁-body simulations have been arXiv:2109.10035v1 [astro-ph.CO] 21 Sep 2021 (for a review see Bullock & Boylan-Kolchin 2017). used to investigate structure formation within the standard ΛCDM model (e.g. Springel et al. 2005; Boylan-Kolchin et al. 2009; Klypin Many solutions have been proposed to solve these problems, in- et al. 2011; Pillepich et al. 2017; Hopkins et al. 2018). These sim- cluding alternative models for DM (e.g. Dodelson & Widrow 1994; ulations agree remarkably well with observations of the large-scale Hu et al. 2000). One class of these models assumes that DM is not structure (e.g. Springel et al. 2006) and allow us to explain the forma- fully collisionless but that DM particles scatter off each other with tion of a wide range of cosmic objects. In the cosmological standard a non-zero cross-section (Spergel & Steinhardt 2000). Such mod- model, today’s Universe mainly consists of dark energy and dark els with self-interacting dark matter (SIDM) provide a promising matter (DM), which is thought to be collisionless and cold (cold dark solution to the small-scale crisis (for a review see Tulin & Yu 2018). matter - CDM). The underlying nature of those two components is A variety of SIDM models exist for which the self-interactions largely unknown. can be velocity-independent or velocity-dependent as several studies Although the large-scale structure of our Universe is reproduced have assumed (e.g. Colin et al. 2002; Vogelsberger et al. 2012; Vo- gelsberger & Zavala 2013; Vogelsberger et al. 2014; Robertson et al. 2017b; Banerjee et al. 2020; Nadler et al. 2020; Turner et al. 2021). ★ E-mail: moritz.fi[email protected] (UHH) The self-interactions could be elastic or inelastic (e.g. Essig et al.

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

Detection of Small Scale Components in Power Law Spectra

Tim Ruhe푎,∗ and Wolfgang Rhode푎 푎TU Dortmund University, Otto Hahn-Straße 4a, 44227 Dortmund, Germany E-mail: [email protected], [email protected]

Spectra in astroparticle physics are commonly approximated by simple power laws. The steeply falling nature of these power laws, however, makes the detection of additional components rather challenging. This holds true especially, if the additional components are small compared to the established ones. Energy spectra of muon neutrinos are an interesting example of such a scenario, where the conventional and astrophysical components to the spectra have been established by the use of different analysis methods, such as likelihood fits or spectral deconvolution. The prompt component, although expected from theoretical models, has not yet been experimentally observed. Furthermore, the extraction of physics parameters is challenged by the large systematic uncertainties, especially at high energies. This contribution presents a different approach to the analysis of power-law spectra, which is based on functional data analysis. The method itself and its implications are discussed using muon and neutrino energy spectra as an example. arXiv:2109.10032v1 [astro-ph.HE] 21 Sep 2021

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

∗Presenter

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

Detection of Small Scale Components Tim Ruhe

10 2 1 r s

1 s

2

m 10 3 c

V e G

/ 7 . 2 E Hoerandel 2003 10 4 Gaisser 2012 IC79 Unfolding Antares

102 103 104 105 106 107 Energy/GeV

Figure 1: Muon neutrino energy spectra measured by Antares (red) and IceCube (green), compared to theoretical fluxes obtained with MCEq for the cosmic ray models by Gaisser and Hoerandel using SYBILL- 2.3c as a hadronic interaction model.

1. Introduction

The measurement of lepton energy spectra via the use of deconvolution techniques has become a routine task for state-of-the-art neutrino telescopes (cf. [1–3]). An observed flattening of the spectrum for energies exceeding 105 GeV has been attributed to a diffuse flux of high energy neutrinos from astrophysical sources [2]. Flux measurements therefore not only complement other analysis techniques, but also have the capability to confirm results in more model-independent analyses. Atmospheric neutrino spectra are expected to consist of three components, the first one being conventional atmospheric neutrinos, which originate from the decay of pions and kaons [4–6]. Due to their relatively long lifetime (휏 ≈ 10−8 s [7]), the parent mesons can lose energy in collision before decaying into neutrinos. Hence, the flux of conventional neutrinos is one power steeper 푑Φ −3.7 (approximately 푑퐸 ∝ 퐸 ), compared to the spectrum of the incident cosmic rays. So-called prompt neutrinos, which are expected from the decay of charmed mesons, but have not yet been experimentally detected, are the second component [4,5,8]. A diffuse flux of neutrinos from astrophysical sources, contributes to the spectrum as a third component(cf. [9, 10]). Figure1 shows muon neutrino energy spectra, obtained with Antares [3] and IceCube-79 [2]. In addition, theoretical predictions obtained with MCEq [11] for the cosmic ray models by Hoeran- del [12] and Gaisser [13] are depicted. For both models, the zenith-averaged muon neutrino fluxes were computed, using SYBILL-2.3c [14] as a hadronic interaction model. The depicted fluxes also include a diffuse flux of astrophysical neutrinos, which is modelled as a power law of the form:

  훾astro 푑Φ 퐸휈 = Φ ,astro × , (1) 푑퐸 0 100 TeV

2 J. Astrophys. Astr. (0000) 000: #### DOI

Wave Effects in Double-Plane Lensing

RAHUL RAMESH1, ASHISH KUMAR MEENA1,2 and JASJEET SINGH BAGLA1*

1Department of Physical Sciences, IISER Mohali, Sector 81, SAS Nagar, Punjab, India - 140306. 2Physics Department, Ben-Gurion University of the Negev, P.O. Box 653, Be’er-Sheva 8410501, Israel

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

MS received xx yy 2021; accepted xx yy 2021

Abstract. We discuss the wave optical effects in gravitational lens systems with two point mass lenses in two different lens planes. We identify and vary parameters (i.e., lens masses, related distances, and their alignments) related to the lens system to investigate their effects on the amplification factor. We find that due to a large number of parameters, it is not possible to make generalized statements regarding the amplification factor. We conclude by noting that the best approach to study two-plane and multi-plane lensing is to study various possible lens systems case by case in order to explore the possibilities in the parameter space instead of hoping to generalize the results of a few test cases. We present a preliminary analysis of the parameter space for a two-lens system here.

Keywords. gravitational waves—microlensing—multi-plane lensing.

1. Introduction objects, as the wavelength of GW signal in this case is comparable to the Schwarzschild radius of the lens In the last few decades, gravitational lensing has be- (RSch ∼ λGW). As a result, wave effects becomes sig- come a useful tool to probe properties of the Uni- nificant, and this leads to frequency dependent ampli- verse (Bartelmann, 2010; Kneib & Natarajan, 2011; fication of GWs. These wave effects due to an iso- Umetsu, 2020). However, most of its applications are lated single/double point mass lens (e.g., Takahashi & discussed in the context of gravitational lensing of Elec- Nakamura, 2003; Christian et al., 2018) or a point mass tromagnetic (EM) radiation originating from a distant with external effects (e.g., Diego et al., 2019; Meena & source. Recent detection of gravitational wave (GW) Bagla, 2020; Mishra et al., 2021) have been discussed signals from merging compact objects by Laser Inter- in recent studies and it has been shown that if the point ferometer Gravitational-Wave Observatory (LIGO; Ab- mass lens is embedded in an external environment then bott et al., 2019, 2020) has opened a new avenue for the wave effects may obtain a significant boost. application of gravitational lensing. In the above studies, only a single lens plane was GWs are subject to deflection due to intervening considered. Another interesting scenario is the lens- matter as they travel from the source to the observer ing due to a double plane lens. Kochanek & Apos- (Ohanian, 1974; Baraldo et al., 1999) in a manner sim- tolakis(1988) first explored the framework of double ilar to EM waves. As the GW signal in these direct plane lensing due to galaxy scale lenses and pointed observations originates from sources at cosmological out that the corresponding population is expected to be distances, the possibility of strong lensing due to an approximately 1% − 10% of all lensed systems. Al-

arXiv:2109.09998v1 [astro-ph.CO] 21 Sep 2021 intervening galaxy or galaxy cluster is significant (Li though the gravitational lensing of gravitational waves et al., 2018; Smith et al., 2018; Broadhurst et al., 2019, due to galaxy scale lenses in double plane configura- 2018, 2020). For gravitational lensing of GW signals tion can be studied using the conventional geometric in the LIGO frequency range, lensing due to a galaxy approach, the microlenses residing in these lens galax- or a galaxy cluster can be described in the same man- ies can lead to wave effects. As these microlenses are ner as that of EM radiation, i.e., in the geometric op- lying in two different lens planes, the corresponding tics limit. This is possible as the Schwarzschild ra- wave effects can be significantly different from the case dius of the lens is much larger in comparison to the of one or two-point mass lenses in single plane. Erdl wavelength of the GW signal (RSch  λGW). How- & Schneider(1993) presented a detailed study of two- ever, the same is not true for gravitational lensing of point mass lens in double plane lensing in the geometric GWs (in LIGO frequency range) due to stellar mass optics limit and pointed out various differences between

© Indian Academy of Sciences 1 MNRAS 000,1–17 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

The two phases of core formation - orbital evolution in the centres of ellipticals with supermassive black hole binaries

M. Frigo1,2?, T. Naab1, A. Rantala1, P. H. Johansson3, B. Neureiter4, J. Thomas4, and F. Rizzuto1 1Max-Planck-Institut f¨urAstrophysik, Karl-Schwarzschild-Straße 1, 85748 Garching, Germany 2Excellence Cluster ORIGINS, Boltzmannstraße 2, 85748 Garching, Germany 3Department of Physics, University of Helsinki, Gustaf H¨allstr¨omin katu 2, 00014 Helsinki, Finland 4Max-Planck-Institut f¨urextraterrestrische Physik, Gießenbachstraße 1, 85748 Garching, Germany

Accepted 2021 September 18. Received 2021 September 17; in original form 2021 May 25

ABSTRACT The flat stellar density cores of massive elliptical galaxies form rapidly due to sinking supermassive black holes (SMBH) in gas-poor galaxy mergers. After the SMBHs form a bound binary, gravitational slingshot interactions with nearby stars drive the core regions towards a tangentially biased stellar velocity distribution. We use collisionless galaxy merger simulations with accurate collisional orbit integration around the central SMBHs to demonstrate that the removal of stars from the centre by slingshot kicks accounts for the entire change in velocity anisotropy. The rate of strong (unbinding) −1 kicks is constant over several hundred Myr at ∼ 3 M yr for our most massive SMBH 10 binary (MBH = 1.7 × 10 M ). Using a frequency-based orbit classification scheme (box, x-tube, z-tube, rosette) we demonstrate that slingshot kicks mostly affect box orbits with small pericentre distances, leading to a velocity anisotropy of β . −0.6 within several hundred Myr as observed in massive ellipticals with large cores. We show how different SMBH masses affect the orbital structure of the merger remnants and present a kinematic tomography connecting orbit families to integral field kinematic features. Our direct orbit classification agrees remarkably well with a modern triaxial Schwarzschild analysis applied to simulated mock kinematic maps. Key words: galaxies: kinematics and dynamics – galaxies: formation – galaxies: super massive black holes – methods: numerical

1 INTRODUCTION show little net rotation at high velocity dispersions (Em- sellem et al. 2007; Cappellari 2016; Ene et al. 2020). In gen- Many properties of massive early-type galaxies correlate eral, photometric, stellar population, and kinematic prop- with the mass of their central super-massive black holes erties of massive ETGs are consistent with a ‘two-phase’ (SMBH), most notably galactic bulge mass and stellar ve- formation process, consisting of a rapid formation at high locity dispersion (Dressler 1989; Gebhardt et al. 2000). This redshift by in-situ star formation, followed by mass assem- suggests a connection between SMBHs and their massive bly through mergers at redshifts z 2 (e.g. Naab et al. arXiv:2109.09996v1 [astro-ph.GA] 21 Sep 2021 host galaxies, however, one that is not yet completely un- . 2009; Oser et al. 2010; Johansson et al. 2012; Naab et al. derstood (Kormendy & Ho 2013; Saglia et al. 2016). The 2014; Penoyre et al. 2017; Rodriguez-Gomez et al. 2016; Xu most massive and luminous (M < −20.5) early-type galax- B et al. 2019; Moster et al. 2020). Here, feedback from accret- ies (ETGs) in the Universe are particularly interesting. They ing SMBHs is thought to be of fundamental importance for have stellar masses of order 1012 M or larger, typically show setting the old ages of these galaxies by suppressing or ter- flat stellar density profiles in their cores and host the most minating star formation (see e.g. Somerville & Dav´e 2015; massive SMBHs that have been observed. The measured 9 Naab & Ostriker 2017, for reviews). SMBH masses are of order 10 M or in some cases even 1010 M (e.g. McConnell et al. 2011; Kormendy & Ho 2013; SMBHs are also thought to be responsible for the ob- Thomas et al. 2016; Mehrgan et al. 2019). served flat stellar density distributions in the centres of mas- Stellar kinematics in most of these very old galaxies sive ellipticals (Lauer 1985). The ejection of stars during the sinking and coalescence of binary SMBHs in mergers of ? E-mail: [email protected] (MF) massive ETGs is a very likely formation mechanism of such

© 2021 The Authors arXiv:2109.09987v1 [astro-ph.SR] 21 Sep 2021 Abstract email: Astrophysics, India. of 034, Institute Bengaluru-560 Indian at working Formerly India. 111, Bengaluru-560 Bikasipura, © Javaraiah J. of Groups Velocity Sunspot Meridional Between and Correlation Residual in Rotation Variations of Study A 10.1007/ DOI: Physics Solar eoiyo uso rusi ieet5 different in t groups and sunspot rate Photoheliograph rotation of residual Debrecen velocity mean yearly and the determined (GPR) and (DPD) Results Photoheliograpic wich [email protected] aee-oeec,adwvltpaedffrneaaysso t of analyrses cr done difference slope also wavelet-phase have We the and 1920–1940. of during wavelet-coherence, weak spectrum very wavelet was it Morlet but window, The the years. of that 9 existence suggests and the (SN) 4 number about 3-year by sunspot in SN data mean the yearly equat from and toward determined intervals slopes transport ind the momentum slope conventio between angular the correlation sign strong of the value a negative per of significant existence As the deviation). motions, standard meridional the its than larger times n h eiinlmto.Tersda oainrt sfudt be to found is rate rotation residual residual the The both motion. in meridional dependence the latitude–time and considerable a exists m ini oiatmsl rudtemxm fslrcce ihvlct 8 velocity with cycles solar of maxima the s equat around The m mostly motions. dominant meridional exist is poleward There tion minima. and during equatorward of of that bands than maxima during weaker extent u iharltvl ekvlct,ol – s m 4–6 only velocity, weak relatively a with but etfi sfudt eraoal od(slope, good reasonably be to found linear- corresponding is The rotation. best-fit residual b and exists motion anticorrelation significant meridional within migra statistically poleward rotation A or fr average equatorward motions. clear meridional than years no equatorw suggests faster the of analysis and This to existence slower belt. the according of suggests folded bands maxima alternate are of that epochs 12–24 respective Cycles Solar during Springer s − − 1 1 [email protected];[email protected] hra h oeadmto sdmnn otyaon h minima the around mostly dominant is motion poleward the whereas , o8 m 80 to •••• eaaye h obnd12yassnptgopdt rmGre from data sunspot-group years 142 combined the analyzed We ••••• ≈ - s ••• − 1ya eidct ntesoeams hogottedata the throughout almost slope the in periodicity 11-year 1 - nalrenme fslrcce h oaini some is rotation the cycles solar of number large a In . ••• - •••• - • ◦ ; aiueitras efidta there that find We intervals. latitude − − 0 1 . h nlsso h data the of analysis The . 028 ± 0 . 0,i bu 3.5 about is 008, esoeadSN. and slope he r h cross- The or. r migrating ard emeridional he igbnsof bands ting least-squares oigtime moving oss-wavelet, h activity the radmo- orward lp leads slope tenthe etween sdfor used n ctsthe icates alternate mtheir om suggests rotation cData ic − 120 –12 en- Inverse mass cascade of self-gravitating collisionless flow and effects on halo mass functions

Zhijie (Jay) Xu1,a

1. Computational Mathematics Group, Physical and Computational Sciences Directorate, Pacific

Northwest National Laboratory, Richland, WA 99352, USA

Abstract:

This paper proposes the inverse mass cascade for the statistically steady state of self-gravitating

collisionless flow (SG-CFD). The continuous mass transfer from small to large mass scales

(inverse cascade) is formulated and direct effect on halo mass functions is discussed with new

mass function proposed. Mass cascade is local, two-way, and asymmetric in mass space. Halos

inherit/pass their mass mostly from/to halos of similar size. The net mass transfer proceeds in a

“bottom-up” fashion. Two distinct ranges can be identified, i.e. a propagation range with a scale-

independent rate of mass transfer and a deposition range with cascaded mass consumed to grow

halos. Dimensional analysis leads to a simple power-law mass function in propagation range with

(λ-1) -1 an exponent (λ) dependent on halo geometry. A fundamental merging frequency f0~mp a (same as the redshift of photon frequency) is identified, where a is the scale factor. The particle mass mp

-1 can be determined if frequency is known. The rate of mass transfer εm~a is independent of halo

mass (a key feature of propagation range) and decreasing with a. Typical halos of mass grow as

3/2 -λ mh~a and the waiting time (halo lifespan) scales as mh . Chain reaction can be used to describe

the mass cascade as a sequence of mass accretion and provides non-equilibrium systems a

mechanism to continuously release energy and maximize entropy when evolving toward the

limiting equilibrium. A continuous injection of mass (“free radicals”) at the smallest scale is

required to sustain the everlasting inverse mass cascade such that the total halo mass Mh increases

a) Electronic mail: [email protected]; [email protected]

1

arXiv:2109.09981v1 [astro-ph.CO] 21 Sep 2021 nwta h pndrcino h akmte halos matter dark the of direction spin the matter that dark know the the unclear and between is galaxy relation halo the the of Namely, momentum the angular for halos. predictions matter have dark models matter, baryonic theoretical the most observe while we that is challenges main theoretically the alignments predict to momentum underway angular been lensing have galaxies’ the efforts contaminate many could data, correlations such Since in known cor- simulations is the the universe data the and hand, both of for momentum (LSS) other parameters structure angular the large-scale of galaxy On between handful relation a cosmology. down precision nail align- can galaxy surveys to sensitive ment are surveys galaxy-lensing oetmo aaisi seta o ohgalaxy- both for data galaxy-lensing essential and is theories evolution galaxies of momentum oi tal. et Codis ( 2008 nwn h rgnadeouino h angular the of evolution and origin the Knowing yee sn L using 2021 Typeset 22, September version Draft ); Schaefer Keywords: at flawlessly almost 5pret eso httdltru hoyfist rdc corr predict to fails theory tenso torque tidal tidal the that show of We eigendirection percent. middle 15 the with aligned preferentially ihtemdl iedrcino h ia esr oevr efind we Moreover, sp the tensor. that tidal predicts the theory of ex torque the eigendirection tidal that middle that find the at show We with We percent th 10 axis. use median. about rotation we is the halos’ speed prediction, the and this of spin examine between alignment To bet the alignment for halos. partial look high-speed predicts large theory the locking for tidal m T angular the ha While their that matter change (TLT). argues drastically dark Our could Theory the haloes ages. Locking of between non-linear encounters during momentum Tidal evolution angular momentum it angular initial the dubbed explains the we for responsible locking; is tidal to similar ekl&Schaefer & Merkel eivsiaetesi lgmn ftedr atrhlsb conside by halos matter dark the of alignment spin the investigate We hn tal. et Zhang ( hsr tal. et Chisari 2015a A in tal. et Jiang T ( 1. E 2009 X pnainet—TdlTru hoy—Dr atrHalos Matter Dark — Theory Torque Tidal — alignment spin ); INTRODUCTION twocolumn ); agee al. et Laigle ( 2015 oi tal. et Codis ( ( ( 2015 z 2019 2013 ); pnAineto akMte ao:MdHalos Mad Halos: Matter Dark of Alignment Spin 10. = rli tal. et Kraljic tl nAASTeX631 in style ); .Mroe,tdywe today Moreover, ). .I atclr the particular, In ). ha Ebrahiman Ehsan ( oi tal. et Codis 2015 ( 2015b .Oeo the of One ). 1 et tal. et Berta ( .These ). hrfUiest fTechnology of University Sharif 2021 ( 2012 and ) ABSTRACT 1 ern Iran Tehran, ). zd Ave. Azadi and z o athle—ihvlcte agrta twice than larger velocities haloes—with fast for 0 = labrAbolhasani Aliakbar igehl,adb sn h edvc approximation Zeldovich the using Zeldovich by and halo, single a infalling the on of torque tensor a inertial matter. exerts the that matter and infalling tidal matter the the ambient between the misalignment of the mo- tensor is angular halos the for of explanation mentum perhaps primary the and that unnecessary argue is symmetry matter misleading spherical collapsing assuming the that for later out and Dorshkevich pointed stage. White linear the ac- momentum during angular quisition estimate to theory perturbation its changes that matter infalling momentum angular the field on tidal torque gravitational a applies background form momentum. to the collapses galaxy/halo, angular matter a the the when of that suggested origin Hoyle the work. for this explanation in problem latter the on focus mainly predictions will the from deviates e nua oetmaqie yaclasn spheri- region collapsing overdense a cal by acquired momentum angular net ob oepeie aealteprilsta n pin up end that particles the all take precise, more be To well-settled viable a has (TTT) Theory Torque Tidal ( c lgmn at alignment ect 1970 Doroshkevich mnu.Tecretmanuscript current The omentum. ensedadtesi direction spin the and speed ween no aotnst ealigned be to tends halo a of in ,wietepyia oiino them, of position physical the write ), rvoswr hwdta close that showed work previous iha xespoaiiyof probability excess an with r lutiTGsmlto and simulation IllustrisTNG e o,teTdllcigTheory locking Tidal the los, httehlsat halos the that espoaiiyo alignment of probability cess igamcaimsomewhat mechanism a ring 1 Hoyle Peebles ( ( 1970 1951 dlTru Theory Torque idal z ( hl tworks it while 0 = 1969 ´ pze al. L´opez et ); .Pelsetmtdthe estimated Peebles ). White .H sdtelinear the used He ). z 0are 10 = ( 1984 ( 2021 .They ). .We ). r , arXiv:2109.09980v1 [astro-ph.SR] 21 Sep 2021 itneidctr,wihldt h icvr fteaccel ( the Universe of the discovery of the expansion to ing cosm led accurate which as used indicators, been distance inte have Ia in SNe res (WDs) systems. Ia) dwarfs binary (SNe acting white supernovae of Ia explosions Type thermonuclear that from accepted generally is It INTRODUCTION 1. hne l 2018 2021 al. et Shen ai sapoiigmdlt xli omlSeI (e.g., Ia SNe normal explain to model promising a is nario D oacmlt elyro t ufc.I h eshell He the If surface. its on He layer a or He star a He-burning accumulate a be to could WD, which star, companion rich oe o xliigSeI ( 2014 Ia SNe explaining for un- posed an remain Di however, problem. mechanism, explosion solved their and Ia SNe 1998 etme 2 2021 22, September Astronomy bevtoa intrso h uvvn oo tri the in star donor surviving the of signatures Observational h u-hnrska asdul-eoain(Dt sce- (DDet) double-detonation mass sub-Chandrasekhar The .I hsseai,aW crtsmtra rmisHe- its from material accretes WD a scenario, this In ). ). ; e words. Key star. donor surviving the of motion current the iayt rvla eoiyo bu 0 ms km 400 about of velocity a at travel to binary h u-hnrska asdul-eoain(Dt scen (DDet) double-detonation mass sub-Chandrasekhar The 2021 22, September S78 nta,i S78i nedtesriigH-trdon He-star surviving the indeed is 708 US if Instead, 708. US eim(e trcmainwti h Dtseai.W a t map We scenario. DDet the within companion three-dimension star perform (He) consistently helium explosio we the time, survive first to the expected is For scenario DDet the in star ytmwt uhamsieH-trdnrcno e ls eno close get cannot on donor Based He-star 708. massive US a of such properties with observed system the explain to needed tntrlcniaefra jce oo enn faDe SN DDet a of remnant donor ejected an for candidate natural it yevlct trU 0,wihi h ats non star unbound fastest the is which 708, US not star are hypervelocity properties observational its equilibrium, thermal opnosbcm infiatyoelmnu o bu 10 about for overluminous significantly become companions bevtos uhsriigcmain r jce ihh with ejected are D companions surviving of Such companions observations. He-star su surviving the of of evolution signatures long-term observable the follow and codes evolution 9 8 7 6 4 5 3 1 2 elutre l 1999 al. et Perlmutter aoaoy atLnig I48412,USA 48824-1321, MI Superc Lansing, National East Laboratory, Astrophysics, Nuclear for Institute Joint etro xelnefrAtohsc nTreDmnin ( Dimensions Three in Astrophysics for Excellence of Center utainRsac oni eteo xelnefrGravi for Excellence of Centre Council Research Australian colo hsc n srnm,Mns nvriy 9Rain 19 University, Monash Astronomy, and Physics of School etu ü srnmedrUiesttHiebr,Instit Heidelberg, Universität der Astronomie für Zentrum edlegrIsiu ü hoeiceSuin Schloss- Studien, Theoretische für Institut Heidelberger Germany nvriyo hns cdm fSine ejn 002 P 100012, Beijing Science, of Academy Chinese of University unnOsraois hns cdm fSine,36Yan 396 Sciences, e-mail: of Academy Chinese Observatories, Yunnan e aoaoyfrteSrcueadEouino Celestial of Evolution and Structure the for Laboratory Key & [email protected] Astrophysics hn-e Liu Zheng-Wei ; tr:sproa:gnrl–bnre:close binaries: – general supernovae: stars: onlye l 2019 al. et Townsley ff rn rgntrseaishv enpro- been have scenarios progenitor erent eoainmdlo yeI supernovae Ia Type of model detonation .Teseicpoeio ytm of systems progenitor specific The ). aucitn.ms no. manuscript ilbad ta.2013 al. et Hillebrandt is ta.1998 al. et Riess 1 , 2 , 3 rerc .Röpke K. Friedrich , ; rnwe l 2020 al. et Gronow − ; 1 tcud o xml,hv eneetdfo lblrclu globular a from ejected been have example, for could, It . cmd tal. et Schmidt ; aze al. et Maoz estv otedtiso h jcadnritrcin W interaction. ejecta-donor the of details the to sensitive ofbunne 5 91 edleg Germany Heidelberg, 69118 35, Wolfsbrunnenweg l(D yrdnmclsmltoso h neato fSN of interaction the of simulations hydrodynamical (3D) al n ob jce ttehg ria eoiyo opc bi compact a of velocity orbital high the at ejected be to and n g eoiiso pt bu 3 ms km 930 about to up of velocities igh nutn ylto aoaoy ihgnSaeUniversit State Michigan Laboratory, Cyclotron onducting erat- ABSTRACT 6 e N a hc ilspottesac o uhcompanions such for search the support will which Ia, SNe Det 4 vvn esa opnos u anga st rvd h p the provide to is goal main Our companions. He-star rviving roi otmoaymdlfrTp asproa SeIa). (SNe supernovae Ia Type for model contemporary a is ario tfrTertsh srpyi,Pioohne 2 69120 12, Philosophenweg Astrophysik, Theoretische für ut , rdrn h hra eeulbainpae fe h sta the After phase. re-equilibration thermal the during yr norGlx,taeln ihavlct faot120k s km 1,200 about of velocity a with travelling Galaxy, our in ult 5 ainlWv icvr OGa) lyo,VC30,Austr 3800, VIC Clayton, (OzGrav), Discovery Wave tational bet,CS umn 526 ..China P.R. 650216, Kunming CAS, Objects, ic r- ata Zeng Yaotian , SR-D,Srmo C 61 Australia 2611, ACT Stromlo, ASTRO-3D), , g ttemmn fS xlso oepantehg velocity high the explain to explosion SN of moment the at ugh R China .R. a efidta esa oo iha nta asof mass initial an with donor He-star a that find We Ia. ro DtS a twudrqieteetr r-Nprogenit pre-SN entire the require would it Ia, SN DDet a of or eotoe f3 matsmltosit n-iesoa s one-dimensional into simulations impact 3D of outcomes he u ealdbnr vlto acltos oee,tep the however, calculations, evolution binary detailed our oetWl,VC30,Australia 3800, VIC Walk, forest fnwn,Gad itit umn 526 ..China P.R. 650216, Kunming District, Guandu gfangwang, ece rtclmass, critical a reaches oes al Vfls intrso jcaitrcinwith interaction ejecta of signatures ( companion flash UV Early models. 2021 2013 Woosley & Moll 1995 Arnett mate a & undergoes Livne core then the (e.g., WD of explosion mass thermonuclear detonation sub-Chandrasekhar entire second The a rial. ignites that detonation 2016 al. et Neunteufel N a o ntne h ako msini NI spectra Ia SN in emission H of lack observati the ( current instance, explaining For Ia. in SNe features attractive its of aini S mgn ( imaging HST in panion uain hwta h Dtseai a xli ag fra large a explain ca can synthesis scenario population DDet the addition, that In show explosions. culations the SN at close of very generally ment are systems binary because nario enr 2007 Leonard .TeDe cnrohsrcnl andatninbecause attention gained recently has scenario DDet The ). 1 ae 2010 Kasen , 2 , n h o-eeto fapeepoincom- pre-explosion a of non-detection the and ) 3 n lxne Heger Alexander and , ; ; ; ike l 2007 al. et Fink ≃ oi ta.2019 al. et Polin hudb eysali h Dtsce- DDet the in small very be should ) rnwe l 2020 al. et Gronow − 0 1 ie l 2011 al. et Li . efidta u uvvn He-star surviving our that find We . 02–0 am1980 Taam . M 2 ⊙ ril ubr ae1o 11 of 1 page number, Article tri h ieto of direction the in ster ( , pl u eut to results our apply e ose ae 2011 Kasen & Woosley r neett DDet to inherent are ) 2010 ,i rgesa initial an triggers it ), 6 ; , 7 ose ta.1986 al. et Woosley , , re-establishes r 8 ,1Cyclotron 1 y, ; , 2021 jcawt a with ejecta aysystem. nary 9 & i ta.2010 al. et Sim Heidelberg, − ost-impact h donor The 1 double 0 rogenitor making , nfuture in © alia . ; M 5 S 2021 ESO ose al. et Boos tellar ⊙ or of n for ons is mo- c- l- a - ; ; ; Draft version September 22, 2021 Typeset using LATEX default style in AASTeX631

Detecting preheating in proto-clusters with Lyman-α Forest Tomography

Robin Kooistra,1 Shigeki Inoue,2, 3, 4 Khee-Gan Lee,1 Renyue Cen,5 and Naoki Yoshida6, 1, 7

1Kavli IPMU (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 2Center for Computational Sciences, University of Tsukuba, Ten-nodai, 1-1-1 Tsukuba, Ibaraki 305-8577, Japan 3Chile Observatory, National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 4Department of Cosmosciences, Graduates School of Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan 5Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 6Department of Physics, School of Science, The University of Tokyo, Bunkyo, Tokyo 113-0033, Japan 7Institute for Physics of Intelligence, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan

ABSTRACT Studies of low redshift galaxy clusters suggest the intra-cluster medium (ICM) has experienced non-gravitational heating during the formation phase of the clusters. Using simple phenomenological heating prescriptions, we simulate the effect of this preheating of the nascent ICM in galaxy proto- clusters and examine its effect on Lyman-α (Lyα) forest tomographic maps. We analyse a series of cosmological zoom-in simulations of proto-clusters within the framework of the Lyα transmission-dark matter (DM) density distribution. We find that the more energy is injected into the proto-ICM at z = 3, the more the distribution at high DM density tilts towards higher Lyα transmission. This effect has been confirmed in both low-resolution simulations adopting a preheating scheme based on entropy floors, as well as in higher-resolution simulations with another scheme based on energy floors. The evolution of the slope of this distribution is shown to vary with redshift. The methodology developed here can be applied to current and upcoming Lyα forest tomographic survey data to help constrain feedback models in galaxy proto-clusters.

Keywords: Galaxy clusters (584) — Intracluster medium(858) — Lyman alpha forest(980)

1. INTRODUCTION According to the hierarchical structure formation scenario, galaxy clusters are expected to form through gravita- tional collapse, where galaxies, galaxy groups and small clusters merge to form increasingly more massive structures. However, even in early X-ray observations of galaxy clusters, it was found that the slope of the X-ray luminosity- temperature relation can be much steeper than the predicted slope of 2 expected from clusters evolving solely through self-similarity (Kaiser 1986; Edge & Stewart 1991; Markevitch 1998). In order to explain this steeper slope, self- similarity needs to be broken and the intracluster medium (ICM) gas needs to be heated early on by non-gravitational sources (e.g. Evrard & Henry 1991). This has become known in the field as ‘cluster preheating’. Support for the preheating scenario has additionally been found in measurements of the entropy in cluster cores (Pratt et al. 2010; Chaudhuri et al. 2013; Iqbal et al. 2017), suggesting that something is heating up the intra-cluster gas in the early arXiv:2109.09954v1 [astro-ph.GA] 21 Sep 2021 stages of the cluster formation, likely during the ‘Cosmic Noon’ epoch of peak star-formation and AGN activity at z ∼ 3.

The origin of the preheating of the ICM is still the subject of debate, and several possibilities have been explored. These include energy injection by star formation through supernova feedback, although the amount of energy super- novae can transport is likely insufficient (Loewenstein 2000; Menci & Cavaliere 2000; Bower et al. 2001). Therefore, the major contributor to the preheating is thought to be black hole feedback from active galactic nuclei (AGN; e.g. Valageas & Silk 1999; Wu et al. 2000; Chaudhuri et al. 2012). The preheating scenario has also been implemented in multiple numerical simulations of clusters (e.g. Navarro et al. 1995; Borgani et al. 2001), as well as the intergalactic medium (IGM; Borgani & Viel 2009), helping to narrow down the possible mechanisms that can lie behind the non-gravitational heat injection. Recently, a radio AGN has been found in a proto-cluster at redshift z = 3.3 (Shen et al. 2021), which adds to the growing consensus that black hole feedback already plays a significant role in the early A Horizon Study for Cosmic Explorer Science, Observatories, and Community arXiv:2109.09882v1 [astro-ph.IM] 20 Sep 2021

Contents

1 Executive Summary1

2 Purpose and Scope6

Science Objectives8

3 Overview9

4 Status of Ground-Based Gravitational-Wave Observatories 14

5 Key Science Questions 16 5.1 Black Holes and Neutron Stars Throughout Cosmic Time 16 5.1.1 Remnants of the First Stars 16 5.1.2 Seed Black Holes and Galaxy Formation 17 5.1.3 Formation and Evolution of Compact Objects 18 5.2 Dynamics of Dense Matter 19 5.2.1 Neutron Star Structure and Composition 20 5.2.2 New Phases in Quantum Chromodynamics 22 5.2.3 Chemical Evolution of the Universe 23 5.2.4 Gamma-Ray Burst Jet Engine 24 5.3 Extreme Gravity and Fundamental Physics 24 5.3.1 Nature of Strong Gravity 26 5.3.2 Unusual and Novel Compact Objects 27 5.3.3 Dark Matter and Dark Energy 27 5.4 Discovery Potential 29 5.4.1 Quantum Gravity 32 5.4.2 New Particles and Fields 33 5.4.3 Stochastic Gravitational-Wave Backgrounds 33

Observatories 36

6 A Science-Driven Design for Cosmic Explorer 37 6.1 Design Concept for Cosmic Explorer 37 6.2 Technology survey 38

iii The cosmic coincidences of primordial-black-hole dark matter

Yi-Peng Wua, Elena Pinettiab, and Joseph Silkcde aLaboratoire de Physique Th´eorique et Hautes Energies (LPTHE), UMR 7589 CNRS & Sorbonne Universit´e,4 Place Jussieu, F-75252, Paris, France bDipartimento di Fisica, Universit´adi Torino & INFN, Sezione di Torino, via P. Giuria 1, I-10125 Torino, Italy cInstitut d’Astrophysique de Paris, UMR 7095 CNRS & Sorbonne Universit´e,98 bis boulevard Arago, F-75014 Paris, France dDepartment of Physics and Astronomy, The Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, U.S.A. and eBeecroft Institute for Particle Astrophysics and Cosmology, University of Oxford, Keble Road, Oxford OX1 3RH, U.K. (Dated: September 22, 2021) If primordial black holes (PBHs) contribute more than 10 percent of the dark matter (DM) density, their energy density today is of the same order as that of the baryons. Such a cosmic coincidence might hint at a mutual origin for the formation scenario of PBHs and the baryon asymmetry of the Universe. Baryogenesis can be triggered by a sharp transition of the rolling rate of inflaton from slow-roll to (nearly) ultra-slow-roll phases that produces large curvature perturbations for PBH formation in single-field inflationary models. We show that the baryogenesis requirement drives the PBH contribution to DM, along with the inferred PBH mass range, the resulting stochastic gravitational wave background frequency window, and the associated cosmic microwave background tensor-to-scalar ratio amplitude, into potentially observable regimes.

I. INTRODUCTION II. PBHS FROM INFLATION

Primordial black holes (PBHs) are one of the most in- Let us focus on single-field inflation for PBH forma- teresting dark matter (DM) candidates that have been tion [39–55, 63–66]. The generic assumption is that the inflaton φ receives a sudden deceleration on comoving severely constrained by joint astrophysical and cosmo- 12 15 −1 scales k0 ∼ 10 −10 Mpc which leads to a sharp de- logical observations [1–12]. The asteroid-mass window 2 −16 −12 crease of the first slow-roll parameter H ≡ −H/H˙ and MPBH/M ∼ 10 − 10 for PBHs to be all DM [13, 14] will be tested in the near future by femto-lensing thus largely enhances the power spectrum of the curva- of gamma ray bursts [15–17], microlensing of X-ray pul- ture perturbation Pζ [56]. Such an enhancement is due sars [18], primary photons measured by next-generation to the temporal dominance of the entropy mode in the MeV detectors [7, 19], neutron star disruption [20], white curvature perturbation [54, 57, 68]. dwarf explosions [21], MeV photons from the Galactic The key parameter for realizing the enhancement of Pζ Center [22] and radio emission measurable with the next- is the rate-of-rolling: generation radio telescopes [23]. Even if PBHs only oc- cupy a small fraction of the DM density today, their exis- φ¨ δ = < −3/2, t0 < t < t∗, (1) tence could be (or could have been [24]) probed by grav- Hφ˙ itational wave observations through binary mergers, ei- ther mutual or with neutron stars (for recent reviews, see where δ → 0 is the standard slow-roll inflation and δ → [25, 26]). −3 is the so-called ultra-slow-roll (USR) limit [58–62]. If PBHs really constitute more than 10 percent of With δ < −3, Pζ can exhibit a spiky peak at a desired the DM density today, the PBH density up to matter- scale for producing a nearly monochromatic distribution arXiv:2109.09875v1 [astro-ph.CO] 20 Sep 2021 radiation equality is of the same order as that of the of PBH mass. baryons, namely ΩPBHeq/ΩBeq ∼ O(1). Such a cosmic We consider that δ takes constant values in different coincidence would hint at a mutual origin for PBHs and phases of inflation (see [46, 54] for some physical mod- the baryon asymmetry of the Universe. In this work, we els). In terms of the e-folding numbers, N ≡ ln a, where N∗−N0 argue that the PBH-baryon density coincidence could be k0/k∗ ≈ a(t0)/a(t∗) = e with N0 ≡ 0, the duration a natural consequence due to baryogenesis triggered by N∗ (for δ having a negative value) is the fundamental pa- inflation models for PBH formation. We show that the rameter that controls the spectral amplitude of Pζ (and abundance of PBHs and baryons are indirectly correlated thus the PBH abundance). to each other via the dynamics of inflation, and, thus, the The analytic structure of the USR power spectrum at scenario is very different from previous suggestions that the end of inflation (N = Nend) has been intensively in- existing PBHs create baryon asymmetry [27–35] or ac- vestigated [63–68]. Dilatation symmetry of the de Sitter count for cosmic coincidence [36–38]. background requires the momentum scaling at each phase Mode-by-mode Relative Binning: Fast Likelihood Estimation for Gravitational Waveforms with Spin-Orbit Precession and Multiple Harmonics

Nathaniel Leslie∗ and Liang Dai† Department of Physics University of California, Berkeley, CA 94720, USA

Geraint Pratten‡ School of Physics and Astronomy & Institute for Gravitational Wave Astronomy, University of Birmingham Birmingham, B15 2TT, UK (Dated: September 22, 2021) Faster likelihood evaluation enhances the efficiency of gravitational wave signal analysis. We present Mode-by-mode Relative Binning (MRB), a new method designed for obtaining fast and ac- curate likelihoods for advanced waveform models that include spin-orbit precession effects and multi- ple radiation harmonics from compact binary coalescence. Leveraging the “twisting-up” procedure of constructing precessing waveform modes from non-precessing ones, the new method mitigates degrade of relative binning accuracy due to interference from superimposed modes. Additionally, we supplement algorithms for optimizing the choice of frequency bins specific to any given strain signal under analysis. Using the new method, we are able to evaluate the likelihood with up to an order of magnitude reduction in the number of waveform model calls per frequency compared to the previously used relative binning scheme, and achieve better likelihood accuracy than is sufficient for obtaining source parameter posterior distributions that are indistinguishable from the exact ones.

I. INTRODUCTION time of a likelihood evaluation scales nearly linearly with the number of queried frequencies [1]. Astrophysical information about gravitational wave A number of methods have been developed to approx- sources is commonly extracted via parameter infer- imate the likelihood with fewer frequency evaluations, ence within the Bayesian framework. This process re- trading an insignificantly amount of accuracy degrade for quires thorough sampling of a high-dimensional param- substantial speedup. Reduced order quadrature (ROQ) eter space (e.g. with a total of 10–15 source parame- [1] is one such method, which pre-computes a reduced ters), and hence a huge number (typically on the order basis for a given waveform model, and approximates any of 107–108) of likelihood evaluations. To map either the physical waveform with an interpolant in the linear space likelihood function or the posterior distribution, it is im- spanned by this basis. This method only needs as many portant to accurately evaluate the likelihood function on independent evaluations as there are basis elements, and manageable timescales. has been applied to reduce the computational time of a Brute-force calculation of the likelihood is often costly. variety of waveform models [2–5]. In order to be compatible with Fast Fourier Trans- Relative binning [6–8], originally developed for non- form (FFT), it requires evaluating the frequency-domain precessing waveforms in the dominant (2, 2) radiation model waveform on a uniformly sampled, sufficiently fine mode, is another method to achieve the same goal. In grid that has ∼ T fmax frequencies, where fmax is the relative binning, the ratio between a given waveform and highest frequency that includes information and T is the an appropriately chosen fiducial waveform is considered. duration of the gravitational wave signal in the sensi- If the fiducial waveform is a decent fit of the data, only tive frequency band of the detector. The computational waveforms that resemble the fiducial one need to be ex- cost can be especially prohibitive if T is long or if fmax is amined in parameter inference, for which the waveform large. Even non-uniform frequency sampling still requires ratio conveniently has a milder frequency dependence. arXiv:2109.09872v1 [astro-ph.IM] 20 Sep 2021 a large number of grid points, as the waveform phase For this reason, the waveform ratio can be well approx- typically evolves over a large number of cycles across the imated as a piecewise-linear function in the frequency sensitive frequency band. For GW parameter inference domain, which only requires waveform model calls at the with sophisticated waveform models, model calls at indi- edges of a reduced number of frequency bins. For sim- vidual frequencies often dominate the runtime. For many ple waveform models with only one harmonic such as frequency domain waveform models, the computational IMRPhenomD [9], this method has been shown to require fewer waveform evaluations than ROQ for the same error tolerance [6], and has been routinely used in independent analysis of public LIGO/Virgo data [10–14]. An earlier ∗ nathaniel [email protected] and somewhat similar method referred to as heterodyne, † [email protected] which also samples the ratio between two waveforms us- ‡ [email protected] ing fewer grid points, was proposed in Ref. [15]. For submission to ApJ Letters FERMILAB-PUB-21-449-AE Draft version September 22, 2021 Typeset using LATEX default style in AASTeX631

C/2014 UN271 (Bernardinelli-Bernstein): the nearly spherical cow of comets Pedro H. Bernardinelli, Gary M. Bernstein, Benjamin T. Montet, Robert Weryk, Richard Wainscoat, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, S. Avila, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, R. Cawthon, C. Conselice, M. Costanzi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, H. T. Diehl, S. Everett, I. Ferrero, B. Flaugher, J. Frieman, J. Garc´ıa-Bellido, E. Gaztanaga, D. W. Gerdes, D. Gruen, R. A. Gruendl, J. Gschwend, G. Gutierrez, S. R. Hinton, D. L. Hollowood, K. Honscheid, D. J. James, K. Kuehn, N. Kuropatkin, O. Lahav, M. A. G. Maia, J. L. Marshall, F. Menanteau, R. Miquel, R. Morgan, R. L. C. Ogando, F. Paz-Chinchon,´ A. Pieres, A. A. Plazas Malagon,´ M. Rodriguez-Monroy, A. K. Romer, A. Roodman, E. Sanchez, M. Schubnell, S. Serrano, I. Sevilla-Noarbe, M. Smith, M. Soares-Santos, E. Suchyta, M. E. C. Swanson, G. Tarle, C. To, M. A. Troxel, T. N. Varga, A. R. Walker, Y. Zhang, (The DES Collaboration)

ABSTRACT

C/2014 UN271 (Bernardinelli-Bernstein) is a comet incoming from the Oort cloud which is remarkable in having the brightest (and presumably largest) nucleus of any well-measured comet, and having been discovered at heliocentric distance rh ≈ 29 au, farther than any Oort-cloud member. We describe in this work the discovery process and observations, and the properties that can be inferred from images recorded until the first reports of activity in June 2021. The orbit has i = 95◦, with perihelion of 10.97 au to be reached in 2031, and previous aphelion at 40, 400 ± 260 au. Backwards integration of the orbit under a standard Galactic tidal model and known stellar encounters suggests this is a pristine new comet, with a perihelion of q ≈ 18 au on its previous perihelion passage 3.5 Myr ago. The photometric data show an unresolved nucleus with absolute magnitude Hr = 8.0, colors that are typical of comet nuclei or Damocloids, and no secular trend as it traversed the range 34–23 au. For −0.5 r-band geometric albedo pr, this implies a diameter of 150(pr/0.04) km. There is strong evidence of brightness fluctuations at ±0.2 mag level, but no rotation period can be discerned. A coma, nominally consistent with a “stationary” 1/ρ surface-brightness distribution, grew in scattering cross-section at an exponential rate from Afρ ≈ 1 m to ≈ 150 m as the comet approached from 28 to 20 au. The activity rate is consistent with a very simple model of sublimation of a surface species in radiative equilibrium with the Sun. The inferred enthalpy of sublimation matches those of CO2 and NH3. More-volatile species such as N2,CH4, and CO must be far less abundant on the sublimating surfaces.

1. INTRODUCTION Our knowledge of the content of the Oort cloud is highly fragmentary—all inferences are based upon the small subset of its members that are torqued into orbits with perihelia q . 10 au, and until recently only the subset of these which develop comae bright enough to be noticed as comets. The cometary activity makes the objects easier to find and arXiv:2109.09852v1 [astro-ph.EP] 20 Sep 2021 makes it easier to identify the composition of the surface volatiles, but it can also obscure the properties of the nuclear body. The diversity of Oort cloud bodies has only recently begun to be explored, with the discovery of objects having varying levels of activity beyond the water frost line at ≈ 5 au (Meech et al. 2009; S´arneczkyet al. 2016; Jewitt et al. 2017, 2021; Meech et al. 2017; Hui et al. 2018, 2019). The discovery of C/2014 UN271 (Bernardinelli-Bernstein) (BB hereafter, for brevity) has expanded this known diversity substantially: as we will elaborate below, it is probably the largest Oort body ever found (indeed the largest of any kind of comet), and the first high-quality observations were taken when BB was at heliocentric distance rh ≈ 29 au in 2014, well before the first announced detection of coma

[email protected]

[email protected] Prepared for submission to JCAP

Light bending by the cosmological constant

Lingyi Hua Alan Heavensa David Baconb

aImperial Centre for Inference and Cosmology (ICIC), Imperial College, London, SW7 2AZ, U.K. bInstitute of Cosmology and Gravitation, University of Portsmouth, 1-8 Burnaby Road, Portsmouth PO1 3FX, U.K. E-mail: [email protected], [email protected], [email protected]

Abstract. We revisit the question of whether the cosmological constant Λ affects the cosmo- logical gravitational bending of light, by numerical integration of the geodesic equations for a Swiss cheese model consisting of a point mass and a compensated vacuole, in a Friedmann- Robertson-Walker background. We find that there is virtually no dependence of the light bending on the cosmological constant that is not already accounted for in the angular diam- eter distances of the standard lensing equations, plus small modifications that arise because the bending is restricted to a finite region covered by the hole. The residual Λ dependence for 13 7 a 10 M lens is at the level of 1 part in 10 , and even this might be accounted for by small changes in the hole size evolution as the photon crosses. We therefore conclude that there is no need for modification of the standard cosmological lensing equations in the presence of a cosmological constant. arXiv:2109.09785v1 [astro-ph.CO] 20 Sep 2021 DES-2015-0048 FERMILAB-PUB-21-385-SCD

Mon. Not. R. Astron. Soc. 000, 1–18 (0000) Printed 22 September 2021 (MN LATEX style file v2.2)

Finding quadruply imaged quasars with machine learning. I. Methods

A. Akhazhanov,1,2? A. More,3,4 A. Amini,5 C. Hazlett,5,6 T. Treu,7† S. Birrer,8 A. Shajib,9‡ P. Schechter,10 C. Lemon,11 B. Nord,12,13,9 M. Aguena,14 S. Allam,12 F. Andrade-Oliveira,15,14 J. Annis,12 D. Brooks,16 E. Buckley-Geer,9,12 D. L. Burke,17,18 A. Carnero Rosell,14 M. Car- rasco Kind,19,20 J. Carretero,21 A. Choi,22 C. Conselice,23,24 M. Costanzi,25,26,27 L. N. da Costa,14,28 M. E. S. Pereira,29,30 J. De Vicente,31 S. Desai,32 J. P. Dietrich,33 P. Doel,16 S. Everett,34 I. Ferrero,35 D. A. Finley,12 B. Flaugher,12 J. Frieman,12,13 J. Garc´ıa-Bellido,36 D. W. Gerdes,37,29 D. Gruen,38 R. A. Gruendl,19,20 J. Gschwend,14,28 G. Gutierrez,12 S. R. Hinton,39 D. L. Hollowood,34 K. Honscheid,22,40 D. J. James,41 A. G. Kim,42 K. Kuehn,43,44 N. Kuropatkin,12 O. Lahav,16 M. Lima,45,14 H. Lin,12 M. A. G. Maia,14,28 M. March,46 F. Menanteau,19,20 R. Miquel,47,21 R. Morgan,48 A. Palmese,12,13 F. Paz-Chinch´on,19,49 A. Pieres,14,28 A. A. Plazas Malag´on,50 E. Sanchez,31 V. Scarpine,12 S. Serrano,51,52 I. Sevilla-Noarbe,31 M. Smith,53 M. Soares-Santos,29 E. Suchyta,54 M. E. C. Swanson,19 G. Tarle,29 C. To,55,17,18 T. N. Varga,56,38 and J. Weller56,38 (DES Collaboration)

Accepted . Received

ABSTRACT Strongly lensed quadruply imaged quasars (quads) are extraordinary objects. They are very rare in the sky – only a few tens are known to date – and yet they provide unique information about a wide range of topics, including the expansion history and the composition of the Universe, the distribution of stars and dark matter in galaxies, the host galaxies of quasars, and the stellar initial mass function. Finding them in astronomical images is a classic “needle in a haystack” problem, as they are outnumbered by other (contaminant) sources by many orders of magnitude. To solve this problem, we develop state-of-the-art deep learning methods and train them on realistic simulated quads based on real images of galaxies taken from the Dark Energy Survey, with realistic source and deflector models, including the chromatic effects of microlensing. The performance of the best methods on a mixture of simulated and real objects is excellent, yielding area under the receiver operating curve in the range 0.86 to 0.89. Recall is close to 100% down to total magnitude i ∼ 21 indicating high completeness, while precision declines from 85% to 70% in the range i ∼ 17 − 21. The methods are extremely fast: training on 2 million samples takes 20 hours on a GPU machine, and 108 multi-band cutouts can be evaluated per GPU-hour. The speed and performance of the method pave the way to apply it to large samples of astronomical sources, bypassing the need for photometric pre-selection that is likely to be a major arXiv:2109.09781v1 [astro-ph.CO] 20 Sep 2021 cause of incompleteness in current samples of known quads. Key words: gravitational lensing: strong – methods: statistical – astronomical data bases: catalogs

1 INTRODUCTION

Strong gravitational lenses are extremely valuable sources of information about the Universe. They provide unique infor- ? e-mail:[email protected] mation about the expansion rate of the Universe, the prop- † Packard fellow erties of distant sources that would be too faint (compact) ‡ NHFP Einstein Fellow to be detected (resolved), and about the distribution of mass

© 0000 RAS KCL-PH-TH-2021-68

Footprints of population III stars in the gravitational-wave background

Katarina Martinovic,1 Carole P´erigois,2 Tania Regimbau,2 and Mairi Sakellariadou1 1Theoretical Particle Physics and Cosmology Group, Physics Department, King’s College London, University of London, Strand, London WC2R 2LS, UK 2LAPP, CNRS, 9 Chemin de Bellevue, 74941 Annecy-le-Vieux, France (Dated: September 2021) We investigate detection prospects of the gravitational-wave background (GWB) that originates from the merging of compact objects formed by the collapse of population III stars. Younger population I/II stars lead to a GWB in the LIGO/Virgo frequency band at the inspiral phase, while population III stars would likely show up at the later merger and ringdown phases. We show that, using a network of third-generation detectors, we may be able to separate a population I/II signal from a population III one, provided we can subtract individual coalescence events. A detection of a population III GWB could reveal important information, such as the average redshifted total mass.

Introduction—We have witnessed a rapid expansion signal due to different mass and redshift distributions [25]. of gravitational-wave (GW) astrophysics in the last decade The GWB is comprised of many sources, of astrophysical due to the success of the Advanced LIGO and Virgo or cosmological origin, but we expect the CBC signal to GW detectors [1, 2] at uncovering signals from numerous be the foreground to all sources [26]. In this study, for the compact binary coalescences (CBCs) [3]. The most recent first time, we consider the possibility of separating pop Advanced LIGO/Virgo observing run, O3, presented us I/II and pop III GWB contributions. Numerous models with dozens of new merger events and has significantly suggest the total CBC background is dominated by pop expanded the stellar graveyard [4]. Despite the increase I/II. However, pop III can be uncovered using subtrac- in the number of detections, we are yet to observe with tion techniques and studying the residual backgrounds confidence an event that would suggest that the progenitor [27, 28]. As the sensitivity of detectors increases and GW compact objects are remnants of the oldest stars in the interferometers see more individual CBC events, a pop Universe [5] - the theoretically-postulated population III III residual background emerges as the dominant signal (pop III) stars (GW190521 could be a potential candidate over pop I/II residual background. We first study how to [6]). Pop III stars are thought to have formed at high detect the GWB from pop III stars, and in the case of a redshifts and as such have low metallicity compared to successful detection, we explore subsequent implications, the more recently formed, population I/II (pop I/II) stars namely information about masses and redshifts of the [7–9]. These old stars have hitherto evaded sky surveys [10– population. 12], and their detection remains an objective for upcoming The paper is organised as follows: in Sec. II we present experiments, such as the James Webb Space Telescope the pop III models and their resulting GWB in different [13]. detector networks, highlighting characterisation of the re- Pop III stars may solve some of the puzzles in black sidual backgrounds. We then introduce Bayesian analysis hole formation, as well as help understanding the early used in Sec. III, and discuss search filters we consider for epochs of the Universe such as reionisation and galaxy pop III stars. Sec. IV is an implications study in case of evolution [14–16]. Numerical simulations show that these a detection of a pop III signal. We use StarTrack (ST) primordial stars could have led to the formation of super- simulation data [29] and apply our detection methods, ul- massive black holes at high redshifts [17–19]. Mergers timately showing consistency of our implications analysis of such heavy remnants would appear in the millihertz with the underlying population. We select the ST data frequency range explored by future space-based detect- since it is the most recent extensive catalogue of merging arXiv:2109.09779v1 [astro-ph.SR] 20 Sep 2021 ors such as LISA. The scope of this study, however, is binaries from pop I, II and III stars that lead to a GWB detection prospects of terrestrial detector networks, and in the LIGO/Virgo frequency range [30]. we therefore focus on models that predict a pop III signal Population III GWB— The GWB is defined as the in the LIGO/Virgo frequency range. The contribution superposition of GWs from all unresolved sources. It 1 to the gravitational-wave background (GWB) from a su- is characterised by the dimensionless parameter ΩGW perposition of unresolved pop III-seeded CBCs has been [31], expressed as the ratio of GW energy density per explored in several studies [21–25]. They show significant logarithmic frequency bin dρGW/d ln(f), normalised deviation of a pop III star signals from a pop I/II stars by the critical energy density of the Universe, ρc = 2 2 (3H0 c )/(8πG):

1 Often referred to as the stochastic gravitational-wave background 1 dρGW(f) [20]. ΩGW(f) = , (1) ρc dln(f) Draft version September 22, 2021 Typeset using LATEX twocolumn style in AASTeX631

Tidal Distortions in NGC1052-DF2 and NGC1052-DF4: Independent Evidence for a Lack of Dark Matter

Michael A. Keim ,1 Pieter van Dokkum ,1 Shany Danieli ,2, 3, ∗ Deborah Lokhorst ,4, 5 Jiaxuan Li (李嘉轩) ,2 Zili Shen ,1 Roberto Abraham ,4, 5 Seery Chen ,4, 5 Colleen Gilhuly ,4, 5 Qing Liu (刘青) ,4, 5 Allison Merritt ,6 Tim B. Miller ,1 Imad Pasha ,1 and Ava Polzin 1

1Department of Astronomy, Yale University, PO Box 208101, New Haven, CT 06520-8101, USA 2Department of Astrophysical Sciences, 4 Ivy Lane, Princeton University, Princeton, NJ 08544 3Institute for Advanced Study, 1 Einstein Drive, Princeton, NJ 08540, USA 4Department of Astronomy & Astrophysics, University of Toronto, 50 St. George St., Toronto, ON M5S 3H4, Canada 5Dunlap Institute for Astronomy and Astrophysics, University of Toronto, Toronto ON, M5S 3H4, Canada 6Max-Planck-Institut f¨urAstronomie, K¨onigstuhl17, D-69117 Heidelberg, Germany

ABSTRACT Two ultra diffuse galaxies in the same group, NGC1052-DF2 and NGC1052-DF4, have been found to have little or no dark matter and to host unusually luminous globular cluster populations. Such low mass diffuse objects in a group environment are easily disrupted and are expected to show evidence of tidal distortions. In this work we present deep new imaging of the NGC1052 group obtained with the Dragonfly Telephoto Array to test this hypothesis. We find that both galaxies show strong position angle twists and are significantly more elongated in their outskirts than in their interiors. The group’s central massive elliptical NGC1052 is the most likely source of these tidal disturbances. The observed distortions imply that the galaxies have a very low total mass; otherwise, they must be very close to NGC1052. Taking into account a recently derived lower limit on the relative line-of-sight distance between NGC1052-DF2 and NGC1052-DF4 from the tip of the red giant branch stars, we infer that the dark matter halo masses of these galaxies cannot be much greater than their stellar masses. Our findings provide strong evidence, completely independent of kinematic constraints, that both galaxies are indeed dark matter deficient. While tidal distortions do not directly provide new information on the formation of NGC1052-DF2 and NGC1052-DF4, the similarity of their tidal features suggests that they are at comparable distances from NGC1052.

Keywords: Dark matter (353) — Galaxy evolution (594) — Galaxy structure (622) — Low surface brightness galaxies (940) — Tidal distortion (1697)

1. INTRODUCTION with the massive elliptical NGC1052 (Shen et al. 2021b) The galaxies NGC1052-DF2 and NGC1052-DF4 both and fall into the category of ultra diffuse galaxies; large, have little to no dark matter based on the velocity dis- faint objects with the sizes of L∗ galaxies but the stellar persion of kinetic tracers (van Dokkum et al. 2018a; masses of dwarf galaxies (van Dokkum et al. 2015). arXiv:2109.09778v1 [astro-ph.GA] 20 Sep 2021 Wasserman et al. 2018; van Dokkum et al. 2019; Danieli Both their luminous globular clusters and low dark et al. 2019; Emsellem et al. 2019) and host a unique matter content set NGC1052-DF2 and NGC1052-DF4 population of globular clusters nearly two magnitudes apart from almost all other galaxies known. The near- brighter than expected (van Dokkum et al. 2018b; Ma universal globular cluster luminosity function has a et al. 2020; Shen et al. 2021a). They are likely associated peak at MV ≈ −7.5 mag (Rejkuba 2012), whereas the NGC1052-DF2 and NGC1052-DF4 distributions peak at MV ≈ −9 mag (Shen et al. 2021a). Galaxies with the Corresponding author: Michael A. Keim stellar mass of NGC1052-DF2 and NGC1052-DF4 (ap- [email protected] 8 8 proximately 2.0×10 M and 1.5×10 M , respectively; ∗ NASA Hubble Fellow van Dokkum et al. 2018a, 2019) ought to have a >300 arXiv:2109.09776v1 [astro-ph.EP] 20 Sep 2021 ★ © ( ytm ( systems MNRAS zkl tal. et Czekala 2020 al. et Kraus etayselrcmain()a rjce itneo ab of a distance and projected au a 1 about at ( of (C) au 8 separation companion a stellar at tertiary (A,B) binary spectroscopic discovered. been has orbit circumtriple a in planet no olntr ici hc h uti bevdt xedt ab to extend to observed is dust the which in disc toplanetary 2 of masses 2014a ( Bi Jiaqing ob 0 be to eeyL Smallwood, L. Jeremy planets and rings circumtriple Ori: GW f1 of ar hl bu 0 r bevdi rpeo higher-order or triple in observed ( are b tems a 20% in about are stars while of 40%-50% pair, than mult more to that estimated leading is It naturally tems. clusters, stellar relativ within in form regions stars of majority the that show Observations INTRODUCTION 1 cetdXX eevdYY noiia omZZZ form original in YYY; Received XXX. Accepted 2 1 3 6 5 4 7 avte l 2004 al. et Calvet 2020 al. et Collaboration Gaia eateto hsc n srnm,Uiest fLeicest of University Astronomy, and Physics of Department eateto hsc n srnm,Uiest fNevada, of University Astronomy, and Physics of Department etefrEolnt n aiaiiy nvriyo Warw of University Habitability, and Exoplanets for Centre colo hsc n srnm,Mns nvriy Clayto University, Monash Astronomy, and Physics of School eateto hsc n srnm,Uiest fVictori of University Astronomy, and Physics of Department eateto hsc,Uiest fWrik oetyCV4 Coventry Warwick, of University Physics, of Department nv rnbeAps NS PG -80 rnbe France Grenoble, F-38000 IPAG, CNRS, Alpes, Grenoble Univ. 00TeAuthors The 2020 -al [email protected] E-mail: . M 2 WOii n uhheacia rpesa ytm iha with system, star triple hierarchical such one is Ori GW ahe ta.1991 al. et Mathieu , uuno ao 1991 Mayor & Duquennoy . b 2( 22 .T ae hr aebe lnt on n3 rpestar triple 32 in found planets been have there date, To ). 000 ⊙ uet ta.2018 al. et Busetti n h cetiiyo h ue opno sestimated is companion outer the of eccentricity the and zkl ta.2017 al. et Czekala . , M 8 1 – ( 16 2017 ⊙ 5 .Tesse sa itneo 408 of distance a at is system The ). (2020) ubn Dong Ruobing , n 1 and .Tedrvdtil trprmtr from parameters star triple derived The ). r hw nTable in shown are ) . M 7 ; egre l 2011 al. et Berger ; ⊙ rgoee l 2019 al. et Fragione epciey hl trChsamass a has C star while respectively, , .Tesse a icmrpepro- circumtriple a has system The ). n a naeo 1 of age an has and ) ; ahvne l 2010 al. et Raghavan 1 asdb neetdpaes—tefis lnts nacircum disc a words: in the pla Key planet(s) that massive first conclude the a — We of planets disc. undetected presence outer by caused the and inner that the suggest We separates doe system. disc star the the discs, protoplanetary of parameters simulations motivated hydrodynamic 3-dimensional and integrations niiul WOri GW individual: n eoemslge.W netgt h rgn ftedus pr the nodally rings of the origins th broken, the planet is investigate disc a We the or misaligned. Once system become disc. star and in the triple break in the A gap from rings. a torque the ri the misaligned of by dust each a either between three with misalignments identified system and have au star 100 observations triple ALMA hierarchical Recent a disc. is Ori GW ABSTRACT ★ eec Nealon Rebecca 2 5 tr n have B and A Stars . ; n hitpePinte Christophe and zkl ta.2017 al. et Czekala .Hwvr syet as However, ). crto,aceindss–hdoyais–pae-ici planet-disc – hydrodynamics – discs accretion accretion, . c,Cvnr V A,UK 7AL, CV4 Coventry ick, 0 ,Vcoi,B 8 A,Canada 1A1, V8P BC Victoria, a, r nvriyRa,LietrLI7H UK 7RH, LEI Leicester Road, University er, ± ; i 80 Australia 3800, Vic n ± a ea,40 ot ayadPrwy a ea,N 89154, NV Vegas, Las Parkway, Maryland South 4505 Vegas, Las Tokovinin -trsys- i-star A,UK 7AL, 0parsec 10 l dense ely 0 . rpit2 etme 01Cmie sn NA L MNRAS using Compiled 2021 September 22 Preprint Myr 1 inary sys- out out ; 2 , 3 , 4 hn Chen Cheng , ne n idernsi h ds raig hnmnn wh phenomenon, breaking’ ‘disc the is rings middle and inner h ig n r iaindbtcoet olnrwt eac with coplanar to close but misaligned are wh 3 coplanar other. being and to 2 close are rings plane the binary sug AB-C results the and their 1 Therefore, ring plane. orbital the to relative 245 n eemslge y11 by misaligned were 3 and orbit binary AB-C the of 150 2011 node be to ascending found of inclinat was The longitude kinematics. the CO and and continuum from dust plane orbital of star modelling triple and ring 74 dust be individual each to estimated were 3 and 2, 1, dust The rings disc. the protoplanetary of a in found ever ring widest the ( au 1300 about to farther extends 2017 gas the while au, 400 rjce ria direction, orbital projected ai faot4,18 n 3 ufo h rpesa.Testr The star. Table in triple outlined the are from rings au the dust 338 d of and circumtriple 188, the 46, the in about of rings radii dust di three of circumtriple identified the They of Ori. observations GW emission gas molecular (ALMA) and continuum Array ter/submillimeter ± oeta xlnto o h ag iainetbetween misalignment large the for explanation potential A .Asmn htteetr ichstesm lcws on-sk clockwise same the has disc entire the that Assuming ). al. et Bi ). 6 8M 28 , 7 ⊕ ( ie l 2020 al. et Bi ( 2020 1 rsne tcm ag Millime- Large Atacama presented ) ± eec .Martin G. Rebecca , 7 ° n 28 and ie al. et Bi o ra u otetru from torque the to due break observationally- not for s that find We . ± .Mslgmnsaepeetbetween present are Misalignments ). 6 rpeorbit. triple ° 1 35 , ti asv nuht carve to enough massive is at h ailwdho utrn is 3 ring dust of width radial The . raigi WOii likely is Ori GW in breaking ig ymasof means by rings t g ihapoietgpat gap prominent a with ngs cs ndffrn timescales different on ecess ± icmrpeprotoplanetary circumtriple e o lnt)i h disc the in planets) (or net h a icmyb driven be may disc gas the ± ( 9 2020 ° 5 epciey( respectively , ° n 40 and USA on htrns1 2, 1, rings that found ) trcin stars: – nteractions ± ± A ,168 8, T 5 E ° tl l v3.0 file style X respectively, , egre al. et Berger 1 etdthat gested age al. et Fang , ± # visibility 9 and 19, masses –body ucture s at isc cof sc ion ere the ile h y Research in Astron. Astrophys. 2018 Vol. X No. XX, 000–000 http://www.raa-journal.org http://www.iop.org/journals/raa Research in (LATEX: msRAA-2019-0285.tex; printed on September 22, 2021; Astronomy and 0:35) Astrophysics

The Multiplicity of Early-type Stars from LAMOST Medium-resolution Spectroscopic Survey

Yanjun Guo1,2, Jiao Li3, Jianping Xiong3, Jiangdan Li1,2, Luqian Wang1, Heran Xiong4, Feng Luo3, Yonghui Hou5,6, Chao Liu2,3, Zhanwen Han1,2,7 and Xuefei Chen1,2,7

1 Yunnan observatories, Chinese Academy of Sciences, P.O. Box 110, Kunming, 650011, China 2 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing, 100049, People’s Republic of China 3 Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, People’s Republic of China 4 Research School of Astronomy and Astrophysics, Mount Stromlo Observatory, The Australian National University, ACT 2611, Australia 5 Nanjing Institute of Astronomical Optics, & Technology, National Astronomical Observatories, Chinese Academy of Sciences, Nanjing 210042, China 6 School of Astronomy and Space Science, University of Chinese Academy of Sciences 7 Center for Astronomical Mega-Science, Chinese Academy of Sciences, Beijing 100012, China

Abstract Massive binaries play significant roles in many fields. Identification of massive stars, particularly massive binaries, is of great importance. In this paper, by adopting the technique of measuring the equivalent widths of several spectral lines, we identified 9,382 early-type stars from LAMOST medium-resolution survey and divided the sample into four groups, T1 (∼O-B4), T2 (∼B5-B7), T3 (∼B5-B7), and T4 (∼B8-A). The relative radial velocities RVrel were calculated using the Maximum Likelihood Estimation, and the stars with significant changes of RVrel were identified as spectroscopic binaries. We found that the observed spectroscopic binary fractions for the four groups are 24.6% ± 0.5%, 20.8%±0.6%, 13.7%±0.3%, and 7.4%±0.3%, respectively. Assuming that orbital period (P) and mass ratio (q) have intrinsic distributions as f (P) ∝ Pπ (1

Key words: Stars: early-type — Stars: spectroscopic binary — Stars: statistics — Catalogs MNRAS 000,1–21 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

NGC 5746: formation history of a massive disc-dominated galaxy

Marie Martig,1★, Francesca Pinna2, Jesús Falcón-Barroso3,4, Dimitri A. Gadotti5, Bernd Husemann2, Ivan Minchev6, Justus Neumann7, Tomás Ruiz-Lara8, Glenn van de Ven9 1Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK 2Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117 Heidelberg, Germany 3Instituto de Astrofísica de Canarias, Calle Via Láctea s/n, 38200 La Laguna, Tenerife, Spain 4Depto. Astrofísica, Universidad de La Laguna, Calle Astrofísico Francisco Sánchez s/n, 38206 La Laguna, Tenerife, Spain 5 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei München, Germany 6Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany 7 Institute of Cosmology and Gravitation, University of Portsmouth, Burnaby Road, Portsmouth, PO1 3FX, UK 8Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands 9Department of Astrophysics, University of Vienna, Türkenschanzstrasse 17, 1180 Wien, Austria

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT The existence of massive galaxies lacking a classical bulge has often been proposed as a challenge to ΛCDM. However, recent simulations propose that a fraction of massive disc galaxies might have had very quiescent merger histories, and also that mergers do not necessarily build classical bulges. We test these ideas with deep MUSE observations of NGC 5746, a massive 11 (∼ 10 M ) edge-on disc galaxy with no classical bulge. We analyse its stellar kinematics and stellar populations, and infer that a massive and extended disc formed very early: 80% of the galaxy’s stellar mass formed more than 10 Gyr ago. Most of the thick disc and the bar formed during that early phase. The bar drove gas towards the center and triggered the formation of the nuclear disc followed by the growth of a boxy/peanut-shaped bulge. Around ∼ 8 Gyr ago, a ∼1:10 merger happened, possibly on a low-inclination orbit. The satellite did not cause significant vertical heating, did not contribute to the growth of a classical bulge, and did not destroy the bar and the nuclear disc. It was however an important event for the galaxy: by depositing its stars throughout the whole galaxy it contributed ∼ 30% of accreted stars to the thick disc. NGC 5746 thus did not completely escape mergers, but the only relatively recent significant merger did not damage the galaxy and did not create a classical bulge. Future observations will reveal if this is representative of the formation histories of massive disc galaxies. Key words: galaxies: evolution, galaxies: formation, galaxies: individual: NGC 5746

1 INTRODUCTION Milky Way). One way to form thick discs is through mergers: mergers can thicken pre-existing discs (Quinn et al. 1993; Read et al. 2008; Disc galaxies host a variety of components, and the relation between Villalobos & Helmi 2008) and/or add accreted stars in a thick con- their detailed morphology and their formation history is complex, figuration (Abadi et al. 2003). Alternatively, discs can form already and still not fully understood. The general picture of disc galaxy thick from turbulent gas, either following gas-rich mergers (Brook formation starts with an early violent phase, with mergers or violent et al. 2004) or violent disc instabilities in gas-rich discs (Bournaud internal instabilities creating thick discs and bulges. Then, at lower arXiv:2109.09770v1 [astro-ph.GA] 20 Sep 2021 et al. 2009). Thick discs typically have older, more metal-poor pop- redshift, as the merger rate decreases, thinner discs can settle, form ulations than thin discs, although there exists some diversity (Mould bars and spiral arms and are dominated by secular evolution (slow 2005; Rejkuba et al. 2009; Yoachim & Dalcanton 2008b; Comerón changes driven by internal processes, see Kormendy & Kennicutt et al. 2015, 2016; Kasparova et al. 2016, 2020; Pinna et al. 2019a,b; 2004; Sellwood 2014). Those mechanisms all act to shape galaxies, Scott et al. 2021). In particular, some thick discs are complex com- and create the diversity of morphologies observed along the Hubble ponents, with their outer parts being formed from flared younger Sequence. populations (Minchev et al. 2015; García de la Cruz et al. 2021), Most spiral galaxies, including the Milky Way, contain both a thin creating radial age gradients like the one observed in the Milky Way and thick disc (Gilmore & Reid 1983; Dalcanton & Bernstein 2002; (Martig et al. 2016). Whether this picture applies to thick discs in Yoachim & Dalcanton 2006; Comerón et al. 2018), although it is general is unknown, but it is likely that thick discs could arise from not always clear if they are clearly meaningful distinct components several different formation channels. (see discussions in Bovy et al. 2012a; Mackereth et al. 2017 for the Bulges are also complex, and are usually classified in two main categories: classical and pseudobulges (see reviews by Kormendy & ★ E-mail: [email protected] Kennicutt 2004; Kormendy 2013; Fisher & Drory 2016). Classical

© 2021 The Authors DRAFTVERSION SEPTEMBER 22, 2021 Typeset using LATEX twocolumn style in AASTeX63

Hyper Suprime-Cam Subaru Strategic Program: A Mass-Dependent Slope of the Galaxy Size-Mass Relation at z < 1∗

† ‡ LALITWADEE KAWINWANICHAKIJ,1 , JOHN D.SILVERMAN,1, 2 XUHENG DING,1 ANGELO GEORGE,3 IVANA DAMJANOV,3 , ‡ MARCIN SAWICKI,3 , MASAYUKI TANAKA,4, 5 DAN S.TARANU,6 SIMON BIRRER,7 SONG HUANG,8 JUNYAO LI,9, 10, 1 MASATO ONODERA,11, 12 TAKATOSHI SHIBUYA,13 AND NAOKI YASUDA1

1Kavli Institute for the Physics and Mathematics of the Universe, The University of Tokyo, Kashiwa, Japan 277-8583 (Kavli IPMU, WPI) 2Department of Astronomy, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 3Institute for Computational Astrophysics and Department of Astronomy & Physics, Saint Mary’s University, 923 Robie Street, Halifax, Nova Scotia, B3H 3C3, Canada 4National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 5Department of Astronomy, School of Science, Graduate University for Advanced Studies (SOKENDAI), 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 6Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA 7Kavli Institute for Particle Astrophysics and Cosmology and Department of Physics, Stanford University, Stanford, CA 94305, USA 8Department of Astronomy and Astrophysics, University of California Santa Cruz, 1156 High St., Santa Cruz, CA 95064, USA 9CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei 230026, China 10School of Astronomy and Space Science, University of Science and Technology of China, Hefei 230026, China 11Subaru Telescope, National Astronomical Observatory of Japan, National Institutes of Natural Sciences (NINS), 650 North A’ohoku Place, Hilo, HI 96720, USA 12Department of Astronomical Science, The Graduate University for Advanced Studies, SOKENDAI, 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan 13Kitami Institute of Technology, 165 Koen-cho, Kitami, Hokkaido 090-8507, Japan

ABSTRACT

We present the galaxy size-mass (Re −M∗) distributions using a stellar-mass complete sample of 1.5 million 2 ∼ galaxies, covering 100 deg , with log(M∗/M ) > 10.2 (9.2) over the redshift range 0.2 < z < 1.0 (z < 0.6) ∼ from the second public data release of the Hyper Suprime-Cam Subaru Strategic Program. We confirm that, at fixed redshift and stellar mass over the range of log(M∗/M ) < 11, star-forming galaxies are on average larger than quiescent galaxies. The large sample of galaxies with accurate size measurements, thanks to the excellent imaging quality, also enables us to demonstrate that the Re − M∗ relations of both populations have a form of broken power-law, with a clear change of slopes at a pivot stellar mass Mp. For quiescent galaxies, 0.1 below an (evolving) pivot mass of log(Mp/M ) = 10.2 − 10.6 the relation follows Re M∗ ; above Mp the 0.6−0.7 ∝ relation is steeper and follows Re M∗ . For star-forming galaxies, below log(Mp/M ) 10.7 the relation 0.2 ∝ 0.3−0.6 ∼ follows Re M ; above Mp the relation evolves with redshift and follows Re M . The shallow power- ∝ ∗ ∝ ∗ law slope for quiescent galaxies below Mp indicates that large low-mass quiescent galaxies have sizes similar to those of their counterpart star-forming galaxies. We take this as evidence that large low-mass quiescent galaxies have been recently quenched (presumably through environment-specific process) without significant structural transformation. Interestingly, the pivot stellar mass of the Re − M∗ relations for both populations also coincides with mass at which half of the galaxy population is quiescent, implied that the pivot mass represents the transition of galaxy growth from being dominated by in-situ star formation to being dominated by (dry) arXiv:2109.09766v1 [astro-ph.GA] 20 Sep 2021 mergers.

Keywords: galaxies: evolution – galaxies: morphology

1. INTRODUCTION

Corresponding author: Lalitwadee Kawinwanichakij The size distribution and its evolution with cosmic time [email protected] provide important clues about assembly history of galaxies and the relationships with the dark matter halos in which ∗ Released on September 22, 2021 they reside (e.g., Mo et al. 1998; Kravtsov 2013; Shibuya † LSSTC Data Science Fellow ‡ Canada Research Chair et al. 2015; Huang et al. 2017; Somerville et al. 2018). The Draft version September 22, 2021 Typeset using LATEX default style in AASTeX63

On the Three-Dimensional Structure of Local Molecular Clouds Catherine Zucker,1 Alyssa Goodman,1 Joao˜ Alves,2 Shmuel Bialy,1 Eric W. Koch,3, 1 Joshua S. Speagle,4, 5, 6, 1 Michael M. Foley,1 Douglas Finkbeiner,1 Reimar Leike,7, 8 Torsten Enßlin,7, 8 Joshua E.G. Peek,9, 10 and Gordian Edenhofer7, 8

1Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., Cambridge, MA, USA 02138 2University of Vienna, Department of Astrophysics, T¨urkenschanzstraße 17, 1180 Vienna, Austria 3University of Alberta, Department of Physics, 4-183 CCIS, Edmonton AB T6G 2E1, Canada 4Department of Statistical Sciences, University of Toronto, Toronto, ON M5S 3G3, Canada 5David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, Toronto, ON M5S 3H4, Canada 6Dunlap Institute for Astronomy & Astrophysics, University of Toronto, Toronto, ON M5S 3H4, Canada 7Max Planck Institute for Astrophysics, Karl-Schwarzschildstraße 1, 85748 Garching, Germany 8Ludwig-Maximilians-Universit¨at,Geschwister-Scholl Platz 1, 80539 Munich, Germany 9Space Telescope Science Institute, 3700 San Martin Dr, Baltimore, MD 21218, USA 10Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA

Submitted to ApJ

ABSTRACT We leverage the 1 pc spatial resolution of the Leike et al.(2020) 3D dust map to characterize the three-dimensional structure of nearby molecular clouds (d . 400 pc). We start by “skeletonizing” the clouds in 3D volume density space to determine their “spines,” which we project on the sky to constrain cloud distances with ≈ 1% uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well-described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is ≈ 3 : 1. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, n, falling off like n−2 at large radii. Using Bayesian model selection, we find that parameterizing the clouds’ profiles using a single Gaussian is disfavored. We compare our results with 2D dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction (AV & 2−3 mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium.

Keywords: molecular clouds, solar neighborhood, interstellar reddening, astronomy data visualization arXiv:2109.09765v1 [astro-ph.GA] 20 Sep 2021 1. INTRODUCTION Much of our understanding of the structure and properties of molecular clouds is obtained via either two-dimensional (2D) observations of dust in emission or extinction (see Andr´eet al. 2010; Lombardi 2009; Planck Collaboration et al. 2011) or three-dimensional (3D) spectral-line observations of gas (see Dame et al. 2001), where the third dimension is not distance, but radial velocity. Until very recently, any insights into the true 3D spatial structure of molecular clouds were available only in numerical simulations, and were limited by the extent to which these simulations could

Corresponding author: Catherine Zucker [email protected] MNRAS 000,1–22 (2021) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0

Outflows from starburst galaxies with various driving mechanisms and their X-ray properties

B. P. Brian Yu1,2★, Ellis R. Owen1,3, Kuo-Chuan Pan1,3, Kinwah Wu2,4 and Ignacio Ferreras5,6,7 1Institute of Astronomy, National Tsing Hua University, Hsinchu, Taiwan (ROC) 2Mullard Space Science Laboratory, University College London, Dorking, Surrey, RH5 6NT, United Kingdom 3Center for Informatics and Computation in Astronomy, National Tsing Hua University, Hsinchu, Taiwan (ROC) 4Research Center for Astronomy, Astrophysics and Astrophotonics, Macquarie University, Sydney, NSW 2019, Australia 5Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom 6Instituto de Astrofísica de Canarias, C/Vía Láctea, s/n, E38205 La Laguna, Tenerife, Spain 7Departamento de Astrofísica, Universidad de La Laguna (ULL), E-38206 La Laguna, Tenerife, Spain

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Outflows in starburst galaxies driven by thermal-mechanical energy, cosmic rays and their mix are investigated with 1D and 2D hydrodynamic simulations. We show that these outflows could reach a stationary state, after which their hydrodynamic profiles asymptotically approach previous results obtained semi-analytically for stationary outflow configurations. The X-rays from the simulated outflows are computed, and high-resolution synthetic spectra and broadband light curves are constructed. The simulated outflows driven by thermal mechanical pressure and CRs have distinguishable spectral signatures, in particular, in the sequence of the keV K훼 lines of various ions and in the L-shell Fe emission complex. We demonstrate that broadband colour analysis in X-rays is a possible alternative means to probe outflow driving mechanisms for distant galaxies, where observations may not be able to provide sufficient photons for high-resolution spectroscopic analyses. Key words: galaxies: starburst – hydrodynamics – ISM: jets and outflows – X-rays: galaxies – cosmic rays – methods: numerical

1 INTRODUCTION et al. 2002a; Rupke et al. 2005b; Rubin et al. 2014). Outflow extents may reach a few tens of kpc (Veilleux et al. 2005; Bland-Hawthorn Large-scale outflows (winds) in galaxies are often initiated by star- et al. 2007; Bordoloi et al. 2011; Martin et al. 2013; Rubin et al. forming activities, e.g. in the nearby starburst galaxies Arp 220 (e.g. 2014; Bordoloi et al. 2016), often with X-ray emitting gas at lower Varenius et al. 2016; Barcos-Muñoz et al. 2018), M82 (e.g. Shopbell altitudes bound by a cap a few kpc from the starburst nucleus (Strick- & Bland-Hawthorn 1998; Bland & Tully 1988) and NGC 253 (e.g. land et al. 2000; Cecil et al. 2002b,a; Devine & Bally 1999; Tsuru Walter et al. 2017; Krieger et al. 2019), while more distant, younger et al. 2007). The gas temperature in some outflows may reach 107K galaxies are also known to host such phenomena (Frye et al. 2002; (McKeith et al. 1995; Shopbell & Bland-Hawthorn 1998). There

arXiv:2109.09764v1 [astro-ph.GA] 20 Sep 2021 Rubin et al. 2014; Rupke et al. 2005a; Arribas et al. 2014). These is evidence of multi-phase structure in galactic winds, with cooler outflows develop from the confluence of stellar winds from young dense clumps, of temperature 102 − 104 K (Strickland et al. 1997; stars, and energetic particles and gas from supernova events. Galac- Lehnert et al. 1999), entrained within the hot gas (see Heckman tic outflows tend to have a bi-conical structure, extending above and 2003; Zhang et al. 2014; Walter et al. 2017; Wu et al. 2020). The below the plane of their host galaxy, following the path of least re- presence of this multi-phase structure implies the coexistence of sistance through the interstellar environment (Veilleux et al. 2005). material with very different structural and local thermal conditions. The structures of outflows, however, vary among galaxies, show- The diversity of the physical conditions within an outflow naturally ing a variety of subtleties in morphology and dynamics. Velocities gives rise to a multitude of radiative processes, and this complexity of galactic outflows are mostly between 10 − 1000 km s−1 (Cecil is not often present in more uniform astrophysical systems. Non-AGN associated outflows are a consequence of star- ★ E-mail: [email protected] (BPBY), [email protected] forming activities in galaxies. Such outflows carry away energy, (ERO), [email protected] (KCP), [email protected] gas and metals from star-forming regions to the circum-galactic (KW) medium of their host galaxy. This regulates the inflow of circum-

© 2021 The Authors Preprint typeset using LATEX style emulateapj v. 12/16/11

THE PER-TAU SHELL: A GIANT STAR-FORMING SPHERICAL SHELL REVEALED BY 3D DUST OBSERVATIONS

SHMUEL BIALY1?,CATHERINE ZUCKER1,ALYSSA GOODMAN1,2,MICHAEL M.FOLEY1,JOAO˜ ALVES2,3,VADIM A.SEMENOV1, ROBERT BENJAMIN4,REIMAR LEIKE5,6,TORSTEN ENßLIN5,6

ABSTRACT A major question in the field of star formation is how molecular clouds form out of the diffuse Interstellar Medium (ISM). Recent advances in 3D dust mapping are revolutionizing our view of the structure of the ISM. Using the highest-resolution 3D dust map to date, we explore the structure of a nearby star-forming region, which includes the well-known Perseus and Taurus molecular clouds. We reveal an extended near-spherical shell, 156 pc in diameter, hereafter the “Per-Tau Shell”, in which the Perseus and Taurus clouds are embedded. We also find a large ring structure at the location of Taurus, hereafter, the “Tau Ring”. We discuss a formation scenario for the Per-Tau Shell, in which previous stellar and supernova (SN) feedback events formed a large expanding shell, where the swept-up ISM has condensed to form both the shell and the Perseus and Taurus molecular clouds within it. We present auxiliary observations of HI, Hα, 26Al, and X-rays that further support this scenario, and estimate Per-Tau Shell’s age to be ≈ 6 − 22 Myrs. The Per-Tau shell offers the first three- dimensional observational view of a phenomenon long-hypothesized theoretically, molecular cloud formation and star formation triggered by previous stellar and SN feedback. Subject headings: Interstellar medium (847) – Molecular clouds (1072) – Solar neighborhood (1509) – Stellar feedback (1602) – Superbubbles (1656) – Astronomy data visualization (1968)

1. INTRODUCTION et al. 2009; Lee & Chen 2009; Dawson et al. 2011, 2015; The conversion of diffuse gas into a dense-cold phase is the Mackey et al. 2017, see reviews by Elmegreen 2011; Daw- first step and a potential bottleneck for star formation. The son 2013). However, a major limitation with 2D extinction balance between heating and cooling processes in the ISM re- maps and PPV cubes is that they do not include information sults in a multiphase gas, consisting of the cold/warm neutral on the distance along the line-of-sight (LOS), and the data are media (CNM/WNM; Field et al. 1969; Wolfire et al. 2003; limited to the 2D “plane-of-the sky.” This results in blending Seifried et al. 2010; Saury et al. 2014; Bialy & Sternberg of structures along the LOS and a limited access to the real 2019). The CNM, being colder and denser, is the phase sus- 3D structure of gas in the ISM. This is especially a problem ceptible to molecule formation, gravitational collapse and star for structures that are only marginally denser than the sur- formation. Numerical simulations have shown that the WNM- rounding ISM, such as old SN remnants, which may be major to-CNM conversion takes place in converging flows where drivers of WNM-CNM conversion, molecule formation and WNM is shock compressed and then radiatively cools, form- star formation triggering (see Fig. 2 in Inutsuka et al. 2015). ing CNM (Hennebelle & Perault´ 1999; Koyama & Inutsuka Since the advent of Gaia (Brown et al. 2018), rapid progress 2000, 2002; Audit & Hennebelle 2005;V azquez-Semadeni´ has been made in the field of 3D dust mapping, which lever- et al. 2006; Heitsch et al. 2006; Inoue & Inutsuka 2008, 2009; ages stellar distance and extinction estimates made by model- Hennebelle & Inutsuka 2019). These flows naturally occur in ing photometry and astrometry to chart the 3D distribution of expanding superbubbles (SB) powered by SN and winds from dust (Lallement et al. 2019; Green et al. 2019; Chen et al. massive stars, introducing a positive feedback-loop for star 2019; Leike & Enßlin 2019). The basic idea is intuitive: formation: [star formation] → [expanding shell] → [gas com- starlight extinction measures the dust column density whereas pression] → [gravitational collapse] → [star formation] (e.g., a parallax measurement constrains the star’s distance. Synthe- Hartmann et al. 2002; Hosokawa & Inutsuka 2006; Ntormousi sizing column densities and distances towards a large num- et al. 2011; Palousˇ 2014). ber of stars provides information on the dust 3D density. Re- The formation of dense gas and stars in expanding shells cently, Leike et al.(2020, herefter L20) extended this frame- arXiv:2109.09763v1 [astro-ph.GA] 20 Sep 2021 has been observed in CO and HI spectral position-position- work, combining metric Gaussian variational inference with velocity (PPV) cubes and in 2D dust extinction maps (e.g. Su Gaussian Processes to construct a 3D dust map of the solar neighborhood’s ISM at 1 pc resolution.

? In this paper, we use the L20 data to study the 3D den- [email protected] sity structure of the ISM in the Taurus-Perseus star-forming 1 Center for Astrophysics | Harvard & Smithsonian, 60 Garden st., Cambridge, MA, 02138 region (see Fig.1). This region is amongst the most thor- 2 Radcliffe Institute for Advanced Study, Harvard University, 10 Garden oughly observed and studied star-forming regions in the Solar Street, Cambridge, MA 02138, USA Neighborhood. When large-scale mapping of molecular gas 3 University of Vienna, Department of Astrophysics, first became feasible, Ungerechts & Thaddeus(1987, here- Turkenschanzstrasse¨ 17, 1180 Wien, Austria 12 4 Department of Physics, University of Wisconsin–Whitewater, 800 after UT87) mapped the Perseus-Taurus region in CO, and West Main Street, Whitewater, WI 53190 included a section in their paper entitled “Are the Taurus and 5 Max Planck Institute for Astrophysics, Karl-Schwarzschildstraße IC 348 (a prominent part of Perseus) Clouds Connected?.” In 1,85748 Garching, Germany velocity space, the CO data show a smooth connecting bridge 6 Ludwig-Maximilians-Universitat,¨ Geschwister-SchollPlatz1, 80539 between Taurus and Perseus, suggesting that a physically con- Munich, Germany MNRAS 000, 000–000 (0000) Preprint 22 September 2021 Compiled using MNRAS LATEX style file v3.0 First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations Philip F. Hopkins1, Iryna S. Butsky2, Georgia V. Panopoulou1, Suoqing Ji1, Eliot Quataert3, Claude-André Faucher-Giguère4, Dušan Kereš5 1 TAPIR, Mailcode 350-17, California Institute of Technology, Pasadena, CA 91125, USA. E-mail:[email protected] 2 Astronomy Department, University of Washington, Seattle, WA 98195, USA 3 Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA 4 Department of Physics and Astronomy and CIERA, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA 5 Department of Physics, Center for Astrophysics and Space Science, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA

ABSTRACT We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV-TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live kinetic-MHD galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighborhood (LISM) observations (spectra, B/C, e+/e−, p¯/p, 10Be/9Be, ionization). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circum-galactic medium in realistic galaxies necessarily imposes a 10kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary∼ spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms domi- nated by 0.1 1kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from∼ local− structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov-Taylor stages of SNe remnants. Key words: cosmic rays — plasmas — methods: numerical — MHD — galaxies: evolution — ISM: structure

1 INTRODUCTION space and time (even at a given galacto-centric radius) as CRs prop- agate through the ISM. Likewise, the injection itself being propor- The propagation and dynamics of cosmic rays (CRs) in the interstel- tional to e.g. SNe rates is strongly clustered in both space and time lar medium (ISM) and circum/inter-galactic medium (CGM/IGM) and specifically related to certain ISM phases (see Evans 1999; is an unsolved problem of fundamental importance for space plasma Vázquez-Semadeni et al. 2003; Mac Low & Klessen 2004; Walch physics as well as star and galaxy formation and evolution (see re- et al. 2015; Fielding et al. 2018), and other key loss terms depend views in Zweibel 2013, 2017; Amato & Blasi 2018; Kachelrieß & on e.g. local ionized vs. neutral fractions, magnetic and radiation Semikoz 2019). For decades, the state-of-the-art modeling of Galac- energy densities – quantities that can vary by ten orders of mag- tic (Milky Way; MW) CR propagation has largely been dominated nitude within the MW (Wolfire et al. 1995; Evans 1999; Draine by idealized analytic models, where a population of CRs is propa- 2011). And these static models cannot, by construction, capture gated through a time-static MW model, with simple or freely-fit as- non-linear effects of CRs actually modifying the galaxy/ISM struc- sumptions about the “halo” or thick disk around the galaxy and no ture through which they propagate. This in turn means that most in- appreciable circum-galactic medium (CGM)1 with “escape” (as a ferred physical quantities such as CR diffusivities. residence times, leaky box or flat halo-diffusion type model) outside of some radius re-acceleration efficiencies, and “convective” speeds (let alone their (Blasi & Amato 2012a; Strong & Moskalenko 2001; Vladimirov dependence on CR energy or ISM properties) are potentially sub- et al. 2012; Gaggero et al. 2015; Guo et al. 2016; Jóhannesson et al. ject to order-of-magnitude systematic uncertainties. That is not to 2016; Cummings et al. 2016; Korsmeier & Cuoco 2016; Evoli et al. say these static-Galaxy models are simple, however: their complex- 2017). ity focuses on evolving an enormous range of CR energies from These calculations generally ignore phase structure or inho- MeV to PeV, including a huge number of different species, and mogeneity in the ISM/CGM, magnetic field structure (anisotropic . & incorporating state-of-the-art nuclear networks for detailed spalla- arXiv:2109.09762v1 [astro-ph.HE] 20 Sep 2021 CR transport), streaming, complicated inflow/outflow/fountain and tion, annihilation, and other reaction rates (recently, see Liu et al. turbulent motions within the galaxy, and time-variability of galactic 2018; Amato & Blasi 2018). structure and ISM phases (although see e.g. Blasi & Amato 2012b; Jóhannesson et al. 2016; Liu et al. 2018; Giacinti et al. 2018), even Meanwhile, simulations of galaxy structure, dynamics, evolu- though, for example, secondary production rates depend on the lo- tion, and formation have made tremendous progress incorporating cal gas density which varies by several orders of magnitude in both and reproducing detailed observations of the time-dependent, multi- phase complexity of the ISM and CGM (Hopkins et al. 2012a; Kim & Ostriker 2017; Grudic´ et al. 2019; Benincasa et al. 2020; Keating 1 The term “halo” is used differently in CR and galaxy literature. In most CR literature, the “halo” is generally taken to have a size ∼ 1 − 10 kpc, cor- et al. 2020; Gurvich et al. 2020), galaxy inflows/outflows/fountains responding to the “thick disk” or “disk-halo interface” region in galaxy for- (Narayanan et al. 2006; Hayward & Hopkins 2017; Muratov et al. mation/structure terminology. In the galaxy community, the gaseous “halo” 2017; Anglés-Alcázar et al. 2017; Hafen et al. 2019b,a; Hopkins usually refers to the circum-galactic medium (CGM), with scale-lengths et al. 2020a; Ji et al. 2020b), and turbulent motions (Hopkins ∼ 20 − 50kpc and extent ∼ 200 − 500kpc (Tumlinson et al. 2017). 2013a,b; Guszejnov et al. 2017b; Escala et al. 2018; Guszejnov et al.

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

Simulating the complexity of the dark matter sheet II: halo and subhalo mass functions for non-cold dark matter models

Jens Stücker,1? Raul E. Angulo1,2 Oliver Hahn,3,4 and Simon D.M. White5 1Donostia International Physics Centre (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastian, Spain. 2IKERBASQUE, Basque Foundation for Science, E-48013, Bilbao, Spain. 3Department of Astrophysics, University of Vienna, Türkenschanzstrasse 17, 1180 Vienna, Austria 4Department of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria 5Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany.

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We present “sheet+release” simulations that reliably follow the evolution of dark matter structure at and below the dark matter free-streaming scale, where instabilities in traditional N-body simulations create a large population of spurious artificial haloes. Our simulations sample a large range of power-spectrum cutoff functions, parameterized through the half-mode scale khm and a slope parameter β. This parameter space can represent many non-cold dark matter models, including thermal relic warm dark matter, sterile-neutrinos, fuzzy dark matter, and a significant fraction of ETHOS models. Combining these simulations with additional N-body simulations, we find the following results. (1) Even after eliminating spurious haloes, the halo mass function in the strongly suppressed regime (nX/nCDM < 5%) remains uncertain because it depends strongly on the definition of a halo. At these mass scales traditional halo finders primarily identify overdensities that are unbound, highly elongated, dominated by tidal fields, or far from virialized. (2) The regime where the suppression is smaller than a factor of 20 is quite robust to these uncertainties, however, and can be inferred reliably from suitable N-body simulations. (3) Parameterizing the suppression in the halo- and subhalo mass functions through the scales where the suppression reaches 20%, 50% and 80%, we provide simple formulae which enable predictions for many non-cold dark matter models. (4) The halo mass-concentration relations in our sheet+release simulations agree well with previous results based on N-body simulations. (5) In general, we confirm the validity of previous N-body studies of warm dark matter models, largely eliminating concerns about the effects of artificial haloes. Key words: methods: numerical – cosmology: large-scale structure of Universe – dark matter

1 INTRODUCTION ferent properties will in general cause different behaviour for the detailed small-scale physics. However, at first order, the main effect While the cosmological parameters of our Universe have been mea- of these properties onto the formation of structures is the suppres- sured to an astonishing accuracy, some key questions about our sion of the power spectrum on small scales (Murgia et al. 2017). Universe remain unanswered. One of those key questions is: What is

arXiv:2109.09760v1 [astro-ph.CO] 20 Sep 2021 The usually considered canonical case for a dark matter candidate the nature of Dark Matter? So far, dark matter seems observationally with such a small-scale suppression is thermal relic warm dark mat- consistent with a perfectly cold, collisionless fluid. However, most ter. However, for example, also fuzzy dark matter (e.g. Niemeyer particle dark matter candidates are expected to deviate from this 2020) and many ETHOS models (Cyr-Racine et al. 2016) exhibit fiducial scenario at some scale. Such deviations can for example a small-scale suppression, just with a different precise shape than be caused by the warmth1 of dark matter, its self-interaction or its canonical warm dark matter. This suppression will already go a long quantum nature. way for explaining a lot of the difference to cold dark matter scenar- Dark matter cannot only be investigated by possible interac- ios – for example in the halo mass function. Here, we will refer to tions with known species in a particle physics lab, but also by the any dark matter model with a small-scale cutoff as “non-cold” dark imprints that its properties leave on the formation of structures. Dif- matter, and we aim to understand the general impact of a small-scale suppression in the power spectrum.

? E-mail: [email protected] There exist several promising probes that can be used to con- 1 more accurately: the primordial velocity dispersion strain (or possibly detect) non-cold dark matter candidates obser-

© 2021 The Authors Draft version September 22, 2021 Preprint typeset using LATEX style emulateapj v. 12/16/11

ALMA REVEALS LARGE OVERDENSITY AND STRONG CLUSTERING OF GALAXIES IN QUASAR ENVIRONMENTS AT Z ∼ 4 Cristina Garc´ıa-Vergara1,Matus Rybak1,2,Jacqueline Hodge1,Joseph F. Hennawi3,Roberto Decarli4,Jorge Gonzalez´ -Lopez´ 5,6,Fabrizio Arrigoni-Battaia7,Manuel Aravena5, and Emanuele P. Farina7 Draft version September 22, 2021

ABSTRACT We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4–3) line emitting galax- ies in 17 quasar fields at z ∼ 4 aimed to perform the first systematic search of dusty galaxies in high−z quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N≥ 5.6 within a projected radius of R . 1.5 h−1 cMpc and a velocity range of ∆ v = ±1000 km s−1 around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity +11.9 of 17.6−7.6 in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross- +2.42 −1 correlation length of r0,QG = 8.37−2.04 h cMpc, assuming a fixed slope γ = 1.8. This contradicts the reported mild overdensities (x1.4) of Lyman alpha emitters (LAEs) in the same fields at scales of R . 7 h−1 cMpc which is well described by a cross-correlation length 3 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at z ∼ 4, finding −1 an auto-correlation length of r0,CO = 3.14 ± 1.71 h cMpc (with γ = 1.8). Our work together with the previous study of LAEs around quasars traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multi-wavelength studies, and highlighting important questions about galaxy properties in high−z dense environments. Keywords: Quasars (1319), Quasar-galaxy pairs (1316), High-redshift galaxies (734), CO line emission (262), Clustering (1908), Large-scale structure of the universe (902), Submillimeter astronomy (1647)

1. INTRODUCTION the early universe are quasars, and thus they represent key Theoretical studies and simulations of the formation and objects to test structure formation models at early epochs. evolution of large-scale structure in the universe have been Quasars are accreting supermassive black holes (SMBHs), largely developed in the last decades, but observational con- and are some of the most luminous sources in the universe. straints in this field remain one of the major challenges for They are observed to peak at z ∼ 2 − 3 (e.g. Richards et al. astrophysicists. In the standard current paradigm, the cos- 2006) and have been detected up to z = 7.5 (Banados˜ et al. mological model is parametrized by the Lambda cold dark 2018; Yang et al. 2020). The quasar auto-correlation has been matter (ΛCDM) model, and the structures grow hierarchi- measured as a function of redshift up to z ∼ 4, and it has cally through gravitational instability (e.g. Dodelson 2003; been found to rise steadily over the redshift range 0 . z . 2 Padmanabhan 2006; Schneider 2015), forming a filamentary (Croom et al. 2005; Myers et al. 2006; Porciani & Norberg structure of dark matter and galaxies. In this framework, the 2006; Shen et al. 2007; da Angelaˆ et al. 2008) with a cluster- most massive galaxies in the early universe are placed in the ing similar to that of local galaxies, and to strongly increase most massive dark matter halos (Springel et al. 2005; Angulo between z ∼ 2 and z ∼ 4 (Shen et al. 2007; White et al. 2012; et al. 2012), tracing the most massive structures. Since the Eftekharzadeh et al. 2015; Timlin et al. 2018). dark matter halo mass can be directly related to the cluster- Additionally, the quasar auto-correlation function has been −1 ing of a population (Cole & Kaiser 1989; Mo & White 1996; shown to get progressively steeper on small (R . 40 h pkpc) Sheth & Tormen 1999), the generic prediction is that massive scales at z < 3 (Hennawi et al. 2006), and the small-scale galaxies should have a large auto-correlation function. clustering has been shown to also increase from z ∼ 3 to arXiv:2109.09754v1 [astro-ph.GA] 20 Sep 2021 The population with the largest auto-correlation function in z ∼ 4 (Shen et al. 2010). The measurements of z ∼ 4 quasar clustering result in an auto-correlation length r0 = 24.3 ± 2.4 h−1cMpc (at a fixed slope γ = 2.0) implying that 1 Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA quasars are the most clustered population at this redshift and Leiden, the Netherlands × 12 −1 2 Faculty of Electrical Engineering, Mathematics and Computer Sci- that they inhabit very massive (Mhalo > 6 10 M h ) dark ence, TU Delft, the Netherlands matter halos (Shen et al. 2007), in agreement with theoretical 3 Department of Physics, University of California, Santa Barbara, CA predictions (Springel et al. 2005; Angulo et al. 2012). These 93106, USA results imply that z ∼ 4 quasars are tracers of massive struc- 4 INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Gobetti 93/3, 40129 Bologna, Italy tures, and thus they should be surrounded by a large overden- 5 Nucleo´ de Astronom´ıa de la Facultad de Ingenier´ıa y Ciencias, Uni- sity of galaxies. versidad Diego Portales, Av. Ejercito´ Libertador 441, Santiago, Chile Several efforts have been made to detect such overdensities 6 Las Campanas Observatory, Carnegie Institution of Washington, in quasars environments, but the results reveal a contradic- Casilla 601, La Serena, Chile tory and confusing picture. Most of the quasar environment 7 Max-Planck-Institut fur¨ Astrophysik, Karl-Schwarzschild-Str 1, 85748 Garching bei Munchen,¨ Germany studies at z & 4 have been performed at optical wavelengths, Draft version September 22, 2021 Typeset using LATEX twocolumn style in AASTeX631

SDSS J1059+4251, a highly magnified 풛 ∼ 2.8 star-forming galaxy: ESI observations of the rest-frame UV spectrum

Annalisa Citro,1 Dawn K. Erb,1 Max Pettini,2 Matthew W. Auger,2 George D. Becker,3 and Bethan L. James4

1The Leonard E. Parker Center for Gravitation, Cosmology and Astrophysics, Department of Physics, University of Wisconsin-Milwaukee, 3135 N Maryland Avenue, Milwaukee, WI 53211, USA 2Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK 3Department of Physics & Astronomy, University of California, Riverside, CA 92521, USA 4AURA for ESA, Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218

ABSTRACT Detailed analyses of high-redshift galaxies are challenging due to their faintness, but this difficulty can be overcome with gravitational lensing, in which the magnification of the flux enables high signal-to-noise ratio (S/N) spectroscopy. We present the rest-frame ultraviolet (UV) Keck Echellette Spectrograph and Imager (ESI) spectrum of the newly discovered 푧 = 2.79 lensed galaxy SDSS J1059+4251. With an observed magnitude F814W = 18.8 and a magnification factor 휇 = 31 ± 3, J1059+4251 is both highly magnified and intrinsically ∗ 10 luminous, about two magnitudes brighter than 푀UV at 푧 ∼ 2–3. With stellar mass 푀∗ = (3.22±0.20) ×10 M , −1 star formation rate SFR = 50 ± 7 M yr , and stellar metallicity 푍∗ ' 0.15 − 0.5 푍 , J1059+4251 is typical of bright star-forming galaxies at similar redshifts. Thanks to the high S/N and the spectral resolution of the ESI spectrum, we are able to separate the interstellar and stellar features and derive properties that would be inaccessible without the aid of the lensing. We find evidence of a gas outflow with speeds up to −1000 km s−1, and of an inflow that is probably due to accreting material seen along a favorable line of sight. We measure relative elemental abundances from the interstellar absorption lines and find that 훼-capture elements are overabundant compared to iron-peak elements, suggestive of rapid star formation. However, this trend may also be affected by dust depletion. Thanks to the high data quality, our results represent a reliable step forward in the characterization of typical galaxies at early cosmic epochs.

Keywords: Galaxy evolution (594), High-redshift galaxies (734), Strong gravitational lensing (1643)

1. INTRODUCTION 2016, Rigby et al. 2018a, 2018b), another effective way to In the last two decades or so, the number of galaxies discov- overcome these difficulties is to study high redshift gravita- ered at redshifts 푧 ' 1–10 has increased spectacularly thanks tionally lensed galaxies, where the magnification furnished by largely to the development of appropriate colour selection cri- the gravitational lensing provides high S/N spectra which can teria, pioneered by Steidel et al.(1996) (see also Adelberger be analyzed in greater detail. In particular, rest-frame ultravi- et al. 2004), and reviewed by Shapley(2011) (see also Madau olet wavelengths are especially interesting, since they encode & Dickinson 2014). However, galaxies at high redshifts are meaningful information about the star formation processes

arXiv:2109.09748v1 [astro-ph.GA] 20 Sep 2021 ∗ that took place in the early Universe. UV stellar continua faint (푚R = 24.4 at 푧 = 2 − 3; Steidel et al. 1999; Reddy et al. 2008) and the signal-to-noise ratio (S/N) needed to en- can be used to infer the galaxy luminosity, star-formation rate able detailed spectroscopic studies is hard to achieve, until and dust extinction (e.g. Kennicutt 1998, Calzetti et al. 2000, the next generation of 30+ m optical/IR telescopes comes on Salim et al. 2007, Wilkins et al. 2012), while individual spec- line. Besides long integration times (e.g. Erb et al. 2010) tral features (in absorption and in emission) provide valuable and stacking techniques (e.g. Zhu et al. 2015, Steidel et al. information on both the galaxy stellar populations and the interstellar medium (ISM) (e.g. Shapley et al. 2003, Berry et al. 2012, Talia et al. 2012, Steidel et al. 2016, Du et al. Corresponding author: Annalisa Citro 2018, Cullen et al. 2019), including the chemical composi- [email protected] tion of young OB stars and the ISM from which these stars have recently formed. Furthermore, the relative abundances Multifield Cosmology with Artificial Intelligence

Francisco Villaescusa-Navarroa,b,1, Daniel Anglés-Alcázarc,b, Shy Genelb,d, David N. Spergelb,a, Yin Lib, Benjamin Wandelte,b, Andrina Nicolaa, Leander Thielef, Sultan Hassanb,g, Jose Manuel Zorrilla Matillaa, Desika Narayananh,i, Romeel Davej,g,k, and Mark Vogelsbergerl

aDepartment of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ, 08544, USA; bCenter for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, New York, NY, 10010, USA; cDepartment of Physics, University of Connecticut, 196 Auditorium Road, U-3046, Storrs, CT, 06269, USA; dColumbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY, 10027, USA; eSorbonne Universite, CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis boulevard Arago, 75014 Paris, France; fDepartment of Physics, Princeton University, Princeton, NJ 08544, USA; gDepartment of Physics & Astronomy, University of the Western Cape, Cape Town 7535, South Africa; hDepartment of Astronomy, University of Florida, 211 Bryant Space Sciences Center, Gainesville, FL, USA; iUniversity of Florida Informatics Institute, 432 Newell Drive, CISE Bldg E251, Gainesville, FL, USA; jInstitute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK; kSouth African Astronomical Observatories, Observatory, Cape Town 7925, South Africa; lKavli Institute for Astrophysics and Space Research, Department of Physics, MIT, Cambridge, MA 02139, USA

This manuscript was compiled on September 22, 2021

Astrophysical processes such as feedback from supernovae and ac- In the case of Gaussian density fields, such as the tempera- tive galactic nuclei modify the properties and spatial distribution of ture anisotropies of the early Universe, the power spectrum dark matter, gas, and galaxies in a poorly understood way. This un- is the summary statistic that completely characterizes the certainty is one of the main theoretical obstacles to extract informa- statistical properties of those fields. This means that, for these tion from cosmological surveys. We use 2,000 state-of-the-art hy- cases, the optimal procedure to extract the maximum amount drodynamic simulations from the CAMELS project spanning a wide of information is known: compute the power spectrum from variety of cosmological and astrophysical models and generate hun- the data and perform a maximum likelihood estimation of the dreds of thousands of 2-dimensional maps for 13 different fields: parameters given the theory predictions. from dark matter to gas and stellar properties. We use these maps to This procedure is routinely used to perform parameter train convolutional neural networks to extract the maximum amount inference from observables such as the spatial distribution of of cosmological information while marginalizing over astrophysical galaxies on sufficiently large scales, since their distribution effects at the field level. Although our maps only cover a small area resembles a Gaussian density field. Unfortunately, even if of (25 h−1Mpc)2, and the different fields are contaminated by as- the early Universe was a Gaussian density field, non-linear trophysical effects in very different ways, our networks can infer the gravitational evolution would transform it into a non-Gaussian values of Ωm and σ8 with a few percent level precision for most of density field on small scales and at low redshift. Numerous the fields. We find that the marginalization performed by the network studies have shown that extracting information using the power retains a wealth of cosmological information compared to a model spectrum in this regime is suboptimal (e.g.1–19). trained on maps from gravity-only N-body simulations that are not For generic non-Gaussian fields, the summary statistic(s) contaminated by astrophysical effects. Finally, we train our networks that fully characterizes their statistical properties is unknown. on multifields – 2D maps that contain several fields as different col- In other words, the procedure to extract the maximum infor- ors or channels – and find that not only they can infer the value of all mation from these fields remains a mystery. This represents a parameters with higher accuracy than networks trained on individual major challenge that limits our ability to extract the maximum fields, but they can constrain the value of Ωm with higher accuracy information from cosmological surveys. A second challenge than the maps from the N-body simulations.

cosmology, astrophysics, deep learning, simulations Significance Statement osmology is becoming a precise branch of astronomy. A wealth of information from cosmological surveys remains CQuantities characterizing critical properties of our Uni- out of reach owing to (1) the lack of an optimal method to ex- verse, such as its age and geometry, are known with relatively tract that information and (2) the contamination introduced by high precision. On the other hand, some fundamental ques- uncertain astrophysical processes such as the impact of super- tions remain unanswered: What is the nature of dark energy? nova explosions and growing supermassive black holes. This What are the neutrino masses? Is General Relativity correct hinders our ability to learn about fundamental physics behind arXiv:2109.09747v1 [astro-ph.CO] 20 Sep 2021 on large scales? Data from upcoming cosmological surveys the laws and constituents of our Universe from cosmological such as Euclid∗, eRosita†, DESI‡, PFS§, Rubin observatory¶, data. We show for the first time that neural networks trained Roman observatory‖, SKA∗∗, and Simons Observatory†† are on thousands of galaxy formation simulations with varying cos- expected to shed light on these questions. Unfortunately, the mological and astrophysical models can help overcome these procedure needed to extract the maximum amount of cosmo- challenges by extracting the relevant cosmological informa- logical information from such surveys is currently unknown. tion directly from multiple raw data fields while simultaneously marginalizing over uncertainties in astrophysical effects. ∗https://www.euclid-ec.org † https://www.mpe.mpg.de/eROSITA Author contributions: FVN, DAA, SG, DNS, YL, and BW designed research; FVN performed re- ‡ https://www.desi.lbl.gov search; FVN, DAA, SG, YL, DN, and RM contributed software; FVN, DAA, SG, DNS, YL, BW, AN, §https://pfs.ipmu.jp LT, SH, JMZM, analyzed data; FVN, DAA, SG, YL, BW, AN, LT, SH, JMZM, DN, and RM wrote the paper ¶https://www.lsst.org ‖https://roman.gsfc.nasa.gov The authors declare no conflict of interest. ∗∗https://www.skatelescope.org 1 ††https://simonsobservatory.org E-mail: [email protected]

www.pnas.org/cgi/doi/10.1073/pnas.XXXXXXXXXX PNAS | September 22, 2021 | vol. XXX | no. XX | 1–11

arises because astrophysical processes, such as feedback from on each channel, outperforming the results from maps not supernovae and active galactic nuclei (AGN), are expected to contaminated by astrophysical effects. Finally, we discuss how affect the distribution and properties of matter and galaxies in our results can be used for forward models and the challenges a poorly understood manner; for a generic statistic, the effect associated with this method before they can be applied to of baryonic physics is largely unknown. data from cosmological surveys. One possibility to tackle these problems is to use artificial neural networks to extract information from the field while Results simultaneously marginalizing over uncertain baryonic effects. We now present the main results of this work. We first il- We have recently demonstrated using Gaussian density fields lustrate in detail our method with one particular field: gas that neural networks can achieve both goals (20). While temperature. We then describe the results obtained for all many works have used neural networks to extract cosmological other fields. Next, we present the results for the multifield information from gravity-only N-body simulations (21–29), maps and the N-body simulations. Unless otherwise stated, to our knowledge no such analysis has been performed using we only show results obtained from maps from the 1,000 Illus- state-of-the-art cosmological hydrodynamic simulations with trisTNG hydrodynamic simulations belonging to the CAMELS full galaxy formation physics and leveraging the information project. Results are very similar if we instead use maps from content of multiple baryonic fields simultaneously. the 1,000 SIMBA hydrodynamic simulations from CAMELS, In this work, we make use of data from the Cosmology and which are run with a completely different code and utilize a Astrophysics with MachinE Learning Simulations (CAMELS) very distinct subgrid physics implementation. We emphasize project (30) to perform likelihood-free inference of cosmo- that unless otherwise stated, we show results when training logical parameters from 2D projected distributions of a vari- networks on maps from the IllustrisTNG simulations and test- ety of fields evolved in the simulations. CAMELS contains ing the model on maps from the same set of simulations‡‡. In more than 4,000 cosmological N-body and state-of-the-art the discussion section we describe the plausible outcome of hydrodynamic simulations, where each simulation represents training the networks on maps from one simulation suite and a Universe with a different cosmological and astrophysical test them on maps from the other simulation suite. model. The hydrodynamic simulations have been run with two different codes and baryonic subgrid models: IllustrisTNG Gas temperature. In the upper panels of Fig.1 we show exam- (31, 32) and SIMBA (33) simulations, that include radiative ples of one of the 13 fields considered in this paper: gas tem- cooling and heating, star formation, stellar evolution, galac- perature. Each map represents a region of 25 × 25 (h−1Mpc)2 tic winds, the formation and growth of supermassive black and contains 256 × 256 pixels; the value of the cosmological holes, and feedback from AGN; this set of input physics can and astrophysical parameters is different for each map. In yield observationally-concordant galaxy populations in current terms of the astrophysical parameters, ASN1 and ASN2 control models (Somerville+Dave 2015). The value of the cosmologi- the efficiency of galactic winds, while AAGN1 and AAGN2 char- cal and astrophysical parameters in the simulations is varied acterize AGN feedback. Since each simulation was run with a across a wide range, so that posterior location and width from different value of the random seed used to generate the initial simulated data from the middle of the range are not affected conditions, each map also represents a different realization of by priors. the cosmic web. The regions where the gas is colder are shown From the hydrodynamic simulations, we generate tens of in black, and they tend to correspond to cosmic voids. In white thousands of 2D maps for 13 different fields (from dark matter we show the regions where the gas is hotter, corresponding to to gas and stellar properties), which are affected by astrophys- places where dark matter halos hosting galaxies are located. ical effects at very different levels. These maps represent a These regions are typically surrounded by hot bubbles whose small area of 25 × 25 (h−1Mpc)2 and have a spatial resolution morphology depends on cosmology but also on the efficiency of ' 100 h−1kpc. We then feed these maps to convolutional of the astrophysical processes. These maps clearly illustrate neural networks that are trained to perform likelihood-free the diversity of cosmological and astrophysical models present inference on the value of the cosmological parameters while in the CAMELS simulations. marginalizing over the astrophysical effects at the field level. Our goal in this paper is to infer cosmological and astro- We show that neural networks not only can marginalize physical parameters directly from these maps. To achieve over baryonic effects at the field level, but they also achieve this, we take 13,500 of these maps (training set) and train a percent level constraints on the value of the cosmological pa- neural network to predict the mean and standard deviation rameters for most of the fields. This illustrates the enormous of the marginal posterior for each of the six cosmological and amount of cosmological information that is embedded in small astrophysical parameters (see the Neural networks subsection scales in the fully non-linear regime. We compare our results below for details). Once the network is trained, we input to to those obtained from maps that are not contaminated by it 40 maps that it has never seen before (belonging to the astrophysical effects (the maps from the gravity-only N-body test set) and show the results in the bottom panels of Fig. simulations), and we show that the network does not lose much 1. We find that the network is able to infer the value of the information when performing the marginalization over the as- cosmological parameters with a surprisingly high accuracy: trophysical effects. In addition, we train neural networks to hδΩm/Ωmi ' 3.5% and hδσ8/σ8i ' 5.3%. We emphasize that perform likelihood-free inference on multifields, defined as 2D these constraints are obtained from the maps themselves (i.e. maps that contain multiple fields simultaneously as different without using intermediate summary statistics) and incorpo- colors or channels. We show for the first time that the multi- rate the marginalization over astrophysical effects as well as field configuration provides the tightest constraints on most the effects of cosmic variance. parameters despite the widely different astrophysical effects ‡‡The maps in the testing set have not been seen by the neural network though.

2 | www.pnas.org/cgi/doi/10.1073/pnas.XXXXXXXXXX Francisco Villaescusa-Navarro et al. Draft version September 22, 2021 Typeset using LATEX twocolumn style in AASTeX63

Optical to X-ray Signatures of Dense Circumstellar Interaction in Core-Collapse Supernovae

Ben Margalit ,1, ∗ Eliot Quataert ,2 and Anna Y. Q. Ho 3, 4

1Astronomy Department and Theoretical Astrophysics Center, University of California, Berkeley, Berkeley, CA 94720, USA 2Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA 3Department of Astronomy, University of California, Berkeley, 501 Campbell Hall, Berkeley, CA, 94720, USA 4Miller Institute for Basic Research in Science, 468 Donner Lab, Berkeley, CA 94720, USA

ABSTRACT Progenitors of core-collapse supernovae (SNe) can shed significant mass to circumstellar material (CSM) in the months–years preceding core-collapse. The ensuing SN explosion launches ejecta that may subsequently collide with this CSM, producing shocks that can power emission across the electro- magnetic spectrum. In this work we explore the thermal signatures of dense CSM interaction, when the CSM density profile is truncated at some outer radius. CSM with optical depth > c/v (where v is the shock velocity) will produce primarily ∼blackbody optical/UV emission whereas lower optical- depth CSM will power bremsstrahlung X-ray emission. Focusing on the latter, we derive light-curves and spectra of the resulting X-ray transients, that include a detailed treatment of Comptonization. Due to strong photoelectric absorption, the X-ray light-curve is dominated by the ‘post-interaction’ phase that occurs after the shock reaches the CSM truncation radius. We treat this regime here for the first time. Using these results, we present the phase-space of optical, UV, and X-ray transients as a function of CSM properties, and discuss detectability prospects. We find that ROSAT would not have been sensitive to CSM X-ray transients but that eROSITA is expected to detect many such events. Future wide-field UV missions such as ULTRASAT will dramatically enhance sensitivity to large optical-depth CSM configurations. Finally, we present a framework within which CSM properties may be directly inferred from observable features of X-ray transients. This can serve as an important tool for studying stellar mass loss using SN X-ray detections.

Keywords: Time domain astronomy (2109) — X-ray transient sources (1852) — High energy as- trophysics (739) — Ultraviolet transient sources (1854) — Core-collapse supernovae (304)

1. INTRODUCTION will ultimately shock the surrounding CSM. This CSM Towards the end of their lives, massive stars can shed interaction is observable in a myriad of ways. significant mass through winds or eruptions, thus pol- Type IIn SNe exhibit narrow emission lines indicative luting their immediate environment with dense circum- of ionization of the surrounding CSM, and their optical stellar material (CSM). The death of these stars results light-curves have long been modelled as being powered in a violent supernova (SN) explosion, whose observa- by CSM interaction (e.g. Smith 2017). Radio SNe are arXiv:2109.09746v1 [astro-ph.HE] 20 Sep 2021 tional features may be markedly shaped by such CSM. another well-studied class of events, whose bright radio First light produced by the SN (at shock breakout from luminosity are powered by non-thermal emission in SN- the stellar surface) can flash-ionize nearby CSM and pro- ejecta–CSM shocks (Chevalier 1982; Weiler et al. 1986; duce narrow emission lines (e.g. Yaron et al. 2017; Bruch Chevalier 1998; Weiler et al. 2002; Chevalier & Fransson et al. 2021). As material ejected from the SN expands, it 2017). More recently, wide-field optical surveys are discov- ering rare optical transients whose light-curves cannot Corresponding author: Ben Margalit be powered by radioactivity (Quimby et al. 2007; Drout [email protected] et al. 2014). CSM-interaction is an appealing alternative ∗ NASA Einstein Fellow energy source. One class of such optical transients are superluminous SNe (SLSNe; Gal-Yam 2012). These SNe A new type of exoplanet direct imaging search: the SCExAO/CHARIS survey of accelerating stars∗

Thayne Curriea,b,c, Timothy D. Brandtd, Masayuki Kuzuharae, Jeffrey Chilcotef, Edward Cashmang,c, R. Y. Liuh,c, Kellen Lawsoni, Taylor Tobinf, G. Mirek Brandtd, Olivier Guyonb,e,j,k, Julien Lozib, Vincent Deob, Sebastien Vievardb, Kyohoon Ahnb, and Nour Skafb,l,m aNASA-Ames Research Center, Moffett Field, California, USA bSubaru Telescope, 650 N. Aohoku Pl., Hilo, Hawai’i, USA cEureka Scientific, Oakland, California, USA dUniversity of California-Santa Barbara, Santa Barbara, CA, USA eAstrobiology Center of NINS, 2-21-1, Osawa, Mitaka, Tokyo, Japan fUniversity of Notre Dame, South Bend, IN, USA gUniversity of Hawai‘i-Hilo, Hilo, HI, USA hColumbia University, New York, NY, USA iUniversity of Oklahoma, Norman, Oklahoma, USA jSteward Observatory, University of Arizona, Tucson, AZ, USA kJames C. Wyant College of Optical Sciences, University of Arizona, Tucson, AZ, USA 1LESIA, Observatoire de Paris, Universit´ePSL, CNRS, Sorbonne Universit´e,Universit´ede Paris, 5 place Jules Janssen, 92195 Meudon, France mDepartment of Physics and Astronomy, University College London, London, United Kingdom

ABSTRACT We present first results from a new exoplanet direct imaging survey being carried out with the Subaru Coron- agraphic Extreme Adaptive Optics project (SCExAO) coupled to the CHARIS integral field spectrograph and assisted with Keck/NIRC2, targeting stars showing evidence for an astrometric acceleration from the Hippar- cos and Gaia satellites. Near-infrared spectra from CHARIS and thermal infrared photometry from NIRC2 constrain newly-discovered companion spectral types, temperatures, and gravities. Relative astrometry of com- panions from SCExAO/CHARIS and NIRC2 and absolute astrometry of the star from Hipparcos and Gaia together yield direct dynamical mass constraints. Even in its infancy, our survey has already yielded multiple discoveries, including at least one likely jovian planet. We describe how our nascent survey is yielding a far higher detection rate than blind surveys from GPI and SPHERE, mass precisions reached for known companions, and the path forward for imaging and characterizing planets at lower masses and smaller orbital separations than previously possible. Keywords: adaptive optics, extrasolar planets, infrared

arXiv:2109.09745v1 [astro-ph.EP] 20 Sep 2021 1. INTRODUCTION Over the past decade, facility and extreme adaptive optics (AO) systems on 8-10m class ground based-telescopes have provided the first direct detections of superjovian extrasolar planets.1–8 The vast majority of imaged exoplanets orbit at wide separations from their host star, typically in systems with ages between 1 and 100 Myr. Follow-up near to mid infrared (IR) photometry and spectroscopy have revealed key insights into the clouds, chemistries, and gravities of young gas giants.9–11 Dedicated exoplanet direct imaging surveys have shed light

∗ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. Further author information: (Send correspondence to T. Currie) T. Currie: E-mail: [email protected] DRAFTVERSION SEPTEMBER 22, 2021 Typeset using LATEX twocolumn style in AASTeX631

Bumpy Declining Light Curves are Common in Hydrogen-poor Superluminous Supernovae

GRIFFIN HOSSEINZADEH,1, 2 EDO BERGER,2 BRIAN D.METZGER,3, 4 SEBASTIAN GOMEZ,5, 2 MATT NICHOLL,6 AND PETER BLANCHARD7

1Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 2Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138-1516, USA 3Department of Physics and Columbia Astrophysics Laboratory, Columbia University, Pupin Hall, New York, NY 10027, USA 4Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, New York, NY 10010-5902, USA 5Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 6Birmingham Institute for Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 7Center for Interdisciplinary Exploration and Research in Astrophysics and Department of Physics and Astronomy, Northwestern University, 1800 Sherman Avenue, 8th Floor, Evanston, IL 60201, USA

(Received 2021 September 20) Submitted to The Astrophysical Journal

ABSTRACT Recent work has revealed that the light curves of hydrogen-poor (Type I) superluminous supernovae (SLSNe), thought to be powered by magnetar central engines, do not always follow the smooth decline predicted by a simple magnetar spin-down model. Here we present the first systematic study of the prevalence and properties of “bumps” in the post-peak light curves of 34 SLSNe. We find that the majority (44–76%) of events cannot be explained by a smooth magnetar model alone. We do not find any difference in the supernova properties between events with and without bumps. By fitting a simple Gaussian model to the light curve residuals, we characterize each bump with an amplitude, temperature, phase, and duration. We find that most bumps correspond with an increase in the photospheric temperature of the ejecta, although we do not see drastic changes in spectroscopic features during the bump. We also find a moderate correlation (ρ ≈ 0.5; p ≈ 0.01) between the phase of the bumps and the rise time, implying that such bumps tend to happen at a certain “evolutionary phase,” (3.7 ± 1.4)trise. Most bumps are consistent with having diffused from a central source of variable luminosity, although sources further out in the ejecta are not excluded. With this evidence, we explore whether the cause of these bumps is intrinsic to the supernova (e.g., a variable central engine) or extrinsic (e.g., circumstellar interaction). Both cases are plausible, requiring low-level variability in the magnetar input luminosity, small decreases in the ejecta opacity, or a thin circumstellar shell or disk.

Unified Astronomy Thesaurus concepts: Circumstellar matter (241), Circumstellar shells (242), Magne- tars (992), Supernovae (1668)

1. INTRODUCTION tral energy injection by a newly-born magnetar (Kasen & arXiv:2109.09743v1 [astro-ph.HE] 20 Sep 2021 In the decade since their discovery (Chomiuk et al. 2011; Bildsten 2010; Woosley 2010), central energy injection from Quimby et al. 2011; Gal-Yam 2012), a major focus in the fallback accretion onto a newly-born black hole (Dexter study of hydrogen-poor (Type I) superluminous supernovae & Kasen 2013), radioactive decay of an extremely large 56 (SLSNe) has been to identify the power source(s) behind mass of Ni (as expected for a pair-instability supernova; their extreme and long-lived luminosities. Four major mech- Kasen et al. 2011; Dessart et al. 2012), and/or kinetic en- anisms have been proposed, alone or in combination: cen- ergy released through collision with circumstellar material (CSM; as expected for a pulsational pair-instability super- nova; Woosley et al. 2007). Each of these proposals has Corresponding author: Griffin Hosseinzadeh strengths and weaknesses, but in general, the magnetar model griffi[email protected] has had the most success in explaining Type I SLSN light Draft version September 22, 2021 Typeset using LATEX twocolumn style in AASTeX631

Hard X-ray Irradiation Potentially Drives Negative AGN Feedback by Altering Molecular Gas Properties

Taiki Kawamuro,1, 2, ∗ Claudio Ricci,1, 3, 4 Takuma Izumi,2, 5, † Masatoshi Imanishi,2, 5 Shunsuke Baba,2, ‡ Dieu D. Nguyen,2 and Kyoko Onishi6

1Nu´cleo de Astronom´ıade la Facultad de Ingenier´ıa,Universidad Diego Portales, Av. Ej´eercito Libertador 441, Santiago, Chile 2National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan 3Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, People’s Republic of China 4George Mason University, Department of Physics & Astronomy, MS 3F3, 4400 University Drive, Fairfax, VA 22030, USA 5Department of Astronomy, School of Science, Graduate University for Advanced Studies (SOKENDAI), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 6Department of Earth and Space Sciences, Chalmers University of Technology, Onsala Space Observatory, 439 94 Onsala, Sweden

ABSTRACT To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and ALMA CO(J=2–1) data for 26 ultra-hard X-ray (> 10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray- irradiated gas via Fe-Kα emission, the CO(J=2–1) observations reveal that of cold molecular gas. At 00 00 high resolutions . 1 , we derive Fe-Kα and CO(J=2–1) maps for the nuclear 2 region, and for the external annular region of 200–400, where 200 is ∼ 100–600 pc for most of our AGNs. First, focusing on the external regions, we find the Fe-Kα emission for six AGNs above 2σ. Their large equivalent widths (& 1 keV) suggest a fluorescent process as their origin. Moreover, by comparing 6–7 keV/3– 6 keV ratio, as a proxy of Fe-Kα, and CO(J=2–1) images for three AGNs with the highest significant Fe-Kα detections, we find a possible spatial separation. These suggest the presence of X-ray-irradiated ISM and the change in the ISM properties. Next, examining the nuclear regions, we find that (1) The 20–50 keV luminosity increases with the CO(J=2–1) luminosity. (2) The ratio of CO(J=2–1)-to- HCN(J=1–0) luminosities increases with 20–50 keV luminosity, suggesting a decrease in the dense gas fraction with X-ray luminosity. (3) The Fe-Kα-to-X-ray continuum luminosity ratio decreases with the molecular gas mass. This may be explained by a negative AGN feedback scenario: the mass accretion rate increases with gas mass, and simultaneously, the AGN evaporates a portion of the gas, which possibly affects star formation.

Keywords: galaxies: active – X-rays: galaxies – submm/mm: galaxies

1. INTRODUCTION SMBHs are expected to have obtained a significant frac- Since the discovery of a potential link between the tion of their mass by gas accretion (e.g., Soltan 1982; spheroidal components of galaxies and their central su- Marconi et al. 2004; Ueda et al. 2014). Thus, in any permassive black holes (SMBHs; e.g., Magorrian et al. scenario related to the growth of galaxies and SMBHs 1998; Gebhardt et al. 2000; Marconi & Hunt 2003; (e.g., Di Matteo et al. 2005; Croton et al. 2006), it is arXiv:2109.09742v1 [astro-ph.GA] 20 Sep 2021 G¨ultekinet al. 2009; Kormendy & Ho 2013), the co- fundamental to understand detailed conversion process evolution of galaxies and SMBHs has been one of the of gas into key components (e.g., stars and SMBHs). most debated topics in astronomy. Generally, galaxies The properties of the ISM around an accreting SMBH, have grown by converting gas into stars, whereas the or an active galactic nucleus (AGN), likely differ from those in pure star-forming regions because the X-ray emission of an AGN is much stronger than that of star- Corresponding author: Taiki Kawamuro forming regions (e.g., Elvis et al. 1994; Mineo et al. [email protected] 2012; Hickox & Alexander 2018). Given that hard X-ray ( 1 keV) emission has high penetrating power, it can ∗ FONDECYT postdoctoral fellow & † deeply pierce through the ISM. Accordingly, the physi- NAOJ fellow cal and chemical properties of the ISM are expected to ‡ JSPS fellow (PD) be widely altered (e.g., Krolik & Kallman 1984; Mal- DRAFTVERSION WEDNESDAY 22ND SEPTEMBER, 2021 Preprint typeset using LATEX style emulateapj v. 01/23/15

LIKELIHOOD-FREE FORWARD MODELING FOR CLUSTER WEAK LENSING AND COSMOLOGY

SUT-IENG TAM1 ,KEIICHI UMETSU1 , AND ADAM AMARA2 Draft version Wednesday 22nd September, 2021

ABSTRACT Likelihood-free inference provides a rigorous approach to preform Bayesian analysis using forward simu- lations only. The main advantage of likelihood-free methods is its ability to account for complex physical processes and observational effects in forward simulations. Here we explore the potential of likelihood-free forward modeling for Bayesian cosmological inference using the redshift evolution of the cluster abundance combined with weak-lensing mass calibration. We use two complementary likelihood-free methods, namely Approximate Bayesian Computation (ABC) and Density-Estimation Likelihood-Free Inference (DELFI), to develop an analysis procedure for inference of the cosmological parameters (Ωm, σ8) and the mass scale of the survey sample. Adopting an eROSITA-like selection function and a 10% scatter in the observable–mass relation in a flat ΛCDM cosmology with Ωm = 0.286 and σ8 = 0.82, we create a synthetic catalog of observable-selected NFW clusters in a survey area of 50 deg2. The stacked tangential shear profile and the number counts in redshift bins are used as summary statistics for both methods. By performing a series of forward simulations, we obtain convergent solutions for the posterior distribution from both methods. We find that ABC recovers broader posteriors than DELFI, especially for the Ωm parameter. For a weak-lensing sur- −2 0.3 vey with a source density of ng = 20 arcmin , we obtain posterior constraints on S8 = σ8(Ωm/0.3) of 0.836 ± 0.032 and 0.810 ± 0.019 from ABC and DELFI, respectively. The analysis framework developed in this study will be particularly powerful for cosmological inference with ongoing cluster cosmology programs, such as the XMM-XXL survey and the eROSITA all-sky survey, in combination with wide-field weak-lensing surveys. Subject headings: cosmology: theory — dark matter — galaxies: clusters: general — gravitational lensing: weak

α 1. INTRODUCTION 2017), where S8 = σ8(Ωm/0.3) with α being a constant As the largest bound objects formed in the universe, galaxy that depends on the degree of parameter degeneracy (typi- clusters play a fundamental role in testing models of back- cally, α = 0.3–0.5; see Section5). In general, there are ground cosmology and structure formation. In the standard various challenging issues associated with cosmological tests picture of hierarchical structure formation, the abundance of using low-redshift probes, especially galaxy clusters, which cluster halos as a function of mass and redshift is sensitive to involve complex measurement processes and modeling in the the amplitude and growth rate of density fluctuations and the highly nonlinear regime of structure formation coupled with cosmic volume–redshift relation (e.g., Haiman et al. 2001). baryonic physics (Pratt et al. 2019). To obtain robust cosmo- Cluster number counts measured over a wide range in mass logical constraints from clusters in the present era of precision and redshift can thus provide powerful cosmological con- cosmology, one needs to conduct accurate statistical inference straints especially on the matter density parameter Ω and the accounting for various observational and instrumental effects m in modeling processes. amplitude of linear density fluctuations σ8 (defined in detail at the end of this section) (e.g., Mantz et al. 2015). In this con- Accurate calibration of cluster mass measurements is an- text, recent and ongoing cluster surveys covering a significant other critical ingredient of cluster cosmology (Pratt et al. fraction of the sky allow us to place stringent constraints on 2019). In cluster surveys, different observational techniques the cosmological parameters (e.g. de Haan et al. 2016; Schel- are employed to define an observable-selected cluster sam- lenberger & Reiprich 2017; Pacaud et al. 2018; Bocquet et al. ple using a low-scatter proxy that correlates with the underly- 2019; Costanzi et al. 2021; To et al. 2021; Chiu et al. 2021). ing cluster mass. With the assumption of hydrostatic equilib- arXiv:2109.09741v1 [astro-ph.CO] 20 Sep 2021 Cosmological parameters in the standard Λ cold dark mat- rium or virial theorem, these mass proxies can provide cluster ter (ΛCDM) model derived from low-redshift cosmological mass estimates, which however are expected to be biased by probes, such as galaxy clusters and cosmic shear, are often the presence of merging substructures, non-gravitational pro- in tension with those from observations of cosmic microwave cesses, or instrumentation effects (Nagai et al. 2007; Donahue background (CMB) anisotropies (Planck Collaboration et al. et al. 2014). Consequently, cosmological cluster studies of- ten require an external mass calibration of the survey sample 2020). Such apparent discrepancies in terms of Ωm and σ8 are often referred to as the “S tension” (e.g., Hildebrandt et al. using direct mass measurements (Planck Collaboration et al. 8 2016; Pacaud et al. 2018). Weak gravitational lensing offers a direct probe of the to- [email protected] [email protected] tal mass distribution around galaxy clusters projected along [email protected] the line of sight, irrespective of their dynamical state (e.g. von 1 Academia Sinica Institute of Astronomy and Astrophysics (ASIAA), der Linden et al. 2014; Umetsu et al. 2014; Hoekstra et al. No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan 2015; Okabe & Smith 2016; Medezinski et al. 2018; Diet- 2 Institute of Cosmology & Gravitation, University of Portsmouth, Den- rich et al. 2019; Herbonnet et al. 2020; Tam et al. 2020; Chiu nis Sciama Building, Burnaby Road, Portsmouth PO1 3FX, UK