Phonon thermal conductivity in nanolaminated composite metals via molecular dynamics

We use nonequilibrium molecular dynamics to characterize the phonon contribution to thermal conduction of Al nanostructures and the role of interfaces in metallic nanocomposites. We characterize the lattice thermal conductivity of pure Al samples as a function of size and temperature from which we obtain, using kinetic theory, the temperature dependence of the phonon mean free path. We also calculated the thermal conductivity of AlAl* and AlNi nanolaminate composites (where Al* differs from Al only in its mass) for various periodic sizes and compositions as well as the associated interfacial thermal resistivities (ITRs). We find that simple, additive models provide good estimates of the thermal conductivities of the nanocomposites in terms of those of the individual components and interfaces if size effects on the behavior of the individual components are considered. The additive models provide important insight to the decrease in thermal conductivity of the nanolaminates as their periodicity (thickness of a bilayer) is reduced to a size comparable with the phonon mean free path and break down when this characteristic size is reduced further. At this point the system can be regarded as homogeneous and the conductivity increases with decreasing periodicity of the laminates. We also observe that the ITR depends on the direction of the heat flux; this is the first molecular level characterization of such thermal diode behavior in a realistic three dimensional material.     

A new approach to lattice thermal conductivity: Application to Ge

The new expression τ _3ph ^?1 =g(ω) (B_N + B_UE^?Θ/αT)T^m is proposed for the three-phonon scattering relaxation rate, considering contributions due to three-phonon normal and umklapp processes, which give a new approach to the lattice thermal conductivity. With use of the above expression, the lattice thermal conductivity of Ge has been calculated in the entire temperature range 2–1000 K: good agreement is found between the experimental and calculated values of the phonon conductivity in the entire temperature range of investigation. Analytical expressions are also obtained to calculate an approximate value of the lattice thermal conductivity. The role of four-phonon processes is also included in the present study.     

Thermal conductivity of molten alkali halides from equilibrium molecular dynamics simulations

The thermal conductivity of molten sodium chloride and potassium chloride has been computed through equilibrium molecular dynamics Green-Kubo simulations in the microcanonical ensemble (N,V,E). In order to access the temperature dependence of the thermal conductivity coefficient of these materials, the simulations were performed at five different state points. The form of the microscopic energy flux for ionic systems whose Coulombic interactions are calculated through the Ewald method is discussed in detail and an efficient formula is used by analogy with the methods used to evaluate the stress tensor in Coulombic systems. The results show that the Born-Mayer-Huggins-Tosi-Fumi potential predicts a weak negative temperature dependence for the thermal conductivity of NaCl and KCl. The simulation results are in agreement with part of the experimental data available in the literature with simulation values generally overpredicting the thermal conductivity by 10%-20%.     

Thermal transport coefficients for liquid and glassy water computed from a harmonic aqueous glass

We compute thermal transport coefficients for liquid and glassy water in terms of the vibrations of the quenched liquid. The thermal conductivity and thermal diffusivity are computed for H(2)O and D(2)O at densities from 0.93 to 1.2 g cm(-3). The computed thermal diffusivity of liquid water is in reasonable agreement with measured values and is found to increase with increasing temperature due largely to the thermal accessibility of delocalized librational modes. The influence of structure and density on the thermal conductivity of amorphous ices is investigated. The calculations reveal that density alone is unable to explain the measured thermal conductivity of amorphous ices, particularly low-density amorphous ices, for which the thermal conductivity decreases with increasing temperature near 100 K. To investigate the influence of structure on thermal transport in amorphous ices we have computed the thermal transport coefficients for low-density amorphous ices prepared in two different ways, one formed by quenching the liquid at 0.93 g cm(-3) and the other by distortion of cubic ice at the same density. The computed thermal conductivity of the latter is higher, but the structures of both forms are too disordered for the thermal conductivity to exhibit the unusual variation observed experimentally.     

Transport Properties of Simple Fluids from the Soft‐Core Double Yukawa Potential

The present work is concerned with checking a new and simple pair potential function (soft‐core double Yukawa potential) for noble gases by calculation of the transport properties. The viscosity, thermal conductivity and self diffusion coefficient in dilute gas limit in the temperature range of 298‐1400 K are calculated and agreement with the measurements is, in general, within the experimental error. A comparison of the calculated and experimental values of the viscosity, thermal conductivity and the diffusion coefficients yields an average absolute deviation of 0.5%, 1.5% and 1.2%, respectively. Also, the calculated transport properties from this potential have been compared with those calculations via the accurate experimental potential and also the corresponding state.     

Thermal conductivity of wurtzite and zinc blende cubic phases of BeO from ab initio calculations

The structural, mechanical, thermal and thermodynamic properties of Beryllium oxide (BeO) in the zinc blende (ZB) and wurtzite (WZ) form have been calculated using the density functional theory (DFT) in the general gradient approximation (GGA). The ground state structural and elastic properties of wurtzite BeO (w-BeO) is calculated using the new GGA ultrasoft pseudopotentials for solids (pbesol); the simulated results have shown excellent agreement with the experiments. The thermodynamic properties are studied using quasi-harmonic approximation (QHA), and the predicted properties agree well for the WZ phase for which the experimental data are available, while for ZB phase it remains to be validated with future experiments. Both Boltzmann transport equation (BTE) and Slack model were used to calculate the lattice thermal conductivity of wurtzite BeO (w-BeO). Furthermore, the thermal conductivity along the crystallographic ‘a’ and ‘c’ axis of wurtzite BeO is investigated using BTE. Our calculation of w-BeO agrees well with the available experimental measurements. Apart from these studies on w-BeO, we have also compared the mechanical, structural and phonon dispersions of z-BeO with previously reported theoretical studies. Additionally we report the volume thermal expansion and the heat capacity at constant pressure of z-BeO for the first time and the bulk thermal conductivity of zinc blende BeO (z-BeO) using BTE.     

Elastic, electronic and thermal properties of YSZ from first principles

First principles calculations were performed to investigate the elastic, electronic and thermal properties of 14% cubic yttria-stabilized zirconia (YSZ) using the pseudo potential plane-wave method within the gradient generalized approximation (GGA) for the exchange and correlation potential. Computed lattice constant parameters are in good agreement with the available experimental results. The three independent elastic constants were computed by means of the stress-strain method, indicating that 14% cubic YSZ is a mechanically stable structure. From the knowledge of the elastic constants, a set of related properties, namely bulk, shear modulus, Young’s modulus, sound velocity, Debye temperature, thermal capacity and minimum thermal conductivity are numerically estimated in the frame work of the Voigt-Reuss-Hill approximation for YSZ polycrystalline. The calculated bulk modulus, shear modulus, Young’s modulus, sound velocity, Debye temperature, thermal capacity and minimum thermal conductivity are in reasonable agreement with the available experimental and theory data. Density of states, charge density and Mulliken population analysis show that the 14% cubic YSZ is covalent and possess ionic character.     

Phonon Scattering by Externded Defects in Polycrystalline Graphene

Phonon scattering by static stress fields created by grain boundaries (GBs) in graphene was studied with the deformation potential method. This method provides exact analytical expressions for phonon mean free paths as a result of scattering by GBs with arbitrary geometries in the Born approximation. Two regimes, k^–1 and k^–3, of the mean free path behavior are discovered in the limit of small momenta. The first (dislocation) regime is realized for open configurations in polycrystalline graphene, while the second regime is realized for closed configurations and Stone-Wales defects in pure graphene. The mean free path does not depend on k for any GB in the short-wave limit (large k). Thermal conductivity in graphene with GBs was calculated using the Callaway model which considers both normal phonon processes and umklapp processes. The contribution from these defects to the thermal conductivity is shown to exceed substantially the contribution from point defects and vacancies in a wide range of temperatures.     

Tailoring thermal conductivity of bulk graphene oxide by tuning the oxidation degree

The thermal conductivity of graphene oxides can be tailored by tuning oxidation degree due to the introduction of atomic- and nano-scale phonon scattering centers.     

Thermodynamic and transport properties of carbon dioxide from molecular simulation

Monte Carlo and molecular dynamics simulations have been used in order to test the ability of a three center intermolecular potential for carbon dioxide to reproduce literature experimental thermophysical values. In particular, both the shear viscosity under supercritical conditions and along the phase coexistence line, as well as the thermal conductivity under supercritical conditions, have been calculated. Together with the already reported excellent agreement for the phase coexistence densities, the authors find that the agreement with experimental values is, in general, good, except for the thermal conductivity at low density. Although extended versions of the model were employed, which include an explicit account of bending and vibrational degrees of freedom, a significant difference was still found with respect to the reported experimental value.     

Energy decay mechanisms and anharmonic lattice dynamics: The case of solid nitrogen

The anharmonic frequencies and linewidths of the lattice phonons in -N₂ are calculated on the basis of three different intermolecular potentials which include atom-atom and electrostatic interactions. The distinction between stationary anharmonicity and decay anharmonicity is stressed and the mechanism of energy transfer between the optical lattice phonons and the two-phonon manifold of the crystal is discussed in detail. The temperature dependence of the phonon self-energy is also considered. The results thus obtained for -N₂ are compared with predictions from previous lattice dynamics. SCP and molecular dynamics calculations. The calculated anharmonic effects are substantially independent of the adopted potential: the agreement with experimental data is reasonably good as far as the linewidths are concerned, while the anharmonic deformation of the potential wells (and thus the frequency shifts) is overestimated. We suggest that, while higher orders in the diagram expansion are necessary for a proper account of the stationary anharmonicity, the decay anharmonicity limits its effectiveness to two-phonon processes, thus allowing proper predictions of the phonon lifetimes by using the lowest-order diagrams. Finally, -N₂ is compared to -CO, and the role played by the translation-rotation coupling is discussed.     

氮化硼纳米管热输运性能的分子动力学模拟

采用基于声子散射理论的Boltzmann-Peierls声子传输方程(BTE)和非平衡态分子动力学模拟(NEMD)方法研究了氮化硼纳米管(BNNT)的热输运性能.分析了BNNT的热力耦合效应,通过BTE与NEMD两种方法相结合,分析了温度和长度对BNNT热输运性能的影响,并应用量子修正扩大了NEMD的研究范围.结果表明:随着拉伸或压缩应变的增加,BNNT热输运性能均呈降低的趋势.通过计算声子态密度(PDOS)在理论上分析了以上结果,发现在拉伸状态下,声子模式的变化是决定BNNT热输运性能变化的主要因素;在压缩状态下,热导率变化是由于模型发生明显的屈曲变形引起的.在低温段,BNNT的热输运性能受量子效应影响最初有一个线性增加的过程,当温度超过一定值时,其开始显著地降低;当BNNT长度小于120nm时,随着长度的增加,其弹道性能逐渐减弱,但仍主要体现为弹道特征,其热导率(κ)与长度(L)基本满足κ∝Lα这一关系.     

Experimental and theoretical study of two new pyrazoline derivatives based on dibenzofuran

Two novel pyrazoline derivatives, named 2,8-bis(1,3-diphenyl-pyrazoline-5-yl)dibenzofuran (A) and 2,8-bis(1-(4-bromophenyl)-3-phenyl-pyrazoline-5-yl)dibenzofuran (B), were synthesized and characterized by elemental analysis, NMR, MS and thermogravimetric analysis. The absorption and emission spectra of them were determined by experimental methods in different polar solvents and were computed using the density functional theory (DFT) and the time-dependent density functional theory (TDDFT) at the same time. The calculated absorption and emission wavelengths are in good agreement with the experimental data. The fluorescence quantum yields and fluorescence lifetimes of them in different polar solvents were studied by means of steady state and time resolved fluorescence. The calculated reorganization energy for hole and electron indicates that the two compounds are in favor of hole transport than electron transport. The results show the two compounds present high fluorescence quantum yields and excellent thermal stability. It makes them of great interest as novel fluorescent probes and optoelectronic materials.     

Relativistic effective valence shell Hamiltonian method: excitation and ionization energies of heavy metal atoms

The relativistic effective valence shell Hamiltonian H(v) method (through second order) is applied to the computation of the low lying excited and ion states of closed shell heavy metal atoms/ions. The resulting excitation and ionization energies are in favorable agreement with experimental data and with other theoretical calculations. The nuclear magnetic hyperfine constants A and lifetimes tau of excited states are evaluated and they are also in accord with experiment. Some of the calculated quantities have not previously been computed.     

Modifying thermal transport in electrically conducting polymers: effects of stretching and combining polymer chains

If their thermal conductivity can be lowered, polyacetylene (PA) and polyaniline (PANI) offer examples of electrically conducting polymers that can have potential use as thermoelectrics. Thermal transport in such polymers is primarily influenced by bonded interactions and chain orientations relative to the direction of heat transfer. We employ molecular dynamics simulations to investigate two mechanisms to control the phonon thermal transport in PANI and PA, namely, (1) mechanical strain and (2) polymer combinations. The molecular configurations of PA and PANI have a significant influence on their thermal transport characteristics. The axial thermal conductivity increases when a polymer is axially stretched but decreases under transverse tension. Since the strain dependence of the thermal conductivity is related to the phonon scattering among neighboring polymer chains, this behavior is examined through Herman's orientation factor that quantifies the degree of chain alignment in a given direction. The conductivity is enhanced as adjacent chains become more aligned along the direction of heat conduction but diminishes when they are orthogonally oriented to it. Physically combining these polymers reduces the thermal conductivity, which reaches a minimum value for a 2:3 PANI/PA chain ratio.     

Transport and Helfand moments in the Lennard-Jones fluid. II. Thermal conductivity

The thermal conductivity is calculated with the Helfand-moment method in the Lennard-Jones fluid near the triple point. The Helfand moment of thermal conductivity is here derived for molecular dynamics with periodic boundary conditions. Thermal conductivity is given by a generalized Einstein relation with this Helfand moment. The authors compute thermal conductivity by this new method and compare it with their own values obtained by the standard Green-Kubo method. The agreement is excellent.     

Extensive theoretical studies on the low-lying electronic states of indium monochloride cation, InCl+

The global potential energy curves for the 14 low-lying doublet and quartet Lambda-S states of InCl+ are calculated at the scalar relativistic MR-CISD+Q (multireference configuration interaction with single and double excitations, and Davidson's correction) level of theory. Spin-orbit coupling is accounted for via the state interaction approach with the full Breit-Pauli Hamiltonian, which leads to 30 Omega states. The computed spectroscopic constants of nine bound Lambda-S states and 17 bound Omega states are in good agreement with the available experimental data. The transition dipole moments and Franck-Condon factors of selected transitions are also calculated, from which the corresponding radiative lifetimes are derived.     

Ab initio vibrational predissociation dynamics of He-I2(B) complex

Three-dimensional quantum mechanical calculations on the vibrational predissociation dynamics of HeI2 B state complex are performed using a potential energy surface accurately fitted to unrestricted open-shell coupled cluster ab initio data, further enabling extrapolation for large I2 bond lengths. A Lanczos iterative method with an optimized complex absorbing potential is used to determine energies and lifetimes of the vibrationally predissociating He,I2(B,v') complex for v'相似文献