Golden-rule treatment on the ClO/ClO+ electron-transfer system

The structures, properties and electron transfer reactivity of the ClO/ClO+ coupling system are studied in this paper at ab initio (HF and MP2) levels and the density functional theory (DFT: B3LYP, B3P86, B3PW91) levels employing 6311 + G(3df) basis set and on the basis of the golden-rule of the time-dependent perturbation theory. Investigations indicate that the results got from the B3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated by using the B3LYP/6-311 + G(3df) method, and then the electron transfer rates are determined at this level. The electronic coupling matrix element of EC.6 is very small, only 0.03 kcal/mol, while that of EC.7 is the biggest, being 12.41 kcal/mol, the corresponding electron transfer rate is also the fastest among these seven encounter complexes. The averaged electron transfer rate is about 1.672 X 10(11) M-1 s(-1). It is indicated that the structures optimized by B3LYP method are more reliable than the results got from the other four methods. It also testified that the electronic coupling matrix element is the vital factor that significantly affects the electron transfer rate. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical study on the ClO/ClO- system electron-transfer reactivity by the golden-rule

The structures, properties and electron transfer reactivity of the ClO/ClO- coupling system are studied in this paper at ab initio (UHF and UMP2) levels and the Density Functional Theory (DFT: UB3LYP, UB3P86, UB3PW91) levels employing 6311 + G(3df) basis set and on the basis of the Golden-rule of the time-dependent perturbation theory. Investigations indicate that the results obtained using the UB3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. For this coupling system, six stable coupling modes have been found which correspond to six different encounter complexes and denote six different electron transfer mechanism: four O-O directly linked structures (one collinear: D-h, one anti-parallel: C-s, two twist: C-2) and two Cl-O linked structures (cis- and anti- C-s structures). The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated for the electron transfer reactions via these six different mechanism at the UB3LYP/6-311 + G(3df) level, and then the electron transfer rates are determined at the same level. The most favorable coupling mode to the electron transfer is the anti-parallel mechanism. The averaged electron transfer rate is about 5.58 X 10(11) M-1 s(-1). It is also implied that the B3LYP method can give more reasonable results for the electron transfer reactivity of this system. (C) 2003 Elsevier B.V. All rights reserved.

Collisional quantum interference effect on rotational energy transfer in an atom-diatom system

The theoretical model of collisional quantum interference (CQI) in intramolecular rotational energy transfer is described in an atom-diatom system, based on the first Born approximation of time-dependent perturbation theory and considering a long-range interaction potential. The relation between differential and integral interference angles is obtained. For the CO A(1)Pi (v = 0)/e(3)Sigma (-)(v = 1)-He collision system, the calculated integral interference angles are consistent with the experimental values. The physical significance of interference angle and the essential factors it depends on as well as the influence of the short-range interaction on CQI are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Multipolar model for collisional quantum interference on rotational energy transfer

A theoretical model of collisional quantum interference (CQI) is developed in a diatom-diatom system based on the first-order Born approximation of time-dependent perturbation theory and the multipolar interaction potential. The transition cross section is obtained. The relations between the differential and integral interference angles are discussed. The key factors on the determination of the differential and integral interference angles are obtained. The changing tendency of the interference angles with the experimental temperatures is obtained.

基于扰动方法的广义最佳鉴别矢量集求解的一种迭代算法

本文对广义最佳鉴别矢量的求解方法进行研究 ,根据矩阵的扰动理论 ,改进了郭提出的求解广义最佳鉴别矢量的一种迭代算法 ,提出了求解广义最佳鉴别矢量的一种新的迭代算法 .本文算法的一个突出优点是随着类别数目的增加 ,计算时间反而缩短 ;而老算法随着类别数目的增加计算时间随着增加 ;不仅如此 ,新算法的识别率不劣于老算法 .在 ORL人脸数据库的数值实验 ,验证了上述论断的正确性

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

Energy Conversion in Shape Memory Alloy Heat Engine Part I: Theory

Shape Memory Alloy (SMA) can be easily deformed to a new shape by applying a small external load at low temperature, and then recovers its original configuration upon heating. This unique shape memory phenomenon has inspired many novel designs. SMA based heat engine is one among them. SMA heat engine is an environment-friendly alternative to extract mechanical energy from low-grade energies, for instance, warm wastewater, geothermal energy, solar thermal energy, etc. The aim of this paper is to present an applicable theoretical model for simulation of SMA-based heat engines. First, a micro-mechanical constitutive model is derived for SMAs. The volume fractions of austenite and martensite variants are chosen as internal variables to describe the evolution of microstructure in SMA upon phase transition. Subsequently, the energy equation is derived based on the first thermodynamic law and the previous SMA model. From Fourier’s law of heat conduction and Newton’s law of cooling, both differential and integral forms of energy conversion equation are obtained.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

A general theory on media with randomly distributed inclusions .1. The average field behaviors

In this paper, a theory is developed to calculate the average strain field in the materials with randomly distributed inclusions. Many previous researches investigating the average field behaviors were based upon Mori and Tanaka's idea. Since they were restricted to studying those materials with uniform distributions of inclusions they did not need detailed statistical information of random microstructures, and could use the volume average to replace the ensemble average. To study more general materials with randomly distributed inclusions, the number density function is introduced in formulating the average field equation in this research. Both uniform and nonuniform distributions of inclusions are taken into account in detail.

Kinetic theory for particle concentration distribution in two phase flow

The results of experiments in open channels and closed pipelines show two kinds of patterns for the vertical distribution of particle concentration (i.e., pattern I and pattern II). The former shows a pattern of maximum concentration at some location above the bottom and the downward decay of the concentration below the location. The latter always shows an increase of the particle concentration downward over the whole vertical, with the maximum value at the bottom. Many investigations were made on the pattern II, but few were made on pattern I. In this paper, a particle velocity distribution function is first obtained in the equilibrium state or in dilute steady state for the particle in two-phase flows, then a theoretical model for the particle concentration distribution is derived from the kinetic theory. More attention is paid to the predictions of the concentration distribution of pattern I and comparisons of the present model are made with the data measured by means of laser doppler anemometry (LDA). Very good agreements are obtained between the measured and calculated results.

A closure theory of intermittency of turbulence

The variational approach to the closure problem of turbulence theory, proposed in an earlier article [Phys. Fluids 26, 2098 (1983); 27, 2229 (1984)], is extended to evaluate the flatness factor, which indicates the degree of intermittency of turbulence. Since the flatness factor is related to the fourth moment of a turbulent velocity field, the corresponding higher-order terms in the perturbation solution of the Liouville equation have to be considered. Most closure methods discard these higher-order terms and fail to explain the intermittency phenomenon. The computed flatness factor of the idealized model of infinite isotropic turbulence ranges from 9 to 15 and has the same order of magnitude as the experimental data of real turbulent flows. The intermittency phenomenon does not necessarily negate the Kolmogorov k−5/3 inertial range spectrum. The Kolmogorov k−5/3 law and the high degree of intermittency can coexist as two consistent consequences of the closure theory of turbulence. The Kolmogorov 1941 theory [J. Fluid Mech. 62, 305 (1974)] cannot be disqualified merely because the energy dissipation rate fluctuates.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

Singularity criteria for perturbation series

This paper is concerned with some mathematical aspects of the Van Dyke method inperturbation theory, i.e. the singularity criteria of perturbation series. The author suggestsa sign criterion and a Domb-syke plot for the cases with complex conjugate singularities, thussucceeding in extending the conclusions of Van Dyke's. Subsequently. effects of singularitiesof the lower order upon the criteria are taken into account. In addition, a method of locat-ing singular points is developed by analysing the new perturbation series derived by the Eulertransformation.

On the non-linear stability theory of spinning solid bodies with liquid-filled cavities and its application to the columbus problem

In this paper, we first present a system of differential-integral equations for the largedisturbance to the general case that any arbitrarily shaped solid body with a cavity contain-ing viscous liquid rotates uniformly around the principal axis of inertia, and then develop aweakly non-linear stability theory by the Lyapunov direct approach. Applying this theoryto the Columbus problem, we have proved the consistency between the theory and Kelvin'sexperiments.

Stocking models of Chinese mitten crab (Eriocheir japonica sinensis) in Yangtze lakes

The fanning of Chinese mitten crab, a quality aquatic product in China and neighbouring Asian countries, has been developing rapidly in China since last decade. It reached a total yield of 3.4 X 10(5) tonnes in 2002. Due to the successive over-stocking year after year, many lakes in the mid-lower Yangtze Basin, the main farming area, are under deterioration, leading to a reduction of crab yield and quality, and, subsequently, a loss of fanning profits. Aiming at a normal development of crab culture and the sustainable use of lakes, an annual investigation dealing with lake environmental factors in relation to stocked crab populations was carried out at 20 farms in 4 lakes. The results show that the submersed macrophyte biomass (B-Mac) is the key factor affecting annual crab yield (CY). Using the ratio of Secchi depth to mean depth (Z(SD)/Z(M)), an easily measured parameter closely correlated to BMac, as driving variable, 10 regression models of maximal crab yields were generated (r(2) ranging 0.49-0.81). Based on the theory of MSY (Maximum Sustainable Yield), in combination with body-weight (BW) and recapture rate (RR) of adult crabs, a general optimal stocking model was eventually formulated. All models are simple and easy to operate. Comments on their applications and prospects are given in brief. (c) 2006 Elsevier B.V. All rights reserved.