# Spectral and transport properties of time-periodic PT-symmetric tight-binding lattices

## Abstract

We investigate the spectral properties and dynamical features of a time-periodic -symmetric Hamiltonian on a one-dimensional tight-binding lattice. It is shown that a high-frequency modulation can drive the system under a transition between the broken- and the unbroken- phases. The time-periodic modulation in the unbroken- regime results in a significant broadening of the quasi-energy spectrum, leading to a hyper-ballistic transport regime. Also, near the -symmetry breaking the dispersion curve of the lattice band becomes linear, with a strong reduction of quantum wave packet spreading.

###### pacs:

03.65.-w, 11.30.Er, 72.10.Bg, 42.82.Et,## I Introduction

Non-Hermitian Hamiltonians play an important role in physics (1). They are introduced, for example, to model open systems and dissipative phenomena in a wide variety of physical contexts (see, for instance, (1); (2)). A special class of non-Hermitian Hamiltonians is provided by complex potentials having parity-time () symmetry (3); (4); (5), that is invariance under simultaneous parity transform (: , , where and stand for momentum and position operators, respectively) and time reversal (: , , ). An important property of Hamiltonians is to admit of an entirely real-valued energy spectrum below a phase transition symmetry-breaking point. Unbroken -symmetric Hamiltonians possess a complete set of eigenvectors and have originally attracted great attention for the possibility to formulate a consistent quantum mechanical theory in a non-Hermitian framework (5); (4); (6); (7). The appeal of such an idea originated in the years-long emphasized relevance of systems with -symmetry in quantum field theory (see e.g. (8); (9)), but have found interest and applications in several other physical fields, such as in classical mechanics (10), magnetohydrodynamics (11), cavity quantum electrodynamics (12), and electronics (13). Very recently, the peculiar features of -symmetry have found important applications in optics, with the theoretical proposal and experimental demonstration of innovative materials and optical devices based on -invariant complex dielectric functions (14); (15); (17); (16); (19); (18); (20); (21); (22); (23); (24); (25). optical media offer rather unique transport properties, including, as example, double refraction, power oscillations, and nonreciprocal diffraction patterns (17), unidirectional Bragg scattering and invisibility (18); (27); (24); (28); (29); (30); (25), giant Goos-Hänchen shift (32), and simultaneous perfect absorption and laser behaviour (33); (34).

Though the physics of -symmetric Hamiltonians has been widely explored, most of the attention has been devoted to static (i.e. time-independent) potentials. It is only very recently that the study of time-dependent -symmetric potentials has attracted an increasing attention (35); (36); (37); (38); (39); (40); (41); (42). The time-dependent part of the Hamiltonian considered in most of such previous studies is an Hermitian driving term acting on a static non-Hermitian potential, whereas a genuine (and eventually strong) non-Hermitian driving has been investigated so far solely for a simple low-dimensional (dimeric) -symmetric system (40); (41) and for the -symmetric kicked rotor model (37) in connection to dynamical localization and quantum chaos. On the other hand, time-dependent Hermitian Hamiltonians are of fundamental importance in many branches of physics, with particular interest to time-periodic Hamiltonians entailing a plethora of phenomena ranging from Rabi oscillations (see e.g. (43)) and Autler-Townes splitting (44), to dynamic localization (45) and coherent destruction of tunneling (46); (47), just to mention a few (see also (48) for a review on driven quantum tunneling). Therefore, it is of great interest to extend the theoretical investigations of coherently-driven quantum systems within the context of -symmetric Hamiltonians, especially for spatially-extended systems, like in complex crystals, where symmetry can strongly affect the transport properties.

In the present paper we consider a time-dependent -symmetric potential defined on a one-dimensional tight-binding lattice with nearest neighbor hopping rate, in which a time-periodic complex modulation of the potential acts as a non-Hermitian periodic driving of a static Hermitian lattice. We show that in the time-periodic system the frequency and amplitude of the driving can control the symmetry phase transition, and provide a detailed investigation of the spectral properties and dynamical features of the driven lattice below the symmetry breaking point. A comparison with the un-driven Hamiltonian and with Hermitian periodic driving (that is a real potential modulation) is given to highlight the peculiar effects of non-Hermitian driving on the quasi-energy spectrum of the system and related transport properties. In particular, we show that the time-periodic modulation in the unbroken regime results in a significant broadening of the quasi-energy spectrum, leading to a hyper-ballistic transport regime and to a strong reduction of quantum wave packet spreading near the -symmetry breaking threshold.

The paper is organized as follows. In Sec. II we introduce the time-periodic -symmetric lattice model and investigate the effects of the time-periodic modulation on the -symmetry breaking point of the system. In Sec. III we discuss the peculiar spectral features and transport properties of the lattice under strong high frequency modulation, showing the appearance of enhanced ballistic transport and reduction of wave packet dispersion as the -symmetry breaking boundary is approached. Finally, in Sec. IV the main conclusions and future developments are outlined.

## Ii A time-periodic -symmetric tight-binding lattice

We consider here a -invariant tight-binding Hamiltonian describing the hopping dynamics of a single particle on a one-dimensional binary superlattice, with a superposed modulation of the ’imaginary’ part of the site energies (loss or gain term) with amplitudes at alternating lattice sites. The Hamiltonian of the superlattice reads:

(1) |

where

(2) |

(3) |

In above Eqs. (2)-(3), is the hopping rate between adjacent lattice sites and is the Wannier state localized at the -th site. Note that the resulting total Hamiltonian is -invariant, regardless the temporal dependence detailed by . For a static (i.e. time-independent) real-valued the -invariant lattice model Eqs.(1-3) was previously investigated in Refs.(18); (49) and was shown to be always in the broken phase for any infinitesimally-small value of . For a time-periodic but imaginary , corresponding to Hermitian ac-driving, the lattice model Eqs.(1-3) was previously investigated in Refs.(50); (51) and it was shown to realize, under certain modulation conditions, approximate suppression of particle spreading (dynamic localization).

The state vector of the system in Wannier basis representation can be written as:

(4) |

where is the complex amplitude for occupation of the Wannier state. Given above decomposition of the state vector, standard projection technique provides the following evolution equations for the amplitudes :

(5) |

where the dot stands for the derivative with respect to time . Let us assume that the driving amplitude is a time-periodic function of period with zero mean value, i.e. . The spectral and transport properties of the ac-driven lattice can be throughly investigated by computing the quasi-energy spectrum of , where the quasi-momentum varies in the interval (see e.g. (48) and references therein). In the following analysis, we will consider a square-wave ac modulation for the sake of definiteness, however similar results are obtained for different types of modulation, such as for a sinusoidal modulation. For a square-wave driving , where is the modulation frequency and an initial phase. In this case the monodromy matrix (connecting the solution to the time-periodic system of Eq. (5) in Fourier space over one oscillation cycle) can be analytically determined as , with:

(6a) | |||

(6b) |

where , , and . Details are given in the Appendix.

Once the complex eigenvalues of are computed as a function of the quasi-momentum , the real part and imaginary part of the quasi-energy spectrum can be evaluated as follows (see the Appendix for details):

(7a) | |||

(7b) |

The nature of the -symmetry in the plane of parameters can be throughly determined by inspecting (Fig. 1). We find that the in the low-frequency regime, the system exhibits a broken -symmetry phase, with complex quasi-energies [c.f. left inset in Fig. 1]. Above a minimum value of the frequency , the system can be driven to the unbroken -symmetric phase (shaded area), having a purely real quasi-energy spectrum [c.f. right inset in Fig. 1], provided that the driving amplitude is lower than a threshold value which is an almost linear function of the modulation frequency [solid line in Fig. 1]. Therefore, contrary to the static case (18); (49), a broad range of modulation parameters exists where the driven Hamiltonian defined by Eqs. (1-3) has a real-valued energy spectrum and can exhibit potentially interesting novel transport properties.

It is worth mentioning briefly that the time-periodic -symmetric Hamiltonian above discussed can find a physical implementation as example in optical systems. Optical realizations of -symmetric photonic lattices, based on light propagation in coupled waveguides with alternating regions of optical gain and absorption, have been discussed in many previous works (see, e.g., (15); (17); (16); (18)). In Fig. 2 we show a sketch of a possible implementation of the one-dimensional tight-binding lattice with symmetric time-periodic driving, in the form of a linear array of tunneling-coupled optical waveguides with an alternation of loss and gain sections along propagation axis . This represents somehow a generalization of the original approach to symmetric photonic lattices, recently demonstrated experimentally for a static -symmetric directional coupler (21). Another possibility to implement our Hamiltonian is to exploit the recent proposal of mimicking -symmetric optical crystals with a network of coupled fiber loops (24). This would represent a time-discrete version of the photonic lattice of Fig. 2, where time evolution is implemented as a repeated circulation in the fiber loops (instead of space propagation along the optical axis of a photonic structure), and the loops can provide on demand gain or loss in equal amount by means of optical amplifiers and amplitude modulators.

## Iii Spectral and transport properties

The quasi-energy spectrum in a time-periodic lattice is known to govern the transport properties of one particle in the lattice. It is worth noting that ac-driving is a powerful and ubiquitous technique of quantum control within Hermitian systems, allowing one to shrink the quasi-energy spectrum and eventually collapse the quasi-energies, giving rise to dramatic effects on transport, like coherent suppression of tunneling and dynamic localization (see e.g. (48) and references therein). In this section we investigate the transport properties of non-Hermtian ac-driving of the -symmetric lattice, highlighting the appearance of a novel transport regime, referred to as hyper-ballistic motion, and the possibility to strongly reduce wave packet spreading.

In Fig. 3 we compare the quasi-energy spectra of the Hamiltonian considered in the previous section for a fixed value of the modulation frequency and different values of the driving amplitude . The energy spectrum of the undriven Hamiltonian is also shown (dotted line). We found that as is increased, the quasi-energy spectrum is progressively broadened, contrary to the case of Hermitian driving (also reported in dashed line for comparison) that always gives rise to a shrink of the spectrum (50). The non-Hermitian nature of the periodic driving (even in the unbroken -symmetric phase) is thus qualitatively distinct from an Hermitian driving acting on the same static Hamiltonian. Moreover, as the -symmetry breaking is approached, the quasi-energy spectrum attains a progressive linear dispersion relation, culminating in a zero group-velocity dispersion in the whole Brillouin zone, with the exception of the zone center and edges.

The impact of the spectral features above detailed on the transport properties of the system is twofold: (i) the broadening of the quasi-energies is expected to enhance the ballistic motion of the particle, i.e. to introduce a novel regime (with no counterpart in Hermtian lattices) that we can refer to as hyper-ballistic motion: the group velocity in the lattice of a particle wave packet is not limited by the hopping rate between adjacent sites; (ii) the linearization of the dispersion curve near the symmetry breaking point results in a strong reduction of the spreading (diffraction) for a particle wavepacket.

We have checked the above predictions by numerical simulations of Eqs. (5) for either narrow (single-site) and broad wave packet excitation. As is well known in the Hermitian ac-driving case, the phase (i.e. time origin) of the modulation has a noticeable impact onto the wave packet dynamics (see, for instance, (52)). For example, in the high-frequency limit, the intial phase is responsible to an additional phase term in the renormalized tunneling rate, as discussed in (52). In the following analysis we will limit to consider a modulation function satisfying the constraint , where . This condition is satisfied, for instance, for a sinusoidal or square-wave modulation, under an appropriate choice of . In this case, for the Hermitian ac-driving problem the renormalization of the tunneling rate affects the amplitude (but not the phase) of the effective hopping rate (52).

The hyper-ballistic behavior of the system when driven close to the -symmetry breaking can be investigated by solving Eqs. (5) under initial excitation of an individual site (e.g. the site). In Fig. 4 we show, as an example, the numerically-computed evolution of obtained in a lattice with sites for the static Hamiltonian [panel (a)], for the dynamic Hamiltonian under driving at frequency and amplitude [panel (b)] corresponding to an Hermitian driving, and for driving at the same frequency and absolute value of the amplitude but in the non-Hermitian case, i.e. [panel (c)]. The square-wave modulation is shown in the inset of Fig. 4(b)-(c). In the non-Hermitian case, the norm of the vector state is bounded and does not diverge in the unbroken phase, however it is not conserved, as shown in the inset of Fig. 4(c) (53).

Most interestingly, the non-Hermitian driving results in a strong enhancement of the ballistic spreading in the lattice which can be quantitatively estimated by computing the variance of the probability distribution as a function of time [Fig. 4(d)], defined as , where is the average site occupation at time . It is seen that the ballistic spreading velocity is about times larger as compared to the undriven Hamiltonian (limited to ) for non-Hermitian driving close to the -symmetry breaking [cf. solid and dotted lines in Fig. 4(d)]. This enhancement factor is in agreement with an analytical estimation of the effective hopping rate in the driven lattice based on asymptotic analysis of Eq. (5). Actually, in the high frequency limit, the driven lattice described by Eq. (5) turns out to be equivalent to a static effective lattice with renormalized hopping rate given by (see e.g. (48)) (54):

(8) |

For a square-wave modulation, one readily obtains , i.e. for the parameters considered in our simulation. Therefore, the increase in the effective hopping rate corresponds to approximately the same increase in the ballistic rate. Note that Hermitian driving with similar parameters strongly reduces the effective hopping rate by a factor of [cf. dotted and dashed lines in Fig. 4(d)]. Indeed, in the Hermitian driving the ratio is always smaller than one.

Finally, we checked the strong reduction of wave packet spreading due to quantum diffraction resulting from the linearization of the quasi-energy spectrum close to the -symmetry breaking. To this aim, we numerically solved Eqs. (5) under initial excitation with a Gaussian probability distribution of full-width at of about sites and carrier momentum . The results for a lattice with sites are reported in Fig. 5(a) for the static Hamiltonian and in Fig. 5(b) for the dynamic Hamiltonian with a modulation frequency and , that is close to the -symmetry breaking. Normalized probability distributions at time are also shown in Fig. 5(c), aside with the initial wavepacket distribution (dashed line) for comparison. The spreading of the wavepacket turned out to be dramatically reduced [solid line in Fig. 5(c)] as compared to the spreading in the static lattice [dotted line in Fig. 5(c)].

## Iv Conclusions

We reported on a theoretical investigation of the spectral features and transport properties of a one-dimensional tight-binding superlattice under a time-periodic -symmetric modulation of the complex site energies. We observed a -symmetry phase transition of the system above a minimum value of the modulation frequency, provided that the amplitude of the modulation reaches a certain threshold value. The threshold turns out to increase almost linearly with the modulation frequency. Interestingly, the unbroken -symmetric phase of the driven system exhibits novel spectral features and transport behavior, with no counterpart in Hermitian time-periodic lattices. In particular, as the -symmetry breaking is approached, the quasi-energy spectrum is progressively broadened and linearized. As opposed to the ac-driven control scheme of Hermitian lattices, one of the main property of the complex ac modulation of the lattice is to enhance (rather than reduce or even suppress) the ballistic motion of the particle, giving rise to an hyper-ballistic regime, where the speed of propagation in the lattice in not limited by the hopping rate of the static lattice. Moreover, a nearly non-dispersive propagation regime is observed, as confirmed by direct numerical simulations.

We believe that our results can disclose novel opportunities for coherent control of transport in non-Hermitian lattices by exploiting the unique features of -symmetric Hamiltonians in the unbroken phase. In particular, it is envisioned that our approach can find a physical implementation in active optical systems based on synthetic photonic lattices or coupled fiber loops.

###### Acknowledgements.

This work was supported by the Fondazione Cariplo (Grant No. 2011-0338).## Appendix A Quasi-energy spectrum and the monodromy matrix

In this appendix we provide analytical derivation of the monodromy matrix for the system of Eqs. (5), leading to the quasi-energies of the time-periodic -symmetric superlattice. Let us first introduce two new complex variables, for the even-site and odd-site Wannier basis elements respectively:

(9a) | |||

(9b) |

Eqs. (5) are thus split into two sets of coupled equations:

(10a) | |||

(10b) |

Since Wannier basis and Bloch basis are related by Fourier transformation, we can expand the Wannier coefficients and in terms of the Bloch spectrum of the vector state defined by , with , i.e.:

(11a) | |||

(11b) |

Substitution of Eqs. (11) into Eqs. (10) yields the following set of coupled equations for the evolution of the Bloch spectral functions and :

(12a) | |||

(12b) |

When is periodic of period , for any given value of the quasi-momentum , Eqs. (12) represent a homogeneous system of ordinary differential equations with time-periodic coefficients. According to Floquet theory, the solution of this system has the general form (see for instance (55)):

(13) |

where is the initial value at time , is a time-periodic matrix of period with being the identity matrix, and is a time-independent matrix whose eigenvalues are the quasi-energies (so-called Floquet exponents) of the system. The real and imaginary parts of the quasi-energies can be computed from the eigenvalues of the monodromy matrix using Eqs.(7) given in the text. The monodromy matrix connects the solution to Eqs. (A4) over one oscillation cycle (i.e. from to ) and is given by . The columns of are the two solutions and of Eqs. (A4) obtained under initial conditions (at time ) and respectively. Since for a square wave modulation the coefficients in Eqs. (12) are constants over a semi-cycle, and can be analytically determined over each semi-cyle, thus factorizes into the product of the two matrices (one for each semi-cycle) and given by Eqs. (6).

### References

- N. Moiseyev, Non-Hermitian Quantum Mechanics (Cambridge University Press, London, Cambridge, 2011).
- N. Moiseyev, Phys. Rep. 302, 212 (1998); J.G. Muga, J.P. Palao, B. Navarro, and I.L. Egusquiza Phys. Rep. 395, 357 (2004); I. Rotter, J. Phys. A: Math. Theor. 42, 153001 (2009).
- V. Buslaev and V. Grecchi, J. Phys. A: Math. Gen. 26 5541 (1993).
- C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998).
- C. M. Bender, Rep. Prog. Phys. 70, 957 (2007).
- A. Mostafazadeh, J. Phys. A: Math. Gen. 36, 7081 (2003).
- As pointed out by Mostafazadeh [A. Mostafazadeh, J. Math. Phys. 43, 205 (2002)], PT-symmetric Hamiltonians are a special case of pseudo-Hermitian Hamiltonians, which are Hamiltonians that have an antilinear symmetry.
- C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. Lett. 93, 251601 (2004).
- A. M. Shalaby, Phys. Rev. D 80, 025006 (2009).
- C. M. Bender, S. Boettcher, and P. N. Meisinger, J. Math. Phys. 40, 2201 (1999).
- U. Guenther, F. Stefani, and M. Znojil, J. Math. Phys. 46, 063504 (2005); U. Guenther, B. F. Samsonov, and F. Stefani, J. Phys. A: Math. Theor. 40, F169 (2007).
- M. B. Plenio, P. L. Knight, Rev. Mod. Phys. 70, 101 (1998).
- J. Schindler, A. Li, M.C. Zheng, F.M. Ellis, and T. Kottos, Phys. Rev. A 84, 040101 (2011); H. Ramezani, J. Schindler, F. M. Ellis, U. Günther, and T. Kottos, Phys. Rev. A 85, 062122 (2012); Z. Lin, J. Schindler, F.M. Ellis, and T. Kottos, Phys. Rev. A 85, 050101 (2012).
- A. Ruschhaupt, F. Delgado, and J.G. Muga, J. Phys. A: Math. Gen. 38, L171 (2005).
- R. El-Ganainy, K.G. Makris, D. N. Christodoulides, and Z. H. Musslimani, Opt. Lett. 32, 2632 (2007).
- S. Klaiman, U. Günther, and N. Moiseyev, Phys. Rev. Lett. 101, 080402 (2008).
- K. G. Makris, R. El-Ganainy, D.N. Christodoulides, and Z. H. Musslimani, Phys. Rev. Lett. 100, 103904 (2008).
- S. Longhi, Phys. Rev. Lett. 103, 123601 (2009).
- A. Mostafazadeh, Phys. Rev. Lett. 102, 220402 (2009).
- A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. Christodoulides, Phys. Rev. Lett. 103, 093902 (2009).
- C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev and D. Kip, Nat. Phys. 6, 192 (2010).
- L. Feng, M. Ayache, J. Huang, Y.L. Xu, M.H. Lu, Y.F. Chen, Y. Fainman, and A. Scherer, Science 333, 729 (2011).
- S. V. Dmitriev, A. A. Sukhorukov, and Yu. S. Kivshar, Opt. Lett. 35, 2976(2010); A.A. Sukhorukov, Z. Xu, and Y.S. Kivshar, Phys. Rev. A 82, 043818 (2010); N. V. Alexeeva, I. V. Barashenkov, A.A. Sukhorukov, and Y.S. Kivshar, Phys. Rev. A 85, 063837 (2012).
- A. Regensburger, C. Bersch, M.-A. Miri, G. Onishchukov, D. N. Christodoulides, and U. Peschel, Nature 488, 167 (2012).
- L. Feng, Y.-L. Xu, W.S. Fegadolli, M.-H. Lu, J.E.B. Oliveira, V.R. Almeida, Y.-F. Chen, and A. Scherer, Nature Mat. (2012) doi:10.1038/nmat3495.
- H. Ramezani, T. Kottos, R. El-Ganainy, and D.N. Christodoulides, Phys. Rev. A 82, 043803 (2010).
- Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, Phys. Rev. Lett. 106, 213901 (2011).
- S. Longhi, Phys. Rev. A 81, 022102 (2010).
- S. Longhi, J. Phys. A: Math. Theor. 44, 485302 (2011).
- E.-M. Graefe and H. F. Jones, Phys. Rev. A 84, 013818 (2011).
- A. Mostafazadeh, Phys. Rev. A 87, 012103 (2013).
- S. Longhi, G. Della Valle and K. Staliunas, Phys. Rev. A 84, 042119 (2011).
- S. Longhi, Phys. Rev. A 82, 031801(R) (2010).
- Y. D. Chong, L. Ge, and A. D. Stone, Phys. Rev. Lett. 106, 093902 (2011).
- C. F. M. Faria, and A, Fring, J. Phys. A: Math. Gen. 39, 9269 (2006).
- C. F. M. Faria, and A. Fring, Laser. Phys. 17, 424 (2007).
- C.T. West, T. Kottos, and T. Prosen, Phys. Rev. Lett. 104, 054102 (2010).
- J. Wu, and X.-T. Xie, Phys. Rev. A 86, 032112 (2012).
- S. Longhi, Phys. Rev. B 80, 235102 (2009).
- N. Moiseyev, Phys. Rev. A 83, 052125 (2011).
- R. El-Ganainy, K.G. Makris, and D.N. Christodoulides, Phys. Rev. A 86, 033813 (2012).
- J. Gong, and Q.-H. Wang, arXiv:1210.5344v1 (2012).
- M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge U. Press, Cambridge, 1997).
- S. H. Autler, and C. H. Townes C. H., Phys. Rev. 100, 703 (1955).
- D. H. Dunlap and V. M. Kenkre, Phys. Rev. B 34, 3625 (1986).
- F. Grossmann, T. Dittrich, P. Jung, and P. Han̈ggi, Phys. Rev. Lett. 67, 516 (1991).
- G. Della Valle, M. Ornigotti, E. Cianci, V. Foglietti, P. Laporta, and S. Longhi, Phys. Rev. Lett. 98, 263601 (2007).
- M. Grifoni and P. Hänggi, Phys. Rep. 304, 229 (1998).
- Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. N. Christodoulides, J. Phys. A 41, 244 019 (2008).
- S. Longhi and K. Staliunas, Opt. Comm. 281, 4343 (2008).
- A. Szameit, Y. V. Kartashov, F. Dreisow, M. Heinrich, T. Pertsch, S. Nolte, A. Tünnermann, V. A. Vysloukh, F. Lederer, and L. Torner, Phys. Rev. Lett. 102, 153901 (2009).
- C. E. Creffield and F. Sols, Phys. Rev. Lett. 100, 250402 (2008).
- From our choice of the initial phase [cf. inset in Fig. 4(b)-(c)] and according to Eqs. (5), during the initial steps of the dynamics (more precisely, for ) even sites exhibit loss whereas odd sites exhibit gain. This results in a norm of the state vector for initial excitation of an even site [as shown in the inset of Fig. 4(c)] and for initial excitation of an odd site (not considered in our simulations).
- The equality of the two integrals on the right hand sides of Eq.(8) is ensured by the constraint we imposed on , that is , which for the square wave function (as well as for a sinusoidal modulation) is satisfied provided that .
- G. Teschl, Ordinary Differential Equations and Dynamical Systems (American Mathematical Society, Providence USA, 2010).