Thomas Watson

Benchmark studies on the building blocks of DNA. 1. Superiority of coupled cluster methods in describing the excited states of nucleobases in the Franck-Condon region.

Equation of motion excitation energy coupled-cluster (EOMEE-CC) methods including perturbative triple excitations have been used to set benchmark results for the excitation energy and oscillator strength of the building units of DNA, i.e., cytosine, guanine, adenine and thymine. In all cases the lowest twelve transitions have been considered including valence and Rydberg ones. Triple-ζ basis sets with diffuse functions have been used and the results are compared to CC2, CASPT2, TDDFT, and DFT/MRCI results from the literature. The results clearly show that it is only the EOMEE-CCSD(T) that is capable of providing accuracy of about 0.1 eV. EOMEE-CCSD systematically overshoots the energy of all types of transitions by 0.1-0.3 eV, whereas CC2 is surprisingly accurate for ππ* transitions but fails (often badly) for nπ* and Rydberg transitions. DFT and CASPT2 seem to give reliable results for the lowest transition, but the error increases fast with the excitation level. The differences in the excitation energies often change the energy ordering of the states, which should even influence the conclusions of excited state dynamics obtained with these approximate methods. The results call for further benchmark calculations on larger building blocks of DNA (nucleosides, basis pairs) at the CCSD(T) level.

Benchmarking for Perturbative Triple-Excitations in EE-EOM-CC Methods

Perturbative triples corrections ((T) and (T̃)) to the equation of motion coupled cluster singles and doubles (EOM-CCSD) are rederived and implemented in a pilot parallel code. The vertical excitation energies of molecules in the test set of Sauer et al. [J. Chem. Theor. Comput. 2009, 5, 555] are reported and compared to the iterative EOM-CCSDT-3 method. The average absolute deviations of EOM-CCSD(T) and EOM-CCSD(T̃) from EOM-CCSDT-3 over this wide test set are 0.06 and 0.18 eV, respectively. The poor performance of the latter suggests misbalanced handling of the (T̃) terms. Scaling curves for EOM-CCSD(T) are also presented to demonstrate the performance across multiple compute nodes, thus enabling the routine and accurate study of excited states for ever larger molecular systems.

RDX Geometries, Excited States, and Revised Energy Ordering of Conformers via MP2 and CCSD(T) Methodologies: Insights into Decomposition Mechanism

The geometries, harmonic frequencies, elec-tronic excitation levels, and energetic orderings of various conformers of RDX have been computed at the ab initio MP2 and CCSD(T) levels, providing more reliable results than have been previously obtained. We observe that the various local minimum-energy conformers are all competitive for being the absolute minimum and that, at reasonable temperatures, several conformers will appreciably contribute to the population of RDX. As a result, we have concluded that any mechanistic study to investigate thermal decomposition can reasonably begin from any one of the cyclohexane conformers of RDX. As such, it is necessary to consider the transition states for each RDX conformer to gauge what the activation energy is. Homolytic bond dissociation has long been speculated to be critical to detonation; we report here the most accurate estimates of homolytic BDEs yet calculated, likely to be accurate within 3 kcal mol(-1). The differences in energy for homolytic BDEs among all the possible RDR conformers are again small, such that most all of the conformers can reasonably be speculated as the next step in the mechanism starting from the RDR radical.

Benchmark studies on the building blocks of DNA. 2. Effect of biological environment on the electronic excitation spectrum of nucleobases.

In the first paper of this series (Szalay; et al. J. Phys. Chem. A, 2012, 116, 6702) we have investigated the excited states of nucleobases. It was shown that it is only the equation of motion excitation energy coupled-cluster (EOMEE-CC) methods, which can give a balanced description for all type of the transitions of these molecules; if the goal is to obtain accurate results with uncertainty of about 0.1 eV only, triples corrections in the form of, e.g., the EOMEE-CCSD(T) method need to be included. In this second paper we extend this study to nucleobases in their biological environment, considering hydration, glycoside bond, and base pairing. EOMEE-CCSD and EOMEE-CCSD(T) methods are used with aug-cc-pVDZ basis. The effect of surrounding water was systematically investigated by considering one to five water molecules at different positions. It was found that hydration can modify the order of the excited states: in particular, nπ* states get shifted above the neighboring ππ* ones. The glycoside bonds effect is smaller, as shown by our calculations on cytidine and guanosine. Here the loss of planarity causes some intensity shift from ππ* to nπ* states. Finally, the guanine-cytosine (GC) Watson-Crick pair was studied; most of the states could be identified as local excitations on one of the bases, but there is also a low-lying charge-transfer state. Significant discrepancy with earlier CASPT2 and TDDFT studies was found for the GC pair and triples effects seem to be essential for all of these systems.

Benchmark studies on the building blocks of DNA. 3. Watson-Crick and stacked base pairs.

Excited states of stacked adenine-thymine and guanine-cytosine pairs as well as the Watson-Crick pair of guanine-thymine have been investigated using the equation of motion coupled-cluster (EOM-CC) method with single and double as well as approximate triple excitations. Transitions have been assigned, and the form of the excitations has been analyzed. The majority of the excitations could be classified as localized on the nucleobases, but for all three studied systems, charge-transfer (CT) transitions could also be identified. The main aim of this study was to compare the performance of lower-level methods (ADC(2) and TDDFT) to the high-level EOM-CC ones. It was shown that both ADC(2) and TDDFT with long-range correction have nonsystematic error in excitation energies, causing alternation of the energetic ordering of the excitations. Considering the high costs of the EOM-CC calculations, there is a need for reliable new approximate methods.

Conformers of CL-20 explosive and ab initio refinement using perturbation theory: implications to detonation mechanisms.

We have identified the major conformers of CL-20 explosive, otherwise known as 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane, more formally known as 2,4,6,8,10,12-hexanitrohexaazatetracyclo[5.5.0.0]-dodecane, via Monte Carlo search in conformational space through molecular mechanics and subsequent quantum mechanical refinement using perturbation theory. Our search produced enough conformers to account for all of the various forms of CL-20 found in crystals. This suggests that our methodology will be useful in studying the conformational landscape of other nitramines. The energy levels of the conformers found are all within 0.25 eV of one another based on MBPT(2)/6-311G(d,p); consequently, without further refinement from a method such as coupled cluster theory, all conformers may reasonably be populated at STP in the gas phase. We also report the harmonic vibrational frequencies of conformers, including the implications on the mechanism of detonation. In particular, we establish that the weakest N-N nitramine of CL-20 is the cyclohexane equatorial nitramine. This preliminary mapping of the conformers of CL-20 makes it possible to study the mechanism of detonation of this explosive rigorously in future work.

The great diversity of HMX conformers: probing the potential energy surface using CCSD(T).

The octahydro-1,3,5,7-tetranitro-1,3,5,7-tetraazocine (HMX) molecule is a very commonly studied system, in all 3 phases, because of its importance as an explosive; however, no one has ever attempted a systematic study of what all the major gas-phase conformers are. This is critical to a mechanistic study of the kinetics involved, as well as the viability of various crystalline polymorphs based on the gas-phase conformers. We have used existing knowledge of basic cyclooctane chemistry to survey all possible HMX conformers based on its fundamental ring structure. After studying what geometries are possible after second-order many-body perturbation theory (MBPT(2)) geometry optimization, we calculated the energetics using coupled cluster singles, doubles, and perturbative triples (CCSD(T))/cc-pVTZ. These highly accurate energies allow us to better calculate starting points for future mechanistic studies. Additionally, the plethora of structures are compared to existing experimental data of crystals. It is found that the crystal field effect is sometimes large and sometimes small for HMX.

Pragmatic ab initio prediction of enthalpies of formation for large molecules: accuracy of MP2 geometries and frequencies using CCSD(T) correlation energies

We have addressed the accuracy of calculating the enthalpy of formation of an arbitrary single reference molecule using practical ab initio methodologies. It is known that MP2 geometries with a triple zeta basis set are almost as reliable as CCSD(T) geometries. It is also known that CCSD(T) correlation energies, with basis extrapolation, feature chemical accuracy for single-reference molecules. We investigate what accuracy one might expect in enthalpies of formation from a MP2 geometry, MP2 harmonic vibrational frequencies, a CCSD(T) correlation energy using triple zeta basis sets. It is far from obvious, a priori, as to which error source contributes most significantly. We observe that the accuracy in calculating enthalpies of formation of single-reference molecules with this protocol is 4 kcal mol-1; our error analysis shows this comes almost exclusively from the correlation energy basis extrapolation, rather than errors intrinsic to MP2.