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

Eric Polizzi

Eric Polizzi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Numerical Analysisphysics.comp-phmath.NAMathematical Softwarecond-mat.mtrl-sciphysics.chem-phcond-mat.mes-hallComputational Engineering, Finance, and Sciencecond-mat.str-elmath.DSphysics.atm-clusquant-ph

Trust snapshot

Quick read

Trust 21 - Emerging
15works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

15 published item(s)

preprint2020arXiv

An Iterative Method for Contour-Based Nonlinear Eigensolvers

Contour integration techniques have become a popular choice for solving the linear and non-linear eigenvalue problems. They principally include the Sakurai-Sugiura methods, the Beyn's algorithm, the FEAST/NLFEAST algorithms and other rational filtering techniques. While these methods can result in effective 'black-box' approach for solving linear eigenvalue problems, they still present several shortcomings for addressing nonlinear eigenvalue problems which are both mathematically and practically far more challenging. In this paper, we introduce a new hybrid algorithm that advantageously combines the iterative nature of NLFEAST with the effectiveness of Beyn's approach to deal with general non-linearity. In doing so, this NLFEAST-Beyn hybrid algorithm can overcome current limitations of both algorithms taken separately. After presenting its derivation from both a Beyn's and NLFEAST's perspective, several numerical examples are discussed to demonstrate the efficiency of the new technique.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

FEAST Eigenvalue Solver v4.0 User Guide

The FEAST library package represents an unified framework for solving various family of eigenvalue problems and achieving accuracy, robustness, high-performance and scalability on parallel architectures. Its originality lies with a new transformative numerical approach to the traditional eigenvalue algorithm design - the FEAST algorithm. The algorithm gathers key elements from complex analysis, numerical linear algebra and approximation theory, to construct an optimal subspace iteration technique using approximate spectral projectors. FEAST can be used for solving both standard and generalized forms of the Hermitian or non-Hermitian problems (linear or non-linear), and it belongs to the family of contour integration eigensolvers. FEAST's main computational task consists of a numerical quadrature computation that involves solving independent linear systems along a complex contour, each with multiple right hand sides. In v4.0, FEAST has been reimplemented using an inverse residual iteration algorithm which enables the linear systems to be solved with very low accuracy (in single precision) with no impact on the FEAST double precision convergence rate. As a result, v4.0 is on average 3-4 times faster than v2.1 and v3.0 using new default optimization parameters (v2.1 has been featured as Intel-MKL's principal HPC eigensolver since 2013). v4.0 also implements new important features such as IFEAST (using Inexact Iterative solver), Non-linear polynomial FEAST, and PFEAST with its 3-MPI levels of parallelism. FEAST is both a comprehensive library package, and an easy to use software. It includes flexible reverse communication interfaces and ready to use driver interfaces for dense, banded and sparse systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enhancing the Performance and Robustness of the FEAST Eigensolver

The FEAST algorithm is a subspace iteration method that uses a spectral projector as a rational filter in order to efficiently solve interior eigenvalue problems in parallel. Although the solutions from the FEAST algorithm converge rapidly in many cases, convergence can be slow in situations where the eigenvalues of a matrix are densely populated near the edges of the search interval of interest, which can be detrimental to parallel load balancing. This work introduces two methods that allow one to improve the convergence robustness of the FEAST algorithm in these situations without having to increase the amount of computation. Selected numerical examples are presented and discussed

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

FEAST Eigensolver for non-Hermitian Problems

A detailed new upgrade of the FEAST eigensolver targeting non-Hermitian eigenvalue problems is presented and thoroughly discussed. It aims at broadening the class of eigenproblems that can be addressed within the framework of the FEAST algorithm. The algorithm is ideally suited for computing selected interior eigenvalues and their associated right/left bi-orthogonal eigenvectors,located within a subset of the complex plane. It combines subspace iteration with efficient contour integration techniques that approximate the left and right spectral projectors. We discuss the various algorithmic choices that have been made to improve the stability and usability of the new non-Hermitian eigensolver. The latter retains the convergence property and multi-level parallelism of Hermitian FEAST, making it a valuable new software tool for the scientific community.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

FEAST Eigenvalue Solver v3.0 User Guide

The FEAST eigensolver package is a free high-performance numerical library for solving the Hermitian and non-Hermitian eigenvalue problems, and obtaining all the eigenvalues and (right/left) eigenvectors within a given search interval or arbitrary contour in the complex plane. Its originality lies with a new transformative numerical approach to the traditional eigenvalue algorithm design - the FEAST algorithm. The FEAST eigensolver combines simplicity and efficiency and it offers many important capabilities for achieving high performance, robustness, accuracy, and scalability on parallel architectures. FEAST is both a comprehensive library package, and an easy to use software. It includes flexible reverse communication interfaces and ready to use predefined interfaces for dense, banded and sparse systems. The current version v3.0 of the FEAST package can address both Hermitian and non-Hermitian eigenvalue problems (real symmetric, real non-symmetric, complex Hermitian, complex symmetric, or complex general systems) on both shared-memory and distributed memory architectures (i.e contains both FEAST-SMP and FEAST-MPI packages). This User's guide provides instructions for installation setup, a detailed description of the FEAST interfaces and a large number of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Universal nature of collective plasmonic excitations in finite 1-D carbon-based nanostructures

Tomonaga-Luttinger (T-L) theory predicts collective plasmon resonances in 1-D nanostructure conductors of finite length, that vary roughly in inverse proportion to the length of the structure. In-depth quantitative understanding of such resonances which have not been clearly identified in experiments so far, would be invaluable for future generations of nano-photonic and nano-electronic devices that employ 1-D conductors. Here we provide evidence of the plasmon resonances in a number of representative 1-D finite carbon-based nanostructures using first-principle computational electronic spectroscopy studies. Our special purpose real-space/real-time all-electron Time-Dependent Density-Functional Theory (TDDFT) simulator can perform excited-states calculations to obtain correct frequencies for known optical transitions, and capture various nanoscopic effects including collective plasmon excitations. The presence of 1-D plasmons is universally predicted by the various numerical experiments, which also demonstrate a phenomenon of resonance splitting. For the metallic carbon nanotubes under study, the plasmons are expected to be related to the T-L plasmons of infinitely long 1-D structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A New Highly Parallel Non-Hermitian Eigensolver

Calculating portions of eigenvalues and eigenvectors of matrices or matrix pencils has many applications. An approach to this calculation for Hermitian problems based on a density matrix has been proposed in 2009 and a software package called FEAST has been developed. The density-matrix approach allows FEAST's implementation to exploit a key strength of modern computer architectures, namely, multiple levels of parallelism. Consequently, the software package has been well received and subsequently commercialized. A detailed theoretical analysis of Hermitian FEAST has also been established very recently. This paper generalizes the FEAST algorithm and theory, for the first time, to tackle non-Hermitian problems. Fundamentally, the new algorithm is basic subspace iteration or Bauer bi-iteration, except applied with a novel accelerator based on Cauchy integrals. The resulting algorithm retains the multi-level parallelism of Hermitian FEAST, making it a valuable new tool for large-scale computational science and engineering problems on leading-edge computing platforms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient estimation of eigenvalue counts in an interval

Estimating the number of eigenvalues located in a given interval of a large sparse Hermitian matrix is an important problem in certain applications and it is a prerequisite of eigensolvers based on a divide-and-conquer paradigm. Often an exact count is not necessary and methods based on stochastic estimates can be utilized to yield rough approximations. This paper examines a number of techniques tailored to this specific task. It reviews standard approaches and explores new ones based on polynomial and rational approximation filtering combined with a stochastic procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

FEAST as a Subspace Iteration Eigensolver Accelerated by Approximate Spectral Projection

The calculation of a segment of eigenvalues and their corresponding eigenvectors of a Hermitian matrix or matrix pencil has many applications. A new density-matrix-based algorithm has been proposed recently and a software package FEAST has been developed. The density-matrix approach allows FEAST's implementation to exploit a key strength of modern computer architectures, namely, multiple levels of parallelism. Consequently, the software package has been well received, especially in the electronic structure community. Nevertheless, theoretical analysis of FEAST has lagged. For instance, the FEAST algorithm has not been proven to converge. This paper offers a detailed numerical analysis of FEAST. In particular, we show that the FEAST algorithm can be understood as an accelerated subspace iteration algorithm in conjunction with the Rayleigh-Ritz procedure. The novelty of FEAST lies in its accelerator which is a rational matrix function that approximates the spectral projector onto the eigenspace in question. Analysis of the numerical nature of this approximate spectral projector and the resulting subspaces generated in the FEAST algorithm establishes the algorithm's convergence. This paper shows that FEAST is resilient against rounding errors and establishes properties that can be leveraged to enhance the algorithm's robustness. Finally, we propose an extension of FEAST to handle non-Hermitian problems and suggest some future research directions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Zolotarev Quadrature Rules and Load Balancing for the FEAST Eigensolver

The FEAST method for solving large sparse eigenproblems is equivalent to subspace iteration with an approximate spectral projector and implicit orthogonalization. This relation allows to characterize the convergence of this method in terms of the error of a certain rational approximant to an indicator function. We propose improved rational approximants leading to FEAST variants with faster convergence, in particular, when using rational approximants based on the work of Zolotarev. Numerical experiments demonstrate the possible computational savings especially for pencils whose eigenvalues are not well separated and when the dimension of the search space is only slightly larger than the number of wanted eigenvalues. The new approach improves both convergence robustness and load balancing when FEAST runs on multiple search intervals in parallel.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Discretization of Time-Dependent Quantum Systems: Real-Time Propagation of The Evolution Operator

We discuss time dependent quantum systems on bounded domains. Our work may be viewed as a framework for several models, including linear iterations involved in time dependent density functional theory (TDDFT), the Hartree-Fock model, or other quantum models. A key aspect of the analysis of the algorithms is the use of time-ordered evolution operators, which allow for both a well-posed problem and its approximation. The approximation theorems obtained for the time-ordered evolution operators complement those in the current literature. We discuss the available theory at the outset, and proceed to apply the theory systematically in later sections via approximations and a global existence theorem for a nonlinear system, obtained via a fixed point theorem for the evolution operator. Our work is consistent with first-principle real time propagation of electronic states, aimed at finding the electronic responses of quantum molecular systems and nanostructures. We present two full 3D quantum atomistic simulations using the finite element method for discretizing the real-space, and the FEAST eigenvalue algorithm for solving the evolution operator at each time step. These numerical experiments are representative of the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-linear eigensolver-based alternative to traditional SCF methods

The self-consistent procedure in electronic structure calculations is revisited using a highly efficient and robust algorithm for solving the non-linear eigenvector problem i.e. H({ψ})ψ = Eψ. This new scheme is derived from a generalization of the FEAST eigenvalue algorithm to account for the non-linearity of the Hamiltonian with the occupied eigenvectors. Using a series of numerical examples and the DFT-Kohn/Sham model, it will be shown that our approach can outperform the traditional SCF mixing-scheme techniques by providing a higher converge rate, convergence to the correct solution regardless of the choice of the initial guess, and a significant reduction of the eigenvalue solve time in simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

FEAST fundamental framework for electronic structure calculations: Reformulation and solution of the muffin-tin problem

In a recent article [1], the FEAST algorithm has been presented as a general purpose eigenvalue solver which is ideally suited for addressing the numerical challenges in electronic structure calculations. Here, FEAST is presented beyond the "black-box" solver as a fundamental modeling framework which can naturally address the original numerical complexity of the electronic structure problem as formulated by Slater in 1937 [2]. The non-linear eigenvalue problem arising from the muffin-tin decomposition of the real-space domain is first derived and then reformulated to be solved exactly within the FEAST framework. This new framework is presented as a fundamental and practical solution for performing both accurate and scalable electronic structure calculations, bypassing the various issues of using traditional approaches such as linearization and pseudopotential techniques. A finite element implementation of this FEAST framework along with simulation results for various molecular systems are also presented and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral-based Propagation Schemes for Time-Dependent Quantum Systems with Application to Carbon Nanotubes

Effective modeling and numerical spectral-based propagation schemes are proposed for addressing the challenges in time-dependent quantum simulations of systems ranging from atoms, molecules, and nanostructures to emerging nanoelectronic devices. While time-dependent Hamiltonian problems can be formally solved by propagating the solutions along tiny simulation time steps, a direct numerical treatment is often considered too computationally demanding. In this paper, however, we propose to go beyond these limitations by introducing high-performance numerical propagation schemes to compute the solution of the time-ordered evolution operator. In addition to the direct Hamiltonian diagonalizations that can be efficiently performed using the new eigenvalue solver FEAST, we have designed a Gaussian propagation scheme and a basis transformed propagation scheme (BTPS) which allow to reduce considerably the simulation times needed by time intervals. It is outlined that BTPS offers the best computational efficiency allowing new perspectives in time-dependent simulations. Finally, these numerical schemes are applied to study the AC response of a (5,5) carbon nanotube within a 3D real-space mesh framework.

0 citations0 reviewsNo code linkedOpen access