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.     

Synthesis and characterization of oxidized cellulose

Samples of oxidized cellulose (OC) with various carboxyl contents and degrees of crystallinity were obtained by the oxidation of native and mercerized cellulose with a solution of nitrogen(IV) oxide in CCl₄. A detailed characterization of these OC samples was performed. The effect of oxidation conditions (concentration of N₂O₄ in the solution and oxidation time) and starting cellulose material on OC characteristics (carboxyl, carbonyl and nitrogen content, degree of crystallinity and polymerization, surface area and swelling, and acidic properties) was investigated. Reactivity in the oxidation process was higher in mercerized cellulose than in native cellulose. The action of dilute solutions (10–15%) of N₂O₄ did not affect the degree of crystallinity of cellulose samples. Under these conditions, the oxidation took place mainly in amorphous regions and on the surface of crystallites. Oxidation in a concentrated (40%) N₂O₄ solution led to the destruction of crystallites, which increased the surface area and swelling of cellulose in water. The surface area and the swelling of OC samples increased with a decrease in the index of crystallinity. The acidic properties of OC were shown to increase with an increase of swelling in water. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4785–4791, 2004     

A note on the necessary conditions for the algebraic Riccati equation

We derive some useful and easily computable necessary conditionsfor the existence of a positive semi-definite solution to thealgebraic Ricatti equation (ARE). A motivating example is givento highlight the usefulness of the conditions for controllerand observer designs for nonlinear systems. Further, an upperbound on the trace of the solution to the ARE is also derived.     

On weak spreadability of distributed-parameter systems and its achievement via linear--quadratic control techniques

In this paper, the concept of weak spreadability is introduced.It constitutes an extension of the idea developed by El Jai& Kassara. In the case of linear distributed systems weconsider quadratic control techniques with a conveniently penalizedcriterion which makes the system weakly spreadable. The approachis outlined for a convection—diffusion system, and theresults of a numerical study are also included to illustratethe main features of the considered problem.     

Aqueous solution and solid state interactions of lanthanide ions with a methacrylic ester polymer bearing pendant 15‐crown‐5 moieties

The interaction between trivalent lanthanide ions and poly(1,4,7,10,13‐pentaoxacyclopentadecan‐2‐yl‐methyl methacrylate), PCR5, in aqueous solution and in the solid state have been studied. In aqueous solution, evidence of a weak interaction between the lanthanides and PCR5 comes from the small red shift of the Ce(III) emission spectra and the slight broadening of the Gd(III) EPR spectra. From the Tb(III) lifetimes in the presence of H₂O and D₂O the loss of one or two water coordinated molecules is confirmed when Tb(III) is bound to PCR5. An association constant of the order of 200 M^?1 was obtained for a 1:1 (lanthanide:15‐crown‐5) complex from the shift of the polymer NMR signals induced by Tb(III). A similar association constant is obtained from the differences of the molar conductivity of Ce(III) solution at various concentrations in presence and absence of PCR5. When Tb(III) is adsorbed on PCR5 membranes, lifetime experiments in H₂O and D₂O confirm the loss of 5 or 6 water coordinated molecules indicating that in solid state the lanthanide(III)‐PCR5 interaction is stronger than in solution. The adsorption of Ce(III) in PCR5 membranes shows a Langmuir type isotherm, from which an equilibrium constant of 39 M^?1 has been calculated. SEM shows that the membrane morphology is not much affected by lanthanide adsorption. Support for lanthanide ion–crown interactions comes from ab initio calculations on 15‐crown‐5/La(III) complex. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 1788–1799, 2007     

A microspectroscopic study of the electronic homogeneity of ordered and disordered Sr2FeMoO6

Disordered Sr₂FeMoO₆ shows a drastic reduction in saturation magnetization compared to highly ordered samples, moreover magnetization as a function of the temperature for different disordered samples shows qualitatively different behaviours. We investigate the origin of such diversity by performing spatially resolved photoemission spectroscopy on various disordered samples. Our results establish that e_xtensive electronic inhomogeneity, arising most probably from an underlying chemical inhomogeneity in disordered samples, is responsible for the observed magnetic inhomogeneity. It is further pointed out that these inhomogeneities are connected with composition fluctuations of the type Sr₂Fe_1+x Mo_1-x O₆ with Fe-rich (x > 0) and Mo-rich (x < 0) regions. Dedicated to Prof J Gopalakrishnan on his 62nd birthday.     

Differential–algebraic equations in multibody system modeling

Numerical methods for the efficient integration of both stiff and nonstiff equations of motion of multibody systems having the form of differential-algebraic equations (DAE) of index 3 are discussed. Linear multi-step ABM and BDF methods are considered for the non-iterational integration of nonstiff DAE. The Park method is proposed for integration of stiff equations. This revised version was published online in June 2006 with corrections to the Cover Date.     

Influence of concentration on the activation energy for diffusion in polymer-solvent systems

The self-diffusion of benzene, toluene, and ethylbenzene in polystyrene have been analyzed using the Vrentas/Duda free-volume diffusion model. Diffusion coefficient predictions suggest an exponential concentration dependence of the activation energy required to overcome attractive forces, E. Without the use of any diffusion data approximating E as zero over the entire concentration range yields self-diffusion coefficient predictions which are in good agreement with experimental data. © 1992 John Wiley & Sons, Inc.     

A note on an SOR-like method for augmented systems

Golub et al. (2001, BIT, 41, 71–85) gave a generalizedsuccessive over-relaxation method for the augmented systems.In this paper, the connection between the SOR-like method andthe preconditioned conjugate gradient (PCG) method for the augmentedsystems is investigated. It is shown that the PCG method isat least as accurate (fast) as the SOR-like method. Numericalexamples demonstrate that the PCG method is much faster thanthe SOR-like method.