Effective Mass Calculation Using Vasp

Effective Mass Calculator for VASP Users

A precision tool for calculating carrier effective mass from VASP band structure data.

Carrier Mass Calculator

Energy E₁ (at k – Δk)

Energy of the first point in electronvolts (eV).

Energy E₂ (at band extremum k)

Energy of the band minimum or maximum (eV). This is your reference point.

Energy E₃ (at k + Δk)

Energy of the third point in electronvolts (eV).

k-point step (Δk)

The distance between k-points in reciprocal Ångströms (Å⁻¹).

Please ensure all inputs are valid numbers and Δk is not zero.

Calculation Results

Effective Mass
1.000 m₀

d²E/dk²
0.00 eV·Å²

Carrier Type
Electron

Mass in kg
0.00e+0 kg

Parabolic Fit Check
Symmetric

Visualization of the E-k data points and the fitted parabolic band.

What is an effective mass calculation using VASP?

An effective mass calculation using VASP data is a fundamental procedure in computational materials science. It determines how a charge carrier (an electron or a hole) accelerates in a crystal lattice under an external electric field. It’s not the particle’s true rest mass but a value that represents its inertia within the periodic potential of a solid. VASP (Vienna Ab initio Simulation Package) provides the necessary electronic band structure (the E-k diagram), from which effective mass is derived.

This value is crucial for predicting a material’s electronic properties, such as carrier mobility and electrical conductivity. A low effective mass implies the carrier can move easily through the crystal, leading to high mobility, which is desirable for applications like transparent conductors or high-frequency transistors. Understanding the semiconductor properties is a key part of this analysis.

The Effective Mass Formula and Explanation

The effective mass (m*) is inversely proportional to the curvature of the energy band. For a simple parabolic band, it is defined by the second derivative of energy (E) with respect to the wave vector (k).

m* = ħ² / (d²E / dk²)

Since VASP outputs discrete E-k points rather than a continuous function, we approximate the second derivative using the finite difference method. For three equally spaced points, this approximation is:

d²E/dk² ≈ (E(k-Δk) – 2·E(k) + E(k+Δk)) / (Δk)²

The accuracy of this effective mass calculation using VASP data depends heavily on the density of the k-point mesh around the band extremum. A related concept you might explore is the density of states.

Formula Variables
Variable	Meaning	Unit (in this calculator)	Typical Range
m*	Effective Mass	Ratio to electron rest mass (m₀)	0.01 – 10 m₀
ħ	Reduced Planck Constant	eV·s	6.582 x 10⁻¹⁶ eV·s
E	Energy Eigenvalue	eV	-10 to +10 eV (relative to Fermi level)
k	Wave Vector	Å⁻¹	0 to ~1 Å⁻¹ (within first B.Z.)
Δk	k-point spacing	Å⁻¹	0.001 – 0.05 Å⁻¹

Practical Examples

Example 1: Electron Effective Mass in Silicon (Si)

Let’s calculate the longitudinal effective mass of an electron at the conduction band minimum (CBM) of Silicon. We use hypothetical but realistic data points from a VASP calculation.

Input E₁: 0.0098 eV
Input E₂ (CBM): 0 eV
Input E₃: 0.0098 eV
Input Δk: 0.01 Å⁻¹

Calculation:

d²E/dk² = (0.0098 – 2*0 + 0.0098) / (0.01)² = 196 eV·Å²

m*/m₀ = 7.61996 / 196 ≈ 0.98 m₀. This is a classic result for the longitudinal mass in Si.

Example 2: Heavy Hole Effective Mass in GaAs

Now, let’s calculate the mass for a hole at the valence band maximum (VBM). The curvature is negative, indicating a maximum. The mass itself is positive.

Input E₁: -0.0076 eV
Input E₂ (VBM): 0 eV
Input E₃: -0.0076 eV
Input Δk: 0.02 Å⁻¹

Calculation:

d²E/dk² = (-0.0076 – 2*0 – 0.0076) / (0.02)² = -38 eV·Å²

m*/m₀ = 7.61996 / |-38| ≈ 0.50 m₀. This corresponds to the heavy-hole band in GaAs.

How to Use This Effective Mass Calculator

Follow these steps to perform an accurate effective mass calculation using VASP results:

Perform VASP Calculation: First, run a high-resolution band structure calculation in VASP along the desired high-symmetry lines.
Identify Band Extremum: Plot the band structure and visually identify the Conduction Band Minimum (CBM) for electrons or the Valence Band Maximum (VBM) for holes. Note the energy of this point.
Extract Data Points: From your `EIGENVAL` file or plotted data, find two additional points, one on each side of the extremum, that are equally spaced in k-space. Record their energy values (E₁ and E₃).
Determine Δk: Find the distance in reciprocal space between your central point (E₂) and the outer points (E₁ or E₃). This value can be found from your KPOINTS file or OUTCAR. This is your Δk.
Enter Values: Input E₁, E₂, E₃, and Δk into the calculator fields above.
Interpret Results: The calculator automatically provides the effective mass as a ratio to the free electron mass (m₀), the carrier type, and a visualization of the parabolic fit. Proper interpretation is key, much like understanding band gap analysis.

Key Factors That Affect Effective Mass

Several physical and computational factors influence the final calculated value:

Crystal Symmetry: The effective mass is a tensor, meaning its value depends on the direction of measurement in k-space (e.g., Γ→X vs. Γ→L). Our calculator handles one direction at a time.
Chemical Composition: Alloying or doping a material directly alters the crystal potential, which reshapes the energy bands and changes their curvature.
Lattice Strain: Applying tensile or compressive strain distorts the unit cell, breaks symmetries, and can drastically change the effective mass, a technique used in “strain engineering.”
Quantum Confinement: In nanomaterials like quantum wells or nanowires, the band structure is modified, leading to different effective masses compared to the bulk material. For advanced topics, consider exploring the Wannier90 interface.
VASP Functional: The choice of DFT functional (e.g., LDA, GGA, HSE) affects the predicted band gap and band curvature. HSE functionals are often more accurate for this purpose.
k-point Density: The precision of the finite difference approximation is highly dependent on a dense k-point sampling around the band extremum. A coarse mesh will lead to an inaccurate effective mass calculation using VASP.

Frequently Asked Questions (FAQ)

Why is my calculated mass negative?: A negative sign in the calculation arises from a negative curvature (d²E/dk² < 0), which defines the top of a band. This indicates a hole. The calculator automatically reports the carrier type as "Hole" and provides the positive mass value.
What is the unit of the final result?: The primary result is a dimensionless ratio of the calculated mass to the rest mass of a free electron (m*/m₀). This is the standard convention in solid-state physics.
How do I get E and k data from VASP?: Energy eigenvalues are in the `EIGENVAL` file. The corresponding k-points and their distances are detailed in the `OUTCAR` file (search for “reciprocal lattice vectors”) or can be derived from the `KPOINTS` file for a simple path.
What if my three points are not perfectly symmetric?: This calculator assumes a symmetric parabolic shape. If E₁ and E₃ are different, it indicates an asymmetric band. The tool will still provide a fit, but it’s a sign that the band is non-parabolic, and the effective mass may be energy-dependent. For higher accuracy, you may need to use a more sophisticated fitting code. Asymmetry can be explored in phonon calculations as well.
Can this tool calculate the full effective mass tensor?: No. This is a simplified tool that calculates the scalar effective mass along one user-specified direction. Calculating the full tensor requires fitting the band structure in 3D, typically done with specialized scripts.
What is a “heavy” vs. “light” carrier?: These terms refer to carriers with relatively large or small effective masses. For example, in many materials, the valence band splits into “heavy hole” and “light hole” bands, with different curvatures and thus different masses.
How important is the Δk value?: It is critically important. A Δk that is too large will not accurately capture the curvature at the band edge. A Δk that is too small can be subject to numerical noise from the DFT calculation. Choosing a value that corresponds to 2-4 points away from the extremum in a high-resolution calculation is a good starting point.
Why is the E-k chart important?: The chart provides an essential visual check. It confirms that your three chosen points actually form a parabola around a minimum or maximum. If the points do not look parabolic, your data is likely unsuitable for this simple fitting model.

Related Tools and Internal Resources

For a deeper dive into materials simulation, explore these related resources and tools that complement your effective mass calculation using VASP.

Lattice Dynamics and Phonon Frequencies: Understand the vibrational modes of your crystal.
Optical Properties Simulation: Calculate dielectric functions and absorption spectra.
VASP Input File Generator: Tools to help set up complex calculations automatically.