Researcher profile

Tomáš Dohnal

Tomáš Dohnal contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2015arXiv

Bifurcation of Nonlinear Bloch Waves from the Spectrum in the Gross-Pitaevskii Equation

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.

preprint2014arXiv

Ground States of a Nonlinear Curl-Curl Problem in Cylindrically Symmetric Media

We consider the nonlinear curl-curl problem $\nabla\times\nabla\times U + V(x) U= Γ(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 $Γ$ the variational problem can be posed in a symmetric subspace of $H(\mathrm{curl};\mathbb{R}^3)$. For a strongly defocusing case $\mathrm{esssup}\, Γ<0$ with large negative values of $Γ$ at infinity we obtain ground states by the direct minimization method. For the focusing case $\mathrm{essinf}\, Γ>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.

preprint2011arXiv

Families of Surface Gap Solitons and their Stability via the Numerical Evans Function Method

The nonlinear Schrödinger equation with a linear periodic potential and a nonlinearity coefficient $Γ$ 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 $Γ$. Using asymptotics, we show that when the frequency parameter converges to the bifurcation gap edge, the size of the allowed jump in $Γ$ 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üller&#39;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.

preprint2011arXiv

Surface Gap Soliton Ground States for the Nonlinear Schrödinger Equation

We consider the nonlinear Schrödinger equation $(-Δ+V(x))u = Γ(x) |u|^{p-1}u$, $x\in \R^n$ with $V(x) = V_1(x) χ_{\{x_1>0\}}(x)+V_2(x) χ_{\{x_1<0\}}(x)$ and $Γ(x) = Γ_1(x) χ_{\{x_1>0\}}(x)+Γ_2(x) χ_{\{x_1<0\}}(x)$ and with $V_1, V_2, Γ_1, Γ_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 σ(-Δ+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, Γ\equiv Γ_1$ and $V\equiv V_2, Γ\equiv Γ_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)$.

preprint2010arXiv

Coupled Mode Equations and Gap Solitons for the 2D Gross-Pitaevskii equation with a non-separable periodic potential

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ö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.

preprint2009arXiv

Diverging probability density functions for flat-top solitary waves

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ödinger 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 $η$ exhibits loglognormal divergence near the maximum possible amplitude $η_{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 $η$ to $η_{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 $β$ become random. We therefore study analytically and by Monte Carlo simulations the PDF of the parameter $p$, where $p=η/(1-\varepsilon_sβ/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.

preprint2009arXiv

Perfectly Matched Layers for Coupled Nonlinear Schrödinger Equations with Mixed Derivatives

This paper constructs perfectly matched layers (PML) for a system of 2D Coupled Nonlinear Schrödinger 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 $σ$ 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.