Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
JoVE Journal Chemistry
Click here for the English version

Chemistry

Alignement de bande de type II de sonde dans les hétérostructures unidimensionnelles de Van Der Waals utilisant des calculs de premiers principes

Published: October 12, 2019 doi: 10.3791/60180
Automatic Translation
1, 1, 1, 2
1College of Information Science and Engineering, Ocean University of China, 2Department of Physics, Southern University of Science and Technology

Summary

Les calculs effectués par le paquet de simulation de l'initio de Vienne Ab peuvent être utilisés pour identifier les propriétés électroniques intrinsèques des matériaux à l'échelle nanométrique et prédire les photocatalyseurs potentiels de division de l'eau.

Abstract

Les outils informatiques basés sur la théorie de la densité-fonctionnelle (DFT) permettent l'exploration des composés nanométriques qualitativement nouveaux et expérimentalement réalisables pour une application ciblée. Les simulations théoriques permettent de bien comprendre les propriétés électroniques intrinsèques des matériaux fonctionnels. Le but de ce protocole est de rechercher des candidats photocatalyseurs par dissection computationnelle. Les applications photocatalytiques nécessitent des lacunes de bande appropriées, des positions appropriées de bord de bande par rapport aux potentiels redox. Les fonctions hybrides peuvent fournir des valeurs précises de ces propriétés, mais elles sont coûteuses sur le plan du calcul, alors que les résultats au niveau fonctionnel de Perdew-Burke-Ernzerhof (PBE) pourraient être efficaces pour suggérer des stratégies d'ingénierie de la structure de la bande par l'intermédiaire le champ électrique et la contrainte tendue visant à augmenter la performance photocatalytique. Pour illustrer cela, dans le présent manuscrit, l'outil de simulation basé sur le DFT VASP est utilisé pour étudier l'alignement des bandes de nanocomposites dans des combinaisons de nanotubes et de nanoribbons dans l'état du sol. Pour répondre à la durée de vie des trous photogénérés et des électrons dans l'état excité, des calculs de dynamique nonadiabatic sont nécessaires.

Introduction

La demande mondiale d'énergie propre et durable a stimulé la recherche de matériaux prometteurs afin de réduire la dépendance à l'égard des ressources pétrolières limitées. Les simulations sont plus efficaces et plus économiques que les expériences visant à accélérer la recherche de nouveaux matériaux fonctionnels1. La conception des matériaux d'un point de vue théorique2,3,4est maintenant de plus en plus populaire en raison des progrès rapides dans les ressources de calcul et les développements théoriques, ce qui rend les simulations de calcul plus fiables 5 . Les calculs de la théorie fonctionnelle de la densité (DFT) mis en œuvre dans de nombreux codes sont de plus en plus robustes et donnent des résultats reproductibles6.

Le Vienna Ab initio Simulation Package (VASP)7 présente l'un des codes DFT les plus prometteurs pour prédire les propriétés moléculaires et cristallines et plus de 40 000 études faisant usage de ce code ont été publiées. La plupart des travaux sont exécutés au niveau fonctionnel8de Perdew-Burke-Ernzerhof (PBE), qui sous-estime la taille des écarts de bande, mais qui capture les tendances essentielles en matière d'alignement des bandes et de décalages de bande3. Ce protocole vise à décrire les détails de l'étude des profils de bord de bande et des bandgaps de matériaux à l'échelle nanométrique pour l'énergie propre et renouvelable à l'aide de cet outil de calcul. D'autres exemples utilisant VASP sont disponibles à https://www.vasp.at.

Ce rapport présente le criblage computationnel des hétérostructures unidimensionnelles (1D) vdW avec des alignements de bande de type II9 pour une application prometteuse dans le fractionnement photocatalytique d'eau4. Plus précisément, les nanoribbons (NR) encapsulés à l'intérieur des nanotubes (NT) sont examinés à titre d'exemple10. Pour traiter les interactions non covalentes, les corrections vdW utilisant la méthode DFT-D3 sont incluses11. Les calculs DFT dans les étapes 1.2, 2.2, 3.2, 3.5.2, et la section 4 par VASP sont exécutés à l'aide d'un script de système de lots portables (PBS) par les ordinateurs de recherche haute performance dans le système CenTOS. Un exemple de script PBS est affiché dans les Matériaux Supplémentaires. Le post-traitement des données par le logiciel P4VASP à l'étape 3.3 et la parcelle de figure par le logiciel xmgrace à l'étape 3.4 sont effectués sur un ordinateur local (ordinateur portable ou bureau) dans le système Ubuntu.

Subscription Required. Please recommend JoVE to your librarian.

Protocol

1. Optimiser la structure atomique.

  1. Préparer quatre fichiers d'entrée pour le calcul de la relaxation de la structure par VASP: INCAR, POSCAR, POTCAR et KPOINTS.
    REMARQUE : Il y a des paramètres spécifiés dans le fichier INCAR qui définissent le calcul. La ligne "EDIFFG - 0.02" dans le fichier INCAR indique que tous les atomes sont détendus jusqu'à ce que la force sur chaque atome soit 'lt;0.02 eV/' . Le fichier POSCAR contient les informations de géométrie atomique. Les paramètres de treillis initiaux dans le fichier POSCAR peuvent être choisis à partir de référencesthéoriques 3 ou expérimentales12,13. Le fichier KPOINTS définit le maillage k point et POTCAR est le fichier pseudopotentiel. L'ordre des types d'atomes dans POSCAR devrait être le même que celui de POTCAR. Des exemples de fichiers d'entrée pour la relaxation de la structure sont affichés dans les matériaux supplémentaires (à l'exception du fichier pseudopotentiel, qui a besoin d'une licence de VASP).
    1. Générer la structure initiale des nanoribbons de nitride de bore (BN) (NR) pour "POSCAR".
      1. Téléchargez le fichier POSCAR pour l'unité en vrac BN à partir de https://materialsproject.org.
      2. Utilisez v2xsf pour convertir le fichier POSCAR en un fichier en format xsf qui peut être lu par xcrysden. Type v2xsf POSCAR sur le terminal dans le système Ubuntu pour obtenir "POSCAR.xsf.gz". Type gunzip POSCAR.xsf.gz et sortir le fichier POSCAR.xsf.
      3. Utilisez xcrysden pour construire la supercellule BN.
        1. Type xcrysden --xsf POSCAR.xsf sur le terminal du système Ubuntu. Sélectionnez le menu Modifier/Nombre d'unités dessinées et étendre la cellule dans les directions X et Y.
        2. Sélectionnez le menu Fichier/Enregistrer la structure XSF pour exporter la structure supercellulaire, nommée «supercellule».
          REMARQUE : Le nom de la structure est une définition arbitraire.
      4. Utilisez du xmakemol pour ouvrir la supercellule. Type xmakemol-f supercellule sur le terminal dans le système Ubuntu. Sélectionnez le menu Edit/Visible. Cliquez sur Toggle pour supprimer les atomes à l'intérieur de la région et couper le NR à la largeur désirée et la chiralité.
    2. Générer la structure initiale du nanotube BN (NT) pour POSCAR. Télécharger "NanotubeModeler" de http://www.jcrystal.com/products. Ouvrez NanotubeModeler.exe dans le système Windows. Sélectionnez le menu Sélectionnez le type/B-N et spécifiez la chiralité. Sélectionnez le menu Fichier/Enregistrer la table XYZ pour exporter la structure.
    3. Générez la structure initiale du nanocomposite en encapsulant le NR (à partir de l'étape 1.1.1) à l'intérieur du NT (à partir de l'étape 1.1.2).
      REMARQUE: L'encapsulation peut être terminée en ajustant les coordonnées cartésiennes du NR et le NT10,14,15.
    4. Utilisez le logiciel vmd pour vérifier la structure atomique avant de soumettre le travail de calcul.
      1. Type vmd sur le terminal dans le système Ubuntu. Dans la fenêtre principale vmd ouverte, sélectionnez le fichier de menu/nouvelle molécule et trouvez le fichier POSCAR à travers la fenêtre Parcourir. Chargez POSCAR en tapant VASP-POSCAR.
      2. Afficher la structure dans différents styles dans la fenêtre Représentations graphiques/méthode de dessin.
        REMARQUE : Par exemple, une fois que le CPK est choisi, chaque atome (lien) est représenté par une sphère (bâton). Le guide d'installation et le tutoriel complet de vmd sont disponibles à http://www.ks.uiuc.edu/Research/vmd.
  2. Tapez qsub job.pbs sur le terminal du système Linux pour soumettre le travail au cluster informatique.
    REMARQUE: "job.pbs" représente le nom du script PBS. Le nom du script PBS est une définition arbitraire. Les quatre fichiers d'entrée ainsi que le script PBS doivent être à l'annuaire de travail. La commande qsub job.pbs sera utilisée dans les étapes 2.2, 3.2, 3.5.2, et la section 4. Un exemple de script PBS peut être trouvé dans le fichier de codage supplémentaire. Une fois le travail soumis terminé, si «l'exactitude requise - arrêt dela minimisation de l'énergie structurelle » apparaît à la fin du journal de sortie, le résultat convergent est obtenu. Le fichier CONTCAR qui en résultera sera utilisé comme fichier d'entrée POSCAR dans les étapes 2.1, 3.1, 3.5.1, 3.5.3.1, 4.1.1, 4.1.4, et section 4.2.

2. Calculer l'énergie d'encapsulation.

  1. Tapez mkdir nanocomposite isolé-nanoribbon isolé-nanotube pour créer trois dossiers pour le nanocomposite, le NR, et le NT sur un terminal dans un système Linux. Préparer un script PBS "job.pbs" et quatre fichiers d'entrée INCAR, POSCAR, POTCAR et KPOINTS pour le calcul de l'énergie dans chaque dossier.
    REMARQUE: Le fichier d'entrée POSCAR est le fichier nommé CONTCAR avec la structure détendue de l'étape 1. Des exemples de fichiers d'entrée sont donnés dans des documents supplémentaires (à l'exception de POTCAR).
  2. Allez à chaque dossier et tapeqjob.pbs sur le terminal dans le système Linux.
    REMARQUE : Les trois tâches soumises effectueront les calculs statiques d'énergie auto-consistantes pour le nanocomposite, le NR isolé et le NT isolé, respectivement.
  3. Extraire l'énergie totale du fichier OUTCAR pour chaque système après avoir terminé les calculs statiques auto-cohérents. Type grep "énergie libre TOTEN" ./nanocomposite/OUTCAR - queue -n 1, grep "énergie libre TOTEN" ./isolated-nanoribbon/OUTCAR - tail -n 1, et grep "free energy TOTEN" ./isolated-nanotube/OUTCAR - tail -n 1. Définissez les trois valeurs affichées comme ENT-NR, ENR, et ENT, respectivement. Calculez l'énergie d'encapsulation par angstrom : EL (ENT-NR - ENT -ENR)/L14,15.
    REMARQUE : La direction périodique dans chaque système est le long de l'axe Z et L est la constante de treillis de la cellule unitaire le long de l'axe Z. Des calculs d'essai de la dépendance d'énergie sur l'énergie de coupure d'onde de plan et le maillage de point de k sont nécessaires. L'énergie d'encapsulation peut être utilisée comme estimation de la stabilité énergétique du nanocomposite.

3. Extraire les propriétés électroniques de la structure de la bande.

  1. Préparer un script PBS "job.pbs" et six fichiers d'entrée: INCAR, POSCAR, POTCAR, KPOINTS, CHGCAR et CHG pour le calcul de la bande. Définir ICHARG 11 en INCAR.
    REMARQUE : Les fichiers CHGCAR et CHG préconvergés proviennent des calculs statiques auto-cohérents à l'étape 2.2. L'analyse de bande est au niveau PBE. L'échantillonnage de k point dans le fichier KPOINTS est en mode en ligne. Des exemples de fichiers d'entrée pour cette étape peuvent être trouvés dans les matériaux supplémentaires (sauf POTCAR).
  2. Tapez qsub job.pbs sur le terminal du système Linux pour soumettre le travail.
  3. Utilisez P4VASP pour générer la bande projetée.
    1. Chargez "vasprun.xml" en tapant p4v vasprun.xml sur le terminal du système Ubuntu.
      REMARQUE: "p4v" est utilisé pour démarrer P4VASP. Le fichier "vasprun.xml" doit être à l'annuaire de travail.
    2. Sélectionnez le menu Contrôle électronique/Local DOS-bands, puis Sélectionnez/Bandes.
      1. Spécifiez les numéros atomiques de NT dans la section Sélection atome. Obtenez le nombre atomique en pointant vers les atomes correspondants en utilisant vmd comme mentionné à l'étape 1.1.4. Spécifiez la couleur, le type et la taille du symbole pour la structure de bande projetée à travers le menu Symbole et La taille du symbole. Appuyez sur le menu Ajouter une nouvelle ligne.
        REMARQUE : Le graphique montrera la structure de la bande avec les contributions du NT.
      2. Répétez la même procédure après l'étape 3.3.2.1 pour obtenir la bande projetée avec les contributions de la NR.
    3. Sélectionnez le menu Graphique/Exportation. Exportez le graphique dans un fichier avec un format agr (par exemple, comme « 11-4.agr »).
      REMARQUE : Les données de sortie des bandes projetées par P4VASP sont en trois colonnes où la troisième représente la pondération.
  4. Utilisez xmgrace pour modifier la bande projetée.
    1. Type xmgrace 11-4.agr sur le terminal pour démarrer xmgrace dans le système Ubuntu. Sélectionnez le menu Propriétés Plot/Axis pour modifier l'étiquette et la plage de l'axe.
    2. Sélectionnez l'apparence De la parcelle de menu/set pour lire la valeur énergétique au numéro de bande et au point k spécifiés.
      REMARQUE : Le maximum de bande de valence (VBM) et le minimum de bande de conduction (CBM) de NR/NT peuvent être lus de la bande projetée avec des contributions sur NR/NT, respectivement. Selon les alignements de la bande, les hétérostructures peuvent être classées en trois types : type I (VBMNT 'lt;VBMNR 'lt;CBMNR 'lt;CBMNT ou VBMNR 'lt;VBMNT 'lt;CBMNT 'lt;CBMNR), type II (VBMNT 'lt;VBMNR 'lt;CBMNT 'lt;CBMNR ou VBMNR 'lt;VBMNT 'lt;CBMNR 'lt;CBMNT), ou type III (VBMNT 'lt;VBMNT 'lt;CBMNR Lt;CBMNR ou VBMNR 'lt;VBMNR 'lt;CBMNT 'lt;CBMNT)9.
    3. Calculer le décalage de la bande de valence (VBO), la bande de conduction offset (CBO), et l'écart de bande suivant Kang et al.16.
    4. Sélectionnez le menu Fichier/Impression pour exporter le graphique avec le format eps.
  5. Calculez la densité de charge décomposée de la bande pour le VBM et le CBM.
    1. Préparer un script PBS "job.pbs" et sept fichiers d'entrée: INCAR, POSCAR, POTCAR, KPOINTS, WAVECAR, CHGCAR et CHG. Spécifiez les numéros de bande pour le CBM et VBM par l'étiquette IBAND dans INCAR. Utilisez le seul point k correspondant pour chaque bord de bande.
      REMARQUE : Les fichiers CHGCAR, CHG et WAVECAR préconvergés proviennent des calculs statiques auto-cohérents de l'étape 2.2. Des exemples de fichiers d'entrée pour cette étape sont donnés dans les documents supplémentaires (sauf POTCAR).
    2. Tapez qsub job.pbs sur le terminal du système Linux pour soumettre le travail.
    3. Utilisez vmd pour tracer le VBM et CBM dans l'espace réel après le travail est terminé.
      1. Démarrez une session vmd et chargez le fichier POSCAR comme à l'étape 1.1.4.
      2. Sélectionnez le fichier de menu/nouvelle molécule dans la fenêtre principale vmd. Trouvez le fichier PARCHG via la fenêtre Parcourir. Chargez PARCHG en tapant VASP-PARCHG.
      3. Sélectionnez les menus Dessiner/Surface solide et Afficher/Isosurface dans la fenêtre Représentations graphiques. Modifier l'isovalue à une valeur appropriée (par exemple, 0,02). Changer la couleur de l'isosurface à travers le menu Méthode de coloriage.
        REMARQUE : Il s'agit d'une analyse intuitive des types de bandes à l'étape 3.4. En général, la structure atomique est disposée loin de la limite, sinon la densité de charge visualisée n'est pas montrée d'une manière continue. S'il vous plaît voir La figure supplémentaire 1 pour plus de détails.

4. Moduler les propriétés électroniques du nanocomposite (NT encapsulé à l'intérieur de NR) par des champs externes.

  1. Ajouter un champ électrique transversal au nanocomposite17.
    1. Préparer un script PBS "job.pbs" et quatre fichiers d'entrée: INCAR, POSCAR, POTCAR, et KPOINTS.
    2. Définissez la force du champ électrique par l'étiquette «EFIELD» en unités d'eV/MD.
    3. Définir LDIPOL et T. Spécifiez IDIPOL avec une valeur exacte (1, 2 ou 3).
      REMARQUE : Ces deux balises sont ajoutées pour inclure des corrections dipole. Le champ électrique sera appliqué le long de l'axe X, Y ou Z en fixant la valeur d'IDIPOL à 1, 2 ou 3.
    4. Effectuer les calculs statiques auto-cohérents et les calculs de structure de bande suivant les sections 2 et 3 sans optimisation structurelle.
      REMARQUE : Des études antérieures indiquent que les champs électriques de plus de 5 V/MD peuvent être utilisés pour modifier l'écart de bande de BN-NT et de BN-NR sans déformer la structure18,19.
  2. Ajouter une souche tendinudinale longitudinale au nanocomposite.
    1. Modifier les paramètres de treillis le long de la direction périodique pour refléter l'effet de contrainte.
      REMARQUE : Par exemple, le paramètre de treillis optimisé du nanocomposite le long de l'axe Z est de 2,5045 euros. Si une souche de tension uniaxiale de 1 % est appliquée le long de la direction Z, changez le paramètre de treillis dans POSCAR à 2,5045 x 1,01 à 2,529545.
    2. Détendez la structure modifiée après l'article 1.
    3. Effectuer des calculs statiques auto-cohérents et des calculs de structure de bande suivant les sections 2 et 3.

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

Zigzag BN-NRs encapsulé à l'intérieur du fauteuil BN-NTs (11,11) ont été choisis comme exemples représentatifs d'une hétérostructure 1D vdW. Les paramètres de treillis ont été pris de Sahin et al.20. Pour plus de commodité, les NR en zigzag sont abrégés Zn, où n représente les dimers III-V le long de la largeur14. L'énergie d'encapsulation EL de l'étape 2.3 a été utilisée comme estimation approximative de la stabilité énergétique du nanocomposite. Les valeurs EL de Z2, Z3, et Z4 encapsulées à l'intérieur de BN-NT (11,11) étaient de -0,033 eV/MD, -0,068 eV/MD, et -0,131 eV/MD, respectivement10, comme le montre la figure 1. Bien que l'EL ait varié par ordre de grandeur avec la taille bn-NR, les trois nanocomposites ont présenté des structures de bande de type II (à partir de l'étape 3.4) supérieures aux cas tout carbone14, où le type II n'a émergé que pour NR avec une seule taille appropriée inséré dans NT14.

La structure de bande du nanocomposite de l'étape 3.2, BN-NT (11,11) - Z4, est montrée dans la figure 2. VBM/CBM se situe à NT/NR (à partir de l'étape 3.5), respectivement. L'alignement décalé de la bande a été bénéfique pour la récolte de lumière. Le principal mécanisme de transfert de charge est le suivant: la photo génère des électrons et un trou dans Z4 au point X, montré dans la figure 3, puis le trou se dissocie de Z4 (kX) à NT (11,11) (kVBM, le point k de VBM pour ce nanocomposite), montré dans la figure 4. Le VBO calculé (à partir de l'étape 3.4.3) est 317 meV, plus grand que l'énergie thermique à 300 K (KT 30 meV), et diminue effectivement le taux de recombinaison des porteurs photogénérés10.

Pour améliorer la récolte de lumière à travers un large spectre, les champs électriques transversaux et les souches tendineuses tendinitinales sont appliqués à BN-NT (11,11) et Z4. L'évolution des bords de bande par rapport au niveau du vide à partir de l'étape 4 est indiquée dans la figure 5. Une réduction substantielle de l'écart jusqu'à près de 0,95 eV est observée dans ce nanocomposite par des champs externes. Plus important encore, l'alignement décalé de bande est préservé10. Sur la base de ces résultats, un tel système 1D devrait intégrer la production d'hydrogène photocatalytique et le stockage de capsules sûrs21. Les électrons photogénérés pourraient être collectés par NR. Poussés par l'attraction électrostatique, les protons pénètrent à travers le NT pour générer une molécule d'hydrogène. L'hydrogène produit est complètement isolé dans le nanotube pour éviter une réaction ou une explosion inverse indésirable.

Figure 1
Figure 1 : Zigzag BN nanoribbons Z2, Z3, et Z4 encapsulés à l'intérieur d'un nanotube BN (11,11). L'énergie d'encapsulation (EL) est répertoriée sous chaque structure. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Figure 2
Figure 2 : Structure de bande du nanotube de BN (11,11) - BN nanoribbon Z4. Les contributions du nanotube et du nanoribbon aux bandes d'énergie sont représentées dans les sphères rouge et bleue, respectivement. Les insets de gauche montrent les distributions de densité de charge du CBM et des états VBM (isovalue 0,02e/3). Ce chiffre a été adapté de Gong et coll.10 avec la permission de la Royal Society of Chemistry. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Figure 3
Figure 3 : La photo génère des électrons et un trou dans le nanoribbon Z4 du BN au point X. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Figure 4
Figure 4 : Le trou se dissocie du nanoribbon BN Z4 (kX) au nanotube BN (11,11) (kVBM, le point k de VBM pour ce nanocomposite). Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Figure 5
Figure 5 : Modulation des bords de bande du nanotube BN (11,11) et du nanoribbon Z4 du BN par les champs externes. Evolution des bords de bande par rapport au niveau du vide sous (A) un champ électrique et (B) souche de tension uniaxiale. La direction négative du champ électrique est dénotée de l'atome de bord inférieur B à l'atome de bord supérieur N de Z4. Le potentiel de réduction de H/H2 et le potentiel d'oxydation de O2/H2O sont de -4,44 eV et de -5,67 eV au pH et 0, respectivement. Le pH no 7 déplace les potentiels redox de l'eau (par pH x 0,059 eV) à -4,027 eV et -5,257 eV, respectivement, indiqués sous forme de lignes pointillées bleues. Ce chiffre a été reproduit à partir de Gong et coll.10 avec la permission de la Royal Society of Chemistry. Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Supplemental Figure 1
Figure supplémentaire 1 : (A) Structure atomique d'un nanotube BN (11,11) - BN nanoribbon Z4 arrangé loin de la frontière et de sa bande de conduction correspondante minimum (B). (C) Structure atomique d'un nanotube BN (11,11) et BN nanoribbon Z4 aligné avec une limite et sa bande de conduction correspondante minimum (D). Veuillez cliquer ici pour voir une version plus grande de ce chiffre.

Fichier de codage supplémentaire: S'il vous plaît cliquez ici pour voir ce fichier (Clic droit pour télécharger).

Subscription Required. Please recommend JoVE to your librarian.

Discussion

Les calculs des propriétés électroniques dans les sections 2, 3 et 4 seraient similaires parmi divers matériaux à l'échelle nanométrique. Le modèle atomique initial à l'étape 1 doit être soigneusement conçu pour extraire des informations significatives. Par exemple, le facteur de sélection du modèle pourrait être la taille ou la chiralité des matériaux. En outre, le modèle atomique initial à l'étape 1.1 devrait être raisonnablement préparé pour la relaxation de la structure à faible coût. Prenant le nanocomposite dans le protocole comme exemple, le NR devrait être encapsulé à l'intérieur du NT d'une manière symétrique. Dans le cas contraire, il sera long de rechercher la structure optimisée par VASP.

Pour tenir compte de l'effet d'un champ électrique, une feuille de dipole artificielle est ajoutée au milieu de la partie vide dans la cellule unitaire périodique dans VASP22. La région de vide ne devrait pas être trop large et le champ électrique devrait être assez faible pour éviter l'émission artificielle de champ23.

Alors que l'effet de la souche peut être simplement réalisé en changeant le paramètre de treillis dans POSCAR, dans le nanocomposite la situation serait plus complexe. Les réponses élastiques du NR et du NT peuvent être différentes les unes des autres, subissant la même force. Cela conduira à une structure disproportionnée. Par exemple, lorsque la souche de tension uniaxiale est appliquée le long de la direction périodique, le paramètre de treillis optimisé du NT et du NR le long de cette direction passe d'un initial de 1,8 à 2,0 à 2,2 euros, respectivement. De grandes supercellules sont nécessaires à la modélisation : au moins 11 cellules unitaires de NT et 10 cellules unitaires de NR dans ce cas (11 x 2,0 ' 10 x 2,2 '22 ').

Tandis que les propriétés électroniques d'état de sol des matériaux peuvent être déterminées par VASP tout à fait bien, pour adresser la durée de vie des trous photogénérés et des électrons existants dans un état excité, il est meilleur d'exécuter le calcul nonadiabatic de dynamique24. Ceci est important pour concevoir des photocatalyseurs avec des porteurs à long terme4.

Le rôle de l'approche computationnelle effectuée par VASP joue dans la découverte de nouveaux matériaux et la sélection de photocatalyseurs potentiels pour aider les efforts expérimentaux. L'alignement de la bande au niveau pBE dans le fractionnement de l'eau n'est pas aussi convaincant que le travail expérimental quantitatif. Des valeurs plus précises des bords de bande par rapport aux potentiels redox, CBO, et VBO sont nécessaires. Il serait préférable d'utiliser le Heyd-Scuseria-Ernzerhof (HSE) hybride fonctionnel25, mais il est plus long que PBE. Néanmoins, les résultats au niveau PBE pourraient être efficaces pour suggérer des stratégies pour l'amélioration de l'activité photocatalytique.

Il convient de mentionner que la conception informatique par VASP permettra également la prédiction des matériaux de cellules solaires, matériaux thermoélectriques, matériaux de batterie au lithium, matériaux de capture de gaz, etc2. Des calculs à haut débit ont été combinés avec les procédures d'apprentissage automatique pour une meilleure prévision des matériaux et un coût de calcul inférieur26,27.

Subscription Required. Please recommend JoVE to your librarian.

Disclosures

Les auteurs n'ont rien à révéler.

Acknowledgments

Ce travail a été soutenu par la China Postdoctoral Science Foundation (Grant No. 2017M612348), la Qingdao Postdoctoral Foundation (Grant No. 3002000-861805033070) et par le Young Talent Project de l'Ocean University of China (Grant No. 3002000-861701013151). Les auteurs remercient Mlle Ya Chong Li d'avoir préparé la narration.

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/

DOWNLOAD MATERIALS LIST

References

  1. Collins, C., et al. Accelerated discovery of two crystal structure types in a complex inorganic phase field. Nature. 546 (7657), 280-284 (2017).
  2. Jain, A., Shin, Y., Persson, K. A. Computational predictions of energy materials using density functional theory. Nature Reviews Materials. 1 (1), 15004 (2016).
  3. de Jong, M., et al. Charting the complete elastic properties of inorganic crystalline compounds. Scientific Data. 2, 150009 (2015).
  4. 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).
  5. 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).
  6. Lejaeghere, K., et al. Reproducibility in density functional theory calculations of solids. Science. 351 (6280), aad3000 (2016).
  7. 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).
  8. Perdew, J. P., Burke, K., Ernzerhof, M. Generalized Gradient Approximation Made Simple. Physical Review Letters. 77 (18), 3865-3868 (1996).
  9. 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).
  10. 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).
  11. 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).
  12. 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).
  13. Zakutayev, A., et al. An open experimental database for exploring inorganic materials. Scientific Data. 5, 180053 (2018).
  14. 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).
  15. 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).
  16. 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).
  17. Makov, G., Payne, M. C. Periodic boundary conditions in ab initio calculations. Physical Review B. 51 (7), 4014-4022 (1995).
  18. 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).
  19. 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).
  20. 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).
  21. Yang, L., et al. Combining Photocatalytic Hydrogen Generation and Capsule Storage in Graphene Based Sandwich Structures. Nature Communications. 8, 16049 (2017).
  22. 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).
  23. 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).
  24. 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).
  25. Paiera, J., Marsman, M., Hummer, K., Kresse, G. Screened hybrid density functionals applied to solids. Journal of Chemical Physics. 124 (15), 154709 (2006).
  26. 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).
  27. 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).

Tags

Chimie Numéro 152 hétérostructures unidimensionnelles de van der Waals alignement de bande de type II nanoribbon nanotube calculs de premier principe bande de valence maximale bande de conduction minimum
Alignement de bande de type II de sonde dans les hétérostructures unidimensionnelles de Van Der Waals utilisant des calculs de premiers principes
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Hu, H., Lu, D., Dou, K. P., Shi, X.… More

Hu, H., Lu, D., Dou, K. P., Shi, X. Q. Probe Type II Band Alignment in One-Dimensional Van Der Waals Heterostructures Using First-Principles Calculations. J. Vis. Exp. (152), e60180, doi:10.3791/60180 (2019).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter