Known issues

Below we provide an incomplete list of known issues. Please mind the description to see whether the issue has been fixed.

Color legend: Open Resolved Planned Obsolete

Version fixed	Version first noticed	Date	Description
Open	<5.4.4	2024-11-14	There is a bug in the implementation of the van der Waals dDsC method (IVDW=4). The atomic volumes and charges are not calculated correctly, which translates into errors in the total energy and forces.
Open	<6	2024-08-26	NaNs if NSW*NBLOCK, NSW*ML_OUTBLOCK or NBLOCK*KBLOCK $> 2^{31}-1$: Either product must not exceed the largest integer(4) number 2147483647 ($=2^{31}-1$). As a solution split the molecular dynamics run into multiple calculations with smaller values for NSW or KBLOCK. Thanks to Renjie Chen who reported the behavior in this post: Error in MD annealing simulation with ML potential
Open	6.4.3	2024-08-26	Memory estimation in ML_MODE=TRAIN is wrong: The memory estimation for the major arrays (design matrix - FMAT, covariance matrix - CMAT, etc.) can be significantly too small. This problem mainly appears if ML_MODE=TRAIN is selected especially in a continuation run. Until this is not officially fixed don't rely on the memory estimation!
Open	<6	2024-08-21	Dielectric tensor and Born effective charges from density functional perturbation theory are incorrect for non-collinear spin calculations when symmetries are used: The rotation of the spinor part of the derivatives of the wavefunctions with respect to k was missing which leads to incorrect results when using LEPSILON=.TRUE. in combination with LNONCOLLINEAR=.TRUE. and ISYM>=0. The fix for previous versions is to use ISYM=-1.
Open	6.4.3	2024-08-05	ML_MODE=REFIT is broken if VASP is compiled without the precompiler flag -DscaLAPACK: Some arrays are allocated with the wrong size, leading to a crash with unclear exit messages of the code. This bug should not affect many users, since we strongly suggest to run ML_MODE=REFIT with -DscaLAPACK, because otherwise the SVD and related routines are not parallelized. This bug will be fixed in the next release.
Open	<6	2024-05-31	Compiler optimizations with the Fujitsu compiler on A64FX platforms: The bulk_BN_vdW-DF3-opt1 and TiO2_IBRION=2 tests might fail when VASP is compiled with distributed makefile.include.fujitsu_a64fx or makefile.include.fujitsu_a64fx_omp files. It is strongly suggested to modify the following line OBJECTS_O2 += fft3dlib.o nonl.o vdw_nl.o in the above mentioned files. Thanks to Ivan Rostov for the bug report.
Open	<6	2024-05-27	Calculations with LMODELHF=.TRUE. crash if started without WAVECAR file in the directory: The crash occurs because of a division by the screening parameter that is zero during the first few iterations that are done with the functional from the POTCAR file. If a WAVECAR file is present, then these first few iterations are skipped.
Open	6.4.3	2024-05-14	Using LCALCEPS in combination with Hybrid functionals may lead to a crash when running on GPU: VASP may crash when using Hybrid functionals in combination with LCALCEPS when running on GPUs due to an error when distributing the electronic states to be optimized for a batched fft. For the moment to overcome this issue, users can set NBLOCK_FOCK <= number of occupied states in the INCAR file. Thank you Sergey Lisenkov and Francesco Ricci for reporting the bug.
Open	<6	2024-05-13	Reading the file DYNMATFULL may lead to a crash in MPI-parallel calculations : If SCALEE≠1, then the file DYNMATFULL is read if present. This may lead to a crash in MPI-parallel calculations, in particular with the gfortran compiler. Thanks to Vyacheslav Bryantsev for the bug report.
Open	6.4.3	2024-04-10	Compilation error for GCC with ELPA support: Compilation with ELPA support (Makefile.include#ELPA_(optional)) failes for the GNU Fortran compiler because the Fortran standard for c_loc was not strictly followed. Other compilers (e.g. NVIDIAs Fortran compiler) might not enforce the standard in this case and will produce a working binary. Solution: Add the TARGET attribute to the variable declarations of matrices A in line 2236 and Z in line 2252 in src/scala.F. Thanks to user rogeli_grima for the bug report!
Open	6.4.2	2024-04-10	AOCC >= 4.0 does not produce runnable code when compiling without OpenMP-support: The AOCC compiler version >= 4.0 apparently uses a more aggressive optimization on a particular symmetry routine (SGRGEN) when compiling without OpenMP-support. Thus, code produced using arch/makefile.include.aocc_ompi_aocl exits with: VERY BAD NEWS! internal error in subroutine SGRGEN: Too many elements 49 Solution: adapt your makefile.include by adding symlib.o to the OBJECTS_O1 line. This fixes the issue. The other options is to compile with openMP support (using arch/makefile.include.aocc_ompi_aocl_omp). Thanks to users jelle_lagerweij, huangjs, and jun_yin2 for the bug report and investigations.
Open	6.4.3	2024-04-03	-DnoAugXCmeta is broken: We no longer recommend compilation of VASP with this precompiler option since it negatively affects the results of SCAN and SCAN-like meta-GGA calculations. To make matters worse, this feature is broken in VASP.6.4.3. So definitely do not compile VASP.6.4.3 with -DnoAugXCmeta.
Open	6.4.2	2024-03-21	Wannier90 exclude_bands not supported for SCDM method: When using LSCDM together with LWANNIER90 or LWANNIER90_RUN, the use of exclude_bands in the Wannier90 input file is currently not supported.
Planned	5.4.0	2024-10-14	Uninitialized variable IFLAG in ELMIN for ICHARG=5: When running VASP with ICHARG=5, the variable IFLAG is not properly initialized before calling EDDIAG. This caused inconsistent behavior during the EDDIAG call when using the DAV algorithm (each MPI rank has a random value for IFLAG). Depending on the compiler used, this could result in either VASP waiting indefinitely during EDDIAG due to waiting MPI ranks, or led to the absence of the requested preconditioning for the DAV solver, slowing down convergence. While this did not produce incorrect results, VASP could either hang indefinitely without throwing an error or exhibit slower convergence.
6.5.0	6.4.2	2024-11-13	SCDM method gives incorrect results for second spin channel when ISPIN = 2: The SCDM method (LSCDM = True) uses a rank-revealing QR decomposition with column pivoting to select the optimal columns of the density matrix. For the second spin channel, the pivot array is not correctly initialized to zero. This causes the results to be unreliable. Thanks to Patrick J. Taylor who reported the behavior in this forum post.
6.4.3	<6	2024-08-21	Interface to Wannier90 and PEAD calculations lead to incorrect results for non-collinear spin calculations when symmetries are used: The rotation of the spinor part of the wavefunctions was missing which leads to incorrect results when computing the projections and overlaps written to the AMN and MMN files used by Wannier90 when LNONCOLLINEAR=.TRUE. and ISYM>=0 are set in the INCAR file. The fix for previous versions is to use ISYM=-1.
6.4.3	6.4.2	2024-08-16	Hash codes in POSCAR and CONTCAR files: Hash codes are printed out to CONTCAR files. This does not affect the calculation but confuses some users. This has been fixed as of VASP 6.4.3, cf. forum posts: https://www.vasp.at/forum/viewtopic.php?f=4&t=19108 , https://www.vasp.at/forum/viewtopic.php?f=3&t=19113 , and https://vasp.at/forum/viewtopic.php?p=27108#p27108.
6.4.3	6.4.2	2024-02-06	The combination of VCAIMAGES and ISIF = 3 results in non-averaged trajectories: The issue arises because the stress tensor is not averaged over the two images of VASP runs.
6.4.3	6.2.1	2023-10-19	Phonon calculations (IBRION = 6) fails for some trigonal cells with ISIF = 3: VASP prints a bug error message complaining that it could not find some k points of the original mesh in the larger one with reduced symmetry of a distortion. As a fix, you can set KBLOWUP = F to circumvent this error message while we work on a fix. Thanks for barshab for the bug report.
6.4.3	6.4.2	2023-09-20	Specific cases of SAXIS gave unexpected quantization axis: For sx=0 and sy<0, VASP falsely assumes alpha=pi/2. It should correctly yield alpha=-pi/2. This error has probably been there for a long time but on one hand the setting is probably rarely chosen and on the other hand as the treatment is consistent within the calculation, the results should not be affected much.
6.4.3	6.4.2	2023-08-21	Restarting a calculation from vaspwave.h5 when the number of k points changed crashes with a bug message: This can happen, e.g., because ISYM is changed. VASP should behave the same as restarting from WAVECAR.
6.4.3	6.4.0	2023-04-06	LOCPOT file for vasp_ncl is not written correctly: LVTOT=T for vasp_ncl should write the potential in the "density, magnetization" representation, i.e., the scalar potential (v0), and magnetic field (Bx, By, Bz), to the LOCPOT file. However, VASP writes the potential in the (upup, updown, downup, downdown) representation to real numbers, which is incomplete.
6.4.2	6.4.0	2023-05-31	Fast-mode predictions will crash together with finite difference (IBRION=5,6): At the end of the calculation the fast-mode is supposed to deallocate important arrays using NSW. In the finite differences method NSW is not used and the fast-mode can wrongly deallocate at an earlier stage. This results in an error if the code wants to access the deallocated arrays. Until a patch is released we suggest two possible quick fixes: (1) Avoid explicit deallocations at the end of the program and let the compiler deallocate when the code runs out of scope. For that remove lines 568, 569, 570 and 572 in the ml_ff_ff2.F file. (2) Avoid the fast-prediction mode: Retrain the MLFF without support for the fast mode, i.e., use ML_MODE = refitbayesian instead of refit. Use the resulting ML_FFN file for the ML_MODE = run production runs. Unlike suggestion (1) this does not require a recompilation of VASP but unfortunately MD runs will take significantly longer. Thanks for Soungminbae for the bug report.
6.4.2	6.4.0	2023-05-17	Incorrect MLFF fast-mode predictions for some triclinic geometries: Due to an error in the cell list algorithm the MLFF predictions (energy, forces and stress tensor) in the fast-execution mode (ML_MODE = run) may be incorrect for triclinic systems with small or large lattice angles (i.e., large deviations from right angles). Until a patch is released we suggest two possible quick fixes: (1) Avoid using the cell list algorithm for neighbor list builds (recommended): Add this%algo_type = 2 in a new line below line 923 in src/ml_ff_neighbor.F and recompile VASP, or, (2) Avoid the fast-prediction mode: Retrain the MLFF without support for the fast mode, i.e., use ML_MODE = refitbayesian instead of refit. Use the resulting ML_FFN file for the ML_MODE = run production runs. Unlike suggestion (1) this does not require a recompilation of VASP but unfortunately MD runs will take significantly longer. Thanks to Johan for a very detailed bug report.
6.4.2	6.4.1	2023-05-15	Bugs in interface to wannier90: If no projections are supplied (e.g. LOCPROJ, LSCDM) and there are no projections found in the wannier90 input file, VASP does not produce the UNK files. This also leads to a crash if LWANNIER90 RUN is used. The format of the UNK files is broken for the gamma-only version of VASP. Thanks to guyohad for the bug report.
6.4.1	6.4.0	2023-03-07	Output of memory estimate in machine learning force fields is wrong for SVD refitting: The SVD algorithm (ML_IALGO_LINREG=3, 4) uses the design matrix and two helping arrays with the size of the design matrix. In the memory estimates these two helping arrays are not considered correctly. The entry "FMAT for basis" at the beginning of the ML_LOGFILE should be three times larger. The algorithm will be fixed such that it only requires twice the design matrix arrays instead of three times and the outputs for the estimates will contain the correct values.
6.4.1	6.4.0	2023-03-07	Bug in sparsification routine for machine learning force fields: This bug effects more severely calculatoins where the number of local reference configurations is getting close to ML_MB. By setting ML_MB to a high value this bug can be avoided in most cases (there are still some cases, especially where a small number of local reference configurations is picked and the structure contains many atoms per type or ML_MCONF_NEW is set to a high value). This bug can especially affect refitting runs, resulting in no ML_FFN file.
6.4.1	6.4.0	2023-03-07	ML_ISTART=2 on sub element types broken for fast force field: When the force is trained for multiple element types, but the production runs (ML_ISTART=2) are carried out for a subset of types, the code most likely crashes. This bug will be urgently fixed.
6.4.1	6.2.0	2023-02-20	INCAR reader issues: Moving an INCAR file from a system with Windows line endings to a Unix-based system can cause the INCAR reader to fail. As a workaround convert the INCAR file to Unix line endings e.g. by :set ff=unix in vi. Comments lines do not work properly with inline tags separated by semicolon if the comment character occurs before the semicolon but not at the beginning of the line. As a workaround, split the tags over multiple lines so that you can comment out what you want. Please also refer to the OUTCAR file and see whether VASP understood your input.
6.4.1	6.4.0	2023-02-17	Corrupt ML_FFN files on some file systems: Insufficient protection against concurrent write statements may lead to corrupt ML_FFN files on some file systems. The broken files will often remain unnoticed until they are used in a prediction-only run with ML_ISTART=2. Then, VASP is likely to exit with some misleading error message about incorrect types present in the ML_FF file. As a workaround it may help to refit starting from the last ML_AB file with ML_MODE=refit which may generate a working ML_FFN file (this is anyway highly recommended to gain access to the fast execution mode in ML_ISTART=2). Alternatively, there is a patch for VASP.6.4.0 available (see attachment to this forum post). Thanks a lot to xiliang_lian and szurlle for reporting this and testing the patch.
6.4.0	6.3.2	2023-01-18	makefile.include template does not work for AOCC 4.0.0: The flang preprocessor explicitly requires specifying that the code is in free format -ffree-form. In earlier versions of VASP you can add this flag to the CPP rule in the makefile.include. Thanks to liu_jiyuan for reporting this bug.
6.4.0	6.1.0	2022-11-23	Memory leak in MD in OpenMP version compiled with AOCC and NV: This problem originates from the DEFAULT(PRIVATE) clause in SET_DD_PAW in paw.F because the NV and AOCC compilers do not correctly clean up the memory for arrays that were allocated outside the OMP PARALLEL region and used as private inside. We advise against compiling with OpenMP support with the NV and AOCC compilers for vasp <= 6.3.2.
6.3.2	5.4.4	2021-11-12	Ionic contributions to the macroscopic polarization with atoms at the periodic boundary: Removed a section of code from POINT_CHARGE_DIPOL that adds a copy of the atom when it is at the periodic boundary. This can lead to a different value of "Ionic dipole moment: p[ion]" being reported in the OUTCAR with respect to previous versions of VASP. This result, although numerically different is still correct since the polarization is defined up to integer multiples of the polarization quantum. Thanks to Chengcheng Xiao the bug report.
6.3.2	6.3.1	2022-05-11	ML_ISTART=1 fails for some scenarios: Due to a bug in the rearrangement of the structures found on the ML_AB file, restarting the training of a force field by means of ML_ISTART=1 fails in some cases. N.B.: this problem only occurs in a scenario where one repeatedly restarts the training, and returns to training for a structure that was trained on before (that means exactly same element types and number of atoms per element), but not immediately before. Example: one starts training a force field for structure A, follows this by a continuation run to train for structure B, and then restarts a second time returning to training for structure A again.
6.3.1	6.2.0	2022-05-05	Treatment of the Coulomb divergence in hybrid-functional band-structure calculations is only correct for PBE0: The Coulomb divergence correction for states at and near the Γ-point in hybrid-functional band-structure calculations (see HFRCUT) was only correctly implemented for PBE0 and HFRCUT=-1. Note: HSE band-structure calculations are not expected to be (strongly) affected because this hybrid functional only includes “short-range” Fock exchange.
6.3.1	6.2.0	2022-03-14	Bug in interface with Wannier90 for non-collinear spin calculations: The spin axis for non-collinear spin calculations is not correctly read from the wannier90 input file. This is because this line in the mlwf.F file: MLWF%LPRJ_functions(IS)%spin_qaxis = proj_s_qaxisx(3,IS) should instead be: MLWF%LPRJ_functions(IS)%spin_qaxis = proj_s_qaxisx(:,IS). Thanks to Domenico Di Sante for reporting this bug.
6.3.1	6.3.0	2022-02-04	Incompatibility with Fujitsu compiler: Fujitsu's Fortran compiler does not support overloaded internal subroutines. A simple workaround is to compile without machine-learning–force-fileds capabilities. Comment out the macro definition of ML_AVAILABLE in line 626 of src/symbol.inc by adding a ! in front, i.e. it should look like this: !#define ML_AVAILABLE. Then do a complete rebuild of VASP: run make veryclean followed by your desired build command.
6.3.0	6.2.0	2021-05-28	Bug in interface with Wannier90 writing UNK when exclude_bands present: The UNK files generated by VASP include all bands where bands specified by `exclude_bands` should be excluded. The fix is to pass the `exclude_bands` array to `get_wave_functions` in mlwf.F. Thanks to Chengcheng Xiao for reporting this bug.
6.2.0	6.1.0	2022-08-29	Inconsistent energy for fixed electron occupancies: Rickard Armiento pointed out that the HF total energy for fixed electron occupancies was inconsistent when compared to 5.4.4 or older versions. This bug was introduced in 6.1.0 in order to support IALGO=3 in combination with ISMEAR=-2 (for SPHPRO calculations as post-processing step) but broke the CG algorithms (IALGO=53) The fix was added in src/main.F with IF (INFO%LONESW .OR. (INFO%IALGO==3 .AND. KPOINTS%ISMEAR/=-2)) THEN \n IF (INFO%LONESW) W_F%CELTOT = W%CELTOT .
>=6	<6	2023-10-31	For LORBIT >= 11 and ISYM = 2, the partial charge densities are not correctly symmetrized: This can result in different charges for symmetrically equivalent partial charge densities. For older versions of VASP, we recommend a two-step procedure: 1. Self-consistent calculation with symmetry switched on (ISYM=2) 2. Recalculation of the partial charge density with symmetry switched off (ISYM=0) To avoid unnecessary large WAVECAR files, we recommend setting LWAVE=.FALSE. in step 2.
PBE.64	-	2024-10-25	Date in Nd POTCAR of release PBE.64 should be 25 May 2022 (not 2002): There is a typo in the first line of the Nd POTCAR for the PBE version 64 release. It reads "PAW_PBE Nd 25May2002" instead of "PAW_PBE Nd 25May2022".

---