Random Copolymers of 1,5-Dioxepan-2-one

ABSTRACT

Copolymers of 1,5-dioxepan-2-one (DXO) and e-caprolactone (?-CL), δ-valerolactone (δ-VL) or L-lactide (LLA) have been synthesized and characterized. High molecular weight copolymers were obtained using stannous-2-ethyl hexanoate as catalyst in bulk. Reactivity ratios for the copolymerization of DXO and δ-VL were determined at 110°C as r_VL=0.5 and r_DXO=2.3. At high conversion, depolymerization of δ-VL occurred, resulting in lower molecular weight and variations in the copolymer composition.

Physical properties, such as crystallinity and melting temperature of the DXO-copolymers proved to be strongly dependent on the choice of comonomer and on the molar composition of the copolymers. DXO appears to be incorporated into the poly-?-caprolactone (PCL) crystals and to some extent into the poly-δ-valerolactone (PVL) crystals, resulting in a more gradual decrease in crystallinity with increasing amount of DXO.     

Thirteen ways to estimate global error

Summary Various techniques that have been proposed for estimating the accumulated discretization error in the numerical solution of differential equations, particularly ordinary differential equations, are classified, described, and compared. For most of the schemes either an outline of an error analysis is given which explains why the scheme works or a weakness of the scheme is illustrated.This research is partially supported by NSF Grant No. MCS-8107046     

Fusion around the barrier for 7Li+12C

Fusion cross-sections for the ⁷Li + ¹²C reaction have been measured at energies above the Coulomb barrier by the direct detection of evaporation residues. The heavy evaporation residues with energies below 3 MeV could not be separated out from the α-particles in the spectrum and hence their contribution was estimated using statistical model calculations. The present work indicates that suppression of fusion cross-sections due to the breakup of ⁷Li may not be significant for ⁷Li + ¹²C reaction at energies around the barrier.     

Symplectic integrators and the conservation of angular momentum

In this article we observe that generally symplectic integrators conserve angular momentum exactly, whereas nonsymplectic integrators do not. We show that this observation extends to multiple timesteps and to constrained dynamics. Both of these devices are important for efficient molecular dynamics simulations. © 1995 by John Wiley & Sons, Inc.     

Limits of parallelism in explicit ODE methods

Numerical methods for ordinary initial value problems that do not depend on special properties of the system are usually found in the class of linear multistage multivalue methods, first formulated by J.C. Butcher. Among these the explicit methods are easiest to implement. For these reasons there has been considerable research activity devoted to generating methods of this class which utilize independent function evaluations that can be performed in parallel. Each such group of concurrent function evaluations can be regarded as a stage of the method. However, it turns out that parallelism affords only limited opportunity for reducing the computing time with such methods. This is most evident for the simple linear homogeneous constant-coefficient test problem, whose solution is essentially a matter of approximating the exponential by an algebraic function. For a given number of stages and a given number of saved values, parallelism offers a somewhat enlarged set of algebraic functions from which to choose. However, there is absolutely no benefit in having the degree of parallelism (number of processors) exceed the number of saved values of the method. Thus, in particular, parallel one-step methods offer no speedup over serial one-step methods for the standard linear test problem. Although the implication of this result for general nonlinear problems is unclear, there are indications that dramatic speedups are not possible in general. Also given are some results relevant to the construction of methods.Work supported in part by National Science Foundation grants DMS 89 11410 and DMS 90 15533 and by US Department of Energy grant DOE DEFG02-87ER25026. Work of the second author was completed while at the University of Illinois.     

Fast evaluation of polarizable forces

Polarizability is considered to be the single most significant development in the next generation of force fields for biomolecular simulations. However, the self-consistent computation of induced atomic dipoles in a polarizable force field is expensive due to the cost of solving a large dense linear system at each step of a simulation. This article introduces methods that reduce the cost of computing the electrostatic energy and force of a polarizable model from about 7.5 times the cost of computing those of a nonpolarizable model to less than twice the cost. This is probably sufficient for the routine use of polarizable forces in biomolecular simulations. The reduction in computing time is achieved by an efficient implementation of the particle-mesh Ewald method, an accurate and robust predictor based on least-squares fitting, and non-stationary iterative methods whose fast convergence is accelerated by a simple preconditioner. Furthermore, with these methods, the self-consistent approach with a larger timestep is shown to be faster than the extended Lagrangian approach. The use of dipole moments from previous timesteps to calculate an accurate initial guess for iterative methods leads to an energy drift, which can be made acceptably small. The use of a zero initial guess does not lead to perceptible energy drift if a reasonably strict convergence criterion for the iteration is imposed.     

Secondary ion yields produced by keV atomic and polyatomic ion impacts on a self-assembled monolayer surface

A suite of keV polyatomic or 'cluster' projectiles was used to bombard unoxidized and oxidized self-assembled monolayer surfaces. Negative secondary ion yields, collected at the limit of single ion impacts, were measured and compared for both molecular and fragment ions. In contrast to targets that are orders of magnitude thicker than the penetration range of the primary ions, secondary ion yields from polyatomic projectile impacts on self-assembled monolayers show little to no enhancement when compared with monatomic projectiles at the same velocity. This unusual trend is most likely due to the structural arrangement and bonding characteristics of the monolayer molecules with the Au(111). Copyright 1999 John Wiley & Sons, Ltd.     

Variable step size destabilizes the Störmer/leapfrog/verlet method

For fixed step-sizeh the Störmer method is stable for the standard test equationÿ= ² y,>0, if and only ifh<2. We show that for variable step sizeh _n there does not exist a (positive) limit onh which ensures stability. Nor can we guarantee stability if, in addition, we limit the step size ratioh _n/h _n–1.This work was supported in part by National Science Foundation Grant DMS 90 15533.     

Multiple grid methods for classical molecular dynamics

Presented in the context of classical molecular mechanics and dynamics are multilevel summation methods for the fast calculation of energies/forces for pairwise interactions, which are based on the hierarchical interpolation of interaction potentials on multiple grids. The concepts and details underlying multigrid interpolation are described. For integration of molecular dynamics the use of different time steps for different interactions allows longer time steps for many of the interactions, and this can be combined with multiple grids in space. Comparison is made to the fast multipole method, and evidence is presented suggesting that for molecular simulations multigrid methods may be superior to the fast multipole method and other tree methods.     

Corrected potential energy functions for constrained molecular dynamics

Atomic oscillations present in classical molecular dynamics restrict the step size that can be used. Multiple time stepping schemes offer only modest improvements, and implicit integrators are costly and inaccurate. The best approach may be to actually remove the highest frequency oscillations by constraining bond lengths and bond angles, thus permitting perhaps a 4-fold increase in the step size. However, omitting degrees of freedom produces errors in statistical averages, and rigid angles do not bend for strong excluded volume forces. These difficulties can be addressed by an enhanced treatment of holonomic constrained dynamics using ideas from papers of Fixman (1974) and Reich (1995, 1999). In particular, the 1995 paper proposes the use of “flexible” constraints, and the 1999 paper uses a modified potential energy function with rigid constraints to emulate flexible constraints. Presented here is a more direct and rigorous derivation of the latter approach, together with justification for the use of constraints in molecular modeling. With rigor comes limitations, so practical compromises are proposed: simplifications of the equations and their judicious application when assumptions are violated. Included are suggestions for new approaches.