Ab Initio Characterization of C₂H₄N₂ Isomers: Structures, electronic energies, spectroscopic parameters, and formation pathways

This work presents a comprehensive theoretical investigation of key isomers of C2H4N2 using state-of-the-art quantum chemical methods. The objective is to characterize their molecular structures, spectroscopic constants, and electronic energies, and elucidate plausible formation and destruction pathways, providing data critical for astrochemical and atmospheric detection. High-accuracy ab initio methods were employed, notably CCSD(T)-F12/cc-pVTZ-F12 for optimized geometries. Additional calculations were performed at the CCSD(T)/aug-cc-pVTZ, CCSD(T)/cc-pVTZ,MP2/aug-cc-pVTZ, and CIS levels. Intrinsic reaction coordinate (IRC) calculations were performed at the B3LYP/6-31G(d, p) level to explore reaction pathways. The ZPE-corrections were determined for all of the isomers considered. Six low-energy C2H4N2 isomers were identified, all within 1 eV of the global minimum. Among them, methyl-cyanamide (MCA) exhibits the lowest relative energy (∼0.2 eV) and a significant electric dipole moment of 5.00 D, making it a strong candidate for detection in the gas phase environments. The rotational constants for MCA, computed at the level of CCSD(T)-F12/cc-pVTZ-F12, are Ae = 34,932.44 MHz, Be = 4,995.31 MHz and Ce = 4,520.30 MHz. The V3 torsional barrier was found to be 631.19 cm−1. Centrifugal distortion constants were computed up to sextic order for all isomers. Formation pathways for MCA such as CH3N + HCN →CH3NHCN and related isomers were characterized. The combination of large dipole moments and distinct rotational signatures supports the detectability of methylcyanamide and related C2H4N2 isomers via radioastronomy, IR and MW spectroscopy. Isomerization and reaction pathways processes were found to be key to their formation in the gas phase environments. These results offer a foundation for observational efforts.