Papers for Monday, May 26 2025

We present an expanded template library for the supernova identification (SNID) software, along with updated source files that make it easy to merge our templates - and other major SNID libraries - into the base code. This expansion, dubbed 'Super-SNID', increases the number of spectra for under-represented supernova classes (e.g., SNe Ia-02cx, Ibn) and adds new classes (e.g., SLSNe, TDEs, LFBOTs). Super-SNID includes 841 spectral templates for 161 objects, primarily from the Public ESO Spectroscopic Survey of Transient Objects (PESSTO) Data Releases 1-4. The library is available on GitHub with simple installation instructions.

In this paper we use the imaging results of M87* from the EHT to calculate the rotational velocity of the inner edge of the accretion disk and find a value of $\sim$0.14c. We then calculate the dimensionless spin parameter, $a$, of the black hole, obtaining a value of $a \sim 0.8$. We deduce that this is probably a lower limit. We go on to use the results of the EHT polarization study of the magnetic field direction in the accretion disk of M87* as a proxy for the direction of motion of the spiralling accreting matter in this highly ionized disk. This direction is defined by the vector sum of the tangential rotation velocity and the inward radial accretion velocity. We thus calculate the accretion velocity to be $\sim$ (7 $\pm$ 0.7) $\times$ $10^7$~ms$^{-1}$. We go on to estimate a range of values for the accretion rate from the inner disk to be $\sim$4 $\times$ 10$^{-5}$ to $\sim$4 $\times$ 10$^{-1}$ M$_{\odot}$yr$^{-1}$, and a range of values for the accretion power to be $\sim$10$^{34}$ to 10$^{38}$ J/s. This is the same range as the power of the jet, making it consistent with accretion-driven jet models.

In Paper I, we classified a stellar sample from the thin disc with a broad range in metallicity as being churned outward or inward, or blurred/undisturbed. In this paper (Paper II), we delve deeper by analysing our entire metallicity-stratified sample along with their dynamic properties, focusing on the connection between radial migration and Li depletion. We analyse the chemo-dynamics of a set of 1188 thin disc dwarf stars observed by the \textit{Gaia}-ESO survey, previously classified into six metallicity-stratified groups via Hierarchical Clustering (HC). We examine several features, such as effective temperatures, masses, and dynamic properties. We also implement a parametric survival analysis using penalised splines (logistic distribution) to quantify how stellar properties and motion (or migration) direction jointly influence Li depletion patterns. We find that stars in our sample that appear to have churned outward are predominantly Li-depleted, regardless of their metallicities. These stars are also the oldest, coldest, and least massive compared to those in the same HC group that have either churned inward or kept their orbital radii. Our survival analysis confirms temperature as the primary driver of Li depletion, followed by metallicity and age, while migration direction shows negligible influence. The increasing proportion of outward-churned stars with higher metallicity (and older ages) indicates their dominant influence on the overall trend observed in the [Fe/H]-A(Li) space for stellar groups with [Fe/H]>0. The survival model reinforces that the observed Li depletion stems primarily from intrinsic stellar properties (cool temperatures, higher metallicity, old ages) rather than migration history. This suggests the metallicity-dependent depletion pattern emerges through stellar evolution rather than Galactic dynamical processes.

Separating resolved and unresolved sources in large imaging surveys is a fundamental step to enable downstream science, such as searching for extragalactic transients in wide-field time-domain surveys. Here we present our method to effectively separate point sources from the resolved, extended sources in the Dark Energy Spectroscopic Instrument (DESI) Legacy Surveys (LS). We develop a supervised machine-learning model based on the Gradient Boosting algorithm \texttt{XGBoost}. The features input to the model are purely morphological and are derived from the tabulated LS data products. We train the model using $\sim$$2\times10^5$ LS sources in the COSMOS field with HST morphological labels and evaluate the model performance on LS sources with spectroscopic classification from the DESI Data Release 1 ($\sim$$2\times10^7$ objects) and the Sloan Digital Sky Survey Data Release 17 ($\sim$$3\times10^6$ objects), as well as on $\sim$$2\times10^8$ Gaia stars. A significant fraction of LS sources are not observed in every LS filter, and we therefore build a ``Hybrid'' model as a linear combination of two \texttt{XGBoost} models, each containing features combining aperture flux measurements from the ``blue'' ($gr$) and ``red'' ($iz$) filters. The Hybrid model shows a reasonable balance between sensitivity and robustness, and achieves higher accuracy and flexibility compared to the LS morphological typing. With the Hybrid model, we provide classification scores for $\sim$$3\times10^9$ LS sources, making this the largest ever machine-learning catalog separating resolved and unresolved sources. The catalog has been incorporated into the real-time pipeline of the La Silla Schmidt Southern Survey (LS4), enabling the identification of extragalactic transients within the LS4 alert stream.

We present extensive follow-up time-series photometry of WD J0049$-$2525, the most massive pulsating white dwarf currently known with $T_{\rm eff} = 13\, 020\,{\rm K}$ and $\log{\it g} = 9.34$ cm s$^{-2}$. The discovery observations detected only two significant pulsation modes. Here, we report the detection of 13 significant pulsation modes ranging from 170 to 258 s based on 11 nights of observations with the New Technology Telescope, Gemini, and Apache Point Observatory telescopes. We use these 13 modes to perform asteroseismology and find that the best-fitting models (under the assumption of an ONe core composition) have $M_{\star} \approx 1.29~M_\odot$, surface hydrogen layer mass of $\log(M_{\rm H}/M_{\star}) \lesssim -7.5$, and a crystallized core fraction of $>99\%$. An analysis of the period spacing also strongly suggests a very high mass. The asteroseismic distance derived is in good agreement with the distance provided by Gaia. We also find tentative evidence of a rotation period of 0.3 or 0.67 d. This analysis provides the first look at the interior of a $\sim 1.3~M_{\odot}$ white dwarf.

Gamma-ray bursts (GRBs) are generally categorized into long and short bursts based on their duration ($T_{90}$). Recently, it has been proposed that GRBs can also be classified into type I (merger) and type II (collapsar) bursts based on the different origin. From a sample of \textit{Swift} long GRBs~(LGRBs) with a redshift completeness of 60\% and $P \geq 2.6 \, \text{ph} \, \text{cm}^{-2} \, \text{s}^{-1}$, collected through the end of 2023, we identify a pure sample of 146 Type II GRBs. With this sample, we construct the luminosity function (LF) using both the Broken Power Law (BPL) and Triple Power Law (TPL) models. Our results indicate that, similar to LGRBs, a strong redshift evolution in either luminosity or density is necessary to accurately account for the observations, regardless of the specific form of the LF assumed. The LF of LGRBs remains a topic of debate, with some studies suggesting it follows a BPL form, while others advocate for a TPL form. In our study, we find that the LF of Type II GRBs tends to favor a BPL model.

We study the detectability of supermassive black holes (SMBHs) with masses of $M_{\bullet}\gtrsim 10^{9}\,\mathrm{M}_\odot$ displaced by gravitational wave recoil kicks $(v_{\rm kick}=0\mathrm{-}2000\,\mathrm{km\,s}^{-1})$ in simulations of merging massive $(M_{\star}>10^{11}\,\mathrm{M}_\odot)$ early-type galaxies. The used KETJU code combines the GADGET-4 fast multiple gravity solver with accurate regularised integration and post-Newtonian corrections (up to PN3.5) around SMBHs. The ejected SMBHs carry clusters of bound stellar material (black hole recoil clusters, BRCs) with masses in the range of $10^6 \lesssim M_{\text{BRC}} \lesssim 10^7\,\mathrm{M}_\odot$ and sizes of several $10\,\mathrm{pc}$. For recoil velocities up to $60\%$ of the galaxy escape velocity, the BRCs are detectable in mock photometric images at a Euclid-like resolution up to redshift $z \sim 1.0$. By Monte Carlo sampling the observability for different recoil directions and magnitudes, we predict that in $\sim20\%$ of instances the BRCs are photometrically detectable, most likely for kicks with SMBH apocentres less than the galaxy effective radius. BRCs occupy distinct regions in the stellar mass/velocity dispersion vs. size relations of known star clusters and galaxies. An enhanced velocity dispersion in excess of $\sigma \sim 600\,\mathrm{km\,s}^{-1}$ coinciding with the SMBH position provides the best evidence for an SMBH-hosting stellar system, effectively distinguishing BRCs from other faint stellar systems. BRCs are promising candidates to observe the aftermath of the yet-undetected mergers of the most massive SMBHs and we estimate that up to 8000 BRCs might be observable below $z\lesssim 0.6$ with large-scale photometric surveys such as Euclid and upcoming high-resolution imaging and spectroscopy with the Extremely Large Telescope.

Ly$\alpha$ emitters (LAEs) have now been discovered out to redshift $z=13$, and are valuable probes of the reionization history at redshifts beyond the reach of other currently available tracers. Most inferences of the neutral hydrogen fraction from LAE observations rely on one-point and two-point statistics like the luminosity function and the power spectrum. We present here an analysis of the bispectrum of high-redshift LAEs and demonstrate its sensitivity to the Epoch of Reionization. We use the Sherwood-Relics suite of cosmological hydrodynamical simulations post-processed with the GPU-based radiative transfer code ATON-HE to generate realistic LAE mock catalogues for a wide range of reionization models, varying the ionization history and the source populations, including contributions of AGN to hydrogen reionization. We demonstrate that the bispectrum provides greater sensitivity than the power spectrum to both the timing and spatial morphology of reionization. Using reduced-$\chi^2$ analysis we further show that the bispectrum also responds more strongly to variations in source population and AGN contribution, apparently more efficiently capturing morphological features missed by two-point statistics. The redshift evolution of the bispectrum reflects the increased clustering of ionizing sources at earlier epochs. The sensitivity of the bispectrum to peculiar velocities underscores the importance of velocity corrections in comparisons to observations. Our findings demonstrate that the LAE bispectrum is a powerful higher-order statistic for probing reionization through current and future LAE surveys using telescopes such as Subaru and JWST.

Atmospheric escape's potential to shape the exoplanet population motivates detailed observations of systems actively undergoing escape. AU Mic is a young and active M dwarf hosting two close-in transiting sub- to Neptune-sized planets. Atmospheric escape was previously detected on the inner planet b, with radially-blown neutral hydrogen producing ~30% blue-shifted absorption in Lyman-alpha. We obtained one HST/STIS transit of the outer planet c, to search for the planet's escaping atmosphere in transmission at Lyman-alpha and compare with AU Mic b. We detected 6 short-duration flares in Si IV and C IV, of which only one corresponded to a Lyman-alpha flare. We identified longer-duration stellar variability at the tens of percent level for lines less sensitive to stellar activity, including O I, C II and Lyman-alpha, which inhibits detection of an exosphere. We do not report absorption associated with an exosphere containing neutral hydrogen or any metals detectable in the far-ultraviolet, and discuss the implications of the non-detection. This work highlights the importance of 1) careful consideration of stellar variability in atmospheric escape observations, and 2) the dual-influence of photoionization and stellar wind when interpreting and modeling atmospheric escape.

In this work, we reconstruct cosmic history via supervised learning through three methods: Classification and Regression Trees (CART), Multi-layer Perceptron Regressor (MLPR), and Support Vector Regression (SVR). For this purpose, we use ages of simulated galaxies based on 32 massive, early-time, passively evolving galaxies in the range $0.12 < z < 1.85$, with absolute ages determined. Using this sample, we simulate subsamples of 100, 1000, 2000, 3334, 6680 points, through the Monte Carlo Method and adopting a Gaussian distribution centering on a spatially flat $\Lambda$CDM as a fiducial model. We found that the SVR method demonstrates the best performance during the process. The methods MLPR and CART also present satisfactory performance, but their mean square errors are greater than those found for the SVR. Using the reconstructed ages, we estimate the matter density parameter and equation of state (EoS) and our analysis found the SVR with 600 predict points obtains $\Omega_m=0.329\pm{}^{0.010}_{0.010}$ and the dark energy EoS parameter $\omega= -1.054\pm{}^{0.087}_{0.126}$, which are consistent with the values from the literature. We highlight that we found the most consistent results for the subsample with 2000 points, which returns 600 predicted points and has the best performance, considering its small sample size and high accuracy. We present the reconstructed curves of galaxy ages and the best fits cosmological parameters.

Recent observations of convection in the jovian atmosphere have demonstrated that convection is strongly concentrated at specific locations on planet. For instance, observations of lightning show that the cyclonic features (e.g,. belts and folded filamentary regions - FFRs) show increased convective activity compared to anti-cyclonic regions. Meanwhile, the distribution of ammonia and water vapor show a large enrichment near the equator, which is also suggestive of strong upwelling and convective activity. Marrying these different observations is challenging due to a lack of data concerning the characteristics of the deep jovian atmosphere, and a resulting inability to observe the true deep source of the various convective phenomena. To understand the nature of these convective events and \paperedit{the role of the } structure of the deep atmosphere \paperedit{in driving convective events}, we run simulations of cloud formation and convection using the Explicit Planetary hybrid-Isentropic Coordinate General Circulation Model (EPIC GCM). We vary the dynamics of the atmosphere by parameterizing the deep wind shear and studying the resulting effect on the strength, frequency and distribution of convective storms. We find that convection in our model is strongly tied to the local dynamics and the deep wind shear. We further decompose the generation of convective available potential energy (CAPE) into three components (thermal, mechanical, and moist/chemical), and find that the chemical mechanism is the strongest component, working to advect water vapor from moisture-rich regions to moisture-poor regions and to drive convection along a ``moisture front.''

The Milky Way Backup Program (MWBP), a survey currently underway with the Dark Energy Spectroscopic Instrument (DESI) on the Nicholas U. Mayall 4-m Telescope, works at the margins of the DESI Main surveys to obtain spectra of millions of additional stars from the Gaia catalog. Efficiently utilizing twilight times (<18 deg) and poor weather conditions, the MWBP extends the range of stellar sources studied to both brighter magnitudes and lower Galactic latitude and declination than the stars studied in DESI's Main Milky Way Survey. While the MWBP prioritizes candidate giant stars selected from the Gaia catalog (using color and parallax criteria), it also includes an unbiased sample of bright stars (i.e., 11.2 < G < 16 mag) as well as fainter sources (to G < 19 mag). As of March 1, 2025, the survey had obtained spectra of ~7 million stars, approximately 1.2 million of which are included in the DESI Data Release 1. The full survey, when completed, will cover an area of more than 21,000 square degrees and include approximately 10 million Gaia sources, roughly equal to the number of stellar spectra obtained through the DESI Main Survey, while only utilizing <9% of all DESI observing time.

It is unknown if an Oort cloud reaches its maximum mass within its star's birth cluster or millions of years later. Complicating the Oort cloud evolution process is the fact that comets can be stripped from orbit due to perturbations from passing stars. We explore how a star's cluster escape time (t$_{ \rm esc}$) and the time its Oort cloud reaches maximum mass (t$_{ \rm max}$) affect the Oort cloud's ability to survive via $N$-body simulations. In a 14 M$_\odot$/pc$^3$ cluster, we identify 50 stars of 1 M$_\odot$ with a range of t$_{ \rm esc}$ to host Oort clouds, each with 1000 comets at t$_{ \rm max}$. For each host, we consider Oort clouds that reach maximum mass 0, 50, and 250 Myr after the cluster's formation. Each Oort cloud's evolution is simulated in the cluster from t$_{ \rm max}$ to t$_{ \rm esc}$. Only a fraction of comets tend to remain in orbit, with this amount depending on t$_{ \rm max}$ and t$_{ \rm esc}$. We observe that 12%, 22%, and 32% of Oort clouds with a t$_{ \rm max}$ of 0, 50 and 250 Myr retain >50% of their comets at t$_{ \rm esc}$, respectively. We find that the fraction of comets stripped has the relationship, $\rm f=m\log_{10}(\frac{t_{ \rm esc}-t_{ \rm max}}{Myr})$ where m = 0.32$\pm$0.04, indicating that the longer the Oort cloud remains in the cluster, the more comets are stripped, with this fraction increasing logarithmically at approximately the same rate for each t$_{ \rm max}$.

We present a cosmological model-independent determination of the Hubble constant, $H_0$, by combining time-delay measurements from seven TDCOSMO systems, Einstein radius measurements, and Type Ia Supernovae data sourced from the Pantheon+ sample. For each lens of time-delay system, we calculate the angular diameter distance $D_{A_l}$ using the product $D^{\textrm{Obs}}(z_l) \cdot D_{A,\Delta t}^{\textrm{Obs}}(z_l, z_s)$, where $D^{\textrm{Obs}}(z_l)$ is reconstructed via Gaussian Processes from 99 Einstein radius measurements, and $D_{A,\Delta t}^{\textrm{Obs}}(z_l,z_s)$ is the time-delay angular distance. We also reconstruct the unanchored luminosity distance $H_0 D_L(z_l)$ from supernova data. By using the cosmic distance duality relation validity, we anchor $D_{A_l}$ and $H_0 D_L(z_l)$ to infer $H_0 = 70.55 \pm 7.44$ km/s/Mpc (68\% CL). Our result, though not resolving the Hubble tension, offers a cosmological model-independent consistency check and highlights the potential of using strong lensing and supernovae data via the cosmic distance duality relation to constrain $H_0$.

Studying metal-poor galaxies is crucial for understanding the physical mechanisms driving the evolution of galaxies. Most observational works in dwarf galaxies employ integral field unit data to investigate gas physics in the entire galaxy body. However, previous studies have not explored the detailed spatially resolved properties of individual extragalactic HII regions. We study the only known HII region in the SagDIG dIrr, a metal-poor dwarf galaxy of the local universe, using VIMOS-IFU/VLT and EFOSC2/NTT archival data. We aim to probe the structure of the gaseous nebulae to (i) give insights into the physical processes shaping its gas dynamics, and (ii) relate these mechanisms to the low gas-phase metallicity content estimated with the direct method. The H$\beta$ emission line map shows two clumps as a biconic-like structure aligned across the same axis. The [OIII]$\lambda$5007 map, on the other hand, shows a more symmetric structure. Radial flux density profiles exhibit signatures of ionized stratification. By comparing the gas structure with the stellar population in the same region, the young stellar population ($<700$ Myr) are found closer to the edges of the biconic-like structure. We estimate 12+log(O/H) = $7.23\pm0.04$ dex with the direct method, considering the whole nebula. Temperature fluctuations are detected, as variations of $T_{e}$ from $\sim 23400$K to $\sim 22000$K, and oxygen abundances from 7.24$\pm0.05$ to 7.50$\pm0.14$ dex, from the centre to the outskirts, respectively. In addition, SagDIG is in line with the low-mass end of the MZR. The observed stratification, the distribution of young MS stars, and the metal-poor content of the HII region suggest that the evolution of SagDIG seems to be actively shaped by stellar feedback processes, likely as the result of extragalactic HII regions being subject to the same gas physics as HII regions of the Milky-Way.

We present a spectroscopic investigation of 25 objects previously reported as possible Planetary Nebulae (PNe) in recent catalogs to obtain their physical properties and to establish their true nature. We found 11 objects showing intense emission lines, 11 where it was not possible to measure $\mathrm{H{\beta}}$, and three where no lines are present. We have used diagnostic diagrams to confirm the true PN nature for eight objects. We obtained elemental abundances for three objects whose values are in agreement with the PNe mean values for our Galaxy. Four objects show [N II] $\lambda$6583 more intense than $\mathrm{H{\alpha}}$, and for two of them, this can be explained by the presence of shocks in the gas. Finally, we report angular sizes based on $\mathrm{H{\alpha}}$ and [O III] $\lambda$5007 emission.

We present four galaxy scale lenses discovered in two JWST blank-fields: the ~ 54 arcmin^2 of the PEARLS North-Ecliptic-Pole Time-Domain Field (NEP TDF) and in the ~ 90 arcmin^2 of CEERS. We perform the search by visual inspection of NIRCam photometric data, obtaining an initial list of 16 lens candidates. We down-select this list to 4 high-confidence lens candidates, based on lens modelling of the image configuration and photometric redshift measurements for both the source and the deflector. We compare our results to samples of lenses obtained in ground-based and space-based lens searches and theoretical expectations. We expect that JWST observations of field galaxies will yield approximately 1 galaxy scale lens every three to five NIRCam pointings of comparable depth to these observations (~ 9 arcmin^2 each). This shows that JWST, compared to other lens searches, can yield an extremely high number of secure lenses per unit area, with redshift and size distributions complementary to lens samples obtained from ground-based and wide-area surveys. We estimate that a single JWST pure-parallel survey of comparable depth could yield $\sim 70$ galaxy scale lenses, with a third of them having z_lens>1 and z_source>3.

The radio jets of radio galaxies in galaxy clusters are often bent due to the ram pressure of the intracluster medium. In this paper we start with a well-defined sample of galaxy clusters and subsequently identifying tailed radio sources in these known environments. Our sample consists of 81 galaxy clusters from the Planck ESZ cluster sample. We present a catalogue of 127 extended cluster radio sources, including brightest cluster galaxies, obtained by visually inspecting Karl G. Jansky Very Large Array (1-2 GHz) observations. We have determined the bending angle of 109 well-structured sources, and classified them accordingly: 84 narrow-angle tailed sources (NATs), 16 wide-angle tailed sources (WATs), and 9 non-bent radio sources. We find a negative correlation between the bending angle and the distance to the cluster centre (impact radius), and we observe that NATs generally have smaller impact radii than the regular galaxy population and WATs. We present a phase-space diagram of tailed radio galaxy velocities and impact radii and find that NATs have a significant excess in the high-velocity and low-impact radius region of phase space, indicating they undergo the largest amount of ram pressure bending. We compared the results from our sample with those for jellyfish galaxies, and suggest that the mechanism responsible for bending the radio tails is similar to the stripping of gas in jellyfish galaxies, although tailed radio galaxies are more concentrated in the centre of the phase space. Finally, we find that NATs and WATs have the same occurrence ratio in merging and relaxed clusters. However, their distribution in the phase-space is significantly different. We report an excess of NATs in the high-velocity and low-impact-radius phase-space region in merging clusters, and an excess of relaxed clusters in the low-velocity and low-impact-radius region.

This study attempts to establish a basis for understanding how methods used in research in solar physics have evolved since World War II (WWII). The goal is to begin to explore if and how the changing research environment affects the training of young scientists, and the future of solar physics research at our institutions. A strategy based upon a sample of 650 PhD theses is used to seek possible trends over 8 decades, with the aim of uncovering any correlations between methods used and measures of success. Necessarily subjective, results depend on how methods are defined, and how success is measured. Although a brief justification of the choices made is attempted, trying mainly to avoid pitfalls such as counting citations, it is clear that further assessment is required. The statistical analysis is based upon necessarily subjective categorization and the inference of likelihoods of two different distributions being drawn from the same underlying distribution. The statistics seem to reflect historical events, such as the Kennedy Moonshot program and the associated SKYLAB mission, with changes delayed by a few years. The data suggest that impactful advances are becoming more rare. Yet the methods used have changed little barring those related to obvious technological advances (e.g. the advent of spacecraft, adaptive optics). A follow-up study to explore the 100,000+ publications in solar physics through machine learning seems warranted.

We report on the statistical analysis of quasar photoevaporation at $z\sim2.2$ by comparing the density of surrounding Ly$\alpha$ Emitters (LAEs) and continuum-selected galaxies, based on the imaging data of Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) and CFHT Large Area $U$-band Deep Survey (CLAUDS). We select 18 quasars from Sloan Digital Sky Survey (SDSS) in the HSC Deep/UltraDeep fields, normalize the LAE/continuum-selected galaxy distribution around each quasar with the quasar proximity size, stack them, and then measure the average densities of the galaxies. As a result, we find that the density of LAEs is $\gtrsim 5 \sigma$ lower than that of continuum-selected galaxies within the quasar proximity region. Within the quasar proximity region, we find that the LAEs with high Ly$\alpha$ equivalent widths (EWs) are less dense than those with low EWs at the 3$\sigma$ level and that LAEs with EW of $\gtrsim150$ Å(rest-frame) are predominantly scarce. Finally, we find that both LAEs and continuum-selected galaxies have smaller densities when they are closer to quasars. We argue that the photoevaporation effect is more effective for smaller dark matter haloes predominantly hosting LAEs, but that it may also affect larger haloes.

High-mass star formation involves complex processes, with the hot core phase playing a crucial role in chemical enrichment and the formation of complex organic molecules. However, molecular inventories in hot cores remain limited. Using data from the ALMA Three-millimeter Observations of Massive Star-forming regions survey (ATOMS), the molecular composition and evolutionary stages of two distinct millimeter continuum sources in the high-mass star forming region G336.99-00.03 have been characterized. MM1, with 19 distinct molecular species detected, along with 8 isotopologues and several vibrationally/torsionally excited states, has been identified as a hot core. MM2 with only 5 species identified, was defined as a HII region. Isotopic ratios in MM1 were derived, with $^{12}$C/$^{13}$C ranging from 16.0 to 29.2, $^{16}$O/$^{18}$O at 47.7, and $^{32}$S/$^{34}$S at 19.2. Molecular abundances in MM1 show strong agreement with other sources and three-phase warm-up chemical models within an order of magnitude for most species. Formation pathways of key molecules were explored, revealing chemical links and reaction networks. This study provides a detailed molecular inventory of two millimeter continuum sources, shedding light on the chemical diversity and evolutionary processes in high-mass star-forming regions. The derived molecular parameters and isotopic ratios offer benchmarks for astrochemical models, paving the way for further investigation into the formation and evolution of complex organic molecules during the hot core phase.

Using DESI DR2 baryon acoustic oscillation distance measurements and Planck cosmic microwave background distance priors, we have measured the dark energy density $\rho_X(z)$ and dark energy equation of state $w_X(z)$ as free functions of redshift (smoothly interpolated from values at $z=0, 1/3, 2/3, 1, 4/3, 2.33$), and find both to be consistent with a cosmological constant, with only deviations of 1.2$\sigma$ for $\rho_X(z)$ and 1.9$\sigma$ for $w_X(z)$ at $z=2/3$. We confirm our earlier finding in Wang & Freese (2006) that $w_X(z)$ is significantly less constrained by data than $\rho_X(z)$. Our results differ noticeably from those of the DESI Collaboration, in which they used the same DESI DR2 data combined with Planck data and found a 3.1$\sigma$ deviation from a cosmological constant, a finding which is the consequence of their assuming the parametrization $w_X(z)=w_0+w_a(1-a)$. Our results indicate that assuming a linear $w_X(z)$ could be misleading and precludes discovering how dark energy actually varies with time at higher redshifts. In our quest to discover the physical nature of dark energy, the most urgent goal at present is to determine definitively whether dark energy density varies with time. We have demonstrated that it is of critical importance to measure dark energy density as a free function of redshift from data. Future galaxy redshift surveys by Euclid and Roman at higher redshifts will significantly advance our understanding of dark energy.

Accretion onto supermassive black hole (SMBH) releases energy via radiation, jets or winds, providing feedback effects on the circumnuclear gas environment. However, the dynamics of accreting gas on sub-parsec scales near the SMBH remains poorly understood. With high-resolution X-ray spectra of Mrk 3, we detect a fast inflowing ionized absorber characterized by redshifted Fe XXV and Fe XXVI absorption lines. Photoionization modeling reveals the inflowing absorber is located at $\lesssim0.04-0.74\rm~pc$ and is decelerating from $6.1\pm0.5\times10^3\rm~km~s^{-1}$ to $3.4\pm0.3\times10^3\rm~km~s^{-1}$ over 11 years. Only $\sim0.6$\%--$3$\% of the inflowing material is estimated to reach the event horizon. This direct evidence of sub-parsec scale fueling inflow bridges the gap between the torus and the outer accretion disk. Furthermore, a $0.86$-keV gas component with sub-solar metallicity ($Z\sim0.22$), outflowing at a velocity of $\sim330\rm~km~s^{-1}$, is detected in the soft X-ray band. It may corresponds to the [O III] biconical outflow in the narrow-line region. A putative ultra-fast disk wind outside our line-of-sight, or clouds within the broad-line region, are promising candidates for decelerating the inflowing absorber.

In large-scale galaxy surveys, particularly deep ground-based photometric studies, galaxy blending is inevitable and poses a potential primary systematic uncertainty for upcoming surveys. Current deblenders predominantly rely on analytical modeling of galaxy profiles, facing limitations due to inflexible and imprecise models. We present a novel approach using a U-net structured transformer-based network for deblending astronomical images, which we term the CAT-deblender. It was trained using both RGB and grz-band images, spanning two distinct data formats from the Dark Energy Camera Legacy Survey (DECaLS) database, including galaxies with diverse morphologies. Our method requires only the approximate central coordinates of each target galaxy, bypassing assumptions on neighboring source counts. Post-deblending, our RGB images retain a high signal-to-noise peak, showing superior structural similarity to ground truth. For multi-band images, the ellipticity of central galaxies and median reconstruction error for the r-band consistently lie within +/-0.025 to +/-0.25, revealing minimal pixel residuals. In our comparison focused on flux recovery, our model shows a mere 1 percent error in magnitude recovery for quadruply blended galaxies, significantly outperforming SExtractor's higher error rate of 4.8 percent. By cross-matching with publicly accessible overlapping galaxy catalogs from the DECaLS database, we successfully deblended 433 overlapping galaxies. Furthermore, we demonstrated effective deblending of 63,733 blended galaxy images randomly selected from the DECaLS database.

Magnetic flux emergence is traditionally considered a key trigger of solar filament eruptions; however, its role in suppressing filament eruptions remains less understood. Using multi-wavelength observations from the Solar Dynamics Observatory, this study investigates a unique case of flux emergence below a quiescent filament from January 3 to 5, 2016, where the newly emerging magnetic flux suppressed rather than promoted the eruption of the filament. It is found that the emerging magnetic bipole within the filament channel directly interacted and reconnected with the overlying filament magnetic field and produced a series of two-sided coronal jets along the filament axis. Instead of eruption, the filament kept stable but broke into two segments at the reconnection site. Further magnetic cancellation or recession of the emerged bipole allowed the filament to recover its original structure. Our analysis results revealed that the flux emergence suppressed the filament eruption by reducing the upward net force. The formation and evolution of filament fine structures, such as filament threads, are closely linked to the reconnection processes between the emerging bipole and the horizontal magnetic field of the filament. This study provides direct observational evidence for the stabilization of solar filaments driven by flux emergence, offering new insights into the dual role of magnetic emergence in triggering and suppressing solar eruptions.

We present microsecond-resolution, coherently-dedispersed, polarimetric measurements of 35 fast radio bursts (FRBs) detected during the Commensal Real-time ASKAP Fast Transients (CRAFT) incoherent sum (ICS) survey with the Australian Square Kilometre Array Pathfinder (ASKAP). We find a wide diversity of time-frequency morphology and polarisation properties broadly consistent with those of currently known non-repeating FRBs. The high S/N and fine time-resolution of our data however reveals a wealth of new information. Key results include (i) the distribution of scattering timescales, ${\tau}_{obs}$, is limited purely by instrumental effects, with no downturn at high ${\tau}_{obs}$ as expected from a log-normal distribution; (ii) for the 29 FRBs with known redshift, ${\tau}_{obs}$ is uncorrelated with dispersion measure (DM) fluctuations about the Macquart relation, in contrast to expectations from pulsar scattering-DM relations; (iii) all FRBs probably have multiple components, and at least a large fraction have variable PA, the identification of which is limited by scattering; (iv) at least half of all FRBs exhibit PA microstructure at 200 ${\mu}s$-200 ns timescales, with behaviour most closely resembling a sub-category of Crab main pulses; (v) that there is a break in the FRB circular polarisation distribution at Stokes V $\gtrsim$ 20%, which is suggestive of a discrete sub-population.

The data release from the Atacama Cosmology Telescope (ACT) imposes stronger constraints on primordial black hole (PBH) formation in single-field inflation models than the Planck data. In particular, the updated Cosmic Microwave Background (CMB) measurements favor a higher scalar spectral index $n_s$, a positive running $\alpha_s$, and a smaller slow-roll parameter $\eta$. Even in the absence of PBH production, the updated constraints together place many single-field models under tension with the CMB data. To explore this tension, we study PBH formation in an $\alpha$-attractor E-model. We investigate an impact of bending of the inflaton potential plateau toward reconciling the model with the new bounds for CMB observables. We find that attempts to increase $n_s$ by upward bending lead to negative values of $\alpha_s$ that are disfavored by the new bounds. A possible model-building approach to resolve the tension is proposed.

We present our observations and analysis of SN~2023gpw, a hydrogen-rich superluminous supernova (SLSN~II) with broad emission lines in its post-peak spectra. Unlike previously observed SLSNe~II, its light curve suggests an abrupt drop during a solar conjunction between $\sim$80 and $\sim$180~d after the light-curve peak, possibly analogous to a normal hydrogen-rich supernova (SN). Spectra taken at and before the peak show hydrogen and helium `flash' emission lines attributed to early interaction with a dense confined circumstellar medium (CSM). A well-observed ultraviolet excess appears as these lines disappear, also as a result of CSM interaction. The blackbody photosphere expands roughly at the same velocity throughout the observations, indicating little or no bulk deceleration. This velocity is much higher than what is seen in spectral lines, suggesting asymmetry in the ejecta. The high total radiated energy ($\gtrsim9\times10^{50}$~erg) and aforementioned lack of bulk deceleration in SN~2023gpw are difficult to reconcile with a neutrino-driven SN simply combined with efficient conversion from kinetic energy to emission through interaction. This suggests an additional energy source such as a central engine. While magnetar-powered models qualitatively similar to SN~2023gpw exist, more modeling work is required to determine if they can reproduce the observed properties in combination with early interaction. The required energy might alternatively be provided by accretion onto a black hole created in the collapse of a massive progenitor star.

We study a sample of narrow-line transients that share characteristics with the Type IIn classified supernova (SN) 1994W, a prototypical member of this class of events, via investigation of their explosion sites and spectrophotometric data. The normalised cumulative rank (NCR) method was used to compare the explosion sites of 10 events to the star-formation distributions of their host galaxies, and to the sites of different evolved massive stars. The resulting sample mean value of NCR$_{\mathrm{H}\alpha} = 0.170 \pm 0.076$ is low, while the NCR$_{\mathrm{NUV}}$ distribution is flat with a mean value of $0.488 \pm 0.084$. The NCR distribution of SN 1994W-like events is consistent with relatively low-mass red supergiants (RSGs) and, despite the small sample size, inconsistent with high-mass stars such as luminous blue variables. To explain the nature of SN 1994W-like transients, interaction between an expanding ejecta and a relatively massive circumstellar medium is likely required, with the latter possibly having been produced by a H envelope ejection via a nuclear flash event, or a luminous red nova (LRN) from a stellar merger; both channels are consistent with low-mass RSGs suggested by the NCR results. In this context, we find the early $-26$ d spectrum from light curve maximum of SN 2003G to share similarities to those of F8-type supergiant stars and LRNe. Finally, based on late-time HST imaging, we set the deepest limits for the surviving precursor of SN 2011ht to $M_{\mathrm{F438W}} > -3.8$ and $M_{\mathrm{F555W}} > -4.0$ mag. This would exclude most supergiants as a non-terminal progenitor, assuming that such a star is not completely obscured by newly formed dust.

We demonstrate that dark matter interactions can profoundly influence stellar nucleosynthesis in the early universe by altering thermodynamic gradients and modifying nuclear reaction rates within primordial stars. Incorporating a dark matter-modified Fermi-Dirac distribution and accounting for localized energy injection from annihilation heating, our model predicts enhanced production of carbon and nitrogen alongside reduced oxygen synthesis. These compositional shifts significantly reshape stellar structure and produce synthetic spectra that closely reproduce the observed characteristics of carbon-enhanced metal-poor (CEMP) stars. Our findings reveal a direct and previously overlooked role of dark matter in driving the chemical evolution of the early cosmos, offering a plausible link between fundamental particle physics and observable astrophysical signatures.

This paper is based on a dataset of many strongly polarized solar lines belonging to the ''second solar spectrum'', i.e. the spectrum near the limb in linear scattering polarization. The observations were done at the Pic du Midi Turret Dome in 2006. The solar spectra were recorded at high spectral resolution (R = 400000) with the spectrograph slit orthogonal to the solar limb, so that $\mu$ = cos$\theta$ continuously varied from 0 .0 to 0.45. The crystal liquid polarimeter delivered the linear polarization rate (Q/I). Strong lines such as CaII 3934 Å, CaI 4227 Å, SrI 4607 Å, SrII 4078 Å, BaII 4554 Å were studied. We measured the Hanle depolarization with the help of models predicting the polarization envelope with no magnetic field and we got values in the range 13-25 Gauss for the unresolved turbulent magnetic field, and we found that it often decreases towards the limb, revealing an altitude gradient. This present analysis was not yet published and spectra shown here become freely available to the research community.

General-purpose large language models, despite their broad capabilities, often struggle with specialized domain knowledge, a limitation particularly pronounced in more accessible, lower-parameter versions. This gap hinders their deployment as effective agents in demanding fields such as astronomy. Building on our prior work with AstroSage-8B, this study introduces AstroSage-70B, a significantly larger and more advanced domain-specialized natural-language AI assistant. It is designed for research and education across astronomy, astrophysics, space science, astroparticle physics, cosmology, and astronomical instrumentation. Developed from the Llama-3.1-70B foundation, AstroSage-70B underwent extensive continued pre-training on a vast corpus of astronomical literature, followed by supervised fine-tuning and model merging. Beyond its 70-billion parameter scale, this model incorporates refined datasets, judiciously chosen learning hyperparameters, and improved training procedures, achieving state-of-the-art performance on complex astronomical tasks. Notably, we integrated reasoning chains into the SFT dataset, enabling AstroSage-70B to either answer the user query immediately, or first emit a human-readable thought process. Evaluated on the AstroMLab-1 benchmark -- comprising 4,425 questions from literature withheld during training -- AstroSage-70B achieves state-of-the-art performance. It surpasses all other tested open-weight and proprietary models, including leading systems like o3, Gemini-2.5-Pro, Claude-3.7-Sonnet, Deepseek-R1, and Qwen-3-235B, even those with API costs two orders of magnitude higher. This work demonstrates that domain specialization, when applied to large-scale models, can enable them to outperform generalist counterparts in specialized knowledge areas like astronomy, thereby advancing the frontier of AI capabilities in the field.

We present the image maps, data reduction, analysis and the first source counts from the Herschel SPIRE Dark Field. The SPIRE Dark Field is an area of sky near the North Ecliptic Pole observed many times during the calibration phase of the Herschel mission in order to characterise the stability of the SPIRE instrument and is subsequently one of the deepest imaged fields of the Universe at far-infrared-submillimetre wavelengths. The SPIRE dark field is concurrent with the Spitzer IRAC Dark Field used for a similar purpose. The final Dark Field map is comprised of 141 individual SPIRE observations in Small Map and Large Map modes defined by a deep inner region approximately 12' in diameter and a slightly shallower surrounding area of diameter ~30'. The depth of both regions reach well below the confusion limit of the SPIRE instrument at 250 microns, 350 microns and 500 microns. Two independent processes are used to extract sources, a standard map based method using the SUSSEXtractor algorithm and a list driven photometry approach using the XID algorithm with the Spitzer MIPS 24 microns catalogue as an input prior. The resulting source counts detect the turnover in the galaxy population with both methods shown to be consistent with previous results from other Herschel surveys, with the XID process reaching approximately twice as deep compared to traditional map based algorithms. Finally, we compare our results with two contemporary galaxy evolution models, again showing a good general agreement with the modelled counts.

Context. The environment plays a crucial role in galaxy evolution, particularly for galaxies infalling into clusters. Accurately estimating the infall times of galaxies from observations can significantly enhance our understanding of the environmental effects on galaxy evolution. Aims. This paper aims to evaluate existing methods for estimating infall times via the $R-V$ diagram, explore possible strategies to improve accuracy in estimating infall times, and discuss fundamental limitations. Methods. We utilize a TNG300-1 simulation and construct the $R-V$ diagram that is directly comparable to the observations. Using the same dataset, we systematically compare four commonly used methods, including the projected radii, caustic profiles, and two discrete methods. A simple linear partition is also considered as a reference. Results. Each method exhibits distinct characteristics. While the linear partition slightly outperforms other methods, all methods suffer from limited accuracy ($\gtrsim 2.6$ Gyr), constrained by the intrinsic dispersion ($2.53$ Gyr) of infall times in the $R-V$ diagram. Given this limit, we explore two potential approaches that can improve accuracy: (1) the infall time dispersion is smaller in more dynamically relaxed clusters, and (2) employing two estimates of infall times instead of one reduces the dispersion to $\lesssim1.5$ Gyr. We further demonstrate that the intrinsic dispersion primarily arises from orbital overlap: galaxies in different orbital phases overlap with each other in the $R-V$ diagram and thus appear indistinguishable. Conclusions. Orbital overlap fundamentally limits the accuracy of infall time estimation. The linear partition approach could be a simple and robust estimation.

Young massive clusters (YMC) and open clusters (OC) are widely considered as potential environments for assembling merging binary stellar-remnant black holes (BBH) via dynamical interactions. However, such moderate mass systems are susceptible to being disrupted by the external tidal field of their host galaxies, potentially limiting their effectiveness as GW sources. In this study, I investigate the formation of BBH mergers in tidally dissolving star clusters. This is achieved with a newly computed grid consisting of 95 evolutionary model star clusters, where the clusters are subjected to a varying extent of tidal stripping. The cluster evolutions are computed with the direct N-body integrator NBODY7 that includes, among others, treatments for post-Newtonian (PN) effects in compact-binary members, mass loss due to stellar evolution, formation of stellar remnants, and tidal stripping. It is found that even strong tidal stripping does not quench the formation of a black hole (BH) core inside a cluster or the formation of dynamical BBH mergers in the system. The overall properties of BBH mergers, e.g., the form of the distribution of merger delay time, primary mass, and mass ratio, and the redshift evolution of merger rate are not significantly altered by the extent of tidal stripping of the parent cluster population. Furthermore, even strongly tidally stripped clusters are capable of dynamically forming Gaia-BH-like detached BH--main-sequence-star binaries that escape into the galactic field. Limitations of the present study and potential future improvements are discussed.

In this work, we explore the production of primordial black holes (PBHs) within the context of $\alpha$-attractor inflationary models, focusing on the preheating phase following inflation. During this phase, self-resonance instabilities arise due to deviations of the inflationary potential from a quadratic form. PBH formation is analyzed using three criteria: (1) the perturbation must lie within the instability band, (2) its characteristic length must exceed the Jeans length, and (3) it must have sufficient time to collapse based on the estimations of massive scalar field spherical collapse in Einstein-de Sitter universe. Based on these criteria, we calculate the PBH mass fraction using the Press-Schechter (PS) and Khlopov-Polnarev (KP) formalisms. Our results show that the PS formalism tends to overestimate PBH abundance during preheating, as it neglects nonspherical effects. In contrast, the KP formalism yields more realistic predictions by incorporating such effects. We provide a detailed comparison with observational constraints from evaporating PBHs. Notably, the PS formalism is excluded by these constraints, which are based on Hawking radiation, while the KP formalism remains viable. These findings underscore the importance of accounting for nonspherical effects and accurate collapse dynamics in studies of PBH formation during preheating.

We study gravitational wave production at strong first order phase transitions, with large-scale, long-running simulations of a system with a scalar order parameter and a relativistic fluid. One transition proceeds by detonations with asymptotic wall speed $v_\text{w}=0.92$ and transition strength $\alpha_n=0.67$, and the other by deflagrations, with a nominal asymptotic wall speed $v_\text{w}=0.44$ and transition strength $\alpha_n=0.5$. We investigate in detail the power spectra of velocity and shear stress and - for the first time in a phase transition simulation - their time decorrelation, which is essential for the understanding of gravitational wave production. In the detonation, the decorrelation speed is larger than the sound speed over a wide range of wavenumbers in the inertial range, supporting a visual impression of a flow dominated by supersonic shocks. Vortical modes do not contribute greatly to the produced gravitational wave power spectra even in the deflagration, where they dominate over a range of wavenumbers. In both cases, we observe dissipation of kinetic energy by acoustic turbulence, and in the case of the detonation an accompanying growth in the integral scale of the flow. The gravitational wave power approaches a constant with a power law in time, from which can be derived a gravitational wave production efficiency. For both cases this is approximately $\tilde{\Omega}^\infty_\text{gw} \simeq 0.017$, even though they have quite different kinetic energy densities. The corresponding fractional density in gravitational radiation today, normalised by the square of the mean bubble spacing in Hubble units, for flows which decay in much less than a Hubble time, is $\Omega_{\text{gw},0}/(H_\text{n} R_*)^2=(4.8\pm1.1)\times 10^{-8}$ for the detonation, and $\Omega_{\text{gw},0}/(H_\text{n} R_*)^2=(1.3\pm0.2)\times 10^{-8}$ for the deflagration.

Ambient dark matter (DM) around binary black holes can imprint characteristic signatures on gravitational waves emitted from their merger. The exact signature depends sensitively on the DM density profile around the black holes. We run very high resolution cosmological hydrodynamics simulations of first star formation that follow the collapse of a $3\times10^{5}\,M_\odot$ mini-halo from $z=49$ to $z\simeq22$. Our flagship model achieves a DM particle mass of $3.7\times10^{-4}\,M_\odot$ and resolves the inner-most structure down to $0.02\,$pc. We show that the halo experiences a two-stage gravitational collapse, where a rotating, constant-density core with $r\lesssim3\,$pc is formed first, surrounded by an extended outskirts. Baryonic infall toward the center continues to raise the local Keplerian velocity and promotes adiabatic contraction of DM. The resulting density profile has an approximately power-law shape of $\rho_{\rm dm} \propto r^{-0.6}$ inside $\sim\!1\,$this http URL find that a piecewise power-law fit reproduces the simulation result to better than 10\%, and also find numerical convergence down to $\simeq\!0.01\,$pc. The DM density profile is typical for ordinary Pop~III halos, but our additional simulations reveal that inner slope varies significantly with halo-to-halo scatter, and the effect of Lyman-Werner irradiation and of supersonic baryon-DM streaming velocities, implying a wide distribution of slopes rather than a single universal curve. The large variation should be considered when calculating the predicted DM-induced dephasing of gravitational waves by up to an order of magnitude relative to the classical analytic model of the DM spike.

In this manuscript, an improved Seyfert-LINER classification line (= S-L line) is proposed in the [N~{\sc ii}] BPT diagram, based on a sample of 47,968 low redshift narrow emission line galaxies from SDSS DR16, motivated by different S-L lines reported in the [N~{\sc ii}] BPT diagram through different methods. The method proposed by Kewley et al. in 2006 is firstly applied, however, the method cannot lead to an accepted S-L line in the [N~{\sc ii}] BPT diagram. Meanwhile, the S-L lines proposed by Schawinski et al. in 2007 and Cid Fernandes et al. in 2010 in the [N~{\sc ii}] BPT diagram are different from each other. Therefore, it is meaningful to check which proposed S-L line is better or to determine an improved one in the [N~{\sc ii}] BPT diagram by a new method. In this manuscript, Seyferts and LINERs that have already been classified in the [S~{\sc ii}] and/or [O~{\sc i}] BPT diagrams can be visualized in the [N~{\sc ii}] BPT diagram, leading the intersection boundary of the two contour maps to be considered as the S-L line in the [N~{\sc ii}] BPT diagram. Rather than the previously proposed S-L lines, the new S-L line can lead to more efficient and harmonious classifications of Seyferts and LINERs, especially in the composite galaxy region, in the [N~{\sc ii}] BPT diagram. Furthermore, based on the discussed S-L lines, the number ratio of Type-2 Seyferts to Type-2 LINERs differs significantly from that of Type-1 Seyferts to Type-1 LINERs in the [N~{\sc ii}] BPT diagram, suggesting that about 90$\%$ of Type-2 LINERs are non-AGN-related objects, true Type-2 AGNs, or objects exhibiting both Seyfert and LINER characteristics.

Recent observations of the repeating fast radio burst FRB 20201124A by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) revealed a second-scale periodic modulation ($\sim$1.7\,s) in burst activity during two distinct observational windows. We find that these two periodic activity epochs temporally coincide with the transitional states of the source's Faraday rotation measure (RM), and the chance coincidence is only about 0.07$\%$. This correlation is can be understood within the magnetar/Be-star binary system framework. Considering that only the polar cap region can remain stable for such an extended period, we apply a coherent linear periodic evolution model to jointly constrain the initial burst period \( P_0 \) and the period derivative \( \dot{P} \) across both observation windows (MJD 59310 and MJD 59347). We obtain spin parameters consistent with blind search results: an initial spin period $P_0 = 1.7060155$\,s at the reference time and spin period derivative $\dot{P} = 6.1393 \times 10^{-10}$\,s\,s$^{-1}$. We conclude that during these two observational windows, the magnetar was just crossing the disk of the Be star. The disk-magnetar interaction at these two geometric positions interaction may surpress the multi-polar magnetic fields at low latitudes of the magnetar, which enhances the dominance of the polar cap region emissions and makes the periodic activity detectable.

The Barcelona Raman LIDAR (BRL) will provide continuous monitoring of the aerosol extinction profile along the line of sight of the Cherenkov Telescope Array Observatory (CTAO). It will be located at its Northern site (CTAO-N) on the Observatorio del Roque de Los Muchachos. This article presents the performance of the pathfinder Barcelona Raman LIDAR (pBRL), a prototype instrument for the final BRL. Power budget simulations were carried out for the pBRL operating. under various conditions, including clear nights, moon conditions, and dust intrusions. The LIDAR PreProcessing (LPP) software suite is presented, which includes several new statistical methods for background subtraction, signal gluing, ground layer and cloud detection and inversion, based on two elastic and one Raman lines. Preliminary test campaigns were conducted, first close to Barcelona and later at CTAO-N, albeit during moonlit nights only. The pBRL, under these non-optimal conditions, achieves maximum ranges up to about 35 km, range resolution of about 50 m for strongly absorbing dust layers, and 500 m for optically thin clouds with the Raman channel only, leading to similar resolutions for the LIDAR ratios and Angstrom exponents. Given the reasonable agreement between the extinction coefficients obtained from the Raman and elastic lines independently, an accuracy of aerosol optical depth retrieval in the order of 0.05 can be assumed with the current setup. The results show that the pBRL can provide valuable scientific results on aerosol characteristics and structure, although not all performance requirements could be validated under the conditions found at the two test sites. Several moderate hardware improvements are planned for its final upgraded version [truncated].

The recent detection of very high energy (VHE, $>$ 100 GeV) $\gamma$-ray emission from gamma-ray bursts (GRBs) has provided new insights into afterglow physics. Understanding the temporal and spectral evolution of VHE GRBs requires detailed modelling of multiwavelength observations spanning radio to VHE $\gamma$ rays. Previous studies on afterglow emission of VHE GRBs were interpreted using a range of frameworks, including single- and multi-zone jet configurations, synchrotron radiation from forward and reverse shocks, synchrotron self-Compton (SSC) processes, as well as hadronic emission processes. We conducted a detailed multiwavelength modelling of five long-duration VHE GRBs - GRB 180720B, GRB 190114C, GRB 190829A, GRB 201216C and GRB 221009A; using the NAIMA code and modifications to it. The code deals with single-zone synchrotron radiation, SSC and external Compton (EC) radiation with complete Klein Nishina cross-section, extragalactic background light correction, and the Markov chain Monte Carlo techniques. Our analysis constrains key parameters governing the emission and surrounding environments of these GRBs. The results indicate that SSC is the dominant VHE emission mechanism, with negligible contribution from EC. Most VHE GRBs are well described by the forward shock model in a spherical jet configuration, where constant density interstellar medium is preferred over wind medium. Additionally, we find that VHE GRBs tend to occur in environments with lower magnetic fields and higher ambient medium densities. Interestingly, VHE GRBs lie at the edge of the 3-$\sigma$ region of the $E_{\rm k,iso}$ - $\epsilon_B$ correlation observed in other energetic GRBs. Our model slightly over predicts the radio fluxes, indicating that a more complicated modelling might be required in some cases. These findings provide crucial constraints on VHE GRB emission sites and mechanisms (Abridged).

Astrometry, the precise measurement of stellar positions and velocities, offers a promising approach to probing the low-frequency stochastic gravitational wave background (SGWB). Notably, astrometric vector sky maps are sensitive to parity-violating SGWB signals, which cannot be distinguished using pulsar timing array observations in an isotropic SGWB. We present the first astrometric constraints on parity-violating SGWB using quasar catalogs from Gaia DR3 and VLBA data. By analyzing the $EB$ correlation in the two-point correlation function of the proper motions of the quasars, we find 2$\sigma$ constraints on the parity-violating SGWB amplitude $h_{70}^2\Omega_{V} = -0.020 \pm 0.025$ from Gaia DR3 and $h_{70}^2\Omega_{V} = -0.004 \pm 0.010$ from VLBA. These constraints are valid in the frequency range $4.2 \times 10^{-18}\,{\rm Hz} < f < 1.1 \times 10^{-8}\,{\rm Hz}$. Although not currently a tight constraint on theoretical models, this first attempt lays the groundwork for future investigations using more precise astrometric data.

Of the few thousand discovered exoplanets, a significant number orbit in the habitable zone of their star. Many of them are gas giants lacking a rocky surface and solid water reservoirs necessary for life as we know it. The search for habitable environments may extend to the moons of these giant planets. No confirmed exomoon discoveries have been made as of today, but promising candidates are known. Theories suggest that moon formation is a natural process in planetary systems. We aim to study moon formation around giant planets in a phase similar to the final assembly of planet formation. We search for conditions for forming the largest moons with the highest possibility in circumplanetary disks, and investigate whether the resulting moons can be habitable. We determined the fraction of the circumplanetary disk's mass converted into moons using numerical N-body simulations where moon embryos grow via embryo-satellitesimal collisions, investigated in disks around giant planets consisting of 100 fully interacting embryos and 1000 satellitesimals. In fiducial simulations, a 10 Jupiter-mass planet orbited a solar analog star at distances of 1-5 au. To determine the habitability of the synthetic moons, we calculated the stellar irradiation and tidal heating flux on these moons based on their orbital and physical parameters. The individual moon mass is found to be higher when the host planet orbits at a smaller stellar distance. However, moons leave the circumplanetary disk due to the stellar thief effect, which is stronger closer to the star. We find that 32% of synthetic moons can be habitable in the circumstellar habitable zone. Due to the intense tidal heating, the incidence rate of moon habitability is similar at 2 au, and decreasing to 1% at larger distances (<5 au). We conclude that the circumstellar habitable zone can be extended to moons around giant planets.

Modified gravity theories have been suggested to address the limitations of general relativity, each exhibiting differences, particularly in their strong-field limits. Nonetheless, there lacks effective means to distinguish or test these theories through local strong-field measurements. In this work, we define a global Gaussian bending measure over singular spacetime regions, establish a corresponding global theory, and demonstrate its applications in a general stationary spacetime. The global theory is based on differential geometry, rather than on specific gravity theories, allowing it to depict various physics within general relativity and beyond. For example, it can be applied to describe the gravitational bending of massless or massive messengers, such as photons, neutrinos, cosmic rays, and possibly massive gravitational waves predicted in certain theories of gravity. Besides, the global theory is applicable to any stationary spacetime regions outside a rotating black hole. As an instance of its direct applications, we investigate the highly-curved spacetime effects of the black hole in its immediate surrounding regions and design local strong-field experiments involving different shapes of singular lensing patches. New means can be therefore anticipated to be developed according to the global theory to differentiate between different gravity theories and test them in their strong-field regions.

The displacement-noise-free interferometer (DFI) is designed to eliminate all displacement-induced noise while retaining sensitivity to gravitational wave (GW) signals. Ground-based DFIs suffer from physical arm-length limitations, resulting in poor sensitivity at frequencies below 1 kHz. To address this, previous research introduced a neutron-based DFI, which replaces laser light with neutrons and achieves exceptional sensitivity down to a few hertz. In this study, we conducted a proof-of-principle experiment using a pulsed neutron source at the Japan Proton Accelerator Research Complex (J- PARC). Despite practical constraints that led to deviations from the ideal experimental design, we optimized the setup and developed a novel analysis method that successfully cancels displacement noise while preserving simulated GW signals. This work presents the first successful demonstration of a neutron DFI and a neutron interferometer for GW detection.

The Laser Interferometer Space Antenna (LISA) is a planned space-based observatory designed to detect gravitational waves (GWs) within the millihertz frequency range. LISA is anticipated to observe the inspiral of compact objects into black holes at the centers of galaxies, so called extreme-mass-ratio inspirals (EMRIs). However, the extraction of these long-lived complex signals is challenging due to the large size and multimodality of the search space. In this study, we introduce a new search strategy that allows us to find EMRI signals in noisy data from wide priors all the way to performing parameter estimation. This work is an important step in understanding how to extract EMRIs from future LISA data.

Gravitational waves undergo redshift as they propagate through the expanding universe, and the redshift may exhibit time-dependent drift. Consequently, for any isolated gravitational wave sources, the mass parameter $\mathcal{M}$ and the redshift $z$ exhibit an observational degeneracy, typically manifesting in the waveform as the redshifted mass $\mathcal{M}(1+z)$. Matching together the wave propagation and the wave generation solutions, we show that dimensionless source parameters depending on mass $\mathcal{M}$ can break this degeneracy. Notably, the postmerger signal from binary neutron stars contains several dimensionless parameters that satisfy this condition, including the quality factors of different frequency components and their frequency ratios. Considering the observations of solely the postmerger signal by the Neutron star Extreme Matter Observatory or the Einstein Telescope, based on the Fisher analysis, we find that the redshift can be measured with fractional uncertainties of $\sim30\%$ for sources at $0.01<z<0.09$. Additionally, we present a corrected derivation of the waveform phase correction due to the redshift drift effect, rectifying a sign error in previous studies.

We investigate the standard relativistic geometrically thin and optically thick accretion disk in the background of a deformed compact object. The main purpose of this work is to determine whether such a deformed object possesses its own observational fingerprint that can distinguish it from Schwarzschild and Kerr black holes. Our analysis reveals the properties of this relativistic accretion disk model and its dependence on the initial parameters.

During inflation, higher derivative terms in the gravitational action may play a significant role. Building on new stable formulations of four-derivative scalar-tensor theories, we study the impact of these corrections in the case where the inflaton is also the additional scalar degree of freedom of the modified theory. This case is highly restricted by requiring that in the homogeneous limit inflation must still work, and that the initial data must be in the weak coupling limit to respect the validity of the effective theory. In such cases, the non-linear dynamics of large perturbations are very similar to the GR case, with the main deviations captured by the terms relating to the homogeneous Einstein-scalar-Gauss-Bonnet contributions. We show that in principle it is possible to dynamically drive the field out of the weak-coupling regime from a starting point well within it, but that to do so one has to finely tune the setup, so such cases are unlikely to occur generically. This work provides a basis for the study of less restricted models in future, for example those in which the inflaton and scalar degree of freedom are independent, or in which one is not restricted to a weakly coupled regime.

Through the Hamilton-Jacobi equation of classical mechanics, BPS magnetized Baryonic layers (possessing both baryonic charge and magnetic flux) have been constructed in the gauged non-linear sigma model (G-NLSM), which is one of the most relevant effective theories for Quantum Chromodynamics (QCD) in the strongly interacting low-energy limit. Since the topological charge that naturally appears on the right hand side of the BPS bound is a non-linear function of the baryonic charge, the thermodynamics of these magnetized Baryonic layers is highly non-trivial. In this work, using tools from the theory of Casimir effect, we derive analytical relationship between baryonic charge, topological charge, magnetic flux and relevant thermodynamical quantities (such as pressure, specific heat and magnetic susceptibility) of these layers. The critical Baryonic chemical potential is identified. Quite interestingly, the grand canonical partition function can be related with the Riemann zeta function. On the technical side, it is quite a remarkable result to derive explicit expressions for all these thermodynamics quantities of a strongly interacting magnetized system at finite Baryon density. On the physical side, such expressions allow a direct physical interpretation which we discuss.

Cosmic-ray transport in turbulent astrophysical environments remains a multifaceted problem, and despite decades of study the impact of complex magnetic field geometry -- evident in simulations and observations -- has only recently received more focused attention. To understand how ensemble-averaged transport behavior emerges from the intricate interactions between cosmic rays and structured magnetic turbulence, we run test-particle experiments in snapshots of a strongly turbulent magnetohydrodynamics simulation. We characterize particle-turbulence interactions via the gyro radii of particles and their experienced field-line curvatures, which reveals two distinct transport modes: magnetized motion, where particles are tightly bound to strong coherent flux tubes and undergo large-scale mirroring; and unmagnetized motion characterized by chaotic scattering through weak and highly tangled regions of the magnetic field. We formulate an effective stochastic process for each mode: compound subdiffusion with long mean free paths for magnetized motion, and a Langevin process with short mean free paths for unmagnetized motion. A combined stochastic walker that alternates between these two modes accurately reproduces the mean squared displacements observed in the test-particle data. Our results emphasize the critical role of coherent magnetic structures in comprehensively understanding cosmic-ray transport and lay a foundation for developing a theory of geometry-mediated transport.