Elasticity, band structure, and piezoelectricity of BexZn1-xO alloys

Lattice constants, elasticity, band structure and piezoelectricity of hexagonal wideband gap BexZn1-xO ternary alloys are calculatedusing firstprinciples methods. The alloys' lattice constants obey Vegard's law well. As Be concentration increases, the bulk modulus and Young's modulus of the alloys increase, whereas the piezoelectricity decreases. We predict that BexZn1-xO/GaN/substrate (x = 0.022) multilayer structure can be suitable for high-frequency surface acoustic wave device applications. Our calculated results are in good agreement with experimental data and other theoretical calculations. (c) 2008 Elsevier B.V. All rights reserved.

压力和铁含量对铁方镁石(Mg1-xFex)O(x<=0.25)弹性常数影响的第一性原理研究

By first principle methods based on density functional theory (DFT),the equation of state(EOS) and elastic constants of both periclase and ferropericlase are calculated. The pressure and iron doping effects on the elastic constants of ferropericlase are investigated systematically. Firstly, we calculate the elastic constants of periclase and compare the obtained results with experimental data and other theoretical calculations, which shows a encouraging consistence and demonstrates the practicability of first-principle methods. Secondly, by adding iron into periclase crystal model, we build up ferropericlase with iron contents ranging from 0% to 25% mole percent. The corresponding elastic constants are calculated in a large pressure range(0~120GPa). Emphatically, the strong correlation of 3d electrons in transitional elements, such as iron, is difficult to treat in first-principle methods for a long time. The current solution is to make additional correction. During the initial stage of this study, the strong correlation of 3d electrons in iron is not considered, and we observed that addition of iron decreases the volume of ferropericlase, which is totally contradictory to the experimental data. By applying LDA+U approximation in order to solve the strongly correlated 3d electron of iron, we observed the expansion of volume by iron as expected. On the basis of the LDA+U approximation, the elastic constants of ferropericlase are calculated. After a detailed analysis of data obtained from theoretical calculations, we have reached the following conclusions:(1)pressure imposes positive effects on all elastic constants, and the degree of effects is C11>C12>C44. (2) Iron has no distinctive effects on C11 and C12, although some fluctuations are observed around 60GPa. However, iron has obvious softening effects on C44 The softening effects on C44 are intensified as pressure increases. Above the 100GPa, the effects increase greatly, even surpasses the pressure's positive effects in ferropericlase crystal models with iron mole percent of having 12.5%, 18.75% and 25% iron content. (3)As to the modulus deprived from elastic constants, iron has no effect on the adiabatic bulk module BS, only a little fluctuation around 60GPa. We find iron's softening effects on shear modulus G. (4)We find out that, compared with low iron content, elastic constants with iron content approaching 25mole% is consistently fluctuated，which may be caused by the limitations of the LDA+U approximation method itself. (5)We investigate the pressure and Fe doping effects on elastic anisotropy factor（A=(2C44+C12-C11)/C11） of ferropericlase and find out that iron contents will lower the critical isotropic pressure. At the same pressure, when the pressure is below the isotropic pressure, iron softens the anisotropy factor ; when pressure surpasses the isotropic pressure, iron increases the anisotropy factor.

Studies of magnetic interactions in Ni-doped ZnO from first-principles calculations

The magnetic interactions in Ni-doped ZnO are calculated using GGA and GGA + U method of density functional theory. The following three cases: (i) Ni-doped ZnO, (ii) (Ni, Al)-codoped ZnO, and (iii) (Ni, Li)-codoped ZnO are studied. The ferromagnetic ordering is always favorable for the three cases within GGA method. However, the ferromagnetic state is sometimes favorable after treating within the method of GGA + U. The GGA underestimates the correlated interactions especially when the Ni ions align directly to each other. (C) 2007 Elsevier B.V. All rights reserved.

ENERGY-LEVEL CALCULATIONS OF THE RECONSTRUCTED 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme taking into account 10 atomic orbitals (s-, p-, and d-type) applied to a supercell containing 256 atoms is used to calculate the bound states of the reconstructed 90-degrees partial dislocation in Si. The results differ significantly from our earlier calculations on the unreconstructed 90-degrees partial using the same method. We find two bands separate from each other in the entire Brillouin zone and the upper band penetrates deep into the indirect band gap which is in contradiction with the general opinion that core reconstruction clears the band gap of dislocation states.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

Ab initio calculations of differential cross sections for single charge transfer in He-3(2+)+He-4 collisions

The single charge transfer process in He-3(2+)+He-4 collisions is investigated using the quantum-mechanical molecular-orbital close-coupling method, in which the adiabatic potentials and radial couplings are calculated by using the ab initio multireference single- and double-excitation configuration interaction methods. The differential cross sections for the single charge transfer are presented at the laboratorial energies E = 6 keV and 10 keV for the projectile He-3(2+). Comparison with the existing data shows that the present results are better in agreement with the experimental measurements than other calculations in the dominant small angle scattering, which is attributed to the accurate calculations of the adiabatic potentials and the radial couplings.

Multiphase technique improves mud-pulse velocity calculations

Measurement while drilling (MWD) has become a popular survey technology to monitor directional data, drilling data, formation evaluation data and safety data in the world. And closed loop drilling shows promise in recent years. Obviously, the method of tr

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Numerical approximation improves well survey calculation

It is now possible to improve the precision of well survey calculations by order of magnitude with numerical approximation.

Although the most precise method of simulating and calculating a wellbore trajectory generally requires more calculation than other, less-accurate methods, the wider use of computers in oil fields now eliminates this as an obstacle.

The results of various calculations show that there is a deviation of more than 10 m among the different methods of calculation for a directional well of 3,000 m.¹ Consequently, it is important to improve the precision and reliability of survey calculation-the fundamental, necessary work of quantitatively monitoring and controlling wellbore trajectories.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Permeability of Fractal Porous Media by Monte Carlo Simulations

The permeability of the fractal porous media is simulated by Monte Carlo technique in this work. Based oil the fractal character of pore size distribution in porous media, the probability models for pore diameter and for permeability are derived. Taking the bi-dispersed fractal porous media as examples, the permeability calculations are performed by the present Monte Carlo method. The results show that the present simulations present a good agreement compared with the existing fractal analytical solution in the general interested porosity range. The proposed simulation method may have the potential in prediction of other transport properties (such as thermal conductivity, dispersion conductivity and electrical conductivity) in fractal porous media, both saturated and unsaturated.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Quantum and variational monte-carlo interaction potentials for li2 (x1-sigma-g+)

Optimized trial functions are used in quantum Monte Carlo and variational Monte Carlo calculations of the Li₂(X ¹Σ⁺_g) potential curve. The trial functions used are a product of a Slater determinant of molecular orbitals multiplied by correlation functions of electron—nuclear and electron—electron separation. The parameters of the determinant and correlation functions are optimized simultaneously by reducing the deviations of the local energy E_L (E_L  Ψ⁻¹_THΨ_T, where Ψ_T denotes a trial function) over a fixed sample. At the equilibrium separation, the variational Monte Carlo and quantum Monte Carlo methods recover 68% and 98% of the correlation energy, respectively. At other points on the curves, these methods yield similar accuracies.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

Elastic deformation of soft membrane with finite thickness induced by a sessile liquid droplet

In this paper, the role of vertical component of Surface tension of a droplet on the elastic deformation of a finite-thickness flexible membrane was theoretically analyzed using Hankel transformation. The vertical displacement at the Surface was derived and can be reduced to Lester's or Rusanov's solutions when the thickness is infinite. Moreover, some Simulations of the effect of a liquid droplet on a membrane with a finite thickness were made. The numerical results showed that there exists a saturated membrane thickness of the order of millimeter, when the thickness of a membrane is larger than such a value, the membrane can be regarded as a half-infinite body. Further numerical calculations for soft membrane whose thickness is far below the saturated thickness were made. By comparison between the maximum vertical displacement of an ultrathin soft membrane and a half-infinite body, we found that Lester's or Rusanov's solutions for a half-infinite body cannot correctly describe Such cases. In other words, the thickness of a soft membrane has great effect on the surface deformation of the ultrathin membrane induced by a liquid droplet. (C) 2009 Elsevier Inc. All rights reserved.