Hidden Time-Reversal Symmetry in Dissipative Reciprocal Systems
HHidden Time-Reversal Symmetry in Dissipative Reciprocal Systems M ÁRIO G. S ILVEIRINHA , University of Lisbon – Instituto Superior Técnico and Instituto de Telecomunicações, Avenida Rovisco Pais 1, 1049-001 Lisboa, Portugal *[email protected]
Abstract:
It is proven, without using the microscopic reversibility argument of Onsager, that lossy reciprocal systems have a hidden time-reversal symmetry. The key idea is that the dissipation channels of lossy dielectrics can be mimicked by a distributed network of lossless transmission lines. It is highlighted that the reciprocity of lossy dielectrics is fundamentally rooted on the hidden time-reversal invariance and on linearity of the materials. Furthermore, it is demonstrated that the upper-half plane response of dissipative materials can be approximated as much as desired by the response of some lossless material.
1. Introduction
The laws that rule the propagation of light in free space are invariant under the “time reversal” operation [1, 2, 3]. This means that when the arrow of time is flipped – so that the dynamics of the electromagnetic wave is reversed in time similar to a movie played backwards –the time-reversed wave remains compatible with the laws of electromagnetism, i.e., it still satisfies the Maxwell equations. An extraordinary consequence of the time-reversal invariance is that the propagation of light in typical material platforms is inherently bi-directional. For example, if a wave can go through some channel with no back-reflections, then the time-reversed wave can go through the same channel but propagating in the opposite direction. This rather profound and intriguing bi-directional character of photonic systems is usually regarded a consequence of the Lorentz reciprocity law [4-5], which is more general than the time-reversal invariance. Indeed, the Lorentz reciprocity theorem can be applied to macroscopic systems with dissipative elements, which are not invariant under a time-reversal transformation. The reciprocity of dissipative systems is usually justified by the microscopic reversibility of physical processes [2, 6], i.e., relying on the time-reversal invariance of (classical) physics at the microscopic level (the quantum world is also mostly ruled by time-reversal invariant laws, with a few exceptions related to nuclear physics which are unimportant for the effects discussed here). Microscopic reversibility and statistical mechanics were used by Onsager to derive reciprocal relations for irreversible processes [2, 6]. In this article, I show that it is unnecessary to invoke microscopic arguments to establish a direct link between the reciprocity of macroscopic systems and time-reversal invariance. To this end, I prove that dissipative dielectrics have a hidden time-reversal symmetry, as the relevant dissipation channels can be mimicked by a distributed network of infinitely-extended time-reversal invariant lossless transmission lines. Using this result, I prove that the Lorentz reciprocity is ultimately a consequence of the hidden time-reversal invariance, linearity and conservation of energy. Furthermore, I also show that the response of a passive reciprocal dissipative material can be approximated as much as desired (in a given time interval with finite duration) by that of a time-reversal invariant lossless medium. Previous studies of quantum electrodynamics in lossy material platforms [7-12] have revealed a rather profound link between dissipation, “open systems”, and interactions with a “bath” of oscillators. For example, Ref. [11] demonstrates that a lossy electric circuit escribed by some impedance function can be implemented using lossless networks of idealized inductors and capacitors. Conversely, the dynamics of some classical inhomogeneous systems is determined by a “bath” of harmonic oscillators, and due to this reason their response is effectively dissipative, even when they are formed exclusively by lossless materials [13, 14]. This article exploits and extends these known paradigms to unveil the hidden-symmetry of reciprocal materials.
2. Time-reversal symmetry
The propagation of light in composite media is ruled by the macroscopic Maxwell equations: t BE , t DH j . (1) The time-reversal operation transforms the electromagnetic fields , E H as [2]: TR , , t t E r E E r , TR , , t t H r H H r (2) The electric displacement vector D is transformed similarly to E , and the induction field B and the electric current density j in the same manner as H . Thus, the former is said to be even under the time reversal operation, whereas the latter are odd . Similar to the original fields, the time-reversed fields TR TR TR , , ,...
E H j satisfy the Maxwell’s equations (1). A different (pseudo-) time-reversal symmetry for photonic systems was introduced in Ref. [15]. For a idealized lossless dielectric material with an instantaneous response the constitutive relations are D r E and B H . Evidently, the time-reversed fields are linked in the same manner, i.e.,
TR TR D r E and
TR TR0 B H . Thereby, dielectric platforms formed by lossless dielectrics are time-reversal invariant. This property also holds true when the (lossless) materials are dispersive. Indeed, the electrodynamics of lossless materials can be formulated as a Schrödinger-type time evolution problem [16-19], which can be shown to be time-reversal invariant for dispersive isotropic dielectrics. Moreover, typical lossless nonlinear platforms are time-reversal invariant [20]. For linear media, the Maxwell equations admit time-harmonic solutions with a time variation i t e , with the real-valued oscillation frequency. The time-dependent electromagnetic fields and current density are written in terms of complex amplitudes in the usual way, e.g., , Re i t t e E r E r . Evidently, acts on the electric field as * , Re Re i t i t t e e E r E r E r . Thus, in the frequency domain the time-reversal operation is closely related to a complex conjugation [21]. The complex field amplitudes are transformed as:
TR *
E r E E r , TR *
H r H H r , (3a) TR * j r j j r , (3b) The complex amplitudes of D and B are transformed similarly, taking into account the parity (even or odd) of each vector field. In the lossless case, the response of a material is time-reversal invariant if and only if the material is reciprocal; for an explicit proof see Appendix A. In contrast, dissipative dielectrics are not time-reversal invariant, even though they are reciprocal.
3. Hidden symmetry
Realistic materials are lossy from an electromagnetic point of view: as a wave propagates in a medium part of its energy is irreversibly lost in the form of heat. The absorption effect is ncluded in the dielectric response. Due to this reason, the equations of macroscopic electrodynamics are not time-reversal invariant when the materials are lossy. Nevertheless, as is well known, the response of dissipative systems is constrained by reciprocity relations [2, 6]. As mentioned in the Introduction, the equations that rule the propagation of light in natural media are time-reversal invariant at the microscopic level [2, 6]. The time-reversal of some wave process in a dissipative system requires the medium to pump the time-reversed field distribution, to give back all the energy originally dissipated as heat. No matter how strange as this may look, the theory predicts that if the system could be prepared in a suitable microscopic initial state it could be possible to generate the time-reversed wave dynamics, starting with an apparently chaotic thermal bath [1]. In practice, the preparation of the required initial state is unrealistic. The microscopic reversibility is a central argument of Onsager’s theory [2, 6]. Next, I show –without invoking any microscopic arguments– that a reciprocal dissipative macroscopic system has a hidden time-reversal symmetry.
Fig. 1. a) Dissipative RL circuit. b) Equivalent time-reversal invariant lossless circuit where the resistor is implemented with a semi-infinite transmission line.
Consider the circuit of Fig. 1a, which is formed by a resistor and an inductance fed by a voltage generator. The current circulating in the circuit is the solution of the differential equation: di tR i t L v tdt . (4) The voltage and the current are transformed under a time-reversal transformation as: v t v t , i t i t . (5) Thus, the system is time reversal invariant only when R . The transfer function of the system in the frequency domain is / 1 I V R i L , and as expected has a pole in the lower-half frequency plane due to the dissipation in the resistor. From the theory of transmission lines, a semi-infinite lossless transmission line with characteristic impedance Z is equivalent to a resistor R Z . Therefore, the dissipative circuit of Fig. 1a has exactly the same response as the lossless circuit of Fig. 1b. This property was previously discussed in the context of the analysis of dissipative quantum systems [11]. The circuit of Fig. 1b is manifestly time-reversal invariant. Thus, it follows that an open time-reversal invariant system can mimic perfectly the response of a dissipative system. One may say that the circuit of Fig. 1a has a hidden time-reversal symmetry . ikewise, the response of a lossy dielectric may be perfectly reproduced by an open time-reversal invariant lossless system. To show this, I consider without loss of generality that the dielectric response is determined by the Lorentz model:
202 20 i . (6) The polarization vector P of a material with the Lorentz dispersion is described by the differential equation: t t P P P E . (7) The polarization vector P is even under a time-reversal: , , t t P r P r . Thus, the dynamics determined by (7) is not time-reversal invariant when the material is lossy ( ). Similar to the circuit of Fig. 1a, it is possible to find a lossless open system with the same response as that of the lossy material (6). Here is an example inspired by transmission line theory: u t AZ P V P E , (8a) u u
Lu t
V I , u u C A uu t t
I V P . (8b) I introduced two vector fields , u u V I that depend on the spacetime coordinates , t r and on a fictitious 4 th space-coordinate u . The fields , P E depend only on , t r . The vector field , , x y zu u u u V V V V has the units of a voltage and the field , , x y zu u u u I I I I has the units of a current. The parameters , L C play the roles of a distributed inductance and capacitance, respectively, and / Z L C is the characteristic impedance. The parameter A has units of area so that / A t P has units of current. The polarization vector is coupled to the Maxwell equations in the usual way: / t E H , / / t t H j P E . (8c) The system (8) models a three-dimensional system (the material) coupled to a 4 th dimension through a distributed network of 1D transmission lines. The transmission lines are infinitely extended along the 4 th dimension, u . The coupling between the material and the lines is determined by the voltages calculated at the point u . The transmission line associated with , i iu u V I ( i=x,y,z ) is fed by a current generator placed at 0 u with amplitude / i A P t . From (8b) it is clear that / i i iu u I I A P t and
00 0 i i iuu u
V V V . If one imposes that , i iu u V I satisfy radiation boundary conditions so that the energy can only flow away from the generator, it follows that / i iu u V I Z for 0 u and / i iu u V I Z for u . Thus, in these conditions, iu i V Z A P t . Substituting this result into (8a) we recover the dissipative Lorentz model (7). Therefore, the dynamics of (7) is exactly reproduced by the lossless open system described by (8a)-(8b) when radiation boundary conditions are imposed on the transmission lines. Moreover, the system (8) is time-reversal invariant because the dynamics of the time-reversed fields [ , i iu u V I are transformed as in (5)] is described by the same equations [Eq. (8)]. Thus, the dissipative Lorentzian response (6) has a hidden time reversal symmetry, consistent with the microscopic reversibility discussed by Onsager [2, 6]. he enunciated results can be readily extended to other more complex dispersive models (and to the anisotropic case), because the permittivity of a generic dispersive material can be written as a sum of Lorentz poles (see Appendix B).
Consider again the system of Fig. 1a. It is interesting to note that the response of the resistor may also be imitated by a line with finite length l terminated in short-circuit, i.e., by a closed lossless system. Specifically, for an excitation that starts at t =0, the response of the resistor is perfectly mimicked by the line in the time interval p t T with / p p T l v the time of propagation from the beginning to the end of the line. Thus, as l the response of the line becomes coincident with that of the dissipative system for longer and longer time intervals. Note that the impedance of the short-circuited line is purely reactive. When l the line is semi-infinite, the system is open, and its response mimics exactly that of the resistor. Furthermore, for an excitation of the system with a time variation i t t i t e e e with i in the upper-half frequency plane (UHP) ( ), the line terminated with a short circuit may reproduce well the resistor response for any time instant. The reason is that when l is sufficiently large the echoes coming from the transmission line have negligible amplitude as compared to the signal coming from the generator, because the amplitude of the latter grows exponentially with time ( t e with ). Thus the UHP response of a dissipative circuit element (the resistor) can be approximated as much as desired by that of a lossless reactive circuit (line terminated with a short circuit). It is evident from the transmission-line model developed in Sect. 3.1 that the response of a lossy dielectric can be imitated during an arbitrarily large time period by the response of some lossless material. Furthermore, the UHP response of a lossy dielectric can be approximated as much as desired by that of a lossless time-reversal invariant material. An explicit proof of this result (not based on the transmission-line model of Sect. 3.1) is given in Appendix B. In quantum theory it is often useful to model dissipation effects through the interaction of some Hermitian system with a bath of harmonic oscillators [7-12]. The formalism of Appendix B connects closely with such an idea, as it unveils that the dielectric function of a dissipative material is determined by a “bath” (continuum) of lossless Lorentz-type oscillators.
4. Reciprocity theorem
As discussed in Sect. 2, the concept of reciprocity is broader than the concept of time-reversal invariance. This property might suggest that reciprocity is rooted on some special characteristic (e.g., on the mathematical structure) of the Maxwell equations, beyond the time-reversal invariance. Next, I prove that the Lorentz reciprocity is essentially a consequence of the linearity and of the hidden time-reversal invariance of the materials.
To begin with, I focus on lossless time-reversal invariant systems, and I derive the Lorentz reciprocity theorem without invoking directly the Maxwell’s equations or the constitutive relations of the materials. The proof uses simply the conservation of energy and the “linearity” of the equations. To this end, consider a generic solution , E H of the Maxwell’s equations in time-harmonic regime (with real-valued) created by some external current density j . Assuming that the system is lossless, the conservation of energy (Poynting theorem) implies that av ext,av p S , where *av
1/ 2 Re
S E H is the time-averaged energy-density flux Poynting vector) and *ext,av
1/ 2 Re p E j is the time-averaged power extracted from the external current j per unit of volume. The energy balance constraint can be expressed as * *ext Re , , 0
E H E j , (9) where for convenience I introduced the bilinear forms , E H E H and ext , E j E j . The bilinear forms satisfy * * * , , E H E H and * * *ext ext , ,
E j E j . Let , E H and , E H be field distributions created respectively by j and j . Then, by linearity , E E H H is the solution of the Maxwell equations for the excitation j j , and thus it must satisfy (9). This is only possible if: * * * *ext ext
Re , , , , 0
E H E H E j E j . (10) Now, if the system has the time-reversal symmetry * * * , ,
E H j is also a solution of the Maxwell equations in the same physical platform. Thereby, it is possible to replace , ,
E H j by * * * , , i e E H j in Eq. (10), and the resulting equation remains true for arbitrary , E H and , E H . The constant factor i e is introduced by convenience (this is possible due to the linearity of the system). Thus, it follows that for any : ext ext Re , , , , 0 i e E H E H E j E j . (11) If w is a complex number and Re 0 i e w for an arbitrary then w . Thereby, it follows that ext ext , , , , E H E H E j E j , (12) which is exactly the Lorentz reciprocity theorem in the differential form. The above analysis is rather general, and only requires that the system is linear, non-dissipative, time-reversal invariant, and that the energy density flux is determined by E H (local media). Clearly, the derivation can be readily generalized to other physical systems (unrelated to electromagnetism) wherein the energy density flux and the extracted power density are bilinear forms of the relevant fields. This property helps to understand why the reciprocity constraint is so ubiquitous in “wave systems” [2, 6].
The analysis of the previous subsection can be extended to general lossy reciprocal materials with a hidden-time reversal symmetry. For simplicity, the following discussion is focused on the model of Sect. 2.1. For the system described by Eqs. (8), the conservation of energy in time-harmonic regime is expressed by av ext,av d,av p p S with *d,av 0, , 0 u u p V I , Z A , and , 0 0 , 0 , u u u I I I . Proceeding as in Sect. 4.1 and using the time-reversal invariance of the system (8), it is possible to show that: ext extd , 0 , 0 d , 0 , 0 , , , , , , u u u u
E H E H E j E jV I V I (13) ith d , V I V I . As seen in Sect. 4.1, when radiation boundary conditions are enforced on the transmission lines the corresponding voltages and currents are linked as
0, 0 0, / 2 u u Z V I . In these conditions, the term d , 0 , 0 d , 0 , 0 , , u u u u
V I V I vanishes, and Eq. (13) yields the standard reciprocity theorem in a dissipative material system. Therefore, the reciprocity of dissipative systems is deeply rooted on the hidden time reversal invariance.
5. Summary
The response of an idealized “open” lossless time-reversal invariant system can be made exactly coincident with that of a lossy dielectric. Thus, dissipative (reciprocal) dielectrics have a hidden time-reversal symmetry and their response is fundamentally constrained by time-reversal invariance. In particular, it was shown that the reciprocity of dissipative systems is ultimately a consequence of the hidden time-reversal symmetry and of the linearity of the materials. Furthermore, it was demonstrated that the dynamics of a lossy material can be precisely mimicked by a lossless material during an arbitrarily large time interval. In addition, the UHP response of a lossy material can be approximated as much as desired by that of some lossless material, without any restriction on the duration of the excitation.
Appendix A
Suppose that the complex amplitudes of the electromagnetic fields are related by bianisotropic constitutive relations of the type, cc M D EB H . (A1) The material matrix M is written in terms of permittivity, permeability, and magneto-electric coupling tensors, in a standard way. Then, the time-reversed fields TR TR , D B and
TR TR , E H are linked by a material matrix TR M given by: TR * z z
M σ M σ , with z . (A2) Note that because of the reality condition, * * M M , the above formula is equivalent to TR z z M σ M σ . A system formed by bianisotropic materials is time-reversal invariant when TR M M . From Eq. (A2), the general conditions for time-reversal invariance are * , * , and * and * . For example, dissipative isotropic dielectrics are not time reversal invariant because * . For a lossless medium, the material matrix is necessarily Hermitian symmetric, † M M , for real-valued [22]. In these conditions, Eq. (A2) reduces to (the superscript T stands for the matrix transpose): TR Tz z
M σ M σ , for real-valued. (A3) his constraint is equivalent to the well-known Lorentz reciprocity relations: T , T , and T . Thus, lossless time-reversal invariant linear materials are always reciprocal. Appendix B
In this Appendix, I show that a generic passive reciprocal medium can be regarded as the limit of some sequence of time-reversal invariant lossless materials. For simplicity, I restrict the analysis to isotropic dissipative dielectrics described by some dispersive permittivity . The permittivity is assumed to be a meromorphic function of frequency and hence has a partial-fraction expansion of the form: , * *Re 0, , , p n n n nn p n p n p n . (B1) Here, , p n are the poles and n are the residues of the permittivity. The second identity uses the reality-condition ** . The passivity of the material requires that the poles are in the lower-half frequency plane (for a time variation i t e ), i.e., , Im 0 n p n . Furthermore, in order that
Im 0 in the positive real axis, it is necessary that decays as in the upper-half of the imaginary frequency axis [22]. This is only possible if n is a positive real-number. Thus, writing , p n n n i , it follows that the permittivity of a passive material is necessarily of the form: nn n n n n i i . (B2) In the above, , n n determine the resonant frequency and damping rate, respectively, and are positive numbers. As expected, each term of the sum is equivalent to a Lorentz oscillator [16]. For example, when the summation restricted to the n =1 term, one obtains the Lorentzian response (6) with , , and . Now the key observation is that in the upper-half frequency plane (UHP) ( Im 0 ) one has: nn n n n di (cid:0) . (B3) The identity can be verified using the Cauchy theorem. Therefore, the dielectric function in the UHP may be written as: d , n nn n n . (B4) Suppose that the integral is approximated by a discrete summation of the form , with M mm M . The approximation is obtained by sampling M points of the integration domain spaced by . Evidently, one has lim M when is in the UHP. The function is analytic in the omplex plane, with the exception of the real-frequency axis where has poles at m ,
1, 2,... m . Since for real-valued / 0 d d , the function describes the response of a lossless passive system [22]. In addition, the material response is time-reversal invariant due to . In any compact region of the UHP, can be approximated as much as desired by some with the sampling period small enough and M large enough. Thus, it follows that a lossy reciprocal dielectric may always be regarded as the limit of a sequence of time-reversal invariant lossless materials. This further supports that the reciprocal properties of a lossy dielectric are inherited from the hidden time-reversal symmetry. Even though the previous discussion was focused on isotropic dielectrics, the results can be generalized to other reciprocal and passive material platforms. Acknowledgments
This research was funded by the IET under the A F Harvey Engineering Research Prize and by Fundação para Ciência e a Tecnologia (FCT) under projects PTDC/EEITEL/4543/2014 and UID/EEA/50008/2017. The author gratefully acknowledges stimulating discussions with Dimitrios Tzarouchis.
References
1. R. Feynman, R. Leighton, M. Sands,
The Feynman lectures on physics , California Institute of Technology, 1963. 2. H. B. G. Casimir, “Reciprocity theorems and irreversible processes”,
Proc. IEEE , , 1570, (1963). 3. M. G. Silveirinha, “Time-reversal Symmetry in Antenna Theory”, arXiv:1902.0958 (2019). 4. R. J. Potton, “Reciprocity in optics”, Rep. Prog. Phys. , , 717-754, (2004). 5. C. Caloz, A. Alù, S. Tretyakov, D. Sounas, K. Achouri, Z.-L. Deck-Léger, “Electromagnetic Nonreciprocity”, Phys. Rev. Applied , , 047001, (2018). 6. L. Onsager, “Reciprocal relations in irreversible processes I”, Phys. Rev. , , 405, (1931). 7. J. J. Hopfield, “Theory of the Contribution of Excitons to the Complex Dielectric Constant of Crystals”, Phys. Rev. , , 1555, (1958). 8. A. O. Caldeira, A. J. Leggett, “Influence of Dissipation on Quantum Tunneling in Macroscopic Systems”, Phys. Rev. Lett. , , 211, (1981). 9. A. O. Caldeira, A. J. Leggett, “Quantum tunnelling in a dissipative system”, Ann. Phys. , , 374-456, (1983). 10. B. Huttner, S. M. Barnett, “Quantization of the electromagnetic field in dielectrics”, Phys. Rev. A , , 4306, (1992). 11. U. Vool, M. Devoret, “Introduction to quantum electromagnetic circuits”, Int. J. Circ. Theor. Appl. , , 897–934, (2017). 12. W. E. I. Sha, A. Y. Liu , W. C. Chew, “Dissipative Quantum Electromagnetics”, IEEE J. Multiscale Multiphys. Comput. Tech. , , 198, (2018). 13. D. C. Tzarouchis, N. Engheta, A. Sihvola, “Resonant Graded-Index Plasmonic Nanoscatterers: Enabling Unusual Light-Matter Interactions”, (2018). 14. J. B. Pendry, P. A. Huidobro, Y. Luo, E. Galiffi, “Compacted dimensions and singular plasmonic surfaces”, Science , , 915-917, (2017). 15. M. G. Silveirinha, “PTD symmetry protected scattering anomaly in optics”, Phys. Rev. B , , 035153, (2017). 16. A. Raman, S. Fan, “Photonic Band Structure of Dispersive Metamaterials Formulated as a Hermitian Eigenvalue Problem”, Phys. Rev. Lett. , , 087401, (2010). 17. M. G. Silveirinha, “Chern Invariants for Continuous Media”, Phys. Rev. B , , 125153, (2015). 18. M. G. Silveirinha, “Topological classification of Chern-type insulators by means of the photonic Green function”, Phys. Rev. B , , 115146, (2018). 19. M. G. Silveirinha, “Modal expansions in dispersive material systems with application to quantum optics and topological photonics”, chapter to appear in “Advances in Mathematical Methods for Electromagnetics”, (edited by Paul Smith, Kazuya Kobayashi) IET, (available in arXiv:1712.04272). 20. D. E. Fernandes, M. G. Silveirinha, “Asymmetric Transmission and Isolation in Nonlinear Devices: Why They Are Different”, IEEE Antennas Wirel. Propag. Lett. , , 1953, (2018). 21. J. B. Pendry, “Time Reversal and Negative Refraction”, Science ,322