Papers for Thursday, Jul 03 2025

In this paper, we present the first direct measurement of the proper motion of the neutron star (NS) in the supernova remnant (SNR) G18.9-1.1 using a 15-year Chandra baseline. After correcting the observations' astrometric solutions using reference Gaia stars' positions, we measure a total proper motion of 24.7 $\pm$ 6.8 mas yr$^{-1}$ at an angle of $336^\circ \pm$ 16$^\circ$ east of north. Using the distance estimates from literature of 2.1 kpc and 3.8 kpc, this proper motion corresponds to Galactic rotation-corrected transverse velocities of 264d$_{2.1}$ $\pm$ 79 km s$^{-1}$ and 474d$_{3.8}$ $\pm$ 129 km s$^{-1}$, respectively. Our power ratio method analysis of SNR ejecta slightly favors the higher velocity, as multipole moments calculated from the back-evolved center using the farther distance are more consistent with values from other CCSNRs. The NS's motion is directly opposite the motion of bulk ejecta in G18.9$-$1.1, providing yet more evidence that NS kicks are generated via a conservation of momentum-like process between the NS and the ejecta, as has been observed in other SNRs.

The LIGHTS survey is imaging galaxies at a depth and spatial resolution comparable to what the Legacy Survey of Space and Time (LSST) will produce in 10 years (i.e., $\sim$31 mag/arcsec$^2$; 3$\sigma$ in areas equivalent to 10$^{\prime\prime}$$\times$ 10$^{\prime\prime}$). This opens up the possibility of probing the edge of galaxies, as the farthest location of in-situ star formation, with a precision that we have been unable to achieve in the past. Traditionally, galaxy edges have been analyzed in one-dimension through ellipse averaging or visual inspection. Our approach allows for a two-dimensional exploration of galaxy edges, which is crucial for understanding deviations from disc symmetry and the environmental effects on galaxy growth. In this paper, we propose a novel method using the second derivative of the surface mass density map of a galaxy to determine its edges. This offers a robust quantitative alternative to traditional edge-detection methods when deep imaging is available. Our technique incorporates Wiener-Hunt deconvolution to remove the effect of the Point Spread Function (PSF) by the galaxy itself. By applying our methodology to the LIGHTS galaxy NGC 3486, we identify the edge at 205$^{\prime\prime}$ $\pm$ 5$^{\prime\prime}$. At this radius, the stellar surface mass density is $\sim$1 M$_\odot$/pc$^2$, supporting a potential connection between galaxy edges and a threshold for in-situ star formation. Our two-dimensional analysis on NGC 3486 reveals an edge asymmetry of $\sim$5$\%$. These techniques will be of paramount importance for a physically motivated determination of the sizes of galaxies in ultra-deep surveys such as LSST, Euclid and Roman.

The gravitational potential of supermassive black holes is so powerful that it triggers some of the most intense phenomena in the Universe. Accretion onto these objects and relativistic jet emission from their vicinity are observable across a wide range of frequencies and throughout cosmic history. However, despite this wealth of data, many aspects of their underlying mechanisms remain elusive. Investigating this phenomena across all frequencies is crucial, yet some energy windows are still poorly explored. One such window is the MeV energy range: many key signatures related to the emission from the SMBH environment - both in quiescent and active phases - are expected to lie between one and several hundreds MeV. In this work, we explore some of the open questions regarding the behavior and emission processes in the surroundings of SMBHs, and how these questions might be approached. From the elusive nature of Fermi bubbles around our Galactic Centre, to the origin of high-energy neutrinos in the nuclei and jets of Active Galactic Nuclei, to the nature and emission mechanisms of the most powerful blazars, the MeV window stands out as a crucial key to understanding SMBH physics.

The Helmi streams are remnants of a dwarf galaxy that was accreted by the Milky Way and whose stars now form a distinct kinematic and chemical substructure in the Galactic halo. Precisely age-dating these typically faint stars of extragalactic origin has been notoriously difficult due to the limitations of using only spectroscopic data, interferometry, or coarse asteroseismic measurements. Using observations from NASA's Transiting Exoplanet Survey Satellite, we report the detailed asteroseismic modeling of two of the brightest red giants within the Helmi streams, HD 175305 and HD 128279. By modeling the individual oscillation mode frequencies and the spectroscopic properties of both stars, we determine their fundamental properties including mass, radius, and age ($\tau$). We report $\tau = 11.16 \pm 0.91$ Gyr for HD 175305 and $\tau = 12.52 \pm 1.05$ Gyr for HD 128279, consistent with previously inferred star-formation histories for the Helmi streams and the differential chemical abundances between the two stars. With precise ages for individual stream members, our results reinforce the hypothesis that the Helmi streams' progenitor must have existed at least 12 Gyr ago. Our results also highlight that the ages of metal-poor, $\alpha$-enhanced red giants can be severely underestimated when inferred using global asteroseismic parameters instead of individual mode frequencies.

Many theoretical and observational studies have suggested that galaxy mergers may trigger enhanced star formation or active galactic nuclei (AGN) activity. We present an analysis of merging and nonmerging galaxies from $0.2 \leq z \leq 3$ in the IllustrisTNG50 simulation. These galaxies encompass a range of masses ($M_\star > 10^{8}M_\odot$), multiple merger stages, and mass ratios ($\geq1:10$). We examine the effect that galaxy mergers have on star formation and black hole accretion rates in the TNG50 universe. We additionally investigate how galaxy and black hole mass, merger stage, merger mass ratio, and redshift affect these quantities. Mergers in our sample show excess specific star formation rates (sSFR) at $z \leq 3$ and enhanced specific black hole accretion rates (sBHAR) at $z \lesssim 2$. The difference between sSFRs and sBHARs in the merging sample compared to the non-merging sample increases as redshift decreases. Additionally, we show that these enhancements persist for at least $\sim1$ Gyr after the merger event. Investigating how mergers behave in the TNG50 simulation throughout cosmic time enables both a better appreciation of the importance of spatial resolution in cosmological simulations and a better basis to understand our high-$z$ universe with observations from \textit{JWST}.

We investigated universal relations for compact stars rotating at the Keplerian (mass-shedding) limit, which is highly relevant for understanding the rapidly rotating objects formed in the aftermath of a neutron star-neutron star merger. Our analysis is based on a set of nucleonic EoS featuring systematic variations in the symmetry energy slope parameter $L_{\rm sym}$ and the isoscalar skewness parameter $Q_{\rm sat}$, varied within ranges that are broadly consistent with current laboratory and astrophysical constraints. The global observable properties of isolated maximally rotating stars are examined, focusing on the mass-radius relation, moment of inertia, quadrupole moment, and the Keplerian (maximum) rotation frequency, as well as their variations in the $L_{\rm sym}$-$Q_{\rm sat}$ parameter space. Next, we demonstrate that, in the limit of Keplerian rotation, universal relations remain valid across the same set of EoSs characterized by varying $L_{\rm sym}$ and $Q_{\rm sat}$. In particular, we present explicit results for the moment of inertia and quadrupole moment as functions of compactness, as well as for the moment of inertia - quadrupole moment relation. All of these relations exhibit excellent universality, with deviations typically within a few percent and rarely exceeding 10\% across a wide range of parameters. These findings support the applicability of $I$-Love-$Q$-type universal relations in observational modeling of maximally rotating compact stars and the gravitational wave emitted by them, when accounting for significant variation in the symmetry energy and high-density behavior of the nuclear EoS.

The astrophysical origin of elements synthesized through the rapid neutron capture process ($r-$process) is a long standing mystery. The hot and dense environments of core-collapse supernovae have been suggested as potential $r-$process sites, particularly the neutrino-driven wind from the newly-born protoneutron star (PNS). Wind models that neglect the potential effects of strong magnetic fields and/or rapid rotation of the PNS typically fail to achieve the necessary conditions for production of the third $r-$process peak, but robustly produce a limited or weak $r-$process for neutron-rich winds. Axisymmetric magnetohydrodynamic simulations of rotating and non-rotating PNS winds with magnetar-strength fields reveal that high entropy material is quasi-periodically ejected from the equatorial closed zone of the PNS magnetosphere. Here, we post-process tracer particle trajectories from these simulations using a nuclear reaction network in order to explore the resulting nucleosynthesis across a range of PNS magnetic field strengths, rotation rates, and neutrino luminosities (cooling phase after core-bounce). We find that a robust $r-$process up to and beyond the third peak is generic to magnetar birth, even for magnetic fields as weak as $\sim 5\times 10^{14}$ G. Depending on the distribution of magnetic field strengths and rotation at birth, we estimate that magnetized PNS winds could account for $\sim 5-100\%$ of the Galactic $r-$process inventory, extending up to the third peak. The robust $r-$process in our calculations is accompanied by overproduction of elements with mass number $\rm A\lesssim 120$ compared to the Solar abundances. We also find that $^{92}\rm Mo$ (a $p-$isotope) is produced in significant quantities in neutron-rich winds.

Cosmography is a widely applied method to infer kinematics of the Universe at small cosmological scales while remaining agnostic about the theory of gravity at play. Usually cosmologists invoke the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric in cosmographic analyses, however generalised approaches allow for analyses outside of any assumed geometrical model. These methods have great promise to be able to model-independently map the cosmic neighborhood where the Universe has not yet converged to isotropy. In this regime, anisotropies can bias parameter inferences if they are not accounted for, and thus must be included for precision cosmology analyses, even when the principle aim is to infer the background cosmology. In this paper, we develop a method to predict the dipole in luminosity distances that arises due to nearby inhomogeneities. This is the leading-order correction to the standard isotropic distance-redshift law. Within a very broad class of general-relativistic universe models, we provide an interpretation of the dipole in terms of the gradients in expansion rate and density which is free from any underlying background cosmology. We use numerical relativity simulations, with improved initial data methods, alongside fully relativistic ray tracing to test the power of our prediction. We find our prediction accurately captures the dipole signature in our simulations to within ~10% for redshifts $z\lesssim 0.07$ in reasonably smooth simulations. In the presence of more non-linear density fields, we find this reduces to $z\lesssim 0.02$. This represents up to an order of magnitude improvement with respect to what is achieved by naive, local cosmography-based predictions. Our paper thus addresses important issues regarding convergence properties of anisotropic cosmographic series expansions that would otherwise limit their applicability to very narrow redshift ranges.

The escape fraction of Lyman continuum photons (fesc(LyC)) is the last key unknown in our understanding of cosmic reionization. Directly estimating the escape fraction (fesc) of ionizing photons in the epoch of reionization (EoR) is impossible, due to the opacity of the intergalactic medium (IGM). However, a high fesc leaves clear imprints in the spectrum of a galaxy, due to reduced nebular line and continuum emission, which also leads to bluer UV continuum slopes (betaUV). Here, we exploit the large archive of deep JWST/NIRSpec spectra from the DAWN JWST Archive to analyze over 1'400 galaxies at 5 < zspec < 10 and constrain their fesc based on SED fitting enhanced with a picket fence model. We identify 71 high-confidence sources with significant fesc based on Bayes factor analysis strongly favouring fesc > 0 over fesc = 0 solutions. We compare the characteristics of this high-escape subset against both the parent sample and established diagnostics including betaUV slope, O32, and SFR surface density (SigmaSFR). For the overall sample, we find that most sources have a low escape fraction (<1%), however, a small subset of sources seems to emit a large number of their ionizing photons into the IGM, such that the average fesc is found to be ~10%, as needed for galaxies to drive reionization. Although uncertainties remain regarding recent burstiness and the intrinsic stellar ionizing photon output at low metallicities, our results demonstrate the unique capability of JWST/NIRSpec to identify individual LyC leakers, measure average fesc and thus constrain the drivers of cosmic reionization.

High-cadence microlensing observations uncovered a population of very short-timescale microlensing events, which are believed to be caused by the population of free-floating planets (FFP) roaming the Milky Way. Unfortunately, the light curves of such events are indistinguishable from those caused by wide-orbit planets. To properly differentiate both cases, one needs high-resolution observations that would allow resolving a putative luminous companion to the lens long before or after the event. Usually, the baseline between the event and high-resolution observations needs to be quite long ($\sim 10$ yr), hindering potential follow-up efforts. However, there is a chance to use archival data if they exist. Here, we present an analysis of the microlensing event OGLE-2023-BLG-0524, the site of which was captured in 1997 with the Hubble Space Telescope (HST). Hence, we achieve a record-breaking baseline length of 25 years. A very short duration of the event ($t_E = 0.346 \pm 0.008$ d) indicates an FFP as the explanation. We have not detected any potential companion to the lens with the HST data, which is consistent with the FFP origin of the event. Thanks to the available HST data, we are able to reject from 25% to 48% of potential stellar companions depending on the assumed population model. Based on the finite-source effects in the light curve we measure the angular Einstein radius value $\theta_E = 4.78 \pm 0.23 \mu as$, suggesting a super-Earth in the Galactic disk or a sub-Saturn-mass planet in the Galactic bulge. We show that the archival high-resolution images should be available for several microlensing events, providing us with the unprecedented possibility of seeing the lensing system as it was many years before the event.

We present a comprehensive study of five nearby active galaxies featuring large (tens of kpc) extended emission-line regions (EELRs). The study is based on large-format integral field spectroscopic observations conducted with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). The spatially resolved kinematics of the ionized gas and stellar components show signs of rotation, bi-conical outflows, and complex behavior likely associated with past interactions. Analysis of the physical conditions of the EELRs indicates that in these systems, the active galactic nucleus (AGN) is the primary ionization source. Using radiative transfer simulations, we compare the ionization state across the EELRs to estimate the required AGN bolometric luminosities at different radial distances. Then, considering the projected light travel time, we reconstruct the inferred AGN luminosity curves. We find that all sources are consistent with a fading trend in intrinsic AGN luminosity by 0.2--3 dex over timescales of 40,000--80,000 years, with a time dependence consistent with previous studies of fading AGNs. These results support the hypothesis that most AGN undergo significant fluctuations in their accretion rates over multiple timescales ranging from 10,000 to 1,000,000 years, as proposed by existing theoretical models. These results provide new insights into the transient phases of AGN activity at previously unexplored scales and their potential long-term impact on their host galaxies through various feedback mechanisms.

We investigate the star formation process across the disk of M33 using a multiwavelength dataset and disk dynamics. We computed numerically equilibrium values of gas densities and scale heights across the disk, taking into account dark matter and testing several analytic approximations that are often used to estimate these variables and the hydrostatic pressure. Orthogonal regressions and hierarchical Bayesian models, as well as random forest (RF) analyses, were used to establish the fundamental relations at physical scales from 160~pc to 1~kpc. The gas pressure, is the main driver of the star formation rate (SFR) surface density throughout the whole star-forming disk of M33. High-pressure regions enhance the atomic-to-molecular gas conversion, with the molecular hydrogen mass surface density being tightly correlated to pressure and a uniform scaling law throughout the M33 disk. The relation between pressure and SFR surface density differs, showing a change in slope from the inner to the outer disk. Scaling laws do not depend on the physical scale and brings out an intrinsic scatter linked to variations in the efficiency and relative age of the molecular gas-to-stars conversion. In the inner disk, where spiral arms are present and the stellar surface density dominates gravity, the pressure and SFR surface densidy establish an almost linear correlation with a smaller dispersion than that of the molecular gas -- SFR surface density relation. In the atomic gas-dominated outer disk, the SFR density has a steeper dependence on pressure, which we propose could be the result of an increasing fraction of diffuse molecular gas that does not form stars.

Magnetospheric accretion is a key process that shapes the inner disks of T Tauri stars, controlling mass and angular momentum evolution. It produces strong ultraviolet and optical emission that irradiates the planet-forming environment. In this work, we characterize the magnetospheric geometries, accretion rates, extinction properties, and hotspot structures of 67 T Tauri stars in the largest and most consistent study of ultraviolet and optical accretion signatures to date. To do so, we apply an accretion flow model to velocity-resolved H$\alpha$ profiles for T Tauri stars from the HST/ULLYSES program with consistently-derived stellar parameters. We find typical magnetospheric truncation radii to be almost half of the usually-assumed value of 5 stellar radii. We then model the same stars' HST/STIS spectra with an accretion shock model, finding a diverse range of hotspot structures. Phase-folding multi-epoch shock models reveals rotational modulation of observed hotspot energy flux densities, indicative of hotspots that persist for at least 3 stellar rotation periods. For the first time, we perform a large-scale, self-consistent comparison of accretion rates measured using accretion flow and shock models, finding them to be consistent within $\sim$0.16 dex for contemporaneous observations. Finally, we find that up to 50% of the total accretion luminosity is at short wavelengths accessible only from space, highlighting the crucial role of ultraviolet spectra in constraining accretion spectral energy distributions, hotspot structure, and extinction.

Solar flare ribbons are features in the lower atmosphere that show enhanced emission intensity, particularly noticeable in H$\alpha$ and EUV images of the transition region and upper chromosphere. In this study, we provide observational evidence that the flare ribbons are composed of a large number of coherent sub-structures that we call \textit{riblets}. Our observations, taken with the Swedish 1-m Solar Telescope / CRisp Imaging Spectro-Polarimeter (SST/CRISP) instrument, focused on an X1.5-class flare that occurred on June 10, 2014. Using data obtained from SST/CRISP, we utilize CRISPEX, an interactive graphical user interface tool, to individually track and analyze the dynamic evolution of these riblets, gaining initial insights into their characteristics. Our investigation reveals the existence of two distinct categories of riblets within the flare ribbons. The first category comprises riblets that evolve perpendicular to the ribbon towards the solar surface at a constant velocity over their lifespan. We present one such riblet with a velocity of 18.04 km/s, persisting for a duration of 67 seconds. The second category includes riblets with varying velocities, indicating a more complex and non-linear evolution pattern. We present one such riblet with an average acceleration of 0.70 km/s$^2$ and a duration of 22.8 seconds. Riblets, with their distinct characteristics, provide an opportunity to study the chromospheric response to a solar flare in a more detailed manner.

We conducted a multi-wavelength monitoring campaign of the blazar S5 0716+714 from 2022 November 26 to 2023 May 28 using optical telescopes in Egypt and Bulgaria. Data were taken during 84 nights in 11 of which intranight monitoring was performed. We also use optical and gamma-ray survey data. On long-term time-scales, we find a gradual decrease of the S5 0716+714 activity since the beginning of 2020 in both optical and gamma-rays. On short-term time-scales, the individual optical light curves are strongly correlated among each other with no time lags observed. The V-band percentage variability amplitude equals 97.59 +/- 0.02 per cent. We find moderate flatter-when-brighter spectral behaviour with the strength of the 'spectral index - flux' anti-correlation decreasing towards the longer wavelengths. The main feature of the short-term light curves is a transient quasi-periodic oscillation with a period of 43.5 +/- 3.6 d. The V-band light curve is modelled with two helically moving blobs and a synchrotron flare. We estimate the resulting parameters, as well as limits on the radius, magnetic field strength, and electron Lorentz factor of the region responsible for the flare. On intranight time-scales, we find smooth flux variations with no flares and derive a duty cycle in the range ~10-20 per cent. The lack of flares on intranight time-scales could result from a temporarily homogeneous jet flow without formation of turbulent cells in terms of prevented Kelvin-Helmholtz instability. The analysis of the data reveals a low activity of S5 0716+714 on all time-scales during the observation period.

The gravitational wave event GW230529 was the first compact binary merger observation reported in the fourth observing run of the LIGO--Virgo--KAGRA observatory network. The more massive component of the GW230529 source binary was inferred to contain a compact object in the lower mass gap, a purported gap between the most massive neutron stars and least massive black holes based on compact object observations in the Milky Way. However, the properties of the mass-gap object, such as its exact mass and spin, remain unresolved due to statistical uncertainties. To investigate the properties of the GW230529 source binary, we perform parameter estimation on a suite of simulated gravitational wave systems with similar parameters. We vary the intrinsic properties of the simulated systems, the detector noise properties, the signal-to-noise ratios, and the waveform model used in recovery. We find that the low signal-to-noise ratio of GW230529 is the key reason for the ambiguity in determining whether the mass of the primary object in the binary is consistent with a low-mass black hole or a high-mass neutron star, since the priors for the masses and spins have a significant impact on the posterior distribution. The inclusion of tidal effects in the waveform model also contributes to the observed degeneracies in the posteriors. We show that a higher signal-to-noise ratio in future observations of mass-gap gravitational wave sources will increase the precision of the measurements and allow us to determine whether the components are black holes or neutron stars.

The occurrence rate of cold Jupiters was found to depend on stellar mass. The formation environment in the protoplanetary disks regulates core formation and the subsequent gas accretion. In this study, we simulate giant planet formation via pebble accretion accounting for various stellar masses, core formation times, disk turbulent viscosities, and grain opacities. We use a self-consistent formation model that calculates the solid accretion rate and gas accretion rate of growing protoplanets. We investigate how the planetary formation, in particular, the contraction of the envelope, and the formation timescale change under different conditions. We find that to reproduce the observed occurrence rate of cold Jupiters, giant planets must undergo slow envelope contraction after they reach pebble isolation, which lasts for several Myrs. Such a slow contraction phase can be achieved when the grain opacity is assumed to be as high as that of the interstellar medium (ISM). If the grain opacity is smaller than the ISM opacity by a factor of ten or more, the growing protoplanets reach crossover mass within 3 Myrs and form too many cold Jupiters around stars of >0.4Msun. Protoplanets around low-mass stars <0.4Msun take >10 Myrs to reach crossover mass also with low grain opacity. If the grain opacity in the planetary envelope is much lower than that of ISM, other mechanisms, such as atmospheric recycling or planetesimal accretion, is required for cold Jupiter formation. We next explore how the deposition of the accreted heavy elements to the planetary envelope changes the formation timescale. Our model suggests that the formation timescale could be longer due to heavy-element enrichment, resulting from the lower core mass at pebble isolation. We conclude that the details of the formation processes have a significant effect on the planetary growth and therefore, the formation of gaseous planets.

We have carried out a morphological search for molecular clouds possibly associated with 48 Galactic infrared bubbles with angular radii of $>1'$ in the southern Galactic plane of $295^\circ \le l\le 350^\circ$ and $|b|\le 1^\circ$ presented by Hanaoka et al. (2019). 116 molecular clouds in the $(l,b,V_{\rm LSR})$ space are identified from the archival $^{12}$CO~$J$~=~1--0 line data obtained by the Mopra Southern Galactic plane survey, where $V_{\rm LSR}$ is the CO-line radial velocity. The kinematic distances are derived from $V_{\rm LSR}$ using the most accurate rotation curve of the Milky Way. We also present measurements of velocity dispersion, size parameter, molecular mass, and virial mass of 116 molecular clouds. This catalog is presented to investigate star formation and the origin of molecular shells and cavities, possibly associated with infrared bubbles. CO line intensity maps and position-velocity diagrams of the molecular clouds are available online as supplementary data.

Spectral modeling codes that estimate photometric redshifts (photo-$z$) are a powerful and often reliable method for determining redshifts of galaxies. However, there are notable instances where degeneracies in spectral energy distribution (SED) colors lead to `catastrophic' failures. We highlight the case of COSBO-7, a dusty, intermediate-$z$ galaxy that masqueraded as a high-$z$ source, because it demonstrates a unique scenario where photo-$z$ codes run into issues despite extensive multi-wavelength photometry. We advocate that photo-$z$ fitting should aim to: (1) use the entire available SED (UV--radio) whenever possible to help break color degeneracies, (2) allow flexible dust attenuation prescriptions, both in terms of the attenuation curve slope and a varying 2175Å absorption feature, and (3) implement uncertainty floors to account for limitations in spectral models and also on the photometry itself.

Switchbacks (SBs) are localized magnetic field deflections in the solar wind, marked by abrupt changes in the magnetic field direction relative to the ambient solar wind. Observations onboard Parker Solar Probe (PSP) at heliocentric distances below 50 Solar Radii (Rs) showed that within SBs, perturbations in the magnetic field ({\Delta}B) and the bulk solar wind velocity ({\Delta}V) align, i.e., {\Delta}B~{\Delta}V, producing enhanced radial velocity spikes. In this study, we examine the characteristics of SB boundaries, with particular attention to the role of boundary shear flow instabilities (Kelvin-Helmholtz instability - KHI) for surface wave phenomena based on the in situ magnetic field, plasma speed, and plasma density measurements from PSP. The results indicate that SB boundaries can be unstable for generating KHI-driven surface waves, suggesting that the wave activity observed at SB boundaries is caused by shear flow instabilities. In addition, the continued development of KHI may lead to boundary erosion, contributing to the radial evolution of SBs via structural weakening or broadening. However, when {\Delta}B and {\Delta}V are closely aligned, the boundary remains stable unless the velocity shear significantly exceeds the magnetic shear. Since the observed velocity shear typically ranges from 40% to 90% of magnetic shear, the instability condition is generally not satisfied. Thus, the configuration leading to the instability arises from deviations from precise alignment of {\Delta}B and {\Delta}V in the young solar wind, and the release of the KHI presumably leads to the formation of the {\Delta}B and {\Delta}V alignment observed at SB boundaries located at 35-55 Rs.

Due to its proximity, Markarian 421 is one of the most extensively studied jetted active galactic nuclei. Its spectral energy distribution and light curve are widely studied, serving as primary means for understanding jet radiation mechanisms. Numerous intriguing observational results have been discovered, some of which, such as the hard X-ray excess, and the associated variability between X-ray and very-high-energy (VHE) emissions, challenge the commonly adopted one-zone leptonic model. In this work, by establishing a time-dependent leptohadronic model, we explore whether the hard X-ray excess and the associated variability between X-ray and VHE emissions could be interpreted by emission from hadronic interactions. Our modeling finds that for the hard X-ray excess found in 2013, both of the secondary emissions from photohadronic and hadronuclear interactions could be a possible explanation for the hard X-ray excess without introducing a super-Eddington jet power. The emission from the photohadronic interactions contributes only to the hard X-ray band, while the hadronuclear interactions also predict VHE emissions associated with the hard X-rays. While for the hard X-ray excess found in 2016, only the secondary emissions from photohadronic interactions provide an interpretation at the cost of introducing a super-Eddington jet power. For the associated variability between X-ray and VHE emissions in 2017, we find that hadronic interactions fail to provide a possible interpretation.

Fast radio bursts (FRBs) are extragalactic radio transients that offer valuable insight of intergalactic medium (IGM). However, the dispersion measure (DM) contributed by IGM ($\rm{DM_{IGM}}$) is degenerated with that from the host galaxy ($\rm{DM_{host}}$), necessitating calibration of the $\rm{DM_{IGM}}$$-z$ relation for cosmological applications. As $\rm{DM_{host}}$ is expected to correlate with host galaxy properties, it is feasible to estimate $\rm{DM_{host}}$ from observable host characteristics. In this study, we conduct spectral energy distribution (SED) and Sérsic model fittings to derive the parameters of FRB host galaxies. Then, we examine the correlations between the excess dispersion measure ($\rm{DM_{exc}}$) and host galaxy parameters, including star formation rate (SFR), stellar mass, specific star formation rate (sSFR), inclination angle, and projected area. A tight correlation between $\rm{DM_{exc}}$ and sSFR is found. This correlation is utilized to estimate the $\rm{DM_{host}}$ of FRBs, providing a method to calibrate the DM$_{\rm IGM}-z$ relation. This approach leads to a notable improvement in calibration performance.

In this Letter, we simulate the collision between outflows from the tidal disruption of a 1M$_\odot$ main sequence star around a $10^6$M$_\odot$ black hole and an initially spherically symmetric circumnuclear cloud. We launch super-Eddington outflows self-consistently by simulating the disruption of stars on both bound and unbound initial orbits using general relativistic smoothed particle hydrodynamics. We find shocks formed as early as $\sim 10~$days after the initial stellar disruption produce prompt radio emission. The shock radius ($\approx~10^{17}$~cm), velocity ($\sim 0.15$c) and total energy ($\sim 10^{51}$ erg) in our simulations match those inferred from radio observations of tidal disruption events (TDEs). We ray-trace to produce synthetic radio images and spectra to compare with the observations. While the TDE outflow is quasi-spherical, the synchrotron emitting region is aspherical but with reflection symmetry above and below the initial orbital plane. Our synthetic spectra show continuous decay in peak frequency, matching prompt radio TDE observations. Our model supports the hypothesis that synchrotron radio flares from TDEs result from the collision between outflows and the circumnuclear material.

Understanding the relationship between supermassive black holes (SMBHs) and their host galaxies at different redshifts is crucial for unraveling the processes of SMBH-galaxy co-evolution. We present the properties of nine type 1 Active Galactic Nuclei (AGNs) at intermediate redshift ($2<z<4$) using the JWST Advanced Deep Extragalactic Survey (JADES). All of them show the significant $\mathrm{H\alpha}$ broad line and the AGN contribution in spectral energy distribution. Our sample covers SMBH masses of $10^{6.1-8.2}\ M_\odot$ and stellar masses of $10^{9.3-11.0}\ M_\odot$, comparable to those of the AGNs observed in the local universe. In the low-mass SMBH regime ($<10^{8}\ M_\odot$), the BH-to-stellar mass ratios in our sample ($0.01-0.1\%$) differ from those of the AGNs at $z>4$ ($1-10\%$), suggesting that black holes and galaxies may trace different evolutionary pathways at intermediate and high redshift. We also perform 2D image decomposition using GALFIT to constrain the bulge mass by evaluating the bulge contribution in the rest-frame near-infrared flux. We identify the AGNs with low BH-to-bulge mass ratios compared to those observed in the nearby bulge-dominant galaxies. This finding suggests the existence of a galaxy-first evolutionary path, in which bulge formation occurs before substantial gas is efficiently accreted onto the central engine.

Cosmic reionization marks a critical epoch when the first galaxies ionized the intergalactic medium through the escape of Lyman continuum (LyC) radiation. Young, massive star clusters are believed to be the primary LyC sources, yet the physical mechanisms enabling LyC escape remain poorly understood. Most existing studies rely on spatially integrated observations, which lack the resolution to resolve internal galaxy structure and pinpoint where and how LyC photons escape. To address this, we propose a science case for the Habitable Worlds Observatory (HWO) that enables spatially resolved spectroscopy of LyC-emitting star clusters and their environments in low-redshift galaxies. This requires a UV integral field unit (IFU) with coverage down to ~ 900 Angstrom and a spatial resolution of 10-100 pc-capabilities essential for directly detecting LyC escape and mapping the surrounding interstellar medium. With such instrumentation, we will map cluster-scale LyC escape fractions, characterize the physical conditions of the surrounding interstellar medium, and directly observe feedback-driven outflows that facilitate LyC leakage. These observations will enable novel calibrations of indirect LyC indicators at unprecedented spatial resolution and establish direct connections between local LyC processes and those in high-redshift, clumpy star-forming galaxies. In the long run, this program will build the physical framework needed to understand how galaxies reionized the early universe and shaped its subsequent evolution.

We study the possibility that Type Ia supernovae might be produced by binary systems where the companion of the exploding white dwarf is an M-dwarf star. Such companion would appear as a runaway star, retaining its pre-explosion orbital velocity along with a kick imparted by the supernova ejecta. It might be rapidly rotating, from being tidally locked with the white dwarf prior to explosion in a very close binary. For this study, we perform a series of multidimensional hydrodynamic simulations to investigate the interaction between M-dwarf companions and SN ejecta, followed by post-impact stellar evolution modeling using the MESA code. Our initial models in the 3D simulations had high spin angular momenta and the effects of magnetic braking have been included. They very significantly reduce the final rotation. A surviving companion candidate, MV-G272, has recently been discovered in the supernova remnant G272.2-3.2, which is an 8.9$\sigma$ proper motion outlier, although being slowly rotating. Our results show that the properties of this companion (luminosity, effective temperature, surface gravity) can be reproduced by our post-impact M-dwarf models. The slow rotation, which is a common characteristic with several proposed hypervelocity SN companions, can be explained by magnetic braking during the post-impact evolution, thus supporting the possibility that the MV-G272 star is the surviving companion of the Type Ia supernova now found as G272.2-3.2 SNR.

Protoplanets can interact with their natal disks and generate gas and dust substructures such as gaps and rings. However, how these planet-induced substructures affect the disk temperature, and how that in turn influences the substructures, remains unclear. We aim to study disk substructures and the thermal structure self-consistently and explore their impact on volatile distribution. To this end, we perform iterative multi-fluid hydrodynamical and radiative transfer simulations of planet-disk interactions. We find that the temperature in a structured disk deviates significantly from that of a smooth disk due to giant planet formation. In particular, midplane temperatures in gaps can increase by tens of Kelvin, leading to volatile sublimation as well as radial shifts and multiplication of icelines. Comparing our multi-dust models with previous gas-only models, we find that the former produces slightly shallower gaps and temperatures about 10 K ($\sim25\%$) higher. Furthermore, the temperature at dust rings formed by pressure bumps can drop by several Kelvin, creating volatile freeze-out regions. Nevertheless, the overall midplane ice distribution is not strongly sensitive to whether dust is included. We also investigate the effect of varying disk viscosity. Increasing $\alpha$ viscosity from $10^{-4}$ to $10^{-2}$ leads to a roughly 10 K ($\sim25\%$) warmer midplane due to enhanced vertical dust mixing. However, higher viscosity suppresses gap opening and reduces the temperature enhancement within gaps. As a result, iceline locations do not follow a simple trend with viscosity. Finally, we propose an observational strategy using ALMA to test our predicted temperature changes within disk gaps.

We present analyses of the early data from Rubin Observatory's Data Preview 1 (DP1) for the globular cluster 47 Tuc field. The DP1 dataset for 47 Tuc includes four nights of observations from the Rubin Commissioning Camera (LSSTComCam), covering multiple bands (ugriy). We address challenges of crowding near the cluster core and toward the SMC in DP1, and demonstrate improved star-galaxy separation by fitting fifth-degree polynomials to the stellar loci in color-color diagrams and applying multi-dimensional sigma clipping. We compile a catalog of 3,576 probable 47 Tuc member stars selected via a combination of isochrone, Gaia proper-motion, and color-color space matched filtering. We explore the sources of photometric scatter in the 47 Tuc color-color sequence, evaluating contributions from various potential sources, including differential extinction within the cluster. Finally, we recover five known variable stars, including three RR Lyrae and two eclipsing binaries. Although the DP1 lightcurves have sparse temporal sampling, they appear to follow the patterns of densely-sampled literature lightcurves well. Despite some data limitations for crowded-field stellar analysis, DP1 demonstrates the promising scientific potential for future LSST data releases.

Spectroscopic observations are essential for confirming associations, measuring kinematics, and determining stellar populations in dwarf galaxies. Here, we present Keck Cosmic Web Imager (KCWI) spectra for 12 MATLAS survey dwarfs. For 9, we confirm recession velocities consistent with their literature-assumed host galaxies. We propose revisions of the host galaxy associations for MATLAS-631, 1494, and 1938. For MATLAS-1494, our measured redshift reclassifies it from an ultra-diffuse galaxy candidate to a dwarf galaxy that is of smaller physical size and places it in the field. It also appears old and passive, providing a challenge to models that invoke quenching by tidal effects. Additionally, we measure stellar population estimates for 7 of the 12 galaxies, finding a 'mixed bag' of old quenched galaxies and those that are currently forming stars. Compared to the literature we find generally younger ages and higher metallicities. This result may help reconcile the observed offset of MATLAS survey dwarf galaxies from the universal stellar mass-metallicity relationship reported by Heesters et al. (2023).

Some recent pulsar observations cannot naturally fit into the conventional picture of neutron stars: the compact objects associated with HESS J1731-347 and XTE J1814-338 have too small radii at the low-mass regime, while the secondary component of GW190814 is too massive for neutron stars to be compatible with constraints from the GW170817 event. In this study, we demonstrate that all these anomalous observations and tensions, together with other conventional ones such as recent NICER observations of PSR J0740+6620, J0030+0451, and PSR J0437-4715, can be naturally explained simultaneously by a new general type of hybrid stars that are self-bound and radially stable in the slow phase transition context, and by some subsets even in the rapid phase transition context also. As a proof of concept, we use hybrid quark stars, inverted hybrid stars, and hybrid strangeon stars as benchmark examples to explicitly demonstrate the advantage and feasibility of slow stable self-bound hybrid stars in relieving all tensions related to compact stars' masses, radii, and tidal deformabilities.

Photometric redshifts of the source galaxies are a key source of systematic uncertainty in the Rubin Observatory Legacy Survey of Space and Time (LSST)'s galaxy clustering and weak lensing analysis, i.e., the $3\times 2$pt analysis. This paper introduces a Fisher forecast code FisherA2Z for the LSST Year 10 (Y10) $3 \times 2$pt and cosmic shear analyses, utilizing a 15-parameter redshift distribution model, with one redshift bias, variance, and outlier rate per tomographic bin. FisherA2Z employs the Core Cosmology Library CCL to compute the large-scale structure power spectrum and incorporates a four-parameter nonlinear alignment model for intrinsic alignments. We evaluate the impact of marginalizing over redshift distribution parameters on weak lensing, forecast biases in cosmological parameters due to redshift errors, and assess cosmological parameter sensitivity to redshift systematic parameters using decision trees. The sensitivity study reveals that for LSST $3\times2$pt analysis, $S_8$ is most sensitive to the mean redshift of the fourth out of the five source tomographic bins, while other cosmological parameters possess different sensitivities. Additionally, we provide cosmological analysis forecasts based on different scenarios of spectroscopic training datasets. We find that the figures-of-merit for the cosmological results increase with the number of spectroscopic training galaxies, and with the completeness of the training set above $z=1.6$, assuming the redshift information comes solely from the training set galaxies without other external constraints.

The combination of galaxy-galaxy weak lensing and galaxy clustering is a powerful probe of the cosmological model, and exploration of how to best model and extract this information from the signals is essential. We present the measurement of the galaxy-galaxy weak lensing signals using the SDSS DR11 spectroscopic galaxies as lens galaxies, and the HSC Y3 shear catalog as source galaxies, binned into four tomographic bins by their photometric redshift. The SDSS DR11 galaxies, with a redshift range $0.15<z<0.7$, are binned into three redshift bins, each as a probe for measuring the projected correlation function, $w_p(R_p)$. We measure the galaxy-galaxy lensing signal $\Delta \Sigma (R_p)$ in 12 lens-source bin pairs and show that there is no evidence for significant systematic biases in the measurement with null testing. We combine our $w_p(R_p)$ and $\Delta \Sigma (R_p)$ ($2\times2$pt) data vectors and perform likelihood inference with a flat $\Lambda$CDM model. For $\Delta \Sigma (R_p)$, we extend the lower limit of the scale cut compared to previous HSC Y3 analyses to $2 h^{-1}$Mpc by including a point-mass correction term in addition to the minimal bias model. We present various tests to validate our model and provide extended consistency tests. In the $\Lambda$CDM context, our fiducial model yields $S_8 = 0.804^{+0.051}_{-0.051}$. The $2\times2$pt data vector provides redshift parameter constraints for the third and fourth redshift bins $\Delta z_3 = -0.079^{+0.074}_{-0.084}$, and $\Delta z_4 = -0.203^{+0.167}_{-0.206}$, which is consistent with results from the previous tomographic cosmic shear studies, and serves as the foundation for a future $3\times 2$pt analysis.

We present a detailed comparison of $\Lambda$CDM, $\omega$CDM, and $\omega_0\omega_a$CDM cosmologies using the latest redshift-resolved BAO measurements from DESI DR2. By separating $r_d$-dependent and $r_d$-independent observables and applying both fixed and model-derived sound horizon values, we test the consistency of each model across 24 observables spanning redshift 0.2 to 2.3. Our results show that the $\omega$CDM model is statistically favored over $\Lambda$CDM, with lower $\chi^2$, AIC, and BIC values and reduced residuals, especially at $z = 0.706$ and 0.934. In contrast, the $\omega_0\omega_a$CDM model is strongly disfavored, even when using model-specific $r_d$, due to large residuals caused by the unphysical high-redshift behavior.

The combination of galaxy clustering and weak lensing is a powerful probe of the cosmology model. We present a joint analysis of galaxy clustering and weak lensing cosmology using SDSS data as the tracer of dark matter (lens sample) and the HSC Y3 dataset as source galaxies. The analysis divides HSC Y3 galaxies into four tomographic bins for both galaxy-galaxy lensing and cosmic shear measurements, and employs a point-mass correction model to utilize galaxy-galaxy lensing signals down to 2$h^{-1}$Mpc, extending up to 70$h^{-1}$Mpc. These strategies enhance the signal-to-noise ratio of the galaxy-galaxy lensing data vector. Using a flat $\Lambda$CDM model, we find $S_8 = 0.780^{+0.029}_{-0.030}$, and using a $w$CDM model, we obtain $S_8 = 0.756^{+0.038}_{-0.036}$ with $w = -1.176^{+0.310}_{-0.346}$. We apply uninformative priors on the redshift mean-shift parameters for the third and fourth tomographic bins. Leveraging the self-calibration power of tomographic weak lensing, we measure $\Delta z_3 = -0.112^{+0.046}_{-0.049}$ and $\Delta z_4 = -0.185^{+0.071}_{-0.081}$, in agreement with previous HSC Y3 results. This demonstrates that weak lensing self-calibration can achieve redshift constraints comparable to other methods such as photometric and clustering redshift calibration.

Fourier neural operators (FNOs) provide a mesh-independent way to learn solution operators for partial differential equations, yet their efficacy for magnetized turbulence is largely unexplored. Here we train an FNO surrogate for the 2-D Orszag-Tang vortex, a canonical non-ideal magnetohydrodynamic (MHD) benchmark, across an ensemble of viscosities and magnetic diffusivities. On unseen parameter settings the model achieves a mean-squared error of $\approx 6 \times 10^{-3}$ in velocity and $\approx 10^{-3}$ in magnetic field, reproduces energy spectra and dissipation rates within $96\%$ accuracy, and retains temporal coherence over long timescales. Spectral analysis shows accurate recovery of large- and intermediate-scale structures, with degradation at the smallest resolved scales due to Fourier-mode truncation. Relative to a UNet baseline the FNO cuts error by $97\%$, and compared with a high-order finite-volume solver it delivers a $25\times$ inference speed-up, offering a practical path to rapid parameter sweeps in MHD simulations.

In this paper, we present for the first time a comprehensive statistical study between type II radio bursts from the metric (m) to the dekameric-hectometric (DH) domain and their associated solar and space weather (SW) phenomena, namely, solar flares (SFs), sunspot (SN) configurations, filament eruptions, coronal mass ejections (CMEs), their interplanetary (IP) counterparts (ICMEs) and shocks, in situ detected particles and geomagnetic storms (GSs). The m-only and m+DH radio signatures are identified from dynamic spectra provided by the ground-based RSTN stations distributed over the globe together with Wind/WAVES satellite data. The DH-only type IIs are adopted from a ready catalog based on Wind/WAVES spacecraft data. We perform the temporal and spatial association between the radio emission and the listed above activity events during solar cycle (SC) 24, separately for the three sub-categories, m-only, m+DH and DH-only type IIs. A quantitative assessment on the occurrence rates is presented as a function of the strength of the specific SW phenomena: highest rates are obtained with CMEs, SFs, filament eruptions, and SN configurations, whereas a much weaker relationship is found with ICMEs, IP shocks, energetic particles, and GSs. The potential of the obtained rates to be used in empirical or physics-based models for SW forecasting is discussed.

In this paper, we apply Principal Component Analysis (PCA) to experimental data recorded by the KASCADE experiment to reconstruct the mass composition of cosmic rays around the \textit{knee} region. A set of four extensive air shower parameters sensitive to the primary particle mass ($LCm$, $N_{\mu}$, $N_{e}$, and lateral shower $age$) was considered, whose coordinates were transformed into a new orthogonal basis that maximally captures the data variance. Based on the experimental distributions of the first two principal components (PCA0 vs.\ PCA1) and full Monte Carlo simulations of the KASCADE array considering five types of primary particles (p, He, C, Si, and Fe) and three hadronic interaction models (EPOS-LHC, QGSjet-II-04, and SIBYLL~2.3d), we obtained the evolution of the abundance of each primary species as a function of energy, as well as the evolution of the mean logarithmic mass with energy. We found that the reconstruction of the mass composition resulting from this comprehensive analysis significantly reduces dependence on the hadronic interaction model used in the simulation process, even though the initial input parameters are model-dependent. Moreover, the results support the idea that around the \textit{knee} region, the abundance of the light component (protons) decreases, while the heavy component shows a slight increase. The evolution of $\langle \ln A \rangle$ as a function of energy derived from this analysis shows excellent agreement with recent results from the LHAASO--KM2A experiment and aligns very well with the predictions of the data-driven GSF astrophysical model.

We present an investigation of the Extreme-Ultraviolet (EUV) wave linked to the flare that occurred on 28 October 2021, along with the associated coronal loop oscillation and type II radio burst. The EUV wave was observed by multi-viewpoint with Solar Dynamics Observatory and Solar Terrestrial Relations Observatory - A. The associated coronal mass ejection (CME) was observed by Large Angle and Spectrometric Coronagraph (LASCO) as well by COR1 coronagraph. From the multi-view observation, we found that the EUV wave is propagated ahead of the connected CME. The coronal magnetic field measurement was performed by the coronal loop oscillations as well by the associated m-type II radio burst observations. We found the magnetic field strength values computed by both methods are consistence and are in the range of ~ 5 to 10 G.

Solar radio bursts can provide important insights into the underlying physical mechanisms that drive the small and large-scale eruptions on the Sun. Since metric radio observations can give us direct observational access to the inner and middle corona, they are often used as an important tool to monitor and understand the coronal dynamics. While the sizes of the radio sources that can be observed in the solar corona is essential for understanding the nature of density turbulence within the solar corona and its subsequent influence on the angular broadening observed in radio source measurements, the smallest radio sources associated with solar radio bursts have so far been limited by observational techniques and the radio instrument's baselines. We selected three type II bursts that were observed with the LOFAR core and remote stations in the Solar Cycle 24. We estimated the sizes and shapes (ellipticity) of the radio sources between $20-200$ MHz using a two-dimensional Gaussian approximation. Our analysis shows that the smallest radio source size for type II bursts in the solar corona which can be observed in the solar atmosphere at low frequencies is $1.5^\prime \pm 0.5^\prime$ at 150 MHz. However, even though the observations were taken with remote baselines (with a maximum distance of $\sim 85~km$), the effective baselines were much shorter ($\sim 35~km$) likely due to snapshot imaging of the Sun. Our results show that the radio source sizes are less affected by scattering than suggested in previous studies. Our measurements indicate a smaller source sizes at frequencies below 95 MHz compared to previous reports, though some overlap exists with measurements at higher frequencies, using smaller baselines.

If 0.1% of the dark matter in the Coma cluster is constituted by primordial black holes (PBHs) with masses ranging from 10^-19 to 10^-17 solar masses, then the observed GeV {\gamma}-ray emission from the cluster could potentially be attributed to Hawking radiation. The emitted spectrum is inversely proportional to the black hole's mass, meaning lighter PBHs radiate at higher energies, potentially falling within the GeV range. If 0.1% of the Coma cluster's dark matter is PBH in this mass range, a fit to the cluster's GeV emission is obtained. We then investigate the potential for constraining evaporating PBHs through cross-correlations between the Unresolved Gamma-Ray Background and weak gravitational lensing. Utilizing 12 years of Fermi-LAT observations and weak lensing measurements from the Dark Energy Survey Year 3, we assess whether such correlations can reveal a PBH {\gamma}-ray component. While a statistically significant correlation between the UGRB and large-scale structure has been observed, this signal is consistent with emission from clustered astrophysical sources such as blazars. Attributing a measurable fraction of the UGRB to PBH evaporation would require unrealistically large PBH abundances. We also draw attention to a cluster of Coma-like X-ray clusters, designated Draco X, observed at a redshift of z = 0.12. These systems, characterized by their significant X-ray emission from hot, diffuse intracluster gas, represent massive gravitationally bound structures. The existence and properties of such clusters at these redshifts provide crucial cosmological probes, offering insights into the formation and evolution of large-scale structure and the underlying cosmological parameters. Further investigation of Draco X and similar high-redshift clusters would yield additional constraints on large scale structure.

Reproducing the observed activity of comets with thermophysical models remains a primary challenge of cometary science. We use a pebble-based thermophysical model of gas-pressure build-up in the subsurface to reproduce the global emission rates of dust, water, CO$_{2}$, and CO observed by Rosetta at comet 67P/Churyumov-Gerasimenko (hereafter 67P). For sufficiently low diffusivities, the low tensile strength is overcome, leading to the ejection of $\sim$ millimetre- to decimetre-sized dust-particles as well as roughly the correct outgassing rates. All the ejections, and thus the bulk of the outgassing, come from the southern hemisphere during the time that it is strongly illuminated at perihelion. This leads to a 'blow-off' of the dust-crust that otherwise forms: volatiles are much closer to the surface in the south (within the top centimetre) than in the north (10-or-more cm deep), naturally explaining the strong southern water-outgassing expected from 67P's non-gravitational accelerations and torques. We find that low gas-diffusivity, as well as large heat-capacity and steeply decreasing tensile strength with depth or ice-content, are in best agreement with the outgassing data. However, even in these cases, we struggle not to exceed the observed emission rates of dust, CO$_{2}$, and CO. In the south, it is difficult for models to achieve a balance between triggering activity and generating too much of it (with CO$_{2}$ the critical driving-species here); while in the north, it remains challenging to generate activity at all. Strong constraints are placed on the nature of the activity mechanism by the location of dust-ejection and erosion.

The large variety and number of dark energy (DE) theories make it impractical to perform detailed analyses on a case-by-case basis, which has motivated proposals to ``parameterize" theories to reduce the size of theory space. The leading approach to do this is the effective field theory of dark energy (EFTofDE), which can describe general Horndeski-type theories with a small number of observationally accessible time-dependent functions. However, the EFTofDE primarily works for linear perturbations, and extending it to obtain a fully non-linear description of DE theories, which is critical for theories with screening mechanisms, is challenging. In this paper, we present a general method for reconstructing the non-linear DE Lagrangian from the background expansion history and certain linear-perturbation quantities, building upon the EFTofDE framework. Using numerical examples, we demonstrate that this method is applicable to a wide range of single-scalar-field dark energy and modified gravity theories, including quintessence, scalar-tensor theory, $k$-essence, and generalized cubic Galileon with shift symmetry. For each of these theories, we discuss the validity of the method and factors affecting its results. While this method involves solving differential equations, we find that the initial conditions are not important for quintessence, scalar-tensor theory and $k$-essence, while for shift-symmetric cubic Galileon, the generic tracker solution can help transform differential equations into algebraic equations. This offers a useful framework to connect cosmological observations at the background and linear-perturbation levels to the underlying non-linear dynamics of dark energy, and will enable cosmological simulations to analyze and examine DE theories systematically and in much greater detail.

TeV halos are extended very-high-energy gamma-ray sources found around some middle-aged pulsars. The emission spanning several tens of parsecs suggests an efficient confinement of the ultra-relativistic lepton pairs produced by pulsars in their vicinity. The physical mechanism responsible for this suppressed transport has not yet been identified. In some scenarios, pair confinement may be linked to the medium the pulsars are located in. We aim at understanding which type of medium pulsars are probing over their lifetime. We developed a model for the environment probed by moving pulsars, from their birth in core-collapse explosions, where they receive a natal kick, until their entry into the interstellar medium. The model involves: (i) a Monte-Carlo sampling of the properties of the massive-star progenitors of pulsars; (ii) a calculation of the structure of the surrounding medium shaped by these progenitors, for the two cases of isolated stars and star clusters; (iii) a computation of the evolution of supernova remnants in these parent environments. Ultimately, from a distribution of neutron star kick velocities, we assess in which medium pulsars are located as a function of time. We first derive the statistical properties of a fully synthetic Galactic population, and then apply the model to a selection of known pulsars to assess the likely nature of their environment. We show that pulsars escape into the ISM at around 200 kyr, significantly later than the values most commonly used in the literature. The majority of known pulsars with a confirmed TeV halo have high probabilities of still being in their parent environment, which suggests that efficient pair confinement is connected to the region influenced by progenitor stars. In order to test this idea, we provide the probability of still residing in the parent environment for a list of known pulsars.

Accurate identification of meteoroid streams is central to understanding their origins and evolution. However, overlapping clusters and background noise hinder classification, an issue amplified for missions such as ESA's LUMIO that rely on meteor shower observations to infer lunar meteoroid impact parameters. This study evaluates the performance of the Hierarchical Density-Based Spatial Clustering of Applications with Noise (HDBSCAN) algorithm for unsupervised meteoroid stream identification, comparing its outcomes with the established Cameras for All-Sky Meteor Surveillance (CAMS) look-up table method. We analyze the CAMS Meteoroid Orbit Database v3.0 using three feature vectors: LUTAB (CAMS geocentric parameters), ORBIT (heliocentric orbital elements), and GEO (adapted geocentric parameters). HDBSCAN is applied with varying minimum cluster sizes and two cluster selection methods (eom and leaf). To align HDBSCAN clusters with CAMS classifications, the Hungarian algorithm determines the optimal mapping. Clustering performance is assessed via the Silhouette score, Normalized Mutual Information, and F1 score, with Principal Component Analysis further supporting the analysis. With the GEO vector, HDBSCAN confirms 39 meteoroid streams, 21 strongly aligning with CAMS. The ORBIT vector identifies 30 streams, 13 with high matching scores. Less active showers pose identification challenges. The eom method consistently yields superior performance and agreement with CAMS. Although HDBSCAN requires careful selection of the minimum cluster size, it delivers robust, internally consistent clusters and outperforms the look-up table method in statistical coherence. These results underscore HDBSCAN's potential as a mathematically consistent alternative for meteoroid stream identification, although further validation is needed to assess physical validity.

We present a numerical treatment of the annual modulation of relic neutrinos due to the Sun's gravitational influence. Extending our previously developed N-1-body simulation framework from Milky Way scales to solar system dynamics, we model how cosmic neutrino background densities might vary throughout Earth's orbital cycle. We validate our numerical approach against analytical expectations from previous studies that assumed idealized relic neutrino populations. Our results suggest that the prior gravitational history of neutrinos traversing asymmetric dark matter distributions can affect annual modulation patterns. While our simulations reproduce modulation amplitudes similar to analytical predictions for heavier neutrinos, we find that the amplitude can vary considerably depending on the specific morphology of dark matter halos. These findings highlight the importance of incorporating realistic structure formation effects when predicting potential observational relic neutrino signatures.

The orbit of the S2 star around Sagittarius A* provides a unique opportunity to test general relativity and study dynamical processes near a supermassive black hole. Observations have shown that the orbit of S2 is consistent with a Schwarzschild orbit at a 10$\sigma$ confidence level, constraining the amount of extended mass within its orbit to less than 1200 M$_\odot$, under the assumption of a smooth, spherically symmetric mass distribution. In this work we investigate the effects on the S2 orbit of granularity in the mass distribution, assuming it consists of a cluster of equal-mass objects surrounding Sagittarius A*. Using a fast dynamical approach validated by full N-body simulations, we perform a large set of simulations of the motion of S2 with different realizations of the cluster objects distribution. We find that granularity can induce significant deviations from the orbit in case of a smooth potential, causing precession of the orbital plane and a variation of the in-plane precession. Interactions with the cluster objects also induce a sort of "Brownian motion" of Sagittarius A*. Mock data analysis reveals that these effects could produce observable deviations in the trajectory of S2 from a Schwarzschild orbit, especially near apocenter. During the next apocenter passage of S2 in 2026, astrometric residuals in Declination may exceed the astrometric accuracy threshold of GRAVITY of about 30 $\mu as$, as it happens in 35 to 60% of simulations for black holes of 20 to 100 M$_\odot$. This presents a unique opportunity to detect, for the first time, scattering effects on the orbit of S2 caused by stellar-mass black holes, thanks to the remarkable precision achievable with GRAVITY. We also demonstrate that any attempt to constrain the extended mass enclosed within the orbit of S2 must explicitly account for granularity in the stellar-mass black hole population.

Increasing evidence shows that most stars in the Milky Way, including the Sun, are born in star-forming regions containing also high-mass stars, but due to both observational and theoretical challenges, our comprehension of their chemical evolution is far less clear than that of their low-mass counterparts. We present the project "CHemical Evolution of MassIve star-froming COres" (CHEMICO). The project aims at investigating any aspect of the chemical evolution of high-mass star-forming cores by observing representatives of the three main evolutionary categories: high-mass starless cores, high-mass protostellar objects, and ultra-compact HII regions. We carried out an unbiased spectral line survey of the entire bandwidth at 3, 2, and 1.2 mm with the 30m telescope of the Insitut de Radioastronomie millimetrique towards three targets that represent the three evolutionary stages. The number of lines and species detected increases with evolution. In this first work, we derive the temperature structure of the targets through the analysis of the carbon-bearing species C2H, c-C3H, c-C3H2, C4H, CH3CCH, HC3N, CH3CN, and HC5N. The excitation temperature, Tex, increases with evolution in each species, although not in the same way. Hydrocarbons tend to be associated with the smallest Tex values and enhancements with evolution, while cyanides are associated with the highest Tex values and enhancements. In each target, the higher the number of atoms in the molecule, the higher Tex tends to be. The temperature structure evolves from a cold, uniform envelope traced by simple hydrocarbons in the high-mass starless core stage, to a more stratified envelope in the protostellar stage (made by a hot core, a shell with intermediate Tex, and a larger cold envelope), to finally a hot core surrounded only by a cold envelope in the Ultracompact HII stage.

The formation and growth of supermassive black holes (SMBHs) remain a significant unsolved problem in astrophysics, particularly in the low-mass regime, where observations are sparse. The Habitable Worlds Observatory (HWO), with its diffraction-limited imaging and high-resolution UV-optical spectroscopy, presents a unique opportunity to explore the demographics of massive black holes (MBHs) in nearby dwarf galaxies. We propose a program to dynamically detect quiescent MBHs with masses as low as $10^{4.5} \, \rm M_\odot$ in galaxies within $10-30$ Mpc, by resolving stellar velocity dispersions down to 30 $\rm km \, s^{-1}$. This effort will dramatically extend current black hole-host galaxy scaling relations into the dwarf regime, probing the fundamental connection between black hole seeds and their environments. Using a volume-limited sample of $\sim 100$ dwarf galaxies drawn from existing catalogs, HWO will resolve the gravitational sphere of influence of MBHs and enable precision measurements of stellar kinematics via Ca II triplet absorption lines. These observations will test predictions from competing seed formation scenarios -- light seeds from Population III remnants versus heavy seeds from direct collapse -- and clarify whether dwarf galaxies retain imprints of their initial black hole population. The results will offer critical insight into SMBH formation channel and the co-evolution of black holes and galaxies across cosmic time.

The X-ray Integral Field Unit is the X-ray imaging spectrometer on-board one of ESA's next large missions, Athena. Athena is set to investigate the theme of the Hot and Energetic Universe, with a launch planned in the late-2030s. Based on a high sensitivity Transition Edge Sensor (TES) detector array operated at very low temperature (50 mK), X-IFU will provide spatially resolved high resolution spectroscopy of the X-ray sky in the 0.2-12 keV energy band, with an energy resolution goal of 4 eV up to 7 keV [3 eV design goal]. This paper presents the current calibration plan of the X-IFU. It provides the requirements applicable to the X-IFU calibration, describes the overall calibration strategy, and details the procedure and sources needed for the ground calibration of each parameter or characteristics of the X-IFU.

Increasingly large areas in cosmic shear surveys lead to a reduction of statistical errors, necessitating to control systematic errors increasingly better. One of these systematic effects was initially studied by Hartlap et al. in 2011, namely that image overlap with (bright foreground) galaxies may prevent some distant (source) galaxies to remain undetected. Since this overlap is more likely to occur in regions of high foreground density -- which tend to be the regions in which the shear is largest -- this detection bias would cause an underestimation of the estimated shear correlation function. This detection bias adds to the possible systematic of image blending, where nearby pairs or multiplets of images render shear estimates more uncertain and thus may cause a reduction in their statistical weight. Based on simulations with data from the Kilo-Degree Survey, we study the conditions under which images are not detected. We find an approximate analytic expression for the detection probability in terms of the separation and brightness ratio to the neighbouring galaxies. Applying this fitting formula to weak lensing ray tracing through, and the galaxy distribution in the Millennium Simulation, we estimate that the detection bias alone leads to an underestimate of $S_8=\sigma_8\sqrt{\Omega_\mathrm{m}/0.3}$ by almost 2\% and can therefore not be neglected in current and forthcoming cosmic shear surveys.

The $L_{\rm IR}-L_{\rm HCN}$ relation suggests that there is a tight connection between dense gas and star formation. We use data from the Close AGN Reference Survey (CARS) to investigate the dense gas - star formation relation in AGN hosting galaxies, and the use of dense gas as an active galactic nuclei (AGN) diagnostic. Our sample contains five Type-1 (unobscured) AGN that were observed with the Atacama Large Millimeter/submillimeter Array (ALMA) with the aim to detect HCN(4-3), HCO$^+$(4-3) and CS(7-6). We detect the dense gas emission required for this analysis in 3 of the 5 targets. We find that despite the potential impact from the AGN on the line fluxes of these sources, they still follow the $L_{\rm IR}-L_{\rm HCN}$ relation. We then go on to test claims that the HCN/HCO+ and HCN/CS line ratios can be used as a tool to classify AGN in the sub-mm HCN diagram. We produce the classic ionised emission-line ratio diagnostics (the so-called BPT diagrams), using available CARS data from the Multi Unit Spectroscopic Explorer (MUSE). We then compare the BPT classification with the sub-mm classification made using the dense gas tracers. Where it was possible to complete the analysis we find general agreement between optical and sub-mm classified gas excitation mechanisms. This suggests that AGN can contribute to the excitation of both the low density gas in the warm ionised medium and the high density gas in molecular clouds simultaneously, perhaps through X-ray, cosmic ray or shock heating mechanisms.

Determining the physical parameters of pulsating variable stars such as RR Lyrae is essential for understanding their internal structure, pulsation mechanisms, and evolutionary state. In this study, we present a machine learning framework that uses feedforward artificial neural networks (ANNs) to infer stellar parameters-mass ($M$), luminosity (log($L/L_\odot$)), effective temperature (log($T_{\rm eff}$)), and metallicity ($Z$)-directly from Transiting Exoplanet Survey Satellite (TESS) light curves. The network is trained on a synthetic grid of RRab light curves generated from hydrodynamical pulsation models spanning a broad range of physical parameters. We validate the model using synthetic self-inversion tests and demonstrate that the ANN accurately recovers the input parameters with minimal bias. We then apply the trained model to RRab stars observed by the TESS. The observed light curves are phase-folded, corrected for extinction, and passed through the ANN to derive physical parameters. Based on these results, we construct an empirical period-luminosity-metallicity (PLZ) relation: log($L/L_\odot$) = (1.458 $\pm$ 0.028) log($P$/days) + (-0.068 $\pm$ 0.007) [Fe/H] + (2.040 $\pm$ 0.007). This work shows that ANN-based light-curve inversion offers an alternative method for extracting stellar parameters from single-band photometry. The approach can be extended to other classes of pulsators such as Cepheids and Miras.

General Relativity (GR) quantitatively accounts for the anomalous perihelion precession observed in several planets' orbits, with Mercury exhibiting the most significant deviation from Newtonian predictions. In the present work, I carry out an in-depth study comparing methods for calculating perihelion advance of inner planets using two methods, one based on the rotation of the Laplace-Runge-Lentz vector and the other on the evolution of the perihelion longitude. I show that, although the two methods predict almost identical perihelion advances for inner planets according to classical gravitational theory and the general relativity model, they give divergent results for the asteroid Icarus. Comparing numerical calculations of Icarus'perihelion advance with those predicted by Einstein's formula leads to the conclusion that only the method based on the Laplace-Runge-Lenz vector offers a consistent explanation for the behavior of Icarus' perihelion advance. By incorporating the hypothetical planet Vulcan into the Newtonian gravitational model and analyzing its influence on the perihelion advances of the inner planets, I find that carefully selecting Vulcan's mass and orbital parameters yields a delicate balance between its semi-major axis and mass necessary to simultaneously match the observed perihelion advances of the inner planets. Vulcan, orbiting at a semi-major axis of 0.545 AU and possessing about one-third the mass of Mercury, could exert gravitational effects on Earth and Mars comparable to those predicted by GR. This Vulcan would cause an effect on Icarus' perihelion advance approximately nine times greater than the relativistic effect. Thus, precise observational data of Icarus's orbital dynamics offers a definitive means to either rule out this specific Vulcan hypothesis or, conversely, indicate the necessity of further refinement of GR's applications.

The purpose of this paper is to analyze the multi-wavelength and multi-instrument observations of two quiescent filament eruptions as well as the deflection of associated CMEs from the radial direction. The events occurred on 18 October 2017 and 9 May 2021, respectively, in the southern solar hemisphere. Both of them and associated flares were registered by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO) and the Solar Terrestrial Relations Observatory-Ahead (STEREO-A) Observatory in different EUV wavebands. Using data from STEREO A COR1 and COR2 instruments and the Large Angle and Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO), we investigated morphology and kinematics of the eruptions and the latitudinal offset of the related CMEs with respect to the erupting filaments. Our observations provide the evidence that the two filament eruptions were highly non-radial. The observed deviations are attributed to the presence of low-latitude coronal holes.

Context. The HR 8799 system, hosting four giant planets between a warm and cold debris disc, and an extended dusty tail beyond, serves as an ideal laboratory for studying planetary formation and evolution. The debris discs have been observed across various wavelengths, and the planetary properties are well-constrained. Nonetheless, there are still open questions regarding the role of the planets in shaping the debris discs. Aims. We investigate the system's evolution with the aim of understanding how planetary migration shapes its architecture, both in terms of planets and the disc. Methods. We performed $N$-body simulations to model the HR 8799 system. We examined the orbital evolution of the known four super-Jupiter planets through the course of simple, imposed migration in a gaseous disc as they perturb an external, massless planetesimal disc. We also explore the impact of introducing a fifth planet on the dynamical and morphological aspects of the disc. Results. The planets migrate outward as a result of their imposed interactions with the gaseous disc while maintaining their resonant configuration. This outward migration excites the planetesimal disc, producing a transient scattered population. While a four-planet system partially reproduces the observed cavity between the star and the cold debris disc, the inclusion of a fifth low-mass planet appears to be crucial for better reproducing key morphological aspects of the cold debris disc. Conclusions. This model provides a novel explanation for the architecture of HR 8799. Outward planetary migration, combined with mean motion resonant interactions and a fifth, low-mass planet, can effectively replicate the observed planetary architecture and cold debris disc characteristics. Our findings underscore the potential important role of planetary migration in shaping debris discs.

The aim of this work is to construct mock galaxy catalogues that accurately reproduce the redshift evolution of galaxy number density, clustering statistics, and baryonic properties, such as stellar mass for luminous red galaxies (LRGs) and absolute magnitude in the $r$-band for the bright galaxy sample (BGS), based on the first three years of observations from the Dark Energy Spectroscopic Instrument (DESI). To achieve this, we applied the subhalo abundance matching (SHAM) technique to the Uchuu $N$-body simulation, which follows the evolution of 2.1 trillion particles within a volume of $8\,h^{-3}\,\mathrm{Gpc}^{3}$, assuming a Planck base-$\Lambda$CDM cosmology. Using SHAM, we populated Uchuu subhalos with LRGs and BGS-BRIGHT ($r<19.5$) galaxies up to redshift $z=1.1$, assigning stellar masses to LRGs and luminosities to BGS galaxies (up to $M_{\rm r}\leq 20$). Furthermore, we analyzed the clustering dependence on stellar mass and luminosity for each tracer. Our results show that the Uchuu BGS-BRIGHT and LRG mocks accurately reproduce the observed redshift evolution of clustering, with better than 5\% agreement for separations of $1<r<20\,h^{-1}\,\mathrm{Mpc}$ and below 10\% for $0.1<r<1\,h^{-1}\,\mathrm{Mpc}$. For the Uchuu-LRG mock, we successfully captured the stellar mass dependence of clustering, while for the Uchuu-BGS mock, we replicated the clustering for various volume-limited subsamples. We also find good agreement between the data and mocks in the dependence of large-scale bias on luminosity for BGS-BRIGHT galaxies and on stellar mass for LRGs. Altogether, these results equip DESI with robust tools for generating high-fidelity lightcones for the remainder of the survey, thereby enhancing our understanding of the galaxy--halo connection.

We present the first cosmological constraints from a joint analysis of galaxy clustering and galaxy-galaxy lensing using extended SubHalo Abundance Matching (SHAMe). We analyse stellar mass-selected Galaxy And Mass Assembly (GAMA) galaxy clustering and Kilo-Degree Survey (KiDS-1000) galaxy-galaxy lensing and find constraints on $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.793^{+0.025}_{-0.024}$, in agreement with Planck at 1.7$\sigma$, and consistent with previous results. We are able to constrain all 5 SHAMe parameters, which describe the galaxy-subhalo connection. We validate our methodology by first applying it to simulated catalogues, generated from the TNG300 simulation, which mimic the stellar mass selection of our real data. We show that we are able to recover the input cosmology for both our fiducial and all-scale analyses. Our all-scale analysis extends to scales of galaxy-galaxy lensing below $r_\mathrm{p}<1.4\,\mathrm{Mpc}/h$, which we exclude in our fiducial analysis to avoid baryonic effects. When including all scales, we find a value of $S_8$, which is 1.26$\sigma$ higher than our fiducial result (against naive expectations where baryonic feedback should lead to small-scale power suppression), and in agreement with Planck at 0.9$\sigma$. We also find a 21% tighter constraint on $S_8$ and a 29% tighter constraint on $\Omega_\mathrm{m}$ compared to our fiducial analysis. This work shows the power and potential of joint small-scale galaxy clustering and galaxy-galaxy lensing analyses using SHAMe.

The typical amount of molecular hydrogen (${\rm H_2}$) in interstellar ices is not known, but significant freeze-out of ${\rm H_2}$ on dust grains is not expected. However, chemical models ubiquitously predict large amounts of $\rm H_2$ freeze-out in dense cloud conditions, and specialized treatments are needed to control the $\rm H_2$ population on grains. Here we present a numerical desorption model where the effect of weak heating events induced by cosmic rays (CRs) that heat grains to temperatures of a few tens of Kelvin at high frequencies is included, improving upon earlier desorption models that only consider strong heating events (maximum grain temperature close to 100 K) that occur at a low frequency. A temperature of a few tens of Kelvin is high enough to induce efficient desorption of $\rm H_2$, but we find that even the weak heating events do not occur often enough to lead to significant $\rm H_2$ desorption. Taking the weak heating events into account does affect the predicted abundances of other lightly-bound species, but the effect is restricted to low column densities. We make here the canonical assumption that the grains are spherical with a radius of 0.1 $\mu$m. It is conceivable that in the case of a grain size distribution, weak heating events could provide a boost to $\rm H_2$ desorption coming off small grains, which are the most numerous. Further studies are still required to better quantify the role of CRs in the desorption of $\rm H_2$ and other weakly bound species.

We study the protoplanetary disk lifetimes using a large sample of young stellar objects in nearby clusters. To investigate the final phase of disk dissipation, we selected 32 clusters, located within 500 pc and aged between 1 and 100 Myr, with membership determined using Gaia data. The Age and mass information of the sources are obtained through spectral energy distribution (SED) analysis and using evolutionary models of various ages. Using the IR data from 2MASS and WISE catalogues, we employ three methods to identify disks across the different wavelength regimes (1.1- 22 $\mu$m). We find that disk fraction consistently decreases as stellar systems age, a trend observed across all wavelengths included in this study. However, there is an increase in the time scale of disk decay as wavelength increases, with characteristic timescales of $\tau_{\text{short}}$ = 1.6 $\pm$ 0.1 Myr for shorter wavelengths (1.6-4.6 $\mu$m) versus $\tau_{\text{W3}}$ = 4.4 $\pm$ 0.3 Myr for 12 $\mu$m. This supports the idea that outer disk regions evolve more slowly. Notably, we detect infrared excesses at 12 $\mu$m and 22 $\mu$m in relatively older systems ($>$10 Myr), with some disks with estimated ages up to $\sim$ 100 Myr. Among these, we identify a population of full disks that persist beyond the typical dissipation timescale. We also observe that the median mass of disk-hosting stars decreases from 0.62 $M_\odot$ to 0.27 $M_\odot$ in clusters younger and older than 40 Myr, respectively, indicating slower disk dissipation around lower-mass stars. We identify 33 transitional disk candidates using various color-color diagrams. Using LAMOST DR8 optical spectra and H-alpha equivalent widths, we identify possible accretors and estimate their mass accretion rates, finding most are younger than 10 Myr.

In the standard picture of cosmology, the galaxies reside in dark matter (DM) halos. DM halos are distributed in the cosmic web in different environments. The luminosity of the galaxies in different environments can be used as a probe to assess a cosmological model. This study focuses on the properties of galaxies in void regions, where halos typically do not experience extreme conditions. By examining the galaxy luminosity function, we aim to understand the dependence of galaxy properties on their environment and redshift so that later, we could use this as a tool to evaluate cosmological models. We employ the Excursion Set Theory to incorporate parameters related to the number density of DM halos into the luminosity function. Using the Galaxy and Mass Assembly (GAMA) survey and 2dFGRS datasets, we fit our theoretical models to observational data, examining the environmental and redshift dependence of the galaxy luminosity function. Our results indicate that we model the galaxy luminosity function in voids effectively by considering the linear density contrast of the environment and the growth function $D(z)$ for redshift dependence. This study provides a model for the environmental dependence of galaxy luminosity function that offers an improvement in the $\chi ^2$ parameter compared to the previously proposed model in \cite{mcnaught2014galaxy}. Both BIC and AIC tests support the superiority of this model for the void region.

Periodic variations in the Sun and Moon's gravitational pull cause small changes in Earth's rotational axis direction called nutation. Nutation components in the retrograde quasi-diurnal frequency band measured in the terrestrial reference frame are amplified by resonance with the Free Core Nutation (FCN), a rotational mode of Earth's fluid core. Dissipative processes at the core-mantle boundary (CMB) dampen this resonance, contributing to the observed phase lag between tidal forcing and Earth's rotational response. This phase lag is commonly attributed to electromagnetic (EM) coupling between the core and the electrically conducting lower mantle. However, estimates of mantle conductivity and radial magnetic field strength at the CMB suggest these effects are insufficient. We show that the missing dissipation arises naturally from the excitation of internal waves in the fluid core by topographic features at the CMB. Adapting a theoretical framework originally developed for tidal flow over oceanic topography, we compute the form drag and associated power flux induced by CMB topography. Our results are consistent with a CMB topography characterized by a root mean square amplitude of a few kilometers. The model favors weak stratification at the top of the core, though stronger stratification remains compatible with increased topographic amplitude. This mechanism provides independent constraints on CMB topography and stratification, complementing seismological and magnetic observations. Its generality offers a new framework for probing deep-interior dynamics across terrestrial planets.

Many previous studies have imposed stringent constraints on the particle mass of fuzzy dark matter (FDM) by analyzing observations of Galactic satellite galaxies, which show no significant evidence of the heating effect predicted by FDM. However, these analyses have generally neglected the tidal influence of the Milky Way, which can substantially suppress the FDM-induced heating in satellites. This oversight arises from computational challenges of accurately capturing the tidal effects in FDM simulations. In this study, we present a novel simulation framework that, for the first time, enables the simulation of an FDM-stellar system within an observationally motivated gravitational potential of the Milky Way. This framework incorporates the diverse Galactic components, including the gravitational influence of the Large Magellanic Cloud. Using the Fornax dwarf galaxy as a case study, we demonstrate that tidal effects significantly alleviate the tension between observational data and the predicted heating effect for an FDM particle mass of $m_a\sim 10^{-22}$ eV.

Hydrogen recombination lines (HRLs) are valuable diagnostics of the physical conditions in ionized regions around high-mass stars. Understanding their broadening mechanisms and intensity trends can provide insights into the densities, temperatures, and kinematics of HII regions. We investigate the properties of ionized gas around massive protostars by analyzing hydrogen recombination lines (H-alpha and H-beta) in the Q-band. Observations were conducted using the Yebes 40m radio telescope in the Q-band (30.5~50 GHz) toward six high-mass protostars selected from the SOMA Survey (G45.12+0.13, G45.47+0.05, G28.20-0.05, G35.20-0.74, G19.08-0.29, and G31.28+0.06). The line profiles were analyzed to assess broadening mechanisms, from which electron densities and temperatures were derived. We compared our results with Q-band data from the TianMa 65m Radio Telescope (TMRT) and ALMA Band 1 Science Verification observations of Orion KL. A total of eight H-alpha (n = 51 to 58) and ten H-beta (n = 64 to 73) lines were detected toward G45.12+0.13, G45.47+0.05, and G28.20-0.05, with non-detections in the other sources. Electron densities of \sim \mathrm{1\text{--}5 \times 10^{6}~cm^{-3}} and temperatures of 8000-10000 K were derived. Orion KL shows one order of magnitude lower electron density, but a similar temperature. Notably, G45.12 and G28.20 show increasing intensity with frequency for both H-alpha and H-beta, in contrast to the decreasing trend in Orion KL. The observed line widths indicate contributions from both thermal and dynamical broadening, suggesting high-temperature ionized gas affected by turbulence, outflows, rotation, or stellar winds. Pressure broadening may also play a minor role. The contrasting intensity trends likely reflect differences in local physical conditions or radiative transfer effects, warranting further study through higher-resolution observations and modeling.

Context: Weak gravitational lensing is a key cosmological probe for current and future large-scale surveys. While power spectra are commonly used for analyses, they fail to capture non-Gaussian information from nonlinear structure formation, necessitating higher-order statistics and methods for efficient map generation. Aims: To develop an emulator that generates accurate convergence maps directly from an input power spectrum and wavelet l1-norm without relying on computationally intensive simulations. Methods: We use either numerical or theoretical predictions to construct convergence maps by iteratively adjusting wavelet coefficients to match target marginal distributions and their inter-scale dependencies, incorporating higher-order statistical information. Results: The resulting kappa maps accurately reproduce the input power spectrum and exhibit higher-order statistical properties consistent with the input predictions, providing an efficient tool for weak lensing analyses.

IZw18 is one of the lowest-metallicity star-forming galaxies known at z$\sim$0, considered a unique local analogue of the first galaxies. The origin of its hard ionizing continuum, expected to be a common feature in the early Universe and traced by He\textsc{ii} emission lines, remains intensely debated and challenging to explain. Here we combine optical (GTC/MEGARA) and mid-infrared (JWST/MIRI) integral field spectroscopic observations for IZw18 to shed new light on the high-ionization phenomenon. This letter reports the first detection of the high-ionization [Ne\textsc{v}]14.32 $\mu$m line in IZw18. Its emission is spatially extended and coincident with the He\textsc{ii} peak, revealing the presence of highly energetic ionizing sources that surpass mechanisms previously proposed on the basis of He\textsc{ii} alone. Our kinematic studies highlight that the He\textsc{ii}$\lambda$4686-emitting gas displays higher velocity dispersions and a different velocity pattern compared to the H$\beta$ emission, suggesting the presence of energetic processes such as shocks or stellar-driven feedback. Additionally, integrated spectra show asymmetric blueshifted profiles in the He\textsc{ii}$\lambda$4686 line, possibly indicating \textbf{early-stage} stellar-driven outflows potentially facilitating future ionizing photon leakage. Our spatial analysis also reveals differences in structure between the emission of H$\beta$ and He\textsc{ii}$\lambda$4686, with the He\textsc{ii}$\lambda$4686 peak offset by a projected distance of 140 pc from the peak H$\beta$ emission. This indicates distinct locations for the most extreme ionizing sources compared to moderate ionizing sources. Our findings underscore the complex interplay of physical processes in extremely metal-poor environments with \textbf{high-ionized} gas, offering new insights into the conditions prevailing in the early galaxies.

In this paper we revisit our spectropolarimetric and velocimetric analysis of the young M dwarf AU Mic based on data collected with SPIRou at the Canada-France-Hawaii telescope, over a monitoring period of 2041 d from 2019 to 2024. The longitudinal magnetic field, the small-scale magnetic field, and the differential temperature of AU Mic, derived from the unpolarized and circularly-polarized spectra, were clearly modulated with the stellar rotation period, with a pattern that evolved over time. The magnetic modeling with Zeeman-Doppler imaging provides a consistent description of the global field of AU Mic that agrees not only with the Least-Squares Deconvolved profiles of the circularly-polarized and unpolarized spectral lines, but also with the small-scale field measurements derived from the broadening of spectral lines, for each of the 11 subsets of the full data. We find that the large-scale field was mostly poloidal, with a dominant dipole component slightly tilted to the rotation axis which decreased from 1.4 to 1.1 kG before increasing at the end of the campaign. The average small-scale field followed a similar trend, decreasing from 2.8 to 2.6 kG then rising. The long-term magnetic evolution we report for AU Mic suggests that, if cyclic, the cycle period is significantly longer than 6 yr. From velocimetric data, we derived improved mass estimates for the two transiting planets, respectively equal to M_b = 6.3+2.5-1.8 M_earth and M_c = 11.6+3.3-2.7 M_earth, yielding very contrasting densities of 0.32+0.13-0.10 and 2.9+1.1-0.8 g/cm3, and a new 90% confidence upper limit of 4.9 M_earth for candidate planet d (period 12.7 d) suspected to induce the transit-timing variations of b and c. We also confirm our claim regarding candidate planet e orbiting with a period of 33.11+-0.06 d, albeit with a smaller mass of M_e = 21.1+5.4-4.3 M_earth.

We report in this paper circularly and linearly polarized observations of the young active M dwarfs AU Mic and EV Lac with the near-infrared SPIRou spectropolarimeter at the Canada-France-Hawaii Telescope, collected from August to October 2023 over a few rotation cycles of both stars. Applying Least-Squares Deconvolution (LSD) to our spectra, we detected Zeeman signatures in circular (Stokes V) and linear (Stokes QU) polarization, and Zeeman broadening in unpolarized (Stokes I) LSD profiles, all exhibiting clear rotational modulation. Using the stellar surface tomographic technique of Zeeman-Doppler imaging on our sets of observations, along with a simple parametric description of how the small-scale and large-scale fields relate to each other, we recovered the magnetic topologies of AU Mic and EV Lac successively from LSD Stokes V, Stokes IV and Stokes IVQU profiles, to investigate how the reconstructed maps evolve as we provide more information, and ultimately infer reliable magnetic maps of both stars. We find that AU Mic hosts a fairly simple and mostly poloidal large-scale field aligned with the rotation axis within about 10deg, whereas that of EV Lac is more complex, stronger and less axisymmetric. Both stars feature intense small-scale fields, of about 4 kG for AU Mic and 6 kG for EV Lac when averaged over the whole stellar surface. Stokes QU Zeeman signatures allow one to reconstruct stellar magnetic fields more reliably, and are especially useful for stars with more complex fields and low vsini's like EV Lac.

The origin of magnetic fields observed on both astrophysical and cosmological scales is a compelling problem that has the potential to shed light on the early Universe. We analytically investigate inflationary magnetogenesis in scenarios where a brief departure from slow-roll inflation - akin to mechanisms proposed for primordial black hole formation - leads to enhanced magnetic field generation with a growing power spectrum. Focusing on the Ratra model, we derive an analytic bound on the growth of the magnetic field power spectrum in this context, showing that the spectral index can reach $d \ln {\cal P}_B / d \ln k = 4.75$ during the growth phase. This growth enables amplification from CMB-safe large-scale amplitudes to values of astrophysical relevance. We further compute the stochastic gravitational wave background sourced by the resulting magnetic fields, incorporating their rich spectral features. Under suitable conditions, the induced signal exhibits a characteristic frequency dependence and amplitude within reach of future gravitational wave observatories, providing a distinctive signature of this mechanism and a specific class of templates for upcoming gravitational wave searches.

Advances in optical astrometry allow us to infer the non-radial kinematic structure of the Universe directly from observations. Here I use a supervised machine learning neural network method to predict 1.57 million redshifts based on several photometric and metadata classifier parameters from the unWISE mid-infrared database and from Gaia. These estimates are used to divide the sample into three redshift bins: 1-2, 2-3, and $>3$. For each subset, all available Gaia proper motions are used in a global vector spherical harmonic solution to degree 3 (30 fitting vector functions). I find significant differences in a few fitted proper motion patterns at different redshifts. The largest signals are seen in the comparison of the vector spherical harmonic fits for the 1-2 and 2-3 redshift bins. The significant harmonics include a rigid spin, a dipole glide from the north Galactic pole to the south and an additional quadrupole distortion. Validation tests with filtered subsamples indicate that the detected effect can be caused by hidden systematic errors in astrometry. The results are verified by using an independent source of redshifts and computing the observer's Galactocentric acceleration. This study offers a new observational test of alternative cosmological models.

The obliquity between the stellar spin axis and the planetary orbit, detected via the Rossiter-McLaughlin (RM) effect, is a tracer of the formation history of planetary systems. While obliquity measurements have been extensively applied to hot Jupiters and short-period planets, they remain rare for cold and long-period planets due to observational challenges, particularly their long transit durations. We report the detection of the RM effect for the 19-hour-long transit of HIP 41378 f, a temperate giant planet on a 542-day orbit, observed through a worldwide spectroscopic campaign. We measure a slight projected obliquity of 21 $\pm$ 8 degrees and a significant 3D spin-orbit angle of 52 $\pm$ 6 degrees, based on the measurement of the stellar rotation period. HIP 41378 f is part of a 5-transiting planetary system with planets close to mean motion resonances. The observed misalignment likely reflects a primordial tilt of the stellar spin axis relative to the protoplanetary disk, rather than dynamical interactions. HIP 41378 f is the first non-eccentric long-period (P>100 days) planet observed with the RM effect, opening new constraints on planetary formation theories. This observation should motivate the exploration of planetary obliquities across a longer range of orbital distances through international collaboration.

Astronomical objects in our universe that are too faint to be directly detectable exist and are important - an obvious example being dark matter. The same can also apply to very faint baryonic objects, such as low luminosity dwarf galaxies and gravitationally compact objects (e.g., rogue planets, white dwarfs, neutron stars, black holes, dark sirens). While they are very difficult to observe directly, they have locations that are highly important when studying astrophysical phenomena. Here, we use a machine learning algorithm known as symbolic regression to model the probability of a dark object's existence as a function of their separation distances to their closest two ``bright" (directly observable) neighbors, and the distances of these bright objects to each other. An advantage of this algorithm is that it is interpretable by humans and can be used to make reproducible predictions. Galaxies with masses above $10^9 M_{\odot}$ and halos above $10^{12} M_{\odot} $ are the objects that we separate into ``bright" and ``dark" to be used in our analysis. We find that it is possible to predict the density of dark objects using an analytic expression that depends on their distances to their closest bright neighbors in Illustris-TNG galaxy formation simulations, which is significantly better than the (linear) scale-dependent biasing prediction for $k \sim 1.0~ h$Mpc$^{-1}$ (and potentially beyond, if allowed by the resolution). This could potentially open the avenue for finding dark objects based on their vicinity to directly observable bright sources and make future surveys more targeted and efficient.

We present \textbf{GeQuLEP} (Germanium-based Quantum Sensors for Low-Energy Physics), a conceptual design for an advanced quantum sensing platform integrating high-purity germanium (Ge) crystals with engineered phononic crystal cavities. At cryogenic temperatures, these cavities naturally host dipole-bound states, effectively forming quantum dots coupled to radio-frequency quantum point contact (RF-QPC) readout systems. This innovative coupling approach promises ultra-sensitive phonon-mediated charge detection through phonon-induced charge displacement. GeQuLEP is specifically designed to achieve exceptionally low detection thresholds, theoretically enabling single primary phonon sensitivity with anticipated energy depositions as low as \textbf{0.00745~eV}. This unprecedented sensitivity, if realized experimentally, would provide unique access to searches for low-mass dark matter down to the keV/$c^2$ mass range via nuclear and electronic recoils. Additionally, GeQuLEP aims to facilitate the real-time detection of solar \textit{pp} neutrinos through coherent elastic neutrino--nucleus scattering (CE$\nu$NS). By combining phonon-based quantum transduction with quantum-classical hybrid readout schemes, the GeQuLEP architecture represents a scalable, contact-free phonon spectroscopy design that could significantly advance the capabilities of ultra-low-energy rare-event detection at the quantum limit.

Galaxy clusters, the pinnacle of structure formation in our universe, are a powerful cosmological probe. Several approaches have been proposed to express cluster number counts, but all these methods rely on empirical explicit scaling relations that link observed properties to the total cluster mass. These scaling relations are over-parametrised, inducing some degeneracy with cosmology. Moreover, they do not provide a direct handle on the numerous non-gravitational phenomena that affect the physics of the intra-cluster medium. We present a proof-of-concept to model cluster number counts, that bypasses the explicit use of scaling relations. We rather implement the effect of several astrophysical processes to describe the cluster properties. We then evaluate the performances of this modelling for the cosmological inference. We developed an accelerated machine learning baryonic field-emulator, built upon the Lagrangian Deep Learning method and trained on the CAMELS simulations. We then created a pipeline that simulates cluster counts in terms of XMM observable quantities. We finally compare the performances of our model, with that involving scaling relations, for the purpose of cosmological inference based on simulations. Our model correctly reproduces the cluster population from the calibration simulations at the fiducial parameter values, and allows us to constrain feedback mechanisms. The cosmological-inference analyses indicate that our simulation-based model is less degenerate than the approach using scaling relations. This novel approach to model observed cluster number counts from simulations opens interesting perspectives for cluster cosmology. It has the potential to overcome the limitations of the standard approach, provided that the resolution and the volume of the simulations will allow a most realistic implementation of the complex phenomena driving cluster evolution.

Kinetic misalignment, one of the most compelling scenarios for the non-thermal generation of axion dark matter, is generally accompanied by axion fragmentation, a process in which the energy of the axion condensate is transferred to its perturbations. The dynamics of fragmentation, at least in the context of dark matter production, have so far been studied semi-analytically using perturbation theory. In this work, we present the first classical lattice simulation of kinetic axion fragmentation in the context of dark matter production, focusing on parameters relevant to the QCD axion. Our findings indicate that the non-perturbative dynamics captured by the lattice lead to a significantly broader spectrum of axion fluctuations, with a sustained transfer of energy to mildly relativistic modes and with smaller occupation numbers compared to the linear approximation. As a consequence, the final dark matter abundance is typically O(1) lower than in the linear approximation, which is itself O(1) lower than the zero-mode-only prediction. This broadening and suppression of the spectrum could have a significant impact on axion mini-halo formation, one of the main experimental handles on kinetic fragmentation.

Cosmochemical studies have proposed that Earth accreted roughly 5-10% of its mass from carbonaceous (CC) material, with a large fraction delivered late via its final impactor, Theia (the Moon-forming impactor). Here, we evaluate this idea using dynamical simulations of terrestrial planet formation, starting from a standard setup with a population of planetary embryos and planetesimals laid out in a ring centered between Venus and Earth's orbits, and also including a population of CC planetesimals and planetary embryos scattered inward by Jupiter. We find that this scenario can match a large number of constraints, including i) the terrestrial planets' masses and orbits; ii) the CC mass fraction of Earth; iii) the much lower CC mass fraction of Mars, as long as Mars only accreted CC planetesimals (but no CC embryos); iv) the timing of the last giant (Moon-forming) impact; and v) a late accretion phase dominated by non-carbonaceous (NC) bodies. For this scenario to work, the total mass in scattered CC objects must have been ~ 0.2 - 0.3 M$_{\oplus}$ , with an embryo-to-planetesimal mass ratio of at least 8, and CC embryos in the ~ 0.01 - 0.05 M$_{\oplus}$ mass range. In that case, our simulations show there are roughly 50-50 odds of Earth's last giant impactor (Theia) having been a carbonaceous object - either a pure CC embryo or an NC embryo that previously accreted a CC embryo. Our simulations thus provide dynamical validation of cosmochemical studies.

Assuming a minimal $\Lambda$CDM cosmology with three massive neutrinos, the joint analysis of Planck cosmic microwave background data, DESI baryon acoustic oscillations, and distance moduli measurements of Type Ia supernovae from the Pantheon+ sample sets an upper bound on the total neutrino mass, $\sum m_\nu \lesssim 0.06$-$0.07$ eV, that lies barely above the lower limit from oscillation experiments. These constraints are mainly driven by mild differences in the inferred values of the matter density parameter across different probes that can be alleviated by introducing additional background-level degrees of freedom (e.g., by dynamical dark energy models). However, in this work we explore an alternative possibility. Since both $\Omega_\mathrm{m}$ and massive neutrinos critically influence the growth of cosmic structures, we test whether the neutrino mass tension may originate from the way matter clusters, rather than from a breakdown of the $\Lambda$CDM expansion history. To this end, we introduce the growth index $\gamma$, which characterizes the rate at which matter perturbations grow. Deviations from the standard $\Lambda$CDM value ($\gamma \simeq 0.55$) can capture a broad class of models, including non-minimal dark sector physics and modified gravity. We show that allowing $\gamma$ to vary significantly relaxes the neutrino mass bounds to $\sum m_\nu \lesssim 0.13$-$0.2$ eV, removing any tension with terrestrial constraints without altering the inferred value of $\Omega_\mathrm{m}$. However, this comes at the cost of departing from standard growth predictions: to have $\sum m_\nu \gtrsim 0.06$ eV one needs $\gamma > 0.55$, and we find a consistent preference for $\gamma > 0.55$ at the level of $\sim 2\sigma$. This preference increases to $\sim 2.5$-$3\sigma$ when a physically motivated prior $\sum m_\nu \ge 0.06$ eV from oscillation experiments is imposed.

We present an updated measurement of the diffuse astrophysical neutrino flux using Baikal-GVD cascade data collected between April 2018 to March 2024. In this period, the detector grew from 15% to 55% of its baseline cubic kilometer configuration. The diffuse astrophysical neutrino flux is detected with a statistical significance of 5.1 $\sigma$. Assuming a single power law model of the astrophysical neutrino flux with identical contribution from each neutrino flavor, the following best-fit parameter values are found: the spectral index $\gamma_{astro}$ = 2.64$^{+0.09}_{-0.11}$ and the flux normalization $\phi_{astro}$ = 4.42$^{+2.31}_{-1.29}\times10^{-18} \text{GeV}^{-1}\text{cm}^{-2}\text{s}^{-1}\text{sr}^{-1}$ per one flavor at 100 TeV. These results are broadly consistent with IceCube measurements.

Stellar-gas kinematic misalignments are a transient phenomenon observed in $\sim11\%$ of the local galaxy population. According to current models, misaligned gas is expected to lose angular momentum and relax into the galactic plane on timescales of $\sim0.1$ Gyr, driving gas toward the central regions of the galaxy. Recent observational studies have found a higher incidence of active galactic nuclei in misaligned galaxies. We use the EAGLE simulation to explore the connection between stellar-gas misalignments and enhanced central black hole (BH) activity between $0<z<1$. We use a sample of $\sim5600$ galaxies with a stellar mass of $M_{*}\geqslant \mathrm{10^{9.5}}$ M$_\odot$ that feature long-lived stellar-gas alignment, counter-rotation, and unstable misalignments (non-coplanarity). Over time windows of $0.5$ Gyr, we find that galaxies experiencing an unstable misalignment have systematically enhanced BH growth during relaxation. Galaxies with long-term counter-rotation show little difference in BH growth compared to aligned galaxies. We suggest that this enhanced BH growth is driven by loss of angular momentum in unstable misaligned gas discs which is able to drive gas inward toward the vicinity of the BH. At $z\approx0.1$, we find a greater incidence of overmassive BHs in galaxies that have spent a greater fraction of time with unstable stellar-gas kinematic misalignments over the preceding $\approx2$ Gyr compared to control samples of aligned galaxies. In agreement with observations, we conclude that BH activity is enhanced in misaligned systems in EAGLE and suggest that the presence of overmassive BHs may be indicative of a past stellar-gas kinematic misalignment.

Using high spatiotemporal resolution, multi-wavelength observations from the New Vacuum Solar Telescope (NVST) and the Solar Dynamics Observatory (SDO), we present a detailed analysis of a small-scale chromospheric jet driven by plasmoid-mediated magnetic reconnection. Our results reveal that the entire process is governed by the dynamic evolution of photospheric magnetic footpoints, which proceeds in two distinct stages. An initial separating motion of the footpoints corresponds to a mild reconnection phase, characterized by a short current sheet and the eruption of a cool H$\alpha$ jet. Subsequently, a converging motion of the footpoints triggers an intense reconnection phase. During this intense stage, the current sheet rapidly elongates, and the resulting decrease in its aspect ratio initiates a tearing-mode instability, forming a plasmoid. The appearance of this plasmoid mediates the onset of fast magnetic reconnection, which produces a hot EUV jet and is concurrent with significant magnetic flux cancellation. We interpret this cancellation as the submergence of newly formed, post-reconnection loops. Furthermore, we identify a distinct, high-temperature plasma blob in the jet spire, significantly hotter than the surrounding jet plasma. We attribute this feature to a secondary heating process, likely caused by reconnection between the upward-propagating plasmoid and the overlying magnetic cusp structure. These observations provide a comprehensive, observationally driven picture (from the initial photospheric triggers to the multi-stage, plasmoid-mediated reconnection) that forms chromospheric jets, highlighting the critical role of footpoint motions in solar atmospheric dynamics.

We identify a new class of time periodic attractor solutions in scalar field cosmology, which we term Cosmological Frequency Combs (CFC). These solutions arise in exponential quintessence models with a phantom matter background and exhibit coherent phase-locked oscillations in the scalar field's normalized variables. We demonstrate that such dynamics induce modulations in observables like the Hubble parameter and growth rate, offering a dynamical mechanism to even address the H0 tension exactly. Our results uncover a previously unexplored phase of cosmic acceleration, linking the concept of frequency combs to large scale cosmological evolution.

In recent years, large language models (LLMs) have transformed natural language understanding through vast datasets and large-scale parameterization. Inspired by this success, we present SpecCLIP, a foundation model framework that extends LLM-inspired methodologies to stellar spectral analysis. Stellar spectra, akin to structured language, encode rich physical and chemical information about stars. By training foundation models on large-scale spectral datasets, our goal is to learn robust and informative embeddings that support diverse downstream applications. As a proof of concept, SpecCLIP involves pre-training on two spectral types--LAMOST low-resolution and Gaia XP--followed by contrastive alignment using the CLIP (Contrastive Language-Image Pre-training) framework, adapted to associate spectra from different instruments. This alignment is complemented by auxiliary decoders that preserve spectrum-specific information and enable translation (prediction) between spectral types, with the former achieved by maximizing mutual information between embeddings and input spectra. The result is a cross-spectrum framework enabling intrinsic calibration and flexible applications across instruments. We demonstrate that fine-tuning these models on moderate-sized labeled datasets improves adaptability to tasks such as stellar-parameter estimation and chemical-abundance determination. SpecCLIP also enhances the accuracy and precision of parameter estimates benchmarked against external survey data. Additionally, its similarity search and cross-spectrum prediction capabilities offer potential for anomaly detection. Our results suggest that contrastively trained foundation models enriched with spectrum-aware decoders can advance precision stellar spectroscopy.

We present multi-wavelength imaging and analysis of a low surface brightness (LSB) dwarf galaxy in the Extended Chandra Deep Field South (ECDFS), SMDG0333094-280938, with particular emphasis on data from the Euclid space telescope and from the Vera C.\ Rubin Observatory. The galaxy is clumpy and blue, and appears to host globular clusters (GCs), suggesting a distance of ~50-60 Mpc which would make the dwarf an ultra-diffuse galaxy (UDG). We carry out spectral energy distribution (SED) fitting from the far-ultraviolet to the near-infrared, in order to estimate the galaxy age and metallicity. We infer a recent peak of star formation that may have led to the formation of the UDG through feedback-driven expansion. This early analysis illustrates how Euclid and Rubin are poised to identify and characterize many thousands of UDGs and other LSB galaxies in the near future, including their GCs and stellar populations.

We present the first optical-UV spectral systematic analysis of 30 Type 1 AGN selected in the FIR and X-ray in the Lockman-SpReSO Survey. The sample of faint objects (m_B = 19.6-21.8) covers a large redshift range of 0.33 > z > 4.97 with high S/N (~21 on average). A detailed spectral analysis based on the Quasar Main Sequence phenomenology prescription was applied to deblend the principal optical-UV emitting regions. Our sample spans a bolometric luminosity range of 44.85 < log Lbol < 47.87, absolute B-magnitude -20.46 > M_B > -26.14, BH mass of 7.59 < log MBH < 9.80, and Eddington ratio -1.70 < log REdd < 0.56. The analysis shows that 18 high-z objects correspond to Population B, whereas three low-z fall in Populations A2, B1, and B1+. The remaining eight are candidates to be Pop. B and one Pop. A object. None of them are extreme accretors. We looked for tendencies in our sample and compared them with other samples with different selection criteria. Evidence for winds was explored using CIV1549 line half-height centroid cmed finding wind velocities between 941 and -1587 kms-1. This result is consistent with samples with similar ranges of z and M_B. The Baldwin effect showed a slope of -0.23 pm 0.03 dex consistent with previous studies. Spectra from twelve objects in our sample were found in the Sloan Digital Sky Survey Data Release 17 database. We applied the same methodology to compare them to our spectra, finding no evidence of variability.

We show that enhancement of the axion relic abundance compared to the standard misalignment contribution generically leads to the production of nonzero momentum axion modes, resulting in warm dark matter behavior and enhanced isocurvature perturbations. It leads to universal constraints on the axion parameter space that are independent of detailed model assumptions and cosmological history. For models enhancing relic abundance with gradient axion modes, observations of the Lyman-$\alpha$ forest impose a lower bound on the axion decay constant, $f_a \gtrsim 10^{15} {\rm GeV}\,(10^{-18}{\rm eV}/m_a)$, from the free-streaming effect. For models relying on the delay of coherent axion oscillations, we obtain a slightly weaker bound, $f_a \gtrsim 10^{14} {\rm GeV}\,(10^{-18}{\rm eV}/m_a)$. We make relatively conservative choices to establish these universal bounds but also provide scaling parameters that can be calibrated for stronger constraints in concrete models and updated as observations improve.

The varied and dynamic evolutionary histories of galaxies give rise to the stunning diversity in their properties that we observe in the present-day universe. HST, and now JWST, have pioneered the study of resolved individual stars in the Milky Way and other members of the Local Group, uncovering the drivers of their morphological, star formation, and chemical evolution. HWO will constitute a paradigm shift: introducing the ability to panchromatically resolve the main bodies of every galaxy in the Local Volume into their constituent stars. In this science case, we summarize the breakthrough progress that HWO will advance in the field of galaxy evolution through resolved stellar populations. HWO will transform our understanding of galaxies in three distance regimes: (1) in the nearest galaxies ($\sim$5 Mpc), where it will resolve stars below the oldest Main Sequence Turnoff, enabling precision stellar astrophysics and star formation history (SFH) inferences to the earliest cosmic times; (2) in the greater Local Volume ($\sim$20 Mpc), where it will resolve stars below the Red Clump, providing access to accurate SFHs for hundreds of galaxies, spanning the entire Hubble Sequence; and (3) out to cosmological volumes ($\sim$50+ Mpc), providing access to the luminous stellar populations in thousands of galaxies, enabling unprecedented views of their morphology, stellar abundances, and dust content. The principal technological requirement advanced by this science case is a camera with a resolution of $\leqslant$0.015'' that is diffraction-limited, and Nyquist-sampled (0.01'' per pixel), to at least 550 nm $-$ comparable to the High Definition Imager from the LUVOIR concept.

We demonstrate the formation of quasi-stable localized scalar configurations in spontaneously symmetry breaking U(1) model by 3+1-dimensional classical lattice simulations. Such configurations are called PQ-balls, as the primary motivation of this kind of configuration is Peccei-Quinn theory under the kinetic misalignment mechanism. Our numerical simulations demonstrate that they can form if the PQ charge is generated through the coherent rotation of a complex scalar field in the complex plane, via dynamics analogous to the Affleck-Dine mechanism. These configurations subsequently decay due to the U(1)-breaking effect induced by spontaneous symmetry breaking. We also demonstrate the formation and decay of oscillons in a similar setup in a real scalar field theory.

We extend the framework of spontaneous baryogenesis by investigating the generation of baryon asymmetry when the inflaton, $\theta$, is minimally coupled with a complex spectator scalar field $\phi$, as $\theta^2|\phi|^2$. To do so, we also consider $\phi$ non-minimally coupled with the Ricci scalar curvature $R$ through a Yukawa-like interaction. We do not consider further interactions of the spectator field with the fermions of the Standard Model, considering it \emph{de facto} as a dark scalar field. In evaluating the violation of the baryon-number conservation during the reheating epoch, in a perfectly homogeneous and isotropic universe, we follow a semiclassical approach, where $\theta$, $\phi$ and gravity are considered as classical fields, whereas the fermions are quantized. We solve the equations of motion for the inflaton and spectator fields, respectively at first and zero-order in perturbation theory, neglecting at first stage the expansion of the universe. Afterwards, we quantify how the spectator field modifies the inflationary dynamics and thus find the baryon asymmetry produced via the inflaton decays into fermion-antifermion pairs by computing the corresponding decay amplitudes. We therefore obtain small first order correction to standard spontaneous baryogenesis and finally discuss the mass-mixing between fermions. Accordingly, the effects of considering the universe expansion are accounted, showing when the coupling between $\phi$ and $R$ becomes noticeable in altering the overall baryon asymmetry.

We present a comprehensive and gauge-invariant study of the neutrino dipole portal at the energy frontier. Assuming negligible active-sterile mixing, we analyze sterile neutrino production via dimension-6 dipole operators coupling to the electroweak field strengths. The analysis incorporates, for the first time, both single- and double-gauge-boson effective interactions. We investigate novel collider signatures at the HL-LHC, the FCC-hh -- studied here in this context for the first time -- and a 10 TeV muon collider. Particular emphasis is placed on electroweak boson-initiated processes, which dominate in the high-mass regime above $\sim$1 TeV. At the muon collider, these VBF-like topologies enable production even when the dipole couples to non-muonic flavors, offering a unique and sensitive probe for different flavor scenarios. We derive sensitivity projections for various theoretical benchmarks, reaching dipole couplings down to $d_{\gamma}\sim 6\times10^{-7}$ GeV$^{-1}$ at FCC-hh and $d_{\gamma}\sim 2\times10^{-7}$ GeV$^{-1}$ at the muon collider.

We calculate the contribution of extragalactic dark matter to the local dark matter density and flux in the Milky Way. By analyzing the Galactic escape velocity as a function of direction, we establish the criterion for separating dark matter particles bound to the Milky Way and those originating from the Local Group environment. Our analysis reveals that approximately 25% of dark matter particles in the Solar neighborhood have an extragalactic origin, contributing nearly 38% of the total mass flux. The directional dependence of this extragalactic component shows significant anisotropy across the sky, with implications for direct detection experiments. We provide quantitative predictions for detection rates and signatures that could help identify the extragalactic dark matter component in current and future experiments.

We show that a `QCD dilaton' field, whose vacuum expectation value sets the strong coupling, can render the Quantum Chromodynamic (QCD) confinement transition first-order. The QCD dilaton is cosmologically attracted to a false vacuum at weak coupling in the early universe. Quantum tunnelling towards the true vacuum triggers prompt chiral symmetry breaking and confinement of QCD, leading to detonating bubbles of the hadronic phase. We find that plasma sound waves produced by this dilaton-induced, first-order QCD phase transition generate a stochastic gravitational wave signal strikingly similar to the recently detected gravitational wave background from Pulsar Timing Arrays. We briefly comment on how this theory can be probed through collider experiments and cosmology.

The planned space-based gravitational wave detector, LISA, will provide a fundamentally new means of studying the orbital alignment of close white dwarf binaries. However, due to the inherent symmetry of their gravitational wave signals, a fourfold degeneracy arises in the transverse projections of their angular momentum vectors. In this paper, we demonstrate that by incorporating timing information from electromagnetic observations, such as radial velocity modulations and light curves, this degeneracy can be reduced to twofold.

Gravitational-wave astronomy provides a promising avenue for the discovery of new physics beyond general relativity as it probes extreme curvature and ultra-relativistic dynamics. However, in the absence of a compelling alternative to general relativity, it is difficult to carry out an analysis that allows for a wide range of deviations. To that end, we introduce a Gaussian process framework to search for deviations from general relativity in gravitational-wave signals from binary black hole mergers with minimal assumptions. We employ a kernel that enforces our prior beliefs that - if gravitational waveforms deviate from the predictions of general relativity - the deviation is likely to be localised in time near the merger with some characteristic frequency. We demonstrate this formalism with simulated data and apply it to events from Gravitational-Wave Transient Catalog 3. We find no evidence for a deviation from general relativity. We limit the fractional deviation in gravitational-wave strain to as low as 7% (90% credibility) of the strain of GW190701_203306.

We study the axion effects on quark matter and quark-matter cores in strange quark magnetars using a three-flavor Nambu-Jona-Lasinio model to represent the charge-parity violating effects through the axion field. Here, axions decay to gamma rays in a very strong magnetic field, which the Fermi Large Area Telescope (Fermi-LAT), Imaging X-ray Polarimetry Explorer (IXPE), and XMM-Newton will be able to detect.

We revisit the electroweak crossover of the Standard Model (SM) in the early Universe, focusing on the interplay between generalized global symmetries, magnetic flux dynamics, and baryogenesis. Employing the dimensionally reduced 3d effective field theory of the SM at high temperature, we identify the symmetry structure -- including higher-form and magnetic symmetries -- and analyze their spontaneous breaking patterns across the crossover. We further define a gauge-invariant mixing angle that interpolates between $\mathrm{U}(1)_Y$ and $\mathrm{U}(1)_\mathrm{em}$ magnetic fields. Based on this framework, we examine baryogenesis via decaying magnetic helicity and identify three key effects: the baryon asymmetry is modified by an $\mathcal{O}(1)$ factor due to (1) the gauge-invariant definition of the mixing angle and (2) the approximate conservation of the unconfined magnetic flux; (3) a novel non-perturbative process in the presence of magnetic flux, which has been overlooked in previous analyses. Our findings suggest that the previous estimation of baryon asymmetry from the magnetic helicity decay may have sizable uncertainties, and we caution against relying on it, calling for further investigation.

We present a comprehensive study of the quasinormal modes of a new class of nonlocal static and spherically symmetric black hole (BH) solutions within the framework of the revised Deser-Woodard theory of gravity. These solutions are constructed as linear perturbations of the Schwarzschild spacetime and are characterized by an inverse power-law behavior of the lapse metric function. We derive the radial profiles of the effective potentials corresponding to scalar, electromagnetic and axial gravitational fluctuations on the BH background. Using the WKB method, complemented by Padé approximants to regularize the trend of the effective potential near its peak, we compute the complex quasinormal mode frequencies associated with each type of perturbation. Our results are examined from both mathematical and physical perspectives, and are substantiated with references to existing literature. In particular, we compare the numerical outcomes with the predictions of the Schwarzschild metric to quantify deviations from the framework of general relativity. When all types of perturbations are combined, the relative deviations of the fundamental modes can reach up to $\sim 12\%$. Finally, we discuss the possibility to place observational bounds in the BH parameter space, based on the predicted sensitivities of future gravitational wave detectors.

We investigate the possibility of \emph{freeze-in} dark matter production and baryogenesis in an $E_6$ extension of the Standard Model, featuring a residual $U(1)_{\psi'}$ gauge symmetry. This symmetry arises from a linear combination of $U(1)\chi$ and $U(1)_{\psi}$, both of which are subgroups of the $E_6$. The spontaneous breaking of $U(1)_{\psi'}$ symmetry governs the dynamics of a singlet fermion, which naturally serves as a freeze-in dark matter candidate. The dark matter mass arises from dimension-five operators, and a discrete symmetry ensures their stability. We show that freeze-in production from scalar decay can yield the correct relic abundance for dark matter masses between few MeV to a few hundred GeV. Simultaneously, heavy right-handed neutrinos generate light neutrino masses via the type-I seesaw and produce the observed baryon asymmetry via leptogenesis.

We provide a fully analytical approach to calculate the nonlinearities of the gravitational waves in the ringdown of a Kerr black hole in the eikonal limit. The corresponding quasi-normal modes are associated to the orbits of a closed circular null geodesic and the problem can be analyzed by taking the Penrose limit around it. We calculate analytically the amplitude and the phase of the quadratic quasi-normal modes as well as its dependence on the black hole spin.

We show that the ionization of dense molecular clouds can be used to set strong constraints on dark matter models producing UV/X-ray photons in their annihilation or decay. We place robust and competitive constraints on various dark matter models, such as axion-like particles, scalars and sterile neutrinos, for masses between $\sim30$~eV and $10$~keV, and project forecasts to illustrate the potential of this target. We discuss how these constraints can be significantly improved by considering a more refined sample of molecular clouds near the Galactic Center and above the Galactic plane, a detailed modeling of the cosmic-ray ionization contribution and, potentially, a more refined analysis of the gas density in clouds through dust extinction maps. Thus, ionization of molecular clouds emerges as one of the most powerful tools for probing sub-keV dark matter.