Kinetic limit of heterogeneous melting in metals

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.

Real-time (nanoseconds) determination of liquid phase growth during shock-induced melting

Melting of solids is a fundamental natural phenomenon whose pressure dependence has been of interest for nearly a century. However, the temporal evolution of the molten phase under pressure has eluded measurements because of experimental challenges. By using the shock front as a fiducial, we investigated the time-dependent growth of the molten phase in shock-compressed germanium. In situ x-ray diffraction measurements at different times (1 to 6 nanoseconds) behind the shock front quantified the real-time growth of the liquid phase at several peak stresses. These results show that the characteristic time for melting in shock-compressed germanium decreases from ~7.2 nanoseconds at 35 gigapascals to less than 1 nanosecond at 42 gigapascals. Our melting kinetics results suggest the need to consider heterogeneous nucleation as a mechanism for shock-induced melting and provide an approach to measuring melting kinetics in shock-compressed solids.

Ultrafast Energy Transfer Process in Confined Gold Nanospheres Revealed by Femtosecond X-ray Imaging and Diffraction

Femtosecond laser pulses drive nonequilibrium phase transitions via reaction paths hidden in thermal equilibrium. This stimulates interest to understand photoinduced ultrafast melting processes, which remains incomplete due to challenges in resolving accompanied kinetics at the relevant space-time resolution. Here, by newly establishing a multiplexing femtosecond X-ray probe, we have successfully revealed ultrafast energy transfer processes in confined Au nanospheres. Real-time images of electron density distributions with the corresponding lattice structures elucidate that the energy transfer begins with subpicosecond melting at the specimen boundary earlier than the lattice thermalization, and proceeds by forming voids. Two temperature molecular dynamics simulations uncovered the presence of both heterogeneous melting with the melting front propagation from surface and grain boundaries and homogeneous melting with random melting seeds and nanoscale voids. Supported by experimental and theoretical results, we provide a comprehensive atomic-scale picture that accounts for the ultrafast laser-induced melting and evaporation kinetics.

Direct Measurement of Ballistic and Diffusive Electron Transport in Gold

We experimentally show that the ballistic length of hot electrons in laser-heated gold films can exceed ∼150 nm, which is ∼50% greater than the previously reported value of 100 nm inferred from pump-probe experiments. We also find that the mean free path of electrons at the peak temperature following interband excitation can reach upward of ∼45 nm, which is higher than the average value of 30 nm predicted from our parameter-free density functional perturbation theory. Our first-principles calculations of electron-phonon coupling reveal that the increase in the mean free path due to interband excitation is a consequence of drastically reduced electron-phonon coupling from lattice stiffening, thus providing the microscopic understanding of our experimental findings.

Kinetics of laser-induced melting of thin gold film: How slow can it get?

Melting is a common and well-studied phenomenon that still reveals new facets when triggered by laser excitation and probed with ultrafast electron diffraction. Recent experimental evidence of anomalously slow nanosecond-scale melting of thin gold films irradiated by femtosecond laser pulses motivates computational efforts aimed at revealing the underlying mechanisms of melting. Atomistic simulations reveal that a combined effect of lattice superheating and relaxation of laser-induced stresses ensures the dominance of the homogeneous melting mechanism at all energies down to the melting threshold and keeps the time scale of melting within ~100 picoseconds. The much longer melting times and the prominent contribution of heterogeneous melting inferred from the experiments cannot be reconciled with the atomistic simulations by any reasonable variation of the electron-phonon coupling strength, thus suggesting the need for further coordinated experimental and theoretical efforts aimed at addressing the mechanisms and kinetics of laser-induced melting in the vicinity of melting threshold.

Melting at Mg/Al interface in Mg-Al-Mg nanolayer by molecular dynamics simulations

The melting at the magnesium/aluminum (Mg/Al) interface is an essential step during the fabrications of Mg-Al structural materials and biomaterials. We carried out molecular dynamics simulations on the melting at the Mg/Al interface in a Mg-Al-Mg nanolayer via analyzing the changes of average atomic potential energy, Lindemann index, heat capacity, atomic density distribution and radial distribution function with temperature. The melting temperatures (Tm) of the nanolayer and the slabs near the interface are significantly sensitive to the heating rate (vh) over the range ofvh ≤ 4.0 K ps^-1. The distance (d) range in which the interface affects the melting of the slabs is predicted to be (-98.2, 89.9) Å atvh→0,if the interface is put atd = 0 and Mg (Al) is located at the left (right) side of the interface. TheTmof the Mg (Al) slab just near the interface (e.g.d=4.0Å) is predicted to be 926.8 K (926.6 K) atvh→0,with 36.9 K (37.1 K) below 963.7 K for the nanolayer. These results highlight the importance of regional research on the melting at an interface in the nanolayers consisting of two different metals.

Kinetics of solid-liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal-liquid interface featuring an increase up to a maximum at 10-30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson-Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue 'Transport phenomena in complex systems (part 1)'.

Traveling waves of the solidification and melting of cubic crystal lattices

Using the phase field crystal model (PFC model), an analysis of slow and fast dynamics of solid-liquid interfaces in solidification and melting processes is presented. Dynamical regimes for cubic lattices invading metastable liquids (solidification) and liquids propagating into metastable crystals (melting) are described in terms of the evolving amplitudes of the density field. Dynamical equations are obtained for body-centered cubic (bcc) and face-centered cubic (fcc) crystal lattices in one- and two-mode approximations. A universal form of the amplitude equations is obtained for the three-dimensional dynamics for different crystal lattices and crystallographic directions. Dynamics of the amplitude's propagation for different lattices and PFC mode's approximations is qualitatively compared. The traveling-wave velocity is quantitatively compared with data of molecular dynamics simulation previously obtained by Mendelev et al. [Modell. Simul. Mater. Sci. Eng. 18, 074002 (2010)MSMEEU0965-039310.1088/0965-0393/18/7/074002] for solidification and melting of the aluminum fcc lattice.

Formation of Periodic Nanoridge Patterns by Ultrashort Single Pulse UV Laser Irradiation of Gold

A direct comparison of simulation and experimental results of UV laser-induced surface nanostructuring of gold is presented. Theoretical simulations and experiments are performed on an identical spatial scale. The experimental results have been obtained by using a laser wavelength of 248 nm and a pulse length of 1.6 ps. A mask projection setup is applied to generate a spatially periodic intensity profile on a gold surface with a sinusoidal shape and periods of 270 nm, 350 nm, and 500 nm. The formation of structures at the surface upon single pulse irradiation is analyzed by scanning and transmission electron microscopy (SEM and TEM). For the simulations, a hybrid atomistic-continuum model capable of capturing the essential mechanisms responsible for the nanostructuring process is used to model the interaction of the laser pulse with the gold target and the subsequent time evolution of the system. The formation of narrow ridges composed of two colliding side walls is found in the simulation as well as in the experiment and the structures generated as a result of the material processing are categorized depending on the range of applied fluencies and periodicities.

Heterogeneous melting kinetics in polycrystalline aluminum

The heterogeneous melting kinetics of polycrystalline aluminum is investigated by a theoretical model which represents the overall melting rate as a functional of the Weibull grain-size-distribution. It is found that the melting process is strongly affected by the mean-grain-diameter, but is insensitive to the shape parameter of the Weibull distribution. The temperature-time-transformation (TTT) diagrams are calculated to probe dependence of the characteristic timescale of melting on the overheating temperature and the mean-grain-diameter. The model predicts that the heterogeneous melting time of polycrystalline aluminum exponentially depends on temperature in high temperature range and the exponent constant is an intrinsic material constant independent of the mean-grain-diameter. Comparisons between TTT diagrams of heterogeneous melting and homogenous melting are also provided.

Visualization of ultrafast melting initiated from radiation-driven defects in solids

Materials exposed to extreme radiation environments such as fusion reactors or deep spaces accumulate substantial defect populations that alter their properties and subsequently the melting behavior. The quantitative characterization requires visualization with femtosecond temporal resolution on the atomic-scale length through measurements of the pair correlation function. Here, we demonstrate experimentally that electron diffraction at relativistic energies opens a new approach for studies of melting kinetics. Our measurements in radiation-damaged tungsten show that the tungsten target subjected to 10 displacements per atom of damage undergoes a melting transition below the melting temperature. Two-temperature molecular dynamics simulations reveal the crucial role of defect clusters, particularly nanovoids, in driving the ultrafast melting process observed on the time scale of less than 10 ps. These results provide new atomic-level insights into the ultrafast melting processes of materials in extreme environments.

Heterogeneous to homogeneous melting transition visualized with ultrafast electron diffraction

The ultrafast laser excitation of matters leads to nonequilibrium states with complex solid-liquid phase-transition dynamics. We used electron diffraction at mega-electron volt energies to visualize the ultrafast melting of gold on the atomic scale length. For energy densities approaching the irreversible melting regime, we first observed heterogeneous melting on time scales of 100 to 1000 picoseconds, transitioning to homogeneous melting that occurs catastrophically within 10 to 20 picoseconds at higher energy densities. We showed evidence for the heterogeneous coexistence of solid and liquid. We determined the ion and electron temperature evolution and found superheated conditions. Our results constrain the electron-ion coupling rate, determine the Debye temperature, and reveal the melting sensitivity to nucleation seeds.

Interatomic Potential in the Nonequilibrium Warm Dense Matter Regime

We present a new measurement of lattice disassembly times in femtosecond-laser-heated polycrystalline Au nanofoils. The results are compared with molecular dynamics simulations incorporating a highly optimized, embedded-atom-method interatomic potential. For absorbed energy densities of 0.9-4.3  MJ/kg, the agreement between the experiment and simulation reveals a single-crystal-like behavior of homogeneous melting and corroborates the applicability of the interatomic potential in the nonequilibrium warm dense matter regime. For energy densities below 0.9  MJ/kg, the measurement is consistent with nanocrystal behavior where melting is initiated at the grain boundaries.

Localized thin film damage sourced and monitored via pump-probe modulated thermoreflectance

Damage in the form of dewetting and delamination of thin films is a major concern in applications requiring micro- or nano-fabrication. In non-contact nanoscale characterization, optical interrogation must be kept to energies below damage thresholds in order to conduct measurements such as pump-probe spectroscopy. In this study, we show that the thermoreflectance of thin films can indicate the degree of film damage induced by a modulated optical heating source. By adjusting the absorbed power of the pump heating event, we identify the characteristics of the change in the thermoreflectance signal when leading up to and exceeding the damage threshold of gold films of varying thicknesses on glass substrates.

Vapor-deposited glasses of methyl-m-toluate: How uniform is stable glass transformation?

AC chip nanocalorimetry is used to characterize vapor-deposited glasses of methyl-m-toluate (MMT). Physical vapor deposition can prepare MMT glasses that have lower heat capacity and significantly higher kinetic stability compared to liquid-cooled glasses. When heated, highly stable MMT glasses transform into the supercooled liquid via propagating fronts. We present the first quantitative analysis of the temporal and spatial uniformities of these transformation fronts. The front velocity varies by less than 4% over the duration of the transformation. For films 280 nm thick, the transformation rates at different spatial positions in the film differ by about 25%; this quantity may be related to spatially heterogeneous dynamics in the stable glass. Our characterization of the kinetic stability of MMT stable glasses extends previous dielectric experiments and is in excellent agreement with these results.

Modified Z method to calculate melting curve by molecular dynamics

We extend the recently proposed Z method of estimating the melting temperature from a complete liquid and propose a modified Z method to calculate the melting temperature from a solid-liquid coexistence state. With the simulation box of rectangular parallelepiped, an initial structure of perfect lattice can run in the microcanonical ensemble to achieve steady solid-liquid coexistence state. The melting pressure and temperature are estimated from the coexistence state. For the small system with 1280 atoms, the simulation results show that the melting curve of copper has a good agreement with the experiments and is identical in accuracy with the results of the two-phase coexistence method with 24 000 atoms in the literature. Moreover, the method is conceptually simpler than the two-phase coexistence method.

Thermodynamic stability of nanometre-sized Cu₃Au systems

Classical molecular dynamics simulations have been used to investigate the order-to-disorder transition in bulk and nanometre-sized Cu₃Au systems. The Helmholtz free energy difference between ordered and disordered phases was evaluated at different temperatures through the Bennett's method. The change of free energy differences with temperature was employed to identify the transition temperatures. The obtained information was used to study the dynamics of a nanometre-sized Cu₃Au particle in a He gas thermostat at the transition temperature. It is shown that the system underwent rapid fluctuation in the chemical order degree related to the formation of vacancies and of atoms with defective coordination. Additional information on surface energies was also gained through a thermodynamic description of the transition process.

Solidification velocities in deeply undercooled silver

We measure the solidification velocity of pure Ag as a function of undercooling temperature from the melting point (Tm=1235 K) to 0.6Tm using ultrafast, pump-probe laser experiments. The thickness of the liquid layer, while it solidifies, is measured using optical third harmonic generation. We show that velocity reaches a maximum value at 0.85Tm, and then remains nearly constant with additional undercooling. These results contradict simple collision-limited models, but they are in good agreement with molecular dynamics simulations presented here, which show that the crystallization velocity depends weakly on temperature from 0.85Tm to less than approximately 0.1Tm.