Experimental Basicities of Phosphazene ... - ACS Publications

Apr 19, 2016 - Institute of Chemistry, University of Tartu, Ravila 14a Str, 50411 ... $Institute of Computer Sciences, University of Tartu, J. Liivi 2...
1 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCA

Experimental Basicities of Phosphazene, Guanidinophosphazene, and Proton Sponge Superbases in the Gas Phase and Solution Ivari Kaljurand,*,† Jaan Saame,† Toomas Rodima,† Ivar Koppel,$ Ilmar A. Koppel,† Julius F. Kögel,‡,§ Jörg Sundermeyer,‡ Uwe Köhn,# Martyn P. Coles,∥ and Ivo Leito† †

Institute of Chemistry, University of Tartu, Ravila 14a Str, 50411 Tartu, Estonia Institute of Computer Sciences, University of Tartu, J. Liivi 2 Str, 50409 Tartu, Estonia ‡ Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße, 35032 Marburg, Germany # Institut für Organische Chemie und Makromolekulare Chemie, Friedrich-Schiller-Universität Jena, Humboldtstraße 10, 07743 Jena, Germany ∥ School of Chemical and Physical Sciences, Victoria University of Wellington, P.O. Box 600, Wellington, New Zealand $

S Supporting Information *

ABSTRACT: Experimental gas-phase superbasicity scale spanning 20 orders of magnitude and ranging from bicyclic guanidine 7-methyl-1,5,7triazabicyclo[4.4.0]dec-5-ene to triguanidinophosphazenes and P3 phosphazenes is presented together with solution basicity data in acetonitrile and tetrahydrofuran. The most basic compound in the scale triguanidinophosphazene Et−NP[NC(NMe2)2]3has the highest experimental gas-phase basicity of an organic base ever reported: 273.9 kcal mol−1. The scale includes besides the higher homologues of classical superbasic phosphazenes and several guanidino-substituted phosphazenes also a number of recently introduced bisphosphazene and bis-guanidino proton sponges. This advancement was made possible by a newly designed Fourier transform ion cyclotron resonance (ICR) mass spectrometry setup with the unique ability to generate and control in the ICR cell sufficient vapor pressures of two delicate compounds having low volatility, which enables determining their basicity difference. The obtained experimental gas-phase and solution basicity data are analyzed in terms of structural and solvent effects and compared with data from theoretical calculations.



INTRODUCTION

in electrospray ionization (ESI) liquid chromatography (LC) mass spectrometry (MS) methods.36−40 Experimental GB and proton affinity (PA) data of neutral superbasic compounds up to the year 2005 were collected into NIST Webbook.4 Since then some progress has been made in extending the experimental GB scale in the superbasic region,41−49 but in most of the published works bases having GB values below 260 kcal mol−1 have been used. This is mainly due to the experimental difficulties in preparation and handling of superbasic compounds and limited availability of reliable and suitable reference compounds. At the same time the progress in design and solution studies of superbases has been impressive. Many new classes of nonionic superbases have been designed, and many of them have been prepared during the past decade. Among these are guanidinophosphazenes,50−52 various types of proton sponges, 5 3 − 6 9 superbasic guanidine derivatives,11,42,70−76 macrocyclic proton chelators,77,78 chiral phosphazenes,6,79−81 quinones,82 phosphonium ylides,83,84,52 aza-

1

Nonionic superbasic compounds have achieved prominent status in several fields of chemistry and materials science. Commercially available and novel guanidine, phosphazene, carbene, and other types of nonionic superbases have been used for many applications. Best known are applications in synthesis as catalysts.2,6−14 Recently several commercially available amidine, guanidine, and P1 and P2 phosphazene superbases have been used for preparation of thermally stable protic ionic liquids (PILs).15−19 The stability and the physical and electrochemical properties of these PILs were reported to be similar to common aprotic ionic liquids.15,18,20 These and other smaller and cheaper superbases as constituents of PILs may find use in several fields of energy technologies,21,22 reversible absorption of CO2,23,24 cellulose processing,25,26 in synthesis, catalysis, and many other applications.17,27 Superbasic carbenes have been reported to form in aprotic ionic liquids.27−29 Gas-phase (or intrinsic) basicity (GB) of a basic molecule is its important property, largely determining its basicity also in solution and being directly important, for example, in various types of laser desorption mass spectrometry methods,30−35 and © XXXX American Chemical Society

Received: February 15, 2016 Revised: April 1, 2016

A

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A phosphatranes,83,84 carbenes,45,52,85−89 etc.48,90−96 Theoretical calculations predict that many of these bases are more basic than the most basic experimentally measured compoundP2phosphazene EtP2(dma)5 (dma = dimethylamino; see Chart 1 for its structure; GB = 264 kcal mol−1).41 This work was undertaken for four reasons. First, to obtain reliable experimental GB values for a representative set of superbases of diverse chemical nature and to compare the results with GB values from theoretical calculations. Second, to explore various families of superbasic compounds and find promising (in terms of stability, inertness, etc.) superbasic

compound families for further extending the existing GB ladder toward stronger bases and thus forming foundation for future extension of experimental GB scale. Third, to explore structure−basicity relationships of strong bases in the gas phase and compare these with basicity trends in solvents. Fourthly, to develop and test the instrumentation and experimental methodology for measuring GB values of large superbasic molecules of very low volatility. To fulfill this task the samples were chosen from the available set of superbasic compounds presented in Charts 1 and 2 and having calculated or estimated GB values above 255 kcal mol−1 for extension of experimental GB scale by using Fourier transform (FT) ion cyclotron resonance (ICR) mass spectrometry (MS) equilibrium method.

Chart 1. Structures and Abbreviations of P1, P2, P3, and P4 Phosphazenes

Chart 2. Structures and Abbreviations of Guanidinophosphazenes,a Proton Sponges, and Other Bases

a

B

tmg refers to −NC(NMe2)2. DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

region during experiments was evaluated to be 100 °C ± 10 °C (see Supporting Information for the methodology and exceptions). Volatile superbases were introduced as vapors from an in-house built short heated (up to 150 °C) stainless steel inlet line, and nonvolatile superbases were introduced as vapors from a custom-made dual solid sample probe accessory. Both inlet types had the possibility for degassing samples with turbopumping (at pressure down to 1 × 10−8 mbar) before introduction of their vapors into the high vacuum region of the FT-ICR cell. Volatile compound samples were carefully degassed with repeated freeze−pump−thaw cycles by using liquid nitrogen for sample freezing or degassed at room temperature. Dual solid probe accessory consisted of gate valved and turbopumped vacuum lock, fitted for simultaneous use of two solid sample probes (Scientific Instrument Services, NJ, USA). On the tip of each probe was a heatable holder for attaching a small Pyrex glass tube with a small amount of pure compound. Air-sensitive bases were loaded into the tube in inert gas (Ar, 5.0) glovebox,5 taken to the FT-ICR MS laboratory in vial under the same inert gas, mounted to the probe tip, and taken into the vacuum lock in in-house built small mobile glovebox under inert gas (N2, 5.0) flow. Probes were inserted through the vacuum lock so that their tips were located in the vicinity of the FT-ICR cell. The distance of each probe tip from the FT-ICR cell was varied from 1 to 5 cm depending on the needed probe tip temperature. When using higher probe tip temperature then it was kept at a longer distance from the cell in order not to influence significantly the cell temperature. The probe tip with sample in the glass tube was heated electrically with temperature controller or cooled with nitrogen gas flow through built-in solid probe cooling line to achieve a suitable sample vaporization rate and vapor pressure of the base in the FT-ICR cell region. Needed probe tip temperatures for creating appropriate vapor pressure of compound and pressure ratios of compounds are presented in Table S1 in Supporting Information. Because of the low volatility of the used bases, high risk to contaminate the ultrahigh vacuum system with excess of the sample compounds and need of long time for equilibration of the partial pressure of the compound in the vacuum system it took typically several days to introduce new superbase to appropriate vapor pressure level in the FT-ICR cell. The time for stabilization of the partial pressure ratio in the mixture of two low-volatile superbasic compounds in gas phase inside the FT-ICR cell was also long, ranging from several hours to one week. Stabilization of the constitution of gas phase was monitored by following the change of pressure ratio of the bases in constant short reaction time mass spectra over time (see next paragraph). Gas-phase composition was considered stable when the level of volatile impurities (mainly traces of solvents, water, and starting materials and side products of superbase synthesis, monitored by taking short reaction time MS) of sample had become small in the gas phase, and the compound signal (BH+, plus in some cases also B+) had become intense and stable. The same sample introduction procedure was repeated for the second compound. Typical pressure reading (uncorrected) range during experiments was around 2 × 10−8 mbar. Hot-filament electron ionization (rhenium ribbon filament) was used for ionization of small fraction of neutral gas mixture of bases B1 and B2 directly inside the FT-ICR cell and thus introduces cationic fragments, which are able to donate their acidic proton to the neutral superbasic molecules. Following this mobile proton generation step a series of experiments were performed using varied

Gas-phase basicity and proton affinity refer to the following equilibrium: ΔG b , ΔHb

B + H+ XoooooooooY BH+

(1)

GB and PA are defined as follows: GB = −ΔG b

PA = −ΔHb

(2)

The directly experimentally measured quantity is the relative basicity of two bases ΔΔGb: ΔΔG b

B2 + B1H+ XoooooY B2H+ + B1

(3)

where ΔΔG b = ΔG b(B2) − ΔG b(B1) = −RT ln K

K=

p(B1) ·I(B2H+) p(B2) ·I(B1H+)

(4)

(5)

p(B1) and p(B2) are the partial pressures of the respective neutral species. I(B1H+) and I(B2H+) are the ratios of ion intensities in the mass spectra (with necessary corrections, see Experimental Methods) and are used for obtaining the ratio of abundances of the protonated species B1H+ and B2H+. Several experimental techniques have been used for determination of GB and PA values of bases,4 of which FTICR MS97 and high-pressure mass spectrometry (HPMS)98 are historically the best known and considered most reliable. Kinetics methods99,42 on various MS instruments have also been used extensively over the last decades and are demonstrated to provide accurate basicity data of superbases; however, they have some intrinsic limitations: not all compounds form proton-bound dimers, or due to the limitations of used MS instrumentation. In addition, the kinetics method cannot be used for extending the scaleit needs a reliable scale for its operation. Proton exchange equilibrium measurement methods with FT-ICR-MS instrumentation having electron ionization source with tunable electron energy for in-cell ion generation have proved to provide reasonably accurate data with good agreement with theoretical calculations.4,41,100−105 This work utilizes the previously used GB measurement methodology41,104−106 modified and tested with weaker bases107,108 on an FT-ICRMS instrument specifically developed to overcome the limitations on experiments with superacids and superbases of low volatility we have faced on earlier instruments. Most importantly, the instrument features dual direct insertion probes accessory, which enables independent introduction of the vapors of two low-volatility compounds into the ICR cell for relative basicity measurements. Introducing such compounds has traditionally been the bottleneck of GB measurements.



EXPERIMENTAL AND THEORETICAL METHODS Fourier Transform Ion Cyclotron Resonance Mass Spectrometry Experiments. Detailed description of the FTICR-MS instrumentation and the details of sample introduction and performing the gas-phase proton exchange equilibrium reaction experiments are given in the Supporting Information. The measurements were performed on an in-house modified Varian 930 FT-ICR mass spectrometer with horizontal 7 T superconductive magnet (Santa Clara, CA, USA). Varian Omega 9.1.21 software was used for FT-ICR-MS experiments setup and data acquisition. The temperature of the FT-ICR cell C

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A reaction time during which the proton exchange reaction between two superbasic compounds was followed until the system reached equilibrium. Sequence settings for the collection of mass spectra were quite similar to the previous works on another type of FT-ICR mass spectrometer41,105 and recently developed for the used FT-ICR-MS.107,108 For each reaction time, from 2 to 20 mass spectra were collected and averaged to improve the signal-to-noise ratio. Reaction time after ionization pulse was varied from 100 ms to over 1000 s. Typical radiofrequency sweep excitation and transient signal acquisition mass ranges were from 40 to 1000 m/z units (broadband mode), transient signal collection time was 16.384 or 32.768 ms (256 or 512 K points of data, respectively) and analog-to-digital conversion (ADC) rate was 16 MHz. These settings were appropriate for working throughout overall pressure range of 1 × 10−8 mbar and yielded sufficient signalto-noise ratio and resolution for quantitatively studying gasphase ion−molecule reactions on the used FT-ICR mass spectrometer. All experiment voltages, pulse lengths, delay times, and data acquisition settings were optimized, if needed, to give (a) reasonably strong transient signal and (b) ratio of the intensities of the analytical ion (BH+) and its +1 isotope for each protonated superbase in mass spectrum to match the theoretical isotope distribution at particular data acquisition and resolution settings. Double resonance (ion selection) experiments were used for checking the reversibility of proton exchange reaction and for confirmation that true proton exchange equilibrium of each pair of bases took place. As criteria of equilibrium state of proton exchange reaction the following three conditions had to be fulfilled at the same time: (a) within a series of proton transfer experiments at various reaction times at a particular pressure ratio of neutrals at least on two consecutive reaction times the ΔΔGb was not changing by more than 0.1 kcal mol−1, (b) double resonance experiments confirmed the reversibility of proton transfer reaction, and (c) the observed ΔΔGb was independent (within 0.3 kcal mol−1) of pressure ratio of neutrals. Partial Pressure Ratios of Neutral Compounds. Ratios of p(B1)/p(B2) (eq 5) were obtained from mass spectrometric data by using similar approach as was used in ref 102. The initial ionization by electron beam pulse inside the FT-ICR cell will form molecular ion species B+ by electron removal and some amount of smaller reactive fragments FmH+ from neutral bases molecules B1 and B2. Some of these fragments are very acidic and are reacting very fast with mixture of neutral B1 and B2 gases and are acting as protonating reagents of chemical ionization process, and as a result B1H+ and B2H+ are formed. As B1 and B2 are very strong bases then the assumption is made that strong acidic reactive fragment species will protonate the similar neutral base molecules purely statistically upon collisions. The higher the number of the molecules of the particular base per volume unit inside the FT-ICR cell, the higher amount of corresponding BH+ is formed within first hundreds of milliseconds. The bulky superbasic compounds used in this work exchange protons between themselves slowly under given pressures, and the proton exchange reaction between these two superbasic compounds will become dominant and observable in different reaction time mass spectra after first seconds. Thus, it was generally observed that the ratio X (eq 6) of sums of intensities of molecular cation I(B+) and protonated form I(BH+) of the bases B1 and B2, respectively, has a plateau at reaction times somewhere in the range from 500 to 2000 ms (inset of Figure 1).

Figure 1. Example of the change of ratio of protonated species (compounds 26 and 29) during proton exchange equilibration reaction. The system has reached proton exchange equilibrium at reaction times over 320 s. (inset) The short reaction time region with respective intensities for partial pressure ratio estimation term (eq 6).

The numerical value of this plateau was taken as approximation of molecule size-uncorrected partial pressure (pu) ratios. pu (B1) pu (B2)

≈X=

I(B1H+) + I(B1+) I(B2H+) + I(B2+)

(6)

These intensities were corrected for theoretical isotopic distribution according to IUPAC recommendations of isotope abundances,109 also B+ (if present) + 1 isotope overlapping with BH+ was taken into account. Ion gauge sensitivity corrections110 were applied as reported earlier,41,102,106,111 by using hot filament ion gauge sensitivity correction equation from ref 113 and polarizabilities from ref 114 to give sizecorrected partial-pressure ratio of neutrals for eq 5. The ratio was not extrapolated to zero reaction time (as was done for neutral acids in ref 102 for the data obtained on different instrument system), as plateau on this ratio was always observed somewhere between 500 and 2000 ms, and application of current approach yielded good agreement with earlier GB data of compounds (see below). Besides avoiding inaccuracies arising from the distant location of the pressure gauge, an additional advantage of this pressure ratio estimation approach is that it is rather insensitive to partial decomposition of the compounds (which has a dramatic effect on measurement accuracy if pressure gauge readings are used for estimating partial pressures)the partial pressure ratio determined from intensities of molecular cation and protonated form of the bases B1 and B2 is independent of the presence of low molecular mass decomposition products. The drift of this pressure ratio during experiment series was taken into account when necessary, by interpolating the pressure ratio change in time for data point of particular reaction time from the short reaction time mass spectra taken at the beginning and at the end of reaction series. Intensity values of protonated species I(B1H+) and I(B2H+) for calculating intensity ratio in eq 5 were taken from mass spectra, and these were corrected with correction factor for taking into account theoretical isotope distribution as described above. The FT-ICR-MS instrument was not used in its typical D

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Table 1. Experimental and Calculated Gas-Phase Basicities, Basicities in Tetrahydrofuran and Acetonitrile of Superbases from This Work and Literature

a

Value from this work if not noted otherwise. pKa values of TPPN in THF and MeCN are not included in Table 1; these are 24.1 and 32.10 pKa units, respectively. bKilocalories per mole (1 kcal mol−1 = 4.184 kJ mol−1). cBasicity calculations at DFT B3LYP/6-311+G** level if not noted otherwise.118 dReference 41. eReference 61. fReference 80. gReference 41. hReference 50. iReference 104. jReference 73. kReference 5. lReference 60. mTheoretical value. nReference 121. oReference 122. pReference 56. qReference 123. rReference 4, anchor compound. sReference 84. tReference 116. uReference 124. vCalculations by the DFT BP TZVP method (using Turbomole 6.4). xReference 90. yReference 52. E

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

every compound against at least two others at various pressure ratios. The ΔΔGb measurement data and the assigned GB values from this work and ref 41 above the phosphazene base PhP1(dma)3 are presented in Table 1. The absolute GBexp values were obtained from the overlapping relative measurement data by using similar minimization approach as in previous work41 and keeping the reference values of the compounds 26, 30, 37, and 39 (serving as anchor compounds for the values of this work) from the previous work41 constant. The whole scale is anchored to the GB value of MTBD from NIST database.4 The consistency standard deviation of the new GB values on the ladder is 1.8 kcal mol−1 similar to the previous work.41 Theoretical GB values from this work and literature for the bases under study are presented in Table 1. pKa Measurements. The pKa values in THF and MeCN are given in Table 1. Relative pKa measurement data in THF and MeCN are presented in Supporting Information Table S4. Absolute pKa values in THF and in MeCN were obtained by using similar data treatment and minimization approach as in previous works.5,104 The assigned pKa values from this work and from literature are presented in Tables 1 and S4.

high-resolution mode (R > 100.000), which requires ultrahigh vacuum (pressure in 1 × 10−10 mbar range). Instead it was used in low-resolution mode (R ≈ 1.200−3.000) for having good intensity, sensitivity and proper isotope ratios. If necessary, the intensities of B1H+ and B2H+ were corrected for intensities of molecular ion +1 isotopes I(B1+) and I(B2+). Generally, the molecular ions disappeared on reaction times over a couple of seconds. This approach of approximation of partial pressure ratios has been validated for acids in ref 102. For validation of the approach for bases we can compare the GB measurement results from this work and ref 41, where the partial pressure ratios for ΔΔGb calculations were estimated from ion gauge readings. The difference of basicities of PhP2(dma)5 and 2-ClC6H4P2(dma)5 in ref 41 was 1.2 kcal mol−1, in the present work (see Table 1) over 2-CF3-C6H4P2(pyrr)5 the same difference is 1.17 kcal mol−1. The difference of basicities of 2-tertbutylimino-2-diethylamino-1,3-dimethylperhydro-1,3,2-diazaphosphorine (BEMP) and 4-MeO-C6H4P1(pyrr)3 in ref 41 was 1.1 kcal mol −1 , in the present work over 1,8-bis(tetramethylguanidino)naphthalene (TMGN) it is 1.12 kcal mol−1; over TMGN and PhP1(dma)2tmg it is 0.94 kcal mol−1. This good agreement of the results from this work and previous work, where partial pressures ratios of the neutrals were obtained from ion gauge readings, gives us confidence that the results from present work are consistent with our earlier data.41 Ultraviolet−Vis Spectrophotometric pKa Measurements in Tetrahydrofuran and Acetonitrile. Earlier instrumental setup and experimental methodology was used.5,104,115−117 See the Supporting Information for pKa measurements, experiment descriptions, and data treatment. Chemicals. Synthesis and purification of free phosphazene bases is described earlier.5,50,104 Commercial solution of tBuP2(dma)5 in tetrahydrofuran (THF; 2 mol/L) was kept under vacuum (1 × 10−3 mbar) for removal of solvent. Commercial BEMP was used as received. HP1(pyrr-2-CH2−Npyrr)3 was the same as in ref 80, CH2(TBD)2 was the same as in ref 60, and TMGN was the same as in ref 50. TMPN, TiPrPN, TBPN, TcyPPN, HMPN, and 1,8-bis(trispyrrolidinophosphazenyl)naphthalene (TPPN) were the same as in refs 61 and 62, and free bases were prepared from their salts when needed. Purification of solvents for UV−vis pKa measurements, used reference bases and acidic and basic titrants were described earlier.5,50,104,116 Theoretical Methods of Gas-Phase Basicity Calculations. The quantum chemical computations reported in this work were mostly performed using Gaussian09 series of programs.118 Density functional theory (DFT) calculations were performed using the B3LYP hybrid functional. Full geometry optimization and vibrational analyses were performed using the 6-311+G** basis set. All stationary points were found to be true minima (Nimag = 0). Unscaled B3LYP/6-311+G** frequencies were used to calculate the GB values of the neutral bases, taking into account the zero-point frequencies, finite temperature 0−298 K correction, pressure−volume work term, and the entropy term as appropriate. For some molecules also the GB calculations were made using Turbomole 6.4119,120 calculations using TZVP basis set.



DISCUSSION Agreement of GBexp and GBcalc Data. Table 1 presents besides the experimental data also the theoretical GB values from this work and from the literature. A recent100 analysis showed that for smaller bases (including also some superbases: two P1 phosphazenes and MTBD and almost superbasic 1,5diazabicyclo[4.3.0]non-5-ene (DBN) the correlations of the experimental and theoretical GB values from most of the Gaussian-type, complete basis set, and DFT methods give close to unity slope; the standard deviation of correlation is generally in the range from 2.9 to 4 kcal mol−1, and the average absolute error is in the range from 1.8 to 5 kcal mol−1. The majority of the results from this work are in line with these findings (see Table S3 and Figure S1 in Supporting Information): up to the experimental GB value of 272 kcal mol−1 the agreement between theory and computations is very good: the root-meansquare (RMS) difference is 1.9 kcal mol−1. The biggest discrepancies are observed with the methyl-substituted Verkade’s base and TMGN. In both cases the experimental values are higher by 4.1 and 4.9 kcal mol−1, respectively, which are both still tolerable. However, with the three strongest basesEtP1(tmg)3, t-BuP1(tmg)3, and PhP3(pyrr)7 (experimental GB values up from 273 kcal mol−1)there is a systematic discrepancy with the RMS difference of 7.1 kcal mol−1, and importantly, all three bases have been found much stronger by computations than by experiments. The root of the inconsistency seems to be in the disagreement of experimental and theoretical results of the third-order homologues of both studied phenylphosphazene homological rowsPhP3(dma)7 and PhP3(pyrr)7see Scheme 1. Addition of the second phosphazene unit to the respective P1 phosphazenes (52 and 43, respectively) leads to almost the same basicity increase, as shown by measurements and predicted by calculations: 15.6 and 13.9 kcal mol−1 for 52 and 13.8 and 14.4 kcal mol−1 for 43. The effect of addition of the third phosphazene unit, however, is predicted very differently by experiments and computations. For both homological rows the experiment yields to nonadditive increase, 7.9 and 7.5 kcal mol−1. This is analogous to the change in solution basicities and corresponds to the general nonadditive pattern of basicity increase upon increase of molecular size.52 However, calculations predict additive basicity



RESULTS Gas-Phase Basicity Measurements. The experimental GB values of altogether 28 superbases were successfully obtained by measuring relative basicity (ΔΔGb value) of F

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Scheme 1. Structure−Basicity Relationshipsa of Phosphazenes and Guanidinophosphazenes in Gas Phase

a GB values (below compounds) and effects of substitutions (below of or next to the arrows) are in kilocalories per mole. Non-bracketed values are experimental, and values in brackets are from theoretical calculations. Values from this work and from refs 4, 41, 52, 127, and 125.

increase by up to 14.5 kcal mol−1 for PhP3(dma)7 and more than additive basicity increase of 17.9 kcal mol−1 for PhP3(pyrr)7. This disagreement of experiment and calculations is higher than generally observed and accepted100 and demonstrates the need for further examination of these compounds both experimentally and theoretically (e.g., computations at higher level of theory). This work together with experimental and theoretical (B3LYP/6-311+G** or levels of similar accuracy) GB data from refs 4, 41, 42, 100, 105, and 127 forms a set of 167 bases and covers GB range over 78 kcal mol−1 (57 orders of magnitude). The set includes various types of bases ranging from ammonia to the strongest bases of this work including amines, diamines, guanidines, anilines, pyridines, phosphines, etc. Correlation of experimental and theoretical GB values is as follows: a = 1.02(0.01), b = −4.9(1.9) kcal mol−1, R2 = 0.989, S = 2.25 kcal mol−1. These results demonstrate that theoretical calculations and experiments give consistent GB values starting from below 200 kcal mol−1 to up to 270 kcal mol−1 for all of these different base families, including the phosphazene and proton sponge type of superbases. The results from this work fit well into the general picture of correlation of previously published gas-phase basicities above 190 kcal mol−1.

Gas-Phase Basicity Trends in Homologation of P1 to P2 to P3 of Phosphazenes and Substituent Effects. The main structure−basicity relationships of phosphazenes and guanidinophosphazenes are summarized in Scheme 1. As seen above the general agreement of experimental and calculated effect of structure change on basicity is very good, as long as medium-size bases are considered. For P3 phosphazenes the experimental and calculated effects differ more. Calculations suggest that replacement of the second amino (dma, pyrr) group with iminophosphorane group, yielding compounds tBuP3(dma)7, PhP3(dma)7, and PhP3(pyrr)7, is close to additive. At the same time experiment shows nonadditive behavior when moving already from phenyl P2-phosphazenes to P3-phosphazenes. That is, the compression of basicity increase is observed for tris-substituted homologues. In refs 84 and 52 theoretical GB values were obtained also for higher members of several homological series of phosphazenes. In both works significant compression of basicity increase was observed for analogous phosphazene series starting from P4-phosphazenes, and in ref 52 the model of exponential decay of basicity increase in those homological series was presented and validated. For the tBuP4(dma)9 our calculations at B3LYP/6-311+G** level yield GB value of 287.8 kcal mol−1; this agrees very well with the G

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Table 2. Statistical Analysis of Gas-Phase Basicity of Substituted Bases with Field/Inductive (σF), Resonance (σR−), and Polarizability (σα) Substituent Constantsa series b

4-X-DA 4-X-C6H4P1(pyrr)3c 4-X-C6H4P2(pyrr)5d

GB0

s(GB0)

ρF

s(ρF)

ρR−

s(ρR−)

ρα

s(ρα)

S

r2

n

217.82 251.85 265.53

0.43 0.41 0.17

−17.26 −12.95 −9.02

0.84 1.51 0.47

−20.87 −22.11 −13.53

1.46 1.73 0.68

−6.33 −2.23

1.13 1.33

0.58 0.43 0.17

0.993 0.998 0.999

10 5 4

σ values are taken from Table 9 of ref 126. bReference 127, only experimental data of the bases that protonate on NMe2 are included. cData from ref 41. dWithout polarizability term σα. a

literature values (GB = 287.5 and 287.7 kcal mol−1 in refs 125 and 52, respectively). HP1(pyrr-2-CH2−N-pyrr)3 is the strongest amino-substituted P1 phosphazene base in the gas phase. It is by 11.9 kcal mol−1 or 8.7 orders of magnitude (in log K units) more basic than its simpler analogue HP1(pyrr)3. At the same time in MeCN the basicity difference of these two compounds is only 2.3 pKa units. This rather high basicity strengthening effect in the gas phase is most probably caused by arrangement of the highly polarizable pyrr-2-CH2−N-pyrr groups around the protonation center of the protonated form. Comparison of para-ring-substituted phosphazene P1, P2, and P3 homologues reveals that in the gas phase similar gradual saturation of substituent effect in the higher homologues is observed as was reported in ref 127 for the same substituents in acetonitrile and THF. In the gas phase 4-CF3−C6H4P2(pyrr)5 is by 4.8 kcal mol−1 weaker base than PhP2(pyrr)5; the effect of the same substituent in P1 phosphazenes was found 6.3 kcal mol−1 in ref 41. The decrease of substituent effect of 4-MeO substitution is even larger, by 3.2 and 1.5 kcal mol−1 compared to nonsubstituted P1 and P2 homologues, respectively. Further sensitivity decrease upon moving toward higher homologues is observed from the GB data of the same substitution in the PhP3(dma)7 series. 4-CF3−C6H4P3(dma)7 is only by 3.2 kcal mol−1 weaker, and 4-MeO-C6H4P3(dma)7 is only by 1.3 kcal mol−1 stronger base than PhP3(dma)7. Analysis of the experimental GB data of para-ring-substituted phosphazenes with Taft/Topsom substituent constants σF, σR, and σα126 according to eq 7 GB = GB0 + ρF σF + ρR− σR− + ρα σα

groups in phosphazenes. The results of this work show that substitution of the first N,N-dimethylamino group in PhP1(dma)3 (GB = 246.1 kcal mol−1; ref 41) with a tmg unit to form PhP1(dma)2tmg increases its GB by 10.4 kcal mol−1. Replacement of second and third amino units increases product basicity by additional 9.4 and 4.0 kcal mol−1, respectively. The energetic effect of introduction of the second tmg is similar to the first, but with the third tmg unit strong saturation of the effect is observed. The agreement of the net effect of this consecutive substitution from experimental and theoretical data is excellent, 23.8 and 23.5 kcal mol−1, respectively. It is on the same order to the same substitution sequence also when starting from HP1(dma)3 and t-BuP1(dma)3, where the net effect of exchanging the three substituents from the experimental data is 22.7 and 20.7 kcal mol−1, respectively. Theoretical calculations show for these substitution series somewhat higher net effect, more than 26 and 27 kcal mol−1, respectively. Very small basicity increase, 2.2 kcal mol−1, is observed experimentally upon replacement of the third amino group NEt2 with tmg group in t-BuP1(tmg)2NEt2 to give tBuP1(tmg)3. Superbases with the −NP(tmg)3 fragment are more volatile than aryl-P3-phosphazenes. The tmg-substituted phosphazene derivatives are very promising as stronger and relatively volatile superbases for further extension of the experimental GB scale. The unexpected order of basicities of the HP1(tmg)3, t-BuP1(tmg)3, and EtP1(tmg)3 according to our experimental and theoretical basicity data in gas phase is an example of some nonconsistency of experiment and calculations. The experimentally observed GB order is supported by the same basicity order of these compounds in the nonpolar solvent THF50 and by the analogue with experimentally observed order in gas phase as well as in acetonitrile for HP1(pyrr)3, t-BuP1(pyrr)3, and EtP1(pyrr)3.41 Two possible reasons for the stronger basicity-enhancing effect of Et versus tBu substituent are (1) higher steric strain in the case of t-Bu group than with Et group and (2) partial stabilization of the protonated EtP1(tmg)3 via cationic hyperconjugation. These compounds can be used as useful sample compounds for further testing various experimental and calculation methods with simple superbasic molecules. Future work in extending and supplementing the upper part of the experimental GB scale with other compounds may also bring out possible experimental problems of some of these compounds. The volatility and evaporation rate of PhP3(pyrr)7 and its phenyl-substituted derivatives is already very low, even when introduced from the heated solid probe. This suggests that further homologation of aryl-Pn(pyrr)2n+1 phosphazenes is probably not a suitable path for preparing stronger bases for further stepwise extension of the experimental GB scale by small basicity differences. Homologation of R-Pn(dma)2n+1 phosphazenes, guanidinophosphazenes, variation of amino and guanidino groups linked to P atoms of phosphazenes and

(7)

are presented in Table 2. It is observed that correlations are excellent. For both phosphazene series the resonance term has highest contribution, 60% of the total budget, and the field/inductive term contributes 35 to 40%. Polarizability term of the substituent has small contribution only in the case of 4-X-C6H4P1(pyrr)3. For both series the resonance term is on the same order as it was observed for substituted para-substituted N,N-dimethylanilines (DA).127 The field/inductive and polarizability terms are of lower contribution in the phosphazene series than in the reference series. In the 4-X-C6H4P2(pyrr)5 series the sensitivity toward para-substituent parameters decreases. The resonance term contribution fades by more than 1.61 times, and field/ inductive term decreases by 1.56 times, if compared to the 4-XC6H4P1(pyrr)3 series. This can be regarded as evidence of decrease of ylidic structure contribution (as opposed to the ylenic structure)128 in the neutral phosphazene. Effect of Stepwise Substitution of Dimethylamino Units with −N=C(NMe2)2 Units in P1 Phosphazenes. Previosly50 it was shown experimentally (in solution) and theoretically (in the gas phase) that −N=C(NMe2)2 (tmg) groups have much higher base strengthening effect than dma H

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Table 3. Correlation of Experimental Gas-Phase Basicity with Respective Experimental Basicity Values pKa(X) in Solution According to Equation 8 pK(X)

compounds

a

s(a)

b

s(b)

S

R2

n

a/1.364a

pKa (THF)

all basesb bases measured in this work 4-X-C6H4P1(pyrr)3c 4-X-C6H4P2(pyrr)5 4-X-C6H4P3(dma)7 4-X-N,N-dimethylanilined all basesb bases measured in this work 4-X-C6H4P1(pyrr)3c 4-X-C6H4P2(pyrr)5 4-X-C6H4P3(dma)7 4-X-N,N-dimethylanilined

2.72 1.07 4.38 2.89 1.91 3.07 2.72 1.44 3.16 2.16 1.49 3.12

0.13 0.19 0.15 0.17 0.12 0.39 0.11 0.22 0.18 0.05 0.18 0.34

201.6 243.2 181.5 205.0 224.9 201.9 181.3 224.1 181.7 206.0 223.0 181.3

2.2 3.5 2.4 3.5 3.4 1.7 2.5 6.5 4.0 1.4 5.5 3.7

8.8 3.2 0.3 0.2 0.2 1.6 9.6 2.6 0.5 0.1 0.4 1.3

0.812 0.566 0.997 0.997 0.994 0.938 0.807 0.658 0.990 0.999 0.986 0.954

98 25 5 3 3 6 137 24 5 4 3 6

1.99 0.78 3.21 2.12 1.40 2.25 2.00 1.06 2.32 1.58 1.09 2.29

pKa(MeCN)

a The attenuation factor calculated from the slope is a/2.30RT = a/1.364. bAll bases from this work and refs 41, 4, 127, 42, 73, and 129. cReference 41; GB values are taken in kilocalories per mole. dAnalysis of the data from ref 127, only N-bases are included.

Replacement of −NMe2 groups in DMAN with −NP(pyrr)3 to give TPPN will increase the proton sponge basicity even more: in THF by 13 pKa units and in MeCN by 13.5 pKa units (see Table S2 in Supporting Information). Previously, some of us61,62,56 determined pKa values in MeCN for the bisphosphazene proton sponges by 31P NMR titration method and by theoretical calculations. In this work we report pKa values in MeCN and THF, which are obtained with UV−vis spectrophotometric titration method5,50,104,116 (see Table 1 and Table S2 in Supporting Information). The present results in MeCN are in excellent agreement with the earlier experimental NMR data except for TMPN, which is in this work determined to be ca. 1 pKa units stronger base than was observed previously.61 This can be explained by the experimental difficulties reported in the previous work (peak broadening and low solubility of TMPN for the 31P NMR experiments). In this work we used UV−vis spectrophotometric titration at much lower concentration (n × 10−5 M), and true equilibrium was reached during the titration experiments. The recommended pKa(MeCN) value of TMPN is 30.35 obtained in this work. Theoretical calculations (isodensity polarized continuum model method)61 underestimate the pKa of TMPN by 0.7 pKa units, overestimate the pKa of TiPPN by 0.5 pKa units and TcyPPN by 1.7 pKa units, respectively. For TBPN calculations underestimate its pKa by 1.1 pKa units. For TPPN we determined its pKa(MeCN) as 32.1. This is 0.2 pKa units lower than previous experimental value (32.3)62 and 0.9 pKa units lower than obtained by theoretical calculations (33.0).62 The general observation is that theoretical calculations will give more accurate pKa values for bisphosphazene proton sponges with smaller substituents. Superbasic bis-guanidino proton sponge CH2(TBD)2,60 which is 1 to 2 pKa units weaker base in MeCN solution than all the discussed bisphosphazene proton sponges, appears to be in gas phase stronger than TMPN. Correlations between Gas-Phase Basicity and Solution Basicities. The data analysis results for various sets of compounds from this work and literature are presented in Table 3. The correlation equation has the following form

proton sponges will be a more promising option for preparing compounds for further extending the experimental GB scale. Effect of Substituents on Proton Sponges Based on 1,8-Disubstituted Naphthalene. For the first time the experimental GB values and pKa values in THF for six superbasic 1,8-disubstituted naphthalene proton sponges were determined in this work. Earlier predictions of the superbasic nature of TMGN in the gas phase from theoretical calculations54,90 were confirmed by experiment. Our GB experiments yield by 4.9 kcal mol−1 higher basicity for TMGN than was calculated previously.90 TMGN is by ca. 18 kcal mol−1 stronger base than the original DMAN proton sponge (GB = 238.0 kcal mol−1).4 Our experiments in MeCN confirm the previously found experimental value.123 Our experimental results of bisphosphazeno 1,8-naphthalene proton sponges can be compared with the earlier theoretical results56,61 in the gas phase. The order of the effect of the alkyl group R (Me- < i-Pr- < n-Bu- < cy-pentyl) in −NPR 3 of bisphosphazene base GB agrees with earlier calculations; however, our results yield from 1.05 to 1.14 times lower gain in energy units when substituting Me group with bulkier alkyl groups. The present experimental results confirm earlier findings that basicity of the proton sponges in the gas phase is mainly determined by the basicity of the corresponding −NPR3 fragments and not by the steric effects of the R substituents.44 In solution our results display a different picture from the gas phase. In THF solution the replacement of Me groups in TMPN with linear n-Bu groups leads to basicity increase by 0.3 pKa units (TBPN), and the replacement with bulkier i-Pr and cyclopentyl groups makes the base weaker by 1.0 (TiPrPN) and 0.9 (TcyPPN) pKa units, respectively. This is due to the hindered access of the THF molecules to stabilize the protonated proton sponges with bulkier substituents (TiPrPN, TcyPPN) via hydrogen bonding. In MeCN solution, which has much lower hydrogen bond acceptor ability, this effect of hindrance is less pronounced. Only TcyPPN deviates from the gas-phase order of basicity increase. Our experiments and calculations confirm earlier theoretical finding56 that the effect of replacement of two −NMe2 groups in DMAN (GB = 238.0 kcal mol−1, pKa(THF) = 11.1, pKa(MeCN) = 24.1)4,5,116 by two −NP(dma)3 units, to form bisphosphazene proton sponge HMPN introduces increase of basicity, in total by ca. 31 kcal mol−1. In THF the effect is also enormous, 10.8 pKa units and in MeCN 11.3 pKa units.

GBexp = a pK a(X) + b

(8)

where pKa(X) denotes basicity value in a solvent X, a and b are the slope and intercept, respectively. I

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Previously41 it was found that correlation of basicities in the gas phase and in solution of compounds from various base families is poor due to the varied influence of solvent and other factors on the basicities of bases. Correlations including all experimental data available in a particular solvent and in the gas phase yield an attenuation factor close to 2, with low correlation coefficient (R2 = 0.81) for both solvents. Parameters of the correlation of experimental basicity data in the gas phase and in solutions of the set of compounds from this work could not be used for analyzing solvent influence because of poor correlation. However, the picture is much clearer when trends within compound families are examined. Excellent correlations of basicities in the gas phase and in solution were previously observed41 for the para-substituted aryl P1 phosphazene family: 4-X-C6H4P1(pyrr)3, with attenuation factor of 3.21 for gas phase versus THF and 2.32 for gas phase versus MeCN, see Table 3. Data from ref 127 yield for para-substituted DA-s attenuation factor of 2.25 and 2.29 for gas phase versus THF and gas phase versus MeCN, respectively. While, for both DA and P1 phosphazene series the attenuation factors of MeCN are similar, then attenuation factors of THF for these families are very different. When moving toward higher homological series of para-substituted phosphazenes, then both solvents become increasingly better in differentiating basicities. For 4-XC6H4P3(dma)7 MeCN almost reaches the gas phase by its differentiating ability of basicities in these series. In the gas phase t-BuP2(dma)5 is by 4.2 kcal mol−1 stronger base than EtP2(dma)5. In solution the order is reversed; in THF the former base is by 0.4 pKa units and in MeCN 0.55 pKa units weaker base than the latter. For the same alkyl substituents reversed order was previously obtained for RP1(pyrr)3 pair in gas phase41 and in both of the solvents.104,5 Future Perspective of Extending Gas-Phase Basicity Scale and Use of Obtained Data. As was discussed above the aryl-P3 phosphazenes, which were found to be very useful in building basicity scales in MeCN5 and THF,116,104,50 may become unusable for future extension of the GB scale by using FT-ICR-MS equilibrium method used in the present work because of their low volatility and the dramatic basicityreducing effect of the aromatic ring attached to the imino nitrogen. More promising families are alkyl-P 3 and P 4 phosphazenes and P1 and P2 guanidinophosphazenes.50 There are several strategies that yield superbases stronger than EtP1(tmg)3, being at the same time stable and having relatively low molecular mass, for example, bisphosphazene proton sponges,61 BIG bases,75 etc. Another direction would be measurement of the GB of the bases of practical interest, such as polyguanidine bases,68 chiral phosphazenes,79 stable carbenes, etc. The present results will help making justified decisions when choosing or developing superbases having appropriate properties for base-catalyzed reactions and for preparing new materials. Previously some PILs have been prepared from commercial P1 phosphazene, amidine, and guanidine bases.15,18 Some of these showed exceptionally good thermal stability and are expected to have potential use in electrochemical energy storage and transformation devices,21,130−132 in synthesis and catalysis,27 etc. Several of the compounds used in the present work can be used for preparation of PILs. These candidates have simultaneously high basicity, low molecular weight, and good stability toward ambient atmosphere in protonated form. Superbases have several advantages over common weaker bases (amines, imidazoles, triazoles, etc.) used for preparation of PIL

and organic ionic plastic crystals (OIPC). Use of superbases enables using of weaker and thus less hazardous acids (ΔpKa value higher than 16 pKa units in water has been recommended as the threshold of suitable acidity difference of the acid and the conjugate acid of the base20) for preparation of PIL or OIPC with good stability and electrochemical properties. Superbasic compounds can be used also for extending other affinity scales: lithium cation affinity,133 etc and their GB data is useful for testing various theoretical calculation methods.134



CONCLUSIONS Experimental GB scale was extended to cover basicities of P2 and P3 phosphazenes, P1 guanidinophosphazenes, and bisphosphazeno proton sponges. This was made possible by using dual solid probe accessory for introducing vapors of lowvolatility superbases into FT-ICR-MS. The results of this work confirm the good agreement of experimentally obtained and theoretically (at DFT B3LYP/6-311+G** level or on TZVP basis set) calculated GB values of P2 and P3 phosphazene bases and 1,8-phosphazeno-substituted naphthalene proton sponges. However, for stronger guanidino-substituted phosphazenes the agreement of experimental and theoretical GB values is not good. Structure−basicity dependence in gas phase and solvent influence on the basicity of several superbase families were experimentally evidenced. The results form ground for future experimental work in the field and can be used as a tool for designing new superbases and their applications.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpca.6b01552. Details of FT-ICR-MS instrumentation; GB and pKa experiments and calculations; GB data; complete citation of refs 101, 102, 118, 122. (PDF)



AUTHOR INFORMATION

Corresponding Author

*Phone: +372 7375274. Fax: +372 7375264. E-mail: ivari. [email protected]. Present Address §

FB2 Biologie/Chemie, Universität Bremen, Leobener Str. im NW2, 28359 Bremen, Germany. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the Estonian National R & D infrastructure development program of Measure 2.3 “Promotion of development activities and innovation” (Regulation No. 34) funded by the Enterprise Estonia Foundation, by Grant No. 8689 from the Estonian Research Council, by the EU through the European Regional Development Fund (TK141 “Advanced materials and high-technology devices for energy recuperation systems”, 2014-2020.4.01.15-0011), and by the institutional research grant of Ministry of Education and Research of Estonia IUT20-14 (TLOKT14014I). Financial support “Chemiefonds Fellowship for Ph.D. candidates“ by the Fonds der Chemischen Industrie is gratefully acknowledged (doctoral scholarship for J.F.K). J

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A



(20) Miran, M. S.; Yasuda, T.; Susan, M. A. B. H.; Dokko, K.; Watanabe, M. Electrochemical Properties of Protic Ionic Liquids: Correlation between Open Circuit Potential for H2/O2 Cells under Non-Humidified Conditions and ΔpKa. RSC Adv. 2013, 3, 4141− 4144. (21) MacFarlane, D. R.; Tachikawa, N.; Forsyth, M.; Pringle, J. M.; Howlett, P. C.; Elliott, G. D.; Davis, J. H., Jr.; Watanabe, M.; Simon, P.; Angell, C. A. Energy Applications of Ionic Liquids. Energy Environ. Sci. 2014, 7, 232−250. (22) Sun, X.; Luo, H.; Dai, S. Ionic Liquids-Based Extraction: A Promising Strategy for the Advanced Nuclear Fuel Cycle. Chem. Rev. 2012, 112, 2100−2128. (23) Wang, C.; Luo, H.; Jiang, D.; Li, H.; Dai, S. Carbon Dioxide Capture by Superbase-Derived Protic Ionic Liquids. Angew. Chem., Int. Ed. 2010, 49, 5978−5981. (24) Yang, Z.-Z.; Zhao, Y.-N.; He, L.-N. CO2 Chemistry: Taskspecific Ionic Liquids for CO2 Capture/Activation and Subsequent Conversion. RSC Adv. 2011, 1, 545−567. (25) Sixta, H.; Michud, A.; Hauru, L.; Asaadi, S.; Ma, Y.; King, A. W. T.; Kilpeläinen, I.; Hummel, M. Ioncell-F: A High-strength Regenerated Cellulose Fibre. Nord. Pulp Pap. Res. J. 2015, 30, 43−57. (26) Michud, A.; King, A.; Parviainen, A.; Sixta, H.; Hauru, L.; Hummel, M.; Kilpeläinen, I. Process for the Production of Shaped Cellulose Articles, FI Patent Application PCT/FI2014/050238, 2014. (27) Hallett, J. P.; Welton, T. Room-Temperature Ionic Liquids: Solvents for Synthesis and Catalysis. 2. Chem. Rev. 2011, 111, 3508− 3576. (28) Hollóczki, O.; Nyulászi, L. Carbenes from Ionic Liquids. In Electronic Effects in Organic Chemistry; Kirchner, B., Ed. Topics in Current Chemistry Series 351; Springer: Berlin, Germany, 2014; pp 1−24. (29) Hollóczki, O.; Gerhard, D.; Massone, K.; Szarvas, L.; Németh, B.; Veszprémi, T.; Nyulászi, L. Carbenes in Ionic Liquids. New J. Chem. 2010, 34, 3004−3009. (30) Jørgensen, T. J. D.; Bojesen, G.; Rahbek-Nielsen, H. The Proton Affinities of Seven Matrix-assisted Laser Desorption/Ionization Matrices Correlated with the Formation of Multiply Charged Ions. Eur. Mass Spectrom. 1998, 4, 39−45. (31) Zenobi, R.; Knochenmuss, R. Ion Formation in MALDI Mass Spectrometry. Mass Spectrom. Rev. 1998, 17, 337−366. (32) Barylyuk, K.; Fritsche, L.; Balabin, R. M.; Nieckarz, R.; Zenobi, R. Gas-Phase Basicity of Several Common MALDI Matrices Measured by a Simple Experimental Approach. RSC Adv. 2012, 2, 1962−1969. (33) Knochenmuss, R.; Zenobi, R. MALDI Ionization: The Role of In-Plume Processes. Chem. Rev. 2003, 103, 441−452. (34) Grechnikov, A. A.; Borodkov, A. S.; Alimpiev, S. S.; Nikiforov, S. M.; Simanovsky, Y. O. Gas-Phase Basicity: Parameter Determining the Efficiency of Laser Desorption/Ionization from Silicon Surfaces. J. Anal. Chem. 2013, 68, 19−26. (35) Jaskolla, T. W.; Karas, M. Compelling Evidence for Lucky Survivor and Gas Phase Protonation: The Unified MALDI Analyte Protonation Mechanism. J. Am. Soc. Mass Spectrom. 2011, 22, 976− 988. (36) Oss, M.; Kruve, A.; Herodes, K.; Leito, I. Electrospray Ionization Efficiency Scale of Organic Compounds. Anal. Chem. 2010, 82, 2865− 2872. (37) Ehrmann, B. M.; Henriksen, T.; Cech, N. B. Relative Importance of Basicity in the Gas Phase and in Solution for Determining Selectivity in Electrospray Ionization Mass Spectrometry. J. Am. Soc. Mass Spectrom. 2008, 19, 719−728. (38) Ikonomou, M. G.; Blades, A. T.; Kebarle, P. Investigations of the Electrospray Interface for Liquid Chromatography/Mass Spectrometry. Anal. Chem. 1990, 62, 957−967. (39) Amad, M. H.; Cech, N. B.; Jackson, G. S.; Enke, C. G. Importance of Gas-phase Proton Affinities in Determining the Electrospray Ionization Response for Analytes and Solvents. J. Mass Spectrom. 2000, 35, 784−789. (40) Bagal, D.; Zhang, H.; Schnier, P. D. Gas-Phase Proton-Transfer Chemistry Coupled with TOF Mass Spectrometry and Ion Mobility-

REFERENCES 2,3

(1) Although there are other definitions of superbasicity in use, we consider the limit of superbasicity of DBU (1,8-diazabicyclo[5.4.0] undec-7-ene), GB = 242.7 kcal mol−1,4 pKa(MeCN) = 24.35 due to its higher use in practice than, for example, 1,8-bis(dimethylamino) naphthalene. (2) Superbases for Organic Synthesis: Guanidines, Amidines, Phosphazenes and Related Organocatalysts; Ishikawa, T., Ed.; Wiley: Chichester, U.K., 2009. (3) Busca, G. Bases and Basic Materials in Chemical and Environmental Processes. Liquid versus Solid Basicity. Chem. Rev. 2010, 110, 2217−2249. (4) Hunter, E. P. L.; Lias, S. G. Evaluated Gas Phase Basicities and Proton Affinities of Molecules: An Update. J. Phys. Chem. Ref. Data 1998, 27 (3), 413−656. http://webbook.nist.gov/chemistry/ (accessed March 28, 2016). (5) Kaljurand, I.; Kütt, A.; Sooväli, L.; Rodima, T.; Mäemets, V.; Leito, I.; Koppel, I. A. Extension of the Self-Consistent Spectrophotometric Basicity Scale in Acetonitrile to a Full Span of 28 pKa Units: Unification of Different Basicity Scales. J. Org. Chem. 2005, 70, 1019− 1028. (6) Takeda, T.; Terada, T. Development of a Chiral Bis(guanidino)iminophosphorane as an Uncharged Organosuperbase for the Enantioselective Amination of Ketones. J. Am. Chem. Soc. 2013, 135, 15306−15309. (7) Munshi, M. K.; Gade, S. M.; Mane, M. V.; Mishra, D.; Pal, S.; Vanka, K.; Rane, V. H.; Kelkar, A. A. 1,8-Diazabicyclo[5.4.0]undec-7ene (DBU): A Highly Efficient Catalyst in Glycerol Carbonate Synthesis. J. Mol. Catal. A: Chem. 2014, 391, 144−149. (8) Buitrago Santanilla, A.; Christensen, M.; Campeau, L.-C.; Davies, I. W.; Dreher, S. D. P2Et Phosphazene: A Mild, Functional Group Tolerant Base for Soluble, Room Temperature Pd-Catalyzed C−N, C−O, and C−C Cross-Coupling Reactions. Org. Lett. 2015, 17, 3370− 3373. (9) Krawczyk, H.; Dzięgielewski, M.; Deredas, D.; Albrecht, L.; Albrecht, Ł. Chiral Iminophosphoranes  An Emerging Class of Superbase Organocatalysts. Chem. - Eur. J. 2015, 21, 10268−10277. (10) Coles, M. P. Bicyclic-guanidines, -guanidinates and -guanidinium salts: wide ranging applications from a simple family of molecules. Chem. Commun. 2009, 3659−3676. (11) Glasovac, Z.; Trošelj, P.; Jušinski, I.; Margetić, D.; EckertMaksić, M. Synthesis of Highly Basic Hexasubstituted Biguanides by Environmentally Friendly Methods. Synlett 2013, 24, 2540−2544. (12) Glasovac, Z.; Eckert-Maksić, M. N,N′,N″-Tris-(3-dimethylaminopropyl)-guanidine, the Procedure of Preparation from Carbodiimide and Application in Reactions of Transesterification of Oil. EP 1786763B1, Sep. 16, 2009. (13) Hopkinson, M. N.; Richter, C.; Schedler, M.; Glorius, F. An Overview of N-Heterocyclic Carbenes. Nature 2014, 510, 485−496. (14) de Frémont, P.; Marion, N.; Nolan, S. P. Carbenes: Synthesis, Properties, and Organometallic Chemistry. Coord. Chem. Rev. 2009, 253, 862−892. (15) Luo, H.; Baker, G. A.; Lee, J. S.; Pagni, R. M.; Dai, S. Ultrastable Superbase-Derived Protic Ionic Liquids. J. Phys. Chem. B 2009, 113, 4181−418. (16) Dai, S.; Luo, H.; Baker, G. A. Superbase-Derived Protic Ionic Liquids. U.S. Patent 8524388B2, Sep. 3, 2013. (17) Greaves, T. L.; Drummond, C. J. Protic Ionic Liquids: Evolving Structure − Property Relationships and Expanding Applications. Chem. Rev. 2015, 115, 11379−11448. (18) Miran, M. S.; Kinoshita, H.; Yasuda, T.; Susan, M. A. B. H.; Watanabe, M. Physicochemical Properties Determined by ΔpKa for Protic Ionic Liquids Based on an Organic Super-Strong Base with Various Brønsted Acids. Phys. Chem. Chem. Phys. 2012, 14, 5178− 5186. (19) Kuroha, M.; Gotoh, H.; Miran, M. S.; Yasuda, T.; Watanabe, M.; Sakakibara, K. Proton-Conductivity-Enhancing Ionic Liquid Consisting of Guanidine and Excess Trifluoromethanesulfonic Acid. Chem. Lett. 2014, 43, 649−651. K

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A MS for the Facile Analysis of Poly(ethylene glycols) and PEGylated Polypeptide Conjugates. Anal. Chem. 2008, 80, 2408−2418. (41) Kaljurand, I.; Koppel, I. A.; Kütt, A.; Rõõm, E.-I.; Rodima, T.; Koppel, I.; Mishima, M.; Leito, I. Experimental Gas-Phase Basicity Scale of Superbasic Phosphazenes. J. Phys. Chem. A 2007, 111, 1245− 1250. (42) Glasovac, Z.; Pavošević, F.; Štrukil, V.; Eckert-Maksić, M.; Schlangen, M.; Kretschmer, R. Toward Extension of the Gas-Phase Basicity Scale by Novel Pyridine Containing Guanidines. Int. J. Mass Spectrom. 2013, 354−355, 113−122. (43) Glasovac, Z.; Štrukil, V.; Eckert-Maksić, M.; Schröder, D.; Kaczorowska, M.; Schwarz, H. Gas-Phase Proton Affinities of Guanidines with Heteroalkyl Side Chains. Int. J. Mass Spectrom. 2008, 270, 39−46. (44) Chen, M.; Moerdyk, J. P.; Blake, G. A.; Bielawski, C. W.; Lee, J. K. Assessing the Proton Affinities of N,N′-Diamidocarbenes. J. Org. Chem. 2013, 78, 10452−10458. (45) Liu, M.; Chen, M.; Zhang, S.; Yang, I.; Buckley, B.; Lee, J. K. Reactivity of Carbene•Phosphine Dimers: Proton Affinity Revisited. J. Phys. Org. Chem. 2011, 24, 929−936. (46) Polyakova, S. M.; Kunetskiy, R. A.; Schröder, D. Proton Affinities of 2-Iminoimidazolines with Bulky N-Alkyl-Substituents. Int. J. Mass Spectrom. 2012, 314, 13−17. (47) Raczyńska, E. D.; Decouzon, M.; Gal, J.-F.; Maria, P.-C.; Gelbard, G.; Vielfaure-Joly, F. Gas-Phase Structural (Internal) Effects in Strong Organic Nitrogen Bases. J. Phys. Org. Chem. 2001, 14, 25− 34. (48) Duchácǩ ová, L.; Kadlčiková, A.; Kotora, M.; Roithová, J. Oxygen Superbases as Polar Binding Pockets in Nonpolar Solvents. J. Am. Chem. Soc. 2010, 132, 12660−12667. (49) Tian, Z.; Fattahi, A.; Lis, L.; Kass, S. R. Neutral Intramolecular Hydrogen-Bonded Bases. Croat. Chem. Acta 2009, 82, 41−45. (50) Kolomeitsev, A. A.; Koppel, I. A.; Rodima, T.; Barten, J.; Lork, E.; Röschenthaler, G.-V.; Kaljurand, I.; Kütt, A.; Koppel, I.; Mäemets, V.; Leito, I. Guanidinophosphazenes: Design, Synthesis, and Basicity in THF and in the Gas Phase. J. Am. Chem. Soc. 2005, 127, 17656− 17666. (51) Kovače vić, B.; Maksić, Z. B. High Basicity of Tris(tetramethylguanidinyl)-phosphine Imide in the Gas Phase and Acetonitrile  a DFT study. Tetrahedron Lett. 2006, 47, 2553−2555. (52) Leito, I.; Koppel, I. A.; Koppel, I.; Kaupmees, K.; Tshepelevitsh, S.; Saame, J. Basicity Limits of Neutral Organic Superbases. Angew. Chem., Int. Ed. 2015, 54, 9262−9265. (53) Peran, N.; Maksić, Z. B. Polycyclic Croissant-like Organic Compounds are Powerful Superbases in the Gas Phase and Acetonitrile  a DFT study. Chem. Commun. 2011, 47, 1327−1329. (54) Kovačević, B.; Maksić, Z. B. The Proton Affinity of the Superbase 1,8-Bis(tetramethylguanidino)naphthalene (TMGN) and Some Related Compounds: A Theoretical Study. Chem. - Eur. J. 2002, 8, 1694−1702. (55) Margetič, D.; Ishikawa, T.; Kumamoto, T. Exceptional Superbasicity of Bis(guanidine) Proton Sponges Imposed by the Bis(secododecahedrane) Molecular Scaffold: A Computational Study. Eur. J. Org. Chem. 2010, 2010, 6563−6572. (56) Raab, V.; Gauchenova, E.; Merkoulov, A.; Harms, K.; Sundermeyer, J.; Kovač e vić, B.; Maksić, Z. B. 1,8-Bis(hexamethyltriaminophosphazenyl)naphthalene, HMPN: A Superbasic Bisphosphazene “Proton Sponge. J. Am. Chem. Soc. 2005, 127, 15738−15743. (57) Singh, A.; Ganguly, B. Strategic Design of Small and Versatile Bicyclic Organic Superbases: A Density Functional Study. New J. Chem. 2008, 32, 210−213. (58) Singh, A.; Ganguly, B. DFT Studies on a New Class of Cage Functionalized Organic Superbases. New J. Chem. 2009, 33, 583−587. (59) Belding, L.; Dudding, T. Synthesis and Theoretical Investigation of a 1,8-Bis(bis(diisopropylamino)cyclopropeniminyl)naphthalene Proton Sponge Derivative. Chem. - Eur. J. 2014, 20, 1032−1037. (60) Coles, M. P.; Aragón-Sáez, P. J.; Oakley, S. H.; Hitchcock, P. B.; Davidson, M. G.; Maksić, Z. B.; Vianello, R.; Leito, I.; Kaljurand, I.;

Apperley, D. Superbasicity of a Bis-guanidino Compound with a Flexible Linker: A Theoretical and Experimental Study. J. Am. Chem. Soc. 2009, 131, 16858−16868. (61) Kögel, J. F.; Xie, X.; Baal, E.; Gesevicius, D.; Oelkers, B.; Kovačević, B.; Sundermeyer, J. Superbasic Alkyl-Substituted Bisphosphazene Proton Sponges: Synthesis, Structural Features, Thermodynamic and Kinetic Basicity, Nucleophilicity and Coordination Chemistry. Chem. - Eur. J. 2014, 20, 7670−7685. (62) Kögel, J. F.; Oelkers, B.; Kovačević, B.; Sundermeyer, J. A New Synthetic Pathway to the Second and Third Generation of Superbasic Bisphosphazene Proton Sponges: The Run for the Best Chelating Ligand for a Proton. J. Am. Chem. Soc. 2013, 135, 17768−17774. (63) Lo, R.; Singh, A.; Kesharwani, M. K.; Ganguly, B. Rational Design of a New Class of Polycyclic Organic Bases Bearing Two Superbasic Sites and Their Applications in the CO2 Capture and Activation Process. Chem. Commun. 2012, 48, 5865−5867. (64) Radić, N.; Maksic, Z. B. Basicity of Neutral Organic Superbases with Vinamidine Structure in Gas Phase and Acetonitrile: A Density Functional Theory Study. J. Phys. Org. Chem. 2012, 25, 1168−1176. (65) Alder, R. W. Design of C2-Chiral Diamines That Are Computationally Predicted To Be a Million-fold More Basic than the Original Proton Sponges. J. Am. Chem. Soc. 2005, 127, 7924−7931. (66) Galeta, J.; Potáče k, M. Substituted Diazatetracyclo[4.4.0.13,10.15,8]dodecanes as Stable Caged Proton Sponges. J. Org. Chem. 2012, 77, 1010−1017. (67) Bachrach, S. M. Extended Hydrogen Bond Network Enabled Superbases. Org. Lett. 2012, 14, 5598−5601. (68) Štrukil, V.; Lekšić, E.; Meštrović, E.; Eckert-Maksić, M. The Synthesis and Structural Characterization of Peralkylated Triguanide Superbases. Aust. J. Chem. 2014, 67, 1129−1133. (69) Barić, D.; Kovačević, B. Designing a Next Generation of Proton Sponges: Cyclopropeniminophosphazenes as the Strongest Pincer Ligands. Tetrahedron Lett. 2016, 57, 442−445. (70) Kovačević, B.; Maksić, Z. B. Basicity of Some Organic Superbases in Acetonitrile. Org. Lett. 2001, 3, 1523−1526. (71) Kovačević, B.; Glasovac, Z.; Maksić, Z. B. The Intramolecular Hydrogen Bond and Intrinsic Proton Affinity of Neutral Organic Molecules: N,N′,N″-Tris(3-aminopropyl)guanidine and Some Related Systems. J. Phys. Org. Chem. 2002, 15, 765−774. (72) Barić, D.; Dragičević, I.; Kovačević, B. Design of Superbasic Guanidines − the Role of Multiple Intramolecular Hydrogen Bonds. J. Org. Chem. 2013, 78, 4075−4082. (73) Eckert-Maksić, M.; Glasovac, Z.; Trošelj, P.; Kütt, A.; Rodima, T.; Koppel, I.; Koppel, I. A. Basicity of Guanidines with Heteroalkyl Side Chains in Acetonitrile. Eur. J. Org. Chem. 2008, 2008, 5176−5184. (74) Glasovac, Z.; Eckert-Maksić, M. Effect of Intramolecular Hydrogen Bonds on the Gas-Phase Basicity of Guanidines. Aust. J. Chem. 2014, 67, 1056−1062. (75) Vazdar, K.; Kunetskiy, R.; Saame, J.; Kaupmees, K.; Leito, I.; Jahn, U. Very Strong Organosuperbases Formed by Combining Imidazole and Guanidine Bases: Synthesis, Structure, and Basicity. Angew. Chem., Int. Ed. 2014, 53, 1435−1438. (76) Barić, D.; Dragičević, I.; Kovačević, B. Cyclopropenimine as a Hydrogen Bond Acceptor− Towards the Strongest Non-phosphorus Superbases. Tetrahedron 2014, 70, 8571−8576. (77) Despotović, I.; Kovačević, B.; Maksić, Z. B. Hyperstrong Neutral Organic Bases: Phosphazeno Azacalix[3](2,6)pyridines. Org. Lett. 2007, 9, 4709−4712. (78) Despotović, I.; Kovačević, B.; Maksić, Z. B. Pyridine and sTriazine as Building Blocks of Nonionic Organic Superbases  A Density Functional Theory B3LYP Study. New J. Chem. 2007, 31, 447−457. (79) Terada, M.; Goto, K.; Oishi, M.; Takeda, T.; Kwon, E.; Kondoh, A. Design and Synthesis of Helically Chiral Spirocyclic P3 Phosphazenes and Characterization of Their Onium Salts. Synlett 2013, 24, 2531−2534. (80) Köhn, U.; Schulz, M.; Schramm, A.; Günther, W.; Görls, H.; Schenk, S.; Anders, E. A New Class of Chiral Phosphazene Bases: L

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Synthesis and Characterization. Eur. J. Org. Chem. 2006, 2006, 4128− 4134. (81) Kögel, J. F.; Kneusels, N.-J.; Sundermeyer, J. Two C2-symmetric Chelating P2-Bisphosphazene Superbases Connected Via a Binaphthyl Backbone − Synthesis, Structural Features and Preparation of a Cationic Alkyl Aluminum Complex. Chem. Commun. 2014, 50, 4319− 4321. (82) Despotović, I.; Maksić, Z. B.; Vianello, R. Design of Brønsted Neutral Organic Bases and Superbases by Computational DFT Methods: Cyclic and Polycyclic Quinones and [3]Carbonylradialenes. Eur. J. Org. Chem. 2007, 2007, 3402−3413. (83) Windus, T. L.; Schmidt, M. W.; Gordon, M. S. Theoretical Investigation of Azaphosphatrane Bases. J. Am. Chem. Soc. 1994, 116, 11449−11455. (84) Koppel, I. A.; Schwesinger, R.; Breuer, T.; Burk, P.; Herodes, K.; Koppel, I.; Leito, I.; Mishima, M. Intrinsic Basicities of Phosphorus Imines and Ylides: A Theoretical Study. J. Phys. Chem. A 2001, 105, 9575−9586. (85) Magill, A. M.; Cavell, K. J.; Yates, B. F. Basicity of Nucleophilic Carbenes in Aqueous and Nonaqueous Solvents − Theoretical Predictions. J. Am. Chem. Soc. 2004, 126, 8717−8724. (86) Tonner, R.; Heydenrych, G.; Frenking, G. First and Second Proton Affinities of Carbene Bases. ChemPhysChem 2008, 9, 1474− 1481. (87) Aldeco-Perez, E.; Rosenthal, A. J.; Donnadieu, B.; Parameswaran, P.; Frenking, G.; Bertrand, G. Isolation of a C5Deprotonated Imidazolium, a Crystalline “Abnormal” N-Heterocyclic Carbene. Science 2009, 326, 556−559. (88) Hollóczki, O.; Nyulászi, L. Neutral Species from “Non-protic” N-heterocyclic Ionic Liquids. Org. Biomol. Chem. 2011, 9, 2634−2640. (89) Hollóczki, O.; Kelemen, Z.; Nyulászi, L. On the Organocatalytic Activity of N-Heterocyclic Carbenes: Role of Sulfur in Thiamine. J. Org. Chem. 2012, 77, 6014−6022. (90) Despotović, I.; Vianello, R. Engineering Exceptionally Strong Oxygen Superbases with 1,8-Diazanaphthalene Di-N-oxides. Chem. Commun. 2014, 50, 10941−10944. (91) Ç iftcioğlu, G. A.; Trindle, C. Computational Estimates of Thermochemistry and pKa Values of Cyclopropenyl Imine Superbases. Int. J. Quantum Chem. 2014, 114, 392−399. (92) Biswas, A. K.; Lo, R.; Ganguly, B. First Principles Studies toward the Design of Silylene Superbases: A Density Functional Theory Study. J. Phys. Chem. A 2013, 117, 3109−3117. (93) Bayat, M.; Salehzadeh, S. Theoretical Studies on the Proton Affinities of Four Different Series of Nano-size Diamines and Designing Strong Superbases Based on Fullerene(C60) Molecule. J. Mol. Struct.: THEOCHEM 2010, 957, 120−125. (94) Bachrach, S. M.; Wilbanks, C. C. Using the Pyridine and Quinuclidine Scaffolds for Superbases: A DFT Study. J. Org. Chem. 2010, 75, 2651−2660. (95) Alexandrova, A. V.; Mašek, T.; Polyakova, S. M.; Císařová, I.; Saame, J.; Leito, I.; Lyapkalo, I. M. Synthesis of Electron-Rich Sterically Hindered P1 Phosphazene Bases by the Staudinger Reaction. Eur. J. Org. Chem. 2013, 2013, 1811−1823. (96) Kunetskiy, R. A.; Polyakova, S. M.; Vavřík, J.; Císařová, I.; Saame, J.; Nerut, E. R.; Koppel, I.; Koppel, I. A.; Kütt, A.; Leito, I.; Lyapkalo, I. M. A New Class of Organosuperbases, N-Alkyl- and NAryl-1,3-dialkyl-4,5-dimethylimidazol-2-ylidene Amines: Synthesis, Structure, pKBH+ Measurements, and Properties. Chem. - Eur. J. 2012, 18, 3621−3630. (97) Wolf, J. F.; Staley, R. H.; Koppel, I.; McIver, R. T., Jr.; Beauchamp, J. L.; Taft, R. W. Gas Phase Basicities and Relative Proton Affinities of Compounds Between Water and Ammonia from Pulsed Ion Cyclotron Resonance Thermal Equilibriums Measurements. J. Am. Chem. Soc. 1977, 99, 5417−5429. (98) Yamdagni, R.; Kebarle, P. Gas-Phase Basicities and Proton Affinities of Compounds between Water and Ammonia and Substituted Benzenes from a Contiuous Ladder of Proton Transfer Equilibria Measurements. J. Am. Chem. Soc. 1976, 98, 1320−1324.

(99) Cooks, R. G.; Wong, P. S. H. Kinetic Method of Making Thermodynamic Determinations: Advances and Applications. Acc. Chem. Res. 1998, 31, 379−386. (100) Toomsalu, E.; Koppel, I. A.; Burk, P. Critical Test of Some Computational Chemistry Methods for Prediction of Gas-Phase Acidities and Basicities. J. Chem. Theory Comput. 2013, 9, 3947−3958. (101) Koppel, I. A.; Taft, R. W.; Anvia, F.; Zhu, S.-Z.; Hu, L.-Q.; Sung, K.-S.; DesMarteau, D. D.; Yagupolskii, L. M.; Yagupolskii, Y. L.; Ignat’ev, N. V.; et al. The Gas-Phase Acidities of Very Strong Neutral Brønsted Acids. J. Am. Chem. Soc. 1994, 116, 3047−3057. (102) Leito, I.; Raamat, E.; Kütt, A.; Saame, J.; Kipper, K.; Koppel, I. A.; Koppel, I.; Zhang, M.; Mishima, M.; Yagupolskii, L. M.; et al. Revision of the Gas-Phase Acidity Scale below 300 kcal mol−1. J. Phys. Chem. A 2009, 113, 8421−8424. (103) Zhang, M.; Sonoda, T.; Mishima, M.; Honda, T.; Leito, I.; Koppel, I. A.; Bonrath, W.; Netscher, T. Gas-Phase Acidity of Bis[(perfluoroalkyl)sulfonyl]imides. Effects of the Perfluoroalkyl Group on the Acidity. J. Phys. Org. Chem. 2014, 27, 676−679. (104) Kaljurand, I.; Rodima, T.; Pihl, A.; Mäemets, V.; Leito, I.; Koppel, I. A.; Mishima, M. Acid-Base Equilibria in Nonpolar Media. 4. Extension of the Self-Consistent Basicity Scale in THF Medium. GasPhase Basicities of Phosphazenes. J. Org. Chem. 2003, 68, 9988−9993. (105) Rõõm, E.-I.; Kütt, A.; Kaljurand, I.; Koppel, I.; Leito, I.; Koppel, I. A.; Mishima, M.; Goto, K.; Miyahara, Y. Brønsted Basicities of Diamines in the Gas Phase, Acetonitrile, and Tetrahydrofuran. Chem. - Eur. J. 2007, 13, 7631−7643. (106) Mishima, M.; Mustanir; Fujio, M.; Tsuno, Y. Gas-Phase Substituent Effects in Stabilized Benzylic Carbocations. Basicities of Benzaldehydes, Acetophenones, and Methyl Benzoates in the Gas Phase. Bull. Chem. Soc. Jpn. 1996, 69, 2009−2018. (107) Kekišev, O.; Kaljurand, I.; Toom, L.; Lenoir, D.; Burk, P.; Järv, J. Effect of strain on gas-phase basicity of (E)-1-methyl-2-(1-methyl-2adamantylidene)adamantane. J. Phys. Org. Chem. 2015, 28, 447−451. (108) Adamson, A.; Kaljurand, I.; Guillemin, J.-C.; Burk, P. Relative stability and proton transfer reactions of unsaturated isocyanides and cyanides. J. Phys. Org. Chem. Accepted for publication on March 6, 2016, DOI: 10.1002/poc.3562. (109) Berglund, M.; Wieser, M. Isotopic Compositions of the Elements 2009 (IUPAC Technical Report). Pure Appl. Chem. 2011, 83, 397−410. (110) Hot-filament electron ionization source and Bayard-Alpert type hot filament ion gauge are based on the same physical principle. Thermally emitted electrons are accelerated, and these ionize neutral molecules upon collisions. Under constant flow and kinetic energy of electrons number of formed ions is proportional to the gas density and the size of the molecules forming the gas. (111) There is a novel and tested QSAR method112 available for calculating ion gauge sensitivity corrections for molecules having low molecular masses; however, this approach could not be used in the present work because this QSAR method does not include phosphorus atom, which is present in most of the molecules of current work. (112) Tämm, K.; Mayeux, C.; Sikk, L.; Gal, J.-F.; Burk, P. Theoretical Modeling of Sensitivity Factors of Bayard-Alpert Ionization Gauges. Int. J. Mass Spectrom. 2013, 341−342, 52−58. On-line engine for performing Sr calculations can be used on http://www.ut.ee/qsar (accessed March 28, 2016). (113) Bartmess, J. E.; Georgiadis, R. M. Empirical Methods for Determination of Ionization Gauge Relative Sensitivities for Different Gases. Vacuum 1983, 33, 149−153. (114) Miller, K. J. Additivity Methods in Molecular Polarizability. J. Am. Chem. Soc. 1990, 112, 8533−8542. (115) In the case of THF the term pKα was previously116,117,104,50 used for indicating that it is an estimate of pKa obtained from ion-pair measurements using the Fuoss equation.117 The basicity THF values in this work were obtained the same way, but we use the term pKa in this work because there is no intrinsic difference from pKa values. (116) Rodima, T.; Kaljurand, I.; Pihl, A.; Mäemets, V.; Leito, I.; Koppel, I. A. Acid-Base Equilibria in Nonpolar Media. 2.1 SelfConsistent Basicity Scale in THF Solution Ranging from 2M

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Methoxypyridine to EtP1(pyrr) Phosphazene. J. Org. Chem. 2002, 67, 1873−1881. (117) Abdur-Rashid, K.; Fong, T. P.; Greaves, B.; Gusev, D. G.; Hinman, J. G.; Landau, S. E.; Lough, A. J.; Morris, R. H. An Acidity Scale for Phosphorus-Containing Compounds Including Metal Hydrides and Dihydrogen Complexes in THF: Toward the Unification of Acidity Scales. J. Am. Chem. Soc. 2000, 122, 9155−9171. (118) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian 09, Revision C.1; Gaussian Inc: Wallingford, CT, 2009. (119) Turbomole V 6.4 2012, A development of University of Karlsruhe and Forschungzentrum Karlsruhe GmbH 1989−2007, Turbomole GmbH, since 2007, see http://www.turbomole.com. (accessed March 28, 2016). (120) Ahlrichs, R.; Bär, M.; Häser, M.; Horn, H.; Kølmel, C. Electronic Structure Calculations on Workstation Computers: the Program System Turbomole. Chem. Phys. Lett. 1989, 162, 165−169. (121) Schwesinger, R.; Willaredt, J.; Schlemper, H.; Keller, M.; Schmitt, D.; Fritz, H. Novel, Very Strong, Uncharged Auxiliary Bases; Design and Synthesis of Monomeric and Polymer-Bound Triaminoiminophosphorane Bases of Broadly Varied Steric Demand. Chem. Ber. 1994, 127, 2435−2454. (122) Schwesinger, R.; Schlemper, H.; Hasenfratz, C.; Willaredt, J.; Dambacher, T.; Breuer, T.; Ottaway, C.; Fletschinger, M.; Boele, J.; Fritz, H.; et al. Extremely Strong, Uncharged Auxiliary Bases; Monomeric and Polymer-Supported Polyaminophosphazenes (P2− P5). Liebigs Ann. 1996, 1996, 1055−1081. (123) Raab, V.; Kipke, J.; Gschwind, R. M.; Sundermeyer, J. 1,8Bis(tetramethylguanidino)naphthalene (TMGN): A New, Superbasic and Kinetically Active “Proton Sponge. Chem. - Eur. J. 2002, 8, 1682− 1692. (124) Kisanga, P. B.; Verkade, J. G.; Schwesinger, R. pK a Measurements of P(RNCH2CH3)3N. J. Org. Chem. 2000, 65, 5431− 5432. (125) Kaupmees, K.; Trummal, A.; Leito, I. Basicities of Strong Bases in Water: A Computational Study. Croat. Chem. Acta 2014, 87, 385− 395. (126) Hansch, C.; Leo, A.; Taft, R. W. A Survey of Hammett Substituent Constants and Resonance and Field Parameters. Chem. Rev. 1991, 91, 165−195. (127) Kaljurand, I.; Lilleorg, R.; Murumaa, A.; Mishima, M.; Burk, P.; Koppel, I.; Koppel, I. A.; Leito, I. The Basicity of Substituted N,Ndimethylanilines in Solution and in the Gas Phase. J. Phys. Org. Chem. 2013, 26, 171−181. (128) Kaljurand, I.; Rodima, T.; Leito, I.; Koppel, I. A.; Schwesinger, R. Self-Consistent Spectrophotometric Basicity Scale in Acetonitrile Covering the Range between Pyridine and DBU. J. Org. Chem. 2000, 65, 6202−6208. (129) Haav, K.; Saame, J.; Kütt, A.; Leito, I. Basicity of Phosphanes and Diphosphanes in Acetonitrile. Eur. J. Org. Chem. 2012, 2012, 2167−2172. (130) Rana, U. A.; Forsyth, M.; MacFarlane, D. R.; Pringle, J. M. Toward Protic Ionic Liquid and Organic Plastic Crystal Electrolytes for Fuel Cells. Electrochim. Acta 2012, 84, 213−222. (131) Smiglak, M.; Pringle, J. M.; Lu, X.; Han, L.; Zhang, S.; Gao, H.; MacFarlane, D. R.; Rogers, R. D. Ionic Liquids for Energy, Materials, and Medicine. Chem. Commun. 2014, 50, 9228−9250. (132) Pringle, J. M. Recent Progress in the Development and Use of Organic Ionic Plastic Crystal Electrolytes. Phys. Chem. Chem. Phys. 2013, 15, 1339−1351. (133) Mayeux, C.; Burk, P.; Gal, J.-F.; Kaljurand, I.; Koppel, I.; Leito, I.; Sikk, L. Gas-Phase Lithium Cation Basicity: Revisiting the High Basicity Range by Experiment and Theory. J. Am. Soc. Mass Spectrom. 2014, 25, 1962−1973. (134) Maksić, Z. B.; Kovačević, B.; Vianello, R. Advances in Determining the Absolute Proton Affinities of Neutral Organic Molecules in the Gas Phase and Their Interpretation: A Theoretical Account. Chem. Rev. 2012, 112, 5240−5270. N

DOI: 10.1021/acs.jpca.6b01552 J. Phys. Chem. A XXXX, XXX, XXX−XXX