# Tag Info

5

No, this isn't possible for an undoped semiconductor. The reason it's happening is because the Fermi level is not well defined for your semiconductor at finite temperature, and you're trying to read the valence band maximum (VBM) from the graph of your broadened DOS. I recommend taking the VBM directly from the band energies, instead. In a DFT program, the ...

4

This seems more of a thermodynamics question than a simulations question given the comment from May 12 '21. What is set to zero by definition is the Gibbs free energy of formation at standard temperature and pressure (1 atm, 298.15 K) for pure elements in their natural state. This explains the difference between the cited values for the solid and the liquid (...

4

When dealing with monolayers, I have found that using the ISIF=4 tag is better and is less time-consuming. As mentioned in Vaspwiki, it allows the relaxation of both lattice constants and atomic coordinates with a fixed volume. However, what is important here is uncommenting the FCELL(3,3)=SAVE line in constr_cell_relax.F file and recompiling VASP. This ...

4

I could be wrong, but I am fairly certain that line of the CONTCAR technically doesn't matter. If you put a O instead of H, but put the potentials in the right order, you will get the right result even though the element is wrong. For visualization, it will definitely be wrong but I don't think the calculated result is wrong. In tools such as ASE the ...

6

I think this issue is avoided in practical sense. I am not sure any of the psuedopotentials provided with VASP actually include valence electrons of the same angular momentum (so you never have 3s and 4s as valence), since these states are normally part of the core. I don't think there is any reason you cannot do this in theory, but in practice the core is ...

4

If by "standard", you mean an "official" VASP-supported program, then the answer is no. VASP gives you the raw text output for you to deal with on your own. That being said, the vasprun.xml file produced at the end of a run contains many of the important calculation results, so any parser than can read an XML file can also parse this data ...

8

TL;DR If bonds are broken in creating the lower surface of the slab, then passivation is probably a good idea. Let's start with a little bit of context: Slab calculations in computational materials science are often (but not always) designed to study the properties of a surface of a macroscopic crystal. In these cases, an accurate model of physical reality ...

5

I know this was long ago, but for those that find the same problem in the future, the sisl python package might be useful. It has a generic Grid object and it knows how to read the grids from VASP: import sisl grid = sisl.get_sile("path/to/your/CHGCAR").read_grid() Then grid is a sisl Grid object, so: grid.grid contains the numpy array of values ...

4

# # To run VASP this script calls $vasp_std # (or posibly$vasp_gam and/or $vasp_ncl). # These variables can be defined by sourcing vaspcmd . vaspcmd 2> /dev/null As mentioned in the comments, the answer is actually in your question. vaspcmd is just any script you provide that defines the$vasp_std, $vasp_gam and$vasp_ncl variables which give the full ...

7

The community wiki by Tyberius is a bit misleading here, you are not doing molecular dynamics calculations (IBRION=0) but likely a geometry optimization (IBRION=1,2,3). In the case of geometry optimization, the CHGCAR is not the predicted charge density, but is instead the charge density of the last completed step. For this reason, you should find no ...

5

In the passage just above what you quoted, it says: For dynamic simulation (IBRION=0), the charge density on the file is the predicted charge density for the next step: i.e. it is compatible with the CONTCAR file, but incompatible with the last positions in the OUTCAR file. This means while the density in the CHGCAR file is a density for the position given ...

2

With a large enough vacuum, you don't need a dipole correction as badly to get good results. Sometimes the correction can cause convergence issues or problems with single steps giving incorrect energies. The second calculation looks like a geometry optimization, so maybe this was done to improve stability; a dipole correction should be applied after ...

1

This is related to your other question, you are correct VASP only supports fixing atoms along the X, Y, and Z axis. You can use ASE optimizers with the FixLine or FixPlane constraints to do this. You just specify which atoms are constrained and which direction they are constrained to move in. Just be aware, ASE running VASP is slower than VASP running its ...

0

The FixAtoms constraint will be turned into the selective dynamics tags as "F F F". ASE will also convert selective dynamics to this constraint when reading. I think the FixedPlane and FixedLine could in theory also work, but I am not sure they do since they support directions not aligned with the X, Y, and Z axis. If you need this and it doesn't ...

5

ASE in general picks the VASP default, by simply not including the tag if not directly asked for. This is done by setting the tag to "None" which has no meaning in VASP. Then when the INCAR is generated, the tag is simply not output, leaving VASP to fall back on its default behavior. For things such as the functional, ASE does attempt to set up ...

5

I don't know any VASP-specific details, but the basic answer is that the three versions are almost identical, and solve the same equations in the same way. If we take vasp_std as the reference, then the differences are: vasp_gam If your system has time-reversal symmetry (true for most simulations), then $+k$ and $-k$ are symmetry-related and we can restrict ...

3

the example doesn't generate K-points from the slab, it uses the K points extracted from the Brillouin zone of the bulk primitive cell. Not sure you're reading the figure correctly. The bottom of the figure shows the three-dimensional Brillouin zone of the primitive cell, but the top of the figure (which the G-M-X-G path refers to) shows the two-dimensional ...

4

Following the suggestion by Ian, let's first translate SIGMA into electronic structure language: VASP's SIGMA parameter defines the smearing width of the Fermi-Dirac [1] occupation function f(E) in electron volts. There are at least two different reasons for choosing a finite smearing: You want to simulate a finite electronic temperature. You are dealing ...

7

You should take Binh Tien's advice for the continuity part of your problem. It has also been answered somewhat here: How to ensure a smooth band structure? As for the gap issue, adding to Phil Hasnip's answer, a similar question has also had several proposals for answers here : How to generate the high symmetry paths for band structure calculations? These ...

8

This is actually an important question regardless of the code you're using (I use SIESTA and NWChem). To obtain a good phonon spectrum, the first thing you need is a truly relaxed/optimized structure. This requirement entails several parameter optimizations (mesh cutoff, force tolerance, basis set ..). Afterwards, other requirements are also important. I ...

10

1. Is it possible that the energy gap of the material is direct in the bulk but indirect in the slab? Yes, this is possible (as is the reverse). The 2D slab does not have the same symmetry as the full 3D slab, and there is no requirement for the conduction band minima and valence band minima to be at the same points in reciprocal-space. As you make the slab ...

4

The DOSCAR requires another run after a relaxation to be useful because the DOSCAR gives an average of the DOS across the ionic steps. In the first paragraph of the DOSCAR VASP wiki page: For dynamic simulations and relaxations, an averaged DOS and an averaged integrated DOS is written to the file. For a description of how the averaging is done see the tags ...

2

As far as my knowledge, you need to run SCF with ISMEAR = -5 to just get very accurate energies, while for relaxations you need to use ISMEAR = 0 (for semiconductors and insulators). Also, the DOSCAR file is used for the density of states calculations, and the tetrahedron smearing (ISMEAR = -5) must be toggled on to get the correct Fermi energy.

Top 50 recent answers are included