[Pw_forum] wfc files: heavy I/O, handling for restarts
Dear Prof. Paolo, Thanks a lot for your reply. > which QE code and which files are you referring to? < That is phonon code (ph.x) and the files I mentioned before are given below in detail. With the earlier versions up to 4.1.3, phonon code runs fine and I got following files (bold and underlined items are folders) in the tmp directory of a local disk. The nodes which are used by the codes are node186, 036, 139. node186:/tmpscratch/sksct8/tmp$ ls > *pbmno.save* _phpbmno.com1_phpbmno.dwf2 _phpbmno.igk3 > _phpbmno.prd3 > pbmno.wfc1 _phpbmno.com2_phpbmno.dwf3 _phpbmno.mixd1 > _phpbmno.recover > pbmno.wfc2 _phpbmno.com3_phpbmno.ebar1 _phpbmno.mixd2 > _phpbmno.recover2 > pbmno.wfc3 _phpbmno.dvkb31 _phpbmno.ebar2 _phpbmno.mixd3 > _phpbmno.recover3 > _phpbmno.bar1 _phpbmno.dvkb32 _phpbmno.ebar3 *_phpbmno.phsave** * * > _phpbmno.save* > _phpbmno.bar2 _phpbmno.dvkb33 _phpbmno.igk_phpbmno.prd1 > _phpbmno.bar3 _phpbmno.dwf1_phpbmno.igk2 _phpbmno.prd2 > node036:/tmpscratch/sksct8/tmp$ ls pbmno.wfc4 _phpbmno.dwf4 _phpbmno.mixd4 _phpbmno.recover4 pbmno.wfc5 _phpbmno.dwf5 _phpbmno.mixd5 _phpbmno.recover5 _phpbmno.bar4 _phpbmno.igk4 _phpbmno.prd4 _phpbmno.wfc4 _phpbmno.bar5 _phpbmno.igk5 _phpbmno.prd5 _phpbmno.wfc5 node139:/tmpscratch/sksct8/tmp$ ls pbmno.wfc6 _phpbmno.bar8 _phpbmno.igk7 _phpbmno.recover6 _phpbmno.wfc8 pbmno.wfc7 _phpbmno.dwf6 _phpbmno.igk8 _phpbmno.recover7 pbmno.wfc8 _phpbmno.dwf7 _phpbmno.prd6 _phpbmno.recover8 _phpbmno.bar6 _phpbmno.dwf8 _phpbmno.prd7 _phpbmno.wfc6 _phpbmno.bar7 _phpbmno.igk6 _phpbmno.prd8 _phpbmno.wfc7 However, with the new version 4.3.1, for the exactly same input files and job scripts I only get these files and nothing else : node045:/tmpscratch/sksct84/tmp$ ls *pbmno.save* pbmno.wfc1 *_ph0* node111:/tmpscratch/sksct84/tmp$ ls pbmno.wfc2 pbmno.wfc3 node092:/tmpscratch/sksct84/tmp$ ls pbmno.wfc4 pbmno.wfc5 node080:/tmpscratch/sksct84/tmp$ ls pbmno.wfc6 pbmno.wfc7 node072:/tmpscratch/sksct84/tmp$ ls pbmno.wfc8 Note that, the used nodes in this time are node045, node111, node092, node080, node072. So it is clear from the above example that somehow in the new version 4.3.1, _phpbmno.save and _phpbmno.phsave goes inside the directory "_ph0" ; And the same phonon calculation, which was running fine with the earlier version, now stops like this way (a bit abruptly and rudely, with out much informations or error messages) : Electric field: Dielectric constant Born effective charges in two ways Atomic displacements: There are 5 irreducible representations Representation 1 3 modes -T_1u G_15 G_4- To be done Representation 2 3 modes -T_1u G_15 G_4- To be done Representation 3 3 modes -T_1u G_15 G_4- To be done Representation 4 3 modes -T_1u G_15 G_4- To be done Representation 5 3 modes -T_2u G_25 G_5- To be done simply with this error : -- MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD with errorcode 0. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -- -- mpirun has exited due to process rank 2 with PID 8791 on node node111.cvos.cluster exiting without calling "finalize". This may have caused other processes in the application to be terminated by signals sent by mpirun (as reported here). -- [node045:11580] 6 more processes have sent help message help-mpi-api.txt / mpi-abort [node045:11580] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages Hope this email explains much better. Thanks and Regards, Saha SK -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110906/78efc62c/attachment.htm
[Pw_forum] wfc files: heavy I/O, handling for restarts
# Scatter wfc restart files awk '{ files_for[$1] = files_for[$1] " '$basename'.wfc" NR } END { for (host in files_for) print host, files_for[host] }' $PBS_NODEFILE \ | while read host files do ssh -n $host "cd $PBS_O_WORKDIR; mv $files $ESPRESSO_TMPDIR/" done # on master host, copy .save directory as well rsync -a $basename.save $ESPRESSO_TMPDIR mpirun -x ESPRESSO_TMPDIR \ -np $(wc -l < $PBS_NODEFILE) \ -machinefile $PBS_NODEFILE \ pw.x -inp input.txt > output.txt # Gather remote files uniq $PBS_NODEFILE \ | while read host do ssh -n $host "rsync -a $TMPDIR/ $PBS_O_WORKDIR/" done -- E.g. for a job with nodes=3:ppn=4 the scatter part would distribute the existing files pwscf.wfc{1..12} as follows: << Yes, this looks indeed cumbersome. It becomes more painful when one can not know a priori in which nodes his/her job will go, particularly when it is totally decided by the automatic queue decider, depending on the free nodes available. In such a situation, one's restarted job may go to a totally new set of nodes, and phonon calculation can not get necessary files to restart. Then restarting phonon calculation becomes more difficult. It seems, there is a more serious trouble in the recent version of QE. In the version before QE4.2, the QE codes used to replicate the same necessary files to the distributed local disks of all the nodes. In this case, at least phonon calculation can run smoothly instead of crashing. But in the recent version, phonon calculations just stop by complaining that the distributed .wfc files in one node are not visible by another node. If a quick remedy of this problem is not easy, then at least, for the time being, it is better to keep the earlier option of replicating the same .wfc files in all the nodes still working in the version 4.3.1. Other better option can be to implement the "WF_COLLECT" trick, also in phonon code, as it is already there for PW.x. Thanks and regards, Saha SK R&D Assistant JNCASR Bangalore 560064 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110905/ea2b77e7/attachment-0001.htm
[Pw_forum] LDA PZ USPP for Mg
Dear All, Would someone like to share a PZ-LDA USPP for Mg ? regards, Saha SK stewart at cnf.cornell.eduwrote: > Hi all, > > I was wondering if anyone had a reliable ultrasoft > pseudopotential for Mg in > the local density approximation. The periodic table > on the PWscf site has > only norm-conserving and a GGA ultrasoft. > > Thanks, Hi Derek, I can try generate it. My US GGA pseudopotential for Mg works quite well. It will take few days, hopefully, including tests. Bests, Eyvaz. -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110518/30742ec4/attachment.htm
[Pw_forum] Problem in paralization with QE4.3
Dear All, Has anybody implemented "wfcdir" option in the QE_ph code as like QE_pw code ? The exactly same job which was running okay in the older version of QE (4.2) is not running in the newest version of QE (4.3). Here is the error : Calculation of q =0.000 0.000 0.000 %% from diropn : error #21 %% from diropn : error #21 error opening pbmno.bar error opening pbmno.bar % stopping ... %%% from diropn : error #21 error opening pbmno.bar %% stopping ... % from diropn : error #21 error opening pbmno.bar % stopping ... stopping ... % from diropn : error #21 error opening pbmno.bar %% stopping ... -- MPI_ABORT was invoked on rank 3 in communicator MPI_COMM_WORLD with errorcode 0. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -- -- mpirun has exited due to process rank 5 with PID 4397 on Looking forward to your comments. Thanking you and with my best regards, Saha S. K. R&D Assistant JNCASR Bangalore - 560064 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110512/37aa1739/attachment.htm
[Pw_forum] matdyn.x
Dear Stefano, Thanks a lot for your quick reply. Yes, both set of the frequencies are obtained after imposing ASR. All of frequencies are not negative. Two of the modes among total 5 modes (triply generate) are negative, which are not very unusual in a BaTiO3-type cubic perovskite. regards, Saha S. K. On 5/11/11, Stefano de Gironcoli wrote: > are you imposing acoustic sum rule ? this change slightly the frequencies . > isn't strange that your frequencies are all negative (imaginary) ? > is your system stable ? > stefano >
[Pw_forum] matdyn.x
Dear All, I find a bit strange results when using matdyn.x to perform the interpolation in order to get a complete phonon dispersion plot. In case of my system, at Gamma point, ph.x prints out two noticeably distinguishable triply degenerate phonon modes of frequency values as follows : -18.9595 -18.9595 -18.9595 -11.8835 -11.8835 -11.8835 However, after using "matdyn.x" the above phonons become like : -18.6356 -18.6356 -18.6356 -18.6096 -18.6096 -18.6096 This looks quite strange to me. As far as I understand, the phonon interpolation should not change the frequency values at Gamma point or any other high-symmetry point (suppose X in my case) which have been already computed by the QE_phonon code. Interpolation means, it should only generate more frequencies in between Gamma and X point, not changing the frequency values in the exact Gamma and X point, which have been already printed by the ph.x run. Looking forward to your comments. Thanking you and with my best regards, Saha S. K. R&D Assistant, JNCASR Bangalore -560064
[Pw_forum] epsilon_0
Dear QE USERS, What QE calculates is the dielectric tensor at clamped ions (often referred to as \epsilon_\infty). The ionic contribution must be calculated and added to this in order to obtain the static dielectric tensor (\epsilon_0). Does the most recent version of QE print "\epsilon_0" somewhere in the output file ? Thanking you and with my best regards, Saha SK R&D Assistant JNCASR, Jakkur, Bangalore 560064 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110428/b81e5717/attachment.htm
[Pw_forum] LDA-USPP for Ge
Dear QE Friends, Does any one have LDA-USPP for Ge ? Does this scarcity mean that it is not good to use LDA for Ge ? Thanking you all, Saha SK R&D Assistant JNCASR Bangalore-560012
[Pw_forum] Ge LDA USPP
Dear QE Friends, Does any one have LDA-USPP for Ge ? Thanking you, Saha SK RA JNCASR Bangalore-560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110216/bbb31d29/attachment.htm
[Pw_forum] the files necessary for restarting still go into 'outdir'
Dear All, The doc file says : wfcdir : this directory specifies where to store files generated by each processor (*.wfc{N}, *.igk{N}, etc.). The idea here is to be able to separately store the largest files, while the files necessary for restarting still go into 'outdir' (for now only works for stand alone PW ) So here my question goes : Is this nice trick still works only for stand alone PW or in the most recent version, this trick has been extended to other cases, in particular to PH ? Thanking you and with warm regards, Saha SK JNCASR Bangalore 560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110127/be31ff9c/attachment.htm
[Pw_forum] Ferromagnetic Insulator is the Ground State
> Dear Prof. Andrea Dal Corso, > > Thanks a lot for your helpful reply. > > I have further question, what one should do if the ground state is > Ferromagnetic (FM) INSULATOR ? Is there any way ? > If the AFM state can be stabilized by symmetry, starting with opposite starting magnetizations, then the phonon code should work. If you need a constraint to stabilize the AFM ground state, then presently you cannot do a phonon calculation. HTH, Dear Professor, Thanks a lot for your kind reply. I am still a bit confused. Is this your above reply implying that even if I get ferromagnetic insulator as ground-state, then also I need to follow the above method (which is for ANTI-FERROMAGNETIC insulator) to estimate the Born-effective charges ? Thanking you and with my best regards, Saha SK JNCASR Bangalore 560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110119/a1d9e735/attachment.htm
[Pw_forum] Wrong effective charges in AFM insulator using "zue"
Dear Prof. Andrea Dal Corso, Thanks a lot for your helpful reply. I have further question, what one should do if the ground state is Ferromagnetic (FM) INSULATOR ? Is there any way ? Andrea wrote: The phonon with LSDA + constraint is not programmed. So it will not work. If the AFM insulator is the ground state then you should obtain it with noncolin=.true. without constraint. In this case ph.x at gamma should work and you should get the correct effective charges. If not please provide the complete input. * regards, Saha SK JNCASR Bangalore 560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110119/896782b5/attachment.htm
[Pw_forum] spin polarized insulator
Dear All, In the recent version of QE code, is this possible to compute the dielectric constant of a spin-polarized insulator ? Thanking you and with my best regards, Saha SK JNCASR Bangalore
[Pw_forum] HOW TO COMPUTE Effective charge for FM INSULATOR using PWSCF?
Dear all, The fixed occupations and lsda need nelup and neldw, and the phonon code with constrained magnetization is not yet available. Does it imply that the phonon calculation is not possible for fixed occupations and lsda regards, S. K. Saha JNCASR BANGALORE P.S. Dear Mighfar, it is okay, cheers! -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20101125/6ee5bcf9/attachment.htm
[Pw_forum] cell_dofree
Dear All, Is there any option in PWSCF which allow to move x, y and z axis, keeping volume unchanged ??? Thanking you all, Saha S.K. JNCASR Bangalore -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20101125/327fda39/attachment.htm
[Pw_forum] DFPT+U
Dear QE USERS, It is nice to know that the DFPT+U has been recently implemented. I hope, it will be public soon. If not, may I ask a developmental implementation or executables, so that I can also perform some test calculations for my systems ? Sometime carrying out FROZEN PHONON is boring. By the way, thanks a lot for the good news. Regards, Saha SK, R&D Assistant, JNC, Bangalore On Thu, Jul 15, 2010 at 12:14 PM, Andrea Floris wrote: > Dear Ali, > > The DFPT+U has been recently implemented, and I am doing some > calculations on test materials, to see > how it performs. It will be public in the future, both in > norm-conserving and ultrasoft pseudopotentials. > > Regards, > > Andrea > > -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100715/64c375cc/attachment.htm
[Pw_forum] GGA+U
Dear QE USERS, Yes, the code works fine when I change the name of element from "eu" to "Eu" inside the pseudopotential file. Thank you all for your kind reply. One more question : what is Hubbard_alpha ? Is this Lagrange multipliers? What is its typical value? Is this possible to specify J & U both, in the input file of the newest version of QE code or is this always U_eff ? regards, Saha SK -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100714/3c15394b/attachment.htm
[Pw_forum] GGA+U
Dear QE Users, I got the following error while used GGA+U for EuO: %% from set_hubbard_l : error # 1 pseudopotential not yet inserted %% stopping ... Inside source file, I have found that "Eu" is mentioned. However, I do not understand the meaning of inserting a pseudopotential. Thanks in advance for a helpful explanation. With best regards, Saha SK R&D Assistant JNCASR Bangalore 560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100713/ca2b7ad1/attachment.htm
[Pw_forum] PAW
Dear Stefano, Thanks a lot for your helpful reply. It works fine. Now the error is for Li-1s. * Li-1s: (18:05:29)...STOP 2 ld1 failed! Aborting... Should I use rmatch_augfun= 2.2 in this case ??? In the same way, is it possible to generate Eu PAW using ld1.x ??? With best regards, Saha SK -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100708/22cf1d72/attachment.htm
[Pw_forum] PAW
Dear QE USERS, I got the following error in the version 4.2 of QE : Running oxygen test... * generating * testing * testing with spin done Generating PAW datasets library * Be: (16:35:55)... plotting... done * C: (16:35:56)... plotting... done * Cl: (16:35:56)...STOP 2 ld1 failed! Aborting... Is there someone else who also has encountered the same error ? Does anyone has PAW data for other elements, for instance, Eu or Sm ? Is this possible to use ABINIT PAW inside PWSCF ? With best regards, Saha SK R&D Assistant JNCASR India -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100708/fb97a102/attachment.htm
[Pw_forum] Eu or Sm Pseudo Potential
> Dear All, > > Does anybody have Eu or Sm (LDA/GGA and Ultra-soft/Norm-conserving) pseudopotential? > > Looking forward to your kind reply. > > Thanking you and with my best regards, > > Saha SK > R&D assistant > JNCASR > Jakkur 560012 > India -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100625/9a018c2b/attachment.htm
[Pw_forum] Wu-Cohen GGA
Dear All, Is Wu-Cohen GGA implemented in QE? Would anyone then like to share his or her pseudopotentials? With my best regards, Saha S. K. R&D Assistant JNCASR Bangalore- 560012 -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100618/eb2c7330/attachment.htm
[Pw_forum] flow-chart in PWSCF
Dear All, First, let me thank Ngoc Linh Nguyen for his helpful reply. A part of his reply is as follows: In general, the scf processes are similar, but in QE a density convergence check is applied to accelerate convergence in scf caculation, i.e. see in PRB 38, 12807 (1988). This reference is indeed helpful. I have still one confusion: what, in fact, QE starts with as an initial density n(r) at the 0th iteration, i.e. what is n_0(r) Is that simply constant or r dependent or both r and R dependent Then, n0(r, R) = what ? Here, r and R are electronic and nuclei co-ordinates, respectively. Thanking you and with my best regards, Saha SK R&D assistant JNCASR Bangalore -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100408/218f0cb4/attachment.htm
[Pw_forum] flow-chart in PWSCF
Dear All, At first, let me take this opportunity to thank Prof. Paolo for his helpful reply. A part of his reply is given here: >> I mean, PWSCF check energy convergence or >> density scf convergence or both? > the latter. More detials in the QE paper. Now just to be fully sure, I would like to ask, does PWSCF follow the attached flow-chart (Please see the attached file) where density-convergence is being checked On the other hand, in the earlier link: http://cms.mpi.univie.ac.at/vasp/vasp/node199.html energy-convergence had been checked. I believe, both the flowcharts are more or less same within a technical difference. Am I right? Looking forward to your further comments. Thanking you and with my best regards, Saha SK R&D Assistant JNCASR Bangalore -- next part -- A non-text attachment was scrubbed... Name: flowchart.pdf Type: application/pdf Size: 65095 bytes Desc: not available Url : http://www.democritos.it/pipermail/pw_forum/attachments/20100326/6c7086b0/attachment-0001.pdf
[Pw_forum] flow-chart in PWSCF
Dear All, http://cms.mpi.univie.ac.at/vasp/vasp/node199.html This link has a flow chart describing the algorithms used in VASP. Is there such a flow-chart in PWSCF also? If not, then I would like to know whether PWSCF follows the same flow-chart or not? I mean, PWSCF check energy convergence or density scf convergence or both? It would be really nice if somebody send me a full flow chart which has been used in PWSCF. With my best regards, Saha S. K. JNCASR R&D Assistant Bangalore -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20100317/018b91d1/attachment.htm
[Pw_forum] Xcrysden
Dear all, Thanks a lot for your kind and helpful reply. Few things about Xcrysden are yet to be known. 1. Now I can print the thermometer after installing ImageMagick. However, it prints the thermometer separately. Please have a look on the attached file. Is it possible in XCrysden to print the thermometer together with the background "Xcrysden charge-density" plot 2. >No! XSF uses Carthesian coordinates (Angstrom units). Then, is there any other smart way to change Angstrom coordinates to crystal coordinates in Xcrysden??? The usual pen and paper (analytical and manual) method is not always time-saving especially for non-cubic (for instance, hexagonal) systems. 3. Another thing I would like to know, whether it is possible to increase the size of the "XYZ co-ordinate symbol" in Xcrysden??? One can easily zoom the picture of a molecule or a crystal in the Xcrysden and print the larger size picture because Xcrysden's print is an XWD window dump. But what about co-ordinate symbol??? Is it possible to make its size bigger? regards, SKS -- next part -- A non-text attachment was scrubbed... Name: mol-urea2D.jpg Type: image/jpeg Size: 4822 bytes Desc: not available Url : http://www.democritos.it/pipermail/pw_forum/attachments/20090403/fe5fcd38/attachment-0001.jpg
[Pw_forum] Xcrysden
Dear Users, I have a couple of questions regarding Xcrysden. 1. Usually xcrysden files use the following format: PRIMVEC 5.52460788000.000.00 0.003.90648780400.00 0.000.00 19.5324384880 PRIMCOORD 14 1 31 0.3031580950 0.00-0.3114762280 31 2.6436116560 1.9532436370-1.8580462820 31 -0.0518106620 0.00-3.8817233780 31 2.7623039400 1.9532436370-5.8597314410 31 0.00 0.00-7.8129756070 31 2.7623039400 1.9532436370-9.7662192440 33 1.4763662580 1.9532436370 0.3364242780 where atomic positions are not written in crystal co-ordinates. Is it possible to specify atomic positions in crystal co-ordinates in any of .xyz or .xsf xcrysden input files (except .pwi or .pwo files). 2. How to print the thermometer in Xcrysden??? Everytime I get an error saying "could not execute /usr/bin/convert". I'll be highly obliged if anyone kindly discuss the above issues. regards, SKS JNCASR Bangalore
[Pw_forum] DMFT+LDA
Dear all, Recently, in abinit, I have seen some on-going developments on DMFT+LDA. I'm a little curious. Has anyone already started or is going to start the similar kind of developmental work in PWSCF regards, S. K. Saha JNCASR Bangalore
[Pw_forum] USPP
> Dear Prithwish, during next week, I will add to the pseudopotential library (http://www.pwscf.org/pseudo.htm) the Ce pseudopotential we used to study cerium-oxide. As you know, in general, compounds based on rare earths may display quite intriguing properties and our pseudo has not been tested for systems other than cerias. You can refer to me if you have specific questions, I will be happy to contribute as far as I can. Yours sincerely, Stefano Dear all, Quite a long time ago, the above message was posted regarding Ce pseudopotential. Is anyone of you having an LDA version of the Ce pseudo which takes f-electron in account? Looking forward to your kind reply. regards, S. K. Saha JNCASR Bangalore 560012
[Pw_forum] Silent Calculations
Dear Sir, Sorry, for my statements. I did not want to report them as an absolute truth. My comments was only based on the Quantum esspresso code. In other packages, the LRT can be more efficient. I have not checked them yet. Anyway, I would like to thank Prof. Baroni for his excellent reply which will really help to clear people's misconceptions about LRT regards SKS JNCASR BANGALORE -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080811/815d1fbc/attachment.htm
[Pw_forum] Silent Calculations
On Sat, Aug 9, 2008 at 1:52 AM, Paolo Giannozzi wrote: > Only for Gamma point phonon one can use "Frozen > phonon" method because it is much faster. not true in general. Each linear-response calculation costs approx. as much as a self-consistent one, and you need at least as many scf calculations as linear-response calculations. Dear Sir, I agree with you that for small system each linear-response calculation costs approx. as much as a self-consistent one, and we need at least as many scf calculations as linear-response calculations. However, IN CASE OF very large system (say about 64 atoms per unit cell), I found that a single representation takes too much computer time compared to a single scf run. Moreover, in some cluster which have multiple cpu for each node and in some quad or octa core computer, always phonon paralization is much slower than manual scf paralization (I mean manually divide your all scf calculations in different processors). In case of Frozen phonon method, one can also accomplish the same goal by using less memory and disk space. However, for non-zone center phonon Linear response is the best way. There is still a scope in PWSCF to make Phonon calculation much faster by implementing the representation number parallization. regards, SKS JNCASR BANGALORE On Sat, Aug 9, 2008 at 1:52 AM, Paolo Giannozzi wrote: > > On Aug 8, 2008, at 18:41 , S. K. S. wrote: > > > Only for Gamma point phonon one can use "Frozen > > phonon" method because it is much faster. > > not true in general. Each linear-response calculation > costs approx. as much as a self-consistent one, and > you need at least as many scf calculations as > linear-response calculations > > Paolo > --- > Paolo Giannozzi, Dept of Physics, University of Udine > via delle Scienze 208, 33100 Udine, Italy > Phone +39-0432-558216, fax +39-0432-558222 > > > > -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080809/b3409ab5/attachment.htm
[Pw_forum] Silent Calculations
Dear Miguel, Yes, you are right. The linear response calculation needs much time even for calculating phonons at Gamma point. Only for Gamma point phonon one can use "Frozen phonon" method because it is much faster. The calculations can be got stalled for many reasons. You have not mentioned whether at last (I mean staying in the silence mode for few weeks) your jobs properly (without any error) get over or not??? If yes, then it can be due to the problem in your parallel computer - may be out of 16 processors some particular processor was down for some time. If the answer of the above question is no, then the memory of your computer is not sufficient enough to tackle a large system with lot of atoms because it needs large memory and as well as generates a huge tmp file. All input parameters (deguass, la2F, filedvscf ..) also should be given correctly for electron-phonon calculation. regards, SKS JNCASR Bangalore On 8/8/08, Miguel Mart?nez wrote: > > Dear everybody, > > I've been running into a few issues as of late, when running some > calculations, it seems the calculations have stalled for a pretty long > period of time. > > The one that troubles me less, taking into account that it's pretty > expensive, is an electron-phonon calculation. This is a pretty large > system (20 atoms in the primitive cell) and computing the phonons at > gamma took ages (well two days using 4 processors). It would be great > if some output were written just to remind you the code is still > working, maybe something related to convergence. But I find somewhat > disturbing that no line was written since the electron phonon > calculation began. > > My other issue is a very similar stalling in a Raman calculation. No, > it's not the same structure. Here, I lack from experience, since this > is the first Raman calculation I've done outside the examples. > > Self-consistent input is: > &control >calculation='scf' >restart_mode='from_scratch', >prefix='sih4', >tstress = .true. >tprnfor = .true. >pseudo_dir = '/home/wmbmacam/pseudo/', >outdir='/scratch/wmbmacam/fdd2-raman/' > / > &system >ibrav= 10, celldm(1) = 12.65898481, celldm(2) = 1.07100858, > celldm(3) = 0.621208699, nat=10, ntyp= 2, >ecutwfc = 55.0 >nr1= 72, nr2= 72, nr3= 72 > / > &electrons >conv_thr = 1.0d-8 >mixing_beta = 0.7 > / > ATOMIC_SPECIES > Si 28.0855 Si.pbe-rrkj.UPF > H 1.0079 H.pbe-vbc.UPF > ATOMIC_POSITIONS {crystal} > Si 0.614368447799 0.385631552201 0.614368447799 > Si 0.864368447799 0.635631552201 0.864368447799 > H 0.442372877617 0.836174882136 0.364427192728 > H 0.364427192728 0.357025047520 0.442372877617 > H 0.892974952480 0.885572807272 0.413825117864 > H 0.413825117864 0.807627122383 0.892974952480 > H 0.378768583324 0.338682246661 0.003482898532 > H 0.003482898532 0.279066271483 0.378768583324 > H 0.970933728517 0.246517101468 0.911317753339 > H 0.911317753339 0.871231416676 0.970933728517 > K_POINTS {automatic} > 12 12 12 1 1 1 > > > The phonon input is: > > phonons of Sih4 at Gamma > &inputph > tr2_ph=1.0d-16, > prefix='sih4', > epsil=.true., > trans=.true., > lraman=.true., > elop=.true., > amass(1)=28.0855, > amass(2)=1.00794, > outdir='/scratch/wmbmacam/fdd2-raman/', > fildyn='sih4.dynG', > fildrho='sih4.drho', > / > 0.0 0.0 0.0 > > > The calculations stalled right after the effective charges calculations. > After that, it spent 24 hours without writing anything (it needed 10 > hours to compute the effective charges). No convergence, no threshold, > no crashes, no nothing. > > Is this behaviour normal? Or are my compiled binaries just telling me I > should stop bothering them in August? > > Finally, I just wanted to note that the amasses in example 15 are > seemingly wrong. While the masses in the scf run are fine, the ones > used in the phonon calculations correspond to GaAs instead of AlAs. > > Regards, > > Miguel > > -- > > Miguel Mart?nez Canales > Dto. F?sica de la Materia Condensada > UPV/EHU > Facultad de Ciencia y Tecnolog?a > Apdo. 644 > 48080 Bilbao (Spain) > Fax: +34 94 601 3500 > Tlf: +34 94 601 5326 > > > "If you have an apple and I have an apple and > we exchange these apples then you and I will > still each have one apple. But if you have an > idea and I have an idea and we exchange these > ideas, then each of us will have two ideas." > > George Bernard Shaw > > ___ > Pw_forum mailing list > Pw_forum at pwscf.org > http://www.democritos.it/mailman/listinfo/pw_forum > -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080808/77a0fabe/attachment.htm
[Pw_forum] not an easy system to treat
Dear Adriano, Thanks a lot for your kind help. I'll wait for your next reply. regards, SKS JNCASR -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080529/90512d69/attachment.htm
[Pw_forum] nelup and neldw
Dear all, The fixed occupations and lsda need nelup and neldw, and the phonon code with constrained magnetization is not yet available. Does it imply that the phonon calculation is not possible for fixed occupations and lsda regards, S. K. Saha JNCASR BANGALORE -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080529/9b526b28/attachment.htm
[Pw_forum] XCRYSDEN
Dear PWscf Users, Is there any option in XCRYSDEN software where we can increase the size or the font of the "Coordinate System" logo which usually comes below in the left side corner of any XCRYSDEN Picture? regards, S. K. S. JNCASR Bangalore -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080522/5a5411ed/attachment.htm
[Pw_forum] format UPF
Dear All, Is there any already exits converter to convert a pseudopotential, generated by Octopus, into UPF format?? regards, SKS JNCASR India -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080407/a7fb9ae2/attachment.htm
[Pw_forum] Gd LDA pseudo
Dear Sir, Thanks a lot for your kind reply. First time, I'm trying to generate a pseudo myself. Being newbie, some thing I can not understand clearly. As for example, in pwscf suite it's given for Os that * &input title='Os', zed=76., rel=1, iswitch=3, config='[Xe] 4f14.0 5d6.0 6s2.0 6p0.0', dft='LDA' / &inputp lloc=0, pseudotype=3, file_pseudopw='Os.pz-rrkjus.UPF', nlcc=.true., rcore=0.9, author='ADC', / 4 6P 2 1 0.00 0.00 2.60 2.60 1 5D 3 2 6.00 0.00 1.90 2.30 1 5D 3 2 0.00 -0.20 1.90 2.30 1 6S 1 0 2.00 0.00 2.50 2.50 1 * Though I can understand the meaning of the most of the input parameters, however, in last part , why the numbers are like these 6P 0.00 2.60 2.60 1 5D 0.00 1.90 2.30 1 5D -0.20 1.90 2.30 1 6S 0.00 2.50 2.50 1 , I can not understand clearly. I know that for Gd it will be like that 6P 21 0.00 5D 32 1.00 5D 320.00 6S 102.00 but what will be the other numbers for Gd I'll be highly obliged if somebody kindly help me to understand the above points. First, I'll test the generated Gd pseudo and if it works well I'll defintely give it back to PWSCF website. Look forward to your kind help. With my best regards, SKS JNCASR BANGALORE -- Forwarded message -- From: Dal Corso Andrea Date: Feb 27, 2008 8:29 PM Subject: Re: Gd LDA pseudo To: "S. K. S." Dear S.K.S. Presently I do not have a Gd pseudopotential. Please post this question to the pw_forum mailing list, maybe someone will be able to provide it. Best wishes Andrea -- Forwarded message ---------- From: Dal Corso Andrea Date: Feb 27, 2008 8:29 PM Subject: Re: Gd LDA pseudo To: "S. K. S." Dear S.K.S. Presently I do not have a Gd pseudopotential. Please post this question to the pw_forum mailing list, maybe someone will be able to provide it. Best wishes Andrea On Wed, 2008-02-27 at 20:19 +0530, S. K. S. wrote: > > >Respected Sir, > >May I get a Gd pseudo (LDA) for PWSCF from >you or any of your students??? > >I know that the rare earth is poorly described >by the plane waves. However, there are lot of DFT calculations >on Lanthanides and Ceria in literature. > >Looking forward to your kind help. > > With my best regards, > SKS > JNCASR > Bangalore > > P.S. I'll be very happy to get someone as a coauthor if he/she > helps me by providing the above pseudo. -- Andrea Dal CorsoTel. 0039-040-3787428 SISSA, Via Beirut 2/4 Fax. 0039-040-3787528 34014 Trieste (Italy) e-mail: dalcorso at sissa.it On Wed, 2008-02-27 at 20:19 +0530, S. K. S. wrote: > > >Respected Sir, > >May I get a Gd pseudo (LDA) for PWSCF from >you or any of your students??? > >I know that the rare earth is poorly described >by the plane waves. However, there are lot of DFT calculations >on Lanthanides and Ceria in literature. > >Looking forward to your kind help. > > With my best regards, > SKS > JNCASR > Bangalore > > P.S. I'll be very happy to get someone as a coauthor if he/she > helps me by providing the above pseudo. -- Andrea Dal CorsoTel. 0039-040-3787428 SISSA, Via Beirut 2/4 Fax. 0039-040-3787528 34014 Trieste (Italy) e-mail: dalcorso at sissa.it -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080228/7cae5887/attachment-0001.htm
[Pw_forum] Vanderbilt USPP
Dear Users, Can anyone please tell me whether the Vanderbilt USPP for Bi, Sm, Gd and Dy are available or not ? Regards, SKS PhD Jncasr India -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20071214/1aeb2595/attachment.htm
[Pw_forum] Doc/INPUT_PW: A compensating jellium background is inserted
Dear all, It's written in the documentation that a compensating jellium background is inserted. I'll be highly obliged if someone kindly tell me the name of the subroutine where it has been inserted. Yours sincerely, S. K. S R&D Assistant. Jncasr India -- next part -- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20071118/6897937c/attachment.htm