Papers for Thursday, Apr 03 2025

The interstellar medium (ISM) consists of a diversity of structures across a range of spatial scales, intimately tied to galactic evolution. In this work, Fourier analysis is used to characterize the spatial structures of dust and Polycyclic Aromatic Hydrocarbons (PAHs) in the ISM of PHANGS-JWST galaxies, observed in the four photometric mid-infrared (MIR) filters from F770W to F2100W (i.e., 7.7 to 21$\mu$m)). We quantify the abundance of structure on different spatial scales using the power-law slope, $\alpha$, of the spatial power spectra. The distribution of $\alpha$ across all length scales differs significantly between filters, with steeper slopes for PAH-dominated filters ($\alpha_{F770W} = 2.19^{+0.16}_{-0.15}$, $\alpha_{F1130W} = 1.88^{+0.25}_{-0.37}$) and shallower for the dust-continuum ($\alpha_{F1000W} = 1.48^{+0.33}_{-0.47}$, $\alpha_{F2100W} = 0.94^{+0.23}_{-0.28}$). The distribution of $\alpha$ across galaxies is narrower for PAH-dominated than for thermal dust-dominated bands, highlighting that PAHs trace photo-dissociation regions dominated by similar physical processes, whereas dust structures are an integrated property over the diverse evolutionary histories of their host galaxies. Unlike dust structures, PAH-sensitive bands display a break in the power spectrum: below a characteristic scale, $\ell_0=160\mathrm{pc}^{+110\mathrm{pc}}_{-50\mathrm{pc}}$, PAH structures are suppressed.

Galaxy and halo scaling relations, connecting a broad range of parameters, are well established from observations. The origin of many of these relations and their scatter is still a matter of debate. It remains a sizable challenge for models to simultaneously and self-consistently reproduce as many scaling relations as possible. We introduce the Magneticum Pathfinder hydrodynamical cosmological simulation suite, to date the suite that self-consistently covers the largest range in box volumes and resolutions. It is the only cosmological simulation suite that is tuned on the hot gas content of galaxy clusters instead of the stellar mass function. By assessing the successes and shortcomings of tuning to the hot gas component of galaxy clusters, we aim to further our understanding of the physical processes shaping the Universe. We analyze the importance of the hot and cold gas components for galaxy and structure evolution. We analyze 28 scaling relations, covering large-scale global parameters as well as internal properties for halos ranging from massive galaxy clusters down to galaxies, and show their predicted evolution from z=4 to z=0 in comparison with observations. These include the halo-to-stellar-mass and Kennicutt--Schmidt relations, the cosmic star formation rate density as well as the Fundamental Plane. Magneticum Pathfinder matches a remarkable number of the observed scaling relations from z=4 to z=0, including challenging relations like the number density of quiescent galaxies at cosmic dawn, the mass--size evolution, the mass--metallicity relation, the Magorrian relation, and the temperature--mass relation. We compile our data to allow for straightforward future comparisons. Galaxy properties and scaling relations arise naturally and the large scatter in observables at high redshift is crucial to distinguish the various galaxy formation models reproducing the z=0 relations.

Black hole (BH) accretion disks are often coupled to ultramagnetized and tenuous plasma coronae close to their central BHs. The coronal magnetic field can exchange energy between the disk and the BH, power X-ray emission, and lead to jetted outflows. Up until now, the coronal physics of BH accretion has only been studied using fluid modeling. We construct the first model of a BH feeding on a zero-net-flux accretion disk corona based on kinetic plasma physics. This allows us to self-consistently capture how collisionless relativistic magnetic reconnection regulates the coronal dynamics. We present global, axisymmetric, general relativistic particle-in-cell simulation of a BH coupled, via a series of magnetic loops, to a razor-thin accretion disk. We target the jet-launching regime where the loops are much larger than the BH. We ray-trace high-energy synchrotron lightcurves and track the flow of Poynting flux through the system, including along specific field-line bundles. Reconnection on field lines coupling the BH to the disk dominates the synchrotron output, regulates the flux threading the BH, and ultimately untethers magnetic loops from the disk, ejecting them via a magnetically striped Blandford-Znajek jet. The jet is initially Poynting-dominated, but reconnection operates at all radii, depleting the Poynting power logarithmically in radius. Coronal emission and jet launch are linked through reconnection in our model. This link might explain coincident X-ray flaring and radio-jet ejections observed during hard-to-soft X-ray binary state transitions. It also suggests that striped jet launch could be heralded by a bright coronal counterpart. Our synchrotron signatures resemble variability observed from the peculiar changing-look AGN, 1ES 1927+654, and from Sagittarius A*, hinting that processes similar to our model may be at work in these contexts.

Type Iax supernovae (SNe Iax) are thermonuclear explosions of white dwarfs, peculiar and underluminous compared to the normal type Ia supernovae. Observations of SNe Iax provide insight into the physics of white-dwarf explosions and suggest that some may not be terminal events. Late-time photometry ($\sim$1400 days post-peak) of the type Iax SN 2012Z, the only white dwarf supernova with a pre-explosion detection of a progenitor system, revealed a flux excess that may be explained by a gravitationally bound remnant driving a radioactively powered wind. We present further late-time Hubble Space Telescope photometry of SN 2012Z, $\sim$2500 days after the explosion, and find that the SN is still brighter than, but trending towards, the pre-explosion flux. Additionally, we observe that the excess F555W flux seen in previous data has grown more pronounced. The color of the excess flux disfavors a light echo or interaction with the circumstellar material. The decline rate of the excess flux is consistent with energy deposition from $^{55}$Fe, but the luminosity is higher than expected from models of the ejecta, further suggesting evidence for a bound remnant. Combining our data with future observations should allow for the detection of emission from the ejecta shock-heating of the companion helium star seen in the progenitor system.

Context. The Coronagraphic Instrument (CGI) on the Roman Space Telescope aims for unprecedented contrast for direct imaging of exoplanets, serving as a critical tech demo for future missions like the Habitable Worlds Observatory. This requires advanced wavefront sensing and control (WFS&C), including pair-wise (PW) probing for electric field estimation in the focal plane. Optimizing PW probe designs is vital to enhance performance and reduce overheads. Aims. We investigate different probe designs for PW probing in the context of Roman CGI. We compare classic sinc-sinc-sine probes, previously introduced single-actuator probes, and newly proposed sharp sinc probes in terms of effectiveness in focal-plane modulation, resilience to non-linearities, and overall impact on convergence and contrast. Methods. We conducted experiments on the THD2 testbed, configured to emulate Roman CGI with a custom Hybrid Lyot Coronagraph. We evaluated the three probe designs through WFS&C experiments using PW probing for estimation and electric field conjugation for wavefront correction. Simulations and hardware tests assessed contrast convergence and the impact of non-linear terms at varying probe amplitudes. We also explored low-flux scenarios to demonstrate the use of high-amplitude probes in reducing exposure times or closing the loop on faint targets. Results. Single-actuator probes emerged as the most effective, with faster convergence and reduced non-linear effects at high amplitudes. Sharp sinc probes performed moderately well but were less robust than single actuators. High-amplitude single-actuator probes showed advantages in dark-hole digging under low-flux, through faster iterations without significant degradation in contrast. The THD2 testbed, operating at contrasts analogous to Roman CGI, validated our results and underscored its role as a critical platform for advancing WFS&C techniques.

We present deep JWST/NIRCam and MIRI imaging of Ion1, a previously confirmed Lyman Continuum (LyC)-emitting galaxy at $z_{spec}=3.794$. Together with existing HST imaging, these new observations from the JADES program enable a joint analysis of Ion1's LyC, rest-frame UV, stellar, and dust emission with unprecedented detail. We report the first detection of dust emission at rest-frame $\sim3 \mu$m in a high-redshift LyC-emitting galaxy using MIRI/F1500W. Our analysis suggests a porous distribution of dust in Ion1, with regions exhibiting evidence of dust deficit coinciding both with LyC-emitting regions and with the peak of H$\alpha$ emission. Furthermore, multi-band NIRCam imaging reveals a strong FUV-to-optical color gradient, where LyC-emitting regions appear significantly bluer than the rest of Ion1. Spatially resolved SED fitting confirms that this color gradient is primarily driven by spatially varying dust attenuation. Together, these findings suggest that Ion1's LyC emission originates from a compact star-forming complex near its stellar-light centroid, where stellar feedback carves out low HI column density channels, facilitating LyC escape. However, only a fraction of these LyC photons - specifically those along sightlines with minimal HI obscuration - ultimately escape and reach observers. This work underscores the critical role of dust and neutral gas geometry in shaping LyC escape in galaxies at high redshifts. Anisotropic LyC escape may be a common feature in the early Universe, which must be properly incorporated to constrain the Epoch of Reionization.

Through gravitational lensing, galaxy clusters can magnify supernovae (SNe) and create multiple images of the same SN. This enables measurements of cosmological parameters, which will be increasingly important in light of upcoming telescopic surveys. We study the prospects of detecting strongly lensed SNe in cluster fields with the Nancy Grace Roman Space Telescope (Roman)'s High Latitude Time Domain Survey (HLTDS) and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). We employed two approaches: one focusing on known multiply imaged galaxies behind clusters, along with the SN rates specific to those galaxies, and another based on the expected number of lensed SNe exploding in a given volume behind a galaxy cluster. We collected all the clusters in the literature that feature a well-constrained lens model and multiply imaged galaxies behind clusters with high-quality data for the lensed galaxies. This allowed us to determine the SN rate for each galaxy. We provide predictions for 46 clusters visible to the Vera C. Rubin Observatory, as well as for 9 observable by Roman's HLTDS, depending on whether the clusters fall within the survey's observing field. We predict that the number of multiply imaged SNe discovered by LSST in its first three years is $3.95 \pm 0.89$ from the first approach or $4.94 \pm 1.02$ from the second. For the HLTDS, the expected number of multiply imaged SNe ranges from $0.38 \pm 0.15$ to $5.2 \pm 2.2$, depending on the specific cluster observed, however, the fields to be targeted remain a matter of discussion. We conclude that LSST offers great prospects for detecting multiply imaged SNe. Our predictions are effectively lower limits, as we only considered the most massive and well-studied clusters. We provide a recommendation for HLTDS observing field selection, namely: either MACS J0553.4-3342 or Abell 1758a should be observed by the survey.

Giant bulgeless disk galaxies, theoretically expected to be rare in the early Universe, have been confirmed by the James Webb Space Telescope (JWST) to exist as early as 2 billion years after the Big Bang. These morphologically extreme systems offer valuable insights into the physics of disk formation and the interplay between galaxies and their dark-matter halos. Using cosmological simulations, we identify analogs of such galaxies with stellar masses around $10^{11} M_\odot$ and half-light radii up to 6 kpc at $z \sim 3$ and characterize the factors that contribute to their formation. These galaxies form in young cosmic knots, populating host halos of high spin, low concentration, and spherical shapes. They feature dynamically coherent circum-galactic medium, as well as gas-rich, coherent mergers, which preserve their disk morphology and drive their large sizes. Interestingly, all the simulated giant disks harbor a compact, aligned inner disk, marginally resolvable in JWST images with a Sérsic index near unity. These findings highlight the environmental and structural conditions necessary for forming and sustaining giant bulgeless disks and provide a theoretical framework for interpreting JWST observations of extreme disk morphologies in the early Universe.

As powerful gamma-ray engines, blazars -- relativistic plasma jets launched toward Earth from active galactic nuclei -- are manifestly high-energy particle accelerators. Yet, exactly how these jets accelerate particles as well as what they are made of both remain largely mysterious. In this work, we argue that these issues may be linked through the gamma-ray emission for which blazars are renowned. Namely, high-energy photons produced at sites of intense particle acceleration could be absorbed by soft radiation within the jet, enriching it with electron-positron pairs. We explore this possibility in the specific context of particle acceleration by magnetized radiative relativistic turbulence. Using a combination of theory, particle-in-cell simulations, and Fokker-Planck modeling, we identify and characterize a novel pair-production-mediated equilibration mechanism in such turbulence. Initially, turbulent energy injection outpaces radiative cooling, leading to runaway particle acceleration and gamma-ray radiation. Then, gamma-ray absorption begets copious newborn pairs, slowing subsequent particle acceleration. This eventually brings particle acceleration into balance with radiative cooling and shuts down pair production: a pair-enriched final equilibrium. We estimate that this process could significantly load jets of flat-spectrum radio quasars with fresh pairs. These results represent an important connection between particle acceleration and plasma composition in blazar jets.

Pulsar Timing Array (PTA) experiments have the potential to unveil continuous gravitational wave (CGW) signals from individual massive black hole binaries (MBHBs). Detecting them in both gravitational waves (GW) and the electromagnetic (EM) spectrum will open a new chapter in multimessenger astronomy. We investigate the feasibility of conducting multimessenger studies by combining the CGW detections from an idealized 30-year SKA PTA and the optical data from the forthcoming LSST survey. To this end, we employed the $\texttt{L-Galaxies}$ semi-analytical model applied to the $\texttt{Millennium}$ simulation. We generated 200 different all-sky lightcones that include galaxies, massive black holes, and MBHBs whose emission is modeled based on their star formation histories and gas accretion physics. We predict an average of $\approx 33$ CGW detections, with signal-to-noise ratios $ S/N > 5$. The detected MBHBs are typically at $z < 0.5$, with masses of $ \sim 3 \times 10^{9} M_{\odot}$, mass ratios $> 0.6$ and eccentricities $\lesssim 0.2$. In terms of EM counterparts, we find less than 15% of these systems to be connected with an AGN detectable by LSST, while their host galaxies are easily detectable ($ < 23$ mag) massive ($ M_{\star} > 10^{11} M_{\odot}$) ellipticals with typical star formation rates ($10^{-15} yr^{-1} < sSRF < 10^{-10} yr^{-1}$). Although the CGW-EM counterpart association is complicated by poor sky localization (only 35% of these CGWs are localized within $\rm 100\, deg^2$), the number of galaxy host candidates can be considerably reduced (thousands to tens) by applying priors based on the galaxy-MBH correlations. However, picking the actual host among these candidates is highly non-trivial, as they occupy a similar region in any optical color-color diagram. Our findings highlight the considerable challenges entailed in opening the low-frequency multimessenger GW sky.

Galaxies are not isolated systems; they continuously interact with their surroundings by ejecting gas via stellar feedback and accreting gas from the environment. Understanding the interplay between outflows from the disc and the surrounding circumgalactic medium (CGM) is key to learning how star-forming galaxies evolve. Our goal is to understand how gas in the CGM is accreted onto the inner regions of the disc, making it available for the formation of stars, exploring the connection between stellar feedback and gas accretion from the CGM in Milky Way-like galaxies. We focus on the distribution of vertical and radial gas flows to and from the disc as a function of galactocentric radius, and examine the implications of these processes for the evolution of such galaxies. We use the Arepo code coupled with the SMUGGLE sub-grid model to perform hydrodynamic N-body simulations of 9 different galaxies surrounded by a hot CGM. Each simulation features a gaseous disc with different mass and scale length, allowing us to examine how disc structure impacts gas dynamics. We find evidence of a crucial link between stellar feedback and gas accretion from the CGM, which together play an essential role in sustaining ongoing star formation in the disc. In particular, the ejection of gas from the disc plane by stellar feedback leads to the generation of a baryon cycle in which the CGM gas is mainly accreted onto the external regions of the disc ($ \approx 3-10$ M$_{\odot}$ yr$^{-1}$ of gas is accreted into the whole disc). From these regions it is then transported to the centre with radial mass rates $\approx 1-4$ M$_{\odot}$ yr$^{-1}$ on average, owing to angular momentum conservation, forming new stars and starting the whole cycle again. We find that both vertical accretion onto the inner regions of the disc and the radial transport of gas from the disc outskirts are necessary to sustain star formation.

We present the weak-lensing mass calibration of 124 galaxy clusters and groups at redshift $0.1<z<0.8$ in the first Data Release of the $eROSITA$ All-Sky Survey ($eRASS1$) and constrain their BCG stellar-mass-to-halo-mass-and-redshift relation, using the data sets from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The cluster survey is conducted by the $eROSITA$ X-ray telescope aboard the Spectrum-Roentgen-Gamma space observatory. The cluster sample is X-ray-selected and optically confirmed with a negligibly low contamination rate. On a basis of individual clusters, the shear profiles of $96$ clusters are derived using the HSC Three-Year weak-lensing data, while the BCG stellar masses $M_{\star,\mathrm{BCG}}$ of $101$ clusters are estimated using the SED template fitting to the HSC five-band photometry. The count rate ($C_{\mathrm{R}}-M-z$) and BCG stellar mass ($M_{\star,\mathrm{BCG}}-M-z$) relations are simultaneously constrained in a population modelling. The $C_{\mathrm{R}}-M-z$ relation reveals a mass trend of $\propto M^{1.50^{+0.20}_{-0.30}}$, which is steeper than the self-similar prediction, and no deviation from the self-similar redshift scaling. We obtain the power-law indices of the mass ($B_{\mathrm{BCG}} = 0.25^{+0.11}_{-0.15}$) and redshift ($\gamma_{\mathrm{BCG}} = 0.87^{+0.92}_{-0.78}$) scaling for the $M_{\star,\mathrm{BCG}}-M-z$ relation, and find a strong degeneracy between them. By adopting an informative prior on $\gamma_{\mathrm{BCG}}$ to break the degeneracy, we obtain a statistically consistent $M_{\star,\mathrm{BCG}}-M-z$ relation. Our results suggest that the BCG stellar mass at a fixed halo mass has remained stable with a moderate increase since redshift $z\approx0.8$. This finding supports the picture of the ``rapid-then-slow'' BCG formation, where the majority of the stellar mass must have been assembled at much earlier cosmic time.

Bootstrap techniques relying on the constraints imposed by Extended Galilean Invariance (EGI), have proved to be very useful in the context of perturbation theory of the Large Scale Structure (LSS). It has been formulated in both the Eulerian as well as Lagrangian space. While the Eulerian bootstrap formalism has been successfully applied to both tracer and matter kernels, the application of bootstrap methods in Lagrangian space has so far been restricted to matter. Up to third order, it has been shown that implementing EGI constraints in Eulerian space fully reproduces the bias expansion for tracers. Previous studies have demonstrated that time non-locality affects the bias expansion in a non-trivial way starting from fifth order. Motivated by this fact, we extend the bootstrap approach upto fifth order in both Eulerian and Lagrangian space and demonstrate that it fully captures the time non-local effects. For this, we generalize the Lagrangian bootstrap for tracers, and found that it agrees with the corresponding results obtained in Eulerian space. One of the major challenges in implementing EGI constraints in Eulerian space, is to systematically find out all the "spurious poles" and make them vanish. We have proposed a method that bypasses this difficulty making the procedure tractable at higher orders. From Lagrangian perspective, we have identified coefficients in the tracer kernel whose ratios are independent of tracer properties and may serve as direct probes of the underlying cosmology.

A systematic search for HI 21-cm absorption in Quasar-Galaxy Pairs (QGPs) provides a powerful means to map the distribution of cold gas around high-redshift star-forming galaxies. Fiber spectroscopy of high-redshift quasars enables the serendipitous detection of foreground star-forming galaxies at extremely small impact parameters, forming a unique subset of QGPs known as Galaxies On Top Of Quasars (GOTOQs). In this study, we present results from a pilot upgraded Giant Metrewave Radio Telescope (uGMRT) survey of three GOTOQs, where we achieved a remarkable 100\% detection rate of HI 21-cm absorption. By combining our findings with existing literature, we establish that GOTOQs constitute a distinct population in terms of HI 21-cm absorption, with significantly higher detection rates than those observed in Damped Lyman-$\alpha$ (DLA)-based or metal absorption-based searches. For the GOTOQs, we find a strong correlation between the line-of-sight reddening and the HI 21-cm optical depth, characterized by $\int \tau\, dv\, (\rm{km\,s^{-1}}) = 13.58^{+2.75}_{-2.35} E(B-V) + 0.68^{+1.06}_{-1.27}$, consistent with the Milky Way sightlines. We also show that the HI 21-cm detection rates and optical depth declines rapidly with the impact parameter. With upcoming wide-field spectroscopic surveys expected to substantially expand the catalog of known GOTOQs, the success of this pilot survey lays the foundation for constructing a statistically significant sample of intervening HI 21-cm absorbers.

High resolution surveys reveal that the interstellar medium in the Milky Way and nearby galaxies consists of interlinked hierarchies of filamentary structure and superbubbles extending from galactic to subpc scales. The characterization of filament properties across this hierarchy is of fundamental importance for the origin of giant molecular clouds and their star clusters. In this paper we characterize the properties of filaments greater than 25 pc in length that are produced in the multi-scale galactic MHD simulations of Zhao et al. 2024. By adapting the FilFinder algorithm of Koch & Rosolowsky, 2015, we extract over 500 filaments ranging up to 10 kpc scales, to derive the probability distribution functions for filament masses and lengths, magnetic field orientations, and the gravitational stability and fragmentation patterns of filaments. We find power-law distributions for filament masses and lengths. The former has a power law index $\alpha_m = 1.85$ that is nearly identical to that of observed GMC mass functions in extragalactic and Galactic surveys, suggesting that GMC properties are inherited from their host filaments. The fragmentation of magnetized filaments on 200 pc scales or less occurs when they exceed an average critical line mass, as predicted by theory. On larger scales however, kpc filaments form out of the cold neutral medium (CNM) and fragmentation follows local variations in the critical line mass along spiral arms or at the boundaries of superbubbles.

Organosulfur species are potential major carriers of sulfur in the interstellar medium, as well as interesting ingredients in prebiotic chemistry. The most fundamental question regarding these species is under which conditions they reside in the gas versus solid phase. Here, we characterize the thermal desorption kinetics, binding energies, and entrapment of the organosulfur methyl mercaptan (CH$_3$SH, or MeSH) in different ice environments, comparing them with those of methanol (CH$_3$OH, or MeOH) ices. The derived multi-layer (pure MeSH-MeSH) and sub-monolayer (layered MeSH-H$_2$O) binding energies are surprisingly similar, corresponding to snow line locations where the disk midplane temperature is ~105 K. In both H$_2$O-dominated and more realistic H$_2$O:CO$_2$-dominated ices, 100% of the MeSH is entrapped, almost exclusively desorbing at the molecular volcano desorption peak, indicating that MeSH is retained at the water snow line if initially mixed with water ice during formation. Additionally, the presence of MeSH in an ice mixture enhances the entrapment of CO$_2$ and MeOH (up to 100%) until the onset of volcano desorption; without MeSH, both desorb at their respective pure desorption temperatures and also co-desorb with water. Compared to MeOH, MeSH binds less well to water, explaining why MeSH escapes during water ice crystallization rather than co-desorbing with water. These results show the larger relative size of MeSH compared to MeOH significantly impacts its ability to bind to water and its entrapment efficiency. Therefore, molecular size plays an important role in the adsorption and retention of S-bearing organics and, in turn, other volatiles in ices.

We investigate the capabilities of the mid-infrared instrument (MIRI) of James Webb Space Telescope (JWST) to advance our knowledge of AGN dust using the spectral fitting technique on an AGN collection of 21 nearby (z<0.05) AGN (7 type-1 and 14 type-2) observations obtained with the medium resolution spectroscopy (MRS) mode. This collection includes publicly available AGN and data from the collaboration of Galactic Activity, Torus, and Outflow Survey (GATOS). We developed a tool named MRSPSFisol that decomposes MRS cubes into point-like and extended contributions. We found statistically good fits for 12 targets with current AGN dust models. The model that provides good fits (chi2/dof<2) for {these 12 targets} assumes a combination of clumpy and smooth distribution of dust in a flare-disk geometry where the dust grain size is a free parameter. Still, two and one AGN statistically prefer the disk+wind and the classical clumpy torus model, respectively. However, the currently available models fail to reproduce 40% of the targets, likely due to the extreme silicate features not well reproduced by the models and signatures of water-ice and aliphatic hydrocarbon absorption features in most targets. New models exploring, for instance, new chemistry, are needed to explain the complexity of AGN dust continuum emission observed by JWST.

Eclipsing binary stars with Delta Scuti components are exceptional systems for gaining a deeper understanding of stellar systems. WY Leo is one such system, exhibiting Delta Scuti-type oscillations. However, it has not been extensively studied in the literature. WY Leo was observed by Transiting Exoplanet Survey Satellite (TESS), providing two sectors of high-quality photometric data. This study focuses on the photometric analysis of WY Leo. Binary modeling was performed, and the system's fundamental stellar parameters, such as mass and radius, were determined. The pulsational properties of WY Leo were also investigated, revealing that the more luminous star exhibits Delta Scuti type oscillations with a period of 0.052 days. The position of the primary component was examined on the HR diagram, and it was found to lie within the Delta Scuti instability strip.

Polarimetric light curves of Sagittarius A* (Sgr A*) sometimes exhibit loops in the Stokes Q and U plane over time, often interpreted as orbiting hotspot motion. In this work, we apply the differential geometry of planar curves to develop a new technique for estimating polarimetric rotation rates. Applying this technique to 230 GHz light curves of Sgr A*, we find evidence of clockwise motion not only during a post-flare period on 2017 April 11th, as previously discovered, but also during the quiescent days imaged by the Event Horizon Telescope (EHT). The data exhibit a clockwise fraction of $0.67 \pm 0.09$ and an overall Q-U rotation rate of $-3.7 \pm 1.2 \ \mathrm{deg}\,t_g^{-1}$. We analyze a library of General Relativistic Magnetohydrodynamic (GRMHD) simulations and find that face-on, clockwise-rotating models with strong magnetic fields are most likely to be consistent with the observations. These results are consistent with EHT and GRAVITY Collaboration studies, and indirectly support an interpretation in which the polarized image of Sagittarius A* has been rotated by an external Faraday screen. This technique offers a novel probe of event horizon scale dynamics that complements dynamical reconstructions.

The rare association of three persistent radio sources (confirmed PRS1 and PRS2, candidate PRS3) with repeating fast radio bursts (FRB 20121102A, 20190520B, 20201124A) offers a unique probe into their magneto-ionic environments. PRSs are attributed to synchrotron emission from relativistic charged particles of magnetar wind nebula (MWN) powered by spin-down magnetohydrodynamic wind or internal magnetic field decay. Using a multizone hydrodynamic model, we track MWN evolution to constrain magnetar progenitor properties. For PRS1 and PRS2, we find an equipartition radius $R_\mathrm{eq} \sim 0.1 $ pc that is consistent with the radio scintillation estimates ($> 0.03$ pc) and radio imaging limits ($<0.7$ pc). This compact size favors low expansion speeds and large initial spin periods, $P_\mathrm{i} \gtrsim 10$ ms, ruling out millisecond magnetar progenitors. Given $P_\mathrm{i} \gtrsim 10$ ms, a current size of $\sim 0.1$ pc, a supernova kinetic energy $E_\mathrm{ SN} \sim 10^{50}-10^{51}$ erg and an ejecta mass $M \sim 3-10 \; M_{\odot}$, the PRS age is $t \sim 10-10^{2}$ yr. PRSs with $t>20$ years require an internal field ($B_\mathrm{int} \sim 10^{16}-10^{16.5}\;$G) with a decay timescale $t_\mathrm{d} \sim 10-10^{2.5}$ yr. The slowest field decay ($t_\mathrm{d,max} \sim 500$ yr) favors sub-energetic supernovae ($E_\mathrm{SN} \sim 10^{50}$ erg) with massive ejecta ($M \gtrsim 10\; M_{\odot}$) and low ionization fraction ($\sim 3\% $). For the sub-energetic scenario for the confirmed PRSs, we predict a cooling break at $100-150$ gigahertz at $20-40 \; \mu \text{Jy}$ and self-absorption near 200 megahertz at $180\; \mu \text{Jy}$. For PRS3, a rotation-powered MWN is viable only if $t \sim 10$ yr; an inverted spectrum beyond 150 gigahertz would rule out this scenario.

Galaxy groups are essential for studying the distribution of matter on a large scale in redshift surveys and for deciphering the link between galaxy traits and their associated halos. In this work, we propose a widely applicable method for identifying groups through machine learning techniques in real space taking into account the impact of redshift distortion. Our methodology involves two neural networks: one is a classification model for identifying central galaxy groups, and the other is a regression model for predicting the mass of these groups. Both models input observable galaxy traits, allowing future applicability to real survey data. Testing on simulated datasets indicates our method accurately identifies over $92\%$ of groups with $\mathrm{M}_{vir} \geq 10^{11}h^{-1}\mathrm{M}_\odot$, with $80\%$ achieving a membership completeness of at least $80\%$. The predicted group masses vary by less than 0.3 dex across different mass scales, even in the absence of a priori data. Our network adapts seamlessly to expand to sparse samples with a flux limit of $m_{r} < 14$, to high redshift samples at $z=1.08$, and to galaxy samples from the TNG300 hydrodynamical simulation without further training. Furthermore, the framework can easily adjust to real surveys by training on redshift distorted samples without needing parameter changes. Careful consideration of different observational effects in redshift space makes it promising that this method will be applicable to real galaxy surveys.

Astronomy, a captivating field that draws upon science, mathematics, and engineering, has traditionally relied on visual representations to convey the wonders of the cosmos. While this approach effectively engages the sighted population, the use of imagery can exclude individuals with blindness or visual impairment (B/VI). Astronomical research is incorporated into press releases, media, outreach efforts, and educational systems aimed at enhancing public interest and often skill in science, but visual materials can hamper a population with B/VI. This paper explores the potential of 3D printing as an assistive technology providing an alternative to imagery. We produced textured 3D prints of astronomical research data from the Hubble Space Telescope (HST). Useability assessment of materials is an important phase of production before integration into structured programs, and we used a multi-phased approach in our prior research to create and test appropriate textures for 3D astronomical prints. This paper describes the last step of reviewing our 3D prints through informal useability sessions with diverse individuals. The assessment indicated our 3D prints provide reliable, informative representations of astronomical data appropriate for public use especially for public information, outreach programs, and science education for individuals with BVI.

Detecting and characterizing the atmospheres of terrestrial exoplanets is a key goal of exoplanetary astronomy, one that may now be within reach given the upcoming campaign to conduct a large-scale survey of rocky M-dwarf worlds with the James Webb Space Telescope. It is imperative that we understand where known planets sit relative to the cosmic shoreline, the boundary between planets that have retained atmospheres and those that have not. Previous works modeled the historic XUV radiation received by mid-to-late M-dwarf planets using a scaling relation calibrated using more massive stars, but fully convective M dwarfs display unique rotation/activity histories that differ from Sun-like stars and early M dwarfs. We synthesize observations of the active lifetimes of mid-to-late M dwarfs to present an updated estimate of their historic XUV fluence. For known planets of inactive, mid-to-late M dwarfs, we calculate a historic XUV fluence that is 2.1-3.1 times the canonical XUV scaling relation on average, with the larger value including corrections for the pre-main-sequence phase and energetic flares. We find that only the largest terrestrial planets known to orbit mid-to-late M-dwarfs are likely to have retained atmospheres within the cosmic shoreline paradigm. Our calculations may help to guide the selection of targets for JWST and may prove useful in interpreting the results; to this end, we define a novel Atmosphere Retention Metric (ARM) that indicates the distance between a planet and the cosmic shoreline, and tabulate the ARM for known mid-to-late M-dwarf planets.

Using published simulations of the 10-year Legacy Survey of Space and Time (LSST), we forecast its ability to determine the masses of individual main-belt asteroids (MBAs) through precise astrometry of any pairs of the ~1.2 million known MBAs undergoing close gravitational encounters during the survey. The uncertainty $\sigma_I$ on the impulse applied to a tracer asteroid by its deflector is derived from the Fisher matrix of the tracer's astrometric data, including an azimuthal acceleration $A_2$ from the Yarkovsky effect as a free parameter for each tracer. If only LSST observations are available, $\sigma_I \approx7\times10^{-6}\,{\rm m}\,{\rm s}^{-1}$ for MBAs at apparent magnitude $m_V<19.5,$ degrading ~10x for $m_V=23.$ These tracers yield a median uncertainty on the mass of an MBA of $\approx4\times10^{-14}M_\odot,$ with a wide range of variation depending on the ``luck'' of close encounters. Roughly 125 MBAs obtain mass measures with S/N>5. If pre-LSST astrometry yields a strong constraint on the state vector of the tracer MBA at the start of LSST, then these values improve to median $\sigma_M\approx1.3\times10^{-14}M_\odot$ and 310 MBAs at S/N>5, with >1/2 of these having S/N>10. These yields would be a ~10-fold increase in the number of known asteroid masses, including a nearly complete knowledge of MBAs with H<7.5. If pre-LSST data are sufficient to start constraining the Yarkovsky effect, another factor $\sim1.5$ can be gained. Tables of the measurable deflector MBAs and their tracers are provided.

The planetary nebula NGC 6720, also known as the ``Ring Nebula", is one of the most iconic examples of nearby planetary nebulae whose morphologies present a challenge to our theoretical understanding of the processes that govern the deaths of most stars in the Universe that evolve on a Hubble time. We present new imaging with JWST of the central star of this planetary nebula (CSPN) and its close vicinity, in the near- to mid-IR wavelength range. We find the presence of a dust cloud around the CSPN, both from the spectral energy distribution at wavelengths >~5 micron, as well as radially-extended emission in the 7.7, 10 and 11.3 micron images. From modeling of these data, we infer that the CSPN has a luminosity of 310 Lsun and is surrounded by a dust cloud with a size of ~2600 au, consisting of relatively small amorphous silicate dust grains (radius ~0.01 micron) with a total mass of 1.9 x 10^(-6) M_earth. However, our best-fit model shows a significant lack of extended emission at 7.7 micron -- we show that such emission can arise from a smaller (7.3 x 10^(-7) M_earth) but uncertain mass of (stochastically-heated) ionized PAHs. However, the same energetic radiation also rapidly destroys PAH molecules, suggesting that these are most likely being continuously replenished, via the outgassing of cometary bodies and/or the collisional grinding of planetesimals. We also find significant photometric variability of the central source that could be due to the presence of a close dwarf companion of mass < ~0.1 Msun.

Mass loss is a key aspect of stellar evolution, particularly in evolved massive stars, yet episodic mass loss remains poorly understood. To investigate this, we need evolved massive stellar populations across various galactic environments. However, spectral classifications are challenging to obtain in large numbers, especially for distant galaxies. We addressed this by leveraging machine-learning techniques. We combined \textit{Spitzer} photometry and Pan-STARRS1 optical data to classify point sources in 26 galaxies within 5 Mpc, and a metallicity range 0.07-1.36 Z$_\odot$. \textit{Gaia} DR3 astrometry was used to remove foreground sources. Classifications are derived using a machine-learning model developed by Maravelias et al. (2022). We report classifications for 1,147,650 sources, with 276,657 sources ($\sim24\%$) being robust. Among these are 120,479 Red Supergiants (RSGs; $\sim11\%$). The classifier performs well even at low metallicities ($\sim0.1$ Z$_\odot$) and distances under 1.5 Mpc, with a slight decrease in accuracy beyond $\sim3$ Mpc due to \textit{Spitzer}'s resolution limits. We also identified 21 luminous RSGs ($\textrm{log}(L/L_\odot)\ge5.5$), 159 dusty Yellow Hypergiants in M31 and M33, as well as 6 extreme RSGs ($\textrm{log}(L/L_\odot)\ge6$) in M31, challenging observed luminosity limits. Class trends with metallicity align with expectations, though biases exist. This catalog serves as a valuable resource for individual-object studies and \textit{James Webb} Space Telescope target selection. It enables follow-up on luminous RSGs and Yellow Hypergiants to refine our understanding of their evolutionary pathways. Additionally, we provide the largest spectroscopically confirmed catalog of massive stars and candidates to date, comprising 5,273 sources (including $\sim330$ other objects).

We investigate the influence of the cosmic web on galaxy properties in the IllustrisTNG simulations. To disentangle the effects of galaxy groups and cosmic filaments, we classify the cosmic web environment into four categories: group, group-dominated, filament-dominated, and field. By controlling for stellar mass, we reveal evident differences in specific star formation rates (sSFR), quenched fraction, gas fractions, local density, and stellar ages among central galaxies in different cosmic web environments, particularly for lower-mass galaxies. However, these differences largely diminish when the effect of local overdensity is further accounted for, indicating its dominant role. Additionally, we observe distinct differences in these properties among satellite galaxies across environments, mainly driven by stellar mass, halo mass, and overdensity. Notably, residual differences between satellites in field and filament-dominated region persist even after controlling for these factors, suggesting a stronger susceptibility of satellite galaxies to filaments compared to centrals. Our findings highlight the importance of differentiating between central and satellite to accurately assess the environmental effects of the cosmic web. Our analysis suggests that the relationship between galaxy properties and their distance from filaments arises from a combination of factors, including stellar and halo mass, groups, overdensity, and the intrinsic influence of the cosmic web. Additionally, we find that the effect of the cosmic web on galaxy properties is reduced at $z=0.5$, compared to $z=0$. Furthermore, central galaxies near thick filaments tend to exhibit slightly to moderately lower sSFR and cold gas fractions compared to those near thin filaments.

Quasars are effective tracers of the large-scale distribution of galaxies at high redshift thanks to their high luminosity and dedicated surveys. Previous studies have shown that quasars exhibit a bias similar to that of rich groups, indicating that quasar pairs may be linked to higher density environments, serving as protocluster proxies. In this work, we aim to characterize close quasar pairs within the same halo by identifying them in redshift space. We analyze pair-quasar cross correlations and CMB-derived lensing convergence profiles centered in these systems. We have identified 2777 pairs of quasars in the redshift range from 1.2 to 2.8. Quasar pairs exhibit a distribution of relative luminosities that differs from that of random pairs with the same redshift distribution, indicating that quasars in these systems are distinct from isolated ones. The cross-correlation between pairs and quasars shows a larger amplitude than the auto-correlation function of quasars, suggesting these systems are more strongly biased toward the large-scale mass distribution and reside in more massive halos. This is further supported by higher convergence CMB lensing profiles of pairs compared to isolated quasars with a similar redshift distribution. Our results indicate that quasar pairs are suitable precursors to present-day clusters of galaxies, unlike isolated quasars, which are associated with less dense environments.

The discovery of planets around PSR~1257+12 suggests that planetary systems may be detected around the recycled pulsars found in globular clusters. Planetary systems in dense clusters have lifetimes to disruption due to perturbations by passing stars comparable to or shorter than the pulsar lifetime, and observations of planets in the cores of clusters may reveal planetary systems formally dynamically unstable on time scales short compared to the characteristic age, $\tau_c$, of the system. Planets formed around cluster pulsars will most likely be restricted to semi--major axis of $\sim 0.1-1.0\, AU$, while "scavenged" planets may be observed in wider orbits, with no stable systems expected in the densest clusters. Observation is most probable in the cluster rich high--density pre-core collapse clusters such as 47Tuc. Detection of planets around cluster pulsars can constrain planet formation mechanisms, in particular the effects of low metallicity on planet formation in disks.

One of the crucial components in simulating the growth and evolution of galaxies within a cosmological framework is the modeling of star formation (SF) and its corresponding feedback. Traditionally, the implemented SF law follows the empirical Kennicutt-Schmidt relation, which links the SF rate (SFR) in a gas element to the total gas density. More recently, an even stronger correlation has been observed between the SFR and the amount of molecular hydrogen ($\mathrm{H}_2$). This opens up the question of whether molecular hydrogen is a necessary precursor for SF or is instead a tracer of the total amount of gas, both of which would explain the observed correlations. In this study, we examine the impact of using an $\mathrm{H}_2$-based SF law on the formation of Milky Way-mass galaxies using cosmological, hydrodynamical simulations. We find that the galaxy modeled with our approach exhibits a well-defined, disc-like morphology. Compared to the traditional recipe, our model delays the onset of SF by approximately $500 \, \mathrm{Myr}$ resulting in a lower SFR, a smaller disc size, and a higher proportion of neutral to ionized gas within the disc region. These findings highlight the importance of including sophisticated SF models which can be compared to several observations -- including those related to $\mathrm{H}_2$ -- to better understand the processes affecting galaxy formation and evolution.

Studying the atomic-to-molecular transition is essential for understanding the evolution of interstellar medium. The linear edge of Taurus molecular cloud, clearly identified in the $^{13}$CO(1-0) intensity map, serves as an ideal site for investigating this transition. Utilizing the Arizona Radio Observatory Sub-Millimeter Telescope, we obtained mapping observations of CO(2-1), $^{13}$CO(2-1), and CO(3-2) across this linear edge. The intensity ratio between CO(2-1) and $^{13}$CO(2-1) indicates a lower limit of the $^{12}C/^{13}C$ ratio of $54\pm 17$. Based on multi-transition observations of CO and $^{13}$CO, we performed Markov Chain Monte Carlo (MCMC) fit of the physical properties across this edge using non-Local Thermodynamic Equilibrium analysis with the RADEX code, based on the Large velocity Gradient (LVG) assumption. The number density profile exhibits a pronounced jump coinciding with the H$_2$ infrared emission peak. The cold HI gas within the molecular cloud, manifested as HI-Narrow Self-Absorption (HINSA) features, is detected along the cloud edge. Our quantitative comparison with numerical simulations provides tentative evidence that shocks induced by colliding gas flows may contribute to the atomic-to-molecular phase transition observed along the linear edge.

Time-delay distance measurements from strongly lensed quasars provide a robust, independent method for determining the Hubble constant ($H_0$). This approach cross-checks $H_0$ estimates from the distance ladder in the late universe and the cosmic microwave background in the early universe. However, the mass-sheet degeneracy in lensing models introduces systematic uncertainty, limiting precision. Dynamical modeling complements strong lensing by constraining the mass distribution with independent observational data. We develop a methodology and software framework for joint modeling of stellar kinematics and lensing data. Using simulated data for the lensed quasar RXJ1131$-$1131, we demonstrate that high-quality kinematic data can achieve $\sim$4% precision on $H_0$. Through extensive modeling, we examine the impact of the presence of a supermassive black hole in the lens galaxy and potential systematic biases in kinematic data on $H_0$ measurements. Our results show that imposing priors on black hole mass and orbital anisotropy, or excluding central kinematic bins, mitigates biases in $H_0$ estimates. By testing on mock kinematic data with systematic biases, we highlight the need for sub-percent control of kinematic systematics, which is achievable with current technology. Additionally, we leverage GPU parallelization to accelerate Bayesian inference, reducing a previously month-long process by an order of magnitude. This pipeline offers significant potential for advancing cosmological and galaxy evolution studies with large datasets.

The James Webb Space Telescope (JWST) is challenging our understanding of the nature of the very first galaxies in the Universe, having discovered a surprising abundance of very massive galaxies in early cosmic epochs. By applying a model of primeval dust, we estimate the far-infrared (FIR)-continuum luminosities for galaxies of different masses at redshifts $z\gtrsim7$. In particular, we predict observed fluxes for different available bands (3-10) of the Atacama Large Millimeter/sub-mm Array (ALMA), considering typical conservative values expected for the properties of first galaxies (e.g., gas-phase metallicities, dust-to-metal ratio, star formation efficiency). As expected, FIR fluxes increase with stellar mass for all ALMA bands, but with a steeper slope for bands 9 and 10. Encouragingly, first galaxies are affected by a strong negative-K correction, in such a way that sources with similar properties are brighter in bands 3-8 at higher redshifts. Such behavior is stronger for bands 6-7. Although the trends for bands 9-10 are not clear, the highest fluxes for such bands are reached towards extreme $z\gtrsim 15$. Counterintuitively, our results suggest that JWST sources with similar masses and dust properties would be more easily detectable with ALMA if they are located at higher $z$.

The thermal evolution of magma oceans formed by giant impacts is strongly influenced by a proto-atmosphere through its blanketing effect, which suppresses outgoing planetary radiation. While both radiative absorption and Rayleigh scattering by atmospheric species can contribute to this effect, the role of the scattering in suppressing thermal radiation from magma oceans remains unclear. In this study, we developed a 1-D radiative transfer model for planetary and solar radiation in a proto-atmosphere composed of H2O and H2, and a coupled thermal evolution model of a planetary interior and proto-atmosphere, to investigate the scattering blanketing effect on planetary radiation and magma ocean cooling. Our results show that Rayleigh scattering significantly reduces outgoing planetary radiation at wavelengths below ~1 micrometer, particularly in hot, thick atmospheres where scattering is highly effective. Consequently, the planetary outgoing radiation flux decreases by up to about one to two orders of magnitude, and the magma ocean lifetime is prolonged by up to about three times due to the scattering blanketing effect when the total amounts of H2O and H2 are equivalent to or greater than the present-day terrestrial seawater. These findings suggest that the prolonged magma ocean phase facilitated efficient differentiation between compatible and incompatible elements, even in the lower mantle. Furthermore, they imply that sustained magma oceans likely persisted throughout much of the giant impact phase, supporting a magma ocean origin of the Moon consistent with its observed chemical characteristics.

We investigate the intrinsic and observational properties of $z\gtrsim 6$ galaxies hosting coalescing massive black holes (MBHs) that gives rise to gravitational waves (GWs) detectable with the Laser Interferometer Space Antenna (LISA). We adopt a zoom-in cosmological hydrodynamical simulation of galaxy formation and black hole (BH) co-evolution, zoomed-in on a $M_h \sim 10^{12}~\rm M_{\odot}$ dark matter halo at z = 6, which hosts a fast accreting super-massive black hole (SMBH) and a star-forming galaxy. Following the SMBH formation backward in time, we identify the merging events that concurred to its formation and we pick up the ones that are detectable with LISA. Among these LISA detectable events (LDEs), we select those that, based on their intrinsic properties are expected to be bright in one or more electromagnetic (EM) bands. We post-process these events with dust radiative transfer calculations to make predictions about their spectral energy distributions and continuum maps in the JWST to ALMA wavelength range. We compare the spectra arising from galaxies hosting the merging MBHs with those arising from AGN powered by single accreting BHs. We find that it will be impossible to identify an LDE from the continuum SEDs because of the absence of specific imprints from the merging MBHs. We also compute the profile of the H$_{\rm \alpha}$ line arising from LDEs, considering the contribution from their star-forming regions and the accreting MBHs. We find that the presence of two accreting MBHs would be difficult to infer even if both MBHs accrete at super-Eddington rates. We conclude that the combined detection of GW and EM signals from $z\gtrsim 6$ MBHs is challenging not only because of the poor sky-localization provided by LISA, but also because the loudest GW emitters are not massive enough to leave significant signatures in the emission lines arising from the broad line region.

Astrometry plays a crucial role in understanding the structure, dynamics, and evolution of celestial objects by providing precise measurements of their positions and motions. We propose a new approach to wide-field, relative astrometry, Plate Analysis. Plate Analysis is an innovative algorithm that estimates stellar coordinates and corrects geometric distortion based on precise reference sources, without relying on dedicated calibration fields. It is implemented as a probabilistic framework using Stochastic Variational Inference to efficiently optimize the numerous parameters involved in the model. The methodology was tested through a simplified simulation designed after the Galactic center survey of the JASMINE mission. This simulation, called the JASMINE mini-mock survey, covered three years of observation with 100 satellite orbits, providing a comprehensive dataset for evaluating the performance of Plate Analysis. Although the observation model incorporated more than 30,000 parameters, the parameters were efficiently optimized through Plate Analysis. The results showed that the estimated coordinates closely match the expected stellar motions, with an average positional error of about $70\,\mathrm{\mu as}$. These findings validate the potential of Plate Analysis for precise wide-field astrometric applications, offering significant insights into improving the accuracy of stellar dynamics measurements.

We analyse high resolution optical spectra of MWC148 (optical counterpart of the gamma-ray source HESS J0632+057) obtained at Observatoire de Haute Provence and Rozhen Observatory. We measure equivalent widths of 7 diffuse interstellar bands and estimate the interstellar extinction E_{B-V}=0.85 +/- 0.08.

The interstellar medium (ISM) is the material that fills the space between the stars in all galaxies; it is a multi-phase medium in pressure equilibrium, with densities and temperatures covering over 6 orders of magnitude. Although accounting for only a small fraction of the mass of any given galaxy, it is a vital component, since it holds the material responsible for galaxy growth through star formation. Studying the ISM requires careful observations at all wavelengths of the electromagnetic spectrum. This article describes the multi-phase nature of the ISM, and then puts it in the context of galaxy evolution models, emphasising the importance of the cycling of baryons in and out of galaxies. Within this framework, the ISM plays a central role: it connects the physical processes operating on very large physical- and time-scales which control the accretion of gas onto galaxies, and the small scale processes that regulate star formation.

The nature of progenitors of Type Ia supernovae (SNe Ia) and their explosion mechanism remain unclear. It has been suggested that SNe Ia may be resulted from thermonuclear explosions of hybrid carbon-oxygen-neon white dwarfs(CONe WDs) when they grow in mass to approach the Chandrasekhar mass limit by accreting matter from a binary main-sequence (MS) companion. In this work, we combine the results of detailed binary evolution calculations with population synthesis models to investigate the rates and delay times of SNe Ia in the CONe WD + MS channel at low metallicity environment of Z = 0.0001. For a constant star formation rate of 5 M_sun yr-1, our calculations predict that the SN Ia rates in the CONe WD + MS channel at low metallicity of Z = 0.0001 is about 0.11 - 3.89 * 10-4 yr-1. In addition,delay times in this channel cover a wide range of 0.05 - 2.5 Gyr. We further compare our results to those given by previous study for the CONe WD + MS channel with higher metallicity of Z = 0.02 to explore the influence of metallicity on the results. We find that these two metallicity environments give a slight difference in rates and delay times of SNe Ia from the CONe WD + MS channel, although SNe Ia produced at low metallicity environment of Z = 0.0001 have relatively longer delay times.

In various types of supernovae (SNe), strong interaction between the SN ejecta and circumstellar material (CSM) has been reported. This raises questions on their progenitors and mass-loss processes shortly before the explosion. Recently, the bright long-lived Type~II SN 2021irp was proposed to be a standard Type II SN interacting with disk-like CSM. The observational properties suggest that the progenitor was a massive star in a binary system and underwent a mass-ejection process due to the binary interaction just before the explosion. Here, we study the diversity of the observational properties of bright long-lived Type II (21irp-like) SNe. We analyse the diversity of their CSM properties, in order to understand their progenitors and mass-loss mechanisms and their relations with the other types of interacting SNe. We performed photometry, spectroscopy, and/or polarimetry for four 21irp-like SNe. Based on these observations as well as published data of SN~2021irp itself and well-observed bright and long-lived type II SNe including SNe~2010jl, 2015da and 2017hcc, we discuss their CSM characteristics. This sample of SNe shows luminous and long-lived photometric evolution, with some variations in the photometric evolution (from $\sim-17$ to $\sim-20$ absolute mag in the $r$/$o$ band even at $\sim 200$ days after the explosion). They show photospheric spectra characterized mainly by Balmer lines for several hundreds of days, with some variations in the shapes of the lines. They show high polarization with slight variations in the polarization degrees with rapid declines with time (from $\sim3-6$ \% before the peak to $\sim1$ \% at $\sim200$ days after the peak). The observational properties are consistent with the disk-CSM-interaction scenario, i.e., typical Type~II SNe interacting with disk-like CSM.

Warming Mars' surface could be a step toward making it suitable for life, but would represent a major science and engineering challenge. To warm Mars using engineered aerosol, particles released locally must disperse globally. The winds that transport aerosol respond to the aerosol's radiative forcing, implying strong radiative-dynamical feedbacks. Using a plume-tracking climate model without a water cycle, we investigate radiative-dynamical feedback from surface release of two particle compositions: graphene (which attenuates UV) and aluminum. Both compositions can give fast global warming of ~30 K. We infer that 2 liters/second release rate of graphene made from Mars' atmosphere via CO2-electrolysis could double Mars' greenhouse effect (+5 K). Self-lofting helps particles rise and spread. The Hadley cell strengthens under warming, aiding mixing. Warming can be focused in latitude by tuning particle size. Within our model, Mars radiative-dynamical feedbacks enable engineered-aerosol warming. Challenges remain, including functionalization, agglomeration, dry-deposition experiments, and modeling water cycle feedbacks.

In this work an analysis of the intracluster light (ICL) in the galaxy cluster SMACS J0723.3-7327 (hereafter, SMACS J0723) using JWST/NIRCam deep imaging in six filters (F090W to F444W) is presented. The images were processed for low surface brightness (LSB) science, with additional correction for instrumental scattering in the short-wavelength channels, and analysed using wavelet-based decomposition. The ICL, brightest cluster galaxy (BCG), and satellite galaxies were extracted and modelled, with 2D maps for each component. ICL and ICL+BCG fractions, computed across all filters within a 400 kpc radius, exhibit a flat trend with wavelength, averaging 28% and 34%, respectively. Flux ratios between the BCG and the next brightest members (M$_{12}$, M$_{13}$ and M$_{14}$) also display minimal wavelength dependence. These results indicate that SMACS J0723 is a dynamically evolved cluster with a dominant BCG and well-developed ICL. Five prominent ICL substructures are analysed, contributing to 10-12% of the total ICL+BCG flux budget, slightly exceeding simulation predictions. Their short dynamical timescales suggest an instantaneous ICL injection rate of several $10^3 L_{\odot}\,{\rm yr}^{-1}$, consistent with active dynamical assembly. These findings support a scenario where SMACS J0723's ICL growth is currently driven by galaxy mergers involving the BCG and other bright satellites, rather than by the accretion of pre-processed ICL from a recent cluster merger. However, extrapolating the current injection rate to the cluster's lifetime indicates that additional mechanisms are required to match the growth observed in other clusters over cosmic times.

We report our photometric and spectroscopic observations and analysis of the 2024 eclipse of the symbiotic binary V1413 Aquilae. We found the system in a visually bright state and the eclipse time of minimum consistent with the published ephemeris. The eclipse profile showed that the hot component was an extended object rather than an isolated white dwarf. By analyzing the eclipse profile we estimated the orbital inclination to be 67.9°, the radius of the extended hot component surrounding the white dwarf to be 39.3 Rsun, and that the red giant star was probably filling its Roche Lobe. From our flux calibrated spectra, we determined the brightest component of the system to be the hot component whose continuum and emission lines together are responsible for 83% of the V-band light. The circumbinary nebula and its emission lines contribute over 14%, while the red giant is responsible for less than 3%. Our spectra revealed a rich harvest of low ionization emission lines. By measuring how flux in these emission lines varied through the eclipse, we have provided information which should prove useful for future modelling of this symbiotic system.

Using elemental abundances for 1.26 million K giants in the LAMOST DR8 value-added catalog, we analyze the chemical abundances of the Virgo Overdensity (VOD) and Hercules-Aquila Cloud (HAC). We find two distinct chemical populations in both overdensities, which is in disagreement with the mainstream hypothesis that both overdensities are composed of materials from a single merger event, namely Gaia-Sausage-Enceladus (GSE). The two populations show different chemical trends: one exhibits low metallicities and high $\alpha$ abundances, and the other shows high metallicities and low $\alpha$ abundances, which is associated with the recently discovered Nereus and Virgo Radial Merger (VRM) components in the local stellar halo, respectively. The Nereus component in these overdensities uniquely exhibits a decreasing trend in the [Fe/H]-[Mn/Fe] plane. Out of all observed Milky Way dwarf galaxies, this trend is only found in the Sculptor dwarf galaxy, which provides clues for the properties of Nereus progenitor. We also find that the velocity ellipse with high aniostropy parameters that is usually considered to be part of GSE are actually a mix of the two components. Both overdensities are well-mixed in kinematic spaces, confirming recent claims that the debris of merger pairs are kinematically indistinguishable in a recent simulation. We find that the velocity ellipses of the VRM stars in these overdensities have large inclination angles, which may be an indication of the merger time in simulated merger events.

We applied Stars' Galactic Origin (StarGO) algorithm in 7D space (i.e., [Fe/H], [Mg/Fe], [Al/Fe], $L_{z}$, $J_{r}$, $J_{z}$, $E$) to analyze stars in the inner halo with APOGEE DR17. We identified some known substructures, including Gaia-Sausage-Enceladus (GSE), Sagittarius Stream (Sgr), LMS-1 (Wukong), Thamnos, Metal-Weak Thick Disk (MWTD) and Aleph. Additionally, we identified an undefined metal-poor group (UDG, with [Fe/H] $<-0.8$ dex) probably linked to known substructure like Aleph, as well as a high $\alpha$-abundance substructure (HAS) associated with both the Nyx and Nyx-2 streams. Chemical abundance of the HAS supports the argument that Nyx and Nyx-2 share a common origin. We discovered three substructures, which we refer to as new-substructure candidate-1,2,3 (NSTC-1,2,3). Despite exhibiting disk-like dynamics, these NSTCs demonstrate notably low [Mg/Fe] ($<$ 0.2 dex) and [Al/Fe] ($<$ $-0.15$ dex), similar to the properties of dwarf galaxies. Their high orbital energy and low [$\alpha$/Fe] indicate the association with recent accretion events.

A unified explanation of the variety of long-duration gamma-ray burst (GRB) light curves (LCs) is essential for identifying the dissipation mechanism and possibly the nature of their central engines. In the past, a model was proposed to describe GRB LCs as the outcome of a stochastic pulse avalanche process, possibly originating from a turbulent regime, and it was tested by comparing average temporal properties of simulated and real LCs. Recently, we revived this model and optimised its parameters using a genetic algorithm (GA), a machine-learning-based approach. Our findings suggested that GRB inner engines may operate near a critical regime. Here we present an advanced version of the model, which allows us to constrain the peak flux distribution of individual pulses, and evaluate its performance on a new dataset of GRBs observed by the Fermi Gamma-ray Burst Monitor (GBM). After introducing new model parameters and a further comparison metric, that is the observed signal-to-noise (S/N) distribution, we test the new model on three complementary datasets: CGRO/BATSE, Swift/BAT, and Fermi/GBM. As in our previous work, the model parameters are optimised using a GA. The updated sets of parameters achieve a further reduction in loss compared to both the original model and our earlier optimisation. The different values of the parameters across the datasets are shown to originate from the different energy passbands, effective areas, trigger algorithms, and, ultimately, different GRB populations of the three experiments. Our results further underpin the stochastic and avalanche character of the dissipation process behind long GRB prompt emission, with an emphasis on the near-critical behaviour, and establish this new model as a reliable tool for generating realistic GRB LCs as they would be seen with future experiments.

Type-2 quasars (QSO2s) are active galactic nuclei (AGN) seen through a significant amount of dust and gas that obscures the central supermassive black hole and the broad line region. Here we present new mid-infrared spectra of the central kiloparsec of five optically-selected QSO2s at redshift z~0.1 obtained with JWST/MIRI/MRS. These QSO2s belong to the QSOFEED sample and they have log Lbol=45.5-46.0 erg/s, global SFRs that place them above the main sequence, and practically identical optical spectral shape and [OIII] luminosity, but their nuclear mid-infrared spectra exhibit an unexpected diversity of both continua and features. They show: 1) 9.7 micron silicate features going from emission (strength of S9.7=0.5) to relatively strong absorption (S9.7=-1.0) and 18 and 23 micron silicates either in emission or flat. In addition, two of the QSO2s show absorption bands of CO, H2O, and aliphatic grains, indicating different levels of nuclear obscuration across the sample. 2) [NeV]/[NeII] ratios ranging from 0.1 to 2.1 and [NeIII]/[NeII] from 1.0 to 3.5, indicating different coronal line and ionizing continuum strengths. 3) Warm molecular gas masses of 1-4x10^7 Msun and warm-to-cold gas mass ratios of 1-2%, with molecular gas excitation likely due to jet-induced shocks in J1430+1339, and to UV heating and/or turbulence in J1509+0434. 4) PAH emission features with equivalent widths ranging from <0.002 to 0.075 micron, from which we measure a larger contribution from neutral molecules (PAH 11.3/6.2=1.3-3.4) and SFRs<3-7 Msun/yr. This unprecedented dataset allowed us to start exploring the role of various AGN and galaxy properties including ionizing continuum, obscuration, electron density, and jet-ISM interactions on some of the spectral differences listed above, but larger samples are now required to fully understand the diversity of QSO2s' nuclear mid-infrared spectra.

The Mini-SiTian (MST) project is a pathfinder for China's next-generation large-scale time-domain survey, SiTian, aimed at discovering variable stars, transients, and explosive events. MST generates hundreds of thousands of transient alerts every night, approximately 99\% of which are false alarms, posing a significant challenge to its scientific goals. To mitigate the impact of false positives, we propose a deep learning-based solution and systematically evaluate thirteen convolutional neural networks. The results show that ResNet achieves exceptional specificity (99.70\%), EfficientNet achieves the highest recall rate (98.68\%), and DenseNet provides balanced performance with a recall rate of 94.55\% and specificity of 98.66\%. Leveraging these complementary strengths, we developed a bagging-based ensemble classifier that integrates ResNet18, DenseNet121, and EfficientNet\_B0 using a soft voting strategy. This classifier achieved the best AUC value (0.9961) among all models, with a recall rate of 95.37\% and specificity of 99.25\%. It has now been successfully deployed in the MST real-time data processing pipeline. Validation using 5,000 practically processed samples with a classification threshold of 0.798 showed that the classifier achieved 88.31\% accuracy, 91.89\% recall rate, and 99.82\% specificity, confirming its effectiveness and robustness under real application conditions.

This paper provides a comprehensive introduction to the Mini-SiTian Real-Time Image Processing pipeline (STRIP) and evaluates its operational performance. The STRIP pipeline is specifically designed for real-time alert triggering and light curve generation for transient sources. By applying the STRIP pipeline to both simulated and real observational data of the Mini-SiTian survey, it successfully identified various types of variable sources, including stellar flares, supernovae, variable stars, and asteroids, while meeting requirements of reduction speed within 5 minutes. For the real observational dataset, the pipeline detected 1 flare event, 127 variable stars, and 14 asteroids from three monitored sky regions. Additionally, two datasets were generated: one, a real-bogus training dataset comprising 218,818 training samples, and the other, a variable star light curve dataset with 421 instances. These datasets will be used to train machine learning algorithms, which are planned for future integration into STRIP.

This paper outlines the scientific goals and observational strategies of the Mini-SiTian array. Mounted at Xinglong Observatory, the Mini-SiTian array consists of three 30 cm telescopes and has been in operation since 2022. The large field of view, combined with the capability for multi-band photometric observations, enables the Mini-SiTian array to perform rapid follow-up observations to identify optical counterparts of gravitational waves, capture the early light curves of tidal disruption events and supernovae, and monitor stellar flares, Be star outbursts, and cataclysmic variable stars, although its limiting magnitude is not very deep. By collaborating with the Xinglong 2.16-m telescope and leveraging a real-time image processing pipeline, simultaneous photometric and spectroscopic observations could be performed to reveal their underlying physical mechanisms. The observational and research experience provide critical guidance for the implementation of the full-scale SiTian project in the future.

The SiTian project, with its vast field of view, will become an ideal platform for asteroid scientific research. In this study, we develop a pipeline to analyze the photometry of asteroids and derive their periods from the data collected by the SiTian pathfinder project Mini-SiTian (MST). The pipeline is applied to the MST f02 region, a MST test region with a sky area of $2.29^{\circ} \times 1.53^{\circ}$. Rotation periods of 22 asteroids are derived by the obtained light curves analysis. Among them, there are 8 asteroids available in the Asteroid Lightcurve Photometry Database (ALCDEF), and 6 of them with more photometric points ($>$200) have similar period parameters as the ones in ALCDEF. Additionally, the periods for 14 of these asteroids are newly obtained and are not listed in ALCDEF. This study demonstrates the feasibility of asteroid photometric research by the SiTian project. It shows that future observations from the SiTian project will provide even more photometry of asteroids, significantly increasing the number of available light curves. The potential vast photometric data of asteroids will help us to further understand the physics of asteroids, their material composition, and the formation and evolution of the solar system.

The SiTian project, designed to utilize 60 telescopes distributed across multiple sites in China, is a next-generation time-domain survey initiative. As a pathfinder for the SiTian project, the Mini-SiTian (MST) has been proposed and implemented to test the SiTian's brain and data pipeline, and to evaluate the feasibility of its technology and science cases. Mounted at the Xinglong Observatory, the MST project comprises three 30 cm telescopes and has been operated since Nov. 2022. Each telescope of the MST possesses a large field of view, covering $2.29^{\circ}$ $\times$ $1.53^{\circ}$ FOV, and is equipped with $g'$, $r'$ and $i'$ filters, respectively. Acting as the pioneer of the forthcoming SiTian project, the MST is dedicated to the discovery of variable stars, transients, and outburst events, and has already obtained some interesting scientific results. In this paper, we will summarize the first-two-year operation of the MST project.

We simulate the optical searching of gravitational-wave electromagnetic counterpart of the binary neutron star (BNS) merger event (i.e., a kilonova) using the ground based {\it SiTian} project prototype telescope with a 5-min limiting magnitude of 22.0, triggered by LIGO and Virgo gravitational wave detectors during the ongoing O4 run. Our simulations show that an average of 0.17-0.25 kilonova events can be observed over the entire O4 period of $\sim 2$ years in the most optimistic case we set, while no kilonova can be detected in other cases. We note that it is beneficial for {\it SiTian}'s kilonova searching by extending the exposure time to gain deeper limiting magnitude despite the rapid decline of kilonova luminosity.

The Mini-SiTian Array serves as a pathfinder for the SiTian project, which aims to survey the entire sky in $gri$ bands every 30 minutes, reaching a limiting magnitude of 21. This special issue features 11 papers covering the design, operation, data reduction, and early scientific results from two years of Mini-SiTian observations. The insights gained from these pathfinder experiments represent a significant milestone toward the full realization of the SiTian project.

As a pathfinder of the SiTian project, the Mini-SiTian (MST) array, employed three commercial CMOS cameras, represents a next-generation, cost-effective optical time-domain survey project. This paper focuses primarily on the precise data processing pipeline designed for wide-field, CMOS-based devices, including the removal of instrumental effects, astrometry, photometry, and flux calibration. When applying this pipeline to approximately 3000 observations taken in the Field 02 (f02) region by MST, the results demonstrate a remarkable astrometric precision of approximately 70--80\,mas (about 0.1\,pixel), an impressive calibration accuracy of approximately 1\,mmag in the MST zero points, and a photometric accuracy of about 4\,mmag for bright stars. Our studies demonstrate that MST CMOS can achieve photometric accuracy comparable to that of CCDs, highlighting the feasibility of large-scale CMOS-based optical time-domain surveys and their potential applications for cost optimization in future large-scale time-domain surveys, like the SiTian project.

The Mini-SiTian project, which is the pathfinder for the SiTian project, utilizes three 30 cm telescopes equipped with commercial CMOS cameras (ZWO ASI6200MM Pro) to simulate large-area time-domain survey. Due to the avoidance of the traditional mechanical shutter, the CMOS camera is favorable in time-domain survey projects. In the future, the SiTian telescope array will employ a two-by-two scientific-grade mosaic CMOS camera to survey a 10,000-degree square area every 30 minutes. Therefore, the performance of CMOS directly determines the detection capability of SiTian telescopes for transient sources, and a comprehensive understanding of the performance of CMOS cameras is crucial. In this research, laboratory testing was conducted to thoroughly evaluate three cameras by assessing several critical parameters, including bias stability, dark current, pixel anomalies, linearity, gain, and read noise. We find exceptional short-term bias stability with standard deviations below 0.02 ADU, negligible dark current of approximately 0.002 e$^{-}$ pixel$^{-1}$ s$^{-1}$ at $0^\circ\text{C}$, and excellent linearity with nonlinearity consistently below $\pm$ 0.5\%, and a small proportion (0.06\% to 0.08\%) of pixels with anomalous responses. Furthermore, our analysis demonstrates uniform gain values across all cameras, ranging from 0.252 to 0.255 e$^{-}$ ADU$^{-1}$, with low readout noise, measured to be below 1.6 e$^{-}$ using conventional methods. We also propose a novel method for pixel-level gain and read noise calculation for CMOS sensors, which revealed a narrow gain distribution and a low median read noise of 1.028 e$^-$ for one of the cameras. The laboratory testing of the ZWO ASI6200MM Pro cameras indicates their potential to meet the requirements of time-domain surveys for the Mini-SiTian project.

Time-domain astronomy is one of the most important areas. Large sky area, deep-field, and short timescale are the priority of time-domain observations. SiTian is an ambitious ground-based project processing all sky optical monitoring, aiming for sky-survey timescale of less than 1 day. It is developed by the Chinese Academy of Sciences, an integrated network of dozens of 1-m-class telescopes deployed worldwide. The Mini-SiTian Telescope Array is carried out for demonstrations on optical design, group scheduling, and software pipeline developments, to overcome the high technical and financial difficulties of SiTian project. One array contains three 300 mm F/3 telescope, with FOV of 5 degrees over 400-1000 nm wavelength range. The Mini-SiTian Telescope Array is now under commissioning in Xinglong Observatory, and a perfect platform for technical research and educational purposes.

W Ursa Majoris (W UMa)-type contact binary systems (CBs) with period--luminosity (PL) relations are valuable distance indicators. The PL relations of CBs are affected by metallicity. Here, we establish PL relations and period--luminosity--metallicity (PLZ) relations in nine bands from optical to mid-infrared ($BP$, $G$, $RP$, $J$, $H$, $K_S$, $W1$, $W2$, $W3$) and in five Wesenheit bands based on Gaia DR3 parallaxes. The dispersion of PLZ relations gradually decreases from the optical to mid-infrared bands, with the minimum dispersion of 0.138 mag. We fit the best PL relations for three bands ($W1$, $W_{G,BP,RP}$, $W_{W1,BP,RP}$) under different parallax uncertainty criteria and determine a parallax (after correction) zero point of $zp_\varpi=24\pm4$ $\mu$as. After fixing the parallax zero point, we find that the total zero errors of the PL and PLZ relation are smallest when the parallax uncertainty is less than 2\%, resulting in a calibrated CB luminosity with an accuracy of 0.33\%, which is more accurate than the classical Cepheids. Furthermore, by examining the absolute magnitude residuals in different metallicity intervals, we find that the use of a linear metallicity effect is appropriate for CBs with different metallicities. These results indicate that CBs are excellent standard candles.

Writing proposals and job applications is arguably one of the most important tasks in the career of a scientist. The proposed ideas must be scientifically compelling, but how a proposal is planned, written, and presented can make an enormous difference. This Perspective is the third in a series aimed at training the writing skills of professional astronomers. In the first two papers we concentrated on the writing of papers, here we concentrate on how proposals and job applications can be optimally written and presented. We discuss how to select where to propose or apply, how to optimise your writing, and add notes on the potential use of artificial intelligence tools. This guide is aimed primarily at more junior researchers, but we hope that our observations and suggestions may also be helpful for more experienced applicants, as well as for reviewers and funding agencies.

Sonification, or conveying data using non-verbal audio, is a relatively niche but growing approach for presenting data across multiple specialist domains including astronomy, climate science, and beyond. The STRAUSS Python package aims to provide such a tool, which builds upon previous approaches to provide a powerful means to explore different ways of expressing data, with fine control over the output audio and its format. STRAUSS is a free, open source (FOSS) package, designed to allow flexible and effective sonification to be integrated into data workflows, in analogy to widely used visualisation packages. The remit of STRAUSS is broad; it is intended to be able to bridge between ad-hoc solutions for sonifying very particular datasets, and highly technical compositional and sound-design tools that are not optimised for sonification, or may have a steep learning curve. The code offers a range of approaches to sonification for a variety of contexts (e.g. science education, science communication, technical data analysis, etc). To this end, STRAUSS is packaged with a number of examples of different sonification approaches, and preset configurations to support "low-barrier, high-ceiling" approach. STRAUSS has been used to produce both educational resources and analysis tools.

We study the dynamics of Friedmann-Lemaître-Robertson-Walker models where a dark energy component with a quadratic equation of state (EoS) nonlinearly interacts with cold dark matter. Thus, two energy scales naturally come into play: $\rho_*$ is the scale at which the nonlinearity of the EoS becomes relevant; $\rho_i$ is the energy scale around which the interaction starts to play an important role in the dynamics. Our focus is to understand whether there are parameter ranges for this system that can produce non-singular bouncing and emergent cosmologies for any initial condition. We complete a dynamical systems analysis, and find the parameter range such that trajectories always expand from a high energy non-singular de Sitter state. For flat and negative curvature models this de Sitter state is represented by a fixed point, the asymptotic past from which the universe emerges from. We find a subset of positive curvature models that during contraction get arbitrarily close to the de Sitter state, thus having a quasi-de Sitter bounce, then emerge from the bounce and expand, evolving toward spatial flatness. We find that the dimensionless parameter $q\equiv \rho_*/\rho_i$, which measures the relative strength of the nonlinear terms in the system, plays a crucial role in the topology of the phase space. When $q < 3$, some trajectories expand toward a singularity, while others evolve toward a low energy cosmological constant at late-times, with a subset going through a decelerated matter dominate era before the final acceleration. When $q > 3$, all trajectories are non-singular, and evolve toward a late-time cosmological constant. We find a subclass in this case in which all trajectories have at least one decelerated matter dominated phase, and accelerate at late-times.

The standard cosmological model, $\Lambda$CDM, assumes isotropy on large cosmic scales. However, recent studies using galaxy cluster scaling relations reported an apparent $H_0$ anisotropy at $5.4\sigma$ that could be attributed to large bulk flows extending beyond $500\,\mathrm{Mpc}$, in disagreement with $\Lambda$CDM. To quantify the statistical tension of the observational galaxy cluster data used in past studies with $\Lambda$CDM, we utilise the isotropic ($2.8\,\mathrm{Gpc})^3$ run of the FLAMINGO ($\Lambda$CDM) simulations, the largest hydrodynamical cosmological simulation available to date. We create 1728 simulated lightcones and study the apparent level of anisotropy traced by X-ray and thermal Sunyaev-Zeldovich scaling relations in the same cluster sample selection and methodology as in Migkas et al. (2021, arXiv:2103.13904). We find the probability of such apparent anisotropies randomly emerging in cluster scaling relations within a $\Lambda$CDM universe to be $0.12\%\, (3.2\sigma)$. The discrepancy goes up to $\sim 3.6\sigma$ when modelled as a bulk flow at $z < 0.1$. We find that statistical noise accounts for over $80\%$ of the anisotropy amplitude in each lightcone, with large peculiar velocities contributing less than $20\%$. We also show that anisotropy amplitudes are highly sensitive to the intrinsic scatter in the scaling relations, with tighter relations providing stronger constraints. Nevertheless, the tension between Migkas et al. (2021, arXiv:2103.13904) and $\Lambda$CDM persists, however, at a lower significance than previously reported.

Cas OB5 is an OB association located at a distance of 2.5-3 kpc that intercepts the Perseus spiral arm. It carries a moderate amount of reddening ($A_V \sim$ 2-3 mag) and contains several well-known open clusters within its boundaries, such as King 12, NGC 7788, and NGC 7790. The availability of modern clustering algorithms, together with \textit{Gaia} DR3 kinematics and complementary spectroscopic data, makes it a suitable site for studies of Galactic structure. We seek to quantify the spatial scale of star formation in the spiral arms, using Cas OB5 as a pilot target before extending our study to more distant and extinguished regions of the Galaxy. We selected 129,695 candidate OBA stars in a 6x8 deg$^2$ region around Cas OB5. We applied a spectral energy distribution (SED) fitting process to this sample to derive the physical parameters. Through this process, we found 56 379 OBA stars, which we then clustered using HDBSCAN. We identified 17 open clusters inside this area, four of which appear to form a coherent structure that we identify as Cas OB5. Nevertheless, our findings suggest that these clusters belong to two different age groups despite sharing a similar position and kinematics. Spectroscopic observations confirm the youth of NGC 7788 (10-15 Myr) compared to NGC 7790 ($110\pm15\:$Myr). We have determined a spatial scale for star formation of a few tens of pc to a few hundreds of pc, comparing the clustered to the diffuse population of Cas~OB5 across this part of the Perseus arm. A spectroscopic analysis was required to complement the clustering algorithm, so that we could separate younger OCs (tracers of the spiral arm) from older ones. These results highlight the need to combine these techniques to fully disentangle the Milky Way structure.

The distance duality relation (DDR) relates two independent ways of measuring cosmological distances, namely the angular diameter distance and the luminosity distance. These can be measured with baryon acoustic oscillations (BAO) and Type Ia supernovae (SNe Ia), respectively. Here, we use recent DESI DR1, Pantheon+, SH0ES and DES-SN5YR data to test this fundamental relation. We employ a parametrised approach and also use model-independent Generic Algorithms (GA), which are a machine learning method where functions evolve loosely based on biological evolution. When we use DESI and Pantheon+ data without Cepheid calibration or big bang nucleosynthesis (BBN), there is a $2\sigma$ violation of the DDR in the parametrised approach. Then, we add high-redshift BBN data and the low-redshift SH0ES Cepheid calibration. This reflects the Hubble tension since both data sets are in tension in the standard cosmological model $\Lambda$CDM. In this case, we find a significant violation of the DDR in the parametrised case at $6\sigma$. Replacing the Pantheon+ SNe Ia data by DES-SN5YR, we find similar results. For the model-independent approach, we find no deviation in the uncalibrated case and a small deviation with BBN and Cepheids which remains at 1$\sigma$. This shows the importance of considering model-independent approaches for the DDR.

The transiting long-period exoplanets are the most interesting follow-up targets for the next-generation instruments, using the most sophisticated observational techniques. However, the scarcity in their transit events often leads to larger uncertainties in their known ephemeris, also resulting in larger biases in their known physical properties in many cases. In this work, I have used the newer publicly available observations from TESS and CHEOPS for five very interesting long-period transiting exoplanets, i.e., HD95338 b, TOI-2134 c, K2-290 c, TOI-1898 b, and TOI-813 b, combined with the previously reported observations, to reanalyze their transit properties and estimate the updated ephemeris. The analyses also incorporated a critical noise treatment algorithm, which uses well-tested techniques such as wavelet denoising and Gaussian process regression, to effectively reduce the impact of various noise components present in the lightcurves on the estimated parameters. The study has resulted in a more precise estimation of ephemeris for all the targets, with the precision in the estimated periods being better than 5 seconds, except for TOI-813 b, for which the precision in the estimated period is better than 21 seconds. The other transit parameters also got updated, with statistically significant improvements seen in most of the cases, the major notable improvements being in the estimated value of the impact parameter of TOI-1898 b, the orbital semi-major axis of TOI-2134 c, and the radius of HD95338 b. Although long-period exoplanets are expected to show more significant transit timing variations in the presence of other undetected planetary mass objects, no such variations were recorded for these targets.

We use HDGAS hydrodynamic simulations to study the impact of active galactic nucleus (AGN) feedback on the conversion of atomic-gas to molecular-gas within the circumnuclear-disc (CND) of a typical AGN-dominated galaxy. The comparison of CI, CII, and CO line intensities and their ratios in the HDGAS post-processing radiative-transfer analysis reveals the complex interplay between AGN-activity, cold molecular gas properties, and the physical processes governing the evolution of star-formation in galaxies. Our results demonstrate that the CI/CO intensity ratio serves as a reliable indicator of the atomic-to-molecular gas transition. We present the probability distribution function (PDF) and abundance trends of various metal species related to molecular H$2$ gas, highlighting differences in clumpiness and intensity maps between AGN feedback and NoAGN models. The profile of the integrated intensity (moment-0) maps shows that the AGN-feedback model exhibits a lower CI/CO intensity ratio in the vicinity of the supermassive black hole (< 50 pc), indicating a smaller atomic-gas abundance and the presence of positive AGN-feedback. Our simulations have successfully predicted the presence of faint-CO emissions extending to larger radii from the galactic center. We also explore the relationships between CII/CO and CI/CII intensity ratios, as well as the ratios versus CO intensity, which provides insights into the "CO-dark" issues. One notable feature in the later time-scale of the AGN model is the presence of a "CO-dark" region, where the intensity of CO emission ($\rm I_{CO}$) is depleted relative to the H$_2$ column density ($N_{\rm H_2}$) compared to the NoAGN model.

The vertical metallicity gradient of the Galactic disk offers valuable insights into the disk's formation and chemical evolution over time. We utilized the LAMOST-LRS young stellar sample to investigate this gradient and found that it approaches zero as stellar effective temperature (or age) increases (or decreases) across various Galactocentric distances. To validate this result, we analyzed 295 open clusters younger than 3 Gyr and 976 classical cepheids within the Galactic disk. The findings confirmed that, within a given narrow age range, the vertical metallicity gradient is effectively zero. This relationship between metallicity and age supports the ``upside-down'' disk formation theory, as it indicates that the youngest and most metal-rich stars dominate the midplane, while older and more metal-poor stars formed at larger vertical heights and currently tend to be at these heights. Overall, our results align well with theoretical predictions, offering further insight into the chemical evolution and structural properties of the Milky Way.

Galactic diffuse gamma-ray flux measured by the Tibet AS{\gamma} experiment and the total Galactic gamma-ray flux measured by Large High Altitude Air Shower Observatory (LHAASO) are found to be consistent, within the statistical and systematic uncertainties, for the inner Galactic Plane region in the sub-PeV energy range (E > 10^14 eV). The result suggests that the sub-PeV Galactic gamma-ray flux is dominated by the diffuse emission. On the other hand, the LHAASO observations suggest that the sub-PeV gamma-ray sources presented in the first LHAASO catalog possibly give a significant contribution to the total sub-PeV Galactic gamma-ray emission ($\approx$ 60%). However, the estimate must be regarded as a conservative upper limit. In fact, current gamma-ray observations imply that many of the sub-PeV gamma-ray sources detected by LHAASO have a cutoff or significant softening in their energy spectra in the several tens of TeV energy range, and the resolved-source contribution to the total sub-PeV Galactic gamma-ray emission should be much lower than the above estimate. More sophisticated discussion about the origin of the sub-PeV Galactic gamma-ray emission requires detailed spectral studies of the individual gamma-ray sources and an accurate estimate of the contamination of the source fluxes from the diffuse emission.

The nature and formation history of our Galaxy's largest and most enigmatic stellar cluster, known as Omega Centauri (ocen) remains debated. Here, we offer a novel approach to disentangling the complex stellar populations within ocen based on phylogenetics methodologies from evolutionary biology. These include the Gaussian Mixture Model and Neighbor-Joining clustering algorithms applied to a set of chemical abundances of ocen stellar members. Instead of using the classical approach in astronomy of grouping them into separate populations, we focused on how the stars are related to each other. In this way, we could identify stars that likely formed in globular clusters versus those originating from prolonged in-situ star formation and how these stars interconnect. Our analysis supports the hypothesis that ocen might be a nuclear star cluster of a galaxy accreted by the Milky Way with a mass of about 10^9M_sun. Furthermore, we revealed the existence of a previously unidentified in-situ stellar population with a distinct chemical pattern unlike any known population found in the Milky Way to date. Our analysis of ocen is an example of the success of cross-disciplinary research and shows the vast potential of applying evolutionary biology tools to astronomical datasets, opening new avenues for understanding the chemical evolution of complex stellar systems.

External photoevaporation is one of the dominant mechanisms for mass loss from protoplanetary discs. However this mass loss is theoretically expected to depend upon the microphysical properties of protoplanetary discs, which are currently poorly constrained in observations. In this work we explore the impact of microphysics on the bulk evolution of discs. The polycyclic aromatic hydrocarbon (PAH) abundance, and the extent to which grain growth has occurred in the disc have profound effects on the strength of mass loss rates due to external photoevaporation, which in turn can have a significant impact on the disc evolution, impacting disc radii and accretion rates over time. The strongest sensitivity is to whether grain growth has occurred in the disc, which reduces the amount of dust entrained in the wind to shield the disc, thus increasing the rate at which gas is lost. Additionally, larger PAH abundances result in stronger heating and higher mass loss rates, but to a lesser extent than grain growth. We find that plausible variations in the PAH abundance and disc dust evolution can leave observable differences in disc populations. This work highlights the importance of obtaining observational constraints of the microphysical properties of protoplanetary discs. Future observations from JWST should soon be able to provide these constraints.

Acoustic oscillations in stars are sensitive to stellar interiors. Frequency differences between overtone modes -- large separations -- probe stellar density, while differences between low-degree modes -- small separations -- probe the sound speed gradient in the energy-generating core of main sequence Sun-like stars, and hence their ages. At later phases of stellar evolution, characterised by inert cores, small separations are believed to lose much of their power to probe deep interiors and simply become proportional to large separations. Here, we present clear evidence of a rapidly evolving convective zone as stars evolve from the subgiant phase into red giants. By measuring acoustic oscillations in 27 stars from the open cluster M67, we observe deviations of proportionality between small and large separations, which are caused by the influence of the bottom of the convective envelope. These deviations become apparent as the convective envelope penetrates deep into the star during subgiant and red giant evolution, eventually entering an ultra-deep regime that leads to the red giant branch luminosity bump. The tight sequence of cluster stars, free of large spreads in ages and fundamental properties, is essential for revealing the connection between the observed small separations and the chemical discontinuities occurring at the bottom of the convective envelope. We use this sequence to show that combining large and small separations can improve estimations of the masses and ages of field stars well after the main sequence.

Lorentz invariance violation (LIV) has long been recognized as an observable low-energy signature of quantum gravity. In spite of a great effort to detect LIV effects, so far only lower bounds have been derived. The high energy photons from the gamma ray burst GRB 221009A have been detected by the LHAASO collaboration and one at ${\cal E} \simeq 251 \, \rm TeV$ by the Carpet collaboration using a partial data set. Very recently, the Carpet collaboration has completed the full data analysis, reporting further support for their previously detected photon now at ${\cal E} = 300^{+ 43}_{- 38} \, {\rm TeV}$, which manifestly clashes with conventional physics. Taking this result at face value, we derive the first evidence for LIV and we show that such a detection cannot be explained by axion-like particles (ALPs), which allow for the observation of the highest energy photons detected by LHAASO. We also outline a scenario in which ALPs and LIV naturally coexist. If confirmed by future observations our finding would represent the first positive result in quantum gravity phenomenology.

DESTINY+ is an upcoming JAXA Epsilon medium-class mission to flyby multiple asteroids including Phaethon. As an asteroid flyby observation instrument, a telescope mechanically capable of single-axis rotation, named TCAP, is mounted on the spacecraft to track and observe the target asteroids during flyby. As in past flyby missions utilizing rotating telescopes, TCAP is also used as a navigation camera for autonomous optical navigation during the closest-approach phase. To mitigate the degradation of the navigation accuracy, past missions performed calibration of the navigation camera's alignment before starting optical navigation. However, such calibration requires significant operational time to complete and imposes constraints on the operation sequence. From the above background, the DESTINY+ team has studied the possibility of reducing operational costs by allowing TCAP alignment errors to remain. This paper describes an autonomous optical navigation algorithm robust to the misalignment of rotating telescopes, proposed in this context. In the proposed method, the misalignment of the telescope is estimated simultaneously with the spacecraft's orbit relative to the flyby target. To deal with the nonlinearity between the misalignment and the observation value, the proposed method utilizes the unscented Kalman filter, instead of the extended Kalman filter widely used in past studies. The proposed method was evaluated with numerical simulations on a PC and with hardware-in-the-loop simulation, taking the Phaethon flyby in the DESTINY+ mission as an example. The validation results suggest that the proposed method can mitigate the misalignment-induced degradation of the optical navigation accuracy with reasonable computational costs suited for onboard computers.

We present a characterization of the dust environment of long-period comet C/2023 A3 (Tsuchinshan-ATLAS) by analyzing an extensive dataset, including dust tail images and photometric measurements, from 10 au pre-perihelion to 1.6 au post-perihelion, using our forward Monte Carlo code. For this analysis, we combine pre-perihelion images from the Zwicky Transient Facility with post-perihelion images from the Sierra Nevada Observatory (IAA-CSIC, Granada, Spain), along with both amateur and professional photometric measurements, including Af$\rho$ and magnitude data. We find that the dust loss rate increases monotonically from the assumed start of activity at a heliocentric distance of approximately 15 au down to 4 au inbound, where the comet exhibits a notable decrease in dust production by about one order of magnitude. Following this period of reduced activity, the dust production rate rises again toward perihelion, reaching a peak production rate of 10$^5$ kg s$^{-1}$. The size distribution of the particles follows a power law, with its index decreasing toward perihelion, along with a reduction in the minimum particle radius, leading to both brightness and dust mass being dominated by small particles at perihelion. The particle speeds exhibit a dependence on heliocentric distance ($r_h$) that closely follows the classical $r_h^{-0.5}$ law near perihelion but deviates at larger heliocentric distances. We demonstrate that the dependence of particle speeds on the cosine of the solar zenith angle at the emission point plays a significant role in shaping the synthetic dust tails and in the formation of a dark stripe along the tail axis observed in high spatial resolution images near the comet's perihelion.

We present a comprehensive analysis of the MIRI Extremely Red Object Virgil, a Lyman-$\alpha$ emitter at $z_{spec} = 6.6379 \pm 0.0035$ with the photometric properties of a Little Red Dot. Leveraging new JWST/MIRI imaging from the MIDIS and PAHSPECS programs, we confirm Virgil's extraordinary nature among galaxies in JADES/GOODS-South, exhibiting a strikingly red NIRCam-to-MIRI color (F444W $-$ F1500W = $2.84\pm0.04$~mag). Deep NIRSpec/PRISM spectroscopy from the OASIS program offers key insights into the host galaxy, revealing properties of an average star-forming galaxy during Cosmic Reionization, such as a subsolar metallicity, low-to-moderate dust content, and a relatively high ionization parameter and electron temperature. By estimating the star formation rate of Virgil from UV and H$\alpha$, we find evidence that the galaxy is either entering or fading out of a bursty episode. Although line-ratio diagnostics employed at high-$z$ would classify Virgil as an Active Galactic Nucleus (AGN), this classification becomes ambiguous once redshift evolution is considered. Nonetheless, Virgil occupies the same parameter space as recently confirmed AGNs at similar redshifts. The new deep MIRI data at 15~$\mu$m reinforce the AGN nature of Virgil, as inferred from multiple spectral energy distribution (SED) fitting codes. Virgil's rising infrared SED and UV excess resemble those of Dust-Obscured Galaxies (DOGs) studied with Spitzer at Cosmic Noon, particularly blue-excess HotDOGs. Our results highlight the need for a multi-wavelength approach incorporating MIRI to uncover such extreme sources at $z\gtrsim6$ and to shed light on the interplay between galaxy evolution and early black hole growth during Cosmic Reionization.

Broadband photometric reverberation mapping (RM) provides a measure of the size of the continuum-emitting region in active galactic nuclei (AGN). Previous monitoring campaigns of PG 2130+099 disagree as to whether the continuum emitting region size is consistent with that predicted for a standard optically thick geometrically thin accretion disk. We present $\sim$6 months of observations from several robotic telescopes, providing the highest cadence and widest wavelength coverage photometric RM study of PG 2130+099 to date. Our results indicate that inferred size of the continuum-emitting region in PG 2130+099, like many recently observed AGN, is larger than the simplest predictions for an irradiated geometrically thin, optically thick accretion disk. We also perform a flux-flux analysis, finding a variable spectrum broadly consistent with a disk, and a constant component with enhanced $\textit{i}$-band emission, potentially due to H$\alpha$. We find some evidence of increasing lag with luminosity, but previous lag measurements are too uncertain to be definitive.

Context: Using Gaia DR3 data, Jiménez-Arranz et al. (2024a) computed the LMC bar pattern speed using three different methods. One of them suggested that the LMC might be hosting a bar that barely rotates, and is slightly counter-rotating with respect to the LMC disc, with a pattern speed of $\Omega_p = -1.0 \pm 0.5$ km/s/kpc. Aims: To confirm that tidal interactions can trigger the LMC hosting a non-rotating bar due to its interaction with the SMC, which could cause the LMC bar to slow down significantly until it (momentarily) stops. Methods: We analyse a subset of models (K9 and K21) from the KRATOS suite (Jiménez-Arranz et al. 2024b) where we detected non-rotating bars. We make use of two different methods to track the evolution of the bar pattern speed: the program this http URL (Dehnen et al. 2023), and temporal finite-differences of the change rate the bar major axis' phase angle. Results: In the second LMC-SMC-like pericenter passage of K9, the bar of the LMC-like galaxy weakens to almost disappear and regenerates with a pattern speed that suffers a slowdown from $\Omega_p \sim 20$ km/s/kpc to $\Omega_p \sim 0$ km/s/kpc in less than 75 Myr. Then, the bar rotates at less than $\Omega_p \sim 3-5$ km/s/kpc for around 100 Myr, until it recovers the initial (before interaction) pattern speed of $\Omega_p \sim 10$ km/s/kpc. The results for the K21 simulation are comparable. Conclusions: This work is the first direct evidence that galactic bars can be slowed down or even stopped by tidal interaction, which strengthens the possibility of the LMC hosting a non-rotating bar, and can add an alternative formation scenario for observed slow-rotating bars.

We present a concurrent fitting of spectra and light curves of the whole population of detected gamma-ray pulsars. Using a synchro-curvature model we compare our theoretical output with the observational data published in the Third Fermi Pulsar Catalog, which has significantly increased the number of known gamma-ray pulsars. Our model properly fits all the spectra and reproduces well a considerable fraction of light curves. Light curve fitting is carried out with two different techniques, whose strong points and caveats are discussed. We use a weighted reduced \{chi}^2 of light curves in time domain, and the Euclidean distance of the Fourier transform of the light curves, i.e. transforming the light curves to the frequency domain. The performance of both methods is found to be qualitatively similar, but individual best-fit solutions may differ. We also show that, in our model based on few effective parameters, the light curve fitting is basically insensitive to the timing and spectral parameters of the pulsar. Finally, we look for correlations between model and physical parameters, and recover trends found in previous studies but without any significant correlation involving geometrical parameters.

Star clusters are crucial for understanding how stars evolve. Their colour-magnitude diagrams show the effects of stellar evolution of approximately coeval objects with the same chemical composition. Furthermore, the determination of their astrophysical parameters (age, distance, colour excess and metallicity) together with their spatial distribution provides information about the structure and the evolution of the Galaxy itself. Using data from the \textit{Gaia} DR3 and 2MASS catalogues, we develop methodologies for characterizing open clusters. Precise membership lists, mean astrometric parameters and radii are obtained. Using photometric data from both data sources, we carried out new age calibrations that rely on morphological indices based on colour ($\Delta BR$) and magnitude ($\Delta G$) differences between the red clump and the turnoff for a sample of 34 open clusters with ages covering the interval $8.3 < \log[t({\rm yr})] < 9.9$. A set of age calibration functions based on \textit{Gaia} morphological age indices are determined for the first time. We demonstrate their accuracy, obtaining a mean residual of 0.06 dex in $\log[t(yr)]$. Our results also show that stellar evolution models tend to predict the difference $\Delta G$. However, they typically overestimate the difference $\Delta BR$ for objects younger than $\log[t({\rm yr})] = 8.8$.

Wave or fuzzy dark matter produced with high momenta behaves in many ways like hot particle dark matter while also possessing seemingly different phenomenology due to wave interference. We develop wave perturbation theory to show that white noise density fluctuations generated by the interference of high-momenta waves are gravitationally unstable in the usual way during matter domination above the free streaming scale and stabilize below the free streaming scale, much like the analogous effects for massive neutrinos in hot dark matter. We verify and illustrate these effects in the density power spectra of Newtonian Schrödinger-Poisson simulations. In the cosmological context, this would cause a gradual suppression of the initial white noise isocurvature perturbations below the free streaming scale at matter radiation equality, unlike cold dark matter isocurvature fluctuations, and virial stability of dark matter halos.

Halo bias links the statistical properties of the spatial distribution of dark matter halos to those of the underlying dark matter field, providing insights into clustering properties in both general relativity (GR) and modified-gravity scenarios such as $f(R)$ models. While the primary halo mass-dependent bias has been studied in detailed, the secondary bias, which accounts for the additional dependencies on other internal halo properties, can offer a sensitive probe for testing gravity beyond the $\Lambda$CDM model. To quantify any potential deviations between $\Lambda$CDM and $f(R)$ gravity models in halo clustering, at both the primary and secondary level, as well as in the distributions of halo properties in the cosmic web. Using $N$-body simulations of $f(R)$ gravity models, we assess the scaling relations and the primary and secondary bias signals of halo populations on the basis of a halo-by-halo estimator of large-scale effective bias. Our analysis is performed using halo number density as the independent variable. The relative difference in the effective bias between the $f(R)$ models and $\Lambda$CDM is sensitive, albeit slightly, to the power index of modified gravity. The largest deviations from GR are measured for low-mass halos, where the average bias at fixed number density decreases by up to 5\% for fixed scaling indices. We also show that the scaling relations for some environmental properties, including neighbour statistics, Mach number and local overdensity, exhibit small but non-negligible deviations (~3-5\%) from GR for a wide range of number densities. Our results also suggests that the properties of halos in sheets and voids show the largest departures from GR (> 10\% in some cases). In terms of secondary bias, we do not find any statistically significant deviations with respect to $\Lambda$CDM for any of the properties explored in this work.

James Webb Space Telescope (JWST) has opened up a new chapter in infrared astronomy. Besides the discovery and a deeper understanding of various astrophysical sources, JWST can also uncover the non-gravitational nature of dark matter (DM). If DM is QCD axion or an eV-scale Axion-like particle (ALP), it can decay into two photons in the infrared band. This will produce a distinct line signature in the spectroscopic observations made by JWST. Using the latest NIRSpec IFU spectroscopic observations from JWST, we put the strongest bound on the photon coupling for QCD axion/ ALP DM in the mass range between 0.47 and 2.55 eV. In particular, we are able to probe a new mass range for ALP DM between $\sim$ 0.47 eV to 0.78 eV beyond what can be probed by globular cluster observations. We constrain well-motivated and UV complete models of QCD axion and ALP DM, including predictions from some models derived from string theory and/ or various Grand Unification scenarios. Future JWST observations of DM-rich systems with a better understanding of the astrophysical and instrumental backgrounds can thus enable us to potentially discover QCD axion and ALP DM. The datasets used in this work are available at: this https URL

We investigate general relativistic effects on the photon spectrum emitted from decaying (or annihilating) particle dark matter in the halo surrounding a primordial black hole. The spectrum undergoes significant modification due to gravitational redshifts, which induces broadening as a result of the intense gravitational field near the black hole. This characteristic alteration in the photon spectrum presents a unique observational signature. Future observations of such spectral features may provide critical evidence for a mixed dark matter scenario, involving both primordial black holes and particle dark matter.

Steven Weinberg productive scientific life teaches us many things, one of the most important of which is the power of his example. This essay contains personal reminiscences and a speculation about gravitational wave propagation based on one of his very last papers -- as seems fitting, given his penchant for putting interesting physics into the essays he wrote for other luminaries over the years. (See arXiv:2502.10979 [this http URL-ph] for a less personal summary of his main scientific accomplishments.)

This article investigates the noise performance of a two-element phased array and interferometer containing a recently introduced self-interference canceler, which in the context of this work acts as a mutual-coupling canceler. To this end, a general framework is proposed to permit noise analysis of this network and a large variety of other networks. The framework-based numerical analysis for a two-element-phased array shows that the addition of the canceler significantly increases the beam-equivalent noise temperature. For a two-element interferometer used in cosmology, this increase in noise temperature is still acceptable as the sky noise temperature in the 20-to-200 MHz band is high. When used in an interferometer, the canceler provides the ability to null mutual coherence at the interferometer output. The ability to provide matching to reduce the sensitivity of the null in mutual coherence to the phase of the 90deg hybrids in the canceler is discussed.

A perturbative strategy for inflation described by two-inflaton fields is developed using a mathematical analogy with the renormalization-group. Two small quantities, $\alpha$ and $\lambda$, corresponding to standard slow-roll parameters are defined and systematic expansions of all inflationary quantities in terms of powers of $\alpha$ and $\lambda$ are found. No other slow-roll parameters are needed. To illustrate this perturbative method in the multi-field context, we adopt a simple two-inflaton model with quadratic potentials in the parameter range where both fields contribute similarly to the dynamics of inflation. The model, even though it is not a viable alternative for phenomenological description of the inflationary period, nicely illustrates subtleties of the perturbative approach. In particular, this method allows us to derive two independent gauge-invariant scalar perturbations that are conserved in the superhorizon limit, overcoming typical problems that emerge in multi-field inflation. Furthermore, it is possible to perform nonperturbative resummations that allow to study the model in a true multi-field regime. We derive tensor and scalar power spectra to the next-to-next-to-leading and next-to-leading orders, respectively, as well as their spectral indices. The hierarchy between the two scalar perturbations allows us to single out the dominant entry of the scalar power-spectrum matrix. Modifications due to the second inflaton occur already at the leading order. Finally, we explain why the quadratic two-inflaton model is not compatible with the present experimental constraints even though non-trivial corrections to scalar perturbations do emerge.

AMS-02 is a detector currently in operation onboard the International Space Station (ISS). One of the main scientific goals of the spectrometer is the measurement of charged particle fluxes. The detector design makes possible the identification of particles and antiparticles by precise measurement of particle momentum in the AMS-02 Silicon Tracker, and velocity in the Cherenkov (RICH) detector. The RICH is able to measure the isotopic composition of the light elements (up to charge Z=5) in the kinetic energy range from a few GeV/n to about 10 GeV/n. However, the velocity reconstruction for charge 1 particles is particularly challenging due to the low number of photons they produce in the RICH detector which can lead to wrong event reconstruction. In this paper, we show the high potential of the Multilayer Perceptron deep learning model (MLP-BFGS) for identification of signal and the background due to interactions inside the AMS-02 detector, and to significantly improve particle identification by its mass.

If a coupling between the inflaton and the Gauss-Bonnet term is introduced, many models of inflation that were ruled out by the most recent Planck data can be made viable again. The predictions for the scalar spectral index and tensor-to-scalar ratio are typically computed using the slow-roll approximation. In this paper we instead study the full equations of motion and determine the necessary initial conditions for reasonable inflation epoch. We derive the conditions under which the Friedmann equation admits positive solutions for the Hubble parameter. Then we study the possibility of the inflaton becoming trapped in a local potential minimum induced by the Gauss-Bonnet term. Finally we demonstrate the results on monomial potential models with a quadratic and a quartic potential and show that the slow-roll approximation becomes imprecise in the quartic case.

Within the framework of general $U(1)$ scenario, we demonstrate that the ultra high energy neutrinos recently detected by KM3NeT could originate from a decaying right handed neutrino dark matter (DM), with a mass of 440 PeV. Considering DM production via freeze-in, we delineate the parameter space that satisfies the observed relic abundance and also lies within the reach of multiple gravitational wave detectors. Our study provides a testable new physics scenario, enabled by multi-messenger astronomy.

In this paper, we perform a comprehensive analysis of the quasinormal modes in the Joshi-Malafarina-Narayan (JMN-1) naked singularity by investigating its response to linear perturbations, including scalar, electromagnetic, and gravitational perturbations. To analyze the stability of the JMN-1 naked singularity under axial perturbations, we compute the quasinormal mode frequencies using the Wentzel-Kramers-Brillouin method. The quasinormal mode frequencies provides information about the stability of spacetime, with the real part of the frequency determining the oscillation rate and the imaginary part governing the decay or growth of perturbations. Our results indicate that by imposing appropriate boundary conditions, we find that the JMN-1 spacetime remains dynamically stable under axial perturbations.

The Brans-Dicke scalar-tensor cosmological models are studied in both Einstein and Jordan frames, using hydrodynamical and self-interacting scalar field representations of the energy-momentum tensor, leading to the same background solutions. The main features of the corresponding cosmological scenarios are determined. In many cases, we identify a transition from decelerated to accelerated regimes. The properties of the self-interacting scalar field, including the corresponding potential, are determined. Finally, some possible realistic configurations are analyzed.

We present an analysis of gravitational wave polarization modes within Gravitational Quantum Field Theory (GQFT), a unified theoretical framework reconciling general relativity and quantum field theory. Our study focuses on five fundamental polarization states predicted in GQFT: two tensor ($+, \times$), two vector (x, y), and one scalar (breathing) mode, focusing on their distinctive detection signatures in space-based interferometers like LISA and Taiji. Using first-order orbital dynamics in the Solar System Barycenter frame, we identify three novel observational features: (1) characteristic interference patterns between polarization modes, (2) distinctive null-point signatures enabling mode discrimination, and (3) sky-position-dependent optimal detection windows. Our approach provides complete sky coverage through polarization mapping while remaining fully compatible with existing mission designs, notably avoiding the need for challenging direct breathing-mode measurements. The results are presented through comprehensive sky maps, offering both theoretical insights into gravitational wave polarization and practical tools for future detector networks. This work establishes a new paradigm for testing fundamental gravity theories through their unique polarization fingerprints, with particular relevance for upcoming multi-messenger gravitational wave astronomy.

Density Functional Theory (DFT) is a robust framework for modeling interacting many-body systems, including the equation of state (EoS) of dense matter. Many models, however, rely on energy functionals based on assumptions that have not been rigorously validated. We critically analyze a commonly used ansatz for confinement, where the energy functional scales with density as $U \propto n^{\frac{2}{3}}$ . Our findings, derived from a systematic non-local energy functional, reveal that this scaling does not capture the dynamics of confinement. Instead, the energy functional evolves from $n^2$ at low densities to $n$ at high densities, governed by an infrared cutoff. These results suggest that models relying on such assumptions should be revisited to ensure more reliable EoS construction.

We discuss a minimal renormalizable Pati-Salam theory based on the $SU(4)_{\rm C}\,\times\,SU(2)_{\rm L}\,\times\,SU(2)_{\rm R}$ gauge group, with unification scale Higgs multiplets taken as $SU(2)_{\rm L}$ and $SU(2)_{\rm R}$ doublets, which lead to neutrino Dirac picture. Although a number of scalar particles could be light, even lying at the LHC energies, the unification scale is hopelessly out of reach in any foreseeable future. Moreover, phase transition in the early Universe leads to the production of magnetic monopoles and domain walls, both incompatible with the standard cosmological model. A small explicit breaking of the discrete left-right symmetry allows the domain walls to decay, and in the process possibly sweep away the monopoles, analogously to the previously discussed case of $SU(5)$ grand unified theory. This leaves an important imprint of gravitational waves, within the reach of next generation searches, correlated with monopole detection and new light particles at collider energies. The theory has a dark matter candidate in the form of an inert scalar doublet, with a mass below TeV, which can further trigger electroweak baryogenesis.

Correspondences between apparently distant concepts are ubiquitous in theoretical physics. In the context of Black Holes (BHs), Quasi-Normal Modes (QNMs) were shown to be linked to both the shadows, in the so-called eikonal limit, and Gray-Body Factors (GBFs), using the WKB approximation. We test the accuracy of the quasinormal modes-graybody factors correspondence's in the context of the Hawking spectra for static and rotating black hole configurations, with particular attention to the superradiant regime. Our analysis reveals the correspondence failure to accurately reproduce the Hawking spectrum due to divergences. Furthermore, we bridge the gap between black shadows and graybody factors by drawing a correspondence between such quantities in the case of generic static and spherically symmetric spacetime configurations. The shadow-GBF correspondence is tested for some case studies, including regular BHs, and its limitations and applicability are thereof discussed. This study opens new perspectives, by introducing a new correspondence and remarking on the caution needed when considering these connections.