Strain Rate Sensitivities of Face-Centred-Cubic Metals Using Molecular Dynamics Simulation

We use dislocation theory and molecular dynamics （MD） simulations to investigate the effect of atom properties on the macroscopic strain rate sensitivity of f cc metals. A method to analyse such effect is proposed. The stress dependence of dislocation velocity is identified as the key of such study and is obtained via 2-D MD simulations on the motion of an individual dislocation in an fcc metal. Combining the simulation results with Orowan＇s relationship, it is concluded that strain rate sensitivities of fcc metals are mainly dependent on their atomic mass rather than the interatomic potential. The order of strain rate sensitivities of five fcc metals obtained by analysing is consistent with the experimental results available.     

Bridging Atomistic/Continuum Scales in Solids with Moving Dislocations

We propose a multiscale method for simulating solids with moving dislocations. Away from atomistic subdomains where the atomistic dynamics are fully resolved, a dislocation is represented by a localized jump profile, superposed on a defect-free field. We assign a thin relay zone around an atomistic subdomain to detect the dislocation profile and its propagation speed at a selected relay time. The detection technique utilizes a lattice time history integral treatment. After the relay, an atomistic computation is performed only for the defect-free field. The method allows one to effectively absorb the fine scale fluctuations and the dynamic dislocations at the interface between the atomistic and continuum domains. In the surrounding region, a coarse grid computation is adequate.     

A Simple Theoretical Method to Predict the Hardness of Pure Metal Crystals

Design of superhard bulk materials requires predicting their hardness, challenging current theories for material design. By introducing a concept of condensing force （CF）, it is shown via ab initio calculations for fcc （Ni, Cu, Al, Ir, Rh, Au, Ag, Pd） and hcp Re crystals that materials with larger CF can have greater hardness. Since the calculation of CF is easy, this method might prove a convenient way to evaluate the hardness of newly designed materials.     

Theoretical Strength of Face-Centred-Cubic Single Crystal Copper Based on a Continuum Model

The constitutive relation of single crystal copper based on atomistic potential is implemented to capture the nonlinear inter-atomic interactions. Uniaxial loading tests of single crystal copper with inter-atomic potential finite-element model are carried out to determine the corresponding ideal strength using the modified Born stability criteria. Dependence of the ideal strength on the crystallographic orientation is studied, and tension- compression asymmetry in ideal strength is also investigated. The results suggest that asymmetry for yielding strength of nano-materials may result from anisotropic character of crystal instability. Moreover, the results also reveal that the critical resolved shear stress in the direction of slip is not an accurate criterion for the ideal strength since it could not capture the dependence on the loading conditions and hydrostatic stress components for the ideal strength.     

Electronic structure of Sc2AC (A=Al, Ga, In, Tl)

Using ab initio calculations, we have studied Sc₂AC with A=Al, Ga, In and Tl. We show that C 2p and Sc 3d as well as A p and Sc 3d states are hybridized, but the antibonding states in the vicinity of the Fermi level weaken the overall bonding. In terms of the chemical bonding, the influence of the size of the A element is minute. Furthermore, the bulk modulus of the corresponding binary transition metal carbide is not conserved in these phases. Therefore, Sc₂AC can be classified as weakly coupled MAX phases according to Sun and co-workers [Z. Sun, D. Music, R. Ahuja, S. Li, J.M. Schneider, Phys. Rev. B 70 (2004) 092102]. It is our ambition that these calculations will stimulate experimental research on these compounds.     

Antisite Defects of the L12 Structure Determined by the Phase Field Microelasticity Model

A phase field microelasticity simulation is performed to examine the antisite defect of L12-Ni3Al in Ni75Al5.3 V19.7 ternary alloy. Combinimg strain energy with the phase field model leads to an atom configuration change as time proceeds. For the Ni sublattice, the antisite defect AlNi, the equilibrium occupancy probability （OP） of which declines, precedes NiNi and VNi in reaching equilibrium; subsequently, NiNi and VNi present a phenomenon of symmetrical rise and decline individually. Similarly, for the Al sublattice, the antisite defect NiAl, the OP of which eventually rises, takes fewer time steps than AlAl and VAl to attain equilibrium. Thereafter, AlAl rises while VAl declines symmetrically at the axes of the NiAl curve. Furthermore, the OP for the Al sublattice is much more sensitive to strain energy than that for the Ni sublattice.     

Ab initio study of Nb2SC and Nb2S2C: Differences in coupling between the S and Nb-C layers

Using ab initio calculations, we have studied the structurally related compounds Nb₂SC and Nb₂S₂C. In Nb₂S₂C (space group , prototype Bi₂Te₃), S atoms are nearest neighbours, while in Nb₂SC (space group P6₃/mmc, prototype Cr₂AlC) this is not the case. The calculated equilibrium volume for these two phases deviates by 1.6-3.7% to previously-published experimental data and the bulk modulus-to-c₄₄ ratios obtained are 1.5 and 5.9, respectively. These results indicate a resemblance of Nb₂S₂C to hexagonal BN and graphite. Furthermore, we have demonstrated that the uniform compression method is adequate for estimating the elastic properties of Nb₂SC, a so-called MAX phase. It is our ambition that these calculations will stimulate further experimental research on these compounds.     

From Discreteness to Continuity： Dislocation Equation for Two-Dimensional Triangular Lattice

A systematic method from the discreteness to the continuity is presented for the dislocation equation of the triangular lattice. A modification of the Peierls equation has been derived strictly. The modified equation includes the higher order corrections of the discrete effect which are important for the core structure of dislocation. It is observed that the modified equation possesses a universal form which is model-independent except the factors. The factors, which depend on the detail of the model, are related to the derivatives of the kernel at its zero point in the wave-vector space. The results open a way to deal with the complicated models because what one needs to do is to investigate the behaviour near the zero point of the kernel in the wave-vector space instead of calculating the kernel completely.     

Crystal Structure of β-La2Mo2O9 from First Principles Calculation

Arrangements of O ions in β-La2Mo2O9 are studied by first principles calculation with two different calculation schemes. All final structure configurations consist of MoO4-tetrahedra, MoO5-hexahedra, LaO8 and LaO7 polyhedra. Molybdenum polyhedra are isolated from each other, lanthanum polyhedra are connected together by sharing O ions. The occupancies of three crystallographic distinct O sites O（1）, O（2） and O（3） are 100%, 91.7% and 25%, respectively, consistent with experiments. All configurations are related to each other by one of 12 symmetry operations of P213 space group, suggesting that the structure observed experimentally may be interpreted as a time and spatial average of these local or inherent structures.     

A Modified Free Volume Model for Characterizing of Rate Effect in Bulk Metallic Glasses

We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses （BMGs）. A dimensionless Deborah number DeiD = tr/ti is proposed to characterize the rate effect in BMGs, where tr is the structural relaxing characteristic time of BMGs under shear load, ti is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.     

Structural phase transition (zincblende-rocksalt) and elastic properties in AlY (Y=N, P and As) compounds: Pressure-induced effects

The present paper addresses the pressure-induced structural aspects of ZnS-type (B3) to NaCl-type (B1) structure in AlY (Y=N, P, As). An effective-interionic interaction potential (EIoIP) with long-range Coulomb and three-body interactions and the Hafemeister-and-Flygare-type short-range overlap repulsion extended up to the second-neighbour ions and the van der Waals (vdW) interaction is developed. Emphasis has been given on evaluating the vdW coefficients by the Slater-Kirkwood variational method, as both the ions are polarizable. The lattice model calculations have revealed reasonably good agreement with the available experimental data on the phase-transition pressures (P_t=16, 14, 7.5 GPa) and the elastic properties of AlY (Y=N, P, As). The equation of state curves (plotted between V(P)/V(0) and pressure) for both the B3 and B1 structures obtained are in fairly good agreement with the experimental results. The calculated values of the volume collapses [ΔV(P)/V(0)] are also close to their observed data. Further, the variations of the second-order elastic constants with pressure follow a systematic trend that is almost identical to that exhibited by the observed data measured for other semiconducting compounds with B3→B1 structural phase transitions.     

Indentation Load Effect on Young＇s Modulus and Hardness of Porous Sialon Ceramic by Depth Sensing Indentation Tests

Depth sensing indentation （DSI） tests at the range of 200-1800mN are performed on porous sialon ceramic to determine the indentation load on Young＇s modulus and hardness values. The Young modulus and hardness （Dynamic and Martens） values are deduced by analysing the unloading segments of the DSI test load-displacement curves using the Oliver-Pharr method. It is found that Young＇s modulus ET, the dynamic hardness HD and the Martens hardness HM exhibit significant indentation load dependences. The values of Young＇s modulus and hardness decrease with the increasing indentation load, as a result of indentation load effect. The experimental hf /hm ratios lower than the critical value 0.7, with hm being the maximum penetration depth during loading and hf the final unloading depth, indicate that our sample shows the work hardening behaviour.     

Dislocation emission from an elliptically blunted crack tip with surface effects

The explicit expressions for critical stress intensity factors are derived for edge dislocation emission from an elliptically blunt crack with surface effects under mode I and mode II loadings. The influence of surface effects on dislocation emission criterion is analyzed. The result indicates the impact of the surface stress becomes remarkable for nanoscale blunted cracks and some particular materials, which cannot only affect the value of the critical stress intensity factors for dislocation emission, but also alter the emission angle.     

Crystallographic symmetry effect on the nucleation in non-equilibrium aggregation pattern

Some non-equilibrium patterns formed in different systems are composed by crystallites with correlations in their crystallographic orientation. This reflects the fact that crystallographic symmetry controls the successive nucleation in the crystallite aggregating process. Based on the crystallographic symmetry, a new classification to describe the different kinds of anisotropy for non-equilibrium aggregation pattern has been proposed: parallel aggregation, twin aggregation and radial (random) aggregation. This new classification can clarify the different anisotropy related to the crystal point group, which the previous classifications cannot. The nucleation mechanisms for these three kinds of patterns are also discussed based on the symmetric relationship among the crystallites. It is shown that the crystallographic symmetry plays an important role not only on the crystal polyhedral shape formed by growth, but also on the aggregation pattern formed by nucleation.     

Interaction of Point Defects with Twin Boundaries in Copper

The interaction between small vacancy clusters and twin boundaries in copper is studied by using many-body potential developed by Ackland et aL for fcc metals. The interaction energies of single-, di- and tri-vacancy clusters with （111） and （112） twin boundaries are computed using well established simulation techniques. For （111） twins the vacancy clusters are highly repelled when they are on the adjacent planes, and are attracted when they are away from the boundary. In the case of （112） twins, vacancy clusters are more attracted to the boundary when they are near the boundary as compared to away from it. Vacancy clusters on both the sides of the boundary are also investigated, and it is observed that the clusters energetically prefer to lie on the off-mirror sites as compared to the mirror position across the twin.     

First principles calculations of structural phase transformation in CaTe at high pressure

First principles calculations of the total energy of CaTe as a function of unit cell volume have been carried out for the NaCl, MnP and CsCl structures on the basis of density functional theory (DFT). All these calculations are performed with the CRYSTAL06 program package. The sequence of high-pressure phases for CaTe transforms from NaCl phase to an intermediate state with a mixture of NaCl and MnP phases and then to the CsCl phase is obtained, which is in good agreement with the previous experimental results. Several structural properties (equilibrium lattice constant, bulk modulus, etc.) of NaCl structure have been calculated, which are also in agreement with the previous experimental results.     

First-principles electronic-structure calculation on the spin distribution and the half-metallic properties of the mixed valence iron compound Ba2F2Fe1.5S3

The first principles within the full potential linearized augmented plane wave (FP-LAPW) method with the generalized gradient approximation (GGA) approach were applied to study the new mixed valence compound Ba₂F₂Fe_1.5S₃. The density of states, the electronic band structure and the spin magnetic moment are calculated. The calculations reveal that the compound has an antiferromagnetic interaction between the Fe^III and Fe^II ions arising from the bridging S atoms, which validate the experimental assumptions that there is a low-dimensional antiferromagnetic interaction in Ba₂F₂Fe_1.5S₃. The spin magnetic moment mainly comes from the Fe^III and Fe^II ions with smaller contribution from S anion. By analysis of the band structure, we find that the compound has half-metallic property.     

Ab-initio calculation of stability and structural properties of cadmium chalcogenides CdS, CdSe, and CdTe under high pressure

The structural phase transformations of Cd chalcogenide compounds (CdS, CdSe, and CdTe) under high pressure are studied using the local approximation to the density functional theory, and the one-electron equations are solved by means of the full potential linear muffin-tin-orbital method. CdS, CdSe, and CdTe are found to have nearly similar behaviour of the structural properties under high pressure. In the CdS compound, thePmmn phase is predicted to be thermodynamically stable after the rock-salt structure, and theCmcm structure is found to be thermodynamically stable before the rock-salt structure in both CdSe and CdTe. The resulting structural properties of the zinc-blende, wurtzite, cinnabar, rock-salt,Pmmn, CsCl,Cmcm, and β-Sn phases for these considered compounds are found to be in agreement with the existing experimental data.     

Electronic structure and site occupation for the intermediate phase of LaNi4.5Al0.5Hy

The crystal structure, electronic structure and hydrogen site occupancy of LaNi_4.5Al_0.5H_y intermediate phase (y=2.0, 3.0, 4.0) have been investigated using the full-potential linearized augmented plane wave (FLAPW) method. For the first time we analyzed the interstitial site occupation of hydrogen atoms. The H atoms were found to prefer the 6m, 3f and 12n sites, while no 4h sites were occupied. A narrowed Ni-d band is found due to the lattice expansion: the total DOS at E_F increases with y, which indicates that the compounds become less stable. The interaction between Al and Ni, with H plays a dominant role in the stability of LaNi_4.5Al_0.5H_y intermediate phase. The smaller the shift of E_F towards the higher energy region, the more stable the compounds will be. Our results are compared with experimental data and discussed in the light of previous works.     

Ab Initio Study on Hypothetical Silver Nitride

We perform the ab initio calculations based on norm-conserving pseudopotentials and density functional theory to investigate the structural, elastic, and thermodynamical properties for silver nitride （AgN） compound that is a member of the 4d transition metal group and has not been synthesized yet. The obtained results are compared with the other available theoretical data, and the agreement is, generally, quite good. We also present the pressure-dependent behaviour of some mechanical and thermodynamical properties for the same compounds.