The thermal conductivity of an arbitrarily dense plasma

-4L(o)s , - DISCLAIMER Thisreportwssprepsred aaanaccountof work sprsoredbyanagencyof theUnitedStatesGmxsrrnent NeithertheUnitedStatesGovernment noranyagencythereof,norarryof thetremplops, makesany warranty, expressorimplied,orassumes anylegalliabilityorrcsporrsiiility fortheaccuracy, compIetcnesa, orusefulness of anyinformation, apparatus, product,orprcassd&closed, orrepraanta thatitsw wordd notinfringeprivately ownedrights.Referencehereinto artyapeciflccommercial product,prociss,or servicebytradename,trademark, manufacturer, orotherwise, doesnotneccwarily mrrstitute orimplyits endorsement, rewmmendation, orfavoring bytheUnitedStataGovernment oranyagencythereof.The viewsandopintomof authorsexpressed hereindonotnecessarily stateorreflectthoseof theUnited StatesGovernment oranyagencythereof LA-10202-MS UC-34 Issued: Ohol~er 1984 The Thermal Conductivity of an ArbitrarilyDense Plasma - , , —.— - LOSES anrilos LosAlamosNationalLaboratory LosAlamos,NewMexico87545 THE THERMRL CONDUCTIVITY OF RN RRBXTRRRILY DENSE PLRSHR by George Rinker RBSTRRCT reports This is the second in a series of concerning the transport properties of dense plasmas In this aork we use the formalism of Lampe to extend our previous calculations of electrical conductivity to thermal conductivity and the calculation of thermoelectric coefficient Quantitative result are temperatures r nging from M -2 to gi~en for iron at -4 to 10sglcm3 10 eV and for densities from 3x1O Lampe for [19681 has calculated a weakly-coupled Chapman-Enskog Sonine carried is Ct31CUliltiOtl plasma with any by so[ving out method, using electrical and thermal degree of the Fermi-statistical transport electron degeneracy His equation by the Lenard-Ba\escu generalizations coefficients of the first two polynomials The physical electron-electron These are calculated for structure electron approximation approximate fluctuations coupling at high density report [Rinker of the Ziman and are scheme Our valid to The validity formula for much stronger of this for section a method structure avoid for the Born effects of complote using and couplings but of resistivity, arbitrary the shortcomings electron-ion is unknown use and neg\ect factors, mu(tip(e-scattering scheme potentials his the electrical thus logarithms Coulomb Debye-shielded cross and electron-ion from we described realistic compensate with by principally calculations for the expressed scattering 19841 iOl)iC potentials, degeneracy ultimately arises e(ectron-ion analysis for incorporates Born approximation the In a previous self-consistent are using of weak approximation partial-wave he interactions The requirement lattice models of Born Ue use an and density we have obtained good for agreement liquid generally as well only those by reproduce metals near as theoretical et al [19821 equivalent in the weak-coupling results combining approach, we logarithm with electrical Simply a numerical in that case consistent for set conduct,v,ty, these calculat~ons than prev~ous L&mpe’s all to have a ●,der the present highly are surpassed parametrized ●e for to expect the of val,d, ,on,z~t,on ty the formal stateO calculated for are J = electr~c current e = electron charge ~ = applied electric field def,ned ●,th respect to the transport the to be val~d an ,nternally electrical coeff,c,ent ,n temper~ture ow Coulomb expression se obta,n th,s In our Born approx~mathon Thus are explo,t wtth reproduce &nd thermoelectric range we the electron-,on h,s the equat,on Ue expect and dens,ty calculations coeff,c,ents to for coeff4cAents modify cond, t,ons ty work transport not choices calculations ad,us[ed for conduct,v, where IS as calculations thermal and Lenard-Balescu expression Iogarlthm, virtually of our calculations appropriate the thermal that Ue Coulomb are con~t,v,ty Lampe”s value conductivity electron-electron In improved replace m,th and limit electrical of Lampe to obtain that formula Ziman our models that calculations have shosn factor by In fact data structure fact the point pseudopotent,al pseudopotential known conduction-band Boercker the melting equ~t,ons p = ~ ‘ee ‘e ~ ‘ ‘T = electron ‘e ‘ pressure number density T = temperature heat (j= f{ux ~ = mean kinetic The electrical constraint conductivity S,12=S21 Wi(( Lampe’s reproduced generalized electron Uand thermal expressions here In the present for Chebyshev of extreme potential) The combinations of conductivity c (*ith to the the conventional Fulierton Each integral logarithms for the have not has because with has an asymptotic as where S they involve M expansion previously efficient effective is the expressions numerically factors of of be statistics highly are become as many not been available Lampe’s ultimately will electron obtained (WkT>lk14, fail terms and These approximations that include lengthy Coulomb [19821 degeneracy functions coefflc~ent rather account integrals approximations the are which these electron thermoelectric Sij approximations combinations be evaluated for addition application, in cases the integrals expressions 10-decimal be catted explicit Fermi-llirac Flccurate These per J=ldl are The quantity For energy except chemical conta,n unstable the integrals to the form m Ik(Z) &Zk ~diZ-i o i=O where z=u/kT i= the degeneracy parameter The instabilities arise through cancellation of the leading terms in z manipulations can be done analytically we have chosen and equivalent the easier approximations, combining numerically before f!s an computation for of the for 10-2 the ionization accommodate all arise ~S to 10sg/cm3 bound strongly range states are as kl shows shows the ionization state formed to increases apparent and Surface analysis partial-wave to these computed the readily the coupled various Figure and Fig ionization is in at Temperatures surface, thermal structure iron grid 3x10-4 increases p of the results previous 1-8 show numerically sufficient of we consider our IOU densities~ and show the electrical peak near and 3d states metal In other temperature Here regions no and point of compar our mode The states iC)niz~tiOn in the probably not probable errors point scattering, of and p physical of at comparisons value for arises from which kl, the conductivity which the S-l The of the 4s ionization iron makes behavior data our model electrical the of of 138.6 a exist is strictly resistivity the surprisingly a its transition electrical good the The inapplicability solid of of of phase the model of with other the model in the the agreement as experience factor in regions applicable good resutt This y significant, least U in the units normal son is gives kl d-nave for density conductivity at smatl experimental the experimental additional density more nearly is Virtually point normal and the strong conductivity melting coefficients from Figures from very The onset are asymptotic model sharp With because coefficients Figures clear kl-kl shell Instead Fullerton’s procedure on a logarithmic range fit from terms evaluating conductivity 19841 plot ionization irregularities as this as a three-dimensional electrons pressure transport [Rinker Zi as a contour of coefficients and densities state the and series in electrical and densities zero determining application and the surviving method of simply transport state to 104eV approaches present the the remaining temperatures same data sets and densities from of of ic)nizatit)n temperatures values instability example to extract and them, the expansions In principle, [iqUid The only at 118 ti(l~cm is of melting compared satisfying elements indicates Table at room temperature for the solid phase but I gives and at is the readily The apparent experimentally, calculated the resistivity This decreased considered decreases presumably arises S and show the transition-metal peak magnitude of because does not in dramatically from electrical conductivity, is negligible as the additional with temperature at whereas the are whereas temperature the is calculated values available for the cm-is-i scattering, Table conductivity good The decrease Table “ Comparison Of and experimental Phase calculated the the melting point but values transport monotonical decrease processes lye Experimental T (K) — bUi(son — I \ resistiv,ty for ~c (~(’)ocm) solid and liquid ,d? p (g/cm3) nc become phase electrical resistivity gives with experimental as additional II % phase near The in relative electron-electron is liquid units is smatler the solid lowered then rise in for agreement but first of electrical measurements the c - ‘%east is processes not transport conductivity process temperature not conductivity electrical to experimental as the measurements the contribute with deteriorates thermal the additional comparisons active, dependence here Figures which temperature n8 ~cdl iron ~d\ ~dl — Solid 7.86 293 13s 9.7a Solid 7.36 1810 126 127.Sb 1.lXM-S Liquid 7.0S 1810 118 138.6b l.Ox10-s 1.3X1O-S 6.SX10-3 ““0 2.4x10-4 [1983al [196S1 - — - Table Comparison of calculated conductivity Lx for solid conductivity and experimental aUeast thermal II cc(erg~K-ls-lcm-ll iron at zero pressure p (g/cm3) T (K] 7.86 300 0.ssx106 8.03x106 7,76 600 ~.f0xf06 S.47x106 7.66 900 1.68x106 3.80x106 7.S6 1200 2.26x106 2.82x106 7.46 1s00 2.86x106 3.18x106 C1983bl Figures and show the Figures and 10 show the conductive related simply to the thermal thermoelectric opacity conductivity Le in coefficient cm2g-1 Sj2 This by where ~sb/k It is = 4.106g6XIfj11 included here for S-l cm-2 convenient K-3 o comparison with radiative Opacities in cm-1s-1 quantity is Figures 11 electrical result and and are where assumed to 10-2, The ratit) temperature arises in both 1/2 low density from errors with cases the of in and Born of the Hubbard our work Plotted it)niz~tion states temperature and Lampe not between and comparison he regions In Hubbard varies a calculations conductivity to theirs, is set ~}S•P…»€• show their throughout The fact approximation it and the ratio of this where is not the ratio to be valid region uniformly differences Of our in Figs displayed calculation most that is and density consider and Lampe [19691 in high of apparently the structure factor Figures 13-18 electron-electron give compwisons contributions with the Lorentz to thermal conduction gas model and which yields neglects the simple relationships and Figures 13 and 14 show the ratAo of result Deviation ratio from electron COlliSiOIIS, density Figures In ~ddition regions 18 show Large ratio deviations change the above degeneracy the Figures from low temperature These region (low 8S12/3L0 occur for high 19 and 20 display a are thermal important nondegenerate is readi of at the the is to the Lorentz importance intermediate importance temperature which conductivity indicates 1S and 16 show the same for to of high in this the of very and high identically large electrcm- temperature thermoelectric and low coefficient deviations occur Figures density) gas in the Lorentz 17 in and gas model degeneracy the degeneracy system parameter to a degenerate system u/kT at high The dramatic density and ly apparent o c) c c: - m -, w @4 u w .> al c d - L 27 28 - “; d c - d - mu c OL - a l-l IA 907 m m n u 29 30 q: $ w-i o z + 6-E I O“f I 1“2 2“1 E“o 9“o- I y-i s“I - 08 n 37 L“ * % w c E UJ c - - 32 3MlLJU3dW31 307 0) c - c -< I c Ku N OO?? @ n 33 x m $ \ , , , ‘% , , , , , tgi”c - m = 6-C , 0“s 0“s I“z 1“2 Z“l 3M-l&w31 !307 S“o E“o 9“o- 9“o- s“l - 0 E LL A\@ fn Y al g (J n u- 35 e m o $ ,, ,, ‘, $., ~ \ , , , , , , ,., , 0“s ‘ O“f ,., ,t ,$ ‘“ ‘ $ ,’ \ @“c m m G A 6“E z + , ,., \ 1“2 ,, , , 1“2 ,, , 2“1 , Z“I ,* ,,t S“o ,, S“o , 9“o- 9“o- ‘ ” - z 07 El n c1 -1 cm o s“I - i C u Is a N L&N 03”0 O“K” LWIMI 907 z’ 00z’” fill @ I@ n L- n 1! w ml f- 37 * w L-q ,.- , , t , ‘, , 0“s ‘, “ ,.,, ,, ‘, *, ‘., ,, ,,*“.‘ “, ,’ gi”s $ m m , - , ‘ , ‘, ‘, t ‘ ‘, ‘ ‘,“ , *, “ ., ‘ “, “.,, ‘, ‘., I “z ; , \ \, ‘., N \ , \ , 1“z z“ t ’ , ‘“ ‘, ‘ z“ I - \ 3M1.NJt13dw31 907 \ tt ,, , ,, “, ‘, ‘, “ ., t -; \\\\ ‘ ~ ‘“ / “’., ‘, ‘“’., “’ , ‘, “ ‘., ‘, “., ‘, ‘ “ ,‘, , ‘.,,’,,’ ~ ‘, “., , %,‘.,\‘“,‘$,““ ,““ ‘ ., “,’\, ’ , ,‘., “.,,’,., , \ /~ , “, ‘., , ‘, “,, “, ‘, , ““.,’ $,’.,““q $,”,, , * O“i \ ,, ‘, “’., ;$,:., ““.,:(’ , j’, ~ 0“ ,,,’,, ”,, ‘, “ ,,.,, ,, , ‘, \‘ , , , ~~~ t, ,, “.,’’,’”” ,‘., C)” J -., ,, “6”i I f.”o I / I o “ -’* \g I 9“o- , t 9“o- t , # ; , ,’ .- $ m m G + -1 o (SY l-l IL u 03 ml LLl dL o c 1- S“1- g ‘ ‘\ ‘ S 1’- ‘ -,, ,, 1’ ‘ ,, ?= ’”-” , ’ I !4 + 6z , C.o ,,, ,, - -— ,,,, ’”’’”””” I ’’’””” ” , , ~ ““”:o 7“”” ””” / -I , J,,, ,, -,, ~~~‡“’ ,, ,.” ””” ’ ” ‘“ , , .’ ’” ,,, ““””’‘“””’/ ,, ,’ , ,,,.;, ,:,: ;,::;;+ :+:+::= ,’ ’ J , , I ,,t,,, :1 I , ,, , , , ”,, , ;,: :, t ,: ;,, :, , , , ’ , , : t: ,: * t , ’ ; :; ,:: , ,$,,I, ,: 1’ , :;;;:;;;; ’ ,“ #, ;,,,.”/, ,(;.”;,1,’,’ ,1 ~: z ,%”, ”,;;,;;;;;:,:,::;, ‘ - ,.- - y;,,,, , , -’-””/ / -1 s“ ‘“;;:(W , ”,, ””,, ”” ,-/ - ’””””” ~,,,.g,~@,;::::::::-: - ‘E=ii=”””” / -./ -: - u-i , l - / -0 , ,’ ,’ ,.’ - ~, *I, I8 ,, ;,t: ’ ’ , ,: , ,.I , ,,, ~.- /’,,, , ,, , ’ ‘1 ., ,, ,, J,,, ’ , ,’/; :/ ,,, ,, ’, - .,,’ ’ ~- 8:.,, ~~" ““;”” ”’””/”;/:? m y- I i i 3MlllM3dW31 907 ●U.S GOVERNMENT PRINTING OFFICE: 1984-O-576-034/4213 — ,, .,- , ., , -—- n-.- ~-.-: - =:2 + ., + ,.,-.:!., :+.= :=’ ,- - .-’ ‘ , - ’ .- ,., -w ;,::.’::lp!#!g!! > - -, ,U A.= .-, *.?l, n”l , 8— - : — .: , , , ! , — — ““ —.- l$%zamwMrfii ,.$ ,- — — : ,+ ,.- :x., & -: : - - ., ;~ -, - - -+ ’ “l.—— , ,,-.,.-, >\ - , ,- , - —., -“ - y&JL-L;i , -., ,-, — y-u 5==s r , ~cm!!! ii— -“ -J!!!&& - —.4 , .- -— —.+ .- .> -J - -,.-— - : ,+ r,y~- ,:,,,:,.,: , a ,’R\-F +- !, :;- -.;”” -“-”.:” “, .: - “:!’ T- , ,.;, i~~,-~ ,,~=aai ’: Mmisi%az’al K3W.-:< ‘A: : .,: -, ~ ; ( : ., — —: — - ,? - ,=- ,.-: T -— -L — - .— -= z —— -4’”:_ , ,,.> - - -7 : ? ‘awl “’-.,;,; ... aaanaccountof work sprsoredbyanagencyof theUnitedStatesGmxsrrnent NeithertheUnitedStatesGovernment noranyagencythereof,norarryof thetremplops, makesany warranty, expressorimplied,orassumes anylegalliabilityorrcsporrsiiility... stateorreflectthoseof theUnited StatesGovernment oranyagencythereof LA-10202-MS UC-34 Issued: Ohol~er 1984 The Thermal Conductivity of an ArbitrarilyDense Plasma - , , —.— - LOSES anrilos LosAlamosNationalLaboratory... exist is strictly resistivity the surprisingly a its transition electrical good the The inapplicability solid of of of phase the model of with other the model in the the agreement as experience