A new N-body potential and its application

Based on the embedded atom method (EAM) proposed by Daw and Baskes and Johnson's model, this paper constructs a new N-body potential for bcc crystal Mo. The procedure of constructing the new N-body potential can be applied to other metals. The dislocation emission from a crack tip has been simulated successfully using molecular dynamics method, the result is in good agreement with the elastic solution.

双晶晶界在高温下的热稳定性的分子动力学研究

晶界结构在高温下的热稳定性问题是一个长期争论而又未能解决的问题，其争论的焦点是：在远低于熔点的温度下，晶界结构是否发生了可观察到的无序化，即是否存在一个远低于熔点的结构转化温度。为了能澄清这一争论，本文系统地研究了晶界结构的热稳定性。为了消除相互作用势的影响和系统误差，本文首先采用Morse势和经验多体势分别对铝、铜单晶的熔化过程进行了分子动力学模拟。在平衡态下，通过计算表征结构无序化的静态结构因子、径向分布函数和单晶原子位形图，获得了铝、铜单晶的熔点，结果表明：多体势计算的铝和铜的单晶熔点更接近实验值。因此，采用经验多体势应用分子动力学方法分别模拟了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界晶界结构由有序向无序转化的过程，计算了平衡态下的表征结构无序化的静态结构因子、径向分布函数和晶界原子位形图并将多体势获得的铝、铜单晶熔点作为晶界结构转化温度的约化熔点，获得了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界结构的转化温度和熔点，结果表明：1．Σ5、Σ9、Σ11、Σ19、Σ33五种对称倾侧双晶晶界均在远低于单晶熔点温度时，晶界结构发生了可观察到的无序化，而且双晶晶界结构的转变温度相差不大，双晶晶界熔点也低于单晶熔点。2．Σ3晶界在温度远低于熔点时，其晶界结构没有发生可观察到的无序化；Σ3晶界的转化温度与单晶熔点接近。所以，可以认为Σ3晶界不存在转化温度。这是由于Σ3晶界为共格孪晶，具有较低的能量。综上所述，除Σ3共格孪晶外，在远低于熔点温度下，晶界结构发生了可观察到的无序化，即：存在一个远低于熔点的转化温度，此时其静态结构因子约为0.5左右；晶界结构的熔点均低于单晶熔点，此时其静态结构因子约为0.15左右。从全文模拟结果可以看出，静态结构因子、径向分布函数、晶界原子位形图三种方法在确定晶界的结构转化温度和熔点时，静态结构因子是最有效、最准确的定量方法。

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

Three-body force rearrangement contribution to single nucleon potential in nuclear matter

Within the framework of microscopic Brueckner-Hatree-Fock, the contribution of the three-body force (TBF) rearrangement to the. single nucleon potential is calculated. The TBF rearrangement effects on the momentum and the density dependence of the single nucleon potential are investigated. The influence of the TBF rearrangement on the effective mass of nucleon is also discussed. It is shown that the rearrangement contribution of TBF is repulsive and momentum-dependent. The TBF rearrangement effect and its momentum dependence increase rapidly as increasing density and momentum. At high densities and high momenta, the repulsive rearrangement contribution reduces strongly the attraction of the single nucleon potential and enhances considerably the momentum dependence of the single nucleon potential.

CONSISTENT THREE-BODY FORCE WITH BONN B POTENTIAL AND NEUTRON STAR STRUCTURE

Microscopic three-nucleon force consistent with the Bonn B two-nucleon potential is constructed, which includes Delta(1232), Roper, and nucleon-antinucleon excitation contributions. Recent results for the choice of the meson parameters are discussed. The forces are used in Brueckner calculations and the saturation properties of nuclear matter are determined. At the high densities,the nuclear equation of state and the symmetry energy are calculated. The corresponding neutron star mass-radius relations are presented.

Microcystin accumulation in freshwater bivalves from Lake Taihu, China, and the potential risk to human consumption

The potential risk through ingestion of microcystins (MC) in contaminated mollusks has not been well studied. The present paper studied seasonal changes of MC content (determined by liquid chromatography-mass spectrometry) in various organs of three species of bivalves (Cristaria plicata, Hyriopsis cumingii, and Lamprotula leai) in Lake Taihu, China, where toxic cyanobacterial blooms occurred. Coinciding with peaks of seston MC (maximum, 5.7 mu g/L) and MC in cyanobacterial blooms (maximum, 0.534 mg/g), most organs showed sharp MC peaks during the summer, indicating both fast uptake and fast depuration by bivalves. Because hepatopancreas and intestine had considerably higher MC content than other organs, they are the most dangerous for human consumption. Both the present and previous studies show that the hepatopancreatic MC and total tissue MC often are correlated in various aquatic invertebrates. During the peak of the cyanobacterial blooms, C. plicata had higher hepatopancreatic MC content than the other bivalves, whereas H. cumingii had higher intestinal MC content than the other bivalves. Estimated daily intakes for humans from the consumption of whole tissues of the three bivalves were 0.48 to 0.94 mu g MC-LR equivalent/kg body weight (12- to 23.5-fold the tolerable daily intake value proposed by the World Health Organization), which indicates a high risk for humans consuming these bivalves.

Hydroelastic dynamic characteristics of a slender axis-symmetric body

The slender axis-symmetric submarine body moving in the vertical plane is the object of our investigation. A coupling model is developed where displacements of a solid body as a Euler beam (consisting of rigid motions and elastic deformations) and fluid pressures are employed as basic independent variables, including the interaction between hydrodynamic forces and structure dynamic forces. Firstly the hydrodynamic forces, depending on and conversely influencing body motions, are taken into account as the governing equations. The expressions of fluid pressure are derived based on the potential theory. The characteristics of fluid pressure, including its components, distribution and effect on structure dynamics, are analyzed. Then the coupling model is solved numerically by means of a finite element method (FEM). This avoids the complicacy, combining CFD (fluid) and FEM (structure), of direct numerical simulation, and allows the body with a non-strict ideal shape so as to be more suitable for practical engineering. An illustrative example is given in which the hydroelastic dynamic characteristics, natural frequencies and modes of a submarine body are analyzed and compared with experimental results. Satisfactory agreement is observed and the model presented in this paper is shown to be valid.

Microscopic Three-Body Force Effect on Nucleon-Nucleon Cross Sections in Symmetric Nuclear Matter

We provide a microscopic calculation of neutron-proton and proton-proton cross sections in symmetric nuclear matter at various densities, using the Brueckner-Hartee-Fock approximation scheme with the Argonne V-14 potential including the contribution of microscopic three-body force. We investigate separately the effects of three-body force on the effective mass and on the scattering amplitude. In the present calculation, the rearrangement contribution of three-body force is considered, which will reduce the neutron and proton effective mass, and depress the amplitude of cross section. The effect of three body force is shown to be repulsive, especially in high densities and large momenta, which will suppress the cross section markedly.

Three-body force effect on P-3 F-2 neutron superfluidity in neutron matter, neutron star matter, and neutron stars

We investigate the effect of microscopic three-body forces on the P-3 F-2 neutron superfluidity in neutron matter, beta-stable neutron star matter, and neutron stars by using the BCS theory and the Brueckner-Hartree-Fock approach. We adopt the Argonne V18 potential supplemented with a microscopic three-body force as the realistic nucleon-nucleon interaction. We have concentrated on studying the three-body force effect on the P-3 F-2 neutron pairing gap. It is found that the three-body force effect considerably enhances the P-3 F-2 neutron superfluidity in neutron star matter and neutron stars.

(PF2)-P-3 neutron superfluidity in neutron stars and three-body force effect

We investigate the (PF2)-P-3 neutron superfluidity in beta-stable neutron star matter and neutron stars by using the BCS theory and the Brueckner-Hartree-Fock approach. We adopt the Argonne V-18 potential supplemented with a microscopic three-body force as the realistic nucleon-nucleon interaction. We have concentrated on studying the three-body force effect on the (PF2)-P-3 neutron pairing gap. It is found that the three-body force effect is to enhance remarkably the (PF2)-P-3 neutron superfluidity in neutron star matter and neutron stars.

Separation of potential and kinetic electron emission from Si and W induced by multiply charged neon and argon ions

T he total secondary electron emission yields, gamma(T), induced by impact of the fast ions Neq+ (q = 2-8) and Arq+ (q = 3-12) on Si and Neq+ (q = 2-8) on W targets have been measured. It was observed that for a given impact energy, gamma(T) increases with the charge of projectile ion. By plotting gamma(T) as a function of the total potential energy of the respective ion, true kinetic and potential electron yields have been obtained. Potential electron yield was proportional to the total potential energy of the projectile ion. However, decrease in potential electron yield with increasing kinetic energy of Neq+ impact on Si and W was observed. This decrease in potential electron yield with kinetic energy of the ion was more pronounced for the projectile ions having higher charge states. Moreover, kinetic electron yield to energy-loss ratio for various ion-target combinations was calculated and results were in good agreement with semi-empirical model for kinetic electron emission.