Effects of information asymmetry on cooperation in the prisoners＇ dilemma game

Effects of inforination asymmetry on cooperation in the prisoners＇ dilemma game are investigated. The amplitude A is introduced to describe the degree of information asymmetry. It is found that there exists an optimal value of amplitude Aopt at which the fraction of cooperation reaches its maximal value. The reason lies in that cooperators on the two-dimensional grid form large clusters at Aopt. In addition, the theoretical analysis in terms of the mean- field theory is used to understand this kind of phenomenon. It is confirmed that the information asymmetry plays an important role in the dynamics of the dilemma games of spatial prisoners.     

Cooperative Dynamics in Lattice-Embedded Scale-Free Networks

We investigate cooperative behaviors of lattice-embedded scale-free networking agents in the prisoner＇s dilemma game model by employing two initial strategy distribution mechanisms, which are specific distribution to the most connected sites （hubs） and random distribution. Our study indicates that the game dynamics crucially depends on the underlying spatial network structure with different strategy distribution mechanism. The cooperators＇ specific distribution contributes to an enhanced level of cooperation in the system compared with random one, and cooperation is robust to cooperators＇ specific distribution but fragile to defectors＇ specific distribution. Especially, unlike the specific case, increasing heterogeneity of network does not always favor the emergence of cooperation under random mechanism. Furthermore, we study the geographical effects and find that the graphically constrained network structure tends to improve the evolution of cooperation in random case and in specific one for a large temptation to defect.     

Shell Model for Elastic and Thermodynamic Properties of Gallium Nitride with Hexagonal Wurtzite Structure

Shell model molecular dynamic simulation with interatomic pair potential is utilized to investigate the elastic and thermodynamic properties of gallium nitride with hexagonal wurtzite structure （w-GaN） at high pressure. The calculated elastic constants Cij at zero pressure and 300 K agree well with the experimental data and other calculated values. Meanwhile, the dependences of the relative volume V/Vo, elastic constants Cij, entropy S, enthalpy H, and heat capacities Cv and Up on pressure are successfully obtained. From the elastic constants obtained, we also calculate the shear modulus G, bulk modulus B, Young＇s modulus E, Poisson＇s ratio v, Debye temperature ΘD, and shear anisotropic factor Ashear on pressures.     

Enhancement of water permeation across nanochannels by partial charges mimicked from biological channels

In biological water channel aquaporins （AQPs）, it is believed that the bipolar orientation of the single-file water molecules inside the channel blocks proton permeation but not water transport. In this paper, the water permeation and particularly the water-selective behaviour across a single-walled carbon nanotube （SWNT） with two partial charges adjacent to the wall of the SWNT are studied by molecular dynamics simulations, in which the distance between the two partial charges is varied from 0.14 nm to 0.5 nm and the charges each have a quantity of 0.5 e. The two partial charges are used to mimic the charge distribution of the conserved non-pseudoautosomal （NPA） （asparagine/proline/alanine） regions in AQPs. Compared with across the nanochannel in a system with one ＋1 ε charge, the water permeation across the nanochannel is greatly enhanced in a system with two ＋0.5 e charges when charges are close to the nanotube, i.e. the two partial charges permit more rapid water diffusion and maintain better bipolar order along the water file when the distance between the two charges and the wall of SWNT is smaller than about 0.05 nm. The bipolar orientation of the single-file water molecules is crucial for the exclusion of proton transfer. These findings may serve as guidelines for the future nanodevices by using charges to transport water and have biological implications because membrane water channels share a similar single-file water chain and positive charged region at centre and provide an insight into why two residues are necessitated in the central region of water channel protein.     

Atomic diffusion across NisoTiso-Cu explosive welding interface： Diffusion layer thickness and atomic concentration distribution

Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of NisoTis0-Cu （at.%）. By using a hybrid method which combines molecular dynamics simulation and classical diffusion the- ory, the thickness of the diffusion layer and the atomic concentration distribution across the welding interface are obtained. The results indicate that the concentration distribution curves at different times have a geometric similarity. According to the geometric similarity, the atomic concentration distribution at any time in explosive welding can be calculated. NisoTis0- Cu explosive welding and scanning electron microscope experiments are done to verify the results. The simulation results and the experimental results are in good agreement.     

Coalescence between Cu57 and Cu5s clusters at a room temperature： molecular dynamics simulations

Three coalescence processes of Cu57-Cu57, Cu57-Cu58, and Cu5s Cu58 clusters at 300 K are investigated by employing molecular dynamics simulations. According to the evolutions of mean square displacement and local atom packing, the coalescence process can be separated into three stages including an approaching stage, a coalescing stage, and a coalesced stage. The simulations show that the coalescence processes and the formed products are sensitive to the respective initial structures of, and the relative configuration between, the two coalescing icosahedron-based clusters.     

Assessment of delocalized and localized molecular orbitals through electron momentum spectroscopy

Recently, there was a hot controversy about the concept of localized orbitals, which was triggered by Grushow＇s work titled ＂Is it time to retire the hybrid atomic orbital？＂ [J. Chem. Educ. 88, 860 （2011）]. To clarify the issue, we assess the delocalized and localized molecular orbitals from an experimental view using electron momentum spectroscopy. The delocalized and localized molecular orbitals based on various theoretical models for CH4, NH3, and H20 are compared with the experimental momentum distributions. Our results show that the delocalized molecular orbitals rather than the localized ones can give a direct interpretation of the experimental （e, 2e） results.     

Molecular constants of LiCI（X^1∑＋） and elastic collisions of two ground-state Cl and Li atoms at low and ultralow temperatures

Interaction potentials for LiCI（X^1∑＋） are constructed by the highly accurate valence internally contracted multireference configuration interaction in combination with a number of large correlation-consistent basis sets, which are used to determine the spectroscopic parameters （D0, De, Re, ωe, ωeχe, Be and αe）. The potentials obtained at the basis sets, i.e., aug-cc-pV5Z-JKFI for Cl and cc-pV5Z for Li, are selected to study the elastic collision properties of Li and Cl atoms at the impact energies from 1.0 ×10^-12 to 1.0× 10-4 a.u. The derived total elastic cross sections are very large and almost constant at ultralow temperatures, and their shapes are mainly dominated by the s-partial wave at very low impact energies. Only one shape resonance can be found in the total elastic cross sections over the present collision energy regime, which is rather strong and obviously broadened by the overlap contributions of the abundant resonances coming from various partial waves. Abundant resonances exist for the elastic partial-wave cross sections until l= 22 partial waves. The vibrational manifolds of the LiCI（X^1∑＋） molecule, which are predicted at the present level of theory and the basis sets cc-pV5Z for Li and the aug-cc-pV5Z-JKFI for Cl, should achieve much high accuracy due to the employment of the large correlation-consistent basis sets.     

Thermal transport property of investigated by molecular Ge34 and d-Ge dynamics and the Slack＇s equation

In this study, we evaluate the values of lattice thermal conductivity κL of type Ⅱ Ge clathrate （Ge34） and diamond phase Ge crystal （d-Ce） with the equilibrium molecular dynamics （EMD） method and the Slack＇s equation. The key parameters of the Slack＇s equation are derived from the thermodynamic properties obtained from the lattice dynamics （LD） calculations. The empirical Tersoff＇s potential is used in both EMD and LD simulations. The thermal conductivities of d-Ge calculated by both methods are in accordance with the experimental values. The predictions of the Slack＇s equation are consistent with the EMD results above 250 K for both Ge34 and d-Ge. In a temperature range of 200-1000 K, the κL value of d-Ge is about several times larger than that of Ge34.     

Multireference configuration interaction study on spectroscopic properties of low-lying electronic states of As2 molecule

The potential energy curves (PECs) of four electronic states (X1Σ+g , e3△u , a 3 Σ-u , and d 3Πg ) of an As 2 molecule are investigated employing the complete active space self-consistent field (CASSCF) method followed by the valence internally contracted multireference configuration interaction (MRCI) approach in conjunction with the correlation-consistent aug-cc-pV5Z basis set. The effect on PECs by the relativistic correction is taken into account. The way to consider the relativistic correction is to employ the second-order Douglas-Kroll Hamiltonian approximation. The correction is made at the level of a cc-pV5Z basis set. The PECs of the electronic states involved are extrapolated to the complete basis set limit. With the PECs, the spectroscopic parameters (Te , Re , ωe , ωexe , ωeye , αe , βe , γe , and Be ) of these electronic states are determined and compared in detail with those reported in the literature. Excellent agreement is found between the present results and the experimental data. The first 40 vibrational states are studied for each electronic state when the rotational quantum number J equals zero. In addition, the vibrational levels, inertial rotation and centrifugal distortion constants of d 3Πg electronic state are reported which are in excellent agreement with the available measurements. Comparison with the experimental data shows that the present results are both reliable and accurate.     

Influence of Ni on Cu precipitation in Fe Cu Ni ternary alloy by an atomic study

The early aging Cu precipitations in Fe-3%Cu and Fe-3%Cu-4%Ni ternary alloys are investigated by molecular dynamics （MD） simulations. The results show that the average size of Cu clusters in Fe-3%Cu-4%Ni alloy is larger than that in Fe-3%Cu alloy. The diffusion of Cu is accelerated by Ni according to the mean square displacement （MSD）. Furthermore, the whole formation process of Cu-rich clusters is analyzed in detail, and it is found that the presence of Ni promotes small Cu-rich clusters to be combined into big ones. Ni atoms prefer to stay at the combination positions of small clusters energetically due to a large number of the first nearest neighbor Cu-Ni interactions, which is verified by first-principles calculations based on density functional theory （DFT）.     

Temperature dependence of Cu20 orientations in the oxidation of Cu （111）/ZnO （0001） by oxygen plasma

The role of temperature on the oxidation dynamics of Cu20 on ZnO （0001） was investigated during the oxidation of Cu （111）/ZnO （0001） by using oxygen plasma as the oxidant. A transition from single crystalline Cu20 （111） orientation to micro-zone phase separation with multiple orientations was revealed when the oxidation temperature increased above 300 ~ C. The experimental results clearly show the effect of the oxidation temperature with the assistance of oxygen plasma on changing the morphology of Cu （111） film and enhancing the lateral nucleation and migration abilities of cuprous oxides. A vertical top-down oxidation mode and a lateral migration model were proposed to explain the different nucleation and growth dynamics of the temperature-dependent oxidation process in the oxidation of Cu （lll）/ZnO （0001）.     

Transport of ions through a （6,6） carbon nanotube under electric fields

The transport of water and ions through carbon nanotubes （CNTs） is crucial in nanotechnology and biotechnology. Previous investigation indicated that the ions can hardly pass through （6,6） CNTs due to their hydrated shells. In the present study, utilizing molecular dynamics simulation, it is shown that the energy barrier mainly originating from the hydrated water molecules could be overcome by applying an electric field large enough in the CNT axis direction. Potential of mean force is calculated to show the reduction of energy barrier when the electric field is present for （Na＋, K＋, C1 ） ions. Consequently, ionic flux through （6,6） CNTs can be found once the electric field becomes larger than a threshold value. The variation of the coordination numbers of ions at different locations from the bulk to the center of the CNT is also explored to elaborate this dynamic process. The thresholds of the electric field are different for Na＋, K＋, and CI- due to their characteristics. This consequence might be potentially applied in ion selectivity in the future.     

Isospin effect of the in-medium nucleon-nucleon cross section in excited nuclear reactions

Using relativistic mean field theory, the neutron and the proton density distribution of 56Ni nuclei could be obtained in the ground state and the excited state. Based on the framework of the quantum molecular dynamics model, the 56Ni nuclei have been simulated in ground state and in the neutron or proton excited state. We then used the three different states of 56Ni to collide with the 56Ni in the ground state. To discuss the evolution of the nuclear stopping in different reactions, two kinds of different excited nuclear reactions were studied at different reaction energies and at different impact parameters. Studies have shown that the nuclear stopping of an excited nuclear reaction is sensitive to the isospin-dependent in-medium nucleon-nucleon cross section, compared with the response value of the ground state nuclear reaction. So, it is better for the excited nuclei to extract the isospin dependence of nucleon-nucleon cross section information.     

Sigma Meson Production in Nuclear Reactions

The sigma meson production in p ＋ ^12C and p A- ^40Ca reactions at the incident energy Ep = 1.5 GeV is investigated within the Quantum Molecular Dynamics model. The simulation results indicate a distinctive A dependence of the sigma production, that is, the increase of A is followed by an increase of the production cross sections. We find that the σ meson production in proton-induced reactions is strongly medium-dependent, and the produced σ mesons decaying in a denser medium experience a stronger mass shift towards lower masses. This mass shift is an experimentally accessible observable in the final state pion pairs, which do not suffer from reabsorption by the surrounding nucleons. It is pointed out that the ratio of measured sigma cross sections as a function of the sigma invariant-mass from various reactions is a good probe to explore the existence of the σ meson in a dense nuclear environment.     

Mechanical Properties of Ni-Coated Single Graphene Sheet and Their Embedded Aluminum Matrix Composites

The effects of Ni coating on the mechanical behaviors of single graphene sheet and their embedded Al matrix composites under axial tension are investigated using molecular dynamics （MD） simulation method. The results show that the Young＇s moduli and tensile strength of graphene obviously decrease after Ni coating. The results also show that the mechanical properties of Al matrix can be obviously increased by embedding a single graphene sheet. From the simulation, we also find that the Young＇s modulus and tensile strength of the Ni-coated graphene/Al composite is obviously larger than those of the uncoated graphene/Al composite. The increased magnitude of the Young＇s modulus and tensile strength of graphene/Al composite are 52.27% and 32.32% at 0.01 K, respectively, due to Ni coating. By exploring the effects of temperature on the mechanical properties of single graphene sheet and their embedded Al matrix composites, it is found that the higher temperature leads to the lower critical strain and tensile strength.     

Structural, Spectroscopic and Vibrational Properties for Low-Lying Electronic States of LiCl

A multireference configuration interaction(MRCI) study has been carried out on the LiCl molecule. The potential energy has been calculated over a wide range of internuclear separation for the 21 low-lying electronic states of the LiCl molecule dissociating into Li(2S,2P,3S)+Cl(2P). The(4)1Σ+,(3)1Π, 1-33Σ+, 1-33Π,1,3,1,3Σ,(5)1Σ+,(4)3Σ+,(4)1Π,(4)3Π excited states are studied for the first time in theory. Molecular spectroscopic constants(Re, De,ωe, ωeχe, Be and αe) have been derived for the 9 bound states(X1Σ+,(3)1Σ+,(2)3Σ+,1,3,1,3Σ,(4)1Π,(4)3Π) with a regular shape, and the spectroscopic constants of ground states X1Σ+are in good agreement with available experimental and theoretical values. The relative diferences between experimental values and our values for Re, De, ωe, ωeχe, Be andαe are 1.02%, 0.60%, 1.72%, 9.46%, 2.0%, and 0.75%, respectively. Moreover, vibrational levels of 9 bound states, which have not been investigated experimentally, are computed.     

Surface reconstruction on stishovite Si02, HfO2 and rutile TiO2 （001）

This paper systematically investigates the surface reconstruction processes and patterns on stishovite SiO2, HfO2 and rutile TiO2 （001） by using classical molecular dynamics. It is found that these three surfaces relax instead of reconstruction at 0 K, and have little possibility to reconstruct below 40 K. Above 40 K, surface reconstructions take place as collective atomic motion which can be speeded by higher temperature or compressed strain. Several reconstruction patterns with approximate surface energies are found, and electrostatic potentials on them are also provided in comparison with possible microscopic results.     

Molecular dynamics simulation of an argon cluster filled inside carbon nanotubes

The effects of the diameters of single-walled carbon nanotubes （SWCNTs） （7.83A to 27.40A） and temperature （20 K-45 K） on the equilibrium structure of an argon cluster are systematically studied by molecular dynamics simulation with consideration of the SWCNTs to be fixed. Since the diameters of SWCNTs with different chiralities increase when temperature is fixed at 20 K, the equilibrium structures of the argon cluster transform from monoatomic chains to helical and then to multishell coaxial cylinders. Chirality has almost no noticeable influence on these cylindrosymmetric structures. The effects of temperature and a non-equilibrium sudden heating process on the structures of argon clusters in SWCNTs are also studied by molecular dynamics simulation.     

Initial conformation of kinesin's neck linker

How ATP binding initiates the docking process of kinesin＇s neck linker is a key question in understanding kinesin mechanisms. By exploiting a molecular dynamics method, we investigate the initial conformation of kinesin＇s neck linker in its docking process. We find that, in the initial conformation, the neck linker has interactions with /30 and forms a ＇cover-neck bundle＇ structure with/30. From this initial structure, the formation of extra turns and the docking of the cover- neck bundle structure can be achieved. The motor head provides a forward force on the initial cover-neck bundle structure through ATP-induced rotation. This force, together with the hydrophobic interaction of ILE327 with the hydrophobic pocket on the motor head, drives the formation of the extra turn and initiates the neck linker docking process. Based on these findings, a pathway from ATP binding-induced motor head rotation to neck linker docking is proposed.