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

Charles L. Epstein

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

ResearcherUnclaimed source record
math.APmath.NAPopulations and Evolutionmath.PRmath-phmath.MPNumerical Analysismath.SGeess.IVeess.SPhep-exmath.CVmath.DGmath.SPphysics.class-phphysics.plasm-ph

Catalog footprint

What is connected

18works
16topics
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

18 published item(s)

preprint2021arXiv

Debye source representation for type-I superconductors, I

In this note, we analyze the classical magneto-static approach to the theory of type I superconductors, and a Debye source representation that can be used numerically to solve the resultant equations. We also prove that one of the fields, $\boldsymbol{J}^-$, found within the superconductor via the London equations, is the physical current in that the outgoing part of the magnetic field is given as the Biot-Savart integral of $\boldsymbol{J}^{-}$. Finally, we compute the static currents for moderate values of London penetration depth, $λ_L,$ for a sphere, a stellarator-like geometry and a two-holed torus.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometry of the Phase Retrieval Problem

One of the most powerful approaches to imaging at the nanometer or subnanometer length scale is coherent diffraction imaging using X-ray sources. For amorphous (non-crystalline) samples, the raw data can be interpreted as the modulus of the continuous Fourier transform of the unknown object. Making use of prior information about the sample (such as its support), a natural goal is to recover the phase through computational means, after which the unknown object can be visualized at high resolution. While many algorithms have been proposed for this phase retrieval problem, careful analysis of its well-posedness has received relatively little attention. In this paper, we show that the problem is, in general, not well-posed and describe some of the underlying issues that are responsible for the ill-posedness. We then show how this analysis can be used to develop experimental protocols that lead to better conditioned inverse problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recovering missing data in coherent diffraction imaging

In coherent diffraction imaging (CDI) experiments, the intensity of the scattered wave impinging on an object is measured on an array of detectors. This signal can be interpreted as the square of the modulus of the Fourier transform of the unknown scattering density. A beam-stop obstructs the forward scattered wave and, hence, the modulus Fourier data from a neighborhood of k=0 cannot be measured. In this note, we describe a linear method for recovering this unmeasured modulus Fourier data from the measured values and an estimate of the support of the image's autocorrelation function without consideration of phase retrieval. We analyze the conditioning of this problem, which grows exponentially with the modulus of the maximum spatial frequency not measured, and the effects of noise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Eigenfunctions and the Dirichlet problem for the Classical Kimura Diffusion Operator

We study the classical Kimura diffusion operator defined on the n-simplex, $$L^{Kim}=\sum_{1\leq i,j\leq n+1}x_ix_j\partial_{x_i}\partial_{x_j}$$ We give novel constructions for the basis of eigenpolynomials, and the solution to the inhomogeneous Dirichlet problem, which are well adapted to numerical applications. Our solution of the Dirichlet problem is quite explicit and provides a precise description of the singularities that arise along the boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Smoothed corners and scattered waves

We introduce an arbitrary order, computationally efficient method to smooth corners on curves in the plane, as well as edges and vertices on surfaces in $\mathbb R^3$. The method is local, only modifying the original surface in a neighborhood of the geometric singularity, and preserves desirable features like convexity and symmetry. The smoothness of the final surface is an explicit parameter in the method, and the bandlimit of the smoothed surface is proportional to its smoothness. Several numerical examples are provided in the context of acoustic scattering. In particular, we compare scattered fields from smoothed geometries in two dimensions with those from polygonal domains. We observe that significant reductions in computational cost can be obtained if merely approximate solutions are desired in the near- or far-field. Provided that it is sub-wavelength, the error of the scattered field is proportional to the size of the geometry that is modified.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formal properties of the probability of fixation: identities, inequalities and approximations

The formula for the probability of fixation of a new mutation is widely used in theoretical population genetics and molecular evolution. Here we derive a series of identities, inequalities and approximations for the exact probability of fixation of a new mutation under the Moran process (equivalent results hold for the approximate probability of fixation for the Wright-Fisher process after an appropriate change of variables). We show that the behavior of the logarithm of the probability of fixation is particularly simple when the selection coefficient is measured as a difference of Malthusian fitnesses, and we exploit this simplicity to derive several inequalities and approximations. We also present a comprehensive comparison of both existing and new approximations for the probability of fixation, highlighting in particular approximations that result in a reversible Markov chain when used to model the dynamics of evolution under weak mutation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Harnack Inequalities and Heat-kernel Estimates for Degenerate Diffusion Operators Arising in Population Biology

This paper continues the analysis, started in [2, 3], of a class of degenerate elliptic operators defined on manifolds with corners, which arise in Population Biology. Using techniques pioneered by J. Moser, and extended and refined by L. Saloff-Coste, Grigoryan, and Sturm, we show that weak solutions to the parabolic problem defined by a sub-class of these operators, which consists of those that can be defined by Dirichlet forms and have non-vanishing transverse vector field, satisfy a Harnack inequality. This allows us to conclude that the solutions to these equations belong, for positive times, to the natural anisotropic Holder spaces, and also leads to upper and, in some cases, lower bounds for the heat kernels of these operators. These results imply that these operators have a compact resolvent when acting on $C^0$ or $L^2.$ The proof relies upon a scale invariant Poincare inequality that we establish for a large class of weighted Dirichlet forms, as well as estimates to handle certain mildly singular perturbation terms. The weights that we consider are neither Ahlfors regular, nor do they generally belong to the Muckenhaupt class $A_2.$

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Harnack Inequalities for Degenerate Diffusions

We study various probabilistic and analytical properties of a class of degenerate diffusion operators arising in Population Genetics, the so-called generalized Kimura diffusion operators. Our main results is a stochastic representation of weak solutions to a degenerate parabolic equation with singular lower-order coefficients, and the proof of the scale-invariant Harnack inequality for nonnegative solutions to the Kimura parabolic equation. The stochastic representation of solutions that we establish is a considerable generalization of the classical results on Feynman-Kac formulas concerning the assumptions on the degeneracy of the diffusion matrix, the boundedness of the drift coefficients, and on the a priori regularity of the weak solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Regularity for the supercritical fractional Laplacian with drift

We consider the linear stationary equation defined by the fractional Laplacian with drift. In the supercritical case, that is the case when the dominant term is given by the drift instead of the diffusion component, we prove local regularity of solutions in Sobolev spaces em- ploying tools from the theory of pseudo-differential operators. The regularity of solutions in the supercritical case is as expected in the subcritical case, when the diffusion is at least as strong as the drift component, and the operator defined by the fractional Laplacian with drift can be viewed as an elliptic operator, which is not the case in the supercritical regime. We compute the leading singularity for the Green's kernel in the supercritical range, which displays some unusual behavior: it is more singular in the half plane into which the drift vector points, than in the complementary half plane.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Debye Sources, Beltrami Fields, and a Complex Structure on Maxwell Fields

The Debye source representation for solutions to the time harmonic Maxwell equations is extended to bounded domains with finitely many smooth boundary components. A strong uniqueness result is proved for this representation. Natural complex structures are identified on the vector spaces of time-harmonic Maxwell fields. It is shown that in terms of Debye source data, these complex structures are uniformized, that is, represented by a fixed linear map on a fixed vector space, independent of the frequency. This complex structure relates time-harmonic Maxwell fields to constant-k Beltrami fields, i.e. solutions of the equation curl(E) = kE. A family of self-adjoint boundary conditions are defined for the Beltrami operator. This leads to a proof of the existence of zero-flux, constant-k, force-free Beltrami fields for any bounded region in R^3, as well as a constructive method to find them. The family of self-adjoint boundary value problems defines a new spectral invariant for bounded domains in R^3.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the convergence of local expansions of layer potentials

In a recently developed quadrature method (quadrature by expansion or QBX), it was demonstrated that weakly singular or singular layer potentials can be evaluated rapidly and accurately on surface by making use of local expansions about carefully chosen off-surface points. In this paper, we derive estimates for the rate of convergence of these local expansions, providing the analytic foundation for the QBX method. The estimates may also be of mathematical interest, particularly for microlocal or asymptotic analysis in potential theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The inevitability of unconditionally deleterious substitutions during adaptation

Studies on the genetics of adaptation typically neglect the possibility that a deleterious mutation might fix. Nonetheless, here we show that, in many regimes, the first substitution is most often deleterious, even when fitness is expected to increase in the long term. In particular, we prove that this phenomenon occurs under weak mutation for any house-of-cards model with an equilibrium distribution. We find that the same qualitative results hold under Fisher's geometric model. We also provide a simple intuition for the surprising prevalence of unconditionally deleterious substitutions during early adaptation. Importantly, the phenomenon we describe occurs on fitness landscapes without any local maxima and is therefore distinct from "valley-crossing". Our results imply that the common practice of ignoring deleterious substitutions leads to qualitatively incorrect predictions in many regimes. Our results also have implications for the substitution process at equilibrium and for the response to a sudden decrease in population size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A consistency condition for the vector potential in multiply-connected domains

A classical problem in electromagnetics concerns the representation of the electric and magnetic fields in the low-frequency or static regime, where topology plays a fundamental role. For multiply connected conductors, at zero frequency the standard boundary conditions on the tangential components of the magnetic field do not uniquely determine the vector potential. We describe a (gauge-invariant) consistency condition that overcomes this non-uniqueness and resolves a longstanding difficulty in inverting the magnetic field integral equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Debye Sources and the Numerical Solution of the Time Harmonic Maxwell Equations, II

In this paper, we develop a new integral representation for the solution of the time harmonic Maxwell equations in media with piecewise constant dielectric permittivity and magnetic permeability in R^3. This representation leads to a coupled system of Fredholm integral equations of the second kind for four scalar densities supported on the material interface. Like the classical Muller equation, it has no spurious resonances. Unlike the classical approach, however, the representation does not suffer from low frequency breakdown. We illustrate the performance of the method with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Degenerate Diffusion Operators Arising in Population Biology

We analyze a class of partial differential equations that arise as "backwards Kolmogorov operators" in infinite population limits of the Wright-Fisher models in population genetics and in mathematical finance. These are degenerate elliptic operators defined on manifolds with corners. The classical example is the Kimura diffusion operator, which acts on functions defined on the simplex in R^n. We introduce anisotropic Holder spaces, and prove existence, uniqueness and regularity results for the heat and resolvent equations defined by this class of operators. This suffices to prove that the C^0-graph closure generates a strongly continuous semigroup, and that the associated Martingale problem has a unique solution. We give a detailed description of the nullspace of the forward Kolmogorov operator.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Subelliptic Spin_C Dirac operators, II Basic Estimates

We assume that the manifold with boundary, X, has a Spin_C-structure with spinor bundle S. Along the boundary, this structure agrees with the structure defined by an infinite order integrable almost complex structure and the metric is Kahler. The induced CR-structure on bX is integrable and either strictly pseudoconvex or strictly pseudoconcave. We assume that E->X is a complex vector bundle, which has an infinite order integrable complex structure along bX, compatible with that defined along bX. In this paper use boundary layer methods to prove subelliptic estimates for the twisted Spin_C- Dirac operator acting on sections on S\otimes E. We use boundary conditions that are modifications of the classical dbar-Neumann condition. These results are proved by using the extended Heisenberg calculus.

0 citations0 reviewsNo code linkedOpen access