Asymptotic approximations to the nodes and weights of Gauss-Hermite and Gauss-Laguerre quadratures

Asymptotic approximations to the zeros of Hermite and Laguerre polynomials are given, together with methods for obtaining the coefficients in the expansions. These approximations can be used as a standalone method of computation of Gaussian quadratures for high enough degrees, with Gaussian weights computed from asymptotic approximations for the orthogonal polynomials. We provide numerical evidence showing that for degrees greater than $100$ the asymptotic methods are enough for a double precision accuracy computation ($15$-$16$ digits) of the nodes and weights of the Gauss--Hermite and Gauss--Laguerre quadratures.Comment: Submitted to Studies in Applied Mathematic

Computation of the Marcum Q-function

Methods and an algorithm for computing the generalized Marcum $Q-$function ($Q_{\mu}(x,y)$) and the complementary function ($P_{\mu}(x,y)$) are described. These functions appear in problems of different technical and scientific areas such as, for example, radar detection and communications, statistics and probability theory, where they are called the non-central chi-square or the non central gamma cumulative distribution functions. The algorithm for computing the Marcum functions combines different methods of evaluation in different regions: series expansions, integral representations, asymptotic expansions, and use of three-term homogeneous recurrence relations. A relative accuracy close to $10^{-12}$ can be obtained in the parameter region $(x,y,\mu) \in [0,\,A]\times [0,\,A]\times [1,\,A]$, $A=200$, while for larger parameters the accuracy decreases (close to $10^{-11}$ for $A=1000$ and close to $5\times 10^{-11}$ for $A=10000$).Comment: Accepted for publication in ACM Trans. Math. Soft

Fermions on one or fewer Kinks

We find the full spectrum of fermion bound states on a Z_2 kink. In addition to the zero mode, there are int[2 m_f/m_s] bound states, where m_f is the fermion and m_s the scalar mass. We also study fermion modes on the background of a well-separated kink-antikink pair. Using a variational argument, we prove that there is at least one bound state in this background, and that the energy of this bound state goes to zero with increasing kink-antikink separation, 2L, and faster than e^{-a2L} where a = min(m_s, 2 m_f). By numerical evaluation, we find some of the low lying bound states explicitly.Comment: 7 pages, 4 figure

Non-equilibrium dynamics of a Bose-Einstein condensate in an optical lattice

The dynamical evolution of a Bose-Einstein condensate trapped in a one-dimensional lattice potential is investigated theoretically in the framework of the Bose-Hubbard model. The emphasis is set on the far-from-equilibrium evolution in a case where the gas is strongly interacting. This is realized by an appropriate choice of the parameters in the Hamiltonian, and by starting with an initial state, where one lattice well contains a Bose-Einstein condensate while all other wells are empty. Oscillations of the condensate as well as non-condensate fractions of the gas between the different sites of the lattice are found to be damped as a consequence of the collisional interactions between the atoms. Functional integral techniques involving self-consistently determined mean fields as well as two-point correlation functions are used to derive the two-particle-irreducible (2PI) effective action. The action is expanded in inverse powers of the number of field components N, and the dynamic equations are derived from it to next-to-leading order in this expansion. This approach reaches considerably beyond the Hartree-Fock-Bogoliubov mean-field theory, and its results are compared to the exact quantum dynamics obtained by A.M. Rey et al., Phys. Rev. A 69, 033610 (2004) for small atom numbers.Comment: 9 pages RevTeX, 3 figure

Resonantly enhanced pair production in a simple diatomic model

A new mechanism for the production of electron-positron pairs from the interaction of a laser field and a fully stripped diatomic molecule in the tunneling regime is presented. When the laser field is turned off, the Dirac operator has resonances in both the positive and the negative energy continua while bound states are in the mass gap. When this system is immersed in a strong laser field, the resonances move in the complex energy plane: the negative energy resonances are pushed to higher energies while the bound states are Stark shifted. It is argued here that there is a pair production enhancement at the crossing of resonances by looking at a simple 1-D model: the nuclei are modeled simply by Dirac delta potential wells while the laser field is assumed to be static and of finite spatial extent. The average rate for the number of electron-positron pairs produced is evaluated and the results are compared to the single nucleus and to the free cases. It is shown that positrons are produced by the Resonantly Enhanced Pair Production (REPP) mechanism, which is analogous to the resonantly enhanced ionization of molecular physics. This phenomenon could be used to increase the number of pairs produced at low field strength, allowing the study of the Dirac vacuum.Comment: 11 pages, 4 figure

Power Utility Maximization in Discrete-Time and Continuous-Time Exponential Levy Models

Consider power utility maximization of terminal wealth in a 1-dimensional continuous-time exponential Levy model with finite time horizon. We discretize the model by restricting portfolio adjustments to an equidistant discrete time grid. Under minimal assumptions we prove convergence of the optimal discrete-time strategies to the continuous-time counterpart. In addition, we provide and compare qualitative properties of the discrete-time and continuous-time optimizers.Comment: 18 pages, to appear in Mathematical Methods of Operations Research. The final publication is available at springerlink.co

Landau-Zener-St\"uckelberg interferometry in pair production from counterpropagating lasers

The rate of electron-positron pair production in linearly polarized counter-propagating lasers is evaluated from a recently discovered solution of the time-dependent Dirac equation. The latter is solved in momentum space where it is formally equivalent to the Schr\"odinger equation describing a strongly driven two-level system. The solution is found from a simple transformation of the Dirac equation and is given in compact form in terms of the doubly-confluent Heun's function. By using the analogy with the two-level system, it is shown that for high-intensity lasers, pair production occurs through periodic non-adiabatic transitions when the adiabatic energy gap is minimal. These transitions give rise to an intricate interference pattern in the pair spectrum, reminiscent of the Landau-Zener-St\"uckelberg phenomenon in molecular physics: the accumulated phase result in constructive or destructive interference. The adiabatic-impulse model is used to study this phenomenon and shows an excellent agreement with the exact result.Comment: 22 pages, 7 figure

The Schrodinger equation with Hulthen potential plus ring-shaped potential

We present the solutions of the Schr$\ddot{o}$dinger equation with the Hulth$\acute{e}$n potential plus ring-shape potential for $\ell\neq 0$ states within the framework of an exponential approximation of the centrifugal potential.Solutions to the corresponding angular and radial equations are obtained in terms of special functions using the conventional Nikiforov-Uvarov method. The normalization constant for the Hulth$\acute{e}$n potential is also computed.Comment: Typed with LateX,12 Pages, Typos correcte

Topological Phases in Graphitic Cones

The electronic structure of graphitic cones exhibits distinctive topological features associated with the apical disclinations. Aharonov-Bohm magnetoconductance oscillations (period Phi_0) are completely absent in rings fabricated from cones with a single pentagonal disclination. Close to the apex, the local density of states changes qualitatively, either developing a cusp which drops to zero at the Fermi energy, or forming a region of nonzero density across the Fermi energy, a local metalization of graphene.Comment: 4 pages, RevTeX 4, 3 PostScript figure