A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

关于Banach空间中凸泛函的广义次梯度不等式

本文在前人^[1,2]的基础之上，以凸泛函的次梯度不等式为工具，将Jensen不等式推广到Banach空间中的凸泛函，导出了Banach空间中的Bochner积分型的广义Jensen不等式，给出其在Banach空间概率论中某些应用，从而推广了文献[3—6]的工作．     

Desorption behavior of a surfactant in a linear low‐density polyethylene blend at elevated temperatures

The desorption behavior of a surfactant in a linear low‐density polyethylene (LLDPE) blend at elevated temperatures of 50, 70, and 80 °C was studied with Fourier transform infrared spectroscopy. The composition of the LLDPE blend was 70:30 LLDPE/low‐density polyethylene. Three different specimens (II, III, and IV) were prepared with various compositions of a small molecular penetrant, sorbitan palmitate (SPAN‐40), and a migration controller, poly(ethylene acrylic acid) (EAA), in the LLDPE blend. The calculated diffusion coefficient (D) of SPAN‐40 in specimens II, III, and IV, between 50 and 80 °C, varied from 1.74 × 10^?11 to 6.79 × 10^?11 cm²/s, from 1.10 × 10^?11 to 5.75 × 10^?11 cm²/s, and from 0.58 × 10^?11 to 4.75 × 10^?11 cm²/s, respectively. In addition, the calculated activation energies (E_D) of specimens II, III, and IV, from the plotting of ln D versus 1/T between 50 and 80 °C, were 42.9, 52.7, and 65.6 kJ/mol, respectively. These values were different from those obtained between 25 and 50 °C and were believed to have been influenced by the interference of Tinuvin (a UV stabilizer) at elevated temperatures higher than 50 °C. Although the desorption rate of SPAN‐40 increased with the temperature and decreased with the EAA content, the observed spectral behavior did not depend on the temperature and time. For all specimens stored over 50 °C, the peak at 1739 cm^?1 decreased in a few days and subsequently increased with a peak shift toward 1730 cm^?1. This arose from the carbonyl stretching vibration of Tinuvin, possibly because of oxidation or degradation at elevated temperatures. In addition, the incorporation of EAA into the LLDPE blend suppressed the desorption rate of SPAN‐40 and retarded the appearance of the 1730 cm^?1 peak. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1114–1126, 2004     

X波段六腔渡越管振荡器的高频特性研究

 从麦克斯韦方程出发，采用时域有限差分法（FDTD）和离散傅里叶变换（DFFT）相结合的方法，通过数值计算得出了六腔开放式谐振腔中前四个谐振频率和场分布，计算出的谐振频率与实验测量结果基本相同。比较了开放腔和封闭腔谐振频率，验证了TEM波吸收边界条件，并在实际编程计算中得以应用。计算结果为六腔渡越管振荡器的机理研究提供了依据。     

非均匀网格时域伪谱算法在超宽带技术中的应用

 与传统时域有限差分算法相比，采用以伪谱方法离散Maxwell微分方程为基础的时域伪谱(PSTD)算法计算大的电尺度电磁场时域问题，将大大提高计算效率，降低内存需求。为了拓宽PSTD算法的应用，近年来，基于网格插值方法的非均匀时域伪谱算法得到了发展。研究的重点是算法中非均匀网格技术的实现及其在时域瞬态脉冲电磁场模拟和高功率超宽带脉冲技术方面的应用。以高斯脉冲为激励源，用该算法计算了多层介质的反射和透射，并通过超宽带脉冲穿墙实验对这一方法的应用进行了验证。模拟和实验结果具有较好的一致性。     

Factorization of operator valued Lp for 0 ≤ p < 1

In [5], it is proved that a bounded linear operator u, from a Banach space Y into an L_p(S, ν) factors through L_p1 (S, ν) for some p₁ > 1, if Y* is of finite cotype; (S, ν) is a probability space for p = 0, and any measure space for 0 < p < 1. In this paper, we generalize this result to uv, where u : Y → L_p(S, ν) and v : X → Y are linear operators such that v* is of finite Ka?in cotype. This result gives also a new proof of Grothendieck's theorem. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Moments of an exponential functional of random walks and permutations with given descent sets

The exponential functional of simple, symmetric random walks with negative drift is an infinite polynomial Y = 1 + ξ₁ + ξ₁ξ₂ + ξ₁ξ₂ξ₃ + ⋯ of independent and identically distributed non-negative random variables. It has moments that are rational functions of the variables μ _k = E(ξ ^k ) < 1 with universal coefficients. It turns out that such a coefficient is equal to the number of permutations with descent set defined by the multiindex of the coefficient. A recursion enumerates all numbers of permutations with given descent sets in the form of a Pascal-type triangle. This revised version was published online in August 2006 with corrections to the Cover Date.     

A Characterisation of Posets that Are nearly Antichains

We present a characterisation of posets of size n with linear discrepancy n − 2. These are the posets that are “furthest” from a linear order without being an antichain. This revised version was published online in September 2006 with corrections to the Cover Date.     

Deciding confluence of certain term rewriting systems in polynomial time

We present a characterization of confluence for term rewriting systems, which is then refined for special classes of rewriting systems. The refined characterization is used to obtain a polynomial time algorithm for deciding the confluence of ground term rewrite systems. The same approach also shows the decidability of confluence for shallow and linear term rewriting systems. The decision procedure has a polynomial time complexity under the assumption that the maximum arity of a function symbol in the signature is a constant.     

Boundary position feedback control of Kirchhoff's non‐linear strings

This paper is concerned with global stabilization of an undamped non‐linear string in the case where any velocity feedback is not available. The linearized system has an infinite number of poles and zeros on the imaginary axis. In the case where any velocity feedback is not available, a parallel compensator is effective. The stabilizer is constructed for the augmented system which consists of the controlled system and a parallel compensator. It is proved that the string can be stabilized by linear boundary control. Copyright © 2003 John Wiley & Sons, Ltd.