Papers for Thursday, Jun 20 2024

Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the LVK transient catalogue. For our COMPAS models, we calculate the likelihood with and without the dependence on the predicted number of BBH merger events. We find strong correlations between hyper-parameters governing the specific angular momentum (AM) of mass lost during mass transfer, the mass-loss rates of Wolf--Rayet stars via winds and the chemically homogeneous evolution (CHE) formation channel. We conclude that analysing the marginalised and unmarginalised likelihood is a good indicator of whether the population parameters distribution and number of observed events reflect the LVK data. In doing so, we see that the majority of the models preferred in terms of the population-level parameters of the BBHs greatly overpredict the number of events we should have observed to date. Looking at the smaller number of models which perform well with both likelihoods, we find that those with no CHE, AM loss occurring closer to the donor during the first mass-transfer event, and/or higher rates of mass-loss from Wolf--Rayet winds are generally preferred by current data. We find these conclusions to be robust to our choice of selection criteria.

We investigate the effects of magnetic diffusion on the spectrum of ultra-high energy cosmic rays (UHECRs) from a cosmological perspective. To this end, we consider two modified theories of gravity (MTGs), namely, the $f(R)$ gravity and a symmetric teleparallel gravity, also known as $f(Q)$ gravity. Utilizing these two MTGs, we calculate the suppression in the flux of UHECRs for a collection of sources. Non-evolution (NE) and cosmic star formation rate (SFR) scenarios have been considered in our calculation of the suppression factor. This study also includes a mixed composition scenario involving the nuclei upto iron (Fe). Furthermore, we provide a parameterization of the suppression factor for the proton and also for the mixed compositions within the $f(R)$ and $f(Q)$ theories, considering both NE and SFR scenarios. The influence of the turbulent magnetic field on the suppression factor is also incorporated in our work. Comparative analysis of all our results with the standard $\Lambda$CDM model reveals significant effects of MTGs on the suppression factor that the $f(R)$ power-law model predicts the lowest suppression factor, while the $f(Q)$ model predicts the highest, and interestingly the results from the standard model fall within the range predicted by these two cosmological models.

We select blue-horizontal branch stars (BHBs) from the internal data release of the Hyper Suprime-Cam Subaru Strategic Program to reveal the global structure of the Milky Way (MW) stellar halo. The data are distributed over $\sim 1,100$~deg$^2$ area in the range of $18.5<g<24.5$~mag, so that candidate BHBs are detectable over a Galactocentric radius of $r \simeq 36-575$~kpc. In order to select most likely BHBs by removing blue straggler stars and other contaminants in a statistically significant manner, we develop and apply an extensive Bayesian method, as described in \citet{Fukushima2019}. Our sample can be fitted to either a single power-law profile with an index of $\alpha=4.11^{+0.18}_{-0.18}$ or a broken power-law profile with an index of $\alpha_{\rm in}=3.90^{+0.24}_{-0.30}$ at $r$ below a broken radius of $r_{\rm b}=184^{+118}_{-66}$ kpc and a very steep slope of $\alpha_{\rm out}=9.1^{+6.8}_{-3.6}$ at $r>r_{\rm b}$; the statistical difference between these fitting profiles is small. Both profiles are found to show prolate shapes having axial ratios of $q=1.47^{+0.30}_{-0.33}$ and $1.56^{+0.34}_{-0.23}$, respectively. We also find a signature of the so-called "splashback radius" for the candidate BHBs, which can reach as large as $r \sim 575$~kpc, although it is still inconclusive owing to rather large distance errors in this faintest end of the sample. Our results suggest that the MW stellar halo consists of the two overlapping components: the {\it in situ} inner halo showing a relatively steep radial density profile and the {\it ex situ} outer halo with a shallower profile, being characteristic of a component formed from accretion of small stellar systems.

We revisit the role of primordial black holes (PBHs) as potential dark matter (DM) candidates, particularly focusing on light asteroid-mass PBHs. These PBHs are expected to emit particles through Hawking evaporation that can generate cosmic rays (CRs), eventually producing other secondary radiations through their propagation in the Milky Way, in addition to prompt emissions. Here, we perform a comprehensive analysis of CR signals resulting from PBH evaporation, incorporating the full CR transport to account for reacceleration and diffusion effects within the Milky Way. In particular, we revisit the $e^\pm$ flux produced by PBHs, using Voyager 1, and study for the first time the diffuse X-ray emission from the up-scattering of Galactic ambient photons due to PBH-produced $e^\pm$ via the inverse Compton effect using XMM-Newton data, as well as the morphological information of the diffuse 511 keV line measured by INTEGRAL/SPI. In doing so, we provide leading constraints on the fraction of DM that can be in form of PBHs in a conservative way, whilst also testing how different assumptions on spin and mass distributions affect our conclusions.

Galaxy clusters are one of the most powerful probes to study extensions of General Relativity and the Standard Cosmological Model. Upcoming surveys like the Vera Rubin Observatory's Legacy Survey of Space and Time are expected to revolutionise the field, by enabling the analysis of cluster samples of unprecedented size and quality. To reach this era of high-precision cluster cosmology, the mitigation of sources of systematic error is crucial. A particularly important challenge is bias in cluster mass measurements induced by inaccurate photometric redshift estimates of source galaxies. This work proposes a method to optimise the source sample selection in cluster weak lensing analyses drawn from wide-field survey lensing catalogs to reduce the bias on reconstructed cluster masses. We use a combinatorial optimisation scheme and methods from variational inference to select galaxies in latent space to produce a probabilistic galaxy source sample catalog for highly accurate cluster mass estimation. We show that our method reduces the critical surface mass density $\Sigma_{\rm crit}$ modelling bias on the 60-70% level, while maintaining up to 90% of galaxies. We highlight that our methodology has applications beyond cluster mass estimation as an approach to jointly combine galaxy selection and model inference under sources of systematics.

Pulsar timing arrays (PTAs) have made tremendous progress and are now showing strong evidence for the gravitational-wave background (GWB). Further probing the origin and characteristics of the GWB will require more generalized analysis techniques. Bayesian methods are most often used but can be computationally expensive. On the other hand, frequentist methods, like the PTA Optimal Statistic (OS), are more computationally efficient and can produce results that are complementary to Bayesian methods, allowing for stronger statistical cases to be built from a confluence of different approaches. In this work we expand the capabilities of the OS through a technique we call the Per-Frequency Optimal Statistic (PFOS). The PFOS removes the underlying power-law assumption inherent in previous implementations of the OS, and allows one to estimate the GWB spectrum in a frequency-by-frequency manner. We have also adapted a recent generalization from the OS pipeline into the PFOS, making it capable of accurately characterizing the spectrum in the intermediate and strong GW signal regimes using only a small fraction of the necessary computational resources when compared with fully-correlated Bayesian methods, while also empowering many new types of analyses not possible before. We find that even in the strong GW signal regime, where the GWB dominates over noise in all frequencies, the injected value of the signal lies within the 50th-percentile of the PFOS uncertainty distribution in 41-45% of simulations, remaining 3$\sigma$-consistent with unbiased estimation.

It is well established that maximizing the information extracted from upcoming and ongoing stage-IV weak-lensing surveys requires higher-order summary statistics that complement the standard two-point statistics. In this work, we focus on weak-lensing peak statistics to test two popular modified gravity models, $f(R)$ and nDGP, using the FORGE and BRIDGE weak-lensing simulations, respectively. From these simulations we measure the peak statistics as a function of both cosmological and modified gravity parameters simultaneously. Our findings indicate that the peak abundance is sensitive to the strength of modified gravity, while the peak two-point correlation function is sensitive to the nature of the screening mechanism in a modified gravity model. We combine these simulated statistics with a Gaussian Process Regression emulator and a Gaussian likelihood to generate stage-IV forecast posterior distributions for the modified gravity models. We demonstrate that, in a systematics-free case, peak statistics can constrain $\log_{10}(f_{R0}) = -6$ to 2% precision, and $\log_{10}(H_0 r_c) = 0.5$ to 25% precision. Finally, we find that our weakest models, $\log_{10}(f_{R0})$ < 6.17 and $\log_{10}(H_0 r_c)$ > 1, can be ruled out at the two sigma level with an area of $300 \rm{deg}^2$ and $1000 \rm{deg}^2$, respectively, with the upcoming stage-IV data.

Recent observations with JWST and ALMA have unveiled galaxies with regular discs at significantly higher redshifts than previously expected. This appears to be in contrast with constraints on the stellar populations of the Milky Way, suggesting that the bulk of the Galactic thin disc formed after $z=1$, and raises questions about the history, evolution, and survivability of primordial discs. Here, we use GigaEris, a state-of-the-art $N$-body, hydrodynamical, cosmological ``zoom-in'' simulation with a billion particles within the virial radius, to delve into the formation of the early kinematically cold discs (KCDs), defined by their ratio between the mean rotational velocity and the radial velocity dispersion, of a Milky Way-sized galaxy at redshifts $z\gtrsim 4$. Our analysis reveals a primarily inward migration pattern for disc stars formed at $z \gtrsim 6$, turning into a mix of inward and outward migration at later times. Stars migrating outwards undergo minimal kinematic heating, and might be identified as part of the thin disc forming at much later epochs. We find that approximately 76 per cent of all stars formed in the KCD at $z \sim 7$ become part of a pseudo-bulge by $z = 4.4$. This proportion decreases to below 10 per cent for KCD stars formed at $z \lesssim 5$. The inward migration of stars born in our KCDs at $z \gtrsim 4$ deviates from the expected inside-out formation scenario of thin discs at lower redshifts. Our results suggest a novel ``two-phase'' disc formation process, whereby the early disc transforms primarily into the pseudo-bulge within less than a billion years, whereas the present-day disc forms subsequently from higher-angular momentum material accreted at later times.

Given their dominance of the galaxy number density, dwarf galaxies are central to our understanding of galaxy formation. While the incidence of AGN and their impact on galaxy evolution has been extensively studied in massive galaxies, much less is known about the role of AGN in the evolution of dwarfs. We search for radiatively-efficient AGN in the nearby (0.1 < z < 0.3) dwarf (10^8 MSun < M < 10^10 MSun) population, using SED fitting (via Prospector) applied to deep ultraviolet to mid-infrared photometry of 508 dwarf galaxies. Around a third (32 +/- 2 per cent) of our dwarfs show signs of AGN activity. We compare the properties of our dwarf AGN to control samples, constructed from non-AGN, which have the same distributions of redshift and stellar mass as their AGN counterparts. KS tests between the AGN and control distributions indicates that the AGN do not show differences in their distances to nodes, filaments and nearby massive galaxies from their control counterparts. This indicates that AGN triggering in the dwarf regime is not strongly correlated with local environment. The fraction of AGN hosts with early-type morphology and those that are interacting are also indistinguishable from the controls within the uncertainties, suggesting that interactions do not play a significant role in inducing AGN activity in our sample. Finally, the star formation activity in dwarf AGN is only slightly lower than that in their control counterparts, suggesting that the presence of radiatively-efficient AGN does not lead to significant, prompt quenching of star formation in these systems.

Multi-messenger observations of coalescing binary neutron stars (BNSs) are a direct probe of the expansion history of the universe and carry the potential to shed light on the disparity between low- and high-redshift measurements of the Hubble constant $H_0$. To measure the value of $H_0$ with such observations requires pristine inference of the luminosity distance and the true source redshift with minimal impact from systematics. In this analysis, we carry out joint inference on mock gravitational wave (GW) signals and their electromagnetic (EM) afterglows from BNS coalescences and find that the inclination angle inferred from the afterglow light curve and apparent superluminal motion can be precise, but need not be accurate and is subject to systematic uncertainty that could be as large as $1.5\sigma$. This produces a disparity between the EM and GW inferred inclination angles, which if not carefully treated when combining observations can bias the inferred value of $H_0$. We also find that already small misalignments of $3^{\circ}-6^{\circ}$ between the inherent system inclinations for the GW and EM emission can bias the inference by $\mathcal{O}(1-2\sigma)$ if not taken into account. As multi-messenger BNS observations are rare, we must make the most out of a small number of events and harness the increased precision, while avoiding reduced accuracy. We demonstrate how to mitigate these potential sources of bias by jointly inferring the mismatch between the GW- and EM-based inclination angles and $H_0$.

We present the first step towards deriving cosmological constraints through the abundances of galaxy clusters selected in a $510\,\mathrm{deg}^2$ weak-lensing aperture mass map, constructed with the Year-Three shear catalog from the Hyper Suprime-Cam Subaru Strategic Program. We adopt a conservative source galaxy selection to construct a sample of $129$ weak-lensing peaks with a signal-to-noise ratio above $4.7$. We use semi-analytical injection simulations to derive the selection function and the mass--observable relation of our sample. These results take into account complicated uncertainties associated with weak-lensing measurements, such as the non-uniform survey depth and the complex survey geometry, projection effects from uncorrelated large-scale structures, and the intrinsic alignment of source galaxies. We also propose a novel modeling framework to make parts of the mass--observable relation insensitive to assumed cosmological parameters. Such a framework not only offers a great computational advantage to cosmological studies, but can also benefit future astrophysical studies using shear-selected clusters. Our results are an important step towards utilizing these cluster samples that are constructed nearly independent of any baryonic assumptions in upcoming deep-and-wide lensing surveys from the Vera Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $\lambda\lambda$4960,5008 outflow signatures in quasars with and without Broad Absorption Lines (BALs), aiming to test how the broad absorption troughs seen in the rest-frame ultraviolet are linked to the narrow line region outflows seen in the rest-frame optical. We present new near-infrared spectra from Magellan/FIRE which cover [OIII] in 12 quasars with 2.1 < z < 2.3, selected to have strong outflow signatures in CIV $\lambda$1550. Combining with data from the literature, we build a sample of 73 BAL, 115 miniBAL and 125 non-BAL QSOs with 1.5 < z < 2.6. The strength and velocity width of [OIII] correlate strongly with the CIV emission properties, but no significant difference is seen in the [OIII] emission-line properties between the BALs, non-BALs and miniBALs once the dependence on CIV emission is taken into account. A weak correlation is observed between the velocities of CIV BALs and [OIII] emission, which is accounted for by the fact that both outflow signatures correlate with the underlying CIV emission properties. Our results add to the growing evidence that BALs and non-BALs are drawn from the same parent population and are consistent with a scenario wherein BAL troughs are intermittent tracers of persistent quasar outflows, with a part of such outflow becoming optically thick along our line-of-sight for sporadic periods of time within which BALs are observed.

We present cosmological constraints using the abundance of weak-lensing shear-selected galaxy clusters in the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The clusters are selected on the mass maps constructed using the three-year (Y3) weak-lensing data with an area of $\approx500~$deg$^2$, resulting in a sample size of $129$ clusters with high signal-to-noise ratios $\nu$ of $\nu\geq4.7$. Owing to the deep, wide-field, and uniform imaging of the HSC survey, this is by far the largest sample of shear-selected clusters, in which the selection solely depends on gravity and is free from any assumptions about the dynamical state. Informed by the optical counterparts, the shear-selected clusters span a redshift range of $z\lesssim0.7$ with a median of $z\approx0.3$. The lensing sources are securely selected at $z\gtrsim0.7$ with a median of $z\approx1.3$, leading to nearly zero cluster member contamination. We carefully account for (1) the bias in the photometric redshift of sources, (2) the bias and scatter in the weak-lensing mass using a simulation-based calibration, and (3) the measurement uncertainty that is directly estimated on the mass maps using an injection-based method developed in a companion paper (Chen et al. submitted). In a blind analysis, the fully marginalized posteriors of the cosmological parameters are obtained as $\Omega_{\mathrm{m}} = 0.50^{+0.28}_{-0.24}$, $\sigma_8 = 0.685^{+0.161}_{-0.088}$, $\hat{S}_{8}\equiv\sigma_8\left(\Omega_{\mathrm{m}}/0.3\right)^{0.25} = 0.835^{+0.041}_{-0.044}$, and $\sigma_8\sqrt{\Omega_{\mathrm{m}}/0.3} = 0.993^{+0.084}_{-0.126}$ in a flat $\Lambda$CDM model. We compare our cosmological constraints with other studies, including those based on cluster abundances, galaxy-galaxy lensing and clustering, and Cosmic Microwave Background observed by $Planck$, and find good agreement at levels of $\lesssim2\sigma$. [abridged]

Extreme coronal line emitters (ECLEs) are objects showing transient high-ionization lines in the centers of galaxies. They have been attributed to echoes of high-energy flares of ionizing radiation, such as those produced by tidal disruption events (TDEs), but have only recently been observed within hundreds of days after an optical transient was detected. AT 2022upj is a nuclear UV-optical flare at z=0.054 with spectra showing [Fe X] {\lambda}6375 and [Fe XIV] {\lambda}5303 during the optical peak, the earliest presence of extreme coronal lines during an ongoing transient. AT 2022upj is also the second ever ECLE (and first with a concurrent flare) to show broad He II {\lambda}4686 emission, a key signature of optical/UV TDEs. We also detect X-ray emission during the optical transient phase, which may be related to the source of ionizing photons for the extreme coronal lines. Finally, we analyze the spectroscopic evolution of each emission line and find that [Fe X] and [Fe XIV] weaken within 400d of optical peak, while [Fe VII] {\lambda}5720, [Fe VII] {\lambda}6087, and [O III] {\lambda}{\lambda}4959,5007 emerge over the same period. The velocities of the iron lines indicate circumnuclear gas within 0.1pc of the central supermassive black hole (SMBH), while a dust echo inferred from NEOWISE data indicates that circumnuclear dust lies at a minimum of 0.4pc away, providing evidence of stratified material around a SMBH. AT 2022upj is the first confirmed ECLE-TDE with clear signatures of both classes. This event's spectroscopic evolution on a $\sim$year unveils the impact of highly energetic flares such as TDEs on the complex environments around SMBHs.

The impact of viscosity in the Intracluster Medium (ICM) is still an open question in astrophysics. To address this problem, we have run a set of cosmological simulations of three galaxy clusters with a mass larger than $M_{\mathrm{Vir}} > 10^{15} $M$_{\odot}$ at $z=0$ using the SPMHD-code OpenGadget3. We aim to quantify the influence of viscosity and constrain its value in the ICM. Our results show significant morphological differences at small scales, temperature variations, and density fluctuations induced by viscosity. We observe a suppression of instabilities at small scales, resulting in a more filamentary structure and a larger amount of small structures due to the lack of mixing with the medium. The conversion from kinetic to internal energy leads to an increase of the virial temperature of the cluster of $\sim$5% - 10%, while the denser regions remain cold. The amplitude of density fluctuations is found to increase with viscosity, as well as the velocity fluctuations. However, comparison with observational data indicates consistency with observed density fluctuations, challenging the direct constraint of viscosity solely through density fluctuations. Furthermore, the ratio of density to velocity fluctuations remains close to 1 regardless of the amount of viscosity, in agreement with the theoretical expectations. Our results show for the first time in a cosmological simulation of a galaxy cluster the effect of viscosity in the ICM, a study that is currently missing in the literature.

Galactic cosmic rays (CRs) are of unknown origin, even though they have traditionally been connected to supernovae due to energetic arguments. In the last decades, Galactic black holes in X-ray binaries (BHXBs) have been proposed as candidate sources of CRs revising the CR paradigm. BHXBs launch two relativistic jets during outbursts, but recent observations suggest that such jets may be launched even during quiescence. A0620-00 is such a well-studied object that shows hints of jet emission. In this work, we study the simultaneous radio-to-X-ray spectrum detected from this source while in quiescence to better constrain the jet dynamics. Given that the majority of BHXBs spend their lifetime in quiescence (qBHXBs), we use the jet dynamics of A0620-00 to study a population of 100000 such sources distributed throughout the Galactic disc, and a further 10000 located in the Boxy Bulge around the Galactic centre. While the contribution to the CR spectrum is suppressed, we find that the cumulative intrinsic emission of qBHXBs from both the Boxy Bulge and the Galactic disc adds to the diffuse emission that various facilities detect from radio to TeV gamma rays. We examine the contribution of qBHXBs to the Galactic diffuse emission and investigate the possibility of SKA, INTEGRAL and CTAO to detect individual sources in the future. Finally, we compare the predicted neutrino flux to the recently presented Galactic diffuse neutrino emission by IceCube.

The galaxy SDSS1335+0728, which had exhibited no prior optical variations during the preceding two decades, began showing significant nuclear variability in the Zwicky Transient Facility (ZTF) alert stream from December 2019 (as ZTF19acnskyy). Its behaviour suggests that SDSS1335+0728 hosts a $\sim 10^6 M_{\odot}$ black hole (BH) that is currently in the process of `turning on'. We present a multi-wavelength photometric analysis and spectroscopic follow-up performed with the aim of better understanding the origin of the nuclear variations detected in SDSS1335+0728. We used archival photometry and spectroscopic data to study the state of SDSS1335+0728 prior to December 2019, and new observations from Swift, SOAR/Goodman, VLT/X-shooter, and Keck/LRIS taken after its turn-on to characterise its current state. We find that: (a) since 2021, the UV flux is four times brighter than the flux reported by GALEX in 2004; (b) since June 2022, the mid-infrared flux has risen more than two times, and the W1-W2 WISE colour has become redder; (c) since February 2024, the source has begun showing X-ray emission; (d) the narrow emission line ratios are now consistent with a more energetic ionising continuum; (e) broad emission lines are not detected; and (f) the [OIII] line increased its flux $\sim 3.6$ years after the first ZTF alert, which implies a relatively compact narrow-line-emitting region. We conclude that the variations observed in SDSS1335+0728 could be either explained by an AGN that is just turning on or by an exotic tidal disruption event (TDE). If the former is true, SDSS1335+0728 is one of the strongest cases of an AGN observed in the process of activating. If the latter, it would correspond to the longest and faintest TDE ever observed (or another class of still unknown nuclear transient). Future observations of SDSS1335+0728 are crucial to further understand its behaviour.

The Universal Rotation Curve (URC) of disk galaxies was originally proposed to predict the shape and amplitude of any rotation curve (RC) based solely on photometric data. Here, the URC is investigated with an extensive set of spatially-resolved rotation curves drawn from the PROBES-I, PROBES-II, and MaNGA data bases with matching multi-band surface brightness profiles from the DESI-LIS and WISE surveys for 3,846 disk galaxies. Common URC formulations fail to achieve an adequate level of accuracy to qualify as truly universal over fully sampled RCs. We develop neural network (NN) equivalents for the proposed URCs which predict RCs with higher accuracy, showing that URC inaccuracies are not due to insufficient data but rather non-optimal formulations or sampling effects. This conclusion remains even if the total RC sample is pruned for symmetry. The latest URC prescriptions and their NN equivalents trained on our sub-sample of 579 disk galaxies with symmetric RCs perform similarly to the URC/NN trained on the complete data sample. We conclude that a URC with an acceptable level of accuracy ($\Delta V_{\rm circ} \lesssim15$ per cent) at all radii would require a detailed modelling of a galaxy's central regions and outskirts (e.g., for baryonic effects leading to contraction or expansion of any dark-matter-only halo).

Stellar streams, long thin streams of stars, have been used as sensitive probes of dark matter substructure for over two decades. Gravitational interactions between dark matter substructures and streams lead to the formation of low density "gaps" in streams, with any given stream typically containing no more than a few such gaps. Prior models for the statistics of such gaps have relied on several simplifying assumptions for the properties of the subhalo population in the cold dark matter scenario. With the expected forthcoming increase in the number of streams, and gaps, observed, in this work we develop a more detailed model for the statistics of subhalos interacting with streams, and test some of the assumptions made in prior works. Instead of using simple fits to N-body estimates of subhalo population statistics at z = 0 as in previous work, we make use of realizations of time-dependent subhalo populations generated from a fully physical model, incorporating structure formation, and subhalo orbital evolution, including tidal heating and stripping physics, which has been carefully calibrated to match results of cosmological N-body simulations. We find that this model predicts up to 60% more gaps on average in Pal-5-like streams than prior works.

We simulate a variety of optical systematics for Taurus, a balloon-borne cosmic microwave background (CMB) polarisation experiment, to assess their impact on large-scale E-mode polarisation measurements and constraints of the optical depth to reionisation {\tau}. We model a one-month flight of Taurus from Wanaka, New Zealand aboard a super-pressure balloon (SPB). We simulate night-time scans of both the CMB and dust foregrounds in the 150GHz band, one of Taurus's four observing bands. We consider a variety of possible systematics that may affect Taurus's observations, including non-gaussian beams, pointing reconstruction error, and half-wave plate (HWP) non-idealities. For each of these, we evaluate the residual power in the difference between maps simulated with and without the systematic, and compare this to the expected signal level corresponding to Taurus's science goals. Our results indicate that most of the HWP-related systematics can be mitigated to be smaller than sample variance by calibrating with Planck's TT spectrum and using an achromatic HWP model, with a preference for five layers of sapphire to ensure good systematic control. However, additional beam characterization will be required to mitigate far-sidelobe pickup from dust on larger scales.

We investigate the physical properties of Lyman-alpha emitters (LAEs) and non-Lyman-alpha emitters (non-LAEs) at z$\sim$4.8--9.6 via a stacking analysis of 253 JWST/NIRSpec spectra of galaxies observed as part of the JWST Advanced Deep Extragalactic Survey (JADES). We identify a sample of 42 LAEs with the equivalent width of Ly$\alpha$ $\gtrsim$20Åand a sample of 211 non-LAEs, divide each sample further via the median redshift of the LAEs (z~6.3), and create composite spectra using the low and medium resolution spectra from NIRSpec. We estimate physical quantities such as dust extinction, UV continuum slope $\beta$, electron temperatures, ionization parameter, escape fraction of Ly$\alpha$ and Lyman Continuum, and the photon production rate for each bin/stack. The existing dust-extinction laws do not appear to be valid at these epochs. The emission line ratio analyses show that active galactic nuclei might dominate all sub-samples, irrespective of Ly$\alpha$ emission. LAEs show much higher [OIII]/[OII] and low [OII]/H$\delta$ at z$\lesssim$6.3 compared to non-LAEs, but these line ratios are not sufficient to distinguish the two populations at z$>$6.3. However, the LAEs samples show large EW([OIII]4959, 5007) ($>$1000Å) compared to the non-LAEs sample at all redshifts. CIV/Ly$\alpha$ and CIV/CIII] for LAE population at z$\lesssim$6.3 is $\sim$a factor of 5 larger than that for LAE population at z$>$6.3. The ionizing radiation for LAEs is hard, as revealed from several diagnostics, including CIV detection, high [OIII]/[OII] ($>$8), and large values of $\xi^{\star}_{ion}$.

A fraction of galaxy clusters harbor diffuse radio sources known as radio halos. The currently adopted scenario for their formation is based on second-order Fermi re-acceleration of seed electrons that is driven by merger-driven turbulence in the intra-cluster medium. This mechanism is expected to be inefficient, which implies that a significant fraction of halos should have very steep ($\alpha < -1.5$) energy spectra. We start investigating the potential and current limitations of the combination of the two surveys conducted by LOFAR, LoTSS (144 MHz) and LoLSS (54 MHz), to probe the origin of radio halos. We follow up the 20 radio halos detected in the DR1 of LoTSS, which covers the HETDEX field, with the LoLSS survey, and we study their spectral properties between 54 and 144 MHz. After the removal of compact sources, 9 halos were excluded due to unreliable halo flux density measurements at 54 MHz. Our main finding is that 7 out of 11 ($\sim$ 64%) exhibit an ultra-steep spectrum ($\alpha < -1.5$), which is a key prediction of turbulent re-acceleration models. We also note a tentative trend for more massive systems to host flatter halos, although the currently poor statistics does not allow for a deeper analysis. Our sample suffers from low angular resolution at 54 MHz, which limits the accuracy of the compact-sources subtraction. Nevertheless, this study is the first step towards providing compelling evidence for the existence of a large fraction of radio halos with very steep spectrum, which is a fundamental prediction of turbulent re-acceleration models. In this regard, the forthcoming second data release of LoLSS, along with the integration of LOFAR international stations and the instrumental upgrade to LOFAR2.0, will improve both the statistics and the low-frequency angular resolution, allowing to conclusively determine the origin of radio halos in galaxy clusters.

We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May-June. Companion multiwavelength data during the campaign are likewise shown. The 2-8 keV X-rays exhibit a net polarization degree PD=1.99%+/-0.13% (68% confidence). The polarization signal is found to increase with energy across IXPE's 2-8 keV bandpass. The polarized X-rays exhibit an energy-independent polarization angle of PA=-25.7+/-1.8 deg. East of North (68% confidence). This is consistent with being aligned to Cyg X-1's AU-scale compact radio jet and its pc-scale radio lobes. In comparison to earlier hard-state observations, the soft state exhibits a factor of 2 lower polarization degree, but a similar trend with energy and a similar (also energy-independent) position angle. When scaling by the natural unit of the disk temperature, we find the appearance of a consistent trendline in the polarization degree between soft and hard states. Our favored polarimetric model indicates Cyg X-1's spin is likely high (a* above ~0.96). The substantial X-ray polarization in Cyg X-1's soft state is most readily explained as resulting from a large portion of X-rays emitted from the disk returning and reflecting off the disk surface, generating a high polarization degree and a polarization direction parallel to the black hole spin axis and radio jet. In IXPE's bandpass, the polarization signal is dominated by the returning reflection emission. This constitutes polarimetric evidence for strong gravitational lensing of X-rays close to the black hole.

The star formation rate (SFR) is tightly connected to the amount of dense gas in molecular clouds. However, it is not fully understood how the relationship between dense molecular gas and star formation varies within galaxies and in different morphological environments. In this work, we study dense gas and star formation in the nearby spiral galaxy NGC 4321 to test how the amount of dense gas and its ability to form stars varies with environmental properties at 260 pc scales. We present new ALMA observations of HCN(1-0) line emission. Combined with existing CO(2-1) observations from ALMA, and H-alpha from MUSE, as well as F2100W from JWST to trace the SFR, we measure the HCN/CO line ratio, a proxy for the dense gas fraction and SFR/HCN, a proxy for the star formation efficiency of the dense gas. Towards the centre of the galaxy, HCN/CO systematically increases while SFR/HCN decreases, but these ratios stay roughly constant throughout the disc. Spiral arms, interarm regions, and bar ends show similar HCN/CO and SFR/HCN. On the bar, there is a significantly lower SFR/HCN at a similar HCN/CO. We conclude that the centres of galaxies show the strongest environmental influence on dense gas and star formation, suggesting either that clouds couple strongly to the surrounding pressure or that HCN is tracing more of the bulk molecular gas that is less efficiently converted into stars. On the contrary, across the disc of NGC 4321, where the ISM pressure is typically low, SFR/HCN does not show large variations (< 0.3 dex) in agreement with Galactic observations of molecular clouds. Despite the large variations across environments and physical conditions, HCN/CO is a good predictor of the mean molecular gas surface density at 260 pc scales.

The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter Array observations. These estimates are compared to atmospheric abundances of wide separation (> 10 au) gas giants. The resolved disk gas C/O ratios, from seven systems, generally exhibit C/O > 1 with sub-solar, or depleted, carbon content. In contrast, wide separation gas giants have atmospheric C/O ratios that cluster near or slightly above the presumed stellar value with a range of elemental C/H. From the existing disk composition, we infer that the solid-state mm/cm-sized pebbles have a total C/O ratio (solid cores and ices) that is solar (stellar) in content. We explore simple models that reconstruct the exoplanet atmospheric composition from the disk, while accounting for silicate cloud formation in the planet atmosphere. If wide separate planets formed via the core-accretion mechanism, they must acquire their metals from pebble or planetesimal accretion. Further, the dispersion in giant planet C/H content is best matched by a disk composition with modest and variable factors of carbon depletion. An origin of the wide separation gas giants via gravitational instability cannot be ruled out as stellar C/O ratios should natively form in this scenario. However, the variation in planet metallicity with a stellar C/O ratio potentially presents challenges to these models.

IRASF11119 is an ultra-luminous IR galaxy with post-merger morphology, hosting a type-1 QSO at z=0.189. Its 2013 Suzaku spectrum shows a prominent Ultra Fast Outflow (UFO) absorption feature (v_out~0.25c). In 2021, we obtained the first XMM-Newton long look of the target, coordinated with a simultaneous NuSTAR observation. The new high-quality data allow us to detect at P>99.8% c.l. multiple absorption features associated with the known UFO. Furthermore, an emission plus absorption feature at 1.1-1.3 keV reveals the presence of a blueshifted P-Cygni profile in the soft band. We associate the hard band features with blends of FeXXV and FeXXVI He$\alpha$-Ly$\alpha$ and He$\beta$-Ly$\beta$ line pairs and infer a large column (N$_H$~$10^{24}$ cm$^{-2}$) of highly ionized (log$\xi$~5) gas outflowing at v_out=0.27c. The 1 keV feature can be associated with a blend of Fe and Ne transitions, produced by a lower column (N$_H$~$10^{21}$ cm$^{-2}$) and ionization (log$\xi$~2.6) gas component outflowing at the same speed. Using a radiative-transfer disk wind model to fit the highly ionized UFO, we derive a large mass outflow rate, comparable with the mass accretion rate (M$_{out}$=4.25 M$_{Sun}$/yr, ~1.6 M$_{acc}$), and kinetic energy and momentum flux among the highest reported in the literature. We measure an extremely low high-energy cut-off (E$_c$~25 keV). Several other cases in the literature suggest that a steep X-ray continuum may be related to the formation of powerful winds. The lack of a significant momentum boost between the nuclear UFO and the different phases of the large-scale outflow, observed in IRASF11119 and in a growing number of sources with powerful UFOs, can be explained by (i) a momentum-driven expansion, (ii) an inefficient coupling of the UFO with the host ISM, or (iii) by repeated energy-driven expansion episodes with low duty-cycle, that average out on long time-scales.

We introduce a new approach for analyzing the IGM damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($1216\,Å \lesssim \lambda \lesssim 1280\,Å$) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone ($\lambda \lesssim 1216\,Å$), we construct a generative probabilistic model for the entire spectrum accounting for all sources of error - the stochasticity induced by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, the unknown intrinsic spectrum of the quasar, and spectral noise. This principled Bayesian method allows us to marginalize out nuisance parameters associated with the quasar's radiation and its unknown intrinsic spectrum to precisely measure the IGM neutral fraction, $\langle x_{\rm HI}\rangle$. A key element of our analysis is the use of dimensionality reduction (DR) to describe the intrinsic quasar spectrum via a small number of nuisance parameters. Using a large sample of $15,559$ SDSS/BOSS quasars at $z \gtrsim 2.15$ we trained and quantified the performance of six distinct DR methods, and find that a six parameter PCA model (five coefficients plus a normalization) performs best, with complex machine learning approaches providing no advantage. By conducting statistical inference on 100 realistic mock EoR quasar spectra, we demonstrate the reliability of the credibility contours that we obtain on $\langle x_{\rm HI}\rangle$ and the quasar lifetime, $t_{\rm Q}$. The new method introduced here will transform IGM damping wings into a precision probe of reionization, on the same solid methodological and statistical footing as other precision cosmological measurements.

We investigate the precision with which the Lyman-$\alpha$ damping wing signature imprinted on the spectra of high-redshift quasars (QSOs) by the foreground neutral intergalactic medium (IGM) can measure the history of cosmic reionization. We leverage a novel inference pipeline based on a generative probabilistic model for the entire spectrum (both red- and blueward of the Lyman-$\alpha$ line), accounting for all relevant sources of uncertainty - the stochasticity caused by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, it's unknown intrinsic spectrum, and spectral noise. Performing fast JAX-based Hamiltonian Monte-Carlo (HMC) parameter inference, we precisely measure the underlying global IGM neutral fraction as well as the lifetime of the quasar. Running a battery of tests on over a thousand mocks, we find optimal precision when running the pipeline with a six parameter PCA continuum model (five coefficients and a normalization) on $\mathrm{S}/\mathrm{N} \sim 10$ spectra, binned to a $\sim 500\,\mathrm{km}/\mathrm{s}$ velocity pixel scale, and extending at least out to the C IV $\lambda\,1549\,\text{Å}$ emission line. After marginalizing out nuisance parameters associated with the quasar continuum, a single spectrum constrains the IGM neutral fraction to $28.0_{-8.8}^{+8.2}\,\%$ and the quasar lifetime to $0.80_{-0.55}^{+0.22}\,\mathrm{dex}$, improving notably towards spectra with a stronger IGM damping wing imprint. Higher precision can be achieved by averaging over statistical quasar samples. We identify two primary sources of uncertainty that contribute approximately equally to the total error budget: the uncertain quasar continuum model and the stochastic distribution of neutral regions arising from both the reionization topology and the location of the quasar's ionization front.

We have carried out inversions of travel times as measured by Gizon et al. (2020) to infer the internal profile of the solar meridional circulation (MC). A linear inverse problem has been solved by the regularized least-squares method with a constraint that the angular momentum (AM) transport by MC should be equatorward (HK21-type constraint). Our motivation for using this constraint is based on the result by Hotta and Kusano (2021) where the solar equator-fast rotation was reproduced successfully without any manipulation. The inversion result indicates that the MC profile is a double-cell structure if the so-called HK21 regime, in which AM transported by MC sustains the equator-fast rotation, correctly describes the physics inside the solar convective zone. The sum of the squared residuals computed with the inferred double-cell MC profile is comparable to that computed with the single-cell MC profile obtained when we exclude the HK21-type constraint, showing that both profiles can explain the data more or less at the same level. However, we also find that adding the HK21-type constraint degrades the resolution of the averaging kernels. Although it is difficult for us to determine the large-scale morphology of the solar MC at the moment, our attempt highlights the relevance of investigating the solar MC profile from both theoretical and observational perspectives.

The phenomenon of changing-look (CL) behavior in active galactic nuclei (AGN) is characterized by dramatic changes in luminosity and/or emission line profiles over relatively short periods, ranging from months to years. The origin of CL-AGNs remains a mystery, but one proposed explanation involves the response of the inner AGN disk to tidal disruption events (TDEs) around the supermassive black hole (SMBH). In this Letter, we calculate the predicted frequency of AGN TDEs as a function of SMBH mass and compare the results to the observed CL-AGN distribution. We find that if the fraction of CL-AGNs caused by AGN-TDEs is high, then: (1) most SMBHs in CL-AGN are near maximal spin, with the dimensionless spin parameter $a>0.9$; (2) AGN inner disks have a high surface density ($\geq 10^{7}\, {\rm g\, cm^{-2}}$); (3) typical AGN lifetimes are $\sim 10$-$100$ Myr; and (4) a nuclear star cluster initial mass function (IMF) that scales as $\sim m_*^{-1.6}$ is preferred. Future observations of CL-AGN will help constrain the fraction of CL-AGNs caused by AGN-TDEs, SMBH spins, AGN lifetimes, and the nuclear star cluster IMF.

Observations of the cosmic microwave backgroundradiation are described to remarkable accuracy by the six-parameterLambda CDM cosmology. However, the key ingredients of this model, namely dark matter, dark energy and cosmic inflation are not understood at a fundamental level. It is therefore important to investigate tensions between the CMB and other cosmological probes. I will review aspects of tensions with direct measurements of the Hubble constant H_0, measurements of weak gravitational lensing, and the recent hints of evolving dark energy reported by the Dark Energy Spectroscopic Instrument (DESI) collaboration.

Axion quark nuggets (AQNs) are hypothetical objects with a mass greater than a few grams and sub-micrometer size, formed during the quark-hadron transition. Originating from the axion field, they offer a possible resolution of the similarity between visible and dark components of the Universe. These composite objects behave as cold dark matter, interacting with ordinary matter and resulting in pervasive electromagnetic radiation throughout the Universe. This work aims to predict the electromagnetic signature in large-scale structures from the AQN-baryon interaction, accounting for thermal and non-thermal radiations. We use Magneticum hydrodynamical simulations to describe the distribution and dynamics of gas and dark matter at cosmological scales. We calculate the electromagnetic signature from radio, starting at $\nu \sim$ 1 GHz, up to a few keV X-ray energies. We find that the AQNs signature is characterized by monopole and fluctuation signals. The amplitude of both signals strongly depends on the average AQN mass and the ionization level of the baryonic environment. We identify a most optimistic scenario with a signal often near the sensitivity limit of existing instruments, such as FIRAS and the South Pole Telescope for high-resolution. Fluctuations in the Extra-galactic Background Light caused by the AQN can be tested with space-based imagers Euclid and James Webb Space Telescope. We also identify a minimal configuration, still out of reach of existing instruments, but future experiments might be able to pose constraints on the AQN model. We conclude that this is a viable dark matter model, which does not violate the canons of cosmology, nor existing observations. The best chances for testing this model reside in 1) ultra-deep IR and optical surveys, 2) spectral distorsions of the CMB and 3) low-frequency (1 GHz < $\nu $ < 100 GHz) and high-resolution ($\ell > 10^4$) observations.

Zw~049.057 is a moderate mass, dusty, early-type galaxy that hosts a powerful compact obscured nucleus (CON, L$_{FIR,CON} \geq$10$^{11}$~L$_{\odot}$). The resolution of HST enabled measurements of the stellar light distribution and characterization of dust features. Zw~049.057 is inclined with a prominent three zone disk; the R$\approx$ 1kpc star forming inner dusty disk contains molecular gas, a main disk with less dust and an older stellar population, and a newly detected outer stellar region at R$>$6~kpc with circular isophotes. Previously unknown polar dust lanes are signatures of a past minor merger that could have warped the outer disk to near face-on. Dust transmission measurements provide lower limit gas mass estimates for dust features. An extended region with moderate optical depth and M$\geq$ 2$\times$10$^8$~M$_{\odot}$ obscures the central 2~kpc. Optical spectra show strong interstellar Na~D absorption with a constant velocity across the main disk, likely arising in this extraplanar medium. Opacity measurements of the two linear dust features, pillars, give a total mass of $\geq$10$^6$~M$_{\odot}$, flow rates of $\geq$2~M$_{\odot}$~yr$^{-1}$, and few Myr flow times. Dust pillars are associated with the CON and are visible signs of its role in driving large-scale feedback. Our assessments of feedback processes suggest gas recycling sustains the CON. However, radiation pressure driven mass loss and efficient star formation must be avoided for the AGN to retain sufficient gas over its lifespan to produce substantial mass growth of the central black hole.

We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $\Gamma$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overline{\Gamma} = 1.49\pm0.18$, with an intrinsic scatter of $\sigma^{2}_{\Gamma} = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overline{\Gamma}$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($\Gamma = 1.30$) IMF, as well as the Salpeter ($\Gamma = 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields.

4U 1538-522 is a persistent high mass X-ray binary which exhibits secular spin evolution. In 2019, it underwent a torque reversal from spinning up to spinning down. We performed an extensive study using four NuSATR observations to compare temporal and spectral properties during different states. We observed no abrupt change in luminosity associated with the torque reversal. In addition, the pulse profile, the spectral shape and the power spectrum remained unchanged before and after the torque reversal. The orbital and super-orbital modulation profiles also showed no significant changes. We discuss possible mechanisms for the torque reversal and conclude that it is unlikely to be caused by interactions between the accretion disk and the magnetosphere. Instead, the transition of accretion modes in spherical accretion may be a plausible explanation.

We estimate the absolute age of the globular cluster NGC 3201 using $10,000$ sets of theoretical isochrones constructed through Monte Carlo simulation using the Dartmouth Stellar Evolution Program. These isochrones take into consideration of uncertainty introduced by the choice of stellar evolution parameters. We fit isochrones with 3 detached eclipsing binaries and obtained an age independent of distance. We also fit isochrones with differential reddening corrected HST photometry data utilizing two different Hess diagram based fitting methods. Results from 3 different methods analyzing 2 different types of data agree to within $1 \sigma$, and we find the absolute age of NGC 3201 $= 11.85 \pm 0.74$ Gyr. We also perform variable importance analysis to study the uncertainty contribution from individual parameters and we find the distance is the dominance source of uncertainty in photometry based analysis while total metallicity, Helium abundance, $\alpha$-element abundance, mixing length, and treatment of helium diffusion are important source of uncertainties for all 3 methods.

The distribution of small planet radius ($<$4 R$_\oplus$) is an indicator of the underlying processes governing planet formation and evolution. We investigate the correlation between the radius distribution of exoplanets in Kepler multiplanet systems and the system-level complexity in orbital period spacing. Utilizing a sample of 234 planetary systems with three or more candidate planets orbiting FGK main-sequence stars, we measure the gap complexity ($C$) to characterize the regularity of planetary spacing and compare it with other measures of period spacing and spacing uniformity. We find that systems with higher gap complexity exhibit a distinct radius distribution compared to systems with lower gap complexity. Specifically, we find that the radius valley, which separates super-Earths and sub-Neptunes, is more pronounced in systems with lower gap complexity ($C$$<$0.165). Planets in high complexity systems ($C$$>$0.35) exhibit a lower frequency of sub-Earths (2.5 times less) and sub-Neptunes (1.3 times less) and a higher frequency of super-Earths (1.4 times more) than planets in low complexity systems. This may suggest that planetary systems with more irregular spacings are more likely to undergo dynamic interactions that influence planet scattering, composition, and atmospheric retention. The gap complexity metric proves to be a valuable tool in linking the physical characteristics of planets to their orbital configurations.

In the era of state-of-the-art space-borne telescopes, high-resolution ground-based observation has emerged as a crucial method for characterizing exoplanets, providing essential insights into their atmospheric compositions. In the optical and NIR regions, high-resolution spectroscopy has been powerful for hot Jupiters (HJ) and ultra-hot Jupiters (UHJ) during their primary transits, as it can probe molecules with better sensitivity. Here, we focus on a comparative simulation study of WASP-76 b (UHJ) and WASP-77 A b (HJ) for different number of transits, utilizing three ground-based spectrographs (GIANO-B (TNG), CARMENES (CAHA), and ANDES (E-ELT)) with varying instrumental parameters, spectral coverages, and resolutions. We aim to evaluate the feasibility of the upcoming ground-based European Extremely Large Telescope (E-ELT) in probing molecules from planet atmospheres and how it surpasses other ground-based observatories in terms of detectability. With the 1-D model, petitCODE, we have self-consistently simulated the atmospheric pressure-temperature profiles, which are subsequently integrated into the 1-D chemical kinetics model, VULCAN, to evolve the atmospheric chemistry. High-resolution spectra are obtained by performing line-by-line radiative transfer using petitRADTRANS. Finally, we use the resulting spectra to assess the detectability (sigma_det) of molecular bands, employing the ground-based noise simulator SPECTR. Utilizing cross-correlation spectroscopy, we have successfully demonstrated the robust consistency between our simulation study and real-time observations for both planets. ANDES excels overall in molecular detection due to its enhanced instrumental architecture, reinforcing E-ELT's importance for studying exoplanet atmospheres. Additionally, our theoretical simulations predict the detection of CO, NH3, and H2S on WASP-76 b atmosphere with a sigma_det> 3.

Galaxies infalling into clusters undergo both star-formation quenching and morphological transformation due to environmental effects. We investigate these processes and their timescales using a local sample of 20,191 cluster and 11,674 field galaxies from SDSS. By analysing morphology as a function of distance from the star-formation main sequence, we show that environmental influence is especially pronounced for low-mass galaxies, which emerge from the green valley with early-type morphologies before their star formation is fully suppressed. Using the galaxies' positions in the clusters' Projected Phase Space, we examine the evolution of blue cloud, green valley, and red sequence fractions as a function of time since infall. Interestingly, the green valley fraction remains constant with time since infall, suggesting a balanced flow of galaxies in and out of this class. We estimate that galaxies less massive than $10^{10}\rm M_{\odot}$ spend approximately 0.4 Gyr in the green valley. By comparing quenched and early-type populations, we provide further evidence for the ``slow-then-rapid'' quenching model and suggest that it can also be applied to morphological transitions. Our results indicate that morphological transformation occurs at larger radii than complete star-formation quenching. About 75% of galaxies undergoing morphological transition in clusters are spirals evolving into S0s, suggesting that infalling galaxies retain their disks, while massive ellipticals are relics of early merger events. Finally, we show it takes approximately 2.5 and 1.2 Gyr after the delay-time ($\sim 3.8 {\rm Gyr}$) for the population of low mass galaxies in clusters to reach a 50% threshold in quenched and early-type fraction, respectively. These findings suggest morphological transition precedes full star formation quenching, with both processes possibly being causally linked.

Planet formation models are necessary to understand the origins of diverse planetary systems. Circumstellar disc substructures have been proposed as preferred locations of planet formation but a complete formation scenario has not been covered by a single model so far. We aim to study the formation of giant planets facilitated by disc substructure and starting with sub-micron-sized dust. We connect dust coagulation and drift, planetesimal formation, $N$-body gravity, pebble accretion, planet migration, planetary gas accretion and gap opening in one consistent modelling framework. We find rapid formation of multiple gas giants from the initial disc substructure. The migration trap near the substructure allows the formation of cold gas giants. A new pressure maximum is created at the outer edge of the planetary gap, which triggers the next generation of planet formation resulting in a compact chain of giant planets. A high planet formation efficiency is achieved as the first gas giants are effective in preventing dust from drifting further inwards, which preserves materials for planet formation. Sequential planet formation is a promising framework to explain the formation of chains of gas and ice giants.

Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These ``ultra-long period'' (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission at low Galactic latitude as part of a wide-field survey for transient and variable radio sources. ASKAP\,J175534.9$-$252749.1 does not appear to repeat, with only a single intense two-minute $\sim 200$\,mJy burst detected from 60~hours of observations. The burst morphology and polarisation properties are comparable to those of classical pulsars but the duration is more than one hundred times longer, analogous to ULPs. No comparable bursts are detected in the rest of our widefield survey to date. Combined with the existing ULP population, this suggests that these sources have a strong Galactic latitude dependence and hints at an unexplored population of transient and variable radio sources in the thin disk of the Milky Way. The resemblance of this burst with both ULPs and pulsars calls for a unified coherent emission model for objects with spin periods from milliseconds to tens of minutes. However, whether or not these are all neutron stars or have the same underlying power source remains open for debate.

One of the ultimate goals of the ESA Ariel space mission is to shed light on the formation pathways and evolution of planetary systems in the Solar neighbourhood. Such an endeavour is only possible by performing a large chemical survey of not only the planets, but also their host stars, inasmuch as stellar elemental abundances are the cipher key to decode the planetary compositional signatures. This work aims at providing homogeneous abundances of C, N, and O of a sample of 181 stars belonging to the Tier 1 of the Ariel Mission Candidate Sample. We applied the spectral synthesis and the equivalent width methods to a variety of atomic and molecular indicators (C I lines at 5052 and 5380.3 A, [O I] forbidden line at 6300.3 A, C_2 bands at 5128 and 5165 A, and CN band at 4215 A) using high-resolution and high S/N spectra collected with several spectrographs. We provide carbon abundances for 180 stars, nitrogen abundances for 105 stars, and oxygen abundances for 89 stars. We analyse the results in the light of the Galactic chemical evolution, and in terms of the planetary companions properties. Our sample basically follows the typical trends with metallicity expected for the [C/Fe], [N/Fe], and [O/Fe] abundance ratios. The fraction between C and O abundances, both yields of primary production, is consistent with a constant ratio as [O/H] increases, whereas the abundance of N tends to increase with the increasing of the O abundance, supporting the theoretical assumption of a secondary production of nitrogen. The [C/N], [C/O], and [N/O] ratios are also correlated with [Fe/H], which might introduce biases in the interpretation of the planetary compositions and formation histories if host stars of different metallicity are compared. We provide relations that can be used to qualitatively estimate whether the atmospheric composition of planets is enriched or not with respect to the host stars.

RR Lyrae stars are excellent tracers of the old population II due to their period-luminosity (PL) and period-luminosity-metallicity (PLZ) relations. While these relations have been investigated in detail in many photometric bands, there are few comprehensive studies about them in Sloan-like systems. We present PL and PLZ relations (as well as their counterparts in Wesenheit magnitudes) in the Sloan-Pan-STARSS gP1rP1iP1 bands obtained for Galactic RR Lyrae stars in the vincinity of the Sun. The data used in this paper were collected with the network of 40 cm telescopes of the Las Cumbres Observatory, and geometric parallaxes were adopted from Gaia Data Release 3. We derived PL and PLZ relations separately for RRab and RRc-type stars, as well as for the mixed population of RRab+RRc stars. To our knowledge, these are the first PL and PLZ relations in the Sloan bands determined using RR Lyrae stars in the Galactic field.

Binary mass transfer can occur at high rates due to rapid expansion of the donor's envelope. In the case where mass transfer is unstable, the binary can rapidly shrink its orbit and lead to a merger. In this work we consider the appearance of the system preceding merger, specifically for the case of a low-mass ($\approx 2.5$-$3~M_\odot$) helium star with a neutron star (NS) companion. Modeling the mass transfer history as well as the wind launched by super-Eddington accretion onto the NS, we find that such systems can power slowly rising transients with timescales as long as years, and luminosities of $\sim 10^{40}$-$10^{41}$ erg s$^{-1}$ from optical to UV. The final explosion following the merger (or core-collapse of the helium star in some cases) leads to an interaction-powered transient with properties resembling Type Ibn supernovae (SNe), possibly with a bright early peak powered by shock cooling emission for merger-powered explosions. We apply our model to the Type Ibn SN 2023fyq, that displayed a long-term precursor activity from years before the terminal explosion.

Based on the $Insight$-HXMT archival data, we have detected a new atypical low-frequency quasi-periodic oscillation (LFQPO) in the black hole X-ray binary MAXI J1348$-$630. The new LFQPO is detected in all the three instruments of $Insight$-HXMT with a combined significance of 3--5 $\sigma$, covering a wide energy range of 1--100 keV. The fractional root-mean-square (RMS) seems decrease with energy. It exclusively appears in the hard state during both the main and mini outburst, spanning an X-ray intensity range by a factor of 10, and a very narrow hardness range. The frequency of this new type of LFQPO is moderately stable, in the range of 0.08--0.15 Hz. We discussed different models for the LFQPO, and found none is able to explain the observed properties of this new type of LFQPO.

This work is part of VESTALE, a project initiated within the Rubin-LSST Cadence Strategy Optimization Process . Its goal is to explore the potential of Rubin-LSST observations aimed at the Galaxy's bulge (Bulge) for studying RR Lyrae stars (RRL). Observation and analysis of RR Lyrae stars in the Bulge are crucial for tracing the old population of the central part of our galaxy and reconstructing the history of Bulge formation. Based on observations conducted with CTIO/DECam by Saha et al. 2019 towards the Baade Window, our simulations demonstrate that early Rubin-LSST observations will enable the recovery of RR Lyrae light curves at Galactic center distances with sufficient precision. This will allow us to utilize theoretical relations from Marconi et al. 2022 to determine their distances and/or metallicity, following the REDIME algorithm introduced in Bono et al. 2019. We show how reddening and crowding affect our simulations and highlight the importance of considering these effects when deriving pulsation parameters (luminosity amplitudes, mean magnitudes) based on the light curves especially if the goal is to explore the opposite side of the Bulge through the observation of its RRL. The simulations discussed in this investigation were conducted to support the SCOC's decision to observe this important sky region since it has only recently been decided to include part of the Bulge as a target within the LSST main survey.

Primordial black holes (PBHs) of stellar mass or heavier would accrete baryonic gas, which becomes denser and hotter, injecting energetic photons in the cosmological medium soon after cosmic recombination, in the so-called dark ages. The ionisation history of the universe would be altered, an effect cosmic microwave background (CMB) temperature and polarisation anisotropies is sensitive to. The magnitude of the effect depends on the abundance and mass distribution of the PBHs, as well as on still uncertain aspects of accretion luminosity. We review the current understanding of this field, with an emphasis on the peculiarities of the phenomenon in the cosmological context, and present the existing constraints on the PBH abundances.

Context. The current algorithms used for the calibration and analysis of very long baseline interferometry (VLBI) networks that only use linear polarizers (as is the case of the VLBI Global Observing System, VGOS) do not properly account for instrumental and source-intrinsic polarimetry, which can cause errors in geodetic and astronomical products. Aims. We aim to develop a calibration pipeline for VLBI interferometers that observe in a basis of linear polarization, as is the case of VGOS. The products from this pipeline can be used to obtain valuable full-polarization astronomical information from the observed sources, and they can be used to potentially improve the geodetic results. Methods. We used the algorithm PolConvert to write the correlation products in a basis of circular polarization that is compatible with the standard VLBI calibration procedures. In addition to this, we implemented a wide-band global fringe-fitting algorithm that accounts for dispersive effects (ionospheric delay) and allows us to perform full-polarization imaging of all the observed sources, covering the whole frequency band of VGOS. Results. We present the outcome of our pipeline applied to a global IVS VGOS epoch of observations and show example imaging results in total intensity and polarization. We also discuss issues encountered during the analysis and suggest points of improvement in the VGOS system for an optimum geodetic and astronomical exploitation of this interferometer.

Magnetic fields connect an array of planetary processes, from atmospheric escape to interior convection. Despite their importance, exoplanet magnetic fields are largely unconstrained by both theory and observation. In this Letter, we propose a novel method for constraining the B field strength of hot gas giants: comparing the velocities of heavy ions and neutral gas with high-resolution spectroscopy. The core concept of this method is that ions are directly deflected by magnetic fields. While neutrals are also affected by B fields via friction with field-accelerated ions, ionic gas should be more strongly coupled to the underlying magnetic field than bulk neutral flow. Hence, measuring the difference between the two velocities yields rough constraints on the B field, provided an estimate of the stellar UV flux is known. We demonstrate that heavy ions are particularly well suited for this technique, because they are less likely to be entrained in complex hydrodynamic outflows than their lighter counterparts. We perform a proof-of-concept calculation with Ba II, an ion whose velocity has been repeatedly measured at high confidence with high-resolution spectroscopy. Our work shows that a 10G magnetic field would produce ~ km/s ion--neutral velocity differences at a microbar, whereas a 50G magnetic field would produce ~20km/s velocity difference. With new leverage on magnetic fields, we will be able to investigate magnetic field generation in the extreme edge cases of hot gas giants, with wide-ranging consequences for planetary interior structure, dynamo theory, and habitability.

Transformer models have recently become very successful in the natural language domain. Their value as sequence-to-sequence translators there, also makes them a highly interesting technique for learning relationships between astrophysical time series. Our aim is investigate how well such a transformer neural network can establish causal temporal relations between different channels of a single-source signal. We thus apply a transformer model to the two phases of Gamma-Ray Bursts (GRBs), reconstructing one phase from the other. GRBs are unique instances where a single process and event produces two distinct time variable phenomena: the prompt emission and the afterglow. We here investigate if a transformer model can predict the afterglow flux from the prompt emission. If successful, such a predictive scheme might then be distilled to the most important underlying physics drivers in the future. We combine the transformer model with a novel dense neural network setup to directly estimate the starting value of the prediction. We find that the transformer model can, in some instances, successfully predict different phases of canonical afterglows, including the plateau phase. Hence it is a useful and promising new astrophysical analysis technique. For the GRB test case, the method marginally exceeds the baseline model overall, but still achieves accurate recovery of the prompt-afterglow fluence-fluence correlation in reconstructed light curves. Despite this progress, we conclude that consistent improvement over the baseline model is not yet achieved for the GRB case. We discuss the future improvements in data and modeling that are required to identify new physical-relation parameters or new insights into the single process driving both GRB phases.

Excess chromospheric emissions within deep photospheric lines are effective proxies of stellar magnetism for FGK stars. This emission decays with stellar age and is a potential determinant of this important stellar quantity. We report absolutely calibrated H$\alpha$ chromospheric fluxes for 511 solar-type stars in a wide interval of precisely determined masses, $[$Fe/H$]$, ages, and evolution states from high S/N, moderately high$-$resolution spectra. The comparison of H$\alpha$ and H+K chromospheric fluxes reveals a metallicity bias (absent from H$\alpha$) affecting Ca II H+K fluxes thereby metal-rich stars with deep line profiles mimic low chromospheric flux levels, and vice versa for metal-poor stars. This bias blurs the age-activity relation, precluding age determinations for old, inactive stars unless mass and $[$Fe/H$]$ are calibrated into the relation. The H+K lines being the most widely studied tool to quantify magnetic activity in FGK stars, care should be exercised in its use whenever wide ranges of mass and $[$Fe/H$]$ are involved. The H$\alpha$ age-activity-mass-metallicity calibration appears to be in line with the theoretical expectation that (other parameters being equal) more massive stars possess narrower convective zones and are less active than less massive stars, while more metal-rich stars have deeper convective zones and appear more active than metal-poorer stars. If regarded statistically in tandem with other age diagnostics, H$\alpha$ chromospheric fluxes may be suitable to constrain ages for FGK stars with acceptable precision.

We carry out a parameter study of interacting disk galaxies with impact parameters ranging from central collisions to weakly interacting scenarios. The orientations of the disks are also varied. In particular, we investigate how magnetic field amplification depends on the these parameters. We use magnetohydrodynamics for gas disks in combination with live dark matter halos in adaptive mesh refinement simulations. The disks are initialized using a setup for isolated disks in hydrostatic equilibrium. Small-scale filtering of the velocity and magnetic field allows us to estimate the turbulent electromotive force (EMF) and kinetic helicity. Time series of the average magnetic field in central and outer disk regions show pronounced peaks during close encounters and mergers. This agrees with observed magnetic fields at different interaction stages. The central field strength exceeds 10 microgauss (corresponding to an amplification factor of 2 to 3) for small impact parameters. As the disks are increasingly disrupted and turbulence is produced by tidal forces, the small-scale EMF reaches a significant fraction of the total EMF. The small-scale kinetic helicity is initially antisymmetric across the disk plane. Albeit its evolution is sensitive to both the impact parameter and inclinations of the rotation axes with respect to the relative motion of the disks, antisymmetry is generally broken through interactions and the merger remnant has lost most of the initial helicity. The EMF and the magnetic field also decay rapidly after coalescence. The strong amplification during close encounters of the interacting galaxies are mostly driven by helical flows and a mean-field dynamo. The small-scale dynamo contributes significantly in post-interaction phases. However, the amplification of the magnetic field cannot be sustained.

We explored the impact that Doppler dimming and brightening effects from bulk motions of solar prominences have on the formation of Lya, Ha, and MgII h line profiles. We compared two schemes in which these effects manifest; when the prominence is moving radially away from the solar surface (radial case), and when the prominence is moving parallel to the solar surface (horizontal case). To do this, we analysed 13,332 model profiles generated through the use of the 1D NLTE (i.e. departures from Local Thermodynamic equilibrium) radiative transfer (RT) code Promweaver, built on the Lightweaver NLTE RT framework to mimic the behaviour and output of the 1D NLTE RT code PROM. We found that horizontal velocities are just as, or more important than radial velocities. This demonstrates that horizontal velocities need to be accounted for when attempting to do any sort of forward modelling.

We accurately identify and classify the variability of A-F stars in the southern continuous viewing zone of the TESS satellite. The brightness limit was set to 10 mag to ensure the utmost reliability of our results and allow for spectroscopic follow-up observations using small telescopes. We aim to compare our findings with existing catalogues of variable stars. The light curves from TESS and their Fourier transform were used to manually classify stars in our sample. Cross-matching with other catalogues was performed to identify contaminants and false positives. We have identified 1171 variable stars (51 % of the sample). Among these variable stars, 67 % have clear classifications, which includes $\delta$ Sct and $\gamma$ Dor pulsating stars and their hybrids, rotationally variables, and eclipsing binaries. We have provided examples of the typical representatives of variable stars and discussed the ambiguous cases. We found 20 pairs of stars with the same frequencies and identified the correct source of the variations. Additionally, we found that the variations in 12 other stars are caused by the contamination with the light of faint nearby large-amplitude variable stars. To compare our sample with other variable star catalogues, we have defined two parameters reflecting the agreement in identification of variable stars and their classification. This comparison reveals intriguing disagreements in classification ranging from 52 % to 100 %. However, if we assume that stars without specific types are only marked as variable, then the agreement is relatively good, ranging from 57 % to 85 % (disagreement 15-43 %). We have demonstrated that the TESS classification is superior to the classification based on other photometric surveys. The classification of stellar variability is complex and requires careful consideration. Caution should be exercised when using catalogue classifications.

Observations of gravitationally lensed, high-mass stars at redshifts $\gtrsim1$ occasionally reveal spectral energy distributions that contain two components with different effective temperatures. Given that two separate stars are involved, it suggests that both stars have simultaneously reached very high magnification, as expected for two stars in a binary system close to the caustic curve of the foreground galaxy-cluster lens. The inferred effective temperatures and luminosities of these stars are, however, difficult to reconcile with known binaries, or even with isolated stars of the same age. Here, we explore three alternative explanations for these cases: circumstellar dust around the cooler of the two stars; age differences of a few Myr among stars in the same star cluster, and a scenario in which the stars originate in two separate star clusters of different age along the lensing caustic. While all of these scenarios are deemed plausible in principle, dust solutions would require more circumstellar extinction than seen in local observations of the relevant super/hypergiant stars. Hence, we argue that age differences between the two stars are the most likely scenario, given the current data.

Pulsars are highly magnetized rotating neutron stars, emitting in a broad electromagnetic energy range. Reproducing the observed population of pulsars refines our understanding of their formation, evolution as well as radiation processes and geometry. In this paper, we improve our previous pulsar population synthesis by studying the impact of the Galactic gravitational potential and of the death line putting emphasis on the gamma-ray pulsar population. In order to elucidate the necessity of a death line, refined initial distributions for spin periods and radial positions at birth were implemented, elevating the sophistication of our previous simulations to the most recent state-of-the-art. Our pulsar population synthesis takes into account the secular evolution of a force-free magnetosphere and the magnetic field decay. The radio and gamma-ray emissions are modeled respectively by the polar cap geometry and the striped wind model. When simulating a cohort of ten million pulsars akin to previous endeavors, the integration of a death line proves essential for a better agreement with observational trends. Yet, discrepancies emerges when modeling a reduced ensemble of only one million pulsars, endowed with a higher birth rate. This divergence yields a remarkably refined $P-\dot{P}$ diagram, conspicuously with or without any delineating death line. Intriguingly, this observation suggests an overestimation of the ages of detected pulsars in existent data sets, thereby challenging the necessity of death lines in pulsar population synthesis. Validation through Kolmogorov-Smirnov testing underlines the statistical congruence between the spin period and spin period derivative distributions of both observations and simulations. Moreover simulations with a higher gamma-ray instrumental sensitivity show that the GeV excess in the Galactic centre could be of pulsar origin.

We examine the classical theory of stellar wind formation. The theory requires that to form a supersonic stellar wind, a subsonic flow speed must start at a specific initial speed from the coronal base, called eigenspeed, go along a continuous eigenfunction, and reach the sonic point, which is where the flow speed equals the sonic speed, while the critical condition, which is where the effective driving force is zero, is satisfied. Any mismatch between the sonic point and critical condition distances results in either subsonic winds when the initial speed is below the eigenspeed, or no stellar wind when the initial speed is above the eigenspeed. However, because the critical condition is determined by the stellar wind temperature profile, which depends on ionization process at the top of the chromosphere and the heating process around the coronal base but not by the processes at the sonic point, the requirement of sonic point satisfying the critical condition is generally not met and hence the momentum equation in the conventional theory encounters difficulty when the initial speed is above the eigenspeed. To resolve the difficulty, we propose a discontinuity between the sonic point and the critical condition to reach supersonic stellar wind solutions. As a result, supersonic stellar winds can be produced when the initial speed in the coronal base is greater than the eigenspeed. The critical solution or eigen function provided by the conventional stellar wind model describes the condition that separates the supersonic stellar winds from subsonic ones.

Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast HI gas collisions (~100 km/s) as suggested by recent observations and their subsequent evolution into YMCs by using three-dimensional magnetohydrodynamics simulations involving self-gravity and detailed thermal/chemical processes. In particular, the impact of ionization feedback from stellar radiation is included in an approximate fashion where the temperature within the HII regions is elevated to 10,000 K, while supernova feedback is not included. We examine whether the resulting massive clumps can survive this ionization feedback and evolve into YMCs. Our simulations reveal the emergence of gas clumps that do not only possess substantial mass (~10^5 M_sun) but also sufficient compactness (~5 pc). Notably, these clumps exhibit significantly higher escape velocities compared to the sound speed of the HII region, indicating effective gravitational retention of gas against feedback-induced evaporation. Consequently, these conditions foster efficient star formation within the massive gas clumps, ultimately leading to their evolution into YMCs. We also perform simulations involving lower-velocity gas collisions, approximately 15 km/s, typical shock velocities induced by galactic superbubbles.

We present the first spectroscopic estimates of the chemical abundance of M dwarf stars in a globular cluster (GC), namely 47 Tucanae. By exploiting NIRSpec on board the James Webb Space Telescope (JWST) we gathered low-resolution spectra for 28 stars with masses in the range ~0.4-0.5 solar masses. The spectra are strongly affected by the H2O water vapour bands which can be used as indicators of the oxygen abundance. The spectral analysis reveals that the target stars feature a different O abundance, with a difference of ~0.40 dex between first and the most-polluted second population. The observed range is similar to that observed among red giant stars. This result reinforces previous findings based on the analysis of photometric diagrams, including the ``chromosome maps'', providing a first, and more direct, evidence of light element variations in the M dwarfs' mass regime. The observation that the multiple populations, with their variations in light elements, exhibit the same patterns from the lower main sequence all the way to the red giant branch further strengthens the notion that multiple stellar populations in globular clusters formed in a series of bursts of star formation.

Octans is one of the most distant ($d\sim150$pc) young stellar associations of the solar neighbourhood. Its age is still poorly constrained in the literature and requires further investigation. We take advantage of the state-of-the-art astrometry delivered by the third data release of the Gaia space mission combined with radial velocity measurements obtained from high-resolution spectroscopy to compute the 3D positions and 3D spatial velocities of the stars and derive the dynamical traceback age of the association. We performed an extensive traceback analysis using different subsamples of stars, different metrics to define the size of the association, and different models for the Galactic potential to integrate the stellar orbits in the past. We derive a dynamical age of $34^{+2}_{-2}$Myr that is independent from stellar models and represents the most precise age estimate currently available for the Octans association. After correcting the radial velocity of the stars for the effect of gravitational redshift, we obtain a dynamical age of $33^{+3}_{-1}$Myr, which is in very good agreement with our first solution. This shows that the effect of gravitational redshift is small for such a distant young stellar association. Our result is also consistent with the less accurate age estimates obtained in previous studies from lithium depletion (30-40Myr) and isochrones (20-30Myr). By integrating the stellar orbits in time, we show that the members of Octans and Octans-Near had different locations in the past, which indicates that the two associations are unrelated despite the close proximity in the sky. Our results confirm that it is possible to derive precise dynamical ages via the traceback method for $\sim30$Myr old stellar clusters at about $\sim150$pc with the same precision level that has been achieved in other studies for young stellar groups within 50pc of the Sun.

Despite the large number of studies focused on the characterization of Li-rich stars and understanding of the mechanisms leading to such enrichment, their origin remains a mystery. Magnetic activity, in particular the phenomena usually associated with it, such as spots and plages, and the Li abundance (A(Li)) of stars, are in general thought to be connected. However, as of today it is unclear just how. In this work, we study a sample of young but evolved intermediate mass red giants, inhabitants of open clusters where planets have been searched. We aim at using radial velocity (RV) and stellar activity indicator signals to look for relations between Li abundances and stellar activity/variability. We explore how the standard deviation (STD), peak to peak amplitude (PTP), mean and median of typical stellar activity indicators (BIS, FWHM, $T_\textrm{eff}$, and H$\alpha$ index) change as a function of the Li content of 82 red giants. Furthermore, we compute weighted Pearson Correlation Coefficients ($\rho_w$) between time series of RV measurements and the stellar activity indicators for the stars in our sample. To aid our results, we also study Generalized Lomb-Scargle Periodograms (GLSP) to capture possible significant periodic temporal variations in our data. Our analysis indicates that the STD and PTP of BIS and FWHM, the mean and median of the H$\alpha$ index and $v\sin(i)$ increase exponentially with A(Li) in our sample of red giants. Significant temporal variations and correlations between RVs and activity indicators also tend to be found preferentially for stars where high A(Li) is observed. Most of the Li-rich stars in our sample either show strong correlations of RV with at least one of the stellar activity indicators or reveal significant periodic temporal variations in their stellar activity indicators GLSPs consistent with those found for RV.

The need of real-time of monitoring and alerting systems for Space Weather hazards has grown significantly in the last two decades. One of the most important challenge for space mission operations and planning is the prediction of solar proton events (SPEs). In this context, artificial intelligence and machine learning techniques have opened a new frontier, providing a new paradigm for statistical forecasting algorithms. The great majority of these models aim to predict the occurrence of a SPE, i.e., they are based on the classification approach. In this work we present a simple and efficient machine learning regression algorithm which is able to forecast the energetic proton flux up to 1 hour ahead by exploiting features derived from the electron flux only. This approach could be helpful to improve monitoring systems of the radiation risk in both deep space and near-Earth environments. The model is very relevant for mission operations and planning, especially when flare characteristics and source location are not available in real time, as at Mars distance.

Atmospheric spectroscopy provides a window into the properties of exoplanets. However, the physical interpretation of retrieved data and its implications for the internal properties of exoplanets remains nebulous. This letter addresses three misconceptions held by some atmospheric spectroscopists regarding the connection between observed chemical abundances and theory: (1) Whether atmospheric spectroscopy can provide the bulk atmospheric chemistry, (2) whether it can identify if a planet is cloudless, and (3) whether atmospheric evaporation arguments can be used to dismiss certain compositions inferred through spectroscopy. This letter concludes by exploring applications of remote sensing in the quest for the search for life outside of our solar system.

The nearby oxygen-rich AGB star L2 Pup hosts a well-studied nearly edge-on disk. To date, disks around AGB stars have not been chemically studied in detail. By combining a parameterisation commonly used for protoplanetary disks and archival ALMA observations, we retrieved an updated density and temperature structure of this disk. This physical model was then used as input to the first chemical model of an AGB disk. The model shows that the physical structure of the disk has a large impact on its chemistry, with certain species showing large changes in column density relative to a radial outflow, indicating that chemistry could be used as a tracer of disks that cannot be directly imaged. Despite its oxygen-rich nature, the daughter species formed within the disk are surprisingly carbon-rich. Two chemical regimes can be distinguished: cosmic-ray induced chemistry in the midplane and photochemistry induced by the interstellar radiation field in the outer regions. Certain complex organic molecules are formed in the midplane. This occurs via gas-phase chemistry only, as the disk is too warm for dust-gas chemistry. The photochemistry in the outer regions leads to the efficient formation of (long) carbon-chains. The predictions of the model allow us to tentatively put the disk's age $\lesssim 10^5$ yr. Additional observations are necessary to better constrain the physical structure of L2 Pup's disk and are essential to test the predictions made by the chemical model. Our exploratory work paves the way for a more general study of the chemistry of AGB disks.

We identify the bright Am-type star HD 181793 to be a previously-unknown eclipsing, chemically peculiar heartbeat binary, the second of its kind known. The system carries an orbital period of $P = 11.47578275 \pm 0.00000055$ days. We use TESS photometry and LCOGT NRES radial velocity data to build a self-consistent orbital model and determine the fundamental stellar characteristics of the primary. We use a spectral separation method to unveil the secondary and measure the masses of both stars. The radial velocity amplitude of the primary, $K_1 = 47.41+0.13-0.12 km s^{-1}$, gives a mass $M_1 = 1.57 \pm 0.01 $ Msun. The secondary radial velocity amplitude $K_2 = 84.95+0.12-0.09 km s^{-1}$ yields a mass ratio $q = 0.558 \pm 0.002$ and a secondary mass $M_2 = 0.87 \pm 0.01 $ Msun. From the spectral energy distribution and Gaia parallax we find a radius $R_1 = 2.04 \pm 0.05$ Rsun. The grazing transit profile and spectroscopic luminosity ratio indicate $R_2 = 1.04+0.15-0.10$ Rsun, suggesting an early-K spectral type. We show that the heartbeat feature in the TESS light curve can be explained by time-varying ellipsoidal variation, driven by the orbital eccentricity of $e = 0.3056+0.0024-0.0026$, and relativistic beaming of the light of the primary. We find no evidence of tidally-excited oscillations.

The search for habitable conditions and signs of life on exoplanets is a major frontier in modern astronomy. Detecting atmospheric signatures of Earth-like exoplanets is challenging due to their small sizes and relatively thin atmospheres. Recently, a new class of habitable sub-Neptune exoplanets, called Hycean worlds, has been theorized. Hycean worlds are planets with H2-rich atmospheres and planet-wide oceans with thermodynamic conditions similar to those in the Earth's oceans. Their large sizes and extended atmospheres, compared to rocky planets of similar mass, make Hycean worlds significantly more accessible to atmospheric observations. These planets open a new avenue in the search for planetary habitability and life elsewhere using spectroscopic observations with the James Webb Space Telescope (JWST). We observed the transmission spectrum of a candidate Hycean world, K2-18 b, recently with JWST in its first year of operations. The spectrum reveals multiple spectral features of carbon bearing molecules in the planetary atmosphere, leading to the first detections of methane (CH4) and carbon dioxide (CO2) in a habitable-zone exoplanet. We discuss inferences of the atmospheric chemical composition and its implications for the atmospheric, interior and surface conditions on the planet, along with the possibility of a habitable ocean underneath the atmosphere. We discuss new observational and theoretical developments in this emerging frontier and their implications for exoplanetary habitability and search for life elsewhere.

Repeating fast radio bursts can exhibit a wide range of burst repetition rates, from none to hundreds of bursts per hour. Here, we report the detection and characteristics of 57 bursts from the recently discovered FRB 20240114A, observed with GMRT in the frequency ranges 300-500 MHz and 550-750 MHz. Majority of the bursts show narrow emission-bandwidth with $\Delta\nu/\nu \sim$ around 10 %. All of the bursts we detect are faint ($<$10 Jy ms), and thus probe the lower end of the energy distribution. We determine the rate function for FRB 20240114A at 400 MHz, and downward drift rates at 400 and 650 MHz, and discuss our measurements in the context of the repeating FRB population. We observe sudden variations in the burst activity of FRB 20240114A over time. Our data as well as the other publicly available information on other observations of FRB 20240114A so far, there is an indication that FRB 20240114A potentially exhibit a chromaticity in its burst activity. While the burst properties of FRB 20240114A are similar to ther repeating FRBs, the frequency-dependent activity, if established, could provide crucial clues to the origin of repeating FRBs. We also place the most stringent 5$\sigma$ upper limits of 600 $\mu$Jy and 89 $\mu$Jy on any persistent radio source (PRS) associated with FRB 20240114A at 400 MHz and 650 MHz, respectively, and compare these with the luminosity of the known PRSs associated with FRB121102A and FRB190520B.

We investigate the planetary migration of low-mass planets ($M_p\in[1,15]M_\oplus$, here $M_\oplus$ is the Earth mass) in a gaseous disc containing a previously formed gap. We perform high-resolution 3D simulations with the FARGO3D code. To create the gap in the surface density of the disc, we use a radial viscosity profile with a bump, which is maintained during the entire simulation time. We find that when the gap is sufficiently deep, the spiral waves excited by the planet trigger the Rossby wave instability, forming cyclonic (underdense) vortices at the edges of the gap. When the planet approaches the gap, it interacts with the vortices, which produce a complex flow structure around the planet. Remarkably, we find a widening of the horseshoe region of the planet produced by the vortex at the outer edge of the gap, which depending on the mass of the planet differs by at least a factor of two with respect to the standard horseshoe width. This inevitably leads to an increase in the corotation torque on the planet and produces an efficient trap to halt its inward migration. In some cases, the planet becomes locked in corotation with the outer vortex. Under this scenario, our results could explain why low-mass planets do not fall towards the central star within the lifetime of the protoplanetary disc. Lastly, the development of these vortices produces an asymmetric temporal evolution of the gap, which could explain the structures observed in some protoplanetary discs.