Effects of Ultrafast Solvation on the Rate of Adiabatic Outer-Sphere Electron Transfer Reactions

Recent studies have demonstrated that solvation dynamics in many common dipolar liquids contain an initial, ultrafast Gaussian component which may contribute even more than 60% to the total solvation energy. It is also known that adiabatic electron transfer reactions often probe the high-frequency components of the relevant solvent friction (Hynes, J. T. J. Phys. Chem. 1986, 90, 3701). In this paper, we present a theoretical study of the effects of the ultrafast solvent polar modes on the adiabatic electron transfer reactions by using the formalism of Hynes. Calculations have been carried out for a model system and also for water and acetonitrile. It is found that, in general, the ultrafast modes can greatly enhance the rate of electron transfer, even by more than an order of magnitude, over the rate obtained by using only the slow overdamped modes usually considered. For water, this acceleration of the rate can be attributed to the high-frequency intermolecular vibrational and librational modes. For a weakly adiabatic reaction, the rate is virtually indistinguishable from the rate predicted by the Marcus transition state theory. Another important result is that even in this case of ultrafast underdamped solvation, energy diffusion appears to be efficient so that electron transfer reaction in water is controlled essentially by the barrier crossing dynamics. This is because the reactant well frequency is-directly proportional to the rate of the initial Gaussian decay of the solvation time correlation function. As a result, the value of the friction at the reactant well frequency rarely falls below the value required for the Kramers turnover except when the polarizability of the water molecules may be neglected. On the other hand, in acetonitrile, the rate of electron transfer reaction is found to be controlled by the energy diffusion dynamics, although a significant contribution to the rate comes also from the barrier crossing rate. Therefore, the present study calls for a need to understand the relaxation of the high-frequency modes in dipolar liquids.

Ultrafast underdamped solvation: Agreement between computer simulation and various theories of solvation dynamics

A theoretical analysis of the three currently popular microscopic theories of solvation dynamics, namely, the dynamic mean spherical approximation (DMSA), the molecular hydrodynamic theory (MHT), and the memory function theory (MFT) is carried out. It is shown that in the underdamped limit of momentum relaxation, all three theories lead to nearly identical results when the translational motions of both the solute ion and the solvent molecules are neglected. In this limit, the theoretical prediction is in almost perfect agreement with the computer simulation results of solvation dynamics in the model Stockmayer liquid. However, the situation changes significantly in the presence of the translational motion of the solvent molecules. In this case, DMSA breaks down but the other two theories correctly predict the acceleration of solvation in agreement with the simulation results. We find that the translational motion of a light solute ion can play an important role in its own solvation. None of the existing theories describe this aspect. A generalization of the extended hydrodynamic theory is presented which, for the first time, includes the contribution of solute motion towards its own solvation dynamics. The extended theory gives excellent agreement with the simulations where solute motion is allowed. It is further shown that in the absence of translation, the memory function theory of Fried and Mukamel can be recovered from the hydrodynamic equations if the wave vector dependent dissipative kernel in the hydrodynamic description is replaced by its long wavelength value. We suggest a convenient memory kernel which is superior to the limiting forms used in earlier descriptions. We also present an alternate, quite general, statistical mechanical expression for the time dependent solvation energy of an ion. This expression has remarkable similarity with that for the translational dielectric friction on a moving ion.

Application of microwave heating technique for rapid synthesis of ?,?-unsaturated esters

The use of microwave heating technique for the acceleration of ortho ester Claisen rearrangement (a three step transformation) is described. Irradiation of a DMF solution of the allyl alcohol 5, triethyl orthoacetate and propionic acid (catalytic) in an Erlenmeyer flask for 10 minutes in a microwave oven generated the ester 8 in 83% yield. Analogously, ortho ester Claisen rearrangement of a variety of allyl and propargyl alcohols (9, 12-22) were achieved. The formation of the diester 10 from 2-butyne-1,4-diol (9) via the ortho ester Claisen rearrangement of two allyl alcohol moieties (involving sh steps) in 15 minutes, demonstrates the versatility of the microwave heating technique.

Structure and properties of fibre reinforced zn-27% al alloy based cast MMCs

Metal matrix composites (MMCs) based on a zinc-27% aluminium alloy (ZA-27) were produced using a pressure infiltration technique. Preforms of alumina fibres and aluminosilicate fibres were used for reinforcement. Uniform distribution of fibres and satisfactory interfacial bonding were achieved. UTS, specific strength, hardness and wear resistance were improved significantly by the alumina fibre reinforcement, but UTS decreased when using aluminosilicate fibres for reinforcement mainly due to unavoidable clustering of particles in the fibre preforms. Structure-property relations have been analysed in all cases.

Structure-potential relationship and nature of disorder in liquids and glassy solids

A method based on the minimal-spanning tree is extended to a collection of points in three dimensions. Two parameters, the average edge length and its standard deviation characterize the disorder. The structural phase diagram for a monatomic system of particles and the characteristic values for the uniform random distribution of points have been obtained. The method is applied to hard spheres and Lennard-Jones systems. These systems occupy distinct regions in the structural phase diagram. The structure of the Lennard-Jones system approaches that of the defective close-packed arrangements at low temperatures whereas in the liquid regime, it deviates from the close-packed configuration.

Some Basic Aspects of Reaction Engineering of Precipitation Processes

Analysis of precipitation reactions is extremely important in the technology of production of fine particles from the liquid phase. The control of composition and particle size in precipitation processes requires careful analysis of the several reactions that comprise the precipitation system. Since precipitation systems involve several, rapid ionic dissociation reactions among other slower ones, the faster reactions may be assumed to be nearly at equilibrium. However, the elimination of species, and the consequent reduction of the system of equations, is an aspect of analysis fraught with the possibility of subtle errors related to the violation of conservation principles. This paper shows how such errors may be avoided systematically by relying on the methods of linear algebra. Applications are demonstrated by analyzing the reactions leading to the precipitation of calcium carbonate in a stirred tank reactor as well as in a single emulsion drop. Sample calculations show that supersaturation dynamics can assume forms that can lead to subsequent dissolution of particles that have once been precipitated.

Thomas-Fermi Method For Particles Obeying Generalized Exclusion Statistics

We use the Thomas-Fermi method to examine the thermodynamics of particles obeying Haldane exclusion statistics. Specifically, we study Calogero-Sutherland particles placed in a given external potential in one dimension. For the case of a simple harmonic potential (constant density of states), we obtain the exact one-particle spatial density and a {\it closed} form for the equation of state at finite temperature, which are both new results. We then solve the problem of particles in a $x^{2/3} ~$ potential (linear density of states) and show that Bose-Einstein condensation does not occur for any statistics other than bosons.

Modeling of precipitation in reverse micellar systems

A model of the precipitation process in reverse micelles has been developed to calculate the size of fine particles obtained therein. While the method shares several features of particle nucleation and growth common to precipitation in large systems, complexities arise in describing the processes of nucleation, due to the extremely small size of a micelle and of particle growth caused by fusion among the micelles. Occupancy of micelles by solubilized molecules is governed by Poisson statistics, implying most of them are empty and cannot nucleate of its own. The model therefore specifies the minimum number of solubilized molecules required to form a nucleus which is used to calculate the homogeneous nucleation rate. Simultaneously, interaction between micelles is assumed to occur by Brownian collision and instantaneous fusion. Analysis of time scales of various events shows growth of particles to be very fast compared to other phenomena occurring. This implies that nonempty micelles either are supersaturated or contain a single precipitated particle and allows application of deterministic population balance equations to describe the evolution of the system with time. The model successfully predicts the experimental measurements of Kandori ct al.(3) on the size of precipitated CaCO3 particles, obtained by carbonation of reverse micelles containing aqueous Ca(OH)(2) solution.

Microwave preparation of oxide bronzes

Several vanadium, tungsten, and molybdenum oxide bronzes have been prepared using microwave irradiation. Metal oxides and alkali metal iodides were used as starting materials, Intermittent grinding and inert atmosphere were found to be necessary for the synthesis of most of the bronzes, The reaction temperatures are remarkably lower than those employed for conventional synthetic techniques and the microwave assisted reactions proceed at extremely fast rates. The microwave synthesized bronzes consist of particles having long, rectangular rod-like morphology. (C) 1999 Academic Press.

Stability and mobility of sand-bed channels affected by seepage

Seepage effects on the stability, mobility, and incipient motion of sand-bed particles are experimentally investigated. Seepage through a sand bed in a downward direction (suction) reduces the stability of particles, and it can even initiate their movement. The bed erosion is increased with the increased rates of suction. Whereas the seepage in an upward direction (injection) increases the stability of bed particles, it does not aid initiating their movement. The rate of bed erosion is reduced or even stopped by the increased infection rates. Hydrodynamic conditions leading to the so-called "pseudoincipient motion'' with suction (for the initiation of particles movement that are otherwise at rest under no-seepage conditions), and with injection (for only arresting the particles movement that are otherwise moving initially) are evaluated. The conventional Shields curve cannot be used to predict such pseudoincipient motion conditions with seepage. The concepts thus developed are useful for a better understanding of the sediment transport mechanics and in the design of stable alluvial channels affected by seepage.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Glassy behavior in neural network models of associative memory

Neural network models of associative memory exhibit a large number of spurious attractors of the network dynamics which are not correlated with any memory state. These spurious attractors, analogous to "glassy" local minima of the energy or free energy of a system of particles, degrade the performance of the network by trapping trajectories starting from states that are not close to one of the memory states. Different methods for reducing the adverse effects of spurious attractors are examined with emphasis on the role of synaptic asymmetry. (C) 2002 Elsevier Science B.V. All rights reserved.

Melting-freezing cycles in a relatively sheared pair of crystalline monolayers

The nonequilibrium dynamical behaviour that arises when two ordered two-dimensional monolayers of particles are sheared over each other is studied in Brownian dynamics simulations. A curious sequence of nonequilibrium states is observed as the driving rate is increased, the most striking of which is a sliding state with irregular alternation between disordered and ordered states. We comment on possible mechanisms underlying these cycles, and experiments that could observe them.

Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.