Inorganic polysulfides: Quantum chemistry study and biological implications

Mechanisms of the Gewald Synthesis of 2-Aminothiophenes from Elemental Sulfur

The Gewald reaction is a well-established one-pot method to access 2-aminothiophenes from carbonyl compounds, activated acetonitriles, and elemental sulfur. To elucidate the reaction's poorly understood mechanism, with regard to the decomposition of sulfur and polysulfide intermediates, we have performed a comprehensive computational study using density functional theory (DFT) calculations at the M06-2X (or ωB97X-D)/aug-cc-pV(T + d)Z/SMD(C2H5OH) level of theory. The results show that the reaction is initiated by a Knoevenagel-Cope condensation, followed by opening of the elemental sulfur, leading to polysulfide formation. The polysulfide intermediates can interconvert and decompose using various mechanisms including unimolecular cyclization, nucleophilic degradation, and scrambling. Protonation of the polysulfides changes their electrophilic behavior and provides a kinetically favorable pathway for their decomposition. This protonation-induced intermolecular degradation is feasible for polysulfides of all lengths, but unimolecular decomposition is kinetically favored for long polysulfides (≥6 sulfur atoms). None of the pathways provide any thermodynamic benefit due to the lack of resonance-stabilized leaving group, and a complex equilibrium of polysulfides of all lengths is expected in solution. Cyclization of the monosulfide with aromatization to the thiophene product is the only driving force behind the reaction, funneling all of the various intermediates into the observed product in a thermodynamically controlled process.

Benchmark of density functional theory methods for the study of organic polysulfides

Elemental sulfur is often used in organic synthesis as its low cost and high abundance make it a highly desirable source of sulfur atoms. However, sulfur's unpredictable catenation behavior poses challenges to its widespread usage due to difficulties in designing new reactions that can account for its multifaceted reactivity. In order to accurately model sulfur's mechanisms using computational approaches, it is necessary to identify density functional theory (DFT) methods that are accurate on these systems. This study benchmarks 12 well-known DFT functionals that include local, non-local, and hybrid methods against DLPNO-CCSD(T)/aug-cc-pV(Q+d)Z//MP2/aug-cc-pV(T+d)Z/SMD(MeCN) for the accurate treatment of organic polysulfides, taking cyanide as a nucleophile. Our benchmarking results indicate that the M06-2X and B3LYP-D3(BJ) density functionals are the most accurate for calculating reaction energies, while local functionals performed the worst. For activation energies, MN15, MN15-L, M06-2X, and ωB97X-D are the most accurate. Our analysis of structural parameters shows that all functionals perform well for ground state optimizations except B97D3, while MN15-L and M06-2X performed best for transition structure optimizations. Overall, the four hybrid functionals MN15, M06-2X, ωB97X-D, and B3LYP-D3(BJ) appear adequate for studying the reaction mechanisms of polysulfides.