Papers for Friday, Nov 10 2023

Type Ia supernovae (SNe Ia) are standardizable cosmological candles which led to the discovery of the accelerating universe. However, the physics of how white dwarfs (WDs) explode and lead to SNe Ia is still poorly understood. The initiation of the detonation front which rapidly disrupts the WD is a crucial element of the puzzle, and global 3D simulations of SNe Ia cannot resolve the requisite length scales to capture detonation initiation. In this work, we elucidate a theoretical criterion for detonation initiation in the distributed burning regime. We test this criterion against local 3D driven turbulent hydrodynamical simulations within electron-degenerate WD matter consisting initially of pure helium. We demonstrate a novel pathway for detonation, in which strong turbulent dissipation rapidly heats the helium, and forms carbon nuclei sufficient to lead to a detonation through accelerated burning via $\alpha$ captures. Simulations of strongly-driven turbulent conditions lead to detonations at a mean density of $10^6$ g cm$^{-3}$ and mean temperature of $1.4 - 1.8 \times 10^9$ K, but fail to detonate at a lower density of $10^5$ g cm$^{-3}$, in excellent agreement with theoretical predictions.

We use high-resolution ($\simeq$ 35pc) hydrodynamical simulations of galaxy formation to investigate the relation between gas accretion and star formation in galaxies hosted by dark matter haloes of mass $10^{12}$ $\mathrm{M_\odot}$ at $z = 2$. At high redshift, cold-accreted gas is expected to be readily available for star formation, while gas accreted in a hot mode is expected to require a longer time to cool down before being able to form stars. Contrary to these expectations, we find that the majority of cold-accreted gas takes several hundred Myr longer to form stars than hot-accreted gas after it reaches the inner circumgalactic medium (CGM). Approximately 10% of the cold-accreted gas flows rapidly through the inner CGM onto the galactic disc. The remaining 90% is trapped in a turbulent accretion region that extends up to $\sim$ 50 per cent of the virial radius, from which it takes several hundred Myr for the gas to be transported to the star-forming disc. In contrast, most hot shock-heated gas avoids this 'slow track', and accretes directly from the CGM onto the disc where stars can form. We find that shock-heating of cold gas after accretion in the inner CGM and supernova-driven outflows contribute to, but do not fully explain, the delay in star formation. These processes combined slow down the delivery of cold-accreted gas to the galactic disc and consequently limit the rate of star formation in Milky Way mass galaxies at $z > 2$.

We study the role of star formation and stellar feedback in a galaxy being ram pressure stripped on its infall into a cluster. We use hydrodynamical wind-tunnel simulations of a massive galaxy ($M_\text{star} = 10^{11} M_\odot$) moving into a massive cluster ($M_\text{cluster} = 10^{15} M_\odot$). We have two types of simulations: with and without star formation and stellar feedback, SF and RC respectively. For each type we simulate four realisations of the same galaxy: a face-on wind, edge-on wind, $45^\circ$ angled wind, and a control galaxy not subject to ram pressure. We directly compare the stripping evolution of galaxies with and without star formation. We find that stellar feedback has no direct effect on the stripping process, i.e. there is no enhancement in stripping via a velocity kick to the interstellar medium gas. The main difference between RC and SF galaxies is due to the indirect effect of stellar feedback, which produces a smoother and more homogeneous interstellar medium. Hence, while the average gas surface density is comparable in both simulation types, the scatter is broader in the RC galaxies. As a result, at the galaxy outskirts overdense clumps survive in RC simulation, and the stripping proceeds more slowly. At the same time, in the inner disc, underdense gas in the RC holes is removed faster than the smoothly distributed gas in the SF simulation. For our massive galaxy, we therefore find that the effect of feedback on the stripping rate is almost negligible, independent of wind angle.

The Conditional Colour-Magnitude Distribution (CCMD) is a comprehensive formalism of the colour-magnitude-halo mass relation of galaxies. With joint modelling of a large sample of SDSS galaxies in fine bins of galaxy colour and luminosity, Xu et al. (2018) inferred parameters of a CCMD model that well reproduces colour- and luminosity-dependent abundance and clustering of present-day galaxies. In this work, we provide a test and investigation of the CCMD model by studying the colour and luminosity distribution of galaxies in galaxy groups. An apples-to-apples comparison of group galaxies is achieved by applying the same galaxy group finder to identify groups from the CCMD galaxy mocks and from the SDSS data, avoiding any systematic effect of group finding and mass assignment on the comparison. We find an overall nice agreement in the conditional luminosity function (CLF), the conditional colour function (CCF), and the CCMD of galaxies in galaxy groups inferred from CCMD mock and SDSS data. We also discuss the subtle differences revealed by the comparison. In addition, using two external catalogues constructed to only include central galaxies with halo mass measured through weak lensing, we find that their colour-magnitude distribution shows two distinct and orthogonal components, in line with the prediction of the CCMD model. Our results suggest that the CCMD model provides a good description of halo mass dependent galaxy colour and luminosity distribution. The halo and CCMD mock catalogues are made publicly available to facilitate other investigations.

Analysis of AR Sco optical light curves spanning nine years show a secular change in the relative amplitudes of the beat pulse pairs generated by the two magnetic poles of its rotating white dwarf. Recent photometry now shows that the primary and secondary beat pulses have similar amplitudes, while in 2015 the primary pulse was approximately twice that of the secondary peak. The equalization in the beat pulse amplitudes is also seen in the linearly polarized flux. This rapid evolution is consistent with precession of the white dwarf spin axis. The observations imply that the pulse amplitudes cycle over a period of $\gtrsim 40$ yrs, but that the upper limit is currently poorly constrained. If precession is the mechanism driving the evolution, then over the next 10 years the ratio of the beat pulse amplitudes will reach a maximum followed by a return to asymmetric beat pulses.

Young exoplanets trace planetary evolution, particularly the atmospheric mass loss that is most dynamic in youth. However, the high activity level of young stars can mask or mimic the spectroscopic signals of atmospheric mass loss. This includes the activity-sensitive He 10830 \AA\ triplet, which is an increasingly important exospheric probe. To characterize the He-10830 triplet at young ages, we present time-series NIR spectra for young transiting planet hosts taken with the Habitable-zone Planet Finder. The He-10830 absorption strength is similar across our sample, except at the fastest and slowest rotation, indicating that young chromospheres are dense and populate metastable helium via collisions. Photoionization and recombination by coronal radiation only dominates metastable helium population at the active and inactive extremes. Volatile stellar activity, such as flares and changing surface features, drives variability in the He-10830 triplet. Variability is largest at the youngest ages before decreasing to $\lesssim5-10$ m\AA\ (or 3%) at ages above 300 Myr, with 6 of 8 stars in this age range agreeing with no intrinsic variability. He-10830 triplet variability is smallest and age-independent at the shortest timescales. Intrinsic stellar variability should not preclude detection of young exospheres, except at the youngest ages. We recommend out-of-transit comparison observations taken directly surrounding transit and observation of multiple transits to minimize activity's effect. Regardless, caution is necessary when interpreting transit observations in the context of stellar activity, as many scenarios can lead to enhanced stellar variability even on timescales of an hour.

We report spectroscopic identification of the host galaxies of 18 ultra-strong MgII systems (USMgII) at $0.6 \leq z \leq 0.8$. We created the largest sample by merging these with 20 host galaxies from our previous survey within $0.4 \leq z \leq 0.6$. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leq D[kpc] \leq 120$ with a median of 19 kpc) are much larger than expected, and the USMgII host galaxies do not follow the canonical $\rm W_{2796}-D$ anti-correlation. We show that the presence and significance of this anti-correlation may depend on the sample selection. The $\rm W_{2796}-D$ anti-correlation seen for the general MgII absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range $0.4 \leq z \leq 1.5$ with an increase of the impact parameters. Compared to the host galaxies of normal MgII absorbers, USMgII host galaxies are brighter and more massive for a given impact parameter. While the USMgII systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of $m_r < 23.6$, at least $29\%$ of the USMgII host galaxies are isolated, and the width of the MgII absorption in these cases may originate from gas flows (infall/outflow) in isolated halos of massive star-forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21\%$ cases where interactions may cause wide velocity spread.

We are entering an era in which we will be able to detect and characterize hundreds of dwarf galaxies within the Local Volume. It is already known that a strong dichotomy exists in the gas content and star formation properties of field dwarf galaxies versus satellite dwarfs of larger galaxies. In this work, we study the more subtle differences that may be detectable in galaxies as a function of distance from a massive galaxy, such as the Milky Way. We compare smoothed particle hydrodynamic simulations of dwarf galaxies formed in a Local Volume-like environment (several Mpc away from a massive galaxy) to those formed nearer to Milky Way-mass halos. We find that the impact of environment on dwarf galaxies extends even beyond the immediate region surrounding Milky Way-mass halos. Even before being accreted as satellites, dwarf galaxies near a Milky Way-mass halo tend to have higher stellar masses for their halo mass than more isolated galaxies. Dwarf galaxies in high-density environments also tend to grow faster and form their stars earlier. We show observational predictions that demonstrate how these trends manifest in lower quenching rates, higher HI fractions, and bluer colors for more isolated dwarf galaxies.

Beyond-$\Lambda$CDM models, which were proposed to resolve the "Hubble tension", often have an impact on the discrepancy in the amplitude of matter clustering, the "$\sigma_8$-tension". To explore the interplay between the two tensions, we propose a simple method to visualize the relation between the two parameters $H_0$ and $\sigma_8$: For a given extension of the $\Lambda$CDM model and data set, we plot the relation between $H_0$ and $\sigma_8$ for different amplitudes of the beyond-$\Lambda$CDM physics. We use this visualization method to illustrate the trend of selected cosmological models, including non-minimal Higgs-like inflation, early dark energy, a varying effective electron mass, an extra number of relativistic species and modified dark energy models. We envision that the proposed method could be a useful diagnostic tool to illustrate the behaviour of complex cosmological models with many parameters in the context of the $H_0$ and $\sigma_8$ tensions.

The Kepler spacecraft, whose single instrument was a 0.95 m diameter wide-field telescope, operated in a heliocentric orbit for nearly a decade, returning a wealth of data that have revolutionized exoplanet science. Kepler data have been used to discover thousands of planets, including hundreds of multi-planet systems. Kepler discoveries have greatly expanded the diversity of known exoplanets and planetary system properties. Moreover, Kepler has provided the best estimates of exoplanet occurrence rates as functions of planetary sizes, orbital periods and stellar type, with precise values for planets with $P \lesssim 1$ yr. We provide herein an overview of the mission and its major findings regarding the occurrence rates of planets, the mass-radius relationship for exoplanets and the architectures of planetary systems.

The Milky Way (MW) dwarf spheroidal satellite galaxies (dSphs) are particularly intriguing targets to search for gamma rays from Weakly Interacting Massive Particle (WIMP) dark matter (DM) annihilation or decay. They are nearby, DM-dominated, and lack significant emission from standard astrophysical processes. Previous studies using the Fermi-Large Area Telescope (LAT) of DM emission from dSphs have provided the most robust and stringent constraints on the DM annihilation cross section and mass. We report an analysis of the MW dSphs using over 14 years of LAT data and an updated census of dSphs and $J$-factors. While no individual dSphs are significantly detected, we find slight excesses with respect to background at the $\gtrsim 2\,\sigma$ local significance level in both tested annihilation channels ($b\bar{b}$, $\tau^+\tau^-$) for 7 dSphs. We do not find a significant DM signal from a combined likelihood analysis of the dSphs ($s_{global}\sim 0.5\sigma$), yet a marginal local excess relative to background at a $2-3\,\sigma$ level is observed at a DM mass of $M_{\chi}=150-230$ GeV ($M_{\chi}=30-50$ GeV) for annihilation into $b\bar{b}$ ($\tau^+\tau^-$). Given the lack of a significant detection, we place updated constraints on the $b\bar{b}$ and $\tau^+\tau^-$ annihilation channels that are generally consistent with previous recent results. As in past studies, tension is found with some WIMP DM interpretations of the Galactic Center Excess (GCE), though the limits are consistent with other interpretations given the uncertainties of the Galactic DM density profile and GCE systematics. Based on conservative assumptions of improved sensitivity with increased LAT exposure and moderate increases in the sample of dSphs, we project the local $\sim 2\,\sigma$ signal, if real, could approach the $\sim 4\,\sigma$ local confidence level with additional $\sim 10$ years of observation.

The electromagnetic emission from the non-relativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r-process elements in the universe. However, many properties of the ejecta can alter the light-curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r-process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light-curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity and composition. Our comparisons show that uncertainties in the velocity distribution can lead to factor of 2-4 uncertainties in the inferred ejecta mass based on peak infra-red luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass.

COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Carlo simulation results of the alpha decay components of the COSINE-100 NaI(Tl) crystals. The data strongly indicate that the alpha decays probabilistically undergo one of two possible quenching factors, but this phenomenon is not understood. The fitted results are consistent with independent measurements and improve the overall understanding of the COSINE-100 backgrounds.

The Extragalactic Background Light (EBL) is the main radiation field responsible for attenuating extragalactic gamma-ray emission at very high energies, but its precise spectral intensity is not fully determined. Therefore, disentangling propagation effects from the intrinsic spectral properties of gamma-ray sources (such as active galactic nuclei, AGN) is the primary challenge to interpreting observations of these objects. We present a Bayesian and Markov Chain Monte Carlo approach to simultaneously infer parameters characterizing the EBL and the intrinsic spectra in a combined fit of a set of sources, which has the advantage of easily incorporating the uncertainties of both sets of parameters into one another through marginalization of the posterior distribution. Taking a sample of synthetic blazars observed by the ideal CTA configuration, we study the effects on the EBL constraints of combining multiple observations and varying their exposure. We also apply the methodology to a set of 65 gamma-ray spectra of 36 different AGNs measured by current Imaging Atmospheric Cherenkov Telescopes, using Hamiltonian Monte Carlo as a solution to the difficult task of sampling in spaces with a high number of parameters. We find robust constraints in the mid-IR region while simultaneously obtaining intrinsic spectral parameters for all of these objects. In particular, we identify Markarian 501 (Mkn 501) flare data (HEGRA/1997) as essential for constraining the EBL above 30$\mu$m.

We develop a model for galaxy formation and the growth of supermassive black holes (SMBHs), based on the fact that cold dark matter (CDM) halos form their gravitational potential wells through a fast phase with rapid change in the potential, and that the high universal baryon fraction makes cooled gas in halos self-gravitating and turbulent before it can form rotation-supported disks. Gas fragmentation produces sub-clouds so dense that cloud-cloud collision and drag on clouds are not significant, producing a dynamically hot system of sub-clouds that form stars and move ballistically to feed the central SMBH. Active galactic nucleus (AGN) and supernova (SN) feedback is effective only in the fast phase, and the cumulative effects are to regulate star formation and SMBH growth, as well as to reduce the amount of cold gas in halos to allow the formation of globally stable disks. Using a set of halo assembly histories, we demonstrate that the model can reproduce a number of observations, including correlations among SMBH mass, stellar mass of galaxies and halo mass, the number densities of galaxies and SMBH, as well as their evolution over the cosmic time.

We report spectroscopic observations of vB 120 (HD 30712), a 5.7 yr astrometric-spectroscopic binary system in the Hyades cluster. We combine our radial velocities with others from the literature, and with existing speckle interferometry measurements, to derive an improved 3D orbit for the system. We infer component masses of M1 = 1.065 +/- 0.018 MSun and M2 = 1.008 +/- 0.016 MSun, and an orbital parallax of 21.86 +/- 0.15 mas, which we show to be more accurate than the parallax from Gaia DR3. This is the ninth binary or multiple system in the Hyades with dynamical mass determinations, and one of the examples with the highest precision. An analysis of the spectral energy distribution yields the absolute radii of the stars, R1 = 0.968 +/- 0.012 RSun and R2 = 0.878 +/- 0.013 RSun, and effective temperatures of 5656 +/- 56 K and 5489 +/- 60 K for the primary and secondary, respectively. A comparison of these properties with the predictions of current stellar evolution models for the known age and metallicity of the cluster shows only minor differences.

In this paper, we study six slowly rotating mid-to-late M~dwarfs (rotation period $P_{\mathrm{rot}} \approx 40-190\,\mathrm{dy}$) by analysing spectropolarimetric data collected with SPIRou at the Canada-France-Hawaii Telescope as part of the SPIRou Legacy Survey from 2019 to 2022. From $\approx$100--200 Least-Squares-Deconvolved (LSD) profiles of circularly polarised spectra of each star, we confirm the stellar rotation periods of the six M~dwarfs and explore their large-scale magnetic field topology and its evolution with time using both the method based on Principal Component Analysis (PCA) proposed recently and Zeeman-Doppler Imaging. All M~dwarfs show large-scale field variations on the time-scale of their rotation periods, directly seen from the circularly polarised LSD profiles using the PCA method. We detect a magnetic polarity reversal for the fully-convective M~dwarf GJ~1151, and a possible inversion in progress for Gl~905. The four fully-convective M~dwarfs of our small sample (Gl~905, GJ~1289, GJ~1151, GJ~1286) show a larger amount of temporal variations (mainly in field strength and axisymmetry) than the two partly-convective ones (Gl~617B, Gl~408). Surprisingly, the six M~dwarfs show large-scale field strengths in the range between 20 to 200\,G similar to those of M~dwarfs rotating significantly faster. Our findings imply that the large-scale fields of very slowly rotating M~dwarfs are likely generated through dynamo processes operating in a different regime than those of the faster rotators that have been magnetically characterized so far.

Given the recent successful launch of the James Webb Space Telescope, determining robust calibrations of the slopes and absolute magnitudes of the near- to mid-infrared Tip of the Red Giant Branch (TRGB) will be essential to measuring precise extragalactic distances via this method. Using ground-based data of the Large Magellanic Cloud from the Magellanic Clouds Photometric Survey along with near-infrared (NIR) data from 2MASS and mid-infrared (MIR) data collected as a part of the SAGE survey using the Spitzer Space Telescope, we present slopes and zero-points for the TRGB in the optical (VI), NIR (JHK) and MIR ([3.6] & [4.5]) bandpasses. These calibrations utilize stars +0.3 +/- 0.1 mag below the tip, providing a substantial statistical improvement over previous calibrations which only used the sample of stars narrowly encompassing the tip.

We propose an approach to infer large-scale heterogeneities within a small celestial body from measurements of its gravitational potential, provided for instance by spacecraft radio-tracking. The non-uniqueness of the gravity inversion is here mitigated by limiting the solutions to piecewise-constant density distributions, thus composed of multiple regions of uniform density (mass anomalies) dispersed in a background medium. The boundary of each anomaly is defined implicitly as the 0-level surface of a scalar field (called the level-set function), so that by modifying this field the shape and location of the anomaly are varied. The gravitational potential associated with a density distribution is here computed via a line-integral polyhedron method, yielding the coefficients of its spherical harmonics expansion. The density distribution is then adjusted via an iterative least-squares approach with Tikhonov regularization, estimating at every iteration corrections to the level-set function, the density contrast of each anomaly, and the background density, in order to minimize the residuals between the predicted gravity coefficients and those measured. Given the non-convexity of the problem and the lack of prior knowledge assumed (save for the shape of the body), the estimation process is repeated for several random initial distributions, and the resulting solutions are clustered based on global properties independent of the input measurements. This provides families of candidate interior models in agreement with the data, and the spread of the local density values across each family is used to assess the uncertainties associated with the estimated distributions.

Bright quasar samples at high redshift are useful for investigating active galactic nuclei evolution. In this study, we describe XQz5, a sample of 83 ultraluminous quasars in the redshift range $4.5 < z < 5.3$ with optical and near-infrared spectroscopic observations, with unprecendented completeness at the bright end of the quasar luminosity function. The sample is observed with the Southern Astrophysical Research Telescope, the Very Large Telescope, and the ANU 2.3m Telescope, resulting in a high-quality, moderate-resolution spectral atlas of the brightest known quasars within the redshift range. We use established virial mass relations to derive the black hole masses by measuring the observed Mg\,\textsc{ii}$\lambda$2799\AA\ emission-line and we estimate the bolometric luminosity with bolometric corrections to the UV continuum. Comparisons to literature samples show that XQz5 bridges the redshift gap between other X-shooter quasar samples, XQ-100 and XQR-30, and is a brighter sample than both. Luminosity-matched lower-redshift samples host more massive black holes, which indicate that quasars at high redshift are more active than their counterparts at lower-redshift, in concordance with recent literature.

We analyze the energy density spectrum of \acp{SIGW} using the NANOGrav 15-year data set, thereby constraining the primordial non-Gaussian parameter $f_{\mathrm{NL}}$. For the first time, we calculate the seventeen missing two-loop diagrams proportional to $f_{\mathrm{NL}}A_{\zeta}^3$ that correspond to the two-point correlation function $\langle h^{\lambda,(3)}_{\mathbf{k}} h^{\lambda',(2)}_{\mathbf{k}'} \rangle$ for local-type primordial non-Gaussianity. The total energy density spectrum of \acp{SIGW} can be significantly suppressed by these two-loop diagrams. If \acp{SIGW} dominate the \acp{SGWB} observed in \ac{PTA} experiments, the parameter interval $f_{\mathrm{NL}}\in [-5,-1]$ is notably excluded based on NANOGrav 15-year data set. After taking into account abundance of \acp{PBH} and the convergence of the cosmological perturbation expansion, we find that the only possible parameter range for $f_{\mathrm{NL}}$ might be $-1\le f_{\mathrm{NL}}< 0$.

Extreme precision radial velocity (EPRV) measurements contend with internal noise (instrumental systematics) and external noise (intrinsic stellar variability) on the road to 10 cm/s "exo-Earth" sensitivity. Both of these noise sources are well-probed using "Sun-as-a-star" RVs and cross-instrument comparisons. We built the Solar Calibrator (SoCal), an autonomous system that feeds stable, disc-integrated sunlight to the recently commissioned Keck Planet Finder (KPF) at the W. M. Keck Observatory. With SoCal, KPF acquires signal-to-noise ~1200, R = ~98,000 optical (445--870 nm) spectra of the Sun in 5~sec exposures at unprecedented cadence for an EPRV facility using KPF's fast readout mode (<16 sec between exposures). Daily autonomous operation is achieved by defining an operations loop using state machine logic. Data affected by clouds are automatically flagged using a reliable quality control metric derived from simultaneous irradiance measurements. Comparing solar data across the growing global network of EPRV spectrographs with solar feeds will allow EPRV teams to disentangle internal and external noise sources and benchmark spectrograph performance. To facilitate this, all SoCal data products are immediately available to the public on the Keck Observatory Archive. We compared SoCal RVs to contemporaneous RVs from NEID, the only other immediately public EPRV solar dataset. We find agreement at the 30-40 cm/s level on timescales of several hours, which is comparable to the combined photon-limited precision. Data from SoCal were also used to assess a detector problem and wavelength calibration inaccuracies associated with KPF during early operations. Long-term SoCal operations will collect upwards of 1,000 solar spectra per six-hour day using KPF's fast readout mode, enabling stellar activity studies at high signal-to-noise on our nearest solar-type star.

In this manuscript, strong evidence is reported to support unobscured broad line regions (BLRs) in Type-1.9 AGN SDSS J1241+2602 with reliable broad H$\alpha$ but no broad H$\beta$. Commonly, disappearance of broad H$\beta$ can be explained by the AGN unified model expected heavily obscured BLRs in Type-1.9 AGN. Here, based on properties of two kinds of BH masses, the virial BH mass and the BH mass through the \msig relation, an independent method is proposed to test whether are there unobscured central BLRs in a Type-1.9 AGN. By the reliable measurement of stellar velocity dispersion about 110$\pm$12km/s through the host galaxy absorption features in \obj, the BH mass through the \msig relation is consistent with the virial BH mass $(3.43\pm1.25)\times10^7{\rm M_\odot}$ determined through properties of the observed broad H$\alpha$ without considering effects of obscurations in SDSS J1241+2602. Meanwhile, if considering heavily obscured BLRs in SDSS J1241+2602, the reddening corrected virial BH mass is tens of times larger than the \msig expected value, leading SDSS J1241+2602 to be an outlier in the \msig space with confidence level higher than $5\sigma$. Therefore, the unobscured BLRs are preferred in the Type-1.9 AGN SDSS J1241+2602. The results indicate that it is necessary to check whether unobscured central BLRs are common in Type-1.9 AGN, when to test the AGN unified model of AGN by properties of Type-1.9 AGN.

A fundamental goal of modern-day astrophysics is to understand the connection between supermassive black hole (SMBH) growth and galaxy evolution. Merging galaxies offer one of the most dramatic channels for galaxy evolution known, capable of driving inflows of gas into galactic nuclei, potentially fueling both star formation and central SMBH activity. Dual active galactic nuclei (dual AGNs) in late-stage mergers with nuclear pair separations $<10$ kpc are thus ideal candidates to study SMBH growth along the merger sequence since they coincide with the most transformative period for galaxies. However, dual AGNs can be extremely difficult to confirm and study. Hard X-ray ($>10$ keV) studies offer a relatively contamination-free tool for probing the dense obscuring environments predicted to surround the majority of dual AGN in late-stage mergers. To date, only a handful of the brightest and closest systems have been studied at these energies due to the demanding instrumental requirements involved. We demonstrate the unique capabilities of HEX-P to spatially resolve the soft and - for the first time - hard X-ray counterparts of closely-separated ($\sim2''-5''$) dual AGNs in the local Universe. By incorporating state-of-the-art physical torus models, we reproduce realistic broadband X-ray spectra expected for deeply embedded accreting SMBHs. Hard X-ray spatially resolved observations of dual AGNs - accessible only to HEX-P - will hence transform our understanding of dual AGN in the nearby Universe.

The CGM hosts many physical processes with different kinematic signatures that affect galaxy evolution. We address the CGM-galaxy kinematic connection by quantifying the fraction of HI that is aligned with galaxy rotation with the equivalent width co-rotation fraction, $f_{\rm EWcorot}$. Using 70 quasar sightlines having HST/COS HI absorption (${12<\log (N(HI)/{\rm cm}^{-2})<20}$) within $5R_{\rm vir}$ of $z<0.6$ galaxies we find that $f_{\rm EWcorot}$ increases with increasing HI column density. $f_{\rm EWcorot}$ is flat at $\sim0.6$ within $R_{\rm vir}$ and decreases beyond $R_{\rm vir}$ to $f_{\rm EWcorot}$$\sim0.35$. $f_{\rm EWcorot}$ also has a flat distribution with azimuthal and inclination angles within $R_{\rm vir}$, but decreases by a factor of two outside of $R_{\rm vir}$ for minor axis gas and by a factor of two for edge-on galaxies. Inside $R_{\rm vir}$, co-rotation dominated HI is located within $\sim 20$ deg of the major and minor axes. We surprisingly find equal amounts of HI absorption consistent with co-rotation along both major and minor axes within $R_{\rm vir}$. However, this co-rotation disappears along the minor axis beyond $R_{\rm vir}$, suggesting that if this gas is from outflows, then it is bound to galaxies. $f_{\rm EWcorot}$ is constant over two decades of halo mass, with no decrease for log(M$_{\rm h}/M_{\odot})>12$ as expected from simulations. Our results suggest that co-rotating gas flows are best found by searching for higher column density gas within $R_{\rm vir}$ and near the major and minor axes.

The multiphase CGM hosts critical processes that affect galaxy evolution such as accretion and outflows. We searched for evidence of these phenomena by using the EW co-rotation fraction ($f_{\rm EWcorot}$) to study the kinematic connection between the multiphase CGM and host galaxy rotation. We examined 27 systems with absorption lines from HST/COS (including, but not limited to, SiII, CII, SiIII, CIII, and OVI) within $21\leq D\leq~276$ kpc of galaxies. We find the median $f_{\rm EWcorot}$ for all ions is consistent within errors and the $f_{\rm EWcorot}$ increases with increasing N$(\rm HI)$. The $f_{\rm EWcorot}$ of lower ionization gas likely decreases with increasing $D/R_{\rm vir}$ while OVI and HI are more consistent with being flat. The $f_{\rm EWcorot}$ varies minimally as a function of azimuthal angle and is similar for all ions at a fixed azimuthal angle. The larger number of OVI detections enabled us to investigate where the majority of co-rotating gas is found. Highly co-rotating OVI primarily resides along the galaxies' major axis. Looking at the $f_{\rm EWcorot}$ as a function of ionization potential (${d{f_{\rm EWcorot}}}/{d{(eV)}}$), we find a stronger co-rotation signature for lower-ionization gas. There are suggestions of a connection between the CGM metallicity and major axis co-rotation where low-ionization gas with higher $f_{\rm EWcorot}$ exhibits lower metallicity and may trace large-scale filamentary inflows. Higher ionization gas with higher $f_{\rm EWcorot}$ exhibits higher metallicity and may instead trace co-planar recycled gas accretion. Our results stress the importance of comparing absorption originating from a range of ionization phases to differentiate between various gas flow scenarios.

Theoretical studies of angular momentum transport suggest that isolated stellar-mass black holes are born with negligible dimensionless spin magnitudes $\chi \lesssim 0.01$. However, recent gravitational-wave observations indicate $\gtrsim 15\%$ of binary black hole systems contain at least one black hole with a non-negligible spin magnitude. One explanation is that the first-born black hole spins up the stellar core of what will become the second-born black hole through tidal interactions. Typically, the second-born black hole is the ``secondary'' (less-massive) black hole, though, it may become the ``primary'' (more-massive) black hole through a process known as mass-ratio reversal. We investigate this hypothesis by analysing data from the third gravitational-wave transient catalog (GWTC-3) using a ``single-spin'' framework in which only one black hole may spin in any given binary. Given this assumption, we show that at least $28\%$ (90% credibility) of the LIGO--Virgo--KAGRA binaries contain a primary with significant spin, possibly indicative of mass-ratio reversal. We find no evidence for binaries that contain a secondary with significant spin. However, the single-spin framework is moderately disfavoured (natural log Bayes factor $\ln {\cal B} = 3.1$) when compared to a model that allows both black holes to spin. If future studies can firmly establish that most merging binaries contain two spinning black holes, it may call into question our understanding of formation mechanisms for binary black holes or the efficiency of angular momentum transport in black hole progenitors.

We present a polarization analysis of PSR J0941$-$39 and PSR J1107$-$5907, which exhibit transitions between being pulsars and rotating radio transients (RRATs), using the ultra-wide bandwidth low-frequency (UWL) receiver on Murriyang, the Parkes 64\,m radio telescope. The spectral index of each pulsar was measured, revealing distinct variations among different states. By using the rotating vector model (RVM), we determined that the magnetosphere geometry remains consistent between the RRAT state and the pulsar state for PSR J0941$-$39, with emissions originating from the same height in the magnetosphere. The occurrence of the RRAT state could be attributed to variations in currents within the pulsar's magnetosphere. Our results suggest that the emission mechanism of RRAT may share similarities with that of a typical pulsar.

The full disk, full of gas and dust, determines the upper limit of planet masses, and its lifetime is critical for planet formation, especially for giant planets. In this work, we studied the evolutionary timescales of the full disks of T Tauri stars (TTSs) and their relations to accretion. Combined with Gaia EDR3, 2MASS, and WISE data, 1077 disk-bearing TTS candidates were found in LAMOST DR8, and stellar parameters were obtained. Among them, 783 are newly classified by spectra as classical T Tauri stars (CTTSs; 169) or weak-lined T Tauri stars (WTTSs). Based on EW and FWHM of Ha, 157 TTSs in accretion were identified, with ~ 82% also having full disks. For TTSs with M<0.35 Mo, about 80% seem to already lose their full disks at ~ 0.1 Myr, which may explain their lower mass, while the remaining 20% with full disks evolve at similar rates of non-full disks within 5 Myr, possibly suffice to form giant planets. The fraction of accreting TTSs to disk-bearing TTSs is stable at ~10% and can last $\sim$ 5-10 Myr, suggesting that full disks and accretion evolve with similar rates as non-full disks. For TTSs with M>0.35Mo, almost all full disks can survive more than 0.1 Myr, most for 1 Myr and some even for 20 Myr, which implies planets are more likely to be formed in their disks than those of M<0.35 Mo, and thus M dwarfs with M>0.35Mo can have more planets. The fraction of full-disk TTSs to disk-bearing TTSs decreases with age following the relation $f\propto t^{-0.35}$, and similar relations existed in the fraction of accreting TTSs and the fraction of full-disk CTTSs, suggesting faster full disks and accretion evolution than non-full disks. For full disk stars, the ratio of accretion of lower-mass stars is systematically lower than that of higher-mass stars, confirming the dependence of accretion on stellar mass.

High-resolution (HR) simulations in cosmology, in particular when including baryons, can take millions of CPU hours. On the other hand, low-resolution (LR) dark matter simulations of the same cosmological volume use minimal computing resources. We develop a denoising diffusion super-resolution emulator for large cosmological simulation volumes. Our approach is based on the image-to-image Palette diffusion model, which we modify to 3 dimensions. Our super-resolution emulator is trained to perform outpainting, and can thus upgrade very large cosmological volumes from LR to HR using an iterative outpainting procedure. As an application, we generate a simulation box with 8 times the volume of the Illustris TNG300 training data, constructed with over 9000 outpainting iterations, and quantify its accuracy using various summary statistics.

Imaging spectroscopy is intended to be coupled with adaptive optics (AO) on large telescopes, such as EST, in order to produce high spatial and temporal resolution measurements of velocities and magnetic fields upon a 2D FOV. We propose a high spectral resolution slicer (30 m{\AA} typical) for the Multichannel Subtractive Double Pass (MSDP) of the future European Solar Telescope (EST), using a new generation slicer for thin photospheric lines such as FeI (56 channels, 0.13 mm step) which will benefit of AO and existing polarimeters. The aim is to reconstitute cubes of instantaneous data (X, Y, lambda) at high cadence, allowing the study of the photospheric dynamics and magnetic fields.

Knowledge of the spectrograph's instrumental profile (IP) provides important information needed for wavelength calibration and for the use in scientific analyses. This work develops new methods for IP reconstruction in high resolution spectrographs equipped with Laser Frequency Comb calibration (LFC) systems and assesses the impact that assumptions on IP shape have on achieving accurate spectroscopic measurements. Astronomical LFCs produce $\approx10000$ bright, unresolved emission lines with known wavelengths, making them excellent probes of the IP. New methods based on Gaussian Process regression were developed to extract detailed information on the IP shape from this data. Applying them to HARPS, an extremely stable spectrograph installed on the ESO 3.6m telescope, we reconstructed its IP at 512 locations of the detector, covering 60% of the total detector area. We found that the HARPS IP is asymmetric and that it varies smoothly across the detector. Empirical IP models provide wavelength accuracy better than 10 ms$^{-1}$ (5 ms$^{-1}$) with 92% (64%) probability. In comparison, reaching the same accuracy has a probability of only 29% (8%) when a Gaussian IP shape is assumed. Furthermore, the Gaussian assumption is associated with intra-order and inter-order distortions in the HARPS wavelength scale as large as 60ms$^{-1}$. The spatial distribution of these distortions suggests they may be related to spectrograph optics and therefore may generally appear in cross-dispersed echelle spectrographs when Gaussian IPs are used. Methods presented here can be applied to other instruments equipped with LFCs, such as ESPRESSO, but also ANDES and G-CLEF in the future. The empirical IPs will be crucial for obtaining objective and unbiased measurements of fundamental constants from high resolution spectra, as well as measurements of the redshift drift, isotopic abundances, and other science cases.

The theory of Nested Figures of Equilibrium, expanded in Papers I and II, is investigated in the limit where the number of layers of the rotating body is infinite, enabling to reach full heterogeneity. In the asymptotic process, the discrete set of equations becomes a differential equation for the rotation rate. In the special case of rigid rotation (from center to surface), we are led to an Integro-Differential Equation (IDE) linking the ellipticity of isopycnic surfaces to the equatorial mass-density profile. In constrast with most studies, these equations are not restricted to small flattenings, but are valid for fast rotators as well. We use numerical solutions obtained from the SCF-method to validate this approach. At small ellipticities (slow rotation), we fully recover Clairaut's equation. Comparisons with Chandrasekhar's perturbative approach and with Roberts' work based on Virial equations are successful. We derive a criterion to characterize the transition from slow to fast rotators. The treatment of heterogeneous structures containing mass-density jumps is proposed through a modified IDE.

We extend the study of Papaloizou & Savonije of the tidal interactions of close orbiting giant planets with a central solar type star to the situation where the spin axis of the central star and the orbital angular momentum are misaligned. We determine the tidal response taking into account the possibility of the excitation of r modes and the effect of tidal forcing due to potential perturbations which have zero frequency in a non rotating frame. Although there is near resonance with r modes with degree l' = 1 and orders m = 1 or -1 , half widths turn out to be sufficiently narrow so that in practice dissipation rates are found to be similar to those produced by non resonant potential perturbations. We use our results to determine the evolution of the misalignment for the full range of initial inclination angles taking account of the spin down of the central star due to magnetic braking. Overall we find the rate of tidal evolution to be unimportant for a one Jupiter mass planet with orbital period of about 3.7d over a main sequence lifetime. However, it becomes significant for higher mass planets and shorter orbital periods, approximately scaling as the square of the planet mass and the inverse fourth power of the orbital period.

We analyse perturbations of self-interacting, scalar field dark matter that contains modes both in a coherent condensate state and an incoherent particle-like state. Starting from the coupled equations for the condensate, the particles and their mutual gravitational potential, first derived from first principles in earlier work by the authors, we derive a hydrodynamic limit of two coupled fluids and study their linearized density perturbations in an expanding universe. We find that away from the condensate-only or particle-only limits, and for certain ranges of the parameters, such self-interacting mixtures can significantly enhance the density power spectrum above the standard linear $\Lambda$CDM value at localised wavenumbers, even pushing structure formation into the non-linear regime earlier than expected in $\Lambda$CDM for these scales. We also note that such mixtures can lead to degeneracies between models with different boson masses and self-coupling strengths, in particular between self-coupled models and non-coupled Fuzzy Dark Matter made up of heavier bosons. These findings open up the possibility of a richer phenomenology in scalar field dark matter models and could further inform efforts to place observational limits on their parameters.

Neutron stars and black holes in X-ray binaries are observed to host strong collimated jets in the hard spectral state. Numerical simulations can act as a valuable tool in understanding the mechanisms behind jet formation and its properties. Although there have been significant efforts in understanding black-hole jets from general-relativistic magnetohydrodynamic (GRMHD) simulations in the past years, neutron star jets, however, still remain poorly explored. We present the results from three-dimensional (3D) GRMHD simulations of accreting neutron stars with oblique magnetospheres for the very first time. The jets in our simulations are produced due to the anchored magnetic field of the rotating star in analogy with the Blandford-Znajek process. We find that for accreting stars, the star-disk magnetic field interaction plays a significant role and as a result, the jet power becomes directly proportional to ${\Phi^2}_{\rm jet}$, where $\Phi_{\rm jet}$ is the open flux in the jet. The jet power decreases with increasing stellar magnetic inclination and finally for an orthogonal magnetosphere, it reduces by a factor of $\simeq 2.95$ compared to the aligned case. We also find that in the strong propeller regime, with a highly oblique magnetosphere, the disk-induced collimation of the open stellar flux preserves parts of the striped wind resulting in a striped jet.

The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array (ALMA). The dust continuum emission revealed a filamentary structure with a length of $\sim$5000 au and an average H$_2$ volume density of $\sim$6 $\times$ 10$^7$ cm$^{-3}$. At the tip of this filamentary structure, there is a compact object, which we call a ``nucleus", with a radius of $\sim$150--200 au and a mass of $\sim$0.1 M$_{\odot}$. The nucleus has a central density of $\sim$2 $\times$ 10$^9$ cm$^{-3}$ with a radial density profile of $r^{-1.87{\pm}0.11}$. The density scaling of the nucleus is $\sim$3.7 times higher than that of the singular isothermal sphere. This as well as the very low virial parameter of 0.39 suggest that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N$_2$D$^+$ 3--2 emission, but not by the C$^{18}$O 2--1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N$_2$D$^+$ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N$_2$, or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC.

We explore the cosmology of the Dark Dimension scenario taking into account perturbations in the linear regime. In the context of the Dark Dimension scenario, a natural candidate for dark matter in our universe is the excitations of a tower of massive spin-2 KK gravitons. These dark gravitons are produced in the early universe and decay to lighter KK gravitons during the course of cosmological evolution. The decay causes the average dark matter mass to decrease as the universe evolves. In addition, the kinetic energy liberated in each decay leads to a kick velocity for the dark matter particles leading to a suppression of structure formation. Using current CMB (Planck), BAO and cosmic shear (KiDS-1000) data, we put a bound on the dark matter kick velocity today $v_\mathrm{today} \leq 2.2 \times 10^{-4} c$ at 95\% CL. This leads to rather specific regions of parameter space for the dark dimension scenario. The combination of the experimental bounds from cosmology, astrophysics and table-top experiments lead to the range $l_5\sim 1- 10 \, \mu m$ for the size of the Dark Dimension. The Dark Dimension scenario is found to be remarkably consistent with current observations and provides signatures that are within reach of near-future experiments.

We propose a model-independent \textit{B\'ezier parametric interpolation} to alleviate the degeneracy between baryonic and dark matter abundances by means of intermediate-redshift data. To do so, we first interpolate the observational Hubble data to extract cosmic bounds over the (reduced) Hubble constant, $h_0$, and interpolate the angular diameter distances, $D(z)$, of the galaxy clusters, inferred from the Sunyaev-Zeldovich effect, constraining the spatial curvature, $\Omega_k$. Through the so-determined Hubble points and $D(z)$, we interpolate uncorrelated data of baryonic acoustic oscillations bounding the baryon ($\omega_b = h^2_0\Omega_b$) and total matter ($\omega_m = h^2_0\Omega_m$) densities, reinforcing the constraints on $h_0$ and $\Omega_k$ with the same technique. Instead of pursuing the usual treatment to fix $\omega_b$ via the value obtained from the cosmic microwave background to remove the matter sector degeneracy, we here interpolate the acoustic parameter from correlated baryonic acoustic oscillations. The results of our Monte Carlo--Markov chain simulations turn out to agree at $1$--$\sigma$ confidence level with the flat $\Lambda$CDM model. While our findings are roughly suitable at $1$--$\sigma$ with its non-flat extension too, the Hubble constant appears in tension up to the $2$--$\sigma$ confidence level. Accordingly, we also reanalyze the Hubble tension with our treatment and find our expectations slightly match local constraints.

Astrophysical radio signals are excellent probes of extreme physical processes that emit them. However, to reach Earth, electromagnetic radiation passes through the ionised interstellar medium (ISM), introducing a frequency-dependent time delay (dispersion) to the emitted signal. Removing dispersion enables searches for transient signals like Fast Radio Bursts (FRB) or repeating signals from isolated pulsars or those in orbit around other compact objects. The sheer volume and high resolution of data that next generation radio telescopes will produce require High-Performance Computing (HPC) solutions and algorithms to be used in time-domain data processing pipelines to extract scientifically valuable results in real-time. This paper presents a state-of-the-art implementation of brute force incoherent dedispersion on NVIDIA GPUs, and on Intel and AMD CPUs. We show that our implementation is 4x faster (8-bit 8192 channels input) than other available solutions and demonstrate, using 11 existing telescopes, that our implementation is at least 20 faster than real-time. This work is part of the AstroAccelerate package.

Isotope ratios have recently been measured in the atmospheres of directly-imaged and transiting exoplanets from ground-based observations. The arrival of JWST allows us to characterise exoplanetary atmospheres in further detail and opens up wavelengths inaccessible from the ground. In this work we constrain the carbon and oxygen isotopes $^{13}$C, $^{18}$O and $^{17}$O from CO in the atmosphere of the directly-imaged companion VHS 1256 b through retrievals of the $\sim$4.1-5.3 $\mu$m NIRSpec G395H/F290LP observations from the early release science programme (ERS 1386). We detect and constrain $^{13}$C$^{16}$O, $^{12}$C$^{18}$O and $^{12}$C$^{17}$O at 32, 16 and 10$\sigma$ confidence respectively, thanks to the very high signal-to-noise observations. We find the ratio of abundances are more precisely constrained than their absolute values, with $\mathrm{^{12}C/^{13}C=62^{+2}_{-2}}$, in between previous measurements for companions ($\sim$30) and isolated brown dwarfs ($\sim$100). The oxygen isotope ratios are $\mathrm{^{16}O/^{18}O =425^{+33}_{-28}}$ and $\mathrm{^{16}O/^{17}O=1010^{+120}_{-100}}$. All of the ratios are lower than the local inter-stellar medium and Solar System, suggesting that abundances of the more minor isotopes are enhanced compared to the primary. This could be driven by isotope fractionation in protoplanetary disks, which can potentially alter the carbon and oxygen ratios through isotope selective photodissociation, gas/ice partitioning and isotopic exchange reactions. In addition to CO, we constrain $^{1}$H$_2$$^{16}$O and $^{12}$C$^{16}$O$_2$ (the primary isotopologues of both species), but find only upper limits on $^{12}$C$^1$H$_4$ and $^{14}$N$^{1}$H$_3$. This work highlights the power of JWST to constrain isotopes in exoplanet atmospheres, with great promise in determining formation histories in the future.

The redshifted 21-cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study, we use a single night of observations with the New Extension in Nan\c{c}ay Upgrading LOFAR (NenuFAR) to place upper limits on the 21-cm power spectrum from the Cosmic Dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nan\c{c}ay Radio Observatory in France. It is a phased array instrument with a very dense uv-coverage at short baselines, making it one of the most sensitive instruments for 21-cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration, and residual foregrounds are subsequently removed using Gaussian process regression (GPR). The final power spectra are constructed from the gridded residual data cubes in the uv-plane. Signal injection tests are performed at each step of the analysis pipeline, and the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$\sigma$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements in the analysis pipeline, which can potentially mitigate it and consequently allow us to reach the thermal noise sensitivity when multiple nights of observations are processed in the future.

The hydrogen 21-cm signal is predicted to be the richest probe of the young Universe including eras known as the cosmic Dark Ages, the Cosmic Dawn when the first star and black hole formed, and the Epoch of Reionization. This signal holds the key to deciphering processes that take place at the early stages of cosmic history. In this opinion piece, we discuss the potential scientific merit of lunar observations of the 21-cm signal and their advantages over more affordable terrestrial efforts. The moon is a prime location for radio cosmology which will enable precision observations of the low-frequency radio sky. The uniqueness of such observations is that they will provide an unparalleled opportunity to test cosmology and the nature of dark matter using the Dark Ages 21-cm signal. No less enticing is the opportunity to obtain a much clearer picture of Cosmic Dawn than what is achievable from the ground, which will allow us to probe properties of the first stars and black holes.

Studies of the pulsar B0823+26 have been carried out using the Large Phased Array (LPA) radio telescope. At time span of 5.5 years, the amplitudes of the main pulse (MP), postcursor (PC) and interpulse (IP) were evaluated in daily sessions lasting 3.7 minutes. It is shown that the ratio of the average amplitudes of MP in the bright (B) and quiet (Q) modes is 60. For B-mode, the average ratio of MP amplitudes to IP amplitudes is 65, and the ratio of MP amplitudes to PC amplitudes is 28. The number of sessions with a nulling is 4% of the total number of sessions. Structure function (SF) and correlation function analysis of MP, IP and PC amplitude variations of over a long-time interval allowed us to detect typical time scales 37 \pm 5 days and one year. The analysis of time variations shows that the time scale of 37 days is well explained by refraction on inhomogeneities of interstellar plasma, which is distributed mostly quasi-uniformly in the line-of-sight. This scintillation makes the main contribution to the observed variability. Analysis of the structure function showed that there may be a few days variability. This time scale does not have an unambiguous interpretation but is apparently associated with the refraction of radio waves on the interstellar medium. One-year variability time scale has not been previously detected. We associate its appearance with the presence of a scattering layer on a closely located screen at a distance of about 50-100 pc from the Earth.

After 25 years of the prediction of the possibility of observations, and despite the many hundred of well studied transiting exoplanet systems, we are still waiting for the announcement of the first confirmed exomoon. Following the ``cascade'' structure of the Drake-equation, but applied to the existence of an observable exomoon instead the original scope of extraterrestrial intelligence. The scope of this paper is to reveal the structure of the problem, rather than giving a quantitative solution. We identify three important steps that can lead us to the discovery. The steps are the formation, the orbital dynamics and long-term stability, and the observability of a given exomoon in a given system. This way, the question will be closely related to questions of star formation, planet formation, 5 possible pathways of moon formation; long-term dynamics of evolved planet systems involving stellar and planetary rotation and internal structure; and the proper evaluation of the observed data, taken the correlated noise of stellar and instrumental origin and the sampling function also into account. This way, a successful exomoon observation and the interpretations of the expectable further measurements proves to be among the most complex and most interdisciplinary questions in astrophysics.

Theoretical predictions suggest that very massive stars have the potential to form through multiple collisions and eventually evolve into intermediate-mass black holes (IMBHs) within Population III star clusters that are embedded in mini dark matter haloes. In this study, we investigate the long-term evolution of Population III star clusters, including models with a primordial binary fraction of $f_{\rm b}=0$ and 1, using the $N$-body simulation code \textsc{petar}. We comprehensively examine the phenomenon of hierarchical triple black holes in the clusters, specifically focusing on the merging of their inner binary black holes (BBHs), with post-Newtonian correction, by using the \textsc{tsunami} code. Our findings suggest a high likelihood of the inner BBHs containing IMBHs with masses on the order of $\mathcal{O}(100)M_{\odot}$, and as a result, their merger rate could be up to $0.1{\rm Gpc}^{-3}{\rm yr}^{-3}$. In the model with $f_{\rm b}=0$, the evolution of these merging inner BBHs is predominantly driven by their gravitational wave radiation at an early time, but their evolutionary dynamics are dominated by the interaction between them and tertiary BHs in the case with $f_{\rm b}=1$. The orbital eccentricities of some merging inner BBHs oscillate over time periodically, known as the Kozai-Lidov oscillation, due to dynamical perturbations. Furthermore, the merging inner BBHs tend to have highly eccentric orbits at low frequency range, some of them with low redshift would be detected by LISA/TianQin. In the higher frequency range, the merging inner BBHs across a wider redshift range would be detected by DECIGO/ET/CE/LIGO/KAGRA.

In this letter, we measure the rest-frame optical and near-infrared sizes of ten quiescent candidates at 3<z<5, first reported by Carnell et al. (2023a). We use $\textit{James Webb Space Telescope}$ (JWST) Near-Infrared Camera (NIRCam) F277W and F444W imaging obtained through the public CEERS Early Release Science (ERS) program and $\textbf{imcascade}$, an astronomical fitting code that utilizes Multi-Gaussian Expansion, to carry out our size measurements. When compared to the extrapolation of rest-optical size-mass relations for quiescent galaxies at lower redshift, eight out of ten candidates in our sample (80%) are on average more compact by $\sim$40%. Seven out of ten candidates (70%) exhibit rest-frame infrared sizes $\sim$10% smaller than rest-frame optical sizes, indicative of negative color gradients. Two candidates (20%) have rest-frame infrared sizes $\sim$1.4$\times$ larger than rest-frame optical sizes; one of these candidates exhibits signs of ongoing or residual star formation, suggesting this galaxy may not be fully quenched. The remaining candidate is unresolved in both filters, which may indicate an Active Galactic Nuclei (AGN). Strikingly, we observe three of the most massive galaxies in the sample (log(M$_{\star}$/M$_{\odot}$) = 10.74 - 10.95) are extremely compact, with effective radii ${\sim}$0.7 kpc. Our results indicate that quiescent galaxies may be more compact than previously anticipated beyond z>3, even after correcting for potential color gradients. This suggests that the size evolution of quiescent galaxies is steeper than previously anticipated and our current understanding is biased by the limited wavelength capabilities of the $\textit{Hubble Space Telescope}$ (HST) and the presence of negative color gradients in quiescent galaxies.

The Nancy Grace Roman Space Telescope Coronagraph Instrument will enable the polarimetric imaging of debris disks and inner dust belts in the optical and near-infrared wavelengths, in addition to the high-contrast polarimetric imaging and spectroscopy of exoplanets. The Coronagraph uses two Wollaston prisms to produce four orthogonally polarized images and is expected to measure the polarization fraction with measurement errors < 3% per spatial resolution element. To simulate the polarization observations through the Hybrid Lyot Coronagraph (HLC) and Shaped Pupil Coronagraph (SPC), we model disk scattering, the coronagraphic point-response function, detector noise, speckles, jitter, and instrumental polarization and calculate the Stokes parameters. To illustrate the potential for discovery and a better understanding of known systems with both the HLC and SPC modes, we model the debris disks around Epsilon Eridani and HR 4796A, respectively. For Epsilon Eridani, using astrosilicates with 0.37+/-0.01 as the peak input polarization fraction in one resolution element, we recover the peak disk polarization fraction of 0.33+/-0.01. Similarly, for HR 4796A, for a peak input polarization fraction of 0.92+/-0.01, we obtain the peak output polarization fraction as 0.80+/-0.03. The Coronagraph design meets the required precision, and forward modeling is needed to accurately estimate the polarization fraction.

The scalar-induced gravitational waves (SIGW), arising from large amplitude primordial density fluctuations, provide a unique observational test for directly probing the epoch of inflation. In this work, we provide constraints on the SIGW background by taking into account the non-Gaussianity in the primordial density fluctuations, using the third observing run (O3) data of the LIGO-Virgo-KAGRA collaboration. We find that the non-Gaussianity gives a non-negligible effect on the GW energy density spectrum and starts to affect the analysis of the O3 data when the non-Gaussianity parameter is $F_{\rm NL} > 3.55$. Furthermore, the constraints exhibit asymptotic behavior given by $F_{\rm NL} A_g = \rm{const.}$ at large $F_{\rm NL}$ limit, where $A_g$ denotes the amplitude of the curvature perturbations. In the limit of large $F_{\rm NL}$, we placed a 95% confidence level upper limit $F_{\rm NL} A_g \leq 0.13, 0.09, 0.10$ at fixed scales of $10^{16}, 10^{16.5}, 10^{17}~{\rm Mpc}^{-1}$, respectively.

We propose the use of entropy, measured from the spatial and flux distribution of pixels in the residual image, as a potential diagnostic and stopping metric for the CLEAN algorithm. Despite its broad success as the standard deconvolution approach in radio interferometry, finding the optimum stopping point for the iterative CLEAN algorithm is still a challenge. We show that the entropy of the residual image, measured during the final stages of CLEAN, can be computed without prior knowledge of the source structure or expected noise levels, and that finding the point of maximum entropy as a measure of randomness in the residual image serves as a robust stopping criterion. We also find that, when compared to the expected thermal noise in the image, the maximum entropy of the residuals is a useful diagnostic that can reveal the presence of data editing, calibration, or deconvolution issues that may limit the fidelity of the final CLEAN map.

Recently Lian et al. (2023), thanks to Gaia-ESO data, studied the chemical evolution of neutron-capture elements in the regime [Fe/H]>-1. We aim here to complement this study down to [Fe/H]=-3, and focus on Ba, Y, Sr, and abundance ratios of [Ba/Y] and [Sr/Y], which give comprehensive views on s-process nucleosynthesis channels. We measured LTE and NLTE abundances of Ba, Y, and Sr in 323 Galactic metal-poor stars using high-resolution optical spectra with high S/N. We used the spectral fitting code TSFitPy, together with 1D model atmospheres using previously determined LTE and NLTE atmospheric parameters. The NLTE effects are on the order of -0.1 to ~0.2dex depending on the element. T he ratio between heavy and light s-process elements [Ba/Y] varies weakly with [Fe/H] even in the metal-poor regime, consistently with the behavior in the metal-rich regime. The [Ba/Y] scatter at a given metallicity is larger than the abundance measurement uncertainties. Homogeneous chemical evolution models with different yields prescriptions are unable to accurately reproduce the [Ba/Y] scatter at low-[Fe/H]. Adopting the stochastic chemical evolution model by Cescutti & Chaippini (2014) allows to reproduce the observed scatter in the abundance pattern of [Ba/Y] and [Ba/Sr]. With our observations, we rule out the need for an arbitrary scaling of the r-process contribution as previously suggested by the model authors. We have showed how important it is to properly include NLTE effects when measuring chemical abundances, especially in the metal-poor regime. This work shows that the choice of the Galactic chemical evolution model (stochastic vs. 1-zone) is key when comparing models to observations. The upcoming surveys such as 4MOST and WEAVE will deliver high quality spectra of many thousands of metal-poor stars, and this work gives a typical case study of what could be achieved with such surveys.

We search for narrow-line optical emission from dark matter decay by stacking dark-sky spectra from the Dark Energy Spectroscopic Instrument (DESI) at the redshift of nearby galaxies from DESI's Bright Galaxy and Luminous Red Galaxy samples. Our search uses regions separated by 5 to 20 arcsecond from the centers of the galaxies, corresponding to an impact parameter of approximately $50\,\rm kpc$. No unidentified spectral line shows up in the search, and we place a line flux limit of $10^{-19}\,\rm{ergs}/\rm{s}/\rm{cm}^{2}/\rm{arcsec}^{2}$ on emissions in the optical band ($3000\lesssim\lambda\lesssim9000 \,\mathring{\rm A}$), which corresponds to $34$ in AB magnitude in a normal broadband detection. This detection limit suggests that the line surface brightness contributed from all dark matter along the line of sight is two orders of magnitude lower than the measured extragalactic background light (EBL), which rules out the possibility that narrow optical-line emission from dark matter decay is a major source of the EBL.

We report on a campaign on the bright black hole X-ray binary Swift J1727.8$-$1613 centered around five observations by the Imaging X-ray Polarimetry Explorer (IXPE). This is the first time it has been possible to trace the evolution of the X-ray polarization of a black hole X-ray binary across a hard to soft state transition. The 2--8 keV polarization degree slowly decreased from $\sim$4\% to $\sim$3\% across the five observations, but remained in the North-South direction throughout. Using the Australia Telescope Compact Array (ATCA), we measure the intrinsic 7.25 GHz radio polarization to align in the same direction. Assuming the radio polarization aligns with the jet direction (which can be tested in the future with resolved jet images), this implies that the X-ray corona is extended in the disk plane, rather than along the jet axis, for the entire hard intermediate state. This in turn implies that the long ($\gtrsim$10 ms) soft lags that we measure with the Neutron star Interior Composition ExploreR (NICER) are dominated by processes other than pure light-crossing delays. Moreover, we find that the evolution of the soft lag amplitude with spectral state differs from the common trend seen for other sources, implying that Swift J1727.8$-$1613 is a member of a hitherto under-sampled sub-population.

In this study we analyze Type Ia supernovae (SNe Ia) data sourced from the Pantheon+ compilation to investigate late-time physics effects influencing the expansion history, $H(z)$, at redshifts $(z < 2)$. Our focus centers on a time-varying dark energy (DE) model that introduces a rapid transition in the equation of state, at a specific redshift, $z_a$, from the baseline, $\Lambda = -1$, value to the present value, $w_0$, through the implementation of a sigmoid function. The constraints obtained for the DE sigmoid phenomenological parametrization have broad applicability for dynamic DE models that invoke late-time physics. Our analysis indicates that the sigmoid model provides a slightly better, though not statistically significant, fit to the SNe Pantheon+ data compared to the standard LCDM alone. The fit results, assuming a flat geometry and maintaining $\Omega_m$ constant at the 2018-Planck value of $0.3153$, are as follows: $H_0 = 73.3^{+0.2}_{-0.6}$ km s$^{-1}$ Mpc$^{-1}$, $w_{0} = -0.95^{+0.15}_{-0.02}$, $z_a = 0.8 \pm 0.46$. The errors represent statistical uncertainties only. The available SN dataset lacks sufficient statistical power to distinguish between the baseline LCDM and the alternative sigmoid models. A feature of interest offered by the sigmoid model is that it identifies a specific redshift, $z_a = 0.8$, where a potential transition in the equation of state could have occurred. The sigmoid model does not favor a DE in the phantom region ($w_0 < -1$). Further constraints to the dynamic DE model have been obtained using CMB data to compute the distance to the last scattering surface. While the sigmoid DE model does not completely resolve the $H_0$ tension, it offers a transition mechanism that can still play a role alongside other potential solutions.

The IAG solar observatory is producing high-fidelity, ultra-high-resolution spectra (R>500000) of the spatially resolved surface of the Sun using a Fourier Transform spectrometer (FTS). The radial velocity (RV) calibration of these spectra is currently performed using absorption lines from Earth's atmosphere, limiting the precision and accuracy. To improve the frequency calibration precision and accuracy we plan to use a Fabry-Perot etalon (FP) setup that is an evolution of the CARMENES FP design and an iodine cell in combination. To create an accurate wavelength solution, the iodine cell is measured in parallel with the FP. The FP can then be used to transfer the accurate wavelength solution provided by the iodine via simultaneous calibration of solar observations. To verify the stability and precision of the FTS we perform parallel measurements of the FP and an iodine cell. The measurements show an intrinsic stability of the FTS of a level of 1 m/s over 90 hours. The difference between the FP RVs and the iodine cell RVs show no significant trends during the same time span. The RMS of the RV difference between FP and iodine cell is 10.7 cm/s, which can be largely attributed to the intrinsic RV precisions of the iodine cell and the FP (10.2 cm/s and 1.0 cm/s, respectively). This shows that we can calibrate the FTS to a level of 10 cm/s, competitive with current state-of-the-art precision RV instruments. Based on these results we argue that the spectrum of iodine can be used as an absolute reference to reach an RV accuracy of 10 cm/s.

We observed gamma-ray burst (GRB) 221009A using Very Long Baseline Interferomety (VLBI) with the European VLBI Network (EVN) and the Very Long Baseline Array (VLBA), over a period spanning from 40 to 262 days after the initial GRB. The high angular resolution (~mas) of our observations allowed us, for the second time ever after GRB 030329, to measure the projected size, $s$, of the relativistic shock caused by the expansion of the GRB ejecta into the surrounding medium. Our observations support the expansion of the shock with a 3.6 $\sigma$-equivalent significance, and confirm its relativistic nature by revealing an apparently superluminal expansion rate. Fitting a power law expansion model, $s\propto t^a$, to the observed size evolution, we find a slope $a=1.9_{-0.6}^{+0.7}$, which is steeper than expected from either a forward shock (FS) or reverse shock (RS) model, implying an apparent acceleration of the expansion. Fitting the data at each frequency separately, we find different expansion rates, pointing to a frequency-dependent behaviour. We show that the observed size evolution can be reconciled with a RS plus FS in the case of a wind-like circum-burst medium, provided that the two shocks dominate the emission at different frequencies and, possibly, at different times.

An analytical solution of the perturbed equations is obtained, which exists in all ergodic models of collisionless spherical stellar systems with a single length parameter. This solution corresponds to variations of this parameter, i.e., stretching or shrinking the sphere while preserving the total mass. The system remains in an equilibrium state. The simplicity of the solution allows for explicit expressions for the distribution function, potential, and density at all orders of perturbation theory. This, in turn, helps to clarify the concept of perturbation energy, which, being a second-order quantity in amplitude, cannot be calculated in linear theory. It is shown that the correct expression for perturbation energy, constructed taking into account 2nd order perturbations, and the well-known expression for perturbation energy constructed as bilinear form obtained within linear theory from 1st order perturbations do not coincide. However, both of these energies are integrals of motion and differ only by a constant. The obtained solution can be used to control the correctness of codes and the accuracy of calculations in the numerical study of collisionless stellar models.

The Antarctic Impulsive Transient Antenna (ANITA) detector has observed several radio pulses coming from the surface of the ice cap at the South Pole. These pulses were attributed to upward-going atmospheric particle showers instead of the downward-going showers induced by cosmic rays that exhibit a characteristic polarity inversion of the radio signal due to reflection in the ice. Coherent transition radiation from cosmic-ray showers developing in the atmosphere and intercepting the ice surface has been suggested as a possible and alternative explanation of these so-called "anomalous" events. To test this interpretation, we have developed an extension of ZHS, a program to calculate coherent pulses from electromagnetic showers, to deal with showers that transit a planar interface between two homogeneous and dielectric media, including transition radiation. By considering different geometries, it is found that all pulses from air showers intercepting the ice surface and detected at the height of ANITA, display the same polarity as pulses emitted by ultra-high-energy cosmic-ray showers that fully develop in the atmosphere and are reflected on the ice. We find that transition radiation is disfavored as a possible explanation of the anomalous ANITA events.

J191213.72-441045.1 is a binary system composed of a white dwarf and an M-dwarf in a 4.03-hour orbit. It shows emission in radio, optical, and X-ray, all modulated at the white dwarf spin period of 5.3 min, as well as various orbital sideband frequencies. Like in the prototype of the class of radio-pulsing white dwarfs, AR Scorpii, the observed pulsed emission seems to be driven by the binary interaction. In this work, we present an analysis of far-ultraviolet spectra obtained with the Cosmic Origins Spectrograph at the Hubble Space Telescope, in which we directly detect the white dwarf in J191213.72-441045.1. We find that the white dwarf has an effective temperature of 11485+/-90 K and mass of 0.59+/-0.05 solar masses. We place a tentative upper limit on the magnetic field of ~50 MG. If the white dwarf is in thermal equilibrium, its physical parameters would imply that crystallisation has not started in the core of the white dwarf. Alternatively, the effective temperature could have been affected by compressional heating, indicating a past phase of accretion. The relatively low upper limit to the magnetic field and potential lack of crystallisation that could generate a strong field pose challenges to pulsar-like models for the system and give preference to propeller models with a low magnetic field. We also develop a geometric model of the binary interaction which explains many salient features of the system.

Blue supergiants are the brightest stars in their host galaxies and yet their evolutionary status has been a long-standing problem in stellar astrophysics. In this pioneering work, we present a large sample of 59 early B-type supergiants in the Large Magellanic Cloud with newly derived stellar parameters and identify the signatures of stars born from binary mergers among them. We simulate novel 1D merger models of binaries consisting of supergiants with hydrogen-free cores (primaries) and main-sequence companions (secondaries) and consider the effects of interaction of the secondary with the core of the primary. We follow the evolution of the new-born 16--40\,M$_\{odot}$ stars until core-carbon depletion, close to their final pre-explosion structure. Unlike stars which are born alone, stars born from such stellar mergers are blue throughout their core helium-burning phase and reproduce the surface gravities and Hertzsprung-Russel diagram positions of most of our sample. This indicates that the observed blue supergiants are structurally similar to merger-born stars. Moreover, the large nitrogen-to-carbon and oxygen ratios, and helium enhancements exhibited by at least half our data sample are uniquely consistent with our model predictions, leading us to conclude that a large fraction of blue supergiants are indeed products of binary mergers.

We investigate the influence of a specific class of slow Baryon Number Violation (BNV) -- one that induces quasi-equilibrium evolution -- on pulsar spin characteristics. This work reveals how BNV can potentially alter observable parameters, including spin-down rates, the second derivative of spin frequency, and braking indices of pulsars. Moreover, we demonstrate that BNV could lead to anomalies in pulsar timing, along with a wide array of braking indices, both positive and negative. In addition, we examine the possibility of pulsar spin-up due to BNV, which may result in a novel mechanism for the revival of ``dead'' pulsars. We conclude by assessing the sensitivity required for future pulsar timing efforts to detect such BNV effects, thus highlighting the potential for pulsars to serve as laboratories for testing fundamental physics.

The study of core-collapse supernova remnants (SNRs) presents a fascinating puzzle, with intricate morphologies and a non-uniform distribution of stellar debris. Particularly, young remnants (aged less than 5000 years) hold immense value as they can offer crucial insights into the inner processes of the supernova (SN) engine, revealing details about nucleosynthetic yields and large-scale asymmetries arising from the early stages of the explosion. Furthermore, these remnants also bear characteristics that may reflect the nature of their progenitor stars and the interactions between the remnants and the surrounding circumstellar medium (CSM), shaped by the progenitor's mass-loss history. Hence, investigating the connection between young SNRs, parent SNe, and progenitor massive stars can be of paramount importance to delve into the physics of SN engines, and to investigate the final stages of massive star evolution and the elusive mechanisms governing their mass loss. In this contribution, I review recent advances in modeling the path from massive stars to SNe and SNRs achieved by our team. The focus is on investigating the links between the observed physical and chemical properties of SNRs and their progenitor stars and SN explosions. The unraveling of this connection offers us the opportunity to probe the physics of core-collapse SN explosions and the final stages of evolution of massive stars.

Orbital debris presents a growing risk to space operations, and is becoming a significant source of contamination of astronomical images. Much of the debris population is uncatalogued, making the impact more difficult to assess. We present initial results from the first ten nights of commissioning observations with the International Liquid Mirror Telescope, in which images were examined for streaks produced by orbiting objects including satellites, rocket bodies and other forms of debris. We detected 83 streaks and performed a correlation analysis to attempt to match these with objects in the public database. 48\% of these objects were uncorrelated, indicating substantial incompleteness in the database, even for some relatively-bright objects. We were able to detect correlated objects to an estimated magnitude of 14.5 and possibly about two magnitudes greater for the faintest uncorrelated object.

The International Liquid Mirror Telescope (ILMT) project is a scientific collaboration in observational astrophysics between the Li{\`e}ge Institute of Astrophysics and Geophysics (Li{\`e}ge University, Belgium), the Aryabatta Research Institute of observational sciencES (ARIES, Nainital, India) and several Canadian universities (British Columbia, Laval, Montr{\'e}al, Toronto, Victoria and York). Meanwhile, several other institutes have joined the project: the Royal Observatory of Belgium, the National University of Uzbekistan and the Ulugh Beg Astronomical Institute (Uzbekistan) as well as the Pozna{\'n} Observatory (Poland). The Li{\`e}ge company AMOS (Advanced Mechanical and Optical Systems) has fabricated the telescope structure that has been erected on the ARIES site in Devasthal (Uttarakhand, India). It is the first liquid mirror telescope being dedicated to astronomical observations. First light was obtained on 29 April 2022 and commissioning is being conducted at the present time. In this short article, we describe and illustrate the main components of the ILMT. We also highlight the ILMT papers presented during the third BINA workshop, which discuss various aspects of the ILMT science programs.

The 4m International Liquid Mirror Telescope (ILMT) is the first optical survey telescope in India that performs zenithal observations of a 22$'$ wide strip of the sky. To determine the portion of the sky covered by the ILMT during the entire year, we represent the ILMT Field of View (FoV) in three different coordinate systems - galactic, ecliptic, and equatorial. We adopt a constant declination of $+29^{\circ}21'41.4"$ and varying right ascension (RA) ranges corresponding to the Local Sidereal Time (LST). The observations from June to September are hampered due to the monsoon season. The handiness of such representations will allow us to locate a transient event in the ILMT FoV. This will enable prompt follow-up observations with other facilities.

The International Liquid Mirror Telescope (ILMT) is a 4-meter class survey telescope. It achieved its first light on 29$^{\rm th}$ April 2022 and is now undergoing the commissioning phase. It scans the sky in a fixed \ang{;22;} wide strip centred at the declination of $+$\ang{29;21;41.4} and works in \emph{Time Delay Integration (TDI)} mode. We present a full catalog of sources in the ILMT strip derived by crossmatching \textit{Gaia} DR3 with SDSS DR17 and PanSTARRS-1 (PS1) to supplement the catalog with apparent magnitudes of these sources in $g, r$, and $i$ filters. These sources can serve as astrometric calibrators. The release of Gaia DR3 provides synthetic photometry in popular broadband photometric systems, including the SDSS $g, r$, and $i$ bands for $\sim$220 million sources across the sky. We have used this synthetic photometry to verify our crossmatching performance and, in turn, create a subset of the catalog with accurate photometric measurements from two reliable sources.

The 4m International Liquid Mirror Telescope (ILMT) facility continuously scans the same sky strip ($\sim$22$^\prime$ wide) on each night with a fixed pointing towards the zenith direction. It is possible to detect hundreds of supernovae (SNe) each year by implementing an optimal image subtraction technique on consecutive night images. Prompt monitoring of ILMT-detected SNe is planned under the secured target of opportunity mode using ARIES telescopes (1.3m DFOT and 3.6m DOT). Spectroscopy with the DOT facility will be useful for the classification and detailed investigation of SNe. During the commissioning phase of the ILMT, supernova (SN) 2023af was identified in the ILMT field of view. The SN was further monitored with the ILMT and DOT facilities. Preliminary results based on the light curve and spectral features of SN 2023af are presented.

Gravitationally lensed quasars (GLQs) are known to potentially provide an independent way of determining the value of the Hubble-Lema\^{i}tre parameter $H_{0}$, to probe the dark matter content of lensing galaxies and to resolve tiny structures in distant active galactic nuclei. That is why multiply imaged quasars are one of the main drivers for a photometric monitoring with the 4-m International Liquid Mirror Telescope (ILMT). We would like to answer the following questions -- how many multiply imaged quasars should we be able to detect with the ILMT? And how to derive accurate magnitudes of the GLQ images? Our estimation of the possible number of multiply imaged quasars is $15$, although optimistic forecasts predict up to $50$ of them. We propose to use the adaptive PSF fitting method for accurate flux measurements of the lensed images. During preliminary observations in spring 2022 we were able to detect the quadruply imaged quasar - SDSS J1251+2935 in the $\it{i}$ and $\it{r}$ spectral bands.

Nestled in the mountains of Northern India, is a 4-metre rotating dish of liquid mercury. Over a 10-year period, the International Liquid Mirror Telescope (ILMT) will survey 117 square degrees of sky, to study the astrometric and photometric variability of all detected objects. One of the scientific programs will be a survey of variable stars. The data gathered will be used to construct a comprehensive catalog of light curves. This will be an essential resource for astronomers studying the formation and evolution of stars, the structure and dynamics of our Milky Way galaxy, and the properties of the Universe as a whole. This catalog will be an aid in our advance to understanding the cosmos and provide deeper insights into the fundamental processes that shape our Universe. In this work, we describe the survey and give some examples of variable stars found in the early commissioning data from the ILMT.

Low surface brightness (LSB) galaxies make up a significant fraction of the luminosity density of the local universe. Their low surface brightness suggests a different formation and evolution process compared to more-typical high-surface-brightness galaxies. This study presents an analysis of LSB galaxies found in images obtained by the International Liquid Mirror Telescope during the observation period from October 24 to November 1, 2022. 3,092 LSB galaxies were measured and separated into blue and red LSB categories based on their $g'-i'$ colours. In these samples, the median effective radius is 4.7 arcsec, and the median value of the mean surface brightness within the effective radius is 26.1 mag arcsec$^{-2}$. The blue LSB galaxies are slightly brighter than the red LSB galaxies. No significant difference of ellipticity was found between the blue and the red LSB galaxies.

Recent research suggests a correlation between the variability and intrinsic brightness of quasars. If calibrated, this could lead to the use of quasars on the cosmic distance ladder, but this work is currently limited by lack of quasar light curve data with high cadence and precision. The Python photometric data pipeline SunPhot is being developed as part of preparations for an upcoming quasar variability survey with the International Liquid Mirror Telescope (ILMT). SunPhot uses aperture photometry to directly extract light curves for a catalogue of sources from calibrated ILMT images. SunPhot v.2.1 is operational, but the project is awaiting completion of ILMT commissioning.

The present article is based upon an invited talk delivered at the occasion of the inauguration of the 4m International Liquid Mirror Telescope (ILMT) which took place in Devasthal (ARIES, Uttarakhand, India) on 21st of March 2023. We present hereafter a short history of the liquid mirror telescopes and in particular of the 4m ILMT which is the first liquid mirror telescope entirely dedicated to astrophysical observations. We discuss a few preliminary scientific results and illustrate some direct CCD images taken during the first commissioning phase of the telescope. We invite the reader to refer to the series of ILMT poster papers published in these same proceedings of the BINA3 workshop for more details about the instrument, operation, first observations, performance and scientific results.

I consider the Festina Lente Swampland bound and argue taking thermal effects, as for instance occur during reheating, into account significantly strengthens the implications of this bound. I argue that the confinement scale should be higher than a scale proportional to the vacuum energy, while Festina Lente without thermal effects only bounds the confinement scale to be above the Hubble scale. For Higgsing of nonabelian gauge fields, I find that the magnitude of the Higgs mass should be heavier than a bound proportional to the Electroweak scale (or generally the scale set by the Higgs VEV). The measured values of the Higgs in the SM satisfy the bound. A way to avoid the bound being violated during inflation is to have a large number of species becoming light. If one wants the inflationary scale to lie below the species scale in this case, this bounds the inflationary scale to be $\ll 10^5$ GeV. These bounds have phenomenological implications for BSM physics such as GUTs, suggesting for example a weak or absent gravitational wave signature from the GUT Higgsing phase transition.

Direct detection experiments and the interpretation of their results are sensitive to the velocity structure of the dark matter in our galactic halo. In this work, we extend the formalism that deals with such astrophysics-driven uncertainties, originally introduced in the context of dark-matter-nuclear scattering, to include dark-matter-electron scattering interactions. Using mock data, we demonstrate that the ability to determine the correct dark matter mass and velocity distribution is depleted for recoil spectra which only populate a few low-lying bins, such as models involving a light mediator. We also demonstrate how this formalism allows one to test the compatibility of existing experimental data sets (e.g. SENSEI and EDELWEISS), as well as make predictions for possible future experiments (e.g. GaAs-based detectors).

Motivated by the stunning projections for future CMB surveys, we evaluate the amount of dark radiation produced in the early Universe by two-body decays or binary scatterings with thermal bath particles via a rigorous analysis in momentum space. We track the evolution of the dark radiation phase space distribution, and we use the asymptotic solution to evaluate the amount of additional relativistic energy density parameterized in terms of an effective number of additional neutrino species $\Delta N_{\rm eff}$. Our approach allows for studying light particles that never reach equilibrium across cosmic history, and to scrutinize the physics of the decoupling when they thermalize instead. We incorporate quantum statistical effects for all the particles involved in the production processes, and we account for the energy exchanged between the visible and invisible sectors. Non-instantaneous decoupling is responsible for spectral distortions in the final distributions, and we quantify how they translate into the corresponding value for $\Delta N_{\rm eff}$. Finally, we undertake a comprehensive comparison between our exact results and approximated methods commonly employed in the existing literature. Remarkably, we find that the difference can be larger than the experimental sensitivity of future observations, justifying the need for a rigorous analysis in momentum space.

The use of machine learning techniques has significantly increased the physics discovery potential of neutrino telescopes. In the upcoming years, we are expecting upgrade of currently existing detectors and new telescopes with novel experimental hardware, yielding more statistics as well as more complicated data signals. This calls out for an upgrade on the software side needed to handle this more complex data in a more efficient way. Specifically, we seek low power and fast software methods to achieve real-time signal processing, where current machine learning methods are too expensive to be deployed in the resource-constrained regions where these experiments are located. We present the first attempt at and a proof-of-concept for enabling machine learning methods to be deployed in-detector for water/ice neutrino telescopes via quantization and deployment on Google Edge Tensor Processing Units (TPUs). We design a recursive neural network with a residual convolutional embedding, and adapt a quantization process to deploy the algorithm on a Google Edge TPU. This algorithm can achieve similar reconstruction accuracy compared with traditional GPU-based machine learning solutions while requiring the same amount of power compared with CPU-based regression solutions, combining the high accuracy and low power advantages and enabling real-time in-detector machine learning in even the most power-restricted environments.

We report results from an all-sky search of the LIGO data from the third LIGO-Virgo-KAGRA run (O3) for continuous gravitational waves from isolated neutron stars in the frequency band [30, 150] Hz and spindown range of $[{-1} \times 10^{-8}, {+1} \times 10^{-9}]$ Hz/s. This search builds upon a previous analysis of the first half of the O3 data using the same PowerFlux pipeline. We search more deeply here by using the full O3 data and by using loose coherence in the initial stage with fully coherent combination of LIGO Hanford (H1) and LIGO Livingston (L1) data, while limiting the frequency band searched and excluding narrow, highly disturbed spectral bands. We detect no signal and set strict frequentist upper limits on circularly polarized and on linearly polarized wave amplitudes, in addition to estimating population-averaged upper limits. The lowest upper limit obtained for circular polarization is $\sim 4.5 \times 10^{-26}$, and the lowest linear polarization limit is $\sim {1.3} \times 10^{-25}$ (both near 144 Hz). The lowest estimated population-averaged upper limit is $\sim {1.0} \times 10^{-25}$. In the frequency band searched here, these limits improve upon the O3a PowerFlux search by a median factor of $\sim 1.4$ and upon the best previous limits obtained for the full O3 data by a median factor of $\sim 1.1$.

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in non-rotating stars with total relative violation of the divergence-free constraint of $10^{-16}$. It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin-Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above $80\%$ in excess of $10^5$ CPU cores and excellent weak scaling is shown up to $\sim 5 \times 10^5$ CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exascale computers.

In a previous work [1], given a putative vortex, it was determined whether it is non abelian or not by studying its radiation channels. The example considered there was a $SU(2)$ gauge model whose internal orientational modes are described by an sphere $S^2$. The non abelian effects presented in this reference were not very pronounced, due to the compactness of this space. In the present work, this analysis is extended for a vortex whose internal space is non compact. This situation may be realised by semi-local supersymmetric vortices [2]-[9]. As the internal space has infinite volume, a largely energetic perturbation may propagate along the object. A specific configuration is presented, when the internal space is the resolved conifold with its Ricci flat metric. The curious feature about it is that it corresponds to a static vortex, that is, the perturbation is only due to the internal modes. Even being static, the emission of gravitational radiation is in the present case of considerable order. This suggest that the presence of slowly moving objects that can emit a large amount of gravitational radiation is a hint of non abelianity.

In the past decade, dual-phase xenon time projection chambers (Xe-TPCs) have emerged as some of the most powerful detectors in the fields of astroparticle physics and rare-event searches. Developed primarily towards the direct detection of dark matter particles, experiments presently operating deep underground have reached target masses at the multi-tonne scale, energy thresholds around 1\,keV and radioactivity-induced background rates similar to those from solar neutrinos. These unique properties, together with demonstrated stable operation over several years, allow for the exploration of new territory via high-sensitivity searches for a plethora of ultra-rare interactions. These include searches for particle dark matter, for second order weak decays, and the observation of astrophysical neutrinos. We first review some properties of xenon as a radiation detection medium and the operation principles of dual-phase Xe-TPCs together with their energy calibration and resolution. We then discuss the status of currently running experiments and of proposed next-generation projects, describing some of the technological challenges. We end by looking at their sensitivity to dark matter candidates, to second order weak decays and to solar and supernova neutrinos. Experiments based on dual-phase Xe-TPCs are difficult, and, like all good experiments, they are constantly pushed to their limits. Together with many other endeavours in astroparticle physics and cosmology they will continue to push at the borders of the unknown, hopefully to reveal profound new knowledge about our cosmos.

If gravity is fundamentally quantum, any two quantum particles must get entangled with each other due to their mutual interaction through gravity. This phenomenon, dubbed gravity-mediated entanglement, has led to recent efforts of detecting perturbative quantum gravity in table-top experimental setups. In this paper, we generalize this to imagine two idealized massive oscillators, in their ground state, which get entangled due to gravity in an expanding universe, and find that the curvature of the background spacetime leaves its imprints on the resulting entanglement profile. Thus, detecting gravity-mediated entanglement from cosmological observations will open up an exciting new avenue of measuring the local expansion rate of the cosmos.

There has been growing interest in $f(Q)$ gravity, which has led to significant advancements in the field. However, it is important to note that most studies in this area were based on the coincident gauge, thus overlooking the impact of the connection degrees of freedom. In this work, we pay special attention to the connection when studying perturbations in general teleparallel, metric teleparallel, and symmetric teleparallel theories of gravity. We do not just examine perturbations in the metric, but also in the affine connection. To illustrate this, we investigate cosmological perturbations in $f(G)$, $f(T)$, and $f(Q)$ gravity with and without matter in form of an additional scalar field for spatially flat and curved FLRW geometries. Our perturbative analysis reveals that for general $f(Q)$ backgrounds, there are up to seven degrees of freedom, depending on the background connection. This is in perfect agreement with the upper bound on degrees of freedom established for the first time in $\href{https://doi.org/10.1002/prop.202300185}{Fortschr. Phys. 2023, 2300185}$. In $f(G)$ and $f(T)$ gravity theories, only two tensor modes propagate in the gravity sector on generic curved cosmological backgrounds, indicating strong coupling problems. In the context of $f(Q)$ cosmology, we find that for a particular background connection, where all seven modes propagate, there is at least one ghost degree of freedom. For all other choices of the connection the ghost can be avoided at the cost of strong coupling problem, where only four degrees of freedom propagate. Hence, all of the cosmologies within the teleparallel families of theories in form of $f(G)$, $f(T)$, and $f(Q)$ suffer either from strong coupling or from ghost instabilities. A direct coupling of the matter field to the connection or non-minimal couplings might alter these results.

We investigate the prospect of performing a null test of binary black hole (BBH) nature using spin-induced quadrupole moment (SIQM) measurements. This is achieved by constraining a deviation parameter ($\delta\kappa$) related to the parameter ($\kappa$) that quantifies the degree of deformation due to the spin of individual binary components on leading (quadrupolar) spin-induced moment. Throughout the paper, we refer to $\kappa$ as the SIQM parameter and $\delta\kappa$ as the SIQM-deviation parameter. The test presented here extends the earlier SIQM-based null tests for BBH nature by employing waveform models that account for double spin-precession and higher modes. We find that waveform with double spin-precession gives better constraints for $\delta\kappa$, compared to waveform with single spin-precession. We also revisit earlier constraints on the SIQM-deviation parameter for selected GW events observed through the first three observing runs (O1-O3) of LIGO-Virgo detectors. Additionally, the effects of higher-order modes on the test are also explored for a variety of mass-ratio and spin combinations by injecting simulated signals in zero-noise. Our analyses indicate that binaries with mass-ratio greater than 3 and significant spin precession may require waveforms that account for spin-precession and higher modes to perform the parameter estimation reliably.

Gravitational waves (GWs) can be distorted, similarly to electromagnetic (EM) waves, by massive objects through the phenomenon of gravitational lensing. The importance of gravitational lensing for GW astronomy is becoming increasingly apparent in the GW detection era, in which nearly 100 events have already been detected. As current ground-based interferometers reach their design sensitivities, it is anticipated that these detectors may observe a few GW signals that are strongly lensed by the dark halos of intervening galaxies or galaxy clusters. Analyzing strong lensing effects on GW signals is, thus, becoming important to understand the lens' properties and correctly infer the intrinsic GW source parameters. However, one cannot accurately infer lens parameters for complex lens models with only GW observations because there are strong degeneracies between the parameters of lensed waveforms. In this paper, we discuss how to conduct parameter estimation of strongly lensed GW signals and infer the lens parameters with the help of EM observations, such as the shape and the redshift of the lens. We find that for simple spherically symmetric lens models, the lens parameters can be well recovered using only GW information. On the other hand, recovering the lens parameters requires systems in which four or more GW images are detected with additional EM observations for non-axially symmetric lens models. Combinations of GW and EM observations can further improve the inference of the lens parameters.

Recently the Event Horizon Telescope observed black holes at event horizon scales for the first time, enabling us to now test the existence of event horizons. Although event horizons have by definition no observable features, one can look for their non-existence. In that case, it is likely that there is some kind of surface, which like any other surface could absorb (and thermally emit) and/or reflect radiation. In this paper, we study the potential observable features of such rotating reflecting surfaces. We construct a general description of reflecting surfaces in arbitrary spacetimes. This is used to define specific models for static and rotating reflecting surfaces, of which we study the corresponding light paths and synthetic images. This is done by numerical integration of the geodesic equation and by the use of the general relativistic radiative transfer code RAPTOR. The reflecting surface creates an infinite set of ring-like features in synthetic images inside the photon ring. There is a central ring in the middle and higher order rings subsequently lie exterior to each other converging to the photon ring. The shape and size of the ring features change only slightly with the radius of the surface R, spin a and inclination i, resulting in all cases in features inside the 'shadow region'. We conclude that rotating reflecting surfaces have clear observable features and that the Event Horizon Telescope is able to observe the difference between reflecting surfaces and an event horizon for high reflectivities. Such reflecting surface models can be excluded, which strengthens the conclusion that the black hole shadow indeed indicates the existence of an event horizon.