arXiv:1612.01493v4 [physics.class-ph] 27 Apr 2017 ∗ † [email protected] [email protected] xc ouin ote(+)dmninlnnierMaxwel nonlinear (1+1)-dimensional the to solutions Exact ewrs olna awl qain,Hdgahtransfo Hodograph waves of Equations, superposition Maxwell Nonlinear wa of Keywords: in superposition phenomena, nonlinear nonlinear solution and some oscillations, exact nonlinear discuss the we Through exac equations, . the Maxwell arbitrari without find in p media to propagation we neous waves issue paper, electromagnetic unresolved this of an In problem still coordinates. is curvilinear orthogonal it perspec the fact, practical In and theoretical cated. both from interest great hrceiigeetoantcwv rpgto nnonli in propagation wave electromagnetic Characterizing eateto ahmtc,ScunUiest 105 Che 610065, University Sichuan Mathematics, of Department in u ioZag ah Zhao, Dazhi Zhang, Xiao Hu, Liang rhgnlcriiercoordinates curvilinear orthogonal Dtd pi 8 2017) 28, April (Dated: ves. mto,Clnrclsokwvs Nonlinear waves, shock Cylindrical rmation, ie,ee huhi seteeycompli- extremely is it though even tives, ∗ ynnieradpriual inhomoge- particularly and nonlinear ly eeta nltcmto ohnl the handle to method analytic an resent n aKn Luo MaoKang and eradihmgnosmdai of is media inhomogeneous and near fte(+)dmninlnonlinear (1+1)-dimensional the of s ouin otennierwvsin waves nonlinear the to solutions t ldn yidia hc ae,free waves, shock cylindrical cluding gu iha,China Sichuan, ngdu, qain nthe in equations l † 2

I. INTRODUCTION

Along with the discovery of media possessing nonlinear electromagnetic properties, understanding electromagnetic wave propagation in nonlinear media has become a topical problem in physics [1–4]. Given the fact that most of media in nature or artificial are nonlinear and inhomogeneous, it is of special importance to study the behavior of electromagnetic fields in media with these two properties. However, this problem is yet to be solved [5, 6], especially for exact solutions which are essential in understanding the related physical processes and in developing new computational asymptotic methods [7–15]. Traditional methods like the inverse scattering method and B¨acklund transformation are efficient for constructing exaction solutions in the theory of plane waves [16–18], but they turn out to be inapplicable in cylindrical and spherical waves [19]. Meanwhile, finding new and physically meaningful exact solution of nonlinear waves is a highly urgent task, due to its own limitation of coupled-wave equation method [20] and difficulty to acquire the satisfactory description of strongly nonlinear phenomenons (i.e. self-steepening and shock waves) by . Notably, Ref. [5] reports that they can obtain an exact solution to the cylindrical electromagnetic wave propagation problem in a medium which has a special form of nonlinearity - with the function being an exponential function of electric field. Following their lead, a number of related works have been reported recently [6, 19–28]. To the best of our knowledge, however, serious difficulties have been met in studying other nonlinear forms. Here, we propose a method to exactly solve the (1+1)-dimensional Maxwell equations in arbitrarily nonlinear and particularly inhomogeneous media without dispersion. We show that if exact solutions of a linear system are given, then exact solutions of the corresponding two nonlinear systems can also be obtained; and the exact solutions of a can be got from either of the corresponding two linear systems. By usage of this connection, we can investigate some nonlinear phenomena by analytic methods. As illustrations, we analyze cylindrical shock waves and nonlinear oscillations of cylindrical waves. We also show how to deal with two waves propagation in a nonlinear medium.

II. ELECTROMAGNETIC MODEL

Consider electromagnetic fields in nonmagnetic and nondispersive media in the system of orthogonal curvilinear coordinates (u,v,w). If the fields are independent of v and w, and the direction of the electromagnetic waves is along the v axis, viz. v-axis polarization, then the (1+1)-dimensional Maxwell equations can be written as

1 ∂Ev(u,t) 1 ∂lv ∂Hw(u,t) + Ev(u,t)= −µ0 , lu ∂u lulv ∂u ∂t 1 ∂Hw(u,t) 1 ∂lw ∂Ev(u,t) − − Hw(u,t)= ǫ(u, Ev) , lu ∂u lulw ∂u ∂t where ǫ(u, Ev) = dDv/dEv is a dielectric function, and lu, lv and lw are the Lam´ecoefficients [29]. When selecting cylindrical coordinate system (r, φ, z), taking r axis as u axis and considering z-axis polarization, these equations are reduced into a model of cylindrical waves which are studied in Refs. [5, 6, 19–28]. For convenience, we introduce the dimensionless variables ρ = u/u0, τ = ct/u0, E = Ev/(N/C), and H = µ0cHw/(N/C), where u0 is a constant associated with the characteristic spatial scale, c is the speed of in vacuum, and (N/C) is the unit of electric field strength in the international system of units. Then we can get the dimensionless Maxwell equations as follows: ∂E(ρ, τ) ∂H(ρ, τ) g (ρ) + g (ρ)E(ρ, τ)= , 1 ∂ρ 2 ∂τ (1) ∂H(ρ, τ) ∂E(ρ, τ) g (ρ) + g (ρ)H(ρ, τ)= ǫ(ρ, E) , 3 ∂ρ 4 ∂τ where g1(ρ), g2(ρ), g3(ρ) and g4(ρ) are the dimensionless coefficients. Hereafter, dielectric function ǫ(ρ, E) is used to represent the media. Define

−1 g5(ρ) = exp (g1(ρ)) g2(ρ)dρ , Z 

−1 g6(ρ) = exp (g3(ρ)) g4(ρ)dρ , Z  2 g1(ρ)[g6(ρ)] g5(ρ) f(ρ)= 2 ,h = dρ = h(ρ), g3(ρ)[g5(ρ)] Z g1(ρ)g6(ρ) 3 where hereafter the integral constants in the indefinite integrals are always set as zero. We use the following ansatz:

ξ = αg5(ρ)E + βh + γ, η = αg6(ρ)H + βτ + γ, (2) where α, β and γ are arbitrary constants (β 6= 0).

III. METHODS

In contrast with the previous studies starting with a nonlinear system directly [5, 6, 19–28], our work begins with a linear system and then transforms it into a nonlinear system through a hodograph transformation [30–32]. Because there exists some mathematical difficulties in applying the hodograph transformation on nonlinear system directly. a. Firstly We start with a linear system that characterizes the electromagnetic fields in vacuum as follows: ∂E ∂H ∂H ∂E g (ρ) + g (ρ)E = ,g (ρ) + g (ρ)H = , (3) 1 ∂ρ 2 ∂τ 3 ∂ρ 4 ∂τ which is derived from system (1) by choosing ǫ(ρ, E) = 1. Suppose E0 = E0(ρ, τ) and H0 = H0(ρ, τ) is a solution of Eqs. (3). Then by ansatz (2), Eqs. (3) can be reduced to

∂ξ ∂η ∂η ∂ξ = , = f(ρ) . (4) ∂h ∂τ ∂h ∂τ A solution of Eqs. (4) is

ξ = C1g5(ρ)E0 + C2h + C3, η = C1g6(ρ)H0 + C2τ + C3, (5) where C1, C2 and C3 are arbitrary constants. If the Jacobian of Eqs. (4) D(ξ, η)/D(h, τ) = 0, then their solutions must have the form ξ = F1(τ ± f(ρ)dh) and ′ ′ η = F2(τ ± f(ρ)dh), where F1 and F2 are arbitrary functions which satisfy F2(τ ± f(ρ)dh) =R±p f(ρ)F1(τ ± ′ ′ f(ρ)dh)R.p (Here the prime denotes the derivative.) In particular, f(ρ) = ±F2/FR1p= F (τ ± fp(ρ)dh) holds Ronlyp if f(ρ) is constant. Since our intention is transforming Eqs. (4)p into some nonlinear forms, weR p do not need to consider this trivial case. When the Jacobian D(ξ, η)/D(h, τ) 6= 0, the hodograph transformation is applicable, and hence we can view ξ and η as independent variables and obtain ∂h ∂τ ∂τ ∂h = , = f(ρ) . (6) ∂ξ ∂η ∂ξ ∂η

Here formula (5) is still valid. Now we look at system (1) for the case

ǫ(ρ, E)= f(σ)/f(ρ), (7)

−1 −1 −1 where σ = h (αg5(ρ)E + βh + γ) = h (ξ) and h (ξ) is the inverse function of h(ξ). Through ansatz (2) we get ∂ξ ∂η ∂η ∂ξ = , = f(σ) , ∂h ∂τ ∂h ∂τ whose dependent variables (ξ and η) and independent variables (h and τ) share the same relation as nonlinear system (6). By analogy to formula (5) it yields

h =C1g5(σ)E0(σ, η)+ C2ξ + C3, (8) τ =C1g6(σ)H0(σ, η)+ C2η + C3.

Expressing E and H from ansatz (2), using formula (8) and setting C1 = −α/β, C2 =1/β and C3 = −γ/β, we obtain

g5(σ) g6(σ) E(ρ, τ)= E0(σ, η),H(ρ, τ)= H0(σ, η). (9) g5(ρ) g6(ρ)

These expressions give an exact solution of Maxwell equations (1) in inhomogeneous nonlinear and nondispersive media (7). 4

Typically, the Maxwell equations have a set of boundary and initial conditions. For the preassigned conditions of Eqs. (1) and (7), it is not difficult to get the corresponding conditions of linear Eqs. (3). For example, if Eqs. (1) and (7) meet the initial conditions

E(ρ, 0) = ϕ1(ρ),H(ρ, 0) = ϕ2(ρ), then the solution of Eqs. (3) must meet the generalized initial conditions

ϕ1(ρ)g5(ρ) ϕ2(ρ)g6(ρ) E0(σ0, η0)= ,H0(σ0, η0)= , g5(σ0) g6(σ0)

−1 where σ0 = h (αg5(ρ)ϕ1(ρ)+ βh(ρ)+ γ) and η0 = αg6(ρ)ϕ2(ρ)+ γ. b. Secondly Substituting the following ansatz into linear system (3): g (ρ) h¯ = 6 dρ = h¯(ρ), Z g3(ρ)g5(ρ)

ξ¯ = αg5(ρ)E + βτ + γ, η¯ = αg6(ρ)H + βh¯ + γ, we get ∂ξ¯ 1 ∂η¯ ∂η¯ ∂ξ¯ = , = . (10) ∂h¯ f(ρ) ∂τ ∂h¯ ∂τ

Similarly, assuming E1 = E1(ρ, τ) and H1 = H1(ρ, τ) is a solution of system (3), then a solution of system (10) is α 1 γ α 1 γ ξ¯ = − g (ρ)E + τ − , η¯ = − g (ρ)H + h¯ − . β 5 1 β β β 6 1 β β By applying the hodograph transformation to linear system (10), we get the nonlinear system

∂τ 1 ∂h¯ ∂τ ∂h¯ = , = . (11) ∂η¯ f(ρ) ∂ξ¯ ∂ξ¯ ∂η¯ When the function ǫ is chosen such that

ǫ(ρ, E)=1/[f(ρ)f(¯σ)], (12)

−1 −1 −1 whereσ ¯ = h¯ (αg5(ρ)E + βh + γ) = h¯ (ξ), then system (1) share the same form as system (11): ∂hη = [f(¯σ)] ∂τ ξ, ∂hξ = ∂τ η, by ansatz (2). Hence τ and h satisfy α 1 γ τ = − g (¯σ)E (¯σ, η)+ η − , β 5 1 β β (13) α 1 γ h = − g (¯σ)H (¯σ, η)+ ξ − . β 6 1 β β Substituting ansatz (2) into formula (13), we can write an exact solution of the Maxwell equations in inhomogeneous nonlinear and nondispersive media (12) as:

g6(¯σ) g5(¯σ) E(ρ, τ)= H1(¯σ, η),H(ρ, τ)= E1(¯σ, η). (14) g5(ρ) g6(ρ) c. Thirdly We consider a more generic linear system arising from (1), in which ǫ satisfies:

ǫ(ρ, E)= g(h(ρ))/f(ρ), (15) with g being an arbitrary function of h(ρ). If we define

g˜3(ρ)= g3(ρ)f(ρ)/g(h(ρ)), g˜4(ρ)= g4(ρ)f(ρ)/g(h(ρ)), then the linear system becomes ∂E ∂H ∂H ∂E g (ρ) + g (ρ)E = , g˜ (ρ) +˜g (ρ)H = , (16) 1 ∂ρ 2 ∂τ 3 ∂ρ 4 ∂τ 5 which is a case we have just discussed. Hence an exact solution of the following nonlinear system ∂E ∂H g (ρ) + g (ρ)E = , 1 ∂ρ 2 ∂τ ∂H g(ξ) ∂E g˜ (ρ) +˜g (ρ)H = , 3 ∂ρ 4 g(h(ρ)) ∂τ is given by formula (9), where E0(ρ, τ) and H0(ρ, τ) represent a solution of system (16). In other words, if E0 and H0 is a solution of Maxwell equations (1) in linear media (15), then an exact solution of Eqs. (1) in the inhomogeneous nonlinear and nondispersive media

ǫ(ρ, E)= g(ξ)/f(ρ)= g(αg5(ρ)E + βh + γ)/f(ρ), (17) is Eqs. (9). −1 Similarly, defining G(ρ) = g(ρ)dρ andσ ˜ = h (G(ξ)), and assuming E1 = E1(ρ, τ) and H1 = H1(ρ, τ) being a solution of Eqs. (1) in the linearR media ǫ(ρ, E)=1/[f(ρ)g(G−1(h(ρ)))], (18) we can also find an exact solution of the Maxwell equations in media (17) according to formula (14):

g6(˜σ) g5(˜σ) E(ρ, τ)= H1(˜σ, η),H(ρ, τ)= E1(˜σ, η). (19) g5(ρ) g6(ρ)

IV. RESULTS AND DISCUSSION

A. Results

Up to now, we have shown that from the solution of the Maxwell equations (1) in linear and homogeneous media ǫ(ρ, E) = 1, one can obtain the exact solution in nonlinear media (7) through formula (9) and in nonlinear media (12) through formula (14). And the exact solution of nonlinear system (1) and (17) can be got from linear media (15) by formula (9), or from linear media (18) by formula (19). Since g in media (17) is an arbitrary function, we have obtained the exact solutions of electromagnetic wave in arbitrarily nonlinear and particularly inhomogeneous media without dispersion, using the solutions of either of the corresponding two linear systems. For media (17), α and β embody the intensity of the nonlinear and inhomogeneous factor, respectively. The larger value |α| means stronger nonlinearity. In practice, by using different substrate materials and designing different cellular architectures, the opens a door to realize all possible material properties [33–39]. For Eqs. (1) and media (17) with the given definite-solution conditions, we can choose a simpler linear media from (15) and (18), as formulas (9) and (19) return the same result. Taking α = 0, β = 1 and γ = 0, the following equations represent the relationship between linear media (15) and linear media (18):

−1 g6(h (G(h))) −1 E0(ρ, τ)= H1(h (G(h)), τ), g5(ρ) −1 (20) g5(h (G(h))) −1 H0(ρ, τ)= E1(h (G(h)), τ), g6(ρ) which can be verified by straightforward differentiation. Taking β = 1 and γ = 0, with the weakening of nonlinearity, in the limit α → 0, the solution E and H is reduced into the solution E0 and H0 which corresponds to linear media (15). When selecting cylindrical coordinate system (r, φ, z), taking r axis as u axis and considering z-axis polarization, 1 2 we have g1(ρ) = 1, g2(ρ) = 0, g3(ρ) = 1, g4(ρ) = ρ , g5(ρ) = 1, g6(ρ) = ρ, h(ρ) = ln ρ and f(ρ) = ρ . If taking g so that media (15) satisfy g(h(ρ))/f(ρ) = 1, then from formula (9) we get

β αE E(ρ, τ)= E0(ρ e ,αρH + βτ), β−1 αE β αE H(ρ, τ)= ρ e H0(ρ e ,αρH + βτ), which represents the exact solution of the Maxwell equations in nonlinear inhomogeneous media ρ2(β−1) exp(2αE). This is the case discussed in Refs. [5, 6, 19–28]. 6

When selecting cylindrical coordinate system (r, φ, z), taking r axis as u axis and considering φ-axis polarization, 2 1 ρ 1 we have g1(ρ) = −1, g2(ρ) = − ρ , g3(ρ) = −1, g4(ρ) = 0, g5(ρ) = ρ, g6(ρ) = 1, h(ρ) = − 2 and f(ρ) = ρ2 . In this case, the specific expression of formula (9) is

βρ2 − 2αρE E(ρ, τ)= E ( βρ2 − 2αρE, αH + βτ), p 0 ρ p (21) 2 H(ρ, τ)=H0( βρ − 2αρE, αH + βτ), p and the specific expression of formula (19) is

1 β 2 2 ρ −αρE E(ρ, τ)= H1(e , αH + βτ), ρ (22) β 2 β 2 2 ρ −αρE 2 ρ −αρE H(ρ, τ)=e E1(e , αH + βτ).

Note that the φ-axis polarized and the z-axis polarized [5, 6] cylindrical waves are asymmetric in nonlinear media. To the best of our knowledge, the features of φ-axis polarized cylindrical waves in nonlinear media remain unknown. Now, let us discuss some applications of those results.

B. Initial value problem and cylindrical shock waves

Consider a φ-axis polarized cylindrical wave in vacuum with the following form: ∂E E ∂H ∂H ∂E − − = , − = ,ρ> 0, (23a) ∂ρ ρ ∂τ ∂ρ ∂τ

E(ρ, 0)=0,H(ρ, 0) = ϕ(ρ)=1/ 1+ ρ2, (23b) p applying the Hankel transform, we can represent the solution as ∞ E1(ρ, τ)= sJ1(sρ)Ψ(s) sin(sτ)ds, Z0 ∞ (24) H1(ρ, τ)= sJ0(sρ)Ψ(s)cos(sτ)ds, Z0 ∞ where Jν represents the first kind of the order ν of the Bessel function [40] and Ψ(s)= ρJ0(ρs)ϕ(ρ)dρ. According 0 2 ρ R 2 βρ −2αρE to Eq. (7) and Eq. (12), the corresponding two nonlinear media are ǫ(ρ, E)= βρ−2αE and ǫ(ρ, E)= ρ e . As mentioned earlier, α and β embody the intensity of the nonlinear and inhomogeneous factors, respectively. Ferroelectric bulk materials, ferroelectric thin films and superlattices are typical examples of such a medium [41–43]. Note that ferroelectric materials are known to be of great interest for many promising applications [44]. After substituting Eqs. (24) into formula (22), the exact solution of the initial value problem

∂E E ∂H ∂H 2 ∂E − − = , − = ρ2eβρ −2αρE ,ρ> 0, ∂ρ ρ ∂τ ∂ρ ∂τ β 2 0 < β,E(ρ, 0) = ϕ(e 2 ρ −αρE(ρ,0))/ρ,H(ρ, 0)=0, can be written as ∞ s β 2 2 ρ −αρE E(ρ, τ)= J0(se )Ψ(s) cos[s(αH + βτ)]ds Z0 ρ 1 1 = Re{ }, ρ (1 − i(αH + βτ))2 + eβρ2−2αρE ∞ (25) β 2 p β 2 2 ρ −αρE 2 ρ −αρE H(ρ, τ)=e sJ1(se )Ψ(s) sin[s(αH + βτ)]ds Z0 i + (αH + βτ) =Re{ }. (1 − i(αH + βτ))2 + eβρ2−2αρE p 7

=0 =5 1.5 0.2 =-0.5 =-0.5 =0 =0 0.75 0.1 Electric fields Electric fields

0 0 0.5 1.5 2.5 3.5 0.5 1.5 2.5 3.5

=10 *=11.5 0.2 0.2 =-0.5 =-0.5 =0 =0 0.1 0.1

Electric fields Electric fields *=2.15 0 0 0.5 1.5 2.5 3.5 0.5 1.5 2.5 3.5

FIG. 1. Linear (α = 0, dash-dotted curves) and nonlinear (α = −0.5, solid curves) electric fields of cylindrical wave (25) as function of the coordinate ρ at various times.

Here i is the imaginary unit. Note that the function ϕ describing the initial state of the magnetic field in vacuum can be chosen arbitrarily. One should only ensure the convergence of the integrals in Eqs. (24). Figure 1 shows the results of numerical calculation of E using formula (25) in the cases with α = 0, β = 1 and α = −0.5, β = 1. At various τ, the electric fields as function of the coordinate ρ in linear inhomogeneous media 2 2 ǫ(ρ, E) = ρ2eρ and nonlinear inhomogeneous media ǫ(ρ, E) = ρ2eρ +ρE are shown by the dash-dotted and solid curves, respectively. It is worth mentioning that the velocity of the peak of the nonlinear electric field is slower than that of the trailing edge, and thus leads to the trailing edge becomes increasingly steeper with increasing τ, which is called self-steepening [45, 46] and eventually creates a cylindrical electromagnetic shock wave at a point ρ = ρ∗ at time τ = τ ∗. Such phenomena cannot be described satisfactorily by perturbation theory. In the absence of dispersion, the solution described by Eqs. (25) becomes unsuitable when τ > τ ∗ because multiple different values of E and H satisfying Eqs. (25) correspond to the same value of ρ, leading to discontinuities of wave components.

C. Boundary value problem and free nonlinear oscillations

Now we consider the situation that the linear wave Eq. (23a) satisfies boundary conditions

1

E1(ρ, τ)= A[Y0(ωna)J1(ωnρ) − J0(ωna)Y1(ωnρ)] sin(ωnτ),

H1(ρ, τ)= A[Y0(ωna)J0(ωnρ) − J0(ωna)Y0(ωnρ)] cos(ωnτ), where Jν and Yν represent the first and the second kind of the order ν of the Bessel function respectively, ωn is the nth positive root of the equation J0(aω)Y0(bω) − J0(bω)Y0(aω) = 0, and A is an amplitude factor. According to formula (22), the exact solution of the following boundary value problem

∂E E ∂H ∂H 2 ∂E − − = , − = ρ2eβρ −2αρE , 0 < β, ∂ρ ρ ∂τ ∂ρ ∂τ 2 1 2 1 2 1 2 1 ( ln a) 2 ≤ ρ ≤ ( ln b) 2 , E(( ln a) 2 , 0) = E(( ln b) 2 , 0)=0, β β β β 8 which describes the corresponding nonlinear oscillations of cylindrical electromagnetic waves, can be written as

β 2 β 2 2 ρ −αρE 2 ρ −αρE E =[Y0(ωna)J0(ωne ) − J0(ωna)Y0(ωne )] A × cos[ωn(αH + βτ)], ρ (26) β 2 β 2 2 ρ −αρE 2 ρ −αρE H =[Y0(ωna)J1(ωne ) − J0(ωna)Y1(ωne )] β 2 2 ρ −αρE × Ae sin[ωn(αH + βτ)].

=1.1 4

0 E, H

-4 0 20 40 60 80 100

=1.5 4

0 E, H

-4 0 20 40 60 80 100

FIG. 2. Electric (solid curves) and magnetic (dash-dotted curves) fields of cylindrical wave (26) as functions of τ in the n = 1 . mode by taking α = 0.3, β = 1, a = e0 5, b = e2 and A = 3.

Figure 2 shows the results of numerical calculation of the fields using solution (26) in the case with α =0.3, β = 1, 0.5 2 a = e , b = e and A = 3 in the n = 1 mode (ω1 ≈ 0.5336) at two points ρ = 1.1 and ρ = 1.5. We see that the electric field and the magnetic field fluctuate with the same ωn/(2π) at ρ =1.1. However, the electric field varies at frequency ωn/(2π) while the magnetic field at the second harmonic ωn/π at ρ = 1.5, i.e., the second harmonic component which is one of the most intensively studied effects in nonlinear is pronounced. Similar to Ref. [5], it should be noted that the electric field E and magnetic field H determined by (26) become ambiguous when the mode n ≥ n∗, where threshold n∗ is an integer depending on the parameter α, β, a, b and A (e.g., n∗ = 2 for α =0.1, β = 1, a = 1, b = 2 and A = 1). Due to such phenomena are not physically admissible, the cylindrical standing waves solution in form (26) obtained without allowance for dispersion ceases to be suitable for the higher modes.

D. and sum- and difference-frequency generation

Consider two φ-axis polarized cylindrical waves with ̟1 and ̟2 propagating in the semi infinite inho- ρ mogeneous and nonlinear medium ǫ(ρ, E)= βρ−2αE :

∂E E ∂H ∂H ρ ∂E − − = , − = , ∂ρ ρ ∂τ ∂ρ (βρ − 2αE) ∂τ (27) β > 0, 0

According to Eq. (15), one of the corresponding linear system is Eq. (23a), with the boundary conditions

0 < βa ≤ ρ, E0( βa,τ)=0. p p 9

By the superposition principle, we get the solution of the linear problem as follows:

2 E0(ρ, τ)= Ai[Y1(̟i βa)J1(̟iρ) Xi=1 p

− J1(̟i βa)Y1(̟iρ)] sin(̟iτ), 2 p H0(ρ, τ)= Ai[Y1(̟i βa)J0(̟iρ) Xi=1 p

− J1(̟i βa)Y0(̟iρ)] cos(̟iτ). p where Jν and Yν represent the first and the second kind of the order ν of the Bessel function respectively, and A1 and A2 are amplitude factors. According to formula (21), we obtain the solution of Eqs. (27) as follows:

2 1 1 1 E = A (βρ2 − 2αρE) 2 [Y (̟ βa)J (̟ (βρ2 − 2αρE) 2 ) ρ i 1 i 1 i Xi=1 p 1 2 2 − J1(̟i βa)Y1(̟i(βρ − 2αρE) )] sin[̟i(αH + βτ)], (28) 2 p 1 2 2 H = Ai[Y1(̟i βa)J0(̟i(βρ − 2αρE) ) Xi=1 p 1 2 2 − J1(̟i βa)Y0(̟i(βρ − 2αρE) )] cos[̟i(αH + βτ)]. p Here we should properly choose A1 and A2, to make inequality 2 max(αE/ρ) < β hold. Such an explicit expression can easily be extended to deal with the problem of any amount of cylindrical electromagnetic waves propagation in nonlinear media, while applying the traditional methods is very difficult [20].

80 ω 1

60

40

Spectral amplitude 20 ω +ω ω 2 1 ω −ω 2 2 1 0 0 1 2 3 4 5 Frequency

FIG. 3. Frequency spectrum of the electric field of Eqs. (28) when two cylindrical waves propagate in a semi infinite nonlinear − 2π π 9̟1 9̟2 medium. We take ρ = 2, α = 1, β = 9, a = 1, A1 = A2 = 1, ̟1 = 9 and ̟2 = 3 (ω1 = 2π , ω2 = 2π ).

Figure 3 shows the result of the fast Fourier transform of the electric field determined by Eqs. (28) at ρ = 2 by 2π π 9̟1 9̟2 taking α = −1, β = 9, a = 1, A1 = A2 = 1, ̟1 = 9 and ̟2 = 3 (ω1 = 2π , ω2 = 2π ). We can see the sum frequency (ω1 + ω2) and the difference frequency (ω2 − ω1) which are second-order nonlinear optical processes are pronounced.

V. CONCLUSIONS

In conclusion, we have revealed the correlation between the linear and nonlinear (1+1)-dimensional Maxwell equa- tions under the orthogonal curvilinear coordinates, and shown how to acquire the exact solutions of the Maxwell equations in nondispersive media with arbitrary nonlinearity and particular inhomogeneity, which have puzzled re- searchers for a long time. The exact solutions enable us to analytically analyze some nonlinear problems difficultly tackled by traditional methods. As examples, we analyze the initial value problem, boundary value problem and sum- and difference-frequency generation. The results can be applied to studying and other nonlinear electromagnetic phenomena conveniently and effectively, as well as studying the properties of ferroelectric materials and metamaterials. 10

ACKNOWLEDGMENTS

We thank Kaijie Chen and Ding Liu for constructive advice.

[1] G.B.Whitham, Linear and nonlinear waves (Wiley, New York, 1974). [2] Y. R. Shen, Principles of Nonlinear Optics (Wiley, New York, 1984). [3] R. W. Boyd, Nonlinear Optics (Academic, New York, 2008). [4] N. Bloembergen, “Nonlinear optics and spectroscopy,” Rev. Mod. Phys. 54, 685–695 (1982). [5] E. Y. Petrov and A. V. Kudrin, “Exact axisymmetric solutions of the maxwell equations in a nonlinear nondispersive medium,” Phys. Rev. Lett. 104, 190404 (2010). [6] H. Xiong, L.-G. Si, P. Huang, and X. Yang, “Analytic description of cylindrical electromagnetic wave propagation in an inhomogeneous nonlinear and nondispersive medium,” Phys. Rev. E 82, 057602 (2010). [7] A.D. Polyanin, and A.I. Zhurov, “Exact separable solutions of delay reaction-diffusion equations and other nonlinear partial functional-differential equations,” Commun. Nonlinear Sci. 19, 409–416 (2014). [8] A.D. Polyanin, and A.I. Zhurov, “Exact solutions of linear and non-linear differential-difference heat and diffusion equations with finite relaxation time,” Int. J. Nonlin. Mech. 54, 115–126 (2013). [9] A.D. Polyanin, and A.I. Zhurov, “Integration of linear and some model non-linear equations of motion of incompressible fluids,” Int. J. Nonlin. Mech. 49, 77–83 (2013). [10] A.D. Polyanin, and A.I. Zhurov, “New generalized and functional separable solutions to non-linear delay reaction-diffusion equations,” Int. J. Nonlin. Mech. 59, 16–22 (2014). [11] R.O. Popovych, C. Sophocleous, and O.O. Vaneeva, “Exact solutions of a remarkable fin equation,” Appl. Math. Lett. 21, 209–214 (2008). [12] M. Zamboni-Rached, E. Recami, and H. Hern´andez-Figueroa, “New localized superluminal solutions to the wave equations with finite total and arbitrary frequencies,” Eur. Phys. J. D 21, 217–228 (2002). [13] M. Zamboni-Rached, E. Recami, and F. Fontana, “Superluminal localized solutions to maxwell equations propagating along a normal-sized waveguide,” Phys. Rev. E 64, 066603 (2001). [14] M. Zamboni-Rached, “Diffraction-attenuation resistant beams in absorbing media,” Opt. Express 14, 1804–1809 (2006). [15] L. A. Ambrosio and M. Zamboni-Rached, “Analytical approach of ordinary frozen waves for optical trapping and micro- manipulation,” Appl. Opt. 54, 2584–2593 (2015). [16] J. H. He, “Some Asymptotic Methods For Strongly Nonlinear Equations,” Int. J. Mod Phys B 20, 1141–1199 (2006). [17] J. H. He, “Variational iteration method - a kind of non-linear analytical technique: some examples,” Int. J. Nonlin. Mech. 34, 699 – 708 (1999). [18] J. H. He and X. H. Wu, “Exp-function method for nonlinear wave equations,” Chaos Fract. 30, 700 – 708 (2006). [19] A. V. Kudrin and E. Y. Petrov, “Cylindrical electromagnetic waves in a nonlinear nondispersive medium: Exact solutions of the maxwell equations,” JETP 110, 537–548 (2010). [20] H. Xiong, L.-G. Si, X. Yang, and Y. Wu, “Analytic descriptions of cylindrical electromagnetic waves in a nonlinear medium,” Sci. Rep. 5, 11071– (2015). [21] V. A. Es’kin, A. V. Kudrin, and E. Y. Petrov, “Exact solutions for the source-excited cylindrical electromagnetic waves in a nonlinear nondispersive medium,” Phys. Rev. E 83, 067602 (2011). [22] H. Xiong, L.-G. Si, C. Ding, X. Yang, and Y. Wu, “Classical theory of cylindrical nonlinear optics: Sum- and difference- frequency generation,” Phys. Rev. A 84, 043841 (2011). [23] H. Xiong, L.-G. Si, J. F. Guo, X.-Y. L¨u, and X. Yang, “Classical theory of cylindrical nonlinear optics: Second-,” Phys. Rev. A 83, 063845 (2011). [24] E. Y. Petrov and A. V. Kudrin, “Electromagnetic oscillations in a driven nonlinear resonator: A description of complex nonlinear dynamics,” Phys. Rev. E 85, 055202 (2012). [25] H. Xiong, L.-G. Si, C. Ding, X.-Y. L¨u, X. Yang, and Y. Wu, “Solutions of the cylindrical nonlinear maxwell equations,” Phys. Rev. E 85, 016602 (2012). [26] H. Xiong, L.-G. Si, C. Ding, X. Yang, and Y. Wu, “Second-harmonic generation of cylindrical electromagnetic waves propagating in an inhomogeneous and nonlinear medium,” Phys. Rev. E 85, 016606 (2012). [27] S.-Y. Chen, T. Li, J.-B. Xie, H. Xie, P. Zhou, Y.-F. Tian, H. Xiong, and L.-G. Si, “Initial value problems of cylindrical electromagnetic waves propagation in a nonlinear nondispersive medium,” Phys. Rev. E 88, 035202 (2013). [28] M. Ranjbar and A. Bahari, “Investigation of third order nonlinearity in propagation of cylindrical waves in homogeneous nonlinear media,” Opt. Commun. 375, 19 – 22 (2016). [29] L. Dolin, “To the possibility of comparison of three-dimensional electromagnetic systems with nonuniform anisotropic filling,” Izv. Vyssh. Uchebn. Zaved. Radiofizika 4, 964–967 (1961). [30] C. Pallikaros, and C. Sophocleous, “On point transformations of generalized nonlinear diffusion equations,” J. Phys. A: Math. Gen. 28, 6459 (1995). [31] J. G. Kingston, and C. Sophocleous, “On form-preserving point transformations of partial differential equations,” J. Phys. A: Math. Gen. 31, 1597 (1998). [32] C. Sophocleous, “Hodograph-type transformations,” Nonlinear Anal-Theor 55, 441–466 (2003). 11

[33] M. Lapine, I. V. Shadrivov, and Y. S. Kivshar, “Colloquium : Nonlinear metamaterials,” Rev. Mod. Phys. 86, 1093–1123 (2014). [34] T. J. Cui, D. Smith, and R. Liu, Metamaterials: Theory, Design, and Applications (Springer, New York, 2009). [35] J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85, 3966–3969 (2000). [36] J. E. Sipe and R. W. Boyd, “Nonlinear susceptibility of composite optical materials in the maxwell garnett model,” Phys. Rev. A 46, 1614–1629 (1992). [37] R. W. Boyd and D. J. Gauthier, “Controlling the Velocity of Light Pulses,” Science 326, 1074–1077 (2009). [38] D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electro- magnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [39] J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780–1782 (2006). [40] P. Aleahmad, H. M. Cessa, I. Kaminer, M. Segev, and D. N. Christodoulides, “Dynamics of accelerating bessel solutions of maxwell’s equations,” J. Opt. Soc. Am. A 33, 2047–2052 (2016). [41] H. Chen, C. Yang, C. Fu, L. Zhao, and Z. Gao, “The size effect of Ba0.6Sr0.4TiO3 thin films on the ferroelectric properties,” Appl. Surf. Sci. 252, 4171 – 4177 (2006). [42] P. Zubko, J. C. Wojde l, M. Hadjimichael, S. Fernandez-Pena, A. Sen´e, I. Luk’yanchuk, J.-M. Triscone, and J. ´I˝niguez, “Negative capacitance in multidomain ferroelectric superlattices,” Nature 534, 524–528 (2016). [43] D. D. Fong, G. B. Stephenson, S. K. Streiffer, J. A. Eastman, O. Auciello, P. H. Fuoss, and C. Thompson, “Ferroelectricity in ultrathin perovskite films,” Science 304, 1650–1653 (2004). [44] M. Dawber, K. M. Rabe, and J. F. Scott, “Physics of thin-film ferroelectric oxides,” Rev. Mod. Phys. 77, 1083–1130 (2005). [45] G. Rosen, “Electromagnetic shocks and the self-annihilation of intense linearly polarized radiation in an ideal dielectric material,” Phys. Rev. 139, A539–A543 (1965). [46] F. DeMartini, C. H. Townes, T. K. Gustafson, and P. L. Kelley, “Self-steepening of light pulses,” Phys. Rev. 164, 312–323 (1967).