Sech wave packets, their Wigner functions and Bohmian trajectories

We have investigated the dynamics of a quantum particle in the optical lattice potential. Initially, the quantum particle was represented by a Gaussian wave packet, located in the center of the well. The corresponding Schrödinger equation was solved explicitly by the method of the Chebyshev global propagation. Obtained solutions were also used for the construction of the Wigner functions. We found a great number of local abrupt changes of the solution shape. To explain this behavior, we used the fact that structurally stable systems, which form the largest class of the low dimensional dynamical systems, can be modeled and classified according to the catastrophe theory. All important features of the exact solution were explained on the basis of the mathematical properties of the catastrophic model. Such an approach enabled us to extract relevant information out of numerical solutions without employing any kind of approximations. We have investigated the influence of the Wigner catastrophes on the details of the quantum-classical correspondence breakdown. The wave packet was found to expand rapidly, filling the whole classically available area of the phase space. It was found that its self-interference pattern saturates quickly. A region of the phase space emerges in which the Wigner function oscillations transform into the singularity driven fluctuations. Once this region covers the whole area of the phase space, a wave packet dynamics enters into the new regime where its Wigner function fluctuates around the ergodic average. It will be shown that all mentioned processes are caused by the proliferation of the catastrophes and their mutual interactions.

Summary form only given. Wave packets (WPs) has recently emerged as a very active research field in quantum optics. We show here that the evolution of WPs may result in new class of large-scale interference effects: the formation of highly regular spatio-temporal patterns in the QM probability, |/spl Psi/|/sup 2/, with the key mechanism for this being pair interference between the eigen-modes of the system and the degeneracy of this interference. The new patterns have been discovered in the probability density |/spl Psi/|/sup 2/(x,t) of a particle in an infinitely deep square well, by numerical simulations of the evolution of a Gaussian WP; later work by a few groups used other initial conditions and media (e.g. Bose-Einstein condensate and angular WPs). In all these cases, the QM distribution shows a net of canals (minima), and ridges (maxima), running as straight lines in space-time (x,t). Similar patterns in the near field of a diffraction grating, and in matter waves in atom optics reveals a new aspects of the Talbot effect. So far the numerical exploration of these patterns has relied on the specific choice (i) of the initial conditions and (ii) of the square well; no physical explanation of the patterns has been suggested. Recently, we observed a quantum carpet under more general conditions, offered its universal explanation, and predicted distinct curved patterns in general case of an anharmonic potential.

PurposeThe paper explores the spectral solution of a partial differential equation governing nonlinear transmission lines (NLTLs). It aims to present an efficient simulation technique for a nonlinear fractional-order transmission line model and to show the effect of fractional derivatives on the evolution of solitons and wave packets of increasing spatial frequency.Design/methodology/approachThe paper is concerned with the solution of the equation governing NLTLs for both integer and fractional time derivatives. The spectral method shall be considered with Hermite and Malmquist-Takenaka orthogonal functions and integer and fractional implementations of the implicit trapezoidal rule. While the Hermite functions and Malmquist-Takenaka functions are suitable for solitons, the Malmquist-Takenaka functions are selected as these are superior for approximating wave packets of increasing frequency (Iserles, Luong and Webb, 2023). Furthermore, they possess properties that are favourable from a computational viewpoint.FindingsResults will show that the proposed approach is effective for simulating NLTLs with both integer and fractional time derivatives. In particular, fractional derivatives are shown to have a significant effect on the evolution of wave packets on a nonlinear transmission line.Originality/valueThe method is novel in using Malmquist-Takenaka basis functions for a spectral solution of a nonlinear transmission line equation with both integer and fractional time derivatives and for examining the effect of fractional derivatives on wave packets.

Effects of the medium fractionality and oscillating potential profiles on the Superarrivals of the Gaussian wave packets

We derive an analytical solution to the time-dependent Schrödinger equation for the transmission of a Gaussian wave packet through an arbitrary potential of finite range. We consider the situation where the initial Gaussian wave packet is sufficiently broad in momentum space to include the full resonance structure of the system in the dynamical description. We find that the transmitted wave packet may be written as the free evolving Gaussian wave packet solution times a transient term. We demonstrate that both at very large distances and very long times the transient term tends to the transmission amplitude of the system and hence the time-evolving solution reproduces the resonance spectra of the system. We also prove that at a fixed distance and very long times the analytical solution is t−3/2 which extends previous analysis on this issue to arbitrary finite-range potentials. Our results are exemplified for a multibarrier system.

Variational and wave packet perturbation theory methods are implemented to study the effects on photodissociation and scattering cross sections due to nonadiabatic coupling between two electronic potential energy surfaces. The simplest variational method we implemented can handle large coupling strengths between surfaces, but its restriction to one Gaussian wave packet per surface is found to be too severe in some instances. In such situations the wave packet perturbation theory, which coherently superimposes multiple Gaussians on each surface, overcomes the difficulties encountered by the variational method and is seen to give quite accurate photodissociation and Landau–Zener cross sections over a wide energy range, so long as the coupling between surfaces is reasonably small. In this paper we give details of the appropriate computational schemes to adopt for calculating photodissociation and scattering cross sections from propagating wave packets. In addition, we discuss the strengths and weaknesses of Gaussian wave packet methods for studying nonadiabatic effects.

We analyze the semiclassical evolution of Gaussian wave packets in chaotic systems. We show that after some short time a Gaussian wave packet becomes a primitive WKB state. From then on, the state can be propagated using the standard time-dependent WKB scheme. Complex trajectories are not necessary to account for the long-time propagation. The Wigner function of the evolving state develops the structure of a classical filament plus quantum oscillations, with phase and amplitude being determined by geometric properties of a classical manifold.

The exact solution of the Schrodinger equation is constructed for time-dependent Hamiltonian systems involving quadratic, inverse quadratic and (1/x)p+p(1/x) terms. The solution corresponds to a semiclassical Gaussian wave packet centered around the classical guiding trajectory. The behavior of the Gaussian wave packet we obtained is very similar to that of a classical wave packet because its guiding trajectory follows the classical equation of motion. Gaussian wave packets are ubiquitous in quantum mechanics and are fundamental states of many physical systems that exhibit various nonclassical properties.

Bohmian trajectories are considered for a particle that is free (i.e. the potential energy is zero), except for a half-line barrier. On the barrier, both Dirichlet and Neumann boundary conditions are considered. The half-line barrier yields one of the simplest cases of diffraction. Using the exact time-dependent propagator found by Schulman, the trajectories are computed numerically for different initial Gaussian wave packets. In particular, it is found that different boundary conditions may lead to qualitatively different sets of trajectories. In the Dirichlet case, the particles tend to be more strongly repelled. The case of an incoming plane wave is also considered. The corresponding Bohmian trajectories are compared with the trajectories of an oil drop hopping on the surface of a vibrating bath.

In Bohm's causal or trajectory interpretation of quantum mechanics, it is straightforward to determine, from a sufficiently large number of calculated particle trajectories, probability distributions for time delays caused by a potential barrier, or for transmission times through a barrier. We show that these distributions can be calculated directly and more efficiently from probability currents, without the calculation of Bohm trajectories as an intermediate step. The ideas are illustrated for Gaussian wave packets incident on a square potential barrier and used to explain why average causal delay times differ from the average delay times calculated in other approaches.

A high-dimensional quantum wave packet approach based on Gaussian wave packets in Cartesian coordinates is presented. In this method, the high-dimensional wave packet is expressed as a product of time-dependent complex Gaussian functions, which describe the motion of individual atoms. It is applied to the ultrafast geometrical rearrangement dynamics of NO doped cryogenic Ne matrices after femtosecond laser pulse excitation. The static deformation of the solid due to the impurity as well as the dynamical response after femtosecond excitation are analyzed and compared to reduced dimensionality studies. The advantages and limitations of this method are analyzed in the perspective of future applications to other quantum solids.

We investigate the gravitational lensing effects on neutrino oscillations within the framework of Modified Gravity (MOG). By modeling neutrinos as both plane waves and Gaussian wave packets, we analyze how the presence of a massive lensing object between the neutrino source and detector modifies the flavor transition probabilities. In the plane-wave approximation, we derive analytical expressions for lensing-induced interference patterns, showing explicit dependence on the MOG parameter α , the impact parameters, and the neutrino mass hierarchy. Extending the analysis to a wave-packet framework, we account for decoherence effects and quantify the damping of oscillation amplitudes due to finite energy and spatial spreads of the wave packets. We compute the decoherence length and damping factor for a realistic Sun–Earth system and demonstrate how both the absolute neutrino mass scale and the MOG parameter influence coherence loss. Our results indicate that gravitational lensing in alternative gravity theories can leave observable imprints on neutrino oscillation patterns, offering a novel probe of both neutrino properties and modifications of General Relativity.

We find a relationship between the dynamics of the Gaussian wave packet and the dynamics of the corresponding Gaussian Wigner function from the Hamiltonian/symplectic point of view. The main result states that the momentum map corresponding to the natural action of the symplectic group on the Siegel upper half space yields the covariance matrix of the corresponding Gaussian Wigner function. This fact, combined with Kostant’s coadjoint orbit covering theorem, establishes a symplectic/Poisson-geometric connection between the two dynamics. The Hamiltonian formulation naturally gives rise to corrections to the potential terms in the dynamics of both the wave packet and the Wigner function, thereby resulting in slightly different sets of equations from the conventional classical ones. We numerically investigate the effect of the correction term and demonstrate that it improves the accuracy of the dynamics as an approximation to the dynamics of expectation values of observables.

We present results of real-time tracking of atomic two-electron dynamics in an autoionizing transient wave packet in krypton. A coherent superposition of two Fano resonances is excited with a femtosecond extreme-ultraviolet pulse. The evolution of the corresponding wave packet is subsequently probed with a delayed infrared pulse. In our specific case, we get access to the interference between one- and two-electron excitation channels in the launched wave packet, which is superimposed on its decay through autoionization. A simple model is able to account for the observed dynamical evolution of this wave packet.

The complex quantum Hamilton-Jacobi equation-Bohmian trajectories (CQHJE-BT) method is introduced as a synthetic trajectory method for integrating the complex quantum Hamilton-Jacobi equation for the complex action function by propagating an ensemble of real-valued correlated Bohmian trajectories. Substituting the wave function expressed in exponential form in terms of the complex action into the time-dependent Schrödinger equation yields the complex quantum Hamilton-Jacobi equation. We transform this equation into the arbitrary Lagrangian-Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation describing the rate of change in the complex action transported along Bohmian trajectories is simultaneously integrated with the guidance equation for Bohmian trajectories, and the time-dependent wave function is readily synthesized. The spatial derivatives of the complex action required for the integration scheme are obtained by solving one moving least squares matrix equation. In addition, the method is applied to the photodissociation of NOCl. The photodissociation dynamics of NOCl can be accurately described by propagating a small ensemble of trajectories. This study demonstrates that the CQHJE-BT method combines the considerable advantages of both the real and the complex quantum trajectory methods previously developed for wave packet dynamics.