Time propagation of constrained coupled Gaussian wave packets

The dynamics of quantum systems can be approximated by the time propagation of Gaussian wave packets. Applying a time dependent variational principle, the time evolution of the parameters of the coupled Gaussian wave packets can be calculated from a set of ordinary differential equations. Unfortunately, the set of equations is ill behaved in most practical applications, depending on the number of propagated Gaussian wave packets, and methods for regularization are needed. We present a general method for regularization based on applying adequate nonholonomic inequality constraints to the evolution of the parameters, keeping the equations of motion well behaved. The power of the method is demonstrated for a nonintegrable system with two degrees of freedom.     

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.     

Quantum hydrodynamics: capturing a reactive scattering resonance

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.     

A preliminary study on sea wave packet equations on slowly varying topography

There is a common hypothesis for the presently popular mild-slope equations that wave particle motion is irrotational. In this paper, an attempt is made to abandon the irrotational assumption and to set up new sea wave packet equations on slowly varying topography by use of the WKBJ method. To simplify the deduction, the two-dimensional shallow water equations are used to describe the sea wave particle motion in the very shallow nearshore area. The established equations can give some characteristics of wave propagation near shore.     

A preliminary study on sea wave packet equations on slowly varying topography

There is a common hypothesis for the presently popular mild-slope equations that wave particle motion is irrotational. In this paper, an attempt is made to abandon the irrotational assumption and to set up new sea wave packet equations on slowly varying topography by use of the WKBJ method. To simplify the deduction, the two-dimensional shallow water equations are used to describe the sea wave particle motion in the very shallow nearshore area. The established equations can give some characteristics of wave propagation near shore.     

Locally implicit solution of a reaction-diffusion system with stiff kinetics

The Tyson-Fife reaction-diffusion equations are solved numerically using a locally implicit approach. Since the variables evolve at very different time scales, the resulting system of equations is stiff. The reaction term is responsible for the stiffness and the time step is increased by using an implicit method. The diffusion operator is evaluated explicitly and the system of implicit nonlinear equations is decoupled. The method is particularly useful for parameter values in which the equations are very stiff, such as the values obtained directly from the experimental reaction rate constants. Previous efforts modified the parameters on the equations to avoid stiffness. The equations then become a simplified model of excitable media and, for those cases, the locally implicit method gives a faster although less accurate solution. Nevertheless, since the modified equations no longer represent a particular chemical system an accurate solution is not as important. The algorithm is applied to observe the transition from simple motion to compound motion of a spiral tip.     

First principle nonlinear quantum dynamics using a correlation-based von Neumann entropy

A new concept to describe the quantum dynamics in complex systems is suggested. It extends established schemes based on the Dirac-Frenkel variation principle, e.g., the multi-configurational time-dependent Hartree (MCTDH) approach. The concept is based on a correlation-based von Neumann entropy (CvN-entropy) definition measuring the complexity of the wavefunction. Equations of motion are derived using a CvN-entropy constraint in the variational principle and result in a generally applicable effective Hamiltonian. It consists of the standard Hamilton operator and an additional nonlinear operator which limits the complexity of the wavefunction. Effectively, this nonlinear operator absorbs complex structures which are emerging in the wavefunction and allows one to introduce non-norm conserving equations of motion. Important aspects of the new concept are outlined studying the wave packet propagation on the diabatic B(2) potential energy surfaces of NO(2). First, it is demonstrated that during standard wave packet propagation the CvN-entropy increases strongly with time roughly independent of the coordinate systems employed. Second, one finds that employing CvN-entropy constrained MCTDH propagation yields improved wave function accuracy on longer time scales while compromising on the short time accuracy. Third, the loss of the wavefunction's norm is directly related to the overlap with the exact wavefunction. This provides an error estimate available without knowing an exact reference.     

Solving for the time-dependent multi-configuration hartree fock equations based on sine-discrete variable representation

A parallel quantum electrons wave packet computer code has been developed to study laser-atom interaction in the nonperturbative regime with attosecond resolution. The motion equations of the multi-configuration time-dependent hartree fock (MCTDHF) based on a sine discrete variable representation were solved by using an adaptive stepsize Runge-Kutta integrator of eight orders. Some efficient algorithms and strategies to accelerate the calculation velocity are introduced and discussed in details. Some illustrated imaginary time propagation and real time propagation have been respectively done in the paper. Single ionization probabilities calculated by using this one dimension MCTDHF model underestimate the accurate results calculated by solving time-dependent Schrodinger equation directly.     

A WAVE-OPTICS EFFECT WHICH ENHANCES LIGHT ABSORPTION BY CHLOROPHYLL IN VIVO

Abstract Light is widely thought to move through cells and tissues as rays. Actually light is inherently wave-like; ray motion is a special simple case of wave motion. To determine how well rays simulate light propagation and gradients within cells, absorption and scattering by particles were calculated with two sets of equations: one based on wave propagation within the particle and the other based on ray propagation. The ray model is found to substantially describe absorption at all particle size levels; however, it does not correctly describe scattering. Furthermore, for biological cells and structures larger than λ, the wave equations predict absorption cross sections which are 8–10% larger than expected on the basis of ray optics and extraction data: i.e. absorption is increased by a peculiar wave optics effect. This effect must enhance the efficiency with which the plants collect light energy for photosynthesis. The effect should be accounted for in corrections for the sieve effect.     

Influence of particle-particle interactions and particles rotational motion in traveling wave dielectrophoresis

Traveling wave dielectrophoresis provides an interesting method for the controlled movement of microsized particles in suspended mixtures, and as such is a promising tool in microfluidic technology. In this case, the electrostatic force acting on the particles has two components: one due to the spatially varying magnitude of the electric field and the other due to the spatially varying phase. The actual movement of the particle is determined by the combined effect of these two forces and corresponding torques, the viscous drag exerted by the fluid on the particle, and the electrostatic and hydrodynamic particle-particle interactions. This paper presents the first numerical simulations of the motion of particles subjected to all previous forces and torques. Our technique is based on a finite-element scheme in which the particles are moved using a direct simulation scheme respecting the fundamental equations of motion for both the fluid and the solid particles. The fluid-particle motion is resolved by the method of distributed Lagrange multipliers and the electrostatic forces are computed using the point-dipole approximation. Our simulations show that the particle behavior strongly depends on the mismatch of the dielectric properties between the particles and the fluid, and that the particle-particle interaction force as well as particles rotation speeds play crucial roles in the various regimes.     

Multidimensional quantum trajectories: applications of the derivative propagation method

In a previous publication [J. Chem. Phys. 118, 9911 (2003)], the derivative propagation method (DPM) was introduced as a novel numerical scheme for solving the quantum hydrodynamic equations of motion (QHEM) and computing the time evolution of quantum mechanical wave packets. These equations are a set of coupled, nonlinear partial differential equations governing the time evolution of the real-valued functions C and S in the complex action, S=C(r,t) + iS(r,t)/Planck's over 2pi, where Psi(r,t)=exp(S). Past numerical solutions to the QHEM were obtained via ensemble trajectory propagation, where the required first- and second-order spatial derivatives were evaluated using fitting techniques such as moving least squares. In the DPM, however, equations of motion are developed for the derivatives themselves, and a truncated set of these are integrated along quantum trajectories concurrently with the original QHEM equations for C and S. Using the DPM quantum effects can be included at various orders of approximation; no spatial fitting is involved; there is no basis set expansion; and single, uncoupled quantum trajectories can be propagated (in parallel) rather than in correlated ensembles. In this study, the DPM is extended from previous one-dimensional (1D) results to calculate transmission probabilities for 2D and 3D wave packet evolution on coupled Eckart barrier/harmonic oscillator surfaces. In the 2D problem, the DPM results are compared to standard numerical integration of the time-dependent Schrodinger equation. Also in this study, the practicality of implementing the DPM for systems with many more degrees of freedom is discussed.     

Trajectory approach to dissipative quantum phase space dynamics: Application to barrier scattering

The Caldeira-Leggett master equation, expressed in Lindblad form, has been used in the numerical study of the effect of a thermal environment on the dynamics of the scattering of a wave packet from a repulsive Eckart barrier. The dynamics are studied in terms of phase space trajectories associated with the distribution function, W(q,p,t). The equations of motion for the trajectories include quantum terms that introduce nonlocality into the motion, which imply that an ensemble of correlated trajectories needs to be propagated. However, use of the derivative propagation method (DPM) allows each trajectory to be propagated individually. This is achieved by deriving equations of motion for the partial derivatives of W(q,p,t) that appear in the master equation. The effects of dissipation on the trajectories are studied and results are shown for the transmission probability. On short time scales, decoherence is demonstrated by a swelling of trajectories into momentum space. For a nondissipative system, a comparison is made of the DPM with the "exact" transmission probability calculated from a fixed grid calculation.     

A meshless approach for electromagnetic simulation of metallic carbon nanotubes

In this paper, a study on the electromagnetic behaviour of a single wall carbon nanotube model is described. The electrons available for conduction are treated as a thin cylindrical layer fluid and their motion is described by means of classical hydrodynamics equations in linearized form. These equations are solved in time domain using the Smoothed Particle Hydrodynamics method. The method suitably handled runs on GRID environment.     

Thermal Charging Optimization of a Wavy-Shaped Nano-Enhanced Thermal Storage Unit

A wavy shape was used to enhance the thermal heat transfer in a shell-tube latent heat thermal energy storage (LHTES) unit. The thermal storage unit was filled with CuO–coconut oil nano-enhanced phase change material (NePCM). The enthalpy-porosity approach was employed to model the phase change heat transfer in the presence of natural convection effects in the molten NePCM. The finite element method was applied to integrate the governing equations for fluid motion and phase change heat transfer. The impact of wave amplitude and wave number of the heated tube, as well as the volume concertation of nanoparticles on the full-charging time of the LHTES unit, was addressed. The Taguchi optimization method was used to find an optimum design of the LHTES unit. The results showed that an increase in the volume fraction of nanoparticles reduces the charging time. Moreover, the waviness of the tube resists the natural convection flow circulation in the phase change domain and could increase the charging time.     

Multidimensional reactive scattering with quantum trajectories: dynamics with Morse vibrational modes

The reactive scattering of a wave packet is studied by the quantum trajectory method for a model system with up to 25 Morse vibrational modes. The equations of motion are formulated in curvilinear reaction path coordinates with the restriction to a planar reaction path. Spatial derivatives are evaluated by the least squares method using contracted basis sets. Dynamical results, including trajectory evolution and time-dependent reaction probabilities, are presented and analyzed. For the case of one Morse vibrational mode, the results are in good agreement with those derived through direct numerical integration of the time-dependent Schrodinger equation.     

3-D transient electrophoretic motion of a spherical particle in a T-shaped rectangular microchannel

This paper considers the electrophoretic motion of a spherical particle in an aqueous electrolyte solution in a T-shaped rectangular microchannel, where the size of the channel is close to that of the particle. This is a complicated transient process where the electric field, the flow field, and the particle motion are coupled together. A theoretical model was developed to investigate the influences of the applied electric potentials, the zeta potentials of the channel and the particle, and the size of the particle on the particle motion. A direct numerical simulation method using the finite element method is employed. This method employs a generalized Galerkin finite element formulation that incorporates both equations of the fluid flow and equations of the particle motion into a single variational equation where the hydrodynamic interactions are eliminated. The ALE method is used to track the surface of the particle at each time step. The numerical results show that the electric field in the T-shaped microchannel is influenced by the presence of the particle, and that the particle motion is influenced by the applied electric potentials and the zeta potentials of the channel and the particle. The path of the particle motion is dominated by the local electric field and the ratio of the zeta potential of the channel to that of the particle. The particle's velocity is also dependent on its size in a small channel.     

Cage motions induced by electronic and vibrational excitations: Cl2 in Ar

Femtosecond dynamics of molecular vibrations as well as cage motions in the B<--X transition of Cl2 in solid Ar have been investigated. We observed molecular vibrational wave-packet motion in experimental pump-probe spectra and an additional oscillation with a 500 fs period which is assigned to the zone-boundary phonon of the Ar crystal. The cage motion is impulsively driven by the B<--X transition due to the expansion of the electronic cloud of the chromophore. To clarify the underlying mechanism, we performed simulations based on the diatomics-in-molecules method which takes into account the different shapes of the Cl2 electronic wave function in the B and X states as well as the anisotropic interaction with the matrix. The simulation results show that Ar atom motion in the (100) plane is initiated by the electronic transition and that only those Ar atoms oscillate coherently with an approximately 500 fs period which are essentially decoupled from the molecular vibration. Their phase and time evolution are in good agreement with the experimentally observed oscillation, supporting the assignment as a displacive excitation of coherent phonons.     

Calcul du rayonnement acoustique d'une couche de mélange à l'aide des équations d'Euler linéarisées

The noise radiated by a mixing layer is computed by an hybrid approach using linearized Euler's equations (LEE) as a wave operator. The acoustic source terms are obtained by substituting the aerodynamic fluctuating velocity field computed by the Large Eddy Simulation (LES) method into the right-hand side of the LEE. The radiated sound field as calculated in this way is in good agreement with the direct calculation obtained by LES. This study provides support to the validity of the acoustic analogy based on the LEE, both for the noise source terms and the mean flow effects on sound propagation.     

Molecular dynamics integration and molecular vibrational theory. I. New symplectic integrators

New symplectic integrators have been developed by combining molecular dynamics integration with the standard theory of molecular vibrations to solve the Hamiltonian equations of motion. The presented integrators analytically resolve the internal high-frequency molecular vibrations by introducing a translating and rotating internal coordinate system of a molecule and calculating normal modes of an isolated molecule only. The translation and rotation of a molecule are treated as vibrational motions with the vibrational frequency zero. All types of motion are thus described in terms of the normal coordinates. The method's time reversibility requirement was used to determine the equations of motion for internal coordinate system of a molecule. The calculation of long-range forces is performed numerically within the generalized second-order leap-frog scheme, in the same way as in standard second-order symplectic methods. The new methods for integrating classical equations of motion using normal mode analysis allow us to use a long integration step and are applicable to any system of molecules with one equilibrium configuration.