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.     

Fractional Derivatives Applied to Phase-Space Reconstructions

The concept and application of phase-space reconstructions are reviewed. Fractional derivatives are then proposed for the purpose of reconstructing dynamics from a single observed time history. A procedure is presented in which the fractional derivatives of time series data are obtained in the frequency domain. The method is applied to the Lorenz system. The ability of the method to unfold the data is assessed by the method of global false nearest neighbors. The reconstructed data is used to compute recurrences and correlation dimensions. The reconstruction is compared to the commonly used method of delays in order to assess the choice of reconstruction parameters, and also the quality of results.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Tetrafunctional initiators for cationic polymerization of olefins

3,3′,5,5′‐Tetrakis(2‐chloro‐2‐propyl)biphenyl (biphenyl tetracumyl chloride, BPTCC) and 1,3‐bis[3,5‐bis(2‐chloro‐2‐propyl)phenoxy]propane (diphenoxypropane tetracumyl chloride, DPPTCC) were synthesized as initiators for quasiliving cationic polymerization of isobutylene (IB). In the synthesis of BPTCC, tetrafunctionality was achieved via the coupling of dimethyl 5‐bromoisophthalate (DMBI) using nickel dibromide bis(triphenylphosphine) and zinc in the presence of a base; in the synthesis of DPPTCC, two equivalents of dimethyl 5‐hydroxyisophthalate were linked via reaction with 1,3‐dibromopropane in the presence of potassium carbonate. Both initiators were used to initiate the polymerization of IB under quasiliving cationic polymerization conditions. PIB initiated from BPTCC revealed a chain end/molecule value (as determined by ¹H‐NMR) of 3.85, verifying the nearly exclusive production of 4‐arm polyisobutylene (PIB). GPC analysis revealed a narrow peak representing the target four‐arm PIB, with a slight shoulder at high elution volumes (low molecular weights). GPC analysis of the PIB initiated by DPPTCC revealed multimodal distributions, suggesting the formation of two‐, three‐, and four‐arm star polymers during the polymerization. This behavior was attributed to Friedel–Crafts alkylation of the initiator core after the addition of one IB unit, which was activated by the electron‐donating oxytrimethyleneoxy linking moiety. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5942–5953, 2004     

Dielectric relaxation analysis of starch oligomers and polymers with respect to their chain length

Starch belongs to the polyglucan group. This type of polysaccharide shows a broad β-relaxation process in dielectric spectra at low temperatures, which has its molecular origin in orientational motions of sugar rings via glucosidic linkages. This chain dynamic was investigated for α(1,4)-linked starch oligomers with well-defined chain lengths of 2, 3, 4, 6, and 7 anhydroglucose units (AGUs) and for α(1,4)-polyglucans with average degrees of polymerization of 5, 10, 56, 70, and so forth (up to 3000; calculated from the mean molecular weight). The activation energy (E_a) of the segmental chain motion was lowest for dimeric maltose (E_a = 49.4 ± 1.3 kJ/mol), and this was followed by passage through a maximum at a degree of polymerization of 6 (E_a = 60.8 ± 1.8 kJ/mol). Subsequently, E_a leveled off at a value of about 52 ± 1.5 kJ/mol for chains containing more than 100 repeating units. The results were compared with the values of cellulose-like oligomers and polymers bearing a β(1,4)-linkage. Interestingly, the shape of the E_a dependency on the chain length of the molecules was qualitatively the same for both systems, whereas quantitatively the starch-like substances generally showed higher E_a values. Additionally, and for comparison, three cyclodextrins were measured by dielectric relaxation spectroscopy. The ringlike molecules, with 6, 7, and 8 α(1,4)-linked AGUs, showed moderately different types of dielectric spectra. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 188–197, 2004     

Control of self‐organized low‐dimensional morphology in poly(styrene‐b‐4vinylpyridine)/polystyrene blend thin films

The surface morphologies of poly(styrene‐b‐4vinylpyridine) (PS‐b‐P4VP) diblock copolymer and homopolystyrene (hPS) binary blend thin films were investigated by atomic force microscopy as a function of total volume fraction of PS (?_PS) in the mixture. It was found that when hPS was added into symmetric PS‐b‐P4VP diblock copolymers, the surface morphology of this diblock copolymer was changed to a certain degree. With ?_PS increasing at first, hPS was solubilized into the corresponding domains of block copolymer and formed cylinders. Moreover, the more solubilized the hPS, the more cylinders exist. However, when the limit was reached, excessive hPS tended to separate from the domains independently instead of solubilizing into the corresponding domains any longer, that is, a macrophase separation occurred. A model describing transitions of these morphologies with an increase in ?_PS is proposed. The effect of composition on the phase morphology of blend films when graphite is used as a substrate is also investigated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3496–3504, 2004     

Electrochemical reduction of 1-([(4-halophenyl)imino]methyl)-2-naphthols in aprotic media.

The electrochemical reduction of 1-([(4-halophenyl)imino]methyl)-2-naphthols on graphite electrodes was studied using cyclic voltammetry, chronoamperometry, constant-potential coulometry and preparative constant-potential electrolysis techniques. The data revealed that the reduction on graphite was irreversible and followed an EC mechanism. The diffusion coefficients and the number of electrons transferred were determined using the chronoamperometric Cottrell slope and the ultramicro disc Pt-electrode steady-state current. The number of electrons was also determined by bulk electrolysis. The compounds were subjected to constant-potential preparative electrolysis and the electrolysis products were purified and identified by spectroscopic methods. Based on these findings, a mechanism for the electro-reduction process is proposed.