Papers

The space Inv$(j_1,j_2,j_3,j_4)$ of SU(2)-invariant four-valent tensors, also known as intertwiners, can be understood as the quantum states of a tetrahedron in Euclidean space with fixed areas. In loop quantum gravity, they are states of the smallest “atom of space" with non-zero volume. At the same time they correspond to four-party tensor product states invariant under global rotations. We consider the multipartite entanglement of states in Inv$(j_1,j_2,j_3,j_4)$ using the recently proposed entropic fill. Numerically evaluating entropic fill in the case of equal spins between $1/2$ and $11$, we find that the distributions of entanglement are very different for intertwiners as compared to generic tensors, and for coherent intertwiners as compared to generic ones. The peak in the distribution seems to be at the highest entanglement for generic intertwiners and at the lowest for generic tensors, but in terms of average entanglement, the roles are switched: average entanglement is highest in arbitrary tensors and lower in intertwiners, at least in the regime of large $j$. We also consider the dependence of entanglement on the geometric data of coherent intertwiners. Coherent intertwiners labeled by a symmetric tetrahedron have maximal entanglement, while minimal entanglement implies a coherent state labeled by a degenerate tetrahedron. Away from these cases, entanglement depends on the geometry in an intricate way.

The Mollow triplet is the definitive spectral signature of an optically dressed quantum emitter. We predict that for emitters coupled to localized phonons, this signature is not confined to the zero-phonon line. Under a strong resonant drive, we show that Mollow triplets are strikingly replicated on the associated phonon sidebands -a surprising result, given that phonon sidebands are typically viewed as incoherent, inelastic scattering pathways. These vibronic Mollow triplets are a direct fingerprint of dynamically generated dressed states that hybridize the emitter's electronic, photonic, and vibrational degrees of freedom. We develop a scalable analytical formalism to model this effect in complex, multi-mode molecular systems, such as dibenzoterrylene. Our work provides the precise driving conditions for observing these novel spectral features, establishing a new signature of coherence in vibronically coupled systems.

Negating the Measurement Independence assumption (MI) is often referred to as the `third way' to account for the experimental violation of Bell's inequality. However, this route is generally viewed as ludicrously contrived, implying some implausible conspiracy where experimenters are denied the freedom to choose measurement settings as they like. Here, a locally realistic model of quantum physics is developed (Rational Mechanics - RaQM - based on a gravitational discretisation of Hilbert Space) which violates MI without denying free will. Crucially, RaQM distinguishes experimenters' ability to freely choose measurement settings to some nominal accuracy, from an inability to choose exact settings, which were never under their control anyway. In RaQM, Hilbert states are necessarily undefined in bases where squared amplitudes and/or complex phases are irrational numbers. Such `irrational' bases correspond to conceivable but necessarily impossible counterfactual measurements, and are shown to play a ubiquitous role in the analysis of both single- and entangled-particle quantum physics. It is concluded that violation of Bell inequalities can be understood with none of the strange processes historically associated with it. Instead, using concepts from (non-classical) $p$-adic number theory, we relate RaQM to Bohm and Hiley's concept of a holistic Machian-like Undivided Universe. If this interpretation of Bell's Theorem is correct, building more and more energetic particle accelerators to probe smaller and smaller scales, in the search for a theory which synthesises quantum and gravitational physics and hence a Theory of Everything, may be a fruitless exercise.

We experimentally demonstrate the simultaneous estimation of the three parameters characterizing a pair of incoherent optical sources in the sub-Rayleigh regime, enabling super-resolved scene characterization. Using spatial-mode demultiplexing (SPADE) with two demultiplexers–one deliberately shifted–we determine separations well below the diffraction limit and achieve sensitive joint estimation of separation, centroid, and relative brightness over a broad range of scene configurations in a single experimental setting. We benchmark our performance using Fisher-information-based Cramér-Rao bounds, and discuss the corresponding quantum limits. We investigate two complementary scenarios: a realistic case with slightly non-identical sources, and an idealized case of indistinguishable sources.

Nonlocal light-matter interactions with giant atoms in high-dimensional environments are not only fundamentally intriguing for testing quantum electrodynamics beyond the dipole approximation but also crucial for building high-dimensional quantum networks and engineering multipartite entangled states. Given the enigmatic and largely uncharted collective signatures exhibited by multiple giant atoms within two-dimensional optical baths, we delve into their nonperturbative collective dynamics within the single-excitation subspace, focusing on the case where they are coupled to a common two-dimensional photonic reservoir and employing a resolvent operator approach. We demonstrate that precisely engineered atomic arrangements lead to unconventional quantum dynamics, featuring non-Markovianity-induced beats and long-lived bound states in the continuum, thereby providing a versatile platform for implementing two-dimensional quantum memory. Phenomenologically, we observe the emergence of exotic photon emission patterns in both two- and three-dimensional (3D) baths. The emission directions are shown to be precisely controllable on demand through exact phase engineering of the coupling parameters, enabling a highly efficient chiral light-matter interface. Moreover, our generalization to a 3D bath reveals that coherent dipole-dipole interactions can survive despite the coupling to a continuum of modes, a finding that challenges conventional wisdom regarding decoherence.

We analyse the performance of different error handling methods in the quantum data centre paradigm of distributed quantum computing. We compare the impact of quantum error detection, using the three-qubit repetition code and the [[4, 1, 2]] Leung-Nielsen-Chuang-Yamamoto code, on remote gates with that of conventional entanglement distillation techniques. Detailed classical simulation is used to obtain results for realistic near-term hardware.

Robust quantum routing is essential for scalable quantum technologies. This paper investigates the resilience of routing protocols in network architectures designed for perfect, high-fidelity transfer of both classical and quantum information under ideal conditions. We encode information in the position of a quantum walker on a graph, modelling the routing of a generic qubit state from a single input to multiple (orthogonal) outputs. We analyse and assess routing performance in various regimes, evaluating their robustness against static and dynamical noise.

Understanding the thermal behavior of quantum many-body pure states is one of the most fundamental issues in quantum thermodynamics. It is widely known that typical pure states yield vanishing work, just as thermal states do, when one restricts to local operations that cannot access correlations among subsystems. However, it remains unclear whether this equivalence to thermal states persists under LOCC (local operations and classical communication), where classically accessible correlations can be exploited for work extraction. In this work, we establish criteria for determining whether many-body pure states remain equivalent to thermal states even under LOCC, and show that this thermal equivalence is governed by their multipartite quantum correlation structure. We show that states with asymptotically maximal multipartite entanglement, such as Haar-random states, cannot yield extensive work under LOCC, whereas some states with limited multipartite entanglement, such as constant-degree graph states, allow extensive work extraction despite being locally indistinguishable from thermal states. Thus, our work provides a refined operational notion of thermal equivalence beyond the traditional local regime, which is becoming increasingly important due to the recent expansion of experimentally accessible operations.

In this paper, we develop a blended dynamics framework for open quantum networks with diffusive couplings. The network consists of qubits interconnected through Hamiltonian couplings, environmental dissipation, and consensus-like diffusive interactions. Such networks commonly arise in spontaneous emission processes and non-Hermitian quantum computing, and their evolution follows a Lindblad master equation. Blended dynamics theory is well established in the classical setting as a tool for analyzing emergent behaviors in heterogeneous networks with diffusive couplings. Its key insight is to blend the local dynamics rather than the trajectories of individual nodes. Perturbation analysis then shows that, under sufficiently strong coupling, all node trajectories tend to stay close to those of the blended system over time. We first show that this theory extends naturally to the reduced-state dynamics of quantum networks, revealing classical-like clustering phenomena in which qubits converge to a shared equilibrium or a common trajectory determined by the quantum blended reduced-state dynamics. We then extend the analysis to qubit coherent states using quantum Laplacians and induced graphs, proving orbit attraction of the network density operator toward the quantum blended coherent dynamics, establishing the emergence of intrinsically quantum and dynamically clustering behaviors. Finally, numerical examples validate the theoretical results.

Localization in non-Hermitian quasicrystals can differ fundamentally from its Hermitian counterpart when non-reciprocity is spatially disordered. Here we study a one-dimensional non-Hermitian Aubry-André-Harper chain with a Bernoulli imaginary gauge field and quasiperiodic onsite modulation. In the nearest-neighbor limit, we identify an anomalous transition from a fully erratic non-Hermitian skin effect (ENHSE) phase to a fully localized phase. Although the fractal dimension vanishes in both regimes, the Lyapunov exponent and the fluctuation of the eigenstate center of mass sharply distinguish them. For generic finite-size realizations, this transition is further accompanied by a complex-to-real spectral change under periodic boundary conditions and a change of spectral winding from nontrivial to trivial. With weak next-nearest-neighbor hopping, we uncover an anomalous mobility edge at the same location as in the Hermitian generalized Aubry-André-Harper model, but separating Anderson-localized states from ENHSE-type macroscopic-accumulation states rather than extended states. We further show that this anomalous localization structure is reflected in spectral winding and wave-packet dynamics: single realizations exhibit winding-dependent drift, winding-resolved averaging preserves opposite directional responses, and full disorder averaging largely restores Hermitian-like transport. Our results establish practical diagnostics of anomalous localization and mobility edges in non-Hermitian quasicrystals.

We compare different circuits implementing distributed versions of quantum computations, using entangled pairs only, and using distributed fan-out operations (using GHZ states). We highlight the advantages of using distributed fan-out operations in terms of reductions in circuit depth and (possibly) entanglement resources. We note that distributed fan-out operations (or notably, distributed GHZ states) could be a “primitive” building block for distributed quantum operations in the same way as entangled pairs are, if distributed GHZ states could be realized efficiently.

Integrated optical microresonators exploiting either second-order ($\chi^{(2)}$) or third-order ($\chi^{(3)}$) nonlinearities have become key platforms for frequency conversion, low-noise microwave photonics, and quantum entanglement generation. Here, we present an analytic theory of Kerr-enhanced three-wave mixing amplification in an electro-optic microresonator with both $\chi^{(2)}$ and $\chi^{(3)}$ nonlinearities. We demonstrate that Kerr dressing hybridizes the optical sidebands, renormalizing the $\chi^{(2)}$ couplings and detunings. As a result the system exhibits gain in regions where analogous bare $\chi^{(2)}$ or $\chi^{(3)}$ amplifiers are subthreshold. Time-domain Langevin simulations confirm this threshold reduction, mapping a practical design window for experiments.

Fidelity estimation is a critical yet resource-intensive step in testing quantum programs on noisy intermediate-scale quantum (NISQ) devices, where the required number of measurements is difficult to predefine due to hardware noise, device heterogeneity, and transpilation-induced circuit transformations. We present QuFid, an adaptive and noise-aware framework that determines measurement budgets online by leveraging circuit structure and runtime statistical feedback. QuFid models a quantum program as a directed acyclic graph (DAG) and employs a control-flow-aware random walk to characterize noise propagation along gate dependencies. Backend-specific effects are captured via transpilation-induced structural deformation metrics, which are integrated into the random-walk formulation to induce a noise-propagation operator. Circuit complexity is then quantified through the spectral characteristics of this operator, providing a principled and lightweight basis for adaptive measurement planning. Experiments on 18 quantum benchmarks executed on IBM Quantum backends show that QuFid significantly reduces measurement cost compared to fixed-shot and learning-based baselines, while consistently maintaining acceptable fidelity bias.

Sensing networks underpin applications ranging from fundamental physics to real-world engineering. Distributed quantum sensing (DQS) can improve measurement performance, but existing protocols typically require multipartite entanglement, which poses substantial challenges for scalable implementation. Here, we introduce a DQS protocol based on bidirectional causal routing in a cyclic network, where a single probe sequentially interrogates M independent sensors along two opposite causal routes. By exploiting the noncommutativity between inter-sensor propagation and local sensing operations, the protocol turns propagation from a passive transport process into a source of sensing information, yielding an asymptotic 1/M^2 scaling of the estimation precision without multipartite entanglement. We experimentally demonstrate the protocol for distributed beam-tilt sensing in a free-space quantum optical network comprising up to 9 sensors, achieving picoradian-level precision in estimating the average tilt angle. These results identify propagation dynamics and routing geometry as active metrological resources for scalable distributed quantum sensing.

Quantum superposition is often phrased as the ability to add state vectors. In practice, however, the physical quantity is a ray (a rank-one projector), so each input specifies only a projector and leaves a gauge freedom in the phases of its vector representatives. This becomes a real operational barrier when one asks for a device that, given two independently prepared unknown pure states, outputs a coherent state proportional to a prescribed linear combination. We identify the missing ingredient as not probabilistic but phase-like. One needs a physical scenario that fixes a single phase convention on the relevant set of rays, so that the overlaps become well defined complex numbers. Thus, we formalize this through phase conventions and a single notion – dubbed as "overlap-determinability." Our main theorem gives an exact equivalence: A nonzero completely positive trace-nonincreasing map that probabilistically produces superposition on a domain exists if and only if that domain is overlap-determinable. This unifies modern no-superposition results and characterizes the exceptional yes-go protocols, which succeed precisely when side information supplies the required missing resource. We then show that granting universal access to such convention-fixed overlaps destabilizes the familiar foundational and computational constraints. It enables forbidden transformations akin to quantum cloning and yields super-luminal signaling. It would also permit reflections about unknown states, leading to exponentially fast overlap amplification and a collapse of Grover's search lower bound to a logarithmic query complexity.

The temporal degree-of-freedom of light is incredibly powerful for modern quantum technologies, enabling large-scale quantum computing architectures and record key-rates in quantum key distribution. However, the generalized measurement of large and complex quantum superpositions of the time-of-arrival of a photon remains a unique experimental challenge. Conventional methods based on unbalanced Franson-type interferometers scale poorly with dimension, requiring multiple cascaded devices and active phase stabilization. In addition, these are limited by construction to a restricted set of phase-only superposition measurements. Here we show how the coupling of spatial and temporal information inside a single multi-mode fiber can be harnessed to program completely generalized measurements for high-dimensional superpositions of photonic time-bin. Using the multi-spectral transmission matrix of the fiber, we find special sets of spatial modes that experience distinct dispersive delays through the fiber. By exciting coherent superpositions of these spatial modes, we engineer the equivalent of large, unbalanced multi-mode interferometers inside the fiber and use them to perform high-quality measurements of arbitrary time-bin superpositions in up to dimension 11. The single fiber functions as a scalable, common-path interferometer for time-bin qudits that significantly eases the experimental overheads of standard approaches based on unbalanced Franson-type interferometers, serving as an essential tool for quantum technologies that harness the temporal properties of light.

Integrated quantum photonics holds significant promise for scalable photonic quantum information processing, quantum repeaters, and quantum networks, but its development is hindered by the mismatch between materials hosting high-quality quantum emitters and those compatible with mature photonic technologies. Heterogeneous integration offers a potential solution to this challenge, yet practical implementations have been limited by inevitable insertion losses at material interfaces. Here, we present a self-aligned heterogeneous quantum photonic integration approach that can deterministically achieve near-unity coupling efficiency at the interface. To showcase our approach, we demonstrate Purcell enhancement of a silicon vacancy (SiV) center in diamond induced by a heterogeneous photonic crystal cavity defined by titanium dioxide (TiO2), as well as optical spin control and readout via a TiO2 photonic circuit. We further show that, when combined with inverse photonic design, our approach enables efficient and broadband collection of single photons from a color center into a heterogeneous waveguide. Our approach is not restricted to SiV centers or TiO2; it can be broadly applied to integrate diverse solid-state quantum emitters with thin-film photonic devices where conformal deposition is possible. Together, these results establish a practical route to scalable quantum photonic integrated circuits that combine high-quality quantum emitters with technologically mature photonic platforms.

The increasing number of control lines connecting quantum processors to external electronics constitutes a major bottleneck in the realization of large-scale quantum computers. Frequency-division multiplexing is expected to enable control of multiple qubits through a single microwave cable; however, interference from off-resonant microwave tones hinders precise qubit control. Here, we propose an active interference suppression method for frequency-division-multiplexed simultaneous gates on microwave-controlled qubits. We demonstrate that the deliberate incorporation of off-resonant microwave tones improves single-qubit gate fidelity. In particular, the gate infidelity scales inversely with the square of the number of microwave tones when off-resonant orthogonal or quasi-orthogonal tones are incorporated. Furthermore, we show that fast oscillations, neglected under the rotating wave approximation, degrade the gate fidelity, and that this degradation can be mitigated through optimized frequency allocation. The proposed approach is simple and effective for improving the performance of frequency-division-multiplexed quantum gates.

Photon distinguishability is a fundamental property manifested in multiphoton interference and one of the main sources of noise in any photonic quantum information processing. In this work, rather than relying on first-quantization methods, we build on a generating-function framework based on the phase-space formalism to characterize the effects of partial distinguishability on the interference of single-mode squeezed states. Our approach goes beyond commonly used models that represent distinguishability via additional noninterfering modes and captures genuine multiphoton interference effects induced by the overlap of the internal state of the photons. This description provides a clear physical account of how distinguishability gives rise to effective noise in Gaussian boson sampling protocols while enabling a systematic investigation of phase effects arising from the overlap of the internal states.

Quantum Key Distribution (QKD) allows secure communication without relying on computational assumptions, but can currently only be deployed over relatively short distances due to hardware constraints. To extend QKD over long distances, networks of trusted repeater nodes can be used, wherein QKD is executed between neighbouring nodes and messages between non-neighbouring nodes are forwarded using a relay protocol. Although these networks are being deployed worldwide, no protocol exists which provides provable guarantees of integrity against manipulation from both external adversaries and corrupted intermediates. In this work, we present the first protocol that provably provides both confidentiality and integrity. Our protocol combines an existing cryptographic technique, Algebraic Manipulation Detection (AMD) codes, with multi-path relaying over trusted repeater networks. This protocol achieves Information Theoretic Security (ITS) against the detection of manipulation, which we prove formally through a sequence of games.

This paper presents a Quantum Reinforcement Learning (QRL) solution to the dynamic portfolio optimization problem based on Variational Quantum Circuits. The implemented QRL approaches are quantum analogues of the classical neural-network-based Deep Deterministic Policy Gradient and Deep Q-Network algorithms. Through an empirical evaluation on real-world financial data, we show that our quantum agents achieve risk-adjusted performance comparable to, and in some cases exceeding, that of classical Deep RL models with several orders of magnitude more parameters. However, while quantum circuit execution is inherently fast at the hardware level, practical deployment on cloud-based quantum systems introduces substantial latency, making end-to-end runtime currently dominated by infrastructural overhead and limiting practical applicability. Taken together, our results suggest that QRL is theoretically competitive with state-of-the-art classical reinforcement learning and may become practically advantageous as deployment overheads diminish. This positions QRL as a promising paradigm for dynamic decision-making in complex, high-dimensional, and non-stationary environments such as financial markets. The complete codebase is released as open source at: https://github.com/VincentGurgul/qrl-dpo-public

This paper argues that several canonical puzzles in quantum mechanics, including spin measurement, the double slit, entanglement correlations, and Wigner's friend, share a common origin in a semantic error and the illicit promotion of interaction conditional outcomes to intrinsic properties. I introduce four principles that license only configuration relative predication, grounding outcomes in physical measurement geometry while preserving objectivity. Applying these principles uniformly dissolves each puzzle without new physics or ad hoc interpretive machinery. Bell's theorem and the Kochen-Specker theorem are reframed not as dynamical mysteries but as constraints on permissible explanatory structure, evidence that intrinsic-outcome semantics is incompatible with empirical reality. The result is a relational objectivity that avoids both naive property realism and observer subjectivism.

A result by Gioia and Wang [Phys Rev X 12, 031007 (2022)] showed that translationally symmetric states having nonzero momentum are necessarily long range entangled (LRE). Here, we consider the question: can a notion of momentum for non-translation symmetric states directly encode the nature of their entanglement, as it does for translation symmetric states? We show the answer is affirmative for 1D systems, while higher dimensional extensions and topologically ordered systems require further work. While Gioia and Wang's result applies to states connected via finite depth quantum circuits to a translation symmetric state, it is often impractical to find such a circuit to determine the nature of the entanglement of states that break translation symmetry. Here, instead of translation eigenstates, we focus on the many-body momentum distribution and the expectation value of the translation operator in many-body states of systems having broken translation symmetry. We show that in the continuum limit the magnitude of the expectation value of the translation operator $|<T>|$ necessarily goes to $1$ for delocalized states, a proxy for LRE states in 1D systems. This result can be seen as a momentum-space version of Resta's formula for the localization length. We investigate how accurate our results are in different lattice models with and without well-defined continuum limits. To that end, we introduce two models: a deterministic version of the random dimer model, illustrating the role of the thermodynamic and continuum limits for our result at a lattice level, and a simplified version of the Aubry-Andre model, with commensurate hopping for both momentum and position space. Finally, we use the random dimer model as a test case for the accuracy of $|<T>|$ as a localization (and thus entanglement) probe for 1D periodic lattice models without a well-defined continuum limit.

We formulate a deterministic algorithm for preparing a general $U(1)$-eigenstate of a spin-$s$ chain of length $n$. These states consist of linear combinations of computational basis states $|\vec{m}\rangle$ of $n$ qudits, each with $(2s+1)$ levels and $s= 1/2, 1, 3/2, \ldots$, whose ditstrings $\vec{m}$ have a fixed digit sum. Exploiting a Gray code for bounded integer compositions, whose consecutive ditstrings obey the Gray property, the quantum state is prepared by applying corresponding “Gray gates.” We use this algorithm to prepare exact eigenstates of integrable spin-$s$ XXX Hamiltonians. We also consider the preparation of AKLT states and spin-$s$ Dicke states.

A central unresolved question in the theory of superconductivity is why only a small subset of metallic elements exhibit a superconducting state, whereas many others remain strictly normal. Neither the conventional Bardeen Cooper Schrieffer (BCS) framework nor its extensions involving charge density wave (CDW) or pair density wave (PDW) order provide a predictive or material-selective criterion capable of distinguishing superconducting metals from non-superconducting ones. In particular, the persistent absence of superconductivity in simple noble metals with well-defined Fermi surfaces poses a challenge for all traditional approaches. Here we address this problem using the Modulated Electron Lattice (MEL) Ginzburg Landau (GL) framework introduced in our previous work. In this formulation, a coarse-grained MEL charge field $\rho_{\mathrm{MEL}}(\mathbf{r})$ with momentum dependent stiffness $\alpha(q)$ is coupled to the superconducting (SC) order parameter $\psi(\mathbf{r})$. We show that metallic superconductivity emerges only when the system satisfies a specific “MEL enhancement window,” characterized by a negative minimum of $\alpha(q)$ at either a finite modulation wave vector $q^{\ast}$ or at $q=0$, together with sufficiently strong coupling between $\rho_{\mathrm{MEL}}$ and $\psi$. This unified criterion naturally partitions metallic elements into three universal classes: (i) MEL-enhanced superconductors with a finite-$q^{\ast}$ charge mode, (ii) conventional BCS superconductors as the homogeneous $q^{\ast}=0$ limit of the MEL framework, and (iii) metals for which $\alpha(q)$ remains positive for all $q$, suppressing all MEL modes and preventing any superconducting instability. By applying this criterion to simple metallic elements, we identify why some metals develop superconductivity while others do not, possibly resolving a selection problem long open within the BCS paradigm.

Super-resolution (SR) seeks to reconstruct high-resolution (HR) data from low-resolution (LR) observations. Classical deep learning methods have advanced SR substantially, but require increasingly deeper networks, large datasets, and heavy computation to capture fine-grained correlations. In this work, we present the first study to investigate quantum circuits for SR. We propose a framework based on Variational Quantum Circuits (VQCs) with Adaptive Non-Local Observable (ANO) measurements. Unlike conventional VQCs with fixed Pauli readouts, ANO introduces trainable multi-qubit Hermitian observables, allowing the measurement process to adapt during training. This design leverages the high-dimensional Hilbert space of quantum systems and the representational structure provided by entanglement and superposition. Experiments demonstrate that ANO-VQCs achieve up to five-fold higher resolution with a relatively small model size, suggesting a promising new direction at the intersection of quantum machine learning and super-resolution.

Quantum state exclusion is the task of identifying at least one state from a known set that was not used in the preparation of a quantum system. A set of quantum states is said to admit state exclusion if there exists a measurement whose outcomes can be put in one-to-one correspondence with the states in the set, such that each outcome rules out its corresponding state with certainty (while possibly also ruling out other states), and each outcome occurs with nonzero probability for at least one state in the set. State exclusion, however, is not always possible in the single-copy setting. In this paper, we investigate whether access to multiple identical copies of the system enables state exclusion. We prove that for any set of three or more pure states, state exclusion becomes possible with a finite number of copies. Moreover, we show that the number of copies required may be arbitrarily large: in particular, for every natural number $N$, we construct sets of states for which state exclusion remains impossible with $N$ or fewer copies.

We extend the Pauli Propagation framework to imaginary-time evolution and introduce imaginary-time Pauli Propagation (ITPP), an operator-based algorithm for approximating thermal and ground-state properties directly in the Pauli basis. We derive explicit imaginary-time propagation rules for Pauli strings and analyze the main approximation errors arising from Trotterization and Pauli-space truncation, including bounds that account for the non-unitary nature of the evolution. Benchmarking ITPP on the one-dimensional transverse field Ising model, we find that the method performs best in the high-temperature regime, where thermal states admit a compressed Pauli-basis approximation, while the required number of Pauli terms grows rapidly as the system approaches the ground-state limit. We further compare with tensor network simulations and study finite-temperature scalability through the number of retained Pauli terms required to reach a fixed target accuracy. These results establish ITPP as a complementary framework for imaginary-time simulation and suggest that combining imaginary-time and real-time Pauli Propagation could provide a pathway toward simulating more general non-unitary and open quantum system dynamics within a unified Pauli-based framework.

We study parameter estimation for the transport coefficients of the quark-gluon plasma by differentiating open-quantum-system-based Monte Carlo simulations of quarkonium suppression. The underlying simulator requires solving a Lindblad equation in a large Hilbert space, which makes parameter estimation computationally expensive. We approach the problem using gradient-based optimization. Specifically, we apply the score-function gradient estimator to differentiate through discrete jump sampling in the Monte Carlo wave-function algorithm used to solve the Lindblad equation. The resulting stochastic gradient estimator exhibits sufficiently low variance and can still be estimated in an embarrassingly parallel manner, enabling efficient scaling of the simulations. We implement this gradient estimator in the existing open-source quarkonium suppression code QTraj. To demonstrate its utility for parameter estimation, we infer the two transport coefficients $\hat{\kappa}$ and $\hat{\gamma}$ using gradient-based optimization on synthetic nuclear modification factor data.

In quantum field theory, entanglement entropy under spatial bipartitioning serves as a powerful information-theoretic probe of quantum correlations. In this work, we present the first comprehensive numerical study of the dynamical evolution and geometric scaling of entanglement entropy in a nonperturbative, strong-field QED setting – specifically, in the context of the Sauter-Schwinger effect. While the weak-field regime is dominated by area-law states, we show that the entanglement entropy undergoes a transition from area-law to a volume-law scaling for certain strong-field regimes in the pulse-profile parameter space – signaling a fundamental shift in the underlying correlation structure induced by nonperturbative pair production. For intermediate regimes, the scaling is a power-law that interpolates between area- and volume-law behavior. Finally, we provide interpretations based on the behavior of the low-energy pair-creation spectrum and discuss how these insights could inform future investigations of related phenomena.