Search PubMed⌕ Search

Biomedical subjects

John C Tully

Publications and source records attributed to John C Tully.

5 recordsLinked to original sources
PubMed2006

Vibrational lifetimes of molecular adsorbates on metal surfaces.

We report density functional theory calculations of electron-hole pair induced vibrational lifetimes of diatomic molecules adsorbed on metal surfaces. For CO on Cu(100), Ni(100), Ni(111), Pt(100), and Pt(111), we find that the C-O internal stretch and the bending modes have lifetimes in the 1-6 ps range, and that the CO-surface stretch and the frustrated translational modes relax more slowly, with lifetimes >10 ps for all cases except CO on Ni(111). This strong mode selectivity confirms earlier calculations for CO on Cu(100) and demonstrates that the trends carry over to other metal substrates. In contrast, for NO adsorbed on Pt(111), whereas we still find that the bending mode has the shortest lifetime, about 1.3 ps, we predict the other three modes to have almost equal lifetimes of 8-10 ps. Similarly, for CN adsorbed on Pt(111), we calculate that the internal stretching and molecule-surface stretching modes have approximately equal lifetimes of about 15 ps. Our results are in reasonable agreement with experiment, where available. We discuss some of the underlying factors that may contribute to the observed mode selectivity with adsorbed CO and the altered selectivity with NO and CN.

Journal Article↗
PubMed2005

Ab initio Ehrenfest dynamics.

We present an ab initio direct Ehrenfest dynamics scheme using a three time-step integrator. The three different time steps are implemented with nuclear velocity Verlet, nuclear-position-coupled midpoint Fock integrator, and time-dependent Hartree-Fock with a modified midpoint and unitary transformation algorithm. The computational cost of the ab initio direct Ehrenfest dynamics presented here is found to be only a factor of 2-4 larger than that of Born-Oppenheimer (BO) dynamics. As an example, we compute the vibration of the NaCl molecule and the intramolecular torsional motion of H2C=NH2+ by Ehrenfest dynamics compared with BO dynamics. For the vibration of NaCl with an initial kinetic energy of 1.16 eV, Ehrenfest dynamics converges to BO dynamics with the same vibrational frequency. The intramolecular rotation of H2C=NH2+ produces significant electronic excitation in the Ehrenfest trajectory. The amount of nonadiabaticity, suggested by the amplitude of the coherent progression of the excited and ground electronic states, is observed to be directly related to the strength of the electron-nuclear coupling. Such nonadiabaticity is seen to have a significant effect on the dynamics compared with the adiabatic approximation.

Journal Article↗
PubMed2005

Mixed quantum-classical equilibrium.

We present an analysis of the equilibrium limits of the two most widely used approaches for simulating the dynamics of molecular systems that combine both quantum and classical degrees of freedom. For a two-level quantum system connected to an infinite number of classical particles, we derive a simple analytical expression for the equilibrium mean energy attained by the self-consistent-field (Ehrenfest) method and show that it deviates substantially from Boltzmann. By contrast, "fewest switches" surface hopping achieves Boltzmann quantum state populations. We verify these analytical results with simulations.

Journal Article↗
PubMed2004

A washboard with moment of inertia model of gas-surface scattering.

A washboard with moment of inertia (WBMI) model for gas atom scattering from a flexible surface is proposed and applied. This model is a direct extension of the washboard model [J. Chem. Phys. 92, 680 (1990)] proposed for gas atom scattering from relatively rigid, corrugated surfaces. In addition, a moment of inertia is incorporated in the original washboard model to describe the flexibility of softer, more highly corrugated surfaces such as polymer or liquid surfaces. The moment of inertia of the effective surface object introduces a dependence of the efficiency of energy transfer on the position and direction of impact, a feature that has been shown to be critical by molecular dynamics simulations. The WBMI model is solved numerically by Monte Carlo integration, which makes the implementation of multiple impacts between a colliding atom and the surface very efficient. The model is applied to Ne and Ar atoms scattering from an alkylthiolate self-assembled monolayer surface and reproduces the major results obtained by classical trajectory simulation of the same system, i.e., a bimodal translation energy distribution P(E(f)) with the low-energy component well-fit with a Boltzmann distribution, but with a temperature that may (Ar) or may not (Ne) be the same as the surface temperature. This indicates that the WBMI model, with well-motivated physical assumptions and simplified interaction, reveals many of the major aspects of the gas-surface collision dynamics, though it does not take into account the real-time dynamics explicitly.

Journal Article↗