Semiclassical initial value representation for the Boltzmann operator in thermal rate constants
|
|
- Terence Morris
- 5 years ago
- Views:
Transcription
1 JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 21 1 DECEMBER 22 Semiclassical initial value representation for the Boltzmann operator in thermal rate constants Yi Zhao and William H. Miller a) Department of Chemistry and Kenneth S. Pitzer Center for Theoretical Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 9472 Received 11 June 22; accepted 5 September 22 The thermal rate constant for a chemical reaction, k(t), can be expressed as the long time limit of the flux-side correlation C fs (t)tre Ĥ/2 Fˆ e Ĥ/2 e iĥt/ ĥe iĥt/. Previous work has focused on semiclassical SC approximations implemented via an initial value representation IVR for the time evolution operators exp(iĥt/) in the correlation function, and this paper shows how an SC-IVR can also be used to approximate the Boltzmann operators exp(ĥ/2). Test calculations show that over a wide temperature range little error is introduced in the rate constant by this SC approximation for the Boltzmann operator. 22 American Institute of Physics. DOI: 1.163/ I. INTRODUCTION This paper is a continuation of our efforts to develop the semiclassical SC initial value representation IVR Refs. 1 3 into a practical way of adding quantum effects to classical molecular dynamics simulations, here focused on calculating thermal rate constants for chemical reactions. A formally exact quantum mechanical expression for such rate constants is the long time limit of the flux-side correlation function, 4,5 ktq r T 1 lim C fs t, t which is given by C fs ttrfˆ ĥt, where Fˆ () is the Boltzmannized flux operator Fˆ e Ĥ/2 i Ĥ,hsqˆ e Ĥ/2, and ĥ(t) is the time-evolved projection operator onto products, ĥte iĥt/ hsqˆ e iĥt/. 1.4 s(q above is some function of the coordinates of the system that defines a dividing surface via the equation s(q), and Q r is the reactant partition function per unit volume. There is also a Kubo version 6,7 of the Boltzmannized flux operator that is sometimes used, 8 but we have found the above symmetrically split version to behave better in the numerical calculations presented in this paper. Evaluation of rate constants via Eqs thus involves two essential operators, ĥ(t) and Fˆ (), the former involving real time propagators, exp(iĥt/), and the latter the imaginary time propagator or Boltzmann operator, exp(ĥ/2). Previous work in our group and also by others on applying SC-IVR methods to evaluate rate constants via Eqs has concentrated on various ways of representing the real time propagators in ĥ(t), e.g., a linearization 9 of the difference between the two propagators in Eq. 1.4 leads to the classical Wigner model 1 15 and also what Pollak et al refer to as mixed quantum classical rate theory; the forward backward IVR Refs is a more accurate way for combining the two propagators in ĥ(t) into one effective propagator from to t and from t back to ; and the most complete semiclassical treatment is to use the full SC-IVR for both propagators independently. 27 With these various ways of treating the real time propagators in ĥ(t), the imaginary time propagators in Fˆ (), Eq. 1.3, have typically been dealt with by a harmonic approximation 23 about the transition state or essentially exactly by a Feynman path integral representation 28 evaluated by Monte Carlo methods 27,29. This paper is concerned with evaluating the Boltzmann operators in Eq. 1.3 also via SC-IVR. This will be more accurate than a harmonic approximation which has been noted to break down at low temperature and presumably more efficient than a path integral treatment. Also, it seems more esthetically pleasing and consistent to use the same kind of approximation for the operator exp(ĥ/2) as for exp(iĥt/), since they are essentially the same operator. We note that Makri and Miller 3 have recently explored use of a coherent state i.e., Herman Kluk-type 31 IVR for the Boltzmann operator, but in this paper we employ the coordinate space or Van Vleck 32 IVR. As one of us has recently shown, 33 a Gaussian filtering specifically, a modified Filinov transformation of the coordinate space IVR converts it into the coherent state IVR, so the two IVR s are effectively equivalent. We find it simpler, at least for present purposes, to formulate things in terms of the coordinate space IVR, and then to carry out any necessary filtering at the end. We also note that Pollak and Eckhardt 37 have used the SC approxia Electronic mail: miller@cchem.berkeley.edu /22/117(21)/965/6/$ American Institute of Physics
2 966 J. Chem. Phys., Vol. 117, No. 21, 1 December 22 Y. Zhao and W. H. Miller mation for the Boltzmann operator in Eq. 1.3 but incorporated further SC approximations e.g., steepest descent in evaluating the correlation function. Section II first describes the SC-IVR of the Boltzmann operator which is based on earlier work 38 of one of us and how it is used in Eqs Some numerical tests are then presented in Sec. III to show how well the SC approximation for the Boltzmann operator performs, and Sec. IV concludes. II. SEMICLASSICAL BOLTZMANNIZED FLUX OPERATOR A. The SC-IVR of the Boltzmann operator The standard Van Vleck SC approximation 32 for matrix elements of the time evolution operator i.e., real time propagator is q 2 e iĥt/ q 1 2i q F 2 p 1 e is t (q 2,q 1 )/, 2.1a where the sum is over all classical trajectories that go from q 1 at time to q 2 at time t, and S t (q 2,q 1 ) is the classical action integral, S t t dtpt 2 /2mVqt, 2.1b and where we have for simplicity assumed a Cartesian Hamiltonian of F degrees of freedom with all masses scaled to a common value, i.e., Hp,qp 2 /2mVq. This can be written in equivalent IVR form F as follows: 2i e iĥt/ dp dq q t p e is t (q,p )/ q t q, where q t q t (q,p ) is the coordinate at time t that has evolved from initial conditions (q,p ). It was shown some time ago 38 that Eq. 2.1 can be converted into a SC approximation for matrix elements of the Boltzmann operator by the often used replacement ti. The result involves classical trajectories in the pure imaginary time variable it, but in light of Newton s equation of motion, m d2 Vqt qt 2 dt q which thus becomes m d2 Vq q 2 d q, 2.4a. 2.4b This can be viewed as classical motion in the real time variable on the inverted potential energy surface, V(q). By introducing the real momentumlike variable p, which is related to the actual momentum p which is imaginary by p ip, the equations of motion read d d qp /m, d d p Vq, 2.5 q which we note involve only real-valued quantities. The IVR for the Boltzmann operator is then obtained from Eq. 2.3, e dp dq q Ĥ 2 p F e S (q,p )/ q q, 2.6a where q (q,p ) is the value of q() at which evolves from initial conditions (q,p ) via the equations of motion Eq. 2.5, and the action integral S (q,p ) is given by S dp 2 /2mVq. B. The flux-side correlation function. 2.6b For a Cartesian Hamiltonian Eq. 2.2, the flux operator Fˆ, Fˆ i Ĥ,hsqˆ, becomes Fˆ sq sq q q q sq q sq. 2.7a 2.7b Therefore, if one uses the IVR of Eq. 2.6 for the two Boltzmann operators in Eq. 1.3 and makes use of the derivative relations of the imaginary time action, q 2 S q 2,q 1 p 2, q 1 S q 2,q 1 p 1, the symmetry of the Boltzmann operator, q 2 e Ĥ q 1 q 1 e Ĥ q 2 2.8a 2.8b 2.9 and a judicious integration by parts then the flux-side correlation function of Eq. 1.2 becomes C fs t2 dq dp dp p q q p F sq sq q p p exps /2 q,p /S /2 q,p /qĥtq, 2.1 where qq /2 (q,p ) and qq /2 (q,p ) are the values of q() at /2 that evolve via the imaginary time equations of motion Eq. 2.5 with initial conditions (q,p ) and (q,p ), respectively. Equation 2.1 thus requires that one generates two imaginary time trajectories or real time trajectories on the inverted potential energy surface, both
3 J. Chem. Phys., Vol. 117, No. 21, 1 December 22 Boltzmann operator representation 967 starting from the same initial position q, which lies on the dividing surface s(q ) because of the delta function in the integrand. The above SC-IVR approximation for the Boltzmann operators can now be combined with whatever treatment of the real time dynamics in Eq. 2.1 i.e., the matrix elements of ĥ(t)] that is desired. For example, if the forward backward IVR Refs. 21,22 is used for the real time dynamics, the matrix element qĥ(t)q is given by qĥtq dps 2ip s 1 dq dp 2 F C q,p e is (q,p )/ qp q p q q, 2.11a where qp q and qp q are coherent state wave functions, e.g., qp q F/4 exp 2 qq 2 ip qq /, 2.11b (p,q ) are the initial conditions for real time trajectories that evolve from time to t via the normal molecular Hamiltonian, have a momentum jump at time t, p t p t p s sq t q t, 2.11c and then evolve backward in time from t to ; S is the forward backward action S t dtlps sq t t dtl, 2.11d Lp 2 /2mV(q) being the usual Lagrangian; and C is the Herman Kluk prefactor that involves derivatives of the final variables (q,p ) with respect to the initial ones (q,p ). A more approximate but much simpler treatment of the real time dynamics is the linearization approximation, 9 which assumes that the forward ( t) and backward (t ) trajectories are infinitesimally close to one another and which is also equivalent to the classical Wigner model. Though this approximation is not able to describe interference effects between different classical trajectories 26 because the forward and backward trajectories are assumed to be infinitesimally close, it has been seen in earlier work 9,27 to be reasonably accurate for these flux correlation functions. In this approximation the matrix elements of ĥ(t) are given by qĥtq2 F dp e ip (qq)/ hsq t, 2.12 where q t q t (q,p ) is the coordinate at time t that results from a real time trajectory with initial conditions (q,p ) with q (qq)/2. When this is used in Eq. 2.1, the explicit expression for the flux-side correlation function within the linearized approximation for the real time dynamics and the SC-IVR for the Boltzmann operators is FIG. 1. A schematic depiction of the various trajectories involved in Eq. 2.13: two imaginary time trajectories originate at q, with initial momenta p and p, respectively, and are propagated for time /2, arriving at positions q and q. A real time trajectory is then involved with initial point q (qq)/2 and initial momentum p. C fs t2 dq dp dp 2F dp p p q q sq sq q p p hsq t exps /2 q,p /S /2 q,p / ip qq/, 2.13 where qq /2 (q,p ),qq /2 (q,p ),q (qq)/2, and q t q t (q,p ). Thus to evaluate the integrand above, two imaginary time trajectories are begun at position q on the dividing surface with initial momentum variables p and p, respectively, and integrated for time /2; their final values are q and q. A real time trajectory, with initial coordinate q (qq)/2 and initial momentum p, is integrated until it is clear whether it will emerge in the reactant or product region in the former case the integrand is zero. Figure 1 depicts this situation schematically. III. NUMERICAL TESTS In this section we present some numerical results to test the error involved in using the SC Boltzmann operator in the flux-side correlation function discussed above. The example treated is a standard one, the 1d Eckart barrier, for which the potential function is VqV sech 2 aq, 3.1 with the parameters chosen to correspond approximately to the HH 2 reaction: V.425eV,a1.36 a.u. and m 16 a.u. The imaginary frequency is b 2V a 2 /m a.u.1619 cm 1, and the two Eckart parameters are 2V / b 13.29, u b , respectively; the cross over temperature, corresponding to u2, is 37 K. One complicating aspect of the SC Boltzmann operator should be noted: at the lowest temperatures considered below, the monodromy matrix elements, q /2 /p, appearing in Eq. 2.1 are negative for some values of the integration variables. Since the square root of this quantity is involved, this leads to an imaginary contribution to the matrix element of the Boltzmann operator but one knows that
4 968 J. Chem. Phys., Vol. 117, No. 21, 1 December 22 Y. Zhao and W. H. Miller FIG. 2. The flux-side correlation function for the 1d Eckart potential, using the SC-IVR for the Boltzmann operator and the exact quantum matrix element of ĥ(t); upper figure is for T1 K, and the lower one for T2 K. these matrix elements are real-valued! The most consistent way to deal with this is simply to set the integrand to zero for such values, and this is what has been done in the calculations presented below. Appendix A gives a discussion of the stationary phase approximation for the case of real exponents which justifies this procedure. To provide an unambiguous test of the accuracy of the SC Boltzmann operator for this example, we first carried out calculations using the exact quantum matrix elements of ĥ(t); i.e., we solved the time-dependent Schrödinger equation in a discrete variable representation in order to obtain qĥ(t)q exactly. The SC-IVR for the Boltzmann operator is thus the only approximation in this case. Figure 2 shows C fs (t) for a high temperature T1 K, upper figure and a low temperature T2 K, lower figure compared to the exact quantum results i.e., using exact quantum results for the Boltzmann operator and ĥ(t)]; the corresponding values of the rate constant are given in Table I and shown in Fig. 3. The conclusion is that the SC approximation for the Boltzmann operator is quite good for a wide range of temperature. We also carried out calculations within the linearized approximation for the real time dynamics, i.e., using Eqs In the present 1d example it is possible to carry out the integral over p in Eq analytically as well as that over q by virtue of the delta function, to give the following explicit expression: C fs 2 2 dp dp q p q p p p 2m qq exp S /2,p /S /2,p / iqq2m V b V qq 2 /. 3.2 A Filinov filtering of the integrand was necessary to tame the oscillatory character of the above integral particularly for the lowest temperatures, and Appendix B describes this procedure which is a slight generalization of what we have used before 39 in more detail. The Filinov parameter was taken sufficiently large 2 to be essentially in the limit i.e., larger values gave no change in the results; in this limit the Filinov procedure yields the value of the original integral. The integrand in Eq. 3.2 is quite localized in the two-dimensional space (p,p ), so that only a modest number of these values are needed to calculate the integral here by grid methods. For systems with more degrees of freedom, one would expect that Monte Carlo methods would make the multidimensional integrals efficient. The results of this linearized approximation for the real time dynamics are also given in Table I and Fig. 3, and one sees that this approximation degrades the quality of the rate constant, more so the lower the temperature, e.g., at T2 K, where the tunneling correction factor is 2, the result of the linearization approximation is about a factor of 2 too small. This is consistent with our earlier results 8,9 using the linearization approximation in the tunneling regime where the Boltzmann operator was treated exactly, by Monte Carlo path integration. Thus use of the SC Boltzmann operator with the linearized approximation for the real time dynamics is essentially just as good as that obtained with the exact treatment of the Boltzmann operator. IV. CONCLUDING REMARKS It has been shown how a SC-IVR can be utilized for the Boltzmann operator in reactive flux correlation functions. This can be combined with whatever level of treatment of the real time dynamics one chooses. Tests on a simple but illustrative example suggests that little error is introduced by treating the Boltzmann operator semiclassically, even at low temperatures where tunneling corrections are quite large. ACKNOWLEDGMENTS This work has been supported by the Director, Office of Science, Office of Basic Energy Sciences, Chemical Sci- TABLE I. Comparison of tunneling correction factors 2e V k(t) calculated by different methods. T K Semi. F Quan. ĥ(t) a Semi. F Semi. ĥ(t) b Quantum mechanical c a Calculated with the SC-IVR for the Boltzmann operator e Ĥ and the exact quantum result for the projection operator ĥ(t) see Eq b Calculated with the SC-IVR for e Ĥ and the linearized semiclassical approximation for ĥ(t) see Eq c Calculated by the exact quantum mechanical method.
5 J. Chem. Phys., Vol. 117, No. 21, 1 December 22 Boltzmann operator representation I SPA 2 Ax k k 2 e Sx k S(x k )/, 969 A7 where x k are all the roots of Eq. A5 for which Eq. A6 is also true. That is, one does not include all the stationary points of S(x) in the sum in Eq. A7, only those for which S(x k ). Note that if one forgot and included the stationary roots in Eq. A7 for which S(x k ), these would give an imaginary contribution, which is nonsensical since the original integral, Eq. A4, is manifestly real. FIG. 3. An Arrhenius plot of the thermal rate constant for the 1d Eckart barrier. The solid curve is the exact result, and the dashed line the classical result. The triangles and squares are the results of present calculations that use the SC-IVR for the Boltzmann operator, and the exact quantum result for the projection operator ĥ(t) triangles, or the linearized semiclassical approximation for ĥ(t) squares. ences Division of the U. S. Department of Energy under Contract No. DE-AC3-76SF98, and by the National Science Foundation Grant No. CHE Y.Z. gratefully thanks N. Makri for stimulating discussions. APPENDIX A: STATIONARY PHASE WITH REAL EXPONENTS The standard stationary phase approximation for a one dimensional integral of the form is I 1 dxaxe is(x)/ I 1 SPA k Ax k 2i e Sx k is(x k )/, A1 A2 where x k are all the points of stationary phase, i.e., roots of the equation Sx k. A3 In extending this result to the case of real exponents, i.e., for an integral of the form I 2 dxaxe S(x)/ A4 one notes that the dominant contributions to the integral come from regions where S(x) has a local minimum i.e., where the integrand has a local maximum with Sx k Sx k. A5 A6 Expanding S(x) quadratically about each such point and carrying out the usual Gaussian integral gives the analog of Eq. A2, APPENDIX B: A GENERALIZED FILINOV TRANSFORMATION FOR THE SC-IVR A recent paper 39 described the effectiveness of a generalized Filinov filtering procedure for evaluating multidimensional integrals with an oscillatory integrand, i.e., integrals of the form, I dz Rze i(z). B1 The procedure is to insert into the integrand of Eq. B1 the following representation of unity: 1 F e 1/4 T 1 dz e (zz )T (zz ) e it (zz ) B2 expand the exponent (z) quadratically about z, and then perform the Gaussian integral over z, giving the following result: I dz Rz 2 2iz exp 1 4 T z T 2iz 1 z. e i(z ) B3 The generalization in the procedure, compared to the original version, is inclusion of the linear term in the exponent in Eq. B2. In the treatment above, it was assumed that the preexponential factor in the integrand of Eq. B1, R(z), is sufficiently slowly varying that it was replaced by R(z ) in obtaining Eq. B3, but this may not always be the case. For example, for the following integral: I dq q e Aq 2 /2 e i pq, B4 which is in fact characteristic of some of the integrals encountered in SC-IVR applications, Eq. B3 gives the incorrect result I. It is easy to remedy this defect, however, by expanding R(z) to first order about z in the treatment leading up to Eq. B3. The resulting generalization of Eq. B3 is easily found to be
6 961 J. Chem. Phys., Vol. 117, No. 21, 1 December 22 Y. Zhao and W. H. Miller I dz 2 2iz Rz i Rz T 2iz 1 z e i(z ) exp 1 4 T z T 2iz 1 z. B5 It is not hard to show that Eq. B5 gives the exact result for the test case in Eq. B4. We now specialize Eq. B5 to the type of IVR integrals considered in this paper: the complex function (z) is written in terms of its real and imaginary parts, zzi z, B6 and it is assumed that the Hessian of the real part, (z), is negligible. Equation B5 then becomes I dz 2 2z Rz i Rz T 2z 1 z i z expiz z 1 4 T z i z T 2z 1 z i z. B7 The Gaussian width matrix and vector are the parameters that control the accuracy and efficiency of the procedure. The original integral is recovered in the limits and.) This above procedure, Eq. B7, was applied to the linearized SC-IVR expression, Eq. 2.13, in the variable p p p. It is necessary to carry out the Filinov filtering in only this one variable in order to obtain a well-behaved integrand for Monte Carlo integration. The explicit result of carrying out this procedure for Eq is C fs dq dp dp dp sq with 1 sq hsq t RiR2 1 q i2 1 M qp q,p M qp q,p 2 2 expi i 2 1 i], B8 Rp p, R1, S /2 q,q S /2 q,q, B9a B9b B9c 1 2 p /2 M qp q,p p /2 M qp q,p, B9d 1 4 M ppq,p M qp q,p M pp q,p M qp q,p, p qq/, p 2 M qpq,p M qp q,p, where M qp represents a monodromy matrix. B9e B9f B9g 1 W. H. Miller, J. Chem. Phys. 53, D. J. Tannor and S. Garashchuk, Annu. Rev. Phys. Chem. 51, W. H. Miller, J. Phys. Chem. A 15, W. H. Miller, J. Chem. Phys. 61, W. H. Miller, S. D. Schwartz, and J. W. Tromp, J. Chem. Phys. 79, T. Yamamoto, J. Chem. Phys. 33, G. A. Voth, D. Chandler, and W. H. Miller, J. Phys. Chem. 93, X. Sun, H. B. Wang, and W. H. Miller, J. Chem. Phys. 19, H. Wang, X. Sun, and W. H. Miller, J. Chem. Phys. 18, E. J. Heller, J. Chem. Phys. 65, R. C. Brown and E. J. Heller, J. Chem. Phys. 75, H. W. Lee and M. O. Scully, J. Chem. Phys. 73, J. S. Cao and G. A. Voth, J. Chem. Phys. 14, R. E. Cline and P. G. Wolynes, J. Chem. Phys. 88, V. Khidekel, V. Chernyak, and S. Mukamel, in Femtochemistry: Ultrafast Chemical and Physical Processes in Molecular Systems, edited by M. Chergui World Scientific, Singapore, 1996, p J. Shao, J. L. Liao, and E. Pollak, J. Chem. Phys. 18, J. L. Liao and E. Pollak, J. Chem. Phys. 111, Y. Zheng and E. Pollak, J. Chem. Phys. 114, J. L. Liao and E. Pollak, J. Chem. Phys. 116, N. Makri and K. Thompson, Chem. Phys. Lett. 291, W. H. Miller, Faraday Discuss. 11, X. Sun and W. H. Miller, J. Chem. Phys. 11, H. Wang, M. Thoss, and W. H. Miller, J. Chem. Phys. 112, R. Gelabert, X. Giménez, M. Thoss, H. Wang, and W. H. Miller, J. Phys. Chem. A 14, H. Wang, M. Thoss, K. L. Sorge, R. Gelabert, X. Giménez, and W. H. Miller, J. Chem. Phys. 114, R. Gelabert, X. Giménez, M. Thoss, H. Wang, and W. H. Miller, J. Chem. Phys. 114, T. Yamamoto, H. Wang, and W. H. Miller, J. Chem. Phys. 116, P. R. Feynman, Rev. Mod. Phys. 2, E. Jezek and N. Makri, J. Phys. Chem. 15, N. Makri and W. H. Miller, J. Chem. Phys. 116, M. F. Herman and E. Kluk, Chem. Phys. 91, J. H. V. Vleck, Proc. Natl. Acad. Sci. U.S.A. 14, W. H. Miller, Mol. Phys. 1, V. S. Filinov, Nucl. Phys. B 271, N. Makri and W. H. Miller, Chem. Phys. Lett. 139, J. D. Doll, D. L. Freeman, and T. L. Beck, Adv. Chem. Phys. 78, E. Pollak and B. Eckhardt, Phys. Rev. E 58, W. H. Miller, J. Chem. Phys. 55, H. Wang, D. E. Manolopoulos, and W. H. Miller, J. Chem. Phys. 115,
Coherent state semiclassical initial value representation for the Boltzmann operator in thermal correlation functions
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 21 1 JUNE 2002 Coherent state semiclassical initial value representation for the Boltzmann operator in thermal correlation functions Nancy Makri Departments
More informationSemiclassical initial value representation for rotational degrees of freedom: The tunneling dynamics of HCl dimer
JOURNAL OF CHEMICAL PHYSICS VOLUME 18, NUMBER 1 1 JUNE 1998 Semiclassical initial value representation for rotational degrees of freedom: The tunneling dynamics of HCl dimer Xiong Sun and William H. Miller
More informationUsing the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions
THE JOURNAL OF CHEMICAL PHYSICS 125, 224104 2006 Using the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions Jian Liu and William H. Miller
More informationValidating and Implementing Modified Filinov Phase Filtration in Semiclassical Dynamics
AIP/123-QED Validating and Implementing Modified Filinov Phase Filtration in Semiclassical Dynamics Matthew S. Church, 1 Sergey V. Antipov, 2 1, a) and Nandini Ananth 1) Department of Chemistry and Chemical
More informationSubmitted to The Journal of Physical Chemistry
Submitted to The Journal of Physical Chemistry Semiclassical Nonadiabatic Surface Hopping Wave Function Expansion at Low Energies: Hops in the Forbidden Region Journal: Manuscript D: Manuscript Type: Date
More informationQuantum mechanical transition state theory and tunneling corrections
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 9 1 MARCH 1999 Quantum mechanical transition state theory and tunneling corrections Ward H. Thompson Department of Chemistry and Biochemistry, University of
More informationForward-Backward Semiclassical Dynamics with Linear Scaling
J. Phys. Chem. A 1999, 103, 9479-9486 9479 Forward-Backward Semiclassical Dynamics with Linear Scaling Jiushu Shao School of Chemical Sciences, UniVersity of Illinois, 601 S. Goodwin AVenue, Urbana, Illinois
More informationReal time path integrals using the Herman Kluk propagator
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 7 15 FEBRUARY 2002 Real time path integrals using the Herman Kluk propagator John C. Burant and Victor S. Batista Department of Chemistry, Yale University,
More informationAnalytic continuation of real-time correlation functions to obtain thermal rate constants for chemical reaction
9 July 1999 Chemical Physics Letters 37 1999 463 468 www.elsevier.nlrlocatercplett Analytic continuation of real-time correlation functions to otain thermal rate constants for chemical reaction Haoin Wang,
More informationOptical response functions with semiclassical dynamics
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 Optical response functions with semiclassical dynamics W. G. Noid, Gregory S. Ezra, and Roger F. Loring Department of Chemistry and Chemical
More informationSemiclassical Dynamics with Quantum Trajectories: Formulation and Comparison with the Semiclassical Initial Value Representation Propagator
University of South Carolina Scholar Commons Faculty Publications Chemistry and Biochemistry, Department of 2-8-2003 Semiclassical Dynamics with Quantum Trajectories: Formulation and Comparison with the
More informationClassical Approximation to Nonradiative Electronic Relaxation in Condensed Phase Systems
J. Phys. Chem. A 1999, 103, 9539-9544 9539 Classical Approximation to Nonradiative Electronic Relaxation in Condensed Phase Systems Eran Rabani, S. A. Egorov,, and B. J. Berne* Department of Chemistry,
More informationSupporting Informations
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2018 Supporting Informations Fabio Gabas, Giovanni Di Liberto, Riccardo Conte, and Michele Ceotto
More informationSemiclassical molecular dynamics simulations of intramolecular proton transfer in photoexcited 2-2 -hydroxyphenyl oxazole
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 21 1 DECEMBER 2000 Semiclassical molecular dynamics simulations of intramolecular proton transfer in photoexcited 2-2 -hydroxyphenyl oxazole Victor Guallar,
More information1 Imaginary Time Path Integral
1 Imaginary Time Path Integral For the so-called imaginary time path integral, the object of interest is exp( τh h There are two reasons for using imaginary time path integrals. One is that the application
More informationARTICLES. Mixed Quantum and Forward-Backward Semiclassical Dynamics. Ed Bukhman and Nancy Makri*
J. Phys. Chem. A 2009, 113, 7183 7188 7183 ARTICLES Mixed Quantum and Forward-Backward Semiclassical Dynamics Ed Bukhman and Nancy Makri* Department of Chemistry, UniVersity of Illinois, Urbana, Illinois
More informationSC-Corr Code Package: Theory and Reference Manual
SC-Corr Code Package: Theory and Reference Manual Matthew S. Church Ananth Group Cornell University msc336@cornell.edu Version. June 6, 07 This document serves as a reference manual for the SC-Corr Code
More informationSemi-Classical Theory for Non-separable Systems :
Semi-Classical Theory for Non-separable Systems : Construction of Good Action-Angle Variables for Reaction Rate Constants* BY WILLIAM H. MILLER-^ Department of Chemistry, and Materials and Molecular Research
More informationWavepacket Correlation Function Approach for Nonadiabatic Reactions: Quasi-Jahn-Teller Model
Wavepacket Correlation for Nonadiabatic Reactions Bull. Korean Chem. Soc. 04, Vol. 35, No. 4 06 http://dx.doi.org/0.50/bkcs.04.35.4.06 Wavepacket Correlation Function Approach for Nonadiabatic Reactions:
More informationA comparison between different semiclassical approximations for optical response functions in nonpolar liquid solutions
THE JOURNAL OF CHEMICAL PHYSICS 122, 064506 2005 A comparison between different semiclassical approximations for optical response functions in nonpolar liquid solutions Qiang Shi Eitan Geva Department
More informationA quantum symmetry preserving semiclassical method
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 19 15 NOVEMBER 22 A quantum symmetry preserving semiclassical method Dmitri Babikov, a) Robert B. Walker, and Russell T Pack Theoretical Chemistry and Molecular
More informationNonadiabatic photodissociation dynamics of ICN in the à continuum: A semiclassical initial value representation study
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 13 1 APRIL 2000 Nonadiabatic photodissociation dynamics of ICN in the à continuum: A semiclassical initial value representation study Eduardo A. Coronado,
More informationJOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 7 15 FEBRUARY 2004
JOURAL OF CHEMICAL HYSICS VOLUME 20, UMBER 7 5 FEBRUARY 2004 ath integral calculation of thermal rate constants within the quantum instanton approximation: Application to the H CH 4 \H 2 CH 3 hydrogen
More informationSemiclassical molecular dynamics simulations of ultrafast photodissociation dynamics associated with the Chappuis band of ozone
Semiclassical molecular dynamics simulations of ultrafast photodissociation dynamics associated with the Chappuis band of ozone Victor S. Batista and William H. Miller a) Department of Chemistry, University
More informationMatching-pursuit/split-operator-Fourier-transform computations of thermal correlation functions
THE JOURNAL OF CHEMICAL PHYSICS 122, 064102 2005 Matching-pursuit/split-operator-Fourier-transform computations of thermal correlation functions Xin Chen, Yinghua Wu, and Victor S. Batista Department of
More informationAyan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur
Ayan Chattopadhyay Mainak Mustafi 3 rd yr Undergraduates Integrated MSc Chemistry IIT Kharagpur Under the supervision of: Dr. Marcel Nooijen Associate Professor Department of Chemistry University of Waterloo
More informationEnergy Conserving Approximations to the Quantum Potential: Dynamics with Linearized Quantum Force
University of South Carolina Scholar Commons Faculty Publications Chemistry and Biochemistry, Department of 1-15-2004 Energy Conserving Approximations to the Quantum Potential: Dynamics with Linearized
More informationSpecial functions and quantum mechanics in phase space: Airy functions
PHYSICAL REVIEW A VOLUME 53, NUMBER 6 JUNE 996 Special functions and quantum mechanics in phase space: Airy functions Go. Torres-Vega, A. Zúñiga-Segundo, and J. D. Morales-Guzmán Departamento de Física,
More informationVibrationally Mediated Bond Selective Dissociative Chemisorption of HOD on Cu(111) Supporting Information
Submitted to Chem. Sci. 8/30/202 Vibrationally Mediated Bond Selective Dissociative Chemisorption of HOD on Cu() Supporting Information Bin Jiang,,2 Daiqian Xie,,a) and Hua Guo 2,a) Institute of Theoretical
More informationEran Rabani, S. A. Egorov, a) and B. J. Berne Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 15 15 OCTOBER 1998 A comparison of exact quantum mechanical and various semiclassical treatments for the vibronic absorption spectrum: The case of fast vibrational
More informationQuantum dynamics via a time propagator in Wigner s phase space
Quantum dynamics via a time propagator in Wigner s phase space Michael Gro nager and Niels Engholm Henriksen Technical University of Denmark, Chemistry Department B, DTU-207, DK-2800 Lyngby, Denmark Received
More informationTheory of selective excitation in stimulated Raman scattering
Theory of selective excitation in stimulated Raman scattering S. A. Malinovskaya, P. H. Bucksbaum, and P. R. Berman Michigan Center for Theoretical Physics, FOCUS Center, and Department of Physics, University
More informationPath Integrals. Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz Jena
Path Integrals Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz 1 07743 Jena 5. Auflage WS 2008/09 1. Auflage, SS 1991 (ETH-Zürich) I ask readers to report
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [University of Illinois] On: 22 August 28 Access details: Access Details: [subscription number 788772948] Publisher Taylor & Francis Informa Ltd Registered in England and
More informationSC Corr-Code Package: Theory and Reference Manual
SC Corr-Code Package: Theory and Reference Manual Matthew S. Church Ananth Group Cornell University msc336@cornell.edu Version. November, 08 This is the reference manual for the SC Corr-Code Package (FORTRAN)
More informationPath integration in relativistic quantum mechanics
Path integration in relativistic quantum mechanics Ian H. Redmount and Wai-Mo Suen McDonnell Center for the Space Sciences Washington University, Department of Physics St. Louis, Missouri 63130 4899 USA
More informationForward backward semiclassical and quantum trajectory methods for time correlation functions
PCCP Dynamic Article Links Cite this: Phys. Chem. Chem. Phys., 2011, 13, 14442 14452 www.rsc.org/pccp PERSPECTIVE Forward backward semiclassical and quantum trajectory methods for time correlation functions
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationMixed Semiclassical Initial Value Representation Time-Averaging Propagator for Spectroscopic Calculations. Abstract
Mixed Semiclassical Initial Value Representation Time-Averaging Propagator for Spectroscopic Calculations Max Buchholz Institut für Theoretische Physik, Technische Universität Dresden, 6 Dresden, Germany
More informationConductance fluctuations at the integer quantum Hall plateau transition
PHYSICAL REVIEW B VOLUME 55, NUMBER 3 15 JANUARY 1997-I Conductance fluctuations at the integer quantum Hall plateau transition Sora Cho Department of Physics, University of California, Santa Barbara,
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationSemiclassical molecular dynamics simulations of excited state double-proton transfer in 7-azaindole dimers
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 20 22 MAY 1999 Semiclassical molecular dynamics simulations of excited state double-proton transfer in 7-azaindole dimers Victor Guallar, Victor S. Batista,
More informationThermal rate constant calculation using flux flux autocorrelation functions: Application to Cl H 2 HCl H reaction
Thermal rate constant calculation using flux flux autocorrelation functions: Application to ClH 2 HClH reaction Haobin Wang, Ward H. Thompson, a) and William H. Miller Department of Chemistry, University
More informationSemiclassical spin coherent state method in the weak spin-orbit coupling limit
arxiv:nlin/847v1 [nlin.cd] 9 Aug Semiclassical spin coherent state method in the weak spin-orbit coupling limit Oleg Zaitsev Institut für Theoretische Physik, Universität Regensburg, D-934 Regensburg,
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationWe start with some important background material in classical and quantum mechanics.
Chapter Basics We start with some important background material in classical and quantum mechanics.. Classical mechanics Lagrangian mechanics Compared to Newtonian mechanics, Lagrangian mechanics has the
More informationA New Form of Path Integral for the Coherent States Representation and its Semiclassical Limit
Brazilian Journal of Physics, vol. 35, no. 1, March, 5 175 A New Form of Path Integral for the Coherent States Representation and its Semiclassical Limit L. C. dos Santos and M. A. M. de Aguiar Instituto
More informationForward-Backward Quantum Dynamics for Time Correlation Functions
806 J. Phys. Chem. A 2004, 108, 806-812 Forward-Backward Quantum Dynamics or Time Correlation Functions Nancy Makri Departments o Chemistry and Physics, UniVersity o Illinois, 601 South Goodwin AVenue,
More informationSEMICLASSICAL INITIAL VALUE REPRESENTATION FOR COMPLEX DYNAMICS
Fakultät für Mathematik und Naturwissenschaften Fachrichtung Physik, Institut für Theoretische Physik Professur für Theoretische Atom- und Molekülphysik SEMICLASSICAL INITIAL VALUE REPRESENTATION FOR COMPLEX
More informationProbing surface diffuseness of nucleus-nucleus potential with quasielastic scattering at deep sub-barrier energies
PHYSICAL REVIEW C 73, 034607 (2006) Probing surface diffuseness of nucleus-nucleus potential with quasielastic scattering at deep sub-barrier energies K. Washiyama, K. Hagino, and M. Dasgupta 2 Department
More informationDetailed Balance in Ehrenfest Mixed Quantum-Classical Dynamics
J. Chem. Theory Comput. 26, 2, 229-235 229 Detailed Balance in Ehrenfest Mixed Quantum-Classical Dynamics Priya V. Parandekar and John C. Tully* Department of Chemistry, Yale UniVersity, P.O. Box 2817,
More informationSummary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:
The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:
More informationPhysics 582, Problem Set 4 Solutions
Physics 582, Problem Set 4 Solutions TAs: Hart Goldman and Ramanjit Sohal Fall 2018 1. PATH INTEGRAL FOR A PARTICLE IN A DOUBLE POTENTIAL WELL 1. In real time, q 0, T/2 q 0, T/2 = q 0 e i HT q 0 = where
More informationFirst Principles Semiclassical Calculations of Vibrational Eigenfunctions
First Principles Semiclassical Calculations of Vibrational Eigenfunctions The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation
More informationPath integral in quantum mechanics based on S-6 Consider nonrelativistic quantum mechanics of one particle in one dimension with the hamiltonian:
Path integral in quantum mechanics based on S-6 Consider nonrelativistic quantum mechanics of one particle in one dimension with the hamiltonian: let s look at one piece first: P and Q obey: Probability
More informationarxiv: v1 [quant-ph] 15 Dec 2011
Sharp and Infinite Boundaries in the Path Integral Formalism Phillip Dluhy and Asim Gangopadhyaya Loyola University Chicago, Department of Physics, Chicago, IL 666 Abstract arxiv:.3674v [quant-ph 5 Dec
More informationPartially-time-ordered Schwinger-Keldysh loop expansion of coherent nonlinear optical susceptibilities
PHYSICAL REVIEW A 77, 02380 2008 Partially-time-ordered Schwinger-Keldysh loop expansion of coherent nonlinear optical susceptibilities Shaul Mukamel Department of Chemistry, University of California,
More informationA Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space
Proceedings of Institute of Mathematics of NAS of Ukraine 00, Vol. 43, Part, 645 65 A Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space S. NASIRI and H. SAFARI Institute for Advanced
More informationPath integrals in quantum mechanics
Path integrals in quantum mechanics Phys V3500/G8099 handout #1 References: there s a nice discussion of this material in the first chapter of L.S. Schulman, Techniques and applications of path integration.
More informationOn the Adequacy of Mixed Quantum-Classical Dynamics in Condensed Phase Systems
10978 J. Phys. Chem. B 1999, 103, 10978-10991 FEATURE ARTICLE On the Adequacy of Mixed Quantum-Classical Dynamics in Condensed Phase Systems S. A. Egorov Theoretical Chemistry Institute and Department
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationInformation guided noise reduction for Monte Carlo integration of oscillatory functions
Chemical Physics Letters 4 (24) 446 452 www.elsevier.com/locate/cplett Information guided noise reduction for Monte Carlo integration of oscillatory functions Nancy Makri * Departments of Chemistry Physics,
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationFemtosecond photoelectron spectroscopy of the I 2 anion: A semiclassical molecular dynamics simulation method
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 8 22 FEBRUARY 1999 Femtosecond photoelectron spectroscopy of the I 2 anion: A semiclassical molecular dynamics simulation method Victor S. Batista, Martin
More informationUniform regularized semiclassical propagator for the x 2 potential
PHYSICAL REVIEW A VOLUME 53, NUMBER FEBRUARY 1996 Uniform regularized semiclassical propagator for the x potential Robert S. Manning * Center for Applied Mathematics, Cornell University, Ithaca, New York
More informationQuantum Dissipation: A Primer
Quantum Dissipation: A Primer P. Hänggi Institut für Physik Universität Augsburg NOISE-INDUCED ESCAPE Reaction-rate theory: fifty years after Kramers CONTENTS Peter Hanggi Lehrstuhl fur Theoretische
More informationWave Packet Representation of Semiclassical Time Evolution in Quantum Mechanics
Wave Packet Representation of Semiclassical Time Evolution in Quantum Mechanics Panos Karageorge, George Makrakis ACMAC Dept. of Applied Mathematics, University of Crete Semiclassical & Multiscale Aspects
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationThe Euclidean Propagator in Quantum Models with Non-Equivalent Instantons
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 2, 609 615 The Euclidean Propagator in Quantum Models with Non-Equivalent Instantons Javier CASAHORRAN Departamento de Física
More informationQuantum Corrections for Monte Carlo Simulation
Quantum Corrections for Monte Carlo Simulation Brian Winstead and Umberto Ravaioli Beckman Institute University of Illinois at Urbana-Champaign Outline Quantum corrections for quantization effects Effective
More informationCoulomb entangler and entanglement-testing network for waveguide qubits
PHYSICAL REVIEW A 72, 032330 2005 Coulomb entangler and entanglement-testing network for waveguide qubits Linda E. Reichl and Michael G. Snyder Center for Studies in Statistical Mechanics and Complex Systems,
More informationThe Path Integral Formulation of Quantum Mechanics
The Path Integral Formulation of Quantum Mechanics Shekhar Suresh Chandra March 15, 2005 Prerequisite: Knowledge of the Lagrangian Formalism (which can be found on my site also). Thirty-one years ago,
More informationThe remarkable properties of superfluids, in particular the
Simulation of dynamical properties of normal and superfluid helium Akira Nakayama* and Nancy Makri* *Department of Chemistry, University of Illinois at Urbana Champaign, 601 South Goodwin Avenue, Urbana,
More informationQUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology http://www.tut.fi/semiphys CONTENTS MOTIVATION PATH INTEGRAL APPROACH
More informationUnder evolution for a small time δt the area A(t) = q p evolves into an area
Physics 106a, Caltech 6 November, 2018 Lecture 11: Hamiltonian Mechanics II Towards statistical mechanics Phase space volumes are conserved by Hamiltonian dynamics We can use many nearby initial conditions
More informationS.K. Saikin May 22, Lecture 13
S.K. Saikin May, 007 13 Decoherence I Lecture 13 A physical qubit is never isolated from its environment completely. As a trivial example, as in the case of a solid state qubit implementation, the physical
More informationarxiv:cond-mat/ v4 [cond-mat.stat-mech] 14 Jan 2007
Exact analytical evaluation of time dependent transmission coefficient from the method of reactive flux for an inverted parabolic barrier arxiv:cond-mat/79v4 [cond-mat.stat-mech] 4 Jan 7 Raarshi Chakrabarti
More informationPreface. Preface to the Third Edition. Preface to the Second Edition. Preface to the First Edition. 1 Introduction 1
xi Contents Preface Preface to the Third Edition Preface to the Second Edition Preface to the First Edition v vii viii ix 1 Introduction 1 I GENERAL THEORY OF OPEN QUANTUM SYSTEMS 5 Diverse limited approaches:
More information(Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) p j . (5.1) !q j. " d dt = 0 (5.2) !p j . (5.
Chapter 5. Hamiltonian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 5.1 The Canonical Equations of Motion As we saw in section 4.7.4, the generalized
More informationSymmetry and Simplicity
Mobolaji Williams Motifs in Physics April 11, 2017 Symmetry and Simplicity These notes 1 are part of a series concerning Motifs in Physics in which we highlight recurrent concepts, techniques, and ways
More informationPath Integral for Spin
Path Integral for Spin Altland-Simons have a good discussion in 3.3 Applications of the Feynman Path Integral to the quantization of spin, which is defined by the commutation relations [Ŝj, Ŝk = iɛ jk
More informationB = 0. E = 1 c. E = 4πρ
Photons In this section, we will treat the electromagnetic field quantum mechanically. We start by recording the Maxwell equations. As usual, we expect these equations to hold both classically and quantum
More informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
More informationarxiv: v1 [cond-mat.stat-mech] 27 Jul 2011
arxiv:1107.5361v1 [cond-mat.stat-mech] 7 Jul 011 Exact quantum statistics for electronically nonadiabatic systems using continuous path variables Nandini Ananth 1 and Thomas F. Miller III 1 Division of
More information9 Instantons in Quantum Mechanics
9 Instantons in Quantum Mechanics 9.1 General discussion It has already been briefly mentioned that in quantum mechanics certain aspects of a problem can be overlooed in a perturbative treatment. One example
More informationEulerian semiclassical computational methods in quantum dynamics. Shi Jin Department of Mathematics University of Wisconsin-Madison
Eulerian semiclassical computational methods in quantum dynamics Shi Jin Department of Mathematics University of Wisconsin-Madison outline Classical mechanics: Hamiltonian system, discontinuous Hamiltonian,
More informationSpectral Degree of Coherence of a Random Three- Dimensional Electromagnetic Field
University of Miami Scholarly Repository Physics Articles and Papers Physics 1-1-004 Spectral Degree of Coherence of a Random Three- Dimensional Electromagnetic Field Olga Korotkova University of Miami,
More informationAlthough we rarely see a shock wave in everyday life,
Shear Viscosity under Shock-Loading Conditions Zhuohan Liang 1 and Steven M. Valone 1 T-4, Los Alamos Summer School MST-8, Los Alamos National Laboratory Although we rarely see a shock wave in everyday
More informationSimulation of quantum processes using entangled trajectory molecular dynamics
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 10 8 SEPTEMBER 2003 Simulation of quantum processes using entangled trajectory molecular dynamics Arnaldo Donoso Laboratorio de Física Estadística de Sistemas
More informationFragmentation production of Ω ccc and Ω bbb baryons
Physics Letters B 559 003 39 44 www.elsevier.com/locate/npe Fragmentation production of Ω ccc and Ω bbb baryons M.A. Gomshi Nobary a,b a Department of Physics, Faculty of Science, Razi University, Kermanshah,
More informationSemi-Classical Dynamics Using Hagedorn Wavepackets
Semi-Classical Dynamics Using Hagedorn Wavepackets Leila Taghizadeh June 6, 2013 Leila Taghizadeh () Semi-Classical Dynamics Using Hagedorn Wavepackets June 6, 2013 1 / 56 Outline 1 The Schrödinger Equation
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture Molecular Dynamics Simulation David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu MD of hard disks intuitive Review and Preview collision
More informationFeynman Path Integrals in Quantum Mechanics
Feynman Path Integrals in Quantum Mechanics Christian Egli October, 2004 Abstract This text is written as a report to the seminar course in theoretical physics at KTH, Stockholm. The idea of this work
More informationChapter 8. Chemical Dynamics
Chapter 8. Chemical Dynamics Chemical dynamics is a field in which scientists study the rates and mechanisms of chemical reactions. It also involves the study of how energy is transferred among molecules
More informationLikewise, any operator, including the most generic Hamiltonian, can be written in this basis as H11 H
Finite Dimensional systems/ilbert space Finite dimensional systems form an important sub-class of degrees of freedom in the physical world To begin with, they describe angular momenta with fixed modulus
More information