Abstract visibility
Abstract text size

Papers for Thursday, Oct 12 2023

Papers with local authors

Changseok Kim, Jong-Hak Woo, Rongxin Luo, Aeree Chung, Junhyun Baek, Huynh Anh N. Le, Donghoon Son

27 pages, 18 figures, Accepted for publication in ApJ

0 votes
Paper 11 — arXiv:2310.06928
0 votes
Paper 11 — arXiv:2310.06928

We present spatially resolved gas kinematics, ionization, and energetics of 11 type 1 and 5 type 2 active galactic nuclei (AGNs) with strong ionized gas outflows at z $<0.3$ using Gemini Multi-Object Spectrograph Integral Field Unit (GMOS-IFU) data. We find a strongly blueshifted region in [OIII] velocity maps, representing an approaching cone in biconical outflows, and blueshifted and redshifted regions in H$\alpha$ velocity maps, which show gravitationally rotating kinematics. AGN photoionization is dominant in the central region of most targets, and some of them also show ring-like structures of LINER or composite that surround the AGN-dominated center. Following our previous studies, we kinematically determine outflow sizes by the ratio between [OIII] and stellar velocity dispersion. Outflow sizes of type 1 AGNs follow the same kinematic outflow size-[OIII] luminosity relation obtained from the type 2 IFU sample in Kang & Woo and Luo (updated slope $0.29\pm0.04$), while they are limited to the central kpc scales, indicating the lack of global impact of outflows on the interstellar medium. Small mass outflow rates and large star formation rates of the combined sample support that there is no evidence of rapid star formation quenching by outflows, which is consistent with the delayed AGN feedback.

D. Caprioli, G. Zacharegkas, C. Haggerty, S. Gupta, B. Schroer

8 pages, 3 figures, Proceedings of the 38th International Cosmic Ray Conference (ICRC2023), 26 July - 3 August, 2023, Nagoya, Japan

0 votes
Paper 21 — arXiv:2310.07038
0 votes
Paper 21 — arXiv:2310.07038

The non-resonant (Bell) streaming instability driven by energetic particles is crucial for producing amplified magnetic fields that are key to the acceleration of cosmic rays (CRs) in supernova remnants, around Galactic and extra-galactic CR sources, and for the CR transport. We present a covariant theory for the saturation of the Bell instability, substantiated by self-consistent kinetic simulations, that can be applied to arbitrary CR distributions and discuss its implications in several heliospheric and astrophysical contexts.

Na Wang, Qian Xu, Jun Ma, Zhiyong Liu, Qi Liu, Hailong Zhang, Xin Pei, Maozheng Chen, Richard N. Manchester, Kejia Lee, Xingwu Zheng, Hans J. Kärcher, Wulin Zhao, Hongwei Li, Dongwei Li, Martin Süss, Matthias Reichert, Zhongyi Zhu, Congsi Wang, Mingshuai Li, Rui Li, Ning Li, Guljaina Kazezkhan, Wenming Yan, Gang Wu, Lang Cui, Ming Zhang, Haitao Li

12 pages, 11 figures, accepted for publication in Science China Physics, Mechanics & Astronomy

0 votes
Paper 35 — arXiv:2310.07163
0 votes
Paper 35 — arXiv:2310.07163

This study presents a general outline of the Qitai radio telescope (QTT) project. Qitai, the site of the telescope, is a county of Xinjiang Uygur Autonomous Region of China, located in the east Tianshan Mountains at an elevation of about 1800 m. The QTT is a fully steerable, Gregorian type telescope with a standard parabolic main reflector of 110 m diameter. The QTT has adopted an um-brella support, homology-symmetric lightweight design. The main reflector is active so that the deformation caused by gravity can be corrected. The structural design aims to ultimately allow high-sensitivity observations from 150 MHz up to 115 GHz. To satisfy the requirements for early scientific goals, the QTT will be equipped with ultra-wideband receivers and large field-of-view mul-ti-beam receivers. A multi-function signal-processing system based on RFSoC and GPU processor chips will be developed. These will enable the QTT to operate in pulsar, spectral line, continuum and Very Long Baseline Interferometer (VLBI) observing modes. Electromagnetic compatibility (EMC) and radio frequency interference (RFI) control techniques are adopted throughout the system design. The QTT will form a world-class observational platform for the detection of low-frequency (nanoHertz) gravitational waves through pulsar timing array (PTA) techniques, pulsar surveys, the discovery of binary black-hole systems, and exploring dark matter and the origin of life in the universe.

S.-X. Yi, C.-W. Wang, B. Zhang, S.-L. Xiong, S.-N. Zhang, W.-J. Tan, J.-C. Liu, W.-C. Xue, Y.-Q. Zhang, C. Zheng, R. Moradi, Y. Wang, P. Zhang, Z.-H. An, C. Cai, P.-Y. Feng, K. Gong, D.-Y. Guo, Y. Huang, B. Li, X.-B. Li, X.-Q. Li, X.-J. Liu, Y.-Q. Liu, X. Ma, W.-X. Peng, R. Qiao, L.-M. Song, J. Wang, P. Wang, Y. Wang, X.-Y. Wen, S. Xiao, Y.-B. Xu, S. Yang, Q.-B. Yi, D.-L. Zhang, F. Zhang, H.-M. Zhang, J.-P. Zhang, Z. Zhang, X.-Y. Zhao, Y. Zhao, S.-J. Zheng

6 pages and 2 figures in the main text. 19 pages and 5 figures in total

0 votes
Paper 38 — arXiv:2310.07205
0 votes
Paper 38 — arXiv:2310.07205

The origin of the prompt emission of gamma-ray bursts (GRBs) is still subject to debate because of the not-well-constrained jet composition, location of the emission region, and mechanism with which $\gamma$-rays are produced. For the bursts whose emission is dominated by non-thermal radiation, two leading paradigms are internal shock model invoking collisions of matter-dominated shells and internal-collision-induced magnetic reconnection and turbulence (ICMART) model invoking collisions of magnetically-dominated shells. These two models invoke different emission regions and have distinct predictions on the origin of light curve variability and spectral evolution. The second brightest GRB in history, GRB230307A, provides an ideal laboratory to study the details of GRB prompt emission thanks to its extraordinarily high photon statistics and the single broad pulse shape characterized by an energy-dependent fast-rise-exponential-decay (FRED) profile. Here we demonstrate that its broad pulse is composed of many rapidly variable short pulses, rather than being the superposition of many short pulses on top of a slow component. Such a feature is consistent with the ICMART picture, which envisages many mini-jets due to local magnetic reconnection events in a large emission zone far from the GRB central engine, but raises a great challenge to the internal shock model that attributes fast and slow variability components to shocks at different radii with the emission being the superposition of various components. The results provide strong evidence for a Poynting-flux-dominated jet composition of this bright GRB.

Aviad Levis, Andrew A. Chael, Katherine L. Bouman, Maciek Wielgus, Pratul P. Srinivasan
0 votes
Paper 75 — arXiv:2310.07687
0 votes
Paper 75 — arXiv:2310.07687

The interaction between the supermassive black hole at the center of the Milky Way, Sagittarius A$^*$, and its accretion disk, occasionally produces high energy flares seen in X-ray, infrared and radio. One mechanism for observed flares is the formation of compact bright regions that appear within the accretion disk and close to the event horizon. Understanding these flares can provide a window into black hole accretion processes. Although sophisticated simulations predict the formation of these flares, their structure has yet to be recovered by observations. Here we show the first three-dimensional (3D) reconstruction of an emission flare in orbit recovered from ALMA light curves observed on April 11, 2017. Our recovery results show compact bright regions at a distance of roughly 6 times the event horizon. Moreover, our recovery suggests a clockwise rotation in a low-inclination orbital plane, a result consistent with prior studies by EHT and GRAVITY collaborations. To recover this emission structure we solve a highly ill-posed tomography problem by integrating a neural 3D representation (an emergent artificial intelligence approach for 3D reconstruction) with a gravitational model for black holes. Although the recovered 3D structure is subject, and sometimes sensitive, to the model assumptions, under physically motivated choices we find that our results are stable and our approach is successful on simulated data. We anticipate that in the future, this approach could be used to analyze a richer collection of time-series data that could shed light on the mechanisms governing black hole and plasma dynamics.

All other papers

Ramin Barak, Konstantin Belotsky, Ekaterina Shlepkina

12 pages, 3 figures; contents of this paper correspond to the one that has been published in the journal Universe

The origin of an anomalous excess of high-energy (about 100 GeV and higher) positrons in cosmic rays is one of the rare problems in this field, which is proposed to be solved with dark matter (DM). Attempts to solve this problem are faced with the issue of having to satisfy the data on cosmic positrons and cosmic gamma radiation, which inevitably accompanies positron production, such as FSR (final state radiation), simultaneously. We have been trying to come up with a solution by means of two approaches: making assumptions (*) about the spatial distribution of the dark matter and (**) about the physics of its interactions. This work is some small final step of a big investigation regarding the search for gamma suppression by employing the second approach, and a model with a doubly charged particle decaying into two positrons (X++ $\rightarrow$ e+ e+) is suggested as the most prospective one from those considered before.

The formation of galaxies with warm dark matter is approximately adiabatic. The cold dark matter limit is singular and requires relaxation. In these lecture notes we develop, step-by-step, the physics of galaxies with warm dark matter, and their formation. The theory is validated with observed spiral galaxy rotation curves. These observations constrain the properties of the dark matter particles.

Giulia Santucci, Claudia Del P. Lagos, Katherine E. Harborne, Aaron Ludlow, Caro Foster, Richard McDermid, Adriano Poci, Katy L. Proctor, Sabine Thater, Glenn van de Ven, Ling Zhu, Daniel Walo Martin

Submitted to MNRAS, 21 pages

The merger history of a galaxy is thought to be one of the major factors determining its internal dynamics, with galaxies having undergone different types or mergers (e.g. dry, minor or major mergers) predicted to show different dynamical properties. We study the instantaneous orbital distribution of galaxies in the Eagle simulation, colouring the orbits of the stellar particles by their stellar age, in order to understand whether stars form in particular orbits (e.g. in a thin or thick disc). We first show that Eagle reproduces well the observed stellar mass fractions in different stellar orbital families as a function of stellar mass and spin parameter at z = 0. We find that the youngest stars reside in a thin disc component that can extend to the very inner regions of galaxies, and that older stars have warmer orbits, with the oldest ones showing orbits consistent with both hot and counter-rotating classifications, which is consistent with the trend found in the Milky-Way and other disc galaxies. We also show that counter-rotating orbits trace galaxy mergers - in particular dry mergers, and that in the absence of mergers, counter-rotating orbits can also be born from highly misaligned gas accretion that leads to star formation.

Erica J. Nelson, Gabriel Brammer, Clara Gimenez-Arteaga, Pascal A. Oesch, Hannah Ubler, Anna de Graaff, Jasleen Matharu, Rohan P. Naidu, Alice E. Shapley, Katherine E. Whitaker, Emily Wisnioski, Natascha M. Forster Schreiber, Renske Smit, Pieter van Dokkum, John Chisholm, Ryan Endsley, Abigail I. Hartley, Justus Gibson, Emma Giovinazzo, Garth Illingworth, Ivo Labbe, Michael V. Maseda, Jorryt Matthee, Alba Covelo Paz, Sedona H. Price, Naveen A. Reddy, Irene Shivaei, Andrea Weibel, Stijn Wuyts, Mengyuan Xiao, Stacey Alberts, William M. Baker, Andrew J. Bunker, Alex J. Cameron, Stephane Charlot, Daniel J. Eisenstein, Zhiyuan Ji, Benjamin D. Johnson, Gareth C. Jones, Roberto Maiolino, Brant Robertson, Lester Sandles, Katherine A. Suess, Sandro Tacchella, Christina C. Williams, Joris Witstok

Fig. 3 shows the main result

With the remarkable sensitivity and resolution of JWST in the infrared, measuring rest-optical kinematics of galaxies at $z>5$ has become possible for the first time. This study pilots a new method for measuring galaxy dynamics for highly multiplexed, unbiased samples by combining FRESCO NIRCam grism spectroscopy and JADES medium-band imaging. Here we present one of the first JWST kinematic measurements for a galaxy at $z>5$. We find a significant velocity gradient, which, if interpreted as rotation yields $V_{rot} = 240\pm50$km/s and we hence refer to this galaxy as Twister-z5. With a rest-frame optical effective radius of $r_e=2.25$kpc, the high rotation velocity in this galaxy is not due to a compact size as may be expected in the early universe but rather a high total mass, ${\rm log(M}_{dyn}/{\rm M}_\odot)=11.0\pm0.2$. This is a factor of roughly 4x higher than the stellar mass within the effective radius. We also observe that the radial H$\alpha$ equivalent width profile and the specific star formation rate map from resolved stellar population modeling is centrally depressed by a factor of $\sim1.5$ from the center to $r_e$. Combined with the morphology of the line-emitting gas in comparison to the continuum, this centrally suppressed star formation is consistent with a star-forming disk surrounding a bulge growing inside-out. While large, rapidly rotating disks are common to z~2, the existence of one after only 1Gyr of cosmic time, shown for the first time in ionized gas, adds to the growing evidence that some galaxies matured earlier than expected in the history of the universe.

Jindra Gensior, Timothy A. Davis, Martin Bureau, J. M. Diederik Kruijssen, Michele Cappellari, Ilaria Ruffa, Thomas G. Williams

14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of structure within a gas reservoir. Specifically here, we focus on the power spectrum slope and aim to constrain whether the shear induced by a dominant spheroidal potential can induce sufficient turbulence to suppress fragmentation, resulting in the smooth central gas discs observed. We compute surface density power spectra for the nuclear gas reservoirs of fourteen simulated isolated galaxies and twelve galaxies observed as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project. Both simulated and observed galaxies range from disc-dominated galaxies to spheroids, with central stellar mass surface densities, a measure of bulge dominance, varying by more than an order of magnitude. For the simulations, the power spectra steepen with increasing central stellar mass surface density, thereby clearly linking the suppression of fragmentation to the shear-driven turbulence induced by the spheroid. The WISDOM observations show a different (but potentially consistent) picture: while there is no correlation between the power spectrum slopes and the central stellar mass surface densities, the slopes scatter around a value of 2.6. This is similar to the behaviour of the slopes of the simulated galaxies with high central stellar mass surface densities, and could indicate that high shear eventually drives incompressible turbulence.

Massimo Dotti, Fabio Rigamonti, Stefano Rinaldi, Walter Del Pozzo, Roberto Decarli, Riccardo Buscicchio

9 pages, 6 figures, accepted for publication on Astronomy & Astrophysics

We present a new observational test to identify massive black hole binaries in large multi-epoch spectroscopical catalogues and to probe the real nature of already proposed binary candidates. The test is tailored for binaries with separations large enough to allow each black hole to retain its own broad line region. In this limit the fast AGN variability typically observed over months cannot be associated to the much longer binary period, and it is assumed (as for the case of single black holes) to be the consequence of the evolution of the innermost regions of the two accretion discs. A simple analysis of the cross-correlation between different parts of individual broad emission lines can therefore identify the presence of two massive black holes whose continua vary independently of each other. Our analysis indicates that to be less affected by the noise in the spectra the broad lines should be divided in two close-to-equal-flux parts. This ensures that in the single massive black hole scenario the cross-correlation will always be high. With monitoring campaigns similar to those performed for reverberation mapping studies, on the other way, a binary can show any value of the cross-correlation and can therefore be distinguished from a standard AGN. The new test can be performed over timescales orders of magnitude shorter than the alternative tests already discussed in literature, and can be a powerful complement to the massive black hole binary search strategies already in place.

Ritik Sharma, Mahavir Sharma (IIT Bhilai)

9 pages, 6 figures, submitted to MNRAS, comments welcome

We investigate the steady spherically symmetric accretion in the combined potential of a central black hole and a dark matter halo. For the halo, we consider a Hernquist and an NFW potential and calculate the critical points of the flow. We find that the trans-sonic solution to the centre is not possible without a black hole, whereas two types of trans-sonic solutions are possible in its presence. We also derive the mass accretion rate for a black hole at the centre of a dark matter halo. Our results indicate two phases of accretion. The first is an initial phase with a low accretion rate that depends on the black hole mass, followed by a second phase with a high accretion rate that depends on the halo mass. In the second phase, the black hole mass increases rapidly to supermassive scales, which explains the existence of quasars at redshift $z\ge 6$ and also the recently detected supermassive black hole (SMBH) by the James Webb Space Telescope (JWST). Further, we calculate the evolution of the Eddington ratio for growing black holes. The accretion is mostly sub-Eddington except for a short super-Eddington episode when the mass accretion rate transitions from low to high. However, during that episode, the black hole mass is likely inadequate to hinder accretion through radiative feedback.

Jason Young, Alexandra Pope, Anna Sajina, Lin Yan, Thiago S Goncalves, Miriam Eleazer, Stacey Alberts, Lee Armus, Matteo Bonato, Daniel A. Dale, Duncan Farrah, Carl Ferkinhoff, Christopher C. Hayward, Jed McKinney, Eric J. Murphy, Nicole Nesvadba, Patrick Ogle, Leonid Sajkov, Sylvain Veilleux

10 pages, 5 figures, 1 table

We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydrocarbon (PAH) features at 3.3 micron, 6.2 micron, and 11.3 micron, and the warm molecular gas indicators H2S(5) and H2S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the Northern side of the galaxy when compared to the well-studied star-formation tracer [Ne II]. The H2 rotational lines, which are shock indicators, are strongest and most extended on the Southern side of the galaxy. [Ar II] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe which can be detected with JWST out to z=3. Velocity maps of [Ar II] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and ISM conditions in a galaxy halfway back to the peak epoch of galaxy evolution.

H. Rueda, F. Bradascio, J.A. Barrio, J. Biteau, F. Brun, C. Champion, J-F. Glicenstein, D. Hoffmann, P. Jean, J.P. Lenain, F. Louis, A. Pérez, M. Punch, P. Sizun, K-H. Sulanke, L.A. Tejedor, B. Vallage (for the CTA NectarCAM Project)

Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 ( arXiv:2309.08219 )

NectarCAM is a Cherenkov camera that will be installed on Medium-Sized Telescopes of the northern array of the Cherenkov Telescope Array Observatory (CTAO). It is composed of 265 modules, each of which includes 7 photo-multiplier tubes, a Front-End Board and a camera trigger system for data collection. The first NectarCAM unit is currently being integrated at CEA Paris-Saclay in France. Once installed at the CTAO's northern site, the NectarCAM's timing abilities will be crucial for reducing noise in images, improving image cleaning, and distinguishing between gamma-ray photons and cosmic-ray background. Additionally, it will enable coincidence identification with neighboring telescopes for stereoscopic observations. The timing system of NectarCAM has been tested in a dark room with various light sources. The results of the tests, including timing precision and accuracy of the trigger arrival relative to a laser source, and the timing of individual and multiple pixel signals, will be presented.

Georg Feulner (1), Mona Bukenberger (1, 2), Stefan Petri (1) ((1) Potsdam Institute for Climate Impact Research, (2) ETH Zürich)

18 pages, 8 figures. Published in Earth Syst. Dynam. (open access); accepted version before copy-editing

The instability with respect to global glaciation is a fundamental property of the climate system caused by the positive ice-albedo feedback. The atmospheric carbon dioxide concentration at which this Snowball bifurcation occurs changes through Earth's history because of the slowly increasing solar luminosity. Quantifying this critical CO$_2$ level is not only interesting from a climate dynamics perspective, but also a prerequisite for understanding past Snowball Earth events as well as the conditions for habitability on Earth and other planets. Earlier studies are limited to investigations with simple climate models for Earth's entire history, or studies of individual time slices carried out with a variety of more complex models and for different boundary conditions, making comparisons and the identification of secular changes difficult. Here we use a coupled climate model of intermediate complexity to trace the Snowball bifurcation of an aquaplanet through Earth's history in one consistent model framework. We find that the critical CO$_2$ concentration decreases more or less logarithmically with increasing solar luminosity until about 1 billion years ago, but drops faster in more recent times. Furthermore, there is a fundamental shift in the dynamics of critical states about 1.2 billion years ago, driven by the interplay of wind-driven sea-ice dynamics and the surface energy balance: For critical states at low solar luminosities, the ice line lies in the Ferrel cell, stabilised by the poleward winds despite moderate meridional temperature gradients under strong greenhouse warming. For critical states at high solar luminosities on the other hand, the ice line rests at the Hadley-cell boundary, stabilised against the equatorward winds by steep meridional temperature gradients resulting from the increased solar energy input at lower latitudes and stronger Ekman transport in the ocean.

Andreas Schanz, Florian List, Oliver Hahn

9 pages, 8 figures, to be submitted to OJA, comments welcome

In recent years, deep learning models have been successfully employed for augmenting low-resolution cosmological simulations with small-scale information, a task known as "super-resolution". So far, these cosmological super-resolution models have relied on generative adversarial networks (GANs), which can achieve highly realistic results, but suffer from various shortcomings (e.g. low sample diversity). We introduce denoising diffusion models as a powerful generative model for super-resolving cosmic large-scale structure predictions (as a first proof-of-concept in two dimensions). To obtain accurate results down to small scales, we develop a new "filter-boosted" training approach that redistributes the importance of different scales in the pixel-wise training objective. We demonstrate that our model not only produces convincing super-resolution images and power spectra consistent at the percent level, but is also able to reproduce the diversity of small-scale features consistent with a given low-resolution simulation. This enables uncertainty quantification for the generated small-scale features, which is critical for the usefulness of such super-resolution models as a viable surrogate model for cosmic structure formation.

A. Alvarado, T. Capistrán, I. Torres, J. R. SacahuÍ, R. Alfaro (for the HAWC collaboration)

Presented at the 38th International Cosmic Ray Conference (ICRC2023), Nagoya, Japan

HAWC is a ground-based observatory consisting of 300 water Cherenkov detectors, which observes the extensive air showers induced by cosmic rays from some TeV to a few PeV and, in particular, gamma rays from 300 GeV to more than 100 TeV. One of the crucial features required for a detector of extensive air showers is the estimation of the primary energy of the events to study the spectra of cosmic and gamma rays. For HAWC there are currently two gamma-ray energy estimators: one relies on a ground density parameter, while the other utilizes an artificial neural network. For the cosmic ray energy estimation, there is only one estimator based on maximum likelihood procedures and measurements of the lateral charge distribution of the events. It is worthwhile to update the cosmic-ray energy estimator due to recent improvements of the extensive air shower offline-reconstruction techniques in HAWC. Therefore, we implemented an artificial neural network to reconstruct the primary energy of hadronic events trained with several observables that characterize the air showers. We trained several models and evaluated their performance against the existing cosmic ray energy estimator. In this work, we present the features and performance of these models.

We present three new brown dwarf spectral binary candidates: CWISE J072708.09$-$360729.2, CWISE J103604.84$-$514424.4, and CWISE J134446.62$-$732053.9, discovered by citizen scientists through the Backyard Worlds: Planet 9 project. Follow-up near-infrared spectroscopy shows that each of these objects is poorly fit by a single near-infrared standard. We constructed binary templates and found significantly better fits, with component types of L7+T4 for CWISE J072708.09$-$360729.2, L7+T4 for CWISE J103604.84$-$514424.4, and L7+T7 for CWISE J134446.62$-$732053.9. However, further investigation of available spectroscopic indices for evidence of binarity and large amplitude variability suggests that CWISE J072708.09$-$360729.2 may instead be a strong variability candidate. Our analysis offers tentative evidence and characterization of these peculiar brown dwarf sources, emphasizing their value as promising targets for future high-resolution imaging or photometric variability studies.

N. Jannsen, J. De Ridder, D. Seynaeve, S. Regibo, R. Huygen, P. Royer, C. Paproth, D. Grießbach, R. Samadi, D. R. Reese, M. Pertenais, E. Grolleau, R. Heller, S. M. Niemi, J. Cabrera, A. Börner, S. Aigrain, J. McCormac, P. Verhoeve, P. Astier, N. Kutrowski, B. Vandenbussche, A. Tkachenko, C. Aerts

30 pages, 22 figures

PLAnetary Transits and Oscillations of stars (PLATO) is the ESA M3 space mission dedicated to detect and characterise transiting exoplanets including information from the asteroseismic properties of their stellar hosts. The uninterrupted and high-precision photometry provided by space-borne instruments such as PLATO require long preparatory phases. An exhaustive list of tests are paramount to design a mission that meets the performance requirements, and as such, simulations are an indispensable tool in the mission preparation. To accommodate PLATO's need of versatile simulations prior to mission launch - that at the same time describe accurately the innovative but complex multi-telescope design - we here present the end-to-end PLATO simulator specifically developed for the purpose, namely \texttt{PlatoSim}. We show step-by-step the algorithms embedded into the software architecture of \texttt{PlatoSim} that allow the user to simulate photometric time series of CCD images and light curves in accordance to the expected observations of PLATO. In the context of the PLATO payload, a general formalism of modelling, end-to-end, incoming photons from the sky to the final measurement in digital units is discussed. We show the strong predictive power of \texttt{PlatoSim} through its diverse applicability and contribution to numerous working groups within the PLATO Mission Consortium. This involves the on-going mechanical integration and alignment, performance studies of the payload, the pipeline development and assessments of the scientific goals. \texttt{PlatoSim} is a state-of-the-art simulator that is able to produce the expected photometric observations of PLATO to a high level of accuracy. We demonstrate that \texttt{PlatoSim} is a key software tool for the PLATO mission in the preparatory phases until mission launch and prospectively beyond.

G.-Q. Zhang, D. Bégué, A. Pe'er, B.-B. Zhang

29 pages, submitted to ApJ

The study of Magnetically Arrested Disks (MAD) attract strong interest in recent years, as these disk configurations were found to generate strong jets as observed in many accreting systems. Here, we present the results of 14 general relativistic magnetohydrodynamic(GRMHD) simulations of advection dominated accretion flow in the MAD state across black hole spins, carried with cuHARM. Our main findings are as follows. (i) The jets transport a significant amount of angular momentum to infinity in the form of Maxwell stresses. For positive, high spin, the rate of angular momentum transport is about 5 times larger than for negative spin. This contribution is nearly absent for a non-rotating black hole. (ii) The mass accretion rate and the MAD parameter, both calculated at the horizon, are not correlated. However, their time derivatives are anti-correlated for every spin. (iii) For zero spin, the contribution of the toroidal component of the magnetic field to the magnetic pressure is negligible, while for fast spinning black hole, it is in the same order as the contribution of the radial magnetic component. For high positive spin, the toroidal component even dominates. (iv) For negative spins, the jets are narrower than their positive spin counterparts, while their fluctuations are larger. The weak jet from the non-rotating black hole is the widest with the smallest fluctuations. Our results highlight the complex, non-linear connection between the black hole spin and the resulting disk and jet properties in the MAD regime.

Thierry Serafin Nadeau, Marten H. van Kerkwijk, Cees G. Bassa, Ben W. Stappers, Mitchell B. Mickaliger, Andrew G. Lyne

11 pages, 7 figures, submitted to ApJ

Using archival data from the 42 foot telescope at the Jodrell Bank Observatory, we produce daily stacks of aligned giant pulses for the Crab pulsar, to study changes to the daily profiles between April 2012 to December 2016. From these, we identify echoes, where intervening material away from the line of sight causes pulsed emission to be redirected towards the observer, with delay corresponding to the increased distance of travel, resulting in additional profile components. These observations show that such echoes may be far more common than implied by the previous rate of detections. All the observed echoes are consistent with approaching zero-delay at their closest approach to the normal giant pulse emission. This indicates that the structures responsible for producing these events must be highly anisotropic, with typical lengths greater than $\sim 4\textrm{AU}$, typical widths on the sky of $\sim 0.1 \textrm{AU}$ and typical depths of $\sim 5\textrm{AU}$, given the previously observed electron densities of the nebular filaments, on the order of 1000 cm$^{-3}$. This suggests that these inhomogeneities are likely to be offshoot substructure from the larger nebular filaments of the Crab nebula.

Elizabeth J. Iles, Alex R. Pettitt, Takashi Okamoto, Daisuke Kawata

18 pages, 12 figures, submitted to the Monthly Notices of the Royal Astronomical Society

Different mechanisms driving bar structure formation indicate that bar origins should be distinguishable in the stellar populations of galaxies. To study how these origins affect different bar morphologies and impact stellar orbits and migration, we analyse three simulated discs which are representative of bar formation under isolated evolution motivated by disc instability, and interaction driven tidal development. The first isolated disc and the tidally driven disc produce similar bar structure, while the second isolated disc, generated by the tidal initial condition without the companion, is visibly dissimilar. Changes to radial and vertical positions, angular momentum in the disc-plane, orbital eccentricity and the subsequent disc metallicities are assessed, as is the dependence on stellar age and formation radii. Bar origin is distinguishable, with the tidal disc displaying larger migration overall, higher metallicity difference between the inner and outer disc, as well as a population of inner disc stars displaced to large radii and below the disc-plane. The affect of closest approach on populations of stars formed before, after and during this period is evident. However, bar morphology is also found to be a significant factor in the evolution of disc stellar properties, with similar bars producing similar traits in migration tendency with radius, particularly in vertical stellar motion and in the evolution of central metallicity features.

Luc Dessart, Taeho Ryu, Pau Amaro Seoane, Andrew M. Taylor

submitted to A&A

High-velocity stellar collisions driven by a supermassive black hole (BH) or BH-driven disruptive collisions, in dense, nuclear clusters can rival the energetics of supergiant star explosions following gravitational collapse of their iron core. Here, starting from a sample of red-giant star collisions simulated with the hydrodynamics code AREPO, we generate photometric and spectroscopic observables using the nonlocal thermodynamic equilibrium time-dependent radiative transfer code CMFGEN. Collisions from more extended giants or stronger collisions (higher velocity or smaller impact parameter) yield bolometric luminosities on the order of 1e43 erg/s at 1d, evolving on a timescale of a week to a bright plateau at ~1e41 erg/s, before plunging precipitously after 20-40d at the end of the optically-thick phase. This luminosity falls primarily in the UV in the first days, thus when it is at its maximum, and shifts to the optical thereafter. Collisions at lower velocity or from less extended stars produce ejecta that are fainter but may remain optically thick for up to 40d if they have a small expansion rate. These collision debris show a similar spectral evolution as that observed or modeled for blue-supergiant star explosions of massive stars, differing only in the more rapid transition to the nebular phase. Such BH-driven disruptive collisions should be detectable by high-cadence surveys in the UV like ULTRASAT.

The paper presents numerical models of mean-motion resonances detected in the Quadrantid meteoroid stream consisting of particles of mass 0.003 to 0.03 g, which helps to prove the presence of the following mean-motion resonances: the 1:9 resonance with Venus, 1:5 resonance with Earth, 1:3 and 3:8 resonances with Mars, and 2:1,7:3, 9:4 resonances with Jupiter. Resonance particles create dust trails far from the Earth orbit. Indeed, there is no observational support for resonant effects in the observed Quadrantid meteoroid stream.

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) experiment is a pathfinder mission for future space-based instruments targeting the fluxes of Ultra-High Energy Cosmic Rays (UHECR), with energies exceeding 1EeV and very high energy diffuse and transient neutrinos, with energies exceeding 1PeV. Using two telescope designs: the Fluorescence Telescope (FT) and the Cherenkov Telescope (CT), EUSO-SPB2 made novel observations of the backgrounds relevant for space-based detection. EUSO-SPB2 will launch from Wanaka, NZ in Spring of 2023, for a long duration (up to 100d) flight at a nominal float altitude of 33km. In this contribution, we will focus on the CT's capability to measure cosmic rays from above Earth's limb via the Cherenkov emission produced by the resultant Extensive Air Showers (EAS). Using the EASCherSim optical Cherenkov generation code, we provide an updated estimate of the event rate of above-the-limb cosmic rays for the CT, taking into account updated values for the trigger efficiency as determined during the field testing of the instrument. We take particular care to consider the longitudinal development of EAS in rarefied atmosphere, accounting for the energy dependent elongation rate. In addition, we consider improvements to the magnetic field modeling present in EASCherSim and illustrate their impact on the observed events and detection thresholds. Finally, we compare these simulations to preliminary flight data from EUSO-SPB2.

Mohini Pachchigar, Yvette Perrott, Tulasi Parashar

6 pages, 4 figures, to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Mass is the most fundamental property of galaxy clusters. However, measuring it is still a challenge. Calibrating mass from intracluster medium observables such as the Sunyaev-Zel'dovich (SZ) effect is subject to uncertainty and biases because of the hydrostatic equilibrium assumption. On the other hand, merging cluster systems have been shown to exhibit radio emission which implies a link with disturbances from hydrostatic equilibrium. We present work on studying deviations of galaxy cluster gas pressure profile from the average (universal) pressure profile using an example of galaxy cluster Abell 1413 with SZ effect data from the Arcminute Microkelvin Imager and Planck. This cluster has also been observed at low radio frequency with the Murchison Widefield Array allowing the investigation of links between gas pressure profile deviations and the presence of radio emission.

K. Whitley, A. Kuznetsova, K. Gültekin, M. Ruszkowski

12 pages, 8 figures, submitted to MNRAS

We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disk, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of a SMBHB with mass ratio $q=0.01$ and separation $a=100\;R_g$, using a physically self-consistent steady disk model. We estimate the electromagnetic variability from the system by monitoring accretion onto the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disk, previously only observed for larger mass ratios, which we attribute to our disk model being much thinner ($H/R\approx0.01$ near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion onto the secondary SMBH on a timescale matching the orbital period of the binary, $t_{\rm{bin}}\approx0.1\;\rm{yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by $0.43\;t_{\rm{bin}}$. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and UV monitoring of AGN, rather than single-waveband periodicity alone.

Zhang XueGuang (GXU)

17 pages, 2 tables, 9 figures, accepted to be published in MNRAS

In the manuscript, effects of Tidal Disruption Events (TDEs) are estimated on long-term AGN variability, to provide interesting clues to detect probable hidden TDEs in normal broad line AGN with apparent intrinsic variability which overwhelm the TDEs expected variability features, after considering the unique TDEs expected variability patterns. Based on theoretical TDEs expected variability plus AGN intrinsic variability randomly simulated by Continuous AutoRegressive process, long-term variability properties with and without TDEs contributions are well analyzed in AGN. Then, interesting effects of TDEs can be determined on long-term observed variability of AGN. First, more massive BHs, especially masses larger than $10^7{\rm M_\odot}$, can lead to more sensitive and positive dependence of $\tau_{TN}$ on $R_{TN}$, with $\tau_{TN}$ as variability timescale ratio of light curves with TDEs contributions to intrinsic light curves without TDEs contributions, and $R_{TN}$ as ratio of peak intensity of TDEs expected variability to the mean intensity of intrinsic AGN variability without TDEs contributions. Second, stronger TDEs contributions $R_{TN}$ can lead to $\tau_{TN}$ quite larger than 5. Third, for intrinsic AGN variability having longer variability timescales, TDEs contributions will lead $\tau_{TN}$ to be increased more slowly. The results actually provide an interesting forward-looking method to detect probable hidden TDEs in normal broad line AGN, due to quite different variability properties, especially different DRW/CAR process expected variability timescales, in different epochs, especially in normal broad line AGN with shorter intrinsic variability timescales and with BH masses larger than $10^7{\rm M_\odot}$.

Vadim M. Uritsky, James A. Klimchuk

21 pages, 9 figures, 2 tables

We report results of an in-depth numerical investigation of three-dimensional projection effects which could influence the observed loop-like structures in an optically thin solar corona. Several archetypal emitting geometries are tested, including collections of luminous structures with circular cross-sections of fixed and random size, light-emitting structures with highly anisotropic cross-sections, as well as two-dimensional stochastic current density structures generated by fully-developed magnetohydrodynamic turbulence. A comprehensive set of statistical signatures is used to compare the line of sight -integrated emission signals predicted by the constructed numerical models with the loop profiles observed by the extreme ultraviolet telescope onboard the flight 2.1 of the High-Resolution Coronal Imager (Hi-C). The results obtained for the Hi-C loops cannot be attributed to randomly oriented quasi-two dimensional emitting structures such as those constituting "coronal veils" (Malanushenko et al., 2022), and they indicate that typical cross-sectional envelopes of loop emission cannot have high eccentricity. The possibility of apparent loop-like projections of very small (close to the resolution limit) or very large (comparable with the size of an active region) light-emitting sheets remains open, but the intermediate range of scales commonly associated with observed loop systems is most likely filled with true quasi-one dimensional (roughly axisymmetric) structures embedded into the three-dimensional coronal volume.

Graham S. Kerr, Adam F. Kowalski, Joel C. Allred, Adrian N. Daw, Melissa R. Kane

Accepted for publication in Monthly Notices of the Royal Astronomical Society

The bulk of solar flare energy is deposited in the chromosphere. Flare ribbons and footpoints in the chromosphere therefore offer great diagnostic potential of flare energy release and transport processes. High quality observations from the IRIS spacecraft have transformed our view of the Sun's atmospheric response to flares. Since most of the chromospheric lines observed by IRIS are optically thick, forward modelling is required to fully appreciate and extract the information they carry. Reproducing certain aspects of the Mg II lines remain frustratingly out of reach in state-of-the-art flare models, which are unable to satisfactorily reproduce the very broad line profiles. A commonly proposed resolution to this is to assert that very large values of `microturbulence' is present. We asses the validity of that approach by analysing optically thin lines in the flare chromosphere from the X-class flare SOL2014-10-25T17:08:00, using the derived value of nonthermal width as a constraint to our numerical models. A nonthermal width of the order 10~km~s$^{-1}$ was found within the short-lived red wing components of three spectral lines, with relatively narrow stationary components. Simulations of this flare were produced, and in the post-processing spectral synthesis we include within the downflows a microturbulence of 10~km~s$^{-1}$. While we can reproduce the O I 1355.598~\AA\ line rather well, and we can capture the general shape and properties of the Mg II line widths, the synthetic lines are still too narrow.

Mark Bishop, Yvette Perrott, Tulasi Parashar, Sean Oughton

to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Accurate estimation of galaxy cluster masses is a central problem in cosmology. Turbulence is believed to introduce significant deviations from the hydrostatic mass estimates. Estimation of turbulence properties is complicated by projection of the 3D cluster onto the 2D plane of the sky, and is commonly done in the form of indirect probes from fluctuations in the X-ray surface brightness and Sunyaev-Zeldovich effect maps. In this paper, we address this problem using simulations. We examine different methods for estimating the power spectrum on 2D projected fluctuation data, emulating data projected onto a 2D plane of the sky, and comparing them to the original, expected 3D power spectrum. Noise can contaminate the power spectrum of ICM observations, so we also briefly compare a few methods of reducing noise in the images for better spectral estimation.

Andrew Eberhardt, Alvaro Zamora, Michael Kopp, Tom Abel

to be submitted to phys rev d

The classical field approximation is widely used to better understand the predictions of ultra-light dark matter. Here, we use the truncated Wigner approximation method to test the classical field approximation of ultra-light dark matter. This method approximates a quantum state as an ensemble of independently evolving realizations drawn from its Wigner function. The method is highly parallelizable and allows the direct simulation of quantum corrections and decoherence times in systems many times larger than have been previously studied in reference to ultra-light dark matter. Our study involves simulation of systems in 1, 2, and 3 spatial dimensions. We simulate three systems, the condensation of a Gaussian random field in three spatial dimensions, a stable collapsed object in three spatial dimensions, and the merging of two stable objects in two spatial dimensions. We study the quantum corrections to the classical field theory in each case. We find that quantum corrections grow exponentially during nonlinear growth with the timescale being approximately equal to the system dynamical time. In stable systems the corrections grow quadratically. We also find that the primary effect of quantum corrections is to reduce the amplitude of fluctuations on the deBroglie scale in the spatial density. Finally, we find that the timescale associated with decoherence due to gravitational coupling to Baryonic matter is at least as fast as the quantum corrections due to gravitational interactions. These results strongly imply that quantum corrections do not impact the predictions of the classical field theory.

Yvette Perrott

6 pages, 5 figures, to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

At temperatures above ~5 keV, the non-relativistic approximation used to derive the classical thermal Sunyaev-Zel'dovich effect spectrum begins to fail. When relativistic effects are included, the spectrum becomes temperature-dependent. This leads to both a problem and an opportunity: a problem, because when the temperature dependence is not accounted for the Compton-y estimate is biased; and an opportunity, because it represents a new way to measure the temperature of the intracluster medium independently of X-ray observations. This work presents current results from investigating the impact of relativistic effects on Planck cluster observations, and projections for future measurements of cluster temperatures using the Atacama Large Aperture Sub-millimetre Telescope.

Shunke Ai, He Gao, Yong Yuan, Bing Zhang, Lin Lan

14 pages, 10 figures, accepted for publication on MNRAS

The maximum mass of neutron stars ($M_{\rm TOV}$) plays a crucial role in understanding their equation of state (EoS). Previous studies have used the measurements for the compactness of massive pulsars and the tidal deformability of neutron stars in binary neutron star (BNS) mergers to constrain the EoS and thus the $M_{\rm TOV}$. The discovery of the most massive pulsar, PSR J0952-0607, with a mass $\sim 2.35M_{\odot}$, has provided a valuable lower limit for $M_{\rm TOV}$. Another efficient method to constrain $M_{\rm TOV}$ is by examining the type of central remnant formed after a BNS merger. Gravitational wave (GW) data can provide the total mass of the system, while accompanying electromagnetic signals can help infer the remnant type. In this study, we combine all the previous constraints and utilize the observational facts that about $24\%$ of the short gamma-ray bursts are followed by an X-ray internal plateau, which indicate that roughly this fraction of BNS mergers yield supermassive neutron stars, to perform (Markov Chain) Monte Carlo simulations. These simulations allow us to explore the probability density distribution of $M_{\rm TOV}$ and other parameters related to BNS mergers. Our findings suggest that $M_{\rm TOV}$ is likely around $2.49M_{\odot} - 2.52M_{\odot}$, with an uncertainty range of approximately [$-0.16M_{\odot}$, $0.15M_{\odot}$] ([$-0.28M_{\odot}$, $0.26M_{\odot}$]) at $1\sigma$ ($2\sigma$) confidence level. Furthermore, we examine the type of merger remnants in specific events like GW170817 and GW190425 to further constrain $M_{\rm TOV}$ and other relevant parameters, which can help to understand the physical processes involved in BNS mergers.

Steph Sallum, Andrew Skemer, Deno Stelter, Ravinder Banyal, Natalie Batalha, Natasha Batalha, Geoff Blake, Tim Brandt, Zack Briesemeister, Katherine de Kleer, Imke de Pater, Aditi Desai, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Rebecca Jensen-Clem, Isabel Kain, Charlie Kilpatrick, Renate Kupke, Mackenzie Lach, Michael C. Liu, Bruce Macintosh, Raquel A. Martinez, Dimitri Mawet, Brittany Miles, Caroline Morley, Diana Powell, Ramya Sethuram, Patrick Sheehan, Justin Spilker, Jordan Stone, Arun Surya, Sivarani Thirupathi, Athira Unni, Kevin Wagner, Yifan Zhou

10 pages, 16 figures, submitted to Proceedings of the SPIE

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~\mu$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~\mu$m) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES' unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements, and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in $\sim2025$.

J. L. Qiu, H. Tong, H. G. Wang

accepted in ApJ. significantly rewritten compared with the original version

The rotating vector model and radius-to-frequency mapping in the presence of multipole magnetic field in pulsars and magnetars are considered. An axisymmetric potential field is assumed. It is found that: (1) The radiation beam in the case of multipole field is wider than the dipole case. This may account the increasing pulse width at higher frequency of pulsars (anti-radius-to-frequency mapping). (2) The expression for the polarization position angle is unchanged. Only the inclination angle {\alpha} and phase constant {\phi}_0 will change. The angle between the rotational axis and line of sight, and the position angle constant {\psi}_0 will not change. When fitting the varying position angle of magnetars, these constraints should be considered. The appearance and disappearance of multipole field may account for the changing slope of position angle in the radio emitting magnetar Swift J1818.0-1607. Similar but more active process in magnetar magnetospheres may account for the diverse position angle in fast radius bursts.

Mayu Tominaga (1, 2), Masahiro Tsujimoto (1), Ken Ebisawa (1), Teruaki Enoto (3), Kimitake Hayasaki (4, 5) ((1) Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, (2) Department of Astronomy, Graduate School of Science, The University of Tokyo, (3) Division of Physics and Astronomy, Graduate School of Science, Kyoto University, (4) Department of Space Science and Astronomy, Chungbuk National University, (5) Harvard-Smithsonian Center for Astrophysics)

Accepted for publication in ApJ, 16 pages, 7 figures

Circinus X-1 (Cir X-1) is a neutron star binary with an elliptical orbit of 16.6~days. The source is unique for its extreme youth, providing a key to understanding early binary evolution. However, its X-ray variability is too complex to reach a clear interpretation. We conducted the first high cadence (every 4 hours on average) observations covering one entire orbit using the NICER X-ray telescope. The X-ray flux behavior can be divided into stable, dip, and flaring phases. The X-ray spectra in all phases can be described by a common model consisting of a partially covered disk black body emission and the line features from a highly-ionized photo-ionized plasma. The spectral change over the orbit is attributable to rapid changes of the partial covering medium in the line-of-sight and gradual changes of the disk black body emission. Emission lines of the H-like and He-like Mg, Si, S, and Fe are detected, most prominently in the dip phase. The Fe emission lines change to absorption in the course of the transition from the dip phase to the flaring phase. The estimated ionization degree indicates no significant changes, suggesting that the photo-ionized plasma is stable over the orbit. We propose a simple model in which the disk black body emission is partially blocked by such a local medium in the line-of-sight that has spatial structures depending on the azimuth of the accretion disk. Emission lines upon the continuum emission are from the photo-ionized plasma located outside of the blocking material.

Yunpeng Jin, Ling Zhu, Stefano Zibetti, Luca Costantin, Glenn van de Ven, Shude Mao

20 pages, 18 figures

We employ a recently-developed population-orbit superposition technique to simultaneously fit the stellar kinematic and age maps of 82 CALIFA spiral galaxies, and obtain the ages of stars in different dynamical structures. We first evaluate the capabilities of this method on CALIFA-like mock data created from the Auriga simulations. The recovered mean ages of dynamically cold, warm and hot components match the true values well, with up to $20\%$ observational error in the mock age maps. For CALIFA spiral galaxies, we find that the stellar ages of the cold, warm and hot components all increase with galaxies' stellar mass, from $\overline{t_{\rm cold}}\sim2.2$ Gyr, $\overline{t_{\rm warm}}\sim2.3$ Gyr and $\overline{t_{\rm hot}}\sim2.6$ Gyr for galaxies with stellar mass $M_*<10^{10}\,\rm M_{\odot}$, to $\overline{t_{\rm cold}}\sim4.0$ Gyr, $\overline{t_{\rm warm}}\sim5.1$ Gyr and $\overline{t_{\rm hot}}\sim5.9$ Gyr for galaxies with $M_*>10^{11}\,\rm M_{\odot}$. About $80\%$ of the galaxies in our sample have $t_{\rm hot}>t_{\rm cold}$, and the mean values of $t_{\rm hot}-t_{\rm cold}$ also increase with stellar mass, from $0.7_{-0.2}^{+0.6}$ Gyr in low-mass galaxies ($10^{8.9}\,\rm M_{\odot}<M_*\le10^{10.5}\,\rm M_{\odot}$) to $1.7_{-0.2}^{+0.7}$ Gyr in high-mass galaxies ($10^{10.5}\,\rm M_{\odot}<M_*<10^{11.3}\,\rm M_{\odot}$). The average younger stellar age in disks than in bulges suggests that either disks formed later and/or experienced a more prolonged and extensive period of star formation. The lower mass spiral galaxies have younger bulges and younger disks, while higher mass spiral galaxies generally have older bulges and their disks span a wide range of ages. This is consistent with the scenario that the bulges in more massive spirals formed earlier than those in less massive spirals.

Subrata Kumar Panda, Siddharth Dhanpal, Simon J. Murphy, Shravan Hanasoge, Timothy R. Bedding

Asteroseismology is a powerful tool to probe the structure of stars. Space-borne instruments like CoRoT, Kepler and TESS have observed the oscillations of numerous stars, among which {\delta} Scutis are particularly interesting owing to their fast rotation rates and complex pulsation mechanisms. In this work, we inferred model-dependent masses, metallicities and ages of 60 {\delta} Scuti stars from their photometric, spectroscopic and asteroseismic observations using least-squares minimization. These statistics have the potential to explain why only a tiny fraction of {\delta} Sct stars pulsate in a very clean manner. We find most of these stars with masses around 1.6 {M_\odot} and metallicities below Z = 0.010. We observed a bimodality in age for these stars, with more than half the sample younger than 30 Myr, while the remaining ones were inferred to be older, i.e., hundreds of Myrs. This work emphasizes the importance of the large-frequency separation ({\Delta\nu}) in studies of {\delta} Scuti stars. We also designed three machine learning (ML) models that hold the potential for inferring these parameters at lower computational cost and much more rapidly. These models further revealed that constraining dipole modes can help in significantly improving age estimation and that radial modes succinctly encode information pertaining to stellar luminosity and temperature. Using the ML models, we also gained qualitative insight into the importance of stellar observables in estimating mass, metallicity, and age. The effective surface temperature T eff strongly affects the inference of all structure parameters and the asteroseismic offset parameter {\epsilon} plays an essential role in the inference of age.

W.-J. Kim, J. S. Urquhart, V. S. Veena, G. A. Fuller, P. Schilke, K-T Kim

23 pages, 16 figures, accepted for A&A for publication

The application of silicon monoxide (SiO) as a shock tracer arises from its propensity to occur in the gas phase as a result of shock-induced phenomena, including outflow activity and interactions between molecular clouds and expanding HII regions or supernova remnants. We searched for indications of shocks toward 366 massive star-forming regions by observing the ground rotational transition of SiO ($v=0$, $J=1-0$) at 43 GHz with the Korean VLBI Network (KVN) 21 m telescopes to extend our understanding on the origins of SiO in star-forming regions. We detected SiO emission toward 104 regions that consist of 57 IRDCs, 21 HMPOs, and 26 UCHIIs. The determined median SiO column density, $N$(SiO), and abundance, $X$(SiO), relative to $N$(H$_2$) are $8.12\times10^{12}$ cm$^{-2}$ and $1.28\times10^{-10}$, respectively. These values are similar to those obtained toward other star-forming regions and also consistent with predicted values from shock models with low-velocity shocks ($\lesssim$10 - 15 km s$^{-1}$). While the $X$(SiO) does not exhibit any strong correlation with the evolutionary stages of their host clumps, $L_{\rm SiO}$ is highly correlated with dust clump mass, and $L_{\rm SiO}/L_{\rm bol}$ also has a strong negative correlation with $T_{\rm dust}$. This shows that colder and younger clumps have high $L_{\rm SiO}/L_{\rm bol}$ suggestive of an evolutionary trend. This trend is not due to excess emission at higher velocities, such as SiO wing features, as the colder sources with high $L_{\rm SiO}/L_{\rm bol}$ ratios lack wing features. Comparing SiO emission with H$_2$O and Class I CH$_3$OH masers, we find a significant correlation between $L_{\rm SiO}$/$L_{\rm bol}$ and $L_{\rm CH_3OH}/L_{\rm bol}$ ratios, whereas no similar correlation is seen for the H$_2$O maser emission. This suggests a similar origin for the SiO and Class I CH$_3$OH emission in these sources.

Supriya Pan, Weiqiang Yang

Invited chapter for the edited book "Hubble Constant Tension" (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024); comments are welcome

The Hubble constant $H_0$ is one of the important cosmological parameters measuring the expansion rate of our universe at present moment. Over the last couple of years, $H_0$ has created an enormous amount of debates interests in the astrophysical and cosmological communities for its different estimations at many standard deviations by different observational surveys. The recent estimation of $H_0$ from the cosmic microwave background observations by Planck within the $\Lambda$-Cold Dark Matter ($\Lambda$CDM) paradigm is in tension with $\gtrsim 5 \sigma$ confidence with SH0ES (Supernovae and $H_0$ for the Equation of State of dark energy) collaboration. As a result, revision of the standard $\Lambda$CDM model has been suggested in various ways in order to examine whether such scenarios can solve this $H_0$ tension. Among the list of the proposed cosmological scenarios, in this chapter we focus on a generalized cosmological theory in which the dark components of the universe, namely, Dark Matter (DM) and Dark Energy (DE) are allowed to interact with each other in a non-gravitational way, widely known as the Interacting DE or Coupled DE scenarios. These interacting scenarios have received magnificent attention in the scientific community for their appealing consequences. Specifically, in the context of $H_0$ tension, it has been observed that the interacting DE scenarios can lead to higher values of the Hubble constant ($H_0$) value, and consequently, the tension on $H_0$ can be either alleviated or solved. In this chapter we review various interacting DE scenarios and their roles in alleviating the $H_0$ tension.

Giunchi, E, Poggianti, B. M., Gullieuszik, M., Moretti, A., Werle, A., Zanella, A., Vulcani, B., Tonnesen, S., Calzetti, D., Bellhouse, C., Scarlata, C., Bacchini, C

20 pages, 14 figures, 2 tables. Accepted for publication in ApJ

We characterize the morphological properties of a statistically relevant sample of H$\alpha$ and UV young star-forming clumps and optical complexes, observed with the \textit{Hubble Space Telescope} in six galaxies of the GASP sample undergoing ram-pressure stripping. The catalogs comprise 2406 (323 in the tails) H$\alpha$ clumps, 3750 (899) UV clumps and 424 tail optical complexes. About 15-20\% of the clumps and 50\% of the complexes are resolved in size. We find that more than half of the complexes contain no H$\alpha$ clumps, while most of them contain at least one UV clump. The clump number and size increase with the complex size, while the median complex filling factor is larger for UV clumps ($0.27$) than for H$\alpha$ clumps ($0.10$) and does not correlate with almost any morphological property. This suggests that the clumps number and size grow with the complex keeping the filling factor constant. When studying the position of the clumps inside their complexes, H$\alpha$ clumps, and UV clumps to a lesser extent, show a displacement from the complex center of $0.1-1$ kpc and, in $\sim 60$\% of the cases, they are displaced away from the galactic disk. This is in accordance with the fireball configuration, already observed in the tails of stripped galaxies. Finally, the filling factor and the clump radius increase with the distance from the galactic disk, suggesting that the reciprocal displacement of the different stellar generations increases as a consequence of the velocity gradient caused by ram pressure.

Kiran Jain, Sushanta C. Tripathy

6 pages and 5 figures; Proceedings of the IAU Symposium 365 - Dynamics of Solar and Stellar Convection Zones and Atmospheres, 2023 August 21-25, Yerevan, Armenia

We present analysis of the evolution of subsurface flows in and around active regions with peculiar magnetic configurations and compare their characteristics with the normal active regions. We also study the zonal and meridional components of subsurface flows separately in different polarity regions separately to better understand their role in flux migration. We use the techniques of local correlation tracking and ring diagrams for computing surface and subsurface flows, respectively. Our study manifests an evidence that the meridional component of the flows near anti-Hale active regions is predominantly equatorward which disagrees with the poleward flow pattern seen in pro-Hale active regions. We also find clockwise or anti-clockwise flows surrounding the anti-Joy active regions depending on their locations in the Southern or Northern hemispheres, respectively.

Yuchen Liu, Linhua Jiang, Rogier A. Windhorst, Yucheng Guo, Zhenya Zheng

13 pages, 9 figures. Accepted for publication in ApJ

We present a study of Lyman continuum (LyC) emission in a sample of $\sim$150 Ly$\alpha$ emitters (LAEs) at $z\approx3.1$ in the Subaru-XMM Deep Survey field. These LAEs were previously selected using the narrowband technique and spectroscopically confirmed with Ly$\alpha$ equivalent widths (EWs) $\ge45$ \r{A}. We obtain deep UV images using a custom intermediate-band filter $U_{\rm J}$ that covers a wavelength range of $3330 \sim 3650$ \r{A}, corresponding to 810$\sim$890 \r{A} in the rest frame. We detect 5 individual LyC galaxy candidates in the $U_{\rm J}$ band, and their escape fractions ($f_{\rm esc}$) of LyC photons are roughly between 40% and 80%. This supports a previous finding that a small fraction of galaxies may have very high $f_{\rm esc}$. We find that the $f_{\rm esc}$ values of the 5 LyC galaxies are not apparently correlated with other galaxy properties such as Ly$\alpha$ luminosity and EW, UV luminosity and slope, and star-formation rate (SFR). This is partly due to the fact that these galaxies only represent a small fraction ($\sim3$%) of our LAE sample. For the remaining LAEs that are not detected in $U_{\rm J}$, we stack their $U_{\rm J}$-band images and constrain their average $f_{\rm esc}$. The upper limit of the average $f_{\rm esc}$ value is about 16%, consistent with the results in the literature. Compared with the non-LyC LAEs, the LyC LAEs tend to have higher Ly$\alpha$ luminosities, Ly$\alpha$ EWs, and SFRs, but their UV continuum slopes are similar to those of other galaxies.

Weiyang Liu, Linhua Jiang

16 pages, 13 figures, 4 tables (+4 pages, 8 figures in Appendix). Accepted for publication in ApJ

We present a spectroscopic survey of Ly$\alpha$ emitters (LAEs) at $z\sim3.7$ and $z\sim4.8$. The LAEs are selected using the narrowband technique based on the combination of deep narrowband and broadband imaging data in two deep fields, and then spectroscopically confirmed with the MMT multi-fiber spectrograph Hectospec. The sample consists of 71 LAEs at $z\sim3.7$ and 69 LAEs at $z\sim4.8$ over $\sim 1.5$ deg$^2$, making it one of the largest spectroscopically confirmed sample of LAEs at the two redshifts. Their Ly$\alpha$ luminosities are measured using the secure redshifts and deep photometric data, and span a range of $\sim 10^{42.5}$ - $10^{43.6} \,\rm erg\, s^{-1}$, so these LAEs represent the most luminous galaxies at the redshifts in terms of Ly$\alpha$ luminosity. We estimate and correct sample incompletenesses and derive reliable Ly$\alpha$ luminosity function (LF)s at $z\sim3.7$ and 4.8 based on the two spectroscopic samples. We find that our Ly$\alpha$ LFs are roughly consistent (within a factor of $2-3$) with previous measurements at similar redshifts that were derived from either photometric samples or spectroscopic samples. By comparing with previous studies in different redshifts, we find that the Ly$\alpha$ LFs decrease mildly from $z\sim3.1$ to $z\sim5.7$, supporting the previous claim of the slow LF evolution between $z\sim2$ and $z\sim6$. At $z>5.7$, the LF declines rapidly towards higher redshift, partly due to the effect of cosmic reionization.

Masahiro Morikawa, Akika Nakamichi

11 pages, 6 figures, 1 table

We first study the solar flare time sequence based on the GOES16 data. We find that the power spectrum density of the low-energy (E\leq E_{mean}) flare shows 1/f fluctuations, but the high-energy (E>E_{mean}) flare shows a flat spectrum. Further, we found that the flare timing time-sequence shows 1/f fluctuations clearer. These facts indicate that the solar flare 1/f fluctuations are associated with low-energy phenomena. We investigate the origin of this 1/f fluctuations based on our recent proposal: 1/f fluctuations arise from amplitude modulation and demodulation. We speculate that this amplitude modulation is encoded by the resonance with the Solar Five-minute Oscillation (SFO) and demodulated by magnetic reconnection. We partially demonstrate this scenario by analyzing the SFO eigenmodes resolving the frequency degeneracy in the azimuthal order number m by solar rotation and resonance. Since 1/f fluctuation is robust, we speculate that the solar flare 1/f fluctuations may be inherited by the various phenomena around the Sun, such as the sunspot numbers and the cosmic rays. Finally, we compare the solar flares and the earthquakes, both showing 1/f fluctuations. Interestingly, the same analysis for solar flares is possible for earthquakes if we read SFO as Earth's Free Oscillation, and magnetic reconnection as fault rupture. Furthermore, we point out the possibility that the same analysis also applies to the activity of the black hole/disk system if we read SFO as the Quasi-Periodic Oscillation of a black hole.

Kala G. Pradeep, Kulinder Pal Singh, G. C. Dewangan, Elias Aydi, P. E. Barrett, D. A. H. Buckley, V. Girish, K. L. Page, S. B. Potter, E. M. Schlegel

Accepted for publication in MNRAS, 17 Pages, 13 Figures

We present multi-wavelength temporal and spectral characteristics of a magnetic cataclysmic variable (MCV) Swift J0503.7-2819, using far ultraviolet (FUV) and X-ray data from AstroSat, supplemented with optical data from the Southern African Large Telescope and X-ray data from the XMM-Newton and Swift observatories. The X-ray modulations at 4897.6657 s and 3932.0355 s are interpreted as the orbital ($P_{\Omega}$) and spin ($P_{\omega}$) period, respectively, and are consistent with prior reports. With a spin-orbit period ratio of 0.8 and $P_{\Omega}$ falling below the period gap (2-3 hrs) of CVs, Swift J0503.7-2819 would be the newest addition to the growing population of nearly synchronous MCVs, which we call EX Hya-like systems. Hard X-ray luminosity of $<$ $2.5\times10^{32} erg s^{-1}$, as measured with the Swift Burst Alert Telescope, identifies it to be a low-luminosity intermediate polar, similar to other EX Hya-like systems. The phenomenology of the light curves and the spectral characteristics rule out a purely disc-fed/stream-fed model and instead reveal the presence of complex accretion structures around the white dwarf. We propose a ring-like accretion flow, akin to EX Hya, using period ratio, stability arguments, and observational features. An attempt is made to differentiate between the asynchronous polar/nearly-synchronous intermediate polar nature of Swift J0503.7-2819. Further, we note that with the advent of sensitive surveys, a growing population of MCVs that exhibit characteristics of both polars and intermediate polars is beginning to be identified, likely forming a genealogical link between the two conventional classes of MCVs.

Zongnan Li, Zhiyuan Li, Ruben Garcia-Benito, Yifei Jin

27 pages, 16 figures, accepted for publication in ApJ. Comments welcome

The origin of the low-ionization nuclear emission-line region (LINER) prevalent in local galaxies and its relationship with supermassive black holes are debated for decades. We preform a comprehensive evaluation of traditional photoionization models against the circumnuclear ionized gas in M81, for which recent CAHA/PPAK integral-field spectroscopic observations reveal a LINER characteristic out to a galactocentric radius of ~1 kpc. Constructed with the photoionization code cloudy, the models have the novel aspect of their primary parameters being well constrained by extensive observations of a prototypical low-luminosity active galactic nucleus (LLAGN) and an old stellar bulge in M81. Additionally, these models incorporate a reasonably broad range of uncertain nebular properties. It is found that the integrated photoionization by the LLAGN and hot, low-mass stars distributed in the bulge can roughly reproduce the observed radial intensity distributions of the Halpha, Hbeta and [N ii] lines, with the bulge stars dominating the ionizing flux at radii >~200 pc. However, the models generally fail to reproduce a similarly declining profile of the [Oiii] line or an accordingly flat profile of the [O iii]/Hbeta ratio. This clearly points to a deficiency of ionizing photons in the outer regions despite an extended photoionization source. The discrepancy might be alleviated if much of the observed [O iii] line arose from a bulge-filling, low-density gas surrounding a denser, Halpha-emitting disk, or by a higher AGN luminosity in the recent past. The case of M81 has important implications for the ionization mechanism of LINERs and low-ionization emission-line regions in general.

Thomas Pfeil, Til Birnstiel, Hubert Klahr

27 pages, 15 figures, Accepted for publication in The Astrophysical Journal

Protoplanetary disks exhibit a vertical gradient in angular momentum, rendering them susceptible to the Vertical Shear Instability (VSI). The most important condition for the onset of this mechanism is a short timescale of thermal relaxation ($\lesssim 0.1$ orbital timescales). Simulations of fully VSI active disks are characterized by turbulent, vertically extended dust layers. This is in contradiction with recent observations of the outer regions of some protoplanetary disks, which appear highly settled. In this work, we demonstrate that the process of dust coagulation can diminish the cooling rate of the gas in the outer disk and extinct the VSI activity. Our findings indicate that the turbulence strength is especially susceptible to variations in the fragmentation velocity of the grains. A small fragmentation velocity of $\approx 100 \mathrm{\, cm \,s^{-1}}$ results in a fully turbulent simulation, whereas a value of $\approx 400 \mathrm{\, cm \,s^{-1}}$ results in a laminar outer disk, being consistent with observations. We show that VSI turbulence remains relatively unaffected by variations in the maximum particle size in the inner disk regions. However, we find that dust coagulation can significantly suppress the occurrence of VSI turbulence at larger distances from the central star.

Pedro Carrilho, Chiara Moretti, Maria Tsedrik

Contribution to the 34th Rencontres de Blois on Particle Physics and Cosmology (Blois 2023)

The current discrepancy between the CMB and weak lensing measurements of the amplitude of matter fluctuations, the so-called $S_8$ tension, has attracted a great deal of recent attention, as it may show a crack in the $\Lambda$CDM model of cosmology. We review the evidence for this tension and describe potential solutions, focusing on extensions of the standard cosmological model, including interacting dark energy and modified gravity. We present a likelihood analysis of the BOSS DR12 data, probing these alternative models as well as $\Lambda$CDM. From this analysis, we show hints of non-standard cosmology compatible with those seen in weak lensing observations, demonstrating that interacting dark energy or modified gravity can explain them successfully. We then discuss the robustness of these results to analysis choices, as well as future paths to confirm them with additional data and further distinguish between models.

Shulei Ni, Yichao Li, Xin Zhang

10 pages, 10 figures

The cosmic microwave background (CMB) provides an important source for the study of weak gravitational lensing. The lensed CMB can be used to constrain cosmological parameters, but it also has smoothing effect on the original CMB. The angular power spectrum of the unlensed CMB has sharper acoustic peaks and more prominent damping tails, allowing for improved inferences of cosmological parameters affecting these characteristics. Although delensing reduces the B-mode power spectrum, it aids the search for primordial gravitational waves and allows for lower variance reconstruction of lensing and other sources of secondary CMB anisotropies. In this work, we investigate whether the deep learning methods, concretely the U-Net++ algorithm here, can play a crucial role in CMB delensing. We consider three fields, i.e., T, Q and U sky maps, present the angular power spectra of CMB delensed TT, EE and BB, and compare them with the unlensed CMB angular power spectra. The results show that the angular power spectrum of the lensed CMB processed by U-Net++ is in great agreement with the unlensed CMB angular power spectrum. Therefore, the U-Net++ based CMB delensing can efficiently eliminate the effects caused by weak gravitational lensing and shed new light on the CMB power spectrum for future CMB experiments.

Hiroko Miyamoto, Matteo Battisti, Dario Barghini, Alexander Belov, Mario Bertaina, Marta Bianciotto, Francesca Bisconti, Carl Blaksley, Sylvie Blin, Karl Bolmgren, Giorgio Cambiè, Francesca Capel, Marco Casolino, Igor Churilo, Christophe De La taille, Toshikazu Ebisuzaki, Johannes Eser, Francesco Fenu, Geroge Filippatos, Massimo Alberto Franceschi, Christer Fuglesang, Alessio Golzio, Philippe Gorodetzky, Fumioshi Kajino, Hiroshi Kasuga, Pavel Klimov, Viktoria Kungel, Vladimir Kuznetsov, Massimiliano Manfrin, Laura Marcelli, Gabriele Mascetti, Włodzimierz Marszał, Marco Mignone, Alexey Murashov, Tommaso Napolitano, Hitoshi Ohmori, Angela Olinto, Etienne Parizot, Piergiorgio Picozza, Lech Wiktor Piotrowski, Zbigniew Plebaniak, Guillaume Prévôt, Enzo Reali, Marco Ricci, et al. (10 additional authors not shown)

Mini-EUSO is a wide Field-of-View (FoV, 44$^{\circ}$) telescope currently in operation from a nadia-facing UV-transparent window in the Russian Zvezda module on the International Space Station (ISS). It is the first detector of the JEM-EUSO program deployed on the ISS, launched in August 2019. The main goal of Mini-EUSO is to measure the UV emissions from the ground and atmosphere, using an orbital platform. Mini-EUSO is mainly sensitive in the 290-430 nm bandwidth. Light is focused by a system of two Fresnel lenses of 25 cm diameter each on the Photo- Detector-Module (PDM), which consists of an array of 36 Multi-Anode Photomultiplier Tubes (MAPMTs), for a total of 2304 pixels working in photon counting mode, in three different time resolutions of 2.5 ${\mu}$s, 320 ${\mu}$s, 40.96 ms operation in parallel. In the longest time scale, the data is continuously acquired to monitor the UV emission of the Earth. It is best suited for the observation of ground sources and therefore has been used for the observational campaigns of the Mini-EUSO. In this contribution, we present the assembled UV flasher, the operation of the field campaign and the analysis of the obtained data. The result is compared with the overall efficiency computed from the expectations which takes into account the atmospheric attenuation and the parameterization of different effects such as the optics efficiency, the MAPMT detection efficiency, BG3 filter transmittance and the transparency of the ISS window.

Teppei Okumura, Atsushi Taruya, Toshiki Kurita, Takahiro Nishimichi

27 pages, 9 figures, 3 tables

We present an analytic model of nonlinear correlators of galaxy/halo ellipticities in redshift space. The three-dimensional ellipticity field is not affected by the redshift-space distortion (RSD) at linear order, but by the nonlinear one, known as the Finger-of-God effect, caused by the coordinate transformation from real to redshift space. Adopting a simple Gaussian damping function to describe the nonlinear RSD effect and the nonlinear alignment model for the relation between the observed ellipticity and underlying tidal fields, we derive analytic formulas for the multipole moments of the power spectra of the ellipticity field in redshift space expanded in not only the associated Legendre basis, a natural basis for the projected galaxy shape field, but also the standard Legendre basis, conventionally used in literature. The multipoles of the corresponding correlation functions of the galaxy shape field are shown to be expressed by a simple Hankel transform, as is the case for those of the conventional galaxy density correlations. We measure these multipoles of the power spectra and correlation functions of the halo ellipticity field using large-volume N-body simulations. We then show that the measured alignment signals can be better predicted by our nonlinear model than the existing linear alignment model. The formulas derived here have already been used to place cosmological constraints using from the redshift-space correlation functions of the galaxy shape field measured from the Sloan Digital Sky Survey (Okumura and Taruya, 2023).

Hiroko Miyamoto, Matteo Battisti, Alexander Belov, Mario Bertaina, Sylvie Blin, Alxandre Creusot, Johannes Eser, George Filippatos, Pavel Klimov, Massimiliano Manfrin, Marco Mignone, Etienne Parizot, Lech Wiktor Piotrowski, Guillaume Prévôt, Daniil Trofimov (for the JEM-EUSO collaboration)

The second generation Extreme Universe Space Observatory on a Super-Pressure Balloon (EUSO-SPB2) mission is a stratospheric balloon mission developed within the Joint Exploratory Missions for Extreme Universe Space Observatory (JEM-EUSO) program. The Fluorescence Telescope (FT) is one of the two separate Schmidt telescopes of EUSO-SPB2, which aims at measuring the fluorescence emission of extensive air showers from cosmic rays above the energy of 1 EeV, looking downwards onto the atmosphere from the float altitude of 33 km. The FT measures photons with a time resolution of 1.05 $\mu$s in two different modes: single photon counting (PC) and charge integration (KI). In this paper, we describe the latter and report on the measurements of its characteristics. We also present a new trigger based on this channel, the so-called KI trigger, which allows to measure additional types of events, namely very short and intense light pulses. We report on the tests of this trigger mode in the laboratory and at the TurLab facility, and its implementation in the EUSO-SPB2 mission.

C. Hanser, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M.Calvo, A. Catalano, M. De Petris, F.-X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, A. Gomez, J. Goupy, S. Katsioli, F. Kéruzoré, C. Kramer, B. Ladjelate, G. Lagache, S. Leclercq, J.-F. Lestrade, J. F. Macías-Pérez, S. C. Madden, A. Maury, P. Mauskopf, F. Mayet, A. Monfardini, A. Moyer-Anin, M. Muñoz-Echeverría, C. Payerne, L. Perotto, G. Pisano, E. Pointecouteau, N. Ponthieu, G. W. Pratt, V. Revéret, A. J. Rigby, A. Ritacco, C. Romero, H. Roussel, F. Ruppin, K. Schuster, A. Sievers, C. Tucker

to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 38 galaxy clusters covering a wide mass range at intermediate to high redshift. The measured SZ fluxes will be essential to calibrate the SZ scaling relation and the galaxy clusters mean pressure profile, needed for the cosmological exploitation of SZ surveys. We present in this study a method to infer a mean pressure profile from cluster observations. We have designed a pipeline encompassing the map-making and the thermodynamical properties estimates from maps. We then combine all the individual fits, propagating the uncertainties on integrated quantities, such as $R_{500}$ or $P_{500}$, and the intrinsic scatter coming from the deviation to the standard self-similar model. We validate the proposed method on realistic LPSZ-like cluster simulations.

Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito, Katsuaki Asano, Jann Aschersleben, Halim Ashkar, Luiz Augusto Stuani, Dominik Baack, Michael Backes, Csaba Balazs, Matteo Balbo, Andres Baquero Larriva, Victor Barbosa Martins, Ulisses Barres de Almeida, Juan Abel Barrio, Denis Bastieri, Pedro Ivo Batista, Ivana Batković, Rowan William Batzofin, Joshua Ryo Baxter, Geoffrey Beck, Julia Becker Tjus, Lukas Beiske, Daniele Belardinelli, Wystan Benbow, Elisa Bernardini, Juan Bernete Medrano, et al. (521 additional authors not shown)

Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 ( arXiv:2309.08219 )

The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA.

Samuel T. Spencer, Alison M.W. Mitchell, Brian Reville

7 pages, 2 figures. From proceedings of the 38th International Cosmic Ray Conference (ICRC2023)

Recent results from LHAASO and Tibet AS$\gamma$ suggest that the Crab Nebula's gamma-ray spectrum extends to the PeV energy range, however the production mechanisms of this highest energy emission remain unclear. It has been postulated that a secondary component of hadronic emission could explain the highest energy gamma-ray flux points, however the origin and acceleration mechanism for this hadronic population has yet to be explained. We postulate one scenario in which hadrons diffuse over time into the Crab pulsar wind nebula from the surrounding supernova ejecta, and are subsequently re-accelerated by the pulsar wind termination shock. We present results of direct particle transport simulations (including radial evolution) to determine if this scenario is viable over the lifetime of the Crab system.

Matteo Lucchini, Marina Ten Have, Jingyi Wang, Jeroen Homan, Erin Kara, Oluwashina Adegoke, Riley Connors, Thomas Dauser, Javier Garcia, Guglielmo Mastroserio, Adam Ingram, Michiel van der Klis, Ole König, Collin Lewin, Labani Mallick, Edward Nathan, Patrick O'Neill, Christos Panagiotou, Joanna Piotrowska, Phil Uttley

Accepted for publication in ApJ

During the outbursts of black hole X-ray binaries (BHXRBs), their accretion flows transition through several states. The source luminosity rises in the hard state, dominated by non-thermal emission, before transitioning to the blackbody-dominated soft state. As the luminosity decreases, the source transitions back into the hard state and fades to quiescence. This picture does not always hold, as $\approx$ 40$\%$ of the outbursts never leave the hard state. Identifying the physics that govern state transitions remains one of the outstanding open questions in black hole astrophysics. In this paper we present an analysis of archival RXTE data of multiple outbursts of GX 339-4. We compare the properties of the X-ray variability and time-averaged energy spectrum and demonstrate that the variability (quantified by the power spectral hue) systematically evolves $\approx$ 10-40 days ahead of the canonical state transition (quantified by a change in spectral hardness); no such evolution is found in hard state only outbursts. This indicates that the X-ray variability can be used to predict if and when the hard-to-soft state transition will occur. Finally, we find a similar behavior in ten outbursts of four additional BHXRBs with more sparse observational coverage. Based on these findings, we suggest that state transitions in BHXRBs might be driven by a change in the turbulence in the outer regions of the disk, leading to a dramatic change in variability. This change is only seen in the spectrum days to weeks later, as the fluctuations propagate inwards towards the corona.

F. Rösch, P. Benke, M. Kadler, E. Ros, R. Ojha, P. G. Edwards, F. Eppel, J. Heßdörfer, J. Stevens

8 pages, 2 figures, 1 table, Proceedings of the 38th International Cosmic Ray Conference (ICRC2023)

The origin of high-energy cosmic neutrinos detected by the IceCube observatory is a hotly debated topic in astroparticle physics. There is growing evidence that some of these neutrinos can be associated with active galactic nuclei (AGN) and especially with blazars. Several recent studies have revealed a statistical correlation between radio-bright AGN samples and IceCube neutrino event catalogs. In addition, a growing number of individual high-energy neutrinos have been found to coincide with individual radio-flaring blazars. These observational results strongly call for high-quality, high angular-resolution radio observations of such neutrino-associated blazars to study their parsec-scale jet structures. TANAMI is the only large and long-term VLBI monitoring program focused on the Southern sky. Within TANAMI, we put an emphasis on Southern IceCube neutrino candidate blazars at 2.3 GHz and 8.4 GHz. Here we present first results of the first high-quality, high angular-resolution VLBI observations of nine Southern-Hemisphere blazars that were associated to IceCube neutrino hotspots in the Southern sky. In the near future, the rapidly growing KM3NeT will complement IceCube by being sensitive to high-energy neutrinos mainly from the Southern Hemisphere. This will increase the importance of Southern-Hemisphere radio monitoring programs of neutrino-associated blazars, like TANAMI.

Yu-Mei Wu, Zu-Cheng Chen, Yan-Chen Bi, Qing-Guo Huang

6 pages, 1 figure, 1 table

The recently detected stochastic signal by several pulsar timing array collaborations, offers an opportunity to scrutinize the fundamental properties of gravity, including the potential mass of the graviton. In this study, we analyze the NANOGrav 15-year data set to search for a stochastic gravitational wave background with modified Hellings-Downs correlations predicted by massive gravity. While the Bayesian analysis comparing the massive gravity to massless gravity within the effective searchable mass range of $m_g\in [3\times 10^{-25}, 8 \times 10^{-24}]\,\rm{eV}/c^2$ does not yield an explicit upper bound as all the Bayes factors are smaller than $3$, the combined consideration of the minimum frequency inherent in a massive gravity and the observed spectrum leads to an upper limit of $m_g<8.2\times 10^{-24}\,\rm{eV}/c^2$.

William J. Wolf, Pedro G. Ferreira

14 pages, 8 figures. Comments are welcomed!

There is compelling evidence that the Universe is undergoing a late phase of accelerated expansion. One of the simplest explanations for this behaviour is the presence of dark energy. A plethora of microphysical models for dark energy have been proposed. The hope is that, with the ever increasing precision of cosmological surveys, it will be possible to precisely pin down the model. We show that is unlikely and that, at best, we will have a phenomenological description for the microphysics of dark energy. Furthermore, we argue that the current phenomenological prescriptions are ill-equipped for shedding light on the fundamental theory of dark energy.

Sara Santoni, Marco De Petris, Antonio Ferragamo, Gustavo Yepes, Weiguang Cui

to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Galaxy clusters and their filamentary outskirts reveal useful laboratories to test cosmological models and investigate Universe composition and evolution. Their environment, in particular the filaments of the Cosmic Web to which they are connected, plays an important role in shaping the properties of galaxy clusters. In this project, we analyse the gas filamentary structures present in 324 regions of The Three Hundred hydrodynamical simulation extracted with the DisPerSE filament finder. We estimate the number of gas filaments globally connected to several galaxy clusters, i.e. the connectivity k, with a mass range of $10^{13} \leq M_{200} \, h^{-1} \, M_{\odot} \leq 10^{15} $ at redshift $z=0$. We study the positive correlation between the connectivity and mass of galaxy clusters. Moreover, we explore the impact of filaments on the dynamical state of clusters, quantified by the degree of relaxation parameter $\chi$.

Olivier Lai, Mark Chun, Stefan Kuiper, Niek Doelman, Marcel Carbillet, Mamadou N'Diaye, Frantz Martinache, Lyu Abe, Jean-Pierre Rivet, Dirk Schmidt

9 pages, 9 figures, Conference Proceeding AO4ELT7

Adaptive optics is a technique mostly used on large telescopes. It turns out to be challenging for smaller telescopes (0.5~2m) due to the small isoplanatic angle, small subapertures and high correction speeds needed at visible wavelengths, requiring bright stars for guiding, severely limiting the sky coverage. NGS SCAO is ideal for planetary objects but remains limited for general purpose observing. The approach we consider is a compromise between image quality gain and sky coverage: we propose to partially improve the image quality anywhere in the sky instead of providing the diffraction limit around a few thousand bright stars. We suggest a new solution based on multiple AO concepts brought together: The principle is based on a rotating Foucault test, like the first AO concept proposed by H. Babcock in 1953, on the Ground Layer Adaptive Optics, proposed by Rigaut and Tokovinin in the early 2000s, and on the idea of Layer-oriented MCAO and the pupil-plane wavefront analysis by R. Ragazzoni. We propose to combine these techniques to use all the light available in a large field to measure the ground layer turbulence and enable the high angular resolution imaging of regions of the sky (e.g., nebulas, galaxies) inaccessible to traditional AO systems. The motivation to develop compact and robust AO system for small telescopes is two-fold: On the one hand, universities often have access to small telescopes as part of their education programs. Also, researchers in countries with fewer resources could also benefit from reliable adaptive optics system on smaller telescopes for research and education purposes. On the other hand, amateur astronomers and enthusiasts want improved image quality for visual observation and astrophotography. Implementing readily accessible adaptive optics in astronomy clubs would also likely have a significant impact on citizen science.

Yongmin Yoon, Jae-Woo Kim, Jongwan Ko

20 pages, 14 figures, 1 table, accepted for publication in the ApJ

We study the mass-size relation of quiescent galaxies across various environments, with a particular focus on its environmental dependence at the low-mass part of $\log(M_\mathrm{star}/M_{\odot})\lesssim10.0$. Our sample consists of 13,667 quiescent galaxies with $\log(M_\mathrm{star}/M_{\odot})\ge9.4$ and $0.01<z<0.04$ from the Sloan Digital Sky Survey. We find that the mass-size relation of low-mass quiescent galaxies (LQGs) with $\log(M_\mathrm{star}/M_{\odot})\lesssim10.0$ depends on their environment, with LQGs in the highest-density environments exhibiting an average size $\sim70\%$ larger than those in isolated environments. Moreover, the slope of the mass-size relation for LQGs in high-density environments is significantly shallower than that of their counterparts in isolated environments. This is in contrast with high-mass quiescent galaxies with $\log(M_\mathrm{star}/M_{\odot})\gtrsim10.5$ that show a nearly identical mass-size relation across all environments. Combined with additional discoveries that the mass-size relation slopes of LQGs and star-forming galaxies are similar to each other in high-density environments, and that LQGs in higher-density environments exhibit more disk-like structures, our results support the idea that LQGs in high-density environments have evolved from star-forming galaxies through environmental effects, which are capable of causing their quenching and transformation into quiescent galaxies. With the aid of an analysis of merger rates for simulated galaxies from a cosmological galaxy formation simulation, we suggest that the steep slope and low normalization of the mass-size relation of LQGs in the lowest-density environments may originate from recent gas-rich mergers, which occur over 10-30 times more frequently in the progenitors of LQGs in the lowest-density environments than in their counterparts in high-density environments at low redshifts.

P. Mroz, R. Poleski

submitted, invited chapter for the "Handbook of Exoplanets"

The number of exoplanets detected using gravitational microlensing technique is currently larger than 200, which enables population studies. Microlensing is uniquely sensitive to low-mass planets orbiting at separations of several astronomical units, a parameter space that is not accessible to other planet-detection techniques, as well as free-floating planets, not orbiting around any star. In this review, we present the state-of-the-art knowledge on the demographics of exoplanets detected with microlensing, with a particular emphasis on their occurrence rates. We also summarize the current knowledge about free-floating planets, an elusive population of objects that seem to be more common than ordinary, gravitationally bound exoplanets.

Galaxy clusters could produce gamma-rays from inverse Compton scattering of cosmic ray electrons or hadronic interactions of cosmic ray protons with the intracluster medium. It is still an open question on whether gamma-ray emission ($>$ GeV energies) has been detected from galaxy clusters. We carry out a systematic search for gamma-ray mission based on 300 galaxy clusters selected from the 2500 deg.$^2$ SPT-SZ survey after sorting them in descending order of $M_{500}/z^2$, using about 15 years of Fermi-LAT data in the energy range between 1-300 GeV. We were able to detect gamma-ray emission with significance of about $6.1\sigma$ from one cluster, viz SPT-CL J2012-5649. The estimated photon energy flux from this cluster is approximately equal to $1.3 \times 10^{-6}$ MeV cm$^{-2}$ s$^{-1}$. The gamma-ray signal is observed between $1-10$ GeV with the best-fit spectral index equal to $-3.61 \pm 0.33$. However, since there are six radio galaxies spatially coincident with SPT-CL J2012-5649 within the Fermi-LAT PSF, we cannot rule out the possibility this signal could be caused by some of these radio galaxies. Six other SPT-SZ clusters show evidence for gamma-ray emission with significance between $3-5\sigma$. None of the remaining clusters show statistically significant evidence for gamma-ray emission.

Lalitha Sairam, Amaury H.M.J. Triaud, Thomas A. Baycroft, Jerome Orosz, Isabelle Boisse, Neda Heidari, Daniel Sebastian, Georgina Dransfield, David V. Martin, Alexandre Santerne, Matthew R. Standing

17 pages, 22 figures, 2 tables, Accepted for publication in MNRAS

Ongoing ground-based radial-velocity observations seeking to detect circumbinary planets focus on single-lined binaries even though over nine in every ten binary systems in the solar-neighbourhood are double-lined. Double-lined binaries are on average brighter, and should in principle yield more precise radial-velocities. However, as the two stars orbit one another, they produce a time-varying blending of their weak spectral lines. This makes an accurate measure of radial velocities difficult, producing a typical scatter of 10-15m/s. This extra noise prevents the detection of most orbiting circumbinary planets. We develop two new data-driven approaches to disentangle the two stellar components of a double-lined binary, and extract accurate and precise radial-velocities. Both approaches use a Gaussian Process regression, with the first one working in the spectral domain, whereas the second works on cross-correlated spectra. We apply our new methods to TIC 172900988, a proposed circumbinary system with a double-lined binary, and detect a circumbinary planet with an orbital period of 150 days, different than previously proposed. We also measure a significant residual scatter, which we speculate is caused by stellar activity. We show that our two data-driven methods outperform the traditionally used TODCOR and TODMOR, for that particular binary system.

Oscar Fernando Romero Matamala (for the JEM-EUSO Collaboration)

submitted to ICRC 2023 Proceedings

The detection of astrophysical neutrinos by IceCube in the TeV-PeV energy range motivates the development of instruments for observing these particles at higher energies. Moreover, the detection of very-high-energy (VHE) neutrinos could potentially bring constraints on ultra-high energy cosmic rays (UHECRs) source models. Tau neutrinos skimming the Earth under a shallow angle can be detected through the decay of a tau resulting in an extensive air shower (EAS) in the atmosphere. The EAS can be detected by capturing some of the optical Cherenkov signal originating from the EAS particles. To assess the feasibility of the Earth-skimming technique from high altitudes, we developed a Cherenkov telescope which was deployed on the Extreme Universe Space Observatory Super Pressure Balloon 2 (EUSO-SPB2) from Wanaka, NZ, on May $13^{th}$. It is a precursor for the Probe of Extreme Multi-Messenger Astrophysics. The 1m diameter Cherenkov telescope for EUSO-SPB2 had a focal plane comprised of 512 silicon photomultipliers (SiPMs) covering a 6.4 x 12.8 square degree field of view coupled to a 100 MS/s readout based on the GET switch capacitor ring sampler. We discuss the calibration and commissioning of the telescope and its in-flight performance.

Stefan Adelbert, Anna B. T. Penzlin, Christoph M. Schäfer, Wilhelm Kley, Billy Quarles, Rafael Sfair

9 pages, 10 figures, accepted for publication by A&A

In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the parameters to study general n-body stability of planet pairs in eccentric horseshoe orbits around binaries. The stability is tested by evolving the planet orbits for $10^5$ binary periods with varying initial semi-major axes and planet eccentricities. The unstable region of a single circumbinary planet is used as a comparison to the investigated coorbital configurations in this work. We confirm previous findings on the stability of single planets and find a first order linear relation between orbit eccentricity and pericentre to identify stable orbits for various binary configurations. Such a linear relation is also found for the stability of 1:1 resonant planets around binaries. Stable orbits for eccentric horseshoe configurations exist with a pericentre closer than seven binary separations and, in the case of Kepler-413, the pericentre of the first stable orbit can be approximated by $r_{c,peri} = (2.88 e_p + 2.46) a_{bin}$.

J. M. Casas, L. Bonavera, J. González-Nuevo, G. Puglisi, C. Baccigalupi, M. M. Cueli, D. Crespo, C. González-Gutiérrez, F. J. de Cos

11 pages, 10 figures, Submitted to A&A

Recovering the polarized cosmic microwave background (CMB) is crucial for shading light on Cosmic Inflation. Methods with different characteristics should be developed and optimized. We aim to use a neural network called CENN and train it for recovering the E and B modes of the CMB. We train the network with realistic simulations of 256x256 pixel squared patches at 100, 143 and 217 GHz Planck channels, which contain the CMB, thermal dust, synchrotron, PS and noise. We make several trainings sets: 30, 25 and 20 arcmin resolution patches at the same position in the sky. After being trained, CENN is able to recover the CMB signal at 143 GHz in Q and U patches. Then, we use NaMaster for estimating the EE and BB power spectrum for each input and output patches in the test dataset, as well as the difference between input and output power spectra and the residuals. We also test the methodology using a different foreground model at 5 arcmin resolution without noise. We recover the E-mode generally founding residuals bellow the input signal at all scales. In particular, we found a value of about 0.1 muK2 at l<200, decreasing below 0.01 muK2 at smaller scales. For the B-mode, we similarly recover the CMB with residuals between 0.01 and 0.001 muK2. We also train the network with 5 arcmin Planck simulations without noise, obtaining clearly better results with respect the previous cases. For a different foreground model, the recovery is similar, although B-mode residuals increase above the input signal. In general, we found that, the network performs better when training with the same resolution used for testing. Based on the results, CENN seems to be a promising for recovering both E and B modes at sub-degree scales in ground-base experiments such as POLARBEAR, SO and CMB-S4. Once extending its applicability at all sky, it could be an alternative component separation method for LiteBIRD satellite.

G. Pirola, J. Juryšek, M. Balbo, D. Eckert, A. Tramacere, R. Walter

Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 ( arXiv:submit/5166658 ). The paper "Multiwavelength study of the galactic PeVatron candidate LHAASO J2108+5157" is published in A&A ( this https URL )

LHAASO J2108+5157 is a recently discovered source, detected in the Ultra-High-Energy band by the LHAASO collaboration. Two molecular clouds were identified in the direction coincident with LHAASO J2108+5157 and, from the spectra reported by LHAASO, there is no sign of an energy cutoff up to 200 TeV. This source makes a promising galactic PeVatron candidate. In 2021, the Large-Sized Telescope prototype (LST-1) of the Cherenkov Telescope Array (CTA) Observatory, collected about 50 hours of quality-selected data on LHAASO J2108+5157. Through these observations, we managed to compute stringent upper limits on the source emission in the multi-TeV band. Together with the analysis of XMM-Newton data and 12 years of Fermi-LAT data, we performed a multi-wavelength study of the source, investigating different possible scenarios of particle acceleration. In this contribution, we will present the results of the analysis, as well as the multi-wavelength modeling, and consequent interpretation of different possible scenarios of emission.

C.L. Hale, D.J. Schwarz, P.N. Best, S.J. Nakoneczny, D. Alonso, D. Bacon, L. Böhme, N. Bhardwaj, M. Bilicki, S. Camera, C.S. Heneka, M. Pashapour-Ahmadabadi, P. Tiwari, J. Zheng, K.J. Duncan, M.J. Jarvis, R. Kondapally, M. Magliocchetti, H.J.A. Rottgering, T.W. Shimwell

Accepted for publication in MNRAS. 29 pages, 24 figures

Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $\mu$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies within the survey. We discuss systematics in the data and an improved methodology for generating random catalogues, compared to that used for LoTSS-DR1, before presenting the angular clustering for $\sim$900,000 sources $\geq$$1.5$ mJy and a peak signal-to-noise $\geq$$7.5$ across $\sim$$80\%$ of the observed area. Using the clustering we infer the bias assuming two evolutionary models. When fitting {angular scales of $0.5 \leq\theta<5\,\deg$, using a linear bias model, we find LoTSS-DR2 sources are biased tracers of the underlying matter, with a bias of $b_{C}= 2.14^{+0.22}_{-0.20}$ (assuming constant bias) and $b_{E}(z=0)= 1.79^{+0.15}_{-0.14}$ (for an evolving model, inversely proportional to the growth factor), corresponding to $b_E= 2.81^{+0.24}_{-0.22}$ at the median redshift of our sample, assuming the LoTSS Deep Fields redshift distribution is representative of our data. This reduces to $b_{C}= 2.02^{+0.17}_{-0.16}$ and $b_{E}(z=0)= 1.67^{+0.12}_{-0.12}$ when allowing preferential redshift distributions from the Deep Fields to model our data. Whilst the clustering amplitude is slightly lower than LoTSS-DR1 ($\geq$2 mJy), our study benefits from larger samples and improved redshift estimates.

S. J. Nakoneczny, D. Alonso, M. Bilicki, D. J. Schwarz, C. L. Hale, A. Pollo, C. Heneka, P. Tiwari, J. Zheng, M. Brüggen, M. J. Jarvis, T. W. Shimwell

We combine the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) catalogue with gravitational lensing maps from the Cosmic Microwave Background (CMB) to place constraints on the bias evolution of LoTSS radio galaxies, and on the amplitude of matter perturbations. We construct a flux-limited catalogue, and analyse its harmonic-space cross-correlation with CMB lensing maps from Planck, $C_\ell^{g\kappa}$, as well as its auto-correlation, $C_\ell^{gg}$. We explore the models describing the redshift evolution of the large-scale radio galaxy bias, discriminating between them through the combination of both $C_\ell^{g\kappa}$ and $C_\ell^{gg}$. Fixing the bias evolution, we then use these data to place constraints on the amplitude of large scale density fluctuations. We report the significance of the $C_\ell^{g\kappa}$ signal at a level of $26.6\sigma$. We determine that a linear bias evolution of the form $b_g(z) = b_{g,D} / D(z)$, where $D(z)$ is the growth rate, is able to provide a good description of the data, and measure $b_{g,D} = 1.41 \pm 0.06$ for a sample flux-limited at $1.5\,{\rm mJy}$, for scales $\ell < 250$ for $C_\ell^{gg}$, and $\ell < 500$ for $C_\ell^{g\kappa}$. At the sample's median redshift, we obtain $b(z = 0.82) = 2.34 \pm 0.10$. Using $\sigma_8$ as a free parameter, while keeping other cosmological parameters fixed to the Planck values, we find fluctuations of $\sigma_8 = 0.75^{+0.05}_{-0.04}$. The result is in agreement with weak lensing surveys, and at $1\sigma$ difference with Planck CMB constraints. We also attempt to detect the late-time integrated Sachs-Wolfe effect with LOFAR, but with the current sky coverage, the cross-correlation with CMB temperature maps is consistent with zero. Our results are an important step towards constraining cosmology with radio continuum surveys from LOFAR and other future large radio surveys.

M. Rouble, G. Smecher, A. Anderson, P. S. Barry, K. Dibert, M. Dobbs, K. S. Karkare, J. Montgomery

14 pages, 7 figures, appearing in the proceedings of SPIE Astronomical Telescopes and Instrumentation 2022

We present RF-ICE, a novel readout platform for microwave kinetic inductance detectors (MKIDs), optimized for use on millimeter-wavelength telescopes. The RF-ICE system extends ICE, a versatile, mature signal processing platform currently in use on telescopes around the world, into a new operational domain with MKIDs biased with gigahertz carriers. The system couples the FPGA-based ICE motherboard with a radio-frequency digitization daughterboard to enable direct digital synthesis from 0 to 6 GHz without the need for external mixing. The system operates two independent readout modules, each with 1024 frequency-multiplexed readout channels spaced across 500 MHz of carrier bandwidth. The system, which is under active development, is in operation with prototype detector wafers and will be deployed for the upcoming SPT-SLIM and SPT-3G+ experiments.

Anshul Verma, Sanjeet K. Patel, Pavan K. Aluri, Sukanta Panda, David F. Mota

26 pages, 5 figures, 3 tables. Comments are welcome

We study the Bianchi-I cosmological model motivated by signals of statistical isotropy violation seen in cosmic microwave background (CMB) observations and others. To that end, we consider various kinds of anisotropic matter that source anisotropy in our model, specifically cosmic strings, magnetic fields, domain walls and Lorentz violation generated magnetic fields. These anisotropic matter sources, taking one at a time, are studied for their co-evolution with standard model (isotropic) sources viz., dust-like (dark/normal) matter and dark energy modelled as cosmological constant. We constrain the Hubble parameter, density fractions of anisotropic matter, isotropic cold dark matter (CDM) and dark energy ($\Lambda$) in a Bianchi-I universe with planar symmetry i.e., with a global ellipsoidal geometry, and try to find signatures of a cosmic preferred axis if any. Pantheon Type Ia Supernova (SNIa) data is used in our study to obtain constraints on cosmological parameters of this model. In our analysis, we found a cosmic axis of anisotropy using SN1e data. Interestingly this preferred axis lies in the same direction, broadly, as some other cosmic anisotropy axes reported in the literature such as the CMB hemispherical power asymmetry, etc., with some of anisotropic sources considered. We also found a non-zero evidence for cosmic shear and eccentricity for our universe with those sources. To be more conclusive, we require more SNIa host galaxy data for tighter constraints on distance and absolute magnitude calibration which we hope will be available from the future JWST observations and other surveys.

Anowar J. Shajib

Invited review for the Proceedings of IAU Symposium 381: "Strong gravitational lensing in the era of Big Data", H. Stacey, C. Grillo, and A. Sonnenfeld eds. 7 pages, 4 figures

Galaxy-scale strong lensing is a powerful tool in Astrophysics and Cosmology, enabling studies of massive galaxies' internal structure, their formation and evolution, stellar initial mass function, and cosmological parameters. In this conference proceeding, we highlight key findings from the past decade in astrophysical applications of strong lensing at the galaxy scale. We then briefly summarize the present status of discovery and analyses of new samples from recent or ongoing surveys. Finally, we offer insights into anticipated developments in the upcoming era of big data shaping the future of this field, thanks to the Rubin, Euclid, and Roman observatories.

L. J. Janssen, P. Woitke, O. Herbort, M. Min, K. L. Chubb, Ch. helling, L. Carone

The first JWST observations of hot Jupiters showed an unexpected detection of SO2 in their hydrogen-rich atmospheres. We investigate how much sulphur can be expected in the atmospheres of rocky exoplanets and which sulphur molecules can be expected to be most abundant and detectable by transmission spectroscopy. We run thermo-chemical equilibrium models at the crust-atmosphere interface, considering surface temperatures 500 to 5000 K, surface pressures 1 to 100 bar, and various sets of element abundances based on common rock compositions. Between 1000 K and 2000 K, we find gaseous sulphur concentrations of up to 25 percent above the rock in our models. SO2, SO, H2S and S2 are by far the most abundant sulphur molecules. SO2 shows potentially detectable features in transmission spectra at about 4 micron, between 7 and 8 micron, and beyond 15 micron. In contrast, the sometimes abundant H2S molecule is difficult to detect in these spectra, which are mostly dominated by H2O and CO2. Although the molecule PS only occurs with concentrations below 300 ppm, it can cause a strong absorption feature between 0.3 and 0.65 micron in some of our models for high surface pressures. The detection of sulphur molecules would enable a better characterisation of the planetary surface.

Christina Konstantopoulou, Annalisa De Cia, Cédric Ledoux, Jens-Kristian Krogager, Lars Mattsson, Darach Watson, Kasper E. Heintz, Céline Péroux, Pasquier Noterdaeme, Anja C. Andersen, Johan P. U. Fynbo, Iris Jermann, Tanita Ramburuth-Hurt

Re-submitted to A&A. Abstract abridged

The evolution of the cosmic dust content and the cycle between metals and dust in the interstellar medium (ISM) plays a fundamental role in galaxy evolution. The chemical enrichment of the Universe can be traced through the evolution of the dust-to-metals ratio (DTM) and the dust-to-gas ratio (DTG) with metallicity. We use a novel method to determine mass estimates of the DTM, DTG and dust composition based on our previous measurements of the depletion of metals in different environments (the Milky Way, the Magellanic Clouds, and damped Lyman-$\alpha$ absorbers, DLAs, toward quasars and towards gamma-ray bursts, GRBs), which were calculated from the relative abundances of metals in the ISM through absorption-line spectroscopy column densities observed mainly from VLT/UVES and X-shooter, and HST/STIS. We derive the dust extinction from the estimated dust depletion ($A_{V, \rm depl}$) and compare with the $A_{V}$ from extinction. We find that the DTM and DTG ratios increase with metallicity and with the dust tracer [Zn/Fe]. This suggests that grain growth in the ISM is a dominant process of dust production. The increasing trend of the DTM and DTG with metallicity is in good agreement with a dust production and evolution model. Our data suggest that the stellar dust yield is much lower than the metal yield and thus that the overall amount of dust in the warm neutral medium that is produced by stars is much lower. We find that $A_{V,\rm depl}$ is overall lower than $A_{V, \rm ext}$ for the Milky Way and a few Magellanic Clouds lines of sight, a discrepancy that is likely related to the presence of carbonaceous dust. We show that the main elements that contribute to the dust composition are, O, Fe, Si, Mg, C, S, Ni and Al for all the environments. Abundances at low dust regimes suggest the presence of pyroxene and metallic iron in dust.

Jonathan Brande, Ian J. M. Crossfield, Laura Kreidberg, Caroline V. Morley, Travis Barman, Björn Benneke, Jessie L. Christiansen, Diana Dragomir, Jonathan J. Fortney, Thomas P. Greene, Kevin K. Hardegree-Ullman, Andrew W. Howard, Heather A. Knutson, Joshua D. Lothringer, Thomas Mikal-Evans

Submitted to ApJL. 11 pages, 3 figures, 4 tables

Over the last decade, precise exoplanet transmission spectroscopy has revealed the atmospheres of dozens of exoplanets, driven largely by observatories like the \textit{Hubble} Space Telescope. One major discovery has been the ubiquity of atmospheric aerosols, often blocking access to exoplanet chemical inventories. Tentative trends have been identified, showing that the clarity of planetary atmospheres may depend on equilibrium temperature. Previous work has often grouped dissimilar planets together in order to increase the statistical power of any trends, but it remains unclear from observed transmission spectra whether these planets exhibit the same atmospheric physics and chemistry. We present a re-analysis of a smaller, more physically similar sample of 15 exo-Neptune transmission spectra across a wide range of temperatures (200-1000 K). Using condensation cloud and hydrocarbon haze models, we find that the exo-Neptune population is best described by very low cloud sedimentation efficiency ($\mathrm{f_{sed}}\sim0.01$) and high metallicity ($100\times$ Solar). There is an intrinsic astrophysical scatter of $\sim0.5$ scale height, perhaps evidence of stochasticity in these planets' formation processes. Observers should expect significant attenuation in transmission spectra of Neptune-size exoplanets, up to 6 scale heights for equilibrium temperatures between 500 and 800~K. With JWST's greater wavelength sensitivity, colder ($<500$~K) planets should be high-priority targets given their comparative rarity, clearer atmospheres, and the need to distinguish between the ``super-puffs'' and more typical gas-dominated planets.

Brandon Melcher, Arnab Pradhan, Scott Watson

39 pages and seven figures

We propose a cosmological lingering phase for the initial state prior to inflation which would help address the singularity problem of inflation. The universe begins with a constant (Hagedorn) temperature and then transitions into an inflationary universe while preserving the Null Energy Condition (NEC). In such a universe time is presumably emergent, calling the age of the universe into question. We first consider the phase space of positive spatial curvature models within General Relativity and with matter sources that respect the NEC. Depending on the duration of the post-lingering inflation these models can produce a small amount of observable spatial curvature in the Cosmic Microwave Background. We also discuss how lingering can arise with or without spatial curvature in theories of quantum gravity when considering the thermodynamic scaling of particles and its impact on the early universe. The string theory dilaton is essential to the dynamics. There are many open questions that remain.

Victor V. Golovko

accepted for publication in the European Physical Journal C (19 pages, 2 tables, 5 figures)

An evaluation method supported by robust statistical analysis was used to analyze historical measurements of 39Ar half-life. The method, which combines the most frequent value (MFV) approach with bootstrap analysis, provides a more reliable way to estimate the half-life of 39Ar. The results show that the half-life is T1/2(MFV) = 268.2 + (3.1) - (2.9) years, with an uncertainty corresponding to the 68% confidence level. This uncertainty is three times smaller than the most precise re-calculated measurements by Stoenner et al. (1965) and 2.7 times smaller than the adopted half-life value in nuclear data sheets. Recently, the specific activity of the beta decay of 39Ar in atmospheric argon was measured in various underground facilities. Applying the MFV method to these measurements gives a specific activity of SA(39Ar/Ar)(MFV) = 0.966 + (0.010) - (0.018) Bq/kg(atmAr), with an uncertainty corresponding to the 68% confidence level. This paper also discusses the method used to determine the half-life of 39Ar using the specific activity of 39Ar in atmospheric argon.

N. Meyer-Vernet, A. Lecacheux, M. Moncuquet, K. Issautier, W. S. Kurth

7 pages, 4 figures, paper accepted by journal Astronomy & Astrophysics

The quasi-thermal motion of plasma particles produces electrostatic fluctuations, whose voltage power spectrum induced on electric antennas reveals plasma properties. In weakly magnetised plasmas, the main feature of the spectrum is a line at the plasma frequency -- proportional to the square root of the electron density -- whose global shape can reveal the electron temperature, while the fine structure reveals the suprathermal electrons. Since it is based on electrostatic waves, quasi-thermal noise spectroscopy (QTN) provides in situ measurements. This method has been successfully used for more than four decades in a large variety of heliosphere environments. Very recently, it has been tentatively applied in the very local interstellar medium (VLISM) to interpret the weak line discovered on board Voyager 1 and in the context of the proposed interstellar probe mission. The present paper shows that the line is still observed in the Voyager Plasma Wave Science data, and concentrates on the main features that distinguish the plasma QTN in the VLISM from that in the heliosphere. We give several tools to interpret it in this medium and highlight the errors arising when it is interpreted without caution, as has recently been done in several publications. We show recent solar wind data, which confirm that the electric field of the QTN line in a weakly magnetised stable plasma is not aligned with the local magnetic field. We explain why the amplitude of the line does not depend on the concentration of suprathermal electrons, and why its observation with a short antenna does not require a kappa electron velocity distribution. Finally, we suggest an origin for the suprathermal electrons producing the QTN and we summarise the properties of the VLISM that could be deduced from an appropriate implementation of QTN spectroscopy on a suitably designed instrument.

Filippo Camilloni, Gianluca Grignani, Troels Harmark, Marta Orselli, Daniele Pica

26 pages, 13 figures

Binary-black-hole (BBH) mergers can take place close to a supermassive black hole (SMBH) while being in a bound orbit around the SMBH. In this paper, we study such bound triple systems and show that including the strong gravity effects of describing the SMBH with a Kerr metric can significantly modify the dynamics, as compared to a Newtonian point particle description of the SMBH. We extract the dynamics of the system, using a quadrupole approximation to the tidal forces due to the SMBH. We exhibit how the gyroscope precession is built into this dynamics, and find the secular Hamiltonian by both averaging over the inner and outer orbits, the latter being the orbit of the BBH around the SMBH. We study the long-time-scale dynamics, including the periastron precession and GW radiation-reaction of the binary system, finding that the strong gravity effects of the SMBH can enhance the von Zeipel-Lidov-Kozai mechanism, resulting in more cycles, higher maximum eccentricity, and thereby a shorter merger time, particularly when the binary is close to, or at, the innermost stable orbit of the SMBH. We end with an analysis of the peak frequency of the GW emission from the binary system, highlighting possible observable signatures in the LISA and ET frequency bands.

Ryusuke Jinno, Jun'ya Kume, Masaki Yamada

Main text: 5 pages, 2 figures, Appendices: 2 pages, 2 figures

We discuss the unique phenomenology of first-order phase transitions (FOPTs) catalyzed by primordial black holes (BHs). Because of the enhancement of the vacuum nucleation rate around a BH, vacuum nucleation can be triggered only around BHs. If the average number of BHs within one Hubble volume is smaller than unity at the time of bubble nucleation, each bubble catalyzed around them expands to the Hubble size, and the Hubble patches containing BHs transition to the true vacuum, whereas the rest of the universe remains in the false vacuum. Assuming that the vacuum energy difference is negligible in relation to the total energy until a time close to the bubble collision, most false-vacuum patches do not undergo inflation, and the entire universe is eventually filled with true-vacuum bubbles much after BH formation and bubble nucleation. This "super-slow phase transition" scenario predicts enhanced gravitational wave signals from bubble collisions and can be tested in future observations. Moreover, a small fraction of the universe remains in the false vacuum and inflates by chance, and the corresponding region can be seen as a BH from an outside observer. This scenario of "birth of baby BHs from parent BHs" via the catalyzed phase transition can account for the entire dark matter in our universe.

Aleksandr Azatov, Giulio Barni, Rudin Petrossian-Byrne, Miguel Vanvlasselaer

42 pages + appendices

We quantise from first principles field theories living on the background of a bubble wall in the planar limit with particular focus on the case of spontaneous breaking of gauge symmetry. Using these tools, we compute the average momentum transfer from transition radiation: the soft emission of radiation by an energetic particle passing across the wall, with a particular focus on the longitudinal polarisation of vectors. We find these to be comparable to transverse polarisations in symmetry-breaking transitions with mild super-cooling, and dominant in broken to broken transitions with thin wall. Our results have phenomenological applications for the expansion of bubbles during first order phase transitions. Our general framework allows for the calculation of any particle processes of interest in such translation breaking backgrounds.

A. Baha Balantekin, George M. Fuller, Anupam Ray, Anna M. Suliga

10 pages, 2 figures. Comments are welcome!

The neutrinos in the diffuse supernova neutrino background (DSNB) travel over cosmological distances and this provides them with an excellent opportunity to interact with dark relics. We show that a cosmologically-significant relic population of keV-mass sterile neutrinos with strong self-interactions could imprint their presence in the DSNB. The signatures of the self-interactions would be "dips" in the otherwise smooth DSNB spectrum. Upcoming large-scale neutrino detectors, for example Hyper-Kamiokande, have a good chance of detecting the DSNB and these dips. If no dips are detected, this method serves as an independent constraint on the sterile neutrino self-interaction strength and mixing with active neutrinos. We show that relic sterile neutrino parameters that evade X-ray and structure bounds may nevertheless be testable by future detectors like TRISTAN, but may also produce dips in the DSNB which could be detectable. Such a detection would suggest the existence of a cosmologicaly-significant, strongly self-interacting sterile neutrino background, likely embedded in a richer dark sector.

Andrew Steinmetz

PhD thesis; 155 pages, 17 figures. Includes work done in collaboration with Martin Formanek, Cheng Tao Yang, and Johann Rafelski

Magnetism is a rich subject touching all aspects of physics. My goal with this dissertation is to explore spin and magnetic moments in \emph{relativistic} mechanics from both a quantum and classical perspective. We emphasize the special case of gyromagnetic ratio $g\!=\!2$ and its relationship to the algebraic spin structure of wave equations. In relativistic quantum mechanics, we investigate generalizations of the Dirac equation for arbitrary magnetic moments for fermions. We analyze the homogeneous magnetic field case and the Coulomb problem for hydrogen-like atoms with emphasis on the role of the anomalous magnetic moment (AMM). We explore alternative approaches which combine mass and the magnetic moment. Classically, we propose a relativistic covariant model of the Stern-Gerlach force via the introduction of a magnetic four-potential. This model modifies the covariant torque equations and unites the Amp\`erian and Gilbertian models for magnetic moments. We further study (transition) magnetic dipoles in Majorana neutrinos specifically analyzing the relationship between flavor mixing and electromagnetic (EM) fields. We demonstrate EM flavor mixing explicitly in the 2-flavor model and develop a dynamical mass basis with an EM rotation matrix. EM induced neutrino mass splitting is compared to neutrino mass hierarchy. An interesting application of these theoretical developments is to study primordial magnetization in the early universe during the hot dense electron-positron plasma epoch. We propose a model of magnetic thermal matter-antimatter plasmas. We analyze the paramagnetic characteristics of electron-positron plasma when exposed to an external primordial field. Future research outlooks include: Second order equations for anomalous quantum chromodynamic (QCD) moments, neutrino CP violation in strong EM fields, and fifth-dimension spin dynamics in Kaluza-Klein theory.

Elisa Maggio

14 pages, 5 figures, contribution to the book "Modified and Quantum Gravity - From theory to experimental searches on all scales"

Gravitational waves open the possibility to investigate the nature of compact objects and probe the horizons of black holes. Some models of modified gravity predict the presence of horizonless and singularity-free compact objects. Such dark compact objects would emit a gravitational-wave signal which differs from the standard black hole scenario. In this chapter, we overview the phenomenology of dark compact objects by analysing their characteristic frequencies in the ringdown and the emission of gravitational-wave echoes in the postmerger signal. We show that future gravitational-wave detectors will allow us to perform model-independent tests of the black hole paradigm.

Alex Roskill, Marienza Caldarola, Sachiko Kuroyanagi, Savvas Nesseris

8 pages, 4 figures

In this paper, we study the next-to-leading order corrections of gravitational wave production by hyperbolic encounters of compact objects. The signal is a burst event, with the majority of the released energy occurring during the closest approach. We discuss the lowest order term contribution to the emission, which is the quadrupole radiation, and also the next-to-leading terms, i.e., the mass octupole and the quadrupole current radiation, investigating the relative contribution to the power emitted in gravitational waves and total energy emitted, both in the time and frequency domains. We find specific configurations of systems where these corrections could be important and should be taken into account when analysing burst events.

The recent astronomical milestone by the pulsar timing arrays (PTA) has revealed galactic-size gravitational waves (GW) in the form of a stochastic gravitational wave background (SGWB), correlating the radio pulses emitted by millisecond pulsars. This draws the outstanding questions toward the origin and the nature of the SGWB; the latter is synonymous to testing how quadrupolar the inter-pulsar spatial correlation is. In this paper, we tackle the nature of the SGWB by considering correlations beyond the Hellings-Downs (HD) curve of Einstein's general relativity. We put the HD and non-Einsteinian GW correlations under scrutiny with the NANOGrav and the CPTA data, and find that both data sets allow a graviton mass $m_{\rm g} \lesssim 1.04 \times 10^{-22} \ {\rm eV}/c^2$ and subluminal traveling waves. We discuss gravitational physics scenarios beyond general relativity that could host non-Einsteinian GW correlations in the SGWB and highlight the importance of the cosmic variance inherited from the stochasticity in interpreting PTA observation.

Jun-Chen Wang, Zi-Yan Yuwen, Yu-Shi Hao, Shao-Jiang Wang

23 pages, 1 figure

The gravitational-wave energy-density spectra from cosmological first-order phase transitions crucially depend on the terminal wall velocity of asymptotic bubble expansion when the driving force from the effective potential difference is gradually balanced by the backreaction force from the thermal plasma. Much attention has previously focused on the backreaction force acting on the bubble wall alone but overlooked the backreaction forces on the sound shell and shockwave front, if any, which have been both numerically and analytically accomplished in our previous studies but only for a bag equation of state. In this paper, we will generalize the backreaction force on bubble expansion beyond the simple bag model.