Nucleophilic Aromatic Addition in Ionizing Environments: Observation and Analysis of New C-N Valence Bonds in Complexes between Naphthalene Radical Cation and Pyridine.

Radical organic ions can be stabilized by complexation with neutral organics via interactions that can resemble chemical bonds, but with much diminished bond energies. Those interactions are a key factor in cluster growth and polymerization reactions in ionizing environments such as regions of the interstellar medium and solar nebulae. Such radical cation complexes between naphthalene (Naph) and pyridine (Pyr) are characterized using mass-selected ion mobility experiments. The measured enthalpy of binding of the Naph+•(Pyr) heterodimer (20.9 kcal/mol) exceeds that of the Naph+•(Naph) homodimer (17.8 kcal/mol). The addition of 1-3 more pyridine molecules to the Naph+•(Pyr) heterodimer gives 10-11 kcal/mol increments in binding enthalpy. A rich array of Naph+•(Pyr) isomers are characterized by electronic structure calculations. The calculated Boltzmann distribution at 400 K yields an enthalpy of binding in reasonable agreement with experiment. The global minimum is a distonic cation formed by Pyr attack on Naph+• at the α-carbon, changing its hybridization from sp2 to distorted sp3. The measured collision cross section in helium for the Naph+•(Pyr) heterodimer of 84.9 ± 2.5 Å2 at 302 K agrees well with calculated angle-averaged cross sections (83.9-85.1 Å2 at 302 K) of the lowest energy distonic structures. A remarkable 16 kcal/mol increase in the binding energy between Naph+•(Pyr) and Bz+•(Pyr) (Bz is benzene) is understood by energy decomposition analysis. A similar increase in binding from Naph+•(NH3) to Naph+•(Pyr) (as well as between Bz+•(NH3) and Bz+•(Pyr)) is likewise rationalized.

Gas phase hydration of halogenated benzene cations. Is it hydrogen or halogen bonding?

Halogen bonding (XB) non-covalent interactions can be observed in compounds containing chlorine, bromine, or iodine which can form directed close contacts of the type R1-XY-R2, where the halogen X acts as a Lewis acid and Y can be any electron donor moiety including electron lone pairs on hetero atoms such as O and N, or π electrons in olefin double bonds and aromatic conjugated systems. In this work, we present the first evidence for the formation of ionic halogen bonds (IXBs) in the hydration of bromobenzene and iodobenzene radical cations in the gas phase. We present a combined thermochemical investigation using the mass-selected ion mobility (MSIM) technique and density functional theory (DFT) calculations of the stepwise hydration of the fluoro, chloro, bromo, and iodobenzene radical cations. The binding energy associated with the formation of an IXB in the hydration of the iodobenzene cation (11.2 kcal mol-1) is about 20% higher than the typical unconventional ionic hydrogen bond (IHB) of the CHδ+OH2 interaction. The formation of an IXB in the hydration of the iodobenzene cation involves a significant entropy loss (29 cal mol-1 K-1) resulting from the formation of a more ordered structure and a highly directional interaction between the oxygen lone pair of electrons of water and the electropositive region around the iodine atom of the iodobenzene cation. In comparison, the hydration of the fluorobenzene and chlorobenzene cations where IHBs are formed, -ΔS° = 18-21 cal mol-1 K-1 consistent with the formation of less ordered structures and loose interactions. The electrostatic potentials on the lowest energy structures of the hydrated halogenated benzene radical cations show clearly that the formation of an IXB is driven by a positively charged σ-hole on the external side of the halogen atom X along the C-X bond axis. The size of the σ-hole increases significantly in bromobenzene and iodobenzene radical cations which results in strong interaction potentials with the electron lone pairs of the oxygen atom of the water molecules and thus IXBs provide the most stable hydrated structures of the bromobenzene and iodobenzene radical cations. The results clearly distinguish the hydration behaviors resulting from the ionic hydrogen and halogen bonding interactions of fluorobenzene and iodobenzene cations, respectively, and establish the different bonding and structural features of the two interactions.

Observation of covalent and electrostatic bonds in nitrogen-containing polycyclic ions formed by gas phase reactions of the benzene radical cation with pyrimidine.

Polycyclic aromatic hydrocarbons (PAHs) and polycyclic aromatic nitrogen heterocyclics (PANHs) are present in ionizing environments, including interstellar clouds and solar nebulae, where their ions can interact with neutral PAH and PANH molecules leading to the formation of a variety of complex organics including large N-containing ions. Herein, we report on the formation of a covalently-bonded (benzene·pyrimidine) radical cation dimer by the gas phase reaction of pyrimidine with the benzene radical cation at room temperature using the mass-selected ion mobility technique. No ligand exchange reactions with benzene and pyrimidine are observed indicating that the binding energy of the (benzene·pyrimidine)˙+ adduct is significantly higher than both the benzene dimer cation and the proton-bound pyrimidine dimer. The (benzene·pyrimidine)˙+ adduct shows thermal stability up to 541 K. Thermal dissociation of the (C6D6·C4H4N2)˙+ adduct at temperatures higher than 500 K produces C4H4N2D+ (m/z 82) suggesting the transfer of a D atom from the C6D6 moiety to the C4H4N2 moiety before the dissociation of the adduct. Mass-selected ion mobility of the (benzene·pyrimidine)˙+ dimer reveals the presence of two families of isomers formed by electron impact ionization of the neutral (benzene·pyrimidine) dimer. The slower mobility peak corresponds to a non-covalent family of isomers with larger collision cross sections (76.0 ± 1.8 Å2) and the faster peak is consistent with a family of covalent isomers with more compact structures and smaller collision cross sections (67.7 ± 2.2 Å2). The mobility measurements at 509 K show only one peak corresponding to the family of stable covalently bonded isomers characterized by smaller collision cross sections (66.9 ± 1.9 Å2 at 509 K). DFT calculations at the M06-2X/6-311++G** level show that the most stable (benzene·pyrimidine)˙+ isomer forms a covalent C-N bond with a binding energy of 49.7 kcal mol-1 and a calculated collision cross section of 69.2 Å2, in excellent agreement with the value obtained from the faster mobility peak of the (benzene·pyrimidine)˙+ dimer. Formation of a C-N covalent bond displaces a hydrogen atom from a C-H bond of the benzene cation which is transferred to the second pyrimidine nitrogen atom, thus preserving the pyrimidine π system and yielding the most stable (benzene·pyrimidine)˙+ isomer. The calculations also show less stable non-covalent electrostatically bonded perpendicular isomers of the (benzene·pyrimidine)˙+ dimer with a binding energy of 19 kcal mol-1 and a calculated collision cross section of 74.0-75.0 Å2 in excellent agreement with the value obtained from the slower mobility peak of the (benzene·pyrimidine)˙+ dimer.

Unconventional CH(δ+)···N hydrogen bonding interactions in the stepwise solvation of the naphthalene radical cation by hydrogen cyanide and acetonitrile molecules.

Equilibrium thermochemical measurements using the mass-selected ion mobility (MSIM) technique have been utilized to investigate the binding energies and entropy changes of the stepwise association of hydrogen cyanide (HCN) and acetonitrile (CH3CN) molecules with the naphthalene radical cation (C10H8˙(+)) in the gas phase forming the C10H8˙(+)(HCN)n and C10H8˙(+)(CH3CN)n clusters with n = 1-3 and 1-5, respectively. The lowest energy structures of the C10H8˙(+)(HCN)n and C10H8˙(+)(CH3CN)n clusters for n = 1-2 have been calculated using the M062X and ω97XD methods within the 6-311+G** basis set, and for n = 1-6 using the B3LYP method within the 6-311++G** basis set. In both systems, the initial interaction occurs through unconventional CH(δ+)···N ionic hydrogen bonds between the hydrogen atoms of the naphthalene cation and the lone pair of electrons on the N atom of the HCN or the CH3CN molecule. The binding energy of CH3CN to the naphthalene cation (11 kcal mol(-1)) is larger than that of HCN (7 kcal mol(-1)) due to a stronger ion-dipole interaction resulting from the large dipole moment of CH3CN (3.9 D). On the other hand, HCN can form both unconventional hydrogen bonds with the hydrogen atoms of the naphthalene cation (CH(δ+)···NCH), and conventional linear hydrogen bonding chains involving HCN···HCN interactions among the associated HCN molecules. HCN molecules tend to form "externally solvated" structures with the naphthalene cation where the naphthalene ion is hydrogen bonded to the exterior of an HCN···HCN chain. For the C10H8˙(+)(CH3CN)n clusters, "internally solvated" structures are favored where the acetonitrile molecules are directly interacting with the naphthalene cation through CH(δ+)···N unconventional ionic hydrogen bonds. In both the C10H8˙(+)(HCN)n and C10H8˙(+)(CH3CN)n clusters, the sequential binding energy decreases stepwise to about 6-7 kcal mol(-1) by three HCN or CH3CN molecules, approaching the macroscopic enthalpy of vaporization of liquid HCN (6.0 kcal mol(-1)).

Communication: Ion mobility of the radical cation dimers: (Naphthalene)2(+•) and naphthalene(+•)-benzene: Evidence for stacked sandwich and T-shape structures.

Dimer radical cations of aromatic and polycyclic aromatic molecules are good model systems for a fundamental understanding of photoconductivity and ferromagnetism in organic materials which depend on the degree of charge delocalization. The structures of the dimer radical cations are difficult to determine theoretically since the potential energy surface is often very flat with multiple shallow minima representing two major classes of isomers adopting the stacked parallel or the T-shape structure. We present experimental results, based on mass-selected ion mobility measurements, on the gas phase structures of the naphthalene(+⋅) ⋅ naphthalene homodimer and the naphthalene(+⋅) ⋅ benzene heterodimer radical cations at different temperatures. Ion mobility studies reveal a persistence of the stacked parallel structure of the naphthalene(+⋅) ⋅ naphthalene homodimer in the temperature range 230-300 K. On the other hand, the results reveal that the naphthalene(+⋅) ⋅ benzene heterodimer is able to exhibit both the stacked parallel and T-shape structural isomers depending on the experimental conditions. Exploitation of the unique structural motifs among charged homo- and heteroaromatic-aromatic interactions may lead to new opportunities for molecular design and recognition involving charged aromatic systems.

Proton-bound dimers of nitrogen heterocyclic molecules: substituent effects on the structures and binding energies of homodimers of diazine, triazine, and fluoropyridine.

The bonding energies of proton-bound homodimers BH(+)B were measured by ion mobility equilibrium studies and calculated at the DFT B3LYP/6-311++G** level, for a series of nitrogen heterocyclic molecules (B) with electron-withdrawing in-ring N and on-ring F substituents. The binding energies (ΔH°(dissoc)) of the proton-bound dimers (BH(+)B) vary significantly, from 29.7 to 18.1 kcal/mol, decreasing linearly with decreasing the proton affinity of the monomer (B). This trend differs significantly from the constant binding energies of most homodimers of other organic nitrogen and oxygen bases. The experimentally measured ΔH°(dissoc) for (1,3-diazine)2H(+), i.e., (pyrimidine)2H(+) and (3-F-pyridine)2H(+) are 22.7 and 23.0 kcal/mol, respectively. The measured ΔH°(dissoc) for the pyrimidine(·+)(3-F-pyridine) radical cation dimer (19.2 kcal/mol) is signifcantly lower than that of the proton-bound homodimers of pyrimidine and 3-F-pyridine, reflecting the stronger interaction in the ionic H-bond of the protonated dimers. The calculated binding energies for (1,2-diazine)2H(+), (pyridine)2H(+), (2-F-pyridine)2H(+), (3-F-pyridine)2H(+), (2,6-di-F-pyridine)2H(+), (4-F-pyridine)2H(+), (1,3-diazine)2H(+), (1,4-diazine)2H(+), (1,3,5-triazine)2H(+), and (pentafluoropyridine)2H(+) are 29.7, 24.9, 24.8, 23.3, 23.2, 23.0, 22.4, 21.9, 19.3, and 18.1 kcal/mol, respectively. The electron-withdrawing substituents form internal dipoles whose electrostatic interactions contribute to both the decreased proton affinities of (B) and the decreased binding energies of the protonated dimers BH(+)B. The bonding energies also vary with rotation about the hydrogen bond, and they decrease in rotamers where the internal dipoles of the components are aligned efficiently for inter-ring repulsion. For compounds substituted at the 3 or 4 (meta or para) positions, the lowest energy rotamers are T-shaped with the planes of the two rings rotated by 90° about the hydrogen bond, while the planar rotamers are weakened by repulsion between the ortho hydrogen atoms of the two rings. Conversely, in ortho-substituted (1,2-diazine)2H(+) and (2-F-pyridine)2H(+), attractive interactions between the ortho (C-H) hydrogen atoms of one ring and the electronegative ortho atoms (N or F) of the other ring are stabilizing, and increase the protonated dimer binding energies by up to 4 kcal/mol. In all of the dimers, rotation about the hydrogen bond can involve a 2-4 kcal/mol barrier due to the relative energies of the rotamers.

What Is the Structure of the Naphthalene-Benzene Heterodimer Radical Cation? Binding Energy, Charge Delocalization, and Unexpected Charge-Transfer Interaction in Stacked Dimer and Trimer Radical Cations.

The binding energy of the naphthalene(+•)(benzene) heterodimer cation has been determined to be 7.9 ± 1 kcal/mol for C10H8(+•)(C6H6) and 8.1 ± 1 kcal/mol for C10H8(+•)(C6D6) by equilibrium thermochemical measurements using the mass-selected drift cell technique. A second benzene molecule binds to the C10H8(+•)(C6D6) dimer with essentially the same energy (8.4 ± 1 kcal/mol), suggesting that the two benzene molecules are stacked on opposite sides of the naphthalene cation in the (C6D6)C10H8(+•)(C6D6) heterotrimer. The lowest-energy isomers of the C10H8(+•)(C6D6) and (C6D6)C10H8(+•)(C6D6) dimer and trimer calculated using the M11/cc-pVTZ method have parallel stacked structures with enthalpies of binding (-ΔH°) of 8.4 and 9.0 kcal/mol, respectively, in excellent agreement with the experimental values. The stacked face-to-face class of isomers is calculated to have substantial charge-transfer stabilization of about 45% of the total interaction energy despite the large difference between the ionization energies of benzene and naphthalene. Similarly, significant delocalization of the positive charge is found among all three fragments of the (C6D6)C10H8(+•)(C6D6) heterotrimer, thus leaving only 46% of the total charge on the central naphthalene moiety. This unexpectedly high charge-transfer component results in activating two benzene molecules in the naphthalene(+•)(benzene)2 heterotrimer cation to associate with a third benzene molecule at 219 K to form a benzene trimer cation and a neutral naphthalene molecule. The global minimum of the C10H8(+•)(C6H6)2 heterotrimer is found to be the one where the naphthalene cation is sandwiched between two benzene molecules. It is remarkable, and rather unusual, that the binding energy of the second benzene molecule is essentially the same as that of the first. This is attributed to the enhanced charge-transfer interaction in the stacked trimer radical cation.