
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear Yin,<br>

<br>

You already got great suggestions. I would like to add one more option to calculate IR from ab initio MD. You could do scf calculation on each configuration obtained via MD to calculate Maximally Localized Wannier Functions (MLWF). From MLWF you could calculate

 the dipole moment of the simulation cell or even for single molecules inside the cell at each time step. It is straightforward to obtain IR  once you have the dipole moment at each time step.<br>

<br>

An example for this option can be found here:<br>

<br>

<a href="http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.98.247401" target="_blank">http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.98.247401</a><br>

<br>

Note that the authors of this paper implemented Wannier on the fly which I think is not available in Q.E. But the computational cost to re-track the trajectory and calculate MLWF should not be too expensive if one already affords the MD run.<br>

<br>

<br>

Mostafa Youssef<br>

MIT<br>

<br>

<br>

</div>

</body>

</html>