Summary
I calcoli eseguiti dal pacchetto di simulazione initio di Vienna Ab possono essere utilizzati per identificare le proprietà elettroniche intrinseche dei materiali su nanoscala e prevedere i potenziali fotocatalizzatori di scissione dell'acqua.
Abstract
Gli strumenti computazionali basati sulla teoria della densità-funzionale (DFT) consentono l'esplorazione dei composti su nanoscala qualitativamente nuovi e sperimentalmente raggiungibili per un'applicazione mirata. Le simulazioni teoriche forniscono una profonda comprensione delle proprietà elettroniche intrinseche dei materiali funzionali. L'obiettivo di questo protocollo è quello di cercare i candidati fotocatalizzatori per dissezione computazionale. Le applicazioni fotocatalitiche richiedono spazi di banda adeguati, posizioni appropriate del bordo della banda rispetto ai potenziali di redox. I funzionali ibridi possono fornire valori precisi di queste proprietà, ma sono costosi dal punto di vista computazionale, mentre i risultati a livello di funzionalità per Il valore Di S.C.-Bbe (PBE) potrebbero essere efficaci per suggerire strategie per l'ingegneria della struttura delle bande campo elettrico e tensione con l'obiettivo di migliorare le prestazioni fotocatalitiche. Per illustrare questo concetto, nel presente manoscritto, lo strumento di simulazione basato su DFT VASP viene utilizzato per studiare l'allineamento delle bande di nanocompositi in combinazioni di nanotubi e nanoribbons nello stato del suolo. Per affrontare la durata dei fori e degli elettroni fotogenerati nello stato eccitato, sono necessari calcoli di dinamica non adiabatica.
Introduction
La domanda mondiale di energia pulita e sostenibile ha stimolato la ricerca di materiali promettenti per ridurre la dipendenza dalle risorse petrolifere finite. Le simulazioni sono più efficienti ed economiche degli esperimenti per accelerare la ricerca di nuovi materiali funzionali1. La progettazione dei materiali da una prospettiva teorica2,3,4 è ora sempre più popolare grazie ai rapidi progressi nelle risorse computazionali e agli sviluppi teorici, rendendo le simulazioni computazionali più affidabili5 . I calcoli della teoria funzionale della densità (DFT) implementati in molti codici stanno diventando più robusti e producono risultati riproducibili6.
Il Vienna Ab initio Simulation Package (VASP)7 presenta uno dei codici DFT più promettenti per prevedere le proprietà molecolari e cristalline e sono stati pubblicati più di 40.000 studi che fanno uso di questo codice. La maggior parte del lavoro viene eseguita al livello funzionale Perdew-Burke-Ernzerhof (PBE)8, che sottovaluta le dimensioni gap della banda, ma cattura le tendenze essenziali nell'allineamento delle bande e negli offset delle bande3. Questo protocollo ha lo scopo di delineare i dettagli dello studio dei profili di bordo delle bande e delle bande di materiali su nanoscala per l'energia pulita e rinnovabile utilizzando questo strumento computazionale. Altri esempi di utilizzo di VASP sono disponibili presso https://www.vasp.at.
Questo rapporto presenta lo screening computazionale di eterostrutture vdW unidimensionali (1D) con allineamenti di banda di tipo II9 per un'applicazione promettente nella suddivisione dell'acqua fotocatalitica4 . In particolare, i nanoribbon (NR) incapsulati all'interno di nanotubi (NT) vengono esaminati come esempio10. Per risolvere le interazioni non covalenti, le correzioni vdW mediante il metodo DFT-D3 sono incluse11. I calcoli DFT nei passaggi 1.2, 2.2, 3.2, 3.5.2 e section 4 di VASP vengono eseguiti utilizzando uno script PBXS (Portable Batch System) dai computer di ricerca ad alte prestazioni nel sistema CenTOS. Un esempio di script PBS è mostrato nel materiale supplementare. Il postprocessing dei dati da parte del software P4VASP nel passaggio 3.3 e la figura del software xmgrace nel passaggio 3.4 vengono effettuati su un computer locale (laptop o desktop) nel sistema Ubuntu.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. Ottimizzare la struttura atomica.
- Preparare quattro file di input per il calcolo del rilassamento della struttura da PARTE di VASP: INCAR, POSCAR, POTCAR e KPOINTS.
NOTA: nel file INCAR sono presenti parametri specificati che definiscono il calcolo. La riga "EDIFFG - 0,02" nel file INCAR indica che tutti gli atomi sono rilassati fino a quando la forza su ogni atomo non è <0.02 eV/ . Il file POSCAR contiene le informazioni sulla geometria atomica. I parametri primi del reticolo nel file POSCAR possono essere scelti da3 o riferimenti sperimentali teorici12,13. Il file KPOINTS definisce la mesh di punto k e POTCAR è il file pseudopotenziale. L'ordine dei tipi di atomo in POSCAR deve essere lo stesso di quello in POTCAR. Esempi di file di input per il rilassamento della struttura sono mostrati nei Materiali supplementari (ad eccezione del file pseudopotenziale, che necessita di una licenza da VASP).- Generare la struttura iniziale dei nanoribbonri di nitrato di boro (BN) (NR) per "POSCAR".
- Scaricare il file POSCAR per l'unità bulk BN da https://materialsproject.org.
- Utilizzare v2xsf per convertire il file POSCAR in un file in formato xsf che può essere letto da xcrysden. Digitare v2xsf POSCAR sul terminale nel sistema Ubuntu per ottenere "POSCAR.xsf.gz". Digitare gunzip POSCAR.xsf.gz e generare il file POSCAR.xsf.
- Utilizzare xcrysden per costruire la supercella BN.
- Digitare xcrysden --xsf POSCAR.xsf sul terminale nel sistema Ubuntu. Selezionare il menu Modifica/Numero di unità disegnate ed estendere la cella nelle direzioni X e Y.
- Selezionare il menu File/Salva struttura XSF per esportare la struttura della supercella, denominata "supercella".
NOTA: il nome della struttura è una definizione arbitraria.
- Utilizzare xmakemol per aprire la supercella. Tipo xmakemol -f supercella sul terminale nel sistema Ubuntu. Selezionare il menu Modifica/Visibile. Fare clic su Attiva/disattiva per eliminare gli atomi all'interno della regione e tagliare il NR alla larghezza e chiralità desiderate.
- Generare la struttura iniziale del nanotubi BN (NT) per POSCAR. Scarica "NanotubeModeler" da http://www.jcrystal.com/products. Aprire NanotubeModeler.exe nel sistema Windows. Selezionare il menu Selezionare il tipo/B-N e specificare la chiralità. Selezionate il menu File/Salva tabella XY per esportare la struttura.
- Generare la struttura iniziale del nanocomposito incapsulando il NR (dal punto 1.1.1) all'interno del NT (dal passaggio 1.1.2).
NOTA: l'incapsulamento può essere terminato regolando le coordinate cartesiane del NR e NT10,14,15. - Utilizzare il software vmd per controllare la struttura atomica prima di inviare il processo di calcolo.
- Digitare vmd sul terminale nel sistema Ubuntu. Nella finestra principale di vmd aperta, selezionare il menu File/Nuova molecola e trovare il file POSCAR tramite la finestra Sfoglia. Caricare POSCAR digitando VASP_POSCAR.
- Visualizzare la struttura in stili diversi nella finestra Rappresentazioni grafiche/Metodo di disegno.
NOTA: Ad esempio, una volta scelto il CPK, ogni atomo (obbligazione) è rappresentato da una sfera (bastone). La guida all'installazione e l'esercitazione completa di vmd sono disponibili presso http://www.ks.uiuc.edu/Research/vmd.
- Generare la struttura iniziale dei nanoribbonri di nitrato di boro (BN) (NR) per "POSCAR".
- Digitare qsub job.pbs sul terminale nel sistema Linux per inviare il processo al cluster di computer.
NOTA: "job.pbs" rappresenta il nome dello script PBS. Il nome dello script PBS è una definizione arbitraria. I quattro file di input insieme allo script PBS devono trovarsi nella directory di lavoro. Il comando qsub job.pbs verrà utilizzato nei passaggi 2.2, 3.2, 3.5.2 e sezione 4. Un esempio di script PBS è disponibile nel file di codifica supplementare. Al termine del lavoro inviato, se " ha raggiunto laprecisione richiesta - l'interruzione della minimizzazione dell'energia strutturale" appare alla fine del log di output, si ottiene il risultato convergente. Il file CONTCAR risultante verrà utilizzato come file di input POSCAR nei passaggi 2.1, 3.1, 3.5.1, 3.5.3.1, 4.1.1, 4.1.4 e sezione 4.2.
2. Calcolare l'energia di incapsulamento.
- Tipo mkdir nanocomposito isolated-nanoribbon isolated-nanotube per creare tre cartelle per il nanocomposito, il NR, e il NT su un terminale in un sistema Linux. Preparare uno script PBS "job.pbs" e quattro file di input INCAR, POSCAR, POTCAR e KPOINTS per il calcolo dell'energia in ogni cartella.
NOTA: il file di input POSCAR è il file denominato CONTCAR con la struttura rilassata dal passaggio 1. Esempi di file di input sono forniti in Materiali supplementari (ad eccezione di POTCAR). - Vai a ogni cartella e digita qsub job.pbs sul terminale nel sistema Linux.
NOTA: i tre processi inviati eseguiranno i calcoli statici dell'energia autocoerente rispettivamente per il NR nanocomposito, isolato e isolato. - Estrarre l'energia totale dal file OUTCAR per ogni sistema dopo aver completato i calcoli statici auto-coerenti. Digitare grep "free energy TOTEN" ./nanocomposite/OUTCAR - tail -n 1, grep "free energy TOTEN" ./isolated-nanoribbon/OUTCAR , tail -n 1, e grep "free energy TOTEN" ./isolated-nanotube/OUTCAR - tail -n 1. Definire i tre valori visualizzati rispettivamente come ENT, NR, ENRed ENT. Calcolare l'energia di incapsulamento per angstrom: EL - (ENT-NR - ENT -ENR)/L14,15.
NOTA: La direzione periodica in ogni sistema è lungo l'asse z e L è la costante del reticolo della cella unitaria lungo l'asse . Sono necessari calcoli di prova della dipendenza energetica dall'energia di taglio dell'onda del piano e dalla mesh del punto k. L'energia di incapsulamento può essere utilizzata come stima della stabilità energetica del nanocomposito.
3. Estrarre le proprietà elettroniche dalla struttura della banda.
- Preparare uno script PBS "job.pbs" e sei file di input: INCAR, POSCAR, POTCAR, KPOINTS, CHGCAR e CHG per il calcolo della banda. Impostare ICHARGE - 11 in INCAR.
NOTA: i file CHGCAR e CHG preconvergenti provengono dai calcoli statici auto-coerenti nel passaggio 2.2. L'analisi della banda è a livello PBE. Il campionamento di punti k nel file KPOINTS è in modalità linea. Esempi di file di input per questo passaggio sono disponibili nei Materiali supplementari (ad eccezione di POTCAR). - Digitare qsub job.pbs sul terminale nel sistema Linux per inviare il processo.
- Utilizzare P4VASP per generare la banda proiettata.
- Caricare "vasprun.xml" digitando p4v vasprun.xml sul terminale nel sistema Ubuntu.
NOTA: "p4v" viene utilizzato per avviare P4VASP. Il file "vasprun.xml" deve trovarsi nella directory di lavoro. - Selezionate il menu Eopzionele/Locale DOS-bande controllo e quindi Seleziona/Bande.
- Specificare i numeri atomici di NT nella sezione Atom selection. Ottenere il numero atomico puntando agli atomi corrispondenti usando vmd come indicato nel passaggio 1.1.4. Specificare il colore, il tipo e le dimensioni del simbolo per la struttura della banda proiettata tramite il menu Simbolo e Dimensione simbolo. Premere il menu Aggiungi nuova riga.
NOTA: Il grafico mostrerà la struttura della banda con i contributi del NT. - Ripetere la stessa procedura dopo il passaggio 3.3.2.1 per ottenere la banda proiettata con i contributi della NR.
- Specificare i numeri atomici di NT nella sezione Atom selection. Ottenere il numero atomico puntando agli atomi corrispondenti usando vmd come indicato nel passaggio 1.1.4. Specificare il colore, il tipo e le dimensioni del simbolo per la struttura della banda proiettata tramite il menu Simbolo e Dimensione simbolo. Premere il menu Aggiungi nuova riga.
- Selezionare il menu Grafico/Esporta. Esportare il grafico in un file in formato agr (ad esempio, come "11-4.agr").
NOTA: i dati di output delle bande proiettate da P4VASP si trovano in tre colonne in cui la terza rappresenta la ponderazione.
- Caricare "vasprun.xml" digitando p4v vasprun.xml sul terminale nel sistema Ubuntu.
- Utilizzare xmgrace per modificare la banda proiettata.
- Digitare xmgrace 11-4.agr sul terminale per avviare xmgrace nel sistema Ubuntu. Selezionare il menu Proprietà traccia/asse per modificare l'etichetta e l'intervallo dell'asse.
- Selezionare il menu Stampa/Imposta aspetto per leggere il valore energetico in corrispondenza del numero di banda e del punto k specificati.
NOTA: il numero massimo della banda di valenza (VBM) e il minimo della banda di conduzione (CBM) di NR/NT possono essere letti dalla banda proiettata con contributi su NR/NT, rispettivamente. In base agli allineamenti delle bande, le eterostrutture possono essere classificate in tre tipi: tipo I (VBMNTNR ) , Tipo II (VBMNT NT NR - Calcolare l'offset della banda di valenza (VBO), l'offset della banda di conduzione (CBO) e lo spazio della banda dopo Kang et al.16.
- Selezionare il menu File/Stampa per esportare il grafico con il formato eps.
- Calcolare la densità di carica decomposta a banda per VBM e CBM.
- Preparare uno script PBS "job.pbs" e sette file di input: INCAR, POSCAR, POTCAR, KPOINTS, WAVECAR, CHGCAR e CHG. Specificare i numeri di banda per CBM e VBM dal tag IBAND in INCAR. Utilizzare il singolo punto k corrispondente per ogni bordo della banda.
NOTA: i file CHGCAR, CHG e WAVECAR preconvergenti provengono dai calcoli statici auto-coerenti nel passaggio 2.2. Esempi di file di input per questo passaggio sono forniti nei Materiali supplementari (ad eccezione di POTCAR). - Digitare qsub job.pbs sul terminale nel sistema Linux per inviare il processo.
- Utilizzare vmd per stampare il VBM e CBM nello spazio reale al termine del processo.
- Avviare una sessione vmd e caricare il file POSCAR come nel passaggio 1.1.4.
- Selezionare il menu File/Nuova molecola nella finestra principale vmd. Individuare il file PARCHG nella finestra Sfoglia. Caricare PARCHG digitando VASP_PARCHG.
- Selezionate i menu Disegna/Superficie solida e Mostra/Isosurface nella finestra Rappresentazioni grafiche. Modificare il valore isovalue in un valore appropriato (ad esempio, 0,02). Modificare il colore dell'isosuperficie tramite il menu Metodo colorazione.
NOTA: si tratta di un'analisi intuitiva per i tipi di banda rispetto a quella del passaggio 3.4. In genere, la struttura atomica è disposta lontano dal contorno, altrimenti la densità di carica visualizzata non viene visualizzata in modo continuo. Per informazioni dettagliate, vedere la figura supplementare 1.
- Preparare uno script PBS "job.pbs" e sette file di input: INCAR, POSCAR, POTCAR, KPOINTS, WAVECAR, CHGCAR e CHG. Specificare i numeri di banda per CBM e VBM dal tag IBAND in INCAR. Utilizzare il singolo punto k corrispondente per ogni bordo della banda.
4. Modulare le proprietà elettroniche del nanocomposito (NT incapsulato all'interno di NR) da campi esterni.
- Aggiungere un campo elettrico trasversale al nanocomposito17.
- Preparare uno script PBS "job.pbs" e quattro file di input: INCAR, POSCAR, POTCAR e KPOINTS.
- Definire l'intensità del campo elettrico con il tag "EFIELD" in unità di eV / .
- Impostare LDIPOL - T. Specificare IDIPOL con un valore esatto (1, 2 o 3).
NOTA: questi due tag vengono aggiunti per includere le correzioni del dipolo. Il campo elettrico verrà applicato lungo l'asse X, Y o z impostando il valore di IDIPOL su 1, 2 o 3. - Eseguire i calcoli statici auto-coerenti e la struttura delle bande seguendo le sezioni 2 e 3 senza ottimizzazione strutturale.
NOTA: Studi precedenti indicano che i campi elettrici su 5 V / s possono essere utilizzati per modificare lo spazio di banda di BN-NT e BN-NR senza deformare la struttura18,19.
- Aggiungere un ceppo di tensile longitudinale al nanocomposito.
- Modificare i parametri del reticolo lungo la direzione periodica per riflettere l'effetto di deformazione.
NOTA: Ad esempio, il parametro del reticolo ottimizzato del nanocomposito lungo l'asse z è 2,5045. Se lungo la direzione di , il parametro del reticolo in POSCAR viene applicato l'1% di deformazione uniassiale, impostare il parametro del reticolo in POSCAR su 2,5045 x 1,01 x 2,529545. - Rilassare la struttura modificata seguendo la sezione 1.
- Eseguire calcoli statici auto-coerenti e la struttura delle bande seguendo le sezioni 2 e 3.
- Modificare i parametri del reticolo lungo la direzione periodica per riflettere l'effetto di deformazione.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Gli ZOzag BN-NRs incapsulati all'interno della poltrona BN-NT (11,11) sono stati scelti come esempi rappresentativi per un'eterostruttura 1DW. I parametri del reticolo sono stati presi da Sahin et al.20. Per comodità, gli NR a zig-zag sono abbreviati,n, dove n rappresenta i dimeri III-V lungo la larghezza14. L'energia di incapsulamento EL del punto 2.3 è stata utilizzata come stima approssimativa per la stabilità energetica del nanocomposito. I valori EL di n.2,3e4 incapsulati all'interno di BN-NT (11,11) erano -0,033 eV /z, -0,068 eV / , e -0.131 eV / , rispettivamente10, come mostrato nella Figura 1. Anche se EL variava di un ordine di grandezza con dimensioni BN-NR, tutti e tre i nanocompositi presentavano strutture di banda di tipo II (dal passo 3.4) superiori ai casi di carbonio14, dove il tipo II è emerso solo per La NR con una sola dimensione appropriata inserito in NT14.
La struttura a banda del nanocomposito dal passo 3.2, BN-NT (11,11) -4, è mostrata nella Figura 2. VBM/CBM si trova rispettivamente in NT/NR (dal passaggio 3.5). L'allineamento a banda sfalsata è stato utile per la raccolta della luce. Il meccanismo principale di trasferimento di carica è il seguente: la foto genera elettroni e un foro in4 al punto X, mostrato nella Figura 3,e poi il foro dissocia da4 (kX) a NT (11,11) (kVBM, il punto k di VBM per questo nanocomposito), mostrato figura 4. Il VBO calcolato (dal passo 3.4.3) è 317 meV, più grande dell'energia termica a 300 K (KT 30 meV), e riduce efficacemente il tasso di ricombinazione dei vettori fotogenerati10.
Per migliorare la raccolta della luce attraverso un ampio spettro, sia i campi elettrici trasversali che i ceppi ditensile longitudinali vengono applicati a BN-NT (11,11) . L'evoluzione dei bordi delle bande rispetto al livello di vuoto dal punto 4 è illustrata nella Figura 5. Una sostanziale riduzione del divario fino a quasi 0,95 eV è osservata in questo nanocomposito da campi esterni. Ancora più importante, l'allineamento a banda sfalsata viene mantenuto10. Sulla base di questi risultati, un tale sistema 1D dovrebbe integrare la generazione di idrogeno fotocatalitico e lo stoccaggio sicuro di capsule21. Gli elettroni fotogenerati potrebbero essere raccolti da NR. Spinti dall'attrazione elettrostatica, i protoni penetrano attraverso il NT per generare una molecola di idrogeno. L'idrogeno prodotto è completamente isolato all'interno del nanotubi per evitare una reazione o un'esplosione inversa indesiderata.
Figura 1: nanoribbons di BN zigzagn. 2,3e4 incapsulati all'interno di un nanotubi BN (11,11). L'energia di incapsulamento (EL) è elencata in ogni struttura. Fare clic qui per visualizzare una versione più grande di questa figura.
Figura 2: Struttura a nastro del nanotubi BN (11,11) , nanonastro BNn. 4. I contributi del nanotubi e del nanonastro alle bande energetiche sono rappresentati rispettivamente nelle sfere rosse e blu. Gli inset a sinistra indicano le distribuzioni della densità di carica degli stati CBM e VBM (isovalore 0,02 e/3). Questa cifra è stata adattata da Gong et al.10 con il permesso della Royal Society of Chemistry. Fare clic qui per visualizzare una versione più grande di questa figura.
Figura 3: La foto genera elettroni e un foro nel nanonastro BNn. 4 nel punto X. Fare clic qui per visualizzare una versione più grande di questa figura.
Figura 4: Il foro dissocia dal nanoribbon di BN4 (kX)al nanotubi BN (11,11) (kVBM, il punto k di VBM per questo nanocomposito). Fare clic qui per visualizzare una versione più grande di questa figura.
Figura 5: Modulazione del bordo della banda del nanotubi BN (11,11) e del nanoribbon BNn. 4 da campi esterni. Evoluzione dei bordi delle bande rispetto al livello di vuoto sotto (A) un campo elettrico e (B) ceppo di tensile uniaxial. La direzione negativa del campo elettrico è indicata dall'atomo del bordo inferiore B all'atomo del bordo superiore N di4. Il potenziale di riduzione di H:/H2 e il potenziale di ossidazione di O2/H2O sono rispettivamente -4,44 eV e -5,67 eV a pH e 0. Il pH 7 sposta i potenziali di redox dell'acqua (da pH x 0,059 eV) a -4.027 eV e -5.257 eV, rispettivamente, mostrati come linee tratteggiate blu. Questa cifra è stata riprodotta da Gong et al.10 con il permesso della Royal Society of Chemistry. Fare clic qui per visualizzare una versione più grande di questa figura.
Figura supplementare 1: (A) Struttura atomica di un nanotubi BN (11,11) , nanonastro BNn. 4 disposto lontano dal confine e dalla corrispondente banda di conduzione minima (B). (C) Struttura atomica di un nanotubi BN (11,11) e di un nanoribbon BN4 allineata con un limite e il corrispondente minimo della banda di duzione (D). Fare clic qui per visualizzare una versione più grande di questa figura.
File di codifica supplementare: Fare clic qui per visualizzare questo file (fare clic con il pulsante destro del mouse per scaricare).
Subscription Required. Please recommend JoVE to your librarian.
Discussion
I calcoli per le proprietà elettroniche nelle sezioni 2, 3 e 4 sarebbero simili tra i vari materiali su nanoscala. Il modello atomico iniziale nel passaggio 1 deve essere attentamente progettato per estrarre informazioni significative. Ad esempio, il fattore per la selezione del modello potrebbe essere la dimensione o la chiralità dei materiali. Inoltre, il modello atomico iniziale nel passaggio 1.1 dovrebbe essere ragionevolmente preparato per il rilassamento della struttura a basso costo. Prendendo come esempio il nanocomposito nel protocollo, il NR deve essere incapsulato all'interno del NT in modo simmetrico. In caso contrario, la ricerca della struttura ottimizzata da parte di VASP richiederà molto tempo.
Per considerare l'effetto di un campo elettrico, un foglio di dipolo artificiale viene aggiunto al centro della parte del vuoto nella cella dell'unità periodica in VASP22. La regione del vuoto non dovrebbe essere troppo ampia e il campo elettrico dovrebbe essere sufficientemente debole da evitare l'emissione artificiale23.
Mentre l'effetto del ceppo può essere semplicemente realizzato cambiando il parametro del reticolo in POSCAR, nel nanocomposito la situazione sarebbe più complessa. Le risposte elastiche della NR e NT possono essere diverse l'una dall'altra, subendo la stessa forza. Ciò porterà ad una struttura sproporzionata. Ad esempio, quando la deformazione della tensione uniassiale viene applicata lungo la direzione periodica, il parametro del reticolo ottimizzato di NT e NR lungo questa direzione passa rispettivamente da un 1,8 a 2,0 e 2,2. Per la modellazione sono necessarie supercelle di grandi dimensioni: almeno 11 unità di celle di NT e 10 celle di unità di NR in questo caso (11 x 2,0 x 10 x 2,2 x 22).
Mentre le proprietà elettroniche a stato di terra dei materiali possono essere determinate da VASP abbastanza bene, per affrontare la durata dei fori e degli elettroni fotogenerati esistenti in uno stato eccitato, è meglio eseguire il calcolo della dinamica non adiabatica24. Questo è importante per progettare fotocatalizzatori con vettori a lunga vita4.
Il ruolo dell'approccio computazionale svolto da VASP consiste nella scoperta di nuovi materiali e nello screening di potenziali fotocatalizzatori per assistere gli sforzi sperimentali. L'allineamento delle bande a livello PBE nella scissione dell'acqua non è convincente come il lavoro quantitativo sperimentale. Sono necessari valori più precisi dei bordi delle bande rispetto ai potenziali redox, CBO e VBO. Sarebbe meglio usare il Heyd-Scuseria-Ernzerhof (HSE) ibrido funzionale25, ma è più tempo di PBE. Tuttavia, i risultati a livello di PBE potrebbero essere efficienti per suggerire strategie per il miglioramento dell'attività fotocatalitica.
Va ricordato che la progettazione computazionale di VASP consentirà anche la previsione di materiali per celle solari, materiali termoelettrici, materiali per batterie al litio, materiali di cattura di gas, ecc.2. I calcoli ad alta velocità effettiva sono stati combinati con le procedure di apprendimento automatico per una migliore previsione dei materiali e un costo computazionale inferiore26,27.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Gli autori non hanno nulla da rivelare.
Acknowledgments
Questo lavoro è stato sostenuto dalla China Postdoctoral Science Foundation (Grant No. 2017M612348), dalla Qingdao Postdoctoral Foundation (Grant n. 3002000-861805033070) e dal Young Talent Project presso l'Ocean University of China (Grant n. 3002000-861701013151). Gli autori ringraziano la signorina Ya Chong Li per aver preparato la narrazione.
Materials
| Name | Company | Catalog Number | Comments |
| Nanotube Modeler | Developed by Dr. Steffen Weber | NanotubeModeler1.8 | http://www.jcrystal.com/products/wincnt/NanotubeModeler.exe |
| P4VASP | Orest Dubay | p4vasp 0.3.30 | Open source, available at www.p4vasp.at |
| v2xsf | Developed by Dr. Jens Kunstmann | v2xsf | http://theory.chm.tu-dresden.de/~jk/software.html |
| VASP software | Computational Materials Physics, Dept. of Physics, University of Vienna | vasp.5.4.1 | https://www.vasp.at |
| VMD software | Theoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaign | vmd1.9.3 | https://www.ks.uiuc.edu/Research/vmd |
| xcrysden | Dept. of Physical and Organic Chemistry, Jozef Stefan Institute | XCrySDen1.5.60 | http://www.xcrysden.org/ |
| Xmakemol | Developed by M. P. Hodges | xmakemol5.16 | https://www.nongnu.org/xmakemol/XmakemolDownloads.html |
| Xmgrace software | Grace Development Team under the coordination of Evgeny Stambulchik | xmgrace5.1.25 | http://plasma-gate.weizmann.ac.il/Grace/ |
References
- Collins, C., et al. Accelerated discovery of two crystal structure types in a complex inorganic phase field. Nature. 546 (7657), 280-284 (2017).
- Jain, A., Shin, Y., Persson, K. A. Computational predictions of energy materials using density functional theory. Nature Reviews Materials. 1 (1), 15004 (2016).
- de Jong, M., et al. Charting the complete elastic properties of inorganic crystalline compounds. Scientific Data. 2, 150009 (2015).
- Fu, C. F., Wu, X. J., Yang, J. L. Material Design for Photocatalytic Water Splitting from a Theoretical Perspective. Advanced Materials. 30 (48), 1802106 (2018).
- Gu, T., Luo, W., Xiang, H. J. Prediction of two-dimensional materials by the global optimization approach. Wiley Interdisciplinary Reviews-Computational Molecular Science. 7 (2), e1295 (2017).
- Lejaeghere, K., et al. Reproducibility in density functional theory calculations of solids. Science. 351 (6280), aad3000 (2016).
- Kresse, G., Furthmüller, J. Efficient Iterative Schemes for ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Physical Review B. 54 (16), 11169-11186 (1996).
- Perdew, J. P., Burke, K., Ernzerhof, M. Generalized Gradient Approximation Made Simple. Physical Review Letters. 77 (18), 3865-3868 (1996).
- Ozcelik, V. O., Azadani, J. G., Yang, C., Koester, S. J., Low, T. Band alignment of two-dimensional semiconductors for designing heterostructures with momentum space matching. Physical Review B. 94 (3), 035125 (2016).
- Gong, M., et al. Robust staggered band alignment in one-dimensional van der Waals heterostructures: binary compound nanoribbons in nanotubes. Journal of Materials Chemistry C. 7 (13), 3829-3836 (2019).
- Grimme, S., Antony, J., Ehrlich, S., Krieg, H. A Consistent and Accurate ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. Journal of Chemical Physics. 132 (15), 154104 (2010).
- Zhang, L., Chen, Z. Q., Su, J., Li, J. F. Data mining new energy materials from structure databases. Renewable & Sustainable Energy Reviews. 107, 554-567 (2019).
- Zakutayev, A., et al. An open experimental database for exploring inorganic materials. Scientific Data. 5, 180053 (2018).
- Kou, L. Z., Tang, C., Frauenheim, T., Chen, C. F. Intrinsic Charge Separation and Tunable Electronic Band Gap of Armchair Graphene Nanoribbons Encapsulated in a Double-Walled Carbon Nanotube. Journal of Physical Chemistry Letters. 4 (8), 1328-1333 (2013).
- Kou, L. Z., Tang, C., Wehling, T., Frauenheim, T., Chen, C. F. Emergent properties and trends of a new class of carbon nanocomposites: graphene nanoribbons encapsulated in a carbon nanotube. Nanoscale. 5 (8), 3306-3314 (2013).
- Kang, J., Tongay, S., Zhou, J., Li, J. B., Wu, J. Q. Band offsets and heterostructures of two-dimensional semiconductors. Applied Physics Letters. 102 (1), 012111 (2013).
- Makov, G., Payne, M. C. Periodic boundary conditions in ab initio calculations. Physical Review B. 51 (7), 4014-4022 (1995).
- Chen, C., Lee, M., Clark, S. J. Band gap modification of singlewalled carbon nanotube and boron nitride nanotube under a transverse electric field. Nanotechnology. 15 (12), 1837 (2004).
- Zhang, Z. H., Guo, W. L. Energy-gap Modulation of BN Ribbons by Transverse Electric Fields: First-Principles Calculations. Physical Review B. 77 (7), 075403 (2008).
- Sahin, H., et al. Monolayer Honeycomb Structures of Group-IV Elements and III-V Binary Compounds: First-Principles Calculations. Physical Review B. 80 (15), 155453 (2009).
- Yang, L., et al. Combining Photocatalytic Hydrogen Generation and Capsule Storage in Graphene Based Sandwich Structures. Nature Communications. 8, 16049 (2017).
- Neugebauer, J., Scheffler, M. Adsorbate-substrate and adsorbate-adsorbate interactions of Na and K adlayers on Al(111). Physical Review B. 46 (24), 16067 (1992).
- He, W., Li, Z. Y., Yang, J. L., Hou, J. G. Electronic structures of organic molecule encapsulated BN nanotubes under transverse electric field. Journal of Chemical Physics. 124 (15), 154709 (2006).
- Zhang, R. Q., et al. Direct Z-Scheme Water Splitting Photocatalyst Based on Two-Dimensional Van Der Waals Heterostructures. Journal of Physical Chemistry Letters. 9 (18), 5419-5424 (2018).
- Paiera, J., Marsman, M., Hummer, K., Kresse, G. Screened hybrid density functionals applied to solids. Journal of Chemical Physics. 124 (15), 154709 (2006).
- Pyzer-Knapp, E. O., Suh, C., Gómez-Bombarelli, R., Aguilera-Iparraguirre, J., Aspuru-Guzik, A. What Is High-Throughput Virtual Screening? A Perspective from Organic Materials Discovery. Annual Review of Materials Research. 45, 195-216 (2015).
- Cerqueira, T. F. T., et al. Identification of Novel Cu, Ag, and Au Ternary Oxides from Global Structural Prediction. Chemistry of Materials. 27 (13), 4562-4573 (2015).
