Structural properties of charge-stabilized ferrofluids under a magnetic field: a Brownian dynamics study

We present Brownian dynamics simulations of real charge-stabilized ferrofluids, which are stable colloidal dispersions of magnetic nanoparticles, with and without the presence of an external magnetic field. The colloidal suspensions are treated as collections of monodisperse spherical particles, bearing point dipoles at their centers and undergoing translational and rotational Brownian motions. The overall repulsive isotropic interactions between particles, governed by electrostatic repulsions, are taken into account by a one-component effective pair interaction potential. The potential parameters are fitted in order that computed structure factors are close to the experimental ones. Two samples of ferrofluid differing by the particle diameter and consequently by the intensity of the magnetic interaction are considered here. The magnetization and birefringence curves are computed: a deviation from the ideal Langevin behaviors is observed if the dipolar moment of particles is sufficiently large. Structure factors are also computed from simulations with and without an applied magnetic field H: the microstructure of the repulsive ferrofluid becomes anisotropic under H. Even our simple modeling of the suspension allows us to account for the main experimental features: an increase of the peak intensity is observed in the direction perpendicular to the field whereas the peak intensity decreases in the direction parallel to the field.     

Intermolecular multiple quantum coherences at high magnetic field: the nonlinear regime

Experiments have been carried at magnetic-field strengths of 9.4, 14.1, and 17.6 T to explore the evolution of intermolecular multiple quantum coherences in the nonlinear regime where the system evolves for times that are much greater than the characteristic time of action of the long-range dipolar field, tau(d). The results show the expected Bessel function form of the recorded signal as a function of time of evolution, with evident zeros and sign changes. As expected, the rate of signal evolution increases at higher-field strengths as a result of the increased equilibrium magnetization. A numerical method for calculating the evolution of magnetization under the action of the distant dipolar field, relaxation, and diffusion that is based on Fourier analysis of the magnetization distribution has been applied to the correlated two-dimensional spectroscopy revamped by asymmetric z-gradient echo detection sequence in the nonlinear regime and shown to produce results that are in good agreement with experimental data acquired at different magnetic fields and rates of spatial modulation. Experiments and simulations have also been used to explore the evolution of magnetization in a mixture of two interacting spin species in the nonlinear regime.     

Electroptical and dynamic properties of sulfonated polystyrene molecules and their associates in chloroform

The behavior of sulfonated PS containing 0.5, 1.35, 2.6, and 5.8 mol % of sodium sulfonate groups in chloroform solutions has been studied by static and dynamic light scattering, viscometry, and electric birefringence. The molecular mass of ionomers is measured, and their translational diffusion coefficient, intrinsic viscosity, and free relaxation times are estimated. It has been shown that association in solutions of ionomers containing more than 1.35 mol % of sodium sulfonate groups proceeds according to the open association model. Analysis of autocorrelation functions of scattered light intensity and electric birefringence decay makes it possible to determine translational diffusion coefficients and relaxation times for individual ionomer molecules, their pair associates, and higher multiplicity associates. With an increase in the fraction of sodium sulfonate groups, the hydrodynamic radius of individual ionomer molecules decreases from 8 to 5.8 nm, while the ratio between the hydrodynamic radius of pair associates and individual sulfonated PS molecules increases.     

Hydrodynamic interaction in polymer solutions simulated with dissipative particle dynamics

The authors analyzed extensively the dynamics of polymer chains in solutions simulated with dissipative particle dynamics (DPD), with a special focus on the potential influence of a low Schmidt number of a typical DPD fluid on the simulated polymer dynamics. It has been argued that a low Schmidt number in a DPD fluid can lead to underdevelopment of the hydrodynamic interaction in polymer solutions. The authors' analyses reveal that equilibrium polymer dynamics in dilute solution, under typical DPD simulation conditions, obey the Zimm [J. Chem. Phys. 24, 269 (1956)] model very well. With a further reduction in the Schmidt number, a deviation from the Zimm model to the Rouse model is observed. This implies that the hydrodynamic interaction between monomers is reasonably developed under typical conditions of a DPD simulation. Only when the Schmidt number is further reduced, the hydrodynamic interaction within the chains becomes underdeveloped. The screening of the hydrodynamic interaction and the excluded volume interaction as the polymer volume fraction is increased are well reproduced by the DPD simulations. The use of soft interaction between polymer beads and a low Schmidt number do not produce noticeable problems for the simulated dynamics at high concentrations, except for the entanglement effect which is not captured in the simulations.     

Tracer diffusion in colloidal suspensions under dilute and crowded conditions with hydrodynamic interactions

We consider tracer diffusion in colloidal suspensions under solid loading conditions, where hydrodynamic interactions play an important role. To this end, we carry out computer simulations based on the hybrid stochastic rotation dynamics-molecular dynamics (SRD-MD) technique. Many details of the simulation method are discussed in detail. In particular, our choices for the SRD-MD parameters and for the different scales are adapted to simulating colloidal suspensions under realistic conditions. Our simulation data are compared with published theoretical, experimental and numerical results and compared to Brownian dynamics simulation data. We demonstrate that our SRD-MD simulations reproduce many features of the hydrodynamics in colloidal fluids under finite loading. In particular, finite-size effects and the diffusive behavior of colloids for a range of volume fractions of the suspension show that hydrodynamic interactions are correctly included within the SRD-MD technique.     

Electrostatic relaxation and hydrodynamic interactions for self-diffusion of ions in electrolyte solutions

The concentration dependence of self-diffusion of ions in solutions at large concentrations has remained an interesting yet unsolved problem. Here we develop a self-consistent microscopic approach based on the ideas of mode-coupling theory. It allows us to calculate both contributions which influence the friction of a moving ion: the ion atmosphere relaxation and hydrodynamic interactions. The resulting theory provides an excellent agreement with known experimental results over a wide concentration range. Interestingly, the mode-coupling self-consistent calculation of friction reveal a nonlinear coupling between the hydrodynamic interactions and the ion atmosphere relaxation which enhances ion diffusion by reducing friction, particularly at intermediate ion concentrations. This rather striking result has its origin in the similar time scales of the relaxation of the ion atmosphere relaxation and the hydrodynamic term, which are essentially given by the Debye relaxation time. The results are also in agreement with computer simulations, with and without hydrodynamic interactions.     

Dipolar solvation dynamics in room temperature ionic liquids: an effective medium calculation using dielectric relaxation data

Solvation dynamics in four imidazolium cation based room temperature ionic liquids (RTIL) have been calculated by using the recently measured dielectric relaxation data [ J. Phys. Chem. B 2008, 112, 4854 ] as an input in a molecular hydrodynamic theory developed earlier for studying solvation energy relaxation in polar solvents. Coumarin 153 (C153), 4-aminophthalimide (4-AP), and trans-4-dimethylamino-4'-cyanostilbene (DCS) have been used as probe molecules for this purpose. The medium response to a laser-excited probe molecule in an ionic liquid is approximated by that in an effective dipolar medium. The calculated decays of the solvent response function for these RTILs have been found to be biphasic and the decay time constants agree well with the available experimental and computer simulation results. Also, no probe dependence has been found for the average solvation times in these ionic liquids. In addition, dipolar solvation dynamics have been predicted for two other RTILs for which experimental results are not available yet. These predictions should be tested against experiments and/or simulation studies.     

Static and pulsed field gradient nuclear magnetic resonance studies of water diffusion in protein matrices

Static field gradient and pulsed field gradient NMR are used to study the temperature dependence of water diffusion in myoglobin and lysozyme matrices for low hydration levels of about 0.3?g/g. We show that in order to determine reliable self-diffusion coefficients D in a broad temperature range, it is very important to consider an exchange of magnetization between water and protein protons, often denoted as cross relaxation. Specifically, upon cooling, the observed stimulated-echo decays, which reflect water diffusion near ambient temperature, become more and more governed by cross relaxation. We demonstrate that comparison of experimental results for inhomogeneous and homogeneous magnetic fields enables successful separation of diffusion and relaxation contributions to the stimulated-echo decays. Making use of this possibility, we find that in the temperature range 230-300?K, the temperature-dependent diffusivities D exhibit a Vogel-Fulcher-Tammann behavior, where water diffusion in the studied protein matrices is substantially slower than in the bulk. By comparing present and previous data, we discuss relations between translational and rotational motions and between short-range and long-range water dynamics in protein matrices. In addition, we critically examine the significance of results from previous applications of NMR diffusometry to the temperature-dependent water diffusion in protein matrices.     

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard sphere particles in a dense monodisperse fluid. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from molecular dynamics simulations. At distances exceeding ~3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.     

Polarizable particle aggregation under rotating magnetic fields using scattering dichroism

We used scattering dichroism to study the dynamics of dipolar chains induced in magnetorheological suspensions under rotating magnetic fields. Both the dichroism (proportional to the total number of aggregated particles) and the phase lag show different behavior below and above a cross-over frequency. The cross-over frequency depends linearly on both the square of the magnetization and the inverse of the viscosity. The Mason number (ratio of viscous to magnetic forces) governs the dynamics. Therefore, there is a cross-over Mason number below which the dichroism remains almost constant and above which the rotation of the field prevents the particle aggregation process from taking place. Our experimental results have been compared with particle dynamics simulations showing good agreement.     

Activated hopping and dynamical fluctuation effects in hard sphere suspensions and fluids

Single particle Brownian dynamics simulation methods are employed to establish the full trajectory level predictions of our nonlinear stochastic Langevin equation theory of activated hopping dynamics in glassy hard sphere suspensions and fluids. The consequences of thermal noise driven mobility fluctuations associated with the barrier hopping process are determined for various ensemble-averaged properties and their distributions. The predicted mean square displacements show classic signatures of transient trapping and anomalous diffusion on intermediate time and length scales. A crossover to a stronger volume fraction dependence of the apparent nondiffusive exponent occurs when the entropic barrier is of order the thermal energy. The volume fraction dependences of various mean relaxation times and rates can be fitted by empirical critical power laws with parameters consistent with ideal mode-coupling theory. However, the results of our divergence-free theory are largely a consequence of activated dynamics. The experimentally measurable alpha relaxation time is found to be very similar to the theoretically defined mean reaction time for escape from the barrier-dominated regime. Various measures of decoupling have been studied. For fluid states with small or nonexistent barriers, relaxation times obey a simple log-normal distribution, while for high volume fractions the relaxation time distributions become Poissonian. The product of the self-diffusion constant and mean alpha relaxation time increases roughly as a logarithmic function of the alpha relaxation time. The cage scale incoherent dynamic structure factor exhibits nonexponential decay with a modest degree of stretching. A nearly universal collapse of the different volume fraction results occurs if time is scaled by the mean alpha relaxation time. Hence, time-volume fraction superposition holds quite well, despite the presence of stretching and volume fraction dependent decoupling associated with the stochastic barrier hopping process. The relevance of other origins of dynamic heterogeneity (e.g., mesoscopic domains), and comparison of our results with experiments, simulations, and alternative theories, is discussed.     

High resolution NMR study of T(1) magnetic relaxation dispersion. III. Influence of spin 1∕2 hetero-nuclei on spin relaxation and polarization transfer among strongly coupled protons

Effects of spin-spin interactions on the nuclear magnetic relaxation dispersion (NMRD) of protons were studied in a situation where spin [fraction one-half] hetero-nuclei are present in the molecule. As in earlier works [K. L. Ivanov, A. V. Yurkovskaya, and H.-M. Vieth, J. Chem. Phys. 129, 234513 (2008); S. E. Korchak, K. L. Ivanov, A. V. Yurkovskaya, and H.-M. Vieth, ibid. 133, 194502 (2010)], spin-spin interactions have a pronounced effect on the relaxivity tending to equalize the longitudinal relaxation times once the spins become strongly coupled at a sufficiently low magnetic field. In addition, we have found influence of (19)F nuclei on the proton NMRD, although in the whole field range, studied protons and fluorine spins were only weakly coupled. In particular, pronounced features in the proton NMRD were found; but each feature was predominantly observed only for particular spin states of the hetero-nuclei. The features are explained theoretically; it is shown that hetero-nuclei can affect the proton NMRD even in the limit of weak coupling when (i) protons are coupled strongly and (ii) have spin-spin interactions of different strengths with the hetero-nuclei. We also show that by choosing the proper magnetic field strength, one can selectively transfer proton spin magnetization between spectral components of choice.     

Transient electric birefringence of wormlike macromolecules in electric fields of arbitrary strength: a computer simulation study

We have studied the birefringence decay of linear models of macromolecules for two different types of flexibility, the broken-rod chain and the wormlike chain, using a computer simulation of a transient electric birefringence experiment. We have paid particular attention to the influence of the intensity of the orienting field, including two orienting mechanisms, the induced dipole, and the permanent dipole. We have compared wormlike and broken-rod models of the same radius of gyration, finding that they present a different decay curve under the influence of the same intensity of the field. We have seen that these differences are due to the faster relaxation times (smaller in the wormlike chain model) and amplitudes, because, regardless of the type of flexibility, the overall size of a molecule (measured by the radius of gyration) essentially determines the longest relaxation time. We have also analyzed how the relaxation process is affected by the degree of flexibility, the orientation mechanisms, and the intensity of the field. Studying a different aspect, we have paid attention to the deformation of a molecule in a transient electric birefringence experiment as a source of information. In this work we have developed equations to characterize this deformation in terms of one of the components of the gyration tensor, if a dynamic light scattering experiment under the influence of an electric field could be performed. To develop this work we have simulated the Brownian dynamics of the different models, relaxing after the removal of an orienting external electric field of arbitrary strength. A comparison with other methods such a the rigid body treatment or the correlation analysis of Brownian trajectories has also been included. We have seen that differences between the two Brownian dynamics methods are small and that the rigid-body treatment is only an acceptable approximation to obtain the longest relaxation time.     

良溶剂中环形高分子单链的结构与动力学性质的模拟研究

本文采用多粒子碰撞动力学与分子动力学耦合的模拟方法研究了环形高分子单链在良溶剂中的静态与动态性质,并与线形分子进行了对比.研究发现,环形高分子链内粒子之间的平均距离小于线形链,即粒子排列得更加紧密;相应的均方回转半径也小于线形链,线形链与环形链的均方回转半径的比值为1.77;同时,环形链扩散的速度也比线形链快,两者比值为1.10.模拟结果揭示了扩散行为是排斥体积作用和流体力学相互作用耦合的结果,在扩散过程中,流体力学相互作用消减了排斥体积作用对扩散行为的贡献.此外,通过对有和没有流体力学相互作用的多粒子碰撞动力学得到的结果作对比,研究了流体力学相互作用对高分子静态和动态行为的影响,结果表明,流体力学相互作用使高分子链在极稀溶液中的扩散速度变快.     

Rotational effects on the magnetic properties of an ensemble of free metal clusters

The magnetization of an ensemble of free magnetic metal clusters in an inhomogeneous external magnet field is calculated. In particular we have investigated the effects of the combined lattice anisotropy and cluster rotation on the magnetic properties. If weak anisotropy is present, almost superparamagnetic behavior is obtained. For stronger anisotropies deviations from this are calculated as a consequence of spin resonance due to the anisotropy field and the cluster rotation. This was proposed recently by de Heer et al. to explain his experimental data as generally expected, since a rotating cluster in a static magnetic field should behave similarly than a nonrotating one in an oscillating magnetic field. The magnetization depends also sensitively on the relaxation times.     

Low-temperature susceptibility of concentrated magnetic fluids

The initial susceptibility of concentrated magnetic fluids (ferrocolloids) has been experimentally investigated at low temperatures. The results obtained indicate that the interparticle dipole-dipole interactions can increase the susceptibility by several times as compared to the Langevin value. It is shown that good agreement between recent theoretical models and experimental observations can be achieved by introducing a correction for coefficients in the series expansion of susceptibility in powers of density and aggregation parameter. A modified equation for equilibrium susceptibility is offered to sum over corrections made by Kalikmanov (Statistical Physics of Fluids, Springer-Verlag, Berlin, 2001) and by B. Huke and M. Lucke (Phys. Rev. E 67, 051403, 2003). The equation gives good quantitative agreement with the experimental data in the wide range of temperature and magnetic particles concentration. It has been found that in some cases the magnetic fluid solidification occurs at temperature several tens of kelvins higher than the crystallization temperature of the carrier liquid. The solidification temperature of magnetic fluids is independent of particle concentration (i.e., magneto-dipole interparticle interactions) and dependent on the surfactant type and carrier liquid. This finding allows us to suggest that molecular interactions and generation of some large-scale structure from colloidal particles in magnetic fluids are responsible for magnetic fluid solidification. If the magnetic fluid contains the particles with the Brownian relaxation mechanism of the magnetic moment, the solidification manifests itself as the peak on the "susceptibility-temperature" curve. This fact proves the dynamic nature of the observed peak: it arises from blocking the Brownian mechanism of the magnetization relaxation.     

Collective relaxation of protein protons at very low magnetic field: a new window on protein dynamics and aggregation

Since the recent availability of high sensitivity field-cycling relaxometers, it has become possible to measure the protein proton relaxation in millimolar protein solutions as a function of magnetic field. In principle, this provides direct access to the so-called spectral density function of protein protons and, hence, to a full set of dynamic parameters. Understanding the dynamic behavior of biological molecules is increasingly appreciated as crucial to understanding their function. However, theoretical tools to analyze the collective relaxation behavior of protons in solute macromolecules over a wide range of magnetic fields are lacking. A complete relaxation matrix analysis of such behavior is described here. This analysis provides excellent predictions of the experimental proton magnetization decays/recoveries-measured to an unprecedented level of accuracy by a last-generation fast field-cycling relaxometer-of two different globular proteins, hen egg white lysozyme and human serum albumin. The new experimentally validated theoretical model is then used to extract dynamic information on these systems. A "collective" order parameter SC2, different from, but complementary to, that commonly extracted from heteronuclear relaxation measurements at high field, is defined and measured. An accurate estimate of the rotational correlation time is obtained: in the case of lysozyme it agrees very well with theoretical predictions; in the case of serum albumin it provides evidence for aggregation at millimolar concentration.     

Strong effect of hydrodynamic coupling on the electric dichroism of bent rods

The effect of hydrodynamic coupling on the spatial orientation of rigid bent rods in electric fields has been analyzed by Brownian dynamics simulations. Bead models for smoothly bent rods were constructed with dimensions of DNA double helices, and established simulation procedures were used to calculate their diffusion tensor, including the translational-rotational coupling tensor. The electric and optical parameters were assigned on the basis of known properties of double helices. Brownian dynamics simulations of the orientation of these models in electric fields showed that both transients and amplitudes of the calculated dichroism are very strongly dependent on translational-rotational coupling over a wide range of electric field strengths. For example, the stationary dichroism of a smoothly bent 179 bp DNA fragment calculated at low field strengths is positive in the presence and negative in the absence of hydrodynamic coupling. The transients are converted from a biphasic to a monophasic shape, when hydrodynamic coupling is turned off. The large changes resulting from hydrodynamic coupling were controlled by calculations based on analytical expressions derived for electrooptical response curves in the limit of low electric field strengths; the results obtained by this independent approach are in very satisfactory agreement with our Brownian dynamics simulations. The effect is strongly dependent on the electric dipole and on its direction. In the absence of any dipole the coupling effect was not observed. The coupling effect increases with the size of the bent rods. Because most macromolecular structures are known to have induced and/or permanent dipole moments, large effects of hydrodynamic coupling on both the amplitudes and the transients of the electric dichroism/birefringence must be expected in general for structures with nonsymmetric shape.     

Mesoscopic model for diffusion-influenced reaction dynamics

A hybrid mesoscopic multiparticle collision model is used to study diffusion-influenced reaction kinetics. The mesoscopic particle dynamics conserves mass, momentum, and energy so that hydrodynamic effects are fully taken into account. Reactive and nonreactive interactions with catalytic solute particles are described by full molecular dynamics. Results are presented for large-scale, three-dimensional simulations to study the influence of diffusion on the rate constants of the A + C <==> B + C reaction. In the limit of a dilute solution of catalytic C particles, the simulation results are compared with diffusion equation approaches for both the irreversible and reversible reaction cases. Simulation results for systems where the volume fraction phi of catalytic spheres is high are also presented, and collective interactions among reactions on catalytic spheres that introduce volume fraction dependence in the rate constants are studied.