Molecular dynamics simulations of heat conduction in multi-walled carbon nanotubes

Heat conduction in multi-walled carbon nanotubes (MWNTs) was studied using non-equilibrium molecular dynamics simulations. This research focuses on the effects of the multi-wall structure of the MWNTs on the heat conduction. The results show that the thermal conductivity of a MWNT is almost the same as that of the corresponding single-walled carbon nanotubes (SWNTs) rather than much smaller as has been suggested. Thus, the multi-wall structure does not significantly affect the thermal conduction in the MWNTs. Analysis of the temperature profiles and the phonon density of states confirms that there is almost no heat transport between the MWNT layers and that each layer conducts heat nearly independently along parallel channels. This is physically reasonable since the weak inter-wall interactions and large interfacial thermal resistances make the MWNT layers behave like parallel thermal circuits.     

Molecular Dynamics Study of BaB2O4 in Crystalline and Molten States

Abstract

Structural aspects of BaB₂O₄ liquids have been investigated by the molecular dynamics simulation including the determination on the parameters of the interatomic potential applicable to BaB₂O₄ in both crystalline and molten states. The structure and physical properties of BaB₂O₄ crystals were successfully reproduced by the MD simulation for both α and β phases. The simulated interference function in the liquid state was also in good agreement with the experimental one. Several interesting features on the relaxation phenomena just after melting were reproduced by the simulation that the structure factors of simulated liquid maintain the characteristic features of the original crystal structure for more than 40ps after melting, and the variation of the number of rings formed by B-O bondings was found to increase after melting.     

Study of the stability of aluminium trimeric clusters in aqueous solutions

Ab initio molecular dynamics (AIMD) based on density functional theory has been used to study small aluminium–oxygen complexes in water. Such Al–O clusters have been seen in several recent mass spectrometry studies. In this study, we have focused on trimeric Al–O clusters. The initial very compact trimeric Al–O structures opened up and formed linear Al–O chains. The typical Al–O coordination number in these chain structures was 5. We have performed long (up to 200 ps) AIMD simulations and these chain structures are stable on the nanosecond time scale. We have also studied the reactivity of the Al–O dimer and solvated Al. We found a formation path for the trimeric cluster, which has a action barrier (0.04 eV) and a reaction free energy of ? 0.55 eV. This suggests that the association of a dimer and a monomer Al–O species is fast and thermodynamically a very favourable process.     

Mixing of nanofluids: molecular dynamics simulations and modelling

In this paper we study a mixing scheme, which has recently been proposed for microfluids, on the nanoscale. We do this by performing a series of nonequilibrium molecular dynamics simulations. On the nanoscale the chaotic mixing regime is captured. We discover a new phenomenon where the two mixing fluids exchange positions after leaving the mixing intersection. The results from the molecular dynamics simulations also reveal complex spatio-temporal stream velocity profiles generated by the mixing device. We find that these profiles can be modelled through an approximate analytical solution to the Navier–Stokes equation.     

Water transport through carbon nanotubes with defects

Non-equilibrium molecular dynamics simulations are performed to investigate how changing the number of structural defects in the wall of a (7,7) single-walled carbon nanotube (CNT) affects water transport and internal fluid dynamics. Structural defects are modelled as vacancy sites (missing carbon atoms). We find that, while fluid flow rates exceed continuum expectations, increasing numbers of defects lead to significant reductions in fluid velocity and mass flow rate. The inclusion of such defects causes a reduction in the water density inside the nanotubes and disrupts the nearly frictionless water transport commonly attributed to CNTs.     

Ab initio fragment molecular orbital calculations on the specific interactions between human,mouse and rat PPARα and GW409544

To elucidate the specific interactions between the peroxisome proliferator-activated receptor (PPARα) and ligand GW409544 (GW), we obtained the solvated structures of the PPARα+GW complexes for human, mouse and rat by classical molecular mechanics calculations, and investigated their electronic properties by ab initio fragment molecular orbital calculations. The results indicate that the positively charged amino acids (Lys and Arg) of PPARα make a major contribution to the binding between PPARα and GW. In addition, it was clarified that Ser280 and Tyr314 of human and rat PPARα have a large attractive interaction with GW, while Ser280, Tyr314 and His440 of mouse PPARα have large interaction. These results on the difference in specific interactions between human and mouse/rat PPARα will be useful for predicting the effects of new chemicals on the human body based on the biomedical studies for the experimental animals such as mouse and rat.     

Molecular dynamics simulation of diffusion of hydrogen in binary hydrogen–tetrahydrofuran hydrate

The binary structure II hydrogen–tetrahydrofuran (THF) hydrate was studied with molecular dynamics simulation. The simulations were carried out at 300, 310 K and 10.1 MPa, and with various contents of hydrogen and THF. The migrations of hydrogen molecules from cage to cage were observed. The migration process of hydrogen was also analysed, and the diffusion coefficients of hydrogen in the hydrate were calculated. The calculated diffusion coefficients qualitatively agreed with the experimental data. Double and quintet occupancies of hydrogen molecules were observed in the small and large cages, respectively, without changing the hydrate structure.     

Molecular dynamics simulation on the decomposition of type SII hydrogen hydrate and the performance of tetrahydrofuran as a stabiliser

Molecular dynamics simulation is used to study the decomposition and stability of SII hydrogen and hydrogen/tetrahydrofuran (THF) hydrates at 150 K, 220 K and 100 bar. The modelling of the microscopic decomposition process of hydrogen hydrate indicates that the decomposition of hydrogen hydrate is led by the diffusive behaviour of H₂ molecules. The hydrogen/THF hydrate presents higher stability, by comparing the distributions of the tetrahedral angle of H₂O molecules, radial distribution functions of H₂O molecules and mean square displacements or diffusion coefficients of H₂O and H₂ molecules in hydrogen hydrate with those in hydrogen/THF hydrate. It is also found that the resistance of the diffusion behaviour of H₂O and H₂ molecules can be enhanced by encaging THF molecules in the (5¹²6⁴) cavities. Additionally, the motion of THF molecules is restricted due to its high interaction energy barrier. Accordingly, THF, as a stabiliser, is helpful in increasing the stability of hydrogen hydrate.     

Quantum chemical study on influence of substituents and solvents in reaction complexing ethylene with nickel dithiolene

The influences induced by various terminal substituents and solvents on the reaction mechanism and chemical dynamics of complexing ethylene with Ni dithiolene are theoretically studied by using B3LYP method and Onsager model. It is shown that the reaction should be a two-step process, and the first step is the rate-determining step. We find that the rate constant of the rate-determining step becomes small when the electron-donating ability of the substituents is increased, while it becomes large when the electron-withdrawing ability of the substituents is increased. Subsequently, we consider the solvent effects on the reaction adding ethylene to the simplified hydrogen-substituted nickel dithiolene. It is found that the solvents will make slight changes to the geometries of the reactants, transition states, intermediates and products. However, the corresponding molecular dipole moments become large with the increase of the solvent polarity, which is beneficial to accelerate the reaction. Moreover, we show that, as the solvent polarity becomes large, the activation energies of the reaction decrease exponentially, while the reaction rate constants increase exponentially. These results demonstrate that in polar solvents, the reaction complexing ethylene with Ni dithiolene may become easier and faster to occur, and the product rate is improved. We believe that this research can be seen as a reference for complex and solvent selection in olefin separation process.     

Molecular dynamics simulation of film size and vacancy defect effects on the thermal conductivities of argon thin films

It is well known that there is a size effect for the thermal conductivity of thin films and that vacancy defects in film reduce the film's thermal conduction. In this paper, the film size and vacancy defect effects on the thermal conductivities of argon thin films were studied by molecular dynamics simulations. The results show the existence of phonon boundary scattering. The results also confirm that the theoretical model based on the Boltzmann equation can accurately model the thermal conduction of thin argon films. Both the theoretical and MD results illustrate that, although, both the defect and the thickness of the thin film deduce the thermal conductivity, their physical mechanisms differ.