Abstract visibility
Abstract text size

Papers for Thursday, Nov 16 2023

Papers with local authors

Keiya Hirashima, Kana Moriwaki, Michiko S. Fujii, Yutaka Hirai, Takayuki R. Saitoh, Junichiro Makino, Shirley Ho

11 pages, 9 figures, Accepted for the NeurIPS 2023 AI4Science Workshop

0 votes
Paper 10 — arXiv:2311.08460
0 votes
Paper 10 — arXiv:2311.08460

Some stars are known to explode at the end of their lives, called supernovae (SNe). The substantial amount of matter and energy that SNe release provides significant feedback to star formation and gas dynamics in a galaxy. SNe release a substantial amount of matter and energy to the interstellar medium, resulting in significant feedback to star formation and gas dynamics in a galaxy. While such feedback has a crucial role in galaxy formation and evolution, in simulations of galaxy formation, it has only been implemented using simple {\it sub-grid models} instead of numerically solving the evolution of gas elements around SNe in detail due to a lack of resolution. We develop a method combining machine learning and Gibbs sampling to predict how a supernova (SN) affects the surrounding gas. The fidelity of our model in the thermal energy and momentum distribution outperforms the low-resolution SN simulations. Our method can replace the SN sub-grid models and help properly simulate un-resolved SN feedback in galaxy formation simulations. We find that employing our new approach reduces the necessary computational cost to $\sim$ 1 percent compared to directly resolving SN feedback.

Dong Li, Zhenyong Hou, Xianyong Bai, Chuan Li, Matthew Fang, Haisheng Zhao, Jincheng Wang, Zongjun Ning

accepted for publication in the Science China Technological Sciences

0 votes
Paper 36 — arXiv:2311.08767
0 votes
Paper 36 — arXiv:2311.08767

Kink oscillations, which are frequently observed in coronal loops and prominences, are often accompanied by extreme-ultraviolet (EUV) waves. However, much more needs to be explored regarding the causal relationships between kink oscillations and EUV waves. In this article, we report the simultaneous detection of kink oscillations and EUV waves that are both associated with an X2.1 flare on 2023 March 03 (SOL2023-03-03T17:39). The kink oscillations, which are almost perpendicular to the axes of loop-like structures, are observed in three coronal loops and one prominence. One short loop shows in-phase oscillation within the same period of 5.2 minutes at three positions. This oscillation could be triggered by the pushing of an expanding loop and interpreted as the standing kink wave. Time lags are found between the kink oscillations of the short loop and two long loops, suggesting that the kink wave travels in different loops. The kink oscillations of one long loop and the prominence are possibly driven by the disturbance of the CME, and that of another long loop might be attributed to the interaction of the EUV wave. The onset time of the kink oscillation of the short loop is nearly same as the beginning of an EUV wave. This fact demonstrates that they are almost simultaneous. The EUV wave is most likely excited by the expanding loop structure and shows two components. The leading component is a fast coronal wave, and the trailing one could be due to the stretching magnetic field lines.

Elena Gallo, Edmund Hodges-Kluck, Tommaso Treu, Vivienne Baldassare, Anil Seth, Jenny Greene, Fabio Pacucci, Richard Plotkin, Amy Reines, Belinda Wilkes

This White Paper is part of a series commissioned for the AXIS Probe Concept Mission

0 votes
Paper 58 — arXiv:2311.09161
0 votes
Paper 58 — arXiv:2311.09161

The fraction of local dwarf galaxies that hosts massive black holes is arguably the cleanest diagnostic of the dominant seed formation mechanism of today's supermassive black holes. A 5 per cent constraint on this quantity can be achieved with AXIS observations of 3300 galaxies across the mass spectrum through a combination of serendipitous extra-galactic fields plus a dedicated 1 Msec GO program.

Chen-Yu Chuang, Christian Kragh Jespersen, Yen-Ting Lin, Shirley Ho, Shy Genel

18 pages, 7 figures, 3 tables, 4 pages of appendices. Submitted to ApJ

0 votes
Paper 59 — arXiv:2311.09162
0 votes
Paper 59 — arXiv:2311.09162

Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only simulations and galaxies are strongly preferred, as producing mocks from full-fledged hydrodynamical simulations or semi-analytical models is too expensive. Here we present a Graph Neural Network-based model that is able to accurately predict key properties of galaxies such as stellar mass, $g-r$ color, star formation rate, gas mass, stellar metallicity, and gas metallicity, purely from dark matter properties extracted from halos along the full assembly history of the galaxies. Tests based on the TNG300 simulation of the IllustrisTNG project show that our model can recover the baryonic properties of galaxies to high accuracy, over a wide redshift range ($z = 0-5$), for all galaxies with stellar masses more massive than $10^9\,M_\odot$ and their progenitors, with strong improvements over the state-of-the-art methods. We further show that our method makes substantial strides toward providing an understanding of the implications of the IllustrisTNG galaxy formation model.

All other papers

Joshua Goodeve

7 pages, 5 figures. Submitted to Astronomy & Astrophysics, 13 November 2023

Context. The accuracy of photometric calibration has gradually become a limiting factor in various fields of astronomy, limiting the scientific output of a host of research. Calibration using artificial light sources in low earth orbit remains largely unexplored. Aims. We aim to demonstrate that photometric calibration using light sources in low earth orbit is a viable and competitive alternative/complement to current calibration techniques, and explore the associated ideas and basic theory. Methods. We present the publicly-available Python code Streaktools as a means to simulate and perform photometric calibration using real and simulated light streaks. We use Streaktools to perform `pill' aperture photometry on 131 simulated streaks, and MCMC-based PSF model-fitting photometry on 425 simulated streaks in an attempt to recover the magnitude zeropoint of a real exposure of the DECam instrument on the Blanco 4m telescope. Results. We show that calibration using pill photometry is too inaccurate to be useful, but that PSF photometry is able to produce unbiased and accurate ($1\sigma$ error = 3.4mmag) estimates of the zeropoint of a real image in a realistic scenario, with a reasonable light source.

Franck Selsis, Jérémy Leconte, Martin Turbet, Guillaume Chaverot, Émeline Bolmont

23 pages, 3 figures in the main text, 6 figures and 1 table in the Methods. Received 14 July 2022; Accepted 24 May 2023; Published 09 August 2023

Water vapour atmospheres with content equivalent to the Earth's oceans, resulting from impacts or a high insolation, were found to yield a surface magma ocean. This was, however, a consequence of assuming a fully convective structure. Here we report, using a consistent climate model, that pure steam atmospheres are commonly shaped by radiative layers, making their thermal structure strongly dependent on the stellar spectrum and internal heat flow. Yhe surface cooler when an adiabatic profile is not imposed: melting Earth's crust requires an insolation several times higher than today, which will not happen during the main-sequence of the Sun. Venus' surface can solidify before the steam atmosphere escapes, which is opposite to previous works. Around the reddest stars ($T_{eff}<$3000K), surface magma oceans cannot form by stellar forcing alone, whatever the water content. These findings affect observable signatures of steam atmospheres and exoplanet mass-radius relationships, drastically changing current constraints on the water content of Trappist-1 planets. Unlike adiabatic structures, radiative-convective profiles are sensitive to opacities. New measurements of poorly constrained high-pressure opacities, in particular far from the H$_2$O absorption bands, are thus necessary to refine models of steam atmospheres, which are important stages in terrestrial planet evolution.

Fabio Pacucci, Bryan Seepaul, Yueying Ni, Nico Cappelluti, Adi Foord

This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website: this http URL with a mission overview here: arXiv:2311.00780 . Review article, 7 pages, 3 figures

This white paper explores the detectability of intermediate-mass black holes (IMBHs) wandering in the Milky Way (MW) and massive local galaxies, with a particular emphasis on the role of AXIS. IMBHs, ranging within $10^{3-6} \, M_\odot$, are commonly found at the centers of dwarf galaxies and may exist, yet undiscovered, in the MW. By using model spectra for advection-dominated accretion flows (ADAFs), we calculated the expected fluxes emitted by a population of wandering IMBHs with a mass of $10^5 \, M_\odot$ in various MW environments and extrapolated our results to massive local galaxies. Around $40\%$ of the potential population of wandering IMBHs in the MW can be detected in an AXIS deep field. We proposed criteria to aid in selecting IMBH candidates using already available optical surveys. We also showed that IMBHs wandering in $>200$ galaxies within $10$ Mpc can be easily detected with AXIS when passing within dense galactic environments (e.g., molecular clouds and cold neutral medium). In summary, we highlighted the potential X-ray detectability of wandering IMBHs in local galaxies and provided insights for guiding future surveys. Detecting wandering IMBHs is crucial for understanding their demographics, evolution, and the merging history of galaxies.

Alessandro Peca, Nico Cappelluti, Stefano Marchesi, Edmund Hodges-Kluck, Adi Foord

This study presents the capabilities of the AXIS telescope in estimating redshifts from X-ray spectra alone (X-ray redshifts, XZs). Through extensive simulations, we establish that AXIS observations enable reliable XZ estimates for more than 5500 obscured Active Galactic Nuclei (AGN) up to redshift $z\sim 6$ in the proposed deep (7 Ms) and intermediate (375 ks) surveys. Notably, at least 1600 of them are expected to be in the Compton-Thick regime ($\log N_H/\mathrm{cm^{-2}}\geq 24$), underscoring the pivotal role of AXIS in sample these elusive objects that continue to be poorly understood. XZs provide an efficient alternative for optical/infrared faint sources, overcoming the need for time-consuming spectroscopy, potential limitations of photometric redshifts, and potential issues related to multi-band counterpart association. This approach will significantly enhance the accuracy of constraints on the X-ray luminosity function and obscured AGN fractions up to high redshift. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website (this http URL) with a mission overview here: arXiv:2311.00780.

Thomas Connor, Eduardo Bañados, Nico Cappelluti, Adi Foord

7 Pages, 2 Figures. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website ( this http URL ) with a mission overview here: arXiv:2311.00780

Jets powered by AGN in the early Universe ($z \gtrsim 6$) have the potential to not only define the trajectories of the first-forming massive galaxies but to enable the accelerated growth of their associated SMBHs. Under typical assumptions, jets could even rectify observed quasars with light seed formation scenarios; however, not only are constraints on the parameters of the first jets lacking, observations of these objects are scarce. Owing to the significant energy density of the CMB at these epochs capable of quenching radio emission, observations will require powerful, high angular resolution X-ray imaging to map and characterize these jets. As such, \textit{AXIS} will be necessary to understand early SMBH growth and feedback.

Ian W. Stephens, Zhe-Yu Daniel Lin, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Haifeng Yang, Rachel Harrison, Akimasa Kataoka, Carlos Carrasco-Gonzalez, Satoshi Okuzumi, Ryo Tazaki

Published in Nature

Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $\mu$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $\mu$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations.

Conny Aerts, Andrew Tkachenko

Invited review paper (31 pages, 11 figures). This is the submitted version. We welcome comments by the community; those that reach us by 15 December 2023 will be taken into account

Rotation matters for the life of a star. It causes a multitude of dynamical phenomena in the stellar interior during a star's evolution and its effects accumulate until the star dies. All stars rotate at some level but those born with a mass above about 1.3 times the mass of the Sun rotate rapidly during more than 90% of their nuclear lifetime. Internal rotation guides the angular momentum and chemical element transport throughout the stellar interior. These transport processes change over time as the star evolves. The cumulative effects of stellar rotation and its induced transport processes determine the helium content of the core by the time it exhausts its hydrogen isotopes. The amount of helium at that stage also guides the heavy element yields at the end of the star's life. A proper theory of stellar evolution and any realistic models for the chemical enrichment of galaxies, must be based on observational calibrations of stellar rotation and of the induced transport processes. Since a few years, asteroseismology offers such calibrations, for single and binary stars. We review the current status of asteroseismic modelling of rotating stars for different stellar mass regimes, in an accessible way for the non-expert. While doing so, we describe exciting opportunities sparked by asteroseismology for various domains in astrophysics, touching upon topics from exoplanetary science to galactic structure and evolution towards gravitational wave physics. Along the way, we provide ample sneak-previews for future 'industrialised' applications of asteroseismology to slow and rapid rotators, from exploitation of combined Kepler, TESS, PLATO, Gaia, and spectroscopy surveys. We end the review with a list of take away messages and achievements of asteroseismology, which are of relevance for many fields of astrophysics.

Ka-Wah Wong, Helen Russell, Jimmy Irwin, Nico Cappelluti, Adi Foord

7 pages, 3 figures. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website this http URL with a mission overview here arXiv:2311.00780

Hot gas around a supermassive black hole (SMBH) should be captured within the gravitational "sphere of influence", characterized by the Bondi radius. Deep Chandra observations have spatially resolved the Bondi radii of at least five nearby SMBHs. Contrary to earlier hot accretion models that predicted a steep temperature increase within the Bondi radius, none of the resolved temperature profiles exhibit such an increase. The temperature inside the Bondi radius appears to be complex, indicative of a multi-temperature phase of hot gas with a cooler component at about 0.2-0.3 keV. The density profiles within the Bondi regions are shallow, suggesting the presence of strong outflows. These findings might be explained by recent realistic numerical simulations that suggest that large-scale accretion inside the Bondi radius can be chaotic, with cooler gas raining down in some directions and hotter gas outflowing in others. With an angular resolution similar to Chandra and a significantly larger collecting area, AXIS will collect enough photons to map the emerging accretion flow within and around the "sphere of influence" of a sample of active galactic nuclei (AGNs). AXIS will reveal transitions in the inflow that ultimately fuels the AGN, as well as outflows that provide feedback to the environment.

Uri Keshet, Ilya Gurwich, Assaf Lavi, Dina Avitan, Teodor Linnik

6 pages, 3 figures, SI; comments welcome

The full, radio to $\gamma$-ray spectrum of the Fermi bubbles is shown to be consistent with standard strong-shock electron acceleration at the bubble edge, without ad-hoc energy cutoffs, if the ambient interstellar radiation is strong; the $\gamma$-ray cooling break should then have a microwave counterpart, undetected until now. Indeed, a broadband bubble-edge analysis uncovers a pronounced downstream dust component, which masked the anticipated $\sim35$ GHz spectral break, and the same overall radio softening consistent with Kraichnan diffusion previously reported in $\gamma$-rays.

Dongwook Lim, Andreas J. Koch-Hansen, Seungsoo Hong, Sang-Hyun Chun, Young-Wook Lee

16 pages, 9 figures, accepted for publication in AJ

Wide binaries, with separations between two stars from a few AU to more than several thousand AU, are valuable objects for various research topics in Galactic astronomy. As the number of newly reported wide binaries continues to increase, studying the chemical abundances of their component stars becomes more important. We conducted high-resolution near-infrared (NIR) spectroscopy for six pairs of wide binary candidates using the Immersion Grating Infrared Spectrometer (IGRINS) at the Gemini-South telescope. One pair was excluded from the wide binary samples due to a significant difference in radial velocity between its component stars, while the remaining five pairs exhibited homogeneous properties in 3D motion and chemical composition among the pair stars. The differences in [Fe/H] ranged from 0.00 to 0.07 dex for these wide binary pairs. The abundance differences between components are comparable to the previous results from optical spectroscopy for other samples. In addition, when combining our data with literature data, it appears that the variation of abundance differences increases in wide binaries with larger separations. However, the SVO2324 and SVO3206 showed minimal differences in most elements despite their large separation, supporting the concept of multiple formation mechanisms depending on each wide binary. This study is the first approach to the chemical properties of wide binaries based on NIR spectroscopy. Our results further highlight that NIR spectroscopy is an effective tool for stellar chemical studies based on equivalent measurements of chemical abundances from the two stars in each wide binary system.

Zhen Yuan, Chengdong Li, Nicolas F. Martin, Giacomo Monari, Benoit Famaey, Arnaud Siebert, Anke Ardern-Arentsen, Federico Sestito, Guillaume F. Thomas, Vanessa Hill, Rodrigo A. Ibata, Georges Kordopatis, Else Starkenburg, Akshara Viswanathan

submitted to A&A, comments are welcome

The most metal-poor stars (e.g. [Fe/H] $\leq-2.5$) are the ancient fossils from the early assembly epoch of our Galaxy, very likely before the formation of the thick disc. Recent studies have shown that a non-negligible fraction of them have prograde planar orbits, which makes their origin a puzzle. It has been suggested that a later-formed rotating bar could have driven these old stars from the inner Galaxy outward, and transformed their orbits to be more rotation-dominated. However, it is not clear if this mechanism can explain these stars as observed in the solar neighborhood. In this paper, we explore the possibility of this scenario by tracing these stars backwards in an axisymmetric Milky Way potential with a bar perturber. We integrate their orbits backward for 6 Gyr under two bar models: one with a constant pattern speed and another one with a decelerating speed. Our experiments show that, under the constantly-rotating bar model, the stars of interest are little affected by the bar and cannot have been shepherded from a spheroidal inner Milky Way to their current orbits. In the extreme case of a rapidly decelerating bar, some of the very metal-poor stars on planar and prograde orbits can be brought from the inner Milky Way, but $\sim90\%$ of them were nevertheless already rotation-dominated ($J_{\phi}$ $\geq$ 1000 km s$^{-1}$ kpc) 6 Gyr ago. The chance of these stars having started with spheroid-like orbits with small rotation ($J_{\phi}$ $\lesssim$ 600 km s$^{-1}$ kpc) is very low ($<$ 3$\%$). We therefore conclude that, within the solar neighborhood, the bar is unlikely to have shepherded a significant fraction of inner Galaxy spheroid stars to produce the overdensity of stars on prograde, planar orbits that is observed today.

Fabio Vito, Paolo Tozzi, Roberto Gilli, Stefano Marchesi, Nico Cappelluti, Adi Foord

8 pages, 3 figures. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website ( this http URL ) with a mission overview at arXiv:2311.00780

Overdense regions at high redshift ($z \gtrsim 2$) are perfect laboratories to study the relations between environment and SMBH growth, and the AGN feedback processes on the surrounding galaxies and diffuse gas. In this white paper, we discuss how AXIS will 1) constrain the AGN incidence in protoclusters, as a function of parameters such as redshift, overdensity, mass of the structure; 2) search for low-luminosity and obscured AGN in the satellite galaxies of luminous QSOs at $z>6$, exploiting the large galaxy density around such biased objects; 3) probe the AGN feedback on the proto-ICM via the measurement of the AGN contribution to the gas ionization and excitation, and the detection of extended X-ray emission from the ionized gas and from radio jets; 4) discover new large-scale structures in the wide and deep AXIS surveys as spikes in the redshift distribution of X-ray sources. These goals can be achieved only with an X-ray mission with the capabilities of AXIS, ensuring a strong synergy with current and future state-of-the-art facilities in other wavelengths. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at this http URL with a mission overview at https://arxiv.org/abs/2311.00780.

Vikram Khaire, Teng Hu, Joseph F. Hennawi, Joseph N. Burchett, Michael Walther, Frederick Davies

21 pages (including 4 page appendix), Submitted to MNRAS

We explore the potential of using the low-redshift Lyman-$\alpha$ (Ly$\alpha$) forest surrounding luminous red galaxies (LRGs) as a tool to constrain active galactic nuclei (AGN) feedback models. Our analysis is based on snapshots from the Illustris and IllustrisTNG simulations at a redshift of $z=0.1$. These simulations offer an ideal platform for studying the influence of AGN feedback on the gas surrounding galaxies, as they share the same initial conditions and underlying code but incorporate different feedback prescriptions. Both simulations show significant impacts of feedback on the temperature and density of the gas around massive halos. Following our previous work, we adjusted the UV background in both simulations to align with the observed number density of Ly$\alpha$ lines ($\rm dN/dz$) in the intergalactic medium and study the Ly$\alpha$ forest around massive halos hosting LRGs, at impact parameters ($r_{\perp}$) ranging from 0.1 to 100 pMpc. Our findings reveal that $\rm dN/dz$, as a function of $r_{\perp}$, is approximately 1.5 to 2 times higher in IllustrisTNG compared to Illustris up to $r_{\perp}$ of $\sim 10$ pMpc. To further assess whether existing data can effectively discern these differences, we search for archival data containing spectra of background quasars probing foreground LRGs. Through a feasibility analysis based on this data, we demonstrate that ${\rm dN/dz} (r_{\perp})$ measurements can distinguish between feedback models of IllustrisTNG and Illustris with a precision exceeding 12$\sigma$. This underscores the potential of ${\rm dN/dz} (r_{\perp})$ measurements around LRGs as a valuable benchmark observation for discriminating between different feedback models.

M. Béthermin, C. Accard, C. Guillaume, M. Dessauges-Zavadsky, E. Ibar, P. Cassata, T. Dereveaux, A. Faisst, J. Freundlich, G. C. Jones, K. Kraljic, H. Algera, R. O. Amorin, S. Bardelli, M. Boquien, V. Buat, E. Donghia, Y. Dubois, A. Ferrara, Y. Fudamoto, M. Ginolfi, P. Guillard, M. Giavalisco, C. Gruppioni, G. Gururajan, N. Hathi, C. C. Hayward, A. M. Koekemoer, B. C. Lemaux, G. E. Magdis, J. Molina, D. Narayanan, L. Mayer, F. Pozzi, F. Rizzo, M. Romano, L. Tasca, P. Theulé, D. Vergani, L. Vallini, G. Zamorani, A. Zanella, E. Zucca

9 pages, 7 figures, 2 tables, accepted by A&A (letter)

The Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($\Sigma_{\rm gas}$-$\Sigma_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resolution. We obtained ~0.3"-resolution [CII] and continuum maps of our objects, which we then converted into gas and obscured SFR surface density maps. In addition, we produced unobscured SFR surface density maps by convolving Hubble ancillary data in the rest-frame UV. We then derived the average $\Sigma_{\rm SFR}$ in various $\Sigma_{\rm gas}$ bins, and estimated the uncertainties using a Monte Carlo sampling. Our galaxy sample follows the KS relation measured in main-sequence galaxies at lower redshift and is slightly lower than predictions from simulations. Our data points probe the high end both in terms of $\Sigma_{\rm gas}$ and $\Sigma_{\rm gas}$, and gas depletion timescales (285-843 Myr) remain similar to z~2 objects. However, three of our objects are clearly morphologically disturbed, and we could have expected shorter gas depletion timescales (~100 Myr) similar to merger-driven starbursts at lower redshifts. This suggests that the mechanisms triggering starbursts at high redshift may be different than in the low- and intermediate-z Universe.

David Kipping, Jason Wright

Accepted for publication in The Astronomical Journal

The search for extraterrestrial (alien) life is one of the greatest scientific quests yet raises fundamental questions about just what we should be looking for and how. We approach alien hunting from the perspective of an experimenter engaging in binary classification with some true and confounding positive probability (TPP and CPP). We derive the Bayes factor in such a framework between two competing hypotheses, which we use to classify experiments as either impotent, imperfect or ideal. Similarly, the experimenter can be classified as dogmatic, biased or agnostic. We show how the unbounded explanatory and evasion capability of aliens poses fundamental problems to experiments directly seeking aliens. Instead, we advocate framing the experiments as looking for that outside of known processes, which means the hypotheses we test do not directly concern aliens per se. To connect back to aliens requires a second level of model selection, for which we derive the final odds ratio in a Bayesian framework. This reveals that it is fundamentally impossible to ever establish alien life at some threshold odds ratio, $\mathcal{O}_{\mathrm{crit}}$, unless we deem the prior probability that some as-yet-undiscovered natural process could explain the event is less than $(1+\mathcal{O}_{\mathrm{crit}})^{-1}$. This elucidates how alien hunters need to carefully consider the challenging problem of how probable unknown unknowns are, such as new physics or chemistry, and how it is arguably most fruitful to focus on experiments for which our domain knowledge is thought to be asymptotically complete.

I.V. Yanes-Rizo, M. A. P. Torres, J. Casares, M. Monelli, P. G. Jonker, T. Abbot, M. Armas Padilla, T. Muñoz-Darias

7 pages, 5 figures

We present VLT spectroscopy, high-resolution imaging and time-resolved photometry of KY TrA, the optical counterpart to the X-ray binary A 1524-61. We perform a refined astrometry of the field, yielding improved coordinates for KY TrA and the field star interloper of similar optical brightness that we locate $0.64 \pm 0.04$ arcsec SE. From the spectroscopy, we refine the radial velocity semi-amplitude of the donor star to $K_2 = 501 \pm 52$ km s$^{-1}$ by employing the correlation between this parameter and the full-width at half-maximum of the H$\alpha$ emission line. The $r$-band light curve shows an ellipsoidal-like modulation with a likely orbital period of $0.26 \pm 0.01$ d ($6.24 \pm 0.24$ h). These numbers imply a mass function $f(M_1) = 3.2 \pm 1.0$ M$_\odot$. The KY TrA de-reddened quiescent colour $(r-i)_0 = 0.27 \pm 0.08$ is consistent with a donor star of spectral type K2 or later, in case of significant accretion disc light contribution to the optical continuum. The colour allows us to place a very conservative upper limit on the companion star mass, $M_2 \leq 0.94$ M$_\odot$, and, in turn, on the binary mass ratio, $q = M_2/M_1 \leq 0.31$. By exploiting the correlation between the binary inclination and the depth of the H$\alpha$ line trough, we establish $i = 57 \pm 13$ deg. All these values lead to a compact object and donor mass of $M_1 = 5.8^{+3.0}_{-2.4}$ M$_\odot$ and $M_2 = 0.5 \pm 0.3$ M$_\odot$, respectively, thus confirming the black hole nature of the accreting object. In addition, we estimate a distance toward the system of $8.0 \pm 0.9$ kpc.

Catherine A. Clark, Gerard T. van Belle, Elliott P. Horch, Michael B. Lund, David R. Ciardi, Kaspar von Braun, Jennifer G. Winters, Mark E. Everett, Zachary D. Hartman, Joe Llama

24 pages, 11 figures. Accepted for publication in the Astronomical Journal

Stellar multiplicity is correlated with many stellar properties, yet multiplicity measurements have proven difficult for the M dwarfs -- the most common type of star in our galaxy -- due to their faintness and the fact that a reasonably-complete inventory of later M dwarfs did not exist until recently. We have therefore carried out the Pervasive Overview of "Kompanions" of Every M dwarf in Our Neighborhood (POKEMON) survey, which made use of the Differential Speckle Survey Instrument on the 4.3-meter Lowell Discovery Telescope, along with the NN-EXPLORE Exoplanet Stellar Speckle Imager on the 3.5-meter WIYN telescope. The POKEMON sample is volume-limited from M0V through M9V out to 15 pc, with additional brighter targets at larger distances. In total, 1124 targets were observed. New discoveries were presented in the first paper in the series. In this second paper in the series, we present all detected companions, gauge our astrometric and photometric precision, and compare our filtered and filterless speckle observations. We find that the majority (58.9%) of the companions we detect in our speckle images are not resolved in Gaia, demonstrating the need for high-resolution imaging in addition to long-term astrometric monitoring. Additionally, we find that the majority (73.2%) of simulated stellar companions would be detectable by our speckle observations. Specifically within 100 au, we find that 70.2% of simulated companions are recovered. Finally, we discuss future directions of the POKEMON survey.

Miguel Coloma Puga (1, 2) Barbara Balmaverde (2) Alessandro Capetti (2) Francesco Massaro (1) Cristina Ramos Almeida (3, 4) George Miley (5) Roberto Gilli (6) Alessandro Marconi (7, 8) ((1) Department of Physics, Università di Torino Via Pietro Giuria, 1, 10125, Torino, Italy (2) INAF - Osservatorio Astrofisico di Torino Via Osservatorio 20, I-10025 Pino Torinese, Italy (3) Instituto de Astrofisica de Canarias C. Vía Láctea, s/n, 38205 La Laguna, Tenerife, Spain (4) Departamento de Astrofísica, Universidad de La Laguna 38206, La Laguna, Tenerife, Spain (5) Leiden Observatory, Leiden University PO Box 9513, 2300 RA Leiden, The Netherlands (6) Osservatorio di Astrofisica e Scienza dello Spazio di Bologna Via Gobetti 93/3, I-40129 Bologna, Italy (7) Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze Via G. Sansone 1, I-50019, Sesto Fiorentino, Firenze, Italy (8) INAF - Osservatorio Astrofisico di Arcetri Largo E. Fermi 5, I-50125, Firenze, Italy)

Pre-proofs version - Accepted for publication in ApJL

Spatially resolved observations of AGN host galaxies undergoing feedback processes are one of the most relevant avenues through which galactic evolution can be studied, given the long lasting effects AGN feedback has on gas reservoirs, star formation, and AGN environments at all scales. Within this context we report results from VLT/MUSE integral field optical spectroscopy of TN J1049-1258, one of the most powerful radio sources known, at a redshift of 3.7. We detected extended ($\sim$ 18 kpc) Lyman $\alpha$ emission, spatially aligned with the radio axis, redshifted by 2250 $\pm$ 60 km s$^{-1}$ with respect to the host galaxy systemic velocity, and co-spatial with UV continuum emission. This Lyman $\alpha$ emission could arise from a companion galaxy, although there are arguments against this interpretation. Alternatively, it might correspond to an outflow of ionized gas stemming from the radio galaxy. The outflow would be the highest redshift spatially resolved ionized outflow to date. The enormous amount of energy injected, however, appears to be unable to quench the host galaxy's prodigious star formation, occurring at a rate of $\sim$4500 M$_{\odot} yr^{-1}$, estimated using its far infra-red luminosity. Within the field we also found two companion galaxies at projected distances of $\sim$25 kpc and $\sim$60 kpc from the host, which suggests the host galaxy is harbored within a protocluster.

Theodoros Smponias

Published in 2021. arXiv admin note: text overlap with arXiv:2311.07694

Microquasar binary stellar systems emit electromagnetic radiation and high-energy particles over a broad energy spectrum. However, they are so far away that it is hard to observe their details. A simulation offers the link between relatively scarce observational data and the rich theoretical background. In this work, high-energy particle emission from simulated twin microquasar jets is calculated in a unified manner. From the cascade of emission within an element of jet matter to the dynamic and radiative whole jet model, the series of physical processes involved are integrated together. A programme suite assembled around model data produces synthetic images and spectra directly comparable to potential observations by contemporary arrays. The model is capable of describing a multitude of system geometries, incorporating increasing levels of realism depending on need and available computational resources. As an application, the modelling process is applied to a typical microquasar, which is synthetically observed from different angles using various imaging geometries. Furthermore, the resulting intensities are comparable to the sensitivity of existing detectors. The combined background emission from a potential distribution of microquasars is also modelled.

Core Francisco Park, Victoria Ono, Nayantara Mudur, Yueying Ni, Carolina Cuesta-Lazaro

Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter components that cannot be directly observed. The relationship between dark matter density fields and galaxy distributions can be sensitive to assumptions in cosmology and astrophysical processes embedded in the galaxy formation models, that remain uncertain in many aspects. Based on state-of-the-art galaxy formation simulation suites with varied cosmological parameters and sub-grid astrophysics, we develop a diffusion generative model to predict the unbiased posterior distribution of the underlying dark matter fields from the given stellar mass fields, while being able to marginalize over the uncertainties in cosmology and galaxy formation.

Ellen L. Sirks, Richard Massey, Ajay S. Gill, Jason Anderson, Steven J. Benton, Anthony M. Brown, Paul Clark, Joshua English, Spencer W. Everett, Aurelien A. Fraisse, Hugo Franco, John W. Hartley, David Harvey, Bradley Holder, Andrew Hunter, Eric M. Huff, Andrew Hynous, Mathilde Jauzac, William C. Jones, Nikky Joyce, Duncan Kennedy, David Lagattuta, Jason S.-Y. Leung, Lun Li, Stephen Lishman, Thuy Vy T. Luu, Jacqueline E. McCleary, Johanna M. Nagy, C. Barth Netterfield, Emaad Paracha, Robert Purcaru, Susan F. Redmond, Jason D. Rhodes, Andrew Robertson, L. Javier Romualdez, Sarah Roth, Robert Salter, Jurgen Schmoll, Mohamed M. Shaaban, Roger Smith, Russell Smith, Sut Ieng Tam, Georgios N. Vassilakis

12 pages

In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) capsules containing 5 TB solid state data storage, plus a GNSS receiver, Iridium transmitter, and parachute. Data from the telescope were copied to these, and two were dropped over Argentina. They drifted 61 km horizontally while they descended 32 km, but we predicted their descent vectors within 2.4 km: in this location, the discrepancy appears irreducible below 2 km because of high speed, gusty winds and local topography. The capsules then reported their own locations to within a few metres. We recovered the capsules and successfully retrieved all of superBIT's data - despite the telescope itself being later destroyed on landing.

L.Zadorozhna, A.Tugay, O.Prikhodko, D.Malyshev, Y.Sahai, D.Savchenko, N.Pulatova

6 pages, 2 figures, 2 tables

The X-ray spectrum of the Coma galaxy cluster was studied using the data from the XMM-Newton observatory. We combined 7 observations performed with the MOS camera of XMM-Newton in the 40'x 40' region centered at the Coma cluster. The analyzed observations were performed in 2000-2005 and have a total duration of 196 ksec. We focus on the analysis of the MOS camera spectra due to their lower affection by strong instrumental line-like background. The obtained spectrum was fitted with a model including contributions from the Solar system/Milky Way hot plasma and a power law X-ray background. The contribution of the instrumental background was modeled as a power law (not convolved with the effective area) and a number of Gaussian lines. The contribution from the Coma cluster was modeled with a single-temperature hot plasma emission. In addition, we searched for possible non-thermal radiation present in the vicinity of the center of the Coma cluster, originating e.g. from synchrotron emission of relativistic electrons on a turbulent magnetic field. We compared the results with previous works by other authors and spectra obtained from other instruments that operate in the similar energy range of 1-10 keV. Careful and detailed spectrum analysis shall be a necessary contribution to our future work - searching for axion-like particles' manifestations in the Coma cluster.

Ximena C. Abrevaya, Douglas Galante, Paula M. Tribelli, Oscar J. Oppezzo, Felipe Nobrega, Gabriel G. Araujo, Fabio Rodrigues, Petra Odert, Martin Leitzinger, Martiniano M. Ricardi, Maria Eugenia Varela, Tamires Gallo, Jorge Sanz-Forcada, Ignasi Ribas, Gustavo F. Porto de Mello, Florian Rodler, 1 Maria Fernanda Cerini, Arnold Hanslmeier, Jorge E. Horvath

Halite (NaCl mineral) has exhibited the potential to preserve microorganisms for millions of years on Earth. This mineral was also identified on Mars and in meteorites. In this study, we investigated the potential of halite crystals to protect microbial life forms on the surface of an airless body (e.g., meteorite), for instance, during a lithopanspermia process (interplanetary travel step) in the early Solar System. To investigate the effect of the radiation of the young Sun on microorganisms, we performed extensive simulation experiments by employing a synchrotron facility. We focused on two exposure conditions: vacuum (low Earth orbit, 10^{-4}Pa) and vacuum-ultraviolet (VUV) radiation (range 57.6 - 124 nm, flux 7.14 W m^{-2}), with the latter representing an extreme scenario with high VUV fluxes comparable to the amount of radiation of a stellar superflare from the young Sun. The stellar VUV parameters were estimated by using the very well-studied solar analog of the young Sun, k^{1}Cet. To evaluate the protective effects of halite, we entrapped a halophilic archaeon (Haloferax volcanii) and a non-halophilic bacterium (Deinococcus radiodurans) in laboratory-grown halite. Control groups were cells entrapped in salt crystals (mixtures of different salts and NaCl) and non-trapped (naked) cells, respectively. All groups were exposed either to vacuum alone or to vacuum plus VUV. Our results demonstrate that halite can serve as protection against vacuum and VUV radiation, regardless of the type of microorganism. In addition, we found that the protection is higher than provided by crystals obtained from mixtures of salts. This extends the protective effects of halite documented in previous studies and reinforces the possibility to consider the crystals of this mineral as potential preservation structures in airless bodies or as vehicles for the interplanetary transfer of microorganisms.

Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for acceleration within the plasma's coherent structures. Here, we investigate how turbulent intermittency models connect statistical fluctuations in turbulence to regions of high dissipation. This connection is established by employing a generalized She-Leveque model to describe the exponents $\zeta_p$ for the structure functions $S^p \propto l^{\zeta_p}$. The fitting of the scaling exponents provide us with a measure of the co-dimension of the dissipative structures, and we subsequently measure their filling fraction. We perform our analysis for a range of magnetizations $\sigma$ and magnetic field fluctuations ${\delta B_0}/{B_0}$. We find that increasing the values of $\sigma$ and ${\delta B_0}/{B_0}$ allows the cascade to break sheets into smaller regions of dissipation that resemble chains of plasmoids. However, as their dissipation increases, the dissipative regions become less volume filling. With this work we aim to inform future turbulent acceleration theories that incorporate particle energization from interactions with coherent structures within relativistic turbulence.

Jongho Park, Motoki Kino, Hiroshi Nagai, Masanori Nakamura, Keiichi Asada, Minchul Kam, Jeffrey A. Hodgson

18 pages, 18 figures, accepted for publication in Astronomy and Astrophysics

Highly collimated relativistic jets are a defining feature of certain active galactic nuclei (AGN), yet their formation mechanism remains elusive. Previous observations and theoretical models have proposed that the ambient medium surrounding the jets could exert pressure, playing a crucial role in shaping the jets. However, direct observational confirmation of such a medium has been lacking. In this study, we present very long baseline interferometric (VLBI) observations of 3C 84 (NGC 1275), located at the center of the Perseus Cluster. Through monitoring observations with the Very Long Baseline Array (VLBA) at 43 GHz, a jet knot was detected to have been ejected from the sub-parsec scale core in the late 2010s. Intriguingly, this knot propagated in a direction significantly offset from the parsec-scale jet direction. To delve deeper into the matter, we employ follow-up VLBA 43 GHz observations, tracing the knot's trajectory until the end of 2022. We discovered that the knot abruptly changed its trajectory in the early 2020s, realigning itself with the parsec-scale jet direction. Additionally, we present results from an observation of 3C 84 with the Global VLBI Alliance (GVA) at 22 GHz, conducted near the monitoring period. By jointly analyzing the GVA 22 GHz image with a VLBA 43 GHz image observed about one week apart, we generated a spectral index map, revealing an inverted spectrum region near the edge of the jet where the knot experienced deflection. These findings suggest the presence of a dense, cold ambient medium characterized by an electron density exceeding $\sim10^5\ {\rm cm^{-3}}$, which guides the jet's propagation on parsec scales and significantly contributes to the overall shaping of the jet.

Xunchuan Liu, Tie Liu, Lei Zhu, Guido Garay, Hong-Li Liu, Paul Goldsmith, Neal Evans, Kee-Tae Kim, Sheng-Yuan Liu, Fengwei Xu, Xing Lu, Anandmayee Tej, Xiaofeng Mai, Leonardo Bronfman, Shanghuo Li, Diego Mardones, Amelia Stutz, Ken'ichi Tatematsu, Ke Wang, Qizhou Zhang, Sheng-Li Qin, Jianwen Zhou, Qiuyi Luo, Siju Zhang, Yu Cheng, Jinhua He, Qilao Gu, Ziyang Li, Zhenying Zhang, Suinan Zhang, Anindya Saha, Lokesh Dewangan, Patricio Sanhueza, Zhiqiang Shen

9 figures, 4 tables, accepted by RAA

This paper presents an overview of the QUARKS survey, which stands for `Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures'. The QUARKS survey is observing 139 massive clumps covered by 156 pointings at ALMA Band 6 ($\lambda\sim$ 1.3 mm). In conjunction with data obtained from the ALMA-ATOMS survey at Band 3 ($\lambda\sim$ 3 mm), QUARKS aims to carry out an unbiased statistical investigation of massive star formation process within protoclusters down to a scale of 1000 au. This overview paper describes the observations and data reduction of the QUARKS survey, and gives a first look at an exemplar source, the mini-starburst Sgr B2(M). The wide-bandwidth (7.5 GHz) and high-angular-resolution (~0.3 arcsec) observations of the QUARKS survey allow to resolve much more compact cores than could be done by the ATOMS survey, and to detect previously unrevealed fainter filamentary structures. The spectral windows cover transitions of species including CO, SO, N$_2$D$^+$, SiO, H$_{30}\alpha$, H$_2$CO, CH$_3$CN and many other complex organic molecules, tracing gas components with different temperatures and spatial extents. QUARKS aims to deepen our understanding of several scientific topics of massive star formation, such as the mass transport within protoclusters by (hub-)filamentary structures, the existence of massive starless cores, the physical and chemical properties of dense cores within protoclusters, and the feedback from already formed high-mass young protostars.

Kuan Wang (UMich), Philip Mansfield (Stanford), Dhayaa Anbajagane (UChicago), Camille Avestruz (UMich)

12+3 pages, 5+2 figures. Fig. 4 and 5 are the main figures. To be submitted to MNRAS, comments welcome

Anisotropy properties -- halo spin, shape, position offset, velocity offset, and orientation -- are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of haloes. These properties reflect the dynamical state of haloes, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behaviour of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behaviour of halo properties. This study paves the way towards a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of haloes.

The Event Horizon Telescope Collaboration

30 pages, 13 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

We present Event Horizon Telescope (EHT) 1.3 mm measurements of the radio source located at the position of the supermassive black hole Sagittarius A* (Sgr A*), collected during the 2017 April 5--11 campaign. The observations were carried out with eight facilities at six locations across the globe. Novel calibration methods are employed to account for Sgr A*'s flux variability. The majority of the 1.3 mm emission arises from horizon scales, where intrinsic structural source variability is detected on timescales of minutes to hours. The effects of interstellar scattering on the image and its variability are found to be subdominant to intrinsic source structure. The calibrated visibility amplitudes, particularly the locations of the visibility minima, are broadly consistent with a blurred ring with a diameter of $\sim$50 $\mu$as, as determined in later works in this series. Contemporaneous multi-wavelength monitoring of Sgr A* was performed at 22, 43, and 86 GHz and at near infrared and X-ray wavelengths. Several X-ray flares from Sgr A* are detected by Chandra, one at low significance jointly with Swift on 2017 April 7 and the other at higher significance jointly with NuSTAR on 2017 April 11. The brighter April 11 flare is not observed simultaneously by the EHT but is followed by a significant increase in millimeter flux variability immediately after the X-ray outburst, indicating a likely connection in the emission physics near the event horizon. We compare Sgr A*'s broadband flux during the EHT campaign to its historical spectral energy distribution and find both the quiescent and flare emission are consistent with its long-term behaviour.

The Event Horizon Telescope Collaboration

24 pages, 6 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A$^*$), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of $\lambda=1.3\,{\rm mm}$. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of $51.8 \pm 2.3$\,\uas (68\% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A$^*$ are consistent with the expected appearance of a Kerr black hole with mass ${\sim}4 \times 10^6\,{\rm M}_\odot$, which is inferred to exist at this location based on previous infrared observations of individual stellar orbits as well as maser proper motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination ($i > 50^\circ$), as well as non-spinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way galaxy, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of $10^3-10^5$ gravitational radii to event horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87$^*$ shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.

Tyler Schmaltz, Yue Hu, A. Lazarian

11 pages, 11 figures, submitted to MNRAS

Magnetic fields are a defining yet enigmatic aspect of the interstellar medium (ISM), with their three-dimensional mapping posing a substantial challenge. In this study, we harness the innovative Velocity Gradient Technique (VGT), underpinned by magnetohydrodynamic (MHD) turbulence theories, to elucidate the magnetic field structure by applying it to the atomic neutral hydrogen (HI) emission line and the molecular tracer $^{12}$CO. We construct the tomography of the magnetic field in the low-mass star-forming region L1688, utilizing two approaches: (1) VGT-HI combined with the Galactic rotational curve, and (2) stellar polarization paired with precise star parallax measurements. Our analysis reveals that the magnetic field orientations deduced from stellar polarization undergo a distinct directional change in the vicinity of L1688, providing evidence that the misalignment between VGT-HI and stellar polarization stems from the influence of the molecular cloud's magnetic field on the polarization of starlight. When comparing VGT-$^{12}$CO to stellar polarization and Planck polarization data, we observe that VGT-$^{12}$CO effectively reconciles the misalignment noted with VGT-HI, showing statistical alignment with Planck polarization measurements. This indicates that VGT-$^{12}$CO could be integrated with VGT-HI, offering vital insights into the magnetic fields of molecular clouds, thereby enhancing the accuracy of our 3D magnetic field reconstructions.

A. Yu. Kirichenko, S. V. Zharikov, A. V. Karpova, E. Fonseca, D. A. Zyuzin, Yu. A. Shibanov, E. A. López, M. R. Gilfanov, A. Cabrera-Lavers, S. Geier, F. A. Dong, D. C. Good, J. W. McKee, B. W. Meyers, I. H. Stairs, M. A. McLaughlin, J. K. Swiggum

Accepted for publication in MNRAS

PSR J1641+8049 is a 2 ms black widow pulsar with the 2.2 h orbital period detected in the radio and $\gamma$-rays. We performed new phase-resolved multi-band photometry of PSR J1641+8049 using the OSIRIS instrument at the Gran Telescopio Canarias. The obtained data were analysed together with the new radio-timing observations from the Canadian Hydrogen Intensity Mapping Experiment (CHIME), the X-ray data from the Spectrum-RG/eROSITA all-sky survey, and all available optical photometric observations. An updated timing solution based on CHIME data is presented, which accounts for secular and periodic modulations in pulse dispersion. The system parameters obtained through the light curve analysis, including the distance to the source 4.6-4.8 kpc and the orbital inclination 56-59 deg, are found to be consistent with previous studies. However, the optical flux of the source at the maximum brightness phase faded by a factor of $\sim$2 as compared to previous observations. Nevertheless, the face of the J1641+8049 companion remains one of the most heated (8000-9500 K) by a pulsar among the known black widow pulsars. We also report a new estimation on the pulsar proper motion of $\approx$2 mas yr$^{-1}$, which yields a spin down luminosity of $\approx$4.87$\times 10^{34}$ ergs s$^{-1}$ and a corresponding heating efficiency of the companion by the pulsar of 0.3-0.7. The pulsar was not detected in X-rays implying its X-ray-luminosity was <3 $\times$ 10$^{31}$ erg s$^{-1}$ at the date of observations.

Kohei Kurahara, Takuya Akahori, Aika Oki, Yuki Omiya, Kazuhiro Nakazawa

6 pages, 3 figures, Submitted to PASJ as Letter

Non-thermal components in the intra-cluster medium (ICM) such as turbulence, magnetic field, and cosmic rays imprint the past and current energetic activities of jets from active galactic neuclie (AGN) of member galaxies as well as disturbance caused by galaxy cluster mergers. Meter- and centimeter-radio observations of synchrotron radiation allow us to diagnose the nonthermal component. Here we report on our discovery of an unidentified diffuse radio source, named the Flying Fox, near the center of the Abell 1060 field. The Flying Fox has an elongated ring-like structure and a central bar shape, but there is no obvious host galaxy. The average spectral index of the Flying Fox is -1.4, which is steeper than those for radio sources seen at meter wavelength. We discussed the possibilities of radio lobes, phoenixes, radio halos and relics, and Odd Radio Circle (ORC). In conclusion, the Flying Fox is not clearly explained by known radio sources.

The Event Horizon Telescope Collaboration

65 pages, 35 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fraction of the correlated flux density, reaching $\sim$100\% on some baselines. Through an exploration of simple geometric source models, we demonstrate that ring-like morphologies provide better fits to the Sgr A* data than do other morphologies with comparable complexity. We develop two strategies for fitting static geometric ring models to the time-variable Sgr A* data; one strategy fits models to short segments of data over which the source is static and averages these independent fits, while the other fits models to the full dataset using a parametric model for the structural variability power spectrum around the average source structure. Both geometric modeling and image-domain feature extraction techniques determine the ring diameter to be $51.8 \pm 2.3$ $\mu$as (68\% credible intervals), with the ring thickness constrained to have an FWHM between $\sim$30\% and 50\% of the ring diameter. To bring the diameter measurements to a common physical scale, we calibrate them using synthetic data generated from GRMHD simulations. This calibration constrains the angular size of the gravitational radius to be $4.8_{-0.7}^{+1.4}$ \mathrm{\mu as}, which we combine with an independent distance measurement from maser parallaxes to determine the mass of Sgr A* to be $4.0_{-0.6}^{+1.1} \times 10^6$ M$_{\odot}$.

We carried a detailed temporal and spectral study of the BL\,Lac by using the long-term \emph{Fermi}-LAT and \emph{Swift}-XRT/UVOT observations, during the period MJD\,59000-59943. The daily-binned $\gamma$-ray light curve displays a maximum flux of $1.74\pm 0.09\times 10^{-5} \rm ph\,cm^{-2}\,s^{-1}$ on MJD\,59868, which is the highest daily $\gamma$-ray flux observed from BL\,Lac. The $\gamma$-ray variability is characterised by power-spectral-density (PSD), r.m.s-flux relation and flux-distribution study. We find that power-law model fits the PSD with index $\sim 1$, which suggest for long memory process at work. The observed r.m.s.-flux relation exhibits a linear trend, which indicates that the $\gamma$-ray flux distribution follows a log-normal distribution. The skewness/Anderson-Darling test and histogram-fit reject the normality of flux distribution, and instead suggest that the flux distribution is log-normal distribution. The fractional-variability amplitude shows that source is more variable in X-ray band than in optical/UV/$\gamma$-ray bands. In order to obtain an insight into the underlying physical process, we extracted broadband spectra from different time periods of the lightcurve. The broadband spectra are statistically fitted with the convolved one-zone leptonic model with different forms of the particle energy distribution. We found that spectral energy distribution during different flux states can be reproduced well with the synchrotron, synchrotron-self-Compton and external-Compton emissions from a broken power-law electron distribution, ensuring equipartition condition. A comparison between the best fit physical parameters show that the variation in different flux-states are mostly related to increase in the bulk-Lorentz factor and spectral hardening of the particle distribution.

Seo-Won Chang, Christian Wolf, Christopher A. Onken, Michael S. Bessell

8 pages, 7 figures, Submitted to MNRAS

We report the discovery of a new Blue Large-Amplitude Pulsator, SMSS J184506-300804 (SMSS-BLAP-1) in Data Release 2 of the SkyMapper Southern Sky Survey. We conduct high-cadence photometric observations in the $u$ band to confirm a periodic modulation of the lightcurve. SMSS-BLAP-1 has a ~19-min pulsation period with an amplitude of 0.2 mag in u band, and is similar to the classical BLAPs found by OGLE. From spectroscopic observations with the Wide-Field Spectrograph on the ANU 2.3m telescope, we confirm it as a low-gravity BLAP: best-fit parameters from our spectral model grid are estimated as $T_\mathrm{eff}$ = 27,000 K and $\log g$ (cm s$^{-2}$) = 4.4. Remarkably, we find evidence of a periodic signal in the residual lightcurve that could hint at a non-radial pulsation mode, and an excess of Ca II K and Na I D absorption from potential circumstellar material.

Benedikt Eggemeier, Peter Hayman, Jens C. Niemeyer, Richard Easther

11 pages, 8 figures

The post-inflationary Universe can pass through a long epoch of effective matter-dominated expansion. This era may allow for both the parametric amplification of initial fluctuations and the gravitational collapse of inflaton perturbations. We perform first-of-their-kind high-resolution simulations that span the resonant phase and the subsequent gravitational collapse of the inflaton field by seguing from a full Klein-Gordon treatment of resonance to a computationally efficient Schr\"odinger-Poisson description that accurately captures the gravitational dynamics when most quanta are nonrelativistic. We consider a representative example in which resonance generates $\mathcal{O}(10^{-1})$ overdensities and gravitational collapse follows promptly as resonance ends. We observe the formation of solitonic cores inside inflaton halos and complex gravitational dynamics on scales of $10^{-27}\,\mathrm{m}$, greatly extending the possible scope of nonlinear post-inflationary gravitational dynamics.

Adrian J. Barker, Michael Efroimsky, Valeri V. Makarov, Dimitri Veras

9 pages, 2 figures, accepted for publication in MNRAS (14th Nov 2023)

The gas giant Kepler-1658b has been inferred to be spiralling into its sub-giant F-type host star Kepler-1658a (KOI-4). The measured rate of change of its orbital period is $\dot{P}_{\rm orb}=-131^{+20}_{-22}\mathrm{ms/yr}$, which can be explained by tidal dissipation in the star if its modified tidal quality factor is as low as $Q^{\,\prime}\approx 2.50\times {10}^{4}$. We explore whether this could plausibly be consistent with theoretical predictions based on applying up-to-date tidal theory in stellar models (varying stellar mass, age, and metallicity) consistent with our newly-derived observational constraints. In most of our models matching the combined constraints on the stellar effective temperature and radius, the dissipation in the star is far too weak, capable of providing $Q^{\,\prime}\gtrsim 10^9$, hence contributing negligibly to orbital evolution. Using only constraints on the stellar radius, efficient tidal dissipation sufficient to explain observations is possible due to inertial waves in the convective envelope during the sub-giant phase, providing $Q^{\,\prime}\sim 10^4$, but this period in the evolution is very short-lived (shorter than $10^2$ yrs in our models). We show that dissipation in the planet is capable of explaining the observed $\dot{P}_\mathrm{orb}$ only if the planet rotates non-synchronously. Tidally-induced pericentre precession is a viable explanation if the periastron argument is near $3\pi/2$ and the quadrupolar Love number is above 0.26. Further observations constraining the stellar and planetary properties in this system have the exciting potential to test tidal theories in stars and planets.

A. R. Yeates, M. A. Berger

23 pages, 8 figures, to appear as chapter of AGU book "Helicities in Geophysics, Astrophysics and Beyond", published by Wiley, ISBN 1119841682

Field line helicity measures the net linking of magnetic flux with a single magnetic field line. It offers a finer topological description than the usual global magnetic helicity integral, while still being invariant in an ideal evolution unless there is a flux of helicity through the domain boundary. In this chapter, we explore how to appropriately define field line helicity in different volumes in a way that preserves a meaningful topological interpretation. We also review the time evolution of field line helicity under both boundary motions and magnetic reconnection.

Shogo Nishiyama, Tomohiro Kara, Brian Thorsbro, Hiromi Saida, Yohsuke Takamori, Masaaki Takahashi, Takayuki Ohgami, Kohei Ichikawa, Rainer Schödel

15 pages, 6 figures, 10 tables, accepted for publication in Proceedings of the Japan Academy, Ser. B, Physical and Biological Sciences

The tremendous tidal force that is linked to the supermassive black hole (SMBH) at the center of our galaxy is expected to strongly subdue star formation in its vicinity. Stars within 1" from the SMBH thus likely formed further from the SMBH and migrated to their current positions. In this study, spectroscopic observations of the star S0-6/S10, one of the closest (projected distance from the SMBH of about 0.3") late-type stars were conducted. Using metal absorption lines in the spectra of S0-6, the radial velocity of S0-6 from 2014 to 2021 was measured, and a marginal acceleration was detected, which indicated that S0-6 is close to the SMBH. The S0-6 spectra were employed to determine its stellar parameters including temperature, chemical abundances ([M/H], [Fe/H], [alpha/Fe], [Ca/Fe], [Mg/Fe], [Ti/Fe]), and age. As suggested by the results of this study, S0-6 is very old (> ~10 Gyr) and has an origin different from that of stars born in the central pc region.

Jakob Beise (for the IceCube-Gen2 collaboration)

Presented at the XVIII International Conference on Topics in Astroparticle and Underground Physics (TAUP 2023)

Core-collapse supernovae (CCSNe) are among the most energetic processes in our Universe and are crucial for the understanding of the formation and chemical composition of the Universe. The precise measurement of the neutrino light curve from CCSNe is crucial to understanding the hydrodynamics and fundamental processes that drive CCSNe. The IceCube Neutrino Observatory has mass-independent sensitivity within the Milky Way and some sensitivity to the higher mass CCSNe in the Large and Small Magellanic clouds. The envisaged large-scale extension of the IceCube detector, IceCube-Gen2, opens the possibility for new sensor design and trigger concepts that could increase the number of neutrinos detected from a CCSNe burst compared to IceCube. In this contribution, we study how wavelength-shifting technology can be used in IceCube-Gen2 to measure the fast modulations of the neutrino signal due to standing accretion shock instabilities (SASI).

With the detection of black hole mergers by the LIGO gravitational wave telescope, there has been increasing interest in the possibility that dark matter may be in the form of solar mass primordial black holes. One of the predictions implicit in this idea is that compact clouds in the broad emission line regions of high redshift quasars will be microlensed, leading to changes in line structure and the appearance of new emission features. In this paper the effect of microlensing on the broad emission line region is reviewed by reference to gravitationally lensed quasar systems where microlensing of the emission lines can be unambiguously identified. It is then shown that although changes in Seyfert galaxy line profiles occur on timescales of a few years, they are too nearby for a significant chance that they could be microlensed, and are plausibly attributed to intrinsic changes in line structure. In contrast, in a sample of 53 high redshift quasars, 9 quasars show large changes in line profile at a rate consistent with microlensing. These changes occur on a timescale an order of magnitude too short for changes associated with the dynamics of the emission line region. The main conclusion of the paper is that the observed changes in quasar emission line profiles are consistent with microlensing by a population of solar mass compact bodies making up the dark matter, although other explanations like intrinsic variability are possible. Such bodies are most plausibly identified as primordial black holes.

Wanqiu He, Masayuki Akiyama, Motohiro Enoki, Kohei Ichikawa, Kohei Inayoshi, Nobunari Kashikawa, Toshihiro Kawaguchi, Yoshiki Matsuoka, Tohru Nagao, Masafusa Onoue, Taira Oogi, Andreas Schulze, Yoshiki Toba, Yoshihiro Ueda

38 pages, 26 figures. Submitted to ApJ. Data is available in this https URL

We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH masses of quasars are estimated using the width of the CIV 1549{\AA} line and the continuum luminosity at 1350{\AA}. To obtain the intrinsic broad-line AGN BHMF and ERDF, we correct for the incompleteness in the low-mass and/or low-Eddington-ratio ranges caused by the flux-limited selection. The resulting BHMF is constrained down to $\log M_{\rm BH}/M_{\odot}\sim7.5$. In comparison with broad-line AGN BHMFs at z=2 in the literature, we find that the number density of massive SMBHs peaks at higher redshifts, consistent with the "down-sizing" evolutionary scenario. Additionally, the resulting ERDF shows a negative dependence on BH mass, suggesting more massive SMBHs tend to accrete at lower Eddington ratios at z=4. With the derived intrinsic broad-line AGN BHMF, we also evaluate the active fraction of broad-line AGNs among the entire SMBH population at z=4. The resulting active fraction may suggest a positive dependence on BH mass. Finally, we examine the time evolution of broad-line AGN BHMF between z=4 and 6 through solving the continuity equation. The results suggest that the broad-line AGN BHMFs at z=4-6 only show evolution in their normalization, but with no significant changes in their shape.

Cai-Xia Qu, A-Li Luo, Rui Wang, Hugh R. A. Jones, Bing Du, Xiang-Lei Chen, You-Fen Wang

14 pages, 12 figures, accepted by ApJS

We provide a catalogue of atmospheric parameters for 1,806,921 cool dwarfs from Gaia DR3 which lie within the range covered by LAMOST cool dwarf spectroscopic parameters: 3200 K < T_{eff}< 4300 K, -0.8 < [M/H] < 0.2 dex, and 4.5 <log{g} < 5.5 dex. Our values are derived based on Machine Learning models trained with multi-band photometry corrected for dust. The photometric data comprises of optical from SDSS r, i, z bands, near-infrared from 2MASS J, H, K and mid-infrared from ALLWISE W1, W2. We used both random forest and LightGBM machine learning models and found similar results from both with an error dispersion of 68 K, 0.22 dex, and 0.05 dex for T_{eff}, [M/H], and log {g}, respectively. Assessment of the relative feature importance of different photometric colors indicated W1 -- W2 as most sensitive to both T_{eff} and log{g}, with J -- H most sensitive to [M/H]. We find that our values show a good agreement with APOGEE, but are significantly different to those provided as part of Gaia DR3.

Francesco Zagaria, Cathie J. Clarke, Richard A. Booth, Stefano Facchini, Giovanni P. Rosotti

13 pages, 4 figures; accepted for publication on ApJL

We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a factor of two. We combine our results with literature measurements of the dust surface density and grain size to determine the dust-to-gas ratio and Stokes number in the 67 au and 100 au rings. Our estimates suggest that particle clumping is taking place under the effect of streaming instability (SI) in the 100 au ring. Even though in the presence of external isotropic turbulence this process might be hindered, we provide evidence that turbulence is non-isotropic in both rings and likely originating from mechanisms (such as ambipolar diffusion) that could ease particle clumping under SI. Finally, we determine the mass accretion rate under the assumption that the disc is in steady state and turbulence regulates angular momentum transport. Our results are in tension with spectroscopic measurements and suggest that other mechanisms might be responsible for accretion, in qualitative agreement with the detection of a magneto-centrifugal wind in this system. Applying our method to larger samples can be used to statistically assess if SI is a viable mechanism to form planetesimals in bright rings.

Amy M. Jones, Rachael L. Beaton, Brian A. Cherinka, Karen L. Masters, Sara Lucatello, Aleksandar M. Diamond-Stanic, Sarah A. Bird, Michael R. Blanton, Katia Cunha, Emily E. Farr, Diane Feuillet, Peter M. Frinchaboy, Alex Hagen, Karen Kinemuchi, Britt Lundgren, Mariarosa L. Marinelli, Adam D. Myers, Alexandre Roman-Lopes, Ashley J. Ross, Jose R. Sanchez-Gallego, Sarah J. Schmidt, Jennifer Sobeck, Keivan G. Stassun, Jamie Tayar, Mariana Vargas-Magana, J. C. Wilson, Gail Zasowski

30 pages, 9 figures, accepted in PASP

The Sloan Digital Sky Survey (SDSS) is one of the largest international astronomy organizations. We present demographic data based on surveys of its members from 2014, 2015 and 2016, during the fourth phase of SDSS (SDSS-IV). We find about half of SDSS-IV collaboration members were based in North America, a quarter in Europe, and the remainder in Asia and Central and South America. Overall, 26-36% are women (from 2014 to 2016), up to 2% report non-binary genders. 11-14% report that they are racial or ethnic minorities where they live. The fraction of women drops with seniority, and is also lower among collaboration leadership. Men in SDSS-IV were more likely to report being in a leadership role, and for the role to be funded and formally recognized. SDSS-IV collaboration members are twice as likely to have a parent with a college degree, than the general population, and are ten times more likely to have a parent with a PhD. This trend is slightly enhanced for female collaboration members. Despite this, the fraction of first generation college students (FGCS) is significant (31%). This fraction increased among collaboration members who are racial or ethnic minorities (40-50%), and decreased among women (15-25%). SDSS-IV implemented many inclusive policies and established a dedicated committee, the Committee on INclusiveness in SDSS (COINS). More than 60% of the collaboration agree that the collaboration is inclusive; however, collaboration leadership more strongly agree with this than the general membership. In this paper, we explain these results in full, including the history of inclusive efforts in SDSS-IV. We conclude with a list of suggested recommendations based on our findings, which can be used to improve equity and inclusion in large astronomical collaborations, which we argue is not only moral, but will also optimize their scientific output.

Context: Radio wavelengths offer a unique possibility of tracing the total star-formation rate (SFR) in galaxies, both obscured and unobscured. To prob the dust-unbiased star-formation history, an accurate measurement of the radio luminosity function (LF) for star-forming galaxies (SFGs) is crucial. Aims: We make use of an SFG sample (5915 sources) from the Very Large Array (VLA) COSMOS 3 GHz data to perform a new modeling of the radio LF. By integrating the analytical LF we aim to calculate the cosmic SFR density (SFRD) history since $z\sim5$. Methods: For the first time, we use both models of the pure luminosity evolution (PLE) and joint luminosity+density evolution (LADE) to fit the LFs directly to the radio data using a full maximum likelihood analysis, considering the sample completeness correction. We also incorporate the updated observations on local radio LFs and radio source counts into the fitting process to obtain additional constraints. Results: We find that the PLE model is not applicable to describe the evolution of the radio LF at high redshift ($z>2$). By construct, our LADE models can successfully fit a large amount of data on radio LFs and source counts of SFGs from recent observations. Our SFRD curve shows a good fit to the SFRD points derived by previous radio estimates. In view of that our radio LFs are not biased in contrast to previous studies performed by fitting the $1/V_{\rm max}$ LF points, our SFRD results should be an improvement to these previous estimates. Below $z<2$, our SFRD matches well with the multi-wavelength compilation from Madau \& Dickinson, while our SFRD turns over at a slightly higher redshift ($2<z<2.5$) and falls more rapidly out to high redshift.

The proliferation of space debris in LEO has become a major concern for the space industry. With the growing interest in space exploration, the prediction of potential collisions between objects in orbit has become a crucial issue. It is estimated that, in orbit, there are millions of fragments a few millimeters in size and thousands of inoperative satellites and discarded rocket stages. Given the high speeds that these fragments can reach, even fragments a few millimeters in size can cause fractures in a satellite's hull or put a serious crack in the window of a space shuttle. The conventional method proposed by Akella and Alfriend in 2000 remains widely used to estimate the probability of collision in short-term encounters. Given the small period of time, it is assumed that, during the encounter: (1) trajectories are represented by straight lines with constant velocity; (2) there is no velocity uncertainty and the position exhibits a stationary distribution throughout the encounter; and (3) position uncertainties are independent and represented by Gaussian distributions. This study introduces a novel derivation based on first principles that naturally allows for tight and fast upper and lower bounds for the probability of collision. We tested implementations of both probability and bound computations with the original and our formulation on a real CDM dataset used in ESA's Collision Avoidance Challenge. Our approach reduces the calculation of the probability to two one-dimensional integrals and has the potential to significantly reduce the processing time compared to the traditional method, from 80% to nearly real-time.

It has recently been suggested that black holes (BH) may exhibit growth of their mass with time, so that their mass is proportional to the cosmological scale factor to the power $n$, with suggested values $n \sim 3$ for supermassive BHs in elliptical galaxies. Here we test these predictions with stellar mass BHs in X-ray binaries using their masses and ages. We perform two sets of tests to assess the compatible values of $n$. First, we assume that no compact object grows over the Tolman-Oppenheimer-Volkof limit which marks the borderline between neutron stars and BHs. We show that half of the BHs would be born with a mass below this limit if $n=3$ applies. The possibility that all BHs were born above the limit is rejected at $4\,\sigma$ if $n=3$ applies. In the second test, we assume that masses of BHs at their formation stay the same over cosmic history. We compare the mass distribution of the youngest BHs, which could have not grown yet, to their older counterparts. Distributions are compatible for $n = -0.8^{+1.2}_{-4.5}$, with $n=3$ excluded with $87\,\%$ confidence. This result may be biased, as massive BHs tend to have a massive companion. Correcting for this bias yields $n\sim 0$. We conclude that mass and age estimates of stellar mass BHs are incompatible with cosmological growth with $n \sim 3$ and favor their mass not changing with time.

Natalia Inostroza-Pino, Valerio Lattanzi, C. Zachary Palmer, Ryan C. Fortenberry, Diego Mardones, Paola Caselli, Oko E. Godwin, Timothy J. Lee

15 pages, 2 figures and 1 appendix table; published in Molecular Physics (Tim Lee Memorial Issue)

The synergy between high-resolution rotational spectroscopy and quantum-chemical calculations is essential for exploring future detection of molecules, especially when spectroscopy parameters are not available yet. By using highly correlated ab initio quartic force fields (QFFs) from explicitly correlated coupled cluster theory, a complete set of rotational constants and centrifugal distortion constants for D$_2$CS and cis/trans-DCSD isomers have been produced. Comparing our new ab initio results for D$_2$CS with new rotational spectroscopy laboratory data for the same species, the accuracy of the computed B and C rotational constants is within 0.1% while the A constant is only slightly higher. Additionally, quantum chemical vibrational frequencies are also provided, and these spectral reference data and new experimental rotational lines will provide additional references for potential observation of these deuterated sulfur species with either ground-based radio telescopes or space-based infrared observatories.

Soren Ramdhanie, Yjan A. Gordon, Heinz Andernach, Eric J. Hooper, Brianna Sampson

5 pages, 1 Table. Catalog data available from this https URL

Giant Radio Galaxies (GRGs) are Active Galactic Nuclei (AGN) with radio emission that extends over projected sizes $>0.7\,$Mpc. The large angular sizes associated with GRGs complicate their identification in radio survey images using traditional source finders. In this Note, we use DRAGNhunter, an algorithm designed to find double-lobed radio galaxies, to search for GRGs in the Faint Images of the Radio Sky at Twenty cm survey (FIRST). Radio and optical images of identified candidates are visually inspected to confirm their authenticity, resulting in the discovery of $63$ previously unreported GRGs.

Petri J. Käpylä (KIS)

10 pages, 4 figures, to appear in IAU Symposium 365 proceedings (A. Getling & L. Kitchatinov, eds)

The canonical undestanding of stellar convection has recently been put under doubt due to helioseismic results and global 3D convection simulations. This "convective conundrum" is manifested by much higher velocity amplitudes in simulations at large scales in comparison to helioseismic results, and the difficulty in reproducing the solar differential rotation and dynamo with global 3D simulations. Here some aspects of this conundrum are discussed from the viewpoint of hydrodynamic Cartesian 3D simulations targeted at testing the rotational influence and surface forcing on deep convection. More specifically, the dominant scale of convection and the depths of the convection zone and the weakly subadiabatic -- yet convecting -- Deardorff zone are discussed in detail.

Silvia Celli, Andreas Specovius, Stefano Menchiari, Alison Mitchell, Giovanni Morlino

11 pages, 16 figures, comments are welcome

Context. Star clusters constitute a relevant part of the stellar population in our Galaxy. The feedback processes they exert on the interstellar medium impact multiple physical processes, from the chemical to the dynamical evolution of the Galaxy. In addition, young and massive stellar clusters might act as efficient particle accelerators, possibly contributing to the production of cosmic rays. Aims. We aim at evaluating the wind luminosity driven by the young (< 30 Myr) Galactic open stellar clusters observed by the Gaia space mission, which is crucial to determine the energy channeled into accelerated particles. Methods. To this extent, we develop a method relying on the number, magnitude and line-of-sight extinction of the stars observed per cluster. Assuming that the stellar mass function follows a Kroupa mass distribution, and accounting for the maximum stellar mass allowed by both the parent cluster age and mass, we conservatively estimate the mass and wind luminosity of 387 local clusters within the second data release of Gaia. Results. We compare the results of our computation with recent estimations of young cluster masses. With respect to these, we provide a sample three times more abundant, particularly above a few thousand solar masses, which is of the utmost relevance for predicting the gamma-ray emission resulting from the interaction of accelerated particles. In fact, the cluster wind luminosity distribution we obtain is found to extend up to 3 x 10^38 erg/s, a promising feature in terms of potential particle acceleration scenarios.

We propose a novel idea for the coherent intense millisecond radio emission of cosmic fast radio bursts (FRBs), which have recently been identified with flares from a magnetar. Motivated by the conventional paradigm of Type III solar radio bursts, we will explore the emission of coherent plasma line radiation at the electron plasma frequency and its harmonic as potential candidates of the coherent FRB emissions associated with magnetar flares. We discuss the emissions region parameters in relativistic strongly magnetized plasmas consisting of electrons, positrons and protons. The goal is to make observable predictions of this model to confront the multi-wavelength observations of FRBs from magnetars. These results will impact both observational radio astronomy and space-based astrophysics

D. A. Zyuzin, A. Yu. Kirichenko, A. V. Karpova, Yu. A. Shibanov, S. V. Zharikov, M. R. Gilfanov, C. Perez Tórtola

10 pages, 11 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

An unidentified $\gamma$-ray source 4FGL J1838.2+3223 has been proposed as a pulsar candidate. We present optical time-series multi-band photometry of its likely optical companion obtained with the 2.1-m telescope of Observatorio Astron\'omico Nacional San Pedro M\'artir, Mexico. The observations and the data from the Zwicky Transient Facility revealed the source brightness variability with a period of $\approx$4.02 h likely associated with the orbital motion of the binary system. The folded light curves have a single sine-like peak per period with an amplitude of about three magnitude accompanied by fast sporadic flares up to one magnitude level. We reproduce them modelling the companion heating by the pulsar. As a result, the companion side facing the pulsar is strongly heated up to 11300$\pm$400 K, while the temperature of its back side is only 2300$\pm$700 K. It has a mass of 0.10$\pm$0.05 ${\rm M}_\odot$ and underfills its Roche lobe with a filling factor of $0.60^{+0.10}_{-0.06}$. This implies that 4FGL J1838.2+3223 likely belongs to the `spider' pulsar family. The estimated distance of $\approx$3.1 kpc is compatible with Gaia results. We detect a flare from the source in X-rays and ultraviolet using Swift archival data and another one in X-rays with the eROSITA all-sky survey. Both flares have X-ray luminosity of $\sim$10$^{34}$ erg s$^{-1}$ which is two orders of magnitude higher than the upper limit in quiescence obtained from eROSITA assuming spectral shape typical for spider pulsars. If the spider interpretation is correct, these flares are among the strongest flares observed from non-accreting spider pulsars.

William Giarè, Adrià Gómez-Valent, Eleonora Di Valentino, Carsten van de Bruck

35 pages, 16 figures, 5 tables

We study scatter-like interactions between neutrinos and dark matter in light of different combinations of temperature, polarization and lensing data released by three independent CMB experiments - the Planck satellite, the Atacama Cosmology Telescope (ACT), and the South Pole Telescope (SPT) - in conjunction with Baryon Acoustic Oscillation (BAO) measurements. We apply two different statistical methodologies. Alongside the usual marginalization technique, we cross-check all the results through a Profile Likelihood analysis. As a first step, working under the assumption of massless neutrinos, we perform a comprehensive (re-)analysis aimed at assessing the validity of some recent results hinting at a mild preference for non-vanishing interactions from small-scale CMB data. We find compelling resilience in the results already documented in the literature, confirming that interactions with a strength $u_{\nu\rm{DM}} \sim 10^{-5} - 10^{-4}$ appear to be globally favored by ACT (both alone and in combination with Planck). This result is corroborated by the inclusion of additional data, such as the recent ACT-DR6 lensing likelihood, as well as by the Profile Likelihood analysis. Interestingly, a fully consistent preference for interactions emerges from SPT, as well, although it is weaker than the one obtained from ACT. As a second step, we repeat the same analysis considering neutrinos as massive particles. Despite the larger parameter space, all the hints pointing towards interactions are confirmed also in this more realistic case. In addition, we report a very mild preference for interactions in Planck+BAO alone (not found in the massless case) which aligns with small-scale data. While this latter result is not fully confirmed by the Profile Likelihood analysis, the profile distribution does confirm that interactions are not disfavoured by Planck.

Ana Ennis, Johanna Hartke, Magda Arnaboldi, Claudia Pulsoni, Fuyan Bian, Chiara Spiniello

Submitted to the Proceedings of IAU Symposium 384: Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics. Eds: O. De Marco, A. Zijlstra, R. Szczerba

Planetary nebulae (PNe) are essential tracers of the kinematics of the diffuse halo and intracluster light where stellar spectroscopy is unfeasible, due to their strong emission lines. However, that is not all they can reveal about the underlying stellar population. In recent years, it has also been found that PNe in the metal-poor halos of galaxies have different properties (specific frequency, luminosity function), than PNe in the more metal-rich galaxy centers. A more quantitative understanding of the role of age and metallicity in these relations would turn PNe into valuable stellar-population tracers. In order to do that, a full characterization of PNe in regions where the stellar light can also be analysed in detail is necessary. In this work, we make use of integral-field spectroscopic data covering the central regions of galaxies, which allow us to measure both stellar ages and metallicities as well as to detect PNe. This analysis is fundamental to calibrate PNe as stellar population tracers and to push our understanding of galaxy properties at unprecedented galactocentric distances.

Some of the most important information on a radio pulsar is derived from its average pulse profile. Many early pulsar studies were necessarily based on only few such profiles. There, discrete profile components were linked to emission mechanism models for individual stars through human interpretation. For the population as a whole, profiles morphology must reflect the geometry and overall evolution of the radio emitting regions. The problem, however, is that this population is becoming too large for intensive studies of all sources individually. Moreover, connecting profiles from a large collection of pulsars rapidly becomes cumbersome. In this article, we present ToPP, the first-ever unsupervised method to sort pulsars by profile-shape similarity, using graph topology. We apply ToPP to the publicly available European Pulsar Network profile database, providing the first organised visual overview of multi-frequency profiles representing 90 individual pulsars. We find discrete evolutionary tracks, varying from simple, single component profiles at all frequencies, towards diverse mixtures of more complex profiles with frequency evolution. The profile evolution is continuous, extending out to millisecond pulsars, and does not fall in sharp classes. We interpret the profiles as a mixture of pulsar core/cone emission type, spin-down energetics, and the line-of-sight impact angle towards the magnetic axis. We show how ToPP can systematically classify sources into the Rankin empirical profile scheme. ToPP forms one of the key unsupervised methods that will be essential to explore upcoming pulsar census data such as expected by the Square Kilometer Array.

The subtle influence of gravitational waves on the apparent positioning of celestial bodies offers novel observational windows. We calculate the expected astrometric signal induced by an isotropic Stochastic Gravitational Wave Background (SGWB) in the short distance limit. Our focus is on the resultant proper motion of Solar System objects, a signal on the same time scales addressed by Pulsar Timing Arrays (PTA). We derive the corresponding astrometric deflection patterns, finding that they manifest as distinctive dipole and quadrupole correlations, or in some cases, may not be present. Our analysis encompasses both Einsteinian and non-Einsteinian polarisations. We estimate the upper limits for the amplitude of a scale-invariant SGWB that could be obtained by tracking the proper motions of large numbers of solar system objects such as asteroids. With the Gaia satellite and the Vera C. Rubin Observatory poised to track an extensive sample of asteroids-ranging from $O(10^5)$ to $O(10^6)$, we highlight the significant future potential for similar surveys to contribute to our understanding of the SGWB.

Kasturi Shah, Gregory P. Chini, Colm-cille P. Caulfield, Pascale Garaud

25 pages, 1 figure

Quantifying transport by strongly stratified turbulence in low Prandtl number ($Pr$) fluids is critically important for the development of better models for the structure and evolution of stellar interiors. Motivated by recent numerical simulations showing strongly anisotropic flows suggestive of scale-separated dynamics, we perform a multiscale asymptotic analysis of the governing equations. We find that, in all cases, the resulting slow-fast system naturally takes a quasilinear form. Our analysis also reveals the existence of several distinct dynamical regimes depending on the emergent buoyancy Reynolds and P\'eclet numbers, $Re_b = \alpha^2 Re$ and $Pe_b = Pr Re_b$, respectively, where $\alpha$ is the aspect ratio of the large-scale turbulent flow structures, and $Re$ is the outer scale Reynolds number. Scaling relationships relating the aspect ratio, the characteristic vertical velocity, and the strength of the stratification (measured by the Froude number $Fr$) naturally emerge from the analysis. When $Pe_b \ll \alpha$, the dynamics at all scales is dominated by buoyancy diffusion, and our results recover the scaling laws empirically obtained from direct numerical simulations by Cope et al. (2020). For $Pe_b \ge O(1)$, diffusion is negligible (or at least subdominant) at all scales and our results are consistent with those of Chini et al. (2022) for strongly stratified geophysical turbulence at $Pr = O(1)$.Finally, we have identified a new regime for $\alpha \ll Pe_b \ll 1$, in which slow, large scales are diffusive while fast, small scales are not. We conclude by presenting a map of parameter space that clearly indicates the transitions between isotropic turbulence, non-diffusive stratified turbulence, diffusive stratified turbulence and viscously-dominated flows.

Jorge L. Cervantes-Cota, Salvador Galindo-Uribarri, George F. Smoot

64 pages, 14 figures, 3 tables. Submitted to Universe

One of main sources of uncertainty in modern cosmology is the present rate of the Universe's expansion, H0, called the Hubble Constant. Once again different observational techniques bring different results causing a new "Hubble tension". In the present work we review the historical roots of Hubble constant from the beginning of the XX century, when modern cosmology started, to the present. We develop the arguments that gave rise to the importance of measuring the expansion of the Universe, its discovery and describing the different pioneering works to measure it. There has been a long dispute on this matter, even in the present epoch that is marked by high-tech instrumentation and therefore resulting in smaller uncertainties in the relevant parameters. It is again now necessary to conduct a careful and critical revision of the different methods, before one invokes new physics to solve the so-called Hubble tension.

Jarred Novak, Nickolas Solomey, Brooks Hartsock, Brian Doty, Jonathan Folkerts

Provisional Patent submitted and submitted to JINST for publication

Low energy anti-neutrinos detected from reactors or other sources have typically used the conversion of an anti-neutrino on hydrogen, producing a positron and a free neutron. This neutron is subsequently captured on a secondary element with a large neutron capture cross section such as Gadolinium or Cadmium. We have studied the anti-neutrino conversion and suggest other elements that have a comparable cross section for anti-neutrino reactions. With most neutron captures on Gadolinium, it is possible to get two or three delayed gamma signals of known energy to occur. With today's fast electronics, this leads to the possibility of having a triple delayed coincidence using the positron annihilation on atomic shell electrons as the starting signal. We have also found an isotope of Tungsten, $^{183}$W that offers a large anti-neutrino interaction cross section of $1.19 \times 10^{-46}$ m$^2$ and an anti-neutrino threshold energy of 2.094 MeV. This reaction makes a nuclear m1 excited state of $^{183}$Ta$^*$ that emits a signature secondary gamma pulse of 73 keV with a 107 ns half-life. This offers a new delayed coincidence technique that can be used to cleanly identify anti-neutrinos with less shielding with the added advantage of the anti-neutrino threshold shifting down to low energies.

Gravitational waves (GW) from chirping binary black holes (BBHs) provide unique opportunities to test general relativity (GR) in the strong-field regime. However, testing GR can be challenging when incomplete physical modeling of the expected signal gives rise to systematic biases. In this study, we investigate the potential influence of wave effects in gravitational lensing (which we refer to as microlensing) on tests of GR using GWs for the first time. We utilize an isolated point-lens model for microlensing with the lens mass ranging from $10-10^5~$M$_\odot$ and base our conclusions on an astrophysically motivated population of BBHs in the LIGO-Virgo detector network. Our analysis centers on two theory-agnostic tests of gravity: the inspiral-merger-ringdown consistency test (IMRCT) and the parameterized tests. Our findings reveal two key insights: First, microlensing can significantly bias GR tests, with a confidence level exceeding $5\sigma$. Notably, substantial deviations from GR $(\sigma > 3)$ tend to align with a strong preference for microlensing over an unlensed signal, underscoring the need for microlensing analysis before claiming any erroneous GR deviations. Nonetheless, we do encounter scenarios where deviations from GR remain significant ($1 < \sigma < 3$), yet the Bayes factor lacks the strength to confidently assert microlensing. Second, deviations from GR correlate with pronounced interference effects, which appear when the GW frequency ($f_\mathrm{GW}$) aligns with the inverse time delay between microlens-induced images ($t_\mathrm{d}$). These false deviations peak in the wave-dominated region and fade where $f_\mathrm{GW}\cdot t_\mathrm{d}$ significantly deviates from unity. Our findings apply broadly to any microlensing scenario, extending beyond specific models and parameter spaces, as we relate the observed biases to the fundamental characteristics of lensing.

We present two dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein-Greene-Kruskal (BGK) modes, which are exact nonlinear steady-state solutions of the Vlasov-Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 2048^2 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Kim V. Berghaus, Tanvi Karwal, Vivian Miranda, Thejs Brinckmann

27 pages, 16 figures, 3 tables

In this work, we quantify the cosmological signatures of dark energy radiation -- a novel description of dark energy, which proposes that the dynamical component of dark energy is comprised of a thermal bath of relativistic particles sourced by thermal friction from a slowly rolling scalar field. For a minimal model with particle production emerging from first principles, we find that the abundance of radiation sourced by dark energy can be as large as $\Omega_{\text{DER}} = 0.03$, exceeding the bounds on relic dark radiation by three orders of magnitude. Although the background and perturbative evolution of dark energy radiation is distinct from Quintessence, we find that current and near-future cosmic microwave background and supernova data will not distinguish these models of dark energy. We also find that our constraints on all models are dominated by their impact on the expansion rate of the Universe. Considering extensions that allow the dark radiation to populate neutrinos, axions, and dark photons, we evaluate the direct detection prospects of a thermal background comprised of these candidates consistent with cosmological constraints on dark energy radiation. Our study indicates that a resolution of $\sim 6 \, \text{meV}$ is required to achieve sensitivity to relativistic neutrinos compatible with dark energy radiation in a neutrino capture experiment on tritium. We also find that dark matter axion experiments lack sensitivity to a relativistic thermal axion background, even if enhanced by dark energy radiation, and dedicated search strategies are required to probe new parameter space. We derive constraints arising from a dark photon background from oscillations into visible photons, and find that several orders of magnitude of viable parameter space can be explored with planned experimental programs such as DM Radio and LADERA.

Michal Stratený, Georgios Lukes-Gerakopoulos

11 pages, 3 figures, 1 table, Proceedings of RAGtime 23-25. Edited by Z. Stuchl\'ik, G. T\"or\"ok and V. Karas. Institute of Physics in Opava

This work studies the dynamics of geodesic motion within a curved spacetime around a Schwarzschild black hole, perturbed by a gravitational field of a far axisymmetric distribution of mass enclosing the system. This spacetime can serve as a versatile model for a diverse range of astrophysical scenarios and, in particular, for extreme mass ratio inspirals as in our work. We show that the system is non-integrable by employing Poincar\'e surface of section and rotation numbers. By utilising the rotation numbers, the widths of resonances are calculated, which are then used in establishing the relation between the underlying perturbation parameter driving the system from integrability and the quadrupole parameter characterising the perturbed metric. This relation allows us to estimate the phase shift caused by the resonance during an inspiral.

Joaquim M. Gomes, Edward Hardy, Susha Parameswaran

30 pages + appendices (two columns), 9 figures

We analyse theories that do not have a de Sitter vacuum and cannot lead to slow-roll quintessence, but which nevertheless support a transient era of accelerated cosmological expansion due to interactions between a scalar $\phi$ and either a hidden sector thermal bath, which evolves as Dark Radiation, or an extremely-light component of Dark Matter. We show that simple models can explain the present-day Dark Energy of the Universe consistently with current observations. This is possible both when $\phi$'s potential has a hilltop form and when it has a steep exponential run-away, as might naturally arise from string theory. We also discuss a related theory of multi-field quintessence, in which $\phi$ is coupled to a sector that sources a subdominant component of Dark Energy, which overcomes many of the challenges of slow-roll quintessence.

Antonio Ferreiro, Samuel Monin, Francisco Torrenti

22 pages + references, 5 figures

We propose a new scheme to regularize the stress-energy tensor and the two-point function of free quantum scalar fields propagating in cosmological spacetimes. We generalize the adiabatic regularization method by introducing two additional mass scales not present in the standard program. If we set them to the order of the physical scale of the problem, we obtain ultraviolet-regularized quantities that do not distort the amplitude of the power spectra at the infrared momentum scales amplified by the non-adiabatic expansion of the universe. This is not ensured by the standard adiabatic method. We also show how our proposed subtraction terms can be interpreted as renormalization of coupling constants in the Einstein's equations. We finally illustrate our proposed regularization method in two scenarios of cosmological interest: de Sitter inflation and geometric reheating.

Delaney Farrell, Fridolin Weber, Jia Jie Li, Armen Sedrakian

7 pages, 6 figures

Neutron stars may experience differential rotation on short, dynamical timescales following extreme astrophysical events like binary neutron star mergers. In this work, the masses and radii of differentially rotating neutron star models are computed. We employ a set of equations of states for dense hypernuclear and $\Delta$-admixed-hypernuclear matter obtained within the framework of CDF theory in the relativistic Hartree-Fock (RHF) approximation. Results are shown for varying meson-$\Delta$ couplings, or equivalently the $\Delta$-potential in nuclear matter. A comparison of our results with those obtained for non-rotating stars shows that the maximum mass difference between differentially rotating and static stars is independent of the underlying particle composition of the star. We further find that the decrease in the radii and increase in the maximum masses of stellar models when $\Delta$-isobars are added to hyperonuclear matter (as initially observed for static and uniformly rotating stars) persist also in the case of differentially rotating neutron stars.

Sudarshan Ghonge, Joshua Brandt, J. M. Sullivan, Margaret Millhouse, Katerina Chatziioannou, James A. Clark, Tyson Littenberg, Neil Cornish, Sophie Hourihane, Laura Cadonati

22 pages, 17 figures

We investigate the impact of transient noise artifacts, or {\it glitches}, on gravitational wave inference, and the efficacy of data cleaning procedures in recovering unbiased source properties. Due to their time-frequency morphology, broadband glitches demonstrate moderate to significant biasing of posterior distributions away from true values. In contrast, narrowband glitches have negligible biasing effects owing to distinct signal and glitch morphologies. We inject simulated binary black hole signals into data containing three common glitch types from past LIGO-Virgo observing runs, and reconstruct both signal and glitch waveforms using {\tt BayesWave}, a wavelet-based Bayesian analysis. We apply the standard LIGO-Virgo-KAGRA deglitching procedure to the detector data - we subtract the glitch waveform estimated by the joint {\tt BayesWave} inference before performing parameter estimation with detailed compact binary waveform models. We find that this deglitching effectively mitigates bias from broadband glitches, with posterior peaks aligning with true values post deglitching. This provides a baseline validation of existing techniques, while demonstrating waveform reconstruction improvements to the Bayesian algorithm for robust astrophysical characterization in glitch-prone detector data.