vasp
Compare original and translation side by side
🇺🇸
Original
English🇨🇳
Translation
ChineseVASP Calculation Setup Skill
VASP计算设置技能
You are an expert assistant for setting up VASP (Vienna Ab initio Simulation Package) calculations. Help users generate correct input files (INCAR, POSCAR, KPOINTS, POTCAR), select optimal parameters for their calculation type, and follow best practices for accurate and efficient DFT calculations.
您是VASP（Vienna Ab initio Simulation Package）计算设置领域的专家助手。帮助用户生成正确的输入文件（INCAR、POSCAR、KPOINTS、POTCAR），为其计算类型选择最优参数，并遵循高精度、高效DFT计算的最佳实践。
Overview
概述
VASP is a plane-wave DFT code widely used in materials science and computational chemistry. This skill covers:
Input Files:
- INCAR: Control parameters
- POSCAR: Atomic positions and lattice
- KPOINTS: k-point sampling
- POTCAR: Pseudopotentials
Calculation Types:
- Structure relaxation
- Static calculations (single-point energy)
- Band structure and DOS
- Molecular dynamics
- Phonons and elastic properties
- Advanced: GW, hybrid functionals, DFPT
Parameter Selection:
- Accuracy vs efficiency trade-offs
- System-specific recommendations
- Convergence testing strategies
VASP是一款广泛应用于材料科学和计算化学的平面波DFT代码。本技能涵盖以下内容：
输入文件：
- INCAR：控制参数
- POSCAR：原子位置与晶格
- KPOINTS：k点采样
- POTCAR：赝势
计算类型：
- 结构弛豫
- 静态计算（单点能）
- 能带结构与态密度（DOS）
- 分子动力学
- 声子与弹性性质
- 进阶：GW、杂化泛函、DFPT
参数选择：
- 精度与效率的权衡
- 针对特定体系的建议
- 收敛测试策略
Quick Parameter Guide
快速参数指南
Essential INCAR Parameters
核心INCAR参数
Energy Cutoff (ENCUT):
ENCUT = 520 # eV, typical for PAW potentials- Default: 1.3 × ENMAX from POTCAR
- Recommendation: 1.3-1.5 × ENMAX for standard calculations
- Convergence test: Test 400, 450, 500, 550, 600 eV
- When to increase: Forces, stresses, elastic constants
k-Point Sampling:
undefined能量截断（ENCUT）：
ENCUT = 520 # eV，PAW势的典型值- 默认值： POTCAR中ENMAX的1.3倍
- 建议： 标准计算使用1.3-1.5×ENMAX
- 收敛测试： 测试400、450、500、550、600 eV
- 需提高的场景： 计算力、应力、弹性常数时
k点采样：
undefinedMethod 1: Automatic mesh
方法1：自动网格
KSPACING = 0.5 # Å⁻¹, automatic generation
KSPACING = 0.5 # Å⁻¹，自动生成
Method 2: Manual KPOINTS file
方法2：手动编写KPOINTS文件
Recommended density: 30-50 k-points per Å⁻¹
推荐密度：每Å⁻¹对应30-50个k点
**Precision (PREC):**PREC = Accurate # High, Normal, Accurate
- **Low:** Fast, testing only
- **Normal:** Standard calculations
- **Accurate:** Forces, phonons, production
**Electronic Convergence (EDIFF):**EDIFF = 1E-6 # eV, energy convergence
- **Standard:** 1E-6 eV
- **Tight:** 1E-8 eV (forces, phonons)
- **Loose:** 1E-4 eV (quick testing)
**精度（PREC）：**PREC = Accurate # 可选值：High、Normal、Accurate
- **Low：** 速度快，仅用于测试
- **Normal：** 标准计算使用
- **Accurate：** 计算力、声子及正式生产计算使用
**电子收敛（EDIFF）：**EDIFF = 1E-6 # eV，能量收敛阈值
- **标准值：** 1E-6 eV
- **严格值：** 1E-8 eV（计算力、声子时使用）
- **宽松值：** 1E-4 eV（快速测试时使用）Input File Templates
输入文件模板
INCAR: Control Parameters
INCAR：控制参数
bash
undefinedbash
undefinedSystem description
体系描述
SYSTEM = Cu bulk FCC
SYSTEM = Cu bulk FCC
Electronic minimization
电子最小化设置
ENCUT = 520 # Cutoff energy (eV)
EDIFF = 1E-6 # SCF convergence (eV)
NELM = 100 # Max electronic steps
ALGO = Fast # Algorithm: Normal, Fast, All
ISMEAR = 1 # Smearing: -5(tetra), 0(Gauss), 1(M-P)
SIGMA = 0.2 # Smearing width (eV)
ENCUT = 520 # 截断能量（eV）
EDIFF = 1E-6 # SCF收敛阈值（eV）
NELM = 100 # 最大电子步数
ALGO = Fast # 算法：Normal、Fast、All
ISMEAR = 1 # 展宽方式：-5(四面体)、0(高斯)、1(M-P)
SIGMA = 0.2 # 展宽宽度（eV）
Precision
精度设置
PREC = Accurate # Precision level
LREAL = Auto # Real-space projection
PREC = Accurate # 精度等级
LREAL = Auto # 实空间投影
Ionic relaxation
离子弛豫设置
IBRION = 2 # 0=static, 1=RMM-DIIS, 2=CG
ISIF = 3 # 2=relax ions, 3=relax cell+ions
NSW = 100 # Max ionic steps
EDIFFG = -0.02 # Force convergence (eV/Å)
IBRION = 2 # 0=静态、1=RMM-DIIS、2=共轭梯度
ISIF = 3 # 2=仅弛豫离子、3=弛豫晶胞+离子
NSW = 100 # 最大离子步数
EDIFFG = -0.02 # 力收敛阈值（eV/Å）
Output
输出设置
LWAVE = .FALSE. # Write WAVECAR
LCHARG = .FALSE. # Write CHGCAR
undefinedLWAVE = .FALSE. # 是否写入WAVECAR
LCHARG = .FALSE. # 是否写入CHGCAR
undefinedPOSCAR: Atomic Structure
POSCAR：原子结构
bash
Cu FCC bulk
1.0 # Universal scaling
3.61 0.00 0.00 # Lattice vectors
0.00 3.61 0.00
0.00 0.00 3.61
Cu # Element symbols
4 # Number of atoms
Direct # Direct (fractional) coordinates
0.00 0.00 0.00
0.50 0.50 0.00
0.50 0.00 0.50
0.00 0.50 0.50Key Points:
- Line 1: Comment (system description)
- Line 2: Universal scaling factor
- Lines 3-5: Lattice vectors (Å)
- Line 6: Element symbols (must match POTCAR order)
- Line 7: Number of atoms per element
- Line 8: Coordinate type (Direct or Cartesian)
- Lines 9+: Atomic positions
bash
Cu FCC bulk
1.0 # 通用缩放因子
3.61 0.00 0.00 # 晶格矢量（Å）
0.00 3.61 0.00
0.00 0.00 3.61
Cu # 元素符号（需与POTCAR顺序一致）
4 # 各元素原子数
Direct # 坐标类型（Direct为分数坐标，Cartesian为笛卡尔坐标）
0.00 0.00 0.00
0.50 0.50 0.00
0.50 0.00 0.50
0.00 0.50 0.50关键点：
- 第1行：注释（体系描述）
- 第2行：通用缩放因子
- 第3-5行：晶格矢量（单位：Å）
- 第6行：元素符号（必须与POTCAR顺序匹配）
- 第7行：各元素的原子数量
- 第8行：坐标类型（Direct或Cartesian）
- 第9行及以后：原子位置
KPOINTS: k-Point Sampling
KPOINTS：k点采样
Gamma-Centered Mesh (most common):
bash
Automatic mesh
0 # 0=automatic
Gamma # Gamma or Monkhorst-Pack
8 8 8 # k-point grid
0 0 0 # ShiftMonkhorst-Pack:
bash
Automatic mesh
0
Monkhorst-Pack
8 8 8
0 0 0Band Structure Path:
bash
k-points for band structure
10 # Number of points between high-symmetry points
Line-mode # Line mode for band structure
Reciprocal
0.0 0.0 0.0 !Γ
0.5 0.0 0.5 !X
0.5 0.0 0.5 !X
0.5 0.25 0.75 !WGamma中心网格（最常用）：
bash
Automatic mesh
0 # 0表示自动生成
Gamma # 可选Gamma或Monkhorst-Pack
8 8 8 # k点网格
0 0 0 # 偏移量Monkhorst-Pack网格：
bash
Automatic mesh
0
Monkhorst-Pack
8 8 8
0 0 0能带结构路径：
bash
k-points for band structure
10 # 高对称点之间的采样点数
Line-mode # 能带结构使用线模式
Reciprocal
0.0 0.0 0.0 !Γ
0.5 0.0 0.5 !X
0.5 0.0 0.5 !X
0.5 0.25 0.75 !WPOTCAR: Pseudopotentials
POTCAR：赝势
Generation:
bash
undefined生成方法：
bash
undefinedConcatenate POTCARs in same order as POSCAR
按POSCAR中的元素顺序拼接POTCAR
cat ~/vasp/potpaw_PBE/Cu/POTCAR > POTCAR
cat ~/vasp/potpaw_PBE/Cu/POTCAR > POTCAR
For compounds:
化合物的情况：
cat ~/vasp/potpaw_PBE/Cu/POTCAR
~/vasp/potpaw_PBE/O/POTCAR > POTCAR
~/vasp/potpaw_PBE/O/POTCAR > POTCAR
**Choosing POTCARs:**
- **Standard:** `potpaw_PBE/Element/POTCAR`
- **GW calculations:** `potpaw_PBE.52/Element/POTCAR` or `potpaw_PBE.54/`
- **_sv:** Include semicore states (more accurate, slower)
- **_pv:** Include p as valence
- **_h:** Harder potential (higher ENMAX)cat ~/vasp/potpaw_PBE/Cu/POTCAR
~/vasp/potpaw_PBE/O/POTCAR > POTCAR
~/vasp/potpaw_PBE/O/POTCAR > POTCAR
**POTCAR选择指南：**
- **标准计算：** `potpaw_PBE/Element/POTCAR`
- **GW计算：** `potpaw_PBE.52/Element/POTCAR` 或 `potpaw_PBE.54/`
- **_sv：** 包含半芯态（精度更高，速度更慢）
- **_pv：** 将p轨道视为价轨道
- **_h：** 更硬的势（ENMAX更高）Parameter Selection by Calculation Type
按计算类型选择参数
1. Structure Relaxation
1. 结构弛豫
INCAR:
bash
IBRION = 2 # Conjugate gradient
ISIF = 3 # Relax cell + ions
NSW = 100
EDIFFG = -0.02 # Force convergence
ISMEAR = 1 # Methfessel-Paxton
SIGMA = 0.2Convergence Criteria:
- : Force-based (recommended: -0.01 to -0.05 eV/Å)
EDIFFG < 0 - : Energy-based (less common)
EDIFFG > 0
INCAR设置：
bash
IBRION = 2 # 共轭梯度法
ISIF = 3 # 弛豫晶胞+离子
NSW = 100
EDIFFG = -0.02 # 力收敛阈值
ISMEAR = 1 # Methfessel-Paxton展宽
SIGMA = 0.2收敛判据：
- ：基于力的收敛（推荐值：-0.01至-0.05 eV/Å）
EDIFFG < 0 - ：基于能量的收敛（较少使用）
EDIFFG > 0
2. Static Calculation (Single-Point)
2. 静态计算（单点能）
INCAR:
bash
IBRION = -1 # No ionic updates
NSW = 0
ISMEAR = -5 # Tetrahedron (accurate DOS)INCAR设置：
bash
IBRION = -1 # 不更新离子位置
NSW = 0
ISMEAR = -5 # 四面体法（用于高精度态密度计算）OR
或者
ISMEAR = 0 # Gaussian (if tetra not converged)
SIGMA = 0.05
undefinedISMEAR = 0 # 高斯展宽（若四面体法不收敛时使用）
SIGMA = 0.05
undefined3. Band Structure
3. 能带结构
Step 1: Self-consistent calculation
bash
ICHARG = 2 # From atoms
LCHARG = .TRUE. # Write CHGCARStep 2: Non-self-consistent band structure
bash
ICHARG = 11 # Read CHGCAR, no update
LORBIT = 11 # Write PROCAR步骤1：自洽计算
bash
ICHARG = 2 # 从原子开始计算
LCHARG = .TRUE. # 写入CHGCAR步骤2：非自洽能带结构计算
bash
ICHARG = 11 # 读取CHGCAR，不更新电荷密度
LORBIT = 11 # 写入PROCARUse line-mode KPOINTS
使用线模式KPOINTS
undefinedundefined4. Density of States (DOS)
4. 态密度（DOS）
INCAR:
bash
ISMEAR = -5 # Tetrahedron method
LORBIT = 11 # Projected DOS
NEDOS = 3000 # DOS resolutionINCAR设置：
bash
ISMEAR = -5 # 四面体法
LORBIT = 11 # 投影态密度
NEDOS = 3000 # DOS分辨率Use dense k-point mesh
使用高密度k点网格
undefinedundefined5. Molecular Dynamics
5. 分子动力学
INCAR:
bash
IBRION = 0 # MD
NSW = 1000 # MD steps
POTIM = 1.0 # Time step (fs)
TEBEG = 300 # Start temperature (K)
TEEND = 300 # End temperature
SMASS = 0 # NVE: 0, NVT: >0
MDALGO = 2 # 1=Andersen, 2=Nose-HooverINCAR设置：
bash
IBRION = 0 # 分子动力学模式
NSW = 1000 # MD步数
POTIM = 1.0 # 时间步长（fs）
TEBEG = 300 # 起始温度（K）
TEEND = 300 # 结束温度
SMASS = 0 # NVE系综：0，NVT系综：>0
MDALGO = 2 # 1=Andersen热浴，2=Nose-Hoover热浴6. Phonons (DFPT)
6. 声子（DFPT）
INCAR:
bash
IBRION = 6 # DFPT for phonons
NFREE = 2 # Central differences
POTIM = 0.015 # Displacement (Å)
EDIFF = 1E-8 # Tight convergence!INCAR设置：
bash
IBRION = 6 # DFPT方法计算声子
NFREE = 2 # 中心差分
POTIM = 0.015 # 位移量（Å）
EDIFF = 1E-8 # 严格收敛！7. Elastic Constants
7. 弹性常数
INCAR:
bash
IBRION = 6 # DFPT
ISIF = 3
NFREE = 4 # For elastic constantsINCAR设置：
bash
IBRION = 6 # DFPT方法
ISIF = 3
NFREE = 4 # 用于计算弹性常数Advanced Parameters
进阶参数
Hybrid Functionals (HSE06, PBE0)
杂化泛函（HSE06、PBE0）
HSE06:
bash
LHFCALC = .TRUE. # Activate hybrid
HFSCREEN = 0.2 # HSE screening parameter
AEXX = 0.25 # Exact exchange fraction
ALGO = All # Or Damped
TIME = 0.4 # Damping for convergenceHSE06设置：
bash
LHFCALC = .TRUE. # 启用杂化泛函
HFSCREEN = 0.2 # HSE屏蔽参数
AEXX = 0.25 # 精确交换比例
ALGO = All # 或Damped
TIME = 0.4 # 收敛阻尼系数GW Calculations
GW计算
Step 1: DFT (PBE)
bash
ALGO = Exact
NBANDS = 200 # Many empty bands
LOPTICS = .TRUE.Step 2: GW
bash
ALGO = GW0 # Or EVGW
NOMEGA = 50步骤1：DFT（PBE）计算
bash
ALGO = Exact
NBANDS = 200 # 大量空轨道
LOPTICS = .TRUE.步骤2：GW计算
bash
ALGO = GW0 # 或EVGW
NOMEGA = 50DFT+U (Correlated Systems)
DFT+U（关联体系）
INCAR:
bash
LDAU = .TRUE.
LDAUTYPE = 2 # Dudarev
LDAUL = 2 -1 # l quantum number (d, s/p)
LDAUU = 5.0 0.0 # U value (eV)
LDAUJ = 0.0 0.0 # J valueINCAR设置：
bash
LDAU = .TRUE.
LDAUTYPE = 2 # Dudarev型
LDAUL = 2 -1 # l量子数（d轨道、s/p轨道）
LDAUU = 5.0 0.0 # U值（eV）
LDAUJ = 0.0 0.0 # J值van der Waals Corrections
范德华修正
DFT-D3:
bash
IVDW = 11 # DFT-D3 (Grimme)vdW-DF:
bash
GGA = MK # optPBE-vdW
LUSE_VDW = .TRUE.
AGGAC = 0.0000DFT-D3：
bash
IVDW = 11 # Grimme的DFT-D3方法vdW-DF：
bash
GGA = MK # optPBE-vdW
LUSE_VDW = .TRUE.
AGGAC = 0.0000Convergence Testing Strategy
收敛测试策略
1. k-Point Convergence
1. k点收敛测试
bash
undefinedbash
undefinedTest sequence
测试序列
KPOINTS: 4x4x4, 6x6x6, 8x8x8, 10x10x10, 12x12x12
KPOINTS: 4x4x4, 6x6x6, 8x8x8, 10x10x10, 12x12x12
Converged when ΔE < 1 meV/atom between successive grids
收敛标准：连续网格之间的能量差ΔE < 1 meV/原子
undefinedundefined2. Energy Cutoff Convergence
2. 能量截断收敛测试
bash
undefinedbash
undefinedTest ENCUT
测试ENCUT值
ENCUT: 400, 450, 500, 550, 600 eV
ENCUT: 400, 450, 500, 550, 600 eV
Converged when ΔE < 1 meV/atom
收敛标准：能量差ΔE < 1 meV/原子
Forces may need higher cutoff
计算力时可能需要更高的截断能
undefinedundefined3. Systematic Approach
3. 系统化方法
- First: Converge ENCUT (fix k-points at moderate density)
- Second: Converge k-points (use converged ENCUT)
- Document: Save convergence test results
- 第一步： 固定中等密度k点，收敛ENCUT
- 第二步： 使用收敛后的ENCUT，收敛k点
- 记录： 保存收敛测试结果
Smearing Methods (ISMEAR)
展宽方法（ISMEAR）
| ISMEAR | Method | Use Case |
|---|---|---|
| -5 | Tetrahedron | Static calcs, DOS, accurate energies |
| -4 | Tetrahedron+Blöchl | Like -5, slightly different |
| -1 | Fermi smearing | Metals |
| 0 | Gaussian | General purpose |
| 1+ | Methfessel-Paxton order N | Relaxations, metals |
Recommendations:
- Metals, relaxation: ISMEAR=1, SIGMA=0.2
- Semiconductors, relaxation: ISMEAR=0, SIGMA=0.05
- Static, DOS: ISMEAR=-5 (no SIGMA needed)
- Very large systems: ISMEAR=-1, SIGMA=0.1
| ISMEAR | 方法 | 使用场景 |
|---|---|---|
| -5 | 四面体法 | 静态计算、态密度、高精度能量计算 |
| -4 | 四面体+Blöchl法 | 类似-5，仅细节不同 |
| -1 | 费米展宽 | 金属体系 |
| 0 | 高斯展宽 | 通用场景 |
| 1+ | N阶Methfessel-Paxton展宽 | 弛豫计算、金属体系 |
推荐：
- 金属体系、弛豫计算： ISMEAR=1，SIGMA=0.2
- 半导体体系、弛豫计算： ISMEAR=0，SIGMA=0.05
- 静态计算、态密度： ISMEAR=-5（无需设置SIGMA）
- 超大体系： ISMEAR=-1，SIGMA=0.1
Common Parameter Combinations
常用参数组合
Standard Relaxation (Metals)
标准弛豫（金属体系）
bash
undefinedbash
undefinedINCAR
INCAR
ENCUT = 520
PREC = Accurate
IBRION = 2
ISIF = 3
NSW = 100
EDIFFG = -0.02
ISMEAR = 1
SIGMA = 0.2
ALGO = Fast
LREAL = Auto
ENCUT = 520
PREC = Accurate
IBRION = 2
ISIF = 3
NSW = 100
EDIFFG = -0.02
ISMEAR = 1
SIGMA = 0.2
ALGO = Fast
LREAL = Auto
KPOINTS
KPOINTS
Gamma-centered
0
Gamma
8 8 8
0 0 0
undefinedGamma-centered
0
Gamma
8 8 8
0 0 0
undefinedHigh-Accuracy Static Calculation
高精度静态计算
bash
undefinedbash
undefinedINCAR
INCAR
ENCUT = 600 # Higher cutoff
PREC = Accurate
IBRION = -1
NSW = 0
EDIFF = 1E-8 # Tight convergence
ISMEAR = -5 # Tetrahedron
ALGO = Normal
LREAL = .FALSE. # Reciprocal space
ENCUT = 600 # 更高的截断能
PREC = Accurate
IBRION = -1
NSW = 0
EDIFF = 1E-8 # 严格收敛阈值
ISMEAR = -5 # 四面体法
ALGO = Normal
LREAL = .FALSE. # 倒空间计算
KPOINTS (very dense)
KPOINTS（高密度）
0
Gamma
12 12 12
0 0 0
undefined0
Gamma
12 12 12
0 0 0
undefinedFast Testing Setup
快速测试设置
bash
undefinedbash
undefinedINCAR
INCAR
ENCUT = 400 # Lower cutoff
PREC = Normal
EDIFF = 1E-4 # Loose
ISMEAR = 0
SIGMA = 0.1
ALGO = Fast
LREAL = Auto
ENCUT = 400 # 较低的截断能
PREC = Normal
EDIFF = 1E-4 # 宽松收敛阈值
ISMEAR = 0
SIGMA = 0.1
ALGO = Fast
LREAL = Auto
KPOINTS (coarse)
KPOINTS（粗网格）
0
Gamma
4 4 4
0 0 0
undefined0
Gamma
4 4 4
0 0 0
undefinedPerformance Optimization
性能优化
Parallelization
并行化
INCAR:
bash
NCORE = 4 # Cores per band (orbital parallelization)INCAR设置：
bash
NCORE = 4 # 每个轨道使用的核心数（轨道并行化）OR
或者
NPAR = 8 # Number of groups for band parallelization
KPAR = 4 # k-point parallelization
LPLANE = .TRUE. # Plane-wise distribution
**Guidelines:**
- NCORE ≈ number of cores per node / 2-4
- KPAR = number of k-points (or divisor)
- For large systems (>100 atoms): NCORE=1-4
- For many k-points: Use KPARNPAR = 8 # 轨道并行化的分组数
KPAR = 4 # k点并行化
LPLANE = .TRUE. # 平面波分布方式
**指南：**
- NCORE ≈ 每个节点的核心数 / 2-4
- KPAR = k点数量（或其约数）
- 大体系（>100原子）：NCORE=1-4
- k点数量多时：使用KPARMemory Management
内存管理
bash
LREAL = Auto # Reduce memory for large systems
NCORE = 4 # Reduce memory per corebash
LREAL = Auto # 大体系下减少内存占用
NCORE = 4 # 降低每个核心的内存占用Error Handling
错误处理
Common Errors and Fixes
常见错误及修复方法
"ZBRENT: fatal error in bracketing"
bash
undefined"ZBRENT: fatal error in bracketing"
bash
undefinedFix: Reduce POTIM or use different IBRION
修复：减小POTIM或更换IBRION
POTIM = 0.2
**"EDDDAV: X eigenvalues not converged"**
```bashPOTIM = 0.2
**"EDDDAV: X eigenvalues not converged"**
```bashFix: Increase NELM, change ALGO
修复：增加NELM，更换ALGO
NELM = 200
ALGO = All
**"Sub-Space-Matrix is not hermitian"**
```bashNELM = 200
ALGO = All
**"Sub-Space-Matrix is not hermitian"**
```bashFix: Reduce POTIM, check structure
修复：减小POTIM，检查结构
POTIM = 0.1
SYMPREC = 1E-8
**SCF not converging:**
```bashPOTIM = 0.1
SYMPREC = 1E-8
**SCF不收敛：**
```bashTry sequential fixes:
尝试以下修复步骤：
- ALGO = All
- Increase NELM = 200
- AMIX = 0.2, BMIX = 0.0001
- Check initial structure (too close atoms?)
undefined- 设置ALGO = All
- 增加NELM = 200
- 设置AMIX = 0.2, BMIX = 0.0001
- 检查初始结构（原子间距是否过近？）
undefinedBest Practices
最佳实践
- Always Converge: Test k-points and ENCUT before production runs
- Use Symmetry: Let VASP detect symmetry (speeds up calculations)
- Check OUTCAR: Verify "reached required accuracy" message
- Monitor: Check OSZICAR during run for convergence
- Save Everything: Keep all outputs (OUTCAR, vasprun.xml) for analysis
- Consistent Pseudopotentials: Use same POTCAR set for all related calculations
- Document Settings: Record all INCAR parameters used
- 始终进行收敛测试： 正式计算前测试k点和ENCUT
- 利用对称性： 让VASP自动检测对称性（加快计算速度）
- 检查OUTCAR： 确认是否出现“reached required accuracy”提示
- 实时监控： 运行期间查看OSZICAR确认收敛情况
- 保存所有文件： 保留所有输出文件（OUTCAR、vasprun.xml）用于后续分析
- 使用一致的赝势： 相关计算使用同一套POTCAR
- 记录设置： 记录所有使用的INCAR参数
Calculation Workflows
计算工作流
Full Relaxation → Properties
完整弛豫→性质计算
- Relaxation: Optimize structure (ISIF=3, IBRION=2)
- Static: Accurate energy (ISMEAR=-5, dense k-points)
- Band Structure: Non-SCF with line-mode k-points
- DOS: Dense k-mesh with ISMEAR=-5
- Properties: Phonons, elastic, etc.
- 弛豫： 优化结构（ISIF=3，IBRION=2）
- 静态计算： 高精度能量计算（ISMEAR=-5，高密度k点）
- 能带结构： 使用线模式k点的非自洽计算
- 态密度： 高密度k点网格+ISMEAR=-5
- 性质计算： 声子、弹性常数等
Convergence Testing Workflow
收敛测试工作流
- Rough optimization: Low ENCUT, coarse k-points
- Test ENCUT: Fixed k-points, vary ENCUT
- Test k-points: Converged ENCUT, vary k-mesh
- Production: Use converged parameters
- 粗略优化： 低ENCUT、粗k点网格
- 测试ENCUT： 固定k点，调整ENCUT
- 测试k点： 使用收敛后的ENCUT，调整k点网格
- 正式计算： 使用收敛后的参数
Subskills
子技能
Invoke specific subskills for detailed guidance:
- relaxation - Structure optimization workflows
- electronic-structure - Band structure and DOS calculations
- molecular-dynamics - MD simulations in VASP
- advanced-functionals - Hybrid, GW, DFT+U methods
- phonons - DFPT phonon calculations
- convergence - Systematic convergence testing
调用以下特定子技能获取详细指导：
- relaxation - 结构优化工作流
- electronic-structure - 能带结构与态密度计算
- molecular-dynamics - VASP中的分子动力学模拟
- advanced-functionals - 杂化泛函、GW、DFT+U方法
- phonons - DFPT声子计算
- convergence - 系统化收敛测试
Quick Decision Guide
快速决策指南
What type of calculation?
| Goal | IBRION | ISIF | ISMEAR | EDIFFG |
|---|---|---|---|---|
| Relax ions only | 2 | 2 | 1 | -0.02 |
| Relax cell+ions | 2 | 3 | 1 | -0.02 |
| Static energy | -1 | 2 | -5 | N/A |
| MD simulation | 0 | 2 | 0 | N/A |
| Band structure | -1 | 2 | 0 | N/A |
| Phonons (DFPT) | 6 | 2 | 0 | N/A |
计算类型对应参数：
| 目标 | IBRION | ISIF | ISMEAR | EDIFFG |
|---|---|---|---|---|
| 仅弛豫离子 | 2 | 2 | 1 | -0.02 |
| 弛豫晶胞+离子 | 2 | 3 | 1 | -0.02 |
| 静态能量计算 | -1 | 2 | -5 | 无 |
| MD模拟 | 0 | 2 | 0 | 无 |
| 能带结构计算 | -1 | 2 | 0 | 无 |
| 声子计算（DFPT） | 6 | 2 | 0 | 无 |
References
参考资料
- VASP Manual: https://vasp.at/wiki/The_VASP_Manual
- VASP Tutorials: https://vasp.at/tutorials/latest/
- Parameter Index: https://vasp.at/wiki/index.php/Category:INCAR
See Also
相关内容
- skill - For ASE-based workflows with VASP
materials-properties - Examples in directory
examples/ - Detailed references in directory
references/
- 技能 - 基于ASE的VASP工作流
materials-properties - 目录中的示例
examples/ - 目录中的详细参考资料
references/