Unsaturated Amines, Phosphines, and Arsines
J. Am. Chem. Soc., Vol. 121, No. 19, 1999 4655
Table 1. Gas-Phase Acidities (kcal/mol) of Phosphines and
Vinylamine from Proton-Transfer Equilibrium Constant
Determinations and Bracketing Experiments
optimization was carried out at the B3LYP/6-31G* level. For
the anions, both MP2 and B3LYP geometry optimizations were
performed using a 6-31+G(d,p) basis set expansion. As
expected, for neutral and protonated species the inclusion of
diffuse components in the basis set leads to negligible changes
in the optimized geometries. The corresponding harmonic
vibrational frequencies were calculated at the same level of
theory used for the geometry optimization and the corresponding
ZPE corrections were scaled by the empirical factor 0.98.24 The
final energies were obtained using a 6-311+G(3df,2p) basis set
which was found to be adequate to reproduce the basicity and
∆
G°acid
∆∆G°acid
∆G°acid
a
b
(AH)c
AH
RefH
(RefH)
(338 K)
ethylphosphine
MeCO2Me
MeCN
i-PrCN
365.1
365.2
367.9
+0.92 ( 0.23
-0.28 ( 0.14
-1.39 ( 0.03
n-C5H11OH 368.2
CF3CH2OH 354.1
CF3CH2OH 354.1
HCONHMe 354.0
-2.17 ( 0.21 365.9 ( 0.8
-0.11 ( 0.06 354.0
-0.98 ( 0.04
HCONHMe
vinylphosphine
d
-1.22 ( 0.09 353.0 ( 0.3
+0.19 ( 0.07
ethynylphosphine n-PrSH
347.9
348.9
350.6
22
5
acidity of basis containing first- and third-row atoms. All the
ab initio and DFT calculations have been performed with the
EtSH
MeSH
-0.30 ( 0.04
-1.45 ( 0.08
2
5
CF3CH2OH 354.1 ,-3
348.6 ( 0.3
Gaussian-94 series of programs.
vinylamine
c-PrCN
t-BuOH
EtOH
367.8
367.7
371.7
+1.8 ( 0.2
g+2
g-2
Differences in reactivity usually reflect important dissimilari-
ties in charge distributions. We have used two different partition
techniques, namely the atoms in molecules (AIM) theory of
369.6
a
Absolute gas-phase acidities (Gibbs energies at 298.15 K for the
26
Bader and the natural bond orbital (NBO) analysis of Weinhold
-
+
b
reaction RefH f Ref + H ) from ref 6. Gibbs energies for the
27
-
-
et al. to analyze the charge distributions of amines, phosphines,
and arsines. The AIM theory is based in a topological analysis
reaction AH + Ref f A + RefH; quoted uncertainties correspond
to the standard deviation for three to four measurements. Upper and
c
2
lower limits were obtained from bracketing experiments. Absolute gas-
of the electron charge density F(r) and its Laplacian ∇ F(r).
phase acidities; no temperature correction applied, see text; quoted
uncertainties (standard deviations) correspond to the overlap quality,
and indicate the consistency with the existing absolute gas-phase acidity
scale. The uncertainty on this scale ((2 kcal/mol) should be added to
More specifically, we have located the so-called bond critical
points, i.e., points where F(r) is minimum along the bond path
and maximum in the other two directions. In general the values
d
2
evaluate the overall uncertainty on each value. Reevaluation, this work.
of F and ∇ F at these points provide useful information on the
bonding characteristics. In most cases negative values of the
Laplacian are associated with covalent linkages, while positive
values are usually associated with closed-shell interactions as
those found in ionic bonds, van der Waals complexes, and
hydrogen bonds. The NBO formalism permits a description of
the different bonds of the system in terms of the natural hybrid
orbitals centered on each atom and provides also useful
information on the charge distribution of the system.
G2-theory.18 G2 is a composite theory based on the 6-311G-
(d,p) basis set and several basis extensions, where electron
correlation effects are treated at the MP4 and QCISD(T) levels
of theory. The final energies are effectively at the QCISD(T)/
6
-311+G(3df,2p) level, assuming that basis set effects on the
correlation energies are additive. A small empirical correction
HLC) to accommodate remaining deficiencies is finally added
(
as well as the corresponding zero point energy (ZPE) correction,
estimated at the HF/6-31G* level. The reader is addressed to
ref 18 for a complete description of this method. Also, an
assessment of the G2 theory for the computation of enthalpies
of formation has recently been published.19
Since as mentioned above, electron correlation effects are
crucial for the systems under study, both population analyses
were performed at the MP2 level. For arsines the MP2-optimized
geometries reported in ref 5 were used. AIM calculations have
28
been carried out by using the AIMPAC series of programs.
For the particular case of phosphines we have considered it
of interest to investigate the performance of the B3LYP density
functional theory method,20 because it has been shown to yield
not only reliable geometries and vibrational frequencies for a
great variety of systems,21 but also to provide basicities and
acidities in fairly good agreement both with experimental
values,22 when available, and with high-level ab initio calcula-
tions, provided a flexible enough basis set is used. The B3LYP
approach combines the Becke’s three parameter non local hybrid
exchange potential with the non local correlation functional of
Lee, Yang, and Parr.23 In these DFT calculations, the geometry
Results and Discussion
Structures. Protonation and Deprotonation Effects. The
optimized geometries of ethyl-, vinyl-, and ethynylphosphine
and their protonated and deprotonated species are given as
Supporting Information. For the sake of comparison, the C-C
and C-X (X ) N, P, As) MP2(full)/6-31G* optimized bond
distances for the corresponding amines and arsines have been
summarized in Table 2.
It can be observed that the C-X bond length decreases in
the order ethyl > vinyl > ethynyl for the three series of
compounds. Consistently the charge density at the corresponding
bond critical points increases in the reverse order (See Table
(
18) Curtiss, L. A.; Raghavachari, K.; Trucks, G. W.; Pople, J. A. J.
Chem. Phys. 1991, 94, 7221.
19) Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Pople, J. A. J.
Chem. Phys. 1997, 106, 1063.
20) (a) Becke, A. D. J. Chem. Phys. 1993, 98, 5648. (b) Becke, A. D.
J. Chem. Phys. 1992, 96, 2155.
21) Llamas-Saiz, A. L.; Foces-Foces, C.; M o´ , O.; Y a´ n˜ ez, M.; Elguero,
(
(
(24) Bauschlicher Jr., C. W. Chem. Phys. Lett. 1995, 246, 40.
(25) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Gill, P: M.W.;
Johnson, B. J.; Robb, M. A.; Cheeseman, J. R.; Keith, T. A.; Peterson, G.
A.; Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; Zakrzewski,
V. G.; Ortiz, J. V.; Foresman, J. B.; Cioslowski, J.; Stefanow, B. B.;
Nanayaklara, A.; Challacombe, M.; Peng, C. Y.; Ayala, P. Y.; Chen, W.;
Wong, M. W.; Andres, J. L.; Replogle, E. S.; Gomperts, R.; Martin, R. L.;
Fox, D. J.; Binkley, J. S.; Defrees, D. J.; Baker, J.; Stewart, J. P.; Head-
Gordon, M.; Gonzalez, C.; Pople, J. A.; Gaussian 94; Gaussian, Inc.:
Pittsburgh, PA, 1995.
(
J. J. Comput. Chem. 1995, 16, 263. Flori a´ n, J.; Johnson, B. G. J. Phys.
Chem. 1994, 98, 3681. Bauschlicher, C. W. Chem. Phys. Lett. 1995, 246,
4
1
0; Martell, J. M.; Goddard, J. D.; Eriksson, L. A. J. Phys. Chem. A 1997,
01, 1927. Cui, Q.; Musaev, D. G.; Svensson, M.; Seiber, S.; Morokuma,
K. J. Am. Chem. Soc. 1995, 117, 12366. Ricca, A.; Bauschlicher, C. W. J.
Phys. Chem. 1995, 99, 5922. Ziegler, T. Chem. ReV. 1991, 91, 651.
(22) (a) Smith, B. J.; Radom, L. J. Am. Chem. Soc. 1993, 115, 4885. (b)
Smith, B. J.; Radom, L. Chem. Phys. Lett. 1994, 231, 345. (c) Gonz a´ lez,
A. I.; M o´ , O.; Y a´ n˜ ez, M.; Le o´ n, E.; Tortajada, J.; Morizur, J.-P.; Leito, I.;
Maria, P.-C.; Gal, J.-F. J. Phys. Chem. 1996, 100, 10490. (d) Amekraz, B.;
Tortajada, J.; Morizur, J.-P.; Gonz a´ lez, A. I.; M o´ , O.; Y a´ n˜ ez, M.; Leito, I.;
Maria, P.-C.; Gal, J.-F. New J. Chem. 1996, 20, 1011.
(26) Bader, R. F. W., Atoms in Molecules. A Quantum Theory. Oxford
University Press: Oxford, 1990.
(27) Weinhold, F.; Carpenter, J. E.; The Structure of Small Molecules
and Ions; Plenum: New York, 1988.
(28) The AIMPAC programs package has been provided by J. Cheeseman
and R. F. W. Bader.
(23) Lee, C.; Yang, W.; Parr, R. G. Phys ReV. 1988, B37, 785.