Solvation properties of Li+ and CI− in water: molecular dynamics simulation with a non-rigid model

Molecular dynamics (MD) simulations were carried out to investigate the solvation properties of Li⁺ and C1^? ions in water with a relatively accurate but rarely used non-rigid model, RWK2, in this study. A new set of ion-water interaction parameters was evaluated from the experimental data and first principles calculation results of stable clusters, Li⁺(H₂O) _n (n = 1- 6) and CI^?(H₂O) _n (n = 1–4). With the ion-water potential parameters evaluated from the data of the clusters and the water-water potential predetermined from the non-rigid RWK2 model, the structural (radial distribution functions, angular distribution functions, spatial distribution functions, coordination number), dynamical (residence time) and energetic properties of the ionic salvations in bulk water were studied through a comprehensive analysis of our MD simulation outputs. These results not only agree well with experimental data and first principles calculations, but also reveal some new insights into the microscopic ionic salvation processes.     

Distribution and thermophysical properties of n-heptacosane molecules confined in carbon nanotube

Molecular dynamics (MD) simulations were performed to provide an insight about the molecule distribution and thermophysical properties of n-heptacosane confined in the (25, 25) single-walled carbon nanotube (CNT). The results show that an orderly distribution of n-heptacosane molecules along the CNT inner wall is clearly observed. Meanwhile, n-heptacosane confined in CNT exhibits an increased self-diffusion coefficient, a decreased melting point and an enhanced thermal conductivity, compared to the bulk. The simulations reveal that MD is an effective and convenient method to understand the variation characteristics of alkane-based phase change materials confined in CNT on molecular and atomic scale.     

Measuring the implantation depth of silver clusters in graphite

Scanning tunnelling microscopy (STM) and molecular dynamics (MD) simulations have been used to investigate the implantation of Ag₇ ^- clusters into the graphite surface. An experimental measure of the implantation depth of individual clusters is gained via thermal oxidation of the bombarded graphite surfaces. This process results in etching of the cluster-induced defects to form etch pits which grow laterally whilst retaining the depth of the implanted cluster. STM imaging of the etch pits reveals the distribution of implantation depths for deposition energies of 2 keV and 5 keV. Molecular dynamics simulations for clusters of 5 keV energy show that the implantation depth for Ag₇ ^- is largely independent of the impact site on the graphite surface and the cluster orientation. The implantation depth found by MD lies at the upper edge of the experimental depth distribution. Received 30 November 2000     

The melting mechanism in binary Pd0.25Ni0.75 nanoparticles: molecular dynamics simulations

The melting mechanism for Pd_0.25Ni_0.75 alloy nanoparticles (NPs) was investigated using molecular dynamics (MD) simulations with quantum Sutton-Chen many-body potentials. NPs of six different sizes ranging from 682 to 22,242 atoms were studied to observe the effect of size on the melting point. The melting temperatures of the NPs were estimated by following the changes in both the thermodynamic and structural quantities such as the total energy, heat capacity and Lindemann index. We also used a thermodynamics model to better estimate the melting point and to check the accuracy of MD simulations. We observed that the melting points of the NPs decreased as their sizes decreased. Although the MD simulations for the bulk system yielded higher melting temperatures because of the lack of a seed for the liquid phase, the melting temperatures determined for both the bulk material and the NPs are in good agreement with those predicted from the thermodynamics model. The melting mechanism proceeds in two steps: firstly, a liquid-like shell is formed in the outer regions of the NP with increasing temperature. The thickness of the liquid-like shell increases with increasing temperature until the shell reaches a critical thickness. Then, the entire Pd–Ni NP including core-related solid-like regions melts at once.     

Entropy-driven aggregation of adhesion sites of supported membranes

A large number of interesting phenomena related to the insertion of colloidal particles in liquid crystals (LC) have recently been reported. Here, we investigate effects caused by the addition of spherically shaped ferroelectric nanoparticles to a nematic liquid crystal. Using molecular dynamics (MD) simulations, the density of LC molecules, the orientational order parameter, and the polar and azimuthal angle profiles are calculated as functions of the distance to the center of the immersed nanoparticle for different temperatures of the system. We observe that the assembly of ferroelectric nanoparticles enhances the nematic order in the LC medium changing many properties of its host above the nematic-isotropic transition temperature T ^* _NI .     

Analysis of the heterogeneous dynamics of imidazolium-based [Tf2N−] ionic liquids using molecular simulation

The complex dynamic behaviour of the imidazolium-based ionic liquids [C_nmim⁺][Tf₂N^?], n = 4, 8, 12 is examined at various temperatures and at atmospheric pressure using molecular dynamics simulation. An existing all-atom force field is further optimised in order to attain reasonable agreement with experimental data for transport properties, such as self-diffusivities and viscosities. Dynamical heterogeneity phenomena are quantified through the calculation of the non-Gaussian parameter and the deviation of the self-part of the van Hove correlation function from the expected normal distribution. From this analysis, ions that move faster or slower than expected are detected in the system. These subsets of ‘fast’ and ‘slow’ ions form individual clusters consisting of either mobile or immobile ions. Detailed analysis of the ions’ diffusion reveals preferential motion along the direction of the alkyl tail for the cation and along the vector that connects the two sulphur atoms for the anion. For the longest alkyl tails, the heterogeneity in the dynamics becomes more pronounced and is preserved for several nanoseconds, especially at low temperatures.     

DFT study of electronic structure and optical properties of some Ru- and Rh-based complexes for dye-sensitized solar cells

The paper reports Time Dependent Density Functional Theory (TD DFT) calculations providing the structure, electronic properties and spectra of [Ru(II)(bpy)_{3? n} (dcbpy) _n ]²⁺ and [Rh(III)(bpy)_{3? n} (dcbpy) _n ]³⁺ complexes, where bpy?=?2,2′-bipyridyl, dcbpy?=?4,4′-dicarboxy-2,2′-bipyridyl, and n?=?0,?1,?2,?3, studied as possible pigments for dye-sensitized solar cells. The role of the metallic ion and of the COOH groups on the optical properties of these complexes are compared and contrasted and their relevance as dyes for hybrid organic–inorganic photovoltaic cells is discussed. It was found that the optical spectra are strongly influenced by the metallic ion, with visible absorption bands for the Ru(II) complexes and only ultraviolet bands for the Rh(III) complexes. Upon excitation, the extra positive charge of the Rh³⁺ centre tends to draw electrons towards the metal ion, facilitating some charge transfer from the ligand to the metal, whereas in the case of the Ru²⁺ ion the electron transfer is clearly from the metal to the ligand. The carboxyl groups play an important role in strengthening the absorption bands in solution in the visible region. Of the complexes studied, the most suited as pigments for dye-sensitized solar cells are the [Ru(II)(bpy)_{3? n} (dcbpy) _n ]²⁺ complexes with n?=?1 and 2. This is based on the following arguments: (i) their intense absorption band in the visible region, (ii) the presence of the anchoring groups allowing the bonding to the TiO₂ substrate and the charge transfer, and (iii) the good energy level alignment with the conduction band edge of the semiconducting substrate and the redox level of the electrolyte.     

The structure,dynamics and solvation mechanisms of ions in water from long time molecular dynamics simulations: a case study of CaCl2 (aq) aqueous solutions

Systematic long time (5–20 ns) molecular dynamics (MD) simulations have been carried out to study the structural and dynamical properties of CaCl₂ aqueous solutions over a wide range of concentrations (≤9.26 m) in this study. Our simulations reveal totally different structural characteristics of those yielded from short time (≤1 ns) MD simulations [A.A. Chialvo and J.M. Simonson, J. Chem. Phys. 119, 8052 (2003); T. Megyes, I. Bako, S. Balint, T. Grosz, and T. Radnai, J. Mol. Liq. 129, 63 (2006)]. An apparent discontinuity was found at 4–5 m of CaCl₂ in various properties including ion–water coordination number and self-diffusion coefficient of ions, which were first noticed by Phutela and Pitzer in their thermodynamic modelling [R.C. Phutela and K.S. Pitzer, J. Sol. Chem. 12, 201 (1983)]. In this study, residence time was first taken into consideration in the study of Ca²⁺–Cl^? ion pairing, and it was found that contact ion pair and solvent-sharing ion pair start to form at the CaCl₂(aq) concentrations of about 4.5 and 4 m, respectively, which may be responsible for the apparent discontinuity. In addition, the residence time of water molecules around Ca²⁺ or Cl^? showed that the hydration structures of Ca²⁺ and Cl^? are flexible with short residence time (<1 ns). It needs to be pointed out that it takes much longer simulation time for the CaCl₂–H₂O system to reach equilibrium than what was assumed in previous studies.     

Test of a method for sampling the internal degrees of freedom of a flexible solute molecule based on adiabatic decoupling and temperature or force scaling

Simulation of the folding equilibrium of a polypeptide in solution is a computational challenge. Standard molecular dynamics (MD) simulations of such systems cover hundreds of nanoseconds, which is barely sufficient to obtain converged ensemble averages for properties that depend both on folded and unfolded peptide conformations. If one is not interested in dynamical properties of the solute, techniques to enhance the conformational sampling can be used to obtain the equilibrium properties more efficiently. Here the effect on particular equilibrium properties at 298?K of adiabatically decoupling the motion a β-hepta-peptide from the motion of the solvent and subsequently up-scaling its temperature or down-scaling the forces acting on it is investigated. The ensemble averages and rate of convergence are compared to those for standard MD simulations at two different temperatures and a simulation in which the temperature of the solute is increased to 340?K while keeping the solvent at 298?K. Adiabatic decoupling with a solute mass scaling factor s _m ?=?100 and a temperature scaling factor of s _T ?=?1.1 seems to slightly increase the convergence of several properties such as enthalpy of folding, NMR NOE atom–atom distances and ³J-couplings compared to a standard MD simulation at 298?K. Convergence is still slower than that observed at 340?K. The system with a temperature of 340?K for the solute and 298?K for the solvent without scaling of the mass converges fastest. Using a force scaling factor s _V ?=?0.909 perturbs the system too much and leads to a destabilization of the folded structure. The sampling efficiency and possible distortive effects on the configurational distribution of the solute degrees of freedom due to adiabatic decoupling and temperature or force scaling are also analysed for a simpler model, a dichloroethane molecule in water. It appears that an up-scaling of the mass of the solute reduces the sampling more than the subsequent up-scaling of the temperature or down-scaling of the force enhances it. This means that adiabatic decoupling the solute degrees of freedom from the solvent ones followed by an up-scaling of temperature of down-scaling of the forces does not lead to significantly enhanced sampling of the folding equilibrium.     

Melting and Grüneisen parameters of NaCl at high pressure

The Buckingham potential has been employed to simulate the melting and thermodynamic parameters of sodium chloride (NaCl) using the molecular dynamics (MD) method. The constant-volume heat capacity and Grüneisen parameters have been obtained in a wide range of temperatures. The calculated thermodynamic parameters are found to be in good agreement with the available experimental data. The NaCl melting simulations appear to validate the interpretation of superheating of the solid in the one-phase MD simulations. The melting curve of NaCl is compared with the experiments and other calculations at pressure 0－30GPa range.     

Molecular dynamics study of confined ionic liquids in Au nanopore

Molecular dynamics simulations were carried to investigate the structure and dynamics of [BMIM][PF₆] ionic liquid (IL) confined inside a slit-like Au metal nanopore with a pore size of 5.0 nm. The calculations show that the mass and number densities of the confined ILs are oscillatory; the solid-like high density layers are formed in the vicinity of the metal surface. The orientational investigation shows that the imidazolium ring of [BMIM] cations prefers to form a small tilt angle with the pore walls. Furthermore, the mean squared displacement (MSD) calculation indicates that the dynamics of confined ILs are remarkably slower than those observed in bulk systems. Our results suggest that the confinement of the Au nanopore can strongly affect the structural and dynamical properties of the confined ILs.     

Prediction of Structure and Energy of Trans-1,4-Polybutadiene Glassy Surface by Atomistic Simulations of Free-Standing Ultrathin Films

Atomistic modeling of amorphous trans-1,4-polybutadiene (TPBD), using molecular mechanics and molecular dynamics (MD) simulations, is performed to generate three-dimensionally periodic bulk and two-dimensionally periodic thin film condensed phases. The condensed structures are constructed using multiple polymer chains. Structural and energetic relaxations and sampling of properties are performed using MD in the canonical ensemble (NVT) by a procedure that relieves local high-energy spots and brings the system to realistic thermodynamic states. The calculated surface energy for TPBD, 30.72 erg/cm², is in excellent agreement with the reported experimental value of 31 erg/cm². The structure of the surface layers is probed in terms of the atomic mass density variations, bond-bond orientation function profiles, and the distribution of the dihedral angles about the rotatable backbone bonds. The thickness of the surface layer over which the density varies smoothly but rapidly is found to be approximately 15 Å. The level of agreement of the calculated surface energy with the experimental value is superior in comparison to previous investigations in the literature using the atomistic approach for flexible polymers.     

Structure, dynamic and energetic of mixed transition metal clusters

Classical molecular dynamics simulation (MD) with Sutton-Chen potential has been used to generate the minimum energy and to study the thermodynamic and dynamic properties of mixed transition metal cluster motifs of Ag _n Ni_(13?n) for n ?? 13. Literature results of thirteen particle clusters of neat silver and nickel atoms were first reproduced before the successive replacement of the silver atom by nickel. Calculation was repeated for both silver-centred and nickel-centred clusters. It was found that the nickel-centred clusters were more stable than the silver-centred clusters. Heat capacities and hence the melting points of silver and nickel-centred clusters were determined by using the Histogram method. Species-centric order parameters developed by Hewage and Amar were used to understand the dynamic behaviour in the transition of silver-centred clusters to more stable nickel-centred clusters. This species-centric order parameter calculation further confirmed the stability of nickel-centred clusters over those of silver-centred species.     

Solid-liquid phase transition in small water clusters: a molecular dynamics simulation study

Water clusters, (H₂O) _n , of varying sizes (n = 8, 12, 16, 20, 24, 28, 32, 36, and 40) have been studied at different temperatures from 0 to 200 K using molecular dynamics simulations. Transitions between solid and liquid phases were investigated to estimate the melting temperature of the clusters. Although the melting temperatures showed non-monotonic behaviour as a function of cluster size, their general tendency follows the classical relationship T _m ∝ n ^?1/3 to the cluster size n. Moreover, it was observed that the liquid-solid surface tension decreased with the cluster size in a similar way to the liquid-vapour surface tension in bulk water. Upon cooling, ice-like crystals were formed from the smaller clusters with n up to 20, while the larger clusters were transformed to glassy structures. The decrease in the glass transition temperature with the cluster size was observed to be much less than the corresponding melting temperature. The mutual order of the melting and glass-transition temperatures were found to be reversed compared with that observed for bulk water.     

Semiclassical laser theory in the stochastic and thermodynamic frameworks

The thermodynamic properties of the laser distribution in the steadily oscillating state are investigated to determine the minimum characteristic of the entropy production. First, the laser Langevin equation for five random variables is treated in the light of the stochastic calculus to deduce the photon-number rate equationn = – C+(n – n_c) + [A/(1 + sn)](n–n_A), where n_n and n₄ are the two constants of the fluctuation attributed to the noise forces subject to the usual fluctuation-dissipation theorem, withn ₄ < 0 for the inverted atomic population. We then combine the dynamics of the lasing mode with a model open system of the Lebowitz type with two reservoirs for which the entropy production(p) is expressed and made subject to a variational principle: The modified variation scheme, the same as Prigogine's local potential method, is shown to give the exact lasing distributionp as the optimum between two distributions of thermal type with temperatures far from each other.     

BaTiO3晶体结构及弹性的分子动力学模拟

有效的势函数是分子动力学模拟的关键. 引入了一种势函数,该势函数的特点是运用参数r_eff计算原子间的静电作用. 通过分子动力学方法模拟得到了BaTiO₃晶体立方相、四方相结构的对关联函数和X射线衍射谱,计算得出了它们的晶格常数及弹性常数. 模拟结果与实验结果符合较好.该势函数可以有效地模拟BaTiO₃晶体的热学和力学性能. 关键词： 分子动力学模拟 势函数 3铁电晶体')" href="#">BaTiO₃铁电晶体     

The statistics of Curie-Weiss models

LetS _n denote the random total magnetization of ann-site Curie-Weiss model, a collection ofn (spin) random variables with an equal interaction of strength 1/n between each pair of spins. The asymptotic behavior for largen of the probability distribution ofS _n is analyzed and related to the well-known (mean-field) thermodynamic properties of these models. One particular result is that at a type-k critical point (S _n-nm)/n^1–1/2k has a limiting distribution with density proportional to exp[-_s ^2k/(2k)!], wherem is the mean magnetization per site and A is a positive critical parameter with a universal upper bound. Another result describes the asymptotic behavior relevant to metastability.Research supported in part by National Science Foundation Grants MPS 76-06644 (to RSE) and MPS 74-04870 A01 (to CMN).     

Influence of Fe doping on electrical properties of LCMO

The polycrystalline samples La_0.67Ca_0.33Mn_(1?x)Fe _x O₃ (x?=?0.00,?0.01,?0.03, and 0.1) have been grown in single phase by solid state route. The analysis of the reaction has been done by thermogravimetry and differential thermal analysis measurements. DC electrical resistivity measurements have been carried out down to 15?K. The samples with x?=?0.00, 0.01, and 0.03 exhibit metal–insulator (MI) transition at temperatures 221.5?K, 217?K, and 215?K respectively, whereas the sample with x?=?0.1 is insulating in nature for entire temperature range. Interestingly, the electric transport properties of these samples are not consistent with their magnetic phase transitions and the samples show MI transition at a temperature, T _MI, which is significantly lower than the paramagnetic to ferromagnetic transition temperature (T _c). The resistivity data below T _MI has been analyzed using the empirical relation ρ?=?ρ₀?+?ρ₁ T ⁿ and the data above this temperature has been analyzed using two existing models, Mott's variable range hopping model and spin polaronic conduction model.     

Determination of the Solubility Parameter of Epoxidized Soybean Oil by Inverse Gas Chromatography

The thermodynamic parameters of epoxidized soybean oil (ESO) were determined by means of inverse gas chromatography (IGC) in the temperature range of 303.15 K-343.15 K. Two groups of probe solvents with different chemical natures and polarities were used to obtain information about ESO. The thermodynamic parameters—such as molar heat of sorption, weight fraction activity coefficient, Flory–Huggins interaction parameter, and partial molar heat of mixing—were obtained to judge the interactions between ESO and solvents at the studied temperatures. Also, the solubility parameters of ESO were found by plotting the graph of δ₁ ²/(RT) – χ^∞ ₁₂/V₁ vs. solubility parameters, δ₁, of the probes. The results showed that the selected solvent ethyl acetate was a moderate solvent for ESO; n-hexane was a moderate (but close to a better) solvent; while the probes n-heptane, n-octane, n-nonane, and chloroform were excellent solvents. From the extrapolation to 298 K, the solubility parameter value of ESO was 16.70 (J/cm³)^0.5.     

Theoretical investigations of rotational phenomena and dielectric properties in a nematic liquid crystal

The molecular dynamics (MD) simulation, based on a realistic atom-atom interaction potential, was performed on 4-n-pentyl-4'-cyanobiphenyl (5CB) in the nematic phase. The rotational viscosity coefficients (RVCs) γ _i, (i = 1, 2) and the ratio of the RVCs λ = - γ ₂/γ ₁ were investigated. Furthermore, static and frequency-dependent dielectric constants and ε were calculated using parameters obtained from the MD simulation. Time correlation functions were computed and used to determine the rotational diffusion coefficient, D _⊥. The RVCs and λ were evaluated using the existing statistical-mechanical approach (SMA), based on a rotational diffusion model. The SMA rests on a model in which it is assumed that the reorientation of an individual molecule is a stochastic Brownian motion in a certain potential of mean torque. According to the SMA, γ _i are dependent on the orientational order and rotational diffusion coefficients. The former was characterized using: i) orientational distribution function (ODF), and ii) a set of order parameters, both derived from analyses of the MD trajectory. A reasonable agreement between the calculated and experimental values of γ _i and λ was obtained. Received 22 March 2000 and Received in final form 8 October 2000