Matter and Radiation at Extremes, Volume. 7, Issue 3, 036901(2022)

A molecular dynamics study of laser-excited gold

Jacob M. Molinaa) and T. G. White
Author Affiliations
  • University of Nevada, Reno, Nevada 89557, USA
  • show less
    Figures & Tables(9)
    Comparison of the atomic mean square deviation and Debye temperature, TD, produced by interatomic potential developed by Sheng et al.34 and the values obtained from x-ray diffraction measurements at equilibrium temperatures up to melting by Syneček et al.54
    Comparison between experimentally obtained Laue decay curves10 and those given by molecular dynamics simulations. (a)–(c) The rms difference between simulated and experimental decay curves for a range of initial electron temperatures and electron–ion equilibration rates. Here, the rms differences of the (220) and (420) curves have been added in quadrature. The area enclosed by the dashed line encompasses the region in which the simulation results lie within the experimental error bars. (d)–(f) Comparisons, for each energy density, of the experimental and simulated decays of the (220) and (420) Laue peaks for best-fit values of Te0 and gei obtained in (a)–(c). (g)–(i) The corresponding temporal evolutions of the electron and ion temperatures. Here, the area shaded in blue represents the variation in the ion temperature across the sample.
    (a) Evolution of best-fit values of η for several functional forms of the electron–ion equilibration rate.21–23,25,26,29 Here, Medvedev et al.* is used to denote the gei(Te, Ti) prediction that is both electron- and ion-temperature-dependent. In our calculations of the average η for all theoretical models, we exclude the outlier results produced by the solely electron-temperature-dependent model calculated by Medvedev and Milov at a constant ion temperature of 300 K, and the Migdal et al. model. (b) Average value of gei for a given energy density up until melting. The value plotted has been calculated by taking the average over the gei value present out until 50, 25, and 10 ps in the low, intermediate, and high energy density cases, respectively.
    Comparison of simulated and measured spatially resolved diffraction pattern data for the low energy density case. (a)–(c) Spatially resolved data from the simulation for 100, 1000, and 3000 ps. (d) Angularly resolved diffraction data obtained from azimuthal integrals of the data shown in (a)–(c). (e)–(g) Spatially resolved data obtained by Mo et al.10 (h) Angularly resolved diffraction data obtained from azimuthal integrals of the data shown in (e)–(g).
    Same as Fig. 4, but for the intermediate energy density case and times of 20, 100, and 800 ps.
    Same as Fig.4, but for the high-energy density-case and times of −2, 7, and 17 ps.
    Visualization of gold melting in the low [(a)–(c)], intermediate [(d)–(f)], and high [(g)–(i)] energy density cases via the technique of polyhedral template matching70 in OVITO visualization software.69 For each energy density case, the time steps corresponding to the diffraction patterns shown in Figs. 4–6 are presented. The green color identifies the fcc lattice type, whereas the non-green colors identify liquid present in the sample. These results were produced using the smaller simulation geometry.
    Comparison between experimentally obtained Laue decay curves10 and those given by molecular dynamics simulations utilizing temperature-dependent electron–ion equilibration rates21–23,25,26,29 for the low (a), intermediate (b), and high (c) energy density cases. Here, “Medvedev et al.*” is used to denote a gei(T) prediction that is both electron- and ion-temperature-dependent.23 For all three energy density cases, the rms difference between simulated and experimental decay curves for a range of initial electron temperatures is shown for each of the employed gei(T) models. In our analysis, the rms differences of the (220) and (420) curves have been added in quadrature.
    • Table 1. Experimentally measured properties of face-centered cubic (fcc) gold alongside those reproduced in molecular dynamics simulations using the interatomic potentials developed by Sheng et al. and the electron-temperature-dependent (ETD) potential of Norman et al. Here, the values shown for the Norman et al. potential are for the variant of the potential calculated for Te ∼ 0.1 eV.

      View table
      View in Article

      Table 1. Experimentally measured properties of face-centered cubic (fcc) gold alongside those reproduced in molecular dynamics simulations using the interatomic potentials developed by Sheng et al. and the electron-temperature-dependent (ETD) potential of Norman et al. Here, the values shown for the Norman et al. potential are for the variant of the potential calculated for Te ∼ 0.1 eV.

      PropertySheng et al.34Norman et al.49Experiment
      Melt temperature (K)132012081337 (Ref. 50)
      Lattice constant (Å)4.0784.1554.078 (Ref. 51)
      Liquid density (g/cm3)17.1017.7017.10 (Ref. 52)
      Melting enthalpy (kJ/mol)11.1012.80 (Ref. 53)
    Tools

    Get Citation

    Copy Citation Text

    Jacob M. Molina, T. G. White. A molecular dynamics study of laser-excited gold[J]. Matter and Radiation at Extremes, 2022, 7(3): 036901

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Radiation and Hydrodynamics

    Received: Sep. 28, 2021

    Accepted: Mar. 3, 2022

    Published Online: Jan. 11, 2023

    The Author Email: Molina Jacob M. (jmmolina@nevada.unr.edu)

    DOI:10.1063/5.0073217

    Topics