Papers for Tuesday, Jul 08 2025

White dwarfs represent the most common end stage of stellar evolution and are important for a range of astrophysical questions. The high-resolution ultraviolet spectroscopic capability of the Habitable World Observatory (HWO) offers a unique capability to characterize white dwarfs. In this documents, we focus on two specific science cases for HWO -- white dwarfs as probes of extrasolar planet compositions, and fundamental astrophysics. HWO will have the sensitivity to measure a suite of heavy elements, such as S, C, O, Fe, and Si, in a large sample of polluted white dwarfs to constrain the water content and the light elements in the cores of extrasolar planets. HWO can also be used to search for any small variation on the fine structure constant in the presence of strong gravity. Both science cases require a minimum resolving power of 60,000, and a ultraviolet coverage down to at least 900 Angstrom.

The Habitable Worlds Observatory (HWO) offers a unique opportunity to revolutionize our understanding of planetary formation and evolution. The goal of this Science Case Development Document (SCDD) is to investigate the physical and chemical processes that shape the composition and atmospheric mass loss in exoplanets. We review the key observables currently known as diagnostics of mass loss via transit observations, i.e., absorption lines of escaping hydrogen (Lyman-alpha), helium, and metals (Fe, Mg, C, O). We also explore the challenges to infer planetary formation processes based on atmospheric composition characterization. HWO could enable a broad, continuous coverage from far-ultraviolet to near-infrared spectroscopy (~100--1600 nm) at high resolution (R > 60, 000), which is essential to make these measurements, disentangle their planetary origin from stellar activity, and ultimately, contextualize the escape rates by simultaneously characterizing the composition, cloud predominance, and thermal structure of exoplanet atmospheres.

A detailed analysis of the structural, astrophysical, kinematic, and dynamical properties of the open clusters Roslund 3 and Ruprecht1 74 is carried out using CCD UBV photometry in conjunction with astrometric and photometric data from Gaia DR3. Membership probabilities were computed via the UPMASK algorithm applied to Gaia proper motions and trigonometric parallaxes, leading to the identification of 198 likely members for Roslund 3 and 397 for Ruprecht 174. Astrophysical parameters were derived using both the classical approach, where parameters are independently determined, and a MCMC technique, which estimates them simultaneously. The agreement between the results from both methods confirms their reliability and highlights the robustness of the classical method. The reddening values were determined as $E(B-V)=0.410\pm 0.046$ mag for Roslund 3 and $E(B-V)=0.615\pm 0.042$ mag for Ruprecht 174. The estimated distances are $d=1687 \pm 121$ pc for Roslund 3 and $d=2385 \pm 163$ pc for Ruprecht 174. Both clusters exhibit metallicities close to the solar value, with [Fe/H] = $0.030 \pm 0.065$ dex for Roslund 3 and [Fe/H] = $0.041 \pm 0.064$ dex for Ruprecht 174. The corresponding ages were found to be $\tau=60\pm 6$ and $\tau=520\pm 50$ Myr, respectively. The present day mass function slopes were found to be $1.18 \pm 0.13$ for Roslund 3 and $1.53 \pm 0.30$ for Ruprecht 174, consistent with the canonical Salpeter value within uncertainties. Galactic orbital analyses indicate that both clusters are thin-disk members confined within the Solar circle. Additionally, relaxation times and spatial distributions of stars suggest that both clusters have reached dynamical relaxation and exhibit clear signs of mass segregation.

Planets with large bodies of water on their surface will have more temperate and stable climates, and such planets are the ideal places for life-as-we-know-it to arise and evolve. A key science case for the Habitable Worlds Observatory (HWO) is to determine which planets host surface liquid water. Aside from its implications for planetary climate and astrobiology, detecting surface water on terrestrial exoplanets would place important constraints on our theories of planet formation and volatile delivery. Rotational variability in the reflectance of an exoplanet may reveal surface features rotating in and out of view, including oceans. Orbital changes in reflectance and polarization, meanwhile, are sensitive to the scattering phase function of the planetary surface, including specular reflection from large bodies of water. Although these techniques are applicable to all temperate terrestrial exoplanets, we focus in this document on the directly-imaged planets that are more likely to drive the HWO coronagraph design. Identification of water oceans relies on detecting a liquid, and using other lines of evidence to narrow that liquid down to being water. Liquids have smoother surfaces than most solids, and hence exhibit specular reflection instead of diffuse reflection. In practice, this makes lakes and oceans look dark from most illumination angles, but mirror-like at glancing angles. HWO is uniquely capable of identifying surface liquid oceans via their optical properties. Given that discovering an ocean on an exoplanet would confirm its status as a habitable world, this science case is literally the raison d'etre of the Habitable Worlds Observatory.

Modern cosmological inference increasingly relies on differentiable models to enable efficient, gradient-based parameter estimation and uncertainty quantification. Here, we present a novel approach for predicting the abundance of dark matter haloes and their cosmology dependence using a differentiable, field-level neural network (NN) model, and study how well the cosmology dependence is captured by common parametrisations of the halo mass function (HMF), and by our NN-based approach. By training a 3D U-Net on initial density fields from fast N-body simulations with varying cosmological parameters, we enable direct, differentiable mapping from the linear density field to protohalo patches and their mass bins. Our method achieves competitive accuracy in identifying protohalo regions and in capturing the dependence of the HMF on cosmological parameters. Our NN derivatives agree well with finite differences of both analytical and emulated HMFs, at the level of the disagreement among the different models. We further demonstrate how the NN model can additionally be used to investigate the response of the HMF to changes in the initial Gaussian random field. Finally, we also demonstrate that a differentiable model can be used to extrapolate existing models at very high precision.

Recent studies have reported a non evolution of galaxy ages at redshifts higher than $z\sim$ 2.5, as well as galaxies older than the Universe. In this work, a sample of galaxies from JWST and HST was analysed via photometry to further understand this astronomical phenomenon. No prior cosmological parameters were assumed in the analysis, but the spectroscopic redshift. When compared to stellar population synthesis models, the results for mass-weighted galaxy ages indicate that the analysed objects seem to be divided into two subsets. The results for the subset with the majority of objects (60\% assuming a flat-$\Lambda$CDM cosmology) indicate an evolution of galaxy ages within the redshift range $z$=0.1-7.0, in the sense that higher redshift galaxies are younger than the Universe. Sources of systematic errors were discussed drawing into conclusion that degeneracies between reddening-age-metallicity, and/or AGN emission may explain the rest 40\% of the galaxies with ages older than expected from a flat-$\Lambda$CDM cosmology.

GRB 221009A is the brightest gamma-ray burst (GRB) observed to date. Extensive observations of its afterglow emission across the electromagnetic spectrum were performed, providing the first strong evidence of a jet with a nontrivial angular structure in a long GRB. We carried out an extensive observation campaign in very-high-energy (VHE) gamma rays with the first Large-Sized Telescope (LST-1) of the future Cherenkov Telescope Array Observatory (CTAO), starting on 2022 October 10, about one day after the burst. A dedicated analysis of the GRB 221009A data is performed to account for the different moonlight conditions under which data were recorded. We find an excess of gamma-like events with a statistical significance of 4.1$\sigma$ during the observations taken 1.33 days after the burst, followed by background-compatible results for the later days. The results are compared with various models of afterglows from structured jets that are consistent with the published multiwavelength data, but entail significant quantitative and qualitative differences in the VHE emission after one day. We disfavor models that imply VHE flux at one day considerably above $10^{-11}$ erg cm$^{-2}$ s$^{-1}$. Our late-time VHE observations can help disentangle the degeneracy among the models and provide valuable new insight into the structure of GRB jets.

Protoclusters of galaxies are overdense regions of the Universe characterized by large gas reservoirs. Such environments make them perfect laboratories to investigate galaxy-AGN co-evolution and the growth of SMBHs. Galaxies living in such a dense regions are expected to growth efficiently their SMBH, resulting in a higher incidence of AGN than in the field. Some protoclusters exhibit extended Ly$\alpha$ nebulae in their central region, pinpointing the presence of massive gas reservoirs, but whose main powering mechanism is still debated. We aim to investigate the AGN population, and AGN enhancement, in three protoclusters at 2.3 < z < 3.2 which host enormous Ly$\alpha$ nebulae (ELANe). Additionally, we search for the presence of X-ray diffuse emission in the same region of the Ly$\alpha$ nebulae to reveal multi-phase gas in these protoclusters. We use deep (190-270 ks) Chandra observations to identify AGN among the protocluster members and perform X-ray spectral analysis to derive the properties of those sources. We compare the AGN fraction and space density with those observed in other known protoclusters and from the field environment. We find 11 X-ray detected AGN in the three protoclusters. Each structure hosts a central, X-ray powerful (log$(L_{\rm X}/{\rm erg \, s^{-1}}) \sim 45-46$), QSO, while the other X-ray sources are mostly moderately luminous (log$(L_{\rm X}/{\rm erg \, s^{-1}}) \sim 44$) and obscured Compton-Thin AGN. The fraction of AGN in our targets is comparable with estimates for other protoclusters, and significantly higher than what is found for low-redshift clusters. We also find a significant enhancement (2-4 dex) of AGN density with respect to the field and to non active galaxies in the protoclusters. Finally, we find no significant soft X-ray diffuse emissions from the nebulae, thus ruling out gravitational heating as the main powering mechanism of the ELANe.

We discuss ASASSN-24fw, a 13th-magnitude star that optically faded by $\Delta g = 4.12 \pm 0.02$ mag starting in September 2024 after over a decade of quiescence in ASAS-SN. The dimmimg lasted $\sim$8 months before returning to quiescence in late May 2025. The spectral energy distribution (SED) before the event is that of a pre-main sequence or a modestly evolved F star with some warm dust emission. The shape of the optical SED during the dim phase is unchanged and the optical and near-infrared spectra are those of an F star. The SED and the dilution of some of the F star infrared absorption features near minimum suggest the presence of a $\sim$0.25$M_\odot$ M dwarf binary companion. The 43.8 year period proposed by Nair & Denisenko (2024) appears correct and is probably half the precession period of a circumbinary disk. The optical eclipse is nearly achromatic, although slightly deeper in bluer filters, $\Delta (g-z)=0.31\pm0.15$ mag, and the $V$ band emission is polarized by up to 4%. The materials most able to produce such small optical color changes and a high polarization are big ($\sim$20 $\mu$m) carbonaceous or water ice grains. Particle distributions dominated by big grains are seen in protoplanetary disks, Saturn-like ring systems and evolved debris disks. We also carry out a survey of occultation events, finding 42 additional systems, of which only 7 (4) closely match $\varepsilon$ Aurigae (KH 15D), the two archetypes of stars with long and deep eclipses. The full sample is widely distributed in an optical color-magnitude diagram, but roughly half show a mid-IR excess. It is likely many of the others have cooler dust since it seems essential to produce the events.

Seismic structure inversions have been used to study the solar interior for decades. With the high-precision frequencies obtained using data from the Kepler mission, it has now become possible to study other solar-like oscillators using structure inversions, including both main-sequence and subgiant stars. Subgiant stars are particularly interesting because they exhibit modes of mixed acoustic-buoyancy nature, which provide the opportunity to probe the deeper region of stellar cores. This work examines whether the structure inversion techniques developed for the pure acoustic modes of the Sun and other main-sequence stars are still valid for mixed modes observed in subgiant stars. We construct two grids of models: one of main-sequence stars and one of early subgiant stars. Using these grids, we examine two different parts of the inversion procedure. First, we examine what we call the "kernel errors", which measure how well the mode sensitivity functions can recover known frequency differences between two models. Second, we test how these kernel errors affect the ability of an inversion to infer known structure differences. On the main sequence, we find that reliable structure inversion results can be obtained across the entire range of masses and large frequency separations we consider. On the subgiant branch, however, the rapid evolution of mixed modes leads to large kernel errors and hence difficulty recovering known structure differences. Our tests show that using mixed modes to infer the structure of subgiant stars reliably will require improvements to current fitting approaches and modifications to the structure inversion techniques.

We present new methods to quantify the AGN population in terms of a multi-dimensional luminosity function that describes the space density of sources as a function of both X-ray and radio luminosity. We compile a sample of 1891 radio and X-ray detected extragalactic sources from the Boötes and COSMOS fields. First, we investigate the X-ray--radio luminosity correlation in the sample and find that an apparent correlation is introduced due to the sensitivity limits of the surveys; when considering individual redshift bins we find a wide range of radio luminosities associated with a given X-ray luminosity, and vice versa, indicating little direct connection between the emission processes. We then measure the X-ray luminosity function, radio luminosity function and multi-dimensional X-ray--radio luminosity function across redshift ($0<z<6$). We apply luminosity thresholds in X-ray and radio to restrict our sample to those in the AGN-dominated regime and explore how the fraction of radio-selected AGN within the overall X-ray sample varies with increasing X-ray luminosity (and vice versa). We find that towards the highest X-ray and radio luminosities the fraction of sources with both an X-ray and radio detection increases towards 100%, indicating that at the highest luminosities we are more likely to obtain a detection in both bands, though the source will not necessarily be bright in both bands. Thus, the most luminous accretion events are more likely to be associated with the production of a jet, despite the distinct physical structures that produce the emission and likely persist over very different timescales.

Simulation-based inference (SBI) has become an important tool in cosmology for extracting additional information from observational data using simulations. However, all cosmological simulations are approximations of the actual universe, and SBI methods can be sensitive to model misspecification - particularly when the observational data lie outside the support of the training distribution. We present a method to improve the robustness of cosmological analyses under such conditions. Our approach first identifies and discards components of the summary statistics that exhibit inconsistency across related simulators, then learns a transformation that brings the observation back within the support of the training distribution. We apply our method in the context of a recent SimBIG SBI galaxy clustering analysis using the wavelet scattering transform (WST) summary statistic. The original analysis struggled to produce robust constraints for certain subsets of WST coefficients, where the observational data appeared out-of-distribution (OOD) relative to the training data. We show that our method enables robust cosmological inference and resolves OOD issues, while preserving most of the constraining power. In particular, the improved SimBIG WST analysis yields $\Lambda$CDM constraints of $\Omega_m = 0.32^{+0.02}_{-0.02}$ and $\sigma_8 = 0.80^{+0.02}_{-0.02}$, which are respectively $1.4\times$ and $3.1\times$ tighter than those from a standard perturbation-theory-based power spectrum analysis, confirming the significant information gain of WST summary statistics. The proposed method is easily applicable to other cosmological SBI contexts and represents a step toward more robust SBI pipelines.

Machine learning (ML) and deep learning (DL) techniques are increasingly used across astrophysics, enabled by the growing availability of data and improved acquisition methods. These approaches now support tasks from redshift estimation to source classification. In this work, we aim to (i) classify blazars from the Fermi 4LAC-DR3 catalogue, in particular to identify the likely origin of blazars of uncertain type (BCUs), and (ii) investigate the full blazar sample to study their structure and redshift-luminosity evolution. We focus especially on the transition region between Flat Spectrum Radio Quasars (FSRQs) and BL Lacertae objects (BL Lacs), which may yield insights into accretion disk evolution. We examine Changing-Look Blazars (CLBs) as potential intermediates in this transition. We implement a classification pipeline using both a strong benchmark model (XGBoost) and a foundation model pre-trained on millions of tabular datasets (TabPFN). By extracting and reducing the high-dimensional latent space of the best model, we provide a 2D representation of the blazar population. This reveals a continuum between FSRQs and BL Lacs, including CLBs as transitional sources. These results support a scenario of gradual evolution from radiatively efficient (FSRQ-like) to inefficient (BL Lac-like) accretion. Ultimately, we show that a single probability score, combined with the latent space, offers a new framework for interpreting blazar diversity beyond discrete classes.

String theory naturally leads to the expectation that dark energy is not stable, and may be evolving as captured by the Swampland de Sitter conjectures. Moreover, motivated by the distance conjecture a unification of dark sector has been proposed, where the smallness of dark energy leads to one extra dimension of micron size with dark matter being the KK graviton excitations in this extra dimension. We consider the natural possibility that the radius of the dark dimension varies as the dark energy decreases, leading to the variation of the dark matter mass. This correlates the decrease of the dark energy with the variation of the dark matter mass as they depend on the variations of a scalar field $\phi$ controlling the radius of the extra dimension. A simple realization of this idea for small range of $\phi$ is captured by choosing a potential which is locally of the form $V=V_0\ {\rm exp}(-c\phi)$ and dark matter mass $m_{\rm DM}=m_0\ {\rm exp}(-c' \phi)$ where the sign of $\phi$ is chosen such that $c'\geq 0$ while we have two choices for the sign of $c$ depending on whether the dark dimension expands or shrinks when the dark energy dominates. We find excellent agreement with recent experimental data from DESI DR2 combined with SN measurements and reproduces the same significance as CPL parametrization with the added benefit of providing a natural explanation for the apparent phantom behavior ($w<-1$) reported by DESI and DES based on a physical model. Regardless of the SN dataset, there is a preference for non-zero values of $c'$ and $c$ that are in the expected $O(1)$ range in Planck units as suggested by the Swampland criteria. In particular, there is a remarkable consistency with $c'\simeq 0.05 \pm 0.01$ for all dataset combinations including SN, and close to the experimental upper bound of $c'\lesssim 0.2$ demanded by the lack of detection of fifth force in the dark sector.

Large-scale photometric surveys are revolutionizing astronomy by delivering unprecedented amounts of data. The rich data sets from missions such as the NASA Kepler and TESS satellites, and the upcoming ESA PLATO mission, are a treasure trove for stellar variability, asteroseismology and exoplanet studies. In order to unlock the full scientific potential of these massive data sets, automated data-driven methods are needed. In this review, I illustrate how machine learning is bringing asteroseismology toward an era of automated scientific discovery, covering the full cycle from data cleaning to variability classification and parameter inference, while highlighting the recent advances in representation learning, multimodal datasets and foundation models. This invited review offers a guide to the challenges and opportunities machine learning brings for stellar variability research and how it could help unlock new frontiers in time-domain astronomy.

We present stellar atmosphere modeling of JWST NIRCam photometry of nine highly magnified individual stars in a single galaxy at redshift z=0.94 known as the Warhol arc, which is strongly lensed by the galaxy cluster MACSJ0416. Seven of these transients were identified by Yan et al. (2023). The nine sources are all likely red supergiants with temperatures of T~4000K. We present new longslit spectroscopy of the Warhol arc acquired with Keck-I and the Large Binocular Telescope, and use these data to constrain the arc's oxygen abundance to be 12+log(O/H)=8.45+-0.08. We perform a microlensing simulation on synthetic stellar populations using a range of stellar metallicities and initial mass function slopes. The temperature distribution of the simulated detectable stars is sensitive to the choice of stellar metallicity, and setting the stellar metallicity equal to the arc's nebular metallicity (log(Z*/Zsun)=-0.24) produces a simulated temperature distribution that is consistent with the observations, while lower stellar metallicities (log(Z*/Zsun)<-0.75) produce simulated temperatures that are inconsistent with the observations. The expected detection rate is strongly anticorrelated with the IMF slope for {\alpha}>1.2. For the canonical IMF slope alpha=2.35, the simulation yields expected transient detection rates that agree with the observed detection rates in the HST Flashlights filters, but over predicts the detection rate by a factor of ~3-12 (<2sigma tension) in the JWST filters. The simulated detection rate is sensitive to the choice of stellar metallicity, with lower metallicities (log(Z*/Zsun)<-0.75) yielding a significantly lower simulated detection rate that further reduces the modest tension with the observations.

Massive stars at cosmological distances can be individually detected during transient microlensing events, when gravitational lensing magnifications may exceed mu ~ 1000. Nine such sources were identified in JWST NIRCam imaging of a single galaxy at redshift z = 0.94 known as the "Warhol arc,'' which is mirror-imaged by the galaxy cluster MACS J0416.1-2403. Here we present the discovery of two coincident and well-characterized microlensing events at the same location followed by a third event observed in a single filter approximately 18 months later. The events can be explained by microlensing of a binary star system consisting of a red supergiant (T ~ 4000 K) and a B-type (T > 13,000 K) companion star. The timescale of the coincident microlensing events constrains the projected source-plane size to R < 270 AU. The most likely binary configurations consistent with the observational constraints on the temperatures and luminosities of each star are stars with initial masses M1 = 22.5+7.1-5.5 Msun and an initial mass ratio very close to unity. A kinematic model that reproduces the observed light curve in all filters gives a relatively small transverse velocity of ~50 km/s. This requires the dominant velocity component of several hundreds of km/s to be roughly parallel to the microcaustic. An alternative possibility would be that the three microlensing events correspond to unrelated stars crossing distinct microcaustics, but this would imply a highly elevated rate of events at their common position, even though no underlying knot is present at the location.

The Star-Planet Activity Research CubeSat (SPARCS) is a NASA-funded 6U-CubeSat mission designed to monitor ultraviolet (UV) radiation from low-mass stars. These stars' relatively high-frequency and high-energy UV flares significantly affect the atmospheres of orbiting exoplanets, driving atmospheric loss and altering the conditions for habitability. SPARCS aims to capture time-resolved photometric data in the far-UV and near-UV simultaneously to better characterize the flares and detect the strongest and rarest among them. In addition, SPARCS is testing innovative technology, such as delta-doped detectors with near 100% internal quantum efficiency and detector-integrated metaldielectric UV bandpass filters. This mission will increase the technology readiness level of these critical components, positioning them for inclusion in future flagship missions like the Habitable Worlds Observatory. This paper outlines SPARCS' mission goals and provides an update as the spacecraft is completed and awaits its planned late-2025 launch to a sun-synchronous low-Earth orbit. It also highlights the critical role of small missions in providing training and leadership development opportunities for students and researchers, advancing technology for larger observatories, and shares lessons learned from collaborations between academic, government, and industry partners.

Giant radio sources, including galaxies and quasars (hereafter GRGs), are active galactic nuclei (AGN) hosting relativistic jets with source sizes exceeding the projected length of 0.7 Mpc. They are crucial to understanding the evolution of radio sources and their interaction with the surrounding environment. Some of these enigmatic objects, e.g., NGC 315, have also been reported as gamma-ray emitters. Since GRGs are thought to be aligned close to the plane of the sky, they are invaluable targets to explore the radiative mechanisms responsible for the observed gamma-ray emission. We have carried out a systematic search of gamma-ray emitting GRGs using sensitive low-resolution radio surveys, such as by Low Frequency Array, NRAO VLA Sky Survey, and Rapid ASKAP Continuum Survey, and considering the fourth data release of the fourth Fermi-Large Area Telescope gamma-ray source (4FGL-DR4) catalog. By carefully inspecting the radio maps of all AGN included in the 4FGL-DR4 catalog, we have identified 16 gamma-ray emitting GRGs, including 8 of them being reported as GRGs for the first time. Some of their observed parameters, e.g., core dominance, appeared to differ from that found for the non-gamma-ray detected GRG population, possibly due to the relatively small viewing angle of the gamma-ray emitting jet. The observed gamma-ray properties of these objects were found to be similar to non-GRG gamma-ray emitting misaligned AGN. We conclude that the origin of the gamma-ray emission could be similar in both source populations.

We derive one dimensional (1D) analytical solutions for the transport equations of incompressible magnetohydrodynamic (MHD) turbulence developed by Zank et al. [2012], Adhikari et al. [2023], including the Elsässer energies and the correlation lengths. The solutions are suitable for an arbitrary given background convection speed and Alfvén speed profiles but require near equipartition of turbulent kinetic energy and magnetic field energy. These analytical solutions provide a simple tool to investigate the evolution of turbulence and resulting energetic particle diffusion coefficients in various space and astrophysical environments that possess simple geometry.

Follow-up studies of persistent emission from Fast Radio Burst (FRB) sources are critical for understanding their elusive emission mechanisms and the nature of their progenitors. This work presents new observations of the persistent radio source (PRS) associated with \ourfrb. We observe a gradual decay in the PRS brightness, which is punctuated by periods of brightening and dimming at both 1.5 GHz and 3 GHz. Furthermore, our low-frequency ($<1$ GHz) observations-the first for this source-reveal evidence of a spectral break which can be attributed to absorption processes. Interpreted within the framework of the magnetar wind nebula model, our data constrain the age of the magnetar progenitor to $53^{+17}_{-10}$ years, broadly consistent with previous work. Assuming the observed 1.5 GHz variability is driven by scintillation, we derive a conservative lower limit for the source's radius of $>0.52$ pc, for a screen at 10 kpc. The observations presented here challenge the predictions of the previously published best-fit hypernebula model for this source.

We present a sample of 53 galaxy spectra at z_spec ~ 5-12 from the JWST CEERS and RUBIES surveys, combining NIRSpec PRISM spectroscopy with NIRCam photometry. We aim to use these data to establish best practices for measuring the UV spectral slope ($\beta$) in the era of JWST. We adopt power-law fits to the rest-frame UV continuum from the spectroscopic data as our fiducial, or `true', $\beta$ values, and compare them to photometric estimates derived through three methods: (1) photometric power-law fitting, (2) power-law fitting to an SED model fitted to the photometry, (3) single-color fitting near the Lyman break, and (4) single-color fitting at fixed rest-frame wavelengths. We find that photometric power-law fitting most closely recovers the spectroscopic slopes, with minimal bias and scatter. SED fitting performs moderately well, and can be preferable in cases of low signal-to-noise where photometric power-law fitting may become unreliable. Single-color estimates, while commonly used in past studies, show the most significant deviations and are not recommended when more than a single color is available. Our results highlight the limitations and strengths of each approach and provide practical guidance for measuring $\beta$ from photometry when spectra are unavailable or are of insufficient quality.

We present a sample of $z_{\rm phot}\sim10-16$ galaxies by exploiting one of the richest JWST NIRCam imaging data, taken in the CANUCS survey in Cycle 1 and the Technicolor (TEC) survey in Cycle 2. The combination of the CANUCS+TEC provides multi-epoch, deep NIRCam images in all medium bands (MBs) and broad bands (BBs) onboard NIRCam (22 filters in total), over $\sim23\ {\rm arcmin}^2$ in three independent lines of sight. We select high-$z$ galaxy candidates based on photometric redshifts, and obtain eight candidates at $z\sim10-16$, including a very robust candidate at $z\sim15.4$. The ultraviolet (UV) luminosity function (LF) from our sample is consistent with previous JWST studies showing a scatter of $\sim0.6$ dex across the literature, marking the significance of the field-to-field variance in interpreting galaxy abundance measurements at $z>10$. We find that the UV LF moderately evolves at $z>10$, and the LF normalization and the luminosity density decline by a factor of $\sim7$ from $z\sim11$ to $z\sim15$, indicating less steep evolution than $z<11$. We highlight the importance of MB filters, not only to minimize the contamination by low-$z$ interlopers but also to maximize the completeness. In particular, faint and less blue galaxies could be missed when the sample is built solely on BB data. The contamination and incompleteness of BB-only selected samples can bias our views of earliest galaxy evolution at $z>10$, including the UV LF by $\sim0.6$ dex, the size-magnitude relation by $\sim0.6$ dex, and the UV slope-magnitude relation by $\Delta\beta_{\rm UV}\sim-0.3$.

Jellyfish galaxies provide direct evidence of ram pressure stripping in cluster environments. We investigate the role of magnetic fields in the formation of jellyfish galaxies with a multiphase interstellar medium (ISM) using radiation magneto-hydrodynamic simulations. We impose magnetized (MHD) and non-magnetized (HD) winds on the gas-rich dwarf galaxies containing the magnetized or non-magnetized ISM. The MHD winds strip the disk gas more effectively than the HD winds because of the magnetic force acting against the local density gradient, which results in remarkably different ram pressure stripped features. The magnetic fields induced by the MHD winds generate a strong magnetic pressure, which forms smoothed disks and tail gas features. Since the stripped ISM in MHD wind cases travels while being nearly isolated from the intracluster medium (ICM), the stripped ISM mostly forms stars within 20~kpc of the galactic disks. In contrast, non-magnetized winds facilitate the efficient mixing of the stripped ISM with the ICM, resulting in the formation of abundant warm clouds that cool and collapse in the distant ($\sim50-100\,$kpc) tails at times of a few hundred Myr. Consequently, distant tail star formation occurs only in the HD wind runs. Finally, despite the different tail features, the star formation rates in the disk remain similar owing to the interplay between the increased gas stripping and the gas density increase in the disks of the MHD wind runs. These results suggest that the magnetized ICM may have a significant influence on jellyfish galaxies, whereas the magnetized ISM play a minor role.

Investigating how small-scale physical processes shape large-scale astrophysical phenomena is one of the key science themes of the Astro2020 decadal Survey. An example of this interplay is Cosmic Reionization, where ionizing photons from galaxies escaped into the intergalactic medium (IGM), driving its transition from neutral to ionized at $z>6$. Star-forming galaxies are thought to be the dominant sources of reionization. However, an open question remains as to whether bright or faint star-forming galaxies are the primary reionization contributors. This depends on the escape fraction $f_{~esc}^{LyC}$ -- the fraction of ionizing photons that successfully escape the galaxy's interstellar and circumgalactic medium to reach the IGM. Performing direct measurements of $f_{~esc}^{LyC}$ during the reionization epoch is not feasible, due to the near-zero transmission of the IGM at $z\gtrsim4$. However, by calibrating $f_{~esc}^{LyC}$ against indirect indicators (i.e. galaxy properties that correlate with $f_{~esc}^{LyC}$) at lower $z$, we can estimate its value in galaxies during reionization. Currently, $f_{~esc}^{LyC}$ calibrations are mostly limited to $\sim 90$ sources brighter $M_{UV} \sim -18$ at $z\sim0.3$. Yet, recent models suggest that fainter galaxies ($-19<M_{UV}<-13$) may be the major contributors to reionization. Our science goal is to extend $f_{~esc}^{LyC}$ calibrations to such faint magnitudes. With the Habitable Worlds Observatory, we aim to achieve a $\sim$ 20-fold increase in statistical power. Additionally, we plan to reach magnitudes as faint as 1/100 $L^*$ at $z\sim0.1$. Suitable targets will be selected from the upcoming Ultraviolet Explorer (UVEX) survey. Through this science case, we will be able to place robust constraints on the role of faint galaxies in the reionization of the Universe and identify its primary contributors.

We present the science case for characterizing the origin of the heaviest elements on the periodic table, with a focus on those produced by the rapid neutron-capture process (r-process), using the Habitable Worlds Observatory (HWO). High-resolution ultraviolet (UV) spectroscopy can increase the number of r-process elements detectable in cool stars by more than 50% relative to optical and infrared spectra. These elements are key to characterizing the physical conditions that govern the r-process and identify the nature, sites, and environments where r-process events occurred. HWO has the potential to greatly expand the sample of stars where rarely studied heavy elements can be detected beyond the Solar neighborhood to the Galactic halo, globular clusters, and dwarf galaxies.

We present the science case for characterizing the nature of the first stars using the Habitable Worlds Observatory (HWO). High-resolution ultraviolet (UV) spectroscopy with the HWO has the potential to confirm any surviving low-mass zero-metallicity first stars by placing unprecedented low limits on their metal abundances. It also has the potential to substantially increase the number of elements detectable in the spectra of known long-lived low-mass stars, which exhibit extremely low metal abundances that reveal the metals produced by the first stars. Elements important for this science case with UV transitions include C, Mg, Al, Si, P, S, Sc, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn. HWO would expand the discovery space when compared with the Hubble Space Telescope by enabling high-resolution UV spectroscopy for much fainter stars throughout the Milky Way and neighboring stellar systems.

We use intergalactic medium (IGM) metallicity distributions from several state-of-the-art cosmological simulations of Milky Way analogs and a semi-analytic model of ultra-faint dwarf galaxy (UFD) formation to model the stellar metallicities of UFDs in MW-like environments. We study simulations with different treatments of star formation, stellar feedback, and Population III enrichment, and in all cases, we find that only a few percent of the IGM accretable by UFD progenitors is enriched to metallicities $\rm [Fe/H]\ge-4$. When the metallicity of accreted IGM in the semi-analytic galaxy formation model is set using these IGM metallicity distributions, the model underpredicts UFD metallicities and their scatter compared to the observed luminosity--metallicity relation. Our results indicate that IGM enrichment is not the dominant mechanism setting metallicities of UFD stars. Instead, UFD stellar metallicity is determined primarily by the interplay between internal enrichment and metal loss through feedback-driven outflows. We examine models with different values of the maximum outflow mass loading factor $\eta_{\rm max}$ and show that the full range of average stellar metallicities of UFDs at $M_V<-7$ can be reproduced if the maximum mass loading factor varies in the range $200\lesssim\eta_{\rm max}\lesssim 2000$. We also consider stellar metallicity distribution functions (MDFs) within individual model galaxies with different assumptions about IGM enrichment and $\eta_{\rm max}$. We find that all considered models are in reasonable agreement with observed UFD MDFs, with model differences less than the uncertainties of current metallicity measurements.

Multiphase filamentary nebulae are ubiquitous in the brightest cluster galaxies (BCGs) of cool-core clusters, providing key insights into the cycle of baryons and the feeding and feedback of supermassive black holes. However, BCGs account for less than 1% of all early-type galaxies (ETGs). To broaden our understanding of how multiphase filamentary nebulae form in ETGs and connect to the greater picture of galaxy evolution, it is crucial to explore ETGs that are outside of the dense centers of galaxy clusters or groups. We present VLT-MUSE IFU observations of 126 nearby non-central ETGs, detecting warm ionized gas in 54 of them. 35 out of 54 display rotating disks that are morphologically and kinematically aligned with their stellar components, suggesting stellar mass loss as the origin of their warm gas. The remaining 19 host filamentary nebulae that are decoupled from the stellar components, resembling those observed in BCGs. These filamentary sources display unique emission line properties that cannot be fully explained by photoionization from post-asymptotic giant branch stars, active galactic nuclei, or fast gas shocks. For the eight filamentary sources that have been observed with Chandra, their soft X-ray emission indicates the presence of hot gas. We posit that their emission lines may be powered by EUV and X-ray radiation from the cooling of the hot gas, similar to cool-core clusters, but the detailed mechanisms and physical conditions may be different. As a detailed study case, we investigate NGC 4374, a non-central ETG with extensive Chandra observations, and find that its warm filaments are over pressured compared to the hot filaments - opposite with what is observed in cool-core clusters.

In this paper, we summarize the shock physics and the treatment of radiative transfer in two well-established shock codes -- the MAPPINGS code (Dopita1976,Binette1985,Sutherland1993) and the Cox/Raymond code (hereafter CR code) (Cox1972,Raymond1976,Raymond1979). We compare the ionization states, temperatures, electron densities, and the energy transportation of the shock models with shock velocities of 50, 110, 150 and 300 km/s. In summary, both codes adopt the Rankine-Hugoniot flow equation to describe the shock flows, giving the same shock physical properties at the immediate area behind shock fronts. The different treatments of radiative transfer in these two codes leads to somewhat different computation of the ionization and thermal structures of shocks, as well as the emission-line fluxes. This work highlights the importance of the delicate treatment of photoionization in shock models, providing insight of the future development of shock codes, such as the 3D shock codes.

We analyze imaging from Data Preview 1 of the Vera C. Rubin Observatory to explore the performance of early LSST pipelines in the 47 Tucanae field. The coadd-\texttt{object} catalog demonstrates the depth and precision possible with Rubin, recovering well-defined color magnitude diagrams for 47 Tuc Small Magellanic Cloud. Unfortunately, the existing pipelines fail to recover sources within $\sim$28 pc of the cluster center, due to the extreme source density. Using Rubin's forced photometry on stars identified via Difference Imaging, we can recover sources down to $\sim$14 pc from the cluster center, and find 14744 potential cluster members with this extended dataset. While this forced photometry has significant systematics, our analysis showcases the potential for detailed structural studies of crowded fields with the Rubin Observatory.

We investigate the possibility that the recently identified population of high-redshift, obscured quasars - known as "Little Red Dots" (LRDs) - originates from early black hole seed formation driven by ultra-strongly self-interacting dark matter (uSIDM). In this framework, dark matter halos undergo gravothermal core collapse due to large self-interaction cross sections, resulting in the rapid formation of massive black hole (BH) seeds with masses $\gtrsim 10^{5} M_\odot$ at redshifts $z \gtrsim 5$. We develop a semi-analytic model that tracks the evolution of the dark matter halo population, the redshift of collapse $z_{\rm coll}$, and the corresponding BH mass function. Black hole growth is modeled stochastically via a log-normal Eddington ratio distribution and a finite duty cycle. We find that the uSIDM scenario naturally reproduces key observed properties of LRDs, including their abundance, compactness, and characteristic BH masses, while offering a mechanism for early, obscured black hole formation that is difficult to achieve in standard CDM-based models. The predicted SMBH mass function at $z \sim 5$ shows excellent agreement with LRD observational data and SIDM merger-tree simulations, particularly at the high-mass end $(m_{\rm BH} \gtrsim 10^{7} M_\odot)$. These results suggest that LRDs may serve as powerful observational tracers of exotic dark sector physics and that SMBH formation in the early universe could be significantly shaped by non-gravitational dark matter interactions.

Over the past decade, an abundance of information from neutron-star observations, nuclear experiments and theory has transformed our efforts to elucidate the properties of dense matter. However, at high densities relevant to the cores of neutron stars, substantial uncertainty about the dense matter equation of state (EoS) remains. In this work, we present a semiparametric EoS framework aimed at better integrating knowledge across these domains in astrophysical inference. We use a Meta-model at low densities, and Gaussian Process extensions at high densities. Comparisons between our semiparametric framework to fully nonparametric EoS representations show that imposing nuclear theoretical and experimental constraints through the Meta-model up to nuclear saturation density results in constraints on the pressure up to twice nuclear saturation density. We show that our Gaussian Process trained on EoS models with nucleonic, hyperonic, and quark compositions extends the range of EoS explored at high density compared to a piecewise polytropic extension schema, under the requirements of causality of matter and of supporting the existence of heavy pulsars. We find that maximum TOV masses above $3.2 M_{\odot}$ can be supported by causal EoS compatible with nuclear constraints at low densities. We then combine information from existing observations of heavy pulsar masses, gravitational waves from binary neutron star mergers, and X-ray pulse profile modeling of millisecond pulsars within a Bayesian inference scheme using our semiparametric EoS prior. With current astrophysical observations, we find a favored pressure at two times nuclear saturation density of $P(2\rho_{\rm nuc}) = 1.98^{+2.13}_{-1.08}\times10^{34}$ dyn/cm$^{2}$, a radius of a $1.4 M_{\odot}$ neutron star value of $R_{1.4} = 11.4^{+0.98}_{-0.60}$\;km, and $M_{\rm max} = 2.31_{-0.23}^{+0.35} M_{\odot}$ at the 90\% credible level.

The symbiotic channel of Type Ia supernovae progenitors is crucial for explaining the observed circumstellar material in some Type Ia supernovae. While extensive numerical and observational efforts have been dedicated to exploring the progenitor system, limited emphasis has been placed on studying the surviving companions arising from the symbiotic channel. In this paper, we present a numerical study of the symbiotic systems using {\tt MESA} as potential Type Ia supernova progenitors. We conduct 1260 binary stellar evolution simulations, over a wide range of parameters, incorporating the optically thick wind model developed by Hachisu et al., and predict the post-impact evolution of these surviving companions. We classify four types of progenitor systems based on the evolutionary stage of the companion at the onset of the explosion: red giant companions, with or without prior helium flash events, and asymptotic giant branch companions, with or without the thermal pulsing phase. After the SN impact, a blue dwarf star with either a helium or carbon-oxygen core is left behind. However, if a small portion of the envelope ($\gtrsim$ 0.3\%) remains on the core of the surviving companion, the overall post-supernova evolution may remain similar to its pre-explosion state, albeit slightly fainter, making observation a challenging endeavor.

White-light flares are explosive phenomena accompanied by brightening of continuum from near-ultraviolet to optical, which occur on the Sun and stars. In order to investigate the mechanism of white-light flares, we carried out simultaneous optical photometry (TESS : 6000-10000 Å) and spectroscopy (Seimei Telescope : 4100-8900 Å) of a M-dwarf EV Lac on 2019 September 14. We detected a flare with high-time-cadence ($\sim 50$ sec) spectroscopic observation. At the peak, the continuum of the flare component is well fitted by a blackbody spectrum with temperature of $T = 8122 \pm 273$ K, which is comparable with the results of previous studies that reported the spectral energy distribution of near-ultraviolet to optical during the flare could be approximated by single-temperature blackbody radiation at $T \sim 10^{4}$ K. We also estimated the time evolution of the flare temperature during the decay phase. The radiative energy of this flare within the optical range is $4.4 \times 10^{32}$ erg, taking into account the time-dependent variation in the decreasing flare temperature and expanding flare area. Furthermore, we detected a delayed increase in the flux of H$\alpha$ after the photometric flare peak, secondary increase, and gradual increase even after the white-light flare ended. Comparison of our results with light curves obtained by the Sun-as-a-star analysis of solar flares indicates that these signals may be due to postflare loops near the stellar limb. Our result about time evolution of white-light continuum will help to gain more insight into the mechanism of white-light flares both on the Sun and stars. Additionally, since extreme ultraviolet radiation from flare loops plays a key role in planetary atmospheric escape, the existence of postflare loops on stellar flares and its time evolution will help future studies about habitability of close-in planets.

We present optical and near-infrared (NIR) observations of the outbursting, Halley-type comet 12P/Pons-Brooks. Three NIR spectra were obtained during two outbursts in October and November 2023, with the 3-meter Infrared Telescope Facility and the Palomar 200-inch Telescope, respectively. The NIR spectra exhibited absorption features at 1.5 and 2.0 $\mu$m, consistent with the diagnostic absorption bands of water ice, superimposed on a red dust-scattering continuum. We find that the absorption bands and the red continuum can be well explained by micrometer-sized crystalline ice at 140--170 K, along with sub-micrometer-sized refractory grains (e.g., amorphous carbon). In addition, an optical spectrum was obtained with the Lijiang 2.4-meter Telescope during the November 2023 outburst, which exhibited the emission bands of gaseous CN, C$_3$, C$_2$ and NH$_2$. The C$_3$/CN and C$_2$/CN ratios suggest that 12P/Pons-Brooks was ''typical'' in C$_3$ abundance but somewhat depleted in C$_2$. The specific kinetic energy of the 2023 November outburst is estimated to be $\sim8\times10^3$ J kg$^{-1}$, suggesting a likely triggering mechanism similar to 332P/Ikeya--Murakami and 17P/Holmes, i.e., the crystallization of amorphous water ice. A refractory-to-ice ratio of $\sim$1.7--3.2 is derived from the total mass loss of dust and gas, aligning with the lower-end estimates for 67P/Churyumov-Gerasimenko and 1P/Halley. This suggests either a less evolved nucleus or an outburst region enriched in icy materials relative to the bulk nucleus.

Accurate redshift measurements are essential for studying the evolution of quasi-stellar objects (QSOs) and their role in cosmic structure formation. While spectroscopic redshifts provide high precision, they are impractical for the vast number of sources detected in large-scale surveys. Photometric redshifts, derived from broadband fluxes, offer an efficient alternative, particularly when combined with machine learning techniques. In this work, we develop and evaluate a neural network model for predicting the redshifts of QSOs in the Dark Energy Spectroscopic Instrument (DESI) Early Data Release spectroscopic catalogue, using photometry from DESI, the Widefield Infrared Survey Explorer (WISE) and the Galactic Evolution Explorer (GALEX). We compare the performance of the neural network model against a k-Nearest Neighbours approach, these being the most accurate and least resource-intensive of the methods trialled herein, optimising model parameters and assessing accuracy with standard statistical metrics. Our results show that incorporating ultraviolet photometry from GALEX improves photometric redshift estimates, reducing scatter and catastrophic outliers compared to models trained only on near infrared and optical bands. The neural network achieves a correlation coefficient with spectroscopic redshift of $0.9187$ with normalised median absolute deviation of $0.197$, representing a significant improvement over other methods. Our work combines DESI, WISE and GALEX measurements, providing robust predictions which address the difficulties in predicting photometric redshift of QSOs over a large redshift range.

Cosmic rays, originating from stars, supernovae, and other astrophysical sources, are composed of high-energy particles that enter Earths atmosphere. Upon interaction with atmospheric nuclei, these primary cosmic rays generate secondary particles, including neutrons, electrons, and muons, with muons constituting a dominant component at ground level. Muons, due to their relative abundance, stability, and well-characterized energy loss mechanisms, serve as critical probes for investigating the fundamental properties of cosmic rays. Studies of muon energy distribution, diurnal anisotropy, and their modulation by solar activity provide critical insights into the mechanism of particle acceleration in cosmic ray sources and the effects of solar and this http URL study aims to characterize the counting spectra and anisotropic properties of cosmic ray muons by using a plastic scintillator detector system. The experiment was conducted over a three-month period, from December 2023 to February 2024, leveraging long-bar plastic scintillator detectors equipped with dual-end photomultiplier tubes (PMTs) and a high-resolution digital data acquisition system. A dual-end coincidence measurement technique was used to enhance the signal-to-noise ratio by suppressing thermal noise and other background interferences. Diurnal variations in muon count rates exhibit a pronounced pattern, with a systematic reduction occurring between 8:00 AM and 1:00 PM. This phenomenon is attributed to the solar shielding effects, where enhanced solar activity during daytime hours modulates the flux of galactic cosmic rays reaching Earths surface. The study further corroborates these findings through cross-comparisons with data from the Yangbajing Cosmic Ray Observatory. These observations underscore the robustness of the plastic scintillator detector system for capturing detailed muon spectra and anisotropic patterns.

Despite decades of research, the fundamental processes involved in the initiation and acceleration of solar eruptions remain not fully understood, making them long-standing and challenging problems in solar physics. Recent high-resolution observations by the Goode Solar Telescope have revealed small-scale magnetic flux emergence in localized regions of solar active areas prior to eruptions. Although much smaller in size than the entire active region, these emerging fluxes reached strengths of up to 2000 G. To investigate their impact, we performed data-constrained magnetohydrodynamic (MHD) simulations. We find that while the small-scale emerging flux does not significantly alter the pre-eruption evolution, it dramatically accelerates the eruption during the main phase by enhancing the growth of torus instability, which emerges in the nonlinear stage. This enhancement occurs independently of the decay index profile. Our analysis indicates that even subtle differences in the pre-eruption evolution can strongly influence the subsequent dynamics, suggesting that small-scale emerging flux can play a critical role in accelerating solar eruptions.

I explore models of the dust-scattered component of the Cosmic Ultraviolet Background (CUVB) at the North Galactic Pole (NGP) in order to develop a framework for calculating the dust-scattered light as a function of the optical depths. As expected, I find that the dust-scattered emission scales linearly with reddening up to $E(B-V) \approx 0.1$\ mag and derive a parametric model for this dependence. I have applied these models to fit the far-ultraviolet (1350--1800 Å) observations from the \textit{Galaxy Evolution Explorer (GALEX)} finding that the optical constants of the interstellar dust grains -- albedo ($a$) and phase function asymmetry factor ($g$) -- are consistent with predictions from the Astrodust model ($a = 0.33$, $g = 0.68$). I detect an isotropic offset of $267 \pm 7$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ Å$^{-1}$, half of which remains unaccounted for by known Galactic or extragalactic sources. I will now extend my analysis to wider sky regions with the goal of generating high-resolution extinction maps.

It is suggested that the variation of mass accretion rate in accretion disk may be responsible for the occurrence of most changing-look active galactic nuclei (CL AGNs). However, the viscous timescale of a thin disk is far longer than the observed timescale of CL AGNs. Though this problem can be resolved by introducing the large-scale magnetic field, the mechanism for radio-quiet CL AGNs with weak/absent large-scale magnetic field remains a mystery. In this work, we assume that the thin accretion disk is collapsed from the inner advection-dominated accretion flow (ADAF) instead of substituting by the outer thin disk through advection. This idea is tested by comparing the cooling timescale ($t_{\rm cool}$) of an ADAF with the observed timescale ($t_{\rm tran}$) of turn-on CL AGNs. We compile a sample of 102 turn-on CL AGNs from the archived data and calculate the cooling timescale of an ADAF with the critical mass accretion rate based on some conventional assumptions. It is found that $t_{\rm cool}$ is much shorter than $t_{\rm tran}$ in most of the CL AGNs, which validates our assumption though $t_{\rm cool}$ is not consistent with $t_{\rm tran}$ ($t_{\rm cool}<t_{\rm tran}$). However, this is reasonable since most of the CL AGNs were observed only two times, indicating that the observed timescale $t_{\rm tran}$ is the maximum value because the changing-look can indeed happen before the second observation.

Fermi-LAT observations revealed that each GeV phase-folded light-curve (aka. phaseogram) of the Crab, Geminga, Dragonfly and Vela pulsars consists of two pulses (P1 & P2) and a "Bridge" between them. There is clearly a "bump" at the Bridge phase of Vela's pulse profiles, that could also be regarded as the third pulse (P3). Differently, the Crab's, Geminga's & Dragonfly's bridges relatively resemble a "valley floor". Despite such an apparent difference, it is interesting to investigate whether their bridge emissions are still within the same general picture as Vela's. Assuming the north-south symmetry, we would expect the fourth component (Bridge2/P4) to exist as well. However, such a hypothetical Bridge2/P4 is not intuitively identified on $\gamma$-ray phaseograms of the Crab, Geminga, Dragonfly and Vela pulsars. It is also intriguing to hint at the rationale for the non-discovery of Bridge2/P4. Our prototypical toy model is free of assumptions on emission regions or radiation mechanisms. Instead, it assumes a north-south symmetric geometry and one circularly symmetric beam per hemisphere, while taking into account Doppler shifts (the most innovative element), time delays and energy-dependent beam shapes. Tentative compatibility of our fitting results with wind models is reported. Notably, for the Crab pulsar, we found a preliminary qualitative correlation between our model predictions and the IXPE results on X-ray polarisation. The softer $\gamma$-ray pulsation of the Geminga pulsar is found to span over its full phase. Prompted by systematic evaluations, we outline some potential improvements for our toy model.

The cold outer regions of protoplanetary disks are expected to contain a midplane-centered layer in which gas-phase CO molecules freeze-out and their overall abundance is low. The layer then manifests itself as a void in the channel maps of CO rotational emission lines. We explore whether the frozen-out layer can expose the circumplanetary environment of embedded accreting protoplanets to observations. To this end, we performed 3D radiative gas-dust hydrodynamic simulations, directly linking opacities and optical depths to the redistribution of sub-$\mu$m- and mm-sized dust grains. A Jupiter-mass planet with accretion luminosity $\sim$$10^{-3}\,L_{\odot}$ was considered as the nominal case. The accretion heating sustains a warm bubble around the planet, locally increasing the abundance of gas-phase CO molecules. Radiative transfer predictions of the emergent sky images show that the bubble becomes a conspicuous CO emission source in channel maps, appearing as a low-intensity optically thick spot located in-between the `dragonfly wings' that trace the foreside and backside line-forming surfaces. The emission intensity of the bubble is nearly independent of the tracing isotopologue, suggesting a very rich observable chemistry, as long as its signal can be deblended from the extended disk emission. This can be achieved with isotopologues that are optically thin or weakly thermally stratified across the planet-induced gap, such as C$^{18}$O. For those, the bubble stands out as the brightest residual in synthetic ALMA observations after subtraction of axially-averaged channel maps inferred from the disk kinematics, enabling new automatic detections of forming protoplanets. By contrast, the horseshoe flow steadily depletes large dust grains from the circumplanetary environment which becomes unobservable in the sub-mm continuum, in accordance with the scarcity of ALMA detections.

Very massive stars (VMS) are defined as stars with an initial mass in excess of 100 Msun. Because of their short lifetime and the shape of the stellar mass function, they are rare objects. Only about twenty of them are known in the Galaxy and the Large Magellanic Cloud. However VMS are important in several ways. They efficiently spread nucleosynthesis products through their boosted stellar winds, they are predicted to explode as pair-instability supernovae or to form heavy black-holes from direct collapse, and they outshine all other types of stars in the ultraviolet light, thus dominating the integrated light of starbursts. Their presence is indirectly suspected across all redshifts, all the way to cosmic dawn where they may have played a key role in the formation of the first galaxies. Their search and identification is currently hampered by instrumental limitation, especially spatial resolution. An integral field spectrograph working at the diffraction limit of HWO (5mas) and with a spectral resolution of about 2000 would revolutionize the understanding of VMS. We make the case for such an instrument in this contribution.

Spotless days (SLDs) as well as CMEs in the decay phase of the solar cycle are believed to be a good predictor of the forthcoming cycle. A sequential increase in SLDs is observed since cycle 21, and cycle 24 has the highest number of SLDs (since cycle 14), which offers a unique opportunity to probe the CME characteristics that occurred during SLDs (hereafter CME_SLD), in a statistical sense. Here, we investigate the CME_SLD during the descending phases of solar cycles 23 (2004-2008) and 24 (2015-2019). The fraction of CMEs that occurred on SLDs is found to be 14 and 11%, for cycles 23 and 24, respectively, compared to the total CMEs that occurred in the aforementioned durations. CME_SLD that occurred on the visible side of the solar disk are found to be slower, smaller in width, and carrying low Kinetic energy and mass compared to the entire population of CMEs. The distribution of annual evolution of speed, angular width and acceleration of the CMESLD with the CMEs that occurred on the non-SLD days (hereafter CME_Non-SLD) for the descending phases shows that the CME_SLD are different from the CME_non-SLD in terms of characteristics (such as speed, width and acceleration). A comparative analysis of CME_SLD kinematics of cycle 23 and 24 shows that the weakest cycle 24 has wider and more massive events. In contrast, other parameters of CME_SLD such as speed, acceleration and Kinetic energy do not have a disparate nature. CME_SLD in both SC23 and SC24 are statistically similar. Therefore, this investigation suggests that SLDs, and hence the sunspot number, may not be a sufficient candidate to predict the solar eruptive activities (e.g. CMEs). On the other hand, our analysis of the relation between the strength of the geomagnetic storm and probable candidate parameters revealed Dst index to be very well correlated with the product of V and Bz.

Cosmologically coupled black holes (CCBHs) are alternative black hole models whose masses evolve as $M \propto a^3$ on cosmological scales. This characteristic suggests that CCBHs could contribute to the accelerated expansion of the Universe. In this paper, I consider a CCBH model in which the cosmological constant is effectively induced, while the baryonic mass is conserved within conventional black holes. This model is motivated by the theoretical framework of Schwarzschild - de Sitter black holes. Assuming that the accelerated cosmic expansion is caused by CCBHs, I perform a cosmological parameter estimation using datasets including Planck 2018 CMB, CMB lensing, BAO, and supernovae. The analysis reveals notable shifts in cosmological parameters, such as $H_0 = 72.24^{+0.34}_{-0.35} \mathrm{km/s/Mpc}$, compared to the standard $\Lambda \mathrm{CDM}$. My $H_0$ constraint is consistent with the value $H_0 = 73.04 \pm 1.04 \mathrm{km/s/Mpc}$ reported by SH0ES within $1 \sigma$. However, the overall fit to the data worsens, with a total $\chi^2 = 2884.12$ for the CCBH model, compared to $\chi^2 = 2836.12$ for the $\Lambda$CDM model. I show that the effect of cosmological coupling is suppressed by a factor of $10^{-16}$ at $\sim$pc scales, rendering it negligible compared to the standard black hole mass in local astrophysical phenomena, although the CCBH model can explain the accelerated expansion.

As part of the SKYSURF Hubble Space Telescope (HST) Legacy Archival program we present galaxy number counts which yield measurements of the extragalactic background light (EBL) at 15 different wavelengths. We have processed 82,752 HST images across 23 filters into 16,686 mosaics using the same software and processing pipeline throughout. Using 17/23 filters that give reliable galaxy counts, we constrain the integrated galaxy light (IGL) with a 1.5-9\% error between 0.3 and 1.6 $\mu$m in combination with 8 bands from WAVES (Wide Area VISTA Extragalactic Survey) and DEVILS (Deep Extragalactic Visible Legacy Survey). While HST was never intended to undertake large area surveys, through extensive quality control and filtering, we were able to extract a reliable and representative sample of fields distributed across the sky. Our final catalogs cover a combined $\approx 19.6 °^2$, with individual filters covering areas ranging from $\approx 0.16-7.0 °^2$. The combination of numerous independent sight-lines and area coverage allows us to reduce cosmic variance uncertainties in deep number counts to 0.06\%-1.8\%. For the first time we are able to establish a measurement of the IGL, $\mathrm{9.07 \pm 0.35 nW m^{-2} sr^{-1}}$, at 0.59 $\mu$m using HST data. We obtain a cosmic optical background value of $ 24.45 \pm 0.50 \mathrm{nW m^{-2} sr^{-1}}$. Different techniques used to measure the COB, both directly and indirectly, have recently converged indicating that the COB arises almost exclusively from processes within galaxies. This in combination with the recent values reported from New Horizons and very high energy (VHE) constraints leaves very little room for any diffuse emission coming from outside the Milky Way.

We present new JWST/NIRCam observations of the starburst irregular galaxy NGC 4449, obtained in Cycle 1 as part of the Feedback in Emerging extrAgalactic Star clusTers (FEAST) program, which we use to investigate its resolved stellar populations and their spatial distributions. NGC4449 NIR color-magnitude diagrams reveal a broad range of stellar populations, spanning different evolutionary phases, from young main sequence stars, to old red giant branch stars and asymptotic giant branch (AGB) stars. The analysis of their spatial distributions shows that younger (< 10 Myr) populations form an S-shaped distribution aligned with the galaxy's north-south axis, while stars aged 10 - 60 Myr show shifting concentrations from the north to the south, consistent with the possibility that external interactions or tidal effects may have triggered star formation in spatially distinct bursts. Clusters of comparable ages generally follow these distributions, suggesting that cluster and field stars form at the same pace in each galaxy region. Thanks to the unprecedented high-spatial resolution and sensitivity of the JWST data we recover a clear gap between Oxygen-rich and the carbon star branch of the AGB population and the presence of a massive AGB star "finger". The analysis of these stars can provide constraints on AGB evolution models and dust production in this galaxy. These results confirms NGC 4449 status as a compelling example of a local dwarf starburst galaxy undergoing complex and possibly external driven star formation and underscore the power of JWST in probing the full lifecycle of stars in nearby starburst systems.

Based on a 10 years sample of gamma-ray flares of FSRQs collected with FERMI and AGILE, I report on this proceeding the advance on a statistical study of variability for a sample of more than 300 FSRQs. I will focus on waiting time between flares (defined as the time intervals between consecutive activity peaks (Pacciani 2022). The investigation revealed that gamma-ray activity can be modeled with overlapping bursts of flares, with flares uniformly distributed within each burst, and a typical burst rate of 0.6 y$^{-1}$. Moreover, a statistically relevant fast component with timescale of order of days is revealed. From these results, constraints on flares emission mechanisms were derived. I also discuss the preliminary results on an investigation of flares luminosity and duration in gamma-rays. A Simple fitting model is shown, correlating peak luminosity and duration of gamma-ray flares (paper in preparation).

UV auroral emissions from giant planets are produced by extra-atmospheric energetic particles interacting with an atmosphere. They have been observed on Jupiter, Saturn and Uranus and should be present on Neptune. Even if the mechanisms are similar, each planet is unique due to its specific source of magnetospheric plasma and the structure and dynamics of its magnetosphere. How these precipitations modify atmospheric heating, dynamics and chemical balance at local and global atmospheric scale is still poorly known, especially on Uranus and Neptune, and critical to understanding the global atmosphere-magnetosphere system of giant planets and exoplanets. In this manucript we present how future observations by instruments, aboard \textit{the Habitable World Observatory} (HWO) will provide new information to better understand the origin and the atmospheric effects of these precipitations. A major interest is for the distant magnetospheres of Uranus and Neptune, never explored by an orbital spacecraft whose UV auroral emissions remains at (Uranus) or below (Neptune) the HST sensitivity. \textit{Pollux} is one such UV instrument concept, which will enable unprecedented high spectral resolution at fine spatial scale not previously seen and polarimetric observations of the planetary aurorae while \textit{LUMOS}, another UV instrument will image the full auroral regions with a good spectral resolution.

The Solar Ultraviolet Imaging Telescope (SUIT) observes the Sun in the near-ultraviolet regime on board the Aditya-L1 satellite, India's dedicated mission to study the Sun. SUIT will image the Sun in the wavelength range of 200-400 nm using 11 science bandpasses with varying spectral bandwidths between 0.1-58 nm. Within this range, the Sun provides huge incoming solar flux to the telescope that also varies by a factor of ~ 20 from the lower end to the upper end of the wavelength band of interest. Thermal Filter Assembly (TFA) is an optical component at the SUIT entrance aperture, directly facing the Sun. The TFA is used to control the heat load entering the telescope cavity and also to reduce the signal reaching the SUIT camera system and the charge-coupled device (CCD) sensor, which is limited in full-well capacity and the linear operational regime. The TFA is designed to allow only 0.1-0.45% of the incoming flux to pass within 200-400 nm. The choice of materials for substrate and coating for the filter poses several challenges in terms of contamination, corrosion/ oxidation and durability during the manufacturing process. Additionally, long-term exposure to harsh space environments and the formation of pinholes are other concerns. Direct exposure to the sun leads to a strong temperature gradient along the thickness of the filter. The design and assembly of the TFA are performed to avoid any thermo-elastic stress affecting optical performance. Different levels of qualification tests and the operation of SUIT in orbit for more than 14 months have confirmed the perfect working of the TFA. To the best of our knowledge, the design, development, and testing of such a rejection filter is the first of its kind for space telescopes in the near ultraviolet range.

Thermochemical convection in Earth's outer core is driven by the crystallisation of the inner core that releases latent heat and light elements. A key question in core dynamics is whether a stable layer exists just below the core-mantle boundary. Recent core convection simulations, accounting for CMB heterogeneities, propose locally stable regions (or regional inversion lenses, RILs) rather than a global layer, allowing both stable and unstable regions to coexist. In this study, we consider a suite of numerical simulations of thermal, chemical, and thermochemical convection models focussed on Ekman number ($E=10^{-5}$) with thermal and chemical flux Rayleigh numbers $\widetilde{Ra}_T=30-4000$ and $\widetilde{Ra}_C=30-100000$, and thermal and chemical Prandtl numbers $Pr_T=1$ and $Pr_\xi=10$. Analysis of purely chemical models reveals light element accumulation (LEA) below the CMB, resulting in either locally stable regions near the poles or global layers, depending on the strength of chemical forcing. These chemically stratified regions persist in our thermochemical models even if the thermal field is fully destabilising. The addition of a heterogeneous CMB heat flux leads to the formation of RILs driven by thermal stratification. Stable regions in these thermochemical models have varying locations, properties, and morphologies depending on whether thermal or chemical convection dominates. In the investigated parameter range, these RILs are O(100 km) thick, and their strength and thickness generally increase with the strength of thermal driving; they are comparatively less sensitive to the strength of chemical driving. Our simulations reveal a diverse range of possible stable regions and/or a global layer at the top of Earth's core, with a seismically plausible range of thickness and strength, which may also have a signature in geomagnetic observations.

Space-based photometry from missions such as TESS has revealed that many young delta Scuti stars exhibit regular high-frequency pulsation patterns. These pulsations provide a powerful means of inferring stellar properties, particularly ages, for pre-main-sequence and early main-sequence delta Scuti stars, for which traditional age-dating methods are poorly constrained. Realising this potential requires robust theoretical models that capture the complexities of stellar structure and evolution. We present a comprehensive grid of 30 million stellar pulsation models, computed using the MESA stellar evolution code and the GYRE oscillation code, tailored specifically to delta Scuti stars. The grid spans a wide range of masses, metallicities, and surface rotation velocities, and covers evolutionary phases from the early pre-main-sequence through the main sequence and into the post-main-sequence contraction phase. For each model, we compute hundreds of adiabatic pulsation frequencies for degrees l = 0-3, capturing p-modes, g-modes, f-modes, and their interactions through avoided crossings. Our analysis maps the behaviour of asteroseismic observables, including the large frequency separation (Delta nu) and the phase offset parameter (epsilon), across age, mass, and rotation. We investigate how these parameters change with evolutionary stage and revisit the scaling relations applicable to delta Scuti stars. This new model grid, which is publicly available, is designed to support the asteroseismology community in interpreting delta Scuti pulsations and in probing the evolution and internal structure of these stars. These improvements over previous model grids will allow for reliable age estimates and stronger constraints on stellar evolution pathways and the timing of planet formation across A- and F-type stellar populations.

We aim to deepen our understanding of the M-L transition on planetary-mass companions and isolated brown dwarfs, and search for evidence of possible differences between these two populations of objects. To this end, we present a set of 21 VLT/SINFONI K-band observations from five archival programs at a spectral resolution of 4000. We aim to measure atmospheric properties, such as temperature, surface gravity, and metallicity, to understand the similarities and differences between objects ranging from M5 to L5 in spectral type. We extracted the spectra of these targets with the TExTRIS code. Subsequently, we model them using ForMoSA, a Bayesian forward modeling tool for spectral analysis, exploring four families of self-consistent atmospheric models: ATMO, BT-Settl, Exo-REM, and Sonora. Here we present the spectra of our targets and the derived parameters from the atmospheric modeling. We observed a drop in effective temperature of more than 500 K as a function of spectral type at the M/L transition, likely related to limitations in the current atmospheric models. In addition, we report carbon-to-oxygen ratio measurements for three companions (2M 0103 AB b, AB Pic b, and CD-35 2722 b), which contribute to the growing list of exoplanets for which this value has been measured. In conclusion, the VLT/SINFONI Library highlights two key points. First, there is a critical need to further investigate the discrepancies among grids of spectra generated by self-consistent models, as these models yield varying results and do not uniformly explore the parameter space. Second, we do not observe apparent discrepancies in the K-band spectra between companions and isolated brown dwarfs, which potentially suggests that these super-Jupiter objects formed through a similar process; however, this warrants further investigation.

In the late stage of their evolution, low- to intermediate-mass stars pass through the asymptotic giant branch (AGB) phase, characterised by strong mass loss through dust driven winds. High angular resolution observations reveal that these winds harbour strong deviations from spherical symmetry, such as spirals and arcs, believed to be caused by hidden (sub-)stellar companions. Much more often, one observes spectral lines, where the presence of a companion is less clear. We study the impact of a binary companion on low-J CO spectral lines of AGB star outflows. By varying the orbital separation and wind velocity, we aim to find line shapes characteristic of more complex binary-induced morphologies. We generated a grid of nine 3D models of a mass-losing AGB star using the smoothed particle hydrodynamics code Phantom, with three values for both the outflow velocity and orbital separation. Utilising the radiative transfer code Magritte, we created synthetic spectral lines for the low rotational transitions of CO at different inclinations and position angles. Our simulations show a variety of morphologies, always with a pronounced spiral structure arising in the orbital plane, but with varying shapes in the meridional plane, and different degrees of global flattening. We find that the CO line profiles can deviate strongly from the parabolic or flat-topped profiles expected from spherically symmetric outflows. A variety of line shapes emerge, with two peaks near the terminal velocity, and a central bump near the central velocity being the most pronounced. In specific cases, the spectral lines can appear parabolic, hiding the presence of a binary companion. We find the CO spectral lines can serve as a binary diagnostic. The influence of the companion on the line can however also go easily unnoticed, as the features can be concealed by the beam profile and the noise of the observations.

We show that multi-field inflationary models with negligible turning in field space during inflation can lead to an effective sourcing of adiabatic from entropic perturbations {\it after} the end of inflation. We illustrate this general phenomenon with a detailed analysis of an inflationary model whose scalar potential is determined by modular invariance. Its entropic perturbations are frozen during inflation, but instead, they are converted into adiabatic perturbations in the first post-inflationary $e$-folds. The curvature power spectrum, giving rise to CMB fluctuations, reaches a novel and enhanced plateau in this process; we address the implications for the inflationary observables $A_{s}$, $n_{s}$ and $r$.

The processes governing the production of astrophysical high-energy neutrinos are still debated, and the sources originating them remain an open question. Among the putative emitters, active galactic nuclei have gained increasing attention. Blazars, in particular, stand out due to their ability to accelerate particles in environments with external radiation fields. Recent observations suggest they may contribute to the neutrino flux detected by IceCube. We study the physical properties of a subsample of 52 blazars proposed as candidate neutrino emitters, based on a positional cross-correlation analysis between IceCube hotspots and the 5BZCat catalog. We aim to provide a first characterization of their central engines and physical nature, to explore the potential link with neutrino production. We analyze the optical spectroscopic properties of the 52 candidate neutrino-emitter blazars to infer their accretion regime. The study is complemented by radio and $\gamma$-ray data, which trace the intrinsic jet power. We compare the sample to other blazar populations in the literature, perform statistical tests, and explore, through simulations, the applicability of methods that include censored data. Overall, the target sample shows properties compatible with the reference samples. We observe a mild tendency to prefer objects with intense radiation fields, typical of radiatively efficient accretors, and high radio power. Among them, 24 are detected by Fermi-LAT, spanning various $\gamma$-ray luminosities. We also show that statistical tests commonly used in the literature need to be handled with caution, as they are sensitive to the number of censored data and the sample size.

As a phenomenon that occurs on Earth and on Mars, the diameter of a dust devil helps determine the amount of dust the devil injects into the atmosphere for both worlds -- for a given dust flux density (dust lifted per area per time), a wider devil will lift more dust into the air. However, the factors that determine a dust devil's diameter $D$ and how it might relate to ambient conditions have remained unclear. Moreover, estimating the contribution to an atmospheric dust budget from a population of dust devils with a range of diameters requires an accurate assessment of the differential diameter distribution, but considerable work has yet to reveal the best representation or explain its physical basis. In this study, we propose that this distribution follows a power-law $\propto D^{-5/3}$ and provide a simple physical explanation for why the distribution takes this form. By fitting diameter distributions of martian dust devil diameters reported in several studies, we show that the data from several studies support this proposed form. Using a previous model that treats dust devils as thermodynamic heat engines, we also show that the areal density of dust devils (number per unit area) $N_0$ scales with the product of their thermodynamic efficiency $\eta$ and the sensible heat flux $F_{\rm s}$ as $N_0 \propto \eta F_{\rm s}$.

We perform a comprehensive wide-band ($3-100$keV) spectro-temporal analysis of 13 outbursting BH-XRBs, using data (quasi)simultaneous with radio observations to unravel the complex disc-jet connection. RXTE observations are analyzed for XTEJ1859+226, GX339-4 (2002, 2006, and 2010 outbursts), 4U1543-47, H1743-322 (2003 and 2009 outbursts), XTEJ1550-564, XTEJ1752-223, XTEJ1650-500, SwiftJ1753.5-0127, XTEJ1748-288, and GROJ1655-40. For SwiftJ1727.8-1613 and MAXIJ1535-571, we utilize HXMT data, while both AstroSat and HXMT observations are analyzed for SwiftJ1658.2-4242. Type-C QPOs observed in harder states (LHS, HIMS; $F_{nth}\ge0.4$) exhibit positive lag for low-inclination sources ($i<50^{\circ}$), whereas it generally exhibits negative lag for high-inclination sources ($i>60^{\circ}$), except XTEJ1550-564, SwiftJ1727.8-1613, H1743-322 (2003 outburst) and GROJ1655-40. Notably, type-A QPOs exhibit negative lags ($\sim1-10$ms) regardless of source inclination, while type-B QPOs show positive lags in low-inclination sources, and both positive and negative lags ($\sim1-15$ms) in high-inclination sources, typically occurring in SIMS ($F_{nth}\lesssim0.45$). Systematic appearance of type-A QPOs preceding radio flares in several sources suggests that type-A QPOs indicate telltale signs of jet ejection, while type-B QPOs are closely linked with radio flares (i.e., transient jets). Present findings suggest the corona evolves from a radially extended to a vertically elongated structure during the type-C to type-B transition via type-A QPOs, with type-B QPOs linked to radially compact or vertically extended coronal geometries, resembling jet ejection. The strong radio-X-ray luminosity correlation seems to provide compelling evidence of accretion-powered jets. Finally, we find that jets in SIMS are moderately relativistic in nature with velocities $\gtrsim 0.3-0.8c$ in BH-XRBs under consideration.

The gravitational interaction between the Milky Way (MW) and the Large Magellanic Cloud (LMC) perturbs the MW halo's density and kinematics, encoding information about both galaxies' masses and structures. We present a suite of 2,848 high-resolution ($10^7$ particles) N-body simulations that systematically vary the mass and shape of both galaxies' haloes. We model how the mean velocities and velocity dispersions of halo stars (30--120 kpc) depend on system parameters, and forecast constraints achievable with current and future observations. Assuming Gaia DR3-level astrometry, 20 km/s radial velocity precision, 10% distance precision, and a sample of $\sim$4,000 RR Lyrae stars, we achieve 1$\sigma$ uncertainties of $0.11 \times 10^{12} M_\odot$ in MW mass, $2.33 \times 10^{10} M_\odot$ in LMC mass, 2.38 in halo concentration ($c$), and 0.06 in halo flattening ($q$). These correspond to fractional uncertainties of 11%, 16%, 25%, and 6% respectively, relative to fiducial values. Improved Gaia proper motions (DR5) yield modest gains (up to 14%), while adding radial velocities improves constraints by up to 60% relative to using Gaia astrometry alone. Doubling the sample size to $\sim$8,000 stars yields an additional 30% improvement, whereas reducing distance uncertainties has minimal impact ($\le$10%). Mean velocities trace LMC-induced perturbations, while velocity dispersions constrain MW halo properties, jointly breaking degeneracies. Our results demonstrate that combining Gaia astrometry with large spectroscopic surveys will enable precise characterization of the MW-LMC system. This methodology paper establishes the framework for interpreting observations; future work will apply these tools to existing spectroscopic datasets. The full simulation suite, HaloDance, will be made publicly available at: this https URL.

I present an analysis of archival spectra of 200 sources toward the Orion Nebula Cluster (ONC) that were obtained with the Near-Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST). I have used these data to assess cluster membership and measure spectral types for the targets. Fifty-three sources are classified as likely cluster members, 24 of which have spectral types that are suggestive of brown dwarfs ($>$M6). Seven of the NIRSpec targets were previously identified as "Jupiter-mass binary objects" (JuMBOs), all of which are background sources rather than brown dwarfs based on the NIRSpec data. The spectral classifications of those objects are consistent with the results of my recent study of the JWST photometry in the ONC, which found that only a few JuMBO components have the colors expected for brown dwarfs, none of which form pairs that have uniquely wide separations or low masses relative to known binary brown dwarfs.

Recent measurements of cosmogenic $^{10}$Be in deep-ocean ferromanganese crusts from the Central and Northern Pacific have revealed an anomalous concentration between 11.5 and 9.0 Myr ago, peaking at 10.1 Myr. One possible explanation is a nearby supernova (SN) event. Motivated by this and by the proximity of the Solar System to the Orion star-forming region during that period, we estimate the probability that at least one SN occurred between the onset and peak of the anomaly. Using an open cluster catalog based on Gaia DR3, we trace back the orbits of 2725 clusters and the Sun over the past 20 Myr and compute the expected number of SN events. We find 19 clusters with a probability greater than 1% each of producing at least one SN within 100 pc of the Sun in the time interval 11.5-10.1 Myr ago. The total cumulative probability exceeds zero at 35 pc from the Sun and increases rapidly with distance, reaching 68% near 100 pc. Two young clusters dominate the SN probability: ASCC 20 contributes most within 70 pc, while OCSN 61 becomes more significant beyond that distance. Our results support the plausibility of a SN origin for the $^{10}$Be anomaly and highlight the importance of additional $^{10}$Be records from independent terrestrial archives to determine whether the anomaly is of astrophysical or terrestrial origin.

We have located archival observations of the centaur 2015 OU$_{194}$ from 2017 and 2018, which extend its data-arc length from 1.0 to 3.5 years. We show that it is in an outer $3:4$ mean motion resonance with Uranus, henceforth referred to as U$_{3/4}$. The resonance is stable from at least 1000 kyrs in the past till 500 kyrs in the future. We find no mention in the literature of known objects in this resonance, or in any other resonance between the orbits of Uranus and Neptune. Looking for additional candidates, we find that 2013 RG$_{98}$ also stays in U$_{3/4}$ for several hundred kyrs around the present epoch. A third candidate, 2014 NX$_{65}$, is strongly influenced by Neptune.

Context. The dust and gas temperature in proto-planetary disks play critical roles in determining their chemical evolution and influencing planet formation processes. Aims. We attempted an accurate measurement of the dust and CO temperature profile in the edge-on disk of the Flying Saucer. Methods. We used the unique properties of the Flying Saucer, its edge-on geometry and its fortunate position in front of CO clouds with different brightness temperatures to provide independent constraints on the dust temperature. We compared it with the dust temperature derived using the radiative transfer code DiskFit and the CO gas temperature. Results. We find clear evidence for a substantial gas temperature vertical gradient, with a cold (10 K) disk mid-plane and a warmer CO layer where T(r) is 27 K at 100 au, dropping with exponent 0.3. Direct evidence for CO depletion in the mid-plane, below about 1 scale height, is also found. At this height, the gas temperature is 15-20 K, consistent with the expected CO freeze out temperature. The dust disk appears optically thin at 345 GHz, and exhibits moderate settling.

The large-scale structure survey J-PAS is taking data since October 2023. In this work, we present a forecast based on the Fisher matrix method to establish its sensitivity to the sum of the neutrino masses. We adapt the Fisher Galaxy Survey Code (FARO) to account for the neutrino mass under various configurations applied to galaxy clustering measurements. This approach allows us to test the sensitivity of J-PAS to the neutrino mass across different tracers, with and without non-linear corrections, and under varying sky coverage. We perform our forecast for two cosmological models: $\Lambda CDM + \sum m_\nu$ and $w_0w_a CDM + \sum m_\nu$. We combine our J-PAS forecast with Cosmic Microwave Background (CMB) data from the Planck Collaboration and Type Ia supernova (SN) data from Pantheon Plus. Our analysis shows that, for a sky coverage of 8,500 square degrees, J-PAS galaxy clustering data alone will constrain the sum of the neutrino masses to an upper limit of $\sum m_\nu < 0.32$\;eV for the $\Lambda CDM + \sum m_\nu$ model, and $\sum m_\nu < 0.36$\;eV for the $w_0w_a CDM + \sum m_\nu$ model. When combined with Planck data, the upper limit improves significantly. For J-PAS+Planck, we find $\sum m_\nu < 0.061$\;eV for the $\Lambda CDM + \sum m_\nu$ model, and for J-PAS+Planck+Pantheon Plus, we obtain $\sum m_\nu < 0.12$\;eV for the $w_0w_a CDM + \sum m_\nu$ model. These results demonstrate that J-PAS clustering measurements can play a crucial role in addressing challenges in the neutrino sector, including potential tensions between cosmological and terrestrial measurements of the neutrino mass, as well as in determining the mass ordering.

With perihelia well beyond Neptune, but semimajor axes and eccentricities indicative of substantial perturbation, the origins of detached trans-Neptunian objects (TNOs) remain a dynamical puzzle. In particular, detached TNOs with orbital inclinations below ~25 degrees are not easily generated from any known mechanism currently in the modern solar system. One notable hypothesis for the origins of detached TNOs is that a ~Mars- to Earth-mass planetary embryo detached the perihelia of these objects from Neptune during the process of Kuiper belt formation before the embryo itself was ejected. We numerically model this scenario via simulations of Kuiper belt formation from a primordial planetesimal belt that is dispersed through the migration of the giant planets. In addition to ~100,000 Kuiper belt objects, each of our simulations contains a hypothetical population of embryos in the primordial belt. We find that our embryos are unlikely to reach the high-perihelion, large semimajor axis orbit necessary to efficiently detach TNO perihelia from Neptune's influence. Moreover, embryos will typically take at least 100 Myrs to reach these unlikely orbits, at which point most of the primordial belt will have already been ejected by the planets, limiting the available population that can be detached. Finally, the TNOs that our embryos do detach consistently have a semimajor axis distribution that is more biased toward small values than observed detached TNOs have. Thus, we conclude that planetary embryos in the primordial Kuiper belt are not likely to have been the primary mechanism for the origin of detached TNOs.

The gravitational potential of the Milky Way encodes information about the distribution of all matter -- including dark matter -- throughout the Galaxy. Gaia data release 3 has revealed a complex structure that necessitates flexible models of the Galactic gravitational potential. We make use of a sample of 5.6 million upper-main-sequence stars to map the full 3D gravitational potential in a one-kiloparsec radius from the Sun using a data-driven approach called ``Deep Potential''. This method makes minimal assumptions about the dynamics of the Galaxy -- that the stars are a collisionless system that is statistically stationary in a rotating frame (with pattern speed to be determined). We model the distribution of stars in 6D phase space using a normalizing flow and the gravitational network using a neural network. We recover a local pattern speed of $\Omega_p = 28.2\pm0.1\mathrm{\,km/s/kpc}$, a local total matter density of $\rho=0.086\pm0.010\mathrm{\,M_\odot/pc^3}$ and local dark matter density of $\rho_\mathrm{DM}=0.007\pm0.011\mathrm{\,M_\odot/pc^3}$. The full 3D model exhibits spatial fluctuations, which may stem from the model architecture and non-stationarity in the Milky Way.

We investigate whether the latest combination of DESI DR2 baryon acoustic oscillation (BAO) measurements, cosmic microwave background (CMB) data (Planck 2018 + ACT), and Type Ia supernovae (SNe Ia) compilations (Pantheon+, Union3, and DES Y5) favor a dynamical dark energy component, and explore if such a scenario can simultaneously help resolve the Hubble tension. We contrast two frameworks: the widely used phenomenological $w_0 w_a$CDM model, and bimetric gravity, a fundamental modification of general relativity that naturally gives rise to phantom dark energy. The $w_0 w_a$CDM model is moderately preferred over $\Lambda$CDM, at the $2$-$4 \, \sigma$ level, when fitting DESI DR2 + CMB + SNe Ia, but it exacerbates the Hubble tension. By comparison, bimetric gravity provides a modest improvement in fit quality, at the $1 \, \sigma$ level, but, by inferring $H_0 = 69.0 \pm 0.4 \, \mathrm{km/s/Mpc}$, it partially eases the Hubble tension, from a $5 \,\sigma$ discrepancy to a $3.7 \, \sigma$ tension. Including locally calibrated SNe Ia brings the overall preference for the bimetric model over $\Lambda$CDM to the $2 \, \sigma$ level, comparable to that of the $w_0 w_a$CDM model when including the local SN Ia calibration.

The spectral resolution ($R \equiv \lambda / \Delta \lambda$) of spectroscopic data is crucial information for accurate kinematic measurements. In this letter, we present a robust measurement of the spectral resolution for the JWST's Near Infrared Spectrograph (NIRSpec) in all its modes: fixed slit (FS), integral field spectroscopy (IFS), and multi-object spectroscopy (MOS). We modeled H and He lines of the planetary nebula SMP LMC 58 using a Gaussian line spread function (LSF) to estimate the wavelength-dependent resolution for multiple disperser and filter combinations. We corrected for the intrinsic width of the planetary nebula's H and He lines due to its expansion velocity by measuring it from a higher-resolution X-shooter spectrum. We find that NIRSpec's in-flight spectral resolutions exceed the pre-launch estimates provided in the JWST User Documentation by 5-45% across all the configurations and covered wavelengths. We recover the expected trend that the resolution increases with the wavelength within a configuration. The robust and accurate LSFs presented in this letter will enable high-accuracy kinematic measurements using NIRSpec for applications in cosmology and galaxy evolution.

We show that the `shot noise' bias in angular clustering power spectra observed from discrete samples of points is not noise, but rather a known additive contribution that naturally arises due to degenerate pairs of points. In particular, we show that the true shot noise contribution cannot have a `non-Poissonian' value, even though all point processes with non-trivial two-point statistics are non-Poissonian. Apparent deviations from the `Poissonian' value can arise when significant correlations or anti-correlations are localised on small spatial scales. However, such deviations always correspond to a physical difference in two-point statistics, not a difference in noise. In the context of simulations, if clustering is treated as the tracer of a discretised underlying density field, any sub- or super-Poissonian sampling of the tracer induces such small-scale modifications and vice versa; we show this explicitly using recent innovations in angular power spectrum estimation from discrete catalogues. Finally, we show that the full covariance of clustering power spectra can also be computed explicitly: it depends only on the two-, three- and four-point statistics of the point process. The usual Gaussian covariance approximation appears as one term in the shot noise contribution, which is non-diagonal even for Gaussian random fields and Poisson sampling.

We present a science case for the Habitable Worlds Observatory (HWO) to map the circumgalactic medium (CGM) in emission by targeting ultraviolet emission lines, which trace the 10^4 - 10^6 K gas engaged in the feedback and accretion mechanisms driving galaxy evolution. While the CGM-galaxy connection is clearly evident through absorption line experiments and limited work done in optical and radio emission from the ground, the nature of this connection is poorly understood with regard to how these cosmic ecosystems exchange matter and energy. We outline a two-pronged experiment with Habitable Worlds Observatory (HWO) utilizing both multi-object spectroscopy to map kpc-scale CGM structures such as galactic superwinds and integral field spectroscopy to unveil the sub-kpc-scale processes such as thermal instabilities, which are theorized to govern the cool gas reservoirs long detected in pencil-beam absorption surveys. This article is an adaptation of a science case document developed by HWO's CGM/IGM Working Group.

The origins of the colors of Trans-Neptunian Objects (TNOs) represent a crucial unresolved question, central to understanding the history of our Solar System. Recent observational surveys have revealed correlations between the eccentricity and inclination of TNOs and their colors. This has rekindled the long-standing debate on whether these colors reflect the conditions of TNO formation or their subsequent collisional evolution. In this study, we address this question with 98.7% certainty, using a model-agnostic, data-driven approach based on causal graphs. First, as a sanity check, we demonstrate how our model can replicate the currently accepted paradigms of TNOs' dynamical history, blindly and without any orbital modeling or physics-based assumptions. In fact, our causal model (with no knowledge of the existence of Neptune) predicts the existence of an unknown perturbing body, i.e., Neptune. We then show how this model predicts, with high certainty, that the color of TNOs is the root cause of their inclination distribution, rather than the other way around. This strongly suggests that the colors of TNOs reflect an underlying dynamical property, most likely their formation location. Moreover, our causal model excludes formation scenarios that invoke substantial color modification by subsequent irradiation. We therefore conclude that the colors of TNOs are predominantly primordial.

We investigate the role of tidal forces in molecular cloud formation by examining how apparent boundedness, as diagnosed by the classical virial parameter, relates to the actual gravitational state of clouds subject to tidal forces from their environment. Clouds are identified by a dendrogram algorithm in zoom-in regions taken from a simulation of a Milky Way-mass galaxy with the Voronoi mesh code AREPO that resolves star-forming regions at sub-parsec resolution. To look at a range of environments, we use data from three different regions that evolve differently in the center, near the equivalent of the Solar circle, and the outskirts of the modeled galaxy, at three different times, each spaced 2 Myr apart. We compute the importance of tidal forces on all identified clouds. We then compare the boundedness of clouds including only their internal potentials to boundedness also including the external gravitational potential. This comparison shows that tidal forces can unbind apparently bound clouds and bind apparently unbound clouds. We characterize the cloud population by comparing their virial parameters to their surface densities, finding the ratio of the maximum to the minimum eigenvalues of the tidal tensor, and determining the strength of gravitational instability in each examined region. We find that it is necessary to take the total gravitational potential into account rather than just the internal self-gravity of the clouds to have an accurate understanding of cloud dynamics.

The origins and acceleration mechanisms of ultra-high-energy cosmic rays (UHECRs) are unknown. Many models attribute their extreme energies to powerful astrophysical jets. Understanding whether jet geometry -- specifically the opening angle and its orientation relative to Earth -- affects observational signatures, is crucial for interpreting UHECR data. In this work, we perform numerical simulations of UHECR propagation in a magnetized universe to investigate the spectral signatures of jetted and non-jetted astrophysical sources. We demonstrate, for the first time, that under certain conditions, emission geometry can play a decisive role in shaping the observed spectrum of individual UHECR sources. These findings provide new insights into the conditions necessary for detecting UHECRs from jets, and highlight how the interplay between emission geometry and magnetic fields influences observed energy spectra.

Our Earth, being the only living planet that we know, provides us with clues that photosynthetic life-forms may be dominant on other exoplanets for billions of years. Spectropolarimetric signatures of the terrestrial photosynthetic life (PSLife) are well studied in the lab and remotely sensed with space and airborne instrumentation. An astonishing biosignature revealed by these measurements is an extremely strong linear polarization (tens \%) associated with broad absorption bands of biological pigments (biopigments) driving photosynthesis in various organisms. Also, unique circular-polarization signatures are associated with biopigments and other complex macromolecules as a sign of homochirality which is ubiquitous in terrestrial life forms. Thus, low-resolution spectro- or multi-band polarimetry of exoplanets directly imaged at an unprecedented contrast using the HWO coronagraph is a novel opportunity for a robust discovery of life on exoplanets. Here we propose to carry out two surveys and two follow-up observing programs. Survey 1 will identify potentially habitable planets (PHPs) through detection of atmospheres, clouds and liquid surface water (ocean) using linear polarimetry. Survey 2 will identify Living World (LW) candidates among PHPs by searching for strong linear polarization signatures associated with strong and broad absorption bands reminiscent of terrestrial biopigments. Follow-up program 3 will obtain multi-color surface maps of LWs, determine the distribution and abundance of alien photosynthetic organisms with exo-biopigments (exoBPs) and correlate their properties with the atmospheric and surface compositions. Follow-up program 4 will employ circular polarization to verify homochirality of exoBPs. This comprehensive approach aims at providing a quantitative answer to the ultimate question "Are we are alone in the Universe?".

In recent years, a class of stripped-envelope supernovae (SESNe) showing two distinct light-curve peaks has emerged, where the first peak cannot be attributed to shock cooling emission. Such peculiar SNe are often studied individually, explained by a combination of powering mechanisms, but are rarely discussed broadly as a group. In this paper, we attempt to form a picture of the landscape of double-peaked SESNe and their powering mechanisms by adding two more objects -- SN 2021uvy and SN 2022hgk. SN 2021uvy is a broad, luminous SN Ib with an unusually long first peak rise and constant color evolution with rising photospheric temperature during the second peak. Though its first peak resembles SN 2019stc, their second peaks differ, making SN 2021uvy unique. SN 2022hgk shows photometric similarity to SN 2019cad and spectroscopic similarity to SN 2005bf, both proposed to be powered by a double-nickel distribution in their ejecta. We analyze their light curves and colors, compare them with a sample of double-peaked SESNe from the ZTF archive, and analyze the light curve parameters of the sample. We observe a correlation (p-value~0.025) between the peak absolute magnitudes of the first and second peaks. No single definitive powering mechanism applies to the whole sample, as it shows variety in the photometric and spectroscopic properties. However, sub-groups of similarity exist that can be explained by mechanisms like the double-nickel distribution, magnetar central engine, interaction, and fallback accretion. We also map out the duration between the peaks ($\Delta t^{21}$) vs the difference between peak absolute magnitudes ($\Delta M^{21}$) as a phase-space that could potentially delineate the most promising powering mechanisms for the double-peaked SESNe.

A proper modelling of the cosmic-ray ionisation rate within gas clouds is crucial to describe their chemical evolution accurately. However, this modelling is computationally demanding because it requires the propagation of cosmic rays throughout the cloud over time. We present a more efficient approach that simultaneously guarantees a reliable estimate of the cosmic-ray impact on the chemistry of prestellar cores. We introduce a numerical framework that mimics the cosmic-ray propagation within gas clouds and applies it to magnetohydrodynamic simulations performed with the code GIZMO. It simulates the cosmic-ray attenuation by computing the effective column density of H$_2$ that is traversed, which is estimated using the same kernel weighting approach as employed in the simulation. The obtained cosmic-ray ionisation rate is then used in post-processing to study the chemical evolution of the clouds. We found that cosmic-ray propagation affects deuterated and non-deuterated species significantly and that it depends on the assumed cosmic-ray spectrum. We explored correlations between the electron abundance, the cosmic-ray ionisation rate, and the abundance of the most relevant ions (HCO$^+$, N$_2$H$^+$, DCO$^+$, N$_2$D$^+$, and o-H$_2$D$^+$), with the purpose of finding simple expressions that link them. We provide an analytical formula to estimate the ionisation fraction, X(e$^-$), from observable tracers and applied it to existing observations of high-mass clumps. We obtained values of about 10$^{-8}$, which is in line with previous works and with expectations for dense clouds. We also provide a linear fit to calculate the cosmic-ray ionisation rate from the local H$_2$ density, which is to be employed in three-dimensional simulations that do not include cosmic-ray propagation.

Planets in the liquid-water habitable zone of low-mass stars experience large tidal forces, $10^3$ to $10^4$ times those on Earth, due to the small distance between the habitable zone and the host stars. Therefore, interior solid tides, ocean tides and atmospheric tides on these planets could be much stronger than that on Earth, but rare work has been done to explicitly simulate the ocean tides. Here, for the first time, we perform global ocean tide simulations and show that ocean tides on asynchronously rotating planets with large eccentricities can reach $\mathcal{O}(1000)\,\mathrm{m}$ in height and $\mathcal{O}(10)\,\mathrm{m\,s^{-1}}$ in flow speed. Interactions between tide and bottom topography can induce large energy dissipation, $\sim\mathcal{O}(100)\,\mathrm{W\,m^{-2}}$ in global mean. This tidal energy dissipation can strongly accelerate orbital evolution by 1-2 orders of magnitude. However, for planets with small eccentricities, the ocean tides are much weaker but still comparable to that on modern Earth. Our results suggest that ocean tides on eccentric planets orbiting low-mass stars are orders of magnitude more powerful than those on Earth and can dramatically influence surface geography and orbital evolution.

CMB lensing reconstructions are a sensitive probe of the growth of structure across cosmic time and a key tool to sharpen investigations of the very early Universe via delensing. At present, a large fraction of this information is drawn from the temperature anisotropies, which are ultimately also the most informative when reconstructing lenses on arcminute scales and smaller. But extragalactic foreground emission from galaxies and clusters can contaminate these reconstructions, limiting our ability to use information from small-scale temperature anisotropies. We develop analytic predictions of the biases from the thermal Sunyaev-Zeldovich and cosmic infrared background to CMB lensing auto- and cross-correlations with low-redshift matter tracers, as well as B-mode delensing, based on a halo model for the dominant one- and two-halo contributions to the relevant foreground bi- and tri-spectra. The method is flexible enough to allow variations in cosmology, astrophysical modeling, experimental configurations and analysis choices, thus enabling an improved understanding of the uncertainties involved in current mitigation strategies. We find that the shape of the bias relative to the CMB lensing auto-spectrum signal is remarkably insensitive to changes in cosmological and astrophysical parameter values. On the other hand, the shape appears to depend on $\Omega_m$ for cross-correlations with low-redshift galaxies. We also clarify the ranges of redshifts and masses that simulations need to resolve in order to capture these effects accurately. Our code, CosmoBLENDER, is made publicly available.

In this paper, we report auroral radio emission from a magnetic B star HD 142990 using the MeerKAT radio telescope at $900-1670$ MHz. The star is known to produce such emission (observed as periodic radio pulses) via electron cyclotron maser emission (ECME). However, past studies on ECME from this star were confined to observations at specific rotational phase ranges where one expects to see such pulses. We, for the first time, observed the star for its one complete rotation cycle and discovered that the star also produces 'off-pulse' emission, which we term as secondary enhancements. Two such enhancements were observed, one of which is left circularly polarized (LCP) and the other is right circularly polarized (RCP), the latter is confirmed to be persistent. Using simulation, we infer that such pulses are likely related to the large misalignment between the stellar rotation and magnetic dipole axes ($>80^\circ$), leading to the formation of highly complex magnetospheric plasma distribution. In addition, by extracting dynamic spectra for the primary pulses, we discovered prominent fine structures in one of the LCP pulses, with timescales as small as the instrumental time resolution (8 seconds). This is the first time that such structures are seen from a magnetic hot star, and has the potential to reveal detailed information about how the emission is driven, and the nature of the elementary sources of radiation. To pinpoint the origin of these fine structures and their significance, higher time and spectral resolution observations should be conducted in the future.

The standard formation model of close-in low-mass planets involves efficient inward migration followed by growth through giant impacts after the protoplanetary gas disk disperses. While detailed N-body simulations have enhanced our understanding, their high computational cost limits statistical comparisons with observations. In our previous work, we introduced a semi-analytical model to track the dynamical evolution of multiple planets through gravitational scattering and giant impacts after the gas disk dispersal. Although this model successfully reproduced N -body simulation results under various initial conditions, our validation was still limited to cases with compact, equally-spaced planetary systems. In this paper, we improve our model to handle more diverse planetary systems characterized by broader variations in planetary masses, semi-major axes, and orbital separations and validate it against recent planet population synthesis results. Our enhanced model accurately reproduces the mass distribution and orbital architectures of the final planetary systems. Thus, we confirm that the model can predict the outcomes of post-gas disk dynamical evolution across a wide range of planetary system architectures, which is crucial for reducing the computational cost of planet formation simulations.

We report the results of our long-term multiwavelength spectral energy distribution (SED) study on the flat spectrum radio quasar 3C 279 during the $\sim14$ years (2008--2022) of {\sl Fermi}-LAT (Large Area Telescope) observing period. The {\sl Fermi}-LAT data were complemented with data in other wavebands obtained from {\sl Swift}-XRT/UVOT, Whole Earth Blazar Telescope (WEBT), along with other optical and radio data from several observatories. Different activity states were identified from the weekly binned $\gamma$-ray light curve, and it was possible to create 168 high-quality and quasi-simultaneous broadband SEDs. We modeled the SEDs using a one-zone leptonic scenario, including the emission region outside the broad-line region (BLR), involving synchrotron, synchrotron self-Compton, and external Compton mechanisms. Such extensive broadband modeling is essential for constraining the underlying multiwavelength radiative mechanisms in the 3C 279 jet and permits to estimate the physical parameters and explore their evolution in time. Our SED modeling study suggests that the increase in the Doppler beaming factor along with the variation of the emitting electrons is the cause for the flares in this source. The multiwavelength emission of 3C 279 was found to be well explained by the scenario in which the emission region is outside the BLR at a distance of $\sim6.42\times10^{3} R_S$. However, for two of the very bright $\gamma$-ray states, the emission region was found to be close to the outer boundary of the BLR at a distance of $\sim1.28\times10^{3} R_S$ from the central black hole.

The variability of virial factor $f$ is investigated for two active galactic nucleus, NGC 5548 and NGC 4151, which had been previously reverberation mapped (RM) over 20 times in the past 30 years. Using four velocity tracers from the broad H$\beta$ width at half-maximum ($\rm FWHM_{\rm H\beta} $) or the line dispersion from the mean or rms spectra, $f$ for each RM epoch are calculated. Significant correlations are found between $f$ and observational parameters, such as the broad line widths, the Eddington ratios and the line profile shapes. For NGC 5548, $f \propto \rm {FWHM}_{mean}^{-0.70\pm0.13}$ and for NGC 4151, $f \propto \rm {FWHM}_{mean}^{-3.31\pm0.59}$. This suggests that a variable $f$ should be included to weight the virial SMBH mass. Using a simple model of thick-disc broad line regions (BLRs), we show that changes in mean inclination can explain $f$ variation. The inclination range is $14.1-40.6$ deg for NGC 5548 and $14.0-55.1$ deg for NGC 4151. Modeling the light curves of $f$ with a damped random walk process yields mean timescales of 638 and 668 days, consistent with BLR dynamical timescales within uncertainties. This indicates that $f$ variations are linked to BLR dynamics, likely due to changes in geometry or inclination.

In the context of the development of several space mission projects for UV spectropolarimetry at high resolution and over a wide UV wavelength range, such as Arago, Polstar, and Pollux onboard the Habitable Worlds Observatory, we are studying and developing the UV nanosatellite CASSTOR to obtain the very first UV spectropolarimetric observations of hot stars and test several new technologies, in particular a UV polarimeter and a Fine Guiding System. In this paper, we present the work and outcome of the Phase 0 study of CASSTOR.

The Lane-Emden equation, a nonlinear second-order ordinary differential equation, plays a fundamental role in theoretical physics and astrophysics, particularly in modeling the structure of stellar interiors. Also referred to as the polytropic differential equation, it describes the behavior of self-gravitating polytropic spheres. In this study, we present a novel approach to the solution of the eigenvalue problem which arises when considering the Lane-Emden equation for n = 0, 1, 2, 3, 4 using Physics-Informed Neural Networks (PINNs). The novelty of this work is that, we not only solve the Lane-Emden equation via PINNS but we also determine the eigenvalue, r, which is the stellar radius. Hyperparameter tuning was conducted using Bayesian optimization in the Optuna framework to identify optimal values for the number of hidden layers, number of neurons, activation function, optimizer, and learning rate for each value of n. The results show that, for n = 0, 1, PINNs achieve near-exact agreement with theoretical eigenvalues (errors < 0.000806%). While for more nonlinear cases, n = 2, 3 and n=4, PINNs yield errors below 0.0009% and 0.05% respectively, validating their robustness.

This study explores transient quasi-periodic oscillations (QPOs) in the $\gamma$-ray emission of two blazars, PMN J0531$-$4827 and PKS 1502+106, using over a decade of Fermi Large Area Telescope observations. The analysis focuses on identifying QPO signatures in their long-term light curves and interpreting the variability through a curved jet model, which predicts multiplicative oscillations with exponentially decaying amplitudes. We develop an analysis methodology to characterize the QPO and the specific properties of the amplitude of such QPOs. The findings offer insights into the dynamic processes driving relativistic jet evolution and their potential connections to underlying mechanisms, such as binary systems or other phenomena influencing the observed characteristics of these blazars.

The origin of the terrestrial water remains debated, as standard Solar System formation models suggest that Earth formed from dry grains, inside the snowline of the Proto-Solar Nebula (PSN). Here, we revisit this issue through the lens of computational chemistry. While the classically used snowline relies on a single condensation temperature, recent work in quantum chemistry shows that the binding energy of water on icy grains has a gaussian distribution, which implies a gradual sublimation of water rather than a sharp transition. We use the computed distribution of binding energies to estimate the radial distribution of adsorbed ice on the dust grains across the PSN protoplanetary disk. Our model reproduces the full range of estimated water abundances on Earth and matches the hydration trends observed in chondrite groups at their predicted formation distances. Thus, we suggest that a significant fraction of Earth's water may have been acquired locally at early stages of the Solar System formation, without requiring delivery from beyond the classical snowline.

We present a systematic search for spatial association between a high-confidence sample of 3,166 morphologically selected galaxy mergers detected through SDSS imaging and the 4FGL-DR4 catalog of gamma-ray sources detected by the Fermi Large Area Telescope (LAT). Using a conservative 4$\sigma$ positional uncertainty threshold and a Poisson-based statistical framework, we identify 21 statistically significant associations with match probabilities $p < 0.05$. Among these are known classes of gamma-ray emitters such as flat-spectrum radio quasars (FSRQs), BL Lacertae objects (BLL), and radio galaxies bolstering the hypothesis that merger-driven processes can fuel high-energy activity. Intriguingly, five of the associated sources remain unclassified in 4FGL-DR4, hinting at a possible link between galaxy mergers and a hitherto unrecognized population of gamma-ray sources. The dominance of AGN-like matches supports the scenario in which mergers trigger accretion onto central supermassive black holes, initiating AGN activity observable at gamma-ray energies. Moreover, the recurrent presence of unassociated sources among secure matches underscores the potential of merger catalogs as physically motivated priors in gamma-ray source identification efforts. This work constitutes the first dedicated effort to explore associations between gamma-ray sources and a large, morphologically selected sample of galaxy mergers, opening new avenues for understanding the role of interactions in high-energy astrophysics. We additionally examined the sample of 70 galaxy pairs from the Canadian Network for Observational Cosmology Field Galaxy Redshift Survey and found no statistically significant matches, with $p < 0.05$ within the 4$\sigma$ positional uncertainty threshold.

The powerful jets of blazars have been historically considered as likely sites of high-energy cosmic-ray acceleration. However, particulars of the launched jet and the locations of leptonic and hadronic jet loading remain unclear. In the case when leptonic and hadronic particle injection occur jointly, a temporal correlation between synchrotron radiation and neutrino production is expected. We use a first catalog of millimeter (mm) wavelength blazar light curves from the Atacama Cosmology Telescope for a time-dependent correlation with twelve years of muon neutrino events from the IceCube South Pole Neutrino Observatory. Such mm emission is known to trace activity of the bright jet base, which is often self-absorbed at lower frequencies and potentially gamma-ray opaque. We perform an analysis of the population, as well as analyses of individual, selected sources. We do not observe a significant signal from the stacked population. TXS 0506+056 is found as the most significant, individual source, though this detection is not globally significant in our analysis of selected AGN. Our results suggest that the majority of mm-bright blazars are neutrino dim. In general, it is possible that many blazars have lighter, leptonic jets, or that only selected blazars provide exceptional conditions for neutrino production.

One of JWST's most remarkable discoveries is a population of compact red galaxies known as Little Red Dots (LRDs). Their existence raises many questions about their nature, origin, and evolution. These galaxies show a steep decline in number density-nearly two orders of magnitude-from $z=6$ to $z=3$. In this study, we explore their potential evolution by identifying candidate descendants in CEERS, assuming a single evolutionary path: the development of a blue star-forming outskirt around the red compact core. Our color-magnitude selection identifies galaxies as red as LRDs at $z<4$, surrounded by young, blue stellar outskirts. Morphological parameters were derived from single Sérsic profile fits; physical properties were obtained from SED fitting using a stellar-only model. These "post-LRD" candidates show LRD-like features with $M_\ast \sim 10^{10} \ M_\odot $, central densities ($ \Sigma_\ast \sim 10^{11} \ M_\odot \ \text{kpc}^{-2}$ ), compact sizes, and red rest-frame colors, but with an added extended component. Their number density at $z = 3 \pm 0.5$ ( $ \sim 10^{-4.15} \, \text{Mpc}^{-3} $) matches that of LRDs at $5 < z < 7$ , supporting a possible evolutionary link. We observe a redshift-dependent increase in outskirts mass fraction and galaxy size-from $\sim 250$ pc at $ z = 5 $ to $\sim 600$ pc at $ z = 3 $-suggesting global stellar growth. Meanwhile, the core remains red and compact, but the V-shape fades as the outskirts grow. These findings support an evolutionary scenario in which LRDs gradually acquire an extended stellar component over cosmic time by cold accretion. This may explain the apparent decline in their observed number density at lower redshift.

Measuring the merger rate density history of binary neutron stars (BNS) can greatly aid in understanding the history of heavy element formation in the Universe. Currently, second-generation Gravitational Wave (GW) detectors can only measure the BNS merger rate density history at low redshifts ($z$ $\sim$ 0.1). Short gamma-ray bursts (sGRBs) may trace the BNS merger to higher redshifts ($z$ $\sim$ 3). However, not all BNS mergers result in sGRBs, and it is not certain that all sGRBs originate from BNS mergers. In this study, we simultaneously utilize simulated BNS merger GW signals detected by the advanced LIGO design and sGRB signals detected by {\it Fermi}/GBM to constrain the BNS merger rate density history up to $z$ $\sim$ 3. The results indicate that with $\sim$ 8 GWs and 571 sGRBs, the BNS merger rate density can be measured with an accuracy of about 50\% through $z=0$ to $z=1$. The ratio of the jet opening angle-corrected sGRB event rate density to the BNS merger rate density, denoted as $\eta$, can be constrained to a relative uncertainty of 45\%. With $\sim$ 21 GWs and 761 sGRBs, the BNS merger rate density can be measured to approximately 35\% and 40\% at $z=0$ and $z=1$, respectively. Meanwhile, $\eta$ can be constrained to a relative uncertainty of 28\%. Additionally, in our parameterized simulation, we find that at least approximately $\sim$550 sGRBs are needed to constrain the characteristic delay time in the star formation rate model, given a relative error of 50\% in the estimated redshift.

One important property in studying the exoplanet population is the host star metallicity ([M/H]). In this study, we derived stellar metallicities and oxygen abundances for 45 M dwarf stars (with 65 exoplanets) using the near-infrared high-resolution spectra from the SDSS APOGEE survey and synthetic spectra computed in LTE. We also investigated the exoplanetary radii distribution for a larger sample of 246 exoplanets orbiting 188 M dwarf stars. The [M/H] versus [O/M] distribution obtained indicates that our sample is composed mainly of thin disk stars, which follow the behavior of the low-alpha sequence in the Milky Way thin disk. Small planets with radii smaller than 3R$_{\oplus}$ were found around stars with a range of metallicities (-0.6$<$[M/H]$<$+0.3), while larger planets of the sample orbit only stars with [M/H]$>=0.0$. These results indicate that while small planets can form in different environments, larger planets preferentially form in metal-rich protoplanetary disks. Exoplanets with P$_{\rm orb}<$4.3 days orbit on average more metal-rich stars than planets with longer periods. This threshold is smaller than that found for FGK stars (8--10 days) and might be related to M dwarfs having a smaller dust sublimation radius. The distribution of exoplanets with R$_{\rm p}>$4R$_{\oplus}$ shows a concentration at orbital periods between 2 and 5 days, which may result from inward orbital migration. There is also a different behavior between single-detected exoplanets and planets from multiplanetary systems, with the latter being found on average around more metal-poor stars, and with orbital periods roughly up to 3 days.

An analytical treatment of the pulsar radio emission height as a function of phase is presented. Based on previous works, assuming a dipole field geometry, the emission height can be expressed as a function of phase and the impact angle. We found that: (1) The emission height is a quadratic function of the phase, given the magnetospheric geometry. The emission height is higher at the edge of the pulse profile than at the center. (2) The emission height is also quadratic function of the impact angle. This point can be compared directly with the observations of the geodetic precessing pulsar PSR J1906+0745. (3) Larger inclination angle may imply a higher emission height, when other parameters are similar. (4) By assuming curvature radiation, or inverse Compton scattering, the Lorentz factor of the radio emitting particles can be obtained as a function of phase.

The Siberian Radioheliograph (SRH) is a ground-based radio interferometer in Irkutsk, Russia, designed for high-resolution solar observations in the microwave range. It can observe dynamic solar events with spatial resolutions of 7-30 arcseconds and temporal resolution up to 0.1 seconds. Generating solar radio images from the Siberian Radioheliograph (SRH) is a multi-step calibration process that corrects instrumental and atmospheric distortions, using redundancy-based calibration with both adjacent and non-adjacent antenna pairs to address phase and amplitude errors in visibility data. The CLEAN algorithm is then applied to deconvolve the point spread function, reduce sidelobes, and enhance the visibility of solar features, resulting in higher quality and more reliable images. While the calibration process generally improves image quality, it can sometimes result in noisy or spatially shifted images that are not suitable for scientific use. We developed a deep learning approach for automatic image quality classification. The training dataset was prepared using a zero-shot CLIP model and further validated manually. We evaluated four different models: a fine-tuned EfficientNet, two CatBoost variants using embeddings from CLIP and EfficientNet, and an Ensemble model that combined predictions from all three individual models. The Ensemble model achieved the best performance. The SRH daily image classification service has been created and is available online at this https URL along with an API offering IDL and Python examples. Integration of Ensemble model into SRH image generating and calibration workflow can improve image reliability and reduces low-quality entries in SRH data catalog, enhancing solar research outcomes.

We present a comprehensive study of the galaxy size-stellar mass relation (SMR) at low redshift (z <= 0.125), using a large spectroscopic sample from the SDSS-DR13 survey. Our goal is to investigate how environment affects galaxy structural properties across multiple spatial scales. Galaxies are classified by specific star formation rate, optical color, and bulge-to-total light ratio, allowing us to disentangle environmental effects from intrinsic galaxy properties. We examine the SMR in three contexts: (1) comparing galaxy sizes in two extreme environments-dense clusters versus cosmic voids; (2) analyzing cluster galaxies across a range of cluster masses; and (3) studying member galaxies located in different cluster regions, from the core to the infall zone. In all three cases, we find no significant dependence of the SMR on environment at fixed stellar mass and galaxy type. Cluster and void galaxies follow consistent SMR trends, and no measurable variation is observed with cluster mass or cluster-centric distance. We also confirm that early-type galaxies exhibit steeper SMR slopes than late types. Notably, this consistent lack of environmental dependence on the SMR persists even when accounting for the differing galaxy number densities in voids, supporting the universality of this SMR scaling relation across diverse environments.

We present a torsion-based cosmological model within the Einstein-Cartan framework, analyze it using combined DESI, PantheonPlus/DESY5, and CMB datasets, a reveals several interesting features. The torsion parameter $\alpha$ shows a trend of change from $\alpha=0.042^{+0.035}_{-0.031}$ (CMB alone) to $\alpha=-0.0075\pm0.0058$ (full dataset), while remaining consistent with zero within $1.7\sigma$. The model demonstrates potential in addressing cosmological tensions, yielding $H_0 = 68.80\pm0.48~\mathrm{km~s^{-1}~Mpc^{-1}}$ and reducing the $S_8$ tension from $0.73\sigma$ to $0.1\sigma$ through its torsion-modified growth equation. Bayesian evidence analysis shows a clear preference for torsion cosmology in most dataset combinations, with the notable exception of CMB+DESI+DESY5 which slightly favors $\Lambda$CDM. Remarkably, substituting DESY5 with PantheonPlus data reverses this preference, yielding $\Delta\log\mathcal{Z} = -3.59$ - a strong indication of the model's viability. While not conclusive, these findings suggest torsion cosmology merits further investigation as a potential extension to $\Lambda$CDM.

The "Li-Dip" is an unexpected, striking, and highly non-standard anomaly of severe lithium depletion observed in mid-F dwarf stars, which has puzzled astronomers for nearly 40 years. Mechanisms proposed to explain the Li-Dip include effects related to rotation, magnetic fields, diffusion, gravity waves, and mass loss. The critical question became, which, if any, might be realistic? Here we show that mixing due to shear induced by stellar angular momentum loss is the unique mechanism driving the Li depletion. Each mechanism leaves a different signature in the subsurface Li distribution. The deepening surface convection zones of subgiants of NGC 188 evolving out of the Li-Dip dredge up the sub-surface material and thus reveal the signature of the responsible mechanism, rotation. Beryllium and boron data have also favored rotational mixing; however, these elements can be extremely difficult or impossible to observe. Our highly complementary approach provides fresh and very feasible perspectives on using Li to probe poorly understood physical mechanisms acting below the stellar surface, thereby improving fundamental understanding of stellar evolution. Rotational mixing may be the dominant mechanism that depletes Li in a wide range of Solar-type stars, including in the Sun. Possible connections to Big Bang Nucleosynthesis are discussed.

Context: The surface abundances of massive stars show evidence of internal mixing, while asteroseismic data suggest that efficient angular momentum (AM) transport occurs in stellar interiors. It is of interest to find a consistent physical framework that is able to account for both of these effects simultaneously. Aims: We investigate the impact of the Tayler instability on the surface abundance of boron in massive B-type stars as predicted by rotating stellar models accounting for the advective nature of meridional currents. Methods: We used the Geneva stellar evolution code to compute models of 9, 12, and 15 Msun stars at different rotational velocities and with and without magnetic fields. We compared the surface boron abundances predicted by these models with those of observed B-type stars. Results: We find that models with only hydrodynamic transport processes overestimate the amount of boron depletion for stars with high rotation rates, in disagreement with observational constraints. We show that this excessively high mixing efficiency is a consequence of the high degree of differential rotation predicted by purely hydrodynamic models. We thus conclude that surface abundances of boron indicate that a more efficient AM transport is needed in stellar radiative zones. We then studied the impact of the Tayler instability as a possible physical explanation to this issue. Models including this instability are found to be in good agreement with constraints on the surface boron abundances, the evolutionary state, and the projected rotational velocity of moderately and fast-rotating B-type stars. Finally, we note that at low rotational velocities, models with magnetic fields do not predict sufficient depletion to be consistent with the observations. This could suggest that the current prescriptions for the Tayler instability may overestimate the AM transport in slow-rotating B-type stars.

We investigate the wavevector and variance anisotropies in the inertial range of the young solar wind observed by the Parker Solar Probe (PSP). Using the first 19 encounters of PSP measurements, we identify the young solar wind from different source regions: coronal hole (CH) interiors, streamers, and low Mach-number boundary layers (LMBLs), i.e., the peripheral region inside CHs. We assess the wavevector anisotropy with the 2D and slab turbulence model for the CH wind and the streamer wind, and the nearly incompressible (NI) MHD turbulence model for the LMBL wind where Taylor's hypothesis becomes questionable. Unlike the $\sim80\%$ 2D contribution typically reported at 1 au, our results show that only $26\%$ of the inertial range energy is associated with 2D fluctuations in the CH wind, and this fraction increases to $45\%$ in the streamer wind. As a representation of the LMBL wind, similarly, the oblique sub-Alfvénic intervals and the near-subsonic intervals are characterized by the dominance of slab fluctuations. All the results suggest that slab fluctuations are more abundant in the young solar wind below 0.3 au than at 1 au. Furthermore, we find a dependence of the variance anisotropy in the inertial range on proton plasma beta $\beta_p$. The variance anisotropy is the strongest in the LMBL wind with the lowest $\beta_p$, and the weakest in the streamer wind with the highest $\beta_p$. This contrast can be interpreted as the remnant of fluctuations from the coronal sources.

We present the Rossiter-McLaughlin measurement of the sub-Neptune TOI-1759A b with MAROON-X. A joint analysis with MuSCAT3 photometry and nine additional TESS transits produces a sky-projected obliquity of $|\lambda|$= $4^\circ\pm18^{\circ}$. We also derive a true obliquity of $\psi$=24$\pm12^{\circ}$ making this planet consistent with full alignment albeit to $<1\sigma$. With a period of 18.85 days and an $a/R_{*}$ of 40, TOI-1759A b is the longest period single sub-Neptune to have a measured obliquity. It joins a growing number of smaller planets which have had this measurement made and, along with K2-25 b, is the only single, aligned sub-Neptune known to date. We also provide an overview of the emerging distribution of obliquity measurements for planets with R$<8$ R$_{\oplus}$. We find that these types of planets tend toward alignment, especially the sub-Neptunes and super-Earths implying a dynamically cool formation history. The majority of misaligned planets in this category have 4$<$R$\leq$8 R$_{\oplus}$ and are more likely to be isolated than planets rather than in compact systems. We find this result to be significant at the $3\sigma$ level, consistent with previous studies. In addition, we conduct injection and recovery testing on available archival radial velocity data to put limits on the presence of massive companions in these systems. Current archival data is insufficient for most systems to have detected a giant planet.

Long-duration GRB 211211A, which lacked an associated supernova at such a low redshift $z=0.076$, but was associated with a possible kilonova emission, has attracted great attention. The neutron star-white dwarf (NS-WD) merger is proposed as a possible progenitor of GRB 211211A, and it could naturally explain the long duration of the prompt emission. However, the NS-WD merger is not an ideal site for producing heavy elements via r-process nucleosynthesis. In this Letter, we investigate the heavy elements produced in NS-WD mergers based on numerical simulations of nucleosynthesis via SkyNet, and then calculate the resulting kilonova-like emission to compare with the solidly observed case of possible kilonova emission associated with GRB 211211A. By adopting three models (i.e., Model-A, Model-B, and Model-C) from \cite{2023ApJ...956...71K} at different temperatures ($T=4$ GK, 5 GK, and 6 GK), which are treated as free parameters, we find that the mass number of the heaviest element produced in our simulations is less than 90 ($A< 90$). Moreover, by comparing the calculated kilonova-like emission with the afterglow-subtracted observations of the possible kilonova associated with GRB 211211A, it is found that the merger of NS and WD cannot be ruled out as the origin of GRB 211211A to produce the possible kilonova emission if the remnant of the WD-NS merger is a supramassive or stable magnetar. Otherwise, it is difficult to explain the early possible kilonova emission following GRB 211211A by invoking the merger of a WD and an NS.

The Nuclear Stellar Disk (NSD), together with the Nuclear Stellar Cluster and the supermassive black hole Sgr A*, forms the central region of the Milky Way. Galactic X-ray background emission is known to be associated with the old stellar population, predominantly produced by accreting white dwarfs. In this work we characterize the X-ray emission of the Galactic Center (GC) region using wide-field observations with the ART-XC telescope on board the SRG observatory in the 4-12 keV energy band. Our analysis demonstrates that the X-ray emission of the GC at a spatial scale of a few hundred parsecs is dominated by the regularly shaped NSD aligned in the Galactic plane, and characterized by latitudinal and longitudinal scale heights of approximately 20 pc and approximately 100 pc, respectively. The measured flux, 6.8 (+0.1, -0.3) x 10^-10 erg/s/cm^2 in the 4-12 keV band, corresponds to a luminosity of 5.9 (+0.1, -0.3) x 10^36 erg/s, assuming the GC distance of 8.178 kpc. The average mass-normalized X-ray emissivity of the NSD, 5.6 (+0.5, -0.7) x 10^27 erg/s/M_sun, exceeds the corresponding value of the Galactic ridge by a factor of 3.3 (+0.4, -0.5), confirming other studies. We also perform a deprojection of the observed NSD surface brightness distribution in order to construct a three-dimensional X-ray luminosity density model, which can be directly compared to the existing three-dimensional stellar mass models. The emissivity of the NSD as a function of distance to Sgr A* reveals a centrally concentrated maximum, indicating an enhanced radiative output per unit stellar mass in the inner NSD region. Finally, we conclude that the observed spatial properties of the NSD are fully consistent with the stellar mass density distribution, leaving only a small room for a possible diffuse component.

In this manuscript, we propose, for the first time, an oversimplified but potentially effective gravitational acceleration model to interpret the double-peaked narrow emission lines (DPNELs) shifted in the same direction. We adopt the framework of a merging kpc-scale dual-core system in an elliptical orbit, which has an emission-line galaxy with clear narrow line regions (NLRs) merging with a companion galaxy lacking emission line features. Due to gravitational forces induced by both galaxies on the NLRs, the accelerations of the far-side and near-side NLR components may share the same vector direction when projected along the line-of-sight, leading the velocities of the observed DPNELs to shift in the same direction. Our simulations indicate that the probability of producing double-peaked features shifted in the same direction reaches 5.81% in merging kpc-scale dual core systems containing emission-line galaxies. Besides the expected results from our proposed model, we identify a unique galaxy SDSS J001050.52-103246.6, whose apparent DPNELs shifted in the same direction can be plausibly explained by the gravitational acceleration model. This proposed model provides a new path to explain DPNELs shifted in the same direction in the scenario that the two galaxies align along the line-of-sight in kpc-scale dual-core systems.

Radio-detection is now an established technique for studying ultra-high-energy (UHE) cosmic rays with energies exceeding $\sim 10^{17}$ eV. The next generation of radio experiments, such as the Giant Radio Array for Neutrino Detection (GRAND), aims to expand this technique to the observation of Earth-skimming UHE neutrinos, which requires the detection of very inclined extensive air showers (EAS). Currently, GRAND is validating its detection principle -- autonomous radio detection -- in particular through the prototype array GRANDProto300, deployed in the Gobi Desert. In this phase, the array is limited to detecting inclined EAS from cosmic rays. Neutrinos cannot be observed because of the restricted detector size. We present a method to reconstruct the arrival direction and energy of EAS with zenith angles above $60^\circ$, applicable as well to upward-going trajectories. The approach combines a point-source-like description of the radio wavefront with the so-called Angular Distribution Function (ADF), a phenomenological model describing the angular pattern of radio signal amplitudes in the 50--200 MHz band. Applied directly to the voltage traces, this method enables efficient event selection with accurate direction reconstruction and a first-order energy estimate. We validate the approach with both simulations and experimental data, and reconstruct the first cosmic-ray candidates detected by GRANDProto300.

The NUTRIG project is dedicated to the development of advanced radio self-trigger methods for large-scale arrays such as the Giant Radio Array for Neutrino Detection (GRAND). The developed techniques are based on features of the radio emission of air showers to perform an efficient online rejection of background. We first describe a first-level trigger (FLT) template-matching method that uses the shape of transient radio pulses measured at the detection-unit level to target those induced by air showers. We present trigger efficiencies and throughput tests of the template-matching FLT in controlled laboratory conditions. Next, we describe the second-level trigger (SLT), which utilizes the measured FLT times and corresponding voltage amplitudes to construct a trigger at the array level. We present offline performances of the SLT, which performs a coarse reconstruction of timing, direction, and polarization of the air shower.

Interband lags in the optical continua of active galactic nuclei (AGN) have been observed over years of monitoring, yet their physical origins remain unclear. While variable interband lags have been found in a few individual AGN potentially, the temporal behavior of interband lags of an AGN sample has not been explored systematically. Here, we analyze the interband lags of 94 bright AGN at $z<0.8$, using both seasonal one-year and full six-year $gri$-band light curves from Zwicky Transient Facility Data Release 22. We find that more than half of 94 AGN show significant seasonal variations in the interband lags. Besides, the short-term lags, derived by averaging lags inferred from multiple seasonal light curves, are consistently smaller than the long-term lags, which are inferred from the full six-year light curves. This supports recent theoretical simulations where the lag measurement is sensitive to the baseline of light curve and the lag variation could be simply attributed to the inherent randomness of AGN variability. Our findings suggest that the interband lags of AGN are more complex and stochastic than commonly thought, and highlight the importance of high-precision time-domain surveys in uncovering the properties of AGN variability as well as the associated accretion physics.

Measuring ultra-high energy neutrinos, with energies above $10^{16}$ eV, is the next frontier of the emerging multi-messenger era. Their detection requires building a large-scale detector with 10 times the instantaneous sensitivity of current instruments, sub-degree angular resolution, and wide daily field of view. The Hybrid Elevated Radio Observatory for Neutrinos (HERON) is designed to be that discovery instrument. HERON combines the complementary features of two radio techniques being demonstrated by the BEACON and GRAND prototypes. Its preliminary design consists of 24 compact, elevated phased stations with 24 antennas each, embedded in a sparse array of 360 standalone antennas. This setup tunes the energy threshold to below 100 PeV, where the neutrino flux should be high. The sensitivity of the phased stations combines with the powerful reconstruction capacities of the standalone antennas to produce an optimal detector. HERON is planned to be installed at an elevation of 1,000 m across a 72 km-long mountain range overlooking a valley in Argentina's San Juan province. It would be connected to the worldwide network of multimessenger observatories and search for neutrino bursts from candidate sources of cosmic rays, like gamma-ray bursts and other powerful transients. With HERON's deep sensitivity, this strategy targets discoveries that cast new light into the inner workings of the most violent astrophysical sources at uncharted energies. We present the preliminary design, performances, and observation strategy of HERON.

While decades of observations in the TeV gamma-ray band have revealed more than 200~sources with radio or X-ray counterparts, there remain dozens of unidentified TeV sources, which may provide crucial information of cosmic ray (CR) accelerators. HESS J1626$-$490 is an unidentified TeV gamma-ray source but is expected to originate from CRs that escaped from the nearby supernova remnant (SNR) G335.2+0.1 and are interacting with dense interstellar clouds. To test this scenario, we scrutinize the properties of the SNR and search for non-thermal counterparts by analyzing observational data in the radio, X-ray, and GeV gamma-ray bands. From analysis of the H\,{\sc i} and $^{12}$CO ($J{=}1{-}0$) line data, we identify the cloud associated with the SNR and compare the morphologies of the cloud and the gamma-ray emission. The distance and age of the SNR are estimated to be $3.3 \pm 0.6$~kpc and ${\sim}5$~kyr, respectively. From X-ray and GeV gamma-ray data analyses, we find an extended GeV gamma-ray emission overlapping with the SNR and H.E.S.S. source regions but no X-ray emission. The location of the peak of the extended GeV emission changes from near the SNR at $\lesssim 1$~GeV to the H.E.S.S. source at $>10$~GeV. We find a spectral hardening at ${\sim}50$~GeV, which is consistent with the existence of two components in the gamma-ray emission. We find that a combination of emission from the escaped CRs and the SNR itself can reproduce the observed broadband spectrum, on the assumption that the SNR has accelerated protons to ${\gtrsim}100$~TeV in the past.

Studying the habitability, internal structure and composition of exoplanets is crucial for understanding their potential to sustain life beyond our solar system. Characterizing planetary structures and atmospheric evolution provides valuable insights into surface conditions and the long-term habitability of these planets. In this study, we present a comprehensive analysis of exoplanets spanning from super-Earths to mini-Neptunes ($R_{\textrm{p}}$ $\leq$ 4 $R_{\oplus}$ and $M_{\textrm{p}}$ $\leq$ 15 $M_{\oplus}$) located within the extended habitable zone, along with parameterization of their host stars. We find that the planets in our sample orbit M dwarf stars and are tidally locked to them. Using archival photometric data from Gaia, Pan-STARRS1, 2MASS, and WISE, we estimate the atmospheric and physical parameters of the host stars. We also model the interior structure of these planets to infer their possible compositions. Additionally, under the assumption that these exoplanets can accrete a gaseous layer, we model the envelope fraction of the habitable exoplanets. With an Earth-like rocky composition, LHS 1140 b and TOI-1452 b can hold onto negligible amount of their initial gas layer. However, sustaining a sufficient amount of atmosphere over time, the planets LP 791-18 c, LTT 3780 c and K2-18 b are likely to be water worlds. The models suggest a water rich composition for TOI-1266 c without any significant amount of atmosphere. Modeling interior compositions and atmospheric escape scenarios allow us to assess the potential habitability of these planets by evaluating the likelihood of surface liquid water and the retention of stable atmospheres.

The next generation of cosmology surveys will probe the matter distribution of the universe to unparalleled precision. To match this level of precision in cosmological parameter estimation, we need to use information at small scales of $\sim$ 1 Mpc, which requires an accurate model of baryonic feedback. In this paper, we employ the Dark Matter + Baryon (DMB) model, a flexible halo model that is well-fit to various hydrodynamical simulations, to extract information on baryonic feedback from galaxy cluster observables. Using a sample of thermal Sunyaev-Zeldovich (tSZ) selected galaxy clusters from the Atacama Cosmology Telescope (ACT) - with masses calibrated via weak lensing from the Dark Energy Survey (DES) - we develop a robust end-to-end pipeline that directly models the calibrated observables. Our analysis demonstrates that the tSZ Y-M relation can constrain several DMB model parameters, providing key insights into baryonic feedback effects on cosmic shear at the several percent level. We find a preference for intermediate to strong levels of feedback, which is both consistent with several hydrodynamic simulations and competitive with similar analyses performed on complementary probes. Finally, we discuss the implications of our results in the context of current and upcoming cosmic shear surveys.

Single pulse studies offer vital insights into the emission physics of pulsars, particularly in the case of young, nearby sources where intrinsic variability is often pronounced. PSR~B0656+14, known for its sporadic and sometimes intense pulses, provides an excellent opportunity to investigate such behaviour at low radio frequencies. This study aims to characterize the single pulse behaviour of PSR~B0656+14 using low-frequency observations at 110-190 MHz from the Irish LOFAR station. Single-pulse extraction is performed, and individual pulse DMs are estimated to probe pulse-to-pulse dispersion variability. We also perform a wait-time analysis to understand the statistical nature of pulse occurrence, and estimate the spectral index from frequency-resolved flux density measurements. The pulse energy distribution is modelled using a combination of log-normal and power-law components. A total of 41 pulses were detected in a 5-hour observation, allowing a wait-time distribution analysis which is well-modelled by an exponential function, indicative of a Poisson process. Profile stability analysis indicates that a significant number of pulses are required to reach a stable average profile, unusual compared to many other pulsars. The single-pulse spectral index varies significantly from pulse to pulse, with a mean value of $\alpha = -0.5$ and a standard deviation of $\Delta\alpha=1.3$. The pulse energy distribution shows a hybrid behaviour, consistent of a log-normal distribution and a power-law tail. Our results confirm that PSR~B0656+14 exhibits highly variable, memory-less emission at low frequencies, with characteristics that resemble those seen in some rotating radio transients (RRATs). If such variability proves to be widespread among pulsars, population synthesis models and survey yield predictions would need to incorporate this currently overlooked feature to ensure accuracy.

Superfluid vortices in neutron star crusts are thought to be pinned to the lattice of nuclei in the crust. The unpinning of superfluid vortices in spin glitches therefore motivates us to study the vortex-crust interaction explicitly with molecular dynamics. In this work, we present a new molecular dynamics methods to characterize the response of the crust to a rigid vortex. When vortex pinning forces and nearest neighbor Coulomb forces are comparable, we observe a qualitatively new phenomena of lattice entrainment with implications for the elastic evolution of the crust.

The inner regions of the Milky Way are known to contain an enigmatic population of prominent molecular clouds characterised by extremely broad lines. The physical origin of these ''extended velocity features'' (EVFs) is still debated, although a connection with the ''dust lanes'' of the Galactic bar has been hypothesised. In this paper, we search for analogous features in the dust lanes of nearby barred galaxies using the PHANGS-ALMA CO(2-1) survey. We aim to confirm existence of EVFs in other galaxies and to take advantage of the external perspective to gain insight into their origin. We study a sample of 29 barred galaxies and find that 34% contain one or more EVFs, while the remaining lack obvious signs of EVFs. Upon analysing the physical properties of the EVFs, we find they possess large virial parameters, ranging from few hundreds to several thousand, indicating that they are strongly out-of-equilibrium. The most likely explanation for their origin is extreme cloud-cloud collisions with relative velocities in excess of 100km/s in highly non-circular flow driven by the bar. This interpretation is consistent with previous high-resolution observations in Milky Way. Further corroboration of this interpretation comes from the inspection of high-sensitivity infrared observations from the PHANGS-JWST Treasury Survey that reveals streams of gas that appear to be hitting the dust lanes at locations where EVFs are found. We argue that EVFs are the clearest examples of cloud-cloud collisions available in literature and represent a unique opportunity to study cloud collisions and their impact on star formation.

Identification of primary cosmic rays on an event-by-event basis is a much-desired capability of cosmic-ray observatories. Several cosmic-ray air-shower experiments use so-called photon tags for gamma hadron primary particle discrimination. These photon tag variables are derived from the total signals measured by an array of detectors and are correlated with the total number of muons in the air shower. In this work, variables based on time distribution of signals in detectors (trace-based discriminant variables) are studied and compared to total-signal-based variables. This study relies on simulated high-energy cosmic-ray air showers with energies around 10^17.5eV. Since the variables discussed are derived from total signals and their time traces, which can be directly measured in real data, they are suitable for use as discriminant variables in the real ground-based cosmic ray experiments.

The origin of ultra-high-energy cosmic rays (UHECRs) is one of the most intriguing mysteries in astroparticle physics and high-energy physics. Since UHECRs with light mass compositions are less deflected by the Galactic and extragalactic magnetic fields, their arrival directions are more strongly correlated with their origins. Charged-particle astronomy with UHECRs is hence a potentially viable probe of extremely energetic phenomena in the universe. The Global Cosmic Ray Observatory (GCOS) is a proposed next-generation observatory to elucidate these origins through precise measurements of UHECRs with unprecedented exposure and mass identification capabilities. We will focus on the ideas and requirements for GCOS summarized in arXiv:2502.05657 and share the recent advances in detector developments and future perspectives with interdisciplinary research.

Two fast radio bursts, FRB 20121102A and FRB 20180916B, show periodic activity with cycles of 159.3 and 16.33 days, respectively. These cycles consist of active and inactive windows, with peak activity centered within the active phase. For FRB 20180916B, studies have reported a frequency-dependent-or ``chromatic''-behaviour, where the activity window starts earlier and becomes narrower at higher frequencies. The activity across frequencies is typically modeled with a power law. In this work, we develop a simple model that combines the phase and frequency dependence of the activity windows of FRB 20121102A and FRB 20180916B. Our goal is to perform a chromaticity study for FRB 20121102A, incorporating model improvements to account for the cyclic nature of its activity window, and to compare the chromatic behaviour between both periodic FRBs. We standardise the detections from the 425 observing epochs for FRB 20121102A and the 214 epochs for FRB 20180916B to account for differences in radio telescope sensitivity. To the normalised detection rate phase distribution, we fit a von Mises distribution and extract the peak activity phase and activity width. These quantities, as a function of frequency, are then modelled as power-laws to construct the chromatic model. For both sources, the activity window starts earlier at higher frequencies. However, FRB 20121102A shows an activity window broadening at higher frequencies, while FRB 20180916B broadens at lower frequencies. FRB 20121102A appears to remain active during at least 80% of the cycle at C-band. The observed chromatic behaviour of FRB 20180916B is consistent with previous findings. For FRB 20121102A, a chromaticity in its activity window is also seen; however, the source appears to be active for longer at higher frequencies, opposite to the behaviour of FRB 20180916B.

Constraining neutrino mass through cosmological observations relies on precise simulations to calibrate their effects on large scale structure, while these simulations must overcome computational challenges like dealing with large velocity dispersions and small intrinsic neutrino perturbations. We present an efficient N-body implementation with semi-linear neutrino mass response which gives accurate power spectra and halo statistics. We explore the necessity of correcting the expansion history caused by massive neutrinos and the transition between relativistic and non-relativistic components. The above method of including neutrino masses is built into the memory-, scalability-, and precision-optimized parallel N-body simulation code CUBE 2.0. Through a suite of neutrino simulations, we precisely quantify the neutrino mass effects on the nonlinear matter power spectra and halo statistics.

Recent observations reveal small-scale reconnection-driven plasma ejections, often termed nanojets, triggered by magnetic field interactions at slight misalignment angles. These fast, collimated plasma ejections are $\sim$1.5 Mm long and $\sim$0.5 Mm wide. In this study, we analyze two high-resolution extreme ultraviolet imaging datasets from the Extreme Ultraviolet Imager onboard the Solar Orbiter mission, corresponding to an eruptive (M7.6) and a confined (C1.2) flare, to investigate the dynamics of nanoflare ejections and, for the first time, compare their properties in distinct magnetic environments. We identified 59 nanoflare ejections: 44 in the eruptive flare and 15 in the confined flare event. Our analysis reveals that these events form two distinct classes: confined events exhibit lower speeds (41--174 kms$^{-1}$) and lower kinetic energies ($10^{20}$--$10^{22}$ erg), placing them closely in or near the picoflare energy regime, while eruptive events show higher speeds (131--775 kms$^{-1}$) and higher kinetic energies ($10^{22}$--$10^{24}$ erg), falling within the nanoflare regime. Furthermore, magnetic field extrapolations reveal a highly sheared arcade with greater twist and higher magnetic energy density in the eruptive event, compared to the less twisted configuration in the confined event. We infer that this sheared arcade configuration in the eruptive event creates favorable conditions for higher speeds and kinetic energies, unlike the less braided structure in the confined event. Our findings highlight the crucial role of the surrounding magnetic environment in regulating the energetics of nanoflare ejections in the solar atmosphere.

Very Long Baseline Interferometry (VLBI) measures two standard observables: group delay and fringe frequency (delay rate). While group delay is widely used to estimate a broad range of geodetic and astrometric parameters, fringe frequency has, to date, been largely ignored. Here, we show that the fringe frequency is a unique tool for direct estimation of the instantaneous Earth angular rotation velocity, which is not accessible with the group delay alone. We estimate the magnitude of the Earth angular rotation velocity using a 30-year set of VLBI data and obtain daily estimates of X and Y angles linking the Instantaneous Rotation Pole (IRP) and the International Celestial Reference System (ICRS) pole. The plain least-squares method was applied to analyse the fringe frequency available from routine geodetic VLBI observations. We estimate three components of the Earth rotation vector on a daily basis with a formal error of 1 prad/s in relative units, or better, if a large international VLBI network is at work. The newly obtained values can be used to monitor the Earth rotation irregularity in parallel to the traditional length-of-day (LOD) values and to directly assess the modern precession-nutation theory.

The B-mode of polarization of the CMB is a uniquely powerful probe of gravitational waves produced in the very early Universe. But searches for primordial B-mode anisotropies must contend with gravitational lensing, which induces late-time B-modes not associated with gravitational waves. These lensing B-modes can be removed -- i.e., delensed -- using observations of the E-modes and a proxy of the matter fluctuations along the line of sight that caused the deflections. The number density and redshift reach of galaxy surveys such as the upcoming Rubin observatory offer attractive prospects for using them to delens B-mode data from CMB experiments such as the Simons Observatory, LiteBIRD or CMB-S4. However, stochasticity and non-linear galaxy bias may in principle decorrelate the galaxy field from the matter distribution responsible for the lensing effect, thus hindering efforts to delens B-modes. In addition, non-linear gravitational evolution and bias introduce non-Gaussianities in the large-scale structure which further complicate the modelling. We quantify these effects by populating an N-body simulation with a magnitude-limited, photometric sample of galaxies similar to Rubin's gold selection, and using them to delens CMB maps lensed by the same matter distribution. We find that pipelines that treat the galaxy overdensity as a Gaussian field will incur negligible bias on the inferred tensor-to-scalar ratio, r. Moreover, we show that even in a highly conservative scenario where only the linear bias of the galaxies can be determined, the bias on r arising from this simplification is well within statistical uncertainties for a cosmic-variance limited scenario where Rubin-like galaxies are used for delensing.

We present the results of a deep neutral hydrogen (Hi) observation of the early-type galaxy NGC 2768 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Leveraging the high sensitivity of FAST, we discover an extended gas envelope around NGC 2768. The total Hi mass is measured to be 8.1 x 10^8 M_sun , representing a magnitude increase compared to previous Westerbork Synthesis Radio Telescope (WSRT) studies. Position-velocity (PV) diagram indicates the envelope mainly involves two components: an Hi disk of NGC 2768 and a newly discovered satellite galaxy without detectable counterparts in currently deep optical surveys. The center of the gas disk is mis-aligned with the optical disk of NGC 2768, with more gas redshifted, indicating it has been disturbed. Our study indicates NGC 2768 is currently undergoing a transition from a spiral galaxy to an S0. Previous deep WSRT observations reveal two dense clumps (named as Clumps A and Clump B throughout this paper) in the center of the envelope. We find Clump A corresponds to the densest part of the disk, while Clump B might be a newly discovered satellite galaxy which probably collided NGC 2768 about 0.38 Gyr ago. We also find tidal interactions between Clump B and PGC 2599651, NGC 2768 and UGC 4808. Based on these new findings, we finally analyze hierarchical accretion history of NGC 2768.

The corona is an integral component of active galactic nuclei (AGNs) which can produce the X-ray emission. However, many of its physical properties and the mechanisms powering this emission remain a mystery. In this work, we study the coronal X-ray variabilities of 13 AGNs by Gaussian Process. 2-10 keV light curves of 13 AGNs can be successfully described by the damped-random walk (DRW) model. The extracted coronal X-ray timescales range from 3 to 50 days. In the plot of variability timescale versus black hole mass, the coronal X-ray timescales of four sources occupy almost the same region as the optical timescales of the accretion disk, with the latter matching the predicted thermal instability timescale of the disk. In contrast, the X-ray timescales of the remaining sources exhibit a systematic offset toward lower values. We propose that the coronal X-ray variability may be driven by internal processes within the corona itself (such as thermal conduction). On the other hand, it may also be triggered by local thermal instabilities occurring in different regions (close to the central black hole) of the accretion disk, which propagate to the corona via disk-corona coupling.

We investigate whether the properties of turbulent gas motions recently measured via X-ray spectroscopy in the Coma cluster of galaxies by XRISM are in tension with the "classical" fluid picture of the intracluster medium on large scales, as produced by a typical high-resolution cosmological simulation. We use a high-resolution simulation of Coma-like cluster of galaxies and show that the Kolmogorov-like turbulence measured in the simulation yields to velocity structure functions and line-width that fully compatible with those measured by the XRISM observation of Coma. These results highlight the combined biases driven by the inhomogeneity of turbulence in the intracluster medium and by the X-ray weighting in observations, and appear to release the tension between the XRISM data and current numerical simulations, showing that a turbulent spectrum much steeper than Kolmogorov is not required by current observational data.

The design of transfers to periodic orbits in the Earth--Moon system has regained prominence with NASA's Artemis and CNSA's Chang'e programs. This work addresses the problem of linking ballistic capture trajectories - exploiting multi-body dynamics for temporary lunar orbit insertion - with bounded periodic motion described in the circular restricted three-body problem (CR3BP). A unified framework is developed for optimizing bi-impulsive transfers to families of periodic orbits via a high-order polynomial expansion of the CR3BP dynamics. That same expansion underlies a continuous 'abacus' parameterization of orbit families, enabling rapid targeting and analytic sensitivity. Transfers to planar periodic-orbit families (Lyapunov L1 and L2, and distant retrograde orbits) are addressed first, followed by extension to spatial families, such as butterfly and halo L1/L2 orbits, with an emphasis towards Near-Rectilinear Halo Orbits (NRHOs). Numerical results demonstrate low-cost solutions and validate the method's adaptability for the design of lunar missions. The optimized trajectories can inform an established low-energy transfer database, enriching it with detailed cost profiles that reflect both transfer feasibility and underlying dynamical relationships to specific periodic-orbit families. Finally, the proposed transfers provide reliable initial guesses for rapid refinement, readily adaptable for further optimization across mission-specific needs.

The spectral energy distributions (SEDs) of 20 Milky Way (MW), 9 Large Magellanic Cloud (LMC), 7 Small Magellanic Cloud (SMC), 12 M31, and 7 M33 (classical) Cepheids with periods longer than 50 days were constructed using photometric data from the literature and fitted with model atmospheres with the aim of identifying objects with an infrared excess. The SEDs were fitted with stellar photosphere models to derive the best-fitting luminosity and effective temperature; a dust component was added when required. The distance and reddening values were taken from the literature. WISE and IRAC images were inspected to verify whether potential excess emission was related to the central objects. Only one star with a significant infrared (IR) excess was found in the LMC and none in the SMC, M31, and M33, contrary to earlier work on the MW suggesting that IR excess may be more prominent in MW Cepheids than in the Magellanic Clouds. One additional object in the MW was found to have an IR excess, but it is unclear whether it is a classical Cepheid or a type-{\sc ii} Cepheid. The stars were plotted in a Hertzsprung--Russell diagram (HRD) and compared to evolutionary tracks for CCs and to theoretical instability strips. For the large majority of stars, the position in the HRD is consistent with the instability strip. For stars in the MW uncertainties in the distance and reddening can significantly change their position in the HRD.

We evaluate global asymmetry in the luminosities of neutrinos emitted from rapidly-rotating proto-neutron stars (PNS's). We build axisymmetric models of PNS's in mechanical equilibrium with rotation by adding prescribed angular momentum distributions by hand to non-rotational PNS models, which are extracted from a one-dimensional (spherically symmetric) PNS cooling calculation at different times: \(t=2, 6, 10, 20, 30\)s after a supernova explosion. We then conduct two-dimensional (spatially axisymmetric) neutrino transport calculations on top of them with the matter profiles (and the spacetime geometry) fixed. We find for the rapidly-rotating models with \(T/|W|\sim 5\times 10^{-2}\) that the neutrino luminosity changes by \(\sim 3 \% \) depending on the observer position. We give detailed analyses of the neutrino-hemispheres as well as the neutrino luminosities that are defined observer-wise. We also calculate the low-frequency (\(\lesssim 1{\rm Hz}\)) gravitational waves produced by the neutrinos radiated asymmetrically. We find that those gravitational waves, if emitted from the Galactic center, can be detected by planned detectors such as B-DECIGO, DECIGO and AILA. Finally, we look for crossings in the energy-integrated angular distributions in momentum space for the electron neutrino sector, a signature of the fast flavor conversion. We find them near the PNS surface in all models.

Pulsars are typically characterized by their stable, highly magnetized, and fast-rotating nature, which underpins their persistent radio emissions. However, the discovery of prolonged radio-quiet ("off") states in intermittent pulsars, such as PSR B1931+24, has challenged the most fundamental theory of pulsar magnetospheric emission. Despite long-term monitoring with several telescopes, continuous emission during these "off" states had not been detected in 20 years of observations. Fortunately, sensitive observations via Five-hundred-meter Aperture Spherical radio Telescope (FAST) revealed the mysterious weak emission containing occasional bursting dwarf pulses during the "off" states of PSR B1931+24. Along with a substantial decrease in flux density, a significant contraction in the integrated pulse width is measured in the "off" state compared to the radio-loud ("on") state, indicating alterations in the plasma supply and magnetospheric structure. Additionally, a previously unobserved dyssynchronous, nonuniform emission pattern is found in both states, supporting theories of a spatially inhomogeneous pair-cascade mechanism and challenging models of spatially coherent discharge. Furthermore, occasional dwarf pulses detected during the "off" state show flux and width distributions similar to those of the "on" state pulses, suggesting a potential link between the "on" and "off" state emissions of PSR B1931+24. Consequently, dwarf pulses are unlikely to represent a distinct emission mode as previously thought; instead, they appear to be part of a continuum within the pulsar's emission behaviors observed during the "on" state. These findings strongly support the basic theory of the pulsar magnetospheric emission and significantly advance our understanding of pulsar magnetospheric dynamics and their emission mechanisms.

We determined the ortho/para (o/p) ratios of NH2D and NHD2 in the archetypical pre-stellar core L1544. The core was observed in the two lowest rotational lines of ortho- and para-NH2D using the APEX and the IRAM 30 m telescopes. The ground-state lines of ortho- and para-NHD2 were observed with APEX. The distributions of chemical abundances in the core were predicted using a gas-grain chemistry model with two different scenarios concerning proton transfer reactions in the gas. One of the scenarios, the so-called full scrambling (FS), allows protons and deuterons to be completely mixed in the intermediate reaction complex before dissociation, whereas the other describes these reactions as proton or deuteron hops (PH). We also tested assumed abundance profiles independent of the chemistry models. Radiative transfer calculations were used to simulate the observed NH2D and NHD2 lines from the predicted and assumed abundance profiles. Our modelling efforts suggest that the ground-state lines of NH2D and NHD2 at the wavelength 0.9 mm that are observable with the same beam and in the same spectrometer band are the most reliable probes of the o/p ratios. Simulations using the PH reaction scheme show systematically better agreement with the observations than simulations with the FS model. Simulations using a broken power law abundance profile as a function of the gas density give spin ratios that are close to the predictions of the PH scenario: o/p-NH2D=2.85+-0.05, o/p-NHD2=2.10+-0.06 (1 sigma). The o/p ratios predicted by the PH scenario in the gas phase correspond to the nuclear spin statistical weights, that is, o/p-NH2D=3, o/p-NHD2=2. In view of the fact that H and D atom addition reactions on grain surfaces also result in these ratios, it is reasonable to assume that the spin ratios of interstellar ammonia and its deuterated forms are in general equal to their statistical values.

We report unnoticed but intriguing features in the peculiar nuclear transient AT 2022fpx, and investigate its type. These features include the constantly red optical color of $g-r>0$, stable soft X-ray flare ($kT\sim100$ eV) in the past $\sim$550 days, prominent mid-infrared echo peaked at $\sim$$10^{43.3}$ erg s$^{-1}$ and the confirmation of a weak AGN by weak flares in pre-event WISE mid-infrared light curves with no contemporary optical, radio or X-ray counterparts. The combination of the optical red color and possible TDE origin of AT 2022fpx is particularly attractive, as it challenges the most widely accepted and adopted "blue color" criterion for optical TDE selection. Although we still cannot confirm whether the red color is intrinsic, we however find that the "blue color" criterion can filter out normal TDEs whose optical-UV spectral energy distributions (SEDs) are either severely contaminated by prominent emission lines (especially H$\alpha$) or heavily dust-reddened. Hence, its potential selection effect may have been imprinted on the whole optical TDE family. Blackbody fitting on the optical (rest-frame $\sim$$4000-7000$ Å) and optical-UV ($\sim$$2000-7000$ Å) SEDs of four TDEs with high-cadence UV observations shows that $T_\mathrm{bb}$ rise by $\sim$40$-$110 % when the UV bands are included. The power-law models ($f_{\lambda}\propto\lambda^{-\alpha}$ with $\alpha=2-3$) can fit the rest-frame $\sim$$2000-7000$ Å SEDs more consistently, indicating that SEDs should peak at shorter wavelengths, but not simple blackbodies. Hence, the estimated released energy for the optical-UV bright but X-ray faint TDEs based on blackbody SED fitting should be significantly lower than the intrinsic energy.

The work is devoted to the study of periodic structures in the magnetic field of the photosphere. The polarity of the Sun's magnetic fields shows cyclicity with periods of 1-3 years, which is possibly due to quasi-biennial variations that are found in various parameters of solar activity, interplanetary space, and in cosmic rays. In this work, variations in weak magnetic fields of the photosphere were investigated, which required the use of special data processing methods. Synoptic maps of the photospheric magnetic field for the period 1978-2016 (NSO Kitt Peak) were used. To highlight the contribution of weak magnetic fields, the saturation threshold for synoptic maps was set at 5 G. Based on thus transformed synoptic maps, a time-latitude diagram was constructed. For the selected latitude intervals, smoothing, trend removal and Fast Fourier Transform (FFT) were performed. In the time-latitude diagram, 6 spatiotemporal regions with distinct cyclic structure were observed during 4 solar cycles (21-24). The periods ranged from 0.5 to 4 years with a maximum at 1.2 years. The amplitude of ripples is significantly higher for those intervals in which the polar field had positive sign. This effect confirms the connection of ripples variations with the polarity of the 22-year magnetic cycle.

Polarization measurements provide strong constraints on magnetic fields in star-forming systems. While magnetic field estimates of a few kiloGauss (kG) have been obtained near the surface of low-mass protostars, there are no analogous measurements in the immediate vicinity of the surface of massive protostars. We report the measurement of radio continuum circular polarization (CP) towards a massive protostar IRAS 18162-2048 for the first time wielding Karl G. Jansky Very Large Array (VLA) observations. The fractional CP varies between $3-5\%$ across the observed frequency range of $4-6$ GHz. We consider multiple hypotheses for the production of CP and propose (i) gyrosynchrotron emission and (ii) Faraday conversion due to turbulence in the magnetic medium - both driven by mildly relativistic electrons as plausible mechanisms. We estimate, for the first time, a magnetic field $B\gtrsim20-35$ G close to the massive protostar. The Lorentz factor of the low energy electrons is estimated to be in the range $\gamma_{min}\sim5-7$ for gyrosynchrotron emission and $80-100$ for Faraday conversion from our observations. The magnetic field estimate can provide important constraints to the formation models of massive stars.

LiteBIRD, the Lite (Light) satellite for the study of $B$-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission focused on primordial cosmology and fundamental physics. In this paper, we present the LiteBIRD Simulation Framework (LBS), a Python package designed for the implementation of pipelines that model the outputs of the data acquisition process from the three instruments on the LiteBIRD spacecraft: LFT (Low-Frequency Telescope), MFT (Mid-Frequency Telescope), and HFT (High-Frequency Telescope). LBS provides several modules to simulate the scanning strategy of the telescopes, the measurement of realistic polarized radiation coming from the sky (including the Cosmic Microwave Background itself, the Solar and Kinematic dipole, and the diffuse foregrounds emitted by the Galaxy), the generation of instrumental noise and the effect of systematic errors, like pointing wobbling, non-idealities in the Half-Wave Plate, et cetera. Additionally, we present the implementation of a simple but complete pipeline that showcases the main features of LBS. We also discuss how we ensured that LBS lets people develop pipelines whose results are accurate and reproducible. A full end-to-end pipeline has been developed using LBS to characterize the scientific performance of the LiteBIRD experiment. This pipeline and the results of the first simulation run are presented in Puglisi et al. (2025).

The James Webb Space Telescope (JWST) now allows us to observe galaxies at the end of cosmic dawn ($z \sim 10-15$) with unprecedented detail, revealing their morphologies, sizes, and internal structures. These observations offer crucial insights into the physical processes driving early galaxy formation. In this work, we introduce the Phoebos hydrodynamical cosmological simulation, a state-of-the-art 100 Mpc volume designed to study the formation and evolution of galaxies at the end of cosmic dawn and into the epoch of reionization. Phoebos includes a stellar feedback model that is intentionally weak, in order to address the high abundance of massive galaxies seen by JWST at early epochs. At variance with most large cosmological hydrodynamical simulations, we do not employ an effective equation of state model, instead our radiative cooling model allows us to capture the multi-phase nature of the gas inside and around galaxies. Phoebos reproduces key observables of early galaxy formation at $z \gtrsim 8$, including the stellar mass function and the stellar-to-halo mass relation. It also recovers the observed slope of the stellar size-to-mass relation and matches the specific star formation rate remarkably well. These results suggest that highly efficient star formation in the presence of only mild regulation from stellar feedback, drives early galaxy growth, supporting a scenario of rapid stellar mass assembly during cosmic dawn. There are indications in the cosmic star formation density that, at lower redshifts, Phoebos might overpredict the stellar mass within the systems, suggesting that a transition to a stronger stellar feedback may be necessary to reproduce later-time observations. These results highlight the potential of Phoebos to interpret JWST observations and to probe the evolving physical processes that shape galaxy formation.

Single massive satellites are of great observational interest, as they can produce prominent and potentially detectable signatures. For terrestrial planets and super-Earths, giant impacts in the late stages of formation may generate dense self-gravitating disks - favourable environments for the formation of such satellites. Motivated by this, we explore satellite formation in dense solid-particle disks through three-dimensional N-body simulations, focusing on the effects of disk mass and the surface density exponent. Our results reveal significant variability in the masses and configurations of satellites formed under identical disk parameters, highlighting the stochastic nature of the process. Higher disk masses and flatter surface density profiles favour the formation of more massive satellites. Disks with masses above 0.03 planetary masses typically yield a single dominant satellite, while those between 0.003 and 0.03 tend to form two-satellite systems. On average, the mass of the largest satellite scales linearly with the initial disk mass, in agreement with analytical predictions. We estimate that a disk with a minimal mass of 0.03 planetary masses around a 1.6 Earth-mass planet orbiting a Sun-like star could form an Earth-Moon-like system detectable by telescopes with a photometric precision of 10 parts per million - a level achievable by the James Webb Space Telescope.

We attempt to measure the $z = 8.2$ Epoch of Reionization (EoR) 21-cm bispectrum (BS) using Murchison Widefield Array (MWA) $154.2~\mathrm{MHz}$ data. We find that $B(k_{1\perp}, k_{2\perp}, k_{3\perp}, k_{1\parallel}, k_{2\parallel})$ the 3D cylindrical BS exhibits a foreground wedge, similar to $P(k_{1\perp},k_{1\parallel})$ the 21-cm cylindrical power spectrum. However, the BS foreground wedge, which depends on $(k_{1\perp},k_{1\parallel})$, $(k_{2\perp},k_{2\parallel})$ and $(k_{3\perp},k_{3\parallel})$ the three sides of a triangle, is more complicated. Considering various foreground avoidance scenarios, we identify the region where all three sides are outside the foreground wedge as the EoR window for the 21-cm BS. However, the EoR window is contaminated by a periodic pattern of spikes that arises from the periodic pattern of missing frequency channels in the data. We evaluate the binned 3D spherical BS for triangles of all possible sizes and shapes, and present results for $\Delta^3$ the mean cube brightness temperature fluctuations. The best $2\sigma$ upper limits we obtain for the EoR 21-cm signal are $\Delta^3_{\rm UL} = (1.81\times 10^3)^3~\mathrm{mK}^3$ at $k_1 = 0.008~\mathrm{Mpc}^{-1}$ and $\Delta^3_{\rm UL} = (2.04\times 10^3)^3~\mathrm{mK}^3$ at $k_1 = 0.012~\mathrm{Mpc}^{-1}$ for equilateral and squeezed triangles, respectively. These are foreground-dominated, and are many orders of magnitude larger than the predicted EoR 21-cm signal $(\sim 10^3 ~\mathrm{mK}^3)$.

Most direct N-body integrations of planetary systems use a symplectic integrator with a fixed timestep. A large timestep is desirable because simulations run fast. However, simulations yield unphysical results if the timestep is too large. Surprisingly, no systematic convergence study has been performed on long (Gyr) timescales. In this paper we present numerical experiments to determine the minimum timestep one has to use in long-term integrations of the Solar System in order to recover the system's fundamental secular frequencies and instability rate. We find that timesteps of up to 32 days, i.e. a third of Mercury's orbital period, yield physical results in 5 Gyr integrations. We argue that the chaotic diffusion that drives the Solar System's long-term evolution dominates over numerical diffusion and timestep resonances. Our results bolster confidence that most simulations in the literature are indeed converged and provide guidance on how to run time and energy efficient simulations while making sure results can be trusted.

The dynamically-driven double-degenerate double-detonation model has emerged as a promising progenitor candidate for Type Ia supernovae. In this scenario, the primary white dwarf ignites due to dynamical interaction with a companion white dwarf, which may also undergo a detonation. Consequently, two scenarios exist: one in which the secondary survives and another in which both white dwarfs detonate. In either case, substantial departures from spherical symmetry are imprinted on the ejecta. Here, we compute full non local thermodynamic equilibrium nebular-phase spectra in 1D and 3D to probe the innermost asymmetries. Our simulations reveal that the multidimensional structures significantly alter the overall ionisation balance, width and velocity of features, especially when the secondary detonates. In this scenario, some element distributions may produce orientation-dependent line profiles that can be centrally peaked from some viewing-angles and somewhat flat-topped from others. Comparison to observations reveals that both scenarios produce most observed features from the optical to mid-infrared. However, the current model realisations do not consistently reproduce all line shapes or relative strengths, and yield prominent optical Ar III emission which is inconsistent with the data. When the secondary detonates, including 3D effects, improves the average agreement with observations, however when compared to observations, particularly weak optical Co III emission and the presence of optical O I and near-infrared S I challenge its viability for normal Type Ia supernovae. Thus, our comparisons with normal Type Ia's tentatively favour detonation of only the primary white dwarf, but we stress that more model realisations and mid-infrared observations are needed.

The study of star cluster evolution necessitates modeling how their density profiles develop from their natal gas distribution. Observational evidence indicates that many star clusters follow a Plummer density profile. However, most studies have focused on the phase after gas ejection, neglecting the influence of gas on early dynamical evolution. We investigate the development of star clusters forming within gas clouds, particularly those with a centrally concentrated gas profile. Simulations were conducted using the \texttt{Torch} framework, integrating the \texttt{FLASH} magnetohydrodynamics code into \texttt{AMUSE}. This permits detailed modeling of star formation, stellar evolution, stellar dynamics, radiative transfer, and gas magnetohydrodynamics. We study the collapse of centrally concentrated, turbulent spheres with a total mass of $2.5\times 10^3\, M_\odot$, investigating the effects of varying numerical resolution and star formation scenarios. The free-fall time is shorter at the center than at the edges of the cloud, with a minimum value of $0.55\,\mathrm{Myr}$. The key conclusions from this study are: (1) subclusters initially form even in dense gas clouds; (2) the final stellar density profile is consistent with a Plummer profile, but more centrally concentrated than analytically predicted; (3) gas collapses globally toward the center on the central free-fall time scale, contradicting the assumption in analytical models of local fragmentation and star formation; (4) hence, gas is predominantly ejected from the innermost regions; and (5) the final value of the SFE depends on the mass and formation time of the most massive star as well as the overall star formation timescale.

The interstellar comet 2I/Borisov is the first interstellar object where compositional characterisation was possible throughout its entire perihelion passage. We report all 16 epochs of a comprehensive optical observation campaign with ESO VLT's integral field spectrograph MUSE, spanning 126 days from 2019 November 14 to 2020 March 19. The spatial dust emission of 2I/Borisov was predominantly smooth, with no seasonal effect. A jet-like feature was consistently visible. The gas production morphology of its coma was also smooth and similar for C$_2$, NH$_2$, and CN: symmetric around the photocentre. The production rates of these species gently declined into and beyond perihelion, until 2I's outburst and splitting event in early 2020 March. C$_2$, NH$_2$, and CN production rates all increased, with NH$_2$ being the most significant; the dust emission also slightly reddened. 2I/Borisov is a carbon-depleted, relatively NH$_2$-rich comet when compared to those comets yet measured in the Solar System.

We present \texttt{COMPASS}, a novel simulation-based inference framework that combines score-based diffusion models with transformer architectures to jointly perform parameter estimation and Bayesian model comparison across competing Galactic Chemical Evolution (GCE) models. \texttt{COMPASS} handles high-dimensional, incomplete, and variable-size stellar abundance datasets. % Applied to high-precision elemental abundance measurements, \texttt{COMPASS} evaluates 40 combinations of nucleosynthetic yield tables. The model strongly favours Asymptotic Giant Branch yields from NuGrid and core-collapse SN yields used in the IllustrisTNG simulation, achieving near-unity cumulative posterior probability. Using the preferred model, we infer a steep high-mass IMF slope and an elevated Supernova\,Ia normalization, consistent with prior solar neighbourhood studies but now derived from fully amortized Bayesian inference. % Our results demonstrate that modern SBI methods can robustly constrain uncertain physics in astrophysical simulators and enable principled model selection when analysing complex, simulation-based data.

Long-period transients (LPTs) are a new and enigmatic class of objects that produce bright pulsations in the radio, with periods far exceeding those seen in rotationally powered pulsars. The proposed progenitors for LPTs are contested, with white dwarfs or magnetars being likely candidates. Here, we present the discovery of ILT\,J163430+445010, a new LPT detected in a blind search for Stokes\,V transients in the LOFAR Two-Metre Sky Survey. Unusual for LPTs, J1634+44 shows pulses that are 100\% circularly polarised, as well as pulses that are 100\% linearly polarised, with the polarisation state changing from pulse to pulse. We detect 19 pulses in total, each with a total polarisation fraction of $\sim100\%$ and a pulse duration of at most 10\,s. The pulses show a periodicity at $841.24808\pm0.00015$\,s, implying a low duty cycle of $0.012$. J1634+44 has a marginally detected counterpart in the ultraviolet GALEX MIS survey and the ultraviolet/optical UNIONS survey, suggesting that it contains a white dwarf with an effective temperature between 15000\,K and 33000\,K. We do not detect J1634+44 with a deep $J$-band exposure with UKIRT at a $3\sigma$ AB magnitude limit of 24.7, ruling out a main-sequence star or ultracool dwarf with a spectral type earlier than M7. The pulses from J1634+44 follow a particular pattern, with two pulses being produced every five periods after a waiting time of two or three periods. This pattern could be a result of spin-orbit coupling in a binary system with a 5:2 or 5:3 resonance, where a companion induces beamed radio emission on the white dwarf. The companion is most likely an ultracool dwarf or another white dwarf, making J1634+44 unique among the currently known sample of LPTs.

Disc winds play a crucial role in many accreting astrophysical systems across all scales. In accreting white dwarfs (AWDs) and active galactic nuclei (AGN), radiation pressure on spectral lines is a promising wind-driving mechanism. However, the efficiency of line driving is extremely sensitive to the ionization state of the flow, making it difficult to construct a reliable physical picture of these winds. Recently, we presented the first radiation-hydrodynamic (RHD) simulations for AWDs that incorporated detailed, multi-dimensional ionization calculations via fully frequency-dependent radiative transfer, using the Sirocco code coupled to PLUTO. These simulations produced much weaker line-driven winds (Mdot_wind / Mdot_acc < 1e-5 for our adopted parameters) than earlier studies using more approximate treatments of ionization and radiative transfer (which yielded Mdot_wind / Mdot_acc ~ 1e-4). One remaining limitation of our work was the assumption of an isothermal outflow. Here, we relax this by adopting an ideal gas equation of state and explicitly solving for the multi-dimensional temperature structure of the flow. In the AWD setting, accounting for the thermal state of the wind does not change the overall conclusions drawn from the isothermal approximation. Our new simulations confirm the line-driving efficiency problem: the predicted outflows are too highly ionized, meaning they neither create optimal driving conditions nor reproduce the observed ultraviolet wind signatures. Possible solutions include wind clumping on sub-grid scales, a softer-than-expected spectral energy distribution, or additional driving mechanisms. With the physics now built into our simulations, we are well-equipped to also explore line-driven disc winds in AGN.

Dynamical activity attributed to the destruction of minor planets orbiting white dwarfs has now been photometrically monitored in individual systems for up to one decade, long enough to measure significant cessation and re-emergence of transit features. Further, periodicities which hint at the presence of debris orbiting exterior to the white dwarf Roche radius, along with widely varying estimates for debris disc lifetimes (up to Myrs), complicate theories for the formation and dynamical evolution of these systems. Here, we illustrate that minor planets orbiting white dwarfs with periods of $\approx$5-25 hours and longer while completely or partially avoiding tidal disruption satisfy the conditions for the occurrence of the sesquinary catastrophe, a phenomenon that occurs in the solar system when impacts from returning ejecta from a moon are fast enough to be erosional to the point of destruction. We hence find that the region corresponding to $\approx$1-4 white dwarf rubble-pile Roche radii represents a danger zone where the collisional timescale for the sesquinary catastrophe to occur is $\sim 10^2-10^5$ yr, suggesting that debris discs around white dwarfs are in a state of semi-continuous replenishment.

We report the discovery of CHIME J1634+44, a Long Period Radio Transient (LPT) unique for two aspects: it is the first known LPT to emit fully circularly polarized radio bursts, and it is the first LPT with a significant spin-up. Given that high circular polarization ($>90$\%) has been observed in FRB~20201124A and in some giant pulses of PSR~B1937+21, we discuss the implications of the high circular polarization of CHIME J1634+44 and conclude its emission mechanism is likely to be ``pulsar-like''. While CHIME J1634+44 has a pulse period of 841 s, its burst arrival patterns are indicative of a secondary 4206 s period, probably associated with binary activity. The timing properties suggest it has a significantly negative period derivative of $\dot{P}=-9.03(0.11)$ s s$^{-1}$. Few systems have been known to spin-up, most notably transitional millisecond pulsars and cataclysmic binaries, both of which seem unlikely progenitors for CHIME J1634+44. If the period was only associated with the spin of the object, then the spin up is likely generated by accretion of material from a companion. If, however, the radio pulse period and the orbital period are locked, as appears to be the case for two other LPTs, the spin up of CHIME J1634+44 could be driven by gravitational wave radiation.

Free-floating planetary-mass objects (FFPMOs) are known to harbor disks at young ages. Here, we present 1-13 $\mu m$ spectra for eight young FFPMOs with masses of 5-10 M$_\mathrm{Jup}$ (at ages of 1-5 Myr), using the NIRSpec and MIRI instruments on the James Webb Space Telescope. We derive fundamental properties of these targets, and find spectral types of M9.5 to L4, with effective temperatures of 1600-1900 K. The photospheric spectra of our targets show a clear diversity at similar temperatures, especially in the 3-5 $\mu m$ range, unaccounted for by existing atmospheric models. We find a silicate absorption feature in the photosphere of one of our targets, the first such detection in very young FFPMOs, indicating silicate clouds in their cool atmospheres. Six of our objects show mid-infrared excess emission above the photosphere, as well as silicate emission features, demonstrating the presence of disks. The shape and strength of the latter features constitute strong evidence of grain growth and crystallization, similar to what is seen in more massive brown dwarfs and stars. We also detect emission lines from hydrocarbon molecules in the disks of several targets. These are the lowest mass isolated objects found so far with silicate and hydrocarbon emission features arising in their disks. The presence of disks and their characteristics point to the potential for the formation of rocky companions around free-floating planetary-mass objects.

We examine scale and redshift dependence of mass-property relations (MPRs) for five hot gas properties of two large group- and cluster-scale halo samples realized by the IllustrisTNG, TNG-Cluster and FLAMINGO cosmological hydrodynamical simulations. For intrinsic properties of i) hot gas mass ($M_{\rm gas}$), ii) spectroscopic-like temperature ($T_{\rm sl}$), iii) soft-band X-ray luminosity ($L_{\rm X}$), and iv) X-ray ($Y_{\rm X}$) and v) Sunyaev-Zel'dovich ($Y_{\rm SZ}$) thermal energies, we use MPR parameters to infer mass proxy quality (MPQ) -- the implied scatter in total halo mass conditioned on a property -- for halos with $M_{\rm 500c} \geq 10^{13}{\, {\rm M}_\odot}$ at redshifts, $z \in \{0, 0.5, 1, 2\}$. We find: (1) in general, scaling relation slopes and covariance display moderate to strong dependence on halo mass, with redshift dependence secondary, (2) for halos with $M_{\rm 500c} > 10^{14}{\, {\rm M}_\odot}$, scalings of $M_{\rm gas}$ and $Y_{\rm SZ}$ simplify toward self-similar slope and constant intrinsic scatter (5 and 10%, respectively) nearly independent of scale, making both measures ideal for cluster finding and characterization to $z=2$, (3) halo mass-conditioned likelihoods of hot gas mass and thermal energy at fixed halo mass closely follow a log-normal form, (4) despite normalization differences ranging up to $0.4$ dex, there is good qualitative, and often quantitative, agreement between the scale-dependent slopes and property covariance of the two simulations. Slopes show appreciable redshift dependence at the group scale, while redshift dependence of the scatter is exhibited by low mass FLAMINGO halos only, (5) property correlations are largely consistent between the simulations, with values that mainly agree with existing empirical measurements. We close with a literature survey placing our MPR slopes and intrinsic scatter estimates into context.

Polarimetric radio observations of the Sun can provide rich information about emission mechanisms and the propagation medium. For the past five decades, solar polarimetric studies at low radio frequencies have almost always assumed the absence of linear polarization. This has been based on the expectations from coronal propagation effects. Here we present the first robust evidence of linear polarization from solar emissions at meter wavelengths using simultaneous measurements with two telescopes of very different designs separated by thousands of kilometers - the Murchison Widefield Array and the upgraded Giant Metrewave Radio Telescope. Both datasets show consistent linear polarization fractions, confirming this detection. Rapid changes in morphology, as well as the fractional linear polarization at small time and frequency spans, further rule out any possibilities of an instrumental origin. Assuming the absence of linear polarization in solar radio emissions can result in incorrect interpretation of solar observations as well as those of other flare stars, which are often guided by learnings from solar studies. This discovery highlights the need for relaxing this assumption, and is essential for precise estimation of polarization signatures, ultimately leading to a better understanding of the plasma conditions in the Sun and other stars.

The analogues of $\gamma$ Cassiopea are binary early type Be stars which are X-ray bright and have hard thermal X-ray spectra. The nature of low-mass companions in these stars and mechanisms of their X-ray emission remain enigmatic. Among the proposed ideas is the presence of an accretion disc around a white dwarf (WD) companion to the Be star donor. We use modern radiative transfer models accounting for reflection physics in order to calculate the synthetic spectra of such systems, and assume that the hottest plasma is thermal and is located in the accretion disc boundary layer. The models are then used to analyse the archival X-ray observations of the $\gamma$ Cas analogue $\zeta$ Tau (a.k.a. Tiānguān) which were obtained by the XMM-Newton telescope. Comparisons with X-ray-emitting symbiotic systems, particularly $\delta$- and $\beta/\delta$-type systems, support the idea that the hard X-ray emission in $\zeta$ Tau is best explained by a WD accreting wind material of the Be star. The plasma temperature and luminosity of the boundary layer associated with the accretion disc are used to estimate a mass accretion rate of $\dot{M}_\mathrm{acc} \approx 4\times 10^{-10}$ M$_\odot$ yr$^{-1}$, implying a nova recurrence time above 10$^{5}$ yr. Our analysis advances the understanding the production of hard X-ray emission in $\gamma$ Cas analogues, further supporting the idea of accreting WDs as companions of Be-stars in these systems.

Gamma-ray bursts (GRBs) are the most powerful transient objects in the Universe, and they are a primary target for the MAGIC Collaboration. Recognizing the challenges of observing these elusive objects with Imaging Atmospheric Cherenkov Telescopes (IACTs), we implemented a dedicated observational strategy that included an automated procedure for rapid re-pointing to transient sources. Since 2013, this automated procedure has enabled MAGIC to observe GRBs at a rate of approximately ten per year, which led to the successful detection of two GRBs at very high energies (VHE; E > 100 GeV). We present a comprehensive analysis of 42 non-detected GRBs (4 short GRBs) observed by MAGIC from 2013 to 2019. We derived upper limits (ULs) on the observed energy flux as well as on the intrinsic energy flux corrected for absorption by the extragalactic background light (EBL) from the MAGIC observations in selected energy and time intervals. We conducted a comprehensive study of their properties to investigate the reasons for these non-detections, including the possible peculiar properties of TeV-detected GRBs. We find that strong EBL absorption significantly hinders TeV detection for the majority of GRBs in our sample. For a subset of 6 GRBs with redshift z < 2, we compared the UL on the intrinsic flux in the VHE domain with the simultaneous X-ray flux, which is observed to be at the same level in the current population of TeV-detected GRBs. Based on these inferred MAGIC ULs, we conclude that a VHE component with a luminosity comparable to the simultaneously observed X-ray luminosity cannot be ruled out for this sample.

Imaging of optical Galactic cirrus, the spatially resolved form of diffuse Galactic light, provides important insights into the properties of the diffuse interstellar medium (ISM) in the Milky Way. While previous investigations have focused mainly on the intensity characteristics of optical cirrus, their morphological properties remain largely unexplored. In this study, we employ several complementary statistical approaches -- local intensity statistics, angular power spectrum / $\Delta$-variance analysis, and Wavelet Scattering Transform analysis -- to characterize the morphology of cirrus in deep optical imaging data. We place our investigation of optical cirrus into a multi-wavelength context by comparing the morphology of cirrus seen with the Dragonfly Telephoto Array to that seen with space-based facilities working at longer wavelengths (Herschel 250~$\mu m$, WISE 12~$\mu m$, and Planck radiance), as well as with structures seen in the DHIGLS {\HI} column density map. Our statistical methods quantify the similarities and the differences of cirrus morphology in all these datasets. The morphology of cirrus at visible wavelengths resembles that of far-infrared cirrus more closely than that of mid-infrared cirrus; on small scales, anisotropies in the cosmic infrared background and systematics may lead to differences. Across all dust tracers, cirrus morphology can be well described by a power spectrum with a common power-law index $\gamma\sim-2.9$. We demonstrate quantitatively that optical cirrus exhibits filamentary, coherent structures across a broad range of angular scales. Our results offer promising avenues for linking the analysis of coherent structures in optical cirrus to the underlying physical processes in the ISM that shape them. Furthermore, we demonstrate that these morphological signatures can be leveraged to distinguish and disentangle cirrus from extragalactic light.

The interstellar comet 3I/ATLAS is only the third interstellar object to be discovered. Pre-perihelion measurements provide a unique opportunity to study its activity and composition, which may alter as it is heated in the coming months. We provide an initial baseline from optical spectroscopic observations obtained only two days after discovery, using the MUSE instrument on the VLT on 2025 July 3, while 3I was at 4.47 au from the Sun and 3.46 au from the Earth. These observations confirm the cometary nature of 3I, and reveal a red coma with a spectral slope of $(18\pm4)\%/1000$~Å, redder than most Solar System comets but similar to the surface colour of some Solar System Trans-Neptunian Objects or Centaurs. We searched for but did not detect gas emission from C$_2$, NH$_2$, CN, and [OI], which is consistent with volatile non-detections for Solar System comets at this heliocentric distance. At present, the coma appears entirely dusty. Future observations of 3I as it comes closer to the Sun will provide an invaluable opportunity to witness the evolution of its activity, study its composition, test predictions of interstellar object population models, and compare 3I to Solar System comets.

We present the design for HALO7D-X, a survey of the stellar halo to investigate the accretion history of the Milky Way. The survey will use a combination of Hubble Space Telescope (HST) and Gaia data for sky position and proper motions of faint stars (18<G<21.5 mag), while line-of-sight velocity, distance, [Fe/H], and [alpha/Fe] will be measured using follow-up Keck spectroscopy. The survey will cover 30 lines of sight, made up of multiple HST archival fields and optimized for Keck DEIMOS spectroscopy. We use mock survey observations of the Bullock and Johnston stellar halo simulations to investigate the sensitivity of HALO7D-X to constrain the basic parameters of the accretion history of our Galaxy's stellar halo. We find that we are sensitive to the mass distribution and accretion timeline of the stellar halo progenitors, but not their orbital circularity. We find that the simulated halos fall into three different groups based on the similarities in their distributions of the observable dimensions of our survey. These groups are also distinct from each other in the mass distribution and accretion timeline of their progenitor satellites, showing that by using similarities in our observables among halos, we are able to identify similarities in their accretion histories. With HALO7D-X we will compare real Milky Way data with simulated halos and use this connection between observables and progenitor mass and accretion timeline to learn about the formation of our Galaxy's stellar halo.

Accreting supermassive black holes (SMBHs) in galaxy mergers with separations $<$ 1kpc are crucial to our understanding of SMBH growth, galaxy evolution, and the evolution of SMBH binaries. Despite their importance, there are less than a handful known, and most have been discovered serendipitously. In this work, we employ a new selection method to systematically pre-select candidate advanced mergers likely to contain unresolved substructure at sub-arcsecond scales. By exploiting the large survey area and astrometric precision of the Wide-field Infrared Survey Explorer (WISE) and the Sloan Digital Sky Survey (SDSS), we have identified a sample of 48 nearby advanced mergers that have red WISE colors ($W_1-W_2>0.5$) indicative of accretion activity and significant sub-arcsecond offsets between their optical and infrared coordinates as measured by SDSS and WISE. We conducted high resolution adaptive optics (AO) observations of this sample with the Keck NIRC2 camera in the $K_p$ band ($2.124 ~ \mu m$, $\Delta\lambda = 0.351 \mu m$) to search for evidence of previously unresolved substructure suggested by the optical-to-infared offsets. We find that a significant fraction (20/48 or 42%) of the sample shows substructure tracing the SDSS/WISE offset and unresolved by SDSS, demonstrating that our methodology is a promising pathway to find dual AGN in follow-up spectroscopy. Archival optical Hubble Space Telescope (HST) imaging reveals that substructure identified with Keck is often missed in the optical or erroneously identified due to partial obscuration, underscoring the importance of carrying out studies of late-stage mergers in the infrared.

We describe the physical characteristics of interstellar comet 3I\ATLAS, discovered on 2025 July 1 by the Asteroid Terrestrial-impact Last Alert System. The comet has eccentricity, $e$ $\simeq$ 6.08 and velocity at infinity, v$_{\infty}$ $\simeq$ 57 km/s, indicating an interstellar origin. We obtained g, r, i photometry with the Palomar 200-inch Next Generation Palomar Spectrograph on 2025 July 3. We measured colour indices g-r = 0.43$\pm$0.02 mag, r-i = 0.16$\pm$0.02 mag, and g-i = 0.59$\pm$0.03 mag and a spectral slope of 1.3$\pm$0.9 $\%$/100 nm. We calculate the dust cross-section within 10,000 km of the comet to be 230.0$\pm$5.2 km$^2$, assuming an albedo of 0.10. The FWHM of 3I\ATLAS's detection is $\sim$2.2 arcsec as measured in our r-band images and has A(0$^\circ$)f$\rho$ of 287.2$\pm$2.8 cm. Finally, we use the sunward extent of the coma to constrain the dust ejection speed, finding that \textmu m-scale to mm-scale particles have $\sim$0.01-1 m/s, implying the comet's dust mass-loss rate is $\sim$0.1 - 1.0 kg/s.

Gravitational wave (GW) astronomy has opened a new window on the Universe, allowing to obtain constraints on dark energy and gravity independent from other electromagnetic waves observations, such as large scale structure (LSS). For the purpose of investigating the consistency between different observations the effective field theory (EFT) of dark energy is a useful tool, allowing to derive model and parametrization independent consistency relations (CR) between the effective gravitational constant, the slip parameter, the gravitational and electromagnetic luminosity (EM) distance, and the speed of GWs. We test the constant brading CR, which is also satisfied by general relativity, by comparing for the first time the constraints obtained from LSS observations with those from GW events with and without electromagnetic counterparts, confirming the validity of the CR at the current level of experimental uncertainty. The event GW170817 and its electromagnetic counterpart provides a constraint of the effective gravitational constant with an accuracy comparable with LSS constraints, while the analysis of GW events without electromagnetic counterpart are consistent, but do not have a constraining power comparable to LSS observations. Beside allowing to test the consistency between independent observations, the CRs can be used to estimate the effective gravitational coupling with GWs at high redshift, where other observations are not available.

KM3NeT has recently reported an event where a muon of energy $120^{+110}_{-60}$ PeV was observed at its ARCA detector, which can stem from a very high-energy neutrino interaction in the vicinity of the detector. Besides revolutionizing our understanding of high-energy neutrino sources, this event can serve as a valuable probe for studying Beyond the Standard Model (BSM) interactions of neutrinos. In this work, we study the dark matter (DM)-neutrino interaction by assuming the neutrino for the event KM3-230213A is originated from a blazar. The flux of such neutrinos, traveling through DM distributed across astrophysical and cosmological scales, can get attenuated due to DM interactions. The detection of such event by KM3NeT allows us to place constraints on the interaction cross section at highest-ever neutrino energy. We derive both conservative constraints-neglecting flux attenuation from the host halo-and optimistic ones by including host halo contributions. Our results show that the energy-independent constraints are weaker than previous bounds. For energy-dependent case, the extreme energy of the event allows us to set some of the strongest limits on scattering cross sections. In future, more such neutrino events with well-understood origin will be essential in constraining or potentially discovering DM-neutrino interactions.

The rapid expansion of mega-constellations in low Earth orbits has posed significant challenges to space traffic management, necessitating periodic inspections of satellites to ensure the sustainability of the space environment when economically feasible. This study addresses the orbital design challenge associated with inspecting numerous satellites distributed across multiple orbital planes through flybys by proposing an innovative orbital-plane-based inspection strategy. The proposed methodology reformulates the multi-satellite flyby problem into a multi-rendezvous trajectory planning problem by proposing an analytical approach to determine a maneuver-free inspection orbit that enables flyby of all satellites within a specific orbital plane. Additionally, a three-layer global optimization framework is developed to tackle this problem. The first layer establishes an approximate cost evaluation model for orbital plane visitation sequences, utilizing a genetic algorithm to identify the optimal sequence from a vast array of candidate planes, thereby maximizing inspection targets while minimizing fuel consumption. The second layer constructs a mixed-integer programming model to locally refine the rendezvous epochs and orbital parameters of each inspection orbit to reduce the total velocity increment. The third layer accurately computes the optimal impulsive maneuvers and trajectories between inspection orbits. In contrast to traditional low-Earth orbit rendezvous optimization frameworks, the proposed framework fully leverages the adjustable freedom in inclination and right ascension of the ascending node (RAAN) of inspection orbits, significantly reducing the total velocity increment. Simulation results demonstrate that the proposed method can effectively address the trajectory optimization problem associated with constellation inspection for tens of thousands of satellites.

It has recently been revealed that charged scalar clouds, spatially regular matter configurations which are made of linearized charged massive scalar fields, can be supported by spinning and charged Kerr-Newman black holes. Using analytical techniques, we establish a no-short hair theorem for these stationary bound-state field configurations. In particular, we prove that the effective proper lengths of the supported charged massive scalar clouds are bounded from below by the remarkably compact dimensionless relation $\ell/M>\ln(3+\sqrt{8})$, where $M$ is the mass of the central supporting black hole. Intriguingly, this lower bound is universal in the sense that it is valid for all Kerr-Newman black-hole spacetimes [that is, in the entire regime $\{a/M\in(0,1],Q/M\in[0,1)\}$ of the dimensionless spin and charge parameters that characterize the central supporting black holes] and for all values of the physical parameters (electric charge $q$, proper mass $\mu$, and angular harmonic indexes $\{l,m\}$) that characterize the supported stationary bound-state scalar fields.

The fluid behavior of the solar wind is affected by the heat flux carried by the suprathermal electron populations, especially the electron strahl (or beam) that propagates along the magnetic field. In turn, the electron strahl cannot be stable, and in the absence of collisions, its properties are regulated mainly by self-generated instabilities. This paper approaches the description of these heat-flux instabilities in a novel manner using regularized Kappa distributions (RKDs) to characterize the electron strahl. RKDs conform to the velocity distributions with suprathermal tails observed in situ, and at the same time allow for consistent macromodeling, based on their singularity-free moments. In contrast, the complexity of RKD models makes the analytical kinetic formalism complicated and still inaccessible, and therefore, here heat-flux instabilities are resolved using the advanced solver ALPS. Two primary types of instabilities emerge depending on plasma conditions: the whistler and firehose heat-flux instabilities. The solver is successfully tested for the first time for such instabilities by comparison with previous results for standard distributions, such as Maxwellian and Kappa. Moreover, the new RKD results show that idealized Maxwellian models can overrate or underestimate the effects of these instabilities, and also show differences from those obtained for the standard Kappa, which, for instance, underestimate the firehose heat-flux growth rates.

Dynamical systems are ubiquitous within science and engineering, from turbulent flow across aircraft wings to structural variability of proteins. Although some systems are well understood and simulated, scientific imaging often confronts never-before-seen dynamics observed through indirect, noisy, and highly sparse measurements. We present NeuralDMD, a model-free framework that combines neural implicit representations with Dynamic Mode Decomposition (DMD) to reconstruct continuous spatio-temporal dynamics from such measurements. The expressiveness of neural representations enables capturing complex spatial structures, while the linear dynamical modes of DMD introduce an inductive bias that guides training and supports stable, low-dimensional representations and forecasting. We validate NeuralDMD on two real-world problems: reconstructing near-surface wind-speed fields over North America from sparse station observations, and recovering the evolution of plasma near the Galactic-center black hole, Sgr A*. In both cases, NeuralDMD outperforms established baselines, demonstrating its potential as a general tool for imaging dynamical systems across geoscience, astronomy, and beyond.

Several major open problems in cosmology, including the nature of inflation, dark matter, and dark energy, share a common structure: they involve spacetime-filling media with unknown microphysics, and can be probed so far only through their gravitational effects. This observation motivates a systematic open-system approach to cosmology, in which gravity evolves in the presence of a generic, unobservable environment. In this work, we develop a general framework for open gravitational dynamics based on general relativity and the Schwinger-Keldysh formalism, carefully addressing the nontrivial constraints imposed by diffeomorphism invariance. At the quantum level, our path integral formulation computes the gravitational density matrix in perturbation theory around a semi-classical spacetime. As illustrative applications, we study inflation and the propagation of gravitational waves in classical regimes where environmental interactions are non-negligible. In the inflationary context, our framework reproduces the known Open Effective Field Theory of Inflation in the decoupling limit and extends it to include gravitational interactions. For gravitational waves, we derive the most general conservative and dissipative corrections to propagation. Remarkably, we find that the leading-order gravitational birefringence is dissipative in nature, whereas conservative birefringence appears only at higher derivative order, opposite to the electromagnetic case. Our results pave the way to modeling dissipative effects in the late universe.

We study the dynamics and relativistic precessions of massive particles on spherical orbits around Kerr-MOG black holes in scalar-tensor-vector gravity (STVG). By employing the Hamilton-Jacobi formalism, we derive conserved quantities and analyze how the MOG parameter $\alpha$ and orbital tilt angle $\zeta$ influence the innermost stable spherical orbits (ISSOs) and orbital stability. We compute the nodal and periastron precession frequencies, finding that nodal precession increases monotonically with both black hole spin and MOG parameter, while periastron precession exhibits a more complex behavior: MOG amplifies curvature-induced effects, which can be partially counteracted by spin. Furthermore, to complement the orbital analysis, we examine the Lense-Thirring spin precession of a gyroscope and demonstrate its sensitivity to the MOG parameter, spin, and orbital tilt angle. These results reveal distinctive signatures of modified gravity in orbital dynamics and provide a potential observational probe to test deviations from general relativity near rotating black holes.

Ultralight bosonic fields around a rotating black hole can extract energy and angular momentum through the superradiant instability and form a dense cloud. We investigate the scenario involving two scalar fields, $\phi$ and $\chi$, with a coupling term $\frac{1}{2}\lambda\phi\chi^2$, which is motivated by the multiple-axion framework. The ultralight scalar $\phi$ forms a cloud that efficiently produces $\chi$ particles non-perturbatively via parametric resonance, with a large mass hierarchy, $\mu_\chi \gg \mu_\phi$. Rather than escaping the system as investigated by previous studies, these $\chi$ particles remain bound, orbiting the black hole. Moreover, the particle production occurs primarily at the peak of the cloud's profile, allowing $\chi$ particles in quasi-circular orbits to pass repeatedly through resonant regions, leading to sustained amplification. This selective process naturally forms a dense ring of $\chi$ particles, with a mass ratio to the cloud fixed by $(\mu_\phi/\mu_\chi)^2$. Our findings reveal a novel mechanism for generating bound-state particles via parametric resonance, which also impacts the evolution of the cloud. This process can be probed through its imprint on binary dynamics and its gravitational-wave signatures.

Axions that couple to electromagnetism are produced in the early Universe by, among other channels, freeze-in via the Primakoff process. For sufficiently large axion masses, the same coupling causes the axions to decay into two photons, which subsequently ionize the intergalactic medium. If this decay occurs in the redshift range $20 \lesssim z \lesssim 1100$, then the contribution to the cosmic microwave background optical depth $\tau_{\rm reio}$ can lead to a conflict with observations, excluding models with sufficiently strongly coupled, heavy axions and high reheating temperatures, $T_{\rm reh}$. Using large ensembles of explicit type IIB string theory models with up to $h^{1,1} = 100$ axions, we compute the full cosmic reionization history caused by the decays of multiple axions. We compare this to the posterior on the high-$z$ component of $\tau_{\rm reio}$ derived from model-independent constraints on the ionization state of the Universe, obtained in a full \textit{Planck} analysis presented in a companion paper. For $h^{1,1} = 20, 50, 100$, we find that approximately 15\%, 15\%, and 10\% of the models in the ensemble prefer $T_{\rm reh} \lesssim 10^{10}\,\text{GeV}$ at 95\% CL. We provide a publicly available code at:~\href{this https URL}{this http URL}, which computes the reionization history for arbitrary ensembles of decaying axions. Our analysis opens the door for future large-scale work studying the preference for low-temperature reheating in models with multiple axions.

This study designs, implements, and evaluates a teaching strategy based on computational thinking (CT) 27 to strengthen the teaching of stellar astrophysics in secondary education. Thirty-one 11th-grade students 28 from the Escuela Normal Superior Leonor Álvarez Pinzón (Tunja, Colombia) were worked with using an 29 action research approach. The methodology combined documentary analysis and practical application 30 through a teaching sequence of four workshops, integrating CT skills such as abstraction, decomposition, 31 and algorithmic thinking in experimental activities and computational modeling. The evaluation, based on 32 rubrics and performance matrices, showed improvements in the structuring of thinking and conceptual 33 understanding of astrophysical phenomena. The results suggest that CP facilitates the teaching of physics 34 by promoting active and structured learning. It is recommended that the strategy be replicated in different 35 educational contexts to validate its applicability and optimize its impact.

A supercooled phase transition in a nearly conformal dark sector can provide a natural setting for darkogenesis via its out-of-equilibrium dynamics, where a particle-antiparticle number asymmetry in the dark sector can be reprocessed into the visible sector, yielding the observed baryon asymmetry and an asymmetric dark matter. We consider a scenario where the number asymmetry is generated from the decay of a mother particle produced via parametric resonance during the phase transition induced due to its coupling to the dilaton associated with spontaneous breaking of scale invariance. It is shown that the correct baryon asymmetry and dark matter abundance can be realized for a dark phase transition at $\mathcal{O}(1) \, \rm GeV$, which can also explain the nano-Hz gravitational wave signal reported by pulsar timing array experiments. The scenario will be tested further in neutron-antineutron oscillation experiments.

We investigate the impact of dark matter condensates on the emission and thermodynamic properties of accretion disks, in a spherically-symmetric and static background. We focus on a class of models where dark matter originates from a genuine mass mixing among neutrino fields and compute the corrections to the dark matter's potential within galactic halo. We find a corresponding Yukawa correction induced by the dark matter energy-momentum tensor over the Newtonian potential. In so doing, employing Schwarzschild coordinates, and adopting the Novikov-Thorne formalism, we compute the geodesic structure and the corresponding disk-integrated luminosity profiles. Assuming a constant mass accretion rate, constituted solely by baryonic matter, we find non-negligible deviations in both the disk structure and radiative output, as compared to the standard Schwarzschild case. Afterwards, we discuss physical consequences of our Yukawa correction, comparing it with recent literature, predicting similar potentials, albeit derived from extended theories of gravity. Accordingly, we thus speculate to use our results to distinguish among candidates of dark matter. Indeed, our findings suggest that incoming high-precision observations of accretion disk spectra may provide a tool to probe dark matter's nature under the form of particles, extended theories of gravity or condensates.

In the early 1970s, Jacob Bekenstein proposed that black holes possess an entropy proportional to the area of their event horizon, introducing the Generalized Second Law of thermodynamics. Stephen Hawking initially objected to this idea, but his subsequent analysis of quantum field theory in curved spacetime led to the prediction of Hawking radiation and the concept of black hole temperature. This study offers a concise overview of the development of black hole thermodynamics between 1972 and 1975, highlighting the theoretical evolution of both Bekenstein's and Hawking's contributions. The work also reflects on the lasting impact of these developments on modern theoretical physics and quantum gravity.

We investigate a generalized power-law dark energy equation of state of the form \( p = w\rho - \beta\rho^m \) in a flat FLRW universe, analyzing its dynamical stability and thermodynamic consistency. The model exhibits a rich phase space structure, with an effective cosmological constant \(\rho^* = [(1+w)/\beta]^{1/(m-1)}\) emerging as a stable attractor for \((w < -1,~ m > 1)\). Notably, the universe evolves from an early de Sitter phase (\(w \to -1\)) to a late-time de Sitter-like one with phantom crossing (\(w(z) < -1\)), aligning with DESI observations. Dynamical analysis reveals that the \(m > 1\) regime avoids ghost instabilities while accommodating phantom behavior, with \(m = 2\) providing particular theoretical advantages. Thermodynamically, the Generalized Second Law holds when the null energy condition \(\rho + p \geq 0\) is satisfied, which naturally occurs for \(\rho \geq \rho^*\). The model's compatibility with both observational data and fundamental thermodynamic principles suggests it as a viable framework for describing late-time cosmic acceleration, resolving tensions associated with phantom crossing while maintaining entropy dominance of the cosmological horizon.

We have investigated the bright ring-like features and polarization structures in the Kerr-Sen black hole's images based on general relativistic radiative transfer (GRRT) simulations, which are illuminated by the 230 GHz thermal synchrotron emission from the radiatively inefficient accretion flows (RIAF). Using the REx method, we extracted the bright ring from the black hole images and found that both the bright ring and black hole shadow shrink as the dilation parameter increases. Combing with the Event Horizon Telescope (EHT) observational data of SgrA*, we present the allowed ranges of black hole parameters for different observed inclinations and disk thicknesses. Our results show that effects of the disk thickness on the allowed parameter space are more strong than those of the observed inclination. We also probe the rotational symmetry of EVPA in the polarization structure of the black hole images by analyzing the coefficient $\beta_2$. The dilaton parameter results in that the real part $Re\beta_2$ increases and the imaginary part $Im\beta_2$ decreases, but the magnitude $|\beta_2|$ generally exhibits a declining trend. Finally, we find that effects of the disk thickness on $\beta_2$ are much weaker than those from the dilaton parameter.

The introduction of General Relativity (GR) in 1915 revolutionized our understanding of gravity, but over time, its limitations in explaining phenomena like dark energy, dark matter, and quantum gravity have motivated alternative theories. Early modifications, such as Weyl's 1919 proposal, focused on adding higher-order terms to the Einstein-Hilbert action. GR's non-renormalizability further strengthened the case for extending it. A central theme of modern gravity research is modifying the geometric structure, often by changing the gravitational Lagrangian. This leads to theories such as teleparallel and symmetric teleparallel gravity, utilizing torsion or non-Levi-Civita connections, with differential geometry providing the essential framework. This thesis explores several modified gravity models. Chapter 1 introduces necessary mathematical tools. Chapter 2 develops a novel parametrization of the deceleration parameter, constrained using MCMC and observational data, and applies it to f(Q) gravity. Chapter 3 embeds the LambdaCDM model into f(Q, L\_m) gravity with non-minimal coupling, producing analytic solutions and matching observations through cosmographic analysis. Chapter 4 considers Bianchi-I spacetime in f(R, L\_m) gravity with observational constraints to measure anisotropy. Chapter 5 presents wormhole solutions in f(Q, T) gravity with conformal symmetries. Chapter 6 explores wormholes in f(R, L\_m) with non-commutative geometry, analyzing shape functions, energy conditions, and stability. Chapter 7 studies Big Bang Nucleosynthesis in f(T) gravity, constraining hybrid models using early- and late-time data, and validating intermediate epochs via cosmography.

Independently from the formation mechanism of primordial black holes in the early universe, their generation is accompanied by a ringdown phase during which they relax to a stationary configuration and gravitational waves under the form of quasinormal modes are emitted. Such gravitational waves generate an irreducible and unavoidable stochastic background which is testable by current and future experiments. In particular, the presence of a primordial black hole ringdown allows to significantly constrain the idea of extremely supermassive black holes to comprise the entire dark matter of the universe.

We examine the behaviour of the gauge invariant scalar field perturbations in an analytic inflationary model that transitions from slow roll to an ultra-slow roll (USR) phase. We find that the numerical solution of the Mukhanov-Sasaki equation is well described by Hamilton-Jacobi (HJ) theory, as long as the appropriate branches of the Hamilton-Jacobi solutions are invoked: Modes that exit the horizon during the slow roll phase evolve into the USR as described by the first HJ branch, up to a subdominant $\mathcal{O}(k^2/H^2)$ correction to the Hamilton-Jacobi prediction for their final amplitude that we compute, indicating the influence of neglected gradient terms. Modes that exit during the USR phase are described by a separate HJ branch (very close to de Sitter) once they become sufficiently superhorizon, obtained by the shift $\epsilon_2 \rightarrow \tilde{\epsilon}_2\simeq \epsilon_2+6$. This transition is similar to the conveyor belt concept put forward in our previous work [1] and the intuition deriving from it. We further discuss implications for the validity of the stochastic equations arising from the Hamilton-Jacobi formulation. We conclude that if appropriately used, the Hamilton-Jacobi attractors can describe the dynamics of long wavelength inflationary inhomogeneities beyond slow roll.

We explore the dynamics of pure scalar fields rolling on an exponential potential in the absence of any additional background fluid and demonstrate the existence of self-tracking solutions in which the self-perturbations of the scalar field act as an effective radiation background. The validity of these solutions is demonstrated through both analytic techniques and numerical simulations using CosmoLattice. We discuss applications to string cosmologies with significant trans-Planckian field excursions between inflation and BBN, including the required initial level of scalar perturbations to avoid overshoot.

The IceCube Neutrino Observatory recently published evidence for diffuse neutrino emission from the Galactic Plane at $4.5\sigma$ significance. This new source of astrophysical neutrinos provides an exciting laboratory for probing the nature of neutrino masses. In particular, extremely small mass splittings, such as those predicted by quasi-Dirac neutrino mass models, and finite neutrino lifetimes from neutrino decays, would induce effects on the spectra and flavor ratios of neutrinos with TeV-scale energies traversing kiloparsec-scale baselines. Using $\mathtt{TANDEM}$, an upcoming three dimensional galactic neutrino emission model, we explore the sensitivity of IceCube and KM3NeT/ARCA to these ultra-long-baseline phenomena. We find that a combined analysis would be sensitive to quasi-Dirac mass splittings $10^{-14.0}~\mathrm{eV^2} \lesssim \delta m^2 \lesssim 10^{11.6}~\mathrm{eV^2}$ and neutrino lifetimes $m / \tau \gtrsim 10^{-14.1}~\mathrm{eV^2}$ at $> 1\sigma$, both regions constituting as-yet unexplored parameter space. Our results demonstrate the potential that astrophysical neutrino sources and global neutrino telescope networks have in probing new regions of exotic neutrino mass models.

We study the periapsis shift of timelike bound orbits in the Zipoy-Voorhees spacetime, which is an exact, static, axisymmetric, and vacuum solution characterized by the deformation parameter $\gamma$, including the Schwarzschild spacetime as $\gamma=1$. We derive both the exact formula for the periapsis shift of a quasi-circular orbit and the formula for the periapsis shift by the post-Newtonian (PN) expansion. We show that the periapsis shift in the Zipoy-Voorhees spacetime for $1/11 < \gamma \leq 1/5$ is the same as in the Schwarzschild spacetime to the 2PN order if $e=\sqrt{(\gamma^{-1}-5)/6}$, where $e$ is the eccentricity of the orbit. Furthermore, we show that the parameter $\gamma$ and the corresponding quadrupole moment $\tilde{M}_2$ of the supermassive compact object at Sagittarius A* are constrained to $\gamma \gtrsim 1.7 \times 10^{-2}$ and $\tilde{M}_2 \lesssim 1.2 \times 10^3$, respectively, from observational data on S2 using the obtained PN expansion formula. Finally, we derive a new series representation for the periapsis shift in the Zipoy-Voorhees spacetime using a recently proposed prescription, which shows fast convergence not only in the weak-field regime but also for small eccentricity.

We study the cosmology of a modified majoron model motivated by the need to protect a global $U(1)$ symmetry from gravity-induced hard explicit breaking (by $d \leq 4$ operators) at the Planck scale. The model extends the Standard Model by introducing a gauged $U(1)_{B-L}$ and an approximate global $U(1)$ symmetry, each spontaneously broken by a corresponding complex scalar singlet. This setup gives rise to a network of effectively global and local cosmic strings, whose stochastic gravitational wave signals can jointly account for the spectrum observed by the NANOGrav collaboration, particularly for majoron masses $m_{\chi} < 10^{-23}$ eV. Although the fit is not as strong as that from supermassive black hole mergers, the model still provides an alternative explanation rooted in high-energy physics. The model also generates light neutrino masses via the seesaw mechanism and avoids cosmological constraints from $\Delta N_{\text{eff}}$, CMB anisotropies, and isocurvature fluctuations. Although the majoron can contribute to dark matter through thermal, coherent oscillation, and string-induced production mechanisms, its relic abundance remains subdominant in the NANOGrav-compatible region. In contrast, the measured dark matter relic density is achievable at higher $m_\chi$, though at the cost of tension with cosmological bounds. If the NANOGrav fits are viewed as constraints, given their comparatively lower Bayes factors, they yield bounds that are significantly stronger than those imposed by the CMB and other cosmological data.

We study strongly correlated fractional topological phases on a two-sphere threaded by a magnetic dipole field with globally vanishing flux. Solving the Dirac equation in this background produces spheroidal wavefunctions forming a highly degenerate manifold of normalizable zero modes, with degeneracy proportional to the total absolute flux. We introduce a non-Abelian spin gauge field near the equator to hybridize the north and south domain-confined modes, forming a global flat band. Projecting interactions into this band yields Laughlin-type correlated states. The entanglement spectrum shows a chiral tower consistent with a virtual edge, demonstrating bulk-edge correspondence in a closed geometry. This generalizes the zero-flux flat-band construction of \cite{Parhizkar:2024som} to curved backgrounds, with potential applications to synthetic and astrophysical systems.

We present stationary and axially-symmetric black hole solutions to the Einstein field equations sourced by an anisotropic fluid, describing rotating black holes embedded in astrophysical environments. We compute their physical properties, including quantities associated with the circular geodesics of massless and massive particles, analyze their shadows and image features, and energy conditions. Overall, we find that deviations from the Kerr metric grow with spin.

Compact binaries with large mass asymmetries - such as Extreme and Intermediate Mass Ratio Inspirals - are unique probes of the astrophysical environments in which they evolve. Their long-lived and intricate dynamics allow for precise inference of source properties, provided waveform models are accurate enough to capture the full complexity of their orbital evolution. In this work, we develop a multi-parameter formalism, inspired by vacuum perturbation theory, to model asymmetric binaries embedded in general matter distributions with both radial and tangential pressures. In the regime of small deviations from the Schwarzschild metric, relevant to most astrophysical scenarios, the system admits a simplified description, where both metric and fluid perturbations can be cast into wave equations closely related to those of the vacuum case. This framework offers a practical approach to modeling the dynamics and the gravitational wave emission from binaries in realistic matter distributions, and can be modularly integrated with existing results for vacuum sources.

A theoretical model on the basis of fluid-Maxwell equations for an electron-ion plasma is presented which describes the conversion of current-driven Langmuir waves into type III radiation whereby simultaneously an excitation of whistler waves may occur. In contrast to the classical approach of Ginzburg and Zhelezniakov (1958) after which beam-excited Langmuir waves in a two-step process are converted in electromagnetic radiation, the presented mechanism works without any instability and wave coalescence. Rather the electric field oscillations at the electron plasma frequency can be triggered by different realisations of the driving current, e.g. by the (uncompensated) net current of the strahl at t=0 in a core-strahl plasma or by given current variations which may represent different situations in space, as shocks, magnetic switch-backs etc.. A linearized system of equations is used to describe the mode coupling occurring at oblique propagation between the mostly electrostatic Langmuir wave and the adjacent electromagnetic left-hand polarized (L) wave. The simplicity of the fluid model allows without great effort the parameters of the current profiles to be varied and thus to simulate a wide range of possible experimental conditions. Measurements of Langmuir waves, type III radiation and whistler waves on board various satellites in the solar wind, and in particular some of the recent results of the Parker Solar Probe are interpreted in the light of the theoretical model presented. For the case of the uncompensated strahl, the fluid approach is confirmed by fully kinetic PIC simulations. One comparison is shown in the Appendix.

The aim of the paper is to demonstrate that electron current oscillations may generate electromagnetic waves as type III radiation and whistler waves without the involvement of the classical plasma emission via the coalescence of waves. PIC simulation results of an electron-core-strahl plasma without initial current compensation are presented which describe the conversion of current-driven Langmuir oscillations/waves into type III radiation whereby simultaneously whistler waves are excited. In contrast to the classical approach of Ginzburg and Zhelezniakov (1958) after which beam-excited Langmuir waves in a two-step process are converted in electromagnetic radiation, any instability is suppressed by selecting a low strahl velocity. Rather electric field oscillations at the electron plasma frequency are triggered by the initially non-compensated current of the strahl. The arising electromagnetic fields exhibit amplitude oscillations which are caused by the superposition of the two wave modes of mixed polarisation at the point of mode coupling. This basic mechanism of wave generation and transformation has already been described in earlier papers using simple fluid models. It is also the topic of the companion paper. Besides the fundamental electromagnetic radiation, the second harmonic of nearly the same intensity has been obtained which is an indication for nonlinear currents. Measurements of Langmuir waves, type III radiation and whistler waves on board various satellites in the solar wind, in particular Parker Solar Probe (PSP) observations are analysed in the light of our results. Interpretations of earlier PIC simulations are critically reviewed.

The theory controlling the Universe's evolution in the classical regime has to be motivated by particle physics reasoning and should also generate inflation and dark energy eras in a unified way. One such framework is $F(R)$ gravity. In this work we examine a class of exponential deformations of $R^2$ gravity motivated by fundamental physics of scalaron evolution in a de Sitter background. As we show this class of models describe both inflation and the dark energy era in a viable way compatible with the Planck constraints on inflation and the cosmological parameters. Regarding the inflationary era, the exponentially deformed $R^2$ model also yields a rescaled Einstein-Hilbert term which remarkably does not affect the dynamics and the inflationary evolution is identical to that of an $R^2$ model. The dark energy era is also found to be viable and mimics the $\Lambda$-Cold-Dark-Matter model. More importantly, this class of $F(R)$ gravity exponential $R^2$ deformations also has an important characteristic, and specifically it yields total equation of state oscillations deeply in the matter domination era, for redshifts $z\sim 3400$, so near the matter-radiation equality. These total equation of state deformations at such a large redshift may directly affect the energy spectrum of the primordial gravitational waves. Indeed as we show, the effect is measurable and it leads to an enhancement of the tensor perturbations energy spectrum for low frequencies probed by the future LiteBIRD mission. This enhancement might have a measurable effect on the $B$-modes of the Cosmic Microwave Background radiation and thus may be detectable by the LiteBIRD mission. Only a handful of theoretical frameworks can generate the gravitational wave pattern generated by the class of exponentially deformed $R^2$ models we presented.

The sensitivity of current gravitational wave (GW) detectors to transient GW signals is severely affected by a variety of non-Gaussian and non-stationary noise transients, such as the blip and tomte "glitches". These glitches share some time-frequency resemblance with GW signals from massive binary black holes. In earlier works [Joshi et al., Phys. Rev. D 103, 044035 (2021); Choudhary et al., Phys. Rev. D 110, 044051 (2024)], the authors presented a method for constructing a $\chi^2$-distributed optimized statistic, based on the unified formalism of $\chi^2$ discriminators [Dhurandhar et al., Phys. Rev. D 96, 103018 (2017)], to distinguish the blip and tomte glitches from the compact binary coalescence (CBC) signals. Unlike past works, the new $\chi^2$ discriminator is constructed from the most significant singular vectors obtained from the singular-value decomposition (SVD) of glitches in real detector data. We find that the chi-square developed in this work performs as efficiently as in Choudhary et al. [Phys. Rev. D 110, 044051 (2024)], which used sine-Gaussian basis vectors. This result supports past empirical findings that the blips and tomtes are well-modeled by sine-Gaussians. It also introduces a method for constructing signal- and glitch-based $\chi^2$ discriminators by directly using real data containing the glitches and, thus, holds promise for extensions to glitches that are captured less well by sine-Gaussians or other analytical functions.

An argument is developed that the long-standing mystery in nuclear physics of the effective axial-current coupling constant in nuclei, $g_A^{\rm eff}\approx 1$, could be understood in terms of the mechanism referred to as ``pseudo-conformal sound speed" in dense compact-star matter, $v_{\rm pcs}^2/c^2\approx 1/3$. Both pros and cons are presented using an effective field theory anchored on renormalization-group approach to interacting baryons on the Fermi surface that enables one to go beyond Weinberg's highly successful EFT $\chi$EFT$_\pi$ with the pion field only (in nuclear medium) by implementing heavy-meson degrees of freedom. Both hidden local symmetry and hidden scale symmetry, the former involving the vector mesons $\rho$ and $\omega$ and the latter involving the hidden scalar meson, a dilaton $\hat{\sigma}$ ($f_0(500)$), play the crucial role. Going beyond the density regime applicable to normal nuclear matter $n_0$, the notion of ``hadron-quark continuity" is brought in via the topological structure of the nucleon, i.e., skyrmion considered to be valid in QCD at large $N_c$ limit. The new inputs for the argumentation are the large N limit of the Grassmanian model for hidden local symmetry and the IR fixed point in QCD for $N_f \leq 3$ involving ``genuine/QCD-conformal dilaton" for hidden scale symmetry.

The LISA mission is an international collaboration between ESA, its member states, and NASA, for the detection of gravitational waves from space. It was adopted in January 2024 and is scheduled for launch in the mid-2030's. It will be a constellation of three identical spacecraft forming a near-equilateral triangle in an heliocentric orbit, transferring laser beams over $2.5 \cdot 10^6$ km long arms. Laser interferometry is used to track separations between test masses, thus measuring spacetime strain variations as a function of time. LISA Science Objectives tackle many open questions in astrophysics, fundamental physics and cosmology, including ESA's Cosmic Vision questions "What are the fundamental laws of the universe?" and "How did the universe originate and of what is it made?". In this contribution, based on the LISA Red Book, we present a summary of the LISA Science Objectives relevant for the European Strategy for Particle Physics.

Magnetic monopoles with masses up to $10^{14}$ GeV can be accelerated to relativistic velocities in Galactic and intergalactic magnetic fields. The cosmic flux of relativistic monopoles is constrained by various experiments, with the limits given as functions of the monopole velocity (Lorentz factor) at the detectors. The velocity, however, is usually treated as a free parameter due to the ambiguity in the computation of the acceleration before the monopoles arrive at Earth. We explicitly evaluate the velocity by exploiting recent studies on cosmic magnetic fields and the monopole acceleration therein, to recast experimental limits in terms of the mass of monopoles. By applying our method to various terrestrial experiments, including the Pierre Auger Observatory, IceCube, MACRO, and the upcoming Cherenkov Telescope Array Observatory, as well as to astrophysical constraints, we report limits on the flux of monopoles for a wide range of monopole masses. We also highlight the role of monopoles as messengers of cosmic magnetic fields, and discuss the possibility of using monopole experiments to probe intergalactic magnetic fields.

A new class of exact spacetimes in Einstein's gravity, which are Kerr black holes immersed in an external magnetic (or electric) field that is asymptotically uniform and oriented along the rotational axis, is presented. These are axisymmetric stationary solutions to the Einstein-Maxwell equations such that the null directions of the Faraday tensor are not aligned with neither of the two principal null directions of the Weyl tensor of algebraic type D (unlike the Kerr-Melvin spacetime). Three physical parameters are the black hole mass $m$, its rotation $a$, and the external field value $B$. For vanishing $B$ the metric directly reduces to standard Boyer-Lindquist form of the Kerr black hole, while for zero $m$ we recover conformally flat Bertotti-Robinson universe with a uniform Maxwell field. For zero $a$ the spacetime is contained in the Van den Bergh-Carminati solution which can be understood as the Schwarzschild black hole in a magnetic field. Our family of black holes with non-aligned Maxwell hair - that can be called the Kerr-Bertotti-Robinson (Kerr-BR) black holes - may find application in various studies ranging from mathematical relativity to relativistic astrophysics.

Overlapping gravitational wave (GW) signals are expected in the third-generation (3G) GW detectors, leading to one of the major challenges in GW data analysis. Inference of overlapping GW sources is complicated - it has been reported that hierarchical inference with signal subtraction may amplify errors, while joint estimation, though more accurate, is computationally expensive. However, in this work, we show that hierarchical subtraction can achieve accurate results with a sufficient number of iterations, and on the other hand, neural density estimators, being able to generate posterior samples rapidly, make it possible to perform signal subtraction and inference repeatedly. We further develop likelihood-based resampling to accelerate the convergence of the iterative subtraction. Our method provides fast and accurate inference for overlapping GW signals and is highly adaptable to various source types and time separations, offering a potential general solution for overlapping GW signal analysis.