Next Article in Journal
Reaction of Partially Methylated Polygalacturonic Acid with Iron(III) Chloride and Characterization of a New Mixed Chloride–Polygalacturonate Complex
Previous Article in Journal
Synthesis and Dynamic Behavior of Ce(IV) Double-Decker Complexes of Sterically Hindered Phthalocyanines
Previous Article in Special Issue
Exploring Spin Distribution and Electronic Properties in FeN4-Graphene Catalysts with Edge Terminations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery

Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(4), 889; https://doi.org/10.3390/molecules29040889
Submission received: 28 December 2023 / Revised: 7 February 2024 / Accepted: 13 February 2024 / Published: 17 February 2024
(This article belongs to the Special Issue Molecular Simulations of Energy Materials)

Abstract

:
In a previous study, we proposed an open quantum network model of a quantum battery (QB) that possesses dark states owing to its structural exchange symmetries. While in a dark state, the QB is capable of storing an exciton without any environment-induced population losses. However, when the structural exchange symmetry is broken, the QB begins to discharge the exciton towards its exit site. In this article, we start by demonstrating that this QB is not only loss-free with respect to exciton population during the storage phase, but also with respect to the QB energy. We then explore the exciton population and energy transfer dynamics of the QB during the discharge phase over a wide range of site energies, bath temperatures, and bath reorganization energies. Our results shed light on how to optimize the QB’s population and energy transfer dynamics for different purposes.

1. Introduction

A quantum battery (QB) is a quantum system that can store and release energy as needed. Over the past decade, several types of QBs have been proposed, demonstrating advantages with respect to charging speed [1,2,3,4,5,6,7,8,9] and work extraction [5,7,10,11,12,13,14] over the classical analogues. In practice, a QB may interact with an environment, which could lead to losses that negatively impact its performance. Thus, any theoretical study should treat the QB as an open quantum system [15,16,17]. Previously, open quantum network (OQN) models, in which the network sites (representing quantum systems) may be coupled to each other and to dissipative/decohering environments, have been used to study the dynamics of open QBs [5,6,8,18].
Over the years, a number of ways of protecting a quantum system from its environment have been proposed. One well-known approach is based on the use of a decoherence-free subspace (DFS) [19,20]—a subspace of the Hilbert space in which the dynamics is purely unitary. While in a DFS, the system’s dynamics is dissipationless and decoherence-free despite the coupling to an environment. As has been shown in Refs. [21,22,23], if there exists a unitary operator that commutes with all elements in the system’s master equation exists, then the system will possess invariant subspaces. Among such subspaces, the one-dimensional subspaces are DFSs because the dynamics maps them onto themselves. Owing to this property, DFSs have many potential applications in quantum information and quantum computing [24,25,26,27,28,29].
Recently, researchers have used dark states living in DFSs for stabilizing and enhancing the performance of QBs and other OQNs [18,30,31,32]. In particular, we proposed an OQN model of a QB with site exchange symmetries that support the existence of dark states localized on the bulk sites of the network [18]. We showed that it is possible to store an exciton in one of these dark states, without any exciton population transfer to the two surface sites which are connected to thermal baths at equal temperatures. In this way, it was possible to protect the QB from environment-induced excitation energy losses. Moreover, we showed that by attaching an additional bath to break the exchange symmetry of the QB, it is possible to discharge the exciton towards the designated exit site to be ultimately harnessed by a sink.
Although it was shown that our QB model is capable of operating as an excitonic QB (i.e., a QB that stores and discharges excitons), the question of whether or not the model is capable of operating as an energy battery (i.e., a QB whose energy does not decrease significantly during the storage phase and discharges energy to a sink with minimal loss to the baths) remained to be explored. Therefore, in this study, in addition to monitoring the exciton population dynamics, we calculate the various contributions to the total energy of the system during both the storage and discharge phases for a wide range of bath temperature gaps, bath reorganization energies, and site energies. The aim of the study is to identify the conditions which optimize the performance of the QB model for its use as an excitonic QB, energy QB, or both.
The paper is organized as follows. First, in Section 2, we introduce the open QB model. In Section 3, we explain how the dynamics of the model are simulated and how the populations/energies are calculated. In Section 4, we present the time-dependent populations/energies of the model for different bath temperatures, bath reorganization energies, and site energies. The performance of the QB under the various conditions is also discussed. In Section 5, we summarize our findings.

2. Open QB Model

Following Ref. [18], we consider the same six-site para-benzene-shaped OQN (the numbering of the sites is depicted in Figure 1b). Setting = 1 , the Hamiltonian of the closed network is given by
H ^ N = n = 1 6 E n | n n | + h n , m | n m | ,
where | n corresponds to a singly excited state localized on site n, E n is the energy of site n, h is the electronic coupling strength between sites n and m, and n , m denotes that a cyclic summation over nearest-neighbour sites is performed. To construct the OQN, the two para-sites of the network are coupled to thermal baths, each composed of a set of independent harmonic oscillators. The sum of the bath and network–bath coupling Hamiltonians is given by
H ^ B + H ^ N B = 1 2 n S S s j M P ^ n , j 2 + ω n , j 2 R ^ n , j C n , j ω n , j 2 | n n | 2 ,
where M is the number of harmonic oscillators in each bath, P ^ n , j and R ^ n , j are the mass-weighted momentum and position operators, respectively, of the jth oscillator with frequency ω n , j , and C n , j is the network–bath coupling strength between the nth site and the jth oscillator. The sites coupled to the baths are referred to as surface sites (SSs) while the remaining sites are referred to as bulk sites (BSs). In this study, we take the site energies of all BS to be equal to a value E B S , while considering different values of site energies of the SS.
The network defined by H ^ N and the aforementioned site energies possesses the following unitary symmetry operator [33]
Π ^ = | 1 1 | + | 4 4 | + | 2 6 | + | 3 5 | + h . c . ,
which satisfies
[ Π ^ , H ^ N ] = 0 , [ Π ^ , H ^ N B ] = 0 n SSs .
As a result, Π ^ shares the same eigenstates as H ^ N . Due to the existence of this symmetry operator, the system possesses two DFSs with the following DSs ( | ψ α ) and eigenvalues ( u α ) [22,34]
| ψ 1 = 1 2 | 5 + | 6 | 2 | 3 , u 1 = E B S + h | ψ 2 = 1 2 | 3 + | 6 | 2 | 5 , u 2 = E B S h
If the system is initialized in the dark state | ψ α , it will undergo dissipationless dynamics, i.e., the dark state will be invariant under the effect of the evolution operator of the composite system and the site populations will be
P ^ n n ( t ) = P ^ n n ( 0 ) = 0 , n SSs , 1 4 , n BSs ,
where P ^ n n = | n n | is the projection operator corresponding to site n and · denotes an ensemble average. Because an exciton can be stored indefinitely in this state, this phase is termed the storage phase.
When a symmetry-breaking perturbation (SBP) is connected to the OQN, the symmetry operator and DSs no longer persist and the QB may begin to discharge a stored exciton. In this study, this is achieved by attaching a bath of M harmonic oscillators to sites 2 and 3 simultaneously. The sum of the SBP and network-SBP coupling Hamiltonians is
H ^ P + H ^ N P = 1 2 k M p ^ k 2 + Ω k 2 r ^ k γ k Ω k 2 S ^ 2 ,
where p ^ k , r ^ k , Ω k , and γ k are the momentum operator, position operator, frequency, and network–SBP coupling strength of the kth oscillator, respectively, and S ^ = | 2 2 | + | 3 3 | . This phase is termed the discharge phase.

3. Simulation Details

Due to the large number of degrees of freedom in the composite system (i.e., QB, thermal baths, and SBP), a fully quantum dynamical simulation of the composite system would be computationally expensive. Thus, following Ref. [18], we use a mixed quantum–classical dynamics method known as “Deterministic evolution of coordinates with initial decoupled equations” (DECIDE) [35,36], which treats the OQN quantum mechanically and the baths and SBP in a classical-like way. Previously, the DECIDE method has been successfully applied to a host of model systems over a large range of parameter regimes [35,37,38,39,40], and is, therefore, expected to produce reliable results in this study. That being said, DECIDE may yield inaccurate results for systems with very slow thermal baths (i.e., when the bath cut-off frequency is much smaller than the subsystem energy gaps) or with very low bath temperatures, neither of which is the case in the present study.
To apply the DECIDE method, we must first apply the Wigner transform [41] to the bath and SBP degrees of freedom. The resulting partially Wigner-transformed Hamiltonian of the composite system is
H ^ W = H ^ N + H ^ N B ( { R n , j } ) + H B ( { P n , j } , { R n , j } ) + χ H P ( { p k } , { r k } ) + H ^ N P ( { r k } ) ,
where { R n , j , P n , j } and { r k , p k } are the position and momentum variables of the baths and SBP, respectively. The parameter χ is equal to 1 when the SBP is attached to the OQN and 0 otherwise. The coordinates of the OQN are taken to be P ^ n m = | n m | , while the coordinates of the baths and SBP are their positions and momenta. According to the DECIDE method, the coupled equations of motion for all of the coordinates are given by [18]
d d t P n m β β ( t ) = i l = 1 6 V l n P ^ l m ( t ) v = 1 6 V m v P ^ n v ( t ) β β i 2 j C n , j ( R n , j ( t ) P ^ n m ( t ) + P ^ n m ( t ) R n , j ( t ) ) β β ( δ n , 1 + δ n , 4 ) + i 2 j C m , j ( R m , j ( t ) P ^ n m ( t ) + P ^ n m ( t ) R m , j ( t ) ) β β ( δ m , 1 + δ m , 4 ) χ i 2 k C k ( R k ( t ) P ^ n m ( t ) + P ^ n m ( t ) R k ( t ) ) β β ( δ n , 2 + δ n , 3 ) + χ i 2 k C k ( R k ( t ) P ^ n m ( t ) + P ^ n m ( t ) R k ( t ) ) β β ( δ m , 2 + δ m , 3 ) , d d t R n , j β β ( t ) = P n , j β β ( t ) , d d t P n , j β β ( t ) = ω n , j 2 R n , j β β ( t ) + C n , j P n n β β ( t ) ( δ n , 1 + δ n , 4 ) , d d t r k β β ( t ) = p k β β ( t ) , d d t p k β β ( t ) = Ω k 2 r k β β ( t ) + χ γ k ( P ^ 22 + P ^ 33 ) β β ( t ) ,
where V n n = E n + j = 1 M C n , j 2 / ( 2 ω n , j 2 ) ( δ n , 1 + δ n , 4 ) + χ k = 1 M γ k 2 / ( 2 Ω k 2 ) ( δ n , 2 + δ n , 3 ) , and when n m , V n m = h for | n m | = 1 and 0 otherwise. Here, β labels an arbitrary basis state, i.e., the matrix element of P ^ n m is given by P n m β β = β | P ^ n m | β . In this work, { | β } = { | 1 , | 2 , , | 6 } .
We assume the initial state of the composite system to be factorized, i.e., ρ ^ t o t ( 0 ) = ρ ^ N ( 0 ) ρ B , W ( 0 ) ρ P , W ( 0 ) , where ρ ^ N ( 0 ) is the initial density operator of the network, and ρ B , W ( 0 ) and ρ P , W ( 0 ) are the initial Wigner-transformed densities of the heat baths and SBP, respectively (N.B.: ρ P , W ( 0 ) is omitted when the OQN is not attached to the SBP). The initial state of the network is taken to be the dark state ρ ^ N ( 0 ) = | ψ 1 ψ 1 | , where | ψ 1 is defined in Equation (5). The initial values of the OQN coordinates are always taken to be P n m β β = δ β , n δ m , β . The bath oscillators are initialized in the thermal equilibrium state given by (setting k B = 1 ) [42]
ρ B , W ( 0 ) = n = 1 , 4 j = 1 M tanh ( β ω n , j / 2 ) π exp 2 tanh ( β ω n , j / 2 ) ω n , j P n , j 2 2 + ω n , j 2 R n , j 2 2 ,
where β = 1 / k B T is the inverse temperature. The oscillators of the SBP are also initialized in a thermal equilibrium state with an analogous form. The initial positions and momenta of the bath and SBP oscillators are sampled from Equation (10) and its analog for the SBP, respectively. The system–bath and system–SBP couplings are characterized by a Debye–Drude spectral density, i.e., J ( ω n , j ) = 2 λ b ω n , j ω n , j 2 + 1 . In this work, the spectral density is discretized to yield the following expressions for the coupling strengths C n , j and frequencies ω n , j [43,44]:
ω n , j = tan ( j arctan ( ω m a x / ω c ) / M ) ω c ,
C n , j = 2 λ b arctan ( ω m a x / ω c ) / ( π M ) ω n , j ,
where λ b is the bath reorganization energy and ω c is the bath cut-off frequency [45].
Previously, the fourth-order Runge–Kutta method was used to integrate the DECIDE equations of motion in Equation (9), yielding conserved total populations for the system under study [18,35]. However, in this work, we found that a combination of a smaller time step and a higher order integrator is needed for good energy conservation and more accurate calculations of the various contributions to the total energy. High-order methods such as the eighth-order Runge–Kutta method [46] can yield accurate results with a relatively large time step, but it contains many integration stages which increase the simulation time drastically. Considering the trade-off between time step and number of integration steps, we employed the sixth-order Runge–Kutta method in this work. Using this integrator with a time step of 0.16 fs, the total energy drift is less than 10 2 cm 1 over a 1 ps trajectory [see Section S1 of the Supporting Information (SI)].
The time-dependent population of site n is calculated via an ensemble average of the projection operator P ^ n n , viz.,
P ^ n n ( t ) = β β d X ( 0 ) P n n β β ( t ) ρ N β β ( 0 ) ρ E , W ( 0 ) ,
where X ( 0 ) = ( { R n , j } , { P n , j } , { r k } , { p k } ) are the initial coordinates of the baths and SBP. Similarly, the average total energy of the composite system is
E t o t ( t ) = β β d X ( 0 ) H ^ W β β ( t ) ρ N β β ( 0 ) ρ E , W ( 0 ) .
In the above equations, ρ E , W = ρ B , W in the absence of the SBP and ρ E , W = ρ B , W ρ P , W when the SBP is attached to the QB. Using our numerical results, we have verified that n = 1 6 P ^ n n ( t ) 1 (i.e., population conservation) and d d t E t o t ( t ) 0 (i.e., energy conservation). For the purposes of our analysis, the total energy of the composite system may be decomposed into the following contributions: OQN energy ( E N ), bath energy ( E N B + E B ), and SBP energy ( E N P + E P ). All simulation results are averaged over 10,000 trajectories, which ensures that the error bars are much smaller than the symbols in the figures.

4. Results and Discussion

4.1. Exciton Storage and Discharge

We start by considering the time-dependent site populations and energies of the OQN, baths, and SBP in the storage ( χ = 0 ) and discharge ( χ = 1 ) phases, using the parameter values from Ref. [18], viz., E 1 = 250 cm 1 , E i { 2 , 3 , 5 , 6 } = 200 cm 1 , E 4 = 0 cm 1 (N.B.: E 4 is smaller than the energies of the remaining sites because site 4 is designated as the exit site, which in practice would be attached to a sink that captures the exciton/energy), h = 60 cm 1 , T L = T R = 300 K, T p = 300 K, λ b = 35 cm 1 , λ p = 10 cm 1 , ω c = ω p = 106 cm 1 , ω m a x = 50 ω c , and M = 100 .
In Figure 1a, we see that the BS and SS populations remain constant at 0.25 and 0, respectively, during the storage phase. Simultaneously, as seen in Figure 1c, the OQN, left-bath, and right-bath energies remain constant (the minor deviations in the bath energies are attributed to numerical errors). Thus, our QB model is capable of perfectly storing both population and energy during the storage phase, as predicted by the theory. In addition, as shown in Section S3 of the SI, the coherences remain constant during the storage phase.
After attaching the SBP, the populations of sites 1 and 4 increase from 0 while the remaining site populations decrease from 0.25 over the 1 ps time period, as seen in Figure 1b. At t = 1 ps, the population of site 4 (i.e., the exit site) is greater than that of site 1 and the BS. (The corresponding time-dependent coherences are shown in Figure S2 of the SI). In Figure 1d, we plot the changes in the OQN, left-bath, right-bath, and SBP energies with respect to their initial values. As can be seen, the OQN and left-bath energies decrease while the SBP and right-bath energies increase. Thus, during the discharge phase, energy flows from the OQN and left bath into the right bath and SBP, with considerably more energy flowing to the right bath than to the SBP.
To further analyze the OQN energy ( E N ) change in the discharge phase, we decompose it into the on-site energy and exchange energy, corresponding to the terms n = 1 6 E n | n n | and h n , m | n m | , respectively, in Equation (1). The change in the exchange energy is therefore given by
Δ E exch ( t ) = h n , m P ^ n m ( t ) P ^ n m ( 0 ) .
At t = 0 , only P ^ 23 , P ^ 32 , P ^ 56 , and P ^ 65 are non-zero (since the OQN is initialized in the dark state | ψ 1 ). Their values decrease from 0.25 to 0.07 over the course of 1 ps (Figure S4), which causes Δ E exch to increase. For the dynamics displayed in Figure 1b,d, Δ E exch increases by 25.59 cm 1 over the course of 1 ps (as calculated by Equation (15)). In fact, Δ E exch increases for all of the parameter regimes studied herein. Considering the values of P ^ n n ( 0 ) and the fact that n = 1 6 P ^ n n ( t ) = 1 , the on-site energy change is (setting all BS energies to be equal, i.e., E i { 2 , 3 , 5 , 6 } = E B S )
Δ E on - site ( t ) = n = 1 6 E n P ^ n n ( t ) n = 1 6 E n P ^ n n ( 0 ) = ( E 1 E B S ) P ^ 11 ( t ) + ( E 4 E B S ) P ^ 44 ( t ) .
From this expression, we see that the on-site energy change depends only on the time-dependent populations of the SS and the energy differences E 1 E B S and E 4 E B S (rather than the absolute values of E 1 and E 4 ). [Indeed, shifting all of the site energies by the same constant results in the same dynamics during the discharge phase]. When site 4 is chosen to be the exit site (i.e., E 4 < E B S < E 1 ), a higher site 4 population will lead to a more negative on-site energy change. In fact, for the dynamics displayed in Figure 1b,d, Δ E on - site decreases by 50.88 cm 1 over the course of 1 ps (as calculated by Equation (16)), which is greater than the increase in Δ E exch . This explains the decrease in the OQN energy observed in Figure 1d.

4.2. Effect of Bath Temperature

We next investigate the effect of varying the bath temperature gap on the site populations and energy flow in the QB, in an effort to find parameter sets that maximize the population of the exit site and minimize the energy loss of the OQN. To simplify our exploration of the parameter space, we vary one parameter at a time while keeping the remaining parameters fixed. First, we fix the temperature of the left bath at 300 K and vary the temperature of the right bath, while keeping the OQN, SBP, and remaining bath parameters unchanged.
In Figure 2a,c, we see that increasing the right bath temperature leads to a decrease in the population of site 4 and relatively small increases in the populations of the remaining sites. As for the energy changes (see Figure 2b,d), the OQN and left bath lose energy, but the energy loss becomes smaller and remains relatively constant, respectively, with increasing right bath temperature. On the other hand, the SBP and right bath gain energy, but the energy gain remains relatively constant and becomes smaller, respectively, with increasing right-bath temperature. When we fix the temperature of the right bath at 300 K and increase the temperature of the left bath, the behaviours of the site populations are similar to those observed in the case when the temperature of the right bath is increased, except for that of site 4 which now exhibits a substantially smaller decrease (see Figure 2c). As for the energy changes (see Figure 2d), the behaviours of the SBP and OQN energies are similar to those observed in the case when the temperature of the right bath is varied. However, the bath energies are significantly different, with the energy changes in the left and right baths becoming more negative and remaining relatively constant, respectively, with increasing left-bath temperature.
Increasing the bath temperature (or decreasing β ) will increase the width of the initial Wigner distribution in Equation (10). This increase is particularly significant for oscillators with low frequencies. It can be verified both analytically and numerically that the ensemble averages of the initial bath kinetic energy, j M P n , j 2 , and bath potential energy, j M ω n , j 2 R n , j 2 , grow linearly with increasing temperature. For example, increasing the temperature of a bath from 300 K to 600 K to 900 K will increase its kinetic energy from ≈1.43 ×  10 4 eV to ≈2.38 ×  10 4 eV to ≈3.37 ×  10 4 eV.
As seen in Figure 2d, increasing the left bath or right bath temperature will lead to similar amounts of additional energy transfer out of the bath. This temperature-driven increase in energy transfer to the OQN is relatively independent of the other parameters. As seen in Figure 3, for different site 4 energies, increasing the right bath temperature by a given amount leads to roughly the same increase in energy being transferred from the bath to the OQN. In this way, one can decrease the OQN energy loss by increasing the temperature of any bath. Finally, varying the temperature of a given bath has a larger effect on the site connected to it. For example, if we increase the left bath temperature while keeping the right bath temperature constant, there will be a larger change in the site 1 population than the site 4 population (see Figure 2). These results suggest that one could design a QB that minimizes the OQN energy loss while maintaining a relatively large site 4 population.

4.3. Effect of Bath Reorganization Energy

Next, we vary the bath reorganization energy λ b . From Figure 4, we see that increasing the right-bath reorganization energy causes the site 1/4 population to decrease/increase, the right-bath energy to increase, and the system to lose more energy; the populations of the remaining sites, left-bath energy, and SBP energy remain relatively constant. On the other hand, we see that increasing the left-bath reorganization energy does not have a significant impact on the site population and energy changes.
As seen in Equation (12), increasing the bath reorganization energy λ b will increase the coupling strength C n , j . For example, if we increase λ b from 35 cm 1 to 70 cm 1 , the coupling strength for each oscillator in the bath will increase by a factor of 2 . Increasing the coupling strengths increases the magnitudes of the coupling terms in the equations of motion, which translates into faster energy transfer between the system and bath. As seen in Figure 4b, increasing the energy transfer rate of the right bath results in more energy transfer out of the bath after 1 ps. On the other hand, increasing the energy transfer rate of the left bath does not cause any significant changes in the left bath energy and site 1 population. This may be due to the relatively low population at site 1, viz., after 1 ps, the populations at sites 1 and 4 are ≈0.11 and 0.28, respectively (see Figure 1).

4.4. Effect of Site Energy

We now investigate the effects of varying the SS energies on the site populations and energy changes. Time series and results after 1 ps for different combinations of the SS energies are shown in Figure 5 and Table S1 of the SI, respectively. As seen in Figure 5a, setting the SS energies equal to each other leads to roughly equal SS populations, with the SS populations decreasing when increasing from E 1 = E 4 = 100 cm−1 to E 1 = E 4 = 250 cm−1. When the energies of the SS are greater than those of the BSs, the BS are more populated than the SSs after 1 ps. Conversely, when the energies of the SS are lower than those of the BS, the SS are more populated than the BS after 1 ps. When the energies of the SS are equal to those of the BS, we see that the SS and BS populations approach each other over time, becoming almost equal after 1 ps. With regards to the energy changes (Figure 5b), setting the SS energies equal to each other leads to roughly equal changes in the left and right bath energies, with the bath energy changes decreasing and becoming more negative when increasing from E 1 = E 4 = 100 cm−1 to E 1 = E 4 = 250 cm−1. More specifically, there are energy gains in the baths for E 1 = E 4 = 100 cm−1 and energy losses for the larger SS energies, with the loss increasing with increasing SS energy. Conversely, there is an energy loss from the system when E 1 = E 4 = 100 cm−1 and energy gains for the larger SS energies, with the gain increasing with increasing SS energy. The SBP energy remains mostly unchanged for the different SS energies.
In Figure 5c,d, we plot the results for E 1 = 250 cm−1 and different values of E 4 . As can be seen, increasing the energy of site 4 from 0 to 100 cm−1 does not cause a significant change in the site 4 population; however, increasing the energy of site 4 from 100 to 250 cm−1 causes a ≈50% drop in the site 4 population after 1 ps. Moreover, when the energy of site 4 is smaller than those of the remaining sites, site 4 becomes the most populated site after 1 ps. The populations of the remaining sites do not change significantly after going from E 4 = 0 to 100 cm−1, but they each increase by several percent after 1 ps after going from E 4 = 100 to 250 cm−1. As for the energy changes, increasing the energy of site 4 leads to more energy transfer from the right bath to the system, which in turn causes the OQN energy to change from decreasing to increasing. When both E 1 and E 4 are greater than E B S and the SS become populated, Δ E on - site also becomes positive. Thus, when Δ E exch is also positive, the OQN will gain energy from the baths.
Based on the results above, a few general observations may be made. First, higher site energies are associated with lower site populations. In addition, for a higher site energy, more energy is transferred from the bath to the OQN. As seen in Figure 5, for the highest E 4 , the (positive) OQN energy change increases and the (negative) bath energy changes decrease, while the magnitude of the SBP energy remains low and relatively constant, i.e., energy transfer from the bath to the OQN. If one changes the site energy of a particular site, then the population of that site will be mainly affected. If two sites have the same site energy, we expect them to eventually have equal populations. Finally, when all sites have the same site energy, all sites will have the same population after a sufficiently long period of time, despite starting with different initial populations (see Figure 5a).

5. Concluding Remarks

In this paper, we studied the population and energy transfer dynamics of an open quantum battery model, originally proposed in Ref. [18], over a wide range of parameter regimes. In the battery’s storage phase, we demonstrated that, in addition to no population leakage, there is no energy leakage from the battery into the attached baths. During the discharge phase, the changes in the populations and OQN energy are influenced by the bath temperatures, bath reorganization energies, and site energies. When increasing the temperature of one bath (while keeping the temperature of the other bath constant), we observed an increase in the energy transferred from that bath to the OQN. We found that the right bath (i.e., the bath connected to the exit site) exerts a larger influence on the exit site population than that exerted by the left bath (i.e., the bath connected to site 1) on the site 1 population. Moreover, when increasing the reorganization energy of the right bath, we observed an increase in the exit site’s population and a decrease in OQN energy. On the other hand, varying the reorganization energy of the left bath does not have pronounced effect on the population and energy. Regarding the site energies, when the energy of the exit site is lower than those of the BS, the OQN energy decreases in most of the parameter regimes studied. Lowering a given site energy causes the corresponding site population to increase. When the site energies are equal, the site populations reach roughly equal values after 1 ps, despite the different initial populations.
The results of our parameter space exploration show that different parameter sets render the QB conducive to different applications. In practice, this would amount to designing the QB in such a way that its properties are consistent with those of the desired parameter set. For example, for an energy battery, one may desire that the QB gains energy from its environment during the discharge phase. As we have shown, if one sets the site energies of the SS to be larger than those of the BS, then the QB gains energy as the SS populations grow in time. On the other hand, for an excitonic battery, one may desire to maximize the population of the exit site, regardless of the change in the OQN energy. In such a case, one could lower the exit site energy or lower the right-bath temperature while maintaining the left-bath temperature at its original value. If one would like to simultaneously reduce the loss in OQN energy and increase the exit site population, one could either increase the exit site energy or increase the left bath temperature. Overall, our findings shed light on design principles that could be used to construct different types of QBs operating between two thermal reservoirs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040889/s1. Total energy drift, effect of varying surface site energies, time derivative of the bath and network-bath energies, and time-dependent quantum coherences.

Author Contributions

Conceptualization, Z.L. and G.H.; methodology, Z.L. and G.H.; software, Z.L.; validation, Z.L. and G.H.; formal analysis, Z.L.; investigation, Z.L. and G.H.; resources, G.H.; writing—original draft preparation, Z.L.; writing—review and editing, Z.L. and G.H.; visualization, Z.L.; supervision, G.H.; funding acquisition, G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Sciences and Engineering Research Council of Canada (grant number RGPIN-2020-05977).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article and supplementary materials.

Acknowledgments

The authors would like to thank Mahmoud Mahdian for useful discussions during the early stages of this work.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Stefanatos, D. Optimal efficiency of a noisy quantum heat engine. Phys. Rev. E 2014, 90, 012119. [Google Scholar] [CrossRef]
  2. Binder, F.C.; Vinjanampathy, S.; Modi, K.; Goold, J. Quantacell: Powerful charging of quantum batteries. New J. Phys. 2015, 17, 075015. [Google Scholar] [CrossRef]
  3. Campaioli, F.; Pollock, F.A.; Binder, F.C.; Céleri, L.; Goold, J.; Vinjanampathy, S.; Modi, K. Enhancing the Charging Power of Quantum Batteries. Phys. Rev. Lett. 2017, 118, 150601. [Google Scholar] [CrossRef]
  4. Ferraro, D.; Campisi, M.; Andolina, G.M.; Pellegrini, V.; Polini, M. High-Power Collective Charging of a Solid-State Quantum Battery. Phys. Rev. Lett. 2018, 120, 117702. [Google Scholar] [CrossRef]
  5. Barra, F. Dissipative Charging of a Quantum Battery. Phys. Rev. Lett. 2019, 122, 210601. [Google Scholar] [CrossRef]
  6. Farina, D.; Andolina, G.M.; Mari, A.; Polini, M.; Giovannetti, V. Charger-mediated energy transfer for quantum batteries: An open-system approach. Phys. Rev. B 2019, 99, 035421. [Google Scholar] [CrossRef]
  7. Rossini, D.; Andolina, G.M.; Polini, M. Many-body localized quantum batteries. Phys. Rev. B 2019, 100, 115142. [Google Scholar] [CrossRef]
  8. Caravelli, F.; Yan, B.; Garcia-Pintos, L.P.; Hamma, A. Energy storage and coherence in closed and open quantum batteries. Quantum 2021, 5, 505. [Google Scholar] [CrossRef]
  9. Mazzoncini, F.; Cavina, V.; Andolina, G.M.; Erdman, P.A.; Giovannetti, V. Optimal control methods for quantum batteries. Phys. Rev. A 2023, 107, 032218. [Google Scholar] [CrossRef]
  10. Alicki, R.; Fannes, M. Entanglement boost for extractable work from ensembles of quantum batteries. Phys. Rev. E 2013, 87, 042123. [Google Scholar] [CrossRef] [PubMed]
  11. Giorgi, G.L.; Campbell, S. Correlation approach to work extraction from finite quantum systems. J. Phys. B 2015, 48, 035501. [Google Scholar] [CrossRef]
  12. Andolina, G.M.; Keck, M.; Mari, A.; Campisi, M.; Giovannetti, V.; Polini, M. Extractable Work, the Role of Correlations, and Asymptotic Freedom in Quantum Batteries. Phys. Rev. Lett. 2019, 122, 047702. [Google Scholar] [CrossRef]
  13. Shi, H.L.; Ding, S.; Wan, Q.K.; Wang, X.H.; Yang, W.L. Entanglement, Coherence, and Extractable Work in Quantum Batteries. Phys. Rev. Lett. 2022, 129, 130602. [Google Scholar] [CrossRef] [PubMed]
  14. Dou, F.Q.; Lu, Y.Q.; Wang, Y.J.; Sun, J.A. Extended Dicke quantum battery with interatomic interactions and driving field. Phys. Rev. B 2022, 105, 115405. [Google Scholar] [CrossRef]
  15. Breuer, H.P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2007. [Google Scholar] [CrossRef]
  16. Weiss, U. Quantum Dissipative Systems; World Scientific: Singapore, 2012; Volume 13. [Google Scholar]
  17. De Vega, I.; Alonso, D. Dynamics of non-Markovian open quantum systems. Rev. Mod. Phys. 2017, 89, 015001. [Google Scholar] [CrossRef]
  18. Liu, J.; Segal, D.; Hanna, G. Loss-free excitonic quantum battery. J. Phys. Chem. C 2019, 123, 18303–18314. [Google Scholar] [CrossRef]
  19. Lidar, D.A. Review of Decoherence-Free Subspaces, Noiseless Subsystems, and Dynamical Decoupling. In Quantum Information and Computation for Chemistry; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2014; pp. 295–354. [Google Scholar] [CrossRef]
  20. Jamiołkowski, A.; Kamizawa, T.; Pastuszak, G. On invariant subspace in quantum control systems and some concepts of integrable quantum systems. Int. J. Theor. Phys. 2015, 54, 2662–2674. [Google Scholar] [CrossRef]
  21. Buča, B.; Prosen, T. A note on symmetry reductions of the Lindblad equation: Transport in constrained open spin chains. New J. Phys. 2012, 14, 073007. [Google Scholar] [CrossRef]
  22. Thingna, J.; Manzano, D.; Cao, J. Dynamical signatures of molecular symmetries in nonequilibrium quantum transport. Sci. Rep. 2016, 6, 1–11. [Google Scholar] [CrossRef]
  23. Albert, V.V.; Jiang, L. Symmetries and conserved quantities in Lindblad master equations. Phys. Rev. A 2014, 89, 022118. [Google Scholar] [CrossRef]
  24. Suter, D.; Álvarez, G.A. Colloquium: Protecting quantum information against environmental noise. Rev. Mod. Phys. 2016, 88, 041001. [Google Scholar] [CrossRef]
  25. Zanardi, P.; Rasetti, M. Noiseless quantum codes. Phys. Rev. Lett. 1997, 79, 3306. [Google Scholar] [CrossRef]
  26. Lidar, D.A.; Chuang, I.L.; Whaley, K.B. Decoherence-free subspaces for quantum computation. Phys. Rev. Lett. 1998, 81, 2594. [Google Scholar] [CrossRef]
  27. Lidar, D.A.; Birgitta Whaley, K. Decoherence-Free Subspaces and Subsystems. In Irreversible Quantum Dynamics; Benatti, F., Floreanini, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; pp. 83–120. [Google Scholar] [CrossRef]
  28. Beige, A.; Braun, D.; Tregenna, B.; Knight, P.L. Quantum computing using dissipation to remain in a decoherence-free subspace. Phys. Rev. Lett. 2000, 85, 1762. [Google Scholar] [CrossRef] [PubMed]
  29. Xue, P.; Xiao, Y.F. Universal quantum computation in decoherence-free subspace with neutral atoms. Phys. Rev. Lett. 2006, 97, 140501. [Google Scholar] [CrossRef]
  30. Manzano, D.; Hurtado, P.I. Symmetry and the thermodynamics of currents in open quantum systems. Phys. Rev. B 2014, 90, 125138. [Google Scholar] [CrossRef]
  31. Manzano, G.; Giorgi, G.L.; Fazio, R.; Zambrini, R. Boosting the performance of small autonomous refrigerators via common environmental effects. New J. Phys. 2019, 21, 123026. [Google Scholar] [CrossRef]
  32. Gherardini, S.; Campaioli, F.; Caruso, F.; Binder, F.C. Stabilizing open quantum batteries by sequential measurements. Phys. Rev. Res. 2020, 2, 013095. [Google Scholar] [CrossRef]
  33. Tejero, Á.; Thingna, J.; Manzano, D. Comment on “Loss-Free Excitonic Quantum Battery”. J. Phys. Chem. C 2021, 125, 7518–7520. [Google Scholar] [CrossRef]
  34. Manzano, D.; Hurtado, P. Harnessing symmetry to control quantum transport. Adv. Phys. 2018, 67, 1–67. [Google Scholar] [CrossRef]
  35. Liu, J.; Hanna, G. Efficient and deterministic propagation of mixed quantum-classical Liouville dynamics. J. Phys. Chem. Lett. 2018, 9, 3928–3933. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, Z.; Sergi, A.; Hanna, G. DECIDE: A Deterministic Mixed Quantum-Classical Dynamics Approach. Appl. Sci. 2022, 12, 7022. [Google Scholar] [CrossRef]
  37. Liu, J.; Hsieh, C.Y.; Segal, D.; Hanna, G. Heat transfer statistics in mixed quantum-classical systems. J. Chem. Phys. 2018, 149, 224104. [Google Scholar] [CrossRef]
  38. Liu, J.; Segal, D.; Hanna, G. Hybrid quantum-classical simulation of quantum speed limits in open quantum systems. J. Phys. A 2019, 52, 215301. [Google Scholar] [CrossRef]
  39. Carpio-Martínez, P.; Hanna, G. Nonequilibrium heat transport in a molecular junction: A mixed quantum-classical approach. J. Chem. Phys. 2019, 151, 074112. [Google Scholar] [CrossRef]
  40. Carpio-Martínez, P.; Hanna, G. Quantum bath effects on nonequilibrium heat transport in model molecular junctions. J. Chem. Phys. 2021, 154, 094108. [Google Scholar] [CrossRef]
  41. Wigner, E. On the Quantum Correction For Thermodynamic Equilibrium. Phys. Rev. 1932, 40, 749–759. [Google Scholar] [CrossRef]
  42. Imre, K.; Özizmir, E.; Rosenbaum, M.; Zweifel, P. Wigner method in quantum statistical mechanics. J. Math. Phys. 1967, 8, 1097–1108. [Google Scholar] [CrossRef]
  43. Wang, H.; Thoss, M.; Miller, W.H. Systematic convergence in the dynamical hybrid approach for complex systems: A numerically exact methodology. J. Chem. Phys. 2001, 115, 2979–2990. [Google Scholar] [CrossRef]
  44. Huo, P.; Coker, D.F. Semi-classical path integral non-adiabatic dynamics: A partial linearized classical mapping Hamiltonian approach. Mol. Phys. 2012, 110, 1035–1052. [Google Scholar] [CrossRef]
  45. Ishizaki, A.; Fleming, G.R. Theoretical examination of quantum coherence in a photosynthetic system at physiological temperature. Proc. Natl. Acad. Sci. USA 2009, 106, 17255–17260. [Google Scholar] [CrossRef] [PubMed]
  46. Shanks, E.B. Higher Order Approximations of Runge-Kutta Type; Technical Report NASA-TN-D-2920; NASA Marshall Space Flight Center NTRS Document ID: 19650022581; NTRS Research Center: Moffett Field, CA, USA, 1965. [Google Scholar]
Figure 1. Time-dependent site populations (upper panels) and energies (lower panels) of the OQN, SBP, and baths in the storage (left panels) and discharge (right panels) phases. In panel (d), energy changes are calculated by subtracting the initial value of the energy from the value at each time. The results were generated using the following parameter set: E 1 = 250 cm 1 , E i { 2 , 3 , 5 , 6 } = 200 cm 1 , E 4 = 0 cm 1 , h = 60 cm 1 , T L = T R = 300 K, T p = 300 K, λ b = 35 cm 1 , λ p = 10 cm 1 , ω c = ω p = 106 cm 1 , ω m a x = 50 ω c , and M = 100 .
Figure 1. Time-dependent site populations (upper panels) and energies (lower panels) of the OQN, SBP, and baths in the storage (left panels) and discharge (right panels) phases. In panel (d), energy changes are calculated by subtracting the initial value of the energy from the value at each time. The results were generated using the following parameter set: E 1 = 250 cm 1 , E i { 2 , 3 , 5 , 6 } = 200 cm 1 , E 4 = 0 cm 1 , h = 60 cm 1 , T L = T R = 300 K, T p = 300 K, λ b = 35 cm 1 , λ p = 10 cm 1 , ω c = ω p = 106 cm 1 , ω m a x = 50 ω c , and M = 100 .
Molecules 29 00889 g001
Figure 2. Site populations (left panels) and energy changes (right panels) in the discharge ( χ = 1 ) phase for different bath temperature gradients. (a,b) Time-dependent site populations and energy changes for right bath temperatures of 300 K (solid lines) and 900 K (dashed lines), with T L = 300 K. (c,d) Site populations and energy changes after 1 ps for different right bath temperatures and a left bath temperature of 300 K (denoted by right bath gradient), and different left-bath temperatures and a right-bath temperature of 300 K (denoted by left bath gradient).
Figure 2. Site populations (left panels) and energy changes (right panels) in the discharge ( χ = 1 ) phase for different bath temperature gradients. (a,b) Time-dependent site populations and energy changes for right bath temperatures of 300 K (solid lines) and 900 K (dashed lines), with T L = 300 K. (c,d) Site populations and energy changes after 1 ps for different right bath temperatures and a left bath temperature of 300 K (denoted by right bath gradient), and different left-bath temperatures and a right-bath temperature of 300 K (denoted by left bath gradient).
Molecules 29 00889 g002
Figure 3. (a) Site populations and (b) energy changes after 1 ps in the discharge phase for different right-bath temperatures and site 4 energies, with E 1 = 250 cm 1 , E i { 2 , 3 , 5 , 6 } = 200 cm 1 , and T L = 300 K.
Figure 3. (a) Site populations and (b) energy changes after 1 ps in the discharge phase for different right-bath temperatures and site 4 energies, with E 1 = 250 cm 1 , E i { 2 , 3 , 5 , 6 } = 200 cm 1 , and T L = 300 K.
Molecules 29 00889 g003
Figure 4. (a) Site populations and (b) energy changes in the discharge phase after 1 ps for different right-bath reorganization energies and a left-bath reorganization energy of 35 cm 1 (denoted by right-bath variation), and different left-bath reorganization energies and a right-bath reorganization energy of 35 cm 1 (denoted by left-bath variation).
Figure 4. (a) Site populations and (b) energy changes in the discharge phase after 1 ps for different right-bath reorganization energies and a left-bath reorganization energy of 35 cm 1 (denoted by right-bath variation), and different left-bath reorganization energies and a right-bath reorganization energy of 35 cm 1 (denoted by left-bath variation).
Molecules 29 00889 g004
Figure 5. Time-dependent site populations (left panels) and energy changes (right panels) in the discharge phase ( χ = 1 ) for different SS energies. (a,b) Results for E 1 = E 4 = 100 cm−1, E 1 = E 4 = 200 cm−1, and E 1 = E 4 = 250 cm−1. (c,d) Results for E 1 = 250 cm−1, and E 4 = 0 , 100, and 250 cm−1. The values of the remaining parameters are E i { 2 , 3 , 5 , 6 } = 200 cm 1 , T A = T B = 300 K, and λ b = 35 cm 1 .
Figure 5. Time-dependent site populations (left panels) and energy changes (right panels) in the discharge phase ( χ = 1 ) for different SS energies. (a,b) Results for E 1 = E 4 = 100 cm−1, E 1 = E 4 = 200 cm−1, and E 1 = E 4 = 250 cm−1. (c,d) Results for E 1 = 250 cm−1, and E 4 = 0 , 100, and 250 cm−1. The values of the remaining parameters are E i { 2 , 3 , 5 , 6 } = 200 cm 1 , T A = T B = 300 K, and λ b = 35 cm 1 .
Molecules 29 00889 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Z.; Hanna, G. Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery. Molecules 2024, 29, 889. https://doi.org/10.3390/molecules29040889

AMA Style

Liu Z, Hanna G. Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery. Molecules. 2024; 29(4):889. https://doi.org/10.3390/molecules29040889

Chicago/Turabian Style

Liu, Zhe, and Gabriel Hanna. 2024. "Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery" Molecules 29, no. 4: 889. https://doi.org/10.3390/molecules29040889

APA Style

Liu, Z., & Hanna, G. (2024). Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery. Molecules, 29(4), 889. https://doi.org/10.3390/molecules29040889

Article Metrics

Back to TopTop