Questions tagged [vasp]

Questions about or related to the Vienna ab-initio Simulation Package (VASP).

Filter by
Sorted by
Tagged with
12 votes
5 answers
1k views

Create doped structures to POSCAR files for vasp

I have a perfect supercell structure Ni fcc. And I would like to dope another element O or N with a ratio of O:Ni=8:100 to the supercell of 3x3x3 and find the most stable ordering structure. The ...
Binh Thien's user avatar
  • 2,025
29 votes
2 answers
5k views

Quantum ESPRESSO vs VASP

We know that both Quantum ESPRESSO and VASP can perform lots of different kinds of simulations of materials. But what can VASP do that Quantum ESPRESSO can not, and vice-versa?
Anibal Bezerra's user avatar
21 votes
3 answers
4k views

k-points and ENCUT convergence tests before or after relaxation?

Let's say I want to relax a structure using VASP. K-points and ENCUT convergence tests must be conducted before or after relaxation? and which tags should be ...
Chi Kou's user avatar
  • 5,683
15 votes
3 answers
3k views

Geometry optimization for n-layers in VASP

What is the best way to optimize monolayer geometry in VASP? Should the same method be used for an n-layer (n=2,3,4,5)?
M_Fisher's user avatar
  • 601
10 votes
2 answers
409 views

How can I be sure that I have found the conduction band minimum in DFT?

Let's say I want to find the electronic energy bands from DFT calculations in VASP along a chosen path in the Brillouin zone. Now, I want to make sure I choose a path that captures the conduction band ...
B. Brekke's user avatar
  • 375
6 votes
1 answer
515 views

Can a file in the Gaussian CUBE format be voxelized?

My aim is to voxelize a CHGCAR file. I am trying to build up a workflow to do the same. With the aid of the Matter Modeling community members, I have figured out a way to convert a CHGCAR file into ...
Pranoy Ray's user avatar
  • 1,481
28 votes
2 answers
1k views

What are good ways to reduce computing time when working with large systems in VASP?

I know that parallelization is an important factor, however, I wanted to know if there are any other methods of reducing computing time when dealing with large systems without sacrificing accuracy. ...
Daniel M.'s user avatar
  • 821
22 votes
2 answers
7k views

What are some good band-structure/DOS plotting tools/styles?

I am currently using p4vasp for drawing DOS and band-structure plots, but it is not much handy in modifying the OUTPUTS. Could you please suggest any other good software or tools for that?
Chi Kou's user avatar
  • 5,683
21 votes
9 answers
9k views

How to master Fortran with minimal effort?

Fortran language is used in many first-principles matter modeling codes, such as VASP and QE. These codes usually include many ...
Jack's user avatar
  • 14.9k
20 votes
2 answers
5k views

What is the difference between Ultrasoft, ONCV and PAW Pseudopotentials? Which is better for a spin-orbit coupled calculation?

I am trying to do spin-orbit coupled calculations for various topological insulators. I have found papers using Quantum Espresso with ONCV pseudopotentials and papers using VASP with PAW ...
Ashique Lal's user avatar
  • 1,581
19 votes
1 answer
2k views

What are the pre/post processing tools available for VASP?

I am aware of pyPROCAR and VASPKIT. Are there any other tools supporting VASP
Thomas's user avatar
  • 8,922
17 votes
1 answer
1k views

Are there problems with the Tran-Blaha modified Becke-Johnson potential (mBJ, a.k.a. TB09) in Quantum ESPRESSO?

I have used the modified Becke-Johnson (mBJ) exchange potential before, as it was developed and implemented in Wien2K by Fabien Tran and Peter Blaha (1). In 2013, Éric Germaneau et al. implemented mBJ ...
epalos's user avatar
  • 4,719
13 votes
4 answers
3k views

How to make a k-point convergence test for a bulk structure with a lower k-point in the z-direction?

Let's say I want to make k-point convergence test for graphite. And let's say it converges at 12x12x4. Do I first need to do a k-point convergence test for equal k-points on xyz directions then make ...
Taavi's user avatar
  • 131
11 votes
2 answers
2k views

Parallelization in VASP

I have a system that has 1 k-point, 54 atoms, NBANDS=152, and NGZ=80. I'm running an AIMD simulation on VASP version 5.4.4 The ...
DoubleKx's user avatar
  • 1,060
11 votes
1 answer
1k views

Simulated low energy electron diffraction (LEED) patterns

What modern software supports the generation of LEED patterns? It seems most LEED simulation packages are quite old and do not work on modern hardware. Is there something that prevents common ...
Tristan Maxson's user avatar
10 votes
1 answer
2k views

How to carry out HSE06+SOC band structure calculation using VASP?

I could not find any good source for this type of calculation. I'd be very appreciated if someone who is experienced with that to show it step by step
Savir's user avatar
  • 1,041
9 votes
2 answers
1k views

How to process electron-density data from a CHGCAR for better handling?

We all know that Vesta helps us visualize the electron clouds present in the CHGCAR data, but I am looking for a better way to process the electron-density data we get in a CHGCAR file for compact ...
Pranoy Ray's user avatar
  • 1,481
8 votes
1 answer
540 views

How to evaluate the electron relaxation time of nanostructures computationally

I am working with Si nanowires using the SIESTA ab initio code and the BoltzTraP2 Boltzmann transport equation solver. However, the output of the BoltzTraP2 code comes per relaxation time (...
PBH's user avatar
  • 2,603
6 votes
1 answer
622 views

Different between ENCUT and ENAUG in the VASP setting?

I am trying to figure out the difference between ENCUT and ENAUG in VASP settings, I found the following description on the wiki of VASP. ENAUG specifies the cut-off energy of the plane wave ...
Jack's user avatar
  • 1,977
5 votes
1 answer
925 views

Best practice(s) for q-points convergence in phonon density of states calculations

The answer can be a complete description of all the things we must be attentive to while converging q-points. A peculiar thing that I feel contradicted about: Generally, we don't take an equal # of k-...
Hitanshu Sachania's user avatar
4 votes
1 answer
514 views

Calculation of DOS in HSE06

Do I do rest, scf, dos steps when calculating DOS like PBE in HSE06? Can anyone help me with this?
Davron Kurbanov's user avatar
2 votes
0 answers
139 views

How to include the spin-orbit coupling strength of only two or three elements in VASP? [closed]

I am trying to calculate magnetic properties in my material with VASP, much like Fig 3b of this paper. The paper says: "To understand the origin of the SIAin monolayer MnBi2Te4, we performed the ...
liao guo's user avatar
19 votes
1 answer
2k views

Difference between potential energy, free energy and Coulomb energy in solid state physics

I often encounter terms such as (Helmholtz, Gibbs) free energy, potential energy and total energy when describing the energy of a physical system at atomic level. Sometimes I stumble upon Coulomb ...
Achintha Ihalage's user avatar
18 votes
2 answers
483 views

Calculating binding energy between two systems

Regarding the calculation of binding energy of molecular absorption on graphene, for which this has been discussed before, the process of calculating $E_{AB}$ is fairly straight forward, i.e., it is ...
Shahid Sattar's user avatar
17 votes
1 answer
669 views

How to ensure a smooth band structure?

I performed an HSE06 calculations using VASP and the calculated band gap is in good agreement with experiment. However the band is not smooth, with very sharp, zig-zag like features. How to improve ...
nisarmuhammad131's user avatar
17 votes
1 answer
1k views

How does the mBJ potential give an accurate band gap as compared to LDA and GGA?

I am not able to understand the literature on how the modified Becke-Johnson(MBJ)$^{[1]}$ potential gives an accurate bandgap. Can someone please help? Thank you. The formula for the potential can be ...
sonia rani's user avatar
16 votes
2 answers
2k views

How to start with the optical properties of 2D materials with the VASP code?

I hope you are doing well. I am new to VASP, and my project is "Optical properties of 2D materials on VASP" but I have no idea how to run a simple calculation with VASP. I will be very ...
Haris's user avatar
  • 229
15 votes
1 answer
598 views

Standard values for level of convergence

In matter modelling we use self consistent approaches for calculating different properties. Many DFT codes find the ground state of the system by determining total energy self consistently. Going ...
Thomas's user avatar
  • 8,922
14 votes
1 answer
1k views

Is LDA better than GGA for elastic constants?

I have calculated elastic constants for Si, GaAs, and GaN using LDA and GGA in VASP. LDA is giving better results than GGA for elastic constants. Is there is a reason behind the good results by ...
sonia rani's user avatar
14 votes
3 answers
2k views

How can I use a GUI on a supercomputer?

I am trying to make plots based to VASP simulation results, but I am struggling with how to download so much data from the server. A lot of examples on the internet are taught using the local desktop ...
Jack's user avatar
  • 1,977
13 votes
2 answers
3k views

mBJ-LDA band structure calculation using VASP

I tried to create a band structure with the mBJ-LDA method using VASP. I found information that the KPOINTS file should be created according to the "0-weight (...
M_Fisher's user avatar
  • 601
13 votes
1 answer
920 views

Doubt about energy convergence with different k points

During the convergence test of k points, I observed my energy start increasing instead of decreasing. I am doubtful is this correct or I am doing something wrong. can I take the appropriate K point in ...
poonam sharma's user avatar
12 votes
1 answer
1k views

Why is the Nonequilibrium Green’s function (NEGF) method for quantum transport not implemented in plane wave codes like VASP and Quantum Espresso?

DFT codes like Siesta, QuantumATK, Questaal, and Smeagol having localized basis sets have NEGF implementations available. But not in plane wave codes like VASP and Quantum Espresso. Why is it so?
Thomas's user avatar
  • 8,922
12 votes
1 answer
497 views

Dealing with symmetry of ordered primitive cell during DFT structure relaxation

I have a large disordered cubic structure that I would like to do AIMD on to look at cation diffusion. I found the primitive cell and then ordered it to create several configurations. I will then ...
DoubleKx's user avatar
  • 1,060
12 votes
0 answers
146 views

Correlation Functions in the Cluster Variation Method [closed]

I've been referring to Dr Mohri's paper [DOI: 10.1007/s11837-013-0738-5] for the cluster variation method (CVM). I wish to calculate the configurational entropy of a binary FCC system. Cluster ...
Hitanshu Sachania's user avatar
12 votes
1 answer
542 views

Binding energy calculation for Zn atom on armchair graphene slab?

I am trying to repeat the result of this paper for Zn absorbed on the armchair graphene slab. The relaxed top and side structures are the following: The relaxed structure is almost the same as the ...
Jack's user avatar
  • 14.9k
11 votes
1 answer
450 views

How to start with structural defects in monolayer?

First of all, thank you for your help! You are so helpful every time. I would like to calculate the influence of the structural defects on the electronic structure in the HfS2 monolayer. With VESTA ...
M_Fisher's user avatar
  • 601
11 votes
2 answers
729 views

Which properties of a material calculated using density functional theory are affected by smearing?

Smearing (smearing width $\sigma$, to be precise) has always been confusing. I understand what it is but can't pin it down when it comes to DFT calculations. Should it be converged like k-points and ...
Hitanshu Sachania's user avatar
10 votes
1 answer
537 views

How to remove the fake k-points from vasprun.xml file in the calculation of HSE06 band structure?

Is there any script that can remove fake k-points from a vasprun.xml file for HSE06 band structure calculation? Because removing it manually is a time-consuming ...
Chi Kou's user avatar
  • 5,683
10 votes
2 answers
1k views

How to make a good HSE06 bandstructure for TMDC monolayer material (e.g.: WS2) in terms of VASP calculation INPUTS?

I want to perform HSE06 band-structure calculation of a WS2 monolayerwhich has the following POSCAR with ...
Chi Kou's user avatar
  • 5,683
10 votes
1 answer
929 views

How to calculate the formation energy of a monolayer using VASP?

I want to calculate the formation energy of WS2 monolayer unit cell using VASP. I am wandering if you could show me the detailed ...
Chi Kou's user avatar
  • 5,683
9 votes
1 answer
1k views

Electric Field Calculations in VASP

Can someone share appropriate tags and method to apply Electric Field in VASP?
Shahid Sattar's user avatar
9 votes
1 answer
2k views

How can we calculate the Gibbs free energy from first principles?

I am trying to calculate the Gibbs free energy for a bulk material. Since DFT ground state calculations give the energy at 0K, a temperature dependent part has to be added in this to obtain correct ...
Shahid Sattar's user avatar
9 votes
2 answers
493 views

Problem related to the calculation of elastic constant with VASP5.4.4

I am using VASP5.4.4 to calculate the mechanical properties of the cubic phase SnTe. The structure can be found in the materials project website. In detail, I use the following INCAR to calculate the ...
Jack's user avatar
  • 14.9k
9 votes
1 answer
981 views

How to extract lattice constants a and c from relaxed hexagonal structure using VASP?

Taking a bulk WS$_2$ (bilayer) structure and relaxing it with ISIF=3 and IVDW=10 (DFT-D2 method) using the VASP code. Then ...
Chi Kou's user avatar
  • 5,683
9 votes
1 answer
528 views

How does charge mixing work?

When doing DFT calculations on certain magnetic materials, convergence can be difficult. This problem can be solved by changing the mixing parameters (e.g. AMIN and BMIX in the INCAR file of VASP). My ...
user avatar
9 votes
2 answers
1k views

What is the format of the electron-density data in a CHGCAR file?

Rephrasing the above, what is the pattern in which the electron-density data is stored in the CHGCAR file? I understand that the electron density data in a CHGCAR file is in the form of a 3D Fortran ...
Pranoy Ray's user avatar
  • 1,481
9 votes
1 answer
641 views

Work Function Problem

2D materials (such as graphene) are mostly exfoliated from their 3D bulk counterpart. I am wondering how we can calculate work function for a 2D material? Normally electrostatic potential in the out-...
Shahid Sattar's user avatar
8 votes
1 answer
500 views

When to use pseudohydrogens to passivate a slab?

I often see some papers that don't use pseudohydrogen passivating the slab's bottom side when doing DFT simulation on the metal oxide's surface, for example, the following description is from a paper ...
Jack's user avatar
  • 1,977
8 votes
2 answers
484 views

Why don't we use the principal quantum number when building the projected density of state?

I am using VASP to calculate the projected DOS, but I am not sure why we aren't using the principal quantum number when doing the projection. According to the VASP wiki, LORBIT would control VASP to ...
Jack's user avatar
  • 1,977