Papers

We determine the optimal quantum manipulation protocols for implementing high-fidelity, fast single-qubit gates. We demonstrate that the time-optimal pulse sequence is a “bang-bang” sequence: discrete pulses of either positive or negative maximum amplitude or zero. The non-adiabatic bang-bang pulse sequence minimizes gate duration providing a speedup for low-frequency architectures. We first derive the protocol for a transversally driven two-level system, extracting exact analytic expressions for minimum gate times. We then extend this framework to general multilevel architectures, identifying conditions that enable the coherent suppression of leakage errors. Using the fluxonium circuit as a representative case study, we optimize $X/2$ and $Y/2$ gates through a combination of discrete bang sequences and continuous waveform smoothing. This approach preserves near-optimal execution speeds while mitigating transitions outside the computational subspace. Open-system simulations demonstrate that these sequences outperform commensurate and resonant pulse schemes across different fluxonium regimes, achieving low-error manipulation significantly faster than standard resonant control, even in the presence of $1/f$ flux noise and dissipation.

The process-matrix formalism describes quantum processes without assuming a fixed global causal order, but the physical meaning of indefinite causal order remains open. We address this question within histories theory, where temporal ordering and dynamical evolution are distinct structures. We show that the decoherence functional admits coherent superpositions of dynamics. In measurement settings, under suitable factorization conditions, these generate process matrices. Thus, we identify process-matrix indefinite causal order as an operationally restricted realization of the more general phenomenon of superposition of dynamics. Histories theory also admits a distinct kinematical notion of indefinite order, in which the ordering of physical events is itself a history observable. The resulting distinction between event order and intervention order clarifies the relation between indefinite causal order, quantum dynamics, and spacetime structure.

We demonstrate a static, disorder-free spin chain Hamiltonian which, by tuning the coupling to an auxiliary qubit, realizes a crossover between super-ballistic, ancilla-accelerated scrambling and sub-ballistic operator confinement. Our minimal model is the mixed-field Ising chain with a spin-1/2 ancilla coupled to its longitudinal magnetization. The ancilla mediates an effective all-to-all interaction which accelerates operator spreading and entanglement growth when weakly coupled, but rapidly saturates its entanglement and projects spin chain operators into effectively frozen subspaces when the ancilla coupling is strong. We locate this crossover independently through both a divergent peak in the mutual-information saturation time near $\lambda_c N/h \approx 8$ and an exponential suppression of the late-time OTOC growth rate, $\log\alpha \propto 1/\lambda$. Through a Feshbach-Fano projection and Schrieffer-Wolff transformation, we reveal an effective hidden symmetry on the chain which confines operators on the chain for a time exponential in the coupling strength. This reconciles the fast, $\log(N)$ scrambling reported for random-unitary-circuit realizations of the star geometry with the confinement previously found in its time-independent Hamiltonian analog, showing both emerge from a single Hamiltonian family as a function of one dimensionless parameter.

The sensitivity of miniaturized optically pumped magnetometers (OPMs) is limited by weak atom-light coupling, which an optical cavity can enhance. In this work, we model the photon-shot-noise-limited sensitivity of cavity-enhanced OPMs in the regime of strongly collisionally broadened optical transitions, characteristic of buffer-gas-filled miniaturized vapor cells. The cavity enhancement is benchmarked against a single-pass free-induction-decay OPM employing Faraday-rotation readout, with the probe power and detuning jointly optimized using the Cramér-Rao lower bound as a figure of merit. For a Fabry-Pérot cavity, we compare side-of-fringe, homodyne, Pound-Drever-Hall, and Faraday-rotation readout. All four yield an optimal sensitivity enhancement scaling as $\alpha\sqrt{2\mathcal{F}/\pi}$, where $\mathcal{F}$ is the cavity finesse and $0.5\leq \alpha \leq 1$ is a readout-dependent prefactor. The enhancement is maximized at critical coupling, and we quantify its degradation away from this point. We further show that, despite spin-dependent absorption associated with the ensemble's vector polarizability, near-critical coupling can be maintained throughout spin precession at arbitrary finesse by exceeding a derived probe-power threshold and increasing the atomic detuning with finesse. We also establish a limit to the maximum cavity enhancement set by vector light-shift noise.

Atomistic descriptions of hydrogen diffusion and trapping at defects are essential for understanding hydrogen embrittlement. As the lightest solute in metals, hydrogen exhibits nuclear quantum effects that alter these processes even at room temperature. Explicit treatment of such effects is computationally demanding, limiting large-scale simulations of complex environments. Here, we use an Fe-H machine-learning interatomic potential (MLIP) based on the performant implementation of the atomic cluster expansion (PACE), covering diverse Fe-H environments, and relabel the training configurations underpinning its transferability with quantum mean forces from centroid-constrained path-integral molecular dynamics at 300 K. This yields a nuclear-quantum-corrected PACE (NQC-PACE) without additional density functional theory calculations. At parent PACE, NQC-PACE describes nuclear quantum effects on hydrogen trapping at vacancies, dislocations, surfaces and general grain boundaries, H-H interactions, and diffusion in alpha-Fe. Grand-canonical Monte Carlo/molecular dynamics simulations show nuclear quantum effects markedly enhance hydrogen segregation at general grain boundaries and trapping behaviour in closer agreement with experimental trends. This enhancement arises from selective quantum stabilisation of open, anisotropically soft local environments. Our framework uses finite-temperature quantum mean forces to relabel the configurational space covered by an MLIP, enabling large-scale analysis of complex materials where light-element quantum effects matter.

In this paper, we analytically compute the quench dynamics of simple product states that evolve under a spin-1/2 XY-model with power-law interactions in a staggered magnetic field. We use linear spin-wave theory as well as a low-energy theory that accounts for strong phase fluctuations. This allows us to provide an analytic expression for the dynamics of phase-type correlations with short-range interactions, where the linear spin-wave theory is known to fail. For density-type correlations the low-energy theory and the linear spin-wave theory give the same result. Our analysis is valid for arbitrary magnetic field strength and in particular provides an analytical understanding of the quench dynamics near criticality. Our results can thus provide useful tools for the experimental observation of excitation gaps and critical behavior.

While the impossibility of perfectly identifying non-orthogonal states is a cornerstone of quantum information science, they can nevertheless be distinguished with a finite error probability. Here, we propose a two-reservoir quantum machine driven by such non-orthogonal-state discrimination to map its functional boundaries across the parameter space of the state overlap $\mu$ and the Carnot efficiency $\eta_C$. Within this $\eta_C$-$\mu$ plane, the machine exhibits phase-transition-like functional switching among a pure heat-engine phase, a mixed phase, and a dissipative phase. We identify critical thresholds governing these transitions: strong thermal driving ($\eta_C>0.5$) unconditionally guarantees positive work extraction, whereas weak driving ($\eta_C \lesssim 0.13$) induces an anomalous reentrant transition, where increasing $\mu$ unexpectedly restores engine functionality after a purely dissipative regime. Our results explicitly demonstrate how quantum mechanics and thermodynamics jointly constrain information-to-energy conversion.

The corrections to Casimir pressure (force) on two plates of real metal described by the dielectric constant of the Drude, Drude-Lorentz, Drude-Smith and plasma models at the finite temperature and conductivity are considered. We consider the contour integral lost in the Lifshitz formula and discuss inconsistencies for the corrections to the Casimir force given in the literature. Numerical and analytical results in the near and far zones are presented. The two first-order correction obtained due to the actual dielectric constant of the metal coincided with the literature data.

We investigate the emergence of a scarred discrete time crystal (SDTC) phase in a periodically driven dimerized spin chain. While generic interacting Floquet systems are expected to thermalize according to the eigenstate thermalization hypothesis (ETH), we demonstrate that this system hosts quantum many-body scars (QMBS) that induce a regime of weak ergodicity breaking. Through an analysis of Floquet level statistics, entanglement entropy, and eigenstate fidelity, we identify a manifold of low-entanglement states characterized by semi-Poisson statistics embedded within an otherwise thermal spectrum. These scarred states support robust subharmonic oscillations with period doubling, signaling the spontaneous breaking of discrete time-translation symmetry. We show that the SDTC response is robust against a variety of initial state configurations, demonstrating its stability beyond fine-tuned conditions. A finite-size scaling analysis reveals that the time-crystalline lifetime grows with system size within the range accessible to our exact-diagonalization calculations. However, drawing on the general phenomenology of approximate many-body scars, we expect that hybridization between Floquet scars and the thermal continuum will eventually curtail this growth, causing the lifetime to saturate at system sizes beyond our current numerical reach. This characterizes the SDTC as a long-lived metastable dynamical regime rather than a strictly stable thermodynamic phase, providing a comprehensive framework for understanding the interplay between periodic driving and constrained many-body dynamics in disorder-free systems.

Quantum Feedback Network Theory (also known colloquially as the SLH-framework) is a powerful tool for modeling systems in quantum optics. This paper describes a method for discretising SLH-networks such that they can be implemented on a quantum computer. Building on a discretisation theorem for quantum stochastic differential equations, we establish strong convergence, uniformly on compact time intervals, of the discretised system to the continuous SLH-network as the discretisation mesh goes to zero.

Anyons obey fractional statistics that lie between bosonic and fermionic statistics, giving rise to a broad range of intriguing phenomena. However, how anyonic statistics govern quantum-state complexity is still largely unexplored. In this work, we investigate the interplay between the statistical phase and on-site interactions in the anyon-Hubbard model, identifying exact quantum many-body scar eigenstates and novel quench dynamics. The Krylov complexity exhibits perfect periodic revivals independent of the statistical phase in the scarred dynamics, whereas after a quench it depends on both the statistical phase and the interaction strength. In the strong-interaction regime, we find approximate scarred dynamics, while in the weak-interaction regime the state spreads over Krylov space and the complexity ultimately saturates. Moreover, for the bosonic initial state, the complexity of fermions exhibits the lowest saturation value, and vice versa. For fractional statistics, the saturation plateau is minimized when the post-quench statistical phase is close to that of the initial state. Our results demonstrate the central role of the statistical phase in governing many-body dynamics and provide new insights into Krylov complexity and quantum many-body scars.

Experiments have revealed that collective strong coupling in optical cavities can drastically alter the conductivity and dielectric properties of molecular ensembles, or even trigger abrupt phase transitions (in Rayleigh scattering or dispersion force driven conformational equilibrium). Until now, the underlying physical mechanism has remained unknown. A recently proposed theory attributes these collective effects to cavity-induced electron correlations that share properties of the known spherical Sherrington-Kirkpatrick model of a spin glass. As a consequence, entropic effects can trigger a phase transition in the inter-molecular electron correlations that could potentially explain the above experiments. In the following work, we derive analytic expressions for the correlation (free-)energy-induced polarizability changes, which reveal that under certain conditions collective strong coupling can even render our molecular ensemble metallic by entropically ionizing collectively degenerate molecular orbitals. Our results are a first major step towards a holistic theory of cavity-mediated electron correlation effects.

The predictive simulation of molecules and materials has had a broad and significant impact. It nevertheless remains constrained by the cost of accurately treating electronic correlation, excited states, and complex energy landscapes. Quantum computing offers a fundamentally different computational paradigm in which quantum states are encoded and manipulated directly rather than approximated on classical hardware. Here we discuss where this approach may provide a genuine scientific advantage in chemistry, materials science, and biochemistry. Promising directions include the high-accuracy treatment of correlated active spaces, improved excited-state simulations, and accelerated exploration of combinatorial structure spaces. The central challenge is therefore not qubit scaling alone, but demonstrably chemically meaningful gains in predictive reliability. We argue that near-term value is most likely to come from disciplined workflow integration rather than wholesale replacement of classical methods. Noisy physical devices, error-mitigated utility experiments, early fault-tolerant devices, and fully fault-tolerant quantum computers offer different scientific prospects, and claims of usefulness must be tied to the specific regime being discussed. Quantum computing will become scientifically valuable when it demonstrably reduces uncertainty in computed energies, rates, spectra, or materials stability after the full costs of state preparation, measurement, error handling, and coupling to classical simulation are included.

Based on the van Kampen-Schram method, a formula for the Casimir force between two plane-layered structures is obtained, which includes the sum of the Matsubara frequencies and a contour integral. It is shown that the contour integral is not taken into account in the Lifshitz formula. Taking it into account gives a small contribution at high temperatures when the Lifshitz formula is approximately correct. However, at low temperatures, its contribution is significant, and the Lifshitz formula for the finite temperature does not transform into the Lifshitz formula at zero temperature. An example of interaction of metal plates is considered.1

We identify a non-Hermitian counterpart of dynamical tunnelling induced by a localised absorbing region in a symmetric confining potential. The absorbing region, described by an imaginary potential barrier, induces an effective double-well structure with a pair of long-lived states that have little population in the barrier. In contrast to the periodic tunnel oscillations through a real barrier, the system slowly evolves towards equal population of the two sides, while the absorbing region remains only weakly populated. We demonstrate this generic mechanism in two continuous potentials, and capture its essence analytically in a minimal $3\times 3$ matrix model. We propose a realistic optical waveguide setup for the experimental observation of the effect.

Quantum permutations, or magic unitaries, have in recent years been explored in the context of identifying `genuinely quantum' isometries of graphs. Here, we import this tool in physics, showing that quantum permutations yield a generalisation of quantum reference frames in a discrete setting that is reminiscent of the passage from special to general relativity. We show that the typical quantum reference frames framework corresponds to quantum permutations classified as `classical' in the mathematical literature, and which we demonstrate are quantum controlled transformations (superpositions of classical coordinate maps). Genuinely quantum permutations (i) allow to construct non–commuting quantum reference frames (ii) correspond to local, as opposed to global, superpositions of transformations. Strikingly, we find that the non-commutativity of quantum fields, when used as reference systems, is exactly what implies that the change of frame is achieved through a genuine quantum permutation. We illustrate the above with several examples in both first and second quantization formalism, which demonstrate (a) simultaneous control on non–commuting variables, (b) the existence of bipartite states that can be localized with a genuine quantum permutation and cannot be localized with the usual quantum reference frame transformations (without introducing additional degrees of freedom), (c) extension of the Ising model symmetries to genuinely quantum permutations, and (d) extension of the symmetries of a scalar field action on curved spacetime to genuinely quantum permutations. While we have in mind applications in quantum gravity, we expect our formalism to be of interest in a wide range of topics in quantum information.

Two-photon bound states in the continuum provide a stringent setting in which quantum interference must suppress radiative decay despite photon–photon interactions. Here, we prove a no-go theorem showing that a single nonlinear mode linearly coupled to a noninteracting bosonic continuum cannot support an interaction-active two-photon bound state in the continuum, provided the local spectral density vanishes only at isolated energies. Apart from this regularity condition, the result is independent of the continuum dispersion and the frequency-dependent coupling. Although such an exact bound state is forbidden, long-lived two-photon quasi-bound states remain possible. For a giant Kerr cavity nonlocally coupled to a waveguide, we derive the asymptotic scaling of the two-photon decay rate with the coupling-point separation from weak nonlinearity to the hard-core limit. We further identify a regime in which the two-photon resonance is substantially longer lived than the single-photon excitation. Our work establishes the limits of exact multiphoton confinement while opening a route to engineering long-lived interacting few-photon states.

We establish spectral filtering as a control parameter for the conditional momentum and position correlations of SPDC biphotons used in quantum imaging. The conditional momentum uncertainty in spontaneous parametric down-conversion (SPDC) is strongly crystal-length and spectral-filter dependent along the walk-off axis; therefore, the effect is observed only in critically phase-matched (CPM) crystals such as $\beta$-barium borate (BBO) crystals, while quasi-phase-matched (QPM) crystals and the non-walk-off axis of BBO remain scaled strictly according to the standard pump waist size ($w_0$) dependent scaling law $1/w_0$-independent of the filter. In position space, the spectral-filter effect is universal and produces a flat-dip-rise (FDR) profile in every crystal class examined. Although this FDR profile was previously demonstrated in BBO only in the nondegenerate regime, our results establish its generality: the FDR dip is also present at exact degeneracy in QPM crystals, an unexpected feature that was previously thought to deliver a resolution advantage in the nondegenerate regime alone. Our treatment applies to all SPDC-based quantum-imaging regimes (including CPM and QPM crystals, walk-off and non-walk-off axes, degenerate and nondegenerate emission, and signal and idler filtering) and offers enhanced quantum-imaging resolution via the design rules presented-most prominently along the walk-off axis of CPM crystals (in far field) and across all transverse axes in the near field.

Higher-form symmetries generalize conventional global symmetries and act on lower-dimensional submanifolds of a quantum system. While Abelian topological phases can be organized by 1-form symmetries that form a group, non-Abelian topological phases based on finite groups require 1-form symmetry operators governed by non-invertible fusion algebras. In this work, we make this statement precise in quantum double lattice models $\mathcal D(G)$ for finite non-Abelian groups $G$. We construct the electric, magnetic, and dyonic 1-form operators directly at the lattice fixed point and show that together they form a complete nonlocal diagnostic algebra for the topological Hilbert space. Using these operators, we explicitly determine the cylinder and torus ground-state subspaces for arbitrary finite $G$. Furthermore, we calculate the microscopic fusion and gluing of the 1-form symmetries and show that their topological deformation properties emerge after projection to the defect-free topological subspace. Our results establish ground states of non-Abelian quantum double models as a concrete microscopic realization of spontaneous non-invertible 1-form symmetry breaking and provide an operator language that may be useful for characterizing such states in quantum processors.

We present MoMPy, an open-source Python package that provides a unified, declarative construction of moment matrices for semidefinite programming (SDP) hierarchies. The user declares operator labels together with a small set of structural relations, and MoMPy returns a matrix of SDP variable indices in which every implied identification has already been made, ready for CVXPY or any other modelling layer. Internally, the identification problem is recast as a word-rewriting problem on tuples of integers and solved with a memoised breadth-first closure coupled to a disjoint-set forest, so that each distinct monomial is processed exactly once per build, however many matrix entries it eventually labels. The central abstraction is independent of the physical scenario: the same construction handles tracial, state (NPA), and block-valued moments, with no notion of parties, settings or preparations. We demonstrate this generality on a tripartite Mermin inequality, the bipartite CHSH inequality, measurement compatibility in a steering scenario, state discrimination and dimension witnessing in prepare-and-measure scenarios, device-independent randomness certification, and certification of deterministic correlations from known ensembles. The construction is validated against an independent brute-force implementation and benchmarked across eight structurally distinct scenarios, spanning bipartite and tripartite Bell tests, heterogeneous-outcome, steering-type, jointly-measurable and network configurations.

Detector-based first-arrival models and operator-based arrival-time observables provide two conceptually distinct approaches to the quantum time-of-arrival problem. Here we connect these approaches by using the positive unidirectional first-arrival current of Marchewka and Schuss (MS) as the basis for constructing a family of positive arrival-time operators and, after normalization, an arrival-time POVM. The first-arrival character is inherited from the dynamics with a Dirichlet boundary, while the detector response is introduced at the amplitude level through a non-negative spectral function $\lambda(k)$. For a one-dimensional free particle, resolution of the identity uniquely selects $\lambda(k)=\frac{\pi}{4k}$, yielding an arrival-time POVM. The normalized unidirectional current can therefore be interpreted as a generalized Born-rule probability density rather than as a hazard rate, as in the original MS formulation. We further show that the normalized one-sided MS amplitude coincides, up to an irrelevant phase, with the corresponding momentum-sector amplitude in Kijowski's arrival-time distribution. Consequently, the two directional contributions in Kijowski's distribution can be reproduced by two first-arrival problems defined on opposite sides of the boundary. This establishes an equivalence between the resulting arrival-time statistics while emphasizing the different physical interpretations: MS describes a first-arrival process with coherent momentum components, whereas Kijowski's POVM treats the two momentum sectors as separate directional contributions. Thus, the normalized unidirectional current provides a POVM description of first arrival at a point, combining detector-based dynamics with a positive-operator measurement structure.

Quantum position-verification is a cryptographic task wherein a verifier attempts to establish the location in space of a prover. Recent experiments have implemented a well-studied class of position-verification schemes, known as the $f$-BB84 scenario, but their security under realistic loss and imperfect state preparation remains incompletely understood. We give a new lower bound on any attack on this scheme, in particular proving nearly-linear quantum gate lower bounds on the attacker, even when allowing the attacker to declare a transmission loss of up to $50\%$, allowing for the challenges prepared by the referee to be imperfect, and allowing the quantum messages used in the protocol to be arbitrarily slow. Our results are applicable to recent and upcoming experimental implementations of $f$-BB84, and in particular establish their security under a bounded quantum gate assumption on the attacker. The key ingredient is a tight analytic tradeoff for a lossy BB84 monogamy-of-entanglement game, valid without requiring the attackers' outputs to agree, which replaces observations previously made numerically.

We determine the exact strong converse exponent of quantum soft covering under the sandwiched Rényi divergence for all orders $\alpha\in[\frac{1}{2},\infty)$. For $\alpha\in[\frac{1}{2},1)$, the exponent is characterized by the two-parameter club-sandwiched mutual information, whereas for $\alpha\in[1,\infty)$, it is characterized by the order-$\alpha$ sandwiched Rényi mutual information. We also determine the exact strong converse exponent of privacy amplification against quantum side information under the sandwiched Rényi divergence for $\alpha\in(2,\infty)$, expressed in terms of the corresponding order-$\alpha$ sandwiched Rényi conditional entropy. To the best of our knowledge, these results provide the first exact characterization of the strong converse exponent of quantum soft covering and the first precise operational interpretation of the two-parameter club-sandwiched mutual information in the quantum setting. The key ingredient is that we establish the exponential rate of the $K$-functional, which is instrumental in deriving the strong converse exponent of quantum soft covering for $\alpha\in[\frac{1}{2},1)$.

We show that the maximum winning probability in the repeated CHSH game for classical strategies is at most $\left(\frac{1+\sqrt{5}}{4}\right)^n$. Together with a matching lower bound by Barak et al. (FOCS'2008), this determines the asymptotic value of the repeated CHSH game exactly. We also show that, if an XOR game has a gap between classical and quantum values, there is also a gap between the asymptotic values for the repeated game.

Rayleigh-Ritz on a subspace from sample-based quantum diagonalization yields the spectral function of a probe state; we ask what bounds its error. The captured Born weight w does not do so alone: an analytic counterexample forces any bound indexed on it to its trivial value. On Hubbard rings at broadening $\eta=0.18t$, the probe's whole support (w=1) still has relative error 0.35–0.43: Rayleigh-Ritz displaces the poles a subspace keeps. The error is bounded instead by the missed weight and the Hamiltonian coupling across the subspace boundary. Given the exact ground state, we prove for any orthogonal projector that the relative $L_1$ error lies between 1-w and $\min\{1+w,(1+\sqrt{w})(\sqrt{1-w}+\Lambda(\eta)/\eta)\}$, with $\Lambda$, that coupling at resolution $\eta$, computed from the Ritz pairs and the retained probe. We test it on L-site rings with classical Born draws of an exactly time-evolved probe, or their infinite-shot ranking; our device runs are execution records only. On the ranked subspaces the bound is vacuous at every operating fraction of the resource scan, at all five sizes up to a 10306296-determinant sector: its leakage branch exceeds the trivial bound 1+w by factors 1.66–8.61. It is informative only well above those fractions, where it is calibrated: at L=6–12 it first beats the trivial bound at a true error of (1.2–7.3)$\times10^{-3}$, though an unproven fit in 1-w is tighter on 11 of 14 splits. For the Born-ranked $A(k,\omega)$ at 85% of its sector, the bound is 350–900 times the true error. In the determinant basis, zero leakage at full weight needs every symmetry-allowed determinant: the probe's measured Krylov support at L=6–14. Plane-wave orbitals reduce a momentum probe's support to one block, 1/L of an exponentially growing sector. At full weight a second-order certificate follows; below it, where every operating fraction lies, one is open.

The exact diagonalization of the nuclear shell model scales exponentially, leading to severe memory bottlenecks in classical high-performance computing (HPC). While hybrid quantum algorithms like the Variational Quantum Eigensolver (VQE) aim to overcome these limits, their deep quantum circuits and iterative feedback loops are susceptible to substantial noise inherent in the current Noisy Intermediate-Scale Quantum (NISQ) hardware. This noise renders several algorithms, such as the VQE, impractical for large-scale calculations despite sophisticated noise-mitigation techniques. As a pragmatic approach tolerant to these issues, we apply the Sample-based Quantum Diagonalization (SQD) framework to nuclear shell models for the first time. Using $^{38}\text{Ar}$ as a benchmark to confirm the numerical accuracy, we extend SQD to $^{32}\text{Mg}$, solving a nuclear shell-model Hamiltonian whose underlying Hilbert space cannot be directly diagonalized using conventional classical methods in a given HPC system. We present a systematic comparison of SQD with standard variational quantum schemes and exact classical solvers. By leveraging NISQ hardware connected via the cloud to classical HPC clusters, the SQD-based scheme could outperform conventional supercomputers in memory scaling and total execution time, enabling more rigorous large-scale shell model calculations.

The Virtual Quantum Subspace Expansion (VQSE) extends the Variational Quantum Eigensolver (VQE) by leveraging additional measurements on the reference state to capture the influence of excluded virtual orbitals. This makes VQSE attractive for chemical applications where accurate energy differences along potential energy surfaces are crucial for modeling reaction rates and kinetics. In this work, we analyze VQSE performance on H$_2$ dissociation including references that use Hartree–Fock molecular orbitals with broken spin symmetry. We identify two mechanisms which affect accuracy: overlap of the reference state with the exact full configuration interaction (FCI) wavefunction and operator pool expressivity. We show these mechanisms are strongly co-dependent. When operators are restricted to act only from the active to the virtual space, results become highly sensitive to the reference, and enlarging the active space does not guarantee improved accuracy. In this case, prioritizing reference overlap over energy minimization is therefore essential. Adding single excitations and number operators within the active space recovers the accuracy of MR-CISD (multi-reference configuration interaction singles and doubles) regardless of the reference. In our noisy hardware experiments, we achieve chemical accuracy by adding additional operators and using strict regularization. These findings motivate careful co-design of reference fidelity, pool expressivity, and hardware constraints for practical VQSE deployment.

The fluctuation theorem (FT) relates microscopic trajectory distributions to macroscopic averages. Using quasiprobability trajectories, this framework has recently been extended to quantum information dynamics. Despite sharing the form of conventional thermodynamic FTs, information FTs involve quasiprobabilities whose negativity and statistical properties remain poorly understood. We analytically determine the properties of negative distributions based on the two-qubit Bell-diagonal states evolving under local dephasing, depolarizing, and amplitude-damping channels. We show that the occurrence of negative quasiprobabilities is not determined by the entanglement or Bell nonlocality of the initial state. Instead, a diagonal-interference decomposition of the transition quasiprobability provides a sufficient condition for negativity. More importantly, we prove that quasiprobability distributions containing negative weights can violate Horváth's necessary and sufficient criterion for the Jensen-Steffensen inequality to hold for every continuous convex function, even though the specific exponential Jensen-like relation enforced by the FT and the quantum data-processing inequality remains valid. This trajectory-level contrast with classical stochastic descriptions reveals a distinctive feature of quantum information dynamics that is invisible when only the corresponding average-level information inequalities are considered.

Point-particle descriptions of electron-electron interaction omit the coherent longitudinal extent of a free-electron quantum wave packet (QEW). A relativistic two-electron wave-packet theory identifies the longitudinal QEW size as a direct control parameter for entanglement generated by mutual electromagnetic coupling. For two electrons in spatially separated paths, the quadratic interaction phase gives a dimensionless entangling parameter $\Gee$ and the logarithmic negativity $\EN=\operatorname{arsinh}(\Gee)/\ln2$. Full time-dependent Schrödinger equation calculations confirm the scaling for narrow QEWs and reveal higher-order Coulomb effects at larger spatial extent. Free drift increases the interaction-point QEW width while preserving the momentum probability distribution, thereby enhancing the subsequently generated entanglement. The results establish a direct connection between free-electron wave-particle duality and bipartite entanglement.

We investigate how nonequilibrium quantum initial states can reduce the finite-time thermodynamic cost of Landauer erasure. Considering a general finite-dimensional quantum memory coupled to a thermal reservoir via a Davies generator, we show that the dissipated heat at a fixed operational erasure fidelity is largely controlled by the overlap of the initial state with the slowest Liouvillian relaxation mode. We derive a modified finite-time Landauer bound in which the excess dissipation above the quasistatic limit scales quadratically with this slow mode projection, and we prove a sufficient Mpemba Landauer condition under which a hotter state can erase faster and dissipate less heat than a colder preparation, without violating Landauer's principle. A minimal qutrit model illustrates these quantum Mpemba speedups and reveals broad parameter regimes where coherence and Hamiltonian-induced modes conspire to suppress finite-time entropy production. We further identify practical control knobs, including temperature tuning, coherence engineering, and Hamiltonian shaping of Liouvillian spectra, which enable Mpemba-enhanced erasure to be directly tested on platforms such as superconducting circuits, trapped ions, semiconductor quantum dots, and solid-state spins.