Advances in Determining the Absolute Proton Affinities of Neutral Organic Molecules in the Gas Phase and Their Interpretation: A Theoretical Account

A Systematic Study of Nonionic Di- and Multiborane Catalysts for the Oligomerization and Polymerization of Epoxides

Borane catalysis has emerged as a powerful technology in epoxide polymerization. Still, the structure-activity correlations for these catalysts are not fully understood to date, especially regarding compounds with nonionic backbones. Thus, in this work, 13 different borane catalysts of this respective type are described and investigated for their epoxide oligomerization and polymerization performance, using propylene oxide (PO), 1-butylene oxide (BO) and allyl glycidyl ether (AGE) as monomers. Structurally, special emphasis is put on catalysts with different linker lengths and linker flexibilities as well as the introduction of more than two borane functionalities. Importantly, this screening is conducted both under typical polymerization conditions as well as under the chain transfer agent (CTA)-rich conditions relevant for large-scale production. It is found that suitable preorganization of the borane groups, such as present in biphenyl derivatives, offers a simple route to high-performing catalysts and quantitative monomer conversion of the investigated epoxides. Furthermore, it is demonstrated that a diborane-catalyzed oligomerization can be kept active over weeks, whereby repeated addition of monomer batches (14 steps) constantly results in full conversion and well-defined oligoethers, underlining the practical potential of this method. The absence of co-initiating counter ions is suggested as an inherent advantage of nonionic catalysts.

Unraveling redox pathways of the disulfide bond in dimethyl disulfide: Ab initio modeling

CONTEXT

In cellular environments, the reduction of disulfide bonds is pivotal for protein folding and synthesis. However, the intricate enzymatic mechanisms governing this process remain poorly understood. This study addresses this gap by investigating a disulfide bridge reduction reaction, serving as a model for comprehending electron and proton transfer in biological systems. Six potential mechanisms for reducing the dimethyl disulfide (DMDS) bridge through electron and proton capture were explored. Thermodynamic and kinetic analyses elucidated the sequence of proton and electron addition. MD-PMM, a method that combines molecular dynamics simulations and quantum-chemical calculations, was employed to compute the redox potential of the mechanism. This research provides valuable insights into the mechanisms and redox potentials involved in disulfide bridge reduction within proteins, offering an understanding of phenomena that are challenging to explore experimentally.

METHODS

All calculations used the Gaussian 09 software package at the MP2/6-311 + g(d,p) theory level. Visualization of the molecular orbitals and electron densities was conducted using Gaussview6. Molecular dynamics simulations were performed using GROMACS with the CHARMM36 force field. The PyMM program (Python Program for QM/MM Simulations Based on the Perturbed Matrix Method) is used to apply the Perturbed Matrix Method to MD simulations.

Strong Bases and beyond: The Prominent Contribution of Neutral Push-Pull Organic Molecules towards Superbases in the Gas Phase

In this review, the principles of gas-phase proton basicity measurements and theoretical calculations are recalled as a reminder of how the basicity PA/GB scale, based on Brønsted-Lowry theory, was constructed in the gas-phase (PA-proton affinity and/or GB-gas-phase basicity in the enthalpy and Gibbs energy scale, respectively). The origins of exceptionally strong gas-phase basicity of some organic nitrogen bases containing N-sp3 (amines), N-sp2 (imines, amidines, guanidines, polyguanides, phosphazenes), and N-sp (nitriles) are rationalized. In particular, the role of push-pull nitrogen bases in the development of the gas-phase basicity in the superbasicity region is emphasized. Some reasons for the difficulties in measurements for poly-functional nitrogen bases are highlighted. Various structural phenomena being in relation with gas-phase acid-base equilibria that should be considered in quantum-chemical calculations of PA/GB parameters are discussed. The preparation methods for strong organic push-pull bases containing a N-sp2 site of protonation are briefly reviewed. Finally, recent trends in research on neutral organic superbases, leaning toward catalytic and other remarkable applications, are underlined.

Intramolecular chaperone-assisted dual-anchoring activation (ICDA): a suitable preorganization for electrophilic halocyclization

The halocyclization reaction represents one of the most common methodologies for the synthesis of heterocyclic molecules. Many efforts have been made to balance the relationship between structure, reactivity and selectivity, including the design of new electrophilic halogenation reagents and the utilization of activating strategies. However, discovering universal reagents or activating strategies for electrophilic halocyclization remains challenging due to the case-by-case practice for different substrates or different cyclization models. Here we report an intramolecular chaperone-assisted dual-anchoring activation (ICDA) model for electrophilic halocyclization, taking advantage of the non-covalent dual-anchoring orientation as the driving force. This protocol allows a practical, catalyst-free and rapid approach to access seven types of small-sized, medium-sized, and large-sized heterocyclic units and to realize polyene-like domino halocyclizations, as exemplified by nearly 90 examples, including a risk-reducing flow protocol for gram-scale synthesis. DFT studies verify the crucial role of ICDA in affording a suitable preorganization for transition state stabilization and X+ transfer acceleration. The utilization of the ICDA model allows a spatiotemporal adjustment to straightforwardly obtain fast, selective and high-yielding synthetic transformations.

Estimation of thermodynamic and physicochemical properties of the alkali astatides: On the bond strength of molecular astatine (At2 ) and the hydration enthalpy of astatide (At- )

The recent accurate and precise determination of the electron affinity (EA) of the astatine atom At0 warrants a re-investigation of the estimated thermodynamic properties of At0 and astatine containing molecules as this EA was found to be much lower (by 0.4 eV) than previous estimated values. In this contribution we estimate, from available data sources, the following thermodynamic and physicochemical properties of the alkali astatides (MAt, M = Li, Na, K, Rb, Cs): their solid and gaseous heats of formation, lattice and gas-phase binding enthalpies, sublimation energies and melting temperatures. Gas-phase charge-transfer dissociation energies for the alkali astatides (the energy requirement for M+ At- ➔ M0  + At0 ) have been obtained and are compared with those for the other alkali halides. Use of Born-Haber cycles together with the new AE (At0 ) value allows the re-evaluation of ΔHf (At0 )g (=56 ± 5 kJ/mol); it is concluded that (At2 )g is a weakly bonded species (bond strength <50 kJ/mol), significantly weaker bonded than previously estimated (116 kJ/mol) and much weaker bonded than I2 (148 kJ/mol), but in agreement with the finding from theory that spin-orbit coupling considerably reduces the bond strength in At2 . The hydration enthalpy (ΔHaq ) of At- is estimated to be -230 ± 2 kJ/mol (using ΔHaq [H+ ] = -1150.1 kJ/mol), in good agreement with molecular dynamics calculations. Arguments are presented that the largest alkali halide, CsAt, like the smallest, LiF, will be only sparingly soluble in water, following the generalization from hard/soft acid/base principles that "small likes small" and "large likes large."

A fixed-charge model of the N-protomer of 4-aminobenzoic acid to facilitate the study of the unimolecular and bimolecular chemistry of its "neutral" carboxylic acid group

RATIONALE

There are a growing number of examples of protomers formed via electrospray ionization (ESI) that do not fragment under mobile proton conditions, giving rise to distinct tandem mass spectra. To model the N-protomer of 4-aminobenzoic acid, here we study the gas-phase unimolecular and bimolecular chemistry of the 4-(carboxyphenyl)trimethylammonium ion.

METHODS

4-(Carboxyphenyl)trimethylammonium iodide was synthesized, purified via recrystallization and transferred to the gas phase via ESI. 4-(Carboxyphenyl)trimethylammonium ion, 7, was mass selected and subjected to collision-induced dissociation and ion-molecule reactions in a linear ion trap mass spectrometer.

RESULTS

The major fragmentation channel for the fixed-charge cation 7 is methyl radical loss, whereas loss of trimethylamine and CO2 represents minor pathways. The free carboxylic acid functional group of 7 is unreactive toward a number of neutral reagents (methanol, acetone, acetonitrile, and N,N'-diisopropylcarbodiimide). 7 reacts very slowly with trimethylborate via addition-elimination, consistent with density functional theory (DFT) calculations that show this reaction is slightly endothermic. The deuterated cation 7(D) undergoes slow D/H exchange with ethanol, and DFT calculations reveal that a flip-flop mechanism operates.

CONCLUSIONS

The free carboxylic group of 7 is not very reactive toward neutral reagents in the gas phase.

Predicting a New Δ-Proton Sponge-Base of 4,12-Dihydrogen-4,8,12-triazatriangulene through Proton Affinity, Aromatic Stabilization Energy, and Aromatic Magnetism

Herein, we report designing a new Δ (delta-shaped) proton sponge base of 4,12-dihydrogen-4,8,12-triazatriangulene (compound 1) and calculating its proton affinity (PA), aromatic stabilization, natural bond orbital (NBO), electron density ρ(r), Laplacian of electron density ∇2 ρ(r), (2D-3D) multidimensional off-nucleus magnetic shielding (σzz (r) and σiso (r)), and scanning nucleus-independent chemical shift (NICSzz and NICS). Density functional theory (DFT) at B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP were used to compute the magnetic shielding variables. In addition, relevant bases like pyridine, quinoline, and acridine were also studied and compared. The protonation of compound 1 yields a highly symmetric carbocation of three Hückel benzenic rings. Comparing our findings of the studied molecules showed that compound 1 precedes others in PA, aromatic isomerization stabilization energy, and basicity. Therefore, the basicity may be enhanced when a conjugate acid gains higher aromatic features than its unprotonated base. Both multidimensional σzz (r) and σiso (r) off-nucleus magnetic shieldings outperformed electron-based techniques and can visually monitor changes in aromaticity that occur by protonation. The B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP levels showed no significant differences in detailing isochemical shielding surfaces.

Protonation of serine: conformers, proton affinities and gas-phase basicities at the "gold standard" and beyond

The potential energy surfaces (PESs) of serine and its protonated counterparts are investigated to determine the structures of the minima. A total of 95 neutral serine, 15 N-(amino-) and 46 O-(carbonyl-)protonated serine conformers are found. Their relative energies, geometries and harmonic vibrational frequencies are determined at the MP2/aug-cc-pVDZ level of theory. To obtain highly accurate thermodynamic values, further computations are performed: the ten conformers with the lowest relative energies from each molecule type (neutral, N- and O-protonated) are further optimized using the explicitly correlated CCSD(T)-F12a/cc-pVDZ-F12 method (for neutral serine, harmonic vibrational frequencies were also computed). In addition, auxiliary corrections were determined: basis-set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation and second-order Douglas-Kroll relativistic effects along with zero-point energy contributions. Two important thermodynamic parameters (at 298.15 K), proton affinity (PA)/gas-phase basicity (GB) are calculated considering the two different protonation sites: 218.05 ± 0.2/209.86 ± 0.6 kcal mol^-1 and 205.87 ± 0.2/196.36 ± 0.3 kcal mol^-1 for the amino and carbonyl sites, respectively. The uncertainty of the determined values was approximated based on various sources including auxiliary corrections, basis-set effects, harmonic vibrational frequencies.

Mechanistic DFT Study of 1,3-Dipolar Cycloadditions of Azides with Guanidine

Density functional calculations SMD(chloroform)//B3LYP/6-311+G(2d,p) were employed in the computational study of 1,3-dipolar cycloadditions of azides with guanidine. The formation of two regioisomeric tetrazoles and their rearrangement to cyclic aziridines and open-chain guanidine products were modeled. The results suggest the feasibility of an uncatalyzed reaction under very drastic conditions since the thermodynamically preferred reaction path (a), which involves cycloaddition by binding the carbon atom from guanidine to the terminal azide nitrogen atom, and the guanidine imino nitrogen with the inner N atom from the azide, has an energy barrier higher than 50 kcal mol^-1. The formation of the other regioisomeric tetrazole (imino nitrogen interacts with terminal N atom of azide) in direction (b) can be more favorable and proceed under milder conditions if alternative activation of the nitrogen molecule releases (e.g., photochemical activation), or deamination could be achieved because these processes have the highest barrier in the less favorable (b) branch of the mechanism. The introduction of substituents should favorably affect the cycloaddition reactivity of the azides, with the greatest effects expected for the benzyl and perfluorophenyl groups.

The importance of four-membered NHCs in stabilizing Breslow intermediates on benzoin condensation pathway

N-heterocyclic carbenes (NHCs) have been established to be effective organocatalysts for facilitating the benzoin condensation and many other reactions. These reactions involve the formation of a Breslow intermediate (BI), which exhibits umpolung chemistry. To facilitate organocatalysis, several new cyclic carbenes are being introduced, four-membered NHCs are of special interest. Whether these NHCs can exhibit catalytic influence or not, can be evaluated by exploring the potential energy surface (PES) of the benzoin condensation reaction. Quantum chemical analysis has been carried out to compare the PES of these four-membered NHCs with that of standard five-membered NHCs to explore their catalytic ability. The barrier for the first step of the reaction for the formation of BI is comparable in all the cases. But the barrier for the second step of the reaction leading to the benzoin formation from BI is estimated to be very high for the four membered NHCs. These results indicate that the probability of identifying and isolating the BI is very high in comparison to the completion of benzoin condensation reaction in the case of the four-membered NHCs.

Benchmark Ab Initio Determination of the Conformers, Proton Affinities, and Gas-Phase Basicities of Cysteine

A systematic conformational mapping combined with literature data leads to 85 stable neutral cysteine conformers. The implementation of the same mapping process for the protonated counterparts reveals 21 N-(amino-), 64 O-(carbonyl-), and 37 S-(thiol-)protonated cysteine conformers. Their relative energies and harmonic vibrational frequencies are given at the MP2/aug-cc-pVDZ level of theory. Further benchmark ab initio computations are performed for the 10 lowest-lying neutral and protonated amino acid conformers (for each type) such as CCSD(T)-F12a/cc-pVDZ-F12 geometry optimizations (and frequency computations for cysteine) as well as auxiliary correction computations of the basis set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation effects, second-order Douglass-Kroll relativistic effects, and zero-point energy contributions. Boltzmann-averaged 0 (298.15) K proton affinity and [298.15 K gas-phase basicity] values of cysteine are predicted to be 214.96 (216.39) [208.21], 201.83 (203.55) [194.16], and 193.31 (194.74) [186.40] kcal/mol for N-, O-, and S-protonation, respectively, also considering the previously described auxiliary corrections.

Pyridoneimine-catalyzed anomeric aqueous oxa-Michael additions of native mono- and disaccharides

A pyridoneimine-catalyzed oxa-Michael addition of protecting groups-free, native mono- and disaccharides with Michael acceptors in aq. solution is reported. Several mono- and disaccharides are reacted with acceptors, namely, methylvinyl ketone, acrylonitrile and tert-butyl acrylate in aq. solution, the addition catalyzed by n-pentylpyridone imine. The addition occurs site-selectively at the anomeric lactol and the remaining hydroxy functionalities are un-affected. The resulting keto-glycopyranoside products are explored in aldol, allylation and oxime product formation, occurring at either α-methyl moiety or at the keto-moiety, with appropriate synthons. In another direction, the keto-glycopyranoside is functionalized further with amino acids through reductive amination in aq. methanol solution. Formation of hemiacetal anion occurs in the presence of pyridoneimine in aq. Solution, enabling subsequent addition to occur with acceptors. Facile reductive amination of the resulting keto-glycoside provides an avenue for conjugations with amino acids in the present work.

Nitriles with High Gas-Phase Basicity-Part II Transmission of the Push-Pull Effect through Methylenecyclopropene and Cyclopropenimine Scaffolds Intercalated between Different Electron Donor(s) and the Cyano N-Protonation Site

This work extends our earlier quantum chemical studies on the gas-phase basicity of very strong N-bases to two series of nitriles containing the methylenecyclopropene and cyclopropenimine scaffolds with dissymmetrical substitution by one or two electron-donating substituents such as Me, NR2, N=C (NR2)2, and N=P (NR2)3, the last three being strong donors. For a proper prediction of their gas-phase base properties, all potential isomeric phenomena and reasonable potential protonation sites are considered to avoid possible inconsistencies when evaluating the energetic parameters and associated protonation or deprotonation equilibria B + H+ = BH+. More than 250 new isomeric structures for neutral and protonated forms are analyzed. The stable structures are selected and the favored ones identified. The microscopic (kinetic) gas-phase basicity parameters (PA and GB) corresponding to N sites (cyano and imino in the cyclopropenimine or in the substituents) in each isomer are calculated. The macroscopic (thermodynamic) PAs and GBs, referring to the isomeric mixtures of favored isomers, are also estimated. The total (pushing) substituent effects are analyzed for monosubstituted and disubstituted derivatives containing two identical or two different substituents. Electron delocalization is examined in the two π-π conjugated transmitters, the methylenecyclopropene and cyclopropenimine scaffolds. The aromatic character of the three-membered ring is also discussed.

Conformational Changes Induced by Methyl Side-Chains in Protonated Tripeptides Containing Glycine and Alanine Residues

We present a systematic study of the conformational and isomeric populations in gas-phase protonated tripeptides containing glycine and alanine residues using infrared predissociation spectroscopy of cryogenically cooled ions. Specifically, the protonated forms of Gly-Gly-Gly, Ala-Gly-Gly, Gly-Ala-Gly, Gly-Gly-Ala, Ala-Ala-Gly, Ala-Gly-Ala, Gly-Ala-Ala, and Ala-Ala-Ala allow us to sample all permutations of the methyl side-chain position, providing a comprehensive view of the effects of this simple side-chain on the 3-D structure of the peptide. The individual structural populations for all but one of these peptide species are determined via conformer-specific IR-IR double-resonance spectroscopy and comparison with electronic structure predictions. The observed structures can be classified into three main families defined by the protonation site and the number of internal hydrogen bonds. The relative contribution of each structural family is highly dependent on the exact amino acid sequence of the tripeptide. These observed changes in structural population can be rationalized in terms of the electron-donating effect of the methyl side-chain modulating the local proton affinities of the amine and various carbonyl groups in the tripeptide.

Strong Bases Design: Predicted Limits of Basicity

Brønsted superbases have wide applications in organic chemistry due to their ability to activate C-H bonds. The strongest neutral bases to date are substituted aminophosphazenes developed in the late 1980s by Reinhard Schwesinger. Since then, much effort has been expended to create even stronger neutral bases. In this article, the reasons for the instability of very basic compounds are investigated by means of high-level quantum-chemical calculations. Theoretical basicity limits are suggested for solutions as well as for the gas phase. A record-breaking superbase most likely to be synthesizable and stable at ambient conditions is proposed. Hexamethylphosphoramide is considered a reliable ionizing solvent for superbases.

Proton Affinity Values of Fentanyl and Fentanyl Analogues Pertinent to Ambient Ionization and Detection

Proton affinity is a major factor in the atmospheric pressure chemical ionization of illicit drugs. The detection of illicit drugs by mass spectrometry and ion mobility spectrometry relies on the analytes having greater proton affinities than background species. Evaluating proton affinities for fentanyl and its analogues is informative for predicting the likelihood of ionization in different environments and for optimizing the compounds' ionization and detection, such as through the addition of dopant chemicals. Herein, density functional theory was used to computationally determine the proton affinity and gas-phase basicity of 15 fentanyl compounds and several relevant molecules as a reference point. The range of proton affinities for the fentanyl compounds was from 1018 to 1078 kJ/mol. Fentanyl compounds with the higher proton affinity values appeared to form a bridge between the oxygen on the amide and the protonated nitrogen on the piperidine ring based on models and calculated bond distances. Experiments with fragmentation of proton-bound clusters using atmospheric flow tube-mass spectrometry (AFT-MS) provided estimates of relative proton affinities and showed proton affinity values of fentanyl compounds >1000 kJ/mol, which were consistent with the computational results. The high proton affinities of fentanyl compounds facilitate their detection by ambient ionization techniques in complex environments. The detection limits of the fentanyl compounds with AFT-MS are in the low femtogram range, which demonstrates the feasibility of trace vapor drug detection.

Exploring and modeling the ion mobility spectrometry of perindopril: Example of protonation-dissociation reactions in large molecules

The current research is constructed for considering the chemical ionization and dissociation of perindopril in the positive mode of corona discharge ion mobility spectrometry. Four product ion peaks are observed in the ion mobility spectrum of perindopril erbumine at the cell temperature of 473 K. These peaks are assigned through the obtained intensity variation analysis in the ion mobility spectra over the elapsed time accompanied by the calculations backed by the validated density functional theory (DFT). In this regard, the most stable ionic species associated with each peak and the corresponding reliable generation pathways are found by the well-confirmed meta hybrid density functional method, M06-2X. The peaks are assigned to the protonated perindopril and its dissociation products, including counter ion and the related fragment ions. However, the structures of the neutral perindopril in the gas phase are thoroughly assessed to find a more stable one. The predicted chemical ionization products by the theory are in excellent agreement with our presented experiment here. Theoretical evaluations demonstrated that the production of a fragment by dissociation process occurs when perindopril gets a proton from the ionization region. Also, without protons, there is no dissociation process. Therefore, our mechanism investigated here is the proton transfer one. All possible sites of perindopril are considered theoretically for protonation along with their possible reactions. In addition to the computed PES, the assigned ions for obtained spectra are confirmed by the computed equilibrium constants and rate constants. Our theoretical results show that the peak of the main fragment is for M-CH₃ CH₂ OH produced by a reaction pathway involving no barrier. This study opens new perspectives in interpreting large molecules spectra for future studies.

Proton affinity and gas phase basicity of diamandoid molecules: diamantane to C131H116

Calculated proton affinities (PAs) and gas phase basicities (GPBs) are reported for diamantane (C₁₄H₂₀), triamantane (C₁₈H₂₄), 'globular and planar' isomers of tetramantane (C₂₂H₂₈) and pentamantane (C₂₆H₃₂), and for one 'globular' isomer of each of the larger diamondoid molecules: C₅₁H₅₈, C₇₈H₇₂, C₁₀₂H₉₀, and C₁₃₁H₁₁₆. Assuming C_xH_y as the parent diamondoid molecule, we calculated PA and GPB values for a variety of C_xH_y+1⁺ isomers, as well as for the reaction C_xH_y + H⁺ yielding C_xH_y-1⁺ + H₂(g); the latter is slightly favored based on GPB values for diamantane through pentamantane, but less favored compared to certain C_xH_y+1⁺ isomers of C₅₁H₅₈, C₁₀₂H₉₀, and C₁₃₁H₁₁₆. Indeed, the GPB values of C₅₁H₅₈, C₁₀₂H₉₀, and C₁₃₁H₁₁₆ classifiy them as 'superbases'. Calculations that had the initial location of the proton in an interstitial site inside the diamondoid molecule always showed the H having moved to the outside of the diamondoid molecule; for this reason, we focused on testing a variety of initial configurations with the proton placed in an initial position on the surface. Additional protons were added to determine the limiting number that could be, per these calculations, taken up by the diamondoid molecules and the maximum number of protons are shown in parentheses: C₁₄H₂₀(2), C₁₈H₂₄(3), C₂₂H₂₈(3), C₂₆H₃₂(3), C₅₁H₅₈(4). Bader charge distributions obtained for C_xH_y+1⁺ isomers (for diamantane through pentamantane) suggest that the positive charge is essentially completely delocalized over all the H atoms. NMR spectra were calculated for different isomers of C₁₄H₁₉⁺, and compared to the published NMR spectrum for when diamantane was mixed with magic acid and H₂(g) was produced.

Acid and base strength variations: rationalization for cyclic amine bases and acidic aqua cations

This perspective highlights and evaluates recent key developments in the thermodynamic approach used to analyze trends in acid and base strength variation. According to this approach, acid and base strength ranking can be interpreted by using thermodynamic or thermochemical cycles. Each cycle generally consists of three independent but well-defined steps. The modus operandi described here entails the identification of the dominant step and the rationalization of its free energy/enthalpy/energy change along a selected series in terms of known structural chemical concepts. Developments in this approach are described by focusing on two related series of bases and two series of acids. In the case of the former the protonation of a series of N-heterocyclic amine bases together with their methyl-substituted analogs receives particular attention while in the case of acids, the acidic properties of aqua dications of elements in period 4 and group 2 are probed. It is illustrated how significant progress in computational chemistry and mass spectrometric techniques can be employed to compare 'inherent' basicity or acidity in the selected families of compounds by using simple gas-phase energy cycles. Unique, dual functions for both electronegativity (element and orbital) and charge density (for aqua cations) indicators are identified and used to evaluate these cycles. Solvent effects (in aqueous solution) are accommodated by including dehydration and hydration changes in appropriately-extended, three-step free energy cycles. It is further suggested that the dominant step in the extended thermodynamic cycle for monomeric aqua cations is the transfer of M(H₂O)_n²⁺ complex hydrates from the gas-phase to bulk water. Charge density of the aqua cations again features prominently in proposed rationalizations. Finally, this article also sheds light on salient relationships that exist between empirically and quantum-chemically estimated enthalpy and entropy changes for the aforementioned transfer process.

Stable Singlet Carbenes as Organic Superbases

A simple experimental procedure for scaling carbene Brønsted basicity is described. The results highlight the strong basicity of pyrazol-4-ylidenes, a type of mesoionic carbene, also named cyclic-bentallenes (CBA). They are more basic (pK_aH >42.7 in acetonitrile) than the popular proazaphosphatrane Verkade bases, and even the Schwesinger phosphazene superbase P₄ (^t Bu). The basicity of these compounds can readily be tuned, and they are accessible in multigram quantities. These results open new avenues for carbon centered superbases.

Push–Pull Effect on the Gas-Phase Basicity of Nitriles: Transmission of the Resonance Effects by Methylenecyclopropene and Cyclopropenimine π-Systems Substituted by Two Identical Strong Electron Donors

The gas-phase basicity of nitriles can be enhanced by a push–pull effect. The role of the intercalated scaffold between the pushing group (electron-donor) and the pulling (electron-acceptor) nitrile group is crucial in the basicity enhancement, simultaneously having a transmission function and an intrinsic contribution to the basicity. In this study, we examine the methylenecyclopropene and the N-analog, cyclopropenimine, as the smallest cyclic π systems that can be considered for resonance propagation in a push–pull system, as well as their derivatives possessing two strong pushing groups (X) attached symmetrically to the cyclopropene scaffold. For basicity and push–pull effect investigations, we apply theoretical methods (DFT and G2). The effects of geometrical and rotational isomerism on the basicity are explored. We establish that the protonation of the cyano group is always favored. The push–pull effect of strong electron donor X substituents is very similar and the two π-systems appear to be good relays for this effect. The effects of groups in the two cyclopropene series are found to be proportional to the effects in the directly substituted nitrile series X–C≡N. In parallel to the basicity, changes in electron delocalization caused by protonation are also assessed on the basis of aromaticity indices. The calculated proton affinities of the nitrile series reported in this study enrich the gas-phase basicity scale of nitriles to around 1000 kJ mol−1.

Design of Novel Uncharged Organic Superbases: Merging Basicity and Functionality

ConspectusOne of the constant challenges of synthetic chemistry is the molecular design and synthesis of nonionic, metal-free superbases as chemically stable neutral organic compounds of moderate molecular weight, intrinsically high thermodynamic basicity, adaptable kinetic basicity, and weak or tunable nucleophilicity at their nitrogen, phosphorus, or carbon basicity centers. Such superbases can catalyze numerous reactions, ranging from C-C bond formation to cycloadditions and polymerization, to name just a few. Additional benefits of organic superbases, as opposed to their inorganic counterparts, are their solubility in organic reaction media, mild reaction conditions, and higher selectivity. Approaching such superbasic compounds remains a continuous challenge. However, recent advances in synthetic methodology and theoretical understanding have resulted in new design principles and synthetic strategies toward superbases. Our computational contributions have demonstrated that the gas-phase basicity region of 350 kcal mol^-1 and even beyond is easily reachable by organosuperbases. However, despite record-high basicities, the physical limitations of many of these compounds become quickly evident. The typically large molecular weight of these molecules and their sensitivity to ordinary reaction conditions prevent them from being practical, even though their preparation is often not too difficult. Thus, obviously structural limitations with respect to molecular weight and structural complexity must be imposed on the design of new synthetically useful organic superbases, but strategies for increasing their basicity remain important.The contemporary design of novel organic superbases is illustrated by phosphazenyl phosphanes displaying gas-phase basicities (GB) above 300 kcal mol^-1 but having molecular weights well below 1000 g·mol^-1. This approach is based on a reconsideration of phosphorus(III) compounds, which goes along with increasing their stability in solution. Another example is the preparation of carbodiphosphoranes incorporating pyrrolidine, tetramethylguanidine, or hexamethylphosphazene as a substituent. With gas-phase proton affinities of up to 300 kcal mol^-1, they are among the top nonionic carbon bases on the basicity scale. Remarkably, the high basicity of these compounds is achieved at molecular weights of around 600 g·mol^-1. Another approach to achieving high basicity through the cooperative effect of multiple intramolecular hydrogen bonding, which increases the stabilization of conjugate acids, has recently been confirmed.This Account focuses on our efforts to produce superbasic molecules that embody many desirable traits, but other groups' approaches will also be discussed. We reveal the crucial structural features of superbases and place them on known basicity scales. We discuss the emerging potential and current limits of their application and give a general outlook into the future.

Phosphorus-Containing Superbases: Recent Progress in the Chemistry of Electron-Abundant Phosphines and Phosphazenes

The renaissance of Brønsted superbases is primarily based on their pronounced capacity for a large variety of chemical transformations under mild reaction conditions. Four major set screws are available for the selective tuning of the basicity: the nature of the basic center (N, P, …), the degree of electron donation by substituents to the central atom, the possibility of charge delocalization, and the energy gain by hydrogen bonding. Within the past decades, a plethora of neutral electron-rich phosphine and phosphazene bases have appeared in the literature. Their outstanding properties and advantages over inorganic or charged bases have now made them indispensable as auxiliary bases in deprotonation processes. Herein, an update of the chemistry of basic phosphines and phosphazenes is given. In addition, due to widespread interest, their use in catalysis or as ligands in coordination chemistry is highlighted.

Aggregation and Solvation of Sodium Hexamethyldisilazide: Across the Solvent Spectrum

We report solution structures of sodium hexamethyldisilazide (NaHMDS) solvated by >30 standard solvents (ligands). These include: toluene, benzene, and styrene; triethylamine and related trialkylamines; pyrrolidine as a representative dialkylamine; dialkylethers including THF, tert-butylmethyl ether, and diethyl ether; dipolar ligands such as DMF, HMPA, DMSO, and DMPU; a bifunctional dipolar ligand nonamethylimidodiphosphoramide (NIPA); polyamines N,N,N',N'-tetramethylenediamine (TMEDA), N,N,N',N″,N″-pentamethyldiethylenetriamine (PMDTA), N,N,N',N'-tetramethylcyclohexanediamine (TMCDA), and 2,2'-bipyridine; polyethers 12-crown-4, 15-crown-5, 18-crown-6, and diglyme; 4,7,13,16,21,24-hexaoxa-1,10-diazabicyclo[8.8.8]hexacosane ([2.2.2] cryptand); and tris[2-(2-methoxyethoxy)ethyl]amine (TDA-1). Combinations of 1H, 13C, 15N, and 29Si NMR spectroscopies, the method of continuous variations, X-ray crystallography, and density functional theory (DFT) computations reveal ligand-modulated aggregation to give mixtures of dimers, monomers, triple ions, and ion pairs. 15N-29Si coupling constants distinguish dimers and monomers. Solvation numbers are determined by a combination of solvent titrations, observed free and bound solvent in the slow exchange limit, and DFT computations. The relative abilities of solvents to compete in binary mixtures often match that predicted by conventional wisdom but with some exceptions and evidence of both competitive and cooperative (mixed) solvation. Crystal structures of a NaHMDS cryptate ion pair and a 15-crown-5-solvated monomer are included. Results are compared with those for lithium hexamethyldisilazide, lithium diisopropylamide, and sodium diisopropylamide.

Dissociation of Singly and Multiply Charged Nitromethane Cations: Femtosecond Laser Mass Spectrometry and Theoretical Modeling

Dissociation pathways of singly- and multiply charged gas-phase nitromethane cations were investigated with strong-field laser photoionization mass spectrometry and density functional theory computations. There are multiple isomers of the singly charged nitromethane radical cation, several of which can be accessed by rearrangement of the parent CH₃-NO₂ structure with low energy barriers. While direct cleavage of the C-N bond from the parent nitromethane cation produces NO₂⁺ and CH₃⁺, rearrangement prior to dissociation accounts for fragmentation products including NO⁺, CH₂OH⁺, and CH₂NO⁺. Extensive Coulomb explosion in fragment ions observed at high laser intensity indicates that rapid dissociation of multiply charged nitromethane cations produces additional species such as CH₂⁺, H⁺, and NO₂²⁺.  On the basis of analysis of Coulomb explosion in the mass spectral signals and pathway calculations, sufficiently intense laser fields can remove four or more electrons from nitromethane.

Synthesis and Basicity Studies of Quinolino[7,8-h]quinoline Derivatives

Quinolino[7,8-h]quinoline is a superbasic compound, with a pK_aH in acetonitrile greater than that of 1,8-bis(dimethylaminonaphthalene) (DMAN), although its synthesis and the synthesis of its derivatives can be problematic. The use of halogen derivatives 4,9-dichloroquinolino[7,8-h]quinoline (16) and 4,9-dibromoquinolino[7,8-h]quinoline (17) as precursors has granted the formation of a range of substituted quinolinoquinolines. The basicity and other properties of quinolinoquinolines can be modified by the inclusion of suitable functionalities. The experimentally obtained pK_aH values of quinolino[7,8-h]quinoline derivatives show that N⁴,N⁴,N⁹,N⁹-tetraethylquinolino[7,8-h]quinoline-4,9-diamine (26) is more superbasic than quinolino[7,8-h]quinoline. Computationally derived pK_aH values of quinolinoquinolines functionalized with dimethylamino (NMe₂), 1,1,3,3-tetramethylguanidino (N═C(NMe₂)₂) or N,N,N',N',N″,N″-hexamethylphosphorimidic triamido (N═P(NMe₂)₃) groups are significantly greater than those of quinolino[7,8-h]quinoline. Overall, electron-donating functionalities are observed to increase the basicity of the quinolinoquinoline moiety, while the substitution of electron-withdrawing groups lowers the basicity.

Twisted Eigen Can Induce Proton Transfer at a Hydrophobic-Hydrophilic Interface

The investigation of proton localization at a hydrophobic-hydrophilic interface is an important problem in chemical and materials sciences. In this study, protonated benzene (i.e., benzenium ion) and water clusters [BZH⁺W_n (where n = 1-6)] are selected as prototype models to understand the interfacial interactions and proton transfer mechanism between a carbonaceous surface and water molecules. The excess protons can localize in the vicinity of the hydrophobic-hydrophilic interface, and these clusters are stabilized by various kinds of noncovalent interactions. Calculations are carried out using ab initio (MP2) and density functional theory B3LYP methods to shed more light on geometries, energetics, and spectral signatures of the protonated species [H⁺(H₂O)_n] at the interfaces. These calculations revealed few low-lying isomers, which have not been reported earlier. Scrutiny of the results reveals that proton localization in the hydrophilic environment is more stable than the hydrophobic benzene π-cloud. Furthermore, the occurrence of an O-H⁺···π hydrogen bond significantly influences the O-H⁺···O interactions in the water clusters and also intensively affects the vibrational modes of the Eigen cation. Thus, the aromatic π-clouds can stabilize the Eigen cation and at the same time, a twisted form of Eigen (one O-H⁺···π → two O-H⁺···π) can enhance the proton transfer through the water chain via a Grotthuss-type mechanism. The vibrational spectra of these clusters reveal that there is a large red-shifted frequency for the O-H⁺···O, O-H⁺···π, and O-H···π modes of interaction. The energetic values and vibrational frequencies obtained from the B3LYP method are in close agreement with the MP2 level and experimental values, respectively.

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Stable Mesoionic N-Heterocyclic Olefins (mNHOs)

We report a new class of stable mesoionic N-heterocyclic olefins, featuring a highly polarized (strongly ylidic) double bond. The ground-state structure cannot be described through an uncharged mesomeric Lewis-structure, thereby structurally distinguishing them from traditional N-heterocyclic olefins (NHOs). mNHOs can easily be obtained through deprotonation of the corresponding methylated N,N'-diaryl-1,2,3-triazolium and N,N'-diaryl-imidazolium salts, respectively. In their reactivity, they represent strong σ-donor ligands as shown by their coordination complexes of rhodium and boron. Their calculated proton affinities, their experimentally derived basicities (competition experiments), as well as donor abilities (Tolman electronic parameter; TEP) exceed the so far reported class of NHOs.