PubMed Health⌕ Search

Biomedical subjects

M A Robb

Publications and source records attributed to M A Robb.

4 recordsLinked to original sources

A novel algorithm for non-adiabatic direct dynamics using variational Gaussian wavepackets.

In a recent paper (G. Worth, P. Hunt and M. Robb, J. Phys. Chem. A, 2003, 107, 621), we used surface hopping direct dynamics calculations to study the molecular dynamics of the butatriene radical cation in the X/A manifold, which is coupled by a conical intersection. Here, we present the first direct dynamics calculations using a novel algorithm, again using this ideal test system. The algorithm, which is based on the powerful multi-configuration time-dependent Hartree (MCTDH) wavepacket propagation method, uses a variational basis of coupled frozen Gaussian functions that optimally represent the evolving nuclear wavepacket at all times. Each Gaussian function follows a "quantum trajectory", along which the potential surface is evaluated by quantum chemistry calculations. As far fewer Gaussian functions are needed than classical trajectories in a semi-classical method, the number of quantum chemical calculations is drastically reduced. A crucial point in direct dynamics. To validate the method, initial calculations have been made using an analytic model Hamiltonian, where it is shown to reproduce the main features of the state population transfer with 8-16 basis functions per state. Coupled to the GAUSSIAN quantum chemistry program, the method is then shown to provide a feasible direct dynamics algorithm for the description of this non-adiabatic process.

Journal Article↗

The T(1) (3)(pi-pi*)/S(0) intersections and triplet lifetimes of cyclic alpha,beta-enones.

The ground and triplet excited states of cycloheptenone, cyclohexenone, and cyclopentenone have been studied using CASSCF calculations. For these three molecules, the difference in energy (DeltaE) between the twisted T(1) (3)(pi-pi*) minimum and T(1) (3)(pi-pi*)/S(0) intersection increases as the flexibility of the ring decreases. A strong positive correlation between DeltaE and the natural logarithm of the experimentally determined triplet lifetimes (ln tau) is found, suggesting that DeltaE predominantly determines the relative radiationless decay rates of T(1).

Journal Article↗

Potential energy surface crossings and the mechanistic spectrum for intramolecular electron transfer in organic radical cations.

The structure of the potential energy surface for the intramolecular electron transfer (IET) of four different model radical cations has been determined by using reaction path mapping and conical intersection optimization at the ab initio CASSCF level of theory. We show that, remarkably, the calculated paths reside in regions of the ground-state energy surface whose structure can be understood in terms of the position and properties of a surface crossing between the ground and the first excited state of the reactant. Thus, in the norbornadiene radical cation and in an analogue compound formed by two cyclopentene units linked by a norbornyl bridge, IET proceeds along direct-overlap and super-exchange concerted paths, respectively, that are located far from a sloped conical intersection point and in a region where the excited-state and ground-state surfaces are well separated. A second potential energy surface structure has been documented for 1,2-diamino ethane radical cation and features two parallel concerted (direct) and stepwise (chemical) paths. In this case a peaked conical intersection is located between the two paths. Finally, a third type of energy surface is documented for the bismethyleneadamantane radical cation and occurs when there is, effectively, a seam of intersection points (not a conical intersection) which separates the reactant and product regions. Since the reaction path cannot avoid the intersection, IET can only occur nonadiabatically. These IET paths indicate that quite different IET mechanisms may operate in radical cations, revealing an unexpectedly enriched and flexible mechanistic spectrum. We show that the origin of each path can be analyzed and understood in terms of the one-dimensional Marcus-Hush model.

Journal Article↗

Computational evidence in favor of a two-state, two-mode model of the retinal chromophore photoisomerization.

In this paper we use ab initio multiconfigurational second-order perturbation theory to establish the intrinsic photoisomerization path model of retinal chromophores. This is accomplished by computing the ground state (S(0)) and the first two singlet excited-state (S(1), S(2)) energies along the rigorously determined photoisomerization coordinate of the rhodopsin chromophore model 4-cis-gamma-methylnona-2,4,6,8-tetraeniminium cation and the bacteriorhodopsin chromophore model all-trans-hepta-2,4, 6-trieniminium cation in isolated conditions. The computed S(2) and S(1) energy profiles do not show any avoided crossing feature along the S(1) reaction path and maintain an energy gap >20 kcal small middle dotmol(-1). In addition, the analysis of the charge distribution shows that there is no qualitative change in the S(2) and S(1) electronic structure along the path. Thus, the S(1) state maintains a prevalent ionic (hole-pair) character whereas the S(2) state maintains a covalent (dot-dot) character. These results, together with the analysis of the S(1) reaction coordinate, support a two-state, two-mode model of the photoisomerization that constitutes a substantial revision of the previously proposed models.

Animals↗