Pith. sign in

REVIEW 3 major objections 5 minor 43 references

Eliminating nanometer-scale asperities on metallic thin films through plasma modification processes studied by molecular dynamics and AFM

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read For metallic thin films, heavier inert gases flatten nanoscale asperities while etching less material, because momentum exchange in ion–surface collisions depends on gas mass.

desk verdict Solid industry-backed study; the mass-dependent smoothing trend looks real, but the experiment's fixed-power plasma conditions leave the mechanism attribution partly open. read the letter →

arxiv 2502.09031 v1 pith:6PG4HUOE submitted 2025-02-13 cond-mat.mtrl-sci

Tomoyuki Tsuyama , Tatsuki Oyama , Yu Azuma , Haruhisa Ohashi , Masahiro Irie , Ayumi Yamakawa , Shoko Uetake , Takayuki Konno
show 5 more authors
Takahiro Ukai Kohei Ochiai Nobuyuki Iwaoka Atsushi Hashimoto Yoshishige Okuno
This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords moleculardynamicsinertgasionbombardmentsurfaceasperitysmoothingsputteryieldplasmamodificationharddiskmediaatomicforcemicroscopycobaltthinfilm
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to show that the gas species used in a plasma surface-modification process is a control knob for the tradeoff between flattening nanoscale surface asperities and etching away the film: heavier inert gases (krypton, xenon) smooth more efficiently per atom removed, while lighter gases (neon, argon) etch faster. This matters for hard disk drives, because nanometer-scale protrusions on the media can damage read/write heads, and heat-assisted magnetic recording media in particular have dense nano-asperities that an argon treatment fixes only at the cost of thinning the magnetic stack. If the central claim is right, choosing a heavier working gas gives a path to smoother media without sacrificing cobalt layer thickness, and the same mass-ratio rule should help control surface morphology on other metals. The study supports the claim with molecular dynamics simulations of a cobalt slab carrying a semi-spherical asperity, a simplified sputter-yield expression, and atomic force microscopy and X-ray fluorescence measurements on hard disk media treated with argon, krypton, and xenon.

What carries the argument

The carrying mechanism is the mass-dependent momentum exchange in a single ion–surface collision. On impact, a heavy ion undergoes only a small velocity reversal, keeps moving forward, and collides with several slab atoms, redistributing them over the surface; a light ion rebounds, concentrating its impulse on one or a few atoms and knocking them out. The paper formalizes the etching side of this picture with the classical sputter-yield formula, which, with the standard approximation of the dimensionless factor $\alpha \approx 0.15(1+M_s/M_i)$, reduces to $Y \approx \frac{0.45}{\pi^2}\frac{M_s}{M_s+M_i}\frac{V}{E_{sb}}$, so that the yield falls as $M_i$ grows. On the measurement side, the key diagnostic is the AFM kurtosis $R_{\mathrm{Ku}}$, the fourth moment of the surface height distribution divided by the squared second moment, which is high when nanoscale protrusions are dense and decreases as the process removes them.

What would settle it

Measure the sputter yield of a cobalt film for neon, argon, krypton, and xenon at the same incident energy near 100 eV: the identity $Y \propto M_s/(M_s+M_i)$ predicts a strict ordering Ne > Ar > Kr > Xe, so any other ordering would disprove the mechanism. A second check is to repeat the AFM and X-ray fluorescence comparison at a much lower total ion fluence to see whether the heavier-gas benefit in $R_{\mathrm{Ku}}$ per etched nanometer survives outside the simulation's high-flux regime.

Watch

Extended reading notes

Core claim

The paper's central claim is that surface smoothing and material removal in a plasma-based surface modification process are governed by the mass ratio between the working gas ion and the target atom. In molecular dynamics simulations of inert-gas bombardment of a cobalt slab with a semi-spherical asperity, heavier ions (krypton, xenon) shrink the asperity faster while sputtering fewer cobalt atoms, whereas lighter ions (neon, argon) sputter more but flatten less efficiently. Single-ion trajectories show the reason: heavier ions forward-scatter after their first collision and keep transferring momentum deeper into the slab, while lighter ions back-scatter and deliver a sharp local impulse that ejects atoms. The simplified sputter-yield expression $Y \propto M_s/(M_s+M_i)$ reproduces the mass ordering of the etching rate. Atomic force microscopy on hard disk media treated with argon, krypton, and xenon confirms the simulation trend: the kurtosis $R_{\mathrm{Ku}}$ drops from $5.1$ (as-deposited) to $4.0$, $3.3$, and $3.1$, while X-ray fluorescence shows the etching rate falls as gas mass rises.

Load-bearing premise

The load-bearing premise is that mass-dependent momentum exchange seen in the simulation—neutral atoms, one idealized semi-spherical bump, and a flux near $10^{23}\,\mathrm{cm^{-2}s^{-1}}$—is the same mechanism that operates when real ions at a flux near $10^{15}\,\mathrm{cm^{-2}s^{-1}}$ treat a layered FePt/Co disk surface.

Editorial extensions

If this is right

  • Heavier working gases such as xenon let a manufacturer reduce nanoscale asperity density on hard disk media without thinning the underlying cobalt layer as much as an argon process would.
  • Gas identity becomes an independent process parameter: bias voltage changes roughness and thickness together, while switching from argon to krypton to xenon moves the process along a different roughness–thickness tradeoff curve.
  • For applications where the goal is material removal or cleaning, neon and argon remain preferable because their higher sputter yield removes surface material faster.
  • The same mass-ratio scaling should apply to other metallic thin films with comparable surface binding energy, predicting $Y \propto M_s/(M_s+M_i)$ for any target–gas pair.
  • AFM kurtosis, not RMS roughness, is the metric that tracks nanoscale protrusion density; the measurements show RMS with no clear trend while $R_{\mathrm{Ku}}$ falls systematically with gas mass.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the mass-scaling mechanism transfers to production plasmas, the same $M_s/(M_s+M_i)$ identity could be used to select gas mixtures or isotope-enriched gases for a desired balance of smoothing and sputtering on arbitrary metals.
  • The simulation's flux is roughly eight orders of magnitude above the experiment, so an implication left unaddressed is whether cumulative heating or defect buildup at high flux changes the mechanism; this could be probed by comparing single-ion-scattering signatures at low and high repetition rates.
  • Because the AFM images were taken on the carbon overcoat rather than the bare cobalt surface, the larger xenon effect could partly reflect smoothing of the carbon layer instead of the underlying cobalt; measuring the cobalt surface before overcoat deposition would separate the two contributions.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript investigates how the choice of inert-gas species (Ne, Ar, Kr, Xe) affects plasma-assisted smoothing of nanoscale surface asperities on metallic thin films. Molecular dynamics (MD) simulations of ion bombardment of a Co slab with a semi-spherical asperity show that heavier ions reduce the asperity height more efficiently while producing lower etching yields, a trade-off the authors attribute to mass-dependent momentum exchange: heavier ions scatter forward and redistribute material over an extended region, while lighter ions back-scatter and give a stronger impulse to a few surface atoms. This interpretation is supported by single-ion trajectory analyses and a Sigmund-theory scaling argument. AFM measurements on hard disk media processed with Ar, Kr, and Xe show a monotonic decrease in surface kurtosis (RKu) in the order Ar-Kr-Xe, while XRF etch rates on Co films decrease with increasing atomic number; the authors take the experiments to be consistent with the simulations and recommend heavier working gases for smoothing without significant film thinning.

Significance. If it holds, the paper's main contribution is a practical design rule with mechanistic backing: the working-gas species can be used to trade off asperity smoothing against material removal, with Xe as the better smoothing choice and lighter gases appropriate when etching is desired. The strengths of the paper are the focused MD setup at fixed ion energy, which isolates the mass effect; the parameter-free interpretation via the Sigmund yield scaling, which explains the etch-rate ordering without fitting; the reproducible experimental protocol on industrially relevant HDD media; and the explicit acknowledgment of the simulated-vs-experimental flux gap and of the absence of an RMS trend. The experimental validation, however, is weakened by the use of RF bias power as the control variable, which leaves ion energy and flux unmatched across gases, and by the use of a proxy substrate for the etch-rate measurement; these gaps are discussed in the major comments.

major comments (3)
  1. [Description of experimental; AFM measurements; Etching rate (Fig. 8, Table 1)] The headline claim that the experiments are consistent with the simulation results is not fully established because the experiments and simulations are not compared at matched ion conditions. The MD simulations vary the ion kinetic energy explicitly (bias of 50-250 V) and demonstrate the mass ordering at fixed energy, whereas the experiments fix the RF bias power (200 W) and the working pressure (2.0 Pa) without reporting the resulting DC self-bias, ion energy distribution, or ion flux. In a capacitively coupled plasma, both the ion energy and the ion flux are gas-dependent at constant power, so the lower etching rate and the stronger RKu reduction observed for Xe (Fig. 8 and Table 1) could reflect differences in ion energy or flux rather than the mass-dependent momentum transfer identified in the MD simulations. At minimum, the authors should report the plasma parameters for Ar, Kr, and Xe under the process conditions, or soften the claimed consistency accordingly.
  2. [Etching rate (Fig. 8)] The experimental etch-rate measurements are performed on a ~2 nm Co layer deposited on an aluminum substrate, whereas the AFM/RKu characterization and the actual product stack are FePt/Co with a carbon overcoat. Because the central experimental evidence for the etch-rate reduction combines these two different sample geometries, the paper should justify the transferability of the Co etch rate, or of its gas-species ordering, from the Al-substrate proxy to the media stack, for example by adding a measurement on the media or by discussing the relevant differences in morphology, crystallinity, and thermal transport.
  3. [MD simulations; Simulation results (Figs. 3-5)] The continuous-bombardment MD results in Figures 3-5 appear to be single trajectories with no error bars, whereas the single-ion analysis of Figure 6 is averaged over 10 seeds. Because the ion injection is stochastic and the slab is thermostatted at 700 K with thermal fluctuations, the paper should either provide statistics from multiple independent runs, with standard deviations or shaded bands, or explicitly state that the displayed curves are single runs and discuss the expected run-to-run variance. Without this, the quantitative ordering Ne-Ar-Kr-Xe and the distinct-zones claim of Figure 5 do not yet have the statistical support needed to separate the mass effect from run-to-run fluctuations.
minor comments (5)
  1. [Methods: MD simulations] The text states that the middle layer of the thin film is maintained at 700 K, but later states that the 6 ps ion-injection interval is determined by the relaxation time for thermalization to 600 K (Figures S1, S2, and Table S1); these two temperatures should be reconciled or explained.
  2. [AFM measurements (Table 1)] The RMS values in Table 1 do not show the claimed smoothing trend (Ar raises RMS to 0.49 nm relative to 0.42 nm for as-deposited), even though RKu decreases monotonically; the authors note this in the text, but a short explanation of why the fourth-moment metric captures protrusion removal while the second-moment metric does not would help the reader interpret the AFM evidence.
  3. [Abstract / Etching rate] The abstract says that AFM and x-ray fluorescence measurements were performed on hard disk media, but the XRF etch-rate measurements are made on Co deposited on an aluminum substrate; the sample geometry should be stated precisely in the abstract or wherever the etch-rate result is summarized.
  4. [Author contributions] The author-contributions statement credits H.I with designing the experiments, but no author with those initials appears in the author list; this is likely a typographical error and should be corrected.
  5. [Simulation results (Sigmund formula)] The citation Mahneetal. should read Mahne et al.; the accompanying formula Y ≈ (0.45/π^2)·(M_s/(M_s+M_i))·(V/E_sb) is correctly derived from the preceding expressions.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the MD-predicted gas-mass trend is directly simulated and the experimental AFM/XRF trends are measured independently, with no fitted parameters or load-bearing self-citations.

full rationale

The core claim—heavier inert gases reduce nanoscale asperities more efficiently while etching more slowly—is obtained from two independent legs. The MD leg defines asperity height and etched-atom count from explicit atomic trajectories under ZBL/MEAM potentials, and the mass ordering is an emergent simulation outcome, not a fitted input. The experimental leg measures RKu and Co thickness by AFM and XRF on hard disk media after Ar/Kr/Xe treatments; the RKu reduction ordering and etch-rate ordering are raw measurements (Table 1 and Figure 8) and are not constructed from the simulation. The Sigmund sputter-yield expression is used only post hoc to explain the simulated etch-rate ordering and is parameter-free (standard formula, no fitted constant); it does not feed the MD results or the AFM/XRF data. There are no load-bearing self-citations: the cited prior MD works are external groups, and no uniqueness theorem or ansatz is imported from the present authors. The experimental condition of fixed RF bias power rather than fixed ion energy/flux is a potential validity limitation for the claimed experiment-simulation consistency, but it is a confounding-variable concern, not a circularity: the experimental trends are not defined in terms of the simulation outputs, nor vice versa.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The central claim is not obtained by fitting; it is a direct comparison of MD trajectories and AFM/XRF measurements. The free parameters are simulation setup values, and the axioms are standard potentials, standard sputter theory, and the modeling assumption that the idealized bump represents real defects. No new entities are introduced.

free parameters (3)
  • Ion injection interval = 6 ps
    Chosen to give the target thermalization to 600 K; it corresponds to an ion flux about 10^8 times higher than the experiment, and could change the balance of etching vs redistribution.
  • Thermostat temperature = 700 K
    Set to mimic the experimental sputtering temperature; the rate of asperity flattening is temperature-dependent.
  • Asperity geometry = semi-spherical, radius ~2 nm on 2 nm slab
    A single idealized defect; real surfaces have distributions of asperity shapes and sizes, and the choice may influence the relative effectiveness of different ions.
assumptions (5)
  • domain assumption MEAM potential for Co accurately describes surface binding energies
    Invoked in Methods for accurate scattering at <1 keV; from cited refs 28,29.
  • domain assumption ZBL potential describes Co-ion repulsion during high-energy collisions
    Standard for ion-solid MD; cited ref 31.
  • domain assumption Inert gas ions are neutralized by Auger processes immediately at the surface
    Justifies treating ions as neutral atoms; cited refs 32,33.
  • standard math Sigmund sputter-yield formula with Mahne's alpha approximation is valid below 1 keV
    Used to rationalize the etching-rate ordering in Figure 4; from refs 41,42.
  • ad hoc to paper A single artificial semi-spherical asperity represents nanoscale asperities on real media
    Assumed in constructing the MD slab; real surfaces have heterogeneous FePt grains and MgO underlayer roughness.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Eliminating nanometer-scale asperities on metallic thin films through plasma modification processes studied by molecular dynamics and AFM." pith.science (2026). https://pith.science/paper/6PG4HUOE

@misc{pith2026250209031,
  author       = {Pith},
  title        = {Pith review of: Eliminating nanometer-scale asperities on metallic thin films through plasma modification processes studied by molecular dynamics and AFM},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6PG4HUOE}},
  note         = {Machine review of arXiv:2502.09031}
}
read the original abstract

We report the effects of reducing surface asperity size at the nanometer scale on metallic surfaces by plasma-assisted surface modification processes using simulations and experiments. Molecular dynamics (MD) simulations were conducted by irradiating various inert gas ions (Ne, Ar, Kr, and Xe) onto a cobalt slab with nanoscale asperities on the surface. The MD simulations showed that as the atomic number of the inert gas increased the surface asperity size was reduced more efficiently, while the etching rate decreased. The dependencies of the scattering behaviors on the inert gas ions originated from the mass exchange between the working gas ions and the slab atoms. Atomic force microscopy and x-ray fluorescence measurements were performed on hard disk media subjected to the surface modification processes. These measurements experimentally demonstrated that the density of nanoscale asperities was reduced with a lower etching rate as the atomic number of the inert gas increased, consistent with the simulation results. Through this study, we clarified that heavier working gases were more effective in reducing surface asperity size without significantly reducing the thickness of the material, which can contribute to better control of surface morphologies at the nanometer scale.

Figures

Figures reproduced from arXiv: 2502.09031 by the authors.

Figure 1
Figure 1. The initial structures of slabs simulated in this paper are illustrated, from the standard orientation of crystal shape (a), and top (b), and side (c) view. A semi-spherical asperity is artificially set on a slab. They are composed of 16000 of Cobalt (Co) atoms in total. (a) 0 ns (b) 5 ns (c) 10 ns (d) 20 ns Etched Co atom Injected ion [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Snapshots of the etched slab taken at 0 (a), 5 (b), 10 (c), and 20 (d) nanoseconds from the MD simulation, with a 100 V acceleration applied to Ar ions. The size of the protrusion is gradually decreased, and the Co atoms in the slab etched away by the bombardments of inert-gas ions.                          $SSOLHGELDV 9      $VSHULW\KHLJKW QP D E F G 1H $… view at source ↗
Figure 3
Figure 3. The residual height of the asperity on the slab during surface modification process with various inert gases: (a) Ne, (b) Ar, (c) Kr, and (d) Xe. To illustrate the dependence on the applied bias, the MD simulation results under different bias conditions are represented in distinct colors. The definition of the height of the asperity is provided in the main text. 8/11 [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: The number of etched Co atoms as a function of elapsed time for various inert gas species: (a) Ne, (b) Ar, (c) Kr, and (d) Xe.              1H $U .U ;H 7KHQXPEHURIHWFKHG&RDWRPV +HLJKWRIDVSHULW\ QP *DVW\SH [PITH_FULL_IMAGE:figures/…
Figure 5
Figure 5. Figure 5: Plots for the residual height of the asperity and the number of etched atoms during the ion bombardments with various types of inert gases. 9/11 [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: The Z-components of velocities (a), momenta (b), and kinetic energies (c) of a single inert gas ion with applied biases of 100 V are plotted before and after collision with the slab surface.                       …
Figure 7
Figure 7. Figure 7: AFM images of hard disk media measured on COC: (a) as-deposited, (b) etched by Ar, (c) etched by Kr, and (d) etched by Xe, with a substrate RF bias power of 200 W applied to the specimens [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]
Figure 8
Figure 8. Figure 8: The substrate bias-power dependence of etching rate for various ingredient gases. 11/11 [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 37 canonical work pages

  1. [1]

    the internet of things

    Özdemir, V . & Hekim, N. Birth of industry 5.0: Making sense of big data with artificial intelligence, “the internet of things” and next-generation technology policy. OMICS: A J. Integr. Biol. 22, 65–76, DOI: 10.1089/omi.2017.0194 (2018). PMID: 29293405, https://doi.org/10.1089/omi.2017.0194

    work page arXiv 2018

  2. [2]

    & Rydning, J

    Reinsel, D., Gantz, J. & Rydning, J. The evolution of data to l ife-critical. IDC White Pap. 1–25 (2017)

    work page 2017

  3. [3]

    Weller, D. K. et al. A hamr media technology roadmap to an areal density of 4 tb/in 2. IEEE Transactions on Magn. 50, 1–8 (2014)

    work page 2014

  4. [4]

    & Takeo, A

    Roddick, E., Kief, M. & Takeo, A. A new advanced storage resea rch consortium hdd technology roadmap. In 2022 IEEE 33rd Magnetic Recording Conference (TMRC), 1–2, DOI: 10.1109/TMRC56419.2022.9918580 (2022)

    work page arXiv 2022

  5. [5]

    & Santos, T

    Weller, D., Mosendz, O., Parker, G., Pisana, S. & Santos, T. S . L10 feptx-y media for heat-assisted mag- netic recording. physica status solidi (a) 210, 1245–1260, DOI: https://doi.org/10.1002/pssa.201329106 (2013). https://onlinelibrary.wiley.com/doi/pdf/10.1002/pssa.201329106

    work page doi:10.1002/pssa.201329106 2013

  6. [6]

    & Gangopadhyay, S

    Marchon, B., Pitchford, T., Hsia, Y .-T. & Gangopadhyay, S. T he head-disk interface roadmap to an areal density of 4 tbit/in2. Adv. Tribol. 2013, 521086, DOI: https://doi.org/10.1155/2013/521086 (2013). https://onlinelibrary.wiley.com/doi/pdf/10.1155/2013/521086

    work page doi:10.1155/2013/521086 2013

  7. [7]

    C., Tyndall, G

    Strom, B., Lee, S. C., Tyndall, G. W . & Khurshudov, A. Hard dis k drive reliability modeling and failure prediction. In Asia-Pacific Magnetic Recording Conference 2006 , 1–2, DOI: 10.1109/APMRC.2006.365900 (2006)

    work page arXiv 2006

  8. [8]

    & Raman, V

    Khurshudov, A. & Raman, V . Roughness effects on head-disk interface durability and reliability. Tribol. Int. 38, 646–651, DOI: https://doi.org/10.1016/j.triboint.2005.01.042 (2005). Tribology of Information Storage Devices, TISD 200 3

    work page doi:10.1016/j.triboint.2005.01.042 2005

Show all 43 references
  1. [9]

    & Polycarpou, A

    Zhang, Y ., Shakil, A., Humood, M. & Polycarpou, A. A. Finite e lement simulations of sliding contact of the head-disk interface in magnetic storage with lubricant effects. Appl. Surf. Sci. Adv. 6, 100155, DOI: https://doi.org/10.1016/j.apsadv.2021.100155 (2021)

    2021

  2. [10]

    & Polycarpou, A

    Zhang, Y ., Tang, H. & Polycarpou, A. A. High temperature mech anics, friction, wear and adhesion of heat-assisted magnetic recording. Tribol. Lett. 68, 109, DOI: 10.1007/s11249-020-01348-z (2020). 5/11

    2020 doi

  3. [11]

    & Hono, K

    Wang, J., Sepehri-Amin, H., Takahashi, Y ., Ohkubo, T. & Hono, K. Magnetic in-plane components of fept nanogranular film on polycrystalline mgo underlayer for heat-assisted ma gnetic recording media. Acta Materialia 177, 1–8, DOI: https://doi.org/10.1016/j.actamat.2019.07.017 (2019)

    2019 doi

  4. [12]

    Wicht, S. et al. Atomic resolution structure-property relation in highly a nisotropic granular fept-c films with near-stoner-wohlfarth behaviour. J. Appl. Phys. 114, 063906, DOI: 10.1063/1.4817274 (2013). https://pubs.aip.org/aip/jap/article-pdf/doi/10.106 3/1.4817274/15119795/...

    2013 doi

  5. [13]

    Xu, C., V araprasad, B. S. D. C. S., Laughlin, D. E. & Zhu, J.-G. Bias sputtering of granular l10-fept films with hexagonal boron nitride grain boundaries. Sci. Reports 13, 11087, DOI: 10.1038/s41598-023-38106-9 (2023)

    2023 doi

  6. [14]

    & Panjan, M

    Panjan, P ., Drnovšek, A., Gselman, P ., ˇCekada, M. & Panjan, M. Review of growth defects in thin films p repared by pvd techniques. Coatings 10, DOI: 10.3390/coatings10050447 (2020)

    2020 doi

  7. [15]

    & Kim, C.-I

    Joo, Y .-H., Um, D.-S. & Kim, C.-I. Tunable physical properties of al-doped zno thin films by o2 and ar plasma treatments. Mater. Res. Express 8, 126402, DOI: 10.1088/2053-1591/ac3f0a (2021)

    2021 doi

  8. [16]

    Larsen, A. H. et al. The atomic simulation environment—a python library for wor king with atoms. J. Physics: Condens. Matter 29, 273002, DOI: 10.1088/1361-648X/aa680e (2017)

    2017 doi

  9. [17]

    P .et al

    Thompson, A. P .et al. Lammps - a flexible simulation tool for particle-based mater ials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun. 271, 108171, DOI: https://doi.org/10.1016/j.cpc.2021.108171 (2022)

    2022

  10. [18]

    Visualization and analysis of atomistic simu lation data with ovito-the open visualization tool

    Stukowski, A. Visualization and analysis of atomistic simu lation data with ovito-the open visualization tool. Model. Simul. Mater. Sci. Eng. 18, 015012, DOI: 10.1088/0965-0393/18/1/015012 (2009)

    2009 doi

  11. [19]

    & Izumi, F

    Momma, K. & Izumi, F. V esta: a three-dimensional visualizat ion system for electronic and struc- tural analysis. J. Appl. Crystallogr . 41, 653–658, DOI: https://doi.org/10.1107/S0021889808012016 (2008). https://onlinelibrary.wiley.com/doi/pdf/10.1107/S0021889808012016

    2008 doi

  12. [20]

    Sonobe, Y . et al. Thermally stable cgc perpendicular recording media with pt -rich coptcr and thin pt layers. IEEE Transactions on Magn. 38, 2006–2011, DOI: 10.1109/TMAG.2002.801810 (2002)

    2002

  13. [21]

    Papusoi, C. et al. L10 fept films with high tc capping layer for heat assisted mag netic recording (hamr). J. Magn. Magn. Mater. 483, 249–265, DOI: https://doi.org/10.1016/j.jmmm.2019.03.076 (2019)

    2019 doi

  14. [22]

    Berendsen, H. J. C., Postma, J. P . M., van Gunsteren, W . F., Di Nola, A. & Haak, J. R. Molecular dynam- ics with coupling to an external bath. The J. Chem. Phys. 81, 3684–3690, DOI: 10.1063/1.448118 (1984). https://pubs.aip.org/aip/jcp/article-pdf/81/8/3684/18950084/3684_1_o...

    1984 doi

  15. [23]

    & Wang, J.-P

    Zhu, J., Quarterman, P . & Wang, J.-P . Molecular dynamic simu lation study of plasma etching l1 fept me- dia in embedded mask patterning (emp) process. AIP Adv. 7, 056507, DOI: 10.1063/1.4977223 (2017). https://pubs.aip.org/aip/adv/article-pdf/doi/10.1063/1.4977223/12951509/0...

    2017 doi

  16. [24]

    & Wang, J.-P

    Zhu, J. & Wang, J.-P . Study of grain-patterned and highly ord ered l10-fept hamr media us- ing reactive molecular dynamics method. AIP Adv. 14, 025302, DOI: 10.1063/9.0000762 (2024). https://pubs.aip.org/aip/adv/article-pdf/doi/10.1063/9.0000762/19332330/025302_1_9.0000762.pdf

    2024 doi

  17. [25]

    Tinacba, E. J. C., Isobe, M. & Hamaguchi, S. Surface damage fo rmation during atomic layer etching of silicon with chlorine adsorption. J. V ac. Sci. & T echnol. A 39, 042603, DOI: 10.1116/6.0001117 (2021). https://pubs.aip.org/avs/jva/article-pdf/doi/10.1116/6.0001117/1613034...

    2021 doi

  18. [26]

    J., Humbird, D

    Végh, J. J., Humbird, D. & Graves, D. B. Silicon etch by fluoroc arbon and argon plasmas in the pres- ence of fluorocarbon films. J. V ac. Sci. & T echnol. A 23, 1598–1604, DOI: 10.1116/1.2049304 (2005). https://pubs.aip.org/avs/jva/article-pdf/23/6/1598/15653976/1598_1_online.pdf

    2005 doi

  19. [27]

    Athavale, S. D. & Economou, D. J. Molecular dynamics simulat ion of atomic layer etch- ing of silicon. J. V ac. Sci. & T echnol. A 13, 966–971, DOI: 10.1116/1.579659 (1995). https://pubs.aip.org/avs/jva/article-pdf/13/3/966/11804066/966_1_online.pdf

    1995 doi

  20. [28]

    Baskes, M. I. Modified embedded-atom potentials for cubic ma terials and impurities. Phys. Rev. B 46, 2727–2742, DOI: 10.1103/PhysRevB.46.2727 (1992)

    1992 doi

  21. [29]

    & Lee, B.-J

    Dong, W .-P ., Kim, H.-K., Ko, W .-S., Lee, B.-M. & Lee, B.-J. Atomistic modeling of pure co and co-al system. Calphad 38, 7–16, DOI: https://doi.org/10.1016/j.calphad.2012.04.001 (2012)

    2012 doi

  22. [30]

    S., Tucker, O

    Morrissey, L. S., Tucker, O. J., Killen, R. M., Nakhla, S. & Sa vin, D. W . Sputtering of surfaces by ion irradiation: A comparison of molecular dynamics and binary collision appr oximation models to laboratory measurements. J. Appl. Phys. 130, 013302, DOI: 10.1063/5.0051073 (...

    2021 doi

  23. [31]

    Ziegler, J. F. & Biersack, J. P

    "Ziegler, J. F. & Biersack, J. P . "The Stopping and Range of Ions in Matter" , "93–129" (Springer US, "Boston, MA", 1985)

    1985

  24. [32]

    R., Humbird, D

    V ella, J. R., Humbird, D. & Graves, D. B. Molecular dynamics s tudy of silicon atomic layer etching by chlorine gas and argon ions. J. V ac. Sci. & T echnol. B 40, 023205, DOI: 10.1116/6.0001681 (2022). https://pubs.aip.org/avs/jvb/article-pdf/doi/10.1116/6.0001681/19811880/0...

    2022 doi

  25. [33]

    & Ono, K

    Nakazaki, N., Takao, Y ., Eriguchi, K. & Ono, K. Molecular dynamics simulations of si etching in cl- and br-based plasmas: Cl+ and br+ ion incidence in the presence of cl and br neutrals . J. Appl. Phys. 118, 233304, DOI: 10.1063/1.4937449 (2015). https://pubs.aip.org/aip/jap/...

    2015 doi

  26. [34]

    Xu, J. et al. Different etching mechanisms of diamond by oxygen and hydro gen plasma: a reactive molecular dynamics study. The J. Phys. Chem. C 125, 16711–16718, DOI: 10.1021/acs.jpcc.1c03919 (2021)

    2021 doi

  27. [35]

    & Stoneham, A

    Dunn, A., Duffy, D. & Stoneham, A. A molecular dynamics study of diamond exposed to tritium bombardment for fusion applications. Nucl. Instruments Methods Phys. Res. Sect. B: Beam Interact ions with Mater. Atoms 269, 1724–1726, DOI: https://doi.org/10.1016/j.nimb.2010.12.009 (...

    2011 doi

  28. [36]

    Domínguez-Gutiérrez, F. J. et al. Unraveling the plasma-material interface with real time di agnosis of dynamic boron conditioning in extreme tokamak plasmas. Nucl. Fusion 57, 086050, DOI: 10.1088/1741-4326/aa7b17 (2017)

    2017 doi

  29. [37]

    P ., Domínguez-Gutiérrez, F

    Bedoya, F., Allain, J. P ., Domínguez-Gutiérrez, F. J. & Krsti´c, P . S. Effect of deuterium irradiation on graphite boroniz ed in the nstx-u tokamak. Sci. Reports 9, 2435, DOI: 10.1038/s41598-019-38941-9 (2019)

    2019 doi

  30. [38]

    Dwivedi, N. et al. Understanding the role of nitrogen in plasma-assisted surf ace modification of magnetic recording media with and without ultrathin carbon overcoats. Sci. Reports 5, 7772, DOI: 10.1038/srep07772 (2015)

    2015 doi

  31. [39]

    & Bijkerk, F

    Zameshin, A., Medvedev, R., Y akshin, A. & Bijkerk, F. Interface formation in w/si multilayers studied by low energy ion scattering. Thin Solid Films 724, 138569, DOI: https://doi.org/10.1016/j.tsf.2021.138569 (2021)

    2021

  32. [40]

    U., Isobe, M., Karahashi, K

    Tercero, J. U., Isobe, M., Karahashi, K. & Hamaguchi, S. Mole cular dynamics simulations of silicon nitride atomic layer etching with ar, kr, and xe ion irradiations. Jpn. J. Appl. Phys. 63, 07SP03, DOI: 10.35848/1347-4065/ad5d77 (2024)

    2024 doi

  33. [41]

    Theory of sputtering

    Sigmund, P . Theory of sputtering. i. sputtering yield of amo rphous and polycrystalline targets. Phys. Rev. 184, 383–416, DOI: 10.1103/PhysRev.184.383 (1969)

    1969 doi

  34. [42]

    & Panjan, M

    Mahne, N., ˇCekada, M. & Panjan, M. Total and differential sputtering yi elds explored by srim simulations. Coatings 12, DOI: 10.3390/coatings12101541 (2022)

    2022 doi

  35. [43]

    & Eckstein, W

    Behrisch, R. & Eckstein, W . Sputtering by Particle Bombardment (Springer Berlin, Heidelberg, 2007). Author contributions statement T.T and A.H supervised this work. T.T, K.O, and N.I performed the simulations. T.O, Y .A, H.O, H.I, T.U, and A.H designed the experiments. T.T, T...

    2007

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.