Papers

In this paper, we examine measurement using the Q-based, objective-field model for quantum mechanics. Schrodinger considered a microscopic system prepared in a superposition of states which is then coupled to a macroscopic meter. We analyze the entangled meter and system, and measurements on it, by solving forward-backward stochastic differential equations for real amplitudes $x(t)$ and $p(t)$ that correspond to the phase-space variables of the Q function of the system at a time $t$. We model the system and meter as single-mode fields, and measurement of $\hat{x}$ by amplification of the amplitude $x(t)$. Our conclusion is that the outcome for the measurement is determined at (or by) the time $t_{m}$, when the coupling to the meter is complete, the meter states being macroscopically distinguishable. There is consistency with macroscopic realism. By evaluating the distribution of the amplitudes $x$ and $p$ postselected on a given outcome of the meter, we show how the $Q$-based model represents a more complete description of quantum mechanics: The variances associated with amplitudes $x$ and $p$ are too narrow to comply with the uncertainty principle, ruling out that the distribution represents a quantum state. We conclude that the collapse of the wavefunction occurs as a two-stage process: First there is an amplification that creates branches of amplitudes $x(t)$ of the meter, associated with distinct eigenvalues. The outcome of measurement is determined by $x(t)$ once amplified, explaining Born's rule. Second, the distribution that determines the final collapse is the state inferred for the system conditioned on the outcome of the meter: information is lost about the meter, in particular, about the complementary variable $p$.

The uniformly accelerated Unruh-DeWitt detector serves as a fundamental model in relativistic quantum metrology. While previous studies have mainly concentrated on single-parameter estimation via quantum Cramér-Rao bound, the multi-parameter case remains significantly underexplored. In this paper, we investigate the multiparameter estimation for a uniformly accelerated Unruh-DeWitt detector coupled to a vacuum scalar field in both bounded and unbounded Minkowski vacuum. Our analysis reveals that quantum Cramér-Rao bound fails to provide a tight error bound for the two-parameter estimation involving the initial phase and weight parameters. For this reason, we numerically compute two tighter error bounds, Holevo Cramér-Rao bound and Nagaoka bound, based on a semidefinite program. Notably, our results demonstrate that Nagaoka bound yields the tightest error bound among all the considered error bounds, consistent with the general hierarchy of multiparameter quantum estimation. In the case with a boundary, we observe the introduction of boundary systematically reduces the values of both Holevo Cramér-Rao bound and Nagaoka bound, indicating an improvement on the attainable estimation precision. These results offer valuable insights on and practical guidance for advancing multiparameter estimation in relativistic context.

Two-dimensional topological insulators protected by nonlocal symmetries or with fragile topology usually do not admit robust in-gap edge modes due to the incompatibility between the symmetry and the boundary. Here, we show that in a parity-time (PT) symmetric system robust in-gap topological edge modes can be stably induced by non-Hermitian couplings that spontaneously break the PT symmetry of the eigenstates. The topological edge modes traverse the imaginary spectral gap between a pair of fragile topological bands, which is opened by the presence of the non-Hermitian perturbation. We demonstrate that the net number of resulting in-gap modes is protected by an operator version of anomaly cancellation that extends beyond the Hermitian limit. The results imply that loss and gain can in principle drive fragile topological phenomena to stable topological phenomena.

Localizable measurements are joint quantum measurements that can be implemented using only non-adaptive local operations and shared entanglement. We provide a protocol-independent characterization of localizable projection-valued measures (PVMs) by exploiting algebraic structures that any such measurement must satisfy. We first show that a rank-1 PVM on $\mathbb{C}^d\otimes\mathbb{C}^d$ containing an element with the maximal Schmidt rank can be localized using entanglement of a Schmidt number at most $d$ if and only if it forms a maximally entangled basis corresponding to a nice unitary error basis. This reveals strong limitations imposed by non-adaptive local operations, in contrast to the adaptive setting where any joint measurement is implementable. We then completely characterize two-qubit rank-1 PVMs that can be localized with two-qubit entanglement, resolving a conjecture of Gisin and Del Santo, and finally extend our characterization to ideal two-qudit measurements, strengthening earlier results.

Aggregates of light harvesting 2 (LH2) complexes form the major exciton-relaying domain in the photosynthetic unit of purple bacteria. Application of a generalized master equation to pairs of the B850 units of LH2 complexes, where excitons predominantly reside, provide quantitative information on how the inter-LH2 exciton transfer depends on the distance, relative rotational angle, and the relative energies of the two LH2s. The distance dependence demonstrates significant enhancement of the rate due to quantum delocalization of excitons, the qualitative nature of which remains robust against the disorder. The angle dependence reflects isotropic nature of exciton transfer, which remains similar for the ensemble of disorder. The variation of the rate on relative excitation energies of LH2 exhibits resonance peaks, which however is fragile as the disorder becomes significant. Overall, the average transfer times between two LH2s are estimated to be in the range of 4 - 25 ps for physically plausible inter-LH2 distances.

In this paper, we revisit the application of classical advantage distillation (CAD) to the decoy-state BB84 protocol. Prior work has shown that CAD can greatly improve maximal distances and noise tolerances of the practical decoy state protocol. However, past work in deriving key-rate bounds for this protocol with CAD have assumed a trivial bound on the quantum entropy, whenever Alice sends a vacuum state in a CAD block (i.e., the entropy of such blocks is taken to be zero). Since such rounds contribute, negatively, to the error correction leakage, this results in a correct, though sub-optimal bound. Here, we derive a new proof of security for CAD applied to the decoy state BB84 protocol, computing a bound on Eve's uncertainty in all possible single and vacuum photon events. We use this to derive a new asymptotic key-rate bound which, we show, outperforms prior work, allowing for increased distances and noise tolerances.

Coherent control, a central concept in physics and chemistry, has sparked significant interest due to its ability to fine-tune interference effects in atoms and individual molecules for applications ranging from light-harvesting complexes to molecular qubits. However, precise characterization of the system's dissipative dynamics is required for its implementation, especially at high temperature. In a quantum control experiment, this means learning system-bath parameters and driving coupling strengths. Here, we demonstrate how to infer key physical parameters of a single molecule driven by spectrally modulated pulses at room temperature. We develop and compare two computational approaches based on two-photon absorption photoluminescence signals: an optimization-based minimization scheme and a feed-forward neural network. The robustness of our approach highlights the importance of reliable parameter estimation in designing effective coherent control protocols. Our results have direct applications in ultrafast spectroscopy, quantum materials and technology.

Random pulse sequences are a powerful method for qubit noise spectroscopy, enabling efficient reconstruction of sparse noise spectra. Here, we advance this method in two complementary directions. First, we extend the method using a regularizer based on the total generalized variation (TGV) norm, in order to reconstruct a larger class of noise spectra, namely piecewise-linear noise spectra, which more realistically model many physical systems. We show through numerical simulations that the new method resolves finer spectral features, while maintaining an order-of-magnitude speedup over conventional approaches to noise spectroscopy. Second, we simplify the experimental implementation of the method, by introducing Rademacher measurements for reconstructing sparse noise spectra. These measurements use pseudorandom pulse sequences that can be generated in real time from a short random seed, reducing experimental complexity without compromising reconstruction accuracy. Together, these developments broaden the reach of random pulse sequences for accurate and efficient noise characterization in realistic quantum systems.

This paper introduces a new logic structure for reversible quantum circuit synthesis. Our synthesis method aims to minimize the quantum cost of reversible quantum circuits with decoders. In this method, multi-valued input, binary output (MVI) functions are utilized as a mathematical concept only, but the circuits are binary. We introduce the new concept of “Multi-Valued Input Fixed Polarity Reed-Muller (MVI-RM)" forms. Our decoder-based circuit uses three logical levels in contrast to commonly-used methods based on Exclusive-or Sum of Products (ESOP) with two levels (AND-XOR expressions), realized by Toffoli gates. In general, the high number of input qubits in the resulting Toffoli gates is a problem that greatly impacts the quantum cost. Using decoders decreases the number of input qubits in these Toffoli gates. We present two practical algorithms for three-level circuit synthesis by finding the MVI-FPRM: products-matching and the newly developed butterfly diagrams. The best MVI-FPRM forms are factorized and reduced to approximate Multi-Valued Input Generalized Reed-Muller (MVI-GRM) forms.

We define a quantum charger as an interacting quantum system that transfers energy between two drives. The key figure of merit characterizing a charger is its charging power. Remarkably, the presence of long-range interactions within the charger can induce a collective steady-state charging mode that depends superlinearly on the size of the charger, exceeding the performance of noninteracting, parallel units. Using the driven Lipkin-Meshkov-Glick model and power-law interacting spin chains, we show that this effect persists up to a critical system size set by the breakdown of the high-frequency regime. We discuss optimal work output as well as experimentally accessible initial states. The superlinear charging effect can be probed in trapped-ion experiments, and positions interacting Floquet systems as promising platforms for enhanced energy conversion.

Quantum device measurements are powerful tools to probe dark matter interactions. Among these, transmon qubits stand out for their ability to suppress external noise while remaining highly sensitive to tiny energy deposits. Ambient galactic halo dark matter interacting with electrons can deposit energy in the qubit, leading to changes in its decoherence time. Recent measurements of transmons have consistently measured, in various experimental setups, a residual contribution to the decoherence time unexplained by thermal noise or known external sources. We use such measurements to set the most stringent laboratory-based constraints to date on dark matter-electron scattering at the keV scale and competitive constraints on dark photon absorption.

We show that every holographic entropy inequality can be recast in the form: "some entanglement wedges reach deeper in the bulk than some other entanglement wedges." When the inequality is saturated, the two sets of wedges reach equally deep. Because bulk depth geometrizes CFT scales, the inequalities enforce and protect the holographic Renormalization Group.

Time-translation symmetry strongly constrains physical dynamics, yet systematic characterization for continuous-variable systems lags behind its discrete-variable counterpart. We close this gap by providing a rigorous classification of Gaussian quantum operations that are covariant under time translations, termed Gaussian covariant operations. We show that several key results known for discrete-variable covariant operations break down in the Gaussian optical setting: discrepancies arise in physical and thermodynamic implementation, in the extensivity of asymmetry, and in catalytic advantages. Our results provide comprehensive mathematical and operational toolkits for Gaussian covariant operations, including a peculiar pair of asymmetry measures that are completely non-extensive. Our findings also reveal surprising consequences of the interplay among symmetry, Gaussianity, and thermodynamic constraints, suggesting that real-world scenarios with multiple constraints have a rich structure not accessible from examining individual constraints separately.

Gravitational time dilation implies that clocks held at different heights accumulate different proper times. We analyze a memory-assisted quantum clock interferometer in which a frequency-bin photonic clock is stored in two vertically separated quantum memories for a controllable duration, such that the joint state evolves in a quantum superposition of two proper times. After retrieval, the photonic modes interfere in a Hong-Ou-Mandel (HOM) interferometer, for which we derive analytic expressions for the resulting multiphoton detection statistics. Extending this HOM-based scheme from entangled photon pairs to frequency-entangled 2N-photon inputs, we show that the proper-time dependent phase is amplified by a factor N, leading to an N-times faster collapse and revival of the interference signal compared with the two-photon case. Incorporating finite memory efficiency and lifetime, we identify regimes where this modulation remains observable. For parameters compatible with demonstrated Rb and Cs memories and achievable optical frequency separations, the first collapse occurs for height differences in the order of 10-100 m with subsecond to few-second storage times, while suitable rare-earth ion and alkali memory combinations can reduce the required height to the few-metre scale. These results establish near-term laboratory conditions for observing entanglement dynamics driven by gravitational time dilation in a photonic platform.

The two-dimensional $O(N)$ nonlinear sigma model (NLSM) is asymptotically free for $N>2$: it exhibits neither a nontrivial fixed point nor spontaneous symmetry-breaking. Here we show that a nontrivial fixed point generically does exist in the complex coupling plane and is described by a complex conformal field theory (CCFT). This CCFT fixed point is generic in the sense that it has a single relevant singlet operator, and is thus expected to arise in any non-Hermitian model with $O(N)$ symmetry upon tuning a single complex parameter. We confirm this prediction numerically by locating the CCFT at $N = 3$ in two non-Hermitian spin-1 antiferromagnetic Heisenberg chains, and in a non-Hermitian spin-$1/2$ ladder, finding good agreement between the complex central charge and scaling dimensions and those obtained by analytic continuation of real fixed points from $N\leq 2$. We further construct a realistic Lindbladian for a spin-1 chain whose no-click dynamics are governed by the non-Hermitian Hamiltonian realizing the CCFT. Since the CCFT vacuum is the eigenstate with the smallest decay rate, the system naturally relaxes under dissipative dynamics toward a CFT state, thus providing a route to preparing long-range entangled states through engineered dissipation.

Grover's algorithm relies on the superposition and interference of quantum mechanics, which is more efficient than classical computing in specific tasks such as searching an unsorted database. Due to the high complexity of quantum mechanics, the correctness of quantum algorithms is difficult to guarantee through traditional simulation methods. By contrast, the fundamental concepts and mathematical structure of Grover's algorithm can be formalized into logical expressions and verified by higher-order logical reasoning. In this paper, we formally model and verify Grover's algorithm in the HOL Light theorem prover. We focus on proving key properties such as the unitarity of its oracle and diffusion operators, the monotonicity of the success probability with respect to the number of iterations, and an exact expression for the optimal iteration count. By analyzing a concrete application to integer factorization, we demonstrate the practicality and prospects of our work.

A nontrivial peculiarity of general relativity is that when the horizon region of black holes is rendered harmless, the exterior doubles, resulting in a causally disconnected parallel universe. This intricacy plays a central role in 't Hooft's unitarity arguments, emphasising an exact identification between the physical universe and its duplicate on the other side of the horizon. However, it leads to another tension in the form of a factor of two correction in Hawking's temperature. This discrepancy is concerning because the Rindler temperature is universal and complies with the Bekenstein-Hawking entropy. We demonstrate that the mismatch in the Boltzmann factor gets fixed if the state that forms the corresponding density matrix adopts a generalised thermofield double structure. That leaves room for some interesting discussion.

Bell inequality tests based on high-dimensional entanglement usually require measurements that can resolve multiple possible outcomes. However, the implementation of high-dimensional multi-outcome measurements is often only emulated via a collection of “click or no-click” measurements. This reduction of multi-outcome measurements to binary-outcome measurements opens a loophole in high-dimensional tests Bell inequalities which can be exploited by local hidden variable models [Tavakoli et al., Phys. Rev. A 111, 042433 (2025)]. Here, we close this loophole by using four-dimensional photonic path-mode entanglement and multi-outcome detection. We test both the well-known Collins-Gisin-Linden-Massar-Popescu inequality and a related Bell inequality tailored for maximally entangled states in high-dimension. We observe violations that are large enough to also rule out any quantum model based on entanglement of lower dimension, thereby demonstrating genuinely high-dimensional nonlocality free of the binarisation loophole.

We have extended the origin-invariant length gauge (LG(OI)) approach – originally developed by Caricato and co-workers for optical rotation (OR) and electronic circular dichroism (ECD) – to vibrational circular dichroism (VCD). This approach avoids the need for gauge-including atomic orbitals (GIAOs), which are typically required to circumvent the unphysical dependence of the CD rotatory strengths on the arbitrary choice of coordinate origin for length gauge (LG) computations. Benchmark VCD spectra are presented for (P)-hydrogen peroxide, (S)-methyloxirane, (1R, 5R)-{\alpha}-pinene, and (1R, 4R)-camphor using Hartree-Fock (HF) theory and density functional theory (DFT) methods across a range of basis sets and compared to those obtained from LG, velocity-gauge (VG), and GIAO computations. These analyses show that for VCD the LG(OI) approach does not converge to the basis-set limit as rapidly as the GIAO approach, but does yield similar quality spectra as GIAO for all major VCD peaks for quadruple-zeta-quality basis sets. The LG(OI) and VG VCD spectra are less reliable compared to GIAOs for smaller basis sets.

Single-photon emitters (SPEs) are central to quantum communication, computing, and metrology, yet their development remains constrained by trade-offs in purity, indistinguishability, and tunability. This review presents a mechanism-based classification of SPEs, offering a physics-oriented framework to clarify the performance limitations of conventional sources, including quantum emitters and nonlinear optical processes. Particular attention is given to hybrid organic-inorganic perovskite quantum dots (HOIP QDs), which provide size- and composition-tunable emission with narrow linewidths and room-temperature operation. Through comparative analysis of physical mechanisms and performance metrics, we show how HOIP QDs may address key limitations of established SPE platforms. Recognizing the constraints of current deterministic sources, we introduce a performance framework to guide the development of scalable SPEs, and examine the theoretical potential of bright squeezed vacuum (BSV) states, discussing how BSV mechanisms could serve as a promising avenue for multiplexable, high-purity photon generation beyond conventional heralded schemes. The review concludes by outlining future directions for integrating HOIP- and BSV-based concepts into scalable quantum photonic architectures.

The theory of symmetry of quantum mechanical systems is applied to study the structure and properties of several classes of relevant maps in quantum information theory: CPTP, PPT and Schwarz maps. First, we develop the general structure that equivariant maps $\Phi:\mathcal A \to \mathcal B$ between $C^\ast$-algebras satisfy. Then, we undertake a systematic study of unital, Hermiticity-preserving maps that are equivariant under natural unitary group actions. Schwarz maps satisfy Kadison's inequality $\Phi(X^\ast X) \geq \Phi(X)^\ast \Phi(X)$ and form an intermediate class between positive and completely positive maps. We completely classify $U(n)$-equivariant on $M_n(\mathbb C)$ and determine those that are completely positive and Schwarz. Partial classifications are then obtained for the weaker $DU(n)$-equivariance (diagonal unitary symmetry) and for tensor-product symmetries $U(n_1) \otimes U(n_2)$. In each case, the parameter regions where $\Phi$ is Schwarz or completely positive are described by explicit algebraic inequalities, and their geometry is illustrated. Finally, we further show that the $U(n)$-equivariant family satisfies $\mathrm{PPT} \iff \mathrm{EB}$, while the $DU(2)$, symmetric $DU(3)$, $U(2) \otimes U(2)$ and $U(2) \otimes U(3)$, families obey the $\mathrm{PPT}^2$ conjecture through a direct symmetry argument. These results reveal how group symmetry controls the structure of non-completely positive maps and provide new concrete examples where the $\mathrm{PPT}^2$ property holds.

We propose a disorder-free one-dimensional single-particle Hamiltonian hosting an exact mobility edge (ME), placing the system outside the assumptions of no-go theorems regarding unbounded potentials. By applying a linear Stark potential selectively to one sublattice of a dimerized chain, we generate an effective Hamiltonian with unbounded, staggered hopping amplitudes. The unbounded nature of the hopping places the model outside the scope of the Simon-Spencer theorem, while the staggered scaling allows it to evade broader constraints on Jacobi matrices. We analytically derive the bulk spectrum in reciprocal space, identifying a sharp ME where the energy magnitude equals the inter-cell hopping strength. This edge separates a continuum of extended states from two distinct localized branches: a standard unbounded Wannier-Stark ladder and an anomalous bounded branch accumulating at the ME. The existence of extended states is supported by finite-size scaling of the inverse participation ratio up to system sizes $L \sim 10^9$. Furthermore, we propose an experimental realization using photonic frequency synthetic dimensions. Our numerical results indicate that the ME is robust against potential experimental imperfections, including frequency detuning errors and photon loss, establishing a practical path for observing MEs in disorder-free systems.

Quantum adiabatic algorithms are commonly analyzed through local spectral properties of an interpolating Hamiltonian, most notably the minimum energy gap. While this perspective captures an important constraint on adiabatic runtimes, it does not fully describe the global structure of spectral evolution in optimization problems with degenerate solution manifolds. In this work, we show that degeneracy alone imposes unavoidable global constraints on spectral flow, even in instances where adiabatic algorithms succeed with high probability. Focusing on digitized quantum adiabatic evolutions, we analyze the eigenphases of the cumulative unitary operator generated along the interpolation path. By explicitly tracking eigenphase trajectories, we demonstrate that multiple spectral bands are forced to interact, braid, and permute before coalescing into a degenerate manifold at the end of the evolution. This global reordering manifests as persistent spectral congestion and nontrivial band permutations that cannot be removed by increasing evolution time or refining the digitization. Using MaxCut instances with controlled degeneracy as a concrete setting, we extract quantitative diagnostics of spectral congestion and explicitly compute the induced band permutations. Our results show that successful adiabatic optimization can coexist with complex and constrained spectral flow, revealing a form of topological obstruction rooted in the global connectivity of eigenstates rather than in local gap closures. These findings highlight intrinsic limitations of gap-based analyses and motivate spectral-flow-based diagnostics for understanding adiabatic algorithms in degenerate optimization landscapes.

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 $\ell_{2}$-type coherence functional.

Modern cyberattacks are increasingly complex, posing significant challenges to classical machine learning methods, particularly when labeled data is limited and feature interactions are highly non-linear. In this study we investigates the potential of hybrid quantum-classical learning to enhance feature representations for intrusion detection and explore possible quantum advantages in cybersecurity analytics. Using the UNSW-NB15 dataset, network traffic is transformed into structured feature vectors through classical preprocessing and normalization. Classical models, including Logistic Regression and Support Vector Machines with linear and RBF kernels, are evaluated on the full dataset to establish baseline performance under large-sample conditions. Simultaneously, a quantum-enhanced pipeline maps classical features into variational quantum circuits via angle encoding and entangling layers, executed on a CPU-based quantum simulator, with resulting quantum embeddings classified using a classical SVM. Experiments show that while classical models achieve higher overall accuracy with large datasets, quantum-enhanced representations demonstrate superior attack recall and improved class separability when data is scarce, suggesting that quantum feature spaces capture complex correlations inaccessible to shallow classical models. These results highlight the potential of quantum embeddings to improve generalization and representation quality in cybersecurity tasks and provide a reproducible framework for evaluating quantum advantages as quantum hardware and simulators continue to advance.

Quantum computation is inherently hybrid, and fast classical manipulation of qubit operators is necessary to ensure scalability in quantum software. We introduce PauliEngine, a high-performance C++ framework that provides efficient primitives for Pauli string multiplication, commutators, symbolic phase tracking, and structural transformations. Built on a binary symplectic representation and optimized bit-wise operations, PauliEngine supports both numerical and symbolic coefficients and is accessible through a Python interface. Runtime benchmarks demonstrate substantial speedups over state-of-the-art implementations. PauliEngine provides a scalable backend for operator-based quantum software tools and simulations.

The interplay between topology and nonlinearity represents a central challenge in modern physics. Here, we investigate this interplay by considering a synthetic Su-Schrieffer-Heeger lattice with all-to-all nonlocal interactions. We find that the distinctive nonlinearity maintains an effective chiral symmetry and leads to a quantized nonlinear winding and Berry phase, as corroborated by the developed Bogoliubov nonlinear adiabatic theory. Increasing nonlinearity drives a sequence of topological transitions signaled by the appearance of characteristic swallowtail band structures at intermediate interaction strengths and band swapping in the strong nonlinear regime. The band swapping results in quantized fractional windings and double-period Bloch oscillations that are closely related to discrete time crystals. Remarkably, even starting from a topologically trivial linear system, nonlocal nonlinearity can induce an emergent topological phase with fractional windings. Experimentally, our model can be realized using photons in a degenerate optical cavity with Rydberg-mediated interactions. Our results establish a rigorous framework and pave the way for exploring nonlinear topological phenomena and their applications in synthetic quantum platforms.

Quantum computers are inherently noisy, and a crucial challenge for achieving large-scale, fault-tolerant quantum computing is to implement quantum error correction. A promising direction that has made rapid recent progress is to design hardware that has a specific noise profile, leading to a significantly higher threshold for noise with certain quantum error correcting codes. This Perspective focuses on erasure qubits, which enable hardware-efficient quantum error correction, by concatenating an inner code built-in to the hardware with an outer code. We focus on implementations of dual-rail encoded erasure qubits using superconducting qubits, giving an overview of recent developments in theory and simulation, and hardware demonstrators. We also discuss the differences between implementations; near-term applications using quantum error detection; and the open problems for developing this approach towards early fault-tolerant quantum computers.

Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.

We present a microscopic analysis and evaluation of the dielectric susceptibility of a dielectric medium consisting of vector-type two-energy-level atoms responding on a weak probe mode when the atoms are driven by a strong coherent field. Each atom, in an environment of others, exists as a quasiparticle further structuring a bulk medium. In a limit of dilute atomic gas, the dynamics of each atom follows the Mollow-type nonlinear excitation regime, and the medium susceptibility collectivizes the individual atomic responses to the probe mode. We outline how the collective dynamics can be interpolated up to a dense medium, and we argue from general positions that in such a medium the optical nonlinearity and, in particular, its parametric part could be significantly magnified by manipulating both the coherent pump and the sample density. That indicates certain limitations for potential capabilities of quantum communication protocols utilizing the entangled photons, created by a parametric process, as a main resource of quantum correlations.