Delineating the interplay effects of microstructure topology and residual stresses in ultrafast laser irradiated thin films
Pith reviewed 2026-05-22 05:07 UTC · model grok-4.3
The pith
Microstructure configuration ranks highest in controlling laser melting and expansion in gold thin films, ahead of topology, grain size, and orientation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Our results reveal a clear hierarchy of influence on laser-metal interaction: 1.) Microstructure configuration 2.) Topology 3.) Grain Size 4.) Crystallographic orientations. In fine-grained thin films, grain boundaries act as primary melting precursors, while local crystallographic orientation determines the melting extent in coarser grains. Residual tensile stresses contribute to higher melting and greater laser-induced expansion than unstrained films. Conversely, residual compressive stresses resist deformation, as deposited thermal energy is utilized to overcome lattice compression, leading to reduced expansion. We found that microstructure grain topology and size exert a stronger fingerp
What carries the argument
Hybrid Two Temperature Model-Molecular Dynamics simulations on microstructure-informed atomistic models that vary between randomized and equiaxed grain topologies, different grain sizes, and added tensile or compressive residual stress states.
If this is right
- Fine-grained films melt first at grain boundaries rather than inside grains.
- In coarser-grained films the amount of melting depends on the local crystal orientation of each grain.
- Tensile residual stress produces more melting and more overall expansion than films without stress.
- Compressive residual stress reduces expansion because laser energy first works against the lattice compression.
- Changes in grain topology and size affect final film expansion more than the initial density of defects.
Where Pith is reading between the lines
- Fabrication routes that deliberately create compressive stress could be used to limit unwanted expansion during laser processing of thin films.
- The same hierarchy of influences may hold for other metals or laser wavelengths, offering a general rule for choosing film properties in precision machining.
- Running targeted experiments on films with controlled grain topologies would test whether the simulated ranking survives real-world conditions.
- These rankings could help reduce trial-and-error in designing laser steps for microelectronics by predicting which film features to adjust first.
Load-bearing premise
The hybrid simulation method and the randomized versus equiaxed grain structures accurately represent the behavior of real fabricated gold thin films without needing experimental calibration or validation in the stress and size ranges examined.
What would settle it
Direct measurement showing that laser-irradiated gold films with compressive residual stress expand at least as much as tensile-stressed films, or that grain size controls melting more than topology does.
Figures
read the original abstract
Advanced nanodevices require high-precision machining of thin films using ultrafast lasers. However, thin-film fabrications cause variations in microstructure, crystallographic orientation, and residual stresses owing to coating conditions and substrate choice. This work investigates the complex interplay between these factors in ultrafast laser-irradiated gold (Au) thin films using a hybrid Two Temperature Model-Molecular Dynamics simulations. We realized microstructure-informed atomistic models with varying grain topologies (randomized vs. equiaxed), grain sizes, and residual tensile/compressive stress configurations. Our results reveal a clear hierarchy of influence on laser-metal interaction: 1.) Microstructure configuration 2.) Topology 3.) Grain Size 4.) Crystallographic orientations. In fine-grained thin films, grain boundaries act as primary melting precursors, while local crystallographic orientation determines the melting extent in coarser grains. Residual tensile stresses contribute to higher melting and greater laser-induced expansion than unstrained films. Conversely, residual compressive stresses resist deformation, as deposited thermal energy is utilized to overcome lattice compression, leading to reduced expansion. We found that microstructure grain topology and size exert a stronger fingerprint on film expansion than the initial defect density.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses hybrid Two-Temperature Model–Molecular Dynamics (TTM-MD) simulations to examine ultrafast laser irradiation of gold thin films with controlled variations in grain topology (randomized vs. equiaxed), grain size, crystallographic orientation, and imposed residual tensile or compressive stresses. It reports a hierarchy of influence on laser–metal interaction (microstructure configuration > topology > grain size > crystallographic orientation), with grain boundaries as primary melting precursors in fine-grained films, local orientation controlling melting extent in coarser grains, and tensile stresses increasing melting and expansion while compressive stresses reduce expansion relative to unstrained films. Microstructure topology and size are stated to dominate film expansion over initial defect density.
Significance. If the reported hierarchy proves robust under experimental calibration, the work would supply useful atomistic guidance for optimizing laser machining precision in thin-film nanodevices. The microstructure-informed modeling approach and direct comparison of randomized versus equiaxed topologies constitute a clear methodological strength. The quantitative distinctions drawn between stress states and between fine- versus coarse-grained regimes would be of practical value once anchored to measured melting fluences and expansion data for comparable Au films.
major comments (2)
- [Abstract and Results] Abstract and Results section: The central hierarchy (microstructure configuration > topology > grain size > crystallographic orientations) and the statements that grain boundaries are primary melting precursors in fine-grained films while orientation determines melting extent in coarser grains rest entirely on TTM-MD outputs. No comparison is presented between simulated melting fluence, lattice expansion, or defect evolution and experimental values for polycrystalline Au films of similar grain size or residual-stress state. Without this mapping, the ordering remains an untested prediction of the chosen interatomic potential, electron–phonon coupling parameters, and grain-boundary model.
- [Methods] Methods section: The hybrid TTM-MD implementation imposes residual stresses and constructs randomized versus equiaxed grain topologies, yet no sensitivity analysis or calibration against measured laser-ablation thresholds or thermal-expansion coefficients for the specific grain-size and stress regimes is reported. This directly affects the load-bearing claim that topology and size exert a stronger fingerprint on expansion than initial defect density.
minor comments (2)
- [Abstract] Abstract: the phrase 'stronger fingerprint on film expansion than the initial defect density' would benefit from an explicit definition of how defect density was quantified and from a quantitative metric (e.g., percentage difference in expansion) used to establish the comparison.
- [Figures] Figure captions and text: ensure consistent notation for the two grain topologies (randomized vs. equiaxed) and for the stress states (tensile, compressive, unstrained) across all panels and tables.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback and for recognizing the methodological strengths of our microstructure-informed TTM-MD simulations. We address each major comment below, maintaining a focus on the computational nature of the study.
read point-by-point responses
-
Referee: [Abstract and Results] Abstract and Results section: The central hierarchy (microstructure configuration > topology > grain size > crystallographic orientations) and the statements that grain boundaries are primary melting precursors in fine-grained films while orientation determines melting extent in coarser grains rest entirely on TTM-MD outputs. No comparison is presented between simulated melting fluence, lattice expansion, or defect evolution and experimental values for polycrystalline Au films of similar grain size or residual-stress state. Without this mapping, the ordering remains an untested prediction of the chosen interatomic potential, electron–phonon coupling parameters, and grain-boundary model.
Authors: We agree that the reported hierarchy is derived exclusively from controlled TTM-MD simulations and that no direct quantitative mapping to experimental melting fluences or expansion data for comparable polycrystalline Au films is provided. The interatomic potential (EAM), electron-phonon coupling, and grain-boundary construction follow standard literature values for gold, and the hierarchy is obtained by isolating each microstructural variable while holding all other simulation parameters fixed. This approach reveals relative influences that are difficult to disentangle experimentally. We will add a brief limitations paragraph in the Discussion section explicitly stating that the ordering constitutes a model prediction pending experimental calibration. revision: partial
-
Referee: [Methods] Methods section: The hybrid TTM-MD implementation imposes residual stresses and constructs randomized versus equiaxed grain topologies, yet no sensitivity analysis or calibration against measured laser-ablation thresholds or thermal-expansion coefficients for the specific grain-size and stress regimes is reported. This directly affects the load-bearing claim that topology and size exert a stronger fingerprint on expansion than initial defect density.
Authors: The referee is correct that no parameter sensitivity study or calibration to measured ablation thresholds appears in the Methods. Residual stresses were imposed via uniform lattice scaling, and the two topologies were generated with standard Voronoi-based procedures; defect densities were matched across compared configurations to isolate topology and size effects. We will expand the Methods section with additional justification for the chosen potential and coupling parameters drawn from prior Au studies, and we will qualify the expansion claim by noting that it holds within the simulated ensemble rather than as a calibrated prediction. revision: partial
Circularity Check
No circularity: hierarchy and stress effects are direct simulation outputs
full rationale
The paper reports results from hybrid TTM-MD simulations on atomistic models with imposed variations in grain topology (randomized vs. equiaxed), size, crystallographic orientations, and residual tensile/compressive stresses. The claimed hierarchy (microstructure configuration > topology > grain size > orientations), grain-boundary melting precursors in fine grains, orientation effects in coarse grains, and differential expansion under tensile vs. compressive stress are presented as emergent simulation observables rather than quantities derived from equations that reduce to the inputs by construction. No self-definitional loops, fitted parameters renamed as predictions, load-bearing self-citations, or smuggled ansatzes appear in the described methodology or abstract. The work is self-contained as a parameter-sweep computational study whose outputs are independent of any internal redefinition of the target quantities.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
hybrid Two Temperature Model-Molecular Dynamics simulations... hierarchy of influence... grain boundaries act as primary melting precursors... residual tensile stresses contribute to higher melting
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We realized microstructure-informed atomistic models with varying grain topologies (randomized vs. equiaxed), grain sizes, and residual tensile/compressive stress configurations
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Progressinultrafastlaserprocessingandfutureprospects
SugiokaK. Progressinultrafastlaserprocessingandfutureprospects. Nanophotonics, 6(2):393–413, 2017.doi:10.1515/nanoph-2016-0004
-
[2]
Shugaev MV, Wu C, Armbruster O, Naghilou A, Brouwer N, Ivanov DS, Derrien TJ-Y, Bulgakova NM, Kautek W, Rethfeld B, et al. Fundamentals of ultrafast laser–material interaction.MRS Bulletin, 41(12):960–968, 2016.doi:10.1557/mrs.2016.274
-
[4]
GuoB,SunJ,LuY,andJiangL. Ultrafastdynamicsobservationdur- ing femtosecond laser-material interaction.International Journal of Extreme Manufacturing, 1(3):032004, 2019.doi:10.1088/2631-7990/ ab3a24
-
[5]
Wang D, Wei S, Yuan X, Liu Z, Weng Y, Zhou Y, Xiao T-H, Goda K, Liu S, and Lei C. Ultrafast imaging for uncovering laser–material interaction dynamics.International Journal of Mechanical System Dynamics, 2(1):65–81, 2022.doi:10.1002/msd2.12024
-
[6]
Martin AA, Calta NP, Hammons JA, Khairallah SA, Nielsen MH, ShuttlesworthRM,SinclairN,MatthewsMJ,JeffriesJR,WilleyTM, et al. Ultrafast dynamics of laser-metal interactions in additive man- ufacturing alloys captured by in situ X-ray imaging.Materials Today Advances, 1:100002, 2019.doi:10.1016/j.mtadv.2019.01.001
-
[7]
Andreev SV, Aseev SA, Bagratashvili VN, Vorob’ev NS, Ishchenko AA, Kompanets VO, Malinovsky AL, Mironov BN, Timofeev AA, Chekalin SV, et al. Ultrafast transmission electron microscope for studyingthedynamicsoftheprocessesinducedbyfemtosecondlaser beams.QuantumElectronics,47(2):116,2017.doi:10.1070/QEL16276
-
[8]
Miloshevsky G. Ultrafast laser matter interactions: modeling ap- proaches,challenges,andprospects.ModellingandSimulationinMa- terials Science and Engineering, 30(8):083001, 2022.doi:10.1088/ 1361-651X/ac8abc
work page 2022
-
[9]
Zhigilei LV, Lin Z, and Ivanov DS. Atomistic modeling of short pulse laser ablation of metals: connections between melting, spal- lation, and phase explosion.The Journal of Physical Chemistry C, 113(27):11892–11906, 2009.doi:10.1021/jp902294m
-
[10]
Liu H, Xie W, Ding Y, Chen K, Wang S, Huo H, and Yang L. Re- view of molecular dynamics simulations in laser-based micro/nano- fabrication.Nanoscale, 16(46):21189–21215, 2024.doi:10.1039/ D4NR03305A
work page 2024
-
[11]
Song S, Lu Q, Zhang P, Yan H, Shi H, Yu Z, Sun T, Luo Z, and Tian Y. A critical review on the simulation of ultra-short pulse laser- metal interactions based on a two-temperature model (TTM).Optics & Laser Technology, 159:109001, 2023.doi:10.1016/j.optlastec. 2022.109001
-
[12]
YinCP,ZhangST,DongYM,YeQW,andLiQ. Molecular-dynamics study of multi-pulsed ultrafast laser interaction with copper.Ad- vances in Production Engineering & Management, 16(4):457–472, 2021.doi:10.14743/apem2021.4.413
-
[13]
The atomistic perspective of nanoscale laser ablation
Ivanov DS, Terekhin PN, Kudryashov SI, Klimentov SM, Kabashin AV, Garcia ME, Rethfeld B, and Zavestovskaya IN. The atomistic perspective of nanoscale laser ablation. InUltrafast Laser Nanos- tructuring: The Pursuit of Extreme Scales, pages 65–137. Springer, 2023.doi:10.1007/978-3-031-14752-4_2
-
[14]
ValavanisAS,ChenC,GrigoropoulosCP,EliceiriM,LiJ,andZhig- ilei LV. Laser fluence and film thickness dependence of the mecha- nisms of femtosecond laser ablation of ag films from atomistic simu- lationsandopticalimaging. InHigh-PowerLaserAblationVIII,page PC129391E. SPIE, 2024.doi:10.1117/12.3012635
-
[15]
XieL,LiY,WangF,YaoC,andYuanX. Moleculardynamicssimu- lationstudyontheablativecopperformingprocessattheatomiclevel by femtosecond laser.Journal of Manufacturing Processes, 150:38– 47, 2025.doi:10.1016/j.jmapro.2025.06.043
-
[16]
Amoruso S, Nedyalkov NN, Wang X, Ausanio G, Bruzzese R, and AtanasovPA. Ultrashort-pulselaserablationofgoldthinfilmtargets: Theory and experiment.Thin Solid Films, 550:190–198, 2014.doi: 10.1016/j.tsf.2013.10.165
-
[17]
Lian Y, Jiang L, Sun J, Lin G, and Liang M. Atomistic insight on temperature-dependent laser induced ultrafast thermomechanical re- sponse in aluminum film.International Journal of Heat and Mass Transfer, 231:125809, 2024.doi:10.1016/j.ijheatmasstransfer. 2024.125809
-
[18]
Rouleau CM, Shih C-Y, Wu C, Zhigilei LV, Puretzky AA, and Geo- heganDB. Nanoparticlegenerationandtransportresultingfromfem- tosecond laser ablation of ultrathin metal films: Time-resolved mea- surementsandmoleculardynamicssimulations.AppliedPhysicsLet- ters, 104(19), 2014.doi:10.1063/1.4876601
-
[19]
Zhang Z, Yang Z, Wang C, Zhang Q, Zheng S, and Xu W. Mech- anisms of femtosecond laser ablation of Ni3Al: Molecular dynam- ics study.Optics & Laser Technology, 133:106505, 2021.doi: 10.1016/j.optlastec.2020.106505
-
[20]
Hayder MM, Moumita TM, Chowdhury S, and Rahman KA. Effec- tive electronic properties and coupling for two-temperature model- molecular dynamics simulation of ultrafast laser ablation of nickel. Molecular Simulation, 50(14):1140–1151, 2024.doi:10.1080/ 08927022.2024.2385499
-
[21]
GanesanHandSandfeldS.Capturingthin-filmmicrostructurecontri- butions during ultrafast laser-metal interactions using atomistic sim- ulations.Materials & Design, page 114224, 2025.doi:10.1016/j. matdes.2025.114224
work page doi:10.1016/j 2025
-
[22]
Olbrich M, Punzel E, Lickschat P, Weißmantel S, and Horn A. In- vestigation on the ablation of thin metal films with femtosecond to picosecond-pulsed laser radiation.Physics Procedia, 83:93–103, 2016.doi:10.1016/j.phpro.2016.08.017
-
[23]
Ganesan H, Teijeiro C, and Sutmann G. Parallelization comparison and optimization of a scale-bridging framework to model cottrell at- mospheres.Computational materials science, 155:439–449, 2018. doi:10.1016/j.commatsci.2018.08.055
-
[24]
Ganesan H, Scheider I, and Cyron CJ. Quantifying the high- temperature separation behavior of lamellar interfaces in𝛾-titanium aluminide under tensile loading by molecular dynamics.Frontiers in materials, 7:602567, 2021.doi:10.3389/fmats.2020.602567
-
[25]
GanesanHandSutmannG.Modelingsegregatedsolutesinplastically deformedalloysusingcoupledmoleculardynamics-montecarlosim- ulations.Journal of Materials Science & Technology, 213:98–108, 2025.doi:10.1016/j.jmst.2024.06.030
-
[26]
Norman GE, Starikov SV, and Stegailov VV. Atomistic simulation of laser ablation of gold: effect of pressure relaxation.Journal of Experimental and Theoretical Physics, 114:792–800, 2012.doi:10. 1134/S1063776112040115
work page 2012
-
[27]
YaoJ,QiD,LiangH,HeY,YaoY,JiaT,YangY,SunZ,andZhangS. Exploring femtosecond laser ablation by snapshot ultrafast imaging and molecular dynamics simulation.Ultrafast Science, 2022.doi: 10.34133/2022/9754131
-
[28]
ThompsonAP,AktulgaHM,BergerR,BolintineanuDS,BrownWM, CrozierPS,In’tVeldPJ,AKohlmeyer,SGMoore,TDNguyen,etal. LAMMPS-aflexiblesimulationtoolforparticle-basedmaterialsmod- eling at the atomic, meso, and continuum scales.Computer Physics Communications,271:108171,2022.doi:10.1016/j.cpc.2021.108171
-
[29]
RutherfordAMandDuffyDM.Theeffectofelectron–ioninteractions on radiation damage simulations.Journal of Physics: Condensed Matter, 19(49):496201, 2007.doi:10.1088/0953-8984/19/49/496201
-
[30]
Daw MS. Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals.Physical Review B, 29(12):6443, 1984.doi:10.1103/PhysRevB.29.6443
-
[31]
Understanding creep in TiAl alloys on the nanosecond scale by molecular dynamics simulations
Ganesan H, Scheider I, and Cyron CJ. Understanding creep in TiAl alloys on the nanosecond scale by molecular dynamics simulations. Materials & Design, 212:110282, 2021.doi:10.1016/j.matdes.2021. 110282
-
[32]
ChandranA,GanesanH,andCyronCJ. StudyingtheeffectsofNbon high-temperature deformation in TiAl alloys using atomistic simula- tions.Materials&Design,237:112596,2024.doi:10.1016/j.matdes. 2023.112596
-
[33]
Grovenor CRM, Hentzell HTG, and Smith DA. The development of grain structure during growth of metallic films.Acta Metallurgica, 32(5):773–781, 1984.doi:10.1016/0001-6160(84)90150-0
-
[34]
Olbrich M, Pflug T, Wüstefeld C, Motylenko M, Sandfeld S, Rafaja D, and Horn A. Hydrodynamic modeling and time-resolved imaging reflectometryoftheultrafastlaser-inducedablationofathingoldfilm. Optics and Lasers in Engineering, 129:106067, 2020
work page 2020
-
[35]
HirelP. Atomsk: Atoolformanipulatingandconvertingatomicdata files.Computer Physics Communications, 197:212–219, 2015.doi: 10.1016/j.cpc.2015.07.012
-
[36]
IvanovDSandZhigileiLV.Combinedatomistic-continuummodeling ofshort-pulselasermeltinganddisintegrationofmetalfilms.Physical review B, 68(6):064114, 2003.doi:10.1103/PhysRevB.68.064114. Ganesan & Sandfeld:Preprint Page 12 of 13 Delineating the interplay effects of microstructure topology and residual stresses
-
[37]
Iabbaden D, Amodeo J, Fusco C, Garrelieand F, and Colombier J-P. Molecular dynamics simulation of structural evolution in crystalline and amorphous cuzr alloys upon ultrafast laser irradia- tion.Physical Review Materials, 6(12):126001, 2022.doi:10.1103/ PhysRevMaterials.6.126001
work page 2022
-
[38]
Stukowski A. Visualization and analysis of atomistic simulation data withOVITO–theopenVisualizationTool.ModellingandSimulation in Materials Science and Engineering, 18(1):015012, 2009.doi:10. 1088/0965-0393/18/1/015012
work page 2009
-
[39]
Huff M. Residual stresses in deposited thin-film material lay- ers for micro-and nano-systems manufacturing.Micromachines, 13(12):2084, 2022.doi:10.3390/mi13122084
-
[40]
Chen J, Chen W-K, Tang J, and Rentzepis PM. Time-resolved structural dynamics of thin metal films heated with femtosecond op- tical pulses.Proceedings of the National Academy of Sciences, 108(47):18887–18892, 2011.doi:10.1073/pnas.1115237108
-
[41]
WuJ,TangM,ZhaoL,ZhuP,JiangT,ZouX,HongL,LuoS-N,Xi- ang D, and Zhang J. Ultrafast atomic view of laser-induced melting and breathing motion of metallic liquid clusters with mev ultrafast electron diffraction.Proceedings of the National Academy of Sci- ences, 119(4):e2111949119, 2022.doi:10.1073/pnas.2111949119
-
[42]
Liubchenko O, Antonowicz J, Sokolowski-Tinten K, Zalden P, MinikayevR,MilovI,AlbertTJ,BresslerC,ChojnackiM,Dłużewski P, et al. Laser-induced ultrafast structural transformations in thin fe layer revealed by time-resolved x-ray diffraction.arXiv preprint arXiv:2506.18730, 2025.doi:10.48550/arXiv.2506.18730
-
[43]
Zhou S and Shen H. Experimental and numerical studies on micro- bumps without melting of gold films with different thicknesses in- duced by ultrafast laser.Optics Communications, 514:128178, 2022. doi:10.1016/j.optcom.2022.128178. Ganesan & Sandfeld:Preprint Page 13 of 13
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.