Dear Wien2K Users,<div> </div><div> I want to use virtual crystal approximation(VCA) for a very large system, so I applied it on a simple system [Ni(28)] to check whether everything is correct or not. To do that, I did the following steps.</div>
<div> 1. Initialized the calculation with Ni(Z=28).</div><div> 2. Changed Z to 28.15</div><div> 3. Changed case.inst to:</div><div>--------------------------------------------------------------</div><div><div>Ni</div><div>
Ar 3</div><div>3, 2,2.00 N</div><div>3, 2,2.00 N</div><div>3,-3,3.00 N</div><div>3,-3,1.30 N</div><div>4,-1,1.00 N</div><div>4,-1,0.85 N</div><div>****</div><div>****</div><div>---------------------------------------------------------------</div>
<div> 4. Changed NE in case.in2* from 16 to 16.15</div><div> 5. Run the spin polarized calculation.</div><div> I also did a supercell calculation where the 2x2x2 supercell was containing 32 non equivalent Ni atoms. Then, I replaced 5 of them with Cu (Z=29) to make a 15% doped system. Next "x spacegroup" changed the spacegroup and reduced the number of symmetry operations. At last the scf calculation was performed and solution converged.</div>
<div> The problem is, the density of states is reasonably different in these two above described cases. The density obtained from the VCA looks almost like the pure Ni DOS. Whereas that obtained from supercell approach is different. I need to know whether the ways of approach are correct or not in both cases. Should I consider the supercell approach more reliable over the VCA ? While constructing a supercell, is it needed to increase the number of k points also ?</div>
<div> Thank you very much in advance for helping me.</div>
<div><br></div>Prasenjit Roy<br>
</div><div>Radboud University, Nijmegen</div>