SCIENTIFIC REPpRTS

OPEN Strongly-coupled plasmas formed from laser-heated solids

Received: 28 June 2015 Accepted: 01 October 2015 Published: 27 October 2015

M. Lyon1, S. D. Bergeson2, G. Hart2 & M. S. Murillo3

We present an analysis of ion temperatures in laser-produced plasmas formed from solids with different initial lattice structures. We show that the equilibrium ion temperature is limited by a mismatch between the initial crystallographic configuration and the close-packed configuration of a strongly-coupled plasma, similar to experiments in ultracold neutral plasmas. We propose experiments to demonstrate and exploit this crystallographic heating in order to produce a strongly coupled plasma with a coupling parameter of several hundred.

Strong coupling occurs when the inter-particle interaction energy in a system exceeds the random thermal energy. Strong coupling occurs in fields as diverse as quark-gluon plasmas1, ultracold atoms in the BEC-BCS crossover2,3, interactions at extremely high energy density4, quantum dots5, superconductivity6, and ultracold neutral plasmas7. When a system transitions into a regime in which the inter-particle interactions form the dominant energy scale, new physical effects can occur.

In plasma physics, strong coupling is predicted to have significant consequences. Strong coupling leads to long-range spatial ordering, causing the plasma to look less like an ionized gas and more like a Coulomb liquid. Long-range spatial order leads to cancelations in nearest-neighbor force terms. Consequently, ordering plays the role of shielding, leading to the paradoxical conclusion that ions of a given energy can approach one another more closely when the plasma is strongly coupled compared to the case when no order is present8-11. The effects of spatial ordering can enhance ion collision rates exponentially12-15.

In plasma physics, strong coupling is difficult to achieve. Strong (Coulomb) coupling is defined as the ratio of nearest-neighbor electrical potential energy to the kinetic energy,

Z 2e 2

4ne0aws kBT '

where Z is the ionization state, aws = (3/4nn)1/3 is the Wigner-Seitz radius, and n is the ion density. Notice that strong coupling, when r> 1, occurs for high Z, high density, and/or low temperature. Ionization typically occurs when kinetic energies exceed atomic ionization energies. The corresponding high-speed Coulomb collisions are relatively brief, infrequent, and long-range. The kinetic energy is high in most plasmas, resulting in low values of r ^ 1. It is possible that strong coupling is achieved in high-energy-density plasmas16. However, determining the plasma parameters to high precision can be challenging in these rapidly changing complex environments.

One barrier to achieving strong coupling in plasmas is called disorder-induced heating (DIH). This effect has been characterized in a number of ultracold neutral plasma (UNP) experiments3,17-20. In those experiments, laser-cooled atoms are resonantly ionized. As the plasma evolves, the rms ion velocity is measured directly. Although the plasmas are formed from essentially zero-temperature atoms, the initially non-interacting atoms suddenly experience the repulsive Coulomb interaction after ionization. The ions move to minimize their nearest-neighbor potential energy, and the electrical potential energy in the disordered system is converted to thermal energy. In these UNP experiments, the ion temperature is

Joint Quantum Institute and Department of Physics, University of Maryland, College Park, Maryland 20742, USA. 2Department of Physics and Astronomy, Brigham Young University, Provo, UT 84602, USA. 3New Mexico Consortium, Los Alamos, NM. Correspondence and requests for materials should be addressed to S.D.B. (email: scott.bergeson@byu.edu)

Figure 1. A representation of the alternating honeycomb lattice in graphite. Carbon atoms are located at intersections of thick black lines. Each carbon atom has three nearest neighbors located 0.142 nm away. The graphite unit cell is shown in the shaded gray box.

determined by the disorder in the initial system. The DIH process is rich in physics, as non-equilibrium dynamics can play a strong role in plasma evolution21.

This disorder-induced heating barrier could be overcome by ionizing a "pre-ordered" system. In laser-cooled gases, significant efforts are underway to do this using higher ionization states22, the Rydberg excitation blockade23-25, adiabatic expansion26, and other methods. However, avoiding DIH requires the spatial order in the initial system to match exactly the order that will exist in the strongly-coupled plasma state.

Disorder-induced heating can be a problem even in high density, high temperature plasmas18. In this paper we show that a form of DIH occurs in laser-produced plasmas when the crystallographic structure of the initial system does not overlap the close-packed configuration of a strongly coupled plasma27,28. Strictly speaking, this is not disorder-induced heating per se, but rather heating that arises when the pair distribution function in a given crystal type mismatches the close-packed pair distribution function in a strongly-coupled plasma. Even though the initial lattice may be highly ordered, it is "disordered" relative to the strongly-coupled plasma state. We call this heating effect crystal mismatch heating (CMH). We show that high energy-density plasmas can be strongly coupled if the initial system is appropriately chosen. We calculate the heating that results when different crystallographic configurations are impulsively ionized. For an appropriately chosen initial system, the value of the strong coupling parameter r should be high enough to observe Wigner crystallization in the plasma state.

An Example: a Strongly-Coupled Graphite Plasma

In ref. 29, Brown et al. describe an impressive experiment at the Linac Coherent Light Source in strongly-driven carbon. A thin graphite foil was illuminated with a short optical laser pulse. The laser pulse compressed and heated the sample to densities on the order of 1 to 2 x 1023 cm-3 and an estimated ion temperature of 5,000-10,000 K. The ion charge state was estimated to be Z = 4.5. Under these conditions, Eq. 1 gives r = 500. A highly-ordered Coulomb crystal was expected to form30. At a variable time after the optical pulse, a free electron laser x-ray pulse illuminated the melted carbon sample. The crystal structure of the sample was investigated using x-ray scattering31. The scattering data indicated that the plasma was in a partially correlated state but without any evidence of Wigner crystallization.

As we show in the next section, an analysis of this system based on crystallography suggests that the ion temperature is significantly higher than previously estimated. Carbon atoms in graphite are arranged in a honeycomb lattice, with alternating layers (see Fig. 1). This lattice is significantly different from the more energetically favorable configuration in which the ions are as far away from each other as possible, such as in a close-packed (CP) configuration of a face-centered cubic (FCC) or hexagonal close-packed (HCP) lattice27,28.

Estimating the ion temperature. The kinetic energy gained by the ions in the transition from the graphite lattice to the FCC lattice can be estimated using a simple model. Before ionization occurs, the nearest-neighbor potential energy in graphite is calculated using the shortest bond length. After the plasma relaxes, the nearest-neighbor potential energy is calculated using the characteristic FCC atom separation at the same average density, neglecting compression. The effects of compression on the value of the strong coupling parameter are discussed towards the end of this article.

Calculating the change in energy from the initial graphite structure to the final FCC lattice structure requires us to know the ion charge state and electron temperature because electron screening is important. Laser "strong-field" ionization has been studied at length (see, for example, ref. 32). Given the laser intensity reported by Brown et al., and considering the ionization energies of carbon33, the treatment in ref. 32 suggests that the ion charge state is probably closer to Z = 4. At solid densities, the probability of electron scattering during the laser pulse is high. Therefore the electron temperature is equal to the ponderomo-tive energy of a free electron in the electric field of the laser. In this case, Te = 60 eV for a laser intensity of 1015 W/cm2. Comparing the estimated ion charge state and plasma density using averaged-atom codes, the estimated range of temperatures is 50 to 100 eV34. The Debye length is AD = (s0kBTe/nee2)l/2 = 0.085 nm, where Te is the electron temperature and ne is the electron density.

We now estimate the potential energy for screened C4+ ions in the graphite lattice and in the FCC lattice at a density of 1.14 x 1023 cm-3 (solid density). In the graphite lattice, the nearest-neighbor distance

between atoms is a0 = 0.142 nm. One would expect the nearest-neighbor potential energy per atom to be

Z 2e 2

U 0 = "-exp (-K0) = 30eV>

4^0a0 (2)

where the factor of e-re° arises due to Yukawa screening, k0 = a0/AD, and Z = 4.

In the conventional FCC lattice, the unit cell contains 4 atoms. We can denote the lattice parameter, or the side-length of the cubic unit cell, as i. The nearest-neighbor distance between atoms, af, is the distance from the corner of the cell to the center of the cube face, a f = i/-/2 = 0.225nm at a density

of 1.14 x 1023 cm-3. The nearest-neighbor potential energy per atom in the FCC cell is

Z 2e 2

Uf = -exp (— kf) = 7 eV.

4^0af (3)

where Kf = af/AD. At this level of approximation, ignoring second-nearest-neighbor and other correlation effects, one would expect the ion temperature to be roughly the difference between these two energies, or approximately 23 eV. This is more than an order of magnitude higher than estimated by Brown et al. in ref. 29.

For this calculation, we have allowed the strongly-coupled plasma to relax to the FCC configuration27. The thermodynamics of strongly-coupled plasmas suggests that at low values of k, the plasma might relax into a BCC configuration. However, in a strongly-screened system, with k « 2, the Madelung energy is lower in the FCC configuration27,28.

CMH heating calculations. We calculate the ion temperature due to spatial structure mismatch more rigorously using molecular dynamics (MD) simulations35,36. In the simulations, the ions interact via the Yukawa potential, Uj(r) = (Zel4ne0ri) exp(-r ¡/AD), where rj is the distance between ions i and j. Although the number of particles per Debye sphere is small, first-principles calculations have shown that the Yukawa interaction can be appropriately used in this regime37. We initially place ions in the main simulation cell at locations corresponding to the initial crystal lattice. The equations of motion of the N particles are integrated using a second-order symplectic integrator (velocity-Verlet) subject to periodic boundary conditions. We find that N ~ 5000 ions is sufficient for the calculation to converge22.

In the calculation, careful attention to boundary conditions is required. Our computational grid is cubic and slightly mismatches the graphite lattice. In order to prevent this from dominating the temperature calculation, we adjust the ion positions near the boundaries to minimize edge-heating in the calculation. We also vary the size of the lattice in order to verify that this potential edge effect does not influence the result. Calculations of the ion temperature use ions located far from the simulation edges. "Movie"-like visualizations are used to verify that this systematic error is negligibly small.

The symmetry of the graphite lattice leads to exact force cancelation. Without a perturbation, no heating occurs. In our calculation, we use a variety of initial conditions to simulate a realistic laser-plasma experiment. For example, we start the simulation with the ions at zero temperature, but with their initial positions perturbed by a small amount. We also start the simulation with the ions in a perfect lattice, but with a Maxwellian velocity distribution. A range of calculations exploring these and other initial conditions show that the final ion temperature is independent of the initial perturbation.

The result of our calculation is shown in Fig. 2. The horizontal axis is scaled using the plasma period, _1/2

o>_1 = (ne2 / m i£0) = 8.3 fs. For the calculation shown, the ions are launched from a perfect lattice configuration. Their initial velocities form a Maxwellian distribution with an ion temperature of 0.1 eV - a temperature high enough that the ions start to move during the first few time steps in the calculation but small compared to the final ion temperature. The plotted quantity is the average kinetic energy per ion in the plasma. After roughly 10 plasma periods have passed, this average kinetic energy corresponds to the ion temperature. At earlier times, the ion temperature is not strictly defined because the system is out of equilibrium18. An initial period of slower heating is followed by a rapid temperature increase to approximately 14 eV. A modest overshoot is observed in the ion kinetic energy, as has been studied extensively in the ultracold neutral plasma community3,38. The final temperature in the plasma is within

Figure 2. MD calculation of the ion temperature in a graphite laser-produced plasma. This simulation was started with the ions in a perfect graphite crystal, but with a velocity distribution corresponding to an ion temperature of 0.1 eV. Regardless of how the simulation is started, the ion temperature always reaches the same final temperature. This illustrates that the final temperature is determined by the mismatch between the initial graphite lattice structure and the final FCC structure of the plasma.

a factor of two of the estimated 23 eV temperature due to changes in the potential energy described in the simple calculation given above.

Implications from Crystallography

Estimating the ion temperature due to crystal mismatch heating (CMH) can be more generally described. This factor-of-two calculation could be done using the atomic density and a quantity known as the atomic packing factor, A. The atomic packing factor is the fraction of the crystalline unit cell that is occupied by the atoms. The very compact FCC and HCP lattices have AF = 0.74. The slightly less compact body-centered cubic (BCC) lattice has AB = 0.68. The comparatively open diamond structure has AD = 0.34. Graphite, by comparison has AG = 0.17.

The distance between neighboring atoms in the initial crystal lattice is twice the atomic radius, 2ra. In terms of atomic packing factor, this is

3A v 4nn

When the atoms in this lattice are ionized and move into the FCC lattice, the new separation between atoms will be similarly defined. The distance between atoms will be 2rF, with rF = (3AF/4nn)1/3 and the density n equal to the density of the initial lattice.

With these definitions in mind, the difference in electric potential energy between the initial lattice and the final idealized FCC lattice can be written

kB Tion = °

Z 2e 2n1/3

I A i J

- e-KF

where the subscript i refers to the initial configuration and the subscript F refers to the FCC lattice. The numerical factor of 0.8 results from relating the density n to the distance between ions via the atomic packing factor A. When the differences in the screening length between the initial and final configurations are not too great, this simplifies to

kB T ion —

Z 2e 2n1/3 )-î

I A i J

This equation could be used to give factor-of-two estimates for the CMH ion temperature in laser-produced plasmas for any of the 14 Bravais lattices. This estimate is valid for pure elements only. It could be extended to mixed-species plasmas generated from alloys or molecular solids in a straight-forward manner.

r„ =

This treatment of CMH suggests that the major source of ion heating would be eliminated in laser-produced plasmas of atoms that are initially in an FCC or HCP lattice. This includes 90% of all metals and several other elements. This treatment also suggests that high ion temperatures will be generated in carbon graphite or diamond plasmas. Graphite pushers are used in laser-driven fusion work39. Mass effects aside, an aluminum pusher (FCC lattice) would produce a lower temperature ion plasma in the ablator. Substituting ablators of different crystallographic structures may make it possible to study the effects of ion temperature on the plasma evolution.

It may be possible to systematically explore the effects of CMH in laser-produced plasmas of the transition metals. The elements V, Cr, Mn, and Fe are all BCC at room temperature. Neighboring elements in the periodic table of Co, Ni, Cu, and Zn are all FCC or HCP. In this latter group, the effects of CMH will be absent. In the former group, CMH will be present. For an iron plasma with Z = 4 and Te = 60 eV, consistent with a laser-produced plasma using a laser intensity of 1015 W/cm2, Eq. (6) suggests a CMH temperature of 0.6 eV.

In strongly-coupled plasmas, the FCC and BCC structures are thermodynamically similar. Whether one structure or the other is preferred depends on k. Therefore, studying the influence of CMH using the FCC/BCC structure in metals may be challenging. Perhaps the best place to explore the effects of CMH are in carbon. Its different allotropes (graphite, diamond) give dramatically different values of the CMH. Its different geometries (bulk graphite, graphene, bucky balls, nanotubes) make it possible to explore geometric effects as well.

Additional Considerations for Strongly-Coupled Plasmas

In Yukawa systems, where electrons screen the ion-ion interaction, the meaning of the strong coupling parameter r defined in Eq. (1) is debatable because it is technically valid for a one-component-plasma. However, the potential energy is more complicated than this equation suggests, and tacking a screening factor like e-K onto Eq. (1) also fails to capture the essential physics40. If the determination of strong coupling is based on long-range ordering in the ion system, expressions have been derived that predict the electron and ion temperatures and densities needed to achieve a given value of the pair distribution function41.

The value of the strong-coupling parameter is influenced by plasma compression. Adiabatic compression can occur in solid-density laser-produced plasmas. Plasmas in this regime behave as an ideal gas42. Adiabatic compression of an ideal gas leads to the expression Tn-2/3 = constant. Given that Eq. (1) shows that rn1/3 = constant in a strongly coupled plasma, adiabatic compression predicts that Tan-1. As the density increases at constant Z, r decreases. For highly compressed plasmas, the increase of the average ionization state with density also influences r43.

Two additional heating mechanisms need to be considered in dense systems. Electron-ion thermali-zation heats the ions, reducing r further. For carbon at a density of 1023 cm-3 and an electron temperature of 60 eV, the electron-ion thermalization time is 2 ps44. In this short time, the electrons will heat the ions, reducing r even further. Three-body recombination and electron-Rydberg scattering heats the electrons and eventually also the ions45,46. The three-body recombination time in this system is less than 1 ps. Experiments in UNPs indicate that recombination occurs rapidly in strongly-coupled plasmas47, and this effect also needs to be considered in estimates of the time-evolving electron and ion temperatures in these systems.

All these considerations suggest that an experiment designed to explore the effects of strong coupling needs to meet several criteria. The plasma needs to be generated from a close-packed crystal. This will eliminate the influence of CMH. Heating due to electron-ion thermalization and three-body recombination can be minimized if the plasma expands before the electrons and ions reach a global equilibrium and before significant recombination has occurred48. The expansion velocity is vexp = (kBTe/mi)1/2 = 10,000 m/s for Fe ions in a plasma with an electron temperature of 60 eV. Given a three-body recombination time of ~1 ps, an initial plasma size should be small enough so that the plasma doubles in size in just 1 ps, or r0 = (1 ps) (10,000 m/s) = 10 nm. Alternatively, a larger plasma could be studied if fs-laser-based temperature diagnostics were used to evaluate the ion temperature in the sub-ps temperature regime.

Conclusion

We have shown that the initial crystal configuration of a solid can limit the maximum possible value of r in a laser-produced plasma. We have shown how to estimate and accurately calculate the ion temperature resulting from crystal-mismatch heating (CMH). At early enough times, before the ions have been heated by collisions with the electrons and before significant recombination has occurred but after several ion plasma periods, the ion temperature will be determined solely by CMH. In this regime, which is analogous to the many experiments in ultracold neutral plasmas, the value of the strong coupling parameter r should be high enough to observe Wigner crystallization in the plasma state.

References

1. Shuryak, E. Why does the quark-gluon plasma at RHIC behave as a nearly ideal fluidc? Prog. Part. Nucl. Phys. 53, 273-303

(2004).

2. Bourdel, T. et al. Experimental study of the BEC-BCS crossover region in lithium 6. Phys. Rev. Lett. 93, 050401 (2004).

3. Chen, Y. C. et al. Electron screening and kinetic-energy oscillations in a strongly coupled plasma. Phys. Rev. Lett. 93, 265003 (2004).

4. Nuckolls, J., Wood, L., Thiessen, A. & Zimmerman, G. Laser compression of matter to super-high densities: Thermonuclear (CTR) applications. Nature 239, 139-142 (1972).

5. Alivisatos, A. P. Semiconductor clusters, nanocrystals, and quantum dots. Science 271, 933-937 (1996).

6. McMillan, W. L. Transition temperature of strong-coupled superconductors. Phys. Rev. 167, 331-344 (1968).

7. Killian, T., Pattard, T., Pohl, T. & Rost, J. Ultracold neutral plasmas. Phys. Rep. 449, 77-130 (2007).

8. Salpeter, E. E. Electron screening and thermonuclear reactions. Austral. J. Phys. 7, 373-388 (1954).

9. Salpeter, E. E. & Van Horn, H. M. Nuclear Reaction Rates at High Densities. Ap. J. 155, 183 (1969).

10. Ichimaru, S. Plasm. Phys. (Benjamen/Cummings, Menlo Park, 1986).

11. Glinsky, M. E., O'Neil, T. M., Rosenbluth, M. N., Tsuruta, K. & Ichimaru, S. Collisional equipartition rate for a magnetized pure electron plasma. Phys. Fluids B 4, 1156-1166 (1992).

12. Ichimaru, S. Nuclear fusion in dense plasmas. Rev. Mod. Phys. 65, 255-299 (1993).

13. Jensen, M. J., Hasegawa, T., Bollinger, J. J. & Dubin, D. H. E. Rapid heating of a strongly coupled plasma near the solid-liquid phase transition. Phys. Rev. Lett. 94, 025001 (2005).

14. Dubin, D. H. E. Measurement of screening enhancement to nuclear reaction rates using a strongly magnetized and strongly correlated non-neutral plasma. Phys. Rev. Lett. 94, 025002 (2005).

15. Anderegg, F., Dubin, D., O'Neil, T. & Driscoll, C. Measurement of correlation-enhanced collision rates. Phys. Rev. Lett. 102, 185001 (2009).

16. Drake, R. P. High Energ. Dens. Phys. (Springer, Heidelberg, 2006).

17. Simien, C. E. et al. Using Absorption Imaging to Study Ion Dynamics in an Ultracold Neutral Plasma. Phys. Rev. Lett. 92, 143001 (2004).

18. Murillo, M. S. Ultrafast Dynamics of Strongly Coupled Plasmas. Phys. Rev. Lett. 96, 165001 (2006).

19. Morrison, J. P., Rennick, C. J., Keller, J. S. & Grant, E. R. Evolution from a molecular rydberg gas to an ultracold plasma in a seeded supersonic expansion of no. Phys. Rev. Lett. 101, 205005 (2008).

20. Lyon, M. & Bergeson, S. D. The influence of electron screening on disorder-induced heating. J. Phys. B 44, 184014 (2011).

21. Clerouin, J. et al. Evidence for out-of-equilibrium states in warm dense matter probed by X-ray Thomson scattering. Phys. Rev. E 91, 011101 (2015).

22. Lyon, M., Bergeson, S. D., Diaw, A. & Murillo, M. S. Using higher ionization states to increase coulomb coupling in an ultracold neutral plasma. Phys. Rev. E 91, 033101 (2015).

23. Bannasch, G. & Pohl, T. Rydberg-atom formation in strongly correlated ultracold plasmas. Phys. Rev. A 84, 052710 (2011).

24. Bannasch, G., Killian, T. C. & Pohl, T. Strongly coupled plasmas via rydberg blockade of cold atoms. Phys. Rev. Lett. 110, 253003

(2013).

25. Robert-de Saint-Vincent, M. et al. Spontaneous avalanche ionization of a strongly blockaded rydberg gas. Phys. Rev. Lett. 110, 045004 (2013).

26. McQuillen, P., Langin, T., Strickler, T. & Killian, T. Adiabatic expansion cooling of ions in ultracold neutral plasmas. Bull. Am. Phys. Soc. 59 (2014).

27. Farouki, R. T. & Hamaguchi, S. Thermodynamics of strongly-coupled yukawa systems near the one-component-plasma limit. ii. molecular dynamics simulations. J. Chem. Phys. 101, 9885-9893 (1994).

28. Ichimaru, S. Strongly coupled plasmas: high-density classical plasmas and degenerate electron liquids. Rev. Mod. Phys. 54, 1017-1059 (1982).

29. Brown, C. R. D. et al. Evidence for a glassy stste in strongly driven carbon. Sci. Rep. 4, 5214 (2014).

30. Dubin, D. H. E. First-order anharmonic correction to the free energy of a coulomb crystal in periodic boundary conditions. Phys. Rev. A 42, 4972-4982 (1990).

31. Vorberger, J. & Gericke, D. O. Ab initio approach to model x-ray diffraction in warm dense matter. Phys. Rev. E 91, 033112 (2015).

32. Augst, S., Strickland, D., Meyerhofer, D. D., Chin, S. L. & Eberly, J. H. Tunneling ionization of noble gases in a high-intensity laser field. Phys. Rev. Lett. 63, 2212-2215 (1989).

33. Kramida, A., Ralchenko, Y., Reader, J. & NIST ASD Team. NIST Atomic Spectra Database Ionization Energies Form (2015) Available at: http://physics.nist.gov/PhysRefData/ASD/ionEnergy.html (Accessed 22nd June 2015).

34. Murillo, M. S., Weisheit, J., Hansen, S. B. & Dharma-wardana, M. W. C. Partial ionization in dense plasmas: Comparisons among average-atom density functional models. Phys. Rev. E 87, 063113 (2013).

35. Murillo, M. S. Ultrafast dynamics of neutral, ultracold plasmas. Phys. Plasmas 14, 055702 (2007).

36. Whitley, H. D., Sanchez, D. M., Hamel, S., Correa, A. A. & Benedict, L. X. Molecular dynamics simulations of warm dense carbon. Contrib. Plasm. Phys. 55, 390-398 (2015).

37. Gilles, D., Lambert, F., Clerouin, J. & Salin, G. Yukawa monte carlo and orbital free molecular dynamics approaches for the equation of state and structural properties of hot dense matter. High Energ. Dens. Phys. 3, 95-98 (2007).

38. Chen, W.-T., Witte, C. & Roberts, J. Plasma oscillation damping in an ultracold neutral plasma. Bull. Am. Phys. Soc. 59 (2014).

39. MacKinnon, A. J. et al. High-density carbon ablator experiments on the national ignition facility. Phys. Plasmas 21, 056318

(2014).

40. Lyon, M., Bergeson, S. D. & Murillo, M. S. Limit of strong ion coupling due to electron shielding. Phys. Rev. E 87, 033101 (2013).

41. Ott, T., Bonitz, M., Stanton, L. G. & Murillo, M. S. Coupling strength in coulomb and yukawa one-component plasmas. Phys. Plasmas 21, 113704 (2014).

42. Laha, S. et al. Experimental Realization of an Exact Solution to the Vlasov Equations for an Expanding Plasma. Phys. Rev. Lett. 99, 155001 (2007).

43. Murillo, M. S., Weisheit, J., Hansen, S. B. & Dharma-wardana, M. W. C. Partial ionization in dense plasmas: Comparisons among average-atom density functional models. Phys. Rev. E 87, 063113 (2013).

44. Gericke, D. O., Murillo, M. S. & Schlanges, M. Dense plasma temperature equilibration in the binary collision approximation. Phys. Rev. E 65, 036418 (2002).

45. Robicheaux, F. & Hanson, J. D. Simulated expansion of an ultra-cold, neutral plasma. Phys. Plasmas 10, 2217-2229 (2003).

46. Roberts, J. L., Fertig, C. D., Lim, M. J. & Rolston, S. L. Electron Temperature of Ultracold Plasmas. Phys. Rev. Lett. 92, 253003 (2004).

47. Denning, A., Bergeson, S. D. & Robicheaux, F. Measurement and simulation of laser-induced fluorescence from nonequilibrium ultracold neutral plasmas. Phys. Rev. A 80, 033415 (2009).

48. McQuillen, P., Strickler, T., Langin, T. & Killian, T. C. Ion temperature evolution in an ultracold neutral plasma. Phys. Plasmas 22, 033513 (2015).

Acknowledgements

This research is supported in part by the Air Force Office of Scientific Research (Grant No. FA9950-12-0308).

Author Contributions

M.L. and M.M. conceived this project. G.H. provided crystallographic information. S.B. developed the packing fraction model. M.M. carried out the simulations. All authors reviewed and contributed to the manuscript.

Additional Information

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Lyon, M. et al. Strongly-coupled plasmas formed from laser-heated solids. Sci. Rep. 5, 15693; doi: 10.1038/srep15693 (2015).

l/jjv (J) I 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/