Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence of many-body wavefunction expansions using a plane wave basis in the thermodynamic limit
105
0
0.0
(
0
)
Authors
James J. Shepherd
Ask ChatGPT about the research
No Arabic abstract
Basis set incompleteness error and finite size error can manifest concurrently in systems for which the two effects are phenomenologically well-separated in length scale. When this is true, we need not necessarily remove the two sources of error simultaneously. Instead, the errors can be found and remedied in different parts of the basis set. This would be of great benefit to a method such as coupled cluster theory since the combined cost of $n_{occ}^6 n_{virt}^4$ could be separated into $n_{occ}^6$ and $n_{virt}^4$ costs with smaller prefactors. In this Communication, we present analysis on a data set due to Baardsen and coworkers, containing coupled cluster doubles energies for the 2DEG for $r_s=$ 0.5, 1.0 and 2.0 a.u.~at a wide range of basis set sizes and particle numbers. In obtaining complete basis set limit thermodynamic limit results, we find that within a small and removable error the above assertion is correct for this simple system. This approach allows for the combination of methods which separately address finite size effects and basis set incompleteness error.
rate research
Read More
We describe a method for computing near-exact energies for correlated systems with large Hilbert spaces. The method efficiently identifies the most important basis states (Slater determinants) and performs a variational calculation in the subspace spanned by these determinants. A semistochastic approach is then used to add a perturbative correction to the variational energy to compute the total energy. The size of the variational space is progressively increased until the total energy converges to within the desired tolerance. We demonstrate the power of the method by computing a near-exact potential energy curve (PEC) for a very challenging molecule -- the chromium dimer.
The present work proposes to use density-functional theory (DFT) to correct for the basis-set error of wave-function theory (WFT). One of the key ideas developed here is to define a range-separation parameter which automatically adapts to a given basis set. The derivation of the exact equations are based on the Levy-Lieb formulation of DFT, which helps us to define a complementary functional which corrects uniquely for the basis-set error of WFT. The coupling of DFT and WFT is done through the definition of a real-space representation of the electron-electron Coulomb operator projected in a one-particle basis set. Such an effective interaction has the particularity to coincide with the exact electron-electron interaction in the limit of a complete basis set, and to be finite at the electron-electron coalescence point when the basis set is incomplete. The non-diverging character of the effective interaction allows one to define a mapping with the long-range interaction used in the context of range-separated DFT and to design practical approximations for the unknown complementary functional. Here, a local-density approximation is proposed for both full-configuration-interaction (FCI) and selected configuration-interaction approaches. Our theory is numerically tested to compute total energies and ionization potentials for a series of atomic systems. The results clearly show that the DFT correction drastically improves the basis-set convergence of both the total energies and the energy differences. For instance, a sub kcal/mol accuracy is obtained from the aug-cc-pVTZ basis set with the method proposed here when an aug-cc-pV5Z basis set barely reaches such a level of accuracy at the near FCI level.
We present the N-body simulation techniques in EXP. EXP uses empirically-chosen basis functions to expand the potential field of an ensemble of particles. Unlike other basis function expansions, the derived basis functions are adapted to an input mass distribution, enabling accurate expansion of highly non-spherical objects, such as galactic discs. We measure the force accuracy in three models, one based on a spherical or aspherical halo, one based on an exponential disc, and one based on a bar-based disc model. We find that EXP is as accurate as a direct-summation or tree-based calculation, and in some ways is better, while being considerably less computationally intensive. We discuss optimising the computation of the basis function representation. We also detail numerical improvements for performing orbit integrations, including timesteps.
The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded electron densities, the sum of which is the electron density of the system of interest. Coupling between the substrate and embedded systems is achieved via approximate kinetic energy functionals. Bulk aluminium is examined as a test case for which there is a strong interaction between the substrate and embedded systems. A number of approximations to the kinetic-energy functional, both semi-local and non-local, are investigated. It is found that Kohn-Sham results can be well reproduced using a non-local kinetic energy functional, with the total energy accurate to better than 0.1 eV per atom and good agreement between the electron densities.
Since in periodic systems, a given element may be present in different spatial arrangements displaying vastly different physical and chemical properties, an elemental basis set that is independent of physical properties of materials may lead to significant simulation inaccuracies. To avoid such a lack of material specificity within a given basis set, we present a material-specific Gaussian basis optimization scheme for solids, which simultaneously minimizes the total energy of the system and optimizes the band energies when compared to the reference plane wave calculation while taking care of the overlap matrix condition number. To assess this basis set optimization scheme, we compare the quality of the Gaussian basis sets generated for diamond, graphite, and silicon via our method against the existing basis sets. The optimization scheme of this work has also been tested on the existing Gaussian basis sets for periodic systems such as MoS$_2$ and NiO yielding improved results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
