A subscription to JoVE is required to view this content. Sign in or start your free trial.

Method Article

Probe Type II Band Alignment in One-Dimensional Van Der Waals Heterostructures Using First-Principles Calculations

7.3K views

⸱

DOI:

10.3791/60180

⸱

October 12th, 2019

 ,  ,  , 
1College of Information Science and Engineering, Ocean University of China, 2Department of Physics, Southern University of Science and Technology

In This Article

Summary

Calculations performed by the Vienna Ab initio Simulation Package can be used to identify the intrinsic electronic properties of nanoscale materials and predict the potential water-splitting photocatalysts.

Abstract

Computational tools based on density-functional theory (DFT) enable the exploration of the qualitatively new, experimentally attainable nanoscale compounds for a targeted application. Theoretical simulations provide a profound understanding of the intrinsic electronic properties of functional materials. The goal of this protocol is to search for photocatalyst candidates by computational dissection. Photocatalytic applications require suitable band gaps, appropriate band edge positions relative to the redox potentials. Hybrid functionals can provide accurate values of these properties but are computationally expensive, whereas the results at the Perdew-Burke-Ernzerhof (PBE) functional level could be effective for suggesting strategies for band structure engineering via electric field and tensile strain aiming to enhance the photocatalytic performance. To illustrate this, in the present manuscript, the DFT based simulation tool VASP is used to investigate the band alignment of nanocomposites in combinations of nanotubes and nanoribbons in the ground state. To address the lifetime of photogenerated holes and electrons in the excited state, nonadiabatic dynamics calculations are needed.

Introduction

The worldwide demand for clean and sustainable energy has spurred research for promising materials to reduce dependence on finite petroleum resources. Simulations are more efficient and economical than experiments in accelerating the search for new functional materials1. Material design from a theoretical perspective2,3,4 is now more and more popular due to rapid advances in computational resources and theory developments, making computational simulations more reliable5. The density functional theory (DFT) calculations implemented in many codes are becoming more robust and yield reproducible results6.

The Vienna Ab initio Simulation Package (VASP)7 presents one of the most promising DFT codes for predicting molecular and crystalline properties and more than 40,000 studies making use of this code have been published. Most work is performed at the Perdew-Burke-Ernzerhof (PBE) functional level8, which underestimates the band gap sizes, but captures the essential trends in band alignment and band offsets3. This protocol aims to outline the details of investigating the band edge profiles and bandgaps of nanoscale materials for clean and renewable energy using this computational tool. More examples using VASP are available at https://www.vasp.at.

This report presents the computational screening of one-dimensional (1D) vdW heterostructures with type II band alignments9 for a promising application in photocatalytic water splitting4. Specifically, nanoribbons (NRs) encapsulated inside nanotubes (NTs) are examined as an example10. To address noncovalent interactions, vdW corrections using the DFT-D3 method are included11. The DFT calculations in steps 1.2, 2.2, 3.2, 3.5.2, and section 4 by VASP are performed using a Portable Batch System (PBS) script by the high-performance research computers in the CenTOS system. An example of a PBS script is shown in the Supplementary Materials. The data postprocessing by the P4VASP software in step 3.3 and the figure plot by the xmgrace software in step 3.4 are carried on a local computer (laptop or desktop) in the Ubuntu system.

Access restricted. Please log in or start a trial to view this content.

Protocol

1. Optimize the atomic structure.

  1. Prepare four input files for structure relaxation calculation by VASP: INCAR, POSCAR, POTCAR, and KPOINTS.
    NOTE: There are specified parameters in the INCAR file that define the calculation. The line "EDIFFG = 0.02" in the INCAR file indicates that all atoms are relaxed until the force on each atom is <0.02 eV/Å. The POSCAR file contains the atomic geometry information. The initial lattice parameters in the POSCAR file can be chosen from theoretical3 or experimental references12,13. The KPOINTS file defines the k point mesh and POTCAR is the pseudopotential file. The order of atom types in POSCAR should be the same as that in POTCAR. Examples of input files for structure relaxation are shown in the Supplementary Materials (except the pseudopotential file, which needs a license from VASP).
    1. Generate the initial structure of boron nitride (BN) nanoribbons (NR) for "POSCAR".
      1. Download the POSCAR file for the BN bulk unit from https://materialsproject.org.
      2. Use v2xsf to convert the POSCAR file to a file in xsf format that can be read by xcrysden. Type v2xsf POSCAR on the terminal in the Ubuntu system to get "POSCAR.xsf.gz". Type gunzip POSCAR.xsf.gz and output the POSCAR.xsf file.
      3. Use xcrysden to build the BN supercell.
        1. Type xcrysden --xsf POSCAR.xsf on the terminal in the Ubuntu system. Select the menu Modify/Number of Units Drawn and extend the cell in the X and Y directions.
        2. Select the menu File/Save XSF Structure to export the supercell structure, named "supercell".
          NOTE: The name of the structure is an arbitrary definition.
      4. Use xmakemol to open the supercell. Type xmakemol -f supercell on the terminal in the Ubuntu system. Select the menu Edit/Visible. Click Toggle to delete the atoms inside the region and cut the NR to the desired width and chirality.
    2. Generate the initial structure of the BN nanotube (NT) for POSCAR. Download "NanotubeModeler" from http://www.jcrystal.com/products. Open NanotubeModeler.exe in the Windows system. Select the menu Select type/B-N and specify the chirality. Select the menu File/Save XYZ table to export the structure.
    3. Generate the initial structure of the nanocomposite by encapsulating the NR (from step 1.1.1) inside the NT (from step 1.1.2).
      NOTE: The encapsulation can be finished by adjusting the Cartesian coordinates of the NR and the NT10,14,15.
    4. Use the vmd software to check the atomic structure before submitting the calculation job.
      1. Type vmd on the terminal in the Ubuntu system. In the opened vmd main window, select the menu File/New Molecule and find the POSCAR file through the Browse window. Load POSCAR by typing VASP_POSCAR.
      2. Display the structure in different styles in the Graphical Representations/Drawing Method window.
        NOTE: For example, once the CPK is chosen, each atom (bond) is represented by a sphere (stick). The installation guide and full tutorial of vmd are available at http://www.ks.uiuc.edu/Research/vmd.
  2. Type qsub job.pbs on the terminal in the Linux system to submit the job to the computer cluster.
    NOTE: "job.pbs" represents the name of the PBS script. The name of the PBS script is an arbitrary definition. The four input files together with the PBS script should be at the working directory. The command qsub job.pbs will be used in steps 2.2, 3.2, 3.5.2, and section 4. An example of a PBS script can be found in the supplementary coding file. After the submitted job is finished, if "reached required accuracy - stopping structural energy minimization" appears at the end of the output log, the converged result is obtained. The resulting CONTCAR file will be used as the input file POSCAR in steps 2.1, 3.1, 3.5.1, 3.5.3.1, 4.1.1, 4.1.4, and section 4.2.

2. Calculate the encapsulation energy.

  1. Type mkdir nanocomposite isolated-nanoribbon isolated-nanotube to create three folders for the nanocomposite, the NR, and the NT on a terminal in a Linux system. Prepare one PBS script "job.pbs" and four input files INCAR, POSCAR, POTCAR, and KPOINTS for the energy calculation in each folder.
    NOTE: The input file POSCAR is the file named CONTCAR with the relaxed structure from step 1. Examples of input files are given in Supplementary Materials (except POTCAR).
  2. Go to each folder and type qsub job.pbs on the terminal in the Linux system.
    NOTE: The three submitted jobs will perform the static self-consistent energy calculations for the nanocomposite, isolated NR, and isolated NT, respectively.
  3. Extract the total energy from the file OUTCAR for each system after finishing the static self-consistent calculations. Type grep “free energy TOTEN” ./nanocomposite/OUTCAR | tail -n 1, grep “free energy TOTEN” ./isolated-nanoribbon/OUTCAR | tail -n 1, and grep “free energy TOTEN” ./isolated-nanotube/OUTCAR | tail -n 1. Define the three displayed values as ENT+NR, ENR, and ENT, respectively. Calculate the encapsulation energy per angstrom: EL = (ENT+NR - ENT -ENR)/L14,15.
    NOTE: The periodical direction in each system is along the Z axis and L is the lattice constant of the unit cell along the Z axis. Test calculations of the energy dependence on the plane wave cutoff energy and the k point mesh are needed. The encapsulation energy can be used as an estimate for the energetic stability of the nanocomposite.

3. Extract the electronic properties from the band structure.

  1. Prepare one PBS script "job.pbs" and six input files: INCAR, POSCAR, POTCAR, KPOINTS, CHGCAR, and CHG for band calculation. Set ICHARG = 11 in INCAR.
    NOTE: The preconverged CHGCAR and CHG files are from the static self-consistent calculations in step 2.2. The band analysis is at the PBE level. The k point sampling in the KPOINTS file is in line-mode. Examples of input files for this step can be found in the Supplementary Materials (except POTCAR).
  2. Type qsub job.pbs on the terminal in the Linux system to submit the job.
  3. Use P4VASP to generate the projected band.
    1. Load "vasprun.xml" by typing p4v vasprun.xml on terminal in the Ubuntu system.
      NOTE: "p4v" is used to start P4VASP. The file "vasprun.xml" should be at the working directory.
    2. Select the menu Electronic/Local DOS+bands control and then Select/Bands.
      1. Specify the atomic numbers of NT in the section Atom selection. Get the atomic number by pointing to the corresponding atoms using vmd as mentioned in step 1.1.4. Specify the color, type, and size of the symbol for the projected band structure through the menu Symbol and Symbol size. Press the menu Add new line.
        NOTE: The graph will show the band structure with contributions from the NT.
      2. Repeat the same procedure following step 3.3.2.1 to get the projected band with contributions from the NR.
    3. Select the menu Graph/Export. Export the graph into a file with an agr format (for example, as "11-4.agr").
      NOTE: The output data of the projected bands by P4VASP are in three columns where the third one represents the weighting.
  4. Use xmgrace to edit the projected band.
    1. Type xmgrace 11-4.agr on the terminal to start xmgrace in the Ubuntu system. Select the menu Plot/Axis properties to edit the label and range of the axis.
    2. Select the menu Plot/Set appearance to read the energy value at the specified band number and k point.
      NOTE: The valence band maximum (VBM) and conduction band minimum (CBM) of NR/NT can be read from the projected band with contributions on NR/NT, respectively. According to the band alignments, heterostructures can be classified into three types: type I (VBMNT <VBMNR <CBMNR <CBMNT or VBMNR <VBMNT <CBMNT <CBMNR), type II (VBMNT <VBMNR <CBMNT <CBMNR or VBMNR <VBMNT <CBMNR <CBMNT), or type III (VBMNT <VBMNT <CBMNR <CBMNR or VBMNR <VBMNR <CBMNT <CBMNT)9.
    3. Calculate the valence band offset (VBO), conduction band offset (CBO), and the band gap following Kang et al.16.
    4. Select the menu File/Print to export the graph with eps format.
  5. Calculate the band decomposed charge density for the VBM and CBM.
    1. Prepare one PBS script "job.pbs" and seven input files: INCAR, POSCAR, POTCAR, KPOINTS, WAVECAR, CHGCAR, and CHG. Specify the band numbers for the CBM and VBM by the tag IBAND in INCAR. Use the single corresponding k point for each band edge.
      NOTE: The preconverged CHGCAR, CHG, and WAVECAR files are from the static self-consistent calculations in step 2.2. Examples of input files for this step are given in the Supplementary Materials (except POTCAR).
    2. Type qsub job.pbs on the terminal in the Linux system to submit the job.
    3. Use vmd to plot the VBM and CBM in real space after the job is finished.
      1. Start a vmd session and load the POSCAR file as in step 1.1.4.
      2. Select the menu File/New Molecule in the vmd main window. Find the PARCHG file through the Browse window. Load PARCHG by typing VASP_PARCHG.
      3. Select the menus Draw/Solid Surface and Show/Isosurface in the Graphical Representations window. Change the isovalue to an appropriate value (for example, 0.02). Change the color of the isosurface through the menu Coloring Method.
        NOTE: This is an intuitive analysis for band types with respect to that in step 3.4. Generally, the atomic structure is arranged away from the boundary, otherwise the visualized charge density is not shown in a continuous manner. Please see Supplemental Figure 1 for details.

4. Modulate the electronic properties of the nanocomposite (NT encapsulated inside NR) by external fields.

  1. Add a transverse electric field to the nanocomposite17.
    1. Prepare one PBS script "job.pbs" and four input files: INCAR, POSCAR, POTCAR, and KPOINTS.
    2. Define the strength of the electric field by the tag "EFIELD" in units of eV/Å.
    3. Set LDIPOL = T. Specify IDIPOL with an exact value (1, 2, or 3).
      NOTE: These two tags are added to include dipole corrections. The electric field will be applied along the X, Y, or Z axis by setting the value of IDIPOL to 1, 2, or 3.
    4. Perform the static self-consistent calculations and band structure calculations following sections 2 and 3 without structural optimization.
      NOTE: Previous studies indicate that electric fields over 5 V/Å can be used to modify the band gap of BN-NT and BN-NR without deforming the structure18,19.
  2. Add a longitudinal tensile strain to the nanocomposite.
    1. Change the lattice parameters along the periodical direction to reflect the strain effect.
      NOTE: For example, the optimized lattice parameter of the nanocomposite along the Z axis is 2.5045 Å. If 1% uniaxial tensile strain is applied along the Z direction, change the lattice parameter in POSCAR to 2.5045 x 1.01 = 2.529545 Å.
    2. Relax the modified structure following section 1.
    3. Perform static self-consistent calculations and band structure calculations following sections 2 and 3.

Access restricted. Please log in or start a trial to view this content.

Results

Zigzag BN-NRs encapsulated inside armchair BN-NTs (11,11) were chosen as representative examples for a 1D vdW heterostructure. The lattice parameters were taken from Sahin et al.20. For convenience, zigzag NRs are abbreviated Zn, where n represents the III–V dimers along the width14. The encapsulation energy EL from step 2.3 was used as a rough estimate for the energetic stability of the nanocomposite. The EL values of Z2, Z3...

Access restricted. Please log in or start a trial to view this content.

Discussion

The calculations for electronic properties in sections 2, 3, and 4 would be similar among various nanoscale materials. The initial atomic model in step 1 should be carefully designed to extract meaningful information. For example, the factor for selecting the model could be the size or chirality of the materials. Also, the initial atomic model in step 1.1 should be reasonably prepared for low-cost structure relaxation. Taking the nanocomposite in the protocol as an example, the NR should be encapsulated inside the NT in ...

Access restricted. Please log in or start a trial to view this content.

Disclosures

The authors have nothing to disclose.

Acknowledgements

This work was supported from China Postdoctoral Science Foundation (Grant No. 2017M612348), Qingdao Postdoctoral Foundation (Grant No. 3002000-861805033070) and from the Young Talent Project at Ocean University of China (Grant No. 3002000-861701013151). The authors thank Miss Ya Chong Li for preparing the narration.

Access restricted. Please log in or start a trial to view this content.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
Nanotube ModelerDeveloped by Dr. Steffen WeberNanotubeModeler1.8http://www.jcrystal.com/products/wincnt/NanotubeModeler.exe
P4VASPOrest Dubayp4vasp 0.3.30Open source, available at www.p4vasp.at
v2xsfDeveloped by Dr. Jens Kunstmannv2xsfhttp://theory.chm.tu-dresden.de/~jk/software.html
VASP softwareComputational Materials Physics, Dept. of Physics, University of Viennavasp.5.4.1https://www.vasp.at
VMD softwareTheoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaignvmd1.9.3https://www.ks.uiuc.edu/Research/vmd
xcrysdenDept. of Physical and Organic Chemistry, Jozef Stefan InstituteXCrySDen1.5.60http://www.xcrysden.org/
XmakemolDeveloped by M. P. Hodgesxmakemol5.16https://www.nongnu.org/xmakemol/XmakemolDownloads.html
Xmgrace softwareGrace Development Team under the coordination of Evgeny Stambulchikxmgrace5.1.25http://plasma-gate.weizmann.ac.il/Grace/

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).

Access restricted. Please log in or start a trial to view this content.

Reprints and Permissions

Tags

Density Functional TheoryVASP SimulationBand Structure AnalysisEncapsulation EnergyElectronic PropertiesElectric Field ModulationTensile Strain Engineering
Probe Type II Band Alignment in One-Dimensional Van Der Waals Heterostructures Using First-Principles Calculations
Contact Us
Recommend to Library

Privacy

Terms of Use

Policies

Copyright © 2026 MyJoVE Corporation. All rights reserved.