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.

GINZBURG-LANDAU LIKE THEORY OF HIGH TEMPERATURE SUPERCONDUCTIVITY IN THE CUPRATES: EMERGENT d-WAVE ORDER

High temperature superconductivity in the cuprates remains one of the most widely investigated, constantly surprising and poorly understood phenomena in physics. Here, we describe briefly a new phenomenological theory inspired by the celebrated description of superconductivity due to Ginzburg and Landau and believed to describe its essence. This posits a free energy functional for the superconductor in terms of a complex order parameter characterizing it. We propose that there is, for superconducting cuprates, a similar functional of the complex, in plane, nearest neighbor spin singlet bond (or Cooper) pair amplitude psi(ij). Further, we suggest that a crucial part of it is a (short range) positive interaction between nearest neighbor bond pairs, of strength J'. Such an interaction leads to nonzero long wavelength phase stiffness or superconductive long range order, with the observed d-wave symmetry, below a temperature T-c similar to zJ' where z is the number of nearest neighbors; d-wave superconductivity is thus an emergent, collective consequence. Using the functional, we calculate a large range of properties, e. g., the pseudogap transition temperature T* as a function of hole doping x, the transition curve T-c(x), the superfluid stiffness rho(s)(x, T), the specific heat (without and with a magnetic field) due to the fluctuating pair degrees of freedom and the zero temperature vortex structure. We find remarkable agreement with experiment. We also calculate the self-energy of electrons hopping on the square cuprate lattice and coupled to electrons of nearly opposite momenta via inevitable long wavelength Cooper pair fluctuations formed of these electrons. The ensuing results for electron spectral density are successfully compared with recent experimental results for angle resolved photo emission spectroscopy (ARPES), and comprehensively explain strange features such as temperature dependent Fermi arcs above T-c and the ``bending'' of the superconducting gap below T-c.

PAMAM Dendrimer-Drug Interactions: Effect of pH on the Binding and Release Pattern

Understanding the dendrimer-drug interaction is of great importance to design and optimize the dendrimer-based drug delivery system. Using atomistic molecular dynamics (MD) simulations, we have analyzed the release pattern of four ligands (two soluble drugs, namely, salicylic acid (Sal), L-alanine (Ala), and two insoluble drugs, namely, phenylbutazone (Pbz) and primidone (Prim)), which were initially encapsulated inside the ethylenediamine (EDA) cored polyamidoamine (PAMAM) dendrimer using the docking method. We have computed the potential of mean force (PMF) variation with generation 5 (G5)-PAMAM dendrimer complexed with drug molecules using umbrella sampling. From our calculated PMF values, we observe that soluble drugs (Sal and Ala) have lower energy barriers than insoluble drugs (Pbz and Prim). The order of ease of release pattern for these drugs from G5 protonated PAMAM dendrimer was found to be Ala > Sal > Prim > Pbz. In the case of insoluble drugs (Prim and Pbz), because of larger size, we observe much nonpolar contribution, and thus, their larger energy barriers can be reasoned to van der Waals contribution. From the hydrogen bonding analysis of the four PAMAM drug complexes under study, we found intermolecular hydrogen bonding to show less significant contribution to the free energy barrier. Another interesting feature appears while calculating the PMF profile of G5NP (nonprotonated)-PAMAM Pbz and G5NP (nonprotonated)-PAMAM-Sal complex. The PMF was found to be less when the drug is bound to nonprotonated dendrimer compared to the protonated dendrimer. Our results suggest that encapsulation of the drug molecule into the host PAMAM dendrimer should be carried out at higher pH values (near pH 10). When such complex enters the human body, the pH is around 7.4 and at that physiological pH, the dendrimer holds the drug tightly. Hence the release of drug can occur at a controlled rate into the bloodstream. Thus, our findings provide a microscopic picture of the encapsulation and controlled release of drugs in the case of dendrimer-based host-guest systems.

Dependence of Protein Adsorption on Wetting Behavior of UHMWPE-HA-Al2O3-CNT Hybrid Biocomposites

Ultrahigh-molecular-weight polyethylene (UHMWPE) is used as an articulating surface in total hip and knee joint replacement. In order to enhance long-term durability/wear resistance properties, UHMWPE-based polymer-ceramic hybrid composites are being developed. Surface properties such as wettability and protein adsorption alter with reinforcement or with change in surface chemistry. From this perspective, the wettability and protein adsorption behavior of compression-molded UHMWPE-hydroxyapatite (HA)-aluminum oxide (Al2O3)-carbon nanotube (CNT) composites were analyzed in conjunction with surface roughness. The combined effect of Al2O3 and CNT shows enhancement of the contact angle by similar to 37A degrees compared with the surface of the UHMWPE matrix reinforced with HA. In reference to unreinforced UHMWPE, protein adsorption density also increased by similar to 230% for 2 wt.%HA-5 wt.%Al2O3-2 wt.%CNT addition to UHMWPE. An important conclusion is that the polar and dispersion components of the surface free energy play a significant role in wetting and protein adsorption than do the total free energy or chemistry of the surface. The results of this study have major implications for the biocompatibility of these newly developed biocomposites.

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Supersymmetric Chern-Simons theories with vector matter

In this paper we discuss SU(N) Chern-Simons theories at level k with both fermionic and bosonic vector matter. In particular we present an exact calculation of the free energy of the N = 2 supersymmetric model (with one chiral field) for all values of the `t Hooft coupling in the large N limit. This is done by using a generalization of the standard Hubbard-Stratanovich method because the SUSY model contains higher order polynomial interactions.

Thermodynamics of TmRhO3, Phase Equilibria, and Chemical Potentials in the System Tm-Rh-O

Phase equilibria in the system Tm-Rh-O at 1200 K is established by isothermal equilibration of selected compositions and phase identification after quenching to room temperature. Six intermetallic phases (Tm3Rh, Tm7Rh3, Tm5Rh3, Tm3Rh2, TmRh, TmRh2 +/-delta) and a ternary oxide TmRhO3 are identified. Based on experimentally determined phase relations, a solid-state electrochemical cell is devised to measure the standard free energy of formation of orthorhombic perovskite TmRhO3 from cubic Tm2O3 and beta-Rh2O3 in the temperature range from (900 to 1300) K. The results can be summarized as: Delta G(f,ox)(o) +/- 104/J.mol(-1) = -46474 + 3.925(T/K). Invoking the Neumann-Kopp rule, the standard enthalpy of formation of TmRhO3 from its constituent elements at 298.15 K is estimated as -1193.89 (+/- 2.86) kJ.mol(-1). The standard entropy of TmRhO3 at 298.15 K is evaluated as 103.8 (+/- 1.6) J.mol(-1).K-1. The oxygen potential-composition diagram and three-dimensional chemical potential diagram at 1200 K and temperature-composition diagrams at constant partial pressures of oxygen are computed from thermodynamic data. The compound TmRhO3 decomposes at 1688 (+/- 2) K in pure oxygen and at 1583 (+/- 2) K in air at standard pressure.

Characterization and iodide adsorption behaviour of HDPY+ modified bentonite

Owing to its favourable physical, chemical and rheological properties, densely compacted bentonite or bentonite-sand mix is considered as a suitable buffer material in deep geological repositories to store high level nuclear waste. Iodine-129 is one of the significant nuclides in the high level waste owing to its long half life and poor sorption onto most geologic media. Bentonite by virtue of negatively charged surface has negligible affinity to retain iodide ions. As organo-bentonites are known to retain iodide ions, the present study characterizes hexadecylpyridinium chloride (HDPyCl.H2O) treated bentonite from Barmer India (referred as HDPy+B) for physico-chemical properties, engineering properties and the iodide adsorption behavior of the organo clay. Batch experiments revealed that HDPy+ ions are largely retained (94 % retention) via cation exchange; the ion-exchange process neutralizes the negative surface charge and bridges clay particles leading to reduction in Atterberg limits, clay content and sediment volume. The organo clay retains iodide by Coulombic attraction (at primary sites) and anion exchange (at secondary sites). The free-energy change (Delta G (o) = -25.5 kJ/mol) value indicated that iodide retention by organo clay is favored physical adsorption process. Iodide adsorption capacity of organo clay decreased significantly (85-100 %) on dilution with 50-80 % bentonite. On the other hand, dilution of bentonite with 50 % organo clay caused 58 % reduction in swell potential and 21 % reduction in swell pressure.

Functional characterization, homology modeling and docking studies of beta-glucosidase responsible for bioactivation of cyanogenic hydroxynitrile glucosides from Leucaena leucocephala (subabul)

Glycosyl hydrolase family 1 beta-glucosidases are important enzymes that serve many diverse functions in plants including defense, whereby hydrolyzing the defensive compounds such as hydroxynitrile glucosides. A hydroxynitrile glucoside cleaving beta-glucosidase gene (Llbglu1) was isolated from Leucaena leucocephala, cloned into pET-28a (+) and expressed in E. coli BL21 (DE3) cells. The recombinant enzyme was purified by Ni-NTA affinity chromatography. The optimal temperature and pH for this beta-glucosidase were found to be 45 A degrees C and 4.8, respectively. The purified Llbglu1 enzyme hydrolyzed the synthetic glycosides, pNPGlucoside (pNPGlc) and pNPGalactoside (pNPGal). Also, the enzyme hydrolyzed amygdalin, a hydroxynitrile glycoside and a few of the tested flavonoid and isoflavonoid glucosides. The kinetic parameters K (m) and V (max) were found to be 38.59 mu M and 0.8237 mu M/mg/min for pNPGlc, whereas for pNPGal the values were observed as 1845 mu M and 0.1037 mu M/mg/min. In the present study, a three dimensional (3D) model of the Llbglu1 was built by MODELLER software to find out the substrate binding sites and the quality of the model was examined using the program PROCHEK. Docking studies indicated that conserved active site residues are Glu 199, Glu 413, His 153, Asn 198, Val 270, Asn 340, and Trp 462. Docking of rhodiocyanoside A with the modeled Llbglu1 resulted in a binding with free energy change (Delta G) of -5.52 kcal/mol on which basis rhodiocyanoside A could be considered as a potential substrate.

Translocation and encapsulation of siRNA inside carbon nanotubes

We report spontaneous translocation of small interfering RNA (siRNA) inside carbon nanotubes (CNTs) of various diameters and chirality using all atom molecular dynamics simulations with explicit solvent. We use umbrella sampling method to calculate the free energy landscape of the siRNA entry and translocation event. Free energy profiles show that siRNA gains free energy while translocating inside CNT, and barrier for siRNA exit from CNT ranges from 40 to 110 kcal/mol depending on CNT chirality and salt concentration. The translocation time tau decreases with the increase of CNT diameter with a critical diameter of 24 angstrom for the translocation. In contrast, double strand DNA of the same sequence does not translocate inside CNT due to large free energy barrier for the translocation. This study helps in understanding the nucleic acid transport through nanopores at microscopic level and may help designing carbon nanotube based sensor for siRNA. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4773302]

Electrochemical determination of thermodynamic properties of DyRhO3 and phase relations in the system Dy-Rh-O

Thermodynamic properties of Dysprosium rhodite (DyRhO3) are measured in the temperature range from 900 to 1,300 K using a solid-state electrochemical cell incorporating yttria-stabilized zirconia as the electrolyte. The standard Gibbs free energy of formation of DyRhO3 with O-type perovskite structure from its components binary oxides, Dysprosia with C-rare earth structure and beta-Rh2O3 with orthorhombic structure, can be represented by the equation: Delta G(f(OX))(O) (+/- 182)/J mol(-1) = -52710+3.821(T/K). By using the thermodynamic data for DyRhO3 from experiment and auxiliary data for other phases from the literature, the phase relations in the system Dy-Rh-O are computed. Thermodynamic data for intermetallic phases in the binary system Dy-Rh, required for constructing the chemical potential diagrams, are evaluated using calorimetric data available in the literature for three intermetallics and Miedema's model, consistent with the phase diagram. The results are presented in the form of Gibbs triangle, oxygen potential-composition diagram, and three-dimensional chemical potential diagram at 1,273 K. Temperature-composition diagrams at constant oxygen partial pressures are also developed. The decomposition temperature of DyRhO3 is 1,732 (+/- 2.5) K in pure oxygen and 1,624 (+/- 2.5) K and in air at standard pressure.

System Gd-Rh-O: Thermodynamics and phase relations

Thermodynamic properties of GdRhO3 are investigated in the temperature range from 900 to 1300 K by employing a solid-state electrochemical cell, incorporating calcia-stabilized zirconia as the electrolyte. The standard Gibbs free energy of formation of GdRhO3 from component binary oxide Gd2O3 with C-rare earth structure and Rh2O3 with orthorhombic structure can be expressed as; Delta G(f(ox))(o)(+/- 60)/J mol(-1) = -56603 + 3.78(T/K) Based on the thermodynamic information on GdRhO3 from experiment and auxiliary data for binary oxides from the literature and estimated properties of Gd-Rh alloys, phase relations are computed for the system Gd-Rh-O at 1273 K. Gibbs free energies for intermetallic phases in the binary Gd-Rh are evaluated using calorimetric data available in the literature for two compositions and Miedema's model, consistent with the binary phase diagram. Isothermal section of the ternary phase diagram, oxygen potential-composition diagram and a 3-D chemical potential diagram for the system Gd-Rh-O at 1273 K are developed. Phase relations in the ternary Gd-Rh-O are also computed as a function of temperature at constant oxygen partial pressures. The ternary oxide, GdRhO3 decomposes to Gd2O3 with B-rare earth structure, metallic Rh and O-2 at 1759(+/- 2) K in pure O-2 and 1649(+/- 2) K in air at a total pressure P-0 -0.1 MPa. (c) 2012 Elsevier B.V. All rights reserved.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

Molecular dynamics simulations of PPI dendrimer-drug complexes

Dendrimeric nanoparticles are potential drug delivery devices which can enhance the solubility of hydrophobic drugs, thus increasing their bioavailability and sustained release action. A quantitative understanding of the dendrimer-drug interactions can give valuable insight into the solubility and release profile of hydrophobic drug molecules in various solvent conditions. Fully atomistic molecular dynamics (MD) simulations have been performed to study the interactions of G5 PPIEDA (G5 ethylenediamine cored poly(propylene imine)) dendrimer and two well known drugs (Famotidine and Indomethacin) at different pH conditions. The study suggested that at low pH the dendrimer-drug complexes are thermodynamically unstable as compared to neutral and high pH conditions. Calculated Potential of Mean Force (PMF) by umbrella sampling showed that the release of drugs from the dendrimer at low pH is spontaneous, median release at neutral pH and slow release at high pH. In addition, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding free energy calculations were also performed at each umbrella sampling window to identify the various energy contributions. To understand the effect of dendrimer chemistry and topology on the solubility and release profile of drugs, this study is extended to explore the solubility and release profile of phenylbutazone drug complexed with G3 poly(amidoamine) and G4 diaminobutane cored PPI dendrimers. The results indicate that the pH-induced conformational changes in dendrimer, ionization states, dendrimer type and pK(a) of the guest molecules influence the free energy barrier and stability of complexation, and thus regulate drug loading, solubility and release.