Structural Transformations, Composition Anomalies and a Dramatic Collapse of Linear Polymer Chains in Dilute Ethanol-Water Mixtures

Water-ethanol mixtures exhibit many interesting anomalies, such as negative excess partial molar volume of ethanol, excess sound absorption coefficient at low concentrations, and positive deviation from Raoult's law for vapor pressure, to mention a few. These anomalies have been attributed to different, often contradictory origins, but a quantitative understanding is still lacking. We show by computer simulation and theoretical analyses that these anomalies arise from the sudden emergence of a bicontinuous phase that occurs at a relatively low ethanol concentration of x(eth) approximate to 0.06-0.10 (that amounts to a volume fraction of 0.17-0.26, which is a significant range!). The bicontinuous phase is formed by aggregation of ethanol molecules, resulting in a weak phase transition whose nature is elucidated. We find that the microheterogeneous structure of the mixture gives rise to a pronounced nonmonotonic composition dependence of local compressibility and nonmonotonic dependence in the peak value of the radial distribution function of ethyl groups. A multidimensional free energy surface of pair association is shown to provide a molecular explanation of the known negative excess partial volume of ethanol in terms of parallel orientation and hence better packing of the ethyl groups in the mixture due to hydrophobic interactions. The energy distribution of the ethanol molecules indicates additional energy decay channels that explain the excess sound attenuation coefficient in aqueous alcohol mixtures. We studied the dependence of the solvation of a linear polymer chain on the composition of the water-ethanol solvent. We find that there is a sudden collapse of the polymer at x(eth) approximate to 0.05-a phenomenon which we attribute to the formation of the microheterogeneous structures in the binary mixture at low ethanol concentrations. Together with recent single molecule pulling experiments, these results provide new insight into the behavior of polymer chain and foreign solutes, such as enzymes, in aqueous binary mixtures.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Study of molecular structure and vibrational spectra of poly (beta,-1-malic acid) PMLA] using DFT approach

Poly (beta-L-malic acid) (PMLA) is a biodegradable polymer and it has various important applications in the biomedical field. In the present work the structural and spectral characteristics of PMLA have been studied by methods of infrared. Raman spectroscopy and quantum chemistry. Electrostatic potential surface, optimized geometry, harmonic vibrational frequencies, infrared intensities and activities of Raman scattering were calculated by density functional theory (DFT) using oligomeric approach employing B3LYP with complete relaxation in the potential energy surface using 6-311++G (d, p) basis set. Based on results, we have discussed the correlation between the vibrational modes and the structure of the PMLA. A complete analysis of the experimental infrared and Raman spectra has been reported on the basis of wavenumber of the vibrational bands and potential energy distribution. The calculated HOMO and LUMO energies shows that charge transfer occur within the molecule. The calculated infrared and the Raman spectra of the polymer based on DFT calculations show reasonable agreement with the experimental results. (c) 2012 Elsevier Ltd. All rights reserved.

LOCAL TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE COLLISIONLESS MAGNETOROTATIONAL INSTABILITY

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the plasma in the disk is essentially collisionless. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell plasma simulations. We focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field, B, until the Alfven velocity, v(A), is comparable to the speed of light, c (independent of the initial value of v(A)/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the pressure perpendicular to B being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume-averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission-from the radio to the gamma-rays-of systems such as Sgr A*.

Intermolecular charge transfer and vibrational analysis of hydrogen bonding in acetazolamide

In the present work the structural and spectral characteristics of acetazolamide have been studied by methods of infrared, Raman spectroscopy and quantum chemistry. Electrostatic potential surface, optimized geometry, harmonic vibrational frequencies, infrared intensities and activities of Raman scattering were calculated by density functional theory (DFT) employing B3LYP with complete relaxation in the potential energy surface using 6-311++G(d,p) basis set. Based on these results, we have discussed the correlation between the vibrational modes and the structure of the dimers of acetazolamide. The calculated vibrational spectra of three dimers of acetazolamide have been compared with observed spectra, and the assignment of observed bands was carried out using potential energy distribution. The observed spectra agree well with the values computed from the OFT. A comparison of observed and calculated vibrational spectra clearly shows the effect of hydrogen bonding. The frequency shifts observed for the different dimers are in accord with the hydrogen bonding in acetazolamide. Natural bond orbital (NBO) analyses reflect the charge transfer interaction in the individual hydrogen bond units and the stability of different dimers of acetazolamide. (C) 2012 Elsevier B.V. All rights reserved.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

Dynamical Transition of Water in the Grooves of DNA Duplex at Low Temperature

At low temperature (below its freezing/melting temperature), liquid water under confinement is known to exhibit anomalous dynamical features. Here we study structure and dynamics of water in the grooves of a long DNA duplex using molecular dynamics simulations with TIP5P potential at low temperature. We find signatures of a dynamical transition in both translational and orientational dynamics of water molecules in both the major and the minor grooves of a DNA duplex. The transition occurs at a slightly higher temperature (TGL ≈ 255 K) than the temperature at which the bulk water is found to undergo a dynamical transition, which for the TIP5P potential is at 247 K. Groove water, however, exhibits markedly different temperature dependence of its properties from the bulk. Entropy calculations reveal that the minor groove water is ordered even at room temperature, and the transition at T ≈ 255 K can be characterized as a strong-to-strong dynamical transition. Confinement of water in the grooves of DNA favors the formation of a low density four-coordinated state (as a consequence of enthalpy−entropy balance) that makes the liquid−liquid transition stronger. The low temperature water is characterized by pronounced tetrahedral order, as manifested in the sharp rise near 109° in the O−O−O angle distribution. We find that the Adams−Gibbs relation between configurational entropy and translational diffusion holds quite well when the two quantities are plotted together in a master plot for different region of aqueous DNA duplex (bulk, major, and minor grooves) at different temperatures. The activation energy for the transfer of water molecules between different regions of DNA is found to be weakly dependent on temperature.

Gibbs Energy of Formation of Cobalt Divanadium Tetroxide

The Gibbs energy of formation of V2O3-saturated spinel CoV2O4 has been measured in the temperature range 900–1700 K using a solid state galvanic cell, which can be represented as Pt, Co + CoV2O4 + V2O3/(CaO) ZrO2/Co + CoO, Pt. The standard free energy of formation of cobalt vanadite from component oxides can be represented as CoO (rs) + V2O3 (cor) → CoV2O4 (sp), ΔG° = −30,125 − 5.06T (± 150) J mole−1. Cation mixing on crystallographically nonequivalent sites of the spinel is responsible for the decrease in free energy with increasing temperature. A correlation between “second law” entropies of formation of cubic 2–3 spinels from component oxides with rock salt and corundum structures and cation distribution is presented. Based on the information obtained in this study and trends in the stability of aluminate and chromite spinels, it can be deduced that copper vanadite is unstable.

The effect of particle size distribution on TG

The effect of a particle size distribution on the fractional reaction has been analysed. The analysis shows that for non-isothermal TG the activation energy and frequency factor evaluated from the fractional reaction by conventional method depend on the particle size distribution, and this may lead to a kinetic compensating effect. Particle size distribution may also lead to an erroneous conclusion about the change in the mechanism of reaction.

Quantitative photoacoustic tomography from boundary pressure measurements: noniterative recovery of optical absorption coefficient from the reconstructed absorbed energy map

We describe a noniterative method for recovering optical absorption coefficient distribution from the absorbed energy map reconstructed using simulated and noisy boundary pressure measurements. The source reconstruction problem is first solved for the absorbed energy map corresponding to single- and multiple-source illuminations from the side of the imaging plane. It is shown that the absorbed energy map and the absorption coefficient distribution, recovered from the single-source illumination with a large variation in photon flux distribution, have signal-to-noise ratios comparable to those of the reconstructed parameters from a more uniform photon density distribution corresponding to multiple-source illuminations. The absorbed energy map is input as absorption coefficient times photon flux in the time-independent diffusion equation (DE) governing photon transport to recover the photon flux in a single step. The recovered photon flux is used to compute the optical absorption coefficient distribution from the absorbed energy map. In the absence of experimental data, we obtain the boundary measurements through Monte Carlo simulations, and we attempt to address the possible limitations of the DE model in the overall reconstruction procedure.

Free energy of mixing of divalent basic oxides with silica

An expression derived for the free energy of mixing of a divalent basic oxide (MO) with SiO2 based on a model of silicate structure, takes into account the distribution of O2- (from MO) into the silica network, the mixing of silicate ions with O2- and the enthalpy of mixing. The resulting expression is ΔGmix=RT{N11n (2N1-N)2/4N1(1-N)+N21n N 2-N/1-N}, where N={(β+N1)-√(β+N 1)2-8βN1N2}/2β β=characteristic constant for the system N1=mol fraction of silica N2=mol fraction of MO. For the proper choice of β, calculated values of the activity of MO for the system PbO-SiO2, MnO-SiO2, FeO-SiO2 and CaO-SiO2 are in good agreement with experiment. The model predicts that the activity of the basic oxide decreases with increase in temperature.

Potential functions for hydrogen bond interactions - IV. Minimum Energy Conformation of the α-Helical Structure of PoIy-L-Alanine

Making use of the empirical potential functions for peptide NH .. O bonds, developed in this laboratory, the relative stabilities of the rightand left-handed α-helical structures of poly-L-alanine have been investigated, by calculating their conformational energies (V). The value of Vmin of the right-handed helix (αP) is about - 10.4 kcal/mole, and that of the left-handed helix (αM) is about - 9.6 kcal/mole, showing that the former is lower in energy by 0.8 kcal/mole. The helical parameters of the stable conformation of αP are n ∼ 3.6 and h ∼ 1.5 Å. The hydrogen bond of length 2.85 Å and nonlinearity of about 10° adds about 4.0 kcal/ mole to the stabilising energy of the helix in the minimum enregy region. The energy minimum is not sharply defined, but occurs over a long valley, suggesting that a distribution of conformations (φ{symbol}, ψ) of nearly the same energy may occur for the individual residues in a helix. The experimental data of a-helical fibres of poly-L-alanine are in good agreement with the theoretical results for αP. In the case of proteins, the mean values of (φ{symbol}, ψ) for different helices are distributed, but they invariably occur within the contour for V = Vmin + 2 kcal/mole for αP.

Energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals

Thermotropic liquid crystals are known to display rich phase behavior on temperature variation. Although the nematic phase is orientationally ordered but translationally disordered, a smectic phase is characterized by the appearance of a partial translational order in addition to a further increase in orientational order. In an attempt to understand the interplay between orientational and translational order in the mesophases that thermotropic liquid crystals typically exhibit upon cooling from the high-temperature isotropic phase, we investigate the potential energy landscapes of a family of model liquid crystalline systems. The configurations of the system corresponding to the local potential energy minima, known as the inherent structures, are determined from computer simulations across the mesophases. We find that the depth of the potential energy minima explored by the system along an isochor grows through the nematic phase as temperature drops in contrast to its insensitivity to temperature in the isotropic and smectic phases. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures; the inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the isotropic-nematic transition. We find that this breakdown occurs at a temperature below which the system explores increasingly deeper potential energy minima.

Transient signal detection using wigner-ville distribution and wavelet denoising

The problem of detecting an unknown transient signal in noise is considered. The SNR of the observed data is first enhanced using wavelet domain filter The output of the wavelet domain filter is then transformed using a Wigner-Ville transform,which separates the spectrum of the observed signal into narrow frequency bands. Each subband signal at the output of the Wigner-ville block is subjected kto wavelet based level dependent denoising (WBLDD)to supress colored noise A weighted sum of the absolute value of outputs of WBLDD is passed through an energy detector, whose output is used as test statistic to take the final decision. By assigning weights proportional to the energy of the corresponding subband signals, the proposed detector approximates a frequency domain matched filter Simulation results are presented to show that the performance of the proposed detector is better than that of the wavelet packet transform based detector.