Molecular-dynamics investigation on polarization retention of barium titanate nanofilm arising from ordered oxygen vacancy

Molecular-dynamics simulations have been carried out to investigate the electric hysteresis of barium titanate nanofilm containing oxygen vacancy ordering array parallel to the {101} crystal plane. The results obtained show a significant weakening of polarization retention from non-zero value to zero as the size of the array was reduced to a critical level, which was attributed to the formation and motion of head-to-head domain wall structure under external field loading process. By comparing with materials containing isolated oxygen vacancies, it was found that the zero retention was due to the oxygen vacancy ordering array rather than to the concentration of oxygen vacancy. Copyright (C) EPLA, 2010

Small-scale turbulent fluctuations beyond Taylor's frozen-flow hypothesis

The space-time cross-correlation function C-T(r, tau) of local temperature fluctuations in turbulent Rayleigh-Benard convection is obtained from simultaneous two-point time series measurements. The obtained C-T(r, tau) is found to have the scaling form C-T(r(E), 0) with r(E)=[(r-U tau)(2)+ V-2 tau(2)](1/2), where U and V are two characteristic velocities associated with the mean and rms velocities of the flow. The experiment verifies the theory and demonstrates its applications to a class of turbulent flows in which the requirement of Taylor's frozen flow hypothesis is not met.

Mechanical Behavior of Nanometer Ni by MD Simulation

The indention simulation of the crystal Ni is carried out by molecular dynamics technique (MD) to study the mechanical behavior at nanometer scales, the indenter tips with sphere shape is used. Some defects such as dislocations, point defects are observed. It is found that defects (dislocations, amorphous) nucleated is from local region near the pin tip or the sample surface. The temperature distribution of local region is analyzed and it can explain our MD simulation result.

Nano-size Effect of Interface Energy and Its Effect on Interface Fracture

An analytical model about size-dependent interface energy of metal/ceramic interfaces in nanoscale is developed by introducing both the chemical energy and the structure stain energy contributions. The dependence of interface energy on the interface thickness is determined by the melting enthalpy, the molar volume, and the shear modulus of two materials composing the interfaces, etc. The analytic prediction of the interface energy and the atomic scale simulation of the interface fracture strength are compared with each other for Ag/MgO and Ni/Al2O3 interfaces, the fracture strength of the interface with the lower chemical interface energy is found to be larger. The potential of Ag/MgO interface related to the interface energy is calculated, and the interface stress and the interface fracture strength are estimated further. The effect of the interface energy on the interface strength and the behind mechanism are discussed.

Kinetic Studies of Gas Flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

Study on the Mechanism of Detonation to Quasi-detonation Transition

In this paper, the mechanism of detonation to quasi-detonation transition was discussed, a new physical model to simulate quasi-detonation was proposed, and one-dimensional theoretical and numerical simulation was conducted. This study firstly demonstrates that the quasi-detonation is of thermal choking. If the conditions of thermal choking are created by some disturbances, the supersonic flow is then unable to accept additional thermal energy, and the CJ detonation becomes the unstable quasi-detonation precipitately. The kinetic energy loss caused by this transition process is firstly considered in this new physical model. The numerical results are in good agreement with previous experimental observations qualitatively, which demonstrates that the quasi-detonation model is physically correct and the study are fundamentally important for detonation and supersonic combustion research.

Splitting the Fast and Slow Motions in Molecular Dynamics Simulations Based on the Change of Cold Potential Well Bottom

The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Heat Transfer Problems Induced By Multi-shocks Interaction

Numerical simulations of the multi-shock interactions observable around hypersonic vehicles were carried out by solving Navier-Stokes equations with the AUSMPW scheme and the new type of the IV interaction created by two incident shock waves was investigated in detail. Numerical results show that the intersection point of the second incident shock with the bow shock plays important role on the flow pattern, peak pressures and heat fluxes. In the case of two incident shocks interacting with the bow shock at the same position, the much higher peak pressure and more severe heat transfer rate are induced than the classical IV interaction. The phenomenon is referred to as the multi-shock interaction and higher requirements will be imposed on thermal protection systems.

A New High Order Accurate Shock Capture Method with Wave Booster

A new kind of shock capturing method is developed. Before applying the high order accurate traditional scheme which is called as base scheme in this paper the fluid parameters are preconditioned in order to control the group velocity. The newly constructed scheme is high order accurate, simple, has high resolution of the shock, and less computer time consumed.

Parallel Fully-Implicit Computation of Magnetohydrodynamics Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETS,: package. The results of convergence tests arc plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

香港地区湍蛙遗传多样性研究

依据线粒体上ND2和CO1两个变异较大的基因序列分析了香港地区香港湍蛙7种群、华南湍蛙1种群，以及大陆其他地区华南湍蛙7种群，戴云湍蛙1种群，武夷湍蛙1种群的系统发育关系，进而探讨香港湍蛙的遗传多样性、香港湍蛙特有性、如何确定香港湍蛙最佳保护单元以及这四种湍蛙的物种分类地位。 1． 香港湍蛙保护遗传学研究 香港湍蛙核苷酸传多样性较低，从其遗传多样性信息、单倍型网络分析、中性检验值以及岐点分布结果一致显示香港湍蛙很可能经历了瓶颈后的扩张，种群正在由一个较小的有效种群大小迅速增长, 有足够的时间通过变异用于积累单倍型的多态性, 而对于提高核苷酸多样化而言, 时间尚短（Nei M et al，1975，Avise J C，2000；李明等，2003）。 分子变异分析结果显示香港湍蛙种群间存在较多的基因交流，且系统发育树上各种群间交叉在一起，没有形成与地理单元相关的分支，而从其单倍型网络看，他们源于共同的祖先，是一个单系群，与地理单元间没有形成显著的遗传分化。因此应作为一个进化显著单元（ESU）。结合其与其他湍蛙发育关系及遗传距离以及野外采集信息认为香港湍蛙只在香港地区有分布，属于香港特有种。该物种内遗传多样性较低，又属于世界自然保护联盟红皮书中的近危种，同时也是《野生动物保护条例》中的受保护野生动物，且由于香港城市建设等使得其栖息环境受到威胁，因此在香港特别行政区应该受到重点保护。 从单倍型分布和核苷酸多样性可以看出大榄涌种群和城门种群具有较高的单倍型多样性和核苷酸多样性，应该作为保护的重点区域。 2. 华南湍蛙东、南沿海种群间系统关系 华南湍蛙分布广，各种群存在着丰富的遗传多样性信息且中部种群广西龙胜和湖南张家界种群核苷酸多样性明显高于其他边缘种群华南湍蛙。种群间几乎没有基因交流，且各种群间无共享单倍型，可见已形成了显著的遗传分化。各种群间遗传距离都较远，其中广东南昆山种群以及福建三港种群与其他种群距离最远，因此可以推测其他种群（广东深圳、香港大屿山、广西龙胜和防城以及湖南张家界种群）可能为独立进化的种群。但是否是一新种或一隐存种，还需要结合形态学进行更深入的研究。 本研究中无论从系统关系看还是从遗传距离看，大屿山种群与深圳种群最近，支持陈坚峰等将其定为华南湍蛙，即华南湍蛙新增一个分布点：香港大屿山。 系统树上广西防城种群（支B）与龙胜和湖南种群（支A）形成姐妹群。香港大屿山种群与深圳种群先形成姐妹群（支C），但却没有与其距离很近的广东南岭及南昆山种群（支D）形成姐妹群，可能粤北和粤中的环境及气候较复杂因此与粤南其他种群形成了明显的隔离。同时可以看出华南湍蛙种群遗传分化与地理距离没有显著的相关性。 3. 四种湍蛙间的系统关系 根据线粒体CO1基因建立四种湍蛙间的系统关系及其遗传距离，很清楚地看到，香港湍蛙与戴云湍蛙关系很近，而华南湍蛙则与武夷湍蛙较近。然而，戴云湍蛙同一个种群内部共有两个单倍型DY1和DY2，且两个单倍型间遗传距离大于DY1与香港湍蛙间遗传距离，更远远大于香港湍蛙种群内部的距离，即戴云湍蛙内部两个单倍型间遗传距离达到了种级水平，同样在系统发育树上这两个单倍型与香港湍蛙形成并系。但是，戴云湍蛙种内在形态上差异不显著。因此，其是否属于萌芽物种分化形成（budding speciation）或已经完全分化为两个不同的种值得进一步研究？ 与戴云湍蛙香港湍蛙关系类似，从系统树上看华南湍蛙不形成单系，而是分成两个大支，与武夷湍蛙形成并系，且福建和南昆山的华南湍蛙与武夷湍蛙遗传距离远大于武夷湍蛙种内福建种群与浙江种群的遗传距离，达到了种级分化水平。由此，可以推断武夷湍蛙是有效种。系统树上广东深圳、香港大屿山、广西防城和龙胜以及湖南张家界种群与华南湍蛙福建及南昆山各种群间遗传距离已超出了种内各种群间的遗传距离，但是至于这一支是否应为另外一个种，有必要扩大采样，并结合核基因及形态信息进行进一步研究。 MtDNA of ND2 and CO1 gene were used to investigate genetic diversity of Amolops in Hongkong .We collected seven populations of A. hongkongensis,，one population of A.ricketti from Hong Kong and other seven populations of A.ricketti from East and South of Chinese mainland. As well as one population of A. daiyunensis and one population of A.wuyiensis Phylogenetic relationship were analyzed of four species. Discussed whether A.hongkongensis is an endemic species and how can we make the conservation and management decisions. 1. Conservation Genetics of A. hongkongensis A. hongkongensis has a low nucleotide diversity, the results of genetic diversity, haplotype network, neutrality test and the mismatch distributions indicate that A. hongkongensis experienced a recent expansion after a bottle neck. They had enough time to accumulated haplotype diversity, but it’s too short to have a high nucleotide diversity（Nei M et al，1975，Avise J C，2000；Li et al，2003）. The result of AMOVA reveals that it has much gene exchange among the populations of A. hongkongensis. The clades of the phylogenetic tree were mixed together, no significant genetic differentiation among 8 populations and they share the same ancestor from the network analysis, these indicate that they are monophyly and should be protected as one ESU. Combined with the information of relationships of interspecies, genetic distance and distribution investigate, We conclude that A. hongkongensis is an endemic species of Hong Kong. Considering on the status of low genetic diversity in A.hongkongensis, and this species was listed in the IUCN red list as near threatened, as well as listed in the . Furthermore, it’s habitat loss and degradation more rapidly as the human activity got higher and higher. So it’s urgent to protect them in Hong Kong. Our results suggest that Tai Lam Wu and TAI MO Shan -Shing Mun populations have the higher priority to be protected because their higher genetic diversity. 2．Phylogenetic relationships among populations of Amolops ricketti from the Southern and eastern China A. ricketti has the considerable genetic diversity of mitochondrial haplotypes within and among populations, and Mitochondrial DNA diversity was higher in populations at the central area of the present distribution range of the frog，i. e. the Longsheng population and Zhangjiajie population, than at the edges of their distribution range. They have no share haplotype among populations, and have a significant genetic differentiation. Genetic distance is high among the populations, especially the distance of Nankun and Sangang group with others, which suggested that they evolved independently. May be there is a cryptic species or a new species, a further study is needed. The results of gene tree and the genetic distance clearly demonstrate that the population from LanTau island is A. ricketti, so we agree with Chen et al(2005) . That means A.ricketti have a new distribution site: LanTau island, HongKong. Phylogenetic relationships were analyzed through NJ and Mrbayes methods and got a consistent topological structure, the structure indicated that the ingroup were comprised four groups. Populations Longsheng and Zhangjiajie were first clustered as clade A; Populations Fangcheng was clustered together (clade B) as a sister group to clade A；Populations Shenzhen and Lantau island were sister groups and clustered as clade C；Then the clade D included populations Nankunshan and Nanling in Guangdong province and Sangang in Fujian province. 3. Phylogenetic Relationships among these four specises Phylogenetic relationships based on 1503bp CO1 gene and the genetic distance show that A. hongkongensis close to A. daiyunensis whereas A.ricketti near to A.wuyiensis. Nevertheless, there are two haplotypes in A.daiyunensis and the genetic distance between them higher than the distance between DY1 with A. hongkongensis. A. hongkongensis is nested in the paraphyletic ancestral species A. daiyunensis. Without significant difference in the morphological characters, So, we considered both A.daiyunensis and A.hongkongensis are valid species, may be this represents a case of ‘budding speciation’ like Batrachuperus pinchonii(Fu and Zeng,2008) in the population of A. daiyunensis. Just like two species above A. wuyiensis and A. ricketti are not monophyly, instead, A.wuyiensis is nested in the paraphyletic ancestral species A.ricketti. We need do more research to make sure whether they are new species.

High negative ion production yield in 30 keV F2+ + adenine (C5H5N5) collisions

In collisions between slow F2+ ions (30 keV) and molecular targets, adenine, scattered particle production yields have been measured directly by simultaneous detection of neutrals, positive and negative ions. The relative cross-section for a negative ion formation channel was measured to be 1%. Despite a slight decrease compared to a larger target, the fullerene C-60, the measured negative ion formation cross section is still at least one order of magnitude larger than the yield in ion-atom interactions.

Scaling and the thermal conductivity of the Frenkel-Kontorova model

We have studied the dependence of the thermal conductivity kappa on the strength of the interparticle potential lambda and the strength of the external potential beta in the Frenkel-Kontorova model. We found that the functional relation can be expressed in a scaling form, kappa(proportional to) lambda 3/2/beta(2 center dot). This result is first obtained by nonequilibrium molecular dynamics. It is then confirmed by two analytical methods, the self-consistent phonon theory and the self-consistent stochastic reservoirs method. The thermal conductivity kappa is therefore a decreasing functon of beta and an increasing function of lambda.

Ordered FePt Nanoparticle Arrays Prepared by a Micellar Method

City Univ Hong Kong

Study on Broadband Emitting Self-Assembled Quantum-Dot Material and Devices

City Univ Hong Kong