1,971 research outputs found
Three real-space discretization techniques in electronic structure calculations
A characteristic feature of the state-of-the-art of real-space methods in
electronic structure calculations is the diversity of the techniques used in
the discretization of the relevant partial differential equations. In this
context, the main approaches include finite-difference methods, various types
of finite-elements and wavelets. This paper reports on the results of several
code development projects that approach problems related to the electronic
structure using these three different discretization methods. We review the
ideas behind these methods, give examples of their applications, and discuss
their similarities and differences.Comment: 39 pages, 10 figures, accepted to a special issue of "physica status
solidi (b) - basic solid state physics" devoted to the CECAM workshop "State
of the art developments and perspectives of real-space electronic structure
techniques in condensed matter and molecular physics". v2: Minor stylistic
and typographical changes, partly inspired by referee comment
PoisFFT - A Free Parallel Fast Poisson Solver
A fast Poisson solver software package PoisFFT is presented. It is available
as a free software licensed under the GNU GPL license version 3. The package
uses the fast Fourier transform to directly solve the Poisson equation on a
uniform orthogonal grid. It can solve the pseudo-spectral approximation and the
second order finite difference approximation of the continuous solution. The
paper reviews the mathematical methods for the fast Poisson solver and
discusses the software implementation and parallelization. The use of PoisFFT
in an incompressible flow solver is also demonstrated
Exploring a coarse-grained distributive strategy for finite-difference PoissonâBoltzmann calculations
We have implemented and evaluated a coarse-grained distributive method for finite-difference PoissonâBoltzmann (FDPB) calculations of large biomolecular systems. This method is based on the electrostatic focusing principle of decomposing a large fine-grid FDPB calculation into multiple independent FDPB calculations, each of which focuses on only a small and a specific portion (block) of the large fine grid. We first analyzed the impact of the focusing approximation upon the accuracy of the numerical reaction field energies and found that a reasonable relative accuracy of 10â3 can be achieved when the buffering space is set to be 16 grid points and the block dimension is set to be at least (1/6)3 of the fine-grid dimension, as in the one-block focusing method. The impact upon efficiency of the use of buffering space to maintain enough accuracy was also studied. It was found that an âoptimalâ multi-block dimension exists for a given computer hardware setup, and this dimension is more or less independent of the solute geometries. A parallel version of thedistributive focusing method was also implemented. Given the proper settings, the distributive method was able to achieve respectable parallel efficiency with tested biomolecular systems on a loosely connected computer cluster
The CECAM Electronic Structure Library and the modular software development paradigm
First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include ``heavy-duty'' ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers
A bibliography on parallel and vector numerical algorithms
This is a bibliography of numerical methods. It also includes a number of other references on machine architecture, programming language, and other topics of interest to scientific computing. Certain conference proceedings and anthologies which have been published in book form are listed also
Particle-Particle, Particle-Scaling function (P3S) algorithm for electrostatic problems in free boundary conditions
An algorithm for fast calculation of the Coulombic forces and energies of
point particles with free boundary conditions is proposed. Its calculation time
scales as N log N for N particles. This novel method has lower crossover point
with the full O(N^2) direct summation than the Fast Multipole Method. The
forces obtained by our algorithm are analytical derivatives of the energy which
guarantees energy conservation during a molecular dynamics simulation. Our
algorithm is very simple. An MPI parallelised version of the code can be
downloaded under the GNU General Public License from the website of our group.Comment: 19 pages, 11 figures, submitted to: Journal of Chemical Physic
Solution of partial differential equations on vector and parallel computers
The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed
Real-Space Mesh Techniques in Density Functional Theory
This review discusses progress in efficient solvers which have as their
foundation a representation in real space, either through finite-difference or
finite-element formulations. The relationship of real-space approaches to
linear-scaling electrostatics and electronic structure methods is first
discussed. Then the basic aspects of real-space representations are presented.
Multigrid techniques for solving the discretized problems are covered; these
numerical schemes allow for highly efficient solution of the grid-based
equations. Applications to problems in electrostatics are discussed, in
particular numerical solutions of Poisson and Poisson-Boltzmann equations.
Next, methods for solving self-consistent eigenvalue problems in real space are
presented; these techniques have been extensively applied to solutions of the
Hartree-Fock and Kohn-Sham equations of electronic structure, and to eigenvalue
problems arising in semiconductor and polymer physics. Finally, real-space
methods have found recent application in computations of optical response and
excited states in time-dependent density functional theory, and these
computational developments are summarized. Multiscale solvers are competitive
with the most efficient available plane-wave techniques in terms of the number
of self-consistency steps required to reach the ground state, and they require
less work in each self-consistency update on a uniform grid. Besides excellent
efficiencies, the decided advantages of the real-space multiscale approach are
1) the near-locality of each function update, 2) the ability to handle global
eigenfunction constraints and potential updates on coarse levels, and 3) the
ability to incorporate adaptive local mesh refinements without loss of optimal
multigrid efficiencies.Comment: 70 pages, 11 figures. To be published in Reviews of Modern Physic
- âŠ