Applicability of medium-size basis sets in calculations of molecular dynamic polarisabilities

Static and dynamic average polarisabilities and polarisability anisotropies of seven linear non-polar and polar molecules are calculated within the CCS, CC2, and CC3 approximations using a range of medium-sized basis sets: the polarised LPol-n (n = ds, dl, fs, fl), the aug-pc-n (n = 1, 2), the def2-SVPD, and -TZVPD basis sets. Reference values are obtained using a hierarchy of Dunning's (d-)aug-cc-pVXZ (X = D, T, Q, 5) basis sets. The results are discussed together with the available CCSD values in terms of basis set and correlation method errors, and their ratio. Detailed analysis shows that already the def2-SVPD basis set can be used in CCS polarisability calculations. When affordable, the slightly larger aug-pc-1 basis set is recommended, as it leads to significant reduction of basis set error. The def2-TZVPD, LPol-ds, and aug-pc-2 basis sets are optimal choice within the CC2 approximation, with the latter allowing to approach the CC2 basis set limit. The LPol-ds, -dl, and def2-TZVPD sets outperform the aug-cc-pVTZ set in average polarisability CCSD calculations, with the def2-TZVPD being competitive to other reduced-size sets also in determination of polarisability anisotropy. The aug-pc-2 basis is a particularly attractive choice for CCSD, giving the accuracy of aug-cc-pVQZ at a significantly reduced computational cost. The polarisability anisotropy is shown to be more computationally demanding than the average polarisability, in particular with respect to the accuracy of the correlation method and an accurate evaluation of this property requires at least the CCSD model.     

An improved ensemble learning machine for biological activity prediction of tyrosine kinase inhibitors

Boosting is one of the most important strategies in ensemble learning because of its ability to improve the stability and performance of weak learners. It is nonparametric, multivariate, fast and interpretable but is not robust against outliers. To enhance its prediction accuracy as well as immunize it against outliers, a modified version of a boosting algorithm (AdaBoost R2) was developed and called AdaBoost R3. In the sampling step, extremum samples were added to the boosting set. In the robustness step, a modified Huber loss function was applied to overcome the outlier problem. In the output step, a deterministic threshold was used to guarantee that bad predictions do not participate in the final output. The performance of the modified algorithm was investigated with two anticancer data sets of tyrosine kinase inhibitors, and the mechanism of inhibition was studied using the relative weighted variable importance procedure. Investigating the effect of base learner's strength reveals that boosting is only successful using the classification and regression tree method (a weak to moderate learner) and does not have a significant effect using the radial basis functions partial least square method (a strong base learners). Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.     

Search of truncation of (N−1) electron basis containing full connected triple excitations in computing main and satellite ionization potentials via Fock‐space coupled cluster approach

A valence‐universal multireference coupled cluster (VUMRCC) theory, realized via the eigenvalue independent partitioning (EIP) route, has been implemented with full inclusion of triples excitations for computing and analyzing the entire main and several satellite peaks in the ionization potential spectra of several molecules. The EIP‐VUMRCC method, unlike the traditional VUMRCC theory, allows divergence‐free homing‐in to satellite roots which would otherwise have been plagued by intruders, and is thus numerically more robust to obtain more efficient and dependable computational schemes allowing more extensive use of the approach. The computed ionization potentials (IPs) as a result of truncation of the (N−1) electron basis manifold involving virtual functions such as 2h‐p and 3h‐2p by different energy thresholds varying from 5 to 15 a.u. with 1 a.u. intervals as well as thresholds such as 20, 25, and 30 a.u. have been carefully looked into. Cutoff at around 25 a.u. turns out to be an optimal threshold. Molecules such as C₂H₄ and C₂H₂ (X = D,T), and N₂ and CO (X = D,T,Q) with Dunning's cc‐pVXZ bases have been investigated to determine all main and 2h‐p shake‐up and 3h‐2p double shake‐up satellite IPs. We believe that the present work will pave the way to a wider application of the method by providing main and satellite IPs for some problematic N‐electron closed shell systems. © 2013 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.     

一类非线性四阶波动方程初边值问题解的有限时间爆破

研究四阶带有阻尼项的非线性波动方程的解的初边值问题,利用位势井方法,证明了当初值满足一定条件时解发生爆破.将有关该系统爆破性质的研究结果一般化,通过证明得到了该系统较好的性质.     

Spectral analysis and stability of thermoelastic Bresse system with second sound and boundary viscoelastic damping

In this paper, we consider the energy decay rate of a thermoelastic Bresse system with variable coefficients. Assume that the thermo-propagation in the system satisfies the Cattaneo's law, which can eliminate the paradox of infinite speed of thermal propagation in the assumption of the Fourier's law in the classical theory of thermoelasticity. Meanwhile, we also discuss the effect of a boundary viscoelastic damping on the stability of this system. By a detailed spectral analysis, we obtain the expressions of the spectrum and deduce some spectral properties of the system. Then based on the distribution of the spectrum, we prove that the energy of the system with a boundary viscoelastic damping decays exponentially. However, it no longer decays exponentially if there is no boundary damping. Copyright © 2013 John Wiley & Sons, Ltd.