Papers

Long-range Rydberg molecules (LRMs) form when a highly excited Rydberg electron scatters from ground-state atoms inside its orbit, creating oscillatory, long-range potentials. We present a combined theoretical and experimental study of caesium dimers correlated to 402P3/2 Rydberg states, with an emphasis on decay via autoionisation (associative ionisation). Our model includes a relativistic treatment of electron-atom scattering with spin-orbit coupling, the perturber's hyperfine structure, and coupling of vibrational levels to a continuum of short-range decay channels. Calculated potential-energy curves predict two families of wells: outer wells near the classical outer turning point supporting long-lived states, and inner wells at shorter range whose lifetimes are limited by tunneling and subsequent vibronic decay. Using photoassociation in an ultracold Cs gas and an analysis of pulsed-field-ionisation signals which are highly selective for the detection of molecules, we assign resonances by binding energy and measure lifetimes. The measured lifetimes of inner-well states increase systematically with increasing detuning and agree with calculated lifetimes; detection of Cs2+ product ions supports autoionisation as a dominant channel. We show that the lifetimes are strongly reduced by spin-orbit interactions in the transient Cs-collision complex, which lift the near-degeneracy in Omega observed for states in the outer well and control the inner-well binding. The identified states also provide promising pathways to create ultracold molecules in ion-pair states.

Permanents, hafnians, and loop-hafnians are combinatorial matrix functions defined on the perfect matchings of graphs. These matrix functions arise in the quantum amplitudes for output photon configurations in linear optical networks, which has led both to sampling hardness results and to a variety of applications. The same matrix functions also appear in Ising models, since Ising models naturally encode graphs, and the matrix functions emerge from quantum interference. However, previous results have been restricted to balanced bipartite graphs, and hence to the permanent. In this paper, we extend the graph structure of the Ising model so that all three matrix functions can be obtained. In particular, self-loops can be encoded by introducing additional spins, and the loop-hafnian is obtained from a superposition of spin states, which we call the loop-hafnian state. We show how to recover the hafnian and permanent cases as the number of self-loops decreases. Moreover, we show that the Ising dynamics we consider can be simulated on a quantum circuit with \(O(n^2)\) gates for a \(2n \times 2n\) real symmetric matrix. Specifically, we propose a quantum circuit that maps the loop-hafnian state to a basis state with \(O(n^2)\) gates, and show that, by an adaptive measurement, the loop-hafnian states for all subsets of vertices can be measured simultaneously.

Despite the unitary evolution of closed quantum systems, long-time expectation of local observables are well described by thermal ensembles, providing the foundation of quantum statistical mechanics. A promising route to understanding this quantum thermalization is the eigenstate thermalization hypothesis (ETH), which posits that individual energy eigenstates already appear locally thermal. Subsequent studies have extended this concept to the full ETH, which captures higher-order correlations among matrix elements through nontrivial relations. In this work, we perform a detailed exact-diagonalization study of finite-size corrections to these relations in the canonical ensemble. We distinguish two distinct sources of corrections: those arising from energy fluctuations, which decay polynomially with system size, and those originating from fluctuations within each energy window, which decay exponentially with system size. In particular, our analysis resolves the puzzle that, for certain observables, finite-size corrections exhibit anomalous growth with increasing system size even in chaotic systems. Our results provide a systematic and practical methodology for validating the full ETH in quantum many-body systems.

Lindbladian quantum master equations (LEs) are the most popular descriptions for quantum systems weakly coupled to baths. But, recent works have established that in many situations such Markovian descriptions are fundamentally limited: they cannot simultaneously capture populations and coherences even to the leading-order in system-bath couplings. This can cause violation of fundamental properties like thermalization and continuity equations associated with local conservation laws, even when such properties are expected in the actual setting. This begs the question: given a physical situation, how do we know if there exists an LE that describes it to a desired accuracy? Here we show that, for both equilibrium and non-equilibrium steady states (NESS), this question can be succinctly formulated as a semidefinite program (SDP), a convex optimization technique. If a solution to the SDP can be found to a desired accuracy, then an LE description is possible for the chosen setting. If not, no LE description is fundamentally attainable, showing that a consistent Markovian treatment is impossible even at weak system-bath coupling for that particular setting. Considering few qubit isotropic XXZ-type models coupled to multiple baths, we find that in most parameter regimes, LE description giving accurate populations and coherences to leading-order is unattainable, leading to rigorous no-go results. However, in some cases, LE description having correct populations but inaccurate coherences, and satisfying local conservation laws, is possible over some of the parameter regimes. Our work highlights the power of semidefinite programming in the analysis of physically consistent LEs, thereby, in understanding the limits of Markovian descriptions at weak system-bath couplings.

Determining the position of an entity is a fundamental prerequisite for nearly all activities. Classical means, however, have been proven incapable of providing secure position verification, meaning that a prover can mislead verifiers about its actual position. In this work, we propose and experimentally realize a secure position-verification protocol that leverages quantum optics and relativity within an information-theoretic framework. Using phase-randomized weak coherent states, two verifiers separated by 2 km securely verify the prover's position with an accuracy better than 75 meters. These results establish secure position-based authentication as a practical possibility, paving the way for applications in financial transactions, disaster response, and authenticated secure communications.

We investigate the dissipative dynamics of a two-level atom in a weak gravitational field. Using the Feynman–Vernon influence functional formalism, we derive a quantum master equation describing the two-level atom interacting with a scalar field in a Newtonian gravitational field, and compute the energy dissipation rate of the atom. We find that the spontaneous emission rate (the dissipation rate in vacuum) is modified by the gravitational field. Specifically, this modification depends on the atom's dipole, the position of the atom relative to the source of the gravitational field, and the frequency of the scalar radiation emitted by the atom. Furthermore, we identify the parameter regimes in which the spontaneous emission rate is enhanced or suppressed by gravity. We also discuss how the modification arises from time dilation and dipole radiation in a weak gravitational field. These findings provide a theoretical basis for exploring gravitational effects in open quantum systems.

Product-formula (PF) based quantum simulation is a promising approach for simulating quantum systems on near-term quantum computers. Achieving a desired simulation precision typically requires a polynomially increasing number of Trotter steps, which remains challenging due to the limited performance of current quantum hardware. To alleviate this issue, post-processing techniques such as the multi-product formula (MPF) have been introduced to suppress algorithmic errors within restricted hardware resources. In this work, we propose a dual-channel multi-product formula that achieves a two-fold improvement in Trotter error scaling. As a result, our method enables the target simulation precision to be reached with approximately half the circuit depth compared to conventional MPF schemes. Importantly, the reduced circuit depth directly translates into lower physical error mitigation overhead when implemented on real quantum hardware. We demonstrate that, for a fixed CNOT count as a measure of quantum circuit, our proposal yields significantly smaller algorithmic errors, while the sampling error remains essentially unchanged.

Coulomb-Sturmian functions are complete, orthonormal, and include the full spectrum of continuum states. They are restricted to integer values of quantum numbers, as imposed by boundary and orthonormality conditions. Bagci-Hoggan exponential-type orbitals remove this restriction through a generalization to quantum number with fractional order. The differential equations for N-dimensional Bagci-Hoggan orbitals are derived. It is demonstrated that Coulomb-Sturmian functions satisfy a particular case of these equations. Additionally, Guseinov's Psi-alpha-ETOs are identified as N-dimensional Coulomb-Sturmians with a shifted dimensional parameter alpha, rather than representing an independent complete orthonormal sets of basis in a weighted Hilbert space.

Optimal control is a central problem in quantum thermodynamics. When minimizing dissipated work and work fluctuations defined via the two-point measurement scheme in open quantum systems, existing approaches largely focus on the rapid- and slow-driving limits, leaving the behavior at intermediate timescales elusive. In this work, by numerically optimizing the driving protocols, we demonstrate that the open quantum systems exhibit distinct optimal structures not captured by the conventional limits. Specifically, in the coherent spin-boson model, we find that the optimal protocol switches discontinuously between distinct locally optimal solutions as the relative weight between dissipation and fluctuations is varied. Furthermore, for a single-level quantum dot coupled to a fermionic reservoir, the optimized protocol develops a characteristic multi-step structure.

Computing electronic structures of quantum systems is a key task underpinning many applications in photonics, solid-state physics, and quantum technologies. This task is typically performed through iterative algorithms to find the energy eigenstates of a Hamiltonian, which are usually computationally expensive and suffer from convergence issues. In this work, we develop and implement the Quantum Jacobi-Davidson (QJD) method and its quantum diagonalization variant, the Sample-Based Quantum Jacobi-Davidson (SBQJD) method, and demonstrate their fast convergence for ground state energy estimation. We assess the intrinsic algorithmic performance of our methods through exact numerical simulations on a variety of quantum systems, including 8-qubit diagonally dominant matrices, 12-qubit one-dimensional Ising models, and a 10-qubit water molecule (H$_2$O) Hamiltonian. Our results show that both QJD and SBQJD achieve significantly faster convergence and require fewer Pauli measurements compared to the recently reported Quantum Davidson method, with SBQJD further benefiting from optimized reference state preparation. These findings establish the QJD framework as an efficient general-purpose subspace-based technique for solving quantum eigenvalue problems, providing a promising foundation for sparse Hamiltonian calculations on future fault-tolerant quantum hardware.

In finite-dimensional quantum systems, temperature cannot be uniquely defined. This, in turn, implies that there are several ways to define entropy production in finite-dimensional quantum systems, because the classical entropy production depends on temperature. We propose a unified definition of entropy production based on the difference in quantum relative entropy with respect to reference states characterized by effective temperatures. We demonstrate that the proposed definition naturally decomposes into a Clausius-type entropy production and an additional contribution arising from the time dependence of the effective temperature. Furthermore, we show that requiring the entropy production rate to take the conventional form as the sum of the entropy change and the heat flow constrains the effective temperature to be either constant or equal to a specific energy-matching effective temperature. For general initial states, entropy production can become negative, in which case we derive lower bounds on entropy production and establish sufficient conditions for its non-negativity using the trace distance.

Non-Hermitian systems have attracted significant theoretical interest due to their extreme properties. However, realizations have mostly been limited to classical applications or artificial setups. In this study, we focus on the quantum nature inherent in bosonic Bogoliubov-de Gennes (BdG) systems, which from the perspective of spectral theory corresponds to non-Hermiticity. Based on this insight, we propose a steady-state skin effect in quantum condensed matter utilizing such BdG non-Hermiticity. Specifically, we introduce BdG quantum terms arising from parametric pumping to the edge states of an underlying bosonic Hermitian Chern insulator, thereby realizing non-Hermiticity without dissipation. This system design has the advantage of being largely independent of microscopic model details. Through analysis using non-equilibrium Green's functions, we find that under open boundary conditions, a steady state exhibiting the non-Hermitian skin effect is realized. The pronounced corner particle accumulation observed in this steady state shows quadrature anisotropy, which manifests the bosonic quantum nature. Our results bridge the gap between the fascinating mathematics of non-Hermitian matrices and practical quantum physical systems.

Spin relaxometry using solid-state spin defects, such as the diamond nitrogen-vacancy (NV) center, probes dynamical processes by measuring how environmental fluctuations enhance the spin relaxation rate. In the weak-coupling limit, relaxation rates sample the transverse magnetic-noise power spectral density through a sensor-specific filter function, turning the defect into a local, frequency-selective noise spectrometer. This review bridges theory and experiment, clarifying how measured relaxation rates map onto noise spectra and how near-field geometry shapes the response. We highlight representative applications across condensed-matter physics, chemical and biological sensing, and relaxometry-based magnetic-resonance spectroscopy. We conclude with emerging opportunities and key challenges.

In recent years, quantum computing technologies have steadily matured and have begun to find practical applications across various domains. One important area is network communication security, where Quantum Key Distribution (QKD) enables communicating parties to establish a shared secret that can then be used to generate symmetric keys for subsequent encryption and decryption. This study focuses on implementing and comparing two well-known QKD protocols, namely BB84 and E91, within an actual quantum computing environment. It also proposes the use of SX gate operations to generate uniform quantum superposition states. By leveraging the properties of quantum superposition and quantum entanglement, the study illustrates how communicating parties can securely obtain a shared secret while preventing adversaries from intercepting it. The experiments are conducted using the IBM Quantum Platform to demonstrate the feasibility of the BB84 and E91 protocols on actual quantum hardware. The evaluation considers several metrics, including entropy, Independent and Identically Distributed (IID), and error-rate verifications.

In a previous work we were able to define a non-additive measure that can be used to represent both classical and quantum states in physics. We further extended this idea to work on a generic space of statistical ensembles (i.e. an ensemble space) in a way that connects to Choquet theory. The question of which non-additive integral is suitable to generalize the notion of expectation value remains open. In this paper we show that the Sugeno and Choquet integrals are not suitable.

Satellite-based quantum communications represent a critical advancement in the pursuit of secure, global-scale quantum networks. Leveraging the principles of quantum mechanics, these systems offer unparalleled security through Quantum Key Distribution (QKD) and other quantum communication protocols. This review provides a comprehensive overview of the current state of satellite-based quantum communications, focusing on the evolution from terrestrial to space-based systems. We explore the distinct advantages and challenges of discrete-variable (DV) and continuous-variable (CV) quantum communication technologies in the context of satellite deployments. The paper also discusses key milestones such as the successful implementation of quantum communication via the Micius satellite and outlines the primary challenges, including atmospheric turbulence and the development of quantum repeaters, that must be addressed to achieve a global quantum internet. This review aims to consolidate recent advancements in the field, providing insights and perspectives on the future directions and potential innovations that will drive the continued evolution of satellite-based quantum communications.

Self-accelerating Airy matter waves offer a clean setting to access the intrinsic cubic-in-time phase. Here we reconstruct the relative phase of simulated Airy-shaped Bose-Einstein condensates from interference fringes in free space, a regime approached in microgravity. The cubic phase dynamics are quantified via approximately debiased, windowed polynomial fits with systematics-aware uncertainty estimates that account for window-induced correlations. We compare two physically feasible phase-extraction methods, heterodyne-based and density-based, and show that an Airy-Gaussian geometry yields substantially improved robustness to fit-window selection relative to an Airy-Airy collision. In the weakly interacting regime, the extracted cubic coefficient responds linearly to leading order to the effective interaction strength, its shift from the noninteracting value providing a calibrated probe of weak mean-field nonlinearities in self-accelerating condensates.

Non-Hermitian degeneracies of Lindblad generators (Liouvillian exceptional points) can induce non-exponential relaxation and higher-order poles in dynamical response functions. A collective spin coupled to a polarized Markovian bath exhibits an exceptional spectral phase in which defective Liouvillian modes imprint super-Lorentzian features in frequency-resolved spectra. We compute the emission spectrum via the Liouvillian resolvent, identify symmetry-sector selection rules, and demonstrate that exceptional-point signatures are strongly state-dependent: they are suppressed in steady-state fluorescence yet become unambiguous for generic (infinite-temperature or random) initial states. Our results provide an experimentally accessible spectroscopic diagnostic of many-body Liouvillian exceptional phases and clarify when steady-state emission can (and cannot) reveal them.

Optimisation plays a central role in a wide range of scientific and industrial applications, and quantum computing has been widely proposed as a means to achieve computational advantages in this domain. To date, research into the design of noise-resilient quantum algorithms has been dominated by variational approaches, while alternatives remain relatively unexplored. In this work, we introduce a novel class of quantum optimisation heuristics that forgo this variational framework in favour of a hybrid quantum-classical approach built upon Markov Chain Monte Carlo (MCMC) techniques. We introduce Quantum-enhanced Simulated Annealing (QeSA) and Quantum-enhanced Parallel Tempering (QePT), before validating these heuristics on hard Sherrington-Kirkpatrick instances and demonstrate their superior scaling over classical benchmarks. These algorithms are expected to exhibit inherent robustness to noise and support parallel execution across both quantum and classical resources with only classical communication required. As such, they offer a scalable and potentially competitive route toward solving large-scale optimisation problems with near-term quantum devices.

This paper presents an efficient and scalable tensor network framework for quantum kernel circuit simulation, alleviating practical costs associated with increasing qubit counts and data size. The framework enables systematic large-scale evaluation of a linearly mixed quantum-classical dual kernel of up to 784 qubits. Using Fashion-MNIST, the classification performance of the test dataset is compared between a classical kernel, a quantum kernel, and the quantum-classical dual kernel across the feature dimensions from 2 to 784, with a one-to-one mapping between encoded features and qubits. Our result shows that the quantum-classical dual kernel consistently outperforms both single-kernel baselines, remains stable as the dimensionality increases, and mitigates the large-scale degradation observed in the quantum kernel. Analysis of the learned mixing weights indicates that quantum contributions dominate below 128 features, while classical contributions become increasingly important beyond 128, suggesting that the classical kernel provides a stabilizing anchor against concentration effects and hardware noise while preserving quantum gains at lower dimensions.

I review the formalism, Feynman rules, and combinatorics that constrain a field to propagate “classically", strictly in tree diagrams, either by itself, or interacting with other, purely quantum fields. The perturbation theory is reorganized by virtue of the linear terms that introduce the constraints via Lagrange multipliers, generalizing and giving results that cannot be obtained with the standard procedures which start at the quadratic terms. I apply the formalism to a theory of an $O(N)$-symmetric quantum field interacting with a “classical" scalar field via cubic interactions in six spacetime dimensions. Using the renormalization group, I examine the effective potential, symmetry breaking with radiative corrections, the fixed points in $d=6-\epsilon$ dimensions, and compare with other works. Other possible generalizations and applications of the formalism are also discussed.

In this manuscript, we study the quench dynamics of a transverse-field XY model starting from coherent Gibbs states. The results reveal that the initial strength of magnetization plays a crucial role in the emergence of dynamical quantum phase transitions. In concrete terms, when quenching within the same phase, through the properties of observables such as Fisher zeros and magnetization, we show that the stronger the initial magnetization, the more difficult the emergence of dynamical quantum phase transitions. The underlying mechanism is that the strong initial magnetization provides a directional effect, which inhibits the spin flipping in the process of quantum quench, making the dynamical quantum phase transition difficult to emerge. Since dynamical quantum phase transitions can be experimentally realized in various artificial systems, we hope that the physics predicted here can be experimentally verified in tabletop platforms.

According to the island formula, information in the code subspace defined in the black hole interior is embedded in the Hawking radiation after the Page time. At first sight, this embedding suggests that operations acting on the Hawking radiation could modify the information in the code subspace, potentially leading to an apparent violation of causality. Indeed, in previous studies based on the PSSY model, which incorporates only the topological degrees of freedom of gravity, it was shown that when the error is sufficiently large, a violation of causality can arise, as indicated by a nonvanishing mutual information. In this paper, we investigate the situation in which dynamical gravity also acts on the Hawking radiation. In this case, operations on the Hawking radiation induce nontrivial backreaction on the bulk spacetime appearing in the gravitational path integral for the mutual information – an effect that is absent when the Hawking radiation is non-gravitating. We find that this backreaction renders the relevant mutual information vanishing. This result implies that, in theories with dynamical gravity, the apparent violation of causality is resolved.

We develop an information-theoretic framework connecting stability, privacy, and generalization for quantum learning algorithms. Learning procedures are modeled as quantum instruments with classical-quantum outputs, and losses are represented by observables. We prove that under a classical-quantum sub-Gaussian condition, an information-theoretic stability measure controls the expected generalization error. Furthermore, we establish a high-probability generalization bound using quantum Rényi divergences to manage higher-order dependencies under non-commutativity. In the trusted Data Processor setting, quantum differential privacy (QDP) provides a mechanism for stability. We show that one-neighbor QDP strictly bounds the information leaked by the classical-quantum output. Combining this with our stability theorem yields a direct privacy-to-generalization guarantee. We also explore an untrusted Data Processor setting. Here, output privacy alone is insufficient since an adversarial processor could perform a highly informative procedure before applying noisy post-processing. To combat this, we introduce Information-Theoretic Admissibility (ITA), a certification condition ensuring the prescribed procedure is not just a degraded version of a strictly more informative, physically allowed operation on the encoded ensemble. We prove a fundamental separation: while admissibility and privacy are in strong tension in classical models, quantum non-orthogonality makes them compatible. A quantum measurement can be ITA - exhausting all relevant accessible information - without perfectly recovering the classical dataset. We illustrate this separation through a concrete quantum ITA example.

Sample-based quantum diagonalization (SQD) offers a powerful route to accurate quantum chemistry on noisy intermediate-scale quantum (NISQ) devices by combining quantum sampling with classical diagonalization. Here we introduce HSQD, a novel half-qubit SQD approach that halves the qubit requirement for simulating a chemical system and drastically reduces overall circuit depth and gate counts, suppressing hardware noise. When modeling the dissociation of the nitrogen molecule with a (10e, 26o) active space, HSQD matches the accuracy of SQD on IBM quantum hardware using only half the number of qubits and 40% fewer measurements. We further enhance HSQD with a heat-bath configuration interaction (HCI) inspired selection of the samples, forming HCI-HSQD. This yields sub-millihartree accuracy across the N2 potential energy surface and produces subspaces up to 39% smaller than those from classical HCI, showing a significant improvement in the compactness of the ground-state representation. Finally, we demonstrate the scalability of HCI-HSQD using iron-sulfur clusters, reaching active spaces of up to (54e, 36o) while using only half as many qubits as the original SQD. For these systems, HCI-HSQD reduces SQD energy errors by up to 76% for [2Fe-2S] and 26% for [4Fe-4S], while also reducing subspace sizes, halving measurement requirements, and eliminating expensive post-processing. Together, these results establish half-qubit SQD as a noise-resilient and resource-efficient pathway toward practical quantum advantage in strongly correlated chemistry.

We extend shapelets for the analysis of astronomical images to be available in a phase space, introducing Wigner function shapelets(WFSs). Whereas conventional shapelets expand images separately in configuration or Fourier space using Hermite-Gaussian or Laguerre-Gaussian modes, WFSs represent images directly in the four-dimensional phase space with symplectic group $\mathrm{Sp}(4,\mathbb{R})$, which is quantized by a phase-space cell $2\pi\lambdabar$ that determines a resolution limit of a telescope. WFSs consist of a bilinear form of the cross-Wigner function of the Laguerre-Gaussian modes as an orthogonal and complete basis for the Wigner function of an image, carrying out $\mathrm{SU}(2)$ irreducible representations of the phase space with the Hopf tori. We introduce a scalar function $\mathcal{W}_{k\ell} (Q_0,Q_2)$ from the $\mathrm{U}(1)\times \mathrm{U}(1)$ - covariant tori to a two-dimensional space of constants of motion $(Q_0,Q_2)$ – the harmonic energy and axial angular momentum – thereby yielding a natural phase-space “band structure”, given a pair of winding number $(k,\ell) \in \mathbb{Z}^2$. % The WFSs leverage key properties of the Wigner function for image analysis: (i) it encodes full information of an image in a symmetry-preserving way; (ii) its trasport equation naturally involves with a Liouville equation at $\lambdabar \rightarrow 0$; (iii) it admits positive/negative oscillatory patterns on $(Q_0,Q_2)$ plane that can be sensitive spatial coherent structure of galaxy morphology and cosmological imprints; and (iv) systematics and noise can be manipulated as a quantum channel operation. This paper aims to bring together all the formulae related to the Wigner function in the context of astrophysics and cosmology, formally organizing them in terminologies of both astronomy and of quantum information theory.

Decoherence is strongly influenced by environmental criticality, with conventional Hermitian critical points typically enhancing the loss of quantum coherence. Here, we show that this paradigm is fundamentally altered in non-Hermitian environments. Focusing on qubits coupled to non-Hermitian spin chains and interacting ultracold Fermi gases, we find that approaching exceptional points can either enhance or strongly suppress decoherence, depending on the balance between Hermitian and non-Hermitian system-environment couplings. In particular, when these couplings are comparable, decoherence is dramatically suppressed at exceptional transitions. We trace this behavior to the distinct response of the environmental ground state near non-Hermitian degeneracies and demonstrate the robustness of this effect across multiple models. Finally, we show that the predicted suppression of decoherence is directly observable on current digital quantum simulation platforms. Our results establish exceptional points as a concrete mechanism for suppressing decoherence and identify non-Hermitian criticality as a new avenue for coherence control in open quantum systems and quantum technologies.

This paper surveys various results in the field of Quantum Learning theory, specifically focusing on learning quantum-encoded classical concepts in the Probably Approximately Correct (PAC) framework. The cornerstone of this work is the emphasis on query, sample, and time complexity separations between classical and quantum learning that emerge under learning with query access to different labeling oracles. This paper aims to consolidate all known results in the area under the above umbrella and underscore the limits of our understanding by leaving the reader with 23 open problems.

Why does quantum theory need the complex numbers? With a view toward answering this question, I argue that the usual Hilbert-space formalism is a special case of the general method of Markovian embeddings. I then describe the `indivisible interpretation' of quantum theory, according to which a quantum system can be regarded as an `indivisible' stochastic process unfolding in an old-fashioned configuration space, with wave functions and other exotic Hilbert-space ingredients demoted from having an ontological status. The complex numbers end up being necessary to ensure that the Hilbert-space formalism is indeed a Markovian embedding.

We investigate the effects of Hawking radiation on quantum steering and steering asymmetry in a tripartite system embedded in Schwarzschild spacetime. All tripartite steering types were classified,comprising three "1 to 2" and three "2 to 1" steering cases. Through a systematic analysis of all physically relevant scenarios (including accessible and inaccessible modes), we classify three canonical scenarios with one, two and three physically accessible modes. In the scenario of three physically accessible modes, Hawking radiation disrupts quantum steering, with the maximum steering asymmetry during the two-way steering to one-way steering transition precisely demarcating the phase boundary between these regimes. For two physically accessible modes, Hawking radiation exhibits dual behavior: enhancing the steering from Alice and Bob to anti-Charlie under certain parameters while suppressing it under others, while net strengthening other steering types. When considering one physically accessible mode, the Hawking effect of the black hole significantly enhances quantum steering. These findings provide new insights into quantum correlations in curved spacetime and establish observable signatures of Hawking effects in quantum steering phenomena.