Papers

Hydrogen contamination in Josephson junctions may contribute to device-to-device variability in superconducting qubits, yet its incorporation, location, and bonding in tunnel barriers remain poorly understood. We combine molecular dynamics simulations of aluminum oxidation in dense O2 and H2O environments with atomistic quantum transport calculations. The simulations produce amorphous AlOx barriers with experimentally relevant hydrogen concentrations. Analysis of 400 statistically independent samples shows that the number of incorporated H atoms follows a beta-binomial distribution, reflecting correlations arising from self-limiting oxidation. Most H atoms reside near the AlOx surface, primarily in Al-OH and Al-OH-Al motifs. Quantum transport calculations for Al/Al2O3/Al junctions reveal that hydrogen enhances transmission near the Fermi level by creating an interface dipole that shifts the oxide bands relative to the electrode Fermi level. Combining hydrogen-number statistics with calculated transmission coefficients yields a probability distribution for the Josephson energy. At an average hydrogen content of 2.56 at.%, we predict E_J/h = 10.23 GHz with a standard deviation of 3.6 MHz. Hydrogen shifts the mean Josephson energy by approximately 11%, proportionally to its concentration, but contributes far less to device-to-device variability than observed in state-of-the-art fabricated junctions. These findings provide an atomistic understanding of hydrogen incorporation and its impact on Josephson energy.

We study the protection of information in nearly critical topological quantum codes, constructed by perturbing topological stabilizer codes towards continuous quantum phase transitions. Our focus is on the transverse-field toric code subjected to local Pauli decoherence. Despite the strong quantum fluctuations of anyons when the transverse field is tuned infinitesimally close to the critical point, we show that a finite strength of Pauli decoherence remains necessary to irreversibly destroy information encoded in the ground-state manifold. Using a replica statistical physics mapping for the coherent information, we show that decoherence can be understood as introducing a two-dimensional inter-replica defect within a three-dimensional replica statistical physics model. A field theoretical analysis shows that this defect is perturbatively irrelevant to the bulk critical point, and cannot renormalize the transverse field strength, leading to a finite error threshold. We argue that a qualitatively similar conclusion can be drawn for a broad class of nearly critical topological codes, under a variety of decoherence channels.

Phase estimation with potentially large phase values, i.e., with large dynamic range, has many applications in quantum metrology, for example to atomic clocks. A recently proposed phase estimation scheme approaches the Heisenberg scaling in this global setting using sequences of increasingly squeezed Gaussian states as probes and adaptively chosen, potentially mid-circuit, measurements. In this work, we first observe that the pattern of increase in the squeezing of the probes is applicable even to states with some non-Gaussian features. We then propose an experimentally feasible version of this phase estimation scheme, based on the alternating application of one-axis twist (OAT) operations and rotations. Our protocols are explicitly described in terms of multiple OAT angles whose durations decrease polynomially with system size and spin-squeezing parameters that decay as $N^{-\mu}$, with $\mu>2/3$ in most cases. Using numerical computation of the system-size dependence $N^{-\nu}$ of the Bayesian mean-squared error of an estimator, we show that these states are suitable for use in the phase estimation scheme, and highlight the protocols to achieve $\nu=17/9$ and $53/27$ using two and three OAT operations respectively in the last adaptation stage. We also analyze the limited non-Gaussianity of the resulting probe states and discuss the role of non-Gaussianity in this protocol more generally. Finally, we analyze how robust these protocols are with respect to imperfections such as particle number fluctuations and coherent control fluctuations.

This work combines theory and experiment to examine the mechanisms underlying the harmonic generation of magnons. We develop a nonlinear spin-wave framework that is directly analogous to harmonic generation in nonlinear optics, and combine it with scanning nitrogen-vacancy (NV) center magnetometry to image and quantify magnonic harmonic generation in a Ni$_{81}$Fe$_{19}$/Pt microstripe. Within this framework, the harmonic response arises from nonlinear magnetization dynamics localized at strongly inhomogeneous textures, such as the sample edges and domain walls, that act as anharmonic confining potentials. Scanning probe imaging confirms that the harmonic response is correspondingly nonuniform and concentrated near the sample edges. We measure an expected nonlinear power-law scaling, a systematic shift toward larger wavevector excitations at higher harmonic order, and a spin-selective response indicative of an increasingly chiral harmonic stray field. These results provide a microscopic understanding of magnonic harmonic generation and highlight its potential for engineering nonlinear functionality in magnonic systems.

Topological Quantum Computing has largely evolved towards a paradigm of manipulating edge localized Majorana within $p$-wave topological superconducting nanowires. To bridge the gap between physical qubit systems and quantum algorithms, we perform a dynamical analysis to extend what is known in the adiabatic regime, providing time-dependent gate elements for further qubit and algorithm modeling efforts. Our analysis covers dynamical considerations for two methods of shuttling domain edge bound Majoranas in a single nanowire system which both function by applying spatiotemporally dependent onsite and hopping parameters within the system's Hamiltonian. We then complicate this model by converting it into the T-qubit to calculate the finite-time gate representation of the shuttling techniques used in a more practical setting. These contributions provide insight for realistic experimental setups in the next-generation of qubit implementation and will hopefully facilitate fault tolerant scalable systems and universal gate design.

Large-scale MIMO detection remains challenging because exact or near-maximum-likelihood search is difficult to scale, while available quantum resources are insufficient for directly solving full-size detection instances by QAOA. This paper therefore proposes a Block-QAOA-Aware MIMO Detector (BQA-MD), whose primary purpose is to reorganize the detection chain so that it becomes compatible with limited-qubit local quantum subproblems. Specifically, BQA-MD combines block-QAOA-aware preprocessing in the QR domain, a standards-consistent blockwise 5G NR Gray-HUBO interface, an MMSE-induced dynamic regularized blockwise objective, and K-best candidate propagation. Within this framework, fixed-size block construction gives every local subproblem a uniform circuit width and parameter dimension, which in turn enables parameter-transfer QAOA as a practical realization strategy for structurally matched local subproblems. Experiments are conducted on a 16x16 Rayleigh MIMO system with 16QAM using classical simulation of the quantum subroutine. The results show that the regularized blockwise detector improves upon its unregularized counterpart, validating the adopted blockwise objective and the block-QAOA-aware design rationale. They also show that the parameter-transfer QAOA detector nearly matches the regularized blockwise exhaustive reference and clearly outperforms direct-training QAOA in BER, thereby supporting parameter reuse as the preferred QAOA realization strategy within the proposed framework. In the tested setting, MMSE remains slightly better in the low-SNR region, whereas the parameter-transfer QAOA detector becomes highly competitive from the medium-SNR regime onward.

The Alicki-Lindblad-Fannes dynamical (ALF) entropy measures the rate at which new information is gathered about a quantum system by inspecting its long-time evolution. We propose an extension of the ALF entropy to open quantum dynamics as a measure of back-flow of information from the environment. Such a proposal is stronger than the existing ones based only on the open system reduced dynamics. In the case of a qubit collisionally coupled to a classical spin chain, we obtain an exact expression for the $\textit{open-system ALF entropy}$ explicitly depending on the environment correlations. An extreme case shows how the information flow from environment to system corresponds to vanishing entropy production as for reversible finite quantum systems.

In their 1972 study of approach to equilibrium, Lanford and Robinson showed that gauge-invariant quasi-free states of lattice fermions maximize entropy among all translation-invariant states with a fixed two-point function, and suggested that the maximizer is unique. In subsequent work on this topic, the uniqueness question re-emerged, together with the problem of whether such quasi-free states are weak Gibbs states. We provide a positive answer to both questions within a class of states whose momentum-space two-point function $\widehat C$ satisfies $0<\widehat C(k)<1$ and belongs to the Wiener algebra of the Brillouin zone. The proof reveals that both the entropy maximization principle and weak Gibbsianity follow directly from the thermodynamic formalism for lattice fermions.

Accessing the physical mechanisms behind non-Markovian phenomena in open quantum dynamics requires the study of the statistical properties of the joint system-environment dynamics. This is impossible at the level of the reduced dynamics of the open system alone as the latter is obtained by suitably eliminating the environment. The task is instead made possible by considering multi-time correlation functions involving observables of the open system, only: the open system-environment interactions turn them into global ones thus building up correlations between the two systems. Multi-time correlations form the basis of both the theory of quantum stochastic processes and of the Alicki-Lindblad-Fannes dynamical entropy (ALF entropy for short). This latter quantity provides for quantum systems a measure of the dynamical entropy production as the Kolmogorov-Sinai entropy does for classical systems. In the case of a collisional model whereby the dissipative dynamics of a finite-level system is obtained by its coupling to an infinite classical spin chain, the ALF entropy can be explicitly computed. It turns out to depend on the parameters characterizing the statistical properties of the environment and can be related to the activation and super-activation of memory effects in the open quantum system.

We study low-rank tensor methods for the numerical solution of Schrödinger's equation with time-independent and explicitly time-dependent Hamiltonians, motivated by large-scale simulations of many-body quantum systems and quantum computing devices subject to time-dependent control pulses. We outline the recent application of the "basis update and Galerkin" (BUG) method for tensor trains, and describe the established TDVP and TDVP-2 algorithms based on the time-dependent variational principle. For comparison, we also consider the BUG method in the Tucker format. All these approaches enable memory efficient representations of partially entangled quantum states and thereby mitigate the exponential cost of conventional state-vector formulations. The rank-adaptivity relies on the truncated singular value decomposition, in which the rank of a matrix is reduced by setting its smallest singular values to zero, based on a threshold parameter that controls the truncation error. Numerical experiments on representative time-independent and time-dependent Hamiltonian models quantify the tradeoff between accuracy and compression across methods, with particular attention to the interplay between the time-step and the truncation threshold, and how the computational effort scales with the number of sub-systems in the quantum system.

This work studies the decay of three imaginarity measures, namely the $l_1$-norm-based imaginarity, the robustness of imaginarity, and the relative-entropy-based imaginarity, under typical noise channels. For arbitrary single-qubit pure states subject to real quantum channels, we prove an exact attenuation theorem: the $l_1$-norm-based imaginarity and the robustness exhibit identical attenuation governed by a common channel-dependent factor, whereas the relative-entropy-based imaginarity additionally depends on the input-state orientation relative to the reference basis. We extend the investigation to two-qubit entangled and dual-rail states, establishing a quantitative link between photon-loss probability and residual imaginarity for photonic dual-rail encodings. For separable two-qubit states, we define a real-operation resource preorder and show that $\Omega_{++}=|+\rangle\langle+|\otimes|+\rangle\langle+|$ is operationally maximal. We further distinguish operational de-imaginary power from reference-state imaginarity loss. For several real two-qubit channels we derive exact operational de-imaginary powers, while for other channels we provide analytic reference-state losses serving as benchmarks and rigorous lower bounds.

Two of the most common interpretations of quantum measurement disagree about the fate of quantum amplitudes after measurement, yet this disagreement has not previously led to experimentally distinguishable predictions. In the standard collapse picture, commonly linked to the Copenhagen interpretation of quantum mechanics, measurements eliminate unrealized amplitudes without leaving a memory. In contrast, in the unitary theory, the measurement device registers one of the possible outcomes while remaining part of an entangled state that continues to harbor the unrealized amplitudes. This persistence arises naturally under unitary evolution, since a measurement device that is part of an entangled system cannot serve as a faithful probe of the joint quantum state. Using single-photon measurements of a tunable quantum state, we experimentally show that these two theories make different predictions when three or more consecutive measurements are performed on the same quantum system. Analysis of the joint density matrix of the three measurements reveals coherence among them and supports the unitary theory of quantum measurement. When decoherence is explicitly introduced, the joint density matrix of the quantum system of interest and the apparatus becomes consistent with what a collapse theory would predict. This work clarifies the dynamics of consecutive quantum measurements and offers new insights into the interpretation of quantum measurements.

Near term quantum processors operate in a noise dominated regime, motivating error mitigation techniques that recover accurate expectation values without full fault tolerance. Zero Noise Extrapolation (ZNE) is a widely used but biased error mitigation method that lacks rigorous error bounds. Its effective application requires nontrivial technical choices, most notably the selection of noise scaling factors and extrapolation models, making ZNE sensitive to user expertise and often necessitating costly trial and error procedures. Here, we introduce Folding Free Zero Noise Extrapolation (FF-ZNE), a method that removes the need for noise factor selection by achieving effective noise amplification without circuit folding. FF-ZNE exploits isomorphic hardware layouts with distinct native noise profiles, such that executing a fixed circuit across these layouts induces controllable variations in the effective noise strength. Under a depolarizing noise model, we analytically show that the resulting extrapolation admits a fixed linear form, eliminating extrapolator choice and enabling a seamless, user independent mitigation procedure. We further propose two algorithms that identify sets of isomorphic hardware layouts on which a given circuit yields sufficiently distinct expectation values to enable reliable zero noise extrapolation. Experiments on a 133 qubit IBM Quantum device demonstrate that FF-ZNE yields mitigated expectation values with average deviations of ~6% and 4.5% for up to 50 qubit EfficientSU2 (sparse) and Hamiltonian simulation (dense) circuits, respectively. The method is thus scalable and applicable to a broad range of circuits. By eliminating noise factor and extrapolator selection, FF-ZNE transforms zero noise extrapolation from a technique requiring expert tuning into a practical, scalable, and broadly accessible error mitigation method for current quantum hardware.

We propose collective nuclear polaritons formed by hybridizing a 229Th nuclear ensemble with a vacuum-ultraviolet cavity mode generated via four-wave mixing, achieving a collective light-matter coupling that scales as $\sqrt{N}$. In the strong-coupling regime the system displays vacuum Rabi oscillations, indicating the hybridization between cavity photons and nuclear excitations. In the superradiant regime, the stored excitation is released in a cooperative burst with peak intensity scaling as $N^2$. The emission lifetime shrinks from thousands of seconds to the millisecond scale and remains tunable. Detuning sweeps across the polariton avoided crossing allow adiabatic conversion of the photonic excitation into a collective nuclear excitation, enabling reversible quantum storage. Our results demonstrate that cavity-mediated nuclear polaritons enable deterministic lifetime engineering and coherent quantum storage in nuclear systems.

A quantum channel is defined as being completely positive (CP) and trace preserving (TP). While not every quantum channel is invertible or reversible, every quantum channel admits various kinds of generalized inverses such as the Moore-Penrose inverse and the Drazin inverse. A generalized inverse of a quantum channel may not itself be a quantum channel: it often fails to be CP. However, generalized inverses still play an important role in quantum error mitigation. Here, because it is often desirable for the generalized inverse of a quantum channel to be at least TP, the Drazin inverse, which is TP, is favoured over the Moore-Penrose inverse, which is not in general TP. In this paper, we take a categorical perspective on generalized inverses of quantum channels. This allows us to give a simple proof of the fact that the Drazin inverse of a quantum channel is always TP. It also allows us to show that for unital quantum channels, the Drazin inverse is also unital. We then generalize this result to dagger Drazin inverses, which allows us to show that for unital quantum channels, the Moore-Penrose inverse is both TP and unital as well. This opens the door to new applications of both the Drazin inverse and Moore-Penrose inverse in quantum information theory and, in particular, in quantum error mitigation.

We present a perturbative closed-time-path (in-in) formulation of an optomechanical system in which a quantum field interacts with a moving mirror via radiation pressure. We derive the effective action governing the dynamics of the moving mirror, incorporating the full back-reaction of the cavity field. These effects are encoded in fluctuation and dissipation kernels, that we show satisfy fluctuation-dissipation relations, and whose spectral structure reveals a direct connection with the underlying physical mechanism responsible for the back-reaction, that is particle creation by the dynamical Casimir effect. By deriving the semiclassical equations of motion for the moving mirror, and computing the energy radiated into the field within the in-out formalism of quantum field theory, we verify the energy balance between the mechanical energy dissipated by the optical back-reaction forces acting on the mirror and the energy carried by the particles created in the field.

Quantum criticality has received extensive attention due to its ability to significantly enhance quantum sensing. But its realization and control in many-body quantum systems remain challenging. We present an effective scheme to simulate the Lipkin-Meshkov-Glick (LMG) model using a coherent Ising machine (CIM) composed of a network of degenerate optical parametric oscillators (DOPO). In our work, the spin variables of the LMG model are mapped onto the phases of DOPO pulses, and the spin-spin interactions are realized by all-to-all couplings among them. Through our investigation of the critical behavior in the antiferromagnetically coupled LMG model in the thermodynamic limit, i.e., $N\rightarrow\infty$, and its application in quantum sensing near the critical point, we verify that the CIM does not only effectively capture the second-order quantum phase transition (QPT) at the critical point but also reconstructs its complete phase diagram under ferromagnetic coupling. Furthermore, we demonstrate how the critical dynamics of this simulation platform can be utilized for quantum-enhanced metrology, achieving a measurement precision that diverges near the critical point of the LMG model. These results highlight the capability of the CIM as a flexible experimental platform for investigating the QPT in the fundamental quantum magnetic models, providing valuable insights into quantum simulation and critical phenomena.

The Kirkwood-Dirac (KD) quasiprobability distribution is a fundamental representation for quantum states and has been widely applied in quantum metrology, quantum chaos, weak values in recent years. A quantum state is KD-classical if its KD-quasiprobability distribution forms a valid classical probability distribution with respect to two given bases, and KD-nonclassical otherwise, with the latter being closely associated with quantum advantages in various quantum processes. In this work, we investigate the structural characteristics of the KD-classical state set when the transition matrix between two orthonormal bases takes the form of a discrete Fourier transform (DFT) matrix. First, we adopt an alternative analytical approach to prove that the set of KD-classical states in a $p^r$-dimensional Hilbert space is the convex hull of KD-classical pure states–a conclusion that was recently established by De Bièvre et al [Annales Henri Poincaré, 1-20, 2025]. Furthermore, we define a directed graph and use it to characterize KD-classical pure states in a Hilbert space of arbitrary dimension $d$. That is, the convex hull of KD-classical pure states along any path from the start vertex to the end vertex in this directed graph is exactly the intersection of the KD-classical state set and the linear space spanned by these path-associated KD-classical pure states. This general result not only yields the $p^r$-dimensional conclusion in a straightforward manner but also encompasses Theorem 2 in the existing work [J. Phys. A, 57, 435303, 2024], demonstrating its generality and inclusiveness.

Measurement backaction degrades dispersive readout of superconducting qubits even at modest drive strengths, often via the reduction of qubit lifetimes during readout. In this work, we theoretically and experimentally study this degradation and show how it can result from the interplay between detuned two-level systems (TLSs) and a drive-renormalized qubit spectrum. For modest to strong readout, the qubit emission spectrum becomes non-Lorentzian and depends sensitively on the readout drive frequency (even when measurement rate is fixed). We combine the readout-modified qubit emission spectrum with time-dependent perturbation theory to predict qubit lifetimes in the presence of a TLS bath. Master equation simulations and experimental measurements on a frequency-tunable transmon confirm these predictions quantitatively. In particular, we find that driving at the resonator frequency associated with the qubit ground state yields the narrowest qubit emission spectrum and the least lifetime degradation for a fixed measurement rate, providing a practical guideline for optimizing readout protocols in future quantum processors.

Millimeter waves are emerging as an enabling technology for connecting and enhancing different quantum platforms such as Rydberg atoms, optomechanics, and superconducting qubits. In this work, we explore the interaction between millimeter wave photons and conventional transmon qubits, specifically for qubit readout. We study a circuit quantum electrodynamic (cQED) system consisting of a millimeter wave cavity at $\omega_r = 2\pi \times 34.7$ GHz and a transmon qubit at $\omega_q = 2\pi \times 3.1$ GHz coupled at rate $g = 2\pi \times 1.3$. With such a large detuning between cavity and qubit, $\omega_r/\omega_q > 10$, we are able to suppress drive induced unwanted resonant state transitions, enabling strong drives for qubit readout. We measure no resonant state transitions even when driving with more than $100$ photons and read out the qubit state with more than $100$ photons to achieve a measurement fidelity greater than $99\%$ without the aid of a quantum limited amplifier.

For high-bias operation devices such as silicon carbide (SiC) power devices, early detection of failure mechanisms is essential to ensure reliability. This requires a method to map high electric fields with high spatial resolution, which has not been realized until now. Here we report that the silicon vacancy (Vsi) in SiC has outstanding characteristics for detecting electric fields applied in various directions within a high-biased SiC device. Vsi exhibits an equivalent response to electric field components parallel (Epara) and perpendicular (Eperp) to the c-axis, a feature unique among quantum sensors, and the responsiveness to Epara and Eperp enables detection of arbitrary electric fields encountered in cutting-edge SiC power devices. We confirmed high electric field detection of ~2.3 MV/cm, which is ~90% of the breakdown electric field of a 4H-SiC with typical carrier concentration. Selectively formed Vsi enables high-resolution mapping of electric field distribution. Vsi-based quantum sensors bring data-driven research and development methodologies as well as device degradation diagnosis.

Rydberg-atom electric field sensing has shown great potential from near-DC to THz with state-of-the-art measurement metrics realized in sensitivity, phase extraction, multi-band receptivity, etc. While Rydberg-atom sensors have shown exceptional performance in the GHz regime, low-frequency operation has remained challenging because of electric-field-screening in conventional vapor cells, which suppresses externally applied fields. We overcome this limitation by combining auxiliary modulation and lock-in detection with a paraffin-coated vapor cell, and demonstrate an electrode-free, wideband method for sensing frequencies, ranging from 0.5 Hz to 10 kHz. Our work extends Rydberg-atom sensor range to VLF, ULF, SLF, ELF and sub-ELF frequency bands. In our method, high state-of-the-art sensitivities have been achieved - 819 $\mu$V/cm/$\sqrt{\text{Hz}}$ for 1 Hz, 33 $\mu$V/cm/$\sqrt{\text{Hz}}$ for 10 Hz, 10 $\mu$V/cm/$\sqrt{\text{Hz}}$ for 100 Hz and 2 $\mu$V/cm/$\sqrt{\text{Hz}}$ for 1 kHz.

We investigate the Magnus expansion for a generic time-dependent two-level system under single-axis driving.By virtue of the su(2) Lie algebra, the expansion is decomposed into a commutator-free form. To illustrate the usefulness of the gained expression, we then revisit the Landau-Zener-Stückelberg-Majorana model, with a focus on non-adiabatic transitions as well as the Stokes phase. In addition, the semiclassical Rabi model is systematically treated by determining the Floquet quasienergy up to different orders. We demonstrate how to employ suitable picture transformations as well as on how to enforce the symmetry of the underlying model in order to guarantee convergence of the expansion as well as to achieve satisfactory agreement with the exact results. For both models that we studied it turns out that a third order approximation yields results that are in next to perfect agreement with exact analytical ones. Surprisingly, in the case of the semiclassical Rabi model, even the second order Magnus approximation in the adiabatic picture produces almost exact results for a large parameter range.

In this study, we developed a diamond quantum magnetometer based on Ramsey interferometry with a short sensor-to-sample distance. Conventional biomagnetic sensors with ensemble nitrogen-vacancy centers using continuous-wave optically detected magnetic resonance and Ramsey methods typically rely on watt-level lasers to achieve high sensitivity, resulting in thermal issues. In contrast, by employing the light-trapping diamond waveguide technique in a high-pressure and high-temperature diamond sample treated with electron beam irradiation, we obtained a high photon conversion efficiency of 9.5%, enabling us to simultaneously achieve a high sensitivity of 2.93(7) pT/Hz^1/2 in the 100-400 Hz frequency range and a minimal temperature increase of only approximately 13 K at a low laser power of 210 mW. Using a dry phantom designed to mimic magnetoencephalography signals, we measured a weak magnetic field of 77.7(2) pT without signal averaging at a sensor-to-sample distance of 2.5 mm. This short-distance measurement prevents severe spatial signal attenuation, yielding a high signal-to-noise ratio. The development here is crucial for practical biomagnetic applications based on Ramsey interferometry.

In measurement-based quantum computation (MBQC), quantum circuits are implemented using adaptive measurements on an entangled resource state. In practice, the resource state will always be prepared with some noise, and it is crucial to understand the effect of this noise on the operation of MBQC. Typically, one measures the fidelity of the noisy resource state with the assumption that a high fidelity state means a high fidelity computation. However, the precise relationship between these two fidelities is not known. Here, we derive an expression that equates the average fidelity of the MBQC output state to a certain correlation function evaluated on the noisy resource state. Using this expression, we show that state fidelity provides a tight lower bound on average MBQC fidelity. Conversely, we also find that state fidelity can greatly underestimate average MBQC fidelity, implying that state fidelity is not a good indicator of MBQC performance in general. In response, we formulate an efficient method to directly estimate average MBQC fidelity by measuring the aforementioned correlation function. These results therefore improve our ability to characterize noisy resource states in quantum computers and benchmark MBQC performance.

Scaling solid-state architectures to the millions of qubits required for utility-scale quantum computing could benefit from the integration of control electronics in the immediate vicinity of the quantum layer. However, lithographically fabricated solid-state qubits perform best at temperatures well below 1 K, where available cooling power is limited, whereas the control electronics dissipate substantial power and therefore require the higher cooling power available at elevated temperatures. To address this challenge, we propose a cryopackaging concept that uses broadband phononic Distributed Bragg Reflectors (DBRs) as a thermal barrier between cryoelectronics and the qubit chip. As an experimental realization of this concept, we fabricate and characterize Ta/SiO$_2$ DBR structures. In this architecture, the DBR is intended to provide mechanical support for superconducting vias while offering substantially better thermal insulation than typical bulk materials. For a 600-nm-thick DBR consisting of 10 Ta/SiO$_2$ bilayers, we obtain a thermal conduction below 1 mW/cm$^2$ from 1.5 K to 100 mK. In a centimeter-scale architecture, this level of isolation is compatible with Watt-level cooling power for nearby electronics while maintaining a qubit temperature around 100 mK in commercially available dilution refrigerators.

Quantum channels generally reduce the distinguishability of quantum states, thereby constraining information transmission and processing in open quantum systems. While it is known that distinguishability can be partially recovered through suitable post-processing protocols, a systematic characterization of the maximal achievable gain has remained elusive. Here, we establish a general framework to determine and optimize the recovery of distinguishability induced by a quantum channel. We introduce an algorithm that identifies the optimal implementation of a multi-copy coarse-graining protocol applicable to arbitrary channels. Within this framework, we derive a general upper bound on the attainable distinguishability gain and quantify the performance of the protocol through its tightness relative to this bound. Our results show that collective coarse-graining can lead to effective dynamics exhibiting the operational signatures of essential non-Markovianity even when the underlying microscopic dynamics remain weakly non-Markovian. A detailed analysis reveals a nontrivial trade-off between mathematical optimality (bound saturation) and operational optimality, together with a strong dependence on both the input ensemble and the number of copies. Taken together, these findings provide a unified and quantitative framework to assess, optimize, and interpret distinguishability recovery in open quantum systems while showing how collective processing can qualitatively modify the operational memory properties of effective dynamical descriptions.

Quantum systems inside high-Q cavities offer an excellent testbed for the control of emergent symmetries induced by light and their interplay with quantum matter. Recently several developments in cavity experiments with neutral atoms and other quantum objects such as ions motivate the study of their quantum correlated properties and their entanglement to tailor and control the behavior of the system. Using the enhanced coupling between light and interacting matter we explore the properties of emergent superradiant modes using our newly developed Light-Matter DMRG algorithm with strongly interacting spin chains. We explore a experimentally viable generalization of the transverse Ising chain coupled to the cavity light where it is possible to induce multimode structures tailored by the light pumped into the system. We find a plethora of scenarios can be explored with clear and accesible measurable signatures. This allows to study the physics of emergent orders and strong quantum correlations with quantum spins where the local and long range coupling can be efficiently simulated. We find that quantum spin nematic states with long range order and magnon pairs emerge as the transitions to superradiant phases take place. Notably, we show the cavity field allows the optimization of entanglement between spins for different light induced modes which can be used for quantum state engineering of quantum correlated states. Our methods can be used to model other hybrid quantum systems efficiently.

This paper addresses the Bi-Objective Traveling Thief Problem (BI-TTP), a challenging multi-objective optimization problem that requires the simultaneous optimization of travel cost and item profit. Conventional methods for the BI-TTP often face severe scalability issues due to the complex interdependence between routing and packing decisions, as well as the inherent complexity and large problem size. These difficulties render classical computing approaches increasingly inapplicable. To tackle this, we propose an advanced hybrid approach that combines quantum annealing (QA) with the $\varepsilon$-constraint method. Specifically, we reformulate the bi-objective problem into a single-objective formulation by restricting the second objective through adjustable $\varepsilon$-levels, determined within established upper and lower bounds. The resulting subproblem involves a sum of fractional terms, which is reformulated with auxiliary variables into an equivalent form. Subsequently, the equivalent formulation is transformed into a Quadratic Unconstrained Binary Optimization (QUBO) model, enabling direct solution via a quantum annealing (QA) solver. The solutions obtained from the quantum annealer are subsequently refined using a tailored heuristic procedure to further enhance overall performance. By leveraging the flexibility in selecting $\varepsilon$ parameters, our approach effectively captures a broad Pareto front, enhancing solution diversity. Experimental results on benchmark instances demonstrate that the proposed method effectively balances two objectives and outperforms baseline approaches in time efficiency.

Density functional theory (DFT) is a cornerstone of computational chemistry and materials science, but its computational cost limits its use in large-scale and high-throughput applications. While machine learning has accelerated energy prediction for specific molecular classes, transferable prediction of electron density across diverse chemical spaces remains challenging. Here, we present a universal framework based on Fourier Neural Operators (FNOs) that directly learns the mapping from external potentials to electron density distributions. Unlike conventional approaches that rely on explicit atomic orbitals, basis sets, or handcrafted descriptors, the proposed method captures global electronic interactions and long-range correlations through operator learning in the spatial-frequency domain. Trained on datasets spanning multiple elements and molecular geometries, the model achieves zero-shot generalization to entirely unseen molecular systems and accurately predicts their electron densities without retraining. This transferability arises from the intrinsic ability of FNOs to represent global structure in continuous fields. Our work establishes neural operator learning as a promising route for fast, accurate, and transferable electronic structure prediction, with potential applications in high-throughput screening and chemical space exploration.