Article
Inorganic Chemistry, Vol. 50, No. 3, 2011 857
in vacuo. The crude solid was recrystallized by dissolving the
solid in dichloromethane, filtering to remove any insoluble
material, then layering methanol (2-3 volumes) on the solution.
After 2-3 days the dark red crystals were removed by filtration,
washed with a small amount of cold methanol, then diethyl
ether, and dried in vacuo. Yield=0.52 g, 26%. To grow crystals
of [Ni(4-Spy)(triphos)]BPh4 suitable for X-ray crystallography
the recrystallization was performed by dissolving in the mini-
mum of thf, then layering with methanol.
energy surface. Here, the energy was calculated with the struc-
ture optimized alongside, with energy convergence criterion
below 10-9 hartree. In one case MD was performed at constant
ꢀ
temperature using Nose-Hoover-Chain (NHC) thermostat.
All complex structure calculations were based on the structure
of [Ni(2-Spy)(triphos)]þ as determined by X-ray crystallogra-
phy. Conformation and intramolecular H-transfer barrier cal-
culations were based on minimal pyridinethiolate deformations.
ADF (Version 2004.01) is the Amsterdam Density Functional
program21 based on the density functional theory, (DFT), and
implemented through the use of Slater type basis orbitals. We
have employed the TZP (core double-ζ, valence triple-ζ with
polarization) basis set, which is one of the most extended sets of
the method, because it can provide efficiently convergent results.
Technically, the method is based on Vosko-Wilk-Nusair
(VWN) functional and the Perdew-Zunger self-interaction
correction (SIC) which is implemented self-consistently through
the use of the Krieger-Li-Iafrate approximation to the opti-
mized effective potential.22 Difficulties in the convergence of the
whole optimized structures did not allow us to proceed beyond
the Local Density Approximation level of calculation.
X-ray Crystallography. Data were collected on a Bruker
APEX2 diffractometer with synchrotron radiation (at Daresbury
Laboratory SRS station 9.8) for the 2-Spy complex and with a
Bruker SMART 1K diffractometer with MoKR radiation for
the 4-Spy complex. Semiempirical absorption corrections were
applied, based on repeated and symmetry-equivalent data. The
structures were solved by direct methods and refined on all
unique F2 data, with anisotropic displacement parameters for
non-H atoms, and isotropic displacement parameters for riding
H atoms. In the case of the 4-Spy complex, solvent molecules
were identified and refined, without H atoms as these could not
be unambiguously positioned; the asymmetric unit contains two
methanol molecules and half of a thf molecule, which is dis-
ordered over an inversion center such that it is not possible to
distinguish C and O atoms (all were refined as C, with restraints
on geometry and displacement parameters). Programs were
standard Bruker control and integration software, SHELXTL,
and local programs.
CPMD (Version 3.11.1)23 is a program for Car-Parrinello
molecular dynamics method, at the level of DFT through the use
of BLYP24,25 functional that includes a gradient correction
contribution as proposed by Becke24 and the correlation energy
term of Lee et al.25 The core-valence interactions were de-
scribed by semilocal norm-conserving Martins-Troullier pseu-
dopotentials,26 and the plane-wave expansion of the electronic
states included waves up to 100 Ry.27 During the calculations,
first, the energy was minimized for an initial structure, then the
structure was optimized and MD simulation was followed with
attention that the electron temperature would not increase. Only
a simplified structure of the complex was amenable to calcula-
tions. The full structure did not converge.
Kinetic Studies. All kinetic studies were conducted using an
Applied Photophysics SX 18MV stopped-flow spectrophotom-
eter, constructed to handle oxygen-sensitive solutions. The
temperature was maintained at 25.0 ( 0.1 °C using a Grant
LT D6G recirculating thermostatted tank. All solutions were
prepared under an atmosphere of dinitrogen and were trans-
ferred by gastight, all-glass syringes into the stopped-flow
spectrophotometer. The solutions of the Ni complexes and those
comprising mixtures of lut and [lutH]BPh4 (or [lutD]BPh4)
solutions were prepared from freshly made stock solutions
and were used within 2 h.
Acknowledgment. One of the authors (A.D.K.) thanks the
Computer Center of the National and Kapodistrian University
of Athens for their support in the computations. We thank
EPSRC for funding of crystallographic equipment and of the
National Crystallography Service, and STFC for access to
synchrotron facilities at Daresbury Laboratory.
All kinetics were studied under pseudo first order conditions
in acetonitrile. The concentration of the prepared stock solu-
tions of the nickel complexes was always 0.2 mmol dm-3, so that
when studied by stopped-flow spectrophotometry the concen-
tration was 0.1 mmol dm-3. The concentrations of lut and lutHþ
Supporting Information Available: X-raycrystallographic files
for [Ni(2-Spy)(triphos)]BPh4 and [Ni(4-Spy)(triphos)]BPh4
ranged from 1 to 30 mmol dm-3 and from 2 to 12 mmol dm-3
,
3
respectively. The ratio [lutHþ]/[lut] ranged from 0.25 to 5. The
reaction of the [Ni(S-2py)(triphos)]þ with lutHþ in the presence
of lut was studied at λ=385 nm. At this wavelength there is an
increase in absorbance, as shown in Figure 3. Presumably, this
increase in absorbance at this wavelength is a consequence of the
reaction involving the chelate ring-opening of the 2-thiopyridine
ligand with a necessary decrease in coordination number of the
Ni. Interestingly, the reaction of [Ni(SC6H4R-4)(triphos)]þ with
[lutH]þ, at λ = 350 nm, is associated with an absorbance
decrease, but this reaction only involves protonation of the
sulfur and no gross structural changes at Ni.5
2MeOH 0.5thf. Listings of tables containing the observed rate
3
constants for the reactions between [Ni(2-Spy)(triphos)]þ and
lutHþ in the presence of lut. Details of the theoretical calcula-
tions and video films showing calculated movement of proton
between the saddle point and nitrogen or sulfur. This material is
(17) Ahlswede, B.; Jug, K. J. Comput. Chem. 1999, 20, 563.
(18) Ahlswede, B.; Jug, K. J. Comput. Chem. 1999, 20, 572.
(19) Bredow, T.; Geudtner, G.; Jug, K. J. Comput. Chem. 2001, 22, 861.
(20) Steveson, M.; Bredow, T.; Gerson, A. Phys. Chem. Chem. Phys. 2002,
4, 358.
(21) ADF, Amsterdam Density Functional software; Scientific Computing and
Modelling N. V., Faculty of Sciences, Theoretical Chemistry, Vrije Universiteit: De
Boelelaan 1083, 1081 HVAmsterdam, The Netherlands (www.scm.com).
(22) (a) Patchkovskii, S.; Ziegler, T. J. Chem. Phys. 2001, 115, 26.
(b) Patchkovskii, S.; Ziegler, T. J. Chem. Phys. 2002, 116, 7806. (c) Patchkovs-
kii, S.; Ziegler, T. J. Phys. Chem. A 2002, 106, 1088.
(23) CPMD code, Version 3.11.1; Copyright IBM Corp. 1990-2006, Copy-
right MPI fur Festkorperforschung, Stuttgart, Germany, 1997-2001.
(24) Becke, A. D. Phys. Rev. A. 1988, 38, 3098.
(25) Lee, C. T.; Yang, W.; Parr, R. G. Phys. Rev. B. 1988, 37, 785.
(26) Troullier, N.; Martins, J. L. Phys. Rev. B. 1991, 43, 1993.
(27) Laasonen, K.; Pasquarello, A.; Car, R.; Lee, C.; Vanderbilt, D. Phys.
Rev. B. 1993, 47, 10142.
The dependences of kobs on the concentrations of lut and
lutHþ were determined graphically as described in the text. The
values of kobs presented are the average of at least three duplicate
experiments where the values of kobs all agree within 10%, and
this is reflected in the error bars presented in Figures 4, 7,
and 9.
Calculation Methods. MSINDO (Version 3.2.1)17-20 is a
method for quantum molecular dynamics at the level of the
Intermediate Neglect of Differential Overlap (INDO) semiem-
pirical self-consistent field (SCF) Molecular Orbital method
with MD calculations based on the lower Born-Oppenheimer