117,372 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
O(N) methods in electronic structure calculations
Linear scaling methods, or O(N) methods, have computational and memory
requirements which scale linearly with the number of atoms in the system, N, in
contrast to standard approaches which scale with the cube of the number of
atoms. These methods, which rely on the short-ranged nature of electronic
structure, will allow accurate, ab initio simulations of systems of
unprecedented size. The theory behind the locality of electronic structure is
described and related to physical properties of systems to be modelled, along
with a survey of recent developments in real-space methods which are important
for efficient use of high performance computers. The linear scaling methods
proposed to date can be divided into seven different areas, and the
applicability, efficiency and advantages of the methods proposed in these areas
is then discussed. The applications of linear scaling methods, as well as the
implementations available as computer programs, are considered. Finally, the
prospects for and the challenges facing linear scaling methods are discussed.Comment: 85 pages, 15 figures, 488 references. Resubmitted to Rep. Prog. Phys
(small changes
Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development
Atomic partial charges appear in the Coulomb term of many force-field models
and can be derived from electronic structure calculations with a myriad of
atoms-in-molecules (AIM) methods. More advanced models have also been proposed,
using the distributed nature of the electron cloud and atomic multipoles. In
this work, an electrostatic force field is defined through a concise
approximation of the electron density, for which the Coulomb interaction is
trivially evaluated. This approximate "pro-density" is expanded in a minimal
basis of atom-centered s-type Slater density functions, whose parameters are
optimized by minimizing the Kullback-Leibler divergence of the pro-density from
a reference electron density, e.g. obtained from an electronic structure
calculation. The proposed method, Minimal Basis Iterative Stockholder (MBIS),
is a variant of the Hirshfeld AIM method but it can also be used as a
density-fitting technique. An iterative algorithm to refine the pro-density is
easily implemented with a linear-scaling computational cost, enabling
applications to supramolecular systems. The benefits of the MBIS method are
demonstrated with systematic applications to molecular databases and extended
models of condensed phases. A comparison to 14 other AIM methods shows its
effectiveness when modeling electrostatic interactions. MBIS is also suitable
for rescaling atomic polarizabilities in the Tkatchenko-Sheffler scheme for
dispersion interactions.Comment: 61 pages, 12 figures, 2 table
Numerical methods for computing Casimir interactions
We review several different approaches for computing Casimir forces and
related fluctuation-induced interactions between bodies of arbitrary shapes and
materials. The relationships between this problem and well known computational
techniques from classical electromagnetism are emphasized. We also review the
basic principles of standard computational methods, categorizing them according
to three criteria---choice of problem, basis, and solution technique---that can
be used to classify proposals for the Casimir problem as well. In this way,
mature classical methods can be exploited to model Casimir physics, with a few
important modifications.Comment: 46 pages, 142 references, 5 figures. To appear in upcoming Lecture
Notes in Physics book on Casimir Physic
- …