Accurate calculation of resonances in multiple-well oscillators

Quantum--mechanical multiple--well oscillators exhibit curious complex eigenvalues that resemble resonances in models with continuum spectra. We discuss a method for the accurate calculation of their real and imaginary parts

A model for single electron decays from a strongly isolated quantum dot

Recent measurements of electron escape from a non-equilibrium charged quantum dot are interpreted within a 2D separable model. The confining potential is derived from 3D self-consistent Poisson-Thomas-Fermi calculations. It is found that the sequence of decay lifetimes provides a sensitive test of the confining potential and its dependence on electron occupation.Comment: 9 pages, 10 figure

The optimized Rayleigh-Ritz scheme for determining the quantum-mechanical spectrum

The convergence of the Rayleigh-Ritz method with nonlinear parameters optimized through minimization of the trace of the truncated matrix is demonstrated by a comparison with analytically known eigenstates of various quasi-solvable systems. We show that the basis of the harmonic oscillator eigenfunctions with optimized frequency ? enables determination of boundstate energies of one-dimensional oscillators to an arbitrary accuracy, even in the case of highly anharmonic multi-well potentials. The same is true in the spherically symmetric case of V (r) = {\omega}2r2 2 + {\lambda}rk, if k > 0. For spiked oscillators with k < -1, the basis of the pseudoharmonic oscillator eigenfunctions with two parameters ? and {\gamma} is more suitable, and optimization of the latter appears crucial for a precise determination of the spectrum.Comment: 22 pages,8 figure

Eigenvalue bounds for a class of singular potentials in N dimensions

The eigenvalue bounds obtained earlier [J. Phys. A: Math. Gen. 31 (1998) 963] for smooth transformations of the form V(x) = g(x^2) + f(1/x^2) are extended to N-dimensions. In particular a simple formula is derived which bounds the eigenvalues for the spiked harmonic oscillator potential V(x) = x^2 + lambda/x^alpha, alpha > 0, lambda > 0, and is valid for all discrete eigenvalues, arbitrary angular momentum ell, and spatial dimension N.Comment: 10 pages (plain tex with 2 ps figures). J.Phys.A:Math.Gen.(In Press

A basis for variational calculations in d dimensions

In this paper we derive expressions for matrix elements (\phi_i,H\phi_j) for the Hamiltonian H=-\Delta+\sum_q a(q)r^q in d > 1 dimensions. The basis functions in each angular momentum subspace are of the form phi_i(r)=r^{i+1+(t-d)/2}e^{-r^p/2}, i >= 0, p > 0, t > 0. The matrix elements are given in terms of the Gamma function for all d. The significance of the parameters t and p and scale s are discussed. Applications to a variety of potentials are presented, including potentials with singular repulsive terms of the form b/r^a, a,b > 0, perturbed Coulomb potentials -D/r + B r + Ar^2, and potentials with weak repulsive terms, such as -g r^2 + r^4, g > 0.Comment: 22 page

Variational collocation for systems of coupled anharmonic oscillators

We have applied a collocation approach to obtain the numerical solution to the stationary Schr\"odinger equation for systems of coupled oscillators. The dependence of the discretized Hamiltonian on scale and angle parameters is exploited to obtain optimal convergence to the exact results. A careful comparison with results taken from the literature is performed, showing the advantages of the present approach.Comment: 14 pages, 10 table

Convergence of the Gaussian Expansion Method in Dimensionally Reduced Yang-Mills Integrals

We advocate a method to improve systematically the self-consistent harmonic approximation (or the Gaussian approximation), which has been employed extensively in condensed matter physics and statistical mechanics. We demonstrate the {\em convergence} of the method in a model obtained from dimensional reduction of SU($N$) Yang-Mills theory in $D$ dimensions. Explicit calculations have been carried out up to the 7th order in the large-N limit, and we do observe a clear convergence to Monte Carlo results. For $D \gtrsim 10$ the convergence is already achieved at the 3rd order, which suggests that the method is particularly useful for studying the IIB matrix model, a conjectured nonperturbative definition of type IIB superstring theory.Comment: LaTeX, 4 pages, 5 figures; title slightly changed, explanations added (16 pages, 14 figures), final version published in JHE

Spiked oscillators: exact solution

A procedure to obtain the eigenenergies and eigenfunctions of a quantum spiked oscillator is presented. The originality of the method lies in an adequate use of asymptotic expansions of Wronskians of algebraic solutions of the Schroedinger equation. The procedure is applied to three familiar examples of spiked oscillators

Eigenvalue bounds for polynomial central potentials in d dimensions

If a single particle obeys non-relativistic QM in R^d and has the Hamiltonian H = - Delta + f(r), where f(r)=sum_{i = 1}^{k}a_ir^{q_i}, 2\leq q_i < q_{i+1}, a_i \geq 0$, then the eigenvalues E = E_{n\ell}^{(d)}(\lambda) are given approximately by the semi-classical expression E = \min_{r > 0}[\frac{1}{r^2} + \sum_{i = 1}^{k}a_i(P_ir)^{q_i}]. It is proved that this formula yields a lower bound if P_i = P_{n\ell}^{(d)}(q_1), an upper bound if$P_i = P_{n\ell}^{(d)}(q_k) and a general approximation formula if P_i = P_{n\ell}^{(d)}(q_i). For the quantum anharmonic oscillator f(r)=r^2+\lambda r^{2m},m=2,3,... in d dimension, for example, E = E_{n\ell}^{(d)}(\lambda) is determined by the algebraic expression \lambda={1\over \beta}({2\alpha(m-1)\over mE-\delta})^m({4\alpha \over (mE-\delta)}-{E\over (m-1)}) where \delta={\sqrt{E^2m^2-4\alpha(m^2-1)}} and \alpha, \beta are constants. An improved lower bound to the lowest eigenvalue in each angular-momentum subspace is also provided. A comparison with the recent results of Bhattacharya et al (Phys. Lett. A, 244 (1998) 9) and Dasgupta et al (J. Phys. A: Math. Theor., 40 (2007) 773) is discussed.Comment: 13 pages, no figure

The Fokker-Planck equation for bistable potential in the optimized expansion

The optimized expansion is used to formulate a systematic approximation scheme to the probability distribution of a stochastic system. The first order approximation for the one-dimensional system driven by noise in an anharmonic potential is shown to agree well with the exact solution of the Fokker-Planck equation. Even for a bistable system the whole period of evolution to equilibrium is correctly described at various noise intensities.Comment: 12 pages, LATEX, 3 Postscript figures compressed an