Efficient near-field wireless energy transfer using adiabatic system variations

Description

BACKGROUND

The disclosure relates to wireless energy transfer. Wireless energy transfer can for example, be useful in such applications as providing power to autonomous electrical or electronic devices.

Radiative modes of omni-directional antennas (which work very well for information transfer) are not suitable for such energy transfer, because a vast majority of energy is wasted into free space. Directed radiation modes, using lasers or highly-directional antennas, can be efficiently used for energy transfer, even for long distances (transfer distance L_TRANS>>L_DEV, where L_DEVis the characteristic size of the device and/or the source), but may require existence of an uninterruptible line-of-sight and a complicated tracking system in the case of mobile objects. Some transfer schemes rely on induction, but are typically restricted to very close-range (L_TRANS<<L_DEV) or low power (˜mW) energy transfers.

The rapid development of autonomous electronics of recent years (e.g. laptops, cell-phones, house-hold robots, that all typically rely on chemical energy storage) has led to an increased need for wireless energy transfer.

SUMMARY

Disclosed is a method for transferring energy wirelessly. The method includes i) transferring energy wirelessly from a first resonator structure to an intermediate resonator structure, wherein the coupling rate between the first resonator structure and the intermediate resonator structure is κ_1B; ii) transferring energy wirelessly from the intermediate resonator structure to a second resonator structure, wherein the coupling rate between the intermediate resonator structure and the second resonator structure is κ_B2; and iii) during the wireless energy transfers, adjusting at least one of the coupling rates κ_1Band κ_B2to reduce energy accumulation in the intermediate resonator structure and improve wireless energy transfer from the first resonator structure to the second resonator structure through the intermediate resonator structure.

Embodiments of the method may include one or more of the following features.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to minimize energy accumulation in the intermediate resonator structure and cause wireless energy transfer from the first resonator structure to the second resonator structure.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to maintain energy distribution in the field of the three-resonator system in an eigenstate having substantially no energy in the intermediate resonator structure. For example, the adjustment of at least one of the coupling rates κ_1Band κ_B2can further cause the eigenstate to evolve substantially adiabatically from an initial state with substantially all energy in the resonator structures in the first resonator structure to a final state with substantially all of the energy in the resonator structures in the second resonator structure.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to include adjustments of both coupling rates κ_1Band κ_B2during wireless energy transfer.

The resonator structures can each have a quality factor larger than 10.

The first and second resonator structures can each have a quality factor greater than 50.

The first and second resonator structures can each have a quality factor greater than 100.

The resonant energy in each of the resonator structures can include electromagnetic fields. For example, the maximum value of the coupling rate κ_1Band the maximum value of the coupling rate κ_B2for inductive coupling between the intermediate resonator structure and each of the first and second resonator structures can each be larger than twice the loss rate Γ for each of the first and second resonators. Moreover, The maximum value of the coupling rate κ_1Band the maximum value of the coupling rate κ_B2for inductive coupling between the intermediate resonator structure and each of the first and second resonator structures can each be larger than four (4) times the loss rate Γ for each of the first and second resonators.

Each resonator structure can have a resonant frequency between 50 KHz and 500 MHz.

The maximum value of the coupling rate κ_1Band the maximum value of the coupling rate κ_B2can each be at least five (5) times greater than the coupling rate between the first resonator structure and the second resonator structure.

The intermediate resonator structure can have a rate of radiative energy loss that is at least twenty (20) times greater than that for either the first resonator structure or the second resonator structure.

The first and second resonator structures can be substantially identical.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to cause peak energy accumulation in the intermediate resonator structure to be less than five percent (5%) of the peak total energy in the three resonator structures.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to cause peak energy accumulation in the intermediate resonator structure during the wireless energy transfers to be less than ten percent (10%) of the peak total energy in the three resonator structures.

Adjusting at least one of the coupling rates κ_1Band κ_B2can include adjusting a relative position and/or orientation between one or more pairs of the resonator structures. Moreover, adjusting at least one of the coupling rates κ_1Band κ_B2can include adjusting a resonator property of one or more of the resonator structures, such as mutual inductance.

The resonator structures can include a capacitively loaded loop or coil of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

The resonator structures can include an inductively loaded rod of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

The wireless energy transfers are non-radiative energy transfers mediated by a coupling of a resonant field evanescent tail of the first resonator structure and a resonant field evanescent tail of the intermediate resonator structure and a coupling of the resonant field evanescent tail of the intermediate resonator structure and a resonant field evanescent tail of the second resonator structure.

The adjustment of the at least one of the coupling rates can define a first mode of operation, wherein the reduction in the energy accumulation in the intermediate resonator structure is relative to energy accumulation in the intermediate resonator structure for a second mode of operation of wireless energy transfer among the three resonator structures having a coupling rate κ′_1Bfor wireless energy transfer from the first resonator structure to the intermediate resonator structure and a coupling rate κ′_B2for wireless energy transfer from the intermediate resonator structure to the second resonator structure with κ′_1Band κ′_B2each being substantially constant during the second mode of wireless energy transfer, and wherein the adjustment of the coupling rates κ_1Band κ_2Bin the first mode of operation can be selected to κ_1B,κ_B2<√{square root over ((κ′_1B²+κ′_B2²)/2)}. Moreover, the first mode of operation can have a greater efficiency of energy transferred from the first resonator to the second resonator compared to that for the second mode of operation. Further, the first and second resonator structures can be substantially identical and each one can have a loss rate Γ_A, the intermediate resonator structure can have a loss rate Γ_B, and wherein Γ_B/Γ_Acan be greater than 50.

Also, a ratio of energy lost to radiation and total energy wirelessly transferred between the first and second resonator structures in the first mode of operation is less than that for the second mode of operation. Moreover, the first and second resonator structures can be substantially identical and each one can have a loss rate Γ_Aand a loss rate only due to radiation Γ_A,rad, the intermediate resonator structure can have a loss rate Γ_Band a loss rate only due to radiation Γ_B,rad, and wherein Γ_B,rad/Γ_B>Γ_A,rad/Γ_A.

The first mode of operation the intermediate resonator structure interacts less with extraneous objects than it does in the second mode of operation.

During the wireless energy transfer from the first resonator structure to the second resonator structure at least one of the coupling rates can be adjusted so that κ_1B<<κ_B2at a start of the energy transfer and κ_1B>>κ_B2by a time a substantial portion of the energy has been transferred from the first resonator structure to the second resonator structure.

The coupling rate κ_B2can be maintained at a fixed value and the coupling rate κ_1Bis increased during the wireless energy transfer from the first resonator structure to second resonator structure.

The coupling rate κ_1Bcan be maintained at a fixed value and the coupling rate κ_B2is decreased during the wireless energy transfer from the first resonator structure to second resonator structure.

During the wireless energy transfer from the first resonator structure to second resonator structure, the coupling rate κ_1Bcan be increased and the coupling rate κ_B2is decreased.

The method may further include features corresponding to those listed for one or more of the apparatuses and methods described below.

In another aspect, disclosed is an apparatus including: first, intermediate, and second resonator structures, wherein a coupling rate between the first resonator structure and the intermediate resonator structure is κ_1Band a coupling rate between the intermediate resonator structure and the second resonator structure is κ_B2; and means for adjusting at least one of the coupling rates κ_1Band κ_B2during wireless energy transfers among the resonator structures to reduce energy accumulation in the intermediate resonator structure and improve wireless energy transfer from the first resonator structure to the second resonator structure through the intermediate resonator structure.

Embodiments for the apparatus can include one or more of the following features.

The means for adjusting can include a rotation stage for adjusting the relative orientation of the intermediate resonator structure with respect to the first and second resonator structures.

The means for adjusting can include a translation stage for moving the first and/or second resonator structures relative to the intermediate resonator structure.

The means for adjusting can include a mechanical, electro-mechanical, or electrical staging system for dynamically adjusting the effective size of one or more of the resonator structures.

The apparatus may further include features corresponding to those listed for the method described above, and one or more of the apparatuses and methods described below.

In another aspect, a method for transferring energy wirelessly includes i): transferring energy wirelessly from a first resonator structure to a intermediate resonator structure, wherein the coupling rate between the first resonator structure and the intermediate resonator structure is κ_1B; ii) transferring energy wirelessly from the intermediate resonator structure to a second resonator, wherein the coupling rate between the intermediate resonator structure and the second resonator structure is κ_B2; and iii) during the wireless energy transfers, adjusting at least one of the coupling rates κ_1Band κ_B2to cause an energy distribution in the field of the three-resonator system to have substantially no energy in the intermediate resonator structure while wirelessly transferring energy from the first resonator structure to the second resonator structure through the intermediate resonator structure.

Embodiments for the method above can include one or more of the following features.

Having substantially no energy in the intermediate resonator structure can mean that peak energy accumulation in the intermediate resonator structure is less than ten percent (10%) of the peak total energy in the three resonator structures throughout the wireless energy transfer.

Having substantially no energy in the intermediate resonator structure can mean that peak energy accumulation in the intermediate resonator structure is less than five percent (5%) of the peak total energy in the three resonator structures throughout the wireless energy transfer.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to maintain the energy distribution in the field of the three-resonator system in an eigenstate having the substantially no energy in the intermediate resonator structure.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to further cause the eigenstate to evolve substantially adiabatically from an initial state with substantially all energy in the resonator structures in the first resonator structure to a final state with substantially all of the energy in the resonator structures in the second resonator structure.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can include adjustments of both coupling rates κ_1Band κ_B2during wireless energy transfers.

The resonant energy in each of the resonator structures comprises electromagnetic fields. For example, the maximum value of the coupling rate κ_1Band the maximum value of the coupling rate κ_B2for inductive coupling between the intermediate resonator structure and each of the first and second resonator structures can each be larger than twice the loss rate Γ for each of the first and second resonators. Moreover, the maximum value of the coupling rate κ_1Band the maximum value of the coupling rate κ_B2for inductive coupling between the intermediate resonator structure and each of the first and second resonator structures can each be larger than four (4) times the loss rate Γ for each of the first and second resonators.

The resonator structure can have a resonant frequency between 50 KHz and 500 MHz.

The first and second resonator structures can be substantially identical.

Adjusting at least one of the coupling rates κ_1Band κ_B2can include adjusting a relative position and/or orientation between one or more pairs of the resonator structures.

Adjusting at least one of the coupling rates κ_1Band κ_B2can include adjusting a resonator property of one or more of the resonator structures, such as mutual inductance.

The resonator structures can include a capacitively loaded loop or coil of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

The resonator structures can include an inductively loaded rod of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

The wireless energy transfers can be non-radiative energy transfers mediated by a coupling of a resonant field evanescent tail of the first resonator structure and a resonant field evanescent tail of the intermediate resonator structure and a coupling of the resonant field evanescent tail of the intermediate resonator structure and a resonant field evanescent tail of the second resonator structure.

The first and second resonator structures can each have a quality factor greater than 50.

The first and second resonator structures can each have a quality factor greater than 100.

The adjustment of at least one of the coupling rates κ_1Band κ_B2can be selected to cause the energy distribution in the field of the three-resonator system to have substantially no energy in the intermediate resonator structure improves wireless energy transfer between the first and second resonator structures.

The adjustment of the at least one of the coupling rates can be selected to define a first mode of operation, wherein energy accumulation in the intermediate resonator structure during the wireless energy transfer from the first resonator structure to second resonator structure is smaller than that for a second mode of operation of wireless energy transfer among the three resonator structures having a coupling rate κ′_1Bfor wireless energy transfer from the first resonator structure to the intermediate resonator structure and a coupling rate κ′_B2for wireless energy transfer from the intermediate resonator structure to the second resonator structure with κ′_1Band κ′_B2each being substantially constant during the second mode of wireless energy transfer, and wherein the adjustment of the coupling rates κ_1Band κ_B2in the first mode of operation can be selected to satisfy κ^1B,κ_B2<√{square root over ((κ′_1B²+κ′_B2²)/2)}.

The first mode of operation can have a greater efficiency of energy transferred from the first resonator to the second resonator compared to that for the second mode of operation.

The first and second resonator structures can be substantially identical and each one can have a loss rate Γ_A, the intermediate resonator structure can have a loss rate Γ_B, and wherein Γ_B/Γ_Acan be greater than 50.

A ratio of energy lost to radiation and total energy wirelessly transferred between the first and second resonator structures in the first mode of operation can be less than that for the second mode of operation.

The first and second resonator structures can be substantially identical and each one can have a loss rate Γ_Aand a loss rate only due to radiation Γ_A,rad, the intermediate resonator structure can have a loss rate Γ_Band a loss rate only due to radiation Γ_B,radand wherein Γ_B,rad/Γ_B>Γ_A,rad/Γ_A.

The first mode of operation the intermediate resonator structure interacts less with extraneous objects than it does in the second mode of operation.

The coupling rate κ_B2can be maintained at a fixed value and the coupling rate κ_1Bcan be increased during the wireless energy transfer from the first resonator structure to second resonator structure.

The coupling rate κ_1Bcan be maintained at a fixed value and the coupling rate κ_B2can be decreased during the wireless energy transfer from the first resonator structure to second resonator structure.

During the wireless energy transfer from the first resonator structure to second resonator structure, the coupling rate κ_1Bcan be increased and the coupling rate κ_B2can be decreased.

The method may further include features corresponding to those listed for the apparatus and method described above, and one or more of the apparatuses and methods described below.

In another aspect, disclosed is an apparatus including: first, intermediate, and second resonator structures, wherein a coupling rate between the first resonator structure and the intermediate resonator structure is κ_1Band a coupling rate between the intermediate resonator structure and the second resonator structure is κ_B2; and means for adjusting at least one of the coupling rates κ_1Band κ_B2during wireless energy transfers among the resonator structures to cause an energy distribution in the field of the three-resonator system to have substantially no energy in the intermediate resonator structure while wirelessly transferring energy from the first resonator structure to the second resonator structure through the intermediate resonator structure.

Embodiments for the apparatus can include one or more of the following features.

The means for adjusting can be configured to maintain the energy distribution in the field of the three-resonator system in an eigenstate having the substantially no energy in the intermediate resonator structure.

The means for adjusting can include a rotation stage for adjusting the relative orientation of the intermediate resonator structure with respect to the first and second resonator structures.

The means for adjusting can include a translation stage for moving the first and/or second resonator structures relative to the intermediate resonator structure.

The means for adjusting can include a mechanical, electro-mechanical, or electrical staging system for dynamically adjusting the effective size of one or more of the resonator structures.

The resonator structures can include a capacitively loaded loop or coil of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

The resonator structures can include an inductively loaded rod of at least one of a conducting wire, a conducting Litz wire, and a conducting ribbon.

A source can be coupled to the first resonator structure and a load can be coupled to the second resonator structure.

The apparatus may further include features corresponding to those listed for the apparatus and methods described above, and the apparatus and method described below.

In another aspect, disclosed is a method for transferring energy wirelessly that includes: i) transferring energy wirelessly from a first resonator structure to a intermediate resonator structure, wherein the coupling rate between the first resonator structure and the intermediate resonator structure is κ_1B; ii) transferring energy wirelessly from the intermediate resonator structure to a second resonator, wherein the coupling rate between the intermediate resonator structure and the second resonator structure with a coupling rate is κ_B2; and iii) during the wireless energy transfers, adjusting at least one of the coupling rates κ_1Band κ_B2to define a first mode of operation in which energy accumulation in the intermediate resonator structure is reduced relative to that for a second mode of operation of wireless energy transfer among the three resonator structures having a coupling rate κ′_1Bfor wireless energy transfer from the first resonator structure to the intermediate resonator structure and a coupling rate κ′_B2for wireless energy transfer from the intermediate resonator structure to the second resonator structure with κ′_1Band κ′_B2each being substantially constant during the second mode of wireless energy transfer, and wherein the adjustment of the coupling rates κ_1Band κ_B2in the first mode of operation can be selected to satisfy κ_1B,κ_B2<√{square root over ((κ′_1B²+κ′_B2²)/2)}.

The method may further include features corresponding to those listed for the apparatuses and methods described above.

In another aspect, disclosed is an apparatus that includes: first, intermediate, and second resonator structures, wherein a coupling rate between the first resonator structure and the intermediate resonator structure is κ_1Band a coupling rate between the intermediate resonator structure and the second resonator structure is κ_B2; and means for adjusting at least one of the coupling rates κ_1Band κ_B2during wireless energy transfers among the resonator structures to define a first mode of operation in which energy accumulation in the intermediate resonator structure is reduced relative to that for a second mode of operation for wireless energy transfer among the three resonator structures having a coupling rate κ′_1Bfor wireless energy transfer from the first resonator structure to the intermediate resonator structure and a coupling rate κ′_B2for wireless energy transfer from the intermediate resonator structure to the second resonator structure with κ′_1Band κ′_B2each being substantially constant during the second mode of wireless energy transfer, and wherein the adjustment of the coupling rates κ₁₂and κ_B2in the first mode of operation can be selected to satisfy κ_1B,κ_B2<√{square root over ((κ′_1B²+κ′_B2²)/2)}.

The apparatus may further include features corresponding to those listed for the apparatuses and methods described above.

Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present disclosure, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.

The details of one or more embodiments are set forth in the accompanying drawings and the description below, including the documents appended hereto. Other features and advantages will be apparent from this disclosure and from the claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a schematic of an example wireless energy transfer scheme.

FIGS. 2(a)-(b) show the efficiency of power transmission η_pfor (a) U=1 and (b) U=3, as a function of the frequency detuning D_oand for different values of the loading rate U_o.

FIG. 2(c) shows the optimal (for zero detuning and under conditions of impedance matching) efficiency for energy transfer κ_E*and power transmission κ_P*, as a function of the coupling-to-loss figure-of-merit U.

FIG. 3 shows an example of a self-resonant conducting-wire coil.

FIG. 4 shows an example of a wireless energy transfer scheme featuring two self-resonant conducting-wire coils.

FIG. 5 is a schematic of an experimental system demonstrating wireless energy transfer.

FIG. 6 shows a comparison between experimental and theoretical results for the coupling rate of the system shown schematically in FIG. 5.

FIG. 7 shows a comparison between experimental and theoretical results for the strong-coupling factor of the system shown schematically in FIG. 5.

FIG. 8 shows a comparison between experimental and theoretical results for the power-transmission efficiency of the system shown schematically in FIG. 5.

FIG. 9 shows an example of a capacitively-loaded conducting-wire coil, and illustrates the surrounding field.

FIG. 10 shows an example wireless energy transfer scheme featuring two capacitively-loaded conducting-wire coils, and illustrates the surrounding field.

FIG. 11 illustrates an example circuit model for wireless energy transfer.

FIG. 12 shows the efficiency, total (loaded) device Q, and source and device currents, voltages and radiated powers (normalized to 1 Watt of output power to the load) as functions of the resonant frequency, for a particular choice of source and device loop dimensions, wp and N_sand different choices of N_d=1, 2, 3, 4, 5, 6, 10.

FIG. 13 shows the efficiency, total (loaded) device Q, and source and device currents, voltages and radiated powers (normalized to 1 Watt of output power to the load) as functions of frequency and wp for a particular choice of source and device loop dimensions, and number of turns N_sand N_d.

FIG. 14 shows an example of an inductively-loaded conducting-wire coil.

FIG. 15 shows (a) an example of a resonant dielectric disk, and illustrates the surrounding field and (b) a wireless energy transfer scheme featuring two resonant dielectric disks, and illustrates the surrounding field.

FIG. 16 shows a schematic of an example wireless energy transfer scheme with one source resonator and one device resonator exchanging energy indirectly through an intermediate resonator.

FIG. 17 shows an example of a wireless energy transfer system: (a) (Left) Schematic of loops configuration in two-object direct transfer. (Right) Time evolution of energies in the two-object direct energy transfer case. (b) (Left) Schematic of three-loops configuration in the constant-κ case. (Right) Dynamics of energy transfer for the configuration in (b. Left). Note that the total energy transferred E₂is 2 times larger than in (a. Right), but at the price of the total energy radiated being 4 times larger. (c) (Left) Loop configuration at t=0 in the adiabatic-κ scheme. (Center) Dynamics of energy transfer with adiabatically rotating loops. (Right) Loop configuration at t=t_EIT. Note that E₂is comparable to (b. Right), but the radiated energy is now much smaller: In fact, it is comparable to (a. Right).

FIG. 18 shows a schematic of an example wireless energy transfer scheme with one source resonator and one device resonator exchanging energy indirectly through an intermediate resonator, where an adjustment system is used to rotate the resonator structures to dynamically adjust their coupling rates.

FIG. 19 shows an example of a temporal variation of the coupling rates in a wireless energy transfer system as in FIG. 18 to achieve an adiabatic transfer of energy from the source object R₁to the device object R₂.

FIG. 20 shows the energy distribution in a wireless energy transfer system as in FIG. 18 as a function of time when the coupling rates are time-varying, for Γ_A=0, κ/Γ_B=10, κ_1B=κ sin [πt/(2t_EIT)], and κ_B2=κ cos [πt/(2t_EIT)].

FIGS. 21(a)-(f) show a comparison between the adiabatic-κ and constant-κ energy transfer schemes, in the general case: (a) Optimum E₂(%) in adiabatic-κ transfer, (b) Optimum E₂(%) in constant-κ transfer, (c) (E₂)_adiabatic-κ/(E₂)_constant-κ, (d) Energy lost (%) at optimum adiabatic-κ transfer, (e) Energy lost (%) at optimum constant-κ transfer, (f) (E_lost)_constant-κ/(E_lost)_adiabatic-κ.

FIG. 22(a)-(e) show a comparison between radiated energies in the adiabatic-κ and constant-κ energy transfer schemes: (a) E_rad(%) in the constant-scheme for Γ_B/Γ_A=500 and Γ_rad^A=0, (b) E_rad(%) in the adiabatic-κ scheme for Γ_B/Γ_A=500 and Γ_rad^A=0, (c) (E_rad)_constant-κ/(E_rad)_adiabatic-κ for Γ_B/Γ_A=50, (d) (E_rad)_constant-κ/(E_rad)_adiabatic-κ for Γ_B/Γ_A=500, (e) [(E_rad)_constant-κ/(E_rad)_adiabatic-κ] as function of κ/Γ_Band Γ_B/Γ_A, for Γ_rad^A=0.

FIGS. 23(a)-(b) show schematics for frequency control mechanisms.

FIGS. 24(a)-(c) illustrate a wireless energy transfer scheme using two dielectric disks in the presence of various extraneous objects.

DETAILED DESCRIPTION

Efficient wireless energy-transfer between two similar-frequency resonant objects can be achieved at mid-range distances, provided these resonant objects are designed to operate in the ‘strong-coupling’ regime. ‘Strong coupling’ can be realized for a wide variety of resonant objects, including electromagnetic resonant objects such as inductively-loaded conducting rods and dielectric disks. Recently, we have demonstrated wireless energy transfer between strongly coupled electromagnetic self-resonant conducting coils and capacitively-loaded conducting coils, bearing high-Q electromagnetic resonant modes. See, for example, the following commonly owned U.S. patent applications, all of which are incorporated herein by reference: U.S. application Ser. No. 11/481,077, filed on Jul. 5, 2006, and published as U.S. Patent Publication No. US 2007-0222542 A1; U.S. application Ser. No. 12/055,963, filed on Mar. 26, 2008, and published as U.S. Patent Publication No. US 2008-0278264 A1; and U.S. patent application Ser. No. 12/466,065, filed on May 14, 2009, and published as U.S. Patent Publication No. 2009-0284083. In general, the energy-transfer efficiency between similar-frequency, strongly coupled resonant objects decreases as the distance between the objects is increased.

In this work, we explore a further scheme of efficient energy transfer between resonant objects that extends the range over which energy may be efficiently transferred. Instead of transferring energy directly between two resonant objects, as has been described in certain embodiments of the cross-referenced patents, in certain embodiments, an intermediate resonant object, with a resonant frequency equal or nearly-equal to that of the two energy-exchanging resonant objects is used to mediate the transfer. The intermediate resonant object may be chosen so that it couples more strongly to each of the resonant objects involved in the energy transfer than those two resonant objects couple to each other. One way to design such an intermediate resonator is to make it larger than either of the resonant objects involved in the energy transfer. However, increasing the size of the intermediate resonant object may lower its quality factor, or Q, by increasing its radiation losses. Surprisingly enough, this new “indirect” energy transfer scheme may be shown to be very efficient and only weakly-radiative by introducing a meticulously chosen time variation of the resonator coupling rates.

The advantage of this method over the prior commonly owned wireless energy transfer techniques is that, in certain embodiments, it can enable energy to be transferred wirelessly between two objects with a larger efficiency and/or with a smaller radiation loss and/or with fewer interactions with extraneous objects.

Accordingly, in certain embodiments, we disclose an efficient wireless energy transfer scheme between two similar resonant objects, strongly coupled to an intermediate resonant object of substantially different properties, but with the same resonance frequency. The transfer mechanism essentially makes use of the adiabatic evolution of an instantaneous (so called ‘dark’) resonant state of the coupled three-object system. Our analysis is based on temporal coupled mode theory (CMT), and is general. Of particular commercial interest is the application of this technique to strongly-coupled electromagnetic resonators used for mid-range wireless energy transfer applications. We show that in certain parameter regimes of interest, this scheme can be more efficient, and/or less radiative than other wireless energy transfer approaches.

While the technique described herein is primarily directed to tangible resonator structures, the technique shares certain features with a quantum interference phenomenon known in the atomic physics community as Electromagnetically Induced Transparency (EIT). In EIT, three atomic states participate. Two of them, which are non-lossy, are coupled to one that has substantial losses. However, by meticulously controlling the mutual couplings between the states, one can establish a coupled system which is overall non-lossy. This phenomena has been demonstrated using carefully timed optical pulses, referred to as probe laser pulses and Stokes laser pulses, to reduce the opacity of media with the appropriate collection of atomic states. A closely related phenomenon known as Stimulated Raman Adiabatic Passage (STIRAP) may take place in a similar system; namely, the probe and Stokes laser beams may be used to achieve complete coherent population transfer between two molecular states of a medium. Hence, we may refer to the currently proposed scheme as the “EIT-like” energy transfer scheme.

In certain embodiments, we disclose an efficient near-field energy transfer scheme between two similar resonant objects, based on an EIT-like transfer of the energy through a mediating resonant object with the same resonant frequency. In embodiments, this EIT-like energy transfer may be realized using electromagnetic resonators as have been described in the cross-referenced patents, but the scheme is not bound only to wireless energy transfer applications. Rather, this scheme is general and may find applications in various other types of coupling between general resonant objects. In certain embodiments described below, we describe particular examples of electromagnetic resonators, but the nature of the resonators and their coupling mechanisms could be quite different (e.g. acoustic, mechanical, etc.). To the extent that many resonant phenomena can be modeled with nearly identical CMT equations, similar behavior to that described herein would occur.

1. Efficient Energy-Transfer by Two ‘Strongly Coupled’ Resonances

FIG. 1 shows a schematic that generally describes one example of the invention, in which energy is transferred wirelessly between two resonant objects. Referring to FIG. 1, energy is transferred over a distance D, between a resonant source object having a characteristic size r₁and a resonant device object of characteristic size r₂. Both objects are resonant objects. The wireless near-field energy transfer is performed using the field (e.g. the electromagnetic field or acoustic field) of the system of two resonant objects.

The characteristic size of an object can be regarded as being equal to the radius of the smallest sphere which can fit around the entire object. The characteristic thickness of an object can be regarded as being, when placed on a flat surface in any arbitrary configuration, the smallest possible height of the highest point of the object above a flat surface. The characteristic width of an object can be regarded as being the radius of the smallest possible circle that the object can pass through while traveling in a straight line. For example, the characteristic width of a cylindrical object is the radius of the cylinder.

Initially, we present a theoretical framework for understanding near-field wireless energy transfer. Note however that it is to be understood that the scope of the invention is not bound by theory.

Different temporal schemes can be employed, depending on the application, to transfer energy between two resonant objects. Here we will consider two particularly simple but important schemes: a one-time finite-amount energy-transfer scheme and a continuous finite-rate energy-transfer (power) scheme.

1.1 Finite-Amount Energy-Transfer Efficiency

Let the source and device objects be 1, 2 respectively and their resonance modes, which we will use for the energy exchange, have angular frequencies ω_1,2, frequency-widths due to intrinsic (absorption, radiation etc.) losses Γ_1,2and (generally) vector fields F_1,2(r), normalized to unity energy. Once the two resonant objects are brought in proximity, they can interact and an appropriate analytical framework for modeling this resonant interaction is that of the well-known coupled-mode theory (CMT). This model works well, when the resonances are well defined by having large quality factors and their resonant frequencies are relatively close to each other. In this picture, the field of the system of the two resonant objects 1, 2 can be approximated by F(r,t)=a₁(t)F₁(r)+a₂(t)F₂(r), where a_1,2(t) are the field amplitudes, with |a_1,2(t)|²equal to the energy stored inside the object 1, 2 respectively, due to the normalization. Then, using e^−iωttime dependence, the field amplitudes can be shown to satisfy, to lowest order:

$\begin{matrix} \frac{d}{d t} a_{1} (t) = - i (ω_{1} - i Γ_{1}) a_{1} (t) + i κ_{11} a_{1} (t) + i κ_{12} a_{2} (t) \frac{d}{d t} a_{2} (t) = - i (ω_{2} - i Γ_{2}) a_{2} (t) + i κ_{21} a_{1} (t) + i κ_{22} a_{2} (t) & (1) \end{matrix}$

where κ_11,22are the shifts in each object's frequency due to the presence of the other, which are a second-order correction and can be absorbed into the resonant frequencies (eigenfrequencies) by setting ω_1,2→ω_1,2+κ_11,22, and κ_12,21are the coupling coefficients, which from the reciprocity requirement of the system satisfy κ₂₁=κ₁₂≡κ.

The resonant modes of the combined system are found by substituting [a₁(t),a₂(t)]=[A₁,A₂]e^−iωt. They have complex resonant frequencies

ω_±=ω₁₂±√{square root over ((Δω₁₂)²+κ²)} (2a)

where ω₁₂=[(ω₁+ω₂)−i(Γ₁+Γ₂)]/2, Δω₁₂=[(ω₁−ω₂)−i(Γ¹−Γ₂)]/2 and whose splitting we denote as δ_E≡ω₊−ω₋, and corresponding resonant field amplitudes

$\begin{matrix} {\vec{V}}_{\pm} = {[\begin{matrix} A_{1} \\ A_{2} \end{matrix}]}_{\pm} = [\begin{matrix} κ \\ Δ ω_{12} \mp \sqrt{{(Δ ω_{12})}^{2} + κ^{2}} \end{matrix}] . & (2 b) \end{matrix}$

Note that, at exact resonance ω₁=ω₂=ω_Aand for Γ₁=Γ₂=Γ_A, we get Δω₁₂=0, δ_E=2κ, and then

${\overline{ω}}_{\pm} = ω_{A} \pm κ - i Γ_{A}$

${\vec{V}}_{\pm} = {[\begin{matrix} A_{1} \\ A_{2} \end{matrix}]}_{\pm} = [\begin{matrix} 1 \\ \mp 1 \end{matrix}],$

namely we get the known result that the resonant modes split to a lower frequency even mode and a higher frequency odd mode.

Assume now that at time t=0 the source object 1 has finite energy |a₁(0)|², while the device object has |a₂(0)|²=0. Since the objects are coupled, energy will be transferred from 1 to 2. With these initial conditions, Eqs. (1) can be solved, predicting the evolution of the device field-amplitude to be

$\begin{matrix} \frac{a_{2} (t)}{\langle a_{1} (0) \rangle} = \frac{2 κ}{δ_{E}} \sin (\frac{δ_{E} t}{2}) e^{- \frac{Γ_{1} + Γ_{2}}{2}} e^{- i \frac{ω_{1} - ω_{2}}{2}} . & (3) \end{matrix}$

The energy-transfer efficiency will be η_E≡|a₂(t)|²/|a₁(0)|². The ratio of energy converted to loss due to a specific loss mechanism in resonators 1 and 2, with respective loss rates Γ_1,lossand Γ_2,losswill be η_loss,E=∫₀^tdτ[2Γ_1,loss|a₁(τ)|²+2Γ_2,loss|a₂(τ)|²]/|a₁(0)|². Note that, at exact resonance ω₁=ω₂=ω_A(an optimal condition), Eq. (3) can be written as

$\begin{matrix} \frac{a_{2} (t)}{\langle a_{1} (0) \rangle} = \frac{\sin (\sqrt{1 - Δ^{2}} T)}{\sqrt{1 - Δ^{2}}} \cdot e^{- T / U} e^{- i ω_{A} t} & (4) \end{matrix}$

where ≡κt, Δ⁻¹=2κ/(Γ₂−Γ₁) and U=2κ/(Γ₁+Γ₂).

$\begin{matrix} T_{*} = \frac{\tan^{- 1} (U \sqrt{1 - Δ^{2}})}{\sqrt{1 - Δ^{2}}} & (5) \end{matrix}$

resulting in an optimal energy-transfer efficiency

$\begin{matrix} η_{E *} \equiv η_{E} (T_{*}) = \frac{U^{2}}{1 + U^{2} (1 - Δ^{2})} \exp (- \frac{2 \tan^{- 1} (U \sqrt{1 - Δ^{2}})}{U \sqrt{1 - Δ^{2}}}) . & (6 a) \end{matrix}$

which is a monotonically increasing function of the coupling-to-loss ratio U=2κ/(Γ₁+Γ₂) and tends to unity when U>>1 custom character |Δ|⁻¹>>1. Therefore, the energy transfer is nearly perfect, when the coupling rate is much faster than all loss rates (κ/Γ_1,2>1). In FIG. 2(c) we show the optimal energy-transfer efficiency when Γ₁=Γ₂=Γ_AΔ=0:

$\begin{matrix} η_{E} (T_{*}, Δ = 0) = \frac{U^{2}}{1 + U^{2}} \exp (- \frac{2 \tan^{- 1} U}{U}) . & (6 b) \end{matrix}$

In a real wireless energy-transfer system, the source object can be connected to a power generator (not shown in FIG. 1), and the device object can be connected to a power consuming load (e.g. a resistor, a battery, an actual device, not shown in FIG. 1). The generator will supply the energy to the source object, the energy will be transferred wirelessly and non-radiatively from the source object to the device object, and the load will consume the energy from the device object. To incorporate such supply and consumption mechanisms into this temporal scheme, in some examples, one can imagine that the generator is very briefly but very strongly coupled to the source at time t=0 to almost instantaneously provide the energy, and the load is similarly very briefly but very strongly coupled to the device at the optimal time t=t_*to almost instantaneously drain the energy. For a constant powering mechanism, at time t=t_*also the generator can again be coupled to the source to feed a new amount of energy, and this process can be repeated periodically with a period t_*.

1.2 Finite-Rate Energy-Transfer (Power-Transmission) Efficiency

Let the generator be continuously supplying energy to the source object 1 at a rate κ₁and the load continuously draining energy from the device object 2 at a rate κ₂. Field amplitudes s_±1,2(t) are then defined, so that |s_±1,2(t)|²is equal to the power ingoing to (for the + sign) or outgoing from (for the − sign) the object 1, 2 respectively, and the CMT equations are modified to

$\begin{matrix} \frac{d}{d t} a_{1} (t) = - i (ω_{1} - i Γ_{1}) a_{1} (t) + i κ_{11} a_{1} (t) + i κ_{12} a_{2} (t) - κ_{1} a_{1} (t) + \sqrt{2 κ_{1}} s_{+ 1} (t) \frac{d}{d t} a_{2} (t) = - i (ω_{2} - i Γ_{2}) a_{2} (t) + i κ_{21} a_{1} (t) + i κ_{22} a_{2} (t) - κ_{2} a_{2} (t) s_{- 1} (t) = \sqrt{2 κ_{1}} a_{1} (t) - s_{+ 1} (t) s_{- 2} (t) = \sqrt{2 κ_{2}} a_{2} (t) & (7) \end{matrix}$

where again we can set ω_1,2→ω_1,2+κ_11,22and κ₂₁=κ₁₂≡κ.

Assume now that the excitation is at a fixed frequency ω, namely has the form s₊₁(t)=S₊₁e^−iωt. Then the response of the linear system will be at the same frequency, namely a_1,2(t)=A₁₂e^−iωtand s_−1,2(t)=S_−1,2e^−iωt. By substituting these into Eqs. (7), using δ_1,2≡ω−ω_1,2, and solving the system, we find the field-amplitude transmitted to the load (S₂₁scattering-matrix element)

$\begin{matrix} \begin{matrix} S_{21} \equiv \frac{S_{- 2}}{S_{+ 1}} \\ = \frac{2 i κ \sqrt{κ_{1} κ_{2}}}{(Γ_{1} + κ_{1} - i δ_{1}) (Γ_{2} + κ_{2} - i δ_{2}) + κ^{2}} \\ = \frac{2 i U \sqrt{U_{1} U_{2}}}{(1 + U_{1} - i D_{1}) (1 + U_{2} - i D_{2}) + U^{2}} \end{matrix} & (8) \end{matrix}$

and the field-amplitude reflected to the generator (S₁₁scattering-matrix element)

$\begin{matrix} \begin{matrix} S_{11} \equiv \frac{S_{- 1}}{S_{+ 1}} \\ = \frac{(Γ_{1} - κ_{1} - i δ_{1}) (Γ_{2} + κ_{2} - i δ_{2}) + κ^{2}}{(Γ_{1} + κ_{1} - i δ_{1}) (Γ_{2} + κ_{2} - i δ_{2}) + κ^{2}} \\ = \frac{(1 - U_{1} - i D_{1}) (1 + U_{2} - i D_{2}) + U^{2}}{(1 + U_{1} - i D_{1}) (1 + U_{2} - i D_{2}) + U^{2}} \end{matrix} & (9) \end{matrix}$

where D_1,2≡δ_1,2/Γ_1,2, U_1,2≡κ_1,2/Γ_1,2and U≡κ/√{square root over (Γ₁Γ₂)}. Similarly, the scattering-matrix elements S₁₂, S₂₂are given by interchanging 1 custom character 2 in Eqs. (8), (9) and, as expected from reciprocity, S₂₁=S₁₂. The coefficients for power transmission (efficiency) and reflection and loss are respectively η_P≡|S₂₁|²=|S₋₂|²/—S₊₁|²and |S₁₁−²=|S₋₁−²/|S₊₁|²and 1−|S₂₁|²−|S₁₁|²=(2θ₁|A₁|²+2Γ₂|A₂|²)/|S₊₁|².

In some implementations, the parameters D_1,2, U_1,2can be designed (engineered), since one can adjust the resonant frequencies ω_1,2(compared to the desired operating frequency ω) and the generator/load supply/drain rates κ_1,2. Their choice can target the optimization of some system performance-characteristic of interest.

In some examples, a goal can be to maximize the power transmission (efficiency) η_P≡|S₂₁|²of the system, so one would require

η′_P(D_1,2)=η′_P(U_1,2)=0 (10)

Since S₂₁(from Eq. (8)) is symmetric upon interchanging 1 custom character 2, the optimal values for D_1,2(determined by Eqs. (10)) will be equal, namely D₁=D₂≡D_o, and similarly U₁=U₂≡U_o. Then,

$\begin{matrix} S_{21} = \frac{2 i {UU}_{o}}{{(1 + U_{o} - i D_{o})}^{2} + U^{2}} & (11) \end{matrix}$

and from the condition η′_P(D_o)=0 we get that, for fixed values of U and U_o, the efficiency can be maximized for the following values of the symmetric detuning

$\begin{matrix} D_{o} = \begin{matrix} \pm \sqrt{U^{2} - {(1 + U_{o})}^{2}}, & if & U > 1 + U_{o} \\ 0, & if & U \leq 1 + U_{o} \end{matrix}, & (12) \end{matrix}$

which, in the case U>1+U_o, can be rewritten for the two frequencies at which the efficiency peaks as

$\begin{matrix} ω_{\pm} = \frac{ω_{1} Γ_{2} + ω_{2} Γ_{1}}{Γ_{1} + Γ_{2}} \pm \frac{2 \sqrt{Γ_{1} Γ_{2}}}{Γ_{1} + Γ_{2}} \sqrt{κ^{2} - (Γ_{1} + κ_{1}) (Γ_{2} + κ_{2})}, & (13) \end{matrix}$

whose splitting we denote as δ_P≡ω₊−ω₋. Note that, at exact resonance ω₁=ω₂, and for Γ₁=Γ₂≡Γ_oand κ₁=κ₂≡κ_o, we get δ_P=2√{square root over (κ²−(Γ_o+κ_o)²)}<δ_E, namely the transmission-peak splitting is smaller than the normal-mode splitting. Then, by substituting D_ointo η_Pfrom Eq. (12), from the condition η′_P(U_o)=0 we get that, for fixed value of U, the efficiency can be maximized for

$\begin{matrix} U_{o *} = \sqrt{1 + U^{2}} \overset{Eq . (12)}{\Rightarrow} D_{o *} = 0 & (14) \end{matrix}$

which is known as ‘critical coupling’ condition, whereas for U_o<U_o*the system is called ‘undercoupled’ and for U_o>U_o*it is called ‘overcoupled’. The dependence of the efficiency on the frequency detuning D_ofor different values of U_o(including the ‘critical-coupling’ condition) are shown in FIG. 2(a,b). The overall optimal power efficiency using Eqs. (14) is

$\begin{matrix} η_{P *} \equiv η_{P} (D_{o *}, U_{o *}) = \frac{U_{o *} - 1}{U_{o *} + 1} = {(\frac{U}{1 + \sqrt{1 + U^{2}}})}^{2}, & (15) \end{matrix}$

which is again only a function of the coupling-to-loss ratio U=κ/√{square root over (Γ₁Γ₂)} and tends to unity when U>>1, as depicted in FIG. 2(c).

In some examples, a goal can be to minimize the power reflection at the side of the generator |s₁₁|²and the load |s₂₂|², so one would then need

S_11,22=0 custom character (1U₁−iD₁)(1±U₂−iD₂)+U²=0, (16)

The equations above present ‘impedance matching’ conditions. Again, the set of these conditions is symmetric upon interchanging 1 custom character 2, so, by substituting D₁=D₂≡D_oand U₁=U₂≡U_ointo Eqs. (16), we get

(1−iD_o)²−U_o²+U₂=0, (17)

from which we easily find that the values of D_oand U_othat cancel all reflections are again exactly those in Eqs. (14).

It can be seen that, the two goals and their associated sets of conditions (Eqs. (10) and Eqs. (16)) result in the same optimized values of the intra-source and intra-device parameters D_1,2, U_1,2. Note that for a lossless system this would be an immediate consequence of power conservation (Hermiticity of the scattering matrix), but this is not apparent for a lossy system.

Accordingly, for any temporal energy-transfer scheme, once the parameters specific only to the source or to the device (such as their resonant frequencies and their excitation or loading rates respectively) have been optimally designed, the efficiency monotonically increases with the ratio of the source-device coupling-rate to their loss rates. Using the definition of a resonance quality factor Q=ω/2Γ and defining by analogy the coupling factor k≡1/Q_κ=2κ/√{square root over (ω₁ω₂)}, it is therefore exactly this ratio

$\begin{matrix} U = \frac{κ}{\sqrt{Γ_{1} Γ_{2}}} = k \sqrt{Q_{1} Q_{2}} & (18) \end{matrix}$

that has been set as a figure-of-merit for any system under consideration for wireless energy-transfer, along with the distance over which this ratio can be achieved (clearly, U will be a decreasing function of distance). The operating regime U>1 is sometimes called ‘strong-coupling’ regime and is a sufficient condition for efficient energy-transfer. In particular, for U>1 we get, from Eq. (15), η_P*>17%, large enough for many practical applications. Note that in some applications, U>0.1 may be sufficient. In applications where it is impossible or impractical to run wires to supply power to a device, U<0.1 may be considered sufficient. One skilled in the art will recognize that the sufficient U is application and specification dependent. The figure-of-merit U may be called the strong-coupling factor. We will further show how to design systems with a large strong-coupling factor.

To achieve a large strong-coupling factor U, in some examples, the energy-transfer application preferably uses resonant modes of high quality factors Q, corresponding to low (i.e. slow) intrinsic-loss rates Γ. This condition can be satisfied by designing resonant modes where all loss mechanisms, typically radiation and absorption, are sufficiently suppressed.

This suggests that the coupling be implemented using, not the lossy radiative far-field, which should rather be suppressed, but the evanescent (non-lossy) stationary near-field. To implement an energy-transfer scheme, usually more appropriate are finite objects, namely ones that are topologically surrounded everywhere by air, into where the near field extends to achieve the coupling. Objects of finite extent do not generally support electromagnetic states that are exponentially decaying in all directions in air away from the objects, since Maxwell's Equations in free space imply that k²=ω²/c², where k is the wave vector, ω the angular frequency, and c the speed of light, because of which one can show that such finite objects cannot support states of infinite Q, rather there always is some amount of radiation. However, very long-lived (so-called “high-Q”) states can be found, whose tails display the needed exponential or exponential-like decay away from the resonant object over long enough distances before they turn oscillatory (radiative). The limiting surface, where this change in the field behavior happens, is called the “radiation caustic”, and, for the wireless energy-transfer scheme to be based on the near field rather than the far/radiation field, the distance between the coupled objects must be such that one lies within the radiation caustic of the other. One typical way of achieving a high radiation-Q (Q_rad) is to design subwavelength resonant objects. When the size of an object is much smaller than the wavelength of radiation in free space, its electromagnetic field couples to radiation very weakly. Since the extent of the near-field into the area surrounding a finite-sized resonant object is set typically by the wavelength, in some examples, resonant objects of subwavelength size have significantly longer evanescent field-tails. In other words, the radiation caustic is pushed far away from the object, so the electromagnetic mode enters the radiative regime only with a small amplitude.

Moreover, most realistic materials exhibit some nonzero amount of absorption, which can be frequency dependent, and thus cannot support states of infinite Q, rather there always is some amount of absorption. However, very long-lived (“high-Q”) states can be found, where electromagnetic modal energy is only weakly dissipated. Some typical ways of achieving a high absorption-Q (Q_abs) is to use materials which exhibit very small absorption at the resonant frequency and/or to shape the field to be localized more inside the least lossy materials.

Furthermore, to achieve a large strong-coupling factor U, in some examples, the energy-transfer application may use systems that achieve a high coupling factor k, corresponding to strong (i.e. fast) coupling rate κ, over distances larger than the characteristic sizes of the objects.

Since finite-sized subwavelength resonant objects can often be designed to have high Q, as was discussed above and will be seen in examples later on, such objects may typically be chosen for the resonant device-object. In these cases, the electromagnetic field is, in some examples, of a quasi-static nature and the distance, up to which sufficient coupling can be achieved, is dictated by the decay-law of this quasi-static field.

Note that in some examples, the resonant source-object may be immobile and thus less restricted in its allowed geometry and size. It can be therefore chosen to be large enough that the near-field extent is not limited by the wavelength, and can thus have nearly infinite radiation-Q. Some objects of nearly infinite extent, such as dielectric waveguides, can support guided modes, whose evanescent tails are decaying exponentially in the direction away from the object, slowly if tuned close to cutoff, therefore a good coupling can also be achieved over distances quite a few times larger than a characteristic size of the source- and/or device-object.

2 ‘Strongly-Coupled’ Resonances at Mid-Range Distances for Realistic Systems

In the following, examples of systems suitable for energy transfer of the type described above are described. We will demonstrate how to compute the CMT parameters ω_1,2, Q_1,2and k described above and how to choose or design these parameters for particular examples in order to produce a desirable figure-of-merit U=κ/√{square root over (Γ₁Γ₂)}=k√{square root over (Q₁Q₂)} at a desired distance D. In some examples, this figure-of-merit is maximized when ω_1,2are tuned close to a particular angular frequency ω_U.

2.1 Self-Resonant Conducting Coils

In some examples, one or more of the resonant objects are self-resonant conducting coils. Referring to FIG. 3, a conducting wire of length, l, and cross-sectional radius, a, is wound into a helical coil of radius, r, and height, h, (namely with N=√{square root over (l²−h²)}/2πr number of turns), surrounded by air. As described below, the wire has distributed inductance and distributed capacitance, and therefore it supports a resonant mode of angular frequency ω. The nature of the resonance lies in the periodic exchange of energy from the electric field within the capacitance of the coil, due to the charge distribution ρ(x) across it, to the magnetic field in free space, due to the current distribution j(x) in the wire. In particular, the charge conservation equation ∇·j=iωρ implies that: (i) this periodic exchange is accompanied by a π/2 phase-shift between the current and the charge density profiles, namely the energy W contained in the coil is at certain points in time completely due to the current and at other points in time completely due to the charge, and (ii) if ρ_l(x) and I(x) are respectively the linear charge and current densities in the wire, where x runs along the wire,

$q_{o} = \frac{1}{2} \int d x \langle ρ_{l} (x) \rangle$

is the maximum amount of positive charge accumulated in one side of the coil (where an equal amount of negative charge always also accumulates in the other side to make the system neutral) and I_o=max{|I(x)|} is the maximum positive value of the linear current distribution, then I_o=wq_o. Then, one can define an effective total inductance L and an effective total capacitance C of the coil through the amount of energy W inside its resonant mode:

$\begin{matrix} W \equiv \frac{1}{2} I_{o}^{2} L \Rightarrow L = \frac{μ_{o}}{4 π I_{o}^{2}} \int \int d x d x^{'} \frac{j (x) \cdot j (x^{'})}{\langle x - x^{'} \rangle}, & (19) \\ W \equiv \frac{1}{2} q_{o}^{2} \frac{1}{C} \Rightarrow \frac{1}{C} = \frac{1}{4 {πɛ}_{o} q_{o}^{2}} \int \int d x d x^{'} \frac{ρ (x) ρ (x^{'})}{\langle x - x^{'} \rangle}, & (20) \end{matrix}$

where μ_oand ∈_oare the magnetic permeability and electric permittivity of free space.

With these definitions, the resonant angular frequency and the effective impedance can be given by the formulas ω=1/√{square root over (LC)} and Z=√{square root over (L/C)} respectively.

Losses in this resonant system consist of ohmic (material absorption) loss inside the wire and radiative loss into free space. One can again define a total absorption resistance R_absfrom the amount of power absorbed inside the wire and a total radiation resistance R_radfrom the amount of power radiated due to electric- and magnetic-dipole radiation:

$\begin{matrix} P_{abs} \equiv \frac{1}{2} I_{o}^{2} R_{abs} \Rightarrow R_{abs} \approx ζ_{c} \frac{l}{2 π a} \cdot \frac{I_{rms}^{2}}{I_{o}^{2}} & (21) \\ P_{rad} \equiv \frac{1}{2} I_{o}^{2} R_{rad} \Rightarrow R_{rad} \approx \frac{ζ_{o}}{6 π} [{(\frac{ω \langle p \rangle}{c})}^{2} + {(\frac{ω \sqrt{\langle m \rangle}}{c})}^{4}], & (22) \end{matrix}$

where c=1/√{square root over (μ_o∈_o)} and ζ_o=√{square root over (μ_o/∈_o)} are the light velocity and light impedance in free space, the impedance ζ_cis ζ_c=1/σζ=√{square root over (μ_oω/2σ)} with σ the conductivity of the conductor and δ the skin depth at the frequency ω,

$I_{rms}^{2} = \frac{1}{l} \int d x {\langle I (x) \rangle}^{2}, p = \int d x r ρ_{l} (x)$

is the electric-dipole moment of the coil and

$m = \frac{1}{2} \int d x r \times j (x)$

is the magnetic-dipole moment of the coil. For the radiation resistance formula Eq. (22), the assumption of operation in the quasi-static regime (h,r<<λ=2πc/ω) has been used, which is the desired regime of a subwavelength resonance. With these definitions, the absorption and radiation quality factors of the resonance may be given by Q_abs=Z/R_absand Q_rad=Z/R_radrespectively.

From Eq. (19)-(22) it follows that to determine the resonance parameters one simply needs to know the current distribution j in the resonant coil. Solving Maxwell's equations to rigorously find the current distribution of the resonant electromagnetic eigenmode of a conducting-wire coil is more involved than, for example, of a standard LC circuit, and we can find no exact solutions in the literature for coils of finite length, making an exact solution difficult. One could in principle write down an elaborate transmission-line-like model, and solve it by brute force. We instead present a model that is (as described below) in good agreement (˜5%) with experiment. Observing that the finite extent of the conductor forming each coil imposes the boundary condition that the current has to be zero at the ends of the coil, since no current can leave the wire, we assume that the resonant mode of each coil is well approximated by a sinusoidal current profile along the length of the conducting wire. We shall be interested in the lowest mode, so if we denote by x the coordinate along the conductor, such that it runs from −l/2 to +l/2, then the current amplitude profile would have the form I(x)=I_ocos(πx/l), where we have assumed that the current does not vary significantly along the wire circumference for a particular x, a valid assumption provided a<<r. It immediately follows from the continuity equation for charge that the linear charge density profile should be of the form ρ_l(x)=ρ_osin(πx/l), and thus q_o=∫₀^l/2dxρ_o| sin(πx/l)|=ρ_ol/π. Using these sinusoidal profiles we find the so-called “self-inductance” L_sand “self-capacitance” C_sof the coil by computing numerically the integrals Eq. (19) and (20); the associated frequency and effective impedance are ω_sand Z_srespectively. The “self-resistances” R_sare given analytically by Eq. (21) and (22) using

$I_{rms}^{2} = \frac{1}{l} \int_{- l / 2}^{l / 2} d x {\langle I_{o} \cos (π x / l) \rangle}^{2} = \frac{1}{2} I_{o}^{2}, \langle p \rangle = q_{o} \sqrt{{(\frac{2}{π} h)}^{2} + {(\frac{4 N \cos (π N)}{(4 N^{2} - 1) π} r)}^{2}}$

$and$

$\langle m \rangle = I_{o} \sqrt{{(\frac{2}{π} N π r^{2})}^{2} + {(\frac{\cos (π N) (12 N^{2} - 1) - \sin (π N) π N (4 N^{2} - 1)}{(16 N^{4} - 8 N^{2} + 1) π})}^{2}},$

and therefore the associated Q_sfactors can be calculated.

The results for two examples of resonant coils with subwavelength modes of λ_s/r≧70 (i.e. those highly suitable for near-field coupling and well within the quasi-static limit) are presented in Table 1. Numerical results are shown for the wavelength and absorption, radiation and total loss rates, for the two different cases of subwavelength-coil resonant modes. Note that, for conducting material, copper (σ=5.998·10^−7 S/m) was used. It can be seen that expected quality factors at microwave frequencies are Q_s,abs≧1000 and Q_s,rad≧5000.

TABLE 1

single coil
λ_s/r
f (MHz)
Q_s,rad
Q_s,abs
Q_s

r = 30 cm, h = 20 cm, a = 1 cm,
74.7
13.39
4164
8170
2758

N = 4

r = 10 cm, h = 3 cm, a = 2 mm,
140
21.38
43919
3968
3639

N = 6

Referring to FIG. 4, in some examples, energy is transferred between two self-resonant conducting-wire coils. The electric and magnetic fields are used to couple the different resonant conducting-wire coils at a distance D between their centers. Usually, the electric coupling highly dominates over the magnetic coupling in the system under consideration for coils with h>>2r and, oppositely, the magnetic coupling highly dominates over the electric coupling for coils with h<<2r. Defining the charge and current distributions of two coils 1,2 respectively as ρ_1,2(x) and j_1,2(x), total charges and peak currents respectively as q_1,2and I_1,2, and capacitances and inductances respectively as C_1,2and L_1,2, which are the analogs of ρ(x), j(x), q_o, I_o, C and L for the single-coil case and are therefore well defined, we can define their mutual capacitance and inductance through the total energy:

$\begin{matrix} W \equiv W_{1} + W_{2} + \frac{1}{2} (q_{1}^{*} q_{2} + q_{2}^{*} q_{1}) / M_{C} + \frac{1}{2} (I_{1}^{*} I_{2} + I_{2}^{*} I_{1}) M_{L} \Rightarrow 1 / M_{C} = \frac{1}{4 {πɛ}_{o} q_{1} q_{2}} \int \int d x d x^{'} \frac{ρ_{1} (x) ρ_{2} (x^{'})}{\langle x - x^{'} \rangle} u, M_{L} = \frac{μ_{o}}{4 π I_{1} I_{2}} \int \int d x d x^{'} \frac{j_{1} (x) \cdot j_{2} (x^{'})}{\langle x - x^{'} \rangle} u, where W_{1} = \frac{1}{2} q_{1}^{2} / C_{1} = \frac{1}{2} I_{1}^{2} L_{1}, W_{2} = \frac{1}{2} q_{2}^{2} / C_{2} = \frac{1}{2} I_{2}^{2} L_{2} & (23) \end{matrix}$

and the retardation factor of u=exp(iω|x−x′|/c) inside the integral can been ignored in the quasi-static regime D<<λ of interest, where each coil is within the near field of the other. With this definition, the coupling factor is given by k=√{square root over (C₁C₂)}/M_C+M_L/√{square root over (L₁L₂)}.

Therefore, to calculate the coupling rate between two self-resonant coils, again the current profiles are needed and, by using again the assumed sinusoidal current profiles, we compute numerically from Eq. (23) the mutual capacitance M_C,sand inductance M_L,sbetween two self-resonant coils at a distance D between their centers, and thus k=1/Q_κ is also determined.

TABLE 2

pair of coils
D/r
Q
Q_κ = 1/k
U

r = 30 cm, h = 20 cm,
3
2758
38.9
70.9

a = 1 cm, N = 4
5
2758
139.4
19.8

λ/r ≈ 75
7
2758
333.0
8.3

Q_s^abs≈ 8170, Q_s^rad≈ 4164
10
2758
818.9
3.4

r = 10 cm, h = 3 cm,
3
3639
61.4
59.3

a = 2 mm, N = 6
5
3639
232.5
15.7

λ/r ≈ 140
7
3639
587.5
6.2

Q_s^abs≈ 3968, Q_s^rad≈ 43919
10
3639
1580
2.3

Referring to Table 2, relevant parameters are shown for exemplary examples featuring pairs or identical self resonant coils. Numerical results are presented for the average wavelength and loss rates of the two normal modes (individual values not shown), and also the coupling rate and figure-of-merit as a function of the coupling distance D, for the two cases of modes presented in Table 1. It can be seen that for medium distances D/r=10−3 the expected coupling-to-loss ratios are in the range U˜2-70.

2.1.1 Experimental Results

An experimental realization of an example of the above described system for wireless energy transfer consists of two self-resonant coils, one of which (the source coil) is coupled inductively to an oscillating circuit, and the second (the device coil) is coupled inductively to a resistive load, as shown schematically in FIG. 5. Referring to FIG. 5, A is a single copper loop of radius 25 cm that is part of the driving circuit, which outputs a sine wave with frequency 9.9 MHz. s and d are respectively the source and device coils referred to in the text. B is a loop of wire attached to the load (“light-bulb”). The various κ's represent direct couplings between the objects. The angle between coil d and the loop A is adjusted so that their direct coupling is zero, while coils s and d are aligned coaxially. The direct coupling between B and A and between B and s is negligible.

The parameters for the two identical helical coils built for the experimental validation of the power transfer scheme were h=20 cm, a=3 mm, r=30 cm and N=5.25. Both coils are made of copper. Due to imperfections in the construction, the spacing between loops of the helix is not uniform, and we have encapsulated the uncertainty about their uniformity by attributing a 10% (2 cm) uncertainty to h. The expected resonant frequency given these dimensions is f₀=10.56±0.3 MHz, which is approximately 5% off from the measured resonance at around 9.90 MHz.

The theoretical Q for the loops is estimated to be ˜2500 (assuming perfect copper of resistivity ρ=1/σ=1.7×10⁻⁸Ωm) but the measured value is 950±50. We believe the discrepancy is mostly due to the effect of the layer of poorly conducting copper oxide on the surface of the copper wire, to which the current is confined by the short skin depth (˜20 μm) at this frequency. We have therefore used the experimentally observed Q (and Γ₁=Γ₂=Γ=ω/(2Q) derived from it) in all subsequent computations.

The coupling coefficient κ can be found experimentally by placing the two self-resonant coils (fine-tuned, by slightly adjusting h, to the same resonant frequency when isolated) a distance D apart and measuring the splitting in the frequencies of the two resonant modes in the transmission spectrum. According to Eq. (13) derived by coupled-mode theory, the splitting in the transmission spectrum should be δ_P=2√{square root over (κ²−Γ²)}, when κ_A,Bare kept very small by keeping A and B at a relatively large distance. The comparison between experimental and theoretical results as a function of distance when the two the coils are aligned coaxially is shown in FIG. 6.

FIG. 7 shows a comparison of experimental and theoretical values for the strong-coupling factor U=κ/Γ as a function of the separation between the two coils. The theory values are obtained by using the theoretically obtained κ and the experimentally measured Γ. The shaded area represents the spread in the theoretical U due to the ˜5% uncertainty in Q. As noted above, the maximum theoretical efficiency depends only on the parameter U, which is plotted as a function of distance in FIG. 7. U is greater than 1 even for D=2.4 m (eight times the radius of the coils), thus the system is in the strongly-coupled regime throughout the entire range of distances probed.

The power-generator circuit was a standard Colpitts oscillator coupled inductively to the source coil by means of a single loop of copper wire 25 cm in radius (see FIG. 5). The load consisted of a previously calibrated light-bulb, and was attached to its own loop of insulated wire, which was in turn placed in proximity of the device coil and inductively coupled to it. Thus, by varying the distance between the light-bulb and the device coil, the parameter U_B=κ_B/Γ was adjusted so that it matched its optimal value, given theoretically by Eq. (14) as U_B*=√{square root over (1+U²)}. Because of its inductive nature, the loop connected to the light-bulb added a small reactive component to κ_Bwhich was compensated for by slightly retuning the coil. The work extracted was determined by adjusting the power going into the Colpitts oscillator until the light-bulb at the load was at its full nominal brightness.

In order to isolate the efficiency of the transfer taking place specifically between the source coil and the load, we measured the current at the mid-point of each of the self-resonant coils with a current-probe (which was not found to lower the Q of the coils noticeably.) This gave a measurement of the current parameters I₁and I₂defined above. The power dissipated in each coil was then computed from P_1,2=ΓL|I_1,2|², and the efficiency was directly obtained from η=P_B/(P₁+P₂+P_B). To ensure that the experimental setup was well described by a two-object coupled-mode theory model, we positioned the device coil such that its direct coupling to the copper loop attached to the Colpitts oscillator was zero. The experimental results are shown in FIG. 8, along with the theoretical prediction for maximum efficiency, given by Eq. (15).

Using this example, we were able to transmit significant amounts of power using this setup from the source coil to the device coil, fully lighting up a 60 W light-bulb from distances more than 2 m away, for example. As an additional test, we also measured the total power going into the driving circuit. The efficiency of the wireless power-transmission itself was hard to estimate in this way, however, as the efficiency of the Colpitts oscillator itself is not precisely known, although it is expected to be far from 100%. Nevertheless, this gave an overly conservative lower bound on the efficiency. When transmitting 60 W to the load over a distance of 2 m, for example, the power flowing into the driving circuit was 400 W. This yields an overall wall-to-load efficiency of ˜15%, which is reasonable given the expected ˜40% efficiency for the wireless power transmission at that distance and the low efficiency of the driving circuit.

From the theoretical treatment above, we see that in typical examples it is important that the coils be on resonance for the power transmission to be practical. We found experimentally that the power transmitted to the load dropped sharply as one of the coils was detuned from resonance. For a fractional detuning Δf/f₀of a few times the inverse loaded Q, the induced current in the device coil was indistinguishable from noise.

The power transmission was not found to be visibly affected as humans and various everyday objects, such as metallic and wooden furniture, as well as electronic devices large and small, were placed between the two coils, even when they drastically obstructed the line of sight between source and device. External objects were found to have an effect only when they were closer than 10 cm from either one of the coils. While some materials (such as aluminum foil, styrofoam and humans) mostly just shifted the resonant frequency, which could in principle be easily corrected with a feedback circuit of the type described earlier, others (cardboard, wood, and PVC) lowered Q when placed closer than a few centimeters from the coil, thereby lowering the efficiency of the transfer.

This method of power transmission is believed safe for humans. When transmitting 60 W (more than enough to power a laptop computer) across 2 m, we estimated that the magnitude of the magnetic field generated is much weaker than the Earth's magnetic field for all distances except for less than about 1 cm away from the wires in the coil, an indication of the safety of the scheme even after long-term use. The power radiated for these parameters was ˜5 W, which is roughly an order of magnitude higher than cell phones but could be drastically reduced, as discussed below.

Although the two coils are currently of identical dimensions, it is possible to make the device coil small enough to fit into portable devices without decreasing the efficiency. One could, for instance, maintain the product of the characteristic sizes of the source and device coils constant.

These experiments demonstrated experimentally a system for power transmission over medium range distances, and found that the experimental results match theory well in multiple independent and mutually consistent tests.

The efficiency of the scheme and the distances covered can be improved by silver-plating the coils, which may increase their Q, or by working with more elaborate geometries for the resonant objects. Nevertheless, the performance characteristics of the system presented here are already at levels where they could be useful in practical applications.

2.2 Capacitively-Loaded Conducting Loops or Coils

In some examples, one or more of the resonant objects are capacitively-loaded conducting loops or coils. Referring to FIG. 9 a helical coil with N turns of conducting wire, as described above, is connected to a pair of conducting parallel plates of area A spaced by distance d via a dielectric material of relative permittivity e, and everything is surrounded by air (as shown, N=1 and h=0). The plates have a capacitance C_p=∈_o∈A/d, which is added to the distributed capacitance of the coil and thus modifies its resonance. Note however, that the presence of the loading capacitor may modify the current distribution inside the wire and therefore the total effective inductance L and total effective capacitance C of the coil may be different respectively from L_sand C_s, which are calculated for a self-resonant coil of the same geometry using a sinusoidal current profile. Since some charge may accumulate at the plates of the external loading capacitor, the charge distribution ρ inside the wire may be reduced, so C<C_s, and thus, from the charge conservation equation, the current distribution j may flatten out, so L>L_s. The resonant frequency for this system may be ω=1/√{square root over (L(C+C_p))}<ω_s=1/√{square root over (L_sC_s)}, and I(x)→I_ocos(λx/l) custom character C→C_sω→ω_s,as C_p→0.

In general, the desired CMT parameters can be found for this system, but again a very complicated solution of Maxwell's Equations is required. Instead, we will analyze only a special case, where a reasonable guess for the current distribution can be made. When C_p>>C_s>C, then ω≈1/√{square root over (LC_p)}<<ω_sand Z≈√{square root over (L/C_p)}<<Z_s, while all the charge is on the plates of the loading capacitor and thus the current distribution is constant along the wire. This allows us now to compute numerically L from Eq. (19). In the case h=0 and N integer, the integral in Eq. (19) can actually be computed analytically, giving the formula L=μ_or[ln(8r/a)−2]N². Explicit analytical formulas are again available for R from Eq. (21) and (22), since I_rms=I_o, |p|≈0 and |m|=I_oNπr²(namely only the magnetic-dipole term is contributing to radiation), so we can determine also Q_abs=ωL/R_absand Q_rad=ωL/R_rad. At the end of the calculations, the validity of the assumption of constant current profile is confirmed by checking that indeed the condition C_p>>C_s custom character ω<<ω_sis satisfied. To satisfy this condition, one could use a large external capacitance, however, this would usually shift the operational frequency lower than the optimal frequency, which we will determine shortly; instead, in typical examples, one often prefers coils with very small self-capacitance C_sto begin with, which usually holds, for the types of coils under consideration, when N=1, so that the self-capacitance comes from the charge distribution across the single turn, which is almost always very small, or when N>1 and h>>2Na, so that the dominant self-capacitance comes from the charge distribution across adjacent turns, which is small if the separation between adjacent turns is large.

The external loading capacitance C_Pprovides the freedom to tune the resonant frequency (for example by tuning A or d). Then, for the particular simple case h=0, for which we have analytical formulas, the total Q=ωL/(R_abs+R_rad) becomes highest at the optimal frequency

$\begin{matrix} ω_{Q} = {[\frac{c^{4}}{π} \sqrt{\frac{ɛ_{o}}{2 σ}} \cdot \frac{1}{{aNr}^{3}}]}^{\frac{2}{7}}, & (24) \end{matrix}$

reaching the value

$\begin{matrix} Q_{\max} = \frac{6}{7 π} {(2 π^{2} η_{o} \frac{σ a^{2} N^{2}}{r})}^{\frac{3}{7}} \cdot [\ln (\frac{8 r}{a}) - 2] . & (25) \end{matrix}$

At lower frequencies Q is dominated by ohmic loss and at higher frequencies by radiation. Note, however, that the formulas above are accurate as long as ω_Q<<ω_sand, as explained above, this holds almost always when N=1, and is usually less accurate when N>1, since h=0 usually implies a large self-capacitance. A coil with large h can be used, if the self-capacitance needs to be reduced compared to the external capacitance, but then the formulas for L and ω_Q, Q_maxare again less accurate. Similar qualitative behavior is expected, but a more complicated theoretical model is needed for making quantitative predictions in that case.

The results of the above analysis for two examples of subwavelength modes of λ/r≧70 (namely highly suitable for near-field coupling and well within the quasi-static limit) of coils with N=1 and h=0 at the optimal frequency Eq. (24) are presented in Table 3. To confirm the validity of constant-current assumption and the resulting analytical formulas, mode-solving calculations were also performed using another completely independent method: computational 3D finite-element frequency-domain (FEFD) simulations (which solve Maxwell's Equations in frequency domain exactly apart for spatial discretization) were conducted, in which the boundaries of the conductor were modeled using a complex impedance ζ_c=√{square root over (μ_oω/2σ)} boundary condition, valid as long as ζ_c/ζ_o<<1 (<10⁻⁵for copper in the microwave). Table 3 shows Numerical FEFD (and in parentheses analytical) results for the wavelength and absorption, radiation and total loss rates, for two different cases of subwavelength-loop resonant modes. Note that copper was used for the conducting material (σ=5.998·10⁷S/m). Specific parameters of the plot in FIG. 4 are highlighted in bold in the table. The two methods (analytical and computational) are in good agreement and show that, in some examples, the optimal frequency is in the low-MHz microwave range and the expected quality factors are Q_abs≧1000 and Q_rad≧10000.

TABLE 3

single coil
λ/r
f
Q_rad
Q_abs
Q

r = 30 cm, a = 2 cm

111.4

8.976

29546

4886

4193

ε = 10, A = 138 cm²_,
(112.4)
(8.897)
(30512)
(5117)
(4381)

d = 4 mm

r = 10 cm, a = 2 mm
69.7
43.04
10702
1545
1350

ε = 10, A = 3.14 cm²_,
(70.4)
(42.61)
(10727)
(1604)
(1395)

d = 1 mm

Referring to FIG. 10, in some examples, energy is transferred between two capacitively-loaded coils. For the rate of energy transfer between two capacitively-loaded coils 1 and 2 at distance D between their centers, the mutual inductance M_Lcan be evaluated numerically from Eq. (23) by using constant current distributions in the case ω<<ω_s. In the case h=0, the coupling may be only magnetic and again we have an analytical formula, which, in the quasi-static limit r<<D<<λ and for the relative orientation shown in FIG. 10, is M_L≈πμ_o/2·(r₁r₂)²N₁N₂/D³, which means that k∝(√{square root over (r₁r₂)}/D)³may be independent of the frequency ω and the number of turns N₁, N₂. Consequently, the resultant coupling figure-of-merit of interest is

$\begin{matrix} U = k \sqrt{Q_{1} Q_{2}} \approx {(\frac{\sqrt{r_{1} r_{2}}}{D})}^{3} \cdot \frac{π^{2} η_{o} \frac{\sqrt{r_{1} r_{2}}}{λ} \cdot N_{1} N_{2}}{\prod_{j = 1, 2} {(\frac{\sqrt{{πη}_{o}}}{λσ} \cdot \frac{r_{j}}{a_{j}} N_{j} + \frac{8}{3} π^{5} {η_{o} (\frac{r_{j}}{λ})}^{4} N_{j}^{2})}^{\frac{1}{2}}}, & (26) \end{matrix}$

which again is more accurate for N₁=N₂=1.

From Eq. (26) it can be seen that the optimal frequency ω_U, where the figure-of-merit is maximized to the value U_max, is close to the frequency ω_Q₁_Q₂at which Q₁Q₂is maximized, since k does not depend much on frequency (at least for the distances D<<λ of interest for which the quasi-static approximation is still valid). Therefore, the optimal frequency ω_U≈ω_Q₁_Q₂may be mostly independent of the distance D between the two coils and may lie between the two frequencies ω_Q₁and ω_Q₂at which the single-coil Q₁and Q₂respectively peak. For same coils, this optimal frequency is given by Eq. (24) and then the strong-coupling factor from Eq. (26) becomes

$\begin{matrix} U_{\max} = {kQ}_{\max} \approx {(\frac{r}{D})}^{3} \cdot \frac{3}{7} {(2 π^{2} η_{o} \frac{σ a^{2} N^{2}}{r})}^{\frac{3}{7}} . & (27) \end{matrix}$

In some examples, one can tune the capacitively-loaded conducting loops or coils, so that their angular resonant frequencies are close to ω_Uwithin Γ_U, which is half the angular frequency width for which U>U_max/2.

Referring to Table 4, numerical FEFD and, in parentheses, analytical results based on the above are shown for two systems each composed of a matched pair of the loaded coils described in Table 3. The average wavelength and loss rates are shown along with the coupling rate and coupling to loss ratio figure-of-merit U=κ/Γ as a function of the coupling distance D, for the two cases. Note that the average numerical Γ_radshown are slightly different from the single-loop value of FIG. 3. Analytical results for Γ_radare not shown but the single-loop value is used. (The specific parameters corresponding to the plot in FIG. 10 are highlighted with bold in the table.) Again we chose N=1 to make the constant-current assumption a good one and computed M_Lnumerically from Eq. (23). Indeed the accuracy can be confirmed by their agreement with the computational FEFD mode-solver simulations, which give κ through the frequency splitting of the two normal modes of the combined system (δ_E=2κ from Eq. (4)). The results show that for medium distances D/r=10−3 the expected coupling-to-loss ratios are in the range U˜0.5-50.

TABLE 4

pair of coils
D/r
Q^rad
Q = ω/2Γ
Q_κ = ω/2κ
κ/Γ

r = 30 cm, a = 2 cm
3
30729
4216
62.6 (63.7)
67.4 (68.7)

ε = 10, A = 138 cm²,
5
29577
4194
235 (248)
17.8 (17.6)

d = 4 mm

7

29128

4185

589 (646)

7.1 (6.8)

λ/r ≈ 112
10
28833
4177
1539 (1828)
2.7 (2.4)

Q^abs≈ 4886

r = 10 cm, a = 2 mm
3
10955
1355
85.4 (91.3)
15.9 (15.3)

ε = 10, A = 3.14 cm²,
5
10740
1351
313 (356)
4.32 (3.92)

d = 1 mm
7
10759
1351
754 (925)
1.79 (1.51)

λ/r ≈ 70
10
10756
1351
1895 (2617)
0.71 (0.53)

Q^abs≈ 1546

2.2.1 Derivation of Optimal Power-Transmission Efficiency

Referring to FIG. 11, to rederive and express Eq. (15) in terms of the parameters which are more directly accessible from particular resonant objects, such as the capacitively-loaded conducting loops, one can consider the following circuit-model of the system, where the inductances L_s,L_drepresent the source and device loops respectively, R_s,R_dtheir respective losses, and C_s,C_dare the required corresponding capacitances to achieve for both resonance at frequency ω. A voltage generator V_gis considered to be connected to the source and a load resistance R_lto the device. The mutual inductance is denoted by M.

Then from the source circuit at resonance (wL_s=1/wC_s):

$\begin{matrix} V_{g} = I_{s} R_{s} - j ω {MI}_{d} \Rightarrow \frac{1}{2} V_{g}^{*} I_{s} = \frac{1}{2} {\langle I_{s} \rangle}^{2} R_{s} + \frac{1}{2} j ω {MI}_{d}^{*} I_{s}, & (28) \end{matrix}$

and from the device circuit at resonance (wL_d=1/wC_d)

0=I_d(R_d+R_l)−jwMI_s custom character jwMI_s=I_d(R_d+R_l) (29)

So by substituting Eq. (29) to Eq. (28) and taking the real part (for time-averaged power) we get:

$\begin{matrix} P_{g} = Re {\frac{1}{2} V_{g}^{*} I_{s}} = \frac{1}{2} {\langle I_{s} \rangle}^{2} R_{s} + \frac{1}{2} {\langle I_{d} \rangle}^{2} (R_{d} + R_{l}) = P_{s} + P_{d} + P_{l}, & (30) \end{matrix}$

where we identified the power delivered by the generator P_g=Re{V*_gI_s/2}, the power lost inside the source P_s=|I_s|²R_s/2, the power lost inside the device P_d=|I_d|²R_d/2 and the power delivered to the load P_l=|I_d|₂R_l/2. Then, the power transmission efficiency is:

$\begin{matrix} η_{P} \equiv \frac{P_{l}}{P_{g}} = \frac{R_{l}}{{\langle \frac{I_{s}}{I_{d}} \rangle}^{2} R_{s} + (R_{d} + R_{l})} \overset{(29)}{=} \frac{R_{l}}{\frac{{(R_{d} + R_{l})}^{2}}{{(ω M)}^{2}} R_{s} + (R_{d} + R_{l})} . & (31) \end{matrix}$

If we now choose the load impedance R_lto optimize the efficiency by η′_P(R_l)=0, we get the optimal load impedance

$\begin{matrix} \frac{R_{l^{*}}}{R_{d}} = \sqrt{1 + \frac{{(ω M)}^{2}}{R_{s} R_{d}}} & (32) \end{matrix}$

and the maximum possible efficiency

$\begin{matrix} η_{P^{*}} = \frac{R_{l^{*}} / R_{d} - 1}{R_{l^{*}} / R_{d} + 1} = {[\frac{ω M / \sqrt{R_{s} R_{d}}}{1 + \sqrt{1 + {(ω M / \sqrt{R_{s} R_{d}})}^{2}}}]}^{2} . & (33) \end{matrix}$

To check now the correspondence with the CMT model, note that κ_l=R_l/2L_d, Γ_d=R_d/2L_d, Γ_s=R_s/2L_s, and κ=ωM/2√{square root over (L_sL_d)}, so then U_l=κ_l/Γ_d=R_l/R_dand U=κ/√{square root over (Γ_sΓ_d)}=ωM/√{square root over (R_sR_d)}. Therefore, the condition Eq. (32) is identical to the condition Eq. (14) and the optimal efficiency Eq. (33) is identical to the general Eq. (15). Indeed, as the CMT analysis predicted, to get a large efficiency, we need to design a system that has a large strong-coupling factor U.

2.2.2 Optimization of U

The results above can be used to increase or optimize the performance of a wireless energy transfer system, which employs capacitively-loaded coils. For example, from the scaling of Eq. (27) with the different system parameters, one sees that to maximize the system figure-of-merit U, in some examples, one can:

- Decrease the resistivity of the conducting material. This can be achieved, for example, by using good conductors (such as copper or silver) and/or lowering the temperature. At very low temperatures one could use also superconducting materials to achieve extremely good performance.
- Increase the wire radius a. In typical examples, this action can be limited by physical size considerations. The purpose of this action is mainly to reduce the resistive losses in the wire by increasing the cross-sectional area through which the electric current is flowing, so one could alternatively use also a Litz wire, or ribbon, or any low AC-resistance structure, instead of a circular wire.
- For fixed desired distance D of energy transfer, increase the radius of the loop r. In typical examples, this action can be limited by physical size considerations.
- For fixed desired distance vs. loop-size ratio D/r, decrease the radius of the loop r. In typical examples, this action can be limited by physical size considerations.
- Increase the number of turns N. (Even though Eq. (27) is expected to be less accurate for N>1, qualitatively it still provides a good indication that we expect an improvement in the coupling-to-loss ratio with increased N.) In typical examples, this action can be limited by physical size and possible voltage considerations, as will be discussed in following paragraphs.
- Adjust the alignment and orientation between the two coils. The figure-of-merit is optimized when both cylindrical coils have exactly the same axis of cylindrical symmetry (namely they are “facing” each other). In some examples, particular mutual coil angles and orientations that lead to zero mutual inductance (such as the orientation where the axes of the two coils are perpendicular and the centers of the two coils are on one of the two axes) should be avoided.
- Finally, note that the height of the coil h is another available design parameter, which can have an impact on the performance similar to that of its radius r, and thus the design rules can be similar.

The above analysis technique can be used to design systems with desired parameters. For example, as listed below, the above described techniques can be used to determine the cross sectional radius a of the wire used to design a system including two same single-turn loops with a given radius in order to achieve a specific performance in terms of U=κ/Γ at a given D/r between them, when the loop material is copper (σ=5.998·10⁷S/m):

- D/r=5, U≧10, r=30 cma≧9 mm
- D/r=5, U≧10, r=5 cma≧3.7 mm
- D/r=5, U≧20, r=30 cma≧20 mm
- D/r=5, U≧20, r=5 cma≧8.3 mm
- D/r=10, U≧1, r=30 cma≧7 mm
- D/r=10, U≧1, r=5 cma≧2.8 mm
- D/r=10, U≧3, r=30 cma≧25 mm
- D/r=10, U≧3, r=5 cma≧10 mm

Similar analysis can be done for the case of two dissimilar loops. For example, in some examples, the device under consideration may be identified specifically (e.g. a laptop or a cell phone), so the dimensions of the device object (r_d,h_d,a_d,N_d) may be restricted. However, in some such examples, the restrictions on the source object (r_s,h_s,a_s,N_s) may be much less, since the source can, for example, be placed under the floor or on the ceiling. In such cases, the desired distance between the source and device may be fixed; in other cases it may be variable. Listed below are examples (simplified to the case N_s=N_d=1 and h_s=h_d=0) of how one can vary the dimensions of the source object to achieve a desired system performance in terms of U_sd=κ/√{square root over (Γ_sΓ_d)}, when the material is again copper (σ=5.998·10⁷S/m):

- D=1.5 m, U_sd≧15, r_d=30 cm, a_d=6 mmr_s=1.158 m, a_s≧5 mm
- D=1.5 m, U_sd≧30, r_d=30 cm, a_d=6 mmr_s=1.15 m, a_s≧33 mm
- D=1.5 m, U_sd≧1, r_d=5 cm, a_d=4 mmr_s=1.119 m, a_s≧7 mm
- D=1.5 m, U_sd≧2, r_d=5 cm, a_d=4 mmr_s=1.119 m, a_s≧52 mm
- D=2 m, U^sd≧10, r_d=30 cm, a_d=6 mmr_s=1.518 m, a_s≧7 mm
- D=2 m, U_sd≧20, r_d=30 cm, a_d=6 mmr_s=1.514 m, a_s≧50 mm
- D=2 m, U_sd≧0.5, r_d=5 cm, a_d=4 mmr_s=1.491 m, a_s≧5 mm
- D=2 m, U_sd≧1, r_d=5 cm, a_d=4 mmr_s=1.491 m, a_s≧36 mm
  
  2.2.3 Optimization of k

As described below, in some examples, the quality factor Q of the resonant objects may be limited from external perturbations and thus varying the coil parameters may not lead to significant improvements in Q. In such cases, one can opt to increase the strong-coupling factor U by increasing the coupling factor k. The coupling does not depend on the frequency and may weakly depend on the number of turns. Therefore, in some examples, one can:

- Increase the wire radii a₁and a₂. In typical examples, this action can be limited by physical size considerations.
- For fixed desired distance D of energy transfer, increase the radii of the coils r₁and r₂. In typical examples, this action can be limited by physical size considerations.
- For fixed desired distance vs. coil-sizes ratio D/√{square root over (r₁r₂)}, only the weak (logarithmic) dependence of the inductance remains, which suggests that one should decrease the radii of the coils r₁and r₂. In typical examples, this action can be limited by physical size considerations.
- Adjust the alignment and orientation between the two coils. In typical examples, the coupling is optimized when both cylindrical coils have exactly the same axis of cylindrical symmetry (namely they are “facing” each other). Particular mutual coil angles and orientations that lead to zero mutual inductance (such as the orientation where the axes of the two coils are perpendicular and the centers of the two coils are on one of the two axes) should obviously be avoided.
- Finally, note that the heights of the coils h₁and h₂are other available design parameters, which can have an impact to the coupling similar to that of their radii r₁and r₂, and thus the design rules can be similar.

Further practical considerations apart from efficiency, e.g. physical size limitations, will be discussed in detail below.

2.2.4 Optimization of Overall System Performance

In embodiments, the dimensions of the resonant objects may be determined by the particular application. For example, when the application is powering a laptop or a cell-phone, the device resonant object cannot have dimensions that exceed those of the laptop or cell-phone respectively. For a system of two loops of specified dimensions, in terms of loop radii r_s,dand wire radii a_s,d, the independent parameters left to adjust for the system optimization are: the number of turns N_s,d, the frequency f, the power-load consumption rate κ_l=R_l/2L_dand the power-generator feeding rate κ_g=R_g/2L_s, where R_gis the internal (characteristic) impedance of the generator.

In general, in various examples, the dependent variable that one may want to increase or optimize may be the overall efficiency η. However, other important variables may need to be taken into consideration upon system design. For example, in examples featuring capacitively-loaded coils, the designs can be constrained by, the currents flowing inside the wires I_s,dand other components and the voltages across the capacitors V_s,d. These limitations can be important because for ˜Watt power applications the values for these parameters can be too large for the wires or the capacitors respectively to handle. Furthermore, the total loaded (by the load) quality factor of the device Q_d[l]=ω/2(Γ_d+Γ_l)=ωL_d/(R_d+R_l) and the total loaded (by the generator) quality factor of the source Q_s[g]=ω/2(Γ_s+Γ_g)=ωL_s/(R_s+R_g) are quantities that should be preferably small, because to match the source and device resonant frequencies to within their Q's, when those are very large, can be challenging experimentally and more sensitive to slight variations. Lastly, the radiated powers P_s,radand P_d,radmay need to be minimized for concerns about far-field interference and safety, even though, in general, for a magnetic, non-radiative scheme they are already typically small. In the following, we examine then the effects of each one of the independent variables on the dependent ones.

We define a new variable wp to express the power-load consumption rate for some particular value of U through U_l=κ_l/Γ_d=√{square root over (1+wp·U²)}. Then, in some examples, values which may impact the choice of this rate are: U_l=1 custom character wp=0 to minimize the required energy stored in the source (and therefore I_sand V_s), U_l=√{square root over (1+U²)}>1wp=1 to maximize the efficiency, as seen earlier, or U_l>>1wp>>1 to decrease the required energy stored in the device (and therefore I_dand V_d) and to decrease or minimize Q_d[l]. Similar is the impact of the choice of the power-generator feeding rate U_g=κ_g/Γ_s, with the roles of the source/device and generator/load reversed.

In some examples, increasing N_sand N_dmay increase Q_sand Q_d, and thus U and the efficiency significantly. It also may decrease the currents I_sand I_d, because the inductance of the loops may increase, and thus the energy W_s,d=L_s,d|I_s,d|²/2 required for given output power P_lcan be achieved with smaller currents. However, in some examples, increasing N_dand thus Q_dcan increase Q_d[l], P_d,radand the voltage across the device capacitance V_d. Similar can be the impact of increasing N_son Q_s[g], P_s,radand V_s. As a conclusion, in some examples, the number of turns N_sand N_dmay be chosen large enough (for high efficiency) but such that they allow for reasonable voltages, loaded Q's and/or powers radiated.

With respect to the resonant frequency, again, there may be an optimal one for efficiency. Close to that optimal frequency Q_d[l] and/or Q_s[g] can be approximately maximum. In some examples, for lower frequencies the currents may typically get larger but the voltages and radiated powers may get smaller. In some examples, one may choose the resonant frequency to maximize any of a number of system parameters or performance specifications, such as efficiency.

One way to decide on an operating regime for the system may be based on a graphical method. Consider two loops of r_s=25 cm, r_d=15 cm, h_s=h_d=0, a_s=a_d=3 mm and distance D=2 m between them. In FIG. 12, we plot some of the above dependent variables (currents, voltages and radiated powers normalized to 1 Watt of output power) in terms of frequency f and N_d, given some choice for wp and N_s. FIG. 12 depicts the trend of system performance explained above. In FIG. 13, we make a contour plot of the dependent variables as functions of both frequency and wp but for both N_sand N_dfixed. For example, in embodiments, a reasonable choice of parameters for the system of two loops with the dimensions given above may be: N_s=2, N_d=6, f=10 MHz and wp=10, which gives the following performance characteristics: η=20.6%, Q_d[l]=1264, I_s=7.2 A, I_d=1.4 A, V_s=2.55 kV, V_d=2.30 kV, P_s,rad=0.15 W, P_d,rad=0.006 W. Note that the results in FIGS. 12, 13 and the calculated performance characteristics are made using the analytical formulas provided above, so they are expected to be less accurate for large values of N_s, N_d, but still may give a good estimate of the scalings and the orders of magnitude.

Finally, in embodiments, one could additionally optimize for the source dimensions, if, for example, only the device dimensions are limited. Namely, one can add r_sand a_sin the set of independent variables and optimize with respect to these all the dependent variables of the system. In embodiments, such an optimization may lead to improved results.

2.3 Inductively-Loaded Conducting Rods

A straight conducting rod of length 2h and cross-sectional radius a has distributed capacitance and distributed inductance, and therefore can support a resonant mode of angular frequency ω. Using the same procedure as in the case of self-resonant coils, one can define an effective total inductance L and an effective total capacitance C of the rod through formulas Eqs. (19) and (20). With these definitions, the resonant angular frequency and the effective impedance may be given again by the common formulas ω=1/√{square root over (LC)} and Z=L/C respectively. To calculate the total inductance and capacitance, one can assume again a sinusoidal current profile along the length of the conducting wire. When interested in the lowest mode, if we denote by x the coordinate along the conductor, such that it runs from −h to +h, then the current amplitude profile may have the form I(x)=I_ocos(πx/2h), since it has to be zero at the open ends of the rod. This is the well-known half-wavelength electric dipole resonant mode.

In some examples, one or more of the resonant objects may be inductively-loaded conducting rods. Referring to FIG. 14, a straight conducting rod of length 2h and cross-sectional radius a, as in the previous paragraph, is cut into two equal pieces of length h, which may be connected via a coil wrapped around a magnetic material of relative permeability μ, and everything is surrounded by air. The coil has an inductance L_c, which is added to the distributed inductance of the rod and thus modifies its resonance. Note however, that the presence of the center-loading inductor may modify the current distribution inside the wire and therefore the total effective inductance L and total effective capacitance C of the rod may be different respectively from L_sand C_s, which are calculated for a self-resonant rod of the same total length using a sinusoidal current profile, as in the previous paragraph. Since some current may be running inside the coil of the external loading inductor, the current distribution j inside the rod may be reduced, so L<L_s, and thus, from the charge conservation equation, the linear charge distribution ρ_lmay flatten out towards the center (being positive in one side of the rod and negative in the other side of the rod, changing abruptly through the inductor), so C>C_s. The resonant frequency for this system may be w=1/√{square root over ((L+L_c)C)}<w_s=1/√{square root over (L_sC_s)}, and I(x)→I_ocos(πx/2h) custom character L→L_sω→ω_s, as L_c→0.

In general, the desired CMT parameters can be found for this system, but again a very complicated solution of Maxwell's Equations is generally required. In a special case, a reasonable estimate for the current distribution can be made. When L_c>>L_s>L, then ω≈1/√{square root over (L_cC)}<<ω_sand Z≈√{square root over (L_c/C)}>>Z_s, while the current distribution is triangular along the rod (with maximum at the center-loading inductor and zero at the ends) and thus the charge distribution may be positive constant on one half of the rod and equally negative constant on the other side of the rod. This allows us to compute numerically C from Eq. (20). In this case, the integral in Eq. (20) can actually be computed analytically, giving the formula 1/C=1/(π∈_oh)[ln(h/a)−1]. Explicit analytical formulas are again available for R from Eq. (21) and (22), since I_rms=I_o, |p|=q_oh and |m|=0 (namely only the electric-dipole term is contributing to radiation), so we can determine also Q_abs=1/ωCR_abs, and Q_rad=1/ωCR_rad. At the end of the calculations, the validity of the assumption of triangular current profile may be confirmed by checking that indeed the condition L_c>>L_s custom character ω<<ω_sis satisfied. This condition may be relatively easily satisfied, since typically a conducting rod has very small self-inductance L_sto begin with.

Another important loss factor in this case is the resistive loss inside the coil of the external loading inductor L_cand it may depend on the particular design of the inductor.

In some examples, the inductor may be made of a Brooks coil, which is the coil geometry which, for fixed wire length, may demonstrate the highest inductance and thus quality factor. The Brooks coil geometry has N_Bcturns of conducting wire of cross-sectional radius a_Bcwrapped around a cylindrically symmetric coil former, which forms a coil with a square cross-section of side r_Bc, where the inner side of the square is also at radius r_Bc(and thus the outer side of the square is at radius 2r_Bc), therefore N_Bc≈(r_Bc/2a_Bc)². The inductance of the coil is then L_c=2.0285μ_or_BcN_Bc²≈2.0285μ_or_Bc⁵/8a_Bc⁴and its resistance

$R_{C} \approx \frac{1}{σ} \frac{l_{bc}}{π a_{Bc}^{2}} \sqrt{1 + \frac{μ_{o} ω σ}{2} {(\frac{a_{Bc}}{2})}^{2}},$

where the total wire length is l_Bc≈2≈(3r_Bc/2)N_Bc≈3πr_Bc³/4a_Bc²and we have used an approximate square-root law for the transition of the resistance from the dc to the ac limit as the skin depth varies with frequency.

The external loading inductance L_cprovides the freedom to tune the resonant frequency. For example, for a Brooks coil with a fixed size r_Bc, the resonant frequency can be reduced by increasing the number of turns N_Bcby decreasing the wire cross-sectional radius a_Bc. Then the desired resonant angular frequency ω=1/√{square root over (L_cC)} may be achieved for a_Bc≈(2.0285μ_or_Bc⁵ω²C)^1/4and the resulting coil quality factor may be Q_c≈0.169μ_oσr_Bc²ω/√{square root over (1+ω²μ_oσ√{square root over (2.0285μ_o(r_Bc/4)⁵C)})}. Then, for the particular simple case L_c>>L_s, for which we have analytical formulas, the total Q=1/ωC(R_c+R_abs+R_rad) becomes highest at some optimal frequency w_Q, reaching the value Q_max, that may be determined by the loading-inductor specific design. For example, for the Brooks-coil procedure described above, at the optimal frequency Q_max≈Q_c∞0.8(μ_oσ²r_Bc³/C)^1/4. At lower frequencies it is dominated by ohmic loss inside the inductor coil and at higher frequencies by radiation. Note, again, that the above formulas are accurate as long as ω_Q<<ω_sand, as explained above, this may be easy to design for by using a large inductance.

The results of the above analysis for two examples, using Brooks coils, of subwavelength modes of λ/h≧200 (namely highly suitable for near-field coupling and well within the quasi-static limit) at the optimal frequency w_Qare presented in Table 5.

Table 5 shows in parentheses (for similarity to previous tables) analytical results for the wavelength and absorption, radiation and total loss rates, for two different cases of subwavelength-rod resonant modes. Note that copper was used for the conducting material (σ=5.998·10⁷S/m). The results show that, in some examples, the optimal frequency may be in the low-MHz microwave range and the expected quality factors may be Q_abs>1000 and Q_rad>100000.

TABLE 5

single rod
λ/h
f (MHz)
Q_rad
Q_abs
Q

h = 30 cm, a = 2 cm
(403.8)
(2.477)
(2.72*10⁶)
(7400)
(7380)

μ = 1, r_Bc= 2 cm,

a_Bc= 0.88 mm, N_Bc= 129

h = 10 cm, a = 2 mm
(214.2)
(14.010)
(6.92*10⁵)
(3908)
(3886)

μ = 1, r_Bc= 5 mm,

a_Bc= 0.25 mm, N_Bc= 103

In some examples, energy may be transferred between two inductively-loaded rods. For the rate of energy transfer between two inductively-loaded rods 1 and 2 at distance D between their centers, the mutual capacitance M_Ccan be evaluated numerically from Eq. (23) by using triangular current distributions in the case ω<<ω_s. In this case, the coupling may be only electric and again we have an analytical formula, which, in the quasi-static limit h<<D<<λ and for the relative orientation such that the two rods are aligned on the same axis, is 1/M_C≈1/2π∈_o·(h₁h₂)²/D³, which means that k∝(√{square root over (h₁h₂)}/D)³is independent of the frequency ω. One can then get the resultant strong-coupling factor U.

It can be seen that the frequency ω_U, where the figure-of-merit is maximized to the value U_max, is close to the frequency ω_Q₁_Q₂, where Q₁Q₂is maximized, since k does not depend much on frequency (at least for the distances D<<λ of interest for which the quasi-static approximation is still valid). Therefore, the optimal frequency ω_U≈ω_Q₁_Q₂may be mostly independent of the distance D between the two rods and may lie between the two frequencies ω_Q₁and ω_Q₂at which the single-rod Q₁and Q₂respectively peak. In some typical examples, one can tune the inductively-loaded conducting rods, so that their angular eigenfrequencies are close to ω_Uwithin Γ_U, which is half the angular frequency width for which U>U_max/2.

Referring to Table 6, in parentheses (for similarity to previous tables) analytical results based on the above are shown for two systems each composed of a matched pair of the loaded rods described in Table 5. The average wavelength and loss rates are shown along with the coupling rate and coupling to loss ratio figure-of-merit U=κ/Γ as a function of the coupling distance D, for the two cases. Note that for Γ_radthe single-rod value is used. Again we chose L_c>>L_sto make the triangular-current assumption a good one and computed M_Cnumerically from Eq. (23). The results show that for medium distances D/h=10−3 the expected coupling-to-loss ratios are in the range U˜0.5-100.

TABLE 6

pair of rods
D/h
Q_κ = 1/k
U

h = 30 cm, a = 2 cm
3
(70.3)
(105.0)

μ =1, r_Bc= 2 cm, a_Bc= 0.88 mm, N_Bc= 129
5
(389)
(19.0)

λ/h ≈ 404
7
(1115)
(6.62)

Q ≈ 7380
10
(3321)
(2.22)

h = 10 cm, a = 2 mm
3
(120)
(32.4)

μ =1, r_Bc= 5 mm, a_Bc= 0.25 mm, N_Bc= 103
5
(664)
(5.85)

λ/h ≈ 214
7
(1900)
(2.05)

Q ≈ 3886
10
(5656)
(0.69)

2.4 Dielectric Disks

In some examples, one or more of the resonant objects are dielectric objects, such as disks. Consider a two dimensional dielectric disk object, as shown in FIG. 15(a), of radius r and relative permittivity c surrounded by air that supports high-Q “whispering-gallery” resonant modes. The loss mechanisms for the energy stored inside such a resonant system are radiation into free space and absorption inside the disk material. High-Q_radand long-tailed subwavelength resonances can be achieved when the dielectric permittivity ∈ is large and the azimuthal field variations are slow (namely of small principal number m). Material absorption is related to the material loss tangent: Q_abs□Re{∈}/Im{∈}. Mode-solving calculations for this type of disk resonances were performed using two independent methods: numerically, 2D finite-difference frequency-domain (FDFD) simulations (which solve Maxwell's Equations in frequency domain exactly apart for spatial discretization) were conducted with a resolution of 30 pts/r; analytically, standard separation of variables (SV) in polar coordinates was used.

TABLE 7

single disk
λ/r
Q^abs
Q^rad
Q

Re{ε} = 147.7, m = 2

20.01

10103

1988

1661

(20.00)
(10075)
(1992)
(1663)

Re{ε} = 65.6, m = 3
9.952
10098
9078
4780

(9.950)
(10087)
(9168)
(4802)

The results for two TE-polarized dielectric-disk subwavelength modes of λ/r≧10 are presented in Table 7. Table 7 shows numerical FDFD (and in parentheses analytical SV) results for the wavelength and absorption, radiation and total loss rates, for two different cases of subwavelength-disk resonant modes. Note that disk-material loss-tangent Im{∈}/Re{∈}=10⁻⁴was used. (The specific parameters corresponding to the plot in FIG. 15(a) are highlighted with bold in the table.) The two methods have excellent agreement and imply that for a properly designed resonant low-loss-dielectric object values of Q_rad≧2000 and Q_abs˜10000 are achievable. Note that for the 3D case the computational complexity would be immensely increased, while the physics would not be significantly different. For example, a spherical object of ∈=147.7 has a whispering gallery mode with m=2, Q_rad=13962, and λ/r=17.

The required values of e, shown in Table 7, might at first seem unrealistically large. However, not only are there in the microwave regime (appropriate for approximately meter-range coupling applications) many materials that have both reasonably high enough dielectric constants and low losses (e.g. Titania, Barium tetratitanate, Lithium tantalite etc.), but also ∈ could signify instead the effective index of other known subwavelength surface-wave systems, such as surface modes on surfaces of metallic materials or plasmonic (metal-like, negative-∈) materials or metallo-dielectric photonic crystals or plasmono-dielectric photonic crystals.

To calculate now the achievable rate of energy transfer between two disks 1 and 2, as shown in FIG. 15(b) we place them at distance D between their centers. Numerically, the FDFD mode-solver simulations give κ through the frequency splitting of the normal modes of the combined system (δ_E=2κ from Eq. (4)), which are even and odd superpositions of the initial single-disk modes; analytically, using the expressions for the separation-of-variables eigenfields E_1,2(r) CMT gives κ through

κ=ω₁/2··d³r∈₂(r)E*₂(r)E₁(r)/∫d³r∈(r)|E₁(r)|²,

where ∈_j(r) and ∈(r) are the dielectric functions that describe only the disk j (minus the constant ∈_obackground) and the whole space respectively. Then, for medium distances D/r=10−3 and for non-radiative coupling such that D<2r_c, where r_c=mλ/2π is the radius of the radiation caustic, the two methods agree very well, and we finally find, as shown in Table 8, strong-coupling factors in the range U˜1-50. Thus, for the analyzed examples, the achieved figure-of-merit values are large enough to be useful for typical applications, as discussed below.

TABLE 8

Q =

two disks
D/r
Q^rad
ω/2/Γ
ω/2κ
κ/Γ

Re{ε} = 147.7, m = 2
3
2478
1989
46.9
(47.5)
42.4
(35.0)

λ/r ≈ 20
5
2411
1946
298.0
(298.0)
6.5
(5.6)

Q^abs≈ 10093
7
2196
1804
769.7
(770.2)
2.3
(2.2)

10

2017

1681

1714

(1601)

0.98

(1.04)

Re{ε} = 65.6, m = 3
3
7972
4455
144
(140)
30.9
(34.3)

λ/r ≈ 10
5
9240
4824
2242
(2083)
2.2
(2.3)

Q^abs≈ 10096
7
9187
4810
7485
(7417)
0.64
(0.65)

Note that even though particular examples are presented and analyzed above as examples of systems that use resonant electromagnetic coupling for wireless energy transfer, those of self-resonant conducting coils, capacitively-loaded resonant conducting coils, inductively-loaded resonant conducting rods and resonant dielectric disks, any system that supports an electromagnetic mode with its electromagnetic energy extending much further than its size can be used for transferring energy. For example, there can be many abstract geometries with distributed capacitances and inductances that support the desired kind of resonances. In some examples, the resonant structure can be a dielectric sphere. In any one of these geometries, one can choose certain parameters to increase and/or optimize U or, if the Q's are limited by external factors, to increase and/or optimize for k or, if other system performance parameters are of importance, to optimize those.

Illustrative Example

In one example, consider a case of wireless energy transfer between two identical resonant conducting loops, labeled by L₁and L₂. The loops are capacitively-loaded and couple inductively via their mutual inductance. Let r_Adenote the loops' radii, N_Atheir numbers of turns, and b_Athe radii of the wires making the loops. We also denote by D₁₂the center-to-center separation between the loops.

Resonant objects of this type have two main loss mechanisms: ohmic absorption and far-field radiation. Using the same theoretical method as in previous sections, we find that for r_A=7 cm, b_A=6 mm, and N_A=15 turns, the quality factors for absorption, radiation and total loss are respectively, Q_A,abs=πf/Γ_A,abs=3.19×10⁴, Q_A,rad=πf/Γ_A,rad=2.6×10⁸and Q_A=πf/Γ_A=2.84×10⁴at a resonant frequency f=1.8×10⁷Hz (remember that L₁and L₂are identical and have the same properties). Γ_A,abs, Γ_A,radare respectively the rates of absorptive and radiative loss of L₁and L₂, and the rate of coupling between L₁and L₂is denoted by κ₁₂.

When the loops are in fixed distinct parallel planes separated by D₁₂=1.4 m and have their centers on an axis (C) perpendicular to their planes, as shown in FIG. 17a(Left), the coupling factor for inductive coupling (ignoring retardation effects) is k₁₂≡κ₁₂/πf=7.68×10⁻⁵, independent of time, and thus the strong-coupling factor is U₁₂=k₁₂Q_A=2.18. This configuration of parallel loops corresponds to the largest possible coupling rate κ₁₂at the particular separation D₁₂.

We find that the energy transferred to L₂is maximum at time T_*=κt_*=tan⁻¹(2.18)=1.14 custom character t_*=4775(1/f) from Eq. (5), and constitutes η_E*=29% of the initial total energy from Eq. (6a), as shown in FIG. 17a(Right). The energies radiated and absorbed are respectively η_rad,E(t_*)=7.2% and η_abs,E(t_*)=58.1% of the initial total energy, with 5.8% of the energy remaining in L₁.

We would like to be able to further increase the efficiency of energy transfer between these two resonant objects at their distance D₁₂. In some examples, one can use an intermediate resonator between the source and device resonators, so that energy is transferred more efficiently from the source to the device resonator indirectly through the intermediate resonator.

3. Efficient Energy-Transfer by a Chain of Three ‘Strongly Coupled’ Resonances

FIG. 16 shows a schematic that generally describes one example of the invention, in which energy is transferred wirelessly among three resonant objects. Referring to FIG. 16, energy is transferred, over a distance D_1B, from a resonant source object R₁of characteristic size r₁to a resonant intermediate object R_Bof characteristic size r_B, and then, over an additional distance D_B2, from the resonant intermediate object R_Bto a resonant device object R₂of characteristic size r₂, where D_1B+D_B2=D. As described above, the source object R₁can be supplied power from a power generator that is coupled to the source object R₁. In some examples, the power generator can be wirelessly, e.g., inductively, coupled to the source object R₁. In some examples, the power generator can be connected to the source object R₁by means of a wire or cable. Also, the device object R₂can be connected to a power load that consumes energy transferred to the device object R₂. For example, the device object can be connected to e.g. a resistor, a battery, or other device. All objects are resonant objects. The wireless near-field energy transfer is performed using the field (e.g. the electromagnetic field or acoustic field) of the system of three resonant objects.

Different temporal schemes can be employed, depending on the application, to transfer energy among three resonant objects. Here we will consider a particularly simple but important scheme: a one-time finite-amount energy-transfer scheme

3.1 Finite-Amount Energy-Transfer Efficiency

Let again the source, intermediate and device objects be 1, B, 2 respectively and their resonance modes, which we will use for the energy exchange, have angular frequencies ω_1,B,2, frequency-widths due to intrinsic (absorption, radiation etc.) losses Γ_1,B,2and (generally) vector fields F_1,B,2(r), normalized to unity energy. Once the three resonant objects are brought in proximity, they can interact and an appropriate analytical framework for modeling this resonant interaction is again that of the well-known coupled-mode theory (CMT), which can give a good description of the system for resonances having quality factors of at least, for example, 10. Then, using e^−iωttime dependence, for the chain arrangement shown in FIG. 16, the field amplitudes can be shown to satisfy, to lowest order.

$\begin{matrix} \frac{d}{d t} a_{1} (t) = - i (ω_{1} - i Γ_{1}) a_{1} (t) + i κ_{11} a_{1} (t) + i κ_{1 B} a_{B} (t) \frac{d}{d t} a_{B} (t) = - i (ω_{B} - i Γ_{B}) a_{B} (t) + i κ_{BB} a_{B} (t) + i κ_{B 1} a_{1} (t) + i κ_{B 2} a_{2} (t) \frac{d}{d t} a_{2} (t) = - i (ω_{2} - i Γ_{2}) a_{2} (t) + i κ_{22} a_{2} (t) + i κ_{2 B} a_{B} (t) & (34) \end{matrix}$

where κ_11,BB,22are the shifts in each object's frequency due to the presence of the other, which are a second-order correction and can be absorbed into the eigenfrequencies by setting ω_1,B,2→ω_1,B,2+κ_11,BB,22, and κ_ijare the coupling coefficients, which from the reciprocity requirement of the system must satisfy κ_ij=κ_ji. Note that, in some examples, the direct coupling coefficient κ₁₂between the resonant objects 1 and 2 may be much smaller than the coupling coefficients κ_1Band κ_B2between these two resonant objects with the intermediate object B, implying that the direct energy transfer between 1 and 2 is substantially dominated by the indirect energy transfer through the intermediate object. In some examples, if the coupling rates κ_1Band κ_B2are at least 5 times larger than the direct coupling rate κ₁₂, using an intermediate resonator may lead to an improvement in terms of energy transfer efficiency, as the indirect transfer may dominate the direct transfer. Therefore, in the CMT Eqs. (34) above, the direct coupling coefficient κ₁₂has been ignored, to analyze those particular examples.

The three resonant modes of the combined system are found by substituting [a₁(t) a_B(t),a₂(t)]=[A₁,A_B,A₂]e^−iωt. When the resonators 1 and 2 are at exact resonance ω₁=ω₂=ω_Aand for Γ₁=Γ₂=Γ_A, the resonant modes have complex resonant frequencies

ω_±=ω_AB±√{square root over ((Δω_AB)²+κ_1B²+κ_B2²)} and ω_ds=ω_A−iΓ_A (35a)

where Ω=[(ω_A+ω_B)−i(Γ_A+Γ_B)]/2, Δω_AB=[(ω_A−θ_B)−i(Γ_A−Γ_B)]/2 and whose splitting we denote as {tilde over (δ)}_E≡ω₊−ω₋, and corresponding resonant field amplitudes

$\begin{matrix} {\vec{V}}_{\pm} = {[\begin{matrix} A_{1} \\ A_{B} \\ A_{2} \end{matrix}]}_{\pm} = [\begin{matrix} κ_{1 B} \\ Δ ω_{AB} \mp \sqrt{{(Δ ω_{AB})}^{2} + κ_{1 B}^{2} + κ_{B 2}^{2}} \\ κ_{B 2} \end{matrix}] and {\vec{V}}_{ds} = {[\begin{matrix} A_{1} \\ A_{B} \\ A_{2} \end{matrix}]}_{ds} = [\begin{matrix} - κ_{B 2} \\ 0 \\ κ_{1 B} \end{matrix}] & (35 b) \end{matrix}$

Note that, when all resonators are at exact resonance ω₁=ω₂(=ω_A)=ω_Band for Γ₁=Γ₂(=Γ_A)=Γ_B, we get Δω_AB=0, {tilde over (δ)}_E=2√{square root over (κ_1B²+κ_B2²)}, and then

$\begin{matrix} {\overline{ω}}_{\pm} = ω_{A} \pm \sqrt{κ_{1 B}^{2} + κ_{B 2}^{2}} - i Γ_{A} and {\overline{ω}}_{ds} = ω_{A} - i Γ_{A} & (36 a) \\ {\vec{V}}_{\pm} = {[\begin{matrix} A_{1} \\ A_{B} \\ A_{2} \end{matrix}]}_{\pm} = [\begin{matrix} κ_{1 B} \\ \mp \sqrt{κ_{1 B}^{2} + κ_{B 2}^{2}} \\ κ_{B 2} \end{matrix}] and {\vec{V}}_{ds} = {[\begin{matrix} A_{1} \\ A_{B} \\ A_{2} \end{matrix}]}_{ds} = [\begin{matrix} - κ_{B 2} \\ 0 \\ κ_{1 B} \end{matrix}], & (36 b) \end{matrix}$

namely we get that the resonant modes split to a lower frequency, a higher frequency and a same frequency mode.

Assume now that at time t=0 the source object 1 has finite energy |a₁(0)|², while the intermediate and device objects have |a_B(0)|²=|a₂(0)|²=0. Since the objects are coupled, energy will be transferred from 1 to B and from B to 2. With these initial conditions, Eqs. (34) can be solved, predicting the evolution of the field-amplitudes. The energy-transfer efficiency will be {tilde over (η)}_E≡|a₂(t)|²/|a₁(0)|². The ratio of energy converted to loss due to a specific loss mechanism in resonators 1, B and 2, with respective loss rates Γ_1,loss, Γ_B,lossand Γ_2,losswill be {tilde over (η)}_loss,E=∫₀^tdτ[2Γ_1,loss|a₁(τ)|²+2Γ_B,loss|a_B(τ)|²+2Γ_2,loss|a₂(τ)|²]/|a₁(0)|². At exact resonance ω₁=ω₂(=ωA)=ω_B(an optimal condition) and in the special symmetric case Γ₁=Γ₂=Γ_Aand κ_1B=κ_B2=κ, the field amplitudes are

$\begin{matrix} a_{1} (\tilde{T}) = \frac{1}{2} e^{- i ω_{A} t} e^{- \tilde{T} / \tilde{U}} [\tilde{Δ} \frac{\sin (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T})}{\sqrt{1 - {\tilde{Δ}}^{2}}} + \cos (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T}) + e^{- \tilde{Δ} \tilde{T}}] & (37 a) \\ a_{B} (\tilde{T}) = \frac{1}{2} e^{- i ω_{A} t} e^{- \tilde{T} / \tilde{U}} \frac{\sin (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T})}{\sqrt{1 - {\tilde{Δ}}^{2}}} & (37 b) \\ a_{2} (\tilde{T}) = \frac{1}{2} e^{- i ω_{A} t} e^{- \tilde{T} / \tilde{U}} [\tilde{Δ} \frac{\sin (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T})}{\sqrt{1 - {\tilde{Δ}}^{2}}} + \cos (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T}) + e^{- \tilde{Δ} \tilde{T}}] & (37 c) \end{matrix}$

where {tilde over (T)}≡√{square root over (2)}κt, {tilde over (Δ)}₋₁=2√{square root over (2)}κ/(Γ_B−Γ_A) and Ũ=2√{square root over (2)}κ/(Γ_A+Γ_B).

In some examples, the system designer can adjust the duration of the coupling t at will. In some examples, the duration t can be adjusted to maximize the device energy (and thus efficiency {tilde over (η)}_E). Then, in the special case above, it can be inferred from Eq. (37c) that {tilde over (η)}_Eis maximized for the {tilde over (T)}={tilde over (T)}_*that satisfies

$\begin{matrix} [\tilde{Δ} - \tilde{U} (1 - {\tilde{Δ}}^{2})] \frac{\sin (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T})}{\sqrt{1 - {\tilde{Δ}}^{2}}} + (1 + \tilde{Δ} \tilde{U}) [\cos (\sqrt{1 - {\tilde{Δ}}^{2}} \tilde{T}) - e^{\tilde{Δ} \tilde{T}}] = 0. & (38) \end{matrix}$

In general, this equation may not have an obvious analytical solution, but it does admit a simple solution in the following two cases:

When Γ_A=Γ_B custom character {tilde over (Δ)}=0, Ũ=√{square root over (2)}κ/Γ_B, the energy transfer from resonator 1 to resonator 2 is maximized at

{tilde over (T)}_*({tilde over (Δ)}=0)=2 tan⁻¹Ũ (39)

resulting in an energy-transfer efficiency

$\begin{matrix} {\tilde{η}}_{E} ({\tilde{T}}_{*}, \tilde{Δ} = 0) = {[\frac{{\tilde{U}}^{2}}{1 + {\tilde{U}}^{2}} \exp (- \frac{2 \tan^{- 1} \tilde{U}}{\tilde{U}})]}^{2}, & (40) \end{matrix}$

which has the form of the two-object energy-transfer efficiency of Eq. (6b) but squared.

When Γ_A=0 custom character {tilde over (Δ)}⁻¹=Ũ=2√{square root over (2)}κ/Γ_B, the energy transfer from resonator 1 to resonator 2 is maximized at

$\begin{matrix} {\tilde{T}}_{*} ({\tilde{Δ}}^{- 1} = \tilde{U}) = {\begin{matrix} π \tilde{U} / \sqrt{{\tilde{U}}^{2} - 1}, & \tilde{U} > 1 \\ \infty, & \tilde{U} \leq 1 \end{matrix} & (41) \end{matrix}$

resulting in an energy-transfer efficiency

$\begin{matrix} η_{E} ({\tilde{T}}_{*}, {\tilde{Δ}}^{- 1} = \tilde{U}) = \frac{1}{4} \cdot {\begin{matrix} {[1 + \exp (- π / \sqrt{{\tilde{U}}^{2} - 1})]}^{2}, \tilde{U} > 1 \\ 1, \tilde{U} \leq 1 \end{matrix} . & (42) \end{matrix}$

In both cases, and in general for any {tilde over (Δ)}, the efficiency is an increasing function of Ũ. Therefore, once more resonators that have high quality factors are preferred. In some examples, one may design resonators with Q>50. In some examples, one may design resonators with Q>100.

Illustrative Example

Returning to our illustrative example, in order to improve the ˜29% efficiency of the energy transfer from L₁to L₂, while keeping the distance D₁₂separating them fixed, we propose to introduce an intermediate resonant object that couples strongly to both L₁and L₂, while having the same resonant frequency as both of them. In one example, we take that mediator to also be a capacitively-loaded conducting-wire loop, which we label by L_B. We place L_Bat equal distance (D_1B=D_B2=D₁₂/2=0.7 m) from L₁and L₂such that its axis also lies on the same axis (C), and we orient it such that its plane is parallel to the planes of L₁and L₂, which is the optimal orientation in terms of coupling. A schematic diagram of the three-loops configuration is depicted in FIG. 17b(Left).

In order for L_Bto couple strongly to L₁and L₂, its size needs to be substantially larger than the size of L₁and L₂. However, this increase in the size of L_Bhas considerable drawback in the sense that it may also be accompanied by significant decrease in its radiation quality factor. This feature may often occur for the resonant systems of this type: stronger coupling can often be enabled by increasing the objects' size, but it may imply stronger radiation from the object in question. Large radiation may often be undesirable, because it could lead to far-field interference with other RF systems, and in some systems also because of safety concerns. For r_B=70 cm, b_B=1.5 cm, and N_B=1 turn, we get Q_B,abs=πf/Γ_B,abs=7706, Q_B,rad=πf/Γ_B,rad=400, so Q_B=πf/Γ_B=380, and k_1B≡κ_1B/πf=k_B2=0.0056, so Ũ=2√{square root over (2)}κ/(Γ_A+Γ_B)=5.94 and {tilde over (Δ)}₋₁=2√{square root over (2)}κ/(Γ_B−Γ_A)=6.1, at f=1.8×10⁷Hz. Note that since the coupling rates κ_1Band κ_B2are ≈70 times larger than κ₁₂, indeed we can ignore the direct coupling between L₁and L₂, and focus only on the indirect energy transfer through the intermediate loop L_B, therefore the analysis of the previous section can be used.

The optimum in energy transferred to L₂occurs at time {tilde over (T)}_*=√{square root over (2)}κt_*=3.21 custom character t_*=129(1/f), calculated from Eq. (38), and is equal to {tilde over (η)}_E*=61.5% of the initial energy from Eq. (37c). [Note that, since Q_A>>Q_B, we could have used the analytical conclusions of the case {tilde over (Δ)}⁻¹=Ũ and then we would have gotten a very close approximation of {tilde over (T)}_*=3.19 from Eq. (41).] The energy radiated is {tilde over (η)}_rad,E(t_*)=31.1%, while the energy absorbed is {tilde over (η)}_abs,E(t_*)=3.3%, and 4.1% of the initial energy is left in L₁. Thus, the energy transferred, now indirectly, from L₁to L₂has increased by factor of 2 relative to the two-loops direct transfer case. Furthermore, the transfer time in the three-loops case is now ≈35 times shorter than in the two-loops direct transfer, because of the stronger coupling rate. The dynamics of the energy transfer in the three-loops case is shown in FIG. 17b(Right).

Note that the energy radiated in the three-loop system has undesirably increased by factor of 4 compared to the two-loop system. We would like to be able to achieve a similar improvement in energy-transfer efficiency, while not allowing the radiated energy to increase. In this specification, we disclose that, in some examples, this can be achieved by appropriately varying the values of the coupling strengths between the three resonators.

4. Efficient Energy-Transfer by a Chain of Three Resonances with Adiabatically Varying Coupling Strengths

Consider again the system of a source resonator R₁, a device resonator R₂and an intermediate resonator R_B. For the purposes of the present analysis, R₁and R₂will be assumed to have negligible mutual interactions with each other, while each of them can be strongly coupled to R_B, with coupling rates κ_1Band κ_B2respectively. Note that, in some examples of wireless energy transfer systems, for the resonant object R_Bto be able to have strong coupling with other resonant objects, it may be accompanied with inferior loss properties compared to R₁and R₂, usually in terms of substantially larger radiation losses.

It was seen in a previous section that, when the resonators 1 and 2 are at exact resonance ω₁=ω₂=ω_Aand for Γ₁=Γ₂=Γ_A, the system supports a resonant state (eigenstate) with resonant frequency (eigenfrequency) {tilde over (ω)}_ds=ω_A−iΓ_Aand resonant field amplitudes {right arrow over (V)}_ds=[−κ_B20 κ_1B]^T/√{square root over (κ_1B²+κ_B2²)}, where here we normalized it. This eigenstate {right arrow over (V)}_dswe call the “dark state” in analogy with atomic systems and the related phenomenon of Electromagnetically Induced Transparency (EIT), wherein complete population transfer between two quantum states through a third lossy state, coupled to each of the other two states, is enabled. The dark state is the most essential building block of our proposed efficient weakly-radiative energy-transfer scheme, because it has no energy at all in the intermediate (lossy) resonator R_B, i.e. a_B(t)=0 ∀ t whenever the three-object system is in state {right arrow over (V)}_ds. In fact, if Γ_A→0, then this state is completely lossless, or if Γ_A,rad→0, then this state is completely non-radiative. Therefore, we disclose using predominantly this state to implement the wireless energy transfer. By doing that, the proposed energy transfer scheme can be made completely lossless, in the limit Γ_A→0, no matter how large is the loss rate Γ_B, as shown in FIG. 20, or completely non-radiative, in the limit Γ_A,rad→0, no matter how large is the radiative loss rate Γ_B,rad. Since the energy transfer efficiency increases as the quality factors of the first (source) and second (device) resonances increase, one may design these resonators to have a high quality factor. In some examples, one may design resonators with Q_A>50. In some examples, one may design resonators with Q_A>100.

Let us demonstrate how it is possible to use the dark state for energy transfer with minimal loss. From the expression of {right arrow over (V)}_dsone can see that, if the three-object system is in the state {right arrow over (V)}_ds, then, in general, there is energy in the source resonator with field amplitude proportional to the coupling rate κ_B2between the device resonator and the intermediate resonator, and there is also energy in the device resonator with field amplitude proportional to the coupling rate κ_1Bbetween the source resonator and the intermediate resonator. Then, κ_1B=0 corresponds to all the system's energy being in R₁, while κ_B2=0 corresponds to all the system's energy being in R₂.

So, the important considerations necessary to achieve efficient weakly-radiative energy transfer, consist of preparing the system initially in state {tilde over (V)}_dsand varying the coupling rates in time appropriately to evolve this state {tilde over (V)}_dsin a way that will cause energy transfer. Thus, if at t=0 all the energy is in R₁, then one should have κ_1B(t=0)<<κ_B2(t=0), for example κ_1B(t=0)=0 and κ_B2(t=0)≠0. In order for the total energy of the system to end up in R₂, at a time t_EITwhen the full variation of the coupling rates has been completed, we should have κ_1B(t=t_EIT)>>κ_B2(t=t_EIT), for example κ_1B(t=t_EIT)≠0 and κ_B2(t=t_EIT)=0. This ensures that the initial and final states of the three-object system are parallel to {tilde over (V)}_ds. However, a second important consideration is to keep the three-object system at all times in {tilde over (V)}_ds(t), even as κ_1B(t) and κ_B2(t) are varied in time. This is crucial in order to prevent the system's energy from getting into any of the two other eigenstates {tilde over (V)}_± and thus getting into the intermediate object R_B, which may be highly radiative or lossy in general, as in the example of FIG. 17. This consideration requires changing κ_1B(t) and κ_B2(t) slowly enough so as to make the entire three-object system adiabatically follow the time evolution of {tilde over (V)}_ds(t). The criterion for adiabatic following can be expressed, in analogy to the population transfer case as

$\begin{matrix} \langle 〈 {\vec{V}}_{\pm} ❘ \frac{d {\vec{V}}_{ds}}{d t} 〉 \rangle << \langle {\overline{ω}}_{\pm} - {\overline{ω}}_{ds} \rangle & (43) \end{matrix}$

where {tilde over (V)}_± are the remaining two eigenstates of the system, with corresponding eigenvalues ω_± as shown earlier. Note that any functional dependence of the coupling rates with time (with duration parameter t_EIT) will work, provided it satisfies the adiabaticity criterion Eq. (43) above. The time functional can be linear, sinusoidal (as in the illustrative example to follow) or the temporal analog of a Butterworth (maximally flat) taper, a Chebyshev taper, an exponential taper and the like.

Referring to FIG. 18, an example coupling rate adjustment system 300 for adjusting coupling rates for the one or more of the resonator structures R₁, R₂, or R_Bis shown. As described, the coupling rates between the first and intermediate resonator structures and the intermediate and second resonator structures are characterized by κ_1Band κ_B2respectively. These coupling rates, κ_1Band κ_B2, are several times (e.g., 70 times) greater than the coupling rate κ₁₂between the first and second resonator structure. In some examples, the coupling rate adjustment system can be a mechanical, electrical, or electro-mechanical system for dynamically adjusting, e.g., rotating, or effecting a translational movement, of the one or more resonator structures with respect to each other.

In some examples, the coupling rate κ_1Bis much smaller than the coupling rate the coupling rate κ_B2at the beginning of the energy transfer. By the end, i.e., when a substantial amount of energy has been transferred from the first resonator structure R₁to the second resonator structure, R₂, the coupling rate κ_1Bis much greater than the coupling rate κ_B2. In some examples, the coupling rate κ_1Bcan be set to a fixed value while the coupling rate κ_B2is being varied from its maximum to its minimum value. In some examples, the coupling rate κ_B2can be set to a fixed value while the coupling rate κ_1Bis being varied from its minimum to its maximum value. In some examples, the coupling rate κ_1Bcan be varied from a minimum to a maximum value while the coupling rate κ_B2is being varied from its maximum to minimum value.

Referring now to FIG. 19, a graph for implementing an example coupling rate adjustment system 300 is shown. As shown, in some examples, the coupling rate κ_1Bis set at its minimum value at time, t=0, and increased as a function of time (see, for example, equation 44), while the coupling rate κ_B2is at its maximum value at t=0 and decreased as a function of time (see, for example, equation 45). Accordingly, at the beginning (t=0), the value of κ_1Bis much smaller than the value of κ_B2. In some examples, the value of κ_1Bcan be selected to be any value much smaller than the value of κ_B2. During the wireless energy transfer, the value of κ_1Bis increased, while the value of κ_B2is decreased. After a predetermined amount of time t_EIThas elapsed (e.g., after a substantial amount of energy has been transferred to the second resonator), the value of κ_1Bbecomes much greater than the value of κ_B2.

In some implementations, the coupling rate adjustment system 300 can effect an adjustment of coupling rates between the resonator structures by changing a relative orientation of one or more of the resonator structures with respect to each other. For example, referring again to FIG. 18, the first and second resonator structures, R₁and R₂, can be rotated about their respective axes (e.g., varying angles θ₁and θ₂), with respect to the intermediate resonator structure R_Bto simultaneously change κ_1Band κ_B2. Alternatively, the orientation of the intermediate resonator structure, R_B, can be adjusted, e.g., rotated about an axis, with respect to the first and second resonator structures to simultaneously change κ_1Band κ_B2. Alternatively, the orientation of only the first resonator structure R₁can be rotated about its respective axis to change κ_1B, while the orientations of R₂and R_Bare fixed and thus κ_B2is fixed to a value intermediate between the minimum and maximum values of κ_1B. Alternatively, the orientation of only the second resonator structure R₂can be rotated about its respective axis to change κ_B2, while the orientations of R₁and R_Bare fixed and thus κ_1Bis fixed to a value intermediate between the minimum and maximum values of κ_B2.

In some implementations, the coupling rate adjustment system 300 can effect an adjustment of coupling rates between the resonator structures by translationally moving one or more of the resonator structures with respect to each other. For example, the positions of the first and second resonator structures, R₁and R₂, can be adjusted, e.g., moved along an axis, with respect to the intermediate resonator structure R_Bto simultaneously change κ_1Band κ_B2. Alternatively, a position of the intermediate resonator structure, R_B, can be adjusted, e.g., moved along an axis, with respect to the first and second resonator structures to simultaneously change κ_1Band κ_B2. Alternatively, a position of only the first resonator structure, R₁, can be adjusted, e.g., moved along an axis, with respect to the intermediate R_Bresonator structure to change κ_1B, while the positions of R₂and R_Bare fixed and thus κ_B2is fixed to a value intermediate between the minimum and maximum values of κ_1B. Alternatively, a position of only the second resonator structure, R₂, can be adjusted, e.g., moved along an axis, with respect to the intermediate R_Bresonator structure to change κ_B2, while the positions of R₁and R_Bare fixed and thus κ_1Bis fixed to a value intermediate between the minimum and maximum values of κ_B2.

In some examples, the coupling rate adjustment system 300 can dynamically adjust an effective size of the resonator objects to effect adjustments in the coupling rates κ_1Band κ_B2similar to that described above. The effective size can be adjusted by changing a physical size of the resonator objects. For example, the physical size can be adjusted by effecting mechanical changes in area, length, or other physical aspect of one or more of the resonator structures. The effective size can also be adjusted through non-mechanical changes, such as, but not limited to, applying a magnetic field to change the permeability of the one or more of the resonator structures.

In principle, one would think of making the transfer time t_EITas long as possible to ensure adiabaticity. However there is limitation on how slow the transfer process can optimally be, imposed by the losses in R_land R₂. Such a limitation may not be a strong concern in typical atomic EIT case, because the initial and final states there can be chosen to be non-lossy ground states. However, in our case, losses in R₁and R₂are not avoidable, and can be detrimental to the energy transfer process whenever the transfer time t_EITis not less than 1/Γ_A. This is because, even if the three-object system is carefully kept in {tilde over (V)}_dsat all times, the total energy of the system will decrease from its initial value as a consequence of losses in R_land R₂. Thus the duration of the transfer may be a compromise between these two limits: the desire to keep t_EITlong enough to ensure near-adiabaticity, but short enough not to suffer from losses in R₁and R₂.

Given a particular functional variation of the coupling rates with time with variation duration parameter t_EIT, one may calculate the optimal energy transfer efficiency in the following way: First, for each t_EIT, determine the optimal time t_*, at which the energy at R₂is maximized and the transfer process may be terminated. Then find the optimal variation time t_EIT*based on the counteracting mechanisms discussed above. The optimal efficiency of energy transfer {tilde over (η)}_EIT,E*can then be calculated. In most cases, this procedure may need to be done numerically using the CMT Eqs. (34) as analytical solutions may not be possible. With respect to optimizing the functional dependence of the coupling rates with time, one may choose one that minimizes the coupling of energy to the eigenstates {tilde over (V)}_± for a given t_EIT, which may lead to the temporal analog of a Chebyshev taper.

In some examples, the optimal t_EITmay not be long enough for the adiabadicity criterion of Eq. (43) to be always satisfied. In those cases, some energy may get into at least one of the lossy states {tilde over (V)}_±. Still significant improvement in efficiency and radiation loss may be achieved by the mode of operation where the coupling rates are variable, compared to the mode of operation where the coupling rates are constant, provided the maximum energy that enters the states {tilde over (V)}_± is much less in the variable rate scenario than in the constant rate scenario. In examples, using the proposed scheme of time-varying coupling rates may be advantageous as long as the maximum energy stored in the intermediate resonator is substantially small. In some examples, substantially small may be at most 5% of the peak total energy of the system. In some examples, substantially small may be at most 10% of the peak total energy of the system.

We can now also see why the mode of operation of the system where the coupling rates are kept constant in time may cause a considerable amount of lost (and especially radiated) energy, compared to the proposed mode of operation where the coupling rates are varied adiabatically in time. The reason is that, when κ_1B=κ_B2=const, the energies in R₁and R₂will always be equal to each other if the three-object system is to stay in {tilde over (V)}_ds. So one cannot transfer energy from R₁to R₂by keeping the system purely in state {tilde over (V)}_ds; note that even the initial state of the system, in which all the energy is in R₁and there is no energy in R₃, cannot be solely in {tilde over (V)}_dsfor fixed nonzero κ_1Band κ_B2, and has nonzero components along the eigenstates {tilde over (V)}_± which implies a finite energy will build up in R_B, and consequently result in an increased radiation, especially if Γ_B,rad>>Γ_A,rad, which may be the case if the resonator R_Bis chosen large enough to couple strongly to R₁and R₂, as explained earlier.

Illustrative Example

The previous analysis explains why a considerable amount of energy was radiated when the inductive coupling rates of the loops were kept constant in time, like in FIG. 17b. Let us now consider the modifications necessary for an adiabatically-varied-κ three-loops indirect transfer scheme, as suggested in the previous section, aiming to reduce the total radiated energy back to its reasonable value in the two-loops direct transfer case (FIG. 17a), while maintaining the total energy transfer at level comparable to the constant-K three-loops indirect transfer case (FIG. 17b). In one example, shown in FIG. 17c(Left and Right), we will keep the orientation of L_Bfixed, and start initially (t=0) with L₁perpendicular to L_Bfor κ_1B=0 and L₂parallel to L_Bfor κ_B2=max, then uniformly rotate L₁and L₂, at the same rates, until finally, at (t=t_EIT), L₁becomes parallel to L_Bfor κ_1B=max and L₂perpendicular to L_Bfor κ_B2=0, where we stop the rotation process. In this example, we choose a sinusoidal temporal variation of the coupling rates:

κ_1B(t)=κ sin(πt/2t_EIT) (44)
κ_B2(t)=κ cos(πt/2t_EIT) (45)

for 0<t<t_EIT, and k_B≡κ_1B/πf=k_B2=0.0056 as before. By using the same CMT analysis as in Eq. (34), we find, in FIG. 17c(Center), that for an optimal t_EIT*=1989(1/f), an optimum transfer of {tilde over (η)}_EIT,E*=61.2% can be achieved at t_*=1796(1/f), with only 8.2% of the initial energy being radiated, 28.6% absorbed, and 2% left in L₁. This is quite remarkable: by simply rotating the loops during the transfer, the energy radiated has dropped by factor of 4, while keeping the same 61% level of the energy transferred.

This considerable decrease in radiation may seem surprising, because the intermediate resonator L_B, which mediates all the energy transfer, is highly radiative (650 times more radiative than L₁and L₂), and there is much more time to radiate, since the whole process lasts 14 times longer than in FIG. 17b. Again, the clue to the physical mechanism behind this surprising result can be obtained by observing the differences between the curves describing the energy in R_Bin FIG. 17b and FIG. 17c. Unlike the case of constant coupling rates, depicted in FIG. 17b, where the amount of energy ultimately transferred to L₂goes first through the intermediate loop L_B, with peak energy storage in L_Bas much as 30% of the peak total energy of the system, in the case of time-varying coupling rates, shown in FIG. 17c, there is almost little or no energy in L_Bat all times during the transfer. In other words, the energy is transferred quite efficiently from L₁to L₂, mediated by L_Bwithout considerable amount of energy ever being in the highly radiative intermediate loop L_B. (Note that direct transfer from L₁to L₂is identically zero here since L₁is always perpendicular to L₂, so all the energy transfer is indeed mediated through L_B). In some examples, improvement in efficiency and/or radiated energy can still have been accomplished if the energy transfer had been designed with a time t_EITsmaller than its optimal value (perhaps to speed up the process), if the maximum energy accumulated inside the intermediate resonator was less than 30%. For example, improvement can have been achieved for maximum energy accumulation inside the intermediate resonator of 5% or even 10%.

An example implementation of the coupling rate adjustment system 300 is described below, where the resonators are capacitively-loaded loops, which couple to each other inductively. At the beginning (t=0), the coupling rate adjustment system 300 sets the relative orientation of the first resonator structure L₁to be perpendicular to the intermediate resonator structure L_B. At this orientation, the value of the coupling rate κ_1Bbetween the first and intermediate resonator structure is at its minimum value. Also, the coupling rate adjustment system 300 can set the relative orientation of the second resonator structure L₂to be parallel to the intermediate resonator structure L_B. At this orientation, the value of the coupling rate κ_B2is at a maximum value. During wireless energy transfer, the coupling rate adjustment system 300 can effect the rotation of the first resonator structure L₁about its axis so that the value of κ_1Bis increased. In some examples, the coupling rate adjustment system 300 can also effect the rotation of the second resonator structure, L₂, about its axis so that the value of κ_B2is decreased. In some examples, a similar effect can be achieved by fixing L₁and L₂to be perpendicular to each other and rotating only L_Bto be parallel to L₂and perpendicular to L₁at t=0 and parallel to L₁and perpendicular to L₂at t=t_EIT. In some examples, a similar effect can be achieved by fixing L_Band one of L₁and L₂(e.g., L₁) at a predetermined orientation (e.g. at 45 degrees with respect to the intermediate resonator L_B) and rotating only the other of L₁and L₂(e.g., L₂from parallel to L_Bat t=0 to perpendicular to L_Bat t=t_EIT).

Similarly, in some implementations, at the beginning (t=0), the coupling rate adjustment system 300 can set the position of the first resonator structure L₁at a first large predetermined distance from the intermediate resonator structure L_Bso that the value of the coupling rate κ_1Bis at its minimum value. Correspondingly, the coupling rate adjustment system 300 can set the position of the second resonator structure L₂at a second small predetermined distance from the intermediate resonator structure L_Bso that the value of the coupling rate κ_B2between the first and intermediate resonator structure is at its maximum value. During wireless energy transfer, the coupling rate adjustment system 300 can affect the position of the first resonator structure L₁to bring it closer to L_Bso that the value of κ_1Bis increased. In some examples, the coupling rate adjustment system 300 can also effect the position of the second resonator structure, L₂, to take it away from L_Bso that the value of κ_B2is decreased. In some examples, a similar effect can be achieved by fixing L₁and L₂to be at a fixed distance to each other and effecting the position of only L_Bto be close to L₂and away from L₁at t=0 and close to L₁and away from L₂at t=t_EIT. In some examples, a similar effect can be achieved by fixing L_Band one of L₁and L₂(e.g., L₁) at a predetermined (not too large but not too small) distance and effecting the position only the other of L₁and L₂(e.g., L₂from close to L_Bat t=0 to away from L_Bat t=t_EIT).

5. Comparison of Static and Adiabatically Dynamic Systems

In the abstract case of energy transfer from R₁to R₂, where no constraints are imposed on the relative magnitude of κ, Γ_rad^A, Γ_rad^B, Γ_abs^A, and Γ_abs^B, it is not certain that the adiabatic-κ (EIT-like) system will always perform better than the constant-κ one, in terms of the transferred and radiated energies. In fact, there could exist some range of the parameters (κ, Γ_rad^A, Γ_rad^B, Γ_abs^A, Γ_abs^B), for which the energy radiated in the constant-K transfer case is less than that radiated in the EIT-like case. For this reason, we investigate both the adiabatic-κ and constant-κ transfer schemes, as we vary some of the crucial parameters of the system. The percentage of energies transferred and lost (radiated+absorbed) depends only on the relative values of κ, Γ_A=Γ_rad^A+Γ_abs^Aand Γ_B=Γ_rad^B+Γ_rad^B. Hence we first calculate and visualize the dependence of these energies on the relevant parameters κ/Γ_Band Γ_B/Γ_A, in the contour plots shown in FIG. 21.

The way the contour plots are calculated is as follows. For each value of (κ/Γ_B, Γ_B/Γ_A) in the adiabatic case, where κ_1B(t) and κ_B2(t) are given by Eq. (44)-(45), one tries range of values of t_EIT. For each t_EIT, the maximum energy transferred E₂(%) over 0<t <t_EIT, denoted by max(E₂, t_EIT), is calculated together with the total energy lost at that maximum transfer. Next the maximum of max(E₂, t_EIT) over all values of t_EITis selected and plotted as single point on the contour plot in FIG. 21a. We refer to this point as the optimum energy transfer (%) in the adiabatic-κ case for the particular (κ/Γ_B, Γ_B/Γ_A) under consideration. We also plot in FIG. 21d the corresponding value of the total energy lost (%) at the optimum of E₂. We repeat these calculations for all pairs (κ/Γ_B, Γ_B/Γ_A) shown in the contour plots. In the constant-κ transfer case, for each (κ/Γ_B, Γ_B/Γ_A), the time evolution of E₂(%) and E_lostare calculated for 0<t<2/κ, and optimum transfer, shown in FIG. 21b, refers to the maximum of E₃(t) over 0<t<2/κ. The corresponding total energy lost at optimum constant-transfer is shown in FIG. 21e. Now that we calculated the energies of interest as functions of (κ/Γ_B, Γ_B/Γ_A), we look for ranges of the relevant parameters in which the adiabatic-κ transfer has advantages over the constant-κ one. So, we plot the ratio of (E₂)_adiabatic-κ/(E₂)_constant-κ in FIG. 21c, and (E_lost)_constant-κ/(E_lost)_adiabatic-κ in FIG. 21f. We find that, for Γ_B/Γ_A>50, the optimum energy transferred in the adiabatic-κ case exceeds that in the constant-κ case, and the improvement factor can be larger than 2. From FIG. 21f, one sees that the adiabatic-κ scheme can reduce the total energy lost by factor of 3 compared to the constant-κ scheme, also in the range Γ_B/Γ_A>50. As in the constant-κ case, also in the adiabatic-κ case the efficiency increases as the ratio of the maximum value, κ, of the coupling rates to the loss rate of the intermediate object (and thus also the first and second objects for Γ_B/Γ_A>1) increases. In some examples, one may design a system so that κ is larger than each of Γ_Band Γ_A. In some examples, one may design a system so that κ is at least 2 times larger than each of Γ_Band Γ_A. In some examples, one may design a system so that κ is at least 4 times larger than each of Γ_Band Γ_A.

Although one may be interested in reducing the total energy lost (radiated+absorbed) as much as possible in order to make the transfer more efficient, the undersirable nature of the radiated energy may make it important to consider reducing the energy radiated. For this purpose, we calculate the energy radiated at optimum transfer in both the adiabatic-κ and constant-κ schemes, and compare them. The relevant parameters in this case are κ/Γ_B, Γ_B/Γ_A, Γ_rad^A/Γ_A, and Γ_rad^B/Γ_B. The problem is more complex because the parameter space is now 4-dimensional. So we focus on those particular cross sections that can best reveal the most important differences between the two schemes. From FIGS. 21c and 21f, one can guess that the best improvement in both E₂and E_lostoccurs for Γ_B/Γ_A≧2500. Moreover, knowing that it is the intermediate object R_Bthat makes the main difference between the adiabatic-κ and constant-κ schemes, being “energy-empty” in the adiabatic-κ case and “energy-full” in the constant-κ one, we first look at the special situation where Γ_rad^A=0. In FIG. 22a and FIG. 22b, we show contour plots of the energy radiated at optimum transfer, in the constant-κ and adiabatic-K schemes respectively, for the particular cross section having Γ_B/Γ_A=500 and Γ_rad^A=0. Comparing these two figures, one can see that, by using the adiabatic-κ scheme, one can reduce the energy radiated by factor of 6.3 or more.

To get quantitative estimate of the radiation reduction factor in the general case where Γ_A,rad≠0, we calculate the ratio of energies radiated at optimum transfers in both schemes, namely,

$\begin{matrix} \frac{{(E_{rad})}_{constant - κ}}{{(E_{rad})}_{adiabatic - κ}} = \frac{2 \int_{0}^{t_{*}^{constant - κ}} {\begin{matrix} \frac{Γ_{rad}^{B}}{Γ_{rad}^{A}} {\langle a_{B}^{constant - κ} (t) \rangle}^{2} + \\ {\langle a_{1}^{constant - κ} (t) \rangle}^{2} + {\langle a_{2}^{adiabatic - κ} (t) \rangle}^{2} \end{matrix}}}{2 \int_{0}^{t_{*}^{adiabatic - κ}} {\begin{matrix} \frac{Γ_{rad}^{B}}{Γ_{rad}^{A}} {\langle a_{B}^{adiabatic - κ} (t) \rangle}^{2} + \\ {\langle a_{1}^{adiabatic - κ} (t) \rangle}^{2} + {\langle a_{2}^{adiabatic - κ} (t) \rangle}^{2} \end{matrix}}} & (46) \end{matrix}$

which depends only on Γ_rad^B/Γ_rad^A, the time-dependent mode amplitudes, and the optimum transfer times in both schemes. The latter two quantities are completely determined by κ/Γ_Band Γ_B/Γ_A. Hence the only parameters relevant to the calculations of the ratio of radiated energies are Γ_rad^B/Γ_rad^A, κ/Γ_Band Γ_B/Γ_A, thus reducing the dimensionality of the investigated parameter space from down to 3. For convenience, we multiply the first relevant parameter Γ_rad^B/Γ_rad^Aby Γ_B/Γ_Awhich becomes (Γ_rad^B/Γ_B)/(Γ_rad^A/Γ_A), i.e. the ratio of quantities that specify what percentage of each object's loss is radiated. Next, we calculate the ratio of energies radiated as function of (Γ_rad^B/Γ_B)/(_rad^A/Γ_A) and κ/Γ_Bin the two special cases Γ_B/Γ_A=50, and Γ_B/Γ_A=500, and plot them in FIG. 22c and FIG. 22d, respectively. We also show, in FIG. 22e, the dependence of (E_rad)_constant-κ/(E_rad)_EIT-likeon κ/Γ_Band Γ_B/Γ_A, for the special case Γ_rad^A=0. As can be seen from FIG. 22c-22d, the adiabatic-κ scheme is less radiative than the constant-κ scheme whenever Γ_rad^B/Γ_Bis larger than Γ_rad^A/Γ_A, and the radiation reduction ratio increases as Γ_B/Γ_Aand κ/Γ_Bare increased (see FIG. 22e). In some examples, the adiabatic-κ scheme is less radiative than the constant-κ scheme whenever Γ_rad^B/Γ_rad^Ais larger than about 20. In some examples, the adiabatic-κ scheme is less radiative than the constant-κ scheme whenever Γ_rad^B/Γ_rad^Ais larger than about 50.

It is to be understood that while three resonant objects are shown in the previous examples, other examples can feature four or more resonant objects. For example, in some examples, a single source object can transfer energy to multiple device objects through one intermediate object. In some examples, energy can be transferred from a source resonant object to a device resonant object, through two or more intermediate resonant objects, and so forth.

6. System Sensitivity to Extraneous Objects

In general, the overall performance of an example of the resonance-based wireless energy-transfer scheme depends strongly on the robustness of the resonant objects' resonances. Therefore, it is desirable to analyze the resonant objects' sensitivity to the near presence of random non-resonant extraneous objects. One appropriate analytical model is that of “perturbation theory” (PT), which suggests that in the presence of an extraneous perturbing object p the field amplitude a₁(t) inside the resonant object 1 satisfies, to first order:

$\begin{matrix} \frac{d a_{1}}{d t} = - i (ω_{1} - i Γ_{1}) a_{1} + i (δ κ_{11 (p)} + i δ Γ_{1 (p)}) a_{1} & (47) \end{matrix}$

where again ω₁is the frequency and Γ₁the intrinsic (absorption, radiation etc.) loss rate, while δκ_11(p)is the frequency shift induced onto 1 due to the presence of p and δΓ_1(p)is the extrinsic due to p (absorption inside p, scattering from p etc.) loss rate. δΓ_1(p)is defined as δΓ_1(p)≡Γ_1(p)−Γ₁, where Γ_1(p)is the total perturbed loss rate in the presence of p. The first-order PT model is valid only for small perturbations. Nevertheless, the parameters δκ_11(p), δΓ_1(p)are well defined, even outside that regime, if a₁is taken to be the amplitude of the exact perturbed mode. Note also that interference effects between the radiation field of the initial resonant-object mode and the field scattered off the extraneous object can for strong scattering (e.g. off metallic objects) result in total Γ_1,rad(p)that are smaller than the initial Γ_1,rad(namely δΓ_1,rad(p)is negative).

It has been shown that a specific relation is desired between the resonant frequencies of the source and device-objects and the driving frequency. In some examples, all resonant objects must have the same eigenfrequency and this must be equal to the driving frequency. In some implementations, this frequency-shift can be “fixed” by applying to one or more resonant objects and the driving generator a feedback mechanism that corrects their frequencies. In some examples, the driving frequency from the generator can be fixed and only the resonant frequencies of the objects can be tuned with respect to this driving frequency.

The resonant frequency of an object can be tuned by, for example, adjusting the geometric properties of the object (e.g. the height of a self-resonant coil, the capacitor plate spacing of a capacitively-loaded loop or coil, the dimensions of the inductor of an inductively-loaded rod, the shape of a dielectric disc, etc.) or changing the position of a non-resonant object in the vicinity of the resonant object.

In some examples, referring to FIG. 23a, each resonant object is provided with an oscillator at fixed frequency and a monitor which determines the eigenfrequency of the object. At least one of the oscillator and the monitor is coupled to a frequency adjuster which can adjust the frequency of the resonant object. The frequency adjuster determines the difference between the fixed driving frequency and the object frequency and acts, as described above, to bring the object frequency into the required relation with respect to the fixed frequency. This technique assures that all resonant objects operate at the same fixed frequency, even in the presence of extraneous objects.

In some examples, referring to FIG. 23(b), during energy transfer from a source object to a device object, the device object provides energy or power to a load, and an efficiency monitor measures the efficiency of the energy-transfer or power-transmission. A frequency adjuster coupled to the load and the efficiency monitor acts, as described above, to adjust the frequency of the object to maximize the efficiency.

In other examples, the frequency adjusting scheme can rely on information exchange between the resonant objects. For example, the frequency of a source object can be monitored and transmitted to a device object, which is in turn synched to this frequency using frequency adjusters, as described above. In other embodiments the frequency of a single clock can be transmitted to multiple devices, and each device then synched to that frequency using frequency adjusters, as described above.

Unlike the frequency shift, the extrinsic perturbing loss due to the presence of extraneous perturbing objects can be detrimental to the functionality of the energy-transfer scheme, because it is difficult to remedy. Therefore, the total perturbed quality factors Q_(p)(and the corresponding perturbed strong-coupling factor U_(p)should be quantified.

In some examples, a system for wireless energy-transfer uses primarily magnetic resonances, wherein the energy stored in the near field in the air region surrounding the resonator is predominantly magnetic, while the electric energy is stored primarily inside the resonator. Such resonances can exist in the quasi-static regime of operation (r<<λ) that we are considering: for example, for coils with h<<2r, most of the electric field is localized within the self-capacitance of the coil or the externally loading capacitor and, for dielectric disks, with ∈>>1 the electric field is preferentially localized inside the disk. In some examples, the influence of extraneous objects on magnetic resonances is nearly absent. The reason is that extraneous non-conducting objects p that could interact with the magnetic field in the air region surrounding the resonator and act as a perturbation to the resonance are those having significant magnetic properties (magnetic permeability Re{μ}>1 or magnetic loss Im{μ}>0). Since almost all every-day non-conducting materials are non-magnetic but just dielectric, they respond to magnetic fields in the same way as free space, and thus will not disturb the resonance of the resonator. Extraneous conducting materials can however lead to some extrinsic losses due to the eddy currents induced inside them or on their surface (depending on their conductivity). However, even for such conducting materials, their presence will not be detrimental to the resonances, as long as they are not in very close proximity to the resonant objects.

The interaction between extraneous objects and resonant objects is reciprocal, namely, if an extraneous object does not influence a resonant object, then also the resonant object does not influence the extraneous object. This fact can be viewed in light of safety considerations for human beings. Humans are also non-magnetic and can sustain strong magnetic fields without undergoing any risk. A typical example, where magnetic fields B˜1T are safely used on humans, is the Magnetic Resonance Imaging (MRI) technique for medical testing. In contrast, the magnetic near-field required in typical embodiments in order to provide a few Watts of power to devices is only B˜10⁻⁴T, which is actually comparable to the magnitude of the Earth's magnetic field. Since, as explained above, a strong electric near-field is also not present and the radiation produced from this non-radiative scheme is minimal, the energy-transfer apparatus, methods and systems described herein is believed safe for living organisms.

An advantage of the presently proposed technique using adiabatic variations of the coupling rates between the first and intermediate resonators and between the intermediate and second resonators compared to a mode of operation where these coupling rates are not varied but are constant is that the interactions of the intermediate resonator with extraneous objects can be greatly reduced with the presently proposed scheme. The reason is once more the fact that there is always a substantially small amount of energy in the intermediate resonator in the adiabatic-κ scheme, therefore there is little energy that can be induced from the intermediate object to an extraneous object in its vicinity. Furthermore, since the losses of the intermediate resonator are substantially avoided in the adiabatic-κ case, the system is less immune to potential reductions of the quality factor of the intermediate resonator due to extraneous objects in its vicinity.

6.1 Capacitively-Loaded Conducting Loops or Coils

In some examples, one can estimate the degree to which the resonant system of a capacitively-loaded conducting-wire coil has mostly magnetic energy stored in the space surrounding it. If one ignores the fringing electric field from the capacitor, the electric and magnetic energy densities in the space surrounding the coil come just from the electric and magnetic field produced by the current in the wire; note that in the far field, these two energy densities must be equal, as is always the case for radiative fields. By using the results for the fields produced by a subwavelength (r<<λ) current loop (magnetic dipole) with h=0, we can calculate the ratio of electric to magnetic energy densities, as a function of distance D_pfrom the center of the loop (in the limit r<<D_p) and the angle θ with respect to the loop axis:

$\begin{matrix} \begin{matrix} \frac{w_{e} (x)}{w_{m} (x)} = \frac{ɛ_{o} {\langle E (x) \rangle}^{2}}{μ_{o} {\langle H (x) \rangle}^{2}} \\ = \frac{(1 + \frac{1}{x^{2}}) \sin^{2} θ}{(\frac{1}{x^{2}} + \frac{1}{x^{4}}) 4 \cos^{2} θ + (1 - \frac{1}{x^{2}} + \frac{1}{x^{4}}) \sin^{2} θ}; \end{matrix} x = 2 π \frac{D_{p}}{λ} \Rightarrow \frac{\underset{S_{p}}{∯} w_{e} (x) ⅆ S}{\underset{S_{p}}{∯} w_{m} (x) ⅆ S} = \frac{1 + \frac{1}{x^{2}}}{1 + \frac{1}{x^{2}} + \frac{3}{x^{4}}}; x = 2 π \frac{D_{p}}{λ}, & (48) \end{matrix}$

where the second line is the ratio of averages over all angles by integrating the electric and magnetic energy densities over the surface of a sphere of radius D_p. From Eq. (48) it is obvious that indeed for all angles in the near field (x<<1) the magnetic energy density is dominant, while in the far field (x>>1) they are equal as they should be. Also, the preferred positioning of the loop is such that objects which can interfere with its resonance lie close to its axis (θ=0), where there is no electric field. For example, using the systems described in Table 4, we can estimate from Eq. (48) that for the loop of r=30 cm at a distance D_p=10r=3 m the ratio of average electric to average magnetic energy density would be ˜12% and at D_p=3r=90 cm it would be ˜1%, and for the loop of r=10 cm at a distance D_p=10r=1 m the ratio would be ˜33% and at D_p=3r=30 cm it would be ˜2.5%. At closer distances this ratio is even smaller and thus the energy is predominantly magnetic in the near field, while in the radiative far field, where they are necessarily of the same order (ratio→1), both are very small, because the fields have significantly decayed, as capacitively-loaded coil systems are designed to radiate very little. Therefore, this is the criterion that qualifies this class of resonant system as a magnetic resonant system.

To provide an estimate of the effect of extraneous objects on the resonance of a capacitively-loaded loop including the capacitor fringing electric field, we use the perturbation theory formula, stated earlier, δΓ_1,abs(p)=₁/4·∫d³r Im{∈_p(r)}|E₁(r)|²/W with the computational FEFD results for the field of an example like the one shown in the plot of FIG. 5 and with a rectangular object of dimensions 30 cm×30 cm×1.5 m and permittivity ∈=49+16i (consistent with human muscles) residing between the loops and almost standing on top of one capacitor (˜3 cm away from it) and find δQ_abs(human)˜10⁵and for ˜10 cm away δQ_abs(human)˜5·10⁵. Thus, for ordinary distances (˜1 m) and placements (not immediately on top of the capacitor) or for most ordinary extraneous objects p of much smaller loss-tangent, we conclude that it is indeed fair to say that δQ_abs(p)→∞. The only perturbation that is expected to affect these resonances is a close proximity of large metallic structures.

Self-resonant coils can be more sensitive than capacitively-loaded coils, since for the former the electric field extends over a much larger region in space (the entire coil) rather than for the latter (just inside the capacitor). On the other hand, self-resonant coils can be simple to make and can withstand much larger voltages than most lumped capacitors. Inductively-loaded conducting rods can also be more sensitive than capacitively-loaded coils, since they rely on the electric field to achieve the coupling.

6.2 Dielectric Disks

For dielectric disks, small, low-index, low-material-loss or far-away stray objects will induce small scattering and absorption. In such cases of small perturbations these extrinsic loss mechanisms can be quantified using respectively the analytical first-order perturbation theory formulas

[δQ_1,rad(p)]⁻¹≡2δΓ_1,rad(p)/ω₁∝∫d³r[Re{∈_p(r)}|E₁(r)|]²/w
[δQ_1,abs(p)]⁻¹≡2δΓ_1,rad(p)/ω₁=∫d³rIm{∈_p(r)}|E₁(r)|²/2w

where W=∫d³r∈(r)|E₁(r)|²/2 is the total resonant electromagnetic energy of the unperturbed mode. As one can see, both of these losses depend on the square of the resonant electric field tails E1 at the site of the extraneous object. In contrast, the coupling factor from object 1 to another resonant object 2 is, as stated earlier,

k₁₂=2κ₁₂/√{square root over (ω₁ω₂)}≈∫d³r∈₂(r)E*₂(r)E₁(r)/∫d³r∈(r)|E₁(r)|²

and depends linearly on the field tails E₁of 1 inside 2. This difference in scaling gives us confidence that, for, for example, exponentially small field tails, coupling to other resonant objects should be much faster than all extrinsic loss rates (κ₁₂>>δΓ_1,2(p)), at least for small perturbations, and thus the energy-transfer scheme is expected to be sturdy for this class of resonant dielectric disks.

However, we also want to examine certain possible situations where extraneous objects cause perturbations too strong to analyze using the above first-order perturbation theory approach. For example, we place a dielectric disk close to another off-resonance object of large Re{∈}, Im{∈} and of same size but different shape (such as a human being h), as shown in FIG. 24a, and a roughened surface of large extent but of small Re{∈}, Im{∈} (such as a wall w), as shown in FIG. 24b. For distances D_h,w/r=10−3 between the disk-center and the “human”-center or “wall”, the numerical FDFD simulation results presented in FIGS. 24a and 24b suggest that, the disk resonance seems to be fairly robust, since it is not detrimentally disturbed by the presence of extraneous objects, with the exception of the very close proximity of high-loss objects. To examine the influence of large perturbations on an entire energy-transfer system we consider two resonant disks in the close presence of both a “human” and a “wall”. Comparing Table 8 to the table in FIG. 24c, the numerical FDFD simulations show that the system performance deteriorates from U˜1-50 to U_(hw)˜0.5-10, i.e. only by acceptably small amounts.

In general, different examples of resonant systems have different degree of sensitivity to external perturbations, and the resonant system of choice depends on the particular application at hand, and how important matters of sensitivity or safety are for that application. For example, for a medical implantable device (such as a wirelessly powered artificial heart) the electric field extent must be minimized to the highest degree possible to protect the tissue surrounding the device. In such cases where sensitivity to external objects or safety is important, one should design the resonant systems so that the ratio of electric to magnetic energy density w_e/w_mis reduced or minimized at most of the desired (according to the application) points in the surrounding space.

7 Applications

The non-radiative wireless energy transfer techniques described above can enable efficient wireless energy-exchange between resonant objects, while suffering only modest transfer and dissipation of energy into other extraneous off-resonant objects. The technique is general, and can be applied to a variety of resonant systems in nature. In this Section, we identify a variety of applications that can benefit from or be designed to utilize wireless power transmission.

Remote devices can be powered directly, using the wirelessly supplied power or energy to operate or run the devices, or the devices can be powered by or through or in addition to a battery or energy storage unit, where the battery is occasionally being charged or re-charged wirelessly. The devices can be powered by hybrid battery/energy storage devices such as batteries with integrated storage capacitors and the like. Furthermore, novel battery and energy storage devices can be designed to take advantage of the operational improvements enabled by wireless power transmission systems.

Devices can be turned off and the wirelessly supplied power or energy used to charge or recharge a battery or energy storage unit. The battery or energy storage unit charging or recharging rate can be high or low. The battery or energy storage units can be trickle charged or float charged. It would be understood by one of ordinary skill in the art that there are a variety of ways to power and/or charge devices, and the variety of ways could be applied to the list of applications that follows.

Some wireless energy transfer examples that can have a variety of possible applications include for example, placing a source (e.g. one connected to the wired electricity network) on the ceiling of a room, while devices such as robots, vehicles, computers, PDAs or similar are placed or move freely within the room. Other applications can include powering or recharging electric-engine buses and/or hybrid cars and medical implantable devices. Additional example applications include the ability to power or recharge autonomous electronics (e.g. laptops, cell-phones, portable music players, house-hold robots, GPS navigation systems, displays, etc), sensors, industrial and manufacturing equipment, medical devices and monitors, home appliances (e.g. lights, fans, heaters, displays, televisions, counter-top appliances, etc.), military devices, heated or illuminated clothing, communications and navigation equipment, including equipment built into vehicles, clothing and protective-wear such as helmets, body armor and vests, and the like, and the ability to transmit power to physically isolated devices such as to implanted medical devices, to hidden, buried, implanted or embedded sensors or tags, to and/or from roof-top solar panels to indoor distribution panels, and the like.

A number of examples of the invention have been described. Nevertheless, it will be understood that various modifications can be made without departing from the spirit and scope of the invention.

Claims

1. A method comprising: transferring energy from a first resonant object to a second resonant object using a third resonant object to mediate the energy transfer, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, more strongly than the first resonant object is capable of resonantly coupling with the second resonant object; andmaintaining a reduced level of energy accumulation in the third resonant object by varying over time a coupling rate between the third resonant object and (i) the first resonant object, (ii) the second resonant object, or (iii) both, based on a loss rate of (i) the first resonant object, (ii) the second resonant object, or (iii) both;wherein the reduced level of energy accumulation is reduced relative to energy accumulation that would occur in the third resonant object without the varying of the coupling rate over time.
2. The method of claim 1, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, at least five times more strongly than the first resonant object is capable of resonantly coupling with the second resonant object.
3. The method of claim 1, wherein varying the coupling rate comprises changing a relative orientation between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
4. The method of claim 1, wherein varying the coupling rate comprises changing a relative position between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
5. The method of claim 1, wherein varying the coupling rate comprises changing a size of at least one of the first resonant object, the second resonant object, or the third resonant object.
6. A method comprising: transferring energy from a first resonant object to a second resonant object using a third resonant object to mediate the energy transfer, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, more strongly than the first resonant object is capable of resonantly coupling with the second resonant object; andmaintaining a reduced level of energy accumulation in the third resonant object by varying over time a coupling rate between the third resonant object and (i) the first resonant object, (ii) the second resonant object, or (iii) both, based on a loss rate of (i) the first resonant object, (ii) the second resonant object, or (iii) both;wherein:transferring the energy comprises wirelessly transferring energy from the first resonant object to the second resonant object using the third resonant object to mediate the wireless energy transfer; andmaintaining the reduced level of energy accumulation in the third resonant object comprises varying both (i) a coupling rate between the first resonant object and the third resonant object and (ii) a coupling rate between the third resonant object and the second resonant object, over time in accordance with an adiabatic time evolution of an eigenstate for the three resonant objects.
7. The method of claim 6, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, at least five times more strongly than the first resonant object is capable of resonantly coupling with the second resonant object.
8. The method of claim 6, wherein the first resonant object has a loss rate unequal to that of the second resonant object.
9. The method of claim 6, wherein each of the first resonant object, the second resonant object, and the third resonant object has a resonant frequency between 50 kHz and 500 MHz.
10. The method of claim 6, wherein varying the coupling rate comprises changing a relative orientation between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
11. The method of claim 6, wherein varying the coupling rate comprises changing a relative position between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
12. The method of claim 6, wherein varying the coupling rate comprises changing a resonator property of at least one of the first resonant object, the second resonant object, or the third resonant object.
13. The method of claim 6, wherein varying the coupling rate comprises changing a size of at least one of the first resonant object, the second resonant object, or the third resonant object.
14. A system comprising: a first resonant object to provide energy;a second resonant object to receive energy;a third resonant object to mediate energy transfer between the first resonant object and the second resonant object, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, more strongly than the first resonant object is capable of resonantly coupling with the second resonant object; anda control mechanism configured and arranged to maintain a reduced level of energy accumulation in the third resonant object by varying over time a coupling rate between the third resonant object and (i) the first resonant object, (ii) the second resonant object, or (iii) both, based on a loss rate of (i) the first resonant object, (ii) the second resonant object, or (iii) both;wherein the reduced level of energy accumulation is reduced relative to energy accumulation that would occur in the third resonant object without the varying of the coupling rate over time.
15. The system of claim 14, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, at least five times more strongly than the first resonant object is capable of resonantly coupling with the second resonant object.
16. The system of claim 14, wherein the control mechanism is configured and arranged to change a relative orientation between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
17. The system of claim 14, wherein the control mechanism is configured and arranged to change a relative position between one or more pairs of the first resonant object, the second resonant object, and the third resonant object.
18. The system of claim 14, wherein the control mechanism is configured and arranged to change a size of at least one of the first resonant object, the second resonant object, or the third resonant object.
19. A system comprising: a first resonant object to provide energy;a second resonant object to receive energy;a third resonant object to mediate energy transfer between the first resonant object and the second resonant object, wherein the third resonant object is capable of resonantly coupling with (i) the first resonant object, (ii) the second resonant object, or (iii) both, more strongly than the first resonant object is capable of resonantly coupling with the second resonant object; anda control mechanism configured and arranged to improve energy transfer by varying over time a coupling rate between the third resonant object and (i) the first resonant object, (ii) the second resonant object, or (iii) both, based on a loss rate of (i) the first resonant object, (ii) the second resonant object, or (iii) both;wherein the first resonant object is a first wireless resonator, the second resonant object is a second wireless resonator, the third resonant object is a third wireless resonator, and the control mechanism is configured and arranged to vary both (i) a coupling rate between the first wireless resonator and the third wireless resonator and (ii) a coupling rate between the third wireless resonator and the second wireless resonator, over time in accordance with an adiabatic time evolution of an eigenstate for the three wireless resonators.
20. The system of claim 19, wherein the third wireless resonator is capable of resonantly coupling with (i) the first wireless resonator, (ii) the second wireless resonator, or (iii) both, at least five times more strongly than the first wireless resonator is capable of resonantly coupling with the second wireless resonator.
21. The system of claim 19, wherein the first wireless resonator has a loss rate unequal to that of the second wireless resonator.
22. The system of claim 19, wherein each of the first wireless resonator, the second wireless resonator, and the third wireless resonator has a resonant frequency between 50 kHz and 500 MHz.
23. The system of claim 19, wherein the control mechanism is configured and arranged to change a relative orientation between one or more pairs of the first wireless resonator, the second wireless resonator, and the third wireless resonator.
24. The system of claim 19, wherein the control mechanism is configured and arranged to change a relative position between one or more pairs of the first wireless resonator, the second wireless resonator, and the third wireless resonator.
25. The system of claim 19, wherein the control mechanism is configured and arranged to change a resonator property of at least one of the first wireless resonator, the second wireless resonator, or the third wireless resonator.
26. The system of claim 19, wherein the control mechanism is configured and arranged to change a size of at least one of the first wireless resonator, the second wireless resonator, or the third wireless resonator.
27. An apparatus comprising: an intermediary wireless resonator to mediate energy transfer between a source wireless resonator and a destination wireless resonator, wherein the intermediary wireless resonator is capable of resonantly coupling with (i) the source wireless resonator, (ii) the destination wireless resonator, or (iii) both, more strongly than the source wireless resonator is capable of resonantly coupling with the destination wireless resonator; anda control mechanism configured and arranged to maintain a reduced level of energy accumulation in the intermediary wireless resonator by varying over time a coupling rate between the source wireless resonator and the intermediary wireless resonator and between the intermediary wireless resonator and the destination wireless resonator;wherein the reduced level of energy accumulation is reduced relative to energy accumulation that would occur in the intermediary wireless resonator without the varying of the coupling rate over time.
28. The apparatus of claim 27, wherein the control mechanism is configured and arranged to change an orientation, a resonator property, or a size of the intermediary wireless resonator.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation, and claims benefit, of U.S. application Ser. No. 13/677,407, filed on Nov. 15, 2012, which is a continuation of, and claims benefit, of U.S. application Ser. No. 12/571,949, filed on Oct. 1, 2009, now U.S. Pat. No. 8,362,651 issued on Jan. 29, 2013, which claims priority to U.S. Provisional Application Ser. No. 61/101,809, filed Oct. 1, 2008. The contents of the prior applications are incorporated herein by reference in their entirety.

Government Interests

This invention was made with government support under grant number DE-FG02-99ER45778 awarded by the Department of Energy, grant number DMR-0213282 awarded by the National Science Foundation, and grant number W911NF-07-D-0004 awarded by the Army Research Office. The Government has certain rights in this invention.

US Referenced Citations (485)

Number	Name	Date	Kind
645576	Tesla	Mar 1900	A
649621	Tesla	May 1900	A
787412	Tesla	Apr 1905	A
1119732	Tesla	Dec 1914	A
2133494	Waters	Oct 1938	A
2266262	Polydoroff	Dec 1941	A
3517350	Beaver	Jun 1970	A
3535543	Dailey	Oct 1970	A
3780425	Penn et al.	Dec 1973	A
3871176	Schukei	Mar 1975	A
4088999	Fletcher et al.	May 1978	A
4095998	Hanson	Jun 1978	A
4180795	Matsuda et al.	Dec 1979	A
4280129	Wells	Jul 1981	A
4450431	Hochstein	May 1984	A
4588978	Allen	May 1986	A
4621243	Harada	Nov 1986	A
5027709	Slagle	Jul 1991	A
5033295	Schmid et al.	Jul 1991	A
5034658	Hiering et al.	Jul 1991	A
5053774	Schuermann et al.	Oct 1991	A
5070293	Ishii et al.	Dec 1991	A
5118997	El-Hamamsy	Jun 1992	A
5216402	Carosa	Jun 1993	A
5229652	Hough	Jul 1993	A
5287112	Schuermann	Feb 1994	A
5293308	Boys et al.	Mar 1994	A
5341083	Klontz et al.	Aug 1994	A
5367242	Hulman	Nov 1994	A
5374930	Schuermann	Dec 1994	A
5408209	Tanzer et al.	Apr 1995	A
5437057	Richley et al.	Jul 1995	A
5455467	Young et al.	Oct 1995	A
5493691	Barrett	Feb 1996	A
5522856	Reineman	Jun 1996	A
5528113	Boys et al.	Jun 1996	A
5541604	Meier	Jul 1996	A
5550452	Shirai et al.	Aug 1996	A
5565763	Arrendale et al.	Oct 1996	A
5630835	Brownlee	May 1997	A
5631660	Higashiguchi et al.	May 1997	A
5697956	Bornzin	Dec 1997	A
5703461	Minoshima et al.	Dec 1997	A
5703573	Fujimoto et al.	Dec 1997	A
5710413	King et al.	Jan 1998	A
5742471	Barbee, Jr. et al.	Apr 1998	A
5821728	Schwind	Oct 1998	A
5821731	Kuki et al.	Oct 1998	A
5864323	Berthon	Jan 1999	A
5898579	Boys et al.	Apr 1999	A
5903134	Takeuchi	May 1999	A
5923544	Urano	Jul 1999	A
5940509	Jovanovich et al.	Aug 1999	A
5957956	Kroll et al.	Sep 1999	A
5959245	Moe et al.	Sep 1999	A
5986895	Stewart et al.	Nov 1999	A
5993996	Firsich	Nov 1999	A
5999308	Nelson et al.	Dec 1999	A
6012659	Nakazawa et al.	Jan 2000	A
6047214	Mueller et al.	Apr 2000	A
6066163	John	May 2000	A
6067473	Greeninger et al.	May 2000	A
6108579	Snell et al.	Aug 2000	A
6127799	Krishnan	Oct 2000	A
6130591	Tsuzuki	Oct 2000	A
6176433	Uesaka et al.	Jan 2001	B1
6184651	Fernandez et al.	Feb 2001	B1
6207887	Bass et al.	Mar 2001	B1
6225800	Zhang et al.	May 2001	B1
6232841	Bartlett et al.	May 2001	B1
6238387	Miller, III	May 2001	B1
6252762	Amatucci	Jun 2001	B1
6262639	Shu et al.	Jul 2001	B1
6300760	Schubert et al.	Oct 2001	B1
6407470	Seelig	Jun 2002	B1
6436299	Baarman et al.	Aug 2002	B1
6450946	Forsell	Sep 2002	B1
6452465	Brown et al.	Sep 2002	B1
6459218	Boys et al.	Oct 2002	B2
6473028	Luc	Oct 2002	B1
6483202	Boys	Nov 2002	B1
6515878	Meins et al.	Feb 2003	B1
6533178	Gaul et al.	Mar 2003	B1
6535133	Gohara	Mar 2003	B2
6561975	Pool et al.	May 2003	B1
6563425	Nicholson et al.	May 2003	B2
6597076	Scheible et al.	Jul 2003	B2
6609023	Fischell et al.	Aug 2003	B1
6631072	Paul et al.	Oct 2003	B1
6650227	Bradin	Nov 2003	B1
6664770	Bartels	Dec 2003	B1
6673250	Kuennen et al.	Jan 2004	B2
6683256	Kao	Jan 2004	B2
6696647	Ono et al.	Feb 2004	B2
6703921	Wuidart et al.	Mar 2004	B1
6731071	Baarman	May 2004	B2
6749119	Scheible et al.	Jun 2004	B2
6772011	Dolgin	Aug 2004	B2
6798716	Charych	Sep 2004	B1
6803744	Sabo	Oct 2004	B1
6806649	Mollema et al.	Oct 2004	B2
6812645	Baaman	Nov 2004	B2
6825620	Kuennen et al.	Nov 2004	B2
6831417	Baarman	Dec 2004	B2
6839035	Addonisio et al.	Jan 2005	B1
6844702	Giannopoulos et al.	Jan 2005	B2
6856291	Mickle et al.	Feb 2005	B2
6858970	Malkin et al.	Feb 2005	B2
6906495	Cheng et al.	Jun 2005	B2
6917163	Baarman	Jul 2005	B2
6917431	Soljacic et al.	Jul 2005	B2
6937130	Scheible et al.	Aug 2005	B2
6960968	Odendaal et al.	Nov 2005	B2
6961619	Casey	Nov 2005	B2
6967462	Landis	Nov 2005	B1
6975198	Baarman	Dec 2005	B2
6988026	Breed et al.	Jan 2006	B2
7027311	Vanderelli et al.	Apr 2006	B2
7035076	Stevenson	Apr 2006	B1
7042196	Ka-Lai et al.	May 2006	B2
7069064	Govorgian et al.	Jun 2006	B2
7084605	Mickle et al.	Aug 2006	B2
7116200	Baarman et al.	Oct 2006	B2
7118240	Baarman et al.	Oct 2006	B2
7126450	Baarman et al.	Oct 2006	B2
7127293	MacDonald	Oct 2006	B2
7132918	Baarman et al.	Nov 2006	B2
7147604	Allen et al.	Dec 2006	B1
7180248	Kuennen et al.	Feb 2007	B2
7191007	Desai et al.	Mar 2007	B2
7193418	Freytag	Mar 2007	B2
7212414	Baarman	May 2007	B2
7233137	Nakamura et al.	Jun 2007	B2
7239110	Cheng et al.	Jul 2007	B2
7248017	Cheng et al.	Jul 2007	B2
7251527	Lyden	Jul 2007	B2
7288918	DiStefano	Oct 2007	B2
7375492	Calhoon et al.	May 2008	B2
7375493	Calhoon et al.	May 2008	B2
7378817	Calhoon et al.	May 2008	B2
7382636	Baarman et al.	Jun 2008	B2
7385357	Kuennen et al.	Jun 2008	B2
7443135	Cho	Oct 2008	B2
7462951	Baarman	Dec 2008	B1
7466213	Lobl et al.	Dec 2008	B2
7474058	Baarman	Jan 2009	B2
7492247	Schmidt et al.	Feb 2009	B2
7514818	Abe et al.	Apr 2009	B2
7518267	Baarman	Apr 2009	B2
7525283	Cheng et al.	Apr 2009	B2
7545337	Guenther	Jun 2009	B2
7599743	Hassler, Jr. et al.	Oct 2009	B2
7615936	Baarman et al.	Nov 2009	B2
7639514	Baarman	Dec 2009	B2
7741734	Joannopoulos et al.	Jun 2010	B2
7795708	Katti	Sep 2010	B2
7825543	Karalis et al.	Nov 2010	B2
7843288	Lee et al.	Nov 2010	B2
7863859	Soar	Jan 2011	B2
8022576	Joannopoulos et al.	Sep 2011	B2
8076800	Joannopoulos et al.	Dec 2011	B2
8076801	Karalis et al.	Dec 2011	B2
8084889	Joannopoulos et al.	Dec 2011	B2
8097983	Karalis et al.	Jan 2012	B2
8131378	Greenberg et al.	Mar 2012	B2
8178995	Amano et al.	May 2012	B2
8362651	Hamam et al.	Jan 2013	B2
8395282	Joannopoulos et al.	Mar 2013	B2
8395283	Joannopoulos et al.	Mar 2013	B2
8400018	Joannopoulos et al.	Mar 2013	B2
8400019	Joannopoulos et al.	Mar 2013	B2
8400020	Joannopoulos et al.	Mar 2013	B2
8400021	Joannopoulos et al.	Mar 2013	B2
8400022	Joannopoulos et al.	Mar 2013	B2
8400023	Joannopoulos et al.	Mar 2013	B2
8400024	Joannopoulos et al.	Mar 2013	B2
8760007	Joannopoulos et al.	Jun 2014	B2
8760008	Joannopoulos et al.	Jun 2014	B2
8766485	Joannopoulos et al.	Jul 2014	B2
8772971	Joannopoulos et al.	Jul 2014	B2
8772972	Joannopoulos et al.	Jul 2014	B2
8791599	Joannopoulos et al.	Jul 2014	B2
8836172	Hamam et al.	Sep 2014	B2
9065286	Joannopoulos et al.	Jun 2015	B2
20020032471	Loftin et al.	Mar 2002	A1
20020105343	Scheible et al.	Aug 2002	A1
20020118004	Scheible et al.	Aug 2002	A1
20020130642	Ettes et al.	Sep 2002	A1
20020167294	Odaohhara	Nov 2002	A1
20030038641	Scheible	Feb 2003	A1
20030062794	Scheible et al.	Apr 2003	A1
20030062980	Scheible et al.	Apr 2003	A1
20030071034	Thompson et al.	Apr 2003	A1
20030124050	Yadav et al.	Jul 2003	A1
20030126948	Yadav et al.	Jul 2003	A1
20030160590	Schaefer et al.	Aug 2003	A1
20030199778	Mickle et al.	Oct 2003	A1
20030214255	Bauman et al.	Nov 2003	A1
20040000974	Odenaal et al.	Jan 2004	A1
20040026998	Henriott et al.	Feb 2004	A1
20040100338	Clark	May 2004	A1
20040113847	Qi et al.	Jun 2004	A1
20040130425	Dayan et al.	Jul 2004	A1
20040130915	Baarman	Jul 2004	A1
20040130916	Baarman	Jul 2004	A1
20040142733	Parise	Jul 2004	A1
20040150934	Baarman	Aug 2004	A1
20040189246	Bulai et al.	Sep 2004	A1
20040201361	Koh et al.	Oct 2004	A1
20040222751	Mollema et al.	Nov 2004	A1
20040227057	Tuominen et al.	Nov 2004	A1
20040232845	Baarman	Nov 2004	A1
20040233043	Yazawa et al.	Nov 2004	A1
20040267501	Freed et al.	Dec 2004	A1
20050007067	Baarman et al.	Jan 2005	A1
20050021134	Opie	Jan 2005	A1
20050027192	Govari et al.	Feb 2005	A1
20050030251	Okamura et al.	Feb 2005	A1
20050033382	Single	Feb 2005	A1
20050085873	Gord et al.	Apr 2005	A1
20050093475	Kuennen et al.	May 2005	A1
20050104064	Hegarty et al.	May 2005	A1
20050104453	Vanderelli et al.	May 2005	A1
20050116650	Baarman	Jun 2005	A1
20050116683	Cheng et al.	Jun 2005	A1
20050122058	Baarman et al.	Jun 2005	A1
20050122059	Baarman et al.	Jun 2005	A1
20050125093	Kikuchi et al.	Jun 2005	A1
20050127849	Baarman et al.	Jun 2005	A1
20050127850	Baarman et al.	Jun 2005	A1
20050127866	Hamilton et al.	Jun 2005	A1
20050135122	Cheng et al.	Jun 2005	A1
20050140482	Cheng et al.	Jun 2005	A1
20050151511	Chary	Jul 2005	A1
20050156560	Shimaoka et al.	Jul 2005	A1
20050189945	Reiderman	Sep 2005	A1
20050194926	Di Stefano	Sep 2005	A1
20050253152	Klimov et al.	Nov 2005	A1
20050288739	Hassler, Jr. et al.	Dec 2005	A1
20050288740	Hassler, Jr. et al.	Dec 2005	A1
20050288741	Hassler, Jr. et al.	Dec 2005	A1
20050288742	Giordano et al.	Dec 2005	A1
20060001509	Gibbs	Jan 2006	A1
20060022636	Xian et al.	Feb 2006	A1
20060044188	Tsai et al.	Mar 2006	A1
20060053296	Busboom et al.	Mar 2006	A1
20060061323	Cheng et al.	Mar 2006	A1
20060066443	Hall	Mar 2006	A1
20060090956	Peshkovskiy et al.	May 2006	A1
20060132045	Baarman	Jun 2006	A1
20060159392	Popovic	Jul 2006	A1
20060164866	Vanderelli et al.	Jul 2006	A1
20060164868	Weber	Jul 2006	A1
20060181242	Freed et al.	Aug 2006	A1
20060184209	John et al.	Aug 2006	A1
20060184210	Singhal et al.	Aug 2006	A1
20060185809	Elfrink et al.	Aug 2006	A1
20060199620	Greene et al.	Sep 2006	A1
20060202665	Hsu	Sep 2006	A1
20060205381	Beart et al.	Sep 2006	A1
20060214626	Nilson et al.	Sep 2006	A1
20060219448	Grieve et al.	Oct 2006	A1
20060238365	Vecchione et al.	Oct 2006	A1
20060270440	Shearer et al.	Nov 2006	A1
20060281435	Shearer et al.	Dec 2006	A1
20060284708	Reeves	Dec 2006	A1
20070010295	Greene et al.	Jan 2007	A1
20070013483	Stewart	Jan 2007	A1
20070016089	Fischell et al.	Jan 2007	A1
20070021140	Keyes, IV et al.	Jan 2007	A1
20070024246	Flaugher	Feb 2007	A1
20070064406	Beart	Mar 2007	A1
20070069687	Suzuki	Mar 2007	A1
20070096875	Waterhouse et al.	May 2007	A1
20070105429	Kohl et al.	May 2007	A1
20070117596	Greene et al.	May 2007	A1
20070126650	Guenther	Jun 2007	A1
20070145830	Lee et al.	Jun 2007	A1
20070171681	Baarman	Jul 2007	A1
20070176840	Pristas et al.	Aug 2007	A1
20070178945	Cook et al.	Aug 2007	A1
20070182367	Partovi	Aug 2007	A1
20070208263	John et al.	Sep 2007	A1
20070222542	Joannopoulos et al.	Sep 2007	A1
20070247005	Tetlow	Oct 2007	A1
20070267918	Gyland	Nov 2007	A1
20070276538	Kjellsson et al.	Nov 2007	A1
20080012569	Hall et al.	Jan 2008	A1
20080014897	Cook et al.	Jan 2008	A1
20080030415	Homan et al.	Feb 2008	A1
20080067874	Tseng	Mar 2008	A1
20080191638	Kuennen et al.	Aug 2008	A1
20080197710	Kreitz et al.	Aug 2008	A1
20080211320	Cook et al.	Sep 2008	A1
20080265684	Farkas	Oct 2008	A1
20080266748	Lee	Oct 2008	A1
20080278264	Karalis et al.	Nov 2008	A1
20080294208	Willis et al.	Nov 2008	A1
20090010028	Baarmen et al.	Jan 2009	A1
20090015075	Cook et al.	Jan 2009	A1
20090033564	Cook et al.	Feb 2009	A1
20090045772	Cook et al.	Feb 2009	A1
20090051224	Cook et al.	Feb 2009	A1
20090058189	Cook et al.	Mar 2009	A1
20090067198	Graham et al.	Mar 2009	A1
20090072627	Cook et al.	Mar 2009	A1
20090072628	Cook et al.	Mar 2009	A1
20090072629	Cook et al.	Mar 2009	A1
20090079268	Cook et al.	Mar 2009	A1
20090085408	Bruhn	Apr 2009	A1
20090085706	Baaman et al.	Apr 2009	A1
20090096413	Patovi et al.	Apr 2009	A1
20090102292	Cook et al.	Apr 2009	A1
20090108679	Porwal	Apr 2009	A1
20090108997	Patterson et al.	Apr 2009	A1
20090127937	Widmer et al.	May 2009	A1
20090134712	Cook et al.	May 2009	A1
20090146892	Shimizu et al.	Jun 2009	A1
20090153273	Chen	Jun 2009	A1
20090160261	Elo	Jun 2009	A1
20090167449	Cook et al.	Jul 2009	A1
20090174263	Baarman et al.	Jul 2009	A1
20090179502	Cook et al.	Jul 2009	A1
20090189458	Kawasaki	Jul 2009	A1
20090195332	Joannopoulos et al.	Aug 2009	A1
20090195333	Joannopoulos et al.	Aug 2009	A1
20090212636	Cook et al.	Aug 2009	A1
20090213028	Cook et al.	Aug 2009	A1
20090224608	Cook et al.	Sep 2009	A1
20090224609	Cook et al.	Sep 2009	A1
20090224856	Karalis et al.	Sep 2009	A1
20090230777	Baarman et al.	Sep 2009	A1
20090237194	Waffenschmidt et al.	Sep 2009	A1
20090243394	Levine	Oct 2009	A1
20090243397	Cook et al.	Oct 2009	A1
20090251008	Sugaya	Oct 2009	A1
20090267558	Jung	Oct 2009	A1
20090267709	Joannopoulos et al.	Oct 2009	A1
20090267710	Joannopoulos et al.	Oct 2009	A1
20090271047	Wakamatsu	Oct 2009	A1
20090271048	Wakamatsu	Oct 2009	A1
20090273242	Cook	Nov 2009	A1
20090281678	Wakamatsu	Nov 2009	A1
20090284082	Mohammadian	Nov 2009	A1
20090284083	Karalis et al.	Nov 2009	A1
20090284218	Mohammadian et al.	Nov 2009	A1
20090284220	Toncich et al.	Nov 2009	A1
20090284227	Mohammadian et al.	Nov 2009	A1
20090284245	Kirby et al.	Nov 2009	A1
20090284369	Toncich et al.	Nov 2009	A1
20090286470	Mohammadian et al.	Nov 2009	A1
20090286475	Toncich et al.	Nov 2009	A1
20090286476	Toncich et al.	Nov 2009	A1
20090289595	Chen et al.	Nov 2009	A1
20090299918	Cook et al.	Dec 2009	A1
20100017249	Fincham et al.	Jan 2010	A1
20100033021	Bennett	Feb 2010	A1
20100034238	Bennett	Feb 2010	A1
20100036773	Bennett	Feb 2010	A1
20100038970	Cook et al.	Feb 2010	A1
20100045114	Sample et al.	Feb 2010	A1
20100052431	Mita	Mar 2010	A1
20100052811	Smith et al.	Mar 2010	A1
20100065352	Ichikawa	Mar 2010	A1
20100066349	Lin et al.	Mar 2010	A1
20100081379	Cooper et al.	Apr 2010	A1
20100094381	Kim et al.	Apr 2010	A1
20100096934	Joannopoulos et al.	Apr 2010	A1
20100102639	Joannopoulos et al.	Apr 2010	A1
20100102640	Joannopoulos et al.	Apr 2010	A1
20100102641	Joannopoulos et al.	Apr 2010	A1
20100109443	Cook et al.	May 2010	A1
20100109604	Boys et al.	May 2010	A1
20100115474	Takada et al.	May 2010	A1
20100117454	Cook et al.	May 2010	A1
20100117455	Joannopoulos et al.	May 2010	A1
20100117456	Karalis et al.	May 2010	A1
20100117596	Cook et al.	May 2010	A1
20100123353	Joannopoulos et al.	May 2010	A1
20100123354	Joannopoulos et al.	May 2010	A1
20100123355	Joannopoulos et al.	May 2010	A1
20100123452	Amano et al.	May 2010	A1
20100123530	Park et al.	May 2010	A1
20100127573	Joannopoulos et al.	May 2010	A1
20100127574	Joannopoulos et al.	May 2010	A1
20100127575	Joannopoulos et al.	May 2010	A1
20100127660	Cook et al.	May 2010	A1
20100133918	Joannopoulos et al.	Jun 2010	A1
20100133919	Joannopoulos et al.	Jun 2010	A1
20100133920	Joannopoulos et al.	Jun 2010	A1
20100148589	Hamam et al.	Jun 2010	A1
20100148723	Cook et al.	Jun 2010	A1
20100151808	Toncich et al.	Jun 2010	A1
20100156346	Takada et al.	Jun 2010	A1
20100156570	Hong et al.	Jun 2010	A1
20100164295	Ichikawa et al.	Jul 2010	A1
20100171370	Karalis et al.	Jul 2010	A1
20100181844	Karalis et al.	Jul 2010	A1
20100181961	Novak et al.	Jul 2010	A1
20100184371	Cook et al.	Jul 2010	A1
20100187911	Joannopoulos et al.	Jul 2010	A1
20100187913	Sample	Jul 2010	A1
20100190435	Cook et al.	Jul 2010	A1
20100190436	Cook et al.	Jul 2010	A1
20100194206	Burdo et al.	Aug 2010	A1
20100194207	Graham	Aug 2010	A1
20100194334	Kirby et al.	Aug 2010	A1
20100194335	Kirby et al.	Aug 2010	A1
20100201189	Kirby et al.	Aug 2010	A1
20100201201	Mobarhan et al.	Aug 2010	A1
20100201202	Kirby et al.	Aug 2010	A1
20100201204	Sakoda et al.	Aug 2010	A1
20100201205	Karalis et al.	Aug 2010	A1
20100201310	Vorenkamp et al.	Aug 2010	A1
20100201313	Vorenkamp et al.	Aug 2010	A1
20100201316	Takada et al.	Aug 2010	A1
20100201513	Vorenkamp et al.	Aug 2010	A1
20100207458	Joannopoulos et al.	Aug 2010	A1
20100210233	Cook et al.	Aug 2010	A1
20100213770	Kikuchi	Aug 2010	A1
20100213895	Keating et al.	Aug 2010	A1
20100217553	Von Novak et al.	Aug 2010	A1
20100219695	Komiyama et al.	Sep 2010	A1
20100219696	Kojima	Sep 2010	A1
20100222010	Ozaki et al.	Sep 2010	A1
20100225175	Karalis et al.	Sep 2010	A1
20100225270	Jacobs et al.	Sep 2010	A1
20100225271	Oyobe et al.	Sep 2010	A1
20100225272	Kirby et al.	Sep 2010	A1
20100231053	Karalis et al.	Sep 2010	A1
20100231163	Mashinsky	Sep 2010	A1
20100237706	Karalis et al.	Sep 2010	A1
20100237707	Karalis et al.	Sep 2010	A1
20100237708	Karalis et al.	Sep 2010	A1
20100244576	Hillan et al.	Sep 2010	A1
20100244577	Shimokawa	Sep 2010	A1
20100244578	Yoshikawa	Sep 2010	A1
20100244579	Sogabe et al.	Sep 2010	A1
20100244580	Uchida et al.	Sep 2010	A1
20100244581	Uchida	Sep 2010	A1
20100244582	Yoshikawa	Sep 2010	A1
20100244583	Shimokawa	Sep 2010	A1
20100244839	Yoshikawa	Sep 2010	A1
20100248622	Kirby et al.	Sep 2010	A1
20100253152	Karalis et al.	Oct 2010	A1
20100253281	Li	Oct 2010	A1
20100256831	Abramo et al.	Oct 2010	A1
20100259109	Sato	Oct 2010	A1
20100264745	Karalis et al.	Oct 2010	A1
20100264746	Kazama et al.	Oct 2010	A1
20100276995	Marzetta et al.	Nov 2010	A1
20100277003	Von Novak et al.	Nov 2010	A1
20100277004	Suzuki et al.	Nov 2010	A1
20100277005	Karalis et al.	Nov 2010	A1
20100277120	Cook et al.	Nov 2010	A1
20100289341	Ozaki et al.	Nov 2010	A1
20100289449	Elo	Nov 2010	A1
20100295505	Jung et al.	Nov 2010	A1
20100295506	Ichikawa	Nov 2010	A1
20100327660	Karalis et al.	Dec 2010	A1
20100327661	Karalis et al.	Dec 2010	A1
20110012431	Karalis et al.	Jan 2011	A1
20110018361	Karalis et al.	Jan 2011	A1
20110025131	Karalis et al.	Feb 2011	A1
20110043046	Joannopoulos et al.	Feb 2011	A1
20110049996	Karalis et al.	Mar 2011	A1
20110049998	Karalis et al.	Mar 2011	A1
20110074218	Karalis et al.	Mar 2011	A1
20110074347	Karalis et al.	Mar 2011	A1
20110089895	Karalis et al.	Apr 2011	A1
20110140544	Karalis et al.	Jun 2011	A1
20110148219	Karalis et al.	Jun 2011	A1
20110162895	Karalis et al.	Jul 2011	A1
20110169339	Karalis et al.	Jul 2011	A1
20110181122	Karalis et al.	Jul 2011	A1
20110193419	Karalis et al.	Aug 2011	A1
20110198939	Karalis et al.	Aug 2011	A1
20110221278	Karalis et al.	Sep 2011	A1
20110227528	Karalis et al.	Sep 2011	A1
20110227530	Karalis et al.	Sep 2011	A1
20110241618	Karalis et al.	Oct 2011	A1
20120068549	Karalis et al.	Mar 2012	A1
20120228960	Karalis et al.	Sep 2012	A1
20130181541	Karalis et al.	Jul 2013	A1
20150188321	Joannopoulos et al.	Jul 2015	A1

Foreign Referenced Citations (89)

Number	Date	Country
142352	Aug 1912	CA
1309793	Aug 2001	CN
1370341	Sep 2002	CN
1703823	Nov 2005	CN
1993863	Jul 2007	CN
38 24 972	Jan 1989	DE
100 29147	Dec 2001	DE
200 16 655	Mar 2002	DE
102 21 484	Nov 2003	DE
103 04584	Aug 2004	DE
10 2005 036290	Feb 2007	DE
102006044057	Apr 2008	DE
1 296 407	Mar 2003	EP
1335477	Aug 2003	EP
1 521 206	Apr 2005	EP
1 524 010	Apr 2005	EP
2 307 379	May 1997	GB
61-159804	Jul 1986	JP
02-097005	Apr 1990	JP
4-265875	Sep 1992	JP
6-341410	Dec 1994	JP
7-50508	Feb 1995	JP
9-147070	Jun 1997	JP
9-182323	Jul 1997	JP
9-298847	Nov 1997	JP
10-164837	Jun 1998	JP
11-25238	Jan 1999	JP
11-75329	Mar 1999	JP
11-188113	Jul 1999	JP
2001-309580	Nov 2001	JP
2002-10535	Jan 2002	JP
2003-179526	Jun 2003	JP
2004-166459	Jun 2004	JP
2004-201458	Jul 2004	JP
2005-57444	Mar 2005	JP
2005-149238	Jun 2005	JP
2006-074848	Mar 2006	JP
2007-505480	Mar 2007	JP
2007-537637	Dec 2007	JP
2009-501510	Jan 2009	JP
10-2007-0017804	Feb 2007	KR
112842	Jul 2005	SG
WO 9217929	Oct 1992	WO
WO 9323908	Nov 1993	WO
WO 9428560	Dec 1994	WO
WO 9511545	Apr 1995	WO
WO 9602970	Feb 1996	WO
WO 9850993	Nov 1998	WO
WO 0077910	Dec 2000	WO
WO 03036761	Jan 2003	WO
WO 03081324	Oct 2003	WO
WO 03092329	Nov 2003	WO
WO 03096361	Nov 2003	WO
WO 03096512	Nov 2003	WO
WO 2004015885	Feb 2004	WO
WO 2004038888	May 2004	WO
WO 2004055654	Jul 2004	WO
WO 2004073150	Aug 2004	WO
WO 2004073166	Aug 2004	WO
WO 2004073176	Aug 2004	WO
WO 2004073177	Aug 2004	WO
WO 2004112216	Dec 2004	WO
WO 2005024865	Mar 2005	WO
WO 2005060068	Jun 2005	WO
WO 2005019598	Nov 2005	WO
WO 2005109597	Nov 2005	WO
WO 2006011769	Feb 2006	WO
WO 2007008646	Jan 2007	WO
WO 2007020583	Feb 2007	WO
WO 2007042952	Apr 2007	WO
WO 2007084716	Jul 2007	WO
WO 2007084717	Jul 2007	WO
WO 2008109489	Sep 2008	WO
WO 2008118178	Oct 2008	WO
WO 2009009559	Jan 2009	WO
WO 2009018568	Feb 2009	WO
WO 2009023155	Feb 2009	WO
WO 2009023646	Feb 2009	WO
WO 2009033043	Mar 2009	WO
WO 2009070730	Jun 2009	WO
WO 2009126963	Oct 2009	WO
WO 2009140506	Nov 2009	WO
WO 2009149464	Dec 2009	WO
WO 2009155000	Dec 2009	WO
WO 2010030977	Mar 2010	WO
WO 2010039967	Apr 2010	WO
WO 2010090538	Aug 2010	WO
WO 2010090539	Sep 2010	WO
WO 2011062827	May 2011	WO

Non-Patent Literature Citations (190)

Entry
Chinese Office Action, Application No. 201310280724.1; dated Jun. 16, 2015; Applicant: Massachusetts Institute of Technology (7 pages).
Korean Office Action, Application No. 10-2015-7005681; dated May 1, 2015; Applicant: Massachusetts Institute of Technology (6 pages).
Non-Final Office Action for U.S. Appl. No. 14/629,709 dated Aug. 14, 2015 (40 pages).
Non-Final Office Action for U.S. Appl. No. 14/666,683 dated Aug. 17, 2015 (42 pages).
Supplementary European Search Report for European Patent Application No. 09818510 dated Jan. 27, 2017 (6 pages).
Gurhan Alper Kendir et al., “An Efficient Inductive Power Link Design for Retinal Prosthesis”, IEEE, ISCAS 2004, pp. IV-41-IV-44 (2004).
Guoxing Wang et al., “Power Supply Topologies for Biphasic Stimulation in Inductively Powered Implants”, IEEE, pp. 2743-2746 (2005).
Communication from the European Patent Office for Patent Application No. 11 184 066.6 dated Oct. 20, 2014 (7 pages).
Communication from the European Patent Office for Patent Application No. 06 786 588.1 dated Oct. 20, 2014 (6 pages).
Korean Office Action for Patent Application No. 10-2011-7023643 from the Examination Bureau of the Korean Intellectual Property Office dated Sep. 26, 2014 (7 pages).
Jackson, J. D. ,“Classical Electrodynamics”,3rd Edition, Wiley, New York,1999,pp. 201-203.
Tang, S.C et al.,“Evaluation of the Shielding Effects on Printed-Circuit-Board Transformers Using Ferrite Plates and Copper Sheets”,IEEE Transactions on Power Electronics,vol. 17, No. 6,Nov. 2002.,pp. 1080-1088.
Villeneuve, Pierre R. et al.,“Microcavities in photonic crystals: Mode symmetry, tunability, and coupling efficiency”,Physical Review B, vol. 54, No. 11 ,Sep. 15, 1996,pp. 7837-7842.
European Communication for Application No. 06786588.1 dated Aug. 20, 2014 (23 pages).
“Intel CTO Says Gap between Humans, Machines Will Close by 2050”, Intel News Release, (See intel.com/.../20080821comp.htm?iid=S . . . ) (Printed Nov. 6, 2009).
“Intel Moves to Free Gadgets of Their Recharging Cords”, by John Markoff, The New York Times—nytimes.com, Aug. 21, 2008.
“Physics Update, Unwired Energy”, Physics Today, pp. 26, (Jan. 2007) (See http://arxiv.org/abs/physics/0611063.).
“Unwired energy questions asked, answered”, Physics Today, pp. 16-17 (Sep. 2007).
“Wireless Energy Transfer Can Potentially Recharge Laptops, Cell Phones Without Cords”, by Marin Soljacic of Massachusetts Institute of Technology and Davide Castelvecchi of American Institute of Physics (Nov. 14, 2006).
“‘Evanescent coupling’ could power gadgets wirelessly” by Celeste Biever, NewScientistsTech.com, (see http://www.newscientisttech.com/article.ns?id=dn10575&print=true) (Nov. 15, 2006).
“Air Power—Wireless data connections are common—now scientists are working on wireless power”, by Stephen Cass, Sponsored by Spectrum, (See http://spectrum.ieee.org/computing/hardware/air-power) (Nov. 2006).
“Automatic Recharging, From a Distance” by Anne Eisenberg, The New York Times, (see www.nytimes.com/2012/03/11/business/built-in-wireless-chargeing-for-electronic-devices.html?—r=0) (published on Mar. 10, 2012).
“Electro-nirvana? Not so fast”, by Alan Boyle, MSNBC, (Jun. 8, 2007).
“How Wireless Charging Will Make Life Simpler (and Greener)” by David Ferris, Forbes (See forbes.com/sites/davidferris/2012/07/24/how-wireless-charging-will-make-life-simpler-and-greener/print/) (dated Jul. 24, 2012).
“In pictures: A year in technology”, BBC News, (Dec. 28, 2007).
“Lab report: Pull the plug for a positive charge”, by James Morgan, The Herald, Web Issue 2680 (Nov. 16, 2006).
“Look, Ma—no wires!—Electricity broadcast through the air may someday run your home”, by Gregory M. Lamb, Staff writer, The Christian Science Monitor, (See http://www.csmonitor.com/2006/1116/p14s01-stct.html) (Nov. 15, 2006).
“Man tries wirelessly boosting batteries”, by Seth Borenstein, AP Science Writer, Boston.com, (See http://www.boston.com/business/technology/articles/2006/11/15/man—tries—wirelessly—b . . . ) (Nov. 15, 2006).
“Man tries wirelessly boosting batteries”, by Seth Borenstein, The Associated Press, USA Today, (Nov. 16, 2006).
“MIT discovery could unplug your iPod forever”, by Chris Reidy, Globe staff, Boston.com, (See http://www.boston.com/business/ticker/2007/06/mit—discovery—c.html) (Jun. 7, 2007).
“MIT Scientists Pave the Way for Wireless Battery Charging”, by William M. Bulkeley, The Wall Street Journal, (See http://online.wsj.com/article/SB118123955549228045.html?mod=googlenews—wsj) (Jun. 8, 2007).
“MIT's wireless electricity for mobile phones”, by Miebi Senge, Vanguard, (See http://www.vanguardngr.com/articles/2002/features/gsm/gsm211062007.htm) (Jun. 11, 2007).
“Next Little Thing 2010 Electricity without wires”, CNN Money (See money.cnn.com/galleries/2009/smallbusiness/0911/gallery.next—little—thing—2010.smb/) (dated Nov. 30, 2009).
“Outlets Are Out”, by Phil Berardelli, ScienceNOW Daily News, Science Now, (See http://sciencenow.sciencemag.org/cgi/content/full/2006/1114/2) (Nov. 14, 2006).
“Physics Promises Wireless Power” by Jonathan Fildes, Science and Technology Reporter, BBC News, (Nov. 15, 2006).
“Recharging gadgets without cables”, Infotech Online, Printed from infotech.indiatimes.com (Nov. 17, 2006).
“Recharging, The Wireless Way—Even physicists forget to recharge their cell phones sometimes.” by Angela Chang—PC Magazine, ABC News Internet Ventures, (2006).
“Scientists light bulb with ‘wireless electricity’”, www.Chinaview.cn, (See http://news.xinhuanet.com/english/2007-06/08/content—6215681.htm) (Jun. 2007).
“The Big Story for CES 2007: The public debut of eCoupled Intelligent Wireless Power” Press Release, Fulton Innovation LLC, Las Vegas, NV, Dec. 27, 2006.
“The end of the plug? Scientists invent wireless device that beams electricity through your home”, by David Derbyshire, Daily Mail, (See http://www.dailymail.co.uk/pages/live/articles/technology/technology.html?in—article—id=4 . . . ) (Jun. 7, 2007).
“The Power of Induction—Cutting the last cord could resonate with our increasingly gadget-dependent lives”, by Davide Castelvecchi, Science News Online, vol. 172, No. 3, (Week of Jul. 21, 2007).
“The technology with impact 2007”, by Jonathan Fildes, BBC News, (Dec. 27, 2007).
“The vision of an MIT physicist: Getting rid of pesky rechargers” by Gareth Cooks, Globe Staff, Boston.com, (Dec. 11, 2006).
“The world's first sheet-type wireless power transmission system: Will a socket be replaced by e-wall?” Press Release, Tokyo, Japan, Dec. 12, 2006.
“Wireless charging—the future for electric cars?” by Katia Moskvitch, BBC News Technology (See www.bbc.co.uk/news/technology-14183409) (dated Jul. 21, 2011).
“Wireless Energy Lights Bulb from Seven Feet Away—Physicists vow to cut the cord between your laptop battery and the wall socket—with just a simple loop of wire”, by JR Minkel, ScientificAmerican.com, (See http://www.sciam.com/article.cfm?articleid=07511C52-E7F2-99DF-3FA6ED2D7DC9AA2 . . . ) (Jun. 7, 2007).
“Wireless energy promise powers up” by Jonathan Fifties, Science and Technology Report, BBC News, (See http://news.bbc.co.uk/2/hi/technology/6725955.stm) (Jun. 7, 2007).
“Wireless Energy Transfer May Power Devices at a Distance”, ScientificAmerican.com, (Nov. 14, 2006).
“Wireless Energy”, by Clay Risen, The New York Times, (Dec. 9, 2007).
“Wireless power transfer possible”, PressTV, (See http://www.presstv.ir/detail.aspx?id=12754&sectionid=3510208) (Jun. 11, 2007).
“Wireless revolution could spell end of plugs”, by Roger Highfield, Science Editor, Telegraph.co.uk, (See http://www.telegraph.co.uk/news/main.jhtml?xml=/news/2007/06/07/nwireless107.xml) (Jun. 7, 2007).
A. Mediano et al. “Design of class E amplifier with nonlinear and linear shunt capacitances for any duty cycle”, IEEE Trans. Microwave Theor. Tech., vol. 55, No. 3, pp. 484-492, (2007).
Abe et al. “A Noncontact Charger Using a Resonant Converter with Parallel Capacitor of the Secondary Coil”. IEEE, 36(2):444-451, Mar./Apr. 2000.
Altchev et al. “Efficient Resonant Inductive Coupling Energy Transfer Using New Magnetic and Design Criteria”. IEEE, pp. 1293-1298, 2005.
Amnon Yariv et al., “Coupled-resonator optical waveguide: a proposal and analysis”, Optics Letters, vol. 24, No. 11, pp. 711-713 (Jun. 1, 1999).
Andre Kurs et al., “Wireless Power Transfer via Strongly Coupled Magnetic Resonances”, Science vol. 317, pp. 83-86 (Jul. 6, 2007).
Andre Kurs et al., “Simultaneous mid-range power transfer to multiple devices”, Applied Physics Letters, vol. 96, No. 044102 (2010).
Apneseth et al. “Introducing wireless proximity switches” ABB Review Apr. 2002.
Aristeidis Karalis et al., “Efficient Wireless non-radiative mid-range energy transfer”, Annals of Physics, vol. 323, pp. 34-48 (2008).
Baker et al., “Feedback Analysis and Design of RF Power Links for Low-Power Bionic Systems,” IEEE Transactions on Biomedical Circuits and Systems, 1(1):28-38 (Mar. 2007).
Balanis, C.A., “Antenna Theory: Analysis and Design,” 3rd Edition, Sections 4.2, 4.3, 5.2, 5.3 (Wiley, New Jersey, 2005).
Burri et al. “Invention Description” Feb. 5, 2008.
C. Fernandez et al., “A simple dc-dc converter for the power supply of a cochlear implant”, IEEE, pp. 1965-1970 (2003).
Clemens M. Zierhofer et al., “High-Efficiency Coupling-Insensitive Transcutaneous Power and Data Transmission Via an Inductive Link”, IEEE Transactions on Biomedical Engineering, vol. 37, No. 7, pp. 716-722 (Jul. 1990).
European Search Report with regard to Application Serial No. 11184066.6 dated Mar. 20, 2013.
Examination Report for Australia Application No. 2006269374, dated Sep. 18, 2008.
Final Office Action with regard to U.S. Appl. No. 12/639,958 dated Jun. 6, 2013 (18 pages).
Final Office Action with regard to U.S. Appl. No. 12/639,963 dated Jun. 18, 2013 (16 pages).
Final Office Action with regard to U.S. Appl. No. 12/649,635 dated Jun. 20, 2013 (20 pages).
Final Office Action with regard to U.S. Appl. No. 12/649,777 dated Jun. 26, 2013 (17 pages).
Final Office Action with regard to U.S. Appl. No. 12/649,813 dated Jun. 24, 2013 (17 pages).
Final Office Action with regard to U.S. Appl. No. 12/649,852 dated Jun. 27, 2013 (19 pages).
Final Office Action with regard to U.S. Appl. No. 12/649,904 dated Sep. 26, 2013 (23 pages).
Final Office Action with regard to U.S. Appl. No. 12/639,966 dated Oct. 9, 2012 (20 pages).
Final Office Action with regard to U.S. Appl. No. 12/639,967 dated Oct. 5, 2012 (21 pages).
International Preliminary Report on Patentability for International Application No. PCT/US2006/026480, dated Jan. 29, 2008.
International Preliminary Report on Patentability with regard to International Application No. PCT/US2007/070892 dated Sep. 29, 2009.
International Search Report and Written Opinion for International Application No. PCT/US09/43970, dated Jul. 14, 2009.
International Search Report and Written Opinion for International Application No. PCT/US2006/026480, dated Dec. 21, 2007.
International Search Report and Written Opinion for International Application No. PCT/US2007/070892, dated Mar. 3, 2008.
International Search Report and Written Opinion of the International Searching Authority for International Application No. PCT/US2011/027868 dated Jul. 5, 2011.
International Search Report for International Application No. PCT/US09/58499 dated Dec. 10, 2009.
Non-Final Office Action for U.S. Appl. No. 12/639,963 dated Feb. 27, 2014 (19 pages).
Non-Final Office Action for U.S. Appl. No. 12/648,604 dated Dec. 5, 2011.
Non-Final Office Action for U.S. Appl. No. 12/649,635 dated Feb. 27, 2014 (18 pages).
Non-Final Office Action for U.S. Appl. No. 12/649,777 dated Feb. 26, 2014 (16 pages).
Non-Final Office Action for U.S. Appl. No. 12/649,813 dated Feb. 27, 2014 (16 pages).
Non-Final Office Action for U.S. Appl. No. 12/649,852 dated Feb. 27, 2014 (17 pages).
Non-Final Office Action for U.S. Appl. No. 12/726,742 dated May 11, 2012.
Non-Final Office Action for U.S. Appl. No. 13/030,395 dated May 17, 2012.
Non-Final Office Action for U.S. Appl. No. 13/036,177 dated May 15, 2012.
Non-Final Office Action for U.S. Appl. No. 13/040,810 dated May 17, 2012.
Non-Final Office Action for U.S. Appl. No. 13/078,511 dated May 15, 2012.
Non-Final Office Action with regard to U.S. Appl. No. 12/949,580 dated Jun. 17, 2013 (55 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/415,667 dated Oct. 5, 2012 (20 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/639,958 dated Aug. 16, 2012 (21 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/639,963 dated Aug. 31, 2012 (20 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/646,524 dated Oct. 1, 2012 (11 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/649,635 dated Dec. 21, 2012 (41 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/649,777 dated Dec. 24, 2012 (43 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/649,813 dated Dec. 21, 2012 (40 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/649,852 dated Dec. 21, 2012 (41 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/649,904 dated Dec. 28, 2012 (43 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/868,852 dated Oct. 10, 2012 (26 pages).
Non-Final Office Action with regard to U.S. Appl. No. 12/949,544 dated Sep. 5, 2012 (41 pages).
PCT International Search Report and Written Opinion for PCT/US09/59244, dated Dec. 7, 2009, 12 pages.
U.S. Appl. No. 60/908,383, filed Mar. 27, 2007.
Submission of Publication to the Japanese Patent Office for Japanese Application No. 2011-256,729.
Submission of Publication to the Japanese Patent Office for Japanese Application No. 2011-509,705.
Translation of Information Statement by Third Party submitted to the Japanese Patent Office for Japanese Application No. 2011-83009.
D.H.Freedman. “Power on a Chip”. MIT Technology Review, Nov. 2004.
Yates , David C. et al., “Optimal Transmission Frequency for Ultralow-Power Short-Range Radio Links”, IEEE Transactions on Circuits and Systems—1, Regular Papers, vol. 51, No. 7, pp. 1405-1413 (Jul. 2004).
David H. Staelin et al., Electromagnetic Waves, Chapters 2, 3, 4, and 8, pp. 46-176 and 336-405 (Prentice Hall Upper Saddle River, New Jersey 1998).
David Schneider, “A Critical Look at Wireless Power”, IEEE Spectrum, (May 2010).
David Vilkomerson et al., “Implantable Doppler System for Self-Monitoring Vascular Grafts”, IEEE Ultrasonics Symposium, pp. 461-465 (2004).
Electricity Unplugged, Feature: Wireless Energy, Physics World, pp. 23-25 (Feb. 2009).
Esser et al. “A New Approach to Power Supplies for Robots”. IEEE, 27(5):872-875, Sep./Oct. 1991.
European Examination Report dated Jan. 15, 2009 in connection with Application No. 06 786 588.1-1242.
Fenske et al. “Dielectric Materials at Microwave Frequencies”. Applied Microwave & Wireless, pp. 92-100, 2000.
G. Scheible et al., “Novel Wireless Power Supply System for Wireless Communication Devices in Industrial Automation Systems”, IEEE, (2002).
Gary Peterson, “MIT WiTricity Not So Original After All”, Feed Line No. 9, (See http://www.tfcbooks.com/articles/witricity.htm) printed Nov. 12, 2009.
Geyi, Wen. A Method for the Evaluation of Small Antenna Q. IEEE Transactions on Antennas and Propagation, vol. 51, No. 8, Aug. 2003.
H. Sekiya et al. “FM/PWM control scheme in class DE inverter”, IEEE Trans. Circuits Syst. I, vol. 51, No. 7 (Jul. 2004).
Haus, H.A., “Waves and Fields in Optoelectronics,” Chapter 7 “Coupling of Modes—Reasonators and Couplers” (Prentice-Hall, New Jersey, 1984).
Heikkinen et al. “Performance and Efficiency of Planar Rectennas for Short-Range Wireless Power Transfer at 2.45 GHz”. Microwave and Optical Technology Letters, 31(2):86-91, Oct. 20, 2001.
Hirai et al. “Integral Motor with Driver and Wireless Transmission of Power and Information for Autonomous Subspindle Drive”. IEEE, 15(1):13-20, Jan. 2000.
Hirai et al. “Practical Study on Wireless Transmission of Power and Information for Autonomous Decentralized Manufacturing System”. IEEE, 46(2):349-359, Apr. 1999.
Hirai et al. “Study on Intelligent Battery Charging Using Inductive Transmission of Power and Information”. IEEE, 15(2):335-345, Mar. 2000.
Hirai et al. “Wireless Transmission of Power and Information and Information for Cableless Linear Motor Drive”. IEEE 15(1):21-27, Jan. 2000.
J. B. Pendry. “A Chiral Route to Negative Refraction”. Science 306:1353-1355 (2004).
J. C. Schuder et al., “Energy Transport Into the Closed Chest From a Set of Very-Large Mutually Orthogonal Coils”, Communication Electronics, vol. 64, pp. 527-534 (Jan. 1963).
J. Schutz et al., “Load Adaptive Medium Frequency Resonant Power Supply”, IEEE, (2002).
Jackson, J.D., “Classical Electrodynamics,” 3rd Edition, Sections 1.11, 5.5, 5.17, 6.9, 8.1, 8.8, 9.2, 9.3 (Wiley, New York, 1999).
John C. Schuder “Powering an Artificial Heart: Birth of the Inductively Coupled-Radio Frequency System in 1960”, Artificial Organs, vol. 26, No. 11, pp. 909-915 (2002).
John C. Schuder et al., “An Inductively Coupled RF System for the Transmission of 1 kW of Power Through the Skin”, IEEE Transactions on Bio-Medical Engineering, vol. BME-18, No. 4 (Jul. 1971).
Fildes, Jonathan, “Wireless Energy Promise Powers Up”, BBC News, Jun. 7, 2007 (See http://news.bbc.co.uk/2/hi/6725955.stm ).
Joseph C. Stark III, “Wireless Power Transmission Utilizing a Phased Array of Tesla Coils”, Master Thesis, Massachusetts Institute of Technology (2004).
Kawamura et al. “Wireless Transmission of Power and Information Through One High-Frequency Resonant AC Link Inverter for Robot Manipulator Applications”. IEEE, 32(3):503-508, May/Jun. 1996.
Klaus Finkenzeller, “RFID Handbook (2nd Edition)”, The Nikkan Kogyo Shimbun, Ltd., pp. 19, 20, 38, 39, 43, 44, 62, 63, 67, 68, 87, 88, 291, 292 (Published on May 31, 2004).
Finkenzeller, Klaus, RFID Handbook—Fundamentals and Applications in Contactless Smart Cards—, Nikkan Kohgyo-sya, Kanno Taihei, first version, pp. 32-37, 253 (Aug. 21, 2001).
Lee, “Antenna Circuit Design for RFID Applications,” Microchip Technology Inc., AN710, 50 pages (2003).
Lee, “RFID Coil Design,” Microchip Technology Inc., AN678, 21 pages (1998).
Liang et al., “Silicon waveguide two-photon absorption detector at 1.5 μm wavelength for autocorrelation measurements,” Applied Physics Letters, 81(7):1323-1325 (Aug. 12, 2002).
Ahmadian et al., “Miniature Transmitter for Implantable Micro Systems”, Proceedings of the 25th Annual International Conference of the IEEE EMBS Cancun, Mexico, pp. 3028-3031, Sep. 17-21, 2003.
M. V. Jacob et al. “Lithium Tantalate—A High Permittivity Dielectric Material for Microwave Communication Systems”. Proceedings of IEEE TENCON—Poster Papers, pp. 1362-1366, 2003.
Marin Soljacic et al., “Photonic-crystal slow-light enhancement of nonlinear phase sensitivity”, J. Opt. Soc. Am B, vol. 19, No. 9, pp. 2052-2059 (Sep. 2002).
Marin Soljacic, “Wireless nonradiative energy transfer”, Visions of Discovery New Light on Physics, Cosmology, and Consciousness, Cambridge University Press, New York, NY pp. 530-542 (2011).
Microchip Technology Inc., “microID 13.56 MHz Design Guide—MCRF355/360 Reader Reference Design,” 24 pages (2001).
MIT Team Experimentally Demonstrates Wireless Power Transfer, Potentially Useful for Power Laptops, Cell-Phones Without Cords—Goodbye Wires . . . , by Franklin Hadley, Institute for Soldier Nanotechnologies, Massachusetts Institute of Technology (Jun. 7, 2007).
Nikola Tesla, “High Frequency Oscillators for Electro-Therapeutic and Other Purposes”, Proceedings of the IEEE, vol. 87, No. 7, pp. 1282-1292 (Jul. 1999).
Nikola Tesla, “High Frequency Oscillators for Electro-Therapeutic and Other Purposes”, The Electrical Engineer, vol. XXVI, No. 50 (Nov. 17, 1898).
O'Brien et al. “Analysis of Wireless Power Supplies for Industrial Automation Systems”. IEEE, pp. 367-372, 2003.
O'Brien et al. “Design of Large Air-Gap Transfoliners for Wireless Power Supplies”. IEEE, pp. 1557-1562, 2003.
Powercast LLC. “White Paper” Powercast simply wire free, 2003.
Ho et al., “A Comparative Study Between Novel Witricity and Traditional Inductive Magnetic Coupling in Wireless Charging”, IEEE Transactions on Magnetics, vol. 47, No. 5, pp. 1522-1525 (May 2011).
S. Sensiper. Electromagnetic wave propogation on helical conductors. PhD Thesis, Massachusetts Institute of Technology, 1951.
Sakamoto et al. “A Novel Circuit for Non-Contact Charging Through Electro-Magnetic Coupling”. IEEE, pp. 168-174, 1992.
Sekitani et al. “A large-area flexible wireless power transmission sheet using printed plastic MEMS switches and organic field-effect transistors”. [Publication Unknown].
Sekitani et al. “A large-area wireless power-transmission sheet using printed organic transistors and plastic MEMS switches” www.nature.com/naturematerials. Published online Apr. 29, 2007.
Shanhui Fan et al., “Rate-Equation Analysis of Output Efficiency and Modulation Rate of Photomic-Crystal Light-Emitting Diodes”, IEEE Journal of Quantum Electronics, vol. 36, No. 10, pp. 1123-1130 (Oct. 2000).
Soljacic. “Wireless Non-Radiative Energy Transfer—PowerPoint presentation”. Massachusetts Institute of Technology, Oct. 6, 2005.
Someya, Takao. “The world's first sheet-type wireless power transmission system”. University of Tokyo, Dec. 12, 2006.
Splashpower, “Splashpower—World Leaders in Wireless Power,” PowerPoint presentation, 30 pages (Sep. 3, 2007).
T. Aoki et al. Observation of strong coupling between one atom and a monolithic microresonator. Nature 443:671-674 (2006).
Texas Instruments, “HF Antenna Design Notes—Technical Application Report,” Literature No. 11-08-26-003, 47 pages (Sep. 2003).
Thomsen et al., “Ultrahigh speed all-optical demultiplexing based on two-photon absorption in a laser diode,” Electronics Letters, 34(19):1871-1872 (Sep. 17, 1998).
UPM Rafsec, “Tutorial overview of inductively coupled RFID Systems,” 7 pages (May 2003).
Vandevoorde et al. “Wireless energy transfer for stand-alone systems: a comparison between low and high power applicability”. Sensors and Actuators, A 92:305-311, 2001.
Stewart, Will, “The Power to Set you Free”, Science, vol. 317, pp. 55-56 (Jul. 6, 2007).
Yoshihiro Konishi, Microwave Electronic Circuit Technology, Chapter 4, pp. 145-197 (Marcel Dekker, Inc., New York, NY 1998).
Ziaie, Babak et al., “A Low-Power Miniature Transmitter Using a Low-Loss Silicon Platform for Biotelemetry”, Proceedings—19th International Conference IEEE/EMBS, pp. 2221-2224; Oct. 30-Nov. 2, 1997 (4 pages).
Covic et al., “Inductive Power Transfer”, Proceedings of the IEEE, vol. 101, No. 6, pp. 1276-1289 (Jun. 2013).
Covic et al., “Modern Trends in Inductive Power Transfer for Transportation Applications”, IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 1, No. 1, pp. 28-41 (Mar. 2013).
de Beoij et a., “Contactless Energy Transfer to a Moving Load Part I: Topology Synthesis and FEM simulation”, IEEE ISIE, Montreal, Quebec Canada, pp. 739-744 (Jul. 9-12, 2006).
de Beoij et a., “Contactless Energy Transfer to a Moving Load Part II: Simulation of Electrical and Mechanical Transient”, IEEE ISIE, Montreal, Quebec Canada, pp. 745-750 (Jul. 9-12, 2006).
Hui et al., “A Critical Review of Recent Progress in Mid-Range Wireless Power Transfer”, IEEE Transactions on Power Electronics, vol. 29, No. 9, pp. 4500-4511 (Sep. 2014).
Takanashi et al., “A Large Air Gap 3 kW Wireless Power Transfer System for Electric Vehicles”, IEEE, pp. 269-274 (2012).
Wang et al., “Load models and their application in the design of loosely coupled inductive power transfer systems”, IEEE, pp. 1053-1058 (2000).
Australia Patent Examination Report No. 2 for Australian Patent Application No. 2009246310 dated Aug. 21, 2014 (3 Pages).
Chinese Office Action for Chinese Application No. 201310098809.8 dated Mar. 2, 2015 (10 pages)).
Chinese Office Action for Chinese Application No. 201310585104.9 dated Mar. 2, 2015 (10 pages).
Chinese Office Action for Chinese Application No. 201210472059.1 dated Jun. 5, 2014 (54 pages).
Chinese Office Action for Chinese Application No. 201210472059.1 dated Jan. 29, 2015 (13 pages).
Japanese Office Action for Japanese Application No. 2011-509705 dated Jul. 1, 2014 (27 pages).
Korean Office Action for Korean Application No. 10-2009-7022442 dated Oct. 18, 2012 (5 pages).
Korean Office Action for Korean Application No. 10-2011-7023643 dated Oct. 23, 2012 (5 pages).
Non-Final Office Action for U.S. Appl. No. 13/477,459 dated Mar. 12, 2015 (62 pages).
Non-Final Office Action for U.S. Appl. No. 13/789,860 dated Mar. 13, 2015 (51 pages).
Non-Final Office Action for U.S. Appl. No. 14/302,662 dated Mar. 12, 2015 (42 pages).
F. Turki, “A wireless battery charger concept with lightweight and low cost vehicle equipment: eCPS,” in Proc. Conference on Electric Roads & Vehicles, Feb. 2013, pp. 1-21.

Related Publications (1)

	Number	Date	Country
	20140354071 A1	Dec 2014	US

Provisional Applications (1)

	Number	Date	Country
	61101809	Oct 2008	US

Continuations (2)

	Number	Date	Country
Parent	13677407	Nov 2012	US
Child	14458563		US
Parent	12571949	Oct 2009	US
Child	13677407		US

Efficient near-field wireless energy transfer using adiabatic system variations

Information

Patent Number

Date Filed

Date Issued

Inventors

Original Assignees

Examiners

Agents

CPC

Field of Search

CPC

International Classifications

Term Extension

Abstract

Description

Claims

CROSS REFERENCE TO RELATED APPLICATIONS

Government Interests

US Referenced Citations (485)

Foreign Referenced Citations (89)

Non-Patent Literature Citations (190)

Related Publications (1)

Provisional Applications (1)

Continuations (2)