Abstract visibility
Abstract text size

Papers for Friday, Feb 17 2023

Papers with local authors

Duncan Farrah, Kevin S. Croker, Gregory Tarlé, Valerio Faraoni, Sara Petty, Jose Afonso, Nicolas Fernandez, Kurtis A. Nishimura, Chris Pearson, Lingyu Wang, Michael Zevin, David L Clements, Andreas Efstathiou, Evanthia Hatziminaoglou, Mark Lacy, Conor McPartland, Lura K Pitchford, Nobuyuki Sakai, Joel Weiner

10 pages, 3 figures, published in ApJ Letters

0 votes
Paper 3 — arXiv:2302.07878
0 votes
Paper 3 — arXiv:2302.07878

Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of accretion or mergers, in a manner that depends on the black hole's interior solution. We test this prediction by considering the growth of supermassive black holes in elliptical galaxies over $0<z\lesssim2.5$. We find evidence for cosmologically coupled mass growth among these black holes, with zero cosmological coupling excluded at 99.98% confidence. The redshift dependence of the mass growth implies that, at $z\lesssim7$, black holes contribute an effectively constant cosmological energy density to Friedmann's equations. The continuity equation then requires that black holes contribute cosmologically as vacuum energy. We further show that black hole production from the cosmic star formation history gives the value of $\Omega_{\Lambda}$ measured by Planck while being consistent with constraints from massive compact halo objects. We thus propose that stellar remnant black holes are the astrophysical origin of dark energy, explaining the onset of accelerating expansion at $z \sim 0.7$.

Kareem El-Badry, Hans-Walter Rix, Yvette Cendes, Antonio C. Rodriguez, Charlie Conroy, Eliot Quataert, Keith Hawkins, Eleonora Zari, Melissa Hobson, Katelyn Breivik, Arne Rau, Edo Berger, Sahar Shahaf, Rhys Seeburger, Kevin B. Burdge, David W. Latham, Lars A. Buchhave, Allyson Bieryla, Dolev Bashi, Tsevi Mazeh, Simchon Faigler

22 pages, 15 figures. Submitted to MNRAS

0 votes
Paper 4 — arXiv:2302.07880
0 votes
Paper 4 — arXiv:2302.07880

We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. We show that the system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV) follow-up over a 6-month period spans most of the orbit's dynamic range in RV and is in excellent agreement with predictions of the Gaia solution. UV imaging and high-resolution optical spectra rule out all plausible luminous companions that could explain the orbit. The star is a bright ($G=12.3$), slightly metal-poor ($\rm [Fe/H]=-0.22$) low-luminosity giant ($T_{\rm eff}=4600\,\rm K$; $R = 7.9\,R_{\odot}$; $\log\left[g/\left({\rm cm\,s^{-2}}\right)\right] = 2.6$). The binary's orbit is moderately eccentric ($e=0.52$). The giant is strongly enhanced in $\alpha-$elements, with $\rm [\alpha/Fe] = +0.26$, but the system's Galactocentric orbit is typical of the thin disk. We obtained X-ray and radio nondetections of the source near periastron, which support BH accretion models in which the net accretion rate at the horizon is much lower than the Bondi-Hoyle-Lyttleton rate. At a distance of 1.16 kpc, Gaia BH2 is the second-nearest known BH, after Gaia BH1. Its orbit -- like that of Gaia BH1 -- seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve, perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 likely significantly outnumber their close, X-ray bright cousins, but their formation pathways remain uncertain.

Philippe Z. Yao, Eliot Quataert, Andy Goulding

12 pages, 8 figures, 2 tables, submitted to MNRAS

0 votes
Paper 8 — arXiv:2302.07886
0 votes
Paper 8 — arXiv:2302.07886

Many double white dwarf (WD) mergers likely do not lead to a prompt thermonuclear explosion. We investigate the prospects for observationally detecting the surviving remnants of such mergers, focusing on the case of mergers of Carbon-Oxygen WDs. For $\sim 10^4$ yr, the merger remnant is observationally similar to an extreme AGB star evolving to become a massive WD. Identifying merger remnants is thus easiest in galaxies with high stellar masses (high WD merger rate) and low star formation rates (low birth rate of $\sim 6-10 \,{\rm M_{\odot}}$ stars). Photometrically identifying merger remnants is challenging even in these cases because the merger remnants appear similar to He stars and post-outburst classical novae. We propose that the most promising technique for discovering WD merger remnants is through their unusual surrounding photoionized nebulae. We use CLOUDY photoionization calculations to investigate their unique spectral features. Merger remnants should produce weak hydrogen lines and strong carbon and oxygen recombination and fine-structure lines in the UV, optical and IR. With integral field spectrographs, we predict that hundreds of candidates are detectable in M87 and other nearby massive galaxies. Dust somewhat reduces the optical line emission for nebula radii $\lesssim 10^{17}$ cm, but the sources spend most of their time with larger radii when dust is less important. Our models roughly reproduce the WISE nebula surrounding the Galactic WD merger candidate IRAS 00500+6713; we predict detectable [Ne VI] and [Mg VII] lines with JWST but that the mid-IR WISE emission is dominated by dust not fine-structure lines.

Yan Huang, Zhuo Li

14 pages, 9 figures, 3 tables, submitted to MNRAS, comments welcome

0 votes
Paper 42 — arXiv:2302.08230
0 votes
Paper 42 — arXiv:2302.08230

The multiwavelength observations of GRB afterglows, together with some high-performance particle-in-cell simulations, hint that the magnetic field may decay behind the shock front. In this work, we develop a numerical code to calculate the evolution of the accelerated electron distribution, their synchrotron and inverse-Compton (IC) spectra and accordingly the light curves (LCs) under the assumption of decaying microturbulence (DM) downstream of the shock, $\epsilon_B(t_p')\propto t_p'^{\alpha_t}$ with $t_p'$ the fluid proper time since injection. We find: (1) The synchrotron spectrum in the DM model is similar to that in the homogeneous turbulence (HT) model with very low magnetic field strength. However, the difference in the IC spectral component is relatively more obvious between them, due to the significant change of the postshock electron energy distribution with DM. (2) If the magnetic field decay faster, there are less electrons cool fast, and the IC spectral component becomes weaker. (3) The LCs in the DM model decay steeper than in the HT model, and the spectral evolution and the LCs in the DM model is similar to the HT model where the magnetic field energy fraction decreases with observer time, $\epsilon_B(t) \propto t^{5\alpha_t /8}$. (4) The DM model can naturally produce a significant IC spectral component in TeV energy range, but due to the Klein-Nishina suppression the IC power cannot be far larger than the synchrotron power. We apply the DM model to describe the afterglow data of GRB 190114C and find the magnetic field decay exponent $\alpha_t\sim -0.4$ and the electron spectral index $p\sim2.4$. Future TeV observations of the IC emission from GRB afterglows will further help to probe the poorly known microphysics of relativistic shocks.

All other papers

Gregory W. Henry, Jacob L. Bean

arXiv only posting, contact Greg Henry for access to the data

GJ1214b is the highest signal-to-noise sub-Neptune for atmospheric studies. Although most previous transmission spectroscopy measurements have revealed a frustratingly featureless spectrum, JWST observations are expected to give new insights to this benchmark planet. We have performed photometric monitoring of GJ1214 (the host star) to provide context for these observations. We find that GJ1214 entered a period of relatively high brightness during 2021 and 2022. This implies that the JWST MIRI/LRS phase curve observation of GJ1214b in July 2022 was obtained during an epoch of low activity for the spot-dominated host star. Like previous works, we are unable to definitively identify the star's rotation period. Nevertheless, we confirm that it is likely >50 days.

A. Addis, A. Bulgarelli, N. Parmiggiani, J. Rodi, A. Bazzano

4 pages, 4 figures, Astronomical Data Analysis Software and System XXXII (31 October-4 November 2022)

In the multi-messenger era, facilities share their results with the scientific community through networks such as the General Coordinates Network to study transient phenomena (e.g., Gamma-ray bursts) and implement real-time analysis pipelines to detect transient events, reacting to science alerts received from other observatories. The fast analysis of transient events is crucial for detecting counterparts of gravitational waves and neutrino candidate events. In this context, collecting scientific results from different high-energy satellites observing the same transient event represents a key step in improving the statistical significance of the high-energy candidate events. This project aims to develop a system and a web platform to share information and scientific results of transient events between high-energy satellites with INAF participation (AGILE, FERMI, INTEGRAL and SWIFT). The AFISS platform implements the COMET VO- Event broker and provides a web portal where the users visualize the list of transient events detected by multi-messenger facilities and received through the GCN. The web portal could show, for each event, a summary of the scientific results shared by the real-time analysis pipelines and a list of time-correlated transient events. In addition, the platform is ready to receive results from participating facilities on sub-threshold events (STE) that cannot be shared with the community due to the low statistical significance. If the platform finds a time correlation between two or more STEs, it can promote them to science alerts. The web interface shows the list of STEs with possible time correlation with other STEs or science alerts. The platform notifies the users with an email when a new transient event is received.

Andrea Botteon, Maxim Markevitch, Reinout J. van Weeren, Gianfranco Brunetti, Timothy W. Shimwell

10 pages, 5 figures, 1 table (excluding Appendixes). Abstract abridged to meet arXiv requirements. Submitted to A&A

Dynamical motions in the ICM can imprint distinctive features on the X-ray images that map the thermal emission from clusters, such as sharp surface brightness discontinuities due to shocks and cold fronts. The gas dynamics during cluster mergers may also drive large-scale turbulence in the ICM which in turn generates extended synchrontron sources known as radio halos. The presence of surface brightness edges in the thermal gas of clusters has been established by a number of X-ray observations. In contrast, edges in radio halos have been observed only in a handful of cases. Our goal is to search for new radio surface brightness discontinuities in the ICM. We inspected the images of the Bullet Cluster and the other 25 radio halos reported in the MeerKAT Galaxy Cluster Legacy Survey. To aid the identification of surface brightness discontinuities, we applied a gradient filtering edge detection method to the radio images. We found that the adopted filtering technique is helpful to identify surface brightness edges in radio images, allowing us to identify at least one gradient in half of the radio halos studied. For the Bullet Cluster, we found excellent agreement between the locations of the 4 radio discontinuities detected and X-ray edges. This similarity informs us that there is substantial interplay between thermal and non-thermal components in galaxy clusters. This interplay is likely due to the forzen-in ICM magnetic field which mediates the advection of cosmic rays while being dragged by thermal gas flows. We conclude that radio halos are shaped by dynamical motions in the ICM and that they often display surface brightness discontinuities apparently co-located with edges in the thermal gas emission. Our results demonstrate that new and future generations of radio telescopes will provide a complementary approach to X-rays to efficiently detect shocks and cold fronts in the ICM.

Moritz S. Fischer, Nils-Henrik Durke, Katharina Hollingshausen, Claudius Hammer, Marcus Brüggen, Klaus Dolag

14 pages, 10 figures + appendices, submitted to MNRAS

Mergers of galaxy clusters are promising probes of dark matter (DM) physics. For example, an offset between the DM component and the galaxy distribution can constrain DM self-interactions. We investigate the role of the intracluster medium (ICM) and its influence on DM-galaxy offsets in self-interacting dark matter (SIDM) models. To this end, we employ Smoothed Particle Hydrodynamics + N-body simulations to study idealised setups of equal- and unequal-mass mergers with head-on collisions. Our simulations show that the ICM hardly affects the offsets arising shortly after the first pericentre passage compared to DM-only (DMO) simulations. But later on, e.g. at the first apocentre, the offsets can be amplified by the presence of the ICM. Furthermore, we find that cross-sections small enough not to be excluded by measurements of the core sizes of relaxed galaxy clusters have a chance to produce observable offsets. We found that different DM models affect the DM distribution and also the galaxy and and ICM distribution, including its temperature. Potentially, the position of the shock fronts, combined with the brightest cluster galaxies (BCGs), provides further clues to the properties of DM. Overall our results demonstrate that mergers of galaxy clusters at stages about the first apocentre passage could be more interesting in terms of DM physics than those shortly after the first pericentre passage. This may motivate further studies of mergers at later evolutionary stages.

Emily R. Liepold, Chung-Pei Ma, Jonelle L. Walsh

Accepted for publication in ApJL. 15 pages, 8 figures

The three-dimensional intrinsic shape of a galaxy and the mass of the central supermassive black hole provide key insight into the galaxy's growth history over cosmic time. Standard assumptions of a spherical or axisymmetric shape can be simplistic and can bias the black hole mass inferred from the motions of stars within a galaxy. Here we present spatially-resolved stellar kinematics of M87 over a two-dimensional $250\mbox{$^{\prime\prime}$} \times 300\mbox{$^{\prime\prime}$}$ contiguous field covering a radial range of 50 pc to 12 kpc from integral-field spectroscopic observations at the Keck II Telescope. From about 5 kpc and outward, we detect a prominent 25 $\mathrm{km~s}^{-1}$ rotational pattern, in which the kinematic axis (connecting the maximal receding and approaching velocities) is $40^\circ$ misaligned with the photometric major axis of M87. The rotational amplitude and misalignment angle both decrease in the inner 5 kpc. Such misaligned and twisted velocity fields are a hallmark of triaxiality, indicating that M87 is not an axisymmetrically shaped galaxy. Triaxial Schwarzschild orbit modeling with more than 4000 observational constraints enabled us to determine simultaneously the shape and mass parameters. The models incorporate a radially declining profile for the stellar mass-to-light ratio suggested by stellar population studies. We find that M87 is strongly triaxial, with ratios of $p=0.845$ for the middle-to-long principal axes and $q=0.722$ for the short-to-long principal axes, and determine the black hole mass to be $(5.37^{+0.37}_{-0.25}\pm 0.22)\times 10^9 M_\odot$, where the second error indicates the systematic uncertainty associated with the distance to M87.

Luca Tortorelli, Amata Mercurio

Accepted for publication in Frontiers of Astronomy and Space Sciences, 27 pages, 8 figures

In today's modern wide-field galaxy surveys, there is the necessity for parametric surface brightness decomposition codes characterised by accuracy, small degree of user intervention, and high degree of parallelisation. We try to address this necessity by introducing $\mathrm{morphofit}$, a highly parallelisable $\mathrm{Python}$ package for the estimate of galaxy structural parameters. The package makes use of wide-spread and reliable codes, namely $\mathrm{SExtractor}$ and $\mathrm{GALFIT}$. It has been optimised and tested in both low-density and crowded environments, where blending and diffuse light makes the structural parameters estimate particularly challenging. $\mathrm{morphofit}$ allows the user to fit multiple surface brightness components to each individual galaxy, among those currently implemented in the code. Using simulated images of single S\'ersic and bulge plus disk galaxy light profiles with different bulge-to-total luminosity ($\mathrm{B/T}$) ratios, we show that $\mathrm{morphofit}$ is able to recover the input structural parameters of the simulated galaxies with good accuracy. We also compare its estimates against existing literature studies, finding consistency within the errors. We use the package in Tortorelli et al. 2023 to measure the structural parameters of cluster galaxies in order to study the wavelength dependence of the Kormendy relation of early-type galaxies. The package is available on github (https://github.com/torluca/morphofit) and on the Pypi server (https://pypi.org/project/morphofit/).

D. B. Malesani, A. J. Levan, L. Izzo, A. de Ugarte Postigo, G. Ghirlanda, K. E. Heintz, D. A. Kann, G. P. Lamb, J. Palmerio, O. S. Salafia, R. Salvaterra, N. R. Tanvir, J. F. Agüí Fernández, S. Campana, A. A. Chrimes, P. D'Avanzo, V. D'Elia, M. Della Valle, M. De Pasquale, J. P. U. Fynbo, N. Gaspari, B. P. Gompertz, D. H. Hartmann, J. Hjorth, P. Jakobsson, E. Palazzi, E. Pian, G. Pugliese, M. E. Ravasio, A. Rossi, A. Saccardi, P. Schady, B. Schneider, J. Sollerman, R. L. C. Starling, C. C. Thöne, A. J. van der Horst, S. D. Vergani, D. Watson, K. Wiersema, D. Xu, T. Zafar

9 pages, 4 figures, submitted to Astronomy & Astrophysics

Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium.

L. Tortorelli, A. Mercurio, G. Granata, P. Rosati, C. Grillo, M. Nonino, A. Acebron, G. Angora, P. Bergamini, G. B. Caminha, U. Meštrić, E. Vanzella

Submitted to Astronomy and Astrophysics in the form of letter to the Editor, 5 pages, 1 figure, 1 table

The wavelength dependence of the Kormendy relation (KR) is well characterised at low-redshift, but poorly studied at intermediate redshifts. The KR provides information on the evolution of the population of early-type galaxies (ETGs), therefore, by studying it, we may shed light on the assembly processes of these objects and their size evolution. Since studies at different redshifts are generally conducted in different rest-frame wavebands, investigating whether there is a wavelength dependence of the KR is fundamental to interpret the conclusions we might draw from it. We analyse the KRs of the three Hubble Frontier Fields clusters, Abell S1063 (z = 0.348), MACS J0416.1-2403 (z = 0.396), and MACS J1149.5+2223 (z = 0.542), as a function of wavelength. This is the first time the KR of ETGs has been explored consistently in such a large range of wavelength at intermediate redshifts. We exploit very deep HST photometry, ranging from the observed B-band to the H-band, and VLT/MUSE integral field spectroscopy. We improve the structural parameters estimation we performed in a previous work (Tortorelli et al. 2018) by means of a newly developed Python package called morphofit (Tortorelli&Mercurio 2023). With its use on cluster ETGs, we find that the KR slopes smoothly increase with wavelength from the optical to the near-infrared bands in all three clusters, with the intercepts getting fainter at lower redshifts due to the passivisation of the ETGs population. The slope trend is consistent with previous findings at lower redshifts. The slope increase with wavelength implies that smaller size ETGs are more centrally concentrated than larger size ETGs in the near-infrared with respect to the optical regime. Since different bands probe different stellar populations in galaxies, the slope increase also implies that smaller ETGs have stronger internal gradients with respect to larger ETGs.

Cimatti A., Moresco M

10 pages, 4 figures, 2 tables (including appendix). Submitted to ApJL

Despite the tremendous advance of observational cosmology, the value of the Hubble constant ($H_0$) is still controversial (the so called ``Hubble tension'') because of the inconsistency between local/late-time measurements and those derived from the cosmic microwave background. As the age of the Universe is very sensitive to $H_0$, we explored whether the present-day oldest stars could place independent constraints on the Hubble constant. To this purpose, we selected from the literature the oldest objects (globular clusters, stars, white dwarfs, ultra-faint and dwarf spheroidal galaxies) with accurate age estimates. Adopting a conservative prior on their formation redshifts ($11 \leq z_{\rm f} \leq 30$) and assuming $\Omega_{\rm M} = 0.3 \pm 0.02$, we developed a method based on Bayesian statistics to estimate the Hubble constant. We selected the oldest objects ($>13.3$ Gyr), and estimated $H_0$ both for each of them individually and for the average ages of homogeneous subsamples. Statistical and systematic uncertainties were properly taken into account. The constraints based on individual ages indicate that $H_0<70.6$ km/s/Mpc when selecting the most accurate estimates. If the ages are averaged and analyzed independently for each subsample, the most stringent constraints imply $H_0<73.0$ with a probability of 93.2\% and errors around 2.5 km/s/Mpc. We also constructed an ``accuracy matrix'' to assess how the constraints on $H_0$ become more stringent with further improvements in the accuracy of stellar ages and $\Omega_{\rm M}$. The results show the high potential of the oldest stars as independent and competitive cosmological probes.

Renu Malhotra, Zherui Chen

11 pages, many figures. In press at MNRAS

Mean motion resonances are important in the analysis and understanding of the dynamics of planetary systems. While perturbative approaches have been dominant in many previous studies, recent non-perturbative approaches have revealed novel properties in the low eccentricity regime for interior mean motion resonances of Jupiter in the fundamental model of the circular planar restricted three body model. Here we extend the non-perturbative investigation to exterior mean motion resonances in the low eccentricity regime (up to about 0.1) and for perturber mass in the range 5e-5 to 1e-3 (in units of the central mass). Our results demonstrate that first order exterior resonances have two branches at low eccentricity as well as low-eccentricity bridges connecting neighboring first order resonances. With increasing perturber mass, higher order resonances dissolve into chaos whereas low order resonances persist with larger widths in their radial extent but smaller azimuthal widths. For low order resonances, we also detect secondary resonances arising from small integer commensurabilities between resonant librations and the synodic frequency. These secondary resonances contribute significantly to generating the chaotic sea that typically occurs near mean motion resonances of higher-mass perturbers.

We perform a Bayesian search in the latest Pulsar Timing Array (PTA) datasets for a stochastic gravitational wave (GW) background sourced by curvature perturbations at scales $10^5~\text{Mpc}^{-1}\lesssim k\lesssim 10^8~\text{Mpc}^{-1}$. These re-enter the Hubble horizon at temperatures around and below the QCD crossover phase transition in the early Universe. We include a stochastic background of astrophysical origin in our search and properly account for constraints on the curvature power spectrum from the overproduction of primordial black holes (PBHs). We find that the International PTA Data Release 2 significantly favors the astrophysical model for its reported common-spectrum process, over the curvature-induced background. On the other hand, the two interpretations fit the NANOgrav 12.5 years dataset equally well. We then set new upper limits on the amplitude of the curvature power spectrum at small scales. These are independent from, and competitive with, indirect astrophysical bounds from the abundance of PBH dark matter. Upcoming PTA data releases will provide the strongest probe of the curvature power spectrum around the QCD epoch.

The equation of state for hydrogen and helium is fundamental for studying stars and giant planets. It has been shown that because of interactions at atomic and molecular levels, the behaviour of a mixture of hydrogen and helium cannot be accurately represented by considering these elements separately. This paper aims at providing a simple method to account for interactions between hydrogen and helium in interior and evolution models of giant planets. Using on the one hand ab initio simulations that involve a system of interacting hydrogen and helium particles and pure equations of state for hydrogen and helium on the other, we derived the contributions in density and entropy of the interactions between hydrogen and helium particles. We show that relative variations of up to 15% in density and entropy arise when non-ideal mixing is accounted for. These non-ideal mixing effects must be considered in interior models of giant planets based on accurate gravity field measurements, particularly in the context of variations in the helium-to-hydrogen ratio. They also affect the mass-radius relation of exoplanets. We provide a table that contains the volume and entropy of mixing as a function of pressure and temperature. This table is to be combined with pure hydrogen and pure helium equations of state to obtain an equation of state that self-consistently includes mixing effects for any hydrogen and helium mixing ratio and may be used to model the interior structure and evolution of giant planets to brown dwarfs. Non-linear mixing must be included in accurate calculations of the equations of state of hydrogen and helium. Uncertainties on the equation of state still exist, however. Ab initio calculations of the behaviour of the hydrogen-helium mixture in the megabar regime for various compositions should be performed in order to gain accuracy.

S. Agarwal, B. Banerjee, A. Shukla, J. Roy, S. Acharya, B. Vaidya, V. R. Chitnis, S. M. Wagner, K. Mannheim, M. Branchesi

6 pages, 3 figures, 1 table, Accepted for publication in MNRAS-L

The evolution of the spectral energy distribution during flares constrains models of particle acceleration in blazar jets. The archetypical blazar BL Lac provided a unique opportunity to study spectral variations during an extended strong flaring episode from 2020-2021. During its brightest $\gamma$-ray state, the observed flux (0.1-300 GeV) reached up to $2.15\,\times\,10^{-5}\,\rm{ph\,cm^{-2}\,s^{-1}}$, with sub-hour scale variability. The synchrotron hump extended into the X-ray regime showing a minute-scale flare with an associated peak shift of inverse-Compton hump in gamma-rays. In shock acceleration models, a high Doppler factor value $>$100 is required to explain the observed rapid variability, change of state, and $\gamma$-ray peak shift. Assuming particle acceleration in mini-jets produced by magnetic reconnection during flares, on the other hand, alleviates the constraint on required bulk Doppler factor. In such jet-in-jet models, observed spectral shift to higher energies (towards TeV regime) and simultaneous rapid variability arises from the accidental alignment of a magnetic plasmoid with the direction of the line of sight. We infer a magnetic field of $\sim0.6\,\rm{G}$ in a reconnection region located at the edge of BLR ($\sim0.02\,\rm{pc}$). The scenario is further supported by log-normal flux distribution arising from merging of plasmoids in reconnection region.

B. O'Connor, E. Troja, G. Ryan, P. Beniamini, H. van Eerten, J. Granot, S. Dichiara, R. Ricci, V. Lipunov, J. H. Gillanders, R. Gill, M. Moss, S. Anand, I. Andreoni, R. L. Becerra, D. A. H. Buckley, N. R. Butler, S. B. Cenko, A. Chasovnikov, J. Durbak, C. Francile, E. Hammerstein, A. J. van der Horst, M. Kasliwal, C. Kouveliotou, A. S. Kutyrev, W. H. Lee, G. Srinivasaragavan, V. Topolev, A. M. Watson, Y. H. Yang, K. Zhirkov

Submitted version. 53 pages, 9 figures, 6 tables

Long duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Due to its enormous energy ($E_\textrm{iso}\!\approx$10$^{55}$ erg) and proximity ($z\!\approx$0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multi-wavelength observations covering the first three months of its afterglow evolution. The X-ray brightness decays as a power-law with slope $\approx\!t^{-1.66}$, which is not consistent with standard predictions for jetted emission. We attribute this behavior to a shallow energy profile of the relativistic jet. A similar trend is observed in other energetic GRBs, suggesting that the most extreme explosions may be powered by structured jets launched by a common central engine.

M. Lafarga, I. Ribas, M. Zechmeister, A. Reiners, Á. López-Gallifa, D. Montes, A. Quirrenbach, P. J. Amado, J. A. Caballero, M. Azzaro, V. J. S. Béjar, A. P. Hatzes, Th. Henning, S. V. Jeffers, A. Kaminski, M. Kürster, P. Schöfer, A. Schweitzer, H. M. Tabernero, M. R. Zapatero Osorio

Accepted for publication in A&A

Radial velocities (RVs) measured from high-resolution stellar spectra are routinely used to detect and characterise orbiting exoplanet companions. The different lines present in stellar spectra are created by several species, which are non-uniformly affected by stellar variability features such as spots or faculae. Stellar variability distorts the shape of the spectral absorption lines from which precise RVs are measured, posing one of the main problems in the study of exoplanets. In this work we aim to study how the spectral lines present in M dwarfs are independently impacted by stellar activity. We used CARMENES optical spectra of six active early- and mid-type M dwarfs to compute line-by-line RVs and study their correlation with several well-studied proxies of stellar activity. We are able to classify spectral lines based on their sensitivity to activity in five M dwarfs displaying high levels of stellar activity. We further used this line classification to compute RVs with activity-sensitive lines and less sensitive lines, enhancing or mitigating stellar activity effects in the RV time series. For specific sets of the least activity-sensitive lines, the RV scatter decreases by ~ 2 to 5 times the initial one, depending on the star. Finally, we compare these lines in the different stars analysed, finding the sensitivity to activity to vary from star to star. Despite the high density of lines and blends present in M dwarf stellar spectra, we find that a line-by-line approach is able to deliver precise RVs. Line-by-line RVs are also sensitive to stellar activity effects, and they allow for an accurate selection of activity-insensitive lines to mitigate activity effects in RV. However, we find stellar activity effects to vary in the same insensitive lines from star to star.

Jean-Pierre Lasota

Updated version of Chapter 3 in "Active Galactic Nuclei"; F. Combes ed.; iSTE/Wiley 2022; DOI:10.1002/9781394163724. One figure (Fig. 2.) and several references added. A graduate-student level lecture on accretion discs in active galactic nuclei, not a review article. arXiv admin note: substantial text overlap with arXiv:1505.02172

I present the elements of accretion-disc physics applied to active galactic nuclei. The accretion driving mechanisms are discussed, then models of geometrically-thin discs, both stationary and time-dependent, are addressed. Disc's self-gravitation in the AGN context is presented. The shapes of spectral line from accretion discs are explained in both newtonian and relativistic cases. The physics of disc's thermal and viscous instabilities is dicussed in detail and models are applied to AGN discs. Finally, the problems of thick accretion flows (ADAFs and slim discs), disc coronae, winds and jets are discussed in some detail. This is a graduate-student level lecture, not a review article

R. Alfaro, C. Alvarez, J.C. Arteaga-Velázquez, D. Avila Rojas, H.A. Ayala Solares, R. Babu, E. Belmont-Moreno, K.S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, O. Chaparro-Amaro, U. Cotti, J. Cotzomi, E. De la Fuente, R. Diaz Hernandez, B.L. Dingus, M.A. DuVernois, M. Durocher, J.C. Díaz-Vélez, C. Espinoza, K.L. Fan, N. Fraija, J.A. García-González, F. Garfias, M.M. González, J.P. Harding, S. Hernández-Cadena, D. Huang, F. Hueyotl-Zahuantitla, A. Iriarte, V. Joshi, S. Kaufmann, D. Kieda, J. Lee, H. León Vargas, J.T. Linnemann, A.L. Longinotti, G. Luis-Raya, K. Malone, O. Martinez, J. Martínez-Castro, J.A. Matthews, E. Moreno, M. Mostafá, A. Nayerhoda, L. Nellen, N. Omodei, Y. Pérez Araujo, E.G. Pérez-Pérez, et al. (22 additional authors not shown)

22 pages, 11 figures, 3 tables

We present the results of dark matter (DM) searches in a sample of 31 dwarf irregular (dIrr) galaxies within the field of view of the HAWC Observatory. dIrr galaxies are DM dominated objects, which astrophysical gamma-ray emission is estimated to be negligible with respect to the secondary gamma-ray flux expected by annihilation or decay of Weakly Interacting Massive Particles (WIMPs). While we do not see any statistically significant DM signal in dIrr galaxies, we present the exclusion limits ($95\%~\text{C.L.}$) for annihilation cross-section and decay lifetime for WIMP candidates with masses between $1$ and $100~\text{TeV}$. Exclusion limits from dIrr galaxies are relevant and complementary to benchmark dwarf Spheroidal (dSph) galaxies. In fact, dIrr galaxies are targets kinematically different from benchmark dSph, preserving the footprints of different evolution histories. We compare the limits from dIrr galaxies to those from ultrafaint and classical dSph galaxies previously observed with HAWC. We find that the contraints are comparable to the limits from classical dSph galaxies and $\thicksim2$ orders of magnitude weaker than the ultrafaint dSph limits.

Kevin Almeida Cheminant (1), Dariusz Góra (1), Nataliia Borodai (1), Ralph Engel (2), Tanguy Pierog (2), Jan Pękala (1), Markus Roth (2), Jarosław Stasielak (1), Michael Unger (2), Darko Veberič (2), Henryk Wilczyński (1) ((1) Institute of Nuclear Physics PAS, (2) Karlsruhe Institute of Technology Institute for Astroparticle Physics)

8 pages, 5 figures, 2 tables, accepted for publication in the proceedings of the 27th European Cosmic Ray Symposium

Simulations of extensive air showers using current hadronic interaction models predict too small numbers of muons compared to events observed in the air-shower experiments, which is known as the muon-deficit problem. In this work, we present a new method to calculate the factor by which the muon signal obtained via Monte-Carlo simulations must be rescaled to match the data, as well as the beta exponent from the Heitler-Matthews model which governs the number of muons found in an extensive air shower as a function of the mass and the energy of the primary cosmic ray. This method uses the so-called z variable (difference between the total reconstructed and the simulated signals), which is connected to the muon signal and is roughly independent of the zenith angle, but depends on the mass of the primary cosmic ray. Using a mock dataset built from QGSJetII-04, we show that such a method allows us to reproduce the average muon signal from this dataset using Monte-Carlo events generated with the EPOS-LHC hadronic model, with accuracy better than 6%. As a consequence of the good recovery of the muon signal for each primary included in the analysis, also the beta exponent can be obtained with accuracy of less than 1% for the studied system. Detailed simulations show a dependence of the beta exponent on hadronic interaction properties, thus the determination of this parameter is important for understanding the muon deficit problem.

Juan S. Cruz, Steen Hannestad, Emil Brinch Holm, Florian Niedermann, Martin S. Sloth, Thomas Tram

14 pages, 7 figures, 1 table

Recent interest in New Early Dark Energy (NEDE), a cosmological model with a vacuum energy component decaying in a triggered phase transition around recombination, has been sparked by its impact on the Hubble tension. Previous constraints on the model parameters were derived in a Bayesian framework with Markov-chain Monte Carlo (MCMC) methods. In this work, we instead perform a frequentist analysis using the profile likelihood in order to assess the impact of prior volume effects on the constraints. We constrain the maximal fraction of NEDE $f_\mathrm{NEDE}$, finding $f_\mathrm{NEDE}=0.076^{+0.040}_{-0.035}$ at $68 \%$ CL with our baseline dataset and similar constraints using either data from SPT-3G, ACT or full-shape large-scale structure, showing a preference over $\Lambda$CDM even in the absence of a SH0ES prior on $H_0$. While this is stronger evidence for NEDE than obtained with the corresponding Bayesian analysis, our constraints broadly match those obtained by fixing the NEDE trigger mass. Including the SH0ES prior on $H_0$, we obtain $f_\mathrm{NEDE}=0.136^{+0.024}_{-0.026}$ at $68 \%$ CL. Furthermore, we compare NEDE with the Early Dark Energy (EDE) model, finding similar constraints on the maximal energy density fractions and $H_0$ in the two models. At $68 \%$ CL in the NEDE model, we find $H_0 = 69.56^{+1.16}_{-1.29} \text{ km s}^{-1}\text{ Mpc}^{-1}$ with our baseline and $H_0 = 71.62^{+0.78}_{-0.76} \text{ km s}^{-1}\text{ Mpc}^{-1}$ when including the SH0ES measurement of $H_0$, thus corroborating previous conclusions that the NEDE model provides a considerable alleviation of the $H_0$ tension.

Edoardo Altamura, Scott T. Kay, Jens Chluba, Imogen Towler

Submitted to Monthly Notices of the Royal Astronomical Society. Comments and discussions are welcome. Data and codes can be found at this https URL

The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulations. To utilise the rotational kSZ as a cosmological probe, simulations offer some of the most comprehensive data sets that can inform the modeling of this signal. In this work, we use the MACSIS data set to specifically investigate the rotational kSZ effect in massive clusters. Based on these models, we test stacking approaches and estimate the amplitude of the combined signal with varying mass, dynamical state, redshift and map-alignment geometry. We find that the dark matter, galaxy and gas spins are generally misaligned, an effect that can cause a sub-optimal estimation of the rotational kSZ effect when based on galaxy catalogues. Furthermore, we provide halo-spin-mass scaling relations that can be used to build a statistical model of the rotational kSZ. The rotational kSZ contribution, which is largest in massive unrelaxed clusters ($\gtrsim$100 $\mu$K), could be relevant to studies of higher-order CMB temperature signals, such as the moving lens effect. The limited mass range of the MACSIS sample strongly motivates an extended investigation of the rotational kSZ effect in large-volume simulations to refine the modelling, particularly towards lower mass and higher redshift, and provide forecasts for upcoming cosmological CMB experiments (e.g. Simons Observatory, SKA-2) and X-ray observations (e.g. \textit{Athena}/X-IFU).

Niall Whiteford, Alistair Glasse, Katy L. Chubb, Daniel Kitzmann, Shrishmoy Ray, Mark W. Phillips, Beth A. Biller, Paul I. Palmer, Ken Rice, Ingo P. Waldmann, Quentin Changeat, Nour Skaf, Jason Wang, Billy Edwards, Ahmed Al-Refaie

Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$\sigma$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures.

Javier Licandro, John Tonry, Miguel R. Alarcon, Miquel Serra-Ricart, Larry Denneau

5 pages, 8 figures, Conference paper

In this work we present the design of the ATLAS unit (Asteroid Terrestrial-impact Last Alert System) that will be installed at Teide Observatory in Tenerife island (Spain). ATLAS-Teide will be built by the Instituto de Astrofisica de Canarias (IAC) and will be operated as part of the ATLAS network in the framework of an operation and science exploitation agreement between the IAC and the ATLAS team at University of Hawaii. ATLAS-Teide will be the first ATLAS unit based on commercial on the shelf (COTS) components. Its design is modular, each module (building block) consist of four Celestron RASA 11 telescopes that point to the same sky field, equipped with QHY600PRO CMOS cameras on an equatorial Direct Drive mount. Each module is equivalent to a 56cm effective diameter telescope and provides a 7.3 deg^2 field of view and a 1.26 arcsec/pix plate scale. ATLAS-Teide will consist of four ATLAS modules in a roll-off roof building. This configuration allows to cover the same sky area of the actual ATLAS telescopes. The first ATLAS module was installed in November 2022 in an existing clamshell at the TO. This module (ATLAS-P) is being used as a prototype to test the system capabilities, develop the needed software (control, image processing, etc.) and complete the fully integration of ATLAS-Teide in the ATLAS network. The preliminary results of the tests are presented here, and the benefits of the new ATLAS design are discussed.

Benne W. Holwerda, Frank Bigiel, Albert Bosma, Helene M. Courtois, Nathan Deg, Helga Dénes, Ahmed Elagali, Bi-Qing For, Baerbel Koribalski, Denis A. Leahy, Karen Lee-Waddell, Ángel R. López-Sánchez, Se-Heon Oh, Tristan N. Reynolds, Jonghwan Rhee, Kristine Spekkens, Jing Wang, Tobias Westmeier, O. Ivy Wong

17 figures, 12 figures, 7 tables, accepted by MNRAS

Galaxy morphology in atomic hydrogen (HI) and in the ultra-violet (UV) are closely linked. This has motivated their combined use to quantify morphology over the full H i disk for both H i and UV imaging. We apply galaxy morphometrics: Concentration, Asymmetry, Gini, M20 and Multimode-Intensity-Deviation statistics to the first moment-0 maps of the WALLABY survey of galaxies in the Hydra cluster center. Taking advantage of this new HI survey, we apply the same morphometrics over the full HI extent on archival GALEX FUV and NUV data to explore how well HI truncated, extended ultraviolet disk (XUV) and other morphological phenomena can be captured using pipeline WALLABY data products. Extended HI and UV disks can be identified relatively straightforward from their respective concentration. Combined with WALLABY HI, even the shallowest GALEX data is sufficient to identify XUV disks. Our second goal is to isolate galaxies undergoing ram-pressure stripping in the H i morphometric space. We employ four different machine learning techniques, a decision tree, a k-nearest neighbour, a support-vector machine, and a random forest. Up to 80% precision and recall are possible with the Random Forest giving the most robust results.

Pascal M. Keller, Bojan Nikolic, Nithyanandan Thyagarajan, Chris L. Carilli, Gianni Bernardi, Ntsikelelo Charles, Landman Bester, Oleg M. Smirnov, Nicholas S. Kern, Joshua S. Dillon, Bryna J. Hazelton, Miguel F. Morales, Daniel C. Jacobs, Aaron R. Parsons, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Nico Eksteen, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz, Steven R. Furlanetto, Kingsley Gale-Sides, Brian Glendenning, Deepthi Gorthi, Bradley Greig, Jasper Grobbelaar, Ziyaad Halday, Jacqueline N. Hewitt, Jack Hickish, Austin Julius, et al. (36 additional authors not shown)

16 pages, 14 figures, accepted for publication by MNRAS

Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related.

Debora Mroczek, M. Coleman Miller, Jacquelyn Noronha-Hostler, Nicolas Yunes

10 pages, 4 figures, CSQCD IX proceedings

Gaussian processes provide a promising framework by which to extrapolate the equation of state (EoS) of cold, catalyzed matter beyond $1-2$ times nuclear saturation density. Here we discuss how to extend Gaussian processes to include nontrivial features in the speed of sound, such as bumps, kinks, and plateaus, which are predicted by nuclear models with exotic degrees of freedom. Using a fully Bayesian analysis incorporating measurements from X-ray sources, gravitational wave observations, and perturbative QCD results, we show that these features are compatible with current constraints and report on how the features affect the EoS posteriors.

The extremely accurate estimates of stellar variability and radial velocity in the Gaia Data Release 3 (Gaia DR3) have enabled us to examine the close binarity and radial velocity (RV) of central stars (CSs) of planetary nebulae (PNe). This study is twofold: (1) searching for new close binary CSs candidates to better understand how binarity affects the formation and evolution of PNe; and (2) extending the sample size of known RV of PNe in order to understand their kinematics and the dynamics of the Milky Way. As a target sample, we used all true, possible, and likely PNe available in the literature. Then, we looked for their matched Gaia DR3 sources that provide measurements of variability and RV. As a result, we detected the first large collection of trustworthy photometric variability of 26 symbiotic stars (SySts) and 82 CSs. In this CS group, there are 24 sources already classified as true close binary CSs in the literature. Hence, we discovered 58 new close binary CS candidates. This close binary (CB) sample represents more than half of what is currently available in the literature. In addition, we identified the radial velocities for 51 PNe. To our knowledge, 24 of these were measured for the first time. The RV measurements predicted by Gaia, based on the Doppler shift of the CS absorption lines, and those derived from nebular emission lines, show satisfactory agreement except for a few extremely high-velocity PNe.

Longhua Qin, Jiancheng Wang, Quangui Gao, Weiwei Na, Huaizhen Li, Ao Wang, Chuyuan Yang, Jianping Yang

Accepted by MNRAS

The extragalactic background light (EBL) in the IR to UV bands partly absorbs very high energy (VHE, $E \geq$ 100GeV) $\gamma-$ray photons travelling over cosmological distances via pair production. In this paper, to get stronger constraints on EBL, we use the deliberate selection of the EBL model and data of five BL Lacs with better statistics and the harder spectra to limit the EBL density and the radiation mechanism of BL Lacs. We constrain the upper limit of the EBL density by fitting the spectral energy distributions (SEDs) of TeV BL Lacs and find that our results are compatible with the published measurement, reaching 50 $\rm{nW m^{-2} sr^{-1}}$. We also obtain that the EBL is not necessarily transparent to high VHE photons. We fix the intrinsic spectral index $\Gamma_i$ of TeV BL Lacs as 1.0 and 1.5 under observation evidence and model assumption. Comparing the EBL density given by galaxy count and $Spitzer$ observations, we then obtain that 1ES 1101-232 has $\Gamma_i$ $\leq$ 1.0 and 1ES 0229+200 should have $\Gamma_i$ not harder than 1.0. We demonstrate that the common radiation assumption of BL Lacs, in which the $\Gamma_i$ is softer than 1.5, should be revisited. Furthermore, we propose that the upper EBL density could be given by fitting the hardest energy spectra of TeV BL Lacs.

Tessa Vernstrom, Jennifer West, Franco Vazza, Denis Wittor, Christopher John Riseley, George Heald

Published in Science Advances. Supplementary Materials available at this https URL

On the largest scales, galaxies are pulled together by gravity to form clusters, which are connected by filaments making a web-like pattern. Radio emission is predicted from this cosmic web, which should originate from the strong accretion shocks around the cosmic structures. We present the first observational evidence that Fermi-type acceleration from strong shocks surrounding the filaments of the cosmic web, as well as in peripherals of low-mass clusters, is at work in the Universe. Using all-sky radio maps and stacking on clusters and filaments, we have detected the polarization signature of the synchrotron emission with polarization fractions >= 20%, which is best explained by the organization of local magnetic fields by strong shock waves both at the cluster peripheries and between clusters. Our interpretation is well supported by a detailed comparison with state-of-the-art cosmological simulations.

Rebecca Mayes, Michael Drinkwater, Joel Pfeffer, Holger Baumgardt

14 pages, Submitted to MNRAS

We use the hydrodynamic EAGLE simulation to predict the numbers and masses of supermassive black holes in stripped nuclei and compare these to confirmed measurements of black holes in observed UCDs. We find that black holes in stripped nuclei are consistent with the numbers and masses of those in observed UCDs. Approximately 50 per cent of stripped nuclei with $M > 2 \times 10^6 M_\odot$ should contain supermassive black holes. We further calculate a mass elevation ratio, $\Psi$ of the population of simulated stripped nuclei of $\Psi_{sim} = 1.51^{+0.06}_{-0.04}$ for $M > 10^7 M_\odot$ stripped nuclei, consistent with that of observed UCDs which have $\Psi_{obs} = 1.7 \pm 0.2$ above $M > 10^7 M_\odot$. We also find that the mass elevation ratios of stripped nuclei with supermassive black holes can explain the observed number of UCDs with elevated mass-to-light ratios. Finally, we predict the relative number of massive black holes in stripped nuclei and galaxy nuclei and find that stripped nuclei should increase the number of black holes in galaxy clusters by 30-100 per cent, depending on the black hole occupation fraction of low-mass galaxies. We conclude that the population of supermassive black holes in UCDs represents a large and unaccounted-for portion of supermassive black holes in galaxy clusters.

We explore the environment of 351 Grand-design and 541 Flocculent spiral galaxies recently identified employing machine learning techniques from the $17^{th}$ data release of Sloan Digital Sky Survey. We introduce a novel estimator called the Local Geometric Index to quantify the morphology of the local environment of these 892 galaxies. Based on the Local Geometric Index of the galaxies, we identify their local environments to be void, sheet, filament or cluster. We find that the Grand-designs are mostly located in dense environments like clusters and filaments, whereas the Flocculents lie in sparse environments in voids and sheets. A $p$-value $<$ $10 ^{-10}$ from a Kolmogorov-Smirnov test indicates that our results are statistically significant at $99.9\%$ confidence level. Further, we note that dense environments with large tidal flows are dominated by the Grand-designs. Metal-poor environments, such as sheets and voids, with a high abundance of gas clouds, on the other hand, are mostly populated by the Flocculents.

Bao Nguyen, Pierre Christian, Chi-kwan Chan

10 pages, 7 figures

Direct imaging of supermassive black holes at event horizon scale resolutions, as recently done by the Event Horizon Telescope (EHT), allows for testing alternative models to supermassive black holes such as Kerr naked singularities (KNSs). We demonstrate that the KNS shadow can be closed, open, or vanishing, depending on the spins and observational inclination angles. We study the critical parameters where the KNS shadow opens a gap, a distinctive phenomenon that does not happen with the black hole shadow. We show that the KNS shadow can only be closed for dimensionless spin $a \lesssim 1.18$ and vanishing for $a \gtrsim 1.18$ for certain ranges of inclination angles. We also perform numerical general relativistic ray tracing calculations, which reproduce the analytical topological change in the KNS shadow and illustrate other observational features within the shadow due to the lack of an event horizon. By comparing with black hole shadow observations, the topological change in the shadow of KNSs can be used to test the cosmic censorship hypothesis and KNSs as alternative models to SMBHs.

Amal A. Rahman, S. Sahayanathan, Malik Zahoor, P. A Subha

Submitted to MNRAS Main Joural

Detection of hard X-ray spectrum from the kilo-parsec scale jet of active galactic nuclei cannot be accounted to the synchrotron emission mechanism from the electron distribution responsible for the radio/optical emission. Alternate explanations are the inverse Compton scattering of cosmic microwave background photons (IC/CMB) or synchrotron emission from a second electron population. When the X-ray emission is interpreted as IC/CMB process, the Compton spectrum peak at GeV energy and were predicted to be the Fermi candidate sources. The non-detection of significant gamma ray flux from these galaxies by Fermi disfavoured the IC/CMB interpretation of the high energy emission. We extend this study to predict the very high energy (VHE) gamma ray emission due to IC/CMB model which can be investigated by Cherenkov Telescope Array(CTA). The model parameters deciding the broadband spectral energy distribution are estimated using analytical approximation of the emissivity functions. The emission model is extrapolated to VHE energy and then compared with the CTAO sensitivity. Particularly, we selected the sources for which the IC/CMB model is not ruled out by initial Fermi observations.

M. Fenucci, G. F. Gronchi, B. Novaković

Accepted for publication in A&A

Context. From numerical simulations, it is known that some secular resonances may affect the motion of near-Earth objects (NEOs). However, the specific location of the secular resonance inside the NEO region is not fully known, because the methods previously used to predict their location can not be used for highly eccentric orbits and when the NEOs cross the orbits of the planets. Aims. In this paper, we aim to map the secular resonances with the planets from Venus to Saturn in the NEO region, even for high values of the eccentricity. Methods. We used an averaged semi-analytical model that can deal with orbit crossing singularities for the computation of the secular dynamics of NEOs, from which we can obtain suitable proper elements and proper frequencies. Then, we computed the proper frequencies over a uniform grid in the proper elements space. Secular resonances are thus located by the level curves corresponding to the proper frequencies of the planets. Results. We determined the location of the secular resonances with the planets from Venus to Saturn, showing that they appear well inside the NEO region. By using full numerical N-body simulations we also showed that the location predicted by our method is fairly accurate. Finally, we provided some indications about possible dynamical paths inside the NEO region, due to the presence of secular resonances.

Wei Liu, Aoxiang Jiang, Wenjuan Fang

33 pages, 5 + 4 figures, 4 tables. To be submitted to JCAP. Comments welcome. arXiv admin note: text overlap with arXiv:2204.02945

The characteristic signatures of massive neutrinos on large-scale structure (LSS), if fully captured, can be used to put a stringent constraint on their mass sum, $M_{\nu}$. Previous work utilizing N-body simulations has shown the Minkowski functionals (MFs) of LSS can reveal the imprints of massive neutrinos on LSS, provide important complementary information to two-point statistics and significantly improve constraints on $M_{\nu}$. In this work, we take a step forward and apply the statistics to the biased tracers of LSS, i.e. the galaxies, and in redshift space. We perform a Fisher matrix analysis and quantify the constraining power of the MFs by using the Molino mock galaxy catalogs, which are constructed based on the halo occupation distribution (HOD) framework with parameters for the SDSS $M_r < -21.5$ and -22 galaxy samples. We find the MFs give tighter constraints on all of the cosmological parameters that we consider than the power spectrum. The constraints on $\Omega_{\mathrm{m}}, \Omega_{\mathrm{b}}, h, n_s, \sigma_8$, and $M_\nu$ from the MFs are better by a factor of 1.9, 2.9, 3.7, 4.2, 2.5, and 5.7, respectively, after marginalizing over the HOD parameters. Specifically, for $M_{\nu}$, we obtain a 1$\sigma$ constraint of 0.059 eV with the MFs alone for a volume of only $\left(1 h^{-1} \mathrm{Gpc}\right)^3$.

Anna Lia Longinotti (UNAM), Quentin Salomé (FINCA, Turku), Chiara Feruglio (INAF Trieste), Yair Krongold (UNAM), Santiago García-Burillo (OAN Madrid), Marcello Giroletti (INAF Bologna), Francesca Panessa (IAPS Roma), Carlo Stanghellini (INAF Bologna), Olga Vega (INAOE), Victor Manuel Patiño-Álvarez (INAOE), Vahram Chavushyan (INAOE), Mauricio Elías-Chavez (INAOE), Aitor Robleto-Orús (UNAM)

Accepted for publication on MNRAS, 17 pages, 7 figures

The Narrow Line Seyfert 1 Galaxy IRAS17020+4544 is one of the few AGN where a galaxy-scale energy-conserving outflow was revealed. This paper reports on NOEMA observations addressed to constrain the spatial scale of the CO emission in outflow. The molecular outflowing gas is resolved in five components tracing approaching and receding gas, all located at a distance of 2-3~kpc on the West and East side of the active nucleus. This high velocity gas (up to v_out=~1900 km/s) is not coincident with the rotation pattern of the CO gas in the host galaxy disk. The estimated mass outflow rate shows that with a global mass output of $\dot{M}_{H_2}$=~139$\pm$20$~M_\odot$~yr$^{-1}$, this powerful galaxy-scale outflow is consistent with the wind conserving its energy, and with a momentum rate boost of a factor of ~30 compared to the momentum rate of the nuclear X-ray wind. Preliminary results from ancillary X-ray (Chandra) and radio images (e-MERLIN) are reported. While the nature of the radio source is not conclusive, the Chandra image may tentatively trace extended emission, as expected by an expanding bubble of hot X-ray gas. The outcome of the NOEMA analysis and of past and ongoing publications dedicated to the description of the outflow multi-band phenomenology in IRAS17020+4544 concur to provide compelling reasons to postulate that an outflow shocking with the galaxy interstellar medium is driving the multi-phase wind in this peculiar AGN.

S. del Palacio, F. García, M. De Becker, D. Altamirano, V. Bosch-Ramon, P. Benaglia, B. Marcote, G. E. Romero

12 pages, 6 figures, accepted for publication in A&A

Context: Massive colliding-wind binaries (CWBs) can be non-thermal sources. The emission produced in their wind-collision region (WCR) encodes information of both the shocks properties and the relativistic electrons accelerated in them. The recently discovered system Apep, a unique massive system hosting two Wolf-Rayet stars, is the most powerful synchrotron radio emitter among the known CWBs, being an exciting candidate to investigate the non-thermal processes associated with stellar wind shocks. Aims: We intend to break the degeneracy between the relativistic particle population and the magnetic field strength in the WCR of Apep by probing its hard X-ray spectrum, where inverse-Compton (IC) emission is expected to dominate. Methods: We observe Apep with NuSTAR for 60 ks and combine this with a re-analysis of a deep archival XMM-Newton observation to better constrain the X-ray spectrum. We use a non-thermal emission model to derive physical parameters from the results. Results: We detect hard X-ray emission consistent with a power-law component. This is compatible with IC emission produced in the WCR for a magnetic field of 100-160 mG and a fraction of ~1.5e-4 of the total wind kinetic power being converted into relativistic electron acceleration. Conclusions: This is the first time that the non-thermal emission from a CWB is detected both in radio and high energies. This allows us to derive the most robust constraints of the particle acceleration efficiency and magnetic field intensity in a CWB so far, reducing the typical uncertainty of a few orders of magnitude to just within a factor of two. This constitutes an important step forward in our characterisation of the physical properties of CWBs.

Liana Rauf, Cullan Howlett, Tamara M. Davis, Claudia D. P. Lagos

19 pages, 16 figures - Will be submitted to MNRAS

In this work, we explore the connection between the gravitational wave (GW) merger rates of stellar-mass binary black holes (BBH) and galaxy properties. We do this by generating populations of stars using the binary synthesis code COMPAS and evolving them in galaxies from the semi-analytic galaxy formation model {\sc Shark}, to determine the number of mergers occurring in each simulation time-step. We find that large, metal-rich galaxies with high star formation rates are more likely to have gravitational wave (GW) events compared to younger, more metal poor galaxies. Our simulation with the default input parameters predicts a higher local merger rate density compared to the third gravitational wave transient catalogue (GWTC-3) prediction from LIGO, VIRGO and KAGRA. This is due to short coalescence times, low metallicities and an high formation rate of stars at low redshift in the galaxy simulation, which produces more BBHs that merge within the age of the Universe compared to observations. We identify alternate remnant mass models that more accurately reproduce the observed volumetric rate and provide updated fits to the merger rate distribution as a function of redshift. We then investigate the relative fraction of GW events in our simulation volume that are in observable host galaxies from different upcoming photometric and spectroscopic surveys, determining which of those are more ideal for tracing host galaxies with high merger rates. The implications of this work can be utilised for numerous applications, including for constraining stellar evolution models, better informing follow-up programs, and placing more informative priors on potential host galaxies when measuring cosmological parameters such as the Hubble constant.

Georgios Dimitriadis, Kate Maguire, Viraj R. Karambelkar, Ryan J. Lebron, Chang Liu, Alexandra Kozyreva, Adam A. Miller, Ryan Ridden-Harper, Joseph P. Anderson, Ting-Wan Chen, Michael Coughlin, Massimo Della Valle, Andrew Drake, Lluís Galbany, Mariusz Gromadzki, Steven L. Groom, Claudia P. Gutiérrez, Nada Ihanec, Cosimo Inserra, Joel Johansson, Tomás E. Müller-Bravo, Matt Nicholl, Abigail Polin, Ben Rusholme, Steve Schulze, Jesper Sollerman, Shubham Srivastav, Kirsty Taggart, Qinan Wang, Yi Yang, David R. Young

19 pages, 16 figures, accepted for publication in MNRAS

We present a photometric and spectroscopic analysis of the ultra-luminous and slowly evolving 03fg-like Type Ia SN 2021zny. Our observational campaign starts from $\sim5.3$ hours after explosion (making SN 2021zny one of the earliest observed members of its class), with dense multi-wavelength coverage from a variety of ground- and space-based telescopes, and is concluded with a nebular spectrum $\sim10$ months after peak brightness. SN 2021zny displayed several characteristics of its class, such as the peak brightness ($M_{B}=-19.95$ mag), the slow decline ($\Delta m_{15}(B) = 0.62$ mag), the blue early-time colours, the low ejecta velocities and the presence of significant unburned material above the photosphere. However, a flux excess for the first $\sim1.5$ days after explosion is observed in four photometric bands, making SN 2021zny the third 03fg-like event with this distinct behavior, while its $+313$ d spectrum shows prominent [O I] lines, a very unusual characteristic of thermonuclear SNe. The early flux excess can be explained as the outcome of the interaction of the ejecta with $\sim0.04\:\mathrm{M_{\odot}}$ of H/He-poor circumstellar material at a distance of $\sim10^{12}$ cm, while the low ionization state of the late-time spectrum reveals low abundances of stable iron-peak elements. All our observations are in accordance with a progenitor system of two carbon/oxygen white dwarfs that undergo a merger event, with the disrupted white dwarf ejecting carbon-rich circumstellar material prior to the primary white dwarf detonation.

Ziwei Li, Xiangkun Liu, Zuhui Fan

20 pages, 15 figures. Accepted for publication in MNRAS

In weak-lensing cosmological studies, peak statistics is sensitive to nonlinear structures and thus complementary to cosmic shear two-point correlations. In this paper, we explore a new approach, namely, the peak steepness statistics, with the overall goal to understand the cosmological information embedded there in comparison with the commonly used peak height statistics. We perform the analyses with ray-tracing simulations considering different sets of cosmological parameters $\Omega_{\rm m}$ and $\sigma_8$. A theoretical model to calculate the abundance of high peaks based on steepness is also presented, which can well describe the main trend of the peak distribution from simulations. We employ $\Delta\chi^2$ and Fisher analyses to study the cosmological dependence of the two peak statistics using our limited sets of simulations as well as our theoretical model. Within our considerations without including potential systematic effects, the results show that the steepness statistics tends to have higher sensitivities to the cosmological parameters than the peak height statistics and this advantage is diluted with the increase of the shape noise. Using the theoretical model, we investigate the physical reasons accounting for the different cosmological information embedded in the two statistics. Our analyses indicate that the projection effect from large-scale structures plays an important role to enhance the gain from the steepness statistics. The redshift and cosmology dependence of dark matter halo density profiles also contributes to the differences between the two statistics.

Aya Kubota (1), Chris Done (2), Kazuki Tsurumi (1), Ryuki Mizukawa (1) ((1) Shibaura Institute of Technology, (2) University of Durham)

29 pages, 15 figures, submitted to PASJ

Stellar-mass black hole binary systems show a strong quasi-periodic oscillation (QPO) in their Comptonised emission. The frequency of this feature correlates with the ratio of a disc to Comptonised emission. We build a new energy-conserving model of the accretion flow, SSsed, which is based on the previous agnsed with complex Comptonisation. However, it is tuned to suit for stellar black holes by including a colour temperature correction, and allowing a little more freedom in the potential to compensate for uncertainties in relativistic correction depending on the corona shape. The model is applied to hundreds of RXTE data of a black hole binary XTE J1550-564. It constrains the model parameters of the innermost radii of the accretion flow, Rin, the outer radius of a hot inner flow which extends from Rin-Rhot and the outer radius of the softer (probably nonthermal) Comptonisation region extending from Rhot-Rcor. We showed that the outer radius of the combined Comptonised emission regions, Rcor anti-correlates well with the central frequencies of low-frequency QPOs detected during the same observations. The relation is remarkably consistent with the quantitative predictions of Lense-Thirring precession of the entire Comptonisation regions for the assumed system parameters. This strongly supports the scenario that low-frequency QPOs are caused by Lense-Thirring precession.

Diego Turrini

45 pages, 7 figures, 4 tables; accepted as a chapter in the book "Planetary systems now", eds. Luisa M. Lara and David Jewitt, World Scientific Publishing Co Pte Ltd

The great diversity of the thousands of planets known to date is proof of the multitude of ways in which formation and evolution processes can shape the life of planetary systems. Multiple formation and evolution paths, however, can result in the same planetary architecture. Because of this, unveiling the individual histories of planetary systems and their planets can prove a challenging task. The chemical composition of planets provides us with a guiding light for navigate this challenge, but to understand the information it carries we need to properly link it to the chemical composition and characteristics of the environments in which the planets formed. To achieve this goal it is necessary to combine the information and perspectives provided by a growing number of different fields of study, spanning the whole lifecycle of stars and their planetary systems. The aim of this chapter is to provide the unifying perspective needed to understand and connect such diverse information, and illustrate the process through which we can decode the message contained into the composition of planetary bodies.

Andrea Bulgarelli, Antonio Addis, Alessio Aboudan, Ismael Abu, Carla Andreani, Andrea Argan, Riccardo Campana, Paolo Calabretto, Carlotta Pittori, Fabio D'Amico, Imma Donnarumma, Adriano De Rosa, Fabio Fuschino, Giuseppe Gorini, Giuseppe Levi, Nicolò Parmiggiani, Piergiorgio Picozza, Gianluca Polenta, Enrico Preziosi, Roberto Senesi, Alessandro Ursi, Valerio Vagelli, Enrico Virgilli

4 pages, 1 figure, Astronomical Data Analysis Software and System XXXII (2022)

Gamma-Flash is an Italian project funded by the Italian Space Agency (ASI) and led by the National Institute for Astrophysics (INAF), devoted to the observation and study of high-energy phenomena, such as terrestrial gamma-ray flashes and gamma-ray glows produced in the Earth's atmosphere during thunderstorms. The project's detectors and the data acquisition and control system (DACS) are placed at the "O. Vittori" observatory on the top of Mt. Cimone (Italy). Another payload will be placed on an aircraft for observations of thunderstorms in the air. This work presents the architecture of the data acquisition and control system and the data flow.

Shuo Kong, Volker Ossenkopf-Okada, Héctor G. Arce, Ralf S. Klessen, Duo Xu

31 pages, 20 figures, 1 table

In this paper, we carry out a pilot parameter exploration for the collision-induced magnetic reconnection (CMR) mechanism that forms filamentary molecular clouds. Following Kong et al. (2021), we utilize Athena++ to model CMR in the context of resistive magnetohydrodynamics (MHD), considering the effect from seven physical conditions, including the Ohmic resistivity ($\eta$), the magnetic field ($B$), the cloud density ($\rho$), the cloud radius $R$, the isothermal temperature $T$, the collision velocity $v_x$, and the shear velocity $v_z$. Compared to their fiducial model, we consider a higher and a lower value for each one of the seven parameters. We quantify the exploration results with five metrics, including the density probability distribution function ($\rho$-PDF), the filament morphology (250 $\mu$m dust emission), the $B$-$\rho$ relation, the dominant fiber width, and the ringiness that describes the significance of the ring-like sub-structures. The exploration forms straight and curved CMR-filaments with rich sub-structures that are highly variable in space and time. The variation translates to fluctuation in all the five metrics, reflecting the chaotic nature of magnetic reconnection in CMR. A temporary $B\propto\rho$ relation is noticeable during the first 0.6 Myr. Overall, the exploration provides useful initial insights to the CMR mechanism.

Pairs of supermassive black holes (SMBHs) at different stages are natural results of galaxy mergers in the hierarchical framework of galaxy formation and evolution. However, identifications of close binaries of SMBHs (CB-SMBHs) with sub-parsec separations in observations are still elusive. Recently, unprecedented spatial resolutions achieved by GRAVITY/GRAVITY+ onboard The Very Large Telescope Interferometer through spectroastrometry (SA) provide new opportunities to resolve CB-SMBHs. Differential phase curves of CB-SMBHs with two independent broad-line regions (BLRs) are found to have distinguished characteristic structures from a single BLR \citep{songsheng2019}. Once the CB-SMBH evolves to the stage where BLRs merge to form a circumbinary BLR, it will hopefully be resolved by the pulsar timing array (PTA) in the near future as sources of nano-hertz gravitational waves. In this work, we use a parameterized model for circumbinary BLRs to calculate line profiles and differential phase curves for SA observations. We show that both profiles and phase curves exhibit asymmetries caused by the Doppler boosting effect of accretion disks around individual black holes, depending on the orbital parameters of the binary and geometries of the BLR. We also generate mock SA data using the model and then recover orbital parameters by fitting the mock data. Degeneracies between parameters contribute greatly to uncertainties of parameters but can be eased through joint analysis of multiple-epoch SA observations and reverberation mappings.

Sven Kiefer, David Gobrecht, Leen Decin, Christiane Helling

13 pages, 4 figures

Nucleation is considered to be the first step in dust and cloud formation in the atmospheres of asymptotic giant branch (AGB) stars, exoplanets, and brown dwarfs. In these environments dust and cloud particles grow to macroscopic sizes when gas phase species condense onto cloud condensation nuclei (CCNs). Understanding the formation processes of CCNs and dust in AGB stars is important because the species that formed in their outflows enrich the interstellar medium. Although widely used, the validity of chemical and thermal equilibrium conditions is debatable in some of these highly dynamical astrophysical environments. We aim to derive a kinetic nucleation model that includes the effects of thermal non-equilibrium by adopting different temperatures for nucleating species, and to quantify the impact of thermal non-equilibrium on kinetic nucleation. Forward and backward rate coefficients are derived as part of a collisional kinetic nucleation theory ansatz. The endothermic backward rates are derived from the law of mass action in thermal non-equilibrium. We consider elastic collisions as thermal equilibrium drivers. For homogeneous TiO2 nucleation and a gas temperature of 1250 K, we find that differences in the kinetic cluster temperatures as small as 20 K increase the formation of larger TiO2 clusters by over an order of magnitude. An increase in cluster temperature of around 20 K at gas temperatures of 1000 K can reduce the formation of a larger TiO2 cluster by over an order of magnitude. Our results confirm and quantify the prediction of previous thermal non-equilibrium studies. Small thermal non-equilibria can cause a significant change in the synthesis of larger clusters. Therefore, it is important to use kinetic nucleation models that include thermal non-equilibrium to describe the formation of clusters in environments where even small thermal non-equilibria can be present.

E. Antonucci, C. Downs, G. E. Capuano, D. Spadaro, R. Susino, D. Telloni, V. Andretta, V. Da Deppo, Y. De Leo, S. Fineschi, F. Frassetto, F. Landini, G. Naletto, G. Nicolini, M. Pancrazzi, M. Romoli, M. Stangalini, L. Teriaca, M. Uslenghi

The slow solar wind belt in the quiet corona, observed with the Metis coronagraph on board Solar Orbiter on May 15, 2020, during the activity minimum of the cycle 24, in a field of view extending from 3.8 $R_\odot$ to 7.0 $R_\odot$, is formed by a slow and dense wind stream running along the coronal current sheet, accelerating in the radial direction and reaching at 6.8 $R_\odot$ a speed within 150 km s$^{-1}$ and 190 km s$^{-1}$, depending on the assumptions on the velocity distribution of the neutral hydrogen atoms in the coronal plasma. The slow stream is separated by thin regions of high velocity shear from faster streams, almost symmetric relative to the current sheet, with peak velocity within 175 km s$^{-1}$ and 230 km s$^{-1}$ at the same coronal level. The density-velocity structure of the slow wind zone is discussed in terms of the expansion factor of the open magnetic field lines that is known to be related to the speed of the quasi-steady solar wind, and in relation to the presence of a web of quasi separatrix layers, S-web, the potential sites of reconnection that release coronal plasma into the wind. The parameters characterizing the coronal magnetic field lines are derived from 3D MHD model calculations. The S-web is found to coincide with the latitudinal region where the slow wind is observed in the outer corona and is surrounded by thin layers of open field lines expanding in a non-monotonic way.

R. Mckinven, B. M. Gaensler, D. Michilli, K. Masui, V. M. Kaspi, J. Su, M. Bhardwaj, T. Cassanelli, P. Chawla, F. (Adam) Dong, E. Fonseca, C. Leung, E. Petroff, Z. Pleunis, M. Rafiei-Ravandi, I. H. Stairs, S. Tendulkar, D. Z. Li, C. Ng, C. Patel, A. B. Pearlman, M. Rahman, K. R. Sand, K. Shin

Fast radio bursts (FRBs) display a confounding variety of burst properties and host galaxy associations. Repeating FRBs offer insight into the FRB population by enabling spectral, temporal and polarimetric properties to be tracked over time. Here, we report on the polarized observations of 12 repeating sources using multi-year monitoring with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) over 400-800 MHz. We observe significant RM variations from many sources in our sample, including RM changes of several hundred $\rm{rad\, m^{-2}}$ over month timescales from FRBs 20181119A, 20190303A and 20190417A, and more modest RM variability ($\rm{\Delta RM \lesssim}$ few tens rad m$^{-2}$) from FRBs 20181030A, 20190208A, 20190213B and 20190117A over equivalent timescales. Several repeaters display a frequency dependent degree of linear polarization that is consistent with depolarization via scattering. Combining our measurements of RM variations with equivalent constraints on DM variability, we estimate the average line-of-sight magnetic field strength in the local environment of each repeater. In general, repeating FRBs display RM variations that are more prevalent/extreme than those seen from radio pulsars in the Milky Way and the Magellanic Clouds, suggesting repeating FRBs and pulsars occupy distinct magneto-ionic environments.

C. Fian, D. Chelouche, S. Kaspi, C. Sobrino Figaredo, T. Lewis, S. Catalan

arXiv admin note: text overlap with arXiv:2111.07385

We report the results from a photometric reverberation mapping campaign carried out with the C18 telescope at the Wise Observatory from 2019 to 2020, targeting the active galactic nucleus (AGN) MCG 08-11-011. The monitoring was conducted on a daily basis with specially designed narrow-band filters, spanning from optical to near-infrared wavelengths ($\sim4000$ to $8000${\AA}) and avoiding prominent broad emission lines. We aim to measure inter-band continuum time lags, determine the size-wavelength relation, and estimate the host-subtracted AGN luminosity for this system. We used the point-spread function photometry to extract the continuum light curves and measure the inter-band time lags using several methods, including the interpolated cross-correlation function, the z-transformed discrete correlation function, a von Neumann estimator, JAVELIN (in spectroscopic and photometric mode), MICA, and a multivariate correlation function. We find wavelength-dependent lags, $\tau(\lambda)$, up to $\sim$7 days between the multiband light curves of MCG 08-11-011. The observed lags are larger than predictions based on standard thin-disk theory by a factor of $\sim3-7$. We discern a significantly steeper ($\tau \propto \lambda^{4.74}$) size-wavelength relation than the $\tau \propto \lambda^{4/3}$ expected for a geometrically thin and optically thick accretion disk, which may result from the contribution of diffuse continuum emission to the flux. These results are similar to those found by previous continuum reverberation mapping campaigns.

P. Popesso, A. Biviano, E. Bulbul, A. Merloni, J. Comparat, N. Clerc, Z. Igo, A. Liu, S. Driver, M. Salvato, M. Brusa, Y. E. Bahar, N. Malavasi, V. Ghirardini, G. Ponti, A. Robotham, J. Liske, S. Grandis

15 pages, 13 figures, Submitted to MNRAS

The paper presents the analysis of GAMA spectroscopic groups and clusters detected and undetected in the SRG/eROSITA X-ray map of the eFEDS (eROSITA Final Equatorial Depth Survey) area, in the halo mass range $10^{13}-5x10^{14}$ $M_{\odot}$ and at $z < 0.2$. We compare the X-ray surface brightness profiles of the eROSITA detected groups with the mean stacked profile of the undetected low-mass halos. Overall, we find that the undetected groups exhibit less concentrated X-ray surface brightness, dark matter, and galaxy distributions with respect to the X-ray detected halos. Consistently with the low mass concentration, the magnitude gap indicates that these are younger systems. The later assembly time is confirmed by the bluer average color of the BCG and of the galaxy population with respect to the detected systems. They reside with a higher probability in filaments while X-ray detected low-mass halos favor the nodes of the Cosmic Web. Because of the suppressed X-ray central emission, the undetected systems tend to be X-ray under-luminous at fixed halo mass, and to lie below the $L_X-M_{halo}$ relation. Interestingly, the X-ray detected systems inhabiting the nodes scatter the less around the relation, while those in filaments tend to lie below it. We do not observe any strong relation between the properties of detected and undetected systems with the AGN activity. The fraction of optically selected AGN in the galaxy population is consistent in the two samples. More interestingly, the probability that the BCG hosts a radio AGN is lower in the undetected groups. We, thus, argue that the observed differences between X-ray detected and undetected groups are ascribable to the Cosmic Web, and its role in the halo assembly bias. Our results suggest that the X-ray selection is biased to favor the most concentrated and old systems located in the nodes of the Cosmic Web.

Since the first LIGO/Virgo detection, Gravitational Waves (GWs) have been very promising as a new complementary probe to understand our Universe. One of the next challenges of GW search is the detection and characterization of the stochastic gravitational wave background (SGWB), that is expected to open a window on the very early Universe (cosmological background) and to provide us new information on astrophysical source populations (astrophysical background). One way to characterize the SGWB and to extract information about its origin is through the cross-correlation with other cosmological probes. To this aim, in this paper, we explore the cross-correlation between the astrophysical background anisotropies and the Cosmic Microwave Background (CMB) ones. Such a signal is sensitive to primordial non-Gaussianity (nG) through the GW bias. Thus, we study the capability of next generation space-based interferometers to detect such a cross-correlation signal and to constrain primordial nG.

Zsófia V. Kovács-Stermeczky, József Vinkó

17 pages, 12 figures, 2 tables, accepted in PASP

A tidal disruption event (TDE) occurs when a supermassive black hole disrupts a nearby passing star by tidal forces. The subsequent fallback accretion of the stellar debris results in a luminous transient outburst. Modeling the light curve of such an event may reveal important information, for example the mass of the central black hole. This paper presents the TiDE software based on semi-analytic modeling of TDEs. This object-oriented code contains different models for the accretion rate and the fallback timescale $t_{\rm min}$. We compare the resulting accretion rates to each other and with hydrodynamically simulated ones and find convincing agreement for full disruptions. We present a set of parameters estimated with TiDE for the well-observed TDE candidate AT2019qiz, and compare our results with those given by the MOSFiT code. Most of the parameters are in reasonable agreement, except for the mass and the radiative efficiency of the black hole, both of which depend heavily on the adopted fallback accretion rate.

Context. Chromospheric fibrils are thin and elongated structures that connect nearby photospheric magnetic field concentrations of opposite polarities. Aims. We assess the possibilities and drawbacks related to the use of current instrumentation and inversion techniques to infer the thermodynamic structure of chromospheric fibrils. Methods. We employed spectroscopic observations obtained in the Ca ii 854.2 nm line with the CRISP instrument at the Swedish 1-m Solar Telescope and in coordination with observations in the ultraviolet Mg ii h & k lines taken with the IRIS satellite. We studied the temperature sensitivity of these chromospheric lines to properly invert their spectral profiles with the Stockholm inversion Code and determine the temperature, line-of-sight velocity, and microturbulent velocity of manually traced chromospheric fibrils present in the field of view. Results. Fibril-like structures show a very particular dependence of their temperature as a function of the position along their length. Their temperatures at the detected footpoints are, on average, 300 K higher than the temperature at the midpoint. The temperature variation appears to be almost symmetrical in shape, with partially traced fibrils showing a similar trend for the temperature variation. Additionally, the response of the Ca ii 854.2 nm line core to variations of the temperature for the inverted models of the atmosphere in fibril areas seems to be insufficient to properly resolve the aforementioned temperature structure. Only the addition of more temperature sensitive lines such as the Mg ii h & k lines would make it possible to properly infer the thermodynamic properties of chromospheric fibrils. Comparisons between the results obtained here and in previous studies focused on bright Ca ii K fibrils yield great similarities between these structures in terms of their temperature.

Bingjie Wang, Joel Leja, Rachel Bezanson, Benjamin D. Johnson, Gourav Khullar, Ivo Labbe, Sedona H. Price, John R. Weaver, Katherine E. Whitaker

Accepted for publication in ApJL. 13 pages, 6 figures, 2 tables. The code is made publicly available online as a part of Prospector; the version used in this work corresponds to the state of the Git repository at commit 820ad72

The advent of the James Webb Space Telescope (JWST) signals a new era in exploring galaxies in the high-$z$ universe. Current and upcoming JWST imaging will potentially detect galaxies out to $z \sim 20$, creating a new urgency in the quest to infer accurate photometric redshifts (photo-$z$) for individual galaxies from their spectral energy distributions, as well as masses, ages and star formation rates. Here we illustrate the utility of informed priors encoding previous observations of galaxies across cosmic time in achieving these goals. We construct three joint priors encoding empirical constraints of redshifts, masses, and star formation histories in the galaxy population within the \prospector\ Bayesian inference framework. In contrast with uniform priors, our model breaks an age-mass-redshift degeneracy, and thus reduces the mean bias error in masses from 0.3 to 0.1 dex, and in ages from 0.6 to 0.2 dex in tests done on mock JWST observations. Notably, our model recovers redshifts at least as accurately as the state-of-the-art photo-$z$ code \eazy\ in deep JWST fields, but with two advantages: tailoring a model based on a particular survey renders mostly unnecessary given well-motivated priors; obtaining joint posteriors describing stellar, active galactic nuclei, gas, and dust contributions becomes possible. We can now confidently use the joint distribution to propagate full non-Gaussian redshift uncertainties into inferred properties of the galaxy population. This model, ``\prospector-$\beta$'', is intended for fitting galaxy photometry where the redshift is unknown, and will be instrumental in ensuring the maximum science return from forthcoming photometric surveys with JWST. The code is made publicly available online as a part of \prospector.

Katherine de Kleer, Zachariah Milby, Carl Schmidt, Maria Camarca, Michael E. Brown

Published in Planetary Science Journal

The tenuous atmospheres of the Galilean satellites are sourced from their surfaces and produced by a combination of plasma-surface interactions and thermal processes. Though thin, these atmospheres can be studied via their auroral emissions, and most work to date has focused on their aurora at UV wavelengths. Here we present the first detections of Ganymede's and Callisto's optical aurorae, as well detections of new optical auroral lines at Europa, based on observations of the targets over ten Jupiter eclipses from 1998 to 2021 with Keck/HIRES. We present measurements of OI emission at 6300/6364, 5577, 7774, and 8446 A and place upper limits on hydrogen at 6563 A. These constitute the first detections of emissions at 7774 and 8446 A at a planetary body other than Earth. The simultaneous measurement of multiple emission lines provides robust constraints on atmospheric composition. We find that the eclipse atmospheres of Europa and Ganymede are composed predominantly of O2 with average column densities of (4.1 \pm 0.1) x 10^{14} cm^{-2} and (4.7 \pm 0.1) x 10^{14} cm^{-2}, respectively. We find weak evidence for H2O in Europa's bulk atmosphere at an H2O/O2 ratio of $\sim$0.25, and place only an upper limit on H2O in Ganymede's bulk atmosphere, corresponding to H2O/O2 < 0.6. The column density of O2 derived for Callisto is (4.0 \pm 0.9 x 10^{15} cm^{-2} for an assumed electron density of 0.15 cm^{-3}, but electron properties at Callisto's orbit are very poorly constrained.

Julian B. Muñoz

20+2 pages, 21+3 figures. Comments welcome. Code available at this https URL

The 21-cm signal holds the key to understanding the first structure formation during cosmic dawn. Theoretical progress over the last decade has focused on simulations of this signal, given the nonlinear and nonlocal relation between initial conditions and observables (21-cm or reionization maps). Here, instead, we propose an effective and fully analytic model for the 21-cm signal during cosmic dawn. We take advantage of the exponential-like behavior of the local star-formation rate density (SFRD) against densities at early times to analytically find its correlation functions including nonlinearities. The SFRD acts as the building block to obtain the statistics of radiative fields (X-ray and Lyman-$\alpha$ fluxes), and therefore the 21-cm signal. We implement this model as the public Python package Zeus21. This code can fully predict the 21-cm global signal and power spectrum in $\sim 1$ s, with negligible memory requirements. When comparing against state-of-the-art semi-numerical simulations from 21CMFAST we find agreement to 10\% precision in both the 21-cm global signal and power spectra, after accounting for a (previously missed) underestimation of adiabatic fluctuations in 21CMFAST. Zeus21 is modular, allowing the user to vary the astrophysical model for the first galaxies, and interfaces with the cosmological code CLASS, which enables searches for beyond standard-model cosmology in 21-cm data. This represents a step towards bringing 21-cm to the era of precision cosmology.

A large number of galactic binary systems emit gravitational waves (GW) continuously with frequencies below $\sim$10 mHz. The LISA mission could identify the tens of thousands of binaries over years of observation and will be subject to the confusion noise around 1 mHz yielded by the unresolved sources. Beyond LISA, there are several missions have been proposed to observe GWs in the sub-mHz range where the galactic foreground is expected to be overwhelming the instrumental noises. In this study, we investigate the detectability of sub-mHz GW missions to detect the galactic double white dwarf (DWD) binaries and evaluate the confusion noise produced by the undistinguished DWDs. This confusion noise could also be viewed as a stochastic GW foreground and be effectively observed in the sub-mHz band. The parameter determinations for the modeled foreground are examined by employing different detector sensitivities and population models. By assuming the determined foregrounds could be subtracted from the data, we evaluate the residuals which are expected to have power spectral densities two orders of magnitude lower than the originals data.

Hannah Banks, Dorota M. Grabowska, Matthew McCullough

32 pages, 5 figures

Long baseline atom interferometers offer an exciting opportunity to explore mid-frequency gravitational waves. In this work we survey the landscape of possible contributions to the total 'gravitational wave background' in this frequency band and advocate for targeting this observable. Such an approach is complimentary to searches for resolved mergers from individual sources and may have much to reveal about the Universe. We find that the inspiral phases of stellar-mass compact binaries cumulatively produce a signal well within reach of the proposed AION-km and AEDGE experiments. Hypothetical populations of dark sector exotic compact objects, harbouring just a tiny fraction of the dark energy density, could also generate signatures unique to mid- and low-frequency gravitational wave detectors, providing a novel means to probe complexity in the dark sector.

Sulagna Bhattacharya, Basudeb Dasgupta, Ranjan Laha, Anupam Ray

14 pages, 6 figures. Comments welcome

Dark matter from the galactic halo can accumulate in neutron stars and transmute them into sub-2.5 $M_\odot$ black holes if the dark matter particles are heavy, stable, and have interactions with nucleons. We show that non-detection of gravitational waves from mergers of such low-mass black holes can constrain the interactions of asymmetric dark matter particles with nucleons. We find benchmark constraints with LIGO O3 data, viz., $\sigma_{\chi n} \geq {\cal O}(10^{-47})$ cm$^2$ for bosonic DM with $m_\chi\sim$ PeV (or $m_\chi\sim$ GeV, if they can Bose-condense) and $\geq {\cal O}(10^{-46})$ cm$^2$ for fermionic DM with $m_\chi \sim 10^3$ PeV. These bounds depend on the priors on DM parameters and on the currently uncertain binary neutron star merger rate density. However, if null-detection continues with increased exposure over the next decade, LIGO will set remarkable constraints. We find the forecasted sensitivity to heavy asymmetric dark matter to be world-leading, viz., dipping many orders of magnitude below the neutrino floor and completely testing the dark matter solution to missing pulsars in the Galactic center, and demonstrate a windfall science-case for gravitational wave detectors.

In this work we developed a deep learning technique that successfully solves a non-linear dynamic control problem. Instead of directly tackling the control problem, we combined methods in probabilistic neural networks and a Kalman-Filter-inspired model to build a non-linear state estimator for the system. We then used the estimated states to implement a trivial controller for the now fully observable system. We applied this technique to a crucial non-linear control problem that arises in the operation of the LIGO system, an interferometric gravitational-wave observatory. We demonstrated in simulation that our approach can learn from data to estimate the state of the system, allowing a successful control of the interferometer's mirror . We also developed a computationally efficient model that can run in real time at high sampling rate on a single modern CPU core, one of the key requirements for the implementation of our solution in the LIGO digital control system. We believe these techniques could be used to help tackle similar non-linear control problems in other applications.

Kazuyuki Sekizawa, Kentaro Kaba

7 pages, 4 figures (main texts) + 15 pages, 3 figures, 21 tables (supplement) - Submitted to Physical Review Letters

We investigate outer crust compositions for a wide range of magnetic field strengths, up to $B\simeq 4\times10^{18}$ G, employing the latest experimental nuclear masses supplemented with various mass models. The essential effects of the magnetic field are twoholds: 1) Enhancement of electron fraction, which is connected to that of protons via the charge neutrality condition, due to the Landau-Rabi quantization of electron motion perpendicular to the field, namely, neutron-richness is suppressed for a given pressure. As a result, 2) nuclei can exist at higher pressures without dripping out neutrons. By exploring optimal outer-crust compositions from all possible nuclei predicted by theoretical models, we find that neutron-rich heavy nuclei with neutron magic numbers 50, 82, 126, as well as 184, with various proton numbers emerge for $B\gtrsim 10^{17}$ G. Moreover, we show that superheavy nuclei with proton numbers $Z\ge104$, including unknown elements such as $Z=119, 120, 122$, and/or $124$, depending on mass models, may emerge as an equilibrium composition at bottom layers of the outer crust for $B\gtrsim 10^{18}$ G, which are extremely neutron-rich ($N\approx 260$-$287$, i.e., $N/Z\approx 2.2$-$2.4$). We point out that those extremely neutron-rich superheavy nuclei locate around the next neutron magic number $N=258$ after $N=184$, underlining importance of nuclear structure calculations under such really exotic, extreme conditions. We demonstrate how the superstrong magnetic field substantially alters crustal properties of neutron stars, which may have detectable consequences.

Matthias Schartner, Leonid Petrov, Christian Plötz, Frank G. Lemoine, Eusebio Terrazas, Benedikt Soja

8 pages, 5 figures

In this work, we present a status update and results of the designated research and development VLBI Intensive program VGOS-INT-S, observed between MACGO12M and WETTZ13S for the rapid determination of the Earth's phase of rotation, expressed via UT1-UTC. The main novelty of these sessions is the use of a special observation strategy, rapidly alternating between high- and low-elevation scans, enabling an improved determination of delays caused by the neutral atmosphere. Since 2021, 25 Intensive sessions have been observed successfully. In early 2022, VGOS-INT-S was among the most accurate Intensive programs with an average formal error $\sigma_{UT1-UTC}$ of 3.1 $\mu$s and a bias w.r.t. IERS C04 of 1.1 $\mu$s. Later, the session performance decreased due to multiple technical difficulties.

José A. R. Cembranos, Miguel González Ortiz, Prado Martín-Moruno

12 pages, 4 figures

Graviton-photon oscillation is the conversion of gravitational waves to electromagnetic waves and vice versa in the presence of a background electromagnetic field. We investigate this phenomenon in a cosmological scenario considering a background cosmic magnetic field and assuming different gravitational frameworks. We obtain the damping term that characterizes the attenuation of the conversion probability in cosmological backgrounds. This is a general feature that is present even for standard General Relativity. Furthermore, we show that the effects of decoherence, which are due to the interaction with the cosmological expansion and with the additional degrees of freedom of alternative theories of gravity, can be relevant to the phenomenon of graviton-photon mixing.

Piyabut Burikham, Tiberiu Harko, Kulapant Pimsamarn, Shahab Shahidi

20 pages, 2 figures, accepted for publication in PRD

We investigate the possibility that the observed behavior of test particles outside galaxies, which is usually explained by assuming the existence of dark matter, is the result of the dynamical evolution of particles in a Weyl type geometry, and its associated conformally invariant Weyl geometric quadratic gravity. As a first step in our investigations we write down the simplest possible conformally invariant gravitational action, constructed in Weyl geometry, and containing the Weyl scalar, and the strength of the Weyl vector only. By introducing an auxiliary scalar field, the theoretical model can be reformulated in the Riemann geometry as scalar-vector-tensor theory, containing a scalar field, and the Weyl vector, respectively. The field equations of the theory are derived in the metric formalism, in the absence of matter. A specific static, spherically symmetric model, in which the Weyl vector has only a radial component, is considered. In this case, an exact analytic solution of the gravitational field equations can be obtained. The behavior of the galactic rotation curves is also considered in detail, and it is shown that an effective geometric mass term, with an associated density profile, can also be introduced. Three particular cases, corresponding to some specific functional forms of the Weyl vector, are also investigated. A comparison of the model with a selected sample of galactic rotation curves is also performed when an explicit breaking of conformal invariance is introduced, which allows the fix of the numerical values of the free parameters of the model. Our results show that Weyl geometric models can be considered as a viable theoretical alternative to the dark matter paradigm.

The flux-based statistical theory of the non-hierarchical three-body system predicts that the chaotic outcome distribution reduces to the chaotic emissivity function times a known function, the asymptotic flux. Here, we measure the chaotic emissivity function (or equivalently, the absorptivity) through simulations. More precisely, we follow millions of scattering events only up to the point when it can be decided whether the scattering is regular or chaotic. In this way, we measure a tri-variate absorptivity function. Using it, we determine the flux-based prediction for the chaotic outcome distribution over both binary binding energy and angular momentum, and we find good agreement with the measured distribution. This constitutes a detailed confirmation of the flux-based theory, and demonstrates a considerable reduction in computation to determine the chaotic outcome distribution.

Shawn G. Rosofsky, E. A. Huerta

13 pages, 9 figures, 1 table. First application of physics informed neural operators to solve magnetohydrodynamics equations

We present the first application of physics informed neural operators, which use tensor Fourier neural operators as their backbone, to model 2D incompressible magnetohydrodynamics simulations. Our results indicate that physics informed AI can accurately model the physics of magnetohydrodynamics simulations that describe laminar flows with Reynolds numbers $Re\leq250$. We also quantify the applicability of our AI surrogates for turbulent flows, and explore how magnetohydrodynamics simulations and AI surrogates store magnetic and kinetic energy across wavenumbers. Based on these studies, we propose a variety of approaches to create AI surrogates that provide a computationally efficient and high fidelity description of magnetohydrodynamics simulations for a broad range of Reynolds numbers. Neural operators and scientific software to produce simulation data to train, validate and test our physics informed neural operators are released with this manuscript.

Saturons are macroscopic objects with maximal microstate entropy. Due to this property, they can be produced via quantum transitions from a homogeneous thermal bath, bypassing the standard exponential suppression characteristic of ordinary extended objects. In this sense, saturons carry an advantage with respect to other macroscopic objects such as black holes and ordinary solitons. Due to unsuppressed thermal production, saturons can have interesting cosmological implications. In particular they can serve as viable dark matter candidates with some unique features. Unlike ordinary particle dark matter, the superheavy saturons can freeze-in at very low temperatures. A nucleation of a saturon can be described in terms of a saturated instanton. This has implications for various phase transitions.

A non-perturbative framework is provided to connect QCD with nuclear phenomenology in the intermediate and high density regime. Using QCD Sum Rules, in-medium scalar and vector self-energies of nucleons are calculated as functions of the density of an infinite nuclear medium. The self-energies are used in the relativistic mean field theory lagrangian of a high-density nuclear medium to find the binding energy of in-medium nucleons and the value of light quark condensate, $\langle \bar{q} q \rangle_{\rm{vac}} = -~(0.288 ~\rm{GeV})^3$, in the Borel-improved resummation scheme. The critical mass of an ideal neutron star is obtained by coupling a uniform saturation energy density of cold, dense nuclear matter to Einstein equation in hydrostatic equilibrium. Since it is less likely for a neutron star core to avoid deconfinement and enter the rigid vector repulsion phase where the speed of sound can smoothly approach from conformal to causal limit, a gap should exist in the stellar mass spectrum, $[3.48M_\odot, 5.47M_\odot]$, where it would be rare to find any isolated, cold, non-rotating neutron star or a black hole.

Sanjay Bhattarai, Dongming Mei, Rajendra Panth, Mathbar Raut, Kyler Kooi, Hao Mei, Guojian Wang

7 pages, 8 figures

We investigated charge transport in an n-type germanium detector at 5.2 K to explore new technology for enhancing low-mass dark matter detection sensitivity. Calculations of dipole and cluster dipole state binding energies and electric field-dependent trapping cross-sections are critical to developing low-threshold detectors. The detector operates in two modes: depleting at 77K before cooling, or directly cooling to 5.2 K and applying different bias voltages. Results indicated lower binding energy of charge states in the second mode, at zero field and under an electric field, suggesting different charge states formed under different operating modes. Measured cluster dipole and dipole state binding energies at zero field were 7.884$\pm$0.644 meV and 8.369$\pm$0.748 meV, respectively, signifying high low-threshold potential for low-mass dark matter searches in the future.

Hoa Dinh Thi, Chiranjib Mondal, Francesca Gulminelli

Contribution to the "Journees de Rencontre des Jeunes Chercheurs (JRJC) 2021" proceedings

The nuclear matter equation of state is relatively well constrained at sub-saturation densities thanks to the knowledge from nuclear physics. However, studying its behavior at supra-saturation densities is a challenging task. Fortunately, the extraordinary progress recently made in observations of neutron stars and neutron star mergers has provided us with unique opportunities to unfold the properties of dense matter. Under the assumption that nucleons are the only constituents of neutron star cores, we perform a Bayesian inference using the so-called meta-modeling technique with a nuclear-physics-informed prior. The latest information from the GW170817 event by the LIGO-Virgo Collaboration (LVC) and from the radius measurement of the heaviest known neutron star PSR J0740+6620 by the Neutron Star Interior Composition Explorer (NICER) telescope and X-ray Multi-Mirror (XMM-Newton) are taken into account as likelihoods in the analysis. The impacts of different constraints on the equation of state as well as on the predictions of neutron star properties are discussed. The obtained posterior reveals that all the current observations are fully compatible with the nucleonic hypothesis. Strong disagreements between our results with future data can be identified as a signal for the existence of exotic degrees of freedom.