Electrical conductivity of hot expanded aluminum: experimental measurements and ab initio calculations

Atomistic-Continuum Study of an Ultrafast Melting Process Controlled by a Femtosecond Laser-Pulse Train

In femtosecond laser fabrication, the laser-pulse train shows great promise in improving processing efficiency, quality, and precision. This research investigates the influence of pulse number, pulse interval, and pulse energy ratio on the lateral and longitudinal ultrafast melting process using an experiment and the molecular dynamics coupling two-temperature model (MD-TTM model), which incorporates temperature-dependent thermophysical parameters. The comparison of experimental and simulation results under single and double pulses proves the reliability of the MD-TTM model and indicates that as the pulse number increases, the melting threshold at the edge region of the laser spot decreases, resulting in a larger diameter of the melting region in the 2D lateral melting results. Using the same model, the lateral melting results of five pulses are simulated. Moreover, the longitudinal melting results are also predicted, and an increasing pulse number leads to a greater early-stage melting depth in the melting process. In the case of double femtosecond laser pulses, the pulse interval and pulse energy ratio also affect the early-stage melting depth, with the best enhancement observed with a 2 ps interval and a 3:7 energy ratio. However, pulse number, pulse energy ratio, and pulse interval do not affect the final melting depth with the same total energies. The findings mean that the phenomena of melting region can be flexibly manipulated through the laser-pulse train, which is expected to be applied to improve the structural precision and boundary quality.

Average-atom Ziman resistivity calculations in expanded metallic plasmas: Effect of mean ionization definition

We present calculations of electrical resistivity for expanded boron, aluminum, titanium, and copper plasmas using the Ziman formulation in the framework of the average-atom model. Our results are compared to experimental data, as well as to other theoretical calculations, relying on the Ziman and Kubo-Greenwood formulations, and based on average-atom models or quantum-molecular-dynamics simulations. The impact of the definition of ionization, paying particular attention to the consistency between the definition and the perfect free electron gas assumption made in the formalism, is discussed. We propose a definition of the mean ionization generalizing to expanded plasmas the idea initially put forward for dense plasmas, consisting in dropping the contribution of quasibound states from the ionization due to continuum ones. It is shown that our recommendation for the calculation of the quasibound density of states provides the best agreement with measurements.

Consistent approach for electrical resistivity within Ziman's theory from solid state to hot dense plasma: Application to aluminum

The approach presented in this work allows a consistent calculation of electrical conductivity of dense matter from the solid state to the hot plasma using the same procedure, consisting in dropping elastic scattering contributions to solid's and liquid's structure factors in the framework of the Ziman theory. The solid's structure factor was computed using a multiphonon expansion. The elastic part is the zero-phonon term and corresponds to Bragg peaks, thermally damped by Debye-Waller attenuation factors. For the liquid, a similar elastic contribution to the structure factor results from a long-range order persisting during the characteristic electron-ion scattering time. All the quantities required for the calculation of the resistivities are obtained from our average-atom model, including the total hypernetted-chain structure factor used from the liquid state to the plasma. No interpolation between two limiting structure factors is required. We derive the correction to apply to the resistivity in order to account for the transient long-range order in the liquid and show that it improves considerably the agreement with quantum-molecular dynamics simulations and experimental aluminum's isochoric and isobaric conductivities. Our results suggest that the long-range order in liquid aluminum could be a slightly compressed fcc one. Two series of ultrafast experiments performed on aluminum were also considered, the first one by Milchberg et al. using short laser pulses and the second one by Sperling et al. involving x-ray heating and carried out on the Linac Coherent Light Source facility. Our attempts to explain the latter assuming an initial liquid state at an ion temperature much smaller than the electron one suggest that the actual initial state before main heating is neither perfectly solid nor a normal liquid.

Liquid-Liquid Phase Transitions in Silicon

We use computationally simple neutral pseudoatom ("average atom") density functional theory (DFT) and standard DFT to elucidate liquid-liquid phase transitions (LPTs) in liquid silicon. An ionization-driven transition and three LPTs including the known LPT near 2.5  g/cm^{3} are found. They are robust even to 1 eV. The pair distributions functions, pair potentials, electrical conductivities, and compressibilites are reported. The LPTs are elucidated within a Fermi liquid picture of electron scattering at the Fermi energy that complements the transient covalent bonding picture.

Isochoric, isobaric, and ultrafast conductivities of aluminum, lithium, and carbon in the warm dense matter regime

We study the conductivities σ of (i) the equilibrium isochoric state σ_{is}, (ii) the equilibrium isobaric state σ_{ib}, and also the (iii) nonequilibrium ultrafast matter state σ_{uf} with the ion temperature T_{i} less than the electron temperature T_{e}. Aluminum, lithium, and carbon are considered, being increasingly complex warm dense matter systems, with carbon having transient covalent bonds. First-principles calculations, i.e., neutral-pseudoatom (NPA) calculations and density-functional theory (DFT) with molecular-dynamics (MD) simulations, are compared where possible with experimental data to characterize σ_{ic}, σ_{ib}, and σ_{uf}. The NPA σ_{ib} is closest to the available experimental data when compared to results from DFT with MD simulations, where simulations of about 64-125 atoms are typically used. The published conductivities for Li are reviewed and the value at a temperature of 4.5 eV is examined using supporting x-ray Thomson-scattering calculations. A physical picture of the variations of σ with temperature and density applicable to these materials is given. The insensitivity of σ to T_{e} below 10 eV for carbon, compared to Al and Li, is clarified.

Brewster angle and reflectivity of optically nonuniform dense plasmas

We provide theoretical analysis of the reflectance of shock-compressed plasmas and warm dense matter for normal incidence of laser radiation as well as for the dependence of s- and p-polarized reflectivity on the incidence angle. The self-consistent approach for the calculation of the optical and electronic properties of warm dense matter and nonideal plasmas developed in our previous works is extended for the description of normal and polarized reflectivity from the broadened optically nonuniform medium. Two methods are applied for the calculation of the reflectivity from spatially broadened optically nonuniform medium. The first one is based on the solution of the Helmholtz equation for the amplitudes of the electromagnetic field. Another one is based on Drude theory of reflection. It allows us to calculate the ratio of the s- and p-polarized reflectivity if dependence of the dielectric function on distance is known. For the case of the polarized reflectivity, the particular attention is concentrated on the Brewster angle. The calculation results for the dielectric function, obtained within the framework of the density-functional theory with the longitudinal expression for the dielectric tensor, are applied for the calculation of the reflectivity. Comparison with the experimental data for shock-compressed xenon is performed.

dc conductivity of two-temperature warm dense gold

Using recently obtained ac conductivity data we have derived dc conductivity together with free electron density and electron momentum relaxation time in two-temperature warm dense gold with energy density up to 4.1 MJ/kg (0.8×10^{11}J/m^{3}). The derivation is based on a Drude interpretation of the dielectric function that takes into account contributions of intraband and interband transitions as well as atomic polarizability. The results provide valuable benchmarks for assessing the extended Ziman formula for electrical resistivity and an accompanying average atom model.

Ab initio calculation of shocked xenon reflectivity

Reflectivity of shocked compressed xenon plasma is calculated within the framework of the density functional theory approach. Dependencies on the frequency of incident radiation and on the plasma density are analyzed. The Fresnel formula for the reflectivity is used. The longitudinal expression in the long-wavelength limit is applied for the calculation of the imaginary part of the dielectric function. The real part of the dielectric function is calculated by means of the Kramers-Kronig transformation. The results are compared with experimental data. The approach for the calculation of plasma frequency is developed.

X-ray diagnosis of the pressure induced Mott nonmetal-metal transition

The evolution of the K-edge x-ray absorption near-edge spectroscopy (XANES) spectrum is investigated for an aluminum plasma expanding from the solid density down to 0.5  g/cm{3}, with temperatures lying from 5 down to 2 eV. The dense plasma is generated by nanosecond laser-induced shock compression. These conditions correspond to the density-temperature region where a metal-nonmetal transition occurs as the density decreases. This transition is directly observed in XANES spectra measurements through the progressive formation of a preedge structure for densities around 1.6  g/cm{3}. Ab initio calculations based on density functional theory and a jellium model have been efficiently tested through direct comparison with the experimental measurements and show that this preedge corresponds to the relocalization of the 3p atomic orbital as the system evolves from a dense plasma toward a partially ionized atomic fluid.

Experiments and simulations on hot expanded boron

We measured the thermodynamical and transport properties of boron in the warm dense matter regime ( 15000 K<T<25000 K and rho=0.094 g/cm3). Experimental data are compared with quantum molecular dynamics (QMD) simulations. We find a very good agreement between data and calculations, which permits us to transform experimental energies into temperatures, and allows us to compare conductivities with an average atom model coupled with a Kubo-Greenwood calculations. Contact is made between computationally intensive QMD simulation codes and fast average atom models.

Ab initio molecular dynamics simulations of dense boron plasmas up to the semiclassical Thomas-Fermi regime

We build an "all-electron" norm-conserving pseudopotential for boron which extends the use of ab initio molecular dynamics simulations up to 50 times the normal density rho0. This allows us to perform ab initio simulations of dense plasmas from the regime where quantum mechanical effects are important to the regime where semiclassical simulations based on the Thomas-Fermi approach are, by default, the only simulation method currently available. This study first allows one to establish, for the case of boron, the density regime from which the semiclassical Thomas-Fermi approach is legitimate and sufficient. It further brings forward various issues pertaining to the construction of pseudopotentials aimed at high-pressure studies.

Electronic structure and equation of state data of warm dense gold

Equation of state data and electrical resistivity of warm dense gold were measured in the internal energy range 8 - 12 MJ/kg. Experimental results were compared with quantum molecular dynamics simulations. The theoretical results match well the experimental data, allowing a detailed interpretation of the theoretical thermodynamic properties and frequency-dependent conductivities.

Plasma interferometry and how the bound-electron contribution can bend fringes in unexpected ways

Utilizing a new average atom code, we calculate the index of refraction in C, Al, Ti, and Pd plasmas and show many conditions over which the bound-electron contribution dominates the free electrons as we explore photon energies from the optical to 100 eV (12 nm) soft x rays. For decades measurement of the electron density in plasmas by interferometers has relied on the approximation that the index of refraction in a plasma is due solely to the free electrons and therefore is less than 1. Recent measurements of Al plasmas using x-ray laser interferometers observed fringes bending in the opposite direction than expected due to the bound-electron contribution causing the index of refraction to be larger than 1. During the next decade x-ray free-electron lasers and other sources will be available to probe a wider variety of plasmas at higher densities and shorter wavelengths, so understanding the index of refraction in plasmas is essential.

Ab-initio simulations of the optical properties of warm dense gold

Using a combination of classical and ab-initio molecular dynamics simulations, we calculate the structure and the electrical conductivity of warm dense gold during the first picoseconds after a short-pulse laser illumination. We find that the ions remain in their initial fcc structure for several picoseconds, despite electron temperatures ranging from a few to several eV after the laser illumination. The electrical conductivities calculated under these nonequilibrium conditions and using the latter assumption are in remarkable agreement with recent measurements using a short-pulse laser interacting with gold thin films.

Simulations of the optical properties of warm dense aluminum

Using quantum molecular dynamics simulations, we show that the optical properties of aluminum change drastically along the nonmetal metal transition observed experimentally. As the density increases and the many-body effects become important, the optical response gradually evolves from the one characteristic of an atomic fluid to the one of a simple metal. We show that quantum molecular dynamics combined with the Kubo-Greenwood formulation naturally embodies the two limits and provides a powerful tool to calculate and benchmark the optical properties of various systems as they evolve into the warm dense matter regime.