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

Benjamin Stamm

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

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.str-elphysics.comp-phquant-ph

Catalog footprint

What is connected

15works
5topics
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

15 published item(s)

preprint2026arXiv

Certified Model Order Reduction for parametric Hermitian eigenproblems

This article deals with the efficient and certified numerical approximation of the smallest eigenvalue and the associated eigenspace of a large-scale parametric Hermitian matrix. For this aim, we rely on projection-based model order reduction (MOR), i.e., we approximate the large-scale problem by projecting it onto a suitable subspace and reducing it to one of a much smaller dimension. Such a subspace is constructed by means of weak greedy-type strategies. After detailing the connections with the reduced basis method for source problems, we introduce a novel error estimate for the approximation error related to the eigenspace associated with the smallest eigenvalue. Since the difference between the second smallest and the smallest eigenvalue, the so-called spectral gap, is crucial for the reliability of the error estimate, we propose efficiently computable upper and lower bounds for higher eigenvalues and for the spectral gap, which enable the assembly of a subspace for the MOR approximation of the spectral gap. Based on that, a second subspace is then generated for the MOR approximation of the eigenspace associated with the smallest eigenvalue. We also provide efficiently computable conditions to ensure that the multiplicity of the smallest eigenvalue is fully captured in the reduced space. This work is motivated by a specific application: the repeated identifications of the states with minimal energy, the so-called ground states, of parametric quantum spin system models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surrogate models for quantum spin systems based on reduced order modeling

We present a methodology to investigate phase-diagrams of quantum models based on the principle of the reduced basis method (RBM). The RBM is built from a few ground-state snapshots, i.e., lowest eigenvectors of the full system Hamiltonian computed at well-chosen points in the parameter space of interest. We put forward a greedy-strategy to assemble such small-dimensional basis, i.e., to select where to spend the numerical effort needed for the snapshots. Once the RBM is assembled, physical observables required for mapping out the phase-diagram (e.g., structure factors) can be computed for any parameter value with a modest computational complexity, considerably lower than the one associated to the underlying Hilbert space dimension. We benchmark the method in two test cases, a chain of excited Rydberg atoms and a geometrically frustrated antiferromagnetic two-dimensional lattice model, and illustrate the accuracy of the approach. In particular, we find that the ground-manifold can be approximated to sufficient accuracy with a moderate number of basis functions, which increases very mildly when the number of microscopic constituents grows - in stark contrast to the exponential growth of the Hilbert space needed to describe each of the few snapshots. A combination of the presented RBM approach with other numerical techniques circumventing even the latter big cost, e.g., Tensor Network methods, is a tantalising outlook of this work.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Boundary integral equations for isotropic linear elasticity

This articles first investigates boundary integral operators for the three-dimensional isotropic linear elasticity of a biphasic model with piecewise constant Lamé coefficients in the form of a bounded domain of arbitrary shape surrounded by a background material. In the simple case of a spherical inclusion, the vector spherical harmonics consist of eigenfunctions of the single and double layer boundary operators and we provide their spectra. Further, in the case of many spherical inclusions with isotropic materials, each with its own set of Lamé parameters, we propose an integral equation and a subsequent Galerkin discretization using the vector spherical harmonics and apply the discretization to several numerical test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-center decomposition of molecular densities: a mathematical perspective

The aim of this paper is to analyze from a mathematical perspective some existing schemes to partition a molecular density into several atomic contributions, with a specific focus on Iterative Stockholder Atom (ISA) methods. We provide a unified mathematical framework to describe the latter family of methods and propose a new scheme, named L-ISA (for linear approximation of ISA). We prove several important mathematical properties of the ISA and L-ISA minimization problems and show that the so-called ISA algorithms can be viewed as alternating minimization schemes, which in turn enables us to obtain new convergence results for these numerical methods. Specific mathematical properties of the ISA decomposition for diatomic systems are also presented. We also review the basis-space oriented Distributed Multipole Analysis method, the mathematical formulation of which is also clarified. Different schemes are numerically compared on different molecules and we discuss the advantages and drawbacks of each approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Linear Scaling in Accuracy Numerical Method for Computing the Electrostatic Forces in the $N$-Body Dielectric Spheres Problem

This article deals with the efficient and accurate computation of the electrostatic forces between charged, spherical dielectric particles undergoing mutual polarisation. We use the spectral Galerkin boundary integral equation framework developed by Lindgren et al. (J. Comput. Phys. 371 (2018): 712-731) and subsequently analysed in two earlier contributions of the authors to propose a linear scaling in cost algorithm for the computation of the approximate forces. We establish exponential convergence of the method and derive error estimates for the approximate forces that do not explicitly depend on the number of dielectric particles $N$. Consequently, the proposed method requires only $\mathcal{O}(N)$ operations to compute the electrostatic forces acting on $N$ dielectric particles up to any given and fixed relative error.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Integral Equation Formulation of the $N$-Body Dielectric Spheres Problem. Part I: Numerical Analysis

In this article, we analyse an integral equation of the second kind that represents the solution of $N$ interacting dielectric spherical particles undergoing mutual polarisation. A traditional analysis can not quantify the scaling of the stability constants -- and thus the approximation error -- with respect to the number $N$ of involved dielectric spheres. We develop a new a priori error analysis that demonstrates $N$-independent stability of the continuous and discrete formulations of the integral equation. Consequently, we obtain convergence rates that are independent of $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Integral Equation Formulation of the $N$-body Dielectric Spheres Problem. Part II: Complexity Analysis

This article is the second in a series of two papers concerning the mathematical study of a boundary integral equation of the second kind that describes the interaction of $N$ dielectric spherical particles undergoing mutual polarisation. The first article presented the numerical analysis of the Galerkin method used to solve this boundary integral equation and derived $N$-independent convergence rates for the induced surface charges and total electrostatic energy. The current article will focus on computational aspects of the algorithm. We provide a convergence analysis of the iterative method used to solve the underlying linear system and show that the number of liner solver iterations required to obtain a solution is independent of $N$. Additionally, we present two linear scaling solution strategies for the computation of the approximate induced surface charges. Finally, we consider a series of numerical experiments designed to validate our theoretical results and explore the dependence of the numerical errors and computational cost of solving the underlying linear system on different system parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gradient Flow Finite Element Discretizations with Energy-Based Adaptivity for the Gross-Pitaevskii Equation

We present an effective adaptive procedure for the numerical approximation of the steady-state Gross-Pitaevskii equation. Our approach is solely based on energy minimization, and consists of a combination of gradient flow iterations and adaptive finite element mesh refinements. Numerical tests show that this strategy is able to provide highly accurate results, with optimal convergence rates with respect to the number of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters

This paper presents a posteriori error estimates for conforming numerical approximations of eigenvalue clusters of second-order self-adjoint elliptic linear operators with compact resolvent. Given a cluster of eigenvalues, we estimate the error in the sum of the eigenvalues, as well as the error in the eigenvectors represented through the density matrix, i.e., the orthogonal projector on the associated eigenspace. This allows us to deal with degenerate (multiple) eigenvalues within the framework. All the bounds are valid under the only assumption that the cluster is separated from the surrounding smaller and larger eigenvalues; we show how this assumption can be numerically checked. Our bounds are guaranteed and converge with the same speed as the exact errors. They can be turned into fully computable bounds as soon as an estimate on the dual norm of the residual is available, which is presented in two particular cases: the Laplace eigenvalue problem discretized with conforming finite elements, and a Schr{ö}dinger operator with periodic boundary conditions of the form $--$Δ$ + V$ discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems

We present the first systematic work for deriving a posteriori error estimates for general non-polynomial basis functions in an interior penalty discontinuous Galerkin (DG) formulation for solving eigenvalue problems associated with second order linear operators. Eigenvalue problems of such types play important roles in scientific and engineering applications, particularly in theoretical chemistry, solid state physics and material science. Based on the framework developed in [{\it L. Lin, B. Stamm, http://dx.doi.org/10.1051/m2an/2015069}] for second order PDEs, we develop residual type upper and lower bound error estimates for measuring the a posteriori error for eigenvalue problems. The main merit of our method is that the method is parameter-free, in the sense that all but one solution-dependent constants appearing in the upper and lower bound estimates are explicitly computable by solving local and independent eigenvalue problems, and the only non-computable constant can be reasonably approximated by a computable one without affecting the overall effectiveness of the estimates in practice. Compared to the PDE case, we find that a posteriori error estimators for eigenvalue problems must neglect certain terms, which involves explicitly the exact eigenvalues or eigenfunctions that are not accessible in numerical simulations. We define such terms carefully, and justify numerically that the neglected terms are indeed numerically high order terms compared to the computable estimators. Numerical results for a variety of problems in 1D and 2D demonstrate that both the upper bound and lower bound are effective for measuring the error of eigenvalues and eigenfunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part I: Second order linear PDE

We present the first systematic work for deriving a posteriori error estimates for general non-polynomial basis functions in an interior penalty discontinuous Galerkin (DG) formulation for solving second order linear PDEs. Our residual type upper and lower bound error estimates measure the error in the energy norm. The main merit of our method is that the method is parameter-free, in the sense that all but one solution-dependent constants appearing in the upper and lower bound estimates are explicitly computable by solving local eigenvalue problems, and the only non-computable constant can be reasonably approximated by a computable one without affecting the overall effectiveness of the estimates in practice. As a side product of our formulation, the penalty parameter in the interior penalty formulation can be automatically determined as well. We develop an efficient numerical procedure to compute the error estimators. Numerical results for a variety of problems in 1D and 2D demonstrate that both the upper bound and lower bound are effective.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An embedded corrector problem to approximate the homogenized coefficients of an elliptic equation

We consider a diffusion equation with highly oscillatory coefficients that admits a homogenized limit. As an alternative to standard corrector problems, we introduce here an embedded corrector problem, written as a diffusion equation in the whole space in which the diffusion matrix is uniform outside some ball of radius $R$. Using that problem, we next introduce three approximations of the homogenized coefficients. These approximations, which are variants of the standard approximations obtained using truncated (supercell) corrector problems, are shown to converge when $R \to \infty$. We also discuss efficient numerical methods to solve the embedded corrector problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparison of some Reduced Representation Approximations

In the field of numerical approximation, specialists considering highly complex problems have recently proposed various ways to simplify their underlying problems. In this field, depending on the problem they were tackling and the community that are at work, different approaches have been developed with some success and have even gained some maturity, the applications can now be applied to information analysis or for numerical simulation of PDE's. At this point, a crossed analysis and effort for understanding the similarities and the differences between these approaches that found their starting points in different backgrounds is of interest. It is the purpose of this paper to contribute to this effort by comparing some constructive reduced representations of complex functions. We present here in full details the Adaptive Cross Approximation (ACA) and the Empirical Interpolation Method (EIM) together with other approaches that enter in the same category.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A posteriori error estimates for the Electric Field Integral Equation on polyhedra

We present a residual-based a posteriori error estimate for the Electric Field Integral Equation (EFIE) on a bounded polyhedron. The EFIE is a variational equation formulated in a negative order Sobolev space on the surface of the polyhedron. We express the estimate in terms of square-integrable and thus computable quantities and derive global lower and upper bounds (up to oscillation terms).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Locally adaptive greedy approximations for anisotropic parameter reduced basis spaces

Reduced order models, in particular the reduced basis method, rely on empirically built and problem dependent basis functions that are constructed during an off-line stage. In the on-line stage, the precomputed problem-dependent solution space, that is spanned by the basis functions, can then be used in order to reduce the size of the computational problem. For complex problems, the number of basis functions required to guarantee a certain error tolerance can become too large in order to benefit computationally from the model reduction. To overcome this, the present work introduces a framework where local approximation spaces (in parameter space) are used to define the reduced order approximation in order to have explicit control over the on-line cost. This approach also adapts the local approximation spaces to local anisotropic behavior in the parameter space. We present the algorithm and numerous numerical tests.

0 citations0 reviewsNo code linkedOpen access