Papers for Wednesday, Feb 19 2025

The connection between $\gamma$-ray flares and blazars is a topic of active research, with few sources exhibiting distinct enough such outbursts to be able to conclusively connect them to features in their jet morphology. Here we present an investigation of the sole $\gamma$-ray flare of the blazar OJ 248 thus far, in association with its jet structure, as revealed by very long baseline interferometry (VLBI). We find that throughout the course of the $\gamma$-ray flare, the fractional linear polarisation increases in the jet of OJ 248, and the VLBI electric vector position angles (EVPAs) turn perpendicular to the bulk jet flow. We interpret this behaviour as a moving shock, travelling through a recollimation shock and upscattering photons via the inverse Compton scattering process, producing a $\gamma$-ray flare; we discuss possible mechanisms. Our hypothesised shock-shock interaction scenario is a viable mechanism to induce such EVPA rotations in both optical and radio bands.

Atmospheric characterisation of temperate sub-Neptunes is the new frontier of exoplanetary science with recent JWST observations of possible Hycean world K2-18b. Accurate modelling of atmospheric processes is essential to interpreting high-precision spectroscopic data given the wide range of possible conditions in the sub-Neptune regime, including on potentially habitable planets. Notably, convection is an important process which can operate in different modes across sub-Neptune conditions. Convection can act very differently in atmospheres with a high condensible mass fraction (non-dilute atmospheres) or with a lighter background gas, e.g. water convection in a H$_2$-rich atmosphere, and can be much weaker or even shut down entirely in the latter case. We present a new mass-flux scheme which can capture these variations and simulate convection over a wide range of parameter space for use in 3D general circulation models (GCMs). We validate our scheme for two representative cases, a terrestrial-like atmosphere and a mini-Neptune atmosphere. In the terrestrial case, considering TRAPPIST-1e with an Earth-like atmosphere, the model performs near-identically to Earth-tuned models in an Earth-like convection case. In the mini-Neptune case, considering the bulk properties of K2-18b and assuming a deep H$_2$-rich atmosphere, we demonstrate the capability of the scheme to reproduce non-condensing convection. We find convection occurring at pressures greater than 0.3 bar and the dynamical structure shows high-latitude prograde jets. Our convection scheme will aid in the 3D climate modelling of a wide range of exoplanet atmospheres, and enable further exploration of temperate sub-Neptune atmospheres.

The blue supergiant (BSG) domain contains a large variety of stars whose past and future evolutionary paths are still highly uncertain. Since binary interaction plays a crucial role in the fate of massive stars, investigating the multiplicity among BSGs helps shed light on the fate of such objects. We aim to estimate the binary fraction of a large sample of BSGs in the Small Magellanic Cloud within the Binarity at LOw Metallicity (BLOeM) survey. In total, we selected 262 targets with spectral types B0-B3 and luminosity classes I-II. This work is based on spectroscopic data collected by the GIRAFFE instrument, mounted on the Very Large Telescope, which gathered nine epochs over three months. Our spectroscopic analysis for each target includes the individual and peak-to-peak radial velocity measurements, an investigation of the line profile variability, and a periodogram analysis to search for possible short- and long-period binaries. By applying a 20 km s$^{-1}$ threshold on the peak-to-peak radial velocities above which we would consider the star to be binary, the resulting observed spectroscopic binary fraction for our BSG sample is 23 $\pm$ 3$\%$. In addition, we derived reliable orbital periods for 41 spectroscopic binaries and potential binary candidates, among which there are 17 eclipsing binaries, including 20 SB1 and SB2 systems with periods of less than 10 days. We reported a significant drop in the binary fraction of BSGs with spectral types later than B2 and effective temperatures less than 18 kK, which could indicate the end of the main sequence phase in this temperature regime. We found no metallicity dependence in the binary fraction of BSGs, compared to existing spectroscopic surveys of the Galaxy and Large Magellanic Cloud.

This work consists of the first detailed photometric study of Terzan 6, one of the least known globular clusters in the Galactic bulge. Through the analysis of high angular resolution and multi-wavelength data obtained from adaptive optics corrected and space observations, we built deep, optical and near-infrared color-magnitude diagrams reaching $\approx 4$ magnitudes below the main-sequence turnoff. Taking advantage of 4 different epochs of observations, we measured precise relative proper motions for a large sample of stars, from which cluster members have been solidly distinguished from Galactic field interlopers. A non-canonical reddening law (with $R_V=2.85$) and high-resolution differential reddening map, with color excess variations up to $\delta E(B-V) \approx 0.8 $ mag, have been derived in the direction of the system. According to these findings, new values of the extinction and distance modulus have been obtained: respectively, $E(B-V)=2.36\pm0.05$ and $(m-M)_0=14.46 \pm 0.10$ (corresponding to $d=7.8 \pm 0.3$ kpc). We also provide the first determinations of the cluster center and projected density profile from resolved star counts. The center is offset by more than $7$ arcsec to the east from the literature value, and the structural parameters obtained from the King model fitting to the density profile indicate that Terzan 6 is in an advanced stage of its dynamical evolution. We also determined the absolute age of the system, finding $t=13\pm 1 $ Gyr, in agreement with the old ages found for the globular clusters in the Galactic bulge. From the re-determination of the absolute magnitude of the red giant branch bump and the recent estimate of the cluster global metallicity, we find that Terzan 6 nicely matches the tight relation between these two parameters drawn by the Galactic globular cluster population.

We present identifications and kinematic analysis of 7,426 massive ($\mathrm{\geq}8M_{\odot}$) stars in the Small Magellanic Cloud (SMC), using Gaia DR3 data. We used Gaia ($G_\mathrm{BP}-G_\mathrm{RP}$, $G$) color-magnitude diagram to select the population of massive stars, and parallax to omit foreground objects. The spatial distribution of the 7,426 massive star candidates is generally consistent with the spatial distribution of the interstellar medium, such as H$\alpha$ and H i emission. The identified massive stars show inhomogeneous distributions over the galaxy, showing several superstructures formed by massive stars with several hundred parsecs scale. The stellar superstructures defined by the surface density have opposite mean proper motions in the east and west, moving away from each other. Similarly, the mean line-of-sight velocities of the superstructures are larger to the southeast and smaller to the northwest. The different east-west properties of the superstructures' proper motion, line-of-sight velocity indicate that the SMC is being stretched by tidal forces and/or ram pressure from the Large Magellanic Cloud to the southeast, thereby rejecting the presence of galaxy rotation in the SMC.

Recent years have seen a surge of interest in the community studying the effect of ultraviolet radiation environment, predominantly set by OB stars, on protoplanetary disc evolution and planet formation. This is important because a significant fraction of planetary systems, potentially including our own, formed in close proximity to OB stars. This is a rapidly developing field, with a broad range of observations across many regions recently obtained or recently scheduled. In this paper, stimulated by a series of workshops on the topic, we take stock of the current and upcoming observations. We discuss how the community can build on this recent success with future observations to make progress in answering the big questions of the field, with the broad goal of disentangling how external photoevaporation contributes to shaping the observed (exo)planet population. Both existing and future instruments offer numerous opportunities to make progress towards this goal.

In this work, we trace the complete life cycle of individual GMCs in high-resolution Milky Way-mass galaxy simulations to determine how different stellar feedback mechanisms and galactic-scale processes govern cloud lifetimes, mass evolution, and local star formation efficiency (SFE). We identify GMCs in simulated galaxies and track their evolution using cloud evolution trees. Via cloud evolution trees, we quantify the lifetimes and SFE of GMCs. We further apply our diagnostics on a suite of simulations with varying star formation and stellar feedback subgrid models and explore their impact together with galactic environments to the GMC life cycles. Our analysis reveals that GMCs undergo dynamic evolution, characterized by continuous gas accretion, gravitational collapse, and star formation, followed by disruption due to stellar feedback. The accretion process sustains the gas content throughout most of the GMC life cycles, resulting in a positive correlation between GMC lifetimes and their maximum masses. The GMC lifetimes range from a few to several tens of Myr, with two distinct dynamical modes: (1) GMCs near the galactic center experience strong tidal disturbances, prolonging their lifetimes when they remain marginally unbound; (2) those in the outer regions are less affected by tides, remain gravitationally bound, and evolve more rapidly. In all model variations, we observe that GMC-scale SFE correlates with the baryonic surface density of GMCs, consistent with previous studies of isolated GMCs. Additionally, we emphasize the critical role of galactic shear in regulating GMC-scale star formation and refine the correlation between local SFE and surface density by including its effects. These findings demonstrate how stellar feedback and galactic-scale dynamics jointly shape GMC-scale star formation in realistic galactic environments.

Ultra-hot Jupiters, an extreme class of planets not found in our solar system, provide a unique window into atmospheric processes. The extreme temperature contrasts between their day- and night-sides pose a fundamental climate puzzle: how is energy distributed? To address this, we must observe the 3D structure of these atmospheres, particularly their vertical circulation patterns, which can serve as a testbed for advanced Global Circulation Models (GCM) [e.g. 1]. Here, we show a dramatic shift in atmospheric circulation in an ultra-hot Jupiter: a unilateral flow from the hot star-facing side to the cooler space-facing side of the planet sits below an equatorial super-rotational jet stream. By resolving the vertical structure of atmospheric dynamics, we move beyond integrated global snapshots of the atmosphere, enabling more accurate identification of flow patterns and allowing for a more nuanced comparison to models. Global circulation models based on first principles struggle to replicate the observed circulation pattern [3], underscoring a critical gap between theoretical understanding of atmospheric flows and observational evidence. This work serves as a testbed to develop more comprehensive models applicable beyond our Solar System as we prepare for the next generation of giant telescopes.

Transit spectroscopy usually relies on the integration of one or several transits to achieve the S/N necessary to resolve spectral features. Consequently, high-S/N observations of exoplanet atmospheres are essential for disentangling the complex chemistry and dynamics beyond global trends. In this study, we combined two partial 4-UT transits of the ultrahot Jupiter WASP-121 b, observed with the ESPRESSO at the VLT in order to revisit its titanium chemistry. Through cross-correlation analysis, we achieved detections of H I, Li I, Na I, K I, Mg I, Ca I, Ti I, V I, Cr I, Mn I, Fe I, Fe II, Co I, Ni I, Ba II, Sr I, and Sr II. Additionally, narrow-band spectroscopy allowed us to resolve strong single lines, resulting in significant detections of H$\alpha$, H$\beta$, H$\gamma$, Li I, Na I, K I, Mg I, Ca II, Sr I, Sr II, and Mn I. Our most notable finding is the high-significance detection of Ti I ($\sim$ 5$\sigma$ per spectrum, and $\sim$ 19$\sigma$ stacked in the planetary rest frame). Comparison with atmospheric models reveals that Ti I is indeed depleted compared to V I. We also resolve the planetary velocity traces of both Ti I and V I, with Ti I exhibiting a significant blueshift toward the end of the transit. This suggests that Ti I primarily originates from low-latitude regions within the super-rotating jet observed in WASP-121 b. Our observations suggest limited mixing between the equatorial jet and the mid-latitudes, in contrast with model predictions from GCMs. We also report the non-detection of TiO, which we attribute to inaccuracies in the line list that could hinder its detection, even if present. Thus, the final determination of the presence of TiO must await space-based observations. We conclude that the 4-UT mode of ESPRESSO is an excellent testbed for achieving high S/N on relatively faint targets, paving the way for future observations with the ELT.

We present 2-4 GHz observations of polarized radio galaxies towards eight fast radio bursts (FRBs), producing grids of Faraday rotation measure (RM) sources with sky densities of 9-28 polarized sources per square degree. Using a Bayesian interpolation framework, we constrain Galactic RM fluctuations below ~ 1 degree squared angular scales around the FRB positions. Despite the positions of all eight FRBs far from the Galactic plane, we constrain previously unresolved small-scale Galactic RM structures around six of the eight FRBs. In two of these fields, we find potential changes in the sign of the Galactic RM that are not captured by previous, sparsely sampled RM grid observations. Our Galactic RM estimate towards the FRBs differs between a few rad m^-2 up to ~ 40 rad m^-2 from the all-sky Galactic RM map of Hutschenreuter et al. (2022). Extrapolating our results to the known population of polarized FRB sources, we may be incorrectly interpreting the host galaxy RM for ~ 30% of the FRB source population with current RM grid observations. Measuring small-scale Galactic RM variations is crucial for identifying FRBs in low density and weakly magnetized environments, which in turn could serve as potent probes of cosmic magnetism. This framework of reconstructing continuous Galactic RM structure from RM grid observations can be readily applied to FRBs that fall in the sky coverage of upcoming large-sky radio polarization surveys of radio galaxies, such as the Very Large Array Sky Survey (VLASS) and the Polarization Sky Survey of the Universe's Magnetism (POSSUM).

We perform radiative magnetohydrodynamics simulations in general relativity (GRRMHD) of super-Eddington disk accretion onto neutron stars endowed with a magnetic dipole corresponding to surface strengths not exceeding 100 GigaGauss. Accretion is found to power strong outflows which collimate the emergent radiation of the accretion columns, leading to apparent radiative luminosities of $\sim 100$ Eddington, when the true luminosity is a few Eddington units. Surprisingly, the collimation cone/angle widens with increasing magnetic field. Thus, in our simulations the apparent luminosity of the neutron star is substantially larger for the weaker magnetic fields ($10^{10}\,$G) than for the stronger ones ($10^{11}\,$G). We conclude that a super-Eddington accreting neutron star with the dipole magnetic field $10^{10}\,$G is the most likely source of ultraluminous X-rays.

A recent study shows that if the power spectra (PS) of accreting compact objects consist of a combination of Lorentzian functions that are coherent in different energy bands but incoherent with each other, the same is true for the Real and Imaginary parts of the cross spectrum (CS). Using this idea, we discovered imaginary quasi-periodic oscillations (QPOs) in NICER observations of the black hole candidate MAXI J1820+070. The imaginary QPOs appear as narrow features with a small Real and large Imaginary part in the CS but are not significantly detected in the PS when they overlap in frequency with other variability components. The coherence function drops and the phase lags increase abruptly at the frequency of the imaginary QPO. We show that the multi-Lorentzian model that fits the PS and CS of the source in two energy bands correctly reproduces the lags and the coherence, and that the narrow drop of the coherence is caused by the interaction of the imaginary QPO with other variability components. The imaginary QPO appears only in the decay of the outburst, during the transition from the high-soft to the low-hard state of MAXI J1820+070, and its frequency decreases from approximately 5 Hz to around 1 Hz as the source spectrum hardens. We also analysed the earlier observations of the transition, where no narrow features were seen, and we identified a QPO in the PS that appears to evolve into the imaginary QPO as the source hardens. As for the type-B and C QPOs in this source, the rms spectrum of the imaginary QPO increases with energy. The lags of the imaginary QPO are similar to those of the type-B and C QPOs above 2 keV but differ from the lags of those other QPOs below that energy. While the properties of this imaginary QPO resemble those of type-C QPOs, we cannot rule out that it is a new type of QPO.

Colour light curves of resident space objects (RSOs) encapsulate distinctive features that can offer insights into an object's structure and design, making them an invaluable tool for classification and characterisation. We present the results of the first large systematic colour survey of the GEO belt in which we obtain full-night multi-colour light curves for 112 active geostationary objects between April and May 2023. Colour light curve maps were created to compare and contrast the colours between different satellites and bus configurations. We find that satellites with BSS-702 and STAR-2 buses can be effectively distinguished from the colour measurements on these maps, but comparing the average colour of individual satellites within given solar equatorial phase angle ranges shows that it is difficult to distinguish between bus configurations based on colour alone. We also find tentative evidence to suggest that there is a relationship between colour and time spent on orbit for the Eurostar-3000 class satellites, which is unseen behaviour within other bus configuration classes. The satellites in our sample exhibit `redder' colours than the Sun, which is in agreement with previous findings. We found common light curve features such as symmetrical colour changes as well as unique regions of short timescale glinting which are `bluer' than other regimes within the colour light curves. If these features are indeed seasonal, this would be a powerful characterisation tool. We are able to detect and resolve features in the light curve of the LDPE-3A satellite related to manoeuvres being performed. Finally, we measured the solar panel offsets of 54 satellites in our sample and found variation in the type of colour response. The majority of which did not exhibit any colour change across the solar panel glints compared to them shifting towards 'redder' or 'bluer' colours.

Using spatially resolved spectroscopic data from the MaNGA sample, we investigate the parameters influencing the radial gradients of gas-phase metallicity ($\nabla\log(\mathrm{O/H})$), to determine whether disk formation is primarily driven by coplanar gas inflow or by the independent evolution of distinct regions within the disk. Our results show that $\nabla \log(\mathrm{O/H})$ strongly correlates with local gas-phase metallicity at a given stellar mass, with steeper gradients observed in metal-poorer disks. This trend supports the coplanar gas inflow scenario, wherein the gas is progressively enriched by in situ star formation as it flows inward. In contrast, the radial gradient of stellar mass surface density shows very weak correlations with $\nabla \log(\mathrm{O/H})$, which is inconsistent with the independent evolution mode, where gas inflow, star formation, and metal enrichment occur independently within each annulus of the disk. Furthermore, we find that $\nabla \log(\mathrm{O/H})$ is also closely correlated with an indicator of local gas turbulence $\sigma_{\mathrm{gas}}/R_{\mathrm{e}}$, highlighting the competing roles of turbulence and coplanar inflow in shaping metallicity gradients. Our results provide indirect observational evidence supporting coplanar gas inflow as the driving mechanism for disk evolution.

Accurate knowledge of gamma-ray burst (GRB) redshifts is essential for studying their intrinsic properties and exploring their potential application in cosmology. Currently, only a small fraction of GRBs have independent redshift measurements, primarily due to the need of rapid follow-up optical/IR spectroscopic observations. For this reason, many have utilized phenomenological correlations to derive pseudo-redshifts of GRBs with no redshift measurement. In this work, we explore the feasibility of analytically deriving pseudo-redshifts directly from the Amati and Yonetoku relations. We simulate populations of GRBs that (i) fall perfectly on the phenomenological correlation track, and (ii) include intrinsic scatter matching observations. Our findings indicate that, in the case of the Amati relation , the mathematical formulation is ill-behaved so that it yields two solutions within a reasonable redshift range $z \in [0.1, 10] $. When realistic scatter is included, it may result in no solution, or the redshift error range is excessively large. In the case of the Yonetoku relation, while it can result in a unique solution in most cases, the large systematic errors of the redshift calls for attention, especially when attempting to use pseudo redshifts to study GRB population properties.

Circular polarization in the cosmic microwave background (CMB) offers a promising probe of the parity-violating physics of the early universe. In this paper, we propose a novel method to constrain the primordial circular polarization of high-frequency gravitational waves (GW) in the GHz range. An efficient conversion of gravitons to photons in a transverse cosmological magnetic field at the epoch of last scattering can generate excess chiral photons if the GW background is chiral in nature. This excess radiation distorts the CMB thermal black-body spectrum, which can be estimated by measuring the V-Stokes parameter in the CMB polarization. Using current upper limits on the angular power spectrum of circular polarization $C_l^{VV}$ from the CLASS, MIPOL, and SPIDER experiments, we obtain the most stringent constraints on the characteristic strain and circular polarization of the isotropic background of stochastic GWs at ${40\,\rm GHz}$ and ${150\,\rm GHz}$, respectively. Our work, therefore, provides an interesting possibility to constrain the circular polarization of high-frequency GWs using the V-mode polarization measurements of CMB.

We investigate early galaxy evolution by modeling self-consistently their radially-resolved evolution of gas, stars, heavy elements, and dust content. Our model successfully reproduces various observed properties of JWST-identified galaxies at $z > 5$, including sizes, stellar masses, star formation rates (SFR), metallicities, and dust-to-stellar mass ratios. We show that the star formation efficiency (SFE), $f_\ast \equiv {\rm SFR}/(f_{\rm b} \dot{M}_{\rm h})$, is regulated by the global equilibrium between cosmological gas inflows, star formation, and gas outflows. Our model predicts $f_\ast \lesssim 20~\%$ for galaxies with halo masses of $M_{\rm h} \sim 10^{11-12}\, M_\odot$ down to $z = 5$, allowing them to reach intrinsic UV magnitudes of $M_{\rm UV} \lesssim -22~{\rm mag}$; when dust attenuation is ignored, the predicted UV luminosity function (LF) at $z \sim 12$ agrees well with observations. However, our model also suggests that these galaxies would be heavily obscured by dust, with high optical depths at 1500~Å~of $\tau_{1500} \gtrsim 10$, causing the dust-attenuated UV LF to fall significantly below the observed one. This discrepancy highlights the need for mechanisms that mitigate strong dust attenuation, such as dust evacuation from star-forming regions and/or preferential production of large dust grains. Further exploration of these processes is essential for understanding the early stages of galaxy evolution.

We analytically derive self-similar solutions for a time-dependent, one-dimensional, magnetically driven accretion disk wind model derived from the magnetohydrodynamic equations. The model assumes a geometrically thin, gas-pressure dominated accretion disk, and incorporates both magnetic braking and turbulent viscosity through an extended alpha-viscosity prescription in the vertical and radial directions, respectively. The $\alpha$ parameter for the vertical stress is assumed to vary with the disk aspect ratio. We confirm that our self-similar solutions without the wind matches with the classical solution of Cannizzo et al. (1990) that the mass accretion rate follows the power law of time $t^{-19/16}$, which has been used as a good indicator for the mass accretion rate of a tidal disruption event (TDE) disk. In contrast, in the presence of the wind, the mass accretion and loss rates decay more steeply than $t^{-19/16}$. We also confirm that the power-law indices of the mass accretion and loss rates are consistent with those obtained from the numerical simulations of Tamilan et al. (2024) at late times. In particular, we find that magnetic braking leads to a faster decay of the mass accretion rate, mass loss rate, and bolometric luminosity, and they asymptote to $t^{-5/2}$ in the strong poloidal magnetic field. This steep index can serve as evidence for magnetocentrifugally driven winds with a strong poloidal magnetic field in the context of TDEs.

We report the results of high-resolution 5 GHz Very Long Baseline Array and European VLBI Network observations of 36 nearby galaxies, an extension of the Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies (LeMMINGs) survey. Our sample includes 21 low ionization nuclear emission regions (LINERs), 4 Seyferts, 3 absorption line galaxies (ALGs), and 8 HII galaxies. We achieved an unprecedented detection rate, successfully imaging 23 out of 36 sources with a detection threshold of $\sim$20 $\mu$Jy beam$^{-1}$. The radio sizes are typically of $\leq$ 5 pc. Core identification was achieved in 16 sources, while 7 others were identified as core candidates. Radio luminosities of the sample range from 10$\rm ^{34}$ to 10$\rm ^{38}$ erg s$^{-1}$. Our analysis reveals a predominance of compact core structures, with ten sources exhibiting a one-sided core jet morphology and NGC 2146 exhibiting a rare two-sided jet structure. The study advances our understanding of the compactness of radio sources at various scales, indicating a core-dominated nature in all but one galaxy NGC2655. We find moderate to strong correlations between radio luminosity, black hole mass, optical [O III] line luminosity, and hard X-ray luminosity, suggesting a common active galactic nucleus (AGN) core origin. These results provide new insights into the fundamental plane of black hole activity and support the role of the synchrotron process in Low-luminosity AGN (LLAGN) radio emission.

We present the results of blazar 3C 273 obtained from simultaneous observations obtained using XMM-Newton and NuSTAR satellites during the period 2015-2019 in five epochs. When the spectra are modeled with a power-law, significant residuals arise below 2 keV and in the energy range of 30-78 keV in NuSTAR data. Residuals in the lower energy band represent soft X-ray excess while at higher energies it likely represents Compton reflection hump which might be a weak component arising from dense and cold material. The presence of a faint iron line is present in XMM-Newton observations. We interpret such features as attributed to the coronal emission plus those arising from reflection from an accretion disk. We model the SEDs with the single zone inverse Compton jet model based on Synchrotron Self Compton and External Compton phenomena. It is found that a one-zone synchrotron plus IC model explains quite well the SEDs but the jet component alone fails to fit the multiband X-ray emission for the low state of this object in 2018 and 2019 which arises due to spectral flattening at low energy X-rays, indicating that an additional Seyfert-like thermal component must be present at X-rays. This is further supported by a big blue bump present in the optical/ultraviolet band in all SEDs. Finally, we analyzed all the epochs using relxill model to incorporate relativistic reflection to model those residuals of soft excess and Compton hump in the X-ray bands.

Supermassive binary black holes (SMBBHs) are the anticipated byproducts of galaxy mergers and play a pivotal role in shaping galaxy evolution, gravitational wave emissions, and accretion physics. Despite their theoretical prevalence, direct observational evidence for SMBBHs remains elusive, with only a handful of candidates identified to date. This paper explores optimal strategies and key environments for locating SMBBHs, focusing on observational signatures in the broad Balmer lines. We present a preliminary analysis on a flux-limited sample of sources belonging to an evolved spectral type along the quasar main sequence, and we discuss the spectroscopic clues indicative of binary activity and highlight the critical role of time-domain spectroscopic surveys in uncovering periodic variability linked to binary systems.

In this study, we use a flexible parametrization of the equation of state of dark energy to explore its possible evolution with datasets from the Dark Energy Spectroscopic Instrument (DESI), Planck cosmic microwave background, and either the 5-year Dark Energy Survey (DES) or the Pantheon+ (PP) supernova (SN) compilation. This parametrization, called transitional dark energy, allows for rapid changes in the equation of state but also changes like that in the Chevallier-Polarski-Linder parametrization. We find a 3.8{\sigma} preference for evolving dark energy over {\Lambda}CDM with the DES dataset and a weaker 2.4{\sigma} preference when using the PP dataset. This corroborates the finding of the DESI Collaboration, who found that their baryon acoustic oscillation data preferred evolving dark energy when fit with the CPL parametrization of the equation of state. Our analysis reveals no significant outliers in the DESI data around the TDE best-fit, while the data is asymmetrically distributed around the {\Lambda}CDM best-fit model such that the measured distances are on average smaller. The DESI and SN data both prefer an expansion history that implies a higher dark energy density around z=0.5 than in the Planck-{\Lambda}CDM model, with the inferred equation of state being greater than -1 around z=0 and close to or below -1 at z>0.5. We show that when the expansion rate is greater than that in the Planck-{\Lambda}CDM model (around z=0.5), the growth rate calculated assuming General Relativity is suppressed relative to the Planck-{\Lambda}CDM model, and it rebounds as the expansion rate differences between the models become smaller closer to the present time. The resulting flattening of the $f\sigma_8(z)$ curve compared to the {\Lambda}CDM model could be an independent signature of the temporal evolution of dark energy.

We study the morphology of gamma-ray burst (GRB) afterglows viewed off-axis using a simplified analytical model. We consider steep jets, which are expected to be the most common type. These jets, characterized by steep lateral gradients in energy and Lorentz factor, produce highly beamed emission. The observed signal is dominated by their minimum visible angle at any given time. Consequently, the afterglow morphology depends on when this angle begins to decrease, revealing the inner regions of the jet. Depending on whether this decrease occurs before, at, or after the reverse shock crosses the ejecta, three distinct classes of light curves emerge. In the first scenario, the de-beamed emission can produce a rapidly rising signal even prior to the reverse shock crossing. This is expected in GRBs with long duration, low energy, dense circumburst media, or combinations thereof. In some cases, the ejecta shell can be considered as effectively thick in the inner regions and effectively thin in the outer regions. For forward shocks, the temporal slopes in both regimes are identical, which makes it hard to detect the transition. Reverse shocks, however, have distinct temporal slopes, allowing potential detection of the transition in light curves if their emission surpasses that of the forward shock. The characteristic synchrotron frequency of de-beamed emission evolves independently of jet structure for forward shocks but depends on the lateral energy and Lorentz factor gradients for the reverse shock, with slower evolution for steep energy and shallow Lorentz factor gradients.

We provide a review on semi-empirical models of galaxy formation and evolution. We present a brief census of the three main modeling approaches to galaxy evolution, namely hydrodynamical simulations, semi-analytic models, and semi-empirical models (SEMs). We focus on SEMs in their different flavors, i.e. interpretative, descriptive and hybrid, discussing the peculiarities and highlighting virtues and shortcomings for each of these variants. We dissect a simple and recent hybrid SEM from our team to highlight some technical aspects. We offer some outlook on the prospective developments of SEMs. Finally, we provide a short summary of this review.

Electron irradiation of water-rich ices plays a significant role in initiating the chemical and physical processes on the surface of airless icy bodies in radiation environments such as Europa and Enceladus, as well as on the Moon, comets, and asteroids interacting with the solar wind. The sputtering process by electrons and ions leads to chemical modification and outgassing of their icy surfaces and the subsequent formation of a tenuous atmosphere. Though electron-sputtering yields are known to be lower compared to ion-sputtering yields, one needs to also account for their differential fluxes. In our experiments, the electron-induced sputtering yields of all the gaseous species H2, O, OH, H2O, and O2 are investigated for electron energies lower than 2 keV in terms of partial pressure vs. time of irradiation. The effective averaged change in partial pressure of the desorbed species is converted to the number of sputtered atoms or molecules per second per cm2 from the ice, and then to the sputtering yields (number of species sputtered per electron). Our data agrees well with the previously reported data for the sputtering of O2 and H2O yield for the amorphous ice. We also find that crystalline ice shows significantly lower sputtering yields when compared to amorphous ice, in agreement with the observation of similar trends in the literature. Our work indicates that sputtering yields per keV of O, OH, O2, and H2O drop with increasing electron energy from 0.5 keV to 2 keV.

Jupiter's poles feature striking polygons of cyclones that drift westward over time-a motion governed by beta-drift. This study investigates how beta-drift and the resulting westward motion depend on the depth of these cyclones. Counterintuitively, shallower cyclones drift more slowly, a consequence of stronger vortex stretching. By employing a 2D quasi-geostrophic model of Jupiter's polar regions, we constrain the cyclones' deformation radius, a key parameter that serves as a proxy for their vertical extent, required to replicate the observed westward drift. We then explore possible vertical structures and the static stability of the poles by solving the eigenvalue problem that links the 2D model to a 3D framework, matching the constrained deformation radius. These findings provide a foundation for interpreting upcoming Juno microwave measurements of Jupiter's north pole, offering insights into the static stability and vertical structure of the polar cyclones. Thus, by leveraging long-term motion as a novel constraint on vertical dynamics, this work sets the stage for advancing our understanding of the formation and evolution of Jupiter's enigmatic polar cyclones.

We present seeing-limited (0.8 arcsec) near-infrared integral field spectroscopy data of the type-2 quasars (QSO2s) SDSS J135646.10+102609.0 (J1356) and SDSS J143029.89+133912.1 (J1430, the Teacup), both belonging to the Quasar Feedback (QSOFEED) sample. The nuclear K-band spectra (1.95-2.45 \textmu m) of these radio-quiet QSO2s reveal several $H_2$ emission lines, indicative of the presence of a warm molecular gas reservoir (T$\geq$1000 K). We measure nuclear masses of 5.9, 4.1, and 1.5 $\times 10^3~M_{\odot}$ in the inner 0.8 arcsec diameter region of the Teacup, J1356 north (J1356N), and south nuclei, respectively. The total warm $H_2$ mass budget is $\sim 4.5$ and $\sim 1.3 \times 10^4~M_{\odot}$ for the Teacup and J1356N, implying warm-to-cold molecular gas ratios of $10^{-6}$. The warm molecular gas kinematics, traced with the $H_2$1-0S(1) and S(2) emission lines, is consistent with that of the cold molecular phase, traced by ALMA CO emission at higher angular resolution (0.2 and 0.6 arcsec). In J1430, we detect the blue- and red-shifted sides of a compact warm molecular outflow extending up to 1.9 kpc and with velocities of 450 km/s. In J1356 only the red-shifted side is detected, with a radius of up to 2.0 kpc and velocity of 370 km/s. The outflow masses are 2.6 and 1.5 $\times 10^3~M_{\odot}$ for the Teacup and J1356N, and the warm-to-cold gas ratios in the outflows are 0.8 and 1 $\times 10^{-4}$, implying that the cold molecular phase dominates the mass budget. We measure warm molecular mass outflow rates of 6.2 and 2.9 $\times 10^{-4}~M_{\odot}/yr$ for the Teacup and J1356N, approximately 0.001\% of the total mass outflow rate. We find an enhancement of velocity dispersion in the $H_2$1-0S(1) residual dispersion map of the Teacup, both along and perpendicular to the compact radio jet direction. This enhanced turbulence can be reproduced by simulations of jet-ISM interactions.

In his 1856 Adams Prize essay, James Clark Maxwell demonstrated that Saturn's rings cannot be comprised of a uniform rigid body. This is a consequence of the two-body gravitational interaction between a ring and planet resulting in instability. Similarly, it is also known that a so-called Dyson sphere encompassing a single star would be unstable due to Newton's shell theorem. A surprising finding is reported here that both a ring and a sphere (shell) can be stable in the restricted three-body problem. First, if two primary masses are considered in orbit about their common centre of mass, a large, uniform, infinitesimal ring enclosing the smaller of the masses can in principle be stable under certain conditions. Similarly, a Dyson sphere can, be stable, if the sphere encloses the smaller of the two primary masses, again under certain conditions. These findings extend Maxwell's results on the dynamics of rings and have an interesting bearing on so-called Ringworlds and Dyson spheres from fiction. Moreover, the existence of passively stable orbits for such large-scale structures may have implications for so-called techno-signatures in search for extra-terrestrial intelligence studies.

Many problems in astrophysics cover multiple orders of magnitude in spatial and temporal scales. While simulating systems that experience rapid changes in these conditions, it is essential to adapt the (time-) step size to capture the behavior of the system during those rapid changes and use a less accurate time step at other, less demanding, moments. We encounter three problems with traditional methods. Firstly, making such changes requires expert knowledge of the astrophysics as well as of the details of the numerical implementation. Secondly, some parameters that determine the time-step size are fixed throughout the simulation, which means that they do not adapt to the rapidly changing conditions of the problem. Lastly, we would like the choice of time-step size to balance accuracy and computation effort. We address these challenges with Reinforcement Learning by training it to select the time-step size dynamically. We use the integration of a system of three equal-mass bodies that move due to their mutual gravity as an example of its application. With our method, the selected integration parameter adapts to the specific requirements of the problem, both in terms of computation time and accuracy while eliminating the expert knowledge needed to set up these simulations. Our method produces results competitive to existing methods and improve the results found with the most commonly-used values of time-step parameter. This method can be applied to other integrators without further retraining. We show that this extrapolation works for variable time-step integrators but does not perform to the desired accuracy for fixed time-step integrators.

The return of Comet 1P/Halley will promote a wide interest for ground and space observations of a celestial body of outstanding scientific and cultural interest. In addition to remote observations, space will open the possibility of in situ science similarly to the passage of 1986. In this paper, we first discuss the scientific motivations for a rendezvous mission, capable to overcome the limitations of the flyby missions that took place at that time. In the second part, we describe an example of a rendezvous trajectory that can be carried out with existing power and propulsion technologies. The transfer is made possible by the gravitational assistance of a giant planet. The resulting mission will be capable to reach the comet beyond the distance of Saturn, when the sublimation of super-volatile species will be ongoing, and well before the onset of the sublimation of water (4 AU). After rendezvous, the spacecraft will accompany the comet for several years before, around and after perihelion (July 2061). Our concept mission does not foresee the implementation of solar panels. In this way, operations can occur even inside the dense dust coma at short distance from the nucleus. In the third part of the paper, an innovative imaging system is proposed, with a very large field of view (100°) capable to record on the same frame details on the surface and the surrounding space, in order to follow for several degrees the trajectories of chunks and clouds ejected by pits or fractures, crucial to the understanding of the cometary activity. A concerted effort is needed in the current decade to plan and approve a rendezvous mission to 1P. Indeed, the scenario here described requires launching before 2040, less than 15 years from now. Later launches imply a severe loss of scientific knowledge, because the spacecraft will not be able to reach the comet before the onset of water sublimation.

Detecting stochastic gravitational wave background (SGWB) from cosmic strings is crucial for unveiling the evolutionary laws of the early universe and validating non-standard cosmological models. This study presents the first systematic evaluation of the detection capabilities of next-generation ground-based gravitational wave detector networks for cosmic strings. By constructing a hybrid signal model incorporating multi-source astrophysical foreground noise, including compact binary coalescences (CBCs) and compact binary hyperbolic encounters (CBHEs), we propose an innovative parameter estimation methodology based on multi-component signal separation. Numerical simulations using one-year observational data reveal three key findings: (1) The CE4020ET network, comprising the Einstein Telescope (ET-10 km) and the Cosmic Explorer (CE-40 km and CE-20 km), achieves nearly one order of magnitude improvement in constraining the cosmic string tension $G\mu$ compared to individual detectors, reaching a relative uncertainty $\Delta G\mu / G\mu < 0.5$ for $G\mu > 3.5 \times 10^{-15}$ under standard cosmological framework; (2) The network demonstrates enhanced parameter resolution in non-standard cosmological scenarios, providing a novel approach to probe pre-Big Bang Nucleosynthesis cosmic evolution; (3) Enhanced detector sensitivity amplifies CBHE foreground interference in parameter estimation, while precise modeling of such signals could further refine $G\mu$ constraints by $1-2$ orders of magnitude. This research not only quantifies the detection potential of third-generation detector networks for cosmic string models but also elucidates the intrinsic connection between foreground modeling precision and cosmological parameter estimation accuracy, offering theoretical foundations for optimizing scientific objectives of next-generation gravitational wave observatories.

Since dark matter particles have never been directly detected, we do not know how they move, and in particular we do not know how they fall inside gravitational potential wells. Usually it is assumed that dark matter only interacts gravitationally with itself and with particles of the standard model, and therefore that its motion is governed by Euler's equation. In this paper, we test this assumption for the first time at cosmological scales, by combining measurements of galaxy velocities with measurements of gravitational potential wells, encoded in the Weyl potential. We find that current data are consistent with Euler's equation at redshifts $z\in [0.3,0.8]$, and we place constraints on the strength of a potential fifth force, which would alter the way dark matter particles fall. We find that a positive fifth force cannot exceed 7% of the gravitational interaction strength, while a negative fifth force is limited to 21%. The coming generation of surveys, including the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory and the Dark Energy Spectroscopic Instrument (DESI) will drastically improve the constraints, allowing to constrain a departure from pure gravitational interaction at the level of 2%.

We use data from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) to study the most likely formation of a forced reconnection region and associated plasma blobs, triggered by jet-like structures in a prominence segment. Around 05:44 UT on December 16$^{th}$, 2017, hot jet-like structures lifted from a nearby active region and fell obliquely on one side of the prominence segment with velocities of $\approx$45--65 km s$^{-1}$. These eruptions compressed the boundaries of the prominence and flux rope, forming an elongated reconnection region with inflow velocities of 47--52 km s$^{-1}$ and 36--49 km s$^{-1}$ in the projected plane. A thin, elongated reconnection region was formed, with multiple magnetic plasma blobs propagating bidirectionally at velocities of 91--178 km s$^{-1}$. These dense blobs, associated with ongoing reconnection, may also be linked to the onset of Kelvin-Helmholtz (K-H) instability. The blobs are attributed to plasmoids, moving at slower speeds (91--178 km s$^{-1}$) due to the high density in the prominence segment. The dimensionless reconnection rate varied from 0.57--0.28, 0.53--0.26, and 0.41--0.20, indicating reconnection rate enhancement and supporting the forced reconnection scenario. After reconnection, the prominence plasma heated to 6 MK, releasing significant thermal energy ($\approx$5.4$\times$10$^{27}$ erg), which drained cool prominence plasma and heated it to coronal temperatures. The ubiquity of jets and outflows in the solar atmosphere makes the aforementioned of reconnection and possible co-existence of K-H instability potentially important for the magnetic energy release and heating in the solar atmosphere.

With the 4-meter Multi-Object Spectroscopic Telescope (4MOST) expected to provide an influx of transient spectra when it begins observations in early 2026 we consider the potential for real-time classification of these spectra. We investigate three extant spectroscopic transient classifiers: the Deep Automated Supernova and Host classifier (DASH), Next Generation SuperFit (NGSF) and SuperNova IDentification (SNID), with a focus on comparing the efficiency and purity of the transient samples they produce. We discuss our method for simulating realistic, 4MOST-like, host-galaxy contaminated spectra and determining quality cuts for each classifier used to ensure pure SN Ia samples while maintaining efficient classification in other transient classes. We investigate the classifiers individually and in combinations. We find that a combination of DASH and NGSF can produce a SN Ia sample with a purity of 99.9% while successfully classifying 70% of SNe Ia. However, it struggles to classify non-SN Ia transients. We investigate photometric cuts to transient magnitude and transient flux fraction, finding that both can be used to improve transient classification efficiencies by 7--25% depending on the transient subclass. Finally, we present an example classification plan for live classification and the predicted purities and efficiencies across five transient classes: Ia, Ibc, II, superluminous and non-supernova transients.

Galaxy clusters selected based on overdensities of galaxies in photometric surveys provide the largest cluster samples. Yet modeling the selection function of such samples is complicated by non-cluster members projected along the line of sight (projection effects) and the potential detection of unvirialized objects (contamination). We empirically constrain the magnitude of these effects by cross-matching galaxy clusters selected in the Dark Energy survey data with the \rdmpr$\,$ algorithm with significant detections in three South Pole Telescope surveys (SZ, pol-ECS, pol-500d). For matched clusters, we augment the \rdmpr$\,$catalog by the SPT detection significance. For unmatched objects we use the SPT detection threshold as an upper limit on the SZe signature. Using a Bayesian population model applied to the collected multi-wavelength data, we explore various physically motivated models to describe the relationship between observed richness and halo mass. Our analysis reveals the limitations of a simple lognormal scatter model in describing the data. We rule out significant contamination by unvirialized objects at the high-richness end of the sample. While dedicated simulations offer a well-fitting calibration of projection effects, our findings suggest the presence of redshift-dependent trends that these simulations may not have captured. Our findings highlight that modeling the selection function of optically detected clusters remains a complicated challenge, requiring a combination of simulation and data-driven approaches.

We explore the background equivalence among three dark energy models by constructing explicit mappings between dynamical dark energy (DDE), interacting dark energy (IDE), and running vacuum (RV). In our approach, the dark sector functions that characterize each model-such as the equation of state parameter $\bar{w}(a)$ for DDE, the interaction term $Q$ for IDE, and the functional form $\Lambda(H)$ for RV-are transformed into one another under specific assumptions. Extending previous work by von Marttens et al., we demonstrate that running vacuum models, characterized by $\Lambda(H) = a_0 + a_1 \dot{H} + a_2 H^2$, can be reinterpreted as an interacting dark energy model with $Q = 3H\gamma \hat{\rho}_c$, which in turn is equivalent to a dynamic dark energy model with an appropriately defined $\bar{w}(a)$. Using Bayesian analysis with Type Ia supernovae, Baryon Acoustic Oscillations, and Cosmic Chronometers, our observational constraints confirm that these theoretical equivalences hold at the background level. This study underscores the importance of seeking convergence in dark energy models, facilitating a better understanding of the dark sector.

The paper presents the new enhancement to the Data Center named DC1 at the University of Naples {\em Federico II}. The ICSC funds at INFN have allowed to improve the power and cooling subsystems, while other funds from the PNRR (the STILES project\index{STILES Project}) and funds directly from the MUR have allowed to enhance the computing, storage and network equipments. All these resources are in addition to the IBiSCo cluster and equipments described earlier in this book, but all together, thanks to a strong synergy between projects, have leaded to a very powerful Data Center for scientific applications.

We present the analysis of OGLE-2014-BLG-1760, a planetary system in the galactic bulge. We combine Keck Adaptive Optics follow-up observations in $K$-band with re-reduced light curve data to confirm the source and lens star identifications and stellar types. The re-reduced MOA dataset had an important impact on the light curve model. We find the Einstein ring crossing time of the event to be $\sim$ 2.5 days shorter than previous fits, which increases the planetary mass-ratio and decreases the source angular size by a factor of 0.25. Our OSIRIS images obtained 6 years after the peak of the event show a source-lens separation of 54.20 $\pm$ 0.23 mas, which leads to a relative proper motion of $\mu_{\rm rel}$ = 9.14 $\pm$ 0.05 mas/yr, larger than the previous light curve-only models. Our analysis shows that the event consists of a Jupiter-mass planet of $M_{\rm p}$ = 0.931 $\pm$ 0.117 $M_{\rm Jup}$ orbiting a K-dwarf star of $M_*$ = 0.803 $\pm$ 0.097 $M_{\odot}$ with a $K$-magnitude of $K_{\rm L}$ = 18.30 $\pm$ 0.05, located in the galactic bulge or bar. We also attempt to constrain the source properties using the source angular size $\theta_*$ and $K$-magnitude. Our results favor the scenario of the source being a younger star in the galactic disk, behind the galactic center, but future multicolor observations are needed to constrain the source and thus the lens properties.

The Km3NET collaboration has recently reported the detection of a neutrino event with energy in excess of 100 PeV. This detection is in 2.5-3$\sigma$ tension with the upper limit on the neutrino flux at this energy imposed by IceCube and the Pierre Auger Observatory, if the event is considered part of the diffuse all-sky neutrino flux. We explore an alternative possibility that the event originates from a flare of an isolated source. We show that the data of Km3NET, IceCube and the Pierre Auger Observatory are consistent with the possibility of a source flare of duration $T \lesssim 2$ yr with muon neutrino flux $F \approx 3\times 10^{-10}(1\mbox{ yr }/ T)$ erg cm$^{-2}$ s$^{-1}$. Constraints on the neutrino spectrum indicate that the protons responsible for the neutrino emission have a very hard spectrum in the $E_p\gtrsim 10^{19}$ eV energy range, or otherwise that the neutrinos are produced by photohadronic interactions with infrared photons. The all-sky rate of similar neutrino flaring sources is constrained to be $R\lesssim 0.4/$ yr.

Many spiral galaxies host magnetic fields with energy densities comparable to those of the turbulent and thermal motions of their interstellar gas. However, quantitative comparison between magnetic field properties inferred from observation and those obtained from theoretical modeling has been lacking. In Paper I we developed a simple, axisymmetric galactic dynamo model that uses various observational data as input. Here we apply our model to calculate radial profiles of azimuthally and vertically averaged magnetic field strength and pitch angle, gas velocity dispersion and scale height, turbulent correlation time and length, and the sizes of supernova remnants for the galaxies M31, M33, M51, and NGC 6946, using input data collected from the literature. Scaling factors are introduced to account for a lack of precision in both theory and observation. Despite the simplicity of our model, its outputs agree fairly well with galaxy properties inferred from observation. Additionally, we find that most of the parameter values are similar between galaxies. We extend the model to predict the magnetic field pitch angles arising from a combination of mean-field dynamo action and the winding up of the random small-scale field owing to the large-scale radial shear. We find their magnitudes to be much smaller than those of the pitch angles measured in polarized radio and far infrared emission. This suggests that effects not included in our model, such as effects associated with spiral arms, are needed to explain the pitch angle values.

In this paper, we model the plasma environment of Ganymede by means of a collisionless test particle simulation. By coupling the outputs from a DSMC simulation of Ganymede's exosphere (i.e. number density profiles of neutral species such as $\mathrm{H}$, $\mathrm{H_2}$, $\mathrm{O}$, $\mathrm{HO}$, $\mathrm{H_2O}$, $\mathrm{O_2}$ for which we provide parametrisation) with those of a MagnetoHydroDynamic simulation of the interaction between Ganymede and the Jovian plasma (i.e. electric and magnetic fields), we perform a comparison between simulated ion plasma densities and ion energy spectra with those observed in-situ during 6 close flybys of Ganymede by the Galileo spacecraft. We find that not only our test particle simulation sometimes can well reproduce the in-situ ion number density measurement, but also the dominant ion species during these flybys are $\mathrm{H_2^+}$, $\mathrm{O_2^+}$, and occasionally $\mathrm{H_2O^+}$. Although the observed ion energy spectra cannot be reproduced exactly, the simulated ion energy spectra exhibit similar trends to those observed near the closest approach and near the magnetopause crossings but at lower energies. We show that the neutral exosphere plays an important role in supplying plasma to Ganymede's magnetised environment and that additional mechanisms may be at play to energise/accelerate newborn ions from the neutral exosphere.

Machine learning methods are well established in the classification of quasars (QSOs). However, the advent of light curve observations adds a great amount of complexity to the problem. Our goal is to use the Zwicky Transient Facility (ZTF) to create a catalog of QSOs. We process the ZTF DR20 light curves with a transformer artificial neural network and combine the Pan-STARRS (PS), AllWISE, and Gaia surveys with extreme gradient boosting. Using ZTF g-band data with at least 100 observational epochs per light curve, we obtain 97% F1 score for QSOs. We find that with 3 day median cadence, a survey time span of at least 900 days is required to achieve 90% QSO F1 score. However, one can obtain the same score with a survey time span of 1800 days and the median cadence prolonged to 12 days. We find that ZTF classification is superior to the PS static bands, and on par with WISE and Gaia measurements. Additionally, we find that the light curves provide the most important features for QSO classification in the ZTF dataset. We robustly classify objects fainter than the $5\sigma$ SNR limit at $g=20.8$ by requiring $g < \mathrm{n_{obs}} / 80 + 20.375$. For this sample, we run inference with added WISE observations, and find 4,849,574 objects classified as QSOs. For 33% of QZO objects, with available WISE data, we publish redshifts with estimated error $\Delta z/(1 + z) = 0.14$.

The phenomenology of in-vacuo dispersion, an effect such that quantum properties of spacetime slow down particles proportionally to their energies, has been a very active research area since the advent of the Fermi telescope. One of the assumptions made in this 15-year effort is that the phenomenology of in-vacuo dispersion has a particle-energy sweet spot: the energy of the particle should be large enough to render the analysis immune to source-intrinsic confounding effects but still small enough to facilitate the identification of the source of the particle. We use the gigantic energy of KM3-230213A as an opportunity to challenge this expectation. For a neutrino of a few hundred PeVs a transient source could have been observed at lower energies several years earlier, even assuming that the characteristic scale of in-vacuo dispersion be close to the Planck scale. We report that GRB090401B is in excellent directional agreement with KM3-230213A, and we discuss a strategy of in-vacuo-dispersion analysis suitable for estimating the significance of KM3-230213A as a GRB090401B-neutrino candidate. We find significance at the level of a p-value of 0.01, not small enough to warrant any excitement, but small enough to establish the point that a handful of such coincidences would suffice to meaningfully test in-vacuo dispersion.

The cosmological Mass Varying Neutrino model beyond the standard $\Lambda$CDM scenario is considered. The interaction of the fermionic field and the scalar field with the inverse power law Ratra-Peebles potential via the Yukawa coupling is studied in detail. Depending on the model parameter $\alpha$ of the Ratra-Peebles potential of the scalar field, the expansion rate of the universe, the mass equation, the mass of the scalar field, the sum of neutrino masses, the mutual influence of the sum of neutrino masses and the value of the scalar field Ratra-Peebles potential as well as the total density of the thermodynamic potential of the coupled fermionic and scalar fields at the critical point are explored. The values of the sum of neutrino masses $m_\nu(a_0)\leq0.07~{\rm eV}$ calculated for values of the model parameter $\alpha$ of the Ratra-Peebles potential $0<\alpha\leq0.016$ are consistent with the constraint $m_\nu(a_0)<0.071~{\rm eV}$ of the cosmological DESI measurements, while the values $m_\nu(a_0)\leq0.45~{\rm eV}$ for $0<\alpha\leq0.143$ are consistent with the upper limit $m_\nu(a_0)<0.45~{\rm eV}$ obtained in the KATRIN experiment.

The longitudinal components of massive vector fields generated during inflation constitute a well-motivated dark matter candidate, with interesting phenomenological implications. During the epoch of radiation domination following inflation, their spectrum exhibits a peak at small scales, whose amplitude and position are governed by the parameters of the dark matter model. We calculate the stochastic gravitational wave spectrum induced at second order in fluctuations by such a longitudinal vector peak. We demonstrate that the amplitude of the gravitational wave spectrum can, in principle, reach significant values at nano-Hertz frequencies or lower. This result suggests a novel gravitational wave probe to test inflationary vector dark matter scenarios, independent from assumptions on the coupling of dark vectors to the Standard Model. Additionally, we derive new analytical formulas for the longitudinal vector transfer functions during radiation domination, offering a valuable tool for characterising the convolution integrals that govern the properties of the induced gravitational waves.

Ram pressure stripping is a well-known environmental quenching mechanism that removes gas from galaxies infalling into groups and clusters. In some extreme examples of ram pressure stripping, galaxies with extended gas tails show evidence of enhanced star formation prior to quenching. In this work we use a sample of 5277 local satellite galaxies in which a stripped tail of gas has not necessarily been observed, to quantify the strength of ram pressure-enhanced star formation and compare these results to a control sample of 8360 field galaxies. We use u-band imaging from the Ultraviolet-Near Infrared Northern Survey (UNIONS) as a star formation tracer and several metrics to quantify star formation asymmetry. We compare these results to environmental properties of the galaxy, such as their time since infall and host halo mass, to constrain the degree of ram pressure enhanced star formation as a function of environment. We find no significant differences between the satellite and the field samples. We further restrict our sample to galaxies which we most expect to be experiencing significant ram pressure but find no strong evidence of these galaxies having systematically enhanced star formation. Finally, we investigate the properties of the most asymmetric galaxies in our sample and again find no strong evidence of ram pressure-induced star formation enhancement. We conclude that any star formation enhancement must be small for infalling galaxies, suggesting that this effect is either uncommon or short-lived.

This article introduces a nonstatic Reissner-Nordstr$\phi$m metric, a metric that does not emit electromagnetic waves but can emit gravitational waves. We first use the GR theory to study a charged spherically symmetric gravitational source (CSSGS), the obtained results are further improved in comparison with the previous studies. In particular, this article considers that the field is not necessarily static. The metric tensors $ g_{\mu\nu} $ are considered both outside and inside the gravitational source (the results show that in the first case $ g_{\mu\nu} $ are time independent, in the latter case they are time dependent). The gravitational acceleration and the event horizon of a charged black hole are investigated. The results prove that the gravitational field is always attractive. We then use the perturbative $ f(R) $ theory to consider CSSGS. The obtained results not only correct the solution of Einstein's equation in magnitude (this will describe astronomical and cosmological quantities more accurately than Einstein's equation), but also reveal new effects. Outside the gravitational source, the metric tensors can depend on time, this makes it possible for a spherically symmetric gravitational source to emit gravitational waves (Einstein's equation cannot give this effect). However, a spherically symmetric field still does not emit electromagnetic waves. Next we present a new method for embedding the spherically symmetric metrics of a star (or a black hole) in the background of the FLRW cosmological. Finally, we discuss the uniqueness of the solutions of the f(R) theory. The perturbative TOV equation is also found.

Kaon condensation in hyperon-mixed matter [($Y$+$K$) phase], which may be realized in neutron stars, is discussed on the basis of chiral symmetry. With the use of the effective chiral Lagrangian for kaon--baryon and kaon--kaon interactions; coupled with the relativistic mean field theory and universal three-baryon repulsive interaction, we clarify the effects of the $s$-wave kaon--baryon scalar interaction simulated by the kaon--baryon sigma terms and vector interaction (Tomozawa--Weinberg term) on kaon properties in hyperon-mixed matter, the onset density of kaon condensation, and the equation of state with the ($Y$+$K$) phase. In particular, the quark condensates in the ($Y$+$K$) phase are obtained, and their relevance to chiral symmetry restoration is discussed.

The remarkable sensitivity achieved by the planned Laser Interferometer Space Antenna (LISA) will allow us to observe gravitational-wave signals from the mergers of massive black hole binaries (MBHBs) with signal-to-noise ratio (SNR) in the hundreds, or even thousands. At such high SNR, our ability to precisely infer the parameters of an MBHB from the detected signal will be limited by the accuracy of the waveform templates we use. In this paper, we explore the systematic biases that arise in parameter estimation if we use waveform templates that do not model radiation in higher-order multipoles. This is an important consideration for the large fraction of high-mass events expected to be observed with LISA. We examine how the biases change for MBHB events with different total masses, mass ratios, and inclination angles. We find that systematic biases due to insufficient mode content are severe for events with total redshifted mass $\gtrsim10^6\,M_\odot$. We then compare several methods of predicting such systematic biases without performing a full Bayesian parameter estimation. In particular, we show that through direct likelihood optimization it is possible to predict systematic biases with remarkable computational efficiency and accuracy. Finally, we devise a method to construct approximate waveforms including angular multipoles with $\ell\geq5$ to better understand how many additional modes (beyond the ones available in current approximants) might be required to perform unbiased parameter estimation on the MBHB signals detected by LISA.

If there exist unstable but long-lived relics of the early universe, their decays could produce detectable fluxes of gamma rays and neutrinos. In this paper, we point out that the decays of superheavy particles, $m_{\chi} \gtrsim 10^{10} \, \text{GeV}$,would produce an enhanced flux of ultra-high-energy neutrinos through the processes of muon and pion pair production in the resulting electromagnetic cascades. These processes transfer energy from electromagnetic decay products into neutrinos, relaxing the constraints that can be derived from gamma-ray observations, and increasing the sensitivity of high-energy neutrino telescopes to superheavy particle decays. Taking this into account, we derive new constraints on long-lived superheavy relics from the IceCube Neutrino Observatory, and from the Fermi Gamma-Ray Space Telescope. We find that IceCube-Gen2, and other next generation neutrino telescopes, will provide unprecedented sensitivity to the decays of superheavy dark matter particles and other long-lived relics.

It is often stated in the physics literature that maximally-spinning Kerr black-hole spacetimes are characterized by near-horizon co-rotating circular geodesics of radius $r_{\text{circular}}$ with the property $r_{\text{circular}}\to r^+_{\text{H}}$, where $r_{\text{H}}$ is the horizon radius of the extremal black hole. Based on the famous Thorne hoop conjecture, in the present compact paper we provide evidence for the existence of a non-trivial lower bound ${r_{\text{circular}-r_{\text{H}}}\over{r_{\text{H}}}}\gtrsim (\mu/M)^{1/2}$ on the radii of circular orbits in the extremal Kerr black-hole spacetime, where $\mu/M$ is the dimensionless mass ratio which characterizes the composed black-hole-orbiting-particle system.

We propose a novel cogenesis of baryon and dark matter (DM) in the Universe by utilising a first-order phase transition (FOPT) in the dark sector containing an asymmetric Dirac fermion $\chi$. Due to the mass difference of $\chi$ across the bubble walls, it is energetically favourable for $\chi$ to get trapped in the false vacuum leading to the formation of Fermi-ball, which can self-collapse to form primordial black hole (PBH) if $\chi$ has a sufficiently large Yukawa interaction. While such PBH formed out of false vacuum collapse can give rise to the DM in the Universe, a tiny amount of asymmetric $\chi$ leaking into the true vacuum through the bubble walls can transfer the dark asymmetry into the visible sector via decay. The same mass difference of $\chi$ across the two minima which decides the amount of trapping or filtering of $\chi$, also allows $\chi$ decay into visible sector in the true minima while keeping it stable in the false vacuum. Our filtered cogenesis scenario can be probed via FOPT generated stochastic gravitational waves (GW) at near future detectors in addition to the well-known detection aspects of asteroid mass PBH constituting DM in the Universe.

We develop a theory of various kinds of large-scale clustering of inertial particles in a rotating density stratified or inhomogeneous turbulent fluid flows. The large-scale particle clustering occurs in scales which are much larger than the integral scale of turbulence, and it is described in terms of the effective pumping velocity in a turbulent flux of particles. We show that for a fast rotating strongly anisotropic turbulence, the large-scale clustering occurs in the plane perpendicular to rotation axis in the direction of the fluid density stratification. We apply the theory of the large-scale particle clustering for explanation of the formation of planetesimals (progenitors of planets) in accretion protoplanetary discs. We determine the radial profiles of the radial and azimuthal components of the effective pumping velocity of particles which have two maxima corresponding to different regimes of the particle--fluid interactions: at the small radius it is the Stokes regime, while at the larger radius it is the Epstein regime. With the decrease the particle radius, the distance between the maxima increases. This implies that smaller-size particles are concentrated nearby the central body of the accretion disk, while larger-size particles are accumulated far from the central body. The dynamic time of the particle clustering is about $\tau_{\rm dyn} \sim 10^5$--$10^6$ years, while the turbulent diffusion time is about $10^7$ years, that is much larger than the characteristic formation time of large-scale particle clusters ($\sim \tau_{\rm dyn}$).

Although the QCD phase transition is a crossover in the standard model, nonstandard effects such as a large lepton asymmetry are known to make it first order, with possible applications to gravitational wave production. This process is sensitive to the speed of the bubble walls during the phase transition, which is difficult to compute from first principles. We take advantage of recent progress on wall speed determinations to provide a simple estimate which constrains the wall speed to be significantly lower than what has been used in previous literature. This in turn strongly suppresses the production of gravitational waves, to a level that is just out of reach of the most sensitive projected experiment for this signal, $\mu$Ares.

We deepen the analysis of the cosmological acceleration produced by quantum gravity dynamics in the formalism of group field theory condensate cosmology, treated at the coarse-grained level via a phenomenological model, in the language of hydrodynamics on minisuperspace. Specifically, we conduct a detailed analysis of the late-time evolution, which shows a phantom-like phase followed by an asymptotic De Sitter expansion. We argue that the model indicates a recent occurrence of the phantom crossing and we extract a more precise expression for the effective cosmological constant, linking its value to other parameters in the model and to the scale of the quantum bounce in the early universe evolution. Additionally, we show how the phantom phase produced by our quantum gravity dynamics increases the inferred value of the current Hubble parameter based on observed data, indicating a possible quantum gravity mechanism for alleviating the Hubble tension. Our results represent a concrete example of how quantum gravity can provide an explanation for large-scale cosmological puzzles, in an emergent spacetime scenario.

Microwave SQUID Multiplexing (uMUX) is a widely used technique in the low temperature detectors community as it offers high capacity of reading out large scale Transition-Edge Sensor (TES) arrays. In this paper, we propose a Sliding Flux Ramp Demodulation (SFRD) algorithm for uMUX readout system. It can achieve a sampling rate in the order of MHz while maintaining a multiplexing ratio about one thousand. Advancing of this large array readout technique makes it possible to observe scientiffc objects with improved time resolution and event count rate. This will be highly helpful for TES calorimeters in X-ray applications, such as X-ray astrophysics missions.

Meta-instability is an irreversible precursor of earthquakes. To identify the meta-instability precursor of the Yangbi $M_S$ 6.4 earthquake ($99.87^{\circ}\mathrm{E}$, $25.67^{\circ}\mathrm{N}$)that occurred on May 21, 2021, we selected seismic data from the pre-earthquake period between 1 and 21 May. We then calculated the apparent wave velocity ratio and the apparent Poisson\text{'}s ratio within the region of $98.5^{\circ} \mathrm{E}-101^{\circ}\mathrm{E}$, $24.6^{\circ}\mathrm{N}-27.1^{\circ}\mathrm{N}$ and interpolated these values. Our findings revealed that the trends of the fitted straight lines at the maximum and minimum points of the gradient divergences of the apparent wave velocity ratio and apparent Poisson\text{'}s ratio fields are consistent with the source mechanism solution for Sections 1 and 2, respectively. Similarly, the trend of the fitted straight lines at the minimum and maximum points of their values is also consistent with the source mechanism solution for Sections 1 and 2. Positive gradient divergence values indicate energy released, whereas negative values suggest energy absorption. The observed stress state matches the experimentally demonstrated meta-instable state. We propose that this method can be a reference for identifying the meta-instability of strike-slip strong earthquakes with a significant number of foreshocks. For seismically active regions, increasing the number of stations with rich data acquisition will facilitate more convenient stress field analysis.

Wormholes are exotic compact objects characterized by the absence of essential singularities and horizons, acting as slender bridges linking two distinct regions of spacetime. Despite their theoretical significance, they remain however undetected, possibly due to their ability to closely mimic the observational properties of black holes. This study explores whether a static and spherically symmetric wormhole within General Relativity can reproduce the quasi-normal mode spectrum of a Schwarzschild black hole under scalar, electromagnetic, and axial gravitational perturbations, both individually and in combination. To address this, we reformulate the wormhole metric components using a near-throat parametrization. Our analysis concentrates on the fundamental mode and first overtone, estimated via the Wentzel-Kramers-Brillouin method. By employing a customized minimization strategy, we demonstrate that within a specific region of the parameter space, a wormhole can successfully replicate a subset of the black hole quasi-normal mode spectrum.

In the present work we study cosmology in dilatonic $f(R,T)$ gravity to address the inflationary phase of the early Universe. As usual, in dilatonic gravity the scalar potential assumes the exponential form. However, this potential is not good enough to be in accord with the Planck 2018 data. More strikingly, the generalized $\beta$-exponential cannot take this into account either. It is just the presence of the dilatonic sector, in the intermediate coupling regime, that can help the theory to be in full accord with the observational data.

Resonant axion-plasmon conversion in the magnetospheres of magnetars may substantially impact the landscape of dark-matter axion detection. This work explores how resonant axion-plasmon conversion, through a mechanism that is analogous to the Mikheyev-Smirnov-Wolfenstein (NSW) effect in neutrinos, modify the expected radio signals from axion-photon conversions observed on Earth. Critically, the resonant conversion radius lies within the region expected for axion-photon conversion, introducing a nonradiative power loss that diminishes the anticipated photon flux. Our analysis demonstrates that this effect can reduce radio telescope sensitivities, shifting them into regions excluded by previous experiments. These findings compel a reassessment of experimental constraints derived from radio signatures of axion-photon conversions and highlight the necessity of accounting for plasmon effects in astrophysical axion searches. The presented corrections provide critical insights for refining the detection strategies of future telescope-based dark matter axion experiments.

We propose new versions of the slow-roll approximation for inflationary models with nonminimally coupled scalar fields. We derive more precise expressions for the standard slow-roll parameters as functions of the scalar field. To verify the accuracy of the proposed approximations, we consider inflationary models with the induced gravity term and the fourth-order monomial potential. For specific values of the model parameters, this model is the well-known Higgs-driven inflationary model. We investigate the inflationary dynamics in the Jordan frame and come to the conclusion that the proposed versions of the slow-roll approximation are not only more accurate at the end of inflation, but also give essentially more precise estimations for the tensor-to-scalar ratio $r$ and of the amplitude of scalar perturbations $A_s$.