BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeremy L. Marzuola

Jeremy L. Marzuola appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.SPmath-phmath.MPmath.NAnlin.PSNumerical Analysiscond-mat.mtrl-scimath.CAmath.DGmath.DSnlin.SIphysics.flu-dyn

Catalog footprint

What is connected

41works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

41 published item(s)

preprint2022arXiv

Eigenvalue spacing for 1D singular Schrödinger operators

The aim of this paper is to provide uniform estimates for the eigenvalue spacings of one-dimensional semiclassical Schrödinger operators with singular potentials on the half-line. We introduce a new development of semiclassical measures related to families of Schrödinger operators that provides a means of establishing uniform non-concentration estimates within that class of operators. This dramatically simplifies analysis that would typically require detailed WKB expansions near the turning point, near the singular point and several gluing type results to connect various regions in the domain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nodal Set Openings on Perturbed Rectangular Domains

We study the effects of perturbing the boundary of a rectangle on the nodal sets of eigenfunctions of the Laplacian. Namely, for a rectangle of a given aspect ratio $N$, we identify the first Dirichlet mode to feature a crossing in its nodal set and perturb one of the sides of the rectangle by a close to flat, smooth curve. Such perturbations will often "open" the crossing in the nodal set, splitting it into two curves, and we study the separation between these curves and their regularity. The main technique used is an approximate separation of variables that allows us to restrict study to the first two Fourier modes in an eigenfunction expansion. We show how the nature of the boundary perturbation provides conditions on the orientation of the opening and estimates on its size. In particular, several features of the perturbed nodal set are asymptotically independent of the aspect ratio, which contrasts with prior works. Numerical results supporting our findings are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Algorithms for Water Waves with Background Flow over Obstacles and Topography

We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pointwise dispersive estimates for Schrödinger operators on product cones

We investigate the dispersive properties of solutions to the Schrödinger equation with a weakly decaying radial potential on cones. If the potential has sufficient polynomial decay at infinity, then we show that the Schrödinger flow on each eigenspace of the link manifold satisfies a weighted $L^1\to L^\infty$ dispersive estimate. In odd dimensions, the decay rate we compute is consistent with that of the Schrödinger equation in a Euclidean space of the same dimension, but the spatial weights reflect the more complicated regularity issues in frequency that we face in the form of the spectral measure. In even dimensions, we prove a similar estimate, but with a loss of $t^{1/2}$ compared to the sharp Euclidean estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The radiation field on product cones

We consider the wave equation on a product cone and find a joint asymptotic expansion for solutions near null and future infinities. The rates of decay seen in the expansion at future infinity are the resonances of a hyperbolic cone and were computed by the authors in a previous paper. The expansion treats an asymptotic regime not considered in the influential work of Cheeger and Taylor. The main result adapts the blueprint laid out in the asymptotically Minkowski setting with key new elements including propagation estimates near the conic singularities. The proof of the propagation estimates extends prior work of Melrose-Vasy-Wunsch and Gannot-Wunsch.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Symmetry Breaking in Density Functional Theory due to Dirac Exchange for a Hydrogen Molecule

We study symmetry breaking in the mean field solutions to the 2 electron hydrogen molecule within Kohn Sham (KS) local spin density function theory with Dirac exchange (the XLDA model). This simplified model shows behavior related to that of the (KS) spin density functional theory (SDFT) predictions in condensed and molecular systems. The Kohn Sham solutions to the constrained SDFT variation problem undergo spontaneous symmetry breaking as the relative strength of the non-convex exchange term increases. This results in the change of the molecular ground state from a paramagnetic state to an antiferromagnetic ground states and a stationary symmetric delocalized 1st excited state. We further characterize the limiting behavior of the minimizer when the strength of the exchange term goes to infinity. This leads to further bifurcations and highly localized states with varying character. The stability of the various solution classes is demonstrated by Hessian analysis. Finite element numerical results provide support for the formal conjectures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Proximal-Gradient Algorithm for Crystal Surface Evolution

As a counterpoint to recent numerical methods for crystal surface evolution, which agree well with microscopic dynamics but suffer from significant stiffness that prevents simulation on fine spatial grids, we develop a new numerical method based on the macroscopic partial differential equation, leveraging its formal structure as the gradient flow of the total variation energy, with respect to a weighted $H^{-1}$ norm. This gradient flow structure relates to several metric space gradient flows of recent interest, including 2-Wasserstein flows and their generalizations to nonlinear mobilities. We develop a novel semi-implicit time discretization of the gradient flow, inspired by the classical minimizing movements scheme (known as the JKO scheme in the 2-Wasserstein case). We then use a primal dual hybrid gradient (PDHG) method to compute each element of the semi-implicit scheme. In one dimension, we prove convergence of the PDHG method to the semi-implicit scheme, under general integrability assumptions on the mobility and its reciprocal. Finally, by taking finite difference approximations of our PDHG method, we arrive at a fully discrete numerical algorithm, with iterations that converge at a rate independent of the spatial discretization: in particular, the convergence properties do not deteriorate as we refine our spatial grid. We close with several numerical examples illustrating the properties of our method, including facet formation at local maxima, pinning at local minima, and convergence as the spatial and temporal discretizations are refined.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Defect resonances of truncated crystal structures

Defects in the atomic structure of crystalline materials may spawn electronic bound states, known as \emph{defect states}, which decay rapidly away from the defect. Simplified models of defect states typically assume the defect is surrounded on all sides by an infinite perfectly crystalline material. In reality the surrounding structure must be finite, and in certain contexts the structure can be small enough that edge effects are significant. In this work we investigate these edge effects and prove the following result. Suppose that a one-dimensional infinite crystalline material hosting a positive energy defect state is truncated a distance $M$ from the defect. Then, for sufficiently large $M$, there exists a resonance \emph{exponentially close} (in $M$) to the bound state eigenvalue. It follows that the truncated structure hosts a metastable state with an exponentially long lifetime. Our methods allow both the resonance frequency and associated resonant state to be computed to all orders in $e^{-M}$. We expect this result to be of particular interest in the context of photonic crystals, where defect states are used for wave-guiding and structures are relatively small. Finally, under a mild additional assumption we prove that if the defect state has negative energy then the truncated structure hosts a bound state with exponentially-close energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Edge-localized states on quantum graphs in the limit of large mass

In this work, we construct and quantify asymptotically in the limit of large mass a variety of edge-localized stationary states of the focusing nonlinear Schrödinger equation on a quantum graph. The method is applicable to general bounded and unbounded graphs. The solutions are constructed by matching a localized large amplitude elliptic function on a single edge with an exponentially smaller remainder on the rest of the graph. This is done by studying the intersections of Dirichlet-to-Neumann manifolds (nonlinear analogues of Dirichlet-to-Neumann maps) corresponding to the two parts of the graph. For the quantum graph with a given set of pendant, looping, and internal edges, we find the edge on which the state of smallest energy at fixed mass is localized. Numerical studies of several examples are used to illustrate the analytical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Limiting Eigenfunctions of Sturm-Liouville operators Subject to a Spectral Flow

We examine the spectrum of a family of Sturm--Liouville operators with regularly spaced delta function potentials parametrized by increasing strength. The limiting behavior of the eigenvalues under this spectral flow was described in a previor result of the last two authors with Berkolaiko, where it was used to study the nodal deficiency of Laplacian eigenfunctions. Here we consider the eigenfunctions of these operators. In particular, we give explicit formulas for the limiting eigenfunctions, and also characterize the eigenfunctions and eigenvalues for all values for the spectral flow parameter (not just in the limit). We also develop spectrally accurate numerical tools for comparison and visualization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Local well-posedness for a quasilinear Schrödinger equation with degenerate dispersion

We consider a quasilinear Schrödinger equation on $\mathbb R$ for which the dispersive effects degenerate when the solution vanishes. We first prove local well-posedness for sufficiently smooth, spatially localized, degenerate initial data. As a corollary in the focusing case we obtain a short time stability result for the energy-minimizing compact breather.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the behavior of $1$-Laplacian Ratio Cuts on nearly rectangular domains

Given a connected set $Ω_0 \subset \mathbb{R}^2$, define a sequence of sets $(Ω_n)_{n=0}^{\infty}$ where $Ω_{n+1}$ is the subset of $Ω_n$ where the first eigenfunction of the (properly normalized) Neumann $p-$Laplacian $ -Δ^{(p)} ϕ= λ_1 |ϕ|^{p-2} ϕ$ is positive (or negative). For $p=1$, this is also referred to as the Ratio Cut of the domain. We conjecture that, unless $Ω_0$ is an isosceles right triangle, these sets converge to the set of rectangles with eccentricity bounded by 2 in the Gromov-Hausdorff distance as long as they have a certain distance to the boundary $\partial Ω_0$. We establish some aspects of this conjecture for $p=1$ where we prove that (1) the 1-Laplacian spectral cut of domains sufficiently close to rectangles of a given aspect ratio is a circular arc that is closer to flat than the original domain (leading eventually to quadrilaterals) and (2) quadrilaterals close to a rectangle of aspect ratio $2$ stay close to quadrilaterals and move closer to rectangles in a suitable metric. We also discuss some numerical aspects and pose many open questions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spectral band degeneracies of $\fracπ{2}-$rotationally invariant periodic Schrödinger operators

The dynamics of waves in periodic media is determined by the band structure of the underlying periodic Hamiltonian. Symmetries of the Hamiltonian can give rise to novel properties of the band structure. Here we consider a class of periodic Schrödinger operators, $H_V=-Δ+V$, where $V$ is periodic with respect to the lattice of translates $Λ=\mathbb{Z}^2$. The potential is also assumed to be real-valued, sufficiently regular and such that, with respect to some origin of coordinates, inversion symmetric (even) and invariant under $π/2$ rotation. The present results are the $\mathbb{Z}^2-$ analogue of results obtained for conical degenerate points (Dirac points) in honeycomb structures. Our proofs make use of the framework developed by Fefferman-Weinstein and Fefferman-Lee-Thorp-Weinstein.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The profile decomposition for the hyperbolic Schrödinger equation

In this note, we prove the profile decomposition for hyperbolic Schrödinger (or mixed signature) equations on $\mathbb{R}^2$ in two cases, one mass-supercritical and one mass-critical. First, as a warm up, we show that the profile decomposition works for the ${\dot H}^{\frac12}$ critical problem, which gives a simple generalization of for instance one of the results in Fanelli-Visciglia (2013). Then, we give the derivation of the profile decomposition in the mass-critical case by proving an improved Strichartz estimate. We will use a very similar approach to that laid out in the notes of Killip-Visan (2008), but we are forced to do a double Whitney decomposition to accommodate an extra scaling symmetry that arises in the problem with mixed signature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$L^p$-bounds on spectral clusters associated to polygonal domains

We look at the $L^p$ bounds on eigenfunctions for polygonal domains (or more generally Euclidean surfaces with conic singularities) by analysis of the wave operator on the flat Euclidean cone $C(\mathbb{S}^1_ρ) := \mathbb{R}_+ \times \left(\mathbb{R} \big/ 2πρ\mathbb{Z}\right)$ of radius $ρ> 0$ equipped with the metric $h(r,θ) = d r^2 + r^2 \, dθ^2$. Using explicit oscillatory integrals and relying on the fundamental solution to the wave equation in geometric regions related to flat wave propagation and diffraction by the cone point, we can prove spectral cluster estimates equivalent to those in works on smooth Riemannian manifolds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence and Stability Properties of Radial Bound States for Schrödinger-Poisson with an External Coulomb Potential in Three Space Dimensions

We consider radial solutions to the Schrödinger-Poisson system in three dimensions with an external smooth potential with Coulomb-like decay. Such a system can be viewed as a model for the interaction of dark matter with a bright matter background in the non-relativistic limit. We find that there are infinitely many critical points of the Hamiltonian, subject to fixed mass, and that these bifurcate from solutions to the associated linear problem at zero mass. As a result, each branch has a different topological character defined by the number of zeros of the radial states. We construct numerical approximations to these nonlinear states along the first several branches. The solution branches can be continued, numerically, to large mass values, where they become asymptotic, under a rescaling, to those of the Schrödinger-Poisson problem with no external potential. Our numerical computations indicate that the ground state is orbitally stable, while the excited states are linearly unstable for sufficiently large mass.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gross-Pitaevskii vortex motion with critically-scaled inhomogeneities

We study the dynamics of vortices in an inhomogeneous Gross-Pitaevskii equation $i \partial_t u = Δu + {1\over \varepsilon^2} (p_\varepsilon^2(x) - |u|^2)$. For a unique scaling regime $|p_\varepsilon(x) - 1 | = O(|\log \varepsilon|^{-1})$, it is shown that vortices can interact both with the background perturbation and with each other. Results for associated parabolic and elliptic problems are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Manifold decompositions and indices of Schrödinger operators

The Maslov index is used to compute the spectra of different boundary value problems for Schrödinger operators on compact manifolds. The main result is a spectral decomposition formula for a manifold $M$ divided into components $Ω_1$ and $Ω_2$ by a separating hypersurface $Σ$. A homotopy argument relates the spectrum of a second-order elliptic operator on $M$ to its Dirichlet and Neumann spectra on $Ω_1$ and $Ω_2$, with the difference given by the Maslov index of a path of Lagrangian subspaces. This Maslov index can be expressed in terms of the Morse indices of the Dirichlet-to-Neumann maps on $Σ$. Applications are given to doubling constructions, periodic boundary conditions and the counting of nodal domains. In particular, a new proof of Courant's nodal domain theorem is given, with an explicit formula for the nodal deficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spectrally Optimized Pointset Configurations

The search for optimal configurations of pointsets, the most notable examples being the problems of Kepler and Thompson, have an extremely rich history with diverse applications in physics, chemistry, communication theory, and scientific computing. In this paper, we introduce and study a new optimality criteria for pointset configurations. Namely, we consider a certain weighted graph associated with a pointset configuration and seek configurations which minimize certain spectral properties of the adjacency matrix or graph Laplacian defined on this graph, subject to geometric constraints on the pointset configuration. This problem can be motivated by solar cell design and swarming models, and we consider several spectral functions with interesting interpretations such as spectral radius, algebraic connectivity, effective resistance, and condition number. We prove that the regular simplex extremizes several spectral invariants on the sphere. We also consider pointset configurations on flat tori via (i) the analogous problem on lattices and (ii) through a variety of computational experiments. For many of the objectives considered (but not all), the triangular lattice is extremal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Morse index theorem for elliptic operators on bounded domains

Given a selfadjoint, elliptic operator $L$, one would like to know how the spectrum changes as the spatial domain $Ω\subset \mathbb{R}^d$ is deformed. For a family of domains $\{Ω_t\}_{t\in[a,b]}$ we prove that the Morse index of $L$ on $Ω_a$ differs from the Morse index of $L$ on $Ω_b$ by the Maslov index of a path of Lagrangian subspaces on the boundary of $Ω$. This is particularly useful when $Ω_a$ is a domain for which the Morse index is known, e.g. a region with very small volume. Then the Maslov index computes the difference of Morse indices for the "original" problem (on $Ω_b$) and the "simplified" problem (on $Ω_a$). This generalizes previous multi-dimensional Morse index theorems that were only available on star-shaped domains or for Dirichlet boundary conditions. We also discuss how one can compute the Maslov index using crossing forms, and present some applications to the spectral theory of Dirichlet and Neumann boundary value problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Euler Equation on a Rotating Surface

We study 2D Euler equations on a rotating surface, subject to the effect of the Coriolis force, with an emphasis on surfaces of revolution. We bring in conservation laws that yield long time estimates on solutions to the Euler equation, and examine ways in which the solutions behave like zonal fields, building on work of B.~Cheng and A.~Mahalov, examining how such 2D Euler equations can account for the observed band structure of rapidly rotating planets. Specific results include both an analysis of time averages of solutions and a study of stability of stationary zonal fields. The latter study includes both analytical and numerical work.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Structure of Global Attractors for Dissipative Zakharov Systems with Forcing on the Torus

The Zakharov system was originally proposed to study the propagation of Langmuir waves in an ionized plasma. In this paper, motivated by earlier work of the first and third authors, we numerically and analytically investigate the dynamics of the dissipative Zakharov system on the torus in 1 dimension. We find an interesting family of stable periodic orbits and fixed points, and explore bifurcations of those points as we take weaker and weaker dissipation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-trapping and Josephson tunneling solutions to the nonlinear Schrödinger / Gross-Pitaevskii Equation

We study the long-time behavior of solutions to the nonlinear Schrödinger / Gross-Pitaevskii equation (NLS/GP) with a symmetric double-well potential, continuing work of the 2nd and 3rd authors. NLS/GP governs nearly-monochromatic guided optical beams in weakly coupled waveguides with both linear and nonlinear (Kerr) refractive indices and zero absorption. The optical power ($L^2$ norm) is conserved with propagation distance. At low optical power, the beam energy executes beating oscillations between the two waveguides. There is an optical power threshold above which the set of guided mode solutions splits into two families of solutions. One type of solution corresponds to an optical beam which is concentrated in either waveguide, but not both. Solutions in the second family undergo tunneling oscillations between the two waveguides. NLS/GP can also model the behavior of Bose-Einstein condensates. A finite dimensional reduction (system of ODEs) well-approximates the PDE dynamics on long time scales. In particular, we derive this reduction, find a class of exact solutions and prove the very long-time shadowing of these solutions by applying the approach of the 2nd and 3rd authors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strang Splitting Methods for a quasilinear Schrödinger equation - Convergence, Instability and Dynamics

We study the Strang splitting scheme for quasilinear Schrödinger equations. We establish the convergence of the scheme for solutions with small initial data. We analyze the linear instability of the numerical scheme, which explains the numerical blow-up of large data solutions and connects to analytical breakdown of regularity of solutions to quasilinear Schrödinger equations. Numerical tests are performed for a modified version of the superfluid thin film equation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics near a minimal-mass soliton for a Korteweg-de Vries equation

We study soliton solutions to a generalized Korteweg - de Vries (KdV) equation with a saturated nonlinearity, following the line of inquiry of the authors for the nonlinear Schrödinger equation (NLS). KdV with such a nonlinearity is known to possess a minimal-mass soliton. We consider a small perturbation of a minimal-mass soliton and identify a system of ODEs, which models the behavior of the perturbation for short times. This connects nicely to a work of Comech, Cuccagna & Pelinovsky (2007). These ODEs form a simple dynamical system with a single unstable hyperbolic fixed point with two possible dynamical outcomes. A particular feature of the dynamics are that they are non-oscillatory. This distinguishes the KdV problem from the analogous NLS one.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Higher Dimensional Vortex Standing Waves for Nonlinear Schr{ö}dinger Equations

We study standing wave solutions to nonlinear Schr{ö}dinger equations, on a manifold with a rotational symmetry, which transform in a natural fashion under the group of rotations. We call these vortex solutions. They are higher dimensional versions of vortex standing waves that have been studied on the Euclidean plane. We focus on two types of vortex solutions, which we call spherical vortices and axial vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strichartz estimates on exterior polygonal domains

Using a new local smoothing estimate of the first and third authors, we prove local-in-time Strichartz and smoothing estimates without a loss exterior to a large class of polygonal obstacles with arbitrary boundary conditions and global-in-time Strichartz estimates without a loss exterior to a large class of polygonal obstacles with Dirichlet boundary conditions. In addition, we prove a global-in-time local smoothing estimate in exterior wedge domains with Dirichlet boundary conditions and discuss some nonlinear applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The relaxation of a family of broken bond crystal surface models

We study the continuum limit of a family of kinetic Monte Carlo models of crystal surface relaxation that includes both the solid-on-solid and discrete Gaussian models. With computational experiments and theoretical arguments we are able to derive several partial differential equation limits identified (or nearly identified) in previous studies and to clarify the correct choice of surface tension appearing in the PDE and the correct scaling regime giving rise to each PDE. We also provide preliminary computational investigations of a number of interesting qualitative features of the large scale behavior of the models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Behavior of a Model Dynamical System with Applications to Weak Turbulence

We experimentally explore solutions to a model Hamiltonian dynamical system derived in Colliander et al., 2012, to study frequency cascades in the cubic defocusing nonlinear Schrödinger equation on the torus. Our results include a statistical analysis of the evolution of data with localized amplitudes and random phases, which supports the conjecture that energy cascades are a generic phenomenon. We also identify stationary solutions, periodic solutions in an associated problem and find experimental evidence of hyperbolic behavior. Many of our results rely upon reframing the dynamical system using a hydrodynamic formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the regularity of the flow map for the gravity-capillary equations

We prove via explicitly constructed initial data that solutions to the gravity-capillary wave system in $\mathbb{R}^3$ representing a 2d air-water interface immediately fails to be $C^3$ with respect to the initial data if the initial data $(h_0, ψ_0) \in H^{s+\frac12} \otimes H^{s}$ for $s<3$. Similar results hold in $\mathbb{R}^2$ domains with a 1d interface. Furthermore, we discuss the illposedness threshold for the pure gravity water wave system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasilinear Schrödinger equations I: Small data and quadratic interactions

In this article we prove local well-posedness in low-regularity Sobolev spaces for general quasilinear Schrödinger equations. These results represent improvements of the pioneering works by Kenig-Ponce-Vega and Kenig-Ponce-Rolvung-Vega, where viscosity methods were used to prove existence of solutions in very high regularity spaces. Our arguments here are purely dispersive. The function spaces in which we show existence are constructed in ways motivated by the results of Mizohata, Ichinose, Doi, and others, including the authors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasilinear Schrödinger equations II: Small data and cubic nonlinearities

In part I of this project we examined low regularity local well-posedness for generic quasilinear Schrödinger equations with small data. This improved, in the small data regime, the preceding results of Kenig, Ponce, and Vega as well as Kenig, Ponce, Rolvung, and Vega. In the setting of quadratic interactions, the (translation invariant) function spaces which were utilized incorporated an $l^1$ summability over cubes in order to account for Mizohata's integrability condition, which is a necessary condition for the $L^2$ well-posedness for the linearized equation. For cubic interactions, this integrability condition meshes better with the inherent $L^2$ nature of the Schrödinger equation, and such summability is not required. Thus we are able to prove small data well-posedness in $H^s$ spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonconcentration in partially rectangular billiards

In specific types of partially rectangular billiards we estimate the mass of an eigenfunction of energy $E$ in the region outside the rectangular set in the high-energy limit. We use the adiabatic ansatz to compare the Dirichlet energy form with a second quadratic form for which separation of variables applies. This allows us to use sharp one-dimensional control estimates and to derive the bound assuming that $E$ is not resonating with the Dirichlet spectrum of the rectangular part.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On instability of some approximate periodic solutions for the full nonlinear Schrödinger equation

Using the Fermi Golden Rule analysis developed in several results by the first author, we prove asymptotic stability of asymmetric nonlinear bound states bifurcating from linear bound states for a quintic nonlinear Schrödinger operator with symmetric potential. This goes in the direction of proving that the approximate periodic solutions of the NLS shadowed in a recent work of the second author with Michael Weinstein do not persist for the full NLS.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strichartz estimates for the wave equation on flat cones

We consider the solution operator for the wave equation on the flat Euclidean cone over the circle of radius $ρ> 0$, the manifold $\mathbb{R}_+ \times \mathbb{R} / 2 πρ\mathbb{Z}$ equipped with the metric $\g(r,θ) = dr^2 + r^2 dθ^2$. Using explicit representations of the solution operator in regions related to flat wave propagation and diffraction by the cone point, we prove dispersive estimates and hence scale invariant Strichartz estimates for the wave equation on flat cones. We then show that this yields corresponding inequalities on wedge domains, polygons, and Euclidean surfaces with conic singularities. This in turn yields well-posedness results for the nonlinear wave equation on such manifolds. Morawetz estimates on the cone are also treated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Small data scattering and soliton stability in $\dot{H}^{-\frac16}$ for the quartic KdV Equation

In this note we prove scattering for perturbations of solitons in the scaling space appropriate for the quartic nonlinearity, namely $\dot{H}^{-\frac16}$. The article relies strongly on refined estimates for a KdV equation linearized at the soliton. In contrast to the work of Tao (2006), we are able to work purely in the scaling space without additional regularity assumptions, allowing us to prove some results on the existence of inverse wave operators.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral Analysis for Matrix Hamiltonian Operators

In this work, we study the spectral properties of matrix Hamiltonians generated by linearizing the nonlinear Schrödinger equation about soliton solutions. By a numerically assisted proof, we show that there are no embedded eigenvalues for the three dimensional cubic equation. Though we focus on a proof of the 3d cubic problem, this work presents a new algorithm for verifying certain spectral properties needed to study soliton stability. Source code for verification of our comptuations, and for further experimentation, are available at http://www.math.toronto.edu/simpson/files/spec_prop_code.tgz.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Strichartz estimates for the Schrödinger equation on polygonal domains

We prove Strichartz estimates with a loss of derivatives for the Schrödinger equation on polygonal domains with either Dirichlet or Neumann homogeneous boundary conditions. Using a standard doubling procedure, estimates the on polygon follow from those on Euclidean surfaces with conical singularities. We develop a Littlewood-Paley squarefunction estimate with respect to the spectrum of the Laplacian on these spaces. This allows us to reduce matters to proving estimates at each frequency scale. The problem can be localized in space provided the time intervals are sufficiently small. Strichartz estimates then follow from a result of the second author regarding the Schrödinger equation on the Euclidean cone.

0 citations0 reviewsNo code linkedOpen access