The origin of complex behavior of liquid carbon: an insight from computer simulation

Probleemstelling

Het smelten van grafiet is een langdurig onopgelost probleem in de fysica. Hoewel experimenten al sinds het begin van de 20e eeuw zijn uitgevoerd, vertonen de resultaten voor het smeltpunt van grafiet een enorme spreiding, variërend van $T_m = 4000$ K tot $T_m = 6700$ K. Deze discrepantie wordt toegeschreven aan onzekerheden in de thermodynamische parameters die in eerdere studies zijn gebruikt, zoals:

• De kritieke temperatuur ($T_c$) en het kritieke punt van koolstof, die vaak zijn geschat op basis van aannames die later als onjuist zijn gebleken (bijv. een volumejump bij smelten van 20% in plaats van de experimenteel waargenomen 70%).
• Het negeren van de invloed van de gasfase en sublimatie bij hoge temperaturen.
• De aanname dat de warmtecapaciteit en thermische uitzettingscoëfficiënt van vloeibaar koolstof constant zijn, terwijl deze in de buurt van het kritieke punt sterk kunnen variëren.

De auteurs stellen dat de grote verspreiding in experimentele resultaten mogelijk wordt veroorzaakt door het feit dat het smeltpunt van grafiet dicht bij het kritieke punt van koolstof ligt, wat leidt tot sterke thermodynamische fluctuaties die in eerdere analyses niet correct zijn meegenomen.

Methodologie

De auteurs hebben moleculaire dynamica (MD) simulaties uitgevoerd om de thermodynamische eigenschappen van vloeibaar koolstof te bestuderen en de locatie van het kritieke punt te bepalen.

• Potentiaal: Er is gebruikgemaakt van het GAP-20 (Gaussian Approximation Potential) machine-learning potentieel. Dit potentieel is getraind op DFT-berekeningen en is specifiek ontwikkeld voor de condensatiefasen van koolstof.
• Systeem: Een kubus met 4096 koolstofatomen met periodieke randvoorwaarden.
• Ensemble en Instellingen: Simulaties in het canonieke ensemble (NVT) met een tijdstap van 2,0 fs.
  • Evenwichtstelling: 10 ps.
  • Productie: 20 ps.
• Parameterbereik: Dichtheid ($\rho$) variërend van 0,05 tot 2,0 g/ml en temperatuur ($T$) variërend van 4000 K tot 7000 K.
• Data-analyse: De simulatiedata (druk en interne energie) zijn gefit met een polynoomvergelijking ($X = \sum a_{ij} T^i \rho^j$) om de toestandvergelijking (EoS) te verkrijgen. Hieruit zijn thermodynamische responsfuncties afgeleid, waaronder de isobare warmtecapaciteit ($c_P$), de thermische uitzettingscoëfficiënt ($\alpha_P$) en de isotherme compressibiliteit ($\beta_T$). De kritieke punten en de lijnen van maxima (Widom-lijnen) zijn geïdentificeerd via deze afgeleiden.

Belangrijkste Bijdragen

1. Bepaling van het kritieke punt: Voor het eerst is het kritieke punt van koolstof geschat binnen het kader van het GAP-20 model, wat resulteert in een aanzienlijk lagere kritieke temperatuur dan eerdere schattingen.
2. Karakterisering van Widom-lijnen: De auteurs hebben voor het eerst de locaties van de maxima van thermodynamische responsfuncties (de Widom-lijnen) voor koolstof in kaart gebracht.
3. Verklaring voor experimentele inconsistenties: Het artikel biedt een mechanistische verklaring voor de grote spreiding in experimentele smeltpunten: de aanwezigheid van sterke thermodynamische fluctuaties in de buurt van het kritieke punt, die de metingen van warmtecapaciteit en dichtheid beïnvloeden.

Resultaten

• Kritieke Parameters: De simulaties leveren de volgende waarden op voor het kritieke punt:
  • $T_c = 4235$ K
  • $\rho_c = 0,1887$ g/ml
  • $P_c = 106,61$ bar
    Deze temperatuur is aanzienlijk lager dan de historische schatting van Leider et al. ($T_c \approx 6810$ K) en ligt zeer dicht bij het geschatte smeltpunt van grafiet.
• Thermodynamische Respons: In de buurt van het kritieke punt vertonen de responsfuncties sterke maxima:
  • De isobare warmtecapaciteit ($c_P$) bereikt waarden tot $25 k_B$ per deeltje in de directe omgeving van het kritieke punt, wat de aanname van een constante warmtecapaciteit volledig weerlegt.
  • De thermische uitzettingscoëfficiënt ($\alpha_P$) en de compressibiliteit ($\beta_T$) tonen eveneens sterke pieken.
• Spontane Kristallisatie: Bij $T = 4100$ K en $\rho = 1,55$ g/ml werd spontane kristallisatie van vloeibaar koolstof naar de grafietfase waargenomen. Dit suggereert dat het kritieke punt dicht bij de smeltlijn ligt in het dichtheid-temperatuur vlak.
• Supercritische Vloeistof (SCF): De resultaten suggereren dat wanneer grafiet smelt, het systeem zich mogelijk direct in een supercritische toestand bevindt met sterke fluctuaties, in plaats van een stabiele vloeistof-gas tweefasensituatie.

Betekenis en Discussie

• Invloed op Experimenten: De auteurs concluderen dat de grote spreiding in experimentele smeltpunten waarschijnlijk het gevolg is van het niet meenemen van de sterke fluctuaties rond het kritieke punt. Kleine onzekerheden in druk of temperatuur leiden in dit gebied tot enorme fouten in de berekende dichtheid en warmteabsorptie.
• Beperkingen van het Model: Het GAP-20 potentieel is getraind op condensatiefasen en DFT-data. De auteurs benadrukken dat de exacte numerieke waarden voor het kritieke punt mogelijk niet direct overeenkomen met de werkelijkheid, omdat DFT zelf moeite heeft met het nauwkeurig beschrijven van kritieke fluctuaties en verdunde fasen (gassen).
• Kwalitatieve Validiteit: Hoewel de absolute waarden kunnen verschuiven, is het kwalitatieve beeld waarschijnlijk correct: het kritieke punt ligt waarschijnlijk veel lager dan eerder gedacht (dicht bij het smeltpunt), of de simulaties geven een systeem weer dat qua gedrag correct is maar verschoven is naar lagere temperaturen.
• Toekomstperspectief: Er is behoefte aan nieuwe modellen die tegelijkertijd de eigenschappen van vloeistof, gas en grafiet nauwkeurig beschrijven om de discrepanties tussen simulaties en experimenten volledig op te lossen.

Conclusie: Het artikel toont aan dat de complexe thermodynamische gedragingen van vloeibaar koolstof, en de bijbehorende moeilijkheden bij het bepalen van het smeltpunt, nauw verbonden zijn met de nabijheid van het kritieke punt en de daaruit voortvloeiende sterke fluctuaties.