Controlling the Pump-Probe Optical Response in Asymmetric Tunneling-Controlled Double Quantum Dot Molecule—Metal Nanoparticle Hybrids

: In the present work, we investigate the modiﬁed nonlinear pump-probe optical properties due to the excitonic–plasmonic interaction of a double semiconductor quantum dot (SQD) molecule coupled to a metal nanoparticle (MNP). More speciﬁcally, we study the absorption and the dispersion spectra of a weak electromagnetic ﬁeld in a hybrid structure with two counterparts, a molecule of two coupled SQDs, and a spherical MNP driven by a ﬁeld of high intensity. We solve the relevant density matrix equations, calculate the ﬁrst-order optical susceptibility of the probe ﬁeld in the strong pumping regime, and investigate the way in which the distance between the two counterparts modiﬁes the optical response, for a variety of values of the physical constants of the system, including the pump-ﬁeld detuning, the tunnelling rate, and the energy separation gap associated with the excited states of the coupled SQDs.

In some of these studies, the nonlinear optical effects in an SQD-MNP hybrid structure are examined under interaction with a pump field of high intensity and a tunable probe field of much lower intensity [7,18,23,24,26,28,34,36]. A specific SQD nanostructure that exhibits optical properties of high importance is the asymmetric molecule consisting of a couple of SQDs [39][40][41][42][43][44][45][46][47][48][49][50][51][52]. The band structures of the SQDs are considered to be different from each other, while the rate at which the tunneling effects take place is determined by the value of the applied voltage. More recently, the optical properties of a hybrid structure composed of an asymmetric double SQD molecule coupled to a MNP have been analyzed [32,38]. In these studies, the construction of a diffraction grating, the special characteristics of which can be modified due to plasmonic mediated interaction [32], and the control and the energy absorption rate spectrum [38] were explored. Here, we analyze

Methods
The hybrid system studied has two counterparts: a spherical MNP with dielectric function ε m (ω) and a molecule of two coupled SQDs with dielectric constant ε S , as seen in Figure 1. Moreover, we consider that the condition R > a is met, with R and a respectively denoting the interparticle distance and the radius of the MNP, and we assume that the size of the SQD is negligible as compared to the size of the MNP. The structure is surrounded by a material with dielectric constant ε env .
Appl. Sci. 2021, 11, x FOR PEER REVIEW 2 of 15 diffraction grating, the special characteristics of which can be modified due to plasmonic mediated interaction [32], and the control and the energy absorption rate spectrum [38] were explored. Here, we analyze the response of a system of an MNP coupled to a doublet of SQDs interacting with a pump field of high intensity. Specifically, we combine the theoretical approaches presented in Refs. [36,46] and calculate the optical susceptibility of the probe field, while a second field of much higher intensity is used for the pumping of the system. Then, we analyze the spectral profile of the real and the imaginary part of the susceptibility, for different values of the physical parameters of the structure. We find that, at the exact-resonance case with the pump field, the spectra present exact symmetry and are composed of a quartet of resonances. In the off-resonance case, the spectra exhibit six resonances and are not symmetric, for average values of the tunneling coefficient. For high interparticle distances, analytical expressions for the position of the resonances are provided based on the dressed-states analysis, which is valid in the strong intensity regime applied in the present paper. For low tunneling rates, the contribution of MNP outweighs the contribution of the coupled SQDs and it also introduces a non-negligible slope to the spectra of the total susceptibility. However, for higher values of tunneling coupling coefficient, the contribution of the double SQD molecule seems to play a more important role in the spectral profile of the total susceptibility, determining to a substantial extent the shape of its imaginary part. In every case, as the distance between the counterparts decreases, the resonances are transposed away from the center of the spectra.
This study is structured as follows. In the first part, the details of the applied methodologies are presented and the equations of motion are derived in the density matrix formalism, which should be numerically solved, in order to calculate the optical susceptibility of the probe field as it interacts with the strongly pumped hybrid nanostructure, and, in the second part, a dressed state representation of the system is introduced in the analysis. Then, in the third part of this work, the results associated with the behavior of the optical susceptibility are analyzed, for different parameters of the structure. Finally, in the last part of this study, the conclusions of the above analysis are summarized.

Methods
The hybrid system studied has two counterparts: a spherical MNP with dielectric function ( ) m ε ω and a molecule of two coupled SQDs with dielectric constant S ε , as seen in Figure 1. Moreover, we consider that the condition > R a is met, with R and a respectively denoting the interparticle distance and the radius of the MNP, and we assume that the size of the SQD is negligible as compared to the size of the MNP. The structure is surrounded by a material with dielectric constant env ε .  Schematic diagram of the MNP-double-SQD hybrid system, a, b, and R respectively representing the radius of MNP, the radius of the SQDs, and the center-to-center distance between an SQD and the MNP. Between the SQDs, electron tunneling is present.
The hybrid system interacts with two linearly polarized, oscillating electromagnetic fields, where E i and ω i respectively represent the corresponding amplitude and the angular Appl. Sci. 2021, 11, 11714 3 of 14 frequency (i = a, for the pump field and i = b, for the probe field), that excite plasmons on the surface of the MNP, providing a strong, continuous spectral response. An energylevel scheme is presented in Figure 2, where the possible transitions are also depicted. The coupling between the two SQDs with different band structures, which are typically placed in distances between 7 and 10 nm, owes its presence to the electron tunneling. The coupled SQDs are represented as a three-level quantum system [32,[38][39][40][41][42][43][44][45][46][47][48][49][50][51][52]. With |0 we denote the ground state, and |1 (direct exciton) corresponds to the state describing a single electron excitation, in which case the excited electron, as well as the hole that is created, are both confined within the same quantum dot. On the other hand, an electron can be excited in state |2 (indirect exciton) when it is conveyed from the first quantum dot to the second one, due to tunneling effects following the excitation pathway. The incident field induces an oscillating dipole moment and excites electrons from the ground state |0 to the (direct) exciton state, where an electron-hole pair is bound to the same SQD. The direct exciton is also coupled to the states of the surface plasmons due to the long-range Coulomb interaction that permits the exchange of energy between the counterparts [53][54][55][56]. The tunneling barrier can be controlled by properly adjusting the voltage of an external gate electrode structure [39][40][41][42][43][44][45][46][47][48][49][50][51][52].
tunneling is present.
The hybrid system interacts with two linearly polarized, oscillating electromagnetic fields, where i E and ω i respectively represent the corresponding amplitude and the angular frequency ( = i a, for the pump field and = i b , for the probe field), that excite plasmons on the surface of the MNP, providing a strong, continuous spectral response. An energy-level scheme is presented in Figure 2, where the possible transitions are also depicted. The coupling between the two SQDs with different band structures, which are typically placed in distances between 7 and 10 nm, owes its presence to the electron tunneling. The coupled SQDs are represented as a three-level quantum system [32,38,[39][40][41][42][43][44][45][46][47][48][49][50][51][52]. With 0 we denote the ground state, and 1 (direct exciton) corresponds to the state describing a single electron excitation, in which case the excited electron, as well as the hole that is created, are both confined within the same quantum dot. On the other hand, an electron can be excited in state 2 (indirect exciton) when it is conveyed from the first quantum dot to the second one, due to tunneling effects following the excitation pathway. The incident field induces an oscillating dipole moment and excites electrons from the ground state 0 to the (direct) exciton state, where an electron-hole pair is bound to the same SQD. The direct exciton is also coupled to the states of the surface plasmons due to the long-range Coulomb interaction that permits the exchange of energy between the counterparts [53][54][55][56]. The tunneling barrier can be controlled by properly adjusting the voltage of an external gate electrode structure [39][40][41][42][43][44][45][46][47][48][49][50][51][52]. Schematic energy-level configuration of the hybrid system. The MNP is described as a continuum of energy states, while the asymmetric double SQD molecule is depicted as a three-level Λ-type system with a ground state 0 , a direct exciton state 1 (electron-hole pair in the same quantum dot), and an indirect exciton state 2 (due to the transferring of the excited electron to the excited state of the second quantum dot). The application of a threshold voltage to the double SQD induces strong tunneling between the upper levels, while, at the same time, the ground energy levels are effectively decoupled.
In Equation (1), we present the Hamiltonian that describes the interaction of the hybrid system with a pump and a probe field, as depicted in Figure 2, in the dipole approximation.
The first three terms in Equation (1) describe the unperturbed system, with n Ε being the energy of level n (for the ground state, we take 0 0 Ε = ), while the fourth Figure 2. Schematic energy-level configuration of the hybrid system. The MNP is described as a continuum of energy states, while the asymmetric double SQD molecule is depicted as a three-level Λ-type system with a ground state |0 , a direct exciton state |1 (electron-hole pair in the same quantum dot), and an indirect exciton state |2 (due to the transferring of the excited electron to the excited state of the second quantum dot). The application of a threshold voltage to the double SQD induces strong tunneling between the upper levels, while, at the same time, the ground energy levels are effectively decoupled.
In Equation (1), we present the Hamiltonian that describes the interaction of the hybrid system with a pump and a probe field, as depicted in Figure 2, in the dipole approximation.
The first three terms in Equation (1) describe the unperturbed system, with E n being the energy of level |n (for the ground state, we take E 0 = 0), while the fourth term owes its presence to the tunneling effects between the SQDs, with T e expressing the rate at which electrons tunnel within the SQD molecule. The last term in Equation (1) describes the interaction between the electromagnetic fields and the SQD molecule. Here, µ and E SQD respectively represent the dipole moment associated to the excitation |0 ↔ |1 and the total field detected at the position of the SQD molecule, which is given by the formula where and ε e f f S = (ε S + 2ε env )/3ε env . In the above, we have used the quasistatic approximation, where retardation effects are omitted, following several previous works [53][54][55][56]. The first term in the expression of Equation (3) is ascribed to the direct coupling between the applied electromagnetic field and the SQDs. We expect that the MNP interacts with the applied field, as well as the field emitted by the SQDs, and as a result, it emits in its surrounding space. The coupling between the SQDs and the first (second) component of this emitted field is associated with the second term of Equation (3) (the G i parameter). The parameters s a and s n , in the case of an incident field polarized along the direction of the interparticle axis, are respectively equal to 2 and (n + 1) 2 , while, for an incident field that is polarized perpendicularly to the direction of the interparticle axis, we take s a = −1 and s n = n(n + 1)/2. The G i parameter is calculated including the multipole effects for the MNP [55].
Based on the expression of the Hamiltonian that describes the system, in the rotating wave approximation, a set of nine density-matrix differential equations is derived. The number of these equations is eliminated to five, under the conditions ∑ 2 0 ρ ii = 1 and ρ mn = ρ * nm , with n = m: . . . . .
where the parameters δ = ω b − ω a and ∆ = (E 1 − E 0 )/ − ω a respectively denote the difference of the angular frequencies with reference to the pump and the probe field and the detuning from the energy gap between the energy levels |0 and |1 . The |1 ↔ |2 transition corresponds to the energy difference ω 12 . The dephasing rates are symbolized with γ nm , while the population decay rate for the transition |n → |m is indicated with Γ nm .
To calculate the first-order susceptibility of the hybrid structure, aiming at the examination of the pump-probe optical response, we apply to the first-order Taylor expansion of ρ nm , according to the first-order approximation, in the probe field Rabi frequency: nm .
In the expansions presented in Equation (10), the terms ρ (2) nm and ρ (3) nm are considered small, in the sense that ρ nm . Since the intensity of the pump field is assumed to be four orders of magnitude higher than the intensity of the probe field (Ω a Ω b ), we should maintain in our calculations all the orders of its interaction with the hybrid structure. By substituting the expanded density-matrix elements in Equations (5)-(9), the set of Equations (11)-(25) is derived. . . . . . . . . . . . . . .
where we have introduced the simplified notation for the coefficients defined in Equations (10) and G = G a + G b . For a hybrid system that is strongly pumped, the probe-field first-order absorption and the dispersion are respectively equal to −Im[ρ (2) 01 ] and Re[ρ (2) 01 ]. More specifically, the first-order probe-field optical susceptibility, associated with the double SQD molecule, under the interaction with a pump field, is described by the following analytical expression: where the coefficient ρ (2) 01 (T) is calculated at time T, which is a sufficient time interval for our system to reach its steady state. In Equation (27), ε 0 denotes the dielectric constant on Appl. Sci. 2021, 11, 11714 6 of 14 vacuum, and with V and Γ we respectively indicate the SQD volume and the optical confinement factor [7]. To perform this calculation, we numerically solve Equations (11)- (26).
The expression describing the MNP optical susceptibility at the first order [6] is given by the following expression: SQD .
At this point, if we sum the first-order susceptibilities of the individual components, we derive an expression for the optical susceptibility for the total hybrid structure at the first order: The dressed-states picture provides us with an illuminating physical insight for the results that are presented in the next section, in the specific case of a double SQD molecule non-interacting with the MNP (practically, for an interparticle distance above 50 nm). Under this condition, the expression of the Hamiltonian, given by Equation (1), is simplified. Under the assumption of interaction only with a strong pump field, the population amplitudes of the energy levels satisfy the expression If we apply the transformations a 0 (t) = c 0 (t), a 1 (t) = c 1 (t)e i(∆−iE 1 / )t , and a 2 (t) = c 2 (t)e i(∆−ω 12 −iE 2 / )t in Equation (28), the Hamiltonian is written as follows: The solution to the eigenvalue/eigenstate problem |H b − λI| = 0, under the condition ∆ = ω 12 , yields the eigenvalue λ 0 = 0, which is associated with the dark state |c = (T e |0 + Ω a |2 )/ Ω 2 a + T 2 e , as well as the eigenvalues λ ± = ∆ ± ∆ 2 + 4Ω 2 a + 4T 2 e /2, corresponding to the dressed states Thus, in the dressed-state description, the strong pumping splits the two SQD energy levels, creating two triplets, as presented in Figure 3. The energy gap between the initial levels and the new ones is equal to ∆/2 ± ∆ 2 + 4Ω 2 a + 4T 2 e /2. A schematic representation of the transitions taking place between the dressed-state energy levels is presented in Figure 3a, for a positive value of the detuning of the pump field (∆ = ω 12 > 0). These transitions are associated with the following values of the pump-probe field energy mismatch: λ + + |λ − |, λ + , |λ − |, −|λ − |, −λ + , −λ + − |λ − |. In the case of a pump field that is found at exact resonance with the |0 ↔ |1 transition (∆ = ω 12 = 0), the dressed states are equally separated from each other, which is the case depicted in Figure 3b. Here, the energy of the |c ↔ |+ (|z ↔ |− ) transition is equal to the energy of the |− ↔ |c (|+ ↔ |z ) tran-sition, and hence, a doublet of resonances, arisen around a central frequency, is responsible for the two sidebands of the FWM spectra.

Results
Here, we investigate the characteristics associated with the absorption and the dispersion spectra of a probe field applied to an SQD-MNP hybrid structure, under the interaction with a strong pump field, for various interparticle distances. In the following paragraphs, we consider a linearly polarized strong pump field. Its polarization is taken parallel to the interparticle axis ( 2) a s = , while its Rabi frequency a Ω  is set equal to 3 meV . The weak probe field Rabi frequency is considered to be 100 times lower than

Results
Here, we investigate the characteristics associated with the absorption and the dispersion spectra of a probe field applied to an SQD-MNP hybrid structure, under the interaction with a strong pump field, for various interparticle distances. In the following paragraphs, we consider a linearly polarized strong pump field. Its polarization is taken parallel to the interparticle axis (s a = 2), while its Rabi frequency Ω a is set equal to 3 meV. The weak probe field Rabi frequency is considered to be 100 times lower than the pump field Rabi frequency. The SQD dielectric constant is ε S = 6ε 0 , while for the modeling of the MNP we take the experimental values of the dielectric function referring to gold [57]. The radius a of the MNP is 7.5 nm. The hybrid structure is considered to be surrounded by vacuum and hence, ε env = ε 0 . For the SQD, we consider values typical for InAs/GaAs SQDs that have been used in previous studies [38,42,46]. As far as the dephasing and decay rates related to the double-SQD molecule are concerned, we respectively have γ 01 = 0.22 meV, γ 12 = 0.1 meV, γ 02 = 0.001 meV, and Γ 10 = 0.025 meV, Γ 20 = 0.000025 meV, and Γ 12 = 0. Finally, the energy gap associated with the |1 ↔ |2 transition is ω 12 = 1.4 eV and the dipole moment element µ of the |0 ↔ |1 transition is equal to 0.65 · e nm.
More specifically, in Figures 4-10, we present the spectra of the asymmetric double SQD molecule susceptibility, χ SQD , the MNP susceptibility, χ MNP , and the total susceptibility of the hybrid structure, χ (1) tot , as a function of the pump-probe field detuning δ, for several values of the distance R that separates the SQD from the MNP (R = 100 nm: blue solid curve, 20 nm: red dashed-dotted curve, and 15 nm: green dashed-dotted curve). In Figure 4, the dispersion and the absorption spectra of the probe field are presented for ∆ = ω 12 = 0, in the case of a quite low value of the tunneling rate T e = 0.01 meV, as seen in Figure 4a,c,e and b,d,f, respectively. First, we observe that, at exact resonance, the absorption/dispersion associated with the double SQD molecule exhibits a symmetric/antisymmetric profile for every value of the tunneling rate coefficient. In Figure 4, since the tunneling rate tends to zero, the spectra of the SQD optical susceptibility, as seen in Figure 4a,b, are identical to the ones plotted for a single SQD-MNP hybrid [36]. As the value of R decreases, the sideband resonances are transposed away from the spectral center. For the set of parameters used in Figures 4-6 (∆ = ω 12 = 0), the MNP contribution outweighs the contribution of the double SQD molecule, introducing a slope to the spectra of the total susceptibility due to the presence of the first term in the expression of the MNP susceptibility, which is in accordance with the mathematical expression of Equation (28). optical susceptibility, as seen in Figure 4a,b, are identical to the ones plotted for a single SQD-MNP hybrid [36]. As the value of R decreases, the sideband resonances are transposed away from the spectral center. For the set of parameters used in Figures 4-6 12 0) ( ω Δ = = , the MNP contribution outweighs the contribution of the double SQD molecule, introducing a slope to the spectra of the total susceptibility due to the presence of the first term in the expression of the MNP susceptibility, which is in accordance with the mathematical expression of Equation (28).   Figure 4). However, the amplitude of the resonances is slightly decreased as the distance between the particles is increased. Finally, and most importantly, as we note by comparing Figure 4 to Figure 5, when the tunneling rate coefficient increases, the double SQD molecule susceptibility seems to have an even more important contribution to the spectral profile of the total susceptibility (1) χ tot . In Figure 6, we investigate the optical response of the two counterparts, as well as the total optical response of the system, for a high electron-tunneling rate ( 5 meV) e T =  . Here, as the interparticle distance decreases, we note an infinitesimal transposition of the resonances, while their intensity is slightly increased, exhibiting an inverse trend with respect to the one observed in Figure 5  The susceptibility spectra presented in Figure 5 refer to the case of the non-negligible electron-tunneling rate T e = 0.5 meV. Here, we observe a set of inner resonances arising around the positions ± Ω 2 a + T 2 e , as well as a set of outer resonances, placed at ±2 Ω 2 a + T 2 e , which are hardly noticeable due to their low amplitude. Moreover, the transposition of the dominant inner resonances, owing its presence to the alteration of the value of R, is of low importance as compared to the transposition of the resonances observed in the case with T e = 0.01 meV (Figure 4). However, the amplitude of the resonances is slightly decreased as the distance between the particles is increased. Finally, and most importantly, as we note by comparing Figure 4 to Figure 5, when the tunneling rate coefficient increases, the double SQD molecule susceptibility seems to have an even more important contribution to the spectral profile of the total susceptibility χ (1) tot . In Figure 6, we investigate the optical response of the two counterparts, as well as the total optical response of the system, for a high electron-tunneling rate ( T e = 5 meV). Here, as the interparticle distance decreases, we note an infinitesimal transposition of the resonances, while their intensity is slightly increased, exhibiting an inverse trend with respect to the one observed in Figure 5, for T e = 0.5 meV. In addition, the distance between the inner resonances is doubled. Furthermore, the double SQD molecule contribution almost entirely defines the spectral profile of Im[χ (1) tot ], as shown in Figure 6f, while the MNP mainly adds a shift to the spectral profile due to the contribution of the double SQD molecule, both co-shaping the profile of the Re[χ (1) tot ] spectrum (Figure 6e). The rest of the trends that are followed are similar to the ones presented in Figure 4, the most important of them being the significant transposition of the resonances when the center-to-center distance is modified. Figure 6. The dispersion (a,c,e) and the absorption (b,d,f) spectra for the same parameters as in Figure 4, but with T e = 5 meV. Also, R = 100 nm (blue solid curve), 20 nm (red dashed curve), and 15 nm (green dashed-dotted curve).
In the next three figures, we study the off-resonant case, with ∆ = −1 meV, for various values of the tunneling rate coefficient (Figure 7: T e = 0.01 meV, Figure 8: T e = 0.5 meV, and Figure 9: T e = 5 meV), where we take ω 12 = 1 meV. After carefully studying all the possible cases, it becomes evident that, for ∆ = 0, the spectral profiles lose their symmetry. By comparing Figure 7a,b with the corresponding Figure 4 (for the same electron-tunneling rate, yet with ∆ = ω 12 = 0), we observe that, for T e = 0.01 meV, the general profile of the dispersion spectrum of the one figure roughly resembles the profile of the absorption spectrum of the other figure. The rest of the trends that are followed are similar to the ones presented in Figure 4, the most important of them being the significant transposition of the resonances when the center-to-center distance is modified.
In the case of an average electron-tunneling rate ( T e = 0.5 meV), for ∆ = − ω 12 = −1 meV, which is presented in Figure 8, we note that the typical six-resonance pattern arises on the spectra of the SQDs' susceptibility, as seen in Figure 8a,b. These patterns usually appear in three-level-system spectra, unless we are in a dark state with exactly resonant transitions or at two-photon resonance. Furthermore, the profile of the total susceptibility is ascribed both to the contributions of the coupled SQDs, as well to the MNP, as seen in Figure 8c The rest of the trends that are followed are similar to the ones presented in Figure 4, the most important of them being the significant transposition of the resonances when the center-to-center distance is modified. In the case of an average electron-tunneling rate ( 0.5 meV) e T =  , for 12 1 meV , which is presented in Figure 8, we note that the typical six-resonance pattern arises on the spectra of the SQDs' susceptibility, as seen in Figure  8a,b. These patterns usually appear in three-level-system spectra, unless we are in a dark state with exactly resonant transitions or at two-photon resonance. Furthermore, the profile of the total susceptibility is ascribed both to the contributions of the coupled SQDs, as well to the MNP, as seen in Figure 8c,d. It becomes apparent that the contribution of MNP is responsible for the characteristic shift observed in (1) tot χ . 2 a e T ± Ω +  , the outer being negligible regarding its magnitude. Except for the loss of symmetry, which is detected in the spectral response, for this specific set of parameters, the spectra present additional similarities to the ones presented in Figure 6 for the same values of the energy gap 12 ω  and tunneling rate e T , yet with a zero pump field detuning Δ . It is important to note that, even for a comparatively high tunneling rate, the position of the inner resonance arising for 0 δ > moves towards higher values of the detuning parameter as we decrease the interparticle distance, in which case the two counterparts of the hybrid structure interact more strongly with each other. This interesting feature was not observed in the spectra of the susceptibility at exact-resonance for the same value of the electron-tunneling rate ( Figure 6). However, if we further increase the tunneling-rate coefficient ( T e = 5 meV), as seen in Figure 9, two inner resonances and two outer resonances arise again, respectively emerging at the positions ± Ω 2 a + T 2 e and ±2 Ω 2 a + T 2 e , the outer being negligible regarding its magnitude. Except for the loss of symmetry, which is detected in the spectral response, for this specific set of parameters, the spectra present additional similarities to the ones presented in Figure 6 for the same values of the energy gap ω 12 and tunneling rate T e , yet with a zero pump field detuning ∆. It is important to note that, even for a comparatively high tunneling rate, the position of the inner resonance arising for δ > 0 moves towards higher values of the detuning parameter as we decrease the interparticle distance, in which case the two counterparts of the hybrid structure interact more strongly with each other. This interesting feature was not observed in the spectra of the susceptibility at exact-resonance for the same value of the electron-tunneling rate ( Figure 6). Finally, in Figure 10, where we consider exactly the same values for the other parameters as in Figure 7, except for 1 meV Δ =  , we observe that the change of the detuning sign is responsible for the inversion of symmetry with respect to the origin (0,0) , as far as the spectra of the double SQD molecule susceptibility are concerned (Figure 10a,b). Adding an extra shift, due to the MNP contribution (Figure 10c,d), we finally take the spectra of the total susceptibility (Figure 10e,f).

Conclusions
In summary, we perform calculations in order to investigate the optical response of a strongly pumped hybrid structure made of an SQD molecule and a spherical MNP. Finally, in Figure 10, where we consider exactly the same values for the other parameters as in Figure 7, except for ∆ = 1 meV, we observe that the change of the detuning sign is responsible for the inversion of symmetry with respect to the origin (0, 0), as far as the spectra of the double SQD molecule susceptibility are concerned (Figure 10a,b). Adding an extra shift, due to the MNP contribution (Figure 10c,d), we finally take the spectra of the total susceptibility (Figure 10e,f). Finally, in Figure 10, where we consider exactly the same values for the other parameters as in Figure 7, except for 1 meV Δ =  , we observe that the change of the detuning sign is responsible for the inversion of symmetry with respect to the origin (0,0) , as far as the spectra of the double SQD molecule susceptibility are concerned (Figure 10a,b). Adding an extra shift, due to the MNP contribution (Figure 10c,d), we finally take the spectra of the total susceptibility (Figure 10e,f).

Conclusions
In summary, we perform calculations in order to investigate the optical response of a strongly pumped hybrid structure made of an SQD molecule and a spherical MNP.

Conclusions
In summary, we perform calculations in order to investigate the optical response of a strongly pumped hybrid structure made of an SQD molecule and a spherical MNP. Starting with the Hamiltonian, the density matrix equations are derived and the firstorder optical susceptibility of the weak probe field is calculated, for the total system, as well as for its two counterparts, respectively. We find that, in the exact-resonance pump field regime, the spectra of the real and the imaginary part of the optical susceptibility are four-peaked and exhibit a symmetric profile, while, in the off-resonance regime, the spectra exhibit six resonances and are not symmetric for average tunneling rates. For high interparticle distances, analytical expressions for the position of the resonances are provided, based on the dressed states analysis, which is valid within a specific range of values of the Rabi frequency. For low tunneling rates, the contribution of MNP outweighs the contribution associated with the asymmetric double SQD molecule, introducing a non-negligible slope to the spectra of the total susceptibility. However, for a stronger tunneling rate, the contribution of the SQD molecule seems to play an even more important role to the spectral profile of the total susceptibility, determining to a substantial extent the shape of the spectrum of its imaginary part. In every case studied, as the distance between the counterparts is decreased, the resonances are transposed away from the center of the spectra.