BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Jean-Christophe Pain appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We propose a simple derivation of an upper bound for the perimeter of an ellipse. The procedure, which relies on the use of elliptic integrals, consists in introducing, via inequalities and convexity properties, specific integrals which can be calculated analytically.
In the Ziman-Evans formulation, the electrical resistivity involves several quantities like the plasma mean ionization and chemical potential, the electron-ion scattering cross-section, the ion-ion structure factor, and the derivative of the Fermi distribution with respect to electron energy. Therefore, in order to make significant progress while comparing different models or analyzing experimental measurements, it is important to try to get insight into such partial physical quantities. In the present work, we propose to investigate the sensitivity of the resistivity to the different physical quantities required for its computation.
In previous works, we presented series representations for $π^3$ and $π^5$, in which the prefactor depends only on the golden ratio appears. In this article, we derive a general relation involving trigonometric functions and an infinite series. Such an identity is likely to provide many series representations for any positive power of $π$, among them the above mentioned representations for $π^3$ and $π^5$.
In many fields of statistical physics, for instance in the study of the liquid-gas phase transition in finite nuclear matter, the Virial coefficients of the Fermi gas play a major role. In this note, we provide relations, sum rules, analytical formulas and numerical values for such coefficients.
We provide a family of expressions of $π$ in terms of the golden ratio $ϕ$ in the same spirit of the formula obtained by Bailey, Borwein and Plouffe for $π$. Connection with cyclotomic polynomials is outlined.
Although many series exist for $π$ and $π^2$, very few are known for $π^3$. In this article, we derive, using a trigonometric identity obtained by Euler, two representations of $π^3$ involving infinite sums and the golden ratio. The methodology can be generalized in order to obtain further series, relating by the way $π^3$ to other mathematical constants.
In this didactic note, we describe a procedure to derive successive approximations of $π$ using Euler Beta functions. It is an interesting exercise for undergraduate students, since it involves polynomial roots, integral calculations, inequalities and Euler Beta functions.
We present sum rules for Clebsch-Gordan coefficients in the framework of SO(4) group-theoretical description of the hydrogen atom. The main results are obtained using properties of the Runge-Lenz- Pauli vector, in particular expressing the matrix elements of the powers of its last component both in spherical and parabolic basis. Connections with Stark effect and diamagnetism of the hydrogen atom are outlined.
We present recent evolutions of the detailed opacity code SCO-RCG which combines statistical modelings of levels and lines with fine-structure calculations. The code now includes the Partially-Resolved-Transition-Array model, which allows one to replace a complex transition array by a small-scale detailed calculation preserving energy and variance of the genuine transition array and yielding improved high-order moments. An approximate method for studying the impact of strong magnetic field on opacity and emissivity was also recently implemented. The Zeeman line profile is modeled by fourth-order Gram-Charlier expansion series, which is a Gaussian multiplied by a linear combination of Hermite polynomials. Electron collisional line broadening is often modeled by a Lorentzian function and one has to calculate the convolution of a Lorentzian with Gram-Charlier distribution for a huge number of spectral lines. Since the numerical cost of the direct convolution would be prohibitive, we propose, in order to obtain the resulting profile, a fast and precise algorithm, relying on a representation of the Gaussian by cubic splines.
In a recent paper [J.-C. Pain, Opt. Spectrosc. ${\bf 218}$, 1105-1109 (2020)], we discussed the link between expectation values of powers of $r$ and Clebsch-Gordan coefficients. In this short note we provide additional information, reminding that such a connection is a direct consequence of group theory. The hydrogenic radial wavefunctions form bases for infinite dimensional representations of the algebra of the non-compact group $O(2,1)$ and the expectation values $r^p$ and $r^{-p}$ ($p$ being positive) transform as tensors with respect to this algebra. As shown a long time ago by Armstrong [L. Armstrong Jr., J . Phys. (Paris) Suppl. C 4 ${\bf 31}$, 17 (1970)], analysis of matrix elements of $r^p$ and $r^{-p}$ reveals that the Wigner-Eckart theorem is valid for this group and that the corresponding Clebsch-Gordan coefficients are proportional to the usual $SO(3)$ Clebsch-Gordan coefficients. This proportionality provides simple explanations of the selection rules for hydrogenic radial matrix elements pointed out by Pasternack and Sternheimer, and the proportionality of hydrogenic expectation values of $r^p$ and $r^{-p}$ to $3jm$ symbols.
The contribution of Jacques Raynal to angular-momentum theory is highly valuable. In the present article, I intend to recall the main aspects of his work related to Wigner $3j$ symbols. It is well known that the latter can be expressed with a hypergeometric series. The polynomial zeros of the $3j$ coefficients were initially characterized by the number of terms of the series minus one, which is the degree of the coefficient. A detailed study of the zeros of the $3j$ coefficient with respect to the degree $n$ for $J = a + b + c \leq 240$ ($a$, $b$ and $c$ being the angular momenta in the first line of the $3j$ symbol) by Raynal revealed that most zeros of high degree had small magnetic quantum numbers. This led him to define the order $m$ to improve the classification of the zeros of the $3j$ coefficient. Raynal did a search for the polynomial zeros of degree 1 to 7 and found that the number of zeros of degree 1 and 2 are infinite, though the number of zeros of degree larger than 3 decreases very quickly as the degree increases. Based on Whipple's transformations of hypergeometric $_3F_2$ functions with unit argument, Raynal generalized the Wigner $3j$ symbols to any arguments and pointed out that there are twelve sets of ten formulas (twelve sets of 120 generalized $3j$ symbols) which are equivalent in the usual case. In this paper, we also discuss other aspects of the zeros of $3j$ coefficients, such as the role of Diophantine equations and powerful numbers, or the alternative approach involving Labarthe patterns.
We present three explicit formulas for the number of electronic configurations in an atom, i.e. the number of ways to distribute $Q$ electrons in $N$ subshells of respective degeneracies $g_1$, $g_2$, ..., $g_N$. The new expressions are obtained using the generating-function formalism. The first one contains sums involving multinomial coefficients. The second one relies on the idea of gathering subshells having the same degeneracy. A third one also collects subshells with the same degeneracy and leads to the definition of a two-variable generating function, allowing the derivation of recursion relations. Concerning the distribution of population on $N$ distinct subshells of a given degeneracy $g$, analytical expressions for the first moments of this distribution are given. The general case of subshells with any degeneracy is analyzed through the computation of cumulants. A fairly simple expression for the cumulants at any order is provided, as well as the cumulant generating function. Using Gram-Charlier expansion, simple approximations of the analyzed distribution in terms of a normal distribution multiplied by a sum of Hermite polynomials are given. The Edgeworth expansion has also been tested. Its accuracy is equivalent to the Gram-Charlier accuracy when few terms are kept, but it is much more rapidly divergent when the truncation order increases. While this analysis is illustrated by examples in atomic supershells it also applies to more general combinatorial problems such as fermion distributions.
In this short paper, the commutator of monomials of operators obeying constant commutation relations is expressed in terms of anticommutators. The formula involves Bernoulli numbers or Euler polynomials evaluated in zero. The role of Bernoulli numbers in quantum-mechanical identities such as the Baker-Campbell-Hausdorff formula is emphasized and applications connected to ordering problems as well as to Ehrenfest theorem are proposed.
The expression of the electron broadening operator including the effect of penetrating collisions, i.e. for which the incoming electron enters the extent of bound-electron wavefunctions, is rather complicated, even for hydrogen. It involves integrals of special functions, which evaluation deserves scrutiny. We present a simple approximate form of the electron collision operator for hydrogen including penetration effects, both in direct and interference terms. The new expression is accurate and easy to compute. In the Penetration Standard Theory, the collision operator is convergent whatever the value of the maximum impact parameter. However, when penetration theory is not valid anymore, it should be questioned. We discuss the problem of strong collisions when penetration effects are taken into account.
We present analytical formulas for the vacuum-polarization Uehling potential in the case where the finite size of the nucleus is modeled by a Fermi charge distribution. Using a Sommerfeld-type development, the potential is expressed in terms of multiple derivatives of a particular integral. The latter and its derivatives can be evaluated exactly in terms of Bickley-Naylor functions, which connection to the Uehling potential was already pointed out in the pure Coulomb case, and of usual Bessel functions of the second kind. The cusp and asymptotic expressions for the Uehling potential with a Fermi charge distribution are also provided. Analytical results for the higher-order-contribution Källèn-Sabry potential are given.
Laboratory (laser and Z-pinch) opacity measurements of well-characterized plasmas provide data to assist inertial confinement fusion, astrophysics and atomic-physics research. In order to test the atomic-physics codes devoted to the calculation of radiative properties of hot plasmas, such experiments must fulfill a number of requirements. In this work, we discuss some sources of uncertainty in absorption-spectroscopy experiments, concerning areal mass, background emission, intensity of the backlighter and self-emission of the plasma. We also study the impact of spatial non-uniformities of the sample.
In multiply-charged ion plasmas, a significant number of electrons may occupy high-energy orbitals. These "Rydberg" electrons, when they act as spectators, are responsible for a number of satellites of X-ray absorption or emission lines, yielding a broadening of the red wing of the resonance lines. The contribution of such satellite lines may be important, because of the high degeneracy of the relevant excited configurations which give these large Boltzmann weights. However, it is in general difficult to take these configurations into account since they are likely to give rise to a large number of lines. We propose to model the perturbation induced by the spectators in a way similar to the Partially-Resolved-Transition-Array approach recently published by C. Iglesias. It consists in a partial detailed-line-accounting calculation in which the effect of the Rydberg spectators is included through a shift and width, expressed in terms of the canonical partition functions, which are key-ingredients of the Super-Transition-Arrays model. The resulting method can a priori be used in any detailed-configuration/line-accounting opacity code.
The broadening of lines by Stark effect is an important tool for inferring electron density and temperature in plasmas. Stark-effect calculations often rely on atomic data (transition rates, energy levels,...) not always exhaustive and/or valid for isolated atoms. We present a recent development in the detailed opacity code SCO-RCG for K-shell spectroscopy (hydrogen- and helium-like ions). This approach is adapted from the work of Gilles and Peyrusse. Neglecting non-diagonal terms in dipolar and collision operators, the line profile is expressed as a sum of Voigt functions associated to the Stark components. The formalism relies on the use of parabolic coordinates within SO(4) symmetry. The relativistic fine-structure of Lyman lines is included by diagonalizing the hamiltonian matrix associated to quantum states having the same principal quantum number $n$. The resulting code enables one to investigate plasma environment effects, the impact of the microfield distribution, the decoupling between electron and ion temperatures and the role of satellite lines (such as Li-like $1sn\ell n'\ell' - 1s^2n\ell$, Be-like, etc.). Comparisons with simpler and widely-used semi-empirical models are presented.
Absorption and emission spectra of multicharged-ion plasmas contain a huge number of electron configurations and electric-dipolar lines, which can be handled by global methods. However, some transition arrays consist only of a small bunch of lines. For that reason, we developed the hybrid opacity code SCO-RCG combining the (statistical) super-transition-array method and the (detailed) fine-structure calculation (requiring the diagonalization of the Hamiltonian matrix) of atomic structure. In order to decide whether a detailed treatment of lines is necessary and to determine the validity of statistical methods, the code involves criteria taking into account coalescence of lines and porosity (localized absence of lines) in transition arrays. Data required for the calculation of detailed transition arrays (Slater, spin-orbit and dipolar integrals) are provided by the super-configuration code SCO, which takes into account plasma screening effects on wavefunctions. Then, level energies and lines are calculated by Cowan's code. Transition arrays for which a detailed treatment is not required or impossible are described statistically, by UTA (Unresolved Transition Arrays) and STA (Super Transition Arrays) formalisms used in SCO. Recent developments are presented, such as the extension of the PRTA (Partially Resolved Transition Array) model to the hybrid approach, as well as comparisons with experimental spectra (laser or Z-pinch).
The modeling of complex atomic spectra is a difficult task, due to the huge number of levels and lines involved. In the presence of a magnetic field, the computation becomes even more difficult. The anomalous Zeeman pattern is a superposition of many absorption or emission profiles with different Zeeman relative strengths, shifts, widths, asymmetries and sharpnesses. We propose a statistical approach to study the effect of a magnetic field on the broadening of spectral lines and transition arrays in atomic spectra. In this model, the sigma and pi profiles are described using the moments of the Zeeman components, which depend on quantum numbers and Landé factors. A graphical calculation of these moments, together with a statistical modeling of Zeeman profiles as expansions in terms of Hermite polynomials are presented. It is shown that the procedure is more efficient, in terms of convergence and validity range, than the Taylor-series expansion in powers of the magnetic field which was suggested in the past. Finally, a simple approximate method to estimate the contribution of a magnetic field to the width of transition arrays is proposed. It relies on our recently published recursive technique for the numbering of LS-terms of an arbitrary configuration.
Around 1950, Wigner introduced the idea of modelling physical reality with an ensemble of random matrices while studying the energy levels of heavy atomic nuclei. Since then, the field of random-matrix theory has grown tremendously, with applications ranging from fluctuations on the economic markets to complex atomic spectra. The purpose of this short article is to review several attempts to apply the basic concepts of random-matrix theory to the structure and radiative transitions of atoms and ions, using the random matrices originally introduced by Wigner in the framework of the gaussian orthogonal ensemble. Some intrinsic properties of complex-atom physics, which could be enlightened by random-matrix theory, are presented.
In hot plasmas, a temperature of a few tens of eV is sufficient for producing highly stripped ions where multipole transitions become important. At low density, the transitions from tightly bound inner shells lead to electric-quadrupole (E2) lines which are comparable in strength with electric-dipole ones. In this work, we propose analytical formulas for the estimation of the number of E2 lines in a transition array. Such expressions rely on statistical descriptions of electron states and J-levels. A generalized 'J-file' sum rule for E2 lines and the strength-weighted shift and variance of the line energies of a transition array nl^N+1 \rightarrow nl^Nn'l' of inter-configuration E2 lines are also presented.
Absorption and emission spectra of plasmas with multicharged-ions contain transition arrays with a huge number of coalescent electric-dipole (E1) lines, which are well suited for treatment by the unresolved transition array and derivative methods. But, some transition arrays show detailed features whose description requires diagonalization of the Hamiltonian matrix. We developed a hybrid opacity code, called SCORCG, which combines statistical approaches with fine-structure calculations consistently. Data required for the computation of detailed transition arrays (atomic configurations and atomic radial integrals) are calculated by the super-configuration code SCO (Super-Configuration Opacity), which provides an accurate description of the plasma screening effects on the wave-functions. Level energies as well as position and strength of spectral lines are computed by an adapted RCG routine of R. D. Cowan. The resulting code provides opacities for hot plasmas and can handle mid-Z elements. The code is also a powerful tool for the interpretation of recent laser and Z-pinch experimental spectra, as well as for validation of statistical methods.
The calculation of angular-momentum coupling transformation matrices can be very time consuming and alternative methods, even if they apply only in special cases, are helpful. We present a recursion relation for the calculation of particular 9j symbols used in the quantum theory of angular momentum.
In a recent note [High Energy Density Phys. 7, 161 (2011)], B.G. Wilson and V. Sonnad proposed a very useful closed form expression for the efficient generation of analytic log-linear radial meshes. The central point of the note is an implicit equation for the parameter h, involving Lambert's function W[x]. The authors mention that they are unaware of any direct proof of this equation (they obtained it by re-summing the Taylor expansion of h using high-order coefficients obtained by analytic differentiation of the implicit definition using symbolic manipulation). In the present comment, we present a direct proof of that equation.
We present a self-consistent approach to the modeling of dense plasma mixtures in local thermodynamic equilibrium. In each electron configuration the nucleus is totally screened by electrons in a Wigner-Seitz sphere (ion-sphere model). Bound and free electrons are treated quantum-mechanically. The assumption that all species should have the same electronic environment leads to the equality of the electronic pressure for all ions of all elements having therefore different cell volumes. The variation of the average atomic radii of the different elements with respect to temperature is investigated, and the procedure is applied to the determination of pressure in the Sun center.
The method proposed by Pratt to derive recursion relations for systems of degenerate fermions [Phys. Rev. Lett. 84, 4255 (2000), arXiv:nucl-th/9905055] relies on diagrammatic techniques. This efficient formalism assumes no explicit two-body interactions, makes possible the inclusion of conservation laws and requires low computational time. In this brief report, we show that such recursion relations can be obtained from generating functions, without any restriction as concerns the number of conservation laws (e.g. total energy or angular momentum).
In this short paper, the validity of Koopmans' theorem in the Hartree-Fock theory at non-zero temperature (Hartree-Fock statistical theory) is investigated. It is shown that Koopmans' theorem does not apply in the grand-canonical ensemble, due to a missing contribution to the energy proportional to the interaction between two electrons belonging to the same orbital. Hartree-Fock statistical theory has also been applied in the canonical ensemble [Blenski et al., Phys. Rev. E 55, R4889 (1997)] for the purpose of photo-absorption calculations. In that case, the Hartree-Fock self-consistent-field equations are derived in the super-configuration approximation. It is shown that Koopmans' theorem does not hold in the canonical ensemble, but that a restricted version of the theorem can be obtained, by assuming that a particular quantity multiplying the interaction matrix element in the expression of the energy does not change during the removal of an electron.
A study of the effects of non-extensivity on the modelling of atomic physics in hot dense plasmas is proposed within Tsallis' statistics. The electronic structure of the plasma is calculated through an average-atom model based on the minimization of the non-extensive free energy.
In this Brief Report, we propose general expressions for 6j and 9j coefficients arising in the study of the J-pairing interaction of identical fermions in a single-j shell. In the literature such coefficients are available only for a number of given angular momenta.
The impact of the third (skewness) and fourth (kurtosis) reduced centered moments on the statistical modeling of E1 lines in complex atomic spectra is investigated through the use of Gram-Charlier, Normal Inverse Gaussian and Generalized Gaussian distributions. It is shown that the modeling of unresolved transition arrays with non-Gaussian distributions may reveal more detailed structures, due essentially to the large value of the kurtosis. In the present work, focus is put essentially on the Generalized Gaussian, the power of the argument in the exponential being constrained by the kurtosis value. The relevance of the new statistical line distribution is checked by comparisons with smoothed detailed line-by-line calculations and through the analysis of 2p-3d transitions of recent laser or Z-pinch absorption measurements. The issue of calculating high-order moments is also discussed (Racah algebra, Jucys graphical method, semi-empirical approach ...).
Faussurier et al. [Phys. Rev. E 65, 016403 (2001)] proposed to use a variational principle relying on Jensen-Feynman (or Gibbs-Bogoliubov) inequality in order to optimize the accounting for two-particle interactions in the calculation of canonical partition functions. It consists in a decomposition into a reference electron system and a first-order correction. The procedure appears to be very efficient in order to evaluate the free energy and the orbital populations. In this work, we present numerical applications of the method and propose to extend it using a reference energy which includes the interaction between two electrons inside a given orbital. This is possible thanks to our efficient recursion relation for the calculation of partition functions. We also show that a linear reference energy, however, is usually sufficient to achieve a good precision and that the most promising way to improve the approach of Faussurier et al. is to apply Jensen's inequality to a more convenient convex function.
The exploration of atomic properties of strongly coupled partially degenerate plasmas, also referred to as warm dense matter, is important in astrophysics, since this thermodynamic regime is encountered for instance in Jovian planets' interior. One of the most important issues is the need for accurate equations of state and transport coefficients. The Ziman formula has been widely used for the computation of the static (DC) electrical resistivity. Usually, the calculations are based on the continuum wavefunctions computed in the temperature and density-dependent self-consistent potential of a fictive atom, representing the average ionization state of the plasma (average-atom model). We present calculations of the electrical resistivity of a plasma based on the superconfiguration (SC) formalism. In this modeling, the contributions of all the electronic configurations are taken into account. It is possible to obtain all the situations between the two limiting cases: detailed configurations (a super-orbital is a single orbital) and detailed ions (all orbitals are gathered in the same super-orbital). The ingredients necessary for the calculation are computed in a self-consistent manner for each SC, using a density-functional description of the electrons. Electron exchange-correlation is handled in the local-density approximation. The momentum transfer cross-sections are calculated by using the phase shifts of the continuum electron wavefunctions computed, in the potential of each SC, by the Schroedinger equation with relativistic corrections (Pauli approximation). Comparisons with experimental data are also presented.