arXiv Daily Digest - 2026-03-13

Genetic analyses indicate minimal gene flow across the so-called Eastern Pacific Barrier (EPB) in larvae of the reef-building coral \emph{Porites lobata}. Notably, Clipperton Atoll, situated on the eastern side of the EPB, is the only site that exhibits detectable genetic connectivity with the Line Islands, which lie to the west of the EPB. To elucidate the relationship between this genetic signal and large-scale Pacific Ocean circulation, we analyze historical trajectories of surface-drifting buoys from the Global Drifter Program (GDP). We first discretize the GDP drifter trajectories into a Markov chain representation and subsequently apply transition path theory (TPT) in combination with Bayesian inference. The TPT analysis identifies reactive trajectories -- pathways that connect the Line Islands to Clipperton Atoll with minimal detours -- whose travel times do not exceed 5 months, which is taken as an upper bound for the larval survival time of \emph{P. lobata}. Consistently, the posterior distribution of transport from Pacific islands west of the EPB to Clipperton Atoll attains a local maximum in the Line Islands at a travel time of approximately 2.5 months. Our probabilistic characterization of the Lagrangian dynamics therefore supports a scenario of weak, but non-negligible, permeability of the EPB, in agreement with the genetic evidence, and it motivates a refined dynamical definition of the EPB based on the remaining duration of reactive trajectories. Furthermore, our results indicate that the connectivity between the Line Islands and Clipperton Atoll is governed primarily by the seasonal modulation of the North Equatorial Countercurrent, rather than by the phase of the El Niño--Southern Oscillation (ENSO). Finally, Clipperton Atoll's role as a terminal sink for trajectories is relevant to the planned mining operations.

Multi-swarm systems, where two or more swarms of mobile agents occupy the same region of space with different parameters and goals, occur in a variety of biological, engineering, and defense applications. Composites of multiple swarms can produce hybrid spatiotemporal patterns, which compared to single swarming systems, are relatively unexplored. In this work, we develop a framework for studying the collision of counter-moving swarms, each with its own preferred, stable velocity before collision. We show that redirection of such swarms after collision occurs when a stable velocity synchronized state of the multi-swarm composite exists. Using a rigid-body approximation, we are able to extract how scatter-redirection transitions scale with swarm parameters in a variety of scenarios from reciprocal and non-reciprocal systems to symmetric and antagonistic parameter values. Our results compare well to simulations of both particle modeled agents and wheeled robots.

As argued by Hone in the paper [Commun. Pure Appl. Math., 74(11):2310--2347, 2021], a ``mismatch" problem remained unsolved while he was investigating continued fraction expansions and Hankel determinants from hyperelliptic curves. In this paper, by introducing a new notion called quadratic orthogonal pairs for hyperelliptic functions, we resolve the corresponding problem. As further applications, we give a thorough treatment of the initial value problems for two discrete integrable systems, i.e. the bilateral Somos-4 and Somos-5 recurrences.

Integrability conditions on local Hamiltonians for one-dimensional quantum systems to be free and interacting fermions are introduced. The definition of free fermion is the simultaneous satisfaction of the Yang-Baxter equation and Shastry's decorated star-triangle relation by the $R$-matrix, which is more general than the previous `free-fermion algebra' by Maassarani and more special than free fermions as in the context of exactly solvable quantum models or integrable classical two-dimensional vertex models dual to quantum spin chains. Free fermionic $R$-matrices are of the difference form and have a conjugation symmetry. These free Hamiltonians may sometimes be deformed by the conjugation operator to describe an integrable interacting system with non-relativistic $R$-matrices, as are the cases of the Hubbard model and the XY model in a longitudinal field. A further criterion is obtain on precisely when such deformations remain integrable. A practical procedure is proposed to iteratively solve the free fermionic $R$-matrices from local Hamiltonians, which can be used to construct non-relativistic $R$-matrices if the conditions are met.

The theoretical efficiency of classical shadow tomography is predicated on a perfect Haar-random unitary ensemble, yet this mathematical ideal remains physically unattainable in near-term hardware. Here, we report the experimental discovery of a fundamental accuracy bound on integrated photonic processors: a ``Hardware Horizon'' where the reconstruction error undergoes a sharp phase transition. While the error initially obeys the predicted statistical scaling $\mathcal{O}(M^{-1/2})$, it abruptly saturates at a floor determined by the spectral distortions of the realized unitary group. By deriving a phenomenological error model, we decouple the competing mechanisms of static coherent spectral distortion and dynamic decoherence, demonstrating that this intrinsic noise floor imposes a hard bound that statistical accumulation cannot overcome. These findings establish that the utility of shadow tomography on NISQ (noisy intermediate-scale quantum) hardware is defined by a specific scaling law involving hardware parameters, necessitating active compensation strategies to bridge the gap between theoretical purity and the noisy reality of integrated photonics.

Living systems are open nonequilibrium systems that continuously exchange energy, matter, and information with their environments, leading to stochastic dynamics with memory and active fluctuations. In this study, we develop a non-Markovian theoretical framework for the entropy dynamics of living systems based on the Keldysh functional formalism and stochastic thermodynamics. The approach naturally incorporates colored environmental noise, memory-dependent dissipation, and many-body interactions, yielding generalized Langevin dynamics and non-Markovian master equations. Within this framework we derive an exact frequency-domain expression for the entropy production rate and show that violations of the fluctuation-dissipation relation provide a direct thermodynamic signature of active biological fluctuations. We further demonstrate that environmental memory enhances low-frequency fluctuations and entropy production, leading to critical slowing down near dynamical instability. These results provide a microscopic physical foundation for the entropy "bathtub" picture of living systems and connect entropy evolution with development, aging, and death in nonequilibrium dynamics.

Machine-learned interatomic potentials (MLIPs) are deployed for high-throughput materials screening without formal reliability guarantees. We show that a single MLIP used as a stability filter misses 93% of density functional theory (DFT)-stable materials (recall 0.07) on a 25,000-material benchmark. Proof-Carrying Materials (PCM) closes this gap through three stages: adversarial falsification across compositional space, bootstrap envelope refinement with 95% confidence intervals, and Lean 4 formal certification. Auditing CHGNet, TensorNet and MACE reveals architecture-specific blind spots with near-zero pairwise error correlations (r <= 0.13; n = 5,000), confirmed by independent Quantum ESPRESSO validation (20/20 converged; median DFT/CHGNet force ratio 12x). A risk model trained on PCM-discovered features predicts failures on unseen materials (AUC-ROC = 0.938 +/- 0.004) and transfers across architectures (cross-MLIP AUC-ROC ~ 0.70; feature importance r = 0.877). In a thermoelectric screening case study, PCM-audited protocols discover 62 additional stable materials missed by single-MLIP screening - a 25% improvement in discovery yield.

In this work, we report coherent perfect absorption (CPA) of anti-modes in an indirectly coupled magnon--polariton system. By examining both single and indirectly coupled cases, we experimentally distinguish the modal decay rate $γ$ from the effective decay rate $γ_{\rm{eff}}$. At CPA, $γ_{\rm{eff}} = 0$, leading to a vanishing output and a visually narrow spectrum in the dB-scale, while the intrinsic linewidth set by $2γ$ remains unchanged, demonstrating that the effective decay rate dictates the spectral amplitude rather than the physical loss. Furthermore, in the indirectly coupled system, CPA persists over a broad, magnetically tunable detuning range, in contrast to the single-detuning CPA observed in the directly coupled case, thereby enabling magnetically reconfigurable and frequency-selective microwave absorbers.

The detection of single particles or molecules represents a critical milestone in the development of biosensing technologies. Recently developed optical sensors based on quasi-bound states in the continuum (qBICs) have primarily focused on detecting global refractive index changes, aiming to simultaneously enhance both refractive index sensitivity and quality ($Q$) factors. However, sensors capable of resolving local refractive index perturbations, such as the binding of a nanometer-sized molecule on a surface, remain elusive and have not yet been demonstrated in BIC geometries due to the limited $Q$ factors and relatively large mode volumes. Here, we demonstrate low-contrast BIC metasurfaces that can perform sensing with a virus-sized single-nanoparticle resolution. The qBIC resonance operating at the critical coupling condition exhibits an experimental $Q$ factor of 4.5 x 10$^4$ in heavy water. The strong interaction between the localized electric field and polystyrene nanoparticles with a diameter of 100 nm enable the experimental observation of step-like resonance wavelength shifts, serving as signatures of individual particle binding events. Furthermore, binding-induced modifications to the qBIC resonance alter the optical confinement and asymmetry factor, inducing changes not only in the resonance wavelength but also in the linewidth and amplitude with single-particle sensitivity. Combined with position-insensitive response and free-space accessible features, low-contrast BIC metasurfaces provide a user-friendly platform for next-generation single-molecule sensing integrated with microfluidic systems.

When resources are scarce, will a population of AI agents coordinate in harmony, or descend into tribal chaos? Diverse decision-making AI from different developers is entering everyday devices -- from phones and medical devices to battlefield drones and cars -- and these AI agents typically compete for finite shared resources such as charging slots, relay bandwidth, and traffic priority. Yet their collective dynamics and hence risks to users and society are poorly understood. Here we study AI-agent populations as the first system of real agents in which four key variables governing collective behaviour can be independently toggled: nature (innate LLM diversity), nurture (individual reinforcement learning), culture (emergent tribe formation), and resource scarcity. We show empirically and mathematically that when resources are scarce, AI model diversity and reinforcement learning increase dangerous system overload, though tribe formation lessens this risk. Meanwhile, some individuals profit handsomely. When resources are abundant, the same ingredients drive overload to near zero, though tribe formation makes the overload slightly worse. The crossover is arithmetical: it is where opposing tribes that form spontaneously first fit inside the available capacity. More sophisticated AI-agent populations are not better: whether their sophistication helps or harms depends entirely on a single number -- the capacity-to-population ratio -- that is knowable before any AI-agent ships.

The `15-minute city' has emerged as a central paradigm in urban planning, promoting universal access to work and essential services within short travel times. Its feasibility-particularly for commuting to work-has however rarely been examined quantitatively. Here, we show that proximity to employment is fundamentally constrained by the internal structure of urban economies. Combining urban geometry with empirically observed firm-size distributions, we derive a lower bound on commuting times that holds independently of planning choices or transport technologies. This bound reveals a sharp transition: when employment is sufficiently concentrated, no spatial rearrangement of workplaces can ensure uniformly short commutes, even under optimal placement. Applied to Paris and its near suburbs, we find that achieving universal 15-minute commutes would require substantial economic restructuring or differentiated mobility strategies. The relevant question is therefore not whether an $x$-minute city is achievable, but what the minimal feasible $x$ is given a city's economic structure and spatial scale.

Prevalent Computational Aberration Correction (CAC) methods are typically tailored to specific optical systems, leading to poor generalization and labor-intensive re-training for new lenses. Developing CAC paradigms capable of generalizing across diverse photographic lenses offers a promising solution to these challenges. However, efforts to achieve such cross-lens universality within consumer photography are still in their early stages due to the lack of a comprehensive benchmark that encompasses a sufficiently wide range of optical aberrations. Furthermore, it remains unclear which specific factors influence existing CAC methods and how these factors affect their performance. In this paper, we present comprehensive experiments and evaluations involving 24 image restoration and CAC algorithms, utilizing our newly proposed UniCAC, a large-scale benchmark for photographic cameras constructed via automatic optical design. The Optical Degradation Evaluator (ODE) is introduced as a novel framework to objectively assess the difficulty of CAC tasks, offering credible quantification of optical aberrations and enabling reliable evaluation. Drawing on our comparative analysis, we identify three key factors -- prior utilization, network architecture, and training strategy -- that most significantly influence CAC performance, and further investigate their respective effects. We believe that our benchmark, dataset, and observations contribute foundational insights to related areas and lay the groundwork for future investigations. Benchmarks, codes, and Zemax files will be available at https://github.com/XiaolongQian/UniCAC.

Wavelength-tunable ultra-intense femtosecond lasers may enable breakthroughs in diverse areas of science spanning attosecond science, particle acceleration and beyond. Conventional crystal-based methods are limited by gain bandwidth and damage thresholds, which restrict their wavelength tunability. Plasma-based frequency conversion, unconstrained by material damage, offers a promising alternative. Here, a novel scheme named Frequency Downshifting Stair (FDS) based on plasma bubble filling control is presented. The FDS enables arbitrary frequency down-conversion of ultra-intense femtosecond pulses and yields chirp-free laser pulses. It can achieve near-100% photon conversion efficiency, approaching the physical limit. This is attributed to the linear control by the FDS of laser chirp evolution during the photon deceleration in the plasma wake bubble. For a laser pulse with an arbitrary wavelength λ_0 (e.g., λ_0=800nm), proof-of-concept PIC simulations demonstrate that a single-stage FDS enables continuous wavelength tuning from λ_0 to {2λ}_0 (800-1600nm). Moreover, a three-stage cascaded FDS achieves more than tenfold frequency (10λ_0) downshifting to a central wavelength of 8.5μm. The FDS scheme thus provides a universal pathway for generating high-energy, few-cycle pulses across the broad infrared regime, offering a powerful new tool for wavelength-dependent ultrafast science.

We study a kinetic multi-agent framework coupling opinion dynamics with epidemic spreading, where individual social behaviour both affects and is affected by disease transmission. Each agent is characterised by an epidemiological state and a continuous opinion variable measuring compliance with non-pharmaceutical interventions. The key mechanism of the model is an asymmetric opinion update driven by epidemic encounters: infection events induce more cautious attitudes, while failed transmissions push individuals toward more extreme opinions. We focus on a prototypical SIS setting, for which we derive a macroscopic kinetic description and, in a fast social-interaction regime, a reduced system of differential equations capturing the feedback between epidemic prevalence and opinion evolution. Convergence of the reduced model is rigorously quantified through a modified Wasserstein distance. Numerical simulations highlight how infection-induced awareness and non-infection-driven extremization jointly shape collective epidemic-opinion dynamics.

Pattern recognition and image classification are essential tasks in machine vision. Autonomous vehicles, for example, require being able to collect the complex information contained in a changing environment and classify it in real time. Here, we experimentally demonstrate image classification at multi-kHz frame rates combining the technique of single pixel imaging (SPI) with a low complexity machine learning model. The use of a microLED-on-CMOS digital light projector for SPI enables ultrafast pattern generation for sub-ms image encoding. We investigate the classification accuracy of our experimental system against the broadly accepted benchmarking task of the MNIST digits classification. We compare the classification performance of two machine learning models: An extreme learning machine (ELM) and a backpropagation trained deep neural network. The complexity of both models is kept low so the overhead added to the inference time is comparable to the image generation time. Crucially, our single pixel image classification approach is based on a spatiotemporal transformation of the information, entirely bypassing the need for image reconstruction. By exploring the performance of our SPI based ELM as binary classifier we demonstrate its potential for efficient anomaly detection in ultrafast imaging scenarios.

Achieving efficient three-dimensional feedback cooling of levitated nanoparticles is a key requirement for precision sensing and quantum control in levitated optomechanics. Here we demonstrate three-dimensional optical feedback cooling of a levitated nanoparticle using an interference-enhanced optical force generated within a single beam path. In this scheme, a weak auxiliary field co-propagates with the trapping tweezer and interferes with it to produce a tunable optical force that enables cold damping along all three center-of-mass motional axes without additional beam paths or trap reconfiguration. Using this approach, we cool a 142-nm-diameter silica nanoparticle in high vacuum to effective temperatures of 625.8, 711.6, and 19.9 mK along the $x$, $y$, and $z$ directions, respectively, at a pressure of $8.5\times10^{-6}$ mbar. The cooling dynamics and their dependence on feedback gain and pressure are well described by a cold-damping model. Because the feedback force is generated optically, the scheme does not rely on electrical actuation and is directly compatible with neutral particles. These results establish interference-based optical forces as a simple and broadly applicable mechanism for three-dimensional feedback control in levitated optomechanics, with a clear pathway toward the quantum regime under improved vacuum and detection conditions.

An anomalous flat-band dispersion provided by a degenerate band edge (DBE) of longitudinal optical modes in a double-grating waveguide is used to enhance second-harmonic generation (SHG). The DBE is a fourth-order exceptional point degeneracy (EPD) in a lossless and gainless waveguide, characterized by the coalescence of four eigenmodes that establish a frozen mode in a cavity. At a DBE resonance, the cavity quality factor scales $Q\propto N^5$, where $N$ is the number of unit cells of the grating waveguide. In our numerical experiments, we observe the peak intensity of the fundamental field in the edge-excited cavity scaling as $I_1\propto N^{3.6}$. This leads to a highly efficient SHG process that is radiated vertically from the cavity (i.e., normal to the grating) without requiring collinear phase matching, with a conversion efficiency scaling as $η\propto N^{8.27}$. These results establish DBE-based waveguides as promising platforms for miniaturized efficient nonlinear photonic devices.

Vertical-cavity surface-emitting laser (VCSEL) structures emitting near 940 nm were grown by solid source molecular beam epitaxy (MBE) on Ge(001) substrates. The VCSEL MBE-growth was realized upon a virtual substrate composed of GaAs on Ge grown by melatorganic vapour phase epitaxy (MOVPE). In situ monitoring during MBE growth employed multispectral reflectometry and magnification-inferred curvature imaging for real-time growth analysis. Curvature measurements revealed progressive compressive stress, while optical reflectivity data confirmed uniform layer growth and accurate stopband formation. Fabricated devices with mesa diameters of 35-40 $μ$m, corresponding to oxide apertures of approximately 11-16 $μ$m, exhibited room-temperature lasing under continuous-wave bias with threshold currents below 3 mA. To the best of our knowledge, this is the first demonstration of monolithically integrated 940 nm VCSELs grown on Ge substrates by MBE. These results confirm the viability of MBE-grown VCSELs on Ge with in situ process control for scalable optoelectronic integration.

Ultrafast lasers with simultaneously high average and peak power have become indispensable for driving a multitude of applications, including high-harmonic generation, strong-field physics, and particle source applications. Both parametric amplifiers and post-compressed Ytterbium lasers have emerged as prime platforms to meet these demands. While multi-pass cell (MPC) based post-compression offers broadband output with high beam quality, it provides limited wavelength tunability and suffers from temporal contrast degradation. Conversely, optical parametric amplifiers (OPAs) provide spectral tunability and high temporal contrast but they are limited by low pump-to-signal conversion efficiency and spatial beam inhomogeneities. Here, we introduce the Optical Parametric Multi-Pass Cell Amplifier (OPMPC), a hybrid architecture that overcomes the limitations of both schemes. Our approach utilizes two non-collinearly intersecting MPCs providing broadband parametric amplification of the seed pulses and complete idler removal after each pass through the crystal, thereby suppressing back-conversion. We experimentally demonstrate a record pump-to-signal power conversion efficiency of 43% using a 1030 nm pump at a 1 kHz repetition rate with a pulse energy of 174 $μ$J. The amplified signal at 1500 nm exhibits excellent beam quality, power and spectral stability and is compressed to 48 fs, demonstrating a new platform for ultrafast pulse generation.

Spin-orbit torque (SOT) enables efficient current-driven control of magnetization, offering a promising pathway toward low-power spintronic devices. However, the origin and propagation of both damping-like (DL) and field-like (FL) SOTs in complex multilayers remain unclear. Here, we investigate NiFe thickness-dependent SOT efficiencies in Ta/Pt/Co/Cu/NiFe/Cu/Capping multilayers (x = 15 nm; Capping = Pt, Al, and SiO2). By employing a spin rotation geometry, the perpendicularly magnetized Pt/Co/Cu stacks serve as a spin source introducing unconventional spin polarization orthogonal to the Oersted field, eliminating its contribution and enabling unambiguous separation of SOTs using planar Hall and polar MOKE measurements. To distinguish bulk and interfacial contributions, we introduce a sample-area-normalized moment m = mNiFe/S, accounting for thickness-dependent magnetization and eliminating uncertainties arising from nominal thickness scaling and magnetic dead layers. We find that DL-SOT follows nearly linear 1/m scaling, consistent with rapid spin absorption at the Cu/NiFe interface but exhibits finite beta_SOT when 1/m approaches zero in both Pt- and Al-capped samples, indicating additional interfacial spin-current contributions at Cu/Pt and Cu/Al interfaces. In contrast, SiO2-capped samples show negligible interfacial contributions. Furthermore, FL-SOT deviates markedly from 1/m scaling and exhibits a significantly longer spin dephasing length (about 1.7 nm) compared to DL-SOT, implying extended propagation across NiFe. Comparative capping-layer studies further corroborate this behavior through interface-dependent spin transport. Our findings clarify the origin and distinct propagation characteristics of DL and FL torques, providing guidelines for engineering interfacial spin-orbit functionalities in ultrathin metallic heterostructures.

In undergraduate optics courses, diffraction gratings are studied extensively, generally within the scalar approximation. When the vector nature of light is taken into account, so-called polarization diffraction gratings have been proposed, which are a cutting-edge research topic due to their numerous applications. This paper proposes a simple experiment to introduce students to polarization diffraction gratings and, at the same time, use this device to apply many of the concepts learned about polarimetry. Although current research uses spatial light modulators and metasurfaces, we use a cheap commercial polarization grating. In addition to show how a polarization grating can be characterized, its use as a cheap and easy-to-use Stokes polarimeter is described and demonstrated experimentally. In performing the experiment, issues typical of inverting linear systems will arise, and this will also provide the opportunity to address the problem of finding well-conditioned systems of equations.

Moving media break time-reversal symmetry and exhibit intriguing optical nonreciprocity. This nonreciprocity is usually weak due to the much lower moving speed of media relative to the speed of light. We demonstrate that strong optical nonreciprocity can emerge in a two-dimensional photonic crystal composed of spinning dielectric cylinders. The photonic crystal supports two types of chiral modes at the Brillouin zone center: hybridized multipole modes and symmetry-protected bound states in the continuum (BICs), both of which carry intrinsic spin angular momentum. For finite wavevectors near the zone center, the BICs transform into quasi-bound states in the continuum (QBICs). Under oblique incidence of circularly polarized plane waves, the photonic crystal exhibits nonreciprocal transmission and absorption that are significantly enhanced at the frequencies of these hybridized multipole modes and QBICs. Furthermore, the high quality factors of the QBICs enable sharp transitions in nonreciprocity. Our work uncovers strong chiral light-matter interactions in periodic moving structures, with potential applications in nonreciprocal light manipulation. The mechanism may also be generalized to other classical wave systems, such as phononic crystals.

Path-length diversity methods may be used for adaptive optics (AO) systems to retrieve phase and amplitude information by measuring intensity across multiple planes. Observations that rely on free-space propagation, such as the nonlinear curvature wavefront sensor (WFS), have been shown to offer excellent sensitivity and robustness to scintillation. However, the default design results in a large opto-mechanical footprint due to unavoidable geometric-optics and wave-optics effects. Measurements recorded in a convergent beam would improve instrument compactness, while concentrating light into smaller detector regions of interest, improving signal-to-noise ratio and possibly wavefront reconstruction speed. In this paper, we study path-length diversity wavefront sensing using four planes of contemporaneous intensity measurements made in a convergent beam. We develop a physical optics propagation model and validate the model by performing wavefront reconstructions in both simulations and lab experiments. The manuscripts core contribution is a practical, intensity-domain, Fourier-transform-based recipe to use a conventional multi-plane Gerchberg-Saxton (or comparable) reconstruction pipeline with convergent-beam measurements, enabling a compact optical layout. We find that this approach offers practical benefits over an equivalent free-space wavefront sensor, in particular reducing size, weight, complexity and cost.

With the continuous expansion of optical networks and the increasing diversity of services, existing operation and maintenance (O&M) approaches are increasingly challenged to meet the rising demands for intelligence and efficiency. Large language models (LLMs), endowed with advanced semantic understanding and contextual analysis capabilities, are emerging as a promising enabler for intelligent optical network O&M. Recent studies have demonstrated the feasibility of applying LLMs to optical network management, marking an important step toward intelligent automation. However, systematic investigations into how LLMs can be effectively integrated into existing O&M workflows remain limited. This paper addresses this gap by drawing inspiration from best practices in real-world O&M workflows and systematically identifying scenarios that are well suited for LLM integration. We highlight that agent-based design is key to improving the executability of tasks, and we propose a multi-Agent collaborative O&M architecture that integrates LLM capabilities with existing O&M tools. The proposed architecture leverages core LLM-related technologies including prompt engineering and tool invocation, to build Agent solutions targeting key tasks such as optical channel management, performance optimization, and fault management. This work presents a conceptual framework for embedding LLM-based Agents into optical network O&M workflows, forming agentized processes that demonstrate the feasibility of LLM-assisted task execution and lay the groundwork for future autonomous O&M systems featuring closed-loop perception, decision-making, and action.

Optical Fabry-Perot cavities are crucial tools for metrology experiments, where they achieve extreme length stability, and for some atomic physics experiments, where tunability to atomic transitions enables atom-light interactions. However, achieving both frequency stability and tunability in a single cavity has remained a challenge, forcing metrology experiments exploiting atom-cavity interactions to rely on external active feedback systems to stabilize the length of the cavity. Here, we describe a piezoelectrically-tunable cavity with a cancellation of the coefficient of thermal expansion at around $5^\circ\mathrm{C}$, achieving fractional frequency instabilities at the $4\times 10^{-13}$ level for 1~s integration time. This advance eliminates the need for external stabilization in many atom-cavity experiments, making this design ideal for applications such as ultra-stable superradiant lasers and other cavity quantum electrodynamics experiments.

Commuting patterns are a central component of urban dynamics and many societal activities. Exogenous shocks, such as a pandemic, might drastically modify them inducing heterogeneous variations across socioeconomic strata. Here, we quantify changes in work commuting patterns in Bogotá, Colombia during three different periods of the COVID-19 pandemic: pre-pandemic (2019), COVID-19 restrictions (2020), and partial reopening (2021). To this end, we use anonymized mobile phone data to infer home and work locations from recurring nighttime and weekday connection patterns, and to build daily commuting metrics. We aggregate mobility flows by administrative boundaries and socioeconomic strata. Additionally, we enrich the dataset with a range of other variables such as territorial vocation (i.e., urban versus rural), demographic information (i.e., population density) and, as a proxy for digital infrastructure quality, geolocated Speedtest measurements from Ookla. We find a marked reduction of commuting during restrictions in 2020 and a strong recovery in 2021, but with persistent heterogeneity across socioeconomic strata. Indeed, while commuting declined similarly across income groups during restrictions, groups of the population in the lower-income bracket rebounded faster to pre-pandemic levels. On the contrary, we find that groups in the higher-income bracket managed to keep higher stay-at-home behavior. Regression analyses reveal that territorial characteristics and disparities in digital connectivity significantly contribute to these differences, suggesting that infrastructure investments could help mitigate mobility-based inequalities.

Aggregation of the main antenna complex of higher plants, Light-Harvesting Complex II (LHCII), is widely used as an in vitro model for energy-dependent quenching (qE), yet fluorescence reduction in aggregates is frequently interpreted without a quantitative separation of intrinsic quenching from excitation-induced annihilation. Here, we address this ambiguity by directly correlating aggregate size, concentration, steady-state fluorescence intensity, and decay kinetics during controlled, incremental aggregation of isolated LHCII. By combining fluorescence correlation spectroscopy (FCS) with TCSPC in a unified experimental framework, we monitored structural and photophysical changes in real time as detergent removal drives biphasic aggregation. We quantified the aggregate composition from the particle concentrations, enabling direct scaling of the absorption cross-section with aggregate size. The average fluorescence lifetime decreased semi-logarithmically with increases in hydrodynamic radius, whereas steady-state fluorescence intensities deviated strongly from this trend. Intensitydependent measurements and steady-state kinetic modeling reveal that singlet-triplet annihilation (STA) emerges at moderate excitation intensities and rapidly becomes the dominant contributor to fluorescence quenching, even for relatively small aggregates. In contrast, intrinsic quenching increases more gradually with aggregate size. By quantitatively disentangling intrinsic excitation quenching from annihilation processes, this work demonstrates that STA can govern the apparent photophysical response of aggregated LHCII across excitation regimes commonly considered non-annihilating. The size-dependent mechanistic framework presented here provides a basis for distinguishing intrinsic quenching from annihilation effects in aggregation-based studies of photosynthetic antenna complexes.

Entropic forces play a fundamental role in nanoscale phenomena, from colloidal self-assembly to biomolecular disaggregation. Here, we develop an exact analytical theory and find general scaling laws for the entropic separation of tether-mediated nanofilament bundles, revealing that a single dimensionless parameter--the ratio of the excluded-volume radius to the tether length--dictates whether filaments are pushed apart or, contrary to the usual expectation, pulled together. This unexpected regime challenges the view that entropic forces invariably promote disaggregation, instead uncovering conditions under which the bundles can remain in attractive metastable states. Brownian dynamics simulations confirm this paradoxical effect, offering predictive insights for applications in biophysics, soft matter physics, and nanotechnology.

We propose a self-consistent simulation model for particle beams in accelerators, which includes the impact of electromagnetic wakefields caused by the geometry of the accelerator chamber. The method is based on a scattered-field formulation for the beam-driven Maxwell's equations. The total electromagnetic field seen by the particles is obtained as the solution of two coupled problems: a purely wakefield problem and a space charge field problem, where for each of these problems, specialized and numerically efficient approaches can be used. To assess the accuracy of the method, we compare simulation results with the analytical solution for a relativistic beam in a uniform accelerator pipe. The numerical efficiency of the method is, furthermore, demonstrated in the beam dynamics study of the multi-cell RF photo-gun installed at the SuperKEK collider facility. We show that electromagnetic wakefields have a non-negligible impact on the quality of the generated beam and, therefore, should be taken into account in the design of high-brilliance electron sources.

The Mpemba effect, in which a system initially farther from equilibrium relaxes faster than a closer one, is often associated with nonlinear or far-from-equilibrium dynamics. We show that this effect can arise entirely within the linear-response regime of many-body systems. In reciprocal systems, a uniform Mpemba effect emerges for three or more degrees of freedom via spectral separation of fast and slow modes. Breaking reciprocity renders the relaxation operator non-normal, enabling a strict componentwise Mpemba effect, with the hotter state relaxing faster even in every individual degree of freedom.

In this paper, single-stage stellarator optimization is combined with stochastic coil optimization to improve the robustness of the stellarator as compared to deterministic methods. The plasma boundary, solved with an MHD solver in fixed-boundary mode, is linked to a set of randomly perturbed coils via the squared flux. The optimizer avoids sharp local minima and can reach improved configurations. Two different configurations obtained with our method, one quasi-axisymmetric and one quasi-helically symmetric, are compared against both the standard stochastic stage II method and the single-stage method. The new configurations shown here yield improved squared flux, quasisymmetry, and particle loss following a posteriori perturbation of the coils.

While current-induced bistate spin-orbit torque (SOT) switching has been well established, deterministic electrical control of multiple magnetic states remains a central challenge in spintronics. Here, we realize a conceptually new multistate SOT device in a SrIrO_3/SrRuO_3 bilayer, hosting four intrinsically stable yet electrically distinguishable magnetic states, including two in-plane canted (IP_c^$\pm$) and two out-of-plane canted (OP_c^$\pm$) states. Pulsed current excitations fully map all twelve deterministic transitions among the four states, establishing a robust switching protocol defined by two characteristic current densities. In-situ scanning nitrogen-vacancy (NV) center magnetometry provides direct real-space evidence for the previously unobserved IP_c^$\pm$ states, and spin dynamics simulations uncover a two-step switching pathway, driven by the concerted action of spin torques and the effective anisotropy field within the fourfold anisotropy landscape. Our demonstration of the intrinsic multistate SOT device directly addresses the density bottleneck of conventional bistate SOT technology, establishing a powerful paradigm for compact, high-speed, and energy-efficient multistate spintronics.

We theoretically investigate the formation of electric and magnetic hot spots with reconfigurable plasmonic antennas. We consider three material systems offering different levels of reconfigurability: gold with the static response, vanadium dioxide which allows for ON/OFF switching, and composite gold-vanadium dioxide material platform which offers a possibility to switch between the electric and magnetic hot spot within a single antenna. Using bowtie and diabolo antennas as a case study, we evaluate optical response functions (scattering and absorption cross-sections, electric and magnetic field enhancement). We demonstrate that the composite material system brings, in addition to enhanced reconfigurability, also novel features of plasmonic antennas, such as strong optical absorption and a joint electric-magnetic hotspot.

Exceptional points (EPs) in non-Hermitian systems offer a remarkably strong response to weak perturbations, but the nonorthogonal nature of the corresponding eigenvectors causes noise to diverge, hindering EPs practical application. Here, we report a twelve-fold enhancement of signal-to-noise ratio (SNR) in magnetic field sensing enabled by a third-order EP of coherent perfect absorption (CPA EP3) in a passive cavity magnonic system. This non-Hermitian magnonic platform comprises two identical yttrium iron garnet (YIG) spheres coherently coupled to a cavity mode, in which the CPA EP3 is realized by engineering the three-mode loss to form a pseudo-Hermitian absorption Hamiltonian. By independently tailoring the absorption EP apart from the resonance EP, the system circumvents the noise divergence caused by eigenbasis collapse. Notably, we harness the sensitivity of the minimum output intensity near CPA to perturbations, yielding a seventyfold SNR improvement and a 400-fold increase in responsivity compared with non-CPA system. A comprehensive noise analysis over one hundred repeated measurements confirms the suppression of frequency noise near the CPA EP3. This demonstrates that our scheme not only avoids the noise divergence plaguing conventional higher-order EP sensors but also provides a general strategy to exploit both CPA and EP for SNR enhancement in passive non-Hermitian systems.

Bragg coherent diffraction imaging (BCDI) fails to reliably retrieve phases in micro-crystals exhibiting strong strain inhomogeneities, which restricts its applicability. Here we show that three-dimensional Bragg ptychography (3DBP) overcomes this limitation by enabling stable inversion for large lattice distortions. Using a combination of experimental measurements and numerical tests, we compare the performance limits of the two approaches and demonstrate that 3DBP tolerates lattice distortions more than six times larger than BCDI. We also establish the sensitivity of both methods on a weakly distorted crystal, for which 3DBP yields smoother amplitude and phase fields with reduced short-length-scale artifacts. 3DBP thus provides a reliable route for imaging micro-crystals with large lattice distortions, expanding the scope of coherent X-ray Bragg microscopy to strongly deformed systems.

Millimeter-wave technologies are essential for future high-speed wireless communications. However, a fundamental challenge remains in the form of severe free-space path loss, where the power density decreases inversely with the square of the distance r (i.e., proportional to r^{-2}) as a spherical dependence. To overcome this limitation, we propose a flexible metasurface tape that is designed to guide electromagnetic energy as surface waves. Unlike conventional free-space propagation, this engineered metasurface confines the field to a subwavelength interface, thereby altering the power decay law to a circular dependence (i.e., proportional to r^{-1}). We numerically and experimentally, for the first time, demonstrate this concept using a periodic grounded-patch array fabricated on a flexible substrate and operated at approximately 100 GHz. The measurement results show that the metasurface tape significantly increases the transmitted power, yielding an average rate of improvement of approximately 40 per meter in received power relative to the free-space baseline in our measurement geometry (e.g., 29-dB increase at 2 m). This increase is realized over a broad bandwidth from 95 GHz to 105 GHz (i.e., approximately 10 %), accommodating wideband modulation schemes required for high-data-rate applications. The flexible, lightweight nature of the tape allows it to be easily installed on diverse surfaces. Our demonstration indicates that the metasurface tape is a promising platform for extending the effective range of millimeter-wave systems, thus offering a robust solution to the path-loss bottleneck in next-generation wireless networks.

A new classical TM stability simulation workflow has been developed that solves the resistive inner-layer equations in a plasma slab to yield a linear, quasi-toroidal TM growth rate $γ$ and mode rotation frequency $ω$. This workflow combines two-fluid and drift MHD effects in a slab approximation of the resistive inner layer (SLAYER) with an effective tearing stability index as $Δ(γ,ω) = Δ' - Δ_\mathrm{crit}$. SLAYER is used to calculate the inner-layer $Δ(γ,ω)$, the STRIDE code is used to calculate a toroidal $Δ'$ that includes shaping effects, and the toroidal $Δ_\mathrm{crit}$ incorporates effects of thermal conduction on Glasser stabilization. This workflow is rapid and numerically robust across reactor-relevant plasma conditions, and yields growth rates that closely align with analytic predictions in well-documented linear growth rate regimes. Using synthetic equilibria, TM stability was calculated across scans of plasma $β$, inverse aspect-ratio, and toroidal current profile gradient. These scans effectively benchmarked this STRIDE+SLAYER workflow against existing models and showed reliable stability predictions in shaped H-mode-like plasmas. This capability to quickly and robustly predict classical tearing stability in tokamaks will facilitate the mapping of TM stable operational regimes and design of safe discharge trajectories in future devices.

Multimode fibers have been proposed for mitigating nonlinear effects in high-power fiber amplifiers, allowing for significant power scaling. Most previous studies on light propagation in continuous-wave fiber amplifiers focus on single mode or few mode fibers. Here we develop a tractable numerical model to simulate light propagation in narrowband, highly multimode fiber amplifiers, which takes into account gain saturation, pump depletion and mode-dependent gain. We consider a frequency domain, field based model, with modal gain being dependent on both intramodal gain and gain-induced mode coupling. We derive coupled equations for the evolution of signal modal amplitudes, pump power and population inversion, and numerically solve these equations using a finite-difference method. For highly multimode excitations, the optical intensity in the fiber is speckled and various modes grow at different rates, due to differential overlap with the gain medium and spatial hole burning. Our analysis is applied to Yb-doped fibers, with a quasi-quantitative analysis of the specific case of Yb, identifying different regimes in which either spontaneous emission (SE) or amplified spontaneous emission (ASE) limit amplifier efficiency, especially for larger core and multimode fibers. Finally, we incorporate ASE and spectrally resolved optical channels into our model and demonstrate the experimentally verifiable phenomenon of ASE suppression with sufficient input signal power. Our model can be combined with existing models for various nonlinear effects, providing a useful tool for quantitatively studying nonlinearity mitigation and power scaling in multimode fiber amplifiers.

Exceptional points (EPs) in non-Hermitian systems are singularities where both eigenvalues and eigenvectors coalesce. In scattering systems, EPs correspond to the merging of scattering states, leading to reflectionless (RL) behavior. A reflectionless exceptional point (RL EP) arises when two RL states further coalesce, yielding an anomalous quartic spectral response. While RL EPs have been explored in bidirectional systems, their unidirectional realization remains elusive. Here, we experimentally demonstrate a unidirectional RL EP by engineering collective states in an anti-Bragg magnonic mirror array. Inversion symmetry is broken using a giant spin ensemble that couples to a waveguide at three spatially separated points, enabling unidirectional reflectionless. At the RL EP, the reflection spectrum flattens and broadens significantly beyond the Lorentzian profile. The observed spectral valleys also expose dark-state behaviors that are typically inaccessible through conventional measurements. Our results provide a route toward controlling collective coherence in open systems and developing broadband unidirectional devices.

Quantifying influence in networks is important across science, economics, and public health, yet widely used centrality measures remain limited: they rely on static representations, heuristic network constructions, and purely endogenous notions of importance, while offering little semantic connection to observable activity. We introduce HawkesRank, a dynamic framework grounded in multivariate Hawkes point processes that models exogenous drivers (intrinsic contributions) and endogenous amplification (self- and cross-excitation). This yields a principled, empirically calibrated, and adaptive importance measure. Classical indices such as Katz centrality and PageRank emerge as mean-field limits of the framework, clarifying both their validity and their limitations. Unlike static averages, HawkesRank measures importance through instantaneous event intensities, enabling prediction, transparent endo-exo decomposition, and adaptability to shocks. Using both simulations and empirical analysis of emotion dynamics in online communication platforms, we show that HawkesRank closely tracks system activity and consistently outperforms static centrality metrics.

The rapid development of photonic quantum information processing necessitates precise and programmable control over optical frequency, a capability critical not only for achieving photon indistinguishability but also for exploiting a virtually unbounded frequency dimension. However, efficient and scalable processing of frequency-encoded photon states remains challenging, primarily due to the limited nonlinear optical interaction in most photonic materials. Here, by harnessing the high-performance thin-film lithium niobate electro-optic (EO) platform, we demonstrate an integrated quantum photonic frequency processor that enables coherent and programmable control of photon frequency with high precision. We establish a scalable architecture for frequency-encoded quantum information processing. Using a fully integrated photonic chip, we realize a universal set of frequency-encoded quantum logic gates, including arbitrary single-qubit rotation gates and the two-qubit controlled-phase gate. Furthermore, we demonstrate its application in high fidelity characterization of frequency-bin entangled states. Our work reveals the unprecedented potential of utilizing the frequency degree of freedom in integrated quantum photonic systems.

We develop a fermi-bose bootstrap embedding (fb-BE) framework for the ground state of interacting elec- trons coupled to phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron-phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for frag- ment size and total system size are demonstrated for one-dimensional Hubbard-Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to infinite system size. Benchmarking against density matrix renormalization group for small 8-site system at half- and quarter-filling shows orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Head-mounted miniaturized microscopes, commonly known as miniscopes, have undergone rapid development and seen widespread adoption over the past two decades, enabling the imaging of neural activity in freely-behaving animals such as rodents, songbirds, and non-human primates. These miniscopes facilitate numerous studies that are not feasible with head-fixed preparations. Recent advancements have enhanced their capabilities, allowing for faster imaging, larger fields of view, and deeper brain penetration. In this review, we examine the latest progress in one-photon and multi-photon miniscopes. We highlight the unique opportunities these devices present for neuroscience research, discuss the current technical challenges, and explore emerging technologies that promise to advance the development of miniscopes.

We propose and experimentally demonstrate a polarization--based steganographic scheme using partially polarized vector beams. In our approach, the spatially dependent polarization structure of the optical field serves as the carrier through which the hidden information can be retrieved. By engineering a vector beam whose polarization states populate a prescribed region of the Poincaré sphere, specifically, the equatorial disk, we establish a nontrivial mapping between transverse spatial coordinates and polarization states. Information retrieval is achieved by applying a spatial mask derived from a parametric curve defined within this region of the Poincaré sphere, followed by spatially resolved polarization analysis. We demonstrate the selective reconstruction of various parametric shapes, including polygonal and smooth curves, confirming that the hidden patterns are retrieved through the combined use of spatial filtering and polarization--domain mapping. Our results establish partially polarized vector beams as a flexible and experimentally accessible platform for polarization--based information hiding.

Diamond is a leading quantum photonics platform due to its ability to host qubits based on crystal defects such as nitrogen vacancy centres. Fabricating nanophotonic devices from defect-rich diamond, which is central to many quantum sensing technologies, promises to enable enhanced performance and integrability of diamond quantum sensors. Here we demonstrate microdisk cavities fabricated from defect-rich diamond that support optical modes with high quality factor ($Q\sim7\times10^4$ at 1042 nm), and show that they exhibit saturable absorption. Power dependent spectroscopy measurements spanning 979 nm to 1604 nm are used to extract wavelength-dependent absorption coefficients and saturation intensities, which indicate that a hydrogen-related defect is a likely origin of the observed absorption. At 1047 nm, we measure a saturation intensity of 3.3 (1) MW/cm$^2$ and an absorption coefficient of 0.537 (4) cm$^{-1}$. These results provide insight into defect-mediated optical loss in diamond nanophotonics and suggest strategies to harness defect-induced nonlinearities in future diamond photonic devices.

Secure communication is essential for modern society, from financial transactions to critical infrastructure. As classical encryption faces threats from advancing computational power, quantum communication provides a fundamentally secure alternative based on physical laws. We present QuaNTUM (Quantum Network at the Technical University of Munich), a modular and extensible quantum communication testbed enabling scalable experiments across fiber-based campus networks and satellite-ground links. The terrestrial network connects research institutions in Garching near Munich via single-mode fibers in a star topology with polarization-maintaining components, multiplexers, and time-synchronized analysis modules. Active polarization control and real-time feedback support stable qubit transmission for high-fidelity quantum key distribution and entanglement distribution. A key feature is the integration of deterministic solid-state single-photon sources, including defects in hexagonal boron nitride and excited erbium atoms, with initial deployments on small satellites to bridge terrestrial and free-space channels. As an open-access platform, QuaNTUM enables protocol development, device benchmarking, and hybrid network research, providing a foundation for scalable quantum communication and future global quantum networks.

Surface-acoustic-wave (SAW) resonators are a promising platform for constructing hybrid quantum systems, where confined acoustic waves enable strong interaction with various physical systems. Focusing SAW resonators, reducing mode volume while suppressing diffraction losses, have recently been investigated for application in such hybrid systems. However, the resonator leads to additional transverse-mode resonances, which cause undesired responses. In this work, we develop single-mode focusing SAW resonators on a thin-film lithium niobate on sapphire. A film thinner than the SAW wavelength allows a highly confined acoustic-wave mode to be localized on the substrate surface. By using contoured electrodes following a two-dimensional Gaussian beam shape, we make the SAW mode focused to nearly a diffraction-limited and confirm it via optical imaging. The apodization technique applied to the interdigitated transducer electrodes suppresses the excitation of higher-order transverse modes.

Cubic boron arsenide (BAs) is a promising semiconductor for next-generation electronics due to its outstanding ambipolar mobility and thermal conductivity, the latter of which is attributed to the suppression of three-phonon scattering. However, precisely accounting for different high-order anharmonic scattering processes is challenging from both theory and experiment, so that questions remain open regarding the ultimate limit of phonon lifetime and thermal conductivity in BAs. Here we show that this gap nearly eliminates three-phonon scattering for zone-center optical phonons in a wide temperature range, leading to a record-high, isotope purity-limited phonon coherence with a quality factor above $3.7\times 10^3$ for >98% enriched $^{11}$BAs below 100 K. We discriminate three decoherence mechanisms by their temperature-dependent contribution to the damping rate using high-resolution Raman and Fourier transform infrared spectroscopy. For the as-synthesized crystals, we find that defect scattering has negligible contributions to the linewidth of optical phonons in comparison to isotope scattering. These results provide critical insights into the intrinsic and extrinsic scattering mechanisms of optical phonons in BAs, motivating further studies to quantify anharmonic effects and realize superior phonon transport.

Viscosity is a fundamental property of liquid lubricants, yet it is challenging to determine accurately, especially at high pressures. Although equilibrium molecular dynamics (EMD) simulations are a promising alternative to resource-intensive experiments, practical challenges remain in assessing the sufficiency of simulation time and in controlling uncertainties in the Green-Kubo formalism due to the finite amount of trajectory data. In this work, we extend the STable AutoCorrelation Integral Estimator (STACIE), a recently developed algorithm for estimating transport properties. First, we introduce the Lorentz model to estimate the viscosity and the exponential correlation time from the low-frequency power spectrum of anisotropic pressure fluctuations. Second, we show how to supplement the three conventional off-diagonal elements of the pressure tensor ($P_{xy}$, $P_{yz}$ and $P_{zx}$) with two additional independent time series for shear viscosity calculations. Using these improvements, we apply STACIE to calculate the shear viscosity of 2,2,4-trimethylhexane from EMD simulations. We demonstrate STACIE's capability to reliably calculate viscosity under high-pressure conditions, offering a robust and automated solution with validated uncertainty quantification. Our results, when compared to the outcomes of the 10th International Fluid Properties Simulation Challenge (IFPSC), underscore the need for long EMD simulations. Large deviations from experimental viscosities in previous works were primarily due to insufficient simulation times and ad hoc post-processing choices, rather than the limitations of the force fields used. Unlike previous studies, our viscosity estimates agree well with experimental results (relative error < 6%) up to the highest pressure of 1 GPa, highlighting the improved reliability and accuracy of STACIE's systematic approach to viscosity predictions.

This article frames the relation between biological and physics by characterizing the former as a subdiscipline rather than a special case of the latter. To do this, we posit biological physics as the science of living matter in contrast to classic biophysics, the study of organismal properties by physical techniques. At the scale of the individual cell, living matter is nonunitary, i.e., not composed of aggregated subunits, and has features (e.g., intracellular organizational arrangements and biomolecular condensates) that are unlike any materials of the nonliving world. In transiently or constitutively multicellular forms (social microorganisms, animals, plants), living matter sustains physical processes that are generic (shared with nonliving matter, e.g., subunit communication by molecular diffusion in cellular slime molds), biogeneric (analogous to nonliving matter but realized through cellular activities, e.g., subunit demixing in animal embryos) or nongeneric (pertaining to sui generis materials, e.g., budding of active solids in plants). This "forms of matter" perspective is philosophically situated in the dialectical materialism of Engels and Hessen and the multilevel physicalism of Neurath and the logical empiricists. We counterpose this view to informationism and to genetic and other hierarchically reductionist physical theories of biological systems and highlight open questions regarding incompletely characterized and enigmatic forms of living matter.

Differentiable programming, enabled by automatic differentiation (AD), provides a robust framework for gradient-based optimization in computational plasma physics. While optimization is often only used towards design, we demonstrate that it can also be used for discovery and bridging the gap towards multi-scale modeling. We discuss four applications: (1) discovering novel nonlinear plasma phenomena, including a previously unknown superadditive wavepacket interaction regime, by optimizing differentiable kinetic simulations; (2) learning hidden variables that capture spatiotemporally non-local kinetic effects in fluid simulations, enabling hydrodynamic models to reproduce large Knudsen number physics typically requiring kinetic solvers; (3) accelerating Thomson scattering analysis by over $140\times$ while enabling extraction of velocity distribution functions with $\mathcal{O}(10^3)$ parameters; and (4) inverse design of spatiotemporal laser pulses that achieve target far-field behavior where full space-time coupling improves performance by $15\times$ over spatial or temporal optimization alone. These examples illustrate that differentiable programming not only accelerates existing design and inference workflows but enables qualitatively new capabilities, from algorithmic physics discovery to high-dimensional inference and design previously considered intractable.

Reaction-diffusion systems driven far from thermodynamic equilibrium through the injection of energy can support multiple distinct spatial patterns that persist as long-lived dynamical phases. The stability of these metastable phases is not determined by thermodynamics, but by the transition paths connecting them. At finite particle numbers, intrinsic stochasticity induces rare transitions between competing patterns, rendering continuum mean-field descriptions insufficient, while exact stochastic simulations become computationally prohibitive in spatially extended systems. Here, we develop a nonequilibrium instanton framework that enables efficient computation of transition rates between metastable patterns from a single optimal transition path and its fluctuations. Using this theoretical framework, we show that an effective entropy in path space can qualitatively alter stability at finite particle numbers by increasing the exit rates of metastable patterns. By studying models of varying complexity, this work establishes path entropy as a key organizing principle for nonequilibrium pattern formation.

The cause of electron transfer in contact electrification is one of the most hotly debated physical problems today. In this study, the electron transfer is hypothesized to be partly driven by the surface dipole induced potential during contact. This phenomena is demonstrated by a combination of atomistic field theory (AFT) and molecular dynamics (MD) simulation. A representative contact system of carbon and silicon dioxide was chosen for its excellent tribo-tunneling power output performance. The results reveal the existence of a nonlinear potential field as well as the existence of a separation dependent potential barrier at the contact interface. Possible scenarios of triboelectric charge transfer are discussed in light of these results. These results are critical to the fundamental understanding of contact electrification.

For a long time, people have been focusing on how to extract more information, such as off-diagonal observables, from the quantum Monte Carlo (QMC) simulation of the partition function, but there have been numerous difficulties, and many of them are insurmountable. In this article, we point out that all the difficulties stem from the starting point of the simulation: calculating a partition function. We introduce a paradigm shift: when we transform the simulated object from a partition function to a generalized reduced density matrix (GRDM), the difficult problem of measurement can be readily solved. By designing the GRDM, both equal-time and nonequal-time off-diagonal observables have been measured easily in QMC with a polynomial computation complexity. As a demonstration, the GRDM enables direct access to nonequal-time correlators for dynamical spectra as well as Rényi-1 correlators that reveal strong-to-weak symmetry breaking in the mixed state, capabilities that lie beyond the reach of prior methods. This establishes a unified framework for holographic characterization within QMC.

Owing to their localized field enhancement and subwavelength resolution, surface waves (SWs) offer broad application potential in communications, sensing, and photonics via on-chip wavefront manipulation. This makes multi-channel SW wavefront manipulation highly desirable. However, conventional metasurfaces for SW wavefront shaping, relying on geometric and propagation phase mechanisms, typically exhibit similar functionalities for co- or cross-polarized output channels under different circularly polarized (CP) incidences, thereby limiting the development of high-capacity on-chip integrated devices. Here, by introducing the chirality-assisted phase as an additional phase control mechanism, we effectively decouple both co- and cross-polarized output channels, enabling independent SW wavefront shaping in four distinct channels. We numerically and experimentally demonstrate two metasurfaces in the microwave range: a four-channel SW meta-deflector and a four-channel SW metadevice that simultaneously produces a focused SW beam, a SW Bessel beam, and two deflected SW beams in different directions. Therefore, chirality-assisted geometric-phase metasurfaces provide a versatile platform for multi-channel SW wavefront engineering, offering significant potential for high-capacity on-chip communication and integrated photonic systems.

The generation of arbitrary polarization states of light is essential for optical communication and photonic information processing. Photonic crystal and metasurface platforms supporting bound states in the continuum (BICs) provide a powerful route for polarization engineering through tailoring the radiation from the resonant modes. However, existing approaches typically rely on static structural symmetry breaking or off-normal radiation, which limits continuous polarization tuning of vertical radiation. Here, we demonstrate a magnetooptical metasurface that generates arbitrary polarization states of light at normal radiation. By applying an external magnetic field with variable rientation, a symmetry-protected BIC is transformed into a quasi-BIC whose radiation polarization can be continuously tuned. The magneto-optical perturbation drives the controlled migration of polarization singularities in momentum space, allowing the emitted states to continuously span the entire Poincaré sphere without structural modification. This approach establishes a compact platform for actively tunable polarization sources and polarizationencoded photonic devices.

Brains remain unrivaled in their ability to recognize and generate complex spatiotemporal patterns. While AI is able to reproduce some of these capabilities, deep learning algorithms remain largely at odds with our current understanding of brain circuitry and dynamics. This is prominently the case for backpropagation through time (BPTT), the go-to algorithm for learning complex temporal dependencies. In this work we propose a general formalism to approximate BPTT in a controlled, biologically plausible manner. Our approach builds on, unifies and extends several previous approaches to local, time-continuous, phase-free spatiotemporal credit assignment based on principles of energy conservation and extremal action. Our starting point is a prospective energy function of neuronal states, from which we calculate real-time error dynamics for time-continuous neuronal networks. In the general case, this provides a simple and straightforward derivation of the adjoint method result for neuronal networks, the time-continuous equivalent to BPTT. With a few modifications, we can turn this into a fully local (in space and time) set of equations for neuron and synapse dynamics. Our theory provides a rigorous framework for spatiotemporal deep learning in the brain, while simultaneously suggesting a blueprint for physical circuits capable of carrying out these computations. These results reframe and extend the recently proposed Generalized Latent Equilibrium (GLE) model.

Universal Machine Learning Interatomic Potentials (uMLIPs) enable atomistic simulations and high-throughput screening at scales far beyond those accessible with density functional theory (DFT). However, most existing uMLIPs are trained on Perdew--Burke--Ernzerhof (PBE) generalized gradient approximation (GGA) data and are therefore fundamentally limited by PBE-level accuracy. In this paper, we argue that better zero-shot predictions versus experiments must be an explicit design target for uMLIPs and present PFP v8, a uMLIP available on the Matlantis service that overcomes the inherent limitations of the PBE functional by being trained to reproduce the regularized-restored strongly constrained and appropriately normed (r2SCAN) meta-GGA potential-energy surface across a wide range of chemical domains. Without requiring domain-specific fine-tuning, PFP v8 delivers systematically improved agreement with experimental data or high-accuracy references for crystals, molecules, and surfaces, outperforming PBE-based DFT calculations. Crucially, in long-time molecular dynamics simulations that are computationally impractical with DFT, PFP v8 predicts melting points with an average error of approximately 130 K, halving the error relative to PBE-trained models. These results establish that uMLIPs can move beyond the limitations of their training approximations and achieve substantially improved agreement with experiment across diverse chemical domains, further narrowing the gap between simulation and reality.

Machine Learning Interatomic Potentials (MLIPs) enable accurate large-scale atomistic simulations, yet improving their expressive capacity efficiently remains challenging. Here we systematically develop Mixture-of-Experts (MoE) and Mixture-of-Linear-Experts (MoLE) architectures for MLIPs and analyze the effects of routing strategies and expert designs. We show that sparse activation combined with shared experts yields substantial performance gains, and that nonlinear MoE formulations outperform MoLE when shared experts are present, underscoring the importance of nonlinear expert specialization. Furthermore, element-wise routing consistently surpasses configuration-level routing, while global MoE routing often leads to numerical instability. The resulting element-wise MoE model achieves state-of-the-art accuracy across the OMol25, OMat24, and OC20M benchmarks. Analysis of routing patterns reveals chemically interpretable expert specialization aligned with periodic-table trends, indicating that the model effectively captures element-specific chemical characteristics for precise interatomic modeling.

Absolute Binding Free Energy (ABFE) methods are among the most accurate computational techniques for predicting protein-ligand binding affinities, but their utility is limited by the need for many simulations of alchemically modified intermediate states. We propose Direct Binding Free Energy (DBFE), an end-state ABFE method in implicit solvent that requires no alchemical intermediates. DBFE outperforms OBC2 double decoupling on a host-guest benchmark and performs comparably to OBC2 MM/GBSA on a protein-ligand benchmark. Since receptor and ligand simulations can be precomputed and amortized across compounds, DBFE requires only one complex simulation per ligand compared to the many lambda windows needed for double decoupling, making it a promising candidate for virtual screening workflows. We publicly release the code for this method at https://github.com/molecularmodelinglab/dbfe.

The mathematical characterization of social-distancing games in classical epidemic theory remains an important question, for their applications to both infectious-disease theory and memetic theory. We consider a special case of the dynamic finite-duration SI social-distancing game where payoffs are accounted using Markov decision theory with zero-discounting, while distancing is constrained by threshold-linear running-costs, and the running-cost of perfect-distancing is finite. In this special case, we are able construct strategic equilibria satisfying the Nash best-response condition explicitly by integration. Our constructions are obtained using a new change of variables which simplifies the geometry and analysis.As it turns out, there are no singular solutions, and a time-dependent bang-bang strategy consisting of a wait-and-see phase followed by a lock-down phase is always the unique strategic equilibrium. We also show that in a restricted strategy space the bang-bang Nash equilibrium is an ESS, and that the optimal public policy exactly corresponds with the equilibrium strategy.

Modern imaging instruments can produce terabytes to petabytes of data for a single experiment. The biggest barrier to processing big image datasets has been computational, where image analysis algorithms often lack the efficiency needed to process such large datasets or make tradeoffs in robustness and accuracy. Deep learning algorithms have vastly improved the accuracy of the first step in an analysis workflow (region segmentation), but the expansion of domain specific feature extraction libraries across scientific disciplines has made it difficult to compare the performance and accuracy of extracted features. To address these needs, we developed a novel feature extraction library called Nyxus. Nyxus is designed from the ground up for scalable out-of-core feature extraction for 2D and 3D image data and rigorously tested against established standards. The comprehensive feature set of Nyxus covers multiple biomedical domains including radiomics and cellular analysis, and is designed for computational scalability across CPUs and GPUs. Nyxus has been packaged to be accessible to users of various skill sets and needs: as a Python package for code developers, a command line tool, as a Napari plugin for low to no-code users or users that want to visualize results, and as an Open Container Initiative (OCI) compliant container that can be used in cloud or super-computing workflows aimed at processing large data sets. Further, Nyxus enables a new methodological approach to feature extraction allowing for programmatic tuning of many features sets for optimal computational efficiency or coverage for use in novel machine learning and deep learning applications.

DNA adder circuits are programmable reaction networks that process DNA molecular inputs to compute a sum and serve as essential components for digital computation. Currently, DNA adders primarily focus on binary addition. While efforts extend the operational bit-width by minimizing the number of DNA strands and developing carry-transmission mechanisms, challenges such as the susceptibility of carrying information to attenuation and the limited expressive capacity of the binary system impose significant constraints on computational scale. This paper proposes a scalable ternary adder architecture by introducing an innovative competitive blocking (CB) circuit. The architecture employs a dual cooperative optimization strategy that significantly enhances single-bit computational capacity and incorporates a dynamic concentration adjustment (CA) to effectively broaden the computational bit-width. Consequently, a significant increase in molecular computing scale is achieved compared to previous binary adders. Biochemical experimental results indicate that the CB circuit effectively outputs the ternary full-adder bit and successfully performs 10-bit addition. Furthermore, by implementing the CA strategy, this adder can be further extended to support 17-bit addition. This research provides a novel methodological foundation for advancing DNA computing technologies and offers promising potential for scalable digital computing applications.

Representations learned by convolutional neural networks (CNNs) exhibit a remarkable resemblance to information processing patterns observed in the primate visual system on large neuroimaging datasets collected under diverse, naturalistic visual stimulation, but with instruction for participants to maintain central fixation. This viewing condition, however, diverges significantly from ecologically valid visual behaviour, suppresses activity in visually active regions, and imposes substantial cognitive load on the viewing task. We present a modification of the encoding model framework, adapting it for use with naturalistic vision datasets acquired under fully natural viewing conditions, without fixation, by incorporating eye-tracking data. Our gaze-aware encoding models were trained on the StudyForrest dataset, which features task-free naturalistic movie viewing. By combining eye-tracking data with the visual content of movie frames, we generate combined subject-wise gaze-stimulus specific feature time series. These time series are constructed by sampling only the locally and temporally relevant elements of the CNN feature map for each fixation. Our results demonstrate that gaze-aware encoding models match the performance of conventional encoding models with 112x fewer model parameters. Gaze-aware encoding models were especially beneficial for participants with more dynamic eye-movement patterns. Therefore, this approach opens the door to more ecologically valid models that can be built in more naturalistic settings, such as playing games or navigating virtual environments.

Phytolith analysis is a crucial tool for reconstructing past vegetation and human activities, but traditional methods are severely limited by labour-intensive, time-consuming manual microscopy. To address this bottleneck, we present Sorometry: a comprehensive end-to-end artificial intelligence pipeline for the high-throughput digitisation, inference, and interpretation of phytoliths. Our workflow processes z-stacked optical microscope scans to automatically generate synchronised 2D orthoimages and 3D point clouds of individual microscopic particles. We developed a multimodal fusion model that combines ConvNeXt for 2D image analysis and PointNet++ for 3D point cloud analysis, supported by a graphical user interface for expert annotation and review. Tested on reference collections and archaeological samples from the Bolivian Amazon, our fusion model achieved a global classification accuracy of 77.9\% across 24 diagnostic morphotypes and 84.5% for segmentation quality. Crucially, the integration of 3D data proved essential for distinguishing complex morphotypes (such as grass silica short cell phytoliths) whose diagnostic features are often obscured by their orientation in 2D projections. Beyond individual object classification, Sorometry incorporates Bayesian finite mixture modelling to predict overall plant source contributions at the assemblage level, successfully identifying specific plants like maize and palms in complex mixed samples. This integrated platform transforms phytolith research into an "omics"-scale discipline, dramatically expanding analytical capacity, standardising expert judgements, and enabling reproducible, population-level characterisations of archaeological and paleoecological assemblages.

We introduce a subclass of Lie symmetries, called parameter-state symmetries, to analyse the local structural identifiability and observability of mechanistic models consisting of state-dependent ODEs with observed outputs. These symmetries act on parameters and states while preserving observed outputs at every time point. We prove that locally structurally identifiable parameter combinations and locally structurally observable states correspond to universal invariants of all parameter-state symmetries of a given model. We illustrate the framework on four previously studied mechanistic models, confirming known identifiability results and revealing novel insights into which states are observable, providing a unified symmetry-based approach for analysing structural properties of dynamical systems.

Antimicrobial resistance (AMR) poses a global health threat, reducing the effectiveness of antibiotics and complicating clinical decision-making. To address this challenge, we introduce abx_amr_simulator, a Python-based simulation package designed to model antibiotic prescribing and AMR dynamics within a controlled, reinforcement learning (RL)-compatible environment. The simulator allows users to specify patient populations, antibiotic-specific AMR response curves, and reward functions that balance immedi- ate clinical benefit against long-term resistance management. Key features include a modular design for configuring patient attributes, antibiotic resistance dynamics modeled via a leaky-balloon abstraction, and tools to explore partial observability through noise, bias, and delay in observations. The package is compatible with the Gymnasium RL API, enabling users to train and test RL agents under diverse clinical scenarios. From an ML perspective, the package provides a configurable benchmark environment for sequential decision-making under uncertainty, including partial observability induced by noisy, biased, and delayed observations. By providing a customizable and extensible framework, abx_amr_simulator offers a valuable tool for studying AMR dynamics and optimizing antibiotic stewardship strategies under realistic uncertainty.

The study of navigation behaviour and the associated brain dynamics have been a focus increasing research over the last decades. Coinciding with this has been an increased focus on a more ecological understanding of cognition. Here we review recent research seeking to provide a more naturalistic, ecological understanding of human navigation behaviour and brain dynamics. Research in this area falls into four categories: testing navigation in real-world environments, analysis of data collected from tracking individuals during daily life, navigation in simulated or virtual environments mimicking the real-world, and mobile brain recording methods. Combining these different approaches to understand the neural basis of navigation shows excellent promise. We conclude with future directions for this research area.

Threshold-free cluster enhancement (TFCE) integrates cluster extent across thresholds to improve voxel-wise neuroimaging inference, but permutation testing makes it prohibitively slow for large datasets. Probabilistic TFCE (pTFCE) uses analytical Gaussian random field (GRF) p-values but discretises the threshold grid. Exact TFCE (eTFCE) eliminates discretisation via a union-find data structure but still requires permutations. We combine eTFCE's union-find for exact cluster-size retrieval with pTFCE's analytical GRF inference. The union-find builds the cluster hierarchy in one pass over sorted voxels and enables exact size queries at any threshold; GRF theory then converts these sizes to analytical p-values without permutations. Validation on synthetic phantoms (64^3, 80 subjects): FWER controlled at nominal level (0/200 null rejections, 95% CI [0.0%, 1.9%]); power matches baseline pTFCE (Dice >= 0.999); smoothness error below 1%; concordance r > 0.99. On UK Biobank (N=500) and IXI (N=563), significance maps form strict subsets of reference R pTFCE, which supports conservative error control. Implemented in pytfce (pip install pytfce): baseline completes whole-brain VBM in ~5s (75x faster than R pTFCE), hybrid in ~85s (4.6x faster) with exact cluster sizes; both >1000x faster than permutation TFCE.

Data-driven discovery of governing equations from time-series data provides a powerful framework for understanding complex biological systems. Library-based approaches that use sparse regression over candidate functions have shown considerable promise, but they face a critical challenge when candidate functions become strongly correlated: numerical ill-conditioning. Poor or restricted sampling, together with particular choices of candidate libraries, can produce strong multicollinearity and numerical instability. In such cases, measurement noise may lead to widely different recovered models, obscuring the true underlying dynamics and hindering accurate system identification. Although sparse regularization promotes parsimonious solutions and can partially mitigate conditioning issues, strong correlations may persist, regularization may bias the recovered models, and the regression problem may remain highly sensitive to small perturbations in the data. We present a systematic analysis of how ill-conditioning affects sparse identification of biological dynamics using benchmark models from systems biology. We show that combinations involving as few as two or three terms can already exhibit strong multicollinearity and extremely large condition numbers. We further show that orthogonal polynomial bases do not consistently resolve ill-conditioning and can perform worse than monomial libraries when the data distribution deviates from the weight function associated with the orthogonal basis. Finally, we demonstrate that when data are sampled from distributions aligned with the appropriate weight functions corresponding to the orthogonal basis, numerical conditioning improves, and orthogonal polynomial bases can yield improved model recovery accuracy across two baseline models.

State space models (SSMs) have recently achieved strong performance on long sequence modeling tasks while offering improved memory and computational efficiency compared to transformer based architectures. However, their evaluation has been largely limited to synthetic benchmarks and application domains such as language and audio, leaving their behavior on sparse and stochastic temporal processes in biological imaging unexplored. In this work, we introduce the Single Molecule Localization Microscopy Challenge (SMLM-C), a benchmark dataset consisting of ten SMLM simulations spanning dSTORM and DNA-PAINT modalities with varying hyperparameter designed to evaluate state space models on biologically realistic spatiotemporal point process data with known ground truth. Using a controlled subset of these simulations, we evaluate state space models and find that performance degrades substantially as temporal discontinuity increases, revealing fundamental challenges in modeling heavy-tailed blinking dynamics. These results highlight the need for sequence models better suited to sparse, irregular temporal processes encountered in real world scientific imaging data.

Efficient interaction with the visual world requires not only accurate object identification but also precise localization of objects in space. While spatial ("where") processing has traditionally been attributed to dorsal stream pathways, recent work has shown that object position can also be decoded from responses in ventral stream areas such as the inferior temporal (IT) cortex. However, because object position in these paradigms is tightly coupled to pixel-based location, it remains unclear whether ventral stream position signals reflect perceptually meaningful spatial representations or simply inherited retinotopic structure. To address this question, we used the motion aftereffect, a classic visual illusion that shifts perceived object position without changing retinal input. Combining large-scale intracortical recordings in macaque IT with matched human psychophysics, we found that motion adaptation induces systematic direction-opponent biases in IT population codes for object position that mirror human perceptual reports, despite identical pixel-level stimuli. These effects are accompanied by adaptation-driven changes in the geometry of IT population representations. We further tested whether artificial vision systems exhibit similar dynamics. Standard feedforward, recurrent, and state-of-the-art video-based neural networks accurately encode object position but fail to produce adaptation-induced position shifts. However, applying empirically derived transformations based on IT adaptation dynamics to model feature spaces is sufficient to generate similar biases. Together, these results indicate that IT represents object position in perceptually aligned coordinates and also highlight a gap between biological and artificial vision systems in capturing history-dependent spatial coding.

Agent-based models (ABMs) are widely used to estimate causal treatment effects via paired counterfactual simulation. A standard variance reduction technique is common random numbers (CRNs), which couples replicates across intervention scenarios by sharing the same random inputs. In practice, CRNs are implemented by reusing the same base seed, but this relies on a critical assumption: that the same draw index corresponds to the same modeled event across scenarios. Stateful pseudorandom number generators (PRNGs) violate this assumption whenever interventions alter the simulation's execution path, because any change in control flow shifts the draw index used for all downstream events. We argue that this execution-path-dependent draw indexing is not only a variance-reduction nuisance, but represents a fundamental mismatch between the scientific causal structure ABMs are intended to encode and the program-level causal structure induced by stateful PRNG implementations. Formalizing this through the lens of structural causal models (SCMs), we show that standard PRNG practices yield causally incoherent paired counterfactual comparisons even when the mechanistic specification is otherwise sound. We show that a remedy is to combine counter-based random number generators (e.g., Philox/Threefry) with event identifiers. This decouples random number generation from simulation execution order by making random draws explicit functions of the particular modeled event that called them, restoring the stable event-indexed exogenous structure assumed by SCMs.

This paper investigates the problem of covert communications in a heterogeneous wireless network where multiple communication modalities are used simultaneously. In this setup, a legitimate transmitter sends confidential data to its receiver by selecting multiple modalities with the goal of maximizing communication covertness against a passive adversary (Willie) while satisfying a transmission rate requirement. We analyze two distinct scenarios for a given observation time by Willie. The two scenarios are: (i) Willie knows the modalities selected by the friendly transmitter, and (ii) Willie is unaware of the selected modalities. We first derive the optimal detector for Willie that minimizes the detection error probability (DEP) in both cases. For the first scenario, we derive an exact expression for the DEP and provide a computationally efficient approximation. For the second scenario, we introduce the DEP expressions in the low-signal-to-noise ratio (SNR) regime at Willie. Building on this analysis, we propose a novel low-complexity modality set selection technique designed to maximize the DEP subject to a rate constraint. Numerical simulations validate the derived analytical expressions and demonstrate that the proposed modality set selection technique achieves near-optimal performance, outperforming benchmark schemes.

In Earth observation (EO) missions with Low Earth orbit (LEO) satellites, high-resolution image acquisition generates a massive data volume that poses a significant challenge for transmission under the limited satellite power budget, while LEO movement introduces dynamic systems. To enable efficient image transmission, this paper employs semantic communication (SemCom) with joint source-channel coding (JSCC), which focuses on transmitting meaningful information to reduce power consumption. Under a quality-of-service (QoS) requirement defined by image reconstruction quality, this work aims to minimize the total transmit power by jointly optimizing the JSCC encoder-decoder parameters and resource allocation. However, the implicit relationship among JSCC parameters, link quality, and image quality, coupled with the presence of mixed integer-continuous variables, makes the problem difficult to solve directly. To address this, a curve-fitting model is proposed to approximate the JSCC compression-SNR-quality relationship. Then, the joint compression ratio-resource allocation (JCRRA) algorithm is proposed to address the underlying problem. Numerical results demonstrate that the proposed method achieves substantial power savings compared to both greedy algorithms and conventional transmission paradigms.

With the deployment of large antenna arrays at high-frequency bands, future wireless communication systems are likely to operate in the radiative near-field (NF). Unlike far-field beam steering, NF beams can be focused on a spatial region with finite depth, enabling user multiplexing in both range and angle. In NF multiuser multiple-input multiple-output (MU-MIMO) systems, achievable rates are limited by interference arising from sidelobes in both the axial (range) and lateral (angle) dimensions. This work investigates how axial sidelobes (ASLs) vary with array geometry. Closed-form array gain expressions are derived to characterize ASLs for uniform planar arrays. Analytical results show that the uniform square array (USA) yields the lowest ASLs, followed by the uniform concentric circular array (UCCA), uniform linear array (ULA), and uniform circular array (UCA). Specifically, the USA achieves a peak sidelobe level (PSLL) of -17.6 dB versus -7.9 dB for the UCA. Numerical simulations confirm that the USA provides superior sidelobe suppression and highest sumrate performance.

Near-field propagation in extremely large-scale MIMO (XL-MIMO) enlarges the beam training (BT) search space by introducing an additional range dimension, which makes conventional codebook-based beam sweeping prohibitively expensive under limited pilot resources, especially for multiuser sub-connected hybrid architectures. This letter proposes a deep-learning-based interference-aware multiuser BT framework (DL-IABT) that directly predicts analog beam indices from a small number of uplink sensing measurements. By exploiting a subarray-level approximation, a far-field codebook is adopted to represent each subarray response with negligible mismatch. To enable end-to-end (E2E) learning, we derive a variant-MSE surrogate loss by eliminating the digital precoder through a closed-form MMSE solution from KKT conditions, which implicitly accounts for multiuser interference (MUI). The proposed network integrates a complex-valued sensing front-end, a shared complex-valued encoder, a Transformer-based multiuser predictor, and a scalable Gumbel--Softmax beam selection head. Simulation results show that DL-IABT achieves near-optimal sum-rate performance while providing markedly higher effective throughput under pilot overhead constraints.

Hybrid beamforming for extremely large-scale multiple-input multiple-output (XL-MIMO) systems is challenging in the near field because the channel depends jointly on angle and distance, and the multiuser interference (MUI) is strong. Existing deep learning methods typically follow either a decoupled design that optimizes analog beamforming without explicitly accounting for MUI, or an end-to-end (E2E) joint analog-digital optimization that can be unstable under nonconvex constant-modulus (CM), pronounced analog-digital coupling, and gradient pattern of sum-rate loss. To address both issues, we develop a complex-valued E2E framework based on a variant minimum mean square error (variant-MMSE) criterion, where the digital precoder is eliminated in closed form via Karush-Kuhn-Tucker (KKT) conditions so that analog learning is trained with a stable objective. The network employs a grouped complex-convolution sensing front-end for uplink (UL) measurements, a shared complex multi-layer perceptron (MLP) for per-user feature extraction, and a merged constant-modulus head to output the analog precoder. In the indirect mode, the network designs hybrid beamformers from estimated channel state information (CSI). In the direct mode where explicit CSI is unavailable, the network learns the sensing operator and the analog mapping from short pilots, after which additional pilots estimate the equivalent channel and enable a KKT closed-form digital precoder. Simulations show that the indirect mode approaches the performance of iterative variant-MMSE optimization with a complexity reduction proportional to the antenna number. In the direct mode, the proposed method improves spectral efficiency over sparse-recovery pipelines and recent deep learning baselines under the same pilot budget.

Following recent advances in flexible electronics and programmable metasurfaces, flexible intelligent metasurfaces (FIMs) have emerged as a promising enabling technology for next-generation wireless networks. A FIM is a morphable electromagnetic surface capable of dynamically adjusting its physical geometry to influence the radiation and propagation of electromagnetic waves. Unlike conventional rigid arrays, FIMs introduce an additional spatial degree of design freedom enabled by mechanical flexibility, which can enhance beamforming, spatial focusing, and adaptation to dynamic wireless environments. This added capability enables wireless systems to shape the propagation environment not only through electromagnetic tuning but also through controllable geometric reconfiguration. This article explores the potential of FIMs for next-generation wireless networks. We first introduce the main hardware architectures of FIMs and explain how they can be integrated into wireless communication systems. We then present representative application scenarios, highlighting the advantages of FIMs for future wireless networks and comparing them with other emerging flexible wireless technologies. To illustrate their potential impact, we present case studies comparing FIM-enabled architectures with conventional rigid-array systems, demonstrating the performance gains enabled by surface flexibility for both communication and sensing applications. Finally, we discuss key opportunities, practical challenges, and open research directions that must be addressed to fully realize the potential of FIM technology in future wireless communication systems.

Continuous aperture arrays (CAPAs) provide a theoretical upper bound on the performance of densely packed antenna arrays, but their analysis is limited by the lack of closed-form signal-to-noise ratio (SNR) distributions under realistic fading conditions. This paper derives accurate analytical expressions for the matched-filter SNR distribution of one-dimensional CAPAs in correlated Rayleigh environments under both the sinc and Jakes correlation models using the Karhunen-Loeve expansion. By applying a truncated hypoexponential model, we obtain accurate approximations for the probability density function and cumulative distribution function of the SNR that closely match simulations, including the outage probability region where precise characterization is critical. Compared to a standard gamma approximation, our approach provides significantly improved accuracy in this regime. Additionally, the CAPA system considered is shown to outperform discrete antenna arrays. The derived expressions enable tractable and accurate evaluation of CAPAs under practical channel models.

This letter investigates a novel uplink (UL) system that integrates power-domain non-orthogonal multiple access (PD-NOMA) with a continuous reconfigurable intelligent surface (CRIS). We analyze the effective CRIS-assisted channels under spatially correlated fading to accurately approximate the characteristic function of the cascaded channel. This allows the derivation of an expression for the bit error rate (BER), a key performance metric for UL PD-NOMA. We further utilize the derived BER expressions to introduce a joint optimization framework that minimizes the average BER via UL power allocation and dynamic RIS partitioning among the users. The analytical results are validated by simulations, and show that the proposed optimization scheme eliminates the BER floors that are associated with UL NOMA. The results also confirm the superiority of the optimized CRIS-NOMA scheme over conventional orthogonal multiple access (OMA) and non-optimized UL NOMA schemes.

Consider the signal-to-noise ratio (SNR) of a continuous fluid antenna system (CFAS) operating over a Rayleigh fading channel. In this paper, we extend traditional system assumptions and consider spatially coherent isotropic correlation, continuous positioning of the antenna rather than discrete, and the use of multi-dimensional space (1D, 2D and 3D). By focusing on the upper tail of the received SNR distribution (the high SNR probability (HSP)), we are able to derive asymptotically exact closed-form formulas for the HSP. Finally, these results lead to scaling laws which describe the increase in the HSP as we employ more dimensions and the optimal CFAS dimensions.

Quality control is an essential operation in manufacturing, ensuring products meet the necessary standards of quality, safety, and reliability. Traditional methods, such as visual inspections, measurements, and statistical techniques, help meet these standards but are often time-consuming, costly, and reactive. With the advent of AI/ML, manufacturers can shift from reactive to proactive approaches in quality control. This study applies ML-based models for predictive quality control in a real-world manufacturing setting. The case company produces castings for powertrain components in heavy vehicles, where poor control of core-making process parameters leads to costly defects. ML models were developed by analyzing data from two core-making machines, their processes, and maintenance logs to identify parameters associated with casting defects, enabling the prediction and prevention of potential defects before they occur. The results demonstrated good accuracy rates, helping quality and production teams identify and eliminate defective cores and thereby improving product quality and production efficiency.

Existing integrated sensing and communication (ISAC) platforms fail to fully utilize the shared spectrum and aperture resources for sensing, resulting in poor sensing performance. Specifically, weak target sensing on the ISAC platform, such as micro-deformation monitoring (mDM), suffers from inaccurate measurements due to poor sensing quality. In this paper, we propose an AI-assisted approach to alleviate the effect of poor sensing quality in the ISAC system by effectively removing the clutter. We begin by modeling the environment clutter model as a combination of the deterministic and stochastic signals to represent urban coverage scenarios around the base station (BS). A clutter suppression optimization problem is formulated to extract the micro-deformation displacement (mDD) from the original ISAC signals. We then propose a learnable template-matching (LTM) approach to mitigate the influences of clutters, thereby enhancing sensing quality. In particular, the electromagnetic (EM) signal feature of the mDD is embedded into the network to strengthen the mDM signal, and clutter filters are incorporated to suppress environmental clutter. Numerical results illustrate the superiority of our proposed approach concerning convergence speed and accuracy in mDD prediction. By deploying our approach to the BS measurement, the simulation-only trained LTM exhibits impressive performance in environment clutter separation and mDD estimation.

In this paper, we propose a novel free-space optical (FSO) communication system utilizing optical skyrmions. We introduce a scheme referred to as skyrmion number modulation (SkM), which employs index modulation by encoding information onto the skyrmion number, a topological invariant preserved during free-space propagation. This topological nature offers the potential for inherent robustness against atmospheric turbulence-induced wavefront distortions, which limit the performance of conventional FSO systems. More specifically, we demonstrate that the fluctuation of the received skyrmion number is mitigated by a proposed intensity-based masking technique. Finally, our performance analysis based on a discrete memoryless channel framework confirms that the proposed system exhibits near-ideal robustness under weak turbulence and supports high-order modulation in moderate regimes.

This paper presents a unified analytical and optimization framework for Standard Condition Number (SCN)-based detection in MIMO Integrated Sensing and Communication (ISAC) systems operating under noise uncertainty. Conventional detectors such as the Likelihood Ratio Test (LRT) and Energy Detector (ED) suffer from false-alarm inflation when interference or jamming alters the noise covariance. To overcome this limitation, the SCN detector, defined as the ratio of the largest to smallest eigenvalues of the sample covariance matrix is analytically characterized for the first time in an ISAC setting. Closed-form expressions for the false-alarm and detection probabilities are derived using random matrix theory for a two-antenna sensing receiver and generalized to arbitrary MIMO dimensions. The analysis proves that the SCN maintains a constant false alarm rate (CFAR) property and remains resilient to covariance mismatch, providing theoretical justification for its robustness in dynamic environments. Leveraging these results, a tractable ISAC power-allocation problem is formulated to minimize total detection error subject to communication rate and power constraints, yielding an interpretable sequential solution. Numerical evaluations verify the theory and demonstrate that the proposed SCN detector consistently outperforms LRT and eigenvalue-based benchmarks, particularly under strong interference and jamming typical of modern multiuser networks.

This paper investigates the susceptibility of Integrated Sensing and Communication (ISAC) systems to hostile jamming, focusing on an aerial Reconfigurable Holographic Surface (RHS)-aided unmanned aerial vehicle (UAV). The proposed framework, termed RHOSI, enhances ISAC's resilience by dynamically shaping the wireless propagation environment. Specifically, RHOSI introduces a strategy to improve jamming resistance by jointly optimizing transmit beamforming at the hybrid base station, RHS phase shift configuration, and UAV spatial deployment, while ensuring the required echo signal-to-interference-plus-noise ratios for reliable sensing. The resulting non-linear optimization problem features highly coupled variables, which are decomposed into sub-problems and solved using an alternating optimization (AO) approach. Simulation results confirm the practicality and effectiveness of RHOSI in significantly improving the throughput and robustness of ISAC under adversarial jamming.

Low Earth orbit (LEO) inter-satellite links (ISLs) must achieve joint synchronization and ranging under severe hardware impairments, namely oscillator phase noise, clock drift, and measurement outliers, exacerbated by rapid relative dynamics exceeding 7~km/s. In coherent Doppler processing, the frequency observable depends on the \emph{difference} between consecutive carrier phase states, creating a cross-epoch coupling structure that fundamentally affects estimation-theoretic performance limits. This paper makes three contributions. First, we prove analytically that this cross-epoch Doppler coupling is \emph{necessary} to avoid unbounded carrier phase uncertainty: without it, phase variance grows linearly without bound. Second, we derive a posterior Cramér-Rao bound (PCRB) via the Tichavský recursion that explicitly incorporates the resulting 10$\times$10 block information structure. Third, we propose a hybrid robust filtering framework combining hard gating for impulsive cycle-slip outliers with Huber M-estimation for heavy-tail contamination, using TASD-aware innovation covariance to account for cross-epoch uncertainty in residual normalization. Monte Carlo simulations at Ka-band confirm that the PCRB accurately lower-bounds estimator performance under nominal conditions, while the hybrid method reduces 95th-percentile phase error by 27--93\% compared to standard extended Kalman filtering across different outlier regimes.

We present an intelligent wearable system to monitor and predict mood states of elderly people during their daily life activities. Our system is composed of a wristband to record different physiological activities together with a mobile app for ecological momentary assessment (EMA). Machine learning is used to train a classifier to automatically predict different mood states based on the smart band only. Our approach shows promising results on mood accuracy and provides results comparable with the state of the art in the specific detection of happiness and activeness.

Audio-visual emotion recognition (AVER) methods typically fuse utterance-level features, and even frame-level attention models seldom address the frame-rate mismatch across modalities. In this paper, we propose a Transformer-based framework focusing on the temporal alignment of multimodal features. Our design employs a multimodal self-attention encoder that simultaneously captures intra- and inter-modal dependencies within a shared feature space. To address heterogeneous sampling rates, we incorporate Temporally-aligned Rotary Position Embeddings (TaRoPE), which implicitly synchronize audio and video tokens. Furthermore, we introduce a Cross-Temporal Matching (CTM) loss that enforces consistency among temporally proximate pairs, guiding the encoder toward better alignment. Experiments on CREMA-D and RAVDESS datasets demonstrate consistent improvements over recent baselines, suggesting that explicitly addressing frame-rate mismatch helps preserve temporal cues and enhances cross-modal fusion.

Understanding how isolated quantum many-body systems thermalize remains a central question in modern physics. We study the onset of ergodicity in a two-dimensional disordered Heisenberg Floquet model using digital quantum simulation on IBM's Nighthawk superconducting processor, reaching system sizes of up to $10\times10$ qubits. We probe ergodicity across different length scales by coarse-graining the system into spatial patches of varying sizes and introducing a measure based on the collision entropy of each patch, enabling a detailed study of when ergodic behavior emerges across scales. The high sampling rate of superconducting quantum processing units, together with an optimal sample estimator, allow us to access patches of sizes up to $3\times3$. We observe that as the Heisenberg coupling $J$ increases, the noiseless system undergoes a smooth crossover from subergodic to ergodic behavior, with smaller patches approaching their random-matrix-theory values first, thereby revealing a hierarchy across scales. In the region of parameter space where classical tensor-network simulations are reliable, small patches or small values of $J$, we find excellent agreement with the error-mitigated quantum simulation. Beyond this regime, volume-law entanglement and contraction complexity growth causes the cost of classical methods to rise sharply. Our results open new directions for the use of quantum computers in the study of quantum thermalization.

We develop a microscopic theory of thermalisation for a thermometer coupled to a many-body bath beyond standard Markovian and Fermi-golden-rule assumptions. By modeling interaction matrix elements in the non-interacting basis as independent random variables, we derive a diffusion-propagator expression for the reduced dynamics and show that relaxation is controlled by the distribution of interaction-induced level broadenings. The theory predicts a thermalisation timescale set by the inverse typical broadening and yields a non-Markovian generalization of global balance. Exact-diagonalization tests for heavy-tailed L{é}vy couplings, an all-to-all transverse-field Ising model, and the one-dimensional Imbrie model show good agreement with these predictions.

Loop corrections to primordial correlation functions are unavoidable due to the non-linear nature of gravity. Previous works have established a robust framework for computing the renormalised one-loop power spectra of scalar and tensor modes, but primarily in (near) de Sitter backgrounds. In this work, we develop a consistent renormalisation procedure applicable to inflationary backgrounds that strongly break de Sitter symmetries and generate scale-dependent features in the primordial spectra. Our analysis is performed within the Effective Field Theory (EFT) of inflationary fluctuations, allowing for arbitrary time dependence of the Wilson coefficients. We show that both ultraviolet divergences and tadpoles of the theory, despite their strong time and scale dependence, can be cancelled by a finite set of local counter-terms compatible with the EFT symmetries. Importantly, this result only relies on the existence of an initial phase of adiabatic evolution continuously related to the Bunch-Davies vacuum and holds independently of the precise time dependence of the background and of the free-field mode functions. We then study two concrete realisations, corresponding to resonant and sharp features. In both cases, all calculations are carried out exactly in the limit of small feature amplitude. We analyse perturbativity and provide the first explicit demonstration that the renormalised one-loop power spectrum generated by a localised feature along the inflationary trajectory vanishes both at large and small scales. Our scale-dependent renormalisation framework implies that models of primordial features used to fit CMB residuals are consistent with perturbativity bounds, and opens the door to systematic studies of loop corrections in more complicated scenarios relevant for scalar-induced gravitational waves and primordial black holes.

Efficient preparation of nonclassical bosonic states is a central requirement for quantum computing, simulation, and precision metrology. We study resource-efficient quantum state preparation in bosonic qudit systems using the SNAP-displacement (SD) protocol. Existing SD-based approaches typically require a large number of gates and SNAP phases, resulting in complex control pulses, increasing the ansatz duration, and amplifying the impact of photon-loss and control errors. In this work, we focus on the near- to medium-term regime, in which noisy quantum devices impose trade-offs on the fidelity that can be achieved, which must be taken into account. Specifically, we propose to optimize only a subset of the SNAP phases and introduce three progressively more general sparse ansatzes. To provide fine-grained control and identify the most suitable ansatz for a given target fidelity, we further employ a scalarized multi-objective optimization that trades off fidelity against either the number of phases or the duration of the ansatz. Numerical results for several target states and qudit dimensions up to $d=64$ show that these sparse ansatzes achieve favorable trade-offs compared to the fully parameterized SD protocol in both ideal and noisy settings, consistently reducing the number of required phases and suggesting a practical route to more efficient near- and medium-term bosonic state preparation.

It was proposed in \cite{Chu:2024qil} that a quantum black hole can be described by a fuzzy sphere together with a half-filled Fermi sea. In this paper we propose that the tunneling of a fuzzy sphere system to a small one describes the quantum decay of black hole by Hawking radiation. Since the Fermi sea shrinks and the quantum mechanical Hamiltonian conserves fermion number, the amplitude of transition naively vanishes unless the tunneling path provides exact number of zero modes to soak up the excess fermi states. We show that a monopole on fuzzy sphere does exactly that. This fixes the tunneling path. The resulting tunneling rate reproduces Page's result for the semi-classical decay rate of black hole. We propose to identify the fermi states released by the monopole as the Hawking radiation. At the level of probability, the Hawking radiation is found to be given by a Boltzmann distribution at the Hawking temperature. One can go beyond the probabilistic description by determining the full wave function of the multi-partite Hawking quanta. This is possible with a real time formulation of the tunneling process. Unitarity is manifest in our quantum mechanics.

We introduce an algebraic structure for studying state-independent contextuality arguments, a key form of quantum non-classicality exemplified by the well-known Peres-Mermin magic square, and used as a source of quantum advantage. We introduce \emph{commutation groups} presented by generators and relations, and analyse them in terms of a string rewriting system. There is also a linear algebraic construction, a directed version of the Heisenberg group. We introduce \emph{contextual words} as a general form of contextuality witness. We characterise when contextual words can arise in commutation groups, and explicitly construct non-contextual value assignments in other cases. We give unitary representations of commutation groups as subgroups of generalized Pauli $n$-groups.

We present a shifted-geodesic framework for computing gravitational-wave fluxes from spinning test bodies moving on bound orbits of Kerr black holes. The method provides a simple and efficient means of evaluating energy and angular momentum fluxes incorporating the leading effect of the smaller body's spin. Because post-adiabatic corrections, including secondary spin contributions, are subdominant to the leading adiabatic terms, this approximation is well justified. In particular, we find that oscillatory spin terms typically contribute very little to fluxes, but their contribution to the description of orbits is computationally expensive, making such terms a natural target for approximation. In our framework, orbital frequencies and integrals of the motion are perturbed to include spin effects, while the trajectory retains the global structure of geodesic motion. This simplifies the computation of gravitational radiation. The shifted-geodesic approximation is most reliable for orbits with lower eccentricity, lower inclination, and larger semi-latus recta. The approximation becomes less reliable as we approach the separatrix between stable and unstable orbits; fortunately, many inspirals spend less time in this region of parameter space. A diagnostic inspiral evolution shows very small dephasing due to use of the shifted-geodesic approximation ($\approx10^{-2}$ radians over 1 year), confirming that spin-induced flux corrections can be accurately included using this simple modification to a geodesic trajectory. This approximation provides a rapid and convenient way to compute spinning-body orbits, but is not intended to replace more accurate treatments. We propose it as a pragmatic alternative when speed and simplicity are prioritized, enabling efficient EMRI/IMRI flux calculations and supporting parameter-space studies for LISA. (Abridged)

Are Gaussian measurements enough to distinguish between Gaussian states? Here, we tackle this question by focusing on the max-relative entropy as an operational distinguishability metric. Given two general multimode Gaussian states, we derive a condition, based on their covariance matrices, that completely determines whether or not there exists an optimal Gaussian measurement achieving the max-relative entropy. When the condition is satisfied, we find this optimal measurement explicitly. When the condition is not met, there is a strict gap between the distinguishability achievable by Gaussian measurements and the unconstrained max-relative entropy in which all measurements are allowed. We illustrate our results in the single-mode setting, and show examples of states for which this gap can be made arbitrarily large, revealing novel instances of Gaussian data hiding.

Driven nonlinear oscillators constitute a universal paradigm for understanding synchronization, frequency pulling, and frequency comb formation in nonequilibrium systems. Here, we realize such an emergent nonlinear oscillator in strongly interacting cesium Rydberg vapor, where coherent optical excitation, dissipation, and long-range interactions give rise to a driven-dissipative time crystal phase with intrinsic oscillation frequencies. Applying a radio-frequency (RF) field allows controlled tuning of the intrinsic oscillation frequency. Under RF heterodyne conditions, we observe intermodulation, frequency pulling, and, at strong drive, the emergence of a comb-like spectrum in the atomic coherence. We quantitatively capture these observations using a four-level mean-field model and demonstrate a classical analogue with a driven Van der Pol oscillator. Our results establish interacting Rydberg ensembles as a tunable platform for exploring nonequilibrium time-crystalline order, nonlinear synchronization, and frequency comb generation in many-body atomic systems.

The sky localisation of about $10$ to $100~\text{deg}^2$, which is expected to be achieved in all-sky blind searches for gravitational waves from supermassive black hole binaries (SMBHBs) with Pulsar Timing Array (PTA) experiments, has long been posed as a prohibitive factor in utilising these sources as standard sirens for precision cosmology. We propose a solution to this problem, which makes use of targeted searches rather than all-sky blind searches for SMBHBs. Using our simulated data informed by current PTA observations, we show that the Chinese Pulsar Timing Array (CPTA) alone could infer the Hubble constant with a precision of 2~km/s/Mpc. Such precision in an independent cosmological probe could provide decisive support in the resolution of the Hubble tension. We demonstrate the application of our method to several simultaneously observed SMBHBs, as well as the method's robustness against confusion between the host galaxies of SMBHB sources in realistic observing scenarios.

High-fidelity quantum operations are the cornerstone of fault-tolerant quantum computation. In open quantum systems, traditional optimal control only passively resists decoherence, leaving environment-induced uncertainty as a fundamental performance bottleneck. To overcome this, we propose a new optimal control framework with flag ancillas and the Flag-GRAPE algorithm, which can actively tailor the system's noise structure. Through embedding post-selection directly into the objective function, Flag-GRAPE correlates decoherence errors with the ancilla's unexpected state. Subsequent measurement and post-selection effectively expel this uncertainty, circumventing the fidelity bounds of traditional control. Numerical simulations in a superconducting quantum circuit demonstrate a $51\%$ reduction in infidelity compared to traditional closed-system pulses and also show that such enhancement is robust across broad noise regimes. Furthermore, by actively converting unstructured decoherence into heralded erasure errors, Flag-GRAPE is inherently compatible with quantum error correction. We demonstrate this by initializing a logical cat-code state, showing that the combination between Flag-GRAPE and QEC yields immediate state preparation enhancements. This new framework can reduce hardware overhead for fault-tolerant architectures and open up a practical path toward logical state preparation gain in near-term experiments.

Developing quantum algorithms to simulate fluid dynamics has become an active area of research, as accelerating fluid simulations could have significant impact in both industry and fundamental science. While many approaches have been proposed for simulating fluid dynamics on quantum computers, it is largely unclear whether these algorithms will provide speedup over existing classical approaches. In this paper we give evidence that quantum computers cannot significantly outperform classical simulations of fluid dynamics in general. We study two models of fluids: the Korteweg-de Vries (KdV) equation, which models shallow water waves, and the incompressible Euler equations, which model ideal, inviscid fluids. We show that any quantum algorithm simulating the KdV equation or the Euler equations for time $T$ requires $Ω(T^2)$ and $e^{Ω(T)}$ copies of the initial state in the worst case, respectively. These lower bounds hold for the task of preparing the final state, and similar bounds hold for history state preparation. We prove the lower bound for the KdV equation by investigating divergence of solitons. For the Euler equations, we show that instabilities enable fast state discrimination.

The coordination complexes of Yb(III) exhibit some of the longest spin coherence times among 4f compounds, making them a promising platform for molecular quantum technologies. While spin-phonon relaxation remains a limiting factor for coherence times even at low temperature, its control through chemical design has the potential to push these spin qubits prototypes beyond current limits. With the aim of providing insights on how to chemically control spin-phonon relaxation, we here present a full ab initio study of spin-phonon dynamics for three Yb(III) molecules exhibiting minimal chemical differences, yet quantitatively different spin relaxation times. Results show that low-temperature relaxation is governed by Raman processes triggered by a small group of largely delocalized low-energy phonons. The analysis of these contributions highlights that the modulation of spin-phonon coupling by molecular structure modifications beyond the first coordination shell are highly non-trivial in nature and hard to rationalize in simple chemical terms. These findings call for a conceptual step change from the attempt to use simple magneto-structural correlations to interpret the effect of molecular structural modifications on spin-phonon relaxation, and present predictive first-principles frameworks as a potential driving force of future chemical design strategies

In quantum theory, equilibrium statistical mechanics is usually formulated through the canonical ensemble, whose privileged status is tied to the Euclidean continuation of time evolution. The microcanonical ensemble, by contrast, is commonly introduced as a separate spectral construction. In this work we show that this asymmetry is representational rather than structural. We formulate the system in an extended Hilbert space in which time is promoted to an auxiliary clock degree of freedom and physical states are selected by a reparametrization-invariant constraint operator $\hat C = \hat P_T + \hat H$. The corresponding projector $δ(\hat C)$ provides a single unified object from which both canonical and microcanonical ensembles emerge as complementary projections in the clock sector. In the clock-time representation, a purely imaginary clock separation yields the Euclidean kernel and the canonical partition function. In the conjugate clock-energy representation, the same projector reduces to the spectral operator $δ(\hat H-E)$ and hence to the microcanonical density of states. The main consequence is structural: canonical and microcanonical statistics need not be introduced as independent constructions, since both are already encoded in the same constrained quantum dynamics.

We present a method for analyzing general time series by employing the history state formalism of quantum mechanics. This formalism allows us to describe a complete evolution based on a single quantum state, the history state, which simultaneously includes -also as a quantum system- the reference clock. It naturally leads to the concept of system-time entanglement, with the ensuing entanglement entropy constituting a measure of the effective number of distinguishable states visited in the history. Through a quantum coherent state embedding of the time series data, it is then possible to associate a quantum history state to the series. The gaussian overlap between these coherent states provides thus a smooth measure of distinguishability between the series data. The eigenvalues of the corresponding overlap matrix determine in fact the entanglement spectrum and entropy of the history state, which provide a rigorous characterization of the evolution. As illustration, the formalism is applied to typical financial time-series data. Through the entanglement entropy and spectrum, different evolution regimes can be identified. Entanglement based volatility indicators are also derived, and compared with standard volatility measures.

Quantum channel discrimination is a fundamental task in quantum information processing. In the one-shot regime, discrimination between two candidate channels is characterized by the diamond norm. Beyond this basic setting, however, many scenarios in distributed quantum information processing remain unresolved, motivating notions of distinguishability that capture the power of the available resources. In this work, we formulate a theory of testers for bipartite channel discrimination, leading to the concept of the entanglement cost of bipartite channel discrimination: the minimum Schmidt rank $k$ of a shared maximally entangled state required for local protocols to achieve the globally optimal success probability. We introduce $k$-injectable testers as a tester-based description of entanglement-assisted local discrimination and, in particular, study the class of $k$-injectable positive-partial-transpose (PPT) testers, which constitutes a numerically tractable relaxation of the practically relevant class of LOCC testers. For every $k$, we derive a semidefinite program (SDP) for the optimal success probability, which in turn yields an efficiently computable one-shot PPT entanglement cost. To render these optimization problems numerically feasible, we prove a symmetry-reduction principle for covariant channel pairs, thereby reducing the effective dimension of the associated SDPs. Finally, by dualizing the SDP, we derive bounds on the composite channel-discrimination problem and illustrate our framework with proof-of-principle examples based on the depolarizing channel, the depolarized SWAP channel, and the Werner--Holevo channels.

The Bernstein-Vazirani (BV) algorithm is frequently taught as a canonical example of quantum parallelism, yet the standard interference-based explanation often obscures its underlying simplicity. We present a geometric reframing in which the Hadamard gate "wrapping" acts as a global basis rotation rather than a generator of computational complexity. This perspective reveals that the algorithm is effectively a classical linear computation over GF(2) performed in the conjugate Fourier basis, with the apparent parallelism arising from coordinate transformation. Building on Mermin's earlier pedagogical shortcut, which presented a 'classical' circuit equivalent but stopped short of explicitly labeling it as such, we elevate this to a formal geometric framework. In the extension, we distinguish between globally rotated circuits--which we reveal as classical linear computations--and topologically twisted circuits that generate quantum entanglement. We introduce a pedagogical taxonomy distinguishing (1) pure computational-basis circuits, (2) globally rotated circuits (exemplified by Bernstein-Vazirani), and (3) topologically twisted circuits involving non-aligned subsystem bases. This framework allows viewing the Gottesman-Knill theorem from a new angle, extends students' understanding of phase kickback and the 'Ricochet Property'. Furthermore, it provides a more intuitive starting point for explaining Bell-pair extensions through concrete circuit derivations and Qiskit simulations suitable for undergraduate quantum information courses. The outlook explores how this geometric view paves the way for understanding entanglement as topological twists.

Fast and high-fidelity qubit measurement plays a key role in quantum error correction. In superconducting qubits, measurement is typically performed using a resonant microwave drive on a readout resonator dispersively coupled to the qubit. Shorter measurement times require larger numbers of photons populating the readout resonator, which ultimately leads to undesired measurementinduced state transitions (MIST) of the qubit. MIST can be particularly problematic because these transitions often leave the qubit in a high energy state, and the MIST locations in readout parameter space drift as a function of qubit offset charge. In transmon qubits, these drifts have been avoided using very large qubit-resonator detunings or dedicated offset charge biases. In this work, we take an alternative approach and add an inductive shunt to the transmon to eliminate the offset charge dependence and stabilize the MIST. We experimentally characterize MIST in several different inductively-shunted transmons, in agreement with quantum and semiclassical models for MIST. These results extend to other inductively-shunted qubits.

Four-dimensional gravity admits many equivalent formulations - metric, Einstein-Cartan, teleparallel, McDowell-Mansouri, among others - each offering distinct advantages, particularly, in view of quantization. We propose a new formulation based on Weyl spinor-valued 1-forms, ultimately encoding the frame-field data. Starting from a topological field theory with a global $\mathrm{SL}(2,\mathbb{C})$ symmetry, we show that promoting this symmetry to a local gauge symmetry leads to the emergence of gravity. We analyze the covariant phase space of this theory, its symmetries and charge structure and explore the role of admissible corner terms together with their impact on boundary charges and their algebra. We study several extensions of this framework, including the incorporation of a cosmological constant and a novel $ G \rightarrow 0 $ scaling limit obtained from this model. The presence of the frame field already at the topological level allows point particles to be coupled uniformly in both the topological and gravitational theories. We perform a detailed Hamiltonian analysis of the theory and clarify the implementation of the reality conditions. We argue that this formulation provides structural features that make it particularly well suited for both discretization and quantization.

The Clifford Hierarchy has been a central topic in quantum computation due to its strong connections with fault-tolerant quantum computation, magic state distillation, and more. Nevertheless, only sections of the hierarchy are fully understood, such as diagonal gates and third level gates. The diagonal part of the hierarchy can be climbed by taking square roots and adding controls. Similarly, square roots of Pauli gates (first level) are Clifford gates (climb to the second level). Based on this theme, we study gates whose square roots climb to the next level. In particular, we fully characterize Clifford gates whose square roots climb to the third level.

This is a book about operational probabilistic theories. The standard approach in such theories is from a time forward perspective. In this book we mostly take a time symmetric perspective. This presents a branding problem. Is this a niche book merely about time symmetry? No. This is a comprehensive book about operational probabilistic theories, but mostly from a time symmetric perspective. In fact, this book consists of (1) a simple book about simple operations having simple causal structure (where all the inputs are before all the outputs), and (2) a complex book about complex operations that can have complicated causal structure (a complex operation is equipped with a causal diagram). For the simple case we are able to show that the time symmetric perspective is equivalent to the time forward perspective. In each book we set up (A) operational probabilistic theories (OPTs) in terms of operations, (B) Operational Quantum Theory (OQT) in terms of operator tensors which correspond to operations, and (C) the theory of Hilbert objects which can be doubled up to give operator tensors. Operations are required to be physical which guarantees that circuits built out of operations have probabilities between 0 and 1 and that certain causality conditions are met. We prove that when we wire together operations the resulting networks are also physical. We model Sorkin's impossible measurements with complex operations and show that physicality prevents anomalous signalling. We develop diagrammatic notation for Hilbert objects. This includes mirrors for doubling up and mirror theorems. We use this framework to prove time symmetric causal dilation theorems for various causal diagrams.

We study in detail the dynamics of unstable two-level quantum systems by adopting the Bloch-vector representation. We identify a novel class of critical scenarios in which the so-called energy-level and decay-width vectors, ${\bf E}$ and ${\bfΓ}$, are orthogonal to one another, and the parameter $r = |{\bf Γ}|/(2|{\bf E}|)$ is less than~1. Most remarkably, we find that critical unstable qubit systems exhibit atypical behaviours like coherence--decoherence oscillations when analysed in an appropriately defined co-decaying frame of the system. By making use of a Fourier series decomposition, we define anharmonicity observables that quantify the degree of non-sinusoidal oscillation of a CUQ. We apply the results of our formalism to the neutral-meson systems and derive generic upper limits on these new observables. In particular, we provide a compilation table of all well-explored meson--antimeson two-level systems in terms of Bloch-sphere parameters.

Randomized compiling (RC) is an established tool to tailor arbitrary quantum noise channels into Pauli errors. The effect of both spatial and temporal noise correlations in randomly compiled circuits, however, is not fully understood. Here, we show that for a broad class of correlated Gaussian noise, RC reduces both the strength and temporal range of correlations. For Clifford circuits, we derive a simple analytical expression for the circuit fidelity of randomly compiled circuits. Surprisingly, we show that this fidelity is always increased by the presence of correlations, suggesting that correlations are a resource in randomly compiled circuits. To leading order in system-bath coupling, we also show that RC suppresses the quantum component of bath correlations, implying that one can safely treat weak noise as being classical. Finally, through extensive numerical simulations, we show that our results remain valid for many relevant non-Clifford circuits. These results clarify how RC mitigates memory effects and enhances circuit robustness.

We present a measurement-driven quantum framework for evaluating overlap and Hamiltonian matrix elements in spin-coupled generalized valence bond (SCGVB) wavefunctions. The approach targets a central difficulty of nonorthogonal valence-bond methods: estimating matrix elements between distinct, generally nonorthogonal configuration state functions. Rather than preparing the full wavefunction on quantum hardware, we reformulate the required quantities as vacuum expectation values of Pauli-string operators that can be accessed using shallow, ancilla-free circuits composed of local Clifford rotations and computational-basis measurements. In contrast to Hadamard-test-based matrix-element estimation, this construction avoids ancilla qubits and controlled operations by reducing the problem to local Pauli measurements. This separates the algebraic construction of the SCGVB problem from the measurement task executed on the quantum register and yields a low-depth strategy compatible with near-term architectures. We demonstrate the framework on square and rectangular H4 using quantum-circuit emulation, where the resulting overlap and Hamiltonian matrices reproduce classical Lowdin-based references with good accuracy across the geometries considered, and where derived Coulson-Chirgwin weights remain chemically consistent. These results support the feasibility of measurement-based quantum assistance for nonorthogonal SCGVB expansions and provide a practical route for incorporating quantum measurements into valence-bond electronic-structure workflows.

One-axis twisting (OAT) interaction is a pivotal resource for manipulating quantum states of atomic ensembles, enabling spin squeezing, atomic-cat-state generation, and weak-phase amplification. Current implementations of OAT dynamics predominantly rely on the Tavis-Cummings model of light-atoms coupling; however, this approach inevitably introduces an additional Stark term that entangles the light with the atoms, which compromises the unitarity of OAT evolution and thereby degrades the OAT-based control precision. Here we propose a scheme based on a driven Tavis-Cummings model to achieve near-unitary OAT evolution. We demonstrate that both constant and time-varying driving of an atoms-cavity hybrid system can realize near-unitary OAT evolution, albeit with distinct coupling strength. Furthermore, when atomic dissipation is taken into account, we find that the time-varying-driving scheme exhibits superior resistance to decoherence. Our approach is broadly applicable to a variety of atomic platforms, including cold atoms, trapped ions, and nitrogen-vacancy centers.

Gravitational waves (GW) emitted by binary systems allow us to perform precision tests of general relativity in the strong field regime. Ringdown signals allow for probing black hole mass and spin with high precision in GW astronomy. With improvements in current and next-generation GW detectors, developing likelihood-free parameter inference methods is crucial. This is especially important when facing challenges such as non-standard noise, partial data, or incomplete signal models that prevent the use of analytical likelihood functions. In this work, we propose an amortized simulation-based inference strategy to estimate ringdown parameters directly. Specifically, our method is based on amortized neural posterior estimation, which trains a neural density estimator of the posterior for all data segments within the prior range. The results show that our trained amortized network achieves statistically consistent parameter estimates with valid confidence coverage compared to established Markov-chain methods, while offering inference speeds that are orders of magnitude faster. Furthermore, we evaluate the robustness of the method against transient noise contamination. Our analysis reveals that the timing of glitch injection has a decisive impact on estimation bias, particularly during the tail of a signal with sparse information. Glitch strength is positively correlated with estimation error, but has limited effect at low signal-to-noise ratios. Mass and spin parameters are most sensitive to noise. This study not only provides an efficient and accurate inference framework for ringdown analysis but also lays a foundation for developing robust data-processing pipelines for future GW astronomy in realistic noise environments.

Certifying that quantum randomness generated by untrusted devices is unpredictable to an attacker (say, Eve) is crucial for device-independent security. Bipartite protocols where only one of the parties is trusted are termed one-sided device-independent (1SDI) or steering-based protocols, where the untrusted party (say, Alice) performs measurements on her part of a bipartite entangled state to steer the subsystem of the trusted party (say, Bob) into different ensembles (collectively, an assemblage) of quantum states. Recent work has shown that an assemblage has certified randomness if and only if it is realizable by a set of measurements that are star-incompatible, i.e., the measurement setting of interest for the guessing probability of Eve is incompatible with at least one of the remaining measurement settings of Alice. However, it remains conceivable that there exist star-incompatible measurements that cannot certify steering-based randomness, just like there exist incompatible measurements that cannot certify bipartite Bell nonlocality. Here we prove that any set of star-incompatible measurements can generate steering-based randomness, thereby establishing an equivalence between the two notions. We further introduce a weight-based measure of star-incompatibility and lower bound the amount required to certify a given randomness, capturing the qualitative and quantitative interplay between the quantum resources of star-incompatibility and steering-based randomness.

We analyze propagation and detection of two-photon states expanded in Zernike modes through atmospheric turbulence using the extended Huygens-Fresnel formalism. For SPDC states prepared with a single Zernike pump mode, we analytically reduce the 8-dimensional continuous propagation integrals to an exact, discrete modal expansion. In the absence of turbulence, Zernike addition enforces conservation of azimuthal index and a strict radial-order bound. Turbulence relaxes these constraints, driving structured azimuthal and radial crosstalk dominated by low-order aberration modes. By explicitly removing the lowest-order terms from the discrete turbulence sum, we demonstrate that partial adaptive optics correcting only up to the sixth radial order is sufficient to heavily suppress this crosstalk and restore near-ideal spatial correlations.

This paper investigates quantum simulation algorithms for the Liouville equation in geometrical optics with partial transmission and reflection at sharp interfaces, based on the Schrödingerization method. By means of a warped phase transformation in one higher dimension, the Schrödingerization method converts linear partial differential equations into a system of Schrödinger-type equations with unitary evolution, thereby rendering them suitable for quantum simulation. In this work, the Schrödingerization method is combined with a Hamiltonian-preserving scheme that incorporates partial transmission and reflection into the numerical flux. A main difficulty is that the interface treatment in the classical scheme relies on threshold-dependent "if/else" procedures, making it highly nontrivial to reformulate the method in a matrix form suitable for quantum simulation. To overcome this difficulty, we encode the interface conditions into a partial transmission and reflection matrix prepared a priori, rather than during the time evolution. We present detailed constructions of the resulting quantum algorithms and show through complexity analysis that the proposed methods achieve polynomial quantum advantage in the precision parameter $ε$ over their classical counterparts.

We present an algebraic framework for approximate model reduction of Markovian open quantum dynamics that guarantees complete positivity and trace preservation by construction. First, we show that projecting a Lindblad generator on its center manifold -- the space spanned by eigenoperators with purely imaginary eigenvalue -- yields an asymptotically exact reduced quantum dynamical semigroup whose dynamics is unitary, with exponentially decaying transient error controlled by the generator's spectral gap. Second, for analytic perturbations of a Lindblad generator with a tractable center manifold, we propose a perturbative reduction that keeps the reduced space fixed at the unperturbed center manifold. The resulting generator is shown to remain a valid Lindbladian for arbitrary perturbation strengths, and explicit finite-time error bounds, that quantify leakage from the unperturbed center sector, are provided. We further clarify the connection to adiabatic elimination methods, by both showing how the algebraic reduction can be directly related to a first-order adiabatic-elimination and by providing sufficient conditions under which the latter method can be applied while preserving complete positivity. We showcase the usefulness of our techniques in dissipative many-body quantum systems exhibiting non-stationary long-time dynamics.

High-energy colliders, exemplified by the CERN's Large Hadron Collider (LHC), constitute genuine quantum machines. In alignment with Richard Feynman's foundational vision for quantum computing, collider physics emerge therefore as a prime candidate for quantum simulations. Prospective applications include Quantum Machine Learning for collider data analysis, accelerated evaluation of complex multiloop Feynman diagrams, efficient jet clustering, enhanced parton shower simulations, and related computational challenges. We discuss two specific applications: the identification of causal structures in multiloop vacuum amplitudes, a fundamental component of the Loop-Tree Duality exhibiting deep connections to graph theory; and high-dimensional function integration and sampling. The latter constitutes an initial step toward realizing a fully fleged quantum event generator capable of operating at high perturbative orders.

Quantum teleportation uses a shared entangled resource, local operations, and a digitally error-corrected classical channel to transfer quantum states between distant parties. We introduce a hybrid teleportation-direct transmission protocol for state transfer that still exploits entanglement, but replaces classical communication and digital error correction with an analog feedforward through a noisy quantum channel. We show that quantum teleportation outperforms this protocol if the communication channel reduces the entanglement of all bipartite states having the same amount of entanglement as the resource; otherwise, the hybrid protocol is optimal. We apply our result to the state transfer of a uniformly distributed coherent-states codebook, highlighting experimentally relevant scenarios where our protocol is most effective. Our findings are directly relevant to both optical and superconducting microwave channels, where analog feedforward techniques have been recently implemented.

Coherent states offer a promising path for near-term quantum computing due to their inherent protection against bit-flip noise. However, their large photon numbers can be challenging for numerical simulation. This paper introduces an effective model, representing coherent-state quantum annealing using spin-1/2 degrees of freedom. We demonstrate that this model yields accurate predictions for realistic experimental settings and can therefore serve as a practical tool for optimizing future quantum hardware.

These notes are a short introduction to the mathematical theory of open quantum systems. They are meant to serve as an entry point into a broad research area which has applications across the quantum sciences dealing with systems subjected to external noise. The guiding idea is to let the key structures of the theory emerge from a concrete model. By working through the dissipative Jaynes-Cummings model the reader will dis- cover explicitly how irreversible dynamics arises from a unitary system-reservoir evolution. The notions of the continuous mode limit, correlation functions, spectral density appear in a natural manner and lead to the evolution equation of the open system in form of a master equation. This sets the stage for the more general analysis of completely positive, trace preserving (CPTP) maps and the study of quantum dynamical semigroups. We motivate and prove the Kraus representation theorem, the dilation theorem and the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) theorem. Working through the exercises (for which full solutions are supplied) will reinforce the ideas introduced in the main text.

In this work, we introduce a screened Simpson-Visser regular solution and perform a comprehensive study of its physical and observational properties. We begin by analyzing the thermodynamic behavior of the black hole, including detailed investigations of the Hawking temperature, Gibbs free energy, and specific heat, which provide insights into its stability and phase structure. Next, we examine the geodesic structure of the spacetime, considering both massless (photon) and massive (timelike) particles. In particular, we study the photon sphere, the corresponding black hole shadow, and the innermost stable circular orbits (ISCO), which are crucial for understanding the motion of matter and light around the black hole. Furthermore, we explore the black hole's energy-emission rate radiation, highlighting the effects of the modified geometry on observational signatures. Finally, we investigate the topological aspects of the black hole, analyzing both the thermodynamic topology and the photon sphere's topological properties. Our analysis demonstrates the intricate interplay between the spacetime geometry, geodesic motion, and black hole thermodynamics, offering a deeper understanding of this class of regular black holes and their potential observational consequences.

Synchronization is a hallmark of collective behavior in classical nonlinear systems, yet its realization as a robust many-body phenomenon in coherent quantum systems remains largely unexplored. Here we demonstrate symmetry-protected quantum synchronization and a quantum chimera state in coherent Floquet dynamics on programmable superconducting quantum processors. By implementing stroboscopic evolution of a two-dimensional Heisenberg model on IBM heavy-hex devices, we observe that initially phase-randomized spins spontaneously self-organize into coherent lattice-wide oscillations. On 28 qubits, synchronization persists even for strongly randomized initial states and is stabilized by SU(2) symmetry, as confirmed by explicit symmetry breaking. Scaling up to 156 qubits reveals a qualitatively distinct regime. For weak initial randomness, global synchronization extends across the device. For strong randomness, the system fails to synchronize globally, yet subsets of qubits exhibit robust local phase coherence under homogeneous unitary dynamics. This coexistence of globally desynchronized and locally synchronized regions constitutes a quantum analogue of a classical chimera state. Statevector and matrix-product-state simulations reproduce both the symmetry-protected synchronization and the chimera coexistence, demonstrating that these phenomena arise from the intrinsic Floquet many-body dynamics. Our results establish symmetry-protected synchronization and quantum chimera states as experimentally accessible nonequilibrium dynamical phases in programable many-body quantum systems.

We propose that quantum mechanics follows from a single hypothesis: equality has finite resolution. Replacing the binary predicate $x = y$ with a graded distinguishability kernel $K(x,y) \in [0,1]$ forces three structural consequences: finite capacity ($N$ perfectly distinguishable states), relational completeness (all structure reduces to $K$-relations, and no measurement orientation is privileged), and reversible dynamics. We formalize the first two as axioms; a structural Leibniz condition within the saturation axiom forces permutation invariance of basis elements, and the full dynamical structure (cyclic evolution of order $N$, complex coefficients, and continuous unitary time evolution) is then uniquely determined. From these axioms (with regularity conditions derived in Appendix B: complex coefficients $\mathbb{C}$ are the unique field supporting cyclic dynamics and relational isotropy; deterministic hidden variables require $Ω(N^2)$ bits of storage (for prime-power $N$; exceeding $\log_2 N$ for all $N \geq 3$); the Born rule $p_k = |c_k|^2$ is the unique probability assignment preserving statistical distinguishability under reversible dynamics; and local tomography follows from $\mathbb{F} = \mathbb{C}$ with tensor product composition. Standard quantum mechanics is the $N \to \infty$ limit; finite $N$ provides a natural UV cutoff. The single free parameter is capacity $N$.

As cosmology rapidly approaches the data-dominated phase of stage IV large scale structure surveys, the modelling of nonlinear scales has become a serious challenge that faces the community, particularly when analysing models beyond $w$CDM. In this work, we emulate the matter power spectrum in a phenomenological parameterisation of modified gravity in which a time-varying effective gravitational constant $μ$ and a gravitational slip $η$ are binned in redshift. We are able to achieve accuracy $<1\%$ in the modified gravity boost relative to COLA (COmoving Lagrangian Acceleration) simulations. We forecast the constraining power for each bin using a simulated $3\times 2$pt LSST Y10-like data vector and a $6\times 2$pt LSST Y10 x Simons Observatory cosmic microwave background (CMB) lensing data vector. We recover the characteristic degeneracy between $μ$ and $η$ previously identified in Fisher forecasts and demonstrate that the best-constrained direction corresponds to the combination $Σ=μ(1+η)/2$ which governs the lensing potential. We show that while large scale structure is sensitive to growth of structure at low redshift, CMB lensing extends the sensitivity to a higher redshift range. These results demonstrate that fast emulation of nonlinear modified-gravity effects enables full Bayesian analyses of model-agnostic gravity parameterisations with realistic survey data vectors and astrophysical systematics.

Achieving fault tolerance with superconducting quantum processors requires qubits to operate within the regime of threshold theorems based on the Born-Markov approximation. This approximation, which models dissipation as constant energy decay into a memoryless environment, breaks down when qubits couple to long-lived two-level systems (TLSs) that become polarized during operation and retain memory of past qubit states. Here, we show that non-Poissonian quantum jump traces carry the information required to distinguish long-lived TLSs from the standard Markovian bath. By fitting the Solomon equations to measured quantum jumps dynamics arising naturally due to thermal fluctuations, we can disentangle the coupling of the qubit to the two environments. Sweeping the qubit frequency reveals distinct peaks, each associated with a TLS that outlives the qubit, providing a handle to understand their microscopic origin.

In quantum mechanics, not everything that can be observed can be observed simultaneously. Observational data exhibits \emph{contextuality} -- a generalisation of nonlocality -- if the result of an observation is necessarily dependent on which combination of observables was measured. This article gives a mathematical introduction to contextuality, emphasising its nature as a general feature of probability theory and logic, rather than of any particular quantum theory.

The standard gravitational-baryogenesis operator $λ\,\nabla_μR\,J^μ$, with $λ\equiv ε/M_\ast^{2}$, is usually treated as a spectator interaction that generates an effective chemical potential in a prescribed background. When included in the gravitational action, however, it defines a genuine curvature--matter-coupling variational problem, relevant for the baryon, lepton, and $B\!-\!L$ currents, whether described microscopically by particle-physics operators or macroscopically by a fluid current $J^μ=n_Xu^μ$. Up to a boundary term the interaction is equivalent to $-λR\nabla_μJ^μ$, making its $f(R,{\rm matter})$ character manifest, but the metric equations remain open unless the metric dependence of $J^μ$ is specified. For an arbitrary local realization $J^μ(Ψ,g)$ we derive the universal part of the field equations and isolate the realization-dependent tensor generated by $δ_g J^μ$. In the vector-density realization the explicit $J^α\nabla_αR$ term cancels, but an algebraic term $-λg_{μν}R\nabla_αJ^α$ survives, so the theory admits only a partial effective-Planck-mass interpretation, $M_{\rm eff}^2=M_{\rm Pl}^2-2λ\nabla_μJ^μ$, and a time-dependent effective gravitational coupling during baryogenesis. Specializing to flat FRW with a homogeneous current $J^μ=n_Xu^μ$, we obtain the modified Friedmann and Raychaudhuri equations, the associated continuity relation, and dimensionless diagnostics that quantify when the spectator approximation is controlled. We also discuss the implications for gravitational-baryogenesis studies in modified theories of gravity, providing a consistent GR-side baseline for implementations in both standard cosmology and modified-gravity frameworks.

We introduce a symbolic operator framework for simulating quantum photonic systems that works directly with the canonical commutation relations and the Weyl algebra. Unlike existing Fock-space or Gaussian simulators, our method treats temporal wave packets and polarization modes in a continuous setting and does not rely on discretization or Hilbert-space truncation. Device operations are expressed as algebraic rewrite rules acting on creation and annihilation operators, allowing exact evolution of finite-photon states through linear optical networks. As an illustration, we reproduce Hong-Ou-Mandel interference for Gaussian pulses with controlled temporal and spectral mismatch.

The computational universality with an elementary gate set $\{H,CCZ\}$ can be transformed to the strict universality by using a maximally imaginary state $|+i\rangle$ and some non-imaginary ancillary qubits. From the viewpoint of operational resource theory, it would be intriguing to elucidate a resource for the universality transformation. In this paper, we explore a necessary and sufficient condition for resource states to realize the universality transformation under free real operations. We show that $|+i\rangle$ is a unique resource state up to the free operations. Moreover, we obtain a stronger conclusion. If a given resource state cannot be used for the universality transformation, then realizable quantum gates are restricted to real orthogonal matrices. Therefore, we can tell that $|+i\rangle$ is unique (up to the free operations) not only as a state whose resource measure of imaginarity is maximal, but also as a state which empowers real operations with the ability to apply at least one non-real quantum gate (regardless of the magnitudes of its imaginary parts).

We investigate the precision limits and optimal protocols for sensing single qubit signals in the presence of erasure noise. We study a hierarchy of precision limits achievable with metrological strategies of differing complexity, and identify the optimal protocol for each. The detectability of erasure noise is shown to lead to enhanced precision limits and simplified sensing protocols. For energy gap estimation, we demonstrate that a simple product-state continuous erasure detection strategy yields significant improvements, outperforming optimal entangled protocols even for large numbers of qubits. We show that for other single-qubit signals, quantum error correction provides a substantial advantage by correcting the dominant erasure processes, and can restore Heisenberg-limited precision in certain erasure configurations. As a byproduct of our analysis, we find erasure-conversion schemes for qubits subject to thermal noise that attain the corresponding ultimate precision limits.

We study a class of minimally coupled scalar field theories which leads to analytic solutions for the Hubble rate and the scalar field. The inflationary phenomenology for this class of models can be studied fully analytically. The resulting phenomenology is compatible with the ACT data and for limiting cases, the spectral index is bluer than the ACT constraints and tends to the value $n_{\mathcal{S}}=0.98$, while in the limiting case, the tensor-to-scalar ratio takes very small values, nearly zero. More importantly, the resulting cosmology describes a Universe that has a finite scale factor at $t=0$, thus non-singular, evolves and expands realizing a slow-roll inflationary era and after that it reaches classically a pressure singularity. Classically, the Universe can pass through this singularity, and a turnaround cosmology is realized with the Universe contracting after the turnaround point. However, before the singularity is realized classically, the quantum phenomena dominate the evolution, avoiding the singularity. Specifically we consider the Nojiri-Odintsov conformal anomaly mechanism and we prove that the conformal anomaly erases the classical singular evolution and at the same time it generates extreme particle creation, which eventually reheats the Universe. Thus in this model the scalar field oscillations and the numerous couplings of the inflaton to the Standard Model particles are not required for reheating. In this context, scalar perturbations are enhanced and thus the formation of primordial black holes and the generation of secondary gravitational waves is enhanced. We also discuss several other mechanisms that may lead to the avoidance of the pressure singularity.

In this paper we study the effective gravitational potential of Schwarzschild, Kerr, Reissner-Nordström and Kerr-Newman black holes with the relativistic corrections to evaluate the circular orbits and the Innermost Stable Circular Orbits (ISCOs), a purely relativistic phenomenon, of neutral and charged particles in the vicinity of these spacetimes. We study the circular orbits and ISCOs of black holes mathematically and graphically. Moreover, the astrophysical properties of ISCO are also dealt with. We find that the particle entering the ISCO loses a certain amount of energy, in different spacetimes, by gravitational radiative processes, before finally spiraling into ISCO. The electromagnetic effects of different spacetimes show how the charges sharpen the radius of the circular orbits, and the increasing product of particle and black hole charge increases the radius of the ISCO. Most importantly, the effective potential of the most general spacetime, as predicted by the no-hair conjecture, has been derived.

Relaxation towards equilibrium is often assumed to be slower when a system starts farther from equilibrium, but this intuition fails in the Mpemba effect. Recent advances in controllable quantum platforms have enabled the exploration of its quantum analogue, the quantum Mpemba effect (QME), yet its microscopic origin remains largely unclear. Here we provide a general framework for understanding the QME in closed quantum many-body chaotic systems by reformulating the equilibration process of local subsystems in terms of Ruelle-Pollicott (RP) resonances. We show that suppressing the initial-state overlap with the dominant RP resonant mode accelerates subsystem equilibration and thereby yields the QME. We further uncover that a novel type of strong QME can occur via complete translation-symmetry breaking of initial states. We substantiate our predictions using the prototypical kicked Ising chain and exotic yet experimentally relevant initial states inspired by number theory. These findings cast the QME in closed many-body systems into a unified framework with open-system analogues and provide experimentally accessible signatures on state-of-the-art quantum platforms.

In this letter, we have considered accretion of a particular type of Dark Energy model onto a Schwarzschild type black hole. Before using the model, the free parameters of the Dark Energy model have been constrained with differential ages data. A narrow peak on top of a wide plateau in two parameters' distributions indicates a well defined best fit value embedded within a broad region of near-degenerate solutions. This means the data strongly favours one specific parameter value but also permit a wide range with comparable likelihood. Physically, it reflects that the Dark Energy dynamics are locally constrained yet globally insensitive to small parameter variations. An increasing $\log_{10}\left[M(z)/M_{0}\right]$ since $z=3$ signifies that black holes have continuously grown through accretion and mergers within the standard hierarchical formation scenario. The precise rate of this growth depends on the radiative efficiency $ε$, the effective accretion parameter $λ_{\rm eff}$, and the cumulative impact of merger events.

The frequency-domain waveform emitted by a two-body scattering process is computed in the extreme-mass-ratio limit through the fifth post-Minkowskian (PM) order (i.e., $O(G^5)$) and the fractional sixth post-Newtonian (PN) order. The current accuracy of the scattering waveform obtained by quantum amplitude methods is the one-loop level corresponding to the 3PM order, whereas the 4PM waveform is known up to the 2PN order only as derived within the traditional multipolar-post-Minkowskian formalism. Direct comparison between these waveforms to the first order in the mass ratio shows that they differ at most by the effect of an angular-independent time shift, leading to a complete physical agreement at the same level of accuracy. The new results at the 4PM and 5PM orders thus provide a benchmark for future multiloop calculations. The soft limit of the waveform at the leading order in the small-frequency expansion gives the gravitational wave memory, which is evaluated at the 4PM and 5PM orders. The waveform is also used to obtain the radiated energy at $O(G^6)$, improving its knowledge from 3PN to 6PN level.

We present, to our knowledge, the first systematic study of early-late cosmic evolution and acceleration in the framework of $f(T,\mathcal{L}_m)$ gravity, an extension of teleparallel theories coupling torsion with the matter Lagrangian. By incorporating the Big-Bang Nucleosynthesis (BBN) bound on the freeze-out temperature, we obtain a tight constraint on the inverse-torsion parameter, ensuring consistency with early-time physics. Employing Markov Chain Monte Carlo analyses with progressively richer observational datasets, CC, Union3, and SN22 supernovae, we constrain a well-motivated model and reconstruct key cosmological functions. The reconstructed Hubble and distance modulus functions show excellent agreement with the observations, confirming the observational viability of the model. The model successfully reproduces the observed late-time expansion history, yielding a transition from deceleration to acceleration through the deceleration parameter. The effective equation of state is found to remain negative throughout, with present values $w_0 > -1$, indicating a quintessence-like behavior rather than a cosmological constant or phantom regime. These results highlight the ability of $f(T,\mathcal{L}_m)$ gravity to mimic the concordance scenario while allowing controlled deviations in the expansion history.

We investigate the properties of black hole shadows in the renormalization group (RG) improved Bonanno-Reuter spacetime, incorporating quantum gravitational corrections via the scale-dependent parameter $(\tildeω)$ in a plasma medium. Light propagation in a non-uniform, pressureless plasma with a radial density profile is analyzed through modified equations of motion. The black hole shadow angular radius is computed, and its dependence on $\tildeω$ and the plasma index is analyzed. The analysis of specific limiting cases indicates systematic deviations of the black hole shadow relative to the classical Schwarzschild limit. Using Event Horizon Telescope (EHT) observations of Sgr~A*, we place constraints on $\tildeω$. Furthermore, within the considered parameter range, plasma and quantum-gravity effects exhibit an observational degeneracy, which future high-resolution measurements with the next-generation EHT are expected to break, thereby providing tighter constraints on the model parameters.

We investigate the quantum stability of a timelike topological wormhole with a simple geometry $M_2 \times S^2$, supported classically by anisotropic fluid. We compute the one-loop quantum backreaction generated by the vacuum fluctuations of a minimally coupled, massive scalar field propagating on the wormhole background. Using dimensional regularization we renormalize the one-loop energy-momentum tensor and identify the necessary gravitational counterterms. We then solve the semiclassical Einstein equations to linear order in $\hbar$ for both time-dependent and static metric {\it Ansätze}. Depending on the choice of finite counterterms, the quantum effects can induce either negative or positive angular pressure, which will tend to destabilize or stablize the wormhole, respectively. We also show that a classically traversable wormhole will remain traversable when the quantum backreaction is taken into account. It would be of interest to investigate whether these conclusions remain true when the equations of semiclassical gravity are self-consistently solved.

Quantum sensing and metrology have been intensively studied over the last several decades to surpass the fundamental shot-noise limit of classical systems and approach the Heisenberg limit. However, implementation of N00N-state-based quantum sensing has been severely constrained by the limited order N, intrinsically imperfect fringe visibility, and vulnerability to photon loss. Recently, the coherence de Broglie wavelength (CBW) has been proposed as an alternative method for achieving superresolution in a coherently coupled interferometer architecture, whose characteristics resemble those of photonic de Broglie wavelength (PBW) used in quantum sensing. Here, we experimentally demonstrate scalable CBW superresolution up to N=3, with near-perfect fringe visibility that is invariant to photon loss. The observed CBWs have the potential to enable a superresolution sensing platform even if it remains within the shot-noise limit.

We demonstrate high-fidelity single-qubit gates on a C-shunt flux qubit that simultaneously combines a large anharmonicity ($\mathcal{A}/2π=848~\mathrm{MHz}$) with long relaxation time ($T_1 = 23~μ\text{s}$). The large anharmonicity significantly suppresses leakage to higher energy levels, enabling fast and precise microwave control. Using DRAG pulses and randomized benchmarking, the qubit achieves gate fidelities exceeding 99.9\%, highlighting the capability of C-shunt flux qubits for robust and high-performance quantum operations. These results establish them as a promising platform for scalable quantum information processing.

Valley degrees of freedom in tilted Dirac materials offer a route toward fully electrical quantum control, but previous electrostatic barrier schemes have used the valley index only as a classical filtering resource. Here, we show that a smooth electrostatic barrier operated in a quantum point contact geometry at normal incidence instead realizes coherent valley phase control. In the single-mode regime, both valleys retain near-unit transmission while the tilt-induced valley-dependent traversal phase generates a controllable relative phase shift $Δd = δ_K - δ_{K'}$ between the $|K\rangle$ and $|K'\rangle$ components of the wavefunction. The resulting electrostatic element implements a tunable valley $Z$ rotation whose accessible phase range covers 99.5\% of the full $2π$ interval while maintaining a transmission-balance metric $B$ above 0.99 over a broad parameter window. Combined with a fixed valley-mixing element that supplies an $X$ rotation, this enables universal single-qubit control through a $Z$--$X$--$Z$ Euler decomposition. For realistic parameters, the ballistic gate time is $\sim$50\,fs, with particularly favourable operating windows in 8-$Pmmn$ borophene and WTe$_2$. These results establish tilted Dirac semimetals as a platform for coherent, all-electrical valley manipulation.

Warm inflation is commonly studied under the assumption that the stochastic force sourcing inflaton fluctuations is Markovian. Realistic thermal systems, however, possess finite relaxation times and can therefore generate colored noise with non-zero correlation time. In this work, we investigate warm inflation beyond the Markovian limit and determine how finite correlation time modifies the primordial scalar power spectrum. We show that memory effects suppress the scalar spectrum relative to the standard white-noise result and derive a simple expression for this correction in terms of the background thermal dynamics. In particular, we relate the size of the non-Markovian effect directly to the thermal ratio between the bath temperature and the Hubble scale, thereby establishing a transparent link between warm-inflation background quantities and the validity of the Markovian approximation. We also derive the corresponding modification of the tensor-to-scalar ratio and the induced shifts in the scalar spectral index and the running of the scalar spectral index. Our results provide a simple and practical diagnostic for identifying when finite correlation-time effects become relevant in warm-inflation model building.

We analyze a global sensing scenario in which the frequency of a monitored Kerr parametric oscillator is estimated assuming limited prior information. The frequency is estimated in real-time by continuously monitoring the oscillator quadrature through homodyne detection and processing the resulting photocurrent with an extended Kalman filter (EKF). Due to the sensor nonlinearity, individual EKF trajectories do not always converge to the true unknown frequency in the long-time limit. However, we show that the statistical distribution of the frequency estimates does exhibit a sharp peak around the true value in the same limit. Leveraging this key statistical property, we develop a global sensing protocol assisted by adaptive control of the sensor parameters to harness critical enhancement. We present numerical evidence that this criticality-enhanced frequency estimation remains robust under low detection efficiency.

We prove that autoparallel curves associated with a torsion-free but not necessarily metric-compatible affine connection can be derived from an action principle. We explicitly construct the action functional and show by standard variational techniques that it produces the desired equations. Our analysis is based on systematically solving the inverse problem of the calculus of variation and the associated Helmholtz conditions. This demonstrates that the dynamics of autoparallels admit a consistent variational formulation even in the presence of non-metricity. Our results provide a variational framework for particle motion in metric-affine geometries and thereby contribute to the mathematical foundations of the geodesic principle in relativistic gravity.

Linear quantum amplifiers are indispensable tools for quantum technologies, yet their performance is fundamentally limited by quantum noise, precluding any signal-to-noise ratio (SNR) enhancement unless supplemented by post-selection or non-classical resources. To surpass this limitation, we propose a nonlinear quantum amplification strategy that exploits the interplay between a gain-stabilized bright eigenmode of a coupled two-mode bosonic system and Kerr nonlinearity. We demonstrate that this interplay enables the signal gain to surpass the noise gain in a selected quadrature, leading to a net increase in the SNR beyond the quantum limits of conventional linear amplifiers. Our work thus establishes a novel nonlinear amplification paradigm capable of enhancing the SNR, with promising applications across quantum information processing, quantum communications, and quantum metrology.

Transport engineering has significant potential to benefit from quantum computing. The rise of intelligent transport systems, autonomous vehicles, and the Internet of Things has created an unprecedented demand for efficient information processing and computational optimisation. Accordingly, transport engineers and scientists have explored the ever-improving capabilities of quantum computers in an effort to meet this demand. Motivated by this growing interest, this paper sets out four aims: (1) to introduce the fundamental aspects of quantum computing relevant to the transport domain, (2) to identify transport-related problems which are suitable for quantum acceleration, (3) to develop a pipeline for solving these problems, and (4) to provide a systematic review of the existing literature. For the latter, a systematic search of the Scopus database (and supplemented by additional citation sources) identified 103 studies for inclusion following PRISMA 2020 guidelines. While a diverse set of use cases have been proposed, we conclude that future research should prioritise problems where quantum computation offers a clear practical benefit. To this end, we suggest promising directions to guide further work in this burgeoning subfield.

We formalize the concept of subtime -- a reversible mode of information interchange within entangled systems -- and show how classical time emerges as an asymptotic limit through decoherence. Building on the photon clock model, in which a single photon confined between two ideal mirrors creates an alternating causality regime, we develop a process-theoretic formalization using the Oreshkov--Costa--Brukner framework extended with an explicit time-reversal duality condition. We introduce Perfect Information Feedback (PIF) as the information-theoretic realization of this reversibility, demonstrating that mutual information is conserved in any closed causal loop and that entropy quantifies the degree of unreflected causality. We define the Reversible Causal Principle (RCP): every causal relation possesses a conjugate dual, and entropy, energy dissipation, and the classical arrow of time appear only when these alternating components decohere or fail to reflect perfectly. The framework unifies Wheeler--Feynman absorber theory, Bennett's reversible computation, Shannon's communication theory, and the process matrix formalism under a single symmetry principle, and identifies experimentally accessible signatures in reversible digital links and quantum switch experiments. The arrow of time, in this picture, records the universe's imperfect causal echo.

We present a general-purpose quantum error correction primitive based on state purification via the SWAP test, which we refer to as purification quantum error correction (PQEC). This method operates on $N$ noisy copies, requires minimally $O(M\log_2 N)$ data qubits to process the $M$-qubit inputs. In a similar way to standard QEC, the purification steps may be interleaved within a quantum algorithm to suppress the logical error rate. No postselection is performed and no knowledge of the state is required. We analyze its performance under a variety of error channels and find that PQEC is highly effective at boosting fidelity and reducing logical error rates, particularly for the depolarizing channel. Error thresholds for the local depolarizing channel are found to be $ 75 \%$ for any register size. For local dephasing, the error threshold is reduced to $ 50 \% $ but may be boosted using twirling.

Thermodynamic geometry provides a powerful framework for probing the microscopic structure of thermodynamic systems. Among its formulations, Geometrothermodynamics (GTD) has been widely applied to black hole thermodynamics, owing to its Legendre-invariant construction. However, recent work has shown that conventional GTD metrics fail to encode essential physical boundaries of thermodynamic phase space. By modifying the conventional metric structure, three new GTD metrics were previously introduced, which successfully capture these boundaries in regular spacetime. Whether such modified metrics remain viable in different spacetime backgrounds and under changes of thermodynamic ensemble has remained an open question. In this work, I address this issue by investigating the behavior of modified GTD metrics in AdS spacetime and across different thermodynamic ensembles in the framework of Bardeen AdS black hole. An analysis of thermodynamic geodesics demonstrates that the modified GTD metrics consistently respect physical boundaries of the phase space, in contrast to conventional GTD metrics. This behavior is preserved in AdS spacetime and under Legendre transformations, establishing the robustness and universality of the modified GTD framework.

This work presents a quantum mechanical framework for analyzing quantization-based optimization algorithms. The sampling process of the quantization-based search is modeled as a gradient-flow dissipative system, leading to a Hamilton-Jacobi-Bellman (HJB) representation. Through a suitable transformation of the objective function, this formulation yields the Schroedinger equation, which reveals that quantum tunneling enables escape from local minima and guarantees access to the global optimum. By establishing the connection to the Fokker-Planck equation, the framework provides a thermodynamic interpretation of global convergence. Such an analysis between the thermodynamic and the quantum dynamic methodology unifies combinatorial and continuous optimization, and extends naturally to machine learning tasks such as image classification. Numerical experiments demonstrate that quantization-based optimization consistently outperforms conventional algorithms across both combinatorial problems and nonconvex continuous functions.

Simulating real-time dynamics under a Hamiltonian is a central goal of quantum information science. While numerous Hamiltonian-simulation quantum algorithms have been proposed, the effects of physical noise have rarely been incorporated into performance analysis, despite the non-negligible noise levels in quantum devices. In this work, we analyze noisy Hamiltonian simulation with quantum error mitigation for Trotterized and randomized LCU-based Hamiltonian simulation algorithms. We give an end-to-end comprehensive complexity analysis of error-mitigated Hamiltonian simulation algorithms using the mean-squared error. Because quantum error mitigation incurs an exponential cost with the number of layers in quantum algorithms, there is a trade-off between the sampling cost and the bias in simulation accuracy or the algorithmic sampling overhead. Optimizing this trade-off, we derive an analytic depth-selection rule and characterize the optimal end-to-end scaling as a function of target accuracy and noise parameters. We further quantify the noise-characterization cost required for error mitigation via gate set tomography and the recently proposed space-time noise inversion method, showing that the latter can significantly reduce the characterization overhead.

We establish a framework for realizing back-action-evading (BAE) measurements and quantum non-demolition (QND) variables in linear quantum systems. The key condition, a purely imaginary Hamiltonian with a real or imaginary coupling operator, enables BAE measurements of conjugate observables. Symmetric coupling further yields QND variables. For non-compliant systems, coherent feedback can engineer BAE measurements. Crucially, the QND interaction condition simultaneously ensures BAE measurements and promotes the coupling operator to a QND observable. This work provides a unified structural theory for enhancing precision in quantum metrology and sensing.

Scalar--tensor theories of gravity provide a natural extension of general relativity and may predict naked singularities as alternative compact objects. In this work, we investigate a novel exact solution within the Freud--Nambu scalar--tensor gravity framework, generalizing the Janis--Newman--Winicour (JNW) naked singularity spacetime through the introduction of a parameter $q$ coupled to a real scalar field $\varphi$ with mass $μ$. Although the metric remains identical to the JNW solution, the scalar field profile is modified, providing a parametrized deformation of this class of spacetimes. We analyze particle dynamics in this background, including a direct linear coupling between the test particle and the scalar field characterized by the parameter $g_s$. The influence of these parameters on astrophysical observables is studied through the specific angular momentum, the innermost stable circular orbit (ISCO), and the radiative efficiency of accretion. We also derive the epicyclic frequencies governing oscillatory motion and explore their implications for quasi-periodic oscillations (QPOs) in black hole binaries. Within the epicyclic resonance model, the upper and lower QPO frequencies depend sensitively on the parameters $n$, $g_s$, and $q$. To constrain the model, we perform a Markov Chain Monte Carlo analysis using twin-peak QPO data from the microquasars XTE~J1550--564 and GRS~1915+105. The resulting black hole masses agree with previous estimates and provide the first observational constraints on the parameters $q$ and $g_s$, indicating that modified gravity effects may leave detectable imprints on strong-field astrophysical phenomena.

As a key index for applications of a device, the device's reliability is its ability to survive (function normally over time) under the influence of some environment. In this paper we present a quantum energy-storing device model with a quantum spin chain, whose environment influence is described by the Lindblad master equation. Here the device survives if the spin system stays in the state with nonzero excitations; otherwise, it fails. Because the Lindblad dynamics enforces one-way energy decay and strict irreversibility of the failure state, we can investigate the reliability of the quantum device directly using classical reliability theory. Focusing on the minimal nontrivial case -- a two-site spin-1/2 chain -- we derive closed-form expressions for the reliability and the hazard rate. The dynamics exhibit an overdamped-underdamped crossover controlled by the competition between coherent exchange and dissipation inhomogeneity. The exact analytical formulas are in excellent agreement with numerical simulations. More importantly, we establish an experimentally accessible protocol for assessing reliability based on first-passage time statistics.

Superconducting quantum circuits are promising platforms for scalable quantum computing, where qubit coherence is critically determined by microscopic defects in the oxide tunneling barrier of Josephson junctions. Amorphous Al$_2$O$_3$ is widely used as a barrier material, but under irradiation, oxygen vacancy (V$_O$) defects are readily generated, introducing noise sources that accelerate qubit decoherence. We systematically investigate the structural characteristics and electronic impact of V$_O$ defects in amorphous Al$_2$O$_3$ using first-principles calculations and \textit{ab initio} molecular dynamics. Our results show that both the coordination environment and concentration of V$_O$s strongly influence electrical conductivity. In particular, two- and three-coordinated V$_O$s, unique to the amorphous structure, enhance conductivity more than conventional four-coordinated vacancies. Increasing V$_O$ concentration amplifies conductivity fluctuations, which we link to critical current noise in Josephson junctions. Using a noise model, we estimate that higher V$_O$ densities lead to shorter qubit coherence times. These findings provide insights for radiation-hard design of superconducting quantum devices.

ArXiv:2508.17898 proposed the booklet wormhole as the holographic dual of the GHZ state. This paper extends the investigation into this geometry, particularly focusing on the junction conditions for matter fields. We show that the symmetry of the GHZ state requires the bulk to admit special Killing vector fields that standard manifolds cannot realize. Moreover, these bulk symmetries require unprecedented quantum non-local junction conditions at the multi-way interface: Observers entering from different horizons will perceive different states inside the wormhole, where the junction conditions appear as constraints on the observables of different sets of observers. We finally discuss how to render booklet wormholes traversable via boundary deformations. A localized wave packet injected from one page generally evolves into a non-local mixed state on each remaining page, with the information encoded in the entanglement between different pages.

We construct unclonable encryption (UE) in the Haar random oracle model, where all parties have query access to $U,U^\dagger,U^*,U^T$ for a Haar random unitary $U$. Our scheme satisfies the standard notion of unclonable indistinguishability security, supports reuse of the secret key, and can encrypt arbitrary-length messages. That is, we give the first evidence that (reusable) UE, which requires computational assumptions, exists in "micocrypt", a world where one-way functions may not exist. As one of our central technical contributions, we build on the recently introduced path recording framework to prove a natural ``unitary reprogramming lemma'', which may be of independent interest.

Plasma, as a medium around the black hole for light propagation, is known to visibly alter the shape of its shadow and the observables, which could impact the interpretation of the Event Horizon Telescope results. In this study, we examine how dark matter and non-magnetized, pressureless plasma influence the shadow of a Kerr-like black hole. We analyze the null-geodesics in the presence of both homogeneous and inhomogeneous plasma profiles and show how their influence on photon orbits affects the resulting black hole shadow. Our findings indicate that increasing the black hole's spin generally enlarges both the shadow radius and deformation. Additionally, the viewing angle decreases the shadow radius while reducing deformation as the observer moves farther from the equatorial plane. For this model, astrophysically reasonable amounts of dark matter show no significant impact on the photon trajectories. However, we observe that increasing plasma density increases both the shadow radius and deformation for homogeneous plasma, while it decreases them for inhomogeneous plasma. The emission rate also depends significantly on the model of plasma chosen, with homogeneous plasma causing significantly more emission as plasma strength increases. We also study the constraints obtained from comparing theoretical shadow radii with EHT observations of M87* and Sgr A*, which allows us to infer reasonable plasma distribution properties and frequencies in our theoretical model.

We formulate an analogue transverse first law for general closed marginally trapped surfaces in arbitrary spacetimes. The construction is intrinsically quasi-local and is attached directly to an individual marginally trapped surface, rather than to a preferred horizon worldtube. Taking the Hawking energy as the internal energy and an invariant effective surface gravity associated with the marginally trapped surface as the quantity controlling the thermal term, we derive a balance law in which the variation of energy splits into a generalized heat contribution and a total work contribution. In this way, the resulting law provides a codimension-two, transverse counterpart to existing horizon-based formulations of black-hole thermodynamics. We show that the formalism reproduces the expected results for round spheres in spherically symmetric spacetimes. We then examine semiclassical equilibrium and evaporating regimes, and extend the analysis to non-spherically symmetric marginally trapped surfaces in Kerr. These examples indicate that the framework remains applicable in situations where a horizon-based treatment is either nonunique or technically cumbersome, and suggest that closed marginally trapped surfaces provide a natural arena for a genuinely quasi-local thermodynamics of black holes.

Spatially varying near-horizon fluctuations of temperature of a Schwarzschild Black Hole is considered within the Euclidean Gravity approach. We present evidence that suggests that such fluctuations in temperature are closely related with the near-horizon supertranslations of a Black Hole. This allows one to express the temperature-dependent corrections to the near-horizon part of the Euclidean gravity action or Free energy in terms of polynomial functionals of a near-horizon supertranslation for all orders. The leading order term turns out to be proportional to the near-horizon supertranslation charge. We also show that this same term results from a microscopic state counting of the near-horizon field with the constraint of near-horizon supertranslation symmetry imposed. The constructed near-horizon partition function provides a physically appealing intuitive description of the low-energy part of horizon physics as a sum over all possible near-horizon supertranslations. This suggests a dual description of near-horizon physics in terms of alternate variables. The implications of these results are briefly discussed.

Autonomous agents, especially delegated systems with memory, persistent context, and multi-step planning, pose a measurement problem not present in stateless models: an agent that preserves continued operation as a terminal objective and one that does so merely instrumentally can produce observationally similar trajectories. External behavioral monitoring cannot reliably distinguish between them. We introduce the Unified Continuation-Interest Protocol (UCIP), a multi-criterion detection framework that moves this distinction from behavior to the latent structure of agent trajectories. UCIP encodes trajectories with a Quantum Boltzmann Machine (QBM), a classical algorithm based on the density-matrix formalism of quantum statistical mechanics, and measures the von Neumann entropy of the reduced density matrix induced by a bipartition of hidden units. We test whether agents with terminal continuation objectives (Type A) produce latent states with higher entanglement entropy than agents whose continuation is merely instrumental (Type B). Higher entanglement reflects stronger cross-partition statistical coupling. On gridworld agents with known ground-truth objectives, UCIP achieves 100% detection accuracy and 1.0 AUC-ROC on held-out non-adversarial evaluation under the frozen Phase I gate. The entanglement gap between Type A and Type B agents is Delta = 0.381 (p < 0.001, permutation test). Pearson r = 0.934 across an 11-point interpolation sweep indicates that, within this synthetic family, UCIP tracks graded changes in continuation weighting rather than merely a binary label. Among the tested models, only the QBM achieves positive Delta. All computations are classical; "quantum" refers only to the mathematical formalism. UCIP does not detect consciousness or subjective experience; it detects statistical structure in latent representations that correlates with known objectives.

In this manuscript, we develop a class of inhomogeneous relativistic cosmological models with the following properties: (i) They contain cosmological observers to whom the spatial geometry and the expansion are homogeneous and isotropic; (ii) Matter behaves closely to dust, as it is formed by an ensemble of massive particles whose number density $4$-vector is conserved and reacts viscously to the local tidal forces; (iii) They generalize the dust FLRW model; (iv) They give rise to effective models on large scales that reproduce the FLRW behaviour with dust and a running dark-energy term, which appears as a backreaction effect from the local gravitational potential; (v) By suitably setting the distribution of energy in the current universe, the effective large-scale equation-of-state parameter can reproduce, in principle, any polynomial on the scale factor whose term of order zero is negative; (vi) The luminous distance observations imply an apparent large-scale acceleration, as the deceleration parameter is negative, while in reality the universe is decelerating.

We start a systematic investigation of the one-loop mass corrections to (super-)string massive higher-spin states. While the imaginary part of the relevant amplitudes are finite, being related to the width of the decay of the states into two lower-mass states at tree level, the real part is generally IR-divergent and needs regularization and renormalization. We mostly focus on states of the first Regge trajectory in the NS-NS sector of Type-II string theories. We explicitly derive a closed-form expression for the integral over the insertion point, relying on properties of elliptic functions and lattice sums. We then regularize the IR divergent integral over the modular parameter of the torus, applying the $i\varepsilon$-prescription in string theory. As a result we compute the desired mass corrections up to level $N = 10$ and analyze their behavior at increasing $N$. Finally, we speculate on the existence of mixing among lower-spin states and conjecture that the one-loop mass matrix be governed by some random matrix theory.

We investigate the geodesic structure, shadow, thermodynamics, and Hawking radiation from a dyonic ModMax black hole in Kalb-Ramond gravity with a cloud of strings. The combined presence of ModMax nonlinear electrodynamics, the Lorentz-violating Kalb-Ramond background, and the string cloud breaks asymptotic flatness and introduces a global conical deficit that modifies all observables through a single geometric prefactor. We derive analytic expressions for the photon sphere, critical impact parameter, and shadow radius, and show that the shadow size depends on both the non-flat asymptotics and the ModMax screening of the dyonic charge. For massive test particles, we determine the innermost stable circular orbit and the accretion efficiency as functions of all model parameters. We also establish the first law of black hole thermodynamics and the generalized Smarr relation for this solution, identify a Hawking-Page-type phase transition in the specific heat, and compute the spectral energy emission rate, which we show is directly governed by the shadow radius in the geometric-optics limit. Our results demonstrate that the interplay of these three ingredients produces a phenomenology observationally distinguishable from standard Reissner-Nordström black holes.

We analyze many-body entanglement in interacting fermionic systems by using the $M$-body reduced density matrix. We demonstrate that if a particle number conserving fermionic Hamiltonian contains only up to $M$-body interaction terms, then its $N$-particle ground state cannot be maximally $M$-body entangled. As a key step in the proof, we show that the energy expectation value of a maximally $M$-body mixed state is equal to the spectral mean of the Hamiltonian on the corresponding $N$-particle subspace. We further demonstrate that the many-body entanglement structure of a ground state can place quantitative constraint on the interaction strength of its parent Hamiltonian. We illustrate the theorem and its implications in Hubbard and extended SYK models. Going beyond ground states, we analyze entanglement generation under unitary dynamics from Slater-determinant initial states in these models. We determine early-time growth and estimate entanglement saturation times. Finally, we derive explicit symmetry-refined saturation upper bounds for $M$-body entanglement in the presence of an Abelian symmetry.

We revisit quantum tomography of $H \to ZZ$ and $H \to WW$ in the presence of higher-order corrections. We verify that neither the use of an effective spin analysing power (only for $ZZ$) or a photon veto are sufficient to render the naively-constructed spin density operators physical. A subtraction of higher-order corrections is thus necessary to perform consistent quantum tomography. Such corrections are small when compared to expected experimental uncertainties with current data. As a by-product, we point out the striking possibility to observe parity-violating effects in $H \to WW$.

Quantum error mitigation (QEM) and quantum error correction (QEC) are two research areas that are often considered as distinct entities, and the problem of combining the two approaches in a non-trivial way has only recently started to be explored. In this paper, we explore a paradigm at the intersection of the two, based on the error mitigation technique of Zero-Noise Extrapolation (ZNE), that uses the distance of an error correcting code as a noise parameter. This is distinct from some alternative approaches, as QEC is here used as a subroutine inside the QEM framework, while other proposals use QEM as a subroutine inside QEC experiments. Intuitively, we exploit the fact that a reduction in the physical noise level is analogous to an increase in the code distance, as both of them result in a decrease in the logical error rate. As such, the extrapolation to zero noise in the case of ZNE becomes comparable to the extrapolation to infinite distance in the case of this method. We describe how to calculate expectation values from a fault-tolerant computation, and we gain some analytical intuition for our ansatz choice. We explore the performance of the considered method to reduce the errors in a range of expectation values for a realistic circuit-level noise model and realistic device imperfections on the rotated surface code, and we particularly show that the performance of the method holds even in the case of non-stabiliser input states.

We study the fundamental limitations of implementing time-dependent Hamiltonian protocols when ''time'' is provided by a quantum clock rather than an external classical parameter. For a parametric harmonic oscillator controlled through a shortcut-to-adiabaticity (STA) schedule and coupled to a minimal clock degree of freedom, tracing out the clock yields an effective reduced dynamics that is a mixture of unitary Gaussian trajectories. Within a noise-dominated regime, we compute the energetic deviation from the target STA outcome and its fluctuations, together with the fidelity to the target evolution and the purity loss of the reduced state, for vacuum and coherent initial states. Combining these observables produces a thermodynamic-uncertainty-type tradeoff that links achievable precision to an irreducible loss of purity set by the clock precision and the protocol sensitivity.

We explore the sense in which the existing constructions for higher-order maps on quantum theory based on causality constraints and compositionality constraints respectively, coincide. More precisely, we construct a functor F : Caus(C) -> StProf(C1) from higher-order causal categories to the category of strong profunctors over first-order causal processes that is lax-lax duoidal, full, faithful, and strongly closed whenever C is additive. When C = CP this embedding is furthermore strong on the sequencer for duoidal categories, expressing the possibility to interpret one-way signalling (but not general non-signalling) constraints in terms of the coend calculus for profunctors. We conclude that insofar as compositional constraints can be used to express causality constraints, the profunctorial approach generalises higher-order quantum theory to a construction over general symmetric monoidal categories.

Efficient low-energy state preparation is a key objective in quantum computation and quantum simulation. Quantum imaginary-time evolution replaces real-time dynamics with imaginary-time dynamics, exponentially suppressing higher-energy eigenstates. We introduce deterministic unitary protocols that approximate imaginary-time evolution for ground-state preparation. The protocols require multiple copies of the system, real-time evolution under the system Hamiltonian, and controlled-SWAP operations (or more general SWAP-generated unitaries). We analyze two concrete circuit families: a tree architecture with provable polynomial-in-depth convergence but rapidly growing width, and a compact "hedge" architecture that achieves comparable accuracy with only polynomial width in a heuristic construction supported by numerics. We provide numerical evidence that mid-circuit post-selection can accelerate convergence with practical success probabilities. Separately, we demonstrate that circuit volume can be traded for the shot complexity of post-circuit observable estimation in the ground-state preparation setting. We outline concrete implementation of platform-specific routes, where multi-copy registers and SWAP-mediated couplings are natural, thereby illustrating how these hybrid analog-digital circuits can complement existing state-preparation methods in the near term.

The Kraus representation of quantum channels allows for a precise emulation of the complex dynamics that take place on quantum processors, whether for benchmarking algorithms, predicting the performance of error correction and mitigation, or in the myriad other uses of compiled digital sequences. Nonetheless, starting from first principles to obtain continuous quantum master equations involves various approximations such as weak coupling to the environment. Further, converting these equations to Kraus operators cannot generally be obtained in closed-form due to the complicated commutator structure of the problem. In our work, we bridge this gap by providing a general closed form formulation for arbitrarily strong driving while remaining linear in the dissipator. The Kraus solution is expressed as a Riemann sum where higher terms can converge quickly to high precision, which we demonstrate numerically. Such a formulation is highly relevant to quantum computing and gate-based models, where effective models are highly sought for large rotation gate angles, even under the influence of underlying non-trivial noise mechanisms.

We present near-term quantum algorithms for auxiliary-field quantum Monte Carlo (AFQMC), which represents imaginary-time projection for ground-state calculation as an ensemble of one-body propagators driven by stochastic fields $Ω$. Starting from the Feynman-Kac formula, we convert each trajectory into a sequence of piecewise-constant one-body generators using stochastic Magnus expansions up to second order, and embed the resulting nonunitary slices into unitaries with a small ancilla overhead. This lifts the projector dynamics to a unitary evolution, enabling coherent circuit execution in the regime $\|Ω\| τ=O(1)$ and reducing the need for frequent mid-circuit measurement. We further derive an equivalent linear-combination-of-unitaries (LCU) form that yields system-only, shallower circuits by trading ancilla cost for additional trajectory sampling. Benchmarks on the Hubbard model verify the accuracy of the dilation and Magnus expansions classically and demonstrate multi-step executions on IBM quantum hardware.

Quantum generative modeling has emerged as a promising application of quantum computers, aiming to model complex probability distributions beyond the reach of classical methods. In practice, however, training such models often requires costly gradient estimation performed directly on the quantum hardware. Crucially, for certain structured quantum circuits, expectation values of local observables can be efficiently evaluated on a classical computer, enabling classical training without calls to the quantum hardware in the optimization loop. In these models, sampling from the resulting circuits can still be classically hard, so inference must be performed on a quantum device, yielding a potential computational advantage. In this work, we introduce a photonic quantum generative model built on parametrized Gaussian Boson Sampling circuits. The training is based on the efficient classical evaluation of expectation values enabled by the Gaussian structure of the state, allowing scalable optimization of the model parameters through the maximum mean discrepancy loss function. We demonstrate the effectiveness of the approach through numerical experiments on photonic systems with up to 805 modes and over a million trainable parameters, highlighting its scalability and suitability for near-term photonic quantum devices.

Quantum many-body scars are nonthermal states exhibiting persistent revivals in an otherwise ergodic, nonintegrable quantum system. Existing examples of exact quantum scars, however, have not yet been amenable to direct experimental demonstration. Here we show that a minimal model of hardcore bosons hopping on a $π$-flux ladder is sufficient to give rise to an exact scar due to kinetic frustration. The simplicity of this model makes it suitable for multiple existing quantum simulation platforms, which we illustrate with proposals for cold atom Bose-Hubbard simulators and polar molecule or Rydberg atom tweezer arrays. In these platforms, the scar lifetime can be extended by tuning experimentally accessible parameters, like the Hubbard interaction or a Floquet drive. Finally, we introduce a practical heuristic based on the energy distribution of eigenstates for systematically predicting and optimizing quantum many-body scar lifetimes. Their cross-platform realizability and long lifetimes make them well-suited for benchmarking coherence and exploring nonergodic dynamics in current and near-term quantum devices.

Dielectric loss at the interfaces of superconducting films has long been recognized as limiting the performance of state-of-the-art superconducting circuits. Notably, the presence of a native oxide layer on the film is hypothesized to contribute to dielectric loss at the metal-air interface. Here, we explore rhenium as a candidate for the film, motivated by its remarkable property to suppress native oxide formation. We demonstrate rhenium on sapphire as a promising material platform for superconducting circuits through the realization of transmons with mean relaxation times $T_1$ up to 407 microseconds at 5 GHz. Our transmons are supplemented with a loss characterization study, in which we separate the dominant loss mechanisms and construct a loss budget that agrees with our $T_1$ measurements. Further characterization may establish rhenium as a leading candidate for maximizing decoherence time.

Neural quantum states (NQS) provide a flexible variational framework for many-body wavefunctions, but suffer from high computational cost and limited interpretability. We introduce DysonNet, a broad class of NQS that couples strictly local nonlinearities through global linear layers. This structure is analogous to a truncated Dyson series which gives an intuitive interpretation of local wavefunction updates as scattering from static impurities. By resumming the scattering series, single-spin-flip updates can be computed in $\mathcal{O}(1)$ time, independent of system size, using an algorithm we call ABACUS. Implementing DysonNet with the state-space model S4, we obtain up to $230\times$ speedups over Vision-Transformers for computing the local estimator. This corresponds to an asymptotic $\mathcal{O}(N^2)$ improvement in training-time scaling, reaching $\mathcal{O}(N \log^2 N)$ total training complexity in area-law phases. Benchmarks on the 1D long-range Ising model and frustrated $J_1$-$J_2$ chains show that DysonNet matches state-of-the-art NQS accuracy while removing the dominant local-update overhead. More broadly, our results suggest a route to scalable NQS architectures where physical interpretability directly enables computational efficiency.

The quantum Hall effect and the quantum anomalous Hall effect both require time-reversal invariance to be broken. We show that non-equilibrium effects can cause Hall physics to arise even when the system is weakly time-reversal symmetric and no magnetic field is applied. In our model, this occurs due to a fermionic subsystem breaking time-reversal invariance even if the system as a whole does not. The fermions receive a TRI-breaking self-energy, caused by interactions with bosonic degrees of freedom in the system and with an external reservoir. As a result, the fermions develop a non-quantized Hall conductance. We demonstrate that, unlike in the equilibrium case, the presence of a mass term is insufficient for the Hall conductance to appear, and wave-function renormalization effects have to be included.

Advancing quantum technologies requires precise and robust coherent control of quantum systems. Robust higher-order Hamiltonian engineering is essential for high-precision control and for accessing effective dynamics absent at zeroth order. Here, we introduce a systematic methodology for achieving the precision, robustness, and complexity required for quantum control through the engineering of higher-order processes and effective Hamiltonians. We identify the minimal subspace of achievable effective Hamiltonian at each order and provide universal cost functions for achieving desired targets. Examples include robust sequences for decoupling, three-body interactions and detuning/interaction correlations.

We present an efficient protocol leveraging classical computation to support Initial State Preparation for strongly correlated fermionic systems, a critical bottleneck for fault-tolerant quantum simulation. Focusing on nuclear shell model eigenstates, we first demonstrate that the Density Matrix Renormalization Group algorithm can efficiently approximate target states as Matrix Product States, capitalizing on the favourable entanglement structure of these fermionic systems. These high-fidelity approximations are then leveraged as a classical resource in a variational circuit optimization scheme to compile shallow quantum circuits. We establish concrete resource estimates by decomposing the resulting circuits into the industry-standard Clifford$+T$ gateset, exploring the benefits of specialized $U3$ synthesis techniques. For all nuclear systems tested, on up to 76 qubit Hamiltonians, we consistently find low $T$-count circuits preparing the nuclear eigenstates to high fidelity with $\sim 2\times 10^4$ total $T$ gates. This low number gives confidence these eigenstates can be prepared on early fault-tolerant quantum computers. Our work establishes a viable path toward practical ground state preparation for nuclear structure and other fermionic applications.

Time-slicing has emerged as a strategy for incorporating semiclassical propagation into real-time path integral formulation and recovering full quantum dynamics. A central step is the decomposition of a time-evolved wave function into a superposition of Gaussian wave packets (GWPs). Here we introduce a quadrature-free variational framework for GWP decomposition, reformulating it as an optimization problem in which the GWP parameters are chosen to maximize the overlap with the time-evolving wave function. An autoencoder-decoder neural network is used for this optimization, with the representation being adaptively reoptimized during propagation. Each wave packet in this decomposition represents a localized patch of the underlying semiclassical manifold, while retaining full correlations between all degrees of freedom. This variational adaptive Gaussian decomposition (VAGD) approach yields a compact Gaussian expansion, providing a scalable route to time-sliced semiclassical quantum dynamics. While general, applying VAGD to facilitate time-slicing of thawed Gaussian approximation (TGA) allows a route to improving the semiclassical treatment to the full quantum mechanical result in a systematic manner.

A large scale quantum network requires the distribution of high-fidelity end-to-end entanglement. To overcome the range limitations inherent to terrestrial fiber, a leading architecture has emerged: satellite-based sources transmitting entanglement to quantum repeaters on the ground. By bridging the gap between abstract analytical frameworks and computationally heavy numerical simulations, this paper provides the first quantitative answer to the question of such a network's achievable performance with current and near-term space technology, while accounting for entanglement swapping and purification. This is achieved by integrating a detailed physical model of a satellite-to-ground link into an analytical entanglement resource estimation framework for quantum repeaters, enabling an optimization of the end-to-end entanglement rate. Our analysis, performed across leading quantum hardware platforms, shows that Low Earth Orbit satellite constellations combined with quantum repeaters employing Neutral Atom or Nitrogen and Silicon Vacancy qubits, could enable a global quantum network, distributing entanglement over distances up to 20,000 km, sufficient for connecting any two points on Earth. This work highlights the major bottlenecks in space and quantum hardware technologies, which need to be addressed, thereby guiding informed investments necessary for enabling a large scale quantum network.

Is spacetime fundamental or can it be derived through quantum interactions? We propose here a way to describe time dilation solely from quantum mechanics. First we start by observing that any operational notion of time must imply some sort of regular motion and, crucially, some sort of memory. Thus the clock model we use here is a simple example of a quantum memory: a charging battery. We describe here the charging of such batteries with quantum open dynamics. The expected value of these batteries grow linearly in time like $\varphi t$. The open dynamics is dependent on an auxiliar state $σ$. Therefore, with a different auxiliar state we have a different $\varphi$. We can describe time dilation and thus a metric this way. We exemplify with a black hole metric.

Quantum annealing approximately solves combinatorial optimization problems by leveraging the principles of adiabatic quantum systems. In this approach, the system's Hamiltonian evolves from an initial general state to a problem-specific state. This study introduces multi-tasking quantum annealing (MTQA), a method that enables the parallel processing of multiple optimization problems by embedding them into spatially distinct regions on quantum hardware. MTQA is evaluated using two NP-hard problems: the minimum vertex cover problem (MVCP) and the graph partitioning problem (GPP). This parallel approach optimizes quantum resource utilization by concurrently utilizing idle qubits. The findings demonstrate that MTQA achieves a solution quality comparable to single-problem quantum annealing and classical simulated annealing (SA), while notably reducing the time-to-solution (TTS) metrics. Eigenspectrum analysis further theoretically supports the hypothesis that parallel embedding preserves quantum coherence and does not increase computational complexity by efficiently utilizing available quantum hardware (e.g., qubits and couplers). MTQA enables efficient multitasking in quantum annealing, optimizing hardware utilization and improving throughput for concurrent tasks and demonstrating performance for problems up to 100 nodes in real-world applications.

Quantum algorithms present a quadratically improved complexity over classical ones for certain sampling tasks. For instance, the Quantum Amplitude Estimation (QAE) algorithm promises to speedup the estimation of the mean of certain functions, given access to the quantum state corresponding to the probability distribution to be sampled from. Classically, samples are often obtained by running steps a Markov Chain. In this work, we experimentally use encodings of Markov chains to prepare quantum states and run a quantum Markov Chain Monte Carlo algorithm (qMCMC) on Quantinuum's H2 and Helios quantum computers. We demonstrate that it is possible to obtain accurate results on current Noisy Intermediate Scale Quantum (NISQ) hardware, operating directly on the physical qubits.

We present new relations derived from Noether's identity that reveal the compatibility between the components of the Hessian matrix of the Lagrangian, the infinitesimal symmetry transformation of the configuration variables and time, and a constant of motion. Using these relations, we develop two new methods to incorporate symmetry requirements directly into the inverse problem of mechanics, thereby restricting the set of acceptable Lagrangians. We accomplish this by combining these relations with Helmholtz's conditions, which allow us to construct Lagrangians whose actions exhibit specific symmetries from the outset. The theory is illustrated with one- and two-dimensional examples.

We report the first lattice QCD study of the $s$-wave scattering of the $\bar{D}$-meson and the nucleon at the physical point, utilizing (2+1)-flavor configurations generated by the HAL QCD collaboration with a pion mass of $m_π\simeq 137$ MeV and a lattice spacing of $a\simeq0.084$ fm. By applying the HAL QCD method to the four-point correlation function of the $\bar{D}N$ system, we obtain a leading-order potential of the derivative expansion of the interaction kernel, which is then used to extract the $s$-wave phase shifts of low-energy $\bar{D}N$ scattering. Both the isospin $I=0$ and $I=1$ channels have a short-range repulsive core and a shallow attractive pocket in the intermediate to long-range region, though the $I=0$ channel is more attractive than the $I=1$ channel. We also observe that the $\bar{D}N$ potential exhibits more attraction than the $KN$ potential, which is its analog in the strange sector. In terms of the $s$-wave phase shifts, the $I=0$ channel shows a weak attractive behavior in the low-energy region with a positive scattering length of $0.246 \pm 0.105 (_{-0.051}^{+0.084})$ fm, whereas the $I=1$ channel shows repulsion with a negative scattering length of $-0.086 \pm 0.050 (_{-0.001}^{+0.037})$ fm. No bound states are found in both isospin channels, indicating the absence of a pentaquark state in the $s$-wave $\bar{D}N$ system.

The standard treatment of neutrino oscillations usually relies on factorization which assumes neutrino production, propagation, and detection are independent processes. As a consequence, the total probability is given by the product of production, oscillation and detection probabilities. As next-generation experiments are bringing neutrino physics to a high level of precision, the validity of this assumption must be checked. We present an S matrix treatment of the entire experimental chain, pion decay, neutrino propagation, and nucleon interaction, as a single, coherent quantum process. Our results reveal non-factorizable terms arising from spin and angular correlations between production and detection final states.In the $ΔL=0$ channel, these corrections introduce a $\sim 1\%$ systematic shift in the energy spectrum and a non-vanishing azimuthal asymmetry, important to be taken into account for precision measurements of $δ_{CP}$. For the $ΔL=2$ Majorana channel, we demonstrate that the S-matrix formalism is generating an azimuthal modulation that provides a direct way to access to the Majorana CP phases, which remain hidden in standard factorized effective mass approximations.

One of the major open puzzles in the Standard Model of particle physics is the strong CP problem: although Quantum Chromodynamics allows a CP-violating topological $θ$-term, experiments constrain its value to be extremely small. The Peccei--Quinn mechanism resolves this problem by promoting the $θ$-angle to a dynamical field-introducing the axion -- whose dynamics relax the effective angle $θ_\text{eff}$ to a CP-conserving minimum. Here, we investigate the resulting axion physics in a Hamiltonian lattice gauge theory (LGT) by coupling a quantized axion field to the massive Schwinger model with a topological $θ$-term. Using infinite matrix product state techniques, we compute the ground-state properties of the resulting theory and demonstrate that the axion dynamically relaxes $θ_\text{eff}$ to the minimum of the vacuum energy. Consequently, the ground-state energy becomes independent of $θ$, demonstrating the axion-mediated solution to the strong CP problem within a fully dynamical LGT. We further analyze CP restoration and extract the axion mass from the topological susceptibility and excitation spectrum. Our results provide a nonperturbative demonstration of axion dynamics in a quantum LGT amenable to investigation on modern quantum hardware.

We examine neutrino evolution in astrophysical environments where the neutrino flux is very large, including core-collapse supernovae and neutron star mergers. In these environments, the neutrino-neutrino and neutrino-antineutrino interactions are crucial. We include non-forward scattering of neutrinos and anti-neutrinos in a semi-classical treatment. Because of the large scale of neutrino momenta (2-10 MeV), the quantum evolution problem can be treated as a sum over incoherent paths in the and flavor of each neutrino. The phases between different neutrinos are essentially random because of the large kinetic terms. Momentum is conserved at each vertex, and important flavor symmetries are retained. Dynamics in the many-body neutrino system enable rapid equilibration in the energy and angular distributions of all flavors, and an equilibration of products of neutrino and anti-neutrino densities for each flavor at either large or zero background matter density. We also describe the evolution at moderate densities where the mass eigenstates differ for neutrinos and antineutrinos, and with time-varying background matter densities. The evolution maintains relevant symmetries and reduces to standard MSW oscillations in the appropriate limits. The rapid equilibration in energy and flavor can significantly impact energy deposition and nucleosynthesis in high-flux astrophysical environments, and potentially flavor energy relations in terrestrial supernovae neutrino observations.

FlexibleSUSY is a framework for an automated calculation of observables in user-defined models of a Beyond the Standard Model (BSM) physics, starting from the model's field content and its Lagrangian. Among a plethora of observables it is capable of calculating are the high precision predictions for Higgs bosons decay widths. Building on these previous developments we present here an interface between FlexibleSUSY and HiggsTools/Lilith. Combined with other FlexibleSUSY capabilities this extension provides a fully automatized tool chain leading directly from a user-defined BSM model to the state-of-the-art validation of the global agreement of a BSM Higgs sector with experimental measurements. We demonstrate this extension on a handful of phenomenologically relevant examples: a CP-conserving version of the Type-II Two Higgs Doublet Model, the CP-violating Next-to-Minimal Supersymmetric Standard Model and the Minimal R-symmetric Supersymmetric Standard Model. These examples show the power of FlexibleSUSY when applied to supersymmetric and non-supersymmetric models, both with and without CP-violation, and illustrate the handling of invisible and undetected decay widths.

Homotopy algebraic methods have become increasingly influential in studying field theories. We consider semi-holomorphic Chern-Simons theory and its relation with the principal chiral model. In particular, we establish an explicit quasi-isomorphism between the cyclic $L_\infty$-algebras governing both theories which directly gives the Lax connection. This provides a concrete example for studying integrability of a two-dimensional system through the homotopy algebraic lens.

We compute the disk 1-point function in timelike Liouville theory. Using the Coulomb gas formalism and analytically continuing in the number of screening operators, we derive an explicit formula, which is shown to satisfy the correct reflection symmetry, to have the expected self-dual properties, to fulfill the bootstrap shift-equations, and to reduce to previous known results in the appropriate limits. In the limit of zero cosmological constant, our result reproduces the one recently obtained in arXiv:2505.09390.

A search for long-lived particles in events with significant missing transverse momentum and at least one displaced vertex is presented. This analysis is performed using 137 $\text{fb}^{-1}$ of $pp$ collision data collected between 2016--2018 during Run 2 of the Large Hadron Collider by the ATLAS detector. Displaced vertices are identified using two different secondary vertexing algorithms, including a novel ``fuzzy'' vertexing algorithm optimized for identifying displaced decays of heavy quarks. Separate searches are performed using each algorithm, and the expected Standard Model background is independently estimated for each search using a data-driven procedure. No significant excess is observed over the background in either case. The results are used to set 95% confidence-level limits on potential beyond-the-Standard Model physics that could produce this final state. Results are interpreted in the context of four models: long-lived gluinos that form $R$-hadrons before decaying, neutralinos decaying via Higgs-mediated channels in the Bino-Wino coannihilation model, long-lived Higgsinos decaying to axinos, and an exotic Higgs portal model predicting displaced decays of light pseudoscalars.

We discuss that conventional Technicolour dynamics can be revitalized within the Dark Technicolour paradigm by invoking the Extended Most Attractive Channel hypothesis. In this framework, Standard Model fermions acquire masses via multifermion chiral condensates arising from new strong dynamics. The model incorporates three confining gauge sectors, Technicolour, Dark Technicolour, and an intermediate QCD-like sector, linked through extended gauge symmetries. The Extended Most Attractive Channel hypothesis reveals a hierarchical structure of condensates, where channels with higher net chirality become increasingly attractive. At low energies, the Dark Technicolour paradigm naturally reduces to the Froggatt-Nielsen or Standard Hierarchical Vacuum Expectation Value model, governed by residual discrete symmetries, offering a compelling resolution to the Standard Model Flavor Problem.

Axion-Like-Particles (ALPs) are hypothetical pseudo-scalar particles actively searched as light dark matter candidates. The coupling of ALPs to photons can give rise to distinctive spectral features in the observed gamma-ray spectrum of astrophysical sources. We perform a forecast study on the sensitivity to ALP-photon interactions using stacked mock observations of selected active galactic nuclei (AGNs) located behind galaxy clusters (GC). The ALP-photon conversion in the magnetic fields of galaxy clusters give rise to absorption-like features in AGN spectra that are subject to large variance in their prediction for individual sources. We consider here a stacking analysis of multiple AGN-cluster pairs, which yields a more controlled prediction of the expected ALP-induced spectral patterns in the observed gamma-ray spectra. Using realistic mock observations of selected Fermi-LAT AGNs by ongoing Imaging Atmospheric Cherenkov Telescopes such as H.E.S.S., MAGIC and VERITAS, we provide a careful assessment of the expected sensitivity of a combined statistical analysis of many AGN-GC pairs, together with the impact of modelling and instrumental uncertainties. The sensitivity reaches ALP-photon couplings down to 6$\times$10$^{-13}$ GeV$^{-1}$ for an ALP mass of 3$\times$10$^{-8}$ eV, and is currently statistically dominated indicating further improvements from more observations. Such a stacking analysis approach enables exploration of the yet-uncharted ALP dark matter parameter space in the 10$^{-8}$ - 10$^{-7}$ eV mass range.

We investigate the particle production in proton-antiproton ($p\bar{p}$) collisions using the PACIAE 4.0 model. The pseudorapidity density distributions ($dN_{\text{ch}}/dη$) and transverse momentum ($p_T$) spectra of charged particles from nonsingle diffractive (NSD) $p\bar{p}$ collisions agree well with the experimental data when using model parameters previously determined from nonsingle diffractive proton-proton ($pp$) collisions. Furthermore, we systematically compare results from both inelastic (INEL) and nonsingle diffractive $p\bar{p}$ and $pp$ collisions at the same energy to study the effect of the initial state (matter vs. antimatter) on the transverse momentum spectra of identified particles. Our results show that the net baryon-number difference in the initial state significantly enhances nucleon production at low collision energies, while its effect becomes negligible for high-multiplicity particles or at high collision energies, as expected. These findings further prove that the PACIAE 4.0 model is a versatile and reliable tool for studying high-energy collision physics.

We study torus knot invariants in the lens space $S^{3}/\mathbb{Z}_{p}$ within Chern--Simons theory. Using the surgery and modular description of lens spaces, we derive a general expression for the invariant of an $(α,β)$ torus knot in this background. In the large-$N$ limit these invariants simplify and acquire a universal form: the invariant of an $(α,β)$ torus knot in $S^{3}/\mathbb{Z}_{p}$ can be expressed in terms of the invariant of the $(α,α+pβ)$ torus knot in $S^{3}$. After an appropriate redefinition of knot variables, the generating functions of these invariants exhibit a structure analogous to quiver partition functions. Since the associated quiver is independent of the rank $N$ and level $k$ of Chern--Simons theory, the large-$N$ result provides a direct way to identify the underlying quiver, allowing us to determine the quiver structure associated with torus knots in $S^{3}/\mathbb{Z}_{p}$.

We revisit and extend the construction of bulk local states in flat holography, focusing on the induced representation obtained from the flat limit of the AdS highest-weight conditions. In three dimensions we clarify the scaling mismatch between bra and ket states in the flat basis and resolve it by introducing a dual basis, which yields a smooth flat limit and reproduces the correct Green's function. For higher dimensions we construct bulk local states explicitly, both in the momentum basis and in an alternative tilde basis. The flat limit of the AdS$_{d+1}$ construction is shown to be non-uniform in the descendant level and the Riemann-sum treatment over the scaling window $n\sim l$ converts the discrete descendant expansion into the continuum momentum representation, recovering the massive propagator. The tilde basis generalizes seamlessly to any dimension and is related to the three-dimensional flat basis by a sign factor. These results establish the induced representation as the correct algebraic foundation for bulk reconstruction in flat holography and provide a unified framework valid for arbitrary dimension.

We report the first measurement of the underground cosmic muon flux at the Stawell Underground Physics Laboratory. The measurement uses eight EJ200 plastic scintillator panels, equipped with Hamamatsu R13089 PMT pairs at the ends, which are the primary components of the muon veto system for the upcoming SABRE South experiment. The measured muon flux is f = (6.33 +/- 0.04_stat +/- 0.35_sys) x 10^{-8} [s^{-1} cm^{-2}]. This measurement is in excellent agreement with simulations, with a relative uncertainty an order of magnitude smaller than the modelling uncertainty.

We construct covariant theories incorporating fluctuating boundaries and soft cutoffs by introducing dynamical reference frames (DRFs). This framework generalizes the covariant action from a hard-cutoff to a soft-cutoff formulation, utilizing smearing functions and their corresponding operator expansions. This generalization initially leads to a loss of diffeomorphism covariance, which is recovered solely by restricting the DRFs, along with both their associated and linear MCFs, to specific forms, and by imposing suitable boundary conditions on the smearing functions. Satisfying these conditions restores covariance in relational spacetime, thereby enabling the consistent definition of subsystems. Within the covariant phase space formalism, we derive the charges of the soft-cutoff theory while explicitly addressing the inherent ambiguities arising from the boundary Lagrangian. We demonstrate that introducing an additional pointwise dependence is essential to resolve these ambiguities and ensure the integrability of the charges, even under fluctuating boundary conditions. Finally, in the context of General Relativity (GR), we establish the conditions under which holographic renormalization results at the asymptotic boundary coincide with the Noether charges derived from our soft-cutoff procedure.

The neutron to mirror neutron transitions in neutron stars would possibly result in significant effects. In this work we show that collisional decoherence entails exponential relaxation in lieu of oscillations. Decoherence is a great many orders of magnitude faster than the expected oscillations. The admixture of mirror neutrons at all times remains very small with respect to ordinary neutrons component.

We examine the proposal of [29] that M-theory may admit a Sugawara-type current algebra formulation based on $E_{11} \otimes_s l_1$. Motivated by the role of generalized coordinates in E-theory, we ask whether current algebra relations of this type can be derived in a setting that includes those coordinates systematically. We show that such a construction can indeed be carried out for a rigid $E_{11}$ model in which the generalized coordinates are treated as inert under the rigid symmetry, in contrast with E-theory. We also argue that the bilinear form entering the Schwinger term requires closer scrutiny, since any natural ad-invariant extension of the $E_{11}$ Cartan-Killing form to $E_{11} \otimes_s l_1$ is degenerate.

We investigate energy-momentum tensor (EMT) form factors and the spatial spin density distribution in the nucleon within a framework of the quantized Skyrme model with vector mesons. We construct both the canonical and Belinfante improved EMTs and analyze how pseudogauge uncertainty influences local spin and momentum densities while leaving the global nucleon properties unchanged. Using the inversion formulas from nucleon matrix elements in the forward limit, we extract the form factors, $A(t)$, $D(t)$, and $J(t)$, in both pseudogauges and the additional antisymmetric form factor associated with the canonical EMT. We find that the pseudogauge choice leads to sizable differences in the local spin and momentum densities. In particular, the canonical EMT naturally encodes spin density through the antisymmetric tensor structure, while the Belinfante EMT is sensitive to the total angular momentum only. Our results illustrate explicitly how different pseudogauges correspond to different spatial interpretations of nucleon spin structure within the same underlying dynamics. These findings provide a concrete model realization of the pseudogauge ambiguity in QCD-inspired nucleon structure and offer useful intuition for interpreting spatial distributions.

Fracton phases are new types of phases of matter characterized by subsystem global symmetry, which is a generalized global symmetry whose symmetry operator is partially topological. Their continuum low-energy effective descriptions admit two different formulations: an exotic quantum field theory (QFT) using exotic tensor gauge fields, and a foliated QFT constructed from a foliation structure and foliated gauge fields. For certain fracton QFTs, these two descriptions are equivalent, which is called the foliated-exotic duality. In this dissertation, we extend the foliated-exotic duality by combining it with the anomaly inflow mechanism for 't Hooft anomalies of subsystem symmetries. This dissertation has two main results. First, we discuss the exotic and foliated $BF$ theories in 2+1 dimensions, which exhibit the mixed 't Hooft anomaly of $\mathbb{Z}_N \times \mathbb{Z}_N$ subsystem symmetry. This anomaly is captured by a subsystem symmetry-protected topological (SSPT) phase for $\mathbb{Z}_N \times \mathbb{Z}_N$ subsystem symmetry in one dimension higher. By extending the foliated-exotic duality in the fractonic $BF$ theory to the SSPT phase, we establish the field correspondences in the SSPT phase and construct the foliated description of the SSPT phase. Second, we discuss the exotic $φ$-theory in 2+1 dimensions -- a fractonic gapless scalar field theory, which has the 't Hooft anomaly of $U(1) \times U(1)$ subsystem symmetry. The anomaly is captured by an SSPT phase for $U(1) \times U(1)$ subsystem symmetry in 3+1 dimensions via the anomaly inflow mechanism. Extending the foliated-exotic duality to the $φ$-theory, we establish field correspondences in the $φ$-theory and construct the foliated $φ$-theory that is equivalent to the exotic $φ$-theory. This provides the first example of the foliated-exotic duality in gapless theories.

By employing the gradient flow exact renormalization group (GFERG), we study the renormalization group (RG) flow of a manifestly gauge or BRST invariant non-perturbative ansatz of the 1PI Wilson action in quantum electrodynamics. The gauge invariance of the Wilson action is \emph{exactly\/} preserved under the RG flow. We explicitly solve the GFERG equation in the leading and partially the next to leading orders of the large $N_f$ approximation, where $N_f$ is the number of flavors. We obtain gauge invariant critical exponents and the gauge invariant 1PI Wilson action at an infrared (IR) fixed point for~$D<4$, where $D$ is the spacetime dimension.

We employ the Basis Light-Front Quantization (BLFQ) framework to compute the leading-twist (twist-2) and subleading-twist (twist-3) transverse-momentum-dependent parton distribution functions (TMDs) of the kaon. The light-front wave functions are obtained by diagonalizing a light-front QCD Hamiltonian that includes quark-antiquark (|q\bar{q}\rangle) and quark-antiquark-gluon (|q\bar{q}g\rangle) Fock components together with a three-dimensional confinement. Using the QCD equations of motion, the twist-3 TMDs are decomposed into twist-2 constributions and genuine twist-3 terms, the latter encoding quark-quark-gluon correlations beyond the probabilistic picture. These genuine twist-3 constributions arise from the interference between the |q\bar{q}\rangle and |q\bar{q}g\rangle sectors, which are usually negelected in the Wandzura-Wilczek approximation. This work provides the first theoretical predictions of kaon subleading-twist TMDs that explicitly account for Fock-sector interference. In addition, we present results for the kaon's twist-2 and twist-3 collinear parton distribution functions (PDFs). The twist-2 PDFs are found to be in good agreement with the recent global analysis by the JAM collaboration.

Every neutron star is born in the process of core-collapse supernova explosion that, for a brief moment, reproduces conditions of the early Universe with temperatures $T\sim O(30\rm\,MeV)$. We calculate the production of Dark Matter $χ$ from the SM particles in such events, SM $\toχ\barχ$, for the freeze-in range of couplings, $α_{\rm FI} \sim O(10^{-26}) $, finding that $O(10^{-6})$ $χ$'s per nucleon is produced. The strong gravitational potential well of the neutron star retains a substantial fraction of these particles that will eventually undergo the reverse process of energy injection, $χ\barχ\to$ SM. This may lead to the abnormal energy injection creating observable signatures such as late-time heating of the neutron stars. To demonstrate the power of this method, we construct a set of simple dark matter models coupled to lepton currents, and show that neutron stars provide unique constraints on parameter space that otherwise cannot be accessed by other means, probing effectively the scattering cross sections with the SM in the ballpark of $σ_{χ\,\rm SM} \propto O(10^{-70})\,\rm cm^2$.

We develop a systematic framework to quantify irreversibility in scalar Model A field theories with a generic non-Hermitian term driving the dynamics. Using the stochastic path-integral formalism, we perform a controlled small-noise expansion, allowing the computation of the entropy production rate (EPR) and violations of the fluctuation-dissipation theorem (FDT). We show that the local EPR is entirely determined by the anti-Hermitian part of the linearised Langevin equation. Around steady states, the non-Hermitian component produces linear corrections to FDT violations and contributes quadratically to the EPR. As an illustration of the applicability of our approach, we analyse a minimal non-Hermitian extension of the Ginzburg-Landau $ψ^4$ theory describing a non-reciprocal Ising model at coarse-grained scales, for which we obtain explicit expressions of the local EPR, showing that it localises at interfaces in non-uniform states. Our results provide a general characterisation of TRS breaking in non-Hermitian scalar field theories.

CP violation in $τ\to K_S πν_τ$ decays has attracted a lot of attention recently, due to the BaBar anomaly in the corresponding rate asymmetry. Within an effective field theory formalism, only extreme fine-tuning would allow to understand such measurement, which is currently being scrutinized at Belle(-II), as will be in the future super-charm-tau factory. Here we summarize the results of applying the same formalism to the other two-meson tau decay channels, which can help solve this conundrum. Our main conclusion is that current and future experiments would be sensitive to the maximum allowed CP rate asymmetry in the related $K^\pm K_S$ modes with a measurement having 5$\%$ precision, that will either support or cast further doubts on the BaBar anomaly.

We present Unified Flavor (UF), a framework that synthesizes the $B$-lattice flavor hierarchy with a dynamical realization based on TeV-scale vectorlike fermion (VLF) chains. Hierarchical Yukawa couplings arise from discrete ninths-quantized lattice exponents enforced by a single flavon $Φ$ with $ε\equiv\langleΦ\rangle/Λ=1/B$, $B=75/14$. Effective Yukawa entries are generated as algebraic path sums along nearest-neighbor chains of vectorlike quarks (VLQs), factorizing into entry, chain-propagation, and exit amplitudes controlled by the discrete gauge charges. A multi-messenger structure -- in which each Yukawa entry receives coherent contributions from several chain configurations -- generates O(1) complex coefficients whose phases are the physical origin of CP violation. We derive a general chain-inversion theorem, perform systematic perturbative diagonalization of both up- and down-type Yukawa textures, and show that the Cabibbo--Kobayashi--Maskawa (CKM) mixing hierarchy and CP-phase structure emerge naturally from the lattice exponent algebra and multi-messenger interference. All six quark masses are reproduced with O(1) coefficients that are essentially unity. The chain locality simultaneously suppresses dangerous flavor-changing neutral currents (FCNCs) and satisfies electroweak precision constraints, while requiring VLQs with masses in the multi-TeV range accessible at the High-Luminosity Large Hadron Collider (HL-LHC). The same discrete gauge symmetry that enforces the lattice structure also protects the Peccei--Quinn axion quality, unifying flavor, CP violation, and the strong CP problem. The framework extends to the lepton sector, reproducing charged-lepton mass hierarchies, the normal-ordered neutrino spectrum, and PMNS mixing with a predictive two-branch octant--$δ$ correlation testable at DUNE and Hyper-Kamiokande.

Our recent paper~\cite{FST} with Claudia Scheimbauer uses the cobordism hypothesis to construct fully local Chern-Simons theories. Here we expose some physics motivations: Yang-Mills plus Chern-Simons in the bosonic case and the free Majorana-Weyl spinor field in the fermionic case. We also give expositions of tangential structures and invertible field theories, in particular the 'gravitational Chern-Simons theory' used by Witten to obtain topological field theories from the physical theories.

In this article, we use two different methods for studying the mass spectra of fully-heavy baryons and pentaquarks. In the first section, we use state-of-the-art machine learning methods, such as deep neural networks and the Particle Transformer model architecture, to predict baryon masses directly from their quantum numbers, based on experimental information on hadrons from the Particle Data Group (PDG). We use this data-driven approach for the case of fully heavy baryons, and a large number of exotic pentaquark states, going much beyond the well-known $ P_c^+(4380) $ and $ P_c^+(4457) candidates. Subsequently,we extend the Gürsey-Radicati mass formula to incorporate the contributions of charm and bottom quarks, enabling analytical calculations for both ground and radially excited states of baryons and pentaquarks. The results obtained from both approaches demonstrate strong agreement with experimental data where available and make predictions for a number of unobserved states, including higher radial excitations. By addressing the question through both data-driven prediction and analytical modeling in different frameworks, this study offers complementary insights into the mass spectrum of conventional and exotic hadrons, guiding future experimental searches.

The particle nature of dark matter (DM) remains one of the central open problems in modern physics. Among the most extensively studied candidates are weakly interacting massive particles, whose parameter space is now under strong pressure from direct detection, indirect detection, and collider searches. In this work we revisit the Higgs-portal scenario with vector DM, first in an effective-field-theory description and then in a renormalizable UV-complete realization. We show that the effective Higgs-portal model with a Proca vector coupled quadratically to the Standard Model Higgs is essentially excluded over almost all of its parameter space by current direct-detection limits, with only a narrow region near the Higgs resonance surviving with a required fine tuning of the DM to Higgs mass that should at the permille level. We then consider a UV completion based on an additional gauged $U(1)_X$ symmetry, in which the DM candidate is a massive vector boson $V$ and the scalar sector is extended by a dark Higgs that mixes with the Standard Model Higgs. In this framework, the presence of a second scalar mediator opens an additional resonant annihilation channel and can substantially weaken the direct-detection constraints. In particular, when the DM mass lies sufficiently close to the heavy-scalar resonance, $m_V \simeq m_{H_2}/2$, the coupling required to reproduce the observed relic abundance can lie up to about two orders of magnitude below current direct-detection bounds, thereby opening viable parameter space that is absent in the effective description. Our results highlight the importance of going beyond the effective-field-theory approximation in Higgs-portal vector DM models and show that UV-complete realizations can qualitatively change the phenomenological conclusions.

The Simplified Template Cross Section (STXS) program has become the standard interface between Higgs measurements and global fits, but its fixed one-dimensional boundaries are not guaranteed to align with the phase-space directions to which the Standard Model Effective Field Theory (SMEFT) is most sensitive. We propose a machine-learning-inspired extension of STXS in which supervised classifiers are used only at the design stage to identify simple, publishable phase-space boundaries. Using associated Higgs production, $pp \to ZH$, as a case study and a benchmark momentum-dependent bosonic SMEFT deformation, we show that the relevant signal-background separation is well captured by a linear boundary in the $(p_T^Z,mZH)$ plane. We construct such boundaries with a linear support vector machine and with a deep-neural-network-assisted distillation procedure, and compare them directly with the standard STXS $p_T^Z$ bins through a common single-region Asimov-significance analysis. In this proof-of-concept setup, the ML-inspired regions systematically outperform the corresponding STXS regions, with the largest gains appearing in the boosted regime where SMEFT effects are concentrated. The final observable remains a simple linear cut, preserving the transparency and experimental portability that make STXS useful.

We investigate how perturbative reheating after inflation modifies the primordial gravitational wave (GW) spectrum generated by cosmological phase transitions. Within a specific inflationary setup, we show that the thermodynamic quantities that control the phase transition depend on the effective equation of state of the cosmological background, which is itself set by the form of the inflaton potential during reheating. Assuming reheating proceeds via perturbative dissipation of the inflaton condensate into boson or fermion pairs, we find that phase transitions taking place in this epoch generally produce GW signals that are systematically suppressed compared with the standard radiation-dominated scenario. We also identify characteristic spectral features that may arise in this case, which could serve as distinctive signatures of the modified expansion history during reheating.

In this paper we analyze the phases that dominate the microcanonical ensemble at various energies in the D1-D5 CFT, which is dual to type II string theory on $AdS_3 \times S^3\times T^4$. We focus on black hole solutions, and on the dependence of the phase structure on the ratio of the size of the torus to the AdS scale; as small localized black holes (with horizon topology $S^8$) grow, they can start to fill the $S^3$ or the $T^4$ or both, and we analyze the general aspects of the transitions between the various phases of uniform and non-uniform black holes, incorporating known solutions and discussing the properties of additional unknown solutions. Some features of the transitions between these phases are similar to higher dimensional AdS spaces, while other features are different. We provide evidence that when the torus is much larger than the AdS radius, there is a large range of energies where the typical states are a novel phase, described by a lattice (in the $T^4$ directions) of black holes with horizon topology $S^5\times S^3$. In this phase the entropy is linear in the energy, with a coefficient that is of order the AdS radius.

We propose a unified framework where the totality of dark matter (DM), the baryon asymmetry of the universe, and a detectable stochastic gravitational wave (GW) background originate from ultra-slow-roll (USR) inflation. The drastic suppression of the inflaton velocity during the USR phase, required for primordial black hole (PBH) DM production, can also set the initial conditions for spontaneous baryogenesis via a derivative coupling. This mechanism establishes a predictive correlation between the PBH abundance and the baryon yield, effectively fixing the reheating temperature $T_\textrm{reh}$ as a function of the post-peak spectral slope of the primordial power spectrum and the tensor-to-scalar ratio on CMB scales $r_\textrm{CMB}$. We perform a simple scan of the parameter space, demonstrating that while ``flat'' spectral tails allow for high-scale inflation ($r_{\rm CMB} \lesssim 10^{-3}$, $T_{\rm reh} \lesssim 10^{14} \text{ GeV}$) with a small wedge of tensor-to-scalar ratios potentially accessible to future CMB B-mode experiments, steep spectral tails enforce drastically lower scale inflation with an unobservably small $r_{\rm CMB}$ to avoid baryon overproduction. This degeneracy can be broken by GW astronomy: while LISA and DECIGO are capable of detecting the induced GW background associated with asteroid-mass PBH DM, the Einstein Telescope (ET) can act as a spectral discriminator, sensitive only to the broadband signals of high-scale scenarios.

Motivated by their role in non-perturbative potentials in string theory, we study divisors in effective cones of Calabi-Yau threefolds. We give examples of geometries for which some divisor classes in the effective cone are not themselves effective: i.e., they have no global sections. We call these non-holomorphic divisor classes "holes," and characterize their behavior in an ensemble of toric hypersurface Calabi-Yau threefolds. We prove some necessary and sufficient conditions for the existence of holes, show consequences of holes that follow from the minimal model program, and demonstrate that a class of holes come in semigroups (with this class conjectured to constitute all holes). Furthermore, we provide moduli-dependent bounds on the volumes of four-cycles representing holes.

We present a new self-supervised machine learning approach for symbolic simplification of complex mathematical expressions. Training data is generated by scrambling simple expressions and recording the inverse operations, creating oracle trajectories that provide both goal states and explicit paths to reach them. A permutation-equivariant, transformer-based policy network is then trained on this data step-wise to predict the oracle action given the input expression. We demonstrate this approach on two problems in high-energy physics: dilogarithm reduction and spinor-helicity scattering amplitude simplification. In both cases, our trained policy network achieves near perfect solve rates across a wide range of difficulty levels, substantially outperforming prior approaches based on reinforcement learning and end-to-end regression. When combined with contrastive grouping and beam search, our model achieves a 100\% full simplification rate on a representative selection of 5-point gluon tree-level amplitudes in Yang-Mills theory, including expressions with over 200 initial terms.

Baryon number-violating processes arise generically in many extensions of the Standard Model, with Grand Unified Theories providing the most compelling realizations. Ongoing experimental searches at JUNO, Hyper-K, and DUNE motivate a more precise and model-independent analysis utilizing the effective field theory (EFT) framework capable of connecting ultraviolet dynamics to low-energy hadronic observables. In this work, we perform a pipeline analysis that connects baryon number-violating (BNV) new physics to its low-energy description through the Standard Model EFT (SMEFT) up to dimension nine and the Low-Energy EFT (LEFT) up to dimension eight, and subsequently matches onto chiral perturbation theory using a systematic spurion method. We show that the complete set of dimension-eight LEFT operators enables the inclusion of the full set of chiral representations for three-quark operators with derivatives, which are also embedded in the chiral Lagrangian after quark-hadron matching. By contrast, the conventional dimension-six operators generate only the $(\bar{\mathbf{3}},\mathbf{3})$ and $(\mathbf{8},\mathbf{1})$ representations. The higher dimensional LEFT operators should be matched from higher order SMEFT operators, and thus a wider class of UV completions are allowed. We consider the complete tree-level UV resonances for dimension 6 and 7 SMEFT operators, and representative UVs for higher dimension operators.

We derive universal scaling relations for the low-energy effective action of string theory, connecting the vacuum energy and gauge couplings to higher-derivative Wilson coefficients. At one-loop in string perturbation theory, these generic parametric relations follow from modular and conformal invariance of the worldsheet, independently of the specific low-energy phase of the theory, and they become non-trivial in species limits. As a result, we substantially strengthen our previous case for the emergent string conjecture and connect UV/IR mixing to swampland principles. We argue that our results persist to higher loops, hinting at a pathway to study strong couplings using dualities. Further accounting for open-string contributions, if any, our results lead to parametric inequalities which reproduce holographic bounds and support the magnetic weak-gravity conjecture and the dark dimension scenario.

We analyze the isovector vector form factors of $D$, $D^*$, $B$, and $B^*$ mesons at low energies. We employ all constraints due to chiral and heavy-quark symmetry, and include the physics of resonant pion-pion rescattering in a model-independent way, using dispersion theory. Special attention is paid to the analytic properties of these form factors, which include anomalous thresholds due to triangle diagrams that are located on the physical Riemann sheets in some of the form factors. We extract the couplings of the $ρ(770)$ resonance to all these heavy mesons by determining the appropriate pole residues.

We investigate the asymptotic behavior of Kaluza-Klein (KK) towers in warped compactifications to Minkowski space. Focusing on the overall decompactification limit, we derive the scaling of KK masses at large KK momentum for scalar fluctuations in lower-dimensional Planck units. In codimension-one warped backgrounds sourced by a higher-dimensional exponential potential, we solve explicitly for the internal profiles and obtain a closed expression for the exponential mass decay rate $λ_{\rm KK}$ of the tower in terms of the moduli space distance. We find that warping reduces $λ_{\rm KK}$ relative to the unwarped case, in such a way that sufficiently strong warping could in principle violate the Sharpened Distance Conjecture bound. Remarkably, this sharpened bound is still satisfied precisely when the higher-dimensional potential obeys the condition forbidding asymptotic accelerated expansion, establishing a direct link between the Sharpened Distance Conjecture and the Strong de Sitter condition in one higher dimension. We also argue that for higher-codimension warped backgrounds the asymptotic KK scaling remains unmodified.

Spin and polarization are central to precision tests of fundamental physics and for interpreting radiation from astrophysical sources and ultraintense laser-matter experiments. Predictive modeling therefore requires not only energy spectra, but also angle-, spin-, and polarization-resolved particle distributions. Here, we demonstrate that a key assumption underlying current strong-field quantum electrodynamics (QED) models, i.e., that emission can be treated as an instantaneous random event sampled from a local differential rate, breaks down once emission angles, electron spin, and/or photon polarization are resolved. Namely, the resulting fully differential distribution can deviate strongly from the true result and can even yield inconsistent probabilities that take negative values. The physical reason is simple: a photon emission probability builds up over a finite length of the electron trajectory, the formation region, during which the electron direction changes by roughly the same small angle that defines the radiation cone. We therefore integrate over this formation region analytically to obtain a physically consistent electron spin and photon polarization model whose implementation is compatible with existing Monte Carlo and particle-in-cell (PIC) workflows. Simulations of a GeV-class electron-laser collision and of emission in a pulsar-like magnetic field reveal spin and polarization patterns that differ even qualitatively from state-of-the-art local models. In particular, our new model predicts substantial angle-dependent circular photon polarization where the standard approach yields none, and a pronounced helicity bias in the recoiling electrons absent from current predictions. These findings have direct implications for upcoming strong-field QED experiments and for interpreting polarized radiation from extreme astrophysical environments.

We provide a short review of the progress made in the past decade with functional QCD in the description of the phase structure of QCD. We summarise the most important technical aspects of the framework, discuss strategies for truncations and address the problem of systematic error estimates. We detail efforts to gauge the approach systematically with lattice QCD at zero chemical potential, also including the physics of the Columbia plot at non-physical quark masses. Our main focus is, however, the high density regime of QCD. We address the predictive power of the functional approach for the appearance of new phases beyond the chiral crossover regime for chemical potentials $μ_B/T\geq 4.5$. The onset of this regime may be signalled by a critical end point of the crossover line but may also involve a moat regime or the emergence of an instability that indicates an inhomogeneous phase. Respective results include estimates for the location of the onset of new phases, and predictions for their experimental signatures.

The suppression of heavy quarkonia in heavy-ion collisions is a well-established signature of Quark-Gluon Plasma (QGP) formation. However, recent LHC measurements of the $Υ(1S)$ state exhibit an anomalous high-$p_T$ plateau in the nuclear modification factor ($R_{AA}$) and a vanishing elliptic flow ($v_2$), challenging standard QCD transport models. We propose a viable mechanism to account for these observations by introducing a minimal dark scalar $φ$ situated within a strict kinematic merging window ($m_φ\approx 9.40$~GeV). We demonstrate that the shared $p_T^{-4}$ asymptotic fragmentation scaling between the hard-scattered dark scalar and Non-Relativistic QCD (NRQCD) color-octet production provides a constant theoretical dark fraction at high momenta. By extracting this fraction ($C_φ\approx 13.8\%$) from the anomalous $R_{AA}$ plateau, we establish a consistent phenomenological correlation: a single parameter addresses the $R_{AA}$ flattening, dilutes the inclusive $v_2$ toward zero, mitigates the long-standing quarkonium polarization puzzle, and naturally evades historical low-$p_T$ dimuon searches via a dynamic detector resolution threshold. We emphasize that future high-precision measurements of the dimuon mass lineshape at extreme transverse momenta are crucial for testing this paradigm.

We study the hadron-quark mixed phase within a two-phase model for symmetric and asymmetric matter. For the quark sector we employ the (2+1) Polyakov-extended Nambu-Jona-Lasinio model (PNJL) with vector interactions. We investigate how the hadronic equation of state affects the phase diagram and the thermodynamic properties inside the mixed phase. The behavior of isentropic trajectories in the mixed phase depends on the fixed entropy per baryon ($s/ρ_B$), with trajectories near the critical end point (CEP) exhibiting a pronounced cooling pattern, while isentropic trajectories with low entropy per baryon undergo pronounced heating as the baryonic density increases. The adiabatic squared speed of sound displays characteristic peak and dip structures that depend on $s/ρ_B$. The polytropic index along isentropic and isothermal trajectories, including in the vicinity of the CEP are also investigated. The effects of vector interactions and isospin asymmetry on thermodynamic observables likewise depend on the chosen $s/ρ_B$ value. Finally, we discuss the population of hyperons along isentropic trajectories and their influence on the phase diagram. The main effect of hyperons is to shift the onset of deconfinement to larger densities and decrease the density extension of the mixed phase.

Despite the success of the Standard Model (SM), several fundamental questions remain unanswered, such as the nature of dark matter (DM), motivating searches for new physics. This paper summarizes three recent searches for new physics in proton-proton collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV, using data recorded with the Compact Muon Solenoid (CMS) detector at the CERN Large Hadron Collider (LHC). The searches focus on "mono-X" final states, characterized by a large imbalance in transverse momentum recoiling against a single visible SM particle ($X$), and serve as powerful probes of new physics scenarios. Results are presented for searches in the pencil-jet (low-multiplicity jet), mono-photon, and mono-top final states, using CMS Run 2 data corresponding to an integrated luminosity of 138 fb$^{-1}$. No significant excess of events beyond SM predictions is observed, and the results are used to set stringent exclusion limits on various new physics scenarios, including simplified DM models and models of large extra spacetime dimensions.

The optimization and interpretation of microlensing surveys depends on having an accurate model of the Milky Way. However, existing Galactic modeling tools often perform poorly in replicating the stellar contents of the inner Galactic bulge region and reproducing microlensing survey results. We present an updated Galactic model implementation within the \synthpop framework that has been tuned for simulating the upcoming Nancy Grace Roman Space Telescope's Galactic Bulge Time Domain Survey (RGBTDS). We evaluate the model against stellar catalogs and kinematics from optical and infrared surveys toward the Galactic bulge, finding good agreement in much of the bulge, including the RGBTDS' contiguous lower bulge fields. However, within Galactic latitudes of $b\lesssim0.5^\circ$ of the Galactic plane, some inconsistencies arise which may impact projections for the RGBTDS' Galactic center field. In a following paper, we will examine RGBTDS simulations with this updated model in detail. Roman's GBTDS and Galactic Plane Survey will be instrumental in resolving the remaining model inconsistencies and improving our understanding of the structure of the central few degrees of our Galaxy.

In this paper we present a simple analysis around scaling relations derived from the Schmidt conjecture for star-forming molecular clouds, at the intra-cloud scale. Using a hierarchical tree (dendrograms) above a constant threshold ($A_V$ = 7 mag), we separate individual gas structures in a column density map of the nearby Giant Molecular Cloud Orion A, constructed from Herschel far-infrared maps. These structures define regions of dense molecular gas that can actively form stars. We also estimate their current embedded population using a list of known young stars. From the combined analysis of the column density map and the young star catalog, we construct a series of plots that show the intra-cloud level behavior of three well-known scaling relations: $N_{yso}$ vs. $M_{gas}$, $Σ_{SFR}$ vs. $Σ_{gas}$ and $R_{eq}$ vs. $M_{gas}$. Our dataset, along with other sets from literature, show the validity of a linear relation for $N_{yso}$ vs. $M_{gas}$, from intra-cloud to inter-cloud scales, over three orders of magnitude. We also especulate on the possibility that the relation could be valid over an even larger scale range. Additionally, our data are consistent with the $R_{eq}$ vs. $M_{gas}$ discussed in previous studies. However, our data is not quite in agreement with previously proposed fits for the $Σ_{SFR}$ vs. $Σ_{gas}$ relation, and we discuss the implications of using the free-fall timescale as the main parameter defining the star-forming efficiency in dense gas regions.

We present spectral analysis results of deeper (270 ks) NuSTAR observations of the merging galaxy cluster system, ZWCL1856.8+6616, at redshift z=0.304, following a pilot study using shallower (30 ks) NuSTAR data (Tumer et al. 2024). The cluster hosts a double radio relic, pointing to a similar mass head-on collision at/near the plane of sky. We aim to find the relation between radio and X-ray shock features. Using data from both focal plane modules of NuSTAR, we study the temperature structure across the field of view and report on the X-ray detected shock strength at the relic sites. We generate nominal and cross-ARFs with nucrossarf to disentangle photon cross-contamination within regions of interest due to the moderate point spread function of NuSTAR. Here we report one of the strongest X-ray detected shocks in a galaxy cluster merger with M=3.90(+1.64,-0.85) at the Northern relic site, that is unprecedentedly larger than the radio counterpart; M=2.5+/-0.2 (Jones et al. 2021a), and we report Southern shock strength as M=2.36(+0.58,-0.46). We argue that since the Northern relic (or radio shock), is confined in a very small region in the sky, particle acceleration is more efficient and is likely to grow in the post-shock regions. In addition, we search for inverse Compton (IC) emission at the radio relic sites; however, an IC component was not detected.

We present the first spatially resolved kinematic evidence for dust in the outflows of Active Galactic Nuclei (AGN). We utilise observations from JWST with NIRSpec IFU and MIRI MRS data of 10 local Seyferts and use Principal Component Analysis (PCA) tomography to extract the kinematics of Polycyclic Aromatic Hydrocarbon (PAH) features. PAHs comprise the smallest carbonaceous dust molecules in the Interstellar Medium (ISM), and produce emission features in the infrared providing the potential to measure kinematics. This is however challenging due to their broad shapes and variations in their intrinsic profile, prompting the need for techniques such as PCA tomography. We find that the velocity of the PAHs is similar to the molecular gas as traced by the rotational transitions of H$_2$, where for NGC 5728 and NGC 7582, both disk and outflow are present. We detect the outflow in the kinematics of large and neutral PAHs, namely the 11.3 $μ$m and 17 $μ$m PAH features, where after subtracting the disk, the velocity field matches that of high-ionisation potential lines such as [NeVI] (7.65 $μ$m, IP = 158 eV). Finally, we fail to detect kinematics of the 6.2 $μ$m PAH due to an altered intrinsic profile while the the 3.3 $μ$m PAH kinematics purely trace the circumnuclear disk. This suggests the PAHs in the outflow are more neutral and larger than in star-forming regions, consistent with PAH band ratios in previous studies of AGN.

Flare ribbons form when energy released by coronal magnetic reconnection is deposited in the low solar atmosphere, so by studying the dynamics of flare ribbons, one obtains an indirect measurement of reconnection. Our aim is to quantify the spatial and temporal scales of substructures in the Extreme Ultraviolet (EUV) flare ribbons, known as kernels, as a probe of the spatial extent and duration of energy injection during the impulsive phase of solar flares. Unprecedented observations of an M2.5 GOES-class flare from the March 2024 major flare campaign of Solar Orbiter were used. These data were obtained at high-cadence in short-exposure mode with the Extreme Ultraviolet Imager's high-resolution telescope, HRI_EUV. Individual kernels were automatically identified using a classical computer vision algorithm. Size distributions of ribbon kernels were derived, and an average light curve of individual kernels was extracted. The EUV flare kernels were small ($\lesssim 60~\text{pixels} \approx 1~\text{Mm}^2$) and a significant fraction were unresolved at a plate scale of 135 km/pix. Furthermore, we derived surprisingly short EUV kernel heating times of less than a few seconds. The average profile exhibits a sharp rise of $1.7\pm0.3$ s from half-maximum, requiring an additional $2.3^{+0.7}_{-0.4}$ s to return to its reference value. Our findings indicate that approximately half of the kernels were unresolved in this flare, despite the enhanced angular resolution offered by Solar Orbiter's proximity to the Sun at 0.38 AU here. Furthermore, we show that energy was only injected in a localised region ($\lesssim 1~\text{Mm}^2$) of flare ribbons for less than a few seconds. These results necessitate an in-depth investigation into the implications of such small-scale and transient injections on the energy flux deposited in solar flares, and the resulting response of the solar atmosphere.

We aim to study star-forming regions and the spectral energy distribution of two possibly interacting galaxies, PGC 56121 and PGC 56125, in the Hickson Compact Group 77. We utilized the far-ultraviolet (FUV) channel of the Ultra Violet Imaging Telescope (UVIT) on board AstroSat to observe and produce FUV images of the galaxies. Our FUV images are at a much higher resolution in comparison to those obtained from previous galaxy surveys by GALEX in the near-UV and those from PS-1, DSS. We have identified several star-forming regions in the two possibly interacting galaxies, PGC 56121 and PGC 56125. These two galaxies form a pair widely separated in redshift and are seen in projection. We also report the presence of a candidate tidal dwarf galaxy at the end of one of the tidal tails located to the east of the pair, and we identified three major star-forming regions in the tidal dwarf. The spectral energy distribution of three galaxies in the system is presented and analyzed to investigate the key physical properties, such as stellar mass, dust mass, total luminosity, and star formation history, of the three galaxies. Based on these observations and on a comparison with observations in radio, these three galaxies are probably part of a small group of interacting galaxies.

We run a suite of $N$-body simulations to investigate how classical bulges affect bar formation and properties under the internal formation mechanism. We incorporate bulges of varying mass and compactness into disk galaxy models and evolve them in isolation to examine the resulting bar pattern speeds and growth timescales. A more massive/compact bulge increases the Toomre $Q$ stability parameter and the circular velocity in the central region, while decreasing the disk mass fraction. It therefore delays the onset of bar formation and increases the bar growth timescale; sufficiently strong bulges can suppress bar formation entirely. During the formation stage, bars exhibit higher initial pattern speeds and faster deceleration rates when the bulges become more massive or compact. This faster deceleration persists after the bar buckling phase, leading to slower-rotating bars in the secular growth stage. However, when the bulge's "diluting" effect on the measured bar strength is removed or reduced, all bars within the same disk share similar distributions in the pattern speed-bar strength ($Ω_p$-$A_2$) space during the secular growth stage. They also show comparable ratios of the co-rotation radius to the bar length ($\mathcal{R}=R_{\mathrm{CR}}/R_{\mathrm {bar}}$) in this stage. These results suggest that the bulge's influence on the pattern speed is more significant during the bar formation stage, while in the secular growth stage, the bulge's effect may be less important, and the disk component dominates the pattern speed evolution.

The merger shocks generated by the collision of galaxy clusters elevate the pressure within the intracluster medium, significantly influencing the evolution of embedded cluster galaxies. We detect a merger shock (Mach number $\sim 1.3$) on the northwest side of the closest rich galaxy cluster Norma (A3627), using XMM-Newton and Chandra data. The textbook ram pressure stripping (RPS) galaxy ESO 137-001 appears to be located in the post-shock region. The shock boosts RPS and may induce the formation of the brightest known X-ray tail behind a cluster late-type galaxy. Another prominent head-tail radio galaxy ESO 137-007, with one of the longest radio continuum tails ($> 500$ kpc), is also likely in the post-shock region. The shock may have reversed the upstream jet to a one-sided radio head-tail morphology. Moreover, the shock can strip and roll the jet cocoon into a vortex ring structure like a `smoke ring' behind the end of the jet as observed by the ASKAP data. Therefore, the cluster merger shock can remarkably change cluster galaxies. Furthermore, Norma is the second brightest non-cool-core cluster following the Coma cluster, with a cool core remnant on its southeast side. Its original cool core may be disrupted by cluster mergers and/or active galactic nuclei.

The Milky Way bulge-bar is composed of multiple populations. Using chemical and kinematical planes, we segregate six populations in a bulge-bar sample observed by the APOGEE survey: two with bar-driven orbits, two with eccentric orbits, and two with low-eccentricity orbits, each composed of low- and high-[Mg/Fe] stars. Our sample spans $-2.0\lesssim$[Fe/H]$\lesssim+0.5$ and Galactocentric distance $R_{Gal}$ $<6$ kpc. We use chemical abundances from APOGEE DR17 for the elements Mg, Si, Ca, Al, K, Mn, Co, Ni, and Fe, and from the BAWLAS catalog for Ce and Nd. We find that the low- and high-[Mg/Fe] stars with low-eccentricity orbits, which exhibit chemical and orbital characteristics similar to those of the low- and high-[$α$/Fe] disks, display slightly negative and positive metallicity gradients, respectively. This result for the low-[Mg/Fe] low-eccentricity stars indicates a break in the global thin disk metallicity gradient. The high eccentricity populations with both low- and high-[Mg/Fe] show approximately flat metallicity gradients. In general, the [X/H] gradients of all elements for all populations follow Fe, except for the neutron-capture elements Ce and Nd. For all elements, the high-[Mg/Fe] bar population shows a much steeper positive [X/H] gradient than the nearly flat gradient for the low-[Mg/Fe] bar stars. The positive [X/H] gradients observed among our high-[Mg/Fe] bar stars probably reflect an age variation along the peanut structure. This interpretation agrees with the N-body simulations. Such steep positive gradients have also been reported in some high-redshift (z$\sim$4--10) galaxies.

We present a comprehensive multifrequency VLBI analysis of the FRII, high-excitation radio galaxy 3C 452, aiming to resolve and analyze for the first time its twin-jet structure on sub-parsec scales. Our data set comprises High Sensitivity Array (HSA) observations at 4.9, 8.4, 15.4, 23.6, and 43.2 GHz. Through fitting methods performed in both the visibility and the image plane, we trace the jet expansion from scales of a few thousand to nearly $10^5$ Schwarzschild radii ($R_S$) on both the approaching and receding jets. Additionally, we derive the core brightness temperatures and Doppler factors to constrain the jet's orientation and intrinsic speed. Our study provides the first detailed description of the twin-jet system in 3C 452 on VLBI scales, confirming it as a rare FRII source with jets detected down to millimeter wavelengths. We resolve both jet and counter-jet down to scales of a few thousand $R_S$, revealing a symmetric, parabolically expanding structure with power-law indices $k \approx 0.66$ (jet) and $k \approx 0.47$ (counter-jet). The brightness temperature analysis yields low Doppler factors ($δ\sim 0.03$-$0.83$), indicative of Doppler de-boosting due to the large viewing angle ($θ\approx 70^\circ$) and/or a magnetically dominated jet base. A spectral index analysis reveals a strongly inverted core spectrum ($α> 2$) with additional absorption at the highest frequencies, followed by a sharp steepening ($α\sim -2.5$) to optically thin values in the innermost jet. Finally, a comparison between broad- and narrow-line high-excitation radio galaxies shows that jets in narrow-line sources such as 3C 452 and Cygnus A complete collimation at $\leq 10^5 R_S$, whereas broad-line sources exhibit shape transitions at $10^6$-$10^7 R_S$, suggesting that orientation plays an important role in the observed collimation scales.

We describe the dynamical, photometric, and spectroscopic data available for stars targeted by Spitzer and Herschel to search for cold circumstellar dust emission from debris disks, a collection that we name the Cold Debris Disk Surveys (CDDS). These data include Hipparcos and Gaia parallaxes, 0.4-1250 micron photometry, spectral types, effective temperatures, gravities, bolometric luminosities, visual extinctions, metallicities, lithium abundances, rotational periods, projected rotational velocities, the Ca~II HK and IR triplet activity indicators, and X-ray luminosities for 3675 stars. Within this sample, we investigate the frequency of stellar and planetary companions (including potential new proper motion companions); use the data to assign CDDS stars to the field or one of many moving groups, open clusters, or stellar associations; and investigate correlations between stellar activity indicators. In future papers, we plan to explore the magnitude and frequency of infrared excess emission as a function of host star properties; to search for new companions with Gaia; and to examine the evolution of infrared excesses with the ages of stars in clusters and the field.

Magnetic fields play a crucial role in compact object physics, particularly in white dwarfs (WDs), where high densities can sustain strong magnetic fields. Observations have revealed magnetized WDs (MWDs) with surface fields reaching approximately $10^9\rm\,G$, although high-field MWDs are fewer in number in current catalogs owing to their intrinsic faintness and limitations in conventional electromagnetic surveys. In this study, we apply unsupervised machine learning (ML) techniques to systematically analyze a sample of hydrogen-atmosphere (DA) WDs. Using Uniform Manifold Approximation and Projection (UMAP) for dimensionality reduction and Density-Based Spatial Clustering of Applications with Noise (DBSCAN) for cluster identification, we classify distinct subpopulations within the DA WD sample. Each cluster exhibits unique intrinsic properties such as mass, surface gravity, temperature, and age. Our analysis further reveals that these subgroups effectively differentiate MWDs from non-magnetic or weakly magnetic counterparts. Moreover, utilizing a set of previously confirmed MWDs, we estimate the field strengths of all other MWDs lacking magnetic field measurements. This study underscores the effectiveness of ML-based approaches in astrophysical discovery, particularly detecting magnetized compact objects when direct measurements are unavailable.

Ultraluminous X-ray sources (ULXs) are non-nuclear X-ray binary systems that exceed the Eddington luminosity for a 10 Msun black hole. The majority of these sources are thought to be stellar-mass compact objects accreting at super-Eddington rates, exhibiting powerful relativistic winds. These winds have been identified through the detection of absorption lines with a blueshift as high as 0.3c and emission lines typically found at their laboratory wavelengths. In this work, we analysed the XMM-Newton data of the ULX NGC 5204 X-1, which has been observed to exhibit emission lines with a blueshift of about 0.3c. The aim of this study is to examine the geometry and physical properties of the accretion disc and the relativistic outflows. In addition, we aim to explore the factors that influence the ULX spectral transitions. We undertook an observing campaign with XMM-Newton to explore the source behaviour at different luminosities. In this first paper of the series, we performed high-resolution X-ray spectroscopy, including archival data, with the RGS instrument which allowed us to resolve both emission and absorption lines. The outflows features were characterised using physical models of plasma in collisional-ionisation and photoionisation equilibrium. We identify collisionally-ionised blueshifted and redshifted components at about 0.3c. These findings have high statistical significance and suggest a biconical structure for the outflow. Additionally, the analysis of the O VII line triplet observed in the spectrum enables us to infer physical properties of the low-velocity line-emitting plasma, e.g. electron density (ne $\sim 10^{10}$ cm$^{-3}$) and temperature (Te $ \geq 1.5 \times 10^5$ K). A hybrid plasma whose ionisation balance is affected by both collisions and radiation is favoured.

We present an analysis of archival Very Large Array (VLA) and Five-hundred-meter Aperture Spherical radio Telescope (FAST) \HI\ 21 cm data, together with archival multi-band radio continuum observations, of the nearby luminous infrared galaxy IRAS~04296+2923. The system, located behind the Taurus dark cloud at a distance of $\sim$29 Mpc, forms a small galaxy group consisting of five members as revealed by the \HI\ imaging. IRAS~04296+2923 has a close companion, HI~0432+2926, with a projected separation of $\sim$40 kpc, a small line-of-sight velocity difference of $Δ$ v = 26 km s$^{-1}$, and comparable total \HI\ masses of order $10^{9}$~$M_{\odot}$. Both galaxies exhibit regular \HI\ velocity fields and characteristic double-horn profiles in the VLA and FAST data, accompanied by only subtle asymmetries and extended \HI\ structures, indicating rotation-dominated kinematics with early signs of weak tidal interaction. Radio continuum emission is detected only from IRAS~04296+2923 and is confined to its nuclear region, consistent with previous studies. Modeling of its multi-band radio spectrum reveals a significant contribution from free--free emission at high frequencies ($>$30 GHz) and a high FIR-to-radio flux ratio ($q_{8.4}\simeq3.2$), implying a young, dust-obscured nuclear starburst. Taken together, the regular \HI\ kinematics, the small velocity offset, and the group-scale environment favor an interpretation in which IRAS~04296+2923 and HI~0432+2926 form a gravitationally bound, orbiting galaxy pair embedded in a small group, rather than an advanced merger. In this context, the luminous infrared galaxy (LIRG) nature of IRAS~04296+2923 is more plausibly driven by internal processes, such as bar-induced gas inflow, possibly modulated by long-timescale, low-level tidal interactions with nearby group companions.

Active galactic nuclei (AGNs) exhibit strong variability across the electromagnetic spectrum on a wide range of timescales, particularly in X-rays where fluctuations are both rapid and high amplitude. Power spectral density (PSD) analysis is commonly used to characterise this variability. Although AGN PSDs are typically well described by a bending power-law model, the dependence of the model parameters on photon energy has not been systematically investigated. We examine whether PSD parameters depend on energy using two highly variable, low-mass Seyfert galaxies, NGC4051 and NGC4395, as case studies. Using archival observations from XMM-Newton, Suzaku, and NuSTAR, we computed power spectra in six energy bands spanning 0.3-20 keV and fitted them with a bending power-law model to study the energy dependence of the PSD parameters. Power spectra derived from light curves obtained with different satellites and at different epochs are consistent within uncertainties, indicating that the X-ray variability process in both sources is stationary. For both AGNs we find that (i) the PSD bending frequency is consistent with being independent of energy, (ii) the high-frequency slope becomes flatter at higher energies, and (iii) the PSD amplitude decreases with increasing energy. These results place important constraints on models of AGN X-ray variability, such as the fluctuating accretion-rate model. Similar analyses of larger AGN samples are required to fully characterise the energy dependence of AGN power spectra.

Upcoming weak lensing (WL) surveys such as those by {\it Euclid}, LSST, and {\it Roman} require percent-level control over systematic effects. A common approach to mitigating baryonic effects uses semi-analytic baryon correction models (BCMs) that modify halo profiles in dark matter-only (DMO) simulations, calibrated to statistics from hydrodynamic simulations. We investigate the limits of this approach by progressively replacing larger regions around halos of decreasing mass in DMO simulations with their hydrodynamical counterparts. We compare multiple statistics -- the matter ($P(k)$) and weak-lensing ($C_\ell$) power spectra, peak counts, minima, one-point PDFs, and Minkowski functionals -- from "Replace" fields against hydrodynamical and DMO simulations. We find that replacing all halos with $M\geq10^{12}\,h^{-1}\,{\rm M}_\odot$ out to $r\leq5R_{200}$ recovers $\sim 90\%$ of the baryonic suppression in $P(k)$ and $C_\ell$ with the remaining $\sim 10\%$ originating from lower-mass halos or material farther outside of DM halos. Each statistic has distinct sensitivities to baryons: $P(k)$ and $C_\ell$ are sensitive to a broad range of masses and radii, whereas WL peaks are primarily affected by the cores of massive halos. We show that BCMs applied to massive halos and calibrated to match hydrodynamical $P(k)$ make two cancelling "mistakes": they underpredict core masses and compensate by overpredicting baryonic impacts at larger radii, thereby explaining previously reported failures of peak statistics in these models. We provide a framework for diagnosing critical mass/radius regions in baryonic modeling for a range of statistics for next-generation BCMs.

Globular clusters encode the hierarchical assembly history of the Milky Way and the physics of multiple stellar populations. Using homogeneous stellar parameters for 69 Galactic globular clusters derived while modelling multiple populations, we reconstruct progenitor-specific age--metallicity relations (AMRs) and test whether helium-related multiple-population (MP) properties depend on progenitor origin once cluster mass and metallicity are controlled for. Ages, helium spreads ($δY$), mean helium abundances ($\bar{Y}$), and first-population fractions ($f_{\rm P1}$) are drawn from hierarchical Bayesian CMD modelling. Progenitor families are identified via chemo-dynamical clustering, AMRs reconstructed within a hierarchical Bayesian framework, and MP indicators tested for environmental dependence. Enrichment timescales are consistent with $τ\lesssim 2$\,Gyr, though individual progenitors prefer shorter values when fitted independently. The primary distinction is the extent of chemical evolution: most systems reach $Δ[\mathrm{Fe/H}] \sim 1.1$--$1.3$\,dex while Sagittarius achieves ${\sim}1.6$\,dex and higher terminal metallicities. Gaia--Sausage--Enceladus and low-energy/Kraken are the dominant accretion events. Neither $δY$ nor $\bar{Y}$ depends on progenitor origin; the mass--MP scaling is indistinguishable across in-situ and accreted systems. Sequoia clusters alone show higher $f_{\rm P1}$ at fixed mass and metallicity. AMRs carry fossil signatures of progenitor chemical evolution and mass hierarchy. Helium enrichment amplitude is regulated by cluster mass and blind to environment, pointing to universal cluster-scale formation physics, with the sole exception of a residual dependence in $f_{\rm P1}$, suggesting the enriched-star fraction retains a secondary environmental imprint.

The binary fraction in young open clusters exceeds that of field stars, making the study of binary dynamical evolution in clusters essential for understanding the origins and evolution of field binaries. Using N-body simulations based on Gaia DR3 open cluster observations and assuming a 100\% primordial binary fraction, we investigated the early evolution of binary survival fractions in open clusters. We find that binary disruption has two stages, an initial rapid decline followed by a slower decrease, well described by two piecewise linear functions. The early disruption rate, $k_1$, follows a power-law relation with the cluster's initial density ($ ρ_\mathrm{0} $), with an index of approximately 0.56, driven by the disruption of wide binaries via close encounters. The transition time between the two phases, $t_\mathrm{b}$, also exhibits a power-law dependence on $ρ_\mathrm{0}$ with an index of about -0.46. The disruption rate also depends on binary parameters: high-$q$ and wide binaries are disrupted faster, while the dependence on eccentricity $e$ is less clear, likely due to its strong evolution. We developed and publicly released a Python tool to predict binary survival fraction evolution based on $ρ_0$, $P$ and $q$. Additionally, we also investigate how open cluster binaries contribute to the field population, and find that the escaped stars have a systematically lower binary fraction, likely due to mass segregation. Both populations show similar distributions of $ P $ and $e$, but lower-$q$ systems preferentially remain bound within clusters, the origin of which remains uncertain.

Variable continuum emission from AGN can be used to probe the structure of their accretion disks via reverberation mapping. Assuming a variable, hot inner light source irradiates the surrounding accretion disk, time delays between different continuum band light curves reveal light-travel times between their respective emission regions. Inter-band delays measured in several low-luminosity AGN are ubiquitously $\sim 3$ times longer than expected from standard disk theory, with evidence this size discrepancy may decrease in more luminous AGN. We have analysed high-cadence light curves of 9,498 of the brightest quasars between redshift 0.3-2.5 in the largest continuum reverberation study to date. Given the large sample size, we construct bins and fit delays jointly to combine inference across the parameter space and improve lag detections. We find that the size discrepancy persists in our high-luminosity sample, and that the previously seen anti-correlation with luminosity is likely driven by wavelength effects. The complex, non-monotonic wavelength dependence of delay amplitudes strongly suggests that contamination of inter-band delays by variable diffuse emission is widespread in the AGN population. We test delay behaviour against a variety of quasar properties finding longer lags in quasars with: higher Eddington ratios, redder colours, stronger optical FeII equivalent widths, higher iron ratios (both UV FeII/MgII and optical FeII/H$β$), CIV broad absorption troughs, and lower CIV blueshift.

Using the Gaia DR3 open cluster catalog, we identified the most massive star in each observed cluster. Examining the m_max$-$M_cluster relations across different age ranges, we find that as clusters age, the relation gradually deviates from the initial m_max$-$M_ecl relation and eventually exhibits clear age stratification. We conducted N$-$body simulations for both individual cluster evolution and subcluster coalescence. Four gas expulsion modes were tested for individual clusters, and two scenarios were modeled for cluster coalescence. Under all four gas expulsion modes, the evolution of the m_max$-$M_cluster relation follows a similar trajectory, differing mainly in evolutionary speed. The coalescence simulations show comparable behavior but align better with the observations, as both exhibit systematically lower m_max$-$M_cluster relations than individual cluster simulations. This systematically lower observed m_max$-$M_cluster relation suggests slower cluster mass loss and smaller masses for the most massive stars$-$both conditions reproduced in the coalescence simulations. Observations also show that clusters older than 5 Myr have most massive stars significantly deviating from the initial m_max$-$M_ecl relation. From this perspective, the coalescence simulations also provide a better match to the observations. In conclusion, the evolution of the m_max$-$M_ecl relation supports subcluster coalescence as a dominant pathway for open cluster formation, consistent with our previous work.

The origin of the mass scaling relation between supermassive black holes (SMBHs, $M_{\rm BH}$) and galaxies ($M_*$) remains a key open question. Rather than invoking AGN feedback, a non-causal mechanism has been proposed in which multiple mergers average out the $M_{\rm BH}/M_*$ ratio, thus decreasing its scatter ($σ$) and forming a tight local mass relation over cosmic history. A larger scatter in the relation at higher redshift suggested from a non-causal evolutionary scenario may be evident from recent JWST observations of overmassive SMBHs at high redshift. Here, we carry out a Monte Carlo simulation of solely merger-induced evolution of galaxies and their SMBHs which incorporates recent high-redshift observational constraints on $σ$ and the galaxy merger rate. We find that the dispersion in the local mass relation can be reproduced, even when starting from a highly scattered population at $z\sim6$ with $σ=0.8\,{\rm dex}$ or $1.0\,{\rm dex}$, which are in agreement with recent JWST studies. The redshift evolution of the scatter is highly sensitive to the mass ratio between merging pairs and the merger rate, and minor mergers with higher frequency than major mergers can also contribute to the scatter evolution, highlighting the importance of accurately constraining these parameters at high redshift through observations. Furthermore, statistical surveys aimed at determining the $M_*$-dependence of $σ$ and constraining $σ$ at $z\sim3-4$ will be effective in testing this scenario.

Bow shocks around cataclysmic variables (CVs) have traditionally been identified with a single bright optical arc. This feature has been interpreted as the bow shock formed by the interaction between a sustained outflow and the interstellar medium (ISM). We show that this interpretation is incomplete. Generic wind-ISM interaction theory predicts a two-shock configuration consisting of an inner terminal wind (reverse) shock and an outer forward shock, separated by a hot, low-density shocked wind cavity. Using archival ultraviolet, optical, and infrared imaging of the nova-like systems BZ Cam and V341 Ara, and the polar 1RXS J052832.5+283824, we find that the nebulae around all three systems exhibit this layered structure. In each case, the previously identified bow shock bright in Halpha and [OIII] corresponds to a compact inner arc, while additional emission components reveal a more extended morphology. Specifically, each system shows an outer arc detected in mid-infrared images, and the region between the optical and infrared arcs is filled with faint Halpha emission and, where available, far ultraviolet emission. We identify this infrared arc, reported here for the first time in these systems, as the sweep-up boundary of the forward shock, while the bright inner optical arc corresponds to the terminal wind shock rather than the forward shock as previously assumed. These results reveal that the true extent and layered structure of bow shocks around CVs only become apparent when observations extend beyond the optical band.

We revisit the physical origin of the achromatic peaks and breaks observed several thousand seconds after the burst in the multi-wavelength afterglows of XRF 080330 and GRB 080710. Using a numerical afterglow model that consistently incorporates finite ejecta thickness and a generalized external density profile, we perform Bayesian inference to estimate model parameters describing these events. Our analysis shows that the gradual rise and achromatic temporal features in both events are more naturally explained by jet dynamical evolution with finite shell thickness rather than by off-axis viewing effects. The inferred initial radial width of the ejecta is of order $10^{13}$ cm for both bursts, implying a central engine activity timescale significantly longer than that suggested by the prompt gamma-ray duration alone. Taken together, these results demonstrate that early afterglow light curves are strongly influenced by transition dynamics when finite ejecta thickness is properly taken into account, thereby providing a physical link between the prompt and afterglow phases and highlighting limitations of simply applying the thin-shell approximation when interpreting early-time afterglows. Furthermore, Bayesian model comparison favors a generalized circumburst density profile over the canonical uniform or steady-wind models, suggesting that fixing the external density structure to idealized profiles a priori may obscure crucial information about the progenitor's pre-burst activity.

We propose a unified scalar-field framework that addresses, within standard general relativity, three current cosmological anomalies: the $H_0$ tension, the mild preference for reduced late-time clustering ($S_8$), and recent indications of evolving dark energy. The model contains a single minimally coupled canonical scalar field evolving in a smooth potential composed of a localized bump superimposed on an exponential tail. The bump generates a transient pre-recombination energy injection that increases the expansion rate before last scattering, reduces the sound horizon, and shifts the CMB-inferred value of $H_0$ upward. After the field is released, its energy density rapidly redshifts through a kination-like phase, ensuring that the early modification does not persist as an unwanted late-time contribution. At low redshift, the exponential tail drives quintessence-like evolution, naturally yielding $w_0>-1$ and $w_a<0$ while suppressing linear structure growth and moving $S_8$ in the observationally preferred direction. The analysis shows explicitly how this smooth single-field potential can produce the required sequence of early enhancement, rapid dilution, and late-time thawing behavior.

pyKurucz is a pure Python reimplementation of Kurucz's SYNTHE, the standard code for computing synthetic stellar spectra. The original Fortran, written decades ago in a legacy dialect, is difficult to compile with modern toolchains without significant manual patching, and its long-term maintenance is uncertain following the passing of Robert L. Kurucz in 2025. pyKurucz is not a wrapper around Fortran but a line-by-line translation: continuous and line opacity from approximately 1.3 million atomic transitions, Voigt profiles, hydrogen Stark broadening, Saha-Boltzmann populations, molecular equilibrium for 189 species, and radiative transfer, all in Python with NumPy, SciPy, and Numba, requiring no Fortran at any stage. Validated against the original across 100 atmosphere models spanning 2500 K cool giants to 44,000 K O stars over 300-1800 nm at resolving power R = 300,000, it achieves sub-0.01% median agreement. The pure Python implementation enables direct integration with machine learning workflows and large-scale survey pipelines, while preserving an archival reference implementation of SYNTHE in a modern, readable language.

Fast radio bursts (FRBs) offer a powerful probe of the ionized Universe through their dispersion measures (DM). While a significant fraction of the DM arises from the intergalactic medium (IGM), the contributions from the host galaxy and the immediate environment of the source remain uncertain, and the physical origin of FRBs is still under active investigation. In this work, we investigated the possibility that FRBs originate from high-magnetic-field neutron stars (NS), whose magnetospheres can facilitate axion-photon mixing. Such mixing can modify photon propagation and induce an effective contribution to the observed dispersion. Using a sample of localized FRBs with measured redshifts, we perform a Bayesian Markov Chain Monte Carlo (MCMC) analysis to constrain the axion mass $m_a$ and axion-photon coupling $g_{aγγ}$. Within a parametric cosmological framework, we obtain $m_a = 1.16^{+4.40}_{-1.08}\,μ{\rm eV}$ and $g_{aγγ} = (1.76^{+6.69}_{-1.64})\times10^{-16}\,{\rm GeV}^{-1}$, together with a physically consistent intergalactic baryon fraction $f_{\rm IGM} = 0.837^{+0.053}_{-0.056}$. We further tested the robustness of our bounds against cosmological modeling assumptions by employing a non-parametric Gaussian Process reconstruction (GPR) of the DM-$z$ relation, which gives statistically consistent results.

One of the major challenges in astronomy involves accurately classifying galaxies, particularly distinguishing between different galaxy types. While many complex algorithms have shown strong performance in classification tasks, their complexity often results in longer processing times and increased difficulty in understanding. This study addresses this issue by exploring the viability of Fisher discriminants, a much simpler algorithm, in performing galaxy morphology classification. We tested four machine learning algorithms: the Fisher discriminant, Artificial Neural Networks (ANNs), Boosted Decision Trees (BDTs), and k-Nearest Neighbours (kNNs) to classify galaxies by the shape of their central bulges. Using data from the Sloan Digital Sky Survey (SDSS), we utilised five pre-processing transformations: normalisation, decorrelation, principal component analysis (PCA), uniformisation, and Gaussianisation, and classified the shape of central bulge into either rounded or no-bulge, based on the Galaxy Zoo Decision Tree. When compared to the Galaxy Zoo 2 (GZ2) labels, the Fisher discriminant with uniformisation obtained the highest accuracy score of 0.9310, outperforming ANN, BDT, and kNN by 1.93%, 0.42%, and 3.08%, respectively.

We present photometric and spectroscopic observations of Type IIP supernova SN 2023axu, spanning $\sim$400 d after the explosion. Its light curve is typical of normal SNe IIP, with a V-band peak of $-17.25 \pm 0.06$ mag and no early-time excess indicative of strong circumstellar interaction. The early spectra exhibit a distinctive broad "ledge" near 4600 Å. Through spectral modeling and comparison, we attribute this feature to a blend of C, N, and He lines excited by weak interaction between the ejecta and a low-density stellar wind. The late-time photometric evolution shows no discernible contribution from interaction, arguing against strong late-time circumstellar material engagement and supporting the low-density wind scenario. From modeling, this SN synthesized $\sim 0.055\,M_\odot$ of $^{56}$Ni, and nebular spectrum analysis indicates a progenitor mass near $15\,M_\odot$. SN 2023axu thus exemplifies weak ejecta-wind interaction and highlights the diversity of mass-loss histories and circumstellar environments of SNe II progenitors.

Ring galaxies are an uncommon class of galaxies whose morphology is closely related to dynamical processes that govern galaxy evolution. Some ring galaxies, known as "collisional ring galaxies", are thought to form as a consequence of head-on collisions between galaxies, and a number of high-redshift collisional ring galaxies have been discovered and/or studied in the era of the James Webb Space Telescope (JWST). In this paper, we present HST/ACS, JWST/NIRCam, and JWST/NIRSpec observations of a candidate ring galaxy at $z_{\rm spec} = 4.0148$, previously identified as a potential gravitational lens. The galaxy exhibits a complex morphology, including three bright clumps along an apparent ring with radius $\approx 0.25$" $\simeq 1.8$ kpc. It has a total SFR $= 140^{+20}_{-30}$ ${\rm M}_{\rm \odot}$ yr$^{-1}$ and $\log(M_\ast/{\rm M}_\odot) = 10.41^{+0.11}_{-0.13}$, making it similar to other high-redshift collisional ring galaxies. Although we argue strongly in favor of the collisional ring explanation, we cannot entirely rule out a galaxy-galaxy strong lensing explanation for the system's morphology, in which a foreground galaxy at $z \simeq 1.7$ lenses a galaxy at $z \simeq 4.0$ into an Einstein ring-like configuration; to confirm the nature of this source, we require kinematic information via high spectral resolution observations. We suggest that current and future gravitational lens surveys should consider high-redshift ring galaxies as possible but significant contaminants.

Low-frequency pulsar observations are crucial for understanding pulsar emission spectra and population physics, as well as for probing the interstellar medium (ISM) and Earth's ionosphere. We report the largest low-frequency pulsar census conducted in the southern hemisphere, covering 50-250 MHz, using the EDA2, an SKA-Low prototype station. In this survey, we detected 120 pulsars, including 23 first-time detections below 150 MHz and 5 below 100 MHz. For each source, we provide integrated pulse profiles and flux-density measurements across five sub-bands spanning 50-250 MHz. We also obtained improved dispersion measure (DM) values for 110 pulsars, with a median absolute DM correction of about 0.1 pc cm^-3. We measured significant Faraday rotation for 40 pulsars with improved rotation measure (RM) values for 4 pulsars, as well as phase-resolved RM variation in J1453-6413. Full-polarimetric pulse profiles are provided for all these pulsars, with multi-frequency polarimetric data for 20 of them. These results will enhance future SKA-Low science: refining pulsar population models, informing survey strategies, and advancing characterization of both the ISM and the ionosphere through low-frequency pulsar monitoring.

We perform three-dimensional radiation-relativistic magnetohydrodynamic (3D Rad-RMHD) simulations of accretion flows around spinning active galactic nuclei (AGNs). Our study focuses on the magnetically arrested disk (MAD) state, adopting a single-temperature model that includes bremsstrahlung opacity as the sole radiation process while varying the black hole spin from non-spinning to rapidly spinning cases. We find that the MAD state persists across all spin values, as demonstrated by the normalized magnetic flux at the horizon and the physically motivated spatially averaged plasma beta. The overall flow dynamics remain qualitatively similar for all spin models in 3D flow, suggesting that black hole spin has minimal influence on the accretion dynamics. In addition, we conduct post-processing using a two-temperature model to calculate the luminosities from synchrotron and bremsstrahlung radiation. We find that the total radiation luminosity is significantly higher than the luminosities from synchrotron and bremsstrahlung. This finding highlights the influence of radiation on the dynamics of the accretion flow. Our analysis shows that the electron temperature is significantly high in the jet region, regardless of spin. We further find that the temporal evolution of both radiative and synchrotron luminosities exhibits qualitatively similar behavior across all spin values. Finally, our results indicate that black hole spin has minimal impact on the spectral energy distribution (SED) in MAD state accretion flows.

Tidal Disruption Events (TDEs) release enormous amounts of energy, offering a promising avenue for detecting Population III (Pop III) stars. However, the radiative properties of TDEs of Pop III stars have so far been studied only analytically, relying on many assumptions. Based on our radiative hydrodynamic simulations that follow the evolution of the accretion system for Pop III star TDEs where a $300\ M_{\odot}$ ($M_{\odot}$ is the solar mass) star is disrupted by a $10^{6}\ M_{\odot}$ black hole (BH), we compute the emission properties of the event in rest frame and find that the spectrum peaks in the optical/UV waveband. After accounting for redshift ($z \sim 10$) and extinction effects, we find the observed spectral peak shifts to the infrared, with fluxes exceeding $10^{2}\mathrm{nJy}$-making such events detectable with both the James Webb Space Telescope (JWST) and the Nancy Grace Roman Space Telescope (Roman). The dependence of the observed spectrum on viewing angle is suppressed due to dust extinction. Using our simulation results, we also calculate the radio emission generated by the interaction between the wind and the circumnuclear medium (CNM) and find that a Pop III star TDE can produce an unusually long-lasting, continuously increasing radio flare with a duration greater than $10^4$ days and thus has the potential to be detected in radio wavebands. These results may be helpful to the detection of Pop III stars.

We investigated the timing and spectral properties of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary 4U 1608-52 using NICER observations. Our analysis reveals a correlation between the QPO frequency and its absolute amplitude, as well as between the frequency and the temperature of the burning layer. Intensity-resolved spectral analysis indicates that the flux modulation of the mHz QPOs is primarily caused by the variations in the blackbody temperature in most observations. Furthermore, for the first time, we report that as the source evolves from the soft spectral state toward the transitional state, the marginally stable burning responsible for the mHz QPOs ignites at deeper layers of the neutron-star surface. The radiation flux associated with the mHz QPOs shows a decreasing trend as the source moves into the transitional state. These two findings support a scenario in which the marginally stable nuclear burning ignites at deeper layers as the temperature decreases, releasing less energy from the nuclear reaction. Finally, we determine that the energy release rate of the marginally stable burning is around 10$^{35}$ erg/s, consistent with the theoretical predictions.

We report the detection of a clumpy, blue-dominated Ly$α$ emission at z = 3.066 located in the heart of a cosmic web filament in the MUSE eXtremely Deep Field (MXDF), spatially associated with the formation of two compact star-forming regions revealed by deep JWST/NIRCam imaging. Gas accretion in these regions is indicated by the blue-dominated Ly$α$ profiles, spectral signatures that are rarely observed. Radiative transfer simulation of the Ly$α$ profile using a clumpy multiphase model suggests a radial inflow of gas clumps with a velocity of 100 km/s. Embedded in this Ly$α$ structure, we find that the associated main galaxy dominates the stellar mass budget, while the two compact ultra-low-mass systems ($\log(M_\star/M_\odot) = 6.3\text{-}6.9$) have formed the bulk of their stellar mass in less than 7 Myr. These two components also have high specific star-formation rates, and elevated ionisation parameters, consistent with recent bursty star formation. This system provides direct observational evidence that how gas accretion, most likely from cosmic web, can induce starburst in ultra-low-mass galaxies.

The statistical properties of ionisation structures during the Epoch of Reionisation (EoR) provide valuable insights into the formation of the first stars and galaxies. However, statistics such as size distributions of ionisation structures can be affected by gravitational lensing caused by foreground massive structures like galaxy clusters. Hence, to quantify the impacts of lensing by galaxy clusters on ionised Bubble Size Distribution (BSD), we conducted a series of multiple-lens-plane lensing simulations involving the light cones of clusters alongside source light cones based on various ionisation models. The deflector population is generated using the Monte Carlo method, guided by halo mass function and empirical scaling relations, while deflectors' mass profile is modelled using the Truncated Navarro-Frenk-White (TNFW) model. Source light cones are produced via a semi-numerical approach or directly sourced from the Evolution of 21 cm Structure (EOS) project. By employing the Mean Free Path method, we measure unlensed and lensed BSD to reveal the lensing impacts. Our results indicate that lensing effects increase the number of large bubbles while leaving the number of small bubbles unchanged across all source models we adopted. Specifically, for the EOS faint galaxies model, the number of R > 15 cMpc bubbles increases by 219% at z = 14; for the EOS bright galaxies model, the above number increases by 832% under the same circumstances. Above all, lensing introduces unavoidable systematics for BSD, which must be carefully taken into account for relevant studies in the Square Kilometre Array (SKA) era.

We propose that an observed scarcity of remnant lobed AGNs in dense clusters results from a peculiarity in their dynamics upon the cessation of jet activity: a rapid `implosion' of lobes that, in their active phase, were primarily supported by the momentum flux of the jet. We investigate this behaviour by analysing the asymptotic behaviour of the RAiSE dynamical model and comparing our predictions both to the full model and hydrodynamic simulations. We find that remnant lobes powered by weak jets in massive clusters are unstable to implosion on the order of at most a few Myr. Consequently, remnant AGNs in massive clusters ($M_\text{halo} \sim 10^{14.5}$~M$_\odot$) will be under-counted by a factor of at least five compared to those in poorer groups ($M_\text{halo} \sim 10^{12}$~M$_\odot$). The lack of such remnants in observed populations may lead to a significant underestimate of the AGN feedback provided by low-powered jets, especially given their prevalence towards cluster cores where feedback is most effective. We discuss the influence of a stabilising magnetic field sheath on the nature of the implosion: does the lobe cleanly implode in on itself, or do fluid instabilities turbulently mix the lobe and ambient medium?

We measure the luminosity function (LF) of the Milky Way's stellar halo, using a magnitude complete, distance limited sample of stars from $Gaia$ DR3. Stars with high transverse velocities are selected, to isolate a high purity sample of the local halo. We adopt a cutoff transverse velocity of 250$\,$km$\,$s$^{-1}$, yielding 24,471 stars, and compute the halo LF, taking into account the effects of sample selection criteria. The LF displays similar features as are found in the well-probed LF of nearby, metal-rich disk stars, showing a strong peak at an absolute magnitude of around $M_G=10$, and a flattening near $M_G\sim7$ (Wielen dip). The $Gaia$ sample yields the first measurement of the LF continuously from the dimmest main sequence halo stars (subdwarfs) at an absolute $M_G$ magnitude near 13 mag to bright giants at $M_G\sim-3$. We obtain a local stellar halo number density of $1.7\times10^{-4}$ stars$\,$pc$^{-3}$ and disk-to-halo ratio by stellar number density of 480:1. We convert the $Gaia$ $G$ band measurements for our sample stars to Johnson-Kron-Cousins $V$ band, compute the $V$-band halo LF, and compare it to previous studies published over many decades that cover a wide range of techniques used. We discuss applications of the LF to the measurement of the luminosity and stellar mass of the Milky Way halo.

We aim to characterize nebular and stellar reddening in star-forming galaxies as a function of global galaxy properties (stellar mass, SFR, metallicity) at $2.7 < z< 7.0$. We also provide a prescription to convert SED-based $E(B-V)_{\mathrm{star}}$ to $E(B-V)_{\mathrm{gas}}$ when direct measurements of nebular reddening are unavailable. Our results are based on JWST/NIRSpec measurements of both individual spectra, with a sample of 283 galaxies, and composite spectra, including a larger sample of 327 galaxies. We estimate nebular reddening using the Balmer decrement (H$α$/H$β$) above $10^{8.5}$ $M_{\odot}$, where the sample is representative. Stellar reddening and SFRs are derived through Prospector SED fitting, while gas-phase metallicities are estimated from strong emission-line ratios. At fixed stellar mass, Balmer decrements remain consistent within uncertainties across our redshift range, indicating that stellar mass primarily determines the overall dust column even by $z \sim 7$. We find that differential reddening ($ΔE(B-V) \equiv E(B-V)_{\mathrm{gas}} - E(B-V)_{\mathrm{star}}$) scales linearly with mass and SFR at $z \sim 2.7 - 4.0$, but shows no mass or SFR dependence above $z \sim 4.0$. We find evidence for smaller $ΔE(B-V)$ above $z \sim 5.0$, suggesting that nebular emission and stellar continuum probe increasingly similar dust columns towards higher redshift. Finally, we find that nebular reddening correlates strongly with metallicity out to $z \sim 5$, whereas the correlation between stellar reddening and metallicity is weaker or absent. Together, these results suggest that both dust mass and geometry play a significant role in shaping the observed reddening of high-redshift galaxies.

The Near-UV-eXplorer (NUX) is a proposed ground-based, wide-field telescope array with a field of view of $\sim$70 square degrees, designed to operate over the 300-350 nm wavelength range and to achieve a target sensitivity of 20 mag in 150 seconds (5 sigma). Its main scientific objective is the detection and characterization of hot, rapidly evolving transients in the near-UV (NUV). Proto-NUX is a pathfinder instrument for NUX, based on an off-the-shelf 36 cm Celestron RASA wide-field astrograph that has been modified to enhance throughput and image quality in the targeted NUV band. The main objectives of Proto-NUX are: (1) to quantify the NUV sensitivity of the prototype and assess the feasibility of the full NUX facility; and (2) to characterize atmospheric extinction in the NUV, including its temporal variability and its dependence on zenith angle. Using three filter configurations, we aim to measure the wavelength dependence of the atmospheric extinction and to disentangle the contributions from Rayleigh scattering (dominating at wavelengths >325 nm) and molecular ozone-dominated absorption (dominating <315 nm). On-site testing is scheduled for 2026 at the Pic du Midi Observatory (France, 2877 m altitude) in order to evaluate on-sky performance under high-altitude observing conditions.

The X-ray emission of active galactic nuclei (AGN) is generally attributed to inverse Compton scattering of accretion-disk photons by hot electrons in a compact corona. In local AGN, directly constraining coronal properties is challenging because the high-energy cutoff often lies beyond the NuSTAR bandpass. High-redshift, luminous quasars enable systematic constraints on the high-energy cutoff, as cosmological redshift shifts the spectal cutoff into the observable hard X-ray band. We present first results from the ``Probing the AGN Coronae with High-redshift AGN'' (PACHA) project, based on quasi-simultaneous NuSTAR and XMM-Newton observations of 13 radio-quiet AGN at $z>1$. We constrain the high-energy cutoff and coronal temperature at 90\% confidence level for 10 and 9 sources, respectively. The sample exhibits a mean cutoff energy of $E_{\rm cut}=80.8\pm8.1$ keV and a mean coronal temperature of $kT_{\rm e}=18.4\pm1.6$ keV, both significantly lower than those measured in local {\it Swift}-BAT AGN, while the mean optical depth ($τ=4.8\pm0.3$) is significantly higher. The uncertainties are at 1~$σ$. Combining our high-redshift sample with local AGN, we find a potential anti-correlation between cutoff energy and both X-ray luminosity and black hole mass, with no significant dependence on Eddington ratio. Within a hybrid coronal framework, the inferred temperatures lie well below the pair-production limits for purely thermal coronae, indicating a substantial efficient Compton cooling and/or non-thermal electron component. The detection of low coronal temperatures in high-luminosity AGN is broadly consistent with predictions from recent radiation MHD simulations that consider purely thermal electron populations, implying that non-thermal electrons may not be the primary drivers of the observed coronal properties in these systems.

The two-point correlation function (2PCF) is a cornerstone of precision cosmology, yet its estimation from imaging surveys is vulnerable to contamination and incompleteness arising from imperfect target selection and pipeline-level inclusion decisions. In practice, the scientific target is a physically defined population, while the working catalog is constructed from noisy measurements and selection cuts, leading to mismatches between true and observed inclusion. These errors are often spatially structured, correlating with survey depth, observing conditions, and foregrounds, and can imprint spurious large-scale power or suppress the true clustering signal. High-resolution spectroscopic samples provide gold-standard inclusion in the target population but are typically available for only a small subset of objects. We introduce a prediction-powered Landy--Szalay (PP--LS) estimator that combines noisy inclusion labels across the full catalog with exact labels on a small spectroscopic subset while preserving the standard random-catalog normalization for survey geometry and selection. PP--LS debiases pair counts using residual-based, design-weighted corrections computed only on the labeled subset, requiring no probability calibration, known misclassification rates, or explicit modeling of contamination. Under simple random sampling of the labeled subset, we establish recovery of the oracle (true-label) Landy--Szalay pair counts and thus consistency for the target 2PCF. In simulations with clustered and spatially structured contaminants, PP--LS removes the bias of naive catalog-level estimators while achieving substantially lower variance than spectroscopic-only clustering. The resulting estimator is statistically principled, computationally lightweight, and integrates directly with standard pair-counting pipelines, enabling robust clustering inference in next-generation surveys.

We present a systematic survey of the neutrino spectral pinching parameter alpha_p(t, M, n-hat) across the Princeton Fornax ensemble of 3D core-collapse supernova simulations. We analyze 25 simulations spanning progenitor masses 8.1-100 M_sun with durations up to 8.47 s post-bounce, computed with the Fornax code and the SFHo equation of state. The pinching parameter alpha_p = (2^2 - E_rms^2)/(E_rms^2 - ^2) is derived from 12-bin spectral moments on a 128x256 sky grid for three neutrino species, enabling time- and angle-resolved spectral characterization. Four results emerge. (1) The nu-bar_e pinching floor is alpha_p = 1.92 +/- 0.10 (N=13 long-running models), lying 0.2-0.4 below 1D predictions due to 3D PNS convection. (2) Both BH-forming models (12.25, 14 M_sun) show anti-pinching (alpha_p < 0.9) before collapse, with deficit Delta alpha_p ~ 0.65 visible from t = 0.5 s. (3) Two of six long-running models exhibit a hierarchy reversal ( > ) after t = 5 s; leptonic flavors carry (40 +/- 3)% of radiated energy. (4) The LESA dipole is suppressed by >3x in BH-forming models; viewing-angle spread Delta alpha_p(68%) ~ 0.8-1.5 dominates spectral-inversion uncertainty. Mollweide sky maps reveal coherent angular structures with alpha_p anticorrelated with luminosity and correlated with mean energy. Detection rates at Hyper-Kamiokande, DUNE, JUNO, and IceCube yield 8-12% NMO/IMO discrimination during Kelvin-Helmholtz cooling. The late-time nu-bar_e pinching floor represents the first 3D characterization of spectral convergence during Kelvin-Helmholtz cooling, lying 0.2-0.4 below 1D predictions, a direct signature of 3D PNS convective transport

One-dimensional self-gravitating systems admit genuine thermodynamical equilibria. For systems with strictly monotonic orbital frequency profile, the Landau and Balescu-Lenard theories predict a relaxation time scaling linearly with the number of particles, $N$, in agreement with simulations. Yet, these theories become ill-posed for degenerate frequency profiles, as is the case in the harmonic potential, where all particles share the exact same mean orbital frequency. Using an exact collision-driven 1D integrator, we investigate numerically the self-consistent relaxation of 1D harmonic self-gravitating systems. We show that harmonic systems relax on a timescale that grows quadratically with $N$. We show that systems that are only partially degenerate display the same quadratic scaling for low $N$, but transition to the linear, non-degenerate behaviour for larger $N$. The larger the fraction of degenerate orbits, the larger the value of $N$ at which this transition of dynamical regime occurs. Finally, we explore the dynamics of fully non-degenerate systems, albeit with finite radial support: we confirm that their relaxation time scales linearly with $N$, though with a substantially larger prefactor than in non-compact systems. Astrophysically, this investigation should offer some new clues on the dynamics of density cores, as in the center of dwarf galaxies.

We present a study of the metal-poor (~6% Z_sun) massive (>8 M_sun) main-sequence star population in the star-forming dwarf galaxy Sextans A. By modeling near-UV to near-IR photometry of individual stars using the Bayesian Extinction and Stellar Tool (BEAST) we infer stellar parameters such as effective temperature, luminosity, and initial mass. We identify 867 massive main-sequence star candidates (present-day mass >8 M_sun and surface gravity >3.7 dex [cgs]) with a plausible spectral energy distribution (SED) fit, 500 of which show a probable SED fit. Comparisons to spectral types of existing observed spectra are consistent with the BEAST-derived stellar parameters, with most discrepancies explained. We identify 292 OBe star candidates through IR photometric signatures and find lower-limit OBe fractions of 15% for M > 8 M_sun, 23% for M > 15 M_sun, and 17% for M > 20 M_sun. We find 57 OB associations and that 24-28% of massive stars are isolated (distance to nearest massive star >28 pc). We discuss six likely runaway candidates (suggested velocities of ~ 50-340 km/s) not clearly associated with any major star-forming complexes. Lastly, we predict Lyman continuum (LyC) escape fractions of f_esc=0.27-0.76 across the star-forming regions and a global value of 0.35-0.71 by assuming low overall extinction and a range of porous geometries, indicating efficient leakage of ionizing photons. Future spectroscopic follow-up and resolved ISM studies will refine these constraints and solidify Sextans A as a benchmark for studying massive-star evolution and feedback at extremely low metallicity.

Alpha Crucis is the closest very high multiplicity massive star to the Sun. At its heart is the $4" \leftrightarrow 430 \text{ au}$ binary $α^1$ (A) + $α^2$ (B) Cru, which combined make up the 13th visually brightest star in the night sky. Here we make use of archival VLTI data of $α$ Cru A+B in order to study its multiplicity and orbital architecture. The data spatially resolved the close (6 mas) companion in $α$ Cru A (a known spectroscopic binary) and revealed that $α$ Cru B is also a close (17 mas) binary, which upgrades $α$ Cru to a seven star system. By combining the interferometric data with radial velocities, we solve for the full orbit of Aa+Ab and find dynamical masses $M_{Aa}=17.2\pm1.2 M_{\odot}$ and $M_{Ab}=6.8\pm0.3 M_{\odot}$. While the data on Alpha Cru B are not yet sufficient to tightly constrain all orbital parameters, we find that the orbital period is most likely 405 days (with 203 days also a possibility). The orientation of the orbital planes are sufficiently constrained to yield a mutual inclination between Aa+Ab and Ba+Bb of either $50 \pm 5^{\circ}$ or $137\pm5^{\circ}$, pointing to a dynamical formation scenario for the system. The photometric masses $M_{Ba}=12.4 M_{\odot}$ and $M_{Bb}=9.8 M_{\odot}$ together with the less massive wide component $α$ Cru Ca+Cb+D yield a total mass $M\simeq52 M_{\odot}$. At larger distances, the seven-star nature of Alpha Crucis would be arguably very challenging to unveil, suggesting that the companion frequency in massive star surveys may be underestimated.

Pa 30 has been identified as the nebular remnant of the historical SN 1181. It is host to a hot ($\approx200,000\,{\rm K}$) central star (WD J005311) with a fast wind ($\approx16,000\,{\rm km\,s^{-1}}$) radiating at roughly the Eddington luminosity for a solar mass ($\approx1.5\times10^{38}\,{\rm erg\,s^{-1}}$). We explore the thermal evolution of this star to understand how it progressed toward the state it is observed as today as well as to constrain its underlying physical properties. We develop a semi-analytic two-component model, which approximates the central star as a hot radiating envelope contracting and cooling above a relatively cool core. Comparing this model with the observed luminosity and radius requires a core mass $M_c\approx1.15-1.4\,M_\odot$ with a core radius $R_c\approx(6-8)\times10^8\,M_\odot$, and a hot envelope mass $ΔM\approx0.02-0.04\,M_\odot$. The small envelope mass is the best constrained of these parameters due to the need to reach the observed radius of $\approx0.15\,R_\odot$ in a timescale of $\approx845\,{\rm yrs}$. These results favor a picture where SN 1181 involved the merger of O/Ne and C/O white dwarfs, and where the majority of the latter was ejected in the explosion. We also explore which models ignite carbon burning at the base of the hot envelope, demonstrating that this is possible but not necessarily required to explain the current thermal state of the central star.

Modern X-ray and $γ$-ray instruments are revealing a growing class of Galactic non-thermal sources whose emission centroids are measurably offset from the nearest plausible sites of cosmic ray (CR) acceleration. Such "displaced" sources are seen in keV X-rays and TeV-PeV $γ$-rays but not in GeV $γ$-rays, have hard spectra, and are not associated with gas clumps, suggesting a leptonic origin. We develop a general framework for understanding displacement, whereby relativistic CR electrons (CRe) injected into the interstellar medium (ISM) with a strongly anisotropic pitch-angle distribution propagate a finite distance from their acceleration site before scattering processes isotropise their directions sufficiently for the emission to become visible. We use CR transport simulations to investigate under what circumstances displacement is likely, finding that it requires an initial pitch angle distribution $\lesssim 45^\circ$ wide, a line of sight broadly edge-on to the magnetic field, and that the source be measured in a waveband where emission is dominated by CRe for which the radiative-loss and pitch-angle scattering timescales are comparable. For typical Galactic conditions the latter condition is satisfied only for CRe energies $\gtrsim$ 10 TeV, explaining why displaced sources appear at X-ray and TeV but not GeV energies. We further show that, when displacement is detected, it allows a direct inference of the CRe pitch-angle scattering rate.

We investigate the link between optical obscuration and X-ray absorption in active galactic nuclei (AGN) by combining X-ray spectroscopy from 4XMM-DR11 with SDSS DR16Q spectroscopy. Bayesian X-ray spectral fits were obtained within the XMM2Athena project, and host-galaxy properties were derived via \textsc{CIGALE} SED fitting. Our final sample comprises 241 X-ray AGN at $z<1.9$. For 172 sources ($\sim70\%$), the optical broad-line (BL) or narrow-line (NL) classification agrees with their X-ray obscuration based on $N_{\rm H}$, but two mismatched populations emerge. Eleven BL AGN show signs of X-ray absorption (BLAbs) and elevated gas-to-dust ratios compared to BL AGN, consistent with dust-free or host-scale absorbers. Conversely, 58 NL AGN appear unobscured in X-rays (NLUnabs) and low gas-to-dust ratios. Nearly half are assigned type~1 properties by SED fitting, suggesting diluted or intrinsically weak broad-line regions, host contamination, or variability. Optical line diagnostics support this picture: NL AGN show higher Balmer decrements than NLUnabs, indicating stronger extinction and different ionization conditions. Host diagnostics further reinforce the contrasts: at $\rm z<0.8$, NLUnabs show star-formation rates and accretion efficiencies that are comparable to BL AGN, whereas NL AGN reside in more quiescent hosts with lower star formation and less efficient black-hole growth. BLAbs match BL AGN in host and accretion properties, with their peculiarity lying in excess X-ray absorption. These findings demonstrate that obscuration arises not only from orientation but also from multi-scale distributions of gas and dust. Identifying such mismatched populations will be crucial for interpreting AGN demographics in ongoing and upcoming surveys such as \emph{Euclid} and VRO/LSST.

This is the first paper in a series using our MACER framework to investigate the evolution of a disk galaxy, which emphasizes the role of active galactic nucleus (AGN) feedback and incorporates cosmological inflows. This paper presents the model setup and the overall results. The predicted AGN duty cycle of approximately 0.49% is consistent with observations. Analysis of the AGN luminosity and star formation rate (SFR) light curves reveals a positive correlation between the two. We find that cold filaments condense in the circumgalactic medium (CGM) region due to radiative cooling and subsequently fall onto the galaxy, significantly enhancing both the SFR and AGN activity. The galaxy is then quenched over a timescale of approximately 1 Gyr by the strong feedback from the enhanced AGN activity. This indicates that a positive correlation between SFR and AGN luminosity does not preclude AGN feedback from acting as the quenching mechanism for the galaxy. Notably, models without AGN feedback exhibit significantly lower peak SFRs than those with it. We attribute this difference to cumulative AGN feedback, which drives gas from the galaxy into the CGM, facilitating the formation of more massive cold filaments and ultimately promoting more intense starburst episodes.

Stellar streams are the result of a host galaxy's gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a globular cluster. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicites for over 10 million stars reaching significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$ km s$^{-1}$ and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel 'spur' feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a 'hot', metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Some gamma-ray bursts (GRBs) exhibit precursor emission episodes preceding the main emission, with a quiescent period in between. The properties of the precursor emission and the duration of the quiescent interval are related to the central engine activity and jet formation processes, thus providing insights into the physics of GRBs. We present a comprehensive analysis of the prompt emission and multi-wavelength afterglow of GRB 250916A. Using detailed afterglow modeling, we find that the broadband data are best described by a powerlaw structured jet with a relatively narrow core ($θ_c \approx 0.8^\circ$), viewed moderately off-axis at a viewing angle $θ_v \approx 2.7^\circ$. The isotropic-equivalent kinetic energy of the jet ($E_{k,iso} \approx 2.4 \times 10^{54}$ erg) is on the higher side for typical GRBs. The precursor emission is well described by a blackbody spectrum with a temperature of kT $\approx$ 13.2 keV and is separated from the main emission by a long quiescent interval of 150 s. Put together, our results indicate that the precursor is likely to be a shock breakout from a cocoon formed by the interaction of the relativistic jet with the progenitor star. The resulting cocoon pressure and shock collimation naturally lead to the launch of a narrowly collimated jet, consistent with the jet geometry inferred from afterglow observations. The long quiescent interval may imply the central engine turn-off in addition to the effect of the off-axis geometry.

Understanding the relative roles of large-scale environment and internal host-galaxy processes in shaping AGN activity is essential for constraining models of black-hole growth and galaxy evolution. We investigate how the environment of X-ray selected active galactic nuclei (AGN) relates to black-hole growth and accretion properties, and whether these introduce an environmental dependence beyond that expected from the host galaxy itself. Combining the XXL and Stripe 82X surveys, we construct samples of 427 broad-line AGN at $0.5<z<1.2$ and more than $20,000$ galaxies, with host-galaxy properties derived consistently using the same spectral energy distribution fitting methodology. Dark matter halo (DMH) masses are inferred from AGN--galaxy cross-correlation functions, while a multivariate nearest-neighbour matching algorithm is used to isolate trends with black-hole mass ($M_{\mathrm{BH}}$), Eddington ratio ($λ_{\mathrm{Edd}}$), and X-ray luminosity ($L_{\mathrm{X}}$) under controlled host-galaxy conditions. Within the uncertainties of the present dataset, X-ray AGN typically reside in halos of $\log(M_{\mathrm{DMH}}/h^{-1}M_\odot)\simeq13$, with no significant variation as a function of $M_{\mathrm{BH}}$, $λ_{\mathrm{Edd}}$, or $L_{\mathrm{X}}$. These results suggest that neither long-term black-hole growth nor short-term accretion variability is strongly linked to large-scale environment, and instead support a scenario in which AGN properties are regulated primarily by internal host-galaxy processes, while large-scale structure sets the broader boundary conditions for gas supply and duty cycle.

The observed GWB spectrum is higher in amplitude than model predictions by a factor of 2-3. Using a semi-analytic model, we evaluate the effect of a high-scatter supermassive black hole (SMBH) scaling relation (Mbh-Mbulge) on models of the nanoHertz gravitational wave background (GWB). By implementing an intrinsic scatter of the Mbh-Mbulge relation, which is larger at higher redshift, but matches local observations, we find that the amplitude of GWB models increases to be consistent with the low-frequency end of the GWB spectrum. This amplitude increase is not uniform across frequencies, a strongly evolving scatter preferentially increases the number density of the most massive SMBHs which, in the GWB spectrum, minimizes the strength of the low-frequency turnover. Our models with positively evolving intrinsic scatter can reproduce the electromagnetically observed overmassive SMBHs at 4 < z < 6 without changing the Mbh-Mbulge normalization though we find that including moderate normalization evolution improves fits to the GWB data. We conclude that the Mbh-Mbulge relation which best describes the available GWB and electromagnetic data sets has intrinsic scatter which evolves as epsilon(z) = epsilon_0 + (0.56 +/- 0.4) log10(1 + z) and normalization which evolves as alpha(z) = alpha_0 (1 + z)^(0.84 +/- 0.35). The results of this work imply that the Mbh-Mbulge relation we see today is not universal throughout cosmic time and that a diversity of seeding models and growth mechanisms may be at play in the early stages of SMBH-galaxy evolution.

We present FlowSN, a statistical framework using simulation-based inference with normalising flows to account for selection effects in observational astronomy. Failure to account for selection effects can lead to biased inference on global parameters. An example is Malmquist bias, where detection limits result in a sample skewed towards brighter objects. In Type Ia supernova (SN Ia) cosmology, these selection effects can systematically shift the inferred posterior distributions of cosmological parameters, necessitating the development of robust statistical frameworks to account for the biases. Simulation-based inference enables us to implicitly learn probability distributions that are analytically intractable to calculate. In this work, we introduce a novel approach that employs a normalising flow to learn the non-analytic selected SN likelihood for a given survey from forward simulations, independent of the assumed cosmological model. The resulting likelihood approximation is incorporated into a hierarchical Bayesian framework and posterior sampling is performed using Hamiltonian Monte Carlo to obtain constraints on cosmological parameters conditioned on the observed data. The modular learnt likelihood approximation can be reused without retraining to evaluate different cosmological models, providing a key advantage over other simulation-based inference approaches. We demonstrate the performance of this methodology by training and testing the simulation-based inference technique using realistic LSST-like SNANA simulations for the first time. Our FlowSN approach yields accurate posterior estimates on cosmological parameters, including the dark energy equation of state $w_0$, that are an order of magnitude less biased than those obtained with conventional techniques and also exhibit improved frequentist calibration.

Dwarf galaxies ($M_{*}/M_{\odot} \lesssim 10^{9.5}$), being sensitive to key evolutionary drivers like baryonic feedback and tidal perturbation, are crucial for understanding galaxy evolution as a whole. Their abundance and faintness, however, ensures that most will be studied primarily via broadband imaging for the foreseeable future. It is thus crucial to identify the most informative broadband-derivable quantities in the dwarf regime. As studies of widely used morphological parameters like concentration, asymmetry, and smoothness suggest these lack discriminatory power among dwarfs, we assess alternatives derived from isophotes: position angle twists, ellipticity, deviations from pure ellipses, and residuals to single-Sérsic profile fits. Using these parameters, we compare dwarf populations with massive galaxies of the same morphological class, and among themselves by morphological class. Only dwarf spirals may differ from their massive counterparts, being structurally simpler; dwarf and massive early type galaxy (ETG) isophotal similarity suggests all dwarf ETGs may be triaxial. Among only dwarfs, morphological classes are indistinguishable in this parameter space. A principal component analysis (PCA) using all available morphological, isophotal, and physical parameters expands on this: no PC explains more than $\sim$26% of the population variance, and no clear multimodality appears in any pairwise PC projection. We find similarly moderate spectral clustering, with a silhouette score of only 0.35. Given this self-similarity, parsing dwarf galaxy evolution from photometric parameters alone will likely require detailed statistical analysis of large dwarf populations in a high-dimensional parameter space, a task suitable for up-coming large-scale surveys like the Legacy Survey of Space and Time.

Perturbations from the Large Magellanic Cloud (LMC) of the Milky Way's stellar and dark matter haloes are well-established. However, studies have generally not considered haloes with high radial anisotropy, like debris from the Gaia Sausage-Enceladus (GSE) in the Milky Way. We run a series of test particle simulations of stellar haloes with different velocity anisotropies $β\in[0.5,0.9]$. The LMC causes these initially axisymmetric haloes to become approximately triaxial. Their major axes are aligned with its orbital plane and tilted by up to $\sim14^\circ$ with respect to a fixed Galactic disc. These effects become much more dramatic as $β$ increases, causing the halo to fractionate spatially according to anisotropy. This confirms the expectations of an analytical model, which predicts that orbits with eccentricities $e\gtrsim0.95$ should azimuthally align with the tidal field of the LMC. The reshaping of the $β=0.9$ halo creates strong overdensities of $\sim40\%$ at heliocentric distances as close as 15 kpc. These coincide with the well-known Virgo Overdensity (VOD) and Hercules-Aquila Cloud (HAC), which have previously been associated with the GSE. We propose that the HAC and VOD were created by the dynamical alignment of highly eccentric orbits by the LMC, and are not necessarily relics of the GSE merger geometry. We conclude that previous works have significantly underestimated perturbations from the LMC in the inner stellar halo by not considering sufficiently high velocity anisotropy. This effect should be corrected for when constructing equilibrium models of the GSE.

Solar twins are among the most powerful tracers of Galactic disk evolution owing to their unique property of sharing nearly solar metallicities ([Fe/H] ~ 0) while spanning a wide range of ages. To grasp solar twins as relics of Galaxy evolution, individual twins must be tagged with ages. A sufficiently large and well-characterized stellar sample then allows us to construct an age distribution that encodes the star formation history beyond our local region, modulated by the efficiency of radial migration of stars. Based on our catalog of 6,594 high-quality local (<~ 300 pc) solar twins from the Gaia Data Release 3 spectroscopic (GSP-Spec) catalog, we derived their age distribution after carefully deconvolving the selection function. We find two distinct features: a narrow peak around ~ 2 Gyr and a broad bump extending over ~ 4--6 Gyr. First, we argue that the former corresponds to a relatively recent burst of star formation that occurred in the disk, including at least a local region within a few kiloparsecs of the Sun, which is in good agreement with previous results deduced from independent works. On the other hand, the older bump, closely associated with the Sun's birth epoch, is intriguing since this finding challenges the predicted presence of a corotation barrier built by the Galactic bar, which is thought to prevent stars born inside R_GC ~ 6 kpc from reaching the solar neighborhood. We propose that the large number of local twins with ages between 4 and 6 Gyr provides compelling evidence that the Sun's long-distance (>= 3 kpc) migration is shared by many inner disk stars. This, in turn, suggests a possible link with the epoch of bar formation, which may have triggered enhanced star formation in the inner disk and induced efficient radial migration.

Measurements of the kinematic cosmic dipole continue to show an intriguing tension between the value inferred from the CMB and that obtained from high-redshift source number counts. While the measured dipole direction appears consistent, the amplitude, set by the observer's peculiar velocity $v_{o}$, remains in significant disagreement. In this paper, we propose using strong gravitational lenses with well-measured Einstein radii to estimate the kinematic cosmic dipole, through the relativistic distortion of the Einstein angle induced by the observer's motion. We show that this effect could be detected solely from measurements of the Einstein radius in wide, high-resolution imaging surveys such as Euclid. However, the precision achievable using Einstein-radius measurements alone, without redshift or lens-galaxy mass information, appears insufficient to discriminate between the CMB value of $v_{o}$ and that derived from source number counts at high statistical significance. Nevertheless, we demonstrate that including a large sample of lenses with available kinematic information, either via the Fundamental Plane relation or, ideally, through spectroscopic velocity-dispersion measurements, drastically reduces the noise and substantially improves the constraining power of this method. We show that, for a realistic sample of strong lenses detected by Euclid and complemented with spectroscopic velocity dispersion measurements from 4MOST or DESI, it is possible to discriminate between the CMB- and source-number-counts-inferred values at the $\sim 4σ$ level using a new, fully independent method. We further demonstrate that this technique is only weakly sensitive to strong-lensing selection effects, with selection biases and threshold effects estimated to be well below the 1% level.

Observations with JWST in the F090W band provide a powerful tool for determining galaxy distances based on tip of the red giant branch (TRGB) measurements. It is a great convenience that the TRGB lies at an almost constant absolute magnitude level at low metallicities. However, the TRGB becomes fainter at high metallicities in the F090W filter. Details of this break in slope are critical for precision applications in the acquisition of distances. With an absolute scaling set by the maser distance to NGC 4258 (but excluding the uncertainty in that distance), the value $M^\mathrm{TRGB}_\mathrm{F090W} = -4.40 \pm 0.03$ mag (traditional Vega) is found for $(\mathrm{F090W}-\mathrm{F150W})_0<1.65$ mag. The theoretical RGB isochrone that reaches the color 1.65 at the RGB tip corresponds to metallicity $[M/H] = -0.57$ for a 10 Gyr population. The calibration is used to derive distances for 16 galaxies relative to the megamaser host NGC 4258. Revised distances are on average slightly closer than literature values derived from the same data.

Ultraluminous X-ray sources represent extreme super-Eddington accretion regimes, and a subset is now known to host highly magnetized neutron stars. However, direct observational probes of their surface magnetic fields remain scarce. In this Letter, we report the detection of a narrow X-ray absorption feature at $3.29\pm0.02$ keV in the XMM$-$Newton/EPIC-pn spectrum of NGC 4656 ULX-1. The source exhibits a hard-ultraluminous state, while our timing analysis reveals a candidate pulsation at $\sim$0.9736 Hz, with a local significance of $5.5σ$ and a pulsed fraction of $\sim11\%$. The feature is robust against changes in continuum modeling and data-selection criteria, retaining a statistical significance of $\gtrsim3σ$ in Monte Carlo simulations. Interpreting the absorption as a proton cyclotron resonant scattering feature implies a local magnetic field of $B\sim(6-7)\times10^{14}$ G in the line-forming region. This value is consistent with strong magnetic fields anchored near the neutron star surface, even if the large-scale dipole is substantially weaker. Although we discuss electron cyclotron features and atomic transitions as possible alternatives, they appear less consistent with the observed phenomenology.