<html><head></head><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:16px"><div id="yui_3_16_0_ym19_1_1546601289982_7806">Dear users and developers</div><div id="yui_3_16_0_ym19_1_1546601289982_7807"><div id="yui_3_16_0_ym19_1_1546601289982_7820"><br id="yui_3_16_0_ym19_1_1546601289982_7808"></div><div id="yui_3_16_0_ym19_1_1546601289982_7821" dir="ltr">I am new to QE mailing list. I sent this email one day ago. but could not see it is to the list. Thus I am resending<br></div><div id="yui_3_16_0_ym19_1_1546601289982_7838"><br id="yui_3_16_0_ym19_1_1546601289982_7808"></div></div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7809">I
 need to calculate the second/third order force constant for a magnetic 
system, (ferromagnetic). Is it possible with Quantum espresso?  I 
searched in the information for ph.x, and couldn't find any related 
parameters.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7810"><br id="yui_3_16_0_ym19_1_1546601289982_7811"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7812">If I do scf.in including nspin=2 and starting magnetization and then the ph.x as normal way, will this give what I want?</div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7813"><br id="yui_3_16_0_ym19_1_1546601289982_7814"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7815">Thank you for your time</div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7816"><br id="yui_3_16_0_ym19_1_1546601289982_7817"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7818">Shayam</div><div dir="ltr" id="yui_3_16_0_ym19_1_1546601289982_7819">PhD candidate</div>Uof Regina,ca</div></body></html>