Quantum versus classical electron transfer energy as reaction coordinate for the aqueous Ru2+/Ru3+ redox reaction

Applying density functional theory (DFT)-based molecular dynamics simulation methods we investigate the effect of explicit treatment of electronic structure on the solvation free energy of aqueous Ru²⁺ and Ru³⁺.Our approach is based on the Marcus theory of redox half reactions, focussing on the vertical energy gap for reduction or oxidation of a single aqua ion. We compare the fluctuations of the quantum and classical energy gap along the same equilibrium ab initio molecular dynamics trajectory for each oxidation state. The classical gap is evaluated using a standard point charge model for the charge distribution of the solvent molecules (water). The quantum gap is computed from the full DFT electronic ground state energies of reduced and oxidized species, thereby accounting for the delocalization of the electron in the donor orbital and reorganization of the electron cloud after electron transfer (ET). The fluctuations of the quantum ET energy are well approximated by gaussian statistics giving rise to parabolic free energy profiles. The curvature is found to be independent of the oxidation state in agreement with the linear response assumption underlying Marcus theory. By contrast, the diabatic free energy curves evaluated using the classical gap as order parameter, while also quadratic, are asymmetric reflecting the difference in oxidation state. The response of these two order parameters is further analysed by a comparison of the spectral density of the fluctuations and the corresponding reorganization free energies.     

Relative importance of basicity in the gas phase and in solution for determining selectivity in electrospray ionization mass spectrometry

Electrospray ionization mass spectrometry is a critically important technique for the determination of small molecules, but its application for this purpose is complicated by its selectivity. For positive ion ESI-MS analysis of basic analytes, several investigators have pointed to the importance of analyte basicity as a source of selectivity. Currently, however, it is not known whether basicity in the gas phase or in solution is ultimately most important in determining responsiveness. The objective of these studies was to investigate the relative importance of basicity in solution and in the gas phase as factors that predict selectivity in positive ion ESI-MS analysis. ESI-MS response was compared for a diverse series of protonatable analytes in two different solvents, neat methanol and methanol with 0.5% acetic acid. A correlation was observed between analyte pK(b) and electrospray response. However, the response for the analytes with very high pK(b) values was significantly higher than would be expected based on concentration of the protonated form or the analyte in solution, and this higher response did not appear to result from gas-phase proton transfer reactions. Although all of the analytes investigated had higher gas-phase basicities than the solvent, their relative responses were not dictated by gas-phase basicity. Higher response was observed for all of the analytes studied in acidified methanol compared with neat methanol, and this higher response was most pronounced for weakly basic analytes. These findings support the use of analyte pK(b) for rational method development in ESI-MS analysis of small molecules.     

First-principles prediction of the pK(a)s of anti-inflammatory oxicams

The gas- and aqueous-phase acidities of a series of oxicams have been computed by combining M05-2X/6-311+G(3df,2p) gas-phase free energies with solvation free energies from the CPCM-UAKS, COSMO-RS, and SMD solvent models. To facilitate accurate gas-phase calculations, a benchmarking study was further carried out to assess the performance of various density functional theory methods against the high-level composite method G3MP2(+). Oxicams are typically diprotic acids, and several tautomers are possible in each protonation state. The direct thermodynamic cycle and the proton exchange scheme have been employed to compute the microscopic pK(a)s on both solution- and gas-phase equilibrium conformers, and these were combined to yield the macroscopic pK(a) values. Using the direct cycle of pK(a) calculation, the CPCM-UAKS model delivered reasonably accurate results with MAD ~ 1, whereas the SMD and COSMO-RS models' performance was less satisfactory with MAD ~ 3. Comparison with experiment also indicates that direct cycle calculations based on solution conformers generally deliver better accuracy. The proton exchange cycle affords further improvement for all solvent models through systematic error cancellation and therefore provides better reliability for the pK(a) prediction of compounds of these types. The latter approach has been applied to predict the pK(a)s of several recently synthesized oxicam derivatives.     

Determination of the overlapping pKa values of resorcinol using UV-visible spectroscopy and DFT methods

In this paper we determine the overlapping pK(a) values of resorcinol in water, applying a UV-Vis spectroscopic method that uses absorbance diagrams. On the other hand, in order to explain the pK(a) values obtained, we also investigate the molecular conformations and solute-solvent interactions of the resorcinate anions, using ab initio and density functional theory methods. Several ionization reactions and equilibria in protic solvents, which possess a high hydrogen-bond-donor capability, are proposed. The mentioned reactions and equilibria constituted the indispensable theoretical basis to calculate the acidity constants of resorcinol. Basis sets at the HF/6-31 + G(d) and B3LYP/6-31 + G(d) levels of theory were used for calculations. Tomasi's method was used to analyze the formation of intermolecular hydrogen bonds between the resorcinate anions and water molecules. In this way, it was determined that in alkaline aqueous solutions the monoanion and dianion of resorcinol are solvated with two and four molecules of water, respectively. The agreement between the experimentally determined pK(a) values and those reported in the literature demonstrates the applicability and accurateness of the spectroscopic method here used. On the other hand, the agreement between the experimental and theoretically calculated pK(a) values provides solid support for the acid-base reactions proposed in this work.     

The oxidation of tyrosine and tryptophan studied by a molecular dynamics normal hydrogen electrode

The thermochemical constants for the oxidation of tyrosine and tryptophan through proton coupled electron transfer in aqueous solution have been computed applying a recently developed density functional theory (DFT) based molecular dynamics method for reversible elimination of protons and electrons. This method enables us to estimate the solvation free energy of a proton (H(+)) in a periodic model system from the free energy for the deprotonation of an aqueous hydronium ion (H(3)O(+)). Using the computed solvation free energy of H(+) as reference, the deprotonation and oxidation free energies of an aqueous species can be converted to pK(a) and normal hydrogen electrode (NHE) potentials. This conversion requires certain thermochemical corrections which were first presented in a similar study of the oxidation of hydrobenzoquinone [J. Cheng, M. Sulpizi, and M. Sprik, J. Chem. Phys. 131, 154504 (2009)]. Taking a different view of the thermodynamic status of the hydronium ion, these thermochemical corrections are revised in the present work. The key difference with the previous scheme is that the hydronium is now treated as an intermediate in the transfer of the proton from solution to the gas-phase. The accuracy of the method is assessed by a detailed comparison of the computed pK(a), NHE potentials and dehydrogenation free energies to experiment. As a further application of the technique, we have analyzed the role of the solvent in the oxidation of tyrosine by the tryptophan radical. The free energy change computed for this hydrogen atom transfer reaction is very similar to the gas-phase value, in agreement with experiment. The molecular dynamics results however, show that the minimal solvent effect on the reaction free energy is accompanied by a significant reorganization of the solvent.     

Theoretical evaluation of pK(a) in phosphoranes: implications for phosphate ester hydrolysis

Knowledge of the pK(a) of phosphoranes is important for the interpretation of phosphate ester hydrolysis. Calculated pK(a)'s of the model phosphorane, ethylene phosphorane, are reported. The method of calculation is based on the use of dimethyl phosphate as a reference state for evaluating relative pK(a) values, and on the optimization of the oxygen and acidic hydrogen van der Waals radii to give reasonable pK(1)(a), pK(2)(a), and pK(3)(a) for phosphoric acid in solution. Density functional theory is employed to calculate the gas-phase protonation energies, and continuum dielectric methods are used to determine the solvation corrections. The calculated pK(1)(a) and p(2)(a) for the model phosphorane are 7.9 and 14.3, respectively. These values are within the range of proposed experimental values, 6.5-11.0 for pK(1)(a), and 11.3-15.0 for pK(2)(a). The mechanistic implications of the calculated pK(a)'s are discussed.     

Microsolvated transition state models for improved insight into chemical properties and reaction mechanisms

Over the years, several methods have been developed to effectively represent the chemical behavior of solutes in solvents. The environmental effects arising due to solvation can generally be achieved either through inclusion of discrete solvent molecules or by inscribing into a cavity in a homogeneous and continuum dielectric medium. In both these approaches of computational origin, the perturbations on the solute induced by the surrounding solvent are at the focus of the problem. While the rigor and method of inclusion of solvent effects vary, such solvation models have found widespread applications, as evident from modern chemical literature. A hybrid method, commonly referred to as cluster-continuum model (CCM), brings together the key advantages of discrete and continuum models. In this perspective, we intend to highlight the latent potential of CCM toward obtaining accurate estimates on a number of properties as well as reactions of contemporary significance. The objective has generally been achieved by choosing illustrative examples from the literature, besides expending efforts to bring out the complementary advantages of CCM as compared to continuum or discrete solvation models. The majority of examples emanate from the prevalent applications of CCM to organic reactions, although a handful of interesting organometallic reactions have also been discussed. In addition, increasingly accurate computations of properties like pK(a) and solvation of ions obtained using the CCM protocol are also presented.     

Solid phase extraction purification of carboxylic acid products from 96-well format solution phase synthesis with DOWEX 1x8-400 formate anion exchange resin

The anion exchange resin DOWEX 1x8-400 formate has been developed for the isolation or resin capture of carboxylic acids from solution phase reactions in a 96-well format using a batchwise solid phase extraction technique. Eleven different anion exchange resins (formate forms) were evaluated for their efficiency at scavenging aryl and aliphatic carboxylic acids from solution. The model carboxylic acids had pK(a)s ranging from 3.40 to 4.89. Exchange efficiency onto the resin was pK(a) dependent with the carboxylic acids but not with their diisopropylethylammonium salts. Exchange off of the resin also showed pK(a) dependence with the stronger acids requiring more concentrated solvent acid for exchange. DOWEX 1x8-400 formate was determined to have superior capacity and the fastest exchange rate. Solvents suitable for exchanging the acids onto the resin were CH2Cl2, methanol, and various solvent/water mixtures. Solvents suitable for exchanging the carboxylic acids off of the resin were TFA/solvent or HCO2H/solvent mixtures. The resin was found to swell best in CH2Cl2 and in polar protic solvents such as water, alcohols, and acids. Application of this technique to the crude product mixtures from an arrayed reductive amination and an arrayed Stille reaction provided product carboxylic acids in yields averaging 57% and purities averaging 89%.     

Theoretical analysis of mechanistic pathways for hydrogen evolution catalyzed by cobaloximes

The mechanistic pathways for hydrogen evolution catalyzed by cobalt complexes with supporting diglyoxime ligands are analyzed with computational methods. The cobaloximes studied are Co(dmgBF(2))(2) (dmg = dimethylglyoxime) and Co(dpgBF(2))(2) (dpg = diphenylglyoxime) in acetonitrile. The reduction potentials and pK(a) values are calculated with density functional theory in conjunction with isodesmic reactions, incorporating the possibility of axial solvent ligand loss during the reduction process. The solvent reorganization energies for electron transfer between the cobalt complex and a metal electrode and the inner-sphere reorganization energies accounting for intramolecular rearrangements and the possibility of ligand loss are also calculated. The relative reduction potentials agree quantitatively with the available experimental values. The pK(a)s and reorganization energies agree qualitatively with estimates based on experimental data. The calculations suggest that a peak measured at ca. -1.0 V vs SCE in cyclic voltammetry experiments for Co(dmgBF(2))(2) is more likely to correspond to the Co(II)H/Co(I)H reduction potential than the Co(III)H/Co(II)H reduction potential. The calculations also predict pK(a) values of Co-hydride complexes and reduction potentials for both cobaloximes that have not been determined experimentally. The results are consistent with a mechanism in which the Co(III) and Co(II) complexes have two axial solvent ligands and the Co(I) complex has a single axial ligand along the reaction pathway. Analysis of the free energy diagrams generated for six different monometallic and bimetallic hydrogen production pathways identified the most favorable pathways for Co(dmgBF(2))(2) and tosic acid. The thermodynamically favored monometallic pathway passes through a Co(III)H intermediate, and Co(II)H reacts with the acid to produce H(2). The thermodynamically favored bimetallic pathways also pass through the Co(III)H intermediate, but the pathways in which two Co(III)H or two Co(II)H complexes react to produce H(2) are not thermodynamically distinguishable with these methods. On the basis of the electrostatic work term associated with bringing the two cobalt complexes together in solution, the preferred bimetallic pathway involves the reaction of two Co(III)H complexes to produce H(2). This mechanistic insight is important for designing more effective catalysts for solar energy conversion.     

Kinetics and mechanism of the pyridinolysis of S-2,4-dinitrophenyl 4-substituted thiobenzoates

[reaction: see text] The reactions of S-2,4-dinitrophenyl 4-methyl (1), S-2,4-dinitrophenyl 4-H (2), S-2,4-dinitrophenyl 4-chloro (3), and S-2,4-dinitrophenyl 4-nitro (4) thiobenzoates with a structurally homogeneous series of pyridines are subjected to a kinetic investigation in 44 wt % ethanol-water, at 25.0 degrees C and an ionic strength of 0.2 M (KCl). The reactions are studied spectrophotometrically (420 nm) by monitoring the appearance of 2,4-dinitrobenzenethiolate anion. Pseudo-first-order rate coefficients (k(obsd)) are obtained for all the reactions, employing excess of amine. The plots of k(obsd) vs [free pyridine] at constant pH are linear with the slopes (k(N)) independent of pH. The Br?nsted-type plots (log k(N) vs pK(a) of the conjugate acid of the pyridines) are curved for all the reactions. The Br?nsted curves are in accordance with stepwise mechanisms, through a zwitterionic tetrahedral intermediate (T(+/-)), and a change in the rate-limiting step. An equation based on this hypothesis accounts well for the experimental points. The Br?nsted lines were calculated with the following parameters: Reactions of thiolbenzoate 1: beta(1) 0.33 (slope at high pK(a)), beta(2) 0.95 (slope at low pK(a)), and pK(a)(0) = 8.5 (pK(a) at the curvature center); thiolbenzoate 2: beta(1) 0.30, beta(2) 0.88, and pK(a)(0) = 8.9; thiolbenzoate 3: beta(1) 0.33, beta(2) 0.89, and pK(a)(0) = 9.5; thiolbenzoate 4: beta(1) 0.21, beta(2) 0.97, and pK(a)(0) = 9.9. The increase of the pK(a)(0) value with the increase of the electron-withdrawing effect of the acyl substituent is explained by the argument that the rate of pyridine expulsion from T(+/-) (k(-)(1)) is favored over that of 2,4-dinitrobenzenethiolate leaving (k(2)), i.e., k(-)(1)/k(2) increases, as the acyl group becomes more electron withdrawing. The pK(a)(0) values for the title reactions are smaller than those for the reactions of the corresponding 4-nitrophenyl 4-substituted thiolbenzoates with the same pyridine series. This is explained by the larger k(2) value for 2,4-dinitrobenzenethiolate leaving from T(+/-) compared with 4-nitrobenzenethiolate, which results in lower k(-)(1)/k(2) ratios for the dinitro derivatives. The pK(a)(0) value obtained for the pyridinolysis of thiolbenzoate 2 (pK(a)(0) = 8.9) is smaller than that found for the same aminolysis of 2,4-dinitrophenyl benzoate (pK(a)(0) = 9.5). This is attributed to the greater nucleofugality from T(+/-) of 2,4-dinitrobenzenethiolate (pK(a) of conjugate acid 3.4) relative to 2,4-dinitrophenoxide (pK(a) of conjugate acid 4.1). The title reactions are also compared with the aminolysis of similar esters to assess the effect of the amine nature and leaving and acyl groups on the kinetics and mechanism.     

The kinetics and mechanisms of aromatic nucleophilic substitution reactions in liquid ammonia

The rates of aromatic nucleophilic substitution reactions in liquid ammonia are much faster than those in protic solvents indicating that liquid ammonia behaves like a typical dipolar aprotic solvent in its solvent effects on organic reactions. Nitrofluorobenzenes (NFBs) readily undergo solvolysis in liquid ammonia and 2-nitrofluorobenzene is about 30 times more reactive than the 4-substituted isomer. Oxygen nucleophiles, such as alkoxide and phenoxide ions, readily displace fluorine of 4-NFB in liquid ammonia to give the corresponding substitution product with little or no competing solvolysis product. Using the pK(a) of the substituted phenols in liquid ammonia, the Br?nsted β(nuc) for the reaction of 4-NFB with para-substituted phenoxides is 0.91, indicative of the removal of most of the negative charge on the oxygen anion and complete bond formation in the transition state and therefore suggests that the decomposition of the Meisenheimer σ-intermediate is rate limiting. The aminolysis of 4-NFB occurs without general base catalysis by the amine and the second-order rate constants generate a Br?nsted β(nuc) of 0.36 using either the pK(a) of aminium ion in acetonitrile or in water, which is also interpreted in terms of rate limiting breakdown of the Meisenheimer σ-intermediate. Nitrobenzene and diazene are formed as unusual products from the reaction between sodium azide and 4-NFB, which may be due to the initially formed 4-nitroazidobenzene decomposing to give a nitrene intermediate, which may then give diazene or be trapped by ammonia to give the unstable hydrazine which then yields nitrobenzene.     

Amplified potentiometric determination of pK(00), pK(0), pK(l), and pK(2) of hydrogen sulfides with Ag(2)S ISE

The chemical capacitor theory has been applied to accurately determine dissociation constants of H(2)S with the Ag(2)S ion-selective electrode (ISE). The theory's principle is based on the measurement of the change in electrode charge density as a result of protonated or unprotonated sulfide adsorbed on the electrode surface. This charge density is related to the potential. Connection of each individual capacitor in series amplifies the potential according to the equation, E(total)=E(1)+E(2)+E(3)+cdots, three dots, centeredE(n). As the charges of individual capacitors are concentrated to one capacitor area, the charge density rises, and the potential increases. The pK(00), pK(0), pK(1), and pK(2) are reported as 1.8, 2.12, 7.05, and 12.0, respectively. The pK(00) and pK(0) are reported here for the first time. The pK(1) agrees well with the literature values; however, the pK(2) differs from those reported recently under extreme conditions. Reasons for disproving the unreasonably high pK(2)>17-19 values are given based on calculations. Mainly, when pK(2)>17-19, the experimental results do not fit the equilibrium equations, pH=(pK(1)+pK(2))/2, pK(1)=(pK(0)+pK(2))/2, and pH=pK(2)+log(HS(-))/(S(2-)).     

Kinetic investigation of the reactions of S-4-nitrophenyl 4-substituted thiobenzoates with secondary alicyclic amines in aqueous ethanol

The reactions of S-4-nitrophenyl 4-X-substituted thiobenzoates (X = H, Cl, and NO(2): 1, 2, and 3, respectively) with a series of secondary alicyclic amines (SAA) were subjected to a kinetic investigation in 44 wt % ethanol-water, at 25.0 degrees C and an ionic strength of 0.2 M (KCl). The reactions were followed spectrophotometrically by monitoring the release of 4-nitrobenzenethiolate anion at 420-425 nm. Under excess amine, pseudo-first-order rate constants (k(obsd)) are obtained for all reactions. The plots of k(obsd) vs [SAA] at constant pH are linear with the slope (k(N)) independent of pH. The statistically corrected Br?nsted-type plots (log k(N)/q vs pK(a) + log p/q) for the reactions of 1 and 2 are nonlinear with slopes at high pK(a), beta(1) = 0.27 and 0.10, respectively, and slopes at low pK(a), beta(2) = 0.86 and 0.84, respectively. The Br?nsted curvature is centered at pK(a) (pK(a)(0)) 10.0 and 10.4, respectively. The reactions of SAA with 3 exhibit a linear Br?nsted-type plot of slope 0.81. These results are consistent with a stepwise mechanism, through a zwitterionic tetrahedral intermediate (T(+/-)). For the reactions of 1 and 2, there is a change in rate-determining step with amine basicity, from T(+/-) breakdown to products at low pK(a), to T(+/-) formation at high pK(a). For the reactions of 3, breakdown to products of T(+/-) is rate limiting for all the SAA series (pK(a)(0) > 11). The increasing pK(a)(0) value as the substituent in the acyl group becomes more electron withdrawing is attributed to an increasing nucleofugality of SAA from T(+/-). The greater pK(a)(0) value for the reactions of SAA with 1, relative to that found in the pyridinolysis of 2,4-dinitrophenyl benzoate (pK(a)(0) = 9.5), is explained by the greater nucleofugality from T(+/-) of the former amines, compared to isobasic pyridines, and the greater leaving ability from T(+/-) of 2,4-dinitrophenoxide relative to 4-nitrobenzenethiolate.     

A theoretical study on structures, energetics, and spectra of Br-.nCO2 clusters: towards bridging the gap between micro-domain and macro-domain

Structures, energetics, and spectra of Br(-).nCO(2) (n = 1-8) clusters are studied based on ab initio electronic structure theory. The geometry of each size of clusters is evaluated by employing second-order Moller-Plesset (MP2) perturbation theory. It is observed that the solvent CO(2) molecules approach the bromide moiety from one side in an asymmetric fashion except for the Br(-).8CO(2) cluster. Simple electrostatic model for charge-quadrupole interactions is valid for the Br(-).nCO(2) clusters. Reduced variational space based energy decomposition method shows that the electrostatic interaction is the major component and polarization and charge transfer energies are the other significant components of the total interaction energy. Both adiabatic and vertical electron detachment energies and solvation energies are calculated at MP2 level of theory. We have observed an excellent agreement between theory and experiment for the vertical detachment and solvation energies. Calculated quantities based on the analytical expression which connects the finite domain to macroscopic one are found to be very good in agreement with the available experimental results. The present study reveals a 2.6 eV increase in the detachment energy of bromide anion due to the solvation effect of CO(2), which is relatively small compared to that of the corresponding 4.7 eV increase in detachment energy in water.     

Determination of the overlapping pK(a) values of chrysin using UV-vis spectroscopy and ab initio methods

The overlapping pK(a) values of 5,7-dihydroxyflavone (chrysin) in EtOH-water solutions were determined by means of a UV-vis spectroscopic method that uses absorbance diagrams, at constant ionic strength (0.050 M) and temperature (25.0+/-0.1 degrees C). It was observed that the pK(a) values increase when the polarity-polarizability and solvation abilities of the reaction medium decrease. In order to calculate the pK(a1) and pK(a2) of chrysin in pure water, various relationships between the determined pK(a) and properties of solvents (relative permittivity, alpha-parameter of Taft and parameter Acity), are proposed. Moreover, with the aim of explaining the first pK(a1) value obtained, the molecular conformations and solute-solvent interactions of the 7(O(-))chrysinate monoanion were also investigated, using ab initio methods. Several ionization reactions and equilibria in water, which possesses a high hydrogen-bond-donor ability, are proposed. These reactions and equilibria constituted the necessary theoretical basis to calculate the first acidity constant of chrysin. The HF/6-31G(d) and HF/6-31+G(d) methods were used for calculations. Tomasi's method was used to analyze the formation of intermolecular hydrogen bonds between the 7(O(-))chrysinate monoanion and water molecules. It was proposed that in alkaline aqueous solutions the monoanion of chrysin is solvated with one water molecule. The agreement between the experimental and theoretical pK(a1) values provides good support for the acid-base reactions proposed in this paper.     

Dynamical quantum-electrodynamics embedding: Combining time-dependent density functional theory and the near-field method

We develop an approach for dynamical (ω > 0) embedding of mixed quantum mechanical (QM)∕classical (or more precisely QM∕electrodynamics) systems with a quantum sub-region, described by time-dependent density functional theory (TDDFT), within a classical sub-region, modeled here by the recently proposed near-field (NF) method. Both sub-systems are propagated simultaneously and are coupled through a common Coulomb potential. As a first step we implement the method to study the plasmonic response of a metal film which is half jellium-like QM and half classical. The resulting response is in good agreement with both full-scale TDDFT and the purely classical NF method. The embedding method is able to describe the optical response of the whole system while capturing quantum mechanical effects, so it is a promising approach for studying electrodynamics in hybrid molecules-metals nanostructures.     

Kinetic study of the aminolysis and pyridinolysis of O-phenyl and O-ethyl O-(2,4-dinitrophenyl) thiocarbonates. A remarkable leaving group effect

The reactions of a series of secondary alicyclic (SA) amines with O-phenyl and O-ethyl O-(2,4-dinitrophenyl) thiocarbonates (1 and 2, respectively) and of a series of pyridines with the former substrate are subjected to a kinetic investigation in water, at 25.0 degrees C, ionic strength 0.2 M (KCl). Under amine excess over the substrate, all the reactions obey pseudo-first-order kinetics and are first-order in amine. The Br?nsted-type plots are biphasic, with slopes (at high pK(a)) of beta(1) = 0.20 for the reactions of SA amines with 1 and 2 and beta(1) = 0.10 for the pyridinolysis of 1 and with slopes (at low pK(a)) of beta(2) = 0.80 for the reactions of SA amines with 1 and 2 and beta(2) = 1.0 for the pyridinolysis of 1. The pK(a) values at the curvature center (pK(a)(0)) are 7.7, 7.0, and 7.0, respectively. These results are consistent with the existence of a zwitterionic tetrahedral intermediate (T++) and a change in the rate-determining step with the variation of amine basicity. The larger pK(a)(0) value for the pyridinolysis of 1 compared to that for 2 (pK(a)(0) = 6.8) and the larger pK(a)(0) value for the reactions of SA amines with 1 relative to 2 are explained by the greater inductive electron withdrawal of PhO compared to EtO. The larger pK(a)(0) values for the reactions of SA amines with 1 and 2, relative to their corresponding pyridinolysis, are attributed to the greater nucleofugalities of SA amines compared to isobasic pyridines. The smaller pK(a)(0) value for the reactions of SA amines with 2 than with O-ethyl S-(2,4-dinitrophenyl) dithiocarbonate (pK(a)(0) = 9.2) is explained by the greater nucleofugality from T(++) of 2,4-dinitrophenoxide (DNPO(-)) relative to the thio derivative. The stepwise reactions of SA amines with 1 and 2, in contrast to the concerted mechanisms for the reactions of the same amines with the corresponding carbonates, is attributed to stabilization of T(++) by the change of O(-) to S(-). The simple mechanism for the SA aminolysis of 2 (only one tetrahedral intermediate, T(++)) is in contrast to the more complex mechanism (two tetrahedral intermediates, T(++) and T(-), the latter formed by deprotonation of T(++) by the amine) for the same aminolysis of the analogous thionocarbonate with 4-nitrophenoxide (NPO(-)) as nucleofuge. To our knowledge, this is the first example of a remarkable change in the decomposition path of a tetrahedral intermediate T by replacement of NPO(-) with DNPO(-) as the leaving group of the substrate. This is explained by (i) the greater leaving ability from T(++) of DNPO(-) than NPO(-) and (ii) the similar rates of deprotonation of both T(++) (formed with DNPO and NPO).