Papers for Thursday, Oct 17 2024

This paper provides a comprehensive guide to gravitational wave data processing, with a particular focus on signal generation, noise modeling, and optimization techniques. Beginning with an introduction to gravitational waves and the detection techniques used by LIGO and Virgo, the manual covers the essentials of signal processing, including Fourier analysis, filtering, and the generation of quadratic chirp signals. The analysis of colored Gaussian noise and its impact on interferometric detectors like LIGO is explored in detail, alongside signal detection methods such as the Generalized Likelihood Ratio Test (GLRT). The paper also delves into optimization techniques like Particle Swarm Optimization (PSO), which can be applied to improve signal estimation accuracy. By providing MATLAB-based implementations, this manual serves as both a theoretical and practical resource for researchers in the field of gravitational wave astronomy.

Observations reveal that magnetic fields on neutron stars (NSs) are in the range of $10^{8-15}$ G. Apart from being celestial bodies, NSs are normally rotating. In this work, we study the impact of a chaotic magnetic field on the mass-radius relation of the rotating NSs. We employ an equation of state of NSs with the nuclei in the crust and hyperons in the core. We use Hartle-Thorne formalism as an approximation of the rotating NSs. For the magnetic field ansatz, we use the one coupled to the energy density. We find that the magnetic field decreases the mass of NS. In contrast, the increment of the magnetic field increases the deformation of rotating NSs.

The hydrogen recombination lines H30$\alpha$, H40$\alpha$, H42$\alpha$, H50$\beta$ and H57$\gamma$ and the underlying bremsstrahlung continuum emission were detected with ALMA in the bipolar nebula Mz3. The source was not spatially resolved, but the velocity profile of the H30$\alpha$ line shows clear indication of maser amplification, confirming previous reports of laser amplification in the far infrared H recombination lines observed with Herschel Space Observatory. Comparison between the flux densities of the H50$\beta$, H40$\alpha$ and H42$\alpha$ lines show overcooling, or darkness amplification by stimulated absorption (dasar effect) at the LSR velocity of about $-25$ km s$^{-1}$, which constrains the density of the absorbing region to about 10$^3$ cm$^{-3}$. The H30$\alpha$ line, on the other hand, presents maser lines at LSR velocities of $-69$ and $-98$ km s$^{-1}$, which indicates ionized gas with densities close to 10$^7$ cm$^{-3}$. Although the source of emission was not resolved, it was possible to find the central position of the images for each velocity interval, which resulted in a well defined position-velocity distribution.

Chemical evidence indicates that an appreciable fraction of Sun-like stars have engulfed rocky planets during their main-sequence lifetimes. We investigate whether the tidal evolution and destruction of ultra-short-period planets (USPs) can explain this phenomenon. We develop a simple parameterized model for the formation and engulfment of USPs in a population of MS stars. With this model, it is possible to reproduce both the observed occurrence rate of USPs and the frequency of planet-engulfing Sun-like stars for a reasonable range of USP formation rates and tidal decay lifetimes. Our results support a theory of USP formation through gradual inward migration over many Gyr and suggest that engulfment occurs $\sim 0.1$-$1 \, {\rm Gyr}$ after formation. This lifetime is set by tidal dissipation in the USP itself instead of the host star, due to the perturbing influence of external companions. If USP engulfment is the main source of pollution among Sun-like stars, we predict a correlation between pollution and compact multi-planet systems; some $5$-$10\%$ of polluted stars should have a transiting planet of mass $\gtrsim 5 M_{\oplus}$ and period $\sim 4$-$12$ days.

We analyze near-infrared integral field spectropolarimetry of the AB Aurigae protoplanetary disk and protoplanet (AB Aur b), obtained with SCExAO/CHARIS in 22 wavelength channels covering the J, H, and K passbands ($\lambda_{\rm o}$ = 1.1--2.4 $\mu m$) over angular separations of $\rho$ $\approx$ 0.13" to 1.1" ($\sim$20--175 au). Our images resolve spiral structures in the disk in each CHARIS channel. At the longest wavelengths, the data may reveal an extension of the western spiral seen in previous polarimetric data at $\rho$ $<$ 0.3" out to larger distances clockwise from the protoplanet AB Aur b, coincident with the ALMA-detected $CO$ gas spiral. While AB Aur b is detectable in complementary total intensity data, it is a non-detection in polarized light at $\lambda$ $>$ 1.3 $\mu $m. While the observed disk color is extremely red across $JHK$, the disk has a blue intrinsic scattering color consistent with small dust grains. The disk's polarization spectrum is redder than AB Aur b's total intensity spectrum. The polarization fraction peaks at $\sim$ 0.6 along the major disk axis. Radiative transfer modeling of the CHARIS data shows that small, porous dust grains with a porosity of $p$ = 0.6--0.8 better reproduce the scattered-light appearance of the disk than more compact spheres ($p$ = 0.3), especially the polarization fraction. This work demonstrates the utility of integral field spectropolarimetry to characterize structures in protoplanetary disks and elucidate the properties of the disks' dust.

We study the formation of stars with varying amounts of heavy elements synthesized by the rapid neutron-capture process ($r$-process) based on our detailed cosmological zoom-in simulation of a Milky Way-like galaxy with an $N$-body/smoothed particle hydrodynamics code, ASURA. Most stars with no overabundance in $r$-process elements, as well as the strongly $r$-process enhanced $r$-II stars ([Eu/Fe] $>+0.7$), are formed in dwarf galaxies accreted by the Milky Way within the 6 Gyr after the Big Bang. In contrast, over half of the moderately enhanced $r$-I stars ($+0.3 <$ [Eu/Fe] $\leq +0.7$) are formed in the main in-situ disk after 6 Gyr. Our results suggest that the fraction of $r$-I and $r$-II stars formed in disrupted dwarf galaxies is larger the higher their [Eu/Fe] is. Accordingly, the most strongly enhanced $r$-III stars ([Eu/Fe] $> +2.0$) are formed in accreted components. These results suggest that non-$r$-process-enhanced stars and $r$-II stars are mainly formed in low-mass dwarf galaxies that hosted either none or a single neutron star merger, while the $r$-I stars tend to form in the well-mixed in-situ disk. We compare our findings with high-resolution spectroscopic observations of $r$-process-enhanced metal-poor stars in the halo and dwarf galaxies, including those collected by the R-Process Alliance. We conclude that observed [Eu/Fe] and [Eu/Mg] ratios can be employed in chemical tagging of the Milky Way's accretion history.

The announcement in the summer of 2023 about the discovery of evidence for a gravitational wave background (GWB) using pulsar timing arrays (PTAs) has ignited both the PTA and the larger scientific community's interest in the experiment and the scientific implications of its findings. As a result, numerous scientific works have been published analyzing and further developing various aspects of the experiment, from performing tests of gravity to improving the efficiency of the current data analysis techniques. In this regard, we contribute to the recent advancements in the field of PTAs by presenting the most general, agnostic, per-frequency Bayesian search for a low-frequency (red) noise process in these data. Our new method involves the use of a conjugate Jeffrey's-like multivariate prior which allows one to model all unique parameters of the global PTA-level red noise covariance matrix as a separate model parameter for which a marginalized posterior-probability distribution can be found using Gibbs sampling. Even though perfecting the implementation of the Gibbs sampling and mitigating the numerical stability challenges require further development, we show the power of this new method by analyzing realistic and theoretical PTA simulated data sets. We show how our technique is consistent with the more restricted standard techniques in recovering both the auto and the cross-spectrum of pulsars' low-frequency (red) noise. Furthermore, we highlight ways to approximately characterize a GWB (both its auto- and cross-spectrum) using Fourier coefficient estimates from single-pulsar and so-called CURN (common uncorrelated red noise) analyses via analytic draws from a specific Inverse-Wishart distribution.

We characterise the entropy profiles of 32 very high mass ($M_{500}>7.75\times10^{14}~M_{\odot}$) galaxy clusters (HIGHMz), selected from the CHEX-MATE sample, to study the intracluster medium (ICM) entropy distribution in a regime where non-gravitational effects are minimised. Using XMM-Newton measurements, we measure the entropy profiles up to ~$R_{500}$ for all objects. The scaled profiles exhibit large dispersion in the central regions, but converge rapidly to the expectation from pure gravitational collapse beyond the core. We quantify the correlation between the ICM morphological parameters and scaled entropy as a function of radius, showing that morphologically relaxed (disturbed) objects have low (high) central entropy. We compare our data to other observational samples, finding differences in normalisation which are linked to the average mass of the samples in question. We find that a weaker mass dependence than self-similar in the scaling (Am ~ -0.25) allows us to minimise the dispersion in the radial range [0.3-0.8]$R_{500}$ for clusters spanning over a decade in mass. The deviation from self-similarity is radially dependent and is more pronounced at small and intermediate radii than at $R_{500}$. We also investigate the distribution of central entropy $K_0$, finding no evidence for bimodality, and outer slopes $\alpha$, which peaks at ~1.1. Using weak lensing masses, we find indication for a small suppression of the scatter (~30%) beyond the core when using masses derived from Yx in the rescaling. Finally, we compare to recent cosmological numerical simulations from THE THREE HUNDRED and MACSIS, finding good agreement with our observational data. These results provide a robust observational benchmark in the gravity-dominated regime and will serve as a future reference for samples at lower mass, higher redshifts, and for ongoing work using cosmological numerical simulations.

Mounting evidence suggests that the semi-classical description of a black hole breaks down at the latest after losing an O(1) fraction of its mass. As a result, effects such as memory burden can slow down evaporation so that small primordial black holes (PBHs), in particular those in the mass range 10^6 g to 10^9 g, become viable dark matter candidates. In this paper, we investigate the production of PBHs from a prototype model of polynomial inflation with a non-minimal coupling to gravity. We show that a sufficiently small PBH mass alleviates any tension with CMB observations. Moreover, we develop efficient numerical procedures to identify model parameters and evolve Mukhanov-Sasaki modes to place bounds on the scalar-induced stochastic gravitational wave (GW) background. Whilst we identify some prospects for observation with future GW detectors, our results highlight the need to develop new experiments for high-frequency GW detection in the ~kHz to ~MHz range. Finally, we demonstrate that previously-used ansätze for modelling the power spectrum only yield a reliable approximation for the GW signal if some input from inflation is used.

The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. We apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.

Galactic globular clusters contain two main groups of stars, the pristine or 1P stars, and the polluted or 2P stars. The pristine-star fraction in clusters, $F_{1P}$, is a decreasing function of the cluster present-day mass, $m_{prst}$. Paper~I has introduced a model mapping the region of the $(m_{prst},F_{1P})$ space occupied by clusters, with the cluster mass threshold for 2P-star formation a key building-block. We now expand this model to the pristine-star fraction in dependence of the pristine- and polluted-population masses. Milone et al.(2020) found that $F_{1P}$ anticorrelates more tightly with the polluted-population present-day mass, $m_{2P,prst}$, than with the cluster total mass, $m_{prst}$. In contrast, $F_{1P}$ anticorrelates poorly with the pristine-population current mass, $m_{1P,prst}$. We show the loose anticorrelation between $F_{1P}$ and $m_{1P,prst}$ to result from a roughly constant pristine-population mass among clusters as they start their long-term evolution in the Galactic tidal field. As for the tight anticorrelation between $m_{2P,prst}$ and $F_{1P}$, it stems from the initially shallow relation between $m_{2P}$ and $F_{1P}$ (Fig.4). Clusters of the Large and Small Magellanic Clouds (LMC and SMC, respectively) appear to behave unexpectedly with respect to each other. For a given $F_{1P}$, LMC clusters are more massive than SMC clusters despite their enduring a stronger tidal field (Fig.6). This is opposite to how the Galactic outer- and inner-halo clusters behave (bottom panels of Figs 1-3). The explanation may lie in cluster formation conditions. Finally, we wonder whether the single-population clusters NGC~419 and Rup~106 formed as multiple-populations clusters.

Accreting massive black holes (MBHs, with M$_\bullet > 10^3$ M$_{\odot}$) are known for their panchromatic emission, spanning from radio to gamma rays. While MBHs accreting at significant fractions of their Eddington rate are readily detectable, those accreting at much lower rates in radiatively inefficient modes often go unnoticed, blending in with other astrophysical sources. This challenge is particularly relevant for gas-starved MBHs in external galaxies and those possibly wandering in the Milky Way. We present HOLESOM, a machine learning-powered tool based on the Self-Organizing Maps (SOMs) algorithm, specifically designed to identify slowly-accreting MBHs using sparse photometric data. Trained on a comprehensive set of $\sim$ 20, 000 spectral energy distributions (SEDs), HOLESOM can (i) determine if the photometry of a source is consistent with slowly-accreting MBHs and (ii) estimate its black hole mass and Eddington ratio, including uncertainties. We validate HOLESOM through extensive tests on synthetic data and real-world cases, including Sagittarius A* (Sgr A*), demonstrating its effectiveness in identifying slowly-accreting MBHs. Additionally, we derive analytical relations between radio and X-ray luminosities to further constrain physical parameters. The primary strength of HOLESOM lies in its ability to accurately identify MBH candidates, which can then be targeted for follow-up photometric and spectroscopic observations. Fast and scalable, HOLESOM offers a robust framework for automatically scanning large multi-wavelength datasets, making it a valuable tool for unveiling hidden MBH populations in the local Universe.

Owing to their similarities with giant exoplanets, brown dwarf companions of stars provide insights into the fundamental processes of planet formation and evolution. From their orbits, several brown dwarf companions are found to be more massive than theoretical predictions given their luminosities and the ages of their host stars (e.g. Brandt et al. 2021, Cheetham et al. 2018, Li et al. 2023). Either the theory is incomplete or these objects are not single entities. For example, they could be two brown dwarfs each with a lower mass and intrinsic luminosity (Brandt et al. 2021, Howe et al. 2024). The most problematic example is Gliese 229 B (Nakajima et al. 1995, Oppenheimer et al. 1995), which is at least 2-6 times less luminous than model predictions given its dynamical mass of $71.4\pm0.6$ Jupiter masses ($M_{\rm Jup}$) (Brandt et al. 2021). We observed Gliese 229 B with the GRAVITY interferometer and, separately, the CRIRES+ spectrograph at the Very Large Telescope. Both sets of observations independently resolve Gliese 229 B into two components, Gliese 229 Ba and Bb, settling the conflict between theory and observations. The two objects have a flux ratio of $0.47\pm0.03$ at a wavelength of 2 $\mu$m and masses of $38.1\pm1.0$ and $34.4\pm1.5$ $M_{\rm Jup}$, respectively. They orbit each other every 12.1 days with a semimajor axis of 0.042 astronomical units (AU). The discovery of Gliese 229 BaBb, each only a few times more massive than the most massive planets, and separated by 16 times the Earth-moon distance, raises new questions about the formation and prevalence of tight binary brown dwarfs around stars.

We introduce a new, open-source, Python module for the acquisition and processing of archival data from many X-ray telescopes - Democratising Archival X-ray Astronomy (hereafter referred to as DAXA). Our software is built to increase access to, and use of, large archives of X-ray astronomy data; providing a unified, easy-to-use, Python interface to the disparate archives and processing tools. We provide this interface for the majority of X-ray telescopes launched within the last 30 years. This module enables much greater access to X-ray data for non-specialists, while preserving low-level control of processing for X-ray experts. It is useful for identifying relevant observations of a single object of interest but it excels at creating multi-mission datasets for serendipitous or targeted studies of large samples of X-ray emitting objects. The management and organization of datasets is also made easier; DAXA archives can be version controlled and updated if new data become available. Once relevant observations are identified, the raw data can be downloaded (and optionally processed) through DAXA, or pre-processed event lists, images, and exposure maps can be downloaded if they are available. X-ray observations are perfectly suited to serendipitous discoveries and archival analyses, and with a decade-long `X-ray desert' potentially on the horizon archival data will take on even greater importance; enhanced access to those archives will be vital to the continuation of X-ray astronomy.

The Large Magellanic Cloud (LMC) is the nearest massive galaxy to the Milky Way. Its circumgalactic medium is complex and multi-phase, containing both stripped HI structures like the Magellanic Stream and Bridge, and a diffuse warm corona seen in high-ion absorption. We analyze 28 AGN sightlines passing within 35 kpc of the LMC with archival HST/COS spectra to characterize the cool (T\approx10^4$ K) gas in the LMC CGM, via new measurements of UV absorption in six low ions (OI, FeII, SiII, AlII, SII, and NiII) and one intermediate ion (SiIII). We show that a declining column-density profile is present in all seven ions, with the low-ion profiles having a steeper slope than the high-ion profiles in CIV and SiIV reported by Krishnarao et al. 2022. Crucially, absorption at the LMC systemic velocity is only detected (in all ions) out to 17 kpc. Beyond this distance, the gas has a lower velocity and is associated with the Magellanic Stream. These results demonstrate that the LMC's CGM is composed of two distinct components: a compact inner halo extending to 17 kpc, and a more extended stripped region associated with the Stream. The compactness and truncation of the LMC's inner CGM agree with recent simulations of ram-pressure stripping of the LMC by the Milky Way's extended corona.

Over the last decades, evolutionary population synthesis models have powered an unmatched leap forward in our understanding of galaxies. From dating the age of the first galaxies in the Universe to detailed measurements of the chemical composition of nearby galaxies, the success of this approach built upon simple stellar population (SSP) spectro-photometric models is unquestionable. However, the internal constraints inherent to the construction of SSP models may hinder our ability to analyze the integrated spectra of galaxies in situations where the SSP assumption does not sufficiently hold. Thus, here we revisit the possibilities of fitting galaxy spectra as a linear combination of stellar templates without assuming any a priori knowledge on stellar evolution. We showcase the sensitivity of this alternative approach to changes in the stellar population properties, in particular the direct connection to variations in the stellar initial mass function, as well as its advantages when dealing with non-canonical integrated populations and semi-resolved observations. Furthermore, our analysis demonstrates that the absorption spectra of galaxies can be used to independently constrain tellar evolution theory beyond the limited conditions of the solar neighborhood.

Asymptotic Safety is a promising framework towards the understanding, in a non-perturbative way, of Quantum Gravity. It treats the Newton's constant G_N and the cosmological constant \Lambda as running coupling of an effective action. At the phenomenological point of view the values of G_N and \Lambda depend on the energy density of the system under consideration. This fact has interesting astrophysical and cosmological consequences. The present work discusses the effects on the distance measurements from SNIa light curves and on the strong and weak gravitational lensing measurements.

A Galactic core-collapse supernova (CCSN) is likely to be observed in neutrino detectors around the world minutes to hours before the electromagnetic radiation arrives. The SNEWS2.0 network of neutrino and dark matter detectors aims to use the relative arrival times of the neutrinos at the different experiments to point back to the supernova. One of the simplest methods to calculate the CCSN direction is to use the first neutrino events detected through the inverse beta decay (IBD) process, $\overline{\nu}_e p\rightarrow e^+n$. We will consider neutrino detectors sensitive to IBD interactions with low backgrounds. The difference in signal arrival times between a large and a small detector will be biased, however, with the first event at the smaller detector, on average, arriving later than that at the larger detector. This bias can be mitigated by using these first events in a data-driven approach without recourse to simulations or models. In this article, we demonstrate this method and its uncertainty estimate using pairs of detectors of different sizes and with different supernova distances. Finally, we use this method to calculate probability skymaps using four detectors (Super-Kamiokande, JUNO, LVD, and SNO+) and show that the calculated probabilities yield appropriate confidence intervals. The resulting skymaps should be useful for the multi-messenger community to follow up on the SNEWS2.0 Galactic CCSN neutrino alert.

We present two epochs of radial velocities of the first imaged T dwarf Gliese 229B obtained with Keck/NIRSPEC. The two radial velocities are discrepant with one another, and with the radial velocity of the host star, at $\approx$$11\sigma$ significance. This points to the existence of a previously postulated, but as-yet undetected, massive companion to Gl 229B; we denote the two components as Gl 229Ba and Gl 229Bb. We compute the joint likelihood of the radial velocities to constrain the period and mass of the secondary companion. Our radial velocities are consistent with an orbital period between a few days and $\approx$60 days, and a secondary mass of at least $\approx$15\,$M_{\rm Jup}$ and up to nearly half the total system mass of Gl 229B. With a significant fraction of the system mass in a faint companion, the strong tension between Gl 229B's dynamical mass and the predictions of evolutionary models is resolved.

The Wide Field Monitor (WFM) is one of the three instruments on the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X) mission, which was proposed in response to the NASA 2023 call for a probe class mission. The WFM is a coded-mask camera system that would be the most scientifically capable wide-angle monitor ever flown. The field of view covers one third of the sky, to 50 percent mask coding, and the energy sensitivity is 2 to 50 keV. The WFM is designed to identify new X-ray transients and to capture spectral and timing changes in known sources with data of unprecedented quality. Science applications cover diverse classes, in including X-ray bursts that coincide with gravitational wave detections, gamma ray bursts and their transition from prompt emission to afterglow, subluminous GRBs that may signal shock breakout in supernovae, state transitions in accreting compact objects and their jets, bright flares in fast X-ray transients, accretion onset in transitional pulsars, and coronal flares from many types of active stars.

An analysis of the Ultra Long Period Cepheids (ULPs) properties could significantly contribute to understanding the Hubble constant tension, e.g. the current discrepancy between determinations based on local distance indicators and those relying on cosmic microwave background measurements. These highly luminous variables are observable beyond 100 Mpc, so if they were confirmed to behave as standard candles, they would allow us a direct measurement of cosmological distances without any secondary distance indicator, thus reducing potential systematic errors in the calibration of the cosmic distance scale. This paper presents an analysis of the largest known sample of 73 ULPs, including 15 objects in nearby galaxies, with new accurate and homogeneous photometry obtained by Gaia DR3, and a new object, in our Galaxy, identified as Long Period Variable in Gaia DR3, but recently reclassified as ULP. The obtained results suggest that, by improving photometric accuracy, the ULP Period-Wesenheit relation shows a smaller dispersion than that obtained in literature and is in better agreement with the Classical Cepheid (CC) one, supporting the hypothesis that ULPs are the extension of the CCs at higher period, mass and luminosity. However, to reach this aim, it is necessary to enrich the sample with high-quality data. The Rubin-LSST survey offers the possibility to achieve this thanks to its photometric characteristics and time extension. In particular, we will explore the capabilities of the Rubin-LSST survey to recover ULP theoretical light-curves by using a new tool called PulsationStarRecovery, built by our group for this type of analysis.

Adaptive optics (AO) systems are crucial for high-resolution astronomical observations by compensating for atmospheric turbulence. While laser guide stars (LGS) address high-order wavefront aberrations, natural guide stars (NGS) remain vital for low-order wavefront sensing (LOWFS). Conventional NGS-based methods like Shack-Hartmann sensors have limitations in field of view, sensitivity, and complexity. Focal plane wavefront sensing (FPWFS) offers advantages, including a wider field of view and enhanced signal-to-noise ratio, but accurately estimating low-order modes from distorted point spread functions (PSFs) remains challenging. We propose an AI-powered FPWFS method specifically for low-order mode estimation in infrared wavelengths. Our approach is trained on simulated data and validated on on-telescope data collected from the Keck I adaptive optic (K1AO) bench calibration source in K-band. By leveraging the enhanced signal-to-noise ratio in the infrared and the power of AI, our method overcomes the limitations of traditional LOWFS techniques. This study demonstrates the effectiveness of AI-based FPWFS for low-order wavefront sensing, paving the way for more compact, efficient, and high-performing AO systems for astronomical observations.

We use a custom-made calibrator to measure individual detectors' polarization angles of BICEP3, a small aperture telescope observing the cosmic microwave background (CMB) at 95GHz from the South Pole. We describe our calibration strategy and the statistical and systematic uncertainties associated with the measurement. We reach an unprecedented precision for such measurement on a CMB experiment, with a repeatability for each detector pair of $0.02°$. We show that the relative angles measured using this method are in excellent agreement with those extracted from CMB data. Because the absolute measurement is currently limited by a systematic uncertainty, we do not derive cosmic birefringence constraints from BICEP3 data in this work. Rather, we forecast the sensitivity of BICEP3 sky maps for such analysis. We investigate the relative contributions of instrument noise, lensing, and dust, as well as astrophysical and instrumental systematics. We also explore the constraining power of different angle estimators, depending on analysis choices. We establish that the BICEP3 2-year dataset (2017--2018) has an on-sky sensitivity to the cosmic birefringence angle of $\sigma = 0.078°$, which could be improved to $\sigma = 0.055°$ by adding all of the existing BICEP3 data (through 2023). Furthermore, we emphasize the possibility of using the BICEP3 sky patch as a polarization calibration source for CMB experiments, which with the present data could reach a precision of $0.035°$. Finally, in the context of inflation searches, we investigate the impact of detector-to-detector variations in polarization angles as they may bias the tensor-to-scalar ratio r. We show that while the effect is expected to remain subdominant to other sources of systematic uncertainty, it can be reliably calibrated using polarization angle measurements such as the ones we present in this paper.

The first detection of the quasiperiodic ultrafast outflow in the ASASSN-20qc system was reported by Pasham et al. (2024). The outflow is revealed in the soft X-ray spectra as an absorption feature, which is periodically enhanced every $ \sim 8.3$ days. The periodic nature of the ultrafast outflow is best explained by an orbiting massive perturber, most likely an intermediate-mass black hole (IMBH), that is inclined with respect to the accretion flow around the primary supermassive black hole (SMBH). In this way, the orbiting body pushes the disc gas into the outflow funnel, where it is accelerated by the ordered magnetic field (Suková et al. 2021). Quasiperiodic ultrafast outflows (QPOuts) are thus a novel phenomenon that can help reveal new extreme-/intemediate-mass ratio inspiral (EMRI/IMRI) candidates that otherwise do not exhibit significant periodic patterns in the continuum flux density.

The Magellanic Stream, a tidal tail of diffuse gas falling onto the Milky Way, formed by interactions between the Small and Large Magellanic Clouds, is primarily composed of neutral atomic hydrogen (HI). The deficiency of dust and the diffuse nature of the present gas make molecular formation rare and difficult, but if present, could lead to regions potentially suitable for star formation, thereby allowing us to probe conditions of star formation similar to those at high redshifts. We search for HCO$^+$, HCN, HNC, and C$_2$H using the highest sensitivity observations of molecular absorption data from the Atacama Large Millimeter Array to trace these regions, comparing with HI archival data to compare these environments in the Magellanic Stream to the HI column density threshold for molecular formation in the Milky Way. We also compare the line of sight locations with confirmed locations of stars, molecular hydrogen, and OI detections, though at higher sensitivities than the observations presented here. We find no detections to a 3$\sigma$ significance, despite four sightlines having column densities surpassing the threshold for molecular formation in the diffuse regions of the Milky Way. Here we present our calculations for the upper limits of the column densities of each of these molecular absorption lines, ranging from $3 \times 10^{10}$ to $1 \times 10^{13}$ cm$^{-2}$. The non-detection of HCO$^+$ suggests that at least one of the following is true: (i) $X_{HCO^+, \mathrm{MS}}$ is significantly lower than the Milky Way value; (ii) that the widespread diffuse molecular gas observed in the Milky Way's diffuse ISM does not have a direct analog in the MS; (iii) the HI-to-H$_2$ transition occurs in the MS at a higher surface density in the MS than in the LMC or SMC; or (iv) molecular gas exists in the MS, but only in small, dense clumps.

Classical nova eruptions result from thermonuclear-powered runaways in, and ejection of, the hydrogen-rich envelopes of white dwarf stars accreted from their close binary companions. Novae brighten to up to 1,000,000 solar luminosities, and recur thousands of times over their lifetimes spanning several billion years. Between eruptions, mass transfer from the donor star to the white dwarf proceeds via an accretion disk unless the white dwarf possesses a strong magnetic field which can partially or totally disrupt the disk. In that case, accretion is focussed by the white dwarf's magnetic field towards its magnetic poles. Optical spectroscopy and interferometric radio maps demonstrate the presence of bipolar jets, typically arcsec in angular size, and orders of magnitude smaller than one parsec in linear size, during the days to months after nova eruptions. These jets expel collimated matter from the white dwarfs in nova binary stars, but well-resolved images of them are lacking. Here we report the Condor telescope's detection of a hitherto unknown, highly resolved and braided jet, three degrees (at least 25 parsecs) in length. The jet originates at the white dwarf of the old nova GK Persei (nova Per 1901 CE). It precesses on a ~ 3600 yr timescale, and must be at least 7200 years old. Detected across four decades of wavelength, the jet's ultimate energy source is likely the strong accretion shocks near the white dwarf's magnetic poles.

Over a hundred gravitational-wave signals have now been detected from the mergers of black holes and neutron stars, but other sources of gravitational waves have not yet been discovered. Some of the most violent explosive events in the Universe are predicted to emit bursts of gravitational waves, and may result in the next big multi-messenger discovery. Gravitational-wave burst signals often have an unknown waveform shape, and unknown gravitational-wave energy, due to unknown or very complicated progenitor astrophysics. Potential sources of gravitational-wave bursts include core-collapse supernovae, cosmic strings, fast radio bursts, eccentric binary systems, and gravitational-wave memory. In this review, we discuss the astrophysical properties of the main predicted sources of gravitational-wave bursts, and the known features of their gravitational-wave emission. We summarise their future detection prospects, and discuss the challenges of searching for gravitational-wave burst signals and interpreting the astrophysics of the source.

MAVIS (the MCAO-Assisted Visible Imager and Spectrograph), planned for the VLT Adaptive Optics Facility, represents an innovative step in Multi-Conjugate Adaptive Optics (MCAO) systems, particularly in its operation at visible wavelengths and anticipated contributions to the field of astronomical astrometry. Recognizing the crucial role of high-precision astrometry in realizing science goals such as studying the dynamics of dense starfields, this study focuses on the challenges of advancing astrometry with MAVIS to its limits, as well as paving the way for further enhancement by incorporating telemetry data as part of the astrometric analysis. We employ MAVISIM, Superstar, and DAOPHOT to simulate both MAVIS imaging performance and provide a pathway to incorporate telemetry data for precise astrometry with MAVIS. Photometry analyses are conducted using the Superstar and DAOPHOT platforms, integrated into a specifically designed pipeline for astrometric analysis in MCAO settings. Combining these platforms, our research aims to elucidate the impact of utilizing telemetry data on improving astrometric precision, potentially establishing new methods for ground-based AO-assisted astrometric analysis. This endeavor not only sheds light on the capabilities of MAVIS but also paves the way for advancing astrometry in the era of next-generation MCAO-enabled giant telescopes.

The Thirty Meter Telescope will use a sophisticated adaptive optics system called NFIRAOS. This system utilizes two deformable mirrors conjugate to 0 km and 11.2 km to apply a Multi-Conjugate Adaptive Optics (MCAO) correction over a 2 arcminute field of view. DM0 and DM11 have 63 and 75 actuators across their respective diameters. To study the behavior of these mirrors, we have developed a low-cost, very high-order Shack-Hartmann Wavefront Sensor (WFS). We will use our WFS to calibrate the flatness of the DMs and measure the influence functions of the actuators. NFIRAOS is cooled to reduce the thermal emissivity of optical surfaces visible to the science detectors, so we will also measure the behaviour of the DMs in both warm and cold environments. As the cold chamber is prone to vibrations, a WFS is preferred to a phase-shifting interferometer. Our design was driven by the need to be able to evaluate the DM surface between the actuators, which led to the requirement of at least 248 sub apertures across the diameter. The largest commercially available Shack-Hartmann WFS has only 128 sub-apertures across the diameter, which is not enough to properly sample these DMs. Furthermore, the designed sensor is able to record the wavefront at 50 FPS (50 times per second) at full resolution. To fabricate this WFS, we used a commercial off-the-shelf CMOS detector, camera lens, and lens let array, which kept the total cost less than 20K USD. Here we present the design and performance characteristics of this device.

GRB170817A (also GW170817) became the first binary neutron star (BNS) merger event detected via gravitational waves and electromagnetic signals. Over the next 4 years, various multiband observations have led to re-imagine the various short Gamma Ray Burts (sGRB) and interstellar medium interaction models. While these models successfully explain the observed afterglow until ~900 days, a re-brightening or excess flux was observed in the 1keV X-Ray band after ~1000 days. In this study, we re-evaluate the jet parameters using the new observations (until ~1234 days) with a boosted fireball jet model. We study the observable effects of the counter-jet for GRB170817A, using our new afterglow code, Firefly. Our results show that it is indeed possible for the observed excess to coincide with the emissions from the counter-jet (~800 days). We also computed an empirical scaling law between the jet and counter-jet peak emission timescales and the observer angle. The Firefly code can also track the simulated object through the observers' sky and numerically model the apparent motion. The calculated apparent motion (~2.6c) does not match the observed apparent motion (7.5c to 5.2c). Hence we conclude, the excess flux of GRB170817A may not be associated with the counter jet; however, it is not enough to reject this hypothesis from the traditional counter jet visibility time scales, which predicts > 5000 days. The apparent motion, combined with the multi-band lightcurves, is needed to break degeneracy between geometrical parameters and the microphysical parameters of the afterglow.

We present significant upgrades to the VAMPIRES instrument, a visible-light (600 nm to 800 nm) high-contrast imaging polarimeter integrated within SCExAO on the Subaru telescope. Key enhancements include new qCMOS detectors, coronagraphs, polarization optics, and a multiband imaging mode, improving sensitivity, resolution, and efficiency. These upgrades position VAMPIRES as a powerful tool for studying sub-stellar companions, accreting protoplanets, circumstellar disks, stellar jets, stellar mass-loss shells, and solar system objects. The instrument achieves angular resolutions from 17 mas to 21 mas and Strehl ratios up to 60\%, with 5$\sigma$ contrast limits of $10^{\text{-}4}$ at 0.1'' to $10^{\text{-}6}$ beyond 0.5''. We demonstrate these capabilities through spectro-polarimetric coronagraphic imaging of the HD 169142 circumstellar disk, ADI+SDI imaging of the sub-stellar companion HD 1160B, narrowband H$\alpha$ imaging of the R Aqr emission nebula, and spectro-polarimetric imaging of Neptune.

We investigated two consecutive solar eruption events in the solar active region (AR) 12994 at the solar eastern limb on 2022 April 15. We found that the flare loops formed by the first eruption were involved in the second eruption. During the initial stage of the second flare, the middle part of these flare loops (E-loops) erupted outward along with the flux ropes below, while the parts of the flare loops (I-loops1 and I-loops2) on either side of the E-loops first rose and then contracted. Approximately 1 hour after the eruption, the heights of I-loops1 and I-loops2 decreased by 9 Mm and 45 Mm, respectively, compared to before the eruption. Their maximum descent velocities were 30 km/s and 130 km/s, respectively. The differential emission measure (DEM) results indicate that the plasma above I-loops1 and I-loops2 began to be heated about 23 minutes and 44 minutes after the start of the second flare, respectively. Within 20 minutes, the plasma temperature in these regions increased from ~3 MK to 6 MK. We proposed an adiabatic heating mechanism that magnetic energy would be converted into thermal and kinetic energy when the pre-stretched loops contract. Our calculations show that the magnetic energy required to heat the two high-temperature regions are 10^29-10^30 erg, which correspond to a loss of field strength of 2-3 G.

Dense clumps distributed along filaments are the immediate medium for star formation. Kinematic properties of the clumps, such as velocity gradient and angular momentum, combined with filament orientation, provide important clues to the formation mechanism of filament-clump configurations and the role of filaments in star formation. By cross-matching the Milky Way atlas for linear filaments and the Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium (SEDIGISM) 13CO (2-1) data, we aim to derive the velocity gradient and its direction, the specific angular momentum (J/M), and the ratio (\beta) between the rotational energy and gravitational energy of clumps, as well as to investigate the alignment between clump rotation and filament orientation. We found a monotonic increase in J/M as a function of clump size (R), following a power-law relation J/M~\propto~R^{1.5\pm0.2}. The ratio \beta ranges from 1.1~\times~10^{-5} to 0.1, with a median value 1.0~\times~10^{-3}, suggesting that clump rotation provides insignificant support against gravitational collapse. The distribution of the angle between clump rotation and natal filament orientation is random, indicating that the clumps' rotational axes have no discernible correlation with the orientation of their hosting filaments. Counting only the most massive clump in each filament also finds no alignment between clump rotation and filament orientation.

Relativistic accretors are cosmic objects that pull matter from their surroundings at speeds almost equal to the light's speed. Because of the tremendous gravitational force from the accretors and the angular momentum of infalling material, which often result in discs of gas and dust that are heated to extremely high temperatures. We encounter strong radiation throughout the electromagnetic spectrum, including intense X-rays. The X-ray view provides a unique window into the behavior of accretors. In this review, we discuss different accretors, particularly binaries, and their origin, involved mechanisms, and properties, including energy spectra and the variability of X-rays from the accretors. This X-ray perspective gives a unique insight into the evolution and connections of these systems with their environment. Future research in this area is necessary to fully understand the process underlying X-ray emission from relativistic accretors.

Ethanolamine (NH$_2$CH$_2$CH$_2$OH, EA) has been identified in the gas phase of the ISM within molecular clouds. Although EA has not been directly observed in the molecular ice phase, a solid state formation mechanism has been proposed. However, the current literature lacks an estimation of the infrared band strengths of EA ices. We conducted an experimental investigation of solid EA ice at low temperatures to ascertain its infrared band strengths, phase transition temperature, and multilayer binding energy. The commonly used laser interferometry method was not applied. Infrared band strengths were determined using three distinct methods. The obtained lab spectrum of EA was compared with the publicly available MIRI MRS James Webb Space Telescope observations toward a low mass protostar. The phase transition temperature for EA ice falls within the range of 175 to 185 K. Among the discussed methods, the simple pressure gauge method provides a reasonable estimate of band strength. We derive a band strength value of about $1\times10^{-17}$ cm molecule$^{-1}$ for the NH$_2$ bending mode in the EA molecules. Additionally, temperature-programmed desorption analysis yielded a multilayer desorption energy of 0.61$\pm$0.01 eV. By comparing the laboratory data documented in this study with the JWST spectrum of the low mass protostar IRAS 2A, an upper-limit for the EA ice abundances was derived.

Previous research on Lorentz invariance violation in photons from gamma-ray bursts (GRBs) suggested a scenario where multi-GeV photons could be emitted before lower-energy photons at the GRB source frame. This implies the existence of a new preburst phase in addition to the traditionally identified prompt and afterglow stages observed in earlier studies. In this study, we present direct evidence for this novel preburst phase in gamma-ray bursts based on recent observations of GRB 221009A. Our analysis leverages data from the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) detectors of the Fermi Gamma-ray Space Telescope (FGST), as well as data from the KM2A detector of the Large High Altitude Air-shower Observatory (LHAASO).

Reproducing color-magnitude diagrams (CMDs) of star-resolved galaxies is one of the most precise methods for measuring the star formation history (SFH) of nearby galaxies back to the earliest time. The upcoming big data era poses challenges to the traditional numerical technique in its capacity to deal with vast amounts of data, which motivates us to explore the feasibility of employing machine learning networks in this field. In this study, we refine the synthetic CMD method with a state-of-the-art theoretical stellar evolution model to simulate the properties of stellar populations, incorporate the convolutional neural network (CNN) in the fitting process to enhance the efficiency, and innovate the initial stellar mass estimation to improve the flexibility. The fine-tuned deep learning network, named \texttt{SFHNet}, has been tested with synthetic data and further validated with photometric data collected from the Hubble Space Telescope (\textit{HST}). The derived SFHs are largely in accordance with those reported in the literature. Furthermore, the network provides detailed insights into the distribution of stellar density, initial stellar mass, and star formation rate (SFR) over the age-metallicity map. The application of the deep learning network not only measures the SFH accurately but also enhances the synthetic CMD method's efficiency and flexibility, thereby facilitating a more comprehensive and in-depth understanding of nearby galaxies.

Analysing stellar parameters and abundances from nearly one million Gaia DR3 Radial Velocity Spectrometer (RVS) spectra poses challenges due to the limited spectral coverage (restricted to the infrared Ca II triplet) and variable signal-to-noise ratios of the data. To address this, we use The Cannon, a data-driven method, to transfer stellar parameters and abundances from the GALAH Data Release 4 (DR4; R ~ 28,000) catalogue to the lower resolution Gaia DR3 RVS spectra (R ~ 11,500). Our model, trained on 14,484 common targets, predicts parameters such as Teff, log g, and [Fe/H], along with several other elements across approximately 800,000 Gaia RVS spectra. We utilise stars from open and globular clusters present in the Gaia RVS catalogue to validate our predicted mean [Fe/H] with high precision (~0.02-0.10 dex). Additionally, we recover the bimodal distribution of [Ti/Fe] versus [Fe/H], reflecting the high and low alpha-components of Milky Way disk stars, demonstrating The Cannon's capability for accurate stellar abundance determination from medium-resolution Gaia RVS spectra. The methodologies and resultant catalogue presented in this work highlight the remarkable potential of the RVS dataset, which by the end of the Gaia mission will comprise spectra of over 200 million stars.

Supernova blasts envelop many surrounding stellar systems, transferring kinetic energy to small bodies in the systems. Geologic evidence from $^{60}\rm Fe$ points to recent nearby supernova activity within the past several Myr. Here, we model the transfer of energy and resulting orbital changes from these supernova blasts to the Oort Cloud, the Kuiper belt, and Saturn's Phoebe ring. For the Oort Cloud, an impulse approximation shows that a 50 pc supernova can eject approximately half of all objects less than 1 cm while altering the trajectories of larger ones, depending on their orbital parameters. For stars closest to supernovae, objects up to $\sim$100 m can be ejected. Turning to the explored solar system, we find that supernovae closer than 50 pc may affect Saturn's Phoebe ring and can sweep away Kuiper belt dust. It is also possible that the passage of the solar system through a dense interstellar cloud could have a similar effect; a numerical trajectory simulation shows that the location of the dust grains and the direction of the wind (from a supernova or interstellar cloud) has a significant impact on whether or not the grains will become unbound from their orbit in the Kuiper belt. Overall, nearby supernovae sweep micron-sized dust from the solar system, though whether the grains are ultimately cast towards the Sun or altogether ejected depends on various factors. Evidence of supernova-modified dust grain trajectories may be observed by New Horizons, though further modeling efforts are required.

Molecular gas is the key probe for the complex interaction between the accretion of black holes and star formation of the host galaxy of active galactic nuclei (AGN). The molecular gas discovered around the AGN indicates that this gas is providing fuel for the AGN. According to the theoretical model of the relativistic jet, the spin of a black hole enhances the relativistic jet of AGN. The spin of the black hole is used as an indicator of AGN activity. Therefore, we study the relationship between the activity of AGN and molecular gas. We find a significant strong correlation between molecular gas fraction and CO luminosity and black hole spin for the early-type galaxies. However, there is no correlation between molecular gas fraction and CO luminosity and black hole spin for the late-type galaxies. These results indicate that the spin of black holes mainly regulates the accretion of molecular gas in massive early-type galaxies. The activity of AGN depends on the amount of gas.

We investigate strong gravitational lensing by two static black hole models (Model-1 and Model-2) within the Effective Quantum Gravity (EQG) framework, characterized by mass $M$ and parameter $\zeta$. For $\zeta = 0$, they reduce to the Schwarzschild solution, and depending on the parameters, they describe black holes with an event and Cauchy horizon (Model-1), a single horizon (Model-2), or no horizons. Using SMBHs Sgr A* and M87* as lenses and integrating theoretical predictions with recent EHT data, we identify significant differences in lensing signatures due to quantum corrections. For Model-1, the deviations of the lensing observables: $|\delta\theta_{\infty}|$ of black holes in EQG from Schwarzschild black hole, for SMBHs Sgr A* and M87, can reach as much as $1.75~\mu$as and $1.32~\mu$as, while $|\delta s|$ is about $30.12$~nas for Sgr A* and $22.63$~nas for M87*. The flux ratio of the first image to all subsequent packed images indicates that EQG black hole images are brighter than their Schwarzschild counterparts, with a deviation in the brightness ratio $|\delta r_{mag}|$ reaching up to 2.02. The time delays between the second and first images, denoted $|\delta T_{2,1}|$, exhibit substantial deviations from the GR counterpart, reaching up to 1.53 min for Sgr A* and 1159.9 min for M87*. The EHT constraints on $\theta_{sh}$ of Sgr A* and M87* within the $1\sigma$ region limit the parameters $\zeta$. Our analysis concludes that EQG black holes are consistent with the EHT observations within this finite space.

Extreme-mass-ratio inspirals (EMRIs), consisting of a massive black hole and a stellar compact object, are one of the most important sources for space-borne gravitational wave detectors like TianQin. Their population study can be used to constrain astrophysical models that interpret the EMRI formation mechanisms. In this paper, as a first step, we employ a parametrization method to describe the EMRI population model in the loss cone formation channel. This approach, however, can be extended to other models such as the accretion disc driven formation channel. We present the phenomenological characteristic of the MBH mass, spin, and redshift distributions. Then, we investigate the posterior distribution of the hyper-parameters that describe this population model. Our results show that TianQin could recover almost all the posterior of the hyper-parameters within $1\sigma$ confidence interval. With one hundred detectable EMRI events, the hyper-parameters $\alpha_1, \alpha_2, b$, which describe the MBH mass distribution, could be measured with an accuracy of $37\%$, $24\%$, and $3\%$, respectively. The hyper-parameters $\mu_z$, and $\sigma_z$, which describe the redshift distribution, have $\mu_z$ above the detectable range of TianQin, and $\sigma_z$ measured with an accuracy of $14.5\%$. With this estimation accuracy, the EMRI population characteristics can be effectively demonstrated, potentially serving as evidence for EMRI formation in the future studies. Furthermore, with an increasing number of detectable events, the parameter estimation for the hyper-parameters will improve and the confidence intervals will be narrowed.

We use RR Lyrae stars identified by the Gaia third data release (DR3) to explore the outskirts of the Ursa Minor (UMi) dwarf spheroidal galaxy (dSph) and update the census of its variable star population. We adopted different tools based on the Gaia DR3 astrometric and photometric data (proper motions, Period-Wesenheit-Metallicity relations, spatial distribution, colour-magnitude diagram and stellar isochrone fitting) to discriminate between different types of variable stars, and to identify UMi members. We find a total of 129 RR Lyrae stars and Anomalous Cepheids (ACs) that belong to UMi. We report 47 new RR Lyrae stars (46 bona fide and 1 candidate) and 5 new ACs (4 bona fide and 1 candidate), including new possible members in the extreme periphery of the galaxy at a distance of $\sim$ 12 half-light radii. We reclassified 13 RR Lyrae stars identified by the Gaia DR3 Specific Object Study pipeline for Cepheids and RR Lyrae stars (SOS Cep$\&$RRL), using data from the literature and Gaia astrometry and photometry. Specifically, we assigned these 13 DR3 RR Lyrae stars to 10 Anomalous Cepheids and 3 double-mode RR Lyrae (RRd), respectively. From the average luminosity of the RR Lyrae stars we derive for UMi a distance modulus of $(m-M)_{0}=$ 19.23 $\pm$ 0.09 mag in excellent agreement with the literature. Finally, we investigated whether some of UMi's variable stars might be members of the ultra-faint stellar cluster Mu{ñ}oz~1, that lies at a projected distance of 45 arcmin from UMi's centre. Based on the properties of the variable stars (distances, colours and metallicities), we find unlikely that they belong to the cluster.

The properties, impact, and fate of hot stars cannot be understood without considering their winds. Revealed to be an almost ubiquitous phenomenon in the regime of massive stars, the winds of hot stars arise from a complex physical mechanism that still provides a major challenge for our understanding of massive stars. Different flavours of hot stars vary significantly in their winds with current evolution models still having problems to connect the zoo of observed phenomena. Moreover, the driving of hot star winds is inherently connected to the opacities arising from spectral line transitions, making the properties and strength of the winds strongly dependent on metallicity and changing them over cosmic time. In these proceedings, the current status in our understanding of hot star winds is briefly reviewed and recent progress in our perception of hot star winds and the consequences for the evolution of massive stars are presented. A particular emphasis is given on current efforts with hydrodynamically-consistent atmosphere models towards a better description of radiatively-driven mass loss of Wolf-Rayet stars with remaining hydrogen envelopes.

Forthcoming cosmic shear surveys will make precise measurements of the matter density field down to very small scales, scales which are dominated by baryon feedback. The modelling of baryon feedback is crucial to ensure unbiased cosmological parameter constraints; the most efficient approach is to use analytic models, but these are limited by how well they can capture the physics of baryon feedback. We investigate the fitting and residual errors of various baryon feedback models to a suite of hydrodynamic simulations, and propagate these to cosmological parameter constraints for cosmic shear. We present an alternative formalism to binary scale-cuts through the use of a theoretical error covariance, which is a well-motivated alternative using errors in the power spectrum modelling itself. We depart from previous works by modelling baryonic feedback errors directly in the matter power spectrum, which is the natural basis to do so and thus preserves information in the lensing kernels. When including angular multipoles up to $\ell_{\mathrm{max}} = 5000$, and assuming Euclid-like survey properties, we find that even multi-parameter models of baryon feedback can introduce significant levels of bias. In contrast, our theoretical error reduces the bias in $\Omega_{\mathrm{m}}$ and $S_{8}$ to acceptable levels, with only a modest increase in parameter variances. The theoretical error approach bypasses the need to directly determine the per-bin $\ell_{\mathrm{max}}$ values, as it naturally suppresses the biassing small-scale information. We also present a detailed study of how flexible HMCode-2020, a widely-used non-linear and baryonic feedback model, is at fitting a range of hydrodynamical simulations.

PSR J0737$-$3039A/B is unique among double neutron star systems. Its near-perfect edge-on orbit causes the fast spinning pulsar A to be eclipsed by the magnetic field of the slow spinning pulsar B. Using high-sensitivity MeerKAT radio observations combined with updated constraints on the system geometry, we studied the impact of these eclipses on the incident polarization properties of pulsar A. Averaging light curves together after correcting for the rotation of pulsar B revealed enormous amounts of circular polarization and rapid changes in the linear polarization position angle, which occur at phases where emission from pulsar A is partially transmitted through the magnetosphere of pulsar B. These behaviours confirm that the eclipse mechanism is the result of synchrotron absorption in a relativistic pair-plasma confined to the closed-field region of pulsar B's truncated dipolar magnetic field. We demonstrate that changes in circular polarization handedness throughout the eclipses are directly tied to the average line of sight magnetic field direction of pulsar B, from which we unambiguously determine the complete magnetic and viewing geometry of the pulsar.

We consider the final evolutionary stages of a neutron star-black hole pair. According to the current paradigm, such systems eventually coalesce, which in some cases is accompanied by neutron-star tidal disruption. Using analytical methods, we show that the scenario of slow (of the order of several seconds) neutron-star stripping by the black hole is also possible, depending on the system parameters (the initial masses and intrinsic angular momenta of the components, the equation of state for the neutron star). Reaching the lower mass limit (about one tenth of the solar mass), the neutron star explodes to produce a comparatively powerful electromagnetic transient. Our population calculations show that the stripping mechanism is possible in 50-90% of the cases among all coalescing neutron star-black hole pairs, depending on the model assumptions about the evolution of close binary systems (the common-envelope efficiency parameter, the supernova explosion mechanism) and the initial metallicity of the stellar population. Because of the large mass of the ejected material, the kilonova emission in this scenario has good prospects of detection.

The search for accreted satellites in the Galactic disk is a challenging task, to which Gaia plays a crucial role in synergy with ground-based spectroscopic surveys. In 2021, Re Fiorentin et al. discovered five substructures with disk kinematics including Icarus. To gain more insight into the origin of Icarus as a remnant of a dwarf galaxy rather than a signature of secular processes of disk formation, we complement astrometric Gaia DR3 data with spectroscopy from APOGEE DR17 and GALAH DR3, and explore the chemo-dynamical distributions within 3 kpc of the Sun. We select 622 stars in the accreted/unevolved regions of [Mg/Mn]-[Al/Fe] and [Mg/Fe]-[Fe/H], where we identify 81 and 376 stars with $-2<{\rm[Fe/H]}<-0.7$ belonging to Icarus and Gaia-Sausage-Enceladus (GSE), respectively. The revised properties of Icarus are: $\langle V+V_{\rm LSR}\rangle\simeq171~\rm{km~s}^{-1}$, $\sigma_{V}\simeq37\rm{km~s}^{-1}$, $\langle e\rangle\simeq0.36$, $\langle{\rm[Fe/H]}\rangle\simeq-1.35$, $\langle{\rm[Mg/Fe]}\rangle\simeq+0.27$, $\langle{\rm[Al/Fe]}\rangle\simeq-0.13$, and $\langle{\rm[Mn/Fe]}\rangle\simeq-0.39$. From the CMD of its members, Icarus appears older than 12 Gyr. Such age and dynamical properties are reminiscent of the metal-weak thick disk. However, detailed chemical analysis in the diagnostic spaces [Ni/Fe]-[(C+N)/O], [Y/Eu]-[Fe/H], [Eu/Mg]-[Fe/H], [Ba/Y]-[Fe/H], and [Ba/Mg]-[Mg/H] evidences that Icarus and GSE occupy the accreted region, well separated from the bulk of in situ disk stars. Updated comparisons with N-body simulations confirm that Icarus' stars are consistent with the debris of a dwarf galaxy with a stellar mass of $\sim~10^9~M_\odot$ accreted onto a primordial disk on an initial prograde low-inclination orbit.

Asymptotic giant branch (AGB) stars can experience proton ingestion events (PIEs), leading to a rich nucleosynthesis. During a PIE, the intermediate neutron capture process (i-process) develops, leading to the production of trans-iron elements. It is also suggested that lithium is produced during these events. We investigate the production of lithium and trans-iron elements in AGB stars experiencing a PIE with $1<M_{\rm ini}/M_{\odot}< 3$ and $-3< \mathrm{[Fe/H]} <0$. We find that lithium is produced in all PIE models with surface abundances $3<$ A(Li) $<5$. The surface enrichment and overall AGB lithium yield increases with decreasing stellar mass. The lithium enrichment is accompanied by a production of $^{13}$C with $3<^{12}$C/$^{13}$C$<9$ at the surface just after the PIE. AGB stars experiencing PIE may be related to J-type carbon stars whose main features are excesses of lithium and $^{13}$C. In addition to Li and $^{13}$C, heavy elements (e.g., Sr, Ba, Eu, Pb) are significantly produced in low-metallicity stars up to [Fe/H]$\simeq-1$. The yields of our models are publicly available. Additionally, of interest to the Li nucleosynthesis, we provide an updated fitting formula for the $^{7}$Be($e^-,\nu_e$)$^{7}$Li electron capture rate.

A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetization $\sigma_{\rm p}$, which controls the peak energy of the neutrino spectrum, and the Eddington ratio $\lambda_{\rm Edd}$ (defined as the ratio between X-ray luminosity $L_{\rm X}$ and Eddington luminosity $L_{\rm Edd}$), which controls the amount of energy transferred to secondary particles. For sources with $\lambda_{\rm Edd} \gtrsim \lambda_{\rm Edd, crit}$ (where $\lambda_{\rm Edd, crit} \sim 10^{-1}$ for $\sigma_{\rm p}=10^5$ or $\sim 10^{-2}$ for $\sigma_{\rm p}=10^7$), proton-photon interactions and $ \gamma \gamma$ annihilation produce enough secondary pairs to achieve Thomson optical depths $\tau_{\rm T} \sim 0.1-10$. Additionally, we find that the neutrino luminosity scales as $L^2_{\rm X}/L_{\rm Edd}$ for $\lambda_{\rm Edd} \lesssim \lambda_{\rm Edd, crit}$, while it is proportional to $L_{\rm X}$ for higher $\lambda_{\rm Edd}$ values. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations.

Some exoplanets have much higher densities than expected from stellar abundances of planet-forming elements. There are two theories - metal-rich formation hypothesis and naked core hypothesis - that explain how formation and evolution can alter the compositions and structures of rocky planets to diverge from their primordial building blocks. Here, we revisit the naked core hypothesis, which states that high-density planets are remnant cores of giant planets that remain in a fossil-compressed state, even after envelope loss. Using a planetary interior model and assuming energy-limited atmospheric escape, we show that a large fraction, if not all, of the iron-silicate core of a giant planet is molten during the planet's early evolution. Upon envelope loss, molten part of the planets can rapidly rebound due to low viscosity, resulting in a decrease in radius by at most 0.06%, if they had hydrogen/helium envelopes, or by at most 7%, if they had H$_2$O envelopes, compared to self-compressed counterparts with the same core mass fraction. Based on our findings, we reject the hypothesis that all high-density exoplanets are naked cores with Kolmogorov-Smirnov p-value $\ll$ 0.05 for both envelope compositions. We find that some high-density exoplanets can still possibly be naked cores, but the probabilities are lower than $\sim$1/2 and $\sim$1/3 for the ice giant and gas giant scenario, respectively, in 95% of the cases. We conclude that most high-density exoplanets are unlikely to be remnant giant-planets cores.

Stellar active regions like spots and faculae can distort the shapes of spectral lines, inducing variations in the radial velocities that are often orders of magnitude larger than the signals from Earth-like planets. Efforts to mitigate these activity signals have hitherto focused on either the time or the velocity (wavelength) domains. We present a physics-driven Gaussian process (GP) framework to model activity signals directly in time series of line profiles or Cross-Correlation Functions (CCFs). Unlike existing methods which correct activity signals in line profile time series, our approach exploits the time correlation between velocity (wavelength) bins in the line profile variations, and is based on a simplified but physically motivated model for the origin of these variations. When tested on both synthetic and real data sets with signal-to-noise ratios down to $\sim$ 100, our method was able to separate the planetary signal from the activity signal, even when their periods were identical. We also conducted injection/recovery tests using two years of realistically sampled HARPS-N solar data, demonstrating the ability of the method to accurately recover a signal induced by a 1.5-Earth mass planet with a semi-amplitude of 0.3 m/s and a period of 33 days during high solar activity.

We present evidence for $\gamma$-ray emission from a stacked population of 39 high-latitude globular clusters (GCs) not detected in the Fermi Point Source Catalog, likely attributable to populations of millisecond pulsars within them. In this work, we use 13 years of data collected by the Large Area Telescope aboard the Fermi Gamma-Ray Space Telescope to search for a cumulative signal from undetected GCs and compared them to control fields (CFs), selected to match the celestial distribution of the target clusters so as to distinguish the $\gamma$-ray signal from background emission. The joint likelihood distribution of the GCs has a significant separation ($\sim4\sigma$) from that of the CFs. We also investigate correlations between detected cluster luminosities and other cluster properties such as distance, the number of millisecond pulsars associated with each cluster, and stellar encounter rate but find no significant relationships.

Massive main-sequence stars have convective cores and radiative envelopes, and sub-surface convection zones. However, their properties strongly depend on a star's opacity and metallicity. Non-rotating 1D evolution models of main-sequence stars between $7 \leq M \leq 40$ M$_{\odot}$ and metallicity of the SMC suggest tenuous sub-surface convection zones when using the Rayleigh number as a criterion for convection. We test if massive stars of different metallicities inside and outside of stability windows for sub-surface convection exhibit similar properties in their observed SLF variability. We extract customised light curves from the ongoing TESS mission for a sample of massive stars using an effective point spread function (ePSF) method, and compare them using a Gaussian process (GP) regression methodology. We demonstrate that the properties of SLF variability observed in time-series photometry of massive stars are consistent across the metallicity range from the Milky Way down to the SMC galaxy, for stars both inside and outside of the sub-surface stability windows. We conclude that non-rotating 1D stellar structure models of sub-surface convection cannot alone be used to explain the mechanism giving rise to SLF variability in light curves of massive stars. Additionally, the similar properties of SLF variability across a large range in metallicity, which follow the same trends in mass and age in the Hertzsprung-Russell (HR) diagram at both high and low metallicity, support a transition in the dominant mechanism causing SLF variability from younger to more evolved stars. Specifically, core-excited internal gravity waves (IGWs) are more favourable for younger stars that lack substantial sub-surface convection zones, especially at low-metallicity, and sub-surface convection zones are more favourable for the most massive and evolved stars. (abstract abridged for arxiv submission)

Vertical mixing is a crucial disequilibrium process in exoplanet atmospheres, significantly impacting chemical abundance and observed spectra. While current state-of-the-art observations have detected its signatures, the effect of vertical mixing on atmospheric spectra varies widely based on planetary parameters. In this study, we explore the influence of disequilibrium chemistry across a parameter space that includes eddy diffusion, surface gravity, internal and equilibrium temperature, and metallicity. We also assess the effectiveness of retrieval models in constraining the eddy diffusion coefficient. By running numerous 1D chemical kinetics models, we investigate the impact of vertical mixing on the transmission spectrum. We also built a custom fast-forward disequilibrium model, which includes vertical mixing using the quenching approximation and calculates the model abundance orders of magnitude faster than the chemical kinetics model. We coupled this forward model with an open source atmospheric retrieval code and used it on the JWST simulated output data of our chemical kinetics model and retrieved eddy diffusion coefficient, internal temperature and atmospheric metallicity. We find that there is a narrow region in the parameters space in which vertical mixing has a large effect on the atmospheric transmission spectrum. In this region of the parameter space, the retrieval model can put high constraints on the transport strength and provide optimal exoplanets to study vertical mixing. Also, the NH3 abundance can be used to constrain the internal temperature for equilibrium temperature T_equi > 1400 K.

We analyse the stability of the de Sitter equilibria in multi-resonant planetary systems. The de Sitter equilibrium is the dynamical state of the Laplace resonance in which all resonant arguments are librating. The sequence of equilibria exists all along the possible states balancing resonance offsets and forced eccentricities. Possible additional new-de Sitter equilibria may exist when at least one of the forced eccentricities is large (the paradigmatic case is Gliese-876). In the present work, these families of equilibria are traced up to crossing exact commensurability, where approximate first-order solutions diverge. Explicit exact location of the equilibria are determined allowing us to verify the Lyapunov stability of the standard de Sitter equilibrium and of the stable branches of the additional ones.

LMCe055-1 was recently discovered in a survey for WRs in the Large Magellanic Cloud, and classified as a WN4/O4, a lower excitation version of the WN3/O3 class discovered as part of the same survey. Its absolute magnitude precluded it from being a WN4+O4 binary. OGLE photometry show shallow primary and secondary eclipses with a 2.2 day period. The spectral characteristics and short period pointed to a possible origin due to binary stripping. Such stripped WR binaries should be common but have proven elusive to identify conclusively. In order to establish its nature, we obtained HST UV and Magellan optical spectra, along with imaging. Our work shows that the WR emission and He II absorption arise in one star, and the He I absorption in another. The He I contributor is the primary of the 2.2-day system and exhibits ~300 km/sed radial velocity variations on that time scale. However, the WR star shows 30-40 km/sed radial velocity variations, with a likely 35-day period and a highly eccentric orbit. Possibly LMCe055-1 is a physical triple, but that would require the 2.2-day pair to have been captured by the WR star. A more likely explanation is that the WR star has an unseen companion in a 35-day orbit and that the 2.2-day pair is in a longer period orbit about the two. Such examples of multiple systems are well known among massive stars, such as HD 5980. Regardless, we argue that it is highly unlikely that the WR component of the LMCe055-1 system resulted from stripping.

One of the main open questions in the field of luminous ($L_{\rm bol}>10^{47}\,\rm erg\,s^{-1}$) quasars (QSOs) at $z \gtrsim 6$ is the rapid formation ($< 1\,$Gyr) of their supermassive black holes (SMBHs). For this work we analysed the relation between the X-ray properties and other properties describing the physics and growth of both the accretion disc and the SMBH in QSOs at the Epoch of Reionization (EoR). The sample consists of 21 $z>6$ QSOs, which includes 16 sources from the rapidly grown QSOs from the HYPERION sample and five other luminous QSOs with available high-quality archival X-ray data. We discovered a strong and statistically significant ($>3\sigma$) relation between the X-ray continuum photon index ($\Gamma$) and the $\rm C\,IV$ disc wind velocity ($v_{\rm C\,IV}$) in $z>6$ luminous QSOs, whereby the higher the $v_{\rm C\,IV}$, the steeper the $\Gamma$. This relation suggests a link between the disc-corona configuration and the kinematics of disc winds. Furthermore, we find evidence at $>2-3\sigma$ level that $\Gamma$ and $v_{\rm C\,IV}$ are correlated to the growth rate history of the SMBH. Although additional data are needed to confirm it, this result may suggest that, in luminous $z>6$ QSOs, the SMBH predominantly grows via fast accretion rather than via initial high seed BH mass.