| 83 |
|
traditional methods developed for homogeneous systems. |
| 84 |
|
|
| 85 |
|
Experimentally, various interfaces have been investigated for their |
| 86 |
< |
thermal conductance. Wang {\it et al.} studied heat transport through |
| 86 |
> |
thermal conductance. Cahill and coworkers studied nanoscale thermal |
| 87 |
> |
transport from metal nanoparticle/fluid interfaces, to epitaxial |
| 88 |
> |
TiN/single crystal oxides interfaces, to hydrophilic and hydrophobic |
| 89 |
> |
interfaces between water and solids with different self-assembled |
| 90 |
> |
monolayers.\cite{Wilson:2002uq,PhysRevB.67.054302,doi:10.1021/jp048375k,PhysRevLett.96.186101} |
| 91 |
> |
Wang {\it et al.} studied heat transport through |
| 92 |
|
long-chain hydrocarbon monolayers on gold substrate at individual |
| 93 |
|
molecular level,\cite{Wang10082007} Schmidt {\it et al.} studied the |
| 94 |
|
role of CTAB on thermal transport between gold nanorods and |
| 111 |
|
measurements for heat conductance of interfaces between the capping |
| 112 |
|
monolayer on Au and a solvent phase have yet to be studied with their |
| 113 |
|
approach. The comparatively low thermal flux through interfaces is |
| 114 |
< |
difficult to measure with Equilibrium MD or forward NEMD simulation |
| 114 |
> |
difficult to measure with Equilibrium |
| 115 |
> |
MD\cite{doi:10.1080/0026897031000068578} or forward NEMD simulation |
| 116 |
|
methods. Therefore, the Reverse NEMD (RNEMD) |
| 117 |
|
methods\cite{MullerPlathe:1997xw,kuang:164101} would have the |
| 118 |
|
advantage of applying this difficult to measure flux (while measuring |
| 200 |
|
discrete on the two sides of the interface. $G$ can be calculated |
| 201 |
|
using the applied thermal flux $J$ and the maximum temperature |
| 202 |
|
difference measured along the thermal gradient max($\Delta T$), which |
| 203 |
< |
occurs at the Gibbs deviding surface (Figure \ref{demoPic}): |
| 203 |
> |
occurs at the Gibbs deviding surface (Figure \ref{demoPic}). This is |
| 204 |
> |
known as the Kapitza conductance, which is the inverse of the Kapitza |
| 205 |
> |
resistance. |
| 206 |
|
\begin{equation} |
| 207 |
|
G=\frac{J}{\Delta T} |
| 208 |
|
\label{discreteG} |
| 358 |
|
these simulations. The chemically-distinct sites (a-e) are expanded |
| 359 |
|
in terms of constituent atoms for both United Atom (UA) and All Atom |
| 360 |
|
(AA) force fields. Most parameters are from |
| 361 |
< |
Refs. \protect\cite{TraPPE-UA.alkanes,TraPPE-UA.alkylbenzenes,TraPPE-UA.thiols} (UA) and \protect\cite{OPLSAA} (AA). Cross-interactions with the Au atoms are given in Table \ref{MnM}.} |
| 361 |
> |
Refs. \protect\cite{TraPPE-UA.alkanes,TraPPE-UA.alkylbenzenes,TraPPE-UA.thiols} |
| 362 |
> |
(UA) and \protect\cite{OPLSAA} (AA). Cross-interactions with the Au |
| 363 |
> |
atoms are given in Table \ref{MnM}.} |
| 364 |
|
\label{demoMol} |
| 365 |
|
\end{figure} |
| 366 |
|
|
| 495 |
|
A series of different initial conditions with a range of surface |
| 496 |
|
coverages was prepared and solvated with various with both of the |
| 497 |
|
solvent molecules. These systems were then equilibrated and their |
| 498 |
< |
interfacial thermal conductivity was measured with our NIVS |
| 498 |
> |
interfacial thermal conductivity was measured with the NIVS |
| 499 |
|
algorithm. Figure \ref{coverage} demonstrates the trend of conductance |
| 500 |
|
with respect to surface coverage. |
| 501 |
|
|
| 507 |
|
\label{coverage} |
| 508 |
|
\end{figure} |
| 509 |
|
|
| 510 |
+ |
In partially covered surfaces, the derivative definition for $G^\prime$ (Eq. \ref{derivativeG}) becomes difficult to apply, as the location of maximum change of $\lambda$ becomes washed out. The discrete definition (Eq. \ref{discreteG}) is easier to apply, as the Gibbs dividing surface is still well-defined. Therefore, $G$ (not $G^\prime$) was used in this section. |
| 511 |
|
|
| 512 |
< |
In partially covered butanethiol on the Au(111) surface, the |
| 502 |
< |
derivative definition for $G^\prime$ (Eq. \ref{derivativeG}) becomes |
| 503 |
< |
difficult to apply, as the location of maximum change of $\lambda$ |
| 504 |
< |
becomes washed out. The discrete definition (Eq. \ref{discreteG}) is |
| 505 |
< |
easier to apply, as the Gibbs dividing surface is still |
| 506 |
< |
well-defined. Therefore, $G$ (not $G^\prime$) was used in this |
| 507 |
< |
section. |
| 512 |
> |
From Figure \ref{coverage}, one can see the significance of the presence of capping agents. When even a small fraction of the Au(111) surface sites are covered with butanethiols, the conductivity exhibits an enhancement by at least a factor of 3. Cappping agents are clearly playing a major role in thermal transport at metal / organic solvent surfaces. |
| 513 |
|
|
| 514 |
< |
From Figure \ref{coverage}, one can see the significance of the |
| 510 |
< |
presence of capping agents. When even a small fraction of the Au(111) |
| 511 |
< |
surface sites are covered with butanethiols, the conductivity exhibits |
| 512 |
< |
an enhancement by at least a factor of 3. This indicates the important |
| 513 |
< |
role cappping agents are playing for thermal transport at metal / |
| 514 |
< |
organic solvent surfaces. |
| 514 |
> |
We note a non-monotonic behavior in the interfacial conductance as a function of surface coverage. The maximum conductance (largest $G$) happens when the surfaces are about 75\% covered with butanethiol caps. The reason for this behavior is not entirely clear. One explanation is that incomplete butanethiol coverage allows small gaps between butanethiols to form. These gaps can be filled by transient solvent molecules. These solvent molecules couple very strongly with the hot capping agent molecules near the surface, and can then carry away (diffusively) the excess thermal energy from the surface. |
| 515 |
|
|
| 516 |
< |
We note a non-monotonic behavior in the interfacial conductance as a |
| 517 |
< |
function of surface coverage. The maximum conductance (largest $G$) |
| 518 |
< |
happens when the surfaces are about 75\% covered with butanethiol |
| 519 |
< |
caps. The reason for this behavior is not entirely clear. One |
| 520 |
< |
explanation is that incomplete butanethiol coverage allows small gaps |
| 521 |
< |
between butanethiols to form. These gaps can be filled by transient |
| 522 |
< |
solvent molecules. These solvent molecules couple very strongly with |
| 523 |
< |
the hot capping agent molecules near the surface, and can then carry |
| 524 |
< |
(diffusively) the excess thermal energy away from the surface. |
| 516 |
> |
There appears to be a competition between the conduction of the thermal energy away from the surface by the capping agents (enhanced by greater coverage) and the coupling of the capping agents with the solvent (enhanced by interdigitation at lower coverages). This competition would lead to the non-monotonic coverage behavior observed here. |
| 517 |
|
|
| 518 |
< |
There appears to be a competition between the conduction of the |
| 527 |
< |
thermal energy away from the surface by the capping agents (enhanced |
| 528 |
< |
by greater coverage) and the coupling of the capping agents with the |
| 529 |
< |
solvent (enhanced by physical contact at lower coverages). This |
| 530 |
< |
competition would lead to the non-monotonic coverage behavior observed |
| 531 |
< |
here. |
| 518 |
> |
Results for rigid body toluene solvent, as well as the UA hexane, are within the ranges expected from prior experimental work.\cite{Wilson:2002uq,cahill:793,PhysRevB.80.195406} This suggests that explicit hydrogen atoms might not be required for modeling thermal transport in these systems. C-H vibrational modes do not see significant excited state population at low temperatures, and are not likely to carry lower frequency excitations from the solid layer into the bulk liquid. |
| 519 |
|
|
| 520 |
< |
A comparison of the results obtained from the two different organic |
| 534 |
< |
solvents can also provide useful information of the interfacial |
| 535 |
< |
thermal transport process. The deuterated hexane (UA) results do not |
| 536 |
< |
appear to be substantially different from those of normal hexane (UA), |
| 537 |
< |
given that butanethiol (UA) is non-deuterated for both solvents. The |
| 538 |
< |
UA models, even though they have eliminated C-H vibrational overlap, |
| 539 |
< |
still have significant overlap in the infrared spectra. Because |
| 540 |
< |
differences in the infrared range do not seem to produce an observable |
| 541 |
< |
difference for the results of $G$ (Figure \ref{uahxnua}). |
| 520 |
> |
The toluene solvent does not exhibit the same behavior as hexane in that $G$ remains at approximately the same magnitude when the capping coverage increases from 25\% to 75\%. Toluene, as a rigid planar molecule, cannot occupy the relatively small gaps between the capping agents as easily as the chain-like {\it n}-hexane. The effect of solvent coupling to the capping agent is therefore weaker in toluene except at the very lowest coverage levels. This effect counters the coverage-dependent conduction of heat away from the metal surface, leading to a much flatter $G$ vs. coverage trend than is observed in {\it n}-hexane. |
| 521 |
|
|
| 543 |
– |
\begin{figure} |
| 544 |
– |
\includegraphics[width=\linewidth]{uahxnua} |
| 545 |
– |
\caption{Vibrational spectra obtained for normal (upper) and |
| 546 |
– |
deuterated (lower) hexane in Au-butanethiol/hexane |
| 547 |
– |
systems. Butanethiol spectra are shown as reference. Both hexane and |
| 548 |
– |
butanethiol were using United-Atom models.} |
| 549 |
– |
\label{uahxnua} |
| 550 |
– |
\end{figure} |
| 551 |
– |
|
| 552 |
– |
Furthermore, results for rigid body toluene solvent, as well as other |
| 553 |
– |
UA-hexane solvents, are reasonable within the general experimental |
| 554 |
– |
ranges\cite{Wilson:2002uq,cahill:793,PhysRevB.80.195406}. This |
| 555 |
– |
suggests that explicit hydrogen might not be a required factor for |
| 556 |
– |
modeling thermal transport phenomena of systems such as |
| 557 |
– |
Au-thiol/organic solvent. |
| 558 |
– |
|
| 559 |
– |
However, results for Au-butanethiol/toluene do not show an identical |
| 560 |
– |
trend with those for Au-butanethiol/hexane in that $G$ remains at |
| 561 |
– |
approximately the same magnitue when butanethiol coverage differs from |
| 562 |
– |
25\% to 75\%. This might be rooted in the molecule shape difference |
| 563 |
– |
for planar toluene and chain-like {\it n}-hexane. Due to this |
| 564 |
– |
difference, toluene molecules have more difficulty in occupying |
| 565 |
– |
relatively small gaps among capping agents when their coverage is not |
| 566 |
– |
too low. Therefore, the solvent-capping agent contact may keep |
| 567 |
– |
increasing until the capping agent coverage reaches a relatively low |
| 568 |
– |
level. This becomes an offset for decreasing butanethiol molecules on |
| 569 |
– |
its effect to the process of interfacial thermal transport. Thus, one |
| 570 |
– |
can see a plateau of $G$ vs. butanethiol coverage in our results. |
| 571 |
– |
|
| 522 |
|
\subsection{Effects due to Solvent \& Solvent Models} |
| 523 |
< |
In addition to UA solvent/capping agent models, AA models are included |
| 574 |
< |
in our simulations as well. Besides simulations of the same (UA or AA) |
| 575 |
< |
model for solvent and capping agent, different models can be applied |
| 576 |
< |
to different components. Furthermore, regardless of models chosen, |
| 577 |
< |
either the solvent or the capping agent can be deuterated, similar to |
| 578 |
< |
the previous section. Table \ref{modelTest} summarizes the results of |
| 579 |
< |
these studies. |
| 523 |
> |
In addition to UA solvent and capping agent models, AA models have also been included in our simulations. In most of this work, the same (UA or AA) model for solvent and capping agent was used, but it is also possible to utilize different models for different components. We have also included isotopic substitutions (Hydrogen to Deuterium) to decrease the explicit vibrational overlap between solvent and capping agent. Table \ref{modelTest} summarizes the results of these studies. |
| 524 |
|
|
| 525 |
|
\begin{table*} |
| 526 |
|
\begin{minipage}{\linewidth} |
| 527 |
|
\begin{center} |
| 528 |
|
|
| 529 |
< |
\caption{Computed interfacial thermal conductivity ($G$ and |
| 529 |
> |
\caption{Computed interfacial thermal conductance ($G$ and |
| 530 |
|
$G^\prime$) values for interfaces using various models for |
| 531 |
|
solvent and capping agent (or without capping agent) at |
| 532 |
|
$\langle T\rangle\sim$200K. (D stands for deuterated solvent |
| 533 |
|
or capping agent molecules; ``Avg.'' denotes results that are |
| 534 |
< |
averages of simulations under different $J_z$'s. Error |
| 535 |
< |
estimates indicated in parenthesis.)} |
| 534 |
> |
averages of simulations under different applied thermal flux values $(J_z)$. Error |
| 535 |
> |
estimates are indicated in parentheses.)} |
| 536 |
|
|
| 537 |
|
\begin{tabular}{llccc} |
| 538 |
|
\hline\hline |
| 567 |
|
\end{minipage} |
| 568 |
|
\end{table*} |
| 569 |
|
|
| 570 |
< |
To facilitate direct comparison, the same system with differnt models |
| 627 |
< |
for different components uses the same length scale for their |
| 628 |
< |
simulation cells. Without the presence of capping agent, using |
| 629 |
< |
different models for hexane yields similar results for both $G$ and |
| 630 |
< |
$G^\prime$, and these two definitions agree with eath other very |
| 631 |
< |
well. This indicates very weak interaction between the metal and the |
| 632 |
< |
solvent, and is a typical case for acoustic impedance mismatch between |
| 633 |
< |
these two phases. |
| 570 |
> |
To facilitate direct comparison between force fields, systems with the same capping agent and solvent were prepared with the same length scales for the simulation cells. |
| 571 |
|
|
| 572 |
< |
As for Au(111) surfaces completely covered by butanethiols, the choice |
| 636 |
< |
of models for capping agent and solvent could impact the measurement |
| 637 |
< |
of $G$ and $G^\prime$ quite significantly. For Au-butanethiol/hexane |
| 638 |
< |
interfaces, using AA model for both butanethiol and hexane yields |
| 639 |
< |
substantially higher conductivity values than using UA model for at |
| 640 |
< |
least one component of the solvent and capping agent, which exceeds |
| 641 |
< |
the general range of experimental measurement results. This is |
| 642 |
< |
probably due to the classically treated C-H vibrations in the AA |
| 643 |
< |
model, which should not be appreciably populated at normal |
| 644 |
< |
temperatures. In comparison, once either the hexanes or the |
| 645 |
< |
butanethiols are deuterated, one can see a significantly lower $G$ and |
| 646 |
< |
$G^\prime$. In either of these cases, the C-H(D) vibrational overlap |
| 647 |
< |
between the solvent and the capping agent is removed (Figure |
| 648 |
< |
\ref{aahxntln}). Conclusively, the improperly treated C-H vibration in |
| 649 |
< |
the AA model produced over-predicted results accordingly. Compared to |
| 650 |
< |
the AA model, the UA model yields more reasonable results with higher |
| 651 |
< |
computational efficiency. |
| 572 |
> |
On bare metal / solvent surfaces, different force field models for hexane yield similar results for both $G$ and $G^\prime$, and these two definitions agree with each other very well. This is primarily an indicator of weak interactions between the metal and the solvent, and is a typical case for acoustic impedance mismatch between these two phases. |
| 573 |
|
|
| 574 |
+ |
For the fully-covered surfaces, the choice of force field for the capping agent and solvent has a large impact on the calulated values of $G$ and $G^\prime$. The AA thiol to AA solvent conductivities are much larger than their UA to UA counterparts, and these values exceed the experimental estimates by a large measure. The AA force field allows significant energy to go into C-H (or C-D) stretching modes, and since these modes are high frequency, this non-quantum behavior is likely responsible for the overestimate of the conductivity. |
| 575 |
+ |
|
| 576 |
+ |
The similarity in the vibrational modes available to solvent and capping agent can be reduced by deuterating one of the two components. Once either the hexanes or the butanethiols are deuterated, one can see a significantly lower $G$ and $G^\prime$ (Figure \ref{aahxntln}). Compared to the AA model, the UA model yields more reasonable conductivity values with much higher computational efficiency. |
| 577 |
+ |
|
| 578 |
|
\begin{figure} |
| 579 |
|
\includegraphics[width=\linewidth]{aahxntln} |
| 580 |
< |
\caption{Spectra obtained for All-Atom model Au-butanethil/solvent |
| 580 |
> |
\caption{Spectra obtained for all-atom (AA) Au / butanethiol / solvent |
| 581 |
|
systems. When butanethiol is deuterated (lower left), its |
| 582 |
< |
vibrational overlap with hexane would decrease significantly, |
| 658 |
< |
compared with normal butanethiol (upper left). However, this |
| 659 |
< |
dramatic change does not apply to toluene as much (right).} |
| 582 |
> |
vibrational overlap with hexane decreases significantly. Since aromatic molecules and the butanethiol are vibrationally dissimilar, the change is not as dramatic when toluene is the solvent (right).} |
| 583 |
|
\label{aahxntln} |
| 584 |
|
\end{figure} |
| 585 |
|
|
| 586 |
< |
However, for Au-butanethiol/toluene interfaces, having the AA |
| 664 |
< |
butanethiol deuterated did not yield a significant change in the |
| 665 |
< |
measurement results. Compared to the C-H vibrational overlap between |
| 666 |
< |
hexane and butanethiol, both of which have alkyl chains, that overlap |
| 667 |
< |
between toluene and butanethiol is not so significant and thus does |
| 668 |
< |
not have as much contribution to the heat exchange |
| 669 |
< |
process. Conversely, extra degrees of freedom such as the C-H |
| 670 |
< |
vibrations could yield higher heat exchange rate between these two |
| 671 |
< |
phases and result in a much higher conductivity. |
| 586 |
> |
For the Au / butanethiol / toluene interfaces, having the AA butanethiol deuterated did not yield a significant change in the measured conductance. Compared to the C-H vibrational overlap between hexane and butanethiol, both of which have alkyl chains, the overlap between toluene and butanethiol is not as significant and thus does not contribute as much to the heat exchange process. The presence of extra degrees of freedom in the AA force field for toluene yields higher heat exchange rates between the two phases and results in a much higher conductivity than in the UA force field. |
| 587 |
|
|
| 588 |
< |
Although the QSC model for Au is known to predict an overly low value |
| 589 |
< |
for bulk metal gold conductivity\cite{kuang:164101}, our computational |
| 675 |
< |
results for $G$ and $G^\prime$ do not seem to be affected by this |
| 676 |
< |
drawback of the model for metal. Instead, our results suggest that the |
| 677 |
< |
modeling of interfacial thermal transport behavior relies mainly on |
| 678 |
< |
the accuracy of the interaction descriptions between components |
| 679 |
< |
occupying the interfaces. |
| 588 |
> |
\subsubsection{Are electronic excitations in the metal important?} |
| 589 |
> |
Because they lack electronic excitations, the QSC and related embedded atom method (EAM) models for gold are known to predict unreasonably low values for bulk conductivity ($\lambda$).\cite{kuang:164101,ISI:000207079300006,Clancy:1992} If the conductance between the phases ($G$) is governed primarily by phonon excitation (and not electronic degrees of freedom), one would expect a classical model to capture most of the interfacial thermal conductance. Our results for $G$ and $G^\prime$ indicate that this is indeed the case, and suggest that the modeling of interfacial thermal transport depends primarily on the description of the interactions between the various components at the interface. When the metal is chemically capped, the primary barrier to thermal conductivity appears to be the interface between the capping agent and the surrounding solvent, so the excitations in the metal have little impact on the value of $G$. |
| 590 |
|
|
| 591 |
|
\subsection{Effects due to methodology and simulation parameters} |
| 592 |
|
|
| 593 |
< |
We have varied our protocol or other parameters of the simulations in |
| 684 |
< |
order to investigate how these factors would affect the measurement of |
| 685 |
< |
$G$'s. It turned out that while some of these parameters would not |
| 686 |
< |
affect the results substantially, some other changes to the |
| 687 |
< |
simulations would have a significant impact on the measurement |
| 688 |
< |
results. |
| 593 |
> |
START HERE |
| 594 |
|
|
| 595 |
< |
In some of our simulations, we allowed $L_x$ and $L_y$ to change |
| 691 |
< |
during equilibrating the liquid phase. Due to the stiffness of the |
| 692 |
< |
crystalline Au structure, $L_x$ and $L_y$ would not change noticeably |
| 693 |
< |
after equilibration. Although $L_z$ could fluctuate $\sim$1\% after a |
| 694 |
< |
system is fully equilibrated in the NPT ensemble, this fluctuation, as |
| 695 |
< |
well as those of $L_x$ and $L_y$ (which is significantly smaller), |
| 696 |
< |
would not be magnified on the calculated $G$'s, as shown in Table |
| 697 |
< |
\ref{AuThiolHexaneUA}. This insensivity to $L_i$ fluctuations allows |
| 698 |
< |
reliable measurement of $G$'s without the necessity of extremely |
| 699 |
< |
cautious equilibration process. |
| 595 |
> |
We have varied our protocol or other parameters of the simulations in order to investigate how these factors would affect the computation of $G$. |
| 596 |
|
|
| 597 |
+ |
We allowed $L_x$ and $L_y$ to change during equilibrating the liquid phase. Due to the stiffness of the crystalline Au structure, $L_x$ and $L_y$ would not change noticeably after equilibration. Although $L_z$ could fluctuate $\sim$1\% after a system is fully equilibrated in the NPT ensemble, this fluctuation, as well as those of $L_x$ and $L_y$ (which is significantly smaller), would not be magnified on the calculated $G$'s, as shown in Table \ref{AuThiolHexaneUA}. This insensivity to $L_i$ fluctuations allows reliable measurement of $G$'s without the necessity of extremely cautious equilibration process. |
| 598 |
+ |
|
| 599 |
|
As stated in our computational details, the spacing filled with |
| 600 |
|
solvent molecules can be chosen within a range. This allows some |
| 601 |
|
change of solvent molecule numbers for the same Au-butanethiol |
| 781 |
|
\label{fourier} |
| 782 |
|
\end{equation} |
| 783 |
|
|
| 886 |
– |
From Figure \ref{coverage}, one can see the significance of the |
| 887 |
– |
presence of capping agents. Even when a fraction of the Au(111) |
| 888 |
– |
surface sites are covered with butanethiols, the conductivity would |
| 889 |
– |
see an enhancement by at least a factor of 3. This indicates the |
| 890 |
– |
important role cappping agent is playing for thermal transport |
| 891 |
– |
phenomena on metal / organic solvent surfaces. |
| 784 |
|
|
| 785 |
< |
Interestingly, as one could observe from our results, the maximum |
| 894 |
< |
conductance enhancement (largest $G$) happens while the surfaces are |
| 895 |
< |
about 75\% covered with butanethiols. This again indicates that |
| 896 |
< |
solvent-capping agent contact has an important role of the thermal |
| 897 |
< |
transport process. Slightly lower butanethiol coverage allows small |
| 898 |
< |
gaps between butanethiols to form. And these gaps could be filled with |
| 899 |
< |
solvent molecules, which acts like ``heat conductors'' on the |
| 900 |
< |
surface. The higher degree of interaction between these solvent |
| 901 |
< |
molecules and capping agents increases the enhancement effect and thus |
| 902 |
< |
produces a higher $G$ than densely packed butanethiol arrays. However, |
| 903 |
< |
once this maximum conductance enhancement is reached, $G$ decreases |
| 904 |
< |
when butanethiol coverage continues to decrease. Each capping agent |
| 905 |
< |
molecule reaches its maximum capacity for thermal |
| 906 |
< |
conductance. Therefore, even higher solvent-capping agent contact |
| 907 |
< |
would not offset this effect. Eventually, when butanethiol coverage |
| 908 |
< |
continues to decrease, solvent-capping agent contact actually |
| 909 |
< |
decreases with the disappearing of butanethiol molecules. In this |
| 910 |
< |
case, $G$ decrease could not be offset but instead accelerated. [MAY NEED |
| 911 |
< |
SNAPSHOT SHOWING THE PHENOMENA / SLAB DENSITY ANALYSIS] |
| 912 |
< |
|
| 913 |
< |
A comparison of the results obtained from differenet organic solvents |
| 914 |
< |
can also provide useful information of the interfacial thermal |
| 915 |
< |
transport process. The deuterated hexane (UA) results do not appear to |
| 916 |
< |
be much different from those of normal hexane (UA), given that |
| 917 |
< |
butanethiol (UA) is non-deuterated for both solvents. These UA model |
| 918 |
< |
studies, even though eliminating C-H vibration samplings, still have |
| 919 |
< |
C-C vibrational frequencies different from each other. However, these |
| 920 |
< |
differences in the infrared range do not seem to produce an observable |
| 921 |
< |
difference for the results of $G$ (Figure \ref{uahxnua}). |
| 922 |
< |
|
| 923 |
< |
\begin{figure} |
| 924 |
< |
\includegraphics[width=\linewidth]{uahxnua} |
| 925 |
< |
\caption{Vibrational spectra obtained for normal (upper) and |
| 926 |
< |
deuterated (lower) hexane in Au-butanethiol/hexane |
| 927 |
< |
systems. Butanethiol spectra are shown as reference. Both hexane and |
| 928 |
< |
butanethiol were using United-Atom models.} |
| 929 |
< |
\label{uahxnua} |
| 930 |
< |
\end{figure} |
| 931 |
< |
|
| 932 |
< |
Furthermore, results for rigid body toluene solvent, as well as other |
| 933 |
< |
UA-hexane solvents, are reasonable within the general experimental |
| 934 |
< |
ranges\cite{Wilson:2002uq,cahill:793,PhysRevB.80.195406}. This |
| 935 |
< |
suggests that explicit hydrogen might not be a required factor for |
| 936 |
< |
modeling thermal transport phenomena of systems such as |
| 937 |
< |
Au-thiol/organic solvent. |
| 938 |
< |
|
| 939 |
< |
However, results for Au-butanethiol/toluene do not show an identical |
| 940 |
< |
trend with those for Au-butanethiol/hexane in that $G$ remains at |
| 941 |
< |
approximately the same magnitue when butanethiol coverage differs from |
| 942 |
< |
25\% to 75\%. This might be rooted in the molecule shape difference |
| 943 |
< |
for planar toluene and chain-like {\it n}-hexane. Due to this |
| 944 |
< |
difference, toluene molecules have more difficulty in occupying |
| 945 |
< |
relatively small gaps among capping agents when their coverage is not |
| 946 |
< |
too low. Therefore, the solvent-capping agent contact may keep |
| 947 |
< |
increasing until the capping agent coverage reaches a relatively low |
| 948 |
< |
level. This becomes an offset for decreasing butanethiol molecules on |
| 949 |
< |
its effect to the process of interfacial thermal transport. Thus, one |
| 950 |
< |
can see a plateau of $G$ vs. butanethiol coverage in our results. |
| 951 |
< |
|
| 952 |
< |
\subsection{Influence of Chosen Molecule Model on $G$} |
| 953 |
< |
In addition to UA solvent/capping agent models, AA models are included |
| 954 |
< |
in our simulations as well. Besides simulations of the same (UA or AA) |
| 955 |
< |
model for solvent and capping agent, different models can be applied |
| 956 |
< |
to different components. Furthermore, regardless of models chosen, |
| 957 |
< |
either the solvent or the capping agent can be deuterated, similar to |
| 958 |
< |
the previous section. Table \ref{modelTest} summarizes the results of |
| 959 |
< |
these studies. |
| 960 |
< |
|
| 961 |
< |
\begin{table*} |
| 962 |
< |
\begin{minipage}{\linewidth} |
| 963 |
< |
\begin{center} |
| 964 |
< |
|
| 965 |
< |
\caption{Computed interfacial thermal conductivity ($G$ and |
| 966 |
< |
$G^\prime$) values for interfaces using various models for |
| 967 |
< |
solvent and capping agent (or without capping agent) at |
| 968 |
< |
$\langle T\rangle\sim$200K. (D stands for deuterated solvent |
| 969 |
< |
or capping agent molecules; ``Avg.'' denotes results that are |
| 970 |
< |
averages of simulations under different $J_z$'s. Error |
| 971 |
< |
estimates indicated in parenthesis.)} |
| 972 |
< |
|
| 973 |
< |
\begin{tabular}{llccc} |
| 974 |
< |
\hline\hline |
| 975 |
< |
Butanethiol model & Solvent & $J_z$ & $G$ & $G^\prime$ \\ |
| 976 |
< |
(or bare surface) & model & (GW/m$^2$) & |
| 977 |
< |
\multicolumn{2}{c}{(MW/m$^2$/K)} \\ |
| 978 |
< |
\hline |
| 979 |
< |
UA & UA hexane & Avg. & 131(9) & 87(10) \\ |
| 980 |
< |
& UA hexane(D) & 1.95 & 153(5) & 136(13) \\ |
| 981 |
< |
& AA hexane & Avg. & 131(6) & 122(10) \\ |
| 982 |
< |
& UA toluene & 1.96 & 187(16) & 151(11) \\ |
| 983 |
< |
& AA toluene & 1.89 & 200(36) & 149(53) \\ |
| 984 |
< |
\hline |
| 985 |
< |
AA & UA hexane & 1.94 & 116(9) & 129(8) \\ |
| 986 |
< |
& AA hexane & Avg. & 442(14) & 356(31) \\ |
| 987 |
< |
& AA hexane(D) & 1.93 & 222(12) & 234(54) \\ |
| 988 |
< |
& UA toluene & 1.98 & 125(25) & 97(60) \\ |
| 989 |
< |
& AA toluene & 3.79 & 487(56) & 290(42) \\ |
| 990 |
< |
\hline |
| 991 |
< |
AA(D) & UA hexane & 1.94 & 158(25) & 172(4) \\ |
| 992 |
< |
& AA hexane & 1.92 & 243(29) & 191(11) \\ |
| 993 |
< |
& AA toluene & 1.93 & 364(36) & 322(67) \\ |
| 994 |
< |
\hline |
| 995 |
< |
bare & UA hexane & Avg. & 46.5(3.2) & 49.4(4.5) \\ |
| 996 |
< |
& UA hexane(D) & 0.98 & 43.9(4.6) & 43.0(2.0) \\ |
| 997 |
< |
& AA hexane & 0.96 & 31.0(1.4) & 29.4(1.3) \\ |
| 998 |
< |
& UA toluene & 1.99 & 70.1(1.3) & 65.8(0.5) \\ |
| 999 |
< |
\hline\hline |
| 1000 |
< |
\end{tabular} |
| 1001 |
< |
\label{modelTest} |
| 1002 |
< |
\end{center} |
| 1003 |
< |
\end{minipage} |
| 1004 |
< |
\end{table*} |
| 1005 |
< |
|
| 1006 |
< |
To facilitate direct comparison, the same system with differnt models |
| 1007 |
< |
for different components uses the same length scale for their |
| 1008 |
< |
simulation cells. Without the presence of capping agent, using |
| 1009 |
< |
different models for hexane yields similar results for both $G$ and |
| 1010 |
< |
$G^\prime$, and these two definitions agree with eath other very |
| 1011 |
< |
well. This indicates very weak interaction between the metal and the |
| 1012 |
< |
solvent, and is a typical case for acoustic impedance mismatch between |
| 1013 |
< |
these two phases. |
| 1014 |
< |
|
| 1015 |
< |
As for Au(111) surfaces completely covered by butanethiols, the choice |
| 1016 |
< |
of models for capping agent and solvent could impact the measurement |
| 1017 |
< |
of $G$ and $G^\prime$ quite significantly. For Au-butanethiol/hexane |
| 1018 |
< |
interfaces, using AA model for both butanethiol and hexane yields |
| 1019 |
< |
substantially higher conductivity values than using UA model for at |
| 1020 |
< |
least one component of the solvent and capping agent, which exceeds |
| 1021 |
< |
the general range of experimental measurement results. This is |
| 1022 |
< |
probably due to the classically treated C-H vibrations in the AA |
| 1023 |
< |
model, which should not be appreciably populated at normal |
| 1024 |
< |
temperatures. In comparison, once either the hexanes or the |
| 1025 |
< |
butanethiols are deuterated, one can see a significantly lower $G$ and |
| 1026 |
< |
$G^\prime$. In either of these cases, the C-H(D) vibrational overlap |
| 1027 |
< |
between the solvent and the capping agent is removed (Figure |
| 1028 |
< |
\ref{aahxntln}). Conclusively, the improperly treated C-H vibration in |
| 1029 |
< |
the AA model produced over-predicted results accordingly. Compared to |
| 1030 |
< |
the AA model, the UA model yields more reasonable results with higher |
| 1031 |
< |
computational efficiency. |
| 1032 |
< |
|
| 1033 |
< |
\begin{figure} |
| 1034 |
< |
\includegraphics[width=\linewidth]{aahxntln} |
| 1035 |
< |
\caption{Spectra obtained for All-Atom model Au-butanethil/solvent |
| 1036 |
< |
systems. When butanethiol is deuterated (lower left), its |
| 1037 |
< |
vibrational overlap with hexane would decrease significantly, |
| 1038 |
< |
compared with normal butanethiol (upper left). However, this |
| 1039 |
< |
dramatic change does not apply to toluene as much (right).} |
| 1040 |
< |
\label{aahxntln} |
| 1041 |
< |
\end{figure} |
| 1042 |
< |
|
| 1043 |
< |
However, for Au-butanethiol/toluene interfaces, having the AA |
| 1044 |
< |
butanethiol deuterated did not yield a significant change in the |
| 1045 |
< |
measurement results. Compared to the C-H vibrational overlap between |
| 1046 |
< |
hexane and butanethiol, both of which have alkyl chains, that overlap |
| 1047 |
< |
between toluene and butanethiol is not so significant and thus does |
| 1048 |
< |
not have as much contribution to the heat exchange |
| 1049 |
< |
process. Conversely, extra degrees of freedom such as the C-H |
| 1050 |
< |
vibrations could yield higher heat exchange rate between these two |
| 1051 |
< |
phases and result in a much higher conductivity. |
| 1052 |
< |
|
| 1053 |
< |
Although the QSC model for Au is known to predict an overly low value |
| 1054 |
< |
for bulk metal gold conductivity\cite{kuang:164101}, our computational |
| 1055 |
< |
results for $G$ and $G^\prime$ do not seem to be affected by this |
| 1056 |
< |
drawback of the model for metal. Instead, our results suggest that the |
| 1057 |
< |
modeling of interfacial thermal transport behavior relies mainly on |
| 1058 |
< |
the accuracy of the interaction descriptions between components |
| 1059 |
< |
occupying the interfaces. |
| 1060 |
< |
|
| 1061 |
< |
\subsection{Role of Capping Agent in Interfacial Thermal Conductance} |
| 785 |
> |
\subsubsection{The role of specific vibrations} |
| 786 |
|
The vibrational spectra for gold slabs in different environments are |
| 787 |
|
shown as in Figure \ref{specAu}. Regardless of the presence of |
| 788 |
|
solvent, the gold surfaces covered by butanethiol molecules, compared |
| 817 |
|
\label{specAu} |
| 818 |
|
\end{figure} |
| 819 |
|
|
| 820 |
< |
[MAY ADD COMPARISON OF AU SLAB WIDTHS BUT NOT MUCH TO TALK ABOUT...] |
| 820 |
> |
\subsubsection{Isotopic substitution and vibrational overlap} |
| 821 |
> |
A comparison of the results obtained from the two different organic |
| 822 |
> |
solvents can also provide useful information of the interfacial |
| 823 |
> |
thermal transport process. The deuterated hexane (UA) results do not |
| 824 |
> |
appear to be substantially different from those of normal hexane (UA), |
| 825 |
> |
given that butanethiol (UA) is non-deuterated for both solvents. The |
| 826 |
> |
UA models, even though they have eliminated C-H vibrational overlap, |
| 827 |
> |
still have significant overlap in the infrared spectra. Because |
| 828 |
> |
differences in the infrared range do not seem to produce an observable |
| 829 |
> |
difference for the results of $G$ (Figure \ref{uahxnua}). |
| 830 |
|
|
| 831 |
+ |
\begin{figure} |
| 832 |
+ |
\includegraphics[width=\linewidth]{uahxnua} |
| 833 |
+ |
\caption{Vibrational spectra obtained for normal (upper) and |
| 834 |
+ |
deuterated (lower) hexane in Au-butanethiol/hexane |
| 835 |
+ |
systems. Butanethiol spectra are shown as reference. Both hexane and |
| 836 |
+ |
butanethiol were using United-Atom models.} |
| 837 |
+ |
\label{uahxnua} |
| 838 |
+ |
\end{figure} |
| 839 |
+ |
|
| 840 |
|
\section{Conclusions} |
| 841 |
|
The NIVS algorithm we developed has been applied to simulations of |
| 842 |
|
Au-butanethiol surfaces with organic solvents. This algorithm allows |
