Track the frontier, from condensed matter to cosmology.
Search papers, labs, and topics across Lattice.
To overcome the severe exponential bottleneck and truncation errors inherent in modeling hybrid continuous- and discrete-variable (CV-DV) quantum systems, this work establishes a classical simulation framework based on linear combinations of Gaussian (LCoG) functions. The approach decomposes the hybrid density matrix across DV basis blocks to analytically track bosonic states under controlled Gaussian unitaries and state-dependent noise channels. Crucially, the method bypasses Fock-space truncation entirely, achieving polynomial scaling per Gaussian component in the number of modes and enabling high-fidelity simulations of Gottesman--Kitaev--Preskill (GKP) state preparation in cavity-QED.
Hybrid quantum circuits combining qubits and bosonic modes can now be simulated without Fock-space truncation by mapping controlled operations directly onto linear combinations of Gaussians with polynomial per-mode scaling.
Hybrid quantum systems combining discrete-variable (DV) and continuous-variable (CV) subsystems arise across a broad range of physical platforms and enable quantum information processing beyond purely qubit-based designs. Their classical simulation, however, is challenging since computational costs of both CV and DV systems exponentially increase. In this paper, we introduce a simulation framework for CV--DV hybrid systems based on linear combination of Gaussian (LCoG) decomposition. We decompose the hybrid density matrix into blocks in a DV basis and represent each block as an LCoG expansion. Within the representation, Gaussian operations on the bosonic subsystem independently act on each Gaussian function. We derive analytic formulas describing how complex Gaussian functions transform under controlled Gaussian operations. We further extend the framework to controlled Gaussian channels with state-dependent Gaussian noise. The method avoids a Fock-space truncation for the CV subsystem, and the cost of updating each Gaussian term scales polynomially with the number of modes. The framework provides a general computational tool for analyzing hybrid quantum algorithms, quantum simulation protocols, and non-Gaussian state generation schemes. As an example, we simulate Gottesman--Kitaev--Preskill (GKP) state generation based on a cavity-QED system.