skip to main content


Title: Proton relays in anomalous carbocations dictate spectroscopy, stability, and mechanisms: case studies on C 2 H 5 + and C 3 H 3 +
We present a detailed analysis of the anomalous carbocations: C 2 H 5 + and C 3 H 3 + . This work involves (a) probing electronic structural properties, (b) ab initio dynamics simulations over a range of internal energies, (c) analysis of reduced dimensional potential surfaces directed along selected conformational transition pathways, (d) dynamically averaged vibrational spectra computed from ab initio dynamics trajectories, and (e) two-dimensional time–frequency analysis to probe conformational dynamics. Key findings are as follows: (i) as noted in our previous study on C 2 H 3 + , it appears that these non-classical carbocations are stabilized by delocalized nuclear frameworks and “proton shuttles”. We analyze this nuclear delocalization and find critical parallels between conformational changes in C 2 H 3 + , C 2 H 5 + , and C 3 H 3 + . (ii) The vibrational signatures of C 2 H 5 + are dominated by the “bridge-proton” conformation, but also show critical contributions from the “classical” configuration, which is a transition state at almost all levels of theory. This result is further substantiated through two-dimensional time–frequency analysis and is at odds with earlier explanations of the experimental spectra, where frequencies close to the classical region were thought to arise from an impurity. While this is still possible, our results here indicate an additional (perhaps more likely) explanation that involves the “classical” isomer. (iii) Finally, in the case of C 3 H 3 + our explanation of the experimental result includes the presence of multiple, namely, “cyclic”, “straight”, and propargyl, configurations. Proton shuttles and nuclear delocalization, reminiscent of those seen in the case of C 2 H 3 + , were seen all through and have a critical role in all our observations.  more » « less
Award ID(s):
1665336
NSF-PAR ID:
10065845
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Phys. Chem. Chem. Phys.
Volume:
19
Issue:
40
ISSN:
1463-9076
Page Range / eLocation ID:
27801 to 27816
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A remarkable distinction between boron and carbon hydrides lies in their extremely different bonding patterns and chemical reactivity, resulting in diverse areas of application. Particularly, carbon, characterized by classical two‐center – two‐electron bonds, gives rise to organic chemistry. In contrast, boron forms numerous exotic and non‐intuitive compounds collectively called non‐classical structures. It is reasonable to anticipate that other elements of Group 13 exhibit their own unusual bonding patterns; however, our knowledge of the hydride chemistry for other elements in Group 13 is much more limited, especially for the heaviest stable element, thallium. In this work, we performed a conformational analysis of Tl2Hxand Tl3Hy(x=0–6, y=0–5) series via Coalescence Kick global minimum search algorithm, DFT, andab initioquantum chemistry methods; we investigated the bonding pattern using the AdNDP algorithm, thermodynamic stability, and stability toward electron detachment. All found global minimum structures are classified as non‐classical structures featuring at least one multi‐center bond.

     
    more » « less
  2. Abstract

    Organofunctionalized tetranuclear clusters [(MIICl)2(VIVO)2{((HOCH2CH2)(H)N(CH2CH2O))(HN(CH2CH2O)2)}2] (1, M=Co,2: M=Zn) containing an unprecedented oxometallacyclic {M2V2Cl2N4O8} (M=Co, Zn) framework have been prepared by solvothermal reactions. The new oxo‐alkoxide compounds were fully characterized by spectroscopic methods, magnetic susceptibility measurement, DFT and ab initio computational methods, and complete single‐crystal X‐ray diffraction structure analysis. The isostructural clusters are formed of edge‐sharing octahedral {VO5N} and trigonal bipyramidal {MO3NCl} units. Diethanolamine ligates the bimetallic lacunary double cubane core of1and2in an unusual two‐mode fashion, unobserved previously. In the crystalline state, the clusters of1and2are joined by hydrogen bonds to form a three‐dimensional network structure. Magnetic susceptibility data indicate weakly antiferromagnetic interactions between the vanadium centers [Jiso(VIV−VIV)=−5.4(1); −3.9(2) cm−1], and inequivalent antiferromagnetic interactions between the cobalt and vanadium centers [Jiso(VIV−CoII)=−12.6 and −7.5 cm−1] contained in1.

     
    more » « less
  3. The gas-phase reaction of the methylidyne (CH; X 2 Π) radical with dimethylacetylene (CH 3 CCCH 3 ; X 1 A 1g ) was studied at a collision energy of 20.6 kJ mol −1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. The crossed molecular beam experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C 5 H 7 reaction intermediate(s) ultimately dissociating to C 5 H 6 isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylenecyclopropene, 13%) and p8 (1-penten-3-yne, 81%) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products. The dynamics are controlled by addition to the carbon–carbon triple bond with the reaction intermediates eventually eliminating a hydrogen atom from the methyl groups of the dimethylacetylene reactant forming 1-methyl-3-methylenecyclopropene (p28). The dominating pathways reveal an unexpected insertion of methylidyne into one of the six carbon–hydrogen single bonds of the methyl groups of dimethylacetylene leading to the acyclic intermediate, which then decomposes to 1-penten-3-yne (p8). Therefore, the methyl groups of dimethylacetylene effectively ‘screen’ the carbon–carbon triple bond from being attacked by addition thus directing the dynamics to an insertion process as seen exclusively in the reaction of methylidyne with ethane (C 2 H 6 ) forming propylene (CH 3 C 2 H 3 ). Therefore, driven by the screening of the triple bond, one propynyl moiety (CH 3 CC) acts in four out of five trajectories as a spectator thus driving an unexpected, but dominating chemistry in analogy to the methylidyne – ethane system. 
    more » « less
  4. null (Ed.)
    Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a “ring polymer” whose spread is related to the quantum mechanical uncertainty in the particle’s position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions. 
    more » « less
  5. null (Ed.)
    We investigate the excited electron dynamics in [Au 25 (SR) 18 ] −1 (R = CH 3 , C 2 H 5 , C 3 H 7 , MPA, PET) [MPA = mercaptopropanoic acid, PET = phenylethylthiol] nanoparticles to understand how different ligands affect the excited state dynamics in this system. The population dynamics of the core and higher excited states lying in the energy range 0.00–2.20 eV are studied using a surface hopping method with decoherence correction in a real-time DFT approach. All of the ligated clusters follow a similar trend in decay for the core states (S 1 –S 6 ). The observed time constants are on the picosecond time scale (2–19 ps), which agrees with the experimental time scale, and this study confirms that the time constants observed experimentally could originate from core-to-core transitions and not from core-to-semiring transitions. In the presence of higher excited states, R = H, CH 3 , C 2 H 5 , C 3 H 7 , and PET demonstrate similar relaxations trends whereas R = MPA shows slightly different relaxation of the core states due to a smaller gap between the LUMO+1 and LUMO+2 gap in its electronic structure. The S 1 (HOMO → LUMO) state gives the slowest decay in all ligated clusters, while S 7 has a relatively long decay. Furthermore, separate electron and hole relaxations were performed on the [Au 25 (SCH 3 ) 18 ] −1 nanocluster to understand how independent electron and hole relaxations contribute to the overall relaxation dynamics. 
    more » « less