Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMixed ramp-Gaussian basis sets for core-dependent properties: STO-RG and STO-R2G for Li-Ne
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
An extensive analysis has been carried out of the performance of standard families of basis sets with the hierarchy of coupled cluster methods CC2, CCSD, CC3 and CCSDT in computing selected Oxygen, Carbon and Nitrogen K-edge (vertical) core excitation and ionization energies within a core-valence separated scheme in the molecules water, ammonia, and carbon monoxide. Complete basis set limits for the excitation energies have been estimated via different basis set extrapolation schemes. The importance of scalar relativistic effects has been established within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e).
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.
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.
Heterostructures of transition metal oxides (TMO) perovskites represent an ideal platform to explore exotic phenomena involving the complex interplay between the spin, charge, orbital and lattice degrees of freedom available in these compounds. At the interface between such materials, this interplay can lead to phenomena that are present in none of the original constituents such as the formation of the interfacial 2D electron system (2DES) discovered at the LAO3/STO3 (LAO/STO) interface. In samples prepared by growing a LAO layer onto a STO substrate, the 2DES is confined in a band bending potential well, whose width is set by the interface charge density and the STO dielectric properties, and determines the electronic band structure. Growing LAO (2 nm) /STO (x nm)/LAO (2 nm) heterostructures on STO substrates allows us to control the extension of the confining potential of the top 2DES via the thickness of the STO layer. In such samples, we explore the dependence of the electronic structure on the width of the confining potential using soft X-ray ARPES combined with ab-initio calculations. The results indicate that varying the thickness of the STO film modifies the quantization of the 3d t2g bands and, interestingly, redistributes the charge between the dxy and dxz/dyz bands.
We derive distance-dependent estimators for two-center and three-center electron repulsion integrals over a short-range Coulomb potential, $textrm{erfc}(omega r_{12})/r_{12}$. These estimators are much tighter than one based on the Schwarz inequality and can be viewed as a complement to the distance-dependent estimators for four-center short-range Coulomb integrals and for two-center and three-center full Coulomb integrals previously reported. Because the short-range Coulomb potential is commonly used in solid-state calculations, including those with the HSE functional and with our recently introduced range-separated periodic Gaussian density fitting, we test our estimators on a diverse set of periodic systems using a wide range of the range-separation parameter $omega$. These tests demonstrate the robust tightness of our estimators, which are then used with integral screening to calculate periodic three-center short-range Coulomb integrals with linear scaling in system size.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
