## Non-Hermitian spectral changes in the scattering of partially coherent radiation by periodic structures

P. A. BRANDÃO 1,* AND S. B. CAVALCANTI 1

1 Universidade Federal de Alagoas, Instituto de Física, 57072-900, Brazil * Corresponding author: paulo.brandao@fis.ufal.br

Compiled August 5, 2019

The physical aspects of partially coherent radiation interacting with deterministic non-Hermitian periodic materials remain largely unexplored in the statistical optics literature. Here, we consider the scattering of partially coherent radiation by a deterministic periodic medium, symmetric under the simultaneous transformations of parity inversion and time reversal, that is, a parity-time (PT)-symmetric periodic medium. Taking into account light fluctuations, one is able to describe the spectrum changes on propagation and the influence of the coherence-driven angular divergence effect. The far-field spectral density profile is found to depend crucially on the loss/gain properties of the material, giving rise to unexpected and contrasting spectral diffraction profiles when compared to the Hermitian ones. © 2019

Optical Society of America

http://dx.doi.org/10.1364/ao.XX.XXXXXX

A non-Hermitian formulation of quantum mechanics has been initiated by Bender and Boettcher in 1998 [1-5]. Nowadays, is has been well established, that a Hamiltonian which is invariant under the combined action of parity inversion and time reversal symmetries (PT symmetry), may have all-real valued eigenvalues. This is a necessary but not a sufficient condition to guarantee the reality of the spectrum [1]. In a series of interesting papers, Mostafazadeh has shown that operators having PT symmetry belong to the more general class of pseudo-Hermitian operators [6-8]. Regarding experimental results on PT-symmetric systems, optical physics has provided the most accessible route to verify the reality of eigenvalues of PT-symmetric operators along with the symmetry breaking transition point present in non-Hermitian configurations. This was possible due to the analogy between the Schrödinger and the paraxial wave equation [9-11], where the potential function is interpreted as the refractive index. In the last decades, non-Hermitian PT-optical lattices have been extensively investigated exhibiting exotic phenomena such as: non-unitary evolution [12], unidirectional invisibility [13], which was later experimentally evidenced in PT-symmetric syntethic photonic lattices [14] and in a PT-symmetric metamaterial [15]. Spectral singularities of Bragg scattering processes in periodic optical lattices have been associated with exceptional points, exhibiting a secular growth of plane waves that satisfy the Bragg condition, which saturates in case of wave packets with a broad momentum distribution [16, 17]. The effect of broadband quasi-PT symmetry in a finite frequency range on pulse propagation, has also been investigated in [18]. The interest in non-Hermitian structures has grown steadily and, in the last few years various aspects of scattering periodic structures were studied, such as the transmission properties of finite periodic structures [19] as well as resonant phenomena in layered periodic ones [20].

However, up to now all studies on PT-symmetric systems do not take into account fluctuations, could that be of the electromagnetic field or the medium, or both. As pointed out by Wolf [21], who studied the scattering of a partially coherent radiation field by a periodic medium, short wavelength radiation is not complete spatially coherent and therefore its random nature must be taken into account in diffraction studies. Some aspects of diffracted wavefields by non-Hermitian gratings were recently published in the literature [22-24]. However, to achieve some insight about the dynamical behavior of some observables such as the spectral density one must consider the coherence properties of the radiation field, particularly if one is dealing with short wavelength radiation. Within the limits of the Born approximation, the scattering of partially coherent radiation from a random medium, that is, a medium whose dielectric function is a random function of both position and time, has been reported to generate frequency shifts [25]. Later on, a study on the diffraction of partially coherent beams from three dimensional periodic structures has been addressed, demonstrating coherence-induced angular shifts of the diffraction orders due to the finite size of the coherence area [26]. Therefore, inspired by the elegant theory of coherence, developed mainly by Wolf, based on stationary random processes [27], and considering that PT-symmetric structures have been artificially obtained and have been the object of many studies worldwide, in the following we have studied the scattering of a partially coherent radiation field by a periodic PT-symmetric medium, within the first-order Born approximation. In contrast with the Hermitian potential, we find radical changes in the spectra profiles due to the combined effect of PT symmetry and partial coherence on propagation.

Let us begin by representing light as a stationary random process, characterized by the spectral degree of coherence µ ( ) i ( r 1 , r 2 ; ω ) [27]. Next, we suppose that it is diffracted by a periodic medium described by the function V ( r ) = V ( r + ) a ,

which is proportional to the refractive index of the material. After the interaction, in the asymptotic limit, the spectral density of the scattered radiation is given by [27]

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ S ^ { ( \infty ) } ( \mathbf r ; \omega ) & = \frac { S ^ { ( i ) } ( \omega ) } { r ^ { 2 } } \\ & \times \int \int _ { D } d ^ { 3 } r _ { 1 } d ^ { 3 } r _ { 2 } \mu ^ { ( i ) } ( \mathbf r _ { 1 }, \mathbf r _ { 2 } ; \omega ) V ^ { * } ( \mathbf r _ { 1 } ) \quad \text{$(1)}$} \quad \text{$whac} \\ & \times V ( \mathbf r _ { 2 } ) \exp [ - i k s \cdot ( \mathbf r _ { 2 } - \mathbf r _ { 1 } ) ], \quad \text{$speq} \cdot \\ \text{where $k = \omega / c$ with $c$ being the speed of light in vacuum and}$$

where k = ω / c with c being the speed of light in vacuum and the first-order Born approximation was used. The unit vector s is directed from the scatterer to the observation point. Let us assume a one-dimensional periodic medium described by V ( r ) = v x ( ) δ ( y ) δ ( z ) where v x ( ) is a periodic function with period a and δ denotes the Dirac delta function. Following [21], we also consider that the incident spectral degree of coherence has a Gaussian form given by

$$\mu ^ { ( i ) } ( \rho _ { 1 }, \rho _ { 2 } ; \omega ) = \exp \left [ - \frac { ( \rho _ { 2 } - \rho _ { 1 } ) ^ { 2 } } { 2 \sigma ^ { 2 } } \right ],$$

with σ denoting the range of transversal coherence with ( ρ 1 , ρ 2 ) representing points at the cross section of the beam. Let us now choose a one dimensional PT-symmetric material by requiring that ( PT ) V ( r )( PT ) -1 = V ( r ) which, in turn, implies that v x ( ) = v ∗ ( -x ) or, in words, the real (imaginary) part of v x ( ) must be even (odd) under the transformation x → -x . The operators P and T are the parity inversion ( r → -r ) and time reversal operators ( i → -i ), respectively. Thus, to satisfy these conditions, we consider the PT-symmetric medium described by

$$\mu, \quad v ( x ) = \frac { 1 } { 2 } + v _ { r } \cos \left ( \frac { 2 \pi x } { a } \right ) + i v _ { i } \sin \left ( \frac { 2 \pi x } { a } \right ), \quad \text{(3)}$$

where vr and vi are real and positive numbers which may depend on the frequency ω . The function v x ( ) is a truncated Fourier series of the form

$$v ( x ) = \sum _ { n = - \infty } ^ { \infty } c _ { n } \exp \left ( \frac { 2 \pi i x n } { a } \right ) _ { \begin{array} { c } \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots } \end{array} \quad ( 4 )$$

with coefficients cn given by c ± 1 = ( vr ± vi ) /2 and c 0 = 1/2, and cn = 0 for n ≥ 2. The coefficients cn are real numbers according to the condition v x ( ) = v ∗ ( -x ) which implies that cn = c ∗ n , as can be easily verified. In terms of cn , Eq. (1) can be recast into

$$S ^ { ( \infty ) } ( \mathbf r ; \omega ) & = \frac { L S ^ { ( i ) } ( \omega ) } { r ^ { 2 } } \sum _ { n = - \infty } ^ { \infty } | c _ { n } | ^ { 2 } \int _ { L } d x ^ { \prime } \exp \left [ - \frac { ( x ^ { \prime } ) ^ { 2 } } { 2 \sigma ^ { 2 } } \right ] \quad & \quad \text{(5)} \\ & \times \exp [ - i ( k \mathbf s \cdot \mathbf x ) x ^ { \prime } ] \exp \left ( \frac { 2 \pi i n x ^ { \prime } } { a } \right ), \quad & \quad \text{(33)} \\ & \quad. \quad \text{$\mathbf Fio$}.$$

where L is a characteristic (large) length of the structure, x ′ = x 2 -x 1 and x is the unit vector in the x direction. In the particular case of the material described in Eq. (3) this integral is readily evaluated in closed form

$$\text{readily evaluated in closed form} \, \stackrel { \cdot \, \cdot } { \sim } \, \stackrel { \cdot } { \sim } \, \text{ure 1} } { S ^ { ( \infty ) } ( \mathbf r ; \omega ) } & = \frac { L \sigma \sqrt { 2 \pi } S ^ { ( i ) } ( \omega ) } { 4 \mathbf r ^ { 2 } } \left \{ \exp \left [ \frac { - ( k \sigma \cos \theta ) ^ { 2 } } { 2 } \right ] } { 2 } \right ] \, \quad \text{(a) } \, \text{(m) } \, \text{and} \, \, \text{(inner)} } \, \text{$\frac{$-\sigma^{2} ( k \cos \theta - 2 \pi / a )^{2}} { 2 } \right ] \right \}.$$

where θ = arccos ( s · x ) . In what follows, we will assume that the spectral density of the incident wavefield has the following Gaussian profile:

$$S ^ { ( i ) } ( \omega ) = S _ { 0 } \exp \left [ - \frac { ( \omega - \omega _ { 0 } ) ^ { 2 } } { 2 \delta ^ { 2 } } \right ],$$

where ω 0 is the central frequency, δ the bandwidth and S 0 the spectrum amplitude.

Fig. 1. (Color online) Normalized scattered spectral density of partially coherent radiation with σ = 10 × ( 1/ k ) . (a) Hermitian scattering with vi = 0. Non-Hermitian scattering with (b) vi = 0.25 and (c) vi = 0.5. The real part of the potential function is vr = 0.5, the incident spectrum is described by Eq. (7) with ω 0 = 3.54 × 10 15 rad/s ( λ 0 = 532 nm), δ = 0.1 ω 0 and a = 2 λ . The white dashed lines are fixed at θ = ± π 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00645v1-figure-f75cf38b-0bb2-4fef-9727-63e2745888ee.jpg)

To obtain sensible numerical estimates from this model, we assume that the incident spectral density is centered at ω 0 = 3.54 × 10 15 rad/s ( λ 0 = 532 nm) with the bandwidth δ = 0.1 ω 0 . Figure 1 depicts the normalized spectral density ˜ S ( θ , ω ) = 4 r 2 S ( θ , ω ) / √ 2 πσ LS 0 in the ( θ , ω ) plane for an incident light with a high degree of coherence, σ = 10 × ( 1/ k ) . Part (a) of Figure 1 illustrates the Hermitian case where vi = 0. In this situation, we recover the results studied in [21] in which a symmetrical interference pattern is obtained with a central maximum at θ c = ± π /2, highlighted by the two white dashed lines, accompanied by two secondary maxima at angles θ ext = ± 2 π /3 and θ in = ± π /3. The subindex ext (in) indicate the outermost (innermost) maxima relative to θ = 0. The value vr = 0.5 is assumed in all plots. Notice that θ c corresponds to a direction perpendicular to the lattice. Therefore, in this situation, we ob-

tain the usual diffraction pattern present in crystallographic systems. If we now allow the lattice to absorb energy from and to give energy to the wavefield, by adjusting the non-Hermitian parameter to vi = 0.25, the result is shown in part (b) of Figure 1. In this case, the amplitudes of the secondary maxima change according to:

$$\text{ amplitudes at $\theta_{in} \sim (v_{r} + v_{i})^{2}$,} \\ \text{ amplitudes at $\theta_{ext} \sim (v_{r} - v_{i})^{2}$.} \\... \dots \dots.. \cdot \cdot \cdot \dots.. \cdot \dots.. \cdot \dots........................................................................................................................................................................................................$$

The central maximum at θ c remains unchanged as we activate the loss/gain properties of the lattice. The effect of the nonHermiticity in this model changes only the amplitude factors as shown in Eq. (8). Therefore, it is possible to adjust the spectral properties of the wavefield by tuning the loss/gain properties of the lattice. The situation shown in part (c) of Figure 1 was obtained by letting vi = vr . In this case, the secondary maxima at θ ext are zero. This scenario generates an asymmetrical spectral density in relation to the perpendicular direction of the lattice. In non-relativistic quantum mechanics described by a Schrödinger equation with potential Eq. (3), the situation vr = vi represents the symmetry breaking point where the eigenvalues become degenerate. In the present model, the physical effect of vr = vi is the disappearance of the spectral amplitude at θ ext whose energy is transferred to the secondary maxima located at θ in , while the central maximum at θ c remains unchanged. Furthermore, the enhanced maximum become as intense as the central one, so that the non-Hermitian spectrum becomes a twin pair of maxima, in both directions alike.

Let us now decrease the coherence length of the radiation field to σ = 5 × ( 1/ k ) . In this situation, the temporal and spatial variations of the wavefield observed in any two pair of points have smaller statistical correlations than the previous case. Figure 2 shows the normalized spectral density for the same values used in Figure 1, except for the coherence length which is now σ = 5 × ( 1/ k ) . For the Hermitian case, shown in part (a), the two secondary maxima disappear and we no longer obtain the Hermitian interference pattern, only two maxima of equal intensities, at θ = θ c . The situation is quite different when the lattice is equipped with gain and loss regions. Parts (b) and (c) of the same figure show the spectral density for the nonHermitian case. It can be seen from these plots that two maxima are still discernible, even for an incident radiation field in the low coherence regime. Part (c), in particular, represents a nonHermitian interference pattern induced by the complex lattice. Lattice properties are generally obtained by examining the interference pattern they generate and, since in this case this pattern depends on the Hermiticity of the lattice, this result could be useful for effective spectroscopy and grating diffraction research.

As pointed out before, in the particular situation where vi = vr , the eigenvalues of the Schrödinger equation with the potential described by Eq. (3) undergoes a symmetry breaking transition. In the context of optics and the paraxial wave equation, the propagation of a wavefield also exhibits the process of energy transfer between normal modes. However, above this point, that is, when vi &gt; vr , some of the eigenvalues may become complex and are characterized by an uncontrolled increase of the amplitude of the electric field during propagation [28]. In contrast, in the present model we see no indication of a symmetry breaking point characterized by a divergent behavior of the spectral density. It should be noted, however, that the present approach has little resemblance with the non-relativistic

Fig. 2. (Color online) Same as Figure 1 except for the wavefield, which now is assumed within the low-coherence regime with σ = 5 × ( 1/ k ) . (a) Hermitian and (b,c) complex lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00645v1-figure-0de3d174-df8a-4d9c-9cc7-5cc1c6403421.jpg)

Schrödinger equation such that we do not expect all characteristics of the former system to apply here. In particular, we are observing the wavefield far away from the material, working with an asymptotic solution. Therefore, our theory might be extended to include cases where vi /greatermuch vr = 0.5 and, indeed, in this limiting situation the spectral density is given by

$$\mathfrak { n } _ { \cdot } \quad \cdot \quad \sigma \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \text{$\intercal} \\ \mathfrak { n } _ { \cdot } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \text{$\intercal} \\ \mathfrak { n } _ { \cdot } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot. \quad \text{$\intercal} \\ \mathfrak { n } _ { \cdot } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \text{$\intercal} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd� \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdot \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \. \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \cdot \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \intercal \quad \mathfrak { n } _ { \cdot } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \delta \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \cdot \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,\quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \interce \quad \mathfrak { n } _ { \cdot } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \cdot \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, }$$

The interesting thing about this equation is that the central maximum at θ c is now negligible compared to the secondary ones which are, in fact, the only relevant quantities. This produces a symmetrical interference pattern for moderate values of σ . To see the difference between Hermitian and non-Hermitian lattices in this approximation, Figure 3 shows the normalized spectral density for vi = 0 in part (a) and for vi /greatermuch vr in part (b) within the low-coherence regime where σ = 5 × ( 1/ k ) . A clear non-Hermitian diffraction pattern with two secondary maxima is still visible even in the case of a partially coherent incident field, when compared to the Hermitian case. In particular, a comparison between the results presented in part (c) of Figure 2 and part (b) of Figure 3 indicates that a visible interference pattern can be obtained by controlling the loss/gain properties of the lattice with σ fixed. All these results suggest the possibility to tune and to improve the visibility of the spectral density,

by adjusting gain/loss regions of the complex lattice without changing the geometry of the system. We should also mention that there is a small shift in the central frequency of the scattered spectral density compared to the initial frequency ω 0 that is solely due to the intrinsic nature of the correlations existent in the wavefield [25]. This shift is, however, negligible within the range of parameters considered in the plots and therefore is not visible in these pictures. Since the current experimental research interest on non-Hermitian systems is to consider materials with modulated loss, it is important to verify what classical coherence theory says in this particular context. If one assumes that, for example, v x ( ) = vr + iv i sin 2 ( π x / a ) , then the material has no gain, only loss [29]. A simple calculation shows that the complex coefficients cn in the Fourier expansion are c 0 = vr + iv i /2 and c ± 1 = -iv i /4. It is easy to see that none of the above effects are realized in this pure lossy configuration, implying that gain may be a decisive factor for these effects to appear.

Fig. 3. (Color online) Normalized spectral density of the scattered wavefield in the regime where vi /greatermuch vr . (a) Hermitian lattice with vi = 0 and (b) non-Hermitian lattice with vi /greatermuch vr described by Eq. (9). The coherence length σ is the same as in Figure 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00645v1-figure-6faa0ee0-8a82-4810-bb32-5d2d11726831.jpg)

Summing up, considering the importance of the statistical fluctuations on the evolution of the spectra of short wavelength radiation, we have investigated the scattering of partially coherent electromagnetic fields by a periodic PT-symmetric medium. Assuming a stationary random process to represent the incident partially coherent field, and within the first-order Born approximation, we obtain the spectral density of the scattered radiation. Our results have provided us with new physical insights into a class of scattering process of a stationary partially coherent field from a deterministic PT-symmetric periodic medium. Experimentally, artificial optical structures with controlled gain and loss regions have been fully pursued, and one may study various types of potential functions as well, as various types of statistical processes. In this way our results should be quite useful in the understanding of the diffraction patterns obtained by the scattering of beams in any state of partial coherence from any type of PT-symmetric media, not necessarily a periodic medium. Work in this direction is in current progress and will be published elsewhere.

The authors are indebted to the Brazilian agencies CNPq and FAPEAL for partial financial support.

## REFERENCES

- 1. C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80 , 5243 (1998).
- 2. C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. Lett. 89 , 1 (2002).
- 3. C. M. Bender, S. Boettcher, and P . N. Meisinger, J. Math. Phys. 40 , 2201 (1999).
- 4. C. M. Bender, Reports on Prog. Phys. 70 , 947 (2007).
- 5. C. M. Bender, PT Symmetry: In Quantum and Classical Physics (World Scientific Publishing, 2018).
- 6. A. Mostafazadeh, J. Math. Phys. 43 , 205 (2002).
- 7. A. Mostafazadeh, J. Math. Phys. 43 , 2814 (2002).
- 8. A. Mostafazadeh, J. Math. Phys. 43 , 3944 (2002).
- 9. S. Longhi, Epl. 120 (2017).
- 10. A. Ruschhaupt, F. Delgado, and J. G. Muga, J. Phys. A: Math. Gen. 38 (2005).
- 11. C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, Nat. Phys. 6 , 192 (2010).
- 12. M. V. Berry, J. Phys. A: Math. Theor. 41 , 244007 (2008).
- 13. Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, Phys. Rev. Lett. 106 , 213901 (2011).

14.

A.

Regensburger,

C.

Bersch,

M.-A.

Miri,

G.

Onishchukov,

D.

N.

Christodoulides, and U. Peschel, Nature.

488

, 167 (2012).

- 15. L. Feng, Y.-L. Xu, W. S. Fegadolli, M.-H. Lu, J. E. B. Oliveira, V. R. Almeida, Y.-F. Chen, and A. Scherer, Nat. Mater. 12 , 108 (2013).
- 16. S. Longhi, Phys. Rev. A - At. Mol. Opt. Phys. 81 , 1 (2010).
- 17. E.-M. Graefe and H. F. Jones, Phys. Rev. A 84 , 013818 (2011).
- 18. D. M. Tsvetkov, V. A. Bushuev, V. V. Konotop, and B. I. Mantsyzov, Phys. Rev. A 98 , 053844 (2018).
- 19. V. Achilleos, Y. Aurégan, and V. Pagneux, Phys. Rev. Lett. 119 , 243904 (2017).
- 20. O. V. Shramkova and G. P . Tsironis, J. Opt. 18 , 105101 (2016).
- 21. E. Wolf, Opt. Lett. 38 , 4023 (2013).
- 22. F. Gao, Y.-M. Liu, X.-D. Tian, C.-L. Cui, and J.-H. Wu, Opt. Express 26 , 33818 (2018).
- 23. T. Shui, W. X. Y ang, S. Liu, L. Li, and Z. Zhu, Phys. Rev. A 97 , 33819 (2018).
- 24. X. Y. Zhu, Y. L. Xu, Y . Zou, X. C. Sun, C. He, M. H. Lu, X. P . Liu, and Y. F. Chen, Appl. Phys. Lett. 109 (2016).
- 25. E. Wolf, Phys. Rev. Lett. 63 , 2220 (1989).
- 26. M. Dušek, Phys. Rev. E 52 , 6833 (1995).
- 27. E. Wolf, Introduction to the Theory of Coherence and Polarization of Light (Cambridge University Press, 2007).
- 28. P. Brandão and S. Cavalcanti, Phys. Rev. A 96 (2017).
- 29. M. Berry and D. O'Dell, J. Phys. A: Math. Gen. 31 , 2093 (1998).

## FULL REFERENCES

- 1. C. M. Bender and S. Boettcher, 'Real spectra in non-hermitian hamiltonians having PT symmetry,' Phys. Rev. Lett. 80 , 5243-5246 (1998).
- 2. C. M. Bender, D. C. Brody, and H. F. Jones, 'Complex Extension of Quantum Mechanics,' Phys. Rev. Lett. 89 , 1-4 (2002).
- 3. C. M. Bender, S. Boettcher, and P . N. Meisinger, 'PT-symmetric quantum mechanics,' J. Math. Phys. 40 , 2201-2229 (1999).
- 4. C. M. Bender, 'Making sense of non-Hermitian Hamiltonians,' Reports on Prog. Phys. 70 , 947-1018 (2007).
- 5. C. M. Bender, PT Symmetry: In Quantum and Classical Physics (World Scientific Publishing, 2018).
- 6. A. Mostafazadeh, 'Pseudo-Hermiticity versus PT symmetry: The necessary condition for the reality of the spectrum of a non-Hermitian Hamiltonian,' J. Math. Phys. 43 , 205-214 (2002).
- 7. A. Mostafazadeh, 'Pseudo-Hermiticity versus PT-symmetry. II: A complete characterization of non-Hermitian Hamiltonians with a real spectrum,' J. Math. Phys. 43 , 2814-2816 (2002).
- 8. A. Mostafazadeh, 'Pseudo-Hermiticity versus PT-symmetry III: Equivalence of pseudo-Hermiticity and the presence of antilinear symmetries,' J. Math. Phys. 43 , 3944-3951 (2002).
- 9. S. Longhi, 'Parity-time symmetry meets photonics: A new twist in nonHermitian optics,' Epl. 120 (2017).
- 10. A. Ruschhaupt, F. Delgado, and J. G. Muga, 'Physical realization of PT-symmetric potential scattering in a planar slab waveguide,' J. Phys. A: Math. Gen. 38 (2005).
- 11. C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, 'Observation of parity-time symmetry in optics,' Nat. Phys. 6 , 192-195 (2010).
- 12. M. V. Berry, 'Optical lattices with PT symmetry are not transparent,' J. Phys. A: Math. Theor. 41 , 244007 (2008).
- 13. Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, 'Unidirectional invisibility induced by PT -symmetric periodic structures,' Phys. Rev. Lett. 106 , 213901 (2011).
- 14. A. Regensburger, C. Bersch, M.-A. Miri, G. Onishchukov, D. N. Christodoulides, and U. Peschel, 'Parity-time synthetic photonic lattices,' Nature. 488 , 167 (2012).
- 15. L. Feng, Y.-L. Xu, W. S. Fegadolli, M.-H. Lu, J. E. B. Oliveira, V. R. Almeida, Y.-F. Chen, and A. Scherer, 'Experimental demonstration of a unidirectional reflectionless parity-time metamaterial at optical frequencies,' Nat. Mater. 12 , 108-113 (2013).
- 16. S. Longhi, 'Spectral singularities and Bragg scattering in complex crystals,' Phys. Rev. A - At. Mol. Opt. Phys. 81 , 1-6 (2010).
- 17. E.-M. Graefe and H. F. Jones, ' PT -symmetric sinusoidal optical lattices at the symmetry- breaking threshold,' Phys. Rev. A 84 , 013818 (2011).
- 18. D. M. Tsvetkov, V. A. Bushuev, V. V. Konotop, and B. I. Mantsyzov, 'Broadband quasiPT -symmetry sustained by inhomogeneous broadening of the spectral line,' Phys. Rev. A 98 , 053844 (2018).
- 19. V. Achilleos, Y . Aurégan, and V. Pagneux, 'Scattering by finite periodic PT -symmetric structures,' Phys. Rev. Lett. 119 , 243904 (2017).
- 20. O. V. Shramkova and G. P . Tsironis, 'Scattering properties ofPT- symmetric layered periodic structures,' J. Opt. 18 , 105101 (2016).
- 21. E. Wolf, 'Diffraction of radiation of any state of spatial coherence on media with periodic structure,' Opt. Lett. 38 , 4023 (2013).
- 22. F. Gao, Y.-M. Liu, X.-D. Tian, C.-L. Cui, and J.-H. Wu, 'Intrinsic link of asymmetric reflection and diffraction in non-Hermitian gratings,' Opt. Express 26 , 33818 (2018).
- 23. T. Shui, W. X. Y ang, S. Liu, L. Li, and Z. Zhu, 'Asymmetric diffraction by atomic gratings with optical PT symmetry in the Raman-Nath regime,' Phys. Rev. A 97 , 33819 (2018).
- 24. X. Y. Zhu, Y . L. Xu, Y . Zou, X. C. Sun, C. He, M. H. Lu, X. P . Liu, and Y . F . Chen, 'Asymmetric diffraction based on a passive parity-time grating,' Appl. Phys. Lett. 109 (2016).
- 25. E. Wolf, 'Correlation-induced doppler-type frequency shifts of spectral lines,' Phys. Rev. Lett. 63 , 2220-2223 (1989).
- 26. M. Dušek, 'Diffraction of partially coherent beams on threedimensional periodic structures and the angular shifts of the diffraction maxima,' Phys. Rev. E 52 , 6833-6840 (1995).
- 27. E. Wolf, Introduction to the Theory of Coherence and Polarization of Light (Cambridge University Press, 2007).
- 28. P. Brandão and S. Cavalcanti, 'Bragg-induced power oscillations in PT -symmetric periodic photonic structures,' Phys. Rev. A 96 (2017).
- 29. M. Berry and D. O'Dell, 'Diffraction by volume gratings with imaginary potentials,' J. Phys. A: Math. Gen. 31 , 2093 (1998).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00645v1-figure-efd47455-eed6-465d-bfc2-ce2a99c28038.jpg)
Paulo A. Brandão S. B. Cavalcanti
Non-Hermitian spectral changes in the scattering of partially coherent radiation by periodic structures
The physical aspects of partially coherent radiation interacting with
deterministic non-Hermitian periodic materials remain largely unexplored in the
statistical optics literature. Here, we consider the scattering of partially
coherent radiation by a deterministic periodic medium, symmetric under the
simultaneous transformations of parity inversion and time reversal, that is, a
parity-time (PT)-symmetric periodic medium. Taking into account light
fluctuations, one is able to describe the spectrum changes on propagation and
the influence of the coherence-driven angular divergence effect. The far-field
spectral density profile is found to depend crucially on the loss/gain
properties of the material, giving rise to unexpected and contrasting spectral
diffraction profiles when compared to the Hermitian ones.