High resolution finite volume scheme for the quantum hydrodynamic equations

The theory of quantum fluid dynamics (QFD) helps nanotechnology engineers to understand the physical effect of quantum forces. Although the governing equations of quantum fluid dynamics and classical fluid mechanics have the same form, there are two numerical simulation problems must be solved in QFD. The first is that the quantum potential term becomes singular and causes a divergence in the numerical simulation when the probability density is very small and close to zero. The second is that the unitarity in the time evolution of the quantum wave packet is significant. Accurate numerical evaluations are critical to the simulations of the flow fields that are generated by various quantum fluid systems.A finite volume scheme is developed herein to solve the quantum hydrodynamic equations of motion, which significantly improve the accuracy and stability of this method. The QFD equation is numerically implemented within the Eulerian method. A third-order modified Osher–Chakravarthy (MOC) upwind-centered finite volume scheme was constructed for conservation law to evaluate the convective terms, and a second-order central finite volume scheme was used to map the quantum potential field. An explicit Runge–Kutta method is used to perform the time integration to achieve fast convergence of the proposed scheme.In order to meet the numerical result can conform to the physical phenomenon and avoid numerical divergence happening due to extremely low probability density, the minimum value setting of probability density must exceed zero and smaller than certain value. The optimal value was found in the proposed numerical approach to maintain a converging numerical simulation when the minimum probability density is 10^?5 to 10^?12. The normalization of the wave packet remains close to unity through a long numerical simulation and the deviations from 1.0 is about 10^?4.To check the QFD finite difference numerical computations, one- and two-dimensional particle motions were solved for an Eckart barrier and a downhill ramp barrier, respectively. The results were compared to the solution of the Schrödinger equation, using the same potentials, which was obtained using by a finite difference method. Finally, the new approach was applied to simulate a quantum nanojet system and offer more intact theory in quantum computational fluid dynamics.     

Quantum tomography as a new approach to simulating quantum processes

A new method is proposed for ab initio calculations of nonstationary quantum processes on the basis of a probability representation of quantum mechanics with the help of a positive definite function (quantum tomogram). The essence of the method is that an ensemble of trajectories associated with the characteristics of the evolution equation for the quantum tomogram is considered in the space where the quantum tomogram is defined. The method is applied for detailed analysis of transient tunneling of a wave packet. The results are in good agreement with the exact numerical solution to the Schrödinger equation for this system. The probability density distributions are obtained in the coordinate and momentum spaces at consecutive instances. For transient tunneling of a wave packet, the probability of penetration behind the barrier and the time of tunneling are calculated as functions of the initial energy.     

On a Liouville-master equation formulation of open quantum systems

The dynamics of open quantum systems is formulated in terms of a probability distribution on the underlying Hilbert space. Defining the time-evolution of this probability distribution by means of a Liouvillemaster equation the time-dependent wave function of the system becomes a stochastic Markov process in the sense of classical probability theory. It is shown that the equation of motion for the two-point correlation function of the random wave function yields the quantum master equation for the statistical operator. Stochastic simulations of the Liouville-master equation are performed for a simple example from quantum optics and are shown to be in perfect agreement with the analytical solution of the corresponding equation for the statistical operator.     

Nonlinear diffraction in a quantum-dot system with allowance for the Hubbard interaction

The propagation of monochromatic laser radiation in a volume system of quantum dots (QDs) that are tunnel-coupled along one axis is considered. The electron energy spectrum of the QD system is modeled in the tight-binding approximation with allowance for the Coulomb interaction of electrons in the Hubbard model. The electromagnetic field of laser radiation in a QD system is described quasi-classically by Maxwell equations; as applied to this problem, they are reduced to a non-one-dimensional wave equation for the vector potential. As a result of the analysis of the wave equation in the approximation of varying amplitudes and phases, an effective equation describing the electromagnetic field in a QD system is obtained and numerically solved. The influence of the parameters of the system and the amplitude and frequency of the field of incident laser radiation on the character of its propagation is investigated. Nonmonotonic dependences of the factor characterizing the laser beam diffraction spread on the parameters of the electron energy spectrum of the system are obtained.     

An equation for the quantum dissipative system in statistical mechanics

A formalism for describing quantum dissipative systems in statistical mechanics is developed. A new equation of the Lindblad type with a quadratic superoperator consisting of Hermitian dissipative operators is derived from the Bloch equation for temperature density matrix using the Feynman integral over the trajectories with a modified Menskii weight functional. By way of example, this equation is solved for a one-dimensional quantum harmonic oscillator with linear dissipation. Applying the projection operator technique, an integral-differential equation for a reduced temperature statistical operator is obtained, which is analogous to the Zwanzig equation in statistical mechanics, and its formal solution is found as a convergent series. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 12, pp. 30–34, December, 2006.     

Investigation of self-focusing a Gaussian laser pulse in a weakly relativistic and ponderomotive regime inside a collisionless warm quantum plasma

This article investigates nonlinear self-focusing of an intense right hand circularly polarized Gaussian profile laser pulse in a weakly relativistic and ponderomotive regime inside a collisionless and unmagnetized warm quantum plasma. The nonlinear propagation equation for laser pulse in plasma has been derived. Then, the evolution differential equation for laser spot-size was obtained with considering the parabolic equation approach under the Wentzel-Kramers-Brillouin and paraxial ray approximations. This differential equation was solved numerically by fourth-order Runge-Kutta method. It is shown that our solution confirms the results of the self-focusing of the laser pulse in a weakly relativistic ponderomotive regime in cold quantum plasma in extreme conditions. Numerical results indicate that self-focusing of the laser pulse in the presence of relativistic and ponderomotive nonlinearity inside warm quantum plasma is improved in comparison with relativistic and ponderomotive cold quantum plasma.     

On the description of noise in quantum systems

A formalism of probability operators which generalizes the notion of density operator is introduced into the theory of noisy quantum systems. The Markov property and the connexion between Heisenberg and Schrödinger picture for systems undergoing an irreversible change are discussed in detail. The probability-operator treatment of noise is related to the Langevin method discussed byLax through a generalized Einstein-relation. The master equation for the quantum mechanical oscillator with linear damping is written down in a Fokker-Planck-type approximation. By means of the Einstein-relation the coefficients in the Fokker-Planck-equation are related to the parameters in the phenomenological equations.     

Numerical modeling of the photoionization of Rydberg atoms by the field of an electromagnetic wave

The ionization of excited hydrogen-like atoms in a femtosecond laser pulse is studied by direct numerical integration of the time-dependent Schrödinger equation for a quantum system in the field of an electromagnetic wave. Expressions are obtained for the ionization probability of the system as a function of the parameters of the laser pulse and the initial state of the atom. Ionization suppression is found, confirming the basic premises of the theory of interference stabilization of Rydberg atoms.     

Rigorous Solutions for the One-Dimensional Hartmann Potential in the Quantum Phase-Space Representation

The rigorous solutions of the Schrödinger equation with the one-dimensional Hartmann potential for a particle are solved and discussed within the framework of the quantum phase space representation established by Torres-Vega and Frederick. For a simple example, the uncertainty principle for the quantum probability density functions is revealed in phase space representation.     

Possibility of suppressing quantum light fluctuations when excess photon fluctuations occur inside a cavity

Using the optical excitation of a high-Q cavity as an example, it is shown that when light is observed at the output of this cavity, effective suppression of the photocurrent shot noise below the quantum limit is in general independent of the parameters of the stationary state of the field oscillator (in particular, it is independent of the rms photon fluctuations) inside the cavity and can occur not only at any allowed negative value but even at a positive value of the Mandel parameter. It was assumed in solving the problem that the cavity is optically excited by superimposing the radiation of a sub-Poisson laser and a laser with excess photon noise. A formal solution was obtained in terms of the kinetic equation for the density matrix of the actual fields (inside the laser cavities and the empty cavity), which is derived here on the basis of the Heisenberg-Langevin quantum equations, taking into account directed propagation of the field from the laser cavities inside the empty cavity. The resulting kinetic equation can also be used to solve other physical problems, since it is applicable to optical systems that contain, in principle, an arbitrary number of coupled cavities and interference mixers. Zh. éksp. Teor. Fiz. 111, 1579–1600 (May 1997)     

Quantum description of a dye laser threshold region

The renormalized Fokker-Planck equation for a dye laser is derived from the Liouvillevon Neuman equation and its stationary solution, as well as the numerically calculated line shapes and the noise spectra for the quantum threshold region are given. The calculations were performed for values of the laser parameters which enable to compare the semiclassical bistable solutions with the threshold and those without the threshold. A difference in the behavior of the line shapes in the threshold region for both cases are discussed.     

Fluctuations and stability of stationary non-equilibrium systems in detailed balance

A generalized thermodynamic potential for Markoffian systems with detailed balance and far from thermal equilibrium has been derived in a previous paper. It was shown that the principle of detailed balance is equivalent to a set of conditions fulfilled by this potential (“potential conditions”). The properties of this potential allow us to extend the validity of a number of thermodynamic concepts well known for systems in or near thermal equilibrium to stationary states far from thermal equilibrium. The concept of symmetry breaking phase transitions for these systems is introduced in analogy to thermal equilibrium systems by considering the dependence of the stationary probability density of the system on a set of externally controlled parameters {λ}. A functional of the time dependent probability density of the system is defined in close analogy to the Gibb's definition of entropy. This functional has the properties of a Ljapunov functional of the governing Fokker-Planck equation showing the stability of the stationary probability density. The Langevin equations connected with the Fokker-Planck equation are considered. It is shown that, by means of the potential conditions, generalized “thermodynamic” fluxes and forces may be defined in such a way that the smoothly varying part of the Langevin equations (kinetic equations) constitutes a linear relation between fluxes and forces. The matrix of coefficients is given by the diffusion matrix of the Fokker-Planck equation. The symmetry relations which hold for this matrix due to the potential conditions then lead to the Onsager-Casimir symmetry relations extended to systems with detailed balance near stationary states far from thermal equilibrium. Finally it is shown that under certain additional assumptions the generalized thermodynamic potential may be used as a Ljapunov function of the kinetic equations.     

Open quantum Markovian systems and the microreversibility

The concept of microreversibility (detailed balance) as applied to open quantum Markovian systems is introduced. The properties of such systems obeying detailed balance are studied in detail. It is shown how microreversibility may be used to obtain the steady state solution for the density operator. Some properties of the eigenfunctions of the Liouville operator are discussed. The general formalism is applied to Pauli master equation and some other master equations describing the relaxation of oscillators, two-level atoms in contact with heat bath; parametric frequency conversion and a single mode laser.     

基于钠硼铝硅酸盐玻璃的PbSe量子点红外单模光纤激光

根据实验制备的钠硼铝硅酸盐PbSe量子点玻璃及其透射电子显微镜(TEM)图、吸收谱和发射谱,计算机数值模拟了以PbSe量子点作为激活增益介质的红外单模光纤激光。应用遗传算法,通过数值求解粒子数速率方程和激光谐振腔振荡方程,优化计算了量子点光纤激光器(QDFL)的最佳抽运波长、光纤长度、掺杂浓度及出射镜反射率。结果表明:饱和抽运功率为2 W,在1676nm激光波长处,QDFL最大输出功率可达1.36 W,抽运效率达68%。与通常的掺稀土离子(Yb3+、Er3+)的光纤激光器相比,QDFL具有抽运效率高、激励阈值低、掺杂密度可调、光纤饱和长度短等特点。由于量子点辐射波长的尺寸依赖特性,容易形成多波长激射或波长可调的新型激光器。     

Solutions of the density matrix equation with separable interactions in Liouvillespace

The masterequation for the density matrix is solved with new techniques in Liouvillespace for two examples: 1. The stationary equation for the Lasersystem, for which the resolvent of the Liouvillean ?₀ is explicitly known, is solved by a variational procedure and 2. the exact solution of a general quantum system with simple coupling to reservoirs can be obtained, as the interaction ?₁ is separable in this case.     

Analytical expressions for asymmetric double quantum wells and their application to semiconductor heterostructures

Analytical solution of Schrödinger's equation for an asymmetric double quantum well structure is obtained for the first time in terms of a transcendental equation, roots of which give energy eigenstates and enable derivation of corresponding wave functions. Results obtained are in agreement with those reported in the literature by numerical methods and do away with the need of perturbation approximation used in such cases. As an example, GaAs/Al _x Ga_{1? x} As system has been discussed with respect to variations in barrier width, quantum well (QW) width, confining potential and material composition. The analysis shows that the coupled QWs can be matched for a given energy level to allow maximum transport between them by variation of one or more of these parameters so that the energy difference between it and the next nearest state is minimal. This condition shows a complete delocalization and the probability of finding the electron in the two QWs becomes equal, which amounts to an asymmetric system becoming symmetric with respect to probability density.     

Quantum kinetic theory of trapped particles in a strong electromagnetic field

The idea of treating quantum systems by semiclassical representations using effective quantum potentials (forces) has been successfully applied in equilibrium by many authors, see e.g. [D. Bohm, Phys. Rev. 85 (1986) 166 and 180; D.K. Ferry, J.R. Zhou, Phys. Rev. B 48 (1993) 7944; A.V. Filinov, M. Bonitz, W. Ebeling, J. Phys. A 36 (2003) 5957 and references cited therein]. Here, this idea is extended to nonequilibrium quantum systems in an external field. A gauge-invariant quantum kinetic theory for weakly inhomogeneous charged particle systems in a strong electromagnetic field is developed within the framework of nonequilibrium Green’s functions. The equation for the spectral density is simplified by introducing a classical (local) form for the kinetics. Nonlocal quantum effects are accounted for in this way by replacing the bare external confinement potential with an effective quantum potential. The equation for this effective potential is identified and solved for weak inhomogeneity in the collisionless limit. The resulting nonequilibrium spectral function is used to determine the density of states and the modification of the Born collision operator in the kinetic equation for the Wigner function due to quantum confinement effects.     

Two- and Four-Level Systems in Magnetic Fields Restricted in Time

We describe some new exact solutions for two- and four-level systems. In all the cases, external fields have a restricted behavior in time. First, we consider a method to construct new solutions for one-spin equation and give some explicit examples: One of them is in a external magnetic field that acts during a finite time interval. Then we show how these solutions can be used to solve the two-spin equation problem. A solution for two interacting spins is analyzed in the case when the field difference between the external fields in each spin varies adiabatically, vanishing on the time infinity. The latter system can be identified with a quantum gate realized by two coupled quantum dots. The probability of the Swap operation for such a gate can be explicitly expressed in terms of special functions. Using the obtained expressions, we construct plots for the Swap operation for some parameters of the external magnetic field and interaction function.     

Response and stability of a SDOF strongly nonlinear stochastic system with light damping modeled by a fractional derivative

A stochastic averaging procedure for a single-degree-of-freedom (SDOF) strongly nonlinear system with light damping modeled by a fractional derivative under Gaussian white noise excitations is developed by using the so-called generalized harmonic functions. The approximate stationary probability density and the largest Lyapunov exponent of the system are obtained from the averaged Itô stochastic differential equation of the system. It is shown that the approximate stationary solutions obtained by using the stochastic averaging procedure agree well with those from the numerical simulation of original systems. The effects of system parameters on the approxiamte stationary probability density and the largest Lyapunov exponent of the system are also discussed.     

Pauli equation for a particle in metric and electromagnetic fields

A Pauli theory (Pauli equation and definition of probability current and density) for a particle in weak metric and arbitrary electromagnetic fields is treated. To formulate non-relativistic quantum mechanical problems in arbitrary electromagnetic fields and weak metrics (non-inertial systems, gravitational fields which are distant fields of arbitrary distribution of masses, gravitational waves) it is not necessary to make use of the general-relativistic Dirac equation. Close analogies to the known Pauli theory with electromagnetic fields exist. For different metric fields the corresponding Hamiltonians are given. For quantum systems (H-atoms) which are disturbed by a homogeneous gravitational field and a gravitational wave the resulting shift of energy levels and the transition probability is calculated.