Reactive Force Field for P/Sn/I System: Atomistic Insight into the Early Stage of Black Phosphorus and Phosphorene Synthesis Process

Technische Samenvatting: Reactieve Krachtveld voor het P/Sn/I-systeem en Atomistische Inzicht in de Synthese van Zwarte Fosfor

Probleemstelling
Zwarte fosfor (BP) en zijn tweedimensionale tegenhanger, fosforeen, zijn veelbelovende materialen voor de volgende generatie elektronica vanwege hun hoge ladingsmobiliteit en instelbare bandgap. De synthese van hoogwaardige, grootschalige BP-films via Chemische Damptransport (CVT) blijft echter uitdagend. Het proces, dat doorgaans rode fosfor gebruikt met tin (Sn) en jodium (I₂) additieven, is op atomistisch niveau niet volledig begrepen. Onzekerheden blijven bestaan over de specifieke rollen van additieven, de identiteit van stabiele gasfase-soorten, de mechanismen van fosfortransport en de aard van pre-nucleatie intermediairen. Experimentele observatie van deze vroege-fase dynamica is moeilijk vanwege hoge temperaturen en reactieve omgevingen, terwijl standaard Density Functional Theory (DFT) beperkt is tot kleine systemen en statische 0K-configuraties, waardoor het niet in staat is om de dynamische, eindige-temperatuur reactienetwerken te modelleren die vereist zijn voor het begrijpen van groeimechanismen.

Methodologie
Om deze beperkingen aan te pakken, hebben de auteurs een nieuw reactief krachtveld ontwikkeld voor het P/Sn/I-systeem, aangeduid als P/Sn/I-2512. Dit krachtveld is getraind tegen een uitgebreide kwantummechanische (QM) dataset gegenereerd met behulp van DFT (PBE-functionaal met TZ2P basisset en Grimme's D3 dispersiecorrectie). De trainingsset omvatte:

• Binding-dissociatie-energieën voor Sn–Sn, Sn–I, Sn–P, I–P en I–I bindingen.
• Valentiehoek- en torsieprofielen voor diverse moleculaire configuraties.
• Vormingsenthalpieën voor sleutelsoorten (SnI₂, SnI₄, PI₃).
• Toestandsvergelijkingen voor gecondenseerde fase tin ($\alpha$-Sn en $\beta$-Sn) en vormingsenergieën van tinclusters (Sn$_x$, $x=2-6$).

De parameters van het krachtveld werden geoptimaliseerd met behulp van de Covariance Matrix Adaptation Evolutionary Strategy (CMA-ES) binnen de Amsterdam Modeling Suite. Validatie werd uitgevoerd door de ReaxFF-resultaten te vergelijken met DFT-data voor binding-dissociatie, toestandsvergelijkingen en reactie-energieën, waarbij een sterke overeenkomst werd aangetoond.

Vervolgens werden Reactive Molecular Dynamics (ReaxFF-MD) simulaties uitgevoerd met behulp van LAMMPS om de vroege stadia van pre-nucleatieprocessen bij 923 K te onderzoeken. De studie maakte gebruik van het NVT-ensemble met een Nosé-Hoover thermostaat. De simulaties waren gestructureerd om specifieke variabelen te isoleren:

1. Systeemvalidatie: MD-simulaties van Sn/I₂ en P/I₂ systemen bij variërende dichtheden (effectieve drukken van 0,01, 0,1 en 1 atm) om de vorming van sleutelsoorten zoals SnI₂, SnI₄, PI₂ en PI₃ te verifiëren.
2. Rol van Jodium: Drie systemen werden vergeleken (alleen P; P + Sn; P + Sn + I₂) om het effect van jodium op de fosfor-recombinatie te isoleren.
3. Dichtheidseffecten: Het P/Sn/I₂-systeem werd gesimuleerd onder drie dichtheidsregimes (effectieve drukken van 5, 10 en 100 atm) om de overgang van diffuse damp naar gecondenseerde fasen te modelleren, wat de condities nabootst waarin de vorming van BP experimenteel wordt waargenomen te geschieden.

Belangrijkste Resultaten

• Krachtveld Validatie: Het P/Sn/I-2512 krachtveld reproduceerde accuraat de DFT referentiedata voor binding-dissociatie, valentiehoeken, torsiebarrières en de toestandsvergelijkingen voor tinkristallen en clusters. Het voorspelde succesvol de vorming van SnI₂, SnI₄ en PI₃ bij hoge temperaturen.
• Dichtheidsafhankelijkheid in Sn/I₂ en P/I₂: In het Sn/I₂ systeem verschoof een toenemende dichtheid de compositie van elementaire soorten naar complexe coördinatieomgevingen, wat de vorming van SnI₄ en Sn₂I₄ (dimeren van SnI₂) naast SnI₂ bevorderde. In het P/I₂ systeem faciliteerden hogere dichtheden de vorming van PI₃ en polymere PI$_x$ soorten, wat consistent is met de experimentele moeilijkheden bij het isoleren van zuiver PI₃.
• Rol van Jodium in Recombinatie: Simulaties toonden aan dat in de afwezigheid van jodium, fosforatomen grotendeel eveneel atomair bleven of minimale clusters vormden. De toevoeging van I₂ verhoogde de snelheid van P–P recombinatie aanzienlijk, wat leidde tot de vorming van binaire P$_x$I$_y$ motieven en transiënte ternaire Sn$_x$P$_y$I$_z$ verbindingen. Dit bevestigt de rol van jodium bij het oplosbaar maken van fosfor en het faciliteren van transport.
• Pre-nucleatie Pad bij Hoge Dichtheid: Bij lage dichtheden (5 atm) bestond het systeem voornamelijk uit atomaire fosfor en kleine clusters (P₂, P₃). Bij een medium dichtheid (10 atm) begonnen grotere clusters (P₄, P₅) te vormen. Bij een hoge dichtheid (100 atm) condenseerde het systeem snel naar hoog-orde P$_x$ clusters (x = 8–12) en grote ternaire Sn$_x$P$_y$I$_z$ verbindingen.
• Identificatie van Intermediairen: Bij de hoogste dichtheid produceerden de simulaties gecondenseerde, Hittorf-achtige fosforstructuren aan de randen van Sn$_x$P$_y$I$_z$ clusters. Daarnaast werden BP-achtige hexagonale motieven gestabiliseerd door jodium waargenomen, geïdentificeerd als potentiële nucleatiezaden.

Betekenis en Claims
Het artikel claimt een praktisch, atomistisch model te bieden voor het bestuderen van groeicondities en de effecten van additieven in de BP/fosforeen damp-synthese. De primaire bijdragen zijn:

1. Ontwikkeling van een gevalideerd Krachtveld: Het P/Sn/I-2512 krachtveld biedt een betrouwbaar instrument voor het simuleren van complexe reactieve systemen met fosfor, tin en jodium, waarbij de grootte- en tijdschaalbeperkingen van DFT worden overwonnen.
2. Mechanistisch Inzicht: De studie verduidelijkt dat jodium en de systeemdensiteit de controlerende factoren zijn voor fosfor-recombinatie. Het demonstreert dat jodium niet louter een transportagent is, maar actief P–P recombinatie bevordert door de vorming van intermediaire P$_x$I$_y$ en ternaire Sn$_x$P$_y$I$_z$ soorten.
3. Pre-nucleatie Model: De auteurs stellen een specifiek pre-nucleatie pad voor waarbij hoge-densiteit condities, gecombineerd met jodium, de condensatie van fosfor naar Hittorf-achtige intermediairen en BP-achtige hexagonen faciliteren. Deze soorten worden gesuggereerd als precursoren voor de daaropvolgende groei van zwarte fosfor tijdens de afkoelfase.

Het werk stelt geen nieuwe experimentele protocollen voor, maar biedt eerder een theoretisch kader om bestaande experimentele observaties te interpreteren, specifiek de noodzaak van hoge druk (densiteit) en de cruciale rol van jodium in de CVT-synthese van zwarte fosfor.