Search

We present a natural orbital-based implementation of the intermediate Hamiltonian Fock space coupled-cluster method for the (1, 1) sector of Fock space. The use of natural orbitals significantly reduces the computational cost and can automatically choose an appropriate set of active orbitals. The new method retains the charge transfer separability of the original intermediate Hamiltonian Fock space coupled-cluster method and gives excellent performance for valence, Rydberg, and charge-transfer excited states. It offers significant computational advantages over the popular equation of motion coupled cluster method for excited states dominated by single excitations. [ABSTRACT FROM AUTHOR]

This work describes a domain-based local pair natural orbital (DLPNO) implementation of the equation of motion coupled cluster method for the computation of electron affinities (EAs) including single and double excitations. Similar to our earlier work on ionization potentials (IPs), the method reported in this study uses the ground state DLPNO framework and extends it to the electron attachment problem. While full linear scaling could not be achieved as in the IP case, leaving the Fock/Koopmans' contributions in the canonical basis and using a tighter threshold for singles PNOs allows us to compute accurate EAs and retain most of the efficiency of the DLPNO technique. Thus as in the IP case, the ground state truncation parameters are sufficient to control the accuracy of the computed EA values, although a new set of integrals for singles PNOs must be generated at the DLPNO integral transformation step. Using standard settings, our method reproduces the canonical results with a maximum absolute deviation of 49 meV for bound states of a test set of 24 molecules. Using the same settings, a calculation involving more than 4500 basis functions, including diffuse functions, takes four days on four cores, with only 48 min spent in the EA module itself. [ABSTRACT FROM AUTHOR]

We present a benchmark investigation on the performance of various approximate EOM‐CCSD methods for the ionization potential and electron affinity of the genetic material subunits. Our test set consists of 22 molecules, and the canonical EOM‐CCSD method has been used as the reference. Our study shows that the accuracy of approximate methods varies largely for different quantities. However, the pair natural orbital (PNO) based EOM‐CCSD method gives a proper balance between computational cost and accuracy, more suitable than the methods based on perturbative approximation to study properties related to ionization and electron attachment to genetic materials. [ABSTRACT FROM AUTHOR]

In this work, a domain-based local pair natural orbital (DLPNO) version of the equation of motion coupled cluster theory with single and double excitations for ionization potentials (IP-EOM-CCSD) equations has been formulated and implemented. The method uses ground state localized occupied and pair natural virtual orbitals and applies the DLPNO machinery to arrive at a linear scaling implementation of the IP-EOM-CCSD method. The accuracy of the method is controllable using ground state truncation parameters. Using default thresholds, the method predicts ionization potential (IP) values with good accuracy (mean absolute error of 0.08 eV). We demonstrate that our code can be used to compute IP values for systems with more than 1000 atoms and 10 000 basis functions. [ABSTRACT FROM AUTHOR]

We have developed a reduced-cost algebraic diagrammatic construction (ADC) method based on state-specific frozen natural orbital and natural auxiliary functions. The newly developed method has been benchmarked on the GW100 test set for the ionization problem. The use of state-specific natural orbitals drastically reduces the size of the virtual space with a systematically controllable accuracy and offers a significant speedup over the standard ionization potential (IP)-ADC(3) method. The accuracy of the method can be controlled by two thresholds and nearly a black box to use. The inclusion of the perturbative correction significantly improves the accuracy of the calculated IP values, and the efficiency of the method has been demonstrated by calculating the IP of a molecule with 60 atoms and more than 2216 basis functions. [ABSTRACT FROM AUTHOR]

An efficient scheme for the automatic selection of an active space for similarity transformed equations of motion (STEOM) coupled cluster method is proposed. It relies on state averaged configuration interaction singles (CIS) natural orbitals and makes it possible to use STEOM as a black box method. The performance of the new scheme is tested for singlet and triplet valence, charge transfer, and Rydberg excited states. [ABSTRACT FROM AUTHOR]

The use of back-transformed pair natural orbitals in the calculation of excited state energies, ionization potentials, and electron affinities is investigated within the framework of equation of motion coupled cluster theory and its similarity transformed variant. Possible approaches to a more optimal use of pair natural orbitals in these methods are indicated. [ABSTRACT FROM AUTHOR]

In the present paper, the chain of spheres exchange (COSX) approximation is applied to the highest scaling terms in the equation of motion (EOM) coupled cluster equations with single and double excitations, in particular, the terms involving integrals with four virtual labels. It is found that even the acceleration of this single term yields significant computational gains without compromising the desired accuracy of the method. For an excitation energy calculation on a cluster of five water molecules using 585 basis functions, the four virtual term is 9.4 times faster using COSX with a loose grid than using the canonical implementation, which yields a 2.6 fold acceleration for the whole of the EOM calculation. For electron attachment calculations, the four virtual term is 15 times and the total EOM calculation is 10 times faster than the canonical calculation for the same system. The accuracy of the new method was tested using Thiel's test set for excited states using the same settings and the maximum absolute deviation over the whole test set was found to be 12.945 cm−1 (59 μHartree) for excitation energies and 6.799 cm−1 (31 μHartree) for electron attachments. Using MP2 amplitudes for the ground state in combination with the parallel evaluation of the full EOM equations in the manner discussed in this paper enabled us to perform calculations for large systems. Electron affinity values for the two lowest states of a Zn protoporphyrine model compound (224 correlated electrons and 1120 basis functions) were obtained in 3 days 19 h using 4 cores of a Xeon E5-2670 processor allocating 10 GB memory per core. Calculating the lowest two excitation energies for trans-retinal (114 correlated electrons and 539 basis functions) took 1 day 21 h using eight cores of the same processor and identical memory allocation per core. [ABSTRACT FROM AUTHOR]

We propose a new elegant strategy to implement third order triples correction in the light of many-body perturbation theory to the Fock space multi-reference coupled cluster method for the ionization problem. The computational scaling as well as the storage requirement is of key concerns in any many-body calculations. Our proposed approach scales as N6 does not require the storage of triples amplitudes and gives superior agreement over all the previous attempts made. This approach is capable of calculating multiple roots in a single calculation in contrast to the inclusion of perturbative triples in the equation of motion variant of the coupled cluster theory, where each root needs to be computed in a state-specific way and requires both the left and right state vectors together. The performance of the newly implemented scheme is tested by applying to methylene, boron nitride (B2N) anion, nitrogen, water, carbon monoxide, acetylene, formaldehyde, and thymine monomer, a DNA base. [ABSTRACT FROM AUTHOR]

We present a benchmark theoretical investigation on the electronic structure and singlet-triplet(ST) gap of 1- and 2-naphthyl cations using the CCSD(T) method. Our calculations reveal that the ground states of both the naphthyl cations are singlet, contrary to the results obtained by DFT/B3LYP calculations reported in previous theoretical studies. However, the triplet states obtained in the two structural isomers of naphthyl cation are completely different. The triplet state in 1-naphthyl cation is (π,σ) type, whereas in 2-naphthyl cation it is (σ,σ') type. The S-T gaps in naphthyl cations and the relative stability ordering of the singlet and the triplet states are highly sensitive to the basisset quality as well as level of correlation, and demand for inclusion of perturbative triples in the coupled-cluster ansatz. [ABSTRACT FROM AUTHOR]

Spin flip equation of motion coupled cluster (EOM-SF-CC) can correctly treat situations involving electronic degeneracies or near degeneracies, e.g., bond breaking, di- and tri-radicals, etc. However, for large systems EOM-SF-CC (even in single and double excitations) is computationally prohibitively expensive. Therefore, earlier approximations to EOM-SF-CC methods such as spin flip configuration interaction singles with perturbative doubles (SF-CIS(D)) have been proposed. In this work, we present a new perturbative approximation to EOM-SF-CC, which has been found to be more accurate than SF-CIS(D). The capabilities, advantages, and timings of the new approach have been demonstrated considering the singlet-triplet gaps in di- and triradicals as well as bond breaking examples. The method is extended to double spin flip EOM-CC and its capabilities have been tested. [ABSTRACT FROM AUTHOR]

Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost‐effective implementations of correlated wave function methods for both closed‐shell and open‐shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain‐based DLPNO approximation, and the equations‐of‐motion approach for electron affinities, EA‐EOM‐CCSD. In addition, several variants of canonical, spin‐component‐scaled, orbital‐optimized, and explicitly correlated (F12) Møller–Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO‐based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO‐pCCSD, the DLPNO‐CCSD(T0), and the orbital optimized OO‐SCS‐MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost‐effective DLPNO‐based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc. A set of biologically relevant quinones with experimentally known electron affinities was compiled to evaluate the predictive ability of a range of wave function and density functional methods. Particular emphasis is placed on the performance characteristics of modern approximations to coupled cluster theory, which hold the promise of delivering systematically and consistently high accuracy for ever larger and chemically diverse systems. [ABSTRACT FROM AUTHOR]

Abstract: In this article, we investigate the performance of different approximate variants of the EOM‐CCSD method for calculation of ionization potential (IP), as compared to EOM‐CCSDT reference values. None of the lower scaling approximations to the EOM‐CCSD method give a consistent performance for valence, inner valence, and core ionization, favoring one, or the other depending on the nature of the approximation used. The parent EOMIP‐CCSD method gives superior performance for valence IP but can show large errors for inner valence and core ionization. The problem is particularly severe for core‐ionization, where even the EOMIP‐CCSDT method cannot provide quantitative accuracy. [ABSTRACT FROM AUTHOR]

In this paper, we present a chain-of-sphere implementation of the external exchange term, the computational bottleneck of coupled-cluster calculations at the singles and doubles level. This implementation is compared to standard molecular orbital, atomic orbital and resolution of identity implementations of the same term within the ORCA package and turns out to be the most efficient one for larger molecules, with a better accuracy than the resolution-of-identity approximation. Furthermore, it becomes possible to perform a canonical CC calculation on a tetramer of nucleobases in 17 days, 20 hours. [ABSTRACT FROM AUTHOR]

We present a novel and cost-effective approach of using a second similarity transformation of the Hamiltonian to include the missing higher-order terms in the second-order approximate coupled cluster singles and doubles (CC2) model. The performance of the newly developed ST-EOM-CC2 model has been investigated for the calculation of excitation energies of valence, Rydberg, and charge-transfer excited states. The method shows significant improvement in the excitation energies of Rydberg and charge-transfer excited states as compared to the conventional CC2 method while retaining the good performance of the latter for the valence excited state. This method retains the charge-transfer separability of the charge-transfer excited states, which is a significant advantage over the traditional CC2 method. A second order many-body perturbation theory variant of the new method is also proposed. [ABSTRACT FROM AUTHOR]

We have investigated the impact of microsolvation on shape-type resonance states of nucleobases, taking cytosine as a case study. To characterize the resonance position and decay width of the metastable states, we employed the newly developed DLPNO-based EA-EOM-CCSD method in conjunction with the resonance via Padé (RVP) method. Our calculations show that the presence of water molecules causes a redshift in the resonance position and an increase in the lifetime for the three lowest-lying resonance states of cytosine. Furthermore, there are some indications that the lowest resonance state in isolated cytosine may get converted to a bound state in the presence of an aqueous environment. The obtained results are extremely sensitive to the basis set used for the calculations. [ABSTRACT FROM AUTHOR]

A perturbative analysis of the ground-state similarity transformed Hamiltonian and its effect on the various Fock-space coupled cluster (FSCC) sectors is presented through calculation of ionisation potential, electron affinity, excitation energies and response properties. Various truncation schemes of the effective Hamiltonian are presented with explicit form of the defining equations. Based on such a truncation, the approximate methods are labelled as FSCC(n), wherenrepresents the correlation energy of the ionised, electron attached or excited states corrected at least up tonth order within coupled cluster singles and doubles scheme (CCSD). A lower scaling CC2 type of approach (abbreviated as FS-CC2) is compared against the group of FSCC(n) methods for energies. Electric response properties have been compared and contrasted for the two lower scaling methods: FSCC(2) and FS-CC2. The various truncated methods are tested for a number of small molecules. The results obtained from a range of truncated methods are compared against full FSCCSD calculations. [ABSTRACT FROM PUBLISHER]

We report a benchmark theoretical investigation of both vertical and adiabatic electron affinities of DNA and RNA nucleobases: adenine, guanine, cytosine, thymine, and uracil using equation of motion coupled cluster method. The vertical electron affinity (VEA) values of the first five states of the DNA and RNA nucleobases are computed. It is observed that the first electron attached state is energetically accessible in gas phase. Furthermore, an analysis of the natural orbitals exhibits that the first electron attached states of uracil and thymine are valence-bound in nature and undergo significant structural changes on attachment of an extra electron, which reflects in the deviation of the adiabatic electron affinity (AEA) than that of the vertical ones. Conversely, the first electron attached states of cytosine, adenine, and guanine are in the category of dipole-bound anions. Their structure, by and large, remain unaffected on attachment of an extra electron, which is evident from the observed small difference between the AEA and VEA values. VEA and AEA values of all the DNA and RNA nucleobases are found to be negative, which implies that the first electron attached states are not stable rather quasi bound. The results of all previous theoretical calculations are out of track and shows large deviation with respect to the experimentally measured values, whereas, our results are found to be in good agreement. Therefore, our computed values can be used as a reliable standard to calibrate new theoretical methods. © 2015 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

The Fock-space multi-reference coupled cluster method is used for the geometry optimisation of the low-lying excited states of the molecules. Molecular geometries of the carbon monohydride cation (CH)+, water (H2O), ozone (O3) and formaldehyde (HCHO) in their low-lying excited states are optimised. Excited state gradients are calculated using finite field multi-reference coupled cluster method. We compare our results with other theoretical and/or experimental results, wherever available. [ABSTRACT FROM PUBLISHER]

We have presented a mechanism for electron attachment to solvated nucleobases using accurate wave-function based hybrid quantum/classical (QM/MM) simulations and uracil as a test case. The initial electron attached state is found to be localized in the bulk water, and this water-bound state acts as a doorway to the formation of the final nucleobase bound state. The electron transfer from water to uracil takes place because of the mixing of electronic and nuclear degrees of freedom. The water molecules around the uracil stabilize the uracil-bound anion by creating an extensive hydrogen-bonding network and accelerate the rate of electron attachment to uracil. The complete transfer of the electron from water to the uracil occurs in a picosecond time scale, which is consistent with the experimentally observed rate of reduction of nucleobases in the presence of water. The degree of solvation of the aqueous electron can lead to a difference in the initial stabilization of the uracil-bound anion. However, the anions formed due to the attachment of both surface-bound and bulk-solvated electrons behave similarly to each other at a longer time scale. [ABSTRACT FROM AUTHOR]

We report a comparative single-referenceand multireference coupled-cluster investigation on the structure,potential energy surface, and IR spectroscopic properties of the transperoxo nitrate radical, one of the key intermediates in stratosphericNOXchemistry. The previous single-reference ab initiostudies predicted an unbound structure for the trans peroxo nitrateradical. However, our Fock space multireference coupled-cluster calculationconfirms a bound structure for the trans peroxo nitrate radical, inaccordance with the experimental results reported earlier. Further,the analysis of the potential energy surface in FSMRCC method indicatesa well-behaved minima, contrary to the shallow minima predicted bythe single-reference coupled-cluster method. The harmonic force fieldanalysis, of various possible isomers of peroxo nitrate also revealsthat only the trans structure leads to the experimentally observedIR peak at 1840 cm–1. The present study highlightsthe critical importance of nondynamic correlation in predicting thestructure and properties of high-energy stratospheric NOxradicals. [ABSTRACT FROM AUTHOR]

We present the theory and the implementation of a low-cost four-component relativistic equation of motion coupled cluster method for ionized states based on frozen natural spinors. A single threshold (natural spinor occupancy) can control the accuracy of the calculated ionization potential values. Frozen natural spinors can significantly reduce the computational cost for valence and core-ionization energies with systematically controllable accuracy. The convergence of the ionization potential values with respect to the natural spinor occupancy threshold becomes slower with the increase in basis set dimension. However, the use of a natural spinor threshold of 10−5 and 10−6 gives excellent agreement with experimental results for valence and core ionization energies, respectively. [ABSTRACT FROM AUTHOR]

In this article, we have studied the bound anionic states of DNA and RNA nucleobases using the newly developed bt‐PNO‐EOM‐CCSD method and extended basis sets. All of the five nucleobases have a single bound anionic state, which is dipole bound in nature. The electron affinity corresponding to the dipole‐bound state is found to be extremely sensitive to the used basis set, and a proper agreement with the available experimental data requires a large basis set with a sufficient number of diffuse functions. Electron correlation plays a major role in stabilizing the bound anionic states of nucleobases. No valence‐bound anionic states are observed for any of the nucleobases. All the five nucleobases of DNA and RNA give a single bound anionic state on electron attachment, which is dipole bound in nature. It is characterized by electron density away from nuclear framework. Thus, simulated electron affinity values are very sensitive to the number of diffuse functions present in the basis set, and electron correlation plays a crucial role in determining their accuracy. [ABSTRACT FROM AUTHOR]

We present the theory, implementation, and benchmark results for a frozen natural spinors based reduced cost four-component relativistic coupled cluster method. The natural spinors are obtained by diagonalizing the one-body reduced density matrix from a relativistic second-order Møller–Plesset calculation based on a four-component Dirac–Coulomb Hamiltonian. The correlation energy in the coupled cluster method converges more rapidly with respect to the size of the virtual space in the frozen natural spinor basis than that observed in the standard canonical spinors obtained from the Dirac–Hartree–Fock calculation. The convergence of properties is not smooth in the frozen natural spinor basis. However, the inclusion of the perturbative correction smoothens the convergence of the properties with respect to the size of the virtual space in the frozen natural spinor basis and greatly reduces the truncation errors in both energy and property calculations. The accuracy of the frozen natural spinor based coupled cluster methods can be controlled by a single threshold and is a black box to use. [ABSTRACT FROM AUTHOR]

A cost-effective equation of motion coupled cluster method, EOMIP-CCSD(2), is used to investigate vertical and adiabatic ionization potential as well as ionization-induced structural changes of water clusters and compared with CCSD(T), CASPT2, and MP2 methods. The moderate N5 scaling and low storage requirement yields EOMIP-CCSD(2) calculation feasible even for reasonably large molecules and clusters with accuracy comparable to CCSD(T) method at much cheaper computational cost. Our calculations shed light on the authenticity of EOMIP-CCSD(2) results and establish a reliable method to study of ionization energy of molecular clusters. We have further investigated the performance of several classes of DFT functionals for ionization energies of water clusters to benchmark the results and to get a reliable functionals for the same. [ABSTRACT FROM AUTHOR]

Intermolecular and intramolecular Friedel-Crafts acylation reactions of various aliphatic and aromatic esters at room temperature with the use of very simple reagents and activating groups in are described. The products were obtained in good yield (60-85 %). The detailed mechanistic pathway was studied by DFT calculations and supported by experimental evidence. [ABSTRACT FROM AUTHOR]

This method can be used for intermolecular as well as intramolecular Friedel-Crafts acylations, e. [ABSTRACT FROM AUTHOR]