Hostname: page-component-74656f9884-c52kr Total loading time: 0 Render date: 2026-11-06T23:50:58.705Z Has data issue: false hasContentIssue false
  • English
  • Français

Low velocity ion slowing down in a de-mixing binary ionic mixture

Published online by Cambridge University Press:  15 March 2011

C. Deutsch ,
D. Leger  and
B. Tashev
C. Deutsch*
Affiliation:
Université Paris-Sud, LPGP (UMR-CNRS 8578), Orsay, France
D. Leger
Affiliation:
Université Valenciennes-Hainaut Cambresis, Lab., Monthouy, France
B. Tashev
Affiliation:
Kazakh National University, Department of Physics, Almaty, Kazakhstan
*
Address correspondence and reprint requests to: C. Deutsch, Université Paris-Sud, LPGP (UMR-CNRS 8578)Bât. 210, F-91405 Orsay, France. E-mail: claude.deutsch@pgp.u-psud.fr

Abstract

We consider ion projectile slowing down at low velocity Vp < Vthe, target thermal electron velocity, in a strongly coupled and de-mixing H-He ionic mixture. It is investigated in terms of quasi-static and critical charge-charge structure factors. Non-polarizable as well as polarizable partially degenerate electron backgrounds are given attention. The low velocity ion slowing down turns negative in the presence of long wavelength and low frequency hydromodes, signaling a critical de-mixing. This process documents an energy transfer from target ion plasma to the incoming ion projectile.

Keywords

Critical de-mixingBinary ionic mixturesLow velocity ion stopping

Information

Type
Research Article
Information
Laser and Particle Beams , Volume 29 , Issue 1 , March 2011 , pp. 121 - 124
Copyright
Copyright © Cambridge University Press 2011

INTRODUCTION

A sustained and widespread interest is currently documented for low velocity ion slowing down (LVISD) in binary ionic mixtures (BIM) (Tashev et al., Reference Tashev, Baimbetov, Deutsch and Fromy2008; Fromy et al., Reference Fromy, Tashev and Deutsch2010) fully neutralized by an electron fluid jellium (Leger & Deutsch, 1988a, 1988b, Reference Leger and Deutsch1992).

We pay specific attention to the hydrogen-helium strongly coupled mixture of sustained astrophysical concern (Stevenson, Reference Stevenson1975; Stevenson & Salpeter, Reference Stevenson and Salpeter1977; Vorberger et al., Reference Vorberger, Tamblyn, Militzer and Bonev2007; Lorenzen et al., Reference Lorenzen, Holst and Redmer2009). Up to now, most of these studies were conducted through, a weakly coupled BIM immersed into a classical high temperature electron background taken in a dielectric Fried-Conte description (Fried & Conte, Reference Fried and Conte1961; Arista & Brandt, Reference Arista and Brandt1981). Specifying BIM ion coupling on species i (i = 1, 2) with charge Zi one has

(1)
\Gamma _i={{Z_i^2 e^2 } \over {k_B Ta_i }},

where a i = ((4π/3) ni)−1/3 with temperature T = T e = T i. Here, we stress strongly coupled BIM with Γi ≫ 1, able to display a critical de-mixing process (2, 4).

Pertaining thermodynamic and hydrodynamic features have already been investigated at length (Leger & Deutsch, 1988a, 1988b, Reference Leger and Deutsch1992). In this context, a critical de-mixing behavior is observed on the ion-ion structure factor S αβ(q) in the long-wavelength q → 0 limit. As usual, LVISD amounts to evaluating low velocity ion stopping with neglected intra-beam scattering and in-flight correlation, in a low velocity regime V p ≪ V the target electron thermal velocity. On the other hand, the opposite V p > V the situation is now rather well documented in the so-called standard stopping model (Deutsch, Reference Deutsch1986; Deutsch et al., Reference Deutsch, Maynard, Bimbot, Gardes, DellaNegra, Dumail, Kubica, Richard, Rivet, Servajean, Fleurier, Sanba, Hoffmann, Weyrich and Wahl1989) where most of the ion projectile stopping arises from its interaction with target electron fluid. To the contrary, LVISD is essentially monitored by the projectile interaction with target ions. A foremost motivation for the present undertaking is a possible involment of LVISD in BIM in probing and testing the ion de-mixing process. In the past, some of us (Leger & Deutsch, 1988a, 1988b, Reference Leger and Deutsch1992) have already investigated a critical BIM de-mixing signature on target electrical resistivity. Within a dielectric framework for target particles, the nonrelativistic ion stopping thus reads as (Arista & Brandt, Reference Arista and Brandt1981).

(2)
\eqalign{S=-{{dE} \over {dx}}= & - {1 \over {V_p }}\left({{{dE} \over {dt}}} \right)_0 \cr & ={2 \over {\rm \pi} }\left({{{Z_p e} \over {V_p }}} \right)^2 \int_{0}^{\infty} {{{dq} \over q}\int_{o}^{qV_p} {d{\rm \omega} {\rm \omega} {\mathop{\rm Im}} \left({ - {1 \over {\varepsilon \lpar q\comma \; {\rm \omega} \rpar }}} \right)} }\comma}

which can be straightforwardly re-expressed in terms of the ion charge-ion charge structure factor when switching to very low

V_p \leq {\overline {V_{thi}}} =C_1 V_{th1}+C_2 V_{th2}\comma \;

with V thi, thermal velocity and C i, relative concentration of ion i.

Incoming ion projectile could then be able to probe every available ion fluid fluctuations in target by restricting the global dielectrics expression ɛ(q,ω) to its i-component. Such a picture highlights the BIM electron background following the ionic fluctuations, within a polarization concept.

Energy exchange between ion projectile and thermalized medium is expressed in terms of emission and absorption processes, according to the protocol (Pines, Reference Pines1964)

\eqalign{\left({{{dE} \over {dt}}} \right)_0 & ={\displaystyle \int}_{{\rm \omega}\gt 0} {d^3 \vec qN \lpar {\rm \omega} \rpar f \lpar \vec q\comma \; {\rm \omega} \rpar - } {\displaystyle \vint}_{{\rm \omega}\gt 0} {d^3 \vec q\lpar N\lpar {\rm \omega} \rpar +1\rpar f\lpar \vec q\comma \; {\rm \omega} \rpar }}

yielding, up to the first order in the interaction, a quarter involving only spontaneous emission. The present approach stands at variance with the so-called T-matrix one (Gericke et al., Reference Gericke, Schlanges and Kraeft1996; Gericke & Schlanges, Reference Gericke and Schlanges1999) advocating a strong binary interaction between incoming ion projectile and one target electron thus probing mostly short wavelength modes with q → ∞. Then,

(3)
S_{zz} \lpar \vec q\comma \; {\rm \omega} \rpar = {{\hbar q^2 } \over {4{\rm \pi} ^2 e^2 }}N\lpar {\rm \omega} \rpar {\mathop{\rm Im}\nolimits} \left({ - {1 \over {\rm \varepsilon \lpar q\comma \; {\rm \omega} \rpar }}} \right)\comma

where N(ω)−1 = eβħω − 1, β = (kBT)−1 allow to put Eq. (2) under the form

(4)
\eqalign{S & ={2 \over {\rm \pi} }\left({{{Z_p e} \over {V_p }}} \right)^2 {\vint}_0^{\,\infty} {{{dq} \over q}{\vint}_o^{\,qV_p} {d{\rm \omega} {{4{\rm \pi} ^2 e^2 {\rm \omega} } \over {\hbar q^2 }}} } \cr & \times S_{zz} \left({\vec q\comma \; {\rm \omega} } \right)\left({e^{{\rm \beta} {\hbar} {\rm \omega} } - 1} \right)\comma} \;

with the usual electron ɛ(q,ω) now extended to the ion components building up the BIM, and the charge-charge structure factor

(5)
S_{zz} \lpar q\rpar =\sum\limits_{{\rm \alpha}\comma {\rm \beta}=1}^2 {} \lpar c_{\rm \alpha} c_{\rm \beta} \rpar ^{1/2} Z_{\rm \alpha} Z_{\rm \beta} S_{{\rm \alpha} {\rm \beta} } \lpar q\rpar \comma \;

where C1 + C2 = 1. Ci refers to concentration of species i within a BIM built on target ion charges Z1 and Z2.

Focusing attention on the slow and long wavelength hydromodes ω → 0, q → 0 monitoring BIM de-mixing, one can safely restrict Eq. (4) to its static limit ω → 0

(6)
S = {8{\rm \pi} \over 3} \lpar Z_{\,p}e^{2}\rpar ^{2} {\rm \beta} V_p {\vint}_{0}^{\infty} dq S_{zz} \lpar q\rpar

in-terms of $S_{zz} \lpar q\rpar ={\vint}_0^{\infty} {_{}^{} d{\rm \omega} S_{zz} } \lpar q\comma \; {\rm \omega} \rpar $ Eq. (6) also highlights the expected LVISD linear Vp − dependence ≤  $\overline {V_{thi}}$ Eq. (6) implies an average over every ω-fluctuation, available to $S_{zz} \lpar \vec q \comma \; {\rm \omega}\rpar$.

CRITICAL Szz(q)

Mean field classical description of BIM de-mixing could be rather straightforwardly explained with the static charge-charge structure factor (De Gennes & Friedel, Reference De Gennes and Friedel1958; Fisher & Langer, Reference Fisher and Langer1968),

(7)
S_{ZZ}={{\Sigma \left\vert t \right\vert ^{ - \rm \gamma } \lpar qa_i \rpar ^2 } \over {3\bar \Gamma \bar Z^2 \left({\lpar q{\rm\xi} \rpar ^2+1} \right)}} \comma

where t = (T − Tc/Tc), Tc = critical temperature, γ = 1, and ξ = ion-ion correlation length featuring $\mathop{\lim }\nolimits_{\left\vert t \right\vert \to 0} {\rm\xi} \to \infty$. Z̄ = C1Z1 + C2Z2 and Γ̄ = C1Γ1 + C2Γ2. ∑ denotes a constant normalizing factor accessed through the sum rule (q in āi−1)(where āi = C1a1 + C2a2)

(8)
{\int}_{0}^{\infty} {_{}^{} dqq^2 \left[{S_{zz} \lpar q\rpar - \overline {z^2 } } \right]}=- {{3{\rm \pi} } \over 2}\comma \;

where $\overline{z^{2}} = 1 + {\rm C}_1 {\rm C}_2\lpar Z_1 - Z_2\rpar ^2 / {\bar z}^2$.

Eq. (7) mostly emphasizes long distance hydromodes, of significance at critical de-mixing.

Paying attention first to non-polarizable BIM with a fixed and rigid electron background (Leger & Deutsch, 1988a, 1988b), the correlation length reads as

(9)
\left({{{\bar a_i } \over {\rm \xi} }} \right)^2=3\bar \Gamma\, \overline{z}^2 {{D_I } \over {D_R }}\comma \;

where D I and DR, respectively, denote the q → 0 limit of

(10)
\eqalign{D_1 \lpar q\rpar =\bar Z^2 - C_1 C_2 \left[{Z_1^2 {\mathop{C}\limits^{ \frown}}_{22}^{ R} \lpar q\rpar +Z_2^2 {\mathop{C}\limits^{ \frown}}_{11}^{ R} \lpar q\rpar - 2Z_1 Z_2 {\mathop{C}\limits^{ \frown}}_{22}^{ R} \lpar q\rpar } \right]\comma \;}

and

(11)
D_R \lpar q\rpar =1 - c_1 {\mathop{C}\limits^{ \frown}}_{11}^{ R} \lpar q\rpar - c_2 {\mathop{C}\limits^{\frown}}_{22}^{ R} \lpar q\rpar - c_1 c_2 \det \left\vert {{\mathop{C}\limits^{\frown}}_{{\rm \alpha} {\rm \beta}}^{ R} \lpar q\rpar } \right\vert \comma \;

in terms of

(12)
{\mathop{C}\limits^{ \frown}}_{{\rm \alpha} {\rm \beta}}^{ R} \lpar q\rpar ={\mathop{C}\limits^{ \frown}}_{{\rm \alpha} {\rm \beta}}^{ R} \lpar q\rpar +Z_{\rm \alpha} Z_{\rm \beta} \mathop{\nu}^{\frown} \lpar q\rpar \comma \;

${\mathop{C} \limits^{ \frown}}_{{\rm \alpha} {\rm \beta}}^{ R}\lpar q\rpar$ denotes the partial direct correlation function, viewed in the Ornstein-Zernikc (OZ) equations

(13)
\mathop{h}^{\frown}{\!}_{{\rm \alpha} {\rm \beta} } \lpar k\rpar =\mathop{C}^{\frown}{\!}_{{\rm \alpha} {\rm \beta} } \lpar q\rpar +\sum\limits_{\nu=1}^2 {c_\nu } \mathop{h}^{\frown}{\!} _{a\nu } \lpar q\rpar \mathop{C}^{ \frown}{\!} _{\nu {\rm \beta} } \lpar q\rpar . \;

The right-hand side of Eq. (12) also features the dimensionless and screened Coulomb potential

(14)
\mathop{\nu}^{\frown} \lpar q\rpar ={{4{\rm \pi} {\rm \beta} e^2 \bar n_i } \over {k^2 {\rm\varepsilon} \lpar q\rpar }} \equiv {{3\bar \Gamma } \over {\bar a_i^2 }}{1 \over {k^2 {\rm\varepsilon} \lpar q\rpar }}\comma \;

with the dimensionless static electron fluid dielectric function ɛ(q). Close to criticality, one expects a characteristic diverging behavior of the correlation length, so that

(15)
{\rm\xi}={\rm\xi} _0^ \pm \left\vert t \right\vert ^{ - \nu }\comma \; t \to 0

with ν = 0.5, in a standard mean field OZ approximation.

SUPERELASTIC LVISD

The introduction of Eq. (7) into Eq. (6), thus yields for q āi ≤ 1

(16)
\eqalign{S & ={{8{\rm \pi} } \over 3}\left({Z_p e^2 } \right)^2 .{{{\rm \beta} V_p } \over {\bar a_i }}.{{\Sigma \left\vert t \right\vert ^{ - 1} } \over {\bar \Gamma \bar Z^2 }}\left({{\xi \over {\bar a_i }}} \right)^2 \left[{1 - Tan^{ - 1} \left({{\xi \over {\bar a_i }}} \right)} \right]\comma \; }

demonstrating that ξ-diverging behavior (Eq. (15)) is compensated by that in Eq. (7).

Now we pay attention to correlation length (Eq. (9)) estimates by solving simultaneously (Leger & Deutsch, 1988a, 1988b) OZ Eq. (13) with the Hypernetted Chain (HNC) equations

(17)
g_{{\rm \alpha} {\rm \beta} } \lpar r\rpar =\exp \left[{h_{{\rm \alpha} {\rm \beta} } \lpar r\rpar - \lpar C_{{\rm \alpha} {\rm \beta} }^R \lpar r\rpar } \right]\comma \;

valid for any Γi values.

Lindhard screening (Lindhard, Reference Lindhard1954) involves (cf. Fig. 1)

(18)
{\rm\varepsilon} \lpar k\rpar ={{k^2+k_{TF}^2\, g_L \lpar x\rpar } \over {k^2 }}\comma \;

where kTF = (6πnee2/EF)1/2 denotes the Thomas-Fermi wave vector and the function gL(x) depends only on the dimensionless variable x = k/2kF, with

(19)
g_L \lpar x\rpar ={1 \over 2}+{{1 - x^2 } \over {4x}}\ell n\left\vert {{{1+x} \over {1 - x}}} \right\vert \comma \;

and kF = (3π2ne)1/3, Fermi wave number in terms of electron density ne, while EF = ħ F2kF2/2me. Moreover, Hubbard screening (Hubbard, Reference Hubbard1966) with

(20)
g\lpar x\rpar ={{g_L \lpar x\rpar } \over {1 - {{k_{TF}^2 } \over {k^2 }}g_L \lpar x\rpar G\lpar x\rpar }}\comma \;

where $G\lpar x\rpar ={{x^2 } \over {2x^2+1/2}}\comma$ includes exchange-correlation contributions into the jellium background.

Fig. 1.

Plot of the reduced squared correlation length (ξ/a1)2 in terms of r s along a critical and vertical line (34% He, Γ = 60) with Lindhard screening. There ξ2 remains negative.

Present critical conditions are now significantly different. In Figure 2, critical de-mixing occurs at C2 = 0.75 in lieu of C2 = 0.34 in Figure 1 for the same Γ = 60. More importantly, ξ2 is now allowed to increase positive on the largest rs range side.

Fig. 2.

PBIM model B with Hubbard screening. Plot of the reduced squared correlation length (ξ/a1)2 as a function of rs along a critical and vertical line (75% He at Γ = 60) ξ2 is now allowed to change sign for rs ≥ rs 0.

The resulting correlation length is now rising strongly when ∣t∣→ 0, up to the standard means field behavior (ξ/a1)2 ≫ 1.

The given LVISD is now negative featuring a super elastic interaction between the low velocity incoming ion projectile and the PBIMB target.

It can also be appreciated that before turning negative, LVISD (16) vanishes for (ξ/a1)2 ~ 1.65. The prefactor ∑ in Eq. (16) is then straightforwardly derived from the sum rule (8) under the alternative forms,

(21)
\Sigma=\left({ - {{3{\rm \pi} } \over 2}+\overline {z^2 } } \right)^3 {{D_R \left\vert t \right\vert } \over {D_I }}

in terms of BIM thermodynamics and also (cf. Eq. (15))

(22)
\Sigma=q\bar \Gamma ^2 \overline {z^2 } \left({ - {{3{\rm \pi} } \over 2}+\overline {z^2 } } \right)^3 \left({{{\xi _0^ \pm } \over {\bar a_i }}} \right)

FINAL REMARKS

The present developments highlight for the first time the intriguing interplay of a first order de-mixing process in a strongly coupled and binary ionic mixture with a low velocity incoming ion beam. The latter may be envisioned for diagnostics purposes or target conditioning in the subfields of ICF and warm dense matter, for instance.

Within a fundamental statistical physics perspective, it should be appreciated that the above results document unambiguously the potentialities of probing collective very long wavelength phenomena occurring in a plasma target with low velocity ion beams via the evaluation of a transport coefficient, featured in the present context by a stopping power mechanism.

As far as a super elastic projectile-target interaction is witnessed at very low V p, it should be noticed that it also can appear in the T-matrix approach (Gericke et al., 1996).

ACKNOWLEDGMENT

We thank the referee for very significant and constructive remarks.

References

REFERENCES

Arista, N.R. & Brandt, W. (1981). Energy-loss and straggling of charged-particles in plasmas of all degeneracies. Phys. Rev. A 23, 1898–1905.Google Scholar
De Gennes, P.G. & Friedel, J. (1958). Anomalie de resistivité dans certains metaux magnétiques (Resistivity anomalies in magnetic metals). J. Phys. Chem. Solids 4, 71–77.Google Scholar
Deutsch, C. & Leger, D. (1993). Electron transport with ion demixing. Contrib. Plasma Phys. 33, 409–420.CrossRefGoogle Scholar
Deutsch, C. (1986). Inertial confinement fusion driven by intense ion beams. Ann. Phys. Paris 11, 1–111.Google Scholar
Deutsch, C., Maynard, G., Bimbot, R., Gardes, D., DellaNegra, S.,Dumail, M., Kubica, B., Richard, A., Rivet, M.F., Servajean, A., Fleurier, C., Sanba, A., Hoffmann, D.H.H., Weyrich, K. & Wahl, H. (1989). Ion stopping in dense plasmas. A standard model approach. Nucl. Instrum. Meth. A 278, 38–43.CrossRefGoogle Scholar
Fisher, M.E. & Langer, J.S. (1968). Resistive anomalies at magnetic critical points. Phys. Rev. Lett. 20, 665–668.CrossRefGoogle Scholar
Fromy, P., Tashev, B. & Deutsch, C. (2010). Very low velocity ion slowing down in binary ionic mixtures: Charge-and mass-asymmetry effects. Phys. Rev. Special Topics Accelerators and Beams 13, 1013C2–8.Google Scholar
Fried, B.D. & Conte, S.D. (1961). The Plasma Dispersion Function. New York: Academic Press.Google Scholar
Gericke, D.O., Schlanges, M. & Kraeft, W.D. (1996). Stopping power of a quantum plasma T-matrix approximation and dynamical screening. Phys. Lett. A 222, 241–245.Google Scholar
Gericke, D.O. and Schlanges, M. (1999). Beam-plasma coupling effect on the stopping power of dense plasmas. Phys. Rev. E 60, 904–909.CrossRefGoogle ScholarPubMed
Hubbard, W.B. (1966). Studies in stellar evolution. V. Transport coefficients of degenerate stellar matter. Astrophys. J. 146, 858–869.Google Scholar
Leger, D. & Deutsch, C. (1988). Critical demixing in strongly coupled ionic mixtures . I. Weak electron-ion coupling and stability requirements. Phys. Rev. A 37, 4916–4929.CrossRefGoogle ScholarPubMed
Leger, D. & Deutsch, C. (1988). Critical demixing in strongly coupled ionic mixtures. II. Hydrodynamic correlation functions. Phys. Rev. A 37, 4930–4942.CrossRefGoogle ScholarPubMed
Leger, D. & Deutsch, C. (1992). Linear and electronic transport in strongly coupled ionic H-He mixtures. Phys. Fluids B 4, 3162–3184.CrossRefGoogle Scholar
Lindhard, J. (1954). Electron dielectric function. Mat. Fys. Medd. K. Dan. Vidensk. Selsk. 28, 1.Google Scholar
Pines, D. (1964). Elementary Excitations in Solids. New York: Benjamin.Google Scholar
Lorenzen, W., Holst, B. & Redmer, R. (2009). Demixing of hydrogen and helium at Megabars pressure. Phys. Rev. Lett. 102, 115701/4.CrossRefGoogle Scholar
Stevenson, D.J. (1975). The thermodynamics and phase separation of dense fully ionized H-He fluid mixtures. Phys. Rev. B 12, 3999–4007.Google Scholar
Stevenson, D.J. & Salpeter, E.E. (1977). Phase-diagram and transport properties for H-He fluid planets. Astrophys. J. Suppl. 35, 221–235.CrossRefGoogle Scholar
Tashev, B., Baimbetov, F., Deutsch, C. & Fromy, P. (2008). Low ion velocity slowing down in dense binary ionic mixtures. Phys. Plasmas 15, 10270.CrossRefGoogle Scholar
Vorberger, J., Tamblyn, I., Militzer, B. & Bonev, S.A. (2007). H-He mixtures in the interior of giant planets. Phys. Rev. B 75, 024206.Google Scholar
Figure 0

Fig. 1. Plot of the reduced squared correlation length (ξ/a1)2 in terms of rs along a critical and vertical line (34% He, Γ = 60) with Lindhard screening. There ξ2 remains negative.

Figure 1

Fig. 2. PBIM model B with Hubbard screening. Plot of the reduced squared correlation length (ξ/a1)2 as a function of rs along a critical and vertical line (75% He at Γ = 60) ξ2 is now allowed to change sign for rs ≥ rs0.