We present a new implementation of the template localization for the fully automated and parallizable incremental method, which excludes failures of this important step within the domain-specific basis set approach and thus ensures a higher reliability of the scheme, preserving its very high accuracy. Furthermore, we combine our method with an efficient focal-point ansatz to reach the complete basis set limit and carefully assess this approach for the first time with regard to reaction energies. For a test set of 51 reactions the incremental focal-point method with a basis set of moderate size provides a very high accuracy with respect to the complete basis set limit.
View Article and Find Full Text PDFIn this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using T = T with appropriate values of 10 to 10 for reactions and 10 for interaction or binding energies.
View Article and Find Full Text PDFSecond-order rate constants (k) of the reactions of various barbiturate anions such as the parent barbiturate, 1,3-dimethylbarbiturate, 2-thiobarbiturate, and 1,3-diethyl-2-thiobarbiturate with diarylcarbenium ions and Michael acceptors have been determined in dimethyl sulfoxide solution at 20 °C. The reactivity parameters N and s of the barbiturate anions were derived from the linear plots of log k versus the electrophilicity parameters E of these reference electrophiles, according to the linear-free-energy relationship log k (20 °C) = s (E + N). Several reactions of these nucleophiles with benzylidenemalononitriles and quinone methides proceeded with reversible formation of the new C-C-bond followed by rate-determining proton shift.
View Article and Find Full Text PDFStacking interactions between aromatic compounds and nucleobases are crucial in recognition of nucleotides and nucleic acids, but a comprehensive understanding of the strength and selectivity of these interactions in aqueous solution has been elusive. To this end, model complexes have been designed and analyzed by experiment and theory. For the first time, stacking free energies between five nucleobases and anthracene were determined experimentally from thermodynamic double mutant cycles.
View Article and Find Full Text PDFWe present the first implementation of the fully automated incremental scheme for CCSD unrelaxed dipole moments using the domain-specific basis-set approach. Truncation parameters are varied, and the accuracy of the method is statistically analyzed for a test set of 20 molecules. The local approximations introduce small errors at second order and negligible ones at third order.
View Article and Find Full Text PDFMinimal 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.
View Article and Find Full Text PDFThe established quantum cluster equilibrium (QCE) approach is refined and applied to N-methylformamide (NMF) and its aqueous solution. The QCE method is split into two iterative cycles: one which converges to the liquid phase solution of the QCE equations and another which yields the gas phase. By comparing Gibbs energies, the thermodynamically stable phase at a given temperature and pressure is then chosen.
View Article and Find Full Text PDFJ Chem Theory Comput
January 2016
The incremental scheme makes it possible to compute CCSD(T) correlation energies to high accuracy for large systems. We present the first extension of this fully automated black-box approach to open-shell systems using an Unrestricted Hartree-Fock (UHF) wave function, extending the efficient domain-specific basis set approach to handle open-shell references. We test our approach on a set of organic and metal organic structures and molecular clusters and demonstrate standard deviations from canonical CCSD(T) values of only 1.
View Article and Find Full Text PDFJ Chem Theory Comput
August 2015
In this work we analyze the accuracy and the efficiency of different schemes to obtain the complete basis set limit for CCSD(T). It is found that composite schemes using an MP2 increment to reach the basis set limit provide high accuracy combined with high efficiency. In these composite schemes the MP2-F12/cc-pVTZ-F12 method is suitable to compute the MP2 contribution at the basis set limit.
View Article and Find Full Text PDFThe strength of interactions between CO2 and 23 building blocks of metal-organic frameworks are reported in this paper. This theoretical study is based on an incremental, explicitly correlated coupled-cluster scheme with interference effects. This scheme allows the accurate calculation of molecular complexes such as zinc acetate (32 non-hydrogen atoms) at the CCSD(T) level, close to the basis set limit.
View Article and Find Full Text PDFAs a rule, acetylides and sulfonyl azides do not undergo electrophilic azide transfer because 1,2,3-triazoles are usually formed. We show now that treatment of tritylethyne with butyllithium followed by exposure to 2,4,6-triisopropylbenzenesulfonyl azide leads to products that are easily explained through the generation of short-lived tritylethynyl azide and its secondary product cyanotritylcarbene. Furthermore, it is demonstrated that tritylcarbenes generally do not produce triphenylethenes exclusively, as was stated in the literature.
View Article and Find Full Text PDF4,5-Dihydro-1,2,3-oxadiazoles are postulated to be key intermediates in the industrial synthesis of ketones from alkenes, in the alkylation of DNA in vivo, and in the decomposition of N-nitrosoureas; they are also a subject of great interest for theoretical chemists. In the presented report, the formation of 4,5-dihydro-1,2,3-oxadiazole and the subsequent decay into secondary products have been studied by NMR monitoring analysis. The elusive properties evading characterization have now been confirmed by (1) H, (13) C, and (15) N NMR spectroscopy, and relevant 2D experiments at very low temperatures.
View Article and Find Full Text PDFIn this work we present the first implementation of the incremental scheme for coupled cluster linear-response frequency-dependent dipole polarizabilities. The implementation is fully automated and makes use of the domain-specific basis set approach. The accuracy of the approach is determined on the basis of a test suite of 47 molecules and small clusters.
View Article and Find Full Text PDFThe parametrization of the polarized molecular orbital (PMO) method, which is a neglect-of-diatomic-differential-overlap (NDDO) semiempirical method that includes polarization functions on hydrogens, is extended to include the constituents that dominate the nucleation of atmospheric aerosols, including ammonia, sulfuric acid, and water. The parametrization and validation are based mainly on CCSD(T)/CBS results for atmospheric clusters composed of sulfuric acid, dimethylamine, and ammonia and on M06-2X exchange-correlation functional calculations for other constituents of the atmospheric aerosols. The resulting model, called PMO2a, is parametrized for molecules containing any type of H, C, or O, amino or ammonium N, and S atoms bonded to O.
View Article and Find Full Text PDFIt is important to test methods for simulating water, but small water clusters for which benchmarks are available are not very representative of the bulk. Here we present benchmark calculations, in particular CCSD(T) calculations at the complete basis set limit, for water 26-mers drawn from Monte Carlo simulations of bulk water. These clusters are large enough that each water molecule participates in 2.
View Article and Find Full Text PDFIn this work, we use MP2 and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] as well as their corresponding explicitly correlated (F12) counterparts to compute the interaction energies of water icosamers. The incremental scheme is used to compute benchmark energies at the CCSD(T)/CBS(45) and CCSD(T)(F12*)/cc-pVQZ-F12 level of theory. The four structures, dodecahedron, edge sharing, face sharing, and fused cubes, are part of the WATER27 test set and therefore, highly accurate interaction energies are required.
View Article and Find Full Text PDFJ Chem Theory Comput
December 2013
In this work we present a new partitioning scheme for the incremental approach in combination with the efficient (F12*) approximation for explicitly correlated coupled cluster (J. Chem. Phys.
View Article and Find Full Text PDFA series of ferrocenyl (Fc = ferrocenyl; fc = ferrocen-1,1'-diyl) and biferrocenyl (Bfc = 1',1″-biferrocenyl; bfc = 1',1″-biferrocen-1,1‴-diyl) mono- and biscarbene tungsten(0) complexes of the type [(CO)5W═C(OMe)R] (1, R = Fc; 3, R = Bfc) and [(CO)5W═C(OMe)-R'-(OMe)C═W(CO)5] (2, R' = fc; 4, R' = bfc) were synthesized according to the classical synthetic methodology by reacting W(CO)6 with LiR (R = Fc, fc, bfc), followed by a subsequent alkylation using methyl trifluoromethanesulfonate. Electrochemical investigations were carried out on these complexes to get a closer insight into the electronic properties of 1-4. The ferrocenyl and biferrocenyl moieties in 1-4 show reversible one-electron redox events.
View Article and Find Full Text PDFWe present a theoretical study of carbene formation from the 1-ethyl-3-methylimidazolium acetate ionic liquid in the absence and presence of CO2 in gas and liquid phase. Although CO2 physisorption constitutes a precursory step of chemisorption (the CO2's reaction with carbenes, which forms from cations via proton abstraction by anions), it also enables a very stable CO2-anion associate. However, this counteracts the chemical absorption by reducing the basicity of the anion and the electrophilicity of the CO2, which is reflected by charge transfer.
View Article and Find Full Text PDFJ Chem Theory Comput
January 2013
In this work, we apply the recently proposed MP2 correction to incremental energies within the domain-specific basis set approach to incrementally expanded CCSD(T)(F12) energies. The approach is tested for a set of 27 molecules with different electronic structures including water clusters, aqua complexes, aliphatic hydrocarbons, alkenes, alkynes, aromatic systems, and amino acids. The root mean squared deviation of the absolute energies with respect to the standard calculation is 1.
View Article and Find Full Text PDFThe mechanism of catalytic 4-exo cyclizations without gem-dialkyl substitution was investigated by a comparison of cyclic voltammetry, EPR, and computational studies with previously published synthetic results. The most active catalyst is a super-unsaturated 13-electron titanocene(III) complex that is formed by supramolecular activation through hydrogen bonding. The template catalyst binds radicals via a two-point binding that is mandatory for the success of the 4-exo cyclization.
View Article and Find Full Text PDF