Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMaterial specific optimization of Gaussian basis sets against plane wave data
68
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
Calculations of the hyperpolarizability are typically much more difficult to converge with basis set size than the linear polarizability. In order to understand these convergence issues and hence obtain accurate ab initio values, we compare calculations of the static hyperpolarizability of the gas-phase chloroform molecule (CHCl_3) using three different kinds of basis sets: Gaussian-type orbitals, numerical basis sets, and real-space grids. Although all of these methods can yield similar results, surprisingly large, diffuse basis sets are needed to achieve convergence to comparable values. These results are interpreted in terms of local polarizability and hyperpolarizability densities. We find that the hyperpolarizability is very sensitive to the molecular structure, and we also assess the significance of vibrational contributions and frequency dispersion.
We investigate the use of optimized correlation consistent gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second order M{o}ller--Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely calculations in the Kohn--Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn--Sham basis, and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si and MgO, and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.
We introduce a mixed density fitting scheme that uses both a Gaussian and a plane-wave fitting basis to accurately evaluate electron repulsion integrals in crystalline systems. We use this scheme to enable efficient all-electron Gaussian based periodic density functional and Hartree-Fock calculations.
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.
The traditional Gaussian basis sets used in modern quantum chemistry lack an electron-nuclear cusp, and hence struggle to accurately describe core electron properties. A recently introduced novel type of basis set, mixed ramp-Gaussians, introduce a new primitive function called a ramp function which addresses this issue. This paper introduces three new mixed ramp-Gaussian basis sets - STO-R, STO-RG and STO-R2G, made from a linear combination of ramp and Gaussian primitive functions - which are derived from the single-core-zeta Slater basis sets for the elements Li to Ne. This derivation is done in an analogous fashion to the famous STO-$n$G basis sets. The STO-RG basis functions are found to outperform the STO-3G basis functions and STO-R2G outperforms STO-6G, both in terms of wavefunction fit and other key quantities such as the one-electron energy and the electron-nuclear cusp. The second part of this paper performs preparatory investigations into how standard all-Gaussian basis sets can be converted to ramp-Gaussian basis sets through modifying the core basis functions. Using a test case of the 6-31G basis set for carbon, we determined that the second Gaussian primitive is less important when fitting a ramp-Gaussian core basis function directly to an all-Gaussian core basis function than when fitting to a Slater basis function. Further, we identified the basis sets that are single-core-zeta and thus should be most straightforward to convert to mixed ramp-Gaussian basis sets in the future.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
