非线性化学动力学发展的新阶段——浓度场方程及浓度场理论

在综述各种非线性化学动力学研究发展的基础上,重点介绍了最新研究成果浓度场理论的主要内容。该理论根据质量作用定律和广义相对性原理,构建了非线性复杂反应动力学方程即浓度场方程,并给出了扩散、结晶、吸附、传热及相变等9种基本动力学类型的机理指数,解决了热分析动力学积分不收敛、理论基础不完善等问题,合理解释了分形子动力学所谓"记忆效应"和"分数级反应级数"问题,并在通过数学方程和图线全面、直观、定量表达并解释化学振荡、化学分岔、多重定态等各种非线性化学现象的同时,给出了"三级反应的双解性质,是产生各种非线性化学现象的根本原因"的重要结论。浓度场方程及浓度场理论具有很强的理论研究和实际应用价值。

Molecular dynamics simulation of fracture behaviors of <110> tilt grain boundaries in γ-TiAl

Molecular dynamics(MD) simulations were carried out to study the fracture behaviors of several symmetric tilt grain boundaries in γ-Ti Al bicrystals with <110> misorientation axes. Tensile deformation along direction perpendicular to grain boundary was simulated under various strain rates and temperatures. The results indicate that the relative orientation of the grains and the presence of certain atom units are two critical factors of the interface structure affecting the stress required for dislocation nucleation. Dislocations nucleate and extend at or near the symmetric tilt grain boundaries during the tensile deformation of Σ3(111) 109.5°, Σ9(221) 141.1° and Σ27(552) 148.4° interfaces. For Σ27(115) 31.6° and Σ11(113) 50.5° interfaces, the interfaces fractured directly in a cleavage manner due to no dislocation emitted from the boundary. The tensile fracture mechanisms of the bicrystals are that micro-cracks nucleate at the grain boundary and propagate along the interface. The variance of crack propagation is whether there is accommodation of plastic region at the crack tips.

New Hexagonal-rhombic Trilayer Ice Structure Confined between Hydrophobic Plates

We perform molecular dynamics simulations for water confined between two smooth hydrophobic walls and observe two crystalline structures with one being first reported. Both of these structures obey the ice rule. The novel ice phase is a flat hexagonal-rhombic trilayer ice, obtained under 1 GPa load at wall separation of 1.0 nm. In this structure, the water molecules in the two layers next to one of the walls （outer layers） and in the middle layer form hexagonal rings and rhombic rings, respectively. For a molecule in the outer layers, three of its four hydrogen bonds are in the same layer, and the other one hydrogen bond connects to the middle layer. For a molecule in the middle layer, only two of its four hydrogen-bonds are located in the same layer, and the other two connect to two different outer layers. Despite their different motifs, the area densities of the three layers are almost equal. The other structure is a flat hexagonal bilayer ice produced at wall separation of 0.8 nm under lateral pressure of 100 MPa, analogous to a system demonstrated by Koga et al [Phys. Rev. Lett. 79, 5262 （1997）]. Both first-order and continuous phase transitions take place in these simulations.

Nucleation Mechanism of Iron in an External Magnetic Field

In this work,the solidification of liquid iron with or without external magnetic field was investigated by using two molecular dynamics methods,namely direct cooling and two-phase simulation.The influence of external magnetic field on the solidification is characterized by the critical temperature and radial distribution functions.Our computational results show that under external magnetic field,the solidification point tends to decrease significantly.By further analyzing the diffusion coefficients and viscosity,we attribute the effect to the stronger fluctuation of liquid iron atoms driven by the external magnetic field.

Effects of tensile temperatures on phase transformations in zirconium by molecular dynamics simulations

The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental observation under high resolution transmission electron microscopy.The results show that externally applied loading first induced the HCP to body-centered cubic(BCC)phase transition in the Pitsch-Schrader(PS)orientation relationship(OR).Then,the face-centered cubic(FCC)structure transformed from the BCC phase in the Bain path.However,the HCP-to-BCC transition was incomplete at 100 K and 300 K,resulting in a prismatic-type OR between the FCC and original HCP phase.Additionally,at the temperature ranging from 100 K to 600 K,the inverse BCC-to-HCP transition occurred locally following other variants of the PS OR,resulting in a basal-type relation between the newly generated HCP and FCC phases.A higher tensile temperature promoted the amount of FCC phase transforming into the BCC phase when the strain exceeded 45%.Besides,the crystal stretched at lower temperatures exhibits relatively higher strength but by the compromise of plasticity.This study reveals the deformation mechanisms in HCP-Zr at different temperatures,which may provide a better understanding of the deformation mechanism of zirconium alloys under different application environments.

Special nanostructures in Al-Mg alloys subjected to high pressure torsion

Deformation twins and stacking faults were observed in nanostructure A1-Mg alloys subjected to high pressure torsion. These observations are surprising because deformation twinnings have never been observed in their coarse-grained counterparts under normal conditions. Experimental evidences are introduced on non-equilibrium grain boundaries, deformation twinnings and partial dislocation emissions from grain boundaries. Some of these features can be explained by the results reported from molecular-dynamics simulations of pure FCC metals. Special emphasis is laid on the recent observations of high density hexagonal and rhombic shaped nanostructures with an average size of 3 nm in the A1-Mg alloys processed by high pressure torsion. A possible formation process of these nanostructures is proposed based on molecular-dynamics simulations.

Ring Polymer Molecular Dynamics of the C(^(1)D)+H_(2) Reaction on the Most Recent Potential Energy Surfaces

Ring polymer molecular dynamics(RPMD)calculations for the C(^(1)D)+H_(2)reaction are performed on the Zhang-Ma-Bian ab initio potential energy surfaces(PESs)recently constructed by our group,which are unique in very good descriptions of the regions around conical intersections and of van der Waals(vdW)interactions.The calculated reaction thermal rate coefficients are in very good agreement with the latest experimental results.The rate coefficients obtained from the ground˜a^(1)A′ZMB-a PES are much larger than those from the previous RKHS PES,which can be attributed to that the vdW saddles on our PESs have very different dynamical effects from the vdW wells on the previous PESs,indicating that the RPMD approach is able to include dynamical effects of the topological structures caused by vdW interactions.The importance of the excited˜b^(1)A′′ZMB-b PES and quantum effects in the title reaction is also underscored.

Molecular dynamics simulation on plasticity deformation mechanism and failure near void for magnesium alloy

The plastic deformation processes of magnesium alloys near a void at atomic scale level were examined through molecular dynamics(MD)simulation.The modified embedded atom method(MEAM)potentials were employed to characterize the interaction between atoms of the magnesium alloy specimen with only a void.The void growth and crystal failure processes for hexagonal close-packed(hcp)structure were observed.The calculating results reveal that the deformation mechanism near a void in magnesium alloy is a complex process.The passivation around the void,dislocation emission,and coalescence of the void and micro-cavities lead to rapid void growth.

Necessity of Integral Formalism

To describe the physical reality, there are two ways of constructing the dynamical equation of field, differential formalism and integral formalism. The importance of this fact is firstly emphasized by Yang in case of gauge field [Phys. Rev. Lett. 33 （1974） 44fi], where the fact has given rise to a deeper understanding for Aharonov-Bohm phase and magnetic monopole [Phys. Rev. D 12 （1975） 3846]. In this paper we shall point out that such a fact also holds in general wave function of matter, it may give rise to a deeper understanding for Berry phase. Most importantly, we shall prove a point that, for general wave function of matter, in the adiabatic limit, there is an intrinsic difference between its integral formalism and differential formalism. It is neglect of this difference that leads to an inconsistency of quantum adiabatic theorem pointed out by Marzlin and Sanders [Phys. Rev. Lett. 93 （2004） 160408]. It has been widely accepted that there is no physical difference of using differential operator or integral operator to construct the dynamical equation of field. Nevertheless, our study shows that the Schroedinger differential equation （i.e., differential formalism for wave function） shall lead to vanishing Berry phase and that the Schroedinger integral equation （i.e., integral formalism for wave function）, in the adiabatic limit, can satisfactorily give the Berry phase. Therefore, we reach a conclusion： There are two ways of describing physical reality, differential formalism and integral formalism; but the integral formalism is a unique way of complete description.

A molecular dynamics investigation of the deformation mechanism and shape memory effect of epoxy shape memory polymers

Following deformation, thermally induced shape memory polymers(SMPs) have the ability to recover their original shape with a change in temperature. In this work, the thermomechanical properties and shape memory behaviors of three types of epoxy SMPs with varying curing agent contents were investigated using a molecular dynamics(MD) method. The mechanical properties under uniaxial tension at different temperatures were obtained, and the simulation results compared reasonably with experimental data. In addition, in a thermomechanical cycle, ideal shape memory effects for the three types of SMPs were revealed through the shape frozen and shape recovery responses at low and high temperatures, respectively, indicating that the recovery time is strongly influenced by the ratio of E-51 to 4,4'-Methylenedianiline.

Size-dependent transition of the deformation behavior of Au nanowires

Inspired by the controversy over tensile deformation modes of single-crystalline 〈110〉/{111} Au nanowires, we investigated the dependency of the deformation mode on diameters of nanowires using the molecular dynamics technique. A new criterion for assessing the preferred deformation mode-slip or twin propagation--of nanowires as a function of nanowire diameter is presented. The results demonstrate the size-dependent transition, from superplastic deformation mediated by twin propagation to the rupture by localized slips in deformed region as the nanowire diameter decreases. Moreover, the criterion was successfully applied to explain the superplastic deformation of Cu nanowires.

Investigation on structure and dynamic property of liquid Pd-Cu-Ni-P alloys using ab initio molecular dynamics simulation

Based on the short-range order,it is found that the abundance of the P-centered P-transition-metal clusters are the common feature among the liquid Pd-Cu-Ni-P alloys,and hence this feature alone could not uncover the underlying mechanisms of the variation of glass forming ability among the liquid alloys.For the so called similar elements such as Cu and Ni,their behaviors are significantly different when interacting with Pd or P atoms.Cu has weak bonding with both Pd and P while Ni has very strong bonding with P but nearly no bonding with Pd.The different bonding characters thus underlie the phenomenon that in the best glass formers the ratio of the two similar elements often deviates from 1:1.Only if the parameters of chemical short-range order of Cu and Ni around P become closest to each other the best glass forming ability is reached.It is also illustrated that the calculated dynamic properties are very helpful to locate the composition of the best glass former.