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

Gideon Simpson

Gideon Simpson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.APNumerical Analysisnlin.PSBiological Physicscond-mat.stat-mechphysics.chem-phphysics.comp-phphysics.geo-phphysics.optics

Trust snapshot

Quick read

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

preprint2022arXiv

A Numerical Method for a Nonlocal Diffusion Equation with Additive Noise

We consider a nonlocal evolution equation representing the continuum limit of a large ensemble of interacting particles on graphs forced by noise. The two principle ingredients of the continuum model are a nonlocal term and Q-Wiener process describing the interactions among the particles in the network and stochastic forcing respectively. The network connectivity is given by a square integrable function called a graphon. We prove that the initial value problem for the continuum model is well-posed. Further, we construct a semidiscrete (discrete in space and continuous in time) and a fully discrete schemes for the nonlocal model. The former is obtained by a discontinuous Galerkin method and the latter is based on further discretizing time using the Euler-Maruyama method. We prove convergence and estimate the rate of convergence in each case. For the semidiscrete scheme, the rate of convergence estimate is expressed in terms of the regularity of the graphon, Q-Wiener process, and the initial data. We work in generalized Lipschitz spaces, which allows to treat models with data of lower regularity. This is important for applications as many interesting types of connectivity including small-world and power-law are expressed by graphons that are not smooth. The error analysis of the fully discrete scheme, on the other hand, reveals that for some models common in applied science, one has a higher speed of convergence than that predicted by the standard estimates for the Euler-Maruyama method. The rate of convergence analysis is supplemented with detailed numerical experiments, which are consistent with our analytical results. As a by-product, this work presents a rigorous justification for taking continuum limit for a large class of interacting dynamical systems on graphs subject to noise.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Iterate Averaging, the Kalman Filter, and 3DVAR for Linear Inverse Problem

It has been proposed that classical filtering methods, like the Kalman filter and 3DVAR, can be used to solve linear statistical inverse problems. In the work of Iglesias, Lin, Lu, & Stuart (2017), error estimates were obtained for this approach. By optimally tuning a regularization parameter in the filters, the authors were able to show that the mean squared error could be systematically reduced. Building on the aforementioned work of Iglesias, Lin, Lu, & Stuart, we prove that by (i) considering the problem in a weaker norm and (ii) applying simple iterate averaging of the filter output, 3DVAR will converge in mean square, unconditionally on the choice of parameter. Without iterate averaging, 3DVAR cannot converge by running additional iterations with a fixed choice of parameter. We also establish that the Kalman filter's performance in this setting cannot be improved through iterate averaging. We illustrate our results with numerical experiments that suggest our convergence rates are sharp.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sampling from Rough Energy Landscapes

We examine challenges to sampling from Boltzmann distributions associated with multiscale energy landscapes. The multiscale features, or "roughness," corresponds to highly oscillatory, but bounded, perturbations of a smooth landscape. Through a combination of numerical experiments and analysis we demonstrate that the performance of Metropolis Adjusted Langevin Algorithm can be severely attenuated as the roughness increases. In contrast, we prove that Random Walk Metropolis is insensitive to such roughness. We also formulate two alternative sampling strategies that incorporate large scale features of the energy landscape, while resisting the impact of fine scale roughness; these also outperform Random Walk Metropolis. Numerical experiments on these landscapes are presented that confirm our predictions. Open questions and numerical challenges are also highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture

Probability currents are fundamental in characterizing the kinetics of non-equilibrium processes. Notably, the steady-state current $J_{ss}$ for a source-sink system can provide the exact mean-first-passage time (MFPT) for the transition from source to sink. Because transient non-equilibrium behavior is quantified in some modern path sampling approaches, such as the "weighted ensemble" strategy, there is strong motivation to determine bounds on $J_{ss}$ -- and hence on the MFPT -- as the system evolves in time. Here we show that $J_{ss}$ is bounded from above and below by the maximum and minimum, respectively, of the current as a function of the spatial coordinate at any time $t$ for one-dimensional systems undergoing over-damped Langevin (i.e., Smoluchowski) dynamics and for higher-dimensional Smoluchowski systems satisfying certain assumptions when projected onto a single dimension. These bounds become tighter with time, making them of potential practical utility in a scheme for estimating $J_{ss}$ and the long-timescale kinetics of complex systems. Conceptually, the bounds result from the fact that extrema of the transient currents relax toward the steady-state current.

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

Solitary Wave Benchmarks in Magma Dynamics

We present a model problem for benchmarking codes that investigate magma migration in the Earth's interior. This system retains the essential features of more sophisticated models, yet has the advantage of possessing solitary wave solutions. The existence of such exact solutions to the nonlinear problem make it an excellent benchmark problem for combinations of solver algorithms. In this work, we explore a novel algorithm for computing high quality approximations of the solitary waves and use them to benchmark a semi-Lagrangian Crank-Nicholson scheme for a finite element discretization of the time dependent problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Vortex Collapse for the L2-Critical Nonlinear Schrödinger Equation

The focusing cubic nonlinear Schrödinger equation in two dimensions admits vortex solitons, standing wave solutions with spatial structure, Qm(r,theta) = e^{i m theta} Rm(r). In the case of spin m = 1, we prove there exists a class of data that collapse with the vortex soliton profile at the log-log rate. This extends the work of Merle and Raphaël, (the case m = 0,) and suggests that the L2 mass that may be concentrated at a point during generic collapse may be unbounded. Difficulties with m >= 2 or when breaking the spin symmetry are discussed.

0 citations0 reviewsNo code linkedOpen access