Use Of Basis Sets Research Articles

On the utmost importance of the geometry factor of accuracy in the quantum chemical calculations of 31P NMR chemical shifts: New efficient pecG-n (n = 1, 2) basis sets for the geometry optimization procedure.

31P nuclear magnetic resonance (NMR) chemical shifts were shown to be very sensitive to the basis set used at the geometry optimization stage. Commonly used energy-optimized basis sets for a phosphorus atom containing only one polarization d-function were shown to be unable to provide correct equilibrium geometries for the calculations of phosphorus chemical shifts. The use of basis sets with at least two polarization d-functions on a phosphorus atom is strongly recommended. In this paper, an idea of creating the basis sets purposed for the geometry optimization that provide the least possible error coming from the geometry factor of accuracy in the resultant NMR shielding constants is proposed. The property-energy consisted algorithm with the target function in the form of the molecular energy gradient relative to P-P bond lengths was applied to create new geometry-oriented pecG-n (n = 1, 2) basis sets for a phosphorus atom. New basis sets have demonstrated by far superior performance as compared to the other commonly used energy-optimized basis sets in massive calculations of 31P NMR chemical shifts carried out at the gauge-including atomic orbital-coupled cluster singles and doubles/pecS-2 level of the theory by taking into account solvent, vibrational, and relativistic corrections.

Comparison of linear combination modeling strategies for edited magnetic resonance spectroscopy at 3 T

J‐difference‐edited spectroscopy is a valuable approach for the in vivo detection of γ‐aminobutyric‐acid (GABA) with magnetic resonance spectroscopy (MRS). A recent expert consensus article recommends linear combination modeling (LCM) of edited MRS but does not give specific details regarding implementation. This study explores different modeling strategies to adapt LCM for GABA‐edited MRS. Sixty‐one medial parietal lobe GABA‐edited MEGA‐PRESS spectra from a recent 3‐T multisite study were modeled using 102 different strategies combining six different approaches to account for co‐edited macromolecules (MMs), three modeling ranges, three baseline knot spacings, and the use of basis sets with or without homocarnosine. The resulting GABA and GABA+ estimates (quantified relative to total creatine), the residuals at different ranges, standard deviations and coefficients of variation (CVs), and Akaike information criteria, were used to evaluate the models' performance. Significantly different GABA+ and GABA estimates were found when a well‐parameterized MM3cobasis function was included in the model. The mean GABA estimates were significantly lower when modeling MM3co, while the CVs were similar. A sparser spline knot spacing led to lower variation in the GABA and GABA+ estimates, and a narrower modeling range—only including the signals of interest—did not substantially improve or degrade modeling performance. Additionally, the results suggest that LCM can separate GABA and the underlying co‐edited MM3co. Incorporating homocarnosine into the modeling did not significantly improve variance in GABA+ estimates. In conclusion, GABA‐edited MRS is most appropriately quantified by LCM with a well‐parameterized co‐edited MM3cobasis function with a constraint to the nonoverlapped MM0.93, in combination with a sparse spline knot spacing (0.55 ppm) and a modeling range of 0.5–4 ppm.

Charge Conjugation Symmetry in the Finite Basis Approximation of the Dirac Equation

Four-component relativistic atomic and molecular calculations are typically performed within the no-pair approximation where negative-energy solutions are discarded. These states are, however, needed in QED calculations, wherein, furthermore, charge conjugation symmetry, which connects electronic and positronic solutions, becomes an issue. In this work, we shall discuss the realization of charge conjugation symmetry of the Dirac equation in a central field within the finite basis approximation. Three schemes for basis set construction are considered: restricted, inverse, and dual kinetic balance. We find that charge conjugation symmetry can be realized within the restricted and inverse kinetic balance prescriptions, but only with a special form of basis functions that does not obey the right boundary conditions of the radial wavefunctions. The dual kinetic balance prescription is, on the other hand, compatible with charge conjugation symmetry without restricting the form of the radial basis functions. However, since charge conjugation relates solutions of opposite value of the quantum number κ , this requires the use of basis sets chosen according to total angular momentum j rather than orbital angular momentum ℓ. As a special case, we consider the free-particle Dirac equation, where opposite energy solutions are related by charge conjugation symmetry. We show that there is additional symmetry in that solutions of the same value of κ come in pairs of opposite energy.

Intramolecular BSSE and dispersion affect the structure of a dipeptide conformer

ABSTRACTB3LYP and MP2 calculations with the commonly-used 6-31+G(d) basis set predict qualitatively different structures for the Tyr–Gly conformer book1, which is the most stable conformer identified in a previous study. The structures differ mainly in the ψtyr Ramachandran angle (138° in the B3LYP structure and 120° in the MP2 structure). The causes for the discrepant structures are attributed to missing dispersion in the B3LYP calculations and large intramolecular BSSE in the MP2 calculations. The correct ψtyr value is estimated to be 130°. The MP2/6-31+G(d) profile identified an additional conformer, not present on the B3LYP surface, with a ψtyr value of 96° and a more folded structure. This minimum is, however, likely an artefact of large intramolecular BSSE values. We recommend the use of basis sets of at least quadruple-zeta quality in density functional theory (DFT), DFTaugmented with an empirical dispersion term (DFT-D) and second-order Møller-Plesset perturbation theory (MP2 ) calculations in cases where intramolecular BSSE is expected to be large.

Gaussian basis sets for use in correlated molecular calculations. XI. Pseudopotential-based and all-electron relativistic basis sets for alkali metal (K-Fr) and alkaline earth (Ca-Ra) elements.

New correlation consistent basis sets based on pseudopotential (PP) Hamiltonians have been developed from double- to quintuple-zeta quality for the late alkali (K-Fr) and alkaline earth (Ca-Ra) metals. These are accompanied by new all-electron basis sets of double- to quadruple-zeta quality that have been contracted for use with both Douglas-Kroll-Hess (DKH) and eXact 2-Component (X2C) scalar relativistic Hamiltonians. Sets for valence correlation (ms), cc-pVnZ-PP and cc-pVnZ-(DK,DK3/X2C), in addition to outer-core correlation [valence + (m-1)sp], cc-p(w)CVnZ-PP and cc-pwCVnZ-(DK,DK3/X2C), are reported. The -PP sets have been developed for use with small-core PPs [I. S. Lim et al., J. Chem. Phys. 122, 104103 (2005) and I. S. Lim et al., J. Chem. Phys. 124, 034107 (2006)], while the all-electron sets utilized second-order DKH Hamiltonians for 4s and 5s elements and third-order DKH for 6s and 7s. The accuracy of the basis sets is assessed through benchmark calculations at the coupled-cluster level of theory for both atomic and molecular properties. Not surprisingly, it is found that outer-core correlation is vital for accurate calculation of the thermodynamic and spectroscopic properties of diatomic molecules containing these elements.

Approaching the Hartree-Fock Limit through the Complementary Auxiliary Basis Set Singles Correction and Auxiliary Basis Sets.

Auxiliary basis sets for use in the resolution of the identity (RI) approximation in explicitly correlated methods are presented for the elements H-Ar. These extend the cc-pVnZ-F12/OptRI (n = D-Q) auxiliary basis sets of Peterson and co-workers by the addition of a small number of s- and p-functions, optimized so as to yield the greatest complementary auxiliary basis set (CABS) singles correction to the Hartree-Fock energy. The new sets, denoted OptRI+, also lead to a reduction in errors due to the RI approximation and hence an improvement in correlation energies. The atomization energies and heats of formation for a test set of small molecules, and spectroscopic constants for 27 diatomics, calculated at the CCSD(T)-F12b level, are shown to have improved error distributions for the new auxiliary basis sets with negligible additional effort. The OptRI+ sets retain all of the desirable properties of the original OptRI, including the production of smooth potential energy surfaces, while maintaining a compact nature.

Multipole moments for embedding potentials: Exploring different atomic allocation algorithms

Polarizable quantum mechanical (QM)/molecular mechanics (MM)-embedding methods are currently among the most promising methods for computationally feasible, yet reliable, production calculations of localized excitations and molecular response properties of large molecular complexes, such as proteins and RNA/DNA, and of molecules in solution. Our aim is to develop a computational methodology for distributed multipole moments and their associated multipole polarizabilities which is accurate, computationally efficient, and with smooth convergence with respect to multipole order. As the first step toward this goal, we herein investigate different ways of obtaining distributed atom-centered multipole moments that are used in the construction of the electrostatic part of the embedding potential. Our objective is methods that not only are accurate and computationally efficient, but which can be consistently extended with site polarizabilities including internal charge transfer terms. We present a new way of dealing with well-known problems in relation to the use of basis sets with diffuse functions in conventional atomic allocation algorithms, avoiding numerical integration schemes. Using this approach, we show that the classical embedding potential can be systematically improved, also when using basis sets with diffuse functions, and that very accurate embedding potentials suitable for QM/MM embedding calculations can be acquired. © 2016 Wiley Periodicals, Inc.

Optical response of silver clusters and their hollow shells from linear-response TDDFT

We present a study of the optical response of compact and hollow icosahedral clusters containing up to 868 silver atoms by means of time-dependent density functional theory. We have studied the dependence on size and morphology of both the sharp plasmonic resonance at 3–4 eV (originated mainly from sp-electrons), and the less studied broader feature appearing in the 6–7 eV range (interband transitions). An analysis of the effect of structural relaxations, as well as the choice of exchange correlation functional (local density versus generalised gradient approximations) both in the ground state and optical response calculations is also presented. We have further analysed the role of the different atom layers (surface versus inner layers) and the different orbital symmetries on the absorption cross-section for energies up to 8 eV. We have also studied the dependence on the number of atom layers in hollow structures. Shells formed by a single layer of atoms show a pronounced red shift of the main plasmon resonances that, however, rapidly converge to those of the compact structures as the number of layers is increased. The methods used to obtain these results are also carefully discussed. Our methodology is based on the use of localised basis (atomic orbitals, and atom-centered and dominant-product functions), which bring several computational advantages related to their relatively small size and the sparsity of the resulting matrices. Furthermore, the use of basis sets of atomic orbitals also allows the possibility of extending some of the standard population analysis tools (e.g. Mulliken population analysis) to the realm of optical excitations. Some examples of these analyses are described in the present work.

Optimized Basis Sets for the Environment in the Domain-Specific Basis Set Approach of the Incremental Scheme.

Minimal basis sets, denoted DSBSenv, based on the segmented basis sets of Ahlrichs and co-workers have been developed for use as environmental basis sets for the domain-specific basis set (DSBS) incremental scheme with the aim of decreasing the CPU requirements of the incremental scheme. The use of these minimal basis sets within explicitly correlated (F12) methods has been enabled by the optimization of matching auxiliary basis sets for use in density fitting of two-electron integrals and resolution of the identity. The accuracy of these auxiliary sets has been validated by calculations on a test set containing small- to medium-sized molecules. The errors due to density fitting are about 2-4 orders of magnitude smaller than the basis set incompleteness error of the DSBSenv orbital basis sets. Additional reductions in computational cost have been tested with the reduced DSBSenv basis sets, in which the highest angular momentum functions of the DSBSenv auxiliary basis sets have been removed. The optimized and reduced basis sets are used in the framework of the domain-specific basis set of the incremental scheme to decrease the computation time without significant loss of accuracy. The computation times and accuracy of the previously used environmental basis and that optimized in this work have been validated with a test set of medium- to large-sized systems. The optimized and reduced DSBSenv basis sets decrease the CPU time by about 15.4% and 19.4% compared with the old environmental basis and retain the accuracy in the absolute energy with standard deviations of 0.99 and 1.06 kJ/mol, respectively.

Implications of the complete basis set limit in valence bond theory: a case study of molecular hydrogen

The use of larger basis sets to approach the complete basis limit, now common in quantum chemistry, is applied for the first time to a range of valence bond functions for the simplest case of molecular hydrogen. Good convergence of the energy is slow due to difficulty in getting a correct cusp near the nuclei, but it is significant. The form of the orbitals converges much faster, leading to a slight distortion of the valence bond orbitals and an enhanced overlap, irrespective whether the basis set is restricted to basis functions centred on one atom for each valence bond orbital or the full use of the basis set is allowed. This blurs the distinction between these two approaches and shows that basis set restrictions are not tenable in the complete basis set limit. Furthermore, it supports the general use of the full basis as advocated in the spin-coupled and generalised valence bond methods.

Improved Segmented All-Electron Relativistically Contracted Basis Sets for the Lanthanides.

Improved versions of the segmented all-electron relativistically contracted (SARC) basis sets for the lanthanides are presented. The second-generation SARC2 basis sets maintain efficient construction of their predecessors and their individual adaptation to the DKH2 and ZORA Hamiltonians, but feature exponents optimized with a completely new orbital shape fitting procedure and a slightly expanded f space that results in sizable improvement in CASSCF energies and in significantly more accurate prediction of spin-orbit coupling parameters. Additionally, an extended set of polarization/correlation functions is constructed that is appropriate for multireference correlated calculations and new auxiliary basis sets for use in resolution-of-identity (density-fitting) approximations in combination with both DFT and wave function based treatments. Thus, the SARC2 basis sets extend the applicability of the first-generation DFT-oriented basis sets to routine all-electron wave function-based treatments of lanthanide complexes. The new basis sets are benchmarked with respect to excitation energies, radial distribution functions, optimized geometries, orbital eigenvalues, ionization potentials, and spin-orbit coupling parameters of lanthanide systems and are shown to be suitable for the description of magnetic and spectroscopic properties using both DFT and multireference wave function-based methods.

Geminal embedding scheme for optimal atomic basis set construction in correlated calculations

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

Auxiliary Basis Sets for Density Fitting in Explicitly Correlated Calculations: The Atoms H-Ar.

Auxiliary basis sets specifically matched to the correlation consistent cc-pVnZ-F12 and cc-pCVnZ-F12 orbital basis sets for the elements H-Ar have been optimized at the density-fitted second-order Møller-Plesset perturbation theory level of theory for use in explicitly correlated (F12) methods, which utilize density fitting for the evaluation of two-electron integrals. Calculations of the correlation energy for a test set of small to medium sized molecules indicate that the density fitting error when using these auxiliary sets is 2 to 3 orders of magnitude smaller than the F12 orbital basis set incompleteness error. The error introduced by the use of these fitting sets within the resolution-of-the-identity approximation of the many-electron integrals arising in F12 theory has also been assessed and is demonstrated to be negligible and well-controlled. General guidelines are proposed for the optimization of density fitting auxiliary basis sets for use with F12 methods for other elements.

Correlation Consistent Gaussian Basis Sets for H, B-Ne with Dirac-Fock AREP Pseudopotentials: Applications in Quantum Monte Carlo Calculations.

In this paper, we introduce correlation consistent Gaussian-type orbital basis sets for the H and B-Ne atoms for use with the CASINO Dirac-Fock AREP pseudopotentials. These basis sets are tested in coupled cluster calculations on H2, B2, C2, N2, O2, and F2 as well as in quantum Monte Carlo calculations on the water monomer and dimer and the water-benzene complex, where they are found to give low variances in variational Monte Carlo calculations and to lead to reduced time step errors and improved convergence in diffusion Monte Carlo calculations compared to the use of nonoptimized basis sets. The use of basis sets with a large number of contracted s and p primitives is found to be especially important for the convergence of the energy in the diffusion Monte Carlo calculations.

Generally contracted valence–core/valence basis sets for use with relativistic effective core potentials and spin–orbit coupling operators

A procedure for structuring generally contracted valence–core/valence basis sets of Gaussian-type functions for use with relativistic effective core potentials (gcv-c/v-RECP basis sets) is presented. Large valence basis sets are enhanced using a compact basis set derived for outer core electrons in the presence of small-core RECPs. When core electrons are represented by RECPs, and appropriate levels of theory, these basis sets are shown to provide accurate representations of atomic and molecular valence and outer-core electrons. Core/valence polarization and correlation effects can be calculated using these basis sets through standard methods for treating electron correlation. Calculations of energies and spectra for Ru, Os, Ir, In and Cs are reported. Spectroscopic constants for RuO2+, OsO2+, Cs2 and InH are calculated and compared with experiment.

Methods of encrypting monotonic access structures

We will present some ideas about sharing a secret in a monotonic access structure. We show the relations which occur between the method of encrypting a monotonic access structure with the use of basis sets or maximal unprivileged sets, and that based on a logical formula (used by Benaloh and Leichter in [1]). We will also give some facts connected with the problem of security, including the aspects of a hierarchy security in the structure. The method of encrypting a monotonic access structure using a family of basis sets or a family of maximal sets that cannot reconstruct the secret will be described in a general way. Some aspects of using the latter based on a logical formula will be also given. Any (general) access structure can be encrypted by each of them but the way of sharing a secret is quite different and usually a specified method has to be chosen to achieve a desirable level of security and time complexity.

Gaussian basis sets for use in correlated molecular calculations. VII. Valence, core-valence, and scalar relativistic basis sets for Li, Be, Na, and Mg

Correlation consistent basis sets of double-ζ through quintuple-ζ quality for the alkali and alkaline earth metals Li, Be, Na, and Mg have been developed, including the valence (cc-pVnZ), augmented valence (aug-cc-pVnZ), core-valence (cc-pCVnZ), and weighted core-valence (cc-pwCVnZ) basis sets. The basis sets are also re-contracted for Douglas–Kroll scalar relativistic calculations and are found to be superior to non-relativistic basis sets in recovering scalar relativistic effects. CCSD(T) computations have been performed with these basis sets, and a series of properties have been examined, including atomic ionization potentials and electron affinities, optimized molecular geometries, harmonic vibrational frequencies, atomization energies, and enthalpies of formation for the molecules Li2, LiF, BeO, BeF, BeH2, BeF2, Na2, NaF, MgO, MgF, MgH2, and MgF2.

Correlation consistent basis sets for molecular core-valence effects with explicitly correlated wave functions: The atoms B–Ne and Al–Ar

Correlation consistent basis sets have been optimized for accurately describing core-core and core-valence correlation effects with explicitly correlated F12 methods. The new sets, denoted cc-pCVnZ-F12 (n=D, T, Q) and aug-cc-pC(F12)VnZ (n=D, T, Q, 5), were developed by augmenting the cc-pVnZ-F12 and aug-cc-pVnZ families of basis sets with additional functions whose exponents were optimized based on the difference between all-electron and valence-electron correlation energies. The number of augmented functions added is fewer, in general, than in the standard cc-pCVnZ and cc-pwCVnZ families of basis sets. Optimal values of the geminal Slater exponent for use with these basis sets in MP2-F12 calculations are presented and are also recommended for CCSD-F12b calculations. Auxiliary basis sets for use in the resolution of the identity approximation in explicitly correlated calculations have also been optimized and matched to the new cc-pCVnZ-F12 series of orbital basis sets. The cc-pCVnZ-F12 basis sets, along with the new auxiliary sets, were benchmarked in CCSD(T)-F12b calculations of spectroscopic properties on a series of homo- and heteronuclear first and second row diatomic molecules. Comparing the effects of correlating the outer core electrons in these molecules with those from conventional CCSD(T) at the complete basis set limit, which involved calculations with new cc-pCV6Z basis sets for the second row elements that were also developed in the course of this work, it is observed that the F12 values are reasonably well converged already at just the triple-zeta level.

Correlation consistent basis sets for explicitly correlated wavefunctions: valence and core–valence basis sets for Li, Be, Na, and Mg

Correlation consistent basis sets for the alkali and alkaline earth metals Li, Be, Na, and Mg have been optimized for use with explicitly correlated F12 methods. These include orbital basis sets for valence-only (denoted cc-pVnZ-F12) and core-valence (cc-pCVnZ-F12) correlation, along with compact complementary auxiliary basis sets for use in the resolution of the identity approximation. Additional auxiliary basis sets that can be employed in the density fitting of two-electron integrals in both explicitly correlated methods and in more conventional correlated methods, such as density fitted second-order Møller-Plesset perturbation theory (DF-MP2), have also been developed by adding additional diffuse or core-valence functions to the cc-pVnZ/MP2FIT sets of Hättig. Explicitly correlated calculations with the approximate CCSD(T)-F12b method have been carried out with these basis sets on a series of sixteen test molecules to demonstrate their accuracy for optimized geometries, harmonic vibrational frequencies, and atomization energies. Results comparable to conventional CCSD(T) quintuple-zeta, which are near the complete basis set limits for these molecules, are obtained using CCSD(T)-F12b with just triple-zeta quality basis sets. The effects on the spectroscopic constants from correlating the outer core electrons are accurately recovered with just the cc-pCVDZ-F12 basis sets.

Development of basis sets for molecular calculations. Comparison of tetrahedral gaussian lobe functions and cartesian gaussian basis sets

In order to allow development of a range of basis sets for use of floating spherical Gaussian orbitals in molecular calculations, an analysis of numerical errors that result from using off-center functions was carried out. This analysis included both rotational variance and numerical differencing problems that can arise when lobe functions are employed. The analysis and numerical studies indicated that satisfactory accuracy in integral evaluation can be maintained if the ratio of distance from the origin to orbital radius is taken to be in the range 0.002 to 0.005 for computers having 72-bit words. 2 figures, 1 table.