Abstract
A key component for the description of charged particle systems is the screening of the Coulomb interaction between charge carriers. First investigated in the 1920s by Debye and Hückel for electrolytes, charge screening is important for determining the structural and transport properties of matter as diverse as astrophysical and laboratory plasmas, nuclear matter such as quark-gluon plasmas, electrons in solids, planetary cores and charged macromolecules. For systems with negligible dynamics, screening is still mostly described using a Debye–Hückel-type approach. Here, we report the novel observation of a significant departure from the Debye–Hückel-type model in high-energy-density matter by probing laser-driven, shock-compressed plastic with high-energy X-rays. We use spectrally resolved X-ray scattering in a geometry that enables direct investigation of the screening cloud, and demonstrate that the observed elastic scattering amplitude is only well described within a more general approach.
Similar content being viewed by others
Introduction
The Debye–Hückel theory of charge screening1 is one of the seminal results of electrolyte and plasma physics, whereupon the long-range Coulomb force between a pair of charge carriers is replaced by an exponentially decaying, short-range potential to account for the interaction with the surrounding medium. The extension of this classical description to degenerate electrons, Thomas–Fermi screening2,3, plays an important role in the description of electrons in solids4 and warm dense matter (WDM)5. Moreover, the concept of Debye-like screening is also applied to describe a large range of systems usually not associated with particles interacting via Coulomb forces, including quark-gluon matter6,7, ultra-cold systems in traps8 or chemical and biological systems9. For all these cases, the structural, thermodynamic, transport or relaxation properties are determined by effective interactions between localized charges.
Although Debye-like screening has been applied for many applications, it contains a number of inherent restrictions: weakly interacting particles in the screening cloud, negligible dynamic evolution and the long-wavelength limit. The latter requirement is only fulfilled for processes with a wavelength much larger than the screening length, leading to a wavenumber k effectively approaching zero. In contrast, we will refer to static screening models, which are not similarly restricted, as finite-wavelength screening.
To investigate static screening beyond the long-wavelength limit, processes involving large momentum transfers need to be considered. Large-angle collisions are a typical process where finite-wavelength screening could be observed. Although such collisions are highly unlikely in ideal, low-density plasmas, strong scattering is known to modify transport and relaxation properties in dense plasmas10,11. Another possibility to investigate deviations from Debye-like screening is the interaction of X-rays with dense matter under large scattering angles12. Indeed, spectrally resolved X-ray Thomson scattering (XRTS)13 is particularly suited for these investigations, as it simultaneously allows for the determination of the plasma conditions and the study of the screening cloud from a single spectrum.
In the following, we report observations of finite-wavelength screening in dense matter probed via spectrally resolved XRTS on laser-driven, shock-compressed plastic (CH) capsules. The strength of the elastic Rayleigh feature is used to further constrain simultaneous measurements of the electron density, temperature and mean ionization obtained from the inelastic Compton feature. We show that agreement between modelled and measured values for the Rayleigh amplitude can only be obtained if finite-wavelength screening is considered.
Results
X-ray scattering in high-energy-density matter
Investigation of the screening cloud at large wavenumbers requires highly compressed and moderately heated matter to be probed by high-energy X-rays. These target conditions were created at the Omega laser facility by compressing thin CH shells with multiple coalescent shocks driven by intense laser beams. He-α radiation at ∼9 keV from a laser-produced Zn plasma is used as the probe. Measuring the intensity of the X-ray scatter for angles between 120° and 150° accesses the response at the wavenumbers of interest (see Methods).
The spectrum of radiation scattered by a plasma is proportional to the double-differential scattering cross-section14
where σT=6.65 × 10−25 cm2 is the Thomson cross-section. The frequency and wavenumber shifts associated with the scattering are ω=ωi–ωs and k=ki–ks, respectively. is the total dynamic structure factor containing information on the spatio-temporal correlations of the fully coupled electron system. For a composite plasma probed with high-energy X-rays, the latter can be decomposed into distinct contributions from elastic and inelastic scattering15. Both terms depend on the temperature, density and mean charge state of the ions. Thus, in principle, the thermodynamic state of the target may be determined by fitting the measured spectrum with theoretical calculations16,17,18.
It is well known that robust estimates of the electron density and temperature of the samples probed with XRTS can be made by matching to the inelastic Compton scattering feature of the spectrum19,20. Moreover, the relative contributions of free and bound electrons enable the ionization state, and therefore material density, to be measured21. Uncertainties in the models describing this feature are effectively constrained by probing at large k (see Methods). In this work, we obtain conditions suggestive of a WDM state with free-electron densities of ∼1024 cm−3 and temperatures of . Furthermore, we find the mean charge state of the carbon ions to be , consistent with estimates of continuum lowering in WDM.
In contrast, the elastic Rayleigh feature is the result of scattering from electrons that follow the low-frequency density fluctuations of the ions and, thus, is an excellent measure for the static ionic structure and the associated screening properties22,23,24,25. The elastic scattering strength in WDM with multiple ion species is given by26
where is the concentration of ion species a, with na its density, fa(k) is the form factor of bound electrons and qa(k) is the screening cloud13. The latter quantities are Fourier transforms of the respective electron densities around ions of species a. Sab(k) denotes the partial ion–ion structure factors. The form factors for bound electrons are readily obtained from first-principles methods for the tightly bound K-shell electrons in carbon.
For most conditions, the largest theoretical uncertainty in WR is related to the ionic structure factors (see, for example, ref. 27). In the present experiment, the ambiguity with respect to the ionic structure is circumvented by the large scattering angle applied as shown in Fig. 1a. At the resulting large k values (shaded vertical band), all spatial correlations have decayed and we have unity for Saa and negligible values for Sab. This result is independent of the method applied to obtain the ion structure (see Methods). Thus, the only unknown remaining in the description of the elastic scattering strength is the form of the screening cloud qa(k), which can be determined by our measurements.
Theoretical description of the screening cloud
In general, the response of the electrons to the ions may be derived within a quantum statistical framework from the ratio of dynamic structure factors, thereby rigorously incorporating strong interactions. As it is associated with the ions, the frequency dependence of the screening cloud is negligible since the electrons react to the ion distribution almost instantaneously, that is, electronic screening can be treated in the static limit ω→0. For the conditions of interest, electron–ion correlations are expected to be weak by virtue of the high Fermi energy, justifying a linear response approach28
in which the screening is provided by the dielectric function . Strong correlations between electrons can be accounted for by modifying the potential with local field corrections . However, a weak coupling approximation, Gee(k)≈0, may be assumed for the conditions created. The non-interacting density response is then given by the random phase approximation (RPA)5
where we have defined , and β=1/kBT. Furthermore, is the degeneracy parameter, with being the thermal de Broglie wavelength, and ηe=βμe is the dimensionless chemical potential. Note that equation (4) is valid for arbitrary k, making it suitable to describe finite-wavelength screening.
While the idea of such a generalized approach to screening is not new, until recently the vast majority of XRTS data has been evaluated on the basis of the Debye–Hückel or Thomas–Fermi theories, wherein the dielectric function has the form . The inverse screening length κe is chosen to be consistent with the appropriate limit, but may also be derived for arbitrary degeneracies5. One may readily obtain this form of the dielectric function by taking the long-wavelength (k→0) limit of equation (4). Thus, the Debye-like (Debye–Hückel and Thomas–Fermi) theories are restricted to small momentum transfer events, such as long-range interactions between particles.
Taking the electron–ion interaction to be the Coulomb potential, one may assume and subsequently find the well-known expression for the screening cloud
Reviewing equation (5) reveals that the ratio k/κe determines whether the long-wavelength limit is applicable. For collective scattering29,30, the wavelength of density fluctuations is larger than the screening length: α=κe/k>1. Accordingly, the long-wavelength limit is valid. Conversely, we have α<1 in the non-collective scattering regime, and the screening function arising from combining equations (3) and (4) should be used. In this regime, the tail of the screening cloud decays significantly faster than the 1/k2-scaling predicted by the Debye-like theories equation (5) as can be seen in Fig. 1b. Under non-degenerate conditions, however, the Debye-like behaviour is recovered for intermediate k before changing to scale as 1/k4 in the limit k→∞ (Fig. 1c). Note that the finite-wavelength results tend to a constant for large k as the electron density close to the ions becomes independent of the temperature due to quantum degeneracy.
Furthermore, this strong 1/k4 decay directly leads to the well-known Friedel oscillations of the real-space electron density distribution around the ions31, which have been observed in cold solids32. This effect has thus far been hidden in XRTS experiments by probing higher temperatures, use of smaller scattering angles or weak ionization. In our experiments, the large scattering angle combined with strongly driven WDM states results in an intermediate scattering regime with α∼0.3, where the screening cloud exhibits changes due to finite-wavelength effects and also gives a significant contribution to the elastic Rayleigh feature.
Data analysis
To determine the plasma conditions from the measured XRTS spectra, we apply a statistical analysis with respect to the probability density corresponding to χ2 statistics (see, for example, ref. 33) over a wide range of density–temperature space. The best fit is defined by minimizing the mean square deviation between data and model, that is χ2, whereas the quality and robustness of the fit is estimated by the values of χ2 in the density–temperature region around the best fit.
For these calculations, the mean charge state of carbon is held constant at for all parameters considered. This estimate is supported by previous observations of continuum lowering21 based on the shape of the inelastic Compton feature, and also by the excellent agreement we find comparing ab initio simulations based on density functional theory (DFT)-molecular dynamics to the predictions of the classical fluid theory for the partial structure factors (see Fig. 1a). Moreover, we have performed self-consistent calculations of the ionization equilibrium5 that include continuum lowering34 and use temperature and density profiles given by radiation hydrodynamics simulations of the target. Figure 2 clearly shows that the mean charge state in the bulk of the target is dominated by contributions from helium-like carbon.
Before we analyse the data using the screening models, we note that it is possible to obtain the Rayleigh weight directly from the experimental data. This is achieved by decomposing the total scattered power into elastic and inelastic contributions, as per the generalized Chihara formula26, and rearranging to obtain
Here, is the ratio of the amplitudes for elastic (ω=0) and inelastic (ω=ωpk≈ħk2/2me) scattering in the measured spectra, Σ(ω) is shape of the X-ray source function and WC(k, ω) is the calculated total inelastic contribution to the scattered power.
Using equation (6) to analyse the data, the observed shape and amplitude of the Rayleigh peak is perfectly reproduced. Since it does not depend on any models for the elastic contributions, this approach yields the experimentally determined value for the strength of the Rayleigh peak . The resulting conditions of the best fit are then largely determined by the sensitivity of the free-electron response to the density and temperature19. Under partially ionized conditions, one may additionally infer the existence or absence of contributions from particular bound states from the shape of the Compton peak; in particular, our data show the absence of L-shell states for the carbon ions21.
As demonstrated by Fig. 3, the fitting procedure results in a well-defined state with a clear peak in the probability density. The latter has a highly elongated shape along the temperature axis indicating much greater sensitivity to the electron density, as expected for partially degenerate conditions. For a pump–probe delay of 3.4 ns, the best fit occurs, for example, at ne=1.31(±0.37) × 1024 cm−3 and kBT=10.5(+6.5) eV. Extracted conditions for delays of 3.5 and 3.6 ns yield decreasing densities and temperatures consistent with cooling and expanding plasmas in the release states after shock compression. This trend agrees with the predictions of radiation hydrodynamics simulations.
The errors on our measurements are quantified by the contour delineating the region in density–temperature space containing 68.3% of the total probability, that is, the 1-sigma confidence interval. Using equation (6) to obtain , we find typical error bars of ±20–30% for the density, whereas only an upper limit, typically around 60%, can be rigorously deduced for the temperature due to the temperature insensitivity of the theoretical models for colder, strongly degenerate states. The sensitivity of the inelastic scattering spectra due to such errors is shown in Fig. 3d.
Well-defined conditions are also found when the strength of the Rayleigh feature, WR, is calculated using the Debye-like model, as well as for the finite-wavelength approaches for the screening cloud (see Fig. 3b,c). We note that in both cases the shape of the contours indicate roughly twice the density sensitivity and substantially improved accuracy with respect to the temperature as compared with the approach determining . This stems from the additional constraints due to the temperature and density dependence of the screening cloud. However, the different screening models yield quite different conditions. For example, at the time t=3.4 ns, the best fit using finite-wavelength screening is found at ne=1.52(±0.1) × 1024 cm−3 and kBT=9.3(±2.4) eV, whereas for the Debye-like model, the best fit occurs for ne=4.8(±0.76) × 1023 cm−3 and kBT=21(±2) eV. Again, cooling and decompression is predicted for later times.
Applied to the same conditions, the elastic scattering strength predicted by different screening models significantly differs. In general, we find that the Rayleigh peak is too large if the k→0 (Debye-like) limit is applied, whereas the amplitude given by finite-wavelength screening agrees with the measurements within the error bars (Fig. 3e). Figure 4 compares the elastic Rayleigh strength obtained using the three approaches presented versus the ideal electron pressure (combining electron density and temperature into a single parameter). Here, we take the plasma conditions obtained from fitting with equation (6) as the reference states, and the corresponding errors are propagated through the calculation of WR for each model. We find significant deviations between the results using Debye-like screening and the extracted (experimental) values for all time delays. On the other hand, finite-wavelength screening is in much better agreement for 3.4 and 3.5 ns. For 3.6 ns, the plasma has cooled considerably and one expects strong correlations within the electron gas. The reduction of accuracy found for finite-wavelength screening points to an inadequate description within the RPA as applied here.
Up to now the free electrons have been assumed to respond to the ions via a screened Coulomb-like potential determined by the effective charge of the ions. Clearly, this is justified for the fully ionized hydrogen component. However, the electrons still bound to the carbon ions are known to modify the electron–ion interactions. Here, we estimate this effect using a soft-core pseudopotential similar to the well-known empty-core approach for metals near room temperature28. The cutoff radius of the ionic core (rc=29 pm) is taken from the effective size of C+4 ions. The resulting weaker electron–ion interactions lead to oscillations in Fourier space, and negative values of the screening cloud in the range of wavenumbers under study. This effect is exacerbated by the weak decay of the Debye–Hückel/Thomas–Fermi dielectric function in k-space. Indeed, within this pseudopotential approach, we do not find a distinct peak in the probability density.
The XRTS data and their analyses clearly show the inapplicability of the long-wavelength approximation inherent to the Debye–Hückel and Thomas–Fermi theories of screening for the high-density matter under investigation. In contrast, the full k-dependence of the screening function, which shows a stronger decay for larger wavenumbers and Friedel oscillations in real space near the ion, provides better agreement with measurements under weakly coupled conditions. This finding has important implications when modelling processes with large momentum transfers.
Presently, Debye-like screening builds the basis for theoretical investigations of, for example, the phase diagram35, the ion dynamics36, the nucleation times37 and relaxation properties10 of such states. Thus, our understanding of the particle and energy transport in stars or large planets as well as our modelling capabilities for inertial confinement fusion experiments must be improved with respect to better screening models. As the concept of Debye–Hückel-like screening is also used to describe the electron dynamics in solids38, quark-gluon matter6,7 and charged macromolecules8, the application of our finding goes far beyond the physics of high-energy-density matter.
Methods
Experimental details
The experiments were performed at the Omega Laser Facility at LLE in Rochester, USA, which provides multiple high-energy beams of 351-nm ultraviolet laser light39. Eight laser beams with a total energy of 13.5 kJ were used to launch multiple shocks into a 70-μm-thin CH shell. The conditions created by the coalescent shocks were then probed with different pump–probe delays21.
Another set of eight laser beams illuminated a Zn foil, creating a hot plasma that emits He-α radiation with a strong peak around 9 keV. These X-rays were used to probe the compressed CH and recorded using a time-gated microchannel plate coupled to a high-resolution highly-oriented pyrolytic graphite (HOPG) crystal spectrometer under scattering angles of 135°±15°. A gold cone protected the spectrometer from direct illumination from the Zn plasma. For non-relativistic probe energies, the momentum transfer from the photons to the plasma, for example, the k values probed, is given by the incident X-ray energy Ei and the scattering angle θ via k≈(2Ei/ħc)sin(θ/2) (ref. 13). In the present experiment, we have k=8.4–9.2 Å−1. In comparison, the inverse screening length in the plasma is κe∼2.5 Å−1 for the densities and temperatures achieved in the target, yielding α=κe/k∼0.3. Thus, for the relatively large k values probed, we expect non-collective scattering and significant contributions from the screening cloud.
Inelastic scattering contributions
Following the semi-classical approach of Chihara15, the dynamic structure factor is decomposed into three different contributions: the scattering from electrons following the ion motion (see equation (2)), scattering from free electrons in the continuum and scattering related to bound-free transitions driven by the incident X-rays13. The free-free contribution to the total inelastic scattering is given by the fluctuation-dissipation theorem, which has been evaluated within the RPA5.
Considerations beyond the RPA, such as nonlinear coupling between density fluctuations and collisions with the ions, have also been investigated. We estimate strong coupling in the electronic subsystem using static local field corrections5 and appeal to the Born–Mermin ansatz13 to include screened electron–ion collisions in first Born approximation. We find that nonideal effects beyond RPA yield only minor changes to our fitting procedure. Collisions do not noticeably modify the Compton feature owing to strong Pauli blocking of electron–ion-scattering channels under the degenerate conditions of interest. As shown in Fig. 1b, local field effects do not significantly influence the free-electron response at large k.
Finally, bound-free transitions have been treated within the impulse approximation13 due to the large k values considered and result in a low amplitude tail on the red wing of the inelastic feature. Fortunately, there is only a small overlap between the bound-free and free-free components of the inelastic feature for . Thus, the fitting procedure is not strongly sensitive to the bound-free contribution.
Elastic scattering contributions
Besides the screening function, the modelling of elastic scattering requires the form factor for electrons bound to the carbon ions fC(k) and the partial ion–ion structure factors Sab(k). The form factor can be obtained with high precision using DFT calculations, and is also well described (see Fig. 1b) using hydrogenic wavefunctions with tabulated effective screening constants13.
The partial static structure factors of the ion subsystem were calculated from both DFT coupled to molecular dynamics26 and via a multicomponent generalization of the hypernetted-chain equations from the classical fluid theory27. The structure predicted from DFT-molecular dynamics is well reproduced if the ion–ion potential is approximated with a Yukawa-like pseudopotential. We find that finite-wavelength screening of the effective ion–ion interaction is not important, as the inter-ionic separation is sufficiently large at the densities considered for the k→0 approximation to be reasonable. Further considerations such as the additional short-range repulsion resulting from the electrons bound to the carbon ions can be included using a simple extension25. The latter yields negligible improvements owing to the small size of the effective core radius in helium-like carbon. We find very small deviations from the ideal gas values (unity or zero) for the k values probed by our experiment, making our analysis very robust against uncertainties in the ionic structure. Thus, the only significant uncertainty lies in the description of the screening clouds.
Evaluation of modelling uncertainties
As implied by equation (6), the value of extracted from the experimental data is sensitive to the shape of the normalized X-ray source function Σ(ω). We have considered time-dependent profiles extracted from the data21, and also models such as Gaussian, Lorenztian and Voigtian functions. We find that the latter all require stronger elastic scattering to reproduce the amplitude of the observed peak, yielding typical uncertainties of ∼+10%. Moreover, we find that the asymmetric shape of the extracted source function, which implicitly contains the depth-broadening effect of the HOPG crystal40, is crucial to simultaneously fitting all spectral features of the data.
The strong angular divergence of the source can be accounted for by averaging calculations over a weighted distribution of scattering angles. Such a consideration has been shown to be important in the collective scattering regime29. Here, the scattering is non-collective for all scattering angles considered, and the net effect on the spectrum is a small degree of additional blurring in frequency space. A small impact on the Rayleigh weight stems from enhanced inelastic weight at ω=0 and the explicit k-dependence of the terms of equation (2). We find that changes to the functional form and width of the weighting distribution yield uncertainties that are well within an accuracy of approximately −5%.
Finally, recent work41 for highly compressed plastic targets points to attenuation of the probe X-rays through the target as another potential source of error in WR. In the present experiment, the much smaller size and lesser degree of compression of the target at the time of measurement ensures that negligible attenuation. Furthermore, radiation hydrodynamics simulations of the implosion suggest that the driven shell of material yields a fairly uniform density distribution in the region that dominates the scattering21.
In summary, estimates of error sources have been determined to be on the order of ɛmod∼+10%/−5%. These errors are combined with the statistical errors arising from the fitting procedure by summing in quadrature, .
Additional information
How to cite this article: Chapman, D.A. et al. Observation of finite-wavelength screening in high-energy-density matter. Nat. Commun. 6:6839 doi: 10.1038/ncomms7839 (2015).
References
Debye, P. & Hückel, E. Zur Theorie der Elektrolyte. I. Gefrierpunktserniedrigung und verwandte Erscheinungen. Phys. Zeitschr. 24, 185–206 (1923).
Thomas, L. H. The calculation of atomic fields. Proc. Camb. Philos. Soc. 23, 542–548 (1927).
Fermi, E. A statistical method for the determination of some atomic properties and the application of this method to the theory of the periodic system of elements. Z. Phys. 48, 73–79 (1928).
Ashcroft, N. W. & Mermin, N. D. Solid State Physics Holt, Rinehart & Winston (1976).
Kremp, D., Schlanges, M. & Kraeft, W. D. Quantum Statistics of Nonideal Plasmas Springer (2005).
Braaten, E. & Nieto, A. Next-to-leading order Debye mass for the quark-gluon plasma. Phys. Rev. Lett. 73, 2402–2404 (1994).
Heinz, U. Quark-gluon transport theory I. The classical theory. Ann. Phys. 161, 48–80 (1985).
Bannasch, G., Killian, T. C. & Pohl, T. Strongly coupled plasmas via Rydberg blockade of cold atoms. Phys. Rev. Lett. 110, 253003 (2013).
Yu., A., Grosberg, T. T., Nguyen & Shklovskii, B. I. The physics of charge inversion in chemical and biological systems. Rev. Mod. Phys. 74, 329–345 (2002).
Gericke, D. O., Murillo, M. S. & Schlanges, M. Dense plasma temperature equilibration in the binary collision approximation. Phys. Rev. E 65, 036418 (2002).
Gericke, D. O. & Schlanges, M. Beam-plasma coupling effects on the stopping power of dense plasmas. Phys. Rev. E 60, 904–910 (1999).
Glenzer, S. H. et al. Demonstration of spectrally resolved x-ray scattering in dense plasmas. Phys. Rev. Lett. 90, 175002 (2003).
Glenzer, S. H. & Redmer, R. X-ray Thomson scattering in high energy density plasmas. Rev. Mod. Phys. 81, 1625–1663 (2009).
Crowley, B. J. B. & Gregori, G. X-ray scattering in many-particle systems. New J. Phys. 15, 015014 (2013).
Chihara, J. Interaction of photons with plasmas and liquid metals—photoabsorption and scattering. J. Phys. Condens. Matter 12, 231–247 (2000).
Garcia-Saiz, E. et al. Probing warm dense lithium by inelastic x-ray scattering. Nat. Phys. 4, 940–944 (2008).
Lee, H. J. et al. X-ray Thomson scattering measurements of density and temperature in shock-compressed beryllium. Phys. Rev. Lett. 102, 115001 (2009).
Regan, S. P. et al. Inelastic x-ray scattering from shocked liquid deuterium. Phys. Rev. Lett. 109, 265003 (2012).
Kritcher, A. L. et al. In-flight measurements of capsule shell adiabats in laser-driven implosions. Phys. Rev. Lett. 107, 015002 (2011).
Fortmann, C. et al. Measurement of the adiabatic index in Be compressed by counterpropagating shocks. Phys. Rev. Lett. 108, 175006 (2012).
Fletcher, L. B. et al. Observations of continuum depression in warm dense matter with x-ray Thomson scattering. Phys. Rev. Lett. 112, 145004 (2014).
Riley, D. et al. Direct observation of strong coupling in a dense plasma. Phys. Rev. E 66, 046408 (2002).
Barbrel, B. et al. Measurement of short-range correlations in shock-compressed plastic by short-pulse x-ray scattering. Phys. Rev. Lett. 102, 165004 (2009).
Kritcher, A. L. et al. Measurements of ionic structure in shock-compressed lithium hydride from ultrafast x-ray Thomson scattering. Phys. Rev. Lett. 103, 245004 (2009).
Ma, T. et al. X-ray scattering measurements of strong ion-ion correlations in shock-compressed aluminium. Phys. Rev. Lett. 110, 065001 (2013).
Wünsch, K., Vorberger, J., Gregori, G. & Gericke, D. O. X-ray scattering as a probe for warm dense mixtures and high-pressure miscibility. EPL 94, 25001 (2011).
Wünsch, K., Hilse, K. P., Schlanges, M. & Gericke, D. O. Structure of strongly coupled multicomponent plasmas. Phys. Rev. E 77, 056404 (2008).
Gericke, D. O., Vorberger, J., Wünsch, K. & Gregori, G. Screening of ionic cores in partially ionized plasma within linear response. Phys. Rev. E 81, 065401(R) (2010).
Glenzer, S. H. et al. Observation of plasmons in warm dense matter. Phys. Rev. Lett. 98, 065002 (2007).
Neumeyer, P. et al. Plasmons in strongly coupled and shock-compressed matter. Phys. Rev. Lett. 105, 075003 (2007).
Friedel, J. Electronic structure of primary solid solutions in metals. Adv. Phys. 3, 446–507 (1954).
Ono, M. et al. Electrostatic potential screened by a two-dimensional electron system: a real-space observation by scanning-tunnelling spectroscopy. Phys. Rev. Lett. 96, 016801 (2006).
Bevington, P. R. & Robinson, D. K. Data Reduction and Error Analysis for the Physical Sciences McGraw Hill (2003).
Stewart, J. C. & Pyatt, K. D. Jr. Lowering of ionization potentials in plasmas. Astrophys. J. 144, 1203–1211 (1966).
Hamaguchi, S., Faroukia, R. T. & Dunin, D. H. E. Phase diagram of Yukawa systems near the one-component-plasma limit revisited. J. Chem. Phys. 105, 7641–7647 (1996).
Vorberger, J., Donko, Z., Tkachenko, I. M. & Gericke, D. O. Dynamic ion structure factor of warm dense matter. Phys. Rev. Lett. 109, 225001 (2012).
Daligault, J. Crystal nucleation in the one-component plasma. Phys. Rev. E 73, 056407 (2006).
Kaiser, A., Rethfeld, B., Vicanek, M. & Simon, G. Microscopic processes in dielectrics under irradiation by subpicosecond laser pulses. Phys. Rev. B 61, 11437 (2000).
Boehly, T. R. et al. Initial performance results of the OMEGA laser system. Opt. Commun. 133, 495–506 (1997).
Pak, A. et al. X-ray line measurements with high-efficiency Bragg crystals. Rev. Sci. Instrum. 75, 3747–3749 (2004).
Chapman, D. A. et al. Simulating x-ray Thomson scattering signals from high-density, millimetre scale plasmas at the National Ignition Facility. Phys. Plasmas 21, 082709 (2014).
Acknowledgements
D.A.C., R.A.B. and D.O.G. thank the UK’s EPSRC for support. Part of this work was performed by the assistance of Lawrence Livermore National Laboratory under Contract DE-AC52-07 NA27344 and supported by LDRD grant 11-ER-050, as well as the NLUF Grants DE-FG52-07 NA28057 and DE-FG52-09 NA29035. The work by SLAC HED was supported by DOE Office of Science, Fusion Energy Science under FWP 100182. R.W.F. acknowledges support from the NLUF Grant DE-NA0000910 and DOE NNSA award DE-FG52-10 NA29649.
Author information
Authors and Affiliations
Contributions
D.A.C., J.V., R.A.B. and D.O.G. developed and applied the theory presented. J.V. performed the DFT-MD simulations. L.B.F., L.D., T.D., R.W.F., S.H.G., G.G., A.L.K., T.M. and A.E.P. carried out the experiment and the data were analysed by D.A.C., L.B.F. and T.M.G. D.A.C. and D.O.G. wrote the paper. S.H.G., G.G., O.L.L. and D.O.G. provided additional support for the experiment, data analysis and interpretation.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Chapman, D., Vorberger, J., Fletcher, L. et al. Observation of finite-wavelength screening in high-energy-density matter. Nat Commun 6, 6839 (2015). https://doi.org/10.1038/ncomms7839
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms7839
This article is cited by
-
Observing the onset of pressure-driven K-shell delocalization
Nature (2023)
-
A strong diffusive ion mode in dense ionized matter predicted by Langevin dynamics
Nature Communications (2017)
-
Theory of Thomson scattering in inhomogeneous media
Scientific Reports (2016)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.