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.     

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.     

Anisotropic friction model of DNA electrophoresis in polymer solutions: Comparison with direct observations

The electrophoresis of DNA chains in uncrosslinked polymer solutions with a Brownian dynamics simulation with an anisotropic friction tensor was analyzed. According to the degree of anisotropy, three types of migration behavior are obtained: fluctuation without or with periodicity between U‐shaped and compact conformations, or migration with linear conformation. We found good agreement between our simulation results and the direct observations of DNA by fluorescence microscopy. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 1316–1322, 2003     

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.     

Experimental investigation and theoretical modeling of liquid entrainment in a horizontal tee with a vertical-up branch

This article describes a comprehensive literature review of liquid entrainment in horizontal pipes with vertical-up branches. Deficiencies in the available data and correlations were identified. The Air–water Test Loop for Advanced Thermal–hydraulic Studies (ATLATS) was constructed and entrainment onset and steady-state data were collected for a wide range of flow conditions. Using new insights gained from experimental testing, the authors developed a model for predicting the onset of entrainment and steady-state entrainment rate. Previously published correlations, along with the new model, are compared against all available data. The new model shows very good agreement with the onset data, but is not very good at predicting branch quality at high liquid flow rates.     

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.     

Light Diffraction of Aligned Polymer Fibers Periodically Dispersed by Phase Separation of Liquid Crystal and Polymer

We have confirmed light diffraction of aligned polymer fibers obtained by a phase separation of an anisotropic-phase solution of liquid crystal and polymer. He—Ne laser light passing through the polymer fibers was scattered in the axis vertical to the fibers, and had two peaks of light intensity symmetrical to the center of the transmitting laser spot. The two peaks were found to be caused by light diffraction due to the periodic polymer-fiber dispersion because the peaks corresponded to values calculated by intervals between the fibers. The periodical fiber networks are considered to be formed by anisotropic spinodal decomposition. This effect can be used to measure the dispersion order of the polymer fibers. © 2004 The Optical Society of Japan     

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