Next Article in Journal
Biodegradation of Conventional and Emerging Pollutants
Next Article in Special Issue
Photocrosslinking of cDNA Display Molecules with Their Target Proteins as a New Strategy for Peptide Selection
Previous Article in Journal
Combination of High Hydrostatic Pressure and Ultrafiltration to Generate a New Emulsifying Ingredient from Egg Yolk
Previous Article in Special Issue
Fluorescence Differentiation of ATP-Related Multiple Enzymatic Activities in Synovial Fluid as a Marker of Calcium Pyrophosphate Deposition Disease Using Kyoto Green
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nanoscale Quantum Thermal Conductance at Water Interface: Green’s Function Approach Based on One-Dimensional Phonon Model

by
Toshihito Umegaki
* and
Shigenori Tanaka
*
Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(5), 1185; https://doi.org/10.3390/molecules25051185
Submission received: 18 January 2020 / Revised: 17 February 2020 / Accepted: 25 February 2020 / Published: 5 March 2020

Abstract

:
We have derived the fundamental formula of phonon transport in water for the evaluation of quantum thermal conductance by using a one-dimensional phonon model based on the nonequilibrium Green’s function method. In our model, phonons are excited as quantum waves from the left or right reservoir and propagate from left to right of H 2 O layer or vice versa. We have assumed these reservoirs as being of periodic structures, whereas we can also model the H 2 O sandwiched between these reservoirs as having aperiodic structures of liquid containing N water molecules. We have extracted the dispersion curves from the experimental absorption spectra of the OH stretching and intermolecular modes of water molecules, and calculated phonon transmission function and quantum thermal conductance. In addition, we have simplified the formulation of the transmission function by employing a case of one water molecule (N=1). From this calculation, we have obtained the characteristic that the transmission probability is almost unity at the frequency bands of acoustic and optical modes, and the transmission probability vanishes by the phonon attenuation reflecting the quantum tunnel effect outside the bands of these two modes. The classical limit of the thermal conductance calculated by our formula agreed with the literature value (order of 10 10 W/K) in high temperature regime (>300 K). The present approach is powerful enough to be applicable to molecular systems containing proteins as well, and to evaluate their thermal conductive characteristics.

Graphical Abstract

1. Introduction

In recent years, advances in the fabrication and characterization of nanoscale systems allow for a better understanding of the heat flow at the microscopic level [1]. Regarding the researches for inorganic materials, studies on carbon nanotube systems have provided a prototypical example for nanoscale thermal conduction experimentally [2,3,4,5,6] and theoretically [7,8]. Both theoretical and experimental progresses have thus been reported extensively in each of these studies [1]. As for other nanoscale systems, Lervik et al. [9] analyzed via classical molecular dynamics simulations the heat transfer through nanometer-scale interfaces consisting of n-decane (2–12 nm diameter) droplets in water. Tanaka et al. [10,11] studied molecular dynamics of water by microwave heating. In contrast to these classical-mechanics-based approaches, Fisher [12] showed that the Landauer transport formalism can be applied to the formulations of the quantum thermal conductance of heat flux carried by phonons between hot and cold reservoirs on the basis of nanoscale models of inorganic materials. Furthermore, using the Landauer formulation, Rego et al. [13] studied the quantized thermal conductance of dielectric quantum wires at low temperatures.
On the other hand, attention has recently been focused on the nanoscale thermal conduction in biological systems as well. Concerning experimental reports on endogenous heat production in single cells, Baffou et al. [14] argued that the experimentally observed temperature rise of Δ T (1–2 K) of a whole cell [15,16,17] would be much larger than the theoretical Δ T estimated from the relation Δ T = P / κ L , where κ is the thermal conductivity, L is the diameter of the heat source and P is the power. This equation is derived from the macroscopic heat diffusion equation for continuum [18]. Baffou et al. [19] estimated Δ T within a cell to be around 10 5 K. However, the 10 5 gap issue (the difference between the measured temperature increase of 1 K in a single cell using cellular thermometry as compared to the theoretically calculated increase estimate of 10 5 K) may disappear [15] if once the uncertainties of variable heat sources in stimulated cells (10 1 ), length scales (10 1 –10 2 ), and micro- and nanoscale thermal parameters (10 1 –10 2 ) are considered. A forthcoming challenge in the series of debates [14,15] is thus to establish a relevant theoretical framework for thermal conduction by setting up nanoscale models of biological cells containing liquid water and proteins. Regarding this issue, Lervik et al. [20] calculated the thermal conductivity at the protein–water interface in terms of classical molecular dynamics. Pandey and Leitner [21,22] quantum-mechanically evaluated the thermal energy transport through a trehalose layer between water and protein, and between gold, such as a gold nanoparticle, and its cellular environment.
In this paper, we aim at constructing a theoretical formulation to describe the quantum thermal conduction at nanoscale, paying special attention to the roles by water in biological contexts. We here derive the fundamental formula of phonon transport in water for the first step of evaluations of quantum thermal conductances in intracellular molecular environment. In the nanoscale space relevant to intracellular environment, we expect that the quantum characteristics of heat transfer may manifest themselves, which should be compared to the classical descriptions. The model of this study, which is given by a novel combination of the Landauer formulation and spectroscopic data, is powerful enough to be applicable to molecular systems containing proteins, and to quantitatively evaluate the thermal conductive characteristics in realistic systems, since experimental spectroscopic features are taken into account in the model. In the present formulation, we have applied the Green’s function method (GFM) to water with a one-dimensional phonon model. This model then takes into account not only intramolecular vibration modes but also intermolecular ones. We have extracted the dispersion curves from experimental absorption spectrum of water and evaluated quantitatively quantum heat conduction characteristics such as the phonon transmission function and the thermal conductance. Based on the GFM theoretical framework and the extracted parameters, we have calculated the quantum thermal conductance and compared it with the literature values quantitatively. We show detailed mathematical formulations based on equilibrium Green’s functions, non-equilibrium Green’s functions and phonon transmission function in the Appendix A, and describe simplified formulas used in the calculations in the following sections of main text.

2. Extraction of Phonon Dispersion Curves of Liquid Water from Experimental Results

Figure 1 illustrates the definitions of displacements, masses, and spring constants of the elements H 2 and O in the present one-dimensional phonon model. Figure 1 gives a background of the theoretical formulations. The j-th unit cell of water molecule with length l j consists of two elements, H 2 and O, blue and yellow colored circles, which have displacements u j , 1 and u j , 2 , masses M j , 1 and M j , 2 , and spring constants K j , 1 and K j , 2 , respectively. By applying the Bloch theorem u j + 1 , n = exp ( i k l j ) u j , n to Equations (A1a)–(A2b) in Appendix A, we can obtain a phonon dispersion relation of liquid water between the frequency ω and the dimensionless wavenumber:
k l j = arccos 2 M j , 1 1 M j , 2 1 K j , 1 K j , 2 ( M j , 1 1 + M j , 2 1 ) ( K j , 1 + K j , 2 ) ω 2 + ω 4 2 M j , 1 1 M j , 2 1 K j , 1 K j , 2 ,
where j and n in K j , n and M j , n are indices of unit cells and atomic components, respectively. In the model of Figure 2, j and n take values from ( N L + 1 ) to ( N + N R ) and 1, 2, respectively, where N L , N R , and N represent the number of unit cells in the left and right reservoirs and the water, respectively, and n = 1 , 2 mean two-hydrogen (H 2 ) and one-oxygen (O), respectively. The phonon dispersion relation (1) has corresponding dispersion curves with real k branches of acoustic and optical modes. Hereafter, assuming that M j , 1 and M j , 2 do not depend on j, we have approximated them by M 1 and M 2 , respectively, and similarly, K j , 1 and K j , 2 were also approximated by K 1 and K 2 , respectively. If we assume the acoustic and optical modes have angular frequency bands of 0 ω ω r and ω q ω ω p , respectively, we can obtain equations for K 1 , K 2 and ω p as follows:
K 1 = 1 2 M 1 M 2 ( ω q 2 + ω r 2 ) M 1 + M 2 + M 1 M 2 ( ω q 2 + ω r 2 ) M 1 + M 2 2 M 1 M 2 ω q 2 ω r 2 ,
K 2 = 1 2 M 1 M 2 ( ω q 2 + ω r 2 ) M 1 + M 2 M 1 M 2 ( ω q 2 + ω r 2 ) M 1 + M 2 2 M 1 M 2 ω q 2 ω r 2 ,
ω p = M 1 1 K 1 + M 1 1 K 2 + M 2 1 K 1 + M 2 1 K 2 .
As illustrated in Figure 2, we have assumed the model in which phonons are excited as quantum waves from the left or right reservoir and propagate from left to right of water layers or vice versa. These reservoirs have periodic structures of solid in the red and blue areas of Figure 2, whereas water may have an aperiodic structure of liquid in the rainbow colored area of Figure 2; the reason for employing the aperiodic structure is to generally describe disordered liquid water. We have assumed a one-dimensional phonon model with intervals of water molecule as l j in the j-th unit cell, where j = 1 , 2 , N as shown in Figure 2. Each unit cell of water is modelled to consist of two components, H 2 and O. The present model takes into account not only intramolecular vibration modes but also intermolecular ones. Phonon has two kinds of modes, acoustic and optical modes. The components of water molecules at acoustic and optical modes vibrate with the same and inverse directions, respectively. We have extracted the dispersion curves from the experimental absorption spectrum of H 2 O in the liquid phase [24] by the following procedure:
  • Extracting wave numbers k b and k e at yellow and red broken lines of absorbance maxima in the experimental absorption spectrum of water in the liquid phase shown in Figure 3.
  • Calculating angular frequencies ω q and ω r as 2 π k b c and 2 π k e c , respectively, with c being the light velocity, so that the one-dimensional phonon model is consistent with the experimental spectrum.
  • Calculating spring constants K 1 and K 2 by using Equations (2a) and (2b) above.
  • Calculating angular frequency ω p by using Equation (2c) above.
  • Calculating the wave number k of phonon by using the dispersion relation Equation (1) above.
  • Obtaining dispersion curves in the upper part of Figure 4 below with the frequency on the abscissa and the normalized wave number k l j / π on the ordinate.
Figure 5 shows liquid-phase water motions [25,26] of (a) anti-symmetric stretching, (b) symmetric stretching, (c) bending, (d) libration, and (e) intermolecular vibration, respectively. In step 1, we determined wave numbers k b and k e as modes (b) and (e) from the maximum value of the experimental absorption spectrum shown in Figure 3, where we selected k b as maximum value of intramolecular vibration in the experimental spectrum because the absorbance of mode (b) was higher than that of mode (a) [27]. It is also noted that some other modes such as those associated with hydrogen bond stretch may be contained in the “libration” region (d) in Figure 3, while they are supposed to play an insignificant role in the present simplified model. This procedure is applicable not only to liquid water, but also to water in a wide range of solid, liquid, and gas temperatures [28]. Table 1 shows the wavenumbers and angular frequencies of various molecular vibrational modes of liquid water, where k a , k b , k c , k d , and k e are wavenumbers of liquid-phase water motions [25,26] of (a)–(e), respectively. Because of small difference between k a and k b , we have calculated k a as k a k b + k b , where k a and k b are wave numbers of isolated water molecule 3756 cm 1 and 3657 cm 1 , respectively [27].

3. Formulation of the Phonon Transmission Function through H 2 O

To obtain the physical outlook without losing physical essence, we have derived a simple expression for one H 2 O molecule on the basis of the formula derived from the equation of the phonon transmission function for a region consisting of N water molecules, described in Appendix A.
In this section, we derive equations of the phonon transmission function t r ( T ) for H 2 and O components of one water molecule, that is, the N=1 phonon model of Figure 6 for simplicity, and perform test calculations using the formula based on the N=1 water layer model without losing physical essence.
We set N = 1 in Equations (A12), (A14c), (A14d), (A15a), (A15b), (A23a), (A23b), (A24), and obtained the following equations:
t r ( T ) = Γ L , 11 Γ R , 22 G D , 12 2 ,
G D , 12 = γ 0 , 2 ( ε + i η ε 1 , 1 γ 1 , 1 γ 0 , 4 g L , 22 ) ( ε + i η ε 1 , 2 γ N , 4 γ N + 1 , 1 g R , 11 ) γ 1 , 2 γ 1 , 3 ,
Γ L , 11 = 2 γ 1 , 1 γ 0 , 4 Im ( g L , 22 ) ,
Γ R , 22 = 2 γ N , 4 γ N + 1 , 1 Im ( g R , 11 ) ,
g L , 22 = ε + i η ε 1 , 1 + ε 1 , 1 ε 1 exp ( i ( k L + k L ) l 0 ) ε + i η ε 0 , 1 + ε 0 , 1 ε 1 exp ( i ( k L + k L ) l 0 ) ( ε + i η ε 0 , 2 ) γ 0 , 2 γ 0 , 2 + γ 0 , 1 exp ( i k L l 0 ) 1 exp ( i ( k L + k L ) l 0 ) γ 0 , 3 ,
g R , 11 = ε + i η ε N , 2 + ε N , 2 ε 1 exp ( i ( k R + k R ) l N ) ε + i η ε N , 2 + ε N , 2 ε 1 exp ( i ( k R + k R ) l N ) ( ε + i η ε N , 1 ) γ N , 3 γ N + 1 , 3 + γ N + 1 , 4 exp ( i k R + l N ) 1 exp ( i ( k R + k R ) l N ) γ N , 2 ,
where tr means trace of matrix; Γ L , j n , Γ R , j n , G D , j n , g L , j n and g R , j n (j, n = 1, 2) are j n elements of matrices Γ L , Γ R , G D , g L and g R , respectively; ε , ε j , n and γ j , n (n = 1–4) represent ω 2 , squared on-site energy and squared hopping energy, respectively, as defined in Equations (A2c)–(A2i) and (A6a)–(A6c) in Appendix A; Im ( Z ) is the imaginary part of complex number Z, and η is the imaginary part of the energy, which is a positive infinitesimal. T represents the transmission matrix:
T = Γ L G D Γ R G D ,
as defined in Equation (A24) in Appendix A. The complex dispersion curves shown in the upper part of Figure 4 has three modes of ( I ) acoustic wave, ( I I ) attenuation and ( I I I ) optical wave as shown in Table 2. In this study, we have employed an interval of water molecule l 1 as 0.319 nm [29] in Figure 2. We can calculate the transmission function by Equation (3a), and ( I ) , ( I I ) and ( I I I ) modes have the transmission probabilities of approximately 1, 0, 1 over angular frequency ranges of 0 < ω < ω r , ω r < ω < ω q , and ω q < ω < ω p , respectively, as shown in Table 2. Then we can depict the frequency characteristics of the transmission probability as shown in the lower part of Figure 4. The transmission probability takes not only discrete values of 0 or 1, but continuous values between 0 and 1 near the boundary angular frequencies of ω r and ω q .

4. Thermal Conductance in Water

4.1. Validation of Calculated Thermal Conductance

Since phonons are bosons, they accord to the Bose-Einstein distribution as follows:
N j ( ω ) = exp ( ħ ω / k B θ j ) 1 1 ,
where , ω , k B , θ j are reduced Planck constant, angular frequency, Boltzmann constant, and temperature in the jth unit cell of water molecules in the one-dimensional phonon model of Figure 2. We can thus formulate the Bose-Einstein distribution function at the unit cell j in one-dimensional phonon model of reservoir-H 2 O-reservoir structure shown in Figure 2. Moreover, we can define the heat flow J by the following equation [1,23,31]:
J = 0 d ω ħ ω 2 π k = 1 2 j = 1 N T j , k ; j , k N 0 ( ω ) N N + 1 ( ω ) ,
where the functions N 0 ( ω ) and N N + 1 ( ω ) are the Bose-Einstein distributions at left and right reservoirs, respectively. Then, thermal conductance σ is defined by an equation that divides the heat flow J by the temperature difference θ 0 θ N + 1 as follows:
σ = J θ 0 θ N + 1 .
Figure 7 shows temperature dependencies of thermal conductances by this calculation with Equation (7). The temperature dependences show fair agreement with the previous results [32] for homogeneous and heterogeneous atomic chains, where the mass of a ‘Device’ atom in the homogeneous case is 4.6 × 10 26 kg (Si atom assumed), and the masses of ‘Device’ atoms in the two heterogeneous cases are 9.2 and 2.3 × 10 26 kg, respectively. It is observed that the agreement is in particular better when the temperature difference θ 0 θ N + 1 is smaller. It is noted that the conductance curves in Figure 7 are not so sensitive to the atomic masses employed in the calculations.

4.2. Classical Limit of Thermal Conductance

Because the higher-order terms ( ħ ω / k B θ j ) 2 / 2 ! + ( ħ ω / k B θ j ) 3 / 3 ! + can be neglected in the classical limit ħ ω / k B θ j 1 , the Bose-Einstein distribution function N j ( ω ) of Equation (5) takes the following limiting form:
N j ( ω ) = 1 + ( ħ ω / k B θ j ) / 1 ! + ( ħ ω / k B θ j ) 2 / 2 ! + 1 1 k B θ j ħ ω .
For example, in a case of temperature θ j =300K and angular frequency range ω ω r , we have ħ ω / k B θ j 0 . 912 . In this case we may use the approximation of Equation (8). N 0 ( ω ) N N + 1 ( ω ) of Equation (6) can then be approximated to ( k B / ħ ω ) ( θ 0 θ N + 1 ) . The phonon transmission probability k = 1 2 j = 1 N T j , k ; j , k can be approximated to unity in the bands of 0 < ω < ω r and ω q < ω < ω p , and zero for the other bands of ω according to the lower part of Figure 4. From these two approximations, the integral with respect to angular frequency of Equation (6) can be approximated by the following equation,
J = k B 2 π ω r ω q + ω p ( θ 0 θ N + 1 ) k B ω r 2 π ( θ 0 θ N + 1 ) ,
where we have transformed the 3rd expression in Equation (9) by using ω q + ω p ω r . Equation (9) can be used for the quantum heat flux due to phonon propagation at the molecular size in the cell, which may be formulated according to the actual size of cell biology. In the classical limit, thermal conductance σ is thus expressed by Equations (7) and (9) as follows:
σ = k B ω r 2 π .
By using Equation (10), we can evaluate the thermal conductance in the classical limit as σ = 7 . 87 × 10 11 W/K, which agrees well with the values in the high-temperature limit in Figure 7 showing the order of 10 10 W/K; this evaluation is also similar to the estimated value of about 100 MW/K/m 2 by Lervik et al. [20] with inclusion of the effective cross section (1 nm) 2 . The quantum thermal conductance evaluated in the present study thus agrees with the thermal conductance obtained via classical molecular dynamics simulation in the high-temperature, classical limit. One of significant results in the present study is that the magnitude of thermal conductance σ in the classical limit calculated from ω r is similar to that by classical molecular dynamics for protein–water system, whereas the former would give somewhat lower estimate due to the neglect of the contributions from other modes.
Baffou et al. argued [19] that the observed temperature rise of Δ T (1-2 K) of a whole cell [16,17] would be much larger than the theoretical Δ T estimated from the relation Δ T = P / κ L , where L is the diameter of the heat source, P is power, and κ is thermal conductivity estimated to be around 1 W/K/m [19]. They estimated Δ T within a cell to be around 10 5 K, assuming that the heat source size is 10 nm and that the average heat production per hour based on glucose production is 100 pW (360 nJ/h). However, one may expect that once three uncertainties below are supposed, the 10 5 gap issue may disappear [15]:
  • variable heat sources in stimulated cells (order 10 1 );
  • length scales (order 10 1 –10 2 );
  • micro- and nanoscale thermal parameters (order 10 1 –10 2 ).
Corresponding to these issues, we can make use of analytical model dicussed above. In order to assess these issues, in this paper, we constructed the following models in the framework of nonequilibrium Green’s function method:
  • reservoir model representing heat and bath reservoirs;
  • one-dimensional phonon model with water molecules sandwiched between the reservoirs;
  • phonon transports at atomic scale in water.
Thus Equation (7) can provide quantitative suggestions for the discussions on intracellular temperature distributions [14,15,19,33]. At physiological temperature, our quantum-mechanical estimation for the thermal conductance above is consistent with the values evaluated through classical molecular dynamics simulations for protein–water systems [20]. Then, if we would employ the evaluations of thermal conductivity κ by Lervik et al. [20] as well, the values of κ would fall around 0.1–0.3 W/K/m, which could make the discussions [14,15,19,33] more quantitative in nanoscale intracellular conditions. Of course, further theoretical investigations would be required to resolve the gap issue.

5. Conclusions

We have derived the fundamental formula of phonon transport in water for the first step of evaluations of quantum thermal conductance by using a one-dimensional phonon model. We have calculated the dispersion curves reproducing OH stretching mode of water molecule and inter-molecular mode, and reproduced experimental values of the absorption spectra by fitting the parameters of the phonon modes such as ω r , ω q and so on. We also have formulated the phonon transmission function of N water molecules connected to the left and right heat baths by using the non-equilibrium Green’s function method. We have employed a model in which phonons are excited as quantum waves from the left or right reservoir and propagate from left to right of H 2 O layer or vice versa. We have modelled these reservoirs as having periodic structures, whereas we can regard the water layer part as aperiodic structure of liquid. In this formulation, we have applied the Green’s function method to H 2 O with a one-dimensional phonon model. The model of this formulation is powerful enough to be applicable to molecular system containing proteins, and to evaluate thermal conductive characteristics thereof.
We have extracted dispersion curves from experimental absorption spectrum of H 2 O and quantitatively evaluated quantum heat conduction characteristics such as the phonon transmission function and the thermal conductance. In addition, we have simplified the formulation of the transmission function by employing one water molecule (N=1) in the aperiodic layer. From this calculation, we have obtained the characteristic that the transmission probability is almost unity in the frequency bands of acoustic and optical modes, whereas the transmission probability vanishes reflecting the quantum tunnel effect in the frequency band between these two modes. The classical limit of the thermal conductance calculated by our formula and the literature value [32] in high temperature regime (>300 K) have agreed with each other on the order of 10 10 W/K. The present model has also provided the calculated results for thermal conductance that are consistent with those evaluated via classical molecular dynamics simulations for protein–water systems [20] in the classical limit.
As future issues, in order to extract more information from the experimental results of water absorption spectra [34], we will expand our one-dimensional phonon model to higher-dimensional models including other vibration modes [35]. We expect this refinement of our model to provide more accurate reproduction of heat conduction characteristics. We are aiming at multi-scale heat conduction analysis by modeling not only water but also intracellular molecules such as proteins. We would like to apply the present method also to the analysis of heat conduction characteristics of biomolecules containing surrounding molecules such as flexible sugar chains [36], which would give more insights into the thermal conductive phenomena in intracellular crowding environments.

Author Contributions

T.U. and S.T. conceived and designed the theoretical models and formulations; T.U. performed the calculations; T.U. and S.T. analyzed the data; All authors contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JSPS KAHENHI Grant #17H06353 and 18K03825.

Acknowledgments

This work was supported by JSPS KAKENHI Grant #17H06353 and 18K03825.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
GFMGreen’s function method
Trad/sunit of tera radians per second
trtrace of matrix
Imimaginary part of complex number

Appendix A. Phonon Transport Equations in a One-Dimensional Chain Model of Water

We perform derivations in Appendix A.1: transport equations of phonon in water, Appendix A.2: dispersion equations in left and right reservoirs, Appendix A.3: non-equivribrium Green’s functions, and Appendix A.4: quantum thermal conductance in water in the following.

Appendix A.1. Derivation of Transport Equations of Phonon in Water

We can formulate phonon transport equations in a one-dimensional model of water in terms of atomic displacements u j , n ( t ) , j = N L + 1 , 1 , 0 , 1 , N + N R based on Figure 1 as follows:
M j , 1 d 2 u j , 1 ( t ) d t 2 = K j 1 , 2 u j , 1 ( t ) u j 1 , 2 ( t ) K j , 1 u j , 1 ( t ) u j , 2 ( t ) ,
M j , 2 d 2 u j , 2 ( t ) d t 2 = K j , 1 u j , 2 ( t ) u j , 1 ( t ) K j , 2 u j , 2 ( t ) u j + 1 , 1 ( t ) ,
where n = 1 , 2 in u j , n ( t ) mean two-hydrogen (H 2 ) and one-oxygen (O) atomic components, respectively. This model can mimic not only intramolecular vibration modes but also intermolecular ones. We can obtain time-independent equations of Fourier transforms of Equations (A1a) and (A1b) as follows.
γ j , 1 u ^ j 1 , 2 ( ω ) + ε j , 1 u ^ j , 1 γ j , 2 u ^ j , 2 ( ω ) = ε u ^ j , 1 ( ω ) ,
γ j , 3 u ^ j , 1 ( ω ) + ε j , 2 u ^ j , 2 γ j , 4 u ^ j + 1 , 1 ( ω ) = ε u ^ j , 2 ( ω ) ,
ε = ω 2 ,
ε j , 1 = ( K j 1 , 2 + K j , 1 ) / M j , 1 ,
ε j , 2 = ( K j , 1 + K j , 2 ) / M j , 2 ,
γ j , 1 = K j 1 , 2 / M j , 1 ,
γ j , 2 = K j , 1 / M j , 1 ,
γ j , 3 = K j , 1 / M j , 2 ,
γ j , 4 = K j , 2 / M j , 2 ,
where u ^ j , n ( ω ) is the Fourier transform of u j , n (t) formulated as
u ^ j , n ( ω ) = 1 2 π d t u j , n ( t ) exp ( i ω t ) ,
ε j , 1 and ε j , 2 represent squared on-site energy, and γ j , 1 , γ j , 2 , γ j , 3 and γ j , 4 represent squared hopping energy, respectively. By using Equations (A2a) and (A2b) with the left reservoir j = N L + 1 , N L + 2 , 1 , 0 , the H 2 O array j = 1 , 2 , N , and the right reservoir j = N + 1 , N + N R of whole structure shown in Figure 2, we can obtain phonon transport equations in matrix forms as follows:
H L τ L D 0 τ D L H D τ D R 0 τ R D H R u ^ L u ^ D u ^ R = ε u ^ L u ^ D u ^ R ,
H L = H N L + 1 V N L + 1 0 W N L + 2 H N L + 2 V 1 0 W 0 H 0 ,
H D = H 1 V 1 0 W 2 H 2 V N 1 0 W N H N ,
H R = H N + 1 V N + 1 0 W N + N R 1 H N + N R 1 V N + N R 1 0 W N + N R H N + N R ,
τ L D = 0 V 0 ,
τ D R = 0 V N ,
τ D L = W 1 0 ,
τ R D = W N + 1 0 .
Atomic displacement vectors u L , u D , and u R consist of u j of the j-th unit cell of each water molecule in the left reservoir, the H 2 O array, and the right reservoir, respectively.
u ^ L = u ^ N L + 1 u ^ 1 u ^ 0 ,
u ^ D = u ^ 1 u ^ 2 u ^ N ,
u ^ R = u ^ N + 1 u ^ N + 2 u ^ N + N R ,
u ^ j = u ^ j , 1 u ^ j , 2 .
The 2 × 2 matrices above consist of on-site energy and hopping energy of atoms in each water molecule of the j-th unit cell which are defined by following equations:
H j = ε j , 1 γ j , 2 γ j , 3 ε j , 2 ,
V j = 0 0 γ j , 4 0 ,
W j = 0 γ j , 1 0 0 .

Appendix A.2. Derivation of Eigenvalue Equations in Left and Right Reservoirs

We define κ and λ as Fourier transforms of atomic displacements of ± x propagation waves in the left reservoir, and μ and ν as those in the right reservoir, respectively. We assume that the Fourier transforms of atomic displacements u ^ j , n are linear combinations of ± x propagating waves in the left and right reservoirs as follows:
u ^ j , n = T j κ n + R j λ n , N L + 1 j 0 , n = 1 , 2 T j μ n + R j ν n , N + 1 j N + N R , n = 1 , 2
where T j and R j ( N L + 1 j 0 and N + 1 j N + N R ) are amplitudes of + and -x propagating waves, respectively.
Assuming that the phonon wavenumbers in the left and right reservoirs are k L ± and k R ± , respectively, we derive the following equations from the periodic boundary conditions of amplitude based on Bloch’s theorem for the left and right reservoirs:
T 1 = exp ( i k L + l 0 ) T 0 ,
T 1 = exp ( i k L + l 0 ) T 0 ,
R 1 = exp ( i k L l 0 ) R 0 ,
R 1 = exp ( i k L l 0 ) R 0 ,
T N + 2 = exp ( i k R + l N + 1 ) T N + 1 ,
T N = exp ( i k R + l N + 1 ) T N + 1 ,
R N + 2 = exp ( + i k R l N + 1 ) R N + 1 ,
R N = exp ( i k R l N + 1 ) R N + 1 .
In the left reservoir, applying R j = 0 and T j = 0 as the condition of no phonon propagating to ∓x directions, we obtain eigenvalue equations of plane waves of phonon propagating to ±x directions.
ε 0 , 1 ε γ 0 , 2 γ 0 , 1 exp ( i k L + l 0 ) γ 0 , 3 γ 0 , 4 exp ( i k L + l 0 ) ε 0 , 2 ε κ 1 κ 2 = 0 0 ,
ε 0 , 1 ε γ 0 , 2 γ 0 , 1 exp ( i k L l 0 ) γ 0 , 3 γ 0 , 4 exp ( i k L l 0 ) ε 0 , 2 ε λ 1 λ 2 = 0 0 .
Similarly, in the right reservoir, applying R j = 0 and T j = 0 as the condition of no phonon propagating to ∓x direction, we obtain eigenvalue equations of plane waves of phonon propagating to ±x directions.
ε N + 1 , 1 ε γ R 1 + γ R 2 + ε N , 2 ε μ 1 μ 2 = 0 0 ,
ε N + 1 , 1 ε γ R 1 γ R 2 ε N , 2 ε ν 1 ν 2 = 0 0 ,
γ R 1 ± γ N + 1 , 2 + γ N + 1 , 1 exp ( i k R ± l N + 1 ) ,
γ R 2 ± γ N + 1 , 3 + γ N + 1 , 4 exp ( i k R ± l N + 1 ) .
We can calculate dispersion curves as solutions to the eigenvalue Equations (A9a)–(A10b) in the left and right reservoirs, respectively, and these dispersion curves are the same as the solutions of the dispersion relation (1).

Appendix A.3. Non-Equilibrium Green’s Functions

The purpose of this section is to formulate phonon waves propagating in the water excited from left and right reservoirs, where we use the non-equilibrium Green’s functions defined by arranging Equation (A4a) [23]:
G L G L D 0 G D L G D G D R 0 G R D G R = ( ε + i η ) I H L τ L D 0 τ D L ( ε + i η ) I H D τ D R 0 τ R D ( ε + i η ) I H R 1 ,
where G L , G D , and G R are retarded Green’s functions in the left reservoir, water molecules, and right reservoir, respectively; G L D , G D L , G R D , and G D R are retarded Green’s functions at the interfaces between each reservoir and water. η is the imaginary part of the energy, which is a positive infinitesimal. Eliminating G L , G R , G L D , G D L , G R D , and G D R from Equation (A11), we can get an expression of G D as follows:
G D = ( ε + i η ) I H D Σ L Σ R 1 ,
where Σ L and Σ R are self-energies defined as
Σ L = τ D L g L τ L D ,
Σ R = τ D R g R τ R D .
g L and g R are equilibrium Green’s functions defined as follows:
g L = ( ε + i η ) I H L 1 ,
g R = ( ε + i η ) I H R 1 ,
g L = g N L + 1 , N L + 1 g N L + 1 , N L + 2 g N L + 1 , 0 g N L + 2 , N L + 1 g N L + 2 , N L + 2 g N L + 2 , 0 g 0 , N L + 1 g 0 , N L + 2 g 0 , 0 ,
g R = g N + 1 , N + 1 g N + 1 , N + 2 g N + 1 , N + N R g N + 2 , N + 1 g N + 2 , N + 2 g N + 2 , N + N R g N + N R , N + 1 g N + N R , N + 2 g N + N R , N + N R ,
g j , j = g j , 1 ; j , 1 g j , 1 ; j , 2 g j , 2 ; j , 1 g j , 2 ; j , 2 , j , j = N L + 1 , N L + 2 , , 0 , N + 1 , N + 2 , , N + N R .
Substituting Equations (A4e)–(A4h), (A14c) and (A14d) into (A13a) and (A13b), we can obtain the following equations for the self-energies:
Σ L = W 0 g 0 , 0 V 0 0 0 0 ,
Σ R = 0 0 0 V N g N + 1 , N + 1 W N + 1 .
We can define the equilibrium Green’s function equivalent for each unit cell at each atomic position (n, n ) in the left and right heat reservoirs as follows:
g n , n L = g j , n ; j , n , j , j = N L + 1 , N L + 2 , , 0 ,
g n , n R = g j , n ; j , n , j , j = N + 1 , N + 2 , , N + N R .
By executing the multiplications (A6c) j = 1 g 0 , 0 (A6b) j = 0 and (A6b) j = N g N + 1 , N + 1 (A6c) j = N + 1 , non-zero components in the right hand side of (A15a) and (A15b) become
W 1 g 0 , 0 V 0 = γ 1 , 1 g 2 , 2 L γ 0 , 4 0 0 0 ,
V N g N + 1 , N + 1 W N + 1 = 0 0 0 γ N , 4 g 1 , 1 R γ N + 1 , 1 .

Appendix A.4. Quantum Thermal Conductance in Water

Let us formulate equilibrium Green’s functions in reservoir, g L , g R , by the the following procedure consisting of basis transformations and so on:
  • Transforming wave functions consisting of the tight binding (TB) basis into those consisting of the plane wave (PW) basis.
  • Finding the PW basis solutions in semi-infinite ( N L , N R ) heat reservoir at both ends of water with using the periodic characteristics of the heat reservoir in thermal equilibrium.
  • Inverse transformation from PW basis wave functions to TB basis ones.
Based on the components of the j=1st and Nth unit cells of Equation (A4a), we can obtain the following equations governing the plane waves propagating to x direction in the left and right heat reservoirs, respectively:
ε u 1 = W 1 U L P L U L 1 u 1 + H 1 u 1 + V 1 u 0 ,
ε u N = W N u N 1 + H N u N + V N U R P R U R 1 u N .
We define P L and P R in Equations (A18a) and (A18b) as matrices in which (1,1) and (2,2) components of the diagonal terms are lined with propagation constants along x direction in the left and right heat reservoirs, respectively:
P L = exp ( i k L + l 0 ) 0 0 0 ,
P R = 0 0 0 exp ( i k R l N + 1 ) .
Next, we define the matrices U L and U R in Equations (A18a) and (A18b) as unitary matrices consisting of the eigenvectors of the eigenvalue Equations (A9a)– (A10b) as follows:
U L = κ 1 λ 1 κ 2 λ 2 ,
U R = μ 1 ν 1 μ 2 ν 2 .
We can simplify g L and g R as matrices with the same number of components as the number of the lattice points connected to the left and right heat reservoirs by using Equations (A18a) and (A18b), respectively:
g L = g 0 , 0 = ( ε + i η ) I H 1 W 1 U L P L U L 1 1 ,
g R = g N + 1 , N + 1 = ( ε + i η ) I H N V N U R P R U R 1 1 .
We can obtain spectral functions A L and A R proportional to local state densities A L / 2 π , and A R / 2 π in left and right reservoirs, respectively:
A L = G D Γ L G D ,
A R = G D Γ R G D .
We can obtain functions Γ L and Γ R by the following equations, which are matrices representing the phonon lifetime in the left and right reservoirs, respectively:
Γ L = i Σ L Σ L ,
Γ R = i Σ R Σ R .
We can then formulate the transmission matrix T of phonon waves traveling from the left to the right heat reservoirs via the interfaces of the water as follows [1,23,31]:
T = Γ L G D Γ R G D .
We can thus express the heat flow in water J by the following equation by integrating the product of the transmission function and the Bose-Einstein distribution [32]:
J = 0 d ω ħ ω 2 π n = 1 2 j = 1 N T j , n ; j , n N 0 ( ω ) N N + 1 ( ω ) ,
where n = 1 2 j = 1 N T j , n ; j , n is the transmission function expressed by the trace of T; N 0 ( ω ) and N N + 1 ( ω ) are Bose-Einstein distribution functions at j = 0 and N + 1 expressed as follows [37]:
N j ( ω ) = exp ( ħ ω / k B θ j ) 1 1 .
Thermal conductance in water σ can then be expressed as
σ = J θ 0 θ N + 1 .

References

  1. Dubi, Y.; Ventra, M.D. Colloquium: Heat flow and thermoelectricity in atomic and molecular junctions. Rev. Mod. Phys. 2011, 83, 131–155. [Google Scholar] [CrossRef]
  2. Yu, C.; Shi, L.; Yao, Z.; Li, D.; Majumdar, A. Thermal conductivity of individual silicon nanowires. Nano Lett. 2005, 5, 1842–1846. [Google Scholar] [CrossRef] [PubMed]
  3. Chiu, H.-Y.; Deshp, V.V.; Postma, H.W.C.; Lau, C.N.; Mikò, C.; Forrò, L.; Bockrath, M. Ballistic Phonon Thermal Transport in Multiwalled Carbon Nanotubes. Phys. Rev. Lett. 2005, 95, 226101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fujii, M.; Zhang, X.; Xie, H.; Ago, H.; Takahashi, K.; Ikuta, T.; Abe, H.; Shimizu, T. Measuring the Thermal Conductivity of a Single Carbon Nanotube. Phys. Rev. Lett. 2005, 95, 065502. [Google Scholar] [CrossRef] [PubMed]
  5. Yang, D.J.; Zhang, Q.; Chen, G.; Yoon, S.F.; Ahn, J.; Wang, S.G.; Zhou, Q.; Wang, S.G.; Li, J.Q. Thermal conductivity of multiwalled carbon nanotubes. Phys. Rev. B 2002, 66, 165440. [Google Scholar] [CrossRef] [Green Version]
  6. Gonnet, P.; Liang, Z.; Choi, E.S.; Kadambala, R.S.; Zhang, C.; Brooks, J.S.; Ben, W.; Leslie, K. Thermal conductivity of magnetically aligned carbon nanotube buckypapers and nanocomposites. Curr. Appl. Phys. 2006, 6, 119–122. [Google Scholar] [CrossRef]
  7. Berber, S.; Kwon, Y.; Tomànek, D. Unusually High Thermal Conductivity of Carbon Nanotubes. Phys. Rev. Lett. 2000, 84, 4613–4616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mingo, N.; Broido, D. Length Dependence of Carbon Nanotube Thermal Conductivity and the “Problem of Long Waves”. Nano Lett. 2005, 5, 1221–1225. [Google Scholar] [CrossRef]
  9. Lervik, A.; Bresme, F.; Kjelstrup, S. Heat transfer in soft nanoscale interfaces: The influence of interface curvature. Soft Matter 2009, 5, 2407–2414. [Google Scholar] [CrossRef] [Green Version]
  10. Tanaka, M.; Sato, M. Microwave heating of water, ice, and saline solution: Molecular dynamics study. J. Chem. Phys. 2007, 126, 034509. [Google Scholar] [CrossRef] [Green Version]
  11. Weissenborn, R.; Reinhardt, T.; Hansen, V.; Maret, G.; Gisler, T. A broadband waveguide for protein crystallography under intense microwave fields. Rev. Sci. Inst. 2004, 75, 5253–5256. [Google Scholar] [CrossRef] [Green Version]
  12. Fisher, T. Thermal Energy at the Nanoscale (Lessons from Nanoscience: A Lecture Notes Series); World Scientific: Singapore, 2014; Volume 3, Chapter 4. [Google Scholar]
  13. Rego, L.; Kiczenow, G. Quantized Thermal Conductance of Dielectric Quantum Wires. Phys. Rev. Lett. 1998, 81, 232–235. [Google Scholar] [CrossRef] [Green Version]
  14. Baffou, G.; Rigneault, H.; Marguet, D.; Jullien, L. Reply to “Validating subcellular thermal changes revealed by fluorescent thermosensors” and “The 105 gap issue between calculation and measurement in single-cell thermometry”. Nat. Methods 2015, 12, 803. [Google Scholar] [CrossRef] [PubMed]
  15. Kiyonaka, S.; Sakaguchi, R.; Hamachi, I.; Morii, T.; Yoshizaki, T.; Mori, Y. Validating subcellular thermal changes revealed by fluorescent thermosensors. Nat. Methods 2015, 12, 801–802. [Google Scholar] [CrossRef] [PubMed]
  16. Okabe, K.; Inada, N.; Gota, C.; Harada, Y.; Funatsu, T.; Uchiyama, S. Intracellular temperature mapping with a fluorescent polymeric thermometer and fluorescence lifetime imaging microscopy. Nat. Commun. 2012, 3, 705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kiyonaka, S.; Kajimoto, T.; Sakaguchi, R.; Shinmi, D.; Omatsu-Kanbe, M.; Matsuura, H.; Imamura, H.; Yoshizaki, T.; Hamachi, I.; Morii, T.; et al. Genetically encoded fluorescent thermosensors visualize subcellular thermoregulation in living cells. Nat. Methods 2013, 10, 1232–1238. [Google Scholar] [CrossRef]
  18. Kondepudi, D.; Prigogine, I. Modern Thermodynamics: From Heat Engines to Dissipative Structures; Wiley: Hoboken, NJ, USA, 1998; pp. 333–350. [Google Scholar]
  19. Baffou, G.; Rigneault, H.; Marguet, D.; Jullien, L. A critique of methods for temperature imaging in single cells. Nat. Methods 2014, 11, 899–901. [Google Scholar] [CrossRef]
  20. Lervik, A.; Bresme, F.; Kjelstrup, S.; Bedeaux, D.; Rubi, J.M. Heat transfer in protein–water interfaces. Phys. Chem. Chem. Phys. 2010, 12, 1610–1617. [Google Scholar] [CrossRef] [Green Version]
  21. Pandey, H.D.; Leitner, D.M. Small Saccharides as a Blanket around Proteins: A Computational Study. J. Phys. Chem. B 2018, 122, 7277–7285. [Google Scholar] [CrossRef]
  22. Leitner, D.M.; Pandey, H.D.; Reid, K.M.J. Energy Transport across Interfaces in Biomolecular Systems. J. Phys. Chem. B 2019, 123, 9507–9524. [Google Scholar] [CrossRef]
  23. Umegaki, T.; Ogawa, M.; Miyoshi, T. Investigation of electronic transport in carbon nanotubes using Green’s-function method. J. Appl. Phys. 2006, 99, 034307. [Google Scholar] [CrossRef] [Green Version]
  24. Elsaesser, T.; Bakker, H.J. Ultrafast Hydrogen Bonding Dynamics and Proton Transfer Processes in the Condensed Phase; Kluwer: Dordrecht, The Netherlands, 2002. [Google Scholar]
  25. Ashihara, S. Time-Resolved Observation of Molecular Dynamics in Liquid Water by Mid-Infrared Nonlinear Spectroscopy. Jpn. J. Opt. 2011, 40, 409–414. (In Japanese) [Google Scholar]
  26. Walrafen, G. Raman Spectrum of Water: Transverse and Longitudinal Acoustic Modes below 300 cm−1 and Optic Modes above 300 cm−1. J. Phys. Chem. 1990, 94, 2237–2239. [Google Scholar] [CrossRef]
  27. Nishimoto, Y.; Imamura, Z. Quantum Chemistry for Molecular Design; Kohdansya: Tokyo, Japan, 1991; pp. 105–108. [Google Scholar]
  28. Marechal, Y. The molecular structure of liquid water delivered by absorption spectroscopy in the whole IR region completed with thermodynamics data. J. Mol. Struct. 2011, 1004, 146–155. [Google Scholar] [CrossRef]
  29. Suzuki, K. Mizunohanashi; Kagaku Dohjin: Tokyo, Japan, 1997. (In Japanese) [Google Scholar]
  30. Paulsson, M.; Brandbyge, M. Transmission eigenchannels from nonequilibrium Green’s functions. Phys. Rev. B 2007, 76, 115117. [Google Scholar] [CrossRef] [Green Version]
  31. Bell, R. Large Scale Simulations of Conduction in Carbon Nanotube Networks. Ph.D. Thesis, The University of Cambridge, Cambridge, UK, 2014. Chapter 3. [Google Scholar]
  32. Zhang, W.; Fisher, T.S.; Mingo, N. The Atomistic Green’s Function Method: An Efficient Simulation Approach for Nanoscale Phonon Transport. Numer. Heat Transf. Part B 2007, 51, 333–349. [Google Scholar] [CrossRef]
  33. Suzuki, M.; Zeeb, V.; Arai, S.; Oyama, K.; Ishiwata, S. The 105 gap issue between calculation and measurement in single-cell thermometry. Nat. Methods 2015, 12, 802–803. [Google Scholar] [CrossRef]
  34. Ramakrishnan, G.; Gonzàlez-Jimènez, M.; Lapthorn, A.J.; Wynne, K. Spectrum of Slow and Super-Slow (Picosecond to Nanosecond) Water Dynamics around Organic and Biological Solutes. J. Phys. Chem. Lett. 2017, 8, 2964–2970. [Google Scholar] [CrossRef] [Green Version]
  35. Chieffo, L.; Shattuck, J.; Amsden, J.J.; Erramilli, S.; Ziegler, L.D. Ultrafast vibrational relaxation of liquid H2O following librational combination band excitation. Chem. Phys. 2007, 341, 71–80. [Google Scholar] [CrossRef]
  36. Balaji, P.V.; Qasba, P.K.; Rao, V.S.R. Flexibility of Sugar Chains and Molecular Dynamics Simulations. Trends Glycosci. Glycotechnol. 1994, 6, 93–94. [Google Scholar]
  37. Datta, S. Quantum Transport: Atom to Transistor; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
Sample Availability: Not available.
Figure 1. Definitions of displacements, masses, and spring constants of the elements H 2 (blue) and O (yellow) in the one-dimensional phonon model.
Figure 1. Definitions of displacements, masses, and spring constants of the elements H 2 (blue) and O (yellow) in the one-dimensional phonon model.
Molecules 25 01185 g001
Figure 2. One-dimensional phonon model of reservoir-water-reservoir structure. Each unit cell j has length l j , ( j = N L + 1 , N L + 2 , N + N R ) along the x-direction and temperature θ j , ( j = 0 , 1 , 2 , N + 1 ), where θ N L + 1 = θ N L + 2 = = θ 0 and θ N + 1 = θ N + 2 = = θ N + N R [23].
Figure 2. One-dimensional phonon model of reservoir-water-reservoir structure. Each unit cell j has length l j , ( j = N L + 1 , N L + 2 , N + N R ) along the x-direction and temperature θ j , ( j = 0 , 1 , 2 , N + 1 ), where θ N L + 1 = θ N L + 2 = = θ 0 and θ N + 1 = θ N + 2 = = θ N + N R [23].
Molecules 25 01185 g002
Figure 3. Absorption spectrum of water in the liquid phase [24]. Yellow and red broken lines represent the vibration modes (b) and (e), respectively.
Figure 3. Absorption spectrum of water in the liquid phase [24]. Yellow and red broken lines represent the vibration modes (b) and (e), respectively.
Molecules 25 01185 g003
Figure 4. Complex dispersion curves (upper) and transmission probability spectrum of the phonon (lower). The upper right inset plots enlarged dispersion curves in the frequency range of 642–645 Trad/s.
Figure 4. Complex dispersion curves (upper) and transmission probability spectrum of the phonon (lower). The upper right inset plots enlarged dispersion curves in the frequency range of 642–645 Trad/s.
Molecules 25 01185 g004
Figure 5. Modes of H 2 O molecular motions of (a) anti-symmetric stretching; (b) symmetric stretching; (c) bending; (d) libration; and (e) intermolecular vibration.
Figure 5. Modes of H 2 O molecular motions of (a) anti-symmetric stretching; (b) symmetric stretching; (c) bending; (d) libration; and (e) intermolecular vibration.
Molecules 25 01185 g005
Figure 6. A model of phonon transport through a water molecule between two reservoirs.
Figure 6. A model of phonon transport through a water molecule between two reservoirs.
Molecules 25 01185 g006
Figure 7. Temperature dependences of thermal conductance calculated by Equation (7). The parameters 1, 2, ⋯ 10 with unit of degrees Kelvin in the graph mean temperature differences θ 0 θ N + 1 between the left and right reservoirs.
Figure 7. Temperature dependences of thermal conductance calculated by Equation (7). The parameters 1, 2, ⋯ 10 with unit of degrees Kelvin in the graph mean temperature differences θ 0 θ N + 1 between the left and right reservoirs.
Molecules 25 01185 g007
Table 1. Wavenumbers and angular frequencies for various molecular vibration modes of water (c: light velocity).
Table 1. Wavenumbers and angular frequencies for various molecular vibration modes of water (c: light velocity).
Wave NumberAngular Frequency
Mode α k α (cm 1 )2 π k α c (Trad/s)
OH anti-symmetrica3509661.4
stretching
OH symmetricb3410642.8
stretching
OH bendingc1660312.9
Binding rotationd700131.9
(libration)
Intermoleculare19035.8
vibration
Table 2. Each vibrational mode (left), angular frequency range (center) and transmission probability (right) .
Table 2. Each vibrational mode (left), angular frequency range (center) and transmission probability (right) .
ModeAngular FrequencyTransmission
RangeProbability
( I ) acoustic wave 0 < ω < ω r 1
( I I ) attenuation ω r < ω < ω q 0
( I I I ) optical wave ω q < ω < ω p 1
†: (transmission probability) = tr ( T )/(its maximum value over ω ) [30].

Share and Cite

MDPI and ACS Style

Umegaki, T.; Tanaka, S. Nanoscale Quantum Thermal Conductance at Water Interface: Green’s Function Approach Based on One-Dimensional Phonon Model. Molecules 2020, 25, 1185. https://doi.org/10.3390/molecules25051185

AMA Style

Umegaki T, Tanaka S. Nanoscale Quantum Thermal Conductance at Water Interface: Green’s Function Approach Based on One-Dimensional Phonon Model. Molecules. 2020; 25(5):1185. https://doi.org/10.3390/molecules25051185

Chicago/Turabian Style

Umegaki, Toshihito, and Shigenori Tanaka. 2020. "Nanoscale Quantum Thermal Conductance at Water Interface: Green’s Function Approach Based on One-Dimensional Phonon Model" Molecules 25, no. 5: 1185. https://doi.org/10.3390/molecules25051185

APA Style

Umegaki, T., & Tanaka, S. (2020). Nanoscale Quantum Thermal Conductance at Water Interface: Green’s Function Approach Based on One-Dimensional Phonon Model. Molecules, 25(5), 1185. https://doi.org/10.3390/molecules25051185

Article Metrics

Back to TopTop