Communication: Unusual dynamics of hybrid nanoparticles and their binary mixtures

We present the results on the evolution of microscopic dynamics of hybrid nanoparticles and their binary mixtures as a function of temperature and wave vector. We find unexpectedly a nonmonotonic dependence of the structural relaxation time of the nanoparticles as a function of the morphology. In binary mixtures of two of the largest nanoparticles studied, we observe re-entrant vitrification as a function of the volume fraction of the smaller nanoparticle, which is unusual for such high diameter ratio. Possible explanation for the observed behavior is provided. (C) 2010 American Institute of Physics. doi:10.1063/1.3495480]

Inertial effects in solvation dynamics

The dynamics of solvation of newly created charged species in dense dipolar liquids can proceed at a high speed with time constants often in the subpicosecond domain. The motion of the solvent molecules can be in the inertial limit at such short times. In this paper we present a microscopic study of the effects of inertial motion of solvent molecules on the solvation dynamics of a newly created ion in a model dipolar liquid. Interesting dynamical behavior emerges when the relative contribution of the translational modes in the wave-vector-dependent longitudinal relaxation time is significant. Especially, the theory predicts that the time correlation function of the solvation energy can become oscillatory in some limiting situations. In general, the dynamics becomes faster in the presence of the inertial contribution. We discuss the experimental situations where the inertial effects can be noticeable.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Relaxation dominated by inertia: Solvation dynamics of a small ion in a dipolar solvent

It is shown from an analytical theory that the solvation dynamics of a small ion can be controlled largely by the inertial response of the dipolar solvent when the liquid is in the underdamped limit. It is also shown that this inertial response arises primarily from the long wavelength (with wavevector k≃0) processes which have a collective excitation-like behaviour. The long time decay is dominated by the processes occurring at molecular lengthscales. The theoretical results are in good agreement with recent computer simulation results.

Structure of valinomycin by molecular dynamics studies

Valinomycin is an important ionophore which exhibits a high conformational flexibility. The study of various conformations adopted by this molecule together with the study of flexibility in a given conformation can throw light on the ion transport by the ionophore across the membrane. Molecular dynamics (MD) studies are ideal to characterize the flexibility in different parts of the molecule and can also give an idea of various conformations adopted by the molecule at a given temperature. Hence MD studies at 100K have been carried out on the minimized crystal structure of the molecule to scan the possible conformations in the neighbourhood of the well known 'bracelet' like structure of uncomplexed Valinomycin, Properties, like the flexibility, average values, r.m.s. fluctuations of the various intramolecular hydrogen bonds are discussed. Energy minimization has been carried out on selected MD simulated points to analyze the characteristics of the unique conformation adopted by this molecule at this temperature.

Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (T Delta S) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules. (C) 2010 American Institute of Physics. doi:10.1063/1.3494115]

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.

Use of wetland buffer areas to reduce nitrogen transport from forested catchments: Retention capacity, emissions of N2O and CH4 and vegetation composition dynamics

The use of buffer areas in forested catchments has been actively researched during the last 15 years; but until now, the research has mainly concentrated on the reduction of sediment and phosphorus loads, instead of nitrogen (N). The aim of this thesis was to examine the use of wetland buffer areas to reduce the nitrogen transport in forested catchments and to investigate the environmental impacts involved in their use. Besides the retention capacity, particular attention was paid to the main factors contributing to the N retention, the potential for increased N2O emissions after large N loading, the effects of peatland restoration for use as buffer areas on CH4 emissions, as well as the vegetation composition dynamics induced by the use of peatlands as buffer areas. To study the capacity of buffer areas to reduce N transport in forested catchments, we first used large artificial loadings of N, and then studied the capacity of buffer areas to reduce ammonium (NH4-N) export originating from ditch network maintenance areas in forested catchments. The potential for increased N2O emissions were studied using the closed chamber technique and a large artificial N loading at five buffer areas. Sampling for CH4 emissions and methane-cycling microbial populations were done on three restored buffer areas and on three buffers constructed on natural peatlands. Vegetation composition dynamics was studied at three buffer areas between 1996 and 2009. Wetland buffer areas were efficient in retaining inorganic N from inflow. The key factors contributing to the retention were the size and the length of the buffer, the hydrological loading and the rate of nutrient loading. Our results show that although the N2O emissions may increase temporarily to very high levels after a large N loading into the buffer area, the buffer areas in forested catchments should be viewed as insignificant sources of N2O. CH4 fluxes were substantially higher from buffers constructed on natural peatlands than from the restored buffer areas, probably because of the slow recovery of methanogens after restoration. The use of peatlands as buffer areas was followed by clear changes in plant species composition and the largest changes occurred in the upstream parts of the buffer areas and the wet lawn-level surfaces, where the contact between the vegetation and the through-flow waters was closer than for the downstream parts and dry hummock sites. The changes in the plant species composition may be an undesired phenomenon especially in the case of the mires representing endangered mire site types, and therefore the construction of new buffer areas should be primarily directed into drained peatland areas.

Application of an ecological simulation model to Lake Päijänne

Tiivistelmä: Simulointimallin soveltaminen Pohjois-Päijänteellä

Traditionally Protected Forests´ Role within Transforming Natural Resource Management Regimes in Taita Hills, Kenya

In Taita Hills, south-eastern Kenya, remnants of indigenous mountain rainforests play a crucial role as water towers and socio-cultural sites. They are pressurized due to poverty, shortage of cultivable land and the fading of traditional knowledge. This study examines the traditional ecological knowledge of Taitas and the ways it may be applied within transforming natural resource management regimes. I have analyzed some justifications for and hindrances to ethnodevelopment and participatory forest management in light of recently renewed Kenyan forest policies. Mixed methods were applied by combining an ethnographic approach with participatory GIS. I learned about traditionally protected forests and their ecological and cultural status through a seek out the expert method and with remote sensing data and tools. My informants were: 107 household interviewees, 257 focus group participants, 73 key informants and 87 common informants in participatory mapping. Religious leaders and state officials shared their knowledge for this study. I have gained a better understanding of the traditionally protected forests and sites through examining their ecological characteristics and relation to social dynamics, by evaluating their strengths and hindrances as sites for conservation of cultural and biological diversity. My results show that, these sites are important components of a complex socio-ecological system, which has symbolical status and sacred and mystical elements within it, that contributes to the connectivity of remnant forests in the agroforestry dominated landscape. Altogether, 255 plant species and 220 uses were recognized by the tradition experts, whereas 161 species with 108 beneficial uses were listed by farmers. Out of the traditionally protected forests studied 47 % were on private land and 23% on community land, leaving 9% within state forest reserves. A paradigm shift in conservation is needed; the conservation area approach is not functional for private lands or areas trusted upon communities. The role of traditionally protected forests in community-based forest management is, however, paradoxal, since communal approaches suggests equal participation of people, whereas management of these sites has traditionally been the duty of solely accredited experts in the village. As modernization has gathered pace such experts have become fewer. Sacredness clearly contributes but, it does not equal conservation. Various social, political and economic arrangements further affect the integrity of traditionally protected forests and sites, control of witchcraft being one of them. My results suggest that the Taita have a rich traditional ecological knowledge base, which should be more determinately integrated into the natural resource management planning processes.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

A molecular dynamics study of a PbO.SiO2 glass and melt

The structures of a PbO.SiO2 glass and melt have been studied using molecular dynamics simulation employing Born-Mayer-Huggins pair potentials. Various pair distribution functions are presented and discussed. Pb-Pb correlations persist in the melt, in agreement with experimental observations. The calculated and experimental radial distribution functions are compared.