Nucleus-Electron Model for States Changing from a Liquid Metal to a Plasma and the Saha Equation

We extend the quantal hypernetted-chain (QHNC) method, which has been proved to yield accurate results for liquid metals, to treat a partially ionized plasma. In a plasma, the electrons change from a quantum to a classical fluid gradually with increasing temperature; the QHNC method applied to the electron gas is in fact able to provide the electron-electron correlation at arbitrary temperature. As an illustrating example of this approach, we investigate how liquid rubidium becomes a plasma by increasing the temperature from 0 to 30 eV at a fixed normal ion-density $1.03 \times 10^{22}/cm^3$. The electron-ion radial distribution function (RDF) in liquid Rb has distinct inner-core and outer-core parts. Even at a temperature of 1 eV, this clear distinction remains as a characteristic of a liquid metal. At a temperature of 3 eV, this distinction disappears, and rubidium becomes a plasma with the ionization 1.21. The temperature variations of bound levels in each ion and the average ionization are calculated in Rb plasmas at the same time. Using the density-functional theory, we also derive the Saha equation applicable even to a high-density plasma at low temperatures. The QHNC method provides a procedure to solve this Saha equation with ease by using a recursive formula; the charge population of differently ionized species are obtained in Rb plasmas at several temperatures. In this way, it is shown that, with the atomic number as the only input, the QHNC method produces the average ionization, the electron-ion and ion-ion RDF's, and the charge population which are consistent with the atomic structure of each ion for a partially ionized plasma.Comment: 28 pages(TeX) and 11 figures (PS

Pressure formulas for liquid metals and plasmas based on the density-functional theory

At first, pressure formulas for the electrons under the external potential produced by fixed nuclei are derived both in the surface integral and volume integral forms concerning an arbitrary volume chosen in the system; the surface integral form is described by a pressure tensor consisting of a sum of the kinetic and exchange-correlation parts in the density-functional theory, and the volume integral form represents the virial theorem with subtraction of the nuclear virial. Secondly on the basis of these formulas, the thermodynamical pressure of liquid metals and plasmas is represented in the forms of the surface integral and the volume integral including the nuclear contribution. From these results, we obtain a virial pressure formula for liquid metals, which is more accurate and simpler than the standard representation. From the view point of our formulation, some comments are made on pressure formulas derived previously and on a definition of pressure widely used.Comment: 18 pages, no figur

One-parameter extension of the Doi-Peliti formalism and relation with orthogonal polynomials

An extension of the Doi-Peliti formalism for stochastic chemical kinetics is proposed. Using the extension, path-integral expressions consistent with previous studies are obtained. In addition, the extended formalism is naturally connected to orthogonal polynomials. We show that two different orthogonal polynomials, i.e., Charlier polynomials and Hermite polynomials, can be used to express the Doi-Peliti formalism explicitly.Comment: 10 page

Polynomial solutions of nonlinear integral equations

We analyze the polynomial solutions of a nonlinear integral equation, generalizing the work of C. Bender and E. Ben-Naim. We show that, in some cases, an orthogonal solution exists and we give its general form in terms of kernel polynomials.Comment: 10 page

Phase Diagram for Anderson Disorder: beyond Single-Parameter Scaling

The Anderson model for independent electrons in a disordered potential is transformed analytically and exactly to a basis of random extended states leading to a variant of augmented space. In addition to the widely-accepted phase diagrams in all physical dimensions, a plethora of additional, weaker Anderson transitions are found, characterized by the long-distance behavior of states. Critical disorders are found for Anderson transitions at which the asymptotically dominant sector of augmented space changes for all states at the same disorder. At fixed disorder, critical energies are also found at which the localization properties of states are singular. Under the approximation of single-parameter scaling, this phase diagram reduces to the widely-accepted one in 1, 2 and 3 dimensions. In two dimensions, in addition to the Anderson transition at infinitesimal disorder, there is a transition between two localized states, characterized by a change in the nature of wave function decay.Comment: 51 pages including 4 figures, revised 30 November 200

Thermoluminescence of Simulated Interstellar Matter after Gamma-ray Irradiation

Interstellar matter is known to be strongly irradiated by radiation and several types of cosmic ray particles. Simulated interstellar matter, such as forsterite $\rm Mg_{2}SiO_{4}$, enstatite $\rm MgSiO_{3}$ and magnesite $\rm MgCO_{3}$ has been irradiated with the $\rm ^{60}Co$ gamma-rays in liquid nitrogen, and also irradiated with fast neutrons at 10 K and 70 K by making use of the low-temperature irradiation facility of Kyoto University Reactor (KUR-LTL. Maximum fast neutron dose is $10^{17}n_f{\rm /cm^{2}}$). After irradiation, samples are stored in liquid nitrogen for several months to allow the decay of induced radioactivity. We measured the luminescence spectra of the gamma ray irradiated samples during warming to 370K using a spectrophotometer. For the forsterite and magnesite, the spectra exhibit a rather intense peak at about 645 -- 655 nm and 660 nm respectively, whereas luminescence scarcely appeared in olivine sample. The spectra of forsterite is very similar to the ERE of the Red Rectangle

Probing Ion-Ion and Electron-Ion Correlations in Liquid Metals within the Quantum Hypernetted Chain Approximation

We use the Quantum Hypernetted Chain Approximation (QHNC) to calculate the ion-ion and electron-ion correlations for liquid metallic Li, Be, Na, Mg, Al, K, Ca, and Ga. We discuss trends in electron-ion structure factors and radial distribution functions, and also calculate the free-atom and metallic-atom form-factors, focusing on how bonding effects affect the interpretation of X-ray scattering experiments, especially experimental measurements of the ion-ion structure factor in the liquid metallic phase.Comment: RevTeX, 19 pages, 7 figure

A high order qq-difference equation for qq-Hahn multiple orthogonal polynomials

A high order linear $q$-difference equation with polynomial coefficients having $q$-Hahn multiple orthogonal polynomials as eigenfunctions is given. The order of the equation is related to the number of orthogonality conditions that these polynomials satisfy. Some limiting situations when $q\to1$ are studied. Indeed, the difference equation for Hahn multiple orthogonal polynomials given in \cite{Lee} is corrected and obtained as a limiting case