Quantum Approach to Damped Three Coupled Nano-Optomechanical Oscillators

We investigate quantum features of three coupled dissipative nano-optomechanical oscillators. 'e Hamiltonian of the system is somewhat complicated due not only to the coupling of the optomechanical oscillators but to the dissipation in the system as well. In order to simplify the problem, a spatial unitary transformation approach and a matrix-diagonalization method are used. From such procedures, the Hamiltonian is eventually diagonalized. In other words, the complicated original Hamiltonian is transformed to a simple one which is associated to three independent simple harmonic oscillators. By utilizing such a simplification of the Hamiltonian, complete solutions (wave functions) of the Schrödinger equation for the optomechanical system are obtained. We confirm that the probability density converges to the origin of the coordinate in a symmetric manner as the optomechanical energy dissipates. 'e wave functions that we have derived can be used as a basic tool for evaluating diverse quantum consequences of the system, such as quadrature fluctuations, entanglement entropy, energy evolution, transition probability, and the Wigner function.


Introduction
Physical systems in nature do not behave independently in most cases because they are not isolated in usual. e coupling of a system to another one often results in various mutual phenomena, such as energy exchange, dissipation, entanglement, amplitude fluctuations, and decoherence [1,2]. e model of a chain of oscillatory motions can be utilized in analyzing the dynamical characteristics of coupled optomechanical [3][4][5][6], nanoelectromechanical [7,8], chemical [9,10], and biological systems [11]. e mechanical analyses of coupled oscillators have been extensively explored so far through different approaches. Such systems can be usually investigated using invariant operator methods [12,13], Bogoliubov transformation methods [14,15], Jaynes-Cummings approaches [16], path integral methods [17], and adiabatic approaches [2]. Based on these approaches, we can elucidate diverse mechanical properties of the coupled oscillatory systems even when they are described by a time-dependent Hamiltonian.
Coupled harmonic oscillators are important topics especially in optomechanical systems which utilize interaction and entanglement between optical and mechanical modes. Coupling of nano-and micromechanical oscillators with photons [3] and/or other systems such as electrons [18] and atoms [19] provides a substantial tool for realizing nextgeneration quantum technologies. For instance, photonic couplings provide a basic technology platform for manipulation of phononic structures [20], slow/fast-light technology [21], zero-point cooling [22], and mechanicalfrequency shifting [23]. Besides such applicability, coupled optomechanical oscillatory systems are required for highly sensitive measuring quantum devices which are crucial in quantum state tomographies and information processing with quantum states. Moreover, abundant physical phenomena along this line can be utilized in nonclassical state preparation in hybrid quantum-information systems including quantum networks with appropriate quantum protocols.
In order to know underlying quantum properties related to entanglement [24,25], quantum coherence and decoherence [26], nonclassical correlations [27], and concurrence [28], exact quantum theory which explains the evolution of the coupled system is necessary. It is known that such a system can be quantized with the assistance of a dynamical invariant [29] or imposing an adiabatic condition [2]. In this paper, we will show that we can also unfold the quantum dynamics of optomechanical physical systems described by a Hamiltonian of dissipative three coupled oscillators through neither introducing a dynamical invariant nor imposing an adiabatic condition in the dissipation. e complete wave functions of the system will be derived by diagonalizing the Hamiltonian directly, instead of the diagonalization of the dynamical invariant.
Before we carry out the main diagonalization, we first transform the Hamiltonian by a unitary operator into an appropriate form that can be diagonalizable without any approximation and condition. Indeed, unitary transformation procedure is a useful tool that enables us to treat a complicated Hamiltonian in a simplified transformed space.
is greatly helps us to solve the analytical quantum solutions of the mechanical systems. e present paper is organized as follows. Our problem for dissipative nano-optomechanical three coupled oscillators is formulated by basic equations of mechanics in Section 2. In Section 3, we deal with the quantization of the system and its physical states. At first, we simplify the Hamiltonian of the system by using the unitary transformation method. Such a simplified Hamiltonian will be represented as a rigorous matrix form for further investigation. e matrix Hamiltonian will be diagonalized so that we can easily manage it. By taking advantage of such a diagonalization, the quantum solutions of the system will be investigated. e concluding remarks are given in the last section.

Preliminaries
Before the actual investigation of optomechanical oscillatory systems, let us briefly show the method for describing mechanical oscillators subjected to a dissipational nonconservative force. In classical mechanics, the motion of a damped mechanical system can be described by the Newtonian equation e second term on the left-hand side is a dissipative frictional force proportional to velocity, where δ is the damping coefficient. From the theoretical point of view, there is no limitation on the magnitude of δ; i.e., it can be any real number. For a mechanical oscillatory system that we are interested, the term in the right-hand side is of the form (2) Whereas the system is underdamped provided that δ < ω/2 where ω � (k/m) 1/2 , it is well known that the system becomes an overdamped oscillator when δ > ω/2. For the case of the overdamped oscillator, its mathematical treatment is somewhat difficult from quantum mechanical point of view. e Hamiltonian that gives the equation of motion, equation (1) with equation (2), is given by e role of this Hamiltonian is limited to a generator of the classical equation of motion. Regarding this, it is well known that the energy of the system in this case is another problem and it is given by [30][31][32] where 〈· · ·〉 is an expectation value in the Fock state. While equation (4) is represented in terms of the canonical momentum, we can also represent it in terms of the physical momentum P k using the relation that P k � m _ x � Pe −δt : is implies that the quantum energy of the system dissipates according to the decrease of both 〈P 2 k 〉 and 〈x 2 〉 over time.
is basic description of 1D system can be easily extended to coupled oscillatory systems.
In this work, we treat an optomechanical system shown in Figure 1 as a generalization of 1D system described by equation (3). As can be seen from Figure 1, optical fields are coupled with a mechanical mode through a nonlinear interaction caused by radiation pressure. e system can be linearized to some extent [3] and the resultant Hamiltonian can be parameterized by the optical quadratures (x 1 , x 2 ) and a position of the mechanical membrane (x 3 ) [29, [33][34][35]: where m 1 and m 2 are the effective electric permittivities in the cavity 1 and cavity 2, respectively, in Figure 1 (although an electric permittivity is written as ε in many cases, we denote it as m for the convenience of mathematical expressions), m 3 is the effective mass of the nanomembrane, k i (i � 1, 2, 3) are the stiffness constants, and k 12 , k 13 , and k 23 are the coupling constants; the convention of the subscript i (including j) which is given here will also be used subsequently. We manage only the case of underdamped oscillation of the system for simplicity. For a simple case where the system is a coupling of two oscillators (instead of three) with the condition δ ⟶ 0, the system corresponds to that of [3] and the relevant Hamiltonian, equation (6), reduces to equation (1) in the same reference.
In phase space, the coordinates and the conjugate momenta satisfy the commutation relations We are interested in quantum mechanical treatment of the optomechanical system. To study the quantum behavior of our system, we consider the Schrödinger equation of the form e evolution of the state vector associated with this equation and its physical interpretation is usually made within the framework of the Copenhagen interpretation of quantum mechanics. However, because the Hamiltonian in equation (6) is a somewhat complicated form, it may be not an easy task to solve the solutions of equation (8). In order to overcome such a difficulty, we will manage the system based on the unitary transformation method and a diagonalization procedure in the next section.

Unitary Transformations.
Unitary transformation methods have often been used to solve the time-dependent Schrödinger equation (TDSE) [36][37][38][39][40][41][42] for a complicated dynamical system. If we use this method, it is possible to reduce the complicated form of the Schrödinger equation to a simpler one. is may facilitate the derivation of associated quantum solutions.
In order to achieve an appropriate transformation, a suitable choice of a unitary operator depending on the Hamiltonian is necessary. In our case, we choose the unitary operator as [29,36,40,41] is is a product of two unitary operators. e first operator is represented by considering the symmetry between x i and P i . Later on, we will see that the choice of this operator simplifies the Hamiltonian through a transformation. In terms of U(t), the transformation of the wave function can be carried out in a way that where ψ(t) is the transformed wave function. Under this unitary transformation, the Schrödinger equation of the original Hamiltonian system, equation (8), is mapped into where the new Hamiltonian H has the form e transformation of canonical variables using equation From the use of equation (6) with the above relations, we see that equation (12) becomes Figure 1: Schematic of three coupled nano-optomechanical oscillators that we consider. ω 1 and ω 2 are the optical modes, whereas ω 3 is a mechanical mode. It represents the interaction of the cavity fields with a flexible nanomembrane via radiation-pressure forces. Complexity where We see from equation (14) that we have taken the Hamiltonian which is simplified to some extent, thanks to the transformation with the use of U(t) given in equation (9). However, the coupling terms x i x j still remained in the Hamiltonian. We will further simplify the Hamiltonian in the subsequent subsection by removing the terms x i x j using a diagonalization method.

Diagonalization of the Hamiltonian.
Let us now consider diagonalization of the Hamiltonian. To this end, we represent the Hamiltonian H in a matrix form such that where Γ ij are elements of the matrix which corresponds to ith row and jth column for each. If we think that the dimension of the matrix Γ is (3 × 3), Γ is a diagonalizable square matrix [43,44]. e diagonalizability of a matrix is in principle equivalent to the existence of a basis of eigenvectors, which makes it possible to define a diagonalizable endomorphism of a vector space. To diagonalize the Hamiltonian matrix, equation (16), it is necessary to seek its eigenvalues ϖ 2 i and the corresponding eigenvectors V → i . In our case, by solving the secular equation for Γ given in equation (17), the eigenvalues of the matrix Γ are given as follows: where From a standard procedure associated with the eigenvalue problem of a matrix, we can easily see that the corresponding normalized eigenvectors are given by where Note that the vectors V → i (i � 1; ±) are orthonormal to each other because Γ is symmetric. us, Γ can be diagonalized by Γ � RDR −1 , where 4 Complexity Now, we introduce new coordinates q i and p i as By using equation (22), we can easily obtain the full expressions of q i and p i ; we have represented them in Appendix for completeness. Now, from a straightforward algebraic procedure using those expressions (Equation (A.1)-(A.6)), we see that H takes the form us, the transformed system is eventually decoupled (diagonalized); i.e., it is reduced to a collection of three decoupled harmonic oscillators. During the diagonalization procedure, we have imposed no conditions or restrictions. Hence, the consequence (equation (25)) justifies the appropriateness of the chosen unitary operator together with the subsequently performed diagonalization. Based on this, the solutions of our problem can be easily obtained. Indeed, from a mathematical point of view, the diagonalization method plays a major role in the simplification of the problem without modification of the behavior of the physical system.
We note that the coupled optomechanical oscillatory systems can also be diagonalized by using the Bogoliubov transformation or its advanced technics, provided that δ � 0. However, the resulting expression of the Hamiltonian which is diagonalized based on such a method is somewhat different from that obtained using the method used here (see, for example, equations (54-57) in Ref. [15]).
Notice that we have used no approximation and perturbation methods during the derivation of equation (28) with equations (29)- (36). Hence, the wave functions in equation (28) are exact even if they are somewhat complicated. e normalization factor, equation (29), is chosen in a way that Equation (29) involves e 3δt/4 , and this reduces to e δt/4 for a 1D case, which is a previously known consequence [30]. e obtained wave functions can be used as a basic tool for further investigation of the quantum dynamics of the nano-optomechanical system. As mentioned earlier, we considered only the case of underdamped oscillation. For the cases of critically and overdamped systems, the quantum mechanism of the oscillation is somewhat different. en, the wave functions are not described by quantum number n i but are similar to those of an unbound system [30].
We have plotted the time behavior of the probability density |Ψ n 1 ,n 2 ,n 3 (x 1 , x 2 , x 3 , t)| 2 for several different choices of quantum numbers n i and the damping factor δ in Figure 2. e numerical data associated to the graphics in this figure were obtained from the use of Mathematica program (Wolfram Research) with equations (28)- (36). e probability density converges to the origin of the coordinate in a symmetric manner as time goes by according to the dissipation of energy, regardless of the chosen values of n i . e convergence is fast when δ is large as expected.

Conclusions
We have investigated quantum mechanical features of dissipative three coupled nano-optomechanical oscillators. Exact solutions of the TDSE of the system have been derived using the unitary transformation method and a diagonalization procedure. From the unitary transformation at first, the complicated time-dependent original quantal Hamiltonian was transformed to a simple one which has unit masses; however, the coupling terms x i x j in the Hamiltonian still remained even after this transformation. In order to remove the coupling terms x i x j , we used a diagonalization method. rough these procedures, the Hamiltonian eventually diagonalized. e finally transformed Hamiltonian was given in the form associated to three independent harmonic oscillators. As a consequence, the transformed Hamiltonian was easily treated and we have identified the corresponding quantum solutions without difficulty. By inverse transformation of the solutions of the TDSE for the transformed Hamiltonian, we finally had the complete quantum solutions of the original optomechanical system. e probability density, which is the absolute square of the wave function, converges to the origin of the coordinate as the oscillatory energy dissipates. e wave functions, which we obtained here, can be used to evaluate not only the quantum mechanical expectation values of various observables, such as physical momentum and quantum energy, but also probability densities, optomechanical dissipation, and fluctuations of the canonical variables.
One of the important research tasks for coupled optomechanical systems is characterizing entanglement between oscillators based on, for example, entropy. Usually, the derivation of entanglement entropy for coupled nonstationary oscillators along this line was carried out up until now with an implicit assumption that all masses (and electric permittivities) of the system are unity or equal to each other [45][46][47][48][49][50]. Notice that such a limited treatment was entirely due to the difficulty in the associated diagonalization procedure of the Hamiltonian. e significance of this research is that we neither adopted such an unnecessary assumption nor used a mathematical approximation when we unfold our theory starting from equation (6). Such clear treatment was possible thanks to the managing of the system in a hybrid way, i.e., by combining the unitary transformation approach and the matrix-diagonalization method together.
An obvious trend in current electronic and optical science is that the related devices become smaller and smaller towards atomic scale as the technology advances. Notice that quantum effects are prominent in devices miniaturized, especially below the scale of the Fermi wavelength [51]. Consequently, quantum treatment of devices including optomechanical ones is important in such cases, while classical mechanics is inadequate in describing their characteristics. e quantum results of this research may provide a theoretical foundation which enables the investigation of the entanglement problem for coupled optomechanical oscillators without a necessity of certain assumptions. 8 Complexity