MP2-F12 Basis Set Convergence near the Complete Basis Set Limit: Are h Functions Sufficient?
Ontology highlight
ABSTRACT:
SUBMITTER: Mehta N
PROVIDER: S-EPMC9234959 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S0S1) and the multiexcitonic state (1TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc-pVnZ series of basis sets, that both the aug-cc-pVDZ and 6-311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.
| S-EPMC7898305 | biostudies-literature
Project description:Benchmarking molecular properties with Gaussian-type orbital (GTO) basis sets can be challenging, because one has to assume that the computed property is at the complete basis set (CBS) limit, without a robust measure of the error. Multiwavelet (MW) bases can be systematically improved with a controllable error, which eliminates the need for such assumptions. In this work, we have used MWs within Kohn-Sham density functional theory to compute static polarizabilities for a set of 92 closed-shell and 32 open-shell species. The results are compared to recent benchmark calculations employing the GTO-type aug-pc4 basis set. We observe discrepancies between GTO and MW results for several species, with open-shell systems showing the largest deviations. Based on linear response calculations, we show that these discrepancies originate from artifacts caused by the field strength and that several polarizabilies from a previous study were contaminated by higher order responses (hyperpolarizabilities). Based on our MW benchmark results, we can affirm that aug-pc4 is able to provide results close to the CBS limit, as long as finite difference effects can be controlled. However, we suggest that a better approach is to use MWs, which are able to yield precise finite difference polarizabilities even with small field strengths.
| S-EPMC7467643 | biostudies-literature
Project description:In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested.The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies. Graphical Abstract In this study, the performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, and 6) was tested. For the first time, we demonstrated regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit.
| S-EPMC4723629 | biostudies-literature
Project description:A systematic investigation of three different electron-electron entanglement measures, namely the von Neumann, the linear and the occupation number entropy at full configuration interaction level has been performed for the four helium-like systems hydride, helium, Li(+) and Be(2+) using a large number of different basis sets. The convergence behavior of the resulting energies and entropies revealed that the latter do in general not show the expected strictly monotonic increase upon increase of the one-electron basis. Overall, the three different entanglement measures show good agreement among each other, the largest deviations being observed for small basis sets. The data clearly demonstrates that it is important to consider the nature of the chemical system when investigating entanglement phenomena in the framework of Gaussian type basis sets: while in case of hydride the use of augmentation functions is crucial, the application of core functions greatly improves the accuracy in case of cationic systems such as Li(+) and Be(2+). In addition, numerical derivatives of the entanglement measures with respect to the nucleic charge have been determined, which proved to be a very sensitive probe of the convergence leading to qualitatively wrong results (i.e., the wrong sign) if too small basis sets are used.
| S-EPMC3982574 | biostudies-literature
Project description:We compare dispersion and induction interactions for noble gas dimers and for Ne, methane, and 2-butyne with HF and LiF using a variety of functionals (including some specifically parameterized to evaluate dispersion interactions) with ab initio methods including CCSD(T) and MP2. We see that inductive interactions tend to enhance dispersion and may be accompanied by charge-transfer. We show that the functionals do not generally follow the expected trends in interaction energies, basis set superposition errors (BSSE), and interaction distances as a function of basis set size. The functionals parameterized to treat dispersion interactions often overestimate these interactions, sometimes by quite a lot, when compared to higher level calculations. Which functionals work best depends upon the examples chosen. The B3LYP and X3LYP functionals, which do not describe pure dispersion interactions, appear to describe dispersion mixed with induction about as accurately as those parametrized to treat dispersion. We observed significant differences in high-level wavefunction calculations in a basis set larger than those used to generate the structures in many of the databases. We discuss the implications for highly parameterized functionals based on these databases, as well as the use of simple potential energy for fitting the parameters rather than experimentally determinable thermodynamic state functions that involve consideration of vibrational states.
| S-EPMC3477180 | biostudies-other
Project description:The calculation of electron correlation is vital for the description of atomistic phenomena in physics, chemistry, and biology. However, accurate wavefunction-based methods exhibit steep scaling and often sluggish convergence with respect to the basis set at hand. Because of their delocalization and ease of extrapolation to the basis-set limit, plane waves would be ideally suited for the calculation of basis-set limit correlation energies. However, the routine use of correlated wavefunction approaches in a plane-wave basis set is hampered by prohibitive scaling due to a large number of virtual continuum states and has not been feasible for all but the smallest systems, even if substantial computational resources are available and methods with comparably beneficial scaling, such as the Møller-Plesset perturbation theory to second order (MP2), are used. Here, we introduce a stochastic sampling of the MP2 integrand based on Monte Carlo summation over continuum orbitals, which allows for speedups of up to a factor of 1000. Given a fixed number of sampling points, the resulting algorithm is dominated by a flat scaling of ∼O(N2). Absolute correlation energies are accurate to <0.1 kcal/mol with respect to conventional calculations for several hundreds of electrons. This allows for the calculation of unbiased basis-set limit correlation energies for systems containing hundreds of electrons with unprecedented efficiency gains based on a straightforward treatment of continuum contributions.
| S-EPMC7584365 | biostudies-literature
Project description:In this study, we attempted to achieve to date the most comprehensive coverage of the yeast proteome during log-phase growth and in the K11R mutant ubiquitin yeast strain with a combination of (1) extensive gel molecular weight fractionation and use of either of two digestion enzymes to improve coverage by addressing the issue of dynamic range via decreased sample complexity, (2) use of a more sensitive and faster LC-MS/MS platform than available for many of the previous experiments (Orbitrap Velos coupled to a nanoAcquity UPLC). By this approach, we identified with high confidence a new fungal conserved gene expressed in S. cerevisiae and achieved the highest coverage of the yeast proteome (83.5%) to date. Besides that, we also extended our experimental workflow incorporating the same advantages to identify and quantify proteins in a K11R/reference SILAC/control total yeast lysate mixture.
2019-08-26 | PXD001928 | Pride
Project description:Gene set analysis (GSA) aims to evaluate the association between the expression of biological pathways, or a priori defined gene sets, and a particular phenotype. Numerous GSA methods have been proposed to assess the enrichment of sets of genes. However, most methods are developed with respect to a specific alternative scenario, such as a differential mean pattern or a differential coexpression. Moreover, a very limited number of methods can handle either binary, categorical, or continuous phenotypes. In this paper, we develop two novel GSA tests, called SDRs, based on the sufficient dimension reduction technique, which aims to capture sufficient information about the relationship between genes and the phenotype. The advantages of our proposed methods are that they allow for categorical and continuous phenotypes, and they are also able to identify a variety of enriched gene sets.Through simulation studies, we compared the type I error and power of SDRs with existing GSA methods for binary, triple, and continuous phenotypes. We found that SDR methods adequately control the type I error rate at the pre-specified nominal level, and they have a satisfactory power to detect gene sets with differential coexpression and to test non-linear associations between gene sets and a continuous phenotype. In addition, the SDR methods were compared with seven widely-used GSA methods using two real microarray datasets for illustration.We concluded that the SDR methods outperform the others because of their flexibility with regard to handling different kinds of phenotypes and their power to detect a wide range of alternative scenarios. Our real data analysis highlights the differences between GSA methods for detecting enriched gene sets.
| S-EPMC4744442 | biostudies-literature
Project description:The coupling of distinct systems underlies nearly all physical phenomena. A basic instance is that of interacting harmonic oscillators, giving rise to, for example, the phonon eigenmodes in a lattice. Of particular importance are the interactions in hybrid quantum systems, which can combine the benefits of each part in quantum technologies. Here we investigate a hybrid optomechanical system having three degrees of freedom, consisting of a microwave cavity and two micromechanical beams with closely spaced frequencies around 32?MHz and no direct interaction. We record the first evidence of tripartite optomechanical mixing, implying that the eigenmodes are combinations of one photonic and two phononic modes. We identify an asymmetric dark mode having a long lifetime. Simultaneously, we operate the nearly macroscopic mechanical modes close to the motional quantum ground state, down to 1.8 thermal quanta, achieved by back-action cooling. These results constitute an important advance towards engineering of entangled motional states.
| S-EPMC3432470 | biostudies-literature
Project description:MotivationPairing between the target sequence and the 6-8 nt long seed sequence of the miRNA presents the most important feature for miRNA target site prediction. Novel high-throughput technologies such as Argonaute HITS-CLIP afford meanwhile a detailed study of miRNA:mRNA duplices. These interaction maps enable a first discrimination between functional and non-functional target sites in a bulky fashion. Prediction algorithms apply different seed paradigms to identify miRNA target sites. Therefore, a quantitative assessment of miRNA target site prediction is of major interest.ResultsWe identified a set of canonical seed types based on a transcriptome wide analysis of experimentally verified functional target sites. We confirmed the specificity of long seeds but we found that the majority of functional target sites are formed by less specific seeds of only 6 nt indicating a crucial role of this type. A substantial fraction of genuine target sites arenon-conserved. Moreover, the majority of functional sites remain uncovered by common prediction methods.
| S-EPMC3087955 | biostudies-literature