Arxiv:2109.00912V1 [Hep-Th] 2 Sep 2021

Hydrogenated Amorphous Silicon Carbide: A Low-loss Deposited Dielectric for Microwave to Submillimeter Wave Superconducting Circuits

1, 1 1, 2 1 2 B. T. Buijtendorp, ∗ S. Vollebregt, K. Karatsu, D. J. Thoen, V. Murugesan, K. Kouwenhoven,1, 2 S. Hähnle,1, 2 J. J. A. Baselmans,1, 2 and A. Endo1 1Faculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Mekelweg 4, Delft 2628 CD, The Netherlands 2SRON Netherlands Institute for Space Research, Niels Bohrweg 4, Leiden 2333 CA, The Netherlands Low-loss deposited dielectrics will benefit superconducting devices such as integrated superconducting spectrometers, superconducting qubits and kinetic inductance parametric amplifiers. Com- pared with planar structures, multi-layer structures such as microstrips are more compact and elim- inate radiation loss at high frequencies. Multi-layer structures are most easily fabricated with deposited dielectrics, which typically exhibit higher dielectric loss than crystalline dielectrics. We measured the sub-kelvin and low-power microwave and mm-submm wave dielectric loss of hydrogenated amorphous silicon carbide (a-SiC:H), using a superconducting chip with NbTiN/a-SiC:H/NbTiN microstrip resonators. We deposited the a-SiC:H by plasma-enhanced chemical vapor deposition at a substrate temperature of 400°C. The a-SiC:H has a mm-submm loss tangent ranging from 4 4 0.80 0.01 10− to 1.43 0.04 10− in the range of 270–385 GHz. The microwave loss tangent ± × 5 ± × is 3.2 0.2 10− . These are the lowest low-power sub-kelvin loss tangents that have been reported for microstrip± × resonators at mm-submm and microwave frequencies. We observe that the loss tangent increases with frequency. The a-SiC:H films are free of blisters and have low stress: 20 MPa compressive at 200 nm thickness to 60 MPa tensile at 1000 nm thickness. −

Superconducting transmission lines are essential com- crystalline Si (c-Si) has led to the hypothesis that a more ponents of devices such as integrated superconducting constrained lattice (i.e. an increase in coordination num- spectrometers (ISSs) [1, 2], kinetic inductance paramet- ber) is correlated with a decrease in dielectric loss [13]. ric amplifiers (KIPAs) [3], and superconducting qubits With electron-beam evaporated amorphous silicon (a-Si) [4]. These superconducting devices operate at sub-kelvin it has been observed that a lower microwave loss can temperatures where amorphous solids exhibit excess di- be achieved by depositing at elevated substrate temper- electric loss compared to their crystalline counterparts atures (Tsub), and a correlation was found between the due to two-level systems (TLSs) [5, 6]. This is prob- microwave loss and the dangling bond density [18]. The lematic since the loss limits the efficiency of ISSs [7], electron-beam evaporated and sputtered a-Si films, which decreases the gain of KIPAs [8, 9] and is a significant are hydrogen-free, exhibit an order of magnitude larger source of decoherence in superconducting qubits [10]. microwave loss than the a-Si:H [13, 18], suggesting that The problem is often circumvented by avoiding deposited there is a relation between non-passivated dangling bonds dielectrics, for example by employing planar structures. and dielectric loss. However, a low-loss deposited dielectric will enable the Similar to a-Si:H, a-SiC:H is a hydrogenated and use of microstrips, thereby achieving benefits such as four-fold coordinated material. Deposition of low-stress miniaturization and the removal of radiation loss [11]. a-SiC:H films at high Tsub is possible with plasma- With hydrogenated amorphous silicon (a-Si:H) in mi- enhanced chemical vapor deposition (PECVD) [19]. The crostrip resonators, sub-kelvin and low-power loss tan- a-SiC:H has high wet etching selectivity with Si [19], en- 5 gents tan δ of 4–5 10− have been reported at mi- abling an a-SiC:H film to be used as a membrane, for × crowave frequencies [12–15]. Recently for a-Si:H a tan δ example in superconducting photon detectors [20]. The 4 of 2.1 0.1 10− has been reported at 350 GHz [14]. a-SiC:H can be patterned by reactive ion etching (RIE). ± × Although this is an improvement over SiO2 [13, 16] and Although the dielectric constants of a-SiC:H and of a- SiNx [7, 13, 17], a further reduction in tan δ is desirable. Si:H depend on the deposition parameters, a-SiC:H gen- For example, the efficiency of the mm-submm filters of erally has a lower microwave r 7 than the a-Si:H arXiv:2110.03500v1 [cond-mat.mtrl-sci] 7 Oct 2021 an ISS depends on R tan δ [7], where R = f/∆f is the microwave 10. The a-SiC:H≈ films do not exhibit · r resolving power. With the current mm-submm tan δ of blisters, which≈ are a common issue with a-Si:H [21, 22]. 4 2.1 10− of a-Si:H the filters are limited to R 500. A These properties suggest that a-SiC:H is a promising al- × ∼ higher R requires a reduction of tan δ. ternative to a-Si:H and SiNx for superconducting devices. Although the microscopic origin of the TLSs remains However, the sub-kelvin and low-power microwave and unknown [5, 6], a comparison of SiO2, SiNx, a-Si:H and mm-submm wave dielectric loss of a-SiC:H have not yet been reported. Here we present measurements of the sub-kelvin and low-power dielectric loss of a-SiC:H at microwave ∗ [email protected] (7 GHz) and mm-submm (270–385 GHz) frequen-

4 and the µWRs. We measured a room temperature εr of 10 WR1 6.9 at 15–66 THz using Fourier-transform infrared spec- i

WR2 Q troscopy and an εr of 7.5 at 250 THz (1200 nm) using

r o 5 WR3 5 t ellipsometry. We observe a bandgap of 1.8 eV, which we 10 c

10 a

WR4 f obtained by ﬁtting the Tauc-Lorentz dispersion model to y t

) 0.0 i l n B

a the ellipsometry data in the range of 0.8–2.7 eV. a d

t 2.5 u (

q |

l The observed frequency dependence of the dielectric

6 1 5.0 6 a 10 2 10 n S

7.5 r | loss of a-SiC:H, throughout microwave to mm-submm e t

10.0 n 6.520 6.525 I frequencies, is not explained by the STM which assumes

7 f (GHz) 7 a frequency independent TLS density [5]. A frequency- 10 10 2 1 0 1 2 3 4 5 6 7 8 9 10 10 10 10 10 10 10 10 10 10 10 10 dependent mm-submm tan δ has previously been ob- Average number of photons in the resonator N served for SiNx [17], and a larger mm-submm tan δ than microwave tan δ has been observed for both SiNx [13, 17] and a-Si:H [13, 14]. For a-Si:H [14] it has been suggested FIG. 5. Measured microwave tan δ = 1/ (pQi) versus av- that the mm-submm loss can be caused by phonon modes erage number of photons in the resonator, of µWR1–4. The at far-infrared wavelengths [29, 30], and the same could curves are ﬁts to Eq. 1, resulting in an average tan δ0 of be true for a-SiC:H. Alternatively, a generalized TLS 5 3.2 0.2 10− . The inset illustrates how tan δ is obtained at ± × model could account for a frequency dependent tan δ [31]. each power, by ﬁtting the S21 dips using the model in Ref. | | If the mm-submm loss of a-SiC:H is dominated by TLS [28]. loss, it is still an open question if the loss is mainly due to the bulk of the dielectric, or due to TLSs in the interface layers [32, 33]. N is equal to To conclude, PECVD a-SiC:H exhibits the lowest low- power sub-kelvin mm-submm tan δ that has been re- Pint N = , (4) ported to date for microstrip resonators. The microwave hf 2 tan δ is comparable to the best values that have been re-

where Pint is the resonator’s internal power. We fitted the ported for microstrip resonators [12, 14]. Interestingly, S dips of the four µWRs using the model described in we observe a frequency dependence of tan δ which is not | 21| Ref. [28], to obtain the resonators’ Qi. For illustration explained by the STM. In addition to the low tan δ, the we show one of the S dips in the inset of Fig. 5. a-SiC:H films are free of blisters and exhibit low stress. | 21| We obtained tan δ0 and β by fitting Eq. 1 to the mea- These properties make a-SiC:H a promising dielectric for sured tan δ versus N, as shown in Fig. 5. The µWRs microwave to submm wave superconducting circuits in 5 exhibit an average tan δ0 of 3.2 0.2 10− . This is many applications. comparable to what has been reported± × for a-Si:H [12– 14]. We measured β values in the range of 0.5–0.7. We confirm that the low-power loss is dominated by the ACKNOWLEDGEMENTS a-SiC:H by measuring a shunted NbTiN CPW quarter- wavelength resonator (CPWR) that was fabricated di- We thank the staff of the Else Kooi Laboratory and rectly on top of the c-Si substrate, which exhibits a tan δ0 the Kavli Nanolab Delft for their support. A.E. was sup- 6 of 4.1 0.1 10− . ported by the Netherlands Organization for Scientific Re- For± the a-SiC:H× we derive a sub-kelvin dielectric con- search NWO (Vidi grant no. 639.042.423). The contri- stant εr 7.0 for microwaves and mm-submm waves, bution of J.J.A. Baselmans was supported by the ERC based on≈ the resonance frequencies of the FP resonators CoG 648135 MOSAIC.

[1] A. Endo, K. Karatsu, Y. Tamura, T. Oshima, [2] K. S. Karkare, P. S. Barry, C. M. Bradford, S. Chapman, A. Taniguchi, T. Takekoshi, S. Asayama, T. J. L. C. S. Doyle, J. Glenn, S. Gordon, S. Hailey-Dunsheath, Bakx, S. Bosma, J. Bueno, K. W. Chin, Y. Fujii, R. M. J. Janssen, A. Kovács,H. G. LeDuc, P. Mauskopf, K. Fujita, R. Huiting, S. Ikarashi, T. Ishida, S. Ishii, R. McGeehan, J. Redford, E. Shirokoff, C. Tucker, R. Kawabe, T. M. Klapwijk, K. Kohno, A. Kouchi, J. Wheeler, and J. Zmuidzinas, Full-Array Noise Per- N. Llombart, J. Maekawa, V. Murugesan, S. Nakat- formance of Deployment-Grade SuperSpec mm-Wave subo, M. Naruse, K. Ohtawara, A. Pascual Laguna, On-Chip Spectrometers, J. Low Temp. Phys. 199, 849 J. Suzuki, K. Suzuki, D. J. Thoen, T. Tsukagoshi, (2020). T. Ueda, P. J. de Visser, P. P. van der Werf, S. J. C. [3] B. Ho Eom, P. K. Day, H. G. LeDuc, and J. Zmuidzinas, Yates, Y. Yoshimura, O. Yurduseven, and J. J. A. Basel- A wideband, low-noise superconducting amplifier with mans, First light demonstration of the integrated super- high dynamic range, Nat. Phys. 8, 623 (2012). conducting spectrometer, Nat. Astron. 3, 989 (2019). Dynamical Friction From Ultralight Dark Matter

Yourong Wang1, ∗ and Richard Easther1, † 1Department of Physics, The University of Auckland New Zealand 1010 (Dated: October 8, 2021) We simulate the gravitational dynamics of a massive object interacting with Ultralight / Fuzzy Dark Matter (ULDM/FDM), non-relativistic quantum matter described by the Schrödinger-Poisson equation. We first consider a point mass moving in a uniform background, and then a supermassive black hole (SMBH) moving within a ULDM soliton. After replicating simple dynamical friction scenarios to verify our numerical strategies, we demonstrate that the wake induced by a moving mass in a uniform medium may undergo gravitational collapse that dramatically increases the drag force, albeit in a scenario unlikely to be encountered astrophysically. We broadly confirm simple estimates of dynamical friction timescales for a black hole at the center of a halo but see that a large moving point mass excites coherent “breathing modes” in a ULDM soliton. These can lead to “stone skipping” trajectories for point masses which do not sink uniformly toward the center of the soliton, as well as stochastic motion near the center itself. These effects will add complexity to SMBH-ULDM interactions and to SMBH mergers in a ULDM universe.

I. INTRODUCTION test of the code we begin by recovering known analytic solutions (which are analogous to a much older We analyse a point-like massive particle inter- problem in electron propagation [25]) in the limit acting with self-gravitating quantum matter. The where the self-gravity of the quantum matter is ig- overall investigation is motivated by the dynamics nored and the point mass moves with a constant of super-massive black holes (SMBH) moving in- velocity. We study the full dynamics by allowing side Ultralight Dark Matter (ULDM) halos. ULDM, the ULDM to respond to its own gravitational po- also known as Fuzzy Dark Matter, is based on non- tential. A point mass moving in an otherwise undis- interacting particles with de Broglie wavelengths turbed ULDM background leaves a elongated over- long enough to influence galactic dynamics on sub- density in its wake, which eventually undergoes grav- kiloparsec scales [1–17]. Common realisations of itational collapse. The deep potential of the result- this scenario are built on axions with masses in the ing overdensity then brings the moving mass to a range 10−20 ∼ 10−23 eV. ULDM is non-relativistic rapid standstill. quantum matter interacting with its own Newtonian In the case of a black hole interacting with a gravitational potential and is thus governed by the ULDM halo the central soliton is already a collapsed nonlinear Schrödinger-Poisson equation. object. We consider SMBH in initially circular or- On large scales, ULDM resembles cold dark mat- bits and broadly confirm simple estimates of the dy- ter (CDM) but its quantum properties become ap- namical timescales on which they sink to the center parent on smaller scales, modifying the expectations of the soliton. However, a moving black hole ex- for intra-galactic dynamics relative to CDM [18]. cites oscillations in the overall soliton independently Given that conventional CDM faces a number of of the dynamical friction and the resulting dynam- challenges when confronted with the small-scale ics can be complicated and stochastic. In particular, arXiv:2110.03428v1 [gr-qc] 6 Oct 2021 properties of galaxies, these differences are the pri- we see possible evidence that the black hole will be mary motivation for ULDM models, and under- “reheated” as it interacts with a dynamical soliton, standing the detailed dynamics of ULDM will be key increasing the likelihood of core-stalling in SMBH to testing these scenarios. Moreover, in addition to mergers in a ULDM dominated universe. This work ULDM, self-gravitating quantum matter may arise is complementary to that of Lancaster et al. [26] in the very early universe [19, 20], hypothetical bo- which treats the point source in the perturbative son stars [21–23] and QCD axion miniclusters [24], limit and considers extended sources moving in more so the underlying dynamical system is relevant to realistic backgrounds – the distinction here is the fo- wide range of astrophysical systems. cus on soliton-SMBH dynamics. Conversely, Bar-Or This work rests on numerical solutions of the cou- et al. [27] discusses the “heating” of black hole bina- pled Schrödinger-Poisson equation. However, as a ries induced by the granular nature of ULDM halos, which is distinct from the phenomena examined in our simulations. ∗ [email protected] Our numerical simulations are based on PyUl- † [email protected] traLight [28], a pseudo-spectral Schrödinger- Probing light dark matter particles with astrophysical experiments

Tanmay Kumar Poddar푎,푏,∗ 푎Theoretical Physics Division, Physical Research Laboratory, Ahmedabad-380009, India 푏Discipline of Physics, Indian Institute of Technology, Gandhinagar-382355, India E-mail: [email protected]

The evidence of gravitational wave was ﬁrst indirectly conﬁrmed by the orbital period loss of Hulse-Taylor binary system which agrees well with the Einstein’s general relativistic prediction. The perihelion precession of planets, gravitational light bending and Shapiro time delay are other tests of Einstein’s general theory of relativity. However there are small uncertainties in the measurements of those observations from the general relativistic prediction. To account those uncertainties, we propose radiation of ultralight axions and vector gauge boson particles in the context of 푈(1)0 extended beyond standard model scenario. We obtain constraints on ultralight axion parameters (axion mass and decay constant) from the observational uncertainties of orbital period loss of compact binary systems, gravitational light bending, Shapiro time delay

and birefringence phenomena. We also obtain the bounds on ultralight 푈(1)퐿휇−퐿휏 gauge bosons from the orbital period loss of compact binary systems. The uncertainties in the perihelion

precession of planets also put bounds on the 푈(1)퐿푒−퐿휇,휏 light gauge bosons. These light particles can be promising candidates of fuzzy dark matter which can be probed from the above precision measurements. arXiv:2110.03365v1 [hep-ph] 7 Oct 2021

*** The European Physical Society Conference on High Energy Physics (EPS-HEP2021), *** *** 26-30 July 2021 *** *** Online conference, jointly organized by Universität Hamburg and the research center DESY ***

∗Speaker

Probing light dark matter particles with astrophysical experiments Tanmay Kumar Poddar

1. Introduction

The Standard Model (SM) of particle physics is a very successful theory in a sense that it can describe the interaction between all the fundamental particles in nature to a very good accuracy consistent with the experiments. However, the galactic rotation curve and the Bullet cluster experiment by Chandra X ray observatory conﬁrm the presence of non luminous /dark matter (DM) which cannot be explained by SM. The DM constitutes 25% of the total energy budget of the universe and its constituent is still unknown. To explain such beyond standard model (BSM) phenomenon, one can simply add particle or extend the gauge group with the SM. The weakly interacting massive particle (WIMP) model is one of the popular DM models and it can be a favourable DM candidate. Though, direct detection experiment put strong constraints on WIMP DM. Also it cannot explain the small scale structure problem of the universe. To resolve these problems, physicists consider alternative DM models and one of them is the fuzzy dark matter (FDM) model. In the following, we consider axions or axion like particles (ALPs), light gauge bosons as FDM candidates and they can form Bose-Einstein condensate. The de Broglie wavelength of the FDM particles is of the size of a dwarf galaxy and it behaves as a wave dark matter. In the following, we have considered some astrophysical experiments to probe such ultralight particles.

2. Constraints on ultralight axions from compact bianry systems

For Hulse-Taylor (HT) binary system, the general relativistic (GR) value of the orbital pe- ¤ −12 −1 ¤ riod loss is 푃푏퐺푅 = (−2.40263 ± 0.00005) × 10 푠푠 and its observed value is 푃푏표푏푠푒푟 푣푒푑 = −2.423(1) × 10−12푠푠−1. The observed value matches in good agreement with the GR prediction which is the indirect evidence of GW. However, there is less than 1% uncertainty in the observation from the GR prediction. To resolve this uncertainty, we assume that ultralight ALPs can radiate from the binary system and contribute to the orbital period loss. The orbital frequency of the HT binary system is Ω ∼ 10−19eV. Therefore, if particles with this small mass (which is in the ballpark of FDM particles) radiate from the binary systems, then it can contribute to the orbital period loss. If a compact object is immersed in a low mass axionic potential and it interacts with the nucleons of the compact star then it can develop a long range axionic hair outside of that object. In Table 1, we consider four compact binary systems and put bounds on the axion decay constants from observation. In Table1, 훼 denotes the ratio of axion mediated ﬁfth force to the gravitational force.

Compact binary system 푓푎 (GeV) 훼 PSR J0348+0432 . 1.66 × 1011 . 5.73 × 10−10 PSR J0737-3039 . 9.69 × 1016 . 9.21 × 10−3 PSR J1738+0333 . 2.03 × 1011 . 8.59 × 10−10 PSR B1913+16 . 2.07 × 1017 . 3.4 × 10−2

Table 1: Summary of the upper bounds on the axion decay constant 푓푎 from compact binaries. For all the −19 binaries we assume axion mass 푚푎 << 10 eV which is constrained from the orbital frequencies of the binaries.

11 −19 From Table1, we obtain the stronger bound on 푓푎 is . 10 GeV for ALPs of mass 푚푎 10 eV which implies if ALPs are FDM, they do not couple with quarks[1].

2 Prepared for submission to JHEP

Dark gauge boson production from neutron stars via nucleon-nucleon bremsstrahlung

Chang Sub Shina,b , Seokhoon Yunc,d aCenter for Theoretical Physics of the Universe, Institute for Basic Science (IBS), Daejeon, 34126, Korea bDepartment of Physics and Institute of Quantum Systems (IQS), Chungnam National University, Daejeon 34134, Republic of Korea cDipartimento di Fisica e Astronomia, Università degli Studi di Padova, Via Marzolo 8, 35131 Padova, Italy dIstituto Nazionale di Fisica Nucleare (INFN), Sezione di Padova, Via Marzolo 8, 35131 Padova, Italy E-mail: [email protected], [email protected]

Abstract: We discuss the dark gauge boson emission from neutron stars via nucleon- nucleon bremsstrahlung. Through the rigorous treatment of the effective field theory pre- scription and the thermal effect, we derive the relevant couplings of dark gauge bosons to

hadrons in medium. As a speciﬁc example, the U(1)B−L gauge boson scenario is chosen to investigate dark gauge boson emissivities during supernovae and cooling of young neutron stars. From the stellar cooling argument, we obtain the constraints on the B − L gauge coupling for given gauge boson masses in two observations: the duration of the supernova neutrino signal of SN1987A, and the inferred x-ray luminosity of the compact object in the remnant of SN1987A (NS1987A). In particular, the constraint from SN1987A is revisited, which is one order of magnitude enhanced compared to the previous derivation. arXiv:2110.03362v1 [hep-ph] 7 Oct 2021 KEK-TH-2349

SU(N)-natural inﬂation

Tomohiro Fujita,1, 2 Hiromasa Nakatsuka,3 Kyohei Mukaida,4, 5 and Kai Murai3, 6 1Waseda Institute for Advanced Study, Waseda University, Shinjuku, Tokyo 169-8050, Japan 2Research Center for the Early Universe, The University of Tokyo, Bunkyo, Tokyo 113-0033, Japan 3ICRR, University of Tokyo, Kashiwa, 277-8582, Japan 4Theory Center, IPNS, KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 5Graduate University for Advanced Studies (Sokendai), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 6Kavli IPMU (WPI), UTIAS, University of Tokyo, Kashiwa, 277-8583, Japan

We study SU(N) gauge fields that couple to the inflaton through the Chern-Simons term. We provide a general procedure to construct homogeneous, isotropic, and attractor solutions of the gauge fields during inflation. The gauge fields develop various VEVs corresponding to different spontaneous symmetry breaking patterns of SU(N) where embedded SU(2) subgroups are broken with the spatial rotation SO(3) symmetry. As specific examples, we construct the stable solutions for N 3 and 4. We also numerically solve the gauge = field dynamics and confirm that our analytic solutions are complete and attractor. It is straightforward to apply our procedure to the other simple Lie groups.

Introduction.— Cosmological inflation describes the show in this letter, most of homogeneous and isotropic at- early stage of the Universe successfully by solving the tractor solutions in the SU(N)-natural inflation have been horizon/flatness problems and providing the seeds of unexplored. anisotropies observed in the cosmic microwave back- In this letter, we provide a general procedure for find- ground (CMB). After inflation, the hot Universe should be ing homogeneous and isotropic attractor solutions in the generated, which requires a coupling between inflaton and SU(N)-natural inflation under three reasonable assump- other sectors. Nevertheless, the flatness of the inflaton po- tions. We find new solutions whose amplitudes of the gauge tential should be protected from radiative corrections to fields are larger than that in the chromo-natural inflation. have a sufficient duration of inflation. Axion-like particle is In addition, the effective potential implies that the solu- an attractive candidate for the inflaton [1–3] as its (approx- tion with the biggest amplitude is the most stable solution imate) shift symmetry controls the flatness. In addition, its for SU(N) with N 2. As we will see, each solution corre- > derivative couplings can induce rich phenomena, including sponds to a non-equivalent embedding of SU(2) subgroup the amplification of the gauge fields [4,5]. in SU(N), implying different spontaneous symmetry break- In particular, the model where the inflaton couples to ing patterns of SU(N) under a chemical potential of the SU(2) gauge fields through the Chern-Simons term, dubbed Chern-Simons number from a non-vanishing inflaton ve- chromo-natural inflation [6], has attracted much attention. locity. We also perform numerical simulations of the gauge In the presence of inflaton velocity, the gauge fields have field dynamics and confirm that the analytically found solutions are attractor and complete for N 3,4. a homogeneous and isotropic attractor solution [7–10] suit- = able to explain the current homogeneous and isotropic Uni- Model.— We consider a pseudo-scalar inflaton φ cou- a verse. Also, the backreaction from the gauge fields en- pled to SU(N) gauge fields Aµ through the Chern–Simons ables the slow-roll inflation with a sub-Planckian axion de- coupling: cay constant. Such background gauge fields induce linear 1 a aµν 1 µ φ a aµν couplings between the metric and gauge field tensor pertur- L F F ∂µφ∂ φ V (φ) F F˜ . (1) = −4 µν + 2 − + 4f µν bations, and the chiral gravitational waves are significantly a produced [11, 12]. Although the original chromo-natural in- The field strength of the SU(N) gauge fields Fµν and its dual arXiv:2110.03228v1 [hep-ph] 7 Oct 2021 flation is excluded from the CMB observations due to the Fãµν are defined by overproduction of gravitational waves [13,14], this conflict a a a abc b c can be evaded by introducing additional fields [15–20]. F ∂µ A ∂ν A g f A A , (2) µν ≡ ν − µ − µ ν 1 The SU(2) gauge group is the simplest example for the ãµν µνρσ a F ² Fρσ, (3) gauge fields to have a homogeneous and isotropic attrac- ≡ 2p g˜ tor solution. In general, we can extend the chromo-natural − abc inflation model to an SU(N) gauge group with N 2. We re- where f is the structure constant of the SU(N) algebra, > and the superscripts a,b,c run from 1 to N 2 1. g is the fer to this extension as “SU(N)-natural inflation”. The pre- − vious discussion on the SU(N) case was often limited to the gauge coupling constant, while g˜ is the determinant of the simplest choice of SU(2) subgroup, and it was considered spacetime metric. We assume that the background metric that the SU(N) case simply leads to the chromo-natural in- during inflation is well approximated by the de Sitter Uni- verse, ds2 dt 2 a(t)2dx2, with a exponentially growing flation (e.g., see [6]). Alternatively, some studies consider = − scale factor a(t) eHt and a constant Hubble parameter, N/2 copies of SU(2) subgroups in SU(N) and discuss the ∝ H const. production of gravitational waves [21–23]. However, as we = The Life and Science of Thanu Padmanabhan (1957-2021) remembered by The Paddy Gharana, with contributions from

Jasjeet Singh Bagla, Krishnakanta Bhattacharya, Sumanta Chakraborty, Sunu Engineer, Valerio Faraoni, Sanved Kolekar, Dawood Kothawala, Kinjalk Lochan, Sujoy Modak, V. Parameswaran Nair, Aseem Paranjape, Krishnamohan Parattu, Sarada G. Rajeev, Bibhas Ranjan Majhi, Tirthankar Roy Choudhury, Mohammad Sami, Sudipta Sarkar, Sandipan Sengupta, T. R. Se- shadri, S. Shankaranarayanan, Suprit Singh, Tejinder P. Singh, L. Sriramkumar and Urjit Yajnik

https://www.iucaa.in/~paddy/

ABSTRACT Thanu Padmanabhan was a renowned Indian theoretical physicist known for his research in general relativity, cosmology, and quantum gravity. In an extraordinary career spanning forty-two years, he published more than three hundred research articles, wrote ten highly successful technical and popular books, and mentored nearly thirty graduate students and post-doctoral fellows. He is best known for his deep work investigating gravitation as an emergent thermodynamic phenomenon. He was an outstanding teacher, and an indefatigable populariser of science, who travelled very widely arXiv:2110.03208v1 [physics.hist-ph] 7 Oct 2021 to motivate and inspire young students. Paddy, as he was aﬀectionately known, was also a close friend to his students and collaborators, treating them as part of his extended academic family. On September 17, 2021 Paddy passed away very unexpectedly, at the age of sixty-four and at the height of his research career, while serving as a Distinguished Professor at the Inter-University Centre for Astronomy and Astrophysics, Pune. His untimely demise has come as a shock to his family and friends and colleagues. In this article, several of them have come together to pay their tributes and share their fond memories of Paddy.

1 Angular emission patterns of remnant black holes

Xiang Li,1 Ling Sun,2, 3 Rico Ka Lok Lo,2 Ethan Payne,2, 3, 4, 5 and Yanbei Chen1 1Burke Institute for Theoretical Physics, California Institute of Technology, Pasadena, CA 91125, USA 2LIGO Laboratory, California Institute of Technology, Pasadena, CA 91125, USA 3OzGrav-ANU, Centre for Gravitational Astrophysics, College of Science, The Australian National University, ACT 2601, Australia 4School of Physics and Astronomy, Monash University, Clayton, VIC 3800, Australia 5OzGrav: The ARC Centre of Excellence for Gravitational Wave Discovery, Clayton, VIC 3800, Australia (Dated: October 8, 2021) The gravitational radiation from the ringdown of a binary black hole merger is described by the solution of the Teukolsky equation, which predicts both the temporal dependence and the angular distribution of the emission. Many studies have explored the temporal feature of the ringdown wave through black hole spectroscopy. In this work, we further study the spatial distribution, by introducing a global fitting procedure over both temporal and spatial dependences, to propose a more complete test of General Relativity. We show that spin-weighted spheroidal harmonics are the better representation of the ringdown angular emission patterns compared to spin- weighted spherical harmonics. The differences are distinguishable in numerical relativity waveforms. We also study the correlation between progenitor binary properties and the excitation of quasinormal modes, including higher-order angular modes, overtones, prograde and retrograde modes. Specifically, we show that the excitation of retrograde modes is dominant when the remnant spin is anti-aligned with the binary orbital angular momentum. This study seeks to provide an analytical strategy and inspire the future development of ringdown test using real gravitational wave events.

I. INTRODUCTION tivity [32–35], exploring such a problem defines the theme of this paper. The gravitational waves emitted at the final stage of a bi- In gravitational wave observations, the prospective searches nary black hole (BBH) merger – the ringdown stage, consist for ringdown waveforms [36, 37] would enable the test of the of a series of quasinormal modes (QNMs) [1–5]. QNMs are spatial-temporal features. With the rapidly increasing number solutions to the homogeneous Teukolsky equation [6–9], i.e. of binary coalescences observed [38–41] by ground-based de- the linearized Einstein’s equations in the background of a Kerr tectors like Advanced LIGO [42] and Advanced Virgo [43], black hole [10–14]. The foundation for doing so follows mod- events with detectable higher-order modes [44–46] are ob- els that describe stellar collapses [10–13] – the strong-field re- served, e.g., GW190412 [47] and GW190814 [48]. The de- gion “falls down” toward the future horizon of the final black tectability of higher-order modes not only impacts the param- hole, revealing a spacetime region in which perturbations sat- eter estimation [18, 19, 35, 49–51], but can be used to study isfy the homogeneous Teukolsky equation with ingoing con- angular emission as well. Currently, in the ringdown stage, a dition near the horizon, and outgoing condition near infinity. high signal-to-noise ratio (SNR) is difficult to achieve due to The homogeneous Teukolsky equation not only predicts the the lack of post-merger cycles and the degraded detector sensi- temporal dependence of the ringdown waves, in terms of their tivity at high frequency range. However, the sensitivity of the complex spectra, but also their spatial distributions, in terms proposed next-generation detectors, including Einstein Tele- of angular emission patterns. There have been many studies scope [52, 53], Cosmic Explorer [54, 55], and NEMO [56], on black hole ringdown spectroscopy involving multiple an- will be significantly improved [57], especially at the high fre- gular frequencies by modeling it as the superposition of expo- quencies, opening more possibilities in the BBH post-merger nentially damped sinusoids [15–19], or using other methods studies [24]. Although a single event could provide limited arXiv:2110.03116v1 [gr-qc] 7 Oct 2021 of frequency extraction [20, 21]. Recently, it has been shown information about angular dependence, combining multiple that the inclusion of overtones [4, 22] can improve the fitting events and extracting angular-dependent features will become of numerical relativity (NR) waveforms and lead to better es- possible with the expected large number of events in the fu- timation of ringdown model parameters [23, 24], because of ture. That calls for strategic studies of temporal-spatial emis- the better characterization of the post-merger signal from an sion patterns before more events with high ringdown SNR are earlier time. observed. Many phenomenological fitting studies based on NR wave- In this paper, we study ringdown gravitational waves and forms have been done [25–27], while most previous works show that the spin-weighted spheroidal harmonics are es- only focused on the ringdown temporal properties. Our study sential in the faithful representation of the temporal-spatial further includes spatial dependence on different models, ex- ringdown emission pattern [58], further confirming that the plicitly. Specifically, when the final spin is not aligned with Teukolsky equation describes the ringdown dynamics. We fit the initial orbital angular momentum, the spatial properties for the NR simulated strain data of merging binaries with different retrograde excitations [28–31] have not been carefully stud- parameters provided in the Simulating eXtreme Spacetimes ied. As the temporal-spatial consistency check of ringdown (SXS) Collaboration catalog [59], without adding simulated emission can provide a more complete test of General Rela- noise. The exclusion of nonlinear gravitational wave mem- arXiv:2110.03005v1 [gr-qc] 6 Oct 2021 nti es ehv osac o e olna em.Frti pur this For fields terms. magnetic s nonlinear is and new electrodynamics electric Faraday-Maxwell for the vacuum search the to that have [7] we Jackson sense this in ierin linear h efc iegneadtu a eeiiae.We n el w deals one invariant When true the eliminated. of be function can thus nonlinear and the divergence perfect the a fnniermdfiaino h xo lcrdnmc eassume co been form we has the Lagrangian electrodynamics has such axion of Lagrangian example the The of symmetries. mentioned modification nonlinear of way term easm htfrsmall for that assume We sa nee)o h suoclrfield pseudoscalar the of integer) an is ..[–]frdtiso hssoy.O hswyw edt aeaLag a have to need we way this On story). tensor this Maxwell of the details on for [2–6] e.g. ftesmere ftotps h rttp fsmer sconn symm is discrete symmetry the of to type respect first with The invariant is types. elect electrodynamics two axion of axion cosmic symmetries Nonlinear the the of of models theoretical the Extending invariant φ ∗ ˜ ‡ † [email protected] [email protected] [email protected] = h ocp fSmer a lydadcniust lyoutstan play to continues and played has Symmetry of concept The φ φF 2 + mn ∗ φ I πk F n oeigtebsctr ftelna xo lcrdnmc 4 5 [4, electrodynamics axion linear the of term basic the covering and a erdcdt h tnadinvariant standard the to reduced be can ewrs xo;nniereetoyais al Univer early electrodynamics; nonlinear 04.50.Kd axion; 04.40.Nr, Keywords: 04.40.-b, 04.20.-q, numbers: PACS and produ particles, particles. can charged accelerated flares born these electric the these of that acceleration anom expect, significant this We of flare. nonlin behavior axion term of the the the character of of The growth behavior anomalous Universe. the the magnetized is describing fields m result, gravitational and and the main electric apply The the and theory, by extended field. coordinated nonlinearly plane the of the equations in rotation el assoc the the symmetry to to prescribing discrete symmetry, Jackson’s the the first, second, field, symmetries: two for count xoial nue lcrcflrsi h al magnetized early the in flares electric induced Axionically mn sspotdietclyb h hieo i/oiemlilesi rn o front in multipliers sin/cosine of choice the by identically supported is ecnie h olnal xeddEinstein-Maxwell-ax extended nonlinearly the consider We lxne .Balakin, B. Alexander uprstedsrt ymty ic h ettr 2 term rest the since symmetry, discrete the supports F aa eea nvriy rmesaasr 6,Kzn420 Kazan 16a, str. Kremlevskaya University, Federal Kazan aa eea nvriy rmesaasr 6,Kzn420 Kazan 16a, str. Kremlevskaya University, Federal Kazan 1 ik eateto eea eaiiyadGaiain Institu Gravitation, and Relativity General of Department L I t dual its , ( I h arnincnb rnfre into transformed be can Lagrangian the ,where ), E ~ 2 ∗ eateto aipyis nttt fPhysics, of Institute Radiophysics, of Department olna xo electrodynamics: axion Nonlinear F → I ik ∗ L φ 1, n h xo field axion the and , E = ( ~ I ∗ nrdcdb eciadQin[]adidctdltra xo ed(s field axion as later indicated and [1] Quinn and Peccei by introduced cos ssm olna ucino t ruet n h invariant the and argument, its of function nonlinear some is ) 1 4 → I ldmrV Bochkarev, V. Vladimir (cos α f .INTRODUCTION I. + Dtd coe ,2021) 8, October (Dated: ( 4 1 φF B F φ ~ ( sin F mn ∗ mn mn α F F F , mn mn B mn ~ osntstsytemninddsrt ymty and symmetry, discrete mentioned the satisfy not does ) ∗ sin + + → φ se φF nrdcdnniery lal,tewl-nw linear well-known the Clearly, nonlinearly. introduced B ~ mn crdnmc ob nain ihrespect with invariant be to ectrodynamics F φ ∗ cos msino h lcrmgei ae by waves electromagnetic the of emission mercwt epc olna rnfrainof transformation linear to respect with ymmetric oteBaciImdlwt magnetic with model Bianchi-I the to m 2, srce ntewr 8.Tepooe nonlinear proposed The [8]. work the in nstructed L F oyais ehpe ob nprdb beauty by inspired be to happen we rodynamics, igrl nPyiso udmna Interactions. Fundamental of Physics in role ding clyidcdeeti edi h early the in field electric induced ically lu lcrcfil a eidctdby indicated be can field electric alous ae ihtepoete fteaxion the of properties the with iated ce ihteie htteLgaga fthe of Lagrangian the that idea the with ected mn α † eteeeto-oirnpi creation, pair electron-positron the ce o hoy hc sbsdo h ac- the on based is which theory, ion ∗ ( al ope xo,electromagnetic axion, coupled early mn ty soitdwt h shift the with associated etry, gei ed.W eietemaster the derive We fields. agnetic I − oew trc o osdrto h daof idea the consideration for attract we pose n liaF Nizamieva F. Albina and t h olna oicto ftetheory, the of modification nonlinear the ith aga,wihcnan h nainsbased invariants the contains which rangian, F ) πk E oueteLgaga,wihstse both satisfies which Lagrangian, the use to ~ ) I → mn sin √ ) − α ntelna hoy when theory, linear the In . . gF ihtecntn angle constant the with , eo Physics, of te mn ∗ 0,Russia 008, 0,Russia 008, .Cery h iceesymmetry discrete the Clearly, ]. F mn h tnadterms standard the f nteato ucinlis functional action the in universe 1, ‡ φ ˜ = I α φ φ nthe On . is F +2 → mn πk ,the 0, F (2) (1) mn ee, ( k KCL-PH-TH-2021-31

Fuelling the search for light dark matter–electron scattering

Louis Hamaide∗ and Christopher McCabe Theoretical Particle Physics and Cosmology Group, Department of Physics, King’s College London, Strand, London, WC2R 2LS, UK

Dark matter (DM) detectors employing a Spherical Proportional Counter (SPC) have demonstrated a single-electron detection threshold and are projected to have small background rates. We explore the sensitivity to DM–electron scattering with SPC detectors in the context of DarkSphere, a proposal for a 300 cm diameter fully-electroformed SPC. SPCs can run with different gases, so we investigate the sensitivity for five targets: helium, neon, xenon, methane, and isobutane. We use tools from quantum chemistry to model the atomic and molecular systems, and calculate the expected DM induced event rates. We find that DarkSphere has the potential to improve current exclusion limits on DM masses above 4 MeV by up to five orders of magnitude. Neon is the best all-round gas target but using gas mixtures, where methane and isobutane constitute 10% of the gas, can improve the sensitivity, especially when combined with helium. Our study highlights the currently untapped potential of SPCs to search for DM-electron scattering in the MeV-to-GeV DM mass range.

I. INTRODUCTION tion energy or band-gap energy, needs to be transferred to the electron to produce an observable signal. There al- Inferring more about the particle physics properties of ready exist an impressive range of existing constraints on dark matter (DM) remains one of the most important the DM-electron cross section in the MeV-to-GeV mass open challenges in astro-particle physics [1,2]. It is in range, ranging from detectors designed specifically for this context that the experimental program of DM direct this search channel as well as detectors that have been detection operates, since if a signal from DM is observed designed first and foremost for nuclear recoil searches but in a direct detection experiment, we will learn about the have been repurposed to search for DM-electron scatter- value of the DM mass and DM interaction with normal ing [46–54]. matter [3–5]. Traditionally, direct detection experiments It is in the latter context that we come to the main have focused on the search for nuclear recoils induced by topic of this paper. The New Experiments With Spheres- scattering from DM particles with the canonical WIMP Gas (NEWS-G) detectors employ a Spherical Propor- mass range of around 5-1000 GeV (see e.g. [6–8]). How- tional Counter (SPC) [55, 56] and can be filled with ever, the lack of any conclusive DM measurement has different species of atomic or molecular gases. Limits meant that during the past decade, the mass range of in- on nuclear-recoil signals from sub-GeV DM have already terest has broadened to encompass DM particles with a been placed with a 60 cm diameter SPC [57], while a mass down below and above this narrow mass window [9]. 140 cm diameter detector, SnoGlobe, which utilises The shift to searching for DM in the MeV to few-GeV ultra-radiopure copper to reduce the background rate, is mass range is challenging experimentally and a number being commissioned at SNOLAB [58]. Funding has also of different avenues are being pursued. The search for been secured for ECUME, a planned follow-up 140 cm nuclear recoils induced by sub-5-GeV DM has contin- diameter detector that will benefit from a fully electro- ued by using detectors with significantly smaller energy- formed vessel to lower the background rate, and a feasi- detection thresholds, often in combination with the util- bility study of hosting a 300 cm diameter fully electro- isation of lighter target nuclei (see e.g. [10–15]). Alter- formed SPC dubbed DarkSphere at the Boulby Un- natively, observable signals from sub-5-GeV DM can be derground Laboratory is underway [59]. arXiv:2110.02985v1 [hep-ph] 6 Oct 2021 generated from rare processes that involve the emission In its simplest form, the operation of the SPC is il- of photons or ‘Migdal’ electrons from the recoiling atom lustrated in fig.I. The detector consists of a grounded, (see e.g. [16–27]), or from the small flux of boosted DM spherical, metallic vessel filled with a gas mixture and a arising from interactions between DM and cosmic rays, spherical anode of radius ∼ 1 mm is placed at the centre. the Sun or mesons (see e.g. [28–34]). The anode produces an electric field that varies as r−2 A complementary approach is to search for DM- in an ideal SPC and allows the ionisation electrons pro- electron scattering processes that ionise an electron from duced through particle interactions in the gas volume to an atom or molecule, or excite an electron across the drift to the anode. When the electron reaches the an- band-gap in semi-conductors or other quantum materi- ode an avalanche occurs, providing signal amplification, als (see e.g. [35–45]). These searches have the advantage which is read out by a grounded metallic rod. Larger that only a small amount of energy, of order the ionisa- SPCs will employ multiple anodes (‘ACHINOS’ [60, 61]) but the basic principle of operation is similar to the simple SPC in fig.I. Owing to the low capacitance of the central sensor and ∗ [email protected] the high amplification gain near the anode, it has been Dimuons in Neutrino Telescopes: New Predictions and First Candidates in IceCube

Bei Zhou1, ∗ and John F. Beacom2, 3, 4, † 1Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA 2Center for Cosmology and AstroParticle Physics (CCAPP), Ohio State University, Columbus, Ohio 43210, USA 3Department of Physics, Ohio State University, Columbus, Ohio 43210, USA 4Department of Astronomy, Ohio State University, Columbus, Ohio 43210, USA (Dated: October 06, 2021) Neutrino telescopes are powerful probes of high-energy astrophysics and particle physics. Their power is increased when they can isolate different event classes, e.g., by flavor, though that is not the only possibility. Here we focus on a new event class for neutrino telescopes: dimuons, two energetic muons from one neutrino interaction. We make new theoretical and observational contributions. For the theoretical part, we calculate dimuon production cross sections and detection prospects via deep-inelastic scattering (DIS; where we greatly improve upon prior work) and W -boson production (WBP; where we present first results). We show that IceCube should have ' 400 dimuons (' 8 from WBP) in its current data and that IceCube-Gen2, with a higher threshold but a larger exposure, can detect ' 1200 dimuons (' 30 from WBP) in 10 years. These dimuons are almost all produced by atmospheric neutrinos. For the observational part, we perform a simple but conservative analysis of IceCube public data, finding the first candidate dimuon events (19 events). Though some IceCube experts we consulted argue these events cannot be real dimuons, (A) these events match well all aspects of our predictions and (B) no other compelling hypotheses have been raised. Whether these 19 events are real dimuons or some new background (or signal!), it is important to understand them. Here we share full details to help IceCube and to attract scrutiny from the broader community. Together, these theoretical and observational contributions help open a valuable new direction for neutrino telescopes, one especially important for probing high-energy QCD and new physics.

I. INTRODUCTION Coincident New single µ physics Neutrino observatories have been making ground- breaking discoveries [1–6]. For astrophysics, high-energy Detector neutrino telescopes in particular can reveal the elusive µ sites of hadronic cosmic-ray acceleration. IceCube has discovered TeV–PeV diffuse astrophysical neutrinos [6–9] Through-going µ and a few promising sources [10–12]. For particle physics, dimuon high-energy neutrino telescopes reach unprecedented en- νµ ergy scales, detecting some standard-model processes for the first time [13–24]. They also provide unique infor- νµ mation on new physics, including dark matter [25–33], Starting µ novel neutrino interactions [34–43], and other exotic sig- dimuon nals [44–51]. Much more progress is expected with new µ telescopes, including KM3NeT [52], Baikal-GVD [53], P- ONE [54], and especially IceCube-Gen2 [55]. So far, neutrino-telescope analyses are based on stud- FIG. 1. Examples of dimuons and other events with sim- arXiv:2110.02974v1 [hep-ph] 6 Oct 2021 ies of just a few event classes, including tracks from ilar signatures. Dimuon signals may be either starting and muons, showers from other particles, and double bangs throughgoing. Backgrounds from two coincident single muons from taus. To get more from the data, new event classes are rare except for downgoing events. New-physics signals, are needed. Dimuons are an event class that has been e.g., double staus produced in supersymmetric models [44– 50], can have a similar topology. detected in accelerator neutrino experiments at energies from tens to hundreds of GeV [56–62], but never in neutrino telescopes or in any experiment above 1 TeV. A dimuon event is two energetic muons produced by the the strange-quark PDF [63–65] through deep-inelastic same neutrino. Dimuon data are important for quan- scattering (DIS). Dimuons can also come from W-boson tum chromodynamics (QCD), especially for measuring production (WBP) [21], tridents [20, 21, 66–71], and the Glashow resonance [13, 23]. They are a background for new physics searches, e.g., where a pair of staus are produced [44–50]. Fig.1 sketches these possibilities. ∗ [email protected]; orcid.org/0000-0003-1600-8835 In this paper, we predict yields for, search for, and † [email protected]; orcid.org/0000-0002-0005-2631 discuss the physical significance of neutrino-induced Prepared for submission to JHEP DESY 21-153

Fuzzy Dark Matter Candidates from String Theory

Michele Cicolia,b Veronica Guidettic Nicole Righic Alexander Westphalc aDipartimento di Fisica e Astronomia, Università di Bologna, via Irnerio 46, 40126 Bologna, Italy bINFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy cDeutches Electronen-Synchrotron, DESY, Notkestraße 85, 22607 Hamburg, Germany E-mail: [email protected], [email protected], [email protected], [email protected]

Abstract: String theory has been claimed to give rise to natural fuzzy dark matter candidates in the form of ultralight axions. In this paper we revisit this claim by a detailed study of how moduli stabilisation affects the masses and decay constants of different axion fields which arise in type IIB flux compactifications. We find that obtaining a considerable contribution to the observed dark matter abundance without tuning the axion initial mis- alignment angle is not a generic feature of 4D string models since it requires a mild violation of the Sf . MP bound, where S is the instanton action and f the axion decay constant. Our analysis singles out C4-axions, C2-axions and thraxions as the best candidates to realise fuzzy dark matter in string theory. For all these ultralight axions we provide predictions which can be confronted with present and forthcoming observations.

Keywords: Fuzzy dark matter, ultralight axions, 4D string models arXiv:2110.02964v1 [hep-th] 6 Oct 2021 Prepared for submission to JCAP

Exploration of interacting dynamical dark energy model with interaction term including the equation-of-state parameter: alleviation of the H0 tension

Rui-Yun Guo,a,1 Lu Feng,b Tian-Ying Yao,a Xing-Yu Chena aSchool of Sciences, Xi’an Technological University, Xi’an 710021, China bCollege of Physical Science and Technology, Shenyang Normal University, Shenyang 110034, China E-mail: [email protected], [email protected], [email protected], [email protected]

Abstract. We explore a scenario of interacting dynamical dark energy model with the interaction term Q including the varying equation-of-state parameter w. Using the data combination of the cosmic microwave background, the baryon acoustic oscillation, and the type Ia supernovae, to global fit the interacting dynamical dark energy model, we find that adding a factor of the varying w in the function of Q can change correlations between the coupling constant β and other parameters, and then has a huge impact on the fitting result of β. In this model, the fitting value of H0 is lower at the 3.54σ level than the direct measurement value of H0. Comparing to the case of interacting dynamical dark energy model with Q excluding w, the model with Q including the constant w is more favored by the current mainstream observation. To obtain higher fitting values of H0 and narrow the discrepancy of H0 between different observations, additional parameters including the effective number of relativistic species, the total neutrino mass, and massive sterile neutrinos are considered in the interacting dynamical dark energy cosmology. We find that the H0 tension can be arXiv:2110.02536v1 [gr-qc] 6 Oct 2021 further reduced in these models, but is still at the about 3σ level.

1Corresponding author. Prepared for submission to JHEP

Scalar two-point functions at the late-time boundary of de Sitter

Gizem S¸eng¨ora Constantinos Skordisa aCEICO, Institute of Physics of the Czech Academy of Sciences, Na Slovance 1999/2, 182 21 Praha 8, Czechia E-mail: [email protected], [email protected]

Abstract: We calculate two-point functions of scalar ﬁelds of mass m and their conjugate momenta at the late-time boundary of de Sitter with Bunch-Davies boundary conditions, in general d + 1 spacetime dimensions. We perform the calculation using the wavefunction picture and using canonical quantization. With the latter one clearly sees how the late- time ﬁeld and conjugate momentum operators are linear combinations of the normalized

late-time operators αN and βN that correspond to unitary irreducible representations of the de Sitter group with well-deﬁned inner products. The two-point functions resulting from

these two different methods are equal and we find that both the autocorrelations of αN and βN and their cross correlations contribute to the late-time field and conjugate momentum d two-point functions. This happens both for light scalars (m < 2 H), corresponding to d complementary series representations, and heavy scalars (m > 2 H), corresponding to principal series representations of the de Sitter group, where H is the Hubble scale of de

Sitter. In the special case m = 0, only the βN autocorrelation contributes to the conjugate momentum two-point function. arXiv:2110.01635v1 [hep-th] 4 Oct 2021 Noname manuscript No. (will be inserted by the editor)

Transfer design between neighborhoods of planetary moons in the circular restricted three-body problem The Moon-to-Moon Analytical Transfer Method

David Canales · Kathleen C. Howell · Elena Fantino

Received: date / Accepted: date

Abstract Given the interest in future space missions devoted to the exploration of key moons in the solar system and that may involve libration point orbits, an efficient design strategy for transfers between moons is introduced that leverages the dynamics in these multi-body systems. The moon-to-moon analytical transfer (MMAT) method is introduced, comprised of a general methodology for transfer design between the vicinities of the moons in any given system within the context of the circular restricted three-body problem, useful regardless of the orbital planes in which the moons reside. A simplified model enables analytical constraints to efficiently determine the feasibility of a transfer between two different moons moving in the vicinity of a common planet. In particular, connections between the periodic orbits of such two different moons are achieved. The strategy is applicable for any type of direct transfers that satisfy the analytical constraints. Case studies are presented for the Jovian and Uranian systems. The transition of the transfers into higher-fidelity ephemeris models confirms the validity of the MMAT method as a fast tool to provide possible transfer options between two consecutive moons.

David Canales School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 47907 Tel.: +34-628-929-331 E-mail: [email protected] ORCID: 0000-0003-0166-2391 Kathleen C. Howell School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 47907 Tel.: +1-765-494-5786 E-mail: [email protected] ORCID: 0000-0002-1298-5017 Elena Fantino Aerospace Engineering Department, Khalifa University of Science and Technology, P.O. Box 127788, Abu Dhabi, United Arab Emirates Tel.: +971-2-312-4014 arXiv:2110.03683v1 [astro-ph.EP] 7 Oct 2021 E-mail: [email protected] ORCID: 0000-0001-7633-8567 LAPTH-038/21

3D template-based Fermi-LAT constraints on axion-like particles1

Francesca Calore,1, ∗ Pierluca Carenza,2, 3, 4, † Christopher Eckner,1, ‡ Tobias Fischer,5, § Maurizio Giannotti,6, ¶ Joerg Jaeckel,7, ∗∗ Kei Kotake,8 Takami Kuroda,9 Alessandro Mirizzi,3, 4, †† and Francesco Sivo3, ‡‡ 1Univ. Grenoble Alpes, USMB, CNRS, LAPTh, F-74000 Annecy, France 2The Oskar Klein Centre, Department of Physics, Stockholm University, Stockholm 106 91, Sweden 3Dipartimento Interateneo di Fisica “Michelangelo Merlin”, Via Amendola 173, 70126 Bari, Italy. 4Istituto Nazionale di Fisica Nucleare - Sezione di Bari, Via Orabona 4, 70126 Bari, Italy. 5Institute of Theoretical Physics, University of Wroc law, 50-204 Wroc law, Poland 6Physical Sciences, Barry University, 11300 NE 2nd Ave., Miami Shores, FL 33161, USA 7Institut fürtheoretische Physik, UniversitätHeidelberg, Philosophenweg 16, 69120 Heidelberg, Germany 8Faculty of Science, Department of Applied Physics & Research Institute of Stellar Explosive Phenomena, Fukuoka University, Fukuoka, Jonan 8-19-1, 814-0180, Japan 9Max-Planck-Institut fürGravitationsphysik, Am Mühlenberg 1, D-14476 Potsdam-Golm, Germany (Dated: October 8, 2021) Axion-like particles (ALPs) may be abundantly produced in core-collapse (CC) supernovae (SNe), hence the cumulative signal from all past SN events would contain an ALP component and create a diffuse flux peaked at energies of about 50 MeV. We update the calculation of this flux by including a set of CC SN models with different progenitor masses following the expected mass distribution. Additionally, we include the effects of failed CC SNe, which yield the formation of black holes instead of explosions. Relying on the coupling strength of ALPs to photons and the related Primakoff process, the diffuse SN ALP flux is converted into a diffuse gamma-ray flux while traversing the magnetic field of the Milky Way. The spatial morphology of this signal is expected to follow the shape of the Galactic magnetic field lines. We make use of this via a template-based analysis that utilizes 12 years of Fermi-LAT data in the energy range from 50 MeV to 500 GeV. This strategy yields an improvement of roughly a factor of two for the upper limit on the ALP-photon coupling constant gaγ compared to a previous analysis that accounted only for the spectral shape of the signal. While the improved SN modeling leads to a less energetic flux that is harder to detect, the combined effect is still an improvement of the limit and in particular its statistical reliability. We also show that our results are robust against variations in the modeling of high-latitude Galactic diffuse emission and systematic uncertainties of the LAT.

I. INTRODUCTION clusters [4, 6–9]. In the minimal scenario in which ALPs are coupled only with photons, the main channel for their Observing a supernova (SN) provides unique opportu- emissivity in the SN core is the Primakoff process, lead- nities for fundamental physics. In particular, they are ing to an ALP flux peaked at energies of about 50 MeV. ideally suited to probe feebly interacting particles (cf. [2] Conversion of these ALPs into gamma rays in the Milky for a recent review) with masses up to the 100 MeV Way magnetic field can lead to an observable gamma- range. Indeed, large numbers of such particles∼ can be ray burst in coincidence with the SN explosion [6, 7]. At emitted in SN events [3]. An important and theoretically the time of the SN 1987A, the Gamma-Ray Spectrome- interesting instance of this are axions and axion-like par- ter (GRS) on the Solar Maximum Mission (SMM) had ticles (ALPs) [4, 5]. Indeed, SN 1987A has significantly the SN in its field of view and observed no gamma-ray strengthened astrophysical axion bounds in a region of signal at the time of the explosion, this allowed to con- the yet-incompletely known ALP parameter space, com- strain the photon-ALP coupling early on [6, 7]. In a arXiv:2110.03679v1 [astro-ph.HE] 7 Oct 2021 plementary to the one probed by the Sun and the globular more refined, recent analysis, this upper limit is stated −12 −1 −10 as gaγ . 5.3 10 GeV for ma < 4 10 eV [10]. A future Galactic× SN explosion in the field× of view of the Large Area Telescope (LAT) aboard the Fermi satellite −13 −1 1 Parts of this work have already been shown at ICRC 2021 and a would allow us to constrain gaγ . 2.0 10 GeV −9 × short preview has appeared in the associated proceedings [1]. for ma < 10 eV [11]. Furthermore, a search for ∗ [email protected] gamma-ray bursts from extragalactic SNe with Fermi- † [email protected] −11 −1 LAT has yielded the limit gaγ . 2.6 10 GeV for ‡ [email protected] −10 × ma < 3 10 eV [12] (see also [13, 14]). § tobias.fi[email protected] × ¶ [email protected] While single SN explosions are rare and must fall into ∗∗ [email protected] the detector field of view to be observed, there exists a †† [email protected] guaranteed contribution to the gamma-ray diffuse flux ‡‡ [email protected] which originates from ALPs emitted by all past SNe in Variability in Quasar Light Curves: using quasars as standard candles

Rance Solomon∗ and Dejan Stojkovic HEPCOS, Department of Physics, SUNY at Buﬀalo, Buﬀalo, NY 14260-1500, USA

A relation between the variational slope, sF , and the mean ﬂux, hF i, in the light curves of 59 spectroscopically conﬁrmed quasars is measured with a dispersion of 0.14dex ranging over three orders of magnitude in hF i. Assuming it holds for quasars in general, not only does this relation add to our working knowledge of quasar variability but it may also be used to determine the luminosity distance of a quasar in a model independent way. Accurately obtaining the luminosity distance would allow quasars to be added to the cosmic distance ladder, easily extending the ladder out far beyond the redshifts accessible to type Ia supernovae where cosmological parameters can be better constrained.

I. INTRODUCTION in [3–6] through a nonlinear relation between the X-ray and UV luminosities. Although it had been previously known, the L −L relation had too large of dispersion The H0-tension is at the forefront of many cosmologi- X UV cal discussions today. Since the mere ≈ 8.5% difference (∼0.35dex) to make accurate predictions of the luminosity distance. However the authors were able to argue between the early and late time H0 measurements is sit- ting at a 4−6σ discrepancy we could be looking at either that much of the previously observed dispersion was not new physics or systematic errors in multiple observations. intrinsic to the quasars and could be subtracted out by We do not know of course which it is, but, in either case, choosing clean data sets, thereby reducing the dispersion the discrepancy is of great interest to the community and to 0.2dex. Their results show a tendency towards a time its study has instigated many new observational methods varying dark energy equation of state but requires further and theoretical possibilities. In this article we focus on a checking of systematics before concluding [6]. new late time observational method. The method we discuss in this article focuses on the The main approach to late time measurements is by measurement of the variational rate, sF , of the absolute building up the cosmic distance ladder in a similar man- light curves (i.e. the slopes of small time scale variations ner to the original work of Hubble in his calculation of in the rest frame of the quasar). We find that sF shares a the expansion rate. At present there are three common tight correlation with the mean absolute flux, hF i, with rungs on the distance ladder: geometrical distance mea- a strikingly low dispersion of 0.14dex for 59 spectroscop- surements, variable stars, and type 1a supernovae (SNe). ically confirmed quasars from the MACHO dataset. Fur- Each rung is anchored by the previous with SNe measur- thermore, this dispersion results from a minimal selection ing out at the farthest redshifts of z ≈ 0.15 as of the lat- criteria of the quasar light curves. By assuming this re- est SH0ES release [1] and an additional nine SNe in the lation holds for quasars in general, we show that one can range 1.5 < z < 2.3 from the follow up study [2]. These obtain the luminosity distance from the measured ap- distances to SNe combined with the recessional velocity parent light curves. This method originates in principle inferred by their redshift gives a measure of the expan- from [7] which had attempted to use quasars as standard sion rate throughout our recent history. And although clocks. local measurements are not able to directly constrain H0 We address the details of the data used for the study in they are able to directly constrain Ωm, today’s fractional section II while section III discusses the finding methods matter density, thereby reducing the shared parameter used to measure the variational rate which we apply to 59 space with H0. quasars from the MACHO survey. Section IV derives the Improvements on the distance ladder can be made by luminosity distance from the measured apparent fluxes adding more anchors, adding more SNe at greater red- and contains a demonstration using 304 under sampled shifts, or by adding a new rung of astronomical objects quasars from the Sloan Digital Sky Survey (SDSS). We arXiv:2110.03671v1 [astro-ph.CO] 7 Oct 2021 which can, on average, be observed at greater redshifts finish with conclusions and a brief discussion of future than SNe. For this last reason there have been many work in section V. attempts at including quasars into the distance ladder as they are the brightest objects known in the universe making them also among the most distant observable II. THE QUASAR SAMPLES objects. But much is still unknown about the behavior of quasars which has allowed high dispersion (intrinsic or otherwise) to plague their phenomenological uses. For reasons which will be discussed in section III, we Currently the most promising attempt is being made require the quasar light curves to be sampled with a relatively regular and high density. More than one observation per month is optimal for the analysis but sufficient results can still be obtained with less frequent observa- ∗ rancesol (at) buffalo.edu tions. This significantly reduces the number of available Astronomy & Astrophysics manuscript no. aanda ©ESO 2021 October 8, 2021

The Milky Way’s middle aged inner ring Shola M. Wylie1 ?, Jonathan P. Clarke1 , and Ortwin E. Gerhard1

Max-Planck-Institut fur Extraterrestrische Physik, Gießenbachstraße, D-85748 Garching, Germany

Received-; accepted -

ABSTRACT

Aims. We investigate the metallicity, age, and orbital anatomy of the inner Milky Way specifically focussing on the outer bar region. Methods. We integrate a sample of apogee DR16 inner Galaxy stars in a state of the art bar-bulge potential with a slow pattern speed and investigate the link between the resulting orbits and their [Fe/H] and ages. By superimposing the orbits we build density, [Fe/H], and age maps of the inner Milky Way, which we divide further using the orbital parameters eccentricity, |Xmax|, and |Zmax|. Results. We find that at low heights from the Galactic plane the Galactic bar gradually transitions into a radially thick, vertically thin, elongated inner ring with average solar [Fe/H]. This inner ring is mainly composed of stars with AstroNN ages between 4 and 9 Gyr with a peak in age between 6 and 8 Gyr, making the average age of the ring ∼6 Gyr. The vertical thickness of the ring decreases markedly towards younger ages. We also find very large L4 Lagrange orbits that have average solar to super-solar metallicities and intermediate ages. Lastly, we confirm a clear X-shape in the [Fe/H] and density distributions at large Galactic heights. Conclusions. The orbital structure obtained for the apogee stars reveals that the Milky Way hosts an inner ring-like structure between the planar bar and corotation. This structure is on average metal rich, intermediately aged, and enhances the horizontal gradient along the bar’s major axis. Key words. Galaxy: bulge, disk, structure, abundances, kinematics and dynamics, evolution

1. Introduction 2008; Skrutskie et al. 2006) by Wegg & Gerhard (2013); Wegg et al. (2015), and to kinematic data from the BRAVA, ARGOS, From previous Galactic bulge studies, we know that the stars in and OGLE surveys (Kunder et al. 2012; Ness et al. 2013; Rat- the outer regions of the Galactic bar are on average more metal tenbury et al. 2007). The favoured model from P17 had a slow rich and younger in comparison to the stars in the central bulge pattern speed Ωb which puts the corotation radius of the bar at (Bovy et al. 2019; Hasselquist et al. 2020; Queiroz et al. 2020). slightly over 6 kpc. More recent dynamical studies of both the This difference leads to a pronounced horizontal metallicity gra- stellar kinematics in the bar and resonant stars in the solar neigh- dient along the bar’s major axis (Wylie et al. 2021). The Milky bourhood have found similar or slightly slower values of Ωb (e.g. Way (MW) is not unique in this structure, bars with ends that are Bovy et al. 2019; Binney 2020; Chiba & Schönrich 2021; Li more metal rich and/or younger than their central regions have et al. 2021; Clarke & Gerhard 2021). been observed in a few other galaxies as well (Seidel et al. 2016), including M31 (Gajda et al. 2021). As we will see a mechanism to enhance the horizontal gra- This gradient structure has several possible origins, one of dient in the bar is an inner ring between 4-6 kpc along the bar’s which is that the bulge and bar were formed from dynamical major axis. This has also been seen in some galaxy models in instabilities of co-spatial discs with differing scale lengths and cosmological simulations (Fragkoudi et al. 2020). Inner rings are metallicities (Fragkoudi et al. 2018; Wylie et al. 2021). In this observed in a substantial fraction of barred disc galaxies (Kor- scenario the most metal poor discs dominate in the very central mendy 1979; Buta & Combes 1996; Comerón et al. 2014) and regions due to their shorter scale lengths while the more metal tend to be preferentially orientated with the bar. Inner rings ob- rich discs dominate in the outer regions, resulting in a positive served at optical wavelengths generally contain star-forming re- horizontal metallicity gradient. However the MW’s horizontal gions (Buta 1995) and are thought to form from the collection gradient is likely too steep to be fully explained by this scenario of gas around the bar’s 4:1 resonance (Schwarz 1984; Buta & suggesting that additional mechanisms are at work. Combes 1996). Passive rings which no longer form stars can be arXiv:2110.03658v1 [astro-ph.GA] 7 Oct 2021 To investigate this further we build density, [Fe/H], and age thicker than active rings but are generally seen in earlier type ≥ orbital maps for a sample of inner Galaxy stars from apogee galaxies ( Sab, Comerón 2013). DR16, for which these parameters along with their positional This paper is structured as follows: In Sect.2 we explain in and 3D kinematic information are available from the ASPCAP, more detail the sample of stars we used, the potential we inte- AstroNN, and Gaia DR2 catalogues (García Pérez et al. 2016; grated them in, and the assumptions we made. We also compare Gaia Collaboration et al. 2018; Leung & Bovy 2019; Mackereth the dynamical model from which we generate the potential to et al. 2019). We integrate the orbits of these stars in a realis- the apogee stars in heliocentric velocity space to check that they tic MW bar-bulge potential from Portail et al. (2017a, hereafter are consistent. In Sect.3 we show density maps for the model P17) which was fit to MW star count data derived from the VVV, and the apogee stellar orbits, and [Fe/H] and age maps for the UKIDSS, and 2MASS surveys (Saito et al. 2012; Lucas et al. apogee orbits divided by spatial symmetry and orbital eccentricity. Finally, in Sect.4 we end with a discussion of our results and ? [email protected] our conclusions.

Article number, page 1 of 8 Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX62

Massive Star Cluster Formation and Destruction in Luminous Infrared Galaxies in GOALS II: An ACS/WFC3 Survey of Nearby LIRGs

S. T. Linden,1 A. S. Evans,2, 3 K. Larson,4 G. C. Privon,3, 5 L. Armus,6 J. Rich,7 T. D´ıaz-Santos,8 E. J. Murphy,3, 2 Y. Song,2 L. Barcos-Munoz,˜ 3, 2 J. Howell,6 V. Charmandaris,8, 9 H. Inami,10 V. U,11 J. A. Surace,12 J. M. Mazzarella,6 and D. Calzetti1

1Department of Astronomy, University of Massachusetts at Amherst, Amherst, MA 01003, USA. 2Astronomy Department, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904 USA 3National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903 USA 4Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218, USA 5Department of Astronomy, University of Florida, 211 Bryant Space Sciences Center, Gainesville, FL 32611 6IPAC, California Institute of Technology, MS 100-22, Pasadena, CA 91125 USA 7Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91101 8Department of Physics, University of Crete, GR-71003, Heraklion, Greece 9Institute of Astrophysics, FORTH, GR-71110 Heraklion, Greece 10Hiroshima Astrophysical Science Center, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan 11Department of Physics and Astronomy, University of California, Irvine, 4129 Frederick Reines Hall, Irvine, CA 92697, USA 12Eureka Scientiﬁc, Inc. 2452 Delmer Street Suite 100 Oakland, CA 94602-3017 USA

(Dated: October 8, 2021)

ABSTRACT We present the results of a Hubble Space Telescope WFC3 near-UV and ACS/WFC optical study into the star cluster populations of 10 luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with dL < 110 Mpc in order to avoid issues related to cluster blending. The measured cluster age distribution slope of dN/dτ ∝ τ −0.5+/−0.2 is steeper than what has been observed in lower-luminosity star-forming galaxies. Further, differences in the slope of the observed cluster age distribution between inner- (dN/dτ ∝ τ −1.07+/−0.12) and outer-disk (dN/dτ ∝ τ −0.37+/−0.09) star clusters provides evidence of mass-dependent cluster destruction in the central regions of LIRGs driven primarily by the combined effect of strong tidal shocks and encounters with massive GMCs. Excluding the nuclear ring surrounding the Seyfert 1 nucleus in NGC 7469, the derived cluster mass function (CMF: dN/dM ∝ M α) has 6 marginal evidence for a truncation in the power-law (PL) at Mt ∼ 2x10 M for our three most cluster- rich galaxies, which are all classified as early-stage mergers. Finally, we find evidence of a flattening of the CMF slope of dN/dM ∝ M −1.42±0.1 for clusters in late-stage mergers relative to early-stage (α = −1.65 ± 0.02), which we attribute to an increase in the formation of massive clusters over the course of the interaction.

arXiv:2110.03638v1 [astro-ph.GA] 7 Oct 2021 Keywords: galaxies: star clusters: general - galaxies: starburst - galaxies: ISM - infrared: galaxies

6 1. INTRODUCTION (SSCs) with masses which can reach ≥ 10 M ; i.e., sim- Imaging surveys of high-redshift star-forming galaxies ilar to the most massive known globular clusters and 1-2 reveal that these systems tend to display turbulent disks, orders of magnitude more massive than any young clus- which host star-forming clumps with stellar masses of ter or HII region found in nearby normal galaxies like 7−8 the Milky Way (MW) or M51 (Kobulnicky & Johnson ∼ 10 M and sizes of 0.5 − 5 kpc (Elmegreen et al. 2004, 2009; Daddi et al. 2010; Livermore et al. 2015). 1999; Whitmore et al. 2014; Cook et al. 2016; Dessauges- These extreme clumps are ∼ 100x more massive than Zavadsky & Adamo 2018). the typical giant molecular cloud (GMCs) observed in While the total number of massive clusters formed is the local Universe, and give rise to super star clusters largely set by the gas density in the interstellar medium

Asteroid 16 Psyche: Shape, Features, and Global Map

Michael K. Shepard12, Katherine de Kleer3, Saverio Cambioni3, Patrick A. Taylor4, Anne K. Virkki4, Edgard G. Rívera-Valentin4, Carolina Rodriguez Sanchez-Vahamonde4, Luisa Fernanda Zambrano-Marin5, Christopher Magri5, David Dunham6, John Moore7, Maria Camarca3

Abstract We develop a shape model of asteroid 16 Psyche using observations acquired in a wide range of wavelengths: Arecibo S-band delay-Doppler imaging, Atacama Large Millimeter Array (ALMA) plane-of-sky imaging, adaptive optics (AO) images from Keck and the Very Large Telescope (VLT), and a recent stellar occultation. Our shape model has dimensions 278 (–4/+8 km) x 238(–4/+6 km) x 171 km (–1/+5 km), an effective spherical diameter Deff = 222 -1/+4 km, and a spin axis (ecliptic lon, lat) of (36°, -8°) ± 2°. We survey all the features previously reported to exist, tentatively identify several new features, and produce a global map of Psyche. Using 30 calibrated radar echoes, we find Psyche’s overall radar albedo to be 0.34 ± 0.08 suggesting that the upper meter of regolith has a significant metal (i.e., Fe-Ni) content. We find four regions of enhanced or complex radar albedo, one of which correlates well with a previously identified feature on Psyche, and all of which appear to correlate with patches of relatively high optical albedo. Based on these findings, we cannot rule out a model of Psyche as a remnant core, but our preferred interpretation is that Psyche is a differentiated world with a regolith composition analogous to enstatite or CH/CB chondrites and peppered with localized regions of high metal concentrations. The most credible formation mechanism for these regions is ferrovolcanism as proposed by Johnson et al. [Nature Astronomy vol 4, January 2020, 41-44.].

1Corresponding author [email protected] 2Bloomsburg University, 400 E. Second St., Bloomsburg, PA 17815, USA 3California Institute of Technology, Pasadena, CA 91125, USA 4Arecibo Observatory, University of Central Florida, Arecibo, PR 00612, USA 5University of Maine Farmington, Farmington, ME, 04938 USA. 6International Occultation Timing Association, Greenbelt, MD 20768-0423, USA

Astronomy & Astrophysics manuscript no. main ©ESO 2021 October 8, 2021

The traditional approximation of rotation for rapidly rotating stars and planets II. Deformation and differential rotation

H. Dhouib1, V. Prat1, T. Van Reeth2, and S. Mathis1

1 Département d’Astrophysique-AIM, CEA/DRF/IRFU, CNRS/INSU, Université Paris-Saclay, Université Paris-Diderot, Université de Paris, F-91191 Gif-sur-Yvette, France e-mail: [email protected] 2 Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium

October 8, 2021

ABSTRACT

Context. We examine the dynamics of low-frequency gravito-inertial waves (GIWs) in differentially rotating deformed radiation zones in stars and planets by generalising the traditional approximation of rotation (TAR). The TAR treatment was built on the assumptions that the star is spherical (i.e. its centrifugal deformation is neglected) and uniformly rotating. However, it has been generalised in our previous work by including the effects of the centrifugal deformation using a non-perturbative approach. In the meantime, TAR has been generalised in spherical geometry to take the differential rotation into account. Aims. We aim to carry out a new generalisation of the TAR treatment to account for the differential rotation and the strong centrifugal deformation simultaneously. Methods. We generalise our previous work by taking into account the differential rotation in the derivation of our complete analytical formalism that allows the study of the dynamics of GIWs in differentially and rapidly rotating stars. Results. We derived the complete set of equations that generalises the TAR, simultaneously taking the full centrifugal acceleration and the differential rotation into account. Within the validity domain of the TAR, we derived a generalised Laplace tidal equation for the horizontal eigenfunctions and asymptotic wave periods of the GIWs, which can be used to probe the structure and dynamics of differentially rotating deformed stars with asteroseismology. Conclusions. A new generalisation of the TAR, which simultaneously takes into account the differential rotation and the centrifugal acceleration in a non-perturbative way, was derived. This generalisation allowed us to study the detectability and the signature of the differential rotation on GIWs in rapidly rotating deformed stars and planets. We found that the effects of the differential rotation in early-type deformed stars on GIWs is theoretically largely detectable in modern space photometry using observations from Kepler and TESS. Key words. hydrodynamics - waves - stars: rotation - stars: oscillations - methods: analytical - methods: numerical

1. Introduction 2019; Deheuvels et al. 2020). Then, observations of oscillation modes in F- and A-type stars (Kurtz et al. 2014; Saio et al. 2015; Bedding et al. 2015; Keen et al. 2015; Van Reeth et al. 2015, Understanding how angular momentum and chemicals are trans- 2016; Schmid & Aerts 2016; Murphy et al. 2016; Sowicka et al. ported in the interiors of stars (and planets) along their evolution 2017; Guo et al. 2017; Van Reeth et al. 2018; Saio et al. 2018; is one of the key challenges of modern stellar (and planetary) as- Mombarg et al. 2019; Li et al. 2019, 2020; Ouazzani et al. 2020; trophysics. Indeed, rotation modifies their structure, their chem- Saio et al. 2021) and in B-type stars (Pápics et al. 2015; Triana ical stratification, their internal flows and magnetism, and their et al. 2015; Moravveji et al. 2016; Pápics et al. 2017; Kallinger mass losses and winds (e.g. Maeder 2009; Mathis 2013; Aerts et al. 2017; Buysschaert et al. 2018; Szewczuk & Daszynska-´ et al. 2019, and references therein). In this quest, asteroseismol- Daszkiewicz 2018; Pedersen et al. 2021; Szewczuk et al. 2021) arXiv:2110.03619v1 [astro-ph.SR] 7 Oct 2021 ogy has bought a fundamental breakthrough by demonstrating provided us new Rosetta stones to constrain the transport of that all stars are the seat of a strong extraction of angular mo- angular momentum in the whole Hertzsprung-Russell diagram. mentum during their evolution in comparison with the predic- More particularly, this pushes gravity- and gravito-inertial mode tions by stellar models taking the rotation into account follow- pulsators such as γ-Doradus and SPB stars at the forefront of ing the standard rotational transport and mixing theory (Eggen- this research. For instance, recent theoretical developments have berger et al. 2012; Marques et al. 2013; Ceillier et al. 2013; demonstrated how it is possible to probe stellar internal rotation Cantiello et al. 2014; Ouazzani et al. 2019). This was first ob- in γ-Doradus stars from their surface to their convective core tained thanks to mixed pulsation modes splitted by rotation prop- (Ouazzani et al. 2020; Saio et al. 2021). These stars are rapid agating in evolved low- and intermediate-mass stars (Beck et al. rotators for a large proportion of them. Therefore, it is necessary 2012; Mosser et al. 2012; Deheuvels et al. 2012; Beck et al. to study gravito-inertial modes. These modes are gravity modes, 2014; Deheuvels et al. 2014, 2015; Di Mauro et al. 2016; Tri- which propagate only in stably stratified stellar radiation zones ana et al. 2017; Gehan et al. 2018; Beck et al. 2018; Tayar et al.

Article number, page 1 of 13 Anomalous Cosmic Ray Oxygen Observations in to 0.1 au

J. S. Rankin1, D. J. McComas1, R. A. Leske2, E. R. Christian3, C. M. S. Cohen2, A. C. Cummings2, C. J. Joyce1, A. W. Labrador2, R. A. Mewaldt2, N. A. Schwadron1,4, E. C. Stone2, R. D. Strauss5, M. E. Wiedenbeck6

1 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08540, USA 2 California Institute of Technology, Pasadena, CA 91125, USA 3 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 4 University of New Hampshire, Durham, NH 03824, USA 5 Center for Space Research, North-West University, Potchefstroom, South Africa 6 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA Abstract

The Integrated Science Investigation of the Sun instrument suite onboard NASA’s Parker Solar Probe mission continues to measure solar energetic particles and cosmic rays closer to the Sun than ever before. Here, we present the first observations of cosmic rays into 0.1 au (21.5 solar radii), focusing specifically on oxygen from ~2018.7 to ~2021.2. Our energy spectra reveal an anomalous cosmic ray-dominated profile that is comparable to that at 1 au, across multiple solar cycle minima. The galactic cosmic ray-dominated component is similar to that of the previous solar minimum (Solar Cycle 24/25 compared to 23/24) but elevated compared to the past (Solar Cycle 20/21). The findings are generally consistent with the current trend of unusually weak solar modulation that originated during the previous solar minimum and continues today. We also find a strong radial intensity gradient: 49.4 ± 8.0 % au-1 from 0.1 to 0.94 au, for energies of 6.9 to 27 MeV nuc-1. This value agrees with that measured by Helios nearly 45 years ago from 0.3 to 1.0 au (48 ± 12 % au-1; 9 to 29 MeV nuc-1) and is larger than predicted by models. The large ACR gradients observed close to the Sun by the Parker Solar Probe Integrated Science Investigation of the Sun instrument suite found here suggest that intermediate-scale variations in the magnetic field’s structure strongly influences cosmic ray drifts, well inside 1 au. 1. Introduction

The study of the spatial and temporal variations of anomalous cosmic rays (ACRs) has been ongoing for nearly 50 years and continues to provide compelling insight about the nature of our heliosphere and the acceleration and transport of energetic particles. ACRs consist mainly of singly-ionized atomic species of high first-ionization potential elements – such as hydrogen, helium, nitrogen, oxygen, neon, and argon – typically in the ~5 to ~50 MeV nuc-1 energy range (e.g., Garcia-Munoz et al. 1973; Hovestadt et al. 1973; McDonald et al. 1974; Christian et al. 1988; Klecker et al. 1998; Cummings et al. 2002a, 2002b; Potgieter 2013). Originating as interstellar neutrals, these particles drift into the heliosphere as a part of the interstellar wind at speeds of ~25.4 km s-1 (McComas et al. 2015). On their journey inward toward the Sun, some fraction gets ionized (Fisk et al. 1974), forming “pickup ions” which gain energy in the motional electric field as they are carried back outward as a part of the solar wind (e.g., Möbius et al. 1985; Drews et al. 2016). A tiny portion of these ~keV pickup ions are accelerated in the outer heliosphere, where they gain energies of a few to a hundred MeV nuc-1 over timescales of a year or less (Jokipii 1996; Mewaldt et al. 1996; Barghouty et al. 2000; Giacalone et al. 2012). The newly accelerated particles form a new population at the low-energy end of the cosmic ray spectra, known as “anomalous cosmic rays”. The source and mechanism for ACR acceleration is The smearing scale in Laguerre reconstructions of the correlation function

Farnik Nikakhtar1,∗ Ravi K. Sheth1,2, and Idit Zehavi3 1Department of Physics and Astronomy, University of Pennsylvania, 209 S. 33rd St., Philadelphia, PA 19104 – USA 2Center for Particle Cosmology, University of Pennsylvania, 209 S. 33rd St., Philadelphia, PA 19104 – USA and 3Department of Physics, Case Western Reserve University, Cleveland, OH 44106-7079 – USA (Dated: October 8, 2021) To a good approximation, on large cosmological scales the evolved two-point correlation function of biased tracers is related to the initial one by a convolution. For Gaussian initial conditions, the smearing kernel is Gaussian, so if the initial correlation function is parametrized using simple polynomials then the evolved correlation function is a sum of generalized Laguerre functions of half- integer order. This motivates an analytic ‘Laguerre reconstruction’ algorithm which previous work has shown is fast and accurate. This reconstruction requires as input the width of the smearing kernel. We show that the method can be extended to estimate the width of the smearing kernel from the same dataset. This estimate, and associated uncertainties, can then be used to marginalize over the distribution of reconstructed shapes, and hence provide error estimates on the value of the distance scale which are not tied to a particular cosmological model. We also show that if, instead, we parametrize the evolved correlation function using simple polynomials, then the initial one is a sum of Hermite polynomials, again enabling fast and accurate deconvolution. If one is willing to use constraints on the smearing scale from other datasets, then marginalizing over its value is simpler for Hermite reconstruction, potentially providing further speed-up in cosmological analyses.

I. INTRODUCTION to ask if the data constrain this scale independently – i.e., without having to assume a fiducial cosmological model. On large (∼ 100h−1Mpc) cosmological scales, the This is particularly interesting because the uncertainties evolved two-point correlation function, even of unbiased on this estimate can be propagated into uncertainties on tracers of the density field, differs in shape from the the shape of the reconstructed correlation function. In unbiased linear theory correlation function [1, 2]. To turn, this provides a means to marginalize over the value leading order, this change in shape is due to a (three- of the smearing scale when estimating cosmological pa- dimensional) convolution [3]. For Gaussian initial condi- rameters – such as the cosmological distance scale at the tions, the convolution kernel is very well approximated redshift of the survey – that is not tied to a fiducial model. by a Gaussian. Convolving a polynomial of order n with In Section II, we show that this is indeed possible in a Gaussian yields a generalized Laguerre function of or- principle, but whether or not the constraint which results der n/2, and [4] used this to motivate an algorithm for is sufficiently tight to be interesting (better than ∼ 10% estimating the initial shape from measurements of the precision) will depend on the data set (volume and tracer evolved one. Essentially, fitting a series of half-integer number density). We illustrate the method using the real Laguerre functions to the two-point correlation function space dark matter correlation function, before applying of biased tracers of the cosmological density field allows it to the monopole of the redshift space distorted corre- one to perform the deconvolution analytically, even in lation function measured in mock galaxy catalogs. the presence of redshift space distortions [5]. Section III shows how this fiducial model-free determi- This reconstruction of the initial shape requires as in- nation of the smearing scale can be used to estimate the put a guess for the width of the Gaussian smearing kernel cosmological distance scale and, in particular, to provide Σ. In CDM models, the physics which gives rise to the realistic estimates of the precision of this determination. convolution relates the smearing scale to the background If the smearing scale is not well-constrained, then this arXiv:2110.03591v1 [astro-ph.CO] 7 Oct 2021 cosmological model: decreases the precision of the distance scale constraint. Therefore, especially for small surveys, it may be that Z 2 2 other datasets provide more useful determinations of the Σ = dk PL(k, z)/3π , (1) smearing. We argue that this motivates consideration of other parametrizations of the reconstruction problem.

where PL(k, z) is the linear theory power spectrum [1] at E.g., even if useful constraints on the smearing scale must redshift z. So, if the model parameters are suﬃciently come from independent datasets, some parametrizations well-known, then they can be used to provide a useful of the reconstruction problem may allow one to marginal- estimate of the smearing scale. However, it is interesting ize over the value of the smearing scale more easily than others. These alternative parametrizations are the subject of Section IV. A ﬁnal section summarizes our results. We illustrate our arguments using the dark matter ∗Electronic address: [email protected] correlation functions measured in the z = 1 outputs

The Effect of Age on the Grouping of Open Clusters: The Primordial Group Hypothesis

Juan Casado

Facultad de Ciencias, Universidad Autónoma de Barcelona, 08193, Bellaterra, Catalonia, Spain

Email: [email protected]

Abstract

The Primordial group hypothesis states that only young enough open clusters (OCs) can be multiple, and old OCs are essentially single since the gravitational interaction between OCs in primordial groups is very weak, and the probability of gravitational capture of two OCs without disruption is very low. We test such postulate through four different studies using a manual search of Gaia EDR3 and extensive literature. First, we revisit the work of de La Fuente Marcos & de La Fuente Marcos (2009), which states that only ca. 40% of the OCs pairs are of primordial origin. However, no plausible binary system among their proposed OC pairs having at least one member older than 0.1 Gyr has been found. Second, we research among the youngest OCs (age < 0.01 Gyr) in Tarricq et al. (2021) and obtain that ca. 71% of them remain in their primordial groups. Third, a similar study of the oldest OCs (age > 4 Gyr) shows that they are essentially alone. Forth, the well-known case of the double cluster in Perseus is shown also to accommodate the title hypothesis. A simplified bimodal model allows retrieving the overall fraction of linked OCs (around 12-16 %) from our results, assuming that young clusters remain associated ~0.04 Gyr. The obtained results further support that OCs are born in groups (Casado 2021).

Keywords

Binary Open Clusters; Open Cluster gGroups; Open Cluster Formation; Gaia; Manual Search; Primordial Group Hypothesis

1. Introduction

Open clusters (OCs) are born from the gravitational collapse of gas and dust in giant molecular clouds. There is observational evidence that some of them are born in groups (Camargo et al., 2016). Galactic OCs are sometimes found in pairs, and the number of these optical pairs is significantly higher than would be expected if clusters were randomly distributed (e.g. Rozhavskii et al. 1976). Studies of grouping among OCs provide keys to understanding star formation in the Galactic disk and the subsequent dynamical evolution of OCs.

Until recently, h and χ Persei were the only confirmed physical double cluster in our Galaxy (Vázquez et al. 2010), despite the existing literature on OC grouping was already significant (Subramanian et al. 1995; Piskunov et. al. 2006; de La Fuente Marcos & de La Fuente Marcos 2009). Conversely, roughly 10% of the known OCs in the Large Magellanic Cloud (LMC) seem to belong to pairs (Vázquez et al., 2010). Back to the Galactic disk, the first estimations of the fraction of these binary clusters reached a maximum level of 20 % (Rozhavskii et al. 1976). Subramanian et al. (1995) estimated that about 8% of the OCs in the Galaxy appear to be members of binary systems. De La Fuente Marcos & de La Fuente Marcos (2009) argued that the real fraction is ca. 12%, similar to that in the LMC. However, the exact fraction remains an open question at present, which new precision Gaia data will hopefully ascertain in the next future, at least in the solar neighbourhood.

Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX631

Testing Models of Triggered Star Formation with Young Stellar Objects in Cepheus OB4

Abby Mintz,1, 2 Joseph L. Hora,2 and Elaine Winston2

1Department of Astronomy, Yale University, 52 Hillhouse Ave., New Haven, CT 06511, USA 2Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., MS-65, Cambridge MA 02138, USA

ABSTRACT OB associations are home to newly formed massive stars, whose turbulent winds and ionizing flux create H II regions rich with star formation. Studying the distribution and abundance of young stellar objects (YSOs) in these ionized bubbles can provide essential insight into the physical processes that shape their formation, allowing us to test competing models of star formation. In this work, we examined one such OB association, Cepheus OB4 (Cep OB4) – a well-suited region for studying YSOs due to its Galactic location, proximity, and geometry. We created a photometric catalog from Spitzer/IRAC mosaics in bands 1 (3.6 µm) and 2 (4.5 µm). We supplemented the catalog with photometry from WISE, 2MASS, IRAC bands 3 (5.8 µm) and 4 (8.0 µm), MIPS 24 µm, and MMIRS near IR data. We used color-color selections to identify 821 YSOs, which we classified using the IR slope of the YSOs’ spectral energy distributions (SEDs), finding 67 Class I, 103 flat spectrum, 569 Class II, and 82 Class III YSOs. We conducted a clustering analysis of the Cep OB4 YSOs and fit their SEDs. We found many young Class I objects distributed in the surrounding shell and pillars as well as a relative age gradient of unclustered sources, with YSOs generally decreasing in age with distance from the central cluster. Both of these results indicate that the expansion of the H II region may have triggered star formation in Cep OB4.

Keywords: Young stellar objects; Star formation; Star forming regions

1. INTRODUCTION is difficult due to dust extinction, their relatively low Many stars in our galaxy formed in OB associa- numbers and greater distances, and the relative brevity tions, loose unbound systems of massive O and B stars of their formation period. (Briceñoet al. 2007; Lada 1999). The environments of Observing their effects on less-massive stars’ presence OB associations, typically embedded in giant molecular and development in OB associations is more feasible, in clouds (GMCs), are strongly affected by massive stars part due to the creation of ionized H II regions and the longer life cycle of less-massive stars. The stellar winds (> 8 M ) whose turbulent and fast-paced lives dramatically shape their surroundings (Zinnecker & Yorke and ionizing flux from massive stars form H II regions 2007). Such OB associations also host less massive star within GMCs, pushing material from inside the OB as- formation and therefore young, less massive stars. Mas- sociation to the outskirts of the H II regions. These sive stars are an integral part of galactic evolution, pro- ionized bubbles have borders of dense material, nebu- arXiv:2110.03559v1 [astro-ph.SR] 7 Oct 2021 viding heavy metals, massive outflows, and turbulent lous features, and filaments, all of which can be sites winds that greatly influence their environments (Zin- of lower-mass star formation. Many star-forming H II necker & Yorke 2007). Despite their important role in regions with bubble morphologies have been studied pre- shaping astrophysical conditions, the processes that lead viously, including Sh2-236 (Ortega et al. 2020), Sh2-48 to their formation are still not entirely understood. Di- (Ortega et al. 2013), and Vela OB2 (Cantat-Gaudin et rectly observing O and B stars in the midst of formation al. 2019) among others (Xu et al. 2014). The presence, distribution, and location of YSOs provides information about recent star formation in the Corresponding author: Abby Mintz OB associations and can help determine what physical [email protected] processes are at play. The collect and collapse model posits that winds and supernovae caused by massive Review Space-Based Photometry Of Binary Stars: From Voyager To TESS

John Southworth1

1 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK; [email protected]

Abstract: Binary stars are crucial laboratories for stellar physics, so have been photometric targets for space missions beginning with the very first orbiting telescope (OAO-2) launched in 1968. This review traces the binary stars observed and the scientific results obtained from the early days of ultraviolet missions (OAO-2, Voyager, ANS, IUE), through a period of diversification (Hipparcos, WIRE, MOST, BRITE), to the current era of large planetary transit surveys (CoRoT, Kepler, TESS). In this time observations have been obtained of detached, semi-detached and contact binaries containing dwarfs, sub-giants, giants, supergiants, white dwarfs, planets, neutron stars and accretion discs. Recent missions have found a huge variety of objects such as pulsating stars in eclipsing binaries, multi-eclipsers, heartbeat stars and binaries hosting transiting planets. Particular attention is paid to eclipsing binaries, because they are staggeringly useful, and to the NASA Transiting Exoplanet Survey Satellite (TESS) because its huge sky coverage enables a wide range of scientific investigations with unprecedented ease. These results are placed into context, future missions are discussed, and a list of important science goals is presented.

Keywords: stars: binaries; stars: fundamental parameters; stars: oscillations; techniques: photometric

1. Introduction In 1923 Hermann Oberth [1] pointed out that telescopes placed above Earth’s atmosphere would have major advantages versus those that remain on the ground. These include avoiding the blurring effects of the atmosphere and being able to observe at wavelengths at which the atmosphere is opaque. To this we can add the avoidance of scintillation, which causes noise in astronomical photometry [2,3], the ability to observe during daytime, no need to worry about bad weather [4] or adverse ground conditions1, and isolation from wildlife2. Siting a telescope

Citation: Southworth, J. Space-Based in space therefore enables one to obtain data of higher quality, with minimal interruptions, Photometry Of Binary Stars: From and at wavelengths inaccessible from the ground. Voyager To TESS. Preprints 2021, 1, 0. There are obvious disadvantages of space telescopes, including the cost and complexity of https://doi.org/ constructing them and putting them in orbit, and the inability to perform on-site maintenance and technical interventions. Nevertheless, space-based telescopes can be the best tool with Received: which to tackle a range of scientific goals. This is particularly true for observing programs that Accepted: require long-term and/or high-precision photometry of a large number of stars. This review arXiv:2110.03543v1 [astro-ph.SR] 7 Oct 2021 Published: is aimed at summarising past work on binary stars using photometry from space telescopes, with a particular emphasis on eclipsing binaries (EBs) and the scientific results coming from Publisher’s Note: MDPI stays neutral the Transiting Exoplanet Survey Satellite (TESS) mission [5]. with regard to jurisdictional claims in published maps and institutional affil- 1.1. Importance of binary stars iations. The study of EBs is a crucial component in our understanding of stellar physics. It is possible to measure the masses and radii of the components of an EB from directly observed quantities [light curves of eclipses and radial velocity (RV) measurements from spectra]

1 More than one ground-based telescope has been permanently shut down by a lightning strike, and Mt. Stromlo Observatory was destroyed by a bushfire. 2 Ground-based telescopes have been temporarily closed for, among other things, ant infestations, venomous snakes, and mice making bedding out of the cladding of fibre optic cables. MNRAS 000,1–16 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

Relativistic scattering of a fast spinning neutron star by a massive black hole

Kaye Jiale Li1,2★, Kinwah Wu1, Po Kin Leung2 and Dinesh Singh3 1Mullard Space Science Laboratory, University College London, Holmbury St Mary, Surrey, RH5 6NT, United Kingdom 2Department of Physics, Chinese University of Hong Kong, Shatin, NT, Hong Kong SAR, China 3Department of Physics, University of Regina, Regina, SK S4S 0A2, Canada

Accepted 2021 October 5. Received 2021 October 4; in original form 2021 February 3

ABSTRACT The orbital dynamics of fast spinning neutron stars encountering a massive Black Hole (BH) with unbounded orbits are investigated using the quadratic-in-spin Mathisson-Papapetrou-Dixon (MPD) formulation. We consider the motion of the spinning neutron stars with astrophysically relevant speed in the gravity field of the BH. For such slow-speed scattering, the hyperbolic orbits followed by these neutron stars all have near the 푒 = 1 eccentricity, and have distinct properties compared with those of 푒 1. We have found that compared with geodesic motion, the spin-orbit and spin-spin coupling will lead to a variation of scattering angles at spatial infinity, and this variation is more prominent for slow-speed scattering than fast-speed scattering. Such a variation leads to an observable difference in pulse-arrival-time within a few hours of observation, and up to a few days or months for larger BH masses or longer spinning periods. Such a relativistic pulsar-BH system also emits a burst of gravitational waves (GWs) in the sensitivity band of LISA, and for optimal settings, can be seen up to 100 Mpc away. A radio follow up of such a GW burst with SKA or FAST will allow for measuring the orbital parameters with high accuracy and testing the predictions of General Relativity (GR). Key words: black hole physics – gravitation – celestial mechanics – relativistic processes – pulsars general

1 INTRODUCTION important system parameter in the orbital and spin dynamical modelling. Moreover, as the mass of a neutron star is small in comparison Pulsars are either fast spinning young neutron stars with a strong with the massive nuclear BH, the MSP can be treated as a test object. magnetic field (∼ 1012 − 1013 G) or old recycled neutron stars with In the test-object limit, the motion of a spinning secondary in a weaker magnetic field (∼ 109 G) (see e.g. Lorimer 2008). The lat- the gravity field of a massive black hole is governed by the MPD ter are millisecond pulsars (MSPs), whose spin periods range from equations. Most of the current studies (e.g. Semerák 1999; Yunes ∼ 1 − 10ms and are extremely stable, with a drift of roughly 1 pulse et al. 2011; Singh et al. 2014; Kimpson et al. 2020a,b) using MPD period on a Hubble timescale (see e.g. Manchester 2017). MSPs are equations have put focus on bounded systems and their orbital dy- mass gyros, and the stability of their spin periods make them very namics. Also a circular or a quasi-circular orbit approximation is high-precision timing instruments. Relativistic binaries containing often adopted (e.g. Bini et al. 2005; Han 2010; Velandia & Tejeiro an MSP orbiting around a BH, i.e. MSP - BH binaries, are par- 2018; Chen et al. 2019), where the spin-spin and spin-orbit couplings ticularly useful for the investigation of fundamental physics, where are treated as independent separated components (see also the appli- accurate, reliable measurements are essential. The two subclasses cation of same deposition in the study of binary neutron stars around of the MSP - BH binaries, the extreme-mass-ratio binary (EMRB) arXiv:2110.03494v1 [astro-ph.HE] 7 Oct 2021 a black hole in Remmen & Wu 2013). Compared to the bounded and the extreme-mass-ratio-inspiral (EMRI) systems, each contain orbits, hyperbolic orbits admit the advantage of having clean non- a massive BH. These systems are natural GW sources, expected degenerate observables. Among them the most well recognised is the to be detected by the Laser Interferometric Space Antenna (LISA) scattering angle 휒 (≡ Δ휙 in this paper). The correction to the scatter- (see Amaro-Seoane et al. 2007; Gair et al. 2017). As MSPs are radio ing angle due to the spin of the secondary has been widely studied. sources, an “observable” electromagnetic counterpart will be present For example, Bini & Geralico(2017); Bini, Geralico & Vines(2017) if an EMRB/EMRI GW event occurs, allowing high-precision mea- calculated the analytical formula (to the first order in spin) of the surements to be made in the electromagnetic messenger domain and orbit of a spinning particle on the equatorial plane around a mas- in the gravitational messenger domain independently. In the theoret- sive BH using the MPD equations. Bini & Geralico(2018) studied ical perspective, the fact that neutron stars have a relatively narrow the effects of spin and spin-induced quadrupolar moments on the mass range (see Lattimer 2012; Özel & Freire 2016) reduces one scattering angle in the high-energy limit. MPD equations are also adopted in the theoretical modelling of scattering system by Vines, Kunst, Steinhoff & Hinderer(2016); Vines(2018); Vines, Steinhoff ★ E-mail: [email protected] (KJL), [email protected] (KW), pkle- & Buonanno(2019); Antonelli, Kavanagh, Khalil, Steinhoff & Vines [email protected] (PKL), [email protected] (DS) (2020), etc. While analytical analysis of the dynamics of spinning

The Pantheon+ Analysis: Evaluating Peculiar Velocity Corrections in Cosmological Analyses with Nearby Type Ia Supernovae

Erik R. Peterson,1 W. D’Arcy Kenworthy,2 Daniel Scolnic,1 Adam G. Riess,2, 3 Dillon Brout,4, 5 Anthony Carr,6 Hel´ ene` Courtois,7 Tamara Davis,6 Arianna Dwomoh,1 David O. Jones,8, 5 Brodie Popovic,1 Benjamin M. Rose,1 and Khaled Said6

1Department of Physics, Duke University, Durham, NC 27708, USA 2Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 3Space Telescope Science Institute, Baltimore, MD 21218, USA 4Center for Astrophysics, Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 5NASA Einstein Fellow 6School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 7University Lyon 1, IUF, IP2I Lyon, 69622 Villeurbanne cedex, France 8Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 92064, USA

ABSTRACT Separating the components of redshift due to expansion and motion in the nearby universe (z < 0.1) is critical for using Type Ia Supernovae (SNe Ia) to measure the Hubble constant (H0) and the equation-of-state parameter of dark energy (w). Here, we study the two dominant ‘motions’ contributing to nearby peculiar redshifts: large-scale, coherent-flow (CF) motions and small-scale motions due to gravitationally-associated galaxies deemed to be in a galaxy group. We use a set of 585 low-z SNe from the Pantheon+ sample, and evaluate the efficacy of corrections to these motions by measuring the improvement of SN distance residuals. We study multiple methods for modeling the large and small-scale motions and show that while group assignments and CF corrections individually contribute to small improvements in Hubble residual scatter, the greatest improvement comes from the combination of the two (relative standard deviation of the Hubble residuals RSD improves from 0.167 mag to 0.157 mag). We find the optimal flow corrections derived from various local density maps −1 −1 significantly reduce Hubble residuals while raising H0 by ∼ 0.4 km s Mpc as compared to using CMB redshifts, disfavoring the hypothesis that unrecognized local structure could resolve the Hubble tension. We estimate that the systematic uncertainties in cosmological parameters after optimally −1 −1 correcting redshifts are 0.08–0.17 km s Mpc in H0 and 0.02–0.03 in w which are smaller than the −1 −1 statistical uncertainties for these measurements: 1.5 km s Mpc for H0 and 0.04 for w.

Keywords: cosmology: peculiar velocities, coherent-ﬂows, galaxy groups

1. INTRODUCTION Recession velocities, vcosmo, are related to cosmologi- Type Ia Supernovae (SNe Ia) are critical tools for mea- cal redshifts zcosmo by Z zcosmo arXiv:2110.03487v1 [astro-ph.CO] 7 Oct 2021 suring cosmological parameters including the Hubble dz constant H , which parameterizes the local expansion vcosmo = c , (1) 0 0 E(z) rate of the universe (Riess et al. 2016; Freedman et al. 2019), and the equation-of-state of dark energy param- where we have used the dimensionless Hubble parameter eter w (JLA: Betoule et al. 2014; Pantheon: Scolnic E(z) = H(z)/H0 (Harrison 1993; Peebles 1993). zcosmo et al. 2018; DES: Brout et al. 2019). Measurements is not directly measured, but can be determined from the observed redshift corrected to the cosmic microwave of H0 and/or w which make use of SNe at low-redshift (0.01 < z < 0.1) are sensitive to the ‘peculiar velocities’ background (CMB) rest frame zCMB and the peculiar of host galaxies. Peculiar velocities (PVs) are the phys- redshift approximated as zpec ≈ vpec/c (an approxima- ical motion of a galaxy relative to the cosmological rest tion to the special relativistic Doppler shift) such that:

frame. 1 + zCMB 1 + zcosmo = . (2) 1 + zpec Draft version October 8, 2021 Typeset using LATEX default style in AASTeX631

The Pantheon+ Analysis: Dependence of Cosmological Constraints on Photometric-Zeropoint Uncertainties of Supernova Surveys

Sasha Brownsberger ,1 Dillon Brout ,1, 2 Daniel Scolnic ,3 Christopher W. Stubbs ,1, 2 and Adam G. Riess 4, 5

1Department of Physics, Harvard University 17 Oxford Street Cambridge, MA 02138, USA 2Center for Astrophysics, Harvard & Smithsonian 60 Garden Street Cambridge, MA 02138, USA 3Department of Physics, Duke University Science Drive Durham, NC 27710, USA 4Department of Physics and Astronomy, Johns Hopkins University 3400 Charles Street Baltimore, MD 21218, USA 5Space Telescope Science Institute 3700 San Martin Drive Baltimore, MD 21218, USA

ABSTRACT

Type Ia supernovae (SNe Ia) measurements of the Hubble constant, H0, the cosmological mass density, ΩM , and the dark energy equation-of-state parameter, w, rely on numerous SNe surveys using distinct photometric systems across three decades of observation. Here, we determine the sensitivities of the upcoming SH0ES+Pantheon+ constraints on H0,ΩM , and w to unknown systematics in the relative photometric zeropoint calibration between the 17 surveys that comprise the Pantheon+ supernovae data set. Varying the zeropoints of these surveys simultaneously with the cosmological parameters, we determine that the SH0ES+Pantheon+ measurement of H0 is robust against inter-survey photometric miscalibration, but that the measurements of ΩM and w are not. Specifically, we find that miscalibrated inter-survey systematics could represent a source of uncertainty in the measured value of H0 that is 1 1 no larger than 0.2 km s− Mpc− . This modest increase in H0 uncertainty could not account for 1 1 the 7 km s− Mpc− “Hubble Tension” between the SH0ES measurement of H0 and the Planck ΛCDM-based inference of H0. However, we find that the SH0ES+Pantheon+ best-fit values of ΩM and w respectively slip, to first order, by 0.04 and 0.17 per 25 mmag of inter-survey calibration − uncertainty, underscoring the vital role that cross-calibration plays in accurately measuring these parameters. Because the Pantheon+ compendium contains many surveys that share low-z Hubble Flow and Cepheid-paired SNe, the SH0ES+Pantheon+ joint constraint of H0 is robust against inter- arXiv:2110.03486v1 [astro-ph.CO] 7 Oct 2021 survey photometric calibration errors, and such errors do not represent an impediment to jointly using SH0ES+Pantheon+ to measure H0 to 1% accuracy.

Keywords: Type Ia supernovae(1728) — Cosmology(343) — Flux calibration(544) — Hubble constant(758) — Dark energy(351) Astronomy data analysis(1858)

1. INTRODUCTION

[email protected] Astronomy & Astrophysics manuscript no. fdd ©ESO 2021 October 8, 2021

Fourier-domain dedispersion C. G. Bassa, J. W. Romein, B. Veenboer, S. van der Vlugt, and S. J. Wijnholds

ASTRON, Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD, Dwingeloo, The Netherlands

Received October 8, 2021; Accepted October 8, 2021

ABSTRACT

We present and implement the concept of the Fourier-domain dedispersion (FDD) algorithm, a brute-force incoherent dedispersion algorithm. This algorithm corrects the frequency-dependent dispersion delays in the arrival time of radio emission from sources such as radio pulsars and fast radio bursts. Where traditional time-domain dedispersion algorithms correct time delays using time shifts, the FDD algorithm performs these shifts by applying phase rotations to the Fourier-transformed time-series data. Incoherent dedispersion to many trial dispersion measures (DMs) is compute, memory-bandwidth and I/O intensive and dedispersion algorithms have been implemented on Graphics Processing Units (GPUs) to achieve high computational performance. However, time-domain dedispersion algorithms have low arithmetic intensity and are therefore often memory-bandwidth limited. The FDD algorithm avoids this limitation and is compute limited, providing a path to exploit the potential of current and upcoming generations of GPUs. We implement the FDD algorithm as an extension of the dedisp time-domain dedispersion software. We compare the performance and energy-to-completion of the FDD implementation using an NVIDIA Titan RTX GPU against the standard as well as an optimized version of dedisp. The optimized implementation already provides a factor of 1.5 to 2 speedup at only 66% of the energy utilization compared to the original algorithm. We ﬁnd that the FDD algorithm outperforms the optimized time-domain dedispersion algorithm by another 20% in performance and 5% in energy-to-completion when a large number of DMs (& 512) are required. The FDD algorithm provides additional performance improvements for FFT-based periodicity surveys of radio pulsars, as the Fourier transform back to the time domain can be omitted. We expect that this computational performance gain will further improve in the future since the Fourier-domain dedispersion algorithm better matches the trends in technological advancements of GPU development. Key words. methods: data analysis – pulsars: general

1. Introduction In surveys for new radio pulsars or fast radio bursts, the dispersion measure of a new source is a priori unknown and the Radio waves propagating through the ionized intergalactic, inter- input data need to be dedispersed to many, usually thousands, stellar and/or interplanetary medium undergo dispersion, which of trial dispersion measures. The majority of current surveys introduces a frequency-dependent time delay between emission perform incoherent dedispersion on Stokes I time-series, typ- at the source and reception at the telescope (e.g. Lorimer & ically with sample times between 50 to 100 µs, and channel Kramer 2012). Hence, pulsed radio signals from astrophysical sizes of 0.1 to 0.5 MHz (Keith et al. 2010; Stovall et al. 2014; sources like radio pulsars and fast radio bursts (see Kaspi & Lazarus et al. 2015; Bhattacharyya et al. 2016, see also Lyon Kramer 2016; Petroff et al. 2019, for reviews) have a distinctive et al. 2016), as this requires lower data rates and less storage dispersion measure (DM), directly related to the electron column compared to voltage data. Only at the lowest observing frequen- density between the source and the observer. If not properly cor- cies (ν . 300 MHz), voltage data are recorded to perform a com- rected for, the dispersion delays lead to smearing and hence loss bination of coherent and incoherent dedispersion (Bassa et al. of signal-to-noise of the received radio pulses. 2017a,b; Pleunis et al. 2017). Correcting for the dispersion is called dedispersion, and can The computational cost to incoherently dedisperse is signif- 3 be performed using two techniques, depending on the state of the icant, as the input data typically consist of Nν ∼ 10 frequency 6 radio observation. For Nyquist-sampled time-series data of or- channels for Nt ∼ 10 time samples, and needs to be dedispersed 4 thogonal polarizations (voltage data hereafter), the phase infor- to NDM ∼ 10 dispersion measures. Hence, algorithms have been mation can be used to coherently dedisperse the signals through developed to either reduce the computational cost by reducing arXiv:2110.03482v1 [astro-ph.IM] 7 Oct 2021 a convolution with the inverse of the response function of dis- the complexity, by using e.g. a tree algorithm (Taylor 1974) or persion (Hankins 1971; Hankins & Rickett 1975). This tech- a piece-wise linear approximation of the dispersion delay (Ran- nique completely removes the effects of dispersion. With coher- som 2001), a transformation algorithm (Zackay & Ofek 2017) ent dedispersion, dedispersion is performed before squaring the and/or by using hardware acceleration of the direct, brute force, orthogonal polarizations signals to form the Stokes parameters. dedispersion approach. A review of these algorithms, and an im- For signals that have already been squared to form the Stokes pa- plementation of the direct dedispersion algorithm on graphics rameters, the incoherent dedispersion method is used. The time- processing units (GPUs), is presented by Barsdell et al.(2012). series data of finite channels are shifted in time to correct for Sclocco et al.(2016) found that the direct, brute force, the dispersion delay between channels, while dispersive smear- dedispersion algorithm on many-core accelerators is inherently ing within channels remains uncorrected. Hence, the size of fre- memory-bandwidth bound, due to the low arithmetic intensity quency channels and time samples is typically optimized to min- of the required calculations. Hence, the direct dedispersion algo- imize dispersive smearing within channels. rithm does not optimally benefit from advances in the computing

Article number, page 1 of 7 Acceleration of ultra-high energy cosmic rays in the early afterglows of gamma-ray bursts: concurrence of jet’s dynamics and wave-particle interactions

Ze-Lin Zhang,1, 2 Ruo-Yu Liu,1, 2, ∗ and Xiang-Yu Wang1, 2, † 1School of Astronomy and Space Science, Nanjing University, Xianlin Road 163, Nanjing 210023, China 2Key laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, China The origin of ultra-high energy cosmic rays (UHECRs) remains a mystery. It has been suggested that UHE- CRs can be produced by the stochastic acceleration in relativistic jets of gamma-ray bursts (GRBs) at the early afterglow phase. Here, we develop a time-dependent model for proton energization by cascading compressible waves in GRB jets considering the concurrent effect of the jet’s dynamics and the mutual interactions between turbulent waves and particles. Considering fast mode of magnetosonic wave as the dominant particle scatterer and assuming interstellar medium (ISM) for the circumburst environment, our numerical results suggest that protons can be accelerated up to 1019 eV during the early afterglow. An estimation shows ultra-high energy nuclei can easily survive photodisintegration in the external shocks in most cases, thus allowing the acceleration of 1020 eV cosmic-ray nuclei in the proposed frame. The spectral slope can be as hard as dN/dE ∝ E 0, which is consistent with the requirement for the interpretation of intermediate-mass composition of UHECR as measured by the Pierre Auger Observatory.

PACS numbers: 45.50.Dd, 52.35.Ra, 94.05.Pt, 96.50.sb, 98.70.Rz

I. INTRODUCTION is 1044 erg Mpc−3yr−1 [29–31] while the gamma-ray energy production rate of GRBs is 1043 erg Mpc−3yr−1 for a typical 52 −1 Ultra-high energy cosmic rays (UHECRs) at the ankle en- gamma-ray luminosity of 10 erg s and a local GRB rate −3 −1 ergy 1018.5 eV and above are the most energetic particles in of 1 Gpc yr . Given the predicted spectral slope of the ac- nature. The presence of these particles has been known for celerated particles being p & 2 for the relativistic shock acceleration [32–35], the fraction of the energy of CRs accelerated over half a century [1]. However, the sites and mechanisms 18.5 of their production are still open questions [2]. The study of beyond the ankle (10 eV) is only at the level of 10% of the energy spectrum and the mass composition of UHECR the total CR energy. As a result, it would require a baryon helps to reveal their origin. Recently, the results of cosmic- loading factor (defined as the ratio of total energy populated ray anisotropy observed by Pierre Auger Observatory and the in CRs to that in gamma rays) of ∼ 100 to account for the Telescope Array support the hypothesis of an extragalactic required UHECR energy production rate. It is in tension with origin for the UHECR [3,4]. Extragalactic sources, such as the constraint from the non-detection of GRB neutrinos by the active galactic nuclei [5,6], Gamma-ray bursts (GRBs) [7– IceCube neutrino telescopes in some dissipation mechanisms 9], energetic supernovae (such as hypernovae) [10, 11], tidal of GRBs [36]. Furthermore, it has been pointed out that a very disruption events [12–14], galaxy clusters [15–17], as well hard CR injection spectrum with p . 1 is favored in order to as milli-second magnetars [18, 19], have been considered as fit the spectrum and composition of UHECRs measured by plausible candidates of UHECR sources. the Pierre Auger Observatory, where the best-fit index is even As the most powerful and intense explosive events in the p < 0 [37, 38]. universe, GRBs have been studied extensively as the cos- Recently, a stochastic acceleration (SA) model of UHE- mic accelerator of UHECRs [7–9, 20–23]. However, the CRs via turbulence in GRB jets has been proposed to avoid acceleration mechanisms of these particles in GRBs remain the problems mentioned above [22]. The SA can yield a an enigma. The standard scenario adopted to produce non- hard UHECR spectrum with shallow index p . 2, which thermal particles is the particle acceleration at shocks, e.g., has been discussed as a possible charged particle accelera- diffusive shock acceleration [24, 25]. However, particle accel- tion mechanism in astrophysical plasmas [39–41]. Magneto- eration by relativistic shocks with bulk Lorentz factor Γ 1 hydrodynamic (MHD) turbulence is indispensable in various arXiv:2110.03458v1 [astro-ph.HE] 7 Oct 2021 is limited by a series of factors. For example, the relative en- astrophysical processes. As the magnetic scattering centers in ergy gain drops quickly ( from Γ2 to ' 2 ) after the first shock SA scenario, MHD waves mainly consist of three types: in- crossing circle because of particles do not have sufficient time compressible Alfven´ modes and compressible fast and slow to become isotropic upstream before being caught up by the modes [42]. Particle scattering and diffusion largely rely on shock [e.g., 26, 27, see also a recent review by Marcowith et the properties of plasma turbulence. Fast mode waves show al. 28]. Another possible disadvantage of the shock acceler- an isotropic cascade and it could be the most effective scat- ation is the energy budget. The required energy production terers of cosmic-rays [43, 44]. The spectrum of the isotropic −3/2 rate of CRs to explain the measured flux beyond the ankle cascade was claimed to be k [42]. The excitation of turbulence in plasmas generally stems from the anisotropy of particle distributions and MHD instabilities [45]. In our work, we consider the turbulence ∗ [email protected] driven by MHD instabilities induced by the jet’s propagation † [email protected] in the circumburst interstellar medium (ISM), such as Kelvin- Received 6 October, 2021; Revised 6 October, 2021; Accepted 7 October, 2021

DOI: xxx/xxxx

ARTICLE TYPE The gamma-ray emission from young radio galaxies and quasars

1Universitá di Trieste, Dipartimento di Fisica, I-34127 Trieste, Italy Abstract 2Istituto Nazionale di Fisica Nucleare, According to radiative models, radio galaxies are predicted to produce -rays from Sezione di Trieste, I-34127 Trieste, Italy 3Istituto Nazionale di Astrofisica - Istituto di the earliest stages of their evolution onwards. The study of the high-energy emission Radioastronomia, I-40129 Bologna, Italy from young radio sources is crucial for providing information on the most energetic 4Dipartimento di Fisica “M. Merlin” dell’Università e del Politecnico di Bari, processes associated with these sources, the actual region responsible for this emis- I-70126 Bari, Italy sion, as well as the structure of the newly born radio jets. Despite systematic searches 5 Istituto Nazionale di Fisica Nucleare, for young radio sources at -ray energies, only a handful of detections have been Sezione di Bari, I-70126 Bari, Italy reported so far. Taking advantage of more than 11 years of Fermi-LAT data, we inves- Correspondence tigate the -ray emission of 162 young radio sources (103 galaxies and 59 quasars), *Giacomo Principe, Email: [email protected] the largest sample of young radio sources used so far for a -ray study. We analysed the Fermi-LAT data of each source separately to search for a significant detection. In addition, we performed the first stacking analysis of this class of sources in order to investigate the -ray emission of the young radio sources that are undetected at high energies. In this note we present the results of our study and we discuss their implications for the predictions of -ray emission from this class of sources.

KEYWORDS: Galaxies: evolution – galaxies: active – galaxies: jets – radio continuum: galaxies – gamma-rays: galaxies

1 INTRODUCTION determination of kinematic and radiative ages in some of the most compact sources (t ∼ 102 − 105 years) (Phillips & One of the greatest questions investigated by modern astro- Mutel, 1982). Compact radio objects reside either in galaxies physics is understanding the origin of the -ray emission in or quasars. While for the latter the -ray emission is favored radio galaxies and quasars. While the extragalactic -ray sky by their smaller jet-inclination angle and beaming effects, the is dominated by blazars, due to their small jet inclination to origin of -ray emission in galaxies is still a matter of debate. the line of sight and by relativistic beaming, only a few radio In the model proposed by Stawarz et al. (2008), young radio arXiv:2110.03416v1 [astro-ph.HE] 7 Oct 2021 galaxies have been detected so far (4LAC, Ajello et al., 2020). galaxies are expected to produce isotropic high-energy emis- Due to their misaligned jets, they offer a unique tool to investi- sion through IC scattering of the UV photons by the electrons gate some of the non-thermal processes at work in unbeamed in the compact radio lobes. Depending on the sources’ phys- regions in active galactic nuclei (AGN). ical parameters (e.g. jet power, UV luminosity), the model According to the evolutionary scenario, the size of a radio predicted that young radio sources could constitute a class of galaxy is strictly related to its age (Fanti et al., 1995). There- -ray emitters detectable by Fermi-LAT. However, systematic fore extragalactic compact radio objects (i.e. radio objects with searches for rays from young radio sources have so far been projected linear size LS<20 kpc), are expected to be the pro- unsuccessful (D’Ammando, Orienti, Giroletti, & Fermi LAT genitors of extended radio galaxies (Readhead et al., 1996). coll., 2016). Dedicated studies reported a handful of detec- Support to the young nature of these objects is given by the tions: three young radio galaxies (NGC 6328 (Migliori et al., 2016), NGC 3894 (Principe et al., 2020) and TXS 0128+554 MNRAS 000,1–9 (2020) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

Radiative energy from a reconnection region around massive black hole

Tian-Le Zhao,1,2★ and Rajiv Kumar1,2 1CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei 230026, China 2School of Astronomy and Space Sciences, University of Science and Technology of China, Hefei 230026, China

Last updated XX; in original form XX

ABSTRACT In the previous numerical study, we find the blob formation and ejection in the presence of magnetic reconnection in the environment of the hot flow of the accretion disk. Based on those encouraging results, in the present work, we calculate the energy and the spectrum of the emission in the different bands around Sgr A*. We assume the electrons in the magnetic reconnection region are non-thermal and emits synchrotron radiation. The electrons in the other region are thermal, which follows the thermal distribution, and the thermal electron emission mechanism is thermal synchrotron radiation. During the whole process of the magnetic evolution and reconnection, we find two peaks in the temporal light curve in the recently observed Radio frequencies (230 GHz and 43 GHz) and NIR wavelengths (3.8휇푚 and 2.2휇푚). Although, the light curve of the NIR band is most prominent in a single peak. The first peak appears because of the blob in the plasma flow, which is formed due to the magnetic reconnection. The second peak appears due to the production of the non-thermal electrons with the evolution of the magnetic flux. Both peaks reach luminosity more than 1026 erg/s for a single plasmoid/blob. For the NIR band, the highest luminosity can reach more than 1028 erg/s.These luminosities can be high for the large simulation area and the stronger magnetic field with the multiple blobs. We infer that the observed flares are a group of magnetic reconnection phenomena, not a single one. Key words: Galaxies: jets – Physical Data and Processes: MHD, magnetic reconnection, plasmas, radiation mechanisms: synchrotron radiation

CONTENTS

1 Introduction 2 Calculation of radiative emission 2.1 A brief introduction to MHD numerical simulation 2.2 Non thermal Synchrotron radiation 2.3 Thermal synchrotron radiation 3 Results 4 Conclusion arXiv:2110.03409v1 [astro-ph.HE] 7 Oct 2021

★ E-mail:[email protected]

Constraining the properties of the potential embedded planets in the disk around HD 100546 Max Ackermann Pyerin1, Timmy N. Delage1, Nicolás T. Kurtovic1, Matías Gárate1, Thomas Henning1, and Paola Pinilla1, 2

1 Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany 2 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK

Bachelor Project at the University of Heidelberg. Lead Author e-mail: [email protected], Group Leader e-mail: [email protected]

Received _; accepted _

ABSTRACT

Context. The protoplanetary disk around the star HD 100546 displays prominent substructures in the form of two concentric rings. Recent observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) have revealed these features with high angular resolution, and well resolved the faint outer ring. This allows us to study the nature of the system further. Aims. Our aim is to constrain some of the properties of potential planets embedded in the disk, assuming that they are responsible for inducing the observed rings and gaps. Methods. We present the self-calibrated 0.9 mm ALMA observations of the dust continuum emission from the circumstellar disk around HD 100546. These observations reveal substructures in the disk which are consistent with two rings, the outer ring being much fainter than the inner one. We reproduce this appearance closely with a numerical model that assumes two embedded planets. We vary planet and disk parameters in the framework of the planet-disk interaction code FARGO3D, and use the outputs for the gas and dust distribution to generate synthetic observations with the code RADMC-3D. Results. From this comparison, we find that an inner planet located at r1 = 13 au with a mass M1 = 8MJup, and an outer one located at r2 = 143 au with a mass M2 = 3MJup leads to the best agreement between synthetic and ALMA observations (deviation less than 3σ for the normalized radial profiles). To match the very low brightness of the outer structure relative to the inner ring, the initial disk gas surface density profile needs to follow an exponentially tapered power-law (self-similar solution), rather than a simple power-law profile. Key words. protoplanetary disks – methods: numerical – submillimeter: planetary systems – planets and satellites: formation – planet-disk interactions

1. Introduction emission and 230 ± 20 au in continuum emission (Walsh et al. 2014), and have revealed the following substructures: a ring of How planets form in protoplanetary disks is one of the main sub- emission with a radius of 27 au (Hubble Telescope, Ardila et al. jects of study in modern astrophysics. A very intriguing part of 2007), and a second very faint ring at 200 au (Atacama Large this topic is how newborn planets interact and actively shape the Millimeter/sub-millimeter Array - ALMA, Walsh et al. 2014). disk they are embedded in. This idea of planet-induced substruc- While several other mechanisms have been proposed to pro- tures is currently driven by recent observations with state-of- duce ring-like substructures in protoplanetary disks, such as the the-art telescopes, which have revealed that circumstellar disks dead zone outer edge mechanism (Pinilla et al. 2016), parti- can exhibit various prominent substructures, such as rings, gaps, cle growth by condensation near ice lines (Zhang et al. 2016; asymmetries or spirals (e.g., ALMA Partnership et al. 2015; Stammler et al. 2017) and magnetic disk winds (Suriano et al. Long et al. 2018; Andrews et al. 2018; Andrews 2020; Cieza 2019), in this work we assume the planet-disk interaction to be et al. 2021). Some of these substructures are believed to be in- arXiv:2110.03373v1 [astro-ph.EP] 7 Oct 2021 the main driver for the observed substructures of HD 100546. deed created by planet-disk interactions, rendering the observed In this mechanism, an embedded planet creates a pressure trap features as indicators of planet formation (e.g., Pinilla et al. outwards of its orbit which triggers an accumulation of dust by 2012, 2015; Zhu & Stone 2014), and could be used to constrain preventing the material from moving inwards. This has been mo- the properties of an embedded planet candidate, such as the mass tivated in previous works, where simulations were executed to and the orbital separation (e.g., Rice et al. 2006). prove the similarity between structures created by planets and The subject of this study is the circumstellar disk surround- the observed structures of HD 100546 (Walsh et al. 2014; Mon- ing the Herbig star HD 100546. This target is found at a distance tesinos et al. 2015; Quanz et al. 2015; Pinilla et al. 2015; Fedele of 108.1 ± 0.5 pc (Brown et al. 2021), with a right ascension of et al. 2021). For example, Pinilla et al. 2015 estimated that a 11h 33m 25.3s and a declination of −70◦ 110 41.200 (Brown et al. model with two high-mass planets can qualitatively reproduce 2021). Observations at multiple diﬀerent wavelengths detected the two rings observed in this disk, with the inner planet be- this disk spanning 390 ± 20 au outwards from the star in CO ing located at 10 au and the outer one at 70 au. Furthermore,

Article number, page 1 of 12 MNRAS 000,1–13 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style ﬁle v3.0

Understanding the X-ray spectral curvature of Mkn 421 using broadband AstroSat observations

Jyotishree Hota,1? Zahir Shah,2,3† Rukaiya Khatoon,4 Ranjeev Misra,2 Ananta C. Pradhan,1 and Rupjyoti Gogoi,4 1 Department of Physics and Astronomy, National Institute of Technology Rourkela, Odisha ,India 2 Inter-University Center for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune - 411007, India 3 Department of Physics, Central University of Kashmir, India 4 Tezpur University,Napaam-784028, Assam, India.

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We present a time-resolved X-ray spectral study of the high energy peaked blazar Mkn 421 using simultaneous broadband observations from the LAXPC and SXT instruments on-board AstroSat. The ∼ 400 ksec long observation taken during 3–8 Jan- uary, 2017 was divided into segments of 10 ksecs. Each segment was fitted using synchrotron emission from particles whose energy distribution was represented by a log-parabola model. We also considered particle energy distribution models where (i) the radiative cooling leads to a maximum energy (ξmax model), (ii) the system has energy dependent diffusion (EDD) and (iii) has energy dependent acceleration (EDA). We found that all these models describe the spectra, although the EDD and EDA models were marginally better. Time resolved spectral analysis allowed for studying the correlation between the spectral parameters for different models. In the simplest and direct approach, the observed correlations are not compatible with the predictions of the ξmax model. While the EDD and EDA models do predict the correlations, the values of the inferred physical parameters are not compatible with the model assumptions. Thus, we show that spectrally degenerate models, can be distinguished based on spectral parameter correlations (especially those between the model normalization and spectral shape ones) making time-resolved spectroscopy a powerful tool to probe the nature of these systems. Key words: galaxies: active – BL Lacertae objects: general – BL Lacertae objects: individual: Mkn 421 – acceleration of particles – diffusion – X-rays: galaxies

1 INTRODUCTION are two types of emission processes, that have been used to explain the broadband SED viz. synchrotron process and Blazars belong to a subclass of radio loud active galactic the inverse-Compton (IC) process. The low energy hump nuclei (AGN) with a relativistic jet directed towards the is produced by relativistic electrons interacting with ambi- observer (Urry & Padovani 1995; Blandford et al. 2019). arXiv:2110.03344v1 [astro-ph.GA] 7 Oct 2021 ent jet magnetic ﬁeld through synchrotron process, where They are characterized by non-thermal broadband spectra, as the high energy hump is mainly explained by IC scatter- rapid variability, strong radio and optical polarization (Sam- ing of low-energy photons (Urry & Mushotzky 1982; Ghis- bruna 2000; Fan et al. 2008). The spectral energy distribu- ellini et al. 1985; Begelman et al. 1987; Blandford & Levin- tion (SED) of a blazar extends from radio to γ-ray ener- son 1995; Bloom & Marscher 1996; Sokolov et al. 2004). gies (Fossati et al. 1998) and is generally characterized by The seed photons for IC scattering in leptonic model can a double hump structure with the low energy hump peak- be either synchrotron photons (synchrotron self-Compton, ing in the optical to X-ray energies and the high energy (SSC): Jones et al.(1974); Maraschi et al.(1992); Ghisellini hump with peak located in γ-ray band. The high-energy et al.(1993)) or photons external to jet (external Comp- emission can be produced either through leptonic models or ton (EC): Dermer et al.(1992); Sikora et al.(1994); B la˙ze- through hadronic models. Based on leptonic model, there jowski et al.(2000); Shah et al.(2017)). Alternatively, in the hadronic processes, the the higher frequency emission ? E-mail: [email protected] can also be initiated by relativistic protons such as the pro- † [email protected]

© 2021 The Authors arXiv:2110.03340v1 [astro-ph.EP] 7 Oct 2021 ieae l 2010 al. et Rivera 2021 al. et Agol 2018 al. et Delrez htrsnn hissc steoese nTRAPPIST-1 in seen one shown the have as studies such several chains accretion, resonant planetesimal that or tion satelli Jovian the Ganymede. of and resonances configuration Europa, Laplace orbital as the to to referred reference often are resonances body frsnne novn he rmr lnt se e.g., (see, chains planets more in or observed disc three been gaseous involving have a resonances systems in of exoplanetary planets two eral of attribute or- often migration ( common is convergent MMRs) most (hereafter to the resonances not motion ( mean are configuration they bital although systems, etary 2021 2016 ( 5:3, 3:2 8:5, and pair): 4:3, innermost t the to with known close (starting chain planets ratios lowing adjacent resonant between longest ratios the period with in date, possibly are planets the at-iepaesaon nutaco e wr tro ma of star seven dwarf red least ultra-cool at an of around consists planets system Earth-size planetary TRAPPIST-1 The Introduction 1. M coe ,2021 8, October Astronomy nlgoee l 2001 al. et Snellgrove ⋆ ∗ RPIT1 yaia nlsso h rni-iigvari transit-timing the of analysis Dynamical TRAPPIST-1: -al [email protected] E-mail: sn omto oesbsdete npbl accre- pebble on either based models formation Using eoatcngrtoshv enosre nohrplan- other in observed been have configurations Resonant n h frmnindTAPS- ytm.Three- system). TRAPPIST-1 aforementioned the and , ; ≈ ih2pae enmto eoacs hl ogrperiod longer while resonances, Shorte system. mean-motion the 2-planet of with variables dynamical the of of series a variations possess signals (TTV) variation transit-timing eoacs–smlrt h bevdoe–cnb omdwhich formed be can – one observed chai the resonant to the similar of – origin the resonances explore we forces, migration ob narsnn hi rudarddaf eso htalth all that show We resonan 2-planet dwarf. the red which for a pair around innermost chain the from resonant TR apart a the of in solution be best-fit to published recently the analyze We ZZZ form original in YYY; Received XXX. Accepted htwiete4y olce aao bevtoshl e i key hold observations of data collected 4-yr the while that e words. fo Key on constraints additional therefore and planets, of give pair also would observations resonanc Additional 3-body solution. best-fit of signatures provide soon will TRAPPIST-1 aelts yaia vlto n tblt lnt an planets – stability and evolution dynamical satellites: 2 1 il ta.2016 al. et Mills srnm eatetadVrulPaeayLbrtr,Un Laboratory, Planetary M Virtual of and Department Department Astronomy Systems, Complex Belgium for Namur, Institute Namur naXys, 0 . 09 & M ue ta.2017 al. et Luger ⊙ Astrophysics .Oermral etr fTAPS- sthat is TRAPPIST-1 of feature remarkable One ). ( ; ; ilne l 2016 al. et Gillon eeta ehnc lntds neatos–protopla – interactions planet-disc – mechanics celestial a roe ta.2018 al. et Grootel Van oz zesie l 2016 al. Go´zdziewski et arcye l 2014 al. et Fabrycky ; ; e el 2002 Peale & Lee hitasne l 2018 al. et Christiansen aucitn.trappist no. manuscript ). n rgno h eoatchain resonant the of origin and , 2017 .Teyssandier J. .Tefrainof formation The ). .Frhroe sev- Furthermore, ). ; ; ; urte l 2018 al. et Ducrot ue ta.2017 al. et Luger aDnl tal. et MacDonald ; ee tal. et Leleu :,3:2, 3:2, aelts detection satellites: d 1, mto scenarios. rmation omnprosvrigfo ast eae,wihaeas p also are which decades, to days from varying periods common oepeieifraino h euirrsnn dynamics resonant peculiar the on information precise more efol- he e Io, tes ⋆ eid eg,teTV hrceitctmsaeo . yr 1.3 of timescale characteristic TTVs the (e.g., periods r nlscruae ymaso rqec nlss ehigh we analysis, frequency a of means By circulate. angles t ABSTRACT fTAPS- n n htfrpriua icconditions, disc particular for that find and TRAPPIST-1 of n -.Libert A-S. , PIT1sse,wihcnit fsvnErhsz planet Earth-size seven of consists which system, APPIST-1 in , rs rm3pae eoacs yueof use By resonances. 3-planet from arise s s npriua h . n . rproiiisepce f expected periodicities yr 5.1 and 3.3 the particular in es, fraino h -lntrsnn yais ute moni further dynamics, resonant 2-planet the on nformation to ss d vriyo ahntn etl,W 89,USA 98195, WA Seattle, Washington, of iversity ; ; lnt r iutnosyi -lntad3pae reson 3-planet and 2-planet in simultaneously are planets e cuaeyrpoue h bevdTV.Oraayi sugg analysis Our TTVs. observed the reproduces accurately urn oiin,adweemgainfre r nyappl only planet, are t outermost forces outside the migration just to where placed and initially positions, are to current and assumed are formed, planets already where have approach simpler a Using MMR). a aual rs uigteerysae ffrainand formation ( of disc stages early a the in during migration arise naturally can ne h rsneo nenpaes( 2011 t al. planets et order pe Cochran unseen in strictly 2004 useful of the Schneider are presence to variations the respect These with infer case). times Keplerian transit odic the of viations b driven planets star. inner central the the lat of with re-assembling interactions migration tidal groups, outward to distinct thanks two on, in p the migrated that have arguing system, may the of evolution tidal and gration Finally Myr. 0.5 also of timescale early a o the on stable to is opposed by configuration as tions into timescales, a planets long Such TRAPPIST-1 very configuration. the of current capture their the allows migration al. et Coleman 2019 atfrti ok,btn ared pi hr-re reso- third-order in of inner up ratio the period since ends TRAPPIST-1 a pair for has rele- relevant pair no be is but may which (which work), in MMR nance few this 5:3 a the for with as MMRs), vant 4:3 (such and resonance 3:2 second-order 2:1, the (mainly resonances teais nvriyo au,6 u eBuels 5000 Bruxelles, de Rue 61 Namur, of University athematics, eaydss–paesadstlie:frain–paesa planets – formation satellites: and planets – discs netary h RPIT1sse rsnssrn Ts(.. de- (i.e., TTVs strong presents system TRAPPIST-1 The rse ta.2019 al. et Brasser ; oea ta.2019 al. et Coleman 1 n .Agol E. and ilne al. et Gillon ( 2019 ; gle l 2005 al. et Agol ; ode l 2012 al. et Ford first-order in up end pairs of majority the ), ( 2 2017 ; re ta.2017 al. et Ormel un re 2021 Ormel & Huang aaoe al. et Tamayo ∼ ; hc yial hwdinstability showed typically which ) un re 2021 Ormel & Huang 1 . ,adteeoele erte8:5 the near lies therefore and 6, N ; aaozue al. et Papaloizou ; bd iuain with simulations -body omn&Mra 2005 Murray & Holman ril ubr ae1o 11 of 1 page number, Article Ste ( 2017 iad-sué2002 Miralda-Escudé ff ; ne l 2012 al. et en fteinnermost the of connege al. et Schoonenberg r associated are ) hwdta disc that showed ) rtecurrent the or ih htthe that light tde h mi- the studied ) eeti the in resent N appearing s ations hi of chain a bd simula- -body .A oe by noted As ). oigof toring © ances, S 2021 ESO ( 2018 ests nd n to and ) lanets see , heir ied ri- er o y n ; ; Future Constraints on Primordial Black Holes from XGIS-THESEUS

Diptimoy Ghosh,1, ∗ Divya Sachdeva,1, † and Praniti Singh1, ‡ 1Department of Physics, Indian Institute of Science Education and Research Pune, India Current observations allow Primordial Black Holes (PBHs) in asteroid mass range 1017 − 1023 g to constitute the entire dark matter (DM) energy density (barring a small mass range constrained by 21 cm observations). In this work, we explore the possibility of probing PBH with masses 1017−1019 g via upcoming X and Gamma Imaging Spectrometer (XGIS) telescope array on-board the Transient High-Energy Sky and Early Universe Surveyor (THESEUS) mission. While our projected 16 17 limits are comparable with those proposed in the literature for 10 g < MPBH < 10 g, we show that the XGIS-THESEUS mission can potentially provide the strongest bound for 1017 g < 18 MPBH . 3 × 10 g for non-rotating PBHs. The bounds become more stringent by nearly an order 15 19 of magnitude for maximally rotating PBHs in the mass range 5 × 10 g < MPBH . 10 g.

I. INTRODUCTION

There are a large variety of astrophysical as well as cosmological observations in support of the Dark Matter (DM) hypothesis. In the absence of any DM candidate within the SM, a variety of models have been proposed to explain this additional matter content, with masses as small as 10−22 eV (lower bound comes from studies of dwarf galaxy morphology [1]). On the high mass end, PBHs are one of the popular candidates of DM that could have formed as a result of large density perturbations in the early universe [2–5]. The recent detection of BH merger through direct detection of gravitational waves by LIGO-Virgo [6,7] have given a new direction to search for PBH as DM [8–10]. FIG. 1. BH temperature T as function of its mass M for Regardless of how PBHs are produced, they produce BH BH non-rotating (a∗ = 0) and maximally rotating (a∗ = 0.9999) radiation through Hawking evaporation. The energy case. spectrum of this radiation is almost thermal. An un- charged rotating BH of mass MBH and angular momentum JBH radiates with temperature [11–14]: Consequently, PBHs with mass greater than 1015 g have lifetime longer than age of the Universe and can contribute to the DM density. p 2 1 1 − a∗ TBH = (1) p 2 4πGN MBH 1 + 1 − a∗ Based on the Hawking radiation, there exist cornu- copia of constraints from non-observation of emission 14 17 GN denotes Newton’s Gravitational constant and a∗ = products from PBH of mass region ∼ 5 × 10 − 2 × 10 2 JBH/GN MBH is the dimensionless spin parameter of g which are presently evaporating. Previous works have BH. The energy spectrum of the emitted particles peaks constrained this mass range via absence of any excess in around TBH (E ∼ 5.77 TBH for photons [14, 15]) and then electron-positron pairs [16–21], non-observation of neu- becomes exponentially suppressed for energies greater arXiv:2110.03333v1 [astro-ph.CO] 7 Oct 2021 trino flux from currents detectors [17] as well as the than TBH and falls off as a power law for E TBH. upcoming ones [22–24], the γ-rays observations [16, 25– Such evaporation causes mass loss on timescales given 28], X-rays observations [25, 29], observation of radio- by: signals [30, 31], 21-cm observations [32, 33] and heating of interstellar medium [34, 35] (for in-depth discussion G2M 3 M 3 on similar such studies of PBH, refer to Ref. [36–38]). τ(M) ∼ ∼ π × 1017 BH s (2) ~c4 1015 g. Due to strong constraints on DM fraction, fPBH, going as small as 10−8, this mass range cannot contribute to the total DM density. Nonetheless, there exist no constraints in the asteroid ∗ [email protected] mass range ∼ 2 × 1017 − 1023 g of PBHs as they radiate † [email protected] photons and neutrinos below the threshold sensitivity ‡ [email protected] of the present telescopes and detectors. Therefore, this MNRAS 000,1–21 (2015) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

PLATO Hare-and-Hounds exercise: Asteroseismic model ﬁtting of main-sequence solar-like pulsators

M. S. Cunha,1★ I. W. Roxburgh,2,3 V. Aguirre Børsen-Koch,4 W. H. Ball,3,4 S. Basu,5 W. J. Chaplin,3,4 M.-J. Goupil,6 B. Nsamba,7,1,8 J. Ong,5 D. R. Reese,6 K. Verma,4 K. Belkacem,6 T. Campante,1,9 J. Christensen-Dalsgaard,4 M. T. Clara,1,9 S. Deheuvels,10 M. J. P. F. G. Monteiro,1,9 A. Noll,10 R. M. Ouazzani,6 J. L. Rørsted,4 A. Stokholm,4,11 M. L. Winther4 1 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, PT4150-762 Porto, Portugal 2 Astronomy Unit, Queen Mary University of London, Mile End Road, London E1 4NS, UK 3 School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK 4 Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 5 Department of Astronomy, Yale University, PO Box 208101, New Haven, CT 06520-8101, USA 6 LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 Place Jules Janssen, 92195, Meudon, France 7 Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, D-85748 Garching, Germany 8 Kyambogo University, P.O. Box 1, Kyambogo, Kampala - Uganda 9 Departamento de Física e Astronomia, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre, s/n, PT-4169-007 Porto, Portugal 10 IRAP, Université de Toulouse, CNRS, CNES, UPS, Toulouse, France 11 Dipartimento di Fisica e Astronomia, Università degli Studi di Bologna, Via Gobetti 93/2, I-40129 Bologna, Italy

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Asteroseismology is a powerful tool to infer fundamental stellar properties. The use of these asteroseismic-inferred properties in a growing number of astrophysical contexts makes it vital to understand their accuracy. Consequently, we performed a hare-and-hounds exercise where the hares simulated data for 6 artificial main-sequence stars and the hounds inferred their properties based on different inference procedures. To mimic a pipeline such as that planned for the PLATO mission, all hounds used the same model grid. Some stars were simulated using the physics adopted in the grid, others a different one. The maximum relative differences found (in absolute value) between the inferred and true values of the mass, radius, and age were 4.32 per cent, 1.33 per cent, and 11.25 per cent, respectively. The largest systematic differences in radius and age were found for a star simulated assuming gravitational settling, not accounted for in the model grid, with biases of -0.88 per cent (radius) and 8.66 per cent (age). For the mass, the most significant bias (-3.16 per cent) was found for a star with a helium enrichment ratio outside the grid range. Moreover, a ∼7 per cent dispersion in age was found when adopting different prescriptions for the surface corrections or shifting the classical observations by ±1휎. The choice of the relative weight given to the classical and seismic constraints also impacted significantly the accuracy and precision of the results. Interestingly, arXiv:2110.03332v1 [astro-ph.SR] 7 Oct 2021 only a few frequencies were required to achieve accurate results on the mass and radius. For the age the same was true when at least one 푙 = 2 mode was considered. Key words: asteroseismology – stars: fundamental parameters – stars: evolution – stars: oscillations – methods: statistical

1 INTRODUCTION the general astrophysical context. Exoplanet research (e.g. Winn & Fabrycky 2015; Santos & Buchhave 2018) and Galactic archaeol- Stellar characterisation is a matter of fundamental importance in ogy (e.g. Miglio et al. 2017) are examples of areas where studies often rely on the knowledge of fundamental stellar properties, such as the stellar mass, radius, and age. The advent of space-based aster- ★ E-mail: [email protected]

© 2015 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

Rapporteur Talk: Cosmic Ray Direct

Philipp Mertsch푎,∗ 푎Institute for Theoretical Particle Physics and Cosmology (TTK), RWTH Aachen University, 52056 Aachen, Germany E-mail: [email protected]

This is the report on the cosmic ray direct track of the 37th International Cosmic Ray Conference (ICRC 2021), broadly covering the contributions relating to charged cosmic rays (CRs) of Galactic origin. The contributions highlighted here are of both observational and theoretical nature and aim at interpreting the local fluxes of CRs as well as studying the wider dynamical effects of CRs. New data from space-borne experiments on CR electrons and positrons, proton and helium as well as heavier nuclei and their isotopes are reviewed. We cover some models of CR acceleration and their feedback on Galactic scales. Diffuse emission in 훾-rays as far as it concerns the CR spectra in the immediate vicinity of the solar system and CR anisotropies are covered briefly. Some upcoming and proposed experiments are highlighted towards the end. arXiv:2110.03325v1 [astro-ph.HE] 7 Oct 2021

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

∗Presenter

CRD rapporteur Philipp Mertsch

1. Introduction

Cosmic rays (CRs) are messengers of the violent universe and interpretations of CR measurements can be used to infer the properties of such environments as far as overall energetics, structure and dynamics are concerned. In addition, CRs have a dynamical effect in a variety of environments, shaping not only the medium that they themselves travel in, but also affecting other constituents, e.g. magnetic fields, radiation and the thermal gas. These two roles of CRs as spectators and of actors allow us to pose different kinds of questions. For the former role, the most pressing questions are: What are the sources of CRs? Can we find dark matter in CRs? Is there primordial anti-matter? For the latter role, different effects of CRs can be explored: CRs produce diffuse emission, contribute to ionisation and heating, provide gravitational support, drive winds and generate turbulence. Suffice it to say, what constitutes a valid and interesting model various strongly between these two different directions of research. The spectrum of CRs extends over more than 12 orders of magnitude in energy, from tens of MeV to about 1020 eV. This spectrum is largely a featureless power law, but there are two important spectral breaks around ∼ 3 × 1015 eV and ∼ 5 × 1018 eV. Both features have been argued to constitute evidence for the transition from Galactic to extra-galactic sources to take place at either energy. At energies than ∼ 10 MeV/ per nucleon the flux is dominantly of solar origin. At International Cosmic Ray Conferences (ICRCs), the scientific contributions on charged CRs have traditionally been distinguished by the kind of detection technique that is used in the respective energy regimes. To a certain degree, this distinction falls along the lines of Galactic vs. extra-galactic. Up to hundreds of TeV, direct measurements are possible, that is a particle detector, oftentimes consisting of (at least) trackers, calorimeters and scintillation sensors is flown on a balloon or in space. This allows not only for determining the energy of an event, but also of the charge and ideally the mass on an event-by-event basis. Of course, ideally, direct measurement are desirable at all energies. However, given the quick, ∼ 퐸−3 decrease of (differential) intensities with energy 퐸, the number of events for the areas/volumes typically instrumented becomes too low at energies above a PeV and higher. This is the report for the cosmic ray direct (CRD) track of the 2021 edition of the ICRC that was to be held in Berlin, Germany, but was eventually exclusively online. Naturally, this cannot be an exhaustive review of all contributions to the CRD track, and is unavoidably somewhat subjective. Having said that, we tried to cover the most important experimental results and the most controversially discussed topics. We will keep references to works other than contributions from this ICRC to a minimum. Instead, we refer the interested reader to a few of the recent reviews on the field [1–3], in addition to the classical texts, e.g. [4]. Reports on the other tracks are found elsewhere [5–12]. Usually at ICRCs, most contributions are presented in parallel during short oral presentations with limited amount for questions and discussions, in addition to parallel sessions with review and highlight talks. Given the exclusive online nature of this ICRC, the scientific content at this year’s conference was organised differently. Most oral presentations were pre-recorded and made available to all participants ahead of time. This allowed significantly more time for discussion in a number of live discussion sessions. Overall, it is this rapporteur’s opinion that this concept was a great success. This is partly due to the minute planning on part of the local organising committee;

2 arXiv:2110.03322v1 [astro-ph.GA] 7 Oct 2021 Sweden. is h tegho h antcfil smr lsv.Snecharge Since elusive. more is th field determined magnetic reliably the of be strength often the can tics, which temperature, and density mlctoso h o rqec unoe ntespectrum the in turn-over frequency low the of Implications 1 eateto srnm,AbNv nvriyCne,Stockho Center, University AlbaNova Astronomy, of Department antcfilspa nipratrl nmn srpyia phenome astrophysical many in role important an play fields Magnetic hra aito ore 11)—mgei field magnetic — (1119) sources radiation thermal shock headings: the Subject and C knot radio supernovae. of in emission properties Attention emiss radio the phase. radio between pre-acceleration the ne similarities its Hence, that of the below properties acceleration. is the Fermi electrons to order emitting sights first radio pr in the observe is part of the take gas that energy thermal suggested the strengt is the that field It implies of magnetic electrons. the cooling relativistic fo enhance of addition, account further density In to to needed order shocks. pressure in synchr necessary high to the the for due that factor is filling argued vations volume is low It a depend Since implies emission. the index. this of spectral resolved, an synchrotron omission partially the to an is on due with field is together magnetic ma value minimum volume of lower the emission orders earlier, two the roughly The mag of is its timate which estimated. that mG, 10 previously show least to than at used larger be is to This likely is turn-over. strength frequency low a has Tau h ycrto pcrmo ai ntCi h rtselrobje protostellar the in C knot radio of spectrum synchrotron The lsia ar tr 22 tla es(67 non- — (1607) jets Stellar — (252) stars Tauri T Classical ai ntCi GTau DG in C knot radio [email protected] .Introduction 1. .I BjörnssonC.-I. ABSTRACT 1 mUiest,S–0 1Stockholm, 91 SE–106 University, lm atce r accelerated are particles d og pcrldiagnos- spectral rough pcrlindex spectral d sas rw to drawn also is swl sthe as well as h o ie in- gives ion h obser- the r a ncnrs to contrast In na. h source the ei field netic induced gnitude ddto eded neof ence overes- obably tDG ct otron of DRAFTVERSION OCTOBER 8, 2021 Typeset using LATEX RNAAS style in AASTeX631

Strong lensing: A magnifying glass to detect gravitational-wave microlensing

Eungwang Seo,1, ∗ Otto A. Hannuksela,2, 3 and Tjonnie G. F. Li1

1Department of Physics, The Chinese University of Hong Kong, Shatin, N.T., Hong Kong 2Nikhef – National Institute for Subatomic Physics, Science Park, 1098 XG Amsterdam, The Netherlands 3Institute for Gravitational and Subatomic Physics (GRASP), Department of Physics, Utrecht University, Princetonplein 1, 3584 CC Utrecht, The Netherlands

(Dated: October 8, 2021)

ABSTRACT Microlensing imprints by typical stellar-mass lenses on gravitational waves are challenging to identify in the LIGO–Virgo frequency band because such effects are weak. However, stellar-mass lenses are generally embedded in host galaxies so that strong lensing can accompany microlensing. Therefore, events that are strongly lensed in addition to being microlensed may signiﬁcantly improve the inference of the latter. In particular, since a pair of strongly lensed signals have the same underlying gravitational-wave signal, we can use information from one signal to enhance the inference of the microlensing effects of the other signal. This will signiﬁcantly enhance our future ability to detect the weak imprints from stellar-mass objects on gravitational-wave signals from colliding compact objects.

INTRODUCTION In this letter we will show how strong lensing cooperates The general microlensing effect on gravitational wave in improving our ability to detect microlensing effects. We (GW) gives observable beating patterns on the signals in light will get information from the first strongly lensed signal by of current detector sensitivity (Christian et al. 2018; Jung & conducting parameter estimation (PE) to analyze the second Shin 2019; Cao et al. 2014; Diego, J. M. et al. 2019; Diego signal, which is both microlensed and strongly lensed. It is 2020). However, the microlensing signatures due to stellar- feasible since the two signals have the same waveform pa- mass objects are weak (Cheung et al. 2021; Mishra et al. rameters. Finally, we will show that we can detect the mi- 2021), and it is thus challenging to obtain compelling evi- crolensing effects on the second signal at an improved accu- dence of such signatures. racy. Typical stellar-mass objects are predominantly part of METHODOLOGY larger-scale structures, which may induce strong lensing effects on GWs. Therefore, it is realistic that strong lensing In principle, one or both of the strongly lensed gravita- signatures are observed in conjunction with microlensing sig- tional waves in Fig.1 can be microlensed. However, we natures. Suppose that a GW from a distant black hole binary consider the cases where only one and both of them are mi- undergoes strong lensing by a galaxy hosting microlenses. crolensed separately. The reason we consider the scenario The GW is then split into two or more signals arriving min- explicitly where only one of the signals is microlensed is utes to months apart (Oguri 2018; Li et al. 2018; Ng et al. that the scenario is more probable. In particular, even if 2018), where each signal may include different beating pat- both signals have the same probability of being microlensed, arXiv:2110.03308v1 [astro-ph.HE] 7 Oct 2021 terns caused by the field of stellar-mass microlenses along the probability that only one of the signals is microlensed their respective paths (Fig.1, for an illustration). In such sce- is higher than the probability that both of them are. More narios, the magnification by the microlens is much weaker specifically, the typical lensing rate by stellar-mass objects is than magnification by the strong lens (Cheung et al. 2021) small because its Einstein radius, which is related to the lens- so that we can infer both the waveform and the relative im- ing cross-section, is small (Christian et al. 2018). Further- age properties (e.g. relative magnification and phase shift) more, microlensing is less likely to happen for signals with 1 from the two strongly lensed signals (Liu et al. 2021; Lo & negative parity (or saddle points) because the demagnifica- Hernandez 2021; Janquart et al. 2021; Abbott et al. 2021a). tion region by a microlens is much larger than the magnifica-

1 We can determine which signals have negative parity by their time de- ∗ E-mail: [email protected] lays (Wierda et al. 2021). manuscript submitted to The Planetary Science Journal

Carbon cycling and habitability of massive Earth-like exoplanets

Amanda Kruijver1,∗, Dennis H¨oning1,2,∗, Wim van Westrenen1

1Department of Earth Sciences, Vrije Universiteit Amsterdam, The Netherlands 2Origins Center, The Netherlands

∗AK and DH equally contributed to this paper. arXiv:2110.03295v1 [astro-ph.EP] 7 Oct 2021

–1– Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX62

Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes

Jared A. Goldberg,1 Yan-Fei Jiang(姜燕飞),2 and Lars Bildsten1, 3

1Department of Physics, University of California, Santa Barbara, CA 93106, USA 2Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA 3Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA

ABSTRACT We explore the three-dimensional properties of convective, luminous (L 104.5 105L ), Hydrogen- rich envelopes of Red Supergiants (RSGs) based on radiation hydrodynamic≈ simulations− in spherical geometry using Athena++. These computations comprise 30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational≈ potential for the outer 3M ≈ of the simulated M 15M stars. This work reveals a radius, Rcorr, around which the nature of ≈ the convection changes. For r > Rcorr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well-studied in less luminous stars, but in RSGs, the near- (or above-) Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outwards rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and early lightcurves of Type II-P supernovae. Inside of Rcorr, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length-theory (MLT). Radiation pressure provides 1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length≈ of α = 3 4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability− of these 3D models is mostly on the timescale of the convective plume lifetimes ( 300 days), with amplitudes consistent with those observed photometrically. ≈

Keywords: hydrodynamics — radiative transfer — convection — stars: massive — stars: supergiants

1. INTRODUCTION It is a theoretical challenge to realistically model stars, As massive (M > 10M ) stars leave the main or even parts of stars, in 3D. This is especially true when sequence, they expand∼ to become Red Supergiants radiative transfer must also be simultaneously solved (RSGs), reaching radii of 300 1200R and luminosi- through a highly turbulent medium with large density ties of 104 105.5L (e.g.≈ Levesque− et al. 2006; Drout variations over optical depths ranging from far above et al. 2012≈ ; Massey− et al. 2021), approaching the Edding- unity down to the radiating photosphere. This radia- ton limit and receiving increasing hydrostatic support tion hydrodynamic (RHD) challenge has been very well- from radiation pressure. These stars are characterized addressed in cases where this region is close to plane- by diﬀuse convective hydrogen-rich envelopes with large parallel, starting with the fundamental work for the Sun scale heights (H/r 0.3) and sonic convection near (Stein & Nordlund 1989, 1998), and now ranging across their surfaces. They≈ are intrinsically variable and pul- the HR diagram (e.g. Trampedach et al. 2013, 2014a,b; Magic et al. 2013a,b, 2015; Chiavassa et al. 2018a; Sonoi arXiv:2110.03261v1 [astro-ph.SR] 7 Oct 2021 sate in large-amplitude coherent modes (e.g. Kiss et al. 2006; Soraisam et al. 2018; Chatys et al. 2019; Ren et al. 2019), building on earlier 2D RHD work (see Lud- et al. 2019; Dorn-Wallenstein et al. 2020) and their 3D wig et al. 1999 for an excellent summary). These 3D nature is revealed to us in spectro-interferometric obser- computations have yielded a physical understanding of vations of nearby stars (e.g. Arroyo-Torres et al. 2015; the nature of RHD convection in this limit, and pro- Kravchenko et al. 2019, 2021; Montarg`eset al. 2021; vide a quantitative ability to set the outer boundary Norris et al. 2021). condition in 1D stellar models (e.g. Trampedach et al. 2014a; Salaris & Cassisi 2015; Magic 2016; Mosumgaard et al. 2018; Spada et al. 2021) for log g > 1.5, including Corresponding author: J. A. Goldberg for asteroseismic applications (Mosumgaard∼ et al. 2020). [email protected] Such clarity has not yet been reached in stellar convec- X-ray luminosity function of accreting neutron stars and black holes

Konstantin Postnov1, Alexandre Kuranov1, Lev Yungelson2, and Marat Gil’fanov3

1 Sternberg Astronomical Institute, 14 Universitetskij pr., Moscow, Russia [email protected] 2 Institute of Astronomy, RAS, 48 Pyatnitskaya str., Moscow, Russia 3 Space Research Institute, RAS, 84/32 Profsoyuznaya str., Moscow, Russia

Abstract. We model X-ray luminosity functions (XLF) of accreting 35 41 −1 neutron stars and black holes in 10 ≤ LX ≤ 10 erg s range in star- forming galaxies and galaxies with the initial star formation burst. XLFs are obtained by combining a fast generation of compact object+normal star population using the binary population synthesis code BSE and calculation of the subsequent detailed binary evolution by the MESA code. XLFs in the galaxies of both types is broadly reproduced using the standard assumptions of the binary star evolution.

Keywords: neutron stars, black holes, X-ray sources

1 Introduction

Stellar-mass ﬁeld X-ray sources in galaxies are usually identiﬁed with binaries harbouring neutron stars (NS) or black holes (BH) accreting matter from their companions. Here, we consider the X-ray luminosity functions (XLF) of model galaxies with an instantaneous star-formation burst as a proxy for elliptical galaxies and a model with constant star-formation rate (SFR) as a proxy for star- forming galaxies. In §2, we present our assumptions and method of computations, in §3 we show the results of computations and discuss them. arXiv:2110.03255v1 [astro-ph.HE] 7 Oct 2021 2 Method of Computations

The matter from the optical star in a binary system with NS or BH can be transferred to the compact object in different ways – via the Roche lobe overflow (RLOF) of the donor star or by the gravitational capture of the stellar wind from the optical companion. In the case of RLOF onto magnetized NSs, we take into account that the NS magnetic fields of about 1012 ∼< B ∼< 1014 G modify the standard Shakura-Sunyaev disk accretion pattern. This can result in a super-Eddington X-ray luminosity LX from the accreting NS even before the Eddington accretion rate is attained (Kuranov et al. 2020). This effect allows Astronomy & Astrophysics manuscript no. aanda ©ESO 2021 October 8, 2021

The Aperture Array Veriﬁcation System 1: System overview and early commissioning results P. Benthem1, R. Wayth2, E. de Lera Acedo12, K. Zarb Adami11, 13, M. Alderighi5, C. Belli4, P. Bolli4, T. Booler2, J. Borg11, J. W. Broderick2, 1, S. Chiarucci4, R. Chiello13, L. Ciani10, G. Comoretto4, B. Crosse2, D. Davidson2, A. DeMarco11, D. Emrich2, A. van Es14, D. Fierro7, A. Faulkner12, M. Gerbers1, N. Razavi-Ghods12, P. Hall2, L. Horsley2, B. Juswardy2, D. Kenney2, K. Steele2, A. Magro11, A. Mattana3, B. McKinley2, J. Monari3, G. Naldi3, J. Nanni9, P. Di Ninni4, F. Paonessa8, F. Perini3, M. Poloni3, G. Pupillo3, S. Rusticelli3, M. Schiaﬃno3, F. Schillirò6, H. Schnetler15, R. Singuaroli10, M. Sokolowski2, A. Sutinjo2, G. Tartarini9, D. Ung2, J. G. Bij de Vaate1, G. Virone8, M. Walker2, M. Waterson14, S. J. Wijnholds1, and A. Williams2

1 ASTRON - Dwingeloo, The Netherlands e-mail: [email protected] 2 ICRAR/Curtin University - Perth, Australia 3 INAF/IRA - Institute of RadioAstronomy - Bologna, Italy 4 INAF - Osservatorio Astronomico di Arcetri - Florence, Italy 5 INAF - Istituto di Astroﬁsica Spaziale e Fisica Cosmica - Milan. Italy 6 INAF - Osservatorio Astronomico di Catania - Catania, Italy 7 INAF - Headquarter - Rome, Italy 8 CNR-IEIIT - Turin, Italy 9 Bologna University DEI - Bologna, Italy 10 University of Florence DINFO - Florence, Italy 11 Institute of Space Sciences and Astronomy, University of Malta - Msida, Malta 12 University of Cambridge - Cambridge, United Kingdom 13 University of Oxford - Oxford, United Kingdom 14 SKA Organization HQ - Jodrell Bank, Cheshire, United Kingdom 15 ATC/STFC - Edinburgh, Scotland

October 8, 2021

ABSTRACT

The design and development process for the Square Kilometre Array (SKA) radio telescope’s Low Frequency Aperture Array component was progressed during the SKA pre-construction phase by an international consortium, with the goal of meeting requirements for a critical design review. As part of the development process a full-sized prototype SKA Low ‘station’ was deployed – the Aperture Ar- ray Veriﬁcation System 1 (AAVS1). We provide a system overview and describe the commissioning results of AAVS1, which is a low frequency radio telescope with 256 dual-polarisation log-periodic dipole antennas working as a phased array. A detailed system description is provided, including an in-depth overview of relevant sub-systems, ranging from hardware, ﬁrmware, software, calibration, and control sub-systems. Early commissioning results cover initial bootstrapping, array calibration, stability testing, beam-forming, and on-sky sensitivity validation. Lessons learned are presented, along with future developments. Key words. Radio astronomy – Square Kilometre Array – phased-array – telescopes — instrumentation: miscellaneous

1. Introduction niques. Likewise, the core operating concept for SKA-Low con- verged on the aperture-array design; the telescope consists of The Square Kilometre Array (SKA) radio telescope promises hundreds of thousands of ﬁxed antennas arranged in several hun- arXiv:2110.03217v1 [astro-ph.IM] 7 Oct 2021 to be a transformational instrument for astronomy and astro- dred ‘stations’, where each station acts as one or more phased physics in the coming decade (Braun et al. 2015). The SKA tele- arrays. Modern digital electronics allow great ﬂexibility in the scope is split into the Mid-Frequency Array (SKA-Mid) and the beam-forming process. The current design for SKA-Low con- Low Frequency Array (SKA-Low), the latter to be constructed at sists of 512 stations, each of which is composed of 256 antenna the Murchison Radio-astronomy Observatory (MRO) in Western elements (Dewdney 2016). Australia. It will be the most sensitive instrument ever built in its operational frequency range. In the ‘pre-construction’ phase of the SKA, a consortium SKA-Low will operate between 50 and 350 MHz. At these called the Aperture Array Design and Construction (AADC) radio frequencies, it is feasible and practical to build very large Consortium was responsible for progressing the design of SKA- collecting areas with simple dipole-like antennas operating as a Low up to the formal system engineering critical design Review phased array. Many early-generation radio telescopes were con- (CDR) milestone. Within the high-level overall SKA telescope structed with such a design using analogue beam-forming tech- design, the AADC Consortium was responsible for the Low Fre-

Article number, page 1 of 21 arXiv:2110.03204v1 [astro-ph.SR] 7 Oct 2021 [email protected] orsodn uhr bihkRajhans Abhishek author: Corresponding lge yrdnmcsmltossrea necetto o study eg., (see for event, tool heating efficient dependent an time as generic serve simulations hydrodynamic aligned Dmdl) nti okw etitorevst h omnyused commonly the to ourselves restrict Hence (se we ( work loop Loops the this loop. of In the part models). along coronal 0D the much for vary characteristic not hav are codes pressure do 0D parameters issues, physical f rea such simulations more part, overcome aligned To Additionally field Performing expensive. required. ones. computationally lik are monolithic parameters runs simple of to of variation opposed number of effects large the a study function, to needs one where ieeta msinmauei h rniinrgo.Tecrnlb coronal The region. transition the in measure ba emission the at differential velocity and pressure, density, temperature, average bihkRajhans Abhishek h oa ooahsamra flo iesrcue,wihaemagn are which structures, like loop of myriad a has corona Solar The nanutshell, a In yee sn L using 2021 Typeset 8, October version Draft EBTEL ntnaeu ieeta msinmauedsrbto ntetran code, the the in o of version distribution part fluid measure coronal emission the differential over instantaneous averaged pressure, and temperature, eesr oacutfrti ffc.Hr,w upgrade we Here, impu the effect. during this flows for account supersonic to of necessary presence the show solutions the rdcssproi os sn h imthsi h ouin ep we solution, the in which mismatches under the conditions Using the flows. supersonic However, predicts subsonic. are flows h eoiispeitdby predicted velocities the EBTEL2 transit and corona the between int energy space and ( on mass of loops rely exchange of to the evolution is thermal solutions based costly, physics-based enthalpy computationally quick is but structures approximate these of evolution dynamical emi h airSoe qain ecmaetesltosfrom solutions the compare We equation. Navier-Stokes the in term npesr between pressure in lsafildlo tutrsaecmo nteslrcrn.Becaus corona. solar the in common are structures loop filled Plasma 1 ne-nvriyCnr o srnm n srpyis P Astrophysics, and Astronomy for Centre Inter-University see ; 4 2 AAGdadSaeFih etr oa hsc Laboratory Physics Solar Center, Flight Space Goddard NASA A ieUiest,Dprmn fPyisadAtooy 6100 Astronomy, and Physics of Department University, Rice EBTEL swl stesaeo-h-r edaindhdoyaiscode hydrodynamics aligned field state-of-the-art the as well as T E 3 etrfrAstrophysics for Center eal tl nAASTeX63 in style default X lmhke l 2008 al. et Klimchuk , 1 tde h asadeeg xhnebtentecrn n th and corona the between exchange energy and mass the studies ugs Tripathi Durgesh lw nEtap ae hra vlto fLoops of Evolution Thermal Based Enthalpy in Flows EBTEL3 EBTEL2 EBTEL3 EBTEL and EBTEL3 | ; sue hta l tgsteflw r usnc oee,som However, subsonic. are flows the stages all at that assumes , HYDRAD avr mtsna,6 adnS.CmrdeM 23,USA 02138, MA Cambridge St. Garden 60 Smithsonian, & Harvard agl ta.2012b al. et Cargill , 1 a eue osuyflw ntesystem. the in flows study to used be can r ncoeareetwt hs bandwith obtained those with agreement close in are tpe .Bradshaw J. Stephen Rcie;Rvsd Accepted) Revised; (Received; ouin eit usatal from substantially deviate solutions lmhk2006 Klimchuk 1. sbte hnta between that than better is umte oApJ to Submitted INTRODUCTION ABSTRACT EBTEL rmwr sacmol sdmto ostudy to method used commonly a is framework ) ; eo h ooa op.I locmue h instantaneous the computes also It loops. coronal the of se e EBTEL2 anse l 2016 al. et Barnes , agl ta.2012a al. et Cargill 2015 s a ,Gnskid ue410,India 411007, Pune Ganeshkhind, 4, - Bag ost , oplnt,eeg ugtadpol fheating of profile and budget energy length, loop e 2 endvlpd h dai hti h coronal the in that is idea The developed. been e ; ia .Kashyap L. Vinay oe61 rebl,M 07,USA 20771, MD Greenbelt, 671, Code , to op eoiya ooa ae n the and base, coronal at velocity loop, f s falo sdfie stepsto twhich at position the as defined is loop a of ase n h epneo lsai ooa op oa to loops coronal in plasma of response the ing anSre oso,T 70,USA 77004, TX Houston, Street Main el 2014 Reale Dcd:Etap ae hra vlto of Evolution Thermal Based Enthalpy code: 0D h auso vrg est,tmeaueand temperature density, average of values the itcseaisivlemlisrne op as loops multi-stranded involve scenarios listic EBTEL3 rsc opiae u elsi cnro is scenarios realistic but complicated such or o region. ion sv hs fha nu.I sthus is It input. heat of phase lsive necetmto o computing for method efficient an EBTEL2 iinrgo.Tecretsingle- current The region. sition tcflxtbsfildwt lsa Field plasma. with filled tubes flux etic EBTEL3 HYDRAD ooeacieint determine to criterion a ropose yicuigtekntcenergy kinetic the including by gae Dsmltos The simulations. 0D egrated o details). for , o eiw.Hwvri situations in However review). a for , o ealddsuso nvarious on discussion detailed a for , ealdmdln fthe of modeling detailed e and ihtoeotie using obtained those with HYDRAD EBTEL efidta h match the that find We . rniinrgo.I calculates It region. transition e , 3 HYDRAD and swe h latter the when ’s ovsfrdensity, for solves HYDRAD ae .Klimchuk A. James Additionally, . hnthe when etimes 4 MNRAS 000,1–18 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

Epoch of Reionization Power Spectrum Limits from Murchison Wideﬁeld Array Data Targeted at EoR1 Field

M. Rahimi,1,2★ B. Pindor,1,2, J. L. B. Line,3,2, N. Barry,1,2, C. M. Trott,3,2, R. L. Webster,1,2, C. H. Jordan,3,2 M. Wilensky,4, S. Yoshiura,1,5, A. Beardsley,6, J. Bowman,7, R. Byrne,4, A. Chokshi,1,2, B. J. Hazelton,4,8 K. Hasegawa,9, E. Howard,10, B. Greig,1,2, D. Jacobs,7, R. Joseph,3,2, M. Kolopanis,7, C. Lynch,3,2 B. McKinley,3,2, D. A. Mitchell,10,2, S. Murray,7, M. F. Morales,4, J. C. Pober,11, K. Takahashi,12 S. J. Tingay,3,2 R. B. Wayth,3,2, J. S. B. Wyithe,1,2, Q. Zheng,13

1The University of Melbourne, School of Physics, Parkville, VIC 3010, Australia 2ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO-3D) 3International Centre for Radio Astronomy Research, Curtin University, Perth, WA 6845, Australia 4University of Washington, Department of Physics, Seattle, WA 98195, USA 5Mizusawa VLBI Observatory, National Astronomical Observatory Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 6Winona State University, Department of Physics, Winona, MN 55987, USA 7Arizona State University, School of Earth and Space Exploration, Tempe, AZ 85287, USA 8University of Washington, eScience Institute, Seattle, WA 98195, USA 9Graduate School of Science, Nagoya University, Japan 10CSIRO Astronomy and Space Science (CASS), PO Box 76, Epping, NSW 1710, Australia 11Brown University, Department of Physics, Providence, RI 02912, USA 12Faculty of Science, Kumamoto University, 2-39-1 Kurokami, Kumamoto 860-8555, Japan 13Shanghai Astronomical Observatory, China

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is instructive to explore different aspects of existing datasets. One such aspect is the choice of observing field. To date, MWA EoR observations have largely focused on the EoR0 field. In this work, we present a new detailed analysis of the EoR1 field. The EoR1 field is one of the coldest regions of the Southern radio sky, but contains the very bright radio galaxy Fornax-A. The presence of this bright extended source in the primary beam of the interferometer makes the calibration and analysis of EoR1 particularly challenging. We demonstrate the effectiveness of a recently developed shapelet model of Fornax-A in improving the results from this field. We also describe and apply a series of data quality metrics which identify and remove systematically contaminated data. With substantially improved source models, upgraded analysis algorithms and enhanced data quality metrics, we determine EoR power spectrum upper limits based on analysis of the best ∼14-hours data observed during 2015 and 2014 at redshifts 6.5, 6.8 and 7.1, with the lowest 2휎 upper limit at z=6.5 of Δ2 ≤ (73.78 mK)2 at 푘 = 0.13 h Mpc−1,

arXiv:2110.03190v1 [astro-ph.CO] 7 Oct 2021 improving on previous EoR1 measurement results. Key words: cosmology: dark ages, reionization, ﬁrst stars – techniques: interferometric – methods: observational, data analysis

1 INTRODUCTION & Loeb 2012). Our goal is to measure the HI spin temperature and neutral fraction via its 21-cm signal to follow the evolution of neutral Over cosmic time, hydrogen has gone through three major transitions: gas in the IGM during the Epoch of Reionization (EoR). ionized hydrogen post-nucleosynthesis in the primordial plasma, The GMRT (Giant Metrewave Radio Telescope, Paciga et al. neutral hydrogen as the universe cooled by expansion, and ﬁnally (2013)), LOFAR (LOw Frequency Array1, Yatawatta et al.(2013)), reionized hydrogen as a result of the formation of the ﬁrst stars and galaxies (Furlanetto et al. 2006; Morales & Wyithe 2010; Pritchard

★ E-mail: [email protected] 1 http://www.lofar.org

© 2021 The Authors arXiv:2110.03185v1 [astro-ph.GA] 7 Oct 2021 n otglx,aeepce ob oecmo (Bour- common more be to this satellites expected are in between galaxy, effective interactions host more and and and mergers violent minor are sense, though mergers in Even major formation evolution. the chemical star their triggering thus and and galaxies, gas compressing in the role of of evolution one 2001). and (Wyder is formation galaxies popula- (SFH) the dwarf star history these of old tracers formation powerful and star most age are The intermediate there Simi- with young, tions. dIrrs considered episodes. few are formation a dIrrs star all pop- although multiple popu- age larly, to intermediate old due dominant exclusively Whereas a ulation have possess to galaxies. some assumed two lations, often these are from dSphs and (dIrrs) distances galaxy irregulars farther dwarf 31) (M late-type at and Andromeda (dSphs) can (MW), the Way spheroidals of they Milky dwarf satellites Often early-type are that into evolution. classified galaxy stel- be and between populations connection the lar understand to laboratories H SA ETNTLSOEMNTRN UVYO OA GROUP LOCAL OF SURVEY MONITORING TELESCOPE NEWTON ISAAC THE nvriyo acetr 1 P,UK 9PL, M13 Manchester, of University Iran 1458889694, Iran Tehran, (IPM), Sciences damental Iran; Tehran, 1956836613, Box P.O. (IPM), [email protected] Sciences mental aayitrcin n egr lya important an play mergers and interactions Galaxy excellent are (LG) Group Local the in galaxies Dwarf rpittpstuigL using 2021 typeset 8, Preprint October version Draft H TRFRAINHSOYO NRMD EIE RMLN P LONG FROM DERIVED I ANDROMEDA OF HISTORY FORMATION STAR THE la Saremi Elham 5 4 3 2 1 orl akCnr o srpyis lnTrn Building Turing Alan Astrophysics, for Centre Bank Jodrell Tehra Tecnology, of University Sharif Department, Physics Fu in Research for Institute Funda- Observatory, UK National 5BG, in Iranian ST5 University, Keele Research Laboratories, Lennard-Jones for Institute Astronomy, of School M Keywords: tla vlto oes eetmtdteSHo hsglx.Amajor A 0 reaching galaxy. ago, this Gyr 6.6 of around SFH peaking the half- I two And estimated within (And in we candidates I use Andromeda LPV occurred models, 59 for methodology identified evolution it the our derive stellar present Using then we survey. and the papers, stars of of a (LPV) series variable production a period to dust of this the second relate quantify and this distributions, distribution (SFH), energy their history spectral sup multi-wavelength use formation red instability, and star pulsational of stars their the aims (AGB) on main branch globula carrie The based giant (WFC). isolated was evolution asymptotic Camera four evolved galaxies Field most and Wide dwarf the the galaxies with nearby dwarf observed were the 55 (LG) in (INT). Telescope survey Newton monitoring optical An nraei eetsa omto ekn rud80Mrao ee g We this ago. in Myr LPVs 800 ( dusty around ago some peaking Gyr of formation 4 presence star The recent in ago. increase Gyr 1–2 until declining ⊙ scluae ihntohl-ih ai fAdIfracntn metallicit constant a for I And of radii half-light two within calculated is 1 1. tfhJavadi Atefeh , tr:eouin–sas G n P–sas uioiyfunction, luminosity stars: LPV– and AGB stars: – evolution stars: < z asls tr:oclain aais niiul n aais s galaxies: – I And individual: galaxies: – Group oscillations Local galaxies: stars: – mass-loss A INTRODUCTION T E tl mltajv 12/16/11 v. emulateapj style X 0 . ) hc ae n aeqecigdp.Attlselrms (1 mass stellar total A dSph. late-quenching a I And makes which 5), 1 ade Navabi Mahdieh , onriArbab Bojnordi rf eso coe ,2021 8, October version Draft ABSTRACT 1 n- ac h a Loon van Th. Jacco , n , 3,4 anMcDonald Iain , ta.21) ae nasuyo h e horizontal red the of study a on Based (Skillman dSphs MW 2017). 31 the M al. than the of uniform et SFHs more espe- galaxies. are them, massive dSphs of less most the in among found dwarf cially been in not formation have pres- star they the early galaxies, While active indicates in- galaxies. GCs dwarf (2017) of LG luminous ence of al. as SFHs (GCs) et the clusters of Caldwell tracers globular of time. use the that vestigated by formation star > M mass s.Wize l 21a siae htdafgalaxies dwarf that estimated (2014a) al. that ex- with do et concluded trends Weisz and formation, galaxies star ist. stud- dwarf diverse the have (1998) LG they of Mateo of although age Overall, number galaxies. the a these over many in ied 2020). formation seen 2017), al. are star al. Universe et of et (Saremi periods Skillman extended remain 2015; questions 2014a; open al. et Weisz al. galaxy et host (Ruiz-Lara the dynamics of internal build-up its 2020). the and satel- in mass dwarf effective stellar Therefore be content. with can stellar sim- lites coincided orbit Sgr MW to and according the ulations passages in pericentre Sgr formation proposed star episodes that of enhanced found orbit They of the (Sgr). also, galaxy and dwarf 2018) Sagittarius Collaboration Gaia SFH (e.g., detailed vey the studied the (2020) of streams al. Re- stellar et encounters. Ruiz-Lara or satellite cently, minor substructures these of halo confirmations MW are 2010). naud ept eetsuiso Gdafglxe’SH (e.g., SFHs galaxies’ dwarf LG of studies recent Despite < M > 10 . ∼ 0035 0Grao(ro ordhf ) hl hs with those while 2), redshift to (prior ago Gyr 10 7 p rudteSnbsdo h aasur- Gaia the on based Sun the around kpc 2 M 10 ⊙ uioiyadrdu aitos In variations. radius and luminosity ± 5 oetmt F ae nlong- on based SFH estimate to d a xeine ny3%o hi overall their of 30% only experienced had 5 dms osfo oeln the modelling from loss mass nd M 0 rinsa h n-on ftheir of end-point the at ergiants aut unhn iearound time quenching a valuate . 2 06M 0016 ⊙ ai .Khosroshahi G. Habib , ih ai fAdIadPadova and I And of radii light )dafglx sa example an as galaxy dwarf I) vrlmnst oreconstruct to luminosity over lsesi h oa Group Local the in clusters r omdmr hn8%o hi stellar their of 80% than more formed lx orsod oaslight a to corresponds alaxy u ihte25mIsaac 2.5-m the with out d hssre r oidentify to are survey this pc fsa formation star of epoch y ⊙ Z asfnto stars: – function mass yr 0 = − WR AAIS II. GALAXIES. DWARF 1 RO VARIABLES ERIOD . n nyslowly only and 0007. elrcontent tellar 6 ± 7) 1,3 × Behzad , 10 6 MNRAS 000, 000–000 (0000) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

The Parameter-Level Performance of Covariance Matrix Conditioning in Cosmic Microwave Background Data Analyses

L. Balkenhol ? and C. L. Reichardt School of Physics, University of Melbourne, Parkville, VIC 3010, Australia

8 October 2021

ABSTRACT Empirical estimates of the band power covariance matrix are commonly used in cosmic microwave background (CMB) power spectrum analyses. While this approach easily captures correlations in the data, noise in the resulting covariance estimate can systematically bias the parameter fitting. Conditioning the estimated covariance matrix, by applying prior information on the shape of the eigenvectors, can reduce these biases and ensure the recovery of robust parameter constraints. In this work, we use simulations to benchmark the performance of four different conditioning schemes, motivated by contemporary CMB analyses. The simulated surveys measure the TT , TE, and EE power spectra over the angular multipole range 300 ` 3500 in ∆` = 50 wide bins, for temperature map-noise levels of 10, 6.4 and ≤ ≤ 2 µK-arcmin. We divide the survey data into Nreal = 30, 50, or 100 uniform subsets. We show the results of different conditioning schemes on the errors in the covariance estimate, and how these uncertainties on the covariance matrix propagate to the best-fit parameters and parameter uncertainties. The most significant effect we find is an additional scatter in the best-fit point, beyond what is expected from the data likelihood. For a minimal conditioning strategy, Nreal = 30, and a temperature map-noise level of 10 µK-arcmin, we find the uncertainty on the recovered best fit parameter to be 1.3 larger than the apparent posterior width from the likelihood ( 1.2 larger than the uncertainty × × when the true covariance is used). Stronger priors on the covariance matrix reduce the misestimation of parameter uncertainties to < 1%. As expected, empirical estimates perform better with higher Nreal, ameliorating the adverse effects on parameter constraints. Key words: cosmology: cosmic background radiation – methods: data analysis – cosmology: cosmological parameters

1 INTRODUCTION timators are frequently used on simulated or real data to generate a robust estimate of the covariance matrix. While one Current cosmological surveys provide new insight on infla- can use simulations to produce a covariance estimate to any tion, dark matter, dark energy, and neutrino properties (Do- desired precision, this requires a good understanding of the delson et al. 2013; Abazajian et al. 2015; Planck Collabo- noise properties of the experiment, and, perhaps more impor- ration et al. 2018; Dutcher et al. 2021; Aiola et al. 2020). tantly, can be very computationally expensive. Dodelson & Cosmic microwave background (CMB) and large-scale struc- Schneider(2013a) note that modern surveys underestimate ture (LSS) observations complement one another by provid- their parameter uncertainties by ∼ 5-10%, due to only run- arXiv:2110.03180v1 [astro-ph.CO] 7 Oct 2021 ing snapshots of the universe at different times. The interplay ning O(100) simulations (e.g., Louis et al.(2017a); Henning of these two probes allows us to constrain structure formation et al.(2018a); Planck Collaboration et al.(2019)). Shrinkage and modified gravity models (Omori et al. 2019). estimators (Hamimeche & Lewis 2009), emulators (Schnei- A crucial ingredient in many data analyses is the covari- der et al. 2008; Morrison & Schneider 2013), and large-scale ance matrix, which encodes the uncertainties on the mea- mode-resampling methods (Schneider et al. 2011) have been sured data points. An accurate covariance matrix is essential developed to address the increasing computational cost of for producing reliable model constraints; errors in the co- simulations by reducing the number of realisations required variance matrix directly propagate into skewed cosmological to achieve adequate precision. parameters (Hartlap et al. 2007; Hamimeche & Lewis 2009; Taylor et al. 2013; Dodelson & Schneider 2013a). Alternatively, the measured data itself can be used to es- Since analytic calculations of the covariance matrix may timate the covariance matrix (e.g., Lueker et al. 2010a; Das struggle to capture the full complexity of the instrument, eset al. 2014; Crites et al. 2015; Henning et al. 2018b; Dutcher et al. 2021). This approach does not make any assumptions about the noise profile of the experiment and can avoid or re- ? Corresponding author: [email protected] duce the required number of computationally expensive simu-

Dynamics of relativistic radio jets in asymmetric environments

Patrick M. Yates-Jones,1,2,3★ Stanislav S. Shabala,1,2,3 Martin G. H. Krause2,1 1 School of Natural Sciences, Private Bag 37, University of Tasmania, Hobart, TAS 7001, Australia 2 Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatﬁeld, Herts AL10 9AB, UK 3 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D)

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We have carried out relativistic three-dimensional simulations of high-power radio sources propagating into asymmetric cluster environments. We offset the environment by 0 or 1 core radii (equal to 144 kpc), and incline the jets by 0, 15, or 45◦ away from the environment centre. The different environment encountered by each radio lobe provides a unique opportunity to study the effect of environment on otherwise identical jets. We find that the jets become unstable towards the end of the simulations, even with a Lorentz factor of 5; they nevertheless develop typical FR II radio morphology. The jets propagating into denser environments have consistently shorter lobe lengths and brighter hotspots, while the axial ratio of the two lobes is similar. We reproduce the recently reported observational anti-correlation between lobe length asymmetry and environment asymmetry, corroborating the notion that observed large-scale radio lobe asymmetry can be driven by differences in the underlying environment. Key words: hydrodynamics – galaxies: active – galaxies: jets – radio continuum: galaxies

1 INTRODUCTION lobe morphology is edge-darkened (Fanaroff-Riley class I, FR I) or edge-brightened (FR II). As radio telescope sensitivity and resolu- Radio sources are inextricably linked to their host galaxy through tion have increased, this distinction has become increasingly murky: feedback processes (McNamara & Nulsen 2007; Fabian 2012; Hard- compact radio sources with some low-frequency extended emission castle & Croston 2020) which influence black hole growth (Rafferty have been termed FR 0s (Garofalo & Singh 2019; Baldi et al. 2019; et al. 2006; Dubois et al. 2012; Alexander & Hickox 2012), star for- Capetti et al. 2019); low luminosity (퐿 ≤ 1025 W Hz−1) FR II mation (Gaibler et al. 2012; Weinberger et al. 2018), and the small to 150 radio sources have been observed (Mingo et al. 2019); radio sources large-scale environment through gas circulation, heating, and shock which exhibit hybrid FR I/II characteristics suggest that precession driving (Croton et al. 2006; Rafferty et al. 2008; Mukherjee et al. and projection effects can drastically alter the observed morphology 2016; English et al. 2016; Yates et al. 2018). Their influence on (Harwood et al. 2020; Krause et al. 2019; Horton et al. 2020a,b); the galaxy evolution and formation is underscored by the significant consideration of mergers, cluster weather, and remnant or restarting attention paid towards analytic (Scheuer 1974; Begelman & Cioffi radio sources further complicates the observed lobe morphology 1989; Falle 1991; Kaiser & Alexander 1997; Alexander 2002, 2006; (Mahatma et al. 2018; Yates et al. 2018; English et al. 2019; Joshi Shabala & Alexander 2009; Krause et al. 2012; Raouf et al. 2017), et al. 2019; Hardcastle et al. 2019a; O’Neill et al. 2019; Bruni et al. semi-analytic (Turner & Shabala 2015; Hardcastle 2018; Raouf et al. 2020). 2019), and numerical (Krause 2005; Gaibler et al. 2011; Wagner arXiv:2110.03162v1 [astro-ph.HE] 7 Oct 2021 et al. 2012; Hardcastle & Krause 2014; Mukherjee et al. 2016; Bick- Many complex processes are likely at work behind the ob- nell et al. 2018; Yates et al. 2018; Li et al. 2018; Perucho et al. 2019; served morphologies, however, one key factor is the environment Horton et al. 2020b) models to understand this phenomenon. These into which the radio source is expanding. Analytic models predict a radio sources are most commonly studied through observations of strong dependence on lobe size as a function of environment (Kaiser the kpc to Mpc scale radio lobes (Bennett & Simth 1962; Hardcas- & Alexander 1997; Alexander 2006; Krause et al. 2012), which is tle et al. 2019b; Seymour et al. 2020) inflated by the relativistic jet reproduced by simulations (Mendygral et al. 2012; Hardcastle & beam, throughout which a population of relativistic electrons accel- Krause 2013; Bourne & Sĳacki 2017; Yates et al. 2018). Obser- erated at shocks emits through synchrotron and inverse-Compton vationally, the observed length asymmetry in radio lobes pairs was processes (Kaiser et al. 1997; Turner et al. 2018a) recently shown to be linked to the underlying environment asymme- Radio sources have historically been classified into two distinct try, as traced by optical galaxy counts by Rodman et al.(2019), who categories (Fanaroff & Riley 1974) based on whether the large-scale found a statistically significant anti-correlation between the two. At first sight, it might sound trivial that jets propagate more slowly into a higher density environment. However, the environment ★ E-mail: [email protected] sets the pressure level in the radio lobes, which in turn affects

Dissecting Nearby Galaxies with piXedfit: I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

Abdurro’uf ,1 Yen-Ting Lin ,1 Hiroyuki Hirashita ,1 Takahiro Morishita ,2 Sandro Tacchella ,3 Masayuki Akiyama ,4 Tsutomu T. Takeuchi ,5, 6 and Po-Feng Wu 1, 7

1Institute of Astronomy and Astrophysics, Academia Sinica, 11F of AS/NTU Astronomy-Mathematics Building, No.1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan, R.O.C. 2Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 3Department of Physics, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Republic of Korea 4Astronomical Institute, Tohoku University, Aramaki, Aoba, Sendai 980-8578, Japan 5Division of Particle and Astrophysical Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464–8602, Japan 6The Research Center for Statistical Machine Learning, the Institute of Statistical Mathematics, 10-3 Midori-cho, Tachikawa, Tokyo 190—8562, Japan 7National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

(Received ; Revised ; Accepted ) Submitted to ApJ

ABSTRACT We study spatially resolved properties (on spatial scales of ∼ 1 − 2 kpc out to at least 3 effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously. With this self-consistent analysis, we present the spatially resolved version of the IRX–β relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust- to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

Keywords: Galaxy evolution (594) — Spiral galaxies (1560) — Interstellar medium (847) — Gas-to- arXiv:2110.03158v1 [astro-ph.GA] 7 Oct 2021 dust ratio (638)

1. INTRODUCTION complex manner. Stars are formed from the cold molec- A galaxy is a complex system whose baryonic compo- ular gas. During the lives of stars, metals are produced nent is composed of stars, dust, gas, and metals. These in the process of nucleosynthesis. Metals are ejected four key components are correlated with each other in a from stars to the interstellar medium (ISM) through the stellar wind and supernovae (SNe). Metals are the main ingredient for dust formation. The dust grains Corresponding author: Abdurro’uf form through condensation of metals in the stellar ejecta [email protected] (e.g., Dwek 1998; Zhukovska et al. 2008) that include winds of low mass stars in the asymptotic giant branch MNRAS 000,1–11 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style ﬁle v3.0

GMRT unveils steep-spectrum antique ﬁlaments in the galaxy cluster Abell 725

M. B. Pandge,1★ Ruta Kale,2 Pratik Dabhade,3,4 Mousumi Mahato3 and Somak Raychaudhury3,5 1DST-INSPIRE Faculty, Dayanand Science College, Barshi Road, Latur 413512, Maharashtra, India 2National Centre for Radio Astrophysics (TIFR), Pune 411007, Maharashtra, India 3Observatoire de Paris, LERMA, Collège de France, CNRS, PSL University, Sorbonne University, 75014, Paris, France 4Inter-University Centre for Astronomy and Astrophysics, Pune 411007, Maharashtra, India 5School of Physics and Astronomy, University of Birmingham,Birmingham B15 2TT, UK

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT We present original GMRT radio observations of the galaxy cluster Abell 725, at a redshift of 0.09, along with other archival observations. Our GMRT maps reveal two steep-spectrum diffuse filaments in the cluster, along with a previously reported arc-like structure, and a wide-angle tail (WAT) radio source associated with the Brightest Cluster Galaxy (BCG) at the periphery of the cluster. The bent morphology of the WAT indicates that its jets have been swept back by the dynamic pressure resulting from the motion of the BCG through the surrounding intracluster medium. The BCG associated with the WAT hosts a black hole whose mass we 9 0 estimate to be 1.4±0.4 × 10 M . We observe a 2 (195 kpc in projection) offset between the BCG and the X-ray centroid of the galaxy cluster, which, along with other dynamic features, indicates the cluster’s early stage of evolution. The WAT radio galaxy, the arc and the filaments 240 − ± − ± − ± − ± have spectral indices 훼612 = 0.46 0.15, 0.8 0.3, and ( 1.13 0.48, 1.40 0.50), respectively. The WAT and the arc are connected structures, while the filaments are detached from them, but are found to be along the trail of the WAT. Based on the morphology of the components, and the progressive steepening of the components from the core of the WAT to the filaments, we propose that this system is a radio galaxy with trailing antique filaments. Key words: galaxies:general– galaxies:active– X-rays:galaxies:clusters– galaxies:clusters:individual:ACO725

1 INTRODUCTION 2019). Apart from radio halos and relics, the only sources growing to Mpc scales are the giant radio galaxies (sizes > 0.7 Mpc), which Powerful diffuse radio sources are mostly found in major merging are in general not found in clusters (e.g. Dabhade et al. 2020). galaxy clusters, indicating that their origin is likely connected to The narrow-angle tail (NAT) and wide-angle tail (WAT) ra- the process of merger (Cohen & Clarke 2011). Clusters of galax- dio galaxies are generally found in the galaxy cluster environment, ies often show extended synchrotron emission in the intra-cluster arXiv:2110.03153v1 [astro-ph.CO] 7 Oct 2021 where the radio jets and lobes, in these systems, often appear to be medium (Feretti 2005). These radio sources mostly have very low bent in a common direction (Miley 1980; O’Dea & Owen 1985; surface brightness and steep radio spectrum, with their overall ex- Owen & Ledlow 1997; Patra et al. 2019). In cases where the angle tent growing to megaparsec scales. Such radio sources, associated 휙 between two tails of the radio galaxy is 휙 ≤ 90◦, forming a ‘V’ with clusters of galaxies, have been classified as radio halos and shape, it is referred to as NAT (e.g NGC 1262), while if the angle radio vestigials. The radio halos, with low polarisation, are found is 90◦ ≤ 휙 ≤ 180◦, it is known as a WAT and is usually found in to be at the centres of the galaxy clusters, while radio relics, found merging or not-relaxed clusters Burns(1990). The relative veloci- towards the periphery of clusters, with elongated arc-like shapes, ties of these sources, with respect to the ICM of the host cluster, are have emissions with higher levels of polarisation. The exact physics thought to be sufficient for the tails to be bent, due to the action of leading to the formation of these sources is indeed still a matter of ram pressure interactions (O’Dea & Owen 1985). debate, such sources are clearly associated with particle acceleration processes in the context of cluster mergers (e.g. van Weeren et al. Duty cycles of radio galaxies are a probe of the accretion at the supermassive blackhole. In the cases of powerful double radio sources of FR II kind, two or more phases can be identified due to the ★ E-mail: [email protected] (DSCL) well defined double structure formed by the hot-spots. However in

© 2021 The Authors arXiv:2110.03131v1 [astro-ph.HE] 7 Oct 2021 eago potnt o nesadn h uenv explo supernova the understanding for opportunity gravitatio good the s a only a not be detect such simultaneously from could wave one gravitational it, the of detection the galaxy, eeto rmtesproaepoinqiefrfo s Co us. from far quite explosion grav supernova of Sin the energy from radiation star. detection the massive symmetric, of spherically death almost the accompanies which supernova, hruhon conditions. thoroughgoing ∗ strong properties stellar local with determined frequency wave in gravitational kilohertz the of existence the us tell simulations num the and growing is distance observable the sensitivity, det of set limited a only if spatially poor, is detectors determining wave obse for The especially phenomena. astronomy, wave astronomical probing gravitational for tool a as waves rmgei ae n etio ihteosraino GW1 of [ observation merger the star with neutron neutrinos and waves tromagnetic ttesraeo h rtnurnsa PS,teso-called the (PNS), star protoneutron the of surface the at hc r h atpeei h opc iaymre vn,h event, merger binary compact the in fr piece waves last gravit gravitational the these the are O3, and which the O3, During LIG of 2022. by end August (O3) the after run at observing network third collaboration The this sky kilonova. the improve the to of possible location was detect it wave while gravitational observations, the netic by provided area location sky lcrncaddress: Electronic h etfote ngaiainlwv srnm fe the after astronomy wave gravitational in frontier next The h rvttoa ae rmsproaepoin aebee have explosions supernova from waves gravitational The ietdtcino rvttoa ae trsteeao gr of era the starts waves gravitational of detection Direct 4 etrfrCmuainlAtohsc,Ntoa Astronom National Astrophysics, Computational for Center 3 ∼ iiino cec,Ntoa srnmclOsraoyof Observatory Astronomical National Science, of Division 2 ASnmes 44.g 71.j 04.30.-w 97.10.Sj, 04.40.Dg, numbers: PACS observables. the affects t how state discuss of also equation and the models the in PNS of the dependence of independent the almost that is em confirm the we derive we addition, density, In average frequency. the t Using of function models. a PNS the as of expressed equatio be the accor can frequencies simulations and mode numerical mass gravity the progenitor and fundamental the the as that such find evolutio models, time PNS the the hand, on other (af the before On phase oscillations. (later) damental early the s in numerical oscillation the damental) in appearing Cow signals wave relativistic gravitational the the with asteroseismology (PNS) star tron nedsilnr hoeia ahmtclSinePro Science Mathematical & Theoretical Interdisciplinary 1 sn h ueia iuaindt o w-iesoa co two-dimensional for data simulation numerical the Using eodatrcrbuc.Ti apu inli rgnlyc originally is signal up ramp This corebounce. after second [email protected] 1 , 2 .Nw ebcm bet s o nyeetoantcwvsa waves electromagnetic only not use to able become we Now, ]. nvra eainfrsproagaiainlwaves gravitational supernova for relation Universal 1 aieSotani, Hajime srpyia i agLbrtr,RKN atm 351-019 Saitama RIKEN, Laboratory, Bang Big Astrophysical 5 eateto hsc,Thk nvriy edi980-8578 Sendai University, Tohoku Physics, of Department ,2, 1, ∗ ooaTakiwaki, Tomoya Dtd coe ,2021) 8, October (Dated: .INTRODUCTION I. inl,woefeunyicessfo e ude hertz hundred few a from increases frequency whose signals, ydpnso h ailpsto nietesa,tegravit the star, the inside position radial the on depends ly a ae u loteeetoantcwvsadneutrinos. and waves electromagnetic the also but waves nal clOsraoyo aa,22- sw,Mtk,Tko181- Tokyo Mitaka, Osawa, 2-21-1 Japan, of Observatory ical h orepsto,bcueteaglrrslto ihgr with resolution angular the because position, source the e feet sepce osgicnl increase. significantly to expected is events of ber ufc rvt ( gravity surface r culywr rca ofidn h an iooab ele by kilonova faint the finding to crucial were actually ors cosaeue.Frtecs fteG101,sace withi searches GW170817, the of case the For used. are ectors ttoa ae a o es tog hc ed oadifficul a to leads which strong, so be not may waves itational mtomre vnscmoe fanurnsa ihabakho black a with star neutron a of composed events merger two om oaino h rvttoa aeeet suigta tw it that assuming event, wave gravitational the of location se em ob ut ae eetees neoewudde would one once Nevertheless, rare. quite be to seems ystem vttoa aeatooyadad sanwmsegrt ele to messenger new a as adds and astronomy wave avitational inmcaim n ehv owi o uha potnt wi opportunity an such for wait to have we and mechanism, sion toa aedtcosaenwpeaigfrO,wihmyst may which O4, for preparing now are detectors wave ational v lobe eetd[ detected been also ave sdrn h vn aeo uenv,e.g, supernova, of rate event the nsidering vto feetoantcwvsi opeetr importa complementary is waves electromagnetic of rvation n ig a led nse,weeKGAfial joined finally KAGRA where finished, already has Virgo and O anysuisvanmrclsmltos(.. [ (e.g., simulations numerical via studies mainly n e oprdt h opc iaymres h uenv is supernova the mergers, binary compact the to compared ce, 01,wihi h rvttoa aeeetfo h binary the from event wave gravitational the is which 70817, opc iaymresms ete(aatc core-collaps (galactic) the be must mergers binary compact e)teaoddcosn ewe h rvt n fun- and gravity the between crossing avoided the ter) mlto a ewl dnie ihtegaiy(fun- gravity the with identified well be can imulation igapoiain ssoni h rvosstudy, previous the in shown As approximation. ling aa,22- sw,Mtk,Tko1188,Japan 181-8588, Tokyo Mitaka, Osawa, 2-21-1 Japan, fsproagaiainlwvssrnl depends strongly waves gravitational supernova of n rm(TES,RKN atm 5-18 Japan 351-0198, Saitama RIKEN, (iTHEMS), gram epooeto traeaedniyindependently density average star protoneutron he igt h rvttoa aesgasapaigin appearing signals wave gravitational the to ding eclas uenv,w xmn h protoneu- the examine we supernova, re-collapse N ufc est ntePSaeaedensity average PNS the on density surface PNS iia oml o uenv rvttoa wave gravitational supernova for formula pirical edfeetteteto h o-nfr matter non-uniform the of treatment different he ,4 3, fsaefrdnemte.Nvrhls,we Nevertheless, matter. dense for state of n n aieTogashi Hajime and g niee sarsl fteBrunt-Väisälä the of result frequency a as onsidered )md [ mode -) 3 5 .Oigt h eeomn fdetectors’ of development the to Owing ]. , 7 Japan , ,Japan 8, .Mawie ic h Brunt-Väisälä the since Meanwhile, ]. dnurnsbtas gravitational also but neutrinos nd 5 ∼ 1 e etr nour in century per 4 58 Japan 8588, – toa wave ational 19 pt few a to up avitational hsmust This ) These ]). ctromag- sa the at as the n yof ty tin nt tect art le, th c- e The maximum entropy distributions of collisionless particle velocity, speed, and energy for

statistical mechanics of self-gravitating collisionless flow (SG-CFD)

Zhijie Xu1,a

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

Northwest National Laboratory, Richland, WA 99352, USA

Abstract

The halo-mediated inverse mass cascade is a key feature of the intermediate statistically steady

state for self-gravitating collisionless flow (SG-CFD). How the inverse mass cascade maximizes

the system entropy and develops limiting velocity/energy distributions are fundamental questions

to answer. We present a statistical theory concerning the maximum entropy distributions of particle

velocity, speed, and energy for self-gravitating systems involving a power-law long-range

interaction with an arbitrary exponent n. For system with long-range interaction (-2

spectrum of halos and halo groups are necessary to form from inverse mass cascade to maximize

the system entropy. While particle velocity in each halo group is still Gaussian, the velocity

distribution of entire system can be non-Gaussian. With the virial equilibrium for local mechanical

equilibrium of halo groups, the maximum entropy principle is applied for statistical equilibrium of global system to derive the limiting distributions. Halo mass function is not required in this

formulation, but it is a direct result of entropy maximization. The predicted velocity distribution involves a shape parameter α that is dependent on the exponent n. The velocity distribution approaches Laplacian with   0 and Gaussian with  . For intermediate α, the maximum

entropy distribution naturally exhibits a Gaussian core at small velocity and exponential wings at

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

Dodgeball—Can a Satellite Avoid Being Hit by Debris?

J. I. Katz Department of Physics and McDonnell Center for the Space Sciences Washington University, St. Louis, Mo. 63130 October 8, 2021

Abstract Can a satellite dodge a collision with untracked orbiting debris? Can a satellite dodge collision with a tracked object, making only the avoidance manœuvers actually required to avoid collision, despite the uncertainties of predicted conjunctions? Satellite-borne radar may dis- tinguish actual collision threats from the much greater number of near misses because an object on a collision course has constant bearing, which may be determined by interferometric detection of the radar return. A large constellation of such radars may enable the determination of the ephemerides of all cm-sized debris in LEO.

1 Introduction

Large constellations of satellites are planned to be launched in the coming years, many to provide world-wide internet access, with cumulative totals of 100,000 or more. As the number of satellites and orbital debris increases, arXiv:2110.03113v1 [astro-ph.EP] 6 Oct 2021 so does the threat collisions with other orbiting objects. Collisions not only damage or destroy the colliding objects, but also create multiple pieces of debris, increasing future hazards. How may a satellite may manœuver to avoid an imminent collision? Pre- dicted conjunctions, approaches closer than the orbital uncertainties risk collision, are comparatively frequent because these uncertainties are signiﬁcant.

1 Constraining cosmological phase transitions with the Parkes Pulsar Timing Array

Xiao Xue,1, 2, 3 Ligong Bian,4, 5, ∗ Jing Shu,1, 2, 6, 7, 8, 9, † Qiang Yuan,10, 11, 8, ‡ Xingjiang Zhu,12, 13, 14, § N. D. Ramesh Bhat,15 Shi Dai,16 Yi Feng,17 Boris Goncharov,12, 13 George Hobbs,18 Eric Howard,18, 19 Richard N. Manchester,18 Christopher J. Russell,20 Daniel J. Reardon,13, 21 Ryan M. Shannon,13, 21 Renee´ Spiewak,22, 21 Nithyanandan Thyagarajan,23 and Jingbo Wang24 1CAS Key Laboratory of Theoretical Physics, Insitute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 2School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China 3II. Institute of Theoretical Physics, UniversitätHamburg, 22761 Hamburg, Germany 4 Department of Physics, Chongqing University, Chongqing 401331, China 5 Chongqing Key Laboratory for Strongly Coupled Physics, Chongqing 401331, China 6CAS Center for Excellence in Particle Physics, Beijing 100049, China 7School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310024, China 8Center for High Energy Physics, Peking University, Beijing 100871, China 9International Center for Theoretical Physics Asia-Pacific, Beijing/Hanzhou, China 10Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023, China 11School of Astronomy and Space Science, University of Science and Technology of China, Hefei 230026, China 12School of Physics and Astronomy, Monash University, Clayton, VIC 3800, Australia 13OzGrav: The ARC Centre of Excellence for Gravitational Wave Discovery, Hawthorn, VIC 3122, Australia 14Advanced Institute of Natural Sciences, Beijing Normal University at Zhuhai 519087, China 15International Centre for Radio Astronomy Research, Curtin University, Bentley, WA 6102, Australia 16Western Sydney University, Locked Bag 1797, Penrith South DC, NSW 1797, Australia 17University of Chinese Academy of Sciences, Beijing 100049, China 18CSIRO Astronomy and Space Science, P.O. Box 76, Epping, NSW 1710, Australia 19Macquarie University, Department of Physics and Astronomy, Sydney, NSW, 2109, Australia 20CSIRO Scientific Computing, Australian Technology Park, Locked Bag 9013, Alexandria, NSW 1435, Australia 21Centre for Astrophysics and Supercomputing, Swinburne University of Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia 22Jodrell Bank Centre for Astrophysics, University of Manchester, Manchester M13 9PL, UK 23CSIRO Astronomy and Space Science (CASS), P. O. Box 1130, Bentley, WA 6102, Australia 24Xinjiang Astronomical Observatory, Chinese Academy of Sciences, 150 Science 1-Street, Urumqi, Xinjiang 830011, China A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature (T ∼ 1 − 100 MeV) −1 −2 −1 phase transition with a duration (β/H∗) ∼ 10 − 10 and therefore can be used to constrain the dark and QCD phase transitions.

PACS numbers:

I. INTRODUCTION

arXiv:2110.03096v1 [astro-ph.CO] 6 Oct 2021 expected to form during cosmological phase transitions, such as cosmic strings. Since the discovery of gravitational waves in 2015, the study of SGWB has gained increasingly broad in- The stochastic gravitational wave background (SGWB) is terest because of its close connection with particle physics and an important target of gravitational-wave astronomy after the cosmology. Different gravitational wave experiments, such as successful detection of gravitational waves from 50 compact LIGO and the upcoming Laser Interferometer Space Antenna binary mergers by the LIGO and Virgo detectors [1–3]. The (LISA), aim at the detection of SGWB at different frequency SGWB may come from the superposition of numerous un- bands [4–6]. resolvable binary coalescences, or from the cosmological in- The pulsar timing array (PTA) experiments measure the ar- flation, cosmological phase transition, or topological defects rival times of radio pulses from a number of highly stable millisecond pulsars in the Milky Way, which enables the detection of gravitational waves of nanohertz frequencies. The three major PTA experiments, the Parkes PTA (PPTA; [7]), the ∗Electronic address: [email protected] †Electronic address: [email protected] European PTA [8], and NANOGrav [9] have been monitoring ‡Electronic address: [email protected] the times of arrival (ToAs) of dozens of pulsars on weekly to §Electronic address: [email protected] monthly cadence for over 10 years, and provide unique probes Received xx xxxx xxxx; Revised xx xxxx xxxx; Accepted xx xxxx xxxx

DOI: xxx/xxxx

ARTICLE TYPE The gaseous natal environments of GPS and CSS sources with ASKAP – FLASH

James R. Allison1,2 | Elaine M. Sadler2,3,4 | Elizabeth K. Mahony4 | Vanessa A. Moss4,3 | Hyein Yoon2,3

1Sub-Dept. of Astrophysics, Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Rd., Oxford, GPS and CSS sources are thought to represent a young and/or conﬁned sub- OX1 3RH, U.K. population of radio-loud active galactic nuclei (AGN) that are yet to evacuate their 2ARC Centre of Excellence for All-Sky Astrophysics in 3 Dimensions (ASTRO 3D) surrounding ambient interstellar gas. By studying the gaseous environments of these 3Sydney Institute for Astronomy, School of objects we can gain an insight into the inter-dependent relationship between galaxies Physics A28, University of Sydney, NSW and their supermassive black holes (SMBHs). The First Large Absorption Survey in 2006, Australia 4Space & Astronomy, CSIRO, PO Box 76, H I (FLASH) will build a census of the neutral atomic hydrogen (H I) gas in galaxies Epping, NSW 1710, Australia at intermediate cosmological redshifts. FLASH is expected to detect at least several

Correspondence hundred H I absorbers associated with GPS and CSS sources. These absorbers pro- *James R. Allison, Email: vide an important probe of the abundance and kinematics of line-of-sight neutral [email protected] gas towards radio AGN, in some cases revealing gas associated with infalling clouds and outflows. Observations are now complete for the first phase of the FLASH Pilot Survey and early analysis has already yielded several detections, including the GPS source PKS 2311−477. Optical imaging of this galaxy reveals an interacting system that could have supplied the neutral gas seen in absorption and triggered the radio-loud AGN. FLASH will provide a statistically significant sample with which the prevalence of such gas-rich interactions amongst compact radio galaxies can be investigated.

KEYWORDS: galaxies:active, radio continuum:galaxies, radio lines:galaxies, galaxies:ISM, galaxies:individual (PKS 2311−477)

1 INTRODUCTION radio spectrum, is a particularly useful method of studying the arXiv:2110.03046v1 [astro-ph.GA] 6 Oct 2021 natal gaseous environments of young radio galaxies. Detec- The GHz-Peaked Spectrum (GPS) and Compact Steep Spec- tion rates of H I absorption amongst GPS and CSS sources trum (CSS) radio sources are thought to be a young and/or tend to be relatively high (∼ 30 per cent) and provide line-of- confined sub-population, with linear sizes less than an about sight kinematic information on the cool neutral gas towards 20 kpc (O’Dea & Saikia 2021). As such, they provide excellent the source (see Morganti & Oosterloo 2018 for a review). objects with which to study the interaction between radio- However, the historical number of compact sources searched loud active galactic nuclei (AGN) and their host galaxies. In for H I absorption remains insufficient to draw strong conclu- some cases multi-phase gaseous outflows are detected, possi- sions about the population, particularly at cosmological dis- bly pointing to a form of radio-jet feedback in the early stages tances. Larger spectroscopic radio surveys with the new Square of radio AGN evolution (e.g. Holt, Tadhunter, & Morganti Kilometre Array (SKA) pathfinder telescopes are therefore 2008). The H I 21-cm line, when detected as absorption in the required to build a statistically significant sample. Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX62

Constraints on the preferential disruption of moderately massive stars by supermassive black holes

Brenna Mockler,1 Angela A. Twum,1 Katie Auchettl,2, 3, 1, 4 Sierra Dodd,1 K.D. French,5 Jamie A.P. Law-Smith,6 and Enrico Ramirez-Ruiz1

1Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA 2School of Physics, The University of Melbourne, VIC 3010, Australia 3ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) 4OzGrav, University of Melbourne, Parkville, Victoria 3010, Australia 5Department of Astronomy, University of Illinois, 1002 W. Green St., Urbana, IL 61801, USA 6Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., Cambridge, MA 02138, USA

Submitted to ApJ

ABSTRACT Tidal disruption events (TDEs) provide a unique opportunity to probe the stellar populations around supermassive black holes (SMBHs). By combining light curve modeling with spectral line information and knowledge about the stellar populations in the host galaxies, we are able to constrain the properties of the disrupted star for a handful of TDEs. The TDEs in our sample have UV spectra, and measurements of the UV N III to C III line ratios enabled estimates of the nitrogen-to-carbon abundance ratios for these events. We show that the measured nitrogen line widths are consistent with originating from the disrupted stellar material dispersed by the central SMBH. We ﬁnd that these nitrogen-to-carbon abundance ratios necessitate the disruption of moderately massive stars (& 1 2M ). We determine − 2 that these moderately massive disruptions are over-represented by a factor of & 10 when compared to the overall stellar population of the post-starburst galaxy hosts. This implies that SMBHs are preferentially disrupting higher mass stars, possibly due to ongoing top-heavy star formation in nuclear star clusters or to dynamical mechanisms that preferentially transport higher mass stars to their tidal radii.

Keywords: accretion, accretion disks – black hole physics – galaxies: nuclei – stars: individual (iPTF15af, iPTF16fnl, ASASSN-14li, ASASSN-18pg, AT2018dyb), Tidal Disruption Events, SMBHs, Spectroscopy

1. INTRODUCTION the internal structure of the star has to be taken into The gravitational influence of a super-massive black consideration and we must switch to a hydrodynamic hole dominates the kinematics of stars within a nuclear description to follow the evolution of the disrupted stel- cluster. Each star within this dense region traces out lar debris. A fraction of the debris eventually falls back, arXiv:2110.03013v1 [astro-ph.HE] 6 Oct 2021 an intricate trajectory under the combined influence of circularizes, and accretes onto the SMBH. This accretion the SMBH and the other stars. Stellar dynamics is ca- powers a flare which is a clear sign of the presence of an pable of describing the resulting stellar motions as point otherwise quiescent SMBH. As we show here, it can also masses moving under the influence of gravity until these provide a compelling fingerprint of the properties of the encounters move a star onto a nearly radial orbit. If a disrupted star. star wanders close to the tidal radius it is ripped apart by It has only been in the last decade – with the advent the SMBH’s gravity (Rees 1988). When that happens, of numerous wide field transient surveys – that we have started to collect photometric and spectroscopic data on a myriad of tidal disruption events (TDEs; e.g., van [email protected] Velzen et al. 2020, and references therein). The discov- Corresponding author: Brenna Mockler ery of TDE flares has generated widespread excitement, October 8, 2021 0:34 ws-procs961x669 WSPC Proceedings - 9.61in x 6.69in ws-procs961x669 page 1

Testing extended theories of gravity with GRBs

L. Mastrototaro, G. Lambiase Dipartimento di Fisica “E.R Caianiello”, Universit`adegli Studi di Salerno, Via Giovanni Paolo II, 132 - 84084 Fisciano (SA), Italy Istituto Nazionale di Fisica Nucleare - Gruppo Collegato di Salerno, Via Giovanni Paolo II, 132 - 84084 Fisciano (SA), Italy. ∗E-mail: [email protected]

We present our studies on the neutrino pairs annihilation into electron-positron pairs (νν¯ → e−e+) near the surface of a neutron star in the framework of extended theories of gravity. The latter modifies the maximum energy deposition rate near to the pho- tonsphere and it might be several orders of magnitude greater than that computed in the framework of General Relativity. These results provide a rising in the Gamma-Ray Bursts energy emitted from a close binary neutron star system and might be a fingerprint of modified theories of gravity, changing our view of astrophysical phenomena.

Keywords: Extended theories of gravity, GRB, Neutrino energy deposition, Black Holes, Neutron Stars

1. Introduction General Relativity (GR) is without any doubt the best theory of gravitational interaction. Its predictions have been tested with high accuracy on scales of the solar system (for example, the precession of the Mercury perihelion, the photons deviation and the gravitational leasing effect), on astrophysical scales (the gravitational waves), and cosmological scales (the cosmic microwave background radiation (CMBR) and the formation of primordial light elements (Big Bang Nucleosynthe- sis)). Despite these results, there are still open questions that make GR incomplete. The latter arise at short distances and small time scales (black hole and cosmological singularities, respectively), or at large distance scales, the rotational curve of the galaxies and the observed accelerated phase of the present Universe, for which any predictability is lost. To solve these issues, deviations from the GR (hence from the Hilbert-Einstein action on which GR is based) are needed, and new ingredients, such as dark matter and dark energy, are required for fitting the present picture of our Universe 1–6. Indeed, in the last years, several alternative or modified theories of gravity have been proposed, which try to answer all the opened questions of GR and the Cosmological Standard Model. To give an example, higher-order curvature invariants than the arXiv:2110.02996v1 [astro-ph.HE] 6 Oct 2021 simple Ricci scalar R, allow getting inflationary behaviour, removing the primordial singularity, as well as explaining the flatness and horizon problems 7,8( for further applications, see Refs. 9–37). On the other hand, one can tend to preserve the GR and the Hilbert-Einstein action, adding only some two unidentified components: dark matter and dark energy. These two different approaches try to solve the same problems but at the moment there is not a final solution for that (modified gravity Astronomy & Astrophysics manuscript no. aa ©ESO 2021 October 8, 2021

Modeling the Mg i from the NUV to MIR: I. The Solar Case J.I. Peralta1, 2, Vieytes M.C.1, 2, A.M.P. Mendez1, and D.M. Mitnik1, 3

1 Instituto de Astronomía y Física del Espacio, CONICET–Universidad de Buenos Aires, Argentina e-mail: [email protected]; [email protected] 2 Departamento de Ciencia y Tecnología, UNTREF, Argentina 3 Departamento de Física, Universidad de Buenos Aires, Argentina

Received 6 August 2021, Accepted 6 October 2021

ABSTRACT

Context. Semi-empirical models of the solar atmosphere are used to study the radiative environment of any planet in our solar system. There is a need for reliable atomic data for neutral atoms and ions in the atmosphere to obtain improved calculated spectra. Atomic parameters are crucial to computing the correct population of elements through the whole stellar atmosphere. Although there is a very good agreement between the observed and calculated spectra for the Sun, there is still a mismatch in several spectral ranges due to the lack of atomic data and its inaccuracies, particularly for neutrals like Mg i. Aims. To correctly represent many spectral lines of Mg i from the near-UV to the mid-IR is necessary to add and update the atomic data involved in the atomic processes that drive their formation. Methods. The improvements to the Mg i atomic model are as follows: i) 127 strong lines, including their broadening data, were added. ii) To obtain these new lines, we increased from 26 to 85 the number of energy levels. iii) Photoionization cross-section parameters were added and updated. iv) Effective Collision Strengths (Υi j) parameters were updated for the first 25 levels using the existing data from the convergent close-coupling (CCC) calculations. One of the most significant changes in our model is given by the new Υi j parameters for transitions involving levels between 26 and 54, which were computed with a multi–configuration Breit–Pauli distorted–wave (DW) method. For transitions involving superlevels, we calculated the Υi j parameters with the usual semi-empirical formulas of Seaton(1962) and van Regemorter(1962). Results. More than one hundred transitions were added to our calculations, increasing our capability of reproducing important features observed in the solar spectra. We found a remarkable improvement in matching the solar spectra for wavelengths higher than 30 000 Å when our new DW Υi j data was used in the model. Conclusions. Key words. atomic:data – sun:atmospheres – Sun:infrared – line:formation – line:profile

1. Introduction Neutral magnesium is an important chemical element since its lines are strong in the spectra of late-type and even in metal- 1D semi-empirical model calculation of the solar atmosphere poor stars, making it a good tracer of α-element abundances. At continues to be useful to predict and reproduce the Solar Spec- chromospheric temperatures, magnesium is susceptible to devi- tral Irradiance (SSI), from the far-infrared (FIR) to the extreme ations from local thermodynamic equilibrium (LTE). These de- ultraviolet (EUV), received at the top of the Earth’s atmosphere viations from LTE are predicted to be significant, particularly (Fontenla & Landi 2018). Changes in the SSI can produce mod- for metal-poor stars (Zhao et al. 1998; Zhao & Gehren 2000). ifications in the chemistry of our atmosphere. Therefore to study the non-LTE (NLTE) effects it is necessary to account for accurate atomic data. The advantage of these models is that they are built to reproduce a large number of observations in different wavelengths. Past studies have used semi-empirical formulae when com- The equations of statistical equilibrium and radiative transfer are puting electron atomic collision data. The most widely used treated in detail and solved simultaneously for all the species in semi-classical formulas for radiative allowed transitions are the the atmosphere. One of the most critical parameters is the data ones given by Seaton(1962) and van Regemorter(1962) since arXiv:2110.02992v1 [astro-ph.SR] 6 Oct 2021 describing the atomic and molecular species present in the at- they provide estimates for the collision rates based on the oscilla- mosphere. Great efforts have been made to obtain experimental tor strengths and transition energies. However, these approxima- and theoretical data. However, there is still a lack of complete tions can lead to a great number of uncertainties in the synthetic and accurate atomic data, particularly for neutral and low ion- spectra of NLTE atmospheric models, such as spectral lines in ized atoms. These atomic species, such as Mg i, are of particu- emission instead of in absorption as observations reveal. lar interest because they have relevant spectral features used in There are several studies on NLTE calculations of Mg i in the the diagnosis of solar and stellar astrophysics. The influence of Sun and other stars, for instance, the recent works by Alexeeva using a more complete and accurate atomic model for neutral et al.(2018), the paper series by Osorio et al.(2015), Osorio atoms in SSI calculations can be found, for example, in Vieytes & Barklem(2016) and Barklem et al.(2017), and the relevant & Fontenla(2013), were the improvements made in the Ni i works by Bergemann et al.(2015) and Scott et al.(2015). In atomic model led to a better match with observations, mainly all of these studies, the authors stress the need for accurate Mg i in the near-ultraviolet (NUV) region. atomic data.

Article number, page 1 of 14 Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX631

HIFlow: Generating Diverse HI Maps Conditioned on Cosmology using Normalizing Flow

Sultan Hassan ,1, 2 Francisco Villaescusa-Navarro ,3, 1 Benjamin Wandelt,4, 1 David N. Spergel,1, 3 Daniel Angles-Alc´ azar´ ,5, 1 Shy Genel,1, 6 Miles Cranmer,3 Greg L. Bryan,7, 1 Romeel Dave,´ 8, 2, 9 Rachel S. Somerville,1 Michael Eickenberg,10 Desika Narayanan,11, 12 Shirley Ho,1, 3, 13 and Sambatra Andrianomena14, 2

1Center for Computational Astrophysics, Flatiron Institute, 162 5th Ave, New York, NY, 10010, USA 2Department of Physics & Astronomy, University of the Western Cape, Cape Town 7535, South Africa 3Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ, 08544, USA 4Sorbonne Universite, CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis boulevard Arago, 75014 Paris, France 5Department of Physics, University of Connecticut, 196 Auditorium Road, U-3046, Storrs, CT, 06269, USA 6Columbia Astrophysics Laboratory, Columbia University, New York, NY, 10027, USA 7 Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 10027, USA 8 Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 9 South African Astronomical Observatories, Observatory, Cape Town 7925, South Africa 10Center for Computational Mathematics, Flatiron Institute, 162 5th Ave, New York, NY, 10010, USA 11Department of Astronomy, University of Florida, Gainesville, FL, USA 12University of Florida Informatics Institute, 432 Newell Drive, CISE Bldg E251, Gainesville, FL, USA 13Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA 14 South African Radio Astronomy Observatory (SARAO), Black River Park, Observatory, Cape Town, 7925, South Africa

ABSTRACT A wealth of cosmological and astrophysical information is expected from many ongoing and upcoming large-scale surveys. It is crucial to prepare for these surveys now and develop tools that can efficiently extract the maximum amount of information. We present HIFlow: a fast emulator that is able to generate neutral hydrogen (HI) maps conditioned only on cosmology (Ωm and σ8), after training on the state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. HIFlow is designed using a class of normalizing flow models, the Masked Autoregressive Flow (MAF), which we demonstrate are capable of generating realistic maps without explicitly using the 2D structure or accounting for any symmetries. HIFlow is able to generate 14 21 −2 new diverse HI maps in the column density range NHI ∼ 10 − 10 cm at z ∼ 6, and naturally mimic the cosmic variance effects. Remarkably, HIFlow is able to reproduce the CAMELS average 2 and standard deviation HI power spectrum (Pk) within a factor of . 2, scoring a very high R > 90%. HIFlow will enable the testing of Pk pipelines for HI surveys, and assist in computing other statistical properties beyond Pk that require generating new diverse samples of high dimensional datasets, such as the covariance matrix. This new tool represents a first step towards enabling rapid parameter inference at the field level, maximizing the scientific return of future HI surveys, and opening a new avenue to minimize the loss of information due to data compression. arXiv:2110.02983v1 [astro-ph.CO] 6 Oct 2021

Keywords: early universe – dark ages, reionization, ﬁrst stars, methods: statistical

1. INTRODUCTION coming large-scale surveys such as the Square Kilo- Extracting the maximum amount of cosmological and meter Array (SKA, Mellema et al. 2013), the Hy- astrophysical information remains a challenge in up- drogen Epoch of Reionization Array (HERA, DeBoer et al. 2017), the Low Frequency Array (LOFAR, van Haarlem et al. 2013), the Vera C. Rubin Observatory Corresponding author: Sultan Hassan Legacy Survey of Space and Time (LSST, Ivezi´cet al. shassan@ﬂatironinstitute.org 2019), Nancy Grace Roman Space Telescope (Roman, Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX631

Investigating the Accretion Nature of Binary Supermassive Black Hole Candidate SDSS J025214.67-002813.7 Adi Foord,1 Xin Liu,2, 3 Kayhan Gultekin¨ ,4 Kevin Whitley,4 Fangzheng Shi,5, 6 and Yu-Ching Chen2, 3, 7

1Kavli Institute of Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA 2Department of Astronomy, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA 3National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA 4Department of Astronomy and Astrophysics, University of Michigan, Ann Arbor, MI 48109 5School of Astronomy and Space Science, Nanjing University, Nanjing 210023, China 6Key Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Nanjing 210023, China 7Center for AstroPhysical Surveys, National Center for Supercomputing Applications, Urbana, IL, 61801, USA

ABSTRACT We present results on a multi-wavelength analysis of SDSS J025214.67-002813.7, a system which has been previously classified as a binary AGN candidate based on periodic signals detected in the optical light curves. We use available radio X-ray observations of the system to investigate the true − accretion nature. Analyzing new observations from XMM-Newton and NuSTAR, we characterize the X-ray emission and search for evidence of circumbinary accretion. Although the 0.5–10 keV spectrum shows evidence of an additional soft emission component, possibly due to extended emission from hot nuclear gas, we find the spectral shape consistent with a single AGN. Compiling a full multi-wavelength SED, we also search for signs of circumbinary accretion, such as a “notch” in the continuum due to the presence of minidisks. We find that the radio–optical emission agrees with the SED of a standard, radio-quiet, AGN, however there is a large deficit in emission blueward of 1400 A.˚ Although this ∼ deficit in emission can plausibly be attributed to a binary AGN system, we find that the SED of SDSS J0252 0028 is better explained by emission from a reddened, single AGN. However, future studies on − the expected hard X-ray emission associated with binary AGN (especially in the unequal-mass regime), will allow for more rigorous analyses of the binary AGN hypothesis.

Keywords: X-ray active galactic nuclei (2035) — Supermassive black holes (1663) — Accretion (14) — Galaxy mergers (608)

1. INTRODUCTION wave background signal detected by pulsar timing ar- A binary supermassive black hole (SMBH) represents rays (PTAs; Burke-Spolaor et al. 2019), and they are the final stage of a galaxy merger, where the two massive direct precursors to gravitational-wave events detectable host galaxies have likely been interacting for hundreds by future space-based laser interferometers (Sesana et al. of megayears to gigayears (Begelman et al. 1980). The 2007). merging system is classified as a binary when the SMBHs Currently, detections of binary SMBHs are limited to systems where both black holes are actively accret- arXiv:2110.02982v1 [astro-ph.HE] 6 Oct 2021 are gravitationally bound in a Keplerian orbit, and for a wide range SMBH masses and host galaxy environments ing (active galactic nuclei, or AGN) and emitting light this occurs at orbital separations <10 pc (Mayer et al. across the electromagnetic spectrum. Given the small 2007; Dotti et al. 2007; Khan et al. 2012). The fate and physical separations between the two SMBHs, the an- final coalescence of the system strongly depends on the gular resolution afforded by radio interferomerty is re- amount of matter the SMBHs can interact with (Sesana quired to spatially resolve binary AGN systems, such as et al. 2007; Merritt et al. 2007). As the last stage before in the serendiptious discovery of 0402+379 (Rodriguez coalescence, binary SMBHs represent an observable link et al. 2006). However, radio detections of binary AGN between galaxy mergers and gravitational wave events. remain limited to relatively low-redshifts (where separa- They are strong emitters of low-frequency gravitational tions below 1 pc can only be probed up to z = 0.1) and systems with two radio-bright AGN (where only 15% waves, which are expected to dominate the gravitational ∼ arXiv:2110.02979v1 [astro-ph.SR] 6 Oct 2021 [email protected] rmtedul lc oe(H egrevent merger (BH) first hole signals black the (GW) ( double GW150914 detected wave the from gravitational (LIGO) high-frequency Observatory Gravitational- Interferometer Wave Laser advanced the when ( ae,mrigtebgnigo e r fmulti- of era new a entire of ( radio astronomy the beginning messenger to in the also gamma-rays marking but waves, from GWs, spectrum in detected electromagnetic GW170817 only event not merger was (DNS) NS double remarkable ( other 2021 ( h IOSinicClaoaine l 2021 al. et Collaboration Scientific LIGO The 2017a al. et Abbott Collaboration Virgo & Collaboration Scientific LIGO e r fatohsc ea i er ago years six began astrophysics of era new A yee sn L using 2021 Typeset 8, October version Draft ,a ela opeo obenurnsa (NS) star neutron double of couple a as well as ), ∼ of rnrc oeof core iron-rich h nalvlct fteclasn rncr xed 00k s km 1000 exceeds core iron collapsing Ca the via of modelling velocity the infall continue the and star helium main-sequence zero-age Keywords: repro to anticipated is recyc systems mild of LIGO-network. the family the as this well by as conditions, system, initial post-SN the the on of fate the discuss we h usin fteroii n oncin epeettefis 1 first the ultra-stripped present an We until self-consistently connection. MESA calculated and is origin which their system of questions the eetdsoeiso obenurnsa DS egr n ultr and mergers (DNS) star neutron double of discoveries Recent 0dul Hmreshv endetected been have mergers BH double 50 ∼ 2 oe oe fa lr-tipdSproaPoeio faD a of Progenitor Supernova Ultra-stripped an of Model Novel eateto hsc n srnm,Aru nvriy Ny University, Aarhus Astronomy, and Physics of Department 1 bote l 2016 al. et Abbott oesatn rmaps-omnevlp iaycnitn fa1 a of consisting binary envelope post-common a from starting code . 90 3 A colo hsc n lcrclIfrain hnquNor Shangqiu Information, Electrical and Physics of School T 1. E M X -a iaysas(81,Selreouin(1599) evolution Stellar (1811), stars binary X-ray neatn iaysas(0) eto tr 10) Gravita (1108), stars Neutron (801), stars binary Interacting 4 ⊙ INTRODUCTION ogJiang, Long twocolumn , AI,Micd 5-7 aionaIsiueo Technolog of Institute California 350-17, Mailcode TAPIR, ossigo a of consisting , 2020 bote l 2017b al. et Abbott 1 ∼ colo cec,QndoUiest fTcnlg,Qingd Technology, of University Qingdao Science, of School n ie HN mergers BH/NS mixed and ) 1 . 50 M tl nAASTeX63 in style .Sbeunl,an- Subsequently, ). ,2 3 2, 1, ⊙ h eutn eodbr Shsa siae asof mass estimated an has NS second-born resulting The . Rcie coe ,22;Rvsd Accepted) Revised; 2021; 8, October (Received hmsM Tauris, M. Thomas ∼ 0 . 29 ). M ⊙ .The ). erc neoeebdigaC oeof core CO a embedding envelope He-rich ABSTRACT 2 e-ogChen, Wen-Cong egr a erltdt Rsvapouto fsta- of production via FRBs to related be may mergers pro- besides Hence, short mergers). a BH ducing double companion. (unlike mat- baryonic its involve ter DNSs and of events pulsar merger the GW Third, of masses the of determine theories for tests (e.g. important gravity provide parame- several thereby post-Keplerian and of relativistic ters timing measure pulsar can DNSs radio high-precision Second, 2003 2002 al. et well Dewi ( as explosions phases, bina- (SN) (CE) supernova X-ray two envelope in as common transfer more or mass one sur- of ries, stages via: multiple history of evolutionary vival past key their provide and of evolution information stellar important binary are exotic they of First, probes astrophysics. modern of pects ea n oe notcltaset(..a (i.e. transient optical e.g. an power elements and chemical which nucleosynthesis, decay heavy r-process through of iron), production (beyond to lead events N ytm lyavtlrl nabodrneo asof range broad a in role vital a play systems DNS ukgd 2,D-00Aru ,Denmark C, Aarhus DK-8000 120, Munkegade mrte l 2017 al. et Smartt ; a nvriy hnqu460,China 476000, Shangqiu University, mal osalwk ta.2004 al. et Podsiadlowski − 1 ,Psdn,C 12,USA 91125, CA Pasadena, y, h xldn trhsattlmass total a has star exploding The . igo h rtbr S Depending NS. first-born the of ling rmre l 2006 al. et Kramer -tipdsproa Se raise (SNe) supernovae a-stripped rncr olpe eapythe apply We collapse. core iron ,3 1, o262,China 266525, ao γ eB oh-oeoeflwuntil overflow Roche-lobe BB se ueteDSmresdetected mergers DNS the duce ; oe faDSprogenitor DNS a of model D rybrt( burst -ray n i Fuller Jim and ecysie l 2002 al. et Belczynski . inlwv ore (677), sources wave tional 35 ul eto Star Neutron ouble ; M eze 2019 Metzger ⊙ Sada3 a and NS bote l 2017b al. et Abbott ∼ 1 ∼ swl saccurately as well as ) . ; 4 61 a e evl1976 Heuvel den van 1 arse l 2017 al. et Tauris . M 44 ⊙ .Fnly DNS Finally, ). M ; . n an and 20 os&Tauris & Voss ⊙ M and ⊙ kilonova these ) ). , ; MNRAS 000,1–10 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

Protoplanet Collisions: Statistical Properties of Ejecta

Samuele Crespi1,2?, Ian Dobbs-Dixon1,2,3, Nikolaos Georgakarakos1,2, Nader Haghighipour4,5, Thomas I. Maindl6,7, Christoph M. Sch¨afer8, Philip Matthias Winter9.

1New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates 2Center for Astro, Particle and Planetary Physics (CAP3), New York University Abu Dhabi 3Center for Space Science, New York University Abu Dhabi 4Planetary Science Institute, 1700 East Fort Lowell, Tucson, AZ 85719, USA 5Institute for Astronomy, University of Hawaii-Manoa, Honolulu, HI 96822, USA 6Department of Astrophysics, University of Vienna, 1180 Vienna, Austria 7SDB Science-driven Business Ltd, 6025 Larnaca, Cyprus 8Institut fur¨ Astronomie und Astrophysik, Eberhard Karls Universit¨at Tubingen,¨ Auf der Morgenstelle 10, 72076 Tubingen,¨ Germany 9Institute for Machine Learning, Johannes Kepler University Linz, Altenberger Straße 69, 4040 Linz, Austria

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT

The last phase of the formation of rocky planets is dominated by collisions among Moon- to Mars-sized planetary embryos. Simulations of this phase need to handle the difficulty of including the post-impact material without sat- urating the numerical integrator. A common approach is to include the collision-generated material by clustering it into few bodies with the same mass and uniformly scattering them around the collision point. However, this approach oversimplifies the properties of the collision material by neglecting features that can play important roles in the final structure and composition of the system. In this study, we present a statistical analysis of the orbital architecture, mass, and size distributions of the material generated through embryo-embryo collisions and show how they can be used to develop a model that can be directly incorporated into the numerical integrations. For instance, results of our analysis indicate that the masses of the fragments follow an exponential distribution with an exponent of −2.21±0.17 over the range of 10−7 to 2 × 10−2 Earth-masses. The distribution of the post-impact velocities show that a large number of fragments are scattered toward the central star. The latter is a new finding that may be quite relevant to the delivery of material from the outer regions of the asteroid belt to the accretion zones of terrestrial planets. Finally, we present an analytical model for the 2D distribution of fragments that can be directly incorporated into numerical integrations. Key words: planets and satellites: formation – planets and satellites: terrestrial planets

1 INTRODUCTION As the time of N-body integrations increases with the num- arXiv:2110.02977v1 [astro-ph.EP] 6 Oct 2021 ber of bodies, in order to ensure that formation simulations Based on the standard model for planet formation, the last are carried out in a reasonable amount of time, collisions are phase of rocky planet formation is dominated by collisions often treated as perfectly inelastic. That is, the two collid- between Moon- to Mars-sized solid bodies, known as proto- ing bodies merge completely and no fragments and debris planets or planetary embryos. To a first approximation, this are produced (e.g., Raymond et al. 2004, O’Brien et al. 2006 phase can be simulated as the evolution of a system of ob- Morishima 2015). This approach has been proven to be use- jects that interact exclusively via gravity (e.g., Ida & Makino ful in demonstrating the formation process, especially as a 1993, Raymond et al. 2006, Barnes et al. 2009, Morishima proof of concept. However, it tends to overestimate the col- 2015, Clement et al. 2019). These simulations are generally lision efficiency, leading to shorter formation time and more carried out using N-body integrators with the largest con- massive planets. The latter causes the final planetary system straint being the number of bodies, usually kept below a few to have smaller number of bodies with planets being farther hundred to avoid lengthy computations. away from the central star (Dugaro et al. 2020, Burger et al. 2020). While simulating collisions has been an integral part of ? E-mail: [email protected] planetary science for decades (e.g., Dormand & Woolfson

A systematic study of the escape of LyC and Lyα photons from star-forming, magnetized turbulent clouds

Taysun Kimm ,1 Rebekka Bieri ,2 Sam Geen ,3 Joakim Rosdahl ,4 Jer´ emy´ Blaizot ,4 Leo´ Michel-Dansac,4 and Thibault Garel 5

1Department of Astronomy, Yonsei University, 50 Yonsei-ro, Seodaemun-gu, Seoul 03722, Republic of Korea 2Max Planck Institute for Astrophysics, Karl-Schwarzschild-Straße 1, 85748, Garching, Germany 3Anton Pannekoek Institute for Astronomy, Universiteit van Amsterdam, Science Park 904, NL-1098 XH Amsterdam, the Netherlands 4Univ Lyon, Univ Lyon1, Ens de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F-69230 Saint-Genis-Laval, France 5Observatoire de Genevé,Universitède Genevè,51 Ch. des Maillettes, 1290 Versoix, Switzerland

(Received October 8, 2021)

ABSTRACT Understanding the escape of Lyman continuum (LyC) and Lyman α (Lyα) photons from giant molecular clouds (GMCs) is crucial if we are to study the reionization of the Universe and to interpret spectra of observed galaxies at high redshift. To this end, we perform high-resolution, radiation- magneto-hydrodynamic simulations of GMCs with self-consistent star formation and stellar feedback. We find that a significant fraction (15–70%) of ionizing radiation escapes from the simulated GMCs 5 6 with different masses (10 and 10 M ), as the clouds are dispersed within about 2–5 Myr from the onset of star formation. The fraction of LyC and Lyα photons leaked is larger when the GMCs are less massive, metal-poor, more turbulent, and less dense. The most efficient leakage of LyC radiation occurs when the total star formation efficiency of a GMC is about 20%. The escape of Lyα shows a trend similar to that of LyC photons, except that the fraction of Lyα photons escaping from the 0.27 GMCs is larger (fLyα ≈ f900 ). The simulated GMCs show a characteristic velocity separation of ∆v ≈ 120 km s−1 in the time-averaged emergent Lyα spectra, suggesting that Lyα could be useful to infer the kinematics of the interstellar and circumgalactic medium. We show that Lyα luminosities are a useful indicator of the LyC escape, provided the number of LyC photons can be deduced through stellar population modeling. Finally, we find that the correlations between the escape fractions of Lyα, ultraviolet photons at 1500A,˚ and the Balmer α line are weak.

Keywords: Giant molecular clouds (653) — Photoionization (2060) – Reionization (1383) – Lyman- apha galaxies (978)

1. INTRODUCTION intergalactic space (e.g., Madau et al. 1999; Iliev et al. 2006; Mesinger & Furlanetto 2007; Gnedin 2014; Kimm

arXiv:2110.02975v1 [astro-ph.GA] 6 Oct 2021 Strong Lyman continuum (LyC) radiation produced by young massive stars is known to create low-density & Cen 2014; Norman et al. 2015; O’Shea et al. 2015; channels in giant molecular clouds (GMCs) (e.g., Dale Ocvirk et al. 2016; Pawlik et al. 2017; Rosdahl et al. et al. 2012; Geen et al. 2018; Kim et al. 2018; Grudić 2018; Finlator et al. 2018; Doussot et al. 2019; Hutter et al. 2021b). After escaping from their place of ori- et al. 2021). Therefore, understanding the propagation gin, the photons not only regulate the dynamics of the of LyC photons within GMCs is an important first step interstellar medium (ISM), but also drive the reioniza- to build a coherent picture of star formation, galaxy evo- tion of the Universe by ionizing neutral hydrogen in the lution, and reionization. Reionization theory requires a significant fraction (∼ 10%) of LyC photons to escape from their host dark [email protected] matter haloes in order to explain the end of reionization as well as the electron optical depth (Robertson et al. arXiv:2110.02973v1 [astro-ph.GA] 6 Oct 2021 1 Hz W ffadrmi uecn o ogritra (10 interval longer a for quiescent remain and off ai pcr sP,btisedo hi pia spectra, a optical definition their by have on galaxies instead (NLS1) but 1 Seyfert PS, their narrow-line on as based spectra not Classified radio jets. young kinematically have common jets. relativistic in the have of scenarios age kinematical these galaxies young of the radio Both classical 1995). devel- into al., still et grow (Fanti are eventually that will sources and young oping PS of which in population scenario, a age young represent the is option Another 2009). Stawarz, & Nikiel-Wroczyński, Janiuk, Siemiginowska, Czerny, nti ae nemtetatvt ftencesmyturn 10 may for nucleus it time. maintain the over and of same on activity the jet intermittent the remains case, tem- shape this a spectral In sources with the transient sources which be may Saikia, radio in & Others (O’Dea variable spectrum. bag inverted be mixed porarily may a Some is population 2021). PS The trum. ( mi:[email protected] Email: [email protected]; Email: Järvelä E. & Berton *M. Correspondence 4 3 2 1 ekdsucs(S r ls fpwru ( powerful of class a are (PS) sources Peaked Berton Marco story love a 1s: Seyfert narrow-line and sources Peaked TYPE ARTICLE xxx/xxxx DOI: 2016; April 26 Received < aCñd,Mdi,Spain Madrid, Cañada, de la Villanueva 28692 s/n, Castillo del Bajo Camino (ESAC), Centre Astronomy Space European (ESA), Agency Space European Finland Kylmälä, FI-02540 114, Metsähovintie Observatory, Radio Metsähovi University Aalto Finland Turku, of University FI-20014 5, Vesilinnantie Turku, of University (FINCA), ESO with Astronomy for Centre Finnish Chile 19001, Santiago 19, Casilla 3107, Córdova de Alonso (ESO), Observatory Southern European oee,P r o h nycaso ore hc may which sources of class only the not are PS However, 5kc fterrltvsi esadb ekdrdospec- radio peaked a by and jets relativistic their of kpc) 15 INTRODUCTION −1 ai ore hrceie yasallna size linear small a by characterized sources radio ) 1,2,3 eie ue2016; June 6 Revised mlaJärvelä Emilia | 3 -10 4 cetd6Jn 2016 June 6 Accepted r,adte unit turn then and yrs, cieglci uli ekdsucs arwln Seyfe narrow-line sources; peaked nuclei; galactic Active KEYWORDS: r obsolete. particular, in parameter and, radio-loudness PS of knowledge our invisi impact basically could GHz, NLS1s, 10 in jets radio absorbed of al discovery We developed. was unification review orientation-based brief this this why In PS. among hiding be could NLS1s emitting ore a hr ohcasfiain,adta ato h p the of part that and classifications, both i share Nowadays, can ago. sources years twenty already noticed were narro and galaxies (PS) sources peaked between similarities first The 4 logP 1 . GHz 4 4 -10 6 > yrs, 25 eioso G ihhge lc oemse (Mathur, masses pro- hole the be black may higher they with as one. AGN NLS1s, high-mass for of low-mass the case genitors the the accreting than sense, be younger time to could this considered less This reasonable In be spent is phase. can has it active one AGN mass an hole If in hole black cycle. efficiently black its life lower that same the a assume with the has redshift, them and same of structure the physical at located same AGN, two consider 2018; properties). Gallo, their on 2020; reviews for to Foschini, 2018, close Lister, see 2018; accreting Komossa, (e.g., be limit must they Eddington to AGN, AGN. the comparable massive other is more NLS1s of of of that that luminosity bolometric than the lower Since 1999) Wandel, & Peterson eso emte ie stpclyatiue olwrota- (10 low mass with to hole attributed black typically a around is velocity narrow- lines tional the permitted 2008), of Haardt, ness & hypotheses Fontana, alternative Dotti, central (Decarli, Although the 1 exist visible. and Type directly (BLR) 2 instead region are Type are broad-line engine as the they which obscured but in (AGN), not AGN, nuclei are galactic they active lines, permitted narrow s km 2000 than ulwdha afmxmmo h H the of maximum half at full-width hs rprisldt nitrsighptei.Ltus Let hypothesis. interesting an to led properties These −1 Otrrc og,18) ept such Despite 1985). Pogge, & (Osterbrock t1 nfidmodels unified 1; rt rn ouainof population arent l tfeunisbelow frequencies at ble -ieSyet1(NLS1) 1 Seyfert w-line oso o h recent the how show so skonta several that known is t edsrb o and how describe we , ne h ieyused widely the ender

msinln lower line emission 6 -10

8 -ray M ⊙ , Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX631

Why do M dwarfs have more transiting planets?

Gijs D. Mulders,1, 2, 3 Joanna Drążkowska,4 Nienke van der Marel,5 Fred J. Ciesla,6, 3 and Ilaria Pascucci7, 3

1Facultad de Ingeniería y Ciencias, Universidad Adolfo Ibáñez, Av. Diagonal las Torres 2640, Peñalolén, Santiago, Chile 7 @GijsMulders 2Millennium Institute for Astrophysics, Chile 3Alien Earths Team, NASA Nexus for Exoplanet System Science 4University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München, Scheinerstr. 1, 81679 Munich, Germany 5Physics & Astronomy Department, University of Victoria, 3800 Finnerty Road, Victoria, BC, V8P 5C2, Canada 6Department of the Geophysical Sciences, The University of Chicago, 5734 South Ellis Avenue, Chicago, IL 60637 7Lunar and Planetary Laboratory, The University of Arizona, Tucson, AZ 85721, USA

ABSTRACT We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in the outer disk to form first. When those giant planets reach pebble isolation mass they cut off the flow of pebbles to the inner disk and prevent the formation of close-in super-Earths. We apply this model to stars with masses between 0.1 and 2 M and for a range of initial disk masses. We find that the masses of hot super-Earths and of cold giant planets are anti-correlated. The fraction of our simulations that form hot super-Earths is higher around lower-mass stars and matches the exoplanet occurrence rates from Kepler. The fraction of simulations forming cold giant planets is consistent with the stellar mass dependence from radial velocity surveys. A key testable prediction of the pebble accretion hypothesis is that the occurrence rates of super-Earths should decrease again for M dwarfs near the sub-stellar boundary like Trappist-1.

Keywords: Exoplanets (498), Exoplanet formation (492), Planet formation (1241), Planetary system formation (1257), Protoplanetary disks (1300)

1. INTRODUCTION of exoplanet demographics becomes more clear, these Exoplanets provide important insights into the pro- planet formation hypotheses need to be evaluated against cesses that operate during planet formation. The dis- the known exoplanet populations. covery of super-Earths — planets with masses and radii One observed trend in exoplanets demographics that between that of Earth and Neptune and orbiting within has so far eluded a straightforward explanation is the 1 au of their host stars — have spawned several new elevated occurrence rate of super-Earths around low- mass stars (see, e.g. Mulders 2018 for a review). The

arXiv:2110.02971v1 [astro-ph.EP] 6 Oct 2021 hypotheses about how planets are assembled given that these planets are absent in our own solar system. Most occurrence rates of super-Earths observed with Kepler of these scenarios involve either the radial drift of pebble- increase from F stars to M dwarfs by a factor 3 (Howard sized solids or the migration of earth-sized proto-planets – et al. 2012; Mulders et al. 2015a), a result that has been in a contrast with formation models the solar system ter- recently conﬁrmed by radial velocity surveys (Sabotta restrial planets which are mostly “in situ”. As exoplanet et al. 2021). These higher occurrence rates have been surveys continue to discover planets and the overall view show to correspond to a higher fraction of stars with planetary systems (Yang et al. 2020; He et al. 2021). This trend deﬁes the scaling relations between host Corresponding author: Gijs D. Mulders star mass and giant planet occurrence (e.g. Johnson et al. [email protected] 2010; Ghezzi et al. 2018; Fulton et al. 2021), between host star mass and protoplanetary disk dust mass (e.g. Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX63

Internal r-process abundance spread of M15 and a single stellar population model

Yuta Tarumi,1 Naoki Yoshida ,1, 2, 3 and Shigeki Inoue4, 5, 6

1Department of Physics, School of Science, The University of Tokyo, Bunkyo, Tokyo 113-0033, Japan 2Kavli Institute for the Physics and Mathematics of the Universe (WPI), UTIAS, The University of Tokyo, Chiba, 277-8583, Japan 3Research Center for the Early Universe, School of Science, The University of Tokyo, Tokyo, 113-0033, Japan 4Center for Computational Sciences, University of Tsukuba, Ten-nodai, 1-1-1 Tsukuba, Ibaraki 305-8577, Japan 5Chile Observatory, National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 6Department of Cosmosciences, Graduates School of Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan

(Received; Revised; Accepted) Submitted to ApJ

ABSTRACT The member stars in globular cluster M15 show a substantial spread in the abundances of r-process elements. We argue that a rare and proliﬁc r-process event enriched the natal cloud of M15 in an inhomogeneous manner. To critically examine the possibility, we perform cosmological galaxy formation simulations and study the physical conditions for the inhomogeneous enrichment. We explore a large parameter space of the merger event time and the site. Our simulations reproduce the large r-process abundance spread if a neutron-star merger occurs at ∼ 100 pc away from the formation site of the cluster and in a limited time range of a few tens million years before the formation. Interestingly, a bimodal feature is found in the Eu abundance distribution in some cases, similarly to that inferred from recent observations. M15 member stars do not show clear correlation between the abundances of Eu and light elements such as Na that is expected in models with two stellar populations. We thus argue that a majority of the stars in M15 are formed in a single burst. Considering the abundances of the ﬁrst peak r-process elements such as Y and Zr, we conclude that the main r-process event should have a high lanthanide fraction of [Eu/Y] ∼ 1.0, consistent with the so-called r-II stars in the Milky Way and Reticulum II.

Keywords: Will be chosen from UAT thesaurus

1. INTRODUCTION appears that GCs host multiple populations commonly Globular clusters (GCs) are compact, nearly spherical (Bastian & Lardo 2018). clusters with old stellar populations. Typically, more Peculiar features are found in the chemical abundance than 105 stars gather within ∼ 10 pc, making the most patterns of GC stars. The abundances of light metals densely populated environments in the Universe. It was such as sodium (Na) and oxygen (O) show prominent thought that stars in a GC are of single population: all spread within each cluster. Moreover, the abundances arXiv:2110.02970v1 [astro-ph.GA] 6 Oct 2021 the stars are formed at once, and share the same stellar clearly (anti-)correlate with each other, suggesting that age and chemical abundance. Recent studies find sub- high-temperature hydrogen-burning have played a role. stantial scatter in the color-magnitude diagram (CMD) This is in stark contrast with the nearly uniform iron 1 (see Renzini 2008 and references therein), and the spread (Fe) abundance in GCs . is confirmed in almost all GCs older than 2 Gyr. It Elements heavier than Fe are also found in GCs and the abundances are measured. The heavy elements are mainly synthesized in two neutron-capture processes: Corresponding author: Yuta Tarumi [email protected] 1 An exception is ω centauri, which shows internal metallicity spread (e.g. Johnson & Pilachowski 2010). MNRAS 000,1–20 (2021) Preprint 8 October 2021 Compiled using MNRAS LATEX style file v3.0

The thesan project: Lyman-휶 emission and transmission during the Epoch of Reionization

A. Smith,1★ R. Kannan,2† E. Garaldi,3‡ M. Vogelsberger,1 R. Pakmor,3 V. Springel3 and L. Hernquist2 1Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 2Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 3Max-Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany

Accepted XXX. Received YYY; in original form ZZZ

ABSTRACT The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The thesan project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the arepo-rt radiation-magneto-hydrodynamics solver. In this initial study we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a (95.5 cMpc)3 region of the Universe. To avoid numerical artifacts we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp(−휏) behavior leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages. Key words: galaxies: high-redshift – cosmology: dark ages, reionization, first stars – radiative transfer – methods: numerical

1 INTRODUCTION (Planck Collaboration et al. 2020), the declining fraction of Lyman- alpha emitters (LAEs) among the galaxy population at 푧 & 6 (Stark The Epoch of Reionization (EoR) is the time period in the history of et al. 2011; Schenker et al. 2014; Mason et al. 2019), the imprint of the Universe when the radiation from the first stars and galaxies ini- neutral hydrogen in the IGM as a damping wing absorption feature tiated a cosmic phase transition throughout the intergalactic medium on the spectrum of high-redshift quasars (Simcoe et al. 2012; Davies (IGM), which went from being cold and neutral to warm and ionized. et al. 2018), and the spatial fluctuations of the Ly훼 forest transmis- In recent years we have witnessed significant progress in understand- sion (Oñorbe et al. 2017; Kulkarni et al. 2019). We anticipate a more ing the reionization process and advancing the current observational complete understanding of the EoR from upcoming facilities, includ- and computational frontiers. However, many of the central questions ing the James Webb Space Telescope (JWST) for the characterization arXiv:2110.02966v1 [astro-ph.CO] 6 Oct 2021 only have tentative answers that may be revised as the constraints of high-푧 galaxies and the Low-Frequency Array (LOFAR), Hydro- from available data improve. For example, there is significant uncer- gen Epoch of Reionization Array (HERA), and Square Kilometer tainty about the role and contribution of low- and high-mass galaxies, Array (SKA) for 21 cm cosmology measurements to map out the the sources and escape of ionizing photons, and even the timing, du- distribution of neutral hydrogen in the Universe. On the theory side, ration, and morphology of reionization itself (Barkana & Loeb 2001; radiation hydrodynamic (RHD) simulations have played an increas- Loeb & Furlanetto 2013; Dayal & Ferrara 2018; Wise 2019). ingly important role in capturing the nonlinear, multiscale physics There is mounting evidence for so-called late reionization in which of reionization although there are still important computational chal- most of the IGM is rapidly ionized around redshift 푧 ∼ 7–8 (Naidu lenges to overcome in the coming decades as well (Ciardi et al. 2000; et al. 2020). This understanding comes from the Planck measure- Iliev et al. 2006; Gnedin 2014; Ocvirk et al. 2016; Pawlik et al. 2017; ment of a low optical depth for electron scattering of CMB photons Rosdahl et al. 2018). Ultimately, the joint analysis of observational probes sensitive to ★ E-mail: [email protected]; NHFP Einstein Fellow. unique aspects of galaxy formation and IGM properties will provide † E-mail: [email protected] definitive answers to the main questions about the EoR. It is in this ‡ E-mail: [email protected] spirit that we pursue an initial study of Lyman-alpha (Ly훼) emission

Reading the CARDs: the Imprint of Accretion History in the Chemical Abundances of the Milky Way’s Stellar Halo

EMILY C.CUNNINGHAM,1 ROBYN E.SANDERSON,2, 1 KATHRYN V. JOHNSTON,3, 1 NONDH PANITHANPAISAL,2 MELISSA K.NESS,3, 1 ∗ ANDREW WETZEL,4 SARAH R.LOEBMAN,5 IVANNA ESCALA,6, 7, 8 , DANNY HORTA,9 AND CLAUDE-ANDRE´ FAUCHER-GIGUERE` 10

1Center for Computational Astrophysics, Flatiron Institute, 162 5th Ave., New York, NY 10010, USA 2Department of Physics and Astronomy, University of Pennsylvania, 209 S 33rd St., Philadelphia, PA 19104, USA 3Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY, 10027, USA 4Department of Physics & Astronomy, University of California, Davis, CA 95616, USA 5Department of Physics, University of California, Merced, 5200 Lake Road, Merced, CA 95343, USA 6Department of Astronomy, California Institute of Technology, 1200 E California Blvd, Pasadena, CA 91125, USA 7The Observatories of the Carnegie Institution for Science, 813 Santa Barbara St, Pasadena, CA 91101, USA 8Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA 9Astrophysics Research Institute, Brownlow Hill, Liverpool L3 5RF 10Department of Physics and Astronomy and CIERA, Northwestern University, 1800 Sherman Ave, Evanston, IL 60201, USA

ABSTRACT In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies’ assembly histories by modelling the CARDs of the stellar halos of the Latte galaxies as a linear combination of “template” CARDs from disrupted dwarfs, with different stellar masses M? and quenching times t100. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of < 10%. For the fraction of mass accreted as a function of t100, we find residuals of 20–30% for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Though the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos.

Keywords: Milky Way stellar halo, chemical abundances, galaxy chemical evolution

1. INTRODUCTION properties of MW halo stars, we have the unique opportunity In the hierarchical paradigm of galaxy evolution, the Milky study galaxy formation and evolution in the high redshift uni- Way (MW) is more than just one galaxy. As many as thou- verse: we can study the remains of galaxies that did not survive to present day, including populations that are too faint to arXiv:2110.02957v1 [astro-ph.GA] 6 Oct 2021 sands of dwarf galaxies are predicted to have contributed to the assembly of the MW dark matter halo over cosmic time; observe directly at high redshift, on a star by star basis. the remains of these accreted systems make up the compo- In the current era of MW studies, we have unprecedented nent of our Galaxy known as the stellar halo. The forma- knowledge of the motions and chemical properties of the tion conditions of halo stars are imprinted in their chemical stars in our Galaxy. This is in large part due to the data abundances. Therefore, through observations of the chemical releases from the Gaia mission (Gaia Collaboration et al. 2016, 2018, 2021) and as well as complementary wide-ﬁeld spectroscopic programs, such as the GALAH survey (De Corresponding author: Emily C. Cunningham Silva et al. 2015), the H3 survey (Conroy et al. 2019b), ecunningham@ﬂatironinstitute.org APOGEE (Majewski et al. 2017), RAVE (Steinmetz et al. ∗ Carnegie-Princeton Fellow 2006), SEGUE (Yanny et al. 2009), and LAMOST (Cui et al. 2012). Our knowledge of the chemodynamical structure of Draft version October 8, 2021 Typeset using LATEX twocolumn style in AASTeX63

The Criterion for Chaos in Three-Planet Systems

Jeremy Rath,1, 2 Sam Hadden,3, 4 and Yoram Lithwick1, 2

1Dept. of Physics and Astronomy, Northwestern University, 2145 Sheridan Rd., Evanston, IL 60208, USA 2Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA), 1800 Sherman, Evanston, IL 60201, USA 3Canadian Institute for Theoretical Astrophysics, 60 St George St Toronto, ON M5S 3H8, Canada 4Center for Astrophysics | Harvard & Smithsonian, Cambridge, MA 02138, USA

ABSTRACT We establish the criterion for chaos in three-planet systems, for systems similar to those discovered by the Kepler spacecraft. Our main results are as follows: (i) The simplest criterion, which is based on overlapping mean motion resonances (MMR’s), only agrees with numerical simulations at a very crude level. (ii) Much greater accuracy is attained by considering neighboring MMR’s that do not overlap. We work out the width of the chaotic zones around each of the neighbors, and also provide simple approximate expressions for the widths. (iii) Even greater accuracy is provided by the overlap of three-body resonances (3BR’s), which accounts for fine-grained structure seen in maps from N-body simulations, and also predicts the Lyapunov times. Previous studies conflict on whether overlap of MMR’s or of 3BR’s drive interplanetary chaos. We show that both do, and in fact they are merely different ways of looking at the same effect. (iv) We compare both criteria with high-resolution maps of chaos from N-body simulations, and show that they agree at a high level of detail.

Keywords: Celestial mechanics - Orbits - Orbital resonances, Exoplanet dynamics

1. INTRODUCTION the chaos in planetary systems. But which resonances In the Solar System, the planets’ orbits are chaotic overlap depends on the system, as planetary systems (Laskar 1996; Laskar & Gastineau 2009; Holman & Mur- harbor a variety of different kinds of resonances. In ray 1996). Many exoplanetary systems are also chaotic, nearly circular two-planet systems, decreasing the plan- or at least close to the threshold for chaos, e.g., Kepler- ets’ separation drives them toward chaos. There is a 11 (Mahajan & Wu 2014; Yee et al. 2021), Kepler-36 critical separation inside of which the orbits are always (Deck et al. 2012), Kepler-102 (Volk & Malhotra 2020), chaotic, and its value is determined by where neighbor- GJ876 (Go´zdziewskiet al. 2002; Batygin et al. 2015a), ing first-order mean motion resonances (MMR’s) over- and Kepler-24, -85, and -444 (Yee et al. 2021). The lap (Wisdom 1980; Deck et al. 2013). Two-planet sys- chaos cannot be too violent, because otherwise the systems can also be driven towards chaos by increasing tems would not survive for billions of years. That raises their eccentricity, because the widths of MMR’s increase the question of how planets were emplaced into such with increasing eccentricity. The critical eccentricity for delicate configurations. It also opens the door to using chaos has been determined by, e.g., Mardling(2008); observed planetary configurations to learn about their Hadden & Lithwick(2018). It depends not only on first- arXiv:2110.02956v1 [astro-ph.EP] 6 Oct 2021 early history. But our understanding of how chaos works order MMR’s, but on higher-order MMR’s as well. In in planetary systems is not yet sufficient to answer such systems of more than two planets, Tamayo et al.(2021) questions. suggest that chaos is generally attributable to the over- A simple and powerful criterion for the presence lap of two-body MMR’s—similar to the two-planet case. of chaos is resonance overlap (Walker & Ford 1969; Overlap of other kinds of resonances also produces Chirikov 1979). Resonance overlap accounts for much of chaos. In the inner Solar System, chaos is caused by the overlap of secular resonances (Laskar 1996; Lithwick & Wu 2011; Batygin et al. 2015b). And in the outer So- Corresponding author: Jeremy Rath lar System, chaos has been attributed to the overlap of [email protected] three-body resonances (3BR’s; Murray & Holman 1999). Quillen(2011) and Petit et al.(2020) suggest that chaos