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.     

Rotational analysis ofA 1Π-X 1Σ+ bands of P14N and P15N: Perturbatio studies in theA 1Π state

The ultraviolet band systemsA ¹Π-X ¹Σ⁺ of P¹⁴N and P¹⁵N were excited in an electrodeless tube containing traces of phosphorus specpure nitrogen and neon using a microwave discharge (2450 MHz). Bands of the isotopic species, P¹⁵N, were obtained using¹⁵N₂ enriched to 95.5%. Rotational analyses of eleven bands of P¹⁴N and sixteen bands of P¹⁵N were carried out. Three perturbing statese ³Σ⁻,d ³Δ andb ³Π, arising from the lower valence configurations were identified from the observed perturbations in thev′=0–4 levels of theA ¹Π state. Deperturbation studies led to the determination of molecular constants of the perturbing states. Vibrational assignments of the perturbing states were made from isotope shift studies.     

Structure and Structural Nonrigidity of Tetrachlorodibenzo-p-Dioxin Radical Cations

Ab initio calculations showed that the tetrachlorodibenzo-para-dioxin radical cations (TCDD RCs) with a planar structure have two steady states with asymmetric dioxin cycles. The activation barriers between these states are up to 2 kcal/mole, so that the RCs may be regarded as being structurally nonrigid within the dioxin cycle. The 2,3,7,8-TCDD RC is more stable than the 1,4,6,9-TCDD RC, the energy difference being 5.2 kcal/mole. The adiabatic ionization potential of 2,3,7,8-TCDD (7.54 eV) is 0.1 eV smaller than the corresponding potential of 1,4,6,9-TCDD. These factors account for the increased hemoproteide affinity and hence increased biological activity of 2,3,7,8-TCDD.     

An Application of Branch and Cut to Open Pit Mine Scheduling

The economic viability of the modern day mine is highly dependent upon careful planning and management. Declining trends in average ore grades, increasing mining costs and environmental considerations will ensure that this situation will remain in the foreseeable future. The operation and management of a large open pit mine having a life of several years is an enormous and complex task. Though a number of optimization techniques have been successfully applied to resolve some important problems, the problem of determining an optimal production schedule over the life of the deposit is still very much unresolved. In this paper we will critically examine the techniques that are being used in the mining industry for production scheduling indicating their limitations. In addition, we present a mixed integer linear programming model for the scheduling problems along with a Branch and Cut solution strategy. Computational results for practical sized problems are discussed.     

二维椭圆型方程反问题中优化算法的比较

近年来，决定椭圆型方程系数反问题在地磁、地球物理、冶金和生物等实际问题上有着广泛的应用．本文讨论了二维的决定椭圆型方程系数反问题的数值求解方法．由误差平方和最小原则，这个反问题可化为一个变分问题，并进一步离散化为一个最优化问题，其目标函数依赖于要决定的方程系数．本文着重考察非线性共轭梯度法在此最优化问题数值计算中的表现，并与拟牛顿法作为对比．为了提高算法的效率我们适当选择加快收敛速度的预处理矩阵．同时还考察了线搜索方法的不同对优化算法的影响．数值实验的结果表明，非线性共轭梯度法在这类大规模优化问题中相对于拟牛顿法更有效．     

Computations of two passing‐by high‐speed trains by a relaxation overset‐grid algorithm

This paper presents a relaxation algorithm, which is based on the overset grid technology, an unsteady three‐dimensional Navier–Stokes flow solver, and an inner‐ and outer‐relaxation method, for simulation of the unsteady flows of moving high‐speed trains. The flow solutions on the overlapped grids can be accurately updated by introducing a grid tracking technique and the inner‐ and outer‐relaxation method. To evaluate the capability and solution accuracy of the present algorithm, the computational static pressure distribution of a single stationary TGV high‐speed train inside a long tunnel is investigated numerically, and is compared with the experimental data from low‐speed wind tunnel test. Further, the unsteady flows of two TGV high‐speed trains passing by each other inside a long tunnel and at the tunnel entrance are simulated. A series of time histories of pressure distributions and aerodynamic loads acting on the train and tunnel surfaces are depicted for detailed discussions. Copyright © 2004 John Wiley & Sons, Ltd.     

Photorefractive Four-Wave Mixing Switches Utilizing Pockels Effect -Longitudinal and Lateral Switches-

We carried out detailed calculations for photorefractive wave-mixing switches based on one of three crystals with high electro-optic coefficients, namely, BaTiO₃, Strontium Barium Niobate (SBN (0.75)), and Potasium Sodium Strontium Barium Niobate (KNSBN). A comparison of results for the three crystals shows that a 0_-cut BaTiO³ crystal is suitable for a longitudinal switch and requires a voltage of about 80 for a 2-mm-thick crystal to induce sufficient phase mismatch. The electrodes must be transparent for the incident and diffracted beams. A 45_-cut SBN (0.75) crystal, however, is suitable for a lateral switch and requires a voltage of about 150 for a 1-mm-wide crystal. The electrodes do not need to be transparent.     

Measurements of Time-Varying Characteristics of Optical Signals and Inhomogeneities of Tested Substances with Additional Phase Modulation

The problem of investigation of the amplitude and phase structure of a time-varying probing optical signal and the structure of time-varying inhomogeneities of a substance tested by this signal is considered. The analysis is concerned, in particular, with determination of the structure of signals and processes with resolution in the pico- and femtosecond range. The scheme used for the analysis is based on registration of four spatially separated spectra of the studied radiation. The spectra are formed in a four-channel scheme with a twin-wave Michelson interferometer and a spectral device. Modulators based on electrooptical crystals (perovskites) are placed in the channels. The sum spectra are formed: without modulators, with the effect of either of the modulators, and with both of them affecting the radiation. The effect of the studied substance implies either modulating the radiation (in this case it is described by multiplication) or redistributing the radiation (then it is described by convolution).     

Synthesis,characterization, and biological activity of [2‐oxo‐2‐(4‐acetyl) phenyl amino] ethylene methacrylate and its derivatives

A new type of methacrylate monomer, [2‐oxo‐2‐(4‐acetyl) phenyl amino] ethylene methacrylate (APEMA), was synthesized. The oxime, 2,4‐dinitrophenylhydrazone, and thiosemicarbazone derivatives of poly{[2‐oxo‐2‐(4‐acetyl) phenyl amino] ethylene methacrylate} [poly(APEMA)] were prepared with hydroxylamine hydrochloride, 2,4‐dinitrophenylhydrazine, and thiosemicarbazone hydrochloride, respectively. The radical homopolymerization of APEMA was performed at 65 °C in a 1,4‐dioxane solution with benzoyl peroxide as an initiator. The monomer and its homopolymer were characterized with Fourier transform infrared and NMR techniques. The thermal stabilities of poly(APEMA) and its derivatives were investigated with thermogravimetric analysis and differential scanning calorimetry. The ultraviolet stability of the polymers were compared. The solubility and inherent viscosity of the polymers were also determined. The number‐average and weight‐average molecular weights and polydispersity index of the polymers were determined with gel permeation chromatography. The antibacterial and antifungal effects of the monomer and the polymer and its derivatives were also investigated on various bacteria and fungi. The activation energies of the thermal degradation of the polymers were calculated with the Ozawa method. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 3157–3169, 2004