Bipolar Reaction Path Hamiltonian Approach for Reactive Scattering Problems.

In this work we present a method for calculating the stationary state wave functions and reaction probabilities of a multidimensional reactive scattering system. Our approach builds upon the counter-propagating wave methodology (CPWM) developed by Poirier and co-workers for calculating one-dimensional stationary state wave functions. The method involves the formulation of a bipolar decomposition for multidimensional stationary scattering wave functions within the context of a reaction path Hamiltonian, so we refer to this work as the bipolar reaction path Hamiltonian (BRPH) approach. Benchmark calculations are presented for several 2D model scattering systems with linear reaction coordinates. We show that the BRPH approach is competitive with conventional calculations based on discrete variable representation (DVR) methods.

Supersymmetric approach to excited states.

We present here a supersymmetric (SUSY) approach for determining excitation energies within the context of a quantum Monte Carlo scheme. By using the fact that SUSY quantum mechanics gives rises to a series of isospectral Hamiltonians, we show that Monte Carlo ground-state calculations in the SUSY partners can be used to reconstruct accurately both the spectrum and states of an arbitrary Schrodinger equation. Since the ground state of each partner potential is nodeless, we avoid any "node" problem typically associated with the Monte Carlo technique. Although we provide an example of using this approach to determine the tunneling states in a double-well potential, the method is applicable to any 1D potential problem. We conclude by discussing the extension to higher dimensions.

Single-electron counting spectroscopy: simulation study of porphyrin in a molecular junction.

Electron counting of a single porphyrin molecule between two electrodes shows a crossover from sub- to super-Poissonian statistics as the bias voltage is scanned. This is attributed to the simultaneous activation of states with electron transfer rates spanning several orders of magnitude. Time-series analysis of consecutive single-electron transfer events reveals fast and slow transport channels, which are not resolved by the average current alone.

Optical coherence and theoretical study of the excitation dynamics of a highly symmetric cyclophane-linked oligophenylenevinylene dimer.

Optoelectronic properties of a polyphenylenevinylene-based oligomer and its paracylophane-linked dimer are studied using a variety of experimental and theoretical techniques. Despite the symmetrical structure and redshifted absorption of the dimer versus the monomer, an exciton picture is not the most appropriate. Electronic structure calculations establish changes in charge density upon optical excitation and show localized excitations that cannot be accounted for by a simple Frenkel exciton model. Visible frequency pump-probe anisotropy measurements suggest that the dimer should be considered as a three-level system with a fast, approximately 130 fs, internal conversion from the higher to lower energy excited electronic state. Signatures of nuclear relaxation processes are compared for electric field-resolved transient grating and two-dimensional photon echo spectra. These measurements reveal that nuclear relaxation occurs on similar time scales for the monomer and dimer. The connection between the spectral phase of four-wave mixing signals and the time dependent width of a nuclear wave packet is discussed. Semiempirical electronic structure and metropolis Monte Carlo calculations show that the dominant line broadening mechanisms for the monomer and dimer are associated with inter-ring torsional coordinates. Together, the theoretical calculations and electric field-resolved four-wave mixing experiments suggest that while the structure of dimer is more rigid than that of monomer, the difference in their rigidities is not sufficient to slow down excited state relaxation of dimer with respect to the monomer.

Simulation of single molecule inelastic electron tunneling signals in paraphenylene-vinylene oligomers and distyrylbenzene[2.2]paracyclophanes.

Inelastic resonances in the electron tunneling spectra of several conjugated molecules are simulated using the nonequilibrium Greens function formalism. The vibrational modes that strongly couple to the electronic current are different from the infrared and Raman active modes. Spatially resolved inelastic electron tunneling (IET) intensities are predicted. The simulated IET intensities for a large distyrylbenzene paracyclophane molecule are in qualitative agreement with recent experimental results.

All-forward semiclassical simulations of nonlinear response functions.

We propose a quantum trajectory algorithm for computing nonlinear response functions of condensed phase molecular systems based on a time-ordered expansion of the density matrix. The nth-order response function is expressed as a sum of 2(n) impulsive response pathways representing trajectories involving zero, one, and up to n interactions with short external pulses. These are evaluated using a forward propagation algorithm based upon a Liouville space extension of the Bohmian propagation method.

Quantum dissipation in unbounded systems.

In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie--Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wave function of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to nonclassical equations of motion. Until very recently these schemes have been limited to conservative systems. In this paper, we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme that incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy, we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.