2025-11-22T23:31:16.072362

Stainless steel in an electronically excited state

Medvedev
Understanding the non-equilibrium behavior of stainless steel under extreme electronic excitation remains a critical challenge for laser processing and radiation science. We employ a hybrid framework integrating density-functional tight binding, transport Monte Carlo, and Boltzmann equations to model austenitic stainless steel (Fe$_{0.5875}$Cr$_{0.25}$Mn$_{0.09}$Ni$_{0.07}$C$_{0.0025}$) under ultrafast irradiation. The developed approach uniquely bridges atomic-scale electronic dynamics and mesoscale material responses, enabling the quantitative mapping of electron-temperature-dependent properties (electronic heat capacity, thermal conductivity, and electron-phonon coupling) up to the electronic temperatures Te~25,000 K. Two distinct lattice disordering mechanisms are identified: nonthermal melting at Te~10,000 K (the dose ~1.4 eV/atom), where the lattice collapses on sub-picosecond timescales without atomic heating driven by electronic excitation modifying the interatomic potential; and thermal melting (at ~0.45 eV/atom), induced by electron-phonon coupling on picosecond timescales. The derived parameters enable predictive modeling of stainless steel under extreme conditions, with implications for laser machining and radiation-resistant material design.
academic

L'acier inoxydable dans un état électroniquement excité

Informations fondamentales

  • ID de l'article : 2504.19798
  • Titre : Stainless steel in an electronically excited state
  • Auteur : Nikita Medvedev (Institut de physique et Institut de physique des plasmas, Académie des sciences tchèque)
  • Classification : cond-mat.mtrl-sci (Physique de la matière condensée - Science des matériaux)
  • Date de publication : Avril 2025 (prépublication arXiv)
  • Lien de l'article : https://arxiv.org/abs/2504.19798

Résumé

Cette étude aborde le défi clé du comportement hors équilibre de l'acier inoxydable sous excitation électronique extrême dans le traitement au laser et la science des radiations. Elle utilise un cadre hybride intégrant la théorie de la fonctionnelle de la densité liaisons fortes (DFTB), le transport par Monte-Carlo et l'équation de Boltzmann pour modéliser le comportement de l'acier inoxydable austénitique (Fe₀.₅₈₇₅Cr₀.₂₅Mn₀.₀₉Ni₀.₀₇C₀.₀₀₂₅) sous irradiation ultrarapide. Cette approche relie de manière unique la dynamique électronique à l'échelle atomique et la réponse des matériaux à l'échelle mésoscopique, réalisant une cartographie quantitative des propriétés dépendantes de la température électronique (capacité thermique électronique, conductivité thermique et couplage électron-phonon) jusqu'à Te25 000 K. L'étude identifie deux mécanismes distincts de désordre du réseau : (1) la fusion non thermique, survenant à Te11 500 K (dose~1,9 eV/atome), où le réseau s'effondre à l'échelle de la sous-picoseconde sans chauffage atomique, pilotée par la modification du potentiel interatomique par excitation électronique ; (2) la fusion thermique (~0,45 eV/atome), induite par le couplage électron-phonon à l'échelle de la picoseconde.

Contexte et motivation de la recherche

Définition du problème

  1. Problème central : Comprendre les mécanismes de réponse hors équilibre de l'acier inoxydable austénitique sous irradiation laser ultrarapide et dans les environnements de rayonnement extrême, en particulier le comportement des matériaux dans les états d'excitation électronique hautement énergétiques.
  2. Importance :
    • Applications généralisées : L'acier inoxydable possède des applications irremplaçables en médecine (implants, instruments chirurgicaux), génie civil (environnements corrosifs, aérospatiale), chimie pétrolière (tuyauteries, réservoirs) et environnements sévères de rayonnement (réacteurs nucléaires, accélérateurs de particules, détecteurs de faisceaux laser à électrons libres)
    • Besoins pratiques : Le traitement au laser, l'évaluation des dommages de rayonnement et la conception de matériaux radiorésistants nécessitent une compréhension précise de la réponse des matériaux dans les conditions extrêmes
  3. Limitations des méthodes existantes :
    • Modèle à deux températures : La dynamique moléculaire à deux températures traditionnelle ne décrit que les effets thermiques, ignorant l'influence directe de l'excitation électronique sur le potentiel interatomique
    • Paramètres manquants : Absence de paramètres de capacité thermique électronique, conductivité thermique et couplage électron-phonon pour l'acier inoxydable à haute température électronique (Te > 10 000 K)
    • Effets non thermiques : Le phénomène de fusion non thermique dans les métaux est rarement étudié ; l'existence de ce mécanisme dans l'acier inoxydable reste incertaine
  4. Motivation de la recherche :
    • Établir un modèle multi-échelle capable de décrire simultanément les effets thermiques et non thermiques
    • Évaluer quantitativement les propriétés thermodynamiques et de transport de l'acier inoxydable à température électronique extrême
    • Déterminer l'existence d'un mécanisme de fusion non thermique dans l'acier inoxydable et son seuil

Contributions principales

  1. Cadre de modélisation multi-échelle hybride : Développement de la boîte à outils XTANT-3 intégrant DFTB, transport par Monte-Carlo et équation de Boltzmann, capable de traiter simultanément les cascades électroniques, la thermalisation électronique et la dynamique atomique
  2. Base de données de paramètres dépendant de la température électronique : Premier calcul systématique pour l'acier inoxydable austénitique à température électronique extrême (jusqu'à Te~25 000 K) :
    • Capacité thermique électronique
    • Conductivité thermique électronique (contributions électron-phonon et électron-électron)
    • Paramètres de couplage électron-phonon
  3. Découverte de mécanismes de dommages doubles : Identification de deux mécanismes de fusion distincts dans l'acier inoxydable :
    • Fusion non thermique : seuil de dose ~1,9 eV/atome, échelle de temps sous-picoseconde, sans chauffage atomique
    • Fusion thermique : seuil de dose ~0,45 eV/atome, échelle de temps picoseconde, via couplage électron-phonon
  4. Réponse spécifique aux éléments : Découverte que le sous-système manganèse est particulièrement sensible à l'excitation électronique, présentant un comportement d'accélération non thermique
  5. Capacité de modélisation prédictive : Fourniture de courbes de flux de seuil de dommages dépendant de l'énergie des photons, pouvant guider les applications expérimentales et d'ingénierie

Détails méthodologiques

Définition de la tâche

Entrées : Paramètres d'impulsion laser (énergie des photons, durée d'impulsion, dose absorbée) et état initial du matériau Sorties :

  • Propriétés thermodynamiques et de transport dépendant de la température électronique
  • Évolution temporelle résolue des températures atomique et électronique
  • Changements de structure du matériau (fusion, désordre) Contraintes : Applicable aux processus d'irradiation ultrarapide (sous-picoseconde à picoseconde)

Architecture du modèle

L'architecture globale utilise la boîte à outils de simulation hybride XTANT-3, comprenant trois modules couplés :

1. Module de transport électronique

Électrons rapides (E > 10 eV) : Monte-Carlo piloté par événements

  • Section efficace d'absorption optique : base de données EPICS2023
  • Ionisation par collision : forme de fonction diélectrique complexe (approximation monopole)
  • Diffusion quasi-élastique : section efficace de Molière modifiée
  • Décroissance Auger : incluant la dynamique des trous cœur-coquille

Électrons lents (E < 10 eV) : Équation de Boltzmann

  • Thermalisation électron-électron : approximation du temps de relaxation
  • Diffusion électron-phonon : méthode de couplage dynamique

2. Module de structure électronique

Théorie de la fonctionnelle de la densité liaisons fortes (DFTB) :

  • Paramétrisation : ensemble de paramètres liaisons fortes PTBP
  • Base : combinaison linéaire d'orbitales atomiques sp³d⁵ (LCAO)
  • Grille d'espace k : grille de Monkhorst-Pack 7×7×7
  • Sorties : niveaux d'énergie transitoires (structure de bandes), fonctions d'onde, forces atomiques

3. Module de dynamique atomique

Dynamique moléculaire (MD) :

  • Forces interatomiques : calculées à partir du gradient du hamiltonien liaisons fortes
  • Énergie potentielle : dépendante de la population électronique transitoire
  • Algorithme d'intégration : algorithme d'ordre 4 de Martyna-Tuckerman
  • Pas de temps : ≤1 fs

4. Système de simulation

  • Composition chimique : Fe₀.₅₈₇₅Cr₀.₂₅Mn₀.₀₉Ni₀.₀₇C₀.₀₀₂₅
  • Boîte de simulation : 400 atomes (235 Fe, 100 Cr, 36 Mn, 28 Ni, 1 C)
  • Structure cristalline : cubique à faces centrées (fcc)
  • Dimensions de la boîte : 18,45×18,4×14,76 Ų
  • Densité : 7,3 g/cm³ (valeur expérimentale 7,5-7,9 g/cm³)

Points d'innovation technique

1. Calcul de la capacité thermique électronique

Utilisation de la méthode statistique des niveaux d'énergie transitoires :

Ce(Te)=1V0∑i∂fe(Ei)∂Te(Ei−μ(Te))C_e(T_e) = \frac{1}{V_0} \sum_i \frac{\partial f_e(E_i)}{\partial T_e}(E_i - \mu(T_e))

où fe(Ei)f_e(E_i) est la distribution de Fermi-Dirac et μ(Te)\mu(T_e) est le potentiel chimique électronique

2. Calcul de la conductivité thermique électronique

Combinaison par la règle de Matthiessen : κtot(Te)=(1κe−a(Te)+1κe−e(Te))−1\kappa_{tot}(T_e) = \left(\frac{1}{\kappa_{e-a}(T_e)} + \frac{1}{\kappa_{e-e}(T_e)}\right)^{-1}

Contribution électron-phonon : Forme de Kubo-Greenwood κe−a(Te)=L22−TeL122L11\kappa_{e-a}(T_e) = L_{22} - T_e \frac{L_{12}^2}{L_{11}}

Coefficients d'Onsager : Lij=−(−1)i+jV0me∑kdfdEk(Ek−μ)i+j+2∣⟨k∣p∣k′⟩∣2L_{ij} = -\frac{(-1)^{i+j}}{V_0 m_e} \sum_k \frac{df}{dE_k}(E_k - \mu)^{i+j+2}|\langle k|p|k'\rangle|^2

Contribution électron-électron : Basée sur les sections efficaces de diffusion Monte-Carlo

3. Calcul du couplage électron-phonon

Méthode de couplage dynamique non perturbatif : G(Te,Ta)=1V0(Te−Ta)∑i,jEjIe−aijG(T_e, T_a) = \frac{1}{V_0(T_e - T_a)} \sum_{i,j} E_j I_{e-a}^{ij}

Intégrale de diffusion :

undefined