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.

Levitation Effect: Role of Symmetry and Dependence of Diffusivity on the Bond Length of Homonuclear and Heteronuclear Diatomic Species

Molecular dynamics investigation of model diatomic species confined to the alpha-cages of zeolite NaY is reported. The dependence of self-diffusivity on the bond length of the diatomic species has been investigated. Three different sets of runs have been carried out. In the first set, the two atoms of the diatomic molecule interact with the zeolite atoms with equal strength (example, O-2, the symmetric case). In the second and third sets which correspond to asymmetric cases, the two atoms of the diatomic molecule interact with unequal strengths (example, CO). The result for the symmetric case exhibits a well-defined maximum in self-diffusivity for an intermediate bond length. In contrast to this, the intermediate asymmetry leads to a less pronounced maximum. For the large asymmetric case, the maximum is completely absent. These findings are analyzed by computing a number of related properties. These results provide a direct confirmation at the microscopic level of the suggestion by Derouane that the supermobility observed experimentally by Kemball has its origin in the mutual cancellation of forces. The maximum in diffusivity from molecular dynamics is seen at the value predicted by the levitation effect. Further, these findings suggest a role for symmetry in the existence of a diffusivity maximum as a function of diameter of the diffusant often referred to as the levitation effect. The nature of the required symmetry for the existence of anomalous diffusivity is interaction symmetry which is different from that normally encountered in crystallography.

Effect of Alkyl Chain Arrangement on Conformation and Dynamics in a Surfactant Intercalated Layered Double Hydroxide: Spectroscopic Measurements and MD Simulations

The anionic surfactant dodecyl sulfate (DDS) has been intercalated in an Mg-Al layered double hydroxide (LDH). Monolayer and bilayer arrangements of the alkyl chains of the intercalated surfactant can be engineered by tuning the Al/Mg ratio of the LDH. In both arrangements the anionic headgroup of the surfactant is tethered to the LDH sheets, and consequently translational mobility of the chains is absent. The degrees of freedom of the confined alkyl chains are restricted to changes in conformation. The effects of the arrangement of the intercalated surfactant chains on conformational order and dynamics have been,investigated by spectroscopic measurements and molecular dynamics simulations. Infrared, Raman, and C-13 NMR spectroscopies were used to investigate conformation of the alkyl chains in the monolayer and bilayer arrangements and variable contact time cross-polarization magic angle spinning (VCT CPMAS) NMR measurements to probe molecular motion. The alkyl chains in the monolayer arrangement of the intercalated DDS chains showed considerably greater conformational disorder and faster dynamics as compared to chains in the bilayer arrangement, in spite of the fact that the volume available per chain in the monolayer is smaller than that in the bilayer. Atomistic MD simulations of the two arrangements of the intercalated surfactant were carried out using an isothermal-isobaric ensemble. The simulations are able to reproduce the essential results of the experiment-greater conformational disorder and faster dynamics for the alkyl chains in the monolayer arrangement of the intercalated surfactant. The MD simulations show that these results are a consequence of the fact that the nature of conformational disorder in the two arrangements is different. In the monolayer arrangement the alkyl chains can sustain isolated gauche defects, whereas in the bilayer arrangement gauche conformers occur only as part of a kink a gauche(+) trans gauche(-) sequence.

Effect of two-level systems on the spectral density of universal conductance fluctuations in doped Si

We have studied the power spectral density [S(f) = gamma/f(alpha)] of universal conductance fluctuations (UCF's) in heavily doped single crystals of Si, when the scatterers themselves act as the primary source of dephasing. We observed that the scatterers, with internal dynamics like two-level-systems, produce a significant, temperature-dependent reduction in the spectral slope alpha when T less than or similar to 10 K, as compared to the bare 1/f (alphaapproximate to1) spectrum at higher temperatures. It is further shown that an upper cutoff frequency (f(m)) in the UCF spectrum is necessary in order to restrict the magnitude of conductance fluctuations, [(deltaG(phi))(2)], per phase coherent region (L-phi(3)) to [(deltaGphi)(2)](1/2) less than or similar to e(2)/h. We find that f(m) approximate to tau(D)(-1), where tau(D) = L-2/D, is the time scale of the diffusive motion of the electron along the active length (L) of the sample (D is the electron diffusivity).

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Effect of granularity on the insulator-superconductor transition in ultrathin Bi films

We have studied the insulator-superconductor transition (IST) by tuning the thickness in quench-condensed Bi films. The resistive transitions of the superconducting films are smooth and can be considered to represent ''homogeneous'' films. The observation of an IST very close to the quantum resistance for pairs R-square(N) similar to h/4e(2) on several substrates supports this idea. The relevant length scales here are the localization length, and the coherence length. However, at the transition, the localization length is much higher than the superconducting coherence length, contrary to expectation for a "homogeneous" transition. This suggests the invalidity of a purely fermionic model for the transition. Furthermore, the current-voltage characteristics of the superconducting films are hysteretic, and show the films to be granular. The relevant energy scales here are the Josephson coupling energy and the charging energy. However, Josephson coupling energies (E-J) and the charging energies (E-c) at the IST, they are found to obey the relation E-j < E-c. This is again contrary to expectation, for the IST in a granular or inhomogeneous system. Hence, a purely bosonic picture of the transition is also inconsistent with our observations. We conclude that the IST observed in our experiments may be either an intermediate case between the fermioinc and bosonic mechanisms, or in a regime of charge and vortex dynamics for which a quantitative analysis has not yet been done.

Albedo enhancement of marine clouds to counteract global warming: impacts on the hydrological cycle

Recent studies have shown that changes in solar radiation affect the hydrological cycle more strongly than equivalent CO(2) changes for the same change in global mean surface temperature. Thus, solar radiation management ``geoengineering'' proposals to completely offset global mean temperature increases by reducing the amount of absorbed sunlight might be expected to slow the global water cycle and reduce runoff over land. However, proposed countering of global warming by increasing the albedo of marine clouds would reduce surface solar radiation only over the oceans. Here, for an idealized scenario, we analyze the response of temperature and the hydrological cycle to increased reflection by clouds over the ocean using an atmospheric general circulation model coupled to a mixed layer ocean model. When cloud droplets are reduced in size over all oceans uniformly to offset the temperature increase from a doubling of atmospheric CO(2), the global-mean precipitation and evaporation decreases by about 1.3% but runoff over land increases by 7.5% primarily due to increases over tropical land. In the model, more reflective marine clouds cool the atmospheric column over ocean. The result is a sinking motion over oceans and upward motion over land. We attribute the increased runoff over land to this increased upward motion over land when marine clouds are made more reflective. Our results suggest that, in contrast to other proposals to increase planetary albedo, offsetting mean global warming by reducing marine cloud droplet size does not necessarily lead to a drying, on average, of the continents. However, we note that the changes in precipitation, evaporation and P-E are dominated by small but significant areas, and given the highly idealized nature of this study, a more thorough and broader assessment would be required for proposals of altering marine cloud properties on a large scale.

Effect of oxygen pressure on the grain and domain structure of polycrystalline 0.85PbMg(1/3)Nb(2/3)O(3)-0.15PbTiO(3) thin films studied by scanning probe microscopy

0.85PbMg(1/3)Nb(2/3)O(3)-0.15PbTiO(3) ferroelectric-relaxor thin films have been deposited on La(0.5)nSr(0.5)CoO(3)/(1 1 1) Pt/TiO(2)/SiO(2)/Si by pulsed laser ablation at various oxygen partial pressures in the range 0.05 to 0.4 Torr. All the films have a rhombohedral perovskite structure. The grain morphology and orientation are drastically affected by the oxygen pressure, studied by x-ray diffraction and scanning electron microscopy. The domain structure investigations by dynamic contact electrostatic force microscopy have revealed that the distribution of polar nanoregions and their dynamics is influenced by the grain morphology, orientation and more importantly, oxygen vacancies. The correlation length extracted from autocorrelation function images has shown that the polarization disorder decreases with oxygen pressure up to 0.3 Torr. The presence of polarized domains and their electric field induced switching is discussed in terms of internal bias field and domain wall pinning. Film deposited at 0.4 Torr presents a curious case with unique triangular grain morphology and large polarization disorder.

Dependence of climate forcing and response on the altitude of black carbon aerosols

Black carbon aerosols absorb solar radiation and decrease planetary albedo, and thus can contribute to climate warming. In this paper, the dependence of equilibrium climate response on the altitude of black carbon is explored using an atmospheric general circulation model coupled to a mixed layer ocean model. The simulations model aerosol direct and semi-direct effects, but not indirect effects. Aerosol concentrations are prescribed and not interactive. It is shown that climate response of black carbon is highly dependent on the altitude of the aerosol. As the altitude of black carbon increases, surface temperatures decrease; black carbon near the surface causes surface warming, whereas black carbon near the tropopause and in the stratosphere causes surface cooling. This cooling occurs despite increasing planetary absorption of sunlight (i.e. decreasing planetary albedo). We find that the trend in surface air temperature response versus the altitude of black carbon is consistent with our calculations of radiative forcing after the troposphere, stratosphere, and land surface have undergone rapid adjustment, calculated as ``regressed'' radiative forcing. The variation in climate response from black carbon at different altitudes occurs largely from different fast climate responses; temperature dependent feedbacks are not statistically distinguishable. Impacts of black carbon at various altitudes on the hydrological cycle are also discussed; black carbon in the lowest atmospheric layer increases precipitation despite reductions in solar radiation reaching the surface, whereas black carbon at higher altitudes decreases precipitation.

PAMAM Dendrimer-Drug Interactions: Effect of pH on the Binding and Release Pattern

Understanding the dendrimer-drug interaction is of great importance to design and optimize the dendrimer-based drug delivery system. Using atomistic molecular dynamics (MD) simulations, we have analyzed the release pattern of four ligands (two soluble drugs, namely, salicylic acid (Sal), L-alanine (Ala), and two insoluble drugs, namely, phenylbutazone (Pbz) and primidone (Prim)), which were initially encapsulated inside the ethylenediamine (EDA) cored polyamidoamine (PAMAM) dendrimer using the docking method. We have computed the potential of mean force (PMF) variation with generation 5 (G5)-PAMAM dendrimer complexed with drug molecules using umbrella sampling. From our calculated PMF values, we observe that soluble drugs (Sal and Ala) have lower energy barriers than insoluble drugs (Pbz and Prim). The order of ease of release pattern for these drugs from G5 protonated PAMAM dendrimer was found to be Ala > Sal > Prim > Pbz. In the case of insoluble drugs (Prim and Pbz), because of larger size, we observe much nonpolar contribution, and thus, their larger energy barriers can be reasoned to van der Waals contribution. From the hydrogen bonding analysis of the four PAMAM drug complexes under study, we found intermolecular hydrogen bonding to show less significant contribution to the free energy barrier. Another interesting feature appears while calculating the PMF profile of G5NP (nonprotonated)-PAMAM Pbz and G5NP (nonprotonated)-PAMAM-Sal complex. The PMF was found to be less when the drug is bound to nonprotonated dendrimer compared to the protonated dendrimer. Our results suggest that encapsulation of the drug molecule into the host PAMAM dendrimer should be carried out at higher pH values (near pH 10). When such complex enters the human body, the pH is around 7.4 and at that physiological pH, the dendrimer holds the drug tightly. Hence the release of drug can occur at a controlled rate into the bloodstream. Thus, our findings provide a microscopic picture of the encapsulation and controlled release of drugs in the case of dendrimer-based host-guest systems.

Confinement and viscoelastic effects on chain closure dynamics

Chemical reactions inside cells are typically subject to the effects both of the cell's confining surfaces and of the viscoelastic behavior of its contents. In this paper, we show how the outcome of one particular reaction of relevance to cellular biochemistry - the diffusion-limited cyclization of long chain polymers - is influenced by such confinement and crowding effects. More specifically, starting from the Rouse model of polymer dynamics, and invoking the Wilemski-Fixman approximation, we determine the scaling relationship between the mean closure time t(c) of a flexible chain (no excluded volume or hydrodynamic interactions) and the length N of its contour under the following separate conditions: (a) confinement of the chain to a sphere of radius d and (b) modulation of its dynamics by colored Gaussian noise. Among other results, we find that in case (a) when d is much smaller than the size of the chain, t(c) similar to Nd-2, and that in case (b), t(c) similar to N-2/(2 (2H)), H being a number between 1/2 and 1 that characterizes the decay of the noise correlations. H is not known a priori, but values of about 0.7 have been used in the successful characterization of protein conformational dynamics. At this value of H (selected for purposes of illustration), t(c) similar to N-3.4, the high scaling exponent reflecting the slow relaxation of the chain in a viscoelastic medium. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4729041]

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Effect of pressure on the ionic conductivity of Li+ and Cl- ions in water

A molecular dynamics simulation study of aqueous solution of LiCl is reported as a function of pressure. Experimental measurements of conductivity of Li+ ion as a function of pressure shows an increase in conductivity with pressure. Our simulations are able to reproduce the observed trend in conductivity. A number of relevant properties have been computed in order to understand the reasons for the increase in conductivity with pressure. These include radial distribution function, void and neck distributions, hydration or coordination numbers, diffusivity, velocity autocorrelation functions, angles between ion-oxygen and dipole of water as well as OH vector, mean residence time for water in the hydration shell, etc. These show that the increase in pressure acts as a structure breaker. The decay of the self part of the intermediate scattering function at small wave number k shows a bi-exponential decay at 1 bar which changes to single exponential decay at higher pressures. The k dependence of the ratio of the self part of the full width at half maximum of the dynamic structure factor to 2Dk(2) exhibits trends which suggest that the void structure of water is playing a role. These support the view that the changes in void and neck distributions in water can account for changes in conductivity or diffusivity of Li+ with pressure. These results can be understood in terms of the levitation effect. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4756909]

Measurement of temperature and pressure on the surface of a blunt cone using FBG sensor in hypersonic wind tunnel

Measurement of temperature and pressure exerted on the leeward surface of a blunt cone specimen has been demonstrated in the present work in a hypersonic wind tunnel using fiber Bragg grating (FBG) sensors. The experiments were conducted on a 30 degrees apex-angle blunt cone with 51 mm base diameter at wind flow speeds of Mach 6.5 and 8.35 in a 300 mm hypersonic wind tunnel of Indian Institute of Science, Bangalore. A special pressure insensitive temperature sensor probe along with the conventional bare FBG sensors was used for explicit temperature and aerodynamic pressure measurement respectively on the leeward surface of the specimen. computational fluid dynamics (CFD) simulation of the flow field around the blunt cone specimen has also been carried out to obtain the temperature and pressure at conditions analogous to experiments. The results obtained from FBG sensors and the CFD simulations are found to be in good agreement with each other.