<div dir="ltr"><div>Dear all, <br><br></div><div>I found out that it was impossible de make a phonon calculation without calculating the Born effective charge in a semi-conductor even if we set  epsil = .false.<br><br></div><div>This is due to the routine prepare_q.f90 that is called inside do_phonon.f90<br></div><div>In the routine there is<br>     IF ( lgamma ) THEN<br>        !<br>        IF ( .NOT. lgauss ) THEN<br>           !<br>           ! ... in the case of an insulator at q=0 one has to calculate<br>           ! ... the dielectric constant and the Born eff. charges<br>           ! ... the other flags depend on input<br>           !<br>           epsil = .TRUE.<br>           zeu = .TRUE.<br>           zue = .TRUE.<br><br><br></div><div>This means that if we compute q=Gamma and if we do not have gaussian-smearing (i.e. an semi-cond or insulator), then epsil is automatically set to TRUE.<br><br></div><div>I know that physically one should have such LO/TO splitting but the user should be able to choose. <br><br></div><div>Maybe this forcing could be reported in the input variable ? or simply put the default value to .TRUE. instead of false and not enforcing that rule?<br><br></div><div>What do you think?<br><br></div><div>Best, <br><br></div><div>Samuel Ponce, <br>Department of Materials, University of Oxford</div><div><br></div></div>