• The stationary Gross-Pitaevskii equation in one dimension is considered with a complex periodic potential satisfying the conditions of the PT (parity-time reversal) symmetry. Under rather general assumptions on the potentials we prove bifurcations of PT-symmetric nonlinear bound states from the end points of a real interval in the spectrum of the non-selfadjoint linear Schrodinger operator with a complex PT-symmetric periodic potential. The nonlinear bound states are approximated by the effective amplitude equation, which bears the form of the cubic nonlinear Schrodinger equation. In addition we provide sufficient conditions for the appearance of complex spectral bands when the complex $\PT$-symmetric potential has an asymptotically small imaginary part.
  • The dynamics of single carrier wavepackets in nonlinear wave problems over periodic structures can be often formally approximated by the constant coefficient nonlinear Schr\"odinger equation (NLS) as an effective model for the wavepacket envelope. We provide a detailed proof of this approximation result for the Gross-Pitaevskii equation (GP) and a semilinear wave equation, both with periodic coefficients in $\mathbb{N}\ni d$ spatial dimensions and with cubic nonlinearities. The proof is carried out in Bloch expansion variables with estimates in an $L^1$-type norm, which translates to an estimate of the supremum norm of the error. The regularity required from the periodic coefficients in order to ensure a small residual and a small error is discussed. We also present a numerical example in two spatial dimensions confirming the approximation result and presenting an approximate traveling solitary wave in the GP with periodic coefficients.
  • We present a new numerical scheme to solve the Helmholtz equation in a wave-guide. We consider a medium that is bounded in the $x_2$-direction, unbounded in the $x_1$-direction and $\varepsilon$-periodic for large $|x_1|$, allowing different media on the left and on the right. We suggest a new numerical method that is based on a truncation of the domain and the use of Bloch wave ansatz functions in radiation boxes. We prove the existence and a stability estimate for the infinite dimensional version of the proposed problem. The scheme is tested on several interfaces of homogeneous and periodic media and it is used to investigate the effect of negative refraction at the interface of a photonic crystal with a positive effective refractive index.
  • We consider wavepackets composed of two modulated carrier Bloch waves with opposite group velocities in the one dimensional periodic Nonlinear Schroedinger/Gross-Pitaevskii equation. These can be approximated by first order coupled mode equations (CMEs) for the two slowly varying envelopes. Under a suitably selected periodic perturbation of the periodic structure the CMEs possess a spectral gap of the corresponding spatial operator and allow families of exponentially localized solitary waves parametrized by velocity. This leads to a family of approximate solitary waves in the periodic nonlinear Schroedinger equation. Besides a formal derivation of the CMEs a rigorous justification of the approximation and an error estimate in the supremum norm are provided. Several numerical tests corroborate the analysis.
  • We rigorously analyze the bifurcation of stationary so called nonlinear Bloch waves (NLBs) from the spectrum in the Gross-Pitaevskii (GP) equation with a periodic potential, in arbitrary space dimensions. These are solutions which can be expressed as finite sums of quasi-periodic functions, and which in a formal asymptotic expansion are obtained from solutions of the so called algebraic coupled mode equations. Here we justify this expansion by proving the existence of NLBs and estimating the error of the formal asymptotics. The analysis is illustrated by numerical bifurcation diagrams, mostly in 2D. In addition, we illustrate some relations of NLBs to other classes of solutions of the GP equation, in particular to so called out--of--gap solitons and truncated NLBs, and present some numerical experiments concerning the stability of these solutions.
  • Many physical systems can be described by nonlinear eigenvalues and bifurcation problems with a linear part that is non-selfadjoint e.g. due to the presence of loss and gain. The balance of these effects is reflected in an antilinear symmetry, like e.g. the PT-symmetry, of the problem. Under this condition we show that the nonlinear eigenvalues bifurcating from real linear eigenvalues remain real and the corresponding nonlinear eigenfunctions remain symmetric. The abstract results are applied in a number of physical models of Bose-Einstein condensation, nonlinear optics and superconductivity, and further numerical analysis is performed.
  • We consider the nonlinear curl-curl problem $\nabla\times\nabla\times U + V(x) U= \Gamma(x)|U|^{p-1}U$ in $\mathbb{R}^3$ related to the nonlinear Maxwell equations for monochromatic fields. We search for solutions as minimizers (ground states) of the corresponding energy functional defined on subspaces (defocusing case) or natural constraints (focusing case) of $H(\mathrm{curl};\mathbb{R}^3)$. Under a cylindrical symmetry assumption on the functions $V$ and $\Gamma$ the variational problem can be posed in a symmetric subspace of $H(\mathrm{curl};\mathbb{R}^3)$. For a strongly defocusing case $\mathrm{esssup}\, \Gamma <0$ with large negative values of $\Gamma$ at infinity we obtain ground states by the direct minimization method. For the focusing case $\mathrm{essinf}\, \Gamma >0$ the concentration compactness principle produces ground states under the assumption that zero lies outside the spectrum of the linear operator $\nabla \times \nabla \times +V(x)$. Examples of cylindrically symmetric functions $V$ are provided for which this holds.
  • pdepath 2.0 is an upgrade of the continuation/bifurcation package pde2path for elliptic systems of PDEs over bounded 2D domains, based on the Matlab pdetoolbox. The new features include a more efficient use of FEM, easier switching between different single parameter continuations, genuine multi-parameter continuation (e.g., fold continuation), more efficient implementation of nonlinear boundary conditions, cylinder and torus geometries (i.e., periodic boundary conditions), and a general interface for adding auxiliary equations like mass conservation or phase equations for continuation of traveling waves. The package (library, demos, manuals) can be downloaded at www.staff.uni-oldenburg.de/hannes.uecker/pde2path
  • We analyze a homogenization limit for the linear wave equation of second order. The spatial operator is assumed to be of divergence form with an oscillatory coefficient matrix $a^\varepsilon$ that is periodic with characteristic length scale $\varepsilon$; no spatial symmetry properties are imposed. Classical homogenization theory allows to describe solutions $u^\varepsilon$ well by a non-dispersive wave equation on fixed time intervals $(0,T)$. Instead, when larger time intervals are considered, dispersive effects are observed. In this contribution we present a well-posed weakly dispersive equation with homogeneous coefficients such that its solutions $w^\varepsilon$ describe $u^\varepsilon$ well on time intervals $(0,T\varepsilon^{-2})$. More precisely, we provide a norm and uniform error estimates of the form $\| u^\varepsilon(t) - w^\varepsilon(t) \| \le C\varepsilon$ for $t\in (0,T\varepsilon^{-2})$. They are accompanied by computable formulas for all coefficients in the effective models. We additionally provide an $\varepsilon$-independent equation of third order that describes dispersion along rays and we present numerical examples.
  • The paper studies asymptotics of moving gap solitons in nonlinear periodic structures of finite contrast ("deep grating") within the one dimensional periodic nonlinear Schr\"odinger equation (PNLS). Periodic structures described by a finite band potential feature transversal crossings of band functions in the linear band structure and a periodic perturbation of the potential yields new small gaps. Novel gap solitons with O(1) velocity despite the deep grating are presented in these gaps. An approximation of gap solitons is given by slowly varying envelopes which satisfy a system of generalized Coupled Mode Equations (gCME) and by Bloch waves at the crossing point. The eigenspace at the crossing point is two dimensional and it is necessary to select Bloch waves belonging to the two band functions. This is achieved by an optimization algorithm. Traveling solitary wave solutions of the gCME then result in nearly solitary wave solutions of PNLS moving at an O(1) velocity across the periodic structure. A number of numerical tests are performed to confirm the asymptotics.
  • We investigate second order linear wave equations in periodic media, aiming at the derivation of effective equations in $\R^n$, $n \in \{1, 2, 3\}$. Standard homogenization theory provides, for the limit of a small periodicity length $\eps>0$, an effective second order wave equation that describes solutions on time intervals $[0,T]$. In order to approximate solutions on large time intervals $[0,T\eps^{-2}]$, one has to use a dispersive, higher order wave equation. In this work, we provide a well-posed, weakly dispersive effective equation, and an estimate for errors between the solution of the original heterogeneous problem and the solution of the dispersive wave equation. We use Bloch-wave analysis to identify a family of relevant limit models and introduce an approach to select a well-posed effective model under symmetry assumptions on the periodic structure. The analytical results are confirmed and illustrated by numerical tests.
  • We consider the problem of verifying the existence of $H^1$ ground states of the 1D nonlinear Schr\"odinger equation for an interface of two periodic structures: $$-u" +V(x)u -\lambda u = \Gamma(x) |u|^{p-1}u \ {on} \R$$ with $V(x) = V_1(x), \Gamma(x)=\Gamma_1(x)$ for $x\geq 0$ and $V(x) = V_2(x), \Gamma(x)=\Gamma_2(x)$ for $x<0$. Here $V_1,V_2,\Gamma_1,\Gamma_2$ are periodic, $\lambda <\min\sigma(-\tfrac{d^2}{dx^2}+V)$, and $p>1$. The article [T. Dohnal, M. Plum and W. Reichel, "Surface Gap Soliton Ground States for the Nonlinear Schr\"odinger Equation," \textit{Comm. Math. Phys.} \textbf{308}, 511-542 (2011)] provides in the 1D case an existence criterion in the form of an integral inequality involving the linear potentials $V_{1},V_2$ and the Bloch waves of the operators $-\tfrac{d^2}{dx^2}+V_{1,2}-\lambda$. We choose here the classes of piecewise constant and piecewise linear potentials $V_{1,2}$ and check this criterion for a set of parameter values. In the piecewise constant case the Bloch waves are calculated explicitly and in the piecewise linear case verified enclosures of the Bloch waves are computed numerically. The integrals in the criterion are evaluated via interval arithmetic so that rigorous existence statements are produced. Examples of interfaces supporting ground states are reported including such, for which ground state existence follows for all periodic $\Gamma_ {1,2}$ with $\esssup \Gamma_{1,2}>0$.
  • Out-of-plane gap solitons in 2D photonic crystals are optical beams localized in the plane of periodicity of the medium and delocalized in the orthogonal direction, in which they propagate with a nonzero velocity. We study such gap solitons as described by the Kerr nonlinear Maxwell system. Using a model of the nonlinear polarization, which does not generate higher harmonics, we obtain a closed curl-curl problem for the fundamental harmonic of the gap soliton. For gap solitons with frequencies inside spectral gaps and in an asymptotic vicinity of a gap edge we use a slowly varying envelope approximation based on the linear Bloch waves at the edge and slowly varying envelopes. We carry out a systematic derivation of the coupled mode equations (CMEs) which govern the envelopes. This derivation needs to be carried out in Bloch variables. The CMEs are a system of coupled nonlinear stationary Schr\"odinger equations with an additional cross derivative term. Examples of gap soliton approximations are numerically computed for a photonic crystal with a hexagonal periodicity cell and an annulus material structure in the cell.
  • We examine numerically vortex families near band edges of the Bloch wave spectrum in the Gross--Pitaevskii equation with a two-dimensional periodic potential and in the discrete nonlinear Schroedinger equation. We show that besides vortex families that terminate at a small distance from the band edges via fold bifurcations there exist vortex families that are continued all way to the band edges.
  • We consider the nonlinear Schr\"{o}dinger equation $(-\Delta +V(x))u = \Gamma(x) |u|^{p-1}u$, $x\in \R^n$ with $V(x) = V_1(x) \chi_{\{x_1>0\}}(x)+V_2(x) \chi_{\{x_1<0\}}(x)$ and $\Gamma(x) = \Gamma_1(x) \chi_{\{x_1>0\}}(x)+\Gamma_2(x) \chi_{\{x_1<0\}}(x)$ and with $V_1, V_2, \Gamma_1, \Gamma_2$ periodic in each coordinate direction. This problem describes the interface of two periodic media, e.g. photonic crystals. We study the existence of ground state $H^1$ solutions (surface gap soliton ground states) for $0<\min \sigma(-\Delta +V)$. Using a concentration compactness argument, we provide an abstract criterion for the existence based on ground state energies of each periodic problem (with $V\equiv V_1, \Gamma\equiv \Gamma_1$ and $V\equiv V_2, \Gamma\equiv \Gamma_2$) as well as a more practical criterion based on ground states themselves. Examples of interfaces satisfying these criteria are provided. In 1D it is shown that, surprisingly, the criteria can be reduced to conditions on the linear Bloch waves of the operators $-\tfrac{d^2}{dx^2} +V_1(x)$ and $-\tfrac{d^2}{dx^2} +V_2(x)$.
  • The nonlinear Schr\"{o}dinger equation with a linear periodic potential and a nonlinearity coefficient $\Gamma$ with a discontinuity supports stationary localized solitary waves with frequencies inside spectral gaps, so called surface gap solitons (SGSs). We compute families of 1D SGSs using the arclength continuation method for a range of values of the jump in $\Gamma$. Using asymptotics, we show that when the frequency parameter converges to the bifurcation gap edge, the size of the allowed jump in $\Gamma$ converges to 0 for SGSs centered at any $x_c\in \R$. Linear stability of SGSs is next determined via the numerical Evans function method, in which the stable and unstable manifolds corresponding to the 0 solution of the linearized spectral ODE problem need to be evolved. Zeros of the Evans function coincide with eigenvalues of the linearized operator. Far from the SGS center the manifolds are spanned by exponentially decaying/increasing Bloch functions. Evolution of the manifolds suffers from stiffness but a numerically stable formulation is possible in the exterior algebra formulation and with the use of Grassmanian preserving ODE integrators. Eigenvalues with positive real part above a small constant are then detected using the complex argument principle and a contour parallel to the imaginary axis. The location of real eigenvalues is found via a straightforward evaluation of the Evans function along the real axis and several complex eigenvalues are located using M\"{u}ller's method. The numerical Evans function method is described in detail. Our results show the existence of both unstable and stable SGSs (possibly with a weak instability), where stability is obtained even for some SGSs centered in the domain half with the less focusing nonlinearity. Direct simulations of the PDE for selected SGS examples confirm the results of Evans function computations.
  • Gap solitons near a band edge of a spatially periodic nonlinear PDE can be formally approximated by solutions of Coupled Mode Equations (CMEs). Here we study this approximation for the case of the 2D Periodic Nonlinear Schr\"{o}dinger / Gross-Pitaevskii Equation with a non-separable potential of finite contrast. We show that unlike in the case of separable potentials [T. Dohnal, D. Pelinovsky, and G. Schneider, J. Nonlin. Sci. {\bf 19}, 95--131 (2009)] the CME derivation has to be carried out in Bloch rather than physical coordinates. Using the Lyapunov-Schmidt reduction we then give a rigorous justification of the CMEs as an asymptotic model for reversible non-degenerate gap solitons and even potentials and provide $H^s$ estimates for this approximation. The results are confirmed by numerical examples including some new families of CMEs and gap solitons absent for separable potentials.
  • This paper constructs perfectly matched layers (PML) for a system of 2D Coupled Nonlinear Schr\"odinger equations with mixed derivatives which arises in the modeling of gap solitons in nonlinear periodic structures with a non-separable linear part. The PML construction is performed in Laplace Fourier space via a modal analysis and can be viewed as a complex change of variables. The mixed derivatives cause the presence of waves with opposite phase and group velocities, which has previously been shown to cause instability of layer equations in certain types of hyperbolic problems. Nevertheless, here the PML is stable if the absorption function $\sigma$ lies below a specified threshold. The PML construction and analysis are carried out for the linear part of the system. Numerical tests are then performed in both the linear and nonlinear regimes checking convergence of the error with respect to the layer width and showing that the PML performs well even in many nonlinear simulations.
  • We consider the existence of localized modes corresponding to eigenvalues of the periodic Schr\"{o}dinger operator $-\partial_x^2+ V(x)$ with an interface. The interface is modeled by a jump either in the value or the derivative of $V(x)$ and, in general, does not correspond to a localized perturbation of the perfectly periodic operator. The periodic potentials on each side of the interface can, moreover, be different. As we show, eigenvalues can only occur in spectral gaps. We pose the eigenvalue problem as a $C^1$ gluing problem for the fundamental solutions (Bloch functions) of the second order ODEs on each side of the interface. The problem is thus reduced to finding matchings of the ratio functions $R_\pm=\frac{\psi_\pm'(0)}{\psi_\pm(0)}$, where $\psi_\pm$ are those Bloch functions that decay on the respective half-lines. These ratio functions are analyzed with the help of the Pr\"{u}fer transformation. The limit values of $R_\pm$ at band edges depend on the ordering of Dirichlet and Neumann eigenvalues at gap edges. We show that the ordering can be determined in the first two gaps via variational analysis for potentials satisfying certain monotonicity conditions. Numerical computations of interface eigenvalues are presented to corroborate the analysis.
  • We investigate the statistics of flat-top solitary wave parameters in the presence of weak multiplicative dissipative disorder. We consider first propagation of solitary waves of the cubic-quintic nonlinear Schr\"odinger equation (CQNLSE) in the presence of disorder in the cubic nonlinear gain. We show by a perturbative analytic calculation and by Monte Carlo simulations that the probability density function (PDF) of the amplitude $\eta$ exhibits loglognormal divergence near the maximum possible amplitude $\eta_{m}$, a behavior that is similar to the one observed earlier for disorder in the linear gain [A. Peleg et al., Phys. Rev. E {\bf 72}, 027203 (2005)]. We relate the loglognormal divergence of the amplitude PDF to the super-exponential approach of $\eta$ to $\eta_{m}$ in the corresponding deterministic model with linear/nonlinear gain. Furthermore, for solitary waves of the derivative CQNLSE with weak disorder in the linear gain both the amplitude and the group velocity $\beta$ become random. We therefore study analytically and by Monte Carlo simulations the PDF of the parameter $p$, where $p=\eta/(1-\varepsilon_s\beta/2)$ and $\varepsilon_s$ is the self-steepening coefficient. Our analytic calculations and numerical simulations show that the PDF of $p$ is loglognormally divergent near the maximum $p$-value.
  • We demonstrate existence of waves localized at the interface of two nonlinear periodic media with different coefficients of the cubic nonlinearity via the one-dimensional Gross--Pitaevsky equation. We call these waves the surface gap solitons (SGS). In the case of smooth symmetric periodic potentials, we study analytically bifurcations of SGS's from standard gap solitons and determine numerically the maximal jump of the nonlinearity coefficient allowing for the SGS existence. We show that the maximal jump vanishes near the thresholds of bifurcations of gap solitons. In the case of continuous potentials with a jump in the first derivative at the interface, we develop a homotopy method of continuation of SGS families from the solution obtained via gluing of parts of the standard gap solitons and study existence of SGS's in the photonic band gaps. We explain the termination of the SGS families in the interior points of the band gaps from the bifurcation of linear bound states in the continuous non-smooth potentials.
  • We address a two-dimensional nonlinear elliptic problem with a finite-amplitude periodic potential. For a class of separable symmetric potentials, we study the bifurcation of the first band gap in the spectrum of the linear Schr\"{o}dinger operator and the relevant coupled-mode equations to describe this bifurcation. The coupled-mode equations are derived by the rigorous analysis based on the Fourier--Bloch decomposition and the Implicit Function Theorem in the space of bounded continuous functions vanishing at infinity. Persistence of reversible localized solutions, called gap solitons, beyond the coupled-mode equations is proved under a non-degeneracy assumption on the kernel of the linearization operator. Various branches of reversible localized solutions are classified numerically in the framework of the coupled-mode equations and convergence of the approximation error is verified. Error estimates on the time-dependent solutions of the Gross--Pitaevskii equation and the coupled-mode equations are obtained for a finite-time interval.