# Publications

### 2021

A dual formalism for Lagrange multipliers is developed. The formalism is used to minimize an action function $S(q_2,q_1,T)$ without any dynamical input other than that $S$ is convex. All the key equations of analytical mechanics -- the Hamilton-Jacobi equation, the generating functions for canonical transformations, Hamilton's equations of motion and $S$ as the time integral of the Lagrangian -- emerge as simple consequences. It appears that to a large extent, analytical mechanics is simply a footnote to the most basic problem in the calculus of variations: that the shortest distance between two points is a straight line.

We present a new trajectory formulation of high harmonic generation that treats classically allowed and classically forbidden processes within a single dynamical framework. Complex trajectories orbit the nucleus, producing the stationary Coulomb ground state. When the field is turned on, these complex trajectories continue their motion in the field-dressed Coulomb potential and therefore tunnel ionization, unbound evolution and recollision are described within a single, seamless framework. The new formulation can bring mechanistic understanding to a broad range of strong field physics effects.

### 2019

The absorption spectrum of the vibronically allowed S1(1A2) ← S0(1A1) transition of formaldehyde is computed by combining multiplicative neural network (NN) potential surface fits, based on multireference electronic structure data, with the two-layer Gaussian-based multiconfiguration time-dependent Hartree (2L-GMCTDH) method. The NN potential surface fit avoids the local harmonic approximation for the evaluation of the potential energy matrix elements. Importantly, the NN surface can be constructed so as to be physically well-behaved outside the domain spanned by the ab initio data points. A comparison with experimental results shows spectroscopic accuracy of the converged surface and 2L-GMCTDH quantum dynamics.

We propose and demonstrate, numerically and experimentally, use of sparsity as prior information for extending the capabilities and performance of techniques and devices for laser pulse diagnostics. We apply the concept of sparsity in three different applications. First, we improve a photodiode-oscilloscope system’s resolution for measuring the intensity structure of laser pulses. Second, we demonstrate the intensity profile reconstruction of ultrashort laser pulses from intensity autocorrelation measurements. Finally, we use a sparse representation of pulses (amplitudes and phases) to retrieve measured pulses from incomplete spectrograms of cross-correlation frequency-resolved optical gating traces.

### 2018

*Advances in Chemical Physics, volume 163*. Whaley KB.(eds.). p. 273-323 (1Advances in Chemical Physics). Abstract

This chapter provides a simple pedagogical presentation of the discrete variable representation (DVR). It reviews the von Neumann (vN) basis of phase‐space Gaussians which include the Projected von Neumann Basis (PvN) and the Biorthogonal von Neumann Basis (PvB). The chapter includes a variety of interesting formal properties of nonorthogonal bases that are an extension of the DVR presentation and provides insight into the method. It presents an analysis of multidimensional considerations, including details of a highly efficient tensor formulation for performing pruned multidimensional DVR calculations for sparse but unstructured grids. The chapter contains illustrative applications. Pruned phase‐space methods have been successfully used for computing eigenenergies of (ro‐) vibrational systems. The chapter focuses on the applications in the context of solving the time‐dependent Schrodinger equation (TDSE). The efficiency of phase‐space versus coordinate‐space methods will certainly depend on the particular system studied and the strength of coupling between degrees of freedom.

Quantum computation places very stringent demands on gate fidelities, and experimental implementations require both the controls and the resultant dynamics to conform to hardware-specific constraints. Superconducting qubits present the additional requirement that pulses must have simple parameterizations, so they can be further calibrated in the experiment, to compensate for uncertainties in system parameters. Other quantum technologies, such as sensing, require extremely high fidelities. We present a novel, conceptually simple and easy-to-implement gradient-based optimal control technique named gradient optimization of analytic controls (GOAT), which satisfies all the above requirements, unlike previous approaches. To demonstrate GOAT's capabilities, with emphasis on flexibility and ease of subsequent calibration, we optimize fast coherence-limited pulses for two leading superconducting qubits architectures-flux-tunable transmons and fixed-frequency transmons with tunable couplers.

Stokes phenomenon refers to the fact that an asymptotic expansion of complex functions can differ in different regions of the complex plane, and that beyond the so-called Stokes lines the expansion has an unphysical divergence. An important special case is when the Stokes lines emanate from phase space caustics of a complex trajectory manifold. In this case, symmetry determines that to second order there is a double coverage of the space, one portion of which is unphysical. Building on the seminal but laconic findings of Adachi, we show that the deviation from second order can be used to rigorously determine the Stokes lines and therefore the region of the space that should be removed. The method has applications to wavepacket reconstruction from complex valued classical trajectories. With a rigorous method in hand for removing unphysical divergences, we demonstrate excellent wavepacket reconstruction for the Morse, Quartic, Coulomb, and Eckart systems.

Complex-valued classical trajectories in complex time encounter singular times at which the momentum diverges. A closed time contour around such a singular time may result in final values for q and p that differ from their initial values. In this work, we develop a calculus for determining the exponent and prefactor of the asymptotic time dependence of p from the singularities of the potential as the singularity time is approached. We identify this exponent with the number of singularity loops giving distinct solutions to Hamilton's equations of motion. The theory is illustrated for the Eckart, Coulomb, Morse, and quartic potentials. Collectively, these potentials illustrate a wide variety of situations: poles and essential singularities at finite and infinite coordinate values. We demonstrate quantitative agreement between analytical and numerical exponents and prefactors, as well as the connection between the exponent and the time circuit count. This work provides the theoretical underpinnings for the choice of time contours described in the studies of Doll et al. [J. Chem. Phys. 58(4), 1343-1351 (1973)] and Petersen and Kay [J. Chem. Phys. 141(5), 054114 (2014)]. It also has implications for wavepacket reconstruction from complex classical trajectories when multiple branches of trajectories are involved.

### 2017

Complex-valued semiclassical methods hold out the promise of treating classically allowed and classically forbidden processes on the same footing. In addition, they provide a natural way to describe optical excitation with complex fields within the trajectory framework. Despite their promise, these methods have until now been limited to short time propagation, due to the numerical difficulties introduced by the complexification. Using a new Final Value Representation of the Coherent State Propagator (FINCO), combined with an analysis of the complex classical phase space, we achieve accurate wavepacket propagation all the way to the revival time of a strongly anharmonic system. (C) 2017 Elsevier B. V. All rights reserved.

### 2016

We present an efficient implementation of dynamically pruned quantum dynamics, both in coordinate space and in phase space. We combine the ideas behind the biorthogonal von Neumann basis (PvB) with the orthogonalized momentum-symmetrized Gaussians (Weylets) to create a newbasis, projected Weylets, that takes the best from both methods. We benchmark pruned time-dependent dynamics using phase-space-localized PvB, projectedWeylets, and coordinate-space-localized DVR bases, with real-world examples in up to six dimensions. For the examples studied, coordinate-space localization is the most important factor for efficient pruning and the pruned dynamics is much faster than the unpruned, exact dynamics. Phase-space localization is useful for more demanding dynamics where many basis functions are required. There, projected Weylets offer a more compact representation than pruned DVR bases. Published by AIP Publishing.

We describe the mathematical underpinnings of the biorthogonal von Neumann method for quantum mechanical simulations (PvB). In particular, we present a detailed discussion of the important issue of nonorthogonal projection onto subspaces of biorthogonal bases, and how this differs from orthogonal projection. We present various representations of the Schrodinger equation in the reduced basis and discuss their relative merits. We conclude with illustrative examples and a discussion of the outlook and challenges ahead for the PvB representation.

### 2015

Among the major challenges in the chemical sciences is controlling chemical reactions and deciphering their mechanisms. Since much of chemistry occurs in excited electronic states, in the last three decades scientists have employed a wide variety of experimental techniques and theoretical methods to recover excited-state potential energy surfaces and the wavepackets that evolve on them. These methods have been partially successful but generally do not provide a complete reconstruction of either the excited state wavepacket or potential. We have recently proposed a methodology for reconstructing excited-state molecular wavepackets and the corresponding potential energy surface [Avisar and Tannor, Phys. Rev. Lett., 2011, 106, 170405]. In our approach, the wavepacket is represented as a superposition of the set of vibrational eigenfunctions of the molecular ground-state Hamiltonian. We assume that the multidimensional ground-state potential surface is known, and therefore these vibrational eigenfunctions are known as well. The time-dependent coefficients of the basis functions are obtained by experimental measurement of the resonant coherent anti-Stokes Raman scattering (CARS) signal. Our reconstruction strategy has several significant advantages: (1) the methodology requires no a priori knowledge of any excited-state potential. (2) It applies to dissociative as well as to bound excited-state potentials. (3) It is general for polyatomics. (4) The excited-state potential surface is reconstructed simultaneously with the wavepacket. Apart from making a general contribution to the field of excited-state spectroscopy, our method provides the information on the excited-state wavepacket and potential necessary to design laser pulse sequences to control photochemical reactions.

### 2014

We extend the periodic von Neumann basis to non-Cartesian coordinates. The bound states of two isomerizing triatomic molecules, LiCN/LiNC and HCN/HNC, are calculated using the vibrational Hamiltonian in Jacobi coordinates. The phase space localization of the basis functions leads to a flexible and accurate representation of the Hamiltonian. This results in significant savings compared to a basis localized just in coordinate space. The favorable scaling of the method with dimensionality makes it promising for applications to larger systems. (C) 2014 AIP Publishing LLC.

*Advances In Chemical Physics*. Vol. 156. p. 1-34 (1Advances in Chemical Physics). Abstract

One of the striking parallels in the development of the method in quantum mechanics and signal processing is that the method never became mainstream in either community. A key reason is undoubtedly the problems encountered in converging the method, problems reported independently in both fields. This chapter presents a theory that optimizes the phase space localization to obtain the most compact possible representation. It shows applications to shaped femtosecond pulses, to the solution of time‐independent and time‐dependent quantum mechanical problems, and to audio and image compression. The chapter then discusses the application of the periodic von neumann basis with biorthogonal exchange (pvb) method to quantum mechanics. Finally, it discusses the application to the time‐independent Schrödinger equation (TISE) and the application to the time‐dependent Schrödinger equation.

We consider manipulation of the transmission coefficient for a quantum particle moving in one dimension where the shape of the potential is taken as the control. We show that the control landscape, the transmission as a functional of the potential, has no traps, i.e., any maxima correspond to full transmission.

Three new developments are presented regarding the semiclassical coherent state propagator. First, we present a conceptually different derivation of Huber and Heller's method for identifying complex root trajectories and their equations of motion [D. Huber and E. J. Heller, J. Chem. Phys. 87, 5302 (1987)]. Our method proceeds directly from the time-dependent Schrodinger equation and therefore allows various generalizations of the formalism. Second, we obtain an analytic expression for the semiclassical coherent state propagator. We show that the prefactor can be expressed in a form that requires solving significantly fewer equations of motion than in alternative expressions. Third, the semiclassical coherent state propagator is used to formulate a final value representation of the time-dependent wavefunction that avoids the root search, eliminates problems with caustics and automatically includes interference. We present numerical results for the 1D Morse oscillator showing that the method may become an attractive alternative to existing semiclassical approaches. (C) 2014 AIP Publishing LLC.

### 2012

We extend a recently developed quantum trajectory method [Y. Goldfarb, I. Degani, and D. J. Tannor, J. Chem. Phys. 125, 231103 (2006)] to treat non-adiabatic transitions. Each trajectory evolves on a single surface according to Newton's laws with complex positions and momenta. The transfer of amplitude between surfaces stems naturally from the equations of motion, without the need for surface hopping. In this paper we derive the equations of motion and show results in the diabatic representation, which is rarely used in trajectory methods for calculating non-adiabatic dynamics. We apply our method to the first two benchmark models introduced by Tully [J. Chem. Phys. 93, 1061 (1990)]. Besides giving the probability branching ratios between the surfaces, the method also allows the reconstruction of the time-dependent wavepacket. Our results are in quantitative agreement with converged quantum mechanical calculations. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4739845]

We present a complex quantum trajectory method for treating non-adiabatic dynamics. Each trajectory evolves classically on a single electronic surface but with complex position and momentum. The equations of motion are derived directly from the time-dependent Schrodinger equation, and the population exchange arises naturally from amplitude-transfer terms. In this paper the equations of motion are derived in the adiabatic representation to complement our work in the diabatic representation [N. Zamstein and D. J. Tannor, J. Chem. Phys. 137, 22A517 (2012)]. We apply our method to two benchmark models introduced by John Tully [J. Chem. Phys. 93, 1061 (1990)], and get very good agreement with converged quantum-mechanical calculations. Specifically, we show that decoherence (spatial separation of wavepackets on different surfaces) is already contained in the equations of motion and does not require ad hoc augmentation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4739846]

Recently we introduced a phase space approach for solving the time-independent Schrodinger equation using a periodic von Neumann basis with bi-orthogonal exchange (pvb) [A. Shimshovitz and D. J. Tannor, Phys. Rev. Lett. 109, 070402 (2012)]. Here we extend the approach to allow a wavelet scaling of the phase space Gaussians. The new basis set, which we call the wavelet pvb basis, is simple to implement and provides an appealing alternative to other wavelet approaches. For the 1D Coulomb problems tested in this paper, the method reduces the size of the basis relative to the Fourier grid method by a factor of 13-60. The savings in basis set size is predicted to grow steeply as the dimensionality increases. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751484]

We propose a method for solving the time-independent Schrodinger equation based on the von Neumann (vN) lattice of phase space Gaussians. By incorporating periodic boundary conditions into the vN lattice [F. Dimler et al., New J. Phys. 11, 105052 (2009)], we solve a longstanding problem of convergence of the vN method. This opens the door to tailoring quantum calculations to the underlying classical phase space structure while retaining the accuracy of the Fourier grid basis. The method has the potential to provide enormous numerical savings as the dimensionality increases. In the classical limit, the method reaches the remarkable efficiency of one basis function per one eigenstate. We illustrate the method for a challenging two-dimensional potential where the Fourier grid method breaks down.

We propose a phase space method to propagate a quantum wavepacket driven by a strong external field. The method employs the periodic von Neumann basis with biorthogonal exchange recently introduced for the calculation of the energy eigenstates of time-independent quantum systems [A. Shimshovitz and D. J. Tannor, Phys. Rev. Lett. (in press) [e-print arXiv:1201.2299v1]]. While the individual elements in this basis set are time-independent, a small subset is chosen in a time-dependent manner to adapt to the evolution of the wavepacket in phase space. We demonstrate the accuracy and efficiency of the present propagation method by calculating the electronic wavepacket in a one-dimensional soft-core atom interacting with a superposition of an intense, few-cycle, near-infrared laser pulse and an attosecond extreme-ultraviolet laser pulse. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4732306]

We have recently proposed a methodology for reconstructing excited-state (ExS) molecular wavepackets, and the corresponding potential energy surface, from three-pulse resonant coherent anti-Stokes Raman scattering and knowledge of the ground-state potential [Avisar and Tannor, Phys. Rev. Lett. 106, 170405 (2011)]. The methodology is general for polyatomics and applies to any form of ExS potential - bound or dissociative. In our previous work we demonstrated the method on diatomics. Here, we demonstrate the method on the triatomics H2O and HOD, reconstructing the ExS wavepacket and potential in the two bond-stretching coordinates. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4722648]

We show that the second-order traps in the control landscape for a three-level ?-system found in our previous work (Phys. Rev. Lett. 2011, 106, 120402) are not local maxima: there exist directions in the space of controls in which the objective grows. The growth of the objective is slow at best 4th order for weak variations of the control. This implies that simple gradient methods would be problematic in the vicinity of second-order traps, where more sophisticated algorithms that exploit the higher order derivative information are necessary to climb up the control landscape efficiently. The theory is supported by a numerical investigation of the landscape in the vicinity of the e(t)=0 second-order trap, performed using the GRAPE and BFGS algorithms.

We consider the optimal control problem of transferring population between states of a quantum system where the coupling proceeds only via intermediate states that are subject to decay. We pose the question whether it is generally possible to carry out this transfer. For a single intermediate decaying state, we recover the stimulated Raman adiabatic passage process, which we identify as the global optimum in the limit of infinite control time. We also present analytical solutions for the case of transfer that has to proceed via two consecutive intermediate decaying states. We show that in this case, for finite power the optimal control does not approach perfect state transfer even in the infinite time limit. We generalize our findings to characterize the topologies of paths that can be achieved by coherent control under the assumption of finite power. If two or more consecutive states in an N-level chain are subject to decay, complete population transfer with finite-power controls is not possible.

In a recent paper [B. Poirier, Chem. Phys. 370, 4 (2010)], a formulation of quantum mechanics was presented for which the usual wavefunction and Schrodinger equation are replaced with an ensemble of real-valued trajectories satisfying a principle of least action. It was found that the resultant quantum trajectories are those of Bohmian mechanics. In this paper, analogous ideas are applied to Bohmian Mechanics with Complex Action (BOMCA). The standard BOMCA trajectories as previously defined are found not to satisfy an action principle. However, an alternate set of complex equations of motion is derived that does exhibit this desirable property, and an approximate numerical implementation is presented. Exact analytical results are also presented, for Gaussian wavepacket propagation under quadratic potentials.

### 2011

One of the most popular methods for solving numerical optimal control problems is the Krotov method, adapted for quantum control by Tannor and coworkers. The Krotov method has the following three appealing properties: (1) monotonic increase of the objective with iteration number, (2) no requirement for a line search, leading to a significant savings over gradient (first-order) methods, and (3) macrosteps at each iteration, resulting in significantly faster growth of the objective at early iterations than in gradient methods where small steps are required. The principal drawback of the Krotov method is slow convergence at later iterations, which is particularly problematic when high fidelity is desired. We show here that, near convergence, the Krotov method degenerates to a first-order gradient method. We then present a variation on the Krotov method that has all the advantages of the original Krotov method but with significantly enhanced convergence (second-order or quasi-Newton) as the optimal solution is approached. We illustrate the method by controlling the three-dimensional dynamics of the valence electron in the Na atom.

Probing the real time dynamics of a reacting molecule remains one of the central challenges in chemistry. Here we show how the time-dependent wave function of an excited-state reacting molecule can be completely reconstructed from resonant coherent anti-Stokes Raman spectroscopy. The method assumes knowledge of the ground potential but not of any excited potential. The excited-state potential can in turn be constructed from the wave function. The formulation is general for polyatomics and applies to bound as well as dissociative excited potentials. We demonstrate the method on the Li(2) molecule.

There has been great interest in recent years in quantum control landscapes. Given an objective J that depends on a control field epsilon the dynamical landscape is defined by the properties of the Hessian delta(2)J/delta epsilon(2) at the critical points delta J/delta epsilon = 0. We show that contrary to recent claims in the literature the dynamical control landscape can exhibit trapping behavior due to the existence of special critical points and illustrate this finding with an example of a 3-level Lambda system. This observation can have profound implications for both theoretical and experimental quantum control studies.

Recently we introduced the von Neumann representation as a joint time-frequency description for femtosecond laser pulses. Here we show that the von Neumann basis can be implemented into an evolutionary algorithm for adaptive optimization in coherent control. We perform simulations that demonstrate the efficiency compared to other parametrizations in the frequency domain. We also illustrate pulse-shape simplification by basis-function reduction. Essential structures using the von Neumann basis are retained without losing control performance significantly. In an optical demonstration experiment we show the practicality by producing double pulses with a given time separation. Adaptive control in time-frequency space will be especially valuable for quantum systems requiring specific transition frequencies at definite times.

We have recently shown how the excited-state wavepacket of a polyatomic molecule can be completely reconstructed from resonant coherent anti-Stokes Raman spectroscopy [Avisar and Tannor, Phys. Rev. Lett., 2011, 106, 170405]. The method assumes knowledge of the ground-state potential but not of any excited-state potential, however the latter can be computed once the excited-state wavepacket is known. The formulation applies to dissociative as well as bound excited potentials. We demonstrate the method on the Li-2 molecule with its bound first excited-state as well as with a model dissociative excited state potential. Preliminary results are shown for a model two-dimensional molecular system. The calculations assume constant transition dipole moment (Condon approximation), delta-pulse excitation and a single excited-state potential, but we discuss the implications of removing these assumptions.

Path-integral derivations are presented for two recently developed complex trajectory techniques for the propagation of wave packets: complex WKB and Bohmian mechanics with complex action (BOMCA). The complex WKB technique is derived using a standard saddle-point approximation of the path integral, but taking into account the (h) over bar dependence of both the amplitude and the phase of the initial wave function, thus giving rise to the need for complex classical trajectories. The BOMCA technique is derived using a modification of the saddle-point technique, in which the path integral is approximated by expanding around a near-classical path, chosen so that up to some predetermined order there is no need to add any correction terms to the leading-order approximation. Both complex WKB and BOMCA techniques give the same leading-order approximation; in the complex WKB technique higher accuracy is achieved by adding correction terms, while in the BOMCA technique no additional terms are ever added: higher accuracy is achieved by changing the path along which the original approximation is computed. The path-integral derivation of the methods explains the need to incorporate contributions from more than one trajectory, as observed in previous numerical work. On the other hand, it emerges that the methods provide efficient schemes for computing the higher-order terms in the asymptotic evaluation of path integrals. The understanding we develop of the BOMCA technique suggests that there should exist near-classical trajectories that give exact quantum dynamical results when used in the computation of the path integral keeping just the leading-order term. We also apply our path-integral techniques to give a compact derivation of the semiclassical approximation to the coherent-state propagator.

### 2010

The influence of a dissipative environment on scattering of a particle by a barrier is investigated by using the recently introduced Bohmian mechanics with complex action [J. Chem. Phys. 125, 231103 (2006)]. An extension of this complex trajectory based formalism to include the interaction of the tunneling particle with an environment of harmonic oscillators with a continuous spectral density and at a certain finite temperature allows us to calculate transmission probabilities beyond the weak system bath coupling regime. The results display an increasing tunneling probability for energies below the barrier and a decreased transmission above the barrier due to the coupling. Furthermore, we demonstrate that solutions of a Markovian master equation fail to do so in general.

Recently we introduced the von Neumann representation as a joint time-frequency description for femtosecond laser pulses and suggested its use as a basis for pulse shaping experiments. Here we use the von Neumann basis to represent multidimensional molecular control landscapes, providing insight into the molecular dynamics. We present three kinds of time-frequency phase space scanning procedures based on the von Neumann formalism: variation of intensity, time-frequency phase space position, and/or the relative phase of single subpulses. The shaped pulses produced are characterized via Fourier-transform spectral interferometry. Quantum control is demonstrated on the laser dye IR140 elucidating a time-frequency pump-dump mechanism. (C) 2010 American Institute of Physics. [doi:10.1063/1.3495950]

We analyze in detail the so-called pushing gate for trapped ions, introducing a time-dependent harmonic approximation for the external motion. We show how to extract the average fidelity for the gate from the resulting semiclassical simulations. We characterize and quantify precisely all types of errors coming from the quantum dynamics and reveal that slight nonlinearities in the ion-pushing force can have a dramatic effect on the adiabaticity of gate operation. By means of quantum optimal control techniques, we show how to suppress each of the resulting gate errors in order to reach a high fidelity compatible with scalable fault-tolerant quantum computing.

The non-Markovian master equation is applied to the calculation of reaction rates. Starting from the flux-side correlation function form, we treat both the thermal and real time evolution consistently within second order perturbation theory in the system-bath coupling. It is shown that the non-Markovian dynamics enter formally not only in the time propagation but also in the expressions for the initial system-bath correlations. We show that these initial correlations can have a significant effect on the reaction rate. The method presented, although approximate, is an effective way to calculate reaction rates for weakly coupled systems over a wide range of temperatures. As such it provides a complementary approach to the exact treatment based on the ML-MCTDH method of Craig et al. [1], which serves as reference in this work. (C) 2010 Elsevier B.V. All rights reserved.

### 2009

We recently introduced the von Neumann picture, a joint time-frequency representation, for describing ultrashort laser pulses. The method exploits a discrete phase-space lattice of nonorthogonal Gaussians to represent the pulses; an arbitrary pulse shape can be represented on this lattice in a one-to-one manner. Although the representation was originally defined for signals with an infinite continuous spectrum, it can be adapted to signals with discrete and finite spectrum with great computational savings, provided that discretization and truncation effects are handled with care. In this paper, we present three methods that avoid loss of accuracy due to these effects. The approach has immediate application to the representation and manipulation of femtosecond laser pulses produced by a liquid-crystal mask with a discrete and finite number of pixels.

We combine optimal control theory with the multi-configuration time-dependent Hartree-Fock method to control the dynamics of interacting particles. We use the resulting scheme to optimize state-to-state transitions in a one-dimensional (1D) model of helium and to entangle the external degrees-of-freedom of two rubidium atoms in a 1D optical lattice. Comparisons with optimization results based on the exact solution of the Schrodinger equation show that the scheme can be used to optimize even involved processes in systems consisting of interacting particles in a reliable and efficient way.

It is well known that a finite level quantum system is controllable if and only if the Lie algebra of its generators has full rank. When the rank of the Lie algebra is not full, there is a rich mathematical and physical structure to the subalgebra that to date has been analyzed only in special cases. We show that uncontrollable systems can be classified into reducible and irreducible ones. The irreducible class is the more subtle and can be related to a notion of generalized entanglement. We give a general prescription for revealing irreducible uncontrollable systems: the fundamental representation of su(N), where N is the number of levels, must remain irreducible in the subalgebra of su(N). We illustrate the concepts with a variety of physical examples.

### 2008

We present a significant improvement to a complex time-dependent WKB (CWKB) formulation developed by Boiron and Lombardi [J. Chem. Phys. 108, 3431 (1998)] in which the time-dependent WKB equations are solved along classical trajectories that propagate in complex space. Boiron and Lombardi showed that the method gives very good agreement with the exact quantum mechanical result as long as the wavefunction does not exhibit interference effects such as oscillations and nodes. In this paper, we show that this limitation can be overcome by superposing the contributions of crossing trajectories. Secondly, we demonstrate that the approximation improves when incorporating higher order terms in the expansion. Thirdly, equations of motion for caustics and Stokes lines are implemented to help overcome Stokes discontinuities. These improvements could make the CWKB formulation a competitive alternative to current time-dependent semiclassical methods. (c) 2008 American Institute of Physics.

### 2007

In recent years, the use of joint time-frequency representations to characterize and interpret shaped femtosecond laser pulses has proven to be very useful. However, the number of points in a joint time-frequency representation is daunting as compared with those in either the frequency or time representation. In this article we introduce the use of the von Neumann representation, in which a femtosecond pulse is represented on a discrete lattice of evenly spaced time-frequency points using a non-orthogonal Gaussian basis. We show that the information content in the von Neumann representation using a lattice of root N points in time and root N points in frequency is exactly the same as in a frequency ( or time) array of N points. Explicit formulas are given for the forward and reverse transformation between an N-point frequency signal and the von Neumann representation. We provide numerical examples of the forward and reverse transformation between the two representations for a variety of different pulse shapes; in all cases the original pulse is reconstructed with excellent precision. The von Neumann representation has the interpretational advantages of the Husimi representation but requires a bare minimum number of points and is stably and conveniently inverted; moreover, it avoids the periodic boundary conditions of the Fourier representation. (c) 2007 Optical Society of America.

We present a unified derivation of Bohmian methods that serves as a common starting point for the derivative propagation method (DPM), Bohmian mechanics with complex action (BOMCA), and the zero-velocity complex action method (ZEVCA). The unified derivation begins with the ansatz psi = e(iS/h) where the action (S) is taken to be complex, and the quantum force is obtained by writing a hierarchy of equations of motion for the phase partial derivatives. We demonstrate how different choices of the trajectory velocity field yield different formulations such as DPM, BOMCA, and ZEVCA. The new derivation is used for two purposes. First, it serves as a common basis for comparing the role of the quantum force in the DPM and BOMCA formulations. Second, we use the new derivation to show that superposing the contributions of real, crossing trajectories yields a nodal pattern essentially identical to that of the exact quantum wavefunction. The latter result suggests a promising new approach to deal with the challenging problem of nodes in Bohmian mechanics.

In recent years, intensive effort has gone into developing numerical tools for exact quantum mechanical calculations that are based on Bohmian mechanics. As part of this effort we have recently developed as alternative formulation of Bohmian mechanics in which the quantum action S is taken to be complex [Y. Goldfarb et al., J. Chem. Phys. 125, 231103 (2006)]. In the alternative formulation there is a significant reduction in the magnitude of the quantum force as compared with the conventional Bohmian formulation, at the price of propagating complex trajectories. In this paper we show that Bohmian mechanics with complex action is able to overcome the main computational limitation of conventional Bohmian methods-the propagation of wave functions once nodes set in. In the vicinity of nodes, the quantum force in conventional Bohmian formulations exhibits rapid oscillations that present a severe numerical challenge. We show that within complex Bohmian mechanics, multiple complex initial conditions can lead to the same real final position, allowing for the accurate description of nodes as a sum of the contribution from two or more crossing trajectories. The idea is illustrated on the reflection amplitude from a one-dimensional Eckart barrier. We believe that trajectory crossing, although in contradiction to the conventional Bohmian trajectory interpretation, provides an important new tool for dealing with the nodal problem in Bohmian methods. (C) 2007 American Institute of Physics.

We present a new semiclassical method that yields an approximation to the quantum mechanical wavefunction at a fixed, predetermined position. In the approach, a hierarchy of ODEs are solved along a trajectory with zero velocity. The new approximation is local, both literally and from a quantum mechanical point of view, in the sense that neighboring trajectories do not communicate with each other. The approach is readily extended to imaginary time propagation and is particularly useful for the calculation of quantities where only local information is required. We present two applications: the calculation of tunneling probabilities and the calculation of low energy eigenvalues. In both applications we obtain excellent agrement with the exact quantum mechanics, with a single trajectory propagation. (c) 2007 Elsevier B.V. All rights reserved.

*Analysis And Control Of Ultrafast Photoinduced Reactions*. p. 25-152 (1Springer Series in Chemical Physics). Abstract

### 2006

In recent years there has been a resurgence of interest in Bohmian mechanics as a numerical tool because of its local dynamics, which suggest the possibility of significant computational advantages for the simulation of large quantum systems. However, closer inspection of the Bohmian formulation reveals that the nonlocality of quantum mechanics has not disappeared-it has simply been swept under the rug into the quantum force. In this paper we present a new formulation of Bohmian mechanics in which the quantum action, S, is taken to be complex. This leads to a single equation for complex S, and ultimately complex x and p but there is a reward for this complexification-a significantly higher degree of localization. The quantum force in the new approach vanishes for Gaussian wave packet dynamics, and its effect on barrier tunneling processes is orders of magnitude lower than that of the classical force. In fact, the current method is shown to be a rigorous extension of generalized Gaussian wave packet dynamics to give exact quantum mechanics. We demonstrate tunneling probabilities that are in virtually perfect agreement with the exact quantum mechanics down to 10(-7) calculated from strictly localized quantum trajectories that do not communicate with their neighbors. The new formulation may have significant implications for fundamental quantum mechanics, ranging from the interpretation of non-locality to measures of quantum complexity. (c) 2006 American Institute of Physics.

Finding multidimensional nondirect product discrete variable representations (DVRs) of Hamiltonian operators is one of the long standing challenges in computational quantum mechanics. The concept of a "DVR set" was introduced as a general framework for treating this problem by R. G. Littlejohn, M. Cargo, T. Carrington, Jr., K. A. Mitchell, and B. Poirier (J. Chem. Phys. 2002, 116, 8691). We present a general solution of the problem of calculating multidimensional DVR sets whose points are those of a known cubature formula. As an illustration, we calculate several new nondirect product cubature DVRs on the plane and on the sphere with up to I 10 points. We also discuss simple and potentially very useful finite basis representations (FBRs), based on general (nonproduct) cubatures. Connections are drawn to a novel view on cubature presented by L Degani, J. Schiff, and D. J. Tannor (Num. Math. 2005, 101, 479), in which commuting extensions of coordinate matrices play a central role. Our construction of DVR sets answers a problem left unresolved in the latter paper, namely, the problem of interpreting as function spaces the vector spaces on which commuting extensions act.

The central objective in any quantum computation is the creation of a desired unitary transformation; the mapping that this unitary transformation produces between the input and output states is identified with the computation. In [S.E. Sklarz, D.J. Tannor, arXiv:-quant-ph/0404081 (submitted to PRA) (2004)] it was shown that local control theory can be used to calculate fields that will produce such a desired unitary transformation. In contrast with previous strategies for quantum computing based on optimal control theory, the local control scheme maintains the system within the computational subspace at intermediate times, thereby avoiding unwanted decay processes. In [S.E. Sklarz et al.], the structure of the Hilbert space had a direct sum structure with respect to the computational register and the mediating states. In this paper, we extend the formalism to the important case of a direct product Hilbert space. The final equations for the control algorithm for the two cases are remarkably similar in structure, despite the fact that the derivations are completely different and that in one case the dynamics is in a Hilbert space and in the other case the dynamics is in a Lionville space. As shown in [S.E. Sklarz et al.], the direct sum implementation leads to a computational mechanism based on virtual transitions, and can be viewed as an extension of the principles of Stimulated Raman Adiabatic Passage from state manipulation to evolution operator manipulation. The direct product implementation developed here leads to the intriguing concept of virtual entanglement - computation that exploits second-order transitions that pass through entangled states but that leaves the subsystems nearly separable at all intermediate times. Finally, we speculate on a connection between the algorithm developed here and the concept of decoherence free subspaces. (c) 2005 Elsevier B.V. All rights reserved.

### 2005

The calculation of chemical reaction rates in the condensed phase is a central preoccupation of theoretical chemistry. At low temperatures, quantum-mechanical effects can be significant and even dominant; yet quantum calculations of rate constants are extremely challenging, requiring theories and methods capable of describing quantum evolution in the presence of dissipation. In this paper we present a new approach based on the use of a non-Markovian quantum master equation (NM-QME). As opposed to other approximate quantum methods, the quantum dynamics of the system coordinate is treated exactly; hence there is no loss of accuracy at low temperatures. However, because of the perturbative nature of the NM-QME it breaks down for dimensionless frictions larger than about 0.1. We show that by augmenting the system coordinate with a collective mode of the bath, the regime of validity of the non-Markovian master equation can be extended significantly, up to dimensionless frictions of 0.5 over the entire temperature range. In the energy representation, the scaling goes as the number of levels in the relevant energy range to the third power. This scaling is not prohibitive even for chemical systems with many levels; hence we believe that the current method will find a useful place alongside the existing techniques for calculating quantum condensed-phase rate constants. (c) 2005 American Institute of Physics.

We extend a recently introduced mapping model, which explains the bunching phenomenon in an ion beam resonator for two ions (Geyer and Tannor 2004 J. Phys. B: At. Mol. Opt. Phys. 37 73), to describe the dynamics of the whole ion bunch. We calculate the time delay of the ions from a model of the bunch geometry and find that the bunch takes on a spherical form at the turning points in the electrostatic mirrors. From this condition we derive how the observed bunch length depends on the experimental parameters. We give an interpretation of the criteria for the existence of the bunch, which were derived from the experimental observations by Pedersen et al (2002 Phys. Rev. A 65 042704).

Based on a novel point of view on 1-dimensional Gaussian quadrature, we present a new approach to d-dimensional cubature formulae. It is well known that the nodes of 1-dimensional Gaussian quadrature can be computed as eigenvalues of the so-called Jacobi matrix. The d-dimensional analog is that cubature nodes can be obtained from the eigenvalues of certain mutually commuting matrices. These are obtained by extending (adding rows and columns to) certain non-commuting matrices A(1),...,A(d), related to the coordinate operators x(1),...,x(d), in R-d. We prove a correspondence between cubature formulae and "commuting extensions" of A(1),...,A(d), satisfying a compatibility condition which, in appropriate coordinates, constrains certain blocks in the extended matrices to be zero. Thus, the problem of finding cubature formulae can be transformed to the problem of computing (and then simultaneously diagonalizing) commuting extensions. We give a general discussion of existence and of the expected size of commuting extensions and briefly describe our attempts at computing them.

The Jaynes-Cummings model (JCM) is the simplest fully quantum model that describes the interaction between light and matter. We extend a previous analysis by Phoenix and Knight [Ann. Phys. 186, 381 (1988)] of the JCM by considering mixed states of both the light and matter. We present examples of qualitatively different entropic correlations. In particular, we explore the regime of entropy exchange between light and matter, i.e., where the rate of change of the two are anticorrelated. This behavior contrasts with the case of pure light-matter states in which the rate of change of the two entropies are positively correlated and in fact identical. We give an analytical derivation of the anticorrelation phenomenon and discuss the regime of its validity. Finally, we show a strong correlation between the region of the Bloch sphere characterized by entropy exchange and that characterized by minimal entanglement as measured by the negative eigenvalues of the partially transposed density matrix.

### 2004

We study the problem of optimal control of dissipative quantum dynamics. Although under most circumstances dissipation leads to an increase in entropy (or a decrease in purity) of the system, there is an important class of problems for which dissipation with external control can decrease the entropy (or increase the purity) of the system. An important example is laser cooling. In such systems, there is an interplay of the Hamiltonian part of the dynamics, which is controllable, and the dissipative part of the dynamics, which is uncontrollable. The strategy is to control the Hamiltonian portion of the evolution in such a way that the dissipation causes the purity of the system to increase rather than decrease. The goal of this paper is to find the strategy that leads to maximal purity at the final time. Under the assumption that Hamiltonian control is complete and arbitrarily fast, we provide a general framework by which to calculate optimal cooling strategies. These assumptions lead to a great simplification, in which the control problem can be reformulated in terms of the spectrum of eigenvalues of rho, rather than rho itself. By combining this formulation with the Hamilton-Jacobi-Bellman theorem we are able to obtain an equation for the globally optimal cooling strategy in terms of the spectrum of the density matrix. For the three-level Lambda system, we provide a complete analytic solution for the optimal cooling strategy. For this system it is found that the optimal strategy does not exploit system coherences and is a "greedy" strategy, in which the purity is increased maximally at each instant.

We present a two particle model to explain the mechanism that stabilizes a bunch of positively charged ions in an 'ion trap resonator' (Pedersen et al 2001 Phys. Rev. Lett. 87 055001). The model decomposes the motion of the two ions into two mappings for the free motion in different parts of the trap and one for a compressing momentum kick. The ions' interaction is modelled by a time delay, which then changes the balance between adjacent momentum kicks. Through these mappings we identify the microscopic process that is responsible for synchronization and give the conditions for that regime.

### 2002

We demonstrate that the synchronization effect observed [Pedersen et al., Phys. Rev. Lett. 87, 055001 (2001)], when a bunch of ions oscillates between two mirrors in an electrostatic ion beam trap, can be explained as a negative mass instability. We derive simple necessary conditions for the existence of a regime in which this dispersionless behavior occurs and demonstrate that in this regime, the ion trap can be used as a high resolution mass spectrometer.

Cooling of internal atomic and molecular states via optical pumping and laser cooling of the atomic velocity distribution, rely on spontaneous emission. The outstanding success of such examples, taken together with general arguments, has led to the widely held notion that radiative cooling requires spontaneous emission. We here show by specific examples and direct calculation, based primarily on breaking emission-absorption symmetry as in lasing without inversion, that cooling of internal states by external coherent control fields is possible. We also show that such coherent schemes allow us to practically reach absolute zero in a finite number of steps, in contrast to some statements of the third law of thermodynamics.

It has been proposed that the adiabatic loading of a Bose-Einstein condensate (BEC) into an optical lattice via the Mott-insulator transition can be used to initialize a quantum computer [D. Jaksch , Phys. Rev. Lett. 81, 3108 (1998)]. The loading of a BEC into the lattice without causing band excitation is readily achievable; however, unless one switches on an optical lattice very slowly, the optical lattice causes a phase to accumulate across the condensate. We show analytically and numerically that a cancellation of this effect is possible by adjusting the harmonic trap force constant of the magnetic trap appropriately, thereby facilitating quick loading of an optical lattice for quantum computing purposes. A simple analytical theory is developed for a nonstationary BEC in a harmonic trap.

Using a set of general methods developed by Krotov [A. I. Konnov and V. A. Krotov, Automation and Remote Control 60, 1427 (1999)], we extend the capabilities of optimal control theory to the nonlinear Schrodinger equation (NLSE). The paper begins with a general review of the Krotov approach to optimization. Although the linearized version of the method is sufficient for the linear Schrodinger equation, the full flexibility of the general method is required for treatment of the nonlinear Schrodinger equation. Formal equations for the optimization of the NLSE, as well as a concrete algorithm are presented. As an illustration, we consider a Bose-Einstein condensate (BEC) initially at rest in a harmonic trap. A phase develops across the BEC when an optical lattice potential is turned on. The goal is to counter this effect and keep the phase flat by adjusting the trap strength. The problem is formulated in the language of optimal control theory (OCT) and solved using the above methodology. To our knowledge, this is the first rigorous application of OCT to the nonlinear Schrodinger equation, a capability that is bound to have numerous other applications.

We examine, both numerically and analytically, the limits of controllability and the dynamics of population transfer in systems with degenerate target states embedded in a finite manifold of states. An analysis of the dressed eigenstates of coupled field-matter systems is combined with the theory of the controllability of invariant systems on Lie groups to show the existence of a control scheme using a pair of shaped laser pulses acting upon a four-state system with intermediate degenerate states and a diamond-shaped array of transition moments (mu(12),mu(24),mu(13),mu(34)not equal0,mu(14),mu(23)=0). The limits of the controllability of these systems are linked to the symmetry of the transition dipole moments connecting the states. A broken symmetry of the transition dipole moments (mu(12)/mu(24)not equalmu(13)/mu(34)) permits complete population transfer to one of the pair of degenerate states. Numerical simulations of the optimized population transfer to one of a pair of degenerate states in the system agree quantitatively with the analytically predicted limits of population transfer in different regimes in the joint material-field characteristics of the system. The applicability of the four-state model to selective population transfer in vibrationally mediated photodissociation experiments is discussed.

### 2001

Knopp [J. Raman Spectrosc. 31, 51 (2000)] have recently used resonant femtosecond coherent anti-Stokes Raman spectroscopy (CARS) to prepare and probe highly excited vibrational wave packets on the ground electronic potential surface of molecular iodine. The experiment uses a sequence of three resonant femtosecond pulses with two independently variable time delays. The first two pulses act as a pump and dump sequence to create a predefined, highly excited wave packet on the ground electronic state, whose amplitude is optimized by selecting the proper pump-dump (Raman) frequency difference and varying the time delay. The third pulse promotes the pump-dump wave packet to an excited electronic state, resulting in subsequent coherent emission of light at the anti-Stokes frequency. This fully-resonant CARS signal, measured as a function of time delay between the second and third pulses, oscillates at a frequency characteristic of the pump-dump wave packet. Due to anharmonicity, this frequency is a sensitive measure of the amount of vibrational excitation. Knopp observed that under certain conditions the signal exhibits pronounced beating between the pump-dump wave packet frequency and the frequency characteristic of the bottom of the ground state well. In this paper we show that these beats arise only when the final pump-dump-pump wave packet is above the excited state dissociation threshold of the molecule. We derive analytical expressions showing that under these conditions, where the polarization is short-lived, there may be strong interferences between the contributions from molecules originally in different vibrational states of the thermal ensemble. In contrast, the CARS polarization in the below threshold case is long-lived, and these interferences cancel. Numerical evaluation of the CARS signal through vibrational wave packet propagation confirms the predictions of the analytical theory and reproduces the distinctive beating pattern observed in the experiments.

Optimal control theory (OCT) is applied to laser cooling of molecules. The objective is to cool vibrations, using shaped pulses synchronized with the spontaneous emission. An instantaneous in time optimal approach is compared to solution based on OCT. In both cases the optimal mechanism is found to operate by a "vibrationally selective coherent population trapping". The trapping condition is that the instantaneous phase of the laser is locked to the phase of the transition dipole moment of nu = 0 with the excited population. The molecules that reach a = 0 by spontaneous emission are then trapped. while the others are continually repumped. For vibrational cooling to nu = 2 and rotational cooling, a different mechanism operates. The field completely changes the transient eigenstates of the Hamiltonian creating a superposition composed of many states. Finally this superposition is transformed by the field to the target energy eigenstate. (C) 2001 Elsevier Science B.V. All rights reserved.

### 2000

We consider the equilibrium state of a quantum system weakly coupled to a quantum bath within second order perturbation theory. It was previously shown by Romero-Rochin and Oppenheim [Physica A 155, 52 (1989)] that the equilibrium state deviates from the canonical form, e(-beta Hs)/Z(s) (H(s) is the free system Hamiltonian and Z(s) the canonical partition function). We reproduce this result via a different derivation, starting from the non-Markovian, rather than the Markovian, quantum Master equation. Our derivation sheds new light on the mechanism that stabilizes the deviation from the canonical form and shows that it involves an interplay between a static distortion to the equilibrium state and dynamical system-bath correlations. We show that this deviation is a necessary consequence of translational invariance and vanishes when the rotating-wave-approximation is applied. The deviation is also shown to vanish for a two-level system off-diagonally coupled to a heat bath or when the Lamb shifts are neglected. Two ways for numerically evaluating the second order deviations are described. Finally, the deviations from canonical equilibrium are given an illuminating geometrical interpretation in terms of the phase space Wigner distribution. (C) 2000 American Institute of Physics. [S0021-9606(00)50828-9].

Calculation of chemical reaction dynamics is central to theoretical chemistry. The majority of calculations use either classical mechanics, which is computationally inexpensive but misses quantum effects, such as tunneling and interference, or quantum mechanics, which is computationally expensive and often conceptually opaque. An appealing middle ground is the use of semiclassical mechanics. Indeed, since the early 1970s there has been great interest in using semiclassical methods to calculate reaction probabilities. However, despite the elegance of classical S-matrix theory, numerical results on even the simplest reactive systems remained out of reach. Recently, with advances both in correlation function formulations of reactive scattering as well as in semiclassical methods, it has become possible for the first time to calculate reaction probabilities semiclassically. The correlation function methods are contrasted with recent flux-based methods, which, although providing somewhat more compact expressions for the cumulative reactive probability, are less compatible with semiclassical implementation.

### 1999

Optimal control theory (OCT) is applied to the problem of cooling molecular rotations. The optimal field gives rise to a striking behavior, in which there is no noticeable increase in the lowest rotational state population until the last percent or so of the control interval,:at which point the population jumps to 1. Further analysis of the intermediate time interval reveals that cooling is taking place all along, in the sense that the purity of the system, as measured by Tr(rho(2)), is increasing monotonically in time. Once the system becomes almost completely pure, the external control field can transfer the amplitude to the lowest rotational state by a completely Hamiltonian manipulation. This mechanism is interesting because it suggests a possible way of accelerating cooling, by exploiting the cooling induced by spontaneous emission to all the ground electronic state levels, not just the lowest rotational level. However, it also raises a major paradox: it may be shown that external control fields, no matter how complicated, cannot change the value of Tr(rho(2)); changing this quantity requires spontaneous emission which is inherently uncontrollable. What place is there then for control, let alone optimal control, using external fields? We discuss the resolution to this paradox with a detailed analysis of cooling in a two-level system.

We study different mechanisms of adiabatic population transfer in N-level systems by means of optimal control algorithms. Using two-dimensional topographic maps of the yield of population transfer as a function of time delay and intensity of the pulses we analyze the global properties of the schemes and the conditions that lead to optimization. For three-level systems it is shown that the optimal pulse sequence is the well-known STIRAP (stimulated Raman adiabatic passage) scheme. For five-level systems a family of solutions ranging from the alternating STIRAP scheme to the new straddling STIRAP (S-STIRAP) scheme is obtained and the behavior of the solutions is compared. For four-level systems we obtain as optimal a S-STIRAP type sequence that behaves as an effective two-level system. For both odd and even numbers of N-level systems, the crucial role of the straddling pulse in reducing the population of all intermediate levels is demonstrated. [S1050-2947(99)04509-6].

A numerical method is described for integration of the time-dependent Schrodinger equation within the presence of a Coulomb held. Because of the singularity at r = 0, the wave packet has to be represented on a grid with a high density of points near the origin; at the same timet because of the long-range character of the Coulomb potential, the grid must extend to large values of r. The sampling points are chosen, following E. Fattal, R. Baer, and R. Kosloff [Phys. Rev. E 53, 1217 (1996)], using a classical phase space criterion. Following those workers, the unequally spaced grid points are mapped to an equally spaced grid, allowing use of fast Fourier transform propagation methods that scale as N ln N, where N is the number of grid points. As a first test, eigenenergies for the hydrogen atom are extracted from short-time segments of the electronic wave-packet autocorrelation function; high accuracy is obtained by using the filter-diagonalization method,As a second test, the ionization rate of the hydrogen atom resulting from a half-cycle pulse is calculated. These results are in excellent agreement with earlier calculations. [S1063-651X(99)14210-7].

We present an accurate, efficient, and flexible method for propagating spatially distributed density matrices in anharmonic potentials interacting with solvent and strong fields. The method is based on the Nakajima-Zwanzig projection operator formalism with a correlated reference state of the bath that takes memory effects and initial/final correlations to second order in the system-bath interaction into account. A key feature of the method proposed is a special parametrization of the bath spectral density leading to a set of coupled equations for primary and N auxiliary density matrices. These coupled master equations can be solved numerically by representing the density operator in eigenrepresentation or on a coordinate space grid, using the Fourier method to calculate the action of the kinetic and potential energy operators, and a combination of split operator and Cayley implicit method to compute the time evolution. The key advantages of the method are: (1) The system potential may consist of any number of electronic states, either bound or dissociative. (2) The cost for arbitrarily long solvent memories is equal to only N + 1 times that of propagating a Markovian density matrix. (3) The method can treat explicitly time-dependent system Hamiltonians nonperturbatively, making the method applicable to strong field spectroscopy, photodissociation, and coherent control in a solvent surrounding. (4) The method is not restricted to special forms of system-bath interactions. Choosing as an illustrative example the asymmetric two-level system, we compare our numerical results with full path-integral results and we show the importance of initial correlations and the effects of strong fields onto the relaxation. Contrary to a Markovian theory, our method incorporates memory effects, correlations in the initial and final state, and effects of strong fields onto the relaxation; and is yet much more effective than path integral calculations. It is thus well-suited to stud

At the heart of the search for a quantum transition state theory is the partitioning of dynamic from thermal factors in the quantum rate expression. We explore the possibility of achieving an approximate partitioning by using the coherent state basis. The coherent states provide a tetradic representation of both the dynamic and thermal factors; the degree to which these factors partition is tied to the degree to which one or both of the tetradics is diagonal, and hence phase space localized. We find that for the dynamical factor the off-diagonal contributions, are small, except for matrix elements between coherent states positioned anywhere along the stable branch of the classical separatrix. The thermal factor is nearly diagonal at high temperatures, but has significant off-diagonal contributions at low temperatures. As a result, at high temperatures the thermal factor cuts off long range correlation, leading to the classical limit. At low temperatures, there is a subtle interplay of the thermal and dynamical factors, with the long range off-diagonal portions of the thermal factor combining with the long range off-diagonal portions of the dynamical factor. This phase space picture sheds light on the physical assumptions underlying several commonly applied approximations for calculating thermal reaction rates. In particular, by elucidating the subtlety of the contributions to the low temperature rate it becomes clear why a simple, yet accurate, estimate of the rate in, this regime is elusive, if not impossible.

Calculation of chemical reaction rates lies at the very core of theoretical chemistry. The essential dynamical quantity which determines the reaction rate is the energy-dependent cumulative reaction probability, N(E), whose Boltzmann average gives the thermal rate constant, k(T). Converged quantum mechanical calculations of N(E) remain a challenge even for three- and four-atom systems, and a longstanding goal of theoreticians has been to calculate NO accurately and efficiently using semiclassical methods. In this article we present a variety of methods for achieving this goal, by combining semiclassical initial value propagation methods with a reactant-product wavepacket correlation function approach to reactive scattering. The correlation function approach, originally developed for transitions between asymptotic internal states of reactants and products, is here reformulated using wavepackets in an arbitrary basis, so that N(E) can be calculated entirely from trajectory dynamics in the vicinity of the transition state. This is analogous to the approaches pioneered by Miller for the quantum calculation of N(E), and leads to a reduction in the number of trajectories and the propagation time. Numerical examples are presented for both one-dimensional test problems and for the collinear hydrogen exchange reaction.

We present new expressions for the cumulative reaction probability (N(E)), cast in terms of time-correlation functions of reactant and product wave packets. The derivation begins with a standard trace expression for the cumulative reaction probability, expressed in terms of the reactive scattering matrix elements in an asymptotic internal basis. By combining the property of invariance of the trace with a wave packet correlation function formulation of reactive scattering, we obtain an expression for N( E) in terms of the correlation matrices of incoming and outgoing wave packets which are arbitrary in the internal coordinates. This formulation, like other recent formulations of N(E), allows calculation of the quantum dynamics just in the interaction region of the potential, and removes the need for knowledge of the asymptotic eigenstates. However, unlike earlier formulations, the present formulation is fully compatible with both exact and approximate methods of wave packet propagation. We illustrate this by calculating N(E) for the collinear hydrogen exchange reaction, both quantally and semiclassically. These results indicate that the use of wave packet cross-correlation functions, as opposed to a coordinate basis and flux operators, regularizes the semiclassical calculation, suggesting that the semiclassical implementation described here may be applied fruitfully to systems with more degrees of freedom. (C) 1999 American Institute of Physics. [S0021-9606(99)01706-7].

In the last several years we have discovered a variety of remarkable pulse strategies for manipulating molecular motion by employing a design strategy we call "local optimization.'' Here we review the concept of local optimization and contrast it with optimal control theory. By way of background, we give highlights from two recent examples of the method: (1) a strategy for eliminating population transfer to one or many excited electronic states during strong field excitation, an effect we call 'optical paralysis'; (2) a generalization of the counterintuitive STIRAP (stimulated Raman adiabatic passage) pulse sequence from three levels to N levels, a strategy we call 'straddling STIRAP.' We then turn to a third example, which is the main subject of this paper: laser cooling of molecular internal degrees of freedom. We study a model that includes both coherent interaction with the radiation field and spontaneous emission; the latter is necessary to carry away the entropy from the molecule. An optimal control calculation was performed first and succeeded in producing vibrational cooling, but the resulting pulse sequence was difficult to interpret. Local optimization subsequently revealed the cooling mechanism: the instantaneous phase of the laser is locked to the phase of the transition dipole moment between the excited state amplitude and v=0 of the ground state. Thus, the molecules that reach v=0 by spontaneous emission become decoupled from the field, and no longer absorb, while molecules in all other states are continually repumped. The mechanism could be called "vibrationally selective coherent population trapping,'' in analogy to the corresponding mechanism of velocity selective coherent population trapping in atoms for sub-Doppler cooling of translations.

### 1998

In a recent series of papers we showed how a phase space distribution function approach could be combined with the method of reactive flux to obtain the rate of barrier crossing as a function of solvent friction in the intermediate to high friction limit. Those studies dealt with both the Markovian and non-Markovian cases, but were restricted to analytic results for parabolic barriers. Here we extend the approach to anharmonic barriers. The guiding approximation is to assume that the phase space distribution for each initial velocity, starting at the barrier top, remains Gaussian for all time, with Gaussian parameters given by time-dependent mean field equations. We expect this approximation to be accurate for short times, up to the "Ehrenfest" time; if this time exceeds the "plateau" time - the time for the distribution to reach its asymptotic partitioning - the quality of the results should be high. There are no adjustable parameters, although some reasonable criterion is needed for ending the integration of the mean field equations to prevent divergence. Numerical results for the linear cusp and the quartic potential show that the method is quite accurate for dimensionless frictions Il, although the accuracy degrades for higher frictions. (C) 1998 Elsevier Science B.V. All rights reserved.

A correlation function formulation for the state-selected total reaction probability, N-alpha(E), is Suggested. A wave packet, correlating with a specific set of internal reactant quantum numbers, alpha, is propagated forward in time until bifurcation is complete at which time the nonreactive portion of the amplitude is discarded. The autocorrelation function of the remaining amplitude is then computed and Fourier transformed to obtain a reactivity spectrum. Dividing by the corresponding spectrum of the original,unfiltered, wave packet normalizes the reactivity spectrum, yielding the fetal reaction probability from the internal state, alpha. The procedure requires negligible storage and just one time-energy Fourier transform for each initial reactant state; independent of the number of open channels of products. The method is illustrated numerically for the one-dimensional Eckart barrier, using both quantum-mechanical and semiclassical propagation methods. Summing over internal states of reactants gives the cumulative reaction probability, N(E). The relation to the trace formula [W. H. Miller, S. D. Schwartz, J.W. Tromp, J. Chem. Phys. 79, 4889 (1983)], N(E) = 1/2(2 pi (h) over bar)(2) tr((F) over bar delta(H-E)(F) over bar-delta(H-E)), is established, and a new variant of the trace formula is presented. (C) 1998 American Institute of Physics. [S0021-9606(98)01127-1].

Optimal control theory (OCT) applied to driving molecular systems by means of femtosecond pulses is now a mature area, but many of its intricacies are as yet unexplored. As a numerical tool, the many variations on the basic method differ not only in computer efficiency but in the type of solutions obtained. In this paper we survey this diversity, focusing on the use of multiphoton IR laser excitation to control either (1) the state selectivity or (2) the photodissociation in a ID Morse potential. We compare two distinct algorithms, the Krotov method and the gradient method. The former method generates large changes in the field at each iteration, white the latter does not. As a result, the Krotov method virtually always leads to pulses that are very different from the initial guess, while with the gradient method this is not always the case. We then analyze the effect of changing the final time, T, and find that it also can have a profound effect on the nature of the optimal solutions. Finally, we compare the solutions obtained using two different projectors to describe the bond-breaking process: a coordinate projector and a projector over scattering states. Again we observe that the optimal pulses and the dynamics they generate are markedly different in the two cases. This ambiguity in the definition of the optimal pulses may be viewed as a shortcoming of the approach, or alternatively it may be viewed as giving the method extra flexibility.

### 1997

STIRAP (stimulated Raman adiabatic passage) has proven to be an efficient and robust technique for transferring population in a three-level system without populating the intermediate state. Here we show that the counterintuitive pulse sequence in STIRAP, in which the Stokes pulse precedes the pump, emerges automatically from a variant of optimal control theory we have previously called ''local'' optimization. Since local optimization is a well-defined, automated computational procedure, this opens the door to automated computation of generalized STIRAP schemes in arbitrarily complicated N-level coupling situations. If the coupling is sequential, a simple qualitative extension of STIRAP emerges: the Stokes pulse precedes the pump as in the three-level system. But, in addition, spanning both the Stokes and pump pulses are pulses corresponding to the transitions between the N-2 intermediate states with intensities about an order of magnitude greater than those of the Stokes and pump pulses. This scheme is amazingly robust, leading to almost 100% population transfer with significantly less population transfer to the N-2 intermediate states than in previously proposed extensions of STIRAP. [S1050-2947(97)07611-7].

In a recent paper we showed the equivalence, under certain well-characterized assumptions, of Redfield's equations for the density operator in the energy representation with the Gaussian phase space ansatz for the Wigner function of Yan and Mukamel. The equivalence shows that the solutions of Redfield's equations respect a striking degree of classical-quantum correspondence, Here we use this equivalence to derive analytic expressions for the density matrix of the harmonic oscillator in the energy representation without making the almost ubiquitous secular approximation. From the elements of the density matrix in the energy representation we derive analytic expressions for Gamma(1)(n)(1/T-1(n)) and Gamma(2)(nm)(1/T-2(nm)), i.e., population and phase relaxation rates for individual matrix elements in the energy representation. Our results show that Gamma(1)(n)(t) = Gamma(1)(t) is independent of n; this is contrary to the widely held belief that Gamma(1)(n) is proportional to n. We also derive the simple result that Gamma(2)(nm)(t) = \n-m\Gamma(1)(t)/2, a generalization of the two-level system result Gamma(2) = Gamma(1)/2. We show that Gamma(1)(t) is the classical rate of energy relaxation, which has periodic modulations characteristic of the classical damped oscillator; averaged over a period Gamma(t) is directly proportional to the classical friction, gamma. An additional element of classical-quantum correspondence concerns the time rate of change of the phase of the off diagonal elements of the density matrix, omega(nm), a quantity which has received little attention previously. We find that omega(nm) is time-dependent, and equal to \n - m\Omega(t), where Omega(t) is the rate of change of phase space angle in the classical damped harmonic oscillator. Finally, expressions for a collective Gamma(1)(t) and Gamma(2)(t) are derived, and shown to satisfy the relationship Gamma(2) = Gamma(1)/2. This familiar result, when applied to these collective rate constants, is se

Six major theories of quantum dissipative dynamics are compared: Redfield theory, the Gaussian phase space ansatz of Yan and Mukamel, the master equations of Agarwal, Caldeira-Leggett/Oppenheim-Romero-Rochin, and Louisell/Lax, and the semigroup theory of Lindblad. The time evolving density operator from each theory is transformed into a Wigner phase space distribution, and classical-quantum correspondence is investigated via comparison with the phase space distribution of the classical Fokker-Planck (FP) equation. Although the comparison is for the specific case of Markovian dynamics of the damped harmonic oscillator with no pure dephasing, certain inferences can be drawn about general systems, The following are our major conclusions: (1) The harmonic oscillator master equation derived from Redfield theory, in the limit of a classical bath, is identical to the Agarwal master equation. (2) Following Agarwal, the Agarwal master equation can be transformed to phase space, and differs from the classical FP equation only by a zero point energy in the diffusion coefficient. This analytic solution supports Gaussian solutions with the following properties: the differential equations for the first moments in p and q and all but one of the second moments (q(2) and pq but not p(2)) are identical to the classical equations. Moreover, the distribution evolves to the thermal state of the bare quantum system at lone times. (3) The Gaussian phase space ansatz of Yan and Mukamel (YM), applied to single surface oscillator dynamics, reduces to the analytical Gaussian solutions of the Agarwal phase space master equation. It follows that the YM ansatz is also a solution to the Redfield master equation. (4) The Agarwal/Redfield master equation has a structure identical to that of the master equation of Caldeira-Leggett/ Oppenheim-Romero-Rockin, but the two are equivalent only in the high temperature limit. (5) The Louisell/Lax HO master equation differs from the Agarwal/Redfield form

The possibility of using phase coherent optical pulse sequences to generate large-amplitude vibrational motion while locking excited state population has been demonstrated by us previously [Kosloff et al., Phys. Rev. Lett. 69 (1992) 2172], While the approach is compatible in principle with strong fields, in practice strong fields excite many higher electronic states and may produce multiphoton ionization, processes that are frequently neglected in model calculations. Here we demonstrate that it is possible to lock any number of unwanted electronic excitations by a single condition on the instantaneous phase of the pulse sequence. We call this scheme ''optical paralysis''. Since only the phase of the field is determined by this condition, the amplitude of the field is still unspecified and can be chosen to achieve some desired objective, e.g. monotonic increase of the ground-state vibrational energy. The scheme is demonstrated by solving the time-dependent Schrodinger equation for nine coupled electronic states of Na-2, with energy deposited in the ground state and a single excited state while the population in all other excited states is kept locked. (C) 1997 Elsevier Science B.V.

Scattering matrix elements and symmetric transition-state resonances for the collinear H-2 + H --> H + H-2 reaction are obtained using a time-dependent approach. The correlation function between reactant channel wavepackets and product channel wavepackets is used to determine the S-matrix elements. In a similar fashion, autocorrelation functions are used to extract the positions and widths of transition-state resonances. The time propagation of the wavepackets is performed by the improved semiclassical frozen Gaussian method of Herman and Kluk, which is an initial value, uniformly converged method. The agreement between the quantum and semiclassical results is far better than that obtained previously for this system by other semiclassical methods.

Theoretical progress in the cooling of internal degrees of freedom of molecules using shaped laser pulses is reported. The emphasis is on general concepts and universal constraints. Several alternative definitions of cooling are considered, including reduction of the von Neumann entropy, -tr{} and increase of the Renyi entropy, tr{(2)}. A distinction between intensive and extensive considerations is used to analyse the cooling process in open systems. It is shown that the Renyi entropy increase is consistent with an increase in the system phase space density and an increase in the absolute population in the ground state. The limitations on cooling processes imposed by Hamiltonian generated unitary transformations are analyzed. For a single mode system with a ground and excited electronic surfaces driven by an external field it is shown that it is impossible to increase the ground state population beyond its initial value. A numerical example based on optimal control theory demonstrates this result. For this model only intensive cooling is possible which can be classified as evaporative cooling. To overcome this constraint, a single bath degree of freedom is added to the model. This allows a heat pump mechanism in which entropy is pumped by the radiation from the primary degree of freedom to the bath mode, resulting in extensive cooling. (C) 1997 American Institute of Physics.

### 1996

A procedure is described for calculating both the magnitude and phase of an autocorrelation function from the Wigner distribution. The method is fully compatible with both quantum and classical-based procedures for propagating the Wigner distribution. Numerical results an presented for an initial Gaussian in a considerably anhamionic Morse potential, showing that the method is accurate for several vibrational periods.

Scattering matrix elements for the collinear H-2(nu)+H-->H+H-2(nu') reaction are obtained using a recently developed time-dependent approach to scattering. The correlation function between reactant channel wave packet and product channel wave packet is used to determine the S-matrix elements. The time propagation of the reactant wave packet is performed by the semiclassical method of Herman and Kluk, which is an initial value, uniformly converged method. The agreement between the quantum and semiclassical results is far better than that obtained previously for this system by other semiclassical methods.