Papers for Wednesday, Sep 18 2024

The star cluster surrounding the supermassive black hole in the center of Milky Way is probed using the data on the S2 star. The value of precession found at the physics-informed neural networks (PINN) analysis of the S2 data is used to consider the role of the scattering of S2 star on stars of the cluster, described by a random force given by the Holtsmark distribution. The critical value for the star density of the core cluster for which the observed precession value by PINN lies inside 70\% confidence interval (between 15\% and 85\% quantiles) around the median of precession due to scattering, is obtained as n_crit \approx 8.3 10^6 pc^-3, that is at higher star densities the perturbation of the orbit of S2 would exceed the observed one.

The assembly of galaxy clusters is understood to be a hierarchical process with a continuous accretion of galaxies over time, which increases the cluster size and mass. Late-type galaxies that fall into clusters can undergo ram-pressure stripping, forming extended gas tails within which star formation can happen. The number, location, and tail orientations of such galaxies provide clues about the galaxy infall process, the assembly of the cluster over time, and the consequences of infall for galaxy evolution. Here, we utilise the $\sim$ 0.5 degree diameter circular field of the Ultraviolet Imaging Telescope to image six galaxy clusters at z < 0.06 that are known to contain jellyfish galaxies. We searched for stripping candidates in the ultraviolet images of these clusters, which revealed 54 candidates showing signs of unilateral extra-planar emission, due to ram-pressure stripping. Seven candidates had already been identified as likely stripping based on optical B-band imaging. We identified 47 new candidates through UV imaging. Spectroscopic redshift information is available for 39 of these candidate galaxies, of which 19 are associated with six clusters. The galaxies with spectroscopic redshifts that are not part of the clusters appear to be within structures at different redshifts identified as additional peaks in the redshift distribution of galaxies, indicating that they might be ram-pressure stripped or disturbed galaxies in other structures along the line of sight. We examine the orbital history of these galaxies based on their location in the position-velocity phase-space diagram and explore a possible connection to the orientation of the tail direction among cluster member candidates. The tails of confirmed cluster member galaxies are found to be oriented away from the cluster centre.

Primordial gravitational waves could be non-Gaussian, just like primordial scalar perturbations. Although the tensor two-point function has thus-far remained elusive, the three-point function could, in principle, be large enough to be detected in Cosmic Microwave Background temperature and polarization anisotropies. We perform a detailed analysis of tensor and mixed tensor-scalar non-Gaussianity through the Planck PR4 bispectrum, placing constraints on eleven primordial templates, spanning various phenomenological and physical regimes including modifications to gravity, additional fields in inflation, and primordial magnetic fields. All analysis is performed using modern quasi-optimal binned estimators, and yields no evidence for tensor non-Gaussianity, with a maximum detection significance of $1.8\sigma$. Our constraints are derived primarily from large-scales (except for tensor-scalar-scalar models), and benefit greatly from the inclusion of $B$-modes. Although we find some loss of information from binning, mask effects and residual foreground contamination, our $f_{\rm NL}$ bounds improve over those of previous analyses by $40-600\%$, with six of the eleven models being analyzed for the first time. Unlike for scalar non-Gaussianity, future low-noise experiments such as LiteBIRD, the Simons Observatory and CMB-S4, will yield considerable improvement in tensor non-Gaussianity constraints.

Type IIn supernovae (SNeIIn) are a highly heterogeneous subclass of core-collapse supernovae, spectroscopically characterized by signatures of interaction with a dense circumstellar medium (CSM). Here we systematically model the light curves of 142 archival SNeIIn using MOSFiT (the Modular Open Source Fitter for Transients). We find that the observed and inferred properties of SNIIn are diverse, but there are some trends. The typical SN CSM is dense ($\sim$10$^{-12}$gcm$^{-3}$) with highly diverse CSM geometry, with a median CSM mass of $\sim$1M$_\odot$. The ejecta are typically massive ($\gtrsim10$M$_\odot$), suggesting massive progenitor systems. We find positive correlations between the CSM mass and the rise and fall times of SNeIIn. Furthermore there are positive correlations between the rise time and fall times and the $r$-band luminosity. We estimate the mass-loss rates of our sample (where spectroscopy is available) and find a high median mass-loss rate of $\sim$10$^{-2}$M$_\odot$yr$^{-1}$, with a range between 10$^{-4}$--1M$_\odot$yr$^{-1}$. These mass-loss rates are most similar to the mass loss from great eruptions of luminous blue variables, consistent with the direct progenitor detections in the literature. We also discuss the role that binary interactions may play, concluding that at least some of our SNeIIn may be from massive binary systems. Finally, we estimate a detection rate of 1.6$\times$10$^5$yr$^{-1}$ in the upcoming Legacy Survey of Space and Time at the Vera C. Rubin Observatory.

We analyze the evolution of the radial profiles and the azimuthal variations of the stellar metallicities from the Vintergatan simulation of a Milky Way-like galaxy. We find that negative gradients exist as soon as the disk settles at high redshift, and are maintained throughout the long term evolution of the galaxy, including during major merger events. The inside-out growth of the disk and an overall outward radial migration tend to flatten these gradients in time. Major merger events only have a moderate and short-lived imprint on the [Fe/H] distributions with almost no radial dependence. The reason lies in the timescale for enrichment in Fe being significantly longer than the duration of the starbursts episodes, themselves slower than dynamical mixing during typical interactions. It results that signatures of major mergers become undetectable in [Fe/H] only a few Myr after pericenter passages. We note that considering other tracers like the warm interstellar medium, or monitoring the evolution of the metallicity gradient as a single value instead of a radial full profile could lead to different interpretations, and warn against an oversimplification of this complex problem.

Astrophysical explosions that contain dense and ram-pressure-dominated ejecta evolve through an interaction phase, during which a forward shock (FS), contact discontinuity (CD), and reverse shock (RS) form and expand with time. We describe new self-similar solutions that apply to this phase and are most accurate in the limit that the ejecta density is large compared to the ambient density. These solutions predict that the FS, CD, and RS expand at different rates in time and not as single temporal power-laws, are valid for explosions driven by steady winds and homologously expanding ejecta, and exist when the ambient density profile is a power-law with power-law index shallower than $\sim 3$ (specifically when the FS does not accelerate). We find excellent agreement between the predictions of these solutions and hydrodynamical simulations, both for the temporal behavior of the discontinuities and for the variation of the fluid quantities. The self-similar solutions are applicable to a wide range of astrophysical phenomena and -- although the details are described in future work -- can be generalized to incorporate relativistic speeds with arbitrary Lorentz factors. We suggest that these solutions accurately interpolate between the initial ``coasting'' phase of the explosion and the later, energy-conserving phase (or, if the ejecta is homologous and the density profile is sufficiently steep, the self-similar phase described in Chevalier 1982),

The hydrogen Lyman-alpha (Lya) emission line, the brightest spectral feature of a photoionized gas, is considered an indirect tracer of the escape of Lyman continuum (LyC) photons, particularly when the intergalactic medium is too opaque for direct detection. However, resonant scattering complicates interpreting the empirical properties of Lya photons, necessitating radiative transfer simulations to capture their strong coupling with underlying gas kinematics. In this study, we leverage the exceptional spatial resolution from strong gravitational lensing to investigate the connection between Lya line profiles and LyC leakage on scales of a few 100 pc in the Sunburst Arc galaxy at $z\sim2.37$. New optical echelle spectra obtained using Magellan MIKE show that both the LyC leaking and non-leaking regions exhibit a classic double-peak Lya feature with an enhanced red peak, indicating outflows at multiple locations in the galaxy. Both regions also show a central Gaussian peak atop the double peaks, indicating directly escaped Lya photons independent of LyC leakage. We introduce a machine learning-based method for emulating Lya simulations to quantify intrinsic dynamics ($\sigma_{\mathrm{int}}$), neutral hydrogen column density ($N_{\mathrm{HI}}$), outflow velocity ($v_{\mathrm{exp}}$), and effective temperature ($T$) across continuous parameter spaces. By comparing the spatially and spectrally resolved Lya lines in Sunburst, we argue that the directly escaped Lya photons originate in a volume-filling, warm ionized medium spanning $\sim1$ kpc, while the LyC leakage is confined to regions of $\lesssim200$ pc. These sub-kpc variations in Lya profiles highlight the complexity of interpreting integrated properties in the presence of inhomogeneous mixtures of gas and young stars, emphasizing the need for spatially and spectrally resolved observations of distant galaxies.

The diffuse stellar component of galaxy clusters made up of intergalactic stars is termed the intracluster light (ICL). Though there is a developing understanding of the mechanisms by which the ICL is formed, no strong consensus has yet been reached on which objects the stars of the ICL are primarily sourced from. We investigate the assembly of the ICL starting approximately $10$ Gyr before $z=0$ in 11 galaxy clusters (halo masses between $\sim1\times 10^{14}$ M$_{\odot}$ and $\sim7\times 10^{14}$ M$_{\odot}$ at $z\approx0$) in the Horizon-AGN simulation. By tracking the stars of galaxies that fall into these clusters past cluster infall, we are able to link almost all of the $z\approx0$ ICL back to progenitor objects. Satellite stripping, mergers, and pre-processing are all found to make significant contributions to the ICL, but any contribution from in-situ star-formation directly into the ICL appears negligible. Even after compensating for resolution effects, we find that approximately $90$ per cent of the stacked ICL of the 11 clusters that is not pre-processed should come from galaxies infalling with stellar masses above $10^{9}$ M$_{\odot}$, with roughly half coming from infalling galaxies with stellar masses within half a dex of $10^{11}$ M$_{\odot}$. The fact that the ICL appears largely sourced from such massive objects suggests that the ICL assembly of any individual cluster may be principally stochastic.

We explore full-shape analysis with simulation-based priors, which is the simplest approach to galaxy clustering data analysis that combines effective field theory (EFT) on large scales and numerical simulations on small scales. The core ingredient of our approach is the prior density of EFT parameters which we extract from a suite of 10500 galaxy simulations based on the halo occupation distribution (HOD) model. We measure the EFT parameters with the field-level forward model, which enables us to cancel cosmic variance. On the theory side, we develop a new efficient approach to calculate field-level transfer functions using time-sliced perturbation theory and the logarithmic fast Fourier transform. We study cosmology dependence of EFT parameters of dark matter halos and HOD galaxies and find that it can be ignored for the purpose of prior generation. We use neural density estimation to model the measured distribution of EFT parameters. Our distribution model is then used as a prior in a reanalysis of the BOSS full-shape galaxy power spectrum data. Assuming the $\Lambda$CDM model, we find significant ($\approx 30\%$ and $\approx 60\%$) improvements for the matter density fraction and the mass fluctuation amplitude, which are constrained to $\Omega_m= 0.315 \pm 0.010$ and $\sigma_8 = 0.671 \pm 0.027$. The value of the Hubble constant does not change, $H_0= 68.7\pm 1.1$ km/s/Mpc. This reaffirms earlier reports of the structure growth tension from the BOSS data. Finally, we use the measured EFT parameters to constrain galaxy formation physics.

Recent observations by the James Webb Telescope (JWST) have unveiled numerous galaxy candidates between $z \sim 9 - 16.5$, hinting at an over-abundance of the bright-end of the UV Luminosity Function (UV LF) $z \gtrsim 11$. Possible solutions require extremely bursty star formation, these systems being dust-free, an evolving initial mass function or even cosmic variance. In this work, we develop an analytic formalism to study dust enrichment and its impact on the UV luminosity of both main-sequence early galaxies and extremely bursty star formers. Our dust model, including the key processes of dust production in type II Supernovae, dust destruction, ejection, growth and sputtering, is calibrated against the latest datasets from the Atacama Large Millimeter Array (ALMA) at $z \sim 4-7$. The model has only 3 free parameters: (i) the star formation efficiency; (ii) the dust growth timescale; and (iii) the dust distribution radius. Our key results are: (i) explaining the observed UV LF requires an average star formation efficiency that increases with redshift as $f_*(z) = 10^{0.13z-3.5}$ at $z \sim 5-13$ with a number of observations hinting at objects lying a factor 10 above this main-sequence. (ii) The dust enrichment of early systems is driven by dust production in SNII ejecta; growth and sputtering are the second and third most crucial processes, impacting the dust mass by 60% and 40% respectively at $z \sim 7$. (iii) In our model, galaxies at $z \gtrsim 9$ can still host significant amounts of dust reaching average dust-to-stellar mass ratios of 0.19% (0.14%) at $z \sim 9$ ($z \sim 11$). Dust attenuation decreases with redshift due to dust being increasingly more dispersed within the halo. (iv) the galaxies observed by ALMA at $z \sim 7$ comprise a biased sample that is not representative of the average population that makes up the UV LF.

A very pressing question in contemporary physics is the identity of Dark Matter (DM), and one that has not been answered affirmatively to any degree so far. Primordial Black Holes (PBHs) are one of the most well-motivated DM candidates. Light enough PBHs have been constrained by either the non-detection of their Hawking radiation itself, or by the non-observation of any measurable effects of this radiation on astrophysical and cosmological observables. We constrain the PBH density by their Hawking radiation effect on the intergalactic medium (IGM) temperature evolution. We use the latest deductions of IGM temperature from Lyman-$\alpha$ forest observations. We put constraints on the fraction of PBH DM with masses $5 \times 10^{15}$ g - $10^{17}$ g separately for spinning and non-spinning BHs. We derive constraints by dealing with the heating effects of the astrophysical reionization of the IGM in two ways. In one way, we completely neglect this heating due to astrophysical sources, thus giving us weaker constraints, but completely robust to the reionization history of the universe. In the second way, we utilise some modelling of the ionization and temperature history, and use it to derive more stringent constraints. We find that for non-spinning PBHs of mass $10^{16}$ g, the current measurements can constrain the PBH-density to be $\lesssim$ 0.1% of the total DM. We find that these constraints from the latest Lyman-$\alpha$ forest temperature measurements are competitive, and hence provide a new observable to probe the nature of PBH DM. The systematics affecting Lyman-$\alpha$ forest measurements are different from other constraining observations, and thus this is a complementary probe.

Extreme-mass ratio inspirals (EMRIs) of stellar-mass black holes (BHs) are among the main targets for upcoming low-frequency gravitational wave (GW) detectors such as the Laser Interferometer Space Antenna (LISA). In the classical scenario, EMRIs are formed when BHs scatter off each other and are driven onto highly eccentric orbits that gradually inspiral due to GW emission. If the cluster is in a state of strong mass segregation, the BHs are expected to reside in a steep cusp around the central massive black hole (MBH), which would facilitate more efficient EMRI formation. However, strong mass segregation may also lead to an increased rate of ejections due to close encounters between the BHs. Here, we test the relevance of such ejections for EMRI formation by numerically solving a two-dimensional Fokker-Planck equation. Our formalism includes the effects of two-body relaxation, GW dissipation, and ejections. We find that the EMRI formation rate can be suppressed due to ejections by more than an order of magnitude for strongly segregated BH cusps with density index $\gamma\gtrsim 2.25$ around central MBHs of mass $M_{\bullet} \lesssim 10^6 M_\odot $. The EMRI formation rate levels off up to a maximum value of $\simeq 200~{\rm Gyr}^{-1}$ due to ejections, which is roughly an order of magnitude lower than the usual scenarios ignoring ejections for steep BH cusps around low mass MBHs. Our analysis brings forth the significance of strong scatterings for EMRI formation in galactic nuclei.

Context: New spectroscopic surveys will increase the number of astronomical objects requiring characterization by over tenfold.. Machine learning tools are required to address this data deluge in a fast and accurate fashion. Most machine learning algorithms can not estimate error directly, making them unsuitable for reliable science. Aims: We aim to train a supervised deep-learning algorithm tailored for high-resolution observational stellar spectra. This algorithm accurately infer precise estimates while providing coherent estimates of uncertainties by leveraging information from both the neural network and the spectra. Methods: We train a conditional Invertible Neural Network (cINN) on observational spectroscopic data obtained from the GIRAFFE spectrograph (HR10 and HR21 setups) within the Gaia-ESO survey. A key features of cINN is its ability to produce the Bayesian posterior distribution of parameters for each spectrum. By analyzing this distribution, we inferred parameters and their uncertainties. Several tests have been applied to study how parameters and errors are estimated. Results: We achieved an accuracy of 28K in $T_{\text{eff}}$, 0.06 dex in $\log g$, 0.03 dex in $[\text{Fe/H}]$, and between 0.05 dex and 0.17 dex for the other abundances for high quality spectra. Accuracy remains stable with low signal-to-noise ratio spectra. The uncertainties obtained are well within the same order of magnitude. The network accurately reproduces astrophysical relationships both on the scale of the Milky Way and within smaller star clusters. We created a table containing the new parameters generated by our cINN. Conclusion: This neural network represents a compelling proposition for future astronomical surveys. These coherent derived uncertainties make it possible to reuse these estimates in other works as Bayesian priors and thus present a solid basis for future work.

While there has been an increase in interest in the possibility of quasi-periodic oscillations (QPOs) in blazars, the search has hitherto been restricted to sources with well-sampled light curves. Objects with light curves that include gaps have been, to our knowledge, overlooked. Here, we study two such curves, which have the interesting feature of pertaining to relatively high redshift blazars -- FSRQs, PKS 2155-83 and PKS 2255-282 -- observed by Fermi-LAT. Their redshifts border the 'cosmic noon' era of galaxy formation and merging, and their light curves exhibit a distinctive pattern of repetitive high and low (gap dominant) states for $15.6$ years. To accommodate for the gaps in the curves, data is integrated over extended time intervals of 1 month and 2 months. The resulting curves were also examined using methods suitable for sparsely sampled data. This investigation of PKS 2155-83 and PKS 2255-282 suggests QPOs with periods of $4.69\pm0.79$ yr ($3\sigma$) and $6.82\pm2.25$ yr ($2.8\sigma$), respectively. The flux PDFs of the blazars, along with the correlation between their flux and spectral index, were also analyzed. Given the epochs the objects are observed, the plausibility of a binary black hole scenario as an origin of the apparent periodicity was examined. We estimated the prospective parameters of such a system using a simple geometric model. The total masses were estimated, and found to be consistent, in principle, with independent (dynamical) measurements of the central black hole masses in the two host galaxies.

$f(R)$ theories of modified gravity may be compatible with current observations if the deviations from general relativity are sufficiently well screened in dense environments. In recent work [arXiv:2310.19955] we have shown that approximations commonly used to assess whether galaxies are screened, or unscreened, fail to hold in observationally interesting parts of parameter space. One of the assumptions commonly made in these approximations, and more broadly in the study of $f(R)$ models, is that the mass of the scalar mode can be neglected inside a galaxy. In this work we demonstrate that this approximation may fail spectacularly and discuss the implications of this for tests of the theory.

Context. Computational astronomy has reached the stage where running a gravitational N-body simulation of a stellar system, such as a Milky Way star cluster, is computationally feasible, but a major limiting factor that remains is the ability to set up physically realistic initial conditions. Aims. We aim to obtain realistic initial conditions for N-body simulations by taking advantage of machine learning, with emphasis on reproducing small-scale interstellar distance distributions. Methods. The computational bottleneck for obtaining such distance distributions is the hydrodynamics of star formation, which ultimately determine the features of the stars, including positions, velocities, and masses. To mitigate this issue, we introduce a new method for sampling physically realistic initial conditions from a limited set of simulations using Gaussian processes. Results. We evaluated the resulting sets of initial conditions based on whether they meet tests for physical realism. We find that direct sampling based on the learned distribution of the star features fails to reproduce binary systems. Consequently, we show that physics-informed sampling algorithms solve this issue, as they are capable of generating realisations closer to reality.

We present the discovery of a luminous X-ray AGN in the dwarf galaxy merger RGG 66. The black hole is predicted to have a mass of $M_{\rm BH} \sim 10^{5.4} M_\odot$ and to be radiating close to its Eddington limit ($L_{\rm bol}/L_{\rm Edd} \sim 0.75$). The AGN in RGG 66 is notable both for its presence in a late-stage dwarf-dwarf merger and for its luminosity of $L_{\rm 2-10~keV} = 10^{42.2}$ erg s$^{-1}$, which is among the most powerful AGNs known in nearby dwarf galaxies. The X-ray spectrum has a best-fit photon index of $\Gamma = 2.4$ and an intrinsic absorption of $N_H \sim 10^{21}$ cm$^{-2}$. These results come from a follow-up {\it Chandra X-ray Observatory} study of four irregular/disturbed dwarf galaxies with evidence for hosting AGNs based on optical spectroscopy. The remaining three dwarf galaxies do not have detectable X-ray sources with upper limits of $L_{\rm 2-10~ keV} \lesssim 10^{40}$ erg s$^{-1}$. Taken at face value, our results on RGG 66 suggest that mergers may trigger the most luminous of AGNs in the dwarf galaxy regime, just as they are suspected to do in more massive galaxy mergers.

We obtain spatially resolved kinematics with the Keck Cosmic Web Imager (KCWI) integral-field spectrograph for a sample of 14 massive (11 < log10 M* < 12) lensing early-type galaxies (ETGs) at redshifts z=0.15-0.35 from the Sloan Lens ACS (SLACS) survey. We integrate within the galaxy effective radius and examine the rotational and dispersion velocities, showing that 11/14 are quantitatively classified as slow rotators in comparison with local galaxy surveys. Of key interest is the ability of this data to enable the precision required for cosmological inference with lensing time delays on scales of 1-2% uncertainty. The dataset is unprecedented for galaxy-scale lens galaxies, in terms of signal-to-noise ratio, sampling, and calibration. We test sources of systematic error and identify primary contributions from choice of stellar template library and wavelength range of the spectral fit. Systematics are quantified at the spatial bin level, resulting in systematic error at 3% and positive spatial covariance of 2%. We examine the effects of integration of the kinematic maps within circular apertures of different sizes and compare with SDSS single-aperture velocity dispersions. The most recent velocity dispersion estimates from SDSS spectra are found to be biased by a factor of 5.3% with respect to KCWI data, and to underestimate uncertainties. We examine correlations between scaling relations and show the correlations to agree with previous SLACS analysis with no statistically significant disagreement. A follow-up paper will present Jeans modeling and discuss the context of these observations within broader studies of galaxy evolution and cosmology.

The existence of axions or Axion-Like Particles (ALPs) has been predicted by various Beyond Standard Model (BSM) theories, and the proposed photon-ALP interaction is one of the ways to probe them. Such an interaction will lead to photon-ALP resonant conversion in galaxy clusters, resulting in a polarized spectral distortion in the CMB along the cluster line of sight. The estimation of this signal from galaxy clusters requires an estimation of the electron density and magnetic field in galaxy clusters. We have developed a new Bayesian framework \texttt{SpectrAx} that can use observations from different electromagnetic bands such as radio, CMB, optical, and X-ray to infer the astrophysical properties of a galaxy cluster, such as cluster its redshift, electron density and magnetic field, along with the BSM physics such as ALPs. By using the simulated data for upcoming CMB surveys such as Simons Observatory (SO) and CMB-S4 in combination with Square Kilometer Array (SKA) and extended ROentgen Survey with an Imaging Telescope Array (eROSITA) we demonstrate the capability in accurately inferring the ALPs coupling strength along with the radial profile of electron density and magnetic field from galaxy clusters. The application of this framework to the data from future surveys by combining SKA+SO+eROSITA and SKA+CMB-S4+eROSITA will make it possible for the first time to explore both astrophysics and BSM physics from low-redshift galaxy clusters using a multi-band approach.

Meteorites trace planet formation in the Sun's protoplanetary disk, but they also record the influence of the Sun's birth environment. Whether the Sun formed in a region like Taurus-Auriga with ~10^2 stars, or a region like the Carina Nebula with ~10^6 stars, matters for how large the Sun's disk was, for how long and from how far away it accreted gas from the molecular cloud, and how it acquired radionuclides like 26Al. To provide context for the interpretation of meteoritic data, we review what is known about the Sun's birth environment. Based on an inferred gas disk outer radius ~50-90 AU, radial transport in the disk, and the abundances of noble gases in Jupiter's atmosphere, the Sun's molecular cloud and protoplanetary disk were exposed to an ultraviolet flux G0 ~30-3000 during its birth and first ~10 Myr of evolution. Based on the orbits of Kuiper Belt objects, the Solar System was subsequently exposed to a stellar density ~100 Msol/pc^3 for ~100 Myr, strongly implying formation in a bound cluster. These facts suggest formation in a region like the outskirts of the Orion Nebula, perhaps 2 pc from the center. The protoplanetary disk might have accreted gas for many Myr, but a few x10^5 yr seems more likely. It probably inherited radionuclides from its molecular cloud, enriched by inputs from supernovae and especially Wolf-Rayet star winds, and acquired a typical amount of 26Al.

The detection of the gravitational-wave event GW230529, presumably a neutron star-black hole (NSBH) merger, by the LIGO-Virgo-KAGRA (LVK) Collaboration is an exciting discovery for multimessenger astronomy. The black hole (BH) has a high probability of falling within the ''mass gap'' between the peaks of the neutron star (NS) and the BH mass distributions. Because of the low primary mass, the binary is more likely to produce an electromagnetic counterpart than previously detected NSBH mergers. We investigate the possible kilonova (KN) emission from GW230529, and find that if it was an NSBH, there is a $\sim$ 2-41% probability (depending on the assumed equation of state) that GW230925 produced a KN with magnitude peaking at $\sim 1-2$ day post merger at $g \lesssim 23.5$, $i<23$. Hence, it could have been detected by ground-based telescopes. If it was a binary neutron star (BNS) merger, we find $\sim$ 0-12% probability that it produced a KN. Motivated by these numbers, we simulated a broader population of mgNSBH mergers that may be detected in O4, and we obtained a 9-21% chance of producing a KN, which would be detectable with $g\lesssim 25$ and $ i \lesssim 24$, typically fainter than what is expected from GW230529. Based on these findings, DECam-like instruments may be able to detect up to 80% of future mgNSBH KNe, thus up to $\sim1$ multimessenger mgNSBH per year may be discoverable at the current level of sensitivity (O4).

We present new results from the ASTRID simulation from $z=3$ to $z=0.5$, covering the epoch of cosmic noon. The galaxy stellar mass function, as well as the black hole mass and luminosity functions in ASTRID, exhibit good agreement with recent observational constraints. We study the $M_{\rm BH}$-$M_*$ scaling relation and its connections to AGN luminosity, galaxy color, and star formation rate, demonstrating that AGN feedback plays a crucial role in the quenching of massive galaxies ($M_*>10^{10.5} M_{\odot}$). Although AGN feedback suppresses star formation through quenching, AGN-host galaxies still exhibit statistically higher levels of star formation compared to inactive ones, due to the positive correlation between AGN activity and star formation, both fueled by a shared gas reservoir. The fraction of quiescent galaxies in ASTRID increases with both galaxy mass and redshift evolution, aligning well with observational trends. We find that different quenching mechanisms can leave distinct morphological imprints on quenched galaxies. Massive, compact quiescent galaxies typically experience shorter quenching timescales, have younger central regions, and host overmassive black holes. This is usually due to a compaction-like quenching mechanism that funnels gas into the galaxy center, leading to starbursts and triggering AGN kinetic feedback. In contrast, quiescent galaxies with more diffuse morphologies generally experience `inside-out' quenching, which is characterized by older central regions compared to the outskirts. These galaxies typically experience longer quenching timescales due to quenching processes operating on a larger halo scale, which gradually deplete the galactic star-forming gas. Data of the \astrid simulation down to $z=0.5$ is available at \url{this https URL}.

Future large, space-based observatories with starlight suppression technology, e.g., the Habitable Worlds Observatory (HWO), will directly image and characterize nearby Earth-like exoplanets. Prior limits on planet masses and system architectures from radial velocity (RV) measurements of potential exo-Earth hosts are critical to the success of HWO's science goals. Here, we present a uniform analysis of archival RVs from HIRES/Keck and HARPS/ESO of the most promising targets for the HWO exo-Earth survey. We analyze RVs and stellar activity indicators of 90 stars in the NASA ExEP Mission Star List and SPORES-HWO Catalog, finding 33 Keplerian signals associated with known planets and 12 signals associated with stellar activity. We also identify 5 new RV signals that we classify as either planet candidates or sources requiring confirmation, noting that the RV observations are biased toward cooler and less active stars. Assessing the sensitivity of the HIRES and HARPS data, we calculate RV limits ranging from $K_{\rm RV} = 0.6 \,{\rm m\,s}^{-1}$ (HD 10700) to $371 \,{\rm m\,s}^{-1}$ (HD 17925) in the middle of the conservative habitable zone (HZ), corresponding to projected planet masses of $5.4 \,{\rm M_\oplus}$ and $10.6 \,{\rm M_{Jup}}$ for those stars. The median HZ sensitivity limit of our sample is $M_{\rm p} \sin i \simeq 66 \,{\rm M_\oplus}$. This work demonstrates the need for future extreme precision radial velocity (EPRV) monitoring of high-priority targets for the next generation of DI missions that will search for habitable extrasolar systems. We advocate for the use of these results in developing future EPRV strategies.

This paper introduces a novel method for creating spectral time series, which can be used for generating synthetic light curves for photometric classification but also for applications like K-corrections and bolometric corrections. This approach is particularly valuable in the era of large astronomical surveys, where it can significantly enhance the analysis and understanding of an increasing number of SNe, even in the absence of extensive spectroscopic data. methods: By employing interpolations based on optimal transport theory, starting from a spectroscopic sequence, we derive weighted average spectra with high cadence. The weights incorporate an uncertainty factor, for penalizing interpolations between spectra with significant epoch differences and with poor match between the synthetic and observed photometry. results: Our analysis reveals that even with phase difference of up to 40 days between pairs of spectra, optical transport can generate interpolated spectral time series that closely resemble the original ones. Synthetic photometry extracted from these spectral time series aligns well with observed photometry. The best results are achieved in the V band, with relative residuals less than 10% for 87% and 84% of the data for type Ia and II, respectively. For the B, g, R and r bands the relative residuals are between 65% and 87% within the previously mentioned 10% threshold for both classes. The worse results correspond to the i and I bands where, in the case, of SN~Ia the values drop to 53% and 42%, respectively. conclusions: We introduce a new method to construct spectral time series for individual SN starting from a sparse spectroscopic sequence, demonstrating its capability to produce reliable light curves that can be used for photometric classification.

Recent advances in machine learning algorithms have unlocked new insights in observational astronomy by allowing astronomers to probe new frontiers. In this article, we present a methodology to disentangle the intrinsic X-ray spectrum of galaxy clusters from the instrumental response function. Employing state-of-the-art modeling software and data mining techniques of the Chandra data archive, we construct a set of 100,000 mock Chandra spectra. We train a recurrent inference machine (RIM) to take in the instrumental response and mock observation and output the intrinsic X-ray spectrum. The RIM can recover the mock intrinsic spectrum below the 1-$\sigma$ error threshold; moreover, the RIM reconstruction of the mock observations are indistinguishable from the observations themselves. To further test the algorithm, we deconvolve extracted spectra from the central regions of the galaxy group NGC 1550, known to have a rich X-ray spectrum, and the massive galaxy clusters Abell 1795. Despite the RIM reconstructions consistently remaining below the 1-$\sigma$ noise level, the recovered intrinsic spectra did not align with modeled expectations. This discrepancy is likely attributable to the RIM's method of implicitly encoding prior information within the neural network. This approach holds promise for unlocking new possibilities in accurate spectral reconstructions and advancing our understanding of complex X-ray cosmic phenomena.

Approximately one billion years (Gyr) in the future, as the Sun brightens, Earth's carbonate-silicate cycle is expected to drive CO$_2$ below the minimum level required by vascular land plants, eliminating most macroscopic land life. Here, we couple global-mean models of temperature- and CO$_2$-dependent plant productivity for C$_3$ and C$_4$ plants, silicate weathering, and climate to re-examine the time remaining for terrestrial plants. If weathering is weakly temperature-dependent (as recent data suggest) and/or strongly CO$_2$-dependent, we find that the interplay between climate, productivity, and weathering causes the future luminosity-driven CO$_2$ decrease to slow and temporarily reverse, averting plant CO$_2$ starvation. This dramatically lengthens plant survival from 1 Gyr up to $\sim$1.6-1.86 Gyr, until extreme temperatures halt photosynthesis, suggesting a revised kill mechanism for land plants and potential doubling of the future lifespan of Earth's land macrobiota. An increased future lifespan for the complex biosphere may imply that Earth life had to achieve a smaller number of ``hard steps'' (unlikely evolutionary transitions) to produce intelligent life than previously estimated. These results also suggest that complex photosynthetic land life on Earth and exoplanets may be able to persist until the onset of the moist greenhouse transition.

Fast radio bursts (FRBs) are bright extragalactic transients likely produced by magnetars. We study the propagation of FRBs in magnetar winds, assuming that the wind is strongly magnetized and composed of electron-positron pairs. We focus on the regime where the strength parameter of the radio wave, $a_0$, is larger than unity, and the wave frequency, $\omega_0$, is larger than the Larmor frequency in the background magnetic field, $\omega_{\rm L}$. We show that strong radio waves with $a_0>1$ are able to propagate when $\omega_0 > a_0\omega_{\rm L}$, as the plasma current is a linear function of the wave electric field. The dispersion relation is independent of the wave strength parameter when $\omega_0 > a_0\omega_{\rm L}$. Instead, radio waves could be damped when $\omega_0 < a_0\omega_{\rm L}$, as a significant fraction of the wave energy is used to compress the plasma and amplify the background magnetic field. Our results suggest that FRBs should be produced at large distances from the magnetar (i.e., $R>10^{12}{\rm\; cm}$, where the condition $\omega_0 > a_0\omega_{\rm L}$ is satisfied). Alternatively, the structure of the magnetar wind should be strongly modified during a flare to allow the escape of FRBs produced at radii $R<10^{12}{\rm\; cm}$.

We develop a theoretical framework and use two-dimensional hydrodynamical simulations to study the repulsive effect between two close orbiters embedded in an accretion disk. We consider orbiters on fixed Keplerian orbits with masses low enough to open shallow gaps. The simulations indicate that the repulsion is larger for more massive orbiters and decreases with the orbital separation and the disk's viscosity. We use two different assumptions to derive theoretical scaling relations for the repulsion. A first scenario assumes that each orbiter absorbs the angular momentum deposited in its horseshoe region by the companion's wake. A second scenario assumes that the corotation torques of the orbiters are modified because the companion changes the underlying radial gradient of the disk surface density. We find a substantial difference between the predictions of these two scenarios. The first one fails to reproduce the scaling of the repulsion with the disk viscosity and generally overestimates the strength of the repulsion. The second scenario, however, gives results that are broadly consistent with those obtained in the simulations.

New observational facilities are beginning to enable insights into the three-dimensional (3D) nature of exoplanets. Transmission spectroscopy is the most widely used method for characterizing transiting temperate exoplanet's atmospheres, but because it only provides a glimpse of the planet's limb and nightside for a typical orbit, its ability to probe 3D characteristics is still an active area of research. Here, we use the ROCKE-3D general circulation model to test the impact of rotation rate, a ``low-order'' 3D characteristic previously shown to drive differences in planetary phase curves, on the transmission spectrum of a representative super-Earth across temperate-to-warm instellations (S$_p$=0.8, 1, 1.25, 1.66, 2, 2.5, 3, 4, 4.56 S$_\oplus$). We find that different rotation regimes do display differences in their transmission spectra, primarily driven by clouds and humidity, and that the differences shrink or disappear in hotter regimes where water clouds are unable to condense (though our simulations do not consider haze formation). The small size of the feature differences and potential for degeneracy with other properties, like differing water content or atmospheric structure, mean that we do not specifically claim to have identified a single transmission diagnostic for rotation rate, but our results can be used for holistic spectrum interpretation and sample creation, and suggest the need for more modelling in this area.

This study explores the extension of teleparallel gravity within the framework of general relativity, introducing an algebraic function $f(T)$ dependent on the torsion scalar $T$. Motivated by the teleparallel formulation, we investigate cosmological implications, employing the simplest parametrization of the dark energy equation of state. Our chosen $f(T)$ function, $f(T)=\alpha(-T)^n$, undergoes stringent constraints using recent observational data ($H(z)$, SNeIa, BAO, and CMB). The model aligns well with cosmic dynamics, exhibiting quintessence behavior. The evolution of the deceleration parameter, the behavior of dark energy components, and the $Om(z)$ diagnostic further reveal intriguing cosmological phenomena, emphasizing the model's compatibility with quintessence scenarios.

The way pulsars spin down is not understood in detail, but a number of possible physical mechanisms produce a spin-down rate that scales as a power of the rotation rate ($\dot\nu\propto-\nu^n$), with the power-law index $n$ called the braking index. PSR B0540-69 is a pulsar that in 2011, after 16 years of spinning down with a constant braking index of 2.1, experienced a giant spin-down change and a reduction of its braking index to nearly zero. Here, we show that following this episode the braking index monotonically increased during a period of at least four years and stabilised at ~1.1. We also present an alternative interpretation of a more modest rotational irregularity that occurred in 2023, which was modelled as an anomalous negative step of the rotation rate. Our analysis shows that the 2023 observations can be equally well described as a transient swing of the spin-down rate (lasting ~65 days), and the Bayesian evidence indicates that this model is strongly preferred.

Irradiation from the central star controls the temperature structure in protoplanetary discs. Yet simulations of gravitational instability typically use models of stellar irradiation with varying complexity, or ignore it altogether, assuming heat generated by spiral shocks is balanced by cooling, leading to a self-regulated state. In this paper, we perform simulations of irradiated, gravitationally unstable protoplanetary discs using 3D hydrodynamics coupled with live Monte-Carlo radiative transfer. We find that the resulting temperature profile is approximately constant in time, since the thermal effects of the star dominate. Hence, the disc cannot regulate gravitational instabilities by adjusting the temperatures in the disc. In a 0.1 Solar mass disc, the disc instead adjusts by angular momentum transport induced by the spiral arms, leading to steadily decreasing surface density, and hence quenching of the instability. Thus, strong spiral arms caused by self-gravity would not persist for longer than ten thousand years in the absence of fresh infall, although weak spiral structures remain present over longer timescales. Using synthetic images at 1.3mm, we find that spirals formed in irradiated discs are challenging to detect. In higher mass discs, we find that fragmentation is likely because the dominant stellar irradiation overwhelms the stabilising influence of PdV work and shock heating in the spiral arms.

Recent determinations of the total rate of the 12C+12C nuclear reaction show non-negligible differences with the reference reaction rate commonly used in previous stellar simulations. In addition, the current uncertainties in determining each exit channel constitute one of the main uncertainties in shaping the inner structure of super asymptotic giant branch stars that could have a measurable impact on the properties of pulsating ultra-massive white dwarfs (WDs). We explore how new determinations of the nuclear reaction rate and its branching ratios affect the evolution of WD progenitors. We show that the current uncertainties in the branching ratios constitute the main uncertainty factor in determining the inner composition of ultra-massive WDs and their progenitors. We found that the use of extreme branching ratios leads to differences in the central abundances of 20Ne of at most 17%, which are translated into differences of at most 1.3 and 0.8% in the cooling times and size of the crystallized core. However, the impact on the pulsation properties is small, less than 1 s for the asymptotic period spacing. We found that the carbon burns partially in the interior of ultra-massive WD progenitors within a particular range of masses, leaving a hybrid CONe-core composition in their cores. The evolution of these new kinds of predicted objects differs substantially from the evolution of objects with pure CO cores. Differences in the size of the crystallized core and cooling times of up to 15 and 6%, respectively leading to distinct patterns in the period spacing distribution.

The eruptive YSO FU Ori is the eponym of its variable class. FU Ori stars are known to undergo outbursts with amplitudes of $>4$ magnitudes in the $V$ band and durations of several decades. Interaction with a binary companion is one proposed outburst trigger, so understanding both components of the FU Ori system is crucial. A recent HST/STIS observation of the FU Ori system clearly resolves its North and South components. We report here on the spectrum of FU Ori South. We detect NUV continuum emission but no FUV continuum, although several bright emission lines consistent with those seen in T Tauri stars are present. The presence of the C II] 2325 multiplet and many H$_2$ lines indicate active accretion. We estimate the extinction to the source and find that the UV spectrum favors $A_V < 4$, contrary to past estimates based on the NIR spectrum.

Observations show that globular clusters might be among the best places to find millisecond pulsars. However, the globular cluster Terzan 6 seems to be an exception without any pulsar discovered, although its high stellar encounter rate suggests that it harbors dozens of them. We report the discovery of the first radio pulsar, PSR J1751-3116A, likely associated with Terzan 6 in a search of C-band (4-8 GHz) data from the Green Bank Telescope with a spin period of 5.33 ms and dispersion measure, DM$\simeq$383 ${\rm pc~cm^{-3}}$. The mean flux density of this pulsar is approximately 3 ${\rm \mu Jy}$. The DM agrees well with predictions from the Galactic free electron density model, assuming a distance of 6.7 kpc for Terzan 6. PSR J1751-3116A is likely an isolated millisecond pulsar, potentially formed through dynamical interactions, considering the core-collapsed classification and the exceptionally high stellar encounter rate of Terzan 6. This is the highest radio frequency observation that has led to the discovery of a pulsar in a globular cluster to date. While L-band (1-2 GHz) observations of this cluster are unlikely to yield significant returns due to propagation effects, we predict that further pulsar discoveries in Terzan 6 will be made by existing radio telescopes at higher frequencies.

On April 25th, 2019, the LIGO-Virgo Collaboration discovered a Gravitational-wave (GW) signal from a binary neutron star (BNS) merger, i.e., GW190425. Due to the inferred large total mass, the origin of GW190425 remains unclear. We perform detailed stellar structure and binary evolution calculations that take into account mass-loss, internal differential rotation, and tidal interactions between a He-rich star and a NS companion. We explore the parameter space of the initial binary properties, including initial NS and He-rich masses and initial orbital period. We find that the immediate post-common-envelope progenitor system, consisting of a primary $\sim2.0\,M_\odot$ ($\sim1.7\,M_\odot$) NS and a secondary He-rich star with an initial mass of $\sim3.0-5.5\,M_\odot$ ($\sim5.5-6.0\,M_\odot$) in a close binary with an initial period of $\sim0.08-0.5\,{\rm{days}}$ ($\sim 0.08-0.4\,{\rm{days}}$), that experiences stable Case BB/BC mass transfer (MT) during binary evolution, can reproduce the formation of GW190425-like BNS events. Our studies reveal that the secondary He-rich star of the GW190425's progenitor before its core collapse can be efficiently spun up through tidal interaction, finally remaining as a NS with rotational energy even reaching $\sim10^{52}\,{\rm{erg}}$, which is always much higher than the neutrino-driven energy of the supernova (SN) explosion. If the newborn secondary NS is a magnetar, we expect that GW190425 can be the remnant of a magnetar-driven SN, e.g., a magnetar-driven ultra-stripped SN, a superluminous SN, or a broad-line Type Ic SN. Our results show that GW190425 could be formed through the isolated binary evolution, which involves a stable Case BB/BC MT just after the common envelope phase. On top of that, we show the He-rich star can be tidally spun up, potentially forming a spinning magnetized NS (magnetar) during the second SN explosion.

In this article a new method is proposed for estimating the mass composition of cosmic rays in individual events with energies above $1.25 \times 10^{19}$ eV. It is based on a joint analysis of experimental data and simulation results obtained using the QGSJet-II.04 model for muons with threshold energy $E_{\mu} = 1.0 \times \cos\theta$ GeV in air showers with zenith angles up to 60 degrees. The data from ground-based and underground scintillation detectors of the Yakutsk EAS array were used. Separate groups of nuclei and other primary particles were found.

We address the critical need for accurate Rosseland mean gas opacities in high-pressure environments, spanning temperatures from 100 K to 32000 K. Current opacity tables from Wichita State University and AESOPUS 2.0 are limited to $\log(R) \le 1$, where $R=\rho\, T_6^{-3}$ in units of $\mathrm{g}\,\mathrm{cm}^{-3}(10^6\mathrm{K})^{-3}$. This is insufficient for modeling very low-mass stars, brown dwarfs, and planets with atmospheres exhibiting higher densities and pressures ($\log(R) > 1$). Leveraging extensive databases such as ExoMol, ExoMolOP, MoLLIST, and HITEMP, we focus on expanding the AESOPUS opacity calculations to cover a broad range of pressure and density conditions ($-8 \leq \log(R) \leq +6$). We incorporate the thermal Doppler mechanism and micro-turbulence velocity. Pressure broadening effects on molecular transitions, leading to Lorentzian or Voigt profiles, are explored in the context of atmospheric profiles for exoplanets, brown dwarfs, and low-mass stars. We also delve into the impact of electron degeneracy and non-ideal effects such as ionization potential depression under high-density conditions, emphasizing its notable influence on Rosseland mean opacities at temperatures exceeding $10,000$ K. As a result, this study expands AESOPUS public web interface for customized gas chemical mixtures, promoting flexibility in opacity calculations based on specific research needs. Additionally, pre-computed opacity tables, inclusive of condensates, are provided. We present a preliminary application to evolutionary models for very low-mass stars.

In this article, we extend the standard $\Lambda$CDM model by incorporating causal viscous dark matter (vDM) modelled using Maxwell-Cattaneo theory and obtain a novel Hubble parameter solution for the model, which can even be extended to incorporate both radiation and baryonic matter components. A detailed and comprehensive analysis of the model is carried out by deriving theoretical constraints on model parameters, comparing the model with the latest observational data sets and validating the behaviour of the model under different thermodynamic laws. Such an in-depth analysis of the model yielded several intriguing results, like the presence of sign-switching bulk viscous pressure which aids both, the decelerated expansion in the early universe and the accelerated expansion in the late universe, possibility of having negative specific entropy rate in the early Universe while still satisfying the covariant second law of thermodynamics, and a correlation between relaxation time parameter and sign-switching redshift. Finally, we propose a unified dark matter (UDM) interpretation for the dark sector in this model and hence showed that, under the said UDM interpretation, this model can satisfy the much-required near equilibrium condition associated with the background dissipative theory both in the early and late phase of the Universe. Based on the NEC trend showcased by the model, we conjuncture that the late-accelerating behaviour of Universe can be viewed as a means for the UDM fluid to attain a global equilibrium state.

Forward models of the galaxy density field enable simulation based inference as well as field level inference of galaxy clustering. However, these analysis techniques require forward models that are both computationally fast and robust to modeling uncertainties in the relation between galaxies and matter. Both requirements can be addressed with the Effective Field Theory of Large Scale Structure. Here, we focus on the physical and numerical convergence of the LEFTfield model. Based on the perturbative nature of the forward model, we derive an analytic understanding of the leading numerical errors, and we compare our estimates to high-resolution and N-body references. This allows us to derive a set of best-practice recommendations for the numerical accuracy parameters, which are completely specified by the desired order of the perturbative solution and the cut-off scale. We verify these recommendations by an extended set of parameter recovery tests from fully nonlinear mock data and find very consistent results. A single evaluation of the forward model takes seconds, making cosmological analyses of galaxy clustering data based on forward models computationally feasible.

Using a combination of HST, JWST, and ALMA data, we perform spatially resolved spectral energy distributions (SED) fitting of fourteen 4<z<6 UV-selected main-sequence galaxies targeted by the [CII] Resolved ISM in Star-forming Galaxies with ALMA (CRISTAL) Large Program. We consistently model the emission from stars and dust in ~0.5-1kpc spatial bins to obtain maps of their physical properties. We find no offsets between the stellar masses (M*) and star formation rates (SFRs) derived from their global emission and those from adding up the values in our spatial bins, suggesting there is no bias of outshining by young stars on the derived global properties. We show that ALMA observations are important to derive robust parameter maps because they reduce the uncertainties in Ldust (hence Av and SFR). Using these maps we explore the resolved star-forming main sequence for z~5 galaxies, finding that this relation persists in typical star-forming galaxies in the early Universe. We find less obscured star formation where the M* (and SFR) surface densities are highest, typically in the central regions, contrary to the global relation between these parameters. We speculate this could be caused by feedback driving gas and dust out of these regions. However, more observations of infrared luminosities with ALMA are needed to verify this. Finally, we test empirical SFR prescriptions based on the UV+IR and [CII] line luminosity, finding they work well at the scales probed (~kpc). Our work demonstrates the usefulness of joint HST, JWST, and ALMA resolved SED modeling analyses at high redshift.

We present a spatially resolved analysis of four star-forming galaxies at $z = 4.44-5.64$ using data from the JWST PRIMER and ALMA-CRISTAL surveys to probe the stellar and inter-stellar medium properties on the sub-kpc scale. In the $1-5\,\mu{\rm m}$ JWST NIRCam imaging we find that the galaxies are composed of multiple clumps (between $2$ and $\sim 8$) separated by $\simeq 5\,{\rm kpc}$, with comparable morphologies and sizes in the rest-frame UV and optical. Using BAGPIPES to perform pixel-by-pixel SED fitting to the JWST data we show that the SFR ($\simeq 25\,{\rm M}_{\odot}/{\rm yr}$) and stellar mass (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot}) \simeq 9.5$) derived from the resolved analysis are in close ($ \lesssim 0.3\,{\rm dex}$) agreement with those obtained by fitting the integrated photometry. In contrast to studies of lower-mass sources, we thus find a reduced impact of outshining of the older (more massive) stellar populations in these normal $z \simeq 5$ galaxies. Our JWST analysis recovers bluer rest-frame UV slopes ($\beta \simeq -2.1$) and younger ages ($\simeq 100\,{\rm Myr}$) than archival values. We find that the dust continuum from ALMA-CRISTAL seen in two of these galaxies correlates, as expected, with regions of redder rest-frame UV slopes and the SED-derived $A_{\rm V}$, as well as the peak in the stellar mass map. We compute the resolved IRX-$\beta$ relation, showing that the IRX is consistent with the local starburst attenuation curve and further demonstrating the presence of an inhomogeneous dust distribution within the galaxies. A comparison of the CRISTAL sources to those from the FirstLight zoom-in simulation of galaxies with the same $M_{\star}$ and SFR reveals similar age and colour gradients, suggesting that major mergers may be important in the formation of clumpy galaxies at this epoch.

We present here results of numerical simulations of the formation and early evolution of rocky planets through pebble accretion, with an with an emphasis on hydrogen envelope longevity and the composition of the outgassed atmosphere. We model planets with a range in mass from 0.1 to 5 Earth masses that orbit between 0.7 and 1.7 AU. The composition of the outgassed atmosphere is calculated with the partial pressure of free oxygen fit to geophysical models of magma ocean self-oxidation. XUV radiation powered photoevaporation is considered as the main driver of atmospheric escape. We model planets that remain below the pebble isolation mass and hence accrete tenuous envelopes only. We consider slow, medium or fast initial stellar rotation for the temporal evolution of the XUV flux. The loss of the envelope is a key event that allows the magma ocean to crystallise and outgas its bulk volatiles. The atmospheric composition of the majority of our simulated planets is dominated by CO$_2$. Our planets accrete a total of 11.6 Earth oceans of water, the majority of which enters the core. The hydrospheres of planets lighter than the Earth reach several times the mass of the Earth's modern oceans, while the hydrospheres of planets ranging from 1 to 3.5 Earth masses are comparable to those of our planet. However, planets of 4-5 Earth masses have smaller hydrospheres due to trapping of volatiles in their massive mantles. Overall, our simulations demonstrate that hydrogen envelopes are easily lost from rocky planets and that this envelope loss triggers the most primordial partitioning of volatiles between the solid mantle and the atmosphere.

Supernova remnants (SNRs) have long been suspected to be the primary sources of Galactic cosmic rays. Over the past decades, great strides have been made in the modelling of particle acceleration, magnetic field amplification, and escape from SNRs. Yet, while many SNRs have been observed in non-thermal emission in radio, X-rays, and gamma-rays, there is no evidence for any individual object contributing to the locally observed flux. Here, we propose a particular spectral signature from individual remnants that is due to the energy-dependent escape from SNRs. For young and nearby sources, we predict fluxes enhanced by tens of percent in narrow rigidity intervals; given the percent-level flux uncertainties of contemporary cosmic-ray data, such features should be readily detectable. We model the spatial and temporal distribution of sources and the resulting distribution of fluxes with a Monte Carlo approach. The decision tree that we have trained on simulated data is able to discriminate with very high significance between the null hypothesis of a smooth distribution of sources and the scenario with a stochastic distribution of individual sources. We suggest that this cosmic-ray energy-dependent injection time (CREDIT) scenario be considered in experimental searches to identify individual SNRs as cosmic-ray sources.

We present observations of mid-J J=4-3 or J=5-4 carbon monoxide (CO) emission lines and continuum emission from a sample of ten of the most luminous log(L/L_solar)~14 Hot Dust-Obscured Galaxies (Hot DOGs) discovered by the Wide-field Infrared Survey Explorer (WISE) with redshifts up to 4.6. We uncover broad spectral lines (FWHM~400 km/s) in these objects, suggesting a turbulent molecular interstellar medium (ISM) may be ubiquitous in Hot DOGs. A halo of molecular gas, extending out to a radius of 5 kpc is observed in W2305-0039, likely supplied by 940 km/s molecular outflows. W0831+0140 is plausibly the host of a merger between at least two galaxies, consistent with observations made using ionized gas. These CO(4-3) observations contrast with previous CO(1-0) studies of the same sources: the CO(4-3) to CO(1-0) luminosity ratios exceed 300 in each source, suggesting that the lowest excited states of CO are underluminous. These findings show that the molecular gas in Hot DOGs is consistently turbulent, plausibly a consequence of AGN feedback, triggered by galactic mergers.

Abridged: We present the first observation of the HCO+(1-0) and HCN(1-0) emission in the northern filaments of Centaurus A with ALMA. HCO+(1-0) is detected in 9 clumps of the Horseshoe complex, with similar velocities as the CO(1-0) emission. Conversely, the HCN(1-0) is not detected and we derive upper limits on the flux. At a resolution of ~40 pc, the line ratio of the velocity-integrated intensities I_HCO+/I_CO varies between 0.03 and 0.08, while I_HCO+/I_HCN is higher than unity with an average lower limit of 1.51. These ratios are significantly higher than what is observed in nearby star-forming galaxies. Moreover, the ratio I_HCO+/I_CO decreases with increasing CO integrated intensity, contrary to what is observed in the star-forming galaxies. This indicates that the HCO+ emission is enhanced and may not arise from dense gas within the Horseshoe complex. This hypothesis is strengthened by the average line ratio I_HCN/I_CO<0.03 which suggests that the gas density is rather low. Using non-LTE, large velocity gradient modelling with RADEX, we explored two possible phases of the gas, that we call "diffuse" and "dense", and are characterised by a significant difference in the HCO+ relative abundance to CO, respectively N_HCO+/N_CO=10^-3 and 3x10^-5. The average CO(1-0) and HCO+(1-0) integrated intensities and the upper limit on HCN(1-0) are compatible with both "diffuse" and "dense" gas. The spectral setup of the present observations also covers the SiO(2-1). While undetected, the upper limit on SiO(2-1) is not compatible with the RADEX predictions for the "dense" gas. We conclude that the 9 molecular clouds detected in HCO+(1-0) are likely dominated by diffuse molecular gas. While the exact origin of the HCO+(1-0) emission remains to be investigated, it is likely related to the energy injection within the molecular gas that prevents gravitational collapse and star formation.

Context. Stellar activity consists of different phenomena, mainly spots and faculae, and it is one of the main sources of noise in exoplanetary observations because it affects both spectroscopic and photometric observations. If we want to study young active planetary systems we need to model the activity of the host stars in order to remove astrophysical noise from our observational data. Aims. We modelled the contribution of stellar spots in photometric observations. Through the use of multiband photometry, we aim to extract the geometric properties of the spots and constrain their temperature. Methods. We analyzed multiband photometric observations acquired with the 80 cm Marcon telescope of the Osservatorio Polifunzionale del Chianti of V1298 Tau, assuming that the photometric modulation observed in different bands should be due to cold spots. Results. We constrained the effective temperature of the active regions present on the surface of V1298 Tau, which is composed by the contemporary presence of spots and faculae. We tested our hypothesis on solar data, verifying that we measure the size of the dominant active region and its averaged effective temperature.

The Roman Space Telescope will be a critical mission to demonstrate high-contrast imaging technologies allowing for the characterisation of exoplanets in reflected light. It will demonstrate $10^{-7}$ contrast limits or better at 3--9 $\lambda / D$ separations with active wavefront control for the first time in space. The detection limits for the Coronagraph Instrument are expected to be set by wavefront variations between the science target and the reference star observations. We are investigating methods to use the deformablel mirrors to methodically probe the impact of such variations on the coronagraphic PSF, generating a PSF library during observations of the reference star to optimise the starlight subtraction at post-processing. We are collaborating with STScI to test and validate these methods in lab using the HiCAT tested, a high-contrast imaging lab platform dedicated to system-level developments for future space missions. In this paper, we will present the first applications of these methods on HiCAT.

Rocky planets may acquire a primordial atmosphere by outgassing of volatiles from their magma ocean. The distribution of O between H$_2$O, CO and CO$_2$ in chemical equilibrium subsequently changes significantly with decreasing temperature. We explore here two chemical models: one where CH$_4$ and NH$_3$ are assumed to be irrevocably destroyed by photolysis, and one where these molecules persist. In the first case, we show that CO cannot co-exist with H$_2$O, since CO oxidizes at low temperatures to form CO$_2$ and H$_2$. In both cases, H escapes from the thermosphere within a few ten million years by absorption of stellar XUV radiation. This escape drives an atmospheric self-oxidation process whereby rocky planet atmospheres become dominated by CO$_2$ and H$_2$O, regardless of their initial oxidation state at outgassing. HCN is considered a potential precursor of prebiotic compounds and RNA. Our oxidizing atmospheres are inefficient at producing HCN by lightning. Instead, we demonstrate that lightning-produced NO, which dissolves as nitrate in the oceans, and interplanetary dust particles may be the main sources of fixed nitrogen to emerging biospheres. Our results highlight the need for origin-of-life scenarios where the first metabolism fixes its C from CO$_2$, rather than from HCN and CO.

The exoplanet sub-Neptune population currently poses a conundrum. Are small-size planets volatile-rich cores without atmosphere, or are they rocky cores surrounded by H-He envelope? To test the different hypotheses from an observational point of view, a large sample of small-size planets with precise mass and radius measurements is the first necessary step. On top of that, much more information will likely be needed, including atmospheric characterisation and a demographic perspective on their bulk properties. We present the concept and strategy of THIRSTEE, a project which aims at shedding light on the composition of the sub-Neptune population across stellar types by increasing their number and improving the accuracy of bulk density measurements, as well as investigating their atmospheres and performing statistical, demographic analysis. We report the first results of the program, characterising a 2-planet system around the M dwarf TOI-406. We analyse TESS and ground-based photometry, together with ESPRESSO and NIRPS/HARPS RVs to derive the orbital parameters and investigate the internal composition of the 2 planets orbiting TOI-406, which have radii and masses of $R_b = 1.32 \pm 0.12 R_{\oplus}$, $M_b = 2.08_{-0.22}^{+0.23} M_{\oplus}$ and $R_c = 2.08_{-0.15}^{+0.16} R_{\oplus}$, $M_c = 6.57_{-0.90}^{+1.00} M_{\oplus}$, and periods of $3.3$ and $13.2$ days, respectively. Planet b is consistent with an Earth-like composition, while planet c is compatible with multiple internal composition models, including volatile-rich planets without H/He atmospheres. The 2 planets are located in 2 distinct regions in the mass-density diagram, supporting the existence of a density gap among small exoplanets around M dwarfs. With an equilibrium temperature of only 368 K, TOI-406 c stands up as a particularly interesting target for atmospheric characterisation with JWST in the low-temperature regime.

The standard $\rm \Lambda$CDM cosmological model predicts that a large amount of diffuse neutral hydrogen distributes in cosmic filaments, which could be mapped through Lyman-alpha (Ly$\alpha$) emission observations. We use the hydrodynamical simulation Illustris-TNG50 to investigate the evolution of surface brightness and detectability of neutral hydrogen in cosmic filaments across redshifts $z=2-5$. While the HI column density of cosmic filaments decreases with redshift, due to the rising temperature with cosmic time in filaments, the surface brightness of Ly$\alpha$ emission in filaments is brighter at lower redshifts, suggesting that the detection of cosmic filaments is more feasible at lower redshifts. However, most of the Ly$\alpha$ emission from cosmic filaments is around $10^{-21}$ $\rm erg\ s^{-1}cm^{-2}arsec^{-2}$, making it extremely challenging to detect with current observational instruments. We further generate mock images using the Multi-Unit Spectroscopic Explorer (MUSE) spectrograph installed on both the Very Large Telescope (VLT) and the upcoming Extremely Large Telescope (ELT). Our finding indicates that while the VLT can only detect filamentary structures made of dense gas in galactic centers, the ELT is expected to reveal much finer filamentary structures from diffuse neutral hydrogen outside of galaxies. Compared to the VLT, both the number density and the longest length of filaments are greatly boosted with the ELT. Hence the forthcoming ELT is highly promising to provide a clearer view of cosmic filaments in Ly$\alpha$ emission.

The recent discovery of gravitational waves (GWs) has opened a new avenue for investigating the equation of state (EOS) of dense matter in compact stars, which is an outstanding problem in astronomy and nuclear physics. In the future, next-generation (XG) GW detectors will be constructed, deemed to provide a large number of high-precision observations. We investigate the potential of constraining the EOS of quark stars (QSs) with high-precision measurements of mass $m$ and tidal deformability $\Lambda$ from the XG GW observatories. We adopt the widely-used bag model for QSs, consisting of four microscopic parameters: the effective bag constant $B_{\rm eff}$, the perturbative quantum chromodynamics correction parameter $a_4$, the strange quark mass $m_s$, and the pairing energy gap $\Delta$. With the help of hierarchical Bayesian inference, for the first time we are able to infer the EOS of QSs combining multiple GW observations. Using the top 25 loudest GW events in our simulation, we find that, the constraints on $B_{\rm eff}$ and $\Delta$ are tightened by several times, while $a_4$ and $m_s$ are still poorly constrained. We also study a simplified 2-dimensional (2-d) EOS model which was recently proposed in literature. The 2-d model is found to exhibit significant parameter-estimation biases as more GW events are analyzed, while the predicted $m$-$\Lambda$ relation remains consistent with the full model.

The Be/X-ray binary pulsar RX J0440.9+4431 went through a giant outburst in December 2022 with a peak flux of $\sim$2.3 Crab in 15--50 keV. We studied the broad-band timing and spectral properties of RX J0440.9+4431 using four $AstroSat$ observations, where the source transited between subcritical and supercritical accretion regimes. Pulsations were detected significantly above 100 keV. The pulse profiles were found to be highly luminosity- and energy-dependent. A significant evolution in the pulse profile shape near the peak of the outburst indicates a possible change in the accretion mode and beaming patterns of RX J0440.9+4431. The rms pulsed fraction was luminosity- and energy-dependent, with a concave-like feature around 20--30 keV. The depth of this feature varied with luminosity, indicating changes in the accretion column height and proportion of reflected photons. The broad-band continuum spectra were best fitted with a two-component Comptonization model with a blackbody component or a two-blackbody component model with a thermal Comptonization component. A quasi-periodic oscillation at 60 mHz was detected at a luminosity of $2.6 \times 10^{37}$ erg s$^{-1}$, which evolved into 42 mHz at $1.5 \times 10^{37}$ erg s$^{-1}$. The QPO rms were found to be energy dependent with an overall increasing trend with energy. For the first time, we found the QPO frequency varying with photon energy in an X-ray pulsar, which poses a challenge in explaining the QPO with current models such as the Keplarian and beat frequency model. Hence, more physically motivated models are required to understand the physical mechanism behind the mHz QPOs.

Impacts by rocky and icy bodies are thought to have played a key role in shaping the composition of solar system objects, including the Earth's habitability. Hence, it is likely that they play a similar role in exoplanetary systems. We investigate how an icy cometary impact affects the atmospheric chemistry, climate, and composition of an Earth-like, tidally-locked, terrestrial exoplanet, a prime target in the search for a habitable exoplanet beyond our solar system. We couple a cometary impact model which includes thermal ablation and pressure driven breakup with the 3D Earth System Model WACCM6/CESM2, and use this model to investigate the effects of the water and thermal energy delivery associated with an $R=2.5$ km pure water ice cometary impact on an Earth-like atmosphere. We find that water is the primary driver of longer timescale changes to the atmospheric chemistry and composition by acting as a source of opacity, cloud ice, and atmospheric hydrogen/oxygen. The water opacity drives heating at $\sim5\times10^{-4}$ bar, and cooling below, due to a decreased flux reaching the surface. The increase in atmospheric hydrogen and oxygen also drives an increase in the abundance of hydrogen/oxygen rich molecules, with the exception of ozone, whose column density decreases by $\sim10\%$. These atmospheric changes are potentially observable for $\sim$ 1-2 years post-impact, particularly those associated with cloud ice scattering. They also persist, albeit at a much reduced level, to our quasi-steady-state, suggesting that sustained bombardment or multiple large impacts have the potential to shape the composition and habitability of terrestrial exoplanets.

We explore the potential for improving constraints on gravity by leveraging correlations in the dispersion measure derived from Fast Radio Bursts (FRBs) in combination with cosmic shear. Specifically, we focus on Horndeski gravity, inferring the kinetic braiding and Planck mass run rate from a stage-4 cosmic shear mock survey alongside a survey comprising $10^4$ FRBs. For the inference pipeline, we utilise hi_class to predict the linear matter power spectrum in modified gravity scenarios, while non-linear corrections are modelled with HMcode, including feedback mechanisms. Our findings indicate that FRBs can disentangle degeneracies between baryonic feedback and cosmological parameters, as well as the mass of massive neutrinos. Since these parameters are also degenerate with modified gravity parameters, the inclusion of FRBs can enhance constraints on Horndeski parameters by up to $40$ percent, despite being a less significant measurement. Additionally, we apply our model to current FRB data and use the uncertainty in the $\mathrm{DM}-z$ relation to impose limits on gravity. However, due to the limited sample size of current data, constraints are predominantly influenced by theoretical priors. Despite this, our study demonstrates that FRBs will significantly augment the limited set of cosmological probes available, playing a critical role in providing alternative tests of feedback, cosmology, and gravity. All codes used in this work are made publically available.

Vision foundation models, which have demonstrated significant potential in many multimedia applications, are often underutilized in the natural sciences. This is primarily due to mismatches between the nature of domain-specific scientific data and the typical training data used for foundation models, leading to distribution shifts. Scientific data often differ substantially in structure and characteristics; researchers frequently face the challenge of optimizing model performance with limited labeled data of only a few hundred or thousand images. To adapt foundation models effectively requires customized approaches in preprocessing, data augmentation, and training techniques. Additionally, each vision foundation model exhibits unique strengths and limitations, influenced by differences in architecture, training procedures, and the datasets used for training. In this work, we evaluate the application of various vision foundation models to astrophysics data, specifically images from optical and radio astronomy. Our results show that using features extracted by specific foundation models improves the classification accuracy of optical galaxy images compared to conventional supervised training. Similarly, these models achieve equivalent or better performance in object detection tasks with radio images. However, their performance in classifying radio galaxy images is generally poor and often inferior to traditional supervised training results. These findings suggest that selecting suitable vision foundation models for astrophysics applications requires careful consideration of the model characteristics and alignment with the specific requirements of the downstream tasks.

The removal of angular momentum from protostellar systems drives accretion onto the central star and may drive the dispersal of the protoplanetary disk. Winds and jets can contribute to removing angular momentum from the disk, though the dominant process remain unclear. To date, observational studies of resolved disk winds have mostly targeted highly inclined disks. We report the detection of extended H2 and [Ne II] emission toward the young stellar object SY Cha with the JWST Mid-InfraRed Instrument Medium Resolution Spectrometer (MIRI-MRS). This is one of the first polychromatic detections of extended H2 toward a moderately inclined, i=51.1 degrees, Class II source. We measure the semi-opening angle of the H2 emission as well as build a rotation diagram to determine the H2 excitation temperature and abundance. We find a wide semi-opening angle, high temperature, and low column density for the H2 emission, all of which are characteristic of a disk wind. These observations demonstrate MIRI-MRS's utility in expanding studies of resolved disk winds beyond edge-on sources.

Magnetars are neutron stars with superstrong magnetic fields which can exceed 1e15 G. Some magnetars (the so-called soft gamma-repeaters) demonstrate occasionally very powerful processes of energy release, which result in exceptionally strong flares of electromagnetic radiation. It is believed that these flares are associated with the presence of superstrong magnetic fields. Despite many hypotheses, the mechanism of these flares remains a mystery. In afterglows of the flares, one has often observed quasi-periodic oscillations (QPOs) of magnetar emission. They are interpreted as stellar vibrations, excited by the flares, which are useful for exploring the nature of magnetar activity. The incompleteness of theories employed to interpret magnetar QPOs is discussed.

It is generally accepted that all massive galaxies host supermassive black holes (BHs) in their center and that mergers of two galaxies lead to the formation of BH binaries. The most interesting among them comprise the mergers in their final state, that is to say with parsec (3.2 light years) or sub-parsec orbital separations. It is possible to detect these systems with binary self-lensing. Here we report the potential detection of a central supermassive BH binary in the active galaxy (AGN) NGC1566 based on a microlensing outburst. The light curve of the outburst - based on observations with the All Sky Automated Survey for SuperNovae - lasted from the beginning of 2017 until the beginning of 2020. The steep symmetric light curve as well as its shape look very different with respect to normal random variations in AGN. However, the observations could be easily reproduced with a best-fit standard microlensing light curve. Based on the light curve, we derived a characteristic timescale of 155 days. During the outburst, the continuum as well as the broad line intensities varied; however, the narrow emission lines did not. This is an indication that the lensing object orbits the AGN nucleus between the broad line region (BLR) and the narrow line region (NLR), that is, at a distance on the order of 250 light days. The light curve can be reproduced by a lens with a BH mass of 5*10^{5} M_solar. This implies a mass ratio to the central AGN on the order of 1 to 10.

Current images of the supermassive black hole (SMBH) candidates at the center of our Galaxy and M87 have opened an unprecedented era for studying strong gravity and the nature of relativistic sources. Very-long-baseline interferometry (VLBI) data show images consistent with a central SMBH within General Relativity (GR). However, it is essential to consider whether other well-motivated dark compact objects within GR could produce similar images. Recent studies have shown that dark matter (DM) halos modeled as self-gravitating systems of neutral fermions can harbor very dense fermionic cores at their centers, which can mimic the spacetime features of a black hole (BH). Such dense, horizonless DM cores can satisfy the observational constraints: they can be supermassive and compact and lack a hard surface. We investigate whether such cores can produce similar observational signatures to those of BHs when illuminated by an accretion disk. We compute images and spectra of the fermion cores with a general-relativistic ray tracing technique, assuming the radiation originates from standard $\alpha$ disks, which are self-consistently solved within the current DM framework. Our simulated images possess a central brightness depression surrounded by a ring-like feature, resembling what is expected in the BH scenario. For Milky Way-like halos, the central brightness depressions have diameters down to $\sim 35\, \mu$as as measured from a distance of approximately $8\,$kpc. Finally, we show that the DM cores do not possess photon rings, a key difference from the BH paradigm, which could help discriminate between the models.

Despite having different astronomical characteristics, the studies of mira variables and ultra-cool dwarfs frequently show similar red colors, which could cause leading to photometric misclassification. This study uses photometric data from the WISE, 2MASS, and Pan-STARRS surveys to construct color-based selection criteria for red dwarfs, brown dwarfs, and Mira variables. On analyzing the color indices, we developed empirical rules that separate these objects with an overall classification accuracy of approximately 91%-92%. While the differentiation between red dwarfs and both Mira variables and brown dwarfs is effective, challenges remain in distinguishing Mira variables from brown dwarfs due to overlapping color indices. The robustness of our classification technique was validated by a bootstrap analysis, highlighting the significance of color indices in large photometric surveys for stellar classification.

The Spectro-Polarimeter (SP) is a new instrument installed at the upgraded Andrei B. Severny Solar Tower Telescope (STT) at the Crimean Astrophysical Observatory. The instrument is a traditional echelle slit dual-beam spectropolarimeter with temporal modulation of the polarization. STT-SP provides simultaneous spectropolarimetric observations of the Sun within three 15 Angstrom wide spectral ranges around photospheric Fe I 5250, Fe I 5324, and chromospheric Mg I b2 5172 spectral lines. The spectral resolution of the instrument reaches 70,000 with the seeing-constrained slit width of 1 arcsec. The field-of-view of STT-SP is 200 arcsec allowing one to map a moderate size active region within a single raster scan. The instrument will provide new opportunities in the analysis of magnetic fields and thermodynamics of the lower atmosphere of the Sun. In this paper we describe the optical design of STT-SP and present the preliminary results acquired during the commissioning of the instrument.

We study the observable spectral and temporal properties of kilonova remnants analytically, and showcase quantitative differences with respect to supernova remnants. We provide detection prospects of kilonova remnants in the context of ongoing radio surveys. We find that there is a good chance to expect 10s of such objects in future surveys with a flux threshold of $\sim 0.1$ mJy. Kilonova remnants from a postulated population of long lived supramassive neutron star remnants of neutron star mergers are even more likely to be detected as they are extremely bright and peak earlier. For ongoing survey with threshold of $\sim$ mJy, we expect to find 10-100s of such objects if they are a significant fraction of total kilonova population. Considering that there are no promising such kilonovae candidates in current surveys, we constrain the fraction of such extreme kilonova to be no more than 30 percent of the overall kilonovae rate, depending on the details of ejecta mass and external density distribution.

The SST-1M is a Small-Sized Telescope (SST) designed to provide a cost-effective and high-performance solution for gamma-ray astrophysics, particularly for energies beyond a few TeV. The goal is to integrate this telescope into an array of similar instruments, leveraging its lightweight design, earthquake resistance, and established Davies-Cotton configuration. Additionally, its optical system is designed to function without a protective dome, allowing it to withstand the harsh atmospheric conditions typical of mountain environments above 2000 m. The SST-1M utilizes a fully digitizing camera system based on silicon photomultipliers (SiPMs). This camera is capable of digitizing all signals from the UV-optical light detectors, allowing for the implementation of various triggers and data analysis methods. We detail the process of designing, prototyping, and validating this system, ensuring that it meets the stringent requirements for gamma-ray detection and performance. An SST-1M stereo system is currently operational and collecting data at the Ondřejov observatory in the Czech Republic, situated at 500 m. Preliminary results from this system are promising. A forthcoming paper will provide a comprehensive analysis of the performance of the telescopes in detecting gamma rays and operating under real-world conditions.

We examine a century of radial velocity, visual magnitude, and astrometric observations of the nearest red supergiant, Betelgeuse, in order to reexamine the century-old assertion that Betelgeuse might be a spectroscopic binary. These data reveal Betelgeuse varying stochastically over years and decades due to its boiling, convective envelope, periodically with a $ 5.78$~yr long secondary period, and quasi-periodically from pulsations with periods of several hundred days. We show that the long secondary period is consistent between astrometric and RV datasets, and argue that it indicates a low-mass companion to Betelgeuse, less than a solar mass, orbiting in a 2,110 day period at a separation of just over twice Betelgeuse's radius. The companion star would be nearly twenty times less massive and a million times fainter than Betelgeuse, with similar effective temperature, effectively hiding it in plain sight near one of the best-studied stars in the night sky. The astrometric data favor an edge-on binary with orbital plane aligned with Betelgeuse's measured spin axis. Tidal spin-orbit interaction drains angular momentum from the orbit and spins up Betelgeuse, explaining the spin--orbit alignment and Betelgeuse's anomalously rapid spin. In the future, the orbit will decay until the companion is swallowed by Betelgeuse in the next 10,000 years.

Ultra-light primordial black holes (PBHs) with masses $M_{\rm PBH}<5\times 10^8\mathrm{g}$ can dominate transiently the energy budget of the Universe and reheat the Universe through their evaporation taking place before Big Bang Nucleosynthesis. The isocurvature energy density fluctuations associated to the inhomogeneous distribution of a population of such PBHs can induce an abundant production of GWs due to second-order gravitational effects. In this work, we discuss the effect of primordial non-Gaussianity on the clustering properties of PBHs and study the effect of a clustered PBH population on the spectral shape of the aforementioned induced GW signal. In particular, focusing on local-type non-Gaussianity we find a double-peaked GW signal with the amplitude of the low-frequency peak being proportional to the square of the non-Gaussian parameter $\tau_\mathrm{NL}$. Remarkably, depending on the PBH mass $M_{\rm PBH}$ and the initial abundance of PBHs at formation time, i.e. $\Omega_\mathrm{PBH,f}$, this double-peaked GW signal can lie well within the frequency bands of forthcoming GW detectors, namely LISA, ET, SKA and BBO, hence rendering this signal falsifiable by GW experiments and promoting it as a novel portal probing the primordial non-Gaussianity.

After the main sequence phase, stars more massive than 2.5 M$_\odot$ rapidly evolve through the Hertzsprung gap as yellow giants and supergiants (YSG), before settling into the red giant branch. Identifying YSG in nearby galaxies is crucial for pinpointing progenitors of luminous red novae (LRNe) - astrophysical transients attributed to stellar mergers. In the era of extensive transient surveys like the Vera Rubin Observatory's LSST, this approach offers a new way to predict and select common envelope transients. This study investigates potential progenitors and precursors of LRNe by analysing Hubble Space Telescope (HST) photometry of stellar populations in galaxies within 20 Mpc to identify YSG candidates. Additionally, we use ZTF and MeerLICHT/BlackGEM to identify possible precursors, preparing for future observations by the LSST. We compiled a sample of 369 galaxies with HST exposures in the F475W, F555W, F606W, and F814W filters. We identified YSG candidates using MESA stellar evolution tracks and statistical analysis of color-magnitude diagrams (CMDs). Our sample includes 246,573 YSG candidates with masses between 3 and 20 $M_\odot$ and is affected by various contaminants, such as foreground stars and extinguished main-sequence stars. After excluding foreground stars using Gaia proper motions, contamination is estimated at 1.7\% from foreground stars and 20\% from extinction affecting main-sequence stars. Combining our YSG candidates with time-domain catalogs yielded several interesting candidates. Notably, we identified 12 LRN precursor candidates for which followup is encouraged. We highlight the importance of monitoring future transients that match YSG candidates to avoid missing potential LRNe and other rare transients. LSST will be a game changer in the search for LRN progenitors and precursors, discovering over 300,000 new YSG and 100 precursors within 20 Mpc.

We revisit a quiet 14-day period of solar minimum during January 2008 and track sub-streamer propagating disturbances (PDs) from low heights in STEREO/EUVI to the extended corona through STEREO/COR1 and into STEREO/COR2 along nonradial paths that trace the structure of the underlying streamers. Using our recently developed method for generating nonradial Height-Time profiles of outward PDs (OPDs) and inward PDs (IPDs), we obtained their velocities along the radial and position angle directions. Our analysis of 417 unique OPDs revealed two classes: slow and fast OPDs. Slow OPDs form preferentially at $\approx$1.6 $R_\odot$ closer to the streamer boundaries, with asymmetric occurrence rates, and show speeds of $16.4_{-8.4}^{+26.6}km/s$ at 1.5 $R_\odot$ and accelerate up to $200.1_{-57.9}^{+71.1}km/s$ at 7.5 $R_\odot$. Fast OPDs form preferentially at $\approx$ 1.6 $R_\odot$ and at $\approx$3.0 $R_\odot$ both at the streamer boundaries and slightly more often within them. They show speeds of $87.8_{-24.8}^{+59.1}km/s$ at 1.5 $R_\odot$ up to $197.8_{-46.7}^{+61.8}km/s$ at 7.5 $R_\odot$. IPDs are observed forming at $\approx$1.8 $R_\odot$ with speeds of tens of $km/s$, mostly concentrated in the aftermath of a CME eruption. We present an example in which we show that periodic brightness variations related to OPDs remained in the range of 98 to 128 min, down to $\approx$2.0 $R_\odot$, well within the field of view of COR1. The velocity profiles of slow OPDs for heliocentric height below 3.0 $R_\odot$ show good agreement with speeds more closely related to the bulk solar wind obtained via interplanetary scintillation.

Stars play a key role in the evolution of the Universe, as sources of radiation, as dynamical engines, and as chemical factories. Outputs of stellar models are then central to various studies in astrophysics. Stellar physics links fundamental physical aspects to hydrodynamic and magnetohydrodynamic processes, and the validity of stellar models depends directly on the modelling of these complex mechanisms. We describe here the different transport processes at work in stellar interiors and how the modelling of these processes can be improved thanks to the unique ability of asteroseismology, the study of stellar oscillations, to probe the internal structure and dynamics of stars.

Determining the Hubble constant (H0), a fundamental parameter describing cosmic expansion, remains a challenge due to conflicting measurements from the early and late universe. Gravitational wave (GW) observations from binary neutron star (BNS) mergers, with identified host galaxies through electromagnetic (EM) follow-up, offer an independent method to measure H0. However, this requires detection of numerous events, which could take decades with current GW detectors. LIGO-India can dramatically accelerate this effort. With sensitivity comparable to the existing LIGO detectors, its addition to the LIGO-Virgo network could increase detected events by 70%. This improvement nearly doubles when accounting for the detector's 70% duty cycle, increasing the probability of simultaneous operation of three detectors by a factor of ~2. We perform end-to-end simulations to estimate triple-coincidence detection rates and sky localization, considering realistic BNS populations, lightcurves, and EM observatory specifications. Our findings suggest LIGO-India could increase BNS events with observed kilonovae by ~2-7 times. The factor of few improvements in source localization precision with LIGO-India can allow much deeper EM follow-up campaigns (not considered in the simulations), potentially increasing the overall rate of detection of EM counterparts by a factor of ~20, which can have an enormous impact in addressing critical questions in different areas of astronomy. We evaluate the impact of LIGO-India in the context of H0 measurement and argue that it can cut down the required observation time of several decades by a factor of few and possibly to just few years with regular sensitivity upgrades.

Recent mid-infrared observations with JWST/MIRI have resulted in the first direct detections of absorption features from silicate clouds in the transmission spectra of two transiting exoplanets, WASP-17 b and WASP-107 b. In this paper, we measure the mid-infrared ($5-12$ $\mu$m) dayside emission spectrum of the benchmark hot Jupiter HD 189733 b with MIRI LRS by combining data from two secondary eclipse observations. We confirm the previous detection of H$_2$O absorption at 6.5 $\mu$m from Spitzer/IRS and additionally detect H$_2$S as well as an absorption feature at 8.7 $\mu$m in both secondary eclipse observations. The excess absorption at 8.7 $\mu$m can be explained by the presence of small ($\sim$0.01 $\mu$m) grains of SiO$_2$[s] in the uppermost layers of HD 189733 b's dayside atmosphere. This is the first direct detection of silicate clouds in HD 189733 b's atmosphere, and the first detection of a distinct absorption feature from silicate clouds on the day side of any hot Jupiter. We find that models including SiO$_2$[s] are preferred by $6-7\sigma$ over clear models and those with other potential cloud species. The high altitude location of these silicate particles is best explained by formation in the hottest regions of HD 189733 b's dayside atmosphere near the substellar point. We additionally find that HD 189733 b's emission spectrum longward of 9 $\mu$m displays residual features not well captured by our current atmospheric models. When combined with other JWST observations of HD 189733 b's transmission and emission spectrum at shorter wavelengths, these observations will provide us with the most detailed picture to date of the atmospheric composition and cloud properties of this benchmark hot Jupiter.

Recent reports of cosmological parity violation in the 4PCF raises the question of how such violations could be systematically generated. Here we present a constructive procedure to generate arbitrary violations of vectorial and tensorial types on any scale, which is computationally efficient in the squeezed limit. We directly compute their numerical transfer function, and find strong conservation in the linear regime. This procedure spans all squeezed parity violating observables at the 4PCF, following the quadratic estimator classification.

We develop a transformer-based conditional generative model for discrete point objects and their properties. We use it to build a model for populating cosmological simulations with gravitationally collapsed structures called dark matter halos. Specifically, we condition our model with dark matter distribution obtained from fast, approximate simulations to recover the correct three-dimensional positions and masses of individual halos. This leads to a first model that can recover the statistical properties of the halos at small scales to better than 3% level using an accelerated dark matter simulation. This trained model can then be applied to simulations with significantly larger volumes which would otherwise be computationally prohibitive with traditional simulations, and also provides a crucial missing link in making end-to-end differentiable cosmological simulations. The code, named GOTHAM (Generative cOnditional Transformer for Halo's Auto-regressive Modeling) is publicly available at \url{this https URL}.