Influence of shell structure on the equation of state of matter at high pressures* E . BOSCHI ( * * ) Received on February 21st, 1972 SUMMARY. — W e develop an equation of state of matter in the pres- sure range of geophysical interest introducing shell structure in the Thomas- F e r m i - D i r a c atom model. RIASSUNTO. — Viene sviluppata un'equazione di stato della materia nell'intervallo di pressioni d'interesse geofisico, introducendo la struttura a strati nel modello atomico di Tliomas-Fermi-Dirac. The Thomas-Fermi model, in its original formulation, cannot give account for the existence of discrete energy levels for bound electrons or, in other words, is not sensitive to the effects of shell structure (2). I t is also obvious that these effects are the principal cause of the in- applicability of the Thomas-Fermi-Dirac model in the pressure range of geophysical interest (4). Up today, these effects have been considered only from a qualitative point of view and the influence of them on the equation of state has not been determined. W e believe that it is possible to solve consistently this problem by means of an appropriate modification of the normalization of the Thomas-Fermi-Dirac elec- tronic charged cloud. As it is well-known, in the Thomas-Fermi- Dirac theory, the number of the possible states per volume unit ( * ) P a p e r presented at the " F i r s t European Eartli and P l a n e t a r y Physics C o l l o q u i u m " , Reading ( U . K . ) , 30 March-2 A p r i l 1971. ( * * ) Istituto di Fisica, Università degli Studi di Bologna ( I t a l y ) . 3 3 4 E . B O S C H I of the electrons, whose momentum is between p and p + dp, is 8 n p2 dp¡h3. Their density is given by the Fermi-Dirac statistics: , 8 np2 dp dp = . , Ekin — ew 7i exp I IT + f l where Eun is tlie kinetic energy, cp the Thomas-Fermi potential, k the Boltzman constant, T the temperature, and /t the chemical potential. If Ave consider an atom of atomic number Z, the integral of do over the entire phase-space atomic volume must satisfy the con- dition: 3 2 N 2 h3 r2 dr p2 dp exp E tin — e cp kT Z Furthermore the Poisson's equation gives another relation between the electron density and the potential. Now we assume that only a part of the electrons in tlie atom can satisfy to the Thomas-Fermi-Dirac requirements. The assumption means that there is, for a given value of the pressure and of the tem- perature, some "critical radius", Bc, which divides the atom volume into two parts: an outer part containing electrons which can be ac- tually considered as a charge cloud and an inner part containing electrons whose behaviour needs a quantum mechanical description. So, for given external conditions, we have B " b o u n d " electrons and F " f r e e " electrons, where B and F are linked by the obvious relation: F + B = Z. I n more precise words we introduce a cutoff in the mo- mentum phase space of the atom that means to change the lower limit of momentum integration, 0, which appears in tlie normalization equation, to some value p. I n this manner we obtain the actual number of electrons, the so-called free electron, which can be treated statistically. Since the free electrons are contained in a continuum of energy levels, the use of the Thomas-Fermi-Dirac theory appears well justified (3). The way to determine the counting of bound electrons is to use a trial potential V in Sclirodinger's equation, to obtain the energy I N F L U E N C E OF S H E L L S T R U C T U R E ON T I I E E Q U A T I O N E T C . 3 3 5 levels corresponding to y , and calculate tlie number of bound electrons from tlie relations: exp (— Eni/kT + n) + 1 and B = Sn, where n and I are tlie principal and orbital quantum numbers, N„i is the average number of electrons in the subshell (nl) and B is the total number of bound electrons. After these considerations, the Thomas-Fermi-Dirac normalization condition becomes: 32 n- ? F r2 dr p2 dp h3 / / (E*in — eV ~ 6XP I îcT +-") + 1 0 Rc + V 2(21 + 1) = z / ni exp (—Eni/kT + n) + 1 Searching the appropriate potential V which appears in the last equation, it seems natural to take into consideration the ionization phenomena and to represent the passage from bound electrons to free electrons as some kind of ionization induced by pressure. This makes possible the task of finding an appropriate potential function 0): V = R \r 2 r2 2/ ^ F e f r 2 A X 0 r (1 + Ar)2 R \2 R where Ax is chosen to make V continuous at r — Rc: _ Z R 5 , 3 ' 1 ~ Rc (1 + A Rc)2 F Rc 2 " By means of Poisson's equation: 2 Z A Rc Z (1 + A Rc)3 (1 + A Rc)2 8 E . BOSCHI where A is a constant for a given value of Z. Bc divides the atomic volume into two regions: an inner region around the nucleus where the electron density is high and due to bound electrons, and an outer region where the electron density is low and due to free electrons. To perform the calculations, we use an iterative process starting from a trial value of the electron density. The figure shows the results obtained for an atom of atomic number 26 and puts in evidence the discrepancy between the Thomas-Fermi-Dirac model, represented by I N F L U E N C E OF S H E L L S T R U C T U R E ON T H E E Q U A T I O N E T C . 3 3 7 a continuous curve line, and the results of these considerations, re- presented by a broken and dotted line. The dotted line, near to the continuous one, represents a Thomas-Fermi-Dirac calculation with the potential V made to check the internal consistency of the work. ACKNOWLEDGEMENTS I wish to thank Prof. M. Caputo for valuable discussions on this work. R E F E R E N C E S T1) A R M S T R O N G B . H . , J O H N S T O N R . R . , K E L L Y P . S . , 19(37. - Progress in high temperature physics and chemistry, edited b y Carl. A . Rouse, P e r g a m o n Press. ( 2 ) BOSCHI E., CAPUTO M., 1969. - « N u o v o C i m e n t o » , 1, 441. ( 3 ) BOSCHI E., CAPUTO M., 1971. - Mantle and Core in Planetary Physics, edited by J. Coulomb and M. Caputo, Academic Press, Milano. C ) T A K E U C H I I L , K A N A M O R I H . , 1 9 6 6 . - « J o u r n . G - e o p h y s . R e s . » , 7 1 , 3 9 8 5 .