Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.

Energy dispersive backscattering of electrons from surface resonances of a disordered medium and 1/f noise

Anderson localised states in the bulk of a disordered medium appear as sharp resonances near the surface. The resonant backscattering leads to an energy-dependent random time delay for an incident electron. We derive an analytic expression for the delay-time probability distribution at a given energy. This is shown to give a 1/f noise for the surface currents in general.

Surface modifications in single crystal surfaces of YBa2Cu3O7-delta upon high energy ion irradiation

Atomic force microscopy investigations on swift heavy ion (200 MeV An) irradiated surfaces of a high T-c single crystal YBa2Cu3O7-delta are presented. Results obtained revealed an ion-induced erosion/sputtering clearly confirming our earlier observation on grain boundary dominated thin films. Apart from sputtering, notable effects were seen with many defect structures like dikes/hillocks surrounded by craters, dikes, holes, pearl like structures and ripple formation of sub-micron undulations, all in one crystal. Results are discussed in the light of co-operative phenomena of material re-distribution mechanism related to mass transfer and crater formations.

Dynamics of polar solvation: Route to single exponential relaxation via translational diffusion

A microscopic theoretical calculation of time-dependent solvation energy shows that the solvation of an ion or a dipole is dominated by a single relaxation time if the translational contribution to relaxation is significant.

Polarization Caging in Diffusion-Controlled Electron Transfer Reactions in Solution

In some bimolecular diffusion-controlled electron transfer (ET) reactions such as ion recombination (IR), both solvent polarization relaxation and the mutual diffusion of the reacting ion pair may determine the rate and even the yield of the reaction. However, a full treatment with these two reaction coordinates is a challenging task and has been left mostly unsolved. In this work, we address this problem by developing a dynamic theory by combining the ideas from ET reaction literature and barrierless chemical reactions. Two-dimensional coupled Smoluchowski equations are employed to compute the time evolution of joint probability distribution for the reactant (P-(1)(X,R,t)) and the product (p((2))(X,R,t)), where X, as is usual in ET reactions, describes the solvent polarization coordinate and R is the distance between the reacting ion pair. The reaction is described by a reaction line (sink) which is a function of X and R obtained by imposing a condition of equal energy on the initial and final states of a reacting ion pair. The resulting two-dimensional coupled equations of motion have been solved numerically using an alternate direction implicit (ADI) scheme (Peaceman and Rachford, J. Soc. Ind. Appl. Math. 1955, 3, 28). The results reveal interesting interplay between polarization relaxation and translational dynamics. The following new results have been obtained. (i) For solvents with slow longitudinal polarization relaxation, the escape probability decreases drastically as the polarization relaxation time increases. We attribute this to caging by polarization of the surrounding solvent, As expected, for the solvents having fast polarization relaxation, the escape probability is independent of the polarization relaxation time. (ii) In the slow relaxation limit, there is a significant dependence of escape probability and average rate on the initial solvent polarization, again displaying the effects of polarization caging. Escape probability increases, and the average rate decreases on increasing the initial polarization. Again, in the fast polarization relaxation limit, there is no effect of initial polarization on the escape probability and the average rate of IR. (iii) For normal and barrierless regions the dependence of escape probability and the rate of IR on initial polarization is stronger than in the inverted region. (iv) Because of the involvement of dynamics along R coordinate, the asymmetrical parabolic (that is, non-Marcus) energy gap dependence of the rate is observed.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Distribution of ionic charge carriers and migration barriers in binary alkali silicate glasses

Likely spatial distributions of network-modifying (and mobile) cations in (oxide) glasses are discussed here. At very low modifier concentrations, the ions form dipoles with non-bridging oxygen centres while, at higher levels of modification, the cations tend to order as a result of Coulombic interactions. Activation energies for cation migration are calculated, assuming that the ions occupy (face-sharing) octahedral sites. It is found that conductivity activation energy decreases markedly with increasing modifier content, in agreement with experiment.

Non-Marcus energy gap dependence of electron transfer rate in contact ion pairs. Novel interplay between relaxation and reaction in solution

Electron transfer reactions between donor-acceptor pairs in solution and in organized media exhibit diverse behaviour. Recent experiments have indicated an interesting breakdown of the Marcus parabolic energy gap dependence in the normal regime for back electron transfer from contact ion pairs. A novel explanation of this breakdown has recently been proposed (M. Tachiya and S. Murata, J. Am. Chem. Sec., 116(1994) 2434) which attributes the breakdown to the interplay between the relaxation in the reactant well and the reaction. A particularly interesting aspect of the model is that it envisages the electron transfer in the normal regime to take place from a completely non-equilibrium condition. In this article a time dependent solution of the model is presented for the first time, after generalizing it to include a realistic initial population distribution. The decay of the contact ion pair population is completely non-exponential. This can be used to check the validity of the Tachiya-Murata model. The dynamics of electron transfer from the solvent separated ion pair, which seem to obey the Marcus relation, is exponential.

Solubilization of n-decane into n-decyl-polyoxyethylene surfactants: A cavity forming free energy based model

The present paper reports the results of a theoretical study of the forces and factors driving the solubilization of n-alkane solubilizates into the micellar core of some non-ionic surfactants, based on a micellar model which includes the cavity forming free energy as a component of micellization. The solubilizate is n-decane and the non-ionic surfactants considered are n-decyl-polyoxyethylene surfactants. The extent of solubilization, i.e. the mole fraction of the solubilizate within the core has been calculated. The results indicate that the incorporated solubilizate has more translational and rotational degrees of freedom as compared to those of the tail parts of the surfactants present in the core. This drives the total free energy of aggregation after solubilization into a more favourable direction. The results are in fair agreement with the experimental results.

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.

Energy transfer efficiency distributions in polymers in solution during folding and unfolding

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the two ends of a Lennard-Jones polymer chain both at equilibrium and during folding and unfolding has been calculated, for the first time, by Brownian dynamics simulations. The distribution of FRET efficiency becomes bimodal during folding of the extended state subsequent to a temperature quench, with the width of the distribution for the extended state broader than that for the folded state. The reverse process of unfolding subsequent to a upward temperature jump shows different characteristics. The distributions show significant viscosity dependence which can be tested against experiments.

Coupled electro-mechanical response of an electroactive polymer cantilever structure and its application in energy harvesting

Ionic polymer-metal composites (IPMC), piezoelectric polymer composites and nematic elastomer composites are materials, which exhibit characteristics of both sensors and actuators. Large deformation and curvature are observed in these systems when electric potential is applied. Effects of geometric non-linearity due to the chargeinduced motion in these materials are poorly understood. In this paper, a coupled model for understanding the behavior of an ionic polymer beam undergoing large deformation and large curvature is presented. Maxwell's equations and charge transport equations are considered which couple the distribution of the ion concentration and the pressure gradient along length of a cantilever beam with interdigital electrodes. A nonlinear constitutive model is derived accounting for the visco-elasto-plastic behavior of these polymers and based on the hypothesis that the presence of electrical charge stretches/contracts bonds, which give rise to electrical field dependent softening/hardening. Polymer chain orientation in statistical sense plays a role on such softening or hardening. Elementary beam kinematics with large curvature is considered. A model for understanding the deformation due to electrostatic repulsion between asymmetrical charge distributions across the cross-sections is presented. Experimental evidence that Silver(Ag) nanoparticle coated IPMCs can be used for energy harvesting is reported. An IPMC strip is vibrated in different environments and the electric power against a resistive load is measured. The electrical power generated was observed to vary with the environment with maximum power being generated when the strip is in wet state. IPMC based energy harvesting systems have potential applications in tidal wave energy harvesting, residual environmental energy harvesting to power MEMS and NEMS devices.

WEPA: Wind energy potential assessment-spatial decision support system

Spatial Decision Support System (SDSS) assist in strategic decision-making activities considering spatial and temporal variables, which help in Regional planning. WEPA is a SDSS designed for assessment of wind potential spatially. A wind energy system transforms the kinetic energy of the wind into mechanical or electrical energy that can be harnessed for practical use. Wind energy can diversify the economies of rural communities, adding to the tax base and providing new types of income. Wind turbines can add a new source of property value in rural areas that have a hard time attracting new industry. Wind speed is extremely important parameter for assessing the amount of energy a wind turbine can convert to electricity: The energy content of the wind varies with the cube (the third power) of the average wind speed. Estimation of the wind power potential for a site is the most important requirement for selecting a site for the installation of a wind electric generator and evaluating projects in economic terms. It is based on data of the wind frequency distribution at the site, which are collected from a meteorological mast consisting of wind anemometer and a wind vane and spatial parameters (like area available for setting up wind farm, landscape, etc.). The wind resource is governed by the climatology of the region concerned and has large variability with reference to space (spatial expanse) and time (season) at any fixed location. Hence the need to conduct wind resource surveys and spatial analysis constitute vital components in programs for exploiting wind energy. SDSS for assessing wind potential of a region / location is designed with user friendly GUI’s (Graphic User Interface) using VB as front end with MS Access database (backend). Validation and pilot testing of WEPA SDSS has been done with the data collected for 45 locations in Karnataka based on primary data at selected locations and data collected from the meteorological observatories of the India Meteorological Department (IMD). Wind energy and its characteristics have been analysed for these locations to generate user-friendly reports and spatial maps. Energy Pattern Factor (EPF) and Power Densities are computed for sites with hourly wind data. With the knowledge of EPF and mean wind speed, mean power density is computed for the locations with only monthly data. Wind energy conversion systems would be most effective in these locations during May to August. The analyses show that coastal and dry arid zones in Karnataka have good wind potential, which if exploited would help local industries, coconut and areca plantations, and agriculture. Pre-monsoon availability of wind energy would help in irrigating these orchards, making wind energy a desirable alternative.

Activities and their relation to cation distribution in NiAl2O4 MgAl2O4 spinel solid solutions

The activity of NiAl2O4 in NiAl2O4MgAl2O4 solid solutions has been measured by using a solid oxide galvanic cell of the type, Pt, Ni + NiAl2O4 + Al2O3(α)/CaOZrO2/Ni + NixMg1−xAl2O4 + Al2O3(α). Pt, in the temperature range 750–1150°C. The activities in the spinel solid solutions show negative deviations from Raoult's law. The cation distribution in the solid solutions has been calculated using site preference energies independent of composition for Ni2+, Mg2+, and Al3+ ions obtained from crystal field theory and measured cation disorder in pure NiAl2O4 and MgAl2O4, and assumi g ideal mixing of cations on the tetrahedral and octahedral positions. The calculated values correctly predict the decrease in the fraction, α, of Ni2+ ions on tetrahedral sites for 1>x>0.25, observed by Porta et al. [J. Solid State Chem.11, 135 (1974)] but do not support their tentative evidence for an increase in α for x < 0.25. The measured excess free energy of mixing can be completely accounted for by using either the calculated or the measured cation distributions. This suggests that the Madelung energy is approximately a linear function of composition in the solid solutions. The composition of NiOMgO solid solutions in equilibrium with NiAl2O4MgAl2O4 solid solutions has been calculated from the results and information available in literature.

Model projections for household energy use in India

Energy use in developing countries is heterogeneous across households. Present day global energy models are mostly too aggregate to account for this heterogeneity. Here, a bottom-up model for residential energy use that starts from key dynamic concepts on energy use in developing countries is presented and applied to India. Energy use and fuel choice is determined for five end-use functions (cooking, water heating, space heating, lighting and appliances) and for five different income quintiles in rural and urban areas. The paper specifically explores the consequences of different assumptions for income distribution and rural electrification on residential sector energy use and CO(2) emissions, finding that results are clearly sensitive to variations in these parameters. As a result of population and economic growth, total Indian residential energy use is expected to increase by around 65-75% in 2050 compared to 2005, but residential carbon emissions may increase by up to 9-10 times the 2005 level. While a more equal income distribution and rural electrification enhance the transition to commercial fuels and reduce poverty, there is a trade-off in terms of higher CO(2) emissions via increased electricity use. (C) 2011 Elsevier Ltd. All rights reserved.