Full Counting Statistics of Electrons through Interaction of the Single Quantum Dot System with the Optical Field

In this paper, using the particle-number-resolved master equation, the properties of full counting statistics (FCS) are investigated for a single quantum dot (QD) system interacting with optical fields in the thermal state, Fock state, coherent state, and coherent state with random phase. In these diverse quantum states of optical fields, average tunneling currents have different step shoulder heights at a lower bias voltage with the same light intensity, and a staircase-shaped current can be induced unexpectedly in vacuum state optical field. The characteristics of the Fano factor and skewness in the coherent state differ from those in all of the other cases. For avalanche-like transport at a lower bias voltage, the mechanism is a dynamical channel blockade in a moderate electron–photon interaction regime. There is a pronounced negative differential conductance that results from tuning the phase of the coherent state optical field in a symmetric QD system.


Introduction
The full counting statistics (FCS) of charge transfer in quantum dots (QDs), molecules, and nanostructures are an important component in the characterization of microscopic transport processes because the FCS contains complete information about low-frequency current fluctuations [1]. Theoretically, FCS can provide all higher-order moments of current fluctuations and offer deep insight into the nature of transport mechanisms not encoded in the first two moments [2]. In particular, the system parameters dominating the transport can be extracted from current-noise measurements. Experimentally, real-time counting statistics have been calculated for single-electron tunneling transport through QDs, and this is a crucial achievement that enables counting individual electron tunnel events [3]. Until now, FCS have been used to study different systems, such as a mesoscopic superconductor [4], a voltage-biased superconducting nano-junction [5], multichannel chaotic cavities [6], a ballistic chaotic cavity [7], a coupled QD system [8], single-molecule junctions [9], and atomic spin devices [10].
On the other hand, owing to the versatility of nanofabricated circuits, hybrid circuit quantum electrodynamics (QED) can provide numerous diverse situations that are not accessible with the standard cavity QED [11]. The advantage of hybrid circuit QED systems is that the properties of a nanostructure device can be arbitrarily controlled by an externally applied field. QD systems can be used in circuit QED to probe light-matter interactions [12][13][14][15], implement a QD laser [16,17], and engineer new states of matter relevant to solid-state physics [18]. In addition to that, the QD is one of the candidates for use as a quantum information device [19].
In an effort to construct quantum optical and quantum electronic QD devices, photon-assisted electronic transport has been investigated by many scientists [20]. However, to study the current and shot noise, most of the studies have employed the classical treatment of the external field, which introduces a time-dependent oscillating energy level in the QD [20,21]. Furthermore, the results are the same if the intensities of the different external optical fields are the same. The classical treatment for photons is valid at high field intensity and for the case of weak coupling in an equilibrium system, but it does not work in the vacuum state or other quantum states of light [22]. Moreover, there are many papers that discuss the FCS properties of a single-molecule junction, which has a similar Hamiltonian (QD interacting with a phonon) to that presented in this paper [23,24]. In previously published papers, the phonon distribution is either in equilibrium or nonequilibrium. The phonon distribution is a dominant factor in determining the transport properties. Usually, a molecule attached to its substrate exhibits an equilibrium phonon distribution while the phonons of a molecule air-bridged to electrodes have a nonequilibrium distribution [25]. The equilibrium state of the aforementioned phonon is a thermal state. In this paper, we apply a Hamiltonian with the quantum treatment of an optical field to study the influences of different quantum states of optical fields (thermal state, Fock state, coherent state, and coherent state with random phase) on the FCS properties of a QD in which the photon distribution is in equilibrium.
The method used in this paper is the particle-number-resolved master equation [26,27]. An alternative method is to use rate equations to derive the tunneling rate by Fermi's golden rule, in which the term that describes the relaxation of vibrations toward the equilibrium distribution is added by hand [23,24,28]. In rate equations, all off-diagonal components are ignored because they decay rapidly in certain parameter regimes [29].
This paper is organized as follows. In Section 2, we present the physical model and theoretical research by the particle-number-resolved master equation approach. In Section 3, we study the influences of optical fields with different statistical properties on the average current, Fano factor, and skewness of single-electron tunneling for different parameters. The conclusion is given in Section 4. Figure 1 illustrates the system schematic of the single QD with normal electrode leads coupled to a one-mode optical cavity. The QD is modeled as a one-level system. The loss of the optical field is not considered here because our focus is on the influence of different optical fields on the transport through a single QD.

Physical Model and Formalism
The total Hamiltonian of the system can be written as H = H L + H R + H ph + H D + H T . The Hamiltonian for electrons in the left (L) and right (R) electrode leads are (e = = 1) [22][23][24] where c † k,α (c k,α ) is the conduction electron creation (annihilation) operator with wave vector k in electrode lead α, and ε k,α is the single-electron energy. The third term H ph stands for the single-mode optical field, and H ph = ω 0 a † a, where a † (a) is the photon creation (annihilation) operator with frequency ω 0 . The electron Hamiltonian on the QD is where ε denotes the energy level of the dot and can be controlled by modulating the gate voltage V g , and λ is the coupling constant between the dot electron and photon mode. d † (d) is the corresponding creation (annihilation) of an electron on the QD, and we assume the on-site Coulomb repulsion on the QD is infinite. The last term H T describes the conduction electron hopping between the QD and electrode leads, where T k,α is a tunneling matrix element.
aforementioned phonon is a thermal state. In this paper, we apply a Hamiltonian with the quantum treatment of an optical field to study the influences of different quantum states of optical fields (thermal state, Fock state, coherent state, and coherent state with random phase) on the FCS properties of a QD in which the photon distribution is in equilibrium. The method used in this paper is the particle-number-resolved master equation [26,27]. An alternative method is to use rate equations to derive the tunneling rate by Fermi's golden rule, in which the term that describes the relaxation of vibrations toward the equilibrium distribution is added by hand [23,24,28]. In rate equations, all off-diagonal components are ignored because they decay rapidly in certain parameter regimes [29]. This paper is organized as follows. In Section 2, we present the physical model and theoretical research by the particle-number-resolved master equation approach. In Section 3, we study the influences of optical fields with different statistical properties on the average current, Fano factor, and skewness of single-electron tunneling for different parameters. The conclusion is given in Section 4.  Figure 1 illustrates the system schematic of the single QD with normal electrode leads coupled to a one-mode optical cavity. The QD is modeled as a one-level system. The loss of the optical field is not considered here because our focus is on the influence of different optical fields on the transport through a single QD.

Physical Model and Formalism
The total Hamiltonian of the system can be written as H = HL + HR + Hph + HD + HT. The Hamiltonian for electrons in the left (L) and right (R) electrode leads are (e = ℏ = 1) [22][23][24]  It is convenient to eliminate the electron-photon coupling terms in the Hamiltonian by a canonical transformation, i.e., H = e s He −s with s = (a † − a)d † d and = λ/ω 0 , so the transformed Hamiltonian becomes H = H el + H sys , where H sys = ω 0 a † a + εd † d, and To investigate an open quantum system, such as a single QD interacting with an optical field and coupled to a reservoir of electrons in the left and right electrode leads, it is appropriate to employ the reduced density matrix of this system ρ, which is obtained by taking the trace over the degrees of freedom of the reservoirs and optical field [27,30], Then, we can get the Born-Markov master equation as where the notation D is defined as D[A]ρ = AρA † − 1 2 A † Aρ + ρA † A , and the superscript R and I represent the real part and imaginary part, respectively. W i (i = 1, 2, 3, 4) represents the tunneling rate between the right (left) leads and the QD [27,30].
where Γ α = 2πρ α ∑ k |T k,α | 2 , f α is the Fermi distribution function for lead α, and F >(<) (ω) is the Fourier transform of the bosonic correlation function X † (t)X ( X(t)X † ). We employ the standard decoupling approximation, which assumes that the averages of the products of two operators ( X † (t)X and X(t)X † ) can be evaluated in the equilibrium state. This approximation, called the single-particle approximation [31], is valid if ω 0 or λ is large compared with the QD-lead coupling [32]. After expanding the exponents in a power series and applying the Cauchy product to the power series, we can obtain [33,34] for an optical field in a thermal state, which is a field emitted by a classical light source in thermal equilibrium at temperature T [35]; |N is the eigenstate of the photon particle number operator a † a.
for an optical field in a Fock state, |N Fock , which can be generated by a cavity quantum electrodynamics scheme using the strong coupling between excited atoms and a single-mode cavity field [36]. The number of photons is exactly defined for a Fock state.
for an optical field in a coherent state, |α , which is the eigenstate of the photon annihilation operator a and has a more precisely defined phase. A coherent state is the closest quantum approximation to the field generated by a laser. A coherent state does not contain a definite number of photons, and this follows from the fact that a coherent state is made up of a coherent superposition of Fock states, |α = ∑ ∞ n=0 |n n|α . Here, I n is the nth modified Bessel function, and J n is the nth first kind of Bessel function.
if the optical field is in a coherent state with random phase, ρ opt = e −|α| 2 ∑ ∞ n=0 |α| 2n n! |n n| [37]. Here, L 2l−k k represents Laguerre polynomials. If we stay within the coherent length of the laser and work at frequencies for which the laser's intensity noise is at a quantum noise limit, then a laser beam can represent an excellent realization of an ideal coherent state [38]. Otherwise, it represents a random phase coherent state. Because we assume that the Coulomb repulsion on a QD is infinite, there are only two states for a QD: an unoccupied state and an occupied state. Denoting 0|ρ|0 and 1|ρ|1 as ρ 0 and ρ 1 , respectively, then, using Equation (6), one can obtain the particle-number-resolved master equation, Then, the particle-number-resolved master equation can be rewritten as where In order to facilitate the calculation of FCS, we introduce the following matrix S(χ, t), which is defined by S(χ, t) = ∑ n ρ n (t)e inχ , where χ is the counting field [21]. Since the particle-number-resolved master equation has the form given in Equation (13), S(χ, t) satisfies the equation where The cumulant generating function (CGF) is defined as e −F(χ) = Tr[S(χ, t)]. Thus, under the low-frequency limit, F(χ) can be written as F(χ) = −λ(χ)t, where λ(χ) is the eigenvalue of L χ and tends to zero ( λ(χ) → 0 ) as χ → 0 [26]. Using Equation (15), one can get which can be substituted into After some algebra, one can obtain the first three cumulants, C 1 , C 2 , and C 3 , as The first-order cumulant, defined as the average of the electron number distribution C 1 = n(t), is clearly related to the average current, I = eC 1 /t. The second-order cumulant C 2 = n 2 − n 2 is the mean-square deviation and is related to the zero-frequency shot noise, defined as S = 2e 2 C 2 /t. The third cumulant C 3 = (n − n) 3 characterizes the skewness of the distribution. The overbar denotes the statistical average, that is (· · · ) = ∑ n (· · · )P (n, t). In general, the normalized second-order cumulant F a = C 2 /C 1 is called the Fano factor, which is the shot noise, and the third cumulant S k = C 3 /C 1 is used to describe the skewness. The electron-number-resolved statistic is called the super-Poissonian for F a > 1 and sub-Poissonian for F a < 1.

Discussion
In the following discussion, we set the effective energy level to ε = E 0 = E F + V g − δ, with E F = 0, which is the equilibrium chemical potential of the left (right) electrode lead in the absence of gate voltage V g ; µ L (µ R ) is the chemical potential of the left (right) electrode lead, which is taken as µ L(R) = E F ± V bias /2; and V bias is the bias voltage between two electrode leads. We also assume the tunneling rate Γ α as Γ l = Γ r = Γ 0 = 0.1, with the photon frequency ω 0 = 1 (Γ 0 ω 0 ). The energy unit is taken as the photon frequency ω 0 throughout this paper. Now, we begin to discuss the influences of different types of optical fields on the FCS properties of charge transport on QDs. Figure 2a shows the average tunneling current as a function of V bias in the regime of moderate electron-photon coupling (ζ < 1). In Figure 2a, the well-known staircase-shaped current appears for different states of optical fields. As V bias increases to cover the photon sidebands, the photon-assisted resonant tunneling channels open one by one, leading to an increasing number of current steps with decreasing shoulder heights, as presented in Figure 3. This is similar to the results of the electron-phonon coupling case [9,[23][24][25] or the AC-driven nanostructure system [20,21]. The average tunneling current has turning points at bias voltages of V bias = 2nω 0 because the photon sidebands locate at nω 0 (n = ±1, ±2, . . . . . .), which is understandable from Equations (7)- (11). When the light intensity is the same for all of the quantum states of optical fields, the average tunneling currents have the same saturation current at a higher bias voltage V bias , but they have different shoulder heights at a lower bias voltage V bias in Figure 2. Our result is different from that obtained by the classical treatment of the photon field in which currents are the same if the intensities of different external optical fields are identical [20,21]. rate Γ as Γ = = Γ ＝0.1, with the photon frequency = 1 (Γ ≪ ). The energy unit is taken as the photon frequency throughout this paper. Now, we begin to discuss the influences of different types of optical fields on the FCS properties of charge transport on QDs. Figure 2(a) shows the average tunneling current as a function of Vbias in the regime of moderate electron-photon coupling ( < 1). In Figure 2(a), the well-known staircase-shaped current appears for different states of optical fields. As Vbias increases to cover the photon sidebands, the photon-assisted resonant tunneling channels open one by one, leading to an increasing number of current steps with decreasing shoulder heights, as presented in Figure 3. This is similar to the results of the electronphonon coupling case [9,[23][24][25] or the AC-driven nanostructure system [20,21]. The average tunneling current has turning points at bias voltages of Vbias = 2n because the photon sidebands locate at ( = ±1, ±2, … … ), which is understandable from Equations (7)(8)(9)(10)(11). When the light intensity is the same for all of the quantum states of optical fields, the average tunneling currents have the same saturation current at a higher bias voltage Vbias, but they have different shoulder heights at a lower bias voltage Vbias in Figure 2. Our result is different from that obtained by the classical treatment of the photon field in which currents are the same if the intensities of different external optical fields are identical [20,21]. Furthermore, we find that the vacuum state optical field also induces a staircase-shaped average current (cyan dotted line). That means that the photon sidebands still work in the quantum case because of the interaction of the QD's electron with the vacuum fluctuations, which is undiscovered in the case of an AC-driven nanostructure system [20,21]. Comparing the curve of the vacuum state with those of the other states, we can observe that the average current is markedly suppressed as the photon intensity increases from zero to one. Figure 2b demonstrates that the Fano factor-voltage characteristics exhibit a singularity at zero bias voltage, which corresponds to the finite thermal Nyquist-Johnson noise due to the thermal fluctuations in the occupation number in the left and the right leads [2]. It can be seen that the Fano factor in the vacuum state (cyan dotted line) decreases quickly to 0.5 with an increasing bias voltage Vbias. This limit corresponds to the regime of the Pauli blockade and Coulomb blockade, in which electron transport is anticorrelated due to the Coulomb repulsion and the Pauli exclusion principle [39,40]. The Fano factor-voltage characteristics for the thermal state, Fock state, and coherent state Furthermore, we find that the vacuum state optical field also induces a staircase-shaped average current (cyan dotted line). That means that the photon sidebands still work in the quantum case because of the interaction of the QD's electron with the vacuum fluctuations, which is undiscovered in the case of an AC-driven nanostructure system [20,21]. Comparing the curve of the vacuum state with those of the other states, we can observe that the average current is markedly suppressed as the photon intensity increases from zero to one. Figure 2b demonstrates that the Fano factor-voltage characteristics exhibit a singularity at zero bias voltage, which corresponds to the finite thermal Nyquist-Johnson noise due to the thermal fluctuations in the occupation number in the left and the right leads [2]. It can be seen that the Fano factor in the vacuum state (cyan dotted line) decreases quickly to 0.5 with an increasing bias voltage V bias . This limit corresponds to the regime of the Pauli blockade and Coulomb blockade, in which electron transport is anticorrelated due to the Coulomb repulsion and the Pauli exclusion principle [39,40]. The Fano factor-voltage characteristics for the thermal state, Fock state, and coherent state with random phase are almost the same and decrease to 0.5 with two steps with increasing bias voltage V bias . However, for the coherent state, as the bias voltage V bias increases, F a decreases to about 0.3, then rises abruptly to about 0.6, and finally decreases to 0.5. This pattern is different from those in other optical fields. optical field, which is similar to that of the QD system interacting with a phonon [42]. Once an electron jumps to a QD from an electrode lead, it can jump on the photon sideband outside the bias voltage window by absorbing a photon. Because of the Coulomb blockade and Pauli blockade, no other electrons can jump to the QD, and this leads to the dynamical channel blockade due to the small tunneling rate of the occupied channel. As the bias voltage increases, the dynamical channel blockade can be gradually removed if all the photon sidebands enter the bias voltage window. Thus, the Fano factor will walk into the regime dominated by the Pauli blockade and Coulomb blockade from the regime dominated by the dynamical channel blockade with an increasing bias voltage.
From Figure 2b, it can also be found that the Fano factor for the coherent state is less than that of the vacuum state for a small bias voltage. The reason is that F and F (Equation 10) contain contributions from the off-diagonal terms of 〈 | ( ) | 〉 and 〈 | ( )| 〉 ( ≠ ) , respectively, in the photon number representation. This will make the amplitude of the photon sideband of the coherent state fall below that of the zero-field vacuum contribution for certain ranges of parameters. Thus, the coherent state could suppress the tunneling contribution of the vacuum state for certain ranges of parameters [22]. Figure 2c exhibits the skewness-voltage characteristics for the vacuum state, thermal state, Fock state, coherent state, and coherent state with random phase. The Sk for the coherent state decreases to about 0.7, then steps down to a constant value that is the same as that for the other quantum states with an increasing bias voltage Vbias, but it remains a constant value for the other three quantum states as the bias voltage Vbias increases. In order to explain the skewness-voltage characteristics, the curve of − versus the bias voltage of a steady solution is plotted for all four quantum states. The solid dark curves correspond to the same parameters as in Figure 2. Comparing the corresponding curves in Figure 2c and Figure 4, there seems to be a correlation between skewness and − .  The staircase-shaped Fano factor in Figure 2b can be interpreted by the mechanism termed the dynamical channel blockade [41]. In order to explain the principle, we provide a schematic diagram in Figure 3, which gives the density of the state for a QD's electron coupling with the thermal state and coherent state optical fields, respectively. From Figure 3, one can find that there are photon sidebands located at nω 0 (n = ±1, ±2, . . . . . .) resulting from the interaction of an electron on the QD with the optical field, which is similar to that of the QD system interacting with a phonon [42]. Once an electron jumps to a QD from an electrode lead, it can jump on the photon sideband outside the bias voltage window by absorbing a photon. Because of the Coulomb blockade and Pauli blockade, no other electrons can jump to the QD, and this leads to the dynamical channel blockade due to the small tunneling rate of the occupied channel. As the bias voltage increases, the dynamical channel blockade can be gradually removed if all the photon sidebands enter the bias voltage window. Thus, the Fano factor will walk into the regime dominated by the Pauli blockade and Coulomb blockade from the regime dominated by the dynamical channel blockade with an increasing bias voltage.
From Figure 2b, it can also be found that the Fano factor for the coherent state is less than that of the vacuum state for a small bias voltage. The reason is that F < coherent and F > coherent (Equation (10)) contain contributions from the off-diagonal terms of n|X(t)X † |m and n|X † X(t)|m (n = m), respectively, in the photon number representation. This will make the amplitude of the photon sideband of the coherent state fall below that of the zero-field vacuum contribution for certain ranges of parameters. Thus, the coherent state could suppress the tunneling contribution of the vacuum state for certain ranges of parameters [22]. Figure 2c exhibits the skewness-voltage characteristics for the vacuum state, thermal state, Fock state, coherent state, and coherent state with random phase. The S k for the coherent state decreases to about 0.7, then steps down to a constant value that is the same as that for the other quantum states with an increasing bias voltage V bias , but it remains a constant value for the other three quantum states as the bias voltage V bias increases. In order to explain the skewness-voltage characteristics, the curve of dark curves correspond to the same parameters as in Figure 2. Comparing the corresponding curves in Figures 2c and 4, there seems to be a correlation between skewness and ρ 0 − ρ 1 .   Figure 5a, it can be observed that the first step of the current is shifted due to the shift of the effective energy level of the QD. Figure 5a also exhibits that the current-voltage characteristics for the thermal state, Fock state, coherent state, and coherent state with random phase are almost the same as those in Figure 2a. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, there are no average currents when the bias voltage starts to increase. The shoulder height of the first and third steps of the current becomes smaller, and it becomes bigger for the second step for the coherent state.   Figure 2b. The Fano factors are at least several orders of magnitude higher than those in Figure 2b at lower bias voltage, and the Fano factor in the coherent state decreases to 0.5 with one step as bias voltage Vbias increases, and the step is not obvious for the thermal state, Fock state, and coherent state with random phase compared with that in Figure 2b. This super-Poissonian value of the Fano factor (Fa > 1) is due to avalanche-like transport [24]. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, the dynamical channel blockade arises at the beginning of the increase in the bias voltage, which results in avalanche-like transport for the four quantum states of optical field. The mechanism for avalanche-like transport in this paper is the dynamical channel blockade for a moderate electron-photon interaction regime (λ < 1), while the mechanism for avalanche-like transport in Reference [24] is one in which the transitions between   Figure 5 shows the FCS properties of QD charge transport versus the bias voltage V bias for different states of optical field for the same parameters as shown in Figure 2, except the effective energy level of the QD is shifted to 0.5. From Figure 5a, it can be observed that the first step of the current is shifted due to the shift of the effective energy level of the QD. Figure 5a also exhibits that the current-voltage characteristics for the thermal state, Fock state, coherent state, and coherent state with random phase are almost the same as those in Figure 2a. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, there are no average currents when the bias voltage starts to increase. The shoulder height of the first and third steps of the current becomes smaller, and it becomes bigger for the second step for the coherent state.   Figure 5a, it can be observed that the first step of the current is shifted due to the shift of the effective energy level of the QD. Figure 5a also exhibits that the current-voltage characteristics for the thermal state, Fock state, coherent state, and coherent state with random phase are almost the same as those in Figure 2a. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, there are no average currents when the bias voltage starts to increase. The shoulder height of the first and third steps of the current becomes smaller, and it becomes bigger for the second step for the coherent state.   Figure 2b. The Fano factors are at least several orders of magnitude higher than those in Figure 2b at lower bias voltage, and the Fano factor in the coherent state decreases to 0.5 with one step as bias voltage Vbias increases, and the step is not obvious for the thermal state, Fock state, and coherent state with random phase compared with that in Figure 2b. This super-Poissonian value of the Fano factor (Fa > 1) is due to avalanche-like transport [24]. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, the dynamical channel blockade arises at the beginning of the increase in the bias voltage, which results in avalanche-like transport for the four quantum states of optical field. The mechanism for avalanche-like transport in this paper is the dynamical channel blockade for a moderate electron-photon interaction regime (λ < 1), while the mechanism for avalanche-like transport in Reference [24] is one in which the transitions between  Figure 5b demonstrates that the Fano factor-voltage characteristics are different from that in Figure 2b. The Fano factors are at least several orders of magnitude higher than those in Figure 2b at lower bias voltage, and the Fano factor in the coherent state decreases to 0.5 with one step as bias voltage V bias increases, and the step is not obvious for the thermal state, Fock state, and coherent state with random phase compared with that in Figure 2b. This super-Poissonian value of the Fano factor (F a > 1) is due to avalanche-like transport [24]. Because the effective energy level of the QD is shifted away from the main channel and the photon sideband channel, the dynamical channel blockade arises at the beginning of the increase in the bias voltage, which results in avalanche-like transport for the four quantum states of optical field. The mechanism for avalanche-like transport in this paper is the dynamical channel blockade for a moderate electron-photon interaction regime (λ < 1), while the mechanism for avalanche-like transport in Reference [24] is one in which the transitions between low-lying phonon states are exponentially suppressed for strong electron-phonon coupling (λ 1). Moreover, the photon distribution described in this paper is in equilibrium, while the phonon distribution is in a nonequilibrium distribution in Reference [24]. Figure 5c demonstrates that the skewness-voltage characteristic for the coherent state is almost the same as that in Figure 2c, except that it has one more step before it decreases to a constant value. Compared with that in Figure 2c, at a lower bias voltage, the curve in Figure 5c has two steps before it decreases to a constant value for the thermal state, Fock state, and coherent state with random phase. Comparing the curves in Figure 5c with the dashed red curves in Figure 4, the correlation between skewness and ρ 0 − ρ 1 is apparent. Figure 6 shows the FCS properties of QD charge transport versus the bias voltage V bias for different states of optical field for the same parameters as shown in Figure 2 except for the lower temperature. Comparing Figure 6 with Figure 2, it can be observed that the curves in Figure 6 are almost the same as those in Figure 2, except that the curves change more sharply at the turning points in Figure 6 than those in Figure 2. This is because the Fermi distribution functions at lower temperature in Equations (7a)-(7d) change more steeply at the Fermi surface than that at higher temperature. (This is because the Fermi distribution functions in Equations (7a)-(7d) change more steeply at the Fermi surface when the temperature decreases). low-lying phonon states are exponentially suppressed for strong electron-phonon coupling (λ ≫ 1). Moreover, the photon distribution described in this paper is in equilibrium, while the phonon distribution is in a nonequilibrium distribution in Reference [24]. Figure 5c demonstrates that the skewness-voltage characteristic for the coherent state is almost the same as that in Figure 2c, except that it has one more step before it decreases to a constant value. Compared with that in Figure 2c, at a lower bias voltage, the curve in Figure 5c has two steps before it decreases to a constant value for the thermal state, Fock state, and coherent state with random phase. Comparing the curves in Figure 5c with the dashed red curves in Figure 4, the correlation between skewness and − is apparent.  Figure 6 shows the FCS properties of QD charge transport versus the bias voltage Vbias for different states of optical field for the same parameters as shown in Figure 2 except for the lower temperature. Comparing Figure 6 with Figure 2, it can be observed that the curves in Figure 6 are almost the same as those in Figure 2, except that the curves change more sharply at the turning points in Figure 6 than those in Figure 2. This is because the Fermi distribution functions at lower temperature in Equations (7a-7d) change more steeply at the Fermi surface than that at higher temperature. (This is because the Fermi distribution functions in Equations (7a-7d) change more steeply at the Fermi surface when the temperature decreases.)  Figure 7 shows the FCS properties of QD charge transport versus the bias voltage V bias for different states of optical field for the same parameters as shown in Figure 2 except for the stronger electron-photon coupling coefficient. Comparing Figure 7a to Figure 2a, we find that more steps appear before the average currents reach saturation values as the bias voltage V bias increases for all of the different states of the optical field, because there are more photon sidebands entering the bias voltage window when the electron-photon coupling becomes stronger. Figure 7b demonstrates that the evolution of the Fano factor-voltage characteristic for the coherent state shows completely different behavior from that observed in Figure 2b. However, the changes observed in Figure 7b for the other states are similar to those of the corresponding curves in Figure 2b. The Fano factors decrease quickly to 0.5 in one step with increasing bias voltage V bias . The plateau values of the steps are obviously higher than those in Figure 2b, respectively, and avalanche-like transport also exists at a lower bias voltage. This is because more photon sidebands appear with an increasing bias voltage, and this enforces the effect of the dynamical channel blockade.
temperature. Comparing Figure 6 with Figure 2, it can be observed that the curves in Figure 6 are almost the same as those in Figure 2, except that the curves change more sharply at the turning points in Figure 6 than those in Figure 2. This is because the Fermi distribution functions at lower temperature in Equations (7a-7d) change more steeply at the Fermi surface than that at higher temperature. (This is because the Fermi distribution functions in Equations (7a-7d) change more steeply at the Fermi surface when the temperature decreases.)   Figure 7c illustrates that the skewness-voltage characteristics are similar to that in Figure 2c for all of the different states of the optical field, except there are higher plateau values of the steps, and there is one more step for the optical field in the coherent state as bias voltage V bias increases. Calculations also confirm the correlation between skewness and ρ 0 − ρ 1 . Figure 8 shows the FCS properties of QD charge transport versus the bias voltage V bias for different states of optical field with the same parameters as Figure 2 except for the stronger light intensities. It illustrates that the FCS properties of QD charge transport behave similarly to those in Figure 2 besides the different shoulder heights of the steps. the different states of the optical field, because there are more photon sidebands entering the bias voltage window when the electron-photon coupling becomes stronger. Figure 7b demonstrates that the evolution of the Fano factor-voltage characteristic for the coherent state shows completely different behavior from that observed in Figure 2b. However, the changes observed in Figure 7b for the other states are similar to those of the corresponding curves in Figure 2b. The Fano factors decrease quickly to 0.5 in one step with increasing bias voltage Vbias. The plateau values of the steps are obviously higher than those in Figure 2b, respectively, and avalanchelike transport also exists at a lower bias voltage. This is because more photon sidebands appear with an increasing bias voltage, and this enforces the effect of the dynamical channel blockade. Figure 7c illustrates that the skewness-voltage characteristics are similar to that in Figure 2c for all of the different states of the optical field, except there are higher plateau values of the steps, and there is one more step for the optical field in the coherent state as bias voltage Vbias increases. Calculations also confirm the correlation between skewness and − .  Figure 8 shows the FCS properties of QD charge transport versus the bias voltage Vbias for different states of optical field with the same parameters as Figure 2 except for the stronger light intensities. It illustrates that the FCS properties of QD charge transport behave similarly to those in Figure 2 besides the different shoulder heights of the steps.  The skewness curves coincide for the thermal state, Fock state, and coherent state with random phase. Figure 9 exhibits I , F a , and S k versus the bias voltage V bias for the coherent state with different phases. It shows that the phase of the coherent state can evidently change the FCS property: In particular, negative differential conductance appears at a bias voltage of nearly 4 and 6 in Figure 9a. As for the other states, the thermal state and coherent state with random phase are mixed states, and the phase of the Fock state is completely uncertain. Here, only the coherent state optical field has the advantage of tuning the FCS properties by phase. The interaction between the QD and coherent state optical field can lead to the appearance of obvious negative differential conductance in a symmetric QD system, while the pronounced negative differential conductance, which usually appears in a similar QD system assisted by phonons with asymmetric parameters and negative differential conductance, is ordinarily much less pronounced in symmetric junctions [9,40]. skewness curves coincide for the thermal state, Fock state, and coherent state with random phase. Figure 8 shows the FCS properties of QD charge transport versus the bias voltage Vbias for different states of optical field with the same parameters as Figure 2 except for the stronger light intensities. It illustrates that the FCS properties of QD charge transport behave similarly to those in Figure 2 besides the different shoulder heights of the steps.  Figure 2. The curve of π/4 is the same as that of 3π/4. Figure 9 exhibits 〈 〉, Fa, and Sk versus the bias voltage Vbias for the coherent state with different phases. It shows that the phase of the coherent state can evidently change the FCS property: In particular, negative differential conductance appears at a bias voltage of nearly 4 and 6 in Figure 9a. As for the other states, the thermal state and coherent state with random phase are mixed states, and the phase of the Fock state is completely uncertain. Here, only the coherent state optical field has the , and S k (c) versus the bias voltage V bias for the coherent state with different phases for the same parameters as Figure 2. The curve of π/4 is the same as that of 3π/4. Figure 10 exhibits I , F a , and S k versus the phase of the coherent state with fixed bias voltage, V bias = 3. From Figure 10, we can find that I , F a , and S k can be changed by a period of π for the phase of the coherent state. Although the definite phase of the optical field is unmeasurable, this result implies that the FCS properties can be periodically changed by tuning the relative phase of the laser beam. advantage of tuning the FCS properties by phase. The interaction between the QD and coherent state optical field can lead to the appearance of obvious negative differential conductance in a symmetric QD system, while the pronounced negative differential conductance, which usually appears in a similar QD system assisted by phonons with asymmetric parameters and negative differential conductance, is ordinarily much less pronounced in symmetric junctions [9,40].   Figure 10, we can find that 〈 〉, Fa, and Sk can be changed by a period of π for the phase of the coherent state. Although the definite phase of the optical field is unmeasurable, this result implies that the FCS properties can be periodically changed by tuning the relative phase of the laser beam.
Finally, Figure 11 demonstrates the FCS properties of QD charge transport versus the electronphoton coupling coefficient under infinitely high bias voltage. Figure 11 shows exponentially suppressed average currents for the different states of optical field due to the very strong electronphoton interaction. The Fano factor and skewness always stay within the constant value dominated by the Coulomb repulsion and the Pauli exclusion principle, because all the photon sidebands are in the window of bias voltage.  Finally, Figure 11 demonstrates the FCS properties of QD charge transport versus the electron-photon coupling coefficient under infinitely high bias voltage. Figure 11 shows exponentially suppressed average currents for the different states of optical field due to the very strong electron-photon interaction. The Fano factor and skewness always stay within the constant value dominated by the Coulomb repulsion and the Pauli exclusion principle, because all the photon sidebands are in the window of bias voltage. advantage of tuning the FCS properties by phase. The interaction between the QD and coherent state optical field can lead to the appearance of obvious negative differential conductance in a symmetric QD system, while the pronounced negative differential conductance, which usually appears in a similar QD system assisted by phonons with asymmetric parameters and negative differential conductance, is ordinarily much less pronounced in symmetric junctions [9,40].   Figure 10, we can find that 〈 〉, Fa, and Sk can be changed by a period of π for the phase of the coherent state. Although the definite phase of the optical field is unmeasurable, this result implies that the FCS properties can be periodically changed by tuning the relative phase of the laser beam.
Finally, Figure 11 demonstrates the FCS properties of QD charge transport versus the electronphoton coupling coefficient under infinitely high bias voltage. Figure 11 shows exponentially suppressed average currents for the different states of optical field due to the very strong electronphoton interaction. The Fano factor and skewness always stay within the constant value dominated by the Coulomb repulsion and the Pauli exclusion principle, because all the photon sidebands are in the window of bias voltage.

Conclusion
In summary, we employed the particle-number-resolved master equation and took into account the influences of optical fields in the thermal state, Fock state, coherent state, and coherent state with

Conclusion
In summary, we employed the particle-number-resolved master equation and took into account the influences of optical fields in the thermal state, Fock state, coherent state, and coherent state with random phase. Then, we investigated the FCS properties of single-electron tunneling transport on a QD in a moderate electron-photon interaction regime by single-particle approximation. The results show that the interaction with the optical fields in different quantum states can lead to critical effects on the average tunneling current, Fano factor, and skewness of single-electron tunneling transport on the QD.
Firstly, under the same light intensity, the average tunneling currents have different shoulder heights for all of the studied quantum states of optical fields when the bias voltage V bias is lower, while they have the same saturation current as the bias voltage V bias becomes higher. The vacuum state optical field can unexpectedly induce a staircase-shaped current. For the coherent state, the characteristics of the Fano factor and skewness are different from those of all other states.
Secondly, the FCS properties are sensitive to the energy level of the QD, the electron-photon coupling coefficient, and the phase of the coherent state for single-electron tunneling transport on a QD. With regard to avalanche-like transport at a lower bias voltage, the mechanism in a moderate electron-photon interaction regime is the dynamical channel blockade, which is different from that for a strong electron-photon interaction regime. However, the Fano factor voltage characteristics are dominated by the Pauli blockade and Coulomb blockade at higher bias voltages.
Thirdly, pronounced negative differential conductance results from tuning the phase of the coherent state optical field in a symmetric QD system, while the pronounced negative differential conductance usually appears in a similar QD system assisted by phonons with asymmetric parameters. The FCS properties vary with the phase of the coherent state in a π-period.
Finally, under infinitely high bias voltages, the evolution behaviors of the average current diverge at higher electron-photon coupling coefficients for the different states of optical field.
On the whole, for single-electron tunneling transport on a QD, the FCS properties can be strictly regulated by the optical field under different quantum states, especially by the phase of the coherent state and QD energy level adjusted by gate voltage. Thus, the investigation is significant for controlling the transport of a QD electron for nanostructure devices.
Author Contributions: All Authors prepared the manuscript and contributed to the concepts and theoretical results exposed in this article. Numerical calculations were run by W.L. All authors have read and approved the final manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.