Non-adiabatic quantum molecular dynamics: basic formalism and case study

A general formalism is presented that treats selfconsistently and simultaneously classical atomic motion and quantum electronic excitations in dynamical processes of atomic many-body systems (non-adiabatic quantum molecular dynamics). On the basis of time-dependent density functional theory, coupled highly non-linear equations of motion are derived for arbitrary basis sets for the time-dependent Kohn-Sham orbitals. Possible approximations to make the approach practical for large atomic cluster systems are discussed. As a first application of the still exact equations of motion, non-adiabatic effects in the scattering of H⁺+H, as a case study, are investigated.     

Discrete variable representations and sudden models in quantum scattering theory

An exact formalism in which the scattering problem may be described by sets of coupled equations labeled either by basis functions or quadrature points is presented. Use of each frame and the simply evaluated unitary transformation which connects them results in an efficient procedure for performing quantum scattering calculations. Two approximations are compared with the IOS.     

Multidimensional reactive scattering with quantum trajectories: dynamics with 50-200 vibrational modes

The dynamics of ensembles containing thousands of quantum trajectories are studied for multidimensional systems undergoing reactive scattering. The Hamiltonian and equations of motion are formulated in curvilinear reaction path coordinates, for the case of a planar (zero-torsion) reaction path. In order to enhance the computational efficiency, an improved least squares fitting procedure is introduced. This scheme involves contracted basis sets and the use of inner and outer stencils around points where fitting is performed. This method is applied to reactive systems with 50-200 harmonic vibrational modes which are coupled to motion along the reaction coordinate. Dynamical results, including trajectory evolution and time-dependent reaction probabilities, are presented and power law scaling of computation time with the number of vibrational modes is described.     

The approximate conservation of P-helicity in rotational excitation: A new decoupling scheme

A new decoupling scheme for atom — rigid-rotor scattering based on Jacob and Wick's helicity formalism is suggested and found to be satisfactory for a wide range of systems of chemical interest. The number of coupled equations which needs to be solved simultaneously is greatly reduced. Conditions for the validity of this approximation are explored. Very accurate results are obtained for the opacity functions of the individual transitions.     

A simple model for the scattering of electrons by a rotating dipole

The coupled equations describing the interaction of one electron with a dipole and hard sphere are shown to be exactly soluble, even when the energy levels of the dipole are taken into account. This model is used to discuss the critical moment for binding the electron in the dipole field. The condition for the existence of Feshbach resonances is similarly discussed. When the model is applied to calculate scattering phase shifts, shape resonances are found.     

Applications of symmetry in coupled channel atom-surface scattering calculations

It has been shown that a partial decoupling of the coupled channel equations for atom-surface scattering can be obtained by the use of symmetry and this greatly reduces the numerical èffort required to solve the equations. The method is illustrated for the systems Ne/W(110), Ne/LiF(001) in the simplified framework of a recently-proposed sudden approximation.     

Coupled-channel calculations and the accuracy of the sudden approximation for atom—surface scattering

Exact coupled-channel calculations are presented for the scattering of Ne from W(110) and He from LiF(001), using symmetry to partly decouple the scattering equations. The results are used to test the recently proposed sudden approximation. For Ne/W(110), typical of all metals, the sudden approximation gives excellent quantitative accuracy. For the very unfavorable system He/LiF(001) good semiquantitative agreement is found with the exact results. It is concluded that the sudden approximately provides an efficient and accurate tool for atom—surface scattering calculations.     

Time-dependent wave packet propagation using quantum hydrodynamics

A new approach for propagating time-dependent quantum wave packets is presented based on the direct numerical solution of the quantum hydrodynamic equations of motion associated with the de Broglie–Bohm formulation of quantum mechanics. A generalized iterative finite difference method (IFDM) is used to solve the resulting set of non-linear coupled equations. The IFDM is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a “smart” Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is generalized to treat higher-dimensional problems and anharmonic potentials. The method is applied to a one-dimensional Gaussian wave packet scattering from an Eckart barrier, a one-dimensional Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. The 2D scattering results represent the first successful application of an accurate direct numerical solution of the quantum hydrodynamic equations to an anharmonic potential energy surface.     

Extended wave packet dynamics; exact solution for collinear atom,diatomic molecule scattering

The semiclassical wave packet dynamics method of Heller is extended to provide a formally exact theory of quantum mechanical motion for multidimensional anharmonic systems by introducing a complete, orthonormal, time-dependent basis of generalized oscillator functions. The exact wavefunction is expressed in terms of this basis and the expansions are shown to develop according to linear, coupled first-order differential equations. Application to collinear inelastic atom-diatomic molecule scattering demonstrates the feasibility and convergence of the new method.     

Dimensionality control of coupled scattering equations using partitioning techniques

A method of variable reduction of the dimensionality of the coupled equations for inelastic scattering is presented, based upon a projection operator P with a restricted range of orbital angular momentum states. For rotational states in the range O?j ?j^* and total angular momentum large, the coupled equations have dimensionality (j^* + 1) ? N ?(j^* + 1)², where the value of N is controlled by the choice of P. This is in contrast to conventional partitioning techniques which utilize further restrictions on the important molecular rotational states. The equations for the P subspace and its complementary Q subspace are decoupled by an approximation on the equation of motion of Qψ_scat. Information about scattering into the Q subspace is retained, within this degree of approximation, and is reintroduced at the end of the computation with little additional labor. The theory is developed in terms of atom-rigid-rotor scattering, although addition of vibrational modes would not in any way interfere with the basic techniques used.     

A nondiagonal quasidegenerate fourth-order perturbation theory

A canonical quasidegenerate Rayleigh-Schrödinger perturbation theory, correct through fourth order in the energy, is explored for a block-diagonal unperturbed Hamiltonian. The theory is developed completely within a Lie Algebra in Hilbert space. Explicit equations forn-particle transition elements in terms of solutions of simultaneous linear equations are presented. A two-dimensional anisotropic anharmonic oscillator is used to provide numerical results. The perturbation theory is shown to be stable under small separation of model and complement spaces. An iterative variant of the fourth order perturbation theory is introduced; the iterative variant is related to the non-iterative one in much the same way as nondegenerate coupled cluster theories are related to nondegenerate perturbation theory. The quasidegenerate coupled cluster theory appears to be stable in the presence of multiple intruder states.     

Coupled‐perturbed Hartree–Fock theory for quasi–one‐dimensional periodic systems: Calculation of static and dynamic nonlinear optical properties of model systems

An alternative method to solve the coupled‐perturbed Hartree–Fock (CPHF) equations for infinite quasi–one‐dimensional systems is presented. The new procedure follows a proposal made by Langhoff, Epstein, and Karplus to obtain perturbed wavefunctions free from arbitrary phase factors in each order of perturbation. It is based on the intermediate orthonormalization of the perturbed wavefunctions (which is different from the usual one) and a corresponding selection of the Lagrangian multipliers. In this way it is possible to incorporate the orthonormalization conditions into the set of CPHF equations. Moreover, a new, advantageous procedure to determine the derivatives of the wavefunction with respect to the quasimomentum k is presented. We report calculations of the dipole moment, the polarizability α, and the first hyperpolarizability β for different polymers (poly‐HF, poly‐H₂O, trans‐polyacetylene, polyyne, and polycarbonitrile) for different frequencies. These results are extensively compared with oligomer calculations. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 251–268, 2003     

An application of the refined born approximation to the solution of coupled equations involved in electron-ion and electron-atom scattering problems

A new approach to the solution of coupled equations involved in electron-ion and electron-atom scattering problems is proposed. This method is a combination of iteration and variation procedures. The main advantage of this method is that exchange terms can be calculated in a direct and straightforward manner. The method is based on the Lippmann-Schwinger equation and does not require trial functions satisfying appropriate boundary conditions. Using the Volterra formulation one can find the solution on an interval determined by the range of the exchange potential and the long-range potential terms can be taken into account by a projection procedure giving the asymptotic value of the reactance matrix. The method is tested on the case of electron-hydrogen atom scattering in the 1s-2s and 1s-2s-2p approximation.We have adapted the method proposed originally by Rayski to obtain solutions of coupled equations involved in electron-ion and electron-atom scattering. As mentioned in section 1 the construction of the method secures an automatic fulfilment of the boundary conditions. It allows an easy calculation of the exchange potential as well as an estimation of the introduced approximation. It gives also a possibility of detecting any spurious convergence. Moreover, it is important that this formalism can be applied in the case of normalized as well as unnormalized initial integral equations. This fact is of special importance in the case of long-range interactions. When the method is used for unnormalized (Volterra) equations it allows application of a very convenient projection procedure for treating the long-range terms in the direct potential.Electron-hydrogen atom collisions are investigated as a numerical illustration of the method. In the 1s-2s approximation the normalized equations were solved, while in the 1s-2s-2p approximation the solution was obtained with the help of Volterra equations and the long-range terms of the direct potential were taken into account by the projection procedure. In both cases the calculations were performed in the first iteration step and the obtained solutions agree fairly well with the results obtained by a numerical integration. It is not clear which set of results is more accurate. The numbers of parameters needed to obtain these results was not too large (not more than twenty in each channel) and decreased with the increase of the values of angular momentum and energy. The calculations were performed without weight functions, although the use of an appropriate weight function can improve the effectiveness of the method. This effectiveness could also be improved through a more lucky choice of trial functions.     

Ultra-low energy elastic scattering in a system of three He atoms

Differential Faddeev equations in total angular momentum representation are used for the first time to investigate ultra-low energy elastic scattering of a helium atom on a helium dimer. Six potential models of interatomic interaction are investigated. The results improve and extend the Faddeev equations based results known in literature. The employed method can be applied to investigation of different elastic and inelastic processes in three- and four-atomic weakly bounded systems below three-body threshold.     

Precise potential harmonic approach to directly solving Schrodinger equations of helium-like three-body systems

A complete potential harmonic scheme is presented,including the linked coupled hyperradial ordi nary differential equations and the secular equation of eigencnergy It has been used to directly solve the Scchrodinger equations of helium-like three-body systems (nuclear charge Z=1-9),and very accurate ground state eigonenergies as well as low-lying singlet excited state ones have been obtained     

Numerical Modeling of Steady-State Ion Transfer in Electrochemical Systems with Allowance for Migration

Numerical methods that are used for modeling steady-state ion transfer in electrochemical systems and account for the diffusion, migration, and convection are analyzed. An economical method for calculating steady-state processes is designed. The method splits a set of coupled equations of ion transfer and the electroneutrality conditions. The splitting is done in each iteration step by successively calculating distributions of the electric field potential and the electrolyte component concentrations. Techniques for approximate solution of redetermined set of difference equations used for calculating electric fields are analyzed. On the basis of this analysis, a novel technique is put forth, which ensures a minimum deviation from electroneutrality. Results of computational experiments are presented.     

Transition strengths and first-order properties of excited states from local coupled cluster CC2 response theory with density fitting

A new ab initio method for calculating transition strengths and orbital-unrelaxed first-order properties of singlet ground and excited states of extended molecular systems is presented. It is based on coupled cluster response theory at the level of the CC2 model with local approximations introduced to the doubles-excitation part of the wave function. Density fitting is employed for the calculation of the electron repulsion integrals, so that--with the exception of doubles amplitudes--only three-indexed objects do occur in the formalism. The new method was tested by performing calculations for a set of various molecules and excited states and by comparing the results with corresponding canonical (nonlocal) calculations. It turned out that for calculating transition strengths and properties of excited states the ordinary Boughton-Pulay domains are insufficient in numerous cases. To circumvent this problem a new scheme for extending domains is proposed, which is based on the solution of the coupled perturbed localization and Hartree-Fock equations. When such extended domains are used, a satisfactory agreement between canonical and local results is achieved.     

Precise potential harmonic approach to directly solving Schrödinger equations of helium-like three-body systems

A complete potential harmonic scheme is presented, including the linked coupled hyperradial ordinary differential equations and the secular equation of eigenenergy. It has been used to directly solve the Schrödinger equations of helium-like three-body systems (nuclear chargeZ = 1–9). and very accurate ground state eigenenergies as well as low-lying singlet excited state ones have been obtained.     

An apparatus for measuring Thomson scattering from an inductively coupled plasma

An experimental arrangement is presented which has been proven capable of measuring Thomson scattering from an inductively coupled plasma (ICP). This system has been shown to reject stray light to a sufficient extent that useful scattering signals can be measured as near as 0.3 nm from the incident laser wavelength. In this paper, design considerations are discussed and a viable experimental system is described. Preliminary Thomson scattering results are presented.