Papers

Singlet-triplet qubits offer an attractive encoding for semiconductor quantum computing, combining ancilla-free readout, reduced sensitivity to common-mode noise, and baseband voltage control. However, the Zeeman energy difference $\Delta E_\mathrm{Z}$ is typically fixed by local magnetic field gradients or $g$-factor inhomogeneities, leaving the exchange interaction $J$ as the only dynamically tunable parameter. This always-on $\Delta E_\mathrm{Z}$ precludes orthogonal control of the qubit's rotation axes and introduces unwanted state rotations during idling. Here we demonstrate all-electrical orthogonal control of a degenerate singlet-triplet (DST) qubit formed by two hole spins in a germanium double quantum dot. Exploiting the electrically tunable anisotropic $g$-factors of the two spins, we identify a regime where both $\Delta E_\mathrm{Z}$ and $J$ vanish, making the $S$ and $T_0$ states degenerate at the idle point. By applying only baseband voltage pulses, we independently control both $J$ and $\Delta E_\mathrm{Z}$, enabling fully orthogonal $Z$- and $X$-axis rotations. Randomized benchmarking yields an average physical single-qubit gate fidelity of 99.53% for a gate duration of approximately 100 ns. Finally, we electrically tune the degenerate point across a wide range of magnetic field orientations, enabling operation in a regime of enhanced coherence time and offering a route towards multi-qubit scaling under a shared global magnetic field.

We investigate a magnetophotonic crystal formed by the pairs of layers with the opposite out-of-plane magnetization. Despite vanishing total magnetization such antiferromagnetic pattern opens photonic bandgap, gives rise to nontrivial topological phases and enables strong cross-polarized light reflection at frequencies inside the gap which can be harnessed for magnetically tunable polarization-rotating mirrors. We systematically explore optical properties of this structure highlighting fruitful connections to topological photonics and axion electrodynamics.

We develop a structure-preserving neural variational framework for random Young-diagram ensembles, with representations adapted to the structure and scaling of each measure. The method is validated on the Plancherel, uniform, minimal-difference, and fixed-\(q\) \(q\)-Plancherel ensembles, using known asymptotic profiles only for post-training comparison. We then study a quartically deformed hook-length ensemble without assuming an analytical saddle shape. Large-\(n\) neural profiles are compared with finite-size MAP profiles obtained from exact-action searches and with mean profiles obtained from corner-transfer Metropolis–Hastings sampling. Increasing the deformation suppresses the leading rows and broadens the support, while the neural, discrete, and sampled mean profiles agree at the percent level. These results provide numerical evidence for a deformation-dependent macroscopic saddle family.

Simulating an open quantum system on a digital quantum computer often involves the use of Trotter-Suzuki (TS) Product Formulas (PF) to approximate the system's time evolution. Precise estimates for the required number of Trotter steps (and hence the overall gate count) can be crucial for minimising the computational cost of these methods. Building on established theoretical guarantees, we derive analytic bounds for the First- and Second-Order Deterministic and Randomised TS-PF, directly relating the number of Trotter steps to the model parameters, evolution time and precision. These bounds enable concrete resource estimation for each method. We then present a computationally efficient classical algorithm that uses diamond norm estimates of individual Liouvillian terms and a binary search to significantly reduce the Trotter steps required for a target precision. Our numerical results on two prototypical models - an XX-Spin Chain with boundary driving and local dephasing, and a Transverse-Field Ising Model - show that the theoretical (analytic) bounds are often overly conservative, whereas the empirical (optimised) bounds yield a significantly smaller number of Trotter steps for the same precision. Among the methods investigated, the Second-Order Randomised TS-PF typically achieves the lowest resource demands, especially for larger systems. These findings emphasise the significance of empirical bounding strategies in achieving more resource-efficient simulations of Markovian open quantum systems.

The Teukolsky equation possesses discrete symmetries that constrain the properties of black hole perturbations and their quasinormal mode spectrum. In this study, we explore how a class of modifications of the Teukolsky potential can alter the symmetry structure of the equation and break the m = 0 degeneracy of quasinormal modes. We prove this result in frequency domain using the master Teukolsky equation and in time domain via (2+1)-dimensional simulations. We also show that the discrete symmetries can be leveraged for a more efficient characterization of the m = 0 quasinormal modes from time-domain evolutions. As a theory-specific application, we consider the case of higher-derivative theories of gravity, highlighting that time-domain implementations of frequency-domain potentials can give rise to additional non-physical branches of modes.

The outskirts of nearby galaxies preserve long-lived signatures of hierarchical assembly and provide a fossil record of galaxy formation. J-ATLAS, the Javalambre Andromeda and Triangulum Legacy Astrophysical Survey, is a proposed wide-field optical photospectroscopic survey of the M31–M33 system with JST/T250 and JPCam. Using the dense J-PAS filter system, it will deliver homogeneous low-resolution optical SEDs for resolved stars, compact stellar systems, emission-line sources, and background objects across the nearest massive spiral-galaxy environment visible from the northern sky. The survey will map disks, halos, streams, satellites, and the intergalactic region, constraining stellar populations, photometric metallicities, carbon-rich evolved stars, star clusters, planetary nebulae, and low-surface-brightness substructure. We also present the ongoing pilot survey, centred on the Giant Southern Stream and Andromeda~I, whose 92 complete pointings are being analysed to validate crowded-field photometry, calibration, source classification, and photometric metallicity diagnostics before extending the methodology to the full J-ATLAS footprint survey.

We use Diagonal Adaptive Non-local Observables (DANO) as a canonical decomposition for studying Variational Quantum Circuit model evolution. Separating each learned observable into a diagonal spectrum and a unitary basis gives a quasi-polar description: the spectral weights are viewed as radial coordinates, while the unitary circuit serves as angular coordinates through Lie group identifications. This turns the training process into a trajectory in spectral and Lie-algebra space. Experiments on two classification tasks show that DANO radial spectral expansion correlates with accuracy. DANO angle coordinates reveal a dominant accuracy-correlated component. The framework provides a different perspective to characterize quantum model behavior.

Next generation X-ray observatories require fast, low-noise, low-power, multi-megapixel imaging spectrometers. To meet these demands, the X-ray Astronomy and Observational Cosmology (XOC) Group at Stanford, in partnership with the MIT Kavli Institute and MIT Lincoln Laboratory (MIT-LL), is developing multi-channel X-ray charge-coupled devices (CCDs) and fast readout architectures. We report the energy resolution and noise performance achieved with a full-scale (1440x1440-pixel), 16-channel, front-illuminated MIT-LL CCD detector developed for the Advanced X-ray Imaging Satellite (AXIS) concept, the CCID-100, read out using two Multi-Channel Readout Chip (MCRC) V1 application-specific integrated circuit (ASIC) chips in the new Stanford CCID-100 test setup. We describe an automated method for bias optimization on each CCD channel, and integrated debugging features of the front-end ASIC and readout system. The demonstrated performance confirms that these systems can meet the speed and noise requirements of future strategic X-ray missions.

We present a timing-based framework for interpreting braking-index measurements in young pulsars using a four-channel spin-down model that includes particle-wind, magnetic-dipole, mass-quadrupole, and current-quadrupole torques. The observed braking index is a torque-weighted average of the channel exponents, enabling equation-of-state-independent constraints on the torque fraction of a possible r-mode-like current-quadrupole component from timing data alone. For positive, slowly varying secular torques, the allowed range is (1\le n_obs\le 7), with (f_{7,\min}=\max[0,(n_obs-5)/2]) and (f_{7,\max}^phys=\min{1,\max[0,(n_obs-1)/6]}). These bounds rely on non-negative, slowly evolving torque coefficients; if these assumptions fail, a large braking index need not uniquely indicate a current-quadrupole torque. Applied to young pulsars with measured braking indices, the analysis shows that most sources do not require gravitational-wave spin-down and are consistent with wind-plus-dipole braking. Sources with (3<n_obs<5) require an additional higher-order contribution, but timing alone cannot identify it uniquely as an r-mode torque. PSR~J0537$-$6910 is the most suggestive case: its large inter-glitch braking index approaches the (n\simeq7) limit and is consistent, within the restricted secular model, with a strong current-quadrupole-like contribution. However, vortex-creep and superfluid-recovery effects may produce similar inter-glitch behaviour without gravitational-wave emission. We also derive stellar-model-dependent r-mode amplitude bounds, braking-index-corrected ages, and continuous-wave ranking metrics for current and future detectors, including the reduced sensitivity expected for glitch-limited semi-coherent searches.

Quantum interference is conventionally associated with the coherent superposition of alternative paths. We uncover a temporal form of interference in repeatedly monitored quantum systems, appearing directly as fringes in the stochastic first-detection-time distribution. This temporal interference becomes remarkably long lived near a class of topological transitions at which a pair of dark states emerges in Hilbert space. Close to such transitions, the corresponding quasi-dark modes dominate the late-time dynamics and interfere coherently, producing pronounced temporal fringes whose decay becomes arbitrarily slow upon approaching the transition. We identify the symmetry conditions underlying this mechanism and relate the fringe frequency to the energy spectrum of the Hamiltonian. Our results reveal how the interplay of quantum superposition, topological transitions, and dark-state formation generates long-lived temporal interference in monitored quantum dynamics.

Kerr black-hole quasinormal modes are usually computed with numerical methods, but analytic approximations remain useful for identifying the physics that controls different parts of the spectrum. In this paper, we ask whether the divergent, high-order Wentzel-Kramers-Brillouin (WKB) expansion about the peak of the Chandrasekhar-Detweiler potential can be made predictive through Padé and Borel-Padé resummation. We develop two complementary implementations: a semi-analytic slow-rotation expansion in the dimensionless spin $a$ (carried out through 21st WKB order), and a fixed-spin Padé-WKB implementation for the resummed frequency equation (carried out through 41st WKB order). In the slow-rotation regime, the 21st-order resummed expansion is significantly more accurate than the fourth-order approximation found previously. For damped modes at larger spins, the fixed-point iteration agrees well with Leaver's method, reaching fractional errors below $10^{-7}$ in the real part of the fundamental $m=0$ mode at $a=0.99$. The same strategy fails for modes that approach the zero-damped branch near extremality. We trace this breakdown to the near-horizon structure of the Chandrasekhar-Detweiler potential. As the extremal limit is approached, nearby poles produce rapid variation on the throat scale, so a local Taylor expansion about the potential peak no longer uniformly captures the relevant region.

Understanding the boundaries between quantum thermalization and localization in many-body systems remains a central frontier of condensed matter and quantum information science. In this work, we investigate the dynamics and spectral properties of a generic model with long-range three-body-interaction, namely, a system with non-local three-wave-mixing. This model has been realized recently with a microwave Fabry-Perot cavity terminated on one end by a superconducting qubit mirror. Utilizing exact diagonalization techniques, we uncover a striking paradox: the global energy level spacing statistics show integrability, even though all dynamic observables and inverse participation ratios of the eigenstates indicate ergodicity and delocalization. We show that this behavior is a hallmark of strong Hilbert space fragmentation driven by kinematic constraints rather than an explicit global symmetry. Inside these sectors, dynamics scramble rapidly, as evidenced by the out-of-time-ordered correlator (OTOC), while global transport is heavily bottlenecked, resulting in a logarithmic relaxation to equilibrium. This picture is further confirmed by fluctuations in eigenstate entanglement entropy at the same energy. Finally, we demonstrate that the late time OTOC average scales with system size, providing a distinct experimentally accessible signature of the underlying three-body kinetic bottlenecks.

We propose and analyze a novel quantum random-number generator based on a tripartite quantum system in which two subsystems act as detectors. Within a quantum non-demolition measurement scheme, one detector is used to certify the presence of genuine quantum effects in the system's evolution, while the second generates random numbers from a distribution that can be optimized to maximize their entropy. Using one two-level system and two three-level systems, we generate random numbers from a nearly uniform three-outcome distribution, yielding close-to-maximal entropy and therefore near-optimal randomness generation. A key feature of the protocol is that randomness generation and certification occur simultaneously. Moreover, certification does not rely on spacelike separation between detectors, removing a major constraint of device-independent approaches. This property enables practical implementation and facilitates the miniaturization of the device, making the protocol a promising candidate for scalable quantum technologies.

Controlling nuclear fusion is so challenging that for decades, people have been asking: Are we closer to infinite clean energy? Upon the philosophy of beauty and simplicity, the current work points out that there can be a shortcut: the rotating mirror with detached electrodes and all-directional pinch compressions. This is based on the provisional patents filed recently by the University of Texas at Austin. The device combines the steady-state and fast processes in the two main streams of controlled nuclear fusion research: magnetic confinement fusion and inertial confinement fusion. The fuel plasma is preheated in a steady-state process in a rotating mirror with detached electrodes and then the pinch compressions in both radial and longitudinal directions are applied as the fast process. Preheating and longitudinal compression, in addition to the radial compression, significantly boost the nuclear fusion rate. Fast compression after the preheating limits the time for thermalization between ions and electrons and, therefore, minimizes the impact of electron bremsstrahlung radiation loss on ions. Based on the existing experimental results, the current method can be extrapolated to have the potential to reach or exceed the Lawson criterion for the first demonstration of the feasibility of peaceful usage of nuclear fusion energy.

The magnetic response of magnetite as it cycles through the Verwey transition at $T_{\rm V}\approx 120$ K remains incompletely understood, despite numerous studies of both structural and magnetic domains. For example, although the alternating-current susceptibility increases upon heating through the Verwey temperature $T_{\rm V}$, the magnetization may either increase or decrease. Since many intertwined processes occur at $T_{\rm V}$, a multifaceted view of the magnetic changes may help clarify the nature of these processes and establish specific features of the transition as diagnostic tools. To this end, we report measurements of the magnetic response of magnetite upon heating through $T_{\rm V}$ under a wide range of applied-field conditions. We studied synthetic stoichiometric and doped (Zn, Ti, and Al) single crystals, as well as a natural polycrystalline sample, with the aim of explaining how and why the magnetization changes as a function of temperature through $T_{\rm V}$. The magnetic-field-induced reorientation of the monoclinic $c$ axis may broaden the transition, primarily in stoichiometric magnetite. These findings help clarify the mechanism of the transition.

We present a LiDAR-based pipeline for the segmentation and structural characterisation of tall urban vegetation in Bologna. Starting from airborne LiDAR point clouds previously classified as high vegetation by a Random Forest model, we implement and compare two individual-tree segmentation strategies: a watershed algorithm applied to the Canopy Height Model and a point-wise region growing algorithm operating directly on the three-dimensional cloud. Over an area of 12 tiles covering the Talea district, the watershed method detects a larger number of trees (7589), dominated by short trees with a height peak around 5 meters, whereas region growing yields fewer trees (6432) but retains the highest returns, yielding individuals above 40 meters, reflecting the absence of a smoothing step in the latter. We further confront the segmentation results with the municipal Open Data catalogue Alberi in manutenzione: on a sample tile roughly half of the LiDAR-detected trees turn out to be absent from the catalogue, several catalogued positions correspond to locations where no tree is observed, and the majority of the height records date back about two decades, which makes the dataset unsuitable as a ground truth reference. From the segmented trees we extract height and crown radius and use them to compute preliminary vegetation indicators, namely above-ground biomass and carbon storage through allometric relations and annual pollen production through species-specific inflorescence coefficients. All results are collected in an interactive per-tree metadata map. The pipeline is modular and can be recalibrated as improved classifications, field campaigns, or complementary sensing technologies become available, providing a scalable basis for data-driven urban greenery planning.

Spin-dependent transport in ferromagnets underpins the development of high-density spintronic memories. Spin-dependent transport in strong spin-orbit-coupled ferromagnets exhibits a significant anisotropy. Both the overall spin polarization during charge transport and the magnetic spin Hall conductivity are found to exhibit pronounced anisotropy when the magnetization is tilted away from the crystallographic easy axis or when the electric field is rotated relative to the crystal axes. These anisotropic responses originate primarily from spin-orbit coupling, which is identified as the key driver of the large anisotropy observed in ferromagnet. Furthermore, strain tunability of the magnetic spin Hall anisotropy is demonstrated, with tensile strain progressively enhancing the oscillatory amplitude of the spin Hall conductivity. These findings establish strong spin-orbit-coupled ferromagnets as a platform for anisotropic spin-current generation and field-free spintronic devices that exploit intrinsic material anisotropy for improved performance and energy efficiency.

In this paper, we introduce a novel mixed-order Rényi divergence and investigate its fundamental properties. Using this divergence, we define a family of mixed-order order-two Rényi mutual information and Rényi conditional entropy. We derive exact reliability functions of quantum soft covering and privacy amplification under the sandwiched Rényi divergence with order $\alpha\in[2,\infty)$. The former is jointly characterized by the sandwiched and mixed-order order-two Rényi mutual information quantities, while the latter is characterized by the corresponding conditional entropies. These results provide operational interpretations of the proposed mixed-order Rényi divergence. To the best of our knowledge, this is the first exact characterization of the reliability function for quantum soft covering.

We calculate the probabilities $p_{\{n\}}$ of the measured bitstrings $\{n\}$ for the vacuum of Rydberg ladders with $N_q$ atoms. As $N_q$ increases, the $p_{\{n\}}$ decrease but become more dense in the low $p$ region raising the possibility that their smallness could be compensated by their large number. The importance of the low probability states can be estimated from the cumulative probability distribution $\Sigma(p_{\Lambda},N_q)$, which is the probability to observe any state having a probability $p\leq p_{\Lambda}$. For not too large values of $p_{\Lambda}$, it is possible to approximately collapse the $\Sigma(p_{\Lambda},N_q)$ for successive $N_q$ into a function resembling the Fermi function when plotted as a function of $-\ln(p_{\Lambda})$. We show that the number of shots necessary to reduce $\Sigma(p_{\Lambda},N_q)$ to some low enough value grows exponentially with $N_q$. We discuss the implications for calculating observables associated with the vacuum.

Predictive simulations of catalytic interfaces require correlated electronic-structure treatments that describe localized chemical transformations while retaining the influence of the extended metallic environment. We introduce QC-DFET, a quantum-computing density-functional embedding framework that maps surface-reaction active spaces to compact, environment-aware qubit Hamiltonians. A reaction-consistent active-space protocol preserves orbital continuity along reaction coordinates, while quantum-selected configuration interaction based on measurements from the Zuchongzhi superconducting quantum processor and strongly contracted perturbation theory capture static and dynamic correlation. On Cu(111), QC-DFET treats active spaces up to 28 qubits and is validated through a hierarchy of experimentally constrained surface-chemistry challenges. H2 dissociation/desorption tests balanced bond breaking and recombination barriers, CO adsorption tests site selectivity and metal-adsorbate bonding, and formate hydrogenation tests competing hydrogenation branches with different kinetic and thermodynamic signatures. Across these cases, QC-DFET reproduces bidirectional H2 barriers, recovers the observed top-site preference and adsorption strength of CO, and reconciles the experimentally benchmarked H2COO* reverse barrier with the lower forward barrier to HCOOH*. These results establish embedded quantum computing as a practical route to correlated surface-reaction energetics.

On spin manifolds, we give a proof of Brendle and Wang's refined positive mass theorem using the Dirac operator method. In the course of this proof, we discuss the origin of the special coefficient appearing in Brendle and Wang's spectral positivity condition for scalar curvature from the perspective of spin geometry.

In this paper we use multilevel modeling of 10 introductory courses from four STEM disciplines to quantify associations of Learning Assistants (LAs) with the student-level outcomes, pass-rates, and retention. Overall we find that having an LA in a class is associated with improvement of these metrics for all students. Fur- thermore, along with all students showing improvement, we find noticeable demographic differences in the amount of improvement. For instance, LAs seem to be associated with larger improvements in pass-rates for historically marginalized groups, and larger increases in graduation rates for women identifying as belonging to historically marginalized races and ethnicities when compared to students from the same group but with class sections without LAs. This indicates that Learning Assistant programs may be an excellent investment for diverse institutions.

Universal relations provide a particularly useful way to extract physical information from neutron star observables in the presence of various uncertainties by reducing the dependence on uncertain model parameters and microphysical inputs. In this study, we examine the oscillation frequencies of mirror dark matter admixed neutron stars using a two-fluid description, where the outer and inner fluids give rise to two distinct fundamental frequencies. We confirm that the universal relation between the mass-scaled fundamental frequency of the outer-fluid-led mode and the stellar compactness, established previously for self-interacting dark matter admixed neutron stars, also holds in the mirror dark matter scenario. However, this universal relation becomes less robust when metric perturbations are included, compared with the corresponding results in the Cowling approximation. We further find that the inner-fluid-led fundamental frequency can also be expressed as a compactness-dependent relation that is largely independent of the normal matter equation of state, provided that the dark matter mass fraction is fixed. These results suggest that the simultaneous detection of the two fundamental frequencies could provide a way to constrain the dark matter mass fraction, even when the equation of state of normal matter remains uncertain. Finally, we find that the Cowling approximation estimates the fundamental frequency associated with the outer fluid with an accuracy comparable to that found for standard neutron stars without dark matter, while it performs even better for the frequency associated with the inner fluid.

Under the collisionless cold dark matter (CDM) paradigm, galaxies with stellar masses below 10**(5-6) Msun are expected to preserve primordial cuspy dark matter (DM) profiles. Because baryonic feedback should be too weak to transform cusps into cores at these masses, such galaxies provide especially sensitive tests of DM physics. If cores are observed in these systems, they could indicate departures from CDM. To address this problem, Sanchez Almeida et al. (2025) introduced the Eddington Inversion Method tool (EIM-tool), a photometry-based diagnostic that avoids the need for long-term spectroscopic observations and is suitable for the analysis of the forthcoming large photometric surveys of galaxies. The method relies on the fact that embedding a stellar distribution with a central core inside a cuspy NFW halo can require a negative phase-space distribution function f, which is physically impossible. By fitting observed stellar surface densities using f as a free parameter, EIM-tool tests whether an assumed gravitational potential is physically consistent with an observed stellar distribution. This paper validates EIM-tool using globular clusters, dwarf spheroidal galaxies (dSphs), and numerical simulations with independently known potentials. Applied to 21 globular clusters, the method rejected NFW cusps in 71% of cases while consistently favoring cored Schuster-Plummer models, accurately recovering expected core radii. Applied to dSphs, EIM-tool reproduced conclusions from classical dynamical studies: Sculptor and Fornax favored cored halos, whereas Draco remained compatible with a cusp. Tests on EDGE and FIRE simulations further showed that the tool reliably rejects cuspy NFW profiles when the true DM distribution is core-like. Overall, these results establish EIM-tool as a robust and efficient photometric method for distinguishing cored and cuspy DM halos.

Optimization in non-convex neural network models is strongly influenced by the geometry of the solution space: sparse, isolated, point-like clusters are typically algorithmically inaccessible, whereas wide and flat regions can be found efficiently despite being relatively rare. At zero temperature this picture has been formalized in binary perceptrons through the overlap gap property (OGP), which limits algorithmic access to configurations with zero training error above a critical constraint density $\alpha_{\rm OGP}$. Here we extend this description to finite temperature, where a positive training error is allowed and statistically penalized. We first show that the frozen one-step replica-symmetry-breaking solution, dominating the zero temperature equilibrium measure, survives at any finite temperature. We furthermore derive a general criterion, based on the smoothness of the single-pattern Gibbs weight near the decision boundary, that determines when a finite-temperature relaxation of the loss removes freezing. We then extend the OGP construction to finite temperature and show that dense, algorithmically accessible regions of finite-energy configurations persist beyond $\alpha_{\rm OGP}$, up to a threshold $\alpha_{\rm OGP}(\epsilon)$ that grows with the allowed training error $\epsilon$. Finally, in the teacher-student setting, we show that these wide, finite-energy regions still retain good generalization. Using a finite energy message-passing algorithm, we demonstrate numerically that thermal noise enables effective generalization in the regime of constraint densities where both recovering the teacher and finding a zero temperature solution are computationally hard.

Conditioning on both past and future assigns intermediate-time properties a causal observer does not; these time-symmetric assignments obey exact symmetry theorems and are measurable from trajectories. We use two-state generator extended dynamic mode decomposition (gEDMD): because weak values obey $dA_w/dt=i\langle[H,A]\rangle_w$ exactly, generator extraction, with an exact-derivative baseline, applies unchanged to them. First, a reflection involution on the pre-/post-selected ensemble splits every window-fitted friction uniquely as $\gamma_{fwd}=\gamma_A+\gamma_S$: $\gamma_A$, antisymmetric about the midpoint, carries the modes' boundary-condition physics; $\gamma_S$, symmetric, comes from the differencing scheme; both follow from the same data as $(\gamma_{fwd}\pm\gamma_{bwd})/2$. At a fixed inference resolution the arrow of time has two layers: the coherent-mode arrow reverses at the midpoint, the fluctuation-level one does not, $\gamma_S$ dominating $\gamma_A$ at every size and class. The difference is one of degree: $\gamma_S$ is 34 times larger there than at the mode layer, and with the exact derivative both layers reverse: immunity belongs to the inference, not the ensemble. Second, in a lattice interferometer conditioned only at its ports, the quantum Cheshire-cat structure emerges unimposed: particle and polarization obey separate continuity equations, and a local field in the polarization-carrying arm rotates that phase alone, at exactly twice the field strength, entering the generator as a rigid imaginary shift, while the particle's weak density stays invariant to machine precision. We verify the sample-level identity and the two layers from $2^8$ to $2^{20}$ dimensions: $|\gamma_A|/\gamma_S=0.09$ to $0.27$ across five classes; self-averaging makes it insensitive to class among those sharing a boundary modulation, removing the $2^{-N/2}$ overlap obstruction for $N$ qubits.

We consider the necessary and sufficient conditions for the smooth matching of an expanding McVittie spacetime and a spatially flat FLRW spacetime across a general hypersurface: a hypersurface of arbitrary causal character, where the character possibly changes from point to point. We identify and completely characterize a number of special cases. In the general case, we show that the full set of matching conditions is equivalent to a 3-dimensional non-linear system of first order ODEs, subject to a single initial value constraint. Taking the perspective that the McVittie region is specified, we prove a semi-global existence result for this system and interpret these solutions in terms of the physical set-up. Solutions exist and are unique globally to the future. In this direction, the required FLRW spacetime exists and is unique and the matching hypersurface exists almost everywhere and is unique. However, the necessary embedding conditions for the hypersurface necessarily fail at least once - but only at isolated points - along the putative hypersurface. To the past, the matching hypersurface extends to meet the past singularity of the McVittie region. We determine the causal character of the matching hypersurface in these limits, and consider the implications for the existence (or otherwise) of an isotropic source of the McVittie cosmological black hole spacetime. We find that generically, the matching hypersurface has spacelike as well as timelike portions, and so the isotropic region cannot be considered to be an interior for the McVittie exterior. This is in distinction to the Oppenheimer-Snyder model, in which an isotropic dust sphere matches across an everywhere-timelike boundary to the Schwarzschild exterior. On the basis of the sometimes-spacelike nature of the hypersurface, we conclude that an isotropic source of the McVittie spacetime does not exist in a global sense.

Although dense sphere packings serve as the initial state for simulations in powder metallurgy, additive manufacturing and granular physics, the choice of a packing generator is rarely guided by a systematic benchmark. A representative packing must be (i)-(ii) physically admissible (non-overlapping particles in gravitational equilibrium); (iii) faithful to the target particle size distribution (PSD); (iv) representative in relative density Phi; and (v) computationally affordable. Four open-source tools have been benchmarked, meeting (i)-(ii) by construction: the sequential DR (dropping-and-rolling) and its densified variant DR-ME, and the discrete element method (DEM) codes LAMMPS (gravity) and dp3D (isostatic compression). Across four configurations (2D/3D lognormal, 3D bimodal, and a 3D domain-size study), they are compared against an industrial MIM-grade powder, with PSD fidelity measured by the bin-width-independent Hellinger distance and Phi against the feedstock solid loading (phi_exp = 0.62, by Archimedes' method). In 3D, the DEM codes reach the densest packings but run more than three orders of magnitude slower: for approximatively 20 000 particles, DR shows a 9% phi shortfall relative to dp3D while running 1800x faster. These idealised model packings yield application-driven tool-selection guidelines.

Pressure-induced superconductivity in bilayer nickelates provides a platform for investigating intertwined roles of charge/spin orders and electric transport in unconventional superconductivity. However, important quantitative information on the transport, such as the absolute value of the resistivity, the anisotropy, and the scattering rate of carriers, remains insufficient due to the lack of accurate measurements using large single crystals. Here we establish a high-precision pressure-temperature phase diagram of high-quality La3Ni2O7 single crystals, by measuring the in-plane and out-of-plane resistivities. We resolve two distinct anomalies associated with density-wave formation with contrasting pressure dependences. The pressure-induced structural transition enhances not only the resistivity values for both directions, but also its anisotropy at low temperatures, demonstrating a pronounced effect of density-wave order on the charge dynamics. Superconductivity with zero-resistance emerges near the boundary where the density-wave phases are fully suppressed, and above Tc, the resistivity exhibits a temperature-linear dependence over a wide temperature range while the scattering rate falls within a regime of the Planckian limit. Our results show that pressure dramatically changes the anisotropic charge transport via modifying density-wave orders, and eventually produces a pronounced strange-metal state with strong scatterings, from which superconductivity develops. This establishes robust density-wave correlations and Planckian dissipation as remarkable features of La3Ni2O7.

In magnetic materials angular momentum can be mediated by different carriers, including electrons, magnons, and phonons. The magnons can interact with circularly polarized phonons which are close in energy, provided specific symmetry conditions are met. If the interaction is strong enough the phonons and magnons shift in frequency and start to mix to form hybrid magneto-elastic quasi-particles. In this paper, we develop a constrained Hamiltonian framework which incorporates hybrid stiffness matrices and Berry curvatures. We quantify the degree of hybridization between phonons and magnons (up to 8% in CrI$_3$ and 25% in CrBr$_3$) using a decomposition of the total energy, which is a generalization of the norm decomposition for atomic contributions to phonons used in the literature. We also explore how the total angular momentum is conserved but shared between the phononic and magnonic subsystems upon hybridization.