Shear-induced enhancement of self-diffusion in interacting colloidal suspensions

We set up the generalized Langevin equations describing coupled single-particle and collective motion in a suspension of interacting colloidal particles in a shear how and use these to show that the measured self-diffusion coefficients in these systems should be strongly dependent on shear rate epsilon. Three regimes are found: (i) an initial const+epsilon(.2), followed by (ii) a large regime of epsilon(.1/2) behavior, crossing over to an asymptotic power-law approach (iii) D-o - const x epsilon(.-1/2) to the Stokes-Einstein value D-o. The shear dependence is isotropic up to very large shear rates and increases with the interparticle interaction strength. Our results provide a straightforward explanation of recent experiments and simulations on sheared colloids.

Self-diffusion of colloids at freezing

This is an introduction to the theory of interacting Brownian particles, as applied to charge-stabilised colloidal suspensions near their equilibrium liquid-solid transition. The density functional approach to the statics of the transition is reviewed briefly, and the generalised Langevin equation method for the dynamics presented in detail. Work with A.V. Indrani [1] on a self-consistent approach for calculating the excess single-particle friction is presented, which explains the observed [2] ''universal'' suppression of self-diffusion at freezing as a consequence of the universal structure-factor height at this transition. Criticisms, open questions, and challenges to theory are discussed.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Diffusion of Pure CH4 and Its Binary Mixture with CO2 in Faujasite NaY: A Combination of Neutron Scattering Experiments and Molecular Dynamics Simulations

The self-diffusion properties of pure CH4 and its binary mixture with CO2 within MY zeolite have been investigated by combining an experimental quasi-elastic neutron scattering (QENS) technique and classical Molecular dynamics simulations. The QENS measurements carried out at 200 K led to an unexpected self-diffusivity profile for Pure CH4 with the presence of a maximum for a loading of 32 CH4/unit cell, which was never observed before for the diffusion of apolar species in azeolite system With large windows. Molecular dynamics simulations were performed using two distinct microscopic models for representing the CH4/NaY interactions. Depending on the model, we are able to fairly reproduce either the magnitude or the profile of the self-diffusivity.Further analysis allowed LIS to provide some molecular insight into the diffusion mechanism in play. The QENS measurements report only a slight decrease of the self-diffusivity of CH4 in the presence of CO2 when the CO2 loading increases. Molecular dynamics simulations successfully capture this experimental trend and suggest a plausible microscopic diffusion mechanism in the case of this binary mixture.

Influence of the Methane-Zeolite A Interaction Potential on the Concentration Dependence of Self-diffusivity

Studies on the diffusion of methane in a zeolite structure type LTA (as per IZA nomenclature) have indicated that different types of methane zeolite potentials exist in the literature in which methane is treated within the united-atom model. One set of potentials, referred to as model A, has a methane oxygen diameter of 3.14 angstrom, while another set of potential parameters, model B, employs a larger value of 3.46 angstrom. Fritzsche and co-workers (1993) have shown that these two potentials lead to two distinctly different energetic barriers for the passage of methane through the eight-ring window in the cation-free form of zeolite A. Here, we compute the variation of the self-diffusivity (D) with loading (c) for these two types of potentials and show that this slight variation in the diameter changes the concentration dependence qualitatively: thus, D decreases monotonically with c for model A, while D increases and goes through a maximum before finally decreasing for model B. This effect and the surprising congruence of the diffusion coefficients for both models at high loadings is examined in detail at the molecular level. Simulations for different temperatures reveal the Arrhenius behaviour of the self-diffusion coefficient. The apparent activation energy is found to vary with the loading. We conclude that beside the cage-to-cage jumps, which are essential for the migration of the guest molecules, at high concentrations migration within the cage and guest guest interactions with other molecules become increasingly dominant influences on the diffusion coefficient and make the guest zeolite interaction less important for both model A and model B.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.

Creep Behaviour of AE42 Magnesium Alloy and Its Composites Using Impression Creep Technique

The creep behaviour of a creep-resistant AE42 magnesium alloy has been examined in the temperature range of 150 to 240 degrees C at the stress levels ranging from 40 to 120 MPa using impression creep technique. A normal creep behaviour, i.e., strain rate decreasing with strain and then reaching a steady state, is observed at all the temperatures and stresses employed The stress exponent varies from 5.1 to 5.7 and the apparent activation energy varies from 130 to 140 kJ/mol, which suggests the high temperature climb of dislocation controlled by lattice self-diffusion being the dominant creep mechanism in the stress and temperature range employed The creep behaviour of the AE42 alloy has also been compared with its composites reinforced with Saffil short fibres and SiC particles in four combinations. All the composites exhibited a lower creep rate than the monolithic AE42 alloy tested at the same temperature and stress levels and the decrease in creep rate was greater in the longitudinal direction than in the transverse direction, as expected. All the hybrid composites, i.e., the composites reinforced with a combination of Saffil short fibres and SiC particles, exhibited creep rates comparable to the composite reinforced with 20% Saffil short fibres alone at all the temperature and stress levels employed, which is beneficial from the commercial point of view.

Characterisation of dynamic recrystallisation in nickel using processing map for hot deformation

The hot deformation behaviour of polycrystalline nickel has been characterised in the temperature range 750-1200-degrees-C and strain rate range 0.0003-100 s-1 using processing maps developed in the basis of the dynamic materials model. The efficiency of power dissipation, given by [2m/(m + 1)]. where m is the strain rate sensitivity, is plotted as a function of temperature and strain rate to obtain a processing map. A domain of dynamic recrystallisation has been identified, with a peak efficiency of 31% occurring at 925-degrees-C and 1 s-1. The published results are in agreement with the prediction of the processing map. The variations of efficiency of power dissipation with temperature and strain rate in the dynamic recrystallisation domain are identical to the corresponding variation of hot ductility. The stress-strain curves exhibited a single peak in a single peak in the dynamic recrystallisation domain, whereas multiple peaks and 'drooping' stress-strain curves were observed at lower and higher strain rates, respectively. The results are explained on the basis of a simple model which considers dynamic recrystallisation in terms of rates of interface formation (nucleation) and migration (growth). It is shown that dynamic recrystallisation in nickel is controlled by the rate of nucleation, which is slower than the rate of migration. The rate of nucleation itself depends on the process of thermal recovery by climb, which in turn depends on self-diffusion.

Electrical transport in magnesium aluminate

The conductivity of MgAl2O4 has been measured at 1273, 1473 and 1673 K as a function of the partial pressure of oxygen ranging from 105 to 10−14 Pa. The MgAl2O4 pellet, sandwiched between two platinum electrodes, was equilibrated with a flowing stream of either Ar + O2, CO + CO2 or Ar + H2 + H2O mixture of known composition. The gas mixture established a known oxygen partial pressure. All measurements were made at a frequency of 1 kHz. These measurements indicate pressure independent ionic conductivity in the range 1 to 10−14 Pa at 1273 K, 10−1 to 10−12 Pa at 1473 K and 10−1 to 10−4 Pa at 1673 K. The activation energy for ionic conduction is 1·48 eV, close to that for self-diffusion of Mg2+ ion in MgAl2O4 calculated from the theoretical relation of Glyde. Using the model, the energy for cation vacancy formation and activation energy for migration are estimated.

Effect of strain rate on the high-temperature low-cycle fatigue properties of a nimonic PE-16 superalloy

Strain-rate effects on the low-cycle fatigue (LCF) behavior of a NIMONIC PE-16 superalloy have been evaluated in the temperature range of 523 to 923 K. Total-strain-controlled fatigue tests were per-formed at a strain amplitude of +/-0.6 pct on samples possessing two different prior microstructures: microstructure A, in the solution-annealed condition (free of gamma' and carbides); and microstructure B, in a double-aged condition with gamma' of 18-nm diameter and M23C6 carbides. The cyclic stress response behavior of the alloy was found to depend on the prior microstructure, testing temperature, and strain rate. A softening regime was found to be associated with shearing of ordered gamma' that were either formed during testing or present in the prior microstructure. Various manifestations of dynamic strain aging (DSA) included negative strain rate-stress response, serrations on the stress-strain hysteresis loops, and increased work-hardening rate. The calculated activation energy matched well with that for self-diffusion of Al and Ti in the matrix. Fatigue life increased with an increase in strain rate from 3 x 10(-5) to 3 x 10(-3) s-1, but decreased with further increases in strain rate. At 723 and 823 K and low strain rates, DSA influenced the deformation and fracture behavior of the alloy. Dynamic strain aging increased the strain localization in planar slip bands, and impingement of these bands caused internal grain-boundary cracks and reduced fatigue life. However, at 923 K and low strain rates, fatigue crack initiation and propagation were accelerated by high-temperature oxidation, and the reduced fatigue life was attributed to oxidation-fatigue interaction. Fatigue life was maximum at the intermediate strain rates, where strain localization was lower. Strain localization as a function of strain rate and temperature was quantified by optical and scanning electron microscopy and correlated with fatigue life.

The extended Enskog operator for simple fluids with continuous potentials: single particle and collective properties

We generalized the Enskog theory originally developed for the hard-sphere fluid to fluids with continuous potentials, such as the Lennard–Jones. We derived the expression for the k and ω dependent transport coefficient matrix which enables us to calculate the transport coefficients for arbitrary length and time scales. Our results reduce to the conventional Chapman–Enskog expression in the low density limit and to the conventional k dependent Enskog theory in the hard-sphere limit. As examples, the self-diffusion of a single atom, the vibrational energy relaxation, and the activated barrier crossing dynamics problem are discussed.

Molecular dynamics insight to phase transition in n-alkanes with carbon nanofillers

The present work aims to investigate the phase transition, dispersion and diffusion behavior of nanocomposites of carbon nanotube (CNT) and straight chain alkanes. These materials are potential candidates for organic phase change materials(PCMs) and have attracted flurry of research recently. Accurate experimental evaluation of the mass, thermal and transport properties of such composites is both difficult as well as economically taxing. Additionally it is crucial to understand the factors that results in modification or enhancement of their characteristic at atomic or molecular level. Classical molecular dynamics approach has been extended to elucidate the same. Bulk atomistic models have been generated and subjected to rigorous multistage equilibration. To reaffirm the approach, both canonical and constant-temperature, constant-pressure ensembles were employed to simulate the models under consideration. Explicit determination of kinetic, potential, non-bond and total energy assisted in understanding the enhanced thermal and transport property of the nanocomposites from molecular point of view. Crucial parameters including mean square displacement and simulated self diffusion coefficient precisely define the balance of the thermodynamic and hydrodynamic interactions. Radial distribution function also reflected the density variation, strength and mobility of the nanocomposites. It is expected that CNT functionalization could improve the dispersion within n-alkane matrix. This would further ameliorate the mass and thermal properties of the composite. Additionally, the determined density was in good agreement with experimental data. Thus, molecular dynamics can be utilized as a high throughput technique for theoretical investigation of nanocomposites PCMs. (C) 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Natural convection from a vertical cone in a porous medium due to the combined effects of heat and mass diffusion with non-uniform wall temperature/concentration or heat/mass flux and suction/injection

An analysis has been carried out to study the non-Darcy natural convention flow of Newtonian fluids on a vertical cone embedded in a saturated porous medium with power-law variation of the wall temperature/concentration or heat/mass flux and suction/injection with the streamwise distance x. Both non-similar and self-similar solutions have been obtained. The effects of non-Darcy parameter, ratio of the buoyancy forces due to mass and heat diffusion, variation of wall temperature/concentration or heat/mass flux and suction/injection on the Nusselt and Sherwood numbers have been studied.