Thermal Conductivity of Silicon Thin Films Predicted by Molecular Dynamics Simulations and Theoretical Calculation

2011 ◽  
Vol 55-57 ◽  
pp. 1152-1155 ◽  
Author(s):  
Xing Li Zhang ◽  
Zhao Wei Sun

Molecular, dynamics simulation and the Boltzmann transport equation are used respectively to analyze the phonon transport in Si thin film. The MD result is in good agreement with the theoretical analysis values. The results show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced and is almost independent of the temperature at the nanoscale. It was observed from the simulation results that there exists the obvious size effect on the thermal conductivity.

Author(s):  
J. A. Thomas ◽  
R. M. Iutzi ◽  
A. J. H. McGaughey

The effective thermal conductivity of water/carbon nanotube (CNT) composite systems is predicted using molecular dynamics simulation. Both empty and water-filled CNTs with diameters ranging from 0.83 nm to 1.26 nm are considered. Using a direct application of the Fourier law, we explore the transition to diffusive phonon transport with increasing CNT length and identify the correlation between CNT diameter and fully-diffusive thermal conductivity. Using Green-Kubo linear response theory, we explore how the thermal conductivity of water inside CNT varies with tube diameter. We predict the effective thermal conductivity of the composite systems and examine how the phonon modes in the CNT are affected by interactions with the water molecules.


Author(s):  
Toshihiro Kaneko ◽  
Kenji Yasuoka ◽  
Ayori Mitsutake ◽  
Xiao Cheng Zeng

Multicanonical molecular dynamics simulations are applied, for the first time, to study the liquid-solid and solid-solid transitions in Lennard-Jones (LJ) clusters. The transition temperatures are estimated based on the peak position in the heat capacity versus temperature curve. For LJ31, LJ58 and LJ98, our results on the solid-solid transition temperature are in good agreement with previous ones. For LJ309, the predicted liquid-solid transition temperature is also in agreement with previous result.


Author(s):  
Bo Qiu ◽  
Hua Bao ◽  
Xiulin Ruan

In this paper, thermoelectric properties of bulk PbTe are calculated using first principles calculations and molecular dynamics simulations. The Full Potential Linearized Augmented Plane Wave (FP-LAPW) method is first employed to calculate the PbTe band structure. The transport coefficients (Seebeck coefficient, electrical conductivity, and electron thermal conductivity) are then computed using Boltzmann transport equation (BTE) under the constant relaxation time approximation. Interatomic pair potentials in the Buckingham form are also derived using ab initio effective charges and total energy data. The effective interatomic pair potentials give excellent results on equilibrium lattice parameters and elastic constants for PbTe. The lattice thermal conductivity of PbTe is then calculated using molecular dynamics simulations with the Green-Kubo method. In the end, the figure of merit of PbTe is computed revealing the thermoelectric capability of this material, and the multiscale simulation approach is shown to have the potential to identify novel thermoelectric materials.


2021 ◽  
Author(s):  
Zhi Meng Zhang ◽  
Hua Yang ◽  
Jun Xia Shi ◽  
Jia Jun Wang ◽  
Zheng Guo Huang ◽  
...  

Abstract The orientation of hydrocarbon chains plays a key role in the applications of organic materials. And chain folding in the process of molecular orientation is also of great significance for the design of organic molecular thin films. The effect of chain length and simulation temperature on the isothermal orientation of n-alkanes on graphene surface is studied by molecular dynamics simulation in this paper. And the chain folding is also described. The n-alkanes can form perpendicular ordered structure, parallel ordered structure or perpendicular orientation at relative low temperature and parallel orientation at relative high temperature on graphene surface. The chain fold happens when long n-alkanes form perpendicular ordered structure on graphene surface. And the simulation results show the interactions of n-alkane−graphene and n-alkane−n-alkane affect chain fold.


2012 ◽  
Vol 26 (20) ◽  
pp. 1250117 ◽  
Author(s):  
L. T. VINH ◽  
N. V. HUY ◽  
P. K. HUNG

Molecular dynamics simulation is carried out for liquid SiO 2 at pressure ranged from zero to 30 GPa and by using BKS, Born–Mayer type and Morse–Stretch potentials. The constructed models reproduce well the experimental data in terms of mean coordination number, bond angle and pair radial distribution function. Furthermore, the density of all samples can be expressed by a linear function of fractions SiO x. It is found that the topology of units SiO x and linkages OSi y is unchanged upon compression although the liquid undergoes substantial change in its network structure. Consequently, the partial bond angle distribution for SiO x and OSi y is identical for all samples constructed by the same potential. This result allows to establishing a simple expression between total bond angle distribution (BAD) and fraction of SiO x and OSi y. The simulation shows a good agreement between the calculation and simulation results for both total O–Si–O and Si–O–Si BADs. This supports a technique to estimate amount of units SiO x and linkages OSi y on base of total Si–O–Si and O–Si–O BADs measured experimentally.


1998 ◽  
Vol 545 ◽  
Author(s):  
G. Chen ◽  
S. G. Volz ◽  
T. Borca-Tasciuc ◽  
T. Zeng ◽  
D. Song ◽  
...  

AbstractUnderstanding phonon heat conduction mechanisms in low-dimensional structures is of critical importance for low-dimensional thermoelectricity. In this paper, we discuss heat conduction mechanisms in two-dimensional (2D) and one-dimensional (1D) structures. Models based on both the phonon wave picture and particle picture are developed for heat conduction in 2D superlattices. The phonon wave model, based on the acoustic wave equations, includes the effects of phonon interference and tunneling, while the particle model, based on the Boltzmann transport equation, treats the internal as well interface scattering of phonons. For 1D systems, both the Boltzmann transport equation and molecular dynamics simulation approaches are employed. Comparing the modeling results with experimental data suggest that the interface scattering of phonons plays a crucial role in the thermal conductivity of low-dimensional structures. We also discuss the minimum thermal conductivity of low-dimensional structures based on a generalized thermal conductivity integral, and suggest that the minimum thermal conductivities of low-dimensional systems may differ from those of their corresponding bulk materials. The discussion leads to alternative ways to reduce thermal conductivity based on the propagating phonon modes.


2020 ◽  
Vol 142 (3) ◽  
Author(s):  
Mohammad Bagheri Motlagh ◽  
Mohammad Kalteh

Abstract In this paper, molecular dynamics simulation is used to investigate the effect of copper and argon nanochannels size on the thermal conductivity of argon. Thermal conductivity is calculated by nonequilibrium molecular dynamics (NEMD) simulation. Simulations are performed for different distances between the walls. Results for both copper and argon walls are investigated individually. Results show that the existence of argon walls has little effect on the thermal conductivity. However, the amount of it for the argon confined between the copper walls is affected by the distance between the two walls. In the same way, the effect of wall roughness on the thermal conductivity is investigated, which shows that roughness is effective only for low distances between the walls. Also, the thermal conductivity of argon under Poiseuille flow in a nanochannel is studied. The results indicate that by increasing the driving force, the thermal conductivity increases and the increase ratio is higher for larger forces.


Author(s):  
G. T. Gao ◽  
J. D. Schall ◽  
K. Van Workum ◽  
P. T. Mikulski ◽  
J. A. Harrison

A constant tension and constant temperature molecular dynamics simulation method was used in the calculations of the elastic constants of the nanocomposite systems. The nanocomposite systems contain a core of sp3 diamond structure surrounded by an amorphous carbon network. The simulation results show that the elastic properties of nanocomposites of diamond-like carbons are closely related to the size of the sp3 diamond core; the bigger the core, the larger the elastic constants, and the system becomes more anisotropic.


2012 ◽  
Vol 501 ◽  
pp. 64-69 ◽  
Author(s):  
Yan He ◽  
Yuan Zheng Tang ◽  
Man Ding ◽  
Lian Xiang Ma

Normal thermal conductivity of amorphous and crystalline SiO2nano-films is calculated by nonequilibrium molecular dynamics (NEMD) simulations in the temperature range from 100 to 700K and thicknesses from 2 to 6nm. The calculated temperature and thickness dependences of thermal conductivity are in good agreement with previous literatures. In the same thickness, higher thermal conductivity is obtained for crystalline SiO2nano-films. And more importantly, for amorphous SiO2nano-films, thickness can be any direction of x, y, z-axis without effect on the normal thermal conductivity, for crystalline SiO2nano-films, the different thickness directions obtain different thermal conductivity results. The different results of amorphous and crystalline SiO2nano-films simply show that film thickness and grain morphology will cause different effects on thermal conductivity.


Sign in / Sign up

Export Citation Format

Share Document