The influence of muffin-tin approximation on energy band gap

The influence of muffin-tin approximation on energy band gap was studied using LMTO-ASA (Linear Muffin-Tin Orbital-Atomic Sphere Approximation) approach. Since the diverse data are available for LaX(X=N, P, As, Sb), they are presented in our research as an example in order to test the reliability of our results. Four groups of muffin-tin radii were chosen, they were the fitted muffin-tin radii based on the optical properties of the crystals (the first), 1 : 1 for La : X(the second), 1.5 : 1 for La : X(the third), and a group of radii derived by making the charge in the interstitial space to be zero(the fourth). The results show that the fitted muffin-tin radii (the first group) give the best results compared with experimental values, and the predicted energy band gaps are very sensitive to the choice of muffin-tin radius in comparison with the other groups. The second and the third delivered results somewhere in between, while the fourth provided the worst results compared with the other groups. For the same crystal, with the increase of muffin-tin radius of lanthanum, the calculated energy band gaps decreased, going from semi-conductor to semimetal. This again clearly indicated the sensitivity of energy band structure on muffin-tin approximation.

The role of Li-O bonds in calculations of nonlinear optical coefficients of LiXO3-type complex crystals

The second-order nonlinear optical (NLO) tenser coefficients of LiXO3 (X = I; Nb or Ta) type complex crystals have been calculated using the chemical bond theory of complex crystals. Contributions of each type of bond to the total second-order NLO coefficient d(ij) and the linear susceptibility X are quantitatively determined. All tensor values thus calculated are in good agreement with experimental data. The Li-O bonds are found to be an important group in the contributions to the total NLO tenser coefficient, especially for those in LiNbO3 and LiTaO3. The importance of Li-O bonds depends on the environment of Li atom in these crystals.

Calculations of nonlinear optical responses of isomorphous crystals NaClO3 And NaBrO3 with natural optical activity

This work considers the isomorphous optically active crystals NaClO3 and NaBrO3. The connection between their second-order nonlinear optical (NLO) responses and chemical bond structures is established, starting from the experimental optical activities. The calculation reproduces the well-known experimental fact that crystals of NaClO3 and NaBrO3 with similar structures have different signs of optical rotation and of second harmonic generation (SHG). Unlike previous bond charge models, the method may include more than one type of bond in the calculation, and therefore may be used to study the optical activity and nonlinear optical properties of more general crystals. (C) 1998 Elsevier Science B.V. All rights reserved.

CALCULATIONS FOR THE VIBRATION FREQUENCIES OF THE O-H...O BOND IN KH2PO4 BY A SEMIEMPIRICAL METHOD

The frequencies of the stretching vibration and the bending vibration of the 0-H ... 0 bond in potassium dihydrogen phosphate have been calculated by means of two semiempirical formulae with three parameters. The calculated results can give satisfactory explanation for the experimental spectra of the potassium dihydrogen phosphate crystal. The parameters used in the calculations may be related to the chemical bonding and the charge distribution about the two oxygen atoms of the 0-H ... 0 bond system.

SIMULATION OF THE OXYGEN-ATOM VACANCY IN THE HIGH-TEMPERATURE SUPERCONDUCTOR TIBA2CAN-1CUNO2N+2.5(N = 1)

These simulation calculations for the oxygen-atom vacancy in the high temperature superconductor TlBa2Ca(n-1)Cu(n)O2n+2.5(n = 1) have been performed by means of the tight-binding approximation based on the EHMO method. The results indicate that the effect of the oxygen-atom vacancy on the charge distributions at the Tl-, Ba-, Cu- and O-atom sites is appreciably different and that there may exist two kinds of Cu cation with different net charges (approximately + 3.0 or approximately + 1.0) due to the oxygen-atom vacancy in the lattice. The electric field gradient at the site of the oxygen-atom vacancy has been calculated. The position of the oxygen-atom vacancy which favours the high temperature superconductivity of TlBa2Ca(n-1)Cu(n)O2n+2.5(n = 1) has been discussed.

CALCULATIONS OF THE EIGENVALUES OF PURE QUADRUPOLE INTERACTION IN MOSSBAUER-EFFECT STUDIES

Two approximate formulae to calculate the eigenvalues of pure quadrupole interaction in Mossbauer effect studies have been proposed and the eigenvalue coefficients in the formulae have been given for various excited states and ground states of the nucleus with different spin. All the eigenvalues of pure quadrupole interaction between both excited state and ground state of nucleus with spin I = 3/2 divided-by 9/2 and the electric-field gradient with different asymmetry parameter (eta = 0 divided-by 1.0) have been calculated by these formulae. The results show that the accuracies in all the calculations are more satisfactory or same in comparison with those obtained by the formula of Shenoy and Dunlap, especially when the asymmetry parameter of electric-field gradient is larger than 0.8 for the nucleus with spin I = 5/2.

Application of effective medium approximation theory to ocean remote sensing under wave breaking

Based on the effective medium approximation theory of composites, the empirical model proposed by Pandey and Kakar is remedied to investigate the microwave emissivity of sea surface under wave breaking driven by strong wind. In the improved model, the effects of seawater bubbles, droplets and difference in temperature of air and sea interface (DTAS) on the emissivity of sea surface covered by whitecaps are discussed. The model results indicate that the effective emissivity of sea surface increases with DTAS increasing, and the impacts of bubble structures and thickness of whitecaps layer on the emissivity are included in the model by introducing the effective dielectric constant of whitecaps layer. Moreover, a good agreement is obtained by comparing the model results with the Rose's experimental data.

压力和铁含量对铁方镁石(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.

全局优化的Fourier有限差分偏移方法研究

Seismic exploration is the main method of seeking oil and gas. With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in seismic exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which has obtained good effect. However, in complex media with wider angles, the effect of FFD method is not satisfactory. Based on the FFD operator, we extend the two coefficients to be optimized to four coefficients, then optimize them globally using simulated annealing algorithm. Our optimization method select the solution of one-way wave equation as the objective function. Except the velocity contrast, we consider the effects of both frequency and depth interval. The proposed method can improve the angle of FFD method without additional computation time, which can reach 75° in complex media with large lateral velocity contrasts and wider propagation angles. In this thesis, combinating the FFD method and alternative-direction-implicit plus interpolation(ADIPI） method, we obtain 3D FFD with higher accuracy. On the premise of keeping the efficiency of the FFD method, this method not only removes the azimuthal anisotropy but also optimizes the FFD mehod, which is helpful to 3D seismic exploration. We use the multi-parameter global optimization method to optimize the high order term of FFD method. Using lower-order equation to obtain the approximation effect of higher-order equation, not only decreases the computational cost result from higher-order term, but also obviously improves the accuracy of FFD method. We compare the FFD, SAFFD（multi-parameter simulated annealing globally optimized FFD）, PFFD, phase-shift method(PS）, globally optimized FFD （GOFFD）, and higher-order term optimized FFD method. The theoretical analyses and the impulse responses demonstrate that higher-order term optimized FFD method significantly extends the accurate propagation angle of the FFD method, which is useful to complex media with wider propagation angles.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.