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

Michael I. Weinstein

Michael I. Weinstein appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath-phmath.MPcond-mat.mes-hallnlin.PSquant-phphysics.opticscond-mat.mtrl-scimath.OCphysics.geo-phcond-mat.softmath.SPphysics.flu-dyn

Catalog footprint

What is connected

25works
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

25 published item(s)

preprint2022arXiv

Lower Bound on Quantum Tunneling for Strong Magnetic Fields

We consider a particle bound to a two-dimensional plane and a double well potential, subject to a perpendicular uniform magnetic field . The energy difference between the lowest two eigenvalues--the eigenvalue splitting--is related to the tunneling probability between the two wells. We obtain upper and lower bounds on this splitting in the regime where both the magnetic field strength and the depth of the wells are large. The main step is a lower bound on the hopping probability between the wells, a key parameter in tight binding models of solid state physics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tight-Binding Reduction and Topological Equivalence in Strong Magnetic Fields

Topological insulators (TIs) are a class of materials which are insulating in their bulk form yet, upon introduction of an a boundary or edge, e.g. by abruptly terminating the material, may exhibit spontaneous current along their boundary. This property is quantified by topological indices associated with either the bulk or the edge system. In the field of condensed matter physics, tight binding (discrete) approximate models, parametrized by hopping coefficients, have been used successfully to capture the topological behavior of TIs in many settings. However, whether such tight binding models capture the same topological features as the underlying continuum models of quantum physics has been an open question. We resolve this question in the context of the archetypal example of topological behavior in materials, the integer quantum Hall effect. We study a class of continuum Hamiltonians, $H^λ$, which govern electron motion in a two-dimensional crystal under the influence of a perpendicular magnetic field. No assumption is made on translation invariance of the crystal. We prove, in the regime where both the magnetic field strength and depth of the crystal potential are sufficiently large, $λ\gg1$, that the low-lying energy spectrum and eigenstates (and corresponding large time dynamics) of $H^λ$ are well-described by a scale-free discrete Hamiltonian, $H^{\rm TB}$; we show norm resolvent convergence. The relevant topological index is the Hall conductivity, which is expressible as a Fredholm index. We prove that for large $λ$ the topological indices of $H^λ$ and $H^{\rm TB}$ agree. This is proved separately for bulk and edge geometries. Our results justify the principle of using discrete models in the study of topological matter.

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
preprint2019arXiv

Dirac operators and domain walls

We study the eigenvalue problem for a one-dimensional Dirac operator with a spatially varying ``mass'' term. It is well-known that when the mass function has the form of a kink, or \emph{domain wall}, transitioning between strictly positive and strictly negative asymptotic mass, $\pmκ_\infty$, at $\pm\infty$, the Dirac operator has a simple eigenvalue of zero energy (geometric multiplicity equal to one) within a gap in the continuous spectrum, with corresponding \emph{zero mode}, an exponentially localized eigenfunction. We prove that when the mass function has the form of \emph{two} domain walls separated by a sufficiently large distance $2 δ$, the Dirac operator has two real simple eigenvalues of opposite sign and of order $e^{- 2 |κ_\infty| δ}$. The associated eigenfunctions are, up to $L^2$ error of order $e^{- 2 |κ_\infty| δ}$, linear combinations of shifted copies of the single domain wall zero mode. For the case of three domain walls, there are two non-zero simple eigenvalues as above and a simple eigenvalue at energy zero. Our methods are based on a Lyapunov-Schmidt reduction strategy and we outline their natural extension to the case of $n$ domain walls for which the minimal distance between domain walls is sufficiently large. The class of Dirac operators we consider controls the bifurcation of topologically protected ``edge states'' from Dirac points (linear band crossings) for classes of Schrödinger operators with domain-wall modulated periodic potentials in one and two space dimensions. The present results may be used to construct a rich class of defect modes in periodic structures modulated by multiple domain walls.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Wavepackets in inhomogeneous periodic media: propagation through a one-dimensional band crossing

We consider a model of an electron in a crystal moving under the influence of an external electric field: Schroedinger's equation in one spatial dimension with a potential which is the sum of a periodic function $V$ and a smooth function $W$. We assume that the period of $V$ is much shorter than the scale of variation of $W$ and denote the ratio of these scales by $ε$. We consider the dynamics of $\textit{semiclassical wavepacket}$ asymptotic (in the limit $ε\downarrow 0$) solutions which are spectrally localized near to a $\textit{crossing}$ of two Bloch band dispersion functions of the periodic operator $- \frac{1}{2} \partial_z^2 + V(z)$. We show that the dynamics is qualitatively different from the case where bands are well-separated: at the time the wavepacket is incident on the band crossing, a second wavepacket is `excited' which has $\textit{opposite}$ group velocity to the incident wavepacket. We then show that our result is consistent with the solution of a `Landau-Zener'-type model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Edge states in honeycomb structures

An edge state is a time-harmonic solution of a conservative wave system, e.g. Schroedinger, Maxwell, which is propagating (plane-wave-like) parallel to, and localized transverse to, a line-defect or "edge". Topologically protected edge states are edge states which are stable against spatially localized (even strong) deformations of the edge. First studied in the context of the quantum Hall effect, protected edge states have attracted huge interest due to their role in the field of topological insulators. Theoretical understanding of topological protection has mainly come from discrete (tight-binding) models and direct numerical simulation. In this paper we consider a rich family of continuum PDE models for which we rigorously study regimes where topologically protected edge states exist. Our model is a class of Schroedinger operators on $\mathbb{R}^2$ with a background 2D honeycomb potential perturbed by an "edge-potential". The edge potential is a domain-wall interpolation, transverse to a prescribed "rational" edge, between two distinct periodic structures. General conditions are given for the bifurcation of a branch of topologically protected edge states from Dirac points of the background honeycomb structure. The bifurcation is seeded by the zero mode of a 1D effective Dirac operator. A key condition is a spectral no-fold condition for the prescribed edge. We then use this result to prove the existence of topologically protected edge states along zigzag edges of certain honeycomb structures. Our results are consistent with the physics literature and appear to be the first rigorous results on the existence of topologically protected edge states for continuum 2D PDE systems describing waves in a non-trivial periodic medium. We also show that the family of Hamiltonians we study contains cases where zigzag edge states exist, but which are not topologically protected.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wavepackets in inhomogeneous periodic media: effective particle-field dynamics and Berry curvature

We consider a model of an electron in a crystal moving under the influence of an external electric field: Schrödinger's equation with a potential which is the sum of a periodic function and a general smooth function. We identify two dimensionless parameters: (re-scaled) Planck's constant and the ratio of the lattice spacing to the scale of variation of the external potential. We consider the special case where both parameters are equal and denote this parameter $ε$. In the limit $ε\downarrow 0$, we prove the existence of solutions known as semiclassical wavepackets which are asymptotic up to `Ehrenfest time' $t \sim \ln 1/ε$. To leading order, the center of mass and average quasi-momentum of these solutions evolve along trajectories generated by the classical Hamiltonian given by the sum of the Bloch band energy and the external potential. We then derive all corrections to the evolution of these observables proportional to $ε$. The corrections depend on the gauge-invariant Berry curvature of the Bloch band, and a coupling to the evolution of the wave-packet envelope which satisfies Schrödinger's equation with a time-dependent harmonic oscillator Hamiltonian. This infinite dimensional coupled `particle-field' system may be derived from an `extended' $ε$-dependent Hamiltonian. It is known that such coupling of observables (discrete particle-like degrees of freedom) to the wave-envelope (continuum field-like degrees of freedom) can have a significant impact on the overall dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Localized States and Dynamics in the Nonlinear Schroedinger / Gross-Pitaevskii Equation

This article is a review of results on the nonlinear Schroedinger / Gross-Pitaevskii equation (NLS / GP). Nonlinear bound states and aspects of their stability theory are discussed from variational and bifurcation perspectives. Nonlinear bound states, in the particular cases where the potential is a single-well, double-well potential and periodic potential, are discussed in detail. We then discuss particle-like dynamics of solitary wave solutions interacting with a potential over a large, but finite, time interval. Finally we turn to the very long time behavior of solutions. We focus on the important resonant radiation damping mechanism that drives the relaxation of the system to a nonlinear ground state and underlies the phenomena of {\it Ground State Selection} and {\it Energy Equipartition}. We also analyze linear and nonlinear "toy minimal models", which illustrate these mechanisms. Regarding the overall style of this article, we seek to emphasize the key ideas and therefore do not present all technical details, leaving that to the references.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On-Site and Off-Site Bound States of the Discrete Nonlinear Schrödinger Equation and the Peierls-Nabarro Barrier

We construct multiple families of solitary standing waves of the discrete cubically nonlinear Schrödinger equation (DNLS) in dimensions $d=1,2$ and $3$. These states are obtained via a bifurcation analysis about the continuum (NLS) limit. One family consists {\it on-site symmetric} (vertex-centered) states; these are spatially localized solitary standing waves which are symmetric about any fixed lattice site. The other spatially localized states are {\it off-site symmetric}. Depending on the spatial dimension, these may be bond-centered, cell-centered, or face-centered. Finally, we show that the energy difference among distinct states of the same frequency is exponentially small with respect to a natural parameter. This provides a rigorous bound for the so-called {\it Peierls-Nabarro} energy barrier.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topologically Protected States in One-Dimensional Systems

We study a class of periodic Schrödinger operators, which in distinguished cases can be proved to have linear band-crossings or "Dirac points". We then show that the introduction of an "edge", via adiabatic modulation of these periodic potentials by a domain wall, results in the bifurcation of spatially localized "edge states". These bound states are associated with the topologically protected zero-energy mode of an asymptotic one-dimensional Dirac operator. Our model captures many aspects of the phenomenon of topologically protected edge states for two-dimensional bulk structures such as the honeycomb structure of graphene. The states we construct can be realized as highly robust TM- electromagnetic modes for a class of photonic waveguides with a phase-defect.

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
preprint2013arXiv

Long-lived Scattering Resonances and Bragg Structures

We consider a system governed by the wave equation with index of refraction $n(x)$, taken to be variable within a bounded region $Ω\subset \mathbb R^d$, and constant in $\mathbb R^d \setminus Ω$. The solution of the time-dependent wave equation with initial data, which is localized in $Ω$, spreads and decays with advancing time. This rate of decay can be measured (for $d=1,3$, and more generally, $d$ odd) in terms of the eigenvalues of the scattering resonance problem, a non-selfadjoint eigenvalue problem governing the time-harmonic solutions of the wave (Helmholtz) equation which are outgoing at $\infty$. Specifically, the rate of energy escape from $Ω$ is governed by the complex scattering eigenfrequency, which is closest to the real axis. We study the structural design problem: Find a refractive index profile $n_*(x)$ within an admissible class which has a scattering frequency with minimal imaginary part. The admissible class is defined in terms of the compact support of $n(x)-1$ and pointwise upper and lower (material) bounds on $n(x)$ for $x \in Ω$, i.e., $0 < n_- \leq n(x) \leq n_+ < \infty$. We formulate this problem as a constrained optimization problem and prove that an optimal structure, $n_*(x)$ exists. Furthermore, $n_*(x)$ is piecewise constant and achieves the material bounds, i.e., $n_*(x) \in {n_-, n_+} $. In one dimension, we establish a connection between $n_*(x)$ and the well-known class of Bragg structures, where $n(x)$ is constant on intervals whose length is one-quarter of the effective wavelength.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A variational perspective on cloaking by anomalous localized resonance

A body of literature has developed concerning "cloaking by anomalous localized resonance". The mathematical heart of the matter involves the behavior of a divergence-form elliptic equation in the plane, $\nabla\cdot (a(x)\nabla u(x)) = f(x)$. The complex-valued coefficient has a matrix-shell-core geometry, with real part equal to 1 in the matrix and the core, and -1 in the shell; one is interested in understanding the resonant behavior of the solution as the imaginary part of $a(x)$ decreases to zero (so that ellipticity is lost). Most analytical work in this area has relied on separation of variables, and has therefore been restricted to radial geometries. We introduce a new approach based on a pair of dual variational principles, and apply it to some non-radial examples. In our examples, as in the radial setting, the spatial location of the source $f$ plays a crucial role in determining whether or not resonance occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Honeycomb Lattice Potentials and Dirac Points

We prove that the two-dimensional Schroedinger operator with a potential having the symmetry of a honeycomb structure has dispersion surfaces with conical singularities (Dirac points) at the vertices of its Brillouin zone. No assumptions are made on the size of the potential. We then prove the robustness of such conical singularities to a restrictive class of perturbations, which break the honeycomb lattice symmetry. General small perturbations of potentials with Dirac points do not have Dirac points; their dispersion surfaces are smooth. The presence of Dirac points in honeycomb structures is associated with many novel electronic and optical properties of materials such as graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The lifetime of shape oscillations of a bubble in an unbounded, inviscid and compressible fluid with surface tension

General perturbations of a spherical gas bubble in a compressible and inviscid fluid with surface tension were proved in Shapiro and Weinstein (2011), in the linearized approximation, to decay exponentially, $\sim e^{-Γt}, Γ>0$, as time advances. Formal asymptotic and numerical evidence led to the conjecture that $Γ\approx \frac{A}ε \frac{We}{ε^{2}} \exp(-B \frac{We}{ε^2})$, where $0<ε\ll1$ is the Mach number, We is the Weber number, and $A$ and $B$ are positive constants. In this paper, we prove this conjecture and calculate $A$ and $B$ to leading order in $ε$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wave packets in Honeycomb Structures and Two-Dimensional Dirac Equations

In a recent article [10], the authors proved that the non-relativistic Schrödinger operator with a generic honeycomb lattice potential has conical (Dirac) points in its dispersion surfaces. These conical points occur for quasi-momenta, which are located at the vertices of the Brillouin zone, a regular hexagon. In this paper, we study the time-evolution of wave-packets, which are spectrally concentrated near such conical points. We prove that the large, but finite, time dynamics is governed by the two-dimensional Dirac equations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Waves in Honeycomb Structures

We review recent work of the authors on the non-relativistic Schrödinger equation with a honeycomb lattice potential, $V$. In particular, we summarize results on (i) the existence of Dirac points, conical singularities in dispersion surfaces of $H_V=-Δ+V$ and (ii) the two-dimensional Dirac equations, as a large, but finite time, effective description of $e^{-iH_Vt}ψ_0$, for data $ψ_0$, which is spectrally localized at a Dirac point. We conclude with a formal derivation and discussion of the effective large time evolution for the nonlinear Schrödinger - Gross Pitaevskii equation for small amplitude initial conditions, $ψ_0$. The effective dynamics are governed by a nonlinear Dirac system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Broad Band Solitons in a Periodic and Nonlinear Maxwell System

We consider the one-dimensional Maxwell equations with low contrast periodic linear refractive index and weak Kerr nonlinearity. In this context, wave packet initial conditions with a single carrier frequency excite infinitely many resonances. On large but finite time-scales, the coupled evolution of backward and forward waves is governed by nonlocal equations of resonant nonlinear geometrical optics. For the special class of solutions which are periodic in the fast phase, these equations are equivalent to an infinite system of nonlinear coupled mode equations, the so called it extended nonlinear coupled mode equations, xNLCME. Numerical studies support the existence of long-lived spatially localized coherent structures, featuring a slowly varying envelope and a train of carrier shocks. In this paper we explore, by analytical, asymptotic and numerical methods, the existence and properties of spatially localized structures of the xNLCME system, which arises for a refractive index profile consisting of periodic array of Dirac delta functions. We consider the limit of small amplitude solutions with frequencies near a band-edge. In this case, stationary xNLCME is well-approximated by an infinite system of coupled, stationary, nonlinear Schrödinger equations, the extended nonlinear Schrödinger system, xNLS. We embed xNLS in a one-parameter family of equations, xNLS$^ε$, which interpolates between infinitely many decoupled NLS equations ($ε=0$) and xNLS ($ε=1$). Using bifurcation methods we show existence of solutions for a range of $ε\in(-ε_0,ε_0)$ and, by a numerical continuation method, establish the continuation of certain branches all the way to $ε=1$. Finally, we perform time-dependent simulations of truncated xNLCME and find the small-amplitude solitons to be robust to both numerical errors and the NLS approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Band-edge solitons, Nonlinear Schrodinger / Gross-Pitaevskii Equations and Effective Media

We consider a class of nonlinear Schrodinger / Gross-Pitaevskii (NLS/GP) equations with periodic potentials, having an even symmetry. We construct "solitons", centered about any point of symmetry of the potential. For focusing (attractive) nonlinearities, these solutions bifurcate from the zero state at the lowest band edge frequency, into the semi-infinite spectral gap. Our results extend to bifurcations into finite spectral gaps, for focusing or defocusing (repulsive) nonlinearities under more restrictive hypotheses. Soliton nonlinear bound states with frequencies near a band edge are well-approximated by a slowly decaying solution of a homogenized NLS/GP equation, with constant homogenized effective mass tensor and effective nonlinear coupling coefficient, modulated by a Bloch state. For the critical NLS equation with a periodic potential, e.g. the cubic two dimensional NLS/GP with a periodic potential, our results imply that the limiting soliton power, as the spectral band edge frequency is approached, is equal to a constant ζ_* times the minimal mass soliton of the translation invariant critical NLS equation. ζ_* is expressible in terms of the band edge Bloch eigenfunction and the determinant of the effective mass tensor; and 0<ζ_*<1$ for any non-constant potential. The results are confirmed by numerical computation of bound states with frequencies near the spectral band edge. Finally, these results have implications for the control of nonlinear waves using periodic structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coherent Structures and Carrier Shocks in the Nonlinear Periodic Maxwell Equations

We consider the one-dimensional propagation of electromagnetic waves in a weakly nonlinear and low-contrast spatially inhomogeneous medium with no energy dissipation. We focus on the case of a periodic medium, in which dispersion enters only through the (Floquet-Bloch) spectral band dispersion associated with the periodic structure; chromatic dispersion (time-nonlocality of the polarization) is neglected. Numerical simulations show that for initial conditions of wave-packet type (a plane wave of fixed carrier frequency multiplied by a slow varying, spatially localized function) very long-lived spatially localized coherent soliton-like structures emerge, whose character is that of a slowly varying envelope of a train of shocks. We call this structure an envelope carrier-shock train. The structure of the solution violates the oft-assumed nearly monochromatic wave packet structure, whose envelope is governed by the nonlinear coupled mode equations (NLCME). The inconsistency and inaccuracy of NLCME lies in the neglect of all (infinitely many) resonances except for the principle resonance induced by the initial carrier frequency. We derive, via a nonlinear geometrical optics expansion, a system of nonlocal integro-differential equations governing the coupled evolution of backward and forward propagating waves. These equations incorporate effects of all resonances. In a periodic medium, these equations may be expressed as a system of infinitely many coupled mode equations, which we call the extended nonlinear coupled mode system (xNLCME). Truncating xNLCME to include only the principle resonances leads to the classical NLCME. Numerical simulations of xNLCME demonstrate that it captures both large scale features, related to third harmonic generation, and fine scale carrier shocks features of the nonlinear periodic Maxwell equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Emergence of Periodic Structure from Maximizing the Lifetime of a Bound State Coupled to Radiation

Consider a system governed by the time-dependent Schrödinger equation in its ground state. When subjected to weak (size $ε$) parametric forcing by an "ionizing field" (time-varying), the state decays with advancing time due to coupling of the bound state to radiation modes. The decay-rate of this metastable state is governed by {\it Fermi's Golden Rule}, $Γ[V]$, which depends on the potential $V$ and the details of the forcing. We pose the potential design problem: find $V_{opt}$ which minimizes $Γ[V]$ (maximizes the lifetime of the state) over an admissible class of potentials with fixed spatial support. We formulate this problem as a constrained optimization problem and prove that an admissible optimal solution exists. Then, using quasi-Newton methods, we compute locally optimal potentials. These have the structure of a truncated periodic potential with a localized defect. In contrast to optimal structures for other spectral optimization problems, our optimizing potentials appear to be interior points of the constraint set and to be smooth. The multi-scale structures that emerge incorporate the physical mechanisms of energy confinement via material contrast and interference effects. An analysis of locally optimal potentials reveals local optimality is attained via two mechanisms: (i) decreasing the density of states near a resonant frequency in the continuum and (ii) tuning the oscillations of extended states to make $Γ[V]$, an oscillatory integral, small. Our approach achieves lifetimes, $\sim (ε^2Γ[V])^{-1}$, for locally optimal potentials with $Γ^{-1}\sim\mathcal{O}(10^{9})$ as compared with $Γ^{-1}\sim \mathcal{O}(10^{2})$ for a typical potential. Finally, we explore the performance of optimal potentials via simulations of the time-evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Equipartition of Mass in Nonlinear Schrödinger / Gross-Pitaevskii Equations

We study the infinite time dynamics of a class of nonlinear Schrödinger / Gross-Pitaevskii equations. In our previous paper, we prove the asymptotic stability of the nonlinear ground state in a general situation which admits degenerate neutral modes of arbitrary finite multiplicity, a typical situation in systems with symmetry. Neutral modes correspond to purely imaginary (neutrally stable) point spectrum of the linearization of the Hamiltonian PDE about a critical point. In particular, a small perturbation of the nonlinear ground state, which typically excites such neutral modes and radiation, will evolve toward an asymptotic nonlinear ground state soliton plus decaying neutral modes plus decaying radiation. In the present article, we give a much more detailed, in fact quantitative, picture of the asymptotic evolution. Specificially we prove an equipartition law: The asymptotic soliton which emerges has a mass which is equal to the initial soliton mass plus one half the mass contained in the initially perturbing neutral modes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Multiscale Model of Partial Melts 1: Effective Equations

In this paper a model for partial melts is constructed using two-scale homogenization theory. While this technique is well known to the mathematics and materials communities, it is relatively novel to problems in the solid Earth. This approach begins with a grain scale model of the medium, coarsening it into a macroscopic one. The emergent model is in good agreement with previous work, including D. McKenzie's, and serves as verification. This methodology also yields a series of Stokes problems whose solutions provide constitutive relations for permeability and viscosity. A numerical investigation of these relations appears in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Multiscale Model of Partial Melts 2: Numerical Results

In a companion paper, equations for partially molten media were derived using two-scale homogenization theory. One advantage of homogenization is that material properties, such as permeability and viscosity, readily emerge. A caveat is that the dependence of these parameters upon the microstructure is not self-evident. In particular, one seeks to relate them to the porosity. In this paper, we numerically solve ensembles of the cell problems from which these quantities emerge. Using this data, we estimate relationships between the parameters and the porosity. In particular, the bulk viscosity appears to be inversely proportional to the porosity. Finally, we synthesize these numerical estimates with the models. Our hybrid numerical--analytical model predicts that the compaction length vanishes with porosity.

0 citations0 reviewsNo code linkedOpen access