Beyond the benzene dimer: an investigation of the additivity of pi-pi interactions

The benzene dimer is the simplest prototype of pi-pi interactions and has been used to understand the fundamental physics of these interactions as they are observed in more complex systems. In biological systems, however, aromatic rings are rarely found in isolated pairs; thus, it is important to understand whether aromatic pairs remain a good model of pi-pi interactions in clusters. In this study, ab initio methods are used to compute the binding energies of several benzene trimers and tetramers, most of them in 1D stacked configurations. The two-body terms change only slightly relative to the dimer, and except for the cyclic trimer, the three- and four-body terms are negligible. This indicates that aromatic clusters do not feature any large nonadditive effects in their binding energies, and polarization effects in benzene clusters do not greatly change the binding that would be anticipated from unperturbed benzene-benzene interactions, at least for the 1D stacked systems considered. Three-body effects are larger for the cyclic trimer, but for all systems considered, the computed binding energies are within 10% of what would be estimated from benzene dimer energies at the same geometries.     

High accuracy ab initio studies of Li6+, Li6-, and three isomers of Li6

The structures and energetics of Li(6) (+), Li(6) (-) and three isomers of Li(6) are investigated using the coupled-cluster singles, doubles and perturbative triples [CCSD(T)] method with valence and core-valence correlation consistent basis sets of double- to quadruple-zeta quality (cc-pVXZ and cc-pCVXZ, where X=D-Q). These results are compared with qualitatively different predictions by less reliable methods. Our results conclusively show that the D(4h) isomer is the global minimum structure for Li(6). It is energetically favored over the C(5v) and D(3h) structures by about 5.1 and 7.1 kcal mol(-1), respectively, after the inclusion of the zero-point vibrational energy (ZPVE) correction. Our most accurate total atomization energies are 123.2, 117.6, and 115.7 kcal mol(-1) for the D(4h), C(5v), and D(3h) isomers, respectively. Comparison of experimental optical absorption spectra with our computed electronic spectra also indicate that the D(4h) isomer is indeed the most stable structure. The cation, anion, and some higher spin states are investigated using the less expensive cc-pCVDZ basis set. Adiabatic ionization energies and electron affinities are reported and compared with experimental values. Predictions of molecular properties are found to be sensitive to the basis set used and to the treatment of electron correlation.     

Nucleic acid analysis using an expanded genetic alphabet to quench fluorescence

Organic chemistry has made possible the synthesis of molecules that expand on Nature's genetic alphabet. Using the previously described nonstandard DNA base pair constructed from isoguanine and 5-methylisocytosine, we report a highly specific and sensitive method that allows for the fast and specific quantitation of genetic sequences in a closed tube format. During PCR amplification, enzymatic site-specific incorporation of a quencher covalently linked to isoguanine allows for the simultaneous detection and identification of multiple targets. The specificity of method is then established by analysis of thermal denaturation or melting of the amplicons. The appropriate functions of all reactions are further verified by incorporation of an independent target into the reaction mixture. We report that the method is sensitive down to the single copy level, and specificity is demonstrated by multiplexed end-point genotypic analysis of four targets simultaneously using four separate fluorescent reporters. The method is general enough for quantitative and qualitative analysis of both RNA and DNA using previously developed primer sets. Though the method described employs the commonly used PCR, the enzymatic incorporation of reporter groups into DNA site-specifically should find broad utility throughout molecular biology.     

Systematic Study of Selected Diagonalization Methods for Configuration Interaction Matrices

Several modifications to the Davidson algorithm are systematically explored to establish their performance for an assortment of configuration interaction (CI) computations. The combination of a generalized Davidson method, a periodic two‐vector subspace collapse, and a blocked Davidson approach for multiple roots is determined to retain the convergence characteristics of the full subspace method. This approach permits the efficient computation of wave functions for large‐scale CI matrices by eliminating the need to ever store more than three expansion vectors ( b _i) and associated matrix‐vector products ( σ _i), thereby dramatically reducing the I/O requirements relative to the full subspace scheme. The minimal‐storage, single‐vector method of Olsen is found to be a reasonable alternative for obtaining energies of well‐behaved systems to within μE_h accuracy, although it typically requires around 50% more iterations and at times is too inefficient to yield high accuracy (ca. 10^?10 E_h) for very large CI problems. Several approximations to the diagonal elements of the CI Hamiltonian matrix are found to allow simple on‐the‐fly computation of the preconditioning matrix, to maintain the spin symmetry of the determinant‐based wave function, and to preserve the convergence characteristics of the diagonalization procedure. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1574–1589, 2001     

Reduced occupancy of the deeply bound 0d(5/2) neutron state in 32Ar

The 9Be(32Ar, 31Ar)X reaction, leading to the 5/2+ ground state of a nucleus at the proton drip line, has a cross section of 10.4(13) mb at a beam energy of 65.1 MeV/nucleon. This translates into a spectroscopic factor that is only 24(3)% of that predicted by the many-body shell-model theory. We introduce refinements to the eikonal reaction theory used to extract the spectroscopic factor to clarify that this very strong reduction represents an effect of nuclear structure. We suggest that it reflects correlation effects linked to the high neutron separation energy (22.0 MeV) for this state.     

Direct evidence for the breakdown of the N = 8 shell closure in 12Be

We derive the formalism to obtain spatial distributions of collisional correlation times for macroscopic particles undergoing granular flow from pulsed gradient spin echo nuclear magnetic resonance diffusion data. This is demonstrated with an example of axial motion in the shear flow regime of a 3D granular flow in a horizontal rotating cylinder at one rotation rate.     

Assessing the performance of density functional theory for the electronic structure of metal-salens: the M06 suite of functionals and the d⁴-metals

We have systematically investigated the electronic structure of the d? metal-salen complexes including the Cr(II)-, Mn(III)-, Fe(IV)-, Mo(II)-, Tc(III)-, and Ru(IV)-salen complexes. Density functional theory (DFT) has been employed, using the BP86 and B3LYP functionals, and the entire M05 and M06 suites of meta-generalized gradient functionals. These results are compared to robust complete active-space self-consistent field (CASSCF) optimized geometries and complete active-space third-order perturbation theory (CASPT3) energies for the lowest singlet, triplet, and quintet states. Although the M06 and M06-L DFT functionals have been generally recommended for computations on complexes that contain main group and transition metals, none of the M0-functionals appear statistically better than the B3LYP functional for the computation of spin-state energy gaps. DFT- and CASSCF-optimized geometries normally agree to within 0.3 ? least root mean squared deviation (LRMSD) for the singlet and triplet structures and less than 0.1 ? LRMSD for the quintet structures. It can be concluded that no DFT functional tested here can be considered reliable over all metal-salen complexes and it is highly recommended that the accuracy of any given DFT functional should be assessed on a case-by-case basis.     

Origin of the surprising enhancement of electrostatic energies by electron-donating substituents in substituted sandwich benzene dimers

A recent study of substituted face-to-face benzene dimers by Lewis and co-workers [J. Am. Chem. Soc. 2011, 133, 3854-3862] indicated a surprising enhancement of electrostatic interactions for both electron-withdrawing and electron-donating substituents. Here we demonstrate that charge penetration (an attractive electrostatic interaction arising from the overlap of the electron densities on the two monomers) is the cause of this counterintuitive effect. These charge penetration effects are significant at typical π-π interaction distances, and they are not easily described by multipole models. A simple measure of a substituent's electron-donating or electron-withdrawing character, such as the Hammett parameter σ(m), is unlikely to capture subtle charge penetration effects. Indeed, correlation of the relative total energies or relative electrostatic energies with ∑σ(m) breaks down for multiply substituted face-to-face benzene dimers.