Papers for Thursday, Jun 06 2024

Neutron stars in Low Mass X-ray Binaries (LMXBs) can accrete matter onto their surface from the companion star. Transiently accreting neutron stars go through alternating phases of active accretion outbursts and quiescence. X-ray observations during the quiescence phase show a drop in X-ray luminosity with the time in quiescence. This is also inferred as the drop in surface temperature or the cooling of accreting neutron stars in quiescence. Analyzing these cooling curves reveals a great deal of information about the structure and composition of neutron stars. However, model-observation comparisons of such cooling curves are challenging - partly due to observational uncertainties, and partly due to incomplete knowledge of heating mechanisms during accretion outbursts. This situation is further exacerbated by the recent discovery of Urca cooling in the neutron star crust. These are cycles that alternate between electron-capture and beta-decay to produce a large flux of neutrinos and anti-neutrinos. These freely stream out of the star and carry energy with them, essentially cooling the neutron star crust without changing the composition. As a result, it is necessary to accurately quantify the strength of Urca cooling to constrain the heat sources in neutron star crusts and facilitate better model-observation comparisons of the cooling curves.

Radio galaxies harbouring jetted active galactic nuclei are a frequent target of very-long-baseline interferometry (VLBI) because they play an essential role in exploring how jets form and propagate. Hence, only few have not been detected with VLBI yet; Fornax A is one of the most famous examples. Here we present the first detection of the compact core region of Fornax A with VLBI. At 8.4 GHz the faint core is consistent with an unresolved point source. We constrained its flux density to be $S_0 = 47.5-62.3\,\textrm{mJy}$ and its diameter to be $D^\textrm{min}_0 \leq 70\,\mu\textrm{as}$. The high values of the measured brightness temperature ($T_\textrm{B} \gtrsim 10^{11}\,\textrm{K}$) imply that the observed radiation is of non-thermal origin, likely associated with the synchrotron emission from the active galactic nucleus. We also investigated the possibility of a second radio source being present within the field of view. Adding a second Gaussian component to the geometrical model-fit does not significantly improve the quality of the fit and we, therefore, conclude that our detection corresponds to the compact core of Fornax A. Analysis of the non-trivial closure phases provides evidence for the detection of more extended flux density, on the angular scale of $\sim4000\,\mu\textrm{as}$. Finally, the fractional circular polarisation of the core is consistent with zero, with a conservative upper limit being $m_\textrm{circ} \leq 4\%$.

Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenching definitions and find that both simulations significantly underestimate the comoving number density of quenched massive galaxies at $z \gtrsim 3$ compared to JWST observations by up to $\sim 2$ dex. This fact highlights the necessity for improved physical models of AGN feedback in galaxy formation simulations. In both simulations, the high-$z$ quenched massive galaxies often host overmassive central black holes above the standard $M_{BH}-M_\star$ relation, implying that the AGN feedback plays a crucial role in quenching galaxies in the early Universe. The typical quenching timescales for these galaxies are $\sim 200-600$ Myr. IllustrisTNG primarily employs AGN kinetic feedback, while ASTRID relies on AGN thermal feedback, which is less effective and has a longer quenching timescale. We also study the post-quenching evolution of the high-$z$ massive quiescent galaxies and find that many experience subsequent reactivation of star formation, evolving into primary progenitors of $z=0$ brightest cluster galaxies.

Observations with the James Webb Space Telescope have revealed a high abundance of bright galaxies at redshift, $z\gtrsim 12$, which has been widely interpreted as conflicting with the $\Lambda$CDM model. In Cowley et al. (2018) predictions were made - prior to the JWST observations - for the expected abundance of these galaxies using the Durham semi-analytic galaxy formation model, GALFORM, which is known to produce a realistic population of galaxies at lower redshifts including the present day. Key to this model is the assumption of a "top-heavy" initial mass function of stars formed in bursts (required to explain the number counts and redshift distribution of sub-millimetre galaxies). Here, we compare the rest-frame ultraviolet luminosity functions derived from JWST observations with those predicted by the Cowley et al. model up to $z=14$ and make further predictions for $z=16$. We find that below $z\sim 10$, the Cowley et al. predictions agree very well with observations, while agreement at $z\gtrsim12$ requires extending the model to take into account the timescale for the growth of obscuring dust grains and its dependence on gas metallicity. We trace the evolution of these galaxies from $z=14$ to $z=0$ and find that their descendants typically reside in halos with a median mass of $10^{13.6}\,h^{-1}\,\mathrm{M_{\odot}}$. The stellar masses of the descendants range from $10^{7}\,h^{-1}\,\mathrm{M_{\odot}}$ to $10^{11.5}\,h^{-1}\,\mathrm{M_{\odot}}$. Although these galaxies were all central galaxies at $z=14$, nearly half of their descendants end up as satellites in massive halos.

Realization of the stationary integral solutions of steady state transonic accretion flow in spherical symmetry helps to understand accretion phenomena on various astrophysical objects. In recent years, attempts have been made to study accreting black hole systems as an example of autonomous dynamical systems. The fixed point analysis is used to study the transonic properties of accretion flow onto an astrophysical black hole, hence the nature of the phase orbits for the transonic flow solutions can be understood without constructing the integral solutions. Since a large-scale astrophysical fluid flow is vulnerable to external perturbation, one needs to ensure that the stationary accretion solutions are stable under perturbation. By adopting a time-dependent stability analysis scheme for the accretion flow, one demonstrates under which condition the perturbation will not diverge. It has also been observed that a space-time metric, dubbed the sonic metric, can be constructed to describe the propagation of the perturbation embedded within the accreting fluid, which mimics a black hole like spacetime within the accreting fluid, where the transonic surfaces can be identified with a black hole like horizons. Such identification is accomplished using the theory of causal structure - by constructing Carter-Penrose diagrams. An accreting black hole system, thus, can be perceived as a classical analogue gravity model naturally found in the universe. Hence, accretion phenomena onto astrophysical black holes can be looked upon from three apparently non-overlapping perspectives - astrophysical processes, theory of dynamical systems, and emergent gravity (alternatively, the analogue gravity) phenomena, respectively. The present article illustrates, by taking the simplest possible accretion flow model, how one can study astrophysical accretion processes from the aforementioned perspectives (Abridged).

We report the discovery of an extended emission line region (EELR) in MUSE observations of Markarian 950, a nearby ($z=0.01628$) post-starburst (PSB) galaxy that hosted the tidal disruption event (TDE) iPTF-16fnl. The EELR requires a non-stellar ionizing continuum with a luminosity L$_{\rm ion, min} \gtrsim 10^{43}$ erg s$^{-1}$, inconsistent with the current weak state (L$_{\rm IR,AGN} < 2.5 \times 10^{42}$ erg s$^{-1}$) of the galactic nucleus. The ionized gas has low velocity ($\sim$-50 km s$^{-1}$) and low turbulence ($\sigma_{\rm gas} \lesssim$ 50 km s$^{-1}$), and is kinematically decoupled from the stellar motions, indicating that the gas kinematics are not AGN driven. Markarian 950 is the third post-starburst galaxy to host a weak nuclear ionizing source as well as an EELR and a TDE. The overall properties of these three galaxies, including the kinematics and accretion history, are unusual but strikingly similar. We estimate that the incidence of EELRs in PSB-TDE hosts is a factor of $\sim 10 \times$ higher than in other PSB galaxies. This suggests that a gas-rich post-merger environment is a key ingredient in driving elevated TDE rates. Based on the current observations, we cannot rule out that the EELRs may be powered through an elevated TDE rate in these galaxies. If the EELRs are not TDE-powered, the presence of intermittent AGN activity, and in particular the fading of the AGN, may be associated with an increased TDE rate and/or an increased rate of detecting TDEs.

In the past five years, six quasi-periodic X-ray eruption (QPE) sources have been discovered in the nuclei of nearby galaxies. Their origin remains an open question. We present MUSE integral field spectroscopy of five QPE host galaxies to characterize their properties. We find that 3/5 galaxies host extended emission line regions (EELRs) up to 10 kpc in size. The EELRs are photo-ionized by a non-stellar continuum, but the current nuclear luminosity is insufficient to power the observed emission lines. The EELRs are decoupled from the stars both kinematically and in projected sky position, and the low velocities and velocity dispersions ($<$ 100 km s$^{-1}$ and $\lesssim 75$ km s$^{-1}$ respectively) are inconsistent with being AGN- or shock-driven. The origin of the EELRs is likely a previous phase of nuclear activity. The QPE host galaxy properties are strikingly similar to those of tidal disruption events (Wevers et al. submitted). The preference for a very short-lived (the typical EELR lifetime is $\sim$15000 years), gas-rich phase where the nucleus has recently faded significantly suggests that TDEs and QPEs may share a common formation channel, disfavoring AGN accretion disk instabilities as the origin of QPEs. In the assumption that QPEs are related to extreme mass ratio inspiral systems (EMRIs; stellar-mass objects on bound orbits about massive black holes), the high incidence of EELRs and recently faded nuclear activity can be used to aid in the localization of the host galaxies of EMRIs discovered by low frequency gravitational wave observatories.

Debris disks give us the unique opportunity to probe the properties of small $\mu$m-sized particles, allowing us to peer into the constituents of their parent bodies, young analogs of comets and asteroids of our solar system. In the past, studies of the total intensity phase function have proven powerful to constrain the main characteristics of the dust particles in debris disks. Nonetheless, there can remain some degeneracies in the modeling that can be alleviated when considering polarized intensity observations. We obtained new near-IR scattered light observations of four young debris disks which allow us to constrain the degree of linear polarization as a function of the scattering angle. All four debris disks are detected in polarized intensity, and three are also recovered in total intensity. We measured peak degree of polarization of $\lesssim 40$\% for all three disks. We find that the particles must consist of highly refractive and absorbing material. For HD129590, by measuring the polarization fraction beyond the birth ring, we constrain the width of the size distribution to be smaller and smaller, compatible with the effect of radiation pressure. We put these findings to the test and present a self-consistent approach to produce synthetic images, assuming different profiles for the radiation pressure strength, and accounting for the presence of unbound grains. We find the contribution of these grains to be especially critical to reproduce the increasing degree of polarization with stellocentric distances. Some of our results might seem difficult to reconcile with our understanding of cosmic dust but since similar results have been obtained for other disks, we discuss the current limitation of available light scattering models as well as possible avenues to alleviate these unfortunate limitations.

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar~WGD$\,$2038$-$4008~using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{\Delta t} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $\Lambda$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.

The properties of H-rich, type II-plateau supernova (SN II-P) progenitors remain uncertain, and this is primarily due to the complexities associated with red supergiant (RSG) wind-mass loss. Recent studies have suggested that the interaction of the ejecta with a standard RSG wind should produce unambiguous signatures in the optical (e.g., a broad, boxy H${\alpha}$ profile) and in the UV (especially Ly ${\alpha}$ and Mg ii ${\lambda}{\lambda}$ 2795, 2802) a few years following the explosion. Such features are expected to be generic in all SNe II-P and can be utilized to constrain RSG winds. Here, we investigate the possibility of detecting late-time (0.3-10 years since explosion) SNe II-P in the NUV with the China Space Station Telescope (CSST). Convolving the existing model spectra of ejecta-wind interactions in SNe II-P with the transmission functions of the CSST, we calculated the associated multiband light curves, in particular, the NUV (255 nm${\sim}$317 nm) band, as well as the $NUV-r$ color. We find that the CSST will be able to detect the NUV radiation associated with ejecta-wind interaction for hundreds SNe II-P out to a few hundred Mpc over its ten-year main sky survey. The CSST will therefore provide a sizable sample of SNe II-P with the NUV signatures of ejecta-wind interaction. This will be helpful for understanding the mass loss history of SN II-P progenitors and their origins.

Theoretical models of structure formation predict the presence of a hot gaseous atmosphere around galaxies. While this hot circum-galactic medium (CGM) has been observationally confirmed through UV absorption lines, the detection of its direct X-ray emission remains scarce. We investigate theoretical predictions of the intrinsic CGM X-ray surface brightness (SB) using simulated galaxies and connect them to their global properties such as gas temperature, hot gas fraction and stellar mass. We select a sample of galaxies from the ultra-high resolution ($48\ \rm{cMpc\, h^{-1}}$) cosmological volume of the Magneticum Pathfinder set of hydrodynamical cosmological simulations. We classify them as star-forming (SF) or quiescent (QU) based on their specific star-formation rate. For each galaxy we generate X-ray mock data using the X-ray photon simulator PHOX, from which we obtain SB profiles out to the virial radius for different X-ray emitting components, namely gas, active galactic nuclei and X-ray binaries (XRBs). We fit a $\beta$-profile to each galaxy and observe trends between its slope and global quantities of the simulated galaxy. We find marginal differences between the average total SB profile of the CGM in SF and QU galaxies, with the contribution from hot gas being the largest ($>50\%$) at radii $r>0.05\,R_{\rm{vir}}$. The contribution from X-ray binaries (XRBs) equals the gas contribution for small radii and is non-zero for large radii. The galaxy population shows positive correlations between global properties and normalization of the SB profile. The slope of fitted $\beta$-profiles correlates strongly with the total gas luminosity, which in turn shows strong connections to the current accretion rate of the central super-massive black hole (SMBH).

We study the post-inflationary energy transfer from the inflaton ($\phi$) into a scalar field ($\chi$) non-minimally coupled to gravity through $\xi R|\chi|^2$, considering models with inflaton potential $V_{\rm inf} \propto |\phi|^{\,p}$ around $\phi = 0$. This corresponds to the paradigm of {\it geometric preheating}, which we extend to its non-linear regime via lattice simulations. Considering $\alpha$-attractor T-model potentials as a proxy, we study the viability of proper {\it reheating} for $p=2, 4, 6$, determining whether radiation domination (RD) due to energetic dominance of $\chi$ over $\phi$, can be achieved. For large inflationary scales $\Lambda$, reheating is frustrated for $p = 2$, it can be partially achieved for $p = 4$, and it becomes very efficient for $p = 6$. Efficient reheating can be however blocked if $\chi$ sustains self-interactions (unless these are extremely feeble), or if $\Lambda$ is low enough, so that inflaton fragmentation brings the universe rapidly into RD. Whenever RD is achieved, either due to reheating or to inflaton fragmentation, we characterize the energy and time scales of the problem, as a function of $\Lambda$ and $\xi$.

I present a model for light travel time effects for emission from a plasma blob in a blazar jet. This calculation could be incorporated into more complex models with particle acceleration and radiation mechanisms, but as presented here it is a agnostic as to these mechanisms. This model includes light travel time effects for an expanding or contracting blob. As an example, this model is applied to a flare observed by VERITAS and MAGIC from Mrk 421 in 2013; and to a flare observed by the Fermi Large Area Telescope from 3C 454.3 in 2010.

We present stellar parameters and chemical abundances of 47 elements detected in the bright (V = 11.63) very metal-poor ([Fe/H] = -2.20 +- 0.12) star 2MASS J22132050-5137385. We observed this star using the Magellan Inamori Kyocera Echelle spectrograph as part of ongoing work by the R-Process Alliance. The spectrum of 2MASS J22132050-5137385 exhibits unusually strong lines of elements heavier than the iron group, and our analysis reveals that these elements were produced by rapid neutron-capture (r-process) nucleosynthesis. We derive a europium enhancement, [Eu/Fe] = +2.45 +- 0.08, that is higher than any other r-process-enhanced star known at present. This star is only the eighth r-process-enhanced star where both thorium and uranium are detected, and we calculate the age of the r-process material, 13.6 +- 2.6 Gyr, from the radioactive decay of these isotopes. This star contains relatively large enhancements of elements that may be produced as transuranic fission fragments, and we propose a new method using this characteristic to assess the r-process yields and gas dilution in samples of r-process-enhanced stars. We conclude that 2MASS J22132050-5137385 exhibits a high level of r-process enhancement because it formed in an environment where the r-process material was less diluted than average. Assuming a canonical baryonic minihalo mass of 10^6 M_sun and a 1 percent metal retention rate, this star formed in a cloud of only ~ 600 M_sun.

We present the online service cosmICweb (COSMological Initial Conditions on the WEB) - the first database and web interface to store, analyze, and disseminate initial conditions for zoom simulations of objects forming in cosmological simulations: from galaxy clusters to galaxies and more. Specifically, we store compressed information about the Lagrangian proto-halo patches for all objects in a typical simulation merger tree along with properties of the halo/galaxy across cosmic time. This enables a convenient web-based selection of the desired zoom region for an object fitting user-specified selection criteria. The information about the region can then be used with the MUSIC code to generate the zoom ICs for the simulation. In addition to some other simulations, we currently support all objects in the EAGLE simulation database, so that for example the Auriga simulations are easily reproduced, which we demonstrate explicitly. The framework is extensible to include other simulations through an API that can be added to an existing database structure and with which cosmICweb can then be interfaced. We make the web portal and database publicly available to the community.

Recent results from pulsar timing arrays (PTAs) show evidence for a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole (SMBH) binaries (BHBs). While the data do not yet constrain the slope of the spectrum, this appears to flatten at the lowest frequencies, deviating from the power-law shape expected for circular binaries evolving solely due to gravitational wave (GW) emission. Interestingly, such flattening can be explained with a population of eccentric rather than circular binaries. The eccentricity of BHBs is notoriously difficult to predict based simply on the parameters of the host galaxies and the initial galactic orbit, as it is subject to stochastic effects. We study the evolution of the eccentricity of BHBs formed in galactic mergers with cosmological initial conditions from pairing to coalescence, with a focus on potential PTA sources. We select galactic mergers from the IllustrisTNG100-1 simulation and re-simulate them at high resolution with the N-body code Griffin down to binary separations of the order of a parsec. We then estimate coalescence timescales with a semi-analytical model of the evolution under the effects of GW emission and stellar hardening. We find that most mergers in IllustrisTNG100-1 occur on highly eccentric orbits, and that the eccentricity of BHBs at binary formation correlates with the initial eccentricity of the merger, if this is no larger than approximately 0.9. For extremely eccentric mergers, the binaries tend to form with modest eccentricities. We discuss the implications of these results on the interpretation of the observed GWB.

\textit{LiteBIRD}, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, \textit{LiteBIRD} will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$\mu$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $\delta r = 0.001$, including systematic errors and margin. If $r \geq 0.01$, \textit{LiteBIRD} expects to achieve a $>5\sigma$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe \textit{LiteBIRD}'s scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight \textit{LiteBIRD}'s synergies with concurrent CMB projects.

Observations show evidence that a significant fraction of protoplanetary disks contain warps. A warp in a disk evolves in time affecting the appearance of shadows and greatly influencing kinematic signatures. So far, many theoretical studies of warped disks have been conducted using Smoothed Particle Hydrodynamics (SPH) methods. In our approach, we use a grid-based method in spherical coordinates which has notable advantages: the method allows for accurate modelling of low viscosity values and the resolution does not depend on density or mass of the disk, which allows surface structures to be resolved. We perform 3D simulations using FARGO3D to simulate the evolution of a warped disk and compare the results to one-dimensional models using a ring code. Additionally, we extensively investigate the applicability of grid-based methods to misaligned disks and test their dependency on grid resolution as well as disk viscosity. We find that grid-based simulations are capable of simulating disks not aligned to the grid geometry. Our three-dimensional simulation of a warped disk compares well to one-dimensional models in evolution of inclination. However, we find a twist which is not captured in 1D models. After thorough analysis we suspect this to be a physical effect possibly caused by non-linear effects neglected in the one-dimensional equations. Evaluating the internal dynamics, we find sloshing and breathing motions as predicted in local shearing box analysis. They can become supersonic, which may have consequences on kinematic observations of warped disks. Warped disks can be accurately modelled in 3D grid-based simulations when using reasonably good resolution, especially in the $\theta$-direction. We find good agreement with the linear approximation of the sloshing motion which highlights the reliability of 1D models.

Understanding the formation of multiple populations in globular clusters (GCs) represents a challenge for stellar population studies. Nevertheless, the outermost GC regions, which may retain information of the initial configuration of the multiple populations, are poorly studied. We use synthetic spectra that account for the chemical compositions of first- and second-population (1P, 2P) stars of 47 Tucanae to identify the spectral regions that are sensitive to its multiple populations. Hence, we defined new photometric bands that are efficient to disentangle 1P and 2P giant stars from Gaia XP spectra. To test these new filters, we constructed the pseudo two-color diagrams dubbed chromosome maps (ChMs) and identified for the first time, 1P and 2P stars in the outermost cluster regions of 47 Tucanae and outside the tidal radius. We constructed similar diagrams for NGC3201, NGC6121, NGC6752, and NGC6397, thus exploring GCs with different metallicities. The ChMs allowed us to clearly disentangle 1P and 2P stars in the external regions of all clusters, with the exception of the metal-poor NGC6397. Our findings, together with literature results from more-internal regions, show that the 2P stars of 47 Tucanae and NGC 3201 are more-centrally concentrated than the 1P, whereas the multiple populations of NGC 6121, and NGC 6752 share the same radial distributions. These radial behaviors are consistent with the GC formation scenarios where 2P stars originate in the central regions. Noticeably, results on NGC 3201 are in tension with the conclusion from recent work that its 1P is more centrally concentrated than the 2P and might form with more central concentration.

Radio galaxies are a subclass of active galactic nuclei that drive relativistic jets from their center and are observed in radio to very-high-energy gamma rays. The emission mechanisms and regions are still unknown. High-energy gamma rays can be explained by the emission from the magnetically arrested disks (MADs) around the central supermassive black hole, for which the magnetic flux threading the black hole is in a saturation level, although the emission from the MADs does not explain the optical and X-ray data. We construct a multi-wavelength emission model in which the optical and X-ray emission is emitted by jets and the gamma rays by MADs. Our model takes into account the particle injection by the magnetic reconnection at the jet base close to the black hole and particle entrainment from the ambient gas at the jet emission zone. We apply our model to M87 and find that our model can explain the simultaneous multi-wavelength data. In our model, the emission from the jets is the synchrotron radiation of the nonthermal electrons accelerated by magnetic reconnection, and the emission from the MADs is the synchrotron radiation mainly of the nonthermal protons accelerated by turbulence. We also find that the strong plasma entrainment is necessary to explain the multi-wavelength data. Our model will be tested by variability analysis among the multi-wavelength data.

In the local Universe, the ratio between the mass of a central black hole and the stellar mass of its host galaxy is $\sim 0.1 \%$. Recently, JWST discovered numerous galaxies at $z>4$ that seem to deviate from the local relation, with black holes overmassive by $10-100$ times. Similar galaxies were also discovered at cosmic noon. The intrinsic scatter in the relation describes how much the evolutionary histories of the single galaxies deviate from the mean evolutionary pattern of their parent dataset. This Research Note examines whether a cosmic evolution of the intrinsic scatter can be detected by assessing its value for datasets in various redshift ranges. Using data from the local Universe ($z < 0.055$), low ($ 0.4 \leq z \leq 0.9$), intermediate ($ 0.9 \leq z \leq 4$), and high ($z > 4$) redshift, we conclude that there is no statistically significant redshift evolution of the intrinsic scatter.

Reconstructing images from the Event Horizon Telescope (EHT) observations of M87*, the supermassive black hole at the center of the galaxy M87, depends on a prior to impose desired image statistics. However, given the impossibility of directly observing black holes, there is no clear choice for a prior. We present a framework for flexibly designing a range of priors, each bringing different biases to the image reconstruction. These priors can be weak (e.g., impose only basic natural-image statistics) or strong (e.g., impose assumptions of black-hole structure). Our framework uses Bayesian inference with score-based priors, which are data-driven priors arising from a deep generative model that can learn complicated image distributions. Using our Bayesian imaging approach with sophisticated data-driven priors, we can assess how visual features and uncertainty of reconstructed images change depending on the prior. In addition to simulated data, we image the real EHT M87* data and discuss how recovered features are influenced by the choice of prior.

Binaries in the Kuiper Belt are common. Here we present our analysis of the Solar System Origins Legacy Survey (SSOLS) to show that using a PSF-fitting method can roughly double the number of binaries identified in that dataset. Out of 198 Kuiper Belt objects (KBOs) observed by SSOLS, we find 23 to be visually separated binaries, while a further 19 are blended-PSF binaries detectable with the method we present here. This is an overall binary fraction of 21% for the SSOLS dataset of cold classical KBOs. In addition, we tested our fitting methods on synthetic data, and while we were able to show it to be very effective at detecting certain blended-PSF binary KBOs, fainter or closer binary KBOs may easily be missed, suggesting that the close binary KBO fraction could be even higher. These results strongly support the idea that most (if not all) KBOs were formed through the Streaming Instability process, and as a consequence, most KBOs were formed as near-equal mass binaries.

We present results from a five-month-long observing campaign of the unusual transient AT 2022fnm, which displays properties common to both luminous red novae (LRNe) and intermediate-luminosity red transients (ILRTs). Although its photometric evolution is broadly consistent with that of LRNe, no second peak is apparent in its light curve, and its spectral properties are more reminiscent of ILRTs. It has a fairly rapid rise time of 5.3$\pm$1.5 d, reaching a peak absolute magnitude of $-12.7\pm$0.1 (in the ATLAS $o$ band). We find some evidence for circumstellar interaction, and a near-infrared excess becomes apparent at approximately +100 d after discovery. We attribute this to a dust echo. Finally, from an analytical diffusion toy model, we attempted to reproduce the pseudo-bolometric light curve and find that a mass of $\sim$4 M$_\odot$ is needed. Overall, the characteristics of AT 2022fnm are consistent with a weak stellar eruption or an explosion reminiscent of low-energy type IIP supernovae, which is compatible with expectations for ILRTs.

Dynamo processes are ubiquitous in astrophysical systems. In relativistic astrophysical systems, such as accretion disks around black holes or neutron stars, they may critically affect the launching of winds and jets that can power electromagnetic emission. Dynamo processes are governed by several microscopic parameters, one of them being magnetic helicity. As a conserved quantity in nonresistive plasmas, magnetic helicity is transported across the system. One important implication of helicity conservation is, that in the absence of helicity fluxes some mean-field dynamos can be quenched, potentially affecting the large-scale evolution of the magnetic field. One of the major challenges in computing magnetic helicity is the need to fix a meaningful electromagnetic gauge. We here present a fully covariant formulation of magnetic helicity transport in general-relativistic plasmas based on the concept of relative helicity by Berger & Field and Finn & Antonsen. This formulation is separately invariant under gauge-transformation of the Maxwell and Einstein equations. As an application of this new formalism we present the first analysis of magnetic helicity transport in the merger of two neutron stars. We demonstrate the presence of global helicity fluxes into the outer layers of the stellar merger remnant, which may impact subsequent large-scale dynamo amplification in these regions.

We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity.

Galaxy clusters are crucial to understanding role of the environment in galaxy evolution. However, due to their rarity, only a limited number of clusters have been identified at $z\gtrsim2$. In this paper, we report a discovery of seven cluster candidates with massive quiescent galaxies at $z\sim2$ in the $3.5\,\mathrm{deg}^{2}$ area of the XMM-LSS field, roughly doubling the known cluster sample at this frontier redshift if confirmed. We construct a photometric redshift catalog based on deep ($i\sim26$, $K_\mathrm{s}\sim24$) multi-wavelength photometry from $u^*$-band to $K$-band gathered from the Hyper Suprime-Cam Subaru Strategic Program and other collaborative/public surveys. We adopt a Gaussian kernel density estimate with two different spatial scales (10" and 60") to draw a density map of massive ($\log(M_{*}/M_{\odot})>10.5$) and quiescent ($\log(\mathrm{sSFR\, [\mathrm{yr^{-1}}]})<-10$) galaxies at $z\sim2$. Then, We identify seven prominent overdensities. These candidates show clear red sequences in color-magnitude diagrams ($z-H$ vs. $H$). Moreover, one of them shows an extended X-ray emission with $L_\mathrm{X}=(1.46\pm0.35)\times10^{44}$ erg s$^{-1}$, suggesting its virialized nature. There is no clear evidence of enhancement nor suppression of the star formation rate of the main sequence galaxies in the clusters. We find that cluster galaxies have a higher fraction of transition population with $-10.5<\log(\mathrm{sSFR\, [\mathrm{yr^{-1}}]})<-10$ ($12\%$) than the field ($2\%$), which implies the ongoing star formation quenching. The quiescent fraction in the cluster candidates also exceeds that in the field. We confirm that the excess of a quiescent fraction is larger for higher-mass galaxies. This is the first statistical evidence for the mass-dependent environmental quenching at work in clusters even at $z\sim2$.

We present the first fluorine measurements in 12 normal giants belonging to the Galactic thin and thick disks using spectra obtained with the Phoenix infrared spectrometer on the 2.1m telescope at Kitt Peak. Abundances are determined from the (1-0) R9 2.3358 micron feature of the molecule HF. Additionally, sodium abundances are derived in 25 giants in the thin disk, thick disk, and halo using the Na I line at 2.3379 microns. We report fluorine abundances for thin and thick disk stars in the metallicity range -0.7 < [Fe/H] < 0. We add two abundance measurements for stars with [Fe/H] < 0.5 dex which are at a critical metallicity range to constrain models. We find a larger dispersion in fluorine abundances than sodium abundances despite both species having similar overall uncertainties due to atmospheric parameters, suggesting this dispersion is real and not observational. The dispersion is slightly larger in the thick disk than the thin. The thin and thick disk average [F/Fe] for our sample of stars combined with the literature differ by 0.03 dex. The observations are compared to available chemical evolution models.

We study the importance of precise modelling of the photometric redshift error distributions on the estimation of parameters from cross correlation measurements and present a working example of the scattering matrix formalism to correct for the redshift bin mismatch of objects in tomographic cross correlation analysis. We measured the angular galaxy auto-power spectrum and cross-power spectrum in four tomographic bins with redshift intervals $z = [0.0,0.3,0.45,0.6,0.8]$ from the cross correlation of Planck Cosmic Microwave Background lensing potential and photometric galaxy catalogue from the Dark Energy Spectroscopic Instrument Legacy Imaging Survey Data Release 8. We estimated galaxy linear bias and amplitude of cross correlation using maximum likelihood estimation to put constraints on the $\sigma_{8}$ parameter. We show that the modified Lorentzian function used to fit the photometric redshift error distribution performs well only near the peaks of the distribution. We adopt a sum of Gaussians model to capture the broad tails of the error distribution. Our sum of Gaussians model yields $\sim 2-5\,\sigma$ smaller values of cross correlation amplitude compared to the $\Lambda$CDM expectations. We compute the $\sigma_{8}$ parameter after correcting for the redshift bin mismatch of objects following the scattering matrix approach. The $\sigma_{8}$ parameter becomes consistent with $\Lambda$CDM model in the last tomographic bin but shows $\sim 1-3\,\sigma$ tension in other redshift bins.

The surface conditions of terrestrial bodies strongly reflect their geological evolutionary processes and vary among various terrestrial bodies. This diversity is attributed to variations in the timescales of boulder formation through processes such as impact cratering, rockfalls from crater walls, seismic motion, and boulder fragmentation caused by micrometeoroid impacts and thermal stress. In this study, we examined boulders on Ceres using high-resolution images with a resolution of approximately 5 m/px obtained during the Ceres Extended Mission 2 Orbit 7 of the Dawn mission. Almost all boulders were present around impact craters, even at a resolution of 5 m/px, thus indicating that the boulders on Ceres were created by impact cratering alone. The maximum boulder size on Ceres is approximately 200 m, even around large craters, which may indicate the upper size limit determined by the mechanical strength of the boulders, such as the tensile strength, the scale effect, and/or shattering strength. The slope of the size-frequency distribution of boulders on Ceres varied significantly across the range of boulder sizes, thus making it difficult to describe it using a single function of a power-law relationship; in particular, it changed at approximately 100 m, thus indicating that destructive or formation mechanisms may be different for large boulders > 100 m and for small boulders < 100 m. There may also be a subsurface structure that prevents the formation of small boulders, although this is difficult to argue conclusively. We estimated that the lifetime of boulders larger than 50 m was equivalent to or shorter than 100 Myr. This lifetime is consistent with a theoretical estimation assuming that micrometeoroid impacts are the primary destructive mechanism.

The possibility of slow diffusion regions as the origin for extended TeV emission halos around some pulsars (such as PSR J0633+1746 and PSR B0656+14) challenges the standard scaling of the electron diffusion coefficient in the interstellar medium. Self-generated turbulence by electron-positron pairs streaming out of the pulsar wind nebula was proposed as a possible mechanism to produce the enhanced turbulence required to explain the morphology and brightness of these TeV halos. We perform fully kinetic 1D3V particle-in-cell simulations of this instability, considering the case where streaming electrons and positrons have the same density. This implies purely resonant instability as the beam does not carry any current. We compare the linear phase of the instability with analytical theory and find very reasonable agreement. The non-linear phase of the instability is also studied, which reveals that the intensity of saturated waves is consistent with a momentum exchange criterion between a decelerating beam and growing magnetic waves. With the adopted parameters, the instability-driven wavemodes cover both the Alfvénic (fluid) and kinetic scales. The spectrum of the produced waves is non-symmetric, with left-handed circular polarisation waves being strongly damped when entering the ion-cyclotron branch, while right-handed waves are suppressed at smaller wavelength when entering the Whistler branch. The low-wavenumber part of the spectrum remains symmetric when in the Alfvénic branch. As a result, positrons behave dynamically differently compared to electrons. We also observed a second harmonic plasma emission in the wave spectrum. An MHD-PIC approach is warranted to probe hotter beams and investigate the Alfvén branch physics. This work confirms that the self-confinement scenario develops essentially according to analytical expectations [...](abridged)

Coronal mass ejections (CMEs) are major drivers of geomagnetic storms, which may cause severe space weather effects. Automating the detection, tracking, and three-dimensional (3D) reconstruction of CMEs is important for operational predictions of CME arrivals. The COR1 coronagraphs on board the Solar Terrestrial Relations Observatory spacecraft have facilitated extensive polarization observations, which are very suitable for the establishment of a 3D CME system. We have developed such a 3D system comprising four modules: classification, segmentation, tracking, and 3D reconstructions. We generalize our previously pretrained classification model to classify COR1 coronagraph images. Subsequently, as there are no publicly available CME segmentation data sets, we manually annotate the structural regions of CMEs using Large Angle and Spectrometric Coronagraph C2 observations. Leveraging transformer-based models, we achieve state-of-the-art results in CME segmentation. Furthermore, we improve the tracking algorithm to solve the difficult separation task of multiple CMEs. In the final module, tracking results, combined with the polarization ratio technique are used to develop the first single-view 3D CME catalog without requiring manual mask annotation. Our method provides higher precision in automatic 2D CME catalog and more reliable physical parameters of CMEs, including 3D propagation direction and speed. The aforementioned 3D CME system can be applied to any coronagraph data with the capability of polarization measurements.

Prominences and coronal rain are two forms of coronal condensations for which we still lack satisfactory details on the formation pathways and conditions under which the two come to exist. We compared prominences that formed via a steady versus stochastic type of heating. We performed 2.5D simulations using the open-source MPI-AMRVAC code. To further extend the work and allow for future direct comparison with observations, we used Lightweaver to form spectra of the filament view of our steady case prominence. With that, we analysed a reconnection event that shares certain characteristics with nanojets. We show how different forms of localised heating that induce thermal instability result in prominences with different properties. The steady form of heating results in prominence with a clear vertical structure stretching across the magnetic field lines. On the other hand, stochastic heating produces many threads that predominantly have a horizontal motion along the field lines. In the steady heating case, the prominence is relatively static; however, there is evidence of reconnection happening almost the entire time the prominence is present. In the case of stochastic heating, the threads are highly dynamic, with them also exhibiting a form of transverse oscillation (strongly resembling the decayless type). The fact that the threads in the stochastic heating case are constantly moving along the field lines suppresses any conditions for reconnection. It, therefore, appears that, to first order, the choice of heating prescription defines whether the prominence-internal dynamics are oriented vertically or horizontally. We closely inspected a sample reconnection event and computed the synthetic optically thick radiation using the open-source Lightweaver radiative transfer framework.

We establish a cosmological-model-independent method to extract the apparent magnitude and its derivative at different redshifts from the Pantheon+ type Ia supernova sample, and find that the obtained values deviate clearly from the prediction of the $\Lambda$CDM model at the lowest redshift. This deviation can be explained as a result of a transition of the absolute magnitude $M$ in the low redshift region. The observations seem to favor this transition since the minimum values of $\chi^2$ for two ansatzes of a varying $M$ are less than that of a constant $M$. The Hubble constant tension is alleviated from larger than $5\sigma$ to be about $1$ to $2\sigma$ for a varying $M$, and the growth tension can be resolved after attributing the variation of $M$ to a modification of the effective Newton's constant.

Recent high-resolution and sensitivity ALMA observations have unveiled the carbon isotope ratios ($^{12}$C/$^{13}$C) of Complex Organic Molecules (COMs) in a low-mass protostellar source. To understand the $^{12}$C/$^{13}$C ratios of COMs, we investigated the carbon isotope fractionation of COMs from prestellar cores to protostellar cores with a gas-grain chemical network model. We confirmed that the $^{12}$C/$^{13}$C ratios of small molecules are bimodal in the prestellar phase: CO and species formed from CO (e.g., CH$_{3}$OH) are slightly enriched in $^{13}$C compared to the local ISM (by $\sim$ 10 $\%$), while those from C and C$^{+}$ are depleted in $^{13}$C owing to isotope exchange reactions. COMs are mainly formed on the grain surface and in the hot gas ($>$ 100 K) in the protostellar phase. The $^{12}$C/$^{13}$C ratios of COMs depend on which molecules the COMs are formed from. In our base model, some COMs in the hot gas are depleted in $^{13}$C compared to the observations. Thus, We additionally incorporate reactions between gaseous atomic C and H$_{2}$O ice or CO ice on the grain surface to form H$_2$CO ice or \ce{C2O} ice, as suggested by recent laboratory studies. The direct C-atom addition reactions open pathways to form \ce{^13C}-enriched COMs from atomic C and CO ice. We find that these direct C-atom addition reactions mitigate $^{13}$C-depletion of COMs, and the model with the direct C-atom addition reactions better reproduces the observations than our base model. We also discuss the impact of the cosmic ray ionization rate on the $^{12}$C/$^{13}$C ratio of COMs.

Black holes threaded by massive vector fields can be subject to a superradiant instability, growing a cloud of massive vector particles around it. In this work, we consider what happens if such a dark matter candidate field mimicking a dark photon interacts with an accretion flow onto the black hole. By including a kinetic mixing term with the standard model photon, we extend the commonly used equations of general-relativistic magnetohydrodynamics to a dark photon constituent. The coupling to the dark photon then appears as an effective dynamo term together with a dark Lorentz force acting on the accreting matter. We numerically study the interactions between the superradiant dark photon cloud and the inner accretion flow by solving the coupled system in full numerical relativity. By parameterically varying the mixing parameter between dark and standard model sector, we provide a first investigation of how the accretion flow could be modified. Depending on the coupling strength, our solutions exhibit increased wind launching, as well as oscillation modes in the disk.

The stochasticity in galaxy clustering, the mismatch between galaxy and underlying matter distribution, suppresses the matter clustering amplitude reconstructed by the combination of galaxy auto-correlation and galaxy-galaxy lensing cross-correlation. In this work, we solve the stochasticity systematics by parameterizing the cross correlation coefficient $r(k)$ between galaxy and matter. We investigate the performance of 12 kinds of parameterization schemes, against the cosmoDC2 $\&$ TNG300-1 galaxy samples over a wide range of redshift and flux cut. The 2-parameter fits are found to describe the stochasticity up to $k_{\rm max}=0.9\,{\rm Mpc^{-1}}h$, while the best performing quadratic scheme $r^2_s(k) = 1+c_1 k+c_2 k^2$ reaches better than $1\%$ accuracy for both the direct ${r}^2_s(k)$ fit and reconstructing matter clustering. Then, we apply the accurate quadratic scheme to forecast the tomographic matter clustering reconstruction by the combination DESI-like LRG $\times$ CSST-like cosmic shear. Depending on assumption of stochasticity, we find that the neglect of a serious stochasticity would result in significant systematic bias in both the reconstruction and the inferred cosmological parameters, even if we adopt scale cut $k_{\rm max}=0.1\,{\rm Mpc^{-1}}h$. We demonstrate the necessity of including stochasticity in reconstruction, and forecast that the reconstruction alone enables a $S_8$ constraint at about $1.5\%$ precision, free from galaxy bias and stochasticity. We will validate our method for DESI spectroscopic survey, and the analysis is expected to be complementary to DESI cosmological constraint by BAO and RSD.

Based on the slopes between DESI $g,r$ and IRAS 100 $\mu m$ intensities, specifically $k_{g}$ and $k_{r}$, we have constructed a substantial sample of Galactic cirri. This sample covers 561.25 deg$^2$ at high Galactic latitudes (|b| $\geq$ 30$^{\circ}$), allowing for a systematic study of the physical parameters of the Galactic cirrus on a large scale, such as $g-r$ color, dust temperature, asymmetry factor and albedo. The ratio of $k_{g}$ and $k_{r}$ is consistent with the diffuse Galactic starlight model, suggesting that the diffuse starlight within our own Galaxy serves as the primary illumination source for the cirrus. Both $k_{g}$ and $k_{r}$ decrease slowly with increasing Galactic latitudes and IRAS 100 $\mu m$ intensities, while they do not have a correlation with Galactic longitudes. The distribution of $k_{g}$ and $k_{r}$ confirms a significant scatter in the slopes, reaching a factor of 4-5. Such large scatter cannot be explained by the weak correlation between the slopes and Galactic latitudes and IRAS 100 $\mu m$ intensities. Instead, it is attributed to substantial variations in the intrinsic properties of the dust, e.g., asymmetry factor and albedo. We propose that the properties of dust particles play a critical role in the observed scatter in slopes, making them the primary contributing factors. Moreover, the variations in dust properties within the cirrus are localized rather than exhibiting large-scale gradients.

The development of detectors with a high time resolution has been pivotal to our comprehension of neutron stars and the accurate measurement of their properties. While high-time resolution astronomy has become a standard in the radio and the high-/very-high-energy bands, progress in the visible band has been comparatively much slower. SiFAP2 is a high-speed optical photometer mounted at the INAF Telescopio Nazionale Galileo. Its potential emerged with the discovery of the first two optical millisecond pulsars: these are among the most efficient particle accelerators and natural laboratories of fundamental physics. Optical millisecond pulsations challenge the standard pulsar paradigm, requiring innovative solutions. Higher photon counting statistics of optical telescopes, compared to high-energy instruments, attain unprecedented sensitivity for weak pulsed signals from bright accreting neutron stars, which are the best candidates for still undetected continuous gravitational waves.

The Laser Interferometer Space Antenna (LISA) is a gravitational wave detector in space. It relies on a post-processing technique named time-delay interferometry (TDI) to suppress the overwhelming laser frequency noise by several orders of magnitude. This algorithm requires intersatellite-ranging monitors to provide information on spacecraft separations. To fulfill this requirement, we will use on-ground observatories, optical sideband-sideband beatnotes, pseudo-random noise ranging (PRNR), and time-delay interferometric ranging (TDIR). This article reports on the experimental end-to-end demonstration of a hexagonal optical testbed used to extract absolute ranges via the optical sidebands, PRNR, and TDIR. These were applied for clock synchronization of optical beatnote signals sampled at independent phasemeters. We set up two possible PRNR processing schemes: Scheme 1 extracts pseudoranges from PRNR via a calibration relying on TDIR; Scheme 2 synchronizes all beatnote signals without TDIR calibration. The schemes rely on newly implemented monitors of local PRNR biases. After the necessary PRNR treatments (unwrapping, ambiguity resolution, bias correction, in-band jitter reduction, and/or calibration), Scheme 1 and 2 achieved ranging accuracies of 2.0 cm to 8.1 cm and 5.8 cm to 41.1 cm, respectively, below the classical 1 m mark with margins.

We study the propagation properties of slow magneto-acoustic waves in a multi-thermal coronal loop using a 3D MHD model, for the first time. A bundle of 33 vertical cylinders, each of 100{\,}km radius, randomly distributed over a circular region of radius 1{\,}Mm is considered to represent the coronal loop. The slow waves are driven by perturbing the vertical velocity ($v_z$) at the base of the loop. We apply forward modelling to the simulation results to generate synthetic images in the coronal channels of SDO/AIA. Furthermore, we add appropriate data noise to enable direct comparison with the real observations. It is found that the synthetic images at the instrument resolution show non-cospatial features in different temperature channels in agreement with previous observations. Time-distance maps are constructed from the synthetic data to study the propagation properties. The results indicate that the oscillations are only visible in specific channels depending on the temperature range of plasma existing within the loop. Additionally, the propagation speed of slow waves is also found to be sensitive to the available temperature range. Overall, we propose that the cross-field thermal properties of coronal structures can be inferred using a combination of numerical simulations and observations of slow magneto-acoustic waves.

Planets orbiting binary systems are relatively unexplored compared to those around single stars. Detections of circumbinary planets and planetary systems offer a first detailed view into our understanding of circumbinary planet formation and dynamical evolution. The BEBOP (Binaries Escorted by Orbiting Planets) radial velocity survey plays a special role in this adventure as it focuses on eclipsing single-lined binaries with an FGK dwarf primary and M dwarf secondary allowing for the highest-radial velocity precision using the HARPS and SOPHIE spectrographs. We obtained 4512 high-resolution spectra for the 179 targets in the BEBOP survey which we used to derive the stellar atmospheric parameters using both equivalent widths and spectral synthesis. We furthermore derive stellar masses, radii, and ages for all targets. With this work, we present the first homogeneous catalogue of precise stellar parameters for these eclipsing single-lined binaries.

Primordial black holes (PBHs) are an attractive dark matter candidate, particularly if they can explain the totality of it. At PBH masses below $\sim 10^{17}$g and above $\sim 10^{23}$g this possibility is excluded from the variety of arguments and with different confidence. The range in between, often referred to as the "asteroid mass window", currently remains unconstrained. The most promising, in our view, way to probe this mass range is to use stars as the PBH detectors. If a star captures even a single PBH it starts being accreted onto it and eventually gets destroyed -- converted into a sub-solar mass black hole. This process may have a variety of signatures form a mere star disappearance to supernova-type explosions of a new kind. The viability of this approach depends crucially on the probability of PBH capture by stars. In this chapter we summarize the existing capture mechanisms and discuss their implications for constraining the abundance of (or perhaps discovering) PBHs in the asteroid mass window.

Herschel operated as an observatory, therefore it did not cover the whole sky, but still observed ~8% of it. The first version of an overall Herschel/PACS Point Source Catalogue was released in 2017. The data are still unique and are very important for research, especially because no new far-infrared mission is foreseen for at least the next decade. In the framework of the NEMESIS project, we revisited all the photometric observations obtained by the PACS instrument on-board of the Herschel space observatory. We aimed to build the most complete and most accurate Herschel/PACS catalogue to date. Our primary goal was to increase the number of real sources, and decrease the number of spurious sources identified on a strongly variable background. Our goal was to build a blind catalogue, meaning that source extraction is conducted without relying on prior detections at various wavelengths, allowing us to detect sources never catalogued before. We define a hybrid strategy that includes classical and ML source identification and characterisation methods, providing catalogues at much higher completeness levels than before. Quality assessment also involves ML techniques. Our source extraction methodology facilitates a systematic and impartial comparison of sensitivity levels across various Herschel fields, a task that was typically beyond the scope of individual programs. We created a high-reliability and a rejected source catalogue for each PACS passband, i.e., at 70, 100, and 160 {\mu}m. With the high-reliability catalogue, we managed to significantly increase the completeness in all bands. At the same time, while the number of high-reliability detections decreased, the number of sources matching with existing catalogues increased, suggesting that the purity is also higher than before. The photometric accuracy of our pipeline is ~1% based on comparison with the standard star models.

The morphology of the 21-cm signal emitted by the neutral hydrogen present in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) depends both on the properties of the sources of ionizing radiation and on the underlying physical processes within the IGM. Variation in the morphology of the IGM 21-cm signal due to the different sources of the EoR is expected to have a significant impact on the 21-cm bispectrum, which is one of the crucial observable statistics that can evaluate the non-Gaussianity present in the signal and which can be estimated from radio interferometric observations of the EoR. Here we present the 21-cm bispectrum for different reionization scenarios assuming different simulated models for the sources of reionization. We also demonstrate how well the 21-cm bispectrum can distinguish between different IGM 21-cm signal morphologies, arising due to the differences in the reionization scenarios, which will help us shed light on the nature of the sources of ionizing photons. Our estimated large-scale bispectrum for all unique $k$-triangle shapes shows a significant difference in their magnitude and sign across different reionization scenarios. Additionally, our focused analysis of bispectrum for a few specific $k$-triangle shapes (e.g. squeezed-limit, linear, and shapes in the vicinity of the squeezed-limit) shows that the large scale 21-cm bispectrum can distinguish between reionization scenarios that show inside-out, outside-in and a combination of inside-out and outside-in morphologies. These results highlight the potential of using the 21-cm bispectrum for constraining different reionization scenarios.

High-redshift blazars provide valuable input to studies of the evolution of active galactic nuclei (AGN) jets and provide constraints on cosmological models. Detections at high energies ($0.1<\mathrm{E}<100$ GeV) of these distant sources are rare, but when they exhibit bright gamma-ray flares, we are able to study them. However, contemporaneous multi-wavelength observations of high-redshift objects ($z>4$) during their different periods of activity have not been carried out so far. An excellent opportunity for such a study arose when the blazar TXS 1508+572 ($z=4.31$) exhibited a $\gamma$-ray flare in 2022 February in the $0.1-300$ GeV range with a flux 25 times brighter than the one reported in the in the fourth catalog of the \textit{Fermi} Large Area Telescope. Our goal is to monitor the morphological changes, spectral index and opacity variations that could be associated with the preceding $\gamma$-ray flare in TXS 1508+572 to find the origin of the high-energy emission in this source. We also plan to compare the source characteristics in the radio band to the blazars in the local Universe ($z<0.1$). In addition, we aim to collect quasi-simultaneous data to our multi-wavelength observations of the object, making TXS 1508+572 the first blazar in the early Universe ($z>4$) with contemporaneous multi-frequency data available in its high state. In order to study the parsec-scale structure of the source, we performed three epochs of very-long-baseline interferometry (VLBI) follow-up observations with the Very Long Baseline Array (VLBA) supplemented with the Effelsberg 100-m Telescope at 15, 22, and 43 GHz, which corresponds to 80, 117, and 228 GHz in the rest frame of TXS 1508+572. In addition, one 86 GHz (456 GHz) measurement was performed by the VLBA and the Green Bank Telescope during the first epoch.

White dwarfs are the remnants of stars not massive enough to become supernovae. This review explores the concept of strange dwarfs, a unique class of white dwarfs which contain cores of strange quark matter. Strange dwarfs have different sizes, masses, and evolutionary paths with respect to white dwarfs. They might form through the accumulation of normal matter on strange quark stars or by capture of strangelets. The stability of strange dwarfs has been debated, with initial studies suggesting stability, while later analyses indicated potential instability. This review revisits these discussions, focusing on the critical role of boundary conditions between nuclear and quark matter in determining stability. It also offers insights into their formation, structure, and possible detection in the universe.

Galaxy clusters provide an avenue to expand our knowledge of cosmology and galaxy evolution. Because it is difficult to accurately measure the total mass of a large number of individual clusters, cluster samples are typically selected using an observable proxy for mass. Selection effects are therefore a key problem in understanding galaxy cluster statistics. We make use of the $(2.8~\rm{Gpc})^3$ FLAMINGO hydrodynamical simulation to investigate how selection based on X-ray luminosity, thermal Sunyaev-Zeldovich effect or galaxy richness influences the halo mass distribution. We define our selection cuts based on the median value of the observable at a fixed mass and compare the resulting samples to a mass-selected sample. We find that all samples are skewed towards lower mass haloes. For X-ray luminosity and richness cuts below a critical value, scatter dominates over the trend with mass and the median mass becomes biased increasingly low with respect to a mass-selected sample. At $z\leq0.5$, observable cuts corresponding to median halo masses between $M_\text{500c}=10^{14}$ and $10^{15}~\rm{M_{\odot}}$ give nearly unbiased median masses for all selection methods, but X-ray selection results in biased medians for higher masses. For cuts corresponding to median masses $<10^{14}$ at $z\leq0.5$ and for all masses at $z\geq1$, only Compton-Y selection yields nearly unbiased median masses. Importantly, even when the median mass is unbiased, the scatter is not because for each selection the sample is skewed towards lower masses than a mass-selected sample. Each selection leads to a different bias in secondary quantities like cool-core fraction, temperature and gas fraction.

White dwarfs offer a unique opportunity to search nearby stellar systems for signs of life, but the habitable zone around these stars is still poorly understood. Since white dwarfs are compact stars with low luminosity, any planets in their habitable zone should be tidally locked, like planets around M-dwarfs. Unlike planets around M-dwarfs, however, habitable white dwarf planets have to rotate very rapidly, with orbital periods ranging from hours to several days. Here we use the ExoCAM Global Climate Model (GCM) to investigate the inner edge of the habitable zone (HZ) around white dwarfs. Our simulations show habitable planets with ultrashort orbital periods ($P\lesssim$1 day) enter a ``bat rotation" regime, which differs from typical atmospheric circulation regimes around M dwarfs. Bat rotators feature mean equatorial subrotation and a displacement of the surface's hottest regions from the equator towards the midlatitudes. We qualitatively explain the onset of bat rotation using shallow water theory. The resulting circulation shifts increase dayside cloud cover and decrease stratospheric water vapor, expanding the white dwarf habitable zone by $\sim$50\% compared to estimates based on 1D models. The James Webb Space Telescope (JWST) should be able to quickly characterize bat rotators around nearby white dwarfs thanks to their distinct thermal phase curves. Our work underlines that tidally locked planets on ultrashort orbits may exhibit unique atmospheric dynamics, and guides future habitability studies of white dwarf systems.

We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and sensitivity in the infrared, now reveal that there are more than ten discrete regions with outflows, nearly doubling the number of such known regions. Outflows exist near the radio lobes, at the nucleus, in a biconical structure perpendicular to the jet, and in a bubble moving against the disk. In some of them, velocities above escape velocity are observed. Stratification is also observed, with higher ionization or excitation gas attaining higher velocities. More outflows and bow shocks, found further away from the nucleus than the radio lobes, in regions without significant radio emission, reveal the existence of past or weak radio jets that interacted with the interstellar medium. The coincidence of the bow shocks with the optical extended emission line region (EELR) suggests that the jets also contributed to the gas ionization. Maps of the H2 gas excitation temperature, T_ex, indicate that the molecular gas is most excited in regions with radio emission. There, T_ex is more than 100K higher than in the EELR interior. We argue that a combination of jet-related shocks and cosmic rays is likely responsible for this excess molecular gas excitation.

Most of the gamma-ray sources in the Fermi-LAT 14-year Source Catalog are associated with pulsars and blazars. However, unveiling the nature of the still unassociated gamma-ray sources is important for the understanding of high energy emission mechanisms in astrophysical objects. This work presents a comprehensive study toward the region covered by the Fermi source 4FGL J1846.9$-$0227 previously suggested to be a blazar or a massive protostar. Using multiwavelength observations, we analyzed several astrophysical objects in the region as possible counterparts of the Fermi-LAT source. Having discarded most of them after a comprehensive analysis, we suggest that the most likely candidate to be such a counterpart is IRAS 18443$-$0231. We discovered that this source, previously cataloged as a planetary nebula candidate, actually is a proto-planetary nebula. The radio continuum image at 3 GHz associated with such a nebula allowed us to identify a jet-like structure. Additionally, we identified an associated red-shifted CO molecular outflow and a dense molecular clump in which the source is embedded. We obtained a radio spectral index of $-0.47 \pm 0.08$ for the source, indicating syncrothron emission due to accelerated particles. Thus, we suggest that processes such as proton-proton collisions and relativistic Bremsstrahlung are likely to occur. IRAS 18443$-$0231, lying almost at the center of the Fermi confidence ellipse and related to the hard X-ray source 4XMM J184700.4$-$022752, would be the first association between a proto-planeatry nebula and gamma emission.

The phenomenon of magnification bias can induce a non-negligible angular correlation between two samples of galaxies with nonoverlapping redshift distributions. This signal is particularly clear when background submillimeter galaxies are used, and has been shown to constitute an independent cosmological probe. This work extends prior studies on the submillimeter galaxy magnification bias to the massive neutrino scenario, with the aim being to assess its sensitivity as a cosmological observable to the sum of neutrino masses. The measurements of the angular cross-correlation function between moderate redshift GAMA galaxies and high-redshift submillimeter H-ATLAS galaxies are fit to the weak lensing prediction down to the arcmin scale. The signal is interpreted under the halo model, which is modified to accommodate massive neutrinos. We discuss the impact of the choice of cosmological parametrization on the sensitivity to neutrino masses. The currently available data on the magnification bias affecting submillimeter galaxies are sensitive to neutrino masses when a cosmological parametrization in terms of the primordial amplitude of the power spectrum $(A_s$) is chosen over the local root mean square of smoothed linear density perturbations $(\sigma_8$). A clear upper limit on the sum of neutrino masses can be derived if the value of $A_s$ is either fixed or assigned a narrow Gaussian prior, a behavior that is robust against changes to the chosen value.

We present a study of $1347$ galaxies at $z<0.35$ with detected nuclear X-ray emission and optical emission line diagnostics in the Baldwin-Phillips-Terlevich (BPT) diagram. This sample was obtained by cross-matching the X-ray Multi-Mirror Mission Observatory - Newton (XMM-Newton) DR10 catalogue with Sloan Digital Sky Survey (SDSS) DR17 galaxies with well-measured line ratios. The distribution of these sources in the BPT diagram covers all three excitation regimes: Ionized Hydrogen (HII) regions (23\%), `composites' (30\%), and Seyfert galaxies with the low ionization nuclear emission line regions (LINERs) (47\%). In contrast, the fraction of objects classified as active galactic nuclei (AGN) in the SDSS subsample selected for cross-match with XMM-Newton is only 13\%. This fact illustrates that X-ray emission from galaxies commonly points towards the presence of AGN. Our data show, for the first time, a clear dependence of the BPT position on the ratio of the X-ray to $H\alpha$ fluxes. Sources dominated by X-ray emission lie in the Seyfert and LINER regimes of the BPT diagram. Most sources with a low X-ray-to-$H\alpha$-luminosity ratio, $log_{10}(L_X/L_{H\alpha}) < 1.0$, lie in the HII regime. In our sample, there are even 45 galaxies that have $L^{Star}_{XR}/L^{Total}_{Xray}>0.5$. In contrast, the positions of the sample members in the BPT diagram exhibit {no} dependence on the X-ray hardness ratio. Our finding suggests that the X-ray-to-$H\alpha$ ratio can help us to differentiate galaxies whose X-ray flux is dominated by an AGN {from galaxies with} central X-ray binaries and other stellar X-ray sources.

The spins of binary black holes (BBHs) measured from gravitational waves carry notable information of the formation pathways. Here we propose a quantity "dimensionless net spin" ($\chi_{\rm N}$), which is related to the sum of angular momentum of component black holes in the system, to provide a novel perspective to study the origin(s) of BBHs. By performing hierarchical Bayesian inference on $\chi_{\rm N}$, we find strong evidence that the marginal distribution of this quantity can be better fitted by two Gaussian components rather than one: there is a narrow peak at $\chi_{\rm N} \sim 0.15$ and another extended peak at $\chi_{\rm N} \sim 0.47$. We also find that the rapidly spinning systems likely dominate the high-mass end of the population and they evolve with redshift much quicker. These findings bring new challenges to the field binary scenario, and suggest that dynamical process should plays a key role in forming high total mass BBHs.

For the first time, a calibration between the HeI $\lambda5876$/H$\beta$ emission line ratio and the helium abundance $y$=12+log(He/H) for Narrow line regions (NLRs) of Seyfert~2 Active Galactic Nuclei (AGN) is proposed. In this context, observational data (taken from the SDSS-DR15 and from the literature) and direct abundance estimates (via the $T_{\rm e}$-method) for a sample of 65 local ($z \: < \: 0.2$) Seyfert~2 nuclei are considered. The resulting calibration estimates the $y$ abundance with an average uncertainty of 0.02 dex. Applying our calibration to spectroscopic data containing only strong emission lines, it yields a helium abundance distribution similar to that obtained via the $T_{\rm e}$-method. Some cautions must be considered to apply our calibration for Seyfert~2 nuclei with high values of electron temperature ($\gtrsim\: 20\,000$ K) or ionization parameter ($\log U > -2.0$).

Star clusters are among the first celestial objects catalogued by early astronomers. As simple and coeval populations, their study has been instrumental in charting the properties of the Milky Way and providing insight into stellar evolution through the 20th century. Clusters were traditionally spotted as local stellar overdensities in the plane of the sky. In recent decades, for a limited number of nearby clusters, it became possible to identify cluster members through their clustering in proper motion space. With its astrometric data of unprecedented precision, the Gaia mission has completely revolutionised our ability to discover and characterise Milky Way star clusters, to map their large-scale distribution, and to investigate their internal structure. In this review we focus on the population of open clusters, residing in the Galactic disc. We summarise the current state of the Gaia-updated cluster census and studies of young clusters and associations. We discuss recent developments in techniques for cluster detection and age estimation. We also review results enabled by Gaia data concerning the dynamical evolution of gravitationally bound clusters and their stellar inventory.

This paper is focused on the segregation of FGK dwarf and giant stars through narrow-band photometric data using the Spanish Virtual Observatory (SVO) Filter Profile Service and associated photometric tools. We selected spectra from the MILES, STELIB, and ELODIE stellar libraries, and used SVO photometric tools to derive the synthetic photometry in 15 J-PAS narrow filters, which were especially selected to cover spectral features sensitive to gravity changes. Using machine-learning techniques as the Gaussian mixture model and the support vector machine, we defined several criteria based on J-PAS colours to discriminate between dwarf and giant stars. We selected five colour-colour diagrams that presented the most promising separation between both samples. Our results show an overall accuracy in the studied sample of $\sim$0.97 for FGK stars, although a dependence on the luminosity type and the stellar effective temperature was found. We also defined a colour-temperature relation for dwarf stars with effective temperatures between 4\,000 and 7\,000\,K, which allows one to estimate the stellar effective temperature from four J-PAS filters ($J0450$, $J0510$, $J0550$, and $J0620$). Additionally, we extended the study to M-type giant and dwarf stars, achieving a similar accuracy to that for FGK stars.

We report the detection of high-energy gamma-ray emission towards the G305 star-forming region. Using almost 15 years of observation data from {\sl Fermi} Large Area Telescope, we detected an extended gamma-ray source in this region with a significance of $\sim 13 \sigma$. The gamma-ray radiation reveals a clear pion-bump feature and can be fitted with the power law parent proton spectrum with an index of $-2.5$. The total cosmic ray (CR) proton energy in the gamma-ray production region is estimated to be the order of $10^{49}\ \rm erg$. We further derived the CR radial distribution from both the gamma-ray emission and gas distribution and found it roughly obeys the $1/r$ type profile, which is consistent with other similar systems and expected from the continuous injection of CRs by the central powerful young massive star cluster Danks 1 or Danks 2 in this region. Together with former detections of similar gamma-ray structures, such as Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and W40, the detection supports the hypothesis that young massive star clusters are CR accelerators.

The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos and explain how to interpret mass and radius measurements of NSs inferred from electromagnetic radiation if these dark halos exist.

We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment of galaxies via a shared luminosity function (LF). Considering this framework within the context of a Rubin Observatory's Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we first demonstrate the substantial impact of the LF on both source galaxy redshift distributions and the intrinsic alignment contamination. We establish how the individual parameters of a Schechter LF model impact the redshift distribution of a magnitude-limited sample, and we demonstrate the effect of marginalising over the LF parameters as incorporated in the intrinsic alignment modelling of a standard cosmic shear analysis set-up. We forecast the impact of our joint modelling framework on cosmological parameter constraints. Our preliminary results are promising, indicating that this framework can yield cosmological constraints consistent with those expected from standard analyses, enhanced by the flexibility of not fixing LF parameters. We plan to further validate these findings with comprehensive Markov chain Monte Carlo simulations to robustly quantify bias avoidance and underscore the framework's efficacy. Taking advantage of our forecasting results and the parameter degeneracies, we identify the specific impact of the shape of the LF of source galaxies on the cosmic shear data vector. We also discuss the potential of this method in providing a way to model generic selection functions in redshift distribution estimation, as well as its possibilities for extension to a 3x2pt analysis, particularly with respect to incorporating galaxy bias in this luminosity-function-based framework. Although we consider the context of LSST cosmic shear in this work, the proposed joint modelling framework is generically applicable to weak lensing surveys.

Stage IV dark energy surveys such as the Vera C. Rubin Observatory and Euclid present a unique opportunity to shed light on the nature of dark energy. However, the full constraining power of the new data cannot be unlocked unless accurate predictions are available at all observable scales. Currently, only the linear regime is well understood in models beyond $\Lambda$CDM: on the nonlinear scales, expensive numerical simulations become necessary, making their direct use impractical in the analyses of large datasets. Recently, machine learning techniques have shown potential to break this impasse: by training emulators, we can reproduce complex data fields in a fraction of the time it takes to produce them. In this work, we present a field-level emulator capable of turning a $\Lambda$CDM N-body simulation into one evolved under $f(R)$ gravity. To achieve this, we build on the map2map neural network, using the strength of modified gravity $|f_{R_0}|$ as style parameter. We find that our emulator correctly estimates the changes it needs to apply to the positions and velocities of the input N-body particles to produce the target simulation. We test the performance of our network against several summary statistics, finding $1\%$ agreement in the power spectrum up to $k \sim 1$ $h/$Mpc, and $1.5\%$ agreement against the independent boost emulator eMantis. Although the algorithm is trained on fixed cosmological parameters, we find it can extrapolate to models it was not trained on. Coupled with available field-level emulators and simulations suites for $\Lambda$CDM, our algorithm can be used to constrain modified gravity in the large-scale structure using full information available at the field-level.

Anti-de Sitter vacuum, which correspond to a negative cosmological constant (CC), is theoretically important and well-motivated. It is interesting to see whether current data can allow the existence of a negative CC not. In this paper, we perform the MCMC analysis for the $w_0w_a$CDM+CC model using recent DESI BAO measurements combined with Planck CMB and Pantheon Plus dataset. The results reveal that the bestfit energy density of CC is $\Omega_\Lambda\sim -0.3$ and the fitting to DESI is slightly improved, while $\Omega_\Lambda=0$ is also $1\sigma$ consistent.