
The stationary GrossPitaevskii equation in one dimension is considered with
a complex periodic potential satisfying the conditions of the PT (paritytime
reversal) symmetry. Under rather general assumptions on the potentials we prove
bifurcations of PTsymmetric nonlinear bound states from the end points of a
real interval in the spectrum of the nonselfadjoint linear Schrodinger
operator with a complex PTsymmetric 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 GrossPitaevskii 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
waveguide. 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/GrossPitaevskii 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 GrossPitaevskii (GP) equation with a
periodic potential, in arbitrary space dimensions. These are solutions which
can be expressed as finite sums of quasiperiodic 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 outofgap 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 nonselfadjoint 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 PTsymmetry, 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
BoseEinstein condensation, nonlinear optics and superconductivity, and further
numerical analysis is performed.

We consider the nonlinear curlcurl problem $\nabla\times\nabla\times U +
V(x) U= \Gamma(x)U^{p1}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
multiparameter 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.unioldenburg.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 nondispersive 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 wellposed 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 wellposed, 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
Blochwave analysis to identify a family of relevant limit models and introduce
an approach to select a wellposed 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^{p1}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},
511542 (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$.

Outofplane 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 curlcurl 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 GrossPitaevskii equation with a twodimensional 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^{p1}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 / GrossPitaevskii Equation with a nonseparable 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}, 95131
(2009)] the CME derivation has to be carried out in Bloch rather than physical
coordinates. Using the LyapunovSchmidt reduction we then give a rigorous
justification of the CMEs as an asymptotic model for reversible nondegenerate
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
nonseparable 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 halflines. 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 flattop solitary wave parameters in the
presence of weak multiplicative dissipative disorder. We consider first
propagation of solitary waves of the cubicquintic 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 superexponential 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 selfsteepening
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
onedimensional GrossPitaevsky 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
nonsmooth potentials.

We address a twodimensional nonlinear elliptic problem with a
finiteamplitude 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 coupledmode equations to
describe this bifurcation. The coupledmode equations are derived by the
rigorous analysis based on the FourierBloch 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 coupledmode equations is proved under a nondegeneracy assumption
on the kernel of the linearization operator. Various branches of reversible
localized solutions are classified numerically in the framework of the
coupledmode equations and convergence of the approximation error is verified.
Error estimates on the timedependent solutions of the GrossPitaevskii
equation and the coupledmode equations are obtained for a finitetime
interval.