Polymer Conformations in Ionic Microgels in the Presence of Salt: Theoretical and Mesoscale Simulation Results.

We investigate the conformational properties of polymers in ionic microgels in the presence of salt ions by molecular dynamics simulations and analytical theory. A microgel particle consists of coarse-grained linear polymers, which are tetra-functionally crosslinked. Counterions and salt ions are taken into account explicitly, and charge-charge interactions are described by the Coulomb potential. By varying the charge interaction strength and salt concentration, we characterize the swelling of the polyelectrolytes and the charge distribution. In particular, we determine the amount of trapped mobile charges inside the microgel and the Debye screening length. Moreover, we analyze the polymer extension theoretically in terms of the tension blob model taking into account counterions and salt ions implicitly by the Debye⁻Hückel model. Our studies reveal a strong dependence of the amount of ions absorbed in the interior of the microgel on the electrostatic interaction strength, which is related to the degree of the gel swelling. This implies a dependence of the inverse Debye screening length κ on the ion concentration; we find a power-law increase of κ with the Coulomb interaction strength with the exponent 3 / 5 for a salt-free microgel and an exponent 1 / 2 for moderate salt concentrations. Additionally, the radial dependence of polymer conformations and ion distributions is addressed.

Comparison of scalable fast methods for long-range interactions.

Based on a parallel scalable library for Coulomb interactions in particle systems, a comparison between the fast multipole method (FMM), multigrid-based methods, fast Fourier transform (FFT)-based methods, and a Maxwell solver is provided for the case of three-dimensional periodic boundary conditions. These methods are directly compared with respect to complexity, scalability, performance, and accuracy. To ensure comparable conditions for all methods and to cover typical applications, we tested all methods on the same set of computers using identical benchmark systems. Our findings suggest that, depending on system size and desired accuracy, the FMM- and FFT-based methods are most efficient in performance and stability.

Hydrodynamic fluctuations in thermostatted multiparticle collision dynamics.

In this work we study the behavior of mesoscopic fluctuations of a fluid simulated by Multiparticle Collision Dynamics when this is applied together with a local thermostatting procedure that constrains the strength of temperature fluctuations. We consider procedures in which the thermostat interacts with the fluid at every simulation step as well as cases in which the thermostat is applied only at regular time intervals. Due to the application of the thermostat temperature fluctuations are forced to relax to equilibrium faster than they do in the nonthermostatted, constant-energy case. Depending on the interval of application of the thermostat, it is demonstrated that the thermodynamic state changes gradually from isothermal to adiabatic conditions. In order to exhibit this effect we compute from simulations diverse correlation functions of the hydrodynamic fluctuating fields. These correlation functions are compared with those predicted by a linearized hydrodynamic theory of a simple fluid in which a thermostat is applied locally. We find a good agreement between the model and the numerical results, which confirms that hydrodynamic fluctuations in Multiparticle Collision Dynamics in the presence of the thermostat have the properties expected for spontaneous fluctuations in fluids in contact with a heat reservoir.

Two-dimensional chiral model for liquid crystals, bent hard needles: a Monte Carlo simulation.

The liquid crystalline behavior of a two dimensional (2D) model of hard needles bent into a "zigzag shape" is studied. This model, originally designed to study two dimensional chiral segregation, also shows liquid crystalline behavior and has some anomalous features which are contrasted in relation to the following: (i) Most of the microscopical models used to study liquid crystals have a symmetry axis that coincides with a molecular axis; (ii) in three-dimensions, chiral molecules can form cholesteric instead of nematic phases; (iii) the smectic phase is usually found when attractions are present or at least when the molecules have finite volume. Despite the fact that the present 2D model does not have any of these characteristics, numerical evidence is found for the occurrence of nematic and smectic phases. Since these molecules are athermal, infinitely repulsive, and infinitesimally thin, the liquid crystalline characteristics are attributed to excluded volume effects. To determine the mesophases of the model, both nematic and smectic order parameters as well as distribution functions are computed.

Is coulomb explosion a damaging mechanism for (125)IUdR?

PURPOSE: To test the integrity of the thymine molecule that experiences an increasing number of charges due to the loss of Auger electrons emitted by the decay of incorporated 125I. Besides the radiation action of these electrons, Coulomb explosion is suspected to be an additional mechanism responsible for the strong radiotoxic effect of decaying DNA-incorporated 125I. The two-step decay process initiates a first Auger cascade within 10(-16) to 10(-14) s resulting in the release of about 7 electrons on average and a corresponding large positive charge on the 125Te daughter atom. Being part of iododeoxyuridine (125IUdR), the analogue of the DNA base thymine, the base is suddenly confronted with this charge. Experimentally, the situation was investigated with small molecules (CH3(125)I and C2H5(125)I) resulting in ion fragmentation in agreement with a Coulomb explosion model (Carlson and White, 1963, 1966). MATERIALS AND METHODS: Semi-empirical quantum mechanical calculations on the Parametric Method 3 (PM3) level (Stewart, 1989a, 1989b) were performed and geometry optimisation was applied for the identification of stable molecule conformations. Subsequently, semiempirical molecular dynamics simulations allowed changes in the conformations to be studied as a function of time. RESULTS: First results show that there is no stable molecular configuration with a total charge of > or = +5e. PM3 calculations will not converge for such a charge located at the 125I/125Te position. This finding is supported by total energy considerations, which begin to favour a system of isolated atoms versus molecular bound atoms when the molecular charge is greater than +4e. The distribution of the partial charges indicates that most of the charge will remain on the tellurium atom with slight increases of charge at the other molecular partners within 125IUdR. Moreover, the molecular dynamics simulations reveal a breaking of chemical bonds between those atoms with the strongest charge increase. CONCLUSIONS: Coulomb explosion must be taken into account as a possible damaging mechanism following the decay of DNA-incorporated Auger electron emitters. Lobachevsky and Martin (2000) have identified the same mechanism to be responsible for part of strand breakage in oligo-deoxynucleotides. To elucidate a possible link between both damage patterns the molecular mechanics simulations have to be extended to larger parts of the DNA molecule.