Molecular simulation study of dynamical properties of room temperature ionic liquids with carbon pieces

Room temperature ionic liquids (RTILs) with dispersed carbon pieces exhibit distinctive physiochemical properties. To explore the molecular mechanism, RTILs/carbon pieces mixture was investigated by molecular dynamics (MD) simulation in this work. Rigid and flexible carbon pieces in the form of graphene with different thicknesses and carbon nanotubes in different sizes were dispersed in a representative RTIL 1-butyl-3-methyl-imidazolium dicyanamide ([Bmim][DCA]). This study demonstrated that the diffusion coefficients of RTILs in the presence of flexible carbons are similar to those of bulk RTILs at varying temperatures, which is in contrast to the decreased diffusion of RTILs in the presence of rigid carbons. In addition, interfacial ion number density at rigid carbon surfaces was higher than that at flexible ones, which is correlated with the accessible external surface area of carbon pieces. The life time of cation-anion pair in the presence of carbon pieces also exhibited a dependence on carbon flexibility. RTILs with dispersed rigid carbon pieces showed longer ion pair life time than those with flexible ones, in consistence with the observation in diffusion coefficients. This work highlights the necessity of including the carbon flexibility when performing MD simulation of RTILs in the presence of dispersed carbon pieces in order to obtain the reliable dynamical and interfacial structural properties.     

单壁碳纳米管中甲烷吸附的分子模拟

采用巨正则系统MonteCarlo方法研究了甲烷在单壁碳纳米管(Singlewallcarbonnanotube,SWNT)中于低温74.05K下的吸附等温线及吸附机理,发现在两个较小的孔径(1.225nm和1.632nm)下单壁碳纳米管中甲烷的吸附有着明显的微孔所独有的“填充效应”,而在2.04nm以上的孔的吸附中会出现毛细凝聚现象。通过模拟知道发生毛细凝聚的必要条件是孔内能至少容纳下两层粒子,此外还导出在恒定温度下毛细凝聚吸附量与SWNT孔径关系。本文还模拟了常温300K下甲烷在SWNT内的吸附,对比了2.04nm和4.077nm两种孔径的SWNT吸附甲烷的等温线,推荐在4.077nm孔中的适宜吸附存储压力为5.0～6.0MPa,吸附质量分数可达16%～19%.     

Molecular dynamics simulation study of the structural characteristics of water molecules confined in functionalized carbon nanotubes

Molecular dynamics (MD) simulations were performed to study the structural properties of water molecules confined in functionalized carbon nanotubes (CNTs). Four CNTs, two armchair-type (6, 6), (7, 7) and two zigzag-type (10, 0), (12, 0) CNTs, representing different helicities and different diameters, were chosen and functionalized at their open ends by the hydrophilic -COOH and the hydrophobic -CH3 groups. The structural properties of water molecules inside the functionalized CNTs, including the orientation distributions of dipole moment and O-H bonds, the length of the single-file water chain, and the average number of hydrogen bonds, were analyzed during a process of simulations. MD simulation results in this work showed that the -CH3 functional groups exert little special effects on the structural properties of water molecules. It is mainly due to the relatively small size of the -CH3 group and its hydrophobic nature, which is consistent with hydrophobic CNTs. For CNTs functionalized by -COOH groups, the configurations of -COOH groups, incurvature or excurvature, determine whether water molecules can enter the CNTs. The incurvature or excurvature configurations of -COOH groups are the results of synergy effects of the CNTs' helicity and diameter and control the flow direction of water molecules in CNTs.     

Molecular dynamics simulation of single wall carbon nanotubes polymerization under compression

Single wall carbon nanotubes (SWCNTs) often aggregate into bundles of hundreds of weakly interacting tubes. Their cross-polymerization opens new possibilities for the creation of new super-hard materials. New mechanical and electronic properties are expected from these condensed structures, as well as novel potential applications. Previous theoretical results presented geometric modifications involving changes in the radial section of the compressed tubes as the explanation to the experimental measurements of structural changes during tube compression. We report here results from molecular dynamics simulations of the SWCNTs polymerization for small diameter arm chair tubes under compression. Hydrostatic and piston-type compression of SWCNTs have been simulated for different temperatures and rates of compression. Our results indicate that large diameter tubes (10,10) are unlike to polymerize while small diameter ones (around 5 A) polymerize even at room temperature. Other interesting results are the observation of the appearance of spontaneous scroll-like structures and also the so-called tubulane motifs, which were predicted in the literature more than a decade ago.     

Supported ionic liquid phase catalysis on functionalized carbon nanotubes

Highly active rhodium catalysts have been prepared by immobilization of an ionic liquid film on carbon nanotubes functionalized with imidazolium-based ionic moieties.     

Molecular simulation of guanidinium-based ionic liquids

A new systematic all-atom force field was developed for cyclic guanidinium-based ionic liquids (ILs) based on the AMBER force field. Optimized molecular geometries and equilibrium bond lengths and angles were obtained by ab initio calculations, and charges were allocated to each atom center by fitting the ab initio electrostatic potential. Molecular dynamics simulations were performed for eleven kinds of ILs that are comprised of NO3(-) anions and cyclic guanidinium-based cations. Validation was carried out by comparing our simulated densities with experimental and calculated data from the literature. Transport properties such as self-diffusion coefficients, viscosities, and conductivities were calculated by molecular dynamic simulation, and their dependence on the length of the alkyl chains of cyclic guanidinium-based cations are discussed. Radial distribution functions and spatial distribution functions were investigated to depict the microscopic structures of the ILs, and the relationship between their properties and microstructures is also discussed.     

Nonenzymatic hydrogen peroxide sensor based on a polyaniline-single walled carbon nanotubes composite in a room temperature ionic liquid

A film consisting of polyaniline and single-walled carbon nanotubes was electro-polymerized on a platinum electrode in a room temperature ionic liquid. This resulted in a selective and non-enzymatic electrode for sensing hydrogen peroxide (HP). The morphology of the sensor film was characterized by scanning electron microscopy, cyclic voltammetry, and chronoamperometry. The composite film synthesized in an ionic liquid provided an electrode with enhanced selectivity and excellent stability. The electrode exhibited good electrocatalytic properties in terms of detection of HP, with a linear range from 5 µM to 1 mM, a detection limit of 1.2 µM, and a response time of around 4 s.     

DFT study of ionic peapod structures from single-walled carbon nanotubes and Lindqvist tungstates

DFT investigations have been carried out on encapsulation of Lindqvist-type W(6)O(19)(2-) anion inside hydrogenated (n,n) armchair single-walled carbon nanotubes (h-CNTs) with n = 8, 9, 10 to understand the confinement effect of the CNTs on the rotation of W(6)O(19)(2-). The energy-decomposition analysis (EDA) of interaction between W(6)O(19)(2-) and CNTs shows that with the increase of confinement effect from n = 8, 9, to 10, the destabilizing ΔE(Pauli) plays a more important role in the relative orientation of W(6)O(19)(2-) inside CNTs. For W(6)O(19)(2-)@(9,9) h-CNT, the most stable orientation appears at the y/z angle 45°/36°. The confinement effect reduces significantly the energy gap of W(6)O(19)(2-)@(n,n) h-CNT (n = 8, 9, 10) compared with free W(6)O(19)(2-). Electron transfer from the W(6)O(19)(2-) to CNT is observed.     

碳纳米管内N2吸附行为及等量吸附热的GCMC分子模拟研究

碳纳米管及石墨烯具有高比表面积、高化学稳定性以及高耐蚀性等优点,被认为是一种理想的吸附材料。分子模拟技术的发展和应用丰富了人们对吸附机理研究的方式,而简单气体吸附体系的吸附机理研究对吸附理论的发展有着重要的推动作用。本文以单壁碳纳米管（SWCNT）-N₂吸附体系为研究对象,首先通过透射扫描电镜和氮气吸/脱附测试对所选用碳纳米管的微观孔形貌及吸/脱附等温线进行了表征,然后根据对应孔径参数采用巨正则蒙特卡罗方法对该体系的吸附过程进行了分子模拟,并详细研究了碳纳米管孔径和温度对该体系吸附行为的影响。结果显示,SWCNT孔径越小,吸附能力则越强;孔半径为0.746nm的SWCNT的吸附体系发生凝聚相变的临界温度为66K。通过对等量吸附热进行计算发现,孔半径0.746、1.15、1.56和1.83 nm的SWCNT-N₂吸附体系对应的初始固-液等量吸附热分别为10.9、9.2、8.6和8.4 kJ/mol。67.5K时,孔半径1.56和1.83 nm的吸附体系的等量吸附热有热峰出现。     

Green chemical functionalization of single-walled carbon nanotubes in ionic liquids

Single walled carbon nanotubes (SWNTs) are exfoliated and functionalized predominantly as individuals by grinding them for minutes at room temperature with aryldiazonium salts in the presence of ionic liquids (ILs) and K(2)CO(3). This constitutes an extremely rapid and mild green chemical functionalization process for obtaining the individualized nanotube structures. A number of ILs and various reaction conditions were surveyed. Raman, XPS, UV/vis/NIR spectroscopies, thermogravimetric analysis, and atomic force and transmission electron microscopies were used to characterize the products.     

Transformation of ionic liquid into carbon nanotubes in confined nanospace

We have developed a two-step filling process for the nano-reaction of ionic liquid in a tip-closed SWNT, where fullerenes are inserted at the end of the host SWNT as a plug to prevent the leakage of the confined ionic liquid during heat treatment.     

Molecular dynamics simulation on influence of temperature effect on electro-coalescence behavior of nano-droplets

In electric dehydration of crude oil, the dewatering efficiency can be improved by raising emulsion temperature properly which reduces the viscosity of crude oil. However, it should be noticed that the emulsion temperature does not only affect the emulsion viscosity but also nano-droplets dynamics behavior which impacts the coalescence efficiency either. Therefore the influence of temperature effect on the electro-coalescence of nano-droplets is studied by a molecular dynamics method. The results show that the temperature presents an active or negative effect, depending on the competitive relation between electrostatic interaction and thermal motion. Two stages are distinguished according to the dominant mechanism. During stage I, governed by the electrostatic interaction, lower temperature promotes the polarization and leads to an acceleration of the droplets coalescence, but higher temperature restrains the coalescence process due to molecules thermal motion breaking the polarization process. During stage II, governed by the thermal motion, lower temperature improves the coalescence because of a diffusion effect, but higher temperature deteriorates electro-coalescence because of a violent molecular thermal motion. Additionally, hydrogen bond and radial distribution functions are obtained by statistics to describe droplets micromorphology, which explains the reason why the droplet forms longer chain structure at the critical electric field.     

Molecular dynamics study of the hydration of the hydroxyl radical at body temperature

Classical molecular dynamics (MD) simulation of ˙OH in liquid water at 37 °C has been performed using flexible models of the solute and solvent molecules. We derived the Morse function describing the bond stretching of the radical and the potential for ˙OH-H(2)O interactions, including short-range interactions of hydrogen atoms. Scans of the potential energy surface of the ˙OH-H(2)O complex have been performed using the DFT method with the B3LYP functional and the 6-311G(d,p) basis set. The DFT-derived partial charges, ±0.375e, and the equilibrium bond-length, 0.975 ?, of ˙OH resulted in the dipole moment of 1.76 D. The radical-water radial distribution functions revealed that ˙OH is not built into the solvent structure but it rather occupies distortions or cavities in the hydrogen-bonded network. The solvent structure at 37 °C has been found to be the same as that of pure water. The hydration cage of the radical comprises 13-14 water molecules. The estimated hydration enthalpy -42 ± 5 kJ mol(-1) is comparable with the experimental value -39 ± 6 kJ mol(-1) for 25 °C. Inspection of hydrogen bonds showed the importance of short-range interaction of hydrogen atoms and indicated that neglect of the angular condition greatly overestimates the number of the H-acceptor radical-water bonds. The mean number ?n = 0.85 of radical-water H-bonds has been calculated using geometric definition of H-bond and ?n = 0.62 has been obtained when the energetic condition, E(da)≤-8 kJ mol(-1), was additionally considered. The continuous lifetimes of 0.033 ps and 0.024 ps have been estimated for the radical H-donor and the H-acceptor bonds, respectively. Within statistical uncertainty the radical self-diffusion coefficient, (2.9 ± 0.6) × 10(-9) m(2) s(-1), is the same as (3.1 ± 0.5) × 10(-9) m(2) s(-1) calculated for water in solution and in pure solvent. To the best of our knowledge, this is the first study of the ˙OH(aq) properties at a biologically relevant body temperature.     

Computer simulation of amino acid sorption on carbon nanotubes

A computer simulation of complexes of (6,6) open carbon nanotubes (CNTs) with neutral molecules, zwitterions and glycine, alanine, and phenylalanine amino acid anions is performed. In starting structures amino acids are arranged in three types: on the external side face, the open end, and inside CNT. The structure is optimized within the density functional theory with regard to the GD3 dispersion correction with and without taking into account solvation effects. It is found that the greatest CNT–amino acid interaction occurs in the neutral aqueous medium at dissociative chemisorption of the zwitterion (adsorption energy 80-90 kcal/mol) and in the basic medium at anion chemisorption (energy ~48-50 kcal/mol) on the open CNT end.     

碳纳米管与咪唑类离子液体相互作用机制的理论研究

离子液体(IL)与碳纳米管(CNT)形成的复合体系在许多重要领域有潜在的应用价值,然而人们对于其微观结构的理解尚不清晰,对IL与CNT相互作用机理的认识还存有争议.本文以(8,0)CNT与1-丁基-3-甲基咪唑六氟化磷IL为例研究了CNT与IL的相互作用,分别探讨了1-丁基-3-甲基咪唑阳离子([BMIM+])、六氟化...     

Molecular simulation of the hydration Gibbs energy of barbiturates

In the present work, molecular dynamics calculations of the Gibbs energy of hydration of 10 different substituted barbiturates in SPC/E water were performed using thermodynamic integration. Given that experimental determination of the Gibbs hydration energy for this class of compounds is currently unfeasible, computer simulations appear as the only alternative for the estimation of this important quantity. Several simulation parameters are discussed and optimized based on calculations for barbituric acid. It is concluded that accounting for electrostatic interactions with the Reaction-Field method can be up to two times faster than with Particle-Mesh-Ewald method, without loss of accuracy. Different number of solvent molecules and simulation lengths were also tested. Lennard-Jones and electrostatic contributions were scaled down to zero in an independent way. It is shown that the electrostatic contribution is dominant (representing approximately 90% of the total Gibbs energy of hydration) and that barbiturate intra-molecular interactions cannot be neglected. The importance of the electrostatic contribution is attributed to the formation of hydrogen bonds between the barbiturates and water, which play an important role in the solvation process. The influence of the different substituents and their contribution to the Gibbs energy of hydration was assessed. Finally, the Lennard-Jones contributions and the total hydration Gibbs energy can both be correlated against molecular weight or partition coefficient data for mono- and di-substituted barbiturates.