1, Restart the calculation with a pre-converged PBE WAVECAR (ISTART = 1, the default value when WAVECAR exists) 2, Decrease the mixing parameter AMIX from the default 0.4 to 0.2 and to 0.1, increase NELM to 12 0 or more. 3, Try the Anderson mixing by setting IMIX = 1, AMIX = 0.1, BMIX = 0.0001, AMIX_MAG = 2*AMIX, BMIX_MAG = 0.0001 4, The above methods are usually enough, and I would pray for a while if they are not ... 来自: https://sites.google.com/site/haoweipeng/softwares https://sites.google.com/site/haoweipeng/softwares
上一篇探讨了实数收敛的概念,用无穷级数来确定一个数。很容易把它扩展到无穷的函数序列求极限的问题。在初等微积分里,这是对函数变量的每个值逐点来考察,对每个固定的变量值,这无穷序列对应着一个函数值的数列,如果所有点对应的数列都收敛,那就认为这无穷函数序列收敛,它的极限函数在每个点的函数值是相应数列的极限值。这样函数序列的极限称为逐点收敛的。这个方法被牛顿引入后,广泛地应用,它虽然可行,但在微积分进一步研究时又遇到种种麻烦,于是又附加了许多条件,如“一致连续”,“绝对可积”等等,最后弄得微积分繁杂不堪。能不能把整个函数看成一个数学空间里的一个点,把这些条件都看成空间里的性质,从一个统一的角度来研究收敛极限的问题?这便是这一篇要介绍的概念。 现代的数学建立在比实数更加抽象的集合论基础上,应用于更广泛的空间。要将定义在实数上一元函数微积分的本质说清楚,推广到多元函数,函数逼近,泛函,随机过程,乃至各种抽象数学结构的集合上,我们要了解集合元素间联系的结构,这样才可能描述变动个体的走向,空间的性质,进而谈及趋近、收敛和极限。 微积分是基于无穷逼近极限的数学。收敛描述的是变动差别越来越小,直至微不可察的数列表现。收敛极限的存在,取决于实数的完备性。所以无穷逼近过程的含义和结果,依赖于它所在数学空间的性质。传统微积分是建立在实数 空间 $\mathbb{R}$ 的性质上。这些实数性质是人们在习以为常的计算经验中总结“发现”出来,又用一些更简单基本的原理来证明的。 在现代数学之前,数学一直被认为是研究真理的学问,直到百多年前,哲人们还都认为不但如此,而且是绝对的真理。现在看来,这些真理也许不过是:习惯成了自然的经验在逻辑上自洽的表现。几千年实践认知形成有穷数学的直观与那时的理论吻合无间。无穷则无法验证,除了凌空而来的假定和逻辑自洽之外无所凭借,一切概念和引为依据的原理都需要精确的定义引入,一切的争执最终也只能根据这些定义和假设,用逻辑作出判决。直觉、经验、试验和公认之事,对数学都不足为凭。数学的价值在于有用,在虚构的空间里,以严谨的逻辑铸造出一把犀利的思想宝刀,当你用对了,无数的实践难题可以迎刃而解。 要想自如地在在抽象的无穷空间中行走,你需要正确的直观想象。前面两篇只是个导引,让你意识到必须纠正过去的直觉,才能放开心怀接纳无穷空间里的真实。在“无穷”篇里数学归纳法的例子,证明了在有限情况都成立的论断,在无穷时未必成立。要走出有限的误区,才能走出局限。在“收敛”篇里 0.999…=1 的分析,说明了极限的含义。这是微积分数学模型的最基本假设,无穷过程达到彼岸的立足点。所有收敛的极限都是依相同的原理,将差异的不等式在极限处变成了等式。除非你曾经震惊思考改变过,如果初读前面这两个例子,没有经历过困惑,思考到终于理解了,你大约还没有走出有穷世界的藩篱。建议你回去走通了再来,否则学过的微积分知识只是别人说的,不是自己会的,此后三篇让你构筑无穷空间想象的介绍,就可能格格难下不能吸收,无缘看见彼岸。 直觉是模糊的联想,并非精确可靠。别人说的风景和你所想的也许不是同一个地方。只有踏在共同之处才是 真实的 一样。想象必须依附于证实过的判断作为立足点,只有在坚实的桩脚附近延伸,才靠近真实,过度延伸的浮想只是梦呓。可靠的直觉来自大量被证实过的事例,它们编织成了图像,越多的事例,才有越清晰的直觉。在现实有穷的世界,这些事例来自日积月累的学习和实践,判决正误唯有实验。在抽象无穷的空间,可靠的事例来自符合定义的例子和严格证明过的判断,在这里逻辑是唯一的判据。你必须记忆概念的定义,消化例子,纠正限于有穷世界的直觉,才会形成在无穷空间里的正确直观。 这一篇黑体字强调的概念,如开集、闭集、邻域、收敛、稠集、连续等等,你可能已经认识,不过在初等微积分学的是在实数上的定义,这儿是从抽象集合的高处俯视。这篇把点集拓扑的基本概念串起来介绍,有足够的信息让你可以依想象,走通从定义到例子的逻辑。希望你从过去的认知起步,以这里的定义来消除局限,用例子来校正偏差,在走通的过程中重建图像。花费时间逐个概念走过,你将会开始有看清无穷空间景象的基础。 好吧,抽象!讨论问题的空间一切属性都需要定义。忘掉实数、无穷、收敛、极限种种课堂灌给我们前人发现的真理,从简单干净的抽象集合开始。上一篇说,无穷序列作为数学模型来描述收敛的观念,序列中的项以及极限之间的“相像”性质,决定了收敛和极限的含义。好的。对集合中的两个元素,比如说是两个人,你能从其一了解另一吗?不能。因为元素只是集合中抽象的个体,没有它们间关系的信息。如果说知道它们同在某一个子集里,好比说知道两人同是一族,你就能知道,他们比不在这子集里的外族人有某些更相像。在这里用子集表达其中的元素具有某些共同的属性,这样的子集称之为开集,它们的并和交满足某些封闭性。集合上有了这样的一组开集,便能以此判断元素间属性是否相近。可以用来判别相邻性的这组开集,称为集合上的拓扑,它们组成了讨论无穷逼近问题的数学空间。 考虑集合X上的一个子集族τ,其元素是X的子集,如果它们任意多的并,及有限个的交都属于τ,空集和X也属于τ,那么τ中元素称为 开集 ,τ称为X上的 拓扑 (Topology),组合(X,τ)称为 拓扑空间 ,在简略τ不致混淆时,也可直接称X为拓扑空间。 在任何拓扑,空集和 X 都是开集,只具有这两个开集的拓扑,称为 平凡拓扑 。取 X 上的一组子集,包括空集和 X 本身,定义它们为开集,进而将它们任意多的并,及有限的交,都视为这空间的开集,则生成了一个拓扑τ。给几个实例来帮助想象。 例3.1:定义实数上的开区间为开集,在它们间进行任意多的并和有限个交运算,依定义这生成了实数上的拓扑。通常实数空间指的就是以此为拓扑的空间。 例3.2:在闭区间 上所有连续函数的集合C ,对任意的$f\in C ,\;\;\delta 0$,定义集合 $\{ g \in C \; |\; \max_{a\leq t \leq b} |g(t)-f(t)|\delta\}$ 为开集,它们任意多的并和有限个交,生成了C 空间上的一个拓扑。 例3.3:集合X={a,b},定义空集Φ,单点集{a}和X为开集,它们构成X上的一个拓扑。 例3.4:将X上所有子集都定义为开集的拓扑,称为 离散拓扑 。它也可以由定义任何单点集都为开集来生成的。 拓扑空间( X ,τ)中集合 X 称为空间,它的子集简称为集合, X 的元素称为点,点之间的相邻关系,由τ来确定。两个点,如果同在拓扑空间的一个开集里,认为它们是相邻,分属两个不相交的开集,认为是分离;包含了点 x 的开集,称为 x 的 开邻域 ,包住 x 开邻域的集合称为 x 的 邻域 ( neighborhood ),寓意是与这点某些属性相近的点都在这里。 那么,怎么用拓扑来描述趋近状态呢? 对于拓扑空间 X 的一个点 x 和序列( x n ),如果对于 x 的每个邻域 U ,都有对应着它的一个正数 N ,当 n N 时,都有 $x_n \in U$ ,则称( x n ) 收敛 于 x 。考虑到点的邻域含义,这意味着,对这空间中拓扑所考虑的属性上,这序列后面的点终将与点 x 是密不可分的。 对于拓扑空间 X 中一个点 x 和集合 A ,如果 x 的每个邻域都包含有除了 x 外,还有 A 中其他的点,则 x 称为 A 的一个 聚点 ( cluster point )。例如:在实数空间中, 0 和 1 是开区间( 0 , 1 )的聚点;任何一个无理数都是有理数集的聚点。 显然,序列收敛到非序列中的点 x ,是序列所在集合的聚点;但集合 A 的聚点,未必有集合中的序列收敛于它。如果集合中的点都有一组可数的邻域,这点所有邻域都包含这组的某些成员,这样的拓扑空间称为 第一可数的 ( first countable )。 第一可数空间集合中的聚点都有集合中的序列收敛于它 。在分析应用中的聚点都是用序列来逼近的,所以应用中的拓扑空间都是第一可数的,它也称为可数性的第一公理。 开集的补集称为 闭集 ;包住集合 A 的所有闭集的交,是包住 A 的最小闭集,叫做 A 的 闭包 。闭包也是包括了 A 和它所有聚点的集合。实数是有理数的闭包。闭区间 是开区间 (a,b) 的闭包,闭区间 是集合 [1, 2)U(2, 3) 的闭包。 对拓扑空间 X 的子集 A ,如果任何开集都包含有 A 中的点,则称 A 是 X 中的 稠集 ,意味着 X 中的点的任何邻域,都含有 A 中的点,即可以用 A 中的点来逼近。如果 X 有个稠集是可数的,则称这空间是 可分的 , X 中的点都有这稠集中的一个序列收敛于它。有理数是可数的,它是实数上的稠集,所以实数空间是可分的,每个实数都有一个有理数的数列收敛于它。所有实数都是有理数集合的聚点,实数空间的拓扑是第一可数的。 同一个集合 X 上可以定义不同的拓扑。平凡拓扑是最粗的拓扑,离散拓扑是最细的。开集设定了一个点的邻域,它意味着里面其它点与这点相类。邻域之外就有某些的“不相类”,所以空间的拓扑区分的能力也很重要。如果对两个点,总有一个点开邻域不含另一,称为 T 0 空间;它们各自都有不包含对方的开邻域,是 T 1 空间;如果两个点能有分属不相交的开邻域,叫 T 2 空间 ,这也称 Hausdorff 空间 ; T 3 , T 4 则是说明点与闭集的区分能力。后者包括前者。过粗的拓扑,让空间中的点含糊难分。过细的拓扑也孤立难同,则难以找到点邻域里的其它点来推测它的性质。应用中大多数具有丰富性质的拓扑,既不会太粗也不会太细。现代分析是建立在点集拓扑研究的基础上。这样可以从简单拓扑结构开始,层层加细,从不断丰富拓扑空间性质的研究中,探究各种特性的本质。 对于微小的变动过程,如果集合上的函数,对应的值的变化也是微小的,即 $x$ 趋向 $a$ 时, $f(x)$ 趋向 $f(a)$ ,这是一种叫做“连续”的对应关系。用拓扑的语言来描述函数 $f(\cdot)$ 在点 $a$ 处连续,它的定义是函数值 $f(a)$ 的每个邻域,都包含了 $a$ 一个邻域的映像(即这邻域所有函数值的集合)。说函数 $f(\cdot)$ 是 连续的 ,指它对定义域中所有的点都是连续的。 从拓扑空间 X 到 Y 的连续函数意味着, Y 的每个开集 U ,它的原像 f -1 (U) 在 X 中也是开集;闭集 V 的原像 f -1 (V) 也是闭集。连续函数把它对应的两个空间的拓扑结构连系起来了。 两个拓扑空间之间如果有个一一满映射 F , F 和 F -1 都是连续的,则称 F 是同胚映射,这两个拓扑空间称为是 同胚 的( Homeomorphism )。想象拓扑空间 X 就像一个有弹性的几何体,同胚映射 F 就好比将这几何体拉伸压缩变形,形成另外一个几何体拓扑空间 Y 。连续函数 F 在这变形中,让所有点的相邻关系都没变。如果 F 的映射不是满的,则说 F 将 X 嵌入 ( embedding )到 Y 空间中。传统的拓扑学,用同胚映射和同伦映射对空间进行几何分类和代数特征的研究。分析应用上,将研究的拓扑空间同胚映射或嵌入到熟悉的空间中,可以比较直观地了解它的性质。 上面所说的是在集合的基础上,抽象地定义具有收敛概念的数学空间。这里尽量用简洁清晰的语言表达概念的定义,如果仍有未解之处,建议查看点集拓扑的教科书来加深理解。你的努力,将给了你一双能够看到无穷空间的眼睛。 现在检验你对抽象定义的理解力,请思考下面的概念题。 对拓扑空间 X 中的集合 M , X 的所有开集与 M 的交集定义为 M 中的开集,构成 M 的拓扑,称 M 为 X 的 拓扑子空间 。 M 是个拓扑空间,它继承了 X 的拓扑。例如实数 $\mathbb{R}$ 是 3 维欧几里德空间 $\mathbb{R}^3$ 的拓扑子空间。 开区间( 0 , 1 )继承实数空间的拓扑,问在这拓扑子空间里( 0 , 1 )是闭集吗?它还是开集吗?在这里的有界数列都收敛吗?这子空间是完备的吗?这拓扑子空间与实数空间是同胚的吗? (待续) 【扩展阅读】 维基百科,拓扑空间 http://zh.wikipedia.org/wiki/%E6%8B%93%E6%89%91%E7%A9%BA%E9%97%B4 Stephen Willard ,General Topology, Addison-Wesley(1970) Renzo’s Math 490 ,Introduction to Topology http://www.math.colostate.edu/~renzo/teaching/Topology10/Notes.pdf
对于一些磁性体系、镧系和锕系元素及相关化合物的静态计算(电子迭代),经常会遇到“难收敛”的问题。 下面给出几个相关Flag及设置方法:(有些话从说明书上摘录的,不大好翻译,就copy了,请见谅) 1、LMAXMIX Default: LMAXMIX = 2 An additional flag controls up to which l quantum number the onsite PAW charge densities are passed through the charge density mixer. Higher l-quantum numbers are usually not handled by the mixer. In order to obtain fast convergence to the groundstate, you can try the following setting: LMAXMIX = 4 for d elements LMAXMIX = 6 for f elements 这个FLAG对于含d电子和f电子的体系是非常重要的,很大一部分体系的收敛问题可以通过设置合适的LMAXMIX值来解决。 2、ALGO, IALGO, LDIAG If the self-consistency loop does not converge within 40 steps, it will probably not converge at all. In this case you should reconsider the tags IALGO, LDIAG, and the mixing-parameters. 这是说明书上的建议。 一般情况下,或使用IALGO=48时遇到收敛问题的话,可以考虑设IALGO为38(4.5以前的版本可设为8),或设置ALGO=Normal or Fast (in VASP.4.5 and later versions)。 3、NELMDL NELMDL gives the number of non-selfconsistent steps at the beginning; if one initializes the wave functions randomly the initial wave functions are far from anything reasonable. The resulting charge density is also 'nonsense'. Therefore it makes sense to keep the initial Hamiltonian, which corresponds to the superposition of atomic charge densities, fixed during the first few steps. Choosing a 'delay' for starting the charge density update becomes essential in all cases where the SC-convergence is very bad (e.g. surfaces or molecules/clusters chains). Without setting a delay VASP will probably not converge or at least the convergence speed is slowed down. NELMDL might be positive or negative. A positive number means that a delay is applied after each ionic movement — in general not a convenient option. A negative value results in a delay only for the start-configuration. 4、mixing-parameters 对于一些难收敛的体系,可以使用“linear mixing”,具体详见VASP说明书中的“Mixing-tags”。 For an initial linear mixing (BMIX ~ 0) an optimal setting for A(AMIX) can be found easily by setting Aopt=Acurrent*Γmean. For the Kerker scheme either A or q0(i.e. AMIX or BMIX) can be optimized, but we recommend to change only BMIX and keep AMIX fixed (you must decrease BMIX if the mean eigenvalue is larger than one, and increase BMIX if the mean eigenvalue is smaller than one). 尽管VASP说明书中给出了调节AMIX和BMIX的一些较为明确的建议,但是实际去调节的时候,还是挺难的,但原则上说,是可以通过调节这两个Flag来使得收敛问题得以解决的,只是得有耐心。 5、kmesh, SIGMA 收敛问题还跟kmesh及SIGMA(当使用ISMEAR不等于-5 和-4时)的设置有关。要达到同样的精度,较小的SIGMA则需要较大的kmesh;而且,当SIGMA较小时,若kpoints不够多,也会出现难收敛的情况。 P.S. 对于不同系统的计算,问题的原因不一定一样,因而可能解决之一问题的方法也不一定会一样。其实以上大部分都来自VASP的说明书,只是我在遇到收敛问题时试过这些方法. 写得比较仓促,以后再修改及补充。 希望有这方面经验的朋友也把自己的见解提出来一起讨论。 金属表面体系。换算法,改Mixing都不管用。 收敛和k-mesh有关系吗?谢谢了 别改Mixing。 你把IBRION设为3,然后把SMASS设为0.6试试看 不好意思,我没说清楚 IBRION是对应离子步的算法 我的问题是电子步不收敛 dE 一直0.1左右反复, 200步也不收敛 改IALGO为38, 增加NBANDS都不行:( 那你试试不要采用自恰计算。但是仍然采用damping驰豫。 是不是你的初始构型离平衡态太远了。 还没遇到过这种问题,可能比较运气吧。 调一下cutoff试试。 IALGO我一直用48,5x5x1 Kpoints。我做金属的表面体系都是在100步内就OK了。 如果愿意,可以把INCAR贴出来,让大家出出注意,这样更可能有的放矢。 应该不会 有的时候甚至弛豫出来的结构再进行static的计算都不收敛 收敛跟k-mesh是有一定关系的。 你INCAR中没有设置ISMEAR,应该就是使用默认的ISMEAR=1。此时SIGMA应该也是默认值0.2。 你看看OUTCAR中的EENTRO的值,看是不是符合要求(一般要求小于1 meV/atom,这与你想要的精度有关,如果你要求精度不高的话,这个标准可以设大一些)。 有几个方向你可以试试,或许会有用: 1、如果EENTRO的值够小,你可以尝试加大SIGMA。(SIGMA较大时,相对收敛会快些) 2、增加k点多。 其实就是说,要调整k-mesh和SIGMA。一般SIGMA跟k点数会有一定的关系。 若想达到同一精度,SIGMA越小所需要的K点数就越多。如果SIGMA较小,而k点数不足,有可能会出现难收敛的情况,也可能得到不准确的结果。 3、调节AMIX与BMIX 我之前也遇到过难收敛的问题,当时调这两个参数,还是有点用的。具体调节的方法VASP上有讲到一些,但可能还是要多试才知道。我现在一时想不起来细节了,你先试试,如果有情况,再一起讨论。 VASP的说明书也有提到,如果电子迭代超过40次仍未收敛的话,一般就很难收敛了。当然,我想这只是VASP说明书的一个建议,并不是绝对的。一般设置NELM=60(默认)就够了,如果60次都不能收敛,那可能就要找找问题或是改相应的设置了。 有不对或不妥的地方,大家帮忙指正:) 从你的INCAR来看,貌似没什么问题。如果还是不收敛的话,你可以增加一点截断能,比如400 eV,并把POTIM减小至0.1。另外,添加三个tag:NELM=100; NELMDL=6; WEIMIN=0. 大家遇到收敛问题怎么解决啊??? 小弟仔细读过使用手册,并结合一些帖子,可是还是解决不了问题。恳请高手指点。 以下是我所用的方法和遇到的问题: 解决方法一:设置scf参数 a1. Iterations. 加大循环次数。计算过程中默认50次,我的体系无论加多大次数就是不收敛。 可能这不是主要原因。我现在一般设置一个中等大小的值400。 a2. mixing.对强阻尼体系,降低mixing值到0.01-0.05.甚至用mixing=0.001,diis n=0.保存TAPE21文件, 再用mixing=0.03,diis n=10.我都用了,都不收敛。分析结果计算过程中有能量跳跃,因此我分析我的体系就 是较强的阻尼,因此我设置mixing=0.03,这个因该可以吧? a3. lshift, vshift.手册上未解释怎么用。 这个我不懂。请问在什么情况下用?怎么用?大家一般设置值为多少? 解决方法二:用occupations指定电子占据。 b1. 直接指定分子占据,由于小弟的量化基础薄弱,弄不懂怎么指定。 我是按论坛上的方法之一,先算一个单点能,然后得到所有不可约表示的电子占据数。但是, 我计算单点能虽然scf收敛了,但是计算得到的能量明显不对,估计是收敛到激发态上去了。请问怎么指定到基态? 另外,对电子占据大家都有什么具体的着?哪位有这方面的文献和资料?我算单点的输出文件在附件中,请高手帮我 指定一下基态的电子占据。先谢谢了!!! 这是最重要,最有用的方法,可是我不懂,郁闷啊!!! b2. 用keeporbitals或freeze控制能量跳跃,但是我的体系用这个关键词还是无法得到几何结构收敛。我只用过 keeporbitals=20,10,5. 但是scf有几个收敛,可是最后得不到几何结构收敛。 大家一般设置值为多少?freeze又该如何设置? b3. smearq. 手册上说先用smearq=一个值。几何收敛后再用smearq=0, restart 前面的TAPE21文件。 我设置 smearq=0.2,几何结构收敛了,但是我用smearq=0,restart不收敛。我用smearq=0.1计算连几何结构 都不收敛。请问大家设置smearq=多少合适?用smearq且令smearq=0之后的结果可信度有多大啊?我看过几篇文献, 但是都没有什么有用的东西。 我已收到你的邮件了,个人觉得在这么多方法中用smearq 再 occupations比较有效: 先用启用smearing (例如:smearq=0.02), 再根据结果文件指认电子占据一般都能解决问题, 例如: occupations A1 a//b B1 a//b E a//b . . . end. Dear Vasp Community: I'm trying to do some very basic defect calculations in a class of complex oxides with the stoichiometry AB6O12, where A=U,W,Re and B=Y,Yb,Lu. In most cases, especially when A=W, I haven't had any problems. But, I have had problems with convergence in two cases that I haven't been able to make any progress on. When A=U and B=Y, I can converge the perfect structure very nicely. But, when I look at an antisite defect, by replacing one Y with U, the calculation doesn't converge. The electronic iterations get stuck, changing by +/- the same amount. A number of colleagues have suggested I change the smearing, set ISPIN=2, set LMIXMAX=6, etc. Everything I've tried, either doesn't change the behavior or does fix the electronic iterations, but after 3 or 4 ionic iterations, the forces on the ions blow up (become unphyiscally large, like 100 eV/angstrom). The second problem, when A=Re and B=Lu, is that I don't get any kind of convergence, even in the perfect case. I get a warning at the top VERY BAD NEWS! internal error in subroutine IBZKPT: Reciprocal lattice and k-lattice belong to different class of lattices. 168 But all I've done compared to the A=U, B=Y case is change the POTCAR file. If I let it go, I get similar convergence problems as before. Any suggestions on what I might try to fix these problems would be greatly appreciated. I'm using ISYM=0 and a Monkhorst-Pack grid of 2x2x2 for these calculations, along with the GGA PAW pseudopotentials. Thank you, Blas Back to top admin Tue Jan 08 2008, 01:27PM posts 1815 1) for the electronic convergence --) LMAXMIX=6 should be set if you have rare-earth elements in the cell, for d-elements, please set it to 4 --) it may help to decrease the mixing parameters (AMIX, BMIX, and in addition AMIX_MAX and BMIX_MAX for the spin-polarized runs to improve convergence ) --) also, it may help to increase NELMDL (the number of non-selconsistent steps at the beginning) to improve the pre-convergence of the wavefunction --) maybe it also helps to use L(S)DA+U ( please check the electronic structure and compare to experiment, if results are available) 2) for the forces which suddenly increase: --) please check (in OSZICAR) if each ionic step was fully converged electronically up to the requested convergence limit (EDIFF) before the forces are calculated. (i.e. whether the forces were calculated after el. convergence is reached or because the number of electronic convergence steps (NELM=60, by default) was reached without actually being converged. The latter yields unreasonable forces. (in that case, please either increase NELM and/or try one of the above mentioned to improve the electronic convergence) --) if each of the ionic steps was converged, please check the relaxation history (XDATCAR file). Sometimes, especially the first ionic steps 'overshoot' a little. If that is the case, please choose a different relaxation algorithm (IBRION) and/or decrease the step width (POTIM) 3) the VERY BAD NEWS message is due to different lattice types obtained for the unit cell and the symmetrized k-mesh. (triclinic and simple monoclinic as it seems). maybe the axes ratios of the unit cells are not 1:1:1 ? please check. for the symmetry analysis, the Bravais Matrix , atoms' positions, and, if set, magnetic moments and atoms' velocities are taken into account. If nothing has been changed except the atoms in the PP, it seems strange that this warning shows up for some calculations whereas it does not for others. In any case, if ISYM=0 is set, symmetrization is not taken into account in the calculation later on, the only point to be careful with is the axis' ratio then. LDAU = .TRUE. LDAUL = LDAUU = LDAUJ =