Ich versuche, den Phasenübergang für Gold und Silizium getrennt mit LAMMPS zu simulieren. Ich habe den Schmelzpunkt für Gold richtig mit dem folgenden Code.
units metal
atom_style atomic
boundary p p p
variable a equal 4.0782
lattice fcc 4.0782
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box
mass 1 196.97
pair_style eam
pair_coeff * * Au_u3.eam
minimize 1.0e-8 1.0e-8 1000 100000
min_style cg
timestep 0.001
velocity all create 300.0 454883 mom yes rot yes dist gaussian
thermo 50000
thermo_style custom step pe ke etotal temp vol press density atoms
fix 1 all press/berendsen iso 0.0 0.0 100.0
fix 2 all nvt temp 300.00 2400.00 1.0
run 10000000
Wenn ich diesen Code jedoch auf Silizium anwende, erhalte ich nicht die richtigen Ergebnisse. Ich habe das Gefühl, dass ich den Code und die Physik hinter dem Erhalten eines Phasenübergangsdiagramms mit Molekulardynamik nicht verstehe. Ich schätze, meine Bitte hier ist, ob Sie Empfehlungen für einige Artikel oder Bücher haben, die ich lesen sollte, um die Simulation durchführen zu können (Phasenübergang in der Molekulardynamik). Ich habe das Buch "Molecular Dynamics verstehen von Dan Frenkel und Berend Smit" durchgeblättert, aber ich habe immer noch das Gefühl, dass mir etwas fehlt.
bearbeitet: Die Ausgabe ist unten für Gold und Silizium dargestellt.
Metalle wie Gold unterliegen Phasenübergängen mit einer sehr begrenzten Keimbildungsbarriere, während dies bei Silizium nicht der Fall ist.
Wenn Sie beispielsweise periodische Randbedingungen (keine Oberflächen) und keine anderen Defekte haben, kann Silizium sehr leicht überhitzt werden (zumindest auf der typischen Zeitskala einer MD-Simulation).
Wenn Sie den Schmelzpunkt messen möchten, ist es richtig, zunächst zwei Phasen - kristallin und amorph - in Kontakt miteinander zu erzeugen und dann herauszufinden, bei welcher Temperatur die beiden im Gleichgewicht sind. (Bei niedrigeren Temperaturen wächst die Kristallphase, bei höheren Temperaturen wächst die amorphe Phase).
Wenn es jedoch Ihr Ziel ist, das Nukleationsereignis zu simulieren, dann ist dies ein sehr schwieriges Problem und ein sehr aktives Forschungsgebiet. Sie könnten einfach zu einer höheren Temperatur gehen oder Ihre Zellgröße korrigieren (entfernen Sie die Presse / Berendsen-Korrektur) - das führt zu einem Phasenübergang, aber es ist nicht unbedingt ein realistischer Übergang. Dafür gibt es weitaus kompliziertere Verfahren, bei denen Rare-Event-Sampling-Methoden wie Metadynamik [ 1 ] oder Seeding-Methoden [ 2 ] zum Einsatz kommen.
Emil
valerio
Zitrone
phyM