Pressure-induced Raman-active radial breathing mode transition in single-wall carbon nanotubes

Using classical constant-pressure molecular dynamics simulations and the force constants model, radial breathing mode (RBM) transition of single-wall carbon nanotubes under hydrostatic pressure is reported. With the pressure increased, the RBM shifts linearly toward higher frequency, and the RBM transition occurs at the same critical pressure as the structural transition. The group theory indicates that the RBMs are all Raman-active; however, due to the effect of the frequency transition and the electronic structure change for tube radial deformation, the Raman intensity of the modes becomes so weak as not to be experimentally detected, which is in agreement with a recent experiment by S. Lebedkin [Phys. Rev. B 73, 094109 (2006)]. Furthermore, the calculated RBM transition pressure is well fitted to the cube of diameter (similar to 1/d(3)).

Interface energy and its influence on interface fracture between metal and ceramic thin films in nanoscale

A theoretical model about the size-dependent interface energy between two thin films with different materials is developed by considering the chemical bonding contribution based on the thermodynamic expressions and the structure strain contribution based on the mechanical characteristics. The interface energy decreases with reducing thickness of thin films, and is determined by such available thermodynamic and mechanical parameters as the melting entropy, the melting enthalpy, the shear modulus of two materials, etc. The predicted interface energies of some metal/MgO and metal/Al2O3 interfaces based on the model are consistent with the results based on the molecular mechanics calculation. Furthermore, the interface fracture properties of Ag/MgO and Ni/Al2O3 based on the atomistic simulation are further compared with each other. The fracture strength and the toughness of the interface with the smaller structure interface energy are both found to be lower. The intrinsic relations among the interface energy, the interface strength, and the fracture toughness are discussed by introducing the related interface potential and the interface stress. The microscopic interface fracture toughness is found to equal the structure interface energy in nanoscale, and the microscopic fracture strength is proportional to the fracture toughness. (C) 2010 American Institute of Physics. [doi:10.1063/1.3501090]

Precursor Film in Dynamic Wetting, Electrowetting, and Electro-Elasto-Capillarity

Dynamic wetting and electrowetting are explored using molecular dynamics simulations. The propagation of the precursor film (PF) is fast and obeys the power law with respect to time. Against the former studies, we find the PF is no slip and solidlike. As an important application of the PF, the electro-elasto-capillarity, which is a good candidate for drug delivery at the micro- or nanoscale, is simulated and realized for the first time. Our findings may be one of the answers to the Huh-Scriven paradox and expand our knowledge of dynamic wetting and electrowetting.

Dislocation nucleation governed softening and maximum strength in nano-twinned metals

In conventional metals, there is plenty of space for dislocations-line defects whose motion results in permanent material deformation-to multiply, so that the metal strengths are controlled by dislocation interactions with grain boundaries(1,2) and other obstacles(3,4). For nano-structured materials, in contrast, dislocation multiplication is severely confined by the nanometre-scale geometries so that continued plasticity can be expected to be source-controlled. Nano-grained polycrystalline materials were found to be strong but brittle(5-9), because both nucleation and motion of dislocations are effectively suppressed by the nanoscale crystallites. Here we report a dislocation-nucleation-controlled mechanism in nano-twinned metals(10,11) in which there are plenty of dislocation nucleation sites but dislocation motion is not confined. We show that dislocation nucleation governs the strength of such materials, resulting in their softening below a critical twin thickness. Large-scale molecular dynamics simulations and a kinetic theory of dislocation nucleation in nano-twinned metals show that there exists a transition in deformation mechanism, occurring at a critical twin-boundary spacing for which strength is maximized. At this point, the classical Hall-Petch type of strengthening due to dislocation pile-up and cutting through twin planes switches to a dislocation-nucleation-controlled softening mechanism with twin-boundary migration resulting from nucleation and motion of partial dislocations parallel to the twin planes. Most previous studies(12,13) did not consider a sufficient range of twin thickness and therefore missed this strength-softening regime. The simulations indicate that the critical twin-boundary spacing for the onset of softening in nano-twinned copper and the maximum strength depend on the grain size: the smaller the grain size, the smaller the critical twin-boundary spacing, and the higher the maximum strength of the material.

A method to calculate the surface tension of a cylindrical droplet

The history of Laplace's equations for spherical and cylindrical droplets and the concept of dividing surface in Gibbs' thermodynamic theory of capillary phenomena are briefly reviewed. The existing theories of surface tensions of cylindrical droplets are briefly reviewed too. For cylindrical droplets, a new method to calculate the radius and the surface tension of the surface of tension is given. This method is suitable to be used by molecular dynamics simulations.

Surface transformation and inversion domain boundaries in gallium nitride nanorods

Phase transformation and subdomain structure in [0001]-oriented gallium nitride (GaN) nanorods of different sizes are studied using molecular dynamics simulations. The analysis concerns the structure of GaN nanorods at 300 K without external loading. Calculations show that a transformation from wurtzite to a tetragonal structure occurs along {0110} lateral surfaces, leading to the formation of a six-sided columnar inversion domain boundary (IDB) in the [0001] direction of the nanorods. This structural configuration is similar to the IDB structure observed experimentally in GaN epitaxial layers. The transformation is significantly dependent on the size of the nanorods.

Asymmetric heat conduction through a weak link

We study the heat conduction of two nonlinear lattices joined by a weak harmonic link. When the system reaches a steady state, the heat conduction of the system is decided by the tunneling heat flow through the weak link. We present an analytical analysis by the combination of the self-consistent phonon theory and the heat tunneling transport formalism, and then the tunneling heat flow can be obtained. Moreover, the nonequilibrium molecular dynamics simulations are performed and the simulations results are consistent with the analytical predictions.

Bluish-white emission from radical carbonyl impurities in amorphous Al2O3 prepared via the Pechini-type sol-gel process

Many efforts have been devoted to exploring novel luminescent materials that do not contain expensive or toxic elements, or do not need mercury vapor plasma as the excitation source. In this paper, amorphous Al2O3 powder samples were prepared via the Pechini-type sol-gel process. The resulting samples were characterized by X-ray diffraction (XRD), Fourier transform infrared (FT-IR) spectroscopy, field emission scanning electron microscopy (FESEM), photoluminescence (PL) excitation and emission spectra, kinetic decay, and electron paramagnetic resonance (EPR).

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.

高温高压下H2O-CH4体系PVTx性质的分子动力学模拟研究

Geological fluids are important components in the earth system. To study thephysical chemistry properties and the evolution of fluid system turns out to be one of the most challenging issues in geosciences. Besides the conventional experimental approaches and theoretical or semi-theoretical modeling, molecular level computer simulation(MLCS) emerges as an alternative tool to quantificationally study the physico-chemical properties of fluid under extreme conditions in order to find out the characteristics and interaction of geological fluids in and around earth. Based on our previous study of the intermolecular potential for pure H2O and thestrict evaluation of the competitive potential models for pure CH4 and the ab initio fitting potential surface across H2O-CH4 molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH4 and the H2O-CH4 mixtures. Comparison of 1941 simulations with experimental PVT data for pure CH4 shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the H2O-CH4 mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase,indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region. After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the H2O-CH4 system covering 673 to 2573 K and 0.01 to 10 GPa. Isochores for compositions < 4 mol% CH4 up to 773 K and 600 MPa are also determined in this thesis.

NCMAS系统熔体结构和动力学性质压力效应的分子动力学研究

Molecular dynamics simulations were used to study the pressure dependence of the structure and the dynamic properties of forsterite melt (Mg_2SiO_4), diopside melt (CaMgSi_2O_6), anorthite melt (CaAl_2Si_2O_8), jadite melt (NaAlSi_2O_6) and albite melt (NaAlSi3O8) from 0 GPa to 25 GPa at about 2000 K and the following conclusions have been reached. Firstly, the ratio of NBO to T (NBO and T denote the content of non-bridging oxygen and the total content of Si~(4+) and Al~(3+) respectively) is closely related to the pressure and the composition of the melts. It decreases monotonously in forsterite, diopside and anorthite melts while increases at the initial stage and then decreases in jadite and albite melts with increasing pressure. At a fixed pressure, the shear viscosity of the melts decreases with increasing NBO/T and the variation rate is almost 150 times higher in fully polymerized melts than that in de-polymerized melts in comparison with anorthite melts. Secondly, it is generally accepted that the formation of the Si and A1 will promote the diffusion of the network-forming ions. The hypothesis is frequently employed to explain the emergence of the maximum self-diffusion coefficient of the network-forming ions in fully polymerized melts. However, I detected that the pressure corresponding to the peak of the self-diffusion coefficient of the network-forming ions is lower than that corresponding to the maximum content of Si and A1, and that there exists an approximately linear relationship between the self-diffusion coefficient of the ions and the breaking frequency of the bonds under a given pressure, which is different from the present understanding about the mechanism of self-diffusion. Thirdly, the relationship between the self-diffusion coefficient of Si~(4+), Al~(3+) and O~(2-) and the shear viscosity of the melts evolves from the Stokes-Einstein equation and Sutherland-Einstein equation to the Eyring equation with increasing pressure. And the key to obtain self-diffusion coefficient from shear viscosity under difference pressures is to determine A. in the Eyring equation. For Si~(4+) and O~(2-), this could be done using the linear relationship between A, and NBO% in anorthite melts. However, this method is inapplicable in other kinds of melts.

分子水平上研究地质流体的物理化学性质

Geological fluids exist in every geosphere of the Earth and play important roles in many processes of material transformations, energetic interchanges and geochemical interactions. To study the physicochemical properties and geochemical behaviors of geological fluids turn Girt to be one of the challenging issues in geosciences. Compared with conventional approaches of experiments and semi-theoretical modeling, computer simulation on molecular level shows its advantages on quantitative predictions of the physicochemical properties of geological fluids under extreme conditions and emerges as a promising approach to find the characteristics of geological fluids and their interactions in different geospheres of the Earth interior.This dissertation systematically discusses the physicochemical properties of typical geological fluids with state-of-the-art computer simulation techniques. The main results can be summarized as follows: (1) The experimental phase behaviors of the systems CH4-C2H6 and. CO2 have been successfully reproduced with Monte Carlo simulations. (2) Through comprehensive isothermal-isobaric molecular dynamics simulations, the PVT data of water hia^e been extended beyond experimental range to about 2000 K and 20 GPa and an improved equation of state for water has been established. (3) Based on extensive computer simulations, am optimized molecular potential for carbon dioxide have been proposed, this model is expected to predict different properties of carbon dioxide (volumetric properties, phase equilibria, heat of vaporization, structural and dynamical properties) with improved accuracies. (4) On the basis of the above researches of the end-members, a set of parameters for unlike interactions has been proposed by non-linear fitting to the ab initio potential surface of CO2-H2O and is superior to the common used mixing rule and the results of prior workers vs/Ith remarkable accuracies, then a number of simulations of the mixture have been carried out to generate data under high temperatures and pressures as an important complement to the limited experiments. (5) With molecular dynamics simulations, various structural, dynamical and thermodynamical properties of ionic solvations and associations have been oomprehensively analyzed, these results not only agree well with experimental data and first principle calculation results, but also reveal some new insights into the microscopic ionic solvation and association processes.

Study of adsorption of methylene blue and new methylene blue in liquid-solid interface by slab optical waveguide spectroscopy

A slab optical waveguide (SOWG) has been used for study of adsorption of both methylene blue (MB) and new methylene blue (NMB) in liquid-solid interface. Adsorption characteristics of MB and NMB on both bare SOWG and silanized SOWG by octadecyltrichlorosilane (ODS) were compared. Effect of pH on adsorption on MB and NMB was investigated. Binding rate constant analysis showed that both MB and NMB on bare SOWG demonstrates larger association constants than those on ODS-SOWG. Interactions of NIB and NMB on bare SOWG and ODS-SOWG were analyzed by molecular mechanics calculation method. The binding energy change was in the following order: ENMB-bare > EMB-bare > ENMB-ODS > EMB-ODS. (c) 2004 Elsevier B.V. All rights reserved.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

The Uniaxial Tensile Deformation of Ni Nanowire: Atomic-Scale Computer Simulations

The mechanical deformations of nickel nanowire subjected to uniaxial tensile strain at 300 K are simulated by using molecular dynamics with the quantum corrected Sutten-Chen many-body force field. We have used common neighbor analysis method to investigate the structural evolution of Ni nanowire during the elongation process. For the strain rate of 0.1%/ps, the elastic limit is up to about 11% strain with the yield stress of 8.6 GPa. At the elastic stage, the deformation is carried mainly through the uniform elongation of the distances between the layers (perpendicular to the Z-axis) while the atomic structure remains basically unchanged. With further strain, the slips in the {111} planes start to take place in order to accommodate the applied strain to carry the deformation partially, and subsequently the neck forms. The atomic rearrangements in the neck region result in a zigzag change in the stress-strain curve; the atomic structures beyond the region, however, have no significant changes. With the strain close to the point of the breaking, we observe the formation of a one-atom thick necklace in Ni nanowire. The strain rates have no significant effect on the deformation mechanism, but have some influence on the yield stress, the elastic limit, and the fracture strain of the nanowire.