Papers

A leading theory for hot Jupiter formation is high-eccentricity migration, in which planets are born at large separations and excited to near-unity eccentricities, creating extreme tidal dissipation that shrinks and circularizes their orbits. The most direct evidence for this scenario is the detection of highly eccentric planets caught in the act of migrating. Although such planets are rare, the recent discovery of TIC 241249530 b - a super-eccentric ($e = 0.94$) and retrograde hot Jupiter progenitor - presents an ideal case study to explore high-eccentricity migration and the accompanying tidal physics. In this paper, we examine the migration history of TIC 241249530 b and test two different models of tidal dissipation: equilibrium tides and chaotic dynamical tides. We show that TIC 241249530 b's properties can be explained by high-eccentricity migration triggered by von Zeipel-Lidov-Kozai oscillations induced by the observed distant binary star in the system, but only if the dominant tidal dissipation takes the form of equilibrium tides. Chaotic dynamical tides fail to explain the system because they require the planet to have migrated substantially closer to its star.

We construct an $R^2$ cosmology based on scale invariance and unimodular gravity in which the Goldstone boson of dilatations establishes a predictive connection between inflation and dark energy. The conservation of the corresponding Noether current confines the inflationary trajectory to a one-dimensional orbit in field space, freezing the Goldstone direction and rendering the inflationary dynamics effectively single field. After inflation, the system settles on a Minkowski vacuum manifold. The unimodular integration constant, already present in the full theory but negligible during inflation, lifts the flat Goldstone direction, generating an exponential potential for the canonically normalised field. Unlike in phenomenological quintessence models, its slope is not a free DE parameter, but is instead fixed by the field-space geometry inherited from the inflationary attractor. This geometry excludes matter-era tracking and, throughout the inflationary viability region, places the field in the thawing regime compatible with accelerated expansion. For a representative inflationary benchmark, the now pseudo-Goldstone field remains frozen until close to the present epoch, yielding DE equation-of-state parameters $(w_0,w_a)\simeq(-0.992,-0.011)$ in the CPL parametrisation $w(a)=w_0+w_a(1-a)$. More generally, the same parameter controlling the inflationary spectral tilt also fixes the slope of the DE potential, leading to a tight correlation between early- and late-Universe observables: increasing the thawing signal lowers $n_s$, while values of $n_s$ closer to unity drive the model towards $\Lambda$CDM.

JWST surveys reveal a greater space density of high-redshift UV-bright galaxies than predicted by conventional galaxy formation models. We present results from a $L=100$ cMpc cosmological simulation evolved to $z=5$ with a variation of the COLIBRE galaxy formation model that adopts a density-dependent stellar initial mass function (IMF), such that stellar populations formed from dense gas are born with a top-heavy IMF. Crucially, heavy element and dust yields, and supernova feedback energetics, are self-consistently adjusted to the changing IMF. We model UV/optical emission (including nebular emission) from galaxies and its attenuation by dust. By allowing a significant fraction of high-redshift star formation to proceed with a top-heavy IMF, the rest-frame far-UV luminosities of early galaxies are elevated by up to a factor of $\simeq4$ with respect to the fiducial COLIBRE L100m6 simulation, which assumes a universal Chabrier IMF. This enables the formation of galaxies with observed brightness up to $M_{\rm UV} \simeq -20$ at $z=15$ (c.f. $M_{\rm UV} \simeq -18.5$ in the fiducial simulation), illustrating the potential of star formation with a top-heavy IMF to alleviate tensions with JWST data. Later, the boost in far-UV emission is partly offset by attenuation due to increased dust surface densities from i) additional dust grain ejection from core-collapse supernovae and ii) efficient grain growth promoted by more metal-rich interstellar gas. The simulation reproduces the $z=5$ galaxy stellar mass function and rest-frame optical luminosity function with comparable accuracy to the fiducial simulation, and both simulations exhibit UV continuum slopes that are consistent with JWST observations.

We show that no stabilizer state in a discrete realization of a local quantum field theory can flow in the continuum to the vacuum or to any state that resembles the vacuum at short distances. The argument rests on the fact that the entanglement spectrum is flat for stabilizer states but non-flat for cyclic and separating states in a local QFT as a consequence of the type III$_1$ nature of the von Neumann algebras associated with arbitrary subregions. Our result implies that simulating physically relevant QFT states necessarily requires quantum resources beyond stabilizer states and Clifford operations. We comment on the implications for holography.

Free-floating planet-moon pairs (FFPMs) may form when planetary systems experience close stellar encounters. We quantify moon-retention probabilities and ejection cross sections for planets and planet-moon systems subject to close stellar encounters, and to dentify the orbital conditions most favorable for forming FFPMs. We conducted extensive numerical scattering experiments to compute ejection cross-sections on a grid in the orbital separation of the planet ($a_p$) and the moon ($a_m$), marginalizing over all other initial parameters. The simulations tracked both planetary and lunar fates across a wide range of encounter geometries. FFPM cross-sections depend on the moon semi-major axis $a_m$, decreasing gradually until a drop near $a_m\sim0.45\,R_H$. Planet-only ejection cross-sections instead decline steadily with $a_p$. Io-like moons ($a_m = 0.008\,R_H$) survive, but survival falls beyond $0.4\,R_H$. Compared with planet-planet scattering, stellar encounters preserve moons more effectively across all separations. Surviving moons at $a_m\lesssim 0.4\,R_H$ retain near-circular, low-inclination orbits, while wide-orbit moons show stronger dynamical excitation. Applying our cross-sections to the microlensing system MOA-2011-BLG-262Lb suggests possible progenitor semi-major axes between $a_p \sim 1.3$\,au and $\sim 5.9$ au for a solar-mass host. Uncertainties remain large and the probability distribution is flat but we prefer $a_p \sim 5.2$\,au. Stellar-encounter ejections constitute a viable channel for producing FFPMs whose orbital properties differ from those formed by planet-planet scattering. Moon retention and orbital excitation provide promising diagnostics of ejection history. Current and upcoming microlensing and direct-imaging surveys may be capable of detecting Galilean-mass moons around rogue planets, offering new tests of dynamical formation pathways.

We study a dark non-Abelian gauge sector with GeV-scale confinement. The dark sector is assumed to couple only feebly to the Standard Model, while its low-energy spectrum may contain long-lived flavor-diagonal dark pions. Signals of these states are particularly well suited to the Electron-Ion Collider (EIC), where the absence of a hard trigger requirement and the capability to record soft final-state particles offer a complementary probe of dark hadronization dynamics. We present a benchmark portal construction, discuss the mixing between an axion-like mediator and dark pions, and identify the resulting displaced-decay signature.

Real-space geometric curvature in magnetic thin films introduces a controllable mechanism for tailoring the pathways of superconducting steady-state Josephson currents via the proximity effect. We present a generalized Green's function method for calculating diffusive transport in arbitrarily curved surfaces, and show how the competing mechanisms of curvature and distance regulate conversion between different superconducting pairings in a proximity-coupled ferromagnet. We show how this dictates the distribution of current density, with the possibility of curvature-controlled current density manipulation, induced eddy currents and current reversal.

We investigate the production of particles much heavier than the characteristic scale of a cosmological first-order phase transition through collisions of highly boosted bubble walls. Using the scalar order-parameter field, we derive the ultraviolet behavior of its Fourier-space profile for both elastic and inelastic collisions. In the regime $\chi\equiv\omega^2-\mathbf{k}^2\gg M_h^2$, we find the universal result $\tilde{\phi}(\chi) = -2V^\prime(2v_\phi)\chi^{-2}+O(\chi^{-3}),$ implying that the spectral density scales as $F(\chi)\propto [V^\prime(2v_\phi)]^2\chi^{-4}$. Thus, heavy-particle production is localized near the instant of collision and, at leading order, depends on the scalar potential only through $V^\prime(2v_\phi)$. We verify this behavior using high-precision numerical solutions of the trapping equation, carefully suppressing spectral leakage from the finite integration domain, and obtain agreement over a broad ultraviolet range. We then derive analytical production rates for general heavy-particle thresholds and for fermion pairs, together with their cosmological number density and yield. Finally, we extend the analysis to $(3+1)$ dimensions and incorporate the finite bubble radius, finding an order-one suppression relative to the parallel-wall approximation. Our results revise the ultraviolet scaling used in previous treatments and have direct implications for superheavy dark-matter production and baryogenesis from bubble collisions.

Quiescent black hole low-mass X-ray binaries provide a key setting for probing accretion physics at low luminosities. A0620-00, the archetypal system in this class, has remained in X-ray quiescence for decades and exhibits complex optical variability, yet the long-term behaviour of its accretion flow remains poorly understood. Here, we report an analysis of long-term optical monitoring of A0620-00 from ZTF, LCO, and ATLAS. The full dataset spans nearly two decades, with the ZTF light curve providing the primary $\sim 2760$-day baseline for the period analysis. We identify a superorbital cycle with a period of $P = 261.9 \pm 9.4$ d and a peak-to-peak amplitude of $\sim 0.2$ mag. The signal is recovered independently across all three surveys, and red-noise simulations indicate that it is unlikely to arise from stochastic variability alone. Furthermore, the relative occurrence of the passive and active quiescent states displayed by the system seems to depend on the superorbital phase, with passive states concentrated near the cycle minimum and active states more common near maximum. We find that, among the possible interpretations, retrograde nodal precession of a hot inner accretion flow might be able to explain the observed long-term modulation. In this interpretation, the periodic signal may arise from cyclic reorientation of the inner flow, which modulates the photometric contribution from the innermost regions. The inferred modulation period would correspond to a characteristic dynamical radius of $\sim0.13a$ ($\sim10^4~R_{\rm g}$), where $a$ is the binary semi-major axis, broadly consistent with the expected transition between the outer thin disc and the inner hot accretion flow.

Every two-qubit measurement basis composed of maximally entangled eigenstates can be transformed into the Bell basis via local unitary operations. For higher dimensions, in contrast, there exist inequivalent bases composed of maximally entangled eigenstates. Here, we provide a single-parameter family of two-qutrit maximally entangled measurement bases, and demonstrate that none of the bases are equivalent to each other under local unitaries. These bases are constructed from the continuous family of symmetric informationally complete sets of states in dimension three. By studying the local unitary bases that generate the family of two-qutrit maximally entangled measurement bases, we construct the first examples of wild error bases in the smallest dimension where these can exist. Finally, we discuss how distinct measurements in the family lead to differences in performance in several scenarios relevant in quantum information.

Motivated by growing observational indications that dark energy may be dynamical, we introduce Rolling Galileon gravity: a minimal shift-symmetry-breaking extension of the cubic Galileon in which the coupling coefficients are allowed to vary, giving rise to an evolving braiding strength. The full theory space, shown to be closed under field redefinitions, is characterised by two functions. We derive analytical conditions to satisfy three phenomenological requirements: i) a phantom-crossing equation of state at late times, ii) a positive integrated Sachs-Wolfe signature, and iii) absence of pathologies in the screened scalar force in cosmic voids. We show that these conditions are collectively satisfied by an increasing braiding strength relative to the kinetic sector. A Bayesian analysis of minimal Rolling Galileon models finds that they can satisfy the viability requirements i)-iii) whilst providing an acceptable fit to expansion-history data.

Tunable-range random graphs provide a versatile setting for exploring how geometry shapes quantum dynamics, with potential applications to analog quantum simulation, amorphous materials, and quantum information processing. Here, we investigate transport and entanglement growth in quantum many-body systems whose two-body interactions are encoded by such graphs, with sites as nodes and two-site couplings as edges. In particular, we study dynamics for quantum hopping on a random graph model with an effective range parameter. By tuning it, we interpolate between nearest-neighbor and all-to-all graphs, identifying a new localization regime that is robust in both single-particle and finite-density cases. Our results establish tunable-range random graphs as a platform for studying the interplay between connectivity, transport, and entanglement in quantum many-body dynamics.

Hydrodynamics is a macroscopic theory of long-wavelength dynamics around local thermal equilibrium. We develop attractodynamics, the analogous construction around a far-from-equilibrium attractor. Dynamics near a far-from-equilibrium attractor retains some non-hydrodynamic microscopic information, which attractodynamics systematically organizes. We move towards this general structure by starting in 0+1D, and consider a model for which an anisotropic far-from-equilibrium attractor solution is exactly known. This setting provides a clean benchmark in which the ideal attractodynamic equations and a leading transient residual extension can be compared directly with the full kinetic evolution. The resulting hierarchy gives an improvable description of near-attractor dynamics: the ideal theory captures evolution close to the attracting manifold, while retaining the leading off-attractor moments extends the regime of agreement until the finite truncation breaks down. This example identifies the ingredients needed for local 3+1D attractodynamics and for macroscopic attractodynamic theories not derived from an underlying kinetic description.

We introduce the transverse charge distribution as a spin-sensitive charge-flow probe for fragmentation and nucleon tomography. By measuring the angular distribution of net electric charge around a fragmenting quark, this observable relies entirely on the tracking of charged-particle directions and charge signs, strictly bypassing the need for calorimetric energy measurements. At leading twist, the distribution decomposes into an unpolarized charge monopole and a chiral-odd transverse charge dipole. We derive the operator product expansion of these distributions onto charge-weighted collinear moments: a charge monopole and a charge dipole governed by the Collins effect and couples directly to transversity. Applying this formalism to transversely polarized $p^\uparrow p$ collisions at RHIC, we show that charge weighting suppresses the unpolarized monopole background and causes the spin-dependent dipoles from oppositely charged hadrons to add coherently. This coherence strongly enhances the resulting azimuthal asymmetries, establishing a theoretically clean and experimentally precise track-only avenue for extracting transversity.

Floquet driving underlies Hamiltonian and gate engineering, and produces dynamical orders with no equilibrium counterpart. These phenomena are, however, transient in well-isolated systems. We show that statically coupled dissipative auxiliaries cool the system toward low-energy states of the Floquet-engineered Hamiltonian, stabilizing orders in the steady state. The excess energy density is controllably small at high drive frequency and weak coupling and is captured by rate theory based on Fermi's Golden Rule. We numerically confirm robust cooling in a Floquet-engineered transverse-field Ising chain in both phases, and demonstrate discrete time-crystalline order, with a period-doubled magnetization response, in the steady state of a long-range Ising chain. Our results provide a rare analytical handle on the steady states of driven dissipative systems.

Open-shell $3d$ transition-metal complexes pose a stringent challenge for correlated electronic-structure methods because near-degenerate metal $d$ orbitals give rise to competing spin states and oxidation-state-dependent charge transfer, while solvation can reshape the relative stability of the resulting electronic configurations. Building on prior demonstrations of open-shell sample-based quantum diagonalization (SQD) for light-atom systems and, separately, closed-shell SQD coupled to the integral-equation-formalism polarizable continuum model (IEF-PCM), we combine these capabilities through an open-shell SQD–IEF-PCM framework and apply it to the octahedrally coordinated $\mathrm{[Co(H_2O)_5CO_2]^{2+/3+}}$ complex, spanning the Co(III) singlet and quintet and Co(II) doublet and quartet states along a metal–ligand dissociation coordinate. Using quantum samples collected on an IBM Heron processor and active spaces of up to 50 qubits, SQD closely reproduces coupled-cluster and heat-bath configuration-interaction benchmarks within the same active spaces in both gas phase and implicit solvent; for high-spin quintet $\mathrm{[Co(H_2O)_5CO_2]^{3+}}$, it additionally resolves an avoided crossing associated with internal charge transfer that is absent in the corresponding singlet and in the lower oxidation state, while IEF-PCM suppresses this crossover by preferentially stabilizing the dissociation channel leading to neutral CO$_2$. To our knowledge, this work provides the first hardware demonstration of SQD for an open-shell 3d transition-metal complex with discrete metal-centered d-orbital chemistry, and the first application of open-shell SQD in a dielectric continuum, establishing a quantum-centric framework for transition-metal chemistry in regimes where spin-state energetics, charge transfer, and environmental response are strongly coupled.

It is well known that classical analog circuits and systems benefit from the powerful Thevenin and Norton theorems. These theorems enable rigorous and exact mathematical simplification and representation of the effect of the rest of a system on the part we want to focus on. We show that there are corresponding versions of a "Quantum Norton Theorem" that enable exact simplification and reduction for quantum circuits and systems, not just for one port, but also for multiport and even open quantum systems. We demonstrate the method on two level systems, chains, system environment partitions, Lindblad examples, and Grover search. The Grover examples show how the reduced network isolates the bright pole governing ideal search and exposes how diagonal disorder transfers spectral weight into dark poles that degrade performance. By partitioning finite dimensional Schrodinger dynamics into retained and eliminated sectors, we show that Gaussian elimination gives an exact reduced equation in which the eliminated subsystem appears as a dynamical self energy and a source term. In the circuit representation, these terms become the Norton admittance and Norton current source seen by the retained quantum port. The same Schur complement construction extends to multiport reductions, composite system environment partitions, density matrix dynamics, and Lindblad evolution in Liouville space.

We develop a deterministic, scale-resolved formulation of energy transfer in the three-dimensional incompressible Navier-Stokes equations based on an explicit triadic decomposition of the nonlinear term in Fourier space. Using a systematic dyadic localization of the velocity field, we derive an exact representation of the nonlinear energy flux across scales and organize it in terms of interactions between well-defined scale components. Under suitable smoothness assumptions, we obtain an absolutely convergent triadic expansion and quantitative bounds that distinguish local and nonlocal contributions in scale space. This framework provides a transparent and fully explicit description of how energy transfer is mediated by triadic interactions and how scale locality emerges as a structural property of the nonlinearity. Building on this formulation, we revisit the classical inertial-range picture of turbulence from a deterministic perspective. We show that, under a scale-invariant flux assumption, the Kolmogorov $-5/3$ scaling is formally consistent with the triadic energy-transfer mechanism at a structural level. The result does not rely on statistical assumptions, but instead follows from the structural properties of the Navier-Stokes equations combined with a scale-resolved representation of the energy flux. The present work thus provides a coherent synthesis of triadic interaction analysis, dyadic scale decomposition, and classical turbulence phenomenology, offering a deterministic framework that clarifies how Kolmogorov-type scaling constraints arise in the scale-resolved structure of the underlying equations.

Covariance dynamics on the positive-definite cone are commonly described by gradient flows, which encode dissipative relaxation but obscure the underlying phase-space structure. We construct a finite-dimensional Hamiltonian lift of covariance dynamics on Sym$^+_n$ equipped with the Bures–Wasserstein metric. The natural mechanical Lagrangian yields canonical momentum $\Pi=\tfrac12 L_\Sigma[\dot\Sigma]$, where $L_\Sigma$ is the Lyapunov operator, and explicit Hamiltonian $\mathcal{H}(\Sigma,\Pi) = 2{\rm tr}(\Pi\Sigma\Pi)+V(\Sigma)$. Adding Rayleigh dissipation recovers the Bures–Wasserstein gradient flow in the overdamped limit. For a spectral-floor and trace-control potential, the quadratic fluctuation Hamiltonian around the isotropic equilibrium separates trace and traceless modes; the baseline stiffness diverges as $(s-\nu)^{-2}$ as the equilibrium covariance approaches the floor. The construction identifies a conservative parent system for constrained Bures–Wasserstein covariance relaxation and fixes the local stiffness scale induced by the spectral floor.

One of us is a realist and believes that reality fundamentally consists in an external physical world, governed by laws of physics. Observers are either emergent (physicalism) or exist in addition to physical reality. The other defends a version of idealism and believes that reality fundamentally consists in first-person states, unembedded into worlds, on which laws of nature act directly. Shared "physical worlds" are merely emergent. The disagreement ultimately centers on the overall coherence of a formal model, known as Algorithmic Idealism, and its ability to resolve observer paradoxes, such as the Boltzmann brain paradox. Our aim is to confront this issue in the form of an adversarial collaboration, where we avoid misunderstanding each other, so that readers can see the true source of the disagreement, and decide for themselves. Our work represents an attempt to establish a form of exchange that may help overcome the fragmentation of scholarly communities in philosophy, physics, and elsewhere. The collaboration was conducted using a discipline-neutral template for theoretical adversarial collaboration developed in a companion paper [1], which we hope will be useful in other theoretical disciplines.

For centered bosonic Gaussian states, the covariance pullback of the Bures cometric differs from the lift-normalized classical covariance Fisher–Rao cometric by a state-independent bilinear form. Under the canonical identification of symmetric covariance covectors with $\mathfrak{sp}(2N,\mathbb{R})$, this form is the trace form and hence a fixed multiple of the Killing form. We prove that, within the normalized symmetric Petz family, Bures is uniquely selected by requiring such an additive state-independent symplectic correction. We determine the full Williamson-frame spectrum and show that a boundary stratum with $m$ pure modes has an $m^2$-dimensional cometric kernel isomorphic to $\mathfrak u(m)$, while the pure Gaussian orbit remains nondegenerate. For radial one-mode estimation, ideal heterodyne detection accesses the exact fraction $(\nu-\hbar/2)/(\nu+\hbar/2)$ of the SLD quantum Fisher information. Its vanishing boundary limit reflects finite heterodyne information relative to a divergent radial quantum Fisher information, not zero measurement information. Finally, a minimal Schur realization fixes the auxiliary inertia and identifies the second Schur complement as the true covariant vertical block.

Adversarial collaborations have become increasingly prominent in empirical science, where proponents of competing theories jointly design experiments capable of discriminating between them. Much less attention has been given to theoretical disputes that cannot readily be resolved by experiment. We distinguish theoretical adversarial collaboration from its empirical counterpart and propose a five-step template for conducting it. The procedure is designed to move disputants beyond familiar critique-reply-rejoinder exchanges by requiring sustained clarification, supervised steel-manning of the target theory, and the collaborative development and refinement of objections. Its condition for success is that the collaboration produce insights that substantially exceed the existing literature: either novel, non-obvious arguments show that the theory can meet the objections, or the theory's advocate concedes that new conceptual resources or revisions are required. The template was developed through reflection on, and implemented in, our own adversarial collaboration in philosophy and foundations of physics. We offer it here as a general framework for constructive and publishable progress in philosophy, foundational science, and other areas of inquiry in which disputes cannot readily be resolved by experiment.

In Friedmann cosmology with Einstein–Cartan torsion, the homogeneous torsion mode compatible with a separately conserved matter sector scales as $\Phi\propto a^{-3}$ and enters the Friedmann constraint as a stiff component of negative energy density, $-3\Phi^{2}\propto a^{-6}$. It has recently been claimed that this mode rescues the Hubble radius as an infrared cutoff for holographic dark energy, producing late-time acceleration and a phantom-divide crossing of possible relevance to DESI. We show that it cannot. With the Hubble cutoff the holographic density drops out of the deceleration parameter, and the only accelerating regime is the transient window $\bar{a}\leq a<4^{1/3}\bar{a}$ around the torsion bounce at $H(\bar{a})=0$; placing that window at observable redshifts would force the Hubble rate to vanish in our recent past, against the measured expansion history. Requiring a viable history yields nested upper bounds on $\Omega_{\Phi}\equiv(\Phi_{0}/H_{0})^{2}$: fitting the official DESI DR2 BAO likelihood gives $\Omega_{\Phi}<8.7\times10^{-4}$ ($95\%$ CL), the existence of a spectroscopically confirmed galaxy at $z=14.32$ gives $8.4\times10^{-5}$, the CMB gives $3.1\times10^{-10}$, and Big Bang nucleosynthesis gives $5\times10^{-24}$. Today's imprint on the dark energy equation of state, $|1+\omega_{0}|\leq2\Omega_{\Phi}/\Omega_{\Lambda}$, falls short of the DESI preference by two orders of magnitude at best and twenty-two at worst, and on the phantom side. In the Granda–Oliveros cutoff, torsion deepens rather than prevents the big-rip singularity. An appendix derives the Friedmann pair from the Einstein–Cartan field equations and shows that the $a^{-3}$ scaling is the kinematics of a diluting spin fluid; escaping the no-go requires breaking exactly that.

Efficient operation of neural networks has been linked to criticality in their underlying non-equilibrium excitation dynamics. However, obtaining experimental evidence of this conjecture remains challenging due to limited control and undersampling in biological systems. Here, we experimentally explore neural network criticality using an ultracold Rydberg gas as a highly controllable simulator. We highlight the similarity of the excitation spreading via Rydberg facilitation and the synaptic connection of spiking activity of neurons, giving rise to distinct absorbing and active phases. We systematically explore and resolve criticality criteria, including power-law scaling of excitation avalanches and the emergence of universal avalanche shape collapse. Crucially, we implement a controlled gain mechanism to compensate for atom loss, mimicking metabolic resource replenishment and stabilizing the system in a controlled non-equilibrium steady state. We find peak temporal correlations at the critical point and stochastic oscillations with dragon king avalanches in the active phase, consistent with predictions for systems orbiting criticality. Our work establishes facilitated Rydberg gases as a platform for investigating criticality, resource dynamics, and emergent oscillations in neural networks.

On flat superhydrophobic surfaces, droplet rebound is well described by a single inertio-capillary time scale, yielding a contact-time that is independent of impact energy. This single-mode response reflects the radial symmetry of flat-plate impacts. We demonstrate that an anisotropic geometric constraint, imposing a fixed spreading length along one axis, breaks this degeneracy and splits the rebound into a reciprocal pair of inertio-capillary modes. The fixed length also couples the contact-time to the Weber-dependent maximum spread, introducing an impact-energy dependence absent on the flat plate. We realize this constraint with grooved substrates, simulated using a non-ideal, entropic, multiple-relaxation-time lattice Boltzmann method and validated against the experiments of Chantelot et al. Extending their blob model from a single transverse scale to the reciprocal pair, we organize both modes through a geometric blob number and relate their time scales to the Weber number and groove width. We show that on non-wetting grooves the reciprocal modes are recovered directly, and explore the effects of finite wall affinity, using competition between the two modes to explain an observed two-branch structure in the contact-time response on mildly wetting, superhydrophobic grooves. Predictions tied to global energy balance reproduce cleanly across all conditions, while those tied to the details of the droplet's spread morphology are approximate but directionally correct. These results show that anisotropic confinement turns contact-time reduction from a question of accelerating a single rebound mode into one of selecting between conjugate inertio-capillary modes.

We derive rigorous Lorentz-covariant bounds on relaxation spectra directly from the analytic structure of retarded Green's functions in thermal quantum field theory, using only causality, unitarity, the Kubo-Martin-Schwinger condition, and Lorentz covariance, without reference to any specific dynamical model. A single rest-frame quasinormal pole is generically smeared into a continuum of excitations in boosted frames, with width set by the maximal signal velocity. We prove that the non-hydrodynamic gap $\Gamma_{\mathrm{gap}}$ transforms as $\tilde{\Gamma}_{\mathrm{gap}} \geq \Gamma_{\mathrm{gap}}/[\gamma(1 + v v_{\mathrm{max}})]$, and that the convergence radius of the hydrodynamic gradient expansion satisfies $\tilde{k}_c \in [k_c/\gamma(1+v v_s), k_c/\gamma(1-v v_s)]$ under a boost of velocity $v$. We verify the bounds by a numerical quasinormal-mode computation in the ${N}=4$ super-Yang-Mills plasma: the leading boosted pole moves deeper into the complex plane – the observed relaxation rate increases with boost velocity, in sharp contrast to naive time dilation – while respecting the bound throughout. The results apply non-perturbatively to the quark-gluon plasma, neutron star merger dynamics, and quantum critical systems.

Interactions between dark matter and dark energy can lead to an apparent phantom-crossing behavior that mimics the expansion history preferred by the latest cosmological observations from DESI baryon acoustic oscillations (BAO), Cosmic Microwave Background (CMB), and Type Ia supernovae (SNe Ia) data. In a previous paper [Khoury, Lin, and Trodden 2025 arXiv:2503.16415], we proposed a concrete particle physics realization of this idea, consisting of a strongly coupled dark sector in which a dark axion is coupled to dark baryons. In this paper, we investigate this idea further by comparing its predictions to the latest cosmological data. We implement the dark axion-dark baryon interaction model in a Boltzmann code and confront it with CMB, DESI DR2 BAO, and SNe Ia data. For the CMB+DESI DR2+DES-Dovekie combination, the best-fit model improves the fit relative to $\Lambda$CDM by $\Delta\chi^2=-14.48$. The preferred solution exhibits a non-monotonic dark-matter mass evolution: the mass decreases between matter-radiation equality and recombination, while increasing over the BAO/SNe-sensitive epoch, leading to an apparent phantom crossing in an effective dark-energy description. Interestingly, the same dynamics produces an Early Dark Energy-like energy injection near matter-radiation equality, but in the data-preferred region this component is too small to raise $H_0$ enough to substantially reduce the current tension.

Maxwell demons can convert knowledge into thermodynamic advantage by using measurement-acquired information. However, standard formulations of this problem assume that the demon has access to the system's entire phase space, an assumption that fails in many practical applications. Here, I address this problem by deriving fluctuation relations for general feedback-controlled systems, including such singular cases. The central quantity is a generalized notion of unavailable information: a portion of the acquired information that cannot be used to extract useful work. Finally, a Szilard engine with finite resolution illustrates these results. In the high-resolution limit, the information acquired by the demon may diverge, yet the extractable work remains finite. This work shows that this occurs because an equally divergent amount of information becomes unavailable.

To help address the escalating energy and latency demands of machine-learning workloads, we introduce a blueprint for an energy-efficient and fast thermodynamic computing stack that leverages stochastic analog processes in physical hardware. In this work, we focus on energy-based thermodynamic computing where the stochastic process is well described by Langevin dynamics with tunable energy potentials. The implementation of such potentials in physical hardware enables us to generate and sample from basic parameterized energy-based models. We demonstrate how to construct and train popular classes of machine learning models based on these hardware-native energy-based models, using the framework of probabilistic graphical models. We analyze the runtime and energy consumption of different models in this thermodynamic paradigm based on theoretical considerations and numerical studies. As a preliminary experimental realization of such hardware, we present our stochastic analog superconducting circuits driven by thermal noise. Together, these results outline a path toward energy-efficient thermodynamic hardware for probabilistic machine learning.

Geometric scaling was conjectured and observed at the ISR more than 50 years ago. We argue that it still holds at the LHC. We show that the dip-bump structures of the differential elastic cross-sections exhibit a remarkable regularity, namely that the ratio of the bump to dip positions is constant from 23 GeV to 13 TeV. Using crossing symmetry, analicity and the optical theorem, we identify the imaginary and real parts of the scattering amplitude. This allows us to compute the $\rho$ parameter and the ratio of the bump to dip values of the differential $pp$ cross-section. Finally, we discuss the energy dependence of the total elastic cross-section as compared to the total cross-section, and the violation of geometrical scaling outside the dip-bump region at the LHC.