Papers

We propose a new concept for quantum thermal logic gates – analogous to classical electronic logic gates – that exploit the heat current in a coupled quantum-dot system tunnel-coupled to metallic thermal reservoirs for logic operations in quantum circuits. We obtained a remarkable one-to-one correspondence with the structure of classical electronic logic gate circuits. An experimental setup is presented that demonstrates a realizable nano-electronic quantum circuit architecture for implementing such quantum thermal logic operations.

Strontium ruthenate (Sr$_2$RuO$_4$, SRO) has been the subject of intense research as a model quasi-two-dimensional metal with strong electronic correlations and potential exotic multi-component superconductor. Yet the nature of the superconducting state and its emergence remain debated, despite highly detailed knowledge of the normal-state electronic properties. Here we use diffuse neutron and x-ray scattering to uncover self-organized structural fluctuations on the picosecond timescale in SRO. We show that these nanoscale correlations do not originate from extrinsic disorder but rather involve cooperative displacements of oxygen atoms in the quintessential RuO$_2$ planes. Moreover, the observed displacement pattern is consistent with distortions due to incipient orbital order that we obtain in first-principles calculations, which suggests that orbital effects could play a pivotal role in the physics of SRO. Similar dynamic correlations may play a role in the physical properties of a wide range of prominent oxides with closely-related lamellar structures, such as the cuprates and nickelates.

We present a novel approach to eigenvector universality for generalized Wigner matrices. Our main consequences are asymptotic normality of joint eigenvector projections everywhere in the spectrum as well as, under a uniform subexponential decay assumption, convergence of the extremal process of eigenvector entries and the resulting Gumbel law for the largest entry of an eigenvector. In the case of smooth entries, we are able to obtain joint normality of an explicit growing number of eigenvector projections, and we are also able to obtain explicit rates of convergence in Kolmogorov distance, which are nearly optimal for a fixed number of projections. This is based on a new analysis of the Dyson vector flow which does not rely on the eigenvector moment flow.

Polymers are indispensable materials with daily applications ranging from electronics to medicine, owing to their versatility, which can be tailored by adjusting their chemical composition and architecture. The design space for these compounds is vast and governed by factors such as monomer classes, copolymer configurations (e.g., linear, branched, random, and alternating), chain size, stoichiometry, and material properties (e.g., density, refractive index, solubility, and Poisson's ratio). Thus, its exploration requires efficient computational methodologies for polymer science. To address this challenge, this paper introduces PolyGraphPy, an open-source, unified Python framework that integrates atomistic simulations with machine learning for accurate property prediction and property-guided polymer design. The framework automates quantum mechanics calculations using Density Functional Tight Binding (DFTB+) to efficiently construct structured datasets for monomers, homopolymers, and alternating copolymers. For property prediction, PolyGraphPy employs Bayesian Graph Neural Networks (GNNs) utilizing stochastic graph representations to predict target properties-such as static polarizability, while providing robust uncertainty quantification. Furthermore, the platform incorporates two complementary generative models for the de novo design of targeted molecules: a SELFIES-based Generative Pre-trained Transformer (GPT) and a Genetic Algorithm (GA) based on BRICS graph fragmentation. Demonstrated on a dataset of acrylates, PolyGraphPy provides a highly customizable, end-to-end pipeline that reduces computational costs and accelerates data-driven polymer informatics.

We present a comparative study of four supervised machine learning regression algorithms – Support Vector Regression (SVR), Gradient Boosting Regression (GBR), Gaussian Process Regression (GPR), and Multilayer Perceptron (MLP) – for predicting the proton structure function $F_2^p(x, Q^2)$ using high-precision BCDMS experimental data. Unlike conventional methods that solve the DGLAP evolution equations, our data-driven framework directly captures the complex nonlinear dynamics of partonic structure. To ensure statistical robustness, we employ $k$-fold cross-validation and perform thorough hyperparameter optimization. Our results show that the MLP and GPR models achieve superior predictive accuracy. In particular, MLP exhibits the highest sensitivity to nonlinear gradients, while SVR proves most stable against experimental uncertainties. The close convergence of training and validation metrics confirms that the models learn the underlying QCD physics without overfitting to statistical fluctuations. This work highlights the potential of ML-based regression as a complementary tool for structure function analysis and kinematic extrapolation in high-energy physics.

We formulate stationary-density-preserving nonreversible perturbations of Fokker–Planck dynamics as gauge fields that deform relaxation spectra while leaving the invariant state fixed. When detailed balance holds, a similarity transformation maps the reversible Fokker–Planck operator to a Witten-Laplacian-type supersymmetric Hamiltonian; nonreversible gauges then appear as non-Hermitian perturbations that preserve the zero mode but modify the excited spectrum. This operator viewpoint gives a common language for relaxation gaps, circulating probability currents, hypocoercive acceleration, and finite control costs. We represent admissible gauge currents by antisymmetric tensor fields and identify the detailed-balance-violating Ohzeki–Ichiki force as a constant symplectic example whose infinite-strength limit is Hamiltonian dynamics. The continuous-time spectral gap alone does not select a finite gauge strength, so we introduce a finite-time regularized objective and an actor–critic procedure for learning the gauge. An exactly solvable anisotropic Gaussian Ornstein–Uhlenbeck benchmark separates the spectral transition from the finite-time optimum and shows that the learned gauge recovers the Lyapunov-equation optimum. A double-well benchmark then illustrates the same constrained selection in a nonconvex metastable landscape. Stochastic gradient methods enter this framework as physically relevant Fokker–Planck systems: mini-batch noise acts as an effective diffusion tensor, and adaptive methods such as Adam correspond to metric choices with possible nonequilibrium currents.

Ultralight axions (ULAs) are well-motivated dark matter particle candidates that arise in many extensions of the Standard Model of particle physics. ULAs with mass $m_\mathrm{a} \lesssim 10^{-27}$ eV have been strongly constrained by cosmic microwave background (CMB) observations in temperature and polarization. We fit recent measurements of gravitational lensing of the CMB from Planck, the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT-3G) using a state-of-the-art simulation-calibrated nonlinear clustering model for ULAs. We derive the strongest constraints on ULAs in the mass range $10^{-26}\;\mathrm{eV}\leq m_\mathrm{a}\leq 10^{-24.5}\;\mathrm{eV}$. ULAs of this mass have been shown to alleviate tensions between inferences of the matter clustering if they compose a few percent of the total dark matter content of the Universe. We conclude that ULAs with a mass of $10^{-26}$ eV make up less than 1.5% of the dark matter and $10^{-25}$ eV axions make less than 9% (both at 95% confidence level). We identify a slight preference for non-zero axion density at $10^{-24.5}$ eV at $2.1\sigma$. We find that the preference for ULAs is largely driven by a few data points and that further investigation of non-linear ULA physics is needed to confirm or rule out this signal.

While core-collapse supernovae have been widely discussed as potential neutrino sources, definitive observational evidence has remained elusive. In this work, we report evidence of high-energy neutrino emission in the direction of supernova SN 2021foa, which is one of the closest and brightest interacting supernovae observed to date. Using the second data release of muon track data from the IceCube Neutrino Observatory, we conducted a time-dependent analysis and identified a neutrino clustering that temporally coincides with the optical peak of SN 2021foa, occurring approximately 16 to 22 days after the discovery date, with a maximum test statistic of $\sim 28.2$. Monte Carlo simulations indicate that the probability of observing such a neutrino excess by chance is $\sim6.7 \times 10^{-5}$, corresponding to a significance of $\sim4.0\,\sigma$. The spatial and temporal correlation strongly suggests that the neutrinos originate from the supernova. SN 2021foa is a unique "flip-flop" supernova; its spectra repeatedly transitioned between hydrogen-rich (Type IIn) and helium-rich (Type Ibn) phases within 50 days post-peak, reflecting a violent and complex mass-loss history of its progenitor. The inferred neutrino energy exceeds the optical radiative energy and ejecta kinetic energy of the supernova by orders of magnitude, suggesting that the neutrino emission is likely powered by a delayed central engine driving a jet that is choked within the dense circumstellar medium.

This study presents a novel framework for estimating the central sea-level pressure ($P_\mathrm{c}$) of tropical cyclones (TCs) using infrared radiometers. We leverage the long-overlooked combination of high spatial resolution and sounding capability of the Moderate Resolution Imaging Spectroradiometer (MODIS) to measure warm-core anomalies in TC eyes. We develop the MODIS Thermal Infrared Sounding (MOTIS) framework, which performs instrument-specific preprocessing and estimates $P_\mathrm{c}$ using multiple linear regression. MOTIS yields $r^2 = 0.945$ and RMSE = 4.3 hPa for high-intensity TCs with observed clear eyes (mean $P_\mathrm{c} = 937$ hPa), outperforming all existing methods for intense TCs. We construct a dataset of 3288 (1082 clear-eye) MOTIS estimates from 2002 to 2025 and demonstrate its potential to improve the quality of Best Track $P_\mathrm{c}$, roughly halving uncertainties in the absence of pressure observations. Although MODIS is nearing the end of its mission, the MOTIS framework could be extended to next-generation geostationary sounders to provide accurate real-time $P_\mathrm{c}$ estimation for high-intensity TCs.

We propose a resolution of the firewall paradox in JT gravity by incorporating topology change into canonical quantization under relational time evolution. The gravitational Hamiltonian acts on the black hole interior through a pair of pants interaction, mapping between a single interior sector and a connected two interior sector. To describe dynamics in the interior while keeping track of the exterior, we pass to an extended phase space description obtained by splitting the bulk Hilbert space across the event horizon. The split introduces boost edge modes at the horizon, which the Hawking modes become gravitationally dressed to. Covariance of the resulting crossed product algebra provides a precise gravitational realization of the firewall: a one sided boost of the interior edge mode relative to the exterior holding the matter fixed, or equivalently, a relative phase between the interior and exterior Hawking partners holding the edge modes fixed. Although each topology changing transition is exponentially suppressed, evolution over a Page time causes the connected two interior branch to dominate. One of these is the naive semiclassical interior, which we show carries a nontrivial one sided boost upon gluing the interior back to the exterior, and hence a firewall. The other interior is shown to be a zero mode of the one sided boost generator. Gluing the interior back to the exterior quotients by the gravitational constraints, which annihilates the firewall branch. On the surviving branch, we show the horizon vacuum measurement and the early radiation purity measurement become the same Dirac observable. Equivalently, we show that Page time dynamics induces a large diffeomorphism on the connected branch under which the operator algebra of the interior Hawking partner and that of the decoded early radiation are identified.

Lens spaces are useful test examples in spectral geometry because their spin Dirac eigenspaces admit explicit congruence descriptions. We use these descriptions to study equivariant $\eta$ invariants for three-dimensional lens spaces with the round metric and the standard coordinate-torus action, retaining the spin-Fourier character of each eigenspace rather than only the ordinary scalar $\eta$ value. For the square family $L(\ell^2,\ell-1)$ and $L(\ell^2,2\ell-1)$, with $\ell\geq 5$ odd, we obtain a residual-circle equivariant $\eta$ separation: the ordinary $\eta$ values agree, and the first derivative of the residual $\eta$ germ vanishes by symmetry, but the second derivative is nonzero. For $L(25,4)$ versus $L(25,9)$, the normalized second derivative is $-6080$. Thus, the residual-circle equivariant $\eta$ germ detects a distinction invisible to the ordinary $\eta$ invariant. The calculation uses spin-Fourier residues directly; perturbative Hessian signs serve only as motivation and are not part of the invariant.

Let $\mathcal{A}$ be a (not necessarily rational) conformal net with a faithful action of a finite group $G$. Let $\text{Rep}^G(\mathcal{A})$ be the $G$-crossed balanced $\mathrm{W}^*$-tensor category of $G$-twisted representations of $\mathcal{A}$ as introduced in arXiv:2606.03623. We show that there is an equivalence of balanced $\mathrm{W}^*$-tensor categories $(\text{Rep}^G(\mathcal{A}))^G\cong \text{Rep}(\mathcal{A}^G)$ between the $G$-equivariantization of $\text{Rep}^G(\mathcal{A})$ and the category of representations of the fixed-points conformal net $\mathcal{A}^G$. This generalizes to the non-rational case the equivalence of braided tensor categories $(\text{Rep}^G(\mathcal{A}))^G\cong \text{Rep}(\mathcal{A}^G)$ for $\mathcal{A}$ rational appearing (in the language of localized endomorphisms) in arXiv:math/0403322, and it also includes the balances.

We construct a periodic weighted graph whose discrete Laplacian has a spectrum with precisely $n$ gaps. Moreover, we show that by an appropriate choice of the weights, the endpoints of these gaps, as well as the upper edge of the spectrum, attain the prescribed values. The underlying graph has a brush-like geometry: it consists of an infinite chain of vertices, each of which is connected to $n$ additional pendant vertices by extra edges. Semi-explicit formulae for the weight coefficients are provided: some of the coefficients are determined explicitly, while others are given as roots of an explicitly determined polynomial.

We propose a phenomenological mass estimator for an intermediate ultralight axion dark-matter component using thin stellar streams. The central observation is an analytic relation linking three length scales in a fuzzy-dark-matter stream progenitor: the tidal radius of a bound dark core, its gravitational Bohr radius, and the axion de Broglie wavelength evaluated at the stripped-star velocity scale. If the stream width is (w=C_w r_t), with (C_w=O(1)), then [ \frac{r_t}{r_B} =========================== \frac{4\pi^2}{C_w^2} \left(\frac{w}{\lambda_dB}\right)^2 . ] Thus (\lambda_dB\sim w) automatically implies (r_t/r_B\sim25)–(40): the dark wave core survives well inside the tidal boundary while the extended stellar envelope is stripped into a thin stream. This wave–tide locking gives a no-simulation axion-mass estimator, [ m_a \simeq 2.69\times10^{-19},eV \left(\frac{\sqrt{2}}{q_\kappa}\right) \left(\frac{R}{10,kpc}\right) \left(\frac{38,pc}{w}\right)^2 \left(\frac{w}{\ell_ripple}\right) \left(\frac{220,{\rm km,s^{-1}}}{v_c}\right). ] In the simplest locked case (\ell_ripple\sim w), a homogeneous first-pass set of narrow stream widths points to (m_a) of a few (10^{-19},eV) and hidden-core masses of a few (10^3)–(10^4M_\odot). This is not a detection claim; it is a falsifiable phenomenological test for an ultralight axion component that can be sharpened by homogeneous stream catalogs and Schrödinger–Poisson simulations.

We study the robustness of asymptotic entanglement manipulation for Mazzola–Sutter–Renner (MSR) almost-iid sources, which allow a sublinear number of deviations from tensor-power structure. For pure sources along $|{\phi}\rangle_{AB}$, every concentration rate below $S(\phi_A)$ is achievable by a single Schur–Weyl protocol depending only on the rate and local Hilbert spaces. For general sources along $\rho_{AB}$, every distillation rate below $E_D(\rho_{AB})$ is achievable by a single LOCC protocol sequence independent of the particular source and defect size. Its error vanishes exponentially and uniformly for each prescribed defect bound $r_n=o(n)$, yielding an optimal universal distillation rate equal to $E_D(\rho_{AB})$. For dilution, the entanglement cost of every MSR almost-iid target sequence along $\rho_{AB}$ is at most $E_F^\infty(\rho_{AB})$, the regularized entanglement of formation of the reference state. Consequently, every MSR almost-iid sequence along a pure reference state $|{\phi}\rangle_{AB}$ is asymptotically reversible, with distillable entanglement and entanglement cost equal to $S(\phi_A)$, even when its individual states are mixed. The proofs establish structural and entropic properties of MSR almost-iid sources that may be useful in other information-theoretic settings.

Magic-angle twisted bilayer graphene (MATBG) hosts flat bands whose correlated and topological phases depend sensitively on the electronic point-group symmetry, which is set by the device configuration including twist angle, heterostrain, and substrate and cap-layer coupling. Here we report near-infrared photocurrent measurements of a MATBG device encapsulated by an hBN substrate and a WSe$_2$ cap layer, at a twist angle of 1.05{\deg}. Using spatially resolved photocurrent imaging with polarization analysis, we disentangle the photothermoelectric response from the circular photogalvanic effect (CPGE) across all fillings of the moiré flat bands. A finite CPGE at normal incidence, forbidden in high-symmetry structures, directly constrains the point-group symmetry of MATBG and reveals a Berry curvature dipole intrinsic to the MATBG flat bands, activated by the symmetry lowering from heterostrain and the asymmetric hBN/WSe$_2$ environment. Our results establish nonlinear photocurrent as an all-optical probe of symmetry and quantum geometry in proximity-engineered moiré materials.

This work examines the effects of multiple parton radiation on proton-proton $W^\pm$ boson production at the Relativistic Heavy Ion Collider (RHIC), focusing on the ratio of charged lepton pseudorapidity differential cross sections to test the SU(2) flavor symmetry violation in the proton quark sea. I show that the measured ratio is marginally affected by the fiducial constraints imposed on the QCD recoil radiation by the STAR experiment, and that the ratio is also stable with respect to transverse momentum resummation effects. Hence, the STAR data on the $W^+/W^-$ cross section ratio provides a robust discrimination of $\bar{d}(x)$ and $\bar{u}(x)$ parton distribution functions (PDFs) at momentum fractions of order 0.1, shown by comparing predictions with various PDF models and using the $L_2$ sensitivity analysis and reweighting methods.

Whether cold dense QCD matter undergoes a strong first-order phase transition remains an open question. In nature, neutron stars provide the most direct probe of cold dense QCD matter. Theoretically, chiral effective field theory constrains the equation of state of dense matter near nuclear saturation density, while perturbative QCD calculations constrain it at densities well beyond stable neutron-star interiors. We perform Bayesian inference with non-parametric Gaussian-process equation of state for $\beta$-equilibrated neutron-star matter under the assumption with and without a strong first-order phase transition, using the tidal deformability from GW170817, the NICER mass–radius measurements of PSR~J0740$+$6620, PSR~J0030$+$0451, PSR~J0437$-$4715, PSR~J0614$-$3329, chiral effective field theory, and perturbative QCD. Within this Bayesian inference framework, the model comparison mildly to moderately favors a strong first-order phase transition, with its onset most likely lying above the central density of the most massive neutron star. Such an onset reconciles the stiffness required to support massive neutron stars with the softening favored by perturbative QCD from asymptotically high density.

As a binary neutron star inspirals due to the emission of gravitational waves, the rising tidal frequency resonantly excites vibrational modes. These oscillations are seismological probes of the rich stellar interior, yet it remains to be established whether gravitational-wave interferometers can measure them. Here, we present the first fully Bayesian study of the capability of the Einstein Telescope to detect tidal resonances. We simulate one year of observations and analyse the 200 loudest signals. We find that the Einstein Telescope can identify resonant modes and is sensitive to gravitational-wave phase shifts as small as $\Delta \Phi \approx 0.03$ for favourable events. We further show that neglecting resonances can bias the inferred tidal deformabilities. These results establish tidal resonances as a measurable route for asteroseismology with future detectors.

We study the motion of columnar ice crystals that form in clouds over a range of low temperature. Our focus here is on elongated ice crystals, which are smaller than the size of the smallest eddies in the flow, with a moderate aspect ratio comprised between $3$ and $5$. We determine turbulent solutions of the Navier-Stokes equations over a range of turbulent kinetic energy dissipation characteristic of clouds ($4.41\;{\rm cm}^2/{\rm s}^3 \le \varepsilon \le 1120\;{\rm cm}^2/{\rm s}^3$) by using direct numerical simulations, and we follow the motion of crystals using simplified but realistic models for the motion of non-spherical, elongated particles. The influence of the fluid inertia leads to a preferential alignment of the crystals perpendicular to the direction of gravity, the alignment effect being opposed by the turbulent fluctuations. Along with the strong alignment of the crystal axis perpendicular to gravity, we observe only a weak alignment with the vorticity, much weaker than in the absence of gravity. The settling velocity depends only weakly on the orientation of the crystals, but is strongly enhanced when $\varepsilon$ increases, an effect that we attribute to preferential concentration in the flow. As the inertia of the columnar ice crystals considered here is significant, we observe a strong spatial clustering. Finally, we discuss the relevance of the effects identified here on the collision frequency between ice crystals in cloud conditions.

We construct a quantum-statistical analogue of Einstein's fluctuation argument for black-body radiation in the context of causal-diamond geometry. Starting from the phase space of a stretched horizon inside a Minkowski causal diamond, we quantize the Poisson algebra generated by the fields averaged over stretched-horizon time. We then compute the fluctuations of the averaged area density of the transverse two-spheres in a thermal state constructed in analogue with the black-body thermal state. In the null limit, where the stretched horizon approaches the causal-diamond boundary, this yields a thermal fluctuation formula of the boundary area operator that contains two terms, in direct analogue with the black-body radiation. The term quadratic in the expectation value is the “classical” contribution, while the linear term has the Verlinde–Zurek scaling characteristic of independent microscopic constituents. In direct analogue with Einstein's interpretation of black-body energy fluctuations as evidence for light quanta, we interpret the linear area-fluctuation term as a statistical signature of discrete quanta of geometry. This provides bottom-up evidence for quantum area degrees of freedom and supports the embadon picture of null quantum geometry.

This work analyses and evaluates a few realistic compact objects in the presence of a gravitational interaction between two particles with a nonmetricity $Q$. In the $f(Q)$ gravity framework, we have selected the anisotropic equation of motion and have determined $f(Q)$ to be a linear function of nonmetricity $Q$. To evaluate the field equations in our work, we have opted to employ the Krori-Barua metric. We calculated the anisotropic factor for each of the four compact objects and found that the anisotropic component is positive and increases monotonically and interpreted that the nuclear force can oppose the gravitational attraction. At last, the relationship between mass and radius has been determined and illustrated visually. We have noted that the compactness of the pulsars LMC X-4, SMC X-4, Cen X-3, and Vela X-1 is inside the Buchdahl's limit for varying values of $a$. This has led to the interpretation that these pulsars are neutron stars in a modified gravity background of $f(Q)$. In addition, we calculated the model mass and, using thirty distinct choices of $a$, ran the Chi-Square test to see if there was a noticeable difference between the observed and model-generated masses. We have also looked at how the surface redshift has changed over time and whether the compact objects in our model that were previously described are compact.

Fault-tolerant quantum computation is inherently a spacetime problem, requiring not merely good static quantum error-correcting codes but also low-overhead protocols for protecting and manipulating encoded quantum information over time. Fault complexes provide a homological formalism for treating such protocols as single spacetime objects. Here we initiate the study of low-overhead fault complexes by introducing spacetime lifting, a method that constructs fault complexes from symmetry-reduced product structures beyond standard foliation. We show that spacetime lifting yields fault complexes and in particular memory experiments with almost-linear fault distance in the total spacetime cost, which substantially outperforms existing constructions. We further endow fault complexes with the operational interpretation of measurement-based cluster state protocols and identify general conditions under which they realize fault-tolerant logical teleportation, showing that spacetime-lifted constructions combine favorable parameter scaling with operational realizations. Our work establishes a systematic route towards more efficient quantum fault tolerance through general complex constructions.

We study thermalisation and the possible occurrence of the Mpemba effect in a closed quantum setting that mimics the interaction of a system with thermal reservoirs coupled only at its boundaries. Specifically, we consider a tripartite geometry in which a finite chain, initially prepared at a finite temperature, is suddenly connected on both sides to two semi-infinite chains of the same nature held at a different temperature. These outer chains act as thermal baths, while the full system evolves unitarily under a homogeneous Hamiltonian. This setup provides a simple quantum realisation of a temperature quench and closely resembles the original scenario in which the classical Mpemba effect was first observed. We focus on two paradigmatic free-fermion models, the XX chain and the transverse-field Ising chain, which respectively preserve and break the global $U(1)$ particle-number symmetry. As a probe of relaxation, we consider the Frobenius distance between the time-evolved reduced density matrix of the central subsystem and its stationary state, which is the thermal state at the bath temperature. Exploiting the free-fermionic structure of both models, the dynamics remains Gaussian and the Frobenius distance can be expressed exactly in terms of two-point correlation functions. Combining this representation with generalised hydrodynamics, we derive analytical predictions for the Frobenius distance in the hydrodynamic limit, providing a complete characterisation of the thermalisation process. Using these results, we investigate the possible occurrence of the Mpemba effect. We find that, despite the genuine non-equilibrium dynamics displayed by the system, no Mpemba effect arises in this setting. Our analysis identifies a broad class of boundary-driven thermalisation protocols in which relaxation is fully characterised analytically and exhibits no anomalous acceleration of equilibration.

We present $\tt BlackHawk$ $\tt v3.0$, a major update of the public code designed to compute Hawking radiation spectra of black holes. Building upon previous versions, this release considerably extends the range of black hole geometries that can be studied by implementing several new spherically symmetric metrics: the Bardeen and Hayward regular black holes, the Simpson-Visser and Peltola-Kunstatter black-bounces, the D'Ambrosio-Rovelli black hole-to-white hole metric, and the Babichev-Charmousis-Lehébel black hole. For each metric, we compute the corresponding Hawking temperatures and greybody factors, enabling the determination of primary Hawking emission spectra for particles of different spins. The greybody factors are obtained through dedicated numerical routines based on the companion code $\tt GrayHawk$. Additionally, $\tt BlackHawk$ $\tt v3.0$ introduces several technical improvements aimed at enhancing precision and efficiency, providing a highly versatile tool. The code is publicly available at https://blackhawk.hepforge.org/

We construct an explicit equivariant quantum Satake map for Grassmannians, which enables us to express their torus-equivariant quantum cohomology in terms of that of projective space. We then consider the exterior algebra of the latter, which admits a canonical identification with a Clifford algebra. We describe the resulting action in several complementary ways: first, from a geometric perspective via push-pull maps, and second, in terms of the shuffle product, which also arises in the simplest cohomological Hall algebra associated with the $A_1$-quiver. Exploiting the Clifford algebra structure, we derive new recurrence relations among equivariant Gromov-Witten invariants, yielding a new method for their computation in terms of Wick's Theorem. As an application, we provide combinatorial proofs of Graham positivity for both equivariant quantum Pieri rules, and in one case extend these results to quantum triple Schubert calculus.

We study deconfined quantum critical points (DQCP) associated with Lieb-Schultz-Mattis (LSM) anomalies in one-dimensional spin chains. Our starting point is a structural characterization of the LSM anomaly in the Lyndon-Hochschild-Serre spectral sequence: $\omega_{\mathrm{LSM}}\in E_\infty^{1,2}= H^1(\mathbb Z_{\mathrm{trans}},H^2(G_{\mathrm{int}},\mathrm{U}(1)))\subseteq H^3(G_{\mathrm{int}}\rtimes_{\rho}\mathbb Z_{\mathrm{trans}},\mathrm{U}(1))$. Physically, this class decorates a translation defect with a projective representation of the internal symmetry $G_\mathrm{int}$. We show that gauging the internal symmetry in the presence of an $E_\infty^{1,2}$-type anomaly necessarily produces a non-invertible dual symmetry. This gives a general mechanism for type-II DQCP: in contrast to type-I examples with $E_\infty^{2,1}$-type anomalies which are dual to ordinary group-like symmetry breaking, type-II transitions are dual to spontaneous breaking of a non-invertible symmetry. We illustrate the mechanism using a spin-$1/2$ chain with an anomalous $D_8$ LSM symmetry. We construct a dimer-to-ferromagnet DQCP candidate, provide numerical evidence for a critical theory with central charge $c\approx 1$, and show, using both category theory and explicit lattice constructions, that gauging the internal symmetry yields the non-invertible $\mathrm{Rep}(H_8)$ dual symmetry.

The dynamics of open quantum systems is commonly described by quantum master equations derived under the assumption of weak system-bath coupling and a separation of timescales between system and bath. When the system is additionally subjected to a periodic driving, the validity of the resulting Floquet master equations is further restricted to regimes of weak or high-frequency driving. Here, we benchmark a set of commonly used Floquet master equations for a model of two locally driven spins coupled to a shared Ohmic reservoir at finite temperature. We systematically probe the accuracy of the equations as a function of the driving parameters, thus identifying limits of their applicability. Dynamical maps predicted by each master equation are compared against numerically exact non-Markovian simulations, tracking the full relaxation dynamics. We find that the accuracy of each master equation closely reflects the assumptions underlying its derivation. For the Floquet-Lindblad equation, errors can be strongly amplified near resonances where the secular approximation breaks down, while approaches that avoid the secular approximation perform better and exhibit a more systematic dependence of the error on driving frequency and amplitude.

Alfvénic turbulence is vital to powering the solar wind and corona, yet eludes a comprehensive understanding of the kinetic processes by which it dissipates. Minor ions are sensitive tracers of these processes, showing extreme perpendicular temperatures and mass-weighted temperature trends that can either correlate or anticorrelate with mass-to-charge ratio, $A_i/Z_i$. We use a combination of quasilinear theory and 3D hybrid-kinetic simulations to explain these features and their correlations with properties of turbulence in the fast solar wind. When Alfvénic turbulence is imbalanced, its cascade to ion-Larmor scales is throttled by the helicity barrier. This barrier ultimately leads to high-frequency proton-cyclotron waves (PCWs), both oblique and parallel, the latter of which produce very flat electric-energy spectra ($\mathcal{E}_{E_{\perp}}\sim k_\parallel^{-\eta}$ with $\eta<2$) over the range of scales that are cyclotron resonant with minor ions. While steeper spectra lead to a positive correlation of heating with $A_i/Z_i$, the shallower spectra cause the dependence to invert, with $Q_i\propto Q_{\mathrm{p}}A_i(A_i/Z_i)^{\eta-2}$. Six simulations of balanced and imbalanced turbulence spanning $\beta_{\rm p0}=\{1,0.3,1/16\}$ corroborate this prediction, showing minor-ion heating rates that follow $(A_i/Z_i)^a$. Minor-ion heating is strongest and most perpendicular in our lowest $\beta_{\rm p0}=1/16$ simulation of imbalanced turbulence, reaching $T_{\perp{\rm O}^{5+}}/T_{\perp{\rm p}}\approx40$ and $T_{\perp{\rm O}^{5+}}/T_{\parallel{\rm O}^{5+}}\approx10$, consistent with low-coronal observations. Future minor-ion measurements should test whether intervals in which minor-ion thermal speeds decrease with increasing mass-to-charge ratio are associated with a history of large cross helicity, enhanced power in parallel PCWs, and a steep transition-range spectrum.

We enlarged the capabilities of the publicly available Mathematica-based tool $\tt GrayHawk$. This second version expands the spectrum of metrics that can be considered in two distinct and disjoined directions. First, it enables a fully numerical computation of the tortoise coordinates integral, allowing the user to account for many metrics for which an analytic computation was impractical. Second, it extends the scattering problem to wormhole solutions. The pool of pre-loaded metrics is enriched, enabling immediate testing of the new features, and the code's modular structure is maintained to facilitate easy modification. This implementation proves $\tt GrayHawk$ adaptability and makes it an even more powerful tool for studying black holes, wormholes, Hawking radiation, and other features involving field propagation on curved manifolds. The codes described in this work are publicly available at \href{https://github.com/marcocalza89/GrayHawk-v2}{\faGithub}.