The research addresses a puzzle that has long frustrated laser engineers. When a femtosecond pulse strikes a metal surface, the energy is deposited so quickly that ordinary notions of heating and melting break down entirely. The pulse duration is shorter than the time it takes for electrons to hand their energy over to the atomic lattice, meaning the material’s electrons are driven to enormous temperatures while the atoms themselves barely move at first. This state—known as two-temperature behavior—sits at the heart of why ultrashort pulses can ablate material with astonishing precision, leaving minimal heat-affected zones, recast layers, or burrs around the machined feature. Understanding exactly how the lattice then responds, and how material is ejected, requires a simulation approach that can resolve atomic motion in space and time simultaneously.

Molecular dynamics provides exactly that capability. In the study, the researchers modeled a titanium target as a vast ensemble of interacting atoms, governed by an interatomic potential calibrated to reproduce titanium’s structural, thermal, and mechanical properties. The laser pulse was incorporated through a two-temperature model, in which the absorbed optical energy first elevates the electron temperature, and energy then flows into the lattice via electron-phonon coupling. By coupling this continuum description of the electron subsystem to the atomistic dynamics of the lattice, the simulations captured the full sequence of events: from the instant of energy deposition through lattice heating, phase transformation, and the ultimate ejection of material from the irradiated zone.

One of the study’s central achievements is its systematic exploration of how the outcome depends on laser fluence—the energy delivered per unit area. At low fluences, just above the ablation threshold, the simulations reveal a delicate regime in which the near-surface region undergoes photomechanical stress confinement and fails through the generation and relaxation of intense compressive and tensile stress waves. The topmost atomic layers can be removed essentially intact, propelled outward by the release of stored thermoelastic stress, while the underlying crystal remains largely undisturbed. This gentle regime is precisely what practitioners prize for precision micromachining, because it minimizes collateral thermal damage and produces clean, well-defined surfaces.

As the fluence increases, the picture changes dramatically. The absorbed energy density climbs past the point where the lattice can remain a coherent solid, and the simulation shows the near-surface region superheating far beyond its equilibrium melting point. In this regime, the dominant material removal mechanism shifts toward explosive decomposition: the superheated, deeply undercooled liquid and even critical-point phenomena come into play, and the irradiated volume disintegrates into a mixture of vapor, clusters, and droplets. The researchers tracked the emergence of a foamy, low-density transient structure—sometimes called a subsurface bubble or cavitation zone—that expands from the center of the deposit and ultimately fragments, ejecting both atomic and nanocluster debris. These atomistic observations connect directly to experimental signatures such as the characteristic size distributions of nanoparticles collected in laser ablation plumes of titanium and other metals.

The simulations also shed light on the fate of the material left behind. Below the ablated layer, the models show rapid quenching at rates of trillions of kelvin per second, which can freeze in structural signatures quite unlike those of equilibrium titanium. Depending on depth and local energy density, the resolidified region can display amorphous character, disordered polycrystalline grains, or heavily twinned and defective crystal structures. Such subsurface defects influence surface roughness, hardness, residual stress, and even the biological response of titanium implants whose surfaces are laser-textured. By resolving these features at the atomic scale, the computational study provides a mechanistic bridge between processing parameters and the microstructure that ultimately determines device performance.

From an engineering standpoint, the value of this work lies in its ability to map the parameter space of femtosecond machining far more cheaply and comprehensively than experiment alone. Each simulation is, in effect, a virtual experiment in which fluence, pulse duration, number of pulses, and material temperature can be varied systematically, and every atom can be observed at every instant—something no microscope can achieve. The researchers analyzed how peak electron and lattice temperatures, stress profiles, and ablation depths evolve as a function of deposited energy, allowing them to identify thresholds separating stress-driven removal, phase-explosion-dominated ejection, and regimes where material is merely melted and resolidified without net removal. These thresholds correspond closely to the processing windows that laser manufacturers and job shops must navigate when optimizing titanium micromachining protocols.

The study also speaks to a long-standing debate in the ultrafast laser community about the relative importance of thermal and nonthermal mechanisms. In strongly absorbing metals excited below the threshold for nonlinear optical breakdown, the simulations support the conventional two-temperature picture: the energy deposition is thermal at the electron level, but the subsequent lattice response is so rapid and so far from equilibrium that classical thermal concepts such as boiling points lose their ordinary meaning. Instead, material removal is governed by the interplay of electron-phonon coupling strength, thermomechanical stress confinement, and the kinetics of melting and vaporization under extreme superheating. For titanium, whose electron-phonon coupling is comparatively strong, this coupling time is short enough that lattice heating begins within a few picoseconds, shaping the transition between the stress-dominated and thermally dominated ablation regimes.

The broader implications extend well beyond titanium. The methodological framework—combining a two-temperature description of laser energy deposition with large-scale molecular dynamics—is directly transferable to other transition metals, alloys, and multilayer thin films used in electronics, energy storage, and biomedical engineering. As femtosecond lasers move into high-throughput industrial settings, from drilling cooling holes in turbine blades to patterning stents and creating microtextured antibacterial surfaces, the demand for predictive process models is intensifying. Atomistic simulations of the kind reported here can feed mesoscale and continuum models, ultimately enabling digital twins of laser machining processes in which parameters are tuned in silico before a single physical part is machined.

There remain challenges on the path to fully predictive simulation. Molecular dynamics of this scale captures picoseconds of physical time, while real ablation plumes evolve over nanoseconds to microseconds, and multi-pulse processing introduces heat accumulation over far longer intervals. Experimental validation likewise demands ultrafast pump-probe diagnostics capable of watching plumes and surfaces evolve at the same temporal resolution. Nevertheless, this work marks a significant step forward in turning femtosecond laser machining from an empirically optimized craft into a quantitatively understood science. For a metal as strategically important as titanium—central to next-generation aircraft, prosthetic joints, and clean-energy hardware—knowing precisely how its atoms respond to the shortest light pulses humanity can generate is more than an academic curiosity. It is the foundation for manufacturing the components on which modern technology increasingly depends.

Subject of Research: Molecular dynamics simulation of femtosecond laser ablation of titanium

Article Title: Molecular dynamics study of femtosecond laser ablation of titanium

Article References: Parris, G., Goel, S., Nguyen, D. T., Salter, P., & Zhou, X. W. (2026). Molecular dynamics study of femtosecond laser ablation of titanium. npj Advanced Manufacturing. https://doi.org/10.1038/s44334-026-00114-8

Image Credits: AI Generated

DOI: 10.1038/s44334-026-00114-8

Keywords: femtosecond laser ablation, titanium, molecular dynamics, two-temperature model, electron-phonon coupling, phase explosion, ultrafast lasers, laser micromachining, advanced manufacturing, npj Advanced Manufacturing, Molecular, dynamics