Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations.
Ontology highlight
ABSTRACT:
SUBMITTER: Altun A
PROVIDER: S-EPMC8607505 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method allows one to perform single point energy calculations for systems with hundreds of atoms while retaining essentially the accuracy of its canonical counterpart, with errors that are typically smaller than 1 kcal/mol for relative energies. Crucial to the accuracy and efficiency of the method is a proper definition of the virtual space in which the coupled-cluster equations are solved, which is spanned by a highly compact set of pair natural orbitals (PNOs) that are specific for each electron pair. The dimension of the PNO space is controlled by the TCutPNO threshold: only PNOs with an occupation number greater than TCutPNO are included in the correlation space of a given electron pair, whilst the remaining PNOs are discarded. To keep the error of the method small, a conservative TCutPNO value is used in standard DLPNO-CCSD(T) calculations. This often leads to unnecessarily large PNO spaces, which limits the efficiency of the method. Herein, we introduce a new computational strategy to approach the complete PNO space limit (for a given basis set) that consists in extrapolating the results obtained with different TCutPNO values. The method is validated on the GMTKN55 set using canonical CCSD(T) data as the reference. Our results demonstrate that a simple two-point extrapolation scheme can be used to significantly increase the efficiency and accuracy of DLPNO-CCSD(T) calculations, thus extending the range of applicability of the technique.
| S-EPMC7586325 | biostudies-literature
Project description:In recent years, the number of studies using a cluster-randomized design has grown dramatically. In addition, the cluster-randomized crossover design has been touted as a methodological advance that can increase efficiency of cluster-randomized studies in certain situations. While the cluster-randomized crossover trial has become a popular tool, standards of design, analysis, reporting and implementation have not been established for this emergent design. We address one particular aspect of cluster-randomized and cluster-randomized crossover trial design: estimating statistical power. We present a general framework for estimating power via simulation in cluster-randomized studies with or without one or more crossover periods. We have implemented this framework in the clusterPower software package for R, freely available online from the Comprehensive R Archive Network. Our simulation framework is easy to implement and users may customize the methods used for data analysis. We give four examples of using the software in practice. The clusterPower package could play an important role in the design of future cluster-randomized and cluster-randomized crossover studies. This work is the first to establish a universal method for calculating power for both cluster-randomized and cluster-randomized clinical trials. More research is needed to develop standardized and recommended methodology for cluster-randomized crossover studies.
| S-EPMC3338707 | biostudies-literature
Project description:The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].
| S-EPMC4174990 | biostudies-literature
Project description:The SCISSORS method for approximating chemical similarities has shown excellent empirical performance on a number of real-world chemical data sets but lacks theoretically proven bounds on its worst-case error performance. This paper first proves reductions showing SCISSORS to be equivalent to two previous kernel methods: kernel principal components analysis and the rank-k Nystro?m approximation of a Gram matrix. These reductions allow the use of generalization bounds on these techniques to show that the expected error in SCISSORS approximations of molecular similarity kernels is bounded in expected pairwise inner product error, in matrix 2-norm and Frobenius norm for full kernel matrix approximations and in root-mean-square deviation for approximated matrices. Finally, we show that the actual performance of SCISSORS is significantly better than these worst-case bounds, indicating that chemical space is well-structured for chemical sampling algorithms.
| S-EPMC3183166 | biostudies-literature
Project description:We consider the problem of calculating power and sample size for tests based on generalized estimating equations (GEE), that arise in studies involving clustered or correlated data (e.g., longitudinal studies and sibling studies). Previous approaches approximate the power of such tests using the asymptotic behavior of the test statistics under fixed alternatives. We develop a more accurate approach in which the asymptotic behavior is studied under a sequence of local alternatives that converge to the null hypothesis at root-m rate, where m is the number of clusters. Based on this approach, explicit sample size formulae are derived for Wald and quasi-score test statistics in a variety of GEE settings. Simulation results show that in the important special case of logistic regression with exchangeable correlation structure, previous approaches can inflate the projected sample size (to obtain nominal 90% power using the Wald statistic) by over 10%, whereas the proposed approach provides an accuracy of around 2%.
| S-EPMC3903421 | biostudies-literature
Project description:When calculating sample size or power for stepped wedge or other types of longitudinal cluster randomized trials, it is critical that the planned sampling structure be accurately specified. One common assumption is that participants will provide measurements in each trial period, that is, a closed cohort, and another is that each participant provides only one measurement during the course of the trial. However some studies have an "open cohort" sampling structure, where participants may provide measurements in variable numbers of periods. To date, sample size calculations for longitudinal cluster randomized trials have not accommodated open cohorts. Feldman and McKinlay (1994) provided some guidance, stating that the participant-level autocorrelation could be varied to account for the degree of overlap in different periods of the study, but did not indicate precisely how to do so. We present sample size and power formulas that allow for open cohorts and discuss the impact of the degree of "openness" on sample size and power. We consider designs where the number of participants in each cluster will be maintained throughout the trial, but individual participants may provide differing numbers of measurements. Our results are a unification of closed cohort and repeated cross-sectional sample results of Hooper et al (2016), and indicate precisely how participant autocorrelation of Feldman and McKinlay should be varied to account for an open cohort sampling structure. We discuss different types of open cohort sampling schemes and how open cohort sampling structure impacts on power in the presence of decaying within-cluster correlations and autoregressive participant-level errors.
| S-EPMC7217159 | biostudies-literature
Project description:Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O(N4) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.
| S-EPMC10569055 | biostudies-literature
Project description:ObjectivesTo clarify and illustrate sample size calculations for the cross-sectional stepped wedge cluster randomized trial (SW-CRT) and to present a simple approach for comparing the efficiencies of competing designs within a unified framework.Study design and settingWe summarize design effects for the SW-CRT, the parallel cluster randomized trial (CRT), and the parallel cluster randomized trial with before and after observations (CRT-BA), assuming cross-sectional samples are selected over time. We present new formulas that enable trialists to determine the required cluster size for a given number of clusters. We illustrate by example how to implement the presented design effects and give practical guidance on the design of stepped wedge studies.ResultsFor a fixed total cluster size, the choice of study design that provides the greatest power depends on the intracluster correlation coefficient (ICC) and the cluster size. When the ICC is small, the CRT tends to be more efficient; when the ICC is large, the SW-CRT tends to be more efficient and can serve as an alternative design when the CRT is an infeasible design.ConclusionOur unified approach allows trialists to easily compare the efficiencies of three competing designs to inform the decision about the most efficient design in a given scenario.
| S-EPMC4687983 | biostudies-literature
Project description:Herein, we introduce a fragment-based local coupled cluster embedding approach for the accurate quantification and analysis of noncovalent interactions in molecular aggregates. Our scheme combines two different expansions of the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) energy: the many-body expansion (MBE) and the local energy decomposition (LED). The low-order terms in the MBE are initially computed in the presence of an environment that is treated at a low level of theory. Then, LED is used to decompose the energy of each term in the embedded MBE into additive fragment and fragment-pairwise contributions. This information is used to quantify the total energy of the system while providing at the same time in-depth insights into the nature and cooperativity of noncovalent interactions. Two different approaches are introduced and tested, in which the environment is treated at different levels of theory: the local coupled cluster in the Hartree-Fock (LCC-in-HF) method, in which the environment is treated at the HF level; and the electrostatically embedded local coupled cluster method (LCC-in-EE), in which the environment is replaced by point charges. Both schemes are designed to preserve as much as possible the accuracy of the parent local coupled cluster method for total energies, while being embarrassingly parallel and less memory intensive. These schemes appear to be particularly promising for the study of large and complex molecular aggregates at the coupled cluster level, such as condensed phase systems and protein-ligand interactions.
| S-EPMC8190956 | biostudies-literature
Project description:The performance of error correction protocols are necessary for understanding the operation of potential quantum computers, but this requires physical error models that can be simulated efficiently with classical computers. The Gottesmann-Knill theorem guarantees a class of such error models. Of these, one of the simplest is the Pauli twirling approximation (PTA), which is obtained by twirling an arbitrary completely positive error channel over the Pauli basis, resulting in a Pauli channel. In this work, we test the PTA's accuracy at predicting the logical error rate by simulating the 5-qubit code using a 9-qubit circuit with realistic decoherence and unitary gate errors. We find evidence for good agreement with exact simulation, with the PTA overestimating the logical error rate by a factor of 2 to 3. Our results suggest that the PTA is a reliable predictor of the logical error rate, at least for low-distance codes.
| S-EPMC4588586 | biostudies-other