星系团Abell2199的动力学子结构

近年来Sloan巡天的光谱观测使近邻富星系团Abell 2199的成员星系数目增加了近两倍.基于该星系团1.5 h-1Mpc范围内的366颗成员星系的位置和红移信息,对其动力学结构进行了分析.计算出中心cD星系的本动速度为(189±45)km s-1,说明cD星系在星系团完全维力化之前就形成了.从这些星系的空间分布和局域视向速度分布中,发现了三个明显的子团.这些子团的并合现象表明,星系团Abell 2199远未达到动力学平衡状态,从而支持了星系团形成的等级模型.

星系团Abell 1795的形态和动力学研究

基于SDSS和NED关于星系团Abell 1795的1.5 h-1Mpc范围内154个成员星系的位置和红移信息,对该团的动力学结构和形态进行了分析.计算出中心cD星系的本动速度为235±72 km.s-1,说明cD星系在星系团完全维里化之前就已经形成.从这些星系的空间投影分布和局域视向速度分布中,发现了正在并合之中的一个子团,进一步说明了该星系团并未处于动力学平衡状态.该星系团中心区域早型星系的比例显著偏高,说明中心区域的恒星形成率偏低.外围区域晚型星系所占的比例略高,可能预示着气体、尘埃等的存在,表现出恒星形成的活跃.团星系的空间延展方向和中心cD星系的方位角一致,在观测上支持了等级成团理论.

Structure-controlled slow dynamics in Al−Mg melts

Molecular dynamics simulation was employed to investigate the dynamical and structural properties of Al−Mg melts with the Al concentration systematically changed.The results show that the viscosity of Al67Mg33 abnormally surpasses that of Al85Mg15 below 550 K,inconsistent with the tendency at high temperatures.The evolution of the icosahedral order population is found to account for this dynamic behavior.Structural analysis shows a preferred bonding between Al and Mg atoms in the nearest neighbor shells,while a repelling tendency between them in the second shells,leading to the prepeak emergence in the partial static structure factors.The formation of icosahedral clusters is constrained in the Al-rich compositions because of the lack of sufficient Mg atoms to stabilize the clusters geometrically.These results demonstrate the structural consequence through the interplay between geometric packing and chemical interaction.These findings are crucial to understanding the structure−dynamic properties in Al−Mg melts.

用新的测光红移方法研究星系团Abell 2199

Sloan数字巡天第4批数据的释放（SDSSDR4）提供了新的星系团观测数据。采用最新研究的测光红移技术对Abell 2199没有光谱红移的星系进行了测光红移估计，挑选出了一批暗的成员星系候选体，使得样本的数量比单独用NED的样本增加了近1／4。通过对样本中470个成员星系运用k-test所得到的结果分析得出，原来确定的各子团星系的数目有所增加，但biweight位置和尺度与原来研究的各子团biweight位置和尺度值相比没有太大的变化。通过对团星系的空间分布和方位角的估算，可看出从核区到～5．5。的特大尺度，由于引力作用而显示出线性排列方式，揭示了一个延展到8．0h^-1Mpc由星系或暗物质组成的纤维状结构，进一步支持了星系形成的等级结构模型。

Molecular dynamics simulation of relationship between local structure and dynamics during glass transition of Mg_7Zn_3 alloy

The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated. It was found that the Mg-centered FK polyhedron and the Zn-centered icosahedron play a critical role in the formation of Mg7Zn3 metallic glass. The self-diffusion coefficients of Mg and Zn atoms deviate from the Arrhenius law near the melting temperature and then satisfy the power law. According to the time correlation functions of mean-square displacement, incoherent intermediate scattering function and non-Gaussian parameter, it was found that the β-relaxation in Mg7Zn3 supercooled liquid becomes more and more evident with decreasing temperature, and the α-relaxation time rapidly increases in the VFT law. Moreover, the smaller Zn atom has a faster relaxation behavior than the Mg atom. Some local atomic structures with short-range order have lower mobility, and they play a critical role in the appearance of cage effect in theβ-relaxation regime. The dynamics deviates from the Arrhenius law just at the temperature as the number of local atomic structures begins to rapidly increase. The dynamic glass transition temperature （Tc） is close to the glass transition point in structure （TgStr）.

Local structure of calcium silicate melts from classical molecular dynamics simulation and a newly constructed thermodynamic model

The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.

Atomistic simulation of defected magnesium hydroxide as flame retardants

The mechanical properties and the point defect energy of magnesium hydroxide（Mg（OH）2） were studied using the molecular dynamics. Moreover, the microelectronic structure of Mg（OH）2 with point defects in the bulk and on its surface were investigated using the first principles. The simulation results indicate that Mg（OH）2 was easily modified by other cations because of its strong, favorable interstitial and substitution defects via point defect energy calculation. Mg（OH）2 can provide high-efficiency flame retardancy because of the strong OH（OH Schottky defect） or H bond（H Frenkel defect and Schottky defect）. The potential model of Mg（OH）2 was established, and molecular dynamics simulation was used to investigate the relations between the crystal structure and the mechanical properties. Mg（OH）2 with special morphology such as nano-sheets was a prior consideration to maintain the composite mechanical properties. The detailed electronic structures of Mg（OH）2 with defects were determined. This work may provide theoretical guidance for choosing dopant element and reveal the element doping mechanism of Mg（OH）2.

Mesoscopic Simulation of Aggregate Behaviour of Polymers in β-HMX-based PBXs

The mesoscopic structures of β-HMX （octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine）- based PBXs （polymer bonded explosives） at room temperature were investigated using dissipative particle dynamics method. The parameters and repulsive parameters of dif- ferent polymers and β-HMX, the mesoscopic structures of β-HMX-based polymer-bonded explosives at different temperatures have been studied. The results showed that the compat-ibility between β-HMX and vinylidenedifluoride （VDF）,β-HMX and chlorotrifluoroethylene （CTFE）, VDF and CTFE increased with increasing temperature. The temperature and mo-lar ratio of the polymers played an important role in wrapped process. And there exists the optimum temperature and molar ratio.

Simulation study on non-linear effects of initial melt temperatures on microstructures during solidification process of liquid Mg_7Zn_3 alloy

The non-linear effects of different initial melt temperatures on the microstructure evolution during the solidification process of liquid Mg7Zn3 alloys were investigated by molecular dynamics simulation, The microstructure transformation mechanisms were analyzed by several methods. The system was found to be solidified into amorphous structures from different initial melt temperatures at the same cooling rate of 1×10^12 K/s, and the 1551 bond-type and the icosahedron basic cluster （12 0 12 0 ） played a key role in the microstructure transition. Different initial melt temperatures had significant effects on the final microstructures. These effects only can be clearly observed below the glass transition temperature Tg; and these effects are non-linearly related to the initial melt temperatures, and fluctuated in a certain range. However, the changes of the average atomic energy of the systems are still linearly related with the initial melt temperatures, namely, the higher the initial melt temperature is, the more stable the amorphous structure is and the stronger the glass forming ability will be.

Molecular Dynamics Study on Microstructure of Potassium Dihydrogen Phosphates Solution

Molecular dynamics simulations were carried out to study the internal energy and microstructure of potassium dihydrogen phosphates （KDP） solution at different temperatures. The water molecule was treated as a simple-point-charge model, while a seven-site model for the dihydrogen phosphate ion was adopted. The internal energy functions and the radial distribution functions of the solution were studied in detail. An unusually large local particle number density fluctuation was observed in the system at saturation temperature. It has been found that the specific heat of oversaturated solution is higher than that of unsaturated solution, which indicates the solution experiences a crystallization process below saturation temperature. The radial distribution function between the oxygen atom of water and the hydrogen atom of the dihydrogen phosphate ion shows a very strong hydrogen bond structure. There are strong interactions between potassium cation and oxygen atom of dihydrogen phosphate ion in KDP solution, and much more ion pairs were formed in saturated solution.

Molecular Dynamics Simulations of Silica Nanotube： Structural and Vibrational Properties Under Different Temperatures

Four-, six-, and eight-membered ring silica nanotubes at temperatures from 300 K to 1600 K are relaxed by classical molecular dynamics simulations with three potential models. The simulation results indicate that the stability of the end rings of the three silica nanotubes gradually decreases with increase in temperature. The validity of the vibrational features of silica nanotubes is shown by the vibrational density of states. Infrared spectra on the silica nanotubes under different temperatures are investigated. A detailed assignment of each spectral peak to the corresponding vibrational mode of the three nanotubes has been addressed. The results are in good agreement with the other theoretical and experimental

Melting Behaviour of Shell-symmetric Aluminum Nanoparticles： Molecular Dynamics Simulation

Molecular dynamics simulations with embedded atom method potential were carried out for A1 nanoparticles of 561 atoms in three structures： icosahedron, decahedron, and truncated octahedron. The total potential energy and specific heat capacity were calculated to estimate the melting temperatures. The melting point is 540＋10 K for the icosahedral structure, 500±10 K for the decahedral structure, and 520±10 K for the truncated octahedral structure. With the results of mean square displacement, the bond order parameters and radius of gyration are consistent with the variation of total potential energy and specific heat capacity. The relaxation time and stretching parameters in the Kohlraush-William-Watts relaxation law were obtained by fitting the mean square displacement. The results show that the relationship between the relaxation time and the temperatures is in agreement with standard Arrhenius relation in the high temperature range.

SEMI-EMPIRICAL CALCULATION FOR ELECTRONIC STRUCTURE OF C60 CLUSTER BY MOLECULAR DYNAMICS AND MNDO APPROACH

The electronic structure for C 60 was semi empirically investigated by using MD (molecular dynamics) and MNDO (modified neglect of diatomic overlap) approach of quantum chemistry.Especially,taking both σ and π orbitals into account,one electron energy levels,those symmetries and π orbital occupancies as well as electron excitation energies for different select rules,cohesive energy,ionization energies and electronic affinity forces were calculated.The obtained molecular orbital ratio shows a wide separation of σ and π types,and near HOMO and LUMO levels there are π orbitals mainly.The calculated semi empirical calculation results are in good agreement with experimental and ab initio calculation data.

Local structure changes of Cu55 cluster during heating

The structural relaxation of a cluster containing 55 atoms at elevated temperatures is simulated by molecular dynamics. The interatomic interactions are given by using the embedded atom method （EAM） potential. By decomposing the peaks of the radial distribution functions （RDFs） according to the pair analysis technique, the local structural patterns are identified for this cluster. During increasing temperature, structural changes of different shells determined by atom density profiles result in an abrupt increase in internal energy. The simulations show how local structural changes can strongly cause internal energy to change accordingly.

Molecular Dynamics Simulations of Liquid Phosphorus at High Temperature and Pressure

By performing ab initio molecular dynamics simulations, we have investigated the microstructure, dynamical and electronic properties of liquid phosphorus （P） under high temperature and pressure. In our simulations, the calculated coordination number （CN） changes discontinuously with density, and seems to increase rapidly after liquid P is compressed to 2.5 g/cm3. Under compression, liquid P shows the first-order liquid-liquid phase transition from the molecular liquid composed of the tetrahedral P4 molecules to complex polymeric form with three-dimensional network structure, accompanied by the nonmetal to metal transition of the electronic structure. The order parameters Q6 and Q4 are sensitive to the microstructural change of liquid P. By calculating diffusion coefficients, we show the dynamical anomaly of liquid P by compression. At lower temperatures, a maximum exists at the diffusion coefficients as a function of density; at higher temperatures, the anomalous behavior is weakened. The excess entropy shows the same phenomena as the diffusion coefficients. By analysis of the angle distribution functions and angular limited triplet correlation functions, we can clearly find that the Peierls distortion in polymeric form of liquid P is reduced by further compression.

Orientation and Motion of Water Molecules at Air/Water Interface

Here we report a quantitative study of the orientational structure and motion of water molecule at the air/water interface. Analysis of Sum Frequency Generation （SFG） vibrational peak of the free O-H stretching band at 3700 cm^-1 in four experimental configurations showed that orientational motion of water molecule at air/water interface is libratory within a limited angular range. The free OH bond of the interracial water molecule is tilted around 33°from the interface normal and the orientational distribution or motion width is less than 15°. This picture is significantly different from the previous conclusion that the interracial water molecule orientation varies over a broad range within the ultrafast vibrational relaxation time, the only direct experimental study concluded for ultrafast and broad orient, ational motion of a liquid interface by Wei et al. （Phys. Rev. Lett. 86, 4799, （2001）） using single SFG experimental configuration.

Density functional investigation on structural and electronic properties of small bimetallic Pb_(n)Ag_(n)(n=2–12)clusters

Structural and electronic properties of PbnAgn(n=2–12)clusters were investigated by density functional theory with generalized gradient approximation at BLYP level in DMol3 program package.The optimized bimetallic PbnAgn(n=2–12)clusters were viewed as the initial structures,then,those were calculated by ab initio molecular dynamics(AIMD)to search possible global minimum energy structures of PbnAgn clusters,finally,the ground state structures of PbnAgn(n=2–12)clusters were achieved.According to the structural evolution of lowest energy structures,Ag atoms prefer gather in the central sites while Pb atoms prefer external positions in PbnAgn(n=2–12)clusters,which is in excellent agreement with experimental results from literature and the application in metallurgy.The average binding energies,HOMO-LUMO gaps,vertical ionization potentials,vertical electron affinities,chemical hardnessη,HOMO orbits,LUMO orbits and density of states of PbnAgn(n=2–12)clusters were calculated.The results indicate that the values of HOMO-LUMO gaps,vertical ionization potentials,vertical electron affinities and chemical hardnessηshow obvious odd-even oscillations when n≤5,PbnAgn(n=2–12)clusters become less chemically stable and show insulator-to-metal transition with the variation of cluster size n,PbnAgn(n≥9)cluster are good candidates to study the properties of PbAg alloys.Those can be well explained by the density of states(DOS)distributions of Pb atoms and Ag atoms between–0.5 Ha and 0.25 Ha in PbnAgn(n=2–12)clusters.

Theoretical Study of Hydration Effects on the Selectivity of 18-Crown-6 Between K~＋ and Na~＋

A combination of molecular dynamics （MD） and density functional theory （DFT） calculations were used to study the hydration structures of K＋ and Na＋ ions under the confinement of 18-crown-6 in order to identify the role of water in the selectivity of 18-crown-6 towards K＋.The radial distribution functions,coordination num-bers,orientation structures and interaction energies were analyzed to investigate the hydration of K＋ and Na＋ in 18-crown-6/cation complexes.All calculations of K＋ and Na＋ in bulk water were also conducted for comparison.The simulation results show that the orientation distributions of the water molecules in the first coordination shell of K＋ are more sensitive to the confinement of 18-crown-6 than those of Na＋.It is more favorable to confine a K＋ in 18-crown-6 than a Na＋ in terms of interaction energy.Good agreement is obtained between MD results and DFT results.

Molecular Dynamics Simulation of Temperature-dependent Flexibility of Thermophilic Xylose Isomerase

The complex model of Thermus thermophilus xylose isomerase （TtXI） with D-xylose was constructed, and molecular dynamics （MD） simulations were carried out at 300 and 360 K for 10 ns by NAMD2.5. The radius of gyration （Rg）, subunit interactions, and residue flexibility were analyzed. The results show that residues 60-69, 142-148, 169-172, and 332-340 have high flexibility at 300 and 360 K. Residues with higher flexibility at 360 K than that at 300 K can mainly be divided into two groups： one locates in the helix-loophelix region consisting of residues 55-80 in catalytic domain; the other at subunit interfaces. The Rg of catalytic domain at 360 K shows 0.16 A higher than that at 300 K, but Rg of small C-terminal domain has no obvious difference. The results indicate that enhanced Rg of catalytic domain may lead to the intense motion of the active site of TtXI and promote the D-xylose isomization reaction. Eight hydrogen bonds and five ion pairs are reduced at subunit interfaces at 360 K compared with 300 K, that may be the main reason for the decrease in rigidity and increase in activity at high temperature of TtXI. This result also help to explain the cold-adaption phenomenon of TtXI E372G mutant reported previously. Our results reveal the relationship between temperature and structure flexibility of TtXI, and play an important role in understanding the thermostability of thermophile protein with multiple subunits.

Molecular dynamics simulations on the relationship between the elastic parameters and the molecular structures of nano-hybrid POSS materials

To research the relationship between the elastic parameters and the molecular structures of nano hybrid polyhedral oligomeric silsesquioxanes （POSS） materials, the mechanical properties at different temperatures for three POSS polymers with different molecular architectures, polymerlized norbornene POSS homopolymer （PNPOSS, pedant architecture）, γ- （2, 3 glycidoxy） propyl diaminoethane POSS polymer （GPDP, catena architecture） and trimethoxysilylcyelopentyl POSS polymer （ TSCP, cage - cage network architecture） were obtained by molecular dynamics simulations based on the Compass force-field. Results indicate that the moleculax architectures of the POSS polymers have great influence on the reinforced effects. The effect of the cage-cage network architecture is best, while that of the catena architecture takes second place and the pedant architecture has the least influence comparatively. The reinforced effects of the POSS monomers were examined. The influences of the temperatures on these effects were analyzed also. It may provide some basis for the reasonable applications of the excellent mechanical properties of the organic-inorganic nano-hybrid materials. It may also provide references for exploitation and design of the POSS materials.