Search for: All records

Heteroaromatic species are commonly found in complex gaseous mixtures, from tobacco smoke to petroleum and asphaltene combustion products. At high temperatures, C–H bond rupture produces various dehydro radical isomers. We have used the spin–flip formulation of equation-of-motion coupled cluster theory with single and double substitutions (EOM-SF-CCSD) to characterize the energies and wave functions of the lowest lying singlet and triplet states of the diradical (2,3), (2,4), (2,5), and (3,4) di-dehydro isomers of pyrrole, furan, and thiophene. In all cases, these diradicals are minima on the broken-symmetry ωB97X-D/cc-pVDZ potential energy surface. In most cases, the diradical geometries distort to enhance through-space or through-bond coupling in the singlet states and to avoid Coulombic or exchange repulsion in the triplet states. EOM-SF-CCSD results indicate that all diradical isomers are two-configurational, closed shell singlet states. The only exceptions to this are for (2,3) and (2,4) thiophene and (2,3) pyrrole, which each contain more than two configurations. In all cases, the leading term in the multiconfigurational diradical wave function doubly occupies the symmetric radical σ orbital, indicative of either through-space or 1,3 through-bond coupling. We utilized the nucleus-independent chemical shift (NICS) approach to qualitatively assess aromaticity and find that this property varies and may be related to the energetic splittings in these diradical isomers. 

High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π– π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller–Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π– π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol −1 , while significant changes of up to 3 kcal mol −1 of the electrostatic component are seen for the cationic pyridinium–benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little “tuning” of the direct solute–solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute–solute and solute–solvent interactions.