Papers

The fixed-ray escort integral representation expresses $\alpha$–$z$ Rényi divergence as an average over ordinary relative entropy of a family of escort states. Working in the standard UV-regulated density-matrix description of QFT subregions, we use this representation to derive an escort-averaged entanglement first law, an escort-averaged representation of the $\alpha$–$z$ information kernel, and an escort-averaged Bekenstein-type bound for ball-shaped regions in conformal field theories. For the conjectural $\alpha$–$z$ quantum null energy condition (QNEC), we obtain a quasi-local form in which the null energy is evaluated in an escort-averaged state and is corrected by an escort-transport term encoding the failure of escort formation to commute with restriction to a null-deformed region. The $z=\alpha$ specialization gives a similar quasi-local form of the Rényi QNEC for sandwiched Rényi divergence. We explicitly compute the Rényi QNEC, including the explicit escort transport term, for coherent-state excitations in a free scalar field theory. For the same coherent family, we obtain a positive $\alpha$–$z$ null Hessian, verifying the conjectured diagonal $\alpha$–$z$ QNEC for this family.

Rapid advances in agentic artificial intelligence enable scientific systems to interpret open-ended objectives, combine heterogeneous information, invoke specialized tools, and revise experimental strategies as evidence accumulates. However, physical experimentation also contains many tasks for which agentic reasoning provides little advantage and can reduce reliability. Quantitative analysis, optimization, spatial targeting, validation, and instrument execution are often better posed as deterministic or algorithmic operations with explicit objectives and verifiable outputs. Here, we introduce a hierarchical architecture for autonomous experimentation that separates these roles. Agentic components interpret scientific intent, construct task-dependent experimental representations, evaluate accumulated evidence, and select high-level actions, whereas deterministic algorithms perform numerical analysis, coordinate selection, validation, and physical execution. We implement this architecture in piezoresponse force microscopy. Starting from a broad scientific question concerning the relation between local domain structure and polarization switching, the system constructs spatial descriptors from multichannel imaging, selects and analyzes local hysteresis measurements, adapts the spectroscopy waveform, and terminates the experiment when additional measurements cease to provide new evidence. The autonomous trajectory also identifies a confounding relationship between polarization state and domain-wall proximity and recognizes that the requested contrast is not independently represented within the available field of view. These results demonstrate a route toward scientific autonomy in which agents determine what evidence is required while algorithms determine how that evidence is acquired reproducibly and within validated physical constraints.

The Kalnajs matrix method is a widely used framework for studying the global linear stability and possible Landau damping of collisionless stellar systems. However, for radial perturbations ($l=0$) in three-dimensional spherical models with infinite boundaries, standard biorthogonal sets such as the Clutton-Brock basis often converge slowly. This stems from a physical constraint: mass-conserving radial modes force the perturbed potential to decay faster than a point mass at infinity, whereas individual Clutton-Brock elements carry a fictitious net mass and decay only as $\mathcal{O}(1/r)$, producing unphysical asymptotic tails. We construct two new families of potential-density basis pairs that are free of these tails by construction. The first modifies the Clutton-Brock set: a specific linear combination of adjacent elements analytically cancels the leading $\mathcal{O}(1/r)$ term, giving potentials that decay as $\mathcal{O}(1/r^3)$ and densities as $\mathcal{O}(1/r^5)$. This breaks strict biorthogonality but yields a compact tridiagonal Gram matrix that enters the response equation at negligible additional cost. The second family is built from Jacobi polynomials that embed the required $\mathcal{O}(1/r^2)$ potential decay directly into their construction while retaining strict diagonal biorthogonality. Numerical tests demonstrate convergence that is uniform in radius for both expansions: the truncation error decreases exponentially when the asymptotic tail of the expanded potential matches the parity of the basis elements, and algebraically otherwise. Both bases reduce the dimension of the response matrix required for a given accuracy and are suited to studying radial perturbations in open stellar systems.

Modeling complex physical phenomena directly from experimental data poses fundamental challenges: measurement devices introduce noise and biases into the data, the governing equations are often unknown or intractable, and the dynamics observed experimentally may exhibit stochastic behavior. In this paper, we consider capillary wave turbulence—an example of nonlinear wave interactions at a microfluidic interface—measured by ultra-high-speed digital holographic microscopy. With the goal to learn an efficient model directly from these data, we adapt a stochastic extension of Operator Inference to learn low-dimensional stochastic differential equation representations of the microscale wave dynamics from experimental measurements. Each experiment is repeated 16 times, and 10 different conditions are created by changing the nondimensional acoustic capillary number through an excitation device. In addition, we propose a new strategy to select the reduced-order model dimension, based on both mean and covariance errors. We further add Tikhonov regularization to the stochastic Operator Inference framework and show that, beyond its conventional role as a numerical stabilization technique in deterministic settings, it has a physically meaningful interpretation as controlling the spectral content of the learned stochastic dynamics. We demonstrate that the resulting stochastic reduced-order models faithfully capture salient physical features of capillary wave turbulence across a range of experimental conditions.

The number of metastable states of a random corrugated landscape controls what dynamics, inference and learning can reach. In the Sherrington-Kirkpatrick model these states are marginal, and counting them requires breaking the BRST symmetry. So far only the annealed complexity was known, predicting states below equilibrium. We compute the quenched complexity function with a two-group Ansatz and one-step replica symmetry breaking. It branches continuously off the annealed curve at its instability point and vanishes near the Parisi free energy, while the soft mode turns toward orthogonality to the magnetization.

The self-assembly of colloidal particles enables the creation of structured materials with programmable functionalities; however, controlling interaction specificity and aggregate morphology in a reversible and scalable manner remains a major challenge. Here, we investigate the selective depletion-induced self-assembly of 3D-printed flat polygonal colloids, where nanoscale surface topography is engineered through precise modeling in two-photon polymerization. By designing anisotropic lateral surfaces, we direct specific interactions that govern aggregate morphology, yielding dimers, chains, zigzag, and honeycomb structures depending on the surface configuration. The specificity of interaction is tuned by varying the length scale of the topographic surfaces, the depletant concentration and the ionic strength of the solution, revealing a transition from selective to non-selective aggregation regimes. The relative placement of lateral interacting surfaces on the colloids enables assembly into aggregates spanning a broad range of sizes, while tuning the interaction strength selectively stabilizes distinct structural motifs. We demonstrate this interplay between geometric arrangement and interaction energy experimentally and corroborate through both theory and simulations for specifically hexagonal shaped colloids. This study establishes a versatile framework for programming colloidal interactions via micro-architectural design, offering new routes for fabricating reconfigurable and functional soft materials.

We present measurements of the hyperfine structure, transition frequencies and relative intensities in 14 lines of 169Tm+ ions starting from the $4f^{13}({}^2\!F^o_{7/2})6s_{1/2} (7/2, 1/2)_{3,4}^o$ levels. The magnetic dipole hyperfine ($A$) constant of all 15 involved ionic levels are determined, 14 of which have sub-MHz accuracy. For five levels no previous measurements were reported, while all others are in excellent agreement with the literature and improve on the previous accuracy by typically 1 - 2 orders of magnitude. We compare transition strengths and include improved transition frequencies. These results are useful for the quantum computing community, modelling of Tm$^+$ abundances in stellar spectra, as well as for an ongoing measurement program on short-lived Tm isotopes at ISOLDE/CERN.

The nature of the dissipation of turbulent kinetic energy (TKE) is investigated experimentally in the wake of a diameter $D=0.58$m wind turbine exposed to several “flavours” of high-Reynolds-number free-stream turbulence (FST). For low- and moderate-intensity FST, an annular region of elevated normalised dissipation, $C_{\varepsilon}$, develops in the outer wake, coinciding with a ring of enhanced turbulence intermittency at both small ($\ell \leq \lambda$) and large ($\ell \geq D$) scales, where $\lambda$ is the Taylor microscale. In the blade-tip region, $C_{\varepsilon}$ scales with $\sqrt{Re_D}/Re_{\lambda}$, where $Re_D$ is a global Reynolds number and $Re_{\lambda}$ a local turbulent Reynolds number based on $\lambda$-a scaling indicative of dissipation being out of equilibrium with the inter-scale flux of TKE in the inertial range of the energy cascade. This non-equilibrium regime is interpreted in light of the observed intermittency : large-scale intermittent events (i.e., low-wavenumber perturbations/“kicks”), driven by persistent tip-vortex and tip-shear-layer dynamics under low-intensity FST, require a finite time to cascade down to the dissipative scales, thereby introducing an imbalance between the inter-scale energy flux and dissipation. At the wake centreline, by contrast, $C_{\varepsilon}$ remains approximately constant with streamwise distance, reflecting either classical Kolmogorov-type equilibrium turbulence or balanced non-equilibrium turbulence, with intermittency confined to the small scales. Under high-intensity FST, large-scale intermittency is suppressed, consistent with the erosion of tip-vortex structures, and no comparable scaling for $C_{\varepsilon}$ could be identified using a single turbulent Reynolds number for these cases, where two similarly intense streams of turbulence, but of different origins, are adjacent to one another.

In the heavy-fermion metal YbRh2Si2, quantum criticality at a suppressed antiferromagnetic order is governed by the interplay of local magnetic moments and itinerant conduction electrons. We demonstrate how this can be investigated by a new experimental approach that enables the observation of electron spin resonance (ESR) across a broad range of frequencies and fields at very low temperatures. This allowed us to cover a large part of the phase diagram from the paramagnetic Fermi-liquid phase to the phase with antiferromagnetic order and including the quantum-critical regime. Both the ESR g-factor and the linewidth present distinct behaviors in these three regimes, providing further insight into the physics across a quantum critical point. Notably, when cooling down at a field directly towards the quantum critical point, both g-factor and linewidth continuously decrease. Furthermore, we observe a very good matching of the g-factor behavior upon field-tuning and temperature-tuning towards the quantum-critical point. We analyze and discuss the results in the context of present theories on ESR in strongly correlated electron systems.

Spatial confinement and chemical screening fundamentally reshape the non-equilibrium phase behavior of autophoretic active particles. Here, we present a systematic study mapping the collective dynamics of self-propelled particles governed by chemo-attractive translational forces ($\mu_t < 0$) and chemo-repulsive rotational torques ($\mu_r > 0$) across varying screening parameters $\kappa$, torque magnitudes $\mu_r$, and boundary condition coefficients $\Lambda^c$. Beyond standard chemotactic macro-phase separation and dynamic clustering, we report the emergence of novel boundary- and screening-induced bubbly phases, classified into boiling and bursting bubbles. Using a metric triad of steady-state cluster fraction $\langle S \rangle$, temporal fluctuation magnitude $\sigma_S$, and coordination number $\langle Q \rangle$, we draw phase diagrams to demarcate phases for no-flux boundaries ($\Lambda^c = 1$) and chemically permeable interfaces ($\Lambda^c = 0$) . Increasing chemical screening ($\kappa$) systematically suppresses long-range attraction, driving sequential phase transitions from macro-scale collapse toward bubbly states, dynamic micro-clusters, and homogeneous gas phases, while simultaneously inducing aggregate shape anisotropy. These findings provide predictive design rules for controlling active assembly and transport in microfluidic environments.

Strongly anisotropic non-Kramers rare-earth ions combine giant longitudinal g-factors with a vanishing transverse component imposed by time-reversal symmetry, a combination that makes their spin transitions exquisitely sensitive to the orientation of the applied magnetic field. We show that this sensitivity carries a dual identity: it is simultaneously an overlooked decoherence channel and the basis for a spin-coherence-based angular probe. Using pulsed electron paramagnetic resonance at X-band, we report the first measurements of the quantum coherence of Tb$^{3+}$ in a native-doped CaWO$_4$ crystal (15 ppb) and map the Hahn-echo coherence time $T_2$ as a function of temperature (2 to 10 K) and resonant field ($10^3$ to $10^4$ G). A parameter-free model combining spin-lattice relaxation, instantaneous diffusion and spectral diffusion from all independently quantified impurities overestimates $T_2$ by an order of magnitude at low temperature and wrongly predicts the field dependence of $T_2$, inconsistent with the observed monotonic decrease of $T_2$ with $B_r$. A two-parameter extension, including dynamical angular fluctuations of the crystal axis, reproduces the full dataset across multiple setups and laboratories. Two controlled experiments nominally identical except for different mechanical configuration of the setup establish the mechanical origin of the dominant contribution. The two-parameter extension corresponds to an angular amplitude noise spectral density of overall order 36 n{\deg}/$\sqrt{Hz}$ from global external vibrations (ranging from 10 to 66 n{\deg}/$\sqrt{Hz}$ depending on the exact setup mechanical configuration) estimated at $\sim$ 2.5 kHz plus a temperature-dependent contribution assumed to come from local phonon-driven angular jitter. It identifies and highlights a decoherence pathway of practical relevance to any anisotropic solid-state spin system.

Black holes, with their strong gravitational fields, provide an important testing ground for theories of gravity beyond General Relativity. Among the many proposed alternatives, considerable recent work has focused on scalar-tensor theories in which a scalar field couples to higher-curvature terms. The black hole solutions that arise in such theories can differ significantly from the Schwarzschild and Kerr solutions of General Relativity. Characteristic properties of these black holes include instabilities, shadows, and gravitational wave spectra, which can be used to constrain the couplings of the underlying theories.

A gravitational-wave (GW) detector records a linear mixture of local couplings under the assumption that extra polarizations enter geodesic deviation. Vacuum general relativity (GR) admits two transverse-traceless (TT) amplitudes. It remains to determine the largest set of couplings that can sit in that mixture, and whether a tensor-only arm-length result selects GR uniquely. The little group \(E(2)\) of a null four-momentum classifies the strain amplitudes \(\mathbf{p}=(p_{+},p_{\times},p_{x},p_{y},p_{b},p_{\ell})\), commonly written \(h_{P}\), which determine the electric tidal tensor along a ray. Geodesic deviation, recorded as differential arm length, therefore contains only those \(p_{P}\) that enter that tensor. Lorentz mixing at helicity \(\pm1\) supplies a gravito-magnetic (GEM) field that does not enter the tidal tensor. if GEM field is static, that does not propagate as a wave. A time-varying helicity-\(\pm1\) current sources a GEM wave that enters the mixture as a coupling, to be isolated by its measured quantity. For a radiation-zone wave that depends only on retarded time, \(\betag_{\perp}=\hat{\mathbf{k}}\times\partial_{t}(p_{x},p_{y})\), and the eight couplings are \(\Pvec=(\mathbf{p},\betag_{\perp})\). Here we adopt unified field equations on \(\mathbf{p}\) to clarify the origin of each component of \(\Pvec\); the measured quantity of each coupling then isolates the polarizations in that mixture, which favors identification of distinct polarizations and model tests.

High-resolution cross-correlation spectroscopy (HRCCS) is a powerful technique for detecting molecules whose individual spectral lines are too weak to be identified directly, but its sensitivity is limited by the availability and quality of high-resolution opacity data. Many molecules of atmospheric and astrobiological interest lack complete line lists, restricting traditional template-based searches. In this work, we use Titan as a controlled testbed to develop and validate a new cross-section-based methodology for HRCCS template construction. We analyse K-band CRIRES+ observations of Titan (1.99 - 2.48 {\mu}m), and compute cross-correlation functions using both line-by-line and cross-section-based templates. Our analysis recovers known hydrocarbons such as methane (CH4) and acetylene (C2H2), and yields the first HRCCS detection of ethane (C2H6) with a significance of SNRpeak = 5.17 {\pm} 0.07. The ethane detection was made possible exclusively through cross-section-based templates, as no high-resolution line list currently exists for this molecule. These results demonstrate that cross-section-based template construction is a practical and powerful strategy for extending HRCCS to molecules that currently lack reliable line lists, and establish Titan as a benchmark for calibrating molecular detection techniques that can be applied to both solar system and exoplanet atmospheres. Future applications of this approach to other ground-based high-resolution spectrographs, as well as to JWST's highest-resolution modes and next-generation facilities such as the ELT, could significantly expand the inventory of molecules detectable in planetary atmospheres.

Three-dimensional linear stability of a weak, compressible shock layer is analyzed using exact Reynolds-number-independent forms of the governing base flow and linear stability equations, obtained by rescaling the problem onto the shock's own natural viscous length and time scales. For any base flow with zero transverse velocity components, depending on the shock-normal coordinate alone, the resulting eigenvalue problem is proved covariant under rotation of the transverse wavenumber vector: the eigenvalue depends on the transverse wavenumbers only through their magnitude, and every three-dimensional eigenmode, throughout the continuous spectrum, decouples exactly into two independent constituents. A two-dimensional, in-plane acoustic–entropy branch carries the disturbance's entire dilatation, pressure and thermodynamic coupling, and a one-dimensional, out-of-plane purely solenoidal vortical branch is governed by a single scalar shear-diffusion equation. Three-dimensional shock-layer stability is thus resolved exactly into its compressible and vortical constituents, confirmed to machine precision by the Grosch–Salwen far-field companion spectra. Both decoupled operators carry an explicit damping term growing with the total transverse wavenumber, so a three-dimensional, oblique disturbance is never less stable than the two-dimensional disturbance sharing its in-plane wavenumber alone; the vortical branch, whose spectrum is controlled directly by this term, is shown to retreat from the imaginary axis in close proportion to the wavenumber squared. This finding is an exact, Reynolds-free analogue of Squire's classical theorem, with the stabilizing role of Reynolds number played here by the transverse wavenumber magnitude.

We prove that the electric field generated in three-dimensional space by finitely many point charges of one sign has only finitely many equilibrium points, thereby answering a 1969 question of Morse and Cairns (restated by Eremenko in 2008 and, as a conjecture, by Shapiro in 2015). More generally, for nonzero charges $q_i$ of possibly mixed signs at distinct sites $\mathbf a_i\in\mathbf{R}^3$, we show that the Coulomb field has at most $2^{N-4}(N-1)(9N^2+9N+10)$ equilibria in the region where $S(\mathbf x):=\sum_iq_i|\mathbf x-\mathbf a_i|^{-3}$ does not vanish. Of course, for charges of one sign, $S$ is nonzero everywhere. The proof rules out curves of equilibria using algebraic geometry and complex analysis on an associated complex curve, and then applies a Bézout count to obtain a quantitative bound. Well-known examples show that, in the mixed-sign case, the Coulomb field can vanish on curves contained in the zero set of $S$.

The properties of compact stars with a strong first-order phase transition to quark matter and with an additional fluid of self-interacting bosonic dark matter (DM) are studied. We find that the inclusion of DM changes the stability of mass-radius configurations relative to the naive one-fluid criterion. For compact star configurations with similar masses and different radii, so-called twin stars, the presence of DM removes the unstable segment between the hadronic and the hybrid branch, so that the stable mass-radius sequence becomes continuous after the onset of the phase transition to quark matter. We furthermore find stable ultra-compact objects (UCOs), defined by a total compactness $C = M_\text{tot}/R_\text{grav} \ge 1/3$. We observe two distinct classes of UCOs: a DM-halo class with $f_\text{DM} \gtrsim 0.9$, and a DM-core class at $f_\text{DM} \lesssim 0.02$. The two classes can be separated by the surface redshift of the normal matter, which reaches $z = 0.734 - 0.808$ for the DM-core and extreme values of $z = 2 - 3.468$ for the DM-halo class, allowing us to use their redshift as a way to disentangle the two classes. Furthermore, the high redshift of the DM-halo UCOs could serve as a signal of DM since their redshift exceeds the Buchdahl limit of $z = 2$. Finally, we find hybrid star solutions of “ultimate twins" with similar masses and visible radii, but different dark matter content, leading to different tidal deformabilities and surface redshifts. Future X-ray and gravitational measurements of ultra-compact neutron stars with radii and masses outside the allowed neutron star range can thereby probe the presence and the properties of DM in addition to a first-order phase transition to quark matter.

This work investigates the unsteady aerodynamic performance and vortex dynamics of avian-inspired flapping wings using our in-house sharp-interface immersed-boundary solver, VorteXdyn. A falcon-inspired body-wing model based on NACA 4312 profile is employed to examine aerodynamic force generation and vortex evolution during steady forward flapping flight at Reynolds numbers of 2,500, 5000, and 10000 and Strouhal numbers of 0.18, 0.225, and 0.27. The influence of physiologies is examined using three wing configurations: a simplified wing without distinct feather structures, a wing incorporating feather-like structures (serrations) along the trailing edge, and a geometrically detailed wing incorporating multiple feather layers consisting of primary, secondary, and median feathers over its span. The aerodynamic performance of these configurations is quantified using the temporal profiles and time-averaged characteristics of the lift and drag coefficients and the lift-to-drag ratio. The associated vortex dynamics are characterized through the formation and evolution of multiple leading-edge vortices (LEVs), their spanwise coherence, circulation, characteristic size, and persistence over the wings' surfaces. Particular emphasis is placed on the spanwise development of the LEVs from the root to the wingtip, their interactions with the tip vortices, and the resulting wake evolution during the downstroke. Our results demonstrate that increasing geometric fidelity modifies aerodynamic force production, the formation and evolution of multiple LEV structures, vortex-vortex and votex-wing interactions, and wake topology. These findings provide insight into the aerodynamic role of feather morphology in three-dimensional flapping flight for the aerodynamic design of efficient bio-inspired flapping-wing micro air vehicles.

We report the measurement of the critical velocity of a paraxial superfluid of light flowing past a mobile impurity, realized using a dual-wavelength optical scheme in a warm vapor of $^{87}\text{Rb}$. The superflow is two-dimensional by nature, while the impurity features a highly repulsive potential, a large tunable radius, and a finite mass. The critical velocity is evaluated by simultaneously monitoring three complementary observables: the number of vortices nucleated in the wake of the impurity, the hydrodynamic force exerted on it, and its displacement velocity relative to the flow. These independent methods yield consistent values for the critical velocity, which remain systematically below the Landau criterion and decrease with the impurity radius, in qualitative agreement with theoretical expectations. Our results demonstrate the reliability of these combined diagnostics, establishing paraxial superfluids of light as a versatile platform for exploring the impurity problem in quantum hydrodynamics.

Using a comprehensive set of three-dimensional, high-resolution direct numerical simulations, we investigate the existence of a dissipative anomaly in isothermal, homogeneous, isotropic compressible turbulence driven by solenoidal forcing. We find that the total kinetic-energy dissipation rate, as well as its solenoidal and dilatational components, approaches finite asymptotic values with increasing Reynolds number $Re$. The normalized mean dissipation rates collapse onto two distinct branches: one corresponding to the subsonic and transonic regimes, with root-mean-square Mach numbers ($M_{\rm rms}\lesssim 1$), and another to the highly supersonic regime, with ($M_{\rm rms}\ge 3$). This two-branch Mach-number dependence is most pronounced for the total kinetic-energy dissipation. For the solenoidal and dilatational dissipation rate components, the dependence on $Re$ depends in addition on the specific choice of the integral scale and root-mean-square velocity. Despite grid resolutions of up to $2048^3$ points the Reynolds numbers accessible are not sufficiently large to distinguish conclusively between a weak and a strong dissipative anomaly. We substantiate these findings using three complementary approaches: (i) a detailed analysis of the mechanisms responsible for dissipation generation, based on the corresponding dissipation-rate balance equations and their individual production terms; (ii) an investigation of precursors of anomalous dissipation using the Duchon–Robert framework extended to compressible flows; and (iii) a geometrical characterization of regions of intense dissipation. Taken together, these analyses provide consistent evidence for the existence of a dissipative anomaly in isothermal compressible turbulence, while leaving its precise weak or strong character unresolved.

We prove that, for an $n$-qubit system of dimension $d=2^n$, every state satisfying $\operatorname{Tr}(\rho^2)\le 1/(d-a_\ast)$, with $a_\ast=0.458327\cdots$, lies inside the stabilizer polytope and is therefore magic-free. Combining this result with general geometric properties of high-dimensional polytopes, we establish quantitative estimates for the Hilbert–Schmidt inradius and volume radius of the stabilizer polytope, and use them to characterize the typicality of magic in random induced states obtained by tracing out a $k$-dimensional subsystem from a $d\times k$-dimensional Haar-random pure state. We prove a sharp phase transition in the probability of such states having magic, whose transition dimension $k_\star$ is bounded between $\Omega(d^2/\log^2d)$ and $\mathcal{O}(d^2)$. We further prove that the number of facets of the stabilizer polytope lies between $\exp[\Omega(d^2/\log^2 d)]$ and $\exp[\mathcal{O}(d^2\log^2 d)]$ employing a result of Bourgain and Milman in convex geometry, substantially improving upon the previous quasipolynomial lower bound and implying that doubly-exponentially many linear inequalities in the number of qubits are required for an exact description of the magic-free region. Overall, our results reveal the near-extremal geometry of the stabilizer polytope and provide a quantitative foundation for understanding the typicality, robustness, and detectability of magic.

Hydrogen tunneling is central to many chemical and biological processes. Herein, we introduce the vibronic-SHAKE (V-SHAKE) approach to comprehensively sample energy-conserving molecular configurations enabling hydrogen tunneling. A constrained form of nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) dynamics, where the tunneling hydrogen nucleus is quantized, is used to sample geometries in the diabatic seam space corresponding to the intersection of the reactant and product NEO-DFT diabatic vibronic surfaces. V-SHAKE is applied to hydrogen and deuterium tunneling in 4-cyanobutanolate and Z-4-hydroxybut-3-en-2-one. The vibronic coupling is found to vary significantly in the diabatic seam space, mainly due to changes in the donor-acceptor distance. The reaction coordinate and gradient of the vibronic coupling at the minimum energy crossing point are nearly orthogonal and are dominated by motions stabilizing the product relative to the reactant or decreasing the donor-acceptor distance, respectively. V-SHAKE provides fundamental insights and validation for assumptions underlying rate theories.

The exceptional gravitational-wave event GW231123 exhibited unusually high total mass and near-extremal spins, while revealing strong waveform-dependent variations in its inferred source properties. Several follow-up studies favored a lensed interpretation, with reduced waveform systematics, raising the question of whether apparent lensing evidence could arise from waveform modeling systematics. Using GW231123-like numerical-relativity injections, we find that lensing degrees of freedom can reduce waveform discrepancies and yield non-negligible lensing support. However, this behavior is not systematic across injections, and cases where all waveform and lens models favor lensing while reducing waveform discrepancies are not found.

Maximum-likelihood decoding (MLD) achieves the minimum logical error rate of stabilizer codes under known i.i.d. Pauli noise, but its exact evaluation is \#P-hard. Practical pipelines therefore approximate MLD by minimum-weight decoding (MWD), retaining only the lowest-weight recovery per syndrome and discarding the coset degeneracy. The minimum-weight search is in turn implemented by stochastic solvers. We introduce approximate maximum-likelihood decoding (AMLD), a black-box framework that recycles the candidate samples discarded by stochastic inner decoders into a per-class truncated free-energy estimator. For every logical class represented in the candidate pool, the estimator is provably bounded below by the exact free energy and above by the empirical minimum weight. AMLD returns the logical class minimizing the estimated free energy with linear classical overhead. In SA-based Ising-decoder benchmarks, AMLD closes up to $83\%$ of the MWD–MLD threshold gap across the toric and color codes under bit-flip and depolarizing noise. The largest threshold improvement, from $17.28\%$ to $18.62\%$, occurs on the $6.6.6$ color code under depolarizing noise. We further demonstrate AMLD on the $[[144,12,12]]$ bivariate-bicycle code, whose bit-flip decoding problem has a hypergraph structure. This application requires neither matching-based enumeration nor code-specific tensor-network contraction. At $p=0.05$, AMLD reduces the logical error rate by $13\%$ relative to MWD evaluated on the same BP-OSD candidate pool.

Random quantum objects are powerful resources for quantum information processing, yet exact Haar randomness is costly and typically unnecessary. We introduce an explicit sparse commuting circuit ensemble on $n$ qubits that reproduces low-order Haar moments in the stringent relative-error sense. The circuit consists of a sparse Clifford phase layer followed by independent single-qubit Clifford gates. Acting on a simple product state, the resulting ensemble forms $\epsilon$-approximate projective $2$- and $3$-designs in relative error, with the required logarithmic interaction degree being asymptotically optimal within this circuit family. It admits an ancilla-free implementation of quantum depth $O(\log(n/\epsilon))$ on an all-to-all architecture, as well as an adaptive constant-depth implementation—in fact, depth seven—using $O(n\log(n/\epsilon))$ ancilla qubits. Departing from existing shallow-design paradigms, our analysis exploits the intrinsic moment structure of commuting phase circuits; at third order, this requires a new block decomposition and combinatorial analysis that also suggests a route toward higher-order shallow designs. Our results show that precise Haar-like statistics can emerge from sparse commuting dynamics with remarkably low quantum resources, with applications to randomized characterization, quantum metrology, quantum algorithms, and many-body physics.

We identify non-Gaussian corrections to the central limit theorem for the global and local laws, i.e. for Stieltjes transform of the empirical eigenvalue distribution of a large real symmetric or complex Hermitian Wigner matrix. We find that in the real case the rate of convergence in this CLT is substantially slower than in the complex case. For the proofs, we develop new static and dynamic hierarchies of extended cumulants up to any order that can be analyzed by a refinement of the zigzag strategy.

The fusion excitation function for $^{21}$O + $^{12}$C was measured for the first time and compared to the fusion of less neutron-rich isotopes. The impact of valence neutrons in the d$_{5/2}$ shell on the fusion excitation function is examined. The experimental data manifest a clear dependence of the extracted barrier height, V$_B$, and barrier position, R$_B$ on neutron excess. To assess the role of dynamics the experimental data are compared with both density constrained frozen Hartree Fock (DCFHF) and density constrained time-dependent Hartree Fock (DCTDHF) theories.

The localized electromechanical phase transition in dielectric elastomers involves complex moving boundaries and severe electrostatic fringe fields driven by high-curvature interfaces. Traditional phenomenological models fundamentally underestimate the configurational forces by completely ignoring the in-plane electric field components and geometric singularities. Here, we present a asymptotically exact geometric framework to resolve the domain wall propagation. By mapping the highly deformed current configuration to a regular parametric strip via conformal mapping, the field singularities are algebraically eliminated. An exact integration by parts directly translates the higher-order geometric metric into a rigorous topological mass term, projecting the global conformal electrostatics into a non-linear $\sigma$-model Lagrangian density. Incorporating the Gent strain-stiffening model, the continuous translation symmetry yields a Hamiltonian first integral. Eigenvalue analysis and Bogomol'nyi-Prasad-Sommerfield (BPS) bound calculations prove that the domain wall emerges strictly as an asymmetric heteroclinic orbit connecting two saddle points. This geometric resolution provides exact analytical scalings for the localized interface energy, domain wall thickness, and asymptotic decay lengths, eliminating all phenomenological parameters and offering a deterministic paradigm for geometric instabilities in active soft matter.

Coherent light-matter interfaces controllably modifying the state of a photon upon interaction with a stationary qubit are a key resource for implementing deterministic entangling gates for optical quantum technologies. This requires a one-to-one mapping between the state of the scattered photon and that of the embedded qubit. Here, we present an experimental signature of such a bijection by leveraging the spin-induced Kerr rotation present in a low-noise charged quantum dot-micropillar cavity device. Through time-resolved polarization measurements, we project the electron spin to one of its eigenstates with $95\pm2\%$ fidelity with a single reflected photon detection, and follow the subsequent spin relaxation through the detection of a second reflected photon. We demonstrate that, after a transient regime governed by the trion radiative lifetime, two orthogonal polarization states can be produced, each associated to a given spin eigenstate. While the current results are limited by a timescale competition between electron spin relaxation and trion radiative lifetime, they could be improved using hole spins displaying increased relaxation times. Our work paves the way towards deterministic logic gates exploiting this one-to-one mapping between a spin and the polarization of a scattered photon.

We develop a geometric framework for conservative mechanical systems modeled on conic pseudo-Finslerian manifolds, extending classical Riemannian and semi-Riemannian mechanics to anisotropic geometries. In this setting, we define a Jacobi-type pseudo-Finslerian metric and demonstrate that motions of fixed energy correspond, up to reparametrization, to its geodesics. We establish energy and momentum conservation under suitable symmetry assumptions via a generalized Noether's theorem, and provide new analytical criteria for the completeness of both the Jacobi metric and the mechanical system. Furthermore, for mechanical systems with holonomic and non-holonomic constraints, we prove the existence and uniqueness of reaction forces satisfying d'Alembert's principle. These results provide a unified geometric foundation for constrained and unconstrained anisotropic dynamics.