Capturing Nuclear Quantum Effects in Hydrogen Diffusion through MoS2 via Machine-Learning-Enhanced Path-Integral Simulations

Technische Samenvatting: Het vastleggen van nucleaire kwantumeffecten in waterstofdiffusie door MoS₂ via machine-learning-verbeterde path-integraal simulaties

Probleemstelling
Waterstoftransport door gelaagde twee-dimensionale (2D) materialen, zoals molybdeen disulfide (MoS₂), is cruciaal voor toepassingen in waterstofopslag, brandstofcellen en isotopenscheiding. Hoewel experimentele studies de interlaagse diffusie en isotoop-selectief transport in materialen zoals MoS₂ en h-BN hebben aangetoond, blijven nauwkeurige theoretische beschrijvingen een uitdaging. Eerdere theoretische onderzoeken vertrouwden grotendeels op klassieke behandelingen van atomaire kernen, waarbij nucleaire kwantumeffecten (NQEs) werden genegeerd. Echter, lichte atomen zoals waterstof vertonen uitgesproken kwantumgedrag—inclusclusief nulpuntsbeweging (ZPE), tunneling en delokalisatie—zelfs bij kamertemperatuur. Bovendien vereist het nauwkeurig modelleren van deze systemen het vastleggen van zowel van der Waals (vdW) interacties als waterstofbruggen, wat geavanceerde uitwisselings-correlatiefunctionals (bijv. r2SCAN+rVV10) in dichtheidsfunctionaaltheorie (DFT) noodzakelijk maakt. De computationele kosten van het combineren van dergelijke hoogwaardige DFT met Path-Integral Molecular Dynamics (PIMD) en enhanced sampling voor grote, structureel complexe systemen (zoals gedraaide bilayers) is prohibitief.

Methodologie
Om deze computationele uitdagingen aan te pakken, hebben de auteurs een machine-learning interatomair potentiaal (MLIP) ontwikkeld op basis van het MACE-framework, getraind op een compacte dataset gegenereerd met r2SCAN+rVV10 DFT. Deze MLIP werd gevalideerd tegenover DFT-referentieberekeningen voor energetische stabiliteit en fononbandstructuren, waarbij een goede overeenkomst werd getoond.

De studie maakte gebruik van een hybride simulatieframework dat de volgende elementen combineert:

1. Well-Tempered Metadynamics (WTMetaD): Geïmplementeerd via PLUMED om vrije-energiebarrières ($\Delta F$) voor interlaagse diffusie te berekenen.
2. Path-Integral Molecular Dynamics (PIMD): Geïmplementeerd via de i-PI package om NQEs expliciet te incorporeren. Kernen werden gerepresenteerd door ring-polymeren met $P=16$ beads voor waterstof (en $P=16$ voor deuterium, om convergentie te garanderen) bij 300 K.
3. Systeemconfiguraties: De studie onderzocht hoog-symmetrische bulk MoS₂ polytypen ($H_h^h$, $H_X^h$, $R_M^h$ en 3R fasen) en gedraaide bilayer MoS₂ (TBL-MoS₂) met draaihoeken van $\theta = 3.89^\circ$, $21.79^\circ$, $38.21^\circ$, en $56.11^\circ$.
4. Isotoopanalyse: Zowel protium (H) als deuterium (D) diffusie werden gesimuleerd om kinetische isotoopeffecten (KIE) te kwantificeren.

Belangrijkste Bijdragen

• Ontwikkeling van een MLIP-PIMD Framework: De auteurs hebben een robuuste workflow opgezet die de kloof overbrugt tussen de nauwkeurigheid van DFT en de lengte/tijdschalen die vereist zijn voor kwantumnucleaire dynamica in complexe 2D-systemen.
• Kwantificering van NQEs in MoS₂: De studie biedt een systematische beoordeling van hoe kwantumfluctuaties diffusielandschappen over verschillende MoS₂ polytypen veranderen.
• Isotoop-selectief Mechanisme: Het werk biedt atomistische inzichten in de mechanismen die isotopenscheiding aandrijven, specif으로 de verschillen in vrije-energiebarrières tussen H en D kwantificeren.
• Moiré-gestuurde Diffusielandschappen: De studie brengt waterstoftransport in gedraaide bilayers in kaart, waarbij wordt onthuld hoe moiré-superroosters ruimtelijk variërende diffusieomgevingen creëren die worden beheerst door lokale stapeling (staggered versus eclipsed) en soliton-netwerken.

Belangrijkste Resultaten

• Impact van NQEs op Diffusiebarrières: De inclusie van NQEs verlaagt de vrije-energiebarrières voor waterstofdiffusie aanzienlijk vergeleken met klassieke simulaties. De reductie varieert van tientallen meV tot bijna 60 meV, afhankelijk van de stapelingsconfiguratie.
  • Voor de meest stabiele $H_h^h$ fase daalt de barrière van 122 meV (klassiek) naar 64 meV (kwantum).
  • Voor de $H_X^h$ fase daalt de barrière van 45 meV naar 18 meV.
  • Deze verzachting van het vrije-energielandschap leidt tot een significante toename van de zelf-diffusiecoëfficiënt ($D$), waarbij de Quantum Correction Factors (QCF) tot wel 9.43 bereiken voor de $H_h^h$ fase.
• Kinetic Isotope Effect (KIE): Er wordt een uitgesproken KIE waargenomen. Hoewel deuterium hogere barrières vertoont dan waterstof in klassieke simulaties, vergroot de inclusie van NQEs dit verschil. De kwantum vrije-energiebarrière voor protium is ongeveer 35 meV lager dan die voor deuterium (60 meV versus 95 meV voor $H_h^h$), wat wijst op een sterke massa-afhankelijke kwantumdiffusie.
• Gedraaide Bilayer Dynamica:
  • Hoog-hoek Draaiingen ($21.79^\circ, 38.21^\circ$): Deze systemen missen uitgebreide hoog-symmetrische domeinen, wat resulteert in diffusielandschappen die worden gedomineerd door intermediaire configuraties. Bij consequent zijn de diffusiebarrières hoger (~120 meV kwantum) en de diffusiviteiten lager (minstens een orde van grootte lager dan geïsoleerde hoog-symmetrische fasen).
  • Laag-hoek Draaiingen ($3.89^\circ, 56.11^\circ$): Deze systemen vormen duidelijke moiré-domeinen (staggered) en knopen/solitonen (eclipsed). Waterstof verblijft bij voorkeur in staggered domeinen waar de tussenlagenafstand kleiner is en de coördinatie gunstig is.
  • Ruimtelijke Variatie: Diffusiebarrières variëren aanzienlijk over het moiré-superrooster, met verschillen tot 80 meV tussen gunstige staggered domeinen en ongunstige soliton/eclipsed regio's. Het $56.11^\circ$ systeem vertoont lagere gemiddelde barrières (150 meV) vergeleken met het $3.89^\circ$ systeem (197 meV), toegeschreven aan de aanwezigheid van metastabiele knopen en gunstige $H_X^h$ domeinen.

Betekenis
Het artikel concludeert dat NQEs een beslissende rol spelen bij waterstoftransport door gelaagd MoS₂, zelfs bij omgevingstemperatuur. Het negeren van deze effecten leidt tot een overschatting van de diffusiebarrières en een onderschatting van de diffusiviteit. De studie toont aan dat de omvang van de kwantumversterking zeer gevoelig is voor de lokale stapelingsomgeving. Bovendien is de expliciete inclusie van kwantumnucleaire dynamica essentieel voor een nauwkeurige beschrijving van isotoop-selectieve transportmechanismen. Ten slotte vestigt het onderzoek de mogelijkheid dat 'twist engineering' in 2D-materialen kan worden gebruikt om diffusielandschappen af te stemmen, waardoor ruimtelijk selectieve kanalen voor waterstofmigratie worden gecreëerd. Het gepresenteerde MLIP-PIMD framework is gepositioneerd als een overdraagbare benadering voor het onderzoeken van kwantumtransport in complexe van der Waals heterostructuren die relevant zijn voor energietechnologieën.