Cristiano Ciuti

We introduce a variational method for simulating the dynamics of interacting open quantum bosonic systems deep in the quantum regime. The method is based on a multidimensional Wigner phase-space representation and employs a variational multi-Gaussian , whose accuracy is systematically controlled by the number of Gaussian components. The variational equations of motion are derived from the Dirac-Frenkel principle and evaluated efficiently by combining the analytical structure of Gaussian functions with automatic differentiation. As a key physical application, we study a driven-dissipative two-dimensional Bose-Hubbard lattice with two-boson coherent driving and two-body losses. Using our dynamical approach, we compute the finite-size scaling of the Liouvillian spectral gap—extracted from the relaxation dynamics—which vanishes in the thermodynamic limit. Our results reveal critical slowing down with dynamical exponents of the 2D quantum Ising universality class, demonstrating the power of our method to capture complex quantum dynamics in large open systems.

We theoretically investigate the role of spatial dimension and driving frequency in a nonequilibrium phase transition of a driven-dissipative interacting bosonic system. In this setting, spatial dimension is dictated by the shape of the external driving field. We consider both homogeneous driving configurations, which correspond to standard integer-dimensional systems, and fractal driving patterns, which give rise to a noninteger Hausdorff dimension for the spatial density. The onset of criticality is characterized by critical slowing down in the excited density dynamics as the system asymptotically approaches the steady state. By analyzing the system-size dependence of the asymptotic decay rate using numerical simulations of the full multimode dynamics, complemented by an analytical statistical mean-field treatment, we determine the lower critical dimension of the nonequilibrium phase transition. We show that this dimension can be noninteger and fractal in nature, and that it can be tuned continuously via the frequency detuning of the driving field.

We theoretically explore the low-energy behavior of a Josephson tunnel junction coupled to a finite-length, charge-biased transmission line and compare it to its flux-biased counterpart. For transmission lines of increasing length, we show that the low-energy charge-dependent energy bands of the charge-biased configuration can be exactly mapped onto those of the flux-biased system via a well-defined duality transformation of circuit parameters. In the limit of an infinitely long transmission line, the influence of boundary conditions vanishes, and both circuits reduce to a resistively shunted Josephson junction. This convergence reveals the system's intrinsic self-duality and critical behavior. Our exact formulation of charge-flux duality provides a foundation for generalizations to more complex superconductor-insulator phase transitions.

Recent studies on cavity-coupled two-dimensional electron gas demonstrate that vacuum-field engineering can tailor electronic transport properties of materials. By achieving ultra-strong coupling between a terahertz resonator and mesoscopic graphene, we demonstrate that cavity vacuum fields can alter the effective degeneracies of Landau levels, resulting in a nonlinear Landau fan diagram for massless Dirac fermions while preserving quantum-Hall quantization. Specifically, by leveraging graphene's gate-tunability, we observe that quantum-Hall features, minimum longitudinal and quantized Hall conductance for a given filling factor, occur at carrier densities reduced by more than 20 percent compared to systems without cavity. Theoretical analysis attributes this effect to the virtual cavity photon mediated transitions between the non-equidistant Landau levels in graphene, significantly reducing their effective degeneracy. This study paves the way for investigating cavity quantum electrodynamics in highly tunable, atomically thin two-dimensional crystals.

We investigate the dynamics of multimode optical systems driven by two-photon processes and subject to nonlocal losses, incorporating quantum noise at the Gaussian level. Our findings show that the statistics retrieved from a single Gaussian quantum trajectory exhibits emergent thermal equilibrium governed by an Ising Hamiltonian, encoded in the dissipative coupling between modes. The system's effective temperature is set by the driving strength relative to the oscillation threshold. Given the ultrashort timescales typical of all-optical devices, our Letter demonstrates that such multimode optical systems can operate as ultrafast Boltzmann samplers, paving the way toward the realization of efficient hardware for combinatorial optimization, with promising applications in machine learning and beyond.