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.

笼形水簇和溶解甲烷之间平均力势能的约束型分子动力学计算及其在天然气水合物成核研究中的意义

That the dodecahedral water cluster (DWC) can adsorb dissolved methane molecules, an important phenomenon related to the hydrate nucleation study, has been observed through molecular dynamics simulations, but it has not been explained satisfactorily [Guang-Jun Guo; Yi-Gang Zhang; Hua Liu. J. Phys. Chem. C, 2007, 111, 2595]. In order to explain this phenomenon by using the potential of mean force (PMF) between the DWC and the dissolved methane, we perform several series of constrained molecular dynamics simulations in the methane-water system. The distance between the center of DWC and the methane molecule is constrained from 5 Å to 18 Å by adding 0.2 Å every time. For each fixed distance, we perform 20 independent simulations to improve the statistical precision. We first get the constraint force between the DWC and the dissolved methane in each simulation and then calculate the PMF by integrating these forces. Subsequently, the radial distribution function (RDF) is obtained from the PMF through an equation of statistical mechanics. The results show that the RDF has a sharp peak at about 6.2 Å, successfully explaining why the DWC adsorbs dissolved methane molecules. The preferential binding coefficient is a positive value (=2.05±0.5), indicates that the DWC tends to adsorb dissolved methane rather than water molecules in methane aqueous solutions. The curve of PMF for the DWC encaging a methane almost coincides that for the empty DWC, meaning that it is the DWC rather than the encaged methane who could adsorb dissolved methane molecules. By comparing the curves of PMF for different directions of the DWC relative to the dissolved methane, we find that it is the cage face rather than the cage edge or vertex that plays an essential role when the DWC adsorbing dissolved methane. This research sheds light on the driving force for the methane adsorption, and it is helpful in understanding the nucleation process of methane hydrate.

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.

Polyamorphism in aluminate liquids

McMillan, P. F., Wilson, M., Wilding, M. C. (2003). Polyamorphism in aluminate liquids. Journal of Physics: Condensed Matter, 15 (36), 6105-6121 RAE2008

Accommodating curvature in a highly ordered functionalized metal oxide nanofiber: synthesis, characterization and multi-scale modeling of layered nanosheets

A key element in the rational design of hybrid organic-inorganic nanostructures, is control of surfactant packing and adsorption onto the inorganic phase in crystal growth and assembly. In layered single crystal nanofibers and bilayered 2D nanosheets of vanadium oxide, we show how the chemisorption of preferred densities of surfactant molecules can direct formation of ordered, curved layers. The atom-scale features of the structures are described using molecular dynamics simulations that quantify surfactant packing effects and confirm the preference for a density of 5 dodecanethiol molecules per 8 vanadium attachment sites in the synthesised structures. This assembly maintains a remarkably well ordered interlayer spacing, even when curved. The assemblies of interdigitated organic bilayers on V2O5 are shown to be sufficiently flexible to tolerate curvature while maintaining a constant interlayer distance without rupture, delamination or cleavage. The accommodation of curvature and invariant structural integrity points to a beneficial role for oxide-directed organic film packing effects in layered architectures such as stacked nanofibers and hybrid 2D nanosheet systems.

A multi-scale numerical modelling of crack propagation in a 2D metallic plate

A new multi-scale model of brittle fracture growth in an Ag plate with macroscopic dimensions is proposed in which the crack propagation is identified with the stochastic drift-diffusion motion of the crack-tip atom through the material. The model couples molecular dynamics simulations, based on many-body interatomic potentials, with the continuum-based theories of fracture mechanics. The Ito stochastic differential equation is used to advance the tip position on a macroscopic scale before each nano-scale simulation is performed. Well-known crack characteristics, such as the roughening transitions of the crack surfaces, as well as the macroscopic crack trajectories are obtained.

Ice-binding proteins adsorb to their ligand via anchored clathrate waters

The main success of my thesis has been to establish the mechanism by which antifreeze proteins (AFPs) bind irreversibly to ice crystals, and hence prevent their growth. AFPs organize ice-like water on their ice-binding site, which then merges and freezes with the quasi-liquid layer of ice. This was revealed from studying the exceptionally large (ca. 1.5-MDa) Ca 2+-dependent AFP from the Antarctic bacterium Marinomonas primoryensis (MpAFP). The 34-kDa antifreeze- active region of MpAFP was predicted to fold as a novel Ca 2+-binding β-helix. Site-directed mutagenesis confirmed the model and demonstrated that its ice-binding site (IBS) consisted of solvent-exposed Thr and Asx parallel arrays on the Ca 2+-binding turns. The X-ray crystal structure of the antifreeze region was solved to a resolution of 1.7 Å. Two of the four molecules within the unit cell of the crystal had portions of their IBSs freely exposed to solvent. Identical clathrate-like cages of water molecules were present on each IBS. These waters were organized by the hydrophobic effect and anchored to the protein via hydrogen bonds. They matched the spacing of water molecules in an ice lattice, demonstrating that anchored clathrate waters bind AFPs to ice. This mechanism was extended to other AFPs including the globular type III AFP from fishes. Site-directed mutagenesis and a modified ice-etching technique demonstrated this protein uses a compound ice-binding site, comprised of two flat and relatively hydrophobic surfaces, to bind at least two planes of ice. Reinvestigation of several crystal structures of type III AFP identified anchored clathrate waters on the solvent-exposed portion of its compound IBS that matched the spacing of waters on the primary prism plane of ice. Ice nucleation proteins (INPs), which can raise the temperature at which ice forms in solution to just slightly below 0oC, have the opposite effect to AFPs. A novel dimeric β-helical model was proposed for the INP produced by the bacterium Pseudomonas borealis. Molecular dynamics simulations showed that INPs are also capable of ordering water molecules into an ice- like lattice. However, their multimerization brings together sufficient ordered waters to form an ice nucleus and initiate freezing.

Equilibrium polymerization of cyclic carbonate oligomers

A model of the polymerization of ring oligomers of bisphenol A polycarbonate (BPA-PC) is used to investigate the influence of dimensionality (2D or 3D), density and temperature on the size distribution of the polymer chains. The polymerization step is catalyzed by a single active particle, conserves the number and type of the chemical bonds, and occurs without a significant gain in either potential energy or configurational entropy. Monte Carlo and molecular dynamics simulations show that polymerization of cyclic oligomers occurs readily at high density and is driven by the entropy associated with the distribution of interparticle bonds. Polymerization competes at lower densities with long range diffusion, which favors small molecular species, and is prevented if the system is sufficiently dilute. Polymerization occurs in 2D via a weakly first order transition as a function of density and is characterized by low hysteresis and large fluctuations in the size of polymer chains. Polymerization occurs more readily in 3D than in 2D, and is favored by increasing temperature, as expected for an entropy-driven process. (C) 2001 American Institute of Physics.

Concertedness and solvent effects in multiple proton transfer reactions: The formic acid dimer in solution

The issue of multiple proton transfer (PT) reactions in solution is addressed by performing molecular dynamics simulations for a formic acid dimer embedded in a water cluster. The reactant species is treated quantum mechanically, within a density functional approach, while the solvent is represented by a classical model. By constraining different distances within the dimer we analyze the PT process in a variety of situations representative of more complex environments. Free energy profiles are presented, and analyzed in terms of typical solvated configurations extracted from the simulations. A decrease in the PT barrier height upon solvation is rationalized in terms of a transition state which is more polarized than the stable states. The dynamics of the double PT process is studied in a low-barrier case and correlated with solvent polarization fluctuations. Cooperative effects in the motion of the two protons are observed in two different situations: when the solvent polarization does not favor the transfer of one of the two protons and when the motion of the two protons is not synchronized. This body of observations is correlated with local structural and dynamical properties of the solvent in the vicinity of the reactant. (C) 2000 American Institute of Physics. [S0021-9606(00)51121-0].

Simulation of the surface structure of butylmethylimidazolium ionic liquids

Molecular dynamics simulations of the liquid/vacuum surfaces of the room temperature ionic liquids [bmim][PF6], [bmim][BF4] and [bmim][Cl] have been carried out at various temperatures. The surfaces are structured with a top monolayer containing oriented cations and anions. The butyl side chains tend to face the vacuum and the methyl side chains the liquid. However, as the butyl chains are not densely packed, both anions and rings are visible from the vacuum phase. The effects of temperature and the anion on the degree of cation orientation is small, but the potential drop from the vacuum to the interior of the liquid is greater for liquids with smaller anions. We compare the simulation results with a range of experimental observations and suggest that neutron reflection from samples with protiated butyl groups would be a sensitive probe of the structure.