Conspectus
The reactivity
and dynamics of molecular systems can be explored
computationally by classical trajectory calculations. The traditional
approach involves fitting a functional form of a potential energy
surface (PES) to the energies from a large number of electronic structure
calculations and then integrating numerous trajectories on this fitted
PES to model the molecular dynamics. The ever-decreasing cost of computing
and continuing advances in computational chemistry software have made
it possible to use electronic structure calculations directly in molecular
dynamics simulations without first having to construct a fitted PES.
In this “on-the-fly” approach, every time the energy
and its derivatives are needed for the integration of the equations
of motion, they are obtained directly from quantum chemical calculations.
This approach started to become practical in the mid-1990s as a result
of increased availability of inexpensive computer resources and improved
computational chemistry software. The application of direct dynamics
calculations has grown rapidly over the last 25 years and would require
a lengthy review article. The present Account is limited to some of
our contributions to methods development and various applications.
To improve the efficiency of direct dynamics calculations, we developed
a Hessian-based predictor-corrector algorithm for integrating classical
trajectories. Hessian updating made this even more efficient. This
approach was also used to improve algorithms for following the steepest
descent reaction paths. For larger molecular systems, we developed
an extended Lagrangian approach in which the electronic structure
is propagated along with the molecular structure. Strong field chemistry
is a rapidly growing area, and to improve the accuracy of molecular
dynamics in intense laser fields, we included the time-varying electric
field in a novel predictor-corrector trajectory integration algorithm.
Since intense laser fields can excite and ionize molecules, we extended
our studies to include electron dynamics. Specifically, we developed
code for time-dependent configuration interaction electron dynamics
to simulate strong field ionization by intense laser pulses. Our initial
application of ab initio direct dynamics in 1994 was to CH2O → H2 + CO; the calculated vibrational distributions
in the products were in very good agreement with experiment. In the
intervening years, we have used direct dynamics to explore energy
partitioning in various dissociation reactions, unimolecular dissociations
yielding three fragments, reactions with branching after the transition
state, nonstatistical dynamics of chemically activated molecules,
dynamics of molecular fragmentation by intense infrared laser pulses,
selective activation of specific dissociation channels by aligned
intense infrared laser fields, angular dependence of strong field
ionization, and simulation of sequential double ionization.