The anionic form of p-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) has been extensively employed as a model of the chromophore of the green fluorescence protein. The bright S1 excited state HBDIβˆ’ has a measured lifetime of 1.4 ps in the gas-phase and is dominated by two non-radiative decay mechanisms: internal conversion and autodetachment into the neutral continuum. Here, time-resolved photoelectron spectroscopy has been used to determine the yields of these two channels from which the lifetime for autodetachment was found to be ∼30 ps.

The discovery of the green fluorescent protein (GFP) and related fluorescent proteins has revolutionised bio-imaging as a marker for in vivo studies of biological systems.1–4 In many fluorescent proteins, the absorption and fluorescence can be attributed to a well-defined chromophore encapsulated in the protein structure.2,5 For GFP, p-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) has commonly been used as the representative chromophore. Fluorescence originates from the anionic form, HBDIβˆ’, where the phenoxyl group has been deprotonated (see Fig. 1).1,6 The photo-physical properties of HBDIβˆ’ have been extensively studied with the aim of gaining insight into the photophysics of GFP.7–13 From this perspective, studies on the isolated HBDIβˆ’ are particularly valuable as these provide detailed information on the intrinsic photophysics to which theoretical studies can be benchmarked.11,14–28

FIG. 1.

Excitation scheme for HBDIβˆ’ indicating the two prominent decay routes: internal conversion (IC) and autodetachment (AD). Also shown is the structure of HBDIβˆ’.

FIG. 1.

Excitation scheme for HBDIβˆ’ indicating the two prominent decay routes: internal conversion (IC) and autodetachment (AD). Also shown is the structure of HBDIβˆ’.

Close modal

An important difference between HBDIβˆ’ and GFP is that for the latter, the protein environment stabilises the anionic chromophore and increases the detachment energy by several eV.20,29 For HBDIβˆ’ in isolation, a number of photoelectron (PE) spectroscopy17–19 and theoretical studies15,16,30 have shown that the S1 state lies very close in energy to the neutral ground state, D0, although formally, the S1 is bound in the Franck-Condon region. Despite this, PE spectra from our group19 and the Andersen group17 have shown that autodetachment from the S1 state is possible, enabled by the finite temperature (∼300 K) of HBDIβˆ’ in these experiments. Autodetachment is an important channel, not only for the interpretation of the action (absorption) spectra of HBDIβˆ’, but also for any dynamical studies of the S1 excited state. The excited state dynamics have recently been probed using time-resolved PE spectroscopy.14 An overall observed lifetime of 1.4 ps was measured based on the decay of the PE signal following detachment from the S1 state by a probe pulse at 1.55 eV. However, it was not possible to deconvolute the individual contributions arising from internal conversion and autodetachment. In the present work, through careful measurements of the time-resolved PE spectra, the lifetime of both decay pathways are determined using the excitation scheme shown in Fig. 1.

The experimental setup and methodology has previously been described in detail.31–36 HBDI was synthesized based broadly on published methods.37 Isolated deprotonated HBDIβˆ’ ions were produced by electrospray ionisation (∼1 mM in MeOH, adjusted to pH 11 by the addition of NaOH). Ions were then stored in a ring-electrode trap and injected at 50 Hz into a time-of-flight mass-spectrometer, where HBDIβˆ’ ions were mass selected and then irradiated at the centre of a velocity-map-imaging PE spectrometer.31,38 Femtosecond laser pulses at 2.57 eV (482 nm), 2.48 eV (500 nm), and 2.39 eV (519 nm) served as the pump and 1.55 eV (800 nm) femtosecond pulses probed the population from the S1 excited state. The polarisation of all pulses was set parallel to the detector. All laser pulses were derived from a commercial femtosecond Ti:sapphire system. The spectral width of visible pulses was on the order of 60 meV (12 nm) full-width at half-maximum (FWHM) and the cross-correlation between pump and probe pulses was on the order of 120 fs. The resolution of the PE spectrometer is ∼5% and has been calibrated using the spectrum of Iβˆ’.

Figure 2 shows pump-probe PE spectra of HBDIβˆ’ excited at 2.48 eV (500 nm) with pump-probe delays of t = 0 and t = βˆ’2 ps. Similar spectra have been obtained for excitation at 2.57 eV (482 nm) and 2.39 eV (519 nm). Signal levels at 2.39 eV were significantly lower than at higher pump energies. In order to reduce the impact of long-term drifts in the experimental conditions, PE images were collected at each delay for 700 laser shots before toggling between the two delays. Raw PE images were collected until a desired signal level was attained, amounting to ∼4 Γ— 104 laser shots at 2.57 and 2.48 eV and >2 Γ— 105 laser shots at 2.39 eV. Photoelectron images were deconvoluted using the polar-onion peeling algorithm.39Β 

FIG. 2.

Time-resolved photoelectron (PE) spectra taken with 2.48 eV pump and 1.55 eV probe pulses at delays of t = βˆ’2 ps (black) and t = 0 ps (red). (Inset) The first 300 meV of a PE spectrum taken at 2.48 eV at low VMI voltages.

FIG. 2.

Time-resolved photoelectron (PE) spectra taken with 2.48 eV pump and 1.55 eV probe pulses at delays of t = βˆ’2 ps (black) and t = 0 ps (red). (Inset) The first 300 meV of a PE spectrum taken at 2.48 eV at low VMI voltages.

Close modal

The PE spectra at t = βˆ’2 ps are virtually identical to those recorded with the pump pulse only: the 1.55 eV probe energy is insufficient to directly detach from the ground state and is not resonant with any transitions. At higher intensity, some PE signal can be observed at 1.55 eV, which arises from 2-photon excitation. The main feature of the PE spectra at 2.57, 2.48, and 2.39 eV is an intense peak at low electron kinetic energy (eKE). This is highlighted in the inset of Fig. 2, for which the VMI voltage was reduced by a factor of ∼4. This feature has previously been assigned to vibrational autodetachment from the S1 state.17,19 Although the adiabatic detachment energy (∼2.6 eV) lies above the S1 origin, the 0.3 eV of internal energy due to the finite temperature (300 K)33 leads to transitions from vibrationally excited states that can overcome the adiabatic energy.

PEs generated by vibrational autodetachment from the S1 state have low eKE (<0.1 eV). This feature is in fact structured (see Fig. 2 inset), with a maximum at about 30 meV. The eKE of this feature is roughly independent of photon energy.17 At eKE < 10 meV, an additional PE signal is observed, which we had previously interpreted as originating from thermal electron loss from the ground S0 state (thermionic emission).19 The t = βˆ’2 ps PE spectra also show a PE feature centred at eKE ∼ 2 eV, which arises from two-photon, resonance-enhanced absorption.19Β 

The PE spectrum at t = 0 shows a new feature, centred around eKE ∼ 1.2 eV, which can be assigned to detachment by the 1.55 eV probe from the S1 excited state. As shown by experiments performed in our laboratory, this feature evolves in time, reflecting dynamics occurring on the S1 state, and the integrated PE signal decays on a 1.4 ps time scale.14 The statistics, however, were insufficient to characterise any dynamics occurring in the low energy feature, which we have overcome here. At t = 0, there is a depletion in the low energy feature, which is concomitant with the increase of the pump-probe feature at eKE ∼ 1.2 eV. It is clear from Fig. 2 that the total integrated depletion of the low eKE signal is much smaller than the signal arising from the exited state.

The integrated pump-probe feature which peaks at eKE ∼ 1.2 eV is proportional to the population of excited states in HBDIβˆ’, [S1]. Its population decay is susceptible to all available decay processes. Internal conversion and autodetachment are expected to dominate as shown schematically in Fig. 1. The observed rate constant of the S1 state measured in time-resolved experiments can be expressed as a sum of the rate constants for all possible decay routes from the excited state, i: kobs = Ο„obsβˆ’1 = |$\sum\nolimits_i {k_i }$|βˆ‘iki = (1.4 ps)βˆ’1. The low eKE feature is primarily a measure of the S1 excited states that have undergone autodetachment. Therefore, the integrated depletion of the low eKE feature in Fig. 2, [AD], is a measure of the population fraction of [S1] that has undergone autodetachment. The ratio of the rate constant for autodetachment, kAD, to all other decay rates can be expressed as

Given that |$k_{{\rm obs}} = \sum\nolimits_i {k_i } = k_{{\rm AD}} + \sum\nolimits_{i \ne {\rm AD}} {k_i }$|k obs =βˆ‘iki=k AD +βˆ‘iβ‰  AD ki⁠, we can now write that

where we have only considered autodetachment and internal conversion as the decay pathways from the S1 state. Taking [S1] as the integrated PE signal between 0.15 < eKE < 1.50 eV and [AD] as the integrated PE signal between 0 < eKE < 0.1 eV,40 and using the previously determined14Β kobs = (1.4 ps)βˆ’1, the lifetimes for autodetachment and internal conversion can be calculated. At 2.57, 2.48, and 2.39 eV, the autodetachment lifetime is 36, 28, and 33 ps, respectively. These are more than an order of magnitude longer than the observed lifetime so that the true internal conversion lifetime is only slightly longer than 1.4 ps (1.46, 1.47, and 1.46 ps at 2.57, 2.48, and 2.39 eV, respectively).

There are relatively large errors associated with the Ο„AD and Ο„IC lifetimes determined here, which we estimate to be on the order of Β±20%. This is in part due to the increased noise generated during image reconstruction at low eKE.39 More importantly, however, there is likely to be some overlap between the depleted signal and the [S1] signal. Finally, there are additional contributions to the low eKE feature originating from thermionic emission, which we have neglected.

The relative importance of autodetachment versus internal conversion observed in the current experiments is consistent with the recent work of Bochenkova and Andersen.16 In this, particular vibrational modes have been identified that are strongly coupled to specific decay pathways. Specifically, a high-frequency mode was correlated with the promotion of vibrational autodetachment into the continuum, while lower-frequency modes associated with torsional motion were correlated with the promotion of internal conversion. The high-frequency modes are the more Franck-Condon active, but it is the interplay between these two types of modes that leads to a balance between the rate of vibrational autodetachment and internal conversion.16 Using a statistical model, Bochenkova and Andersen16 were also able to estimate relative rates kAD:kIC at 2.48 eV (500 nm) and found that kAD:kIC ∼ 0.1.16 This is in reasonable agreement with the ratio of kAD:kIC ∼ 0.05 ± 0.01 determined here at 2.48 eV. Given the uncertainties in our experimental measurements, the model of Bochenkova and Andersen16 seems to have captured much of the underlying physics occurring on the S1 excited state.

Their model also predicts that, at higher excitation energies, the mode leading to autodetachment is more readily excited and this is expected to result in a sharp increase in the rate of autodetachment. At 2.57 eV (482 nm), the ratio kAD:kIC was expected to increase to 0.2 and at 2.75 eV (450 nm), autodetachment is expected to account for about half of the parent ion depletion.16 Our results do not follow this trend, although the spectral window probed here is probably too narrow to conclusively state that the kAD:kIC is constant as a function of energy. At 2.39 eV, we are exciting well below the S1 ← S0 0-0 transition which is at 2.52 eV (492 nm). At this pump energy, excitation is enabled by the fact that our pulses have a broad spectral profile (FWHM ∼ 12 nm) and by the finite temperature of the ground state ions (∼300 K). If the spectral overlap of the laser with the action spectrum at 300 K is considered, then it is clear that components of the laser pulse that are significantly blue-shifted drive the transition, although the central energy is at 2.39 eV. Additionally, the initial temperature used in the statistical model by Bochenkova and Andersen16 was 300 K. Although we estimate the internal temperature of HBDIβˆ’ in the current experiment to be similar to this, it may be slightly higher and hence a direct comparison must be cautioned. Ideally, we had hoped to perform our experiments at 2.75 eV (450 nm), but this is not possible using our current set-up.

Action spectra from the Jockusch group show that the ratio of electrons formed to fragmentation is ∼1:1 at 2.57 eV (482 nm).22,25 The former may be assigned to autodetachment, while the latter must to be a consequence of internal conversion, thus suggesting that kAD:kIC ∼ 1 at 2.57 eV (482 nm). This is neither consistent with the ratio determined here nor with that predicted by Bochenkova and Andersen.16 In the PE spectra at eKE < 10 meV, a peak can be attributed to thermionic emission from the S0 ground state. In our analysis above, we have neglected to consider the depletion in the peak at <10 meV. Its depletion at t = 0 is negligible compared to (and convoluted by) the depletion in the peak at 0.01 < eKE < 0.10 eV. The small depletion observed indicates that thermionic emission has a much longer lifetime than the vibrational autodetachment and internal conversion, as may be expected.41–43 The discrepancy between our results and the action spectra from the Jockusch group can, in part, be accounted for by this thermionic emission. We note that we can currently only monitor thermionic electrons over a few 100 ns window (limited by our VMI arrangement). If thermionic emission is much slower than this, which it likely will be, then the true yield of thermionic electrons is not accurately represented in our experiments and this feature will be much larger than measured here.

In conclusion, we have experimentally determined the lifetimes of electron autodetachment and internal conversion from the S1 excited state of the HBDI gaseous anion at three excitation wavelengths. The lifetime for autodetachment is on the order of 30 ps following excitation near the maximum of the S1 ← S0 absorption (action) spectrum and is consistent with the calculations by Bochenkova and Andersen.16 We have also noted the importance of thermionic emission as an operative electron loss channel over long time scales.

This work was supported by the EPSRC (EP/D073472) and Durham University. J.R.R.V. thanks the ERC for a Starting Grant.

1.
H.
β€ˆ
Niwa
,
S.
β€ˆ
Inouye
,
T.
β€ˆ
Hirano
,
T.
β€ˆ
Matsuno
,
S.
β€ˆ
Kojima
,
M.
β€ˆ
Kubota
,
M.
β€ˆ
Ohashi
, and
F. I.
β€ˆ
Tsuji
,
Proc. Natl. Acad. Sci. U.S.A.
β€ˆ
93
,
13617
(
1996
).
2.
F.
β€ˆ
Yang
,
L. G.
β€ˆ
Moss
, and
G. N.
β€ˆ
Phillips
,
Nat. Biotechnol.
β€ˆ
14
,
1246
(
1996
).
3.
G. H.
β€ˆ
Patterson
,
S. M.
β€ˆ
Knobel
,
W. D.
β€ˆ
Sharif
,
S. R.
β€ˆ
Kain
, and
D. W.
β€ˆ
Piston
,
Biophys. J.
β€ˆ
73
,
2782
(
1997
).
4.
R. Y.
β€ˆ
Tsien
,
Annu. Rev. Biochem.
β€ˆ
67
,
509
(
1998
).
5.
C. W.
β€ˆ
Cody
,
D. C.
β€ˆ
Prasher
,
W. M.
β€ˆ
Westler
,
F. G.
β€ˆ
Prendergast
, and
W. W.
β€ˆ
Ward
,
Biochemistry
β€ˆ
32
,
1212
(
1993
).
6.
M.
β€ˆ
Chattoraj
,
B. A.
β€ˆ
King
,
G. U.
β€ˆ
Bublitz
, and
S. G.
β€ˆ
Boxer
,
Proc. Natl. Acad. Sci. U.S.A.
β€ˆ
93
,
8362
(
1996
).
7.
J. J.
β€ˆ
van Thor
,
Chem. Soc. Rev.
β€ˆ
38
,
2935
(
2009
).
8.
P. J.
β€ˆ
Tonge
and
S. R.
β€ˆ
Meech
,
J. Photochem. Photobiol., A
β€ˆ
205
,
1
(
2009
).
9.
S. R.
β€ˆ
Meech
,
Chem. Soc. Rev.
β€ˆ
38
,
2922
(
2009
).
10.
L. H.
β€ˆ
Andersen
and
A. V.
β€ˆ
Bochenkova
,
Eur. Phys. J. D
β€ˆ
51
,
5
(
2009
).
11.
A.
β€ˆ
Toniolo
,
S.
β€ˆ
Olsen
,
L.
β€ˆ
Manohar
, and
T. J.
β€ˆ
Martinez
,
Faraday Discuss.
β€ˆ
127
,
149
(
2004
).
12.
L. H.
β€ˆ
Andersen
,
H.
β€ˆ
Bluhme
,
S.
β€ˆ
Boye
,
T. J. D.
β€ˆ
Jorgensen
,
H.
β€ˆ
Krogh
,
I. B.
β€ˆ
Nielsen
,
S. B.
β€ˆ
Nielsen
, and
A.
β€ˆ
Svendsen
,
Phys. Chem. Chem. Phys.
β€ˆ
6
,
2617
(
2004
).
13.
N. M.
β€ˆ
Webber
,
K. L.
β€ˆ
Litvinenko
, and
S. R.
β€ˆ
Meech
,
J. Phys. Chem. B
β€ˆ
105
,
8036
(
2001
).
14.
C. R. S.
β€ˆ
Mooney
,
D. A.
β€ˆ
Horke
,
A. S.
β€ˆ
Chatterley
,
A.
β€ˆ
Simperler
,
H. H.
β€ˆ
Fielding
, and
J. R. R.
β€ˆ
Verlet
,
Chem. Sci.
β€ˆ
4
,
921
(
2013
).
15.
K. B.
β€ˆ
Bravaya
and
A. I.
β€ˆ
Krylov
, β€œ
On the photodetachment from the green fluorescent protein chromophore
,”
J. Phys. Chem. A
(published online).
16.
A. V.
β€ˆ
Bochenkova
and
L. H.
β€ˆ
Andersen
,
Faraday Discuss.
β€ˆ
163
,
297
(
2013
).
17.
Y.
β€ˆ
Toker
,
D. B.
β€ˆ
Rahbek
,
B.
β€ˆ
Klaerke
,
A. V.
β€ˆ
Bochenkova
, and
L. H.
β€ˆ
Andersen
,
Phys. Rev. Lett.
β€ˆ
109
,
128101
(
2012
).
18.
C. R. S.
β€ˆ
Mooney
,
M. E.
β€ˆ
Sanz
,
A. R.
β€ˆ
McKay
,
R. J.
β€ˆ
Fitzmaurice
,
A. E.
β€ˆ
Aliev
,
S.
β€ˆ
Caddick
, and
H. H.
β€ˆ
Fielding
,
J. Phys. Chem. A
β€ˆ
116
,
7943
(
2012
).
19.
D. A.
β€ˆ
Horke
and
J. R. R.
β€ˆ
Verlet
,
Phys. Chem. Chem. Phys.
β€ˆ
14
,
8511
(
2012
).
20.
K. B.
β€ˆ
Bravaya
,
B. L.
β€ˆ
Grigorenko
,
A. V.
β€ˆ
Nemukhin
, and
A. I.
β€ˆ
Krylov
,
Acc. Chem. Res.
β€ˆ
45
,
265
(
2012
).
21.
G.
β€ˆ
Jonasson
,
J.-M.
β€ˆ
Teuler
,
G.
β€ˆ
Vallverdu
,
F.
β€ˆ
Merola
,
J.
β€ˆ
Ridard
,
B.
β€ˆ
Levy
, and
I.
β€ˆ
Demachy
,
J. Chem. Theory Comput.
β€ˆ
7
,
1990
(
2011
).
22.
M. W.
β€ˆ
Forbes
,
A. M.
β€ˆ
Nagy
, and
R. A.
β€ˆ
Jockusch
,
Int. J. Mass. Spectrom.
β€ˆ
308
,
155
(
2011
).
23.
K.
β€ˆ
Chingin
,
R. M.
β€ˆ
Balabin
,
V.
β€ˆ
Frankevich
,
K.
β€ˆ
Barylyuk
,
R.
β€ˆ
Nieckarz
,
P.
β€ˆ
Sagulenko
, and
R.
β€ˆ
Zenobi
,
Int. J. Mass. Spectrom.
β€ˆ
306
,
241
(
2011
).
24.
I. V.
β€ˆ
Polyakov
,
B. L.
β€ˆ
Grigorenko
,
E. M.
β€ˆ
Epifanovsky
,
A. I.
β€ˆ
Krylov
, and
A. V.
β€ˆ
Nemukhin
,
J. Chem. Theory Comput.
β€ˆ
6
,
2377
(
2010
).
25.
M. W.
β€ˆ
Forbes
and
R. A.
β€ˆ
Jockusch
,
J. Am. Chem. Soc.
β€ˆ
131
,
17038
(
2009
).
26.
C.
β€ˆ
Filippi
,
M.
β€ˆ
Ziccheddu
, and
F.
β€ˆ
Buda
,
J. Chem. Theory Comput.
β€ˆ
5
,
2074
(
2009
).
27.
T. J.
β€ˆ
Martinez
,
Acc. Chem. Res.
β€ˆ
39
,
119
(
2006
).
28.
M. E.
β€ˆ
Martin
,
F.
β€ˆ
Negri
, and
M.
β€ˆ
Olivucci
,
J. Am. Chem. Soc.
β€ˆ
126
,
5452
(
2004
).
29.
K. B.
β€ˆ
Bravaya
,
M. G.
β€ˆ
Khrenova
,
B. L.
β€ˆ
Grigorenko
,
A. V.
β€ˆ
Nemukhin
, and
A. I.
β€ˆ
Krylov
,
J. Phys. Chem. B
β€ˆ
115
,
8296
(
2011
).
30.
I.
β€ˆ
Polyakov
,
E.
β€ˆ
Epifanovsky
,
B.
β€ˆ
Grigorenko
,
A. I.
β€ˆ
Krylov
, and
A.
β€ˆ
Nemukhin
,
J. Chem. Theory Comput.
β€ˆ
5
,
1907
(
2009
).
31.
D. A.
β€ˆ
Horke
,
G. M.
β€ˆ
Roberts
,
J.
β€ˆ
Lecointre
, and
J. R. R.
β€ˆ
Verlet
,
Rev. Sci. Instrum.
β€ˆ
83
,
063101
(
2012
).
32.
J.
β€ˆ
Lecointre
,
G. M.
β€ˆ
Roberts
,
D. A.
β€ˆ
Horke
, and
J. R. R.
β€ˆ
Verlet
,
J. Phys. Chem. A
β€ˆ
114
,
11216
(
2010
).
33.
D. A.
β€ˆ
Horke
,
A. S.
β€ˆ
Chatterley
, and
J. R. R.
β€ˆ
Verlet
,
J. Phys. Chem. Lett.
β€ˆ
3
,
834
(
2012
).
34.
A. S.
β€ˆ
Chatterley
,
D. A.
β€ˆ
Horke
, and
J. R. R.
β€ˆ
Verlet
,
Phys. Chem. Chem. Phys.
β€ˆ
14
,
16155
(
2012
).
35.
D. A.
β€ˆ
Horke
and
J. R. R.
β€ˆ
Verlet
,
Phys. Chem. Chem. Phys.
β€ˆ
13
,
19546
(
2011
).
36.
J. R. R.
β€ˆ
Verlet
,
Chem. Soc. Rev.
β€ˆ
37
,
505
(
2008
).
37.
S.
β€ˆ
Kojima
,
H.
β€ˆ
Ohkawa
,
T.
β€ˆ
Hirano
,
S.
β€ˆ
Maki
,
H.
β€ˆ
Niwa
,
M.
β€ˆ
Ohashi
,
S.
β€ˆ
Inouye
, and
F. I.
β€ˆ
Tsuji
,
Tetrahedron Lett.
β€ˆ
39
,
5239
(
1998
).
38.
A.
β€ˆ
Eppink
and
D. H.
β€ˆ
Parker
,
Rev. Sci. Instrum.
β€ˆ
68
,
3477
(
1997
).
39.
G. M.
β€ˆ
Roberts
,
J. L.
β€ˆ
Nixon
,
J.
β€ˆ
Lecointre
,
E.
β€ˆ
Wrede
, and
J. R. R.
β€ˆ
Verlet
,
Rev. Sci. Instrum.
β€ˆ
80
,
053104
(
2009
).
40.
The intergration ranges defining [S1] for excitation at 2.38 and 2.57 eV were 0.15 < eKE < 1.40 eV and 0.15 < eKE < 1.60 eV, respectively, reflecting that fact that the maximum of the pump-probe feature shifts to higher eKE. The integration range for that of the vibrational autodetachment stays roughly unchanged at the different photon energies used.
41.
E. E. B.
β€ˆ
Campbell
and
R. D.
β€ˆ
Levine
,
Annu. Rev. Phys. Chem.
β€ˆ
51
,
65
(
2000
).
42.
J. R. R.
β€ˆ
Verlet
,
A. E.
β€ˆ
Bragg
,
A.
β€ˆ
Kammrath
,
O.
β€ˆ
Cheshnovsky
, and
D. M.
β€ˆ
Neumark
,
J. Chem. Phys.
β€ˆ
121
,
10015
(
2004
).
43.
J. R. R.
β€ˆ
Verlet
,
A.
β€ˆ
Kammrath
,
G. B.
β€ˆ
Griffin
, and
D. M.
β€ˆ
Neumark
,
J. Chem. Phys.
β€ˆ
123
,
231102
(
2005
).