Transitionless quantum driving based wireless power transfer

TL;DR

This paper introduces a transitionless quantum driving method for wireless power transfer that is faster, robust against parameter variations, and resilient to coil losses, improving efficiency over traditional adiabatic approaches.

Contribution

It presents a novel application of transitionless quantum driving to wireless power transfer, enhancing speed and robustness compared to existing methods.

Findings

Transfer efficiency is significantly improved with TQD.

The scheme is robust against coupling variations.

It shows resilience to intrinsic coil losses.

Abstract

Shortcut to adiabaticity (STA) techniques have the potential to drive a system beyond the adiabatic limits. Here, we present a robust and efficient method for wireless power transfer (WPT) between two coils based on the so-called transitionless quantum driving (TQD) algorithm. We show that it is possible to transfer power between the coils significantly fast compared to its adiabatic counterpart. The scheme is fairly robust against the variations in the coupling strength and the coupling distance between the coils. Also, the scheme is found to be reasonably immune to intrinsic losses in the coils.

Equations

$$H(t)=\sum _n|n(t)⟩E_n(t)⟨n(t)|$$

$$|{\Psi }_n(t)⟩=e^{i{\xi }_n(t)}|n(t)⟩,$$

$$H^{ad}(t)=U^{-1}H(t)U-i\hslash U^{-1}\dot{U}$$

$$H_1(t)=i\hslash \sum _m[|{\partial }_{t}n(t)⟩⟨n(t)|-⟨n(t)⟩{\partial }_{t}n(t)|n(t)⟩⟨n(t)|]$$

$$v(t)=L\frac{di(t)}{dt},$$

$$i(t)=-C\frac{dv(t)}{dt}$$

$$\frac{d^{2}v(t)}{d{t}^2}+{\omega }_0^{2}v(t)=0,$$

$$\frac{d{a}_{\pm }(t)}{dt}=\pm j{\omega }_0{a}_{\pm }(t),$$

$$a_{\pm }(t)=\sqrt{\frac{C}{2}}(v(t)\pm j\sqrt{\frac{L}{C}}i(t))$$

$$v(t)=|V|\cos ({\omega }_{0}t),$$

$$i(t)=\sqrt{\frac{C}{L}}|V|\sin ({\omega }_{0}t)$$

$$\frac{da(t)}{dt}=j{\omega }_{0}a(t)-\Gamma a(t),$$

$$\frac{d{a}_s(t)}{dt}=(j{\omega }_s-{\Gamma }_s)a_s(t)+j\kappa a_d(t)$$

$$\frac{d{a}_d(t)}{dt}=(j{\omega }_d-{\Gamma }_d-{\Gamma }_w)a_s(t)+j\kappa a_s(t)$$

$$H(t)=\left(\begin{array}{cc}-{\omega }_s-j{\Gamma }_s& -\kappa (t)\\ -\kappa (t)& -{\omega }_d-j{\Gamma }_d-j{\Gamma }_w\end{array}\right)$$

$$H(t)=\left(\begin{array}{cc}\frac{\Delta (t)}{2}-j{\Gamma }_s& -\kappa (t)\\ -\kappa (t)& -\frac{\Delta (t)}{2}-j{\Gamma }_d-j{\Gamma }_w\end{array}\right)$$

$$B_{+}(t)$$

$$B_{-}(t)$$

$$|⟨{\partial }_t{B}_{\pm }(t)|B_{\mp }(t)⟩|\ll \sqrt{(4\kappa (t{)}^2+\Delta (t{)}^2)}$$

$$\kappa (t)={\kappa }_0;\Delta (t)=\delta +\beta (t-t_0)$$

$$\Gamma =\left(\begin{array}{cc}{\Gamma }_s& 0\\ 0& {\Gamma }_d+{\Gamma }_w\end{array}\right)$$

$$[B_{+}(t),B_{-}(t){]}^T=U(\Theta (t){)}^{†}[b_s(t),b_d(t){]}^T$$

$$H_1(t)=j\left(\begin{array}{cc}0& {\kappa }_a(t)\\ -{\kappa }_a(t)& 0\end{array}\right)$$

$$H(t)+H_{CD}(t)=\left(\begin{array}{cc}(\Delta -\dot{\varphi })/2& \sqrt{{\kappa }^2(t)+{\kappa }_a^2(t)}\\ \sqrt{{\kappa }^2(t)+{\kappa }_a^2(t)}& -(\Delta -\dot{\varphi }/2)\end{array}\right),$$

$${\kappa }_a(t)=\frac{|\dot{\Delta }(t)\kappa (t)-\dot{\kappa }(t)\Delta (t)|}{{\Delta }^2(t)+4{\kappa }^2(t)}$$

$$\dot{\varphi }(t)=\frac{2\Delta (t){\dot{\Delta }}^2(t)}{{\Delta }^2(t)+4{\kappa }^2(t)+{\dot{\Delta }}^2(t)}$$

$$\left(\begin{array}{c}{b}_s^{\prime }\\ b_d^{\prime }\end{array}\right)=\left(\begin{array}{cc}{e}^{-j\varphi /2}& 0\\ 0& e^{j\varphi /2}\end{array}\right)\left(\begin{array}{c}{b}_s\\ b_d\end{array}\right)$$

$$\frac{d\rho }{dt}=-j[H,\rho ]-\frac{1}{2}\{\Gamma ,\rho \}$$

$$\eta =\frac{{\Gamma }_w{\int }_0^T|b_d(t){|}^{2}dt}{{\Gamma }_s{\int }_0^T|b_s(t){|}^{2}dt+({\Gamma }_d+{\Gamma }_w){\int }_0^T|b_d(t){|}^{2}dt}$$

Full text

I Introduction

Techniques for coherent manipulation of quantum states, specially for two-level and three-level quantum systems, are of great interest in modern day quantum optics and atomic physics. In this context, resonant excitation via pulsed radiation and adiabatic evolution by means of frequency sweep, with controlled chirp across resonance, are particularly well known. These two methods have been investigated theoretically and demonstrated experimentally by a wide range of implementations in past few decades 1 . To reach at a desired target state, from a given one, using resonant excitation, a constant frequency $\pi$ pulse is required. One needs to have precise control over the pulse area which makes it sensitive to the fluctuations in various parameters. The adiabatic processes like rapid adiabatic passage (RAP), Stark chirped rapid adiabatic passage (SCRAP), stimulated Raman adiabatic passage (STIRAP) etc. do not require such precise control and are robust against parameter fluctuations 2 . However, for adiabatic evolution, the system must evolve adiabatically with time making those processes relatively slow. Furthermore, no matter how robust the process is, it is necessary to drive a system as fast as possible in order to reduce the effects of definite decoherences present in the system. Efforts have been made to address those issues in recent years. Current developments in this direction are mainly based on the idea of accelerating the evolution beyond the adiabatic regime without violating the principles of state transfer. Few shortcut to adiabaticity (STA) techniques, such as transitionless quantum driving (TQD) 3 , counter diabatic algorithm (CDA) 4 and Lewis-Riesenfeld invariant (LRI) approach 5 are put forward to speed up the time evolution of quantum systems. In past few years, these methods have been explored rigorously across various branches of physics such as waveguide couplers 6 , Bose-Einstein condensates 7 , entangled state preparation 8 , non-hermitian systems 9 and so on. Owing to the ubiquitous nature of adiabatic processes, STA techniques span a broad range of applications. One possible potential use of these techniques may be in the so-called wireless power transfer technology.

In modern age, various wireless technologies play crucial role in our day-to-day life. Since the early days of electromagnetics, significant progress has been made in transferring information via wireless communication. In contrast, wireless power transfer (WPT) gained little progress in the last century. However, recent tide in the use of electronic appliances and requirement for short and mid-range wireless energy exchange has helped WPT getting tremendous attention, and studies on WPT systems has gained momentum in past few years 10 . Modern day WPT techniques are mainly based on the mutual induction between two coils. Two fundamental principles, which are commonly used in such systems, are near-field non-radiative magneto-inductive effects and the resonant coupling between both the emitter and the receiver coils 11 ; 12 ; 13 ; 14 . Like resonant excitation (as discussed above), here also it is extremely important to maintain the resonance in both the coils. Otherwise it may result in decrease in the efficiency11 . Few studies has been put forward to solve this issue 15 ; 16 . WPT is also dependent on the coupling distance between the coils and it has been shown in few articles that in the strong coupling regime, it is possible to enhance efficiency for larger distances 14 .

In this work, we study WPT in a system of two inductively coupled coils by exploiting the adiabatic and the TQD techniques. Using coupled mode theory 17 , we find out the governing equations (which are similar to the Schrödinger equation for two level system) and devise power transfer mechanism which is impervious to the coupling strength and distance between the coils and intrinsic losses present in the system.

The remainder of this article is organized as follows. In Sec.II, at first we present a brief discussion on the TQD algorithm and a detailed study of the coupled LC circuits using the coupled mode theory,followed by application of adiabatic and TQD techniques in order to achieve wireless power transfer. Subsequently, in Sec.III, we have discussed numerical calculations and relevant results followed by the conclusion of our work in Sec. IV.

II Theoretical formulation

II.1 Transitionless quantum driving

Shortcut techniques for adiabatic processes ensued from the limitations of the adiabatic theorem itself 3 . Adiabatic theorem states that if a state $\ket{n(t=0)}$ is an eigenstate of a Hamiltonian $H(t=0)$ and the Hamiltonian evolves very slowly (adiabatically) with time, then $\ket{n(t)}$ will remain in its eigenstate throughout the evolution. However it is possible to inverse engineer the system Hamiltonian from a set of predefined eigenstates to drive the system beyond the adiabatic limit using TQD 3 . According to this method, the transition probability between the instantaneous eigenstates could be made zero to achieve exact evolution of the eigenstates. To understand it, let us consider the Hamiltonian:

[TABLE]

The instantaneous eigenstates of this Hamiltonian are the adiabatic states, given by:

[TABLE]

where $\xi_{n}(t)$ is the adiabatic phase. We can define a unitary transformation to diagonalize $H(t=0)$ in order to obtain the adiabatic Hamiltonian, as follows:

[TABLE]

Here $U$ is a time dependent unitary operator. The first term on the RHS of Eq. (3), has instantaneous eigenvalues and is diagonal itself. So it can drive the system along the adiabatic states $\ket{\Psi_{n}(t)}$. The second term refers to the non-adiabatic correction, and generally off-diagonal. When $H(t)$ changes slowly with time, the off-diagonal part of $H^{ad}(t)$ becomes negligible. Hence the probability for transfer among the instantaneous eigenstates tends to be zero but the evolution becomes slow. Problem arises when $H(t)$ changes fast with time, resulting in the breakdown of the adiabatic condition, as the off-diagonal part of $H^{ad}(t)$ becomes stronger and the system no longer follows the adiabatic path. In such a scenario, the TQD algorithm could be used to negate the effects of the non-adiabatic terms. According to this technique, we need to add an additional coupling Hamiltonian to cancel out the off-diagonal terms from $H^{ad}(t)$, so that it becomes diagonal once and for all and the system can be driven exactly regardless of the speed of the time evolution. In other words, one can drive the system infinitely fast along the adiabatic path. The additional Hamiltonian in the $\ket{n(t)}$ basis, is given by:

[TABLE]

The total effective Hamiltonian, $H_{eff}(t)=H(t)+H_{1}(t)$, drives the system exactly beyond the adiabatic limit. One needs to configure $H_{1}(t)$ appropriately, depending on the system, to drive it in infinitely short time.

II.2 Coupled mode Theory for two interacting coils

Let us consider a loss-less $LC$ circuit, without loss, with current $i(t)$ flowing in it and voltage, $v(t)$, across $L$ and $C$. The equations describing the voltage and the current can be written in terms of the following coupled differential equations:

[TABLE]

Above equations could be rewritten, as a second order differential equation for voltage, as follows:

[TABLE]

where $\omega_{0}^{2}=1/LC$. The coupled equations in Eq. (5) can be expressed by two uncoupled equations for mode amplitudes $a_{+}(t)$ and $a_{-}(t)$:

[TABLE]

Here, $j=\sqrt{-1}$. The solutions for the current and the voltage from Eq. (5) and Eq. (6), when subjected to proper boundary conditions, are given by:

[TABLE]

Here $|V|$ is the peak voltage. Using above solutions, the total energy of the system could be obtained as: $|a_{\pm}|^{2}=\frac{C}{2}|V|^{2}=W$, where $a_{+}$ being the positive frequency component of the mode amplitudes, while $a_{-}$ is its negative counterpart. In the rest of the analysis,we will consider the positive frequency component only and will drop the ‘+’ subscript for simplicity. Taking loss in the system into account, the equation is written in the following modified form:

[TABLE]

where $\Gamma$ is the decay rate due to the dissipation from the coils. In our study we consider two such coils, namely the Source and the Drain (Fig. (1)), which are coupled by mutual inductance between them. These two coils are off resonant to each other having different resonant frequencies $\omega_{S}$ and $\omega_{D}$. The coupling between the coils is given by $\kappa(t)=M(t)\sqrt{\omega_{s}\omega_{d}/L_{s}L_{d}}$, where $M(t)$ is the mutual inductance between these two coils. The mode amplitudes are $a_{s}(t)$ and $a_{d}(t)$ respectively and are coupled to each other. Total energy of the system is $|a_{s}(t)|^{2}+|a_{d}(t)|^{2}$. This system can be expressed using the following set of coupled equations:

[TABLE]

Here $\Gamma_{s}$ and $\Gamma_{d}$ are the intrinsic loss rates of the source and the drain respectively. $\Gamma_{w}$ is the work extracted from the drain coil.

II.3 Energy transfer protocols

In order to design the energy transfer protocols, we discuss two methods here: the adiabatic passage and transitionless quantum driving. From the coupled mode theory, discussed above, we can characterise the system by defining the Hamiltonian in $[a_{s}(t)\quad a_{d}(t)]^{T}$ basis, as follows:

[TABLE]

As the circuits,chosen here, are not resonant to each other, we consider the system in a rotating frame of reference, where the interaction Hamiltonian in the diabatic basis is as follows:

[TABLE]

The diabatic basis are given by, $b_{s,d}(t)=a_{s,d}\exp[-j(\omega_{s}+\omega_{d})t/2]$.In this new basis, one have to consider only the frequency difference between the coils, given by, $\Delta(t)=\omega_{d}(t)-\omega_{s}(t)$.

II.3.1 Adiabatic following

The instantaneous eigenvectors of Eq. (14) without losses are,

[TABLE]

Here $\Theta(t)$ is the angle of mixing,given by $\Theta(t)=\frac{1}{2}\tan^{-1}(2\kappa/\Delta)$. $B_{\pm}(t)$ are also known as adiabatic states. For adiabatic evolution, one needs to vary the Hamiltonian infinitely slowly or adiabatically so that the system always follows a particular state $B_{+}$ or $B_{-}$ during a complete cycle of the time evolution. To achieve this, the evolution ought to follow the adiabatic condition, which is obtained by comparing the non-adiabatic correction terms (Eq. (3)) to the instantaneous eigen-energies. It is expressed as follows:

[TABLE]

The fulfillment of this condition makes the transition probability between $B_{+}(t)$ and $B_{-}(t)$ zero. Moreover, if we assume that the system is initially in the state $B_{+}(t)$ and the power is in the source coil at $t=0$, which refers to $\Theta(0)=0$ , then, by rotating $\Theta$ clockwise to $\frac{\pi}{2}$ , one could arrive at the final state, $b_{d}(t)$. Thus the power ends up in the drain coil. This rotation could be achieved by sweeping $\Delta(t)$ from a large negative value to a large positive value. For our system, we chose $\kappa(t)$ and $\Delta(t)$ according to the well known Landau-Zener scheme 18 , given by:

[TABLE]

where $\delta$ is some arbitrary offset between the frequencies of the two coils and $\beta$ determines the slope in $\Delta(t)$, which eventually controls the speed of the process. Under such choices, large $t_{0}$ is needed to satisfy the adiabatic condition and $t_{0}$ ,effectively, determines the width of the evolution cycle. When the loss rates are non-zero in both the source and the drain coil, the system could be considered as dissipative. These dissipations can be modelled mathematically via the dissipation matrix:

[TABLE]

It should be noted that, for adiabatic evolution, the intrinsic loss rates $\Gamma_{s}$ and $\Gamma_{d}$ should be less than the coupling strength $\kappa_{0}$, otherwise the evolution would not be possible and the power would be lost from the coil itself.

II.3.2 Shortcut to adiabaticity

For shortcut, we first transform our Hamiltonian in Eq. (14), in the adiabatic basis, using Eq. (3),where the basis states are related as:

[TABLE]

The non-adiabatic correction terms are generally negligible under adiabatic approximation. The Landau-Zener formula for the total transition probability between adiabatic states in our case is, $p=\exp[-2\pi\kappa_{0}^{2}/\beta]$. One can obtain $p\simeq 0$ when Eq. (16) is satisfied, which can be written simply as $\beta\ll 8\kappa_{0}^{2}$. However if one wants to drive the evolution faster, non-adiabatic corrections becomes stronger and $p$ no longer remains zero. To avoid such a scenario, the adiabatic Hamiltonian needs to be diagonalized exactly. To serve this cause, we add the additional interaction to the adiabatic Hamiltonian as proposed by Berry 3 .We use Eq. (4) to find the additional interaction. The second term in Eq. (4) vanishes owing to the orthogonality of $B_{\pm}(t)$. The first term can be written in the adiabatic basis as follows:

[TABLE]

The total Hamiltonian, required for the transitionless driving, is given by

[TABLE]

where

[TABLE]

and

[TABLE]

Here $(\Delta-\dot{\phi})/2$ is the effective detuning and $\sqrt{\kappa^{2}(t)+\kappa_{a}^{2}(t)}=\kappa_{eff}$ is the modified coupling between the two coils. $\phi$ characterizes another unitary rotation given by

[TABLE]

With all these unitary transformations, one needs to keep track of all the bases used in the process and keep them consistent. The mixing angle $\Theta$ should be adjusted properly via the boundary conditions, given by $\Theta(0)=0$, $\Theta(T)=\pi/2$ and $\kappa_{a}(0)=\kappa_{a}(T)=0$.

III results and discussion

To envisage the transfer of energy from the source to the drain ,with the effects of intrinsic losses taken into account, we followed the standard density matrix approach and therefore solved the master equation, given by

[TABLE]

Here $\rho(t)$ is the density matrix and $\Gamma$ represents the dissipation matrix as described in Eq. (18). The diagonal elements of $\rho(t)$ are $\rho_{ss}=|b_{s}|^{2}$ and $\rho_{dd}=|b_{d}|^{2}$ which represents the energy of the source and

the drain coil respectively. For adiabatic evolution we chose the interaction Hamiltonian in Eq. (14) and for TQD, $H(t)$ is taken from Eq. (21).

In Fig. (2), we present the results showing the evolution of fractional energies, $|b_{s,d}(t)|^{2}/|b_{s}(0)|^{2}$ , in presence of intrinsic losses, using both the adiabatic and the TQD algorithm. Even when the adiabatic condition is being satisfied, as shown in Fig. (2a) , as a result of adiabatic evolution, the fractional energy attains the value, on the order of $0.25$ or so at the end of the given time window. The requirement of large time, i.e. $t_{0}=10^{-4}s$, so that the adiabaticity condition is maintained, results in energy dissipation due to the intrinsic losses.

However, when TQD is applied, fractional energy of the drain coil almost attains nearly the same value for the same period of time (Fig. (2b)). But when the adiabatic condition is violated i.e. $\beta\geq 8\kappa_{0}^{2}$, enhancement of energy in the drain coil could be observed in Fig. (2c). The affects of intrinsic losses are almost eliminated as the time period becomes shorter and shorter(Fig. (2d)). It is worthwhile to note the effect of $\beta$ on the power transfer mechanism. Physically, $\beta$ determines the slope of $\Delta$ and thereby it controls the time period required to complete a single power transfer cycle. In the adiabatic regime, $\beta$ is relatively small and hence the required period $T_{adiabatic}$ is larger and because of that, power is dissipated from the source coil during the process. When $\beta$ is large, the frequency sweep becomes faster and $\omega_{s}(t)$ becomes steeper (here $\omega_{d}$ is taken as constant) as shown in Fig. (3) which results in squeezing of the period. Thus the power is transferred to the drain in short time with minimum loss from the source. Also it is obvious that one needs to repeat the cycle over and over again to transfer power for practical purposes.

The work efficiency of our system is defined as the ratio between the useful extracted power from the drain, $P_{work}=\Gamma_{w}\int_{0}^{T}|b_{d}(t)|^{2}dt$ to the total time averaged power $P_{total}$ in the system over a particular time period $T$, given by

[TABLE]

We studied efficiency of the system against the variation of the initial frequency difference between the coils for different values of $\kappa_{0}/\Gamma_{s,d}$. The efficiency strongly depends on the coupling strength $\kappa_{0}$ for both the adiabatic and the TQD approach and it increases with increasing $\kappa_{0}$ as depicted in Fig. (4a) . Although this is not that surprising, but in the beyond adiabatic regime (for shorter periods), efficiency for adiabatic method decreases rapidly. However efficiency remains intact for TQD algorithm which can be seen from Fig. 4b.

The coupling between the coils for WPT systems are generally very sensitive to the distance, $d$ , between the coils and it decays very rapidly with larger coupling distances 11 ; 12 . In Fig. 5(a) and 5(b) we depict the dependence of coupling strength $\kappa(d)$ and efficiency $\eta(d)$ respectively on the distance $d$ between the coils. From Fig. (5b) it is clear that, even in the adiabatic regime, the efficiency for the adiabatic case decreases with increasing $d$. $\eta$ tends to zero for $d=2m$ or so as the strength of $\kappa$ goes down. But in the case of TQD, $\eta$ maintains a steady value for large $d$ and that certainly gets enhanced with larger $\beta$ values. The reason behind such behavior could be understood from the formalism of inverse engineering. We can observe from Fig. (5a) that as $\kappa$ decreases with distance, we require an additional coupling $\kappa_{a}$ which increases with distance so that the effective coupling $\kappa_{eff}$ constitutes a reasonable strength for sustained power transfer over a certain range of distance.

We have also studied the efficiency of our scheme against the variations in $\kappa_{0}$ and the intrinsic losses ($\Gamma_{s}=\Gamma_{d}$). From the contour plots in Fig. (6), we observe that $\eta$ is nearly unity for lower losses and higher $\kappa_{0}$ in the adiabatic regime. As we move from Fig. (6a) to Fig. (6c), adiabaticity gradually breaks down and efficiency also decreases gradually. Finally it goes to zero when $t_{0}=10^{-6}$s for any reasonable amount of losses. On the other hand, using TQD, we find that the achievable efficiency is highly robust against variations in $\kappa_{0}$ and $\Gamma_{s,d}$ unlike its adiabatic counterpart. Also the efficiency is found to get enhanced with decreasing time window.

IV Conclusion

In conclusion, we have explored a WPT system in the light of adiabatic method and transitionless quantum driving method. Our findings could be summarized as follows. The adiabatic evolution of power is a useful way to transfer power between two coils. Unlike resonant WPT systems, it uses two off-resonant coils and frequency sweeping is used for power transfer. However, it has to fulfil the adiabatic condition which makes the time required to transfer power in each cycle longer, resulting in dissipation of power from the source coil and hence efficiency decreases. On the other hand, the TQD algorithm enables us to enhance the efficiency of power transfer. The algorithm suggest that, it is possible to decrease the transfer time and increase the efficiency by invoking an additional interaction between the coils.The amount of power dissipated from the source also gets decreased. The WPT method using TQD shows more robustness compared to the adiabatic one against the variations in the parameters like the coupling strength, intrinsic losses and the coupling distance.

ACKNOWLEDGMENTS

K. P. would like to gratefully acknowledge the research fellowship provided by MHRD, Govt. of India.