Campionamento multiscala di un eterogeneo acqua/metallo catalizzatore interfaccia utilizzando la teoria del funzionale della densità e campo di forza dinamica molecolare

I modelli di protocollo, fenomeni termici e quantistici in catalisi eterogenea in fase liquida. È il primo, per quanto ne so, ad incorporare la meccanica quantistica con campionamento completo di un ambiente liquido esplicito. Le configurazioni di molecole liquide generate da questo protocollo rappresentano quelle attese in condizioni di reazione reali e possono essere utilizzate per esplorare fenomeni a livello molecolare che dipendono dalla disposizione spaziale delle molecole.

Le configurazioni di molecole liquide generate da questo protocollo forniscono informazioni sui ruoli che il solvente svolge sulla catalisi eterogenea in fase liquida. Se si sta tentando questo protocollo per la prima volta, è consigliabile assicurarsi innanzitutto di avere accesso e di poter eseguire il software VASP, MCPliQ, VMD e LAMMPS. A dimostrare la procedura sarà Tianjun Xie, uno studente laureato del mio laboratorio.

Dopo aver generato la struttura adsorbato, generare i file di input LAMMPS per una simulazione NPT ed equilibrare il volume della cella utilizzando FFMD. Copiare il file di input LAMMPS nella directory di lavoro. Modificare la variabile di gruppo sulla linea 34 per indicare gli indici di tipo atomico per gli atomi di ossigeno e idrogeno dell'acqua.

Nella variabile di gruppo sulla linea 35, per indicare gli indici di tipo atomico per gli atomi di platino e adsorbato. Sulla riga 17 del file di input, modificare la variabile run step per impostare la durata della simulazione NPT abbastanza a lungo da comprendere un'esecuzione di equilibrazione e un'esecuzione di produzione. Eseguire il software LAMMPS digitando il comando dall'interfaccia della riga di comando, che include le informazioni sul numero di core della CPU da utilizzare e sul nome dell'eseguibile LAMMPS.

La minimizzazione dell'energia affina la configurazione della molecola d'acqua ed è seguita da una simulazione FFMD eseguita a un numero costante di molecole d'acqua, volume e temperatura per portare l'acqua alla temperatura di simulazione. Un'altra simulazione FFMD viene quindi eseguita a un numero costante di molecole d'acqua, pressione e temperatura per determinare l'altezza fisicamente corretta della scatola di simulazione. I file di output verranno utilizzati in un secondo momento.

Dopo la simulazione NPT, tracciare l'altezza della supercella rispetto al tempo. Il punto in cui si stabilizzare a un valore di stato stazionaria è il punto nella simulazione TNP quando la produzione può iniziare. Verificare l'equilibrazione della simulazione del TNP assicurando che le fluttuazioni nell'altezza della supercella siano minime o siano confluita in un valore costante.

Se si verificano grandi fluttuazioni, aprire l'ingresso. equil e ridurre il passo del tempo sulla linea 92 per rigenerare la configurazione della molecola d'acqua ed eseguire il software LAMMPS come prima. Per iniziare, digitare l'interfaccia della riga di comando per eseguire lo script.

Questo script restituisce l'altezza media della supercella dalla parte di esecuzione di produzione della simulazione NPT a un file TXT. Lo script legge la lunghezza della dimensione z della cella a intervalli di 1.000 femtosecondi, che è l'intervallo predefinito per la stampa delle informazioni in LAMMPS. Se si desidera un intervallo di stampa diverso, è possibile modificarlo modificando la riga 20 del get_npt_lz.

py script e la riga 16 dell'ingresso LAMMPS. file equile. Lo script rileva e scarta i primi due nanosecondi di valori lz, in quanto comprendono la parte di equilibrazione della simulazione.

La durata dell'esecuzione dell'equilibrazione può essere modificata modificando la riga 19 del file. I restanti tre nanosecondi comprendono la porzione di produzione, e sono quindi usati per calcolare la lunghezza media della dimensione z. Inoltre, lo script restituisce un altro file TXT che fornisce valori di lz in funzione del passo del tempo e di un file PNG, che traccia gli stessi dati.

Il grafico può essere utilizzato per verificare l'equilibrazione della simulazione del TNP. Per ricostruire la supercella utilizzando l'altezza media determinata in NPT, copiare il file di dati generato in precedenza in una nuova directory di lavoro e rinominarlo come data.myadsorbate. Modificare quindi il nuovo file di dati per modificare zlo in 0.0 e zhi sul valore lz dal valore medio restituito nel file TXT.

Copiare il file di input LAMMPS nella nuova directory di lavoro. Modificare la variabile di gruppo sulla linea 32 per indicare gli indici di tipo atomico per gli atomi di ossigeno e idrogeno dell'acqua e la variabile di gruppo sulla linea 33 per indicare gli indici di tipo atomo per gli atomi di platino e adsorbato. Quindi, sulla riga 16, modificare la variabile runStep in modo che sia abbastanza lunga da comprendere un'esecuzione di equilibrazione e un ciclo di produzione.

Digitare il comando per eseguire LAMMPS nell'interfaccia della riga di comando per eseguire il software LAMMPS. Questo eseguirà una simulazione NVT costante sulle molecole d'acqua e viene generato il file di output chiave. La simulazione NVT comprende una porzione di equilibrazione e una porzione di produzione.

La parte di produzione inizia quando l'energia del sistema viene tracciata rispetto ai livelli di tempo spenti. Per iniziare, aprire lo script di durata del legame idrogeno. Modificate la variabile iniziale effettiva sulla linea 22 per impostare il passo temporale del primo intervallo di tempo.

Modificate la variabile timestep sulla riga 23 per impostare la frequenza con cui i fotogrammi vengono scritti nel file di traiettoria LAMMPS. Modificare N_first e N-ultime variabili nelle righe 24 e 25 per impostare il primo e l'ultimo passo indietro e modificare la variabile nevery sulla riga 26 per impostare se i fotogrammi consecutivi vengono considerati o ignorati. Impostate il numero di linee per sezione cornice del file di traiettoria modificando la variabile della linea del fotogramma sulla linea 27.

Modificare inoltre le linee da 31 a 35, per specificare quali tipi di atomo all'interno dei dati. il file del mioadsorbato appartiene all'adsorbato e quali tipi di atomi appartengono alle molecole d'acqua. Lo script analizza le configurazioni dell'acqua nel ciclo di produzione e determina se alcune molecole d'acqua sono legate all'idrogeno all'adsorbato.

Conta quindi il tempo di simulazione in cui ogni legame idrogeno rimane intatto e riporta queste informazioni come una distribuzione delle vite di legame idrogeno in unità di picosecondi. LAMMPS scrive la configurazione delle molecole d'acqua nel file ogni 1.000 femtosecondi, che è il valore predefinito nel file di input LAMMPS fornito. Rileva e scarta le prime due configurazioni del valore di due nanosecondi nel file in quanto comprendono la parte di equilibrazione della simulazione e utilizza i restanti tre nanosecondi per calcolare le vite del legame idrogeno.

Per eseguire lo script, digitare nell'interfaccia della riga di comando. Genera quindi un file DAT. Tracciare i dati nel file per visualizzare la distribuzione delle vite di legame idrogeno che si sono verificate durante la simulazione NVT.

Per determinare l'incremento di tempo da utilizzare per l'intervallo di campionamento del tempo, utilizzare un incremento di tempo maggiore o uguale alla durata massima del legame idrogeno. Determinare il numero di configurazioni dall'esecuzione di produzione della traiettoria FFMD NVT per campionare in modo che il tempo minimo tra le configurazioni sia uguale o maggiore dell'intervallo di campionamento temporale identificato in precedenza. Sullo script di estrazione dei frame scritto in precedenza, modificare il valore predefinito per il numero di frame variabili sulla riga 21 per specificare il numero di configurazioni da estrarre.

Per eseguire lo script, digitare il nome dello script nell'interfaccia della riga di comando. Verrà prodotto un elenco di tempi di simulazione corrispondenti alle configurazioni che devono essere estratte dal file di simulazione NVT. Queste configurazioni possono essere utilizzate come strutture di partenza nelle simulazioni AIMD o QM.

In questa procedura, FFMD è stato utilizzato per generare la configurazione iniziale delle molecole d'acqua. La simulazione AIMD mostra che una molecola d'acqua originariamente incollata all'idrogeno ad un adsorbato di alcol di zucchero su una superficie platino-111 astrae l'idrogeno dall'adsorbato alcolico e deposita un secondo idrogeno sulla superficie del platino-111. Le strutture delle molecole di acqua liquida dipendono dalle impostazioni di input.

Impostare questi in modo improprio può avere influenze indesiderate sulle strutture dell'acqua. In questa figura, il lato sinistro è la struttura iniziale per un'esecuzione FFMD. E il lato destro è entro un picosecondo dall'inizio della simulazione.

La simulazione FFMD esplode a causa di impostazioni di forza non fisicamente grandi, causando l'allontanamento delle molecole d'acqua dalla superficie. Le configurazioni possono essere utilizzate nella meccanica quantistica, o simulazione QMM, o potrebbero essere utilizzate per analizzare statistiche relative alle posizioni spaziali delle molecole. Questa tecnica apre la strada ai ricercatori per esplorare i ruoli che gli ambienti di reazione liquida hanno sulla catalisi generando configurazioni effettive di molecole liquide alle interfacce catalitiche.