-
Efficient real space formalism for hybrid density functionals
Authors:
Xin Jing,
Phanish Suryanarayana
Abstract:
We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn-Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly co…
▽ More
We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn-Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme, while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with literature. Overall, the current formalism provides an avenue for efficient real space DFT calculations with hybrid density functionals.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
Spectral scheme for atomic structure calculations in density functional theory
Authors:
Sayan Bhowmik,
John E. Pask,
Andrew J. Medford,
Phanish Suryanarayana
Abstract:
We present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadratur…
▽ More
We present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number $Z = 1$ to $83$.
△ Less
Submitted 6 June, 2024; v1 submitted 1 June, 2024;
originally announced June 2024.
-
Overcoming the chemical complexity bottleneck in on-the-fly machine learned molecular dynamics simulations
Authors:
Lucas R. Timmerman,
Shashikant Kumar,
Phanish Suryanarayana,
Andrew J. Medford
Abstract:
We develop a framework for on-the-fly machine learned force field molecular dynamics simulations based on the multipole featurization scheme that overcomes the bottleneck with the number of chemical elements. Considering bulk systems with up to 6 elements, we demonstrate that the number of density functional theory calls remains approximately independent of the number of chemical elements, in cont…
▽ More
We develop a framework for on-the-fly machine learned force field molecular dynamics simulations based on the multipole featurization scheme that overcomes the bottleneck with the number of chemical elements. Considering bulk systems with up to 6 elements, we demonstrate that the number of density functional theory calls remains approximately independent of the number of chemical elements, in contrast to the increase in the smooth overlap of atomic positions scheme.
△ Less
Submitted 12 June, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
Kohn-Sham accuracy from orbital-free density functional theory via $Δ$-machine learning
Authors:
Shashikant Kumar,
Xin Jing,
John E. Pask,
Andrew J. Medford,
Phanish Suryanarayana
Abstract:
We present a $Δ$-machine learning model for obtaining Kohn-Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn-Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations,…
▽ More
We present a $Δ$-machine learning model for obtaining Kohn-Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn-Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations, and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas-Fermi-von Weizs{ä}cker (TFW) orbital-free energies and forces by more than two orders of magnitude, but is also more accurate than MLFFs based solely on Kohn-Sham DFT, while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al$_{0.88}$Si$_{0.12}$, the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn-Sham study performed at an order of magnitude smaller length and time scales.
△ Less
Submitted 10 October, 2023;
originally announced October 2023.
-
Self-consistent convolutional density functional approximations: Application to adsorption at metal surfaces
Authors:
Sushree Jagriti Sahoo,
Qimen Xu,
Xiangyun Lei,
Daniel Staros,
Gopal R. Iyer,
Brenda Rubenstein,
Phanish Suryanarayana,
Andrew J. Medford
Abstract:
The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals is available, but nearly all are based on the hierarchy of inputs commonly referred to as "Jacob's ladder." This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron den…
▽ More
The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals is available, but nearly all are based on the hierarchy of inputs commonly referred to as "Jacob's ladder." This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBE$α$ framework where $α$ is a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.
△ Less
Submitted 26 February, 2024; v1 submitted 9 August, 2023;
originally announced August 2023.
-
Strain engineering of Zeeman and Rashba effects in transition metal dichalcogenide nanotubes and their Janus variants: An ab initio study
Authors:
Arpit Bhardwaj,
Phanish Suryanarayana
Abstract:
We study the influence of mechanical deformations on the Zeeman and Rashba effects in synthesized transition metal dichalcogenide (TMD) nanotubes and their Janus variants from first principles. In particular, we perform symmetry-adapted density functional theory simulations with spin-orbit coupling to determine the variation in the Zeeman and Rashba splittings with axial and torsional deformations…
▽ More
We study the influence of mechanical deformations on the Zeeman and Rashba effects in synthesized transition metal dichalcogenide (TMD) nanotubes and their Janus variants from first principles. In particular, we perform symmetry-adapted density functional theory simulations with spin-orbit coupling to determine the variation in the Zeeman and Rashba splittings with axial and torsional deformations. We find significant splitting in molybdenum and tungsten nanotubes, for which the Zeeman splitting decreases with increase in strain, going to zero for large enough tensile/shear strains, while the Rashba splitting coefficient increases linearly with shear strain, while being zero for all tensile strains, a consequence of the inversion symmetry remaining unbroken. In addition, the Zeeman splitting is relatively unaffected by nanotube diameter, whereas the Rashba coefficient decreases with increase in diameter. Overall, mechanical deformations represent a powerful tool for spintronics in TMD nanotubes as well as their Janus variants.
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
On the Effective Mass of Mechanical Lattices with Microstructure
Authors:
Francesco Fedele,
Phanish Suryanarayana,
Arash Yavari
Abstract:
We present a general formalism for the analysis of mechanical lattices with microstructure using the concept of effective mass. We first revisit a classical case of microstructure being modeled by a spring-interconnected mass-in-mass cell. The frequency-dependent effective mass of the cell is the sum of a static mass and of an added mass, in analogy to that of a swimmer in a fluid. The effective m…
▽ More
We present a general formalism for the analysis of mechanical lattices with microstructure using the concept of effective mass. We first revisit a classical case of microstructure being modeled by a spring-interconnected mass-in-mass cell. The frequency-dependent effective mass of the cell is the sum of a static mass and of an added mass, in analogy to that of a swimmer in a fluid. The effective mass is derived using three different methods: momentum equivalence, action equivalence, and dynamic condensation. These methods are generalized to mechanical systems with arbitrary microstructure. As an application, we calculate the effective mass of a $1$D composite lattice with microstructure modeled by a chiral spring-interconnected mass-in-mass cell. A reduced (condensed) model of the full lattice is then obtained by lumping the microstructure into a single effective mass. A dynamic Bloch analysis is then performed using both the full and reduced lattice models, which give the same spectral results. In particular, the frequency bands follow from the full lattice model by solving a linear eigenvalue problem, or from the reduced lattice model by solving a smaller nonlinear eigenvalue problem. The range of frequencies of negative effective mass falls within the bandgaps of the lattice. Localized modes due to defects in the microstructure have frequencies within the bandgaps, inside the negative-mass range. Defects of the outer, or macro stiffness yield localized modes within each bandgap, but outside the negative-mass range. The proposed formalism can be applied to study the odd properties of coupled micro-macro systems, e.g., active matter.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
Version 2.0.0 -- SPARC: Simulation Package for Ab-initio Real-space Calculations
Authors:
Boqin Zhang,
Xin Jing,
Qimen Xu,
Shashikant Kumar,
Abhiraj Sharma,
Lucas Erlandson,
Sushree Jagriti Sahoo,
Edmond Chow,
Andrew J. Medford,
John E. Pask,
Phanish Suryanarayana
Abstract:
SPARC is an accurate, efficient, and scalable real-space electronic structure code for performing ab initio Kohn-Sham density functional theory calculations. Version 2.0.0 of the software provides increased efficiency, and includes spin-orbit coupling, dispersion interactions, and advanced semilocal as well as hybrid exchange-correlation functionals, where it outperforms state-of-the-art planewave…
▽ More
SPARC is an accurate, efficient, and scalable real-space electronic structure code for performing ab initio Kohn-Sham density functional theory calculations. Version 2.0.0 of the software provides increased efficiency, and includes spin-orbit coupling, dispersion interactions, and advanced semilocal as well as hybrid exchange-correlation functionals, where it outperforms state-of-the-art planewave codes by an order of magnitude and more, with increasing advantages as the number of processors is increased. These new features further expand the range of physical applications amenable to first principles investigation.
△ Less
Submitted 1 May, 2024; v1 submitted 12 May, 2023;
originally announced May 2023.
-
Towards understanding flexoelectricity at the nanoscale
Authors:
David Codony,
Phanish Suryanarayana,
Irene Arias
Abstract:
We review the authors' recent works on flexoelectricity at the nanoscale [arXiv:2010.01747, arXiv:2010.13899], while emphasizing the role of continuum mechanics in interpreting the electromechanical response of quantum mechanical systems under bending.
We review the authors' recent works on flexoelectricity at the nanoscale [arXiv:2010.01747, arXiv:2010.13899], while emphasizing the role of continuum mechanics in interpreting the electromechanical response of quantum mechanical systems under bending.
△ Less
Submitted 2 April, 2023;
originally announced April 2023.
-
Assessing the source of error in the Thomas-Fermi-von Weizsäcker density functional
Authors:
Bishal Thapa,
Xin Jing,
John E. Pask,
Phanish Suryanarayana,
Igor I. Mazin
Abstract:
We investigate the source of error in the Thomas-Fermi-von Weizsäcker (TFW) density functional relative to Kohn-Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn-Sham d…
▽ More
We investigate the source of error in the Thomas-Fermi-von Weizsäcker (TFW) density functional relative to Kohn-Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn-Sham density, the corresponding energy deviates significantly from the Kohn-Sham value. We show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn-Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn-Sham solution.
△ Less
Submitted 28 June, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
-
Version 2.0.0 -- M-SPARC: Matlab-Simulation Package for Ab-initio Real-space Calculations
Authors:
Boqin Zhang,
Xin Jing,
Shashikant Kumar,
Phanish Suryanarayana
Abstract:
M-SPARC is a Matlab code for performing ab initio Kohn-Sham Density Functional Theory simulations. Version 2.0.0 of the software further extends its capability to include relativistic effects, dispersion interactions, and advanced semilocal/nonlocal exchange-correlation functionals. These features significantly increase the fidelity of first principles calculations that can be performed using M-SP…
▽ More
M-SPARC is a Matlab code for performing ab initio Kohn-Sham Density Functional Theory simulations. Version 2.0.0 of the software further extends its capability to include relativistic effects, dispersion interactions, and advanced semilocal/nonlocal exchange-correlation functionals. These features significantly increase the fidelity of first principles calculations that can be performed using M-SPARC.
△ Less
Submitted 14 December, 2022; v1 submitted 12 December, 2022;
originally announced December 2022.
-
Ab initio study on the electromechanical response of Janus transition metal dihalide nanotubes
Authors:
Arpit Bhardwaj,
Phanish Suryanarayana
Abstract:
We study the electronic response of Janus transition metal dihalide (TMH) nanotubes to mechanical deformations using Kohn-Sham density functional theory. Specifically, considering twelve armchair and zigzag Janus TMH nanotubes that are expected to be stable from the phonon analysis of flat monolayer counterparts, we first compute their equilibrium diameters and then determine the variation in band…
▽ More
We study the electronic response of Janus transition metal dihalide (TMH) nanotubes to mechanical deformations using Kohn-Sham density functional theory. Specifically, considering twelve armchair and zigzag Janus TMH nanotubes that are expected to be stable from the phonon analysis of flat monolayer counterparts, we first compute their equilibrium diameters and then determine the variation in bandgap and effective mass of charge carriers with the application of tensile and torsional deformations. We find that the nanotubes undergo a linear and quadratic decrease in bandgap with tensile and shear strain, respectively. In addition, there is a continual increase and decrease in the effective mass of electrons and holes, respectively. We show that for a given strain, the change in bandgap for the armchair nanotubes can be correlated with the transition metal's in-plane $d$ orbital's contribution to the projected density of states at the bottom of the conduction band.
△ Less
Submitted 23 October, 2022;
originally announced October 2022.
-
Roadmap on Electronic Structure Codes in the Exascale Era
Authors:
Vikram Gavini,
Stefano Baroni,
Volker Blum,
David R. Bowler,
Alexander Buccheri,
James R. Chelikowsky,
Sambit Das,
William Dawson,
Pietro Delugas,
Mehmet Dogan,
Claudia Draxl,
Giulia Galli,
Luigi Genovese,
Paolo Giannozzi,
Matteo Giantomassi,
Xavier Gonze,
Marco Govoni,
Andris Gulans,
François Gygi,
John M. Herbert,
Sebastian Kokott,
Thomas D. Kühne,
Kai-Hsin Liou,
Tsuyoshi Miyazaki,
Phani Motamarri
, et al. (16 additional authors not shown)
Abstract:
Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources d…
▽ More
Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
On the bending of rectangular atomic monolayers along different directions: an ab initio study
Authors:
Shashikant Kumar,
Phanish Suryanarayana
Abstract:
We study the bending of rectangular atomic monolayers along different directions from first principles. Specifically, choosing the phosphorene, GeS, TiS$_3$, and As$_2$S$_3$ monolayers as representative examples, we perform Kohn-Sham density functional theory calculations to determine the variation in transverse flexoelectric coefficient and bending modulus with the direction of bending. We find t…
▽ More
We study the bending of rectangular atomic monolayers along different directions from first principles. Specifically, choosing the phosphorene, GeS, TiS$_3$, and As$_2$S$_3$ monolayers as representative examples, we perform Kohn-Sham density functional theory calculations to determine the variation in transverse flexoelectric coefficient and bending modulus with the direction of bending. We find that while the flexoelectric coefficient is nearly isotropic, there is significant and complex anisotropy in bending modulus that also differs between the monolayers, with extremal values not necessarily occurring along the principal directions. In particular, the commonly adopted orthotropic continuum plate model with uniform thickness fails to describe the observed variations in bending modulus for GeS, TiS$_3$, and As$_2$S$_3$. We determine the direction-dependent effective thickness for use in such continuum models. We also show that the anisotropy in bending modulus is not associated with the rehybridization of atomic orbitals.
△ Less
Submitted 28 September, 2022; v1 submitted 29 July, 2022;
originally announced August 2022.
-
A note on transversal flexoelectricity in two-dimensional systems
Authors:
David Codony,
Irene Arias,
Phanish Suryanarayana
Abstract:
In a recent letter [Phys. Rev. Lett. 127, 216801 (2021)], the authors introduced the effective flexoelectric coefficient $μ^\textrm{2D}$ for quantifying the flexoelectric effect in 2D systems, and reported a disagreement with the flexoelectric coefficient $μ_\text{T}$ introduced in Ref. [Phys. Rev. Mat. 5, L030801 (2021)], attributed to the neglect of the metric term $\varphi^\text{M}$ -- quadrupo…
▽ More
In a recent letter [Phys. Rev. Lett. 127, 216801 (2021)], the authors introduced the effective flexoelectric coefficient $μ^\textrm{2D}$ for quantifying the flexoelectric effect in 2D systems, and reported a disagreement with the flexoelectric coefficient $μ_\text{T}$ introduced in Ref. [Phys. Rev. Mat. 5, L030801 (2021)], attributed to the neglect of the metric term $\varphi^\text{M}$ -- quadrupolar moment of the unperturbed charge density -- in the formulation of Ref. [Phys. Rev. Mat. 5, L030801 (2021)]. Here, we show that the model in Ref. [Phys. Rev. Mat. 5, L030801 (2021)] is correct and is in agreement with that in Ref. [Phys. Rev. Lett. 127, 216801 (2021)]. The discrepancies in the numerical values of the coefficients arise due to the difference in their definitions: $μ_\textrm{T}$ measures changes in bending-induced out-of-plane polarization, whereas $μ^\textrm{2D}$ measures changes in bending-induced voltage drop across the 2D system.
△ Less
Submitted 26 April, 2023; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Strain engineering of Janus transition metal dichalcogenide nanotubes: An ab intio study
Authors:
Arpit Bhardwaj,
Phanish Suryanarayana
Abstract:
We study the electromechanical response of Janus transition metal dichalcogenide (TMD) nanotubes from first principles. In particular, considering both armchair and zigzag variants of twenty-seven select Janus TMD nanotubes, we determine the change in bandgap and charge carriers' effective mass upon axial and torsional deformations using density functional theory (DFT). We observe that metallic na…
▽ More
We study the electromechanical response of Janus transition metal dichalcogenide (TMD) nanotubes from first principles. In particular, considering both armchair and zigzag variants of twenty-seven select Janus TMD nanotubes, we determine the change in bandgap and charge carriers' effective mass upon axial and torsional deformations using density functional theory (DFT). We observe that metallic nanotubes remain unaffected, whereas the bandgap in semiconducting nanotubes decreases linearly and quadratically with axial and shear strains, respectively, leading to semiconductor--metal transitions. In addition, we find that there is a continuous decrease and increase in the effective mass of holes and electrons with strains, respectively, leading to n-type--p-type semiconductor transitions. We show that this behavior is a consequence of the rehybridization of orbitals, rather than charge transfer between the atoms. Overall, mechanical deformations form a powerful tool for tailoring the electronic response of semiconducting Janus TMD nanotubes.
△ Less
Submitted 28 September, 2022; v1 submitted 20 January, 2022;
originally announced January 2022.
-
Accurate approximations of density functional theory for large systems with applications to defects in crystalline solids
Authors:
Kaushik Bhattacharya,
Vikram Gavini,
Michael Ortiz,
Mauricio Ponga,
Phanish Suryanarayana
Abstract:
This chapter presents controlled approximations of Kohn-Sham density functional theory (DFT) that enable very large scale simulations. The work is motivated by the study of defects in crystalline solids, though the ideas can be used in other applications. The key idea is to formulate DFT as a minimization problem over the density operator, and to cast spatial and spectral discretization as systema…
▽ More
This chapter presents controlled approximations of Kohn-Sham density functional theory (DFT) that enable very large scale simulations. The work is motivated by the study of defects in crystalline solids, though the ideas can be used in other applications. The key idea is to formulate DFT as a minimization problem over the density operator, and to cast spatial and spectral discretization as systematically convergent approximations. This enables efficient and adaptive algorithms that solve the equations of DFT with no additional modeling, and up to desired accuracy, for very large systems, with linear and sublinear scaling. Various approaches based on such approximations are presented, and their numerical performance demonstrated through selected examples. These examples also provide important insight about the mechanics and physics of defects in crystalline solids.
△ Less
Submitted 11 December, 2021;
originally announced December 2021.
-
Elastic properties of Janus transition metal dichalcogenide nanotubes from first principles
Authors:
Arpit Bhardwaj,
Phanish Suryanarayana
Abstract:
We calculate the elastic properties of Janus transition metal dichalcogenide (TMD) nanotubes using first principles Kohn-Sham density functional theory (DFT). Specifically, we perform electronic structure simulations that exploit the cyclic and helical symmetry in the system to compute the Young's moduli, Poisson's ratios, and torsional moduli for twenty-seven select armchair and zigzag Janus TMD…
▽ More
We calculate the elastic properties of Janus transition metal dichalcogenide (TMD) nanotubes using first principles Kohn-Sham density functional theory (DFT). Specifically, we perform electronic structure simulations that exploit the cyclic and helical symmetry in the system to compute the Young's moduli, Poisson's ratios, and torsional moduli for twenty-seven select armchair and zigzag Janus TMD nanotubes at their equilibrium diameters. We find the following trend in the moduli values: MSSe > MSTe > MSeTe, while their anisotropy with respect to armchair and zigzag configurations has the following ordering: MSTe > MSeTe > MSSe. This anisotropy and its ordering between the different groups is confirmed by computing the shear modulus from the torsional modulus using an isotropic elastic continuum model, and comparing it against the value predicted from the isotropic relation featuring the Young's modulus and Poisson's ratio. We also develop a model for the Young's and torsional moduli of Janus TMD nanotubes based on linear regression.
△ Less
Submitted 28 September, 2022; v1 submitted 11 November, 2021;
originally announced November 2021.
-
Ab-initio investigation of finite size effects in rutile titania nanoparticles with semilocal and nonlocal density functionals
Authors:
Sushree Jagriti Sahoo,
Xin Jing,
Phanish Suryanarayana,
Andrew J. Medford
Abstract:
In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and band structure of rutile titania (TiO$_2$) nanoparticles with particle size. The surface energies and band structures are calculated for cuboidal particles with minimum dimension ranging from 3.7 Å (24 atoms) to 10.3…
▽ More
In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and band structure of rutile titania (TiO$_2$) nanoparticles with particle size. The surface energies and band structures are calculated for cuboidal particles with minimum dimension ranging from 3.7 Å (24 atoms) to 10.3 Å (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on band structure and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Å. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the band structure is highly dependent on the shape and size of these particles. The DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the bandgap play a key role in determining the band structure of nanoparticles and the bandgap does not converge to the bulk limit for the particle sizes investigated.
△ Less
Submitted 7 July, 2023; v1 submitted 5 October, 2021;
originally announced October 2021.
-
Implementation of Perdew-Zunger self-interaction correction in real space using Fermi-Löwdin orbitals
Authors:
Carlos M. Diaz,
Phanish Suryanarayana,
Qimen Xu,
Tunna Baruah,
John Pask,
Rajendra Zope
Abstract:
Most widely used density functional approximations suffer from self-interaction (SI) error, which can be corrected using the Perdew-Zunger (PZ) self-interaction correction (SIC). We implement the recently proposed size-extensive formulation of PZ-SIC using Fermi-Löwdin Orbitals (FLOs) in real space, which is amenable to systematic convergence and large-scale parallelization. We verify the new form…
▽ More
Most widely used density functional approximations suffer from self-interaction (SI) error, which can be corrected using the Perdew-Zunger (PZ) self-interaction correction (SIC). We implement the recently proposed size-extensive formulation of PZ-SIC using Fermi-Löwdin Orbitals (FLOs) in real space, which is amenable to systematic convergence and large-scale parallelization. We verify the new formulation within the generalized Slater scheme by computing atomization energies and ionization potentials of selected molecules and comparing to those obtained by existing FLOSIC implementations in Gaussian based codes. The results show good agreement between the two formulations, with new real-space results somewhat closer to experiment on average for the systems considered. We also obtain the ionization potentials and atomization energies by scaling down the Slater statistical average of SIC potentials. The results show that scaling down the average SIC potential improves both atomization energies and ionization potentials, bringing them closer to experiment. Finally, we verify the present formulation by calculating the barrier heights of chemical reactions in the BH6 dataset, where significant improvements are obtained relative to Gaussian based FLOSIC results.
△ Less
Submitted 2 August, 2021;
originally announced August 2021.
-
Torsional strain engineering of transition metal dichalcogenide nanotubes: An ab initio study
Authors:
Arpit Bhardwaj,
Abhiraj Sharma,
Phanish Suryanarayana
Abstract:
We study the effect of torsional deformations on the electronic properties of single-walled transition metal dichalcogenide (TMD) nanotubes. In particular, considering forty-five select armchair and zigzag TMD nanotubes, we perform symmetry-adapted Kohn-Sham density functional theory calculations to determine the variation in bandgap and effective mass of charge carriers with twist. We find that m…
▽ More
We study the effect of torsional deformations on the electronic properties of single-walled transition metal dichalcogenide (TMD) nanotubes. In particular, considering forty-five select armchair and zigzag TMD nanotubes, we perform symmetry-adapted Kohn-Sham density functional theory calculations to determine the variation in bandgap and effective mass of charge carriers with twist. We find that metallic nanotubes remain so even after deformation, whereas semiconducting nanotubes experience a decrease in bandgap with twist -- originally direct bandgaps become indirect -- resulting in semiconductor to metal transitions. In addition, the effective mass of holes and electrons continuously decrease and increase with twist, respectively, resulting in n-type to p-type semiconductor transitions. We find that this behavior is likely due to rehybridization of orbitals in the metal and chalcogen atoms, rather than charge transfer between them. Overall, torsional deformations represent a powerful avenue to engineer the electronic properties of semiconducting TMD nanotubes, with applications to devices like sensors and semiconductor switches.
△ Less
Submitted 8 June, 2021;
originally announced June 2021.
-
Torsional moduli of transition metal dichalcogenide nanotubes from first principles
Authors:
Arpit Bhardwaj,
Abhiraj Sharma,
Phanish Suryanarayana
Abstract:
We calculate the torsional moduli of single-walled transition metal dichalcogenide (TMD) nanotubes using ab initio density functional theory (DFT). Specifically, considering forty-five select TMD nanotubes, we perform symmetry-adapted DFT calculations to calculate the torsional moduli for the armchair and zigzag variants of these materials in the low-twist regime and at practically relevant diamet…
▽ More
We calculate the torsional moduli of single-walled transition metal dichalcogenide (TMD) nanotubes using ab initio density functional theory (DFT). Specifically, considering forty-five select TMD nanotubes, we perform symmetry-adapted DFT calculations to calculate the torsional moduli for the armchair and zigzag variants of these materials in the low-twist regime and at practically relevant diameters. We find that the torsional moduli follow the trend: MS$_2$ $>$ MSe$_2$ $>$ MTe$_2$. In addition, the moduli display a power law dependence on diameter, with the scaling generally close to cubic, as predicted by the isotropic elastic continuum model. In particular, the shear moduli so computed are in good agreement with those predicted by the isotropic relation in terms of the Young's modulus and Poisson's ratio, both of which are also calculated using symmetry-adapted DFT. Finally, we develop a linear regression model for the torsional moduli of TMD nanotubes based on the nature/characteristics of the metal-chalcogen bond, and show that it is capable of making reasonably accurate predictions.
△ Less
Submitted 15 March, 2021;
originally announced March 2021.
-
Transversal flexoelectric coefficients for fifty select atomic monolayers from first principles
Authors:
Shashikant Kumar,
David Codony,
Irene Arias,
Phanish Suryanarayana
Abstract:
We calculate transversal flexoelectric coefficients along the principal directions for fifty select atomic monolayers using ab initio Density Functional Theory (DFT). Specifically, considering representative materials from each of Groups IV, III-V, V monolayers, transition metal dichalcogenides (TMDs), Group III monochalcogenides, Group IV monochalcogenides, transition metal trichalcogenides (TMTs…
▽ More
We calculate transversal flexoelectric coefficients along the principal directions for fifty select atomic monolayers using ab initio Density Functional Theory (DFT). Specifically, considering representative materials from each of Groups IV, III-V, V monolayers, transition metal dichalcogenides (TMDs), Group III monochalcogenides, Group IV monochalcogenides, transition metal trichalcogenides (TMTs), and Group V chalcogenides, we perform symmetry-adapted DFT simulations to calculate the flexoelectric coefficients at practically relevant bending curvatures. We find that the materials demonstrate linear behavior and have similar coefficients along both principal directions, with values for TMTs being up to a factor of five larger than graphene. In addition, we find electronic origins for the flexoelectric effect, which increases with monolayer thickness, elastic modulus along bending direction, and sum of polarizability of constituent atoms.
△ Less
Submitted 26 October, 2020;
originally announced October 2020.
-
Transversal flexoelectric coefficient for nanostructures at finite deformations from first principles
Authors:
David Codony,
Irene Arias,
Phanish Suryanarayana
Abstract:
We present a novel formulation for calculating the transversal flexoelectric coefficient of nanostructures at finite deformations from first principles. Specifically, we introduce the concept of \emph{radial polarization} to make the coefficient a well-defined quantity for uniform bending deformations. We use the framework to calculate the flexoelectric coefficient for group IV atomic monolayers u…
▽ More
We present a novel formulation for calculating the transversal flexoelectric coefficient of nanostructures at finite deformations from first principles. Specifically, we introduce the concept of \emph{radial polarization} to make the coefficient a well-defined quantity for uniform bending deformations. We use the framework to calculate the flexoelectric coefficient for group IV atomic monolayers using density functional theory. We find that graphene's coefficient is significantly larger than previously reported, with a charge transfer mechanism that differs from other members of its group.
△ Less
Submitted 4 October, 2020;
originally announced October 2020.
-
SPARC: Simulation Package for Ab-initio Real-space Calculations
Authors:
Qimen Xu,
Abhiraj Sharma,
Benjamin Comer,
Hua Huang,
Edmond Chow,
Andrew J. Medford,
John E. Pask,
Phanish Suryanarayana
Abstract:
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 comparabl…
▽ More
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.
△ Less
Submitted 20 May, 2020;
originally announced May 2020.
-
Bending moduli for thirty-two select atomic monolayers from first principles
Authors:
Shashikant Kumar,
Phanish Suryanarayana
Abstract:
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 de…
▽ More
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.
△ Less
Submitted 9 June, 2020; v1 submitted 24 March, 2020;
originally announced March 2020.
-
M-SPARC: MATLAB-Simulation Package for Ab-initio Real-space Calculations
Authors:
Qimen Xu,
Abhiraj Sharma,
Phanish Suryanarayana
Abstract:
We present M-SPARC: MATLAB-Simulation Package for Ab-initio Real-space Calculations. It can perform pseudopotential spin-polarized and unpolarized Kohn-Sham Density Functional Theory (DFT) simulations for isolated systems such as molecules as well as extended systems such as crystals, surfaces, and nanowires. M-SPARC provides a rapid prototyping platform for the development and testing of new algo…
▽ More
We present M-SPARC: MATLAB-Simulation Package for Ab-initio Real-space Calculations. It can perform pseudopotential spin-polarized and unpolarized Kohn-Sham Density Functional Theory (DFT) simulations for isolated systems such as molecules as well as extended systems such as crystals, surfaces, and nanowires. M-SPARC provides a rapid prototyping platform for the development and testing of new algorithms and methods in real-space DFT, with the potential to significantly accelerate the rate of advancements in the field. It also provides a convenient avenue for the accurate first principles study of small to moderate sized systems.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
On preconditioning the self-consistent field iteration in real-space Density Functional Theory
Authors:
Shashikant Kumar,
Qimen Xu,
Phanish Suryanarayana
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 6 November, 2019; v1 submitted 15 October, 2019;
originally announced October 2019.
-
Symmetry-adapted real-space density functional theory for cylindrical geometries: application to large X (X=C, Si, Ge, Sn) nanotubes
Authors:
Swarnava Ghosh,
Amartya S. Banerjee,
Phanish Suryanarayana
Abstract:
We present a symmetry-adapted real-space formulation of Kohn-Sham density functional theory for cylindrical geometries and apply it to the study of large X (X=C, Si, Ge, Sn) nanotubes. Specifically, starting from the Kohn-Sham equations posed on all of space, we reduce the problem to the fundamental domain by incorporating cyclic and periodic symmetries present in the angular and axial directions…
▽ More
We present a symmetry-adapted real-space formulation of Kohn-Sham density functional theory for cylindrical geometries and apply it to the study of large X (X=C, Si, Ge, Sn) nanotubes. Specifically, starting from the Kohn-Sham equations posed on all of space, we reduce the problem to the fundamental domain by incorporating cyclic and periodic symmetries present in the angular and axial directions of the cylinder, respectively. We develop a high-order finite-difference parallel implementation of this formulation, and verify its accuracy against established planewave and real-space codes. Using this implementation, we study the band structure and bending properties of X nanotubes and Xene sheets, respectively. Specifically, we first show that zigzag and armchair X nanotubes with radii in the range 1 to 5 nm are semiconducting. In particular, we find an inverse linear dependence of the bandgap with respect to the radius for all nanotubes, other than the armchair and zigzag type III carbon variants, for which we find an inverse quadratic dependence. Next, we exploit the connection between cyclic symmetry and uniform bending deformations to calculate the bending moduli of Xene sheets in both zigzag and armchair directions. We find Kirchhoff-Love type bending behavior for all sheets, with graphene and stanene possessing the largest and smallest moduli, respectively. In addition, other than graphene, the sheets demonstrate significant anisotropy, with larger bending moduli along the armchair direction. Finally, we demonstrate that the proposed approach has very good parallel scaling and is highly efficient, enabling ab initio simulations of unprecedented size for systems with a high degree of cyclic symmetry. In particular, we show that even micron-sized nanotubes can be simulated with modest computational effort.
△ Less
Submitted 15 September, 2019; v1 submitted 30 April, 2019;
originally announced April 2019.
-
Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
Authors:
Shuai Zhang,
Amy Lazicki,
Burkhard Militzer,
Lin H. Yang,
Kyle Caspersen,
Jim A. Gaffney,
Markus W. Däne,
John E. Pask,
Walter R. Johnson,
Abhiraj Sharma,
Phanish Suryanarayana,
Duane D. Johnson,
Andrey V. Smirnov,
Philip A. Sterne,
David Erskine,
Richard A. London,
Federica Coppari,
Damian Swift,
Joseph Nilsen,
Art J. Nelson,
Heather D. Whitley
Abstract:
The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functio…
▽ More
The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (ACTEX), and all-electron Green's function Korringa-Kohn-Rostoker (MECCA), and compute the pressure and internal energy of BN over a broad range of densities ($ρ$) and temperatures ($T$). Our experiments were conducted at the Omega laser facility and measured the Hugoniot of BN to unprecedented pressures (12--30 Mbar). The EOSs computed using different methods cross validate one another, and the experimental Hugoniot are in good agreement with our theoretical predictions. We assess that the largest discrepancies between theoretical predictions are $<$4% in pressure and $<$3% in energy and occur at $10^6$ K. We find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the warm-dense regime. Moreover, SQ and FOE data have significantly smaller error bars than PIMC, and so represent significant advances for efficient computation at high $T$. We also construct tabular EOS models and clarify the ionic and electronic structure of BN over a broad $T-ρ$ range and quantify their roles in the EOS. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.
△ Less
Submitted 2 February, 2019;
originally announced February 2019.
-
On the calculation of the stress tensor in real-space Kohn-Sham Density Functional Theory
Authors:
Abhiraj Sharma,
Phanish Suryanarayana
Abstract:
We present an accurate and efficient formulation of the stress tensor for real-space Kohn-Sham Density Functional Theory (DFT) calculations. Specifically, while employing a local formulation of the electrostatics, we derive a linear-scaling expression for the stress tensor that is applicable to simulations with unit cells of arbitrary symmetry, semilocal exchange-correlation functionals, and Brill…
▽ More
We present an accurate and efficient formulation of the stress tensor for real-space Kohn-Sham Density Functional Theory (DFT) calculations. Specifically, while employing a local formulation of the electrostatics, we derive a linear-scaling expression for the stress tensor that is applicable to simulations with unit cells of arbitrary symmetry, semilocal exchange-correlation functionals, and Brillouin zone integration. In particular, we rewrite the contributions arising from the self energy and the nonlocal pseudopotential energy to make them amenable to the real-space finite-difference discretization, achieving up to three orders of magnitude improvement in the accuracy of the computed stresses. Using examples representative of static and dynamic calculations, we verify the accuracy and efficiency of the proposed formulation. In particular, we demonstrate high rates of convergence with spatial discretization, consistency between the computed energy and stress tensor, and very good agreement with reference planewave results.
△ Less
Submitted 21 September, 2018;
originally announced September 2018.
-
Discrete discontinuous basis projection method for large-scale electronic structure calculations
Authors:
Qimen Xu,
Phanish Suryanarayana,
John E. Pask
Abstract:
We present an approach to accelerate real-space electronic structure methods several fold, without loss of accuracy, by reducing the dimension of the discrete eigenproblem that must be solved. To accomplish this, we construct an efficient, systematically improvable, discontinuous basis spanning the occupied subspace and project the real-space Hamiltonian onto the span. In calculations on a range o…
▽ More
We present an approach to accelerate real-space electronic structure methods several fold, without loss of accuracy, by reducing the dimension of the discrete eigenproblem that must be solved. To accomplish this, we construct an efficient, systematically improvable, discontinuous basis spanning the occupied subspace and project the real-space Hamiltonian onto the span. In calculations on a range of systems, we find that accurate energies and forces are obtained with 8--25 basis functions per atom, reducing the dimension of the associated real-space eigenproblems by 1--3 orders of magnitude.
△ Less
Submitted 27 April, 2018;
originally announced April 2018.
-
On real-space Density Functional Theory for non-orthogonal crystal systems: Kronecker product formulation of the kinetic energy operator
Authors:
Abhiraj Sharma,
Phanish Suryanarayana
Abstract:
We present an accurate and efficient real-space Density Functional Theory (DFT) framework for the ab-initio study of non-orthogonal crystal systems. Specifically, employing a local reformulation of the electrostatics, we develop a novel Kronecker product formulation of the real-space kinetic energy operator that significantly reduces the number of operations associated with the Laplacian-vector mu…
▽ More
We present an accurate and efficient real-space Density Functional Theory (DFT) framework for the ab-initio study of non-orthogonal crystal systems. Specifically, employing a local reformulation of the electrostatics, we develop a novel Kronecker product formulation of the real-space kinetic energy operator that significantly reduces the number of operations associated with the Laplacian-vector multiplication, the dominant cost in practical computations. In particular, we reduce the scaling with respect to finite-difference order from quadratic to linear, thereby significantly bridging the gap in computational cost between non-orthogonal and orthogonal systems. We verify the accuracy and efficiency of the proposed methodology through selected examples.
△ Less
Submitted 11 March, 2018; v1 submitted 3 February, 2018;
originally announced February 2018.
-
Two-level Chebyshev filter based complementary subspace method: pushing the envelope of large-scale electronic structure calculations
Authors:
Amartya S. Banerjee,
Lin Lin,
Phanish Suryanarayana,
Chao Yang,
John E. Pask
Abstract:
We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (> 1000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to: 1) compute a set…
▽ More
We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (> 1000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to: 1) compute a set of vectors that span the occupied subspace of the Hamiltonian; 2) reduce subspace diagonalization to just partially occupied states; and 3) obtain those states in an efficient, scalable manner via an inner Chebyshev-filter iteration. By reducing the necessary computation to just partially occupied states, and obtaining these through an inner Chebyshev iteration, our approach reduces the cost of large metallic calculations significantly, while eliminating subspace diagonalization for insulating systems altogether. We describe the implementation of the method within the framework of the Discontinuous Galerkin (DG) electronic structure method and show that this results in a computational scheme that can effectively tackle bulk and nano systems containing tens of thousands of electrons, with chemical accuracy, within a few minutes or less of wall clock time per SCF iteration on large-scale computing platforms. We anticipate that our method will be instrumental in pushing the envelope of large-scale ab initio molecular dynamics. As a demonstration of this, we simulate a bulk silicon system containing 8,000 atoms at finite temperature, and obtain an average SCF step wall time of 51 seconds on 34,560 processors; thus allowing us to carry out 1.0 ps of ab initio molecular dynamics in approximately 28 hours (of wall time).
△ Less
Submitted 20 February, 2018; v1 submitted 12 December, 2017;
originally announced December 2017.
-
SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
Authors:
Phanish Suryanarayana,
Phanisri P. Pratapa,
Abhiraj Sharma,
John E. Pask
Abstract:
We present SQDFT: a large-scale parallel implementation of the Spectral Quadrature (SQ) method for $\mathcal{O}(N)$ Kohn-Sham Density Functional Theory (DFT) calculations at high temperature. Specifically, we develop an efficient and scalable finite-difference implementation of the infinite-cell Clenshaw-Curtis SQ approach, in which results for the infinite crystal are obtained by expressing quant…
▽ More
We present SQDFT: a large-scale parallel implementation of the Spectral Quadrature (SQ) method for $\mathcal{O}(N)$ Kohn-Sham Density Functional Theory (DFT) calculations at high temperature. Specifically, we develop an efficient and scalable finite-difference implementation of the infinite-cell Clenshaw-Curtis SQ approach, in which results for the infinite crystal are obtained by expressing quantities of interest as bilinear forms or sums of bilinear forms, that are then approximated by spatially localized Clenshaw-Curtis quadrature rules. We demonstrate the accuracy of SQDFT by showing systematic convergence of energies and atomic forces with respect to SQ parameters to reference diagonalization results, and convergence with discretization to established planewave results, for both metallic and insulating systems. We further demonstrate that SQDFT achieves excellent strong and weak parallel scaling on computer systems consisting of tens of thousands of processors, with near perfect $\mathcal{O}(N)$ scaling with system size and wall times as low as a few seconds per self-consistent field iteration. Finally, we verify the accuracy of SQDFT in large-scale quantum molecular dynamics simulations of aluminum at high temperature.
△ Less
Submitted 19 November, 2017; v1 submitted 25 August, 2017;
originally announced August 2017.
-
On nearsightedness in metallic systems for $\mathcal{O}(N)$ Density Functional Theory calculations: A case study on Aluminum
Authors:
Phanish Suryanarayana
Abstract:
We investigate the locality of electronic interactions in aluminum as a function of smearing/electronic temperature in the context of $\mathcal{O}(N)$ Density Functional Theory calculations. Specifically, we determine the convergence in energy and atomic forces with truncation region size for smearing of $0.001-0.15$ Ha. We find exponential convergence accompanied by a rate that increases sub-line…
▽ More
We investigate the locality of electronic interactions in aluminum as a function of smearing/electronic temperature in the context of $\mathcal{O}(N)$ Density Functional Theory calculations. Specifically, we determine the convergence in energy and atomic forces with truncation region size for smearing of $0.001-0.15$ Ha. We find exponential convergence accompanied by a rate that increases sub-linearly with smearing, with truncation region sizes of $48-64$ Bohr required to achieve chemical accuracy for typical smearing values of $0.001-0.01$ Ha. This translates to $\mathcal{O}(N)$ scaling for systems larger than $\mathcal{O}(1000)$ atoms.
△ Less
Submitted 24 April, 2017; v1 submitted 5 March, 2017;
originally announced March 2017.
-
Cyclic Density Functional Theory : A route to the first principles simulation of bending in nanostructures
Authors:
Amartya S. Banerjee,
Phanish Suryanarayana
Abstract:
We formulate and implement Cyclic Density Functional Theory (Cyclic DFT) -- a self-consistent first principles simulation method for nanostructures with cyclic symmetries. Using arguments based on Group Representation Theory, we rigorously demonstrate that the Kohn-Sham eigenvalue problem for such systems can be reduced to a fundamental domain (or cyclic unit cell) augmented with cyclic-Bloch boun…
▽ More
We formulate and implement Cyclic Density Functional Theory (Cyclic DFT) -- a self-consistent first principles simulation method for nanostructures with cyclic symmetries. Using arguments based on Group Representation Theory, we rigorously demonstrate that the Kohn-Sham eigenvalue problem for such systems can be reduced to a fundamental domain (or cyclic unit cell) augmented with cyclic-Bloch boundary conditions. Analogously, the equations of electrostatics appearing in Kohn-Sham theory can be reduced to the fundamental domain augmented with cyclic boundary conditions. By making use of this symmetry cell reduction, we show that the electronic ground-state energy and the Hellmann-Feynman forces on the atoms can be calculated using quantities defined over the fundamental domain. We develop a symmetry-adapted finite-difference discretization scheme to obtain a fully functional numerical realization of the proposed approach. We verify that our formulation and implementation of Cyclic DFT is both accurate and efficient through selected examples.
The connection of cyclic symmetries with uniform bending deformations provides an elegant route to the ab-initio study of bending in nanostructures using Cyclic DFT. As a demonstration of this capability, we simulate the uniform bending of a silicene nanoribbon and obtain its energy-curvature relationship from first principles. A self-consistent ab-initio simulation of this nature is unprecedented and well outside the scope of any other systematic first principles method in existence. Our simulations reveal that the bending stiffness of the silicene nanoribbon is intermediate between that of graphene and molybdenum disulphide. We describe several future avenues and applications of Cyclic DFT, including its extension to the study of non-uniform bending deformations and its possible use in the study of the nanoscale flexoelectric effect.
△ Less
Submitted 15 March, 2017; v1 submitted 28 May, 2016;
originally announced May 2016.
-
Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations
Authors:
Amartya S. Banerjee,
Phanish Suryanarayana,
John E. Pask
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 January, 2016; v1 submitted 4 December, 2015;
originally announced December 2015.
-
Coarse-graining Kohn-Sham Density Functional Theory
Authors:
Phanish Suryanarayana,
Kaushik Bhattacharya,
Michael Ortiz
Abstract:
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 achiev…
▽ More
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.
△ Less
Submitted 8 August, 2012;
originally announced August 2012.