Unitary quantum time evolution by iterative Lanczos reduction

Park, Tae Jun; Light, J. C.
November 1986
Journal of Chemical Physics;11/15/1986, Vol. 85 Issue 10, p5870
Academic Journal
A general unitary time evolution method for wave packets defined on a fixed L2 basis is developed. It is based on the Lanczos reduction of the full N×N Hamiltonian to a p-dimensional subspace defined by the application of H p-1 times to the initial vector. Unitary time evolution in the subspace is determined by exp{-iHpt}, retaining accuracy for a time interval τ, which can be estimated from the Lanczos reduced Hamiltonian Hp. The process is then iterated for additional time intervals. Although accurate results over long times can be obtained, the process is most efficient for large systems over short times. Time evolution employing this method in one- (unbounded) and two-dimensional (bounded) potentials are done as examples using a distributed Gaussian basis. The one-dimensional application is to direct evaluation of a thermal rate constant for the one-dimensional Eckart barrier.


Related Articles

  • Quantum flux operators and thermal rate constant: Collinear H+H2. Park, Tae Jun; Light, J. C. // Journal of Chemical Physics;4/15/1988, Vol. 88 Issue 8, p4897 

    The exact quantum formulation of the thermal rate constant, k(T), given by Miller et al. [W. H. Miller, J. Chem. Phys. 61, 1823 (1974); W. H. Miller, S. D. Schwartz, and J. W. Tromp, ibid. 79, 4889 (1983)] is evaluated in a localized [ATOTHER]@sL[/ATOTHER] 2 basis (distributed Gaussian basis)...

  • Relativistic Gaussian basis sets for molecular calculations: Cs–Hg. Koga, Toshikatsu; Tatewaki, Hiroshi; Matsuoka, Osamu // Journal of Chemical Physics;10/22/2002, Vol. 117 Issue 16, p7813 

    Relativistic Gaussian basis sets suitable for molecular calculations are presented for the 26 atoms [sub 55]Cs through [sub 80]Hg. This is the second part of our systematic generation of the basis sets, proposed in a previous paper [J. Chem. Phys. 115, 3561 (2001)]. For the 26 atoms, the maximum...

  • Gaussian basis sets for CO2 molecule generated with the molecular improved generator coordinate Hartree-Fock method. Barreto, M.T.; Muniz, E.P.; Jorge, F.E.; Cunha, A.G. // Theoretical Chemistry Accounts: Theory, Computation, & Modeling;Mar2005, Vol. 113 Issue 2, p69 

    The molecular improved generator coordinate Hartree-Fock (MIGCHF) method is used to generate increasing size atom-centered Gaussian basis sets for the CO2 molecule. From these basis sets total HF energies and second-order correlation energies were calculated and compared with results obtained...

  • Gaussian basis sets for the atoms from K through Xe. Centoducatte, R; de Castro, E VR; Jorge, F E // Canadian Journal of Chemistry;Feb2001, Vol. 79 Issue 2, p121 

    An improved generator coordinate Hartree-Fock (IGCHF) method is used to generate Gaussian basis sets for the atoms from K (Z = 19) through Xe (Z = 54). The Griffin-Hill-Wheeler-HF equations are integrated using the integral discretization technique. The ground state HF total energies obtained by...

  • Wave packet evolution of damped oscillator with a time dependent linear potential. Ikot, Akpan N.; Uwah, Edet J.; Akpabio, Louis E.; Akpan, Ita O. // Latin-American Journal of Physics Education;May2010, Vol. 4 Issue 2, p331 

    We discuss the general solution of a time-dependent Schrödinger wave equation (SWE) with time-dependent linear potential within the framework of invariant operators. We obtain the Gaussian wave packet evolution by choosing the ansatz for the weight factor of the form g (λ) = Ae...

  • Bohmian dynamics on subspaces using linearized quantum force. Rassolov, Vitaly A.; Garashchuk, Sophya // Journal of Chemical Physics;4/15/2004, Vol. 120 Issue 15, p6815 

    In the de Broglie–Bohm formulation of quantum mechanics the time-dependent Schrödinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is...

  • Quasi-bound states, resonance tunnelling, and tunnelling times generated by twin symmetric barriers. Maheswari, A. Uma; Prema, P.; Mahadevan, S.; Shastry, C. S. // Pramana: Journal of Physics;Dec2009, Vol. 73 Issue 6, p969 

    In analogy with the definition of resonant or quasi-bound states used in three-dimensional quantal scattering, we define the quasi-bound states that occur in one-dimensional transmission generated by twin symmetric potential barriers and evaluate their energies and widths using two typical...

  • Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. VandeVondele, Joost; Hutter, Jürg // Journal of Chemical Physics;9/21/2007, Vol. 127 Issue 11, p114105 

    We present a library of Gaussian basis sets that has been specifically optimized to perform accurate molecular calculations based on density functional theory. It targets a wide range of chemical environments, including the gas phase, interfaces, and the condensed phase. These generally...

  • Nonequilibrium membrane fluctuations driven by active proteins. Lin, Lawrence C.-L.; Gov, Nir; Brown, Frank L. H. // Journal of Chemical Physics;2/21/2006, Vol. 124 Issue 7, p074903 

    We extend a model for nonthermal membrane undulations driven by active (adenosine triphosphate-dependent or light-harvesting) membrane proteins [N. Gov, Phys. Rev. Lett. 93, 268104 (2004)]. The present model accounts for the fact that proteins can diffuse laterally across the membrane surface...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics