Open AccessArticle

Quantum Photovoltaic Cells Driven by Photon Pulses

Sangchul Oh

^1,*

Jung Jun Park

² and

Hyunchul Nha

Qatar Environment and Energy Research Institute, Hamad Bin Khalifa University, Qatar Foundation, P.O. Box 5825 Doha, Qatar

Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 02455, Korea

Department of Physics, Texas A&M University at Qatar, Education City, P.O. Box 23874 Doha, Qatar

Author to whom correspondence should be addressed.

Entropy 2020, 22(6), 693; https://doi.org/10.3390/e22060693

Submission received: 23 May 2020 / Revised: 12 June 2020 / Accepted: 16 June 2020 / Published: 20 June 2020

Download

Browse Figures

Review Reports Versions Notes

Abstract

We investigate the quantum thermodynamics of two quantum systems, a two-level system and a four-level quantum photocell, each driven by photon pulses as a quantum heat engine. We set these systems to be in thermal contact only with a cold reservoir while the heat (energy) source, conventionally given from a hot thermal reservoir, is supplied by a sequence of photon pulses. The dynamics of each system is governed by a coherent interaction due to photon pulses in terms of the Jaynes-Cummings Hamiltonian together with the system-bath interaction described by the Lindblad master equation. We calculate the thermodynamic quantities for the two-level system and the quantum photocell including the change in system energy, the power delivered by photon pulses, the power output to an external load, the heat dissipated to a cold bath, and the entropy production. We thereby demonstrate how a quantum photocell in the cold bath can operate as a continuum quantum heat engine with a sequence of photon pulses continuously applied. We specifically introduce the power efficiency of the quantum photocell in terms of the ratio of output power delivered to an external load with current and voltage to the input power delivered by the photon pulse. Our study indicates a possibility that a quantum system driven by external fields can act as an efficient quantum heat engine under non-equilibrium thermodynamics.

Keywords:

open quantum system; photovoltaic cell; quantum heat engines; quantum thermodynamics; master equations

1. Introduction

Thermodynamics deals with the evolution of systems, usually in contact with reservoirs, describing the dynamics under universal laws independent of microscopic details. Among its four laws, the second law dictates the total entropy of a closed system can never decrease over time and that the closed system spontaneously evolves toward the state with maximum entropy. One of the possible statements about the second law of thermodynamics is to set the upper bound on the efficiency of heat engines. Heat engines convert heat energy, which typically flows from a hot source to a cold sink, to mechanical energy or chemical energy. The efficiency of energy conversion is defined by the ratio of the work output to the amount of heat energy input. The ultimate efficiency of the heat engine is known in equilibrium thermodynamics to be determined only by temperatures of hot and cold heat baths,

T_{h}

and

T_{c}

, respectively, that is,

η = 1 - T_{c} / T_{h}

, the so-called the Carnot limit.

Photovoltaic cells (or solar cells) and photosynthesis, just like classical heat engines, convert photon energy from the sun into electric energy and chemical energy, respectively. The upper limit of the efficiency of p-n junction solar cells with an energy bandgap is known as the Shockley-Queisser limit [1]. The key assumptions in deriving the Shockley-Queisser limit are (i) photons with energies less than the bandgap are not utilized, (ii) a photon with energy greater than the bandgap produces only one electron-hole pair, and (iii) only the radiative recombination of electron-hole pairs is considered. While the non-radiative loss may be minimized by the manufacturing technology, the radiative recombination is the intrinsic energy loss governed by the law of physics. Assuming the sun and the solar cell are described as black-bodies with temperature

T_{s} = 6000 K

and

T_{c} = 300 K

, respectively, the maximum efficiency is about 30% for a solar cell with a bandgap of

1.137 eV

[1]. Shockley and Queisser [1] also showed that the maximum efficiency of a single p-n junction solar cell would be approximately 44% around 1.137 eV if there is no radiative recombination loss.

Recently, many theoretical studies have suggested that noise-induced quantum coherence [2,3], Fano-induced coherence [4] or delocalized quantum states of interacting dipoles [5,6,7,8] can reduce the radiative recombination loss of a solar cell, thus enhancing the efficiency of solar cells. The same idea is applied to the photosynthetic complex [9,10,11,12]. Most of these studies employ the donor-acceptor quantum photocell model where the donor is in thermal equilibrium with a hot bath, that is, the sun at

5800 K

and the acceptor is at room temperature. The photocell operating while continuously contacting both heat reservoirs is called a continuous quantum heat engine [13,14]. A typical example of the continuous quantum heat engine is Scovil and Schulz-Dubois’ three-level masers whose efficiency achieves the Carnot efficiency [15]. By solving the master equation for the quantum photocell, it was shown that the noise-induced quantum coherence or the dark state of the donor may enhance the power output. However, it was not clear whether the efficiency of the photocell could be enhanced by the quantum effect. Some works assumed that the mean photon number of the hot thermal bath could be

\bar{n} = 60000

, but the mean photon number of the sun at energy

1.8 eV

as a black body at

5800 K

is only about

\bar{n} = 0.037

[16,17]. To address this issue, Reference [16] introduced the pumping term and showed that the dark state could enhance the power output but not its efficiency.

In this paper, we explore another form of a quantum heat engine. We consider two quantum systems, a two-level system and a donor-acceptor quantum photocell, and investigate their quantum dynamics under the coherent driving and the system-bath interaction. Each quantum system is in thermal contact only with a cold reservoir but not with a hot reservoir. Instead, they are driven by a sequence of photon pulses that supplies input energy to the systems, which is conventionally done by a hot reservoir. The photon pulses represent the stream of energy source to the system and may thus remove the unrealistic assumption of the high mean photon number of the sun by the previous works. We solve the time-dependent Markov Lindblad master equation and investigate the thermodynamic quantities such as the change in energy, the heat dissipation to the cold bath, the power delivered by the photon pulse and the entropy generation. Specifically, we introduce the power efficiency of the quantum photocell in terms of the ratio of output power delivered to an external load to input power delivered by the photon pulses.

This paper is organized as follows. In Section 2, we review briefly the quantum dynamics and quantum thermodynamics of an open quantum system based on the master equation approach. In Section 3, we examine a two-level system in a cold bath driven by photon pulses and investigate how the energy, heat current, and entropy change. In Section 4 a quantum photovoltaic cell with donor and acceptor driven by repeated photon pulses is considered. We calculate the quantum thermodynamic quantities and the power output by the sequence of the photon pulses together with engine efficiency. Finally, in Section 5 we summarize our results with some discussion.

2. Quantum Thermodynamics of Open Quantum Systems

We start with a brief review of quantum thermodynamics of an open quantum system that exchanges energy and entropy with its environment. The equations presented in this section will be applied to those examples in the next two Section 3 and Section 4. As usual, we assume the Born-Markov approximations: a weak interaction between an open quantum system and an environment, and the extremely short correlation time of the environment, that is, no memory effect. The density operator

ρ (t)

of the quantum system with a slowly varying time-dependent Hamiltonian obeys the Lindblad-Gorini-Kossakowski-Sudarshan (LGKS) master equation [18,19,20,21,22]

\frac{d}{d t} ρ (t) = L [ρ (t)] = - \frac{i}{ℏ} [H_{S} (t), ρ (t)] + D [ρ (t)],

(1)

where

H_{S} (t) = H_{0} + H_{1} (t)

is the Hamiltonian of the system. Here

H_{0}

represents a time-independent unperturbed Hamiltonian and

H_{1} (t)

an external time-dependent perturbation. The decoherence and the dissipation of the open quantum system due to an environmental interaction are described by the non-unitary operator

D [ρ] = \sum_{k} (2 L_{k} ρ L_{k}^{†} - L_{k}^{†} L ρ - ρ L_{k}^{†} L_{k}),

(2)

where

L_{k}

are the Lindblad operators determined according to the type of interaction.

From the solution

ρ (t)

of Equation (1), one can calculate the quantum thermodynamic quantities. The first law of classical thermodynamics states the energy conservation,

d E = δ Q + δ W

. The time-dependent internal energy of the system is given by

E (t) = tr \{ρ (t) H_{S} (t)\}

. Its derivative with respect to time gives rise to the first law of quantum thermodynamics [23,24,25]

\frac{d}{d t} E (t) = \dot{Q} (t) + \dot{W} (t) = J (t) + P (t) .

(3)

Here

J (t)

is the heat current from the environment into the system

J (t) \equiv \dot{Q} (t) = tr (\frac{d ρ (t)}{d t} H_{S} (t)),

(4)

and

P (t)

is the power delivered to system by external forces,

P (t) \equiv \dot{W} (t) = tr (ρ (t) \frac{d H_{S}}{d t}) .

(5)

Since

H_{0}

is the time-independent Hamiltonian, the power can be written as

P (t) = tr (ρ (t) \frac{d H_{1} (t)}{d t})

. The change in energy of the system for finite time can be obtained by integrating the heat current and the power as

Δ E (t) = Q (t) + W (t) = \int_{0}^{t} J (s) d s + \int_{0}^{t} P (s) d s .

(6)

The second law of thermodynamics describes the irreversibility of dynamics, where the entropy plays a key role. The von Neumann entropy

S (t)

of the system in the state

ρ (t)

is given by

S (t) = - tr \{ρ (t) log ρ (t)\} .

(7)

The thermodynamic entropy

S

is written as

S = k_{B} S (t)

with the Boltzmann constant

k_{B}

. The net change in the entropy

d S_{net}

(entropy production) of the whole system+reservoir can be written in terms of the entropy change of the system,

d S

, and the entropy flow due to heat from an environment to a system,

d S_{e}

, as

d S_{net} = d S - d S_{e} .

(8)

The change in the entropy of the system,

d S

, over time is written as

\begin{matrix} \frac{d S}{d t} & = - tr \{\dot{ρ} (t) log ρ (t)\} \end{matrix}

(9a)

\begin{matrix} = - tr \{L [ρ (t)] log ρ (t)\}, \end{matrix}

(9b)

where

tr \{\dot{ρ} (t)\} = 0

and the quantum Markov master Equation (1) are used. The entropy flow

d S_{e}

per unit time from the environment into the system is written as

\begin{matrix} \frac{d S_{e}}{d t} & \equiv J_{S} = β \dot{Q} (t) \end{matrix}

(10a)

\begin{matrix} = β tr \{\dot{ρ} H_{S} (t)\} = β tr \{L [ρ (t)] H_{S} (t)\}, \end{matrix}

(10b)

where

β = 1 / k_{B} T

and T is the temperature of the environment. The net entropy production rate

σ (t)

of the system is given by

σ (t) = \frac{d S_{net}}{d t} = \dot{S} (t) - β \dot{Q} (t) \geq 0,

(11)

where

σ (t) \geq 0

comes from the Spohn inequality [23,24,25,26]. Equation (11) may be written as

σ \equiv - \frac{d}{d t} S (ρ (t) ‖ ρ_{ss}) \geq 0,

(12)

where

S (ρ (t) ∥ ρ_{ss}) \equiv tr [ρ (t) (log ρ (t) - log ρ_{ss})]

is the relative entropy of

ρ (t)

with respect to the stationary state

ρ_{ss}

, for example, the canonical state of the system,

ρ_{ss} = ρ_{β} = e^{- β H_{s} (t)} / Z

. This is called the second law of nonequilibrium quantum thermodynamics in the weak coupling limit.

As an application of quantum thermodynamics of open quantum systems presented in Section 2, we first consider a two-level quantum system which is in contact with a cold bath at temperature

T_{c}

and driven by repeated photon pulses, as depicted in Figure 1. The hot thermal bath supplying energy does not have direct contact with the quantum system. Its role is here replaced by a sequence of photon pulses to the two-level system. We examine how the two-level system absorbs and dissipates energy and generates entropy during this process in order to gain insight into the nonequilibrium dynamics due to photon pulses.

3. A Two-Level System Driven by Photon Pulses

The unperturbed Hamiltonian of the two-level system with energy levels

E_{0}

and

E_{1}

may be written as

\begin{matrix} H_{0} = - \frac{ℏ ω_{0}}{2} σ_{z}, \end{matrix}

(13)

where

ω_{0} = (E_{1} - E_{0}) / ℏ

and

σ_{z} = | 0 〉 〈 0 | - | 1 〉 〈 1 |

. The interaction between a two-level system and incoming photon pulses is described by the Jaynes-Cummings Hamiltonian [27]

H_{1} (t) = i ℏ [g^{*} (t) \sqrt{γ} σ_{-} - g (t) \sqrt{γ} σ_{+}],

(14)

where

σ_{+} = | 1 〉 〈 0 |

and

σ_{-} = | 0 〉 〈 1 |

are the raising and the lowering operators, respectively. The Jaynes-Cummings Hamiltonian is derived based on the dipole and rotating wave approximations [28]. We set

E_{1} - E_{0} = 1 eV

. Here

γ

is the Weisskopf-Wigner spontaneous decay rate

\begin{matrix} γ = \frac{1}{4 π ϵ_{0}} \frac{4 ω_{0}^{3} d_{01}^{2}}{3 ℏ c^{3}}, \end{matrix}

(15)

where

d_{01}

is the transition dipole moment between the two states

| 0 〉

and

| 1 〉

. A typical value of

γ

for an atom or a quantum dot for visible light emission is in the order of nanoseconds corresponding to

μ eV

, while the energy at visible frequencies is about eV, that is, femtoseconds. As a numerical calculation becomes demanding with a big difference between these time-scales, we use for our study the values of parameters as listed in Table 1. We consider the photon pulses given at peak times

t_{i}

g (t) = α \sum_{i} ξ (t; t_{i})

with coherent states having average photon number

〈 n 〉 = {| α |}^{2}

and a Gaussian pulse shape [29,30,31]

\begin{matrix} ξ (t; t_{i}) \equiv {(\frac{Ω^{2}}{2 π})}^{1 / 4} exp [- \frac{Ω^{2} {(t - t_{i})}^{2}}{4} - i ω_{0} t] . \end{matrix}

(16)

Here

1 / Ω

is the pulse bandwidth.

Under the Born-Markov approximation, the interaction of a two-level system and the thermal photon bath is recast to the dissipative operator

D

acting on the density matrix of the system [21,32]

\begin{matrix} D_{C} [ρ] & = \frac{γ}{2} ({\bar{n}}_{c} + 1) (2 σ_{-} ρ σ_{+} - σ_{+} σ_{-} ρ + ρ σ_{+} σ_{-}) \\ + \frac{γ}{2} {\bar{n}}_{c} (2 σ_{+} ρ σ_{-} - σ_{-} σ_{+} ρ + ρ σ_{-} σ_{+}) . \end{matrix}

(17)

Here

{\bar{n}}_{c}

is the mean photon number of the cold bath at the frequency

ω_{0}

in thermal equilibrium of temperature

T_{c}

\begin{matrix} {\bar{n}}_{c} = \frac{1}{e^{ℏ ω_{0} / k_{B} T_{c}} - 1} . \end{matrix}

(18)

As noted in References [32,33], at optical frequencies and room temperature, the mean photon number

\bar{n}

is very small and negligible while it has a finite value at microwave frequencies and the room temperature. For example, with

ℏ ω_{0} = eV

and

T_{c} = 300 K

one obtains

\bar{n} \approx 6.5 \times 10^{- 31}

. At the optical frequencies,

ℏ ω_{0} = 1.8 eV

and the temperature of the sun as a black body,

T_{s} = 5800 K

, the mean photon number is

\bar{n} \approx 0.0317

[16,33]. With this in mind, Equation (17) reduces to

\begin{matrix} D_{C} [ρ] \approx \frac{γ}{2} (2 σ_{-} ρ σ_{+} - σ_{+} σ_{-} ρ + ρ σ_{+} σ_{-}) . \end{matrix}

(19)

With Equations (13), (17), and (14), the Lindblad equation for the two-level system, in contact with the cold thermal bath and driven by a Gaussian photon pulse, is given by

\frac{d}{d t} ρ (t) = - \frac{i}{ℏ} [H_{0} + H_{1} (t), ρ (t)] + D_{C} [ρ (t)] .

(20)

The quantum dynamics and the quantum thermodynamics of the two-level system are investigated by solving Equation (20) numerically using the Runge-Kutta method [34]. The parameters used in the numerical simulation are shown in Table 1.

Figure 2 and Figure 3 describe the thermodynamic quantities of the system when a sequence of Gaussian photon pulses are applied at a regular interval (

g (t)

: green curves in (a)) with

〈 n 〉 = 1

and

〈 n 〉 = 10

, respectively. Figure 2a shows the time-evolution of the density matrix elements and the sequence of Gaussian photon pulses, which is first applied around the peak time

ω_{0} t / 2 π = 50

with

〈 n 〉 = 1

. The initial state of the two level system is assumed to be in a superposed state

| ψ (0) 〉 = \frac{1}{\sqrt{2}} (| 0 〉 + | 1 〉)

. As shown in Figure 2a, the superposed state decays to the ground state, that is, the system becomes in thermal equilibrium with the cold thermal bath before the photon pulse is applied. When the Gaussian photon pulse is first applied around

ω_{0} t / 2 π = 50

, the system gets excited and then becomes decayed into the ground state after the pulse is gone. This process is repeated according to each Gaussian pulse.

With regard to the first law of quantum thermodynamics, Figure 2b plots the rate of energy change, the power, and the heat current. The heat current

J (t)

from the environment to the system is always negative. This means the excited state of the two-level system releases its energy to the cold bath. In contrast, the power

P (t)

and

\dot{E} (t)

oscillate out of phase while the photon pulse is applied. Figure 2c shows the entropy

S (t)

of the two-level system together with its energy

E (t)

, the work done

W (t)

on the system and the heat transfer

Q (t)

. Figure 2d shows the entropy production

σ (t) = \dot{S} (t) - β \dot{Q}

as a function of time in relation to the second law of thermodynamics. The entropy production

σ (t)

is always positive confirming the second law.

In Figure 3, we see more oscillatory behaviors in the quantities due to a stronger photon pulse with

〈 n 〉 = 10

than those with

〈 n 〉 = 1

in Figure 2. Nevertheless, the overall trend is similar to that explained above for Figure 2.

We now examine how the quantum thermodynamic quantities depend on the temporal shapes of a Gaussian pulse sequence. In Figure 4 and Figure 5, we plot the same quantities as those in Figure 2 and Figure 3, but compare two cases, that is, regularly spaced (Figure 4) and irregularly spaced (Figure 5) sequence of Gaussian pulses with the same mean number

〈 n 〉 = {| α |}^{2} = 1

. As described by the curve

g (t)

, the peak times

t_{i}

’s in Equation (16) are regularly (not regularly) spaced in the left (right) panel. In both cases, the overall trend of the thermodynamic quantities are similar to that explained for Figure 2 while the actual response of the system does depend on the temporal shape of the pulse sequence. Remarkably, we see that the output power

P (t)

(blue curves in Figure 4b and Figure 5b) and the accumulated work

W (t)

(blue curves in Figure 4c and Figure 5c) depend on the temporal shape of the incoming pulses even with the same

{| α |}^{2} = 1

, which can have implications for practical photocell operation. In particular, we find that the case of regular sequencing of pulses yields a higher value of work.

4. Quantum Photocell Driven by Photon Pulses

For a quantum heat engine, let us now consider a quantum photovoltaics cell driven by photon pulses, as shown in Figure 6. The quantum photocell we consider is a 4-level quantum system composed of a donor and an acceptor. In 1959, Scovil and Du-Bois [15] consider the 3-level system as the simplest quantum heat engine where one part of the 3-level system is in thermal equilibrium with a hot bath and the other part with a cold bath. Many previous studies [2,3,4,5,6,7,8,9,10] took a similar assumption that the donor of the quantum photocell is in contact with the hot bath, that is, the sun, and the acceptor is in thermal contact with the cold bath. In contrast to the previous works, we assume that the quantum photocell is in thermal contact only with the cold bath. In our previous work [16], the pumping term was introduced in the Lindblad master equation to describe the energy flow from the hot bath. Here the input energy is supplied by the sequence of incoming photon pulses.

The cyclic operation of the quantum photocell can be performed with the sequence as follows: (i) The donor absorbs incoming photons and the electron becomes excited with the transition from the ground state

| 0 〉

to the excited state

| 1 〉

. (ii) The phonon vibration makes the excited electron at the donor transfer to the acceptor state

| 2 〉

. (iii) The acceptor is coupled to an external load and the current flow (electric work) is represented by the transition decay from the state

| 2 〉

to the state

| 3 〉

. (iv) The electron in the state

| 3 〉

of the acceptor returns to the ground state

| 0 〉

of the donor by a vibrational or non-radiative decay.

The unperturbed Hamiltonian of the quantum photocell with 4-levels is written as

H_{0} = - E_{0} | 0 〉 〈 0 | - E_{1} | 1 〉 〈 1 | - E_{2} | 2 〉 〈 2 | - E_{3} | 3 〉 〈 3 | .

(21)

Similar to Equation (14), the interaction of the donor of the photocell with the incoming photon pulses is again described by the Jaynes-Cummings Hamiltonian

H_{1} (t) = i ℏ [g^{*} (t) \sqrt{γ} σ_{-} - g (t) \sqrt{γ} σ_{+}],

(22)

where

σ_{+} = | 1 〉 〈 0 |

and

σ_{+} = | 0 〉 〈 1 |

. Same as Equation (17), the interaction of the donor of the quantum photocell with the cold thermal bath is represented by the Lindblad operator,

\begin{matrix} D_{c} [ρ] & = \frac{γ_{01}}{2} ({\bar{n}}_{c} + 1) (2 σ_{-} ρ σ_{+} - σ_{+} σ_{-} ρ + ρ σ_{+} σ_{-}) \\ + \frac{γ_{01}}{2} {\bar{n}}_{c} (2 σ_{+} ρ σ_{-} - σ_{-} σ_{+} ρ + ρ σ_{-} σ_{+}) . \end{matrix}

(23)

The electron transfer between the states

| 1 〉

and

| 2 〉

and that between the state

| 3 〉

and

| 0 〉

are described by the Lindblad operator

D_{ph} [ρ] = D_{ph}^{(1, 2)} [ρ] + D_{ph}^{(3, 0)} [ρ]

, where

\begin{matrix} D_{ph}^{(i, j)} [ρ] & = \frac{γ_{i j}}{2} ({\bar{n}}_{ph} + 1) (2 L_{i j} ρ L_{i j}^{†} - L_{i j}^{†} L_{i j} ρ + ρ L_{i j}^{†} L_{i j}) \\ + \frac{γ_{i j}}{2} {\bar{n}}_{ph} (2 L_{i j}^{†} ρ L_{i j} - L_{i j} L_{i j}^{†} ρ + ρ L_{i j} L_{i j}^{†}) . \end{matrix}

(24)

Here

γ_{12}

and

γ_{30}

represents the transition rate between

| 1 〉

and

| 2 〉

and between

| 3 〉

and

| 0 〉

, respectively.

L_{i j} = | i 〉 〈 j |

and

L_{i j}^{†} = | j 〉 〈 i |

are the lowering and raising operators, respectively.

{\bar{n}}_{ph}

is the phonon occupation number at

ℏ ω = E_{1} - E_{2} = E_{3} - E_{0}

and

T_{c} = 300 K

. The work done by the quantum photocell to the external load is described by the ohmic dissipation

D_{ohm} [ρ] = \frac{Γ}{2} (2 L_{3} ρ L_{3}^{†} - L_{3}^{†} L_{3} ρ + ρ L_{3}^{†} L_{3}),

(25)

where

L_{3} = | 3 〉 〈 2 |

. Here

Γ

represent the conductance of the external load and may be changed from zero corresponding to the open circuit and to a big number representing the short-circuit of the quantum photocell. With Equations (21)–(25), the LGKS equation for the quantum photocell is written as

\begin{matrix} \frac{d}{d t} ρ (t) = - \frac{i}{ℏ} [H_{0} + H_{1} (t), ρ (t)] + D_{c} [ρ] + D_{ph}^{(1, 2)} [ρ] + D_{ph}^{(3, 0)} [ρ] + D_{ohm} [ρ] . \end{matrix}

(26)

Since the quantum photocell has no direct interaction Hamiltonian between the donor and the acceptor, the Hamiltonian

H_{S} (t) = H_{0} + H_{1} (t)

of the quantum photocell can be written as the sum of the time-dependent donor Hamiltonian

H_{D}

and the time-independent acceptor Hamiltonian

H_{A}

\begin{matrix} H_{S} (t) = H_{D} (t) + H_{A}, \end{matrix}

(27)

where

H_{D} (t) \equiv - E_{0} | 0 〉 〈 0 | - E_{1} | 1 〉 〈 1 | + H_{1} (t)

and

H_{A} \equiv - E_{2} | 2 〉 〈 2 | - E_{3} | 3 〉 〈 3 |

. This partition makes it possible to express some quantum thermodynamic quantities as the sum of the donor and acceptor parts. The energy of the quantum photocell is given by the sum of the energies of the donor and acceptor

\begin{matrix} E (t) = tr \{ρ (t) H_{S} (t)\} = E_{D} (t) + E_{A} (t), \end{matrix}

(28)

where the donor energy

E_{D} (t)

and the acceptor energy

E_{A} (t)

are given by

\begin{matrix} E_{D} (t) & = {tr}_{D} {ρ_{D} (t) H_{D} (t)}, \end{matrix}

(29a)

\begin{matrix} E_{A} (t) & = {tr}_{A} {ρ_{A} (t) H_{A}}, \end{matrix}

(29b)

respectively. Here

ρ_{A} = {tr}_{D} {ρ}

and

ρ_{D} = {tr}_{A} {ρ}

are the density operators of the donor and acceptor, respectively. Since the photon pulse delivers the power only to the donor, the power

P (t)

is given by the power of the donor

\begin{matrix} P (t) = tr {ρ (t) {\dot{H}}_{S} (t)} = {tr}_{D} {ρ_{D} {\dot{H}}_{1} (t)} \equiv P_{D} (t) . \end{matrix}

(30)

The heat dissipation occurs at the donor and acceptor. Thus the heat current

J (t)

is written as the sum of the two parts

\begin{matrix} J (t) & = {tr}_{D} {{\dot{ρ}}_{D} (t) H_{D} (t)} + {tr}_{A} {{\dot{ρ}}_{A} (t) H_{A}} \\ = J_{D} (t) + J_{A} (t) . \end{matrix}

(31)

Finally, the entropy of the quantum photocell

S (t) = tr {ρ log ρ}

can be written as the sum of the entropies of the donor and acceptor,

S_{D} (t) = - {tr}_{D} {ρ_{D} log ρ_{D}}

and

S_{A} (t) = - {tr}_{A} {ρ_{A} log ρ_{A}}

, too. This is because there is no coherent interaction between the donor and the acceptor so the whole density operator of the system has a structure

ρ = ρ_{D} \oplus ρ_{A}

We calculate the current through the external load as

I = e Γ \cdot ρ_{22},

(32)

and the voltage across the external load

e V = E_{2} - E_{3} + k_{B} T_{C} log (\frac{ρ_{22}}{ρ_{33}}) .

(33)

The latter comes from the relation

ρ_{22} / ρ_{33} = exp (- (E_{2} - E_{3} - e V) / k_{B} T_{C})

. The electric power delivered to the external load by the photocell is written as

P_{out} = I (t) \cdot V (t)

which depends on the external conductance

Γ

. Now that we have the power delivered by the photon pulse, Equation (30) and the electric power output

P_{out} (t)

, we can define the power efficiency of the quantum photocell as

We solve numerically the LGKS Equation (26) using the Runge-Kutta method for different sequences of photon pulses to obtain the quantum thermodynamic quantities. Figure 7a plots the population of each level of the quantum photocell as a function of dimensionless time when the sequences of the Gaussian pulses are applied one immediately after another almost in a continuum limit. Figure 7b shows the change in energies of the donor and acceptor,

{\dot{E}}_{D} (t)

and

{\dot{E}}_{A} (t)

, the power

P_{D} (t)

delivered to the donor by the photon, the power output

P_{out}

, the heat dissipation at the donor and the acceptor,

J_{D} (t)

and

J_{A} (t)

. Note that

P_{D} (t)

becomes equal to

J_{D} (t)

in the steady state. The heat dissipation

J_{D} (t)

of the donor is mainly associated with the transfer of electrons from the donor to the acceptor. The heat flow from the donor to the environment is relatively small because the decay rate

γ_{01}

to the environment is much smaller than the electron transfer rate

γ_{12}

from the donor to the acceptor. We note that only the fraction of the transferred heat energy is converted to the power output

P_{out}

J_{D} (t)

mainly plays the role of populating the acceptor level 2, which generates the acceptor current for power output as indicated by Equation (32). It would be interesting to see whether the coherent transfer between the donor and the acceptor may give rise to a different result.

Figure 7c shows the total entropy

S (t)

of the system and the entropies of the donor and the acceptor,

S_{D} (t)

and

S_{A} (t)

. Our numerical calculation confirms that the total entropy is the sum of those of the donor and the acceptor,

S (t) = S_{D} (t) + S_{A} (t)

, as explained before.

Figure 7d depicts the current

I (t)

, the voltage

V (t)

, and the power efficiency

η

defined by the ratio of the electric power output

P_{out} (t)

and the power delivered by the photon

P_{D} (t)

. From the figures, we see that these quantities initially show an oscillatory behavior then become saturated in the long time limit. In particular, the asymptotic value of power efficiency is as high as

η_{p} \sim 0.36

, which can also be interpreted as the work efficiency, that is, work output divided by energy input, when

η_{p}

is constant.

There are different types of quantum heat engines like continuous engine, two-stroke engine and four-stroke engine. Many other studies considered the quantum photocell as the continuous heat engine where the donor is in thermal contact with the hot reservoir and the acceptor is in the cold bath. Figure 7 demonstrates the photocell as a continuous heat engine, which is not in contact with a hot bath, but is supplied input energy by photon pulses. In our case, we have the flexibility of engineering the input photon pulses as desired. In Figure 8, we further compare the power efficiency between the two cases. Figure 8a is the case where the photon pulses are applied at a finite time interval (discrete mode operation). On the other hand, Figure 8b is the case where the photon pulses are applied almost continuously. Both have the same energy parameter

〈 n 〉 = 1

of incoming Gaussian pulses. In the discrete mode, we see an oscillatory behavior of power efficiency between 0.2 and 0.6. In the continuum mode, the power efficiency does not oscillate but asymptotically approaches the value

η_{p} \sim 0.36

η_{p} = \frac{P_{out} (t)}{P_{D} (t)} .

(34)

Let us compare the two cases in Figure 7 and Figure 8b both dealing with the continuum limit with different energies

〈 n 〉 = 10

and

〈 n 〉 = 1

, respectively. It is interesting to find that the efficiency turns out to be the same

η_{p} \sim 0.36

regardless of the pulse energy of our consideration. Of course, the transient behaviors are different as the high-energy case shows an oscillatory behavior while the low-energy case does not. Aside from details, we see that our model of heat engine offers a possibility to make an efficient quantum engine with a proper design.

5. Summary

In this paper, we studied quantum thermodynamics of two open quantum systems, the two-level system and the quantum photovoltaic model, driven by the Gaussian photon pulses. By solving the master equation with the time-dependent Hamiltonian of the Gaussian photon pulses, we calculated quantum thermodynamic quantities. For the two-level system in the cold bath, we examined the first law of quantum thermodynamics, which relates the energy change of the system, the heat current, and the power. We also illustrated the second law of thermodynamics by confirming that the entropy production is positive.

More importantly, we investigated the quantum photovoltaic cell in the cold bath driven by the sequence of the Gaussian photon pulses. The power efficiency of the quantum photocell was considered as the ratio of the output power delivered to the external load by the photocell to the input power delivered by the photon pulses. We showed that the quantum photocell as a heat engine can operate both in the discrete stroke mode and in the continuous stroke mode by changing the sequence of the photon pulses.

Our model of quantum heat engine based on a driven quantum system in contact with a single bath seems worthwhile to further investigate. In our work, we showed that the efficiency as high as

η_{p} = 0.36

can be achieved, which should be further explored in a broad range of system parameters. The maximum efficiency of the quantum photocell may be obtained by applying the optimal control method [35] or by adopting the quasi-static or quantum adiabatic processes [36,37] There are some meaningful directions to consider. One is to study how the dark state or the quantum coherence can further enhance the performance of the photocell. We also note that recently Chan et al. [31] studied the quantum dynamics of excitons by absorption of single photons in photosynthetic light-harvesting complexes. It would be interesting how the photosynthetic light-harvesting complexes behave when the photon pulses are applied. Moreover, while we considered the Gaussian photon pulses in the current work, other photon pulses, for example, hyperbolic secant, rectangular, or symmetric exponential pulses may be tested to come up with an optimal design [30]. An open problem is how to mimic the thermal photon from the hot thermal bath and to incorporate the thermal photons into the simulation. The quantum photocells considered here can be simulated on quantum computers [38,39] or by using transmon qubits as a working substance [40].

Author Contributions

Conceptualization, S.O.; Investigation, Writing—Review & editing, S.O., J.J.P. and H.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

Shockley, W.; Queisser, H.J. Detailed Balance Limit of Efficiency of p-n Junction Solar Cells. J. Appl. Phys. 1961, 32, 510–519. [Google Scholar] [CrossRef]
Scully, M.O. Quantum Photocell: Using Quantum Coherence to Reduce Radiative Recombination and Increase Efficiency. Phys. Rev. Lett. 2010, 104, 207701. [Google Scholar] [CrossRef]
Scully, M.O.; Chapin, K.R.; Dorfman, K.E.; Kim, M.B.; Svidzinsky, A. Quantum heat engine power can be increased by noise-induced coherence. Proc. Natl. Acad. Sci. USA 2011, 108, 15097–15100. [Google Scholar] [CrossRef] [Green Version]
Svidzinsky, A.A.; Dorfman, K.E.; Scully, M.O. Enhancing photovoltaic power by Fano-induced coherence. Phys. Rev. A 2011, 84, 053818. [Google Scholar] [CrossRef] [Green Version]
Creatore, C.; Parker, M.A.; Emmott, S.; Chin, A.W. Efficient Biologically Inspired Photocell Enhanced by Delocalized Quantum States. Phys. Rev. Lett. 2013, 111, 253601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Zhang, Y.; Oh, S.; Alharbi, F.H.; Engel, G.S.; Kais, S. Delocalized quantum states enhance photocell efficiency. Phys. Chem. Chem. Phys. 2015, 17, 5743–5750. [Google Scholar] [CrossRef] [Green Version]
Fruchtman, A.; Gómez-Bombarelli, R.; Lovett, B.W.; Gauger, E.M. Photocell Optimization Using Dark State Protection. Phys. Rev. Lett. 2016, 117, 203603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Higgins, K.D.B.; Lovett, B.W.; Gauger, E.M. Quantum-Enhanced Capture of Photons Using Optical Ratchet States. J. Phys. Chem. C 2017, 121, 20714–20719. [Google Scholar] [CrossRef] [Green Version]
Dorfman, K.E.; Voronine, D.V.; Mukamel, S.; Scully, M.O. Photosynthetic reaction center as a quantum heat engine. Proc. Natl. Acad. Sci. USA 2013, 110, 2746–2751. [Google Scholar] [CrossRef] [Green Version]
Killoran, N.; Huelga, S.F.; Plenio, M.B. Enhancing light-harvesting power with coherent vibrational interactions: A quantum heat engine picture. J. Chem. Phys. 2015, 143, 155102. [Google Scholar] [CrossRef]
Yamada, Y.; Yamaji, Y.; Imada, M. Exciton Lifetime Paradoxically Enhanced by Dissipation and Decoherence: Toward Efficient Energy Conversion of a Solar Cell. Phys. Rev. Lett. 2015, 115, 197701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Stones, R.; Hossein-Nejad, H.; van Grondelle, R.; Olaya-Castro, A. On the performance of a photosystem II reaction centre-based photocell. Chem. Sci. 2017, 8, 6871–6880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Kosloff, R.; Levy, A. Quantum Heat Engines and Refrigerators: Continuous Devices. Annu. Rev. Phys. Chem. 2014, 65, 365–393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Uzdin, R.; Levy, A.; Kosloff, R. Equivalence of Quantum Heat Machines, and Quantum-Thermodynamic Signatures. Phys. Rev. X 2015, 5, 031044. [Google Scholar] [CrossRef] [Green Version]
Scovil, H.E.D.; Schulz-DuBois, E.O. Three-Level Masers as Heat Engines. Phys. Rev. Lett. 1959, 2, 262–263. [Google Scholar] [CrossRef]
Oh, S. Efficiency and power enhancement of solar cells by dark states. Phys. Lett. A 2019, 383, 125857. [Google Scholar] [CrossRef]
Tomasi, S.; Kassal, I. Classification of Coherent Enhancements of Light-Harvesting Processes. J. Phys. Chem. Lett. 2020, 11, 2348–2355. [Google Scholar] [CrossRef] [Green Version]
Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
Gorini, V.; Kossakowski, A.; Sudarshan, E.C.G. Completely positive dynamical semigroups of N-level systems. J. Math. Phys. 1976, 17, 821–825. [Google Scholar] [CrossRef]
Alicki, R.; Lendi, K. Quantum Dynamical Semigroup and Applications; Springer: Berlin, Germny, 1984. [Google Scholar]
Breuer, H.P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: New York, NY, USA, 2002. [Google Scholar]
Rivas, Á.; Huelga, S.F. Open Quantum Systems, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar] [CrossRef] [Green Version]
Spohn, H. Entropy production for quantum dynamical semigroups. J. Math. Phys. 1978, 19, 1227–1230. [Google Scholar] [CrossRef]
Spohn, H.; Lebowitz, J.L. Irreversible Thermodynamics for Quantum Systems Weakly Coupled to Thermal Reservoirs. In Advances in Chemical Physics; Rice, S.A., Ed.; John Wiley & Sons, Ltd.: New York, NY, USA, 1978; Volume XXXVIII, pp. 109–142. [Google Scholar] [CrossRef]
Alicki, R. The quantum open system as a model of the heat engine. J. Phys. Math. Gen. 1979, 12, L103–L107. [Google Scholar] [CrossRef]
Das, S.; Khatri, S.; Siopsis, G.; Wilde, M.M. Fundamental limits on quantum dynamics based on entropy change. J. Math. Phys. 2018, 59, 012205. [Google Scholar] [CrossRef] [Green Version]
Jaynes, E.T.; Cummings, F.W. Comparison of quantum and semiclassical radiation theories with application to the beam maser. Proc. IEEE 1963, 51, 89–109. [Google Scholar] [CrossRef] [Green Version]
Orszag, M. Quantum Optics; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
Loudon, R. The Quantum Theory of Light; Clarendon Press: Oxford, UK, 1973. [Google Scholar]
Wang, Y.; Minář, J.C.V.; Sheridan, L.; Scarani, V. Efficient excitation of a two-level atom by a single photon in a propagating mode. Phys. Rev. A 2011, 83, 063842. [Google Scholar] [CrossRef] [Green Version]
Chan, H.C.H.; Gamel, O.E.; Fleming, G.R.; Whaley, K.B. Single-photon absorption by single photosynthetic light-harvesting complexes. J. Phys. B Atomic Mol. Opt. Phys. 2018, 51, 054002. [Google Scholar] [CrossRef] [Green Version]
Carmichael, H. An Open Systems Approach to Quantum Optics; Springer: Berlin, Germany, 1993. [Google Scholar]
Mandel, L.; Wolf, E. Optical Coherence and Quantum Optics; Cambridge University Press: New York, NY, USA, 1995. [Google Scholar]
Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes 3rd Edition: The Art of Scientific Computing, 3rd ed.; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
Stefanatos, D. Optimal efficiency of a noisy quantum heat engine. Phys. Rev. E 2014, 90, 012119. [Google Scholar] [CrossRef] [Green Version]
Alvarado Barrios, G.; Peña, F.J.; Albarrán-Arriagada, F.; Vargas, P.; Retamal, J.C. Quantum Mechanical Engine for the Quantum Rabi Model. Entropy 2018, 20, 767. [Google Scholar] [CrossRef] [Green Version]
Peña, F.J.; Zambrano, D.; Negrete, O.; De Chiara, G.; Orellana, P.A.; Vargas, P. Quasistatic and quantum-adiabatic Otto engine for a two-dimensional material: The case of a graphene quantum dot. Phys. Rev. E 2020, 101, 012116. [Google Scholar] [CrossRef] [Green Version]
Potočnik, A.; Bargerbos, A.; Schröder, F.A.Y.N.; Khan, S.A.; Collodo, M.C.; Gasparinetti, S.; Salathé, Y.; Creatore, C.; Eichler, C.; Türeci, H.E.; et al. Studying light-harvesting models with superconducting circuits. Nat. Commun. 2018, 9, 904. [Google Scholar] [CrossRef]
Maslennikov, G.; Ding, S.; Hablützel, R.; Gan, J.; Roulet, A.; Nimmrichter, S.; Dai, J.; Scarani, V.; Matsukevich, D. Quantum absorption refrigerator with trapped ions. Nat. Commun. 2019, 10, 202. [Google Scholar] [CrossRef]
Cherubim, C.; Brito, F.; Deffner, S. Non-Thermal Quantum Engine in Transmon Qubits. Entropy 2019, 21, 545. [Google Scholar] [CrossRef] [Green Version]

Figure 1. A two-level system with energy levels

E_{0}

and

E_{1}

in contact with a cold thermal bath at

T_{c}

is driven by Gaussian photon pulses serving as an energy source in our work.

Figure 1. A two-level system with energy levels

E_{0}

and

E_{1}

in contact with a cold thermal bath at

T_{c}

is driven by Gaussian photon pulses serving as an energy source in our work.

Figure 2. (a) The density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

are plotted over time. (b) The rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

are calculated as functions of time. (c) The energy

E (t)

, the work

W (t)

, the heat transfer

Q (t)

, and the system entropy

S (t)

are plotted as functions of time. (d) The rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted over time. The parameters are taken as

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 2. (a) The density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

are plotted over time. (b) The rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

are calculated as functions of time. (c) The energy

E (t)

, the work

W (t)

, the heat transfer

Q (t)

, and the system entropy

S (t)

are plotted as functions of time. (d) The rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted over time. The parameters are taken as

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 3. (a) The density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

are plotted over time. (b) The rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

are plotted as functions of time. (c) The energy

E (t)

, the work

W (t)

, the heat transfer

Q (t)

, and the system entropy

S (t)

are calculated as functions of time. (d) The rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are shown over time. The parameters are

〈 n 〉 = 10

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 3. (a) The density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

are plotted over time. (b) The rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

are plotted as functions of time. (c) The energy

E (t)

, the work

W (t)

, the heat transfer

Q (t)

, and the system entropy

S (t)

are calculated as functions of time. (d) The rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are shown over time. The parameters are

〈 n 〉 = 10

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 4. When a Gaussian photon pulse is overlapped with the subsequent Gaussian photon pulse and the interval between them is regular, (a) the density matrix elements of the two-level system, (b) the rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, the heat current

J (t)

, (c) energy

E (t)

, work

W (t)

, heat

Q (t)

system entropy

S (t)

, (d) the rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted as a function of time. The parameters are taken as

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

\frac{d E (t)}{d t}

, the power

P (t)

, the heat current

J (t)

, (c) energy

E (t)

, work

W (t)

, heat

Q (t)

system entropy

S (t)

, (d) the rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted as a function of time. The parameters are taken as

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 5. When an irregularly spaced sequence of photon pulses is applied, (a) the density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

, (b) the rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

, (c) the energy

E (t)

, work

W (t)

, and heat

Q (t)

, and the system entropy

S (t)

, (d) the rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted as functions of time. Parameters:

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 5. When an irregularly spaced sequence of photon pulses is applied, (a) the density matrix elements of the two-level system and the sequence of Gaussian photon pulses

g (t)

, (b) the rate of energy change

\frac{d E (t)}{d t}

, the power

P (t)

, and the heat current

J (t)

, (c) the energy

E (t)

, work

W (t)

, and heat

Q (t)

, and the system entropy

S (t)

, (d) the rate of system entropy change

\frac{d S}{d t}

and the entropy production

σ (t)

are plotted as functions of time. Parameters:

〈 n 〉 = 1

γ = 10^{- 2} ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = 1 eV

Figure 6. Schematic diagram of a donor-acceptor photocell.

γ_{01}

is the spontaneous decay due to the coupling with the cold thermal bath.

γ_{21}

and

γ_{03}

are the transfer rate between the donor and the acceptor.

Γ

stands for the external load or electrical resistance.

Figure 6. Schematic diagram of a donor-acceptor photocell.

γ_{01}

is the spontaneous decay due to the coupling with the cold thermal bath.

γ_{21}

and

γ_{03}

are the transfer rate between the donor and the acceptor.

Γ

stands for the external load or electrical resistance.

Figure 7. (a) The diagonal matrix elements of the density operator of the photocell and the pulse profiles are plotted as a function of time. (b) The changes in the energy of the donor and acceptor

{(d E / d t)}_{D} = E_{D}^{'}

and

{(d E / d t)}_{A} = E_{A}^{'}

, the power delivered to the donor

P_{D} (t)

by the photon pulses, the power output

P_{out}

, and the heat currents of the donor and acceptor

J_{D}

and

J_{A}

are plotted as a function of time. (c) The entropy of the quantum photocell,

S (t)

, the entropy of the donor,

S_{D} (t)

, and the entropy of the acceptor,

S_{A} (t)

, are calculated as a function of time. (d) The current

I (t)

, the voltage

V (t)

, and the efficiency

η

are plotted as a function of time. Parameters:

〈 n 〉 = 10

γ_{12} = γ_{03} = 10^{- 2} ω_{0}

γ_{01} = 10^{- 2} ω_{0}

Γ = 0.1 ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = E_{1} - E_{0} = 1.8 eV

{(d E / d t)}_{D} = E_{D}^{'}

and

{(d E / d t)}_{A} = E_{A}^{'}

, the power delivered to the donor

P_{D} (t)

by the photon pulses, the power output

P_{out}

, and the heat currents of the donor and acceptor

J_{D}

and

J_{A}

are plotted as a function of time. (c) The entropy of the quantum photocell,

S (t)

, the entropy of the donor,

S_{D} (t)

, and the entropy of the acceptor,

S_{A} (t)

, are calculated as a function of time. (d) The current

I (t)

, the voltage

V (t)

, and the efficiency

η

are plotted as a function of time. Parameters:

〈 n 〉 = 10

γ_{12} = γ_{03} = 10^{- 2} ω_{0}

γ_{01} = 10^{- 2} ω_{0}

Γ = 0.1 ω_{0}

Ω = ω_{0} / 4 π

, and

ℏ ω_{0} = E_{1} - E_{0} = 1.8 eV

Figure 8. (a) From discrete mode to (b) the continuous mode operation by changing the interval of the pulses. Parameters:

〈 n 〉 = 1

γ_{21} = γ_{03} = 10^{- 3} ω_{0}

γ_{01} = 10^{- 2} ω_{0}

, and

Ω = ω_{0} / 4 π

Figure 8. (a) From discrete mode to (b) the continuous mode operation by changing the interval of the pulses. Parameters:

〈 n 〉 = 1

γ_{21} = γ_{03} = 10^{- 3} ω_{0}

γ_{01} = 10^{- 2} ω_{0}

, and

Ω = ω_{0} / 4 π

Table 1. Typical parameters used in this work.

Energy gap of the two-level system	$E_{1} - E_{0} = ℏ ω_{0} = 1.0$ eV
Energy gap of the donor of the quantum photocell	$E_{1} - E_{0} = ℏ ω_{0} = 1.8$ eV
Energy gap of the acceptor of the quantum photocell	$E_{2} - E_{3} = 1.6$ eV
Weisskopf-Winger constant	$γ / ω_{0} = γ_{01} / ω_{0} = 10^{- 3} \sim 10^{- 6}$
Phonon decay constant	$γ_{12} / ω_{0} = γ_{03} / ω_{0} = 10^{- 2} \sim 10^{- 3}$
Photon number of a pulse	$〈 n 〉 = {\| α \|}^{2} = 1 or 10$
Temperature of the cold bath	$T_{c}$ = 300 K
Width of a Gaussian pulse	$Ω = ω_{0} / 4 π$

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Oh, S.; Park, J.J.; Nha, H. Quantum Photovoltaic Cells Driven by Photon Pulses. Entropy 2020, 22, 693. https://doi.org/10.3390/e22060693

AMA Style

Oh S, Park JJ, Nha H. Quantum Photovoltaic Cells Driven by Photon Pulses. Entropy. 2020; 22(6):693. https://doi.org/10.3390/e22060693

Chicago/Turabian Style

Oh, Sangchul, Jung Jun Park, and Hyunchul Nha. 2020. "Quantum Photovoltaic Cells Driven by Photon Pulses" Entropy 22, no. 6: 693. https://doi.org/10.3390/e22060693

APA Style

Oh, S., Park, J. J., & Nha, H. (2020). Quantum Photovoltaic Cells Driven by Photon Pulses. Entropy, 22(6), 693. https://doi.org/10.3390/e22060693

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Quantum Photovoltaic Cells Driven by Photon Pulses

Abstract

1. Introduction

2. Quantum Thermodynamics of Open Quantum Systems

3. A Two-Level System Driven by Photon Pulses

4. Quantum Photocell Driven by Photon Pulses

5. Summary

Author Contributions

Funding

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI