Papers for Friday, Apr 10 2026

I present here PowerSpectR, an R package for computing and visualizing median-based radial Fourier power spectra from imaging data. Power spectra provide a representation of spatial structure by decomposing contributions across spatial scales, and the resulting slopes can serve as compact, low-dimensional summaries of morphological complexity across images. PowerSpectR provides a workflow for estimating these slopes, combining edge-effect mitigation through Hann windowing, Fourier-domain analysis, and radial binning with azimuthal median statistics. The use of median aggregation helps to reduce sensitivity to bright compact sources, masking artifacts, and other localized features that can bias standard estimators. PowerSpectR is released under the MIT license at \href{this https URL}{this repository}.

MAXI J1744$-$294, likely a low-mass X-ray binary system, is a Galactic-center transient source, detected at radio and X-ray wavelengths, located approximately $19''$ southeast of Sgr A*. We report the first detection of its variable linear polarization in four epochs spanning 2025 Apr 04--09. The normalized 33 and 43 GHz Stokes parameters $q$ and $u$ over the four epochs imply a common Faraday rotation screen with a rotation measure RM $=-63\,606^{+844}_{-861}$ radians m$^{-2}$, the third largest RM detected within the Galaxy. The RM is consistent with that of the Galactic center magnetar PSR J1745$-$2900, giving the first direct evidence that MAXI J1744 lies within the Galactic center region, is bound to Sgr A*, and therefore, is part of the nuclear star cluster. The uniformity in the Galactic center Faraday screen suggests that Sgr A*'s $\approx-10^5$ rad m$^{-2}$ RM is intrinsic rather than originating from an unrelated line-of-sight source. On 2025 Apr 06, we detected a secondary polarized component with an additional RM $\approx-6000$ rad m$^{-2}$, which was not seen at any other epoch. Assuming this secondary component primarily cools by synchrotron radiation, the implied local magnetic field strength is $\sim$15--30 gauss. In the context of a jetted X-ray binary progenitor, the additional RM screen and magnetic field strength are explainable with a short-lived knot in a putative jet.

Photometric missions such as Kepler and TESS have generated millions of light curves covering almost the entire sky, offering unprecedented opportunities to study stellar variability and advance our understanding of the Universe. In this data-rich environment, machine learning has emerged as a powerful tool to efficiently and accurately process and classify light curves according to their type of stellar variability. In this work, we introduce ASTRAFier: a novel Transformer-based model for variability classification that integrates Bidirectional Long Short-Term Memory (BiLSTM) and Convolutional Neural Networks (CNNs). The model operates directly on time series without requiring feature engineering, creating an easy-to-maintain and efficient end-to-end classification framework. We train and validate our model using both Kepler and TESS light curves and, respectively, achieve a classification accuracy of $94.26\%$ on Kepler and $88.22\%$ on TESS. We demonstrate scalability by deploying our model on $\sim 2.8$ million TESS light curves from sectors 14, 15, and 26 (Kepler Field-of-View) delivered by MIT's Quick-look Pipeline (QLP) and release the resulting stellar variability catalog.

Strong gravitational lensing offers a unique probe of dark matter (DM) on sub-galactic scales, where the abundance and distribution of low-mass halos are highly sensitive to the underlying properties of DM particles. In this work, we forecast LSST's sensitivity to DM substructure in galaxy-galaxy strong lenses using simulated samples and neural ratio estimators (NREs). Our simulations include both subhalos within the main deflector and line-of-sight (LOS) halos, with halo masses down to $\sim 10^7 M_\odot$ under the expected LSST ten-year survey imaging quality. We show that the constraining power on halo mass function (HMF) parameters improves significantly with sample size. Analyses based on a few hundred lenses yield broad posteriors comparable with other probes like the Ly-$\alpha$ forest. By contrast, when combining 2500 lenses, $\approx 74\%$ and $\approx 36\%$ of the prior volume considered can be excluded at the $3\sigma$ and $5\sigma$ levels respectively, enabling statistically significant exclusions of non-$\Lambda$CDM scenarios. We further demonstrate that the sensitivity arises not only from the high-mass end of the HMF but also from low-mass halos: masking halos below $\log (m_{\rm halo}/M_\odot) \leq 7.5$ induces a measurable shift in the inferred posteriors. Finally, we find that LOS halos contribute significantly to the constraining power, with increasing importance of LOS halos at higher redshifts. While this analysis assumes perfect knowledge of the data-generating process and cannot be directly applied to data analysis, it quantifies constraints achievable with LSST alone and motivates the development of robust inference methods for real survey data.

Recent JWST observations revealed two massive ($M_{\star} \gtrsim 10^{11}\,M_{\odot}$), unexpectedly large spiral galaxies at $z \sim 3$, both in overdense environments. We focus on one of these, ADF22.1 at $z = 3.09$, which hosts an active galactic nucleus (AGN), exploiting its extended [CII] emission ($\sim$30 kpc in diameter) with high-resolution observations from the Atacama Large Millimetre Array and JWST. We find a flat outer rotation curve reaching $\sim$530 km s$^{-1}$, and perform, for the first time for a system of this type, a rotation-curve decomposition. We infer a dark-matter halo mass of $\log(M_{200}/M_{\odot})=12.9^{+0.4}_{-0.3}$, a baryon-to-halo mass ratio of $0.4^{+0.6}_{-0.3}$ in units of the cosmological baryon fraction, and a ratio between the baryonic and dark-matter halo specific angular momentum of $1.0^{+0.7}_{-0.5}$. Comparing these quantities with those of local galaxies, we find that ADF22.1 is indistinguishable from $z=0$ giant discs, pointing to the inefficiency of AGN feedback in halting disc growth. Using the Mapping Nearby Galaxies at Apache Point Observatory survey, we identify potential $z=0$ descendants of ADF22.1, suggesting it will evolve into an extreme (in either mass or angular momentum) early-type galaxy. Finally, we argue that cold-stream accretion, invoked to explain disc formation at $z > 1$, cannot simultaneously account for its size, dynamical properties, high specific angular momentum, and baryon-to-halo mass ratio. Instead, sustained accretion from the hot circumgalactic medium, either via spontaneous or fountain-driven condensation, offers a more physically plausible formation pathway.

The tip of the red giant branch (TRGB) provides a key standard candle for extragalactic distance measurements and for refining the Hubble constant. We test its robustness by quantifying how metallicity, $\alpha$-element enhancement, age, and initial helium abundance modulate the TRGB luminosity, using synthetic composite color--magnitude diagrams in the $I$ and $F814W$ bands. We find that metallicity and $\alpha$-element enhancement are the primary drivers of TRGB variation, while age introduces only a modest effect and helium abundance is negligible. At fixed age and helium content, increasing the mean metallicity by 0.5 dex or the $\alpha$-element enhancement by 0.3 dex produces the well-known systematic dimming of 0.046 and 0.050 mag, respectively, in $M_I^{\rm TRGB}$, and of 0.093 and 0.044 mag, respectively, in $M_{F814W}^{\rm TRGB}$. By comparison, changes in age of 3~Gyr and in initial helium abundance of 0.10 yield minor luminosity shifts, with average changes of 0.031 and 0.009~mag, respectively, in $M_I^{\rm TRGB}$, and of 0.035 and 0.027 mag, respectively, in $M_{F814W}^{\rm TRGB}$, substantially smaller than those caused by variations in metallicity or $\alpha$-element enhancement. For mixed stellar populations under typical stellar-halo metallicity conditions, the net variation in $M_I^{\rm TRGB}$ arising from each combination of the $\alpha$-element enhancement, age, and initial helium abundance remains below 0.028~mag, well within reported systematic uncertainties. Together, these results reaffirm the TRGB as a highly robust distance indicator and support its continued use as an independent anchor for precision cosmology in the era of the Hubble-tension debate.

Black hole impacts on accretion disks in galactic nuclei can power luminous transients, but predicting their observable signatures is challenging because the post-collision flow is highly time-dependent and inhomogeneous. We present a radiative post-processing framework for relativistic hydrodynamics simulations of black hole-disk collisions. Using physically motivated prescriptions for shock heating, optical depth via an eikonal solver, and photon escape fractions that account for advection trapping and diffusion, we predict light curves and spectral energy distributions over a range of disk densities and collision velocities. Our results indicate that the emission is dominated by the long-lived, highly super-Eddington accretion flow onto the secondary black hole, rather than by cooling of the unbound ejecta. In the parameter range explored, the luminosity can reach several times the Eddington luminosity of the secondary, and the emission is generically dominated by soft X-rays. We find that lower velocity collisions produce brighter flares, while the disk surface density mainly controls spectral evolution: low-density disks typically produce keV-peaked flares with weak spectral evolution, whereas high-density disks show softer early emission and late-time hardening. A depletion-time estimate calibrated to our results suggests characteristic durations of hours to days for intermediate-mass secondaries, and yields $t_{\rm flare} \propto P_{\rm QPE}$. We discuss implications for QPE-like transients and for the SMBH-binary candidate OJ 287.

We present a weak-lensing (WL) analysis of filamentary structures in the merging galaxy cluster Abell 2744 using wide-field Magellan/MegaCam imaging data. We employ two complementary techniques: standard matched-filter scans to identify global orientations, and a new stepwise 2D tracing method to reconstruct locally varying filament orientations. The matched-filter analysis detects coherent filamentary features in the northwest and east directions across both inner (1.0-2.2 Mpc) and outer (2.2-3.4 Mpc) annuli. However, while the northwest filament yields consistent constraints across both regions, parameter inference for the eastern structure remains unstable and radially inconsistent when restricted to global reference-point scans. We demonstrate that re-characterizing the eastern structure using the locally preferred elongation directions from our stepwise tracer significantly resolves these tensions, improving fit quality and bringing inner and outer constraints into agreement. Furthermore, the detected filaments align well with diffuse X-ray structures and previously identified merger axes, supporting their physical connection to the cluster's mass assembly. These results highlight that stepwise 2D tracing is essential for characterizing curved or complex filaments where global reference-point scans are insufficient.

Tracing the compositional link between terrestrial super-Earths and their host stars provides clues to their dominant formation pathway. By constraining the stellar abundances of refractory elements, we can predict the core mass fractions (CMFs) of their super-Earths. The level of agreement between this prediction and the planetary CMF derived from their masses and radii can reveal past formation processes, like mantle stripping and water-rich formation plus sequestration in the planet's core. Here, we present the first results from the Near Infrared Planet Searcher (NIRPS) GTO CMF subprogram: an intensive radial velocity campaign to refine masses and compute host stellar abundances of three hot super- Earths around M dwarfs (GJ 1132 b, GJ 1252 b, and LTT 3780 b), calculating masses of $1.69 \pm 0.15M_\oplus$, $1.54 \pm 0.18M_\oplus$, and $2.34 \pm 0.10M_\oplus$ respectively. We measure the CMFs of these and six further hot super-Earths with precise masses already available in the literature to 10-15% precision. We compare these to CMF predictions made from measuring the Fe, Mg, and Si abundances of their host stars measured from the NIRPS spectra. We find that the CMFs of these planets are smaller than expected from their host stellar abundances, to a statistically significant degree. This discrepancy is suggestive of significant reservoirs of water, and while these planets are too hot to harbor surface water, they likely have interior water mass fractions of $\sim$1%.

The sub-kpc scale gas structure encodes key information of giant molecular cloud (GMC) formation. Therefore, we aim for a quantitative description of molecular gas structure across 150-1000 pc using a sample of 8984 GMCs from 40 galaxies observed by PHANGS-ALMA. We homogenize our data to a fixed resolution of 150 pc and mass sensitivity of 2.5 M$_{\odot}$ pc$^{-2}$ to remove observational bias. We then calculate nearest neighbour distances, neighbour number density, and two-point correlation functions for the catalogued GMCs. When analysing the two-point correlation function, we generate several control samples that reflect different null hypotheses on large spatial scales. We stack integrated intensity CO emission profiles around the position of catalogued GMCs to probe the gas distribution on scales between the resolution limit and the typical GMC-GMC spacing. Our measurements of cloud spacing and number of neighbours show that GMC clustering follows the large-scale gas distribution. Once we account for this contribution, the peak excess clustering in the two-point correlation function drops from 1+$\omega$ of 2.3 to 1.3, with the power-law slope flattened from -0.25 to 0. We show that the stacked CO intensity profiles around CO peaks can be recovered by the "GMC size" measured by CPROPS, with an additional 20% of the flux in an extended component beyond 500 pc. We find that our stacked profiles can be fit with a double Gaussian function plus a constant offset. The broad Gaussian component accounts for 70% of the over-density power above the constant offset, and is stronger around massive and gravitationally bound GMCs. Our results indicate that galactic structure regulates the GMC distribution in galaxy disks, and the formation of massive, gravitationally bound GMCs is related to strong local gas clustering.

Modeling subphotospheric shocks in a gamma-ray burst (GRB) is challenging due to the various timescales that must be resolved, and the fact that the same radiation dynamically mediates the shocks while forming the observed signal. Here, we present the first self-consistent radiation-hydrodynamic simulation of a subphotospheric internal collision, following the system from formation and propagation of forward and reverse radiation-mediated shocks all the way to photon decoupling and free streaming toward the observer. The simulation evolves the plasma and photon field with full Compton coupling, including the feedback on the hydrodynamic flow. As the ejecta expands and the optical depth decreases, both shocks broaden and the radiation field becomes highly non-thermal. Surprisingly, we find that the reverse shock remains completely radiation-mediated down to upstream optical depths of order a few $\times 10^{-1}$, which indicates that Compton coupling is important even in moderately optically thin regions. The photons undergo last scattering over a broad range of radii rather than at a single photospheric surface. The light curve shows a late, quasi-thermal post-cursor produced by photons that decouple upstream of the reverse shock, which could be searched for in observations. The emitted time-integrated spectrum is GRB-like, with a low-energy photon index $\alpha \sim -1$ and a high-energy photon index $\beta \sim -2.5$. These results show how radiation-mediated shocks evolve close to the photosphere and how they shape the emitted photon field.

Gravitational-wave observations have revealed an excess of binary black hole mergers with primary masses near $\sim 35\,M_\odot$. We show that if this feature originates from dynamical formation in dense stellar systems, and if the pair-instability supernova truncates the first-generation black hole mass spectrum, then second-generation mergers inevitably produce a second peak near $\sim 70\,M_\odot$. This structure reflects the suppression of first-generation black holes above a characteristic mass and the accumulation of merger remnants near twice that scale. Its location is robust, whereas its amplitude depends strongly on cluster initial conditions. Using a large suite of cluster population-synthesis models, we show that current gravitational-wave data already constrain the birth properties of globular clusters, irrespective of their overall contribution to the observed population. If clusters dominate mergers above the pair-instability scale, these constraints tighten further and imply a minimum first-generation merger rate of $\mathcal{R}(m_1 \leq 50\,M_\odot) \geq 0.099\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ ($99\%$ confidence). We further show that a drop or gap in the secondary black hole mass spectrum is not a robust signature of a cluster origin for high-mass mergers within the pair-instability mass gap. A confirmed excess near $\sim 70\,M_\odot$ would support a dynamical origin of the $\sim 35\,M_\odot$ feature and provide independent evidence for a pair-instability mass gap with a lower edge at $\lesssim 50M_\odot$.

We investigate how well the Large Interferometer for Exoplanets (LIFE) mission concept can detect habitable conditions on exoplanets through the presence of atmospheric water vapor as a proxy for surface oceans. We model the atmosphere of a pre-biotic Earth-like planet across a range of water concentrations, from water-poor to water-rich, with surface partial pressures from 10$^{-7}$ to 1 bar of H$_2$O. We simulate LIFE-like noise at spectral resolutions R = 50 and 100 using LIFEsim and perform Bayesian atmospheric retrievals to determine the technical requirements for LIFE to confirm habitability. We model three vertical water distributions: a vertically constant profile, a Manabe-Wetherald based Earth-like profile, and a diffusion and photochemistry profile to test how the assumed vertical structure influences the retrieved abundances. Clouds are not modeled. We find the ability for LIFE to detect water strongly depends on the vertical profile assumed. LIFE is unable to constrain the highest water cases and provides upper limits on low water planets. For the highest water abundances, absorption features saturate and reduce sensitivity to characterize precise H$_2$O levels. Water vapor is not detectable in any profile modeled for $\leq10^{-6}$ bar in surface water, comparable to Mars. For an Earth-like profile, LIFE could constrain H$_2$O concentrations from $\sim10^{-3}$ to 1 bar, spanning below and above present-day Earth concentrations of 10$^{-2}$ bar. Detectable atmospheric water may imply surface oceans, as water is highly reactive and rapidly removed by surface mineral reactions. Thus, LIFE can characterize water abundances indicative of habitable surface conditions.

Observations of temperate sub-Neptunes suggest active chemical environments, finding evidence of both water vapor and photochemical hazes in their atmospheres. Hazes formed in water-rich atmospheres are chemically complex, containing molecules relevant to prebiotic chemistry, and their strong optical opacity obscures sought-after gaseous molecular absorption features. While many studies have investigated haze formation and properties across diverse atmospheric conditions, little is known about the evolution of these hazes in their environment once formed. In particular, interactions with water can drive hydrolysis reactions that alter haze composition and optical behavior, affecting our interpretations of habitability and observational spectroscopy. Here, we perform hydrolysis experiments on haze analogs of temperate water-rich exoplanets and measure their optical properties. Transmittance measurements from 0.4 to 28.5 $\mu$m reveal changes in key functional groups after hydrolysis, along with an overall increase in sample absorbance. We report the derived optical constants for use in observational and modeling studies. Through synthetic atmospheric spectra, we demonstrate the need for physically informed haze optical properties in models, consistent with expected planetary conditions. The increased absorptivity and high imaginary refractive index of hydrolyzed hazes almost completely flatten features in model spectra, presenting critical consequences for atmospheric characterization of water-rich sub-Neptunes.

Amidst the exoplanet revolution in which multiple techniques have successfully found planets, the Doppler (Radial Velocity, or "RV") technique is unique in its sensitivity to giant planets at very long orbital periods around Sun-like stars. The upcoming retirement of Keck-HIRES will incur irreversible changes in the continuation of HIRES's decades-long stable RV baseline and with it, the exoplanet community's ability to detect giant exoplanets with periods longer than Jupiter. With the time elapsed from the last HIRES RV for many stars of interest at ~3 years and growing, we tested the impact of a "critical RV", one that would bridge this gap between past HIRES RVs and future stable Keck-KPF RVs, on the recovery of long-period giant exoplanets. We generated 2000 1-planet systems with RVs sampled at a timeseries representative of this situation and used the planet-finding code Octofitter to perform injection-recovery experiments including and omitting this critical RV for each system. For these injected long-period super-Jupiters (~8-55 years, 1-13 $M_J$), including the critical RV induced a $1.5\times$ enhancement in overall planet recovery and a more specific $3.5\times$ enhancement in the recovery of super-Jupiters with Saturn-like periods. These experiments show that gathering a critical RV for stars of interest can help ensure that HIRES's decades-long stable RV baseline in conjunction with future KPF RVs, or indeed that the RV baselines containing an observational gap of any instruments that will undergo an RV zeropoint offset, will continue to be foundational to the discovery of long-period giant exoplanets in years to come.

Short-period white dwarf+main-sequence binaries are Post-Common-Envelope Binaries (PCEB) that have survived a common envelope phase. Such systems, if detached and eclipsing, enable precise measurements of the constituent stars, providing a unique opportunity to probe the effects of the common envelope phase on the system. We report the discovery of one such nearby (57 pc) system, TIC-460388167, using a combination of multi-band photometric light curves and spectroscopic radial velocities. In addition to eclipses, the system exhibits a continuously variable light curve that we model as a combination of ellipsoidal variations and star spots. We determine a period $P$=0.63596258$\pm$0.00000012 d and inclination $i$=89.0$\pm$0.4 deg. The best-fitting model specifies a white dwarf with T$_1$=7607$\pm$127 K and radius R$_1$=0.0131$\pm$0.0003 $R_\odot$, which is eclipsed by a T$_2$=3151 $\pm$ 59 K, R$_2$=0.327$\pm$0.006 $R_\odot$ M dwarf. The white dwarf mass is 0.61$\pm$0.04 M$_\odot$. We present the first velocity resolved profile for a PCEB secondary and show that the rotation of the M-dwarf is synchronous with the orbital period, as expected. We compare the constituent stars to other PCEB systems and find TIC-460388167A is one of the coolest known white dwarfs in such systems. TIC-460388167 is among the longest period eclipsing PCEB systems known.

Understanding the initiation of solar flares and coronal mass ejections (CMEs) is essential for improving forecasts of extreme space weather. Soft X-ray (SXR) and Extreme Ultraviolet (EUV) observations provide critical diagnostics of the highly dynamic solar corona; however, significant measurement gaps persist despite decades of observations. The Solar Extreme Ultraviolet Spectrograph and High-energy Imager (SEUSHI) aims to address these gaps by combining multi-pinhole SXR imaging with grazing-incidence EUV spectroscopy on a shared camera. SEUSHI delivers spatially-resolved temperature and emission measure maps at 1 arcminute resolution and 5 second cadence to identify Hot Onset Precursor Events (HOPEs), which provide early alerts of flares. Additionally, high-cadence (100 Hz) readouts of selected image rows enable photon-counting spectroscopy over 1.1-6.8 keV with approx. 0.08 keV energy resolution, to investigate elemental abundance evolution in active regions, a key diagnostic of coronal heating. SEUSHI also provides high-resolution (approx. 0.2 nm) EUV spectra measurements across the 16.1-33.8 nm range at 5 second cadence for studies of coronal dimming and generation of early alerts for CMEs. SEUSHI is designed with low power, mass, and volume requirements, making it suitable for small satellite platforms. A technology demonstration version of SEUSHI is currently under development for flight aboard the Solar Dynamics Observatory Extreme Ultraviolet Variability Experiment calibration sounding rocket. This paper presents the instrument design, development, and calibration. Successful demonstration on the sounding rocket platform is an important step towards the opportunity to fly SEUSHI on future satellite missions and thus to improve space weather operations.

Magnetic cataclysmic variables provide a natural laboratory for studying how accretion interacts with compact-object magnetospheres and generates stochastic variability. We present an optical variability study of the intermediate-polar candidate EP240309a, an Einstein Probe X-ray transient, using BOOTES photometry, high-cadence TESS light curves, and a SOAR/Goodman optical spectrum. Previous studies found a white-dwarf spin period of 3.97 min (Pspin ~ 238 s) and an orbital period of Porb = 3.7614(4) h. Power spectral densities from the BOOTES data are consistent with single power laws with slopes alpha ~ 1.2-1.8, with no statistically significant evidence for a bend across the sampled frequency range. Using red-noise simulations and injection-recovery tests, we place one-sided constraints on any putative break frequency, which translate, under standard dynamical identifications, into an upper limit on the magnetospheric radius of Rm <= few x 10^10 cm for MWD = 0.8 Msun. In the TESS data, we detect a linear rms-flux relation on hour timescales in three high-cadence sectors, while two other sectors do not show a robust detection, indicating epoch-dependent rms-flux behavior. The SOAR spectrum shows Balmer and He II emission lines with FWHM about 1000-1600 km s^-1; under a Keplerian interpretation, these imply characteristic radii of r about (0.9-3.4) x 10^10 cm, broadly comparable to the timing-based constraints. Overall, the data provide conservative, order-of-magnitude radius constraints consistent with accretion onto a magnetic white dwarf, but they do not establish the detailed accretion geometry or exclude stream-fed or mixed accretion scenarios.

We present a concise multimessenger inference of the neutron-star core sound-speed profile using GW170817 and three \textit{NICER} mass--radius posteriors (PSR J0030$+$0451, PSR J0740$+$6620, and PSR J0437$-$4715). The main result is not only a preference for intermediate-density stiffening within smooth equation-of-state families, but a translation of that inference into the stellar masses that access the relevant density regime. In the baseline smooth peaked family, the posterior probability that $c_s^2 > 1/3$ at $3.5\,n_{\rm sat}$ is $85.4\,\%$, while equal-prior averaging over peaked, monotonic, piecewise, and transition-capable families gives a more conservative $79.0\,\%$. Posterior-resampled exact Tolman--Oppenheimer--Volkoff solutions show that the onset density of the inferred stiffening is typically reached near $1.6\,M_\odot$, whereas the peak region is accessed only near $2.1\,M_\odot$. A J0740-like $2.07\,M_\odot$ pulsar reaches the onset in $91\,\%$ of posterior draws but the peak in only $46\,\%$, showing that current data mainly constrain whether massive stars have entered the stiffening regime rather than traversed its full peak.

Changing-look transitions challenge our understanding of active galactic nuclei (AGN), exhibiting dramatic changes in broad-line emission and continuum flux on timescales of months to years. We present a detailed study of the spectroscopically confirmed changing-look AGN ZTF18abuamgo. Combining photometric survey data with spectroscopy spanning three epochs over 20 years, we identify a turn-on transition from a Type 1.5 to Type 1.2 AGN and estimate the timescale of this change to be as short as four years. Spectral analysis indicates that this transformation is driven by a rapid increase in accretion rate, with the Eddington ratio rising from $0.032 \pm 0.005$ in the dim state to $0.08 \pm 0.01$ in the bright state. For the first time in a changing-look AGN, we apply the Boltzmann plot method to the visible Balmer series emission, deriving broad line region electron temperatures of $11,800 \pm 900$ K and $11,900 \pm 2,400$ K in 2022 and 2024, respectively. Applying single-epoch black hole mass estimation to the brightening H$\alpha$ emission, we find a mass of $(5.0 \pm 0.4) \times 10^7 M_\odot$. The consistency in this estimate across all spectroscopic epochs suggest that even highly variable broad lines in CL-AGN do not bias the results derived using this method. Our results demonstrate that objects like ZTF18abuamgo provide a unique laboratory to study extreme AGN variability, probe the physical conditions in the broad line region, and assess the limitations of widely used black hole mass estimation methods.

We investigate two neighboring clusters in the Cygnus complex, Berkeley 86 and Berkeley 87, with a primary emphasis on the evaluation of extinction in the field of view towards and across the clusters. We also analyze their kinematic behavior in space and time to discern their possible common origin and relation to the Cyg~OB1 association. New CCD photometry in the Vilnius seven-color system, obtained down to V=19.0 mag in the fields of these two clusters, is used to classify stars in terms of spectral and luminosity classes and to determine the individual values of interstellar extinction. The probable cluster members are identified in a 5-parameter space based on Gaia DR3. The cluster ages and stellar masses are derived through the use of the HR diagrams. To obtain the 3D kinematics of the clusters and trace their orbits back in time, we combine the Gaia-based proper-motions and distances with radial velocities from the literature. The estimated cluster properties show that both clusters are almost equidistant (1.7 kpc) and nearly coeval, with average ages of 6.1$\pm$0.5 and 6.5$\pm$0.4 Myr, respectively, and age dispersion of 3 Myr. The nonuniformity of extinction is evident within each cluster, especially pronounced across the face of Berkeley 86 where the most-massive stars show substantial substructure. By extrapolating the observed mass function to a minimum stellar mass, we obtain cluster masses of 519 M(Sun) and 1551 M(Sun) for Berkeley 86 and 87, respectively. Although both clusters share very similar properties, their orbital paths show no indication that they had a common birthplace, however Berkeley 87 and its neighbor NGC 6913 are very likely to have been born in pair.

The forthcoming launch of the Ultraviolet Transient Astronomy Satellite (ULTRASAT) will transform our understanding of the transient ultraviolet sky by increasing our ability to identify transients due to its unprecedented 204 deg2 field of view. While rapid (extragalactic) transients are a priority science area for the mission, flaring stars and AGN can often contaminate searches for such objects. To prepare for these challenges, the Zwicky Transient Facility (ZTF)-ULTRASAT experiment observed five fields at high cadence over three nights, in close proximity to ULTRASAT's three northern high-cadence fields. A real-time filter identified seven transient candidates, of which five were presistent variable sources and two were spurious. Periods and amplitudes derived from the ZRF Source Classifcation Project (SCoPe) showed that three candidates were RR Lyrae stars with short periods and high amplitudes, while the remaining two displayed flairing behavior. We demonstrate that short-timescale, high-amplitude variables can sysematically mimic transient alerts in high-cadence UV surveys, and we provide a concrete startegy to this contamination using pre-existing machine learning catalogs.

Gamma-ray bursts (GRBs) are among the most luminous and rapidly evolving transients in the Universe. While space-based instruments have extended GRB observations up to energies of $\sim$100 GeV, the detection of very-high-energy (VHE; $E>100$ GeV) emission from ground-based telescopes, especially during prompt or/and the early afterglow phase, remains challenging. These difficulties arise from the rapid temporal decay of GRB afterglows, strong attenuation by the extragalactic background light (EBL), observational latency, and the typical poor sky localization provided by MeV-detectors such as Fermi/GBM. In this work, we investigate the prospects for detecting TeV ($\sim$100 GeV--1 TeV) emission from poorly localized GRBs by adopting optimized follow-up strategies based on rapid tiling of large localization regions. We simulate a realistic population of GRBs informed by more than fifteen years of Fermi/GBM and Swift/XRT detections and recent progresses in the afterglow emission modeling. Using these simulations, we evaluate the detectability of GRB early afterglows by the next-generation Imaging Atmospheric Cherenkov Telescopes, equipped with larger field-of-view (FoV), as a function of latency, exposure time, and observational strategy. Our strategy can significantly enhance the detection rate; for instruments such as ASTRI and LACT, it increases by up to a factor of two compared to strategies limited to well-localized (Swift-like) events. For CTAO, our proposed approach provides up to four VHE detections per year.

We present a study of the multiphase gas structure and kinematics of the circumnuclear disk (CND) of NGC 5506, a nearby edge-on Seyfert galaxy, at a spatial resolution of $\sim20$ pc. Observations of [C I](1-0), CO(3-2), and HCO$^{+}$(4-3) obtained with the Atacama Large Millimeter/submillimeter Array reveal the CND dominated by rotational motion on scales of several hundred parsecs. No significant differences in geometrical thickness or velocity structure are found between [C I](1-0) and CO(3-2) across the CND, whereas HCO$^{+}$(4-3) emission is more concentrated toward the disk plane. The ratio of velocity dispersion to rotational velocity, a proxy for disk scale height-to-radius ratio, is high ($\gtrsim0.9$) in the central region ($\lesssim30$ pc) for both [C I](1-0) and CO(3-2), indicating geometrically thick structures in both tracers. Regions where the [C I](1-0)/CO(3-2) ratio exceeds the CND average are spatially correlated with the [O III]$\lambda$5007 bicone observed with the Hubble Space Telescope, suggesting that CO is preferentially dissociated by the AGN-driven biconical ionized outflow. The observed CND scale height and velocity dispersions traced by [C I](1-0) and CO(3-2) are consistent with a model in which supernova-driven turbulence provides the vertical support for the CND.

The central few parsecs of the Milky Way host a complex X-ray-emitting environment in which several extended plasma components are blended along the line of sight, complicating attempts to measure the intrinsic properties of individual components. In particular, the supernova remnant (SNR) Sgr A East is strongly confused with the stellar wind-fed plasma associated with Sagittarius A* and the surrounding nuclear environment. Here we apply Poissonian Generalized Morphological Component Analysis (pGMCA) to deep, stacked Chandra ACIS-I observations of the Galactic Center to disentangle these overlapping X-ray components. By comparing the separated X-ray components with multiwavelength data, we identify the location of the reflected shock in Sgr A East and construct spatially resolved maps of Fe and S/Ar/Ca emission. The Fe emission is centrally concentrated, consistent with the properties of mixed-morphology supernova remnants. Separating the SNR emission from the shocked wind plasma around Sgr A* allows us to recover uncontaminated SNR properties and improve the robustness of the derived parameters. Spectral modeling of the isolated Sgr A East component reveals a lower ionization age and a higher electron density than previously reported, indicating strong interaction with dense surrounding material.

The European Space Agency's Ariel mission will conduct a survey of the atmospheric properties of exoplanets around bright stars. The mission is nominally divided into three Tiers. The Tier 1 survey will consist of low-precision observations of ~1000 planets, with a subset of these included in the higher-precision Tier 2 survey expected to be necessary for atmospheric characterization. Tier 3 will be repeated observations of a small number of benchmark planets. Though previous studies have assessed the ability of Ariel to uncover population-level trends, they have generally presupposed a given Tier. Here we interrogate this assumption and assess the information content of Ariel transmission spectra as a function of Tier for three benchmark planets: a hot-Saturn, warm-Neptune, and temperate sub-Neptune. We simulate a grid of Ariel transit spectra at different Tiers for each target and use retrievals to assess which chemical species are detectable. We find that for giant planets like a hot-Saturn or warm-Neptune, Tier 1-quality observations are sufficient for <1.5dex constraints on H2O and CO2, irrespective of the presence of clouds -- meaning important chemical insights are already obtainable in the Tier 1 survey. Moving to Tiers 2 and 3 result in an incremental increase in precision as well as other molecules becoming detectable in certain scenarios (e.g., H2S, CO). Tier 1 observations are also sufficient to constrain CH4 in a cloud-free, temperate sub-Neptune, whereas observations with at least Tier 2 precision are necessary if the atmosphere is cloudy. The number of transits necessary to reach this precision, however, may be prohibitive for the inclusion of temperate sub-Neptunes in even the Tier 1 survey.

Galaxy environment plays a crucial role in quenching star formation in dwarf galaxies. In Milky Way (MW)-like environments, dwarf satellite quenching is primarily driven by ram pressure stripping (RPS), the direct removal of satellite gas by the host halo gas. Using a suite of 20-pc resolution hydrodynamical wind tunnel simulations, we constrain the satellite mass scale at which the stripping of a dwarf galaxy's interstellar medium (ISM) becomes inefficient in MW-like halos. The simulations include radiative cooling in a multiphase satellite ISM, star formation, and stellar feedback, and vary both satellite masses ($M_{\star}=10^{6.2}, 10^{6.8}, 10^{7.2}\ M_{\odot}$) and host halo gas densities along a first-infall and post-pericentric orbit. We find that the degree of ISM stripping in our dwarf galaxies is consistent with the analytical prediction by McCarthy et al. (2008). Star formation is rapidly quenched when RPS is effective, but can be mildly enhanced or temporarily quenched and subsequently reignited when RPS is incomplete. ISM stripping is efficient for satellites with $M_{\star} \lesssim 10^{7}\ M_{\odot}$ (or $M_{200} \lesssim 10^{10}\ M_{\odot}$) but highly inefficient above this scale. This transitional mass ($M_{\star} \approx 10^{7}\ M_{\odot}$) is 0.5-1 dex lower than that found in observations and cosmological simulations, suggesting that additional mechanisms are needed to quench more massive satellites, such as tidal stripping of the satellite dark matter or RPS from a clumpy gaseous halo.

Globular clusters (GCs) are powerful tracers of galaxy assembly, frequently used to identify accreted substructure and reconstruct hierarchical merger histories. With advances in GC formation models and cosmological simulations, we can now better quantify the information about galaxy evolution encoded in present-day GCs. Here, we investigate how GC kinematics evolve over cosmic time and assess the extent to which GCs retain memory of the past of their host galaxy. Using a GC formation model applied to five Milky Way (MW) analogues from the Latte suite of the FIRE-2 simulations, we track the evolution of kinematic properties. At $z=0$, in-situ and ex-situ GCs exhibit substantial overlap in kinematic space, indicating that these populations are not clearly separable. We find that a subset of kinematic properties evolve in an ordered fashion across both in-situ and ex-situ populations, whereas others are dominated by stochastic variations. As a result, by the present day, most memory of the progenitor of an accreted GC is erased and only a few correlations persist. These correlations link progenitor halo mass to the total mass and number of a GC population, and the galactocentric distance of GC substructure to progenitor maximum circular velocity. These results highlight how both deterministic and stochastic processes driven by galaxy evolution shape GC kinematics and demonstrate the limits of reconstructing the assembly history of a galaxy from present-day GC orbits alone.

Wide-orbit ($>$10 AU) gas giant planets shape the architecture of planetary systems, yet their occurrence rate remains poorly constrained. JWST has obtained the deepest mid-infrared images of nearby stars to date through substantial MIRI time-series observations of transiting planets, providing sensitive probes for wide-orbit companions. Here we leverage 15 micron observations from four programs targeting ten M-dwarf systems to search for such planets. By applying reference differential imaging for precise PSF subtraction, we achieve a median 5$\sigma$ contrast of $8.9 \times 10^{-4} - 6.2 \times 10^{-3}$ (median sensitivity in apparent magnitude of 15.8-16.8 mag) at a separation of 1" and $1.2 -9.1 \times 10^{-4}$ (17.5-19.0 mag) at separations $\gtrsim$3". The sensitivity is converted to planet detection probability for each system as a function of planet mass versus semimajor axis. Assuming solar metallicity and a clear atmosphere, we are sensitive to Jupiter-sized planets with an effective temperature of ${\sim}170$ K at separations beyond 35 AU in systems at 12.5 pc. Additionally, we catalog the nearby sources and estimate their possible impact on future observations assuming they are background sources. Our results demonstrate that archival MIRI time-series imaging data is a powerful window into the population of wide-orbit gas giants around M-dwarfs.

The magnetar 1E 1841-045 exhibited a new active episode starting on August 20, 2024, marked by X-ray bursts and enhanced persistent emission. Using data from the Einstein Probe (EP), we report on the timing and spectral results following the onset of this outburst. The pulse profile displays a multi-peaked structure, with notable phase shifts in the secondary peak. Energy-resolved pulse profile analysis indicates a transition in the dominant peak of the pulse profile above 5.8 keV. The 0.5-10 keV X-ray spectrum is well-modeled by a combined blackbody and power-law (BB+PL) model, showing a $\sim 20\%$ flux increase following the outburst. Phase-resolved spectroscopy indicates a correlation between BB temperature and pulse profile intensity, along with spectral hardening at a specific pulse phase. The high spatial resolution of EP enables effective separation of the supernova remnant emission, which is crucial for measuring the intrinsic pulse emission of the source. These findings underscore the intricate relationship between magnetar outbursts, pulse profile evolution, and spectral characteristics.

We report new results from a nulling study of PSR~B0751+32 (PSR J0754+3231), observed at 1250~MHz with the Five hundred meter Aperture Spherical radio Telescope (FAST). Our analysis confirms the presence of periodic nulling in this pulsar. Using the recently developed mixture model method, we obtained a nulling fraction (NF) of $35.1\% \pm 0.6\%$. Three independent approaches were employed to estimate the nulling periodicity, and the results reveal significant temporal evolution of the modulation both within individual observations and across different \textbf{observing} essions. The pulsar exhibits an asymmetric two-component mean pulse profile, with the leading component brighter and narrower than the trailing one. Pulse energy analysis shows that both components remain stable immediately after the onset of the burst state, but subsequently undergo a progressive decline, with the trailing component most severely affected prior to burst termination. Notably, no evidence of the previously reported subpulse drifting was detected in our data. Our results challenge previous models that ascribed periodic nulling to purely geometric effects.

The enigmatic origin of the diffuse neutrino background detected by IceCube in the energy range from TeV to PeV remains one of the central open problems in high-energy astrophysics, and this puzzle is further deepened by the recent evidence for a spectral break. Could this convex-spectrum background arise predominantly from the evolution of active galactic nuclei (AGNs)? In this work, we claim that the spectral break is naturally predicted when AGN evolution is taken into account, and the diffuse flux can be interpreted as the superposition of contributions from AGNs at different evolutionary phases. We develop a unified framework that incorporates AGN evolution, where cosmic rays (CRs) accelerated during the active phase subsequently diffuse and interact in the host galaxy after the central engine switches off, producing a long-lived hadronic afterglow. Adopting physically plausible parameters, our model successfully reproduces both the spectral features of the diffuse background and the observed neutrino emission from representative sources such as TXS 0506$+$056 and NGC 1068. Our results suggest that AGN host galaxies are more efficient CR reservoirs than previously expected. Moreover, the model favors a lepton-dominated scenario for most AGNs. This conclusion accounts for the relatively low detection rate of point-like sources by IceCube and underscores the need for next-generation neutrino telescopes with larger effective areas and higher sensitivity.

We calculate the gravitational wave signal from the collapse of a rotating 300 $M_\odot$ star at the upper end of the pair-instability regime. The large-scale asymmetries that develop during the collapse produce a strong signal in the deci-Hz range that has a characteristic shape which is likely amenable to a template-based search. The most ambitious designs for deci-Hz detectors could detect such signals out to distances of 200 Mpc, possibly at a rate of 0.5 per year.

Meteor clusters are typically defined as groups of meteors that appear close together in both space and time. To date, only a handful of such events have been recorded instrumentally and analysed in detail. In many documented cases, thermal stress has been identified as the most likely cause of meteoroid fragmentation near Earth. This paper documents two further cases and provides a summary of all currently known clusters. The two clusters that were recorded over Hawaii Island in 2023 and 2024 represent two distinct scenarios. The 2024 meteor cluster was characterised by a dominant mass body and, with the fragments arranged along the antisolar direction according to their mass. Such cases enable us to reliably determine the age of the cluster and identify the most likely formation scenario. This cluster was around three days old, and the thermal stress was the most likely mechanism of its formation. The 2023 cluster was not such a case. It does not contain a mass dominant body, nor are its fragments arranged by their mass. Therefore, it was only possible to estimate its age to be no more than four days. Furthermore, other potential formation mechanisms besides thermal stress cannot be ruled out. This fact was observed in all analysed clusters. All clusters known up to date were formed in close proximity to Earth. The volume of a cluster increases with its age. This means that older clusters, formed by the fragmentation far away from Earth may remain undetected, as their fragments are also dispersed too widely to be observed by local experiment. However, global networks can detect such dispersed clusters.

Chang'E-4 (CE4), the first mission to soft-land on the lunar farside, provides a unique opportunity for astronomical observations from an environment shielded from terrestrial radio interference, and thus serves as pathfinder for lunar farside radio search for extraterrestrial intelligence (SETI) studies. We present a search for periodic technosignatures using low-frequency radio observations from the CE-4 mission, the first radio SETI study based on data from on the observation in lunar farside. We analyze the CE4 dynamic spectra with a component-level framework that combines principal component analysis (PCA), cross-antenna basis alignment, as well as temporal periodicity and frequency comb structure diagnostics. No final periodic candidate signal is found after the selection procedure, and we therefore find no evidence in the present CE4 sample for a credible periodic artificial signal. This study serves as a pathfinder and provides a practical framework for lunar radio SETI analysis. As more future lunar missions begin to incorporate radio instrumentation, lunar farside may become a promising site for expanding radio SETI research.

We present detailed light-curve analyses of ten binary-lens microlensing events observed during the 2023--2025 seasons and selected as candidates for hosting brown-dwarf companions. The sample includes OGLE-2023-BLG-0249, KMT-2023-BLG-1246, OGLE-2023-BLG-0079, KMT-2024-BLG-0072, KMT-2024-BLG-0897, KMT-2024-BLG-1876, KMT-2024-BLG-2379, KMT-2025-BLG-0922, KMT-2025-BLG-1056, and KMT-2025-BLG-2427. For each event, we carry out modeling of the light curve, explore relevant degeneracies, and, when finite-source effects are present, determine the angular Einstein radius. For OGLE-2023-BLG-0249, we additionally measure the microlens parallax, which allows a direct determination of the lens masses and distance. For the remaining events, we estimate the physical lens properties via Bayesian analyses incorporating Galactic priors. The resulting posteriors show that the lens companions in all systems have median masses in the brown-dwarf regime, and the lenses of two events (KMT-2025-BLG-0922 and KMT-2025-BLG-1056) are consistent with binaries in which both lens components fall within the brown-dwarf mass range. Spanning a wide range of projected separations and distances, these detections illustrate the power of high-cadence microlensing surveys to build a census of brown-dwarf companions, including faint and distant systems beyond the reach of flux-limited methods.

The rotation-activity connection is the standard paradigm for interpreting chromospheric and coronal activity in late-type stars, namely, stars with outer convection zones. This paradigm states that activity increases with decreasing rotation period until a saturation limit is reached. By scaling rotation periods with the convective turnover time via the Rossby number, $\text{Ro}$, saturation is expected to occur at a universal value across all spectral types. In our paper, we systematically investigate the relationship between rotation and activity as measured though X-ray emission for a large sample of late-type stars to test the universal applicability of this paradigm. To this end, we utilized TESS short-cadence space photometry to determine the rotation periods for late-type stars identified in the eROSITA all-sky survey. This combined dataset provides rotation and X-ray measurements for 14004 stars, representing a sample size increase of more than an order of magnitude compared to previous studies. We find that the convective turnover times derived from this sample closely agree with theoretical computations, supporting the idea that Rossby number-activity relations hold for all late-type main sequence stars. The lower level of activity in earlier spectral types (e.g., F-type and late A-type stars) is a physical consequence of their short convective turnover times, which prevent them from rotating rapidly enough to ever reach the saturation regime. We demonstrate that a simple model incorporating our derived turnover times versus color can successfully reproduce the observed characteristics of the eROSITA X-ray activity distribution, as measured by the L$_X$/L$_{\text{bol}}$ ratio and {\it Gaia} BP-RP color.

Bayesian methods for the detection of continuous gravitational waves (CGWs) in Pulsar Timing Array (PTA) data incur substantial computational costs that grow rapidly due to the number of noise and signal parameters characterizing the fitted model being proportional to the size of the PTA. This computational burden limits the scalability of these methods for large-scale PTAs comprising hundreds of pulsars anticipated from next-generation radio astronomy facilities. In this work, we introduce a computationally efficient frequentist method designed to circumvent this challenge. This is achieved by combining an adaptive spline fitting algorithm that non-parametrically suppresses red noise, thereby eliminating the need for complex noise modeling inherent to Bayesian methods, with a novel scheme for optimizing the subsets of pulsars included in the search. We quantify the performance of our method on a simulated dataset based on the NANOGrav 15-year data release and find that it achieves a performance comparable to that of Bayesian analysis: for a CGW signal with a signal-to-noise ratio of $\approx 10$, our method yields a relative characteristic strain error of 1.0\% and a frequency error of 0.072\% from the injected values by using the optimal pulsar selections, while the same errors are 1.7\% and 0.16\%, respectively, for the standard Bayesian analysis. At the same time, our analysis completes in less than 5 hours, in contrast to the 1-2 days required by Bayesian methods. This allows us to perform a rigorous study of our method using multiple data realizations and signal parameters, establishing it as an efficient and scalable tool for CGW searches with large-scale PTAs.

According to activity-rotation relations, rapid rotators are expected to show high levels of magnetic activity. However, recent studies with TESS have found Ultra Fast Rotating (UFR) M dwarfs with periods $<1$ d displaying low levels of flaring activity. There have been efforts to explore their magnetic field strengths through spectropolarimetric measurements and to assess the potential for binarity. However, neither could fully explain the lack of observed flaring activity despite their rapid rotation. Another avenue for investigation is to measure their coronal emission for signs of supersaturation: an underluminosity in X-rays observed for some rapidly rotating FGK stars. Therefore, in this study, we utilise X-ray observations from Swift and XMM-Newton of ten M dwarf UFRs with P$_{\rm{rot}}$<1 d to determine their X-ray luminosities. Overall, we do not find evidence for supersaturation amongst our UFR M dwarf stars, instead determining them to be at the saturated level, or perhaps even enhanced. Therefore, supersaturation seems not to be the main driver behind the reduced level of flaring activity observed in these stars, and the mystery behind the magnetic activity of UFR low-mass stars remains. Additionally, we provide an updated analysis on the long term variability within our sample using TESS light curves taken during Cycles 5 and 7. We identify 352 optical flares from our sample with energies between $1.2\times10^{31}$ and $8.7\times10^{34}$ erg. We determine flare rates for each TESS cycle, compare them, identifying variations across a 7 year timespan and attribute this to potential activity cycles.

This work establishes oblique shocks in Massive Star Clusters (MSC) as a primary mechanism for accelerating cosmic rays (CR) up to the knee of the energy spectrum. We develop a model that incorporates the combined contribution of supernova and collective wind shocks, emphasizing the critical role of the shock obliquity angle in determining the maximum particle energy. We illustrate, within our model that oblique shocks can significantly enhance acceleration efficiency, allowing particles to reach multi-PeV energies in a rigidity-dependent manner. Our preferred model, which incorporates oblique shocks, reproduces the all-particle spectrum and composition observed by The Large High Altitude Air Shower Observatory (LHAASO), interpreting the knee as arising from a sequence of rigidity-dependent cutoffs. The model also predicts subdominant but detectable gamma-ray and neutrino emissions. This study provides an attempt at building a unified framework connecting MSC particle acceleration to the observed features of the cosmic-ray knee.

We present the discovery of SN 2025mkn, a gravitationally lensed Type II supernova. First detected as a blue transient in ZTF, 0.83$^{\prime\prime}$ from a $z=0.42$ elliptical galaxy, follow-up SNIFS/UH2.2m and LRIS/Keck spectra revealed absorption lines at $z=1.371$. Later JWST NIRCam imaging shows that the bright transient is a close pair of point sources separated by $\sim 0.07^{\prime\prime}$, and a 30 times fainter counterimage opposite the lens, for which NIRSpec reveals strong H$\alpha$ emission also at $z=1.371$. The light curves and spectra are consistent with the Type II supernova source being magnified $\gtrsim 100$ times, with $\sim 250$ required to reconcile its luminosity with that of nearby events such as SN 2023ixf. Lens models are consistent with such high magnifications, and always show that the faint image arrived first (undetected in earlier ZTF imaging), consistent with the later spectral phase of this fainter image. A fourth image is also predicted and possibly detected in the NIRSpec data. Light-curve-based time-delay measurements are not possible due to the first image being the faintest; however, the resolved NIRSpec spectra offer a future opportunity for time-delay cosmography through supernova phase measurements.

The upper atmospheres of temperate sub-Neptunes are strongly influenced by atmospheric dynamics due to their cool equilibrium temperature and thereby longer chemical timescales than the atmospheric dynamical timescales. In this study, we used a three-dimensional (3D) general circulation model to investigate the transport-induced disequilibrium chemistry and vertical mixing on temperate gas-rich mini-Neptunes, using K2-18b as an example. We model K2-18b assuming 180 times solar metallicity and consider it as either a synchronous or an asynchronous rotator, exploring spin-orbit resonances of 2:1, 6:1, and 10:1. We find that the vertical transport affects the chemical structure significantly, making CO$_2$ and CO more abundant ($\sim$10$^{-3}$) in the upper atmosphere compared to the chemical equilibrium abundance (<10$^{-15}$), and horizontal winds further homogenize the chemical composition zonally in this region. Molecular abundances in the photosphere generally agree across different rotation periods. We employ a passive tracer in the model to estimate the one-dimensional (1D) equivalent eddy-diffusion coefficient ($K_{zz}$) of K2-18b, providing a parameter useful for future 1D atmospheric models. Additionally, synthetic transmission spectra generated from our model are compared with the JWST observations, and we find that our model can provide a comparable fit to the observations. This work offers a 3D perspective on transport-induced chemistry on a temperate sub-Neptune and derives vertical mixing parameters to support 1D modelling.

The mass flux of solar and stellar winds is a key quantity for stellar evolution and space weather, yet its physical regulation mechanism remains an unsolved problem. In particular, conventional Alfvén wave--driven models that self-consistently connect the stellar surface to the stellar wind fail to reproduce the observed scaling between stellar X-ray flux and mass-loss rate, a discrepancy that can be largely attributed to the dissipation of a substantial fraction of the wave energy by chromospheric turbulence. To address this issue, we aim to clarify the role of chromospheric turbulence in regulating the stellar wind mass flux. We perform one-dimensional wave-driven wind simulations, comparing cases with and without chromospheric turbulence suppression to assess its impact on coronal and wind properties. We find that suppressing chromospheric turbulence leads to a systematic increase in the coronal particle flux, and hence the wind mass flux, by up to an order of magnitude, particularly in regions of moderately strong magnetic field. This behavior arises from a combination of changes in the Poynting flux at the coronal base and in the asymptotic wind speed. Furthermore, the model with chromospheric turbulence suppression reproduces the observed empirical scaling between coronal magnetic field strength and mass flux without invoking additional energy input mechanisms such as interchange reconnection. These results identify the chromospheric turbulence as a key factor in regulating stellar wind mass flux and highlight the importance of incorporating its effects in models that connect the stellar surface and the stellar wind.

X-ray spectra of accretion-powered X-ray pulsars can often be described using a power-law continuum with a high-energy cutoff, which might be further modified by additional spectral components. The Be X-ray binary system 4U 0115+63 is well known for having one of the highest numbers of detected harmonics of its cyclotron resonant scattering features (CRSFs), a pronounced spectral component known as the ''10 keV feature,'' and quasiperiodic oscillations (QPOs) with a period of about 500 s during outbursts. The changes in count rate by a factor of two during the approximately 500 s QPOs allow us to probe the variation in the spectral components with flux. We study the ''10 keV feature'' in emission, aiming to disentangle it from the broadband continuum and CRSFs and investigate its origin. We focus on the flux-dependent behavior of the CRSF and its harmonics, and particularly the contribution of the ''10 keV feature,'' as seen in the flux-resolved analysis of two NuSTAR observations of the 2015 outburst. Comparing the flux-resolved spectra of a given observation with the respective total dataset revealed a distinct change in overall spectral shape at the position of the ''10 keV feature'' but no comparable deviation at the energies of the harmonic CRSFs. The change associated with the ''10 keV feature'' does not seem to involve its centroid energy, which remains constant within a given observation. We find indications for an anticorrelation between the continuum flux and the ratio of the ''10 keV feature'' flux to the continuum flux within each observation. The analysis strengthens previous claims that the ''10 keV feature'' shows some independence from the remaining features. This result supports the interpretation that the ''10 keV feature'' has a different formation mechanism than the continuum emission, although its origin lies within the same physical environment.

ESA's Hera space mission is on its way to the mission target, the binary asteroid (65803) Didymos. HyperScout-H, one of the instruments onboard Hera, is a hyperspectral imager operating in the visible and near-infrared regions between 0.65 and 0.95 microns. HyperScout-H will enable a detailed assessment of the composition of both objects, Didymos and its satellite Dimorphos, the characterization of space weathering effects, and the possible presence of exogenous material on their surfaces. To monitor instrument functionality, calibration exposures are acquired regularly. This article describes the in-flight calibrations carried out for HyperScout-H during the commissioning and cruise phases. Bias and dark exposures, as well as stellar field observations, were acquired several times after launch. We update the calibration data and monitor instrument performance in the space environment. In five images, the surface of Mars fills the entire field of view, enabling cross-validation of HyperScout-H results with those reported by other Mars missions. The calibration data indicate that the bias pattern is stable, the dark current remains negligible for short exposures, and the detector response is highly linear. We quantify the field-of-view alignment and geometric distortion, and evaluate the point spread function based on the stellar field observations. Stellar observations and Mars swing-by data provide updated radiometric calibration constants, suggesting that in-flight conditions have slightly modified the detector's spectral response. In-flight calibrations are essential to ensure data quality and reliability. The results obtained for HyperScout-H demonstrate that the instrument can achieve its scientific goals in observations of the Didymos-Dimorphos system.

DustPedia and LVL are two samples representative of the local galaxy population, including in total 1011 unique objects of all morphological types, with a wide range of stellar masses ($M_*$) and star-formation activity, and a spectral coverage from the FUV to the FIR. The purpose of this work is to show that these samples cover two complementary ranges in $M_*$ and morphology, making them an ideal set for constraining the dominant processes in the evolution of the galactic dust content. Using the multiwavelength data provided by the two surveys, we fitted the galaxies' spectral energy distribution and estimated their physical properties, in particular the $M_*$, $sM_\mathrm{dust}=M_\mathrm{dust}/M_*$, and sSFR = SFR$/M_*$. By combining DustPedia and LVL, we highlight that the $\log_{10}(sM_{\rm dust})$-$\log_{10}(M_*)$ trend is not monotonic. Thanks to a large number of objects across a wide range of $M_*$, we have been able to fit two smoothly-joined linear correlations: a positive for $\log_{10}(M_{*}/$M$_\odot)\lesssim9.5$ (mainly LVL late spirals and irregulars), and a negative one for larger-mass, mainly DustPedia spirals (early types are distinct and more dispersed in the same mass regime). For $\log_{10}(M_{*}/$M$_\odot)>9.5$, we confirm a strong sM_{\rm dust}-sSFR correlation; dwarf galaxies, instead, lie below this trend, with a large scatter of $sM_{\rm dust}$, for -10.5<$\log_{10}$(sSFR/yr$^{-1}$)<-9.0. By using chemical evolution models we find that the observed $\log_{10}(sM_{\rm dust})$-$\log_{10}(M_{*})$ and $\log_{10}(sM_{\rm dust})$-$\log_{10}$(sSFR) trends can be interpreted mainly by variations in the initial gas mass budget and the galaxy ages, respectively. Low-mass Sm-Irr galaxies with low $sM_{\rm dust}$ and high sSFR can only be reproduced by the models by assuming high photofragmentation rate of large grains, and/or low grain-growth in clouds.

A population of Galactic gamma-ray binaries is currently emerging due to ever increasing sensitivity of gamma-ray observatories. The detection of very high energy (VHE) photons with energies well above 10 TeV from a dozen of sources and the estimated power of those sources make them potentially interesting cosmic ray accelerators. Multi-wavelength observations of gamma-ray binaries revealed that most of them include a young massive star in pair with a relativistic companion, either a black hole or energetic pulsar. Fast stellar winds interacting with powerful relativistic outflows from pulsars or the black hole jets in microquasars are favorable sites for VHE particle acceleration. To estimate the expected number of gamma-ray binaries, we present results of population synthesis calculations of Galactic binaries in which a young massive OB- or Be-star is accompanied by a pulsar capable of producing a powerful relativistic outflow. The distributions over the binary eccentricities, orbital periods, Be-disk inclinations, and the pulsar braking energy losses are taken into account. Conditions for a binary to accelerate VHE particles, radiate and absorb the non-thermal photons that may reach the observer are discussed. We model the anisotropic structure of the zone of interaction of the relativistic pulsar wind with the strongly magnetized massive star's wind. The stellar winds with strong ($\sim$ Gauss) magnetic fields at $\sim$ AU distances colliding with powerful pulsar outflows are capable of accelerating particles up to PeV energies at some orbital configurations and phases. The strong magnetic field in the interaction region produces a highly anisotropic structure of the particle accelerator and emitter in the pulsar outflow. The anisotropic radiation pattern may affect the gamma-ray photon absorption and the number of the observed gamma-ray loud systems.

We review the current state of knowledge of the long-term evolution of the small bodies that give rise to comets and exocomets, as well as their reservoirs. The active cometary phase is only transitory, and bodies that become comets pass from a source population, such as the Kuiper Belt, Oort Cloud or their extra-solar analogues, through the active cometary phase, to eventual dormancy or destruction. We discuss dynamical delivery channels that can move comets from their source reservoirs to orbits with small periapsides, and the depletion of these reservoirs by dynamical and collisional means. We also discuss the physical evolution of cometary nuclei, especially in light of recent advances from missions to Solar System comets such as Rosetta's visit to 67P. We then describe our current knowledge of interstellar objects, which can originate from the same source regions as exocomets but be amenable to detailed study when they enter the Solar System. We include a summary of stellar winds emanating from different types of stars, which become increasingly strong once stars leave the Main Sequence. This is followed by a description of how small bodies are affected by stellar evolution, and the range of comet-like phenomena observed towards white dwarf stars. Overall, while we have an increasingly good picture of the physical and dynamical evolution of Solar System comets, a number of large gaps remain in our knowledge of the physics of exocomets, related to our inability to directly probe these bodies and many of the planets that might be affecting their orbits.

Quasars exhibit stochastic variability across wavelengths, typically well-described by a Damped Random Walk (DRW). However, extreme luminosity changes, known as quasar flares, represent significant departures from this baseline and offer crucial insights into accretion disc dynamics and the fundamental physics of supermassive black hole fueling. While transient surveys have spurred interest in flare detection, a systematic search within the legacy SDSS Stripe 82 dataset -- containing 9,258 confirmed quasars -- has not yet been performed. The primary statistical challenge lies in distinguishing these rare events from ever-present intrinsic noise. To address this, we present FLARE (Flare detection via physics-informed Learning, Anomaly scoring, and Recognition Engine), a generalized three-stage framework for detecting flares present in quasar data. FLARE operates by modeling baseline DRW behavior, applying anomaly scoring to flag potential flares, and utilizing a recognition engine to verify candidates. For Stripe 82, we implement this framework using a physics-informed probabilistic Gated Recurrent Unit (GRU) for baseline modeling, Extreme Value Theory (EVT) for anomaly detection, and benchmarking various open-weight and proprietary Vision Language Models as recognition engines for final verification. Detection is executed on r-band light curves, with candidates cross-checked against g-band data to definitively rule out instrumental artifacts. Applying this pipeline, we successfully identify 27 quasars exhibiting distinct flaring activity.

Merging luminous galaxies are ideal laboratories to study some of the most extreme astrophysical phenomena. The local (z=0.1232) obscured quasar J1356+1026 has two nuclei, North and South (J1356N and J1356S), but despite numerous efforts, J1356S had not yet been confirmed as an AGN. Thanks to the superb sensitivity and spatial resolution of the MIRI/MRS instrument on board the JWST, we present new evidence suggesting that J1356S may indeed host an AGN with log L$_{\rm bol}=43.4\pm^{0.6}_{0.5} erg s^{-1}$. This is supported by the detection of strong coronal line emission at this location and by a spectral shape that differs from that of J1356N and those of the narrow-line region (NLR). Aided by the spatially resolved information of MIRI/MRS and VLT/SINFONI, we also find that the high ionization gas, traced by the coronal lines [Ne V]$14.3~\mu$m and [Si VI]$1.963 \mu$m, has an extension of ~13-15.5 kpc. This is likely a lower limit of the true extension, as suggested by the comparison with optical imaging from HST. {The extended [Ne V] emission can be accounted for by photoionization from the quasar in J1356N in a relatively low density environment, ranging from $\rm n_e\leq 2000-3800 cm^{-3}$ in J1356N and $\rm n_e\leq 600-1200 cm^{-3}$ in J1356S and the NLR, as measured from the [Ne V]$14.3\mu$m and $24.3~\mu$m lines.

It is known that the external irradiation of protoplanetary disks by nearby massive stars can result in mass loss that impacts the disk evolution, however the dynamical impact of external irradiation upon the disk itself has not been explored in detail. We aim to investigate the dynamical effect of asymmetric external irradiation on the structure of such disks. We perform two-dimensional multi-fluid radiation hydrodynamical simulations of protoplanetary disks subject to external irradiation using the PLUTO code, with external irradiation modeled as a plane-parallel flux and a simplified nonaxisymmetric heating rate corresponding to the thermal reemission from hot material within the region marginally optically thick to the external irradiation. We find that a nearby massive star can, under certain conditions, induce significant dynamical effects on a protoplanetary disk, including a shadowed region, pronounced spiral arms in gas, and rings and gaps in dust. The dynamics are caused by the temperature asymmetry driven and maintained by external irradiation, akin to the well-established mechanism of shadow-induced spirals and rings in disk with shadowing from their inner regions. Our results show that if an external temperature asymmetry can be induced it can have a significant dynamical impact on the disk itself (in addition to the well-studied mass loss and truncation effects due to external irradiation), possibly even driving substructure. This prompts further investigation with detailed, dynamical radiative transfer models.

Small carbon hydride cations, such as the methylidyne ion (CH$^+$), play an important role in the chemistry of the interstellar medium (ISM). They participate in gas-phase reaction networks leading to the formation of hydrocarbon species that act as precursors to more complex organic molecules. CH$^+$ is a highly reactive ion that is rapidly destroyed by H, H$_2$, and free electrons, making its excitation challenging to model. Its level populations depend not only on radiative and inelastic processes but also on chemical formation and destruction rates, a mechanism known as chemical pumping. We investigate this effect using a new set of ab initio state-resolved ro-vibrational (reactive and inelastic) collision data to model the observed CH$^+$ emission. Multiple rotational and ro-vibrational transitions of CH$^+$ detected toward the planetary nebula NGC 7027 are analyzed. The chemical structure of CH$^+$ is modeled with the CLOUDY code using updated reaction rates, providing the temperature and density structure across the nebula. A non-local thermodynamic equilibrium (NLTE) analysis is performed using CLOUDY and the single-zone RADEX code with a comprehensive set of spectroscopic and collisional data. In addition, chemical formation and destruction processes are implemented in RADEX and explored via Markov Chain Monte Carlo sampling. The CLOUDY model reproduces the observed CH$^+$ line fluxes within a factor of 1.3 on average. It indicates that rotational and ro-vibrational lines arise from physically distinct regions, primarily differing in temperature. RADEX models show that chemical pumping significantly enhances populations above ($\upsilon = 0, J = 1$), strongly increasing ro-vibrational emission, especially in the $\upsilon =2 \to 1$ band. Single-zone models remain limited, highlighting the need for full 1D modeling including all excitation processes.

We present an analysis of the parallel peculiar velocity correlation function using data from the Cosmicflows-4 (CF4) survey. CF4 significantly extends the depth of the peculiar velocity measurements, mitigating the impact of observers on the cosmic variance. We examine the distribution of cosmic variance using different velocity correlation estimators. The combination of the large peculiar velocity uncertainties and the anisotropy distribution of the CF4 data across the northern and southern hemispheres results in substantial statistical uncertainties in the velocity correlation function. To address this, we test different weighing schemes in the velocity correlation function and implement a more accurate peculiar velocity estimator that reduces velocity uncertainties, consequently decreasing the statistical uncertainty. Using the CF4 group dataset, we derive a growth rate of $f\sigma_8=0.384^{+0.116}_{-0.194}$ and a local growth rate of $f\sigma_8=0.569^{+0.054}_{-0.06}$ through a Markov Chain Monte Carlo method.

The stochastic gravitational wave background (SGWB) is conventionally treated as a passive relic of its astrophysical and cosmological sources, with negligible back-reaction on the matter content of the Universe. Here we show that this assumption needs to be modified once the SGWB and matter are treated as a dynamically coupled non-equilibrium system. Combining linearized general relativity with the fluctuation-dissipation theorem, we derive a generalized Langevin framework that drives the coupled system toward a dynamical equilibrium, which is characterized by a distinctive strain spectrum with a high-frequency cutoff $\mathcal{W}$, and a scale-dependent coupling parameter that screens gravity progressively for the most massive structures. Three findings support this framework. Fitting the equilibrium spectrum to the NANOGrav 15-year dataset yields a Bayes factor of $48\pm 3.8$ over the supermassive black hole binary baseline, achieved entirely within general relativity and the Standard Model. The PTA-calibrated screening mass scale $m_{c}\sim 10^{12}\text{--}10^{14}\,M_{\odot}$ overlaps, with no free cosmological parameter, the $\Lambda$CDM-derived linear-to-nonlinear transition mass $M_{\rm NL}$ of cosmic structure at $\sim 8\,h^{-1}\,\mathrm{Mpc}$. Most strikingly, promoting this concordance to a structural identification expresses $\mathcal{W}$ entirely in terms of $M_{\rm NL}$, and its inverse acquires a transparent physical reading as a coherence threshold for SGWB-matter coupling. $\mathcal{W}$ is thereby a derived quantity linking nanohertz gravitational-wave observables to the late-time cosmological sector. The framework makes distinctive scale-dependent predictions testable by forthcoming large-scale structure surveys and space-borne gravitational-wave observatories.

Understanding the kinematic behaviour of the Small Magellanic Cloud (SMC) remains a challenge addressed by many authors using diverse approaches. Over time, increasing observational evidence has accumulated for tidal perturbations induced by the Large Magellanic Cloud (LMC) on the SMC, especially in its outer regions. In this study, we adopt star clusters as kinematic tracers of the SMC. We analyse 36 clusters distributed across the galaxy's structural regions (Northern Bridge, Southern Bridge, Wing/Bridge, West Halo, Main Body and Counter-Bridge). From each cluster's proper motions, radial velocity and heliocentric distance we estimate Cartesian velocities \((V_x,\,V_y,\,V_z)\) in the SMC reference frame. We also compute the same velocity components under the assumption that the SMC behaves as a rotating disc. We then define the residual velocity \(\Delta V\) for each cluster as the difference between the two velocities derived. Additionally, we perform a kinematic anisotropy analysis to characterise the distribution of kinetic energy across the SMC. We find that increasing values of \(\Delta V\) correlate with increasing cluster distance from the SMC center, and that \(\Delta V \approx 60\ \mathrm{km\,s^{-1}}\) it appears to be a lower limit that separates, in kinematic terms, the areas of tidal origin from those with the best behavior.

Langevin stochastic differential equations provide a dynamical description of pulsar timing noise and gravitational wave background (GWB) signals. They are also central to state space algorithms that have gained traction in pulsar timing array analysis due to their linear computational scaling with the number of observations. In this work, we utilize established methods in diffusion theory to derive analytical solutions (means, covariances, and probability density functions) to Langevin equations relevant to red noise and the GWB signal in pulsars. The solutions give direct physical insight on the dynamics of pulsar timing signals. As a canonical example, we show that the pulsar spin frequency modeled as an Ornstein-Uhlenbeck process is mathematically inconsistent with a stationary GWB signal when the timing residual is the direct observable. The nonstationarity can be partially dealt with by marginalizing over long time deterministic trends in the data. Then, we show that a random process based on an overdamped harmonic oscillator supports both a stationary spin frequency and phase residuals, consistent with a stationary GWB signal. We also turn our attention to a phenomenological model of a neutron star -- a two-component model with spin wandering -- that has been motivated to explain observed timing noise in radio pulsars. We derive analytical expressions for the means, covariances, and cross-covariances of the crust and superfluid rotational states driven by white noise. The associated constant deterministic torques are linked to the quadratic spin-down of pulsars. The solutions reveal the physical origin of nonstationarity in the residual model: the coexistence of damped and diffusive eigenmodes of the system.

Electromagnetic cascades are initiated by TeV gamma rays propagating through the intergalactic medium (IGM), and they can be used to constrain the weak intergalactic magnetic field (IGMF) in cosmic voids. Primary TeV photons produce electrons and positrons through electromagnetic pair production, which can be deflected out of the line-of-sight to the observer by IGMF. In addition, electron-positron pairs can perturb the IGM, triggering plasma instabilities that can cool down the pairs before they upscatter cosmic background photons to GeV energies via inverse Compton (IC) scattering. We investigate the influence of plasma instabilities on the cascade spectrum by introducing a parameterized model for the instability using a publicly available Monte Carlo framework CRPropa. We use extended-emission observations within the field of view of the observer to constrain the IGMF in the presence of plasma instability cooling. Based on spectral observations of the blazar 1ES 0229+200 from Fermi-LAT, we find the best-fit photon spectrum including the plasma instability and IGMF parameters that reproduces the observational data for different observer field-of-view angles and obtain the IGMF constraint in cosmic voids. We find that plasma instabilities with a characteristic length scale of order $10^{2}~\text{kpc}$ reproduce the observed photon spectrum and imply an IGMF strength of order $10^{-17}~\text{G}$.

Scattering phase functions (SPFs) derived from resolved scattered-light images of debris discs are widely used to infer dust grain properties, often via parametric forms such as the Henyey-Greenstein (HG) phase function. However, it remains unclear to what extent the inferred scattering behaviour reflects intrinsic dust properties rather than projection effects, disc geometry, or methodological choices. We test how reliably SPFs and HG asymmetry parameters can be recovered from scattered-light images and identify regimes where geometric and observational effects introduce significant biases. We use a physically motivated forward-modelling framework combining dust-scattering calculations, grain dynamics, and ray-tracing to generate synthetic total-intensity images. Since the intrinsic SPFs are known a priori, phase functions extracted from the images can be directly compared to the input scattering behaviour. We explore a grid of grain size distributions, disc inclinations, and opening angles, and fit two-component HG functions to evaluate how well the forward-scattering parameter $g_{1}$ traces grain properties. Even under idealised conditions with perfect knowledge of disc geometry, the recovered phase functions can differ substantially from the intrinsic SPFs. Limited scattering-angle coverage is the dominant effect: strong forward-scattering peaks at small angles are typically unobservable, leading to non-monotonic trends of apparent anisotropy with grain size. Projection effects, line-of-sight mixing, and SPF-extraction choices further modify the recovered phase functions, causing the fitted $g_{1}$ to depend strongly on viewing geometry and methodology. We conclude that SPFs and HG parameters derived from scattered-light images should be interpreted as effective, observation-dependent quantities rather than direct proxies for dust properties.

Cold Neptunes and sub-Neptunes are among the most common products of planet formation and likely dominate the angular-momentum budgets in most planetary systems, yet their dynamical impact on planetary architectures remains poorly understood. Using N-body simulations, we investigate the evolution of multi-Neptune systems assembled into resonant chains during the gas-disk phase and later coupled to remnant planetesimal disks. We show that planetesimal disks containing $\simeq 1$-$4\%$ of the planetary mass efficiently disrupt resonant chains and trigger global dynamical instabilities on timescales of $1~\mathrm{Myr}$-$1~\mathrm{Gyr}$, providing a pathway for delayed instability long after gas-disk dispersal, albeit with instability timescales that are highly sensitive to disk mass. The ensuing instability drives large-scale orbital rearrangement and loss of planets through collisions, tidal disruption, and ejections. Notably, in most systems at least one planet is scattered inward to $\sim 0.1~\mathrm{au}$ on $\sim 10$-$100$ Myr timescales (for $\sim 5$-$50\; M_\oplus$ planets) following instability onset, with a substantial fraction undergoing tidal capture or disruption. This tidal capture can provide a natural pathway to hot Neptune formation, while compact inner chains, if present, would be destroyed on $\sim 100~\mathrm{Myr}$ timescales by cold sub-Neptunes, naturally explaining the observed decline in the resonant fraction. We argue that the predictions of our model, which yields mass-segregated planets and corresponding relative abundances of cold, wide-orbit, and free-floating planets, can be tested by ongoing and upcoming microlensing surveys.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) promises to drastically accelerate the discovery of ultracool dwarfs (UCDs) over the course of its 10-year survey of the Southern Hemisphere. With the official start of LSST imminent, we showcase LSST's capabilities for discovering and characterizing UCDs using early commissioning data (Data Preview 1). The LSST photometric system at this stage remains poorly understood for faint UCDs. Thus, we begin by cross-matching Data Preview 1 against known UCD catalogs. We recover 1 known UCD from the Ultracool Sheet, 17 UCDs from the Dark Energy Survey, and 17 low mass stars from the Gaia Catalog of Nearby Stars. Using these known UCDs alongside recent spectroscopically-confirmed Euclid objects, we select 89 ultracool dwarf candidates in LSST fields, 17 of which are unique to this work. We present our candidates, a photometric temperature estimate, and discuss lessons learned from using early LSST data products. Finally, we turn to the future and predict potential UCD counts in upcoming LSST commissioning data (Data Preview 2), which is expected to be available to the Rubin community in 2026. Using synthetic populations of brown dwarfs, we forecast over 17,000 objects may be discovered and characterized in Data Preview 2. We predict that several hundred known objects and thousands of as-of-yet undiscovered UCDs may be detected in Data Preview 2 fields.

Weak gravitational lensing is one of the primary cosmological probes, providing powerful constraints on the cosmological model. As Stage IV surveys are expected to deliver data of unprecedented precision, accurate modeling of weak gravitational lensing observables across both linear and non-linear scales becomes increasingly important. In this work, we investigate weak lensing in modified gravity (MG) models, extensions of the standard $\Lambda$CDM cosmology in which gravity deviates from general relativity, generally introducing modifications to the lensing equation. We parametrize these modifications through the common phenomenological function $\Sigma_\mathrm{mg}$ and apply ray-tracing to the density maps of N-body and hydrodynamical simulations. We model the time dependence of $\Sigma_\mathrm{mg}$ analytically, while we introduce a phenomenological scale dependence to represent the screening mechanisms by which MG models reduce to general relativity in high-density environments. Starting from the output of the FLAMINGO hydrodynamical simulations, we generate fully ray-traced convergence maps using our modified lensing model. We analyze how the parameters of our prescription affect the weak lensing convergence power spectrum and compare these effects to other known sources of variation, in particular cosmological parameters and baryonic feedback. We find that the modifications to the lensing equation deriving from the MG model produce non-negligible signatures in the convergence power spectrum and that, within extensions of the $\Lambda$CDM framework, these effects can be larger than those induced by baryonic physics. Our results indicate that modified lensing should become a standard ingredient of the analysis of modified gravity simulations.

We study the atomic Hydrogen (HI) content of galaxies in the core of the Shapley Supercluster (SSC) at <z> ~ 0.048, using observations from the MeerKAT Galaxy Cluster Legacy Survey and optical data from the Shapley Supercluster Survey (ShaSS) project. Our sample comprises 169 galaxies with HI detections in the dynamically active region of Abell 3558 and SC1329-313. Following the literature, we classify galaxies into star-forming main sequence (SFMS), transition (TZ), and red sequence (RS) populations, and examine how the HI content varies across these populations. Galaxies on the SFMS exhibit an average HI gas fraction offset of 0.038 dex from the gas fraction main sequence, while TZ and RS populations show depleted HI fractions of -0.034 and -0.211 dex. HI depletion timescales span from 6 to 170 Gyr (SFMS-TZ-RS) confirming increasingly inefficient star formation with quenching. Scaling relations between HI mass and stellar mass in the SSC are generally consistent with field samples. The most direct signature of the dense environment of the SSC is the marked predominance of TZ galaxies, in contrast to what is observed in the field-dominated sample of xGASS, where the population is mostly composed of SFMS galaxies. Moreover, the SFMS and RS populations have similar size, again in contrast with field populations. These results suggest that galaxies in the SSC are undergoing environmental quenching through starvation or strangulation, rather than rapid gas stripping. Despite detectable HI reservoirs, many galaxies exhibit long depletion times, indicating reduced gas accretion and inefficient star formation.

Recent astronomical observations, in particular from the Kepler and TESS missions and their related follow-ups, have revealed an abundance of exoplanets in the size range between Neptune (4 Earth radii) and Earth (1 Earth radii ), as well as a low occurrence rate of planets around twice the radius of Earth (2 Earth radii). This paper uses statistical methods, in particular, the survival function analysis, to address the known exoplanet population observed mainly from the Kepler's primary mission, in order to mathematically elucidate the orbital distributions (expressed in either the orbital period P or the orbital semi-major axis a), for each of the host stars, in both a collective way, and also separately for the planets grouped into various radius bins. We uncover a log-uniform distribution for the majority of planets except the giants. Based on the results of the statistics, we then visit several possible formation scenarios and pathways for planets in different size ranges, in order to explain the results from a theoretical point-of-view.

As statistical systems, galaxies exhibit a rich interplay between organized structure and stochastic fluctuations across a broad range of spatial scales. This duality motivates the need for quantitative frameworks capable of capturing their morphological complexity. The ordinal patterns framework, along with its associated statistical measures: permutation entropy ($H$), disequilibrium ($D_E$), statistical complexity ($C$), and ordinal network node entropy, has recently emerged as a powerful tool for analyzing such complexity in physical systems. We apply this framework in a multiwavelength, multiscale analysis of the galaxy NGC 628, utilizing observations in the near-ultraviolet, near-infrared, mid-infrared, and millimeter bands. Our results reveal a characteristic spatial scale of approximately 200 parsecs, marking the transition from small-scale structures influenced by star formation and stellar feedback to larger-scale morphology governed by the galaxy's dynamics. Furthermore, we find that the $C$ vs. $H$ trajectories for all wavelengths converge toward a common attractor curve, consistent with the behavior of isotropic Gaussian random fields. This convergence suggests a universal statistical behavior in galactic structure at large scales, despite the differing physical processes traced by each wavelength.

We combine Gaia DR3 5-parameter astrometry with calibrated radial velocities for 23 nearby (<1 kpc) young (<60 Myr) clusters, with membership lists from Cantat-Gaudin et al. (2020). We characterise the plane-of-sky structure of the clusters using Q-Parameter and Angular Dispersion Parameter (ADP) methods. We measure plane-of-sky expansion using several methods. We determine plane-of-sky orientations along which expansion is maximised. We also estimate expansion timescales and traceback ages and compare to isochronal ages. We then look for correlations between cluster properties and discuss sample-wide trends. We find that most young clusters are more smoothly structured in their centers where the rate of dynamical interactions is highest, while hierarchical structure can survive in the sparse outskirts for >10 Myr. We also find that the majority of nearby young clusters exhibit clear signatures of expansion in the plane-of-sky, which in many cases is significantly anisotropic, even at ages >30 Myr. We find evidence that older clusters tend to have directions of maximum expansion oriented closer to parallel with the Galactic plane. The high degree of spatial structure and significant expansion anisotropy imply that the majority of these young clusters have formed with significant spatial and kinematic substructure and not as dense, monolithic clusters. Kinematic ages estimated from expansion timescales and on-sky traceback are generally in good agreement with estimates inferred from stellar evolution models for clusters <10 Myr old. However, many clusters with older isochronal ages appear to have significantly younger kinematic ages. We discuss potential reasons for this discrepancy, including a prolonged embedded and/or gravitationally bound phase in the early stages of the clusters.

More than one hundred years ago, physics has been revolutionized when people realized that electronic orbitals, or electromagnetic interactions in general, are quantized. Now, in this study, we are presenting evidence of quantization of planet orbits around stars. Confining a wave in spatial dimensions "quantizes" its wave number. Therefore, this study points to the evidence of the existence of long-range standing waves in the proto-planetary disks. Such waves, although being on a much larger scale of few tens of AU, have already been found by ALMA observation-so called ring-like structure. Now we see that it may exist within 1 AU, and may exert its effect on the existence and distribution of planets within this distance range to the host star. Careful analysis has been carried out to compare the results of different surveys.

This paper investigates the shadow and polarization images of a Kerr-Sen black hole illuminated by geometrically thick and optically thin accretion disks. We adopt two classes of accretion models, namely the phenomenological radiatively inefficient accretion flow (RIAF) model and the analytical ballistic approximation accretion flow (BAAF) model. Based on radiative transfer theory, we examine the effects of the spin parameter $a$, black hole charge $Q$, and observer inclination angle $\theta$ on the shadow images. Both models show that, as the charge $Q$ increases, the photon rings and the central dark regions shrink simultaneously. Meanwhile, frame dragging gives rise to a pronounced brightness asymmetry, which becomes more significant with increasing $a$ and $\theta$. The main difference between isotropic and anisotropic radiation is that, in the latter case, the higher order images are brighter in the upper and lower polar regions. For the BAAF model, because the conical approximation renders certain regions geometrically thinner, the spatial extent of the higher order images is narrower than that in the RIAF model, and the separation between the direct image and the higher order images is more distinct. In the polarization images, the spatial distribution of the polarization vector directions is mainly determined by gravitational lensing and frame dragging, whereas the intensity near the photon ring and the scale of the higher order images are significantly influenced by $Q$.

Motivated by the recent Dark Energy Spectroscopic Instrument (DESI) DR2 preference for dynamical dark energy, we study interacting dark energy models in which a canonical quintessence field couples to cold dark matter through a field-dependent mass $m(\phi)$. In such scenarios, the effective equation of state inferred under the assumption of non-interacting dark sectors, $w_{\rm eff}(z)$, can differ from the intrinsic scalar-field equation of state $w_\phi(z)$, making an apparent phantom crossing $w_{\rm eff}<-1$ possible without introducing a phantom scalar. We show that a viable realization of this mechanism requires the scalar field to originate from a frozen phase deep in the radiation era, in order for the effective coupling to remain sufficiently suppressed before recombination to evade cosmic microwave background constraints, and for the late-time evolution to become strong enough to reproduce the apparent behavior of $w_{\rm eff}(z)$ preferred by DESI. We identify the general conditions that allow these requirements to be satisfied simultaneously, and present an illustrative phenomenological realization in which $w_{\rm eff}(z)$ evolves from $w_{\rm eff}\approx -1.2$ at $z \approx 1.0$ to $w_{\rm eff}\approx -0.9$ at $z\approx 0.4$. These conditions and requirements serve as a guide for designing future models of this kind which can safely navigate the phantom divide at $w=-1$ in an effective way without phantom fields.

We present a novel, JAX-powered implementation of a parametric component-separation method for CMB polarization data, explicitly designed to handle spatially varying foreground Spectral Energy Distributions (SEDs). The approach models this variation across the sky by grouping sets of pixels that share common foreground spectral parameters, scanning over thousands of such configurations to evaluate the trade-off between model complexity and residual systematic contamination. Built within the FURAX framework -- a JAX-powered environment for CMB data analysis -- our pipeline extends the fgbuster parametric formalism. It enables fully vectorized, GPU-accelerated evaluation of the spectral likelihood, map reconstruction, and diagnostic metrics across tens of thousands of pixel subset configurations, noise realizations, and sky regions. Our implementation achieves up to $\sim 100\times$ speed-up over the scipy TNC optimizer used in fgbuster when running on GPUs, as well as giving more robust results. When applied to LiteBIRD-like simulations with spatially varying foreground SEDs, our optimized K-means configuration reduces the 68% upper limit on the tensor-to-scalar ratio $r$ by $\approx 30\%$ relative to a fixed, previously derived multi-resolution configuration, while maintaining competitive statistical uncertainties.

One of the main questions regarding planet formation is how to cross the metre-scale barrier. Several theories rely on the formation of dust clumps dense enough to collapse under their own gravity. Vortices are promising candidate sites of clump formation because they can concentrate dust 'laminarly' by capturing particles, and 'turbulently' by creating the ideal conditions for the streaming instability. In this two-part series, we assess the validity of both pathways by investigating the effect of backreacting dust on vortices. This first paper focuses on the laminar pathway. We use multiple timescale analysis to create two models of vortex evolution. They differ in their assumptions regarding how much gas crosses the vortex's boundary: the first one assumes that the vortex's mass is constant, whereas the second one assumes that the gas density is constant. These two options epitomize the two ways in which vortices can respond to dust concentration. Essentially, as dust gets closer to the vortex centre, it loses angular momentum. To compensate, the gas must either move away from the vortex centre or change its vorticity (and therefore its shape). This choice neatly emerges from the conservation of a quantity akin to potential vorticity. Interestingly, we find that vortices that adjust their vorticity all evolve towards elliptically unstable shapes. And since the elliptical instability destroys the vortex, we conclude that dust imposes an upper bound on vortex lifetimes. If vortex destruction happens before the dust reaches the Hill density, the 'laminar' vortex pathway to planetesimal formation fails.

One of the main questions in planet formation theory is how to cross the metre-scale barrier. In this two-part series, we assess the merits of vortex-based theories by investigating the effect of backreacting dust on vortices. Specifically, this second paper focuses on the 'turbulent' vortex theory, according to which the streaming instability (SI) might be active in vortices. We re-purpose the dusty vortex models derived in paper I as background flows for a linear stability analysis. To simplify the perturbation equations, we build an analogue of the shearing box that follows vortex streamlines instead of Keplerian orbits. This allows us to study the evolution of small wavelength perturbations. We find that inertial waves and dust density waves can propagate in vortices, but that they are not sinusoidal in time. We then extend resonant drag instability theory to these non-modal waves. This allows us to demonstrate that a close cousin of the SI remains active in vortices, a result that greatly strengthens the case for vortex-induced planetesimal formation. Our results also complement past simulations - which showed that the dust's backreaction makes vortices unstable - by providing insights into the nature of (some of) the unstable modes. The caveat is that our work is restricted to the limit of dilute well-coupled dust and to the linear phase of the instability. Finally, our 'vortex SI' extends to 2D. We explain the mechanism of this 'zonal flow RDI', but remain unsure whether it is the unknown instability seen in 2D vortex simulations.

We analyse the effects of a non-standard evolution of the Universe during the reheating epoch on the spectrum of scalar-induced gravitational waves (SIGWs) accounting for the presence of primordial non-Gaussianity. We show that given values of $w$ and $c_s^2$ leave characteristic features in the spectrum which can be detectable by third generation interferometers like LISA. In addition, we argue that the specific reheating dynamics can suppress or even enhance the spectrum, with crucial consequences for its detectability. We perform a Fisher forecast for different values of $w$ and different scans to assess the detectability of the signal when different values of the amplitude and central frequency are considered.

We study the non-Gaussian features of the gravitational wave (GW) background generated by a population of inspiraling supermassive black hole (SMBH) binaries. We show that the SMBH GW amplitude distribution (GWAD) features a universal heavy power-law tail $\propto A^{-4}$, while the low-amplitude tail depends on the SMBH merger rate and the energy-loss mechanisms of the binaries. The distribution of the induced timing residuals inherits this heavy tail. As a result, the ensemble averaged statistical moments of order three and higher diverge, limiting their usefulness as measures of non-Gaussianity, and the GW background from SMBH binaries exhibits the single loud source principle, according to which the strongest signals are more likely to be caused by a small number of loud sources. We confirm that the variance-averaged Gaussian approximation accurately describes the timing residual statistics. This approximation justifies a factored likelihood structure that combines standard Gaussian-process PTA posteriors with the non-Gaussian population prior, enabling consistent incorporation of non-Gaussian effects into SMBH model inference. We provide a fast and flexible Python implementation to compute the distribution of timing residuals from a given SMBH merger rate or GWAD.

Recent cosmological data reveal tension between parameters inferred from measurements of the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), and supernovae (SN) under $\Lambda$CDM. Typical dynamical dark energy parameterizations (such as $w_0w_a$) that seek to jointly resolve these tensions have an equation of state parameter that crosses into the phantom regime, leading to potential instabilities for physical models. We show that the BAO (early-time) and SN (late-time) sides of the tension can instead be treated independently. Early dark energy (EDE) can reduce the tension between CMB-BAO data by changing the calibration of the sound horizon at the drag epoch $r_d$, with a $\Delta\chi^2 = -{9.4}$ relative to $\Lambda$CDM, raising $H_0$ to 70.87 $\rm km s^{-1}Mpc^{-1}$. EDE alone cannot bring consistency between CMB, BAO, and SN data, but combining with a thawing-quintessence component of dark energy reduces tensions between the three datasets, with $\Delta\chi^2=-12.6$ relative to $\Lambda$CDM without a phantom component, vs. $\Delta\chi^2=-15.8$ for $w_0 w_a$ with one. We consider different SN datasets, using the most recent DES Dovekie catalog as our default while assessing differences with the original DESY5 and Pantheon+ catalogs. While the significance of adding thawing quintessence changes, the EDE solution to the CMB-BAO tension remains nearly unaffected. Moreover, though we do not include direct Hubble constant measurements in these $\Delta\chi^2$ values, the EDE solution reduces the Hubble tension with the Local Distance Network value from $7\sigma$ in $\Lambda$CDM to $2-3\sigma$ depending on the SN dataset, nominally the equivalent of an extra $\Delta\chi^2 \sim -40$ or more.

We report our new lattice results for the sphaleron rate calculated within a thermal effective field theory of soft SU(N) gluons, where $N=2,3$, for a wide range of temperatures spanning from $0.6$-$10^{15}$ GeV at sufficiently large volumes. Comparing these results with sphaleron rates in a non-thermal SU(N) plasma where the infrared gluons are over-occupied, we estimate the typical thermalization time for these ultra-soft gluons during the early stages of reheating after inflation. We have also calculated the non-perturbative thermal axion production rate using lattice techniques which shows significant deviation from its perturbative estimate even at the electroweak scale.

In dense neutrino gases, which exist for instance in supernovae, the flavour states of different neutrinos may become entangled with one another. The theoretical description of such systems may therefore call for simulations on a quantum computer. Existing quantum simulations of simple toy systems are not optimal in the sense that they do not fully exploit the symmetries of the system. Here, we propose a new class of qubit-efficient algorithms based on Dicke states and the $su(2)$ spin algebra. We demonstrate the excellent performance of these algorithms both on classical and on quantum hardware.

We apply a general relativistic framework to static and rotating black hole solutions in Scalar-Tensor-Vector Gravity or modified gravity (MOG). Our results yield exact analytic, closed-form relations expressing the mass $M$, the MOG coupling parameter $\alpha$, and the distance $D$ of a Schwarzschild-MOG black hole in terms of a minimal set of directly measurable elements of the accretion disk: the total frequency shift, the telescope aperture angle, and the redshift rapidity. The resulting expressions are derived for particles close to the midline and line of sight, where the redshift rapidity is treated as a relativistic invariant encoding the evolution of the frequency shift with respect to the emitter's proper time in MOG spacetime. We further extend the formalism to the rotating Kerr-MOG geometry and obtain corresponding relations that determine the rotation parameter $a$ jointly with $M$, $\alpha$, and $D$ on the midline. In the rotating background, we introduced the redshift acceleration (general-relativistic version of jerk) to disentangle the spacetime parameters. Crucially, the explicit appearance of $\alpha$ in these formulas enables direct empirical estimation of this parameter, thereby providing a means to test for departures from standard general relativity. The previous results obtained in the standard Schwarzschild/Kerr backgrounds are recovered in the limit $\alpha \to 0$. The derived expressions are concise and suitable for incorporation into black hole parameter-estimation pipelines.

Data from gravitational-wave (GW) detectors often contains a high rate of non-Gaussian transient noise, known as glitches. The parameters estimated from GW signals coinciding with detector glitches are occasionally biased away from their true values. During the first part of the fourth LIGO-Virgo-KAGRA (LVK) observing run, 29% of GW candidates had overlapping or nearby glitches in one or more detectors. In the latter part of the fourth observation run, sensitivity improvements have increased the rates of GW detection. Consequently, scenarios in which GW signals and detector glitches overlap in time are more likely. In this study, we quantify shifts in inferred posterior distributions for short-duration compact binary coalescence GW signals interacting with common LIGO glitches as a function of time between the signal merger time and the glitch. We find statistically significant biases in parameter estimation for mass, spin, and sky position for "blip", "thunder", and "fast-scattering" glitches. Using these results, we provide estimates of what parameters are most affected by overlapping noise sources, as well as what constitutes a "safe" time separation between a gravitational wave signal and a glitch, without requiring glitch subtraction for unbiased source property estimation. We find that in a majority of cases, all parameters are susceptible to significant bias due to glitch interference. Additionally, we find that glitches that occur within the time prior of the GW signal cause more extreme biases than glitches outside of the time prior.

We investigate measures of distance and redshift in cosmological space-times that admit a shear-free foliation, which we henceforth refer to as `quasi-Newtonian'. Space expands isotropically in this description, and small-scale gravitational physics has a natural Newtonian limit, which makes it ideal for considering the physics of wide classes of cosmological models. By assuming that the energy-momentum tensor is dominated by rest-mass density, and that the 3-velocity of matter is small in the quasi-Newtonian frame, we derive fundamental results for kinematics and light propagation. Our results provide a new way of formulating general-relativistic cosmologies with non-perturbative structures in terms of quantities that can be understood from cosmological perturbation theory and post-Newtonian expansions, and allow us to quantify departures of observables from the predictions of Friedmann cosmology. It thereby provides a route to understanding inherently relativistic space-time structures, such as those that occur in Lemaître-Tolman-Bondi, Szekeres solutions, and Bianchi cosmologies in terms of Newtonian degrees of freedom. We illustrate our results using the degenerate Kasner solution as an example, and explain how our approach can be used to provide new insights into the current cosmological tensions.

We investigate propagation of generic waves on thin planar domain walls effectively described by the scalar DBI model. We pay a particular attention to the possibility of caustic (shock) formation - the process, which may lead to intensive particle emission by domain walls. It is demonstrated that no singularities arise in DBI in 2D flat spacetime in the hyperbolic case, if one starts from smooth initial conditions. Technically, this happens because the same family characteristics of the relevant PDE remain parallel at all the times, albeit not being straight lines generically. Crucially, characteristic curves cease to be parallel beyond the simplified setup of DBI in 2D flat spacetime. In particular, this is shown to be the case in $D>2$ for spherical waves, in an expanding Universe, and in the case of a minimal deformation of DBI necessary for avoiding the domain wall problem in cosmology. However, we prove that DBI remains shock free in the hyperbolic case in all these physically relevant situations. This strongly suggests that caustics can form on planar domain walls only due to the loss of hyperbolicity, and they have a cusp profile. We demonstrate, how the non-trivial structure of DBI characteristics beyond the 2D flat spacetime setup uncovered in this work can significantly affect cusp formation.

The Majoron is a hypothetical (pseudo) Nambu-Goldstone boson arising from the spontaneous breaking of a global lepton number symmetry, and is known as a candidate for dark matter in our Universe. In this paper, we investigate the possibility of probing the Majoron dark matter with a linear optical cavity used in the interferometric gravitational wave detectors. We consider a scenario in which the Majoron dark matter couples to photons through a QED anomaly, leading to an oscillatory photon birefringence induced by the coherent dark matter background. The anomaly coefficient is fixed by requiring the model to simultaneously reproduce the electroweak Higgs scale and a typical right-handed Majorana neutrino mass scale, and the resulting dark matter-photon coupling naturally falls within the sensitivity range of optical interferometers. By incorporating additional optics to extract the birefringence signal, we find that ground-based laser interferometers such as Advanced LIGO, KAGRA, as well as future detectors, can probe a region of the parameter space of Majoron dark matter.

The aim of this chapter, focused on relativistic mean-field models and part of the Encyclopedia of Nuclear Physics, is to provide an introductory, self-contained discussion accessible to a broad audience, including advanced undergraduate students. The chapter surveys the fundamental ideas, assumptions, and theoretical framework underlying relativistic mean-field models, and illustrates their wide range of applications across nuclear science. Particular emphasis is placed on the central role that these models play in the construction of equations of state for strongly interacting matter, as well as on the intimate connections between nuclear experiments, astrophysical observations, and theoretical modeling. In this context, relativistic mean-field theory is shown to provide a unified description of bulk nuclear properties and dense neutron-rich matter, enabling the interpretation of the remarkable structural and observational properties of neutron stars in the emerging era of multi-messenger astronomy.

The cosmological gravitational wave background provides a powerful window on parity-violating physics at energies far beyond the reach of terrestrial experiments. However, any co-located planar detector network is insensitive to isotropic circular polarization, independent of its relative orien- tation. In this letter, we show that this no-go result can be evaded by a new class of co-located 3D interferometer designs, which we call Pyramid, whose non-coplanar configuration geometrically isolates chirality. This new design is a natural extension of the third generation of gravitational wave detectors. The coplanar correlation channel is blind to circular polarization, whereas the co-located non-coplanar channel is insensitive to the unpolarized background and acquires a response only in the presence of nonzero net helicity. Pyramid interferometers therefore furnish a unique probe of cosmological gravitational-wave chirality, opening a realistic terrestrial pathway to test parity violation and fundamental symmetry breaking in the early Universe.