1 |
\documentclass[11pt]{article} |
2 |
\usepackage{amsmath} |
3 |
\usepackage{amssymb} |
4 |
\usepackage{setspace} |
5 |
\usepackage{endfloat} |
6 |
\usepackage{caption} |
7 |
%\usepackage{tabularx} |
8 |
\usepackage{graphicx} |
9 |
\usepackage{multirow} |
10 |
%\usepackage{booktabs} |
11 |
%\usepackage{bibentry} |
12 |
%\usepackage{mathrsfs} |
13 |
%\usepackage[ref]{overcite} |
14 |
\usepackage[square, comma, sort&compress]{natbib} |
15 |
\usepackage{url} |
16 |
\pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm |
17 |
\evensidemargin 0.0cm \topmargin -21pt \headsep 10pt \textheight |
18 |
9.0in \textwidth 6.5in \brokenpenalty=10000 |
19 |
|
20 |
% double space list of tables and figures |
21 |
\AtBeginDelayedFloats{\renewcommand{\baselinestretch}{1.66}} |
22 |
\setlength{\abovecaptionskip}{20 pt} |
23 |
\setlength{\belowcaptionskip}{30 pt} |
24 |
|
25 |
%\renewcommand\citemid{\ } % no comma in optional referenc note |
26 |
\bibpunct{[}{]}{,}{s}{}{;} |
27 |
\bibliographystyle{aip} |
28 |
|
29 |
\begin{document} |
30 |
|
31 |
\title{Simulating interfacial thermal conductance at metal-solvent |
32 |
interfaces: the role of chemical capping agents} |
33 |
|
34 |
\author{Shenyu Kuang and J. Daniel |
35 |
Gezelter\footnote{Corresponding author. \ Electronic mail: gezelter@nd.edu} \\ |
36 |
Department of Chemistry and Biochemistry,\\ |
37 |
University of Notre Dame\\ |
38 |
Notre Dame, Indiana 46556} |
39 |
|
40 |
\date{\today} |
41 |
|
42 |
\maketitle |
43 |
|
44 |
\begin{doublespace} |
45 |
|
46 |
\begin{abstract} |
47 |
The abstract version 2 |
48 |
\end{abstract} |
49 |
|
50 |
\newpage |
51 |
|
52 |
%\narrowtext |
53 |
|
54 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
55 |
% BODY OF TEXT |
56 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
57 |
|
58 |
\section{Introduction} |
59 |
|
60 |
The intro. |
61 |
|
62 |
\section{Methodology} |
63 |
\subsection{Algorithm} |
64 |
There have been many algorithms for computing thermal conductivity |
65 |
using molecular dynamics simulations. However, interfacial conductance |
66 |
is at least an order of magnitude smaller. This would make the |
67 |
calculation even more difficult for those slowly-converging |
68 |
equilibrium methods. Imposed-flux non-equilibrium |
69 |
methods\cite{MullerPlathe:1997xw} have the flux set {\it a priori} and |
70 |
the response of temperature or momentum gradients are easier to |
71 |
measure than the flux, if unknown, and thus, is a preferable way to |
72 |
the forward NEMD methods. Although the momentum swapping approach for |
73 |
flux-imposing can be used for exchanging energy between particles of |
74 |
different identity, the kinetic energy transfer efficiency is affected |
75 |
by the mass difference between the particles, which limits its |
76 |
application on heterogeneous interfacial systems. |
77 |
|
78 |
The non-isotropic velocity scaling (NIVS)\cite{kuang:164101} approach in |
79 |
non-equilibrium MD simulations is able to impose relatively large |
80 |
kinetic energy flux without obvious perturbation to the velocity |
81 |
distribution of the simulated systems. Furthermore, this approach has |
82 |
the advantage in heterogeneous interfaces in that kinetic energy flux |
83 |
can be applied between regions of particles of arbitary identity, and |
84 |
the flux quantity is not restricted by particle mass difference. |
85 |
|
86 |
The NIVS algorithm scales the velocity vectors in two separate regions |
87 |
of a simulation system with respective diagonal scaling matricies. To |
88 |
determine these scaling factors in the matricies, a set of equations |
89 |
including linear momentum conservation and kinetic energy conservation |
90 |
constraints and target momentum/energy flux satisfaction is |
91 |
solved. With the scaling operation applied to the system in a set |
92 |
frequency, corresponding momentum/temperature gradients can be built, |
93 |
which can be used for computing transportation properties and other |
94 |
applications related to momentum/temperature gradients. The NIVS |
95 |
algorithm conserves momenta and energy and does not depend on an |
96 |
external thermostat. |
97 |
|
98 |
(wondering how much detail of algorithm should be put here...) |
99 |
|
100 |
\subsection{Simulation Parameters} |
101 |
Our simulation systems consists of metal gold lattice slab solvated by |
102 |
organic solvents. In order to study the role of capping agents in |
103 |
interfacial thermal conductance, butanethiol is chosen to cover gold |
104 |
surfaces in comparison to no capping agent present. |
105 |
|
106 |
The Au-Au interactions in metal lattice slab is described by the |
107 |
quantum Sutton-Chen (QSC) formulation.\cite{PhysRevB.59.3527} The QSC |
108 |
potentials include zero-point quantum corrections and are |
109 |
reparametrized for accurate surface energies compared to the |
110 |
Sutton-Chen potentials\cite{Chen90}. |
111 |
|
112 |
Straight chain {\it n}-hexane and aromatic toluene are respectively |
113 |
used as solvents. For hexane, both United-Atom\cite{TraPPE-UA.alkanes} |
114 |
and All-Atom\cite{OPLSAA} force fields are used for comparison; for |
115 |
toluene, United-Atom\cite{TraPPE-UA.alkylbenzenes} force fields are |
116 |
used with rigid body constraints applied. (maybe needs more details |
117 |
about rigid body) |
118 |
|
119 |
Buatnethiol molecules are used as capping agent for some of our |
120 |
simulations. United-Atom\cite{TraPPE-UA.thiols} and All-Atom models |
121 |
are respectively used corresponding to the force field type of |
122 |
solvent. |
123 |
|
124 |
To describe the interactions between metal Au and non-metal capping |
125 |
agent and solvent, we refer to Vlugt\cite{vlugt:cpc2007154} and derive |
126 |
other interactions which are not parametrized in their work. (can add |
127 |
hautman and klein's paper here and more discussion; need to put |
128 |
aromatic-metal interaction approximation here) |
129 |
|
130 |
\section{Computational Details} |
131 |
Our simulation systems consists of a lattice Au slab with the (111) |
132 |
surface perpendicular to the $z$-axis, and a solvent layer between the |
133 |
periodic Au slabs along the $z$-axis. To set up the interfacial |
134 |
system, the Au slab is first equilibrated without solvent under room |
135 |
pressure and a desired temperature. After the metal slab is |
136 |
equilibrated, United-Atom or All-Atom butanethiols are replicated on |
137 |
the Au surface, each occupying the (??) among three Au atoms, and is |
138 |
equilibrated under NVT ensemble. According to (CITATION), the maximal |
139 |
thiol capacity on Au surface is $1/3$ of the total number of surface |
140 |
Au atoms. |
141 |
|
142 |
\cite{packmol} |
143 |
|
144 |
\section{Acknowledgments} |
145 |
Support for this project was provided by the National Science |
146 |
Foundation under grant CHE-0848243. Computational time was provided by |
147 |
the Center for Research Computing (CRC) at the University of Notre |
148 |
Dame. \newpage |
149 |
|
150 |
\bibliography{interfacial} |
151 |
|
152 |
\end{doublespace} |
153 |
\end{document} |
154 |
|