[Pw_forum] problem with supercell

2013-11-03 Thread Lorenzo Paulatto
Finally, I do not understand what you're doing with the magnetization and Hubbard U, but I'm don't know that method very well, hence I'll leave any comment to more expert people. good work! -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 4

[Pw_forum] linking quantum espresso with ACML library

2013-10-28 Thread Lorenzo Paulatto
On 10/28/2013 06:53 PM, Jan Gryko wrote: > I have no problem linking with fftw3, but when I try to link with the ACML, > it cannot find all subroutines. Are you sure you are linking correctly? best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33

[Pw_forum] q2r.x error

2013-10-27 Thread Lorenzo Paulatto
y hand? best > Is there anyway to bypass this issue in existing results? and calculate > lambda and Tc? Send the input files you are using and the header of the dyn files to this list for more accurate answers. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Pari

[Pw_forum] MPI

2013-10-09 Thread Lorenzo Paulatto
t; here? Dear Mike, if you oonly have 32 cpus it is a bad idea to run 4096 parallel processes. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te co

[Pw_forum] Problem of XCrySDen launching

2013-10-09 Thread Lorenzo Paulatto
> As already said, this is not the XCD mailing list. However, this looks like a possible conflict of variables between different versions of xcrysden. Clean up the hidden file $HOME/.bashrc : remove every possible reference to xcrysden, then reboot your computer (or at least the session) and t

[Pw_forum] problem with vaccum

2013-10-01 Thread Lorenzo Paulatto
repared some time ago: <https://www.dropbox.com/s/rl3vtbpajiiilni/qe-belfast13.tgz> If you also do PART1 and 3 it would not hurt... it takes a few hours. good work! -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http:/

[Pw_forum] time of scf and nscf computations

2013-09-28 Thread Lorenzo Paulatto
pression that SCF would diagonalize the empty bands with reduced accuracy, while NSCF would diagonalize all the bands with full accuracy. However, quickly checking the source code I found no trace of my memory; possibly I did not recall correctly. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC --

[Pw_forum] error in installing espresso-5.0.2

2013-07-03 Thread Lorenzo Paulatto
urs if you reuse a 4.3.x make.sys file in 5.x; to fix it either re-configure (but that pretty much kills the point of reusing make.sys) or add this line to make.sys: ELPA_LIBS_SWITCH = disabled best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44

[Pw_forum] error in routine electrons

2013-07-01 Thread Lorenzo Paulatto
On 07/01/2013 12:26 PM, Vincenzo Verdolino wrote: > Here is part of my input file Guess: It could be a combination of smearing set too high with the number of bands set too low. good work -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084

[Pw_forum] HgTe

2013-06-27 Thread Lorenzo Paulatto
y here: http://arxiv.org/abs/1305.5973 cheers -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5 --

[Pw_forum] HgTe

2013-06-26 Thread Lorenzo Paulatto
ation used for the reference isolated atom does not matter. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu

[Pw_forum] increasing the parameter npk

2013-06-26 Thread Lorenzo Paulatto
part of your question: you can modify it yourself in Modules/parameters.f90. Do not forget to recompile the code afterward. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/

[Pw_forum] Floating point exception

2013-06-25 Thread Lorenzo Paulatto
n your input. Why don't you try the latest version 5.0.3 of the code? best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5

[Pw_forum] a question about the adsorption of 3d atoms on a surface with spin -polarized calculations

2013-06-24 Thread Lorenzo Paulatto
. In the latter case, it is important to check that the final state actually retains the magnetization you asked for. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 2

[Pw_forum] Problem with Matdyn.x

2013-06-12 Thread Lorenzo Paulatto
ine K-POINTS crystal I do not know where you found it, but it does not exist in matdyn (scarce) documentation. of course, you can just run matdyn, no need to repeat the entire calculation. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084

[Pw_forum] "Raman activities" do not converge

2013-05-30 Thread Lorenzo Paulatto
non, epsilon, raman...) calculations it's better to use a much stricter conv_thr (i.e. 1.d-12). Sometimes I even reduce the hard-coded threshold (ethr) in PHonon/PH/set_defaults_pw.f90 bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 08

[Pw_forum] cell dm (1) for (8,0) nanotube

2013-05-24 Thread Lorenzo Paulatto
. On the other hand, if ibrav is different from 0, you have to check the corresponding cell definition in PW/Doc/INPUT_PW.txt and set celldm(1) accordingly (in bohr). best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skyp

[Pw_forum] Installation problem

2013-05-20 Thread Lorenzo Paulatto
re pickier than others, thank you for reporting we will fix it for the next release. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115,

[Pw_forum] fft 3D

2013-05-17 Thread Lorenzo Paulatto
wo; the wavefunction grids, each centered around its k point; the spherical grid, for ion-centered quantities that depend only on |G|. I'll go more in depth on these subjects if there is some demand. hth, bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6

[Pw_forum] fft 3D

2013-05-16 Thread Lorenzo Paulatto
the code snip you pasted have nothing to do with the fourier transform; they are a re-shuffling of the order of the plane waves. That said, the noise in "after" looks like a typical example of aliasing. Like here: <http://en.wikipedia.org/wiki/File:SquareWave.gif> bests -- D

[Pw_forum] Fermi Energy of Cu

2013-05-14 Thread Lorenzo Paulatto
aces. good work -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5 -- next part -- An HTML att

[Pw_forum] Details of partial density of states calculation

2013-05-14 Thread Lorenzo Paulatto
hope Stefano can be more helpful > Thanks, Natalie > you're welcome -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5

[Pw_forum] Details of partial density of states calculation

2013-05-13 Thread Lorenzo Paulatto
nd, there is a code that reconstructs the PAW charge (it's pawplot.f90 in the PP/src directory), I do not know if it has anything in common with what you did. I'll have a look at it if you wish, if it improves prowfc it is of course welcome. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC

[Pw_forum] Charge density at any k-point in the Brillouin zone

2013-05-04 Thread Lorenzo Paulatto
ut I cannot be sure either: I have not seen your input nor your output file. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4

[Pw_forum] Charge density at any k-point in the Brillouin zone

2013-05-02 Thread Lorenzo Paulatto
;inputpp is the index of the k-point, as listed in the output of pw.x. If you want to plot a specific point that is not in the list you will have to do a nscf calculation first. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / s

[Pw_forum] Is the results of structural relaxiation by pwscf accurate?

2013-04-29 Thread Lorenzo Paulatto
eliable? It is > And can you provide some suggestions? > Provide a sample case where do you think QE is failing best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-2

[Pw_forum] Not converging scf calculation

2013-04-28 Thread Lorenzo Paulatto
tly high energy every pseudopotential fails. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5

[Pw_forum] Question on the wsweight subroutine

2013-04-20 Thread Lorenzo Paulatto
id breaking symmetry) there are many ways to enforce these 2 simple rules, wsweights is just the most complicated > > I would really appreciate if you could shed some light on these issues. > > Many thanks once again! you can find a longer explanation in appendix C of this pre-prin

[Pw_forum] PW is frozen

2013-04-16 Thread Lorenzo Paulatto
mining a process you can detach the debugger and proceed with the next one. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 5

[Pw_forum] k point grid hybrid functionals

2013-04-15 Thread Lorenzo Paulatto
t confused about this topic, do > you have any paper to recommend me? > I do not have any paper where this is mechanism is explained in detail... maybe someone else who's reading this thread can recommend some reading? bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Unive

[Pw_forum] k point grid hybrid functionals

2013-04-15 Thread Lorenzo Paulatto
nd any {k+q} grid are just shifted and re-indexed w.r.t each other. I hope this helps, it is a bit confusing but it makes sense eventually. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone:+33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~pa

[Pw_forum] projwfc.x

2013-04-12 Thread Lorenzo Paulatto
On 04/12/2013 11:35 AM, Sakhrawi Taoufek wrote: > my input is : > &inputpp > prefix='fe', > outdir='$outdir', > degauss=0.01, > Emin=2.0, Emax=22.0, DeltaE=0.1, smoothing=0.3, > / > Starting from version 5.0 you must repla

[Pw_forum] d3.x on various q-points

2013-04-09 Thread Lorenzo Paulatto
=browse&path=%2F%2Acheckout%2A%2Ftrunk%2Fespresso%2FPHonon%2FD3%2Fd3_setup.f90&revision=10127>> and recompile the D3 code to fix the problem. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.imp

[Pw_forum] d3.x on various q-points

2013-04-09 Thread Lorenzo Paulatto
three inputs (pw, ph, d3) as well as any script you are using to run them and I'll try to reproduce your problem: without it is impossible to say anything definitive. cheers -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www:

[Pw_forum] How to add semicore state electron to valance state?

2013-04-07 Thread Lorenzo Paulatto
its all-electron configuration as [Xe] 4f14 5d3 6s2 or [Kr] 5s2 4d10 5p6 4f14 5d3 6s2 it does not matter. The states that go to valence have to be specified later on. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www:

[Pw_forum] fft_scalar_MOD_cfft3d at fft_scalar.f90:1218 error in vc-relax

2013-04-06 Thread Lorenzo Paulatto
rk properly if you let nrXX alone? best -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05

[Pw_forum] fix ions during relaxation / compile as a library?

2013-04-04 Thread Lorenzo Paulatto
asy to do. I know that some work in this direction has been done, and IS being done right now. i.e. PW can be used as a plugin for PLUMED (or PLUMED as a plugin for PW, depending on the website you're on); so maybe you can study and reuse this interface mecanism. cheers > Kind regards

[Pw_forum] How to set starting_magnetization?

2013-04-02 Thread Lorenzo Paulatto
e two distinct non-degenerate solutions could exist and be separated enough that the initial guess can matter. > > Thank you! > Best regards! you're welcome -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-in

[Pw_forum] Question on the wsweight subroutine

2013-03-29 Thread Lorenzo Paulatto
3 and so on. It takes some thought and some time to understand wsweight; if I remember correctly, Schwarz inequality is used < http://en.wikipedia.org/wiki/Cauchy%E2%80%93Schwarz_inequality> <http://en.wikipedia.org/wiki/Cauchy%E2%80%93Schwarz_inequality> bests -- Dr. Loren

[Pw_forum] Electronic entropy

2013-03-27 Thread Lorenzo Paulatto
entropy? I don't know :) > Thank you for your responses. > > Best regards, tchao -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4

[Pw_forum] Pw_forum Digest, Vol 69, Issue 23

2013-03-25 Thread Lorenzo Paulatto
a specified occupation? Is there a bug in the code? did you not see my answer? bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussi

[Pw_forum] nscf calculation error

2013-03-24 Thread Lorenzo Paulatto
to fixed. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05

[Pw_forum] Pw_forum Digest, Vol 69, Issue 20

2013-03-20 Thread Lorenzo Paulatto
p.s. you can also prepend "auto:" to the value of fildrho, which will automatically generate the name of the file, i.e. fildrho = "auto:drho" you must do it in both input files. -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275

[Pw_forum] Pw_forum Digest, Vol 69, Issue 20

2013-03-20 Thread Lorenzo Paulatto
g was quite high. Now, the stupid ascii file is not written any more, which made you think that the fildrho is not written, and the binary file is always written in the same place (outdir/_ph0). So, please, just believe me: use different names for the two files and everything will work. best

[Pw_forum] Problem in d3.x calculation, espresso-5.0.2

2013-03-20 Thread Lorenzo Paulatto
On 20 March 2013 10:10, Aaditya Manjanath wrote: > fildrho='Si.drho', > fild0rho='Si.drho', Dear Aaditya, you cannot use the same name for both files; they will be written one over the other. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit?

[Pw_forum] Unreadable .ps

2013-03-19 Thread Lorenzo Paulatto
ckerItemEdit&tracker_item_id=38&start=0> Anyway, just replace all the NaN with 0 to get a decent postscript. -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te co

[Pw_forum] Unreadable .ps

2013-03-18 Thread Lorenzo Paulatto
"***" instead of the actual number. Try to find the problem in the ps file, then go back to the fortran to fix it. It is not as hard as it sounds. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc

[Pw_forum] a question on pseudopotential generation

2013-03-12 Thread Lorenzo Paulatto
March 2013 14:46, Yue-Wen Fang wrote: > Dear Dr. Lorenzo Paulatto > Thank you! In the tutorial* INPUT_LD1.html*, it says that > rcore REAL > > Matching radius (a.u.) for the smoothing of the core charge. > If not specified, the matching radius is determined > by the condi

[Pw_forum] a question on pseudopotential generation

2013-03-12 Thread Lorenzo Paulatto
On 12 March 2013 13:14, Yue-Wen Fang wrote: > I really don not know how to understand "push slightly outwards"? make it larger > In addition, why to try rcore=2.0 but not other values? 2.0 is an example, try other values too bests -- Dr. Lorenzo Paulatto IdR @ IMPMC --

[Pw_forum] weighted sum over k-points

2013-03-08 Thread Lorenzo Paulatto
pin degeneracy put somewhere else. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05

[Pw_forum] Error in Execution while linked with Intel MKL - 11.0 64-bit : QE-5.0.2

2013-03-07 Thread Lorenzo Paulatto
ing" module (google it, if you want to knowmore). It is not the default yet, because it is not clear if the new method is supported by all compilers, but when it is supported I find it it more reliable. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone

[Pw_forum] Error in Execution while linked with Intel MKL - 11.0 64-bit : QE-5.0.2

2013-03-07 Thread Lorenzo Paulatto
Try to set MANUAL_DFLAGS = -D__ISO_C_BINDING in make.sys, than make clean and recompile. cheers -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 11

[Pw_forum] Changing nr3

2013-03-05 Thread Lorenzo Paulatto
tion in the parallel case. Maybe this feature got lost at that time, it should not be too difficult to put it back, but it requires time. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/

[Pw_forum] Changing nr3

2013-03-02 Thread Lorenzo Paulatto
You must set all the three dimensions: nr1, nr2 and nr3 for your setting to be taken in account. In the most recent versions of the code there is a warning in the output if you specify only one of the dimensions, I'll add a comment in the manual too. best regards -- Dr. Lorenzo Paulatt

[Pw_forum] gamma_only problem with bands.x

2013-02-28 Thread Lorenzo Paulatto
Thank you for reporting. It is a known bug of v5.0.2 see here for a fix: <http://qe-forge.org/pipermail/pw_forum/2013-January/100812.html> bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.imp

[Pw_forum] (no subject)

2013-02-25 Thread Lorenzo Paulatto
er manual available online < http://www.quantum-espresso.org/wp-content/uploads/Doc/ph_user_guide/> best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16

[Pw_forum] local and global minimums in relaxation

2013-02-15 Thread Lorenzo Paulatto
st). -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05

[Pw_forum] ph.x: Avoiding the recalculation of the band structure in distributed phonon dispersion jobs

2013-02-12 Thread Lorenzo Paulatto
hose that contain the irrep parallelism, but it may be possible to recompose the 2. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place

[Pw_forum] Disparities between HPC's

2013-02-07 Thread Lorenzo Paulatto
any case the result is that we have interference between different augmentation spheres. I've spent an afternoon in a naive effort to reimpose orthonormality on the G-space grid but could not finish before having to move on more urgent matters. If any good-soul has some time to dedicate t

[Pw_forum] error in make pw, after rev. 9674

2013-02-06 Thread Lorenzo Paulatto
c > I'm not sure if there is any simple solution and I cannot help you farther as it depends on details on your system which I do not have. It would be nie to see what "uname -a" returns on your system, to see how the harmless update of autoconf could break it bests -- Dr. Lore

[Pw_forum] negative frequency in high q-mesh

2013-01-31 Thread Lorenzo Paulatto
ultrasoft pseudopotential. Try to read something about this subject i.e. here < http://www.stanford.edu/~hkulik/www/Tutorials/Entries/2012/4/17_Convergence_101.html> (the "dual" is ecutrho/ecutwfc) bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0

[Pw_forum] converting UPF 2.0.1 into old versions

2013-01-29 Thread Lorenzo Paulatto
Dear John, it is not possible, but if you have the input (recent version UPF have the input included in the file) you can generate the pseudopotential with the 4.3.1 version of the ld1.x code bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275

[Pw_forum] changing total energy due to changing ibrav

2013-01-16 Thread Lorenzo Paulatto
ve me > some information on why the total energy changes? > Dear Tram, maybe you did not use 2 sets of k-points that map perfectly from one system to the other. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http:/

[Pw_forum] Relax Calculation for ZnO

2013-01-12 Thread Lorenzo Paulatto
is a semiconductor, not a metal. > > If it is not a metal, you can normally use fixed occupations. >1. After my relax calculation is done, should I run a scf calculation >again with atomic positions I get from relax calculation? > > You only need to for variable cell

[Pw_forum] setting grid size and memory management for pw scf

2013-01-12 Thread Lorenzo Paulatto
nnot just reduce the cutoff, you need to test convergence and use a converged value. If you are running out of memory you can either but a bigger computer or study a smaller system. In some specific case changing the parallelism scheme or some computation parameter can help. bests -- Dr. Lore

[Pw_forum] negative phonon frequencies again :: graphene

2012-12-20 Thread Lorenzo Paulatto
Dear Martin, which version of the code are you using? bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?de

[Pw_forum] questions about computing the third order coupling tensor implemented by d3.x code

2012-12-20 Thread Lorenzo Paulatto
m using a 4.2.1 version. > There was a bug with pools up to a certain version, use the latest one (and apply this fix: < http://qe-forge.org/gf/project/q-e/scmsvn/?action=browse&path=%2Ftrunk%2Fespresso%2FPHonon%2FPH%2Fsgam_ph.f90&r1=9630&r2=9724 >) bests -- Dr. Lorenzo Paul

[Pw_forum] questions about computing the third order coupling tensor implemented by d3.x code

2012-12-19 Thread Lorenzo Paulatto
Dear Tian, please see my post from February: < http://www.democritos.it/pipermail/pw_forum/2012-February/023255.html> best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paul

[Pw_forum] unit of mode vectors in matdyn.modes

2012-12-13 Thread Lorenzo Paulatto
om. If I remember correctly, those in the dyn files do not include it. > * the q-points (I guess fractional w.r.t. reciprocal unit cell) > cartesian axes in units of 2pi/celldm(1) > Thank you. > -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33

[Pw_forum] error: "ortho went bananas" in cp.x run

2012-12-13 Thread Lorenzo Paulatto
do *exactly* as it says on the web page: svn checkout http://qeforge.qe-forge.org/svn/q-e/trunk/espresso bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te c

[Pw_forum] contradictory results of electron-phonon interaction for different k-grids

2012-12-12 Thread Lorenzo Paulatto
d" phonon dispersions, even if they are not properly relaxed; when the symmetry is lower even a small displacement from the equilibrium will give some negative frequency. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype:

[Pw_forum] error while compiling quantum espresso-5.0.1 with ifort 12.0

2012-12-06 Thread Lorenzo Paulatto
84/mkl/include -c hills.c > > Deat Ruchir, your compiling environment seems to be mixed up. Try to restart from scratch, make clean than reconfigure with ./configure [any of your options] MPIF90=mpiifort FC=ifort CC=icc and recompile bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS &

[Pw_forum] K_POINTS automatic

2012-12-03 Thread Lorenzo Paulatto
efunctions from two artificial replicas of the graphene bilayer to hybridize with each other, which is an unphysical effect. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail:

[Pw_forum] Error while running tests and examples

2012-12-03 Thread Lorenzo Paulatto
Here they are (for now): *https://www.dropbox.com/s/2zkiebh2ss1owgx/all_test_pseudos.tgz* some pseudos from the examples may still be missing, but at least all the tests should run. Just uncompress the archive in the pseudo/ directory inside QE bests On 3 December 2012 15:05, Lorenzo Paulatto

[Pw_forum] Error while running tests and examples

2012-12-03 Thread Lorenzo Paulatto
the example and test pseudos somewhere I'll be glad to prepare it (for you, and for me too as I'm going to need sooner or later). bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr

[Pw_forum] K_POINTS automatic

2012-12-03 Thread Lorenzo Paulatto
etail for us to understand. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05 -- next par

[Pw_forum] K_POINTS automatic

2012-12-03 Thread Lorenzo Paulatto
han the band structure should be the same in the two case. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05 ---

[Pw_forum] dummy atom

2012-11-29 Thread Lorenzo Paulatto
On 29 November 2012 09:34, Maryam Jamali wrote: > Please let me know, could I define dummy atom in ATOMIC_POSITIONS in pw.x? > if it could, how? > Thanks in advance, > > What do you mean by "dummy atoms"? bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Unive

[Pw_forum] electron-phonon interaction error.

2012-11-22 Thread Lorenzo Paulatto
> What are these numbers? > > DO you meet this case? Could you tell me the reason or how to correct it? > Did your phonon dispersion look ok in the first place? Did you actually compute the el-ph coupling and found it to be uncorrect? best regards -- Dr. Lorenzo Paulatto IdR @ IMPM

[Pw_forum] binding energy

2012-11-21 Thread Lorenzo Paulatto
SPRESSO is not photoshop, i.e. you cannot just play around with parameters until you get something you like. You should understand the meaning of every variable in the input in order to learn something new. Refer to the user manuals in directory PW/Doc. -- Dr. Lorenzo Paulatto IdR @ IMP

[Pw_forum] Format of .pp postprocessing file

2012-11-14 Thread Lorenzo Paulatto
e files and xcrysden xsf files. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 place Jussieu 75252 Paris C?dex 05 -- next

[Pw_forum] SCF calculation don' converge for surface slab

2012-11-14 Thread Lorenzo Paulatto
On 14 November 2012 08:10, patriot pershing wrote: > Dear all > i have modified the mixing mode and smearging parameters and the total > energy converge but it take many iteration (14-20 iteration by relax step) > > It is normal if you have many atoms. bests -- Dr. Loren

[Pw_forum] SCF calculation don' converge for surface slab

2012-11-12 Thread Lorenzo Paulatto
urface, making convergence unstable. You coudl also try to set mixing_mode="local-TF" in &electrons, it helps sometimes. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/

[Pw_forum] Phonon calculations

2012-11-12 Thread Lorenzo Paulatto
omputes the dynamical matrix and store it in the fildyn files, then it also diagonalize it, but you can ignore this part if you wish. > > Thanks > bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.imp

[Pw_forum] Weird jump in phonon dispersion at Gamma point

2012-11-10 Thread Lorenzo Paulatto
Yes, it was another thread: < http://pwscf.org/pipermail/pw_forum/2012-October/099881.html> (please note that the mailing list archive has moved since october 4th). I'm happy you could sort it out -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33

[Pw_forum] Parallel Phonon calculations

2012-11-09 Thread Lorenzo Paulatto
m mailing list > Pw_forum at pwscf.org > http://pwscf.org/mailman/listinfo/pw_forum > -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 p

[Pw_forum] Weird jump in phonon dispersion at Gamma point

2012-11-09 Thread Lorenzo Paulatto
lo_to_split=.FALSE. ENDIF ! !!! CALL nonanal (nat, nat_blk, itau_blk, epsil, qhat, zeu, omega, dyn) i.e. set lo_to_split to .false. in every case and do not compute the non-analitical long-range term. Then recompile the phonon code best regards -- Dr. Lorenzo

[Pw_forum] Error in compiling espresso-4.0.3

2012-11-06 Thread Lorenzo Paulatto
eems that the latter was the case, maybe you did a small mistake while editing make.sys? Try to compare make.sys with make.sys.in and/or do a "make veryclean" and redo the compilation from scratch. -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 442

[Pw_forum] Relaxtion

2012-11-06 Thread Lorenzo Paulatto
On 6 November 2012 10:11, Giovanni Pizzi wrote: > try also to put the &electrons namelist before ions and cell. This > should solve the problem. > > That is it! I'm not sure it is stated in the documentation, but &electrons must always preced &ions and &cell.

[Pw_forum] pdos files connot bw generated by projwfc.x

2012-11-06 Thread Lorenzo Paulatto
On 6 November 2012 10:10, ??? wrote: > what at bottoom is my input for and thank you for your help!! > > Please post your output too. best regards -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.

[Pw_forum] Relaxtion

2012-11-05 Thread Lorenzo Paulatto
range character in your input (you can check with "cat -v input.file" and see if anything strange appears). In a desperate case, you may try to invert the order of the namelists &cell and &ions, with some compiler one of the two orders is not good, but I do not remember which. go

[Pw_forum] scf_must_converge

2012-11-05 Thread Lorenzo Paulatto
every dataset must be tested separately and its convergence thresholds determined w.r.t. the desired accuracy. -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te cour

[Pw_forum] Phonon DOS per polarization

2012-10-26 Thread Lorenzo Paulatto
not very hard to do, but a bit boring as it requires a rethinking of the code structures. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upmc.fr/~paulatto/ mail: 23-24/4?16 Bo?te courrier 115, 4 p

[Pw_forum] Questions about the DFPT review paper

2012-10-26 Thread Lorenzo Paulatto
gt; If I remeber correctly, replace \Delta\psi from eq 28 in eq 23 > > 2. Why the left side of Eq.25 has a null eigenvalue? > > It may be better explained in Messiah "Quantum mechanics" book (1962), have you tried to have a look there? bests -- Dr. Lorenzo Paulatto IdR @ I

[Pw_forum] Phonon DOS per polarization

2012-10-25 Thread Lorenzo Paulatto
problem at all. Finally, you can think of a specific criteria to discriminate the bands in a specific material, but implementing it is up to you. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www: http://www-int.impmc.upm

[Pw_forum] IFCs and Difference between unit cell and supercell

2012-10-20 Thread Lorenzo Paulatto
comes Gamma (weight 1/4) + X =a/2pi(1 0 0) (weight 3/4) in the FCC cell. In general it is not easy, you will need to think some time (let's say a few hours) about it. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +33 (0)1 44275 084 / skype: paulatz www:

[Pw_forum] The jobs always cut after couple iterations

2012-10-19 Thread Lorenzo Paulatto
528, 220) > Arrays for rho mixing 854.30 Mb (6998400, 8) > > Dear Ihsan, I would bet you are running out of memory. There should be some sign of it in the various output files that the queue system generates. bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Univer

[Pw_forum] Different optical mode frequency at the same q-vector

2012-10-16 Thread Lorenzo Paulatto
On 16 October 2012 05:09, Kyeong-hyun Park wrote: > Dear Dr. Lorenzo Paulatto > > Thank you for your suggestion though, unfortunately commenting out the > statement didn't resolve > my problem. > Well, I think your first problem was solved. However I'm not sure if y

[Pw_forum] Different optical mode frequency at the same q-vector

2012-10-15 Thread Lorenzo Paulatto
On 15 October 2012 07:24, Kyeong-hyun Park wrote: > Dear Dr. Lorenzo Paulatto and other users > > First of all, thank you for your advice. > Hi, in order to properly remove the anisotropy you must also comment out the line CALL nonanal (nat, nat_blk, itau_blk, epsi

[Pw_forum] Different optical mode frequency at the same q-vector

2012-10-15 Thread Lorenzo Paulatto
ode, but maybe you have not (the code is unnecessarily complicated) 4. I do not have a handy example to check out how to enable/disable LO-TO splitting, but if you want to send me the *dyn files I'll have a look bests -- Dr. Lorenzo Paulatto IdR @ IMPMC -- CNRS & Universit? Paris 6 phone: +

<    5   6   7   8   9   10   11   12   13   14   >