Mean spherical model for the square-well potential

The mean spherical model (MSM) equation has been solved for the aquare-well potential. Comparison of the MSM results with the simulation results shows the MSM radial distribution function and energy to be satisfactory in the liquid states.     

Mean spherical model for asymmetric electrolytes

The mean spherical approximation for the primitive model of electrolytic solutions is solved for the most general case of an arbitrary charge and size. It is found that the excess thermodynamic properties are scaled to the charge and the size by means of a rational expression involving a single parameter. This parameter is found by solving an algebraic equation. The explicit form of this equation is obtained for the cases of a binary mixture and in the limit of low concentrations.     

Mean spherical model-structure of liquid argon

A new formulation forS(k), the structure factor, has been obtained, treating the square-well potential as a perturbation on the hard-core in the mean spherical model approximation. The potential parameters have been varied not only to get a satisfactory peak height at the fitted temperature but also over a wide range of temperatures and densities. The agreement between the computed and experimental structure factor values shows that the representation of the attractive tail by a square-well potential yields a satisfactory method in understanding the structure of liquid argon.     

Stability analysis of a spherical electrostatic suspension in a liquid in the induction approximation

In [1] the image force was shown to impose additional conditions for the electrostatic suspension of a sphere without dynamic control of the electrode potential, and the dependence of the critical voltage between the electrodes on the sphere radius was derived experimentally. In this work, this dependence is found analytically by calculating electrical forces in the third-order approximation in shift of the sphere from equilibrium.     

Analytic solution of the mean spherical approximation for ion-dipole model in a neutralizing background

The analytic solution of the mean spherical approximation (MSA) for a multicomponent mixture of hard ions and hard dipoles with arbitrary valences and sizes of particles in a uniform neutralizing background is found. Expressions for the pair correlation functions and thermodynamics in the MSA are obtained.     

Mean spherical model for hard ions and dipoles: Thermodynamics and correlation functions

The solution of the mean spherical model of a mixture of equal-size hard ions and dipoles is reinvestigated. Simple expressions for the coefficients of the Laplace transform of the pair correlation function and the other thermodynamic properties are given.Supported by DOE through the Center of Environment and Energy Research of the University of Puerto Rico and NSF Grant CH 77-14611.     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

Abstract

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

Extended mean spherical approximation for electrolyte solutions

The range of applicability of the mean spherical approximation to a primitive-model electrolyte is extended to allow for stronger Coulomb coupling. The new terms in the approximation are extracted from the HNC formalism, which is adequate for strongly coupled Coulomb systems. The new approximation satisfies the Stillinger-Lovett conditions. The new integral equation thus obtained is presented.     

Scale-free memory model for multiagent reinforcement learning. Mean field approximation and rock-paper-scissors dynamics

A continuous time model for multiagent systems governed by reinforcement learning with scale-free memory is developed. The agents are assumed to act independently of one another in optimizing their choice of possible actions via trial-and-error search. To gain awareness about the action value the agents accumulate in their memory the rewards obtained from taking a specific action at each moment of time. The contribution of the rewards in the past to the agent current perception of action value is described by an integral operator with a power-law kernel. Finally a fractional differential equation governing the system dynamics is obtained. The agents are considered to interact with one another implicitly via the reward of one agent depending on the choice of the other agents. The pairwise interaction model is adopted to describe this effect. As a specific example of systems with non-transitive interactions, a two agent and three agent systems of the rock-paper-scissors type are analyzed in detail, including the stability analysis and numerical simulation. Scale-free memory is demonstrated to cause complex dynamics of the systems at hand. In particular, it is shown that there can be simultaneously two modes of the system instability undergoing subcritical and supercritical bifurcation, with the latter one exhibiting anomalous oscillations with the amplitude and period growing with time. Besides, the instability onset via this supercritical mode may be regarded as “altruism self-organization”. For the three agent system the instability dynamics is found to be rather irregular and can be composed of alternate fragments of oscillations different in their properties.     

Simple pair potential model for real fluids. III. Parameter determination and a revised model for spherical molecules

A method combining macroscopic data with atomic structure constants for determination of the parameters of simple potential models is discussed. A revised hybrid potential model, MMH, whose hard core diameter and location of the potential minimum are determined a priori from the atomic structure is proposed. Calculations of the second and third virial coefficients and the viscosity coefficient indicate that the MMH model may be a reliable potential superior to all other simple potential models.     

Comment on the rank two spherical harmonic model for nematic liquid crystals

Both order parameters of the uniaxial nematic phase are calculated from the conventional rank two spherical harmonic model. The conclusion, based on the calculation and experimental values of the order parameters, is drawn that the rank two model fails more seriously than previously realized to describe the nematic state.     

Critical parameters of asymmetric primitive model electrolytes in the mean spherical approximation

The effect of size and charge asymmetries on the critical parameters of primitive model electrolytes is comprehensively surveyed by using the mean spherical approximation.     

Mean field approximation for noisy delay coupled excitable neurons

Mean field approximation of a large collection of FitzHugh-Nagumo excitable neurons with noise and all-to-all coupling with explicit time-delays, modelled by N?1 stochastic delay-differential equations is derived. The resulting approximation contains only two deterministic delay-differential equations but provides excellent predictions concerning the stability and bifurcations of the averaged global variables of the exact large system.     

Accurate delta potential approximation for a coordinate-dependent potential and its analytical solution

In this Letter, analytical expression in recursive form for reflection amplitude is derived through a transformation which reduces the equation of motion for a coordinate-dependent potential to an equation with a chain of delta function potentials. The formulation provides accurate results of eigenvalues of bound and resonance states generated by a variety of potentials.     

The self-consistent mean spherical approximation for the one-component plasma

The consequences of choosing the adjustable hard-core diameter in the mean spherical approximation for the one-component plasma so as to achieve thermodynamic consistency between the energy and compressibility equations are investigated. Such a choice is found to be possible only for >8.5 and, although the resulting correlation functions are discontinuous, the height of the main peak in the static structure factor is remarkably accurate. Two especially noteworthy aspects of the thermodynamic results are that the compressibility equation is much more accurate than in any previous approximation free of input from computer simulations and that the nonstatic part of the internal energy has a ^1/4 dependence in the strong coupling limit in agreement with Monte Carlo data.     

Density functional theory for a model of nonuniform liquid metal in partially ionized states

Summary The problem of the electronic, ionic and atomic density profiles in a nonuniform liquid metal which can locally exist in partially ionized states is examined using the density functional formalism. Ions and electrons are allowed to bind into atoms through a ?reaction? governed by the mass action law. Formally exact equations for the density profiles are given in terms of the inverse response matrix of the nonuniform system, which consists of two terms: the first corresponding to to a mixture of ions, electrons and structureless atoms and the other to the atomic internal degrees of freedom. Approximate schemes are proposed for both contributions, stressing, in partieular, i) how ionization of the atoms arises with increasing density and the relation with Mott’s criterion for the metal-insulator transition, and ii) the usefulness of a weak-coupling assumption for interspecies correlations. This formalism used together with properly parametrized trial functions for the density profiles should be particularly useful for studying the liquid-vapour interface of alkali metals.

---

Riassunto Il lavoro usa il formalismo del funzionale di densità per discutere la determinazione dei profili di densità in un fluido non omogeneo che può trovarsi localmente in stati di ionizzazione parziale. Con riferimento specifico ai metalli alcalini, la formazione di atomi per legame di ioni ed electtroni è descritta come una ?reazione? governata dalla legge di azione di massa. Le condizioni di equilibrio microscopico del fluido sono espresse tramite una matrice di risposta in cui appaiono due contributi, associati rispettivamente ai gradi di libertà interni degli atomi e ad una miscela di tre specie di particelle senza struttura interna. Si propongono quindi schemi approssimati per il calcolo di questi due contributi, con attenzione particolare all’uso di teorie perturbative e alla relazione tra grado di ionizzazione locale e densità localle, in connessione col criterio di Mott per la transizione isolante-metallo. Il formalismo proposto è in particolare appropriato per lo studio dell’interfaccia liquido-vapore in metalli alcalini.

---

Резюме Используя функционала плотности, исследуется проблема профилей электронной, ионной и атомой ллотностей В неоднодном зидком металле, которий мозет локалщно существоватщ в частично иониэированных состояних. Ионы и электроны могут şвяэыватщся в атомщ череэ ?реакцию, определяемую эаконом дейşтвуюших масс. Приводятся формалщно точные уравнения для профнлеи плотности в термина⇆ обратной матрици отклика для неоднородной системы, которая состоит иэ двух чуленов: первий соответствует ссмеси ионов, электронов и бесстрктурных атмов и второй член соответствует атомным внутренним степеням свобды. Предлагаются приближенные схемы для обоых вкладов, обрашая особое бнмание, в частности, 1) как иониэация атомов воэликает бри увеличени влотности, на свяэы с критерием Мота для берехода ?металл-иеолятор?, и 2) на полеэносты предположения слабой свяэи для межвидовых кореляций. Предложенный Формалиэм, исполыеовалный вместе с параметриэованными пробными функциями для профилей плотнорти, окаэываетця оробенне полеэным для исследования грраницы ?жидкостъ-пар? в шелочлых металлах.

     

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 an 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.03x10(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 that are consistent with the atomic structure of each ion for a partially ionized plasma.