Re: [Wien] wien2venus.py problem
Dear Yasir Ali Please use XCrysden for calculating the 3D electron density distribution in WIEN2k. It can handle spin polarized cases as well. Regards Sharat Chandra On Thu, 18 Jul 2013, Yasir Ali wrote: Date: Thu, 18 Jul 2013 22:16:35 -0700 (PDT) From: Yasir Ali yasiralikhan...@yahoo.com Reply-To: A Mailing list for WIEN2k users wien@zeus.theochem.tuwien.ac.at To: wien wien@zeus.theochem.tuwien.ac.at Subject: [Wien] wien2venus.py problem I am facing problem in wien2venus.py. i used wien2venus.py 50 50 50 to generate case.rho3d file of electron density. When i used to plot this generated file, i.e., case.rho3d in VESTA, it does not give plot but only give an empty cube. So any one can tell me what is the problem? Regards: Yasir Ali ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] Initialization issue with core leak out of sphere And RMT Sums Greater than NN-DIST
Have you given the lattice parameters in Bohr units or in nm units? I Bohr = 0.5291772108 A On Mon, 5 Dec 2011, Qiwen YAO wrote: Date: Mon, 05 Dec 2011 19:24:36 +0900 From: Qiwen YAO Yao.Qiwen at nims.go.jp Reply-To: A Mailing list for WIEN2k users wien at zeus.theochem.tuwien.ac.at To: A Mailing list for WIEN2k users wien at zeus.theochem.tuwien.ac.at Subject: [Wien] Initialization issue with core leak out of sphere And RMT Sums Greater than NN-DIST Dear all, I have been facing a few compounds that I am unable to do the initialization (using w2web) for the calculation. The problem is: I was facing both core leak out of sphere, and RMT Sums Greater than NN-DIST at the same time. A recent example: I am trying to do an calculation initialization on the LuFe2O4 (R-3m group) crystal (cell parameters are from the Pearson's Crystal Data sheet 382676, a=b=0.3435nm, c=2.525nm, Alpha=Beta=90 degree, Theta=120 degree. Atomic positions: O1: 0 0 0.1295 Fe: 0 0 0.2141 O2: 0 0 0.2914 Lu: 0 0 0 ). Using the WIEN default RMT values and the default ENERGY to separate core and valence states (-6) , I was facing core electrons leak out of MT-sphere message for all four atoms. So I go back to the x start step and change the separation energy to the lowest possible value (-10 in this case) - if going any lower than -10, I will be seeing the error message after running x start, like these: forrtl: severe (24): end-of-file during read, unit -4, file stdin error: command /home/wen/WIEN2K11/wien2k11install/lstart lstart.def failed So, the natural step I took next is to increase the RMT values for the atoms correspondingly (according to each atom's core leak value). Then I am facing an another problem: If I am to increase the RMT values too large, then in the x nn step I will be seeing errors like ERROR !!! RMT( 1)=0.53000 AND RMT( 1)=0.53000 SUMS TO 1.06000 GT NNN-DIST= 0.84061 But if the RMT value is not increased enough then the core leaking message will still be there. So how would I solve such problems? I have been facing these for a few compounds - for example, the other compound I tried was BiAlO3 which would lead me to the same bottle-neck problem. Is there an automatic approaching way that I can take to get rid of both error message at the same time? For cases like these, if there is no practical solution to the initialization would that means the unit cell parameters were out of touch (reality)? Thank you for reading. I have been searching the archive for a solution but no luck. Kind regards, Wen ** Dr QiWen YAO JSPS Fellow Multifunctional Materials Group Optical and Electronic Materials Unit Environment and Energy Materials Research Division National Institute for Materials Science 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan Phone: +81-29-851-3354, ext. no. 6482, Fax: +81-29-859-2501 ** ___ Wien mailing list Wien at zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien -- Dr. Sharat Chandra %%%_%%___%%% Home Address %%% Scientist Asst. Prof. HBNI (_) ___ ___ / _ \ ___ 117, Fifth Avenue Materials Science Group, | |/ _ \ / __\| |_| |' _ \ DAE Township IGCAR, Kalpakkam,| | (_) | |__,| _ | | \| Kalpakkam 603102, Tamil Nadu, INDIA|_|\__/ |\___/|_| |_|_|603102 TN INDIA Ph: +91 44 27480500*22167 / /| | +91 44 27481521 email: sharat at igcar.gov.in | |_| | sharat.c at gmail.com \___/%%%
[Wien] About optic code in WIEN2k_11.1 (Release 5/4/2011)
Hi I have a report to make about the optic code included in WIEN2k_11.1 (Release 5/4/2011). I am using the WIEN2k on a Fedora 12 system with ifort (IFORT) 10.1 20080602. I am calculating the optical constants for Si diamond structure crystal. (lattice parameter 10.270700 atu, Si at (1/8,1/8,1/8). My Si.inop file is given below (just the default file from case.inop, where I have put the proper number of bands in line 2 for NBvalMAX): --- 9 1 number of k-points, first k-point -5.0 3.0 4Emin, Emax for matrix elements, NBvalMAX 2 number of choices (columns in *outmat): 1 Re xx 3 Re zz OFF ON/OFF writes MME to unit 4 Choices: 1..Re xx 2..Re yy 3..Re zz 4..Re xy 5..Re xz 6..Re yz 7..Im xy 8..Im xz 9..Im yz - When I run the command: x optic, I get the following error: emin,emax,nbvalmax -5.00 3.00 4 creating ust kmin,kmax -16 -16 -16 16 16 16 forrtl: severe (174): SIGSEGV, segmentation fault occurred Image PCRoutineLineSource optic 00428975 planew_ 166 planew_tmp.f optic 0043AAA7 mom_mat_ 588 sph-UP_tmp.f optic 0047 MAIN__447 opmain.f optic 0040328E Unknown Unknown Unknown libc.so.6 0030BE61EB1D Unknown Unknown Unknown optic 004031A9 Unknown Unknown Unknown 8.090u 0.232s 0:08.33 99.8% 0+0k 0+5840io 0pf+0w error: command /usr/local/wien2k111/optic optic.def failed -- This input file runs properly with the optic package included in WIEN2k_09.2 (Release 29/9/2009). I am not able to understand where is the error coming from as the codes in both the versions are quite different. Regards Sharat Chandra
[Wien] charge density of a single atom
Hi I am interested in obtaining the variation in the charge density with respect to r, for isolated Y atom. For this I have seen the convergence of total energy with respect to the unit cell length, RMT*kmax and Gmax and I use a single k-point for all the calculations. I obtain 15 a.u. as the converged unit cell parameter. I have put the atom in the centre of the unit cell. scf shows that there are 11 electrons in the valence band, ie., 4s2, 4p6, 4d1, 5s2. The input for lapw2 in electron density calculation is the range of the valence electrons energy obtained in the scf file. I carry out the calculation of the 3D charge density using XCrysden, using RHO, ANG, VAL, NODEBUG as the inputs for lapw5 program. I calculate the the distribution of the charge density from the (1/2,1/2,1/2) to (1,1/2,1/2) direction along the single line to get the rho(r) versus r on a 300 point grid from the xsf data save by XCrysden. When I perform an integration of this data (by integrating the r versus 4*pi*r^2*rho(r) data) to obtain the total charge, I get ~20 electrons. I should be getting either 11 electrons, or 39 electrons, which is the atomic number of Y, what I get is neither here nor there. Can you please tell me what I am doing wrong? Is it that I have to use AIM to get the real charge density, or any other method is to be used? Regards Sharat Chandra
[Wien] mpich2 with Wien2k
Dear Lawrence We have been able to get the WIEN working only with the mpich1 till now and we too get the same error as you have reported. The error is mainly due to the error in the devices that get made during the compilation of mpich2. However, WIEN works properly with lammpi package. We have the mpich1/2 and lam on our 64 node cluster machine with Redhat 4E and Intel ifort 9/mkl8 compilers. Till now we have found the mpich1 or lam mpi to be compatible with most of the dft packages. Hope this helps you. Regards Sharat Chandra On Fri, 18 Apr 2008, Laurence Marks wrote: Has anyone got this working? I'm getting an error in CALL SL_INIT(ICTXTALL, 1, NPE) in lapw1 of Fatal error in MPI_Comm_size: Invalid communicator, error stack: MPI_Comm_size(110): MPI_Comm_size(comm=0x5b, size=0xa4ee68) failed MPI_Comm_size(69).: Invalid communicator (it might be some bad hardware/drivers)
[Wien] Where is the MKL?
Hi You have to run the proper mklvars*.sh or mklvars*.csh file depending on your shell (whether bash or csh/tcsh) at the login time (put it in the .bashrc or .tcshrc). Regards Sharat Chandra