Critical micelle concentration from a lattice gas model

The temperature dependence of the critical micelle concentration (CMC) and a closed-loop coexistence curve are obtained, via Monte Carlo simulations, in the water surfactant limit of a two-dimensional version of a statistical mechanical model for micro-emulsions, The CMC and the coexistence curve reproduce various experimental trends as functions of the couplings. In the oil-surfactant limit, there is a conventional coexistence cure with an upper consolute point that allows for a region of three-phase coexistence between oil-rich, water-rich and microemulsion phases.

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some of the current issues of the solar simulations are put forward.

Estimating Frequencies of Minority Nevirapine-Resistant Strains in Chronically HIV-1-Infected Individuals Naive to Nevirapine by Using Stochastic Simulations and a Mathematical Model

Nevirapine forms the mainstay of our efforts to curtail the pediatric AIDS epidemic through prevention of mother-to-child transmission of HIV-1. A key limitation, however, is the rapid selection of HIV-1 strains resistant to nevirapine following the administration of a single dose. This rapid selection of resistance suggests that nevirapine-resistant strains preexist in HIV-1 patients and may adversely affect outcomes of treatment. The frequencies of nevirapine-resistant strains in vivo, however, remain poorly estimated, possibly because they exist as a minority below current assay detection limits. Here, we employ stochastic simulations and a mathematical model to estimate the frequencies of strains carrying different combinations of the common nevirapine resistance mutations K103N, V106A, Y181C, Y188C, and G190A in chronically infected HIV-1 patients naive to nevirapine. We estimate the relative fitness of mutant strains from an independent analysis of previous competitive growth assays. We predict that single mutants are likely to preexist in patients at frequencies (similar to 0.01% to 0.001%) near or below current assay detection limits (>0.01%), emphasizing the need for more-sensitive assays. The existence of double mutants is subject to large stochastic variations. Triple and higher mutants are predicted not to exist. Our estimates are robust to variations in the recombination rate, cellular superinfection frequency, and the effective population size. Thus, with 10(7) to 10(8) infected cells in HIV-1 patients, even when undetected, nevirapine-resistant genomes may exist in substantial numbers and compromise efforts to prevent mother-to-child transmission of HIV-1, accelerate the failure of subsequent antiretroviral treatments, and facilitate the transmission of drug resistance.

The Mechanics and Statistics of Active Matter

Active particles contain internal degrees of freedom with the ability to take in and dissipate energy and, in the process, execute systematic movement. Examples include all living organisms and their motile constituents such as molecular motors. This article reviews recent progress in applying the principles of nonequilibrium statistical mechanics and hydrodynamics to form a systematic theory of the behavior of collections of active particles-active matter-with only minimal regard to microscopic details. A unified view of the many kinds of active matter is presented, encompassing not only living systems but inanimate analogs. Theory and experiment are discussed side by side.

Disordered matter studied by x-ray Raman scattering

Spectroscopy can provide valuable information on the structure of disordered matter beyond that which is available through e.g. x-ray and neutron diffraction. X-ray Raman scattering is a non-resonant element-sensitive process which allows bulk-sensitive measurements of core-excited spectra from light-element samples. In this thesis, x-ray Raman scattering is used to study the local structure of hydrogen-bonded liquids and solids, including liquid water, a series of linear and branched alcohols, and high-pressure ice phases. Connecting the spectral features to the local atomic-scale structure involves theoretical references, and in the case of hydrogen-bonded systems the interpretation of the spectra is currently actively debated. The systematic studies of the intra- and intermolecular effects in alcohols, non-hydrogen-bonded neighbors in high-pressure ices, and the effect of temperature in liquid water are used to demonstrate different aspects of the local structure that can influence the near-edge spectra. Additionally, the determination of the extended x-ray absorption fine structure is addressed in a momentum-transfer dependent study. This work demonstrates the potential of x-ray Raman scattering for unique studies of the local structure of a variety of disordered light-element systems.

Decomposition of Scots pine fine woody debris in boreal conditions

"Litter quality and environmental effects on Scots pine (Pinus sylvestris L.) fine woody debris (FWD) decomposition were examined in three forestry-drained peatlands representing different site types along a climatic gradient from the north boreal (Northern Finland) to south (Southern Finland) and hemiboreal (Central Estonia) conditions. Decomposition (percent mass loss) of FWD with diameter <= 10 mm (twigs) and FWD with diameter > 10 mm (branches) was measured using the litter bag method over 1-4-year periods. Overall, decomposition rates increased from north to south, the rate constants (k values) varying from 0.128 to 0.188 year(-1) and from 0.066 to 0.127 year(-1) for twigs and branches, respectively. On average, twigs had lost 34%, 19% and 19%, and branches 25%, 17% and 11% of their initial mass after 2 years of decomposition at the hemiboreal, south boreal and north boreal sites, respectively. After 4 years at the south boreal site the values were 48% for twigs and 42% for branches. Based on earlier studies, we suggest that the decomposition rates that we determined may be used for estimating Scots pine FWD decomposition in the boreal zone, also in upland forests. Explanatory models accounted for 50.4% and 71.2% of the total variation in FWD decomposition rates when the first two and all years were considered, respectively. The variables most related to FWD decomposition included the initial ash, water extractives and Klason lignin content of litter, and cumulative site precipitation minus potential evapotranspiration. Simulations of inputs and decomposition of Scots pine FWD and needle litter in south boreal conditions over a 60-year period showed that 72 g m(-2) of organic matter from FWD vs. 365 g m(-2) from needles accumulated in the forest floor. The annual inputs varied from 5.7 to 15.6 g m(-2) and from 92 to 152 g m(-2) for FWD and needles, respectively. Each thinning caused an increase in FWD inputs, Up to 510 g m(-2), while the needle inputs did not change dramatically. Because the annual FWD inputs were lowered following the thinnings, the overall effect of thinnings on C accumulation from FWD was slightly negative. The contribution of FWD to soil C accumulation, relative to needle litter, seems to be rather minor in boreal Scots pine forests. (C) 2008 Elsevier B.V. All rights reserved."

Redundancy resolution using tractrix -- Simulations and Experiments

Hyper-redundant robots are characterized by the presence of a large number of actuated joints, many more than the number required to perform a given task. These robots have been proposed and used for many applications involving avoiding obstacles or, in general, to provide enhanced dexterity in performing tasks. Making effective use of the extra degrees of freedom or resolution of redundancy has been an extensive topic of research and several methods have been proposed in literature. In this paper, we compare three known methods and show that an algorithm based on a classical curve called the tractrix leads to a more 'natural' motion of the hyper-redundant robot, with the displacements diminishing from the end-effector to the fixed base. In addition, since the actuators nearer the base 'see' a greater inertia due to the links farther away, smaller motion of the actuators nearer the base results in better motion of the end-effector as compared to other two approaches. We present simulation and experimental results performed on a prototype eight link planar hyper-redundant manipulator.

Large-amplitude chirped coherent phonons in tellurium mediated by ultrafast photoexcited carrier diffusion

We report femtosecond time-resolved reflectivity measurements of coherent phonons in tellurium performed over a wide range of temperatures (3-296 K) and pump-laser intensities. A totally symmetric A(1) coherent phonon at 3.6 THz responsible for the oscillations in the reflectivity data is observed to be strongly positively chirped (i.e., phonon time period decreases at longer pump-probe delay times) with increasing photoexcited carrier density, more so at lower temperatures. We show that the temperature dependence of the coherent phonon frequency is anomalous (i.e, increasing with increasing temperature) at high photoexcited carrier density due to electron-phonon interaction. At the highest photoexcited carrier density of (1.4 x 10(21) cm(-3) and the sample temperature of 3 K, the lattice displacement of the coherent phonon mode is estimated to be as high as similar to 0.24 angstrom. Numerical simulations based on coupled effects of optical absorption and carrier diffusion reveal that the diffusion of carriers dominates the nonoscillatory electronic part of the time-resolved reflectivity. Finally, using the pump-probe experiments at low carrier density of 6 x 10(18) cm(-3), we separate the phonon anharmonicity to obtain the electron-phonon coupling contribution to the phonon frequency and linewidth.

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

Soft matter electrolytes based on polymethylmetacrylate dispersions in lithium bis(trifluoromethanesulfonyl)imide/1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ionic liquids

Ion transport in a polymer-ionic liquid (IL) soft matter composite electrolyte is discussed here in detail in the context of polymer-ionic liquid interaction and glass transition temperature The dispersion of polymethylmetacrylate (PMMA) in 1-butyl-3-methylimidazolium hexafluorophosphate (BMIPF6) and 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide (BMITFSI) resulted in transparent composite electrolytes with a jelly-like consistency The composite ionic conductivity measured over the range -30 C to 60 C was always lower than that of the neat BMITFSI/BMIPF6 and LiTFSI-BMITFSI/LiTFSI-BMIPF6 electrolytes but still very high (>1 mS/cm at 25 degrees C up to 50 wt% PMMA) While addition of LiTFSI to IL does not influence the glass T-g and T-m melting temperature significantly dispersion of PMMA (especially at higher contents) resulted in increase in T-g and disappearance of T-m In general the profile of temperature-dependent ionic conductivity could be fitted to Vogel-Tamman-Fulcher (VTF) suggesting a solvent assisted ion transport However for higher PMMA concentration sharp demarcation of temperature regimes between thermally activated and solvent assisted ion transport were observed with the glass transition temperature acting as the reference point for transformation from one form of transport mechanism to the other Because of the beneficial physico-chemical properties and interesting ion transport mechanism we envisage the present soft matter electrolytes to be promising for application in electrochromic devices (C) 2010 Elsevier Ltd All rights reserved

Cell-dynamical simulation of magnetic hysteresis in the two-dimensional Ising system

We present results from numerical simulations using a ‘‘cell-dynamical system’’ to obtain solutions to the time-dependent Ginzburg-Landau equation for a scalar, two-dimensional (2D), (Φ2)2 model in the presence of a sinusoidal external magnetic field. Our results confirm a recent scaling law proposed by Rao, Krishnamurthy, and Pandit [Phys. Rev. B 42, 856 (1990)], and are also in excellent agreement with recent Monte Carlo simulations of hysteretic behavior of 2D Ising spins by Lo and Pelcovits [Phys. Rev. A 42, 7471 (1990)].

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.