Enhancement of quantum battery performance is a popular subject in quantum thermodynamics. An interesting phenomenon is the quick charging effect [Phys. Rev. Res. 6, 023136 (2024)], which has been explored by utilizing a quantum interferometric technique known as superposition of trajectories. A similar technique used to boost quantum battery performance is indefinite causal order. Here, we propose a new charging protocol that utilizes cyclic indefinite causal order, whereby N charging sequences are superposed when utilizing N chargers. We observe charging efficiency bursts when implementing our cyclic indefinite charging protocol. The duration of these bursts increase with N. Additionally, we present a circuit model to implement our charging protocol for the two-charger scenario and perform proof-of-concept demonstrations on IonQ, Quantinuum and IBMQ quantum processors. The results validate the existence of charging efficiency bursts as shown by our theoretical analysis and numerical simulations.

The Hierarchical equations of motion (HEOM) method is an important non-perturbative technique, allowing numerically exact treatment of open quantum systems with strong coupling and non-Markovian memory. However, its encoding of bath memory into auxiliary density operators often limits direct access to detailed bath information. In contrast, the reaction-coordinate (RC) mapping allows direct and transparent access to the dominant collective bath mode, but its perturbative and often Markovian treatment of the residual bath restricts its reliability. In this work, we introduce RC-HEOM, a hybrid method that unifies the strengths of both approaches by combining RC mapping with a fully non-perturbative HEOM description of the residual bath. RC-HEOM simultaneously retains exact non-Markovian memory and access to the RC mode, which enables analysis of system-RC information. Applying this method to the Anderson impurity models, we directly track the emergence of the Kondo singlet from the growth of the Kondo resonance and uncover a nontrivial RC-mediated coherence revival. These results demonstrate that RC-HEOM is a promising method for characterizing open quantum systems in regimes that are difficult to access with conventional master-equation methods.

High-fidelity and rapid qubit readout is essential for superconducting quantum processors, typically realized through the quantum non-demolition (QND) dispersive interaction within a qubit-resonator architecture. However, the achievable readout speed and fidelity are fundamentally limited by measurement-induced state transitions (MIST). For a transmon qubit, MIST is highly sensitive to the offset charge $n_g$ due to the charge dispersion of its higher-lying energy levels. In this work, we systematically investigate $n_g$-dependent MIST dynamics governed by the diabaticity and symmetry of pulse shaping within a charge-sensitive transmon architecture. We engineer fast-load and fast-clear pulses that effectively suppress resonator photon overshoots, thereby demonstrating a highly practical strategy to mitigate MIST without requiring complex waveforms or real-time feedback. Utilizing active gate-voltage control and rapid feedback, the measurement-induced transition probability is precisely mapped against $n_g$ and the steady-state resonator photon number, exhibiting strong agreement with numerical Floquet branch analysis. Ultimately, we evaluate the $n_g$-averaged total error probabilities for both readout and post-readout stages, verifying that a straightforward three-step pulse scheme consistently minimizes overall readout errors. Within the framework of large-scale superconducting quantum processors, this practical, hardware-free approach inherently offers a better trade-off between the readout signal-to-noise ratio and QND preservation.

Quantum steering, measurement incompatibility, and instrument incompatibility have recently been recognized as unified manifestations of quantum incompatibility. Building on this perspective, we develop a general framework for constructing optimization-free, nonlinear incompatibility witnesses based on convex functionals, valid in arbitrary dimensions. We prove that these witnesses are nontrivial precisely when the underlying functional is non-affine on extremal points (e.g., pure states for ensembles). For pure bipartite states, the witnesses yield lower bounds on entanglement measures, thereby outperforming most linear steering inequalities in the pure-state regime. Moreover, the construction extends in full generality to certify measurement and instrument incompatibility, where the witnesses act as genuine incompatibility monotones. We demonstrate the versatility of our approach with two operationally relevant functionals: the Wigner-Yanase skew information and an ℓ_2-type coherence functional.

Spatio-temporal quantum steering provides a framework for benchmarking the nonclassicality of general quantum state transfer processes. A central diagnostic is the no-signaling-in-time (NSIT) condition, whose violation can indicate basis-dependent hardware errors. However, finite measurement statistics may also yield apparent violations, thereby obscuring the detection of basis-dependent hardware errors. To address this, we construct a statistical hypothesis test under the null hypothesis that NSIT violations arise solely from statistical fluctuations. Combining the statistical properties of NSIT violation under the null hypothesis with Chebyshev's inequality, we obtain a distribution-free upper bound on the $p$-value without parametric assumptions. We apply this method to two examples. For a single-qubit state-transfer experiment on a superconducting processor, we observe several instances that the NSIT violation is observed and the null hypothesis is simultaneously rejected by a small $p$-value, providing statistical evidence of basis-dependent hardware errors. For a seven-qubit Hayden-Preskill teleportation protocol on IonQ devices, the null hypothesis is also rejected even when the average fidelity exceeds the classical threshold, while the associated nonclassicality measure vanishes. Our results highlight the necessity of statistical hypothesis testing for detecting basis-dependent errors in near-term quantum devices.

Quantum information scrambling (QIS) describes the rapid spread of initially localized information across an entire quantum many-body system through entanglement generation. Once scrambled, the original local information becomes encoded globally, inaccessible from any single subsystem. In this work, we introduce a circuit-based decoding protocol. By utilizing the concept of postselected closed timelike curves (PCTCs), we demonstrate how postselection allows us to interpret an ordinary quantum experiment as an example of a paradox-free trajectory, simulating a consistent time loop and reliable information recovery. Specifically, when conditioned on a final postselected outcome, this experiment can be interpreted as decoding the scrambled information even before the original information is generated. Furthermore, the success probability of the PCTC is governed by out-of-time-ordered correlations, which is a standard measure of QIS. We experimentally implement our protocol on cloud-based Quantinuum and IBM quantum processors. Our approach illuminates a unique quantum task under postselection: the causally consistent simulation of future-to-past scrambled information retrieval.

The quantum kernel method, a promising quantum machine learning algorithm, possesses substantial potential for demonstrating quantum advantage. Although the majority of the quantum kernel is constructed in the context of gate-based quantum circuits, inspired by the idea of analog quantum computing, here we construct an analog quantum kernel and a hybrid quantum kernel, and show their competitiveness against other kernel methods in a benchmarking task and the practical problem of estimating non-Markovianity from sparse data. Additionally, we also incorporate operational noise into the quantum kernels. Our results reveal that the presence of operational noise can be beneficial to the performance of the developed quantum kernels. We attribute this counterintuitive noise-enhanced performance to the improved expressivity and higher model complexity induced by noise. These results pave the way for practical implementations of quantum kernel methods and provide an efficient approach for estimating non-Markovianity with reduced experimental demands.

We consider the problem of modeling a single qubit in contact with a one-dimensional waveguide beyond the standard perturbative and Markovian approximations. Using the recently developed input-output hierarchical equations of motion (io-HEOM), we investigate multiple examples of such waveguides, characterized by different spectral densities. Our examples highlight that the io-HEOM method can accurately capture non-Markovianity in waveguide QED from two distinct origins. The first source of non-Markovianity is spatially non-local coupling between the qubit and the waveguide. By examining two examples with non-local coupling, we show how the coupling function affects the steady-state bound photons, and demonstrate the release of these photons when the qubit energy is quenched. The second source of non-Markovianity is non-linear dispersion. We illustrate this scenario using the example of a cavity array with point-like coupling, where the non-linear dispersion leads to persistent oscillations due to Van Hove singularities in the spectral density.

One of the remarkable aspects of quantum steering is its ability to violate local uncertainty complementarity relations. In this vein of study, various steering witnesses have been developed. Here, we introduce a novel complementarity relation between the system’s quantum and classical uncertainties corresponding to the distillable coherence and the von Neumann entropy, respectively. We show that the proposed complementarity relation is tighter than the entropic uncertainty relation (EUR). Leveraging this result, we propose a steering witness that is more efficient than the EUR. From the operational perspective, the steering witness quantifies the amount of extra distillable coherence facilitated by quantum steerability. Notably, the proposed steering witness serves as a full entanglement measure for pure bipartite states–an ability that the EUR lacks. We also experimentally validate such a property through a photonic system. Furthermore, a deeper connection to the uncertainty principle is revealed by showcasing the steering-induced distillable coherence can quantify measurement incompatibility and quantum steerability under genuine incoherent operations. Our work establishes a clear quantitative and operational link between coherence and steering, which are vital resources of quantum technologies, and underscores our efforts in bridging the uncertainty principle with quantum coherence.

Temporal quantum correlations provide an intriguing way of testing quantumness at the macroscopic level, with a logical hierarchy present among the quantum correlations associated with nonmacrorealism, temporal steering, and temporal inseparability. By manipulating the dynamics of a superconducting qubit, we observe the full hierarchy of temporal quantum correlations. Moreover, we show that the rich dynamics of the temporal quantum correlations, such as sudden death or revival of temporal steering, provides a useful and unique measure for benchmarking qubits on a realistic circuit. Our work finds applications in identifying the casual structure in a quantum network, the non-Markovianity of open quantum systems, and the security bounds of quantum key distribution. As an example, we demonstrate the non-Markovianity of a single superconducting qubit on the quantum circuit.

This work addresses the critical challenge of dye molecule oxidation and its impact on device stability by investigating the suppression of photobleaching through engineering the chromophore's surrounding environment. We induce strong coupling with confined light modes [optical cavities or localized surface plasmons (LSPs)] to reduce the triplet state population, thereby mitigating photo-oxidation. Utilizing the hierarchical-equations-of-motion (HEOM) approach to capture non-Markovian and non-perturbative effects, we analyze both cavity-chromophore and LSP-chromophore systems. Our analysis reveals that the optimal antioxidation performance depends on the competition between cavity-chromophore coupling and cavity-bath (dissipation) interaction. Importantly, in the strong cavity-chromophore coupling regime, increasing cavity dissipation can enhance antioxidation through quantum coherence-induced population transfer. Conversely, in the weak cavity-chromophore coupling regime, increasing cavity dissipation can counterintuitively reduce the antioxidation capability due to a quantum Zeno-like effect. Furthermore, in the LSP-chromophore system, engineering the LSP structure, particularly the LSP dissipation rate, can similarly be used to optimize the antioxidation effect. These findings, complemented by analytical results for finding optimal system parameters, provide practical guidelines for designing photostable organic materials with enhanced performance in various optoelectronic applications.

We investigate non-Markovian transport dynamics and signatures of the Kondo effect in a single quantum dot (QD) model. The QD is coupled to a left lead non-Markovian bath and weakly coupled to a right lead, acting as a detector. We calculate the waiting time distribution (WTD) of electrons tunneling into the detector using a combination of the hierarchical equations of motion (HEOM) approach and the Born-Markov (BM) approximation. Oscillations emerge in the short-time WTD, becoming more pronounced with stronger left-lead coupling. Fourier analysis reveals a blueshift in the oscillation frequency as coupling increases, indicating enhanced system-bath hybridization. Crucially, comparison with a master equation confirms that these oscillations are a direct consequence of non-Markovian system-bath correlations. By introducing a toy model, we analyze the role of different parameters, such as coupling strength and bandwidth, in significantly influencing the oscillatory behavior in WTD. In addition, we examine the Kondo effect's influence on these oscillations by varying the coupling strength and the bandwidth of the non-Markovian bath. Decreasing the bandwidth and increasing the coupling strength enhance the WTD oscillations, while also enhancing the Kondo resonance in the quantum dot's density of states. Our results demonstrate that WTD oscillations offer a valuable tool for probing non-Markovian system-bath interactions and the emergence of Kondo correlations within QD systems.

We present QuantumToolbox.jl, an open-source Julia package for simulating open quantum systems. Designed with a syntax familiar to users of QuTiP (Quantum Toolbox in Python), it harnesses Julia's high-performance ecosystem to deliver fast and scalable simulations. The package includes a suite of time-evolution solvers supporting distributed computing and GPU acceleration, enabling efficient simulation of large-scale quantum systems. We also show how QuantumToolbox.jl can integrate with automatic differentiation tools, making it well-suited for gradient-based optimization tasks such as quantum optimal control. Benchmark comparisons demonstrate substantial performance gains over existing frameworks. With its flexible design and computational efficiency, QuantumToolbox.jl serves as a powerful tool for both theoretical studies and practical applications in quantum science.

This work addresses the critical challenge of dye molecule oxidation and its impact on device stability by investigating the suppression of photobleaching through engineering the chromophore's surrounding environment. We induce strong coupling with confined light modes [optical cavities or localized surface plasmons (LSPs)] to reduce the triplet state population, thereby mitigating photo-oxidation. Utilizing the hierarchical-equations-of-motion approach to capture non-Markovian and non-perturbative effects, we analyze both cavity-chromophore and LSP-chromophore systems. Our analysis reveals that the optimal antioxidation performance depends on the competition between cavity-chromophore coupling and cavity-bath (dissipation) interaction. Importantly, in the strong cavity-chromophore coupling regime, increasing cavity dissipation can enhance antioxidation through quantum coherence-induced population transfer. Conversely, in the weak cavity-chromophore coupling regime, increasing cavity dissipation can counterintuitively reduce the antioxidation capability due to a quantum Zeno-like effect. Furthermore, in the LSP-chromophore system, engineering the LSP structure, particularly the LSP dissipation rate, can similarly be used to optimize the antioxidation effect. These findings, complemented by analytical results for finding optimal system parameters, provide practical guidelines for designing photostable organic materials with enhanced performance in various optoelectronic applications.

Certifying nonclassical correlations typically requires access to all subsystems, presenting a major challenge in open quantum systems coupled to inaccessible environments. Recent works have shown that, in autonomous pure dephasing scenarios, quantum discord with the environment can be certified from system-only dynamics via the Hamiltonian ensemble formulation. However, this approach leaves open whether stronger correlations, such as entanglement, can be certified. Moreover, its reliance on Fourier analysis requires full-time dynamics, which is experimentally resource-intensive and provides limited information about when such correlations are established during evolution. In this work, we present a method that enables the certification of system-environment quantum entanglement solely from the reduced dynamics of the system. The method is based on the theory of mixed-unitary channels and applies to general non-autonomous pure dephasing scenarios. Crucially, it relaxes the need for full-time dynamics, offering a resource-efficient approach that also reveals the precise timing of entanglement generation. We experimentally validate this method on a Quantinuum trapped-ion quantum processor with a controlled-dephasing model. Finally, we highlight its potential as a tool for certifying gravitationally induced entanglement.

The Liouvillian skin effect and the non-Hermitian skin effect have both been used to explain the localization of eigenmodes near system boundaries, though the former is arguably more accurate in some regimes due to its incorporation of quantum jumps. However, these frameworks predominantly focus on weak Markovian interactions, neglecting the potentially crucial role of memory effects. To address this, we investigate, utilizing the powerful hierarchical equations of motion method, how a non-Markovian environment can modify the Liouvillian skin effect. We demonstrate that a non-Markovian environment can induce a “thick skin effect,” where the skin mode broadens and shifts into the bulk. We further identify that the skin-mode quantum coherence can only be generated when the coupling contains counter-rotating terms, leading to the coherence delocalization and oscillatory relaxation with a characteristic linear scaling with system size. Remarkably, both the skin-mode and steady-state coherence exhibit resistance to decoherence from additional environmental noise. These findings highlight the profound impact of system-bath correlations on relaxation and localization, revealing unique phenomena beyond conventional Markovian approximations.

Exceptional points (EPs) are singularities in the spectra of non-Hermitian operators where eigenvalues and eigenvectors coalesce. Open quantum systems have recently been explored as EP testbeds due to their non-Hermitian nature. However, most studies focus on the Markovian limit, leaving a gap in understanding EPs in the non-Markovian regime. This work addresses this gap by proposing a general framework based on two numerically exact descriptions of non-Markovian dynamics: the pseudomode equation of motion (PMEOM) and the hierarchical equations of motion (HEOM). The PMEOM is particularly useful due to its Lindblad-type structure, aligning with previous studies in the Markovian regime while offering deeper insights into EP identification. This framework incorporates non-Markovian effects through auxiliary degrees of freedom, enabling the discovery of additional or higher-order EPs that are inaccessible in the Markovian regime. We demonstrate the utility of this approach using the spin-boson model and linear bosonic systems. The full description of exceptional points in the non-Markovian regime is complicated by the unclear applicability of standard techniques such as spectral analysis. Here, the authors fill this gap by proposing a general framework that combines the pseudomode equation of motion and the hierarchical equations of motion.

To unequivocally distinguish genuine quantumness from classicality, a widely adopted approach focuses on the negative values of a quasi-distribution representation as compelling evidence of nonclassicality. Prominent examples include the dynamical process nonclassicality characterized by the canonical Hamiltonian ensemble representation (CHER) and the nonclassicality of quantum states characterized by the Wigner function. However, to construct a multivariate joint quasi-distribution function with negative values from experimental data is typically highly cumbersome. Here we propose a computational approach utilizing a deep generative model, processing three marginals, to construct the bivariate joint quasi-distribution functions. We first apply our model to tackle the challenging problem of the CHERs, which lacks universal solutions, rendering the problem ground-truth (GT) deficient. To overcome the GT deficiency of the CHER problem, we design optimal synthetic datasets to train our model. While trained with synthetic data, the physics-informed optimization enables our model to capture the detrimental effect of the thermal fluctuations on nonclassicality, which cannot be obtained from any analytical solutions. This underscores the reliability of our approach. This approach also allows us to predict the Wigner functions subject to thermal noises. Our model predicts the Wigner functions with a prominent accuracy by processing three marginals of probability distributions. Our approach also provides a significant reduction of the experimental efforts of constructing the Wigner functions of quantum states, giving rise to an efficient alternative way to realize the quantum state tomography.

Photon-mediated quantum networks generally consist of quantum channels, repeaters, and end nodes. Remote state preparation (RSP) enables one of the end nodes to prepare the states of the other end nodes remotely. RSP also serves as a deterministic single-photon source for networking communications. Herein, we theoretically and experimentally investigate how networking RSP surpasses any classical emulation without entanglement and qubit unitaries. We introduce a new type of quantum resource, which we refer to as RSP capability, to validate all the static and dynamic elements required for nonclassical state preparation and transmission, such as quantum channels and repeaters. This goes beyond the static resources of quantum correlations. We experimentally demonstrate the RSP capability measurement of the photon pairs created by a polarization Sagnac interferometer, including the transition between classical and nonclassical RSP depending on the photon-pair qualities. Our results help reveal the quantum advantages arising when networking RSP plays a role.

Studies have shown that the Hilbert spaces of non-Hermitian systems require nontrivial metrics. Here, we demonstrate how evolution dimensions, in addition to time, can emerge naturally from a geometric formalism. Specifically, in this formalism, Hamiltonians can be interpreted as a Christoffel symbol-like operators, and the Schroedinger equation as a parallel transport in this formalism. We then derive the evolution equations for the states and metrics along the emergent dimensions and find that the curvature of the Hilbert space bundle for any given closed system is locally flat. Finally, we show that the fidelity susceptibilities and the Berry curvatures of states are related to these emergent parallel transports.