Papers for Tuesday, Jul 22 2025

This work concerns a detailed review of data analysis methods used for remotely sensed images of large areas of the Earth and of other solid astronomical objects. In detail, it focuses on the problem of inferring the materials that cover the surfaces captured by hyper-spectral images and estimating their abundances and spatial distributions within the region. The most successful and relevant hyper-spectral unmixing methods are reported as well as compared, as an addition to analysing the most recent methodologies. The most important public data-sets in this setting, which are vastly used in the testing and validation of the former, are also systematically explored. Finally, open problems are spotlighted and concrete recommendations for future research are provided.

The Milky Way is surrounded by large amounts of hot gas at temperatures T>10^6 K, which represents a major baryon reservoir. We here explore the prospects of studying the hot coronal gas in Milky Way halo by analyzing the highly forbidden optical coronal lines of [FeX] and [FeXIV] in absorption against bright extragalactic background sources. We use a semi-analytic model of the Milky Way's coronal gas distribution together wih HESTIA simulations of the Local Group and observational constraints to predict the expected FeX and FeXIV column densities as well as the line shapes and strengths. We predict column densities of log N(FeX)=15.40 and log N(FeXIV)=15.23 in the Milky Way's hot halo and estimate that a minimum S/N of 50,000 (25,000) is required to detect [FeX] l6374.5 ([FeXIV] l5302.9) absorption at a 3sigma level. Using archical optical data from an original sample of 739 high resolution AGN spectra from VLT/UVES and KECK/HIRES, we generate a stacked composite spectrum to measure an upper limit for the column densities of FeX and FeXIV in the Milky Way's coronal gas. No [FeX] and [FeXIV] is detected in our composite spectrum, which achieves a maximum S/N= 1,240 near 5300 A. We derive 3sigma upper column-density limits of log N(FeX)<16.27 and log N(FeXIV)<15.85, in line with the above-mentioned predictions. While [FeX] and [FexIX] absorption is too weak to be detected with current optical data, we outline how up-coming extragalactic spectral surveys with millions of medium- to high-resolution optical spectra will provide the necessary sensitivity and spectral resolution to measure velocity-resolved [FeX] and [FeXIV] absorption in the Milky Way's coronal gas (and beyond). This gives the prospect of opening a new window for studying the dominant baryonic mass component of the Milky Way in the form of hot coronal gas via optical spectroscopy.

We demonstrate a novel setup for hybrid particle-in-cell simulations designed to isolate the physics of the shock precursor over long time periods for significantly lower computational cost than previous methods. This is achieved using a "faux-shock" or shock-like boundary condition on one edge of our simulation domain such that particles that interact with the boundary either pass through it or are reflected off of it with a change in momentum that mimics scattering in the downstream. We show that our faux-shock setup reproduces the same fluid quantities and phase spaces as traditional shock simulations, including those which could otherwise only be done in 3D, with higher particle resolution and for reduced computational cost. While the method involves an assumed boundary condition, it nonetheless captures the essential physics of interest, establishing it as a reliable and efficient tool for future self-consistent studies of instabilities driven by cosmic rays in a shock upstream medium.

We study the observational signatures from the interactions between a newly born neutron star and a companion star that is impacted by the supernova ejecta. We focus on the cases with bound post-explosion orbits, where the neutron star may periodically gravitationally capture gas from the companion. We find that neutron star accretion must occur if the pre-supernova binary separation is less than about 20 Rsun. This is because the stellar radius expands beyond this radius before the shock-inflated envelope undergoes Kelvin-Helmholtz contraction back to the main sequence. We then consider the internal shocks formed between adjacent episodes of disk wind. The shocks efficiently convert the wind kinetic energy into radiation (due to inverse-Compton cooling), which heats up the supernova ejecta located at much larger radii. The extra heating powers bright optical emission that is periodically modulated on the orbital timescale. The shocks also accelerate non-thermal particles which produce gamma-ray and neutrino emission from 100 MeV to 10 PeV via hardronic pp collisions. The high-energy photons leak out of the supernova ejecta after a delay of several months to one year. Photo-ionization of the slowest parts of the disk wind produces hydrogen recombination lines. We then use the model to explain the puzzling Type Ic supernova SN2022jli which shows a double-peaked optical lightcurve along with many peculiar properties, including delayed onset and rapid shutoff of the second peak, periodic modulations, delayed GeV emission, and narrow Balmer lines. Under this model, SN2022jli had a close-by companion at a pre-supernova binary separation of 10 to 20 Rsun, likely due to an earlier phase of common-envelope evolution.

We present the joint modelling of weak-lensing and count measurements of the galaxy clusters detected with the AMICO code, in the fourth data release of the Kilo Degree Survey (KiDS-1000). The analysed sample comprises about 8000 clusters, covering an effective area of 839 deg$^{2}$ and extending up to a redshift of $z = 0.8$. Stacked cluster weak-lensing and count measurements have been derived in bins of redshift and intrinsic richness, $\lambda^*$. Based on self-organising maps, we reconstructed the true redshift distributions of the background galaxy samples. We accounted for the systematic uncertainties arising from impurities in the background and cluster samples, biases in the cluster $z$ and $\lambda^*$, projection effects, halo orientation and miscentring, truncation of cluster halo mass distributions, matter correlated with cluster haloes, multiplicative shear bias, baryonic matter, geometric distortions in the lensing profiles, uncertainties in the theoretical halo mass function, and super-sample covariance. We also employed a blinding strategy based on perturbing the cluster sample completeness. The improved statistics and photometry compared to the previous KiDS data release, KiDS-DR3, have led to a halving of the uncertainties on $\Omega_{\rm m}$ and $\sigma_8$, as we obtained $\Omega_{\rm m}=0.22\pm0.02$ and $\sigma_8=0.86\pm0.03$. The constraint on $S_8 \equiv \sigma_8(\Omega_{\rm m}/0.3)^{0.5}$, $S_8=0.74\pm0.03$, is in excellent agreement with recent cluster count and KiDS-1000 cosmic shear analyses, while it shows a $2.8\sigma$ tension with Planck cosmic microwave background results. The constraints on the $\log\lambda^*-\log M_{200}$ relation imply a mass precision of 8%, on average. In addition, the result on the intrinsic scatter of the $\log\lambda^*-\log M_{200}$ relation, $\sigma_{\rm intr}=0.05\pm0.02$, confirms that $\lambda^*$ is an excellent mass proxy.

Gamma-ray bursts (GRBs) are singular outbursts of high-energy radiation with durations typically lasting from milliseconds to minutes and, in extreme cases, a few hours. They are attributed to the catastrophic outcomes of stellar-scale events and, as such, are not expected to recur. Here, we present observations of an exceptional GRB\,250702BDE which triggered the {\em Fermi} gamma-ray burst monitor on three occasions over several hours, and which was detected in soft X-rays by the \textit{Einstein Probe} a day before the $\gamma$-ray triggers (EP250702a). We present the discovery of an extremely red infrared counterpart of the event with the VLT, as well as radio observations from MeerKAT. Hubble Space Telescope observations pinpoint the source to a non-nuclear location in a host galaxy with complex morphology, implying GRB 250702BDE is an extragalactic event. The multi-wavelength counterpart is well described with standard afterglow models at a relatively low redshift $z \sim 0.2$, but the prompt emission does not readily fit within the expectations for either collapsar or merger-driven GRBs. Indeed, a striking feature of the multiple prompt outbursts is that the third occurs at an integer multiple of the interval between the first two. Although not conclusive, this could be indicative of periodicity in the progenitor system. We discuss several possible scenarios to explain the exceptional properties of the burst, which suggest that either a very unusual collapsar or the tidal disruption of a white dwarf by an intermediate-mass black hole are plausible explanations for this unprecedented GRB.

We present AstroWISP: a collection of image processing tools for source extraction, background determination, point spread function/pixel response function fitting, and aperture photometry. AstroWISP is particularly well-suited for working with detectors featuring a Bayer mask (an array of microfilters applied to each detector pixel to allow color photography), such as consumer DSLR cameras. Such detectors pose significant challenges for existing tools while offering a much cheaper alternative to specialized devices. As a result, consumer DSLR cameras with Bayer masks are often underutilized for precision photometry. \astrowisp{} addresses this limitation in an effort to democratize precision photometry and support broader community participation in research. We demonstrate that our tools produce high-precision photometry from such images, enabling the use of such devices for detecting exoplanet transits. We package our tools for all major operating systems to ensure accessibility for amateur astronomers.

A popular model for Type Ia supernovae (SNe Ia) is the detonation of a CO white dwarf (WD) that is triggered by the prior detonation of a thin surface layer of helium, known as a double detonation (DD). We explore the unique early electromagnetic signatures that are expected from collision of the CO detonation with the He detonation. The three features are (1) a shock breakout flash, (2) a stage of planar shock breakout cooling, and finally (3) shock cooling emission from the thermal energy released by the collision. The planar phase is unique to the unusual density profile of the He-detonated layer in comparison to the steep profile at a stellar edge as is usually considered for shock breakout. The shock cooling emission can be modified by recombination, and we explore these effects. All together, we expect an initial flash dominated by the planar phase of $\sim6\times10^{43}\,{\rm erg\,s^{-1}}$, which lasts ~5 s in the soft X-rays. This is followed by ~12-24 hrs of shock cooling at a luminosity of $3-10\times10^{40}\,{\rm erg\,s^{-1}}$ in the optical/UV. We discuss prospects for detection of this early DD emission with current and upcoming surveys.

The broad emission lines (BELs) emitted by Active Galactic Nuclei respond to variations in the ionizing continuum emission from the accretion disk surrounding the central supermassive black hole (SMBH). This reverberation response provides insights into the structure and dynamics of the Broad Line Region (BLR). In Rosborough et al., 2024, we introduced a new forward-modeling tool, the Broad Emission Line MApping Code (BELMAC), which simulates the velocity-resolved reverberation response of the BLR to an input light curve. In this work, we describe a new version of BELMAC, which uses photoionization models to calculate the cloud luminosities for selected BELs. We investigated the reverberation responses of H$\alpha$, H$\beta$, MgII$\lambda$2800 and CIV$\lambda$1550 for models representing a disk-like BLR with Keplerian rotation, radiatively driven outflows, and inflows. The line responses generally provide a good indication of the respective luminosity-weighted radii. However, there are situations when the BLR exhibits a negative response to the driving continuum, causing overestimates of the luminosity-weighted radius. The virial mass derived from the models can differ dramatically from the actual SMBH mass, depending mainly on the disk inclination and velocity field. In single zone models, the BELs exhibit similar responses and profile shapes; two-zone models, such as a Keplerian disk and a biconical outflow, can reproduce observed differences between high- and low-ionization lines. Radial flows produce asymmetric line profile shapes due to both anisotropic cloud emission and electron scattering in an inter-cloud medium. These competing attenuation effects complicate the interpretation of profile asymmetries.

Context. Galaxy clusters provide key insights into cosmic structure formation, galaxy formation and are essential for cosmological studies. Aims. We present a catalog of galaxy clusters detected in the Kilo-Degree Survey (KiDS-DR4) optimized for cosmological analyses and investigations of cluster properties. Each detection includes probabilistic membership assignments for the KiDS-DR4 galaxies within the magnitude range 15<r'<24. Methods. Using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm, we identified 23965 clusters over an effective area of about 839 deg2 in the redshift range $0.1\le z \le0.9$, with a signal-to-noise ratio $S/N>3.5$. The sample is highly homogeneous across the entire survey thanks to the restrictive galaxy selection criteria we adopted. Spectroscopic data from the GAMA survey were used to calibrate the clusters photometric redshift and assess their uncertainties. We introduced algorithmic enhancements to AMICO to mitigate border effects among neighbor tiles. Quality flags are also provided for each cluster detection. The sample purity and completeness assessments have been estimated using the SinFoniA data driven approach, thus avoiding strong assumptions embedded in numerical simulations. We introduced a blinding scheme of the selection function meant to support the cosmological analyses. Results. Our cluster sample includes 321 cross-matches with the X-ray eRASS1 "primary" sample and 235 matches with the ACT-DR5 cluster sample. We derived a mass-proxy scaling relation based on intrinsic richness, $\lambda_*$, using masses from the eRASS1 catalog. Conclusions. The KiDS-DR4 cluster catalog provides a valuable data set for investigating galaxy cluster properties and contributes to cosmological studies by offering a large, well-characterized cluster sample.

We present a crossmatch between a combined catalog of X-ray sources and the Vera C. Rubin Observatory Data Preview 1 (DP1) to identify optical counterparts. The six fields targeted as part of DP1 include the Extended Chandra Deep Field South (E-CDF-S), the Euclid Deep Field South (EDF-S), the Fornax Dwarf Spheroidal Galaxy (Fornax dSph), 47 Tucanae (47 Tuc) and science validation fields with low galactic and ecliptic latitude (SV\_95\_-25 and SV\_38\_7, respectively). We find matches to 2314 of 3830 X-ray sources. We also compare our crossmatch to DP1 in the E-CDF-S field to previous efforts to identify optical counterparts. The probability of a chance coincidence match varies across each DP1 field, with overall high reliability in the E-CDF-S field, and lower proportion of high-reliability matches in the other fields. The majority of previously known sources that we detect are, unsurprisingly, active galaxies. We plot the X-ray-to-optical flux ratio against optical magnitude and color in an effort to identify Galactic accreting compact objects using a {\em Gaia} color threshold transformed to LSST $g$--$i$, but do not find any strong candidates in these primarily extragalactic counterparts. The DP1 dataset contains high-cadence photometry collected over a number of nights. We calculate the Stetson \( J \) variability index for each object under the hypothesis that X-ray counterparts tend to exhibit higher optical variability; however, the evidence is inconclusive whether our sample is more variable over DP1 timescales when compared to field objects.

Tidal disruption events (TDEs) as excellent beacons to black hole (BH) accreting systems have been studied for more than five decades with a single star tidally disrupted by a central massive BH. However, if considering two stars passing through a central BH and being tidally disrupted in a short period, so-called double TDEs could be expected and lead to unique variability features very different from features from standard TDEs. Here, we report such oversimplified double TDEs in the known changing-look AGN Mrk1018, of which 15years-long light curve with plateau features can be described by two main-sequence stars tidally disrupted by the central supermassive BH. Meanwhile, the BH mass determined by the double TDEs is consistent with the M-sigma relation determined value by the measured stellar velocity dispersions in Mrk1018. The results indicate tight connections between the TDEs and the changing-look properties in Mrk1018.

We present the analysis of the extended Main Sequence Turn-Off (eMSTO) in the open cluster NGC\,6067. We derive the projected rotational velocity, \textit{v}sin\textit{i}, of the stars belonging to the eMSTO region of the main sequence (MS) utilizing \textit{Gaia}-ESO spectra. Our results reveal a positive correlation between \textit{v}sin\textit{i} and the color of eMSTO stars, where fast-rotating stars predominantly occupy the red part of the MS while slow-rotating ones prefer a bluer side of the MS. The gravity-darkening effect might be a reason for this correlation. We find that most of the close binaries present in the eMSTO population would be slow-rotating due to the tidal-locking phenomenon. We identify four double-lined spectroscopic binaries (SB2) featuring slow-rotating companions, further supporting this tidal-locking hypothesis. However, the spatial distribution and the cumulative radial distribution indicate a higher concentration of red eMSTO stars in the cluster's central region than their bluer counterparts. This suggests that tidal locking is less likely to be the cause of the observed spread in rotation rates among eMSTO stars. Instead, we propose that star-disk interactions during the pre-main-sequence phase might have played a crucial role in spreading the rotation rates of stars, leading to the eMSTO phenomenon in NGC\,6067.

High-mass stars are notable for several reasons: they are characterized by strong winds, which inject momentum and enriched material into their surroundings, and die spectacularly as supernovae, leaving behind compact remnants and heavy elements (such as those that make life on Earth possible). Despite their relative rarity, they play a disproportionate role in the evolution of the galaxies that host them, and likely also played a significant role in the early days of the Universe. A subset ($\sim$10\%) of these stars was also found to host magnetic fields on their surface. These fields impact their evolution, and may lead to exotic physics (e.g., heavy stellar-mass black holes, pair-instability supernovae, magnetars, etc.). However, the detection and measurement of magnetic fields is limited, due to current instrumentation, to nearby massive stars in the Milky Way. To truly understand how magnetism arises in massive stars, and what role it might have played in earlier stages of our Universe, we require next-generation hardware, such as the proposed near-infrared-to-ultraviolet spectropolarimeter Pollux, on the Habitable Worlds Observatory (HWO). In this contribution, we detail how Pollux @ HWO will enable new frontiers in the study of magnetic massive stars, delivering results that will profoundly impact the fields of stellar formation, stellar evolution, compact objects, and stellar feedback.

We present a comprehensive analysis of the HL Tau dust disk by modeling its intensity profiles across six wavelengths (0.45 to 7.9 mm) with a resolution of 0.05 arcsec ($\sim7$ au). Using a Markov Chain Monte Carlo (MCMC) approach, we constrain key dust properties including temperature, surface density, maximum grain size, composition, filling factor, and size distribution. The full fitting, with all parameters free, shows a preference for organics-rich dust with a low filling factor in the outer region ($r \gtrsim 40$ au), where the spectral index is $\sim3.7$, but amorphous-carbon-rich dust also reasonably reproduces the observed intensity profiles. Considering the scattering polarization observed at 0.87 mm, compact, amorphous-carbon-rich dust is unlikely, and moderately porous dust is favored. Beyond 40 au, the maximum dust size is likely $\sim100~{\rm \mu m}$ if dust is compact or amorphous-carbon rich. However, if the dust is moderately porous and organics-rich, both the predicted dust surface density and dust size can be sufficiently large for the pebble accretion rate to reach $\sim10M_{\oplus}~{\rm Myr^{-1}}$ in most regions, suggesting that pebble accretion could be a key mechanism for forming planets in the disk. In contrast, if the dust is amorphous-carbon-rich, forming a giant planet core via pebble accretion is unlikely due to the combined effects of low dust surface density and small dust size required to match the observed emission, suggesting other mechanisms, such as disk fragmentation due to gravitational instability, may be responsible for planet formation in the HL Tau disk.

We report the discovery of several new pulses from the source ASKAP J175534.9-252749.1 (J1755-2527), originally identified from a single 2-min long pulse, confirming it as a long period transient (LPT) with a period of ~1.16 hours. The pulses are significantly scattered, consistent with Galactic electron density models. Two of the new pulses also had measurable polarisation, but unlike the originally detected pulse, the polarisation angle does not behave as expected from the rotating vector model. We interpret historical non-detections of J1755-2527 as an intrinsic intermittency that occurs on month-long timescales, and discuss possible causes. We conjecture that, like some other LPTs with periods >~ 1 hour, J1755-2527 may host a white dwarf in a binary orbit, but note that its period is marginally shorter than the canonical orbital period minimum of cataclysmic variables. Our work highlights the importance of additional observations in establishing the nature of unusual radio-emitting objects.

Carbon-rich Wolf-Rayet (WR) stars are significant contributors of carbonaceous dust to the galactic environment, however the mechanisms and conditions for formation and subsequent evolution of dust around these stars remain open questions. Here we present JWST observations of the WR+WR colliding-wind binary Apep which reveal an intricate series of nested concentric dust shells that are abundant in detailed substructure. The striking regularity in these substructures between successive shells suggests an exactly repeating formation mechanism combined with a highly stable outflow that maintains a consistent morphology even after reaching 0.6 pc (assuming a distance of 2.4 kpc) into the interstellar medium. The concentric dust shells show subtle deviations from spherical outflow, which could reflect orbital modulation along the eccentric binary orbit or some mild degree of non-sphericity in the stellar wind. Tracking the evolution of dust across the multi-tiered structure, we measure the dust temperature evolution that can broadly be described with an amorphous carbon composition in radiative thermal equilibrium with the central stars. The temperature profile and orbital period place new distance constraints that support Apep being at a greater distance than previously estimated, reducing the line-of-sight and sky-plane wind speed discrepancy previously thought to characterise the system.

ALMA-IMF is a large program of the Atacama Large Millimeter/submillimeter Array (ALMA) that aims to determine the origin of the core mass function (CMF) of 15 massive Galactic protoclusters (~ 1.0-25.0 x 10^3 solar mass within ~ 2.5 x 2.5 pc^2 ) located towards the Galactic plane. In addition, the objective of the program is to obtain a thorough understanding of their physical and kinematic properties. Here we study the turbulence in these protoclusters with C18O (2-1) emission line using the sonic Mach number analysis (M_s ) and the size-linewidth relation. The probability distribution functions (PDFs) for M_s show a similar pattern, exhibiting no clear trend associated with evolutionary stage, peaking in the range between 4 and 7, and then extending to ~ 25. Such values of M_s indicate that the turbulence in the density regime traced by the C18O line inside the protoclusters is supersonic in nature. In addition, we compare the non-thermal velocity dispersions (sigma_nth, C18O) obtained from the C18O(2-1) line with the non-thermal line widths (sigma_nth, DCN ) of the cores obtained from the DCN (3-2) line. We observe that, on average, the non-thermal linewidth in cores is half that of the gas surrounding them. This suggests that turbulence diminishes at smaller scales or dissipates at the periphery of the cores. Furthermore, we examine the size-linewidth relation for the structures we extracted from the position-position-velocity C18O(2-1) line emission cube with dendrogram algorithm. The power-law index (p) obtained from the size-linewidth relation is between 0.41 and 0.64, steeper than the Kolmogorov law of turbulence, as expected for compressible media. In conclusion, this work is one of the first to carry out such a statistical study of turbulence for embedded massive protoclusters.

Human activities degrade the Earth environment at an unprecedented scale and pace, threatening Earth-system stability, resilience and life-support functions. We can of course deny the facts, get angry about them, try to bargain, or fall into deep depression. Or we may overcome these stages of grief and move towards accepting that human activities need to change, including our own ones. The purpose of this paper is to support astronomers in this transition, by providing insights into the origins of environmental impacts in astronomical research and proposing changes that would make the field sustainable. The paper focuses on the environmental impacts of research infrastructures, since these are the dominant sources of greenhouse gas emissions in astronomy, acknowledging that impact reductions in other areas, for example professional air travelling, need also to be achieved.

Atmospheric muons produced in cosmic-ray air showers are classified as conventional muons from pion and kaon decays and prompt muons from heavy hadron decays. Conventional muons dominate at lower energies, and the prompt component becomes dominant at PeV energies and above. Precisely measuring the atmospheric muon flux from a few GeV to several PeV is valuable for advancing our understanding of cosmic-ray interactions and testing hadronic interaction models. Low-energy muons that stop within the IceCube in-ice array provide valuable information about the energy spectrum of muons from a few hundred GeV up to 10 TeV. Machine learning techniques are employed to enhance event reconstruction and selection to provide insights into the conventional and prompt components. This contribution presents the unfolding of the energy spectrum of stopping muons in IceCube as well as the unfolding of high-energy muons to probe the prompt component.

Observations of Milky Way stars by multiplexed spectroscopic instruments and of gas in nearby galaxies using integral field units have made it possible to measure the abundances of multiple elements in both the interstellar medium and the stars that form out of it. These observations have revealed complex correlations between elemental abundances, but thus far there has been no theoretical framework to interpret these data. In this paper we extend the simple stochastically-forced diffusion model of Krumholz & Ting (2018), which has proven successful at explaining the spatial abundance patterns of single elements, to multiple elements, clarifying why elements are correlated and what controls their degree of correlation, and making quantitative predictions for the degree of correlation in both gas and young stars. We show that our results are qualitatively consistent with observed patterns, and point out how application of this theory to measured correlations should enable determination of currently unknown parameters describing r-process nucleosynthesis.

Context. Molecular hydrogen and its cation H+2 are among the first species formed in the early Universe, and play a key role in the thermal and chemical evolution of the primordial gas. In molecular clouds, H+2 ions formed through ionization of H2 by particles react rapidly with H2 to form H+3 , triggering the formation of almost all detected interstellar molecules. Aims. We present a new set of cross sections and rate coefficients for state-to-state ro-vibrational transitions of the H+2 and HD+ ions, induced by low-energy electron collisions. Study includes the major electron-impact processes relevant for low-metallicity astrochemistry: inelastic and superelastic scattering, and dissociative recombination. Methods. The electron-induced processes involving H+2 and HD+ are treated using the multichannel quantum defect theory. Results. The newly calculated thermal rate coefficients show significant differences compared to those used in previous studies. When introduced into astrochemical models, particularly for shock-induced chemistry in metal-free gas, the updated dissociative recombination rates produce substantial changes in the predicted molecular abundances. Conclusions. These data provide updated and improved input for the modeling of hydrogen-rich plasmas in environments

The Calcium-K line profiles as functions of solar latitude and time were obtained through our observations from the Kodaikanal Solar Observatory using the solar tunnel telescope and spectrograph with a CCD detector. Observations were conducted on all days with favourable sky conditions. We analysed the data collected over a period of about ten years to study the variations in the Ca II K line profiles recorded between 2015 and 2024, of which 709 days of data were found useful. The temporal and time-averaged latitudinal variations of the K$_{1}$ width, K$_{2}$ width, K$_{3}$ intensity and the intensity ratios of K$_{2v}$/K$_{2r}$ and K$_{2v}$/K$_{3}$ were computed using a semi-automated program. The parameters showed asymmetric increases towards the higher latitudes, with the rates of increase being higher in the southern hemisphere. The temporal plots for K$_{1}$ width and K$_{3}$ intensity showed positive correlations with the plage and spot filling factors, whereas the temporal plots for K$_{2}$ width, K$_{2v}$/K$_{2r}$ and K$_{2v}$/K$_{3}$ intensity ratios showed negative correlations. The time-averaged latitudinal plot for K$_{1}$ width has small peaks near 25°N and 20°S. The K$_{2}$ width has a small peak at 0°. The K$_{3}$ intensity has peaks at 20°N and 15°S. The K$_{2v}$/K$_{2r}$ intensity ratio shows peaks at 50$\degree$N, 0° and 40°S. The K$_{2v}$/K$_{3}$ intensity ratio shows peaks at 60°N, 0° and 60°S. Slope profiles show spectral response to magnetic activity peaks near K$_{3}$ with north-south asymmetries. Such variations in the line profiles are important in the studies of solar irradiance, surface flux transport and solar dynamo.

Much of the carbonaceous dust observed in the early universe may originate from colliding wind binaries (CWBs) hosting hot, luminous Wolf-Rayet (WR) stars. Downstream of the shock between the stellar winds there exists a suitable environment for dust grain formation, and the orbital motion of the stars wraps this dust into richly structured spiral geometries. The Apep system is the most extreme WR-CWB in our Milky Way: two WR stars produce a complex spiral dust nebula, whose slow expansion has been linked to a gamma-ray burst progenitor. It has been unclear whether the O-type supergiant 0.7" distant from the WR+WR binary is physically associated with the system, and whether it affects the dusty nebula. Multi-epoch VLT/VISIR and JWST/MIRI observations show that this northern companion star routinely carves a cavity in the dust nebula - the first time such an effect has been observed in a CWB - which unambiguously associates the O star as a bound component to the Apep system. These observations are used together with a new geometric model to infer the cavity geometry and the orbit of the WR+WR binary, yielding the first strong constraints on wind and orbital parameters. We confirm an orbital period of over 190 years for the inner binary - nearly an order of magnitude longer than the next longest period dust-producing WR-CWB. This, together with the confirmed classification as a hierarchical triple, cements Apep as a singular astrophysical laboratory for studying colliding winds and the terminal life stages of the most massive star systems.

We present a comprehensive analysis of 32 type II supernovae (SNe II) with plateau phase photometry and late phase ($nebular$) spectroscopy available, aiming to bridge the gap between the surface and core of their red supergiant (RSG) progenitors. Using \texttt{MESA}\,+\texttt{STELLA}, we compute an extensive grid of SN II light curve models originating from RSG with effective temperatures $T_{\rm eff}$ around 3650\,K and hydrogen-rich envelopes artificially stripped to varying degrees. These models are then used to derive the hydrogen-rich envelope masses $M_{\rm Henv}$ for SNe II from their plateau phase light curves. Nebular spectroscopy further constrains the progenitor RSG's luminosity log\,$L_{\rm prog}$, and is employed to remove the degeneracies in light curve modeling. The comparison between log\,$L_{\rm prog}$-$M_{\rm Henv}$ reveals that $M_{\rm Henv}$ spans a broad range at the same log\,$L_{\rm prog}$, and almost all SNe II have lower $M_{\rm Henv}$ than the prediction of the default stellar wind models. We explore alternative wind prescriptions, binary evolution models, and the possibility of more compact RSG progenitors. Although binary interaction offers a compelling explanation for the non-monotonicity and large scatter in the log\,$L_{\rm prog}$-$M_{\rm Henv}$ relation, the high occurrence rate of partially-stripped RSGs cannot be accounted for by stable binary mass transfer alone without fine-tuned orbital parameters. This highlights that, despite being the most commonly observed class of core-collapse SNe, SNe II likely originate from a variety of mass-loss histories and evolutionary pathways that are more diverse and complex than typically assumed in standard stellar evolution models.

Dwarf galaxies dominate the galaxy population across all redshifts, and majority of mergers are expected to occur between them. However, the impact of dwarf-dwarf mergers on star formation (SF) is less understood. In this context, we study SF in 6,155 isolated dwarf galaxies (no massive galaxy, $M_{*} > 10^{10} M_{\odot}$, within 1 Mpc$^3$), including 194 post-mergers and 5,961 non-interacting galaxies, spanning stellar masses of $10^7$ - $10^{9.6}$ M$_{\odot}$ and redshifts of 0.01 - 0.12. Post-mergers sample is from a previous study, which used deep optical images to identify dwarf galaxies with merger signatures. Using GALEX FUV data, we estimate star formation rates (SFRs) and find difference in log(SFR) between a post-merger galaxy and the median of its corresponding control sample, matched in stellar mass and redshift. SFR offsets range from -2 to +2 dex, indicating both enhancement and suppression of SF in our sample with 67% of post-mergers showing enhancement. The median SFR is found to be elevated by 0.24 dex ($\sim$1.73 times) in post-mergers, comparable to enhancements seen in massive galaxies. SF is found to be similarly enhanced in both the central (6" diameter region) and outer regions of post-mergers with respect to their non-interacting counterparts. This is in contrast to what is observed in massive galaxies, where the merger-triggered SF is more significant in the central regions. In the given small range of redshift, post-merger dwarfs exhibit a higher median specific SFR compared to their non-interacting counterparts. About 33% of post-mergers are quenched, possibly reflecting a later stage of the post-merger regime, where quenching can happen as observed in massive galaxies. This study suggests that dwarf-dwarf mergers can affect SF in the local universe. A more comprehensive study of post-merger dwarfs is required to understand their evolution.

The dawning of the Space Age marked the start of an ongoing relationship between the professional astronomical community and both state and non-state actors that launch and operate spacecraft in near-Earth orbital space. While the Cold War heated up in the late 1950s, military uses of outer space quickly came into conflict with the priorities of astronomers then building ever-bigger ground-based telescopes and envisioning the first generation of space telescopes. As the threat of global thermonuclear war loomed, the United States carried out Project West Ford, which tested an 'artificial ionosphere' for microwave radio propagation by placing several hundred million tiny copper dipoles into a belt orbiting the Earth. While the test was ultimately successful, it ignited a firestorm of concern and criticism among astronomers and ultimately influenced the framing of the United Nations Outer Space Treaty. Here we examine the history of Project West Ford as it prompted astronomers to react, comparing it with the ongoing problem of the potential impact of large satellite constellations on astronomical research.

Searching for signs of life is a primary goal of the Habitable Worlds Observatory (HWO). However, merely detecting oxygen, methane, or other widely discussed biosignatures is insufficient evidence for a biosphere. In parallel with biosignature detection, exoplanet life detection additionally requires characterization of the broader physicochemical context to evaluate planetary habitability and the plausibility that life could produce a particular biosignature in a given environment. Life detection further requires that we can confidently rule out photochemical or geological phenomena that can mimic life (i.e. "false positives"). Evaluating false positive scenarios may require different observatory specifications than biosignature detection surveys. Here, we explore the coronagraph requirements for assessing habitability and ruling out known false positive (and false negative) scenarios for oxygen and methane, the two most widely discussed biosignatures for Earth-like exoplanets. We find that broad wavelength coverage ranging from the near UV (0.26 $\mu$m) and extending into the near infrared (1.7 $\mu$m), is necessary for contextualizing biosignatures with HWO. The short wavelength cutoff is driven by the need to identify Proterozoic-like biospheres via O$_3$, whereas the long wavelength cutoff is driven by the need to contextualize O$_2$ and CH$_4$ biosignatures via constraints on C-bearing atmospheric species. The ability to obtain spectra with signal-to-noise ratios of 20-40 across this 0.26-1.7 $\mu$m range (assuming R=7 UV, R=140 VIS, and R=70 NIR) is also required. Without sufficiently broad wavelength coverage, we risk being unprepared to interpret biosignature detections and may ultimately be ill-equipped to confirm the detection of an Earth-like biosphere, which is a driving motivation of HWO..

A primary objective in contemporary low background physics is the search for rare and novel phenomena beyond the Standard Model of particle physics, e.g. the scattering off of a potential Dark Matter particle or the neutrinoless double beta decay. The success of such searches depends on a reliable background prediction via Monte Carlo simulations. A widely used toolkit to construct these simulations is Geant4, which offers the user a wide choice of how to model the physics of particle interactions. For example, for electromagnetic interactions, Geant4 provides pre-defined sets of models: physics constructors. As decay products of radioactive contaminants contribute to the background mainly via electromagnetic interactions, the physics constructor used in a Geant4 simulation may have an impact on the total energy deposition inside the detector target. To facilitate the selection of physics constructors for simulations of experiments that are using CaWO$_\mathrm{4}$ and Ge targets, we quantify their impact on the total energy deposition for several test cases. These cases consist of radioactive contaminants commonly encountered, covering energy depositions via $\alpha$, $\beta$, and $\gamma$ particles, as well as two examples for the target thickness: thin and bulky. We also consider the computing performance of the studied physics constructors.

This work describes new dynamical simulations of terrestrial planet formation. The simulations started at the protoplanetary disk stage, when planetesimals formed and accreted into protoplanets, and continued past the late stage of giant impacts. We explored the effect of different parameters, such as the initial radial distribution of planetesimals and Type-I migration of protoplanets, on the final results. In each case, a thousand simulations were completed to characterize the stochastic nature of the accretion process. In the model best able to satisfy various constraints, Mercury, Venus, and Earth accreted from planetesimals that formed early near the silicate sublimation line near 0.5 au and migrated by disk torques. For Venus and Earth to end up at 0.7-1 au, Type-I migration had to be directed outward, for example as the magnetically driven winds reduced the surface gas density in the inner part of the disk. Mercury was left behind near the original ring location. We suggest that Mars and multiple Mars-sized protoplanets grew from a distinct outer source of planetesimals at 1.5-2 au. While many migrated inwards to accrete onto the proto-Earth, our Mars was the lone survivor. This model explains: (1) the masses and orbits of the terrestrial planets, (2) the chemical composition of the Earth, where ~70% and ~30% come from reduced inner-ring and more-oxidized outer-ring materials, and (3) the isotopic differences of the Earth and Mars. It suggests that the Moon-forming impactor Theia plausibly shared a similar isotopic composition and accretion history with that of the proto-Earth.

The landmark detection of both gravitational waves (GWs) and electromagnetic (EM) radiation from the binary neutron star merger GW170817 has spurred efforts to streamline the follow-up of GW alerts in current and future observing runs of ground-based GW detectors. Within this context, the radio band of the EM spectrum presents unique challenges. Sensitive radio facilities capable of detecting the faint radio afterglow seen in GW170817, and with sufficient angular resolution, have small fields of view compared to typical GW localization areas. Additionally, theoretical models predict that the radio emission from binary neutron star mergers can evolve over weeks to years, necessitating long-term monitoring to probe the physics of the various post-merger ejecta components. These constraints, combined with limited radio observing resources, make the development of more coordinated follow-up strategies essential -- especially as the next generation of GW detectors promise a dramatic increase in detection rates. Here, we present RADAR, a framework designed to address these challenges by promoting community-driven information sharing, federated data analysis, and system resilience, while integrating AI methods for both GW signal identification and radio data aggregation. We show that it is possible to preserve data rights while sharing models that can help design and/or update follow-up strategies. We demonstrate our approach through a case study of GW170817, and discuss future directions for refinement and broader application.

Moreton waves are widely regarded as the chromospheric counterpart of extreme ultraviolet (EUV) waves propagating in the corona. However, direct observational evidence confirming their simultaneous propagation across multiple atmospheric layers from the corona through the transition region to the chromosphere has been lacking. In this study, we present comprehensive observational evidence of a three-dimensional (3D) fast-mode wave propagating from the corona through the transition region into the chromosphere, exhibiting a gradual deceleration. Additionally, this wave interacts with three filaments (F1, F2, and F3) along its path, inducing oscillation with multiple amplitudes: Filaments F1 and F2 exhibit simultaneous horizontal and vertical large-scale oscillations ($\sim$\speed{20}), while Filament F3 only exhibits vertical small-scale oscillation ($\sim$\speed{4}). Interestingly, F1 displays a similar oscillation period of about 500\,s in both horizontal and vertical directions, whereas F2 shows significantly different periods in these two dimensions (1100\,s and 750\,s), and F3 exhibits only a vertical oscillation with a period of about 450\,s. Based on this kinematic behavior, we propose that their oscillations were likely triggered by compression from the flanks of the dome-shaped wavefront. We further estimate the magnetic fields of the filaments. The radial (axial) magnetic fields for F1 and F2 are estimated to be 14.9\,G (28.6\,G) and 9.9\,G (18.6\,G), respectively. For F3, we estimate its radial magnetic field to be 16.6\,G.

The increasing field of view of radio telescopes and improved data processing capabilities have led to a surge in the detection of Fast Radio Bursts (FRBs). The discovery rate of FRBs is already a few per day and is expected to increase rapidly with new surveys coming online. The growing number of events necessitates prioritized follow-up due to limited multi-wavelength resources, requiring rapid and automated classification. In this study, we introduce Frabjous, a deep learning framework for an automated morphology classifier with an aim towards enabling the prompt follow-up of anomalous and intriguing FRBs, and a comprehensive statistical analysis of FRB morphologies. Deep learning models require a large training set of each FRB archetype, however, publicly available data lacks sufficient samples for most FRB types. In this paper, we build a simulation framework for generating realistic examples of FRBs and train a network based on a combination of simulated and real data, starting with the CHIME/FRB catalog. Applying our framework to the first CHIME/FRB catalog, we achieve an overall classification accuracy of approximately 55%, well over a random multiclass classification rate of 20 % with five balanced classes during training. While this falls short of desirable performance, we critically discuss the limitations of our approach and propose potential avenues for improvement. Future work should explore strategies to augment training datasets and broaden the scope of FRB morphological studies, aiming for more accurate and reliable classification results.

Planets orbiting young, solar-type stars are embedded in a more energetic environment than that of the solar neighbourhood. They experience harsher conditions due to enhanced stellar magnetic activity and wind shaping the secular evolution of a planetary atmosphere. This study is dedicated to the characterisation of the magnetic activity of eleven Sun-like stars, with ages between 0.2 and 6.1 Gyr and rotation periods between 4.6 and 28.7 d. Based on a sub-sample of six stars, we aim to study the large-scale magnetic field, which we then use to simulate the associated stellar wind and environment. Finally, we want to determine the conditions during the early evolution of planetary habitability. We analysed high-resolution spectropolarimetric data collected in 2018 and 2019 with Narval. We computed activity diagnostics from chromospheric lines such as CaII H&K, H$\alpha$, and the CaII infrared triplet, as well as the longitudinal magnetic field from circularly polarised least-squares deconvolution profiles. For six stars exhibiting detectable circular polarisation signals, we reconstructed the large-scale magnetic field at the photospheric level by means of Zeeman-Doppler imaging (ZDI). In agreement with previous studies, we found a global decrease in the activity indices and longitudinal field with increasing age and rotation period. The large-scale magnetic field of the six sub-sample stars displays a strength between 1 and 25 G and reveals substantial contributions from different components such as poloidal (40-90 %), toroidal (10-60 %), dipolar (30-80 %), and quadrupolar (10-40 %), with distinct levels of axisymmetry (6-84 %) and short-term variability of the order of months. Ultimately, this implies that exoplanets tend to experience a broad variety of stellar magnetic environments after their formation.

Ultraluminous and hyperluminous X-ray sources (ULXs and HLXs) are among the brightest astrophysical objects in the X-ray sky. While ULXs most likely host stellar-mass compact objects accreting at super-Eddington rates, HLXs are compelling candidates for accreting intermediate-mass black holes. Our goal is to produce a clean sample of HLXs by removing possible contaminants and characterise the spectral properties of the remaining population. Starting with a set of 115 HLXs detected by XMM-Newton, we identified and removed contaminants (AGNs, X-ray diffuse emission detected as point-like, and tidal disruption event candidates) and retrieved 40 sources for which XMM-Newton spectra are available. We fitted them with an absorbed power law model and determined their unabsorbed luminosities and hardness ratios. We constructed the hardness-luminosity diagram, compared the results with the spectral properties of the HLX prototype, ESO 243-49 HLX-1, and conducted a deeper analysis on a few promising candidates. The resulting HLX population spans a luminosity range from $1\times10^{41}$ erg s$^{-1}$ to nearly $10^{43}$ erg s$^{-1}$ and is homogeneously spread in hardness between 0.5 and 5. Half of the population has hardness ratios higher than a typical AGN, and could be considered the extension of the ULX population at higher energies. We found four very soft outliers, which are characterised by steep power law spectra and no X-ray emission above 1$-$2 keV, similarly to ESO 243-49 HLX-1. Those with multi-epoch archival data show changes in luminosity up to almost two orders of magnitudes. We show that sources currently identified as HLXs can be more diverse than ULXs and disentangling between different types of objects is not trivial with currently available data. New observations would be beneficial to expand the current sample and uncover the true nature of many objects of this class.

Recent 3D dust maps of the local Milky Way are revolutionizing our understanding of the Sun's Galactic neighborhood, bringing much-needed insight into the large-scale organization of the interstellar medium. Focusing on the largest scales in $\textit{Gaia}$-based 3D dust maps, we find a pattern of seven highly elongated, mostly parallel structures in the local $5\,\mathrm{kpc}^2$, five of which were previously unknown. These structures show pitch angles of $33.5 \pm 4.0 ^\circ$, and masses ranging from $10^5$ to $10^6$ $M_\odot$. We refer to these structures as superclouds. Nearly all known star-forming regions in the solar neighborhood lie within the superclouds, primarily along their central axes, supporting the idea that they are precursors to giant molecular clouds. All but one of the seven superclouds show an underlying undulation, indicating that this is not a property unique to the Radcliffe wave. We find that while the superclouds have line masses that vary by about a factor of four, their volume densities vary by about 10$\%$ only. This suggests that superclouds self-regulate their physical sizes and internal structure to maintain pressure equilibrium with their environment. These findings establish a new framework for understanding how large-scale Galactic structures shape the conditions for star formation in the solar vicinity, and likely in galaxies like the Milky Way.

Ethylene glycol ($\mathrm{(CH_2OH)_2}$, hereafter EG) and Glycolonitrile ($\mathrm{HOCH_2CN}$, hereafter GN) are considered molecular precursors of nucleic acids. EG is a sugar alcohol and the reduced form of Glycolaldehyde ($\mathrm{CH_2(OH)CHO}$, hereafter GA). GN is considered a key precursor of adenine formation (nucleotide) and can be a precursor of glycine (amino acid). Detections of such prebiotic molecules in the interstellar medium are increasingly common. How much of this complexity endures to the planet formation stage, and thus is already present when planets form, remains largely unknown. Here we report Atacama Large Millimeter/sub-millimeter Array (ALMA) observations in which we tentatively detect EG and GN in the protoplanetary disk around the outbursting protostar V883 Ori. The observed EG emission is best reproduced by a column density of $\mathrm{3.63^{+0.11}_{-0.12} \times 10^{16} \; cm^{-2}}$ and a temperature of at least 300 K. The observed GN emission is best reproduced by a column density of $\mathrm{3.37^{+0.09}_{-0.09} \times 10^{16} \; cm^{-2}}$ and a temperature of $88^{+1.2}_{-1.2}$ K. Comparing the abundance of EG and GN relative to methanol in V883 Ori with other objects, V883 Ori falls between hot cores and comets in terms of increasing complexity. This suggests that the build up of prebiotic molecules continues past the hot core phase into the epoch of planet formation. Nascent planets in such environments may inherit essential building blocks for life, enhancing their potential habitability. Further observations of this protoplanetary disk at higher spectral resolution are required to resolve blended lines and to confirm these tentative detection.

We present optical and near-infrared spectroscopy of the interstellar object 3I/ATLAS, obtained with Gemini-S/GMOS and NASA IRTF/SpeX on 2025 July 5 and July 14. The optical spectrum shows a red slope of approximately 10% per 1000 angstroms between 0.5 and 0.8 microns, closely resembling that of typical D-type asteroids. At longer wavelengths, the near-infrared spectrum flattens significantly to approximately 3% per 1000 angstroms from 0.9 to 1.5 microns, consistent with the spectral behavior of large water ice grains in the coma. Spectral modeling with an areal mixture of 70% Tagish Lake meteorite and 30% 10-micron-sized water ice successfully reproduces both the overall continuum and the broad absorption feature near 2.0 microns. The 1.5-micron water ice band, however, is not detected, likely due to the limited signal-to-noise of the IRTF data and dilution by refractory materials. The ~30% ice fraction should be interpreted as an approximate, order-of-magnitude estimate of the coma composition. The agreement between the GMOS and SpeX spectra, taken nine days apart, suggests short-term stability in the coma's optical properties. Our observations reveal that 3I/ATLAS is an active interstellar comet containing abundant water ice, with a dust composition more similar to D-type asteroids than to ultrared trans-Neptunian objects.

Despite growing efforts to find the sources of high energy neutrinos measured by IceCube, the bulk of the neutrinos remain with unknown origins. We aim to constrain the emissivity of cosmic high-energy neutrinos from extragalactic sources through their correlation with the large-scale structure. We use cross-correlations between the IceCube 10-year dataset and tomographic maps of the galaxy overdensity to place constraints on the bias-weighted high-energy neutrino emissivity out to redshift $z\sim3$. We test two different models to describe the evolution of neutrino emissivity with redshift, a power law model $\propto (1+z)^a$, and a model tracking the star formation history, assuming a simple power law model for the energy injection spectrum. We also consider a non-parametric reconstruction of the astrophysical neutrino emissivity as a function of redshift. We do not find any significant correlation, with our strongest results corresponding to a $1.9 \sigma$ deviation with respect to a model with zero signal. We use our measurements to place upper bounds on the bias-weighted astrophysical high-energy neutrino emission rate as a function of redshift for different source models. This analysis provides a new probe to test extragalactic neutrino source models. With future neutrino and galaxy datasets we expect the constraining and detection power of this type of analysis analysis to increase.

A spherical system with arbitrary mean velocities in the radial and transverse directions, but with zero dispersion of radial velocities, is considered as a model describing the violent relaxation of star clusters and galaxies. The model allows to reduce the infinite chain of Jeans' equations of moments to just three hydrodynamic equations. The found system of equations is consistent with the law of conservation of energy and the virial equation. This part of the work is devoted to a general description of the model, then the calculation results will be presented.

A majority of JWST/NIRSpec/IFU studies at high redshifts to date have focused on UV-bright or massive objects, while our understanding of low-mass galaxies at early cosmic times remains limited. In this work, we present NIRSpec/IFS high-resolution observations of two low-mass ($M_* < 10^9 \ M_\odot$), low-metallicity ($[12 + \log(\text{O/H})] < 8$) galaxies at $z \sim 7.66$ with evidence of hosting AGN. Using spatially-resolved maps of rest-frame optical emission lines, we find flat metallicity profiles, indicative of ISM redistribution by outflows or past merging. We identify kinematical components decoupled from galactic rotation with velocities of $\sim 250 - 500 \ \text{km} \ \text{s}^{-1}$. We argue that these components are likely tracing outflows, possibly AGN-driven, for which we infer outflow rates of $\sim 21 - 40 \ M_\odot \ \text{yr}^{-1}$, suggesting they may suppress future star formation. We compare our observational results to those from the new large-volume AESOPICA simulations, which fully incorporate different models of black hole growth and AGN feedback. We find that our observational results of $v_\text{out}/v_\text{esc}$ and $\dot{M}_\text{out}$/SFR are consistent with the AGN population in these simulations, hinting that AGN-driven feedback may contribute to quenching both in our systems and in a wider population of low-mass galaxies in the early Universe. This novel study demonstrates the necessity of deep IFU observations to decompose the complex kinematics and morphology of high-$z$ galaxies, trace outflows, and constrain the effect of feedback in these low-mass systems.

Gamma-ray bursts observed in high-energies allow the investigation of the emission processes of these still puzzling events. In this study, we perform general relativistic magnetohydrodynamic (GRMHD) simulations to investigate GRB 090510, a peculiar short GRB detected by Fermi-LAT. Our primary goal is to model the energetics, jet structure, variability, and opening angle of the burst to understand its underlying physical conditions. We tested the 2D and 3D models and estimated the time scale of variability. The predicted energetics and the jet opening angle reconcile with the observed ones with 1$\sigma$ when considering that the jet opening angles also evolve with redshift. Furthermore, we extend our analysis by incorporating dynamical ejecta into selected models to study its impact on jet collimation at smaller distances. In addition, we investigated a suite of models exhibiting a broad range of observable GRB properties, thereby extending our understanding beyond this specific event.

While astrophysical observations imply that 85% of the matter content is unaccounted for, the nature of this dark matter (DM) component remains unknown. Weakly Interacting Massive Particles (WIMPs) - DM particles that interact at or below the weak interaction scale - could naturally explain this missing matter. These interactions with the Standard Model (SM) allow them to be gravitationally captured in celestial bodies like the Sun. Trapped DM in the solar core could subsequently annihilate, producing stable SM particles, of which only neutrinos can escape the Sun's dense interior. Therefore, an excess of neutrinos originating from the direction of the Sun would serve as evidence of DM. The IceCube Upgrade, a dense infill of the IceCube Neutrino Observatory, will lower the energy threshold and improve sensitivity in the range from 1 to 500 GeV, thereby enhancing IceCube's ability to detect GeV-scale DM. In this contribution, I present projections of the IceCube Upgrade's sensitivity to the DM-proton scattering cross section for DM masses between 3 GeV and 500 GeV. These sensitivities position IceCube as the most sensitive indirect detection experiment for DM in the mass range from 3 GeV to 10 TeV.

Aims. Particles ejected from the lunar surface via hypervelocity impacts form a torus between the Earth and the Moon. According to our previous study (Yang et al., A\&A, 659, A120), among them about $2.3\times10^{-4}\,\mathrm{kg/s}$ particles impact the Earth after long-term orbital evolution. We mainly focus on these Earth impactors, analyze their orbital element distribution, and estimate their influence on Earth-based observations. Methods. In previous work we simulated the long-term orbital evolution of particles ejected from the lunar surface, and obtained their steady-state spatial distribution in the Earth-Moon system. In this work, we analyze the simulation results about the Earth impactors, including the fraction of impactors with different initial parameters among all impactors, the orbital element distribution, and the projection of particles onto several Earth-based observatories. Results. Particles ejected from the lunar surface are more likely to impact the Earth within a certain range of initial parameters. Most of these lunar-ejected impactors ($\sim70\%$) reach the Earth within one year, while most of the small ones ($87.2\%$ of $0.2\,\mathrm{\mu m}$ particles and $64.6\%$ of $0.5\,\mathrm{\mu m}$ particles) reach the Earth within one week. A large proportion of lunar-ejected Earth impactors can be distinguished from interplanetary dust particles according to the differences in their orbital distributions. Besides, lunar-ejected particles may exhibit distinct configurations and orientations from the perspectives of different Earth-based observatories.

We present an in-depth photometric and spectroscopic study of the contact binary KIC 7766185. Spectroscopic observations were conducted on the Mayall 4-m telescope at Kitt Peak National Observatory and used to extract radial velocities. Using the radial velocity measurements, Kepler photometry, and Gaia multi-color photometry, binary analysis was performed in PHOEBE to determine orbital and stellar parameters. The results classify KIC 7766185 as an A-type W UMa system and reveal one of the largest and most massive secondary stars in the sample of well studied W UMa systems. The two stars are in a shallow contact state, with a small fillout factor of $0.029\pm0.002$. Along with the binary analysis, we conduct a period study that reveals evidence for a cyclic variation in the eclipse timings.

The search for habitable, Earth-like exoplanets faces major observational challenges due to their small size and faint signals. M-dwarf stars offer a promising avenue to detect and study smaller planets, especially sub-Neptunes-among the most common exoplanet types. K2-18 b, a temperate sub-Neptune in an M-dwarf habitable zone, has been observed with HST and JWST, revealing an H2-rich atmosphere with CH4 and possible CO2. Conflicting interpretations highlight the importance of non-equilibrium chemistry, which is critical for constraining atmospheric parameters like metallicity, C/O ratio, and vertical mixing (Kzz). This study explores the parameter space of metallicity, C/O ratio, and Kzz for K2-18 b using the non-equilibrium chemistry model FRECKLL and JWST data. We generated spectra from a 3D grid of models and compared them to observations to refine atmospheric constraints. A fixed pressure-temperature profile was used to capture first-order chemical trends, acknowledging some uncertainties. Our best-fit model favors high metallicity (266^{+291}_{-104} at 2 sigma) and high C/O ratio (C/O > 2.1 at 2 sigma). CH4 is robustly detected (log10[CH4] = -0.3^{+0.1}_{-1.7} at 1 mbar), while CO2 remains uncertain due to spectral noise. Kzz has no clear impact on the fit and remains unconstrained. Non-equilibrium models outperform flat spectra at > 4 sigma confidence, confirming atmospheric features. Minor species, such as H2O and NH3, may be present but are likely masked by dominant absorbers. Our results highlight the limits of constant-abundance retrievals. The atmosphere has a high C/O ratio suggesting possible aerosol formation. Better constraints require higher-precision data. Future JWST NIRSpec G395H and ELT/ANDES observations will be critical for probing habitability and refining models.

Thin and thick discs are prominent components in the Milky Way and other spiral galaxies, but their formation histories are not yet understood, which is particularly true for thick discs. This is partially due to the fact that we lack sufficient understanding of thin and thick discs in other galaxies to determine how common features of the Milky Way's discs are. Here, we conduct an integral field spectroscopy study of the thin and thick discs in the edge-on Milky Way analogue IC 2531, using observational data from the WiFeS IFU spectrograph on the ANU 2.3m telescope. We provide spectral analysis of IC 2531's kinematics and stellar populations in the thin and thick disc regions by conducting pPXF fitting with Vazdekis models. We found that IC 2531's disc above the dust plane generally has chemical properties between the Milky Way's chemical thin and thick disc. IC 2531's thick disc is somewhat similar to the Milky Way's in terms of stellar populations: on average, it is older, more metal-poor than its corresponding thin disc, but similarly alpha-enhanced in general. But we do find a clearly alpha-rich thick disc bin though one of our "thick disc" bins may be dominated by a warped thin disc. These results may help to constrain formation theories of thick discs.

Blazars, including BL Lacs and FSRQs, are the most luminous extragalactic {\gamma}-ray sources. They account for about 56% of the sources listed in the recent Fermi-LAT catalog (4FGL-DR4). The optical and UV spectra of BL Lacs are nearly featureless, making it difficult to precisely determine their redshifts. Consequently, nearly half of the {\gamma}-ray BL Lacs lack reliable redshift measurements. This poses a major challenge, since redshift is crucial for studying the cosmic evolution of the blazar population and {\gamma}-ray propagation studies such as indirect evidence of EBL, placing constraints on IGMF and searches for LIV and ALPs. This paper is the fourth in a series dedicated to determining the redshift of a sample of blazars identified as key targets for future observations with the Cherenkov Telescope Array Observatory (CTAO). We performed Monte Carlo simulations to select {\gamma}-ray blazars detected by Fermi-LAT with hard spectra, that lack redshift measurements. These blazars are expected to be detectable by CTAO within 30 hours or less of exposure assuming an average flux state. In this fourth paper, we report the results of detailed spectroscopic observations of 29 blazars using the ESO/VLT, Keck II, and SALT telescopes. Our analysis involved a thorough search for spectral lines in the spectra of each blazar, and when features of the host galaxy were identified, we modeled its properties. We also compared the magnitudes of the targets during the observations to their long-term light curves. In the sample studied, 9 of 29 sources were observed with a high signal-to-noise ratio (S/N > 100), while the remaining 20 were observed with a moderate or low S/N. We successfully determined firm redshifts for 12 blazars, ranging from 0.1636 to 1.1427, and identified two lower limit redshifts at z > 1.0196 and z > 1.4454. The remaining 15 BL Lac objects exhibited featureless spectra.

It is projected that more than 100,000 communication satellites will be deployed in Low-Earth Orbit (LEO) over the next decade. These LEO satellites (LEOsats) will be captured frequently by the survey camera onboard the China Space Station Telescope (CSST), contaminating sources in the images. As such, it is necessary to assess the impact of LEOsats on CSST survey observations. We use the images taken by the Hubble Space Telescope (HST) in its F814W band to simulate $i$-band images for the CSST. The simulation results indicate that LEOsats at higher altitudes cause more contamination than those at lower altitudes. If 100,000 LEOsats are deployed at altitudes between 550 km and 1200 km with a 53-degree orbital inclination, the fraction of contaminated sources in a 150-s exposure image would remain below 0.50%. For slitless spectroscopic images, the contaminated area is expected to be below 1.50%. After removing the LEOsat trails, the residual photon noise contributes to relative photometric errors that exceed one-tenth of the total error budget in approximately 0.10% of all sources. Our investigation shows that even though LEOsats are unavoidable in CSST observations, they only have a minor impact on samples extracted from the CSST survey.

The WALOPControl software is designed to facilitate comprehensive control and operation of the WALOP (Wide Area Linear Optical Polarimeter) polarimeters, ensuring safe and concurrent management of various instrument components and functionalities. This software encompasses several critical requirements, including control of the filter wheel, calibration half-wave plate, calibration polarizer, guider positioning, focusers, and 4 concurrent CCD cameras. It also manages the host telescope and dome operations while logging operational parameters, user commands, and environmental conditions for troubleshooting and stability. It provides a user-friendly graphical user interface, secure access control, a notification system for errors, and a modular configuration for troubleshooting are integral to the software's architecture. It is accessible over the internet with the backend developed using NodeJS and ExpressJS, featuring a RESTful API that interacts with a MongoDB database, facilitating real-time status updates and data logging. The frontend utilizes the this http URL framework, with Redux for state management and Material UI for the graphical components. The system also allows for automatic observations based on user-defined schedules. A Continuous Integration and Continuous Deployment (CI CD) pipeline ensures the software's reliability through automated testing and streamlined deployment. The WALOPControl software is a key component of the PASIPHAE (Polar-Areas Stellar Imaging in Polarimetry High Accuracy Experiment) project, which aims to study the dust and magnetic field of the Milky Way by observing the polarization of starlight.

Detecting continuous gravitational waves (CWs) is challenging due to their weak amplitude and high computational demands, especially with poorly constrained source parameters. Stochastic gravitational-wave background (SGWB) searches using cross-correlation techniques can identify unresolved astrophysical sources, including CWs, at lower computational cost, albeit with reduced sensitivity. This motivates a hybrid approach where SGWB algorithms act as a first-pass filter to identify CW candidates for follow-up with dedicated CW pipelines. We evaluated the discovery potential of the SGWB analysis tool PyStoch for detecting CWs, using simulated signals from spinning down NSs. We then applied the method to data from the third LIGO-Virgo-KAGRA observing run (O3), covering the (20-1726) Hz frequency band, and targeting four supernova remnants: Vela Jr., G347.3-0.5, Cassiopeia A, and the NS associated with the 1987A supernova remnant. If necessary, significant candidates are followed up using the 5-vector Resampling and Band-Sampled Data Frequency-Hough techniques. However, since no interesting candidates were identified in the real O3 analysis, we set 95\% confidence-level upper limits on the CW strain amplitude $h_0$. The most stringent limit was obtained for Cassiopeia A, and is $h_0 = 1.13 \times 10^{-25}$ at $201.57$ Hz with a frequency resolution of $1/32$ Hz. As for the other targets, the best upper limits have been set with the same frequency resolution, and correspond to $h_0 = 1.20 \times 10^{-25} $ at $202.16$ Hz for G347.3-0.5, $1.20 \times 10^{-25}$ at $217.81$ Hz for Vela Jr., and $1.47 \times 10^{-25}$ at $186.41$ Hz for the NS in the 1987A supernova remnant.

The discovery of fast and variable coherent signals in a handful of ultraluminous X-ray sources (ULXs) testifies to the presence of super-Eddington accreting neutron stars, and drastically changed the understanding of the ULX class. Our capability of discovering pulsations in ULXs is limited, among others, by poor statistics. However, catalogues and archives of high-energy missions contain information which can be used to identify new candidate pulsating ULXs (PULXs). The goal of this research is to single out candidate PULXs among those ULXs which have not shown pulsations due to an unfavourable combination of factors. We applied an AI approach to an updated database of ULXs detected by XMM-Newton. We first used an unsupervised clustering algorithm to sort out sources with similar characteristics into two clusters. Then, the sample of known PULX observations has been used to set the separation threshold between the two clusters and to identify the one containing the new candidate PULXs. We found that only a few criteria are needed to assign the membership of an observation to one of the two clusters. The cluster of new candidate PULXs counts 85 unique sources for 355 observations, with $\sim$85% of these new candidates having multiple observations. A preliminary timing analysis found no new pulsations for these candidates. This work presents a sample of new candidate PULXs observed by XMM-Newton, the properties of which are similar (in a multi-dimensional phase space) to those of the known PULXs, despite the absence of pulsations in their light curves. While this result is a clear example of the predictive power of AI-based methods, it also highlights the need for high-statistics observational data to reveal coherent signals from the sources in this sample and thus validate the robustness of the approach.

A subclass of intermediate mass variables Delta Scuti stars, known as High-amplitude Delta Scuti (HADS) stars, exhibits pronounced radial pulsations with high amplitudes. The ground-based and space-based observations of the HADS star EH Lib are used to help making asteroseismological analysis of this pulsating star. Following the reduction of the light curves, the frequency analysis reveals the fundamental frequency as $f_0=11.3105$ c day$^{-1}$ and two more significant frequencies $f_1$ and $f_2$, in addition to the harmonics of $f_0$ and a linear combination. The period change rate is determined as $(1/P_0)(dP_0/dt)=(5.4\pm0.5)\times10^{-9}$ yr$^{-1}$ derived from an O-C diagram, which is constructed from 342 times of maximum light spanning over 70 years. Using these observational constraints, along with the metallicity reported in the literature, we construct theoretical models using the stellar evolution code MESA and calculate the theoretical frequencies of the eigen modes using the oscillation code GYRE. The appropriate models are selected by matching both $f_0$ and $(1/P_0)(dP_0/dt)$ within their respective uncertainties. The results indicate that the observed period change of EH Lib can be attributed to stellar evolutionary effects. The stellar parameters of EH Lib are derived as: the mass of $1.715\pm0.065$ M$_{\odot}$, the luminosity of log $(L/L_{\odot})=1.38\pm0.06$, and the age of $(1.14\pm0.13)\times10^{9}$ years. EH Lib is classified as a single-mode HADS star, locating currently in the Hertzsprung gap, with a helium core and a hydrogen-burning shell. This work expands the asteroseismological sample of HADS stars and establishes a foundation for future investigations into their commonalities and specific properties, thereby advancing our understanding of these variables.

With the same general purposes as Part I of this monograph, we analyze here major events in the history of the Earth, such as the formation of the Earth itself, the origin of life, the great glaciations and the mass extinctions of species, and we also analyze the astronomical context in which they occurred. We argue that the Sun was captured around 500-700 Myr ago by a supercloud, of which the Orion arm and Gould's belt are currently part, and that this event marked the history of the Earth. With this fact, we associate a massive capture of comets by the Sun, the great glaciations of the Cryogenic period, the emergence of complex life in the Cambrian and the recurrent mass extinctions of species.

We analyze in a systematic way the implications for the classes of hilltop and hilltop-squared inflation of the recent DR6 observations by the Atacama Cosmology Telescope collaboration. We find that the reported shift in the spectral index leads to parameter ranges for these models that are significantly reduced when compared to the results obtained from the {\sc Planck} observations. We mainly focus on the more dramatic implications for the hilltop-squared class, but along the way we also highlight the milder impact on the class of hilltop models.

We investigate the Type I migration of planets in low-density cavities and inner discs of strongly magnetized young stars using global three-dimensional (3D) magnetohydrodynamic (MHD) simulations, where the strong magnetic field carves the low-density cavity. Simulations show that a planet in the cavity migrates inwards up to the radius at which the outer Lindblad resonances are inside the cavity. At smaller radii, the migration stalls. The migration is faster if a star accretes in the unstable regime where the temporary tongues penetrate the magnetosphere. If a planet is in a highly inclined orbit, it interacts with the disc, and the eccentricity increases due to the Kozai-Lydov mechanism. A planet may stop or reverse its migration in the inner disc before entering the cavity. The magnetosphere interacts with the inner disc, changing its density distribution such that migration slows down or is even reversed. A tilted magnetosphere also excites density and bending waves in the disc, which may slow down migration and also increase the inclination and eccentricity of the planet. When a planet reaches the disc-cavity boundary, it is typically trapped at the boundary by asymmetric corotation torque. A planet moves together with the boundary when the cavity expands. Overall, a magnetized star provides an environment for slow or reverse migration.

Open clusters serve as important tools for accurately studying the chemical evolution of the Milky Way. By combining precise chemical data from high-resolution spectra with information on their distances and ages, we can effectively uncover the processes that have shaped our Galaxy. This study aims to derive NLTE atmospheric parameters and chemical abundances for approximately one hundred giant stars across 33 open clusters with near-solar metallicity. The clusters span a wide range of ages, enabling an assessment of the presence and extent of any age-related abundance gradients. In the Stellar Population Astrophysics (SPA) project, we acquired new high-resolution spectra of open clusters using the HARPS-N echelle spectrograph at the Telescopio Nazionale Galileo. We chemically characterized nine open clusters for the first time and reanalyzed previously studied SPA clusters, resulting in a consistent and homogeneous sample. We determined NLTE atmospheric parameters using the equivalent width method and derived NLTE chemical abundances through spectral synthesis for various elements, including alpha elements (Mg, Si, and Ti), light odd-Z elements (Na, Al), iron-peak elements (Mn, Co, and Ni), and neutron-capture elements (Sr, Y, and Eu). Our findings are compared with the existing literature, revealing good agreement. We examine the trends of [X/Fe] versus age, confirming previous observations and the enrichment patterns predicted by nucleosynthesis processes. Positive correlations with age are present for Mg, Si, Ti, Al, Mn, Co, Ni, and Sr, while Na and Y and Eu show a negative trend. This study emphasizes the significance of NLTE corrections and reinforces the utility of open clusters as tracers of Galactic chemical evolution. Furthermore, we provide a benchmark sample of NLTE abundances for upcoming open cluster surveys within large-scale projects such as 4MOST and WEAVE.

We analyze bars formed in $N$-body simulations to investigate two key aspects of stellar kinematic structure of barred galaxies: the angular distributions of the radial and azimuthal components of stellar velocities, and the impact of bars on rotation curves. We find that stars on bar-supporting $x_1$-like orbits exhibit characteristic sawtooth-like radial velocity patterns and arch-like tangential velocity patterns as a function of azimuth. In contrast, stars on box and disk orbits show little azimuthal variation, effectively smoothing the overall velocity distribution. When averaged over all orbital families, the resulting kinematics are broadly consistent with the bisymmetric model of Spekkens & Sellwood, with the amplitudes of bar-induced velocity perturbations increasing with bar strength. In addition, bars amplify the radial pressure gradient associated with enhanced random stellar motions, leading to a noticeable reduction in the mean rotational velocity. This effect becomes more pronounced with increasing bar strength, resulting in a shallower rotation curve within the bar region. We discuss our results in the context of the kinematic properties of observed barred galaxies.

We present an optimised multipole algorithm for computing the three-point correlation function (3PCF), tailored for application to large-scale cosmological datasets. The algorithm builds on a $in\, situ$ interpretation of correlation functions, wherein spatial displacements are implemented via translation window functions. In Fourier space, these translations correspond to plane waves, whose decomposition into spherical harmonics naturally leads to a multipole expansion framework for the 3PCF. To accelerate computation, we incorporate density field reconstruction within the framework of multiresolution analysis, enabling efficient summation using either grid-based or particle-based schemes. In addition to the shared computational cost of reconstructing the multipole-decomposed density fields - scaling as $\mathcal{O}(L^2_{\text{trun}} N_g \log N_g)$ (where $N_g$ is the number of grids and $L_{\text{trun}}$ is the truncation order of the multipole expansion) - the final summation step achieves a complexity of $\mathcal{O}(D^6_{\text{sup}} N_g)$ for the grid-based approach and $\mathcal{O}(D^3_{\text{sup}} N_p)$ for the particle-based scheme (where $D_{\text{sup}}$ is the support of the basis function and $N_p$ is the number of particles). The proposed $in\, situ$ multipole algorithm is fully GPU-accelerated and implemented in the open-source $Hermes$ toolkit for cosmic statistics. This development enables fast, scalable higher-order clustering analyses for large-volume datasets from current and upcoming cosmological surveys such as Euclid, DESI, LSST, and CSST.

Extending the investigation of the presumed primordial comet as part of continuing work on a new model of the Kreutz sungrazer system, I confront a previously derived set of orbital elements with Aristotle's remarks in his Meteorologica to test their compatibility and determine the comet's perihelion time. The two translations of the treatise into English that I am familiar with differ at one point substantially from each other. Unambiguously, the year and season of the comet's appearance was early 372 BC (or -371). From Aristotle's constraint on the comet's setting relative to sunset, I infer that the probable date of perihelion passage was January 20, a date also consistent with the vague remark on frosty weather. On the day that Aristotle claims the comet was not seen, its head may have been hidden behind the Sun's disk or in contact with it. The observation that the `comet receded as far as Orion's belt, where it dissolved' is being satisfied by the tested orbit if the perihelion was reached between January 20 and February 10. Aristotle's third statement, which describes the tail as a streak 60 degrees in length, suggests a plasma feature stretching in space over 0.8 AU. The dust tail was developing more gradually and it was all that could be seen from the comet when it was approaching Orion's belt in early April. The comet was seen over a period of more than 10 weeks. The results of this study strengthen the notion that Aristotle's comet indeed was the gigantic progenitor of Kreutz sungrazers.

The prompt $\gamma$-rays of gamma-ray bursts (GRBs) may originate from the photosphere of a relativistic jet. However, only a few GRBs have been observed with evident blackbody-like emission, for example, GRB~090902B. It has been demonstrated that internal dissipation processes, such as magnetic reconnection, can occur within the relativistic jet and thereby drive violent turbulence in the dissipation region. In this paper, we study the photospheric emission of a jet with turbulence below its photosphere. Here, the turbulence is modeled phenomenological under the assumption that the four-velocity of its eddies follows a Gaussian distribution in the jet's co-moving frame. It is found that the turbulence scatters photons to high energies and thus intensifies the emission in the high-energy regime. The corresponding distortion of the radiation field can be preserved if and only if the turbulence occurs in a region with incomplete photon-electron coupling. Consequently, the observed radiation spectrum can be reshaped into a Band-like spectrum.

Stellar-mass black holes ($3$ $M_\odot \lesssim M_{\rm BH} \lesssim 150$ $M_\odot$) are the natural product of the evolution of heavy stars ($M_{\rm star} \gtrsim 20$ $M_\odot$). In our Galaxy, we expect $10^8$-$10^9$ stellar-mass black holes formed from the gravitational collapse of heavy stars, but currently we know less than 100 objects. We also know $\sim 100$ stellar-mass black holes in other galaxies, most of them discovered by gravitational wave observatories in the past 10 years. The detection of black holes is indeed extremely challenging and possible only in very special cases. This article is a short review on the physics and astrophysics of stellar-mass black holes.

In this work, we present a study on the H$\alpha$ emission line flux concentration of 3098 low-redshift star-forming galaxies (SFGs) using the MaNGA data available in the Data Release 17 from the Sloan Digital Sky Survey. We define the H$\alpha$ flux concentration index ($C_{\rm H\alpha}$) as $C_{\rm H\alpha}=F_{\rm H\alpha,0.8~R_e}/F_{\rm H\alpha,1.5~R_e}$, where $F_{\rm H\alpha,0.8~R_e}$ and $F_{\rm H\alpha,1.5~R_e}$ are the cumulative H$\alpha$ flux inside $0.8$ and $1.5$ $r-$band effective radius, respectively. We find that $C_{\rm H\alpha}$ is strongly correlated with the luminosity weighted stellar age gradient. $C_{\rm H\alpha}$ is also sensitive to environmental effects, in the sense that low-mass satellite galaxies below the star formation main sequence tend to have higher $C_{\rm H\alpha}$. For central galaxies, we find that massive disk galaxies with enhanced star formation rate tend to have higher $C_{\rm H\alpha}$, while such a phenomenon is not seen in the low-mass regime. We interpret this as evidence that compaction events more frequently occur in the high-mass regime, which eventually resulting in the buildup of prominent bulges in massive SFGs. Implications of these findings on galaxy structure formation are discussed.

Gamma-ray bursts are expected to be generated by structured jets, whose profiles significantly impact their afterglow emission. Previously, we developed a numerical code jetsimpy, to model the afterglow of jets with arbitrary angular profiles. In this study, we extend the code to incorporate a stratified radial profile, enabling it to model jets with arbitrary axisymmetric two-dimensional structures. The radial profile leads to the formation of a reverse shock. We modeled the shock system using an energy conservation prescription, which differs from the pressure balance approach. This leads to remarkably different predictions for reverse shock emission. In particular, we find that the reverse shock emission in the thin shell case is significantly overestimated in analytic models. We also explore the off-axis reverse shock emission from structured jets, where the cores belong to thick shell cases and the wings belong to thin shell cases. We have confirmed the prediction that off-axis observers may see a thin-to-thick transition, but we find that the light curve morphology is hard to distinguish from pure thin or thick shell cases. A radial profile also introduces hydrodynamic energy injection. As such, our code can naturally apply to refreshed shock cases, where the modeling of kilonova afterglows is demonstrated as an example. To validate our method, we fit the optical flash of GRB 990123, showing good agreement with the data. The upgraded jetsimpy provides unprecedented flexibility in modeling the afterglow emission of jets with various profiles, including those derived from general relativistic magnetohydrodynamic simulations.

Despite the ever-increasing number of known exoplanets, no uncontested detections have been made of their satellites, known as exomoons. The quest to find exomoons is at the forefront of exoplanetary sciences. Certain space-born instruments are thought to be suitable for this purpose. We show the progress made with the CHaracterizing ExOPlanets Satellite (CHEOPS) in this field using the HD 95338 planetary system. We present a novel methodology as an important step in the quest to find exomoons. We utilize ground-based spectroscopic data in combination with Gaia observations to obtain precise stellar parameters. These are then used as input in the analysis of the planetary transits observed by CHEOPS and the Transiting Exoplanet Survey Satellite (TESS). In addition, we search for the signs of satellites primarily in the form of additional transits in the Hill sphere of the eccentric Neptune-sized planet HD 95338b in a sequential approach based on four CHEOPS visits. We also briefly explore the transit timing variations of the planet. We present refined stellar and planetary parameters, narrowing down the uncertainty on the planet-to-star radius ratio by a factor of $10$. We also pin down the ephemeris of HD 95338b. Using injection/retrieval tests, we show that a $5 \sigma$ detection of an exomoon would be possible at $R_{\rm Moon} = 0.8$~$R_\oplus$ with the methodology presented here. We exclude the transit of an exomoon in the system with $R_{\rm Moon} \approx 0.6$~$R_\oplus$ at the $1\sigma$ level. The algorithm used for finding the transit-like event can be used as a baseline for other similar targets, observed by CHEOPS or other missions.

Long-period radio transients (LPTs) represent a recently uncovered class of Galactic radio sources exhibiting minute-to-hour periodicities and highly polarised pulses of second-to-minute duration. Their phenomenology does not fit exactly in any other class, although it might resemble that of radio magnetars or of white dwarf (WD) radio emitting binary systems. Notably, two LPTs with confirmed multi-wavelength counterparts have been identified as synchronised white dwarf-M dwarf binaries (polars). Meanwhile, systems such as AR Scorpii and J1912-44 exhibit short-period pulsations in a hr-tight orbit, with polarized radio emission proposed to be generated by the interaction of the WD magnetosphere with the low-mass companion wind. Here we demonstrate that both LPTs and WD binary pulsars can be explained within a single geometric model in which radio emission is triggered when the magnetic pole of a rotating white dwarf intersects its companion's wind in the binary orbital plane. We use a 36-year timing baseline to infer the orbital period and binary geometry solely from radio data of GPM J1839-10, the longest-active LPT known. The model naturally predicts its intermittent emission and double-pulse structure. Crucially, we show that the beat period between the spin and the orbit matches the observed pulse substructure and polarisation signatures, providing strong support for the model. Applying this model to the WD binary pulsar J1912-44, it can also reproduce the system emission and geometry. Our results place GPM J1839-10, and other LPTs in general and radio emitting WD binaries, at different stages of a continuum between intermediate and synchronised polars, suggesting a unified population of magnetic WD binaries driving coherent radio emission.

In the era of photometry with space-based telescopes, such as CHEOPS (CHaracterizing ExOPlanets Satellite), JWST (James Webb Space Telescope), PLATO (PLAnetary Transits and Oscillations of stars), and ARIEL (Atmospheric Remote-sensing Infrared Exoplanet Large-survey), the road has opened for detecting subtle distortions in exoplanet transit light curves -- resulting from their non-spherical shape. We investigate the prospects of retrieval of rotational flatness (oblateness) of exoplanets at various noise levels. We present a novel method for calculating the transit light curves based on the Gauss-Legendre quadrature. We compare it in the non-rotating limit to the available analytical models. We conduct injection-and-retrieval tests to assess the precision and accuracy of the retrievable oblateness values. We find that the light curve calculation technique is about $25$\% faster than a well-known analytical counterpart, while still being precise enough. We show that a $3 \sigma$ oblateness detection is possible for a planet orbiting bright enough stars, by exploiting a precise estimate on the stellar density obtained e.g. from asteroseismology. We also show that for noise levels $\geq 256$ ppm (expressed as point-to-point scatter with a $60$~s exposure time) detection of planetary oblateness is not reliable.

To understand the chemical origin of the Solar system, the chemical evolution along the star/planet formation is a key issue. Extensive observational studies have demonstrated a chemical diversity in young low-mass protostellar sources so far. Furthermore, chemical differentiations in the vicinity of the protostars have recently been reported. This suggests that molecular distribution is sensitive to a change in the physical conditions associated with disk formation. Some kinds of molecular lines, especially Sulfur-bearing species, are therefore prospected to work as molecular markers to highlight particular structures of disk-forming regions. Conversely, detailed physical characterization is essential for elucidating the chemical evolution occurring there. Machine learnings may help us to disentangle the observed structures. Angular momentum of the gas is the key topic to understand the structure formation, which is also essential to the integration of the chemical and physical characterization.

The evolution of the supermassive Black Hole (BH) population across cosmic times remains a central unresolved issue in modern astrophysics, due to the many noticeable uncertainties in the involved physical processes. Here we tackle the problem via a semi-empirical approach with minimal assumptions and data-driven inputs. This is firmly grounded on a continuity plus Smoluchowski equation framework that allows to unitarily describe the two primary modes of BH growth: gas accretion and binary mergers. Key quantities related to the latter processes are incorporated through educated parameterizations, and then constrained in a Bayesian setup from joint observational estimates of the local BH mass function, of the large-scale BH clustering, and of the nano-Hz stochastic gravitational wave (GW) background measured from Pulsar Timimg Array (PTA) experiments. We find that the BH accretion-related parameters are strongly dependent on the local BH mass function determination: higher normalizations and flatter high-mass slopes in the latter imply lower radiative efficiencies and mean Eddington ratios with a stronger redshift evolution. Additionally, the binary BH merger rate is estimated to be a fraction $\lesssim 10^{-1}$ of the galaxy merger rate derived from galaxy pairs counts by JWST, and constrained not to exceed the latter at $\gtrsim 2\sigma$. Relatedly, we highlight hints of a possible tension between current constraints on BH demographics and the interpretation of the nano-Hz GW background as predominantly caused by binary BH mergers. Specifically, we bound the latter's contribution to $\lesssim 30-50\%$ at $\sim 3\sigma$, suggesting that either systematics in the datasets considered here have been underestimated so far, or that additional astrophysical/cosmological sources are needed to explain the residual part of the signal measured by PTA experiments.

The launch of the Imaging X-ray Polarimetry Explorer (IXPE), the first space-based polarimeter since 1978, offers a two order of magnitude improvement to the measurement of X-ray polarisation than its predecessor OSO-8, offering unprecedented precision for the measurement of polarisation degree and polarisation angle of X-ray sources. This advancement lends itself to the birth of a number of contemporary techniques to study Galactic compact objects, including X-ray polarimetry-timing, the study of how polarisation properties evolve over short timescales. However, the statistical nature of polarisation measurements poses a challenge for studies on arbitrarily short timescales, as a large number of photons are required to achieve statistically significant measurements of polarisation degree and angle for time-resolved analyses. Furthermore, if the polarisation variability is stochastic, then phase-folding techniques introduce systematic errors in the phase assignment of photons. Ingram and Maccarone presented a model independent Fourier-based technique that circumvents these issues. It can be used on arbitrarily short timescales for any kind of variability, whether aperiodic, quasi-periodic or purely periodic. Here we implement this method on real IXPE data. We address several instrumental effects and test the technique on X-ray pulsars, RX-J0440.9+4431 and Hercules X-1 . We verify that our technique recovers the polarisation variability signal that we already know to be there from typical phase-folding techniques. It will now be possible to study fast stochastic polarisation variability of X-ray sources, with applications including quasi-periodic oscillations, mass accretion rate fluctuations, and reverberation mapping.

Gamma-ray emission in the GeV-TeV range from the solar disk is likely to arise from collisions of galactic cosmic rays (GCRs) with solar atmospheric plasma. In a previous study, we demonstrated that closed turbulent magnetic arcades trap efficiently GCRs leading to a gamma-ray flux consistent with the Fermi-HAWC observations (from $\sim 0.1$ GeV to $\sim 1$ TeV). Here, we model a synthetic magnetic field with a static, laminar structure of open field lines in the chromosphere increasingly braiding near the solar surface, with a scale height of $\sim 10^{-2} R_\odot$. The height-dependent increase in magnetic field line braiding is modulated by an exponential scalar function, mimicking the bending of the photo- and chromo-spheric magnetic field revealed by polarimetric observations and reproduced by MHD simulations. Employing 3D test-particle numerical simulations, we investigate how distorted magnetic field lines affect the gamma-rays production by injecting GeV-TeV protons into both magnetically laminar and braided regions. We find that with the chosen spatial resolution this synthetic magnetic field can account for the $> 10$ GeV gamma-ray spectrum observed by Fermi-LAT/HAWC. A rebrightening between approximately $30$ and $100$ GeV (following a $\sim 30$ GeV spectral dip), suggests an enhanced confinement within the photo-/chromospheric layer by a stronger braiding.

Context. The origin and nature of spiral arms remain unclear. Star forming regions and young stars are generally strongly associated to the spiral structure, but there are few quantitative predictions from simulations about the involvement of stars of different ages. Aims. We aim to quantify the interplay between spiral arms and different populations. Methods. We use a hydrodynamical simulation of an isolated disc galaxy displaying a dynamic multi-armed spiral structure. Inspired by cosmological structure metrics, we develop a new method, the local dimension, that robustly delineates arms across populations and through space and time. Results. We find that all stars, including those as old as 11Gyr, support the arms. The spiral strength decreases with stellar age up to 2Gyr-old stars and remains nearly constant for older stars. However, the scaling between arm strength and age (or velocity dispersion) depends on the strength of the global spiral structure at each time. Almost all stars formed in arms remain within them for no more than 140-180Myr, whereas old stars leave arms about three times faster. Even if the youngest populations dominate in the production of the spiral torques at early times, all populations contribute equally at later times. Conclusions. Our results highlight the power of the local dimension for studying complex spiral structures and show that all stellar populations in the disc partake in the arms. Since in our model we see spiral arms in populations with velocity dispersions up to 90km/s, which are comparable to those of the Milky Way, we predict that old Galactic populations could also exhibit spiral structure.

We report the detection and characterization of ultrafast outflows (UFOs) in the X-ray spectra of the tidal disruption event (TDE) AT2020afhd, based on observations from NICER, Swift, and XMM. Prominent blueshifted absorption features were detected exclusively during the intermediate phase of the event, occurring between days 172 and 212 within the first 300 days post-discovery. During this period, the UFO appeared no earlier than day 74, strengthened between days 172 and 194, and disappeared after day 215. This marks the first time that the full evolutionary sequence of X-ray outflows has been observed in a TDE. Moreover, the outflows exhibited a dramatic deceleration from ~0.19c to ~0.0097c over a span of approximately 10 days. Photoionization spectral analysis reveals an inverse correlation between outflow velocity and ionization parameter, in contradiction to the predictions from radiation pressure-driven wind. Eventually, we propose that the delayed onset of the outflows may result from an increase in the wind opening angle and/or metal enrichment, particularly iron and oxygen, during the disk formation phase.

Surveys for exoplanets indicate that the occurrence rate of gas giant planets orbiting late-type stars in orbits with periods shorter than 1000 days is lower than in the case of Sun-like stars. This is in agreement with planet formation models based on the core or pebble accretion paradigm. The CARMENES exoplanet survey has been conducting radial-velocity observations of several targets that show long-period trends or modulations that are consistent with the presence of giant planets at large orbital separations. We present an analysis of five such systems that were monitored with the CARMENES spectrograph, as well as with the IRD spectrograph. In addition, we used archival data to improve the orbital parameters of the planetary systems. We improve the parameters of three previously known planets orbiting the M dwarfs GJ 317, GJ 463, and GJ 3512. We also determine the orbital parameters and minimum mass of the planet GJ 3512 c, for which only lower limits had been given previously. Furthermore, we present the discovery of two new giant planets orbiting the stars GJ 9733 and GJ 508.2, although for the second one only lower limits to the orbital properties can be determined. The new planet discoveries add to the short list of known giant planets orbiting M-dwarf stars with subsolar metallicity at long orbital periods above 2000 days. These results reveal that giant planets appear to form more frequently in wide orbits than in close-in orbits around low-mass and lower metallicity stars.

Inflation with an inflection point potential is a popular model for producing primordial black holes. The potential near the inflection point is approximately flat, with a local maximum next to a local minimum, prone to eternal inflation. We show that a sufficient condition for eternal inflation is $\lambda_1 \leq 3$, where $\lambda_1$ is the index of the `exponential tail,' the lowest eigenvalue of the Fokker--Planck equation over a bounded region. We write $\lambda_1$ in terms of the model parameters for linear and quadratic regions. Wide quadratic regions inflate eternally if the second slow-roll parameter $\eta_V \geq -6$. We test example models from the literature and show this condition is satisfied; we argue eternal inflation is difficult to avoid in inflection point PBH models. Eternally inflating regions correspond to type II perturbations and form baby universes, hidden behind black hole horizons. These baby universes are inhomogeneous on large scales and dominate the multiverse's total volume. We argue that, if volume weighting is used, eternal inflation makes inflection point primordial black hole models incompatible with large-scale structure observations.

We investigate the observational signatures of outflow rotation in protostellar systems using magnetohydrodynamics simulations of protostellar evolution with radiative transfer and synthetic observation. The velocity gradient perpendicular to the outflow axis indicates outflow rotation. The rotation signature is clearly seen in the moment 1 map and a position-velocity (PV) diagram across an outflow lobe made from our model with an inclination angle of i>~85{degree sign}, as in observational studies of protostellar outflows. Velocity projection with lower inclinations distorts the moment 1 map because the outflow vertical (propagation) velocity contributes more to the line-of-sight velocity, leading to an incorrect outflow axis direction. The PV diagram adopting the incorrect outflow axis shows no clear velocity gradient. These effects may prevent us from identifying outflow rotation. Our analysis implies that rotational signatures can be obscured in ~2/3 to ~4/5 of the total outflow population (i<70{degree sign}-80{degree sign}), regardless of the evolutionary stage. Complicated structures in observed outflows make it difficult to determine the outflow, which may result in the apparent non-detection of outflow rotation.

Recently discovered supermassive black holes with masses of $\sim10^8\,M_\odot$ at redshifts $z\sim9$-$11$ in active galactic nuclei (AGN) pose severe challenges to our understanding of supermassive black hole formation. One proposed channel are rapidly accreting supermassive PopIII stars (SMSs) that form in large primordial gas halos and grow up to $<10^6\,M_\odot$. They eventually collapse due to the general relativistic instability and could lead to supernova-like explosions. This releases massive and energetic ejecta that then interact with the halo medium via an optically thick shock. We develop a semi-analytic model to compute the shock properties, bolometric luminosity, emission spectrum and photometry over time. The initial data is informed by stellar evolution and general relativistic SMS collapse simulations. We find that SMS explosion light curves reach a brightness $\sim10^{45\mathrm{-}47}\,\mathrm{erg/s}$ and last $10$-$200$ years in the source frame - up to $250$-$3000$ years with cosmic time dilation. This makes them quasi-persistent sources which vary indistinguishably to little red dots and AGN within $0.5$-$9\,(1+z)$ yrs. Bright SMS explosions are observable in long-wavelength JWST filters up to $z\leq20$ ($24$-$26$ mag) and pulsating SMSs up to $z\leq15$. EUCLID and the Roman space telescope (RST) can detect SMS explosions at $z<11$-$12$. Their deep fields could constrain the SMS rate down to $10^{-11}$Mpc$^{-3}$yr$^{-1}$, which is much deeper than JWST bounds. Based on cosmological simulations and observed star formation rates, we expect to image up to several hundred SMS explosions with EUCLID and dozens with RST deep fields.

As part of Saudi Arabia Vision 2030 and under the guidance of the Royal Commission for AlUla (RCU), efforts are underway to establish AlUla Manara as the Kingdom first major astronomical observatory. This study presents a preliminary assessment of the site based on ECMWF ERA5 reanalysis data to evaluate its suitability for hosting a 4m-class optical-IR telescope. AlUla Manara is located on a remote plateau 74 km north of the historical town of AlUla and was recently designated as an International Dark Sky Park. The analysis focuses on key astro-meteorological parameters such as seeing, temperature regimes, wind patterns, cloud cover and precipitable water vapor (PWV). Results show a median nighttime seeing of 1.5 arcsec, a median cold season PWV of 3.2 mm, and over 79% of nighttime hours with clear sky conditions. Wind regimes are generally mild, posing no constraints on infrastructure. The analysis includes three further sites in the Kingdom, namely Volcanic Top, Ward Mountain, and Dubba Mountain. These sites exhibit better turbulence conditions, but are located outside RCU jurisdiction. Nevertheless, AlUla Manara remains a competitive candidate thanks to its alignment with broader regional development goals. To validate these preliminary results, a dedicated Astronomical Site Monitor has been deployed on site to support the design and operational planning of the observatory.

In recent years, with an increasing number of type Ia supernovae (SNe Ia) discovered soon after their explosions, a non-negligible fraction of SNe Ia with early-excess emissions (EExSNe Ia) have been confirmed. In this letter, we present a total of \textbf{67} early-phase normal SNe Ia from published papers and ongoing transient survey projects to systematically investigate their photometric behaviors from very early time. We found that EExSNe Ia in our sample have longer rise and brighter peak luminosities compared to those of non-EExSNe Ia. Moreover, EExSNe Ia commonly have ``red-bump" features in the early $B-V$ color while non-EExSNe Ia show blueward evolution from the very beginning. Here, we propose that the thin-helium double-detonation scenario can phenomenologically explain the photometric diversities of normal SNe Ia considering different white dwarf-He-shell mass combinations and the viewing-angle effect, implying a unified explosion mechanism of normal-type SNe Ia. To further testify the possible common origin of normal SNe Ia, systematical studies of multiband photometric and spectral properties of early-phase SNe Ia through the new generation wide-field time-domain survey facilities and global real-time follow-up networks are highly demanded.

The ASTRI Mini-Array is an international collaboration led by the Italian National Institute for Astrophysics. The project aims to construct and operate an array of nine Imaging Atmospheric Cherenkov Telescopes to study gamma-ray sources at very high energy (TeV) and perform stellar intensity interferometry observations. We describe the updated Online Observation Quality System (OOQS) software architecture. The OOQS is one of the subsystems of the Supervisory Control and Data Acquisition (SCADA) system. It aims to execute real-time data quality checks on the data acquired by the Cherenkov cameras and intensity interferometry instruments and provide feedback to both SCADA and the Operator about abnormal conditions detected. The data quality results are stored in the Quality Archive for further investigation and sent to the Operator Human Machine Interface (HMI) through Kafka.

Downward Terrestrial Gamma-ray Flashes (TGFs) are sub-millisecond bursts of MeV gamma rays produced in thunderclouds. According to the Relativistic Runaway Electron Avalanche model, gamma rays are produced, via bremsstrahlung, from electron cascades activated by a relativistic "seed" electron. It is not clear what mechanism is responsible for the acceleration of electrons to relativistic energies in electric discharges. To better understand the acceleration sites and the TGF production mechanisms, it is critically important to identify the TGF source position and geometry in the atmosphere and to study the gamma emission characteristics. The Surface Detector of the Pierre Auger Observatory, with its 1600 water-Cherenkov detectors very sensitive to high-energy photons and with a very fine time-sampling, is a valuable instrument to study downward TGFs. The possibility to analyze the radiation emission in detail led to the observation of the first TGFs with an asymmetric azimuthal structure, suggesting a complex source different from the initially hypothesized downward beam. We report on these observations and the new perspectives which may open with the incorporation of new instruments at the Auger site to study lightning development alongside gamma emission, and the increasingly detailed data provided by satellites and global lightning networks.

The mass distribution of binary black holes inferred from gravitational wave measurements is expected to shed light on their formation scenarios. An emerging structure in the mass distribution indicates the presence of multiple peaks around chirp masses of $8M_\odot$, $14M_\odot$, and $27M_\odot$. In particular, there is a lack of observations between chirp masses of 10 and 12 $M_\odot$. In this letter, we report that observations significantly favour the model supporting suppression of the rate in a narrow chirp mass range compared to the model that doesn't include suppression at a confidence greater than 99.5\%. Using another test, which measures the deviation between the inferred chirp mass distributions from the two models, we conservatively estimate a 95\% confidence in the presence of a feature. A lack of confidence has been reported in the presence of a gap around a comparable location in the component mass distribution. The differing conclusions are due to a unique correlation between the primary~(heavier of the two masses) and the secondary~(lighter of the two masses) masses of binary black holes. This correlation results in increased clustering of measured chirp masses around specific values.

Environmental effects within cosmological overdensities, such as galaxy groups and clusters, have been shown to impact galaxies and their cold gas reservoirs and thereby provide constraints on galaxy evolution models. Galaxy groups foster frequent galaxy-galaxy interactions, making them rich environments in which to study galaxy transformation. In this work, we study a serendipitously discovered large overdensity of HI galaxies at z~0.04. The galaxies appear to lie in a filamentary-like structure of megaparsec scale. MeerKAT's angular resolution and field of view allow us to spatially resolve the HI galaxies while simultaneously probing large-scale structure. The HI and sub-arcsec Dark Energy Survey (DES) imaging reveal a large number of interacting galaxies in this collective group. MeerKAT data enables us to derive HI masses and investigate interacting galaxies. We use DES and Wide-field Infrared Survey Explorer (WISE) data to quantify the star formation rates, stellar masses, and stellar morphologies of member galaxies and compare these with field scaling relations. To place this discovery and the environmental effects in context, we use the SIMBA simulation to investigate the prevalence of qualitatively similar HI overdensities and their large-scale morphological properties. This enables a prediction of how frequently such structures might be serendipitously discovered with MeerKAT and SKA-Mid HI observations in comparable observation time. The combination of spatially resolved HI data and optical imaging reveals a group rich in interactions, suggesting environmental processes are already shaping galaxy properties within the structure. More of these serendipitous discoveries are expected and, alongside ongoing targeted programmes, these will provide a rich sample to study galaxy transformation and enable a MeerKAT HI perspective on large-scale structure, including filaments.

We investigate the possibility that the recently reported GW231123 event, with component masses $M_1=137^{+22}_{-17}\,M_\odot$, $M_2=103^{+20}_{-52}\,M_\odot$ and a local merger rate $R_{\mathrm{local}}=0.08^{+0.19}_{-0.07}\,\mathrm{Gpc^{-3}\,yr^{-1}}$, originates from primordial black holes (PBHs) formed during an early matter-dominated era. We compute the PBH mass function, abundance, spin distribution and the merger rate density and find a set of choices for the parameters to reproduce the key properties of GW231123. We also show that the resulting PBH abundance, $f_{\mathrm{pbh}}= 1.64^{+5.00}_{-1.59}\times10^{-1}$, is at the border of the exclusion regions, but still remains self-consistent. Finally, we estimate the scalar-induced gravitational waves (SIGWs) that are inevitably generated during PBH formation. PBHs that interpret GW231123 are accompanied by negligible SIGWs in the nano-hertz band, indicating no conflict with current pulsar timing arrays data.

Simulated maps of the microwave background (CMB) radiation are generally created using one of two methods: all-sky simulations use the spherical harmonic transform, while maps covering small areas approximate the sky as flat, allowing the use of fast Fourier transforms (FFTs). Current and near-future experiments, particularly ones like CMB S4, will cover areas too large for the flat-sky approximation but significantly less than the full sky. In this regime, it can be more efficient to simulate maps in a 3-D box using FFTs, and then sample onto the observed part of the celestial sphere. We present a method for performing such simulations and show that it can be more efficient than full-sky simulations. We develop the method for scalar maps, but we expect it to be applicable to higher-spin (e.g., polarization) simulations as well.

To explore the formation and properties of Saturn's G ring, we study the dynamics of micron-sized dust particles originating from the arc of debris near the inner edge of the ring. The dynamical evolution of particles due to various perturbation forces and the plasma sputtering that erodes the particles is simulated by a well-tested numerical code. Based on the simulation results, the normal $I/F$ of the G ring observed by the Cassini spacecraft can be explained by dust particles originating from the arc. Other properties of the G ring are also estimated, including the steady-state size distribution and the number density of ring particles, the geometric optical depth, the apparent edge-on thickness, the age and the remaining lifetime of the G ring. We find that the particle size distribution of the G ring follows a power law with an exponent of 2.8, and dust particles in the size range of $[5, 10]\,{\mu}$m are dominant within the ring. The average number density of particles of the G ring in the radial direction is about $10^{-3}$-$10^{-2}\,\mathrm{m}^{-3}$. The peak value of the edge-on geometric optical depth of the G ring is about $3.9\times10^{-2}$. The maximum apparent edge-on thickness of the G ring with the geometric optical depth larger than $1\times10^{-8}$ is approximately $9,000\,\mathrm{km}$. The age of the G ring is estimated to be $10^{6}$-$10^{7}\,\mathrm{years}$, and the remaining lifetime of the ring is on the order of $10^{4}\,\mathrm{years}$.

Betelgeuse -- the closest M-supergiant to the Sun -- has recently been predicted to host a lower-mass stellar companion that orbits the primary with a period of $\sim 6$ years. The putative stellar companion is thought to cause the long photometric modulation observed in Betelgeuse, which cannot be explained by stellar pulsations. Additionally, radial velocity and astrometric data also point to a stellar companion. Here we present diffraction-limited optical speckle imaging observations obtained on the 8.1-meter Gemini North telescope in 2020 and 2024. The 2020 observations were taken during the Great Dimming event, and at a time when the stellar companion was predicted to be unobservable because it was directly in-line with Betelgeuse itself. The 2024 observations were taken three days after the predicted time of greatest elongation for the companion. A comparison of the 2020 and 2024 data reveal no companion in 2020 (as expected) and the probable detection of a companion in 2024. The presumed stellar companion has an angular separation and position angle of 52 mas and $115^\circ$ east of north, respectively, which is in excellent agreement with predictions from dynamical considerations. The detected companion is roughly six magnitudes fainter than Betelgeuse at 466 nm.\bf{While this is only a 1.5$\sigma$ detection, the results are in reasonable agreement with the predictions: the appearance of the companion at quadrature; the angular separation from Betelgeuse; the position angle with respect to Betelgeuse; the magnitude difference; and the estimated mass of the companion.

The Iron Stars XX Oph and AS 325, both binaries consisting of a Be + late K,M II star present remarkable optical spectra. Low velocity, dense winds from the late type star collide with high velocity, optically thin material expelled from the hot Be star producing a plethora of emission lines and complex P Cygni absorption profiles. The members of the American Association of Variable Star Observers have faithfully observed XX Oph for 85 years and AS 325 for 32 years. This long term photometric monitoring has revealed AS 325 to be an eclipsing system while XX Oph shows a complex light curve behavior. We present archival and recent photometry, new high-resolution optical imaging, and new high-resolution optical spectroscopy of the two stars. The orbital period of AS 325 is refined as 512.943 days and the stellar components are Be+K2.5 II. XX Oph is shown to be non-eclipsing and consists of Be+M6II stars.

We investigate the feasibility of a spacecraft mission to conduct a flyby of 3I/ATLAS, the third macroscopic interstellar object discovered on July 1 2025, as it traverses the Solar System. There are both ready-to-launch spacecraft currently in storage on Earth, such as Janus, and spacecraft nearing the end of their missions at Mars. We calculate minimum $\Delta V$ single-impulse direct transfer trajectories to 3I/ATLAS both from Earth and from Mars. We consider launch dates spanning January 2025 through March 2026 to explore obtainable and hypothetical mission scenarios. Post-discovery Earth departures require a challenging $\Delta V\gtrsim24$ km s$^{-1}$ to fly by 3I/ATLAS. By contrast, Mars departures from July 2025 - September 2025 require $\Delta V\sim5$ km s$^{-1}$ to achieve an early October flyby -- which is more feasible with existing propulsion capabilities. We discuss how existing spacecraft could be used to observe 3I/ATLAS and how spacecraft at other locations in the Solar System could be repurposed to visit future interstellar objects on short notice.

We revisit the theoretical modeling and simulation of a Gaussian stochastic gravitational wave background (SGWB) signal in a pulsar timing array (PTA). We show that the correlation between Fourier components of pulsar timing residuals can be expressed using transfer functions; that are indicative of characteristic temporal correlations in a SGWB signal observed in a finite time window. These transfer functions, when convolved with the SGWB power spectrum and spatial correlation (Hellings \& Downs curve), describe the variances and correlations of the pulsar timing residuals' Fourier coefficients. The convolutions are the exact frequency- and Fourier-domain representations of the time-domain covariance function. We derive explicit forms for the transfer functions for unpolarized and circularly polarized SGWB signals. We validate our results by comparing Gaussian theoretical expectation values with standard simulations based on point sources and our own covariance-matrix-based approach. The unified frequency- and Fourier-domain formalism provides a robust foundation for future PTA precision analyses and highlights the importance of temporal correlations in interpreting GW signals.

The case of Ross 176 is a late K-type star that hosts a promising water-world candidate planet. The star has a radius of $R_*$=0.569$\pm$0.020$R_{\odot}$ and a mass of $M_{\star}$ = 0.577 $\pm$ 0.024 $M_{\odot}$. We constrained the planetary mass using spectroscopic data from CARMENES, an instrument that has already played a major role in confirming the planetary nature of the transit signal detected by TESS. We used Gaussian Processes (GP) to improve the analysis because the host star has a relatively strong activity that affects the radial velocity dataset. In addition, we applied a GP to the TESS light curves to reduce the correlated noise in the detrended dataset. The stellar activity indicators show a strong signal that is related to the stellar rotation period of $\sim$ 32 days. This stellar activity signal was also confirmed on the TESS light curves. Ross 176b is an inner hot transiting planet with a low-eccentricity orbit of $e = 0.25 \pm 0.04$, an orbital period of $P \sim 5$ days, and an equilibrium temperature of $T_{eq}\sim 682K$. With a radius of $R_p = 1.84\pm0.08R_{\oplus}$ (4% precision), a mass of $M_p = 4.57^{+0.89}_{-0.93} M_{\oplus}$ (20% precision), and a mean density of $\rho_p = 4.03^{+0.49}_{-0.81} g cm^{-3}$, the composition of Ross 176b might be consistent with a water-world scenario. Moreover, Ross 176b is a promising target for atmospheric characterization, which might lead to more information on the existence, formation and composition of water worlds. This detection increases the sample of planets orbiting K-type stars. This sample is valuable for investigating the valley of planets with small radii around this type of star. This study also shows that the dual detection of space- and ground-based telescopes is efficient for confirm new planets.

A fraction of the Ultra Luminous X-ray (ULX) sources are known to be accreting neutron stars as they show coherent X-ray pulsations with pulse periods ranging from ~1-30 seconds. While initially thought to host intermediate-mass black holes, ULXs have since been recognized as a diverse class of objects, including ULX pulsars. These pulsars require models specifically tailored to account for their unique accretion physics, distinct from those used for Galactic black hole binaries. The X-ray spectra of all Galactic accreting X-ray pulsars (including sources in the Magellanic Clouds) are dominated by a high energy cut-off power-law and some of the sources show a soft excess, some emission lines, cyclotron absorption features, etc. In this work, we undertake a comprehensive analysis of the broadband X-ray spectra of five ULX pulsars using simultaneous XMM-Newton and NuSTAR observations and show that their X-ray spectra can be effectively described by spectral models, similar to those used for the local accretion-powered X-ray pulsars. A soft excess is detected in all the sources which is also consistent with the local X-ray pulsars that have low absorption column density. We have marginal detection or low upper limit on the presence of the iron K-alpha emission line from these sources, which is a key difference of the ULX pulsars with the local accreting X-ray pulsars. We discuss the implication of this on the nature of the binary companion and the accretion mechanism in the ULX pulsars.

The origin of fast radio bursts (FRBs), the brightest cosmic radio explosions, is still unknown. Bearing critical clues to FRBs' origin, the long-term evolution of FRBs has yet to be confirmed, since the field is still young and most FRBs were seen only once. Here we report clear evidence of decadal evolution of FRB~20121102A, the first precisely localized repeater. In conjunction with archival data, our FAST and GBT monitoring campaign since 2020 reveals a significant 7% decline of local dispersion measure (DM). The rotation measure (RM) of 30,755$\pm$16 $\mathrm{rad\,m^{-2}}$ detected in the last epoch represents a 70% decrease compared to that from December 2016. The $\sigma_{RM}$ parameter, which describes the complexity of the magneto-ionic environment surrounding the source, was shown to have decreased by 13%. These general trends reveal an evolving FRB environment, which could originate from an early-phase supernova associated with an enhanced pair wind from the FRB central engine.

The detection of the Stochastic Gravitational Wave Background (SGWB) is one of the most challenging tasks for both current and next-generation detectors. Successfully distinguishing the SGWB from instrumental noise and environmental effects requires accurate and flexible analysis tools capable of detecting the signal and determining its origin. In this paper, we introduce a unified framework and a user-friendly tool for SGWB characterization: \texttt{GWBird} (Gravitational Wave Background Inventory of Response functions for Detectors). This code enables the computation of overlap reduction functions (ORFs), power-law integrated sensitivity curves (PLS), angular response functions, and angular PLS (APLS). It supports the full range of gravitational wave polarization modes (tensor, scalar, and vector), allowing for the characterization of both isotropic and anisotropic SGWB components for all the polarizations. Additionally, the code includes functions for circular polarization characterization, which is particularly relevant for probing parity-violating signals. The framework integrates analyses for ground-based, space-based, and Pulsar Timing Array (PTA) detectors, offering a versatile framework for SGWB analysis. The \texttt{GWBird} code is publicly available at:~\github{this https URL}

In our previous work, we identified $\sim100,000$ metal-poor stars ([Fe/H] $<$ -1.0) from the LAMOST Survey. This work estimates their chemical abundances and explores the origin and evolution of the Galactic metal-poor disc. Our chemo-dynamical analysis reveals four main populations within the metal-poor disc: (1) a primordial disc older than 12 Gyr with [Fe/H] $>$ -1.5; (2) debris stars from the progenitor galaxy of Gaia-Sausage-Enceladus (GSE), but now residing in the Galactic disc; (3) the metal-poor tail of the metal-rich, high-$\alpha$ disc formed 10-12 Gyr ago, with metallicity lower limit extending to -2.0; (4) the metal-poor tail of the metal-rich, low-$\alpha$ disc younger than 8 Gyr, reaching a lower metallicity limit of -1.8. These results reveal the presence of a primordial disc and show that both high-$\alpha$ and low-$\alpha$ discs reach lower metallicities than previously thought. Analysis of merger debris reveals that Wukong, with extremely low metallicity, likely originate from merger events distinct from GSE. Additionally, three new substructures are identified: ShangGu-1, characterized by unusual [Fe/H]-eccentricity correlations; ShangGu-2, possibly heated disc stars; and ShangGu-3, which can be divided into four subgroups based on differing orbital directions, with two aligning with the previously known Nyx and Nyx-2.

Recent developments in the study of pulsar radio emission revealed that the microphysics of quantum electrodynamic (QED) pair cascades at pulsar polar caps may be responsible for generating the observed coherent radio waves. However, modeling the pair cascades in the polar cap region poses significant challenges, particularly under conditions of high plasma multiplicity. Traditional Particle-in-Cell (PIC) methods often face rapidly increasing computational costs as the multiplicity grows exponentially. To address this issue, we present a new simulation code using the Vlasov method, which efficiently simulates the evolution of charged particle distribution functions in phase space without a proportional increase in computational expense at high multiplicities. We apply this code to study $e^\pm$ pair cascades in 1D, incorporating key physical processes such as curvature radiation, radiative cooling, and magnetic pair production. We study both the Ruderman-Sutherland (RS) and the Space-charge-limited Flow (SCLF) regimes, and find quasiperiodic gap formation and pair production bursts in both cases. These features produce strong electric field oscillations, potentially enabling coherent low-frequency radio emission. We construct a unified analytic model that describes the key features of the polar cap cascade, which can be used to estimate the return current heating rate that can be used to inform X-ray hotspot models. Spectral analysis shows that a significant amount of energy is carried in superluminal modes -- collective excitations that could connect to observed radio features. Our results align with previous PIC studies while offering enhanced fidelity in both dense and rarefied regions.

Astrotourism has emerged as a powerful cross sectoral tool to promote science education, sustainable economic development, and cultural exchange. Recognising its potential, the International Astronomical Union's Office of Astronomy for Development (IAU OAD) has developed a suite of openly accessible resources to support individuals and institutions interested in implementing astrotourism initiatives globally. These resources also encourage individuals and existing businesses to broaden their offerings to include activities that use the night sky as a backdrop, such as food experiences, wellness practices, and cultural exploration. This paper offers a comprehensive summary of these resources, available on the OAD's Astrotourism Portal, and situates them within the broader context of astronomy for development work. The paper is targeted at educators, policymakers, tourism operators, grassroots organisers, and entrepreneurs, providing guidance on how they can foster inclusive, locally grounded, and sustainable astrotourism efforts, particularly in underresourced or emerging contexts.

We present new, deep narrow band imagery and discuss the nature of SDSO1, the large [O III]-emitting nebula centered 1.5 degrees SE of M31. We find strong evidence to support the hypothesis that SDSO1 is unrelated to M31 and is instead a bow shock driven by a faded, giant (D = 20 pc), 400 kyr-old, ghost planetary nebula (GPN) expelled by the symbiotic WD binary star EG Andromedae. Because EG And lags the rotation of the Milky Way, its hypersonic velocity of 107 km/s drives a shock into the local interstellar medium. SDSO1 also sports a 45-pc long turbulent tail that crosses M31. We establish the shock-powered GPN phase as a new phase of PN evolution, and identify seven more candidate GPNe by their large size and shock-tail morphology. This includes several giant halos of younger planetary nebulae, possibly expelled by now degenerate binary companions. The interaction of an old, fast-moving GPN with the ISM generates a shock that remains visible long after the photoionized PN shell has faded below the limit of detectability. The SDSO1 GPN has reached a phase in its evolution where its outward expansion has been slowed considerably by the work it has done on the interstellar medium. Once the GPN outflow expends the rest of its kinetic energy, it will be stripped from EG And and mix into the surrounding ISM.

Aerosols is an old topic in the young field of exoplanet atmospheres. Understanding what they are, how they form, and where they go has long provided a fertile playground for theorists. For observers, however, aerosols have been a multi-decade migraine, as their chronic presence hides atmospheric features. For hot Jupiters, the large day-night temperature contrast drives inhomogeneous thermal structures and aerosol distribution, leading to different limb properties probed by transit spectra. We present JWST NIRISS/SOSS spectra of morning and evening limbs for nine gas giants with equilibrium temperatures of ~800-1700 K. By measuring feature size of the 1.4 $\mu$m water band for both limbs, we found three planets (WASP-39 b, WASP-94 Ab, and WASP-17 b) show prominent ($>$5$\sigma$) limb-limb atmospheric opacity difference with muted morning and clear evening limbs. The heavily muted water features on morning limbs indicate high-altitude (0.1 to 0.01 mbar) aerosols. To simultaneously have clear evening limbs requires processes with timescales ($\sim$day) comparable to advection to remove these lofted grains, and we found that both downwelling flow and dayside cloud evaporation could be plausible mechanisms. We hypothesize an empirical boundary--termed the "asymmetry horizon"--in temperature-gravity space that marks the transition where inhomogeneous aerosol coverage begins to emerge. Heterogeneous aerosol coverage is common among hot Jupiters. If unrecognized, limb averaging suppresses spectral features, mimicking high-mean-molecular-weight atmospheres, inflating inferred metallicity by up to 2 dex, and underestimating limb temperatures by as much as half. Finally, we introduce the Limb Spectroscopy Metric (LSM) to predict limb spectral feature size based on planet parameters.

The intracluster medium (ICM), composed of hot plasma, dominates the baryonic content of galaxy clusters and is primarily observable in X-rays. Its thermodynamic properties, pressure, temperature, entropy, and electron density, offer crucial insight into the physical processes shaping clusters, from accretion and mergers to radiative cooling and feedback. We investigate the thermodynamic properties of galaxy clusters in the Simulating the LOcal Web (SLOW) constrained simulations, which reproduce the observed large-scale structure of the local Universe. We assess how well these simulations reproduce observed ICM profiles and explore the connection between cluster formation history and core classification. Three-dimensional thermodynamic profiles are extracted and compared to deprojected X-ray and Sunyaev - Zel'dovich (SZ) data for local clusters classified as solid cool-core (SCC), weakly cool-core (WCC), and non-cool-core (NCC) systems. We also examine the mass assembly history of the simulated counterparts to link their formation to present-day ICM properties. The simulations reproduce global thermodynamic profiles for clusters such as Perseus, Coma, A85, A119, A1644, A2029, A3158, and A3266. Moreover, they show that CC clusters typically assemble their mass earlier, while NCC systems grow through more extended, late-time merger-driven histories. WCC clusters show intermediate behavior, suggesting an evolutionary transition. Our results demonstrate that constrained simulations provide a powerful tool for linking cluster formation history to present-day ICM properties and point to possible refinements in subgrid physics as well as in resolution that could improve the agreement in cluster core regions.

In this work, we study some characteristics and gravitational signatures of the Schwarzschild black hole immersed in a Hernquist dark matter halo (SBH-HDM). We determine the black hole's remnant radius and mass, which provide useful residual information at the end of its evaporation. We then explore the luminosity of the accretion disk from the SBH-HDM model. In this way, we determine the key orbital parameters of the test particles within the accretion disk, such as angular velocity, angular momentum, energy, and the radius of the innermost stable circular orbit, based on the dark matter model parameters. We also numerically estimate the accretion disk's efficiency in converting matter into radiation. We also demonstrate that dark matter, which significantly alters the geometry surrounding a Schwarzschild black hole, influences the accretion disk's radiative flux, temperature, differential luminosity, and spectral luminosity. The stability of a black hole spacetime is determined in the eikonal regime. The Lyapunov exponent is also analyzed to quantify the stability of the particle regime and to demonstrate the infall into or escape from the black hole to infinity, as well as the quasi-normal modes. Finally, some properties of black holes are studied from a topological perspective.

We present predictions for the zero-temperature equation of state at finite isospin density using the Nambu--Jona-Lasinio (NJL) model within the Medium Separation Scheme (MSS) -- a scheme that explicitly disentangles medium effects from the ultraviolet divergent vacuum terms. Recent lattice QCD results reveal a nonmonotonic speed of sound ($c_s^2$) as a function of isospin chemical potential ($\mu_I$), exhibiting explicit violation of the conformal bound $c_s^2 = 1/3$. These findings have attracted significant theoretical interest, as established models -- including the NJL model -- failed to anticipate this behavior prior to lattice simulations. Conventional NJL implementations yield unphysical artifacts, often attributed to regularization scale sensitivity stemming from nonrenormalizability. However, in this work, we demonstrate that the standard NJL framework combined with MSS quantitatively reproduces state-of-the-art lattice data for isospin QCD.

Astrophysical compact objects are studied in the context of quadratic non-metricity gravity. The solutions to the gravitational field equations, which include fluid components, are analyzed to investigate the density and pressure properties of radio pulsars. It is explicitly demonstrated that the theoretically stable models are consistent with astronomical data, due to the geometric features of the quadratic component. Furthermore, it is shown that, in contrast to the compactness limits of black holes in general relativity, the core density can significantly exceed the density at which nuclear saturation occurs, and the surface density can also surpass the value of nuclear saturation. Additionally, it is found that the radial sound speed remains below the conformal upper bound for sound velocity established by perturbative quantum chromodynamics.

The intense magnetic fields inferred from magnetars suggest they may be strong gravitational-wave emitters. Although emissions due to hydromagnetic deformations are more promising from a detection standpoint, exterior fields also contribute a strain. However, numerical evidence suggests that the free energy of stable magnetospheric solutions cannot exceed a few tens of percent relative to the potential state, implying that the gravitational-wave luminosity cannot differ significantly between models. This prompts `universality', in the sense that the strain provides a direct probe of the near-surface field without being muddied by magnetospheric currents. Using a suite of three-dimensional, force-free, general-relativistic solutions for dipole and dipole-plus-quadrupole fields, we find that space-based interferometers may enable marginal detections out to $\lesssim$ kpc distances for slowly-rotating magnetars with fields of $\gtrsim 10^{15}$ G.

We study the Cowling approximation by analytical means as applied to a system of linear differential equations arising from models of non-radial stellar pulsation. We consider various asymptotic cases, including those of high harmonic degree and high oscillation frequency. Our methods involve a reformulation of the system in terms of an integro-differential equation for which certain Hilbert-space methods apply. By way of a more complete asymptotic study, we extend our results to certain fundamental solution sets, characterized according to certain multi-point boundary-value problems: Such asymptotics further enable us to produce sharp estimates as confirmation of our general results.

Wormholes are non-trivial topological structures that arise as exact solutions to Einstein's field equations, theoretically connecting distinct regions of spacetime via a throat-like geometry. While static traversable wormholes necessarily require exotic matter that violates the classical energy conditions, subsequent studies have sought to minimize such violations by introducing time-dependent geometries embedded within cosmological backgrounds. This review provides a comprehensive survey of evolving wormhole solutions, emphasizing their formulation within both general relativity and alternative theories of gravity. We explore key developments in the construction of non-static wormhole spacetimes, including those conformally related to static solutions, as well as dynamically evolving geometries influenced by scalar fields. Particular attention is given to the wormholes embedded into Friedmann-Lemaître-Robertson-Walker (FLRW) universes and de Sitter backgrounds, where the interplay between the cosmic expansion and wormhole dynamics is analyzed. We also examine the role of modified gravity theories, especially in hybrid metric-Palatini gravity, which enable the realization of traversable wormholes supported by effective stress-energy tensors that do not violate the null or weak energy conditions. By systematically analyzing a wide range of time-dependent wormhole solutions, this review identifies the specific geometric and physical conditions under which wormholes can evolve consistently with the null and weak energy conditions. These findings clarify how such configurations can be naturally integrated into cosmological models governed by general relativity or modified gravity, thereby contributing to a deeper theoretical understanding of localized spacetime structures in an expanding universe.

In this work, we calculate the 2025 impact factors for 16 journals in Particle Physics (and related areas such as Nuclear Physics) using citations collated by NASA/ADS (Astrophysics Data System). We then compare them to the official impact factors calculated by Clarivate. We also compare these impact factors to the median-based impact factors introduced in a previous work. We do not find any systematic bias between the ADS and official impact factors. We find the maximum relative difference between official and median-based impact factors for PTEP and Annual Review of Nuclear and Particle Science. For PTEP, this difference is due to one outlier, viz Particle Data Group with over 1400 citations in the last two years. Similarly for Annual Review of Nuclear and Particle Science, the difference is due to a review paper on Hubble tension and early dark energy with over 100 citations in the past two years. Both these papers have significantly elevated the official impact factors above its median value.

The persistent discrepancy between the experimental measurement and the Standard Model (SM) prediction of the muon's anomalous magnetic moment $(g-2)_\mu$ remains one of the most intriguing hints of physics beyond the SM. A well-motivated explanation involves a light $Z'$ gauge boson associated with a broken $U(1)_{L_\mu - L_\tau}$ symmetry. Such a boson not only resolves the $(g-2)_\mu$ anomaly, but also induces resonant interactions between high-energy cosmic neutrinos and the cosmic neutrino background (C$\nu$B), potentially shaping the observable neutrino flux at Earth. In this work, we explore the implications of such interactions for the cosmic propagation of high-energy neutrinos. We compute the optical depth for neutrino attenuation via $Z'$-mediated scattering, accounting for neutrino masses, hierarchies, and thermal distributions. We delineate the regions in $(m_{Z'}, m_\nu)$ space where the optical depth exceeds unity, defining a "neutrino cosmic horizon" beyond which high-energy neutrinos are significantly attenuated. We confront these results with the parameter space required to simultaneously explain the muon $g-2$ anomaly and ease the Hubble tension via an additional contribution to the effective number of relativistic degrees of freedom, $\Delta N_{\mathrm{eff}} \simeq 0.2-0.5$. Our analysis reveals a consistent region in parameter space where all three phenomena - $(g-2)_\mu$, $N_{\mathrm{eff}}$, and high-energy neutrino attenuation-can be explained by the same light mediator. These findings motivate future searches for spectral features in IceCube and its next-generation successors as indirect probes of new physics in the neutrino sector.

This paper presents an adaptive symplectic integrator, SQQ-PTQ, developed on the basis of the fixed-step symplectic integrator SQQ. To mitigate the Runge phenomenon, SQQ-PTQ employs Chebyshev interpolation for approximating the action, enhancing both the precision and stability of the interpolation. In addition, to reduce the computational cost of evaluating interpolation functions, SQQ-PTQ introduces a projection method that improves the efficiency of these computations. A key feature of SQQ-PTQ is its use of the time transformation to implement an adaptive time step. To address the challenge of computing complicated Jacobian matrices attributed to the time transformation, SQQ-PTQ adopts a quasi-Newton method based on Broyden's method. This strategy accelerates the solution of nonlinear equations, thereby improving the overall computational performance. The effectiveness and robustness of SQQ-PTQ are demonstrated via three numerical experiments. In particular, SQQ-PTQ demonstrates adaptability in handling close-encounter problems. Moreover, during long-term integrations, SQQ-PTQ maintains the energy conservation, further confirming its advantages as a symplectic algorithm.

The electromagnetic (EM) helicity flux density and the magnetic helicity are related by topological Chern-Simons terms. We show that the helicity flux density is distinguished from magnetic helicity by analysing Hopf solitons. We find the helicity flux density for a point charge moving with an acceleration, extending the Liénard-Wiechert angular distribution of radiant power. We also derive the multipole expansion of the helicity flux density, generalizing the Larmor's formula for the radiant power. These formulae have been applied to discuss the helicity flux density in several toy models such as circular and helical motion as well as soft bremsstrahlung. We also comment on the potential applications of the EM helicity flux density to pulsar systems.

We propose a new, general form for a pseudo-Newtonian gravitational potential (PNP), expressed as a series of Paczyński-Wiita-like functions with the addition of increasing negative powers of $r$ with arbitrary coefficients. We present a procedure for determining these coefficients to construct a custom PNP that replicates key features of Schwarzschild geodesics for a test particle near a black hole. As an example, we construct potentials set to reproduce (I) the presence of an innermost stable circular orbit at the $r=6$ (geometric units), with the correct infall velocity for small deviations (on the geodesic universal infall), (II) the periapsis advance at large distances, and (III) the presence of a marginally bound circular orbit with specific angular momentum $L=4$, and the periapsis advance of parabolic orbits close to it. We compare the performance of our examples against the Paczyński-Wiita potential and other existing potentials. Finally, we discuss the limitations and advantages of our formulation.

Ultralight scalars can form superradiant clouds around rotating black holes. These may alter the dynamics of compact binaries and the ensuing waveform through orbital resonances and cloud ionization. We re-examine resonances involving states with nonzero decay width, deriving an effective treatment for resonances that are wider than the binary's frequency chirp. We demonstrate the utility of this approach by calculating an upper bound for the cloud's mass surviving up to the latest stages of the inspiral. Next, we study the accumulation of resonances with high-energy bound states. When these infinitely many, increasingly weak resonances are properly taken into account, they smooth out the "sharp features" in the binary's evolution that had been attributed to the ionization of the cloud. We compare our Newtonian results with recent relativistic calculations, highlighting common features as well as discrepancies. Our conclusions emphasize the need to carefully incorporate resonances in boson cloud waveform modeling.

A method is presented to use a fiber-optic device known as a photonic lantern to generate a reconfigurable custom wavefront for a null test of spherical, aspheric, and freeform optical surfaces. By modulating input intensity and phases at single mode fiber input ports, the wavefront of the output light field from the multimode end can be controlled to generate a custom nulling phase function. Generation of a desired wavefront is demonstrated by simulating a nineteen-port non-mode-selective photonic lantern. Using a linear response-matrix approach, a phase function with an RMS error of 44 nm from the target was generated in simulation. A compact form-factor noninterferometric null test for freeform optical surfaces is then described utilizing the photonic lantern as both a reconfigurable nulling source and a wavefront sensor.

We propose a new dark energy (DE) model from four parameter generalized entropy function of apparent horizon in a spatially flat universe. Such kind of generalized entropy is able to generalize all the known entropies proposed so far, for suitable representations of the entropic parameters. It turns out that the scenario can describe the correct thermal history of the universe, with the sequence of matter and dark energy epochs. Comparing with the $\Lambda$CDM model, the proposed generalized entropic DE model provides a higher value of present Hubble parameter for certain range of entropic parameter(s) leading to a possible resolution of Hubble tension issue. We confront the scenario with CC, PantheonPlus+SH0ES, DESI DR1 and compressed Planck likelihood datasets, which clearly depicts the phenomenological viability of the present model for some best fitted values of entropic parameter(s) that are indeed consistent with the resolution of Hubble tension.

We present a multi-scale optimal control framework for active seismic isolation in the Einstein Telescope, a third-generation gravitational-wave observatory. Our approach jointly optimizes feedback and blending filters in a cross-coupled opto-mechanical system using a unified cost function based on the "acausal optimum," which quantifies sensor signal-to-noise ratios across frequencies. This method enables efficient re-optimization under varying sensor configurations and environmental conditions. We apply the framework to two candidate sensing systems: OmniSens-a six-degree-of-freedom inertial isolation system-and BRS-T360, which combines Beam Rotation Sensor (BRS) as an inertial tilt sensor with T360 as a horizontal seismometer. We demonstrate superior low-frequency isolation with OmniSens, reducing platform motion by up to two orders of magnitude near the microseism. The framework allows for ready optimization and projection of sensor noise to metrics relevant to the performance of the instrument, aiding the design of the Einstein Telescope.

Accurate and reliable gravitational waveform models are crucial in determining the properties of compact binary mergers. In particular, next-generation gravitational-wave detectors will require more accurate waveforms to avoid biases in the analysis. In this work, we extend the recent NRTidalv3 model to account for higher-mode corrections in the tidal phase contributions for binary neutron star systems. The higher-mode, multipolar NRTidalv3 model is then attached to several binary-black-hole baselines, such as the phenomenological IMRPhenomXHM and IMRPhenomXPHM models, and the effective-one-body-based model SEOBNRv5HM_ROM. We test the performance and validity of the newly developed models by comparing them with numerical-relativity simulations and other tidal models. Finally, we employ them in parameter estimation analyses on simulated signals from both comparable-mass and high-mass-ratio systems, as well as on the gravitational-wave event GW170817, for which we find consistent results with respect to previous analyses.

The study of solar neutrinos presents significant opportunities in astrophysics, nuclear physics, and particle physics. However, the low-energy nature of these neutrinos introduces considerable challenges to isolate them from background events, requiring detectors with low-energy threshold, high spatial and energy resolutions, and low data rate. We present the study of solar neutrinos with a kiloton-scale liquid argon detector located underground, instrumented with a pixel readout using the Q-Pix technology. We explore the potential of using volume fiducialization, directional topological information, light signal coincidence and pulse-shape discrimination to enhance solar neutrino sensitivity. We find that discriminating neutrino signals below 5 MeV is very difficult. However, we show that these methods are useful for the detection of solar neutrinos when external backgrounds are sufficiently understood and when the detector is built using low-background techniques. When building a workable background model for this study, we identify {\gamma} background from the cavern walls and from capture of {\alpha} particles in radon decay chains as both critical to solar neutrino sensitivity and significantly underconstrained by existing measurements. Finally, we highlight that the main advantage of the use of Q-Pix for solar neutrino studies lies in its ability to enable the continuous readout of all low-energy events with minimal data rates and manageable storage for further offline analyses.

There is growing interest in leveraging LLMs to aid in astronomy and other scientific research, but benchmarks for LLM evaluation in general have not kept pace with the increasingly diverse ways that real people evaluate and use these models. In this study, we seek to improve evaluation procedures by building an understanding of how users evaluate LLMs. We focus on a particular use case: an LLM-powered retrieval-augmented generation bot for engaging with astronomical literature, which we deployed via Slack. Our inductive coding of 368 queries to the bot over four weeks and our follow-up interviews with 11 astronomers reveal how humans evaluated this system, including the types of questions asked and the criteria for judging responses. We synthesize our findings into concrete recommendations for building better benchmarks, which we then employ in constructing a sample benchmark for evaluating LLMs for astronomy. Overall, our work offers ways to improve LLM evaluation and ultimately usability, particularly for use in scientific research.

We present GravLensX, an innovative method for rendering black holes with gravitational lensing effects using neural networks. The methodology involves training neural networks to fit the spacetime around black holes and then employing these trained models to generate the path of light rays affected by gravitational lensing. This enables efficient and scalable simulations of black holes with optically thin accretion disks, significantly decreasing the time required for rendering compared to traditional methods. We validate our approach through extensive rendering of multiple black hole systems with superposed Kerr metric, demonstrating its capability to produce accurate visualizations with significantly $15\times$ reduced computational time. Our findings suggest that neural networks offer a promising alternative for rendering complex astrophysical phenomena, potentially paving a new path to astronomical visualization.