BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jeremy L. Marzuola

Jeremy L. Marzuola contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.SPmath-phmath.MPmath.NANumerical Analysiscond-mat.mtrl-scimath.CAmath.DGnlin.PSnlin.SIphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
12topics
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

14 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