Non-monotonic, distance-dependent relaxation of water in reverse micelles: Propagation of surface induced frustration along hydrogen bond networks

Layer-wise, distance-dependent orientational relaxation of water confined in reverse micelles (RM) is studied using theoretical and computational tools. We use both a newly constructed ``spins on a ring'' (SOR) Ising-type model (with Shore-Zwanzig rotational dynamics) and atomistic simulations with explicit water. Our study explores the effect of reverse micelle size and role of intermolecular correlations, compromised by the presence of a highly polar surface, on the distance (from the interface) dependence of water relaxation. The ``spins on a ring'' model can capture some aspects of distance dependence of relaxation, such as acceleration of orientational relaxation at intermediate layers. In atomistic simulations, layer-wise decomposition of hydrogen bond formation pattern clearly reveals that hydrogen bond arrangement of water at a certain distance away from the surface can remain frustrated due to the interaction with the polar surface head groups. This layer-wise analysis also reveals the presence of a non-monotonic slow relaxation component which can be attributed to this frustration effect and which is accentuated in small to intermediate size RMs. For large size RMs, the long time component decreases monotonically from the interface to the interior of the RMs with slowest relaxation observed at the interface. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4732095]

Dynamics of the Indian-Ocean oxygen minimum zones

In the Indian Ocean, mid-depth oxygen minimum zones (OMZs) occur in the Arabian Sea and the Bay of Bengal. The lower part of the Arabian-Sea OMZ (ASOMZ; below 400 m) intensifies northward across the basin; in contrast, its upper part (above 400 m) is located in the central/eastern basin, well east of the most productive regions along the western boundary. The Bay-of-Bengal OMZ (BBOMZ), although strong, is weaker than the ASOMZ. To investigate the processes that maintain the Indian-Ocean OMZs, we obtain a suite of solutions to a coupled biological/physical model. Its physical component is a variable-density, 6 1/2-layer model, in which each layer corresponds to a distinct dynamical regime or water-mass type. Its biological component has six compartments: nutrients, phytoplankton, zooplankton, two size classes of detritus, and oxygen. Because the model grid is non-eddy resolving (0.5 degrees), the biological model also includes a parameterization of enhanced mixing based on the eddy kinetic energy derived from satellite observations. To explore further the impact of local processes on OMZs, we also obtain analytic solutions to a one-dimensional, simplified version of the biological model. Our control run is able to simulate basic features of the oxygen, nutrient, and phytoplankton fields throughout the Indian Ocean. The model OMZs result from a balance, or lack thereof, between a sink of oxygen by remineralization and subsurface oxygen sources due primarily to northward spreading of oxygenated water from the Southern Hemisphere, with a contribution from Persian-Gulf water in the northern Arabian Sea. The northward intensification of the lower ASOMZ results mostly from horizontal mixing since advection is weak in its depth range. The eastward shift of the upper ASOMZ is due primarily to enhanced advection and vertical eddy mixing in the western Arabian Sea, which spread oxygenated waters both horizontally and vertically. Advection carries small detritus from the western boundary into the central/eastern Arabian Sea, where it provides an additional source of remineralization that drives the ASOMZ to suboxic levels. The model BBOMZ is weaker than the ASOMZ because the Bay lacks a remote source of detritus from the western boundary. Although detritus has a prominent annual cycle, the model OMZs do not because there is not enough time for significant remineralization to occur.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.

Impact dynamics of high Weber number drops on chemically modified metallic surfaces

Hydrophobic/superhydrophobic metallic surfaces prepared via chemical treatment are encountered in many industrial scenarios involving the impingement of spray droplets. The effectiveness of such surfaces is understood through the analysis of droplet impact experiments. In the present study, three target surfaces with aluminum (Al-6061) as base material-acid-etched, Octadecyl Trichloro Silane (OTS) coated, and acid-etched plus OTS-coated-were prepared. Experiments on the impact of inertia dominated water drops on these chemically modified aluminum surfaces were carried out with the objective to highlight the effect of chemical treatment on the target surfaces on key sub-processes occurring in drop impact phenomenon. High speed videos of the entire drop impact dynamics were captured at three Weber number (We) conditions representative of high We (We > 200) regime. During the early stages of drop spreading, the drop impact resulted in ejection of secondary droplets from spreading drop front on the etched surfaces resembling prompt splash on rough surfaces whereas no such splashing was observable on untreated aluminum surface. Prominent development of undulations (fingers) were observed at the rim of drop spreading on the etched surfaces; between the etched surfaces the OTS-coated surface showed a subdued development of fingers than the uncoated surface. The impacted drops showed intense receding on OTS-coated surfaces whereas on the etched surface a highly irregular receding, with drop liquid sticking to the surface, was observed. Quantitative analyses were performed to reveal the effect of target surface characteristics on drop impact parameters such as temporal variation of spread factor of drop lamella, temporal variation of average finger length during spreading phase, maximum drop spreading, time taken to attain maximum spreading, sensitivity of maximum spreading to We, number of fingers at maximum spreading, and average receding velocity of drop lamella. Existing models for maximum drop spreading showed reasonably good agreement with the experimental measurements on the target surfaces except the acid-etched surface. (C) 2014 Elsevier B.V. All rights reserved.

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Driving force of water entry into hydrophobic channels of carbon nanotubes: entropy or energy?

Spontaneous entry of water molecules inside single-wall carbon nanotubes (SWCNTs) has been confirmed by both simulations and experiments. Using molecular dynamics simulations, we have studied the thermodynamics of filling of a (6,6) carbon nanotube in a temperature range from 273 to 353K and with different strengths of the nanotube-water interaction. From explicit energy and entropy calculations using the two-phase thermodynamics method, we have presented a thermodynamic understanding of the filling behaviour of a nanotube. We show that both the energy and the entropy of transfer decrease with increasing temperature. On the other hand, scaling down the attractive part of the carbon-oxygen interaction results in increased energy of transfer while the entropy of transfer increases slowly with decreasing the interaction strength. Our results indicate that both energy and entropy favour water entry into (6,6) SWCNTs. Our results are compared with those of several recent studies of water entry into carbon nanotubes.

Cell Dynamics of Model Proteins

We design rapidly folding sequences by assigning the strongest couplings to the contacts present in a target native state in a two dimensional model of heteropolymers. The pathways to folding and their dependence on the temperature are illustrated via a mapping of the dynamics into motion within the space of the maximally compact cells.

Dynamics of ribonuclease A and ribonuclease S: Computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNaseS and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model, According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Dynamics of a dilute sheared inelastic fluid. II. The effect of correlations

The effect of correlations on the viscosity of a dilute sheared inelastic fluid is analyzed using the ring-kinetic equation for the two-particle correlation function. The leading-order contribution to the stress in an expansion in epsilon=(1-e)(1/2) is calculated, and it is shown that the leading-order viscosity is identical to that obtained from the Green-Kubo formula, provided the stress autocorrelation function in a sheared steady state is used in the Green-Kubo formula. A systemmatic extension of this to higher orders is also formulated, and the higher-order contributions to the stress from the ring-kinetic equation are determined in terms of the terms in the Chapman-Enskog solution for the Boltzmann equation. The series is resummed analytically to obtain a renormalized stress equation. The most dominant contributions to the two-particle correlation function are products of the eigenvectors of the conserved hydrodynamic modes of the two correlated particles. In Part I, it was shown that the long-time tails of the velocity autocorrelation function are not present in a sheared fluid. Using those results, we show that correlations do not cause a divergence in the transport coefficients; the viscosity is not divergent in two dimensions, and the Burnett coefficients are not divergent in three dimensions. The equations for three-particle and higher correlations are analyzed diagrammatically. It is found that the contributions due to the three-particle and higher correlation functions to the renormalized viscosity are smaller than those due to the two-particle distribution function in the limit epsilon -> 0. This implies that the most dominant correlation effects are due to the two-particle correlations.

Anomalous thermal dynamics of Bragg gratings inscribed in germanosilicate optical fiber

An interesting, periodic appearance of a new peak has been observed in the reflected spectrum of a Fiber Bragg Grating (FBG) inscribed in a germanosilicate fiber during thermal treatment. The new peak occurs on the longer wavelength side of the spectrum during heating and on the shorter wavelength side during cooling, following an identical reverse dynamics. Comparison with a commercial grating with 99.9% reflectivity shows a similar decay dynamics. It is proposed that the distortion due to simultaneous erasure and thermal expansion of the index modulation profile may be responsible for the observed anomaly. The reported results help us in understanding the thermal behavior of FBGs and provide additional insights into the mechanisms responsible for the photosensitivity in germanosilicate fibers.

The dynamics of polymerized carbon nanotubes in semiconductor polymer electronics and electro-mechanical sensing

Polymerized carbon nanotubes (CNTs) are promising materials for polymer-based electronics and electro-mechanical sensors. The advantage of having a polymer nanolayer on CNTs widens the scope for functionalizing it in various ways for polymer electronic devices. However, in this paper, we show for the first time experimentally that, due to a resistive polymer layer having carbon nanoparticle inclusions and polymerized carbon nanotubes, an interesting dynamics can be exploited. We first show analytically that the relative change in the resistance of a single isolated semiconductive nanotube is directly proportional to the axial and torsional dynamic strains, when the strains are small, whereas, in polymerized CNTs, the viscoelasticity of the polymer and its effective electrical polarization give rise to nonlinear effects as a function of frequency and bias voltage. A simplified formula is derived to account for these effects and validated in the light of experimental results. CNT–polymer-based channels have been fabricated on a PZT substrate. Strain sensing performance of such a one-dimensional channel structure is reported. For a single frequency modulated sine pulse as input, which is common in elastic and acoustic wave-based diagnostics, imaging, microwave devices, energy harvesting, etc, the performance of the fabricated channel has been found to be promising.

Programming an FPGA to Emulate the Dynamics of DC Machines

This paper presents the programming an FPGA (Field Programmable Gate Array) to emulate the dynamics of DC machines. FPGA allows high speed real time simulation with high precision. The described design includes block diagram representation of DC machine, which contain all arithmetic and logical operations. The real time simulation of the machine in FPGA is controlled by user interfaces they are Keypad interface, LCD display on-line and digital to analog converter. This approach provides emulation of electrical machine by changing the parameters. Separately Exited DC machine implemented and experimental results are presented.

Dynamics of the ammonium ion in the ferroelectric transition in sodium ammonium selenate, NaNH4SeO4·2H2O-an E.S.R. study

E.S.R. investigations of γ-irradiated ferroelectric Sodium ammonium selenate, NaNH4SeO4•2H2O and its deuteriated analogue in powder and single crystal forms have led to a deeper understanding of the nature of the ferroelectric transition of 180 K. A number of paramagnetic species formed due to γ-irradiation have been identified on the basis of their g-factors and hyperfine features from 77Se. The radical SeO4 has been used as a microprobe in studying the phase transition.

Classical Langevin dynamics of a charged particle moving on a sphere and diamagnetism: A surprise

It is generally known that the orbital diamagnetism of a classical system of charged particles in thermal equilibrium is identically zero —the Bohr-van Leeuwen theorem. Physically, this null result derives from the exact cancellation of the orbital diamagnetic moment associated with the complete cyclotron orbits of the charged particles by the paramagnetic moment subtended by the incomplete orbits skipping the boundary in the opposite sense. Motivated by this crucial but subtle role of the boundary, we have simulated here the case of a finite but unbounded system, namely that of a charged particle moving on the surface of a sphere in the presence of an externally applied uniform magnetic field. Following a real space-time approach based on the classical Langevin equation, we have computed the orbital magnetic moment that now indeed turns out to be non-zero and has the diamagnetic sign. To the best of our knowledge, this is the first report of the possibility of finite classical diamagnetism in principle, and it is due to the avoided cancellation.