Papers for Thursday, Mar 05 2026

Stellar masses are a fundamental property to understand models of pre-main sequence evolution, but their values derived from Hertzsprung-Russell (HR) diagrams are strongly model dependent. We benchmark pre-main sequence stellar evolutionary tracks using stellar masses dynamically estimated by fitting a parametric model to ALMA observations of the $^{12}$CO $(J=3-2)$ line transition emitted by the disks orbiting 20 sources in the old ($4-14$ Myr) Upper Scorpius star forming region. We derive stellar masses from HR diagram fitting for ten different stellar evolutionary models, which we then compare with their stellar dynamical masses for comparison in the stellar mass range $0.1-1.3 \> M_\odot$. Models with a moderate-to-low fraction of cold stellar spots ($f=17\%$) most accurately reproduce the dynamical stellar masses ($100\%$ of the targets agree within $\pm1\sigma$). While a higher spot coverage ($f=34\%$) provides similar stellar mass predictions similar to magnetic equipartition models, larger fractions ($f\geq51\%$) significantly disagree with dynamical masses. Magnetic equipartition models overestimate stellar masses up to a factor $\sim20\%$, whereas non-magnetic models underestimate them up to $\sim12\%$. For some models, there is evidence that the stellar mass discrepancies are anticorrelated with dynamical stellar masses. When stellar dynamical mass priors are considered in HR diagram fitting, the median age of a single source can change up to $\sim25\%$, while the median ages inferred across different tracks become consistent, with the age scatter decreasing by $\gtrsim77\%$. These results provide strong empirical constraints for testing and developing evolutionary models of pre-main sequence stars.

We present a JWST analysis of intracluster light (ICL) in XLSSC 122 at z = 1.98, currently the most distant known strong lensing galaxy cluster with an evolved member population. Using deep JWST imaging complemented by HST data and careful control of systematics, we robustly detect diffuse emission extending to several hundred kpc from the brightest cluster galaxy (BCG) down to about 29 mag arcsec^-2. Multi component PSF convolved Sersic modeling separates the surface brightness profiles into three components: a BCG core, a BCG envelope, and an ICL component, with stable Sersic indices across wavelengths. Nearly flat color profiles indicate minimal radial variation in the stellar populations of the BCG envelope and the ICL. The median ICL fraction measured across seven bands is about 17 percent, demonstrating that the buildup of intracluster stars in massive halos was already well underway by z about 2. The ICL fraction peaks near 5000 Angstrom in the rest frame, resembling the behavior observed in dynamically active clusters. We also detect a southern excess of ICL relative to the best fit Sersic model and quantify it using wavelet based modeling, providing additional support that this system is dynamically active. The BCG + ICL light distribution and strong lensing mass map show strong morphological agreement within about 100 kpc. These findings establish the ICL as an early forming and dynamically informative component of massive halos.

We present the first systematic census of quenching mechanisms using kinematic asymmetries in a large sample of $\sim$6,700 galaxies from the MaNGA survey, providing a unified view of what halts star formation in the local Universe ($z<0.2$). We quantify stellar and nebular gas disturbances through the higher-order terms of a Fourier series expansion. These asymmetries serve as powerful diagnostics, as different quenching mechanisms leave distinct kinematic signatures on gas and stars. Our analysis reveals that the most effective quenching pathways leave minimal kinematic imprints by the time galaxies are fully quenched. This "kinematic regularity" points toward slow-acting processes (>3 Gyr) such as starvation and maintenance feedback. A striking finding emerges from our mass-matched analysis: quenched symmetric satellites are significantly more compact than their asymmetric counterparts ($3.4\sigma$), a trend that is even more pronounced for symmetric centrals ($12.3\sigma$). Our results suggest that environment drives the dominant satellite quenching pathway through rapid gas stripping followed by long-term starvation. These compact, kinematically undisturbed satellites (the most representative case within our sample) have undergone intense gas stripping and central compaction, creating bulge-like structures with old, metal-rich stellar populations. Combined with halo gas cut-off and the prevention of cosmological accretion due to starvation, this creates an irreversible quenching path. Conversely, the larger sizes of disturbed, quenched centrals are consistent with merger-driven growth. Internal processes, likely driven by the AGN cycle over 1-3 Gyr that prevents hot halo gas cooling, sustain quenching maintenance in this population. The absence of asymmetric satellites in the star-forming regime suggests environmental quenching operates without significant kinematic perturbation.

Galaxy mergers have long been proposed as a mechanism for funneling gas toward galactic centres, potentially triggering accretion onto supermassive black holes (SMBHs) and igniting active galactic nuclei (AGN). While simulations often support this scenario, observational studies have yielded conflicting results regarding the AGN-merger connection. In this study, we analyze 31 galaxies from cosmological zoom-in simulations spanning redshifts $0.5 < z < 3$. We identify mergers using detailed merger trees based on six-dimensional dark matter particle information and identify AGN activity through SMBH accretion histories. To bridge the gap between simulations and observations, we generate mock JWST-like images and extract non-parametric morphological parameters. Employing a $k$-nearest neighbours (KNN) classifier in a five-dimensional space (four morphological parameters and redshift), we identify mergers in the mock-observed dataset. Our analysis reveals a statistically significant enhancement of AGN activity in merging systems, particularly at lower redshifts ($0.5 < z < 0.9$), where central gas reservoirs are more depleted. This supports the view that mergers contribute more significantly to AGN triggering in environments with low internal gas reservoirs, while their impact may be less pronounced in gas-rich systems. However, when relying solely on morphological classifications from mock observations, the observed AGN-merger connection weakens, especially at higher redshifts. This underscores the challenges in detecting merger-induced AGN activity observationally and highlights the importance of combining simulations with realistic mock observations to fully understand the AGN-merger relationship.

Galaxy clusters offer powerful insights into the large-scale structure of the Universe and the physics of baryons in hot state. Their scientific exploitation, however, hinges on our ability to accurately measure key thermodynamic properties. In this work, we aim to assess the reliability of current analysis techniques in reconstructing these properties, with particular focus on samples similar to those observed in the Cluster HEritage project with XMM-Newton (CHEX-MATE). We develop a suite of dedicated end-to-end simulations of CHEX-MATE-like clusters selected from large scale hydrodynamical simulations, and processed through a newly developed realistic XMM-Newton simulator. We apply a full X-ray data analysis pipeline to the mock datasets, including imaging, spectral fitting, and profile reconstruction. The gas density profiles can be robustly recovered across a wide radial range, when using azimuthal mean surface brightness profiles. Our reconstruction techniques are able to reproduce the intrinsic density profile with the correct scatter, with deviations of at most 10% between 0.1 and 1xR500c. The gas mass is reconstructed with better than 1% accuracy. Accurate measurement of temperature profiles is more challenging and possibly subject to biases, particularly in the presence of azimuthal variations and multi-temperature gas along the line of sight, which dominate over projection effects. Our results highlight the need for caution in interpreting cluster temperature measurements and underscore the value of tailored mock observations for understanding observational systematics. These findings also suggest that biases in X-ray temperature measurements may alter the interpretation of the thermodynamical state of the intra-cluster medium, an outlook particularly relevant in light of recent low velocity measurements from the XRISM mission.

Stage-IV photometric galaxy surveys are designed to measure the position and shapes of billions of galaxies. Their aim is to characterise the large-scale distribution of matter in the Universe using galaxy clustering and weak gravitational lensing. As a byproduct, stage-IV surveys are expected to detect more than a hundred thousand strong gravitational lenses. In this article, we propose the use of weak-lensing perturbations to strong lenses, specifically their line-of-sight (LOS) shear, as a cosmological probe. This new observable allows us to define three new correlation functions: the LOS shear with itself, with galaxy positions, and with galaxy shapes, thereby promoting the standard $3\times 2$pt correlation method to a $6\times 2$pt scheme. We design estimators for these new correlation functions and determine their expectation values as a function of the matter power spectrum. We then derive the analytical expression for the full covariance matrix of the $6\times 2$pt correlation scheme. Considering various scenarios for the stage-IV strong-lensing samples, we demonstrate that the cosmological information carried by the LOS shear of strong lenses will be detectable with a very high signal-to-noise ratio, even in the most pessimistic of cases. Strong lenses are thus extremely promising cosmological probes, whose synergy with galaxy positions and shapes should also contribute to mitigating systematics in stage-IV surveys.

The classical picture that planet formation occurs in protoplanetary disks that are isolated from their environment is undergoing a major shift toward a more connected picture. An increasing amount of evolved disks are found to be actively interacting with their environment, often showing various types of spiral structures. In this work, we aim to investigate if these spirals can be a direct result of ongoing late infall using the grid-based 3D hydrodynamics code FARGO3D. We perform a detailed analysis of the spiral properties and appearance in scattered light and CO line emission using the radiative transfer code RADMC3D. In scattered light, we find both well-defined spirals with few arms (m=2) and more flocculent structures: The gradual accretion of gas remnants after a major accretion event has the most success in the former, whereas active accretion via streamers favors the latter. The m=2 spirals we find have a very low pattern speed, making them easily discernible from spirals caused by a perturber. We also find spiral patterns in the $^{12}$CO residual motions, but their morphology does not match the one found in scattered light. The disk perturbations are strongest in the upper layers (z>4H), which is reflected by the reduced amplitude of the residual motions in the more optically thin $^{13}$CO emission. Moreover, we find that the formation of m=2 spirals is not promoted in disks with lower mass, despite being more susceptible to deeper kinematic perturbations. While the late-infall streamers impact planet formation directly through the delivery of fresh material, we show that the midplane remains unperturbed unless the infalling mass is of the same order of magnitude as the disk mass. Planet formation can therefore only be impacted by late infall through secondary mechanisms that lead to dust trapping or the generation of turbulence starting from surface-level perturbations.

Analyses of baryon acoustic oscillations (BAO) commonly employ template-based methods to extract compressed parameters from the clustering of dark-matter tracers, which are then interpreted in terms of ratios of the sound-horizon scale and cosmological distances relative to a fiducial cosmology. A small mismatch between the sound-horizon scale derived from the standard analytic formulation (integral over the sound speed) and the effective scale imprinted in clustered matter can, however, introduce a systematic bias in cosmological inference. We extend previous work to a broader class of cosmological models, quantify this bias for surveys with DESI-like precision, and propose strategies to correct for the effect. We find that the induced bias becomes a significant fraction of the statistical uncertainty for deviations from the fiducial cosmology, at the level of $|\Delta \Omega_m| = 0.03$ and $|\Delta N_\mathrm{eff}| = 0.3$, and for very precise data corresponding to a forecasted Year-5 DESI survey (or other stage IV dark energy galaxy surveys). We present several ways to correct for this effect, suitable for a variety of applications. We therefore recommend that analyses exploring such parameter regimes either apply the proposed corrections or include an appropriate systematic error budget.

The precision era of multi-messenger astronomy, together with modern astrophysical, cosmological, and gravitational wave observations, increasingly points toward the existence of a ``dark" sector that cannot be explained within the framework of the Standard Model of particle physics and General Relativity. In this review, we explore extensions of standard physics and examine how observational data can be used to probe new particles and interactions. We consider a wide range of scales, from Solar System tests to galactic and cosmological observations, and investigate both conventional dark matter candidates, such as weakly interacting massive particles, and alternative scenarios including ultralight fields and primordial black holes. We discuss constraints derived from compact objects such as neutron stars, black holes, pulsars, and magnetars observations as well as from high-energy astrophysical phenomena. In addition, we analyze extensions of General Relativity involving additional scalar fields and their impact on gravitational wave signals and stochastic backgrounds from primordial black holes. We also study the capture and accumulation of dark matter in compact objects, which can alter properties such as mass, radius, and tidal deformability, and consider scenarios in which dark matter decays into Standard Model particles. While current observations already place significant limits on dark matter and modified-gravity models, upcoming experiments and observatories are expected to further probe or discover such new physics by improving constraints on particle masses and interaction strengths.

The STIS team re-derived on-orbit sensitivity curves for the echelle modes for post-servicing mission 4 observations using the standard DA white dwarf G 191-B2B. These new updates relied on the recent CALSPECv11 models, which introduced improvements in the fluxes of the primary standard stars of the order of~1-3% depending on the wavelength of interest. As part of this effort, the team also released new blaze shift coefficients and echelle ripple tables. We present a detailed description of the procedures followed in the derivation of these new throughputs and the accompanying updates.

Constraining stellar models using asteroseismic and spectroscopic observations is a powerful method for precisely determining the fundamental properties of stars in different kinematic components of our galaxy. We use spectroscopy and individual oscillation mode frequencies to perform a homogeneous modeling study of eight evolved metal-poor stars enhanced in $\alpha$-elements. We compare a full treatment of $\alpha$-enhancement against an ad hoc correction to the total metallicity and show that the stellar properties inferred from asteroseismic modeling using both sets of models agree with each other. Additionally, we find that the uncertainties on stellar parameters derived from the both $\alpha$-enhanced modeling methods are comparable. This is in qualitative disagreement with existing works showing red-giant ages constrained by only the global asteroseismic parameters to depend strongly on the opacities and abundances assumed in 1D modeling. We also show that the observed frequency of maximum oscillation power ($\nu_{\text{max}}$) is larger than the value predicted from applying the $\nu_{\text{max}}$ scaling relation to the masses, radii, and temperatures inferred from the detailed modeling. This discrepancy is pronounced at low metallicities, consistent with recent findings indicating a breakdown of the $\nu_{\text{max}}$ scaling relation for metal-poor stars. Understanding the extent to which the $\nu_{\text{max}}$ scaling relation fails for low-metallicity solar-like oscillators through detailed modeling will enable more accurate mass and age determinations for hundreds of giant stars in the Galactic Halo for which only global asteroseismic parameters are available.

We report the discovery of significant X-ray polarization from the dynamically confirmed black hole X-ray binary (BHXB) GS 1354-64 during its 2025-2026 outburst, obtained with the Imaging X-ray Polarimetry Explorer (IXPE). The observation, obtained shortly after a bright X-ray flare, captures the source in an intermediate state following a stalled (failed) state transition. We discover significant 2-8 keV polarization at the ~4% level with high statistical support--14-sigma significance from frequentist analysis and log Bayes Factor 283+/-1 in a Bayesian framework--measuring PD 4.0+/-0.2% and PA-1+/-2 degrees (90% credible interval). The PD exhibits a statistically significant increasing trend with energy--the strongest such increase yet observed by IXPE in a BHXB--going from 2.1+/-0.3% in the 2-3 keV band to 11+/-3% in the 6.5-7 keV band, while the PA appears stable across both energy and time to within statistical uncertainties. Timing analysis of the IXPE data reveals a ~5 Hz Type-C quasi-periodic oscillation. IXPE + NuSTAR spectropolarimetric modeling suggests that the data can be described by polarized thermal disk and Comptonized components with PAs differing by ~90 degrees, or by a dominant Comptonized polarized component whose effective PD increases across the IXPE bandpass--the inferred component-level polarization levels are therefore model-dependent. In either picture, GS 1354-64 retains a strong coronal component during the transitional period observed by IXPE. These results illustrate how X-ray polarimetry can provide a sensitive diagnostic of the accretion state and geometry in black hole X-ray binary accretion flows, exploring a liminal phase at the cusp of state transition.

In this work, we present comprehensive observational constraints on the time-varying vacuum $\Lambda(t)$CDM cosmology using the latest baryon acoustic oscillation (BAO) data from Dark Energy Spectroscopic Instrument (DESI) Data Release 2 measurements in combination with cosmic chronometer $H(z)$ data. We explicitly quantify the impact of vacuum dynamics on the expansion history, the total effective equation of state parameter of the unified cosmic fluid, the effective dark energy equation of state parameter, and the deceleration parameter in the spatially flat $\Lambda(t)$CDM model. We perform a full Markov Chain Monte Carlo (MCMC) analysis and statistical model comparison, providing a consistent assessment of the $\Lambda(t)$CDM model relative to the standard $\Lambda$CDM framework. Our results demonstrate that $H(z)$ and BAO observations strongly constrain deviations from the $\Lambda$CDM model, driving the vacuum dynamics parameter toward $\alpha \simeq 0$, while significantly reducing parameter degeneracies and alleviating the Hubble tension.

The noise in bias frames for all four readout amplifiers in the Advanced Camera for Surveys (ACS) Wide Field Channel (WFC) is dependent on row number. This is because dark current accumulated during readout increases across the detector, influencing and increasing the read noise as a function of row number. In this report, we investigate bias frames taken with the ACS/WFC to explore the column dependence of read noise for each of the amplifiers for different anneal periods. Analyzing the data, we find that there is no column dependence of read noise and that the read noise values for the physical pre-scans are approximately 0.5 e$^-$ lower than in the science arrays because there is no readout dark accumulated in this area. We further investigate 1) the evolution of read noise over an anneal period, 2) a linear decrease in read noise within the initial columns per amplifier, and 3) pixels in elevated read noise columns. We conclude that 1) there is no visual trend of read noise over an anneal period, 2) amplifiers A and C have an initial linear decrease of read noise in the science arrays, and 3) masking unstable hot pixels in a column will decrease its read noise values.

We present the spatial part of the point source signal extraction strategy for the upcoming CHORD galaxy survey. CHORD, the Canadian Hydrogen Observatory and Radio-transient Detector, is an under-construction drift-scanning compact interferometric radio telescope. CHORD comprises 512 six meter dishes and observes in the 300 to 1500 MHz frequency range. One of its science goals is producing a catalogue of galaxies detected by the neutral hydrogen (HI) 21 cm emission line. CHORD's highly redundant dish layout creates the problem of spatial aliasing, the effect where the same signal could be feasibly produced from sources at multiple locations on the sky. The search will be done with a matched filter in the visibility plane. This paper presents the search strategy and a prediction tool that can quickly estimate the matched filter response at a given sky position, allowing a prediction of alias locations and severity. This tool confirms that although aliases are impossible to distinguish in a single snapshot, they become possible to distinguish when combining data over a period of time. It predicts that aliases will be harder to distinguish for observations closer to the celestial equator, but that scanning with offset adjacent strips can remove this degeneracy. It predicts that the optimal strategy for a single offset to disambiguate aliases is to re-point the array in declination by about two degrees. A future paper will combine these findings with realistic noise estimates and galaxy population statistics to make forecasts of the population of galaxies that CHORD will detect.

We test whether parsec-scale stellar systems in the Milky Way follow the galactic radial acceleration relation (RAR) or the baryonic TullyFisher relation (BTFR). We analyse 5646 Gaia DR3 open clusters from the Hunt \& Reffert catalogue. Observed accelerations are derived from velocity dispersions and characteristic radii, and baryonic accelerations from stellar masses and characterisitc radii. The clusters are placed on the RAR and BTFR planes and compared with Newtonian and MOND expectations. Approximately 90 per cent of open clusters (those with $N_\star \leq 250$) lie close to the RAR, albeit with significant scatter. In a first-of-its-kind test, a smaller fiducial sample is consistent with a best-fitting acceleration scale $g_\dagger \approx 1.2 \times 10^{-10}\ \mathrm{m\,s^{-2}} \pm 0.5$ dex, compatible with canonical MOND values. More massive clusters approach the Newtonian virial expectation. No correlations are found between RAR residuals and galactocentric radii, distance to the Galactic disk midplane, age, or morphology. Tidal effects and unresolved binaries are insufficient to reproduce the observations without fine-tuning. Interpreted within a MOND framework, the alignment of most open clusters with the RAR and BTFR suggests that low-acceleration dynamics operate on parsec scales within the Milky Way. This implies that the Galactic gravitational field is not smooth on these scales and may include regions where the total gravitational acceleration falls below $a_0$, partially mitigating the external field effect, thereby motivating higher-resolution modelling of the Galactic potential and informing other small-scale gravity tests within the Galaxy.

We present a systematic analysis of 2544 cataclysmic variable systems and related candidates observed by the Transiting Exoplanet Survey Satellite (TESS), with the aim of compiling a comprehensive catalogue of orbital periods. Using 2-minute photometric time-series data, we applied an automated algorithm to generate Lomb-Scargle periodograms and identify the most significant coherent periodic signals, which were subsequently verified through visual inspection. This process yielded a confident sample of 910 sources exhibiting at least one periodic signal, hereafter referred to as the Cataclysmic Variable Confident Catalogue (CCC). For each object, we report the most likely orbital period together with additional periodic features such as spin modulations and/or superhump signals when present. To assess consistency with previously published measurements, we cross-matched the CCC with the Ritter \& Kolb catalogue, identifying 300 overlapping systems, of which 215 showed full agreement with the R\&K orbital periods, while 39 displayed discrepancies for which the R\&K values were revised based on our TESS measurements and supporting evidence from the literature. Overall, the CCC provides a means to validate known orbital periods, propose corrections where necessary, and offer new determinations for systems with previously unknown periods, thereby supporting a more complete and reliable characterisation of the cataclysmic variable population.

Thermal imaging of debris disks is resolving the vertical height in an increasing number of systems, enabling the use of the vertical structure to decode the dynamical state of the planetary system. In this study, we examine the multi-wavelength structure of the archetypical edge-on debris disk of $\beta$ Pic, extensive imaging of which across mid-infrared to millimeter wavelengths makes it the prime system to study the vertical height across different grain size populations. We non-parametrically modelled the radial profiles and constrained the vertical height at each wavelength while taking into account the vertical warping, finding the disk to be on average 1.5 times thicker vertically in the mid-infrared compared to the millimeter and the scale height to be relatively constant across radius. The decreasing scale height with wavelength is in contrast to predictions from collisional damping, and could be a result of the combined effect of radiation pressure and random collisions. We also show that the disk is warped at millimeter wavelengths and find tentative evidence for clumps in ALMA images which will require follow-up observations to confirm. The millimeter vertical warping is consistent with findings in scattered light and the secular perturbation interpretation due to the inner giant planets, which could also explain the relatively constant apparent scale height across radius, and potentially earlier findings of a non-Gaussian vertical profile which this study confirms.

A recently proposed generalized entropy by Luciano and Saridakis extends the standard Boltzmann-Gibbs and Bekenstein-Hawking framework through a microscopically motivated construction involving two independent entropic exponents. When applied within the gravity-thermodynamics correspondence, this entropy leads to a modified cosmological dynamics that can be interpreted as an effective dark energy sector of entropic origin, while recovering $\Lambda$CDM in appropriate limits. In this work, we perform the first observational confrontation of the resulting entropic cosmology at the background level. Focusing on the case $\alpha_\delta=0$, we constrain the model using Cosmic Chronometers, Pantheon$^+$ Type Ia supernovae calibrated with SH0ES, BAO measurements from DESI DR2 and compressed Planck 2018 CMB information. We find that the model yields a statistically robust fit to the combined data sets and can simultaneously satisfy Pantheon$^+$, SH0ES and CMB shift-parameter constraints, unlike $\Lambda$CDM. Although the entropic parameters remain close to their standard values, the $\Lambda$CDM limit is excluded at the $2\sigma$ level within the restricted parameter space considered. These results indicate that the Luciano-Saridakis entropic cosmology offers a viable extension of the standard model with the potential to alleviate the Hubble tension at the background level.

We report the discovery of a Lyman-alpha emitter (LAE) protocluster at z = 4.90 in the COSMOS field, comprising four distinct overdensity peaks spanning ~65 x 36 cMpc^2, with the primary concentration exhibiting a 4-fold surface density enhancement relative to the field within a 1.5 proper Mpc (pMpc) radius. Using SILVERRUSH narrowband survey data combined with JWST COSMOS-Web imaging, we perform a first systematic rest-frame optical and UV morphological comparison of protocluster versus field LAEs at this redshift using JWST NIRCam rest-frame UV (F150W, ~2540 Angstrom) and optical (F277W, ~4700 Angstrom) imaging. Sersic profile fitting for 16 protocluster members and 23 field LAEs reveals a significant size difference: protocluster LAEs are ~40% larger in rest-optical (median R_e = 0.81+0.26/-0.04 kpc vs. 0.58+0.11/-0.04 kpc, p = 0.041) with no significant difference in rest-UV (p = 0.51) or Sersic index. At fixed stellar mass, protocluster LAEs are offset by +0.12 dex (~31%) in rest-optical size from the field size-mass relation (68% CI: [+0.08, +0.21]; Mann-Whitney p = 0.033), with 75% exhibiting positive size residuals compared to 44% of field LAEs. This wavelength-dependent environmental signature suggests that protocluster environments at z ~ 5 preferentially affect extended stellar populations, possibly through tidal interactions or an earlier onset of star formation in the dense environment, with no significant environmental difference detected in rest-UV sizes, providing observational evidence for environmental influences on the structure of LAEs during the early build-up phase of cosmic star formation.

Accretion is the primary driver of protostellar evolution, regulating mass assembly and shaping the physical and chemical environments of young stellar objects. Quantifying accretion in the Class 0 protostellar phase is particularly important, yet remains observationally challenging due to high extinction toward the central protostars. In this paper, we present JWST NIRSpec and MIRI/MRS IFU data towards the Class 0 protostar L1527 IRS. We extract one-dimensional spectra and find emission from atomic and molecular hydrogen, water, OH, and several ionic species. The atomic hydrogen lines, Br$\alpha$, Pf$\alpha$, and Pf$\gamma$ are the most critical to this study since they can be used as accretion diagnostics. The existence of these atomic hydrogen lines viewed in scattered light indicates that accretion is likely occurring magnetospherically rather than through a boundary layer. Moment 0 emission maps show that the hydrogen emission is co-spatial with the scattered light continuum with a strong east-west asymmetry which is not due to outflow shocks. We additionally present moment 0 maps of other detected species and discuss their emission morphology. By primarily analyzing the Br$\alpha$ line, the strongest of our detected atomic hydrogen lines, we characterize the accretion onto L1527 IRS by estimating the accretion luminosity to be $0.4~\text{L}_\odot$ and the accretion rate to be around $1\times10^{-7}~ \text{M}_\odot \text{yr}^{-1}$. We lastly discuss the implications of our results with respect to both non-steady and asymmetric accretion possibly occurring in L1527 IRS.

Modern datasets have the precision necessary to uncover new information by including higher-order, non-Gaussian information into cosmological inference. The marked power spectrum offers access to such information while preserving the structure of two-point correlators. This approach to higher-order statistics has the advantage that many modeling questions can directly benefit from progress already made in standard cosmological analyses using the power spectrum and correlation function, while increasing the data vector size negligibly and retaining much of the degeneracy-breaking power of the bispectrum. In this work, we first restructure the marked power spectrum to isolate its higher-order information and demonstrate its ability to break parameter degeneracies. We then investigate the effect of survey geometry on the marked power spectrum and find that a treatment similar to that of the power spectrum is sufficient. Additionally, we investigate the perturbative modeling and covariance structure of the marked power spectrum, shedding light on its degeneracy breaking power and cross-covariance with the power spectrum. Finally, we demonstrate that the cosmology dependence of the marked power spectrum is smooth, indicating that cosmological inference is possible by modeling the cosmology dependence through interpolation rather than analytical modeling.

Cosmological reionization was a highly out-of-equilibrium event that affected every parcel of the intergalactic medium, making it a candidate for astrophysical generation of intergalactic magnetic fields. During reionization, the first stars and galaxies ionized the surrounding, largely neutral, medium in ever-expanding envelopes. Photoionization from sources on one side of the front, combined with the quadrupolar angular dependence of the photoionization cross section, leads to an anisotropic electron velocity distribution. We investigate instabilities in these reionization fronts as a mechanism to generate seed magnetic fields. The Weibel instability has the potential to create a magnetic field from these anisotropies. We calculate the magnitude of the isotropic and anisotropic distribution within a simulated reionization front. We find that the fractional anisotropy can grow to $6\times 10^{-3}$ toward the middle of the ionization front. We show that the linear growth timescale of the Weibel instability is fast compared to the crossing time of the ionization front ($\sim 2\times 10^5$ seconds). We briefly speculate on the possible non-linear evolution of the instability and the implications for cosmological magnetogenesis.

M dwarfs are prime targets for exoplanet searches due to their low masses and radii, which enable the detection of small planets in their habitable zones (HZs). However, the magnetic activity of M dwarfs can introduce signals in radial velocity measure- ments that may be mistaken for planetary signatures, making the understanding of stellar activity cycles crucial for accurate planet detection and characterisation. We aim to identify and characterise long-term magnetic activity cycles in M dwarfs using a homogeneous and extensive spectroscopic dataset in order to better understand their magnetic variability and its implications for exoplanet detection. We analysed 13 years of high-resolution spectra obtained with the SOPHIE spectrograph for two early M dwarfs known to host exoplanets. We simultaneously monitored chromospheric activity using two indicators, the H{\alpha} index and the Mount Wilson S -index. Long-term trends were modelled using both sinusoidal and low-order polynomial fits to robustly identify stellar activity cycles. For GJ 617A, we report a cycle of approximately 4.8 years, while for GJ 411, we find several characteristic timescales of variability of about 4.9 years. In addition, TESS photometric data reveal signs of short-term variability in GJ617A. The periods of the long-term variability detected for GJ 617A and GJ 411 do not coincide with any of the planetary signals previously reported, which reinforces the hypothesis that they are of magnetic origin. If indeed the variability is due to activity, the cycles detected would not be driven by the same mechanism: The cycle in GJ 617A is consistent with a solar-like dynamo, while the rotation seems to play a different role in the long-term cycles detected in GJ 411.

Neutrino astronomy provides another window to exploring the Universe, exemplified by the detection of a megaelectronvolt neutrino burst from the core-collapse supernova (CCSN) SN~1987A (refs.~\citenum{hir+87,bio+87}). Commonly discussed theories suggest that some CCSNe could produce neutrinos with energies a thousand times more than those of SN~1987A \cite{tm18}, which has been probed with new-generation facilities \cite{abb+12,aar+15,abb+23}. The interaction of SN ejecta with a dense circumstellar medium (CSM) or a jet, launched in a CCSN, being choked in the stellar envelope of the progenitor or an outside CSM are both well-accepted scenarios for the high-energy neutrino production. Here we report the detection of a high-energy neutrino flare at a 3.9$\sigma$ significance from SN~2017hcd, made by our analysis of the public track-like neutrino data taken by the IceCube Neutrino Observatory \cite{IceCube17}. A Type IIn SN with optical emissions arising from the ejecta--CSM interaction, SN~2017hcd's neutrino flare lasted $\sim$1--2 month, with its central time $\sim$14-day prior to the SN's optical discovery time. Its estimated isotropic neutrino energy (all flavors) is approximately two orders of magnitude higher than the energy ($\sim 10^{50}$\,erg) carried in the SN's ejecta, too high to be explained with the ejecta--CSM scenario. Thus, a choked jet may be the source of the neutrino flare.

Based on the Ca\,II H and K lines observed by LAMOST, we employ the photospheric ($R'_{\rm HK}$) and basal ($R^+_{\rm HK,L}$) flux calibrated chromospheric activity indices to examine the relationship between chromospheric activity and the stellar rotation rate. We identify the rotation periods of 11,108 stars observed by Kepler and TESS by cross-matching our chromospheric activity catalog with previous studies. Our statistical results show that chromospheric activity increases with the rotation rate until it reaches a saturation level. As the stellar effective temperature increases from 4950 to 5850 K, the saturation values of the rotation period ($P_{\rm rot}$) vary correspondingly from 4.38 to 1.23 days for $R'_{\rm HK}$ and from 9.88 to 1.33 days for $R^+_{\rm HK,L}$. Similarly, the corresponding saturation Rossby number Ro ranges from 0.200 to 0.032 for $R'_{\rm HK}$ and from 0.302 to 0.107 for $R^+_{\rm HK,L}$. The saturation is also found to be significant in stars with thick convective zones, whereas it is less apparent in stars with higher effective temperatures. For solar-like stars in the $T_{\rm eff}$ range of 4800 to 6000 K, The values of chromospheric activity indicators are saturated when $P_{\rm rot}<1.45 $ days (Ro$<$0.100) and $P_{\rm rot}<2.85 $ days (Ro$<$0.097) for $R'_{\rm HK}$ and $R^+_{\rm HK,L}$, respectively.

A machine-learning-based method is developed to identify objects with unusual stellar spectra. The method employs an autoencoder, a neural network trained to compress spectral data into a low-dimensional representation and subsequently reconstruct it. Spectra that deviate significantly from the dominant patterns in the training dataset are identified using the reconstruction error as an anomaly score. The models are applied to selected datasets from the MaNGA Stellar Library, an empirical library of stellar spectra. Several spectra are flagged as anomalous: an object with likely instrumental and/or reduction issues, two carbon stars, and an oxygen-rich thermally pulsating asymptotic giant branch star. The sources of the large reconstruction errors are examined, and the effectiveness and limitations of autoencoder-based approaches for detecting anomalous stellar spectra are discussed.

We investigate the accretion dynamics of the black hole X-ray binary Swift J1727.8-1613 during its $2023-2024$ discovery outburst that lasted for $\sim10$ months. Insight-HXMT monitored the rising phase of the outburst of Swift J1727.8-1613 roughly continuously from 2023 Aug 25 to 2023 Oct 05. Strong signatures of type-C Quasi-Periodic Oscillations (QPOs) are observed during this phase of the outburst. In our recent paper, nature of the QPOs are studied with the propagating oscillatory shock (POS) model. In this paper, we report on the observation of both positive (or hard) and negative (or soft) time-lags in the $4-10$ keV (LE), $10-30$ keV (ME), and $30 -150$ keV (HE) bands, computed with respect to the $2-4$ keV reference band. We detect a clear transition from hard to soft lags as the outburst evolves. We show the evolution of QPOs and associated time-lags between different X-ray energy bands, correlated with changes in the QPO frequency, spectral state, and the size of the Comptonizing region. Our analysis reveals strong anti-correlations between the time-lags and both QPO frequency and photon index, and a strong positive correlation with the shock location. These evolving lag characteristics and their correlations provide crucial insights into the changing accretion geometry and the interplay of radiative processes, further supporting dynamic models like the POS in explaining the coupled spectro-temporal evolution in black hole X-ray binaries.

Clouds are the largest source of uncertainty in climate simulations. For exoplanets, cloud simulation is particularly challenging because of the lack of observational data to tune parameterized cloud models. Here we apply Community Aerosol and Radiation Model for Atmospheres (CARMA), a size-resolved bin cloud microphysics model, to the atmospheric global climate model Community Atmosphere Model (CAM6) and simulate exoplanets with a range of planetary rotation rates. CARMA produces fewer liquid clouds than the native CAM6 parameterized cloud microphysics scheme (Morrison-Gettelman two-moment microphysics, MG), more ice clouds, and a significantly different ice cloud size distribution. Overall, this leads to a decrease in the magnitude of the net CRE by 4-10 $W/m^2$, which is unlikely to change the determination of habitability from a climate perspective in most cases. The difference in ice cloud size distribution is likely to strongly affect transmission spectral retrievals. Our work confirms that the MG parameterized cloud microphysics scheme can produce reasonable climate simulation when extrapolated to some exoplanet contexts and highlights the value of resolved cloud microphysics for evaluating parameterized schemes and for interpreting observations.

We measure galaxy stellar mass functions (SMFs) for nine of the most massive galaxy clusters in the local universe ($0.07 < z < 0.11$) using deep and complete spectroscopy from the MAssive Cluster Survey with Hectospec (MACH). We construct the cluster SMFs down to $\log(M_*/M_\odot) \gtrsim 8.5$. For comparison, we measure the SMF for field galaxies, complete to $\log(M_*/M_\odot) \gtrsim 10.5$, based on Sloan Digital Sky Survey (SDSS) spectroscopy over the same redshift range. The mean MACH SMF shows a shape similar to that of the field SMF but with a significantly higher amplitude at $\log(M_*/M_\odot) < 11.4$. At $\log(M_*/M_\odot) > 11.4$, the MACH SMF shows a clear excess, indicating the contribution of massive galaxies, including Brightest Cluster Galaxies (BCGs). Based on homogeneous MACH spectroscopy, we compare SMF shapes for quiescent and star-forming members as a function of cluster-centric distance. The quiescent SMFs display a curved shape with a peak at $\log(M_*/M_\odot) \approx 10.5$; the star-forming SMFs decline monotonically with increasing stellar mass. We further compare the mean MACH SMF with SMFs derived from similarly massive clusters in the IllustrisTNG-300 simulations. The shape of the observed and simulated SMFs agree well overall. However, the MACH clusters contain roughly a factor of two more galaxies at $9.0 < \log(M_*/M_\odot) < 10.5$. These results demonstrate that constructing cluster SMFs from complete spectroscopic samples can test simulations and provide powerful constraints on galaxy formation and evolution in dense environments.

We study the impact of supermassive black hole (SMBH) growth, $\langle \dot{M}_\mathrm{SMBH}\rangle$, major and minor galaxy mergers, and gas processes, on the average gas metallicity of galaxies, with the aim to uncover which of these processes drive the scatter in the gas metallicity-stellar mass relation (MZR) at different redshifts in nodes, filaments and voids. At $z=5$, minor mergers produce the largest differential in $\log[Z_g/Z_\odot]$ for all environments, where the node population displays a maximum $0.38$ dex increase in the average $\log[Z_g/Z_\odot]$ compared to non-merging galaxies. The node population also displays a consistent $0.1$ dex reduction in $\delta \log[Z_g/Z_{\odot}]$ across all redshifts, whilst filament and void galaxies show a lower magnitude of reduction. Major mergers show little influence on these same properties. This suggests minor mergers regulate metallicity and contribute to over galaxy mass growth concurrently, accelerating chemical evolution post merger. Between $z=1-3$, a high $\langle \dot{M}_\mathrm{SMBH}\rangle$ leads to a reduction in $\delta \log[Z_g/Z_{\odot}]$ for all environments. Here, node galaxies show the largest reduction of approximately $0.25$ dex, suggesting that metal-rich outflows strongly drive the MZR at intermediate times. Finally, galaxies with low $M_{gas}/{M_{tot}}$ show increased $\delta \log[Z_g/Z_{\odot}]$ across all redshifts and environments, again a $0.25$ dex maximum for node galaxies. These galaxies also spike in $\delta \log[Z_g/Z_{\odot}]$ at late times, below $z=1$. At this time, galaxies in the nodes show negative $\langle \dot{M}_\mathrm{gas} \rangle$ whilst also showing the largest $\delta \log[Z_g/Z_{\odot}]$ values we observe of $0.2$ dex, suggesting the importance of the balance between gas accretion and starvation in driving MZR scatter at low redshifts.

We present a timing analysis of the high-cadence optical light curve of the high-redshift flat-spectrum radio quasar PKS\,0805$-$07 obtained during \textit{TESS} Sector~34 (MJD $59230.90$--$59239.90$). We search for short-timescale quasi-periodic oscillations (QPOs) using complementary time-series techniques, including the Lomb--Scargle periodogram (LSP) and the weighted wavelet $Z$-transform (WWZ), and evaluate their significance against red-noise variability using Monte Carlo simulations. The LSP reveals a dominant modulation at $f \approx 0.597\,\mathrm{d^{-1}}$ ($P \approx 1.7\,\mathrm{d}$) exceeding the $99.99\%$ confidence level, while the WWZ independently recovers a consistent timescale at the $\sim 99.9\%$ level and shows that the signal is temporally localized rather than persistent across the full light curve. The modulation spans $\sim 5$ coherent cycles, indicating a transient quasi-periodic feature. We discuss possible physical interpretations of the detected modulation. In a disk-based scenario, orbital motion of a hotspot near the innermost stable circular orbit implies a black hole mass of $M_{\rm BH}\sim7.2\times10^{8}\,M_\odot$, consistent with typical FSRQ values. Alternatively, magnetohydrodynamic kink instabilities in the relativistic jet can naturally produce day-scale variability for standard blazar parameters and account for the transient character of the signal. We conclude that the observed modulation is consistent with a compact, short-lived structure embedded within stochastic jet variability.

Sagittarius A* (Sgr A*) is a compact radio source at the Galactic center. Observations have confirmed that its mass is approximately (4.1)*10$^{6}$ M$_{\odot}$, and Sgr A* is generally believed to be powered by gas accretion onto a supermassive black hole. Multifrequency radio observations of the pulsar J1745-2900, about 0.12 pc away from Sgr A*, reveal an unusually large Faraday rotation. Combined with X-ray observations, this indicates that there is a strong magnetic field (greater than 8 mG) leading to a low $\beta$ plasma at large this http URL show that the gas starts to be captured by the black hole below tens of thousands of the Schwarzschild radii $r_S$, where the gas pressure starts to dominate. Assuming that the accretion rate along magnetic fields at large scales decreases with the distance to the black hole following a power law, it is shown that, with an accretion disk below tens of $r_S$, as revealed with the EHT observations, there should be a supersonic wind above such a small accretion disk, and the accretion flow may be convection-dominated from tens of $r_S$ to tens of thousands of $r_S$. Detailed modeling is warranted.

Ridges in galaxy density fields measured by photometric surveys are 2D projections of filaments in the cosmic web, and so should lens light from background galaxies. We report on a detection of this effect in Dark Energy Survey Year 3 data at high significance, though not independently of galaxy-galaxy lensing. We describe improvements to the existing subspace-constrained mean shift algorithm to locate these ridges efficiently at scale, and examine the dependence of the signal in simulations on cosmological and algorithmic parameters. We find that it depends primarily on $S_8=\sigma_8 \left( \Omega_m / 0.3 \right)^{1/2}$, and discuss improvements to our methodology that would be needed to allow precision parameter estimation.

Galaxy morphology encodes key information about formation and evolution. Large imaging surveys require automated, reproducible methods beyond visual inspection. Non--parametric indices provide an useful framework, but their performance must be quantitatively assessed. We present a homogeneous catalog of non--parametric morphological indices for DECaLS galaxies with effective radii larger than 2 arcsec. Our goal is to evaluate the reliability of indices in separating spirals and ellipticals, test their consistency with existing classification schemes, and establish their applicability for the upcoming surveys focused in the southern hemisphere. We developed galmex, a modular Python package for preprocessing images and measuring a variety of non--parametric indices. Using bona-fide spirals and ellipticals as control samples, we assessed the discriminatory power of each index, and compared them with CNN-based T-Types and Galaxy Zoo DECaLS labels. We use the indices as input for a Light Gradient Boosted Machine (LightGBM) to obtain probabilistic classifications. Concentration is the most reliable parameter from the Concentratiom + Asymmetry + Smoothness system (CAS), while asymmetry--based indices (A and S) are limited to detecting disturbed morphologies. MEGG indices (M20, Entropy, Gini, G2) provide stronger separation and trace a gradient with T--Type. By using a simple binary (0/1) label for ellipticals/spirals, classifiers trained on non--parametric indices achieve high accuracy and well--calibrated probabilities, dominated by entropy, concentration, and Gini. We release the first public catalog of CA[A_S]S+MEGG indices for DECaLS, together with galmex. We combine the non-parametric indices with machine learning framework to derive spiral/elliptical separation for galaxies below z~0.15 through a probabilistic approach.

This work presents a description of the scientific goals and objectives of OPTIMus (OPTical, Infrared, Millimeter survey of massive star-forming regions), a survey of massive star-forming regions in the optical, infrared, and millimeter wavelengths. The survey is aimed at constructing a comprehensive characterization of the multicomponent and structurally complex interstellar medium in the vicinity of young massive stars, combining both observational and theoretical aspects. Using multi-wavelength observational data, we will reconstruct the three-dimensional structure and determine the physical parameters of HII regions, photodissociation regions, and the surrounding molecular clouds. The paper describes the observational data obtained with the BTA 6-m and Zeiss-1000 telescopes of the Special Astrophysical Observatory of the Russian Academy of Sciences, the 2.5-m telescope of the Caucasian Mountain Observatory of the Sternberg Astronomical Institute of Moscow State University, and the 20-m telescope of the Onsala Space Observatory.

Gas accretion process can fuel both star formation and black hole activity, playing a critical role in galaxy evolution. The counter-rotating structures are believed to originate from gas accretion, serving as an ideal laboratory for studying its impact on galaxy evolution. Based on the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we built a sample of 147 galaxies with counter-rotating stellar disks (CRDs). This is the largest CRD sample to date, accounting for $\sim$1.5% of the MaNGA survey. For a subset of 138 CRDs, global stellar mass ($M_\ast$) and star formation rate (SFR) were measured in reference. We constructed a control sample with similar $M_\ast$ and SFR but lacking counter-rotating structures. The CRDs relatively exhibit more bulge-dominated morphology, lower molecular gas mass fraction and reside in less dense environment, supporting the hypothesis that they primarily originate from gas accretion. We classified 96 out of 138 CRDs into four types based on their stellar and gas kinematics following the criteria from Bao et al. (2022). There are two additional CRD types: 8 CRDs show misalignment between both stellar disks and gas disk, indicating multiple gas accretion events with differing angular momentum directions; 34 CRDs lack ionized gas emission, showing the highest $M_\ast$ among all the CRD types, which may represent a final stage of CRD evolution. We compared the radial gradients of gas-phase metallicity and stellar population properties between CRD types, and found that the impact of gas accretion on galaxy evolution primarily depends on the abundance of pre-existing gas in progenitors.

We present a novel gradient regularization to completely eliminate the magnetic divergence error in meshless magnetohydrodynamics (MHD), which offers a high spatial resolution and conservative advantage, due to its Lagrangian nature. Comparing with the counterpart of constrained-gradient (CG) technique, we reform $\nabla \cdot \mathbf{B}=0$ by an implicit projection method to modify the magnetic-field gradients. The accuracy of modified-gradient (MG) method is verified and it achieves exact divergence-free results with round-off precision, by using tests of shock tube, 2D and 3D vortex, magneto-rotational instability, and especially, advection experiment, compared with CG method and the GIZMO code. It leads to noticeable improvement in pattern, amplitude and numerical dissipation of divergence error of magnetic field.

At the end of the Cassini mission, Saturn's rings have been claimed to be spectacularly young compared to the age of the Solar System: their unusual ice-rich composition corresponds to initially pure ice rings polluted by interplanetary dust particles for 100 to 400 Myr. Since then, this exposure age has been commonly accepted as the real age of the rings. In this paper, we review the processes that are involved in determining the exposure age. We aim to see how the exposure age depends on various parameters and how relevant it is to define the real rings age. First, a new expression for the gravitational focusing onto planar rings, important parameter but crudely defined in the literature, is derived. Then, an analytical formula describing how the dust fraction varies with time in static or viscously evolving rings is provided, including possible vaporisation at impact. Finally, we introduce a cleaning process from space weathering to possibly alter dust and reduce its amount to make rings look younger than they are. We first found that the gravitational focusing is 5 times less important than previously thought, which automatically increases the exposure age from 0.5 to 2 Gyr. Moreover, depending on the impact properties (vaporisation rate, space weathering efficiency), several billion years can easily be reached. Finally, we find that the dust fraction in the rings converges towards a finite value which, in particular with an efficient space weathering mechanism, can be close to the observed one in the current rings. In this case, neither the age nor the initial composition of the rings can be derived, and the measure of the dust fraction and bombardment rate only constrains the physical parameters of the impacts and the efficiency of the space weathering. As long as the latter parameters are not known, the exposure age argument in favour of young rings is completely undercut.

The New Robotic Telescope (NRT) is an international collaboration to build and operate a 4 m diameter fully robotic telescope. The telescope will take advantage of the superb atmospheric conditions at the Observatory of the Roque de los Muchachos (ORM). In conjunction with a large aperture, entirely robotic operation, quick response, and a set of versatile instrumentation in the optical and near-infrared this guarantees a high scientific impact focused mainly in the area of time domain astronomy. This paper presents the scientific motivation and the status of the project, discussing possible technical solutions under evaluation for the optics, mechanics and control system.

The supernova remnant (SNR) G315.4$-$2.3 (MSH 14$-$63 or RCW 86) exhibits strong emission across the electromagnetic spectrum. Radio polarization observations probe magnetic fields and will help to understand the evolution of the SNR. We aim to investigate the radio spectrum and magnetic field properties of the SNR. We observed G315.4$-$2.3 using the Australia Telescope Compact Array (ATCA), covering the frequency range of 1.1-3.1 GHz. We performed rotation measure (RM) synthesis on the $Q$ and $U$ frequency cubes to obtain polarized intensity and RM. The regular component of the line-of-sight magnetic field was estimated from RM. The fractional polarization versus wavelength squared was used to constrain the properties of the turbulent magnetic field. We obtained image cubes of Stokes $I$, $Q$, and $U$, along with images of polarized intensity $P$, RM, and fractional polarization $p$. The radio spectra are very similar for different areas of the SNR. The foreground RM was estimated to be 55 rad m$^{-2}$, and the internal RM of most SNR areas is less than about 50 rad m$^{-2}$. The regular magnetic field along the line of sight was estimated to be about 1.4 $\mu$G in the southwest, much smaller than the total magnetic field. For most parts of the southwest and northeast, $p$ is less than 8% and is nearly constant with $\lambda^2$. We estimated the ratio of turbulent to regular magnetic field to be larger than about 3. The scale of the turbulent magnetic field for some area in the northwest might be smaller than about 0.4 pc. The radio characteristics, including spectrum and turbulent magnetic field, are very similar in the northeast and southwest, even though the evolution is quite different for these two regions based on the current models. These should be taken into account for future modeling of the evolution of the SNR.

We explore what fraction of delta sunspots in which the polarity inversion line (PIL) is sharp in photospheric magnetograms are made from a writhe kink in an emerging twisted flux rope. We searched simultaneous full-disk magnetograms and continuum images from Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO) to find 28 random sharp-PIL delta sunspots that are born well on the disk. Only one of these is made from a single newly emerged bipolar magnetic region (BMR) and therefore is a candidate for being made from a single emerging writhe-kinked flux rope. That outcome indicates that few, if any, sharp-PIL delta sunspots are made by a single emerging writhe-kinked flux rope. That is the main new finding of this paper. Each of the other 27 is made by merging of two or more emerging or emerged BMRs. We name delta-sunspot genesis from a single BMR Type I genesis. We identify another three genesis types among the other 27 delta sunspots: Type II, Type III, and Type IV. We present an observed example genesis for each of the four genesis types, and for each example present schematic drawings depicting our scenario(s) for the cause of that example genesis. The core idea of each scenario is that the delta sunspot is made by packing opposite-polarity magnetic flux together by advection into a convection downflow.

NRT is an international project to build and operate the world's largest robotic telescope. The telescope will have a segmented primary mirror with an equivalent diameter of 4 m, a set of simultaneously mounted optical and near-infrared instruments, and a response time of less than 30 seconds. The project builds on the experience gained with the successful twenty-year operation of the Liverpool telescope, and with the GTC optics and control system. All of the above together with the excellent conditions for astronomical observation of La Palma, represents a solid base and guarantees that NRT will be one of the leading facilities in the field of time domain astronomy. This contribution will analyze the current status of the project with special emphasis on the development of its optics, and the plans for its construction and operation.

Alfvén waves are widely believed to play an important role in the transport of energy from the solar photosphere to the corona through the partially ionized chromosphere. In previous work, the properties of torsional Alfvén waves were theoretically studied using a multi-fluid model. Here, we compare those multi-fluid results with those obtained using the single-fluid magnetohydrodynamic approximation, as a way to assess the performance of the latter in the context of Alfvénic waves in the lower solar atmosphere. We consider a broadband photospheric driver that excites torsional Alfvén waves with frequencies ranging from 0.1 mHz to 300 mHz. These waves propagate upwards to the corona along a magnetic flux tube expanding with height. For both models, we compare the energy flux, chromospheric reflection, transmission and absorption coefficients, and the associated heating rates. In general, the results are almost identical in the two models, with the exception of two minor differences: (1) the net energy flux reaching the corona is approximately 5% larger in the single-fluid model, mainly owing to the higher reflectivity found in the multi-fluid model for wave frequencies exceeding 10 mHz; and (2) in a narrow region around 500 km above the photosphere, the single-fluid model underestimates the plasma heating rate due to ion-neutral damping by about a factor of two compared with the multi-fluid model. Both discrepancies arise from the approximate treatment of the ion-neutral drift in the single-fluid model and are expected to have a very limited impact on practical applications.

In the excursion-set formalism, the mass distribution of primordial black holes (PBHs) is derived from the first-passage time of a random walk describing the density contrast as the coarse-graining scale varies. We address two recent criticisms that have been raised about this approach. First, it was argued that the random walks are subject to colored (i.e. correlated over time) noise, making the first-passage-time problem cumbersome. We show that this arises from an incorrect separation of drift and noise when sampling on the Hubble-crossing surface: if Fourier modes are uncorrelated, the noise is strictly white. Moreover, sampling along the Hubble-crossing surface precludes using the density dispersion as a time variable, explaining the reported pathologies. Sampling instead on a synchronous surface removes both issues. This requires solving a first-passage-time problem with a moving barrier, for which we provide an efficient numerical framework. Second, it was suggested that cloud-in-cloud (i.e. that large black holes may engulf smaller ones) is irrelevant for PBHs and that the excursion set is therefore not needed. While valid for widely separated scales, this statement fails for broad power spectra with enhanced continua of modes. We further show that Press-Schechter estimates neglecting boundary evolution can break down even without cloud-in-cloud effects. Our results establish the robustness and necessity of the excursion-set formalism in realistic PBH formation scenarios.

Surface runoff shapes planetary landscapes, but global hydrological models often lack the resolution and flexibility to simulate dynamic surface water bodies beyond Earth. Recent studies of Mars have revealed abundant geological and mineralogical evidence for past surface water, including valley networks, crater lakes, deltas and possible ocean margins dating from late Noachian to early Hesperian times. These features suggest that early Mars experienced periods allowing liquid water stability, runoff and sediment transport. To investigate where surface water could accumulate and how it may have been redistributed, we developed a global high-resolution (km-scale) surface hydrological model. The model uses a pre-computed hydrological database that maps topographic depressions, their spillover points, hierarchical connections between basins, and lake volume-area-elevation relationships. This database approach greatly accelerates simulations by avoiding repeated geomorphic processing. The model dynamically forms, grows, merges and dries lakes and putative seas without prescribing fixed coastlines, by transferring water volumes between depressions according to their storage capacities and overflow rules. We explore model behavior over the present-day Mars' topography measured by MOLA (Mars Orbiter Laser Altimeter) topography for a range of evaporation rates (from 0.1 m/yr to 10 m/yr) and total water inventories expressed as Global Equivalent Layer (from 1 mGEL to 1000 mGEL). 48 Simulations are iterated to reach the steady state. The model outputs the extent and depth of surface water bodies and identifies main drainage pathways using overflow fluxes as runoff indicators. Results show a transition toward a contiguous northern ocean between low (1-10 m) GEL values and increasing concentration of water in northern lowlands and major impact basins at higher GEL.

We present the first spatially resolved at $\sim20$ pc scale application of AGN-specific metallicity diagnostics for the nearby Compton-thick Seyfert 2 galaxy Mrk 573 ($z = 0.017$). We use Hubble Space Telescope narrow-band imaging, MUSE integral-field spectroscopy and apply AGN strong-line metallicity diagnostics based on [O III], [S II], H$\beta$, H$\alpha$, and [N II] emission lines. We construct maps of $12 + \log$(O/H) for two different metallicity calibrations and two different N/O-O/H scaling relations out to $\sim1$ kpc and down to $\sim20$ pc scales. Our analysis reveals metallicity enhancement in AGN-dominated regions, with oxygen abundances reaching up to few times Solar. The metallicity shows a patchy spatial distribution, varying on $\sim100$ pc scales, appears to trace the high Seyfert/LINER index (SLI) value regions and the VLA 6 cm jet/radio lobe emission. These spatial correspondences and the lack of evidence for star formation in the bicone region suggest that the enrichment originates from metals transported from the nuclear AGN regions by winds, outflows, or jets.

Since 2019, three scintillator detectors of the EEE collaboration have been continuously measuring cosmic muon rates at 78.9°N at the Ny-Ålesund Research Station (Svalbard). The resulting six-year time series reveals a pronounced annual modulation, driven primarily by seasonal atmospheric variations. Utilizing routine radiosonde profiles collected above the same site, we applied several established techniques --along with a tailored analysis approach-- to investigate the relationship between muon rate and atmospheric temperature. The temperature-corrected muon-rates are analysed using the Lomb-Scargle periodogram technique in order to investigate the presence of remaining periodic structures. Finally, the temperature corrections coefficients of our analysis are compared with measurements in other stations located at lower latitudes.

Wave-optics effects in strongly lensed gravitational waves (GWs) provide a new interferometric probe of dark matter substructure. We compute the full diffraction integral for GWs propagating through statistically generated cold dark matter subhalo populations and quantify the resulting frequency-dependent amplification in the Laser Interferometer Space Antenna (LISA) band. We show that realistic galaxy-scale lenses generically produce percent-level amplitude and phase distortions in strongly magnified images, primarily induced by subhalos in the mass range $10^4$-$10^7\,M_{\odot}$. These signatures arise naturally within the standard cold dark matter paradigm and should be detectable in high signal-to-noise LISA events. Strongly lensed GWs thus offer a direct and complementary window on dark matter structure at subgalactic mass scales inaccessible to electromagnetic measurements.

We present a novel method for the tomographic reconstruction of the bias-weighted mean electron pressure $\langle bP_e \rangle$ and star formation rate density $\langle b\rho_{\mathrm{SFR}} \rangle$, by simultaneously modelling the contribution from the thermal Sunyaev-Zel'dovich (tSZ) effect and the Cosmic Infrared Background (CIB) to the cross-correlation between photometric galaxy samples and multi-frequency Cosmic Microwave Background (CMB) maps. The resulting measurements are independent of the galaxy clustering properties and robust against cross-contamination between tSZ and CIB. Applying this method to publicly available data, we reconstruct the cosmic evolution of $\langle bP_e \rangle$ and $\langle b\rho_{\mathrm{SFR}} \rangle$ out to $z\sim1$, making our measurements publicly available. Our measurements of both quantities are broadly compatible with predictions from the fiducial FLAMINGO hydrodynamical simulation, although we observe a lower gas pressure at low redshifts, in agreement with other measurements.

Strong gravitational lensing provides an independent and powerful probe of cosmic expansion by directly linking observables to cosmological distances. Upcoming surveys such as LSST will discover large number of galaxy-galaxy strong lensing systems, offering a new route to precise cosmological constraints. In this paper, we propose a Fisher-like sensitivity factor to map how the cosmological information of strong-lensing distances changes across the lens-source redshift plane. Applying such factor to the distance ratio $D_{ls}/D_s$, the time-delay distance $D_{\Delta t}$, and the double-source-plane ratio, we determine the ``sensitivity valleys'' where an observable becomes insensitive to a given parameter. The realistically simulated LSST lens population, which largely lies outside the distance-ratio valleys, covers the most sensitive region for $(w_0,w_a)$ parameter space. We then develop a new hierarchical framework, which could calibrate the redshift evolution of lens mass-density slopes and constrain cosmological parameters simultaneously. Focusing on the LSST mock data, we demonstrate that ignoring mass-profile evolution can bias $\Omega_m$ by up to $\sim 10\sigma$, while modeling the lens evolution could perfectly recovers the fiducial cosmology and yield stringent cosmological constraints (e.g., $\Delta\Omega_m \simeq 0.01$ and $\Delta w \simeq 0.1$ for $\sim 10^4$ lenses).

We present an alternative measurement of the Baryon Acoustic Oscillation (BAO) using Legendre multipole representation of the Ly$\alpha$ forest correlation functions from the second data release (DR2) of the Dark Energy Spectroscopic Instrument survey. Compressing the auto- and cross-correlation functions into Legendre multipoles yields a positive-definite covariance matrix without any smoothing -- unlike the baseline DR2 analysis -- thanks to a significantly reduced data vector size. We introduce the statistical corrections required to debias the finite-sample covariance matrix estimate and demonstrate that monopole and quadrupole terms for both auto- and cross-correlations can be used even when the correlation functions are distorted by continuum errors and contaminated by metals. This formalism has slightly diminished the constraining power of the BAO scale, while considerably weakening constraints on nuisance parameters. We measure the isotropic BAO scale with $0.96\%$ precision at $z_\mathrm{eff}=2.35$, the Hubble parameter $H(z_\mathrm{eff})=(238.7\pm3.4)~(147.09~\mathrm{Mpc}/r_d) ~\mathrm{km~s}^{-1}~\text{Mpc}^{-1}$, and the transverse comoving distance $D_M(z_\mathrm{eff})=(5.79 \pm 0.10)~(r_d/147.09~\mathrm{Mpc})$~Gpc for a given value of the sound horizon ($r_d$). Our BAO results are entirely consistent with the baseline DR2 analysis.

We present near-infrared (NIR) observations of the 2021 eruption of the recurrent nova RS Ophiuchi. The dataset provides both pre- and post-eruption perspectives on the eruption, as well as NIR spectra at high cadence. The spectrum obtained in 2020 June (14.3 years after the 2006 eruption, and 428.1 days before the 2021 eruption), is that of the red giant secondary, on which are superimposed several emission lines which most likely arise in the red giant wind. Spectra obtained during the eruption consist of emission (including coronal) lines, superimposed on a bremsstrahlung continuum at 8900K. The temperature of the coronal gas is estimated to be $10^{6.0}$K on day 11.7, and $10^{5.9}$K on day 31.7. The high cadence observations, obtained on day 31.7 of the eruption, provide no conclusive evidence for rapid ($<\sim1$~minute) variations in the HeI 1.0833$\,\mu$m line. Data obtained about one year after the eruption show that there may have been changes in the spectral type of the secondary after the 2021 eruption.

Small-scale vortical motions in the upper solar atmosphere are abundant and occupy about 2.8% of the photosphere at any given time. Although considerable work has focused on the detection and analysis of individual solar vortices, the interconnected and multi-scale behaviour of these coherent structures remains largely unexplored. We present a methodology for studying this behaviour through vortex interactions, to improve our understanding of how small- and large-scale photospheric flows contribute to energy transfer into the upper solar atmosphere and to the driving of solar activity. We represent vortices as a network of interacting structures. We apply a community detection algorithm to derive an optimal reduced network composed of highly interconnected vortex groups. From the interaction patterns and group structure, we define three roles within each community: peripheral, connector and hub. We then track both vortex communities and their member vortices from the photosphere into the chromosphere and across their lifetimes. On average, vortices assigned to these roles persist to greater heights in the chromosphere and have longer lifetimes than unclassified vortices. This shows that community detection can identify vortices with greater dynamical influence on the upper solar atmosphere. We also find that 32% to 58.6% of vortex communities exhibit global periodic behaviour following a helical path. This collective vortical motion may indicate an enhanced mechanism for wave excitation. Solar vortical community detection, therefore, offers a new framework for studying solar vortices and a new perspective on the importance of collective vortex dynamics.

Methanol, the simplest complex organic molecule found in space, is considered a key compound necessary for the formation of chemical species of prebiotic interest. Methanol detections in protoplanetary disks remain scarce, even though it is frequently detected in the material surrounding other Young Stellar Objects. We investigate the presence of methanol in the protoplanetary disk around the HL Tau protostar, motivated by the detection of spatially resolved warm water emission. Given the similar volatility of methanol and water, thermally desorbed gas-phase methanol is expected to emit from the same region of the HL Tau disk where water vapour has been observed. Accordingly, we selected and imaged the most promising ALMA archival observations to search for rotational methanol lines. We found no methanol emission in the analysed archival datasets. Assuming optically thin emission and LTE, we derive stringent upper limits on the methanol column density for different excitation temperatures: < 7.2 x 10^(14) cm^(-2) at 100 K and < 1.8 x 10^(15) cm^(-2) at 200 K, assuming a circular emitting region with a radius of 17 au (~ 0.12''). Furthermore, we obtain a stringent upper limit on the methanol-to-water column density ratio (< 0.55 x 10^(-3) at 100 K and < 1.4x 10^(-3) at 200 K), which is, on average, an order of magnitude lower than the values measured for other Young Stellar Objects and Solar System comets. We argue that the most likely explanation for the methanol non-detection in HL Tau is the presence of optically thick dust in the central region of the disk, which obscures part of the methanol emission. The upper limit on the methanol-to-water ratio in the HL Tau disk is at least an order of magnitude smaller than most clouds, YSOs and comets, possibly due to radiative transfer and/or excitation effects, or due to a different chemical evolution compared to the other sources.

We present an automated framework to search for optical counterparts of LIGO-Virgo-KAGRA (LVK) gravitational wave (GW) superevents using public Zwicky Transient Facility (ZTF) alerts processed through the ALeRCE broker. The goal is to filter and identify optical transients potentially associated with binary black hole (BBH) mergers during the LVK O4a and O4b observing runs. Using the Automatic Learning for the Rapid Classification of Events (ALeRCE) infrastructure, we spatially query ZTF alerts within GW localization regions and apply machine learning classifiers, host-galaxy crossmatching, and temporal cuts within 200 days post-merger to isolate plausible candidates. Our search yielded one candidate in O4a and four in O4b, several consistent with the supernova or tidal disruption event regime. This work demonstrates that public alert brokers can establish a robust baseline for systematic searches for electromagnetic counterparts to GW superevents in current and future observing runs. Our algorithm provides a systematic approach to search for BBH counterparts for all significant LVK GW superevents using survey telescope alerts. The search, together with the accompanying analysis, demonstrates the significance of the counterpart candidates, with one candidate ultimately identified as a transient event consistent with a Bowen fluorescence flare in a now discarded active galactic nucleus (AGN).

The large population of broad-line Active Galactic Nuclei (AGN) observed with the James Webb Space Telescope (JWST) at $z \gtrsim 4$ opens a new window onto the black hole-galaxy connection in the first Gyr of cosmic history. We use the JADES survey-level dataset and develop a forward-modeling Bayesian framework that explicitly accounts for broad H$\alpha$ detectability, ensuring that selection effects are incorporated into the likelihood function. With this approach, we constrain the black hole-stellar mass ($M_{\mathrm{BH}}$-$M_\star$) relation to be $\log M_{\rm BH} = -4.06^{+0.50}_{-0.51} + 1.17^{+0.06}_{-0.06}\,\log M_\star$, with an intrinsic orthogonal scatter of $\sigma_{\rm int} = 0.63^{+0.14}_{-0.11}$ dex. The slope and normalization are consistent with local determinations, indicating that the average scaling was already established by $z \sim 4$-6. This suggests that the primary evolution of the relation occurs in its dispersion rather than in its mean normalization. In contrast, the substantially larger intrinsic scatter relative to the nearby Universe reveals a wider diversity of black hole-galaxy growth histories, likely driven by bursty accretion, delayed feedback, and differences in merger or seeding histories. Future JWST samples will be crucial to test whether this increased scatter is a persistent feature of the high-redshift Universe.

The study of Stark broadening of neutral helium lines, despite significant advances over recent decades, has not led to updated large grids of helium line profiles relevant to the spectroscopic study of helium-rich stars. While the semi-analytical approach based on the standard Stark broadening theory is efficient for generating such grids, it presents challenges in incorporating additional physical effects into the model. Motivated by recent studies that highlight potential issues with line profiles in the context of white dwarf stars, this paper leverages advances in computer simulations to create a new grid of line profiles for 13 neutral helium lines in the optical range. These profiles cover densities ranging from 10^14 to 6 x 10^17 cm^-3 and temperatures from 10,000 K to 40,000 K, with the exception of the narrower He I 4713 line, for which the profile grid begins at 10^15.5 cm^-3. The primary goal of this research is to present the new grid and compare it with both the semi-analytical approach and other simulation results. By doing so, corrections to the previous grid will be explored, providing a foundation for future studies that utilize this updated grid. We also examine the impact of these new profiles on the determination of physical parameters for a range of astrophysical objects, including DB white dwarfs and other helium-rich stars.

We present a novel method for constructing high-accuracy, time-varying disk-integrated stellar absorption line profiles that isolate the effects of granulation alone. This framework provides an effectively unlimited supply of physically consistent training data, offering a unique opportunity to study granulation-driven velocity variability with no contamination from other stellar processes or instrumental systematics. Our interpolation scheme enables accurate profile generation at arbitrary limb angles and successfully reproduces observed disk integrated solar bisector shapes from IAG spectra. Using four Fe I lines (525.0, 615.2, 617.3, and 627.1 nm), we produce 1000 model star disk-integrated realisations per line and find an isolated granulation-induced RV scatter of 0.16-0.21 m s^-1. Using our synthetic profiles and assuming infinite signal-to-noise, we find strong correlations between various line-shape metrics and convective blueshift, demonstrating that line-shape diagnostics can, in principle, trace granulation effects. Equivalent width proves the strongest diagnostic, achieving up to 60% scatter reduction. However, the strength of all simple line shape diagnostics rapidly diminishes once photon noise is injected. Even when artificially boosting the signal to represent a spectrum containing ~1000 spectral lines, the achievable improvement with these metrics remains below 10% at typical signal-to-noise ratios. Our results highlight the need for more robust, noise-resilient diagnostics and position our synthetic dataset as a valuable testbed for developing and benchmarking such methods.

The discovery rate of optical transients will explode to 10 million public alerts per night once the Vera C. Rubin Observatory's Legacy Survey of Space and Time comes online, overwhelming the traditional physics-based inference pipelines. A continuous-time forecasting AI model is of interest because it can deliver millisecond-scale inference for thousands of objects per day, whereas legacy MCMC codes need hours per object. In this paper, we propose SELDON, a new continuous-time variational autoencoder for panels of sparse and irregularly time-sampled (gappy) astrophysical light curves that are nonstationary, heteroscedastic, and inherently dependent. SELDON combines a masked GRU-ODE encoder with a latent neural ODE propagator and an interpretable Gaussian-basis decoder. The encoder learns to summarize panels of imbalanced and correlated data even when only a handful of points are observed. The neural ODE then integrates this hidden state forward in continuous time, extrapolating to future unseen epochs. This extrapolated time series is further encoded by deep sets to a latent distribution that is decoded to a weighted sum of Gaussian basis functions, the parameters of which are physically meaningful. Such parameters (e.g., rise time, decay rate, peak flux) directly drive downstream prioritization of spectroscopic follow-up for astrophysical surveys. Beyond astronomy, the architecture of SELDON offers a generic recipe for interpretable and continuous-time sequence modeling in any time domain where data are multivariate, sparse, heteroscedastic, and irregularly spaced.

A complete picture of the gas thermodynamics around galaxies is imprinted on the cosmic microwave background (CMB). Indeed, the thermal, kinematic, and relativistic Sunyaev-Zel'dovich effects (tSZ, kSZ, rSZ) measure the gas density, temperature, pressure of baryonic feedback and bulk velocity around galaxies, along with the gravitational potential it sits in. This full thermodynamic picture promises to constrain galaxy formation models and gas related uncertainties in the impact on galaxy lensing. Recent kSZ measurements around galaxies suggest that the gas may be more extended than anticipated, pointing to powerful feedback processes and large baryonic corrections to lensing. How robust are these conclusions about the galaxy-halo connection, including satellite fraction and high-mass outliers, or to 2D projection effects and large-scale velocity modes? In this paper, we give estimates for these effects using a simulated sample of DESI-like luminous red galaxies within the IllustrisTNG hydrodynamical simulation and the Abacus N-body simulation. We show that analyzing projected 2D profiles can lead to biases when computing quantities like the gas fraction. We also find that in the absence of spatial filtering, the 2-halo term is non-negligible for kSZ even at the smaller radii where the 1-halo term dominates. We show that a 1% uncertainty in the satellite fraction of galaxies can propagate into uncertainties of $\pm$1%, $\pm$3% and $\pm$5% in the 1-halo terms of the kSZ, tSZ, and rSZ signals, respectively. We show that masking the 2% most massive objects in the sample reduces the profile amplitudes by up to 10%, 40%, and 75% for the kSZ, tSZ, and rSZ signals, respectively. Finally, we show that naive simulations of the kSZ effect can be biased by an artificial Doppler term, which is automatically removed when high-pass or compensated aperture filtering is applied.

Among the most compelling cold dark matter candidates, the axion has recently been subject to a wide range of astrophysical studies aiming to constraints its properties. We present updated bounds on the isocurvature fraction, $f_{\rm{iso}}$, which parameterizes the contribution of isocurvature perturbations induced by post-inflationary produced axion-like particles (ALPs) to the ordinary power spectrum. We use new simulations based on the Sherwood-Relics suite to fit high-resolution Lyman-$\alpha$ forest flux power spectrum data. With the published noise model of the Lyman-$\alpha$ forest data, we find a tentative detection of $f_{\rm{iso}}$ = ${0.0064^{+0.0012}_{-0.0014}}$ (68% C.L), after accounting for the degenerate effect of IGM thermal evolution. With a more conservative modelling of the residual noise in the data, the upper bound is weakened to $f_{\rm{iso}}< 0.0084$ (95% C.L), which translates into an ALP temperature-independent mass $m_a > 1.73 \times 10^{-18}$eV. Our constraints are stronger than bounds derived from large-scale structure probes at higher and lower redshifts and are competitive with those derived from UV luminosity function data. Interestingly, the best current Lyman-$\alpha$ forest data prefers a non-zero contribution from isocurvature modes.

Carbonaceous Chondrites have special significance in the stellar evolution and in particular in the evolution of life on earth. The carbonaceous meteorite that fell in Mukundpura village, Jaipur, Rajasthan on 6th June 2017 is one such rare CM2 (Carbonaceous Chondrite) carbonaceous meteorite. We carried out high resolution scanning and transmission electron microscopic (TEM) studies on typical thin sections, showing abundant grains of iridium (Ir), pentlandite (NiS), and more interestingly crystalline carbon (C). These crystallite carbon grains resemble nanodiamond like signature in the freshest Mukundpura meteorite. The high-resolution Raman spectroscopic measurements are carried out on the crystalline carbon grains, showing well resolved three distinct peaks with a vibrational mode at 1315 cm-1, with the onset of a weak vibrational mode at 1150 cm^-1, substantiating the observation of nanocrystalline diamond in Mukundpura meteorite. The broad peak centered at 1360 cm^-1 and 1575 cm^-1 (as an average), suggest the presence of graphitic carbon as well together with apparent presence of nanocrystalline diamond. The average size of nanocrystalline diamond is ~ 3-5 nm. High iridium content in this meteorite supports the meteoric impact related iridium anomaly in geological stratigraphic boundaries (this http URL-Tertiary boundary) that has caused mass extinction of flora and fauna.

Because of the special angular distribution of excited electrons by the photoelectric effect, the Gas Pixel Detector (GPD) is effective in measuring keV X-ray polarization of astrophysical events (e.g. gamma-ray bursts), by capturing ionization tracks of excited electrons as polarized images. Traditionally, the emission angles of photoelectrons are extracted from polarized images first, and statistics are then performed on these angles to infer the polarization direction and intensity. However, observation with the wide field of view requires the incident angle of X-rays not directly attainable through the traditional analysis process. In this paper, we propose using the generalized sliced Wasserstein (GSW) distance, projected by neural networks with random weights, as a completely data-driven approach to analyze X-ray polarization based on two-dimensional polarized images. We find the structures of the randomized neural networks matter when focusing on different aspects of the polarized images, and take advantage of the discrimination abilities by different neural network structures. The proposed method, named the structured GSW distance, successfully distinguishes polarized images with different configurations of incident angles and polarization directions. Furthermore, we build a simplified statistical model based on the von Mises distribution and the circular Wasserstein distance and compare the model against the proposed method, showing their high consistency. The computational method reported in this paper may benefit GPD-based polarimetry in astroparticle experiments and also pattern analysis on raw data from pixel detectors.

We will address the existence of a new symmetry for an imperfect fluid by introducing local four-velocity gauge-like transformations for the case when there is vorticity. A similar tetrad formulation as to the Einstein-Maxwell spacetimes formalism presented in previous manuscripts will be developed in this manuscript for the imperfect fluids. The four-velocity curl and the metric tensor will be invariant under these kind of four-velocity gauge-like local transformations. While the Einstein-Maxwell stress-energy tensor is locally gauge invariant under electromagnetic gauge transformations, the perfect fluid stress-energy tensor will not be invariant under four-velocity gauge-like local transformations. We will dedicate our analysis to the imperfect fluid stress-energy tensor that will be invariant under local four-velocity gauge-like transformations when additional transformations are introduced for several variables included in the stress-energy tensor itself. We will also pay special attention to the construction of a vorticity stress-energy tensor invariant under local four-velocity gauge-like transformations. An application on neutron stars will be developed in order to show the simplifications brought about by these new tetrads.

The non-linearly coupled Einstein-matter field equations predict the existence of shadows with well-defined boundaries around black holes. We prove that, in spherically symmetric hairy black-hole spacetimes whose matter fields satisfy the weak energy condition, the radii of these shadows are bounded from below by the dimensionless relation $r_{\text{sh}}/r_{\text{H}}\geq 3\sqrt{3}/2$, where $r_{\text{H}}$ is the horizon radius of the central hairy black hole. The characteristic shadow of the (bald) Schwarzschild black-hole spacetime saturates the analytically derived lower bound.

Typical gravitational production of relics from amplification of inflationary perturbations assumes Bunch-Davies initial conditions, i.e. a vacuum with initially no particles. In this paper we investigate the impact of non Bunch-Davies initial conditions to the final abundance of relics, with particular attention to the parameter space where the total dark matter abundance is reproduced. We present a general framework for any initial condition, through which we show their non-trivial effect on both spectrum and late-time abundance. We argue that for particles whose source of conformal symmetry breaking comes only from a mass term (spin-1/2 fermions and conformally coupled scalars), the choice of initial conditions has little impact on the mass range relevant to dark matter. For other particles, e.g. the longitudinal mode of spin-1, we see a large deviation from the standard computation. We exemplify and quantify our results with an initial thermal state and a two-stage inflation scenario, highlighting that the total dark matter can be obtained for a wide range of masses.

Dark matter direct detection experiments involving electron recoils are beginning to test highly-predictive, thermal-relic milestones for sub-GeV dark matter models. Due to the Lee-Weinberg bound, thermal dark matter candidates in this mass range necessarily require comparably-light mediator particles to achieve a suitably large annihilation cross section. Here we present new thermal-relic milestones for sub-GeV dark matter candidates that couple to vector mediators. In these models, the mediators are massive gauge bosons of anomaly-free abelian extensions to the Standard Model, including the dark photon, gauged $L_i - L_j, B-L$, and $B-3L_i$ models, where $B$ is the baryon number, $L$ is the lepton number, and $i,j$ index the lepton families. Since the same interactions that govern cosmological production also govern electron scattering, the targets we present are firmly predictive and allow for these models to be robustly discovered or falsified. Furthermore, since the mediators we study exhaust the minimal anomaly-free U(1) extensions to the Standard Model, our results offer a complete list of predictive milestones for sub-GeV dark matter coupled to vector mediators.

Low-thrust electric propulsion missions are often designed under simplifying assumptions such as constant thrust or fixed specific impulse, neglecting the strong coupling between trajectory dynamics, spacecraft power availability, and propulsion performance. In deep-space environments with reduced solar irradiance, these assumptions can lead to suboptimal or infeasible designs, underscoring the need to simultaneously optimize the trajectory and power subsystem. This paper presents a multidisciplinary design optimization (MDO) framework for the simultaneous design of low-thrust trajectories and spacecraft power systems, with explicit coupling to electric propulsion performance. The framework incorporates a high-fidelity variable-specific impulse model of the SPT-140 Hall thruster, in which thrust and efficiency are directly constrained by time-varying solar power availability and solar array degradation, rather than treated as fixed parameters. The coupled problem is posed as a time-optimal control problem and addressed using a framework built on top of OpenMDAO and Dymos toolchains, where Dymos employs a collocation-based direct-transcription approach for trajectory optimization. OpenMDAO provides accurate analytic partial derivatives, enabling efficient gradient-based optimization. A Fast Fourier Series shape-based method is used to generate dynamically feasible initial guess trajectories, and the resulting nonlinear programming problem is solved using IPOPT. The proposed framework is demonstrated through a low-thrust orbit insertion scenario around asteroid 16-Psyche, a regime in which reduced solar irradiance makes power-aware trajectory design particularly critical. Simulation results demonstrate the framework's ability to capture key power-propulsion-trajectory trade-offs, highlighting the importance of integrated power optimization for realistic electric propulsion mission design.

Recently, an exact rotating black hole solution in a parity-violating theory of gravity was obtained via a conformal transformation of the Kerr solution in general relativity, with parity-violating effects encoded in the conformal factor. We study the quasinormal modes (QNMs) of a test scalar field minimally coupled to gravity on this conformal Kerr background, treating the parity-violating effects perturbatively while allowing for arbitrary black hole spin, from the non-rotating case to the near-extremal regime. For low spin, we derive a perturbative formula for the QNM frequencies that includes the leading-order parity-violating correction. For high spin, particularly in the near-extremal regime, we find sizable deviations from the Kerr QNM frequencies. Our results point to a new avenue for probing parity-violating physics in the strong-gravity regime through black hole QNMs.

We investigate the principles of quantum field theory using a stiff de Sitter space. We demonstrate that a non-unitary Lagrangian on a Euclidean AdS geometry can produce the perturbative expansion of late-time correlation functions to all orders. This discovery greatly simplifies perturbative computations while also allowing us to prove fundamental features of these correlators, which are part of a Euclidean CFT. This allows us to construct an OPE expansion, limit the operator spectrum, and deduce the analytic structure of the spectral density that captures the conformal partial wave expansion of a late-time four-point function. In general, the standard CFT concept of unitarity does not apply to dimensions and OPE coefficients. Rather, the positivity of the spectral density represents the unitarity of the de Sitter theory. This assertion is non-perturbative and does not depend on the use of Euclidean AdS Lagrangians. In a scalar theory, we compute tree-level and entire one-loop-resummed exchange diagrams to demonstrate and verify these characteristics. In the spectrum density, an exchanged particle shows up as a resonant characteristic that may be helpful in experimental searches.

Modeling inviscid compressible flows with shocks and vortex dominated dynamics remains challenging for particle grid methods due to moving discontinuities, cell crossing noise, and quadrature degradation under strong deformation. Building on a FLIP/APIC framework with vorticity aware tensor artificial viscosity, we identify a long time RTI failure mode: particle depletion at spike heads degrades quadrature and particle grid coupling, producing nonphysical, void-like dents. Standard mitigations (CPDI lite and subcell-jittered seeding) reduce but do not eliminate this artifact. We therefore add two sampling-aware controls: (i) conservative split resampling that replenishes depleted cells while exactly conserving mass, momentum, and internal energy; and (ii) a soft-switch that attenuates only the APIC affine term when local support is insufficient. Tests on the Sod shock tube and single/multi mode RTI show that the method removes spike head voids in long-time RTI while preserving vortex roll up, and matches reference Euler growth metrics

A new phenomenological dark energy model, originally associated to the large-scale structure formation and considered as a solution to the fine-tuning and coincidence problems related to the cosmological constant, was analyzed within the framework of General Relativity in a Friedman-Robertson-Walker spacetime and its model parameters were estimated using cosmic chronometers and recent DESI data. It turns out that the proposed model can serve as an alternative evolving dark energy model with a novel equation of state function, apart from other popular propositions in the literature. Due to the form of this phenomenological energy density ansatz, which starts to rise with the nonlinear structure growth in the universe and falls with the domination of cosmic voids, we prefer to call it structure-induced dark energy. Observational constraints show that it is not only a suitable solution for the fundamental problems such as coincidence or fine-tuning problems, it gives flexibility, when considering the cosmic tensions and presents a new perspective on the evolving dark energy models.

In this letter, we explore the phenomenological impact of inflationary gravitational particle production in the physics of Dark Matter (DM). Large-scale DM fluctuations generated during inflation behave as gravitational particles upon their post-inflationary horizon reentry and alter the conventional Boltzmann dynamics of DM with a non-conserving source term, thereby producing significant phenomenological consequences. Within this framework, we analyze four distinct types of DM classified according to their production mechanisms. Dark matter may be completely non-interacting with the thermal bath, behaving as Inert Dark Matter. Alternatively, depending on the strength of its interactions with bath particles, DM may exhibit WIMPy, UFOy, or FIMPy behavior, sharing characteristics with their conventional counterparts. The late-time enhancement of the DM number density, driven by the successive horizon reentry of gravitationally produced low-momentum modes, enlarges the viable parameter space for both thermal and non-thermal DM scenarios. Remarkably, this expanded parameter space remains consistent with current constraints from $\Delta N_{\rm eff}$ and Lyman-$\alpha$ bound.

We study slow-roll inflation driven by a scalar field non-minimally coupled to gravity within the effective framework of Loop Quantum Cosmology (LQC), including inverse-volume corrections. We consider two physically motivated classes of potentials, a Higgs-like quartic potential $V\propto\phi^{4}$ and string-inspired fractional monomial potentials $V\propto\phi^{p}$ with $p<1$. Working at first order in the slow-roll expansion, we derive analytic expressions for the inflationary observables, namely the scalar spectral index $n_s$, the tensor-to-scalar ratio $r$, and the running $\alpha_s\equiv dn_s/d\ln k$, and then solve the corrected background dynamics numerically to obtain quantitative predictions. Confronting these results with current observational constraints from Planck 2018 and ACT DR6, we find that the model can lie within the allowed region of the $(n_s,r,\alpha_s)$ parameter space, including a mild preference for slightly larger $n_s$, as suggested by recent ground-based measurements. We also compute the probability of achieving sufficient slow-roll inflation in this setting. Although effective LQC replaces the initial singularity with a nonsingular quantum bounce, the likelihood of a sufficiently long inflationary phase depends on the pre-inflationary dynamics and on the inflaton potential. Using the canonical Liouville measure on the effective phase space, we determine the fraction of post-bounce trajectories that yield sufficient inflation and find that the non-minimal coupling parameter $\xi$ substantially enlarges the phase-space volume of favorable initial conditions relative to the minimally coupled case, exhibiting an attractor-like enhancement that saturates at large $\xi$.

Two-dimensional particle-in-cell simulations with a realistic mass ratio reveal the generation mechanisms of the out-of-plane magnetic field in magnetic islands/flux ropes during magnetic reconnection. In the absence of an initial guide field, reconnection produces a large electron temperature anisotropy (around 4.5) inside magnetic islands that drives the Weibel instability. Strong out-of-plane magnetic fields (Bz/B0 around 0.4, greatly exceeding the Hall field) with a regular bipolar structure grow inside islands. A space-time analysis reveals a one-to-one correspondence between the temperature anisotropy and the development of the Weibel magnetic field. The instability relaxes the anisotropy, but island merging leads to anisotropy reemergence and re-excitation. In the presence of a strong ambient guide field (Bg/B0 = 0.5), the electron outflow from the X-point deflects along the separatrices and forms a circular current loop wrapping the flux ropes. This flux-rope separatrix current generates an out-ofplane magnetic field that reinforces the ambient guide field, reaching Bz/B0 around 1.4. The current can, in some cases, drive the electron Kelvin-Helmholtz instability, which produces electron vortices and strengthens the magnetic field. Mergers significantly broaden the islands and further strengthen the field. These self-generated out-of-plane magnetic fields scatter electrons and reduce their temperature anisotropy, which can potentially affect electron heating via Fermi reflection. The simulation results are supported by spacecraft observations suggesting that ambient guide fields can be enhanced within flux ropes in Earth's magnetotail.

Hermiticity is usually treated as a foundational axiom of quantum mechanics, guaranteeing real spectra and unitary time evolution. In this work we argue that Hermiticity is more naturally understood as a symmetry law arising from the global conservation of an inner product current. We show that in spacetimes admitting complete Cauchy surfaces without boundary flux this conservation reduces to the familiar Hermiticity condition of the canonical inner product. However, in the presence of causal horizons, most strikingly in black hole geometries, this conservation law becomes obstructed for restricted observers. Tracing over inaccessible degrees of freedom then inevitably yields completely positive trace preserving dynamics with an effective non-Hermitian generator. Using quantum thermodynamics and the monotonicity of relative entropy, we demonstrate that the generalized second law may be reinterpreted as an entropy balance that compensates precisely for the flux of inner product charge through the horizon. The structure of Einstein equations, through the Bianchi identity and the Raychaudhuri focusing equation, provides the geometric mechanism underlying this balance. We also show that black hole ringdown can serve as a realistic observational probe of this idea and may provide quantitative upper bounds on the strength of horizon-induced inner product flux. In this way gravity, entropy production, and effective non-Hermiticity are unified under a single structural principle, with Hermiticity emerging as the special case of globally conserved inner product symmetry.