6,290 research outputs found
A spectral scheme for Kohn-Sham density functional theory of clusters
Starting from the observation that one of the most successful methods for
solving the Kohn-Sham equations for periodic systems -- the plane-wave method
-- is a spectral method based on eigenfunction expansion, we formulate a
spectral method designed towards solving the Kohn-Sham equations for clusters.
This allows for efficient calculation of the electronic structure of clusters
(and molecules) with high accuracy and systematic convergence properties
without the need for any artificial periodicity. The basis functions in this
method form a complete orthonormal set and are expressible in terms of
spherical harmonics and spherical Bessel functions. Computation of the occupied
eigenstates of the discretized Kohn-Sham Hamiltonian is carried out using a
combination of preconditioned block eigensolvers and Chebyshev polynomial
filter accelerated subspace iterations. Several algorithmic and computational
aspects of the method, including computation of the electrostatics terms and
parallelization are discussed. We have implemented these methods and algorithms
into an efficient and reliable package called ClusterES (Cluster Electronic
Structure). A variety of benchmark calculations employing local and non-local
pseudopotentials are carried out using our package and the results are compared
to the literature. Convergence properties of the basis set are discussed
through numerical examples. Computations involving large systems that contain
thousands of electrons are demonstrated to highlight the efficacy of our
methodology. The use of our method to study clusters with arbitrary point group
symmetries is briefly discussed.Comment: Manuscript submitted (with revisions) to Journal of Computational
Physic
Model ab initio study of charge carrier solvation and large polaron formation on conjugated carbon chains
Using long C_{N}H_{2} conjugated carbon chains with the polyynic structure as
prototypical examples of one-dimensional (1D) semiconductors, we discuss
self-localization of excess charge carriers into 1D large polarons in the
presence of the interaction with a surrounding polar solvent. The solvation
mechanism of self-trapping is different from the polaron formation due to
coupling with bond-length modulations of the underlying atomic lattice
well-known in conjugated polymers. Model ab initio computations employing the
hybrid B3LYP density functional in conjunction with the polarizable continuum
model are carried out demonstrating the formation of both electron- and
hole-polarons. Polarons can emerge entirely due to solvation but even larger
degrees of charge localization occur when accompanied by atomic displacements
A Real-Space Full Multigrid study of the fragmentation of Li11+ clusters
We have studied the fragmentation of Li11+ clusters into the two
experimentally observed products (Li9+,Li2) and (Li10+,Li) The ground state
structures for the two fragmentation channels are found by Molecular Dynamics
Simulated Annealing in the framework of Local Density Functional theory.
Energetics considerations suggest that the fragmentation process is dominated
by non-equilibrium processes. We use a real-space approach to solve the
Kohn-Sham problem, where the Laplacian operator is discretized according to the
Mehrstellen scheme, and take advantage of a Full MultiGrid (FMG) strategy to
accelerate convergence. When applied to isolated clusters we find our FMG
method to be more efficient than state-of-the-art plane wave calculations.Comment: 9 pages + 6 Figures (in gzipped tar file
QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo
We review recent advances in the capabilities of the open source ab initio
Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for
greater efficiency and reproducibility. The auxiliary field QMC (AFQMC)
implementation has been greatly expanded to include k-point symmetries,
tensor-hypercontraction, and accelerated graphical processing unit (GPU)
support. These scaling and memory reductions greatly increase the number of
orbitals that can practically be included in AFQMC calculations, increasing
accuracy. Advances in real space methods include techniques for accurate
computation of band gaps and for systematically improving the nodal surface of
ground state wavefunctions. Results of these calculations can be used to
validate application of more approximate electronic structure methods including
GW and density functional based techniques. To provide an improved foundation
for these calculations we utilize a new set of correlation-consistent effective
core potentials (pseudopotentials) that are more accurate than previous sets;
these can also be applied in quantum-chemical and other many-body applications,
not only QMC. These advances increase the efficiency, accuracy, and range of
properties that can be studied in both molecules and materials with QMC and
QMCPACK
Two-level Chebyshev filter based complementary subspace method: pushing the envelope of large-scale electronic structure calculations
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).Comment: Resubmitted version (version 2
Electrons in High-Tc Compounds: Ab-Initio Correlation Results
Electronic correlations in the ground state of an idealized infinite-layer
high-Tc compound are computed using the ab-initio method of local ansatz.
Comparisons are made with the local-density approximation (LDA) results, and
the correlation functions are analyzed in detail. These correlation functions
are used to determine the effective atomic-interaction parameters for model
Hamiltonians. On the resulting model, doping dependencies of the relevant
correlations are investigated. Aside from the expected strong atomic
correlations, particular spin correlations arise. The dominating contribution
is a strong nearest neighbor correlation that is Stoner-enhanced due to the
closeness of the ground state to the magnetic phase. This feature depends
moderately on doping, and is absent in a single-band Hubbard model. Our
calculated spin correlation function is in good qualitative agreement with that
determined from the neutron scattering experiments for a metal.Comment: 21pp, 5fig, Phys. Rev. B (Oct. 98
- …