Molecular assemblies for solar energy conversion - Biomimetic approach

Excitation energy migration followed by electron transfer forms the main components of natural photosynthesis. An understanding of these aspects is essential to reenact the primary processes in laboratory under conditions that are precisely repeatable. Here we describe the state of understanding of the natural process and several model systems designed to harvest solar energy and conversion to useful form of chemical energy. The molecular assemblies constituting the model systems offer a great advantage in terms of better comprehension of the mechanistic aspects and yield valuable information on the design of molecular photonic devices.

Probing potential energy surfaces in confined systems: behavior of mean-square displacement in zeolites

Time evolution of mean-squared displacement based on molecular dynamics for a variety of adsorbate-zeolite systems is reported. Transition from ballistic to diffusive behavior is observed for all the systems. The transition times are found to be system dependent and show different types of dependence on temperature. Model calculations on a one-dimensional system are carried out which show that the characteristic length and transition times are dependent on the distance between the barriers, their heights, and temperature. In light of these findings, it is shown that it is possible to obtain valuable information about the average potential energy surface sampled under specific external conditions.

Ricin-membrane interaction: membrane penetration depth by fluorescence quenching and resonance energy transfer

The entry of the plant toxin ricin and its A- and B-subunits in model membranes in the presence as well as absence of monosialoganglioside (GM(1)) has been studied. Dioleoylphosphatidylcholine and 5-, 10-, and 12-doxyl- or 9,10-dibromophosphatidylcholines serve as quenchers of intrinsic tryptophan fluorescence of the proteins. The parallax method of Chattopadhyay and London [(1987) Biochemistry 26, 39-45] has been employed to measure the average membrane penetration depth of tryptophans of ricin and its B-chain and the actual depth of the sole Trp 211 in the A-chain. The results indicate that both of the chains as well as intact ricin penetrate the membrane deeply and the C-terminal end of the A-chain is well inside the bilayer, especially at pH 4.5. An extrinsic probe N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl) ethylenediamine (I-AEDANS) has been attached to Cys 259 of the A-chain, and the kinetics of penetration has been followed by monitoring the increase in AEDANS fluorescence at 480 nm. The insertion follows first-order kinetics, and the rate constant is higher at a lower pH. The energy transfer distance analysis between Trp 211 and AEDANS points out that the conformation of the A-chain changes as it inserts into the membrane. CD studies indicate that the helicity of the proteins increases after penetration, which implies that some of the unordered structure in the native protein is converted to the ordered form during this process. Hydrophobic forces seem to be responsible for stabilizing a particular protein conformation inside the membrane.

Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Monte-Carlo estimation of time-dependent statistical characteristics of random dynamical systems

The problem of estimating the time-dependent statistical characteristics of a random dynamical system is studied under two different settings. In the first, the system dynamics is governed by a differential equation parameterized by a random parameter, while in the second, this is governed by a differential equation with an underlying parameter sequence characterized by a continuous time Markov chain. We propose, for the first time in the literature, stochastic approximation algorithms for estimating various time-dependent process characteristics of the system. In particular, we provide efficient estimators for quantities such as the mean, variance and distribution of the process at any given time as well as the joint distribution and the autocorrelation coefficient at different times. A novel aspect of our approach is that we assume that information on the parameter model (i.e., its distribution in the first case and transition probabilities of the Markov chain in the second) is not available in either case. This is unlike most other work in the literature that assumes availability of such information. Also, most of the prior work in the literature is geared towards analyzing the steady-state system behavior of the random dynamical system while our focus is on analyzing the time-dependent statistical characteristics which are in general difficult to obtain. We prove the almost sure convergence of our stochastic approximation scheme in each case to the true value of the quantity being estimated. We provide a general class of strongly consistent estimators for the aforementioned statistical quantities with regular sample average estimators being a specific instance of these. We also present an application of the proposed scheme on a widely used model in population biology. Numerical experiments in this framework show that the time-dependent process characteristics as obtained using our algorithm in each case exhibit excellent agreement with exact results. (C) 2010 Elsevier Inc. All rights reserved.

Effect of dynamical asymmetry on the viscosity of a random copolymer melt

The variation of the viscosity as a function of the sequence distribution in an A-B random copolymer melt is determined. The parameters that characterize the random copolymer are the fraction of A monomers f, the parameter lambda which determines the correlation in the monomer identities along a chain and the Flory chi parameter chi(F) which determines the strength of the enthalpic repulsion between monomers of type A and B. For lambda>0, there is a greater probability of finding like monomers at adjacent positions along the chain, and for lambda<0 unlike monomers are more likely to be adjacent to each other. The traditional Markov model for the random copolymer melt is altered to remove ultraviolet divergences in the equations for the renormalized viscosity, and the phase diagram for the modified model has a binary fluid type transition for lambda>0 and does not exhibit a phase transition for lambda<0. A mode coupling analysis is used to determine the renormalization of the viscosity due to the dependence of the bare viscosity on the local concentration field. Due to the dissipative nature of the coupling. there are nonlinearities both in the transport equation and in the noise correlation. The concentration dependence of the transport coefficient presents additional difficulties in the formulation due to the Ito-Stratonovich dilemma, and there is some ambiguity about the choice of the concentration to be used while calculating the noise correlation. In the Appendix, it is shown using a diagrammatic perturbation analysis that the Ito prescription for the calculation of the transport coefficient, when coupled with a causal discretization scheme, provides a consistent formulation that satisfies stationarity and the fluctuation dissipation theorem. This functional integral formalism is used in the present analysis, and consistency is verified for the present problem as well. The upper critical dimension for this type of renormaliaation is 2, and so there is no divergence in the viscosity in the vicinity of a critical point. The results indicate that there is a systematic dependence of the viscosity on lambda and chi(F). The fluctuations tend to increase the viscosity for lambda<0, and decrease the viscosity for lambda>0, and an increase in chi(F) tends to decrease the viscosity. (C) 1996 American Institute of Physics.

Investigation of the characteristics of an SF6 rotating arc by a mathematical model

A mathematical model has been developed for predicting the performance of rotating arcs in SF6 gas by considering the energy balance and force balance equations. The finite difference technique has been adopted for the computer simulation of the arc characteristics. This method helps in considering the spatial variation of the transport and radiative properties of the arc. All the three heat loss mechanisms-conduction, convection, and radiation-have been considered. Results obtained over a 10 ms (half cycle of 50 Hz wave) current flow period for 1.4 kA (peak) and 4.2 kA (peak), show that the proposed arc model gives the expected behavior of the arc over the range of currents studied.

Model Hamiltonians for nonlinear optical properties of conjugated polymers

Quantum cell models for delocalized electrons provide a unified approach to the large NLO responses of conjugated polymers and pi-pi* spectra of conjugated molecules. We discuss exact NLO coefficients of infinite chains with noninteracting pi-electrons and finite chains with molecular Coulomb interactions V(R) in order to compare exact and self-consistent-field results, to follow the evolution from molecular to polymeric responses, and to model vibronic contributions in third-harmonic-generation spectra. We relate polymer fluorescence to the alternation delta of transfer integrals t(1+/-delta) along the chain and discuss correlated excited states and energy thresholds of conjugated polymers.

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.

The structures of bacteriorhodopsin with different retinal-Schiff base orientations--computer modeling and energy minimization studies

Bacteriorhodopsin has been the subject of intense study in order to understand its photochemical function. The recent atomic model proposed by Henderson and coworkers based on electron cryo-microscopic studies has helped in understanding many of the structural and functional aspects of bacteriorhodopsin. However, the accuracy of the positions of the side chains is not very high since the model is based on low-resolution data. In this study, we have minimized the energy of this structure of bacteriorhodopsin and analyzed various types of interactions such as - intrahelical and interhelical hydrogen bonds and retinal environment. In order to understand the photochemical action, it is necessary to obtain information on the structures adopted at the intermediate states. In this direction, we have generated some intermediate structures taking into account certain experimental data, by computer modeling studies. Various isomers of retinal with 13-cis and/or 15-cis conformations and all possible staggered orientations of Lys-216 side chain were generated. The resultant structures were examined for the distance between Lys-216-schiff base nitrogen and the carboxylate oxygen atoms of Asp-96 - a residue which is known to reprotonate the schiff base at later stages of photocycle. Some of the structures were selected on the basis of suitable retinal orientation and the stability of these structures were tested by energy minimization studies. Further, the minimized structures are analyzed for the hydrogen bond interactions and retinal environment and the results are compared with those of the minimized rest state structure. The importance of functional groups in stabilizing the structure of bacteriorhodopsin and in participating dynamically during the photocycle have been discussed.

Structured systems methodology for evaluation of random interruptions in continuous process type manufacturing systems

A structured systems methodology was developed to analyse the problems of production interruptions occurring at random intervals in continuous process type manufacturing systems. At a macro level the methodology focuses on identifying suitable investment policies to reduce interruptions of a total manufacturing system that is a combination of several process plants. An interruption-tree-based simulation model was developed for macroanalysis. At a micro level the methodology focuses on finding the effects of alternative configurations of individual process plants on the overall system performance. A Markov simulation model was developed for microlevel analysis. The methodology was tested with an industry-specific application.

The energy balance in the Ramdas layer

On calm clear nights, air at a height of a few decimetres above bare soil can be cooler than the surface by several degrees in what we shall call the Ramdas layer (Ramdas and Atmanathan, 1932). The authors have recently offered a logical explanation for such a lifted temperature minimum, together with a detailed numerical model. In this paper, we provide physical insight into the phenomenon by a detailed discussion of the energy budget in four typical cases, including one with a lifted minimum. It is shown that the net cooling rate near ground is the small difference between two dominant terms, representing respectively radiative upflux from the ground and from the air layers just above ground. The delicate energy balance that leads to the lifted minimum is upset by turbulent transport, by surface emissivity approaching unity, or by high ground cooling rates. The rapid variation of the flux emissivity of humid air is shown to dominate radiative transport near the ground.

Nonequilibrium phase transition in the kinetic Ising model: Divergences of fluctuations and responses near the transition point

The nonequilibrium dynamic phase transition in the kinetic Ising model in the presence of an oscillating magnetic field is studied by Monte Carlo simulation. The fluctuation of the dynamic older parameter is studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined ''susceptibility'' is also studied near the dynamic transition point. Similarly, the fluctuation of energy and appropriately defined ''specific heat'' is studied as a function of temperature near the dynamic transition point. In both cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (in power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).

Decentralized Sustainable Energy Planning of Tumkur District, India

Present work shows the feasibility of decentralized energy options for the Tumkur district in India. Decentralized energy planning (DEP) involves scaling down energy planning to subnational or regional scales. The important aspect of the energy planning at decentralized level would be to prepare an area-based DEP to meet energy needs and development of alternate energy sources at least-cost to the economy and environment. The geographical coverage and scale reflects the level at which the analysis takes place, which is an important factor in determining the structure of models. In the present work, DEP modeling under different scenarios has been carried out for Tumkur district of India for the year 2020. DEP model is suitably scaled for obtaining the optimal mix of energy resources and technologies using a computer-based goal programming technique. The rural areas of the Tumkur district have different energy needs. Results show that electricity needs can be met by biomass gasifier technology, using biomass feedstock produced by allocating only 12% of the wasteland in the district at 8 t/ha/yr of biomass productivity. Surplus electricity can be produced by adopting the option of biomass power generation from energy plantations. The surplus electricity generated can be supplied to the grid. The sustainable development scenario is a least cost scenario apart from promoting self-reliance, local employment, and environmental benefits. (C) 2010 American Institute of Chemical Engineers Environ Prog, 30: 248-258, 2011

Phase separation in a simple model with dynamical asymmetry

We perform computer simulations of a Cahn-Hilliard model of phase separation that has dynamical asymmetry between the two coexisting phases. The dynamical asymmetry is incorporated by considering a mobility function that is order parameter dependent. Simulations of this model reveal morphological features similar to those observed in viscoelastic phase separation. In the early stages, the minority phase domains form a percolating structure that shrinks with time, eventually leading to the formation of disconnected regions that are characterized by the presence of random interfaces as well as isolated droplets. The domains grow as L(t)similar to t(1/3) in the very late stages. Although dynamical scaling is violated in the area shrinking regime, it is restored at late times. However, the form of the scaling function is found to depend on the extent of dynamical asymmetry. [S1063-651X(99)12101-9].