Researchers led by Bomi Kim of the Korea Institute of Geoscience and Mineral Resources, together with colleagues at Kongju National University, turned their attention to the Myeonsan titanium deposit in the Taebaeksan Basin of eastern South Korea. The deposit, hosted by the Cambrian Myeonsan Formation, is interpreted as a paleo-placer: titanium-bearing minerals such as rutile, ilmenite, and hematite were eroded from Precambrian basement rocks, transported by ancient rivers, and concentrated in conglomeratic sandstone during sedimentary deposition. Because the ore horizons are surrounded by a diverse suite of lithologies, including limestone, clastic sedimentary rocks, intrusive dykes, and granite, the site provides an ideal natural laboratory for testing which physical measurements can reliably separate ore-bearing rock from barren host rock.

The team selected 39 representative drill-core samples from borehole KDMS 2023-10, which penetrates the mineralized interval and its surroundings. The samples were divided into five groups: nine Myeonsan ore specimens, twelve limestones, four dyke rocks, eleven clastic sedimentary rocks, and three granites. Each specimen was prepared as a polished cylinder, a geometry that promotes uniform current distribution and simplifies the calculations needed for electrical measurements. Where core material was scarce, half-core semi-cylindrical specimens were prepared for the titanium ores, ensuring that every sample could pass through the full battery of tests.

That battery was deliberately broad. The researchers measured density and porosity using the Archimedes buoyancy method, following International Society for Rock Mechanics procedures: specimens were dried at 110 degrees Celsius for 24 hours, vacuum-saturated with water, and weighed in both saturated and submerged states. Magnetic susceptibility was recorded with a Bartington MS-3 meter at six positions on each sample to average out mineralogical anisotropy. Direct-current electrical resistivity, time-domain induced polarization (TDIP), and electrochemical impedance spectroscopy (EIS), known in geophysics as spectral induced polarization or SIP, were measured on fully saturated specimens using Solartron instruments with copper-mesh and filter-paper electrodes. Portable X-ray fluorescence (p-XRF) provided elemental compositions, and petrographic thin sections revealed the mineralogical texture behind the numbers.

The results deliver a clear hierarchy of diagnostic power. Density proved the most straightforward discriminator: the titanium ores averaged 3.17 grams per cubic centimeter, well above limestone at 2.74, dyke rocks at 2.73, clastic sediments at 2.84, and granite at 2.68. The elevated density reflects the abundance of heavy oxide minerals such as rutile and ilmenite packed into the ore horizons. The larger standard deviation among ore samples also revealed heterogeneous mineral enrichment within the mineralized interval, a useful warning that ore grades can vary significantly over short distances.

Chargeability, the capacity of a rock to store and slowly release electrical charge, produced an even more dramatic contrast. The ore samples averaged 91.57 millivolts per volt, roughly sixteen times higher than limestone at 5.64, dyke rocks at 3.41, clastic sediments at 4.62, and granite at 6.56. This signal arises because semiconductive ore minerals accumulate charge at mineral-fluid interfaces when an electric current is applied, then release it gradually after the current stops. In practical terms, an induced-polarization survey over this deposit should light up the ore zones while leaving most host rocks nearly invisible, making TDIP one of the most valuable exploration tools for sedimentary titanium deposits.

Not every measurement proved useful, and the negative results are as informative as the positive ones. Porosity values, ranging from 0.36 to 1.19 percent across groups, showed no meaningful separation between ore and host rocks, reflecting instead the compaction and cementation history of the sedimentary sequence. Magnetic susceptibility failed as a discriminator because the ore assemblage itself is magnetically quiet: rutile is only weakly magnetic and hematite is antiferromagnetic, so the intrusive dyke rocks, averaging 10.00 millisiemens per meter, outmagnetized the ores at 2.82. Electrical resistivity was similarly ambiguous, controlled mainly by pore structure and water saturation rather than ore content, with dyke and granite samples showing lower resistivities than the ores themselves.

The frequency-domain measurements added a powerful new dimension. In EIS spectra, several ore specimens, including MO1, MO2, MO5, and MO9, displayed strong phase responses at low frequencies below about 100 hertz, a signature of interfacial polarization between semiconductive ore minerals and pore fluids. These same samples showed the highest TDIP chargeability and elevated titanium contents above 3.11 weight percent and iron above 8.82 weight percent. Other ore samples and most host rocks instead showed phase values that increased toward higher frequencies, a pattern associated with dielectric polarization in insulating silicate and carbonate minerals. The two spectral fingerprints give geophysicists a way to distinguish genuine ore polarization from background dielectric effects, sharpening the interpretation of field survey data.

Petrography tied the physics back to mineralogy. Under the microscope, the ore specimens were dominated by opaque rutile, ilmenite, and hematite with minor quartz, while limestones were composed of calcite, dyke rocks of quartz, feldspar, and amphibole, clastic sediments of alternating coarse and fine silicate grains, and granites of quartz, feldspar, mica, and tourmaline. Sericitization of feldspar in the dykes and granites pointed to hydrothermal alteration. These observations directly explain the measured contrasts: dense, polarizable oxide minerals produce high density, high chargeability, and strong low-frequency phase responses, whereas insulating silicate and carbonate host minerals generate only weak, high-frequency dielectric behavior.

The study also exposed practical limits. One ore specimen, MO3, contained the highest titanium (9.85 weight percent) and iron (24.27 weight percent) of all samples yet showed weaker polarization than its peers, likely because p-XRF was measured at a single spot on a heterogeneous specimen. Because the core samples had to be preserved intact for the mining project, destructive analyses such as X-ray diffraction and ICP-MS were not possible. The authors conclude that spot geochemistry is best used to flag the presence of metallic minerals, while volume-integrating measurements such as density, TDIP, and EIS provide the more robust bulk characterization. Together, the integrated dataset establishes quantitative criteria for distinguishing titanium-bearing rocks from their hosts and demonstrates how combined petrophysical, geochemical, and mineralogical measurements can sharpen the search for a critical mineral whose secure supply increasingly depends on finding it faster and smarter.

Subject of Research: Integrated physical property measurements of titanium-bearing ores and host rocks at the Myeonsan deposit for mineral exploration

Article Title: Integrated physical property measurements of titanium-bearing and host rocks for mineral exploration

Article References: Kim, B., Shin, S., Shin, D., & Lee, G. (2026). Integrated physical property measurements of titanium-bearing and host rocks for mineral exploration. Environmental Earth Sciences, 85(16), Article 405. https://doi.org/10.1007/s12665-026-13140-8

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13140-8

Keywords: titanium, mineral exploration, Myeonsan deposit, induced polarization, spectral induced polarization, density, chargeability, paleo-placer deposit, petrophysics, critical minerals, South Korea, rutile

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