Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIterative solutions to the steady state density matrix for optomechanical systems
106
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a sparse matrix permutation from graph theory that gives stable incomplete Lower-Upper (LU) preconditioners necessary for iterative solutions to the steady state density matrix for quantum optomechanical systems. This reordering is efficient, adding little overhead to the computation, and results in a marked reduction in both memory and runtime requirements compared to other solution methods, with performance gains increasing with system size. Either of these benchmarks can be tuned via the preconditioner accuracy and solution tolerance. This reordering optimizes the condition number of the approximate inverse, and is the only method found to be stable at large Hilbert space dimensions. This allows for steady state solutions to otherwise intractable quantum optomechanical systems.
rate research
Read More
We discuss the numerical solution methods available when solving for the steady-state density matrix of a time-independent open quantum optical system, where the system operators are expressed in a suitable basis representation as sparse matrices. In particular, we focus on the difficulties posed by the non-Hermitian structure of the Lindblad super operator, and the numerical techniques designed to mitigate these pitfalls. In addition, we introduce a doubly iterative inverse-power method that can give reduced memory and runtime requirements in situations where other iterative methods are limited due to poor bandwidth and profile reduction. The relevant methods are demonstrated on several prototypical quantum optical systems where it is found that iterative methods based on iLU factorization using reverse Cuthill-Mckee ordering tend to outperform other solution techniques in terms of both memory consumption and runtime as the size of the underlying Hilbert space increases. For eigenvalue solving, Krylov iterations using the stabilized bi-conjugate gradient method outperform generalized minimal residual methods. In contrast, minimal residual methods work best for solvers based on direct LU decomposition. This work serves as a guide for solving the steady-state density matrix of an arbitrary quantum optical system, and points to several avenues of future research that will extend the applicability of these classical algorithms in absence of a quantum computer.
We study the physical properties of double-cavity optomechanical system in which the mechanical resonator interacts with one of the coupled cavities and another cavity is used as an auxiliary cavity. The model can be expected to achieve the strong optomechanical coupling strength and overcome the optomechanical cavity decay, simultaneously. Through the coherent auxiliary cavity interferences, the steady-state squeezing of mechanical resonator can be generated in highly unresolved sideband regime. The validity of the scheme is assessed by numerical simulation and theoretical analysis of the steady-state variance of the mechanical displacement quadrature. The scheme provides a platform for the mechanical squeezing beyond the resolved sideband limit and addresses the restricted experimental bounds at present.
We study self-oscillations of an optomechanical system, where coherent mechanical oscillations are induced by a driven optical or microwave cavity, for the case of an anharmonic mechanical oscillator potential. A semiclassical analytical model is developed to characterize the limit cycle for large mechanical amplitudes corresponding to a weak nonlinearity. As a result, we predict conditions to achieve subpoissonian phonon statistics in the steady state, indicating classically forbidden behavior. We compare with numerical simulations and find very good agreement. Our model is quite general and can be applied to other physical systems such as trapped ions or superconducting circuits.
Continuing our work on the nature and existence of fluctuation-dissipation relations (FDR) in linear and nonlinear open quantum systems [1-3], here we consider such relations when a linear system is in a nonequilibrium steady state (NESS). With the model of two-oscillators (considered as a short harmonic chain with the two ends) each connected to a thermal bath of different temperatures we find that when the chain is fully relaxed due to interaction with the baths, the relation that connects the noise kernel and the imaginary part of the dissipation kernel of the chain in one bath does not assume the conventional form for the FDR in equilibrium cases. There exists an additional term we call the `bias current that depends on the difference of the baths initial temperatures and the inter-oscillator coupling strength. We further show that this term is related to the steady heat flow between the two baths when the system is in NESS. The ability to know the real-time development of the inter-heat exchange (between the baths and the end-oscillators) and the intra-heat transfer (within the chain) and their dependence on the parameters in the system offers possibilities for quantifiable control and in the design of quantum heat engines or thermal devices.
Variational approaches for the calculation of vibrational wave functions and energies are a natural route to obtain highly accurate results with controllable errors. However, the unfavorable scaling and the resulting high computational cost of standard variational approaches limit their application to small molecules with only few vibrational modes. Here, we demonstrate how the density matrix renormalization group (DMRG) can be exploited to optimize vibrational wave functions (vDMRG) expressed as matrix product states. We study the convergence of these calculations with respect to the size of the local basis of each mode, the number of renormalized block states, and the number of DMRG sweeps required. We demonstrate the high accuracy achieved by vDMRG for small molecules that were intensively studied in the literature. We then proceed to show that the complete fingerprint region of the sarcosyn-glycin dipeptide can be calculated with vDMRG.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
