广义非简谐振子的多尺度微扰理论

应用多尺度微扰理论到广义非简谐振子, 得到了一阶经典和量子微扰解. 特别是 我们的量子解在极限条件下能方便地转变为经典解, 并且坐标和动量算符的对易 关系的简化十分自然. 与Taylor级数解相比较, 无论是在经典还是在量子解 中频率移动都出现在各阶振动表达式中, 所以多尺度微扰解是弱耦合非简谐振动的较好解法.     

Integral representations of nonnegative solutions for parabolic equations and elliptic Martin boundaries

We consider nonnegative solutions of a parabolic equation in a cylinder D×(0,T), where D is a noncompact domain of a Riemannian manifold. Under the assumption [IU] (i.e., the associated heat kernel is intrinsically ultracontractive), we establish an integral representation theorem: any nonnegative solution is represented uniquely by an integral on (D×{0})∪(M_∂D×[0,T)), where M_∂D is the Martin boundary of D for the associated elliptic operator. We apply it in a unified way to several concrete examples to explicitly represent nonnegative solutions. We also show that [IU] implies the condition [SP] (i.e., the constant function 1 is a small perturbation of the elliptic operator on D).     

Maximum-norm error analysis of a non-monotone FEM for a singularly perturbed reaction-diffusion problem

A non-monotone FEM discretization of a singularly perturbed one-dimensional reaction-diffusion problem whose solution exhibits strong layers is considered. The method is shown to be maximum-norm stable although it is not inverse monotone. Both a priori and a posteriori error bounds in the maximum norm are derived. The a priori result can be used to deduce uniform convergence of various layer-adapted meshes proposed in the literature. Numerical experiments complement the theoretical results. AMS subject classification (2000)  65L10, 65L50, 65L60     

圆柱型正交各向异性圆形薄板的非线性非对称弯曲问题（Ⅱ）

本文利用“正则摄动法”［２］研究文献［１］中由“两变量法”［３］所得到的关于Ｗｎｍ，和的递推方程和递推边界条件的求解问题。求得了本问题的一致有效渐近解，最后作为实例，我们利用“混合摄动法”［４］研究了轴对称线性弯曲问题，并将所得结果与文献［５］中给出的精确解相比较，两者基本上是一致的。     

两种群竞争模型的奇摄动群

本文应用多重尺度法构造出非线性微分方程组的解的渐近展开式。并用微分不等式的技巧，证明原问题的解的存在性，且给出解的一致有效渐近估计．     

Interval parameter perturbation method for evaluating the bounds on natural frequencies of structures with interval parameters

For structural parameters with uncertainties, interval mathematics can, in the case where the probabilistic distribution density of uncertain variables is unavailable, deal with the influence of uncertainties in structural parameters on the response of structures. In order to evaluate the region containing natural frequencies of structures with interval parameters, the interval parameter perturbation method is presented in this paper. The advantage of the present method is its computational efficiency in evaluating the region containing natural frequencies. A numerical example is used to illustrate the efficiency of the method proposed. The project is supported by National Youth Natural Science Foundation of China and National Outstanding Youth Science Foundation of China.     

Application of symplectic algorithms to QCT calculation:H+H_2 system

The basic theory of symplectic algorithm was introduced. A comparison between Runge-Kutta method and symplectic integration method was preformed in the simulation of the long time behavior of H + H2 system on BKMP potential energy surface. Our results reveal a dis-sipative behavior in the integral of ordinary differential equation by the fourth order Runge-Kutta method, which causes incorrect simulation results in QCT calculations. However, when the symplectic integration method is applied, the dissipative behavior is not found in the same system. When the initial state is the same, the energy deviation of fourth order symplectic integral method is almost one percent of that of fourth order Runge-Kutta method in a 60000-step simulation, and that of sixth order symplectic integral method is much less. These results show that the symplectic integral methods are always the better choice in the integral calculation of the long time behavior in maintaining energy conservation.     

Free energy perturbation approach to the critical assessment of selective cyclooxygenase-2 inhibitors

The discovery of selective cyclooxygenase-2 (COX-2) inhibitors represents a major achievement of the efforts over the past few decades to develop therapeutic treatments for inflammation. To gain insights into designing new COX-2-selective inhibitors, we address the energetic and structural basis for the selective inhibition of COX isozymes by means of a combined computational protocol involving docking experiment, force field design for the heme prothetic group, and free energy perturbation (FEP) simulation. We consider both COX-2- and COX-1-selective inhibitors taking the V523I mutant of COX-2 to be a relevant structural model for COX-1 as confirmed by a variety of experimental and theoretical evidences. For all COX-2-selective inhibitors under consideration, we find that free energies of binding become less favorable as the receptor changes from COX-2 to COX-1, due to the weakening and/or loss of hydrogen bond and hydrophobic interactions that stabilize the inhibitors in the COX-2 active site. On the other hand, COX-1-selective oxicam inhibitors gain extra stabilization energy with the change of residue 523 from valine to isoleucine because of the formations of new hydrogen bonds in the enzyme-inhibitor complexes. The utility of the combined computational approach, as a valuable tool for in silico screening of COX-2-selective inhibitors, is further exemplified by identifying the physicochemical origins of the enantiospecific selective inhibition of COX-2 by -substituted indomethacin ethanolamide inhibitors.     

Exponentially fitted symplectic methods for the numerical integration of the Schrödinger equation

The computation of the energy eigenvalues of the one-dimensional time-independent Schrödinger equation is considered. Exponentially fitted and trigonometrically fitted symplectic integrators are obtained, by modification of the first and second order Yoshida symplectic methods. Numerical results are obtained for the one-dimensional harmonic oscillator and Morse potential.AMS subject classification: 65L15Funding by research project 71239 of Prefecture of Western Macedonia and the E.U. is gratefully acknowledged.     

On the relationship betweenT 1 andT 2 for stochastic relaxation models

We consider the relaxation dynamics of two quantum levels coupled to a stochastic bath. We emphasize that even if the matrix elements of the fluctuating Hamiltonian are Gaussian, a second-order cumulant truncation is not exact. For various stochastic models, including the case of a spin-1/2 particle in a fluctuating magnetic field, we calculate 1/T ₁, the population relaxation rate, and 1/T ₂, the phase relaxation rate, up to fourth order in perturbation theory. We show that unlike the commonly accepted second-order result that 1/T ₂1/2T ₁, when fourth-order terms are included, in some instances 1/T ₂<1/2T ₁.