## Thermal conductivity of porous jennite by molecular dynamics method

Song-Nam Hong a , Chol-Jun Yu a, GLYPH&lt;3&gt; , Un-Song Hwang a , Chung-Hyok Kim a , Byong-Hyok Ri b a Chair of Computational Materials Design (CMD), Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People's Republic of Korea b Changhun-Kumsu Agency for New Technology Exchange, Mangyongdae District, Pyongyang, Democratic People's Republic of Korea

## Abstract

The thermal conductivity of porous jennite, as the major component of cement paste, and its porosity and temperature dependences are simulated by molecular dynamics methods using ClayFF force field. The porous jennite models with di GLYPH&lt;11&gt; erent porosities are created by removing atoms within the sphere from bulk jennite model. The thermal conductivity elements of bulk jennite calculated by equilibrium Green-Kubo method are almost identical to those by non-equilibrium M¨ uller-Plathe (MP) method at 300 K. The volumetric thermal conductivity of porous jennite is found to decrease from 1.141 W m K to 0.144 W m K as increasing the porosity from 0% to 72.22% at 300 K, / GLYPH&lt;1&gt; / GLYPH&lt;1&gt; following the empirical coherent potential model when the pore is assumed to be filled with air. As increasing temperature, the thermal conductivity is observed to increase from 240 K to 560 K and gradually decrease until 1100 K for porous jennites with the porosities of 0%, 15.99% and 32.79%.

Keywords: Jennite, Calcium Silicate Hydrate (C-S-H), Porous material, Thermal conductivity, Molecular dynamics

GLYPH&lt;3&gt; Corresponding author

Email address: cj.yu@ryongnamsan.edu.kp (Chol-Jun Yu)

## 1. Introduction

Building materials have always been indispensable for survival of human beings, and among them, concrete has been the most widely used since the invention of modern cement. While concrete consists of cement matrix and aggregates such as sand and gravel or industrial slag, its main binding phase is known to be calcium silicate hydrates (CaO GLYPH&lt;0&gt; SiO 2 GLYPH&lt;0&gt; H2 O; C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H) gel produced by cement hydration [1, 2]. The C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; Hgel is responsible for various important properties of concrete, including strength, toughness, and durability. Therefore, considerably mass of experimental and theoretical research has been carried out in revealing and improving mechanical properties of concrete and C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H gels [3, 4, 5, 6, 7, 8].

Although atomic modeling is widely accepted as a powerful tool to provide invaluable insights into the fundamental behavior of materials, the modeling of real C S GLYPH&lt;0&gt; GLYPH&lt;0&gt; H gel is known to be exceedingly challenging. This is partly associated with the di GLYPH&lt;14&gt; culties in identifying the atomistic structures of C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H gel by material characterization experiments, such as X-ray di GLYPH&lt;11&gt; raction (XRD), transition or scanning electron microscopy (TEM or SEM) and Raman spectroscopy [9]. In fact, the C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H gel in concrete is imperfect crystalline, in which the calcium to silicon (C S) ratio / is changeable from 1.2 to 2.1 according to the hydration process condition [10]. In such circumstances, they are often plausibly modeled by tobermorite and jennite minerals based on their structural similarities observed by XRD [11], although more realistic model with a C S ratio of 1.65 for C / GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H has been suggested by randomly removing charge-neutralized SiO 2 units from tobermorite 11 Å [12, 13, 14, 15, 16]. The mineral tobermorite has much smaller C S ratio of / GLYPH&lt;24&gt; 0.83 compared with experiment, whereas for jennite it is reasonable to be about 1.5 due to the corrugation of Ca GLYPH&lt;0&gt; O sheets [17]. So far, jennite and tobermorite family have been intensively studied as substitutes for the C-S-H gels. Bonaccorsi et al . have done pioneering work on characterizing the crystalline structures of jennite [18], tobermorite 11 Å [19] and tobermorite 14 Å [20] through XRD experiments. The structural and mechanical properties of jennite and tobermorite have been investigated using firstprinciples calculations, [17, 21, 22, 23, 24, 25, 26] and molecular dynamics (MD) sim-

ulations [27, 28, 29, 30]. However, there is only a limited number of studies on thermal properties of jennite and tobermorite as well as C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H gel [31, 32, 33, 34, 35], in spite of its importance in design and construction of energy e GLYPH&lt;14&gt; cient buildings with concrete [36]. In particular, studies of revealing the thermal conductivity of jennite at molecular scale are truly scarce.

Recently, there is a growing interest in high thermal insulation materials in construction to save energy use in buildings [37, 38]. In general, thermal insulation materials are porous [39, 40]. The thermal conductivity of such porous materials decreases as their porosity increases due to an increasing amount of inside small-sized air pores. Porous concrete is in particular advantageous due to its high thermal insulation, low density, fire-proofing and acoustic insulation [41]. Qomi et al . [34] reported the thermal properties of C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H, such as specific-heat capacity, thermal expansion coe GLYPH&lt;14&gt; cient and thermal conductivity, by means of MD method using core-only CSH force field potential [42]. Jin et al . [43] proposed the fractal model for predicting the thermal conductivity of porous concrete based on their measurements of autoclaved aerated concrete with various porosities.

In this work, we aim to determine the thermal conductivity of jennite and its dependence on the porosity and temperature by using MD simulations. We organized the paper as follows. In Section 2, we described computational methods, including the employed force field potential, jennite models and MD methods for predicting the thermal conductivity. Section 3 shows the results and discussion, where comparisons with the previous MD works, experimental results, and predictions based on the continuum models for heat transfer in porous materials, were provided. Temperature dependence of thermal conductivity was also presented together with a discussion in comparison with the available literature data. Finally, we provided the conclusions of this work in Section 4.

## 2. Methods

## 2.1. Force field

The accuracy of classical MD simulations strongly depends on the selected interatomic force field potential. To date, many force fields have been developed for

cementitious materials [44, 45]. They are for example ClayFF [46], IFF (interface force field) [47], CementFF [48], ReaxFF [49], and CSH-FF [42]. Among them, the ClayFF potential was selected in this work, since it has been completely tested through lots of studies on various properties of clays and cement materials [30, 50, 51, 52]. As Mutisya et al . [29] pointed out, moreover, the ClayFF is able to capture thermodynamic properties of both hydrous and anhydrous cement phases.

In the ClayFF formalism, the total energy of a molecular system is represented as follows [44, 46],

$$E _ { \text{total} } = E _ { \text{Coul} } + E _ { \text{vdW} } + E _ { \text{bond} } + E _ { \text{angle} }$$

where E Coul is the Coulomb interaction energy, E vdW is the short-range van der Waals (vdW) energy, E bond is the bond stretching term, and E angle is the angle bending term, respectively. The electrostatic interaction energy between the charged particles is described by the Coulomb's law as,

$$E _ { \text{Coul} } = \frac { e ^ { 2 } } { 4 \pi \varepsilon _ { 0 } } \sum _ { i \neq j } \frac { q _ { i } q _ { j } } { r _ { i j } }$$

where qi is the atomic partial charge of i th particle, ri j is the distance between the i th and j th particles, e is the elementary charge, and " 0 is the dielectric permittivity of vacuum. The attractive vdW dispersion energy is written as [44, 46],

$$E _ { \text{vd} } = \sum _ { i \neq j } D _ { i j } ^ { 0 } \left [ \left ( \frac { R _ { i j } ^ { 0 } } { r _ { i j } } \right ) ^ { 1 2 } - 2 \left ( \frac { R _ { i j } ^ { 0 } } { r _ { i j } } \right ) ^ { 6 } \right ]$$

where D ij 0 = q D i 0 D j 0 and R i j 0 = ( R i 0 + R j 0 ) = 2. D i 0 , R i 0 are the empirical parameters derived from the experimental data for physical property. The bonded interactions, which are used to describe the O GLYPH&lt;0&gt; H bonding in water molecules and hydroxyl groups, and the covalent bonding in polyatomic species, are simplified as harmonic terms as follows [44, 46],

$$E _ { \text{bond} } = \sum _ { i \neq j } k _ { 1 } ( r _ { i j } - r _ { 0 } ) ^ { 2 }$$

$$E _ { \text{angle} } = \sum _ { i \neq j \neq k } k _ { 2 } ( \theta _ { i j k } - \theta _ { 0 } ) ^ { 2 }$$

where k 1 and k 2 are the harmonic force constant parameters, and r 0 and GLYPH&lt;18&gt; 0 are the equilibrium values of bond length and bond angle respectively. The parameters qi ,

D i 0 , R i 0 and k for Ca, Si, O and H atoms considered in this work were generated by using the msi 2 lmp tool of Large Scale Atomic Molecular Massively Parallel Simulator / (LAMMPS) [53], and given in Table S1 GLYPH&lt;0&gt; S3 (Supplementary data).

## 2.2. Models for porous jennite

Figure 1: Polyhedral view of atomistic structures of jennite models. (a) Unit cell of jennite crystal, b) supercell for bulk jennite model with 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 unit cells and c) its equilibrated phase, and porous jennite models created by removing atoms within sphere with radii of (d) 5 Å, (e) 7 Å and (f) 9 Å. Blue, yellow, red, white, and green balls indicate Ca, Si, O, H, and interlayer Ca atoms, and yellow polyhedra represent the SiO 4 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-c11f9b43-73ee-4823-a402-5f51fbf3340a.jpg)

The chemical formula of jennite is Ca Si O 9 6 18 (OH) 6 GLYPH&lt;1&gt; 8H 2 O, and its crystalline lattice is triclinic with P ¯ 1 space group [17, 18]. Fig. 1(a) shows the unit cell of jennite (see jennite-unitcell.cif in Supplementary material). For MD simulations, we constructed a supercell for bulk jennite model by using 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 unit cells, which contains 2484 atoms, and subsequently carried out constant pressure and temperature ( NPT ) simulation with a time step of 0.1 fs for 100 ps at 300 K and atmospheric pressure, then it's size was about 32 6 : GLYPH&lt;2&gt; 29 6 : GLYPH&lt;2&gt; 32 8 Å, as shown in Fig. 1(b). : At this NPT simulation the cuto GLYPH&lt;11&gt; radii for the LJ and Coulombic interactions were set to be 15 Å. The periodic boundary condition was applied in all three crystallographic directions, and the

pppm method was used for calculation of electrostatic term. Nose-Hoover thermostat and barostat were used to control temperature and pressure of the system, and the velocity Verlet time integration algorithm was used to update the position and velocity of atoms. Fig. 1(c) shows the equilibrated atomistic structure of the jennite model. Table 1 presents the determined lattice parameters, which are in good agreement with experiment with a maximal relative error of 2.73% [18].

To verify the validity of jennite model and the force field, we also evaluated the radial distribution function (RDF) of Ca GLYPH&lt;0&gt; O in comparison with available experiment and previous MD simulations. According to the di GLYPH&lt;11&gt; erence of local structures, calcium atoms are classified into sheet (Cas) and interlayer (Caw) atoms, while oxygen atoms are divided into surrounding (Os) and water (Ow) oxygen atoms. Therefore, we considered four di GLYPH&lt;11&gt; erent RDFs for Ca GLYPH&lt;0&gt; O bonds (see Fig. S1). For the Cas GLYPH&lt;0&gt; O bonds, the first peaks were found at 2.41 and 2.59 Å, which are slightly larger than the previous MDresults of 2.37 and 2.57 Å [27] while in reasonable agreement with the experimental value of 2.45 GLYPH&lt;6&gt; 0.06 Å [54]. The distances between Caw and Ow or Oh (hydroxyl groups) atoms were estimated to be 2.35 or 2.50 Å, being in good agreement with the XRD experimental data of 2.36 or 2.50 Å [18]. These indicate that our jennite model can be used safely and the ClayFF employed in this work is reasonably accurate and transferable for jennite.

Then, the porous jennite models were created by removing all the atoms within sphere with a certain radius located at the center of perfect jennite model. Note that if one atom of a polyatomic molecule (e.g., CaO, SiO , H O) belonged to the spherical 2 2 region, the molecule was totally removed. We increased the radius gradually from 3 Å to 11 Å, and estimated the corresponding mass density and porosity of each porous model. The porosity can be estimated by using the mass densities of perfect model GLYPH&lt;26&gt; 0

Table 1: Lattice parameters of unit cell of jennite crystal, equilibrated through NPT run, in comparison with experiment.

|          |   a (Å) |   b (Å) |   c (Å) |   GLYPH<11> ( GLYPH<14> ) |   GLYPH<12> ( GLYPH<14> ) |   GLYPH<13> ( GLYPH<14> ) |
|----------|---------|---------|---------|---------------------------|---------------------------|---------------------------|
| ClayFF   |  10.865 |   7.402 |  10.946 |                    102.06 |                     94.57 |                    110.33 |
| Exp [18] |  10.576 |   7.265 |  10.931 |                    101.3  |                     96.98 |                    109.65 |

Table 2: Sphere radius ( r ), total number of atoms ( N tot ), number of deleted atoms ( N del ), mass density ( GLYPH&lt;26&gt; ), and porosities estimated using number of atoms ( PN ) and mass density ( P GLYPH&lt;26&gt; ) in perfect and porous jennite models.

|   model |   r (Å) |   N tot |   N del |   GLYPH<26> (g / cm 3 ) |   P N (%) |   P GLYPH<26> (%) |
|---------|---------|---------|---------|-------------------------|-----------|-------------------|
|       1 |     0   |    2484 |       0 |                   2.205 |      0    |              0    |
|       2 |     3   |    2346 |     138 |                   2.099 |      5.56 |              4.79 |
|       3 |     5   |    2208 |     276 |                   1.976 |     11.11 |             10.39 |
|       4 |     5.5 |    2070 |     414 |                   1.852 |     16.67 |             15.99 |
|       5 |     6   |    1863 |     621 |                   1.667 |     25    |             24.39 |
|       6 |     7   |    1656 |     828 |                   1.482 |     33.33 |             32.79 |
|       7 |     9   |    1173 |    1311 |                   1.05  |     52.78 |             52.39 |
|       8 |    11   |     690 |    1794 |                   0.612 |     72.22 |             72.22 |

and porous model GLYPH&lt;26&gt; p as follows,

$$P = \left ( 1 - \frac { \rho _ { p } } { \rho _ { 0 } } \right ) \times 1 0 0 \%$$

For reference, we also provided the porosity calculated using the number of atoms instead of density. Table 2 lists the characteristics of porous jennite models with various porosities from 0 to 72.22%. Fig. 1(d), (e) and (f) show the equilibrated atomistic structures of porous jennite models with the porosities of 10.39%, 32.79% and 52.39% respectively.

## 2.3. Molecular dynamics methods for predicting the thermal conductivity

The thermal conductivity K of insulating solid is defined as a second-order tensor coe GLYPH&lt;14&gt; cient relating the temperature gradient r T to the heat flux q through Fourier's law, q = K GLYPH&lt;1&gt; r T . Typically, two di GLYPH&lt;11&gt; erent methods have been developed to calculate K via MD simulation: (1) equilibrium Green-Kubo (GK) method [55, 56] and (2) nonequilibrium M¨ uller-Plathe (MP) method [57]. The former method does not impose temperature gradients and yields the full thermal conductivity tensor via one time simulation, while the latter approach imposes a heat flux on the simulation box and gives the heat conductivity separately in each direction for an anisotropic material. In this work, we applied these two approaches in order to guarantee an accuracy of calculation, relying mostly on the MP method.

In the GK approach, once the molecular system is equilibrated through NPT simulation, the net heat flux fluctuates around zero during NVE run. The dissipation rate of heat flux fluctuations is evaluated, and the thermal conductivity tensor for

an anisotropic materials is measured at equilibrium by using the following GK relation [34, 58],

$$K _ { \alpha \beta } = \frac { V } { k _ { B } T ^ { 2 } } \int _ { 0 } ^ { \infty } \left \langle q _ { \alpha } ( t ) \otimes q _ { \beta } ( 0 ) \right \rangle d t$$

where GLYPH&lt;10&gt; is the dyadic product in tensor notation, kB the Boltzmann constant, V the volume of system, and GLYPH&lt;11&gt; , GLYPH&lt;12&gt; the components of coordinate system. In this equation, GLYPH&lt;30&gt; GLYPH&lt;11&gt;GLYPH&lt;12&gt; ( ) t = h q GLYPH&lt;11&gt; ( ) t GLYPH&lt;10&gt; q GLYPH&lt;12&gt; (0) i = h q GLYPH&lt;11&gt; (0) GLYPH&lt;10&gt; q GLYPH&lt;12&gt; (0) i is the second-order tensor known as the heat flux autocorrelation function (HFACF), and thus, the thermal conductivity can be calculated from integration of HFACF. The heat flux in a polyatomic system with a number of atom N is written as follows [34],

$$q _ { \alpha } ( t ) = \frac { 1 } { V } \frac { d } { d t } \sum _ { i = 1 } ^ { N } r _ { i, \alpha } H _ { i }$$

where r i is the position vector of the i th atom, and Hi is the Hamiltonian, which is composed of kinetic and potential (Eq. (1)) energies. The GK approach has been applied to C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H [34] and amorphous silicon dioxide thin films [59]. In this work, the equilibrium GK approach was applied to the perfect jennite crystal with the 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 supercell using LAMMPS [53]. After equilibration of the system through NPT run, we subsequently performed NVE simulation with a time step of 0.1 fs for 20 ns (2 GLYPH&lt;2&gt; 10 5 correlation time steps).

In the non-equilibrium MP method, thermal conductivity can be calculated by imposing a heat flux on the simulation box and measuring the resultant temperature gradient [57, 60, 61, 62]. If the heat propagates along GLYPH&lt;11&gt; direction in the simulation box, the thermal conductivity in this direction is calculated by

$$K _ { \alpha } = - \frac { q _ { \alpha } } { \partial T / \partial r _ { \alpha } } = - \frac { q _ { \alpha } } { \Delta T / \Delta r _ { \alpha } }$$

GLYPH&lt;11&gt;

As implemented in the LAMMPS [53], the simulation box is divided into an even number of parallel slices, and the temperature of each slice is evaluated at every time step and then averaged over the whole time steps [61]. The heat flux is imposed through velocity rescaling, i.e., by exchanging velocities between atoms with the largest kinetic energy of the first slice and atoms with the lowest kinetic energy of the middle slice at a certain simulation time step. After reaching steady state by imposing heat flux

for su GLYPH&lt;14&gt; ciently long time, a temperature profile T r ( GLYPH&lt;11&gt; ) is evaluated and thus the thermal conductivity is ready for estimation by using Eq. (9). For an anisotropic material, this procedure has to be performed separately to obtain the thermal conductivity in each direction of coordinate system. The MP method is easy to understand and has moderate simulation time, so it has been used in many studies of porous materials [59, 61, 62].

In this work, we applied the MP method mostly to the porous jennite models with various porosities for estimating their thermal conductivities at di GLYPH&lt;11&gt; erent temperatures, using the LAMMPS program [53]. The porous jennite models, constructed by removing atoms within the spherical region as described above, were equilibrated by performing constant volume and constant temperature ( NVT ) simulations with a time step of 0.1 fs for 100 ps. Then, the simulation boxes were created by extending these equilibrated porous jennite models by from two to ten times in [100], [010] and [001] crystallographic directions, and also equilibrated by NVT run. To calibrate the temperature profile and the thermal conductivity, the extended simulations boxes are divided into 20 plane-parallel slices, and micro canonical ensemble ( NVE ) simulations were performed with a time step of 0.1 fs during 200 ps. Like the application of GK approach, the Nose-Hoover thermostat and barostat, and velocity Verlet time integration algorithm were adopted.

## 3. Results and discussion

## 3.1. Thermal conductivity of jennite by Green-Kubo method

We first presented the thermal conductivity of perfect jennite calculated by GK approach using the 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 supercell model. Figure 2 shows the HFACF along the [100], [010] and [001] directions, and thus the principal elements of thermal conductivity tensor at temperature of 300 K. After approximately 200 fs of NVE simulation at equilibrium state, the HFACF was found to fluctuate around zero while the thermal conductivity to converge into a certain value. The principal values were determined to be K 1 = 1 125, : K 2 = 1 393 and : K 3 = 0 921 W m K along the [100], [010] and [001] di-: / GLYPH&lt;1&gt; rections respectively, and the average volumetric value was estimated to be K v = 1 146 : WmK. As expected earlier, the thermal conductivity of jennite was confirmed to be / GLYPH&lt;1&gt; anisotropic like K 2 &gt; K 1 &gt; K 3 , which is di GLYPH&lt;11&gt; erent from K 1 &gt; K 2 &gt; K 3 for C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H

predicted by MD [34], possibly due to di GLYPH&lt;11&gt; erences in model and force field. The reason of such anisotropy is that fairly heat conductive Si GLYPH&lt;0&gt; O bonds attribute to higher thermal conduction along the [010] direction (the so-called Si GLYPH&lt;0&gt; O dreierketten chains run along the [010] direction [18]), whereas in the [001] direction water molecules and looser Ca GLYPH&lt;0&gt; O bonds scatter phonons, resulting in lower thermal conduction in this direction [27, 34].

## 3.2. Influence of system length on thermal conductivity

In the previous works, it has been found that in the MP method the accuracy of thermal conductivity calculation severely depends on the system length, whereas the e GLYPH&lt;11&gt; ect of its cross-section area is independent for the case of amorphous nanoporous silica. [61]. Therefore, we tested the influence of simulation box length on the thermal conductivity of bulk jennite to determine the proper length. In this test, the perfect

⋅Κ)

Figure 2: (a) The heat flux autocorrelation function (HFACF) and (b) the thermal conductivity along the crystallographic directions as functions of correlation time. Insets show the fluctuations of HFACF and thermal conductivities with dashed lines indicating average values estimated from 1000 step.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-eb1f61c4-47b7-44ba-a07f-25b20d7fe868.jpg)

jennite model was used. Due to the anisotropy of thermal conductivity, the jennite model was repeated from two to ten times in the [100], [010] and [001] directions, yielding the simulation boxes with a range of system length from GLYPH&lt;24&gt; 60 to GLYPH&lt;24&gt; 300 Å, cross-section area of about 30 GLYPH&lt;2&gt; 30 Å 2 , and number of atoms from 4968 to 24 840.

As depicted in Fig. 3(b), the simulation boxes were divided into 20 plane slices perpendicular to the repeating direction (heat flux propagation direction), of which the 11th slice plays the role of heat source and the first slice is for heat sink. It is worth noting that this configuration is due to the periodic boundary condition in all the three directions, while under the fixed boundary condition the first and last slices can be used for heat source and sink [62]. Temperature was fixed at 300 K through NVT simulation, and subsequently NVE runs were performed with a time step of 0.1 fs during a total time of 200 ps. Figures 3(c)-(e) show temperature profiles in the simulation boxes with increasing lengths along the heat flux directions of [100], [010] and [001]. It was observed that, when the system length was not enough long, the temperature reduction going from the heat source to the heat sink slices did not follow the straight line but rather meandering path, raising di GLYPH&lt;14&gt; culty and ambiguity in measuring the temperature gradient. For the system lengths over GLYPH&lt;24&gt; 150 Å (extending supercell model five times), the temperature profile became almost the straight line.

In Table 3, we summarize the principal and volumetric averaged thermal conductivities of jennite crystal as increasing the system length. As can be seen in Fig. 4, the principal thermal conductivities converged into the reliable values over the system length of GLYPH&lt;24&gt; 150 Å, which is comparable with the previous value of 100 Å for amorphous silica [61]. For the sake of high reliability, we selected the appropriate system

Table 3: Principal ( K 1, K 2, K 3) and volumetric averaged ( K v) thermal conductivities of jennite crystal as increasing the system length, i.e. number of repeating cells N Re.

|   N Re |   K 1 (W / mK) GLYPH<1> |   K 2 (W / mK) GLYPH<1> |   K 3 (W / mK) GLYPH<1> |   K v (W / mK) GLYPH<1> |
|--------|-------------------------|-------------------------|-------------------------|-------------------------|
|      2 |                   1.521 |                   1.357 |                   0.566 |                   1.148 |
|      3 |                   1.239 |                   1.506 |                   0.748 |                   1.164 |
|      4 |                   1.182 |                   1.339 |                   0.886 |                   1.136 |
|      5 |                   1.14  |                   1.378 |                   0.917 |                   1.145 |
|      8 |                   1.1   |                   1.354 |                   0.967 |                   1.141 |
|     10 |                   1.085 |                   1.368 |                   1.003 |                   1.152 |

length as about 240 Å (repeating jennite model 8 times, 19 872 atoms) by considering the computational cost and statistical error.

When compared with our previous result by GK approach, the obtained thermal conductivities, K 1 = 1 100 : GLYPH&lt;6&gt; 0 025 W m K, : / GLYPH&lt;1&gt; K 2 = 1 354 : GLYPH&lt;6&gt; 0 034 W m K, : / GLYPH&lt;1&gt; K 3 = 0 967 : GLYPH&lt;6&gt; 0 020 W m K and : / GLYPH&lt;1&gt; K v = 1 141 : GLYPH&lt;6&gt; 0 026 W m K, are in good agreement each : / GLYPH&lt;1&gt; other. Due to a lack of experimental data for jennite, we compared our results with the previous MD works for similar materials. For C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H, Qomi et al . [34] obtained 0 98 : GLYPH&lt;6&gt; 0 2 W m K, which is comparable with ours. : / GLYPH&lt;1&gt; Wu et al . [63] also reported the similar value of concrete hydration product as 1.015 W m K through the volumetric / GLYPH&lt;1&gt; weighted average of the alite and belite. These indicate that our models, adopted force field, and computational method are su GLYPH&lt;14&gt; ciently accurate for study of porous jennite in the following subsections.

## 3.3. Porosity dependence of thermal conductivity

We further investigated the porous jennite with porosities from 0% to 72.22% at fixed temperature of 300 K, comparing with the di GLYPH&lt;11&gt; erent models for describing the e GLYPH&lt;11&gt; ective thermal conductivity of porous material and with the experimental data. It should be noted that introducing pores into the perfect crystal or media reduces the thermal conductivity substantially, and thus, porous media acts as a thermal insulator or a barrier to heat removal. As mentioned above, the simulation boxes extended by 8 times in each direction were employed (see Fig. 3(a) for porous jennite with a porosity of 24.39%). However, for the case of 72.22% porosity, the number of atoms was too small (5 520 atoms) and thus the width and height of simulation box were again doubled, keeping the length unchanged, so that the number of atoms became enough large (22 080).

Before presenting the results of this work, let us consider the previously suggested models that relate the thermal conductivity of porous media to the porosity. The Garnett model was first derived to determine the electrical conductivity of a binary composite, in which spherical inclusions with random size were dispersed, as follows [64],

$$K ^ { G } = K _ { b } \frac { 2 K _ { b } + K _ { p } - 2 p ( K _ { b } - K _ { p } ) } { 2 K _ { b } + K _ { p } + p ( K _ { b } - K _ { p } ) }$$

where Kb and Kp are the thermal conductivities of bulk and pore, and p is the porosity. Landauer [65] proposed the coherent potential (CP) model to estimate the conductivity of a binary composite made of matrix and spherical inclusions as follows,

$$K ^ { \text{CP} } = \frac { 1 } { 4 } \left \{ ( 3 p - 1 ) K _ { p } + ( 2 - 3 p ) K _ { b } + \\ \left [ \left ( ( 3 p - 1 ) K _ { p } + ( 2 - 3 p ) K _ { b } \right ) ^ { 2 } + 8 K _ { b } K _ { p } \right ] ^ { \frac { 1 } { 2 } } \right \} \ \ ( 1 1 )$$

By using the CP model, Cahill and Allen [66] well predicted the thermal conductivity of porous material with a porosity of 30%. The Russel model was developed to derive the e GLYPH&lt;11&gt; ective thermal conductivity of a dry porous material, in which cubical inclusions with uniform pores were evenly distributed [67],

$$K ^ { R } = K _ { b } \frac { \ v p ^ { 2 / 3 } + 1 - p ^ { 2 / 3 } } { v ( p ^ { 2 / 3 } - p ) + 1 + p - p ^ { 2 / 3 } }$$

where v = Kp = Kb . The Russell model assumed a parallel heat flow, whereas the Frey model treated a series heat flow [68],

$$K ^ { F } = K _ { b } \frac { v ( 1 - p ^ { 1 / 3 } + p ) + p ^ { 1 / 3 } - p } { v ( 1 - p ^ { 1 / 3 } ) + p ^ { 1 / 3 } }$$

The Parallel and Series models were developed to give the maximum and minimum of thermal conductivity of a porous material because planes of constituent phases are in parallel and in series with the heat flow direction, as follows,

$$K ^ { P } & = K _ { b } ( 1 - p ) + K _ { p } p & ( 1 4 )$$

$$K ^ { S } = \frac { 1 } { ( 1 - p ) / K _ { b } + p / K _ { p } }$$

Meanwhile, Qomi et al . [34] derived the self-consistent (SC) model to estimate the homogenized thermal conductivity of matrix-inclusion system, as follows,

$$K _ { v } ^ { S C } & = \frac { \sum _ { s = 1 } ^ { n _ { p } } f _ { s } K _ { v } ^ { s } B _ { s } ^ { s p h } } { \sum _ { s = 1 } ^ { n _ { p } } f _ { s } B _ { s } ^ { s p h } }$$

where K s v is the volumetric thermal conductivity of spherical inclusion, f s is the volume fraction between the matrix and inclusion, and Bs sph = 3 Kv SC = (2 Kv SC + K s v ) is the spherical localization factor of the s th phase.

In Table 4, we present the thermal conductivities of the porous jennite with different porosities from 0 to 72.22% calculated by applying the MP method at 300 K. Figure 5 shows normalized thermal conductivity of the porous jennite as a function of porosity, comparing with the e GLYPH&lt;11&gt; ective thermal conductivities of porous media derived from the various models described above. As expected, the thermal conductivity of porous jennite was observed to decrease as increasing the porosity, since the thermal conductivity of air (0.026 W m K [69]) is much smaller than that of jennite. / GLYPH&lt;1&gt; As can be seen in Fig. 5, although all models exhibit decreasing tendency of e GLYPH&lt;11&gt; ective thermal conductivity as increasing the porosity, the CP and SC models (they are almost identical each other) were found to be in quantitative agreement with our MD simulation results. It should be noted that an iterative procedure is required in the calculation with the SC model, and thus, in order to calculate the thermal conductivity of porous jennite by using the SC model we made use of the result from the Parallel model as the initial value and iterated the calculation for 10 times during the self-consistent cycle, reaching convergence after 7 or 8 times. For the case of amorphous SiO , the CP 2 model has also been found to reasonably produce its thermal conductivity for a range of porosities such as 10 GLYPH&lt;24&gt; 35% [61] and 0.54 GLYPH&lt;24&gt; 24.15% [62]. Furthermore, Santos [70] revealed that the thermal conductivity of conventional aluminous refractory concrete exhibited the decreasing tendency as increasing the porosity from 0 to GLYPH&lt;24&gt; 36.37% at various temperatures, nearly according to that from the CP model. Meanwhile, Coquil et al . [61] found a slight overestimation of the CP model compared with the MD result, whereas Zhu et al . [62] found an underestimation when the porosity was 24.15%. But

Table 4: Principal ( K 1, K 2, K 3) and volumetric averaged ( K v) thermal conductivities of the porous jennite with various porosities ( P ) at 300 K.

GLYPH&lt;26&gt;

|   P GLYPH<26> (%) |   K 1 (W / mK) GLYPH<1> |   K 2 (W / mK) GLYPH<1> |   K 3 (W / mK) GLYPH<1> |   K v (W / mK) GLYPH<1> |
|-------------------|-------------------------|-------------------------|-------------------------|-------------------------|
|              0    |                   1.1   |                   1.354 |                   0.967 |                   1.141 |
|              4.79 |                   0.996 |                   1.253 |                   0.867 |                   1.039 |
|             10.39 |                   0.932 |                   1.159 |                   0.806 |                   0.966 |
|             15.99 |                   0.828 |                   0.986 |                   0.732 |                   0.849 |
|             24.39 |                   0.66  |                   0.77  |                   0.671 |                   0.701 |
|             32.79 |                   0.525 |                   0.684 |                   0.545 |                   0.585 |
|             52.39 |                   0.2   |                   0.384 |                   0.338 |                   0.307 |
|             72.22 |                   0.013 |                   0.273 |                   0.146 |                   0.144 |

in our study the MD results were in good agreement with those from the CP model over the porosity range from 0 to 52.39% and not ruled out for the porosity of 72.22%. Based on this analysis, it can be concluded that the CP model is able to well predict the thermal conductivity of not only amorphous porous silica but also porous cement based materials.

Next, we made a careful analysis of our MD result in comparison with those from the CP model, in which the pores were assumed to be filled with air or water, and from experiment. Whether the pores of porous cement hydration product would be filled with air or water (or even vacuum) depends on its hardening process condition. If the pores are assumed to be filled with water, the porous jennite has the maximal thermal conductivity, whereas filling with air will give the minimal values, as obtained by the CPmodel shown in Fig. 6. It is worth noting that for the case of water filling the thermal conductivity of porous jennite should decrease with increasing the porosity as well, due to smaller conductivity of water (0.6062 W m K [69]) than jennite. / GLYPH&lt;1&gt; According to the work of Constantinides and Ulm [74], the porosities of high and low density C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H were found to be 24 and 36% respectively, and therefore, if hydrated cement paste would be composed of only jennite, its thermal conductivity should be ranged from 0.553 W m K with air filling to 0.992 W m K with water filling on the CP model lines. / GLYPH&lt;1&gt; / GLYPH&lt;1&gt;

In our previous work [77], we demonstrated that when applying the multi-functional composite coating material [78] made from silica alumina-rich minerals to the concrete / wall, high dense ceramic layer with a thickness of GLYPH&lt;24&gt; 1 mm and a porosity of 10.05% can

Table 5: Thermal conductivity ( K ) of the cement hydration product with di GLYPH&lt;11&gt; erent water to cement ratios (w c) /

|                       | K (W / mK) GLYPH<1>        | w / c    |
|-----------------------|----------------------------|----------|
| This work             | 1.030                      | -        |
| Yoon et al . [31]     | 1.28                       | 0.4      |
| Xu et al . [32]       | 0.53                       | 0.35     |
| Bentz [33]            | 0.90 GLYPH<24> 1.05        | 0.3, 0.4 |
| Hansen et al . [75]   | 0.88 (early)               | 0.5      |
| Mounanga et al . [76] | 1.0 (fresh) 1.07 (28 days) | 0.348    |
| Wu et al . [63]       | 0 : 9568 GLYPH<6> 0 : 0257 | 0.45     |

In this work, laser flash method was used to measure thermal conductivity.

Data of Wu et al . [63] was from calculation, and others were from experiment.

be formed on the wall surface by two-step hydration of cement and mineral powders. By using the laser flash method, we measured the thermal conductivity of this ceramic layer to be 1.030 W m K at 223 K (see Table S2), which is placed between the red and / GLYPH&lt;1&gt; blue lines as shown in Fig. 6. In Table 5, we summarize the thermal conductivity of cement hydration products with di GLYPH&lt;11&gt; erent water to cement ratios (w c), measured from / experiment or calculated from theory. Wu et al . [63] calculated the thermal conductivity of hardened cement paste to be 0 96 : GLYPH&lt;6&gt; 0 03 W m K when its hydration degree was : / GLYPH&lt;1&gt; 0.95 and w c / = 0.45. When compared with the available experimental data of cement paste with di GLYPH&lt;11&gt; erent w c ratios such as 0 90 / : GLYPH&lt;24&gt; 1 05 W m K with the hydration degree : / GLYPH&lt;1&gt; of about 0 3 : GLYPH&lt;24&gt; 0 8, 1.28 W m K [31], 0.78 and 0.88 W m K for early and late stage of : / GLYPH&lt;1&gt; / GLYPH&lt;1&gt; hydration process [75], 0.53 W m K [32], and 1.0 and 1.07 W m K for fresh and after / GLYPH&lt;1&gt; / GLYPH&lt;1&gt; 28 days [76], our predicted values with MD simulation can be said to be well suited with them.

## 3.4. Temperature dependence of thermal conductivity

At the final stage, we investigated the thermal conductivity of porous jennite as increasing temperature from 240 K to 560 K with a step of 60 K for details and from 300 K to 1300 K with a step of 200 K. For the sake of simplicity without losing generality, we selected the three di GLYPH&lt;11&gt; erent typical porosities of 0%, 15.99% and 32.79%. We followed the above-mentioned procedure to calibrate the thermal conductivity: doing NPT simulation of perfect 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 jennite model for 100 ps, creating the pore with a radius corresponding to the porosity and running NVT for 400 ps, constructing simulation box by repeating 8 times the model and running NVT for 100 ps, and finally performing NVE run for 200 ps with imposing the heat flux, all these with a time step of 0.1 fs. These processes were repeated at di GLYPH&lt;11&gt; erent temperatures.

Table 6 lists the unit cell parameters and density of perfect jennite crystal obtained by doing NPT simulations as increasing temperature from 300 K to 1300 K. It was observed that as increasing temperature the unit cell volume increases monotonically while the mass density decreases, although the lattice constants varies irregularly. It should be noted that at 1300 K the density was too small to think of crystal; it was actually transformed to liquid state. Moreover, as revealed by Yu and Kirkpatrick [79],

the jennite crystal is known to be transformed to wollastonite and larnite above 1073 K. That is why we did not perform further simulation at 1300 K. Another point is that since the volume of jennite model increases with increasing temperature, the radius of spherical region for creating the pore should be adjusted to meet the target porosity.

Figure 7 shows the temperature profiles in the simulation boxes of perfect jennite along the principal crystallographic directions. At moderate temperatures, one can see a good linearity of temperature change along the slices, which makes it easy to measure the temperature gradients, while an inappreciable dispersion of temperatures was found at 1100 K.

In Fig. 8, we show the temperature dependence of thermal conductivities of porous jennite with porosities of 0, 15.99 and 32.79%. For the three cases of porosities, the principal thermal conductivities, especially volumetric averaged conductivities, increase to 500 K and decrease after that. Also, it was found that the [100] element of thermal conductivity became the biggest after 900 K in all the three cases. Figure 8(d) shows the volumetric thermal conductivities calculated in the detailed temperature range from 240 K to 560 K as well as in the temperature range from 300 K to 1100 K in comparison with our experimental values measured for the ceramic coating layer on the concrete wall by using the laser flash method. We note that our multifunctional composite building material was developed for use of thermal insulation, moisture barrier, sterilization of buildings [78]. Our experimental values were found to mostly locate between the porosity of 15.99% and 32.79%. Since the porosity of this ceramic coating layer was measured to be about 10% [77], these are slightly deviated

Table 6: The unit cell parameters and density of jennite crystal at di GLYPH&lt;11&gt; erent temperatures, obtained by MD simulation.

|                         |   300 K |   500 K |   700 K |   900 K |   1100 K |   1300 K |
|-------------------------|---------|---------|---------|---------|----------|----------|
| a (Å)                   |  10.865 |  11.6   |  12.691 |  13.666 |   12.647 |   14.238 |
| b (Å)                   |   7.402 |   7.627 |   7.188 |   7.246 |    6.183 |    7.312 |
| c (Å)                   |  10.946 |  10.474 |  11.014 |  10.832 |   15.031 |   25.811 |
| GLYPH<11> ( GLYPH<14> ) | 102.06  | 103.66  | 110.64  |  98.35  |   89.83  |   90.78  |
| GLYPH<12> ( GLYPH<14> ) |  94.57  |  93.95  |  84.53  |  93.47  |   99.12  |   89.64  |
| GLYPH<13> ( GLYPH<14> ) | 110.33  | 111.55  | 110.63  | 114.55  |  104.31  |  106.85  |
| V (Å 3 )                | 796.1   | 824.9   | 879.4   | 956.6   | 1123.6   | 2571.8   |
| GLYPH<26> (g / cm 3 )   |   2.218 |   2.141 |   2.008 |   1.846 |    1.633 |    0.687 |

from the predicted values, possibly due to not a pure jennite. Based on our calculation result, we suggested the 3rd-order polynomial relation between the thermal conductivity of porous jennite and temperature by

$$K ( T ) = c _ { 0 } + c _ { 1 } T + c _ { 2 } T ^ { 2 } + c _ { 3 } T ^ { 3 }$$

where ci is the coe GLYPH&lt;14&gt; cient. Table 7 shows them resulting from interpolation process and in Fig. 8(d) we see that the 3rd-order interpolation fits well with our calculation data.

Table 7: Coe GLYPH&lt;14&gt; cients of the 3rd-order polynomial relation between the thermal conductivity of porous jennite and temperature.

|   P (%) |   c 0 |   c 1 ( GLYPH<2> 10 GLYPH<0> 3 ) | c 2 ( GLYPH<2> 10 GLYPH<0> 7 )   | c 3 ( GLYPH<2> 10 GLYPH<0> 10 )   |
|---------|-------|----------------------------------|----------------------------------|-----------------------------------|
|    0    | 1.095 |                            0.145 | 7.510                            | GLYPH<0> 11.220                   |
|   15.99 | 0.445 |                            2.217 | GLYPH<0> 30.583                  | 9.635                             |
|   32.79 | 0.313 |                            1.535 | GLYPH<0> 22.750                  | 8.029                             |

We compared our calculation results with other experimental data for cement paste and concrete. Ukrainczyk et al . [80] conducted experiments with calcium aluminate cement with w c / = 0.3, 0.4, finding that its thermal conductivity increases with increasing temperature from 293 K to 353 K due to loss of water from jennite at 343 GLYPH&lt;24&gt; 363 K. At this temperature, the lattice constant c of jennite crystal was found to shrink and simultaneously the jennite transformed from Ca [Si O 9 6 18 (OH) ] 8H O to 6 GLYPH&lt;1&gt; 2 Ca 9 [Si 6 O16 (OH) ](OH) 2 8 GLYPH&lt;1&gt; 2H 2 O, called metajennite [79]. Xianzhi et al . [81] demonstrated that the thermal conductivity of cement paste decreased like 0.82, 0.78, 0.60 WmKat 333, 363 and 393 K. Wang / GLYPH&lt;1&gt; et al . [82] reported that the fly ash concrete had decreasing thermal conductivity from 1.69 to 0.95 W m K as increasing temperature / GLYPH&lt;1&gt; from 293 K to 793 K. As they pointed out, the change of microstructure was found to start at 573 K from SEM image observation, and Yu [79] stated that jennite lost the interlayer water molecules and thus basic layer structure upon heating to 623 K. It can be said that our calculation results are similar to these experimental data. On the contrary, some irregular changing tendency was observed from experiment for di GLYPH&lt;11&gt; erent kinds of concrete. For example, Gencel et al . [73] found that the thermal conductivity of concrete, made using vermiculite as dominant component, decreased from 293 K to 1193 K and again increased to 1373 K. Khaliq and Kodur [83] studied the temperature e GLYPH&lt;11&gt; ect on the thermal properties for self-consolidating concrete, finding that

its thermal conductivity decreased up to 693 K, marginally increased to 793 K, and again decreased to 1093 K. For aluminous refractory concrete, Santos [70] reported the increasing tendency from 300 K to 333 K, significantly decreasing to 793 K, and again gradual increasing to 1273 K. For explanation of such complicated behavior, the moisture loss by evaporation of free and pore water with the increase of temperature was considered as the main cause, which is di GLYPH&lt;11&gt; erent from our calculation models with the pore filled with air.

## 4. Conclusions

In this work, we have investigated the thermal conductivity of porous jennite by means of molecular dynamics method, aiming at clarifying the thermal behavior of cement hydration product C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H gel. The porous jennite models with di GLYPH&lt;11&gt; erent porosities have been created by removing atoms within spherical region with a radius corresponding to the porosity, using 3 GLYPH&lt;2&gt; 4 GLYPH&lt;2&gt; 3 supercell model for perfect jennite. We have employed two di GLYPH&lt;11&gt; erent MD methods, namely equilibrium Green-Kubo and nonequilibrium M¨ uller-Plathe methods, to calculate the thermal conductivity of jennite along the crystallographic [100], [010] and [001] directions. We found that the thermal conductivity of perfect jennite calculated by GK method using the supercell model were almost identical to those by MP method, exhibiting the crystallographic anisotropy like K 2 &gt; K 1 &gt; K 3 . The non-equilibrium MP method has been mostly applied in this work to investigate the porosity and temperature dependences of thermal conductivity of porous jennite, using the simulation boxes extended the porous jennite model 8 times along the heat flux direction and divided into 20 plane slices. At room temperature of 300 K, the principal and volumetric thermal conductivities of bulk jennite were determined to be K 1 = 1 100 : GLYPH&lt;6&gt; 0 025, : K 2 = 1 354 : GLYPH&lt;6&gt; 0 034, : K 3 = 0 967 : GLYPH&lt;6&gt; 0 020 and : K v = 1 141 : GLYPH&lt;6&gt; 0 026 W m K in reasonable agreement with the experimental values for : / GLYPH&lt;1&gt; C S GLYPH&lt;0&gt; GLYPH&lt;0&gt; H. We have found the decreasing tendency of e GLYPH&lt;11&gt; ective thermal conductivity of porous jennite as increasing the porosity, which was confirmed to follow the empirical coherent potential or self-consistent model proposed for porous materials among other kinds of models. Finally we have investigated the temperature dependence of thermal conductivity of porous jennite with the porosities of 0, 15.99 and 32.79%, finding its

gradual increase from 240 K to 560 K and decrease from 560 K to 1100 K. By comparing with our own measurements and available experimental data for cement paste and concrete, the calculated values were said to be high reasonable and reliable. Therefore, we believe this work will be helpful for a design of high thermal insulating materials with cement and concrete for energy e GLYPH&lt;14&gt; cient buildings.

## Acknowledgments

This work belongs to the state research project 'Design of Innovative Functional Materials for Energy and Environmental Application' (No. 2016-20), supported partially by the State Committee of Science and Technology, Democratic People's Republic of Korea. The simulations have been performed on the HP Blade System C7000 (HP BL460c) that is managed by Faculty of Materials Science, Kim Il Sung University.

## Notes

Declarations of interest: none.

## References

- [1] I. G. Richardson, The calcium silicate hydrates, Cem. Concr. Res. 38 (2008) 137158.
- [2] D. Gastaldi, F. Canonico, L. Capelli, E. Boccaleri, M. Milanesio, L. Palin, G. Croce, F. Marone, K. Mader, M. Stampanoni, In situ tomographic investigation on the early hydration behaviors of cementing systems, Constr. Build. Mater. 29 (2012) 284-290.
- [3] D. Fan, S. Yang, Mechanical properties of C-S-H globules and interfaces by molecular dynamics simulation, Construct. Build. Mater. 176 (2018) 573-582.
- [4] A. Al-Ostaz, W. Wu, A. H. D. Cheng, C. R. Song, A molecular dynamics and microporomechanics study on the mechanical properties of major constituents of hydrated cement, Composites Part B 41 (2010) 543-549.

- [5] A. Zaoui, Insight into elastic behavior of calcium silicate hydrated oxide (C-S-H) under pressure and composition e GLYPH&lt;11&gt; ect, Cem. Concr. Res. 42 (2012) 306-312.
- [6] S. Hajilar, B. Shafei, Nano-scale investigation of elastic properties of hydrated cement paste constituents using molecular dynamics simulations, Comput. Mater. Sci. 101 (2015) 216-226.
- [7] L. B. Skinner, S. R. Chae, C. J. Benmore, H. R. Wenk, P. J. M. Monteiro, Nanosturucture of calcium silicate hydrates in cements, Phys. Rev. Lett. 104 (2010) 195502.
- [8] E. Masoero, E. D. Gado, R. J. M. Pellenq, F. J. Ulm, S. Yip, Nanostructure and Nanomechanics of Cement: Polydisperse Colloidal Packing, Phys. Rev. Lett. 109 (2012) 155503.
- [9] I. G. Richardson, J. Skibsted, R. J. Kirkpatrick, Characterisation of cement hydrate phases by TEM, NMR and Raman spectroscopy, Adv. Cem. Res. 22 (2010) 233-248.
- [10] I. G. Richardson, The nature of C GLYPH&lt;0&gt; S GLYPH&lt;0&gt; H in hardened cements, Cem. Concr. Res. 29 (1999) 1131-1147.
- [11] H. F. W. Taylor, Proposed structure for calcium silicate hydrate gel, J. Am. Ceram. Soc. 69 (2005) 464-467.
- [12] R. J. M. Pellenq, A. Kushima, R. Shahsavari, K. J. V. Vliet, M. J. Buehler, S. Yip, et al., A realistic molecular model of cement hydrates, Proc. Natl. Acad. Sci. 106 (2009) 16102-16108.
- [13] A. J. Allen, J. J. Thomas, H. M. Jennings, Composition and density of nanoscale calcium-silicate-hydrate in cement, Nat. Mater. 6 (2007) 311-316.
- [14] D. Hou, T. Zhao, P. Wang, Z. Li, J. Zhang, Molecular dynamics study on the mode I fracture of calcium silicate hydrate under tensile loading, Eng. Fract. Mech. 131 (2014) 557-569.

- [15] D. Hou, Z. Li, T. Zhao, P. Zhang, Water transport in the nano-pore of the calcium silicate phase: reactivity, structure and dynamics, Phys. Chem. Chem. Phys. 17 (2015) 1411-1423.
- [16] M. J. A. Qomi, K. J. Krakowiak, M. Bauchy, K. L. Stewart, R. Shahsavari, D. Jagannathan, D. B. Wrommer, A. Baronnet, M. J. Buehler, S. Yip, F. J. Ulm, K. J. V. Vliet, R. J. M. Pellenq, Combinatorial molecular optimization of cement hydrates, Nat. Commun. 5 (2014) 4960.
- [17] A. Vidmer, G. Sclauzero, A. Pasquarello, Infrared spectra of jennite and tobermorite from first-principles, Cem. Concr. Res. 60 (2014) 11-23.
- [18] E. Bonaccorsi, S. Merlino, H. F. W. Taylor, The crystal structure of jennite, Ca9Si6(OH)6 GLYPH&lt;1&gt; 8H2O, Cem. Concr. Res. 34 (2004) 1481-1488.
- [19] S. Merlino, E. Bonaccorsi, T. Armbruster, The real structure of tobermorite 11 Å normal and anomalous forms, OD character and polytypic modifications, Eur. J. Mineral. 13 (2001) 577-590.
- [20] E. Bonaccorsi, S. Merlino, A. R. Kampf, The crystal structure of 14 Å(plombierite), a C-S-H phase, J. Am. Ceram. Soc. 88 (2005) 505-516.
- [21] C. Dharmawardhana, A. Misra, S. Aryal, P. Rulis, W. Ching, Role of interatomic bonding in the mechanical anisotropy and interlayer cohesion of CSH crystals, Cem. Concr. Res. 52 (2013) 123-130.
- [22] S. Churakov, Hydrogen bond connectivity in jennite from ab initio simulations, Cem. Concr. Res. 38 (2008) 1359-1364.
- [23] R. Shahsavari, M. J. Buehler, R. J. M. Pellenq, F. J. Ulm, First-principles study of elastic constants and interlayer interactions of complex hydrated oxides: case study of tobermorite and jennite, J. Am. Ceram. Soc. 92 (2009) 2323-2330.
- [24] J. Moon, S. Yoon, P. J. Monteiro, Mechanical properties of jennite: A theoretical and experimental study, Cem. Concr. Res. 71 (2015) 106-114.

- [25] D. Tunega, A. Zaoui, Understanding of bonding and mechanical characteristics of cementitious mineral tobermorite from first-principles, J. Comput. Chem. 32 (2011) 306-314.
- [26] P. Rejmak, J. S. Dolado, M. J. Stott, A. Ayuela, 29 Si NMR in Cement: A Theoretical Study on Calcium Silicate Hydrates, J. Phys. Chem. C 116 (2012) 9755-9761.
- [27] D. Hou, T. Zhao, Z. Jin, Z. Li, Structure, reactivity and mechanical properties of water ultra-confined in the ordered crystal: A case study of jennite, Micropor. Mesopor. Mater. 204 (2015) 106-114.
- [28] H. Manzano, J. S. Dolado, A. Ayuela, Elastic properties of the main species present in Portland cement pastes, Acta Mater. 57 (2009) 1666-1674.
- [29] S. M. Mutisya, J. M. de Almeida, C. R. Miranda, Molecular simulations of cement based materials: A comparison between first principles and classical force field calculations, Comput. Mater. Sci. 138 (2017) 392-402.
- [30] D. Hou, Z. Li, Molecular dynamics study of water and ions transport in nano-pore of layered structure: A case study of tobermorite, Micropor. Mesopor. Mater. 195 (2014) 9-20.
- [31] S. Yoon, D. E. Macphee, M. S. Imbabi, Estimation of the thermal properties of hardened cement paste on the basis of guarded heat flow meter measurements, Thermochim. Acta 588 (2014) 1-10.
- [32] Y. Xu, D. D. L. Chung, E GLYPH&lt;11&gt; ect of sand addition on the specific heat and thermal conductivity of cement, Cem. Concr. Res. 30 (2000) 59-61.
- [33] D. Bentz, Transient plane source measurements of the thermal properties of hydrating cement pastes, Mater. Construc. Mater. Structures 40 (2007) 1073-1080.
- [34] M. J. A. Qomi, F.-J. Ulm, R. J. M. Pellenq, Physical Origins of Thermal Properties of Cement Paste, Phys. Rev. Applied 3 (2015) 064010.

- [35] Z. Misri, M. H. W. Ibrahim, A. S. M. A. Awal, M. S. M. Desa, N. S. Ghadzali, Review on factors influencing thermal conductivity of concrete incorporating various type of waste materials, IOP Conf. Ser.: Earth Environ. Sci. 140 (2018) 012141.
- [36] G. M. Sabnis (Ed.), Green building with concrete, CRC Press (Taylor &amp; Francis Group), Boca Raton, 2012.
- [37] K. Gobakis, D. Kolokotsa, N. Maravelaki-Kalaitzaki, V. Perdikatsis, M. Santamouris, Development and analysis of advanced inorganic coatings for buildings and urban structures, Energy Build. 89 (2015) 196-205.
- [38] A. M. Papadopoulos, State of the art in thermal insulation materials and aims for future developments, Energy Build. 37 (2005) 77-86.
- [39] J. Jambor, Pore structure and strength development of cement composites, Cem. Concr. Res. 20 (1990) 948-954.
- [40] M. A. R. Bhutta, K. Tsuruta, J. Mirza, Evaluation of high-performance porous concrete properties, Construct. Build. Mater. 31 (2012) 67-73.
- [41] S. Hatanaka, N. Mishima, A. Maegawa, E. Sakamoto, Fundamental Study on Properties of Small Particle Size Porous Concrete, J. Adv. Concr. Technol. 12 (2014) 24-33.
- [42] R. Shahsavari, R. J. M. Pellenq, F. J. Ulm, Empirical force fields for complex hydrated calcio-silicate layered materials, Phys. Chem. Chem. Phys. 13 (2011) 1002-1011.
- [43] H.-Q. Jin, X.-L. Yao, L.-W. Fan, Z.-T. Yu, Experimental determination and fractal modeling of the e GLYPH&lt;11&gt; ective thermal conductivity of autoclaved aerated concrete: E ects of moisture content, Int. J. Heat Mass Transfer 92 (2016) 589-602. GLYPH&lt;11&gt;
- [44] R. K. Mishra, A. K. Mohamed, D. Geissb¨ uhler, H. Manzano, T. Jamil, R. Shahsavari, A. G. Kalinichev, S. Galmarini, L. Tao, H. Heinz, R. Pellenq, A. C. T. van

- Duin, S. C. Parker, R. J. Flatt, P. Bowen, cem GLYPH&lt;11&gt; : A force field database for cementitious materials including validations, applications and opportunities, Cem. Concr. Res. 102 (2017) 68-89.
- [45] D. Lau, W. Jian, Z. Yu, D. Hui, Nano-engineering of construction materials using molecular dynamics simulations: Prospects and challenges, Composits Part B 143 (2018) 282-291.
- [46] R. T. Cygan, J. J. Liang, A. G. Kalinichev, Molecular Models of Hydroxide, Oxyhydroxide, and Clay Phases and the Development of a General Force Field, J. Phys. Chem. B 108 (2004) 1255-1266.
- [47] H. Heinz, T. J. Lin, R. K. Mishra, F. S. Emami, Thermodynamically consistent force fields for the assembly of inorganic, organic, and biological nanostructures: the INTERFACE force field, Langmuir 29 (2013) 1754-1765.
- [48] C. L. Freeman, J. H. Harding, D. J. Cooke, J. A. Elliott, J. S. Lardge, D. M. Du GLYPH&lt;11&gt; y, New forcefields for modeling biomineralization processes, J. Phys. Chem. C 111 (2007) 11943-11951.
- [49] A. C. V. Duin, S. Dasgupta, F. Lorant, W. A. Goddard, ReaxFF: a reactive force field for hydrocarbons, J. Phys. Chem. A 105 (2001) 9396-9409.
- [50] A. G. Kalinichev, J. W. Wang, R. J. Kirkpatrick, Molecular dynamics modeling of the structure, dynamics and energetics of mineral-water interfaces: application to cement materials, Cem. Concr. Res. 37 (2007) 337-347.
- [51] A. G. Kalinichev, R. J. Kirkpatrick, Molecular dynamics simulation of cationic complexation with natural organic matter, Eur. J. Soil Sci 58 (2007) 909-917.
- [52] D. Tavakoli, A. Tarighat, Molecular dynamics study on the mechanical properties of portland cement clinker phases, Comput. Mater. Sci. 119 (2016) 65-73.
- [53] S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, J. Comput. Phys. 117 (1995) 1-19.

- [54] G. M. Bowers, R. J. Kirkpatrick, Natural abundance 43 Ca NMR spectroscopy of tobermorite and jennite: Model compounds for C-S-H, J. Am. Ceram. Soc. 92 (2009) 545-548.
- [55] M. S. Green, Marko GLYPH&lt;11&gt; random processes and the statistical mechanics of timedependent phenomena. II. Irreversible processes in fluids, J. Chem. Phys. 22 (1954) 398-413.
- [56] R. Kubo, M. Yokota, S. Nakajima, Statistical-mechanical theory of irreversible processes. II. Response to thermal disturbance, J. Phys. Soc. Jpn 12 (1957) 12031211.
- [57] F. M¨ uller-Plathe, A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity, J. Chem. Phys. 106 (1997) 6082-6085.
- [58] D. Frenkel, B. Smit, Understanding Molecular Simulation From Algorithms to Applications, Academic Press, Florida, USA, 2002.
- [59] H. Gu, H. Wang, E GLYPH&lt;11&gt; ect of strain on thermal conductivity of amorphous silicon dioxide thin films: A molecular dynamics study, Comput. Mater. Sci. 144 (2018) 133-138.
- [60] S. S. Mahajan, G. Subbarayan, B. G. Sammakia, Estimating thermal conductivity of amorphous silica nanoparticles and nanowires using molecular dynamics simulations, Phys. Rev. E 76 (2007) 056701.
- [61] T. Coquil, J. Fang, L. Pilon, Molecular dynamics study of the thermal conductivity of amorphous nanoporous silica, Int. J. Heat Mass Transfer 54 (2011) 45404548.
- [62] W. Zhu, G. Zheng, S. Cao, H. He, Thermal conductivity of amorphous SiO2 thin film: A molecular dynamics study, Sci. Rep. 8 (2018) 10537.
- [63] T. Wu, I. Temizer, P. Wriggers, Computational thermal homogenization of concrete, Cem. Concr. Compos. 35 (2013) 59-70.

- [64] J. C. M. Garnett, Colours in metal glasses and in metallic films, Phil. Trans. Roy. Soc. Lond. A 203 (1904) 385-420.
- [65] R. Landauer, The Electrical Resistance of Binary Metallic Mixtures, A. Appl. Phys. 23 (1952) 779.
- [66] D. G. Cahill, T. H. Allen, Thermal conductivity of sputtered and evaporated SiO2 and TiO2 optical coatings, Appl. Phys. Lett. 65 (1994) 309-311.
- [67] H. W. Russell, Principles of Heat Flow in Porous Insulators, J. Am. Ceram. Soc. 18 (1935) 1.
- [68] S. Frey, Uber Die Elektrische Leitfahigkeit Binarer Aggregate, Zeitschrift Fur Elektrochemie 38 (1932) 260-274.
- [69] D. R. Lide, ed., CRC Handbook of Chemistry and Physics, CRC Press Taylor and / Francis, Boca Raton, FL., 2009.
- [70] W. N. D. Santos, E GLYPH&lt;11&gt; ect of moisture and porosity on the thermal properties of a conventional refractory concrete, J. Eur. Ceram. Soc. 23 (2003) 745-755.
- [71] H. Awang, etc., E ect of additives GLYPH&lt;11&gt; on mechanical and thermal properties of lightweight foamed concrete, Adv. Appl. Sci. Res. 3 (2012) 3326-3338.
- [72] F. Batool, E GLYPH&lt;11&gt; ect of Microstructure on Thermal Conductivity of Cement-Based Foam, 2015.
- [73] O. Gencel, F. Koksal, M. Sahin, M. Y. Durgun, H. E. H. Lobland, W. Brostow, Modeling of Thermal Conductivity of Concrete with Vermiculite using by Artificial Neural Networks Approaches, Exp. Heat Transf. 26 (2013) 360-383.
- [74] G. Constantinides, F.-J. Ulm, The nanogranular nature of C-S-H, J. Mech. Phys. Solids 55 (2007) 64-90.
- [75] P. F. Hansen, J. Hansen, K. Hougaard, E. J. Pedersen, Thermal properties of hardening cement paste, in: Proceeding of RILEM international conference on concrete at early ages, RILEM, Paris, 23-36, 1982.

- [76] P. Mounanga, A. Khelid, G. Bastian, Experimental study and modeling approaches for the thermal conductivity evolution of hydrating cement paste, Adv. Cem. Res. 16 (2004) 95-103.
- [77] C.-J. Yu, B.-H. Ri, C.-H. Kim, U.-S. Hwang, K.-C. Ri, C.-J. Song, U.-C. Kim, Formation and characterization of ceramic coating from alumino silicate mineral powders in the matrix of cement composite on the concrete wall, Mater. Chem. Phys. 227 (2019) 211-218.
- [78] B.-H. Ri, K. K. Ri, D. K. Ri, Multifunctional composite building materials and construction thereof, WIPO Patent WO 2018 124314 A1, 2018. /
- [79] P. Yu, R. Kirkpatrick, Thermal dehydration of tobermorite and jennite, Concr. Sci. Eng. 1 (1999) 185-191.
- [80] N. Ukrainczyk, T. Matusinovi´ c, Thermal properties of hydrating calcium aluminate cement pastes, Cem. Concr. Res. 40 (2010) 128-136.
- [81] S. Xianzhi, Z. Rui, L. Ruixia, L. Gensheng, S. Baojiang, S. Yu, W. Gaosheng, Z. Shijie, Study on thermal conductivity of cement with thermal conductive materials in geothermal well, Geothermics 81 (2019) 1-11.
- [82] W. Wang, C. Lu, Y. Li, Q. Li, An investigation on thermal conductivity of fly ash concrete after elevated temperature exposure, Construct. Build. Mater. 148 (2017) 148-154.
- [83] W. Khaliq, V. Kodur, Thermal and mechanical properties of fi ber reinforced high performance self-consolidating concrete at elevated temperatures, Cem. Concr. Res. 41 (2011) 1112-1122.

Figure 3: (a) The simulation box created by repeating the supercell (porosity 24.4%) eight times in [100] direction. (b) Sketch of the simulation box divided into 20 slices, of which the first and 11th slices are for heat sink and source. Temperature profiles as increasing the system length (number of repeated jennite cells N Re) in (c) [100], (d) [010], and (e) [001] directions at 300 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-00b06aef-1c44-496f-a88a-65ace3abd280.jpg)

Figure 4: Thermal conductivity of perfect jennite as increasing the system length, obtained by applying the non-equilibrium M¨ uller-Plathe method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-4fffc457-4093-411a-b3d7-9e780a31e191.jpg)

Figure 5: Normalized thermal conductivity of porous jennite as a function of porosity obtained from the suggested models and MD simulation in this work.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-b85f2217-3c04-4ce8-80e3-b7fffcf5908f.jpg)

Figure 6: Thermal conductivity of cement hydration products including jennite with their pores filled with air or water, comparing with experimental data. Exp1, 2, 3, 4 are from Ref. [43, 71, 72, 73]. They are well placed between the blue line for air filling and the red line for water filling estimated by using the CP model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-b00085e2-dd5d-4248-80a2-4d657ebe6f5d.jpg)

Figure 7: Temperature profiles in simulation boxes of perfect jennite crystal along (a) [100], (b) [010] and (c) [001] directions as increasing temperature from 300 K to 1100 K with a step of 200 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-684e293f-5dd4-4352-af4f-97ffa5fb07ff.jpg)

Figure 8: Thermal conductivities of porous jennite with various porosities of (a) 0%, (b) 15.99%, and (c) 32.79% along the [100], [010] and [001] directions as increasing temperature from 300 K to 1100 K with a step of 200 K, where Vol means the volumetric averaged value. (d) Volumetric thermal conductivities as functions of temperature, merged those from 240 K to 560 K with a step of 60 K and those from 300 K to 1100 K with a step of 200 K, in comparison with our measurement (Exp1, 2) performed by using the laser flash method. Solid lines indicate the 3rd-order polynomial interpolation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00402v1-figure-c42780c2-6a2a-4767-8abf-18a095a07184.jpg)
Song-Nam Hong Chol-Jun Yua Un-Song Hwang Chung-Hyok Kim Byong-Hyok Ri
Thermal conductivity of porous jennite by molecular dynamics method
The thermal conductivity of porous jennite, as the major component of cement
paste, and its porosity and temperature dependences are simulated by molecular
dynamics methods using ClayFF force field. The porous jennite models with
different porosities are created by removing atoms within the sphere from bulk
jennite model. The thermal conductivity elements of bulk jennite calculated by
equilibrium Green-Kubo method are almost identical to those by non-equilibrium
M\"{u}ller-Plathe (MP) method at 300 K. The volumetric thermal conductivity of
porous jennite is found to decrease from 1.141 W/m$\cdot$K to 0.144 W/m$\cdot$K
as increasing the porosity from 0\% to 72.22\% at 300 K, following the
empirical coherent potential model when the pore is assumed to be filled with
air. As increasing temperature, the thermal conductivity is observed to
increase from 240 K to 560 K and gradually decrease until 1100 K for porous
jennites with the porosities of 0\%, 15.99\% and 32.79\%.