Maximum entropy and the problem of moments: A stable algorithm

We present a technique for entropy optimization to calculate a distribution from its moments. The technique is based upon maximizing a discretized form of the Shannon entropy functional by mapping the problem onto a dual space where an optimal solution can be constructed iteratively. We demonstrate the performance and stability of our algorithm with several tests on numerically difficult functions. We then consider an electronic structure application, the electronic density of states of amorphous silica and study the convergence of Fermi level with increasing number of moments.Comment: 4 pages including 3 figure

Entanglement and nonclassicality for multi-mode radiation field states

Nonclassicality in the sense of quantum optics is a prerequisite for entanglement in multi-mode radiation states. In this work we bring out the possibilities of passing from the former to the latter, via action of classicality preserving systems like beamsplitters, in a transparent manner. For single mode states, a complete description of nonclassicality is available via the classical theory of moments, as a set of necessary and sufficient conditions on the photon number distribution. We show that when the mode is coupled to an ancilla in any coherent state, and the system is then acted upon by a beamsplitter, these conditions turn exactly into signatures of NPT entanglement of the output state. Since the classical moment problem does not generalize to two or more modes, we turn in these cases to other familiar sufficient but not necessary conditions for nonclassicality, namely the Mandel parameter criterion and its extensions. We generalize the Mandel matrix from one-mode states to the two-mode situation, leading to a natural classification of states with varying levels of nonclassicality. For two--mode states we present a single test that can, if successful, simultaneously show nonclassicality as well as NPT entanglement. We also develop a test for NPT entanglement after beamsplitter action on a nonclassical state, tracing carefully the way in which it goes beyond the Mandel nonclassicality test. The result of three--mode beamsplitter action after coupling to an ancilla in the ground state is treated in the same spirit. The concept of genuine tripartite entanglement, and scalar measures of nonclassicality at the Mandel level for two-mode systems, are discussed. Numerous examples illustrating all these concepts are presented.Comment: Latex, 46 page

Coherent states of a charged particle in a uniform magnetic field

The coherent states are constructed for a charged particle in a uniform magnetic field based on coherent states for the circular motion which have recently been introduced by the authors.Comment: 2 eps figure

Total Widths And Slopes From Complex Regge Trajectories

Maximally complex Regge trajectories are introduced for which both Re $\alpha(s)$ and Im $\alpha(s)$ grow as $s^{1-\epsilon}$ ($\epsilon$ small and positive). Our expression reduces to the standard real linear form as the imaginary part (proportional to $\epsilon$) goes to zero. A scaling formula for the total widths emerges: $\Gamma_{TOT}/M\to$ constant for large M, in very good agreement with data for mesons and baryons. The unitarity corrections also enhance the space-like slopes from their time-like values, thereby resolving an old problem with the $\rho$ trajectory in $\pi N$ charge exchange. Finally, the unitarily enhanced intercept, $\alpha_{\rho}\approx 0.525$, \nolinebreak is in good accord with the Donnachie-Landshoff total cross section analysis.Comment: 9 pages, 3 Figure

Lyapunov exponent and natural invariant density determination of chaotic maps: An iterative maximum entropy ansatz

We apply the maximum entropy principle to construct the natural invariant density and Lyapunov exponent of one-dimensional chaotic maps. Using a novel function reconstruction technique that is based on the solution of Hausdorff moment problem via maximizing Shannon entropy, we estimate the invariant density and the Lyapunov exponent of nonlinear maps in one-dimension from a knowledge of finite number of moments. The accuracy and the stability of the algorithm are illustrated by comparing our results to a number of nonlinear maps for which the exact analytical results are available. Furthermore, we also consider a very complex example for which no exact analytical result for invariant density is available. A comparison of our results to those available in the literature is also discussed.Comment: 16 pages including 6 figure

Function reconstruction as a classical moment problem: A maximum entropy approach

We present a systematic study of the reconstruction of a non-negative function via maximum entropy approach utilizing the information contained in a finite number of moments of the function. For testing the efficacy of the approach, we reconstruct a set of functions using an iterative entropy optimization scheme, and study the convergence profile as the number of moments is increased. We consider a wide variety of functions that include a distribution with a sharp discontinuity, a rapidly oscillatory function, a distribution with singularities, and finally a distribution with several spikes and fine structure. The last example is important in the context of the determination of the natural density of the logistic map. The convergence of the method is studied by comparing the moments of the approximated functions with the exact ones. Furthermore, by varying the number of moments and iterations, we examine to what extent the features of the functions, such as the divergence behavior at singular points within the interval, is reproduced. The proximity of the reconstructed maximum entropy solution to the exact solution is examined via Kullback-Leibler divergence and variation measures for different number of moments.Comment: 20 pages, 17 figure

Linked-cluster expansion for the Green\u27s function of the infinite-U Hubbard model

We implement a highly efficient strong-coupling expansion for the Green\u27s function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green\u27s function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Padé approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green\u27s function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order

Vector Continued Fractions using a Generalised Inverse

A real vector space combined with an inverse for vectors is sufficient to define a vector continued fraction whose parameters consist of vector shifts and changes of scale. The choice of sign for different components of the vector inverse permits construction of vector analogues of the Jacobi continued fraction. These vector Jacobi fractions are related to vector and scalar-valued polynomial functions of the vectors, which satisfy recurrence relations similar to those of orthogonal polynomials. The vector Jacobi fraction has strong convergence properties which are demonstrated analytically, and illustrated numerically.Comment: Published form - minor change

Two-band random matrices

Spectral correlations in unitary invariant, non-Gaussian ensembles of large random matrices possessing an eigenvalue gap are studied within the framework of the orthogonal polynomial technique. Both local and global characteristics of spectra are directly reconstructed from the recurrence equation for orthogonal polynomials associated with a given random matrix ensemble. It is established that an eigenvalue gap does not affect the local eigenvalue correlations which follow the universal sine and the universal multicritical laws in the bulk and soft-edge scaling limits, respectively. By contrast, global smoothed eigenvalue correlations do reflect the presence of a gap, and are shown to satisfy a new universal law exhibiting a sharp dependence on the odd/even dimension of random matrices whose spectra are bounded. In the case of unbounded spectrum, the corresponding universal `density-density' correlator is conjectured to be generic for chaotic systems with a forbidden gap and broken time reversal symmetry.Comment: 12 pages (latex), references added, discussion enlarge

Generating Bounds for the Ground State Energy of the Infinite Quantum Lens Potential

Moment based methods have produced efficient multiscale quantization algorithms for solving singular perturbation/strong coupling problems. One of these, the Eigenvalue Moment Method (EMM), developed by Handy et al (Phys. Rev. Lett.{\bf 55}, 931 (1985); ibid, {\bf 60}, 253 (1988b)), generates converging lower and upper bounds to a specific discrete state energy, once the signature property of the associated wavefunction is known. This method is particularly effective for multidimensional, bosonic ground state problems, since the corresponding wavefunction must be of uniform signature, and can be taken to be positive. Despite this, the vast majority of problems studied have been on unbounded domains. The important problem of an electron in an infinite quantum lens potential defines a challenging extension of EMM to systems defined on a compact domain. We investigate this here, and introduce novel modifications to the conventional EMM formalism that facilitate its adaptability to the required boundary conditions.Comment: Submitted to J. Phys.