2-DELTA-FUNCTION GENERALIZED PLASMA POLE MODEL FOR 1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES IN MANY-ELECTRON SYSTEMS

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-delta-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single delta-function GPP model has the drawback of over simplifications and for the crystals without the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-delta GPP model is an optimization of the two approaches. We analyze the two-delta GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these delta-functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.

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.

Friction phenomena and phase transition in the underdamped two-dimensional Frenkel-Kontorova model

Locked-to-sliding phase transition has been studied in the driven two-dimensional Frenkel-Kontorova model with the square symmetric substrate potential. It is found that as the driving force increases, the system transfers from the locked state to the sliding state where the motion of particles is in the direction different from that of driving force. With the further increase in driving force, at some critical value, the particles start to move in the direction of driving force. These two critical forces, the static friction or depinning force, and the kinetic friction force for which particles move in the direction of driving force have been analyzed for different system parameters. Different scenarios of phase transitions have been examined and dynamical phases are classified. In the case of zero misfit angle, the analytical expressions for static and kinetic friction force have been obtained.

Restudy on the wave functions of Gamma(nS) states in the light-front quark model and the radiative decays of Gamma(nS) -> eta(b) + gamma

The light-front quark model has been applied to calculate the transition matrix elements of heavy hadron decays. However, it is noted that using the traditional wave functions of the light-front quark model given in the literature, the theoretically determined decay constants of the Gamma(nS) obviously contradict the data. This implies that the wave functions must be modified. Keeping the orthogonality among the nS states and fitting their decay constants, we obtain a series of the wave functions for Gamma(nS). Based on these wave functions and by analogy with the hydrogen atom, we suggest a modified analytical form for the Gamma(nS) wave functions. Using the modified wave functions, the obtained decay constants are close to the experimental data. Then we calculate the rates of radiative decays of Gamma(nS) -> eta(b) + gamma. Our predictions are consistent with the experimental data on decays Gamma(3S) -> eta(b) + gamma within the theoretical and experimental errors.

Improvement of a Fission-Like Model for Nuclear alpha Decay

The alpha decay constant is the product of the penetrability P and assault frequency nu(0) in the fission-like model. An effective assault frequency P-nu replacing the previous assault frequency nu(0) is introduced for improvement of a fission-like model named the generalized liquid drop model (GLDM) to describe the nuclear alpha decay process more accurately. Two analytical formulae are proposed for the effective assault frequency due to experimental data within the GLDM. The improved model can be used to give accurate calculations for alpha decay half-lives.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

A collision model for DNA separation by capillary electrophoresis in dilute polymer solution

A theoretical description. based on chemical kinetics and electrochemistry, is given of DNA separation in dilute polymer solution by capillary electrophoresis. A self-consistent model was developed leading to predictions of the DNA electrophoretic velocity as a function of the experimental conditions - polymer concentration, temperature, and electric field strength. The effect of selected experimental variables is discussed. The phenomena discussed are illustrated for the example of 100 bp DNA ladder separation in dilute HPMC solution by capillary electrophoresis. This model is the first single model that can fully explain the dependence of DNA electrophoretic velocity on electrophoretic conditions.