Papers for Friday, Jan 17 2025

Planetesimal belts are ubiquitous around nearby stars, and their spatial properties hold crucial information for planetesimal and planet formation models. We present resolved dust observations of 74 planetary systems as part of the REsolved ALMA and SMA Observations of Nearby Stars (REASONS) survey and archival reanalysis. We uniformly modelled interferometric visibilities for the entire sample to obtain the basic spatial properties of each belt, and combined these with constraints from multi-wavelength photometry. We report key findings from a first exploration of this legacy dataset: (1) Belt dust masses are depleted over time in a radially dependent way, with dust being depleted faster in smaller belts, as predicted by collisional evolution. (2) Most belts are broad discs rather than narrow rings, with much broader fractional widths than rings in protoplanetary discs. We link broad belts to either unresolved substructure or broad planetesimal discs produced if protoplanetary rings migrate. (3) The vertical aspect ratios (h = H/R) of 24 belts indicate orbital inclinations of 1-20 degrees, implying relative particle velocities of 0.1-4 km/s, and no clear evolution of heights with system age. This could be explained by early stirring within the belt by large bodies (with sizes of at least 140 km to the size of the Moon), by inheritance of inclinations from the protoplanetary disc stage, or by a diversity in evolutionary pathways and gravitational stirring mechanisms. We release the REASONS legacy multidimensional sample of millimetre-resolved belts to the community as a valuable tool for follow-up multi-wavelength observations and population modelling studies.

Galaxy groups with $M_{tot} \lesssim 10^{14}$ $M_\odot$ and up to a few tens of members are the most common galaxy environment, marking the transition between field and massive clusters. Identifying groups plays a crucial role in understanding structure formation and galaxy evolution. Modern deep surveys allow us to build well-characterized samples of groups up to the regime where structures were taking shape. We aimed to build the largest deep catalog of galaxy groups to date over the COSMOS-Web field effective area of 0.45 deg$^2$, leveraging the deep high quality data of the new COSMOS-Web photometric catalog resulted from the James Webb Space Telescope observations of the COSMOS-Web field. We performed the group search with the AMICO algorithm, a linear matched filter based on an analytical model for the group signal. AMICO has already been tested in wide and deep field surveys, including COSMOS data up to $z=2$. In this work, we tested the algorithm performances at even higher redshift and searched for protocluster cores at $z>2$. We compiled a list of known protoclusters in COSMOS at $2 \leq z \leq 3.7$, matched them with our detections and studied the clustering of the detected cores. We estimated purity and completeness of our sample by creating data-driven mocks with the SinFoniA code and linked signal-to-noise to purity. We detected 1678 groups in the COSMOS-Web field up to $z=3.7$, including lists of members extending nearly two magnitudes deeper than the previous AMICO-COSMOS catalog. 756 groups were detected with purity of 80\%. More than 500 groups have their redshift confirmed by assigning spectroscopic counterparts. This group catalog offers a unique opportunity to explore galaxy evolution in different environments spanning $\sim$12 Gyr and to study groups, from the least rich population to the formation of the most massive clusters.

Quiescent galaxies over $3<z<6$ are rare and puzzling. They formed and quenched within two billion years and simulations routinely struggle to predict their observed abundances. Developing a robust identification technique for these galaxies is crucial for constraining galaxy evolution models. Traditional rest-frame color-color selection techniques for quiescent galaxies are known to break down or require adjustments at $z\gtrsim3$. Recently, observed-frame color-color criteria have been established with JWST/NIRCam colors that efficiently pre-select high-redshift quiescent galaxies using only $\lesssim1\%$ of a given sample. In this work, the Uniform Manifold Approximation and Projection machine-learning technique is applied to pre-select quiescent galaxies over $3<z<6$ using observed NIRCam colors. From a parent sample of 43,926 galaxies in JADES, we ultimately find 44 quiescent candidates from a pool of $\approx2,300$ galaxies. This is about five times fewer galaxies than what would be pre-selected using color-color criteria. Two-thirds of these candidates can be pre-selected from a pool as small as 247, which is about twice as efficient as existing observed-frame color selection techniques. Nearly two-thirds of the candidates are new discoveries and include quiescent galaxies with mass-weighted ages as young as $\lesssim300$ Myr. We obtain number densities in agreement with the literature at $z<4$ and find generally higher abundances at $z>4$, although our measurements are consistent within errors. This technique may be applied to other JWST surveys.

The galaxy cluster PSZ2 G181.06+48.47 ($z=0.234$) is a post-merging system that exhibits symmetric double radio relics separated by ~2.7 Mpc. We present the first weak-lensing analysis of PSZ2 G181.06+48.47 and propose possible merging scenarios using numerical simulations. Our analysis with Subaru Hyper Suprime-Cam imaging identifies a binary dark matter structure consisting of northern and southern components, separated by ~500 kpc. Assuming Navarro-Frenk-White (NFW) halos, the masses for the northern and southern subclusters are $M_{200c}^{N} = 0.88_{-0.30}^{+0.35} \times 10^{14} M_{\odot}$ and $M_{200c}^{S} = 2.71_{-0.48}^{+0.51} \times 10^{14} M_{\odot}$, respectively. By superposing the two NFW halos, we determine the total mass of the cluster to be $M_{200c} = 4.22_{-1.00}^{+1.10} \times 10^{14} M_{\odot}$ ($M_{500c} = 2.90_{-0.69}^{+0.75} \times 10^{14} M_{\odot}$). Our mass estimate suggests that the two relics are located around the cluster $R_{200c}$, where the density of the intracluster medium is very low. Our idealized simulations find that an off-axis collision of a 3:1 major merger can simultaneously reproduce the observed relic and dark matter halo separations. From these findings, we suggest that the system is observed ~0.9 Gyr after the first pericenter passage and is returning from the first apocenter.

Merging binary black holes (BBHs) formed dynamically in dense star clusters are expected to have uncorrelated spin--orbit orientations since they are assembled through many random interactions. However, measured effective spins in BBHs detected by LIGO/Virgo/KAGRA hint at additional physical processes that may introduce anisotropy. Here we address this question by exploring the impact of stellar collisions, and accretion of collision debris, on the spin--orbit alignment in merging BBHs formed in dense star clusters. Through hydrodynamic simulations, we study the regime where the disruption of a massive star by a BBH causes the stellar debris to form individual accretion disks bound to each black hole. We show that these disks, which are randomly oriented relative to the binary orbital plane after the initial disruption of the star, can be reoriented by strong tidal torques in the binary near pericenter passages. Following accretion by the BHs on longer timescales, BBHs with small but preferentially positive effective spin parameters ($\chi_{\rm eff} \lesssim 0.2$) are formed. Our results indicate that BBH collisions in young massive star clusters could contribute to the observed trend toward small positive $\chi_{\rm eff}$, and we suggest that the standard assumption often made that dynamically assembled BBHs should have isotropically distributed BH spins is not always justified.

We use NIRSpec MSA spectroscopy and NIRCam Photometry to explore the properties of JADES-GS8-RL-1, a rapidly quenched, $z=8.5$ galaxy with a stellar mass of $10^{8.9}M_\odot$, a steep blue UV slope, a Balmer break, and no sign of strong emission lines. With a $\beta_{UV}$=-2.8$\pm 0.2$, as measured from the NIRSpec spectrum, JADES-GS8-RL-1 is consistent with negligible dust attenuation and little to no contribution from the nebular continuum alongside a probable high escape fraction. The $\beta_{UV}$ slope measured from photometry varies from -3.0 in the central regions to -2.2 at the outskirts suggesting possible regional differences in the escape fraction. There are no high-ionisation emission lines, only a tentative 2.9\sig detection of [OII]. Using photometry, this emission appears to be extended, possibly corresponding to weakly ionised gas expelled during or after the quenching process. JADES-GS8-RL-1 is spatially resolved with a half-light radius of 240 pc and has an exponential, disc-like morphology. It appears to have formed all its stars in a short burst within the past 100 Myr with a formation time of $\approx$70 Myr and a quenching time of $\approx$30 Myr. This quenching would have occurred rapidly, making it a more distant example of the kind of low-mass "mini-quenched" galaxies previously observed at high-z. Due to the extremely blue $\beta_{UV}$ slope, our best-fit model predicts a high value for \fesc of >10\%, consistent with the value derived from the $\beta_{UV}$ slope, which when combined with our extraordinarily low O32 upper limit suggests JADES-GS8-RL-1 is a fascinating example of a high-z "remnant leaker" in one of its earliest phases, deep in the epoch of reionisation.

The James Webb Space Telescope, launched in 2021, is an infrared observatory of novel design: deployable, with active optics, fully open to space for radiative cooling and orbiting the Lagrange point no. 2. This article explains the rationale leading to this specific design and describes the various other architectures that were considered along the way: from a monolithic 10-meter telescope in geosynchronous orbit to a 6-meter one in High Earth Orbit, then a 16-meter observatory on the Moon, a 4- or 6-meter one in an elliptical heliocentric orbit, and a segmented 8-meter one passively cooled to 50 K at L2, which was finally descoped to 6.6 meters. It also addresses the optimization for scientific performance, the challenge of dealing with such an ultra-low operating temperature, cost issues, supporting technology, modifications made during final design and, finally, how the architecture performs on orbit.

GX 339$-$4 is one of the prototypical black hole X-ray transients, exhibiting recurrent outbursts that allow detailed studies of black hole accretion and ejection phenomena. In this work we present four epochs of optical and near-infrared spectroscopy obtained with X-shooter at the Very Large Telescope. The dataset includes two hard state spectra, collected during the 2013 and 2015 outbursts, and two soft state spectra observed during the 2021 outburst. Strong Balmer, Paschen, He I and He II emission lines are consistently observed in all spectra, while Brackett transitions and the Bowen blend are only prominent in the soft state. Although P-Cygni profiles are not identified, the presence of wind signatures, such as extended emission wings, flat-top and asymmetric red-skewed profiles, is consistently observed through most emission lines, suggesting the presence of wind-type ejecta. These features are particularly evident in the hard state, but they are also observed in the soft state, especially in the near-infrared. This strengthens the case for state-independent winds in black hole transients and increases the evidence for wind signatures in low-to-intermediate orbital inclination systems. We also study the spectral energy distribution, which provides evidence for the presence of synchrotron emission during the hard state. The jet significantly affects the near-infrared continuum, greatly diluting the emission features produced in the accretion flow. The simultaneous identification of both jet and wind signatures during the hard state reinforces the idea of a complex outflow scenario, in which different types of ejecta coexist.

Recently, a double-lined binary (SB2) classical Cepheid, OGLE-LMC-CEP-1347, was discovered, with the orbital period (P$_{\rm orb} = 59$ days) five times shorter than of any binary Cepheid known before. The expected mass of the Cepheid was below $3.5$ M$_\odot$, which, if confirmed, would also probe the uncharted territory. The system configuration also pointed to the Cepheid as a merger. We present a novel method for determining precise physical parameters of binary Cepheids using both theory and observations. This q-PED method combines the measured mass ratio (q), pulsation (P), and evolutionary (E) models, and the known distance (D) supplemented with multi-band photometry. Applying it, we determined the mass of the Cepheid of $3.42 \pm 0.09$ M$_\odot$, its radius of $13.65 \pm 0.27$ R$_\odot$, the companion mass of $1.89 \pm 0.04$ M$_\odot$ and radius of $12.5 \pm 0.62$ R$_\odot$. With the current configuration, the apparent evolutionary age difference of almost $1$ Gyr between the components strongly favors the Cepheid merger origin scenario. If so, the actual age of the Cepheid would be $1.1$ Gyr, on the edge of Population II stars, indicating a significant fraction of Cepheids may be much older than typically assumed. We also applied our method to an eclipsing binary Cepheid with accurately determined physical parameters, obtaining a close agreement, which confirmed our method's reliability.

We conduct a study of the gas kinematics of two quasar host galaxies at $z\gtrsim6$ traced by the [CII] emission line using ALMA. By combining deep observations at both low and high resolution, we recover the diffuse emission, resolve its structure, and measure the rotation curves from the inner region of the galaxy to its outskirts using DysmalPy and 3DBarolo. Assuming that both galaxies exhibit disk rotation driven by the gravitational potential of the galaxy, we find that the best-fit disk models have a $V_{\rm rot}/\sigma \approx 2$ and inferred circular velocities out to $\sim$6-8 kpc scales, well beyond the likely stellar distribution. We then determine the mass profiles of each component (stars, gas, dark matter) with priors on the baryon and dark matter properties. We find relatively large dark matter fractions within their effective radii ($f_{\rm DM}(R<R_e)$ = $0.61_{-0.08}^{+0.08}$ and $0.53_{-0.23}^{+0.21}$, respectively), which are significantly larger than those extrapolated from lower redshift studies and remain robust under different input parameters verified by Monte-Carlo simulations. The large $f_{\rm DM}(R<R_e)$ corresponds to halo masses of $\sim 10^{12.5}-10^{12.8}\, M_\odot$, thus representative of the most massive halos at these redshifts. Notably, while the masses of these SMBHs are approximately 1 dex higher than the low-redshift relationship with stellar mass, the closer alignment of SMBH and halo masses with a local relationship may indicate that the early formation of these SMBHs is linked to their dark matter halos, providing insights into the co-evolution of galaxies and black holes in the early universe.

Observations of the circumgalactic medium (CGM) often display coincident absorption from species with widely varying ionization states, providing direct evidence for complex, multiphase interactions. Motivated by these measurements, we perform a series of cloud-crushing simulations that model cold clouds traveling through the hot CGM. We analyze the ion distributions of these clouds, generate mock absorption spectra, and study their implications on quasar (QSO) absorption observations. Our results show interesting multiphase features, in which ions with significantly different ionization potentials exist in the same absorber and share similar spectral features. However, our simulations are unable to explain high ions like O \textsc{vi} and their coexistence with lower ions that appear in many observed QSO absorption systems.

Recent years have seen a growing sample of TeV emission detections in gamma-ray burst afterglows, as well as an increasing role for structured jets in afterglow modelling. Using a kinetic approach, we show that the structure of an afterglow jet impacts its TeV emission, with jets where the energy falls off more sharply with angle showing a decrease in Inverse Compton (IC) peak flux relative to synchrotron peak flux. We use a modified version of the code katu, to which we have added adiabatic expansion and a fully self-consistent treatment of IC cooling both for the electron and photon populations. We compare our results to the semi-analytical code afterglowpy, finding a good agreement with our model except at early times off axis where the effects of baryon loading are important. We compare electron cooling in the cases where there is no IC cooling, Thomson cooling and an inclusion of Klein-Nishina effects, finding that the spectra can only be distinguished if the Compton potential is significantly increased. The smooth and gradual transition of our KN cooling also leads to a disparity in the cooling between our results and semi-analytical solutions based on asymptotic limits. Finally, we combine best fit parameters from afterglowpy to reproduce the light curves of GRB 170817A in our model. For our choice of parameters, we find that GRB 170817 would not have been detected in the TeV domain if seen on-axis, even by the upcoming Cherenkov Telescope Array Observatory.

We confirm the planetary nature of a pair of transiting sub-Neptune exoplanets orbiting the bright F-type sub-giant star TOI-6054 ($V=8.02$, $K=6.673$) as a part of the OrCAS radial velocity survey using WIYN/NEID observations. We find that TOI-6054b and TOI-6054c have radii of $2.65 \pm 0.15$ $R_{\oplus}$ and $2.81 \pm 0.18$ $R_{\oplus}$, respectively, and masses of $12.4 \pm 1.7$ $M_{\oplus}$ and $9.2 \pm 2.0$ $M_{\oplus}$. The planets have zero-albedo equilibrium temperatures of $1360 \pm 33$ K and $1144 \pm 28$ K. The host star has expanded and will evolve off of the Main Sequence within the next $\sim$500 Myr, and the resulting increase in stellar luminosity has more than doubled the stellar flux the two planets receive compared to the start of the host star's main sequence phase. Consequently, TOI-6054b may be losing some of its primordial H/He atmosphere -- if it has one. Based on dynamical simulations performed using the orbital parameters of the two planets, TOI-6054b, and TOI-6054c are very likely in a 5:3 mean motion resonance. The TOI-6054 system thus has the potential to be an excellent candidate for future atmospheric follow-up observations, with two similarly sized sub-Neptunes around a bright star. We also estimate that if TOI-6054b is currently losing its H/He atmosphere this should be observable from space and from the ground.

Observations of massive supermassive black holes (SMBHs) in the early universe challenge existing black hole formation models. We propose that soliton cores in fuzzy dark matter (FDM) offer a potential solution to this timing problem. Our FDM cosmological zoom-in simulations confirm that for a particle mass $m_{\rm FDM}\sim 10^{-22}~{\rm eV}$, solitons are well developed at redshift $z \sim 7$ with masses of $\sim10^9~M_\odot$, comparable to the observed SMBHs. We then demonstrate using hydrodynamic simulations that, compared to cold dark matter, these high-$z$ massive FDM solitons with mass $M_s$ can provide additional gravitational potential to accrete gas and boost the Bondi accretion rate of a growing black hole seed with mass $M_{\rm BH}$ by up to two to four orders of magnitude, in the regime of efficient cooling and negligible radiation pressure. This accretion boosting mechanism is effective for $10^{-22}~{\rm eV} \lesssim m_{\rm FDM} \lesssim 10^{-20}~{\rm eV}$ and potentially beyond as long as $M_s > M_{\rm BH}$. The simulation code GAMER is accessible at this https URL.

This study explores the origins of cosmic rays and their secondary messengers, focusing on the potential role of four BL Lacs W Comae, 1ES 1959+650, PKS 2005-489, and PKS 2155-304 as potential sources of astrophysical neutrinos and gamma rays. We analyzed a single-zone model to understand the interactions between high-energy protons and ambient photons within blazar jets, leading to neutrino production observables and gamma-ray emission. This modeling contextualizes the emissions within multiwavelength observations and evaluates the capabilities of the next-generation Cherenkov Telescope Array Observatory (CTAO) in detecting these emissions. Our estimations suggest that these sources could be effective emitters of CRs, highlighting the need for future multimessenger observations to further investigate and constrain this class of sources.

We present Mantis Shrimp, a multi-survey deep learning model for photometric redshift estimation that fuses ultra-violet (GALEX), optical (PanSTARRS), and infrared (UnWISE) imagery. Machine learning is now an established approach for photometric redshift estimation, with generally acknowledged higher performance in areas with a high density of spectroscopically identified galaxies over template-based methods. Multiple works have shown that image-based convolutional neural networks can outperform tabular-based color/magnitude models. In comparison to tabular models, image models have additional design complexities: it is largely unknown how to fuse inputs from different instruments which have different resolutions or noise properties. The Mantis Shrimp model estimates the conditional density estimate of redshift using cutout images. The density estimates are well calibrated and the point estimates perform well in the distribution of available spectroscopically confirmed galaxies with (bias = 1e-2), scatter (NMAD = 2.44e-2) and catastrophic outlier rate ($\eta$=17.53$\%$). We find that early fusion approaches (e.g., resampling and stacking images from different instruments) match the performance of late fusion approaches (e.g., concatenating latent space representations), so that the design choice ultimately is left to the user. Finally, we study how the models learn to use information across bands, finding evidence that our models successfully incorporates information from all surveys. The applicability of our model to the analysis of large populations of galaxies is limited by the speed of downloading cutouts from external servers; however, our model could be useful in smaller studies such as generating priors over redshift for stellar population synthesis.

Long-period radio transients are a novel class of astronomical objects characterised by prolonged periods ranging from 18 minutes to 54 minutes. They exhibit highly polarised, coherent, beamed radio emission lasting only 10--100 seconds. The intrinsic nature of these objects is subject to speculation, with highly magnetised white dwarfs and neutron stars being the prevailing candidates. Here we present ASKAP J183950.5-075635.0 (hereafter, ASKAP J1839-0756), boasting the longest known period of this class at 6.45 hours. It exhibits emission characteristics of an ordered dipolar magnetic field, with pulsar-like bright main pulses and weaker interpulses offset by about half a period are indicative of an oblique or orthogonal rotator. This phenomenon, observed for the first time in a long-period radio transient, confirms that the radio emission originates from both magnetic poles and that the observed period corresponds to the rotation period. The spectroscopic and polarimetric properties of ASKAP J1839-0756 are consistent with a neutron star origin, and this object is a crucial piece of evidence in our understanding of long-period radio sources and their links to neutron stars.

Cosmological information is encoded in the structure of galaxy clusters. In Universes with less matter and larger initial density perturbations, clusters form earlier and have more time to accrete material, leading to a more extended infall region. Thus, measuring the mean mass distribution in the infall region provides a novel cosmological test. The infall region is largely insensitive to baryonic physics, and provides a cleaner structural test than other measures of cluster assembly time such as concentration. We consider cluster samples from three publicly available galaxy cluster catalogues: the Spectrsopic Identification of eROSITA Sources (SPIDERS) catalogue, the X-ray and Sunyaev-Zeldovich effect selected clusters in the meta-catalogue M2C, and clusters identified in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Survey. Using a preliminary shape catalogue from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), we derive excess surface mass density profiles for each sample. We then compare the mean profile for the DESI Legacy sample, which is the most complete, to predictions from a suite of simulations covering a range of $\Omega_{\rm m}$ and $\sigma_8$, obtaining constraints of $\Omega_{\rm m}=0.29\pm 0.05$ and $\sigma_8=0.80 \pm 0.04$. We also measure mean (comoving) splashback radii for SPIDERS, M2C and DESI Legacy Imaging Survey clusters of $1.59^{+0.16}_{-0.13} {\rm cMpc}/h$, $1.30^{+0.25}_{-0.13} {\rm cMpc}/h$ and $1.45\pm0.11 {\rm cMpc}/h$ respectively. Performing this analysis with the final UNIONS shape catalogue and the full sample of spectroscopically observed clusters in DESI, we can expect to improve on the best current constraints from cluster abundance studies by a factor of 2 or more.

Context. During the primary Kepler mission, between 2009 and 2013, about 150,000 pre-selected targets were observed with a 29.42 minute-long cadence. However, a survey of background stars that fall within the field of view (FOV) of the downloaded apertures of the primary targets has revealed a number of interesting objects. In previous papers we have presented surveys of short period Eclipsing Binaries and RR Lyrae stars. Aims. The current survey of the Kepler background is concentrated on identifying longer-period eclipsing binaries and pulsating stars. These will be the subject of later papers. In the course of this survey, in addition to eclipsing binaries and pulsating stars, seven exoplanet candidates have been uncovered and in this paper we report on these candidates. Methods. We use Lomb-Scargle (LS), light curve transit search and Phase Dispersion Minimisation (PDM) methods to reveal pixels that show significant periodicities, resulting in the identification of the seven exoplanet candidates. We have prepared the light curves for analysis using Pytransit software and cross matched the pixel coordinates with Gaia and other catalogues to identify the sources. this http URL identify seven Hot Jupiter exoplanet candidates with planet radii ranging from 0.8878 to 1.5174 RJup and periods ranging from 2.5089 to 4.7918 days.

We present a new code and approach, JERALD -- JAX Enhanced Resolution Approximate Lagrangian Dynamics -- , that improves on and extends the Lagrangian Deep Learning method of Dai & Seljak (2021), producing high-resolution dark matter, stellar mass and neutral hydrogen maps from lower-resolution approximate $N$-body simulations. The model is trained using the Sherwood-Relics simulation suite (for a fixed cosmology), specifically designed for the intergalactic medium and the neutral hydrogen distribution in the cosmic web. The output is tested in the redshift range from $z=5$ to $z=0$ and the generalization properties of the learned mapping is demonstrated. JERALD produces maps with dark matter, stellar and neutral hydrogen power spectra in excellent agreement with full-hydrodynamic simulations with $8\times$ higher resolution, at large and intermediate scales; in particular, JERALD's neutral hydrogen power spectra agree with their higher-resolution full-hydrodynamic counterparts within 90% up to $k\simeq1\,h$Mpc$^{-1}$ and within 70% up to $k\simeq10\,h$Mpc$^{-1}$. JERALD provides a fast, accurate and physically motivated approach that we plan to embed in a statistical inference pipeline, such as Simulation-Based Inference, to constrain dark matter properties from large- to intermediate-scale structure observables.

The origin of heavy elements synthesized through the rapid neutron capture process ($r$-process) has been an enduring mystery for over half a century. Cehula et al. (2024) recently showed that magnetar giant flares, among the brightest transients ever observed, can shock-heat and eject neutron star crustal material at high velocity, achieving the requisite conditions for an $r$-process. Patel et al. (in prep.) confirmed an $r$-process in these ejecta using detailed nucleosynthesis calculations. Radioactive decay of the freshly synthesized nuclei releases a forest of gamma-ray lines, Doppler broadened by the high ejecta velocities $v \gtrsim 0.1c$ into a quasi-continuous spectrum peaking around 1 MeV. Here, we show that the predicted emission properties (light-curve, fluence, and spectrum) match a previously unexplained hard gamma-ray signal seen in the aftermath of the famous December 2004 giant flare from the magnetar SGR 1806-20. This MeV emission component, rising to peak around 10 minutes after the initial spike before decaying away over the next few hours, is direct observational evidence for the synthesis of $\sim 10^{-6}M_{\odot}$ of $r$-process elements. The discovery of magnetar giant flares as confirmed $r$-process sites, contributing at least $\sim 1$-$10\%$ of the total Galactic abundances, has implications for the Galactic chemical evolution, especially at the earliest epochs probed by low-metallicity stars. It also implicates magnetars as potentially dominant sources of heavy cosmic rays. Characterization of the $r$-process emission from giant flares by resolving decay line features offers a compelling science case for NASA's forthcoming COSI nuclear spectrometer, as well as next-generation MeV telescope missions.

Ultra-soft narrow line Seyfert 1 (US-NLSy) are a poorly observed class of active galactic nuclei characterized by significant flux changes and an extreme soft X-ray excess. This peculiar spectral shape represents a golden opportunity to test whether the standard framework commonly adopted for modelling local AGN is still valid. We thus present the results on the joint XMM-Newton and HST monitoring campaign of the highly accreting US-NLSy RBS 1332. The optical-to-UV spectrum of RBS 1332 exhibits evidence of both a stratified narrow-line region and an ionized outflow, that produces absorption troughs over a wide range of velocities (from ~1500 km s-1 to ~1700 km s-1) in several high-ionization transitions (Lyalpha, N V, C IV). From a spectroscopic point of view, the optical/UV/FUV/X-rays emission of this source is due to the superposition of three distinct components which are best modelled in the context of the two-coronae framework in which the radiation of RBS 1332 can be ascribed to a standard outer disk, a warm Comptonization region and a soft coronal continuum. The present dataset is not compatible with a pure relativistic reflection scenario. Finally, the adoption of the novel model reXcor allowed us to determine that the soft X-ray excess in RBS 1332 is dominated by the emission of the optically thick and warm Comptonizing medium, and only marginal contribution is expected from relativistic reflection from a lamppost-like corona.

Fast radio bursts (FRBs) are millisecond-duration pulses occurring at cosmological distances with a mysterious origin. Observations show that at least some FRBs are produced by magnetars. All magnetar-powered FRB models require some triggering mechanisms, among which the most popular is the crust cracking of a neutron star, which is called starquake. However, so far there has been no decisive evidence for this speculation. Here we report the energy functions of the three most active repeating FRBs, which show a universal break around $10^{38}$ erg. Such a break is similar to that of the frequency-magnitude relationship of earthquakes. The break and change of the power-law indices below and above it can be well understood within the framework of FRBs triggered by starquakes in the magnetar models. The seed of weak FRBs can grow both on the magnetar surface and in the deeper crust. In contrast, the triggering of strong FRBs is confined by the crustal thickness and the seed of strong FRBs can only grow on the surface. This difference in dimensionality causes a break in the scaling properties from weak to strong FRBs, occurring at a point where the penetration depth of starquakes equals the crustal thickness. Our result, together with the earthquake-like temporal properties of these FRBs, strongly supports that FRBs are triggered by starquakes, providing a new opportunity to study the physical properties of the neutron star crust.

Recent baryon acoustic oscillation (BAO) measurements from the Dark Energy Spectroscopic Instrument (DESI) collaboration, combined with the cosmic microwave background (CMB) and type Ia supernovae (SNe Ia) observations, suggest a preference for dynamical dark energy (DDE) with $w_0>-1$ and $w_a<0$. Given the cosmological origin of fast radio bursts (FRBs), the combination of their dispersion measures and host galaxy redshifts makes localized FRBs a valuable tool for probing cosmology. Using an updated sample of 92 localized FRBs, along with DESI BAO, PlantheonPlus and CMB data, we constrain the dark energy (DE) equation of state (EoS) under the Chevallier-Polarski-Linder (CPL) parameterization. We find that even without incorporating CMB data, DDE remains preferred with $w_0 = -0.855 ^{+0.084}_{-0.084}$ and $w_a = -1.174^{+0.462}_{-0.491}$ at a confidence level of $\sim2.5 \sigma$. A joint analysis constrains these to be $w_0 = -0.784^{+0.064}_{-0.064}$ and $w_a = -0.872^{+0.269}_{-0.278}$, showing a discrepancy with $\Lambda$CDM at a $\sim3.1\sigma$ level. Furthermore, using localized FRBs alone, we estimate the Hubble constant $H_0$ to be $69.04^{+2.30}_{-2.07}$ and $75.61^{+2.23}_{-2.07} \, \rm km \, s^{-1} \, Mpc^{-1}$, assuming the Galactic electron density models to be NE2001 (Cordes \& Lazio) and YMW16 (Yao et al.), respectively. Thus, accurate accounting of the Galactic dispersion measure is crucial for resolving the Hubble tension with FRBs. Future BAO measurements, next-generation CMB experiments, and more localized FRBs will further constrain the DE EoS and the cosmological parameters.

We report a search for high-energy astrophysical neutrino multiplets, detections of multiple neutrino clusters in the same direction within 30 days, based on an analysis of 11.4 years of IceCube data. A new search method optimized for transient neutrino emission with a monthly time scale is employed, providing a higher sensitivity to neutrino fluxes. This result is sensitive to neutrino transient emission, reaching per-flavor flux of approximately $10^{-10}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm sec}^{-1}$ from the Northern sky in the energy range $E\gtrsim 50$~TeV. The number of doublets and triplets identified in this search is compatible with the atmospheric background hypothesis, which leads us to set limits on the nature of neutrino transient sources with emission timescales of one month.

We present an extensive temporal and spectral study of the Seyfert 1 AGN Mrk 50 using 15 years (2007-2022) of multiwavelength observations from XMM-Newton, Swift, and NuSTAR for the first time. From the timing analysis, we found that the source exhibited variability of $\sim$20 % during the 2007 observation, which reduced to below 10 % in the subsequent observations and became non-variable in the observations from 2010 onward. From the spectral study, we found that the spectra are nearly featureless. Non-detection of absorption in the low-energy domain during the 15 years of observation infers the absence of obscuration around the central engine, rendering the nucleus a `bare' type. A prominent soft X-ray excess below 2 keV was detected in the source spectrum during the observations between 2007 and 2010, which vanished during the later observations. To describe the nature of the soft excess, we use two physical models, such as warm Comptonization and blurred reflection from the ionized accretion disk. Both the physical models explain the nature and origin of the soft excess in this source. Our analysis found that Mrk~50 accretes at sub-Eddington accretion rate ($\lambda_{Edd}=0.13-0.02$) during all the observations used in this work.

We study the nucleosynthesis in a core-collapse supernova model including newly calculated neutrino-induced reaction rates with both collective and Mikheyev-Smirnov-Wolfenstein (MSW) neutrino-flavor oscillations considered. We show that the measurement of a pair of $^{11}$B/$^{10}$B and $^{138}$La/$^{139}$La or $^6$Li/$^7$Li and $^{138}$La/$^{139}$La in presolar grains that are inferred to have originated from core-collapse supernovae could constrain the neutrino mass hierarchy. The new shell-model and the model of quasi-particle random phase approximation in the estimate of three important neutrino-induced reactions, $\nu+^{16}$O, $\nu+^{20}$Ne, and $\nu+^{138}$Ba are applied in our reaction network. The new rates decrease the calculated $^{7}$Li/$^{6}$Li ratio by a factor of five compared with the previous study. More interestingly, these new rates result in a clear separation of the isotopic ratio of $^{11}$B/$^{10}$B between normal and inverted mass hierarchies in the O/Ne, O/C, and C/He layers where $^{138}$La abundance depends strongly on the mass hierarchy. In these layers, the sensitivity of the calculated abundances of $^{10,11}$B and $^{6,7}$Li to the nuclear reaction uncertainties is also tiny. Therefore, we propose that the $^{11}$B/$^{10}$B vs. $^{138}$La/$^{139}$La and $^6$Li/$^7$Li vs. $^{138}$La/$^{139}$La in type X silicon carbide grains sampled material from C/He layer can be used as a new probe to constrain the neutrino mass hierarchy.

We present the results from long-term simultaneous monitoring observations of SiO and H2O masers toward the Mira variable star WX Serpentis. This study has been conducted with 21m single-dish radio telescopes of the Korean VLBI Network from 2009 June to 2021 June. Five maser lines were considered: SiO v=1, 2, J=1-0; SiO v=1, J=2-1, 3-2, and H2O 6(1,6)-5(2,3) transitions, with the SiO maser lines distributed near the stellar velocity and the H2O maser exhibiting an asymmetric line profile with five to six peaked components. Intense H2O maser emissions suddenly appeared in 2019 September, indicating flaring. The intensity variations of SiO and H2O masers are strongly correlated with the optical light curve (OLC) of the central star, with individual phase lags; the phase lag of the H2O maser relative to the OLC is larger than that of the SiO masers. The consequent phase difference between the SiO masers and the H2O maser likely indicates that their formation regions and main driving mechanisms are different from each other. The SiO masers in WX Ser exhibit a dominant single-peak velocity distribution, similar to other Mira variable stars. However, the H2O maser displays distinct morphological features, showing a radial acceleration and preferential intensity dominance at blueshifted velocities. This suggests that the H2O maser clouds of WX Ser are moving outward, thereby developing an asymmetric outflow owing to nonuniform material ejection from the stellar atmosphere. The findings confirm that an initial asymmetric outflow structure emerged during the thermally pulsing asymptotic giant branch phase, specifically in the Mira variable star stage.

Turbulent small-scale structures in the envelopes and winds of massive stars have long been suggested as the cause for excessive line broadening that could not be explained by other mechanisms such as thermal broadening. However, the origin of these structures, particularly in the envelope, has not been extensively studied. We study the origin of structures seen in 2D unified stellar atmosphere and wind simulations of O stars and Wolf-Rayet (WR) stars. Particularly, we study whether the structure growth in the simulations is consistent with sub-surface convection, as is commonly assumed to be the origin of this turbulence. Using a linear stability analysis of the optically thick envelopes of massive stars, we identified multiple instabilities that could drive structure growth. We quantified the structure growth in the non-linear simulations of O stars and WR stars by computing density power spectra and tracking their temporal evolution. Then, we compared these results to the analytical results from the stability analysis. The stability analysis leads to two possible instabilities: the convective instability and an acoustic instability. Analytic expressions for the growth rates of these different instabilities are found. In particular, strong radiative diffusion damps the growth rate $\omega$ of the convective instability leading to a distinct $\omega \sim 1/k^2$ dependence on wavenumber $k$. From our power spectra analysis of the simulations, however, we find that structure growth rather increases with $k$ - tentatively as $\omega \sim \sqrt{k}$. Our results suggest that structures in luminous O and WR star envelopes do not primarily develop from the sub-surface convective instability. Rather the growth seems compatible with either the acoustic instability in the radiation-dominated regime or with Rayleigh-Taylor type instabilities, although the exact origin remains inconclusive for now.

We assess the detectability of tidal disruption events (TDEs) using mock observations from the Mini SiTian array. We select 100 host galaxy samples from a simulated galaxy catalog based on specific criteria such as redshift, BH mass, and event rate. Taking into account the site conditions and survey strategy, we simulate observations over a 440 deg$^2$ field. The results indicate that $0.53\pm 0.73$ TDEs can be detected per year when observing in both $g$ and $r$ bands with 300-second exposures every 3 days. Applying this method to the SiTian project, we expect to discover approximately 204 TDEs annually, heralding a new era in TDE science.

We present a numerical proof of the concept that the void spin distributions can in principle provide a tight constraint on the amplitude of matter density fluctuation on the scale of $8\,h^{-1}{\rm Mpc}$ ($\sigma_{8}$) without being severely deteriorated by the degeneracies of $\sigma_{8}$ with cold dark matter density parameter multiplied by the dimensionless Hubble parameter square ($\Omega_{\rm cdm}h^{2}$), total neutrino mass ($M_{\nu}$) and dark energy equation of state ($w$). Applying the Void-Finder algorithm to a total of $15$ AbacusSummit $N$-body simulations of $15$ different cosmological models, we identify the giant voids to measure their spins, defined as the magnitudes of rescaled specific angular momenta of void halos. The $15$ cosmologies include the Planck $\Lambda$CDM and $14$ non-Planck models, each of which differs among one another only in one of $\{\sigma_{8},\ \Omega_{\rm cdm}h^{2},\ M_{\nu},\ w\}$. The probability density distribution of void spins is determined for each model and found to be well approximated by the generalized Gamma distribution with two characteristic parameters, $k$ and $\theta$. It turns out that the best-fit values of $k$ and $\theta$ exhibit very sensitive dependences only on $\sigma_{8}$, being almost insensitive to $\Omega_{\rm cdm}h^{2}$, $M_{\nu}$, $w$. This exclusive $\sigma_{8}$-dependence of the void spin distributions is confirmed to be robust against the variation of the mass and number cuts of void halos. We also test an observational feasibility of estimating the void spins from real data on the galaxy redshifts.

Medium-timescale (minutes to hours) radio transients are a relatively unexplored population. The wide field-of-view and high instantaneous sensitivity of instruments such as MeerKAT provides an opportunity to probe this class of sources, using image-plane detection techniques. We aim to systematically mine archival synthesis imaging data in order to search for medium-timescale transients and variables that are not detected by conventional long-track image synthesis techniques. We deploy a prototype blind transient and variable search pipeline named TRON. This processes calibrated visibility data, constructs high-time cadence images, performs a search for variability on multiple timescales, and extracts lightcurves for detected sources. As proof of concept, we apply it to three MeerKAT observations of globular clusters, known to host transient or variable sources. We detect a previously known eclipsing MSP suspected to be a `black widow' system, in the globular cluster Omega Centauri, with a light curve confirming the eclipsing nature of the emission. We detect a previously known `red back' eclipsing MSP in the globular cluster Terzan 5. Using observations of the globular cluster 47 Tucanae, we detect two known millisecond pulsars (MSPs), and one previously reported MSP candidate, with hints of eclipsing behaviour.

Medium-timescale (minutes to hours) radio transients are a relatively unexplored population. The wide field-of-view and high instantaneous sensitivity of instruments such as MeerKAT provides an opportunity to probe this class of sources, using image-plane detection techniques. The previous letter in this series describes our project and associated TRON pipeline designed to mine archival MeerKAT data for transient and variable sources. In this letter, we report on a new transient, a radio flare, associated with Gaia DR3 6865945581361480448, a G type star, whose parallax places it at a distance of 1334 pc. Its duration and high degree of circular polarization suggests electron cyclotron maser instability as the mechanism, consistent with an RS CVn variable.

The dynamics of the origin of gamma-ray emissions in gamma-ray bursts (GRBs) remains an enigma. Through a joint analysis of GRB 180427A, observed by the Fermi Gamma-ray Space Telescope and AstroSat's Cadmium Zinc Telluride Imager, we identify emissions from two distinct regions with varying polarisation properties. Time-resolved polarisation analysis reveals a synchronous evolution of the polarisation angle (PA) and fraction (PF) with two emission pulses, peaking with a delay of $4.82 \pm 0.12\, \mathrm{s}$. Spectral analysis indicates that the first pulse is dominated by blackbody radiation, while the second pulse exhibits a non-thermal spectrum (power law with an exponential cutoff). Using a bottom-to-top approach through simulations, we decouple the polarisation properties of the individual spectral components, revealing polarisation fractions of 25\% - 40\% for the blackbody spectrum and 30\% - 60\% for the non-thermal spectrum. At a redshift of $z \sim 0.05$, the blackbody emission originates from the jet photosphere at $10^{11}\, \mathrm{cm}$, whereas the non-thermal emission arises from an optically thin region at $10^{15}\, \mathrm{cm}$. The changing dominance of these emissions explains the observed PA shift of $60^\circ \pm 22.3^\circ$. The spectral cutoff at 1 MeV suggests pair opacity due to the jet's low bulk Lorentz factor ($\Gamma \sim$ tens). The high polarisation and hard spectral slopes ($\alpha > -0.5$) imply a top-hat jet structure observed off-axis, near the jet's edge. This off-axis viewing introduces anisotropy in the radiation within the viewing cone ($1/\Gamma$), accounting for the observed polarisation.

The classical picture of our Solar System being the archetypal outcome of planet formation has been rendered obsolete by the astonishing diversity of extrasolar-system architectures. From rare hot-Jupiters to abundant super-Earths and sub-Neptunes, most detected exoplanets have no analogs in our system, and their interior and atmospheric compositions remain largely unknown. Fortunately, new methodologies enable us to analyze exoplanet atmospheres, inferring their compositions, temperatures, dynamics, and even formation pathways. Specifically, ground-based high-resolution Doppler spectroscopy (HRDS) can disentangle spectral-line profiles of weak exo-atmospheric signals from the dominating features of Earth's atmosphere in the observed flux. For over a decade, HRDS has focused on hot Jupiters (close-orbiting gas giants) due to their high signal-to-noise ratio, which makes them ideal laboratories for advancing our knowledge. However, there have been concerns regarding potential biases in exo-atmospheric-detection methods, hindering comparative planetology. Here we propose a modeling framework based on extensive simulations of HRDS exo-atmospheric observations to systematically explore in-silico underlying biases in commonly-used pipelines, particularly under the presence of observational noise. Our findings show that exo-atmospheric detection-significances are highly contingent on details of the analysis-pipeline used, with different techniques responding differently to noise: A given technique may fail to recover a true-signal that is detected by another. Noise effects in the computed significances are non-trivial and pipeline-dependent. Statistical analyses provide a complementary tool to contextualize signal-significances, which will gain in relevance as we move towards studying the atmospheres of smaller, potentially habitable exoplanets with even weaker signals.

The origin of the correlation between the effective spins ($\chi_{\rm eff}$) and mass ratios ($q$) of LIGO-Virgo-KAGRA's binary black holes (BBHs) is still an open question. Motivated by recent identification of two subpopulations of the BBHs, in this work we investigate the potential $\chi_{\rm eff}-q$ correlation for each subpopulation. Surprisingly, the $\chi_{\rm eff}$-$q$ correlation vanishes for the low-mass subpopulation if we introduce a second $\chi_{\rm eff}$ distribution for the high-mass subpopulation likely originating from hierarchical mergers. The first subpopulation has a narrow $\chi_{\rm eff}$ distribution peaking at $\sim0.05$, whose primary-mass function cuts off at $\sim 45M_{\odot}$, in agreement with first-generation BBHs. The second $\chi_{\rm eff}$ distribution is broad and peaks at $\mu_{\chi,2}=0.35^{+0.18}_{-0.22}$, consistent with the expectation of hierarchical mergers formed in the disks of active galactic nucleus (AGNs). We infer $\mu_{\chi,2}>0$ at 98.7\% credible level, and a symmetric $\chi_{\rm eff}$ distribution for the second subpopulation is disfavored by $\mathcal{B}\sim5$. However, negative values of $\chi_{\rm eff}$ are also measured, indicating that the hierarchical mergers may take place within both star clusters and AGN disks. We find a Bayes factor of $\ln\mathcal{B}=5.2$ for two distinct $\chi_{\rm eff}$ distributions relative to single $\chi_{\rm eff}$ distribution that conditioned on mass ratios. Therefore we conclude that the $\chi_{\rm eff}$-$q$ correlation in the entire population can be explained as the superposition of two subpopulations. Additionally, we suggest to use a flexible mass function to reduce the bias in $\chi_{\rm eff}$-$q$ correlation that may be introduced by model mis-specification.

Ultraluminous X-ray sources (ULXs) with neutron star (NS) accretors challenge traditional accretion models, and have sparked a debate regarding the role of geometrical beaming and strong magnetic fields (B). The reduction of the Thomson cross-section in the presence of strong B, leads to a modification of the Eddington limit, and therefore is expected to affect significantly the observational appearance of NS-ULXs. We investigate the role of this modification using population synthesis models, and explore its effects on the X-ray luminosity functions, spin-up rates, and outflow energetics of the observed NS-ULXs. Our results show that the new prescription allows NS-ULXs to achieve super-Eddington luminosities with milder beaming compared to before, improving the agreement with observations. In addition, it broadens the range of spin-up rates allowing for more diverse conditions in NS-ULXs in terms of accretion rates and magnetic fields. More importantly, the reduced beaming increases the likelihood of observing the NS-ULXs within wind-powered nebulae such as NGC 5907 ULX-1. Our findings highlight the necessity of taking into account B effects independently of the approach: geometrical beaming or strong B, and call for magnetospheric accretion prescriptions that can be integrated in population synthesis codes.

An avenue for understanding cosmological galaxy formation is to compare morphometric parameters in observations and simulations of galaxy assembly. In this second paper of the ASymba: Asymmetries of HI in SIMBA Galaxies series, we measure atomic gas HI asymmetries in spatially-resolved detections from the untargetted WALLABY survey, and compare them to realizations of WALLABY-like mock samples from the SIMBA cosmological simulations. We develop a Scanline Tracing method to create mock galaxy HI datacubes which minimizes shot noise along the spectral dimension compared to particle-based methods, and therefore spurious asymmetry contributions. We compute 1D and 3D asymmetries for spatially-resolved WALLABY Pilot Survey detections, and find that the highest 3D asymmetries A3D>0.5 stem from interacting systems or detections with strong bridges or tails. We then construct a series of WALLABY-like mock realizations drawn from the SIMBA 50 Mpc simulation volume, and compare their asymmetry distributions. We find that the incidence of high A3D detections is higher in WALLABY than in the SIMBA mocks, but that difference is not statistically significant (p-value = 0.05). The statistical power of quantitative comparisons of asymmetries such as the one presented here will improve as the WALLABY survey progresses, and as simulation volumes and resolutions increase.

Multiple populations are ubiquitous in the old massive globular clusters (GCs) of the Milky Way. It is still unclear how they arose during the formation of a GC. The topic of iron and metallicity variations has recently attracted attention with the measurement of iron variations among the primordial population (P1) stars of Galactic GCs. We use the spectra of more than 8000 RGB stars in 21 Galactic GCs observed with MUSE to derive individual stellar metallicities [M/H]. For each cluster, we use the HST photometric catalogs to separate the stars into two main populations (P1 and P2). We measure the metallicity spread within the primordial population of each cluster by combining our metallicity measurements with the stars $\Delta_{\rm F275W,F814W}$ pseudo-color. We also derive metallicity dispersions ($\sigma_{\rm [M/H]}$) for the P1 and P2 stars of each GC. In all but three GCs, we measure a significant correlation between the metallicity and the $\Delta_{\rm F275W,F814W}$ pseudo-color of the P1 stars such that stars with larger $\Delta_{\rm F275W,F814W}$ have higher metallicities. We measure metallicity spreads that range from 0.03 to 0.24 dex and correlate with the GC masses. As for the intrinsic metallicity dispersions, when combining the P1 and P2 stars, we measure values ranging from 0.02 dex to 0.08 dex that correlate very well with the GC masses. We compared the metallicity dispersion among the P1 and P2 stars and found that the P2 stars have metallicity dispersions that are smaller or equal to that of the P1 stars. We find that both the metallicity spreads of the P1 stars (from the $\Delta_{\rm F275W,F814W}$ spread in the chromosome maps) and the metallicity dispersions ($\sigma_{\rm [M/H]}$) correlate with the GC masses, as predicted by some theoretical self-enrichment models presented in the literature.

Galaxy clusters are currently the endpoint of the hierarchical structure formation; they form via the accretion of dark matter and cosmic gas from their local environment. In particular, filaments contribute grandly by accreting gas from cosmic matter sheets and underdense regions and feeding it to the galaxy clusters. Along the way, the gas in filaments is shocked and heated, which, together with the velocity structure within the filament, induces swirling and, thus, turbulence. In this work, we study a constrained hydrodynamical simulation replica of the Virgo cluster to characterise the velocity field in the two cosmic filaments connected to the cluster with unprecedented high resolution. First, we conduct a qualitative examination of slices extracted from the simulation. We study the temperature, the velocity field, and derived quantities in longitudinal cuts to study the general structure of the filaments and in transverse cuts to study their inner organisation and connection to cosmic matter sheets and underdense regions. Then, we conduct a quantitative study of velocities in Virgo's filaments by computing the 2D energy spectrum from 1 and 5~Mpc square maps extracted from the slices and centred on the core of the filaments. We show that the velocity field goes from mostly compressive far in the filaments to mostly solenoidal in Virgo's core. Moreover, we observe that the total energy spectrum in the filaments gains in amplitude and steepens towards Virgo.

We explore spectroscopic and photometric methods for identifying high-redshift galaxies containing an Active Galactic Nucleus (AGN) with JWST observations. After demonstrating the limitations of standard optical methods, which appear ineffective in the low-metallicity environment of the early universe, we evaluate alternative diagnostic techniques using the current JWST observational capabilities. Our analysis focuses on line ratios and equivalent widths (EWs) of UV emission lines: CIV, HeII $\lambda$1640, OIII] $\lambda$1665, and CIII], and the faint optical line, HeII $\lambda$4686. We find that the most valuable diagnostic quantities for finding AGN are the line ratios: (CIII] + CIV) / HeII $\lambda$1640 and CIII] / HeII $\lambda$1640, as well as the EW of HeII $\lambda$1640. For more reliable AGN identification, the HeII $\lambda$1640 and OIII] $\lambda$1665 lines would need to be detected separately. We show that the HeII $\lambda$1640/H$\beta$ ratio effectively separates AGN from star-forming galaxies, though it is contingent on a low dust content. We also show that in order to effectively use these diagnostics, future observations require longer exposure times, especially for galaxies at $z > 6$. Subsequently, we plot three real high-redshift sources on these diagrams which present strong UV emission lines. However, in order to classify them as strong AGN candidates, further study is needed due to the blending of HeII + OIII] and unreliable optical lines. Lastly, we carry out a selection process using spectral energy distribution (SED) fitting with EAZY to identify strong AGN candidates in the JADES NIRCam photometry. One galaxy in our sample emerged as a strong AGN candidate, supported by both photometric selection and strong UV emission. We present a sample of similar AGN candidates in the JADES data based on this method.

The white dwarf mass distribution has been studied primarily at two extremes: objects that presumably evolved as single stars and members of close binaries that likely underwent substantial interaction. This work considers the intermediate separation regime of ~1 au and demonstrates how binary interaction affects white dwarf masses. The binary mass ratio distribution is utilized for this purpose. Modeled as a truncated Pareto profile, this distribution provides insights into the populations' properties and evolutionary history. When applied to homogeneous samples of binaries with giant primaries of similar age, the distribution's shape constrains the fraction of white dwarf companions, the white dwarf mass distribution, and the properties of their progenitors. As a test case, this method is applied to a small spectroscopic sample of binaries in open clusters with red giant primaries and orbital periods between 0.5 and 20 years. The analysis reveals that white dwarfs in these systems are ~20% less massive than their isolated counterparts, with a typical mass of ~0.55 Msun. Their progenitors likely lost 80-85% of their mass, with binary interactions enhancing mass loss by an additional ~0.2 Msun. These findings highlight the utility of this approach for studying binary evolution and improving population models, particularly with future datasets from Gaia and other large-scale surveys.

Gamma-ray bursts (GRBs) are widely suggested as potential sources of ultrahigh-energy cosmic rays (UHECRs). The kinetic energy of the jets dissipates, leading to the production of an enormous amount of $\gamma$-ray photons and possibly also the acceleration of protons. The accelerated protons will interact with the radiation of the GRB via the photomeson and Bethe-Heitler processes, which can initiate electromagnetic cascades. This process can give rise to broadband radiation up to the GeV-TeV $\gamma$-ray regime. The expected $\gamma$-ray flux from cascades depends on properties of the GRB jet, such as the dissipation radius $R_{\rm diss}$, the bulk Lorentz factor $\Gamma$, and the baryon loading factor $\eta_p$. Therefore, observations of Fermi-LAT can impose constraints on these important parameters. In this study, we select 12 GRBs of high keV-MeV fluence and constrain the baryon loading factor, under different combinations of the bulk Lorentz factor and the dissipation radius based on Fermi-LAT's measurements. Our findings indicate a strong constraint of $\eta_p<10$ for most selected GRBs over a large parameter space except for large dissipation radii ($\gtrsim 10^{15}\rm cm$) and high bulk Lorentz factors ($\gtrsim 600$). The constraint is comparable to, and in some GRBs even stronger than, that from high-energy neutrinos for stacked GRBs. Our results suggest that for typical bulk Lorentz factor of several hundreds, the dissipation radii of GRBs need be large to avoid overshooting the GeV gamma-ray flux during the prompt emission phase of GRBs, which can be used to constrain GRBs.

Changing-look active galactic nuclei (CLAGNs) are known to change their spectral type between 1 and 2 (changing-state) or change their absorption between Compton-thick and Compton-thin (changing-obscuration) on timescales of years or less. The physical mechanism and possible connection between the two types of CLAGNs are still unclear. We explore the evolution of the broadband X-ray spectra from Nuclear Spectroscopic Telescope Array (\nustar\,) and column density in five CLAGNs with moderate inclination viewing angles, which have shown significant variations of both optical types and X-ray absorption. Based on a phenomenological and two clumpy torus models, we find that the X-ray photon index ($\Gamma$) and the Eddington-scaled X-ray $2-10$ keV luminosity ($L_{\rm X}/L_{\rm Edd}$) are positively correlated for the five sources, which are similar to other bright AGNs and optical CLAGNs at type 1 phase. We find a significant negative correlation between log$N_\mathrm{H,los}$ and log$L_{\rm X}/L_{\rm Edd}$ except for ESO 362-G18. Similar to changing-state AGNs, changing-obscuration AGNs may be also triggered by the evolution of the accretion disc. Our results support the disc wind scenario, where the disc wind proportional to the accretion rate and formed at moderate inclination angles would push the obscuration material further away and decrease the column density from the line of sight observed in the changing-look AGNs.

The Chern-Simons gravitational term during inflation is usually coupled to the inflaton field. The resulting theory suffers from ghost-field formation in the tensor sector, which limits the observational effects of P-violation on cosmological correlators. In this work, we consider the Chern-Simons term coupled to an isocurvature component in a multi-field model of inflation. Since the resulting theory does not affect the quadratic action of tensor perturbations, ghost fields do not appear. This operator provides (P-violating) interactions between the isocurvature perturbation and the curvature and tensor perturbations. We show that combining these couplings with interactions between the curvature and isocurvature components coming from a turning trajectory, the resulting $\langle sst \rangle_{PV}$ non-Gaussianities can reach $f^{sst, PV}_{\rm NL}=B_{PV}^{\zeta\zeta h}(k,k,k)/P^2_{\zeta}(k)\sim \mathcal O(1)$ within the parameter space of the theory. Our result motivates the systematic study of the Chern-Simons gravitational term coupled to isocurvature fields in multi-field models of inflation with couplings between the curvature and isocurvature fields or other mechanisms that transfer effects on the isocurvature field into the curvature field.

Tidal locking of planets to their host stars results in an atmospheric circulation with a hotspot fixed to the frame of reference of the planet. On the other hand, asynchronously rotating planets feature moving hotspots either lagging or leading the corresponding substellar point as it translates along the surface. We show that a planet falling in the latter category could mimic the circulation of tidally synchronous planets under the influence of time-varying instellation, possibly provided by pulsating or multiple star systems. This happens when the planets diurnal period is in resonance with the period of instellation variation, leading to a planet-frame-fixed hotspot. Slight differences in the above periods lead to East-West or West-East creeping hotspots with a period significantly longer than both. The rate of hotspot motion is given by the difference between the diurnal and instellation variation rates, similar to the lower envelope frequency of beat patterns formed by two superposed waves in linear wave theory. We call this phenomenon beating. A combination of the radiative, rotational, wave propagation, and drag timescales establishes dynamical constraints on beating. Based on this we classify a set of Kepler and TESS circumbinary planets with two candidates exhibiting climatic departures from the no-variation scenario. In general, hotter and faster-spinning planets are more susceptible to climatic departures. Beating, if it occurs, may additionally create optimistic extensions of habitable zones for corresponding systems.

We introduce the Serendipitous H-ATLAS-fields Observations of Radio Extragalactic Sources (SHORES) multiple pencil beam survey that observed at 2.1 GHz with the Australia Telescope Compact Array (ATCA) 29 fields in total intensity and polarization within the Herschel-ATLAS Southern Galactic Field. This paper presents the observations, calibration and analysis of the 27 shallow fields that cover an overall area of $\sim 26$ square degree with increasing sensitivity towards the phase centers of each pointing according to the ATCA 22 m dish response function, down to $\sigma\lesssim 33\, \mu$Jy. Two additional (deep) fields have been observed to even higher sensitivity. All the SHORES observations have been calibrated to account also for linear polarization. Polarization and deeper field analysis will be presented in future papers. The SHORES shallow-field sample considered in the present paper counts $2294$ sources detected with BLOBCAT to signal-to-noise ratio $SNR\gtrsim 4.5$. Simulations determined that our procedure and final catalog is 95% reliable above $497.5\, \mu$Jy and $95\%$ complete to the $SNR\gtrsim 4.5$ significance level. By exploiting ATCA E-W 6 km configuration we reached resolutions of $3.2\times 7.2$ arcsec, to which level $81\%$ of our sources are unresolved. We determined source counts down to the $150\, \mu$Jy level. For the sources with a counterpart in H-ATLAS, the FIR-radio correlation is calculated and discussed.

We conduct three-dimensional (3D) hydrodynamical simulations of common envelope evolution (CEE) of a neutron star (NS) that launches jets as it spirals in inside the envelope of a rotating red supergiant (RSG) stellar envelope and find that Rayleigh-Taylor instabilities form filamentary ejecta. We first study the 3D RSG envelope properties before we launch the jets. Adding envelope rotation causes the RSG envelope to expand in the equatorial plane and contract along the poles, leading to non-radial oscillations that decay after two oscillation periods, like the radial oscillation of the non-rotating model. In addition, the envelope becomes convective with large vortices, as in the non-rotating case. Since RSG stars oscillate and have envelope convection, we strengthen the claim that there is no need to relax one-dimensional stellar models of cool giant stars when transporting them to 3D grids. When adding jets, the 3D simulations that include pre-set envelope rotation show that envelope rotation leads to more prominent spiral structures of the ejecta than in the non-rotating case. We map the envelope zones that are Rayleigh-Taylor unstable and conclude that this instability forms the filamentary ejecta, with and without envelope rotation. The jet-inflated high-pressure volumes around the NS accelerate the envelope, a process prone to Rayleigh-Taylor instability.

Low-energy (<300 keV) protons entering the field of view of XMM-Newton are observed in the form of a sudden increase in the background level, the so-called soft proton flares, affecting up to 40% of the mission observing time. In-flight XMM-Newton's observations of soft protons represent a unique laboratory to validate and improve our understanding of their interaction with the mirror, optical filters, and X-ray instruments. At the same time, such models would link the observed background flares to the primary proton population encountered by the telescope, converting XMM-Newton into a monitor for soft protons. We built a Geant4 simulation of XMM-Newton, including a verified mass model of the X-ray mirror, the focal plane assembly, and the EPIC MOS and pn-CCDs. We encoded the energy redistribution and proton transmission efficiency into a redistribution matrix file (RMF) and an auxiliary response file (ARF). For the validation, three averaged soft proton spectra, one for each filter configuration, were extracted from a collection of 13 years of MOS observations of the focused non X-ray background and analysed with Xspec. The best-fit model is in agreement with the power-law distribution predicted from independent measurements for the XMM-Newton orbit, spent mostly in the magnetosheath and nearby regions. For the first time we are able to link detected soft proton flares with the proton radiation environment in the Earth's magnetosphere, while proving the validity of the simulation chain in predicting the background of future missions. Benefiting from this work and contributions from the Athena instrument consortia, we also present the response files for the Athena mission and updated estimates for its focused charged background.

Sgr A* is currently very faint. However, X-ray radiation reflected by the Sgr A complex, a group of nearby molecular clouds, suggests that it went through one or more periods of high activity some hundreds of years ago. We aim to determine whether previously proposed physical scenarios are consistent with the observed X-ray variability over the past 25 years, and to characterize the spatial distribution, shape, and internal structure of the clouds. We exploit the full set of XMM-Newton observations, extending the previously studied dataset on variability by at least 12 years. Starting from the recent IXPE result that places the so-called Bridge cloud 26 pc behind Sgr A*, we reconstruct the LOS position of the other clouds, assuming that they were illuminated by a single flare. Additionally, we derive the probability density function (PDF) of the molecular density. We also study the 3D geometry of the complex in case two flares illuminate the clouds. As of spring 2024, the lightfront is still illuminating the Sgr A complex, with the Bridge currently being the brightest cloud. The other clouds in the complex have faded significantly. In the single flare scenario, the Sgr A complex is located $\simeq$ 25 pc behind Sgr A*. In the past 25 years, the illuminated region spans 10-15 pc along the LOS. The derived PDF is roughly log-normal, consistent with previous Chandra results, with a possible high-density excess. Both a single and a multiple flares scenario can explain the observed X-ray variability. Previous concerns about the single flare scenario, raised by shorter monitoring, are now overcome in the 25 years of monitoring. If two flares illuminate the clouds, they must be separated by at least $\sim$ 30 years. We speculate that these clouds are closer to Sgr A* than the nuclear molecular ring at $\simeq$ 100-200 pc and possibly drifting from the ring to the inner region of the Galaxy.

The physical and orbital parameters of Trans-Neptunian Objects (TNOs) provide valuable information about the Solar System's formation and evolution. In particular, the characterization of binaries provides insights into the formation mechanisms that may be playing a role at such large distances from the Sun. Studies show two distinct populations, and (38628) Huya occupies an intermediate position between the unequal-size binaries and those with components of roughly equal sizes. In this work, we predicted and observed three stellar occultation events by Huya. Huya and its satellite - S/2012 (38628) 1 - were detected during occultations in March 2021 and again in June 2023. Additionally, an attempt to detect Huya in February 2023 resulted in an additional single-chord detection of the secondary. A spherical body with a minimum diameter of D = 165 km can explain the three single-chord observations and provide a lower limit for the satellite size. The astrometry of Huya's system, as derived from the occultations and supplemented by observations from the Hubble Space Telescope and Keck Observatory, provided constraints on the satellite orbit and the mass of the system. Therefore, assuming the secondary is in an equatorial orbit around the primary, the limb fitting was constrained by the satellite orbit position angle. The system density, calculated by summing the most precise measurement of Huya's volume to the spherical satellite average volume, is $\rho_{1}$ = 1073 $\pm$ 66 kg m$^{-3}$. The density that the object would have assuming a Maclaurin equilibrium shape with a rotational period of 6.725 $\pm$ 0.01 hours is $\rho_{2}$ = 768 $\pm$ 42 kg m$^{-3}$. This difference rules out the Maclaurin equilibrium assumption for the main body shape.

In recent years, numerical simulations have become indispensable for addressing complex astrophysical problems. The MagnetoHydroDynamics (MHD) framework represents a key tool for investigating the dynamical evolution of astrophysical plasmas, which are described as a set of partial differential equations that enforce the conservation of mass, momentum, and energy, along with Maxwell's equation for the evolution of the electromagnetic fields. Due to the high nonlinearity of the MHD equations (regardless of their specifications, e.g., classical/relativistic or ideal/resistive), a general analytical solution is precluded, making the numerical approach crucial. Numerical simulations usually end up producing large sets of data files and their scientific analysis leans on dedicated software designed for data visualization. However, in order to encompass all of the code output features, specialized tools focusing on the numerical code may represent a more versatile and built-in tool. Here, we present PyPLUTO, a Python package tailored for efficient loading, manipulation, and visualization of outputs produced with the PLUTO code (Mignone et al., 2007; Mignone et al., 2012). PyPLUTO uses memory mapping to optimize data loading and provides general routines for data manipulation and visualization. PyPLUTO also supports the particle modules of the PLUTO code, enabling users to load and visualize particles, such as cosmic rays (Mignone et al., 2018), Lagrangian (Vaidya et al., 2018), or dust (Mignone et al., 2019) particles, from hybrid simulations. A dedicated Graphical User Interface (GUI) simplifies the generation of single-subplot figures, making PyPLUTO a powerful yet user-friendly toolkit for astrophysical data analysis.

We use the Tolman metric to describe gravitational collapse of a sphere of a fluid without pressure in spacetime with the Hubble parameter $H$ related to the cosmological constant. We show that the largest radius of a galaxy formed from such a fluid with mass $M$ is given by $(GM/H^2)^{1/3}$.

This paper explores the idea that information is an essential and distinctive feature of living systems. Unlike non-living systems, living systems actively acquire, process, and use information about their environments to respond to changing conditions, sustain themselves, and achieve other intrinsic goals. We discuss relevant theoretical frameworks such as ``semantic information'' and ``fitness value of information''. We also highlight the broader implications of our perspective for fields such as origins-of-life research and astrobiology. In particular, we touch on the transition to information-driven systems as a key step in abiogenesis, informational constraints as determinants of planetary habitability, and informational biosignatures for detecting life beyond Earth. We briefly discuss experimental platforms which offer opportunities to investigate these theoretical concepts in controlled environments. By integrating theoretical and experimental approaches, this perspective advances our understanding of life's informational dynamics and its universal principles across diverse scientific domains.

Einstein Telescope (ET) is a proposed next-generation Gravitational Wave (GW) interferometer designed to detect a large number of astrophysical and cosmological sources with unprecedented sensitivity. A key target for ET is the detection of a stochastic gravitational-wave background (SGWB), a faint signal from unresolved GW sources. In its proposed triangular configuration, correlated Newtonian noise of seismic origin poses some challenges for the SGWB detection. We study the impact of correlated noise on the SGWB detection and relative parameter estimation for ET in the triangular configuration, comparing it to a 2L configuration with two separated L-shaped detectors. We perform a Bayesian analysis on simulated data, which shows that accurate reconstruction of the SGWB parameters and instrumental noise is achievable if the noise is properly modeled. We illustrate that neglecting correlated noise leads to significant biases in the parameter reconstruction. Our results show that while the 2L configuration provides slightly better parameter estimation precision, mainly due to its longer arm length, the triangular configuration remains competitive when accurate noise modeling is provided.

Dark Matter can interact with electroweak gauge bosons via higher-dimensional operators, in spite of being neutral under gauge interactions, much like neutral atoms interact with photons through Rayleigh scattering. This study explores effective interactions between a real scalar Dark Matter particle, singlet under the SM gauge group, and electroweak gauge bosons. We present a comprehensive analysis of current constraints and projected sensitivities from both lepton and hadron colliders as well as direct and indirect detection experiments in testing Rayleigh Dark Matter interactions. We find that, thanks to the complementarity between collider experiments and cosmological probes, thermally produced Rayleigh Dark Matter at the hundreds of GeV scale can be thoroughly tested with the next generation of experiments. For lighter candidates, upcoming forecasts will explore uncharted parameter space, significantly surpassing the thermal Dark Matter benchmark.

We provide an accessible review to the eta problem in brane-antibrane inflation and a recently proposed solution. The description of the antibrane by means of nonlinearly realized supersymmetry suggests a simple stabilization mechanism for the volume mode, based on the most general Kahler and superpotentials compatible with the symmetries of the problem. This work is a contribution to the proceedings of the second general meeting of the COST Action CA21106 (Cosmic WISPers) and is based on arXiv:2410.00097 .

Late-decaying particles naturally arise in many extensions of the Standard Model, directly impacting key cosmological processes in the early universe, such as Big Bang Nucleosynthesis (BBN). BBN studies often consider electromagnetic energy injection episodes only, but in practice long-lived particles are also amenable to hadronic decays. The latter can greatly alter the predicted abundances of light elements such as $\mathrm{D}/\mathrm{H}$, $Y_p$, ${}^3\mathrm{He}/\mathrm{D}$, and ${}^7\mathrm{Li}/\mathrm{H}$. Incorporating up-to-date measurements, we place constraints on the primordial abundance of long-lived particles as a function of their lifetime. Lastly, we apply our results to the gravitino problem and set bounds on the reheating temperature, which controls the gravitino primordial abundance.

Energy transfer in the dark sector of the universe gives rise to new phenomena of special interest in modern cosmology. When dark energy is modeled as a phantom scalar field, interactions become crucial to avoid Big Rip singularities. In this work, we revisit the phase-space analysis of the field equations by introducing a new set of dimensionless variables distinct from the traditional Hubble normalization approach. These new variables define a compactified phase space for the evolution of physical parameters. We demonstrate that these compactified variables offer fresh insights into the phase-space analysis in gravitational theories, particularly when the dark energy fluid is allowed to possess a negative kinetic energy density.

We investigate the impact of dark matter halos on the gravitational lensing produced by electrically charged, spherically symmetric black holes in the strong-field regime. The study focuses on two dark matter models: the Universal Rotation Curve Model and the cold dark matter model. We derive the coefficients for the strong deflection limit and numerically analyze the deflection angle variations. Graphical representations of the results show that the strong deflection angle, $\alpha_D$ , increases with the charge parameter $Q$ in the presence of a dark matter halo. We explore the astrophysical consequences for the supermassive black holes $M87^*$ and $SgrA^*$ , comparing the results with standard Reissner-Nordström and Schwarzschild black holes via strong gravitational lensing observations. Our findings suggest that charged black holes with dark matter halos can be differentiated from standard black holes. We constrain the charge parameter $Q$ using observational data from the Event Horizon Telescope Collaboration. For $M87^*$ , we find $0 \leq |Q| \leq 0.366M$ with the Universal Rotation Curve model and $0 \leq |Q| \leq 0.364M$ with the cold dark matter model. For $SgrA^*$ , the constraints are $0 \leq |Q| \leq 0.586M$ and $0 \leq |Q| \leq 0.584M$, respectively. These results suggest that charged black holes with dark matter halos satisfy the Event Horizon Telescope constraints, offering potential for future identification in observational campaigns.

We investigate the breaking of dark $SU(2)_d$ symmetry at different temperature scales, occurring after Peccei-Quinn symmetry breaking or following QCD symmetry breaking. We focus on assessing the potential of the hidden monopoles generated during this process to serve as dark matter candidate. Additionally, we examine the impact of axion-monopole interactions on the axion mass. When the phase transition occurs at extremely high temperature ($\sim 10^8 \mathrm{GeV}$), the contribution of monopoles to the axion mass through witten effect becomes non-negligible, playing a crucial role in accurately determining the axion relic density. Moreover, the stochastic gravitational wave background generated by dark phase transition and axionic domain wall collapse may offer a potential explanation for the low-frequency gravitational wave signals observed in PTA experiments.

The formation of the cosmic structures in the late Universe is considered using Vlasov kinetic approach. The crucial point is the use of the gravitational potential with repulsive term of the cosmological constant which provides a solution to the Hubble tension, that is the Hubble parameter for the late Universe has to differ from its global cosmological value. This also provides a mechanism of formation of stationary semi-periodic gravitating structures of voids and walls, so that the cosmological constant has a role of the scaling and hence can be compared with the observational data for given regions. The considered mechanism of the structure formation in late cosmological epoch then is succeeding the epoch described by the evolution of primordial density fluctuations.

Direct detection experiments have started to explore dark matter scattering off electrons and nucleons through light mediators. Mediators with sub-keV masses are efficiently produced in the Sun and can be absorbed in the same detectors that probe dark matter scattering. We investigate the interplay of dark matter scattering and mediator absorption for two models with a dark photon as mediator. For Dirac dark matter, we find that scattering and absorption can be simultaneously observed at direct detection experiments in the near future. For atomic dark matter, we predict additional signals due to scattering of both dark atoms and constituents from ionized dark atoms. In both models, we determine the parameter space that respects bounds from cosmology and astrophysics, where the strongest constraints come from dark matter self-interactions. In this way, we identify viable targets for dark matter with light mediators at upcoming direct detection experiments. Distinguishing between the various signals, for instance by measuring energy distributions, will be crucial to reveal the underlying model in case of a discovery.

We propose a dark matter (DM) model with a complex scalar charged under a hidden gauge symmetry, denoted as $U(1)_D$. The scalar field is the DM candidate while the $U(1)_D$ gauge field $A'$ plays the role of a mediator, which connects the dark sector to the standard model (SM) sector via a tiny kinetic mixing. We find that both the secluded and catalyzed annihilation scenarios can be realized in this model. The phenomenology of DM, including relic density, indirect detection (Fermi-LAT), and CMB (Planck) constraints, is discussed. We also extend our discussion to DM with other spins, including Dirac fermion and vector boson. Our analysis is carried out in two models, denoted as $U(1)_D \times U(1)_Y$ and $U(1)_D \times U(1)_{L_\mu-L_\tau}$, with the former corresponding to $A'$ kinetically mixing with the $U(1)_Y$ gauge field $B$ and the latter corresponding to $A'$ mixing with the $U(1)_{L_\mu-L_\tau}$ gauge field $Z'$. We find that, in previous studies, the indirect detection limits were overly restrictive because they only considered the simplified $2\mathrm{DM} \to 2\mathrm{SM}$ annihilation channel. In contrast, by performing a complete calculation of the gamma-ray and CMB constraints from the process $2\mathrm{DM} \to 2A' \to 4\mathrm{SM}$ in the models we consider, we observe weaker constraints in both the $U(1)_D \times U(1)_Y$ and $U(1)_D \times U(1)_{L_\mu-L_\tau}$ models, with the $U(1)_D \times U(1)_{L_\mu-L_\tau}$ model being subject to the weakest constraints overall since it involves less hadronic decay processes.