多点边值问题的插值矩阵法及误差分析

插值矩阵法是求解多点边值问题的数值法。本文给出的该法的误差分析,论证了插值矩阵法解得的y(x),y'(x),…,y^(m)(x)有相同的精度,并对二阶方程,给出该法的稳定性证明和收敛阶。     

块多分裂方法与预条件子空间迭代方法

提出一种块多分裂并行PE迭代算法(MPPE),可以克服M^-1r^(s)并行化处理的困难。这种算法格式简单明了,收敛速度快。并证明了当矩阵A是M-阵和H-阵时,该算法是收敛的。同时把这种分裂作为预处理矩阵,对子空间方法类进行了预处理,并给出的计算实例显示该算法很有效,对子空间方法类的余量光滑和加速都起到了比较好的作用。     

模M之m阶单位根的计数与确定

要根据[1]中给出的MCS序列x_n≡αx_n-1(modM)具有周期m的充要条件,确定其相应全部参数α,可归结为模M之m阶单位根的计数与确定问题,而此问题对一般复合模M至今未见解决。本文解决了这个问题,从而也解决了MCS序列具有周期m的全部参数α的计数和确定问题。     

计算流体中一个改进的强隐式格式及迭代的收敛性

本文分析了Stone引入辅助矩阵Ã将五对角阵A变为七对角阵A+Ã的过程,指出Stone引入的辅助阵Ã不具备对称性,由此会给迭代收敛的分析工作带来困难。在此分析的基础上构造了一个对称的辅助矩阵,当A为对称时可使A+Ã对称正定。给出选取最优迭代参数ω的公式和办法。通过典型算例对这种改进的强隐式格式的数值稳定性作了验证。     

可向量计算的块预条件迭代算法

提出一种类似于PE算法的实用并行迭代算法(VPE),可以克服M^-1r^(s)向量或并行化处理的困难.这种算法格式简单明了,收敛速度快.并证明了当矩阵A是M-阵和H-阵时,该算法是收敛的。计算实例显示该算法很有效.     

计算自旋-s算子幺正演化矩阵ds(t)的新方法及其应用

提出了一种严格求解任意自旋-s算子幺正演化矩阵的方法,该方法不同于群论的方法和直接计算的方法,是一种间接的算法.方法的核心是利用两个系统表示的等价性：即自旋-s算子Hamiltonian量H_s=S_x与Heisenberg XX开链带相互作用J_n=n(N-n)的Hamiltonian量的等价性,由于存在这种等价性,自旋-s算子幺正演化矩阵的计算可通过Heisenberg XX开链中态的演化来实现.采用该方法计算了s=3/2,s=2和s=5/2时对应的幺正演化矩阵.由于初始态|sm〉在算子e^-itS_x下的演化实质上相当于对态|sm〉进行一个绕x轴转角为β＝t的转动,演化矩阵元d^s_m'm(t)=〈sm′|e^-itS_x|sm〉就是转动后的态e^-itS_x|sm〉在|sm′〉态上的投影值,所以在t＝π时刻的演化矩阵刚好对应Heisenberg XX开链上量子态的理想传输. 关键词： s算子')" href="#">自旋-s算子 幺正演化矩阵 量子态传输     

色散方程的一类具任意稳定性的显格式

本文对色散方程u_r=au_xxx构造了一类中间层包含四个结点,带两个参数m和θ的三层显式差分格式。当m和θ满足一定的关系时,其稳定性条件为|γ|≤(m+1)/(4(m-1))(|m|>1),从而当取m充分接近1时,可得到任意大的稳定性条件,并且保持截断误差阶不变。数值例子验证了理论分析的结果。     

Jacobi矩阵特征值反问题

研究如下一类Jacobi矩阵特征值反问题:问题IEP:给定两个互异实数λ,μ(λ<μ)和两个n维非零实向量x,y,求n阶Jacobi矩阵J,使得(λ,x),(μ,y)分别恰是J的第i,j(i≠j)个特征对。还分析了Jacobi矩阵的特征性质,给出了一个特征对恰是Jacobi矩阵J的第i个特征对的充分必要条件,由此导出了问题IEP有解的充分必要条件。     

二维非定常Sine-Gordon方程辛算法及其孤子数值模拟

在矩形域[-a,a]×[-a,a]内对微分算子L=(ə²)/(əx²)+(ə²)/(əy²)用5点差分格式将二维非定常Sine Gordon方程离散化为一个2×799²阶非线性Hamilton系统.对该系统使用Euler中心格式,得到一个非线性方程组.对此方程组建立迭代解法并给出了这个迭代方法的收敛条件和收敛速度.Sine Gordon方程单孤子和双孤子的数值模拟试验显示该辛算法是有效的.     

波动方程辛算法的迭代求解

对(∂²)/(∂x²)利用中心差商算子,对expt作对角Padé逼近,由波动偏微分方程可得到两类具有O(Δx²+Δt^2l)和O(Δx⁴+Δt^2l)精度的辛格式.对由此类辛格式产生的线性方程组构造了两种迭代解法,并对l=1,2,3,4给出了它们的收敛条件.并进行了数值实验.     

Singularities and accumulation of singularities of πN scattering amplitudes

It is demonstrated that for the isospin I = 1/2 πN scattering amplitude, T^I=1/2(s, t), $s={\left({m}_{N}^{2}-{m}_{\pi }^{2}\right)}^{2}/{m}_{N}^{2}$ and $s={m}_{N}^{2}+2{m}_{\pi }^{2}$ are two accumulation points of poles on the second sheet of complex s plane, and are hence accumulation of singularities of T^I=1/2(s, t). For T^I=3/2(s, t), $s={\left({m}_{N}^{2}-{m}_{\pi }^{2}\right)}^{2}/{m}_{N}^{2}$ is the accumulation point of poles on the second sheet of the complex s plane. The proof is valid up to all orders of chiral expansions.     

使用迭代法的终止检验

为了保证使用迭代法解偏微分方程数值解时的精度,终止检验是一个非常重要的问题。当线性代数方程组的系数矩阵是严格对角优势或不可约弱对角优势时,本文给出了使用迭代法u^s+1=Gu^s+K解偏微分方程时的终止检验公式。并且还给出一些数值例子。     

Analytical result for dimensionally regularized massless master non-planar double box with one-leg off shell

The dimensionally regularized massless non-planar double box Feynman diagram with powers of propagators equal to one, one leg off the mass shell, i.e., with p₁²=q²≠0, and three legs on shell, p_i²=0, i=2,3,4, is analytically calculated for general values of q² and the Mandelstam variables s,t and u (not necessarily restricted by the physical condition s+t+u=q²). An explicit result is expressed through (generalized) polylogarithms, up to the fourth order, dependent on rational combinations of q²,s,t and u, and simple finite two- and three-fold Mellin–Barnes integrals of products of gamma functions which are easily numerically evaluated for arbitrary non-zero values of the arguments.     

A consistent next-to-leading-order QCD calculation of hadronic diffractive scattering

We calculate the order _s² and order _s³ QCD contributions to colour-singlet exchange in the leading log s approximation. We implement the resulting amplitude at the hadronic level and thus construct the QCD pomeron and odderon to this order of perturbation theory. We show that the structure of the hadronic form factors provides a natural mechanism through which the odderon gets suppressed at t = 0 whereas it dominates the elastic cross section at large t. We also demonstrate that the inclusion of nonperturbative effects through a modification of the gluon propagator accelerates greatly the convergence of the log s expansion, although not enough to provide agreement with the data.     

Exact solutions to the angular Teukolsky equation with s ≠ 0

We first convert the angular Teukolsky equation under the special condition of τ ≠ 0, s ≠ 0, m=0 into a confluent Heun differential equation (CHDE) by taking different function transformation and variable substitution. And then according to the characteristics of both CHDE and its analytical solution expressed by a confluent Heun function (CHF), we find two linearly dependent solutions corresponding to the same eigenstate, from which we obtain a precise energy spectrum equation by constructing a Wronskian determinant. After that, we are able to localize the positions of the eigenvalues on the real axis or on the complex plane when τ is a real number, a pure imaginary number, and a complex number, respectively and we notice that the relation between the quantum number l and the spin weight quantum number s satisfies the relation l=∣s∣+ n, n=0, 1, 2···. The exact eigenvalues and the corresponding normalized eigenfunctions given by the CHF are obtained with the aid of Maple. The features of the angular probability distribution (APD) and the linearly dependent characteristics of two eigenfunctions corresponding to the same eigenstate are discussed. We find that for a real number τ, the eigenvalue is a real number and the eigenfunction is a real function, and the eigenfunction system is an orthogonal complete system, and the APD is asymmetric in the northern and southern hemispheres. For a pure imaginary number τ, the eigenvalue is still a real number and the eigenfunction is a complex function, but the APD is symmetric in the northern and southern hemispheres. When τ is a complex number, the eigenvalue is a complex number, the eigenfunction is still a complex function, and the APD in the northern and southern hemispheres is also asymmetric. Finally, an approximate expression of complex eigenvalues is obtained when n is greater than ∣s∣.     

Neutrino masses, mixing and hierarchy

We build a model to describe neutrinos based on strict hierarchy, incorporating as much as possible, the latest known data, for Δ_sol and Δ_atm, and for the mixing angles determined from neutrino oscillation experiments, including that from KamLAND. Since the hierarchy assumption is a statement about mass ratios, it lets us obtain all three neutrino masses. We obtain a mass matrix, M_ν and a mixing matrix, U, where both M_ν and U are given in terms of powers of Λ, the analog of the Cabibbo angle λ in the Wolfenstein representation, and two parameters, ρ and κ, each of order one. The expansion parameter, Λ, is defined by

Full-size image

, and ρ expresses our ignorance of the lightest neutrino mass m₁, (m₁=ρΛ⁴m₃), while κ scales s₁₃ to the experimental upper limit, s₁₃=κΛ²≈0.16κ. These matrices are similar in structure to those for the quark and lepton families, but with Λ about 1.6 times larger than the λ for the quarks and charged leptons. The upper limit for the effective neutrino mass in double β-decay experiments is 4×10⁻³ eV if s₁₃=0 and 6×10⁻³ eV if s₁₃ is maximal. The model, which is fairly unique, given the hierarchy assumption and the data, is compared to supersymmetric extension and texture zero models of mass generation.     

长程表面等离子体的增强效应

研究了双层金属薄膜构型中构型参数对长程表面等离子体的影响,并发现了衰减全反射激发方法下长程表面等离子体的增强效应.以特征矩阵算法为基础,通过数值计算构型的反射谱,研究构型参数的变化对反射谱的影响.发现由于衰减全反射激发方法中耦合器的存在导致的非对称特性,会使双层金属薄膜构型中的长程表面等离子体拥有本征模式特性以外的有趣特性,如长程模式得到增强而另一支受到抑制,从而使能量更为集中在希望被激发的一支.研究结果对非对称激发构型中的长程表面等离子体研究具有启发意义.     

Preferential attachment network model with aging and initial attractiveness

In this paper, we generalize the growing network model with preferential attachment for new links to simultaneously include aging and initial attractiveness of nodes. The network evolves with the addition of a new node per unit time, and each new node has m new links that with probability Π_i are connected to nodes i already present in the network. In our model, the preferential attachment probability Π_i is proportional not only to k_i + A, the sum of the old node i's degree k_i and its initial attractiveness A, but also to the aging factor ${\tau }_{i}^{-\alpha }$, where τ_i is the age of the old node i. That is, ${{\rm{\Pi }}}_{i}\propto ({k}_{i}+A){\tau }_{i}^{-\alpha }$. Based on the continuum approximation, we present a mean-field analysis that predicts the degree dynamics of the network structure. We show that depending on the aging parameter α two different network topologies can emerge. For α < 1, the network exhibits scaling behavior with a power-law degree distribution P(k) ∝ k^−γ for large k where the scaling exponent γ increases with the aging parameter α and is linearly correlated with the ratio A/m. Moreover, the average degree k(t_i, t) at time t for any node i that is added into the network at time t_i scales as $k({t}_{i},t)\propto {t}_{i}^{-\beta }$ where 1/β is a linear function of A/m. For α > 1, such scaling behavior disappears and the degree distribution is exponential.     

On the Riemann–Hilbert problem of a generalized derivative nonlinear Schrödinger equation

In this work,we present a unified transformation method directly by using the inverse scattering method for a generalized derivative nonlinear Schr?dinger(DNLS)equation.By establishing a matrix Riemann-Hilbert problem and reconstructing potential function q(x,t)from eigenfunctions{Gj(x,t,η)}3/1 in the inverse problem,the initial-boundary value problems for the generalized DNLS equation on the half-line are discussed.Moreover,we also obtain that the spectral functions f(η),s(η),F(η),S(η)are not independent of each other,but meet an important global relation.As applications,the generalized DNLS equation can be reduced to the Kaup-Newell equation and Chen-Lee-Liu equation on the half-line.