%% ****** Start of file apstemplate.tex ****** %
%%
%%
%%   This file is part of the APS files in the REVTeX 4 distribution.
%%   Version 4.1r of REVTeX, August 2010
%%
%%
%%   Copyright (c) 2001, 2009, 2010 The American Physical Society.
%%
%%   See the REVTeX 4 README file for restrictions and more information.
%%
%
% This is a template for producing manuscripts for use with REVTEX 4.0
% Copy this file to another name and then work on that file.
% That way, you always have this original template file to use.
%
% Group addresses by affiliation; use superscriptaddress for long
% author lists, or if there are many overlapping affiliations.
% For Phys. Rev. appearance, change preprint to twocolumn.
% Choose pra, prb, prc, prd, pre, prl, prstab, prstper, or rmp for journal
%  Add 'draft' option to mark overfull boxes with black boxes
%  Add 'showpacs' option to make PACS codes appear
%  Add 'showkeys' option to make keywords appear

\documentclass[reprint,amsmath,amssymb,aps,prb,superscriptaddress]{revtex4-1}
%\documentclass[aps,prl,preprint,superscriptaddress]{revtex4-1}
%\documentclass[aps,prl,reprint,groupedaddress]{revtex4-1}

% You should use BibTeX and apsrev.bst for references
% Choosing a journal automatically selects the correct APS
% BibTeX style file (bst file), so only uncomment the line
% below if necessary.
%\bibliographystyle{apsrev4-1}
\usepackage{textcomp,pxfonts}
\usepackage{amsmath}
\usepackage{graphicx}
\usepackage{multirow}
%\usepackage{authblk}

\usepackage{color}
%\newcommand{\hilight}[1]{\colorbox{yellow}{#1}}



%\usepackage{natbib}
\begin{document}

% Use the \preprint command to place your local institutional report
% number in the upper righthand corner of the title page in preprint mode.
% Multiple \preprint commands are allowed.
% Use the 'preprintnumbers' class option to override journal defaults
% to display numbers if necessary
%\preprint{}

%Title of paper

\title{Anharmonic lattice dynamics of Ag$_2$O studied by inelastic neutron scattering and first principles molecular dynamics simulations\\
Supplemental Material}

% repeat the \author .. \affiliation  etc. as needed
% \email, \thanks, \homepage, \altaffiliation all apply to the current
% author. Explanatory text should go in the []'s, actual e-mail
% address or url should go in the {}'s for \email and \homepage.
% Please use the appropriate macro foreach each type of information

% \affiliation command applies to all authors since the last
% \affiliation command. The \affiliation command should follow the
% other information
% \affiliation can be followed by \email, \homepage, \thanks as well.

\author{Tian Lan} 
\email[]{tianlan@caltech.edu}
\affiliation{\mbox{Department of Applied Physics and Materials Science, \nolinebreak California Institute of Technology, Pasadena, California 91125, USA}}
\author{Chen W. Li}
\affiliation{Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA}
\author{J. L. Niedziela}
\affiliation{Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA}
\author{Hillary Smith}
\affiliation{\mbox{Department of Applied Physics and Materials Science, \nolinebreak California Institute of Technology, Pasadena, California 91125, USA}}
\author{Douglas L. Abernathy}
\affiliation{Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA}
\author{George R. Rossman}
\affiliation{Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California 91125, USA} 
\author{Brent Fultz}
\affiliation{\mbox{Department of Applied Physics and Materials Science, \nolinebreak California Institute of Technology, Pasadena, California 91125, USA}}





\date{\today}


% insert suggested PACS numbers in braces on next line



%\pacs{78.70.Nx, 78.30.-j, 63.20.kg, 63.20.Ry}

% insert suggested keywords - APS authors don't need to do this
%\keywords{}

%\maketitle must follow title, authors, abstract, \pacs, and \keywords
\maketitle

% body of paper here - Use proper section commands
% References should be done using the \cite, \ref, and \label commands


To separate the effects of quasiharmonicity and explicit anharmonicity,
the mode frequency $\omegaup_j=\omegaup_j(V,T)$ is expressed as a function 
of volume and temperature
\begin{equation}
\label{eq:weight}
\left(\frac{\partial \ln \omegaup_j}{\partial T}\right)_P = - \gammaup_j\betaup +
\left(\frac{\partial \ln \omegaup_j}{\partial T}\right)_V
\end{equation}
where $j$ is the phonon mode index, $\betaup$ is the volume thermal expansivity 
and $\gammaup_j$ is the mode Gr\"{u}neisen parameter. 
The left-hand side %, as obtained in our Raman spectra shown in Fig. \ref{fig:Raman}, 
is the temperature-dependent  frequency shift at constant pressure, 
and includes contributions from both  quasiharmonicity and explicit anharmonicity. 
The first term on the right-hand side, the isothermal frequency shift as a function of pressure, 
is the quasiharmonic contribution to the frequency shift. 
The second term on the right is the pure temperature 
contribution to the frequency shift from the explicit anharmonicity.
From the difference of the isobaric and isothermal frequency shifts, 
the explicit anharmonicity can be determined.

In a molecular dynamics simulation, the quasiharmonic contribution can be evaluated explicitly by turning off the temperature-dependent anharmonicity. In principle, this method is equivalent to the QHA method implemented with  self-consistent lattice dynamics,\cite{ScF3} 
and in practice we have found this to be true.
For example, we performed simulations at 40\,K for volumes characteristic of 400\,K, 
which produced a pressure of 0.45\,GPa. This calculation therefore removed the temperature effect while preserving the quasiharmonic volume effect at 400 K. The corresponding phonon DOS curves from MD calculations are shown in Fig.~\ref{DOSSM}. By comparing the phonon spectrum of a simulation at 40 K and 0.45 GPa with a simulation at 400 K, the pure temperature
dependence is identified.
From the spectra of Fig.~\ref{DOSSM}, it is found that the explicit anharmonicity dominates
the softening and broadening of the phonon spectra. 
All features in the phonon spectra with the QHA 
showed little change with temperature, except for small stiffenings at high energies. 
%Of course,  the thermal lifetime broadening of phonon linewidths cannot be  
%predicted with the QHA because it assumes noninteracting normal modes.

\begin{figure}[b]
\includegraphics[width=0.8\columnwidth]{DOS1SM}
\caption{Neutron weighted phonon DOS of Ag$_2$O with the cuprite structure from MD simulations. 
The green spectrum is the MD simulated phonon DOS at 40\,K and 0.45\,GPa. Vertical dashed lines are aligned to the major peak centers at 40\,K and labeled by numbers. 
The incident energy was 100\,meV for panel (a), and 30\,meV for panel (b). The spectra were convoluted with the resolution function characteristic of ARCS for the different energies of the incident neutron beam. }
\label{DOSSM}
\end{figure} 
  




\bibliography{ref} 
\end{document}