Bloch Oscillations in the Chains of Artiﬁcial Atoms Dressed with Photons

: We present a model of one-dimensional chain of two-level artiﬁcial atoms driven with DC ﬁeld and quantum light simultaneously in a strong coupling regime. The interaction of atoms with light leads to electron-photon entanglement (dressing of the atoms with light). The driving via dc ﬁeld leads to the Bloch oscillations (BO) in the chain of dressed atoms. We consider the mutual inﬂuence of dressing and BO and show that scenario of oscillations dramatically differs from predicted by the Jaynes-Cummings and Bloch-Zener models. We study the evolution of the population inversion, tunneling current, photon probability distribution, mean number of photons, and photon number variance, and show the inﬂuence of BO on the quantum-statistical characteristics of light. For example, the collapse-revivals picture and vacuum Rabi-oscillations are strongly modulated with Bloch frequency. As a result, quantum properties of light and degree of electron-photon entanglement become controllable via adiabatic dc ﬁeld turning. On the other hand, the low-frequency tunneling current depends on the quantum light statistics (in particular, for coherent initial state it is modulated accordingly the collapse-revivals picture). The developed model is universal with respect to the physical origin of artiﬁcial atom and frequency range of atom-light interaction. The model is adapted to the 2D-heterostructures (THz frequencies), semiconductor quantum dots (optical range), and Josephson junctions (microwaves). The data for numerical simulations are taken from recently published experiments. The obtained results open a new way in quantum state engineering and nano-photonic spectroscopy. The model is universal with respect to the physical origin of artiﬁcial atoms and frequency ranges of atom-light interaction. The model was adapted to the semiconductor 2D-heterostructures (THz frequencies), semiconductor quantum dots (optical range), and Josephson junctions (microwaves). The initial data for numerical simulations are taken from recently published experiments. The dynamical equations have been studied both analytically and numerically. The idea for analytical simpliﬁcation is based on the separate description of electrons dressing with photons and quasi-classical motion of their geometrical center driven by a dc ﬁeld. The analytical solutions are in good qualitative agreement with numerical simulations. Our model is based on such conventional simpliﬁcations as neglect of damping and RWA. Their validity for considered systems was supported by numerical estimations taken from experimental data.


Introduction
The early quantum theory of electrical conductivity in crystal lattices by Bloch, Zener and Wannier [1][2][3][4] led to the prediction that a homogeneous dc field induces an oscillatory rather than uniform motion of the electrons. These so-called Bloch oscillations (BO) have been observed in bulk N-doped GaAs at a lattice temperature of 300 K in high fields up to 300 kV/cm [5], as well as different types of artificial systems such as semiconductor superlattices [6], interacting atoms in optical lattices [7,8], ultracold atoms [9][10][11][12][13][14], light intensity oscillations in waveguide arrays [15][16][17][18][19][20][21], acoustic waves in layered and elastic structures [22], and atomic oscillations in Bose-Einstein condensates [23] among others. Several recent studies have investigated the dynamics of cold atoms in optical lattices subject to ac forcing; the theoretically predicted renormalization of the tunneling amplitudes has been verified experimentally. The recent observations include global motion of the atom cloud, such as giant "super-Bloch oscillations" [24]. As a result, BO transformed from the specific contra intuitive model to the general experimentally-supported physical concept of oscillatory motion of wave packets placed in a periodic potential when driven by a constant force [8,25].
Rabi oscillations (RO) are periodical transitions of a two-level quantum system between its stationary states under the action of an AC driving field [26,27]. The phenomenon was theoretically The paper is organized as follows. In Section 2, we review the model and basic assumptions, formulate Hamiltonian at the Wannier-Fock basis and obtain equations of motion for probability amplitudes. In Section 3, we obtain the approximate analytical solution of equations of motion basing on quasi-classical concept. In Section 4, we describe and discuss the results of numerical calculations for electron Gaussian wave packets and different initial states of light (coherent state, vacuum field, double-Fock state entangled with Gaussian wave packet). In Section 5, we analyze the potential implementations of future experiments. The main results of the work and some promising tasks for future activity are formulated in Section 6.

Physical System and Model
Let us consider a one-dimensional (1D) structure of identical atoms placed over a line with period a (see for example where N is a number of atoms. We assume the tunneling to be the predominant mechanism of interatomic coupling and neglect other ones (such as Förster and the radiation field transfer). As it was shown in [41], such an assumption can be justified for a wide range of realistic parameter values. Before starting the consideration, we will discuss the potential degree of correspondence of our model to the different types of chains. Three rather conventional simplifications have been done: Before starting the consideration, we will discuss the potential degree of correspondence of our model to the different types of chains. Three rather conventional simplifications have been done: (i) We neglect all types of damping; (ii) we use so called rotating-wave approximation (RWA) [26]; and (iii) we assume the atomic chain to be infinitely long and perfectly periodic.
Assumption (i) relates to the case for scattering and radiation times strongly exceed the Bloch and Rabi periods. Because of RO for the quantum light are non-periodic, we imply here the value t R ∼ = Ω −1 n , as the Rabi-period ( n is the mean number of photons). These conditions are met with a large margin in superconductor junctions (artificial fluxonium atoms) [53][54][55][56]. The damping of BO and RO in heterostructures is defined by the electron scattering and decoherence, respectively. The recent progress in molecular beam epitaxy allowed achieving the values of scattering time τ ∼ = 10 −10 s in the ultra-high-quality AlGaAs/GaAs heterostructures [57]. For this case BO may be considered as ballistic, while dephasing becomes the dominant component of the damping for coherent inter-sub-band THz transitions (the typical values of characteristic times are 100-300 fs [50]). Such values are comparable with typical Bloch and Rabi periods; therefore, the damping does not manifest itself. In [50], the authors experimentally observed the manifestation of RO for the pulse with duration 200 fs, which does not exceed the dephasing time and is comparable with Rabi period (1-2 RO cycles over the pulse). For validity of our model to this case, we will consider interaction with rather short pulse, of which the duration, however, strongly exceeds the period of high-frequency filling and guaranties RWA validity.
The perfect periodicity means the identity and a rather large number of atoms in the chain. Its implementation is rather simple for 2D-heterostructures (for example, a chain of 51 element was used in experiments [50]). The problem is not so easy for QD-chains because of QDs comprise hundreds or thousands of real atoms, with inevitable variations in size and shape and, consequently, unavoidable variability in their energies and relaxation times [58,59]. One of the most promising technologies for applications in quantum photonics is the embedding of QDs within nanowires [59]. The QDs form at the apex of a GaAs/AlGaAs interface, are highly stable, and can be positioned with nanometer precision relative to the nanowire centre. As it was found [59], there is a chain of bright, nanoscale emitters in the red. QD-in-nanowire mimics very closely a two-level atom with high associated lifetime ( ∼ =450 ps). One more way is using a scanning tunneling microscope to create QDs with identical, deterministic sizes [58]. As it was mentioned in [58], the reached digital fidelity opens the door to QD architectures free of intrinsic broadening. This makes it reasonable to recommend our model of 1D QD-chain for visible light applications.
There exist different ways to theoretically describe complex quantum systems strongly interacted with quantum light. In quantum optics an "all-matter" picture is widely used, where the dynamics of light is integrated out (for example, in the optical Bloch equations [60,61]). In Refs. [62,63] the light-matter interaction is treated in an "all-light" picture (Lippmann-Schwinger equation approach). We use as the theoretical approach the probability amplitude method, generalized for the case of 1D chains driven via homogeneous dc field. It consists of solving of the Schrödinger equation for wave function, which is the superposition of various Wannier-Fock states for atom-light system.
To solve equations of motion for probability amplitudes analytically, we use the quasi-classical concept. The dressed electron is described as a wavepacket prepared with a well-defined quasi-momentum. The motion of quasiparticle center of mass is governed by the Newton's law, while the internal degrees of freedom have been described by means of quantum theory, based on the concept of electron-photon entanglement widely used in quantum optics [26,27]. Thus, the position of the quasi-particle center periodically evolves with a Bloch frequency corresponding to the quasi-momentum scanning a complete Brillouin zone. Such an approach was used for BO of conductive electrons in [8,64] and will be adopted for the case of dressed electrons in this paper.
In the general case, the equations of motion were integrated numerically with simplification through RWA. However, the use of the RWA may not describe the atom-light interaction when the coupling becomes sufficiently strong (ultrastrong coupling) [65]. Both lower and upper fundamental limitations of ultrastrong coupling of atoms and light were recently formulated in [66]. We bounded our consideration by the strong coupling regime assuming the lower limitation in [66] not to be reached.     The inductance L 1 is a tool of control of the "junction-resonator" coupling (it is assumed L 2 >> L 1 ).
Josephson energy E J strongly exceeds charging energy E C (E J >> E C ). We model the resonator as an LC-circuit with effective inductance L r and effective capacitance C r . The single fluxonium plays role of atomic chain due to the interaction with inductive environment [53][54][55][56] (BO take place in the space of quasi-charge).

Hamiltonian in Wannier-Fock Basis
Let us denote with a p and b p , the Wannier wave functions centered at p-th atom in the excited and ground states, respectively ( Figure 1). The two neighboring atoms are coupled via the electron tunneling, such that only intra-band tunnel transitions are permitted. It means that the electron due to the tunneling can go from state a p to state a p+1 and from state b p to state b p+1 only. The transitions between ground and excited states of different atoms are forbidden: a p±1 b p ≈ 0.
Let the chain be exposed to a single-mode quantum light, which electric field operator iŝ E =Ê (+) +Ê (−) = e √ ω/2ε 0 V(â +â + ), where V is the normalizing volume, e is the unit polarization vector,Ê (+) = e( ω/2ε 0 V) 1/2â ,Ê (−) = e( ω/2ε 0 V) 1/2â + are positive and negative frequency components of field operator,â,â + are creation-annihilation operators, respectively (an exp(−iωt) time dependence of the light is implicit). The field is assumed to be homogeneous over the chain axis. Such assumption corresponds to the chain excited by the normally incident wide laser beam, or specially-symmetric eigen-mode of microcavity, photonic crystal; etc. The chain is driven by an electrostatic (dc) field directed along the axis with responsibility for BO. We will consider the case of dipole interaction in the regime of strong coupling and assume the resonant condition ω 0 ≈ ω to be fulfilled. The system under consideration exhibits complex single-particle BO of the electrons entangled with photons, for which the theoretical framework will be introduced.
The system is described by the total Hamiltonian (see Appendix A): Here, the first term represents the free motion of the chain with absence of both light and dc field. It is given byĤ 0 = ( ω 0 /2)∑ pσzp , whereσ zp = a p a p − b p b p . The second termĤ ph = ωâ +â is the Hamiltonian of free electromagnetic field. The term is Hamiltonian of electron tunneling, t a,b are the penetration energies of potential barrier at the excited and ground states, respectively. The componentĤ describes the atom-light interaction, where g = − √ ω/2 ε 0 V(e · d ab ) is the interaction constant, d ab is the dipole moment, H.c. means Hermitian conjugation. The transition dipole moments in the chain is assumed to be vectors of identical values and orientations. The operatorsσ + p = a p b p , σ − p = b p a p are creation-annihilation operators of excited state in the p-th atom. The Hamiltonian (3) is written in the RWA form [26].
The last term,Ĥ describes the driving via dc field E dc and is responsible for BO.

Equations of Motion
The evolution of the system in the interaction picture is described by the Schrödinger equation i ∂ t |Ψ =V|Ψ with the interaction Hamiltonian given by: The state vector of the "atomic chain + light" system is: Here, b p , n = b p ⊗ |n , a p , n = a p ⊗ |n , are Fock states with n photons and Wannier states centered around p-th atom in ground and excited state, respectively, a p,n and b p,n are the unknown probability amplitudes. From the Schrodinger equation we get following equations for the probability amplitudes: (7) and (8) goes to the corresponding equations obtained in [41] for Rabi-waves in QD chains.

Studying Observable Values
We will study the values of two types. The first ones directly characterize the spatial-temporal behavior of electronic component and therefore are averaged over photonic state distribution. They are position-dependent, thus we will address to their densities per unit cell of the chain. From among these values we introduce: (i) inversion density, given by and (ii) tunneling current density (see Appendix B) The values of the second type directly characterize the quantum-statistical properties of light and therefore are spatially averaged. From among these values we will use: (iv) mean number of photons where p(n, t) is given by (11), (v) photonic number variance where n(t) is given by relation (12), and (vi) von Neumann quantum entropy where p j are the weights of the various Fock states of the statistical light distribution [44].

Preliminaries
The simplest and conventional model of BO is based on a quasi-classical approximation. It results from Equations (7) and (8) with g = 0. Electron-photon interaction disappears (subscript n becomes unnecessary) and the equations become separated. The system reads, and the same for b p . (the value δε becomes arbitrary and is taken δε = 0). In the absence of a dc field (E dc = 0), the eigenstate of the system corresponds to the electrons tunneling in the periodic potential and has a form of Bloch-wave a p+1 (t) = a p e i(ϕ−νt) , where ϕ is a phase shift per unit cell, coupled with the eigen frequency by ν = 2 −1 t a cos ϕ. The phase shift is coupled with the continuous quasi-momentum via relation h = ϕ/a, conventionally restricted to the first Brillouin zone [−π/a; π/a].
Under the influence of a dc field, a given Bloch-state |Ψ(h 0 ) = N −1/2 ∑ p a p e i(h 0 pa−ν(h 0 )t) evolves up to a phase factor into the state |Ψ(h(t)) with h(t) variation according to (16) or h(t) = −e −1 E dc t + h 0 . Thus, this evolution is periodic with a Bloch frequency corresponding to the time required for the quasi-momentum to scan a full Brillouin zone. It is described by the substitution The obtained solution relates also to the rather wide wave-packet with well-defined quasi-momentum. The periodic motion of its geometrical center corresponds to the quasi-classical model of BO. The main idea of the analytical solution considered next is based on hypothesis that the total dynamics of quasiparticles described by Equations (8) and (9) are reducible to the superposition of two interacting partial motions: (i) internal motion, which is dictated by dressing; and (ii) external motion of the quasiparticle as a whole. The internal motion is of completely quantum origin and doesn't allow any classical interpretation. It corresponds to the Rabi-wave and plays the role of the tunneling in the simplest case of BO considered before. The external motion is a motion of geometric center and will be described quasi-classically, following Equation (16).

Details
We consider the partial solutions of System (7) and (8) in the absence of a dc field (E dc = 0) and the case of zero detuning (ω = ω 0 ). We use the ansatz a p+1,n (t) = a p,n (t)e iϕ , b p+1,n (t) = b p,n (t)e iϕ , a p,n (t) = a pn e −iνt , b p,n (t) = b pn e −iνt , where ϕ is a given phase shift per unit cell, ν is the unknown eigen frequency, a pn , b pn are unknown constant coefficients. Making use of Equations (7) and (8) we obtain for them the matrix equation The eigenfrequencies are found from its characteristic equation and given by where g n = g √ n + 1. The correspondent eigenstates are where The States (14) and (15) describe the travelling of transitions between the states a p , n and b p , n + 1 along the chain (so called Rabi-waves [41]).
The periodicity of the lattice leads to a band structure of the energy spectrum of the dressed electrons in the entangled eigenstates |Ψ α,n (h) with the corresponding eigenenergies ν 1,2 (h, n). They are labeled by the discrete number of photons n and the continuous quasi-momentum h. Under the influence of a dc field, weak enough not to induce inter-band transitions, the state |Ψ α,n (h 0 ) evolves up to a phase factor into the state |Ψ α,n (h(t)) with h(t) variation according to Equation (16). This evolution is periodic with a Bloch frequency ω B corresponding to the time required for the quasi-momentum to scan a full Brillouin zone. It leads to the exchange, The distribution of photonic probabilities for the state |Ψ n,1 is given by The average number of photons is where the value ∆ n (t) is given by Relation (23). The similar elementary calculations for the state |Ψ n,2 give For photonic number variances we have Here, von Neumann entropy [44] accumulated at the given states of light is equal for both states and given by Due to BO, characteristics of quantum light, such as the photonic probabilities, the average number of photons, and the Neiman entropy are oscillatory functions with Bloch frequency ω B . The amplitude of these non-monochromatic oscillations is strongly dependent on the relation between penetrations of potential barrier in the ground and excited states. In particular, as obeyed from Equation (17), for the case of identical tunneling penetrations at the ground and excited states ( t a → t b ) the value ∆ n becomes constant, and oscillations vanish. The oscillating effect is suppressed likewise with the increasing of photons number n. As a consequence, a wave packet of entangled electron-photon prepared with a well-defined quasi-momentum will also oscillate in position with amplitude 2 (t a + t b )/eE dc .
Finally, we note that there is an interesting ability for quantum state control via adiabatically turning the dc field on and off. Starting from E dc = 0 we proceed to turn on E dc until its maximal value. Thereafter we slowly turn E dc off and end it at the moment t = T. This process corresponds to the turning the phase from initial value ϕ = ϕ 0 until the final one given by The optimal choice of the turning time T makes the arbitrary value of the phase turn reachable. As a result, the DC field via BO becomes an effective tool of control of quantum light statistics. In particular, the photonic probabilities, degree of electron-photon entanglement, and Neumann entropy may be adiabatically changed from minimal to maximal values (and vice versa) via adiabatic turning of the DC field.

Numerical Modelling and Discussion
In this section, we will show and discuss the numerical results for the different types of initial states. The System (8) and (9) was solved with the Crank-Nicolson numerical integration technique [67]. We have used periodic Born-Von Karman relations a 1,n = a N,n , b 1,n = b N,n as boundary conditions [41]. Let us note that the concrete form of boundary conditions has no physical importance in our case because the area of oscillations is placed rather distantly from the ends of the chain. For simplicity, we will limit our consideration by the case of zero detuning (ω = ω 0 ).
We will consider three types of initial conditions. For two of them the electron and photon subsystems are assumed to be initially non-interacting. The electronic subsystem is prepared as a coherent superposition of Gaussian wave packets: where a 0 , b 0 are arbitrary complex values that satisfy the normalization condition, u a,b , σ a,b , k a,b are the position of Gaussian center, effective Gaussian width, and the initial value of quasi-momentum at the ground and excited states, respectively. The photonic subsystem is prepared in the next types of state: (i) Coherent initial state; It is given by Poisson distribution p(n) = n n/2 e − n /2 / √ n!, where n is the mean number of photons. The total wavefunction at the moment t = 0 is given by (ii) Vacuum initial state; It is given by the total wave function where a 0 , b 0 are the arbitrary coefficients satisfying normalization condition; (iii) Entangled photon-electron initial state; The electron-photon entanglement prepared initially before the driving fields would be switching on. The initial condition reads where C 0 is normalization constant. The value of inversion for the State (34) is equal to zero.
In Figure 4, we plotted the temporal dynamics of inversion. The case of zero dc field and coherent initial light (Figure 4a) corresponds to the well-known scenario of collapse-revivals given by the Jaynes-Cummings model for the single atom [26]. The initial spatial distribution is not disturbed over the motion via inter-atomic tunneling. For the case of dc field existence and coherent initial light (Figure 4b), the collapse-revivals picture drifts due to BO: The collapses and revivals appear at the different segments of the atomic chain. The inversion behavior dramatically changes for the vacuum state of light (Figure 4c). Here, BO leads to the separation of areas of positive and negative values of inversion. The particle in the excited state (positive inversion) exhibits the single oscillation of rather high amplitude in the course a half of Bloch circle. During the other part of Bloch circle, the particle is in rest, which means its appearance in the ground state (negative inversion). The reason for this is the weak value of transparence of potential barrier at the ground level used in the shown calculations. In contrast with single atom scenario, vacuum RO between maximal and minimal values of inversion occurs with Bloch frequency ω B instead of Rabi frequency Ω 0 = 2g.
In Figure 5, we plotted the time dependence of the tunneling current. For the absence of light (Figure 5a), the current exhibits the periodic dynamics with Bloch frequency. This picture is dramatically changed due to the light-chain coupling both for coherent and vacuum states of light (Figure 5b,c), respectively. For the coherent light, the tunnel current is modulated in exact compliance with the collapse-revivals picture. This result is rather counterintuitive; it demonstrates the influence of the quantum-light statistics on the low-frequency motion of the charges driven separately by a dc field. It makes the spectra of the tunneling current more wide and various wherewith in the case of ordinary BO. It opens a new promising avenue in the spectroscopy of nano-circuits and nano-devices based on the synthesis of quantum optical and dc tools of control [51,52]. For the vacuum state of light, the tunnel current periodically evolves through one pulse after another with frequency ω B instead of continuous periodic behavior in the case of ordinary BO [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25]. Such evolution is accompanied by the modulation with Rabi frequency. The results demonstrate high and various mutual interactions of RO and BO. comparable with RO components. As a result, the mean number of photons is modulated with Bloch frequency B ω . Therefore, the dc field via BO opens one more way of engineering of the quantum light states. This effect depends on the large number of physical factors, such as dc field value, energy and dipole moment of quantum transition for the atom, tunnel coupling, etc. Concluding, the results of numerical modeling agree with the simple analytical model, developed in Section 3. They are promising for different applications in quantum computing, quantum informatics, etc.     In Figure 6, we plot the dynamics of mean number of photons. Again, we see the mutual influence of RO and BO. For dc field absence (Figure 6b), the collapse-revivals scenario corresponds the Jaynes-Cummings model for single atom [26]. It represents the quantum interference of the spectral components with different Rabi-frequencies Ω n and different relative amplitudes [26]. For dc field appearance (Figure 6a), there is one more type of interference components with frequencies mω B (m is integer value). Among them, the term with m = 1 is dominant, and its amplitude is comparable with RO components. As a result, the mean number of photons is modulated with Bloch frequency ω B . Therefore, the dc field via BO opens one more way of engineering of the quantum light states. This effect depends on the large number of physical factors, such as dc field value, energy and dipole moment of quantum transition for the atom, tunnel coupling, etc. Concluding, the results of numerical modeling agree with the simple analytical model, developed in Section 3. They are promising for different applications in quantum computing, quantum informatics, etc.      Figure 7b displays the evolution of photonic variance. As it was mentioned above, the electrons and photons in the initial state are uncorrelated. Thus, electron-photon entanglement is poor at the intervals of motion between the stopping points, whereby the variance oscillates only due to photon number fluctuations of rather small amplitude with vacuum Rabi-frequency. The processes of photons emission-absorption are accompanied the strong narrow peaks of variance. Thus, the vicinities of stopping points are characterized by the strong non-classicality of light and high degree of its entanglement with electronic sub-system. This scenario strongly differs from the case of the harmonic oscillations of variance with Rabi-frequency for the single atom [26]. The obtained results show that such manifestations of non-classicality of light, such as electron-photon entanglement and accumulated entropy become controllable by dc field via BO. As it is well known [26], the physical origin of vacuum RO is a manifestation of spontaneous emission of the single exited state in strong coupling regime. The obtained dynamics qualitatively support the conclusion of an early paper [68], that BO dramatically changes the physical picture of spontaneous emission in the atomic chains.  To illustrate the mutual influence of high-frequency and low-frequency motions, we show the spatial-temporal evolution of inversion for initial condition (34) (Figure 8). One finds out of the phase oscillations of inversion with Bloch frequency. The two halves of every Bloch period correspond to To illustrate the mutual influence of high-frequency and low-frequency motions, we show the spatial-temporal evolution of inversion for initial condition (34) (Figure 8). One finds out of the phase oscillations of inversion with Bloch frequency. The two halves of every Bloch period correspond to the excited and ground state, respectively. Here, BO of inversion is modulated with RO. Such dynamics are similar to the shown for the vacuum initial state at Figure 4c, however, the modulation depth strongly increases due to the influence of initial entanglement (the inversion evolves over the single Rabi-circle between zero and ±1). In Figure 9, we show the tunnel current for initial state (34), which as well as inversion exhibits the high level of RO modulation. Comparing it with Figure 5c, one finds again the strong influence of initial electron-photon correlations. Figure 10 displays the quantum-statistical properties of light. The evolution of mean number of photons (Figure 10a) may be considered as a periodic system of step-like beatings (dashed line), modulated with high-frequency oscillations (solid curve). The value n averaged with respect to the step-period is approximately equal to 0.5, which is equal to the main photon number for both states |Ψ n1,2 given by Equations (20) and (21). The general physical interpretation of Figure 10a may be done from the quasi-classical point of view as an interference of this states (of course, we speak here about moving wavepackets, but not about perfect travelling waves). Again, the photon emission-absorption with Bloch-period takes place in the vicinity of stopping points. One can see, comparing Figure 10a with Figure 9, that the moments of shut-down are agreed with zeros of inversion. Translatory motion of the wavepacket is accompanied by RO between eigenstates, which correspond to high-frequency fluctuations with Rabi-frequency. In Figure 10b we plot the dynamics of photonic variance. In contrast with vacuum photonic state, the variance is maximal at the intervals of translatory motion and strongly decreases in the vicinities of stopping points. The reason for it is that translatory motion corresponds to the statement in one of the Eigenstates (20) and (21). Every one of these states is characterized by the maximal level of electron-photon entanglement. Photons emission-absorption means the mutual transformation of the States (20) and (21). It leads to creation the coherent superposition of these states for a short time with comparable probability amplitudes. It leads to the breaking of electron-photon correlations. Again, RO-fluctuations are imposed to this ideal picture similar to the case of the vacuum photonic state. The photonic variance in these oscillations reaches the value δn = 0.25, which corresponds to the maximal one for two States (20) and (21). emission-absorption with Bloch-period takes place in the vicinity of stopping points. One can see, comparing Figure 10a with Figure 9, that the moments of shut-down are agreed with zeros of inversion. Translatory motion of the wavepacket is accompanied by RO between eigenstates, which correspond to high-frequency fluctuations with Rabi-frequency. In Figure 10b we plot the dynamics of photonic variance. In contrast with vacuum photonic state, the variance is maximal at the intervals of translatory motion and strongly decreases in the vicinities of stopping points. The reason for it is that translatory motion corresponds to the statement in one of the Eigenstates (20) and (21). Every one of these states is characterized by the maximal level of electron-photon entanglement. Photons emission-absorption means the mutual transformation of the States (20) and (21). It leads to creation the coherent superposition of these states for a short time with comparable probability amplitudes. It leads to the breaking of electron-photon correlations. Again, RO-fluctuations are imposed to this ideal picture similar to the case of the vacuum photonic state. The photonic variance in these oscillations reaches the value 0.25 n δ = , which corresponds to the maximal one for two States (20) and (21).  t ω = .

Semiconductor 2D-Heterostructures
The next promising candidate for the experimental implementation of the model considered above is the low-dimensional semiconducting heterostructure (see Figure 2) under femtosecond inter-sub-band excitation [50]. For example, the sample consisting of 51 GaAs quantum wells

Semiconductor 2D-Heterostructures
The next promising candidate for the experimental implementation of the model considered above is the low-dimensional semiconducting heterostructure (see Figure 2) under femtosecond inter-sub-band excitation [50]. For example, the sample consisting of 51 GaAs quantum wells separated by barriers was implemented in [50]. A coherent excitation of the sample was created by a femtosecond pulse with a center frequency resonant to the inter-sub-band transition. The resonant line at 30 THz was homogeneously broadened with coherence time 320 Fs. The amplitude of incident field was varied inside of area (5-50) kV/cm. As a result, there were observed coherent sub-picosecond RO and manipulated in a wide range by varying the strength of the coherent driving field. The measurements [50] qualitatively agreed with the predictions of the simplest model based on Maxwell-Bloch equations for non-interactive two-level systems. Conventional simplifications, such as RWA and omitting of all types of damping have been used. The simulation of the pulse driving field was done by varying the Rabi frequency Ω R = Ω R (t).
Here, we add the dc field for manipulation by atomic chain via BO. For the application of our model to the potential experiments with 2D-heterostructures, it is necessary to describe the driving quantum light as a short transient normally incident to the heterostructure. The driving process in this case is described by the instantaneous coupling coefficient g. The field quantization should be modified following Appendix C. The quantum properties of light in the wavepacket are dictated by the special pair of bosonic creation-annihilation operators, which allows us to rewrite the dynamic Equations (7) and (8) as: i ∂a p,n ∂t = δε − E dc · eR p a p,n + t a a p+1,n + t * a a p−1,n − g(t) where g(t) = d ab · u 0 (t)/ , u 0 (t) is a slow envelope of the driving pulse. The driving with electromagnetic pulse leads to chain dynamics qualitatively different from the case of monochromatic light. Figures 11 and 12 depict the temporal-spatial evolution of inversion for different relations between barrier penetrations at 1S and 2S levels. Figure 11 shows the case of identical penetrations. One finds BO over the pulse duration. The oscillations are modulated with Rabi frequency and stop with the pulse disappearing at zero inversion. The behavior becomes quite different for the case of 2S penetration strongly exceed 1S one. (Figure 12). No oscillations exist before the pulse appearance (t < 0.9). The reason for this is a weak value of penetration at 1S level. The switching of pulse creates both RO and BO. The RO existence is limited by the pulse duration (t < 1.2), while BO continues even after the pulse switches off. It is a result of tunneling at the 2S-level on which the penetration is much higher. The charge existence at the 2S-level is a result of resonant pumping by the optical pulse. Thus, we can consider the resonant pumping as a tool of BO triggering and speak about photon-assisted BO.
( ) where ( ) ( ) is a slow envelope of the driving pulse.
The driving with electromagnetic pulse leads to chain dynamics qualitatively different from the case of monochromatic light. Figures 11 and 12 depict the temporal-spatial evolution of inversion for different relations between barrier penetrations at 1S and 2S levels. Figure 11 shows the case of identical penetrations. One finds BO over the pulse duration. The oscillations are modulated with Rabi frequency and stop with the pulse disappearing at zero inversion. The behavior becomes quite different for the case of 2S penetration strongly exceed 1S one. (Figure 12). No oscillations exist before the pulse appearance ( 0.9 t < ). The reason for this is a weak value of penetration at 1S level. The switching of pulse creates both RO and BO. The RO existence is limited by the pulse duration ( 1 .2 t < ), while BO continues even after the pulse switches off. It is a result of tunneling at the 2Slevel on which the penetration is much higher. The charge existence at the 2S-level is a result of resonant pumping by the optical pulse. Thus, we can consider the resonant pumping as a tool of BO triggering and speak about photon-assisted BO.

Josephson Junction
The next promising candidate for experiment implementation is a single Josephson junction embedded in an inductive environment. As it was shown in [53][54][55][56], the voltage-biased Josephson junction exhibits the temporal dynamics equal to the motion of a fictitious charged particle with

Josephson Junction
The next promising candidate for experiment implementation is a single Josephson junction embedded in an inductive environment. As it was shown in [53][54][55][56], the voltage-biased Josephson junction exhibits the temporal dynamics equal to the motion of a fictitious charged particle with inertia provided by the inductance. This particle moves in a potential, which also has a periodic part and a linear tilt. On the other hand, a Josephson junction strongly coupled to the microwave resonator ( Figure 3) in the resonant regime may be considered as an artificial two-level atom ("fluxonium" [54]). It is described by the Jaynes-Cummings model [54] and exhibits RO dynamics [54]. Here, we will focus on the case of Josephson junction simultaneously driven by dc voltage and quantum light in the regime of small detuning. We will show that the model investigated in this work is applicable to this type of physical system. Therefore, we will start our analysis by identifying the parameters of Josephson junction in terms of this model.
The temporal dynamics of entire circuit may be described using so-called quasi-charge representation [53,55]. Let us first consider the single junction, of which the standard form readŝ H 0 = Q 2 /2C J − E J cosφ J , whereQ = (2e/i)∂/∂φ J ,φ J are operators of junction charge and phase difference across the junction, respectively. It describes the particle of mass C J moving in the periodic potential [55] (C J is a junction capacitance). The eigenfunctions of this Hamiltonian are Floquet-Bloch modes with condition of quasi-periodicity Ψ α,Q (ϕ) = P α,Q (ϕ)e iϕQ/2e , where P α,Q (ϕ) = P α,Q (ϕ + 2π) is periodic function, α is a number of the gap. These functions describe the motion of the quasiparticles, which is similar to the motion of electrons in crystal. The phase ϕ and quasi-charge Q play the role of the conjugated spatial coordinate and quasi-momentum for electrons in crystals, respectively. The energies of the two first modes ε 0,1 (Q), are separated by the gap, which in the regime of interest here (E J >> E C J , E C J = e 2 /2C J ) is of the order of plasma frequency ω p = 8E J E C J / . Taking E J /h = 30 GHz, and E C J /h = 3, we have ω p = 150 GHz. Next, in the limit E J >> E C J the energy bands ε 0,1 (Q) are purely sinusoidal [53] with the bandwidths of the order of [55] (the typical values are ∆ 0 ≈ E C /2 = 1, 5 GHz). The situation is the same in the case of atomic chain considered above (see Equation (13) in the case g n = 0). Thus, we are able to establish the correspondence 2t a,b ↔ ∆ 0 .
The additional series inductance and switching on the bias voltage transforms the Hamiltonian H 0 toĤ where E L J = ( /2e) 2 L −1 2 is the constant of inductance energy, φ b is the total phase difference across the circuit, which is determined by the bias voltage V b , such that .
The case under consideration has a simple physical interpretation based on the charge-phase duality [53]. The voltage-biased junction with inductive environment and the current-biased junction shunted by a capacitor C J are dual to each other when exchanging the role of quasi-charge and phase. The dynamics of the circuit, shown in Figure 3, are equal to the intra-band motion of a fictitious charged particle in the tilted periodic potential. The tilt leads to BO of the particle in the quasi-charge space, which allows for the identification of the Bloch-frequency as dual plasma frequency The inductive strong coupling of a Josephson junction with a microwave resonator leads to its dressing with microwave photons. It corresponds to the periodic emission-absorption of photons and inter-band transitions by the fictitious charged particle (RO). Thus, the qualitative picture of junction behavior is similar to the BO of dressed electrons considered above: The periodic inter-band transitions (RO of fluxonium) are accompanied by the charge oscillations through the bias-voltage driving. For formal identification of the coupling parameters of the system, "junction + microwave resonator", it is necessary to analyze in more detail the total Hamiltonian. It readŝ H =Ĥ A +Ĥ ph +Ĥ I,ph +Ĥ I,dc , whereĤ A is a Hamiltonian of the Josephson junction in an inductive environment,Ĥ ph is Hamiltonian of microwave photons in the resonator,Ĥ I,ph is a Hamiltonian of junction-cavity interaction, andĤ I,dc is a support of biased voltage. For further analysis, it is convenient to use the secondary quantization technique. The Hamiltonian of photons readŝ whereQ r = i 2Z r (â −â + ),Φ r = i Z r 2 (â +â + ),â,â + are creation-annihilation operators for microwave photons, which satisfy the conventional bosonic commutative relations. By following the concept of quasi-charge-phase duality, the Fourier-Bloch states Ψ α,Q (ϕ) exchanged by the Wannier functions (see Appendix A) and the Hamiltonian is written in the Wannier-basis. We denote the two first states with quasi-charge Q as e Q and g Q and introduce the raising operatorsσ Q = e Q g Q andσ + Q = g Q e Q . The atomic Hamiltonian may be rewritten asĤ A = ω p∑QσQz /2, whereσ Qz is Pauli inversion matrix.
For junction-cavity interaction we have [54] H where G ≈ ( /2e) 2 L 1 (L 2 L r ) −1 , φ eg = g Q σ + Q e Q is the matrix element of effective dipole moment. The Hamiltonian is written in rotating-wave approximation. It is equal to the Hamiltonian (3) if the coupling coefficient identified with the value where R 0 is the quantum of resistance. The validity of rotating-wave approximation implies all relevant energies to be smaller than the plasmonic frequency of the Josephson junction and resonant frequency of resonator. In particular, we impose the characteristic energy, E L J , associated with the inductance L 2 to be smaller than ω p . It corresponds to the values L 2 ≈ 200 nH, which are quite reachable [56]. Second, we need also g << ω p . For standard value Z r = 50 kΩ [56] we have L 1 /L 2 ≈ 10 −3 for this condition been imposed. The coherence time and energy relaxation time both are of the order ∼ = µs [56]. Thus, the number of BO flops is high enough for their experimental observation in the fluxonium atom in microwave frequency range. The real cases may be described by lossless model considered above.

Conclusion and Outlook
In this paper, we developed a model of the chain of two-level artificial atoms manipulated simultaneously by a dc field and single-mode quantum light in the strong coupling regime. Atom-light interactions were assumed to be resonant. The cases of monochromatic light and light pulse have been considered. Both, electronic and photonic characteristics of the oscillation process, have been studied, such as: (i) inversion density; (ii) tunneling current density; (iii) distribution of photon probabilities; (iv) mean number of photons; (v) photon number variance; and (vi) Neumann entropy of light. The different types of the initial states of the light have been considered: (i) coherent state; (ii) vacuum photonic state; and (iii) Fock-state entangled with the electron wave packet. The Gaussian wavepackets were chosen as initial states of atoms.
The model is universal with respect to the physical origin of artificial atoms and frequency ranges of atom-light interaction. The model was adapted to the semiconductor 2D-heterostructures (THz frequencies), semiconductor quantum dots (optical range), and Josephson junctions (microwaves). The initial data for numerical simulations are taken from recently published experiments.
The dynamical equations have been studied both analytically and numerically. The idea for analytical simplification is based on the separate description of electrons dressing with photons and quasi-classical motion of their geometrical center driven by a dc field. The analytical solutions are in good qualitative agreement with numerical simulations. Our model is based on such conventional simplifications as neglect of damping and RWA. Their validity for considered systems was supported by numerical estimations taken from experimental data.
The following conclusions are emerged from our studies: (1) The case of initial coherent state of the light exhibits the collapse-revival picture, which drifts over the chain (collapses and revivals placed in the different spatial areas). In contrast with Jaynes-Cummings model, the collapse-revival picture is modulated with Bloch frequency; (2) In the case of initial vacuum state of light, the photon emission and absorption occurs with Bloch frequency, instead of Rabi-frequency in Jaynes-Cummings model. The photonic probabilities are mainly modulated with Bloch frequency, while the contribution of Rabi-components is rather slight. BO strongly squeezes the vacuum state of light entangled with electronic wavepackets; (3) The electron-photon entanglement dramatically modifies the tunnel current behavior. It becomes modulated agreed with the collapse-revival picture for the case of coherent state of light, and periodically modulated by RO for the case of initial photonic vacuum; To conclude, the main result of the paper is the novel effect of the influence of BO and quantum statistical properties of light on each other. It is counterintuitive because of the strongly different frequency ranges for such types of oscillations existence. The reason for it is an entanglement of electronic and photonic states in the system been considered.
The obtained results allow for the control of quantum-statistical properties of light via adiabatically turning the dc field. They are promising for applications in quantum optics, quantum informatics, and quantum computing. The mutual modulation of low-frequency BO and high-frequency inter-band transitions produces the new types of spectral lines in tunneling current and optical polarization. It opens the promising ways in spectroscopy of nanodevices in THz and optical frequency ranges. These problems may be considered as directions for future research activity.
Author Contributions: Developments of the physical models, derivation of the basis equations, interpretation of the physical results and righting the paper have been done by I.L. and G.S. jointly. The Numerical Python calculations and Matplotlib figures were produced by I.L.