-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Orbital-enriched Flat-top Partition of Unity Method for the Schrödinger Eigenproblem
Authors:
Clelia Albrecht,
Constanze Klaar,
John E. Pask,
Marc Alexander Schweitzer,
N. Sukumar,
Albert Ziegenhagel
Abstract:
Quantum mechanical calculations require the repeated solution of a Schrödinger equation for the wavefunctions of the system. Recent work has shown that enriched finite element methods significantly reduce the degrees of freedom required to obtain accurate solutions. However, time to solution has been adversely affected by the need to solve a generalized eigenvalue problem and the ill-conditioning…
▽ More
Quantum mechanical calculations require the repeated solution of a Schrödinger equation for the wavefunctions of the system. Recent work has shown that enriched finite element methods significantly reduce the degrees of freedom required to obtain accurate solutions. However, time to solution has been adversely affected by the need to solve a generalized eigenvalue problem and the ill-conditioning of associated systems matrices. In this work, we address both issues by proposing a stable and efficient orbital-enriched partition-of-unity method to solve the Schrödinger boundary-value problem in a parallelepiped unit cell subject to Bloch-periodic boundary conditions. In our proposed PUM, the three-dimensional domain is covered by overlapping patches, with a compactly-supported, non-negative weight function, that is identically equal to unity over some finite subset of its support associated with each patch. This so-called flat-top property provides a pathway to devise a stable approximation over the whole domain. On each patch, we use $p$-th degree orthogonal polynomials that ensure $p$-th order completeness, and in addition include eigenfunctions of the radial solution of the Schrödinger equation. Furthermore, we adopt a variational lumping approach to construct a block-diagonal overlap matrix that yields a standard eigenvalue problem and demonstrate accuracy, stability and efficiency of the method.
△ Less
Submitted 1 March, 2018; v1 submitted 29 January, 2018;
originally announced January 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.
-
Chebyshev polynomial filtered subspace iteration in the Discontinuous Galerkin method for large-scale electronic structure calculations
Authors:
Amartya S. Banerjee,
Lin Lin,
Wei Hu,
Chao Yang,
John E. Pask
Abstract:
The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis (ALB) set to solve the Kohn-Sham equations of density functional theory (DFT) in a discontinuous Galerkin framework. The adaptive local basis is generated on-the-fly to capture the local material physics, and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. A…
▽ More
The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis (ALB) set to solve the Kohn-Sham equations of density functional theory (DFT) in a discontinuous Galerkin framework. The adaptive local basis is generated on-the-fly to capture the local material physics, and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. A central issue for large-scale calculations, however, is the computation of the electron density (and subsequently, ground state properties) from the discretized Hamiltonian in an efficient and scalable manner. We show in this work how Chebyshev polynomial filtered subspace iteration (CheFSI) can be used to address this issue and push the envelope in large-scale materials simulations in a discontinuous Galerkin framework. We describe how the subspace filtering steps can be performed in an efficient and scalable manner using a two-dimensional parallelization scheme, thanks to the orthogonality of the DG basis set and block-sparse structure of the DG Hamiltonian matrix. The on-the-fly nature of the ALBs requires additional care in carrying out the subspace iterations. We demonstrate the parallel scalability of the DG-CheFSI approach in calculations of large-scale two-dimensional graphene sheets and bulk three-dimensional lithium-ion electrolyte systems. Employing 55,296 computational cores, the time per self-consistent field iteration for a sample of the bulk 3D electrolyte containing 8,586 atoms is 90 seconds, and the time for a graphene sheet containing 11,520 atoms is 75 seconds.
△ Less
Submitted 29 September, 2016; v1 submitted 10 June, 2016;
originally announced June 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.
-
Lithium Ion Solvation and Diffusion in Bulk Organic Electrolytes from First Principles and Classical Reactive Molecular Dynamics
Authors:
Mitchell T. Ong,
Osvalds Verners,
Erik W. Draeger,
Adri C. T. van Duin,
Vincenzo Lordi,
John E. Pask
Abstract:
Lithium-ion battery performance is strongly influenced by the ionic conductivity of the electrolyte, which depends on the speed at which Li ions migrate across the cell and relates to their solvation structure. The choice of solvent can greatly impact both solvation and diffusivity of Li ions. We use first principles molecular dynamics to examine the solvation and diffusion of Li ions in the bulk…
▽ More
Lithium-ion battery performance is strongly influenced by the ionic conductivity of the electrolyte, which depends on the speed at which Li ions migrate across the cell and relates to their solvation structure. The choice of solvent can greatly impact both solvation and diffusivity of Li ions. We use first principles molecular dynamics to examine the solvation and diffusion of Li ions in the bulk organic solvents ethylene carbonate (EC), ethyl methyl carbonate (EMC), and a mixture of EC/EMC. We find that Li ions are solvated by either carbonyl or ether oxygen atoms of the solvents and sometimes by the PF$_6^-$ anion. Li$^+$ prefers a tetrahedrally-coordinated first solvation shell regardless of which species are involved, with the specific preferred solvation structure dependent on the organic solvent. In addition, we calculate Li diffusion coefficients in each electrolyte, finding slightly larger diffusivities in the linear carbonate EMC compared to the cyclic carbonate EC. The magnitude of the diffusion coefficient correlates with the strength of Li$^+$ solvation. Corresponding analysis for the PF$_6^-$ anion shows greater diffusivity associated with a weakly-bound, poorly defined first solvation shell. These results may be used to aid in the design of new electrolytes to improve Li-ion battery performance.
△ Less
Submitted 26 February, 2015;
originally announced February 2015.
-
Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations
Authors:
Yunfeng Cai,
Zhaojun Bai,
John E. Pask,
N. Sukumar
Abstract:
The iterative diagonalization of a sequence of large ill-conditioned generalized eigenvalue problems is a computational bottleneck in quantum mechanical methods employing a nonorthogonal basis for {\em ab initio} electronic structure calculations. We propose a hybrid preconditioning scheme to effectively combine global and locally accelerated preconditioners for rapid iterative diagonalization of…
▽ More
The iterative diagonalization of a sequence of large ill-conditioned generalized eigenvalue problems is a computational bottleneck in quantum mechanical methods employing a nonorthogonal basis for {\em ab initio} electronic structure calculations. We propose a hybrid preconditioning scheme to effectively combine global and locally accelerated preconditioners for rapid iterative diagonalization of such eigenvalue problems. In partition-of-unity finite-element (PUFE) pseudopotential density-functional calculations, employing a nonorthogonal basis, we show that the hybrid preconditioned block steepest descent method is a cost-effective eigensolver, outperforming current state-of-the-art global preconditioning schemes, and comparably efficient for the ill-conditioned generalized eigenvalue problems produced by PUFE as the locally optimal block preconditioned conjugate-gradient method for the well-conditioned standard eigenvalue problems produced by planewave methods.
△ Less
Submitted 11 August, 2013;
originally announced August 2013.
-
Linear scaling solution of the all-electron Coulomb problem in solids
Authors:
J. E. Pask,
N. Sukumar,
S. E. Mousavi
Abstract:
We present a linear scaling formulation for the solution of the all-electron Coulomb problem in crystalline solids. The resulting method is systematically improvable and well suited to large-scale quantum mechanical calculations in which the Coulomb potential and energy of a continuous electronic density and singular nuclear density are required. Linear scaling is achieved by introducing smooth, s…
▽ More
We present a linear scaling formulation for the solution of the all-electron Coulomb problem in crystalline solids. The resulting method is systematically improvable and well suited to large-scale quantum mechanical calculations in which the Coulomb potential and energy of a continuous electronic density and singular nuclear density are required. Linear scaling is achieved by introducing smooth, strictly local neutralizing densities to render nuclear interactions strictly local, and solving the remaining neutral Poisson problem for the electrons in real space. While the formulation includes singular nuclear potentials without smearing approximations, the required Poisson solution is in Sobolev space $H^1$, as required for convergence in the energy norm. We employ enriched finite elements, with enrichments from isolated atom solutions, for an efficient solution of the resulting Poisson problem in the interacting solid. We demonstrate the accuracy and convergence of the approach by direct comparison to standard Ewald sums for a lattice of point charges, and demonstrate the accuracy in all-electron quantum mechanical calculations with an application to crystalline diamond.
△ Less
Submitted 10 August, 2011; v1 submitted 11 April, 2010;
originally announced April 2010.
-
Half metallic digital ferromagnetic heterostructure composed of a $δ$-doped layer of Mn in Si
Authors:
M. C. Qian,
C. Y. Fong,
Kai Liu,
Warren E. Pickett,
J. E. Pask,
L. H. Yang
Abstract:
We propose and investigate the properties of a digital ferromagnetic heterostructure (DFH) consisting of a $δ$-doped layer of Mn in Si, using \textit{ab initio} electronic-structure methods. We find that (i) ferromagnetic order of the Mn layer is energetically favorable relative to antiferromagnetic, and (ii) the heterostructure is a two-dimensional half metallic system. The metallic behavior is…
▽ More
We propose and investigate the properties of a digital ferromagnetic heterostructure (DFH) consisting of a $δ$-doped layer of Mn in Si, using \textit{ab initio} electronic-structure methods. We find that (i) ferromagnetic order of the Mn layer is energetically favorable relative to antiferromagnetic, and (ii) the heterostructure is a two-dimensional half metallic system. The metallic behavior is contributed by three majority-spin bands originating from hybridized Mn-$d$ and nearest-neighbor Si-$p$ states, and the corresponding carriers are responsible for the ferromagnetic order in the Mn layer. The minority-spin channel has a calculated semiconducting gap of 0.25 eV. Analysis of the total and partial densities of states, band structure, Fermi surfaces and associated charge density reveals the marked two-dimensional nature of the half metallicity. The band lineup is found to be favorable for retaining the half metal character to near the Curie temperature ($T_{C}$). Being Si based and possibly having a high $T_{C}$ as suggested by an experiment on dilutely doped Mn in Si, the heterostructure may be of special interest for integration into mature Si technologies for spintronic applications.
△ Less
Submitted 28 November, 2005; v1 submitted 10 August, 2005;
originally announced August 2005.
-
Structural, Electronic, and Magnetic Properties of MnO
Authors:
J. E. Pask,
D. J. Singh,
I. I. Mazin,
C. S. Hellberg,
J. Kortus
Abstract:
We calculate the structural, electronic, and magnetic properties of MnO from first principles, using the full-potential linearized augmented planewave method, with both local-density and generalized-gradient approximations to exchange and correlation. We find the ground state to be of rhombohedrally distorted B1 structure with compression along the [111] direction, antiferromagnetic with type-II…
▽ More
We calculate the structural, electronic, and magnetic properties of MnO from first principles, using the full-potential linearized augmented planewave method, with both local-density and generalized-gradient approximations to exchange and correlation. We find the ground state to be of rhombohedrally distorted B1 structure with compression along the [111] direction, antiferromagnetic with type-II ordering, and insulating, consistent with experiment. We show that the distortion can be understood in terms of a Heisenberg model with distance dependent nearest-neighbor and next-nearest-neighbor couplings determined from first principles. Finally, we show that magnetic ordering can induce significant charge anisotropy, and give predictions for electric field gradients in the ground-state rhombohedrally distorted structure.
△ Less
Submitted 2 March, 2001; v1 submitted 18 December, 2000;
originally announced December 2000.
-
Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element approach
Authors:
J. E. Pask,
B. M. Klein,
C. Y. Fong,
P. A. Sterne
Abstract:
We present an approach to solid-state electronic-structure calculations based on the finite-element method. In this method, the basis functions are strictly local, piecewise polynomials. Because the basis is composed of polynomials, the method is completely general and its convergence can be controlled systematically. Because the basis functions are strictly local in real space, the method allow…
▽ More
We present an approach to solid-state electronic-structure calculations based on the finite-element method. In this method, the basis functions are strictly local, piecewise polynomials. Because the basis is composed of polynomials, the method is completely general and its convergence can be controlled systematically. Because the basis functions are strictly local in real space, the method allows for variable resolution in real space; produces sparse, structured matrices, enabling the effective use of iterative solution methods; and is well suited to parallel implementation. The method thus combines the significant advantages of both real-space-grid and basis-oriented approaches and so promises to be particularly well suited for large, accurate ab initio calculations. We develop the theory of our approach in detail, discuss advantages and disadvantages, and report initial results, including the first fully three-dimensional electronic band structures calculated by the method.
△ Less
Submitted 20 December, 2000; v1 submitted 20 March, 1999;
originally announced March 1999.