Nanomechanically Induced Transparency

Abstract

In this paper, we investigate a nanomechanically induced transparency (NIT) effects that arises from the coupling of a nanoelectromechanical system and a trapped ion. By confining the ion in mesoscopic traps and capacitively coupling it with a nanoelectromechanical system suspended as electrodes, the research is intricately focussed on the implications of including the ion’s degrees of freedom. The Lamb-Dicke approximation is crucial to understanding the effects of phonon exchange with electronic qubits and revealing transparency phenomena in this unique coupling. The results underline the importance of the Lamb-Dicke approximation in modelling the effects of transparency windows in nanoelectromechanical systems.

I Introduction

Contemporary scientific research has witnessed several advances, especially in the field of nanotechnology and electromechanical systems. The Colombian interaction between nanoelectromechanical systems and trapped ions stands out even more for its intrinsic rewriting algebra of Quantum Electrodynamics (QED)[1, 2]. This convergence of phenomena offers a unique perspective that could revolutionise technology and open up new frontiers in the understanding of quantum state engineering [3, 4, 5].

Nanoelectromechanical systems (NEM) represent an emerging frontier in the miniaturisation of mechanical devices [6]. With dimensions in the order of nanometres, these systems exhibit unique properties, such as high sensitivity and responsiveness at tiny scales. Trapped ions, on the other hand, are charged particles confined in electromagnetic fields, usually generated by ion traps. These traps offer precise control over the movements and quantum states of the ions, making them ideal components for studying quantum systems and quantum computing applications. Colombian interaction [7, 8], in this context, refers to the electrostatic attraction between a charged particle and a charged plate, i.e. a crucial phenomenon when dealing with trapped ions and a NEM. In the context of the Columbian interaction between NEM and trapped ions, transparency can be achieved by precisely controlling the oscillation of the NEM, resulting in a controlled manipulation of the mechanical response of the materials involved [9].

From an experimental point of view, investigating the presence and role of quantum coherence in a wide variety of systems has attracted considerable interest over the past years [10]. In view of these stimulating results, here we examine the quantum coherence response for a less explored experiment: an ion trap coupled to NEM.

The first proposal of such a system was made by Tian and Zoller [11], with ions in mesoscopic traps, while suspended nanomechanical resonators played the role of tiny electrodes of the trap, acting as nanomechanical resonators with their own degrees of freedom.

Additionally, Hensinger et al. [12] showed that a single trapped ion could be used to test the quantum nature of a mesoscopic mechanical oscillator. The investigation of quantum coherence in mechanical interactions has great relevance for quantum information processing, with proof of experimental feasibility and generation of entanglement [13, 14, 15].

The article is organised as follows. First, we review the NEM-ion trap coupling, modelling and formulating the interaction Hamiltonian for dynamical calculations of a single-phonon system on command. We then analyse the open quantum dynamics of the system that effectively describe the NEM-ion trap coupling. Finally, we obtain the spectra and analyse the induced transparency effects, thereby closing our results and presenting our conclusions.

II The NEM-Ion Trap Coupling

The interaction of a single ion trap and a NEM is based on the experiment proposed in Ref. [12]; see, e.g., an illustration of this Ion-NEM coupling in Fig. 1.

Given this, the corresponding electrostatic energy of such an Ion-NEM system reads

$$V_{c}=k\frac{V_{0}C_{0}}{d+X(t)-x(t)},$$ (1)

where $C_{0}$ denotes the capacitance of the gate, $V_{0}$ represents the voltage bias, $d$ corresponds the separation of the center of mass of the ion and NEM from equilibrium, and $x(X)$ is the position of the ion (NEM). For small oscillations ($X,x<<d$), expanding up to second order, we have [12]

$$V_{c}=-\chi X(t)x(t),$$χかい italic_X ( italic_t ) italic_x ( italic_t ) , (2)

where $\chi=2keV_{0}C_{0}$χかい = 2 italic_k italic_e italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and the first-order term is neglected due to the rapid oscillations.

Figure 1: Schematic model of experiment: vibrational mode of trapped Ion interacting with the electric field of a laser, and electrostatically with a nanoelectromechanical system.

Now, defining $X=\sqrt{\frac{\hbar}{2M\omega}}\left(a+a^{{\dagger}}\right)$ωおめが end_ARG end_ARG ( italic_a + italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ), and $x=\sqrt{\frac{\hbar}{2m\nu}}\left(b+b^{{\dagger}}\right)$νにゅー end_ARG end_ARG ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ), with $M(m)$ being the mass and $\omega(\nu)$ωおめが ( italic_νにゅー ) being the frequency of NEM(Ion) respectively. The Hamiltonian of the whole system can be written as

$$H=\hbar\omega a^{\dagger}a+\hbar\nu b^{\dagger}b+\hbar\frac{\omega_{a}}{2}\sigma_{z}-\hbar\lambda(a+a^{\dagger})(b+b^{\dagger})+H_{e},$$ωおめが italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + roman_ℏ italic_νにゅー italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b + roman_ℏ divide start_ARG italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - roman_ℏ italic_λらむだ ( italic_a + italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) + italic_H start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , (3)

where $\lambda=keV_{0}C_{0}(mM\nu\omega)^{-1/2}/d^{3}$λらむだ = italic_k italic_e italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_m italic_M italic_νにゅー italic_ωおめが ) start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT / italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT is the controllable coupling between the two phononic modes, $\sigma_{z}=|e\rangle\langle e|-|g\rangle\langle g|$σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = | italic_e ⟩ ⟨ italic_e | - | italic_g ⟩ ⟨ italic_g | is the free qubit operator; $a,(a^{{\dagger}})$ are the anihilation (creation) operators of the NEM; and $b,(b^{{\dagger}})$ that of the anihilation (creation) operators of mechanical oscillation of the ion; $H_{e}$ describes the interaction of an external laser with the ion

$$H_{e}=-\vec{\mu}\cdot\vec{E}(x,t),$$μみゅー end_ARG ⋅ over→ start_ARG italic_E end_ARG ( italic_x , italic_t ) ,

where $\vec{\mu}$μみゅー end_ARG is the electric dipole operator for the internal transition. The electric field describing the laser can be written as

$$E(x,t) =A(e^{i(k_{l}x-\omega_{l}t)}+e^{-i(k_{l}x-\omega_{l}t)})$$ ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ) end_POSTSUPERSCRIPT + italic_e start_POSTSUPERSCRIPT - italic_i ( italic_k start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_x - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ) end_POSTSUPERSCRIPT )

$$=A(e^{i[\eta(b+b^{\dagger})-\omega_{l}t]}+e^{-i[\eta(b+b^{\dagger})-\omega_{l}t]})$$ηいーた ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ] end_POSTSUPERSCRIPT + italic_e start_POSTSUPERSCRIPT - italic_i [ italic_ηいーた ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ] end_POSTSUPERSCRIPT )

where $k_{l}$ and $\omega_{l}$ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT are the wave vector and the frequency of the laser, and $\eta=k_{l}\sqrt{\frac{\hbar}{2m\nu}}$ηいーた = italic_k start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT square-root start_ARG divide start_ARG roman_ℏ end_ARG start_ARG 2 italic_m italic_νにゅー end_ARG end_ARG is the Lamb-Dicke parameter. The Hamiltonian $H_{e}$ can be written as:

$$H_e=\mathrm{\hslash }\mathrm{<ruby><rb>\Omega </rb><rt>おめが</rt></ruby>}[{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{+}{e}^{i(\mathrm{<ruby><rb>\eta </rb><rt>いーた</rt></ruby>}(b+b^{†})-{\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_{l}t)}+{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}{e}^{-i(\mathrm{<ruby><rb>\eta </rb><rt>いーた</rt></ruby>}(b+b^{†})-{\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_{l}t)}].$$Ωおめがdelimited-[]subscript𝜎superscript𝑒𝑖𝜂𝑏superscript𝑏†subscript𝜔𝑙𝑡subscript𝜎superscript𝑒𝑖𝜂𝑏superscript𝑏†subscript𝜔𝑙𝑡H_{e}=\hbar\Omega[\sigma_{+}e^{i(\eta(b+b^{\dagger})-\omega_{l}t)}+\sigma_{-}e% ^{-i(\eta(b+b^{\dagger})-\omega_{l}t)}].italic_H start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = roman_ℏ roman_Ωおめが [ italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_ηいーた ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ) end_POSTSUPERSCRIPT + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ηいーた ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t ) end_POSTSUPERSCRIPT ] . (4)

Where $\mathrm{<ruby><rb>\Omega </rb><rt>おめが</rt></ruby>}=epA/\mathrm{\hslash }$Ωおめが𝑒𝑝𝐴Planck-constant-over-2-pi\Omega=epA/\hbarroman_Ωおめが = italic_e italic_p italic_A / roman_ℏ is the effective transition Rabi frequency and $p$ is the dipole matrix element. The qubit dipole transition operators are $\sigma_{-}=|g\rangle\langle e|$σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT = | italic_g ⟩ ⟨ italic_e | and $\sigma_{+}=\sigma_{-}^{{\dagger}}$σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT = italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT.

In the Lamb-Dicke limit, we assuming the ion localized in a region much smaller than the wavelength of the laser $\eta=k_{l}\sqrt{\frac{\hbar}{2m\nu}}\ll 1$ηいーた = italic_k start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT square-root start_ARG divide start_ARG roman_ℏ end_ARG start_ARG 2 italic_m italic_νにゅー end_ARG end_ARG ≪ 1, we can express the dipolar interaction between the ion and the field as

$$H_{e}=\hbar\Omega[\sigma_{+}e^{-i\omega_{l}t}+i\eta\sigma_{+}(b+b^{\dagger})e^{-i\omega_{l}t}+H.c.],$$Ωおめが [ italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT + italic_i italic_ηいーた italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT - italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT + italic_H . italic_c . ] , (5)

The Hamiltonian (3) can be written as:

$$H=H_{o}+H_{I}+H_{e},$$ (6)

where

$$H_{0} =\hbar\omega a^{\dagger}a+\hbar\nu b^{\dagger}b+\hbar\frac{\omega_{a}}{2}\sigma_{z},$$ ωおめが italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + roman_ℏ italic_νにゅー italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b + roman_ℏ divide start_ARG italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ,

$$H_{I} =-\hbar\lambda(a+a^{\dagger})(b+b^{\dagger}),$$ λらむだ ( italic_a + italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) ,

$$H_{e} =\hbar\Omega[\sigma_{+}e^{-i\omega_{l}t}+\sigma_{-}e^{i\omega_{l}t}]$$ Ωおめがdelimited-[]subscript𝜎superscript𝑒𝑖subscript𝜔𝑙𝑡subscript𝜎superscript𝑒𝑖subscript𝜔𝑙𝑡\displaystyle=\hbar\Omega[\sigma_{+}e^{-i\omega_{l}t}+\sigma_{-}e^{i\omega_{l}% t}]= roman_ℏ roman_Ωおめが [ italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT ]

$$+\hbar[ig\sigma_{+}(b+b^{\dagger})e^{-i\omega_{l}t}-ig^{*}\sigma_{-}(b+b^{\dagger})e^{i\omega_{l}t}].$$σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT - italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT - italic_i italic_g start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ( italic_b + italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT ] .

with $g=\mathrm{<ruby><rb>\eta </rb><rt>いーた</rt></ruby>}\mathrm{<ruby><rb>\Omega </rb><rt>おめが</rt></ruby>}$Ωおめがg=\eta\Omegaitalic_g = italic_ηいーた roman_Ωおめが. In the interaction picture, with $U(t)=e^{-iH_{0}t/\hbar}$, we have

$$\mathcal{V}(t) =U^{\dagger}(t)\big{(}H_{I}+H_{e}\big{)}U(t)$$

$$=-\hbar\lambda\big{(}abe^{-i(\omega+\nu)t}+ab^{\dagger}e^{-i(\omega-\nu)t}$$λらむだ ( italic_a italic_b italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT + italic_a italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが - italic_νにゅー ) italic_t end_POSTSUPERSCRIPT

$$+a^{\dagger}be^{i(\omega-\nu)t}+a^{\dagger}b^{\dagger}e^{i(\omega+\nu)t}\big{)}$$ωおめが - italic_νにゅー ) italic_t end_POSTSUPERSCRIPT + italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_ωおめが + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT )

$+\mathrm{\hslash }\mathrm{<ruby><rb>\Omega </rb><rt>おめが</rt></ruby>}[{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{+}{e}^{-i({\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_l-{\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_a)t}+{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}{e}^{i({\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_l-{\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_a)t}]$Ωおめがdelimited-[]subscript𝜎superscript𝑒𝑖subscript𝜔𝑙subscript𝜔𝑎𝑡subscript𝜎superscript𝑒𝑖subscript𝜔𝑙subscript𝜔𝑎𝑡\displaystyle+\hbar\Omega[\sigma_{+}e^{-i(\omega_{l}-\omega_{a})t}+\sigma_{-}e% ^{i(\omega_{l}-\omega_{a})t}]+ roman_ℏ roman_Ωおめが [ italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_t end_POSTSUPERSCRIPT + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_t end_POSTSUPERSCRIPT ]

$$+\hbar\big{[}ig\sigma_{+}be^{i(\omega_{a}-\omega_{l}-\nu)t}+ig\sigma_{+}b^{\dagger}e^{i(\omega_{a}-\omega_{l}+\nu)t}$$σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b italic_e start_POSTSUPERSCRIPT italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT - italic_νにゅー ) italic_t end_POSTSUPERSCRIPT + italic_i italic_g italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT

$$-ig^{*}\sigma_{-}be^{-i(\omega_{a}-\omega_{l}+\nu)t}-ig^{*}\sigma_{-}b^{\dagger}e^{-i(\omega_{a}-\omega_{l}-\nu)t}\big{]}.$$σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT - italic_i italic_g start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT - italic_νにゅー ) italic_t end_POSTSUPERSCRIPT ] .

Assuming that the frequencies of the NEM and the mechanical oscillation of the ion are nearly in resonance ($\omega+\nu\gg\left|\omega-\nu\right|$ωおめが + italic_νにゅー ≫ | italic_ωおめが - italic_νにゅー |), and adjusting the laser frequency for $\omega_{l}=\omega_{a}-\nu$ωおめが start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT = italic_ωおめが start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_νにゅー,

$$\mathcal{V}(t) =-\hbar\lambda\big{(}abe^{-i(\omega+\nu)t}+ab^{\dagger}$$ λらむだ ( italic_a italic_b italic_e start_POSTSUPERSCRIPT - italic_i ( italic_ωおめが + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT + italic_a italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT

$$+a^{\dagger}b+a^{\dagger}b^{\dagger}e^{i(\omega+\nu)t}\big{)}$$ωおめが + italic_νにゅー ) italic_t end_POSTSUPERSCRIPT )

$+\mathrm{\hslash }\mathrm{<ruby><rb>\Omega </rb><rt>おめが</rt></ruby>}[{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{+}{e}^{-i\mathrm{<ruby><rb>\nu </rb><rt>にゅー</rt></ruby>}t}+{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}{e}^{i\mathrm{<ruby><rb>\nu </rb><rt>にゅー</rt></ruby>}t}]$Ωおめがdelimited-[]subscript𝜎superscript𝑒𝑖𝜈𝑡subscript𝜎superscript𝑒𝑖𝜈𝑡\displaystyle+\hbar\Omega[\sigma_{+}e^{-i\nu t}+\sigma_{-}e^{i\nu t}]+ roman_ℏ roman_Ωおめが [ italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_νにゅー italic_t end_POSTSUPERSCRIPT + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_νにゅー italic_t end_POSTSUPERSCRIPT ]

$$+\hbar\big{[}ig\sigma_{+}b+ig\sigma_{+}b^{\dagger}e^{2i\nu t}$$σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b + italic_i italic_g italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_i italic_νにゅー italic_t end_POSTSUPERSCRIPT

$$-ig^{*}\sigma_{-}be^{-2i\nu t}-ig^{*}\sigma_{-}b^{\dagger}\big{]}.$$σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b italic_e start_POSTSUPERSCRIPT - 2 italic_i italic_νにゅー italic_t end_POSTSUPERSCRIPT - italic_i italic_g start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ] .

In the RWA we have

$$\mathcal{H}\approx\hbar\frac{\omega}{2}\sigma_{z}+\hbar\omega a^{{\dagger}}a+\hbar\omega b^{{\dagger}}b-\hbar\lambda\left(ab^{\dagger}+ba^{\dagger}\right)+\hbar g\left(\sigma_{+}b+\sigma_{-}b^{\dagger}\right).$$ωおめが end_ARG start_ARG 2 end_ARG italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + roman_ℏ italic_ωおめが italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + roman_ℏ italic_ωおめが italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b - roman_ℏ italic_λらむだ ( italic_a italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_b italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) + roman_ℏ italic_g ( italic_σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) . (7)

where $\omega=\nu$ωおめが = italic_νにゅー, $ib\to b$ ($-ib^{{\dagger}}\to b^{{\dagger}}$) with Ramsey rotation.

Adding a drive with amplitude $|\epsilon|$ and frequency $\omega_{p}$ωおめが start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT to NEM

$$H_{d}=\hbar(\epsilon a^{{\dagger}}e^{-i\omega_{p}t}+\epsilon^{*}ae^{i\omega_{p}t}),$$ωおめが start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT + italic_ϵ start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT italic_a italic_e start_POSTSUPERSCRIPT italic_i italic_ωおめが start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT ) ,

the full Hamiltonian results

$$\mathcal{H}_{Full} = \hbar\frac{\Delta_{p}}{2}\sigma_{z}+\hbar\Delta_{p}a^{{\dagger}}a+\hbar\Delta_{p}b^{{\dagger}}b-\hbar\lambda\left(ab^{\dagger}+ba^{\dagger}\right)$$ Δでるた𝑝2subscript𝜎𝑧Planck-constant-over-2-pisubscriptΔでるた𝑝superscript𝑎†𝑎Planck-constant-over-2-pisubscriptΔでるた𝑝superscript𝑏†𝑏Planck-constant-over-2-pi𝜆𝑎superscript𝑏†𝑏superscript𝑎†\displaystyle\hbar\frac{\Delta_{p}}{2}\sigma_{z}+\hbar\Delta_{p}a^{{\dagger}}a% +\hbar\Delta_{p}b^{{\dagger}}b-\hbar\lambda\left(ab^{\dagger}+ba^{\dagger}\right)roman_ℏ divide start_ARG roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + roman_ℏ roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_a + roman_ℏ roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_b - roman_ℏ italic_λらむだ ( italic_a italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_b italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) (8)

$$+ \hbar g\left(\sigma_{+}b+\sigma_{-}b^{\dagger}\right)+\hbar(\epsilon a^{{\dagger}}+\epsilon^{*}a),$$ σしぐま start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_b + italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) + roman_ℏ ( italic_ϵ italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + italic_ϵ start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT italic_a ) ,

where ${\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}_p=\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}-{\mathrm{<ruby><rb>\omega </rb><rt>おめが</rt></ruby>}}_d$Δでるた𝑝𝜔subscript𝜔𝑑\Delta_{p}=\omega-\omega_{d}roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = italic_ωおめが - italic_ωおめが start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT.

The presence of dissipation in the circuit dynamics can be described by the master equation [16]

$$\partial_{t}\hat{\rho}={-i\hbar}\left[H_{Full},\hat{\rho}\right]+\frac{\gamma}{2}\mathcal{D}[\sigma_{-}]\hat{\rho}+\frac{\gamma_{\phi}}{2}\mathcal{D}[\sigma_{z}]\hat{\rho}+\sum_{\alpha=a,b}\frac{\kappa_{\alpha}}{2}\mathcal{D}[\alpha]\hat{\rho}$$ρろー end_ARG = - italic_i roman_ℏ [ italic_H start_POSTSUBSCRIPT italic_F italic_u italic_l italic_l end_POSTSUBSCRIPT , over^ start_ARG italic_ρろー end_ARG ] + divide start_ARG italic_γがんま end_ARG start_ARG 2 end_ARG caligraphic_D [ italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ] over^ start_ARG italic_ρろー end_ARG + divide start_ARG italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG caligraphic_D [ italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ] over^ start_ARG italic_ρろー end_ARG + ∑ start_POSTSUBSCRIPT italic_αあるふぁ = italic_a , italic_b end_POSTSUBSCRIPT divide start_ARG italic_κかっぱ start_POSTSUBSCRIPT italic_αあるふぁ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG caligraphic_D [ italic_αあるふぁ ] over^ start_ARG italic_ρろー end_ARG (9)

where $\mathcal{D}[\beta]\hat{\rho}=2\beta\hat{\rho}\beta^{{\dagger}}-\left\{\beta^{{\dagger}}\beta,\hat{\rho}\right\}$βべーた ] over^ start_ARG italic_ρろー end_ARG = 2 italic_βべーた over^ start_ARG italic_ρろー end_ARG italic_βべーた start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - { italic_βべーた start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_βべーた , over^ start_ARG italic_ρろー end_ARG }, for $\beta=\sigma_{-},\sigma_{z},\alpha$βべーた = italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT , italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_αあるふぁ, with $\kappa_{a}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT, $\kappa_{b}$κかっぱ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT and $\gamma$γがんま being the NEM and QuBit relaxation rates, respectively. The $\gamma_{\phi}$γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT is the QuBit decoherence rate. With this we can obtain the following equations of motion

$$⟨\dot{a}⟩=-i{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a⟨a⟩-iϵ+i\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}⟨b⟩,$$Δでるた𝑎delimited-⟨⟩𝑎𝑖italic-ϵ𝑖𝜆delimited-⟨⟩𝑏\left\langle\dot{a}\right\rangle=-i\bar{\Delta}_{a}\left\langle a\right\rangle% -i\epsilon+i\lambda\left\langle b\right\rangle,⟨ over˙ start_ARG italic_a end_ARG ⟩ = - italic_i over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ⟨ italic_a ⟩ - italic_i italic_ϵ + italic_i italic_λらむだ ⟨ italic_b ⟩ , (10a)

$$⟨\dot{b}⟩=-i{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_b⟨b⟩+i\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}⟨a⟩-ig⟨{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}⟩,$$Δでるた𝑏delimited-⟨⟩𝑏𝑖𝜆delimited-⟨⟩𝑎𝑖𝑔delimited-⟨⟩subscript𝜎\left\langle\dot{b}\right\rangle=-i\bar{\Delta}_{b}\left\langle b\right\rangle% +i\lambda\left\langle a\right\rangle-ig\left\langle\sigma_{-}\right\rangle,⟨ over˙ start_ARG italic_b end_ARG ⟩ = - italic_i over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ⟨ italic_b ⟩ + italic_i italic_λらむだ ⟨ italic_a ⟩ - italic_i italic_g ⟨ italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ⟩ , (10b)

$$⟨{\dot{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}}_{-}⟩=-i{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q⟨{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}⟩+ig⟨b{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_z⟩,$$Δでるた𝑞delimited-⟨⟩subscript𝜎𝑖𝑔delimited-⟨⟩𝑏subscript𝜎𝑧\left\langle\dot{\sigma}_{-}\right\rangle=-i\bar{\Delta}_{q}\left\langle\sigma% _{-}\right\rangle+ig\left\langle b\sigma_{z}\right\rangle,⟨ over˙ start_ARG italic_σしぐま end_ARG start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ⟩ = - italic_i over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ⟨ italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ⟩ + italic_i italic_g ⟨ italic_b italic_σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ⟩ , (10c)

where $\kappa_{q}/2=\gamma_{\phi}+\gamma/2$κかっぱ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT / 2 = italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT + italic_γがんま / 2 and ${\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_j={\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}_p-i{\mathrm{<ruby><rb>\kappa </rb><rt>かっぱ</rt></ruby>}}_j/2$Δでるた𝑗subscriptΔでるた𝑝𝑖subscript𝜅𝑗2\bar{\Delta}_{j}=\Delta_{p}-i\kappa_{j}/2over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT - italic_i italic_κかっぱ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT / 2, for $j=a,b,q$. In single phonon limit, we have $\left\langle b\sigma_{z}\right\rangle=-\left\langle b\right\rangle$σしぐま start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ⟩ = - ⟨ italic_b ⟩. With this we have the following solutions in the steady state

$${⟨a⟩}_{ss}=\frac{ϵ\left({\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_b{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q-g^2\right)}{g^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a+{\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}}^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q-{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_b{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q},$$Δでるた𝑏subscript¯Δでるた𝑞superscript𝑔2superscript𝑔2subscript¯Δでるた𝑎superscript𝜆2subscript¯Δでるた𝑞subscript¯Δでるた𝑎subscript¯Δでるた𝑏subscript¯Δでるた𝑞\left\langle a\right\rangle_{ss}=\frac{\epsilon\left(\bar{\Delta}_{b}\bar{% \Delta}_{q}-g^{2}\right)}{g^{2}\bar{\Delta}_{a}+\lambda^{2}\bar{\Delta}_{q}-% \bar{\Delta}_{a}\bar{\Delta}_{b}\bar{\Delta}_{q}},⟨ italic_a ⟩ start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT = divide start_ARG italic_ϵ ( over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT - italic_g start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_g start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT - over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG , (11a)

$${⟨b⟩}_{ss}=\frac{ϵ\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q}{g^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a+{\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}}^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q-{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_b{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q},$$Δでるた𝑞superscript𝑔2subscript¯Δでるた𝑎superscript𝜆2subscript¯Δでるた𝑞subscript¯Δでるた𝑎subscript¯Δでるた𝑏subscript¯Δでるた𝑞\left\langle b\right\rangle_{ss}=\frac{\epsilon\lambda\bar{\Delta}_{q}}{g^{2}% \bar{\Delta}_{a}+\lambda^{2}\bar{\Delta}_{q}-\bar{\Delta}_{a}\bar{\Delta}_{b}% \bar{\Delta}_{q}},⟨ italic_b ⟩ start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT = divide start_ARG italic_ϵ italic_λらむだ over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG start_ARG italic_g start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT - over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG , (11b)

$${⟨{\mathrm{<ruby><rb>\sigma </rb><rt>しぐま</rt></ruby>}}_{-}⟩}_{ss}=\frac{ϵ\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}g}{{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_b{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q-g^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_a-{\mathrm{<ruby><rb>\lambda </rb><rt>らむだ</rt></ruby>}}^2{\overline{\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}}_q},$$Δでるた𝑎subscript¯Δでるた𝑏subscript¯Δでるた𝑞superscript𝑔2subscript¯Δでるた𝑎superscript𝜆2subscript¯Δでるた𝑞\left\langle\sigma_{-}\right\rangle_{ss}=\frac{\epsilon\lambda g}{\bar{\Delta}% _{a}\bar{\Delta}_{b}\bar{\Delta}_{q}-g^{2}\bar{\Delta}_{a}-\lambda^{2}\bar{% \Delta}_{q}},⟨ italic_σしぐま start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ⟩ start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT = divide start_ARG italic_ϵ italic_λらむだ italic_g end_ARG start_ARG over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT - italic_g start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_λらむだ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG roman_Δでるた end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG , (11c)

For we investigate the phenomenon of induced transparency in this paper, we’ll call it nanomechanically induced transparency (NIT). It is important to highlight that the occurrence of the NIT phenomenon requires that $\kappa_{a}\gg\kappa_{b},\gamma,\gamma_{\phi}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ≫ italic_κかっぱ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , italic_γがんま , italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT. To plot the graphs we can use experimental data recorded in the literature, given with $\lambda\approx 300$kHzきろへるつitalic_λらむだ ≈ 300 italic_k italic_H italic_z [12], and $g\approx 500$kHzきろへるつitalic_g ≈ 500 italic_k italic_H italic_z [17] for $Cd^{+}$, and $\kappa_{q}=11.9$kHzきろへるつitalic_κかっぱ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = 11.9 italic_k italic_H italic_z [18]. The other parameters can be modified by the manufacturing technology of the experimental elements. Here we will only show their implications. Next, we will investigate the effect of transmission and absorption due to the oscillatory mechanical parameters analogous to the electromagnetic parameters of Ref. [9, 19].

Figure 2: Absorption $Im\left\langle a\right\rangle_{ss}$ (red dashed line) and dispersion $Re\left\langle a\right\rangle_{ss}$ (blue solid line) of the cavity mode $a$ when coupled to trapped ion and cavity mode $b$ as a function of the normalized detuning ${\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}_p/{\mathrm{<ruby><rb>\kappa </rb><rt>かっぱ</rt></ruby>}}_a$Δでるた𝑝subscript𝜅𝑎\Delta_{p}/\kappa_{a}roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT / italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT. In this case we use the parameters: a) $\left|\epsilon\right|=0.03\kappa_{a},\lambda=g=0.5\kappa_{a},\kappa_{b}=10^{-3}\kappa_{a},\gamma=\gamma_{\phi}=10^{-3}\kappa_{a}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_λらむだ = italic_g = 0.5 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_κかっぱ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま = italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT; b) $\left|\epsilon\right|=0.03\kappa_{a},\lambda=1.0\kappa_{a},g=0.15\kappa_{a},\gamma=\gamma_{\phi}=10^{-3}\kappa_{a}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_λらむだ = 1.0 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_g = 0.15 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま = italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT.

Figure 3: Absorption $Im\left\langle a\right\rangle_{ss}$ (red dashed line) and dispersion $Re\left\langle a\right\rangle_{ss}$ (blue solid line) of the cavity mode $a$ when coupled to trapped ion and cavity mode $b$ as a function of the normalized detuning ${\mathrm{<ruby><rb>\Delta </rb><rt>でるた</rt></ruby>}}_p/{\mathrm{<ruby><rb>\kappa </rb><rt>かっぱ</rt></ruby>}}_a$Δでるた𝑝subscript𝜅𝑎\Delta_{p}/\kappa_{a}roman_Δでるた start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT / italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT. In this case we use the parameters: a) $\left|\epsilon\right|=0.03\kappa_{a},\lambda=g=0.5\kappa_{a},\kappa_{b}=10^{-3}\kappa_{a},\gamma=10^{-2}\kappa_{a},\gamma_{\phi}=10^{-1}\kappa_{a}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_λらむだ = italic_g = 0.5 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_κかっぱ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま = 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT; b) $\left|\epsilon\right|=0.03\kappa_{a},\lambda=g=0.5\kappa_{a},\gamma=10^{-1}\kappa_{a},\gamma_{\phi}=1.0\kappa_{a}$κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_λらむだ = italic_g = 0.5 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま = 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_γがんま start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT = 1.0 italic_κかっぱ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT.

III Results and Discussion

In Fig. (2), we show the behavior of the real and imaginary parts of the operator $\left\langle a\right\rangle_{ss}$ given by the Eq.(IIa). As can be seen in Fig (2.a), when the couplings are equal, the transparency windows have a symmetrical profile. However, when there is a discrepancy between the coupling values, the transparency windows can be significantly affected, as can be seen in Fig. (2.b). This fact is discussed in Ref. [19] and is related to the interference between the different absorption paths involving the closest energy levels that are in the linewidth of the system. We emphasize that the width and height of the central peak in Fig. (2) are strongly dependent on the transition rates between the ground state and the first excited states. Specifically, these rates can be equal at certain values of g. In such cases, selecting values of g where these rates are equal can result in different absorption profiles for the system. For instance, by choosing $\left(g\sim 0.15\right)$, we obtain the absorption profile depicted in Fig. (2.b). Indeed, increasing the value of g leads to an increase in both the height and width of the central peak, as observed in Fig. (2.a). In this way, we have a system that can be controlled by adjusting the system coupling parameters. Furthermore, in this system it is possible to have Fano-type resonance by adjusting the coupling parameters. This type of resonance arise of the interference of eigenstates whose energy levels are partially inside and outside the linewidth of the system. For example, for the case ($g=3.0\lambda$λらむだ), the absorption exhibits an asymmetric shape. In comparison to the study described in Ref.[20], we show that it is possible to obtain a broader scenario of the system’s absorption profile turning on the coupling parameters of the system. We also observed that dissipative ion channels can substantially modify the absorption profile of the system as shown in Fig. (3.a). As can be seen in Fig. (3.b), the dephasing dissipative channel is what most damages the system’s absorption. In fact, it destroys the central absorption peak, providing a strong indication of a destructive interference effect on the system.

IV Summary

We show how the effects of ion degrees of freedom in the Lamb-Dicke approximation can be relevant in the emergence of quantum interference phenomena in the context of nanoelectromechanical systems. The nanomechanically induced transparency (NIT) phenomenon arises due to interference of the absorption paths that are in the linewidth of the system. We also show that the profile of the absorption spectrum can be controlled by tuning the coupling parameters between the ion and the nanoelectromechanical resonators modes. The results presented can be pave the way for the intense investigations of the analog electromagnetically induced transparency in nanoscale.

Acknowledgments. E. C. Diniz acknowledge the UNEMAT (Campus Tangará da Serra) for the hospitality.