[Pw_forum] Phonon DOS

Subhradip Ghosh subhra at physics.rutgers.edu
Fri May 30 14:50:34 CEST 2003


Hi,

I am trying to compute Phonon DOS using the code matdyn.f90. I found that
it worked for an insulator but in my case the system is metallic and while
I ran the code there is some output generated alongwith an error message

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from matdyn : error #        -1
     non-analytic term for q=0 missing !

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I traced it back to the code and found that this error message pops up in
case we have no dielectric tensor(probably) written in the real-space
force constant
file. There is a logical switch in the force constant file which comes out
to be F in my case and therefore there are no entries of the variable
epsil. I have no idea what the remedy is? Is it a serious error or should
I ignore it?

Thanks in advance,

Subhradip

PS: My system is not a simple metal, rather a metallic alloy.

############################################################################
		     Dr	Subhradip Ghosh
		     Post Doctoral Fellow
		Condensed matter theory group
	        Department of Physics and Astronomy
############################################################################
Office:                                     | Residence:
--------------------------------------------|-------------------------------
Serin E255                                  |132, South Adelaide Avenue
Department of Physics and Astronomy         |Apartment# 2A
Rutgers, The State University of New Jersey |Highland Park
136 Frelinghuysen Road                      |NJ 08904
Piscataway, NJ 08854-8019                   |USA
USA                                         |
Email address: subhra at physics.rutgers.edu,  |Telephone # :(732)4489266
               subhradip at yahoo.com          |
Telephone #  : (732)445-4604                |
##############################################################################




More information about the users mailing list