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

Phanish Suryanarayana

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

ResearcherUnclaimed source record
physics.comp-phcond-mat.mtrl-sciphysics.chem-phphysics.app-ph

Catalog footprint

What is connected

10works
4topics
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

10 published item(s)

preprint2021arXiv

Real-space density kernel method for Kohn-Sham density functional theory calculations at high temperature

Kohn-Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. We present a density matrix based method for Kohn-Sham calculations at high temperature that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann-Feynman forces, and Hellmann-Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw-Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials; and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116,045 K from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bending moduli for thirty-two select atomic monolayers from first principles

We calculate bending moduli along the principal directions for thirty-two select atomic monolayers using ab initio Density Functional Theory (DFT). Specifically, considering representative materials from each of Groups IV, V, III-V monolayers, transition metal dichalcogenides, Group III monochalcogenides, Group IV monochalcogenides, and transition metal trichalcogenides, we utilize the recently developed Cyclic DFT method to calculate the bending moduli in the practically relevant but previously intractable low-curvature limit. We find that the moduli generally increase with thickness of the monolayer and that structures with a rectangular lattice are prone to a higher degree of anisotropy relative to those with a honeycomb lattice. We also find that exceptions to these trends are generally a consequence of unusually strong/weak bonding and/or significant structural relxation related effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for $\mathcal{O}(N)$ Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to $\mathcal{O}(N)$ methods. We then describe its evaluation within the $\mathcal{O}(N)$ infinite-cell Clenshaw-Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic as well as insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results, and show convergence with respect to mesh size to established planewave results. We employ the new formulation to compute the viscosity of hydrogen at a million kelvin from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SPARC: Simulation Package for Ab-initio Real-space Calculations

We present SPARC: Simulation Package for Ab-initio Real-space Calculations. SPARC can perform Kohn-Sham density functional theory calculations for isolated systems such as molecules as well as extended systems such as crystals and surfaces, in both static and dynamic settings. It is straightforward to install/use and highly competitive with state-of-the-art planewave codes, demonstrating comparable performance on a small number of processors and increasing advantages as the number of processors grows. Notably, SPARC brings solution times down to a few seconds for systems with $\mathcal{O}(100-500)$ atoms on large-scale parallel computers, outperforming planewave counterparts by an order of magnitude and more.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On preconditioning the self-consistent field iteration in real-space Density Functional Theory

We present a real-space formulation for isotropic Fourier-space preconditioners used to accelerate the self-consistent field iteration in Density Functional Theory calculations. Specifically, after approximating the preconditioner in Fourier space using a rational function, we express its real-space application in terms of the solution of sparse Helmholtz-type systems. Using the truncated-Kerker and Resta preconditioners as representative examples, we show that the proposed real-space method is both accurate and efficient, requiring the solution of a single linear system, while accelerating self-consistency to the same extent as its exact Fourier-space counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations

Pulay's Direct Inversion in the Iterative Subspace (DIIS) method is one of the most widely used mixing schemes for accelerating the self-consistent solution of electronic structure problems. In this work, we propose a simple generalization of DIIS in which Pulay extrapolation is performed at periodic intervals rather than on every self-consistent field iteration, and linear mixing is performed on all other iterations. We demonstrate through numerical tests on a wide variety of materials systems in the framework of density functional theory that the proposed generalization of Pulay's method significantly improves its robustness and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters

As the first component of SPARC (Simulation Package for Ab-initio Real-space Calculations), we present an accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory (DFT) for isolated clusters. Specifically, utilizing a local reformulation of the electrostatics, the Chebyshev polynomial filtered self-consistent field iteration, and a reformulation of the non-local component of the force, we develop a framework using the finite-difference representation that enables the efficient evaluation of energies and atomic forces to within the desired accuracies in DFT. Through selected examples consisting of a variety of elements, we demonstrate that SPARC obtains exponential convergence in energy and forces with domain size; systematic convergence in the energy and forces with mesh-size to reference plane-wave result at comparably high rates; forces that are consistent with the energy, both free from any noticeable `egg-box' effect; and accurate ground-state properties including equilibrium geometries and vibrational spectra. In addition, for systems consisting up to thousands of electrons, SPARC displays weak and strong parallel scaling behavior that is similar to well-established and optimized plane-wave implementations, but with a significantly reduced prefactor. Overall, SPARC represents an attractive alternative to plane-wave codes for practical DFT simulations of isolated clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory

We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state as well as forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Augmented Lagrangian formulation of Orbital-Free Density Functional Theory

We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic orbital-free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a sequence of minimization problems without any constraint, thereby making it amenable to powerful unconstrained optimization algorithms. Further, we develop a parallel implementation of this approach for the Thomas-Fermi-von Weizscaker (TFW) kinetic energy functional in the framework of higher-order finite-differences and the conjugate gradient method. With this implementation, we establish that the Augmented Lagrangian approach is highly competitive compared to the penalty and Lagrange multiplier methods. Additionally, we show that higher-order finite-differences represent a computationally efficient discretization for performing OF-DFT simulations. Overall, we demonstrate that the proposed formulation and implementation is both efficient and robust by studying selected examples, including systems consisting of thousands of atoms. We validate the accuracy of the computed energies and forces by comparing them with those obtained by existing plane-wave methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coarse-graining Kohn-Sham Density Functional Theory

We present a real-space formulation for coarse-graining Kohn-Sham Density Functional Theory that significantly speeds up the analysis of material defects without appreciable loss of accuracy. The approximation scheme consists of two steps. First, we develop a linear-scaling method that enables the direct evaluation of the electron density without the need to evaluate individual orbitals. We achieve this by performing Gauss quadrature over the spectrum of the linearized Hamiltonian operator appearing in each iteration of the self-consistent field method. Building on the linear-scaling method, we introduce a spatial approximation scheme resulting in a coarse-grained Density Functional Theory. The spatial approximation is adapted so as to furnish fine resolution where necessary and to coarsen elsewhere. This coarse-graining step enables the analysis of defects at a fraction of the original computational cost, without any significant loss of accuracy. Furthermore, we show that the coarse-grained solutions are convergent with respect to the spatial approximation. We illustrate the scope, versatility, efficiency and accuracy of the scheme by means of selected examples.

0 citations0 reviewsNo code linkedOpen access