Improved efficiency of focal point conformational analysis with truncated correlation consistent basis sets

It has been suggested that the computational cost of correlated ab initio calculations could be reduced efficiently by using truncated basis sets on hydrogen atoms (Mintz et al., J Chem Phys 2004, 121, 5629). We now explore this proposal in the context of conformational analysis of small molecules, such as hydrogen peroxide, dimethyl ether, ethyl methyl ether, formic acid, methyl formate, and several small alcohols. It is found that truncated correlation consistent basis sets that lack certain higher angular momentum functions on hydrogen atoms offer accuracy similar to traditional Dunning's basis sets for conformational analysis. Combination of such basis sets with the basis set extrapolation technique to estimate Hartree-Fock and M?ller-Plesset second order energies provides composite extrapolation model chemistries that are significantly more accurate and faster than analogous single point calculations with traditional correlation consistent basis sets. Root mean square errors of best composite extrapolation model chemistries on the used set of molecules are within 0.03 kcal/mol of traditional focal point conformational energies. The applicability of composite extrapolation methods is illustrated by performing conformational analysis of tert-butanol and cyclohexanol. For comparison, conformational energies calculated with popular molecular mechanics force fields are also given.     

Accurate Gaussian basis sets for the H2O molecule generated with the molecular improved generator coordinate Hartree–Fock method

The molecular improved generator coordinate Hartree–Fock (MIGCHF) method is used to generate accurate basis sets of primitive Gaussian-type functions for the H₂O molecule. Sequences of increasing size atom centered basis sets are employed to explore the accuracy that can be achieved with this method. Using the O(24s14p8d5f2g1h);H(22s9p5d2f1g) basis set, the HF and second-order electron correlation energies of the H₂O ground state at the experimental geometry are computed as −76.0674680 and −0.3491935 hartree, respectively. The HF energy is in error by 20 μhartree and the second-order correlation energy corresponds to 96.5% of an estimate of the limiting value. The relevance of the present calculations is to show the accuracy that can be achieved in studies of small polyatomic molecules with the MIGCHF method.     

Polarization functions for the modified m6‐31G basis sets for atoms Ga through Kr

The 2df polarization functions for the modified m6‐31G basis sets of the third‐row atoms Ga through Kr (Int J Quantum Chem, 2007, 107, 3028; Int J. Quantum Chem, 2009, 109, 1158) are proposed. The performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets were examined in molecular calculations carried out by the density functional theory (DFT) method with B3LYP hybrid functional, Møller‐Plesset perturbation theory of the second order (MP2), quadratic configuration interaction method with single and double substitutions and were compared with those for the known 6‐31G basis sets as well as with the other similar 641 and 6‐311G basis sets with and without polarization functions. Obtained results have shown that the performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets are better in comparison with the performances of the known 6‐31G, 6‐31G(d,p) and 6‐31G(2df,p) basis sets. These improvements are mainly reached due to better approximations of different electrons belonging to the different atomic shells in the modified basis sets. Applicability of the modified basis sets in thermochemical calculations is also discussed. © 2013 Wiley Periodicals, Inc.     

Gaussian basis sets for molecular applications

The choice of basis set in quantum chemical calculations can have a huge impact on the quality of the results, especially for correlated ab initio methods. This article provides an overview of the development of Gaussian basis sets for molecular calculations, with a focus on four popular families of modern atom‐centered, energy‐optimized bases: atomic natural orbital, correlation consistent, polarization consistent, and def2. The terminology used for describing basis sets is briefly covered, along with an overview of the auxiliary basis sets used in a number of integral approximation techniques and an outlook on possible future directions of basis set design. © 2012 Wiley Periodicals, Inc.     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

Efficient polarized basis sets for evaluation of static and dynamic molecular electric properties

New medium size Gaussian‐type basis set R‐ORP for evaluation of static and dynamic electric properties in molecular systems is presented. It is obtained in a close resemblance to the original ORP basis set, from the source basis set through addition of two first‐order polarization functions whose exponent values are optimized with respect to the finite field restricted open‐shell Hartree–Fock (ROHF) atomic polarizabilities. As the source set the VTZ basis set of Ahlrichs and coworkers, augmented with additional diffuse functions and contracted to the form [6s/3s] for hydrogen and [11s7p/4s3p] for carbon through fluorine, is chosen. The resulting basis set is of the form [6s2p/3s2p] for hydrogen and [11s7p2d/4s3p2d] for other atoms. Presented basis set is next tested in the CCSD static and dynamic molecular polarizability and hyperpolarizability calculations for a set of ten and four test molecules, respectively, for which very accurate reference data exist. Additionally, the recently developed ORP basis set is employed in the calculations to examine the limits of its applicability. Results are compared to the literature data obtained in both, large and diffuse, as well as reduced‐size basis sets. In the case of polarizability calculations, the aug‐pc‐1 and R‐ORP are the optimal choices among the investigated smaller basis sets, with the overall performance of the aug‐pc‐1 set being better. Among the larger sets, the ORP performs better in the case of average polarizability, while the RMSE values for polarizability anisotropy are practically identical for d‐aug‐cc‐pVDZ and ORP sets. Finally, the R‐ORP and ORP basis sets compete other small bases in the evaluation of the first hyperpolarizability in investigated systems. © 2016 Wiley Periodicals, Inc.     

Oxidative addition of the ethane C-C bond to Pd. An ab initio benchmark and DFT validation study

We have computed a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the ethane C-C bond to the palladium atom and have used this to evaluate the performance of 24 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing this reaction. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods [HF, MP2, CCSD, CCSD(T)] in combination with a hierarchical series of five Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -10.8 (-11.3) kcal/mol for the formation of the reactant complex, 19.4 (17.1) kcal/mol for the activation energy relative to the separate reactants, and -4.5 (-6.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.5 to 2.5 kcal/mol and errors in activation energies ranging from -0.2 to -3.2 kcal/mol. Interestingly, the well-known BLYP functional compares very reasonably with a slight underestimation of the overall barrier by -0.9 kcal/mol. For comparison, with B3LYP we arrive at an overestimation of the overall barrier by 5.8 kcal/mol. On the other hand, B3LYP performs excellently for the central barrier (i.e., relative to the reactant complex) which it underestimates by only -0.1 kcal/mol.     

Irregular 1-sets on the graphs of continuous functions

One can define in a natural way irregular 1-sets on the graphs of several fractal functions, like Takagi’s function, Weierstrass-Cellerier type functions and the typical continuous function. These irregular 1-sets can be useful during the investigation of level-sets and occupation measures of these functions. For example, we see that for Takagi’s function and for certain Weierstrass-Cellerier functions the occupation measure is singular with respect to the Lebesgue measure and for almost every level the level set is finite.     

Rejoinder on “Statistical reasoning with set-valued information: Ontic vs. epistemic views”

This note replies to comments made on our contribution to the Low Quality Data debate.     

Large and Moderate Deviations of Random Upper Semicontinuous Functions

We firstly discuss the topological properties of the space of upper semicontinuous functions, and then we obtain large deviation principles for random upper semicontinuous functions under various topologies. Finally, we prove moderate deviation principles for random sets and random upper semicontinuous functions.