We study the problem of a target search by a Brownian particle subject to stochastic resetting to a pair of sites. The mean search time is minimized by an optimal resetting rate which does not vary smoothly, in contrast with the well-known single site case, but exhibits a discontinuous transition as the position of one resetting site is varied while keeping the initial position of the particle fixed, or vice-versa. The discontinuity vanishes at a "liquid-gas" critical point in position space. This critical point exists provided that the relative weight $m$ of the further site is comprised in the interval $[2.9028...,8.5603...]$. When the initial position follows the resetting point distribution, a discontinuous transition also exists for the optimal rate as the distance between the resetting points is varied, provided that $m$ exceeds the critical value $m_c=6.6008...$ This setup can be mapped onto an intermittent search problem with switching diffusion coefficients and represents a minimal model for the study of distributed resetting.

Gamma-ray bursts are the most luminous electromagnetic events in the universe. Their prompt gamma-ray emission has typical durations between a fraction of a second and several minutes. A rare subset of these events have durations in excess of a thousand seconds, referred to as ultra-long gamma-ray bursts. Here, we report the discovery of the longest gamma-ray burst ever seen with a ~25,000 s gamma-ray duration, GRB 250702B, and characterize this event using data from four instruments in the InterPlanetary Network and the Monitor of All-sky X-ray Image. We find a hard spectrum, subsecond variability, and high total energy, which are only known to arise from ultrarelativistic jets powered by a rapidly-spinning stellar-mass central engine. These properties and the extreme duration are together incompatible with all confirmed gamma-ray burst progenitors and nearly all models in the literature. This burst is naturally explained with the helium merger model, where a field binary ends when a black hole falls into a stripped star and proceeds to consume and explode it from within. Under this paradigm, GRB 250702B adds to the growing evidence that helium stars expand and that some ultra-long GRBs have similar evolutionary pathways as collapsars, stellar-mass gravitational wave sources, and potentially rare types of supernovae.

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel'dovich effect (kSZ). We precisely measure this effect ($S/N\approx 10$) via velocity stacking with more than 800,000 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over $\geq$ 4,000 deg$^2$. We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift, but clear trends with stellar mass and absolute magnitude in $g$, $r$, and $z$ bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5\% confidence, i.e. PTE = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris $z=0.5$, PTE = 0.37) over weaker-feedback cases (IllustrisTNG $z=0.8$, PTE = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around DESI Y1 bright galaxy sample (BGS) and emission-line galaxies (ELG), whose features match qualitative expectations. Finally, we forecast $S/N \sim 50$ for future stacked kSZ measurements using data from ACT, DESI Y3, and Rubin Observatory. These measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

We present the Baryon Acoustic Oscillation (BAO) measurements with the Lyman-alpha (LyA) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the auto-correlation of the LyA forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of two larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped LyA absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at $z_{eff} = 2.33$. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in LyA BAO analysis for the first time. We measure the ratios $D_H(z_{eff})/r_d = 8.632 \pm 0.098 \pm 0.026$ and $D_M(z_{eff})/r_d = 38.99 \pm 0.52 \pm 0.12$, where $D_H = c/H(z)$ is the Hubble distance, $D_M$ is the transverse comoving distance, $r_d$ is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.

Baryon Acoustic Oscillation (BAO) measurements play a key role in ruling out post-recombination solutions to the Hubble tension. However, because the data compression leading to these measurements assumes a fiducial $\Lambda$CDM cosmology, their reliability in testing late-time modifications to $\Lambda$CDM has at times been called into question. We play devil's advocate and posit that fiducial cosmology assumptions do indeed affect BAO measurements in such a way that low-redshift acoustic angular scales (proportional to the Hubble constant $H_0$) are biased low, and test whether such a rescaling can rescue post-recombination solutions. The answer is no. Firstly, strong constraints on the shape of the $z \lesssim 2$ expansion history from unanchored Type Ia Supernovae (SNeIa) prevent large deviations from $\Lambda$CDM. In addition, unless $\Omega_m$ is significantly lower than $0.3$, the rescaled BAO measurements would be in strong tension with geometrical information from the Cosmic Microwave Background. We demonstrate this explicitly on several dark energy (DE) models ($w$CDM, CPL DE, phenomenologically emergent DE, holographic DE, $\Lambda_s$CDM, and the negative cosmological constant model), finding that none can address the Hubble tension once unanchored SNeIa are included. We argue that the $\Lambda_s$CDM sign-switching cosmological constant model possesses interesting features which make it the least unpromising one among those tested. Our results demonstrate that possible fiducial cosmology-induced BAO biases cannot be invoked as loopholes to the Hubble tension "no-go theorem", and highlight the extremely important but so far underappreciated role of unanchored SNeIa in ruling out post-recombination solutions.

The distribution of close-in exoplanets is shaped by the interplay between atmospheric and dynamical processes. The Neptunian Desert, Ridge, and Savanna illustrate the sensitivity of these worlds to such processes, making them ideal to disentangle their roles. Determining how many Neptunes were brought close-in by early disk-driven migration (DDM; maintaining primordial spin-orbit alignment) or late high-eccentricity migration (HEM; generating large misalignments) is essential to understand how much atmosphere they lost. We propose a unified view of the Neptunian landscape to guide its exploration, speculating that the Ridge is a hot spot for evolutionary processes. Low-density Neptunes would mainly undergo DDM, getting fully eroded at shorter periods than the Ridge, while denser Neptunes would be brought to the Ridge and Desert by HEM. We embark on this exploration via ATREIDES, which relies on spectroscopy and photometry of 60 close-in Neptunes, their reduction with robust pipelines, and their interpretation through internal structure, atmospheric, and evolutionary models. We carried out a systematic RM census with VLT/ESPRESSO to measure the distribution of 3D spin-orbit angles, correlate its shape with system properties and thus relate the fraction of aligned-misaligned systems to DDM, HEM, and atmospheric erosion. Our first target, TOI-421c, lies in the Savanna with a neighboring sub-Neptune TOI-421b. We measured their 3D spin-orbit angles (Psib = 57+11-15 deg; Psic = 44.9+4.4-4.1 deg). Together with the eccentricity and possibly large mutual inclination of their orbits, this hints at a chaotic dynamical origin that could result from DDM followed by HEM. ATREIDES will provide the community with a wealth of constraints for formation and evolution models. We welcome collaborations that will contribute to pushing our understanding of the Neptunian landscape forward.

Understanding whether cosmic acceleration arises from a cosmological constant or a dynamical component is a central goal of cosmology, and the Dark Energy Spectroscopic Instrument (DESI) enables stringent tests with high-precision distance measurements. We analyze baryon acoustic oscillation (BAO) measurements from DESI Data Release 1 (DR1) and Data Release 2 (DR2), combined with Type Ia supernovae and a cosmic microwave background (CMB) distance prior. With the larger statistical power and wider redshift coverage of DR2, the preference for dynamical dark energy does not diminish relative to DR1. Using both a shape-function reconstruction and non-parametric approaches with a Horndeski-motivated correlation prior, we find that the dark-energy equation of state $w(z)$ varies with redshift. BAO data alone yield modest constraints, but in combination with independent supernova compilations and the CMB prior they strengthen the evidence for dynamics. Bayesian model comparison shows moderate support for departures from $\Lambda$CDM when multiple degrees of freedom in $w(z)$ are allowed, corresponding to $\approx3\sigma$ tension with $\Lambda$CDM (and higher for some data sets). Despite methodological differences, our results are consistent with companion DESI papers, underscoring the complementarity of approaches. Possible systematics remain under study; forthcoming DESI, \emph{Euclid}, and next-generation CMB data will provide decisive tests.

Recently the Dark Energy Spectroscopic Instrument (DESI) provided constraints on the expansion history from their Data Release 2 (DR2). The DESI baryon acoustic oscillation (BAO) measurements are well described by a flat $\Lambda$CDM model, but the preferred parameters are in mild ($2.3\sigma$) tension with those determined from the cosmic microwave background (CMB). The DESI collaboration has already explored a variety of solutions to this tension relying on variations in the late-time evolution of dark energy. Here we test an alternative -- the introduction of an ``early dark energy'' (EDE) component. We find that EDE models can alleviate the tension, though they lead to differences in other cosmological parameters that have observational implications. Particularly the EDE models that fit the acoustic datasets prefer lower $\Omega_m$, higher $H_0$, $n_s$ and $\sigma_8$ in contrast to the late-time solutions. We discuss the current status and near-future prospects for distinguishing amongst these solutions.

We define strict and lax orthogonal factorization systems on double categories. These consist of an orthogonal factorization system on arrows and one on double cells that are compatible with each other. Our definitions are motivated by several explicit examples, including factorization systems on double categories of spans, relations and bimodules. We then prove monadicity results for orthogonal factorization systems on double categories in order to justify our definitions. For fibrant double categories we discuss the structure of the double orthogonal factorization systems that have a given orthogonal factorization system on the arrows in common. Finally, we study the interaction of orthogonal factorization systems on double categories with double fibrations.

We present new deep imaging and high-resolution spectroscopy of the extreme-abundance-discrepancy planetary nebula Ou 5, together with photoionization modelling aimed at probing its unusual thermal and chemical structure. The nebula exhibits a nested bipolar morphology, including inner and outer shells, faint outer lobes, and polar knots. Remarkably, all these components share a dynamical age of order 10,000 yr. Thermal broadening of the H alpha line relative to heavier ions implies a hydrogen temperature 3000 K to 6000 K, in contrast to the ~ 10,000 K derived from collisionally excited line diagnostics. This provides independent support for the presence of at least two distinct temperature/metallicity phases, as previously proposed to explain extreme abundance discrepancies. Photoionization models with sinusoidally varying metallicity successfully reproduce the observed nebular spectrum and morphology. A mixture of fluctuations with both extreme and moderate metallicity contrasts is required to simultaneously fit the O II and the [O III] observations. The nebular He II emission demands a hotter and more luminous central star than previously inferred, consistent with a ~ 0.58 solar mass post-AGB progenitor evolving toward a CO white dwarf. Ou 5 thus reinforces the link between close-binary nuclei and extreme abundance discrepancies, and provides a valuable benchmark for understanding how common-envelope ejections give rise to the thermal and abundance inhomogeneities observed in planetary nebulae.

Baryon acoustic oscillation data from the first year of the Dark Energy Spectroscopic Instrument (DESI) provide near percent-level precision of cosmic distances in seven bins over the redshift range $z=0.1$-$4.2$. We use this data, together with other distance probes, to constrain the cosmic expansion history using some well-motivated physical classes of dark energy. In particular, we explore three physics-focused behaviors of dark energy from the equation of state and energy density perspectives: the thawing class (matching many simple quintessence potentials), emergent class (where dark energy comes into being recently, as in phase transition models), and mirage class (where phenomenologically the distance to CMB last scattering is close to that from a cosmological constant $\Lambda$ despite dark energy dynamics). All three classes fit the data at least as well as $\Lambda$CDM, and indeed can improve on it by $\Delta\chi^2\approx -5$ to $-17$ for the combination of DESI BAO with CMB and supernova data, while having one more parameter. The mirage class does essentially as well as $w_0w_a$CDM while having one less parameter. These classes of dynamical behaviors highlight worthwhile avenues for further exploration into the nature of dark energy.

We apply the exponential operator method to derive the propagator for a fermion immersed within a rigidly rotating environment with cylindrical geometry. Given that the rotation axis provides a preferred direction, Lorentz symmetry is lost and the general solution is not translationally invariant in the radial coordinate. However, under the approximation that the fermion is completely dragged by the vortical motion, valid for large angular velocities, translation invariance is recovered. The propagator can then be written in momentum space. The result is suited to be used applying ordinary Feynman rules for perturbative calculations in momentum space.

Support points summarize a large dataset through a smaller set of representative points that can be used for data operations, such as Monte Carlo integration, without requiring access to the full dataset. In this sense, support points offer a compact yet informative representation of the original data. We build on this idea to introduce a generative modeling framework based on random weighted support points, where the randomness arises from a weighting scheme inspired by the Dirichlet process and the Bayesian bootstrap. The proposed method generates diverse and interpretable sample sets from a fixed dataset, without relying on probabilistic modeling assumptions or neural network architectures. We present the theoretical formulation of the method and develop an efficient optimization algorithm based on the Convex--Concave Procedure (CCP). Empirical results on the MNIST and CelebA-HQ datasets show that our approach produces high-quality and diverse outputs at a fraction of the computational cost of black-box alternatives such as Generative Adversarial Networks (GANs) or Denoising Diffusion Probabilistic Models (DDPMs). These results suggest that random weighted support points offer a principled, scalable, and interpretable alternative for generative modeling. A key feature is their ability to produce genuinely interpolative samples that preserve underlying data structure.

We investigate the stochastic gravitational wave background (SGWB) generated by primordial black holes (PBHs) in the dense cores of dwarf galaxies (DGs), considering both hierarchical binary black hole (BBH) mergers and close hyperbolic encounters (CHEs). Extending our previous merger framework, we incorporate up to four successive generations of PBHs within a Hubble time and quantify the GW emission from both channels. Our results show that while BBHs dominate the total emission, CHEs occur earlier, provide the first GW signals, and contribute a continuous though subdominant background that becomes relatively more significant once the initial PBH population is depleted and binary formation is suppressed. We compute the resulting SGWB spectra, demonstrating that BBHs and CHEs imprint distinct frequency dependencies consistent with analytical expectations. We then compare the predicted signals with the sensitivity of observatories such as LISA, DECIGO, ET, IPTA, and SKA. The numerical implementation is publicly available at \href{this https URL}.

We report observations of the ultra-high-energy gamma-ray source LHAASO J2108$+$5157, utilizing VERITAS, HAWC, Fermi-LAT, and XMM-Newton. VERITAS has collected $\sim$ 40 hours of data that we used to set ULs to the emission above 200 GeV. The HAWC data, collected over $\sim 2400$ days, reveal emission between 3 and 146 TeV, with a significance of $7.5~\sigma$, favoring an extended source model. The best-fit spectrum measured by HAWC is characterized by a simple power-law with a spectral index of $2.45\pm0.11_{stat}$. Fermi-LAT analysis finds a point source with a very soft spectrum in the LHAASO J2108+5157 region, consistent with the 4FGL-DR3 catalog results. The XMM-Newton analysis yields a null detection of the source in the 2 - 7 keV band. The broadband spectrum can be interpreted as a pulsar and a pulsar wind nebula system, where the GeV gamma-ray emission originates from an unidentified pulsar, and the X-ray and TeV emission is attributed to synchrotron radiation and inverse Compton scattering of electrons accelerated within a pulsar wind nebula. In this leptonic scenario, our X-ray upper limit provides a stringent constraint on the magnetic field, which is $\lesssim 1.5\ \mu$G.

We investigate how differences in the stellar feedback produce disks with different morphologies in Milky Way-like progenitors over 1 $\leq z \leq 5$, using eight state-of-the-art cosmological hydrodynamics simulation codes in the \textit{AGORA} project. In three of the participating codes, a distinct, rotation-dominated inner core emerges with a formation timescale of $\lesssim 300$ Myr, largely driven by a major merger event, while two other codes exhibit similar signs of wet compaction -- gaseous shrinkage into a compact starburst phase -- at earlier epochs. The remaining three codes show only weak evidence of wet compaction. Consequently, we divide the simulated galaxies into two groups: those with strong compaction signatures and those with weaker ones. Galaxies in these two groups differ in size, stellar age gradients, and disk-to-total mass ratios. Specifically, codes with strong wet compaction build their outer disks in an inside-out fashion, leading to negative age gradients, whereas codes with weaker compaction feature flat or positive age gradients caused primarily by outward stellar migration. Although the stellar half-mass radii of these two groups diverge at $z \sim 3$, the inclusion of dust extinction brings their sizes and shapes in mock observations closer to each other and to observed galaxies. We attribute the observed morphological differences primarily to variations in the stellar feedback implementations -- such as delayed cooling timescales, and feedback strengths -- that regulate both the onset and duration of compaction. Overall, our results suggest that disk assembly at high redshifts is highly sensitive to the details of the stellar feedback prescriptions in simulations.

Turbulence in curved spacetimes in general, and in the vicinity of black holes (BHs) in particular, represents a poorly understood phenomenon that is often analysed employing techniques developed for flat spacetimes. We here propose a novel approach to study turbulence in strong gravitational fields that is based on the computation of structure functions on generic manifolds and is thus applicable to arbitrary curved spacetimes. In particular, we introduce, for the first time, a formalism to compute the characteristic properties of turbulence, such as the second-order structure function or the power spectral density, in terms of proper lengths and volumes and not in terms of coordinate lengths and volumes, as customarily done. By applying the new approach to the turbulent rest-mass density field from simulations of magnetised disc accretion onto a Kerr BH, we inspect in a rigorous way turbulence in regions close to the event horizon, but also in the disc, the wind, and in the jet. We demonstrate that the new approach can capture the typical behavior of an inertial-range cascade and that differences up to $40-80\%$ emerge in the vicinity of the event horizon with respect to the standard flat-spacetime approach. While these differences become smaller at larger distances, our study highlights that special care needs to be paid when analysing turbulence in strongly curved spacetimes.