Molecular probe dynamics reveals suppression of ice-like regions in strongly confined supercooled water

The structure of the hydrogen bond network is a key element for understanding water's thermodynamic and kinetic anomalies. While ambient water is strongly believed to be a uniform, continuous hydrogen-bonded liquid, there is growing consensus that supercooled water is better described in terms of distinct domains with either a low-density ice-like structure or a high-density disordered one. We evidenced two distinct rotational mobilities of probe molecules in interstitial supercooled water of polycrystalline ice Banerjee D, et al. (2009) ESR evidence for 2 coexisting liquid phases in deeply supercooled bulk water. Proc Natl Acad Sci USA 106: 11448-11453]. Here we show that, by increasing the confinement of interstitial water, the mobility of probe molecules, surprisingly, increases. We argue that loose confinement allows the presence of ice-like regions in supercooled water, whereas a tighter confinement yields the suppression of this ordered fraction and leads to higher fluidity. Compelling evidence of the presence of ice-like regions is provided by the probe orientational entropy barrier which is set, through hydrogen bonding, by the configuration of the surrounding water molecules and yields a direct measure of the configurational entropy of the same. We find that, under loose confinement of supercooled water, the entropy barrier surmounted by the slower probe fraction exceeds that of equilibrium water by the melting entropy of ice, whereas no increase of the barrier is observed under stronger confinement. The lower limit of metastability of supercooled water is discussed.

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.

Dual Role of Acetate as a Nucleophile and as an Internal Base in Cycloplatination Reaction of sym-N,N `,N `'-Triarylguanidines

Reaction of cis-Cl2Pt(S(O)Me-2)(2)] with 1 equiv of sym-N,N',N `'-triarylguanidines, ArN=C(NHAr)(2) (sym = symmetrical; Ar = 2-MeC6H4 (LH22-tolyl), 2-(MeO)C6H4 (LH22-anisyl), 4-MeC6H4 (LH24-tolyl), 2,5-Me2C6H3 (LH22,5-xylyl), and 2,6-Me2C6H3 (LH22,6-xylyl)) in toluene under reflux condition for 3 h afforded cis- or trans-Cl2Pt(S(O)Me-2)(ArN=C(NHAr)(2))] (Ar = 2-MeC6H4 (1), 2-(MeO)C6H4 (2), 4-MeC6H4 (3), 2,5-h Me2C6H3 (4), and 2,6-Me2C6H3 (5), respectively) in 83-96% yield. Reaction of cis-Cl2Pt(S(O)Me-2)(2)] with 1 equiv of LH22-tolyl and LH24-tolyl in the presence of 1 equiv of NaOAc in methanol under reflux condition for 3 h afforded acetate-substituted products, cis-(AcO)ClPt(S(O)Me-2)(ArN=C(NHAr)(2))] (Ar = 2-MeC6H4 (6) and 4-MeC6H4 (7)) in 83% and 84% yields, respectively. Reaction of cis-Cl2Pt(S(O)Me-2)(2)] with 1 equiv of LH22-anisyl and LH22-tolyl in the presence of 1 equiv of NaOAc in methanol under reflux condition for 3 and 12 h afforded six-membered C,N] platinacycles, Pt{kappa(2)(C,N)-C6H3R-3(NHC(NHAr)(=NAr))-2}Cl(S(O)Me-2)] (Ar = 2-RC6H4; R = OMe (8) and Me (9)), in 92% and 79% yields, respectively. The new complexes have been characterized by analytical and spectroscopic techniques, and further the molecular structures of 1, 2, 4, 5, 6, and 8 have been determined by single-crystal X-ray diffraction. The platinum atom in 1, 4, and 5 exhibited the trans configuration, while that in 2, 6, and 8 exhibited the cis configuration. Complex 6 is shown to be the precursor for 9, and the former is suggested to transform to the latter possibly via an intramolecular C-H activation followed by elimination of AcOH. The solution behavior of new complexes has been studied by multinuclear NMR (H-1, Pt-195, and C-13) spectroscopy. The new complexes exist exclusively as a single isomer (trans (1 and 5) and cis (6 and 7)), a mixture of cis and trans isomers with the former isomer being predominant in the case of 2 and the latter isomer being predominant in the case of 3. Complex 5 in the trans form revealed the presence of one isomer at 0.007 mM concentration and two isomers in about 1.00:0.12 ratio at 0.154 mM concentration as revealed by H-1 NMR spectroscopy, and this has been ascribed to the restricted Pt-S bond rotation at higher concentration. Platinacycle 8 exists as one isomer, while 9 exists as a mixture of seven isomers in solution. The influence of steric factor, pi-acceptor property of the guanidine, subtle solid-state packing forces upon the configuration of the platinum atom, and the number of isomers in solution have been outlined. Factors that accelerate or slow down the cycloplatination reaction, the role of NaOAc, and a plausible mechanism of this reaction have been discussed.

Classification of Vehicles using Magnetic Dipole Model

This paper presents an efficient approach to the modeling and classification of vehicles using the magnetic signature of the vehicle. A database was created using the magnetic signature collected over a wide range of vehicles(cars). A vehicle is modeled as an array of magnetic dipoles. The strength of the magnetic dipole and the separation between the magnetic dipoles varies for different vehicles and is dependent on the metallic composition and configuration of the vehicle. Based on the magnetic dipole data model, we present a novel method to extract a feature vector from the magnetic signature. In the classification of vehicles, a linear support vector machine configuration is used to classify the vehicles based on the obtained feature vectors.

An estimate of equilibrium sensitivity of global terrestrial carbon cycle using NCAR CCSM4

Increasing concentrations of atmospheric CO2 influence climate, terrestrial biosphere productivity and ecosystem carbon storage through its radiative, physiological and fertilization effects. In this paper, we quantify these effects for a doubling of CO2 using a low resolution configuration of the coupled model NCAR CCSM4. In contrast to previous coupled climate-carbon modeling studies, we focus on the near-equilibrium response of the terrestrial carbon cycle. For a doubling of CO2, the radiative effect on the physical climate system causes global mean surface air temperature to increase by 2.14 K, whereas the physiological and fertilization on the land biosphere effects cause a warming of 0.22 K, suggesting that these later effects increase global warming by about 10 % as found in many recent studies. The CO2-fertilization leads to total ecosystem carbon gain of 371 Gt-C (28 %) while the radiative effect causes a loss of 131 Gt-C (10 %) indicating that climate warming damps the fertilization-induced carbon uptake over land. Our model-based estimate for the maximum potential terrestrial carbon uptake resulting from a doubling of atmospheric CO2 concentration (285-570 ppm) is only 242 Gt-C. This highlights the limited storage capacity of the terrestrial carbon reservoir. We also find that the terrestrial carbon storage sensitivity to changes in CO2 and temperature have been estimated to be lower in previous transient simulations because of lags in the climate-carbon system. Our model simulations indicate that the time scale of terrestrial carbon cycle response is greater than 500 years for CO2-fertilization and about 200 years for temperature perturbations. We also find that dynamic changes in vegetation amplify the terrestrial carbon storage sensitivity relative to a static vegetation case: because of changes in tree cover, changes in total ecosystem carbon for CO2-direct and climate effects are amplified by 88 and 72 %, respectively, in simulations with dynamic vegetation when compared to static vegetation simulations.

An approach for assessing the accuracy of accelerometer responses during axial impact tests

The primary objective of the present study is to show that for the most common configuration of an impactor system, the accelerometer cannot exactly reproduce the dynamic response of a specimen subjected to impact loading. An equivalent Lumped Parameter Model (LPM) of the given impactor set-up has been formulated for assessing the accuracy of an accelerometer mounted in a drop-weight impactor set-up for an axially loaded specimen. A specimen under the impact loading is represented by a non-linear spring of varying stiffness, while the accelerometer is assumed to behave in a linear manner due to its high stiffness. Specimens made of steel, aluminium and fibre-reinforced composite (FRC) are used in the present study. Assuming the force-displacement response obtained in an actual impact test to be the true behaviour of the test specimen, a suitable numerical approach has been used to solve the governing non-linear differential equations of a three degrees-of-freedom (DOF) system in a piece-wise linear manner. The numerical solution of the governing differential equations following an explicit time integration scheme yields an excellent reproduction of the mechanical behaviour of the specimen, consequently confirming the accuracy of the numerical approach. However, the spring representing the accelerometer predicts a response that qualitatively matches the assumed force-displacement response of the test specimen with a perceptibly lower magnitude of load.

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.

On the nature of meandering of the springtime western boundary current in the Bay of Bengal

We present evidence that the springtime western boundary current (WBC) in the Bay of Bengal is a continuous northward-flowing current from about 12 degrees N to 17 degrees N, which then separates from the coast at around 18 degrees N. We first revisit a hydrographic data set collected in 1987 from a potential vorticity perspective, and then analyze absolute dynamic height maps from satellite altimeters during the period 2000-2010. The altimetric maps suggest that the mean configuration of the WBC is that of an intense current with two anticyclonic eddies on the offshore side, which are part of the basin-wide anticyclonic circulation. The WBC consistently separates from the coast at around 18 degrees N in all years between 2000 and 2010. The path of the eastward-flowing mean stream after separation appears to be consistent with isolines of f/H and with Ertel's potential vorticity, based on an analysis of the hydrographic data from 1987.

Evaluation of the accuracy of an impact load cell using lumped parameter modeling and analysis

The objective of the current study is to evaluate the fidelity of load cell reading during impact testing in a drop-weight impactor using lumped parameter modeling. For the most common configuration of a moving impactor-load cell system in which dynamic load is transferred from the impactor head to the load cell, a quantitative assessment is made of the possible discrepancy that can result in load cell response. A 3-DOF (degrees-of-freedom) LPM (lumped parameter model) is considered to represent a given impact testing set-up. In this model, a test specimen in the form of a steel hat section similar to front rails of cars is represented by a nonlinear spring while the load cell is assumed to behave in a linear manner due to its high stiffness. Assuming a given load-displacement response obtained in an actual test as the true behavior of the specimen, the numerical solution of the governing differential equations following an implicit time integration scheme is shown to yield an excellent reproduction of the mechanical behavior of the specimen thereby confirming the accuracy of the numerical approach. The spring representing the load cell, however,predicts a response that qualitatively matches the assumed load-displacement response of the test specimen with a perceptibly lower magnitude of load.

Evaluation of the accuracy of an accelerometer response generated by axial impact loading

The basic objective in the present study is to show that for the most common configuration of an impactor system, an accelerometer cannot exactly reproduce the dynamic response of a specimen subject to impact loading. Assessment of the accelerometer mounted in a drop-weight impactor setup for an axially loaded specimen is done with the aid of an equivalent lumped parameter model (LPM) of the setup. A steel hat-type specimen under the impact loading is represented as a non-linear spring of varying stiffness, while the accelerometer is assumed to behave in a linear manner due to its high stiffness. A suitable numerical approach has been used to solve the non-linear governing equations for a 3 degrees-of-freedom system in a piece-wise linear manner. The numerical solution following an explicit time integration scheme is used to yield an excellent reproduction of the mechanical behavior of the specimen thereby confirming the accuracy of the numerical approach. The spring representing the accelerometer, however, predicts a response that qualitatively matches the assumed load–displacement response of the test specimen with a perceptibly lower magnitude of load.

Small-Angle Neutron Scattering Studies of Hemoglobin Confined Inside Silica Tubes of Varying Sizes

In addition to the chemical nature of the surface, the dimensions of the confining host exert a significant influence on confined protein structures; this results in immense biological implications, especially those concerning the enzymatic activities of the protein. This study probes the structure of hemoglobin (Hb), a model protein, confined inside silica tubes with pore diameters that vary by one order of magnitude (approximate to 20-200 nm). The effect of confinement on the protein structure is probed by comparison with the structure of the protein in solution. Small-angle neutron scattering (SANS), which provides information on protein tertiary and quaternary structures, is employed to study the influence of the tube pore diameter on the structure and configuration of the confined protein in detail. Confinement significantly influences the structural stability of Hb and the structure depends on the Si-tube pore diameter. The high radius of gyration (R-g) and polydispersity of Hb in the 20 nm diameter Si-tube indicates that Hb undergoes a significant amount of aggregation. However, for Si-tube diameters greater or equal to 100 nm, the R-g of Hb is found to be in very close proximity to that obtained from the protein data bank (PDB) reported structure (R-g of native Hb=23.8 angstrom). This strongly indicates that the protein has a preference for the more native-like non-aggregated state if confined inside tubes of diameter greater or equal to 100 nm. Further insight into the Hb structure is obtained from the distance distribution function, p(r), and ab initio models calculated from the SANS patterns. These also suggest that the Si-tube size is a key parameter for protein stability and structure.

Curvature based mobility analysis and form closure of smooth planar curves with multiple contacts

This paper presents a simple second-order, curvature based mobility analysis of planar curves in contact. The underlying theory deals with penetration and separation of curves with multiple contacts, based on relative configuration of osculating circles at points of contact for a second-order rotation about each point of the plane. Geometric and analytical treatment of mobility analysis is presented for generic as well as special contact geometries. For objects with a single contact, partitioning of the plane into four types of mobility regions has been shown. Using point based composition operations based on dual-number matrices, analysis has been extended to computationally handle multiple contacts scenario. A novel color coded directed line has been proposed to capture the contact scenario. Multiple contacts mobility is obtained through intersection of the mobility half-spaces. It is derived that mobility region comprises a pair of unbounded or a single bounded convex polygon. The theory has been used for analysis and synthesis of form closure configurations, revolute and prismatic kinematic pairs. (C) 2013 Elsevier Ltd. All rights reserved.

Effect of Substrate Temperature on the Microstructure Variation in Sputter-deposited Hydrogenated Amorphous Silicon Thinfilms

Vacancy, void incorporation and Si-H-x configuration in hydrogenated amorphous silicon (a-Si:H) thin films was studied. Films were grown by Direct Current (DC), pulsed DC and Radio Frequency (RF) magnetron sputtering. Fourier Transform Infrared (FTIR) spectroscopic analysis has been carried out on the films and found that, the a-Si: H films grown by DC magnetron sputtering are of good quality compared to pulsed DC and RF deposited films. The effect of Substrate temperature (T-S) on the total hydrogen concentration (C-H), configuration of hydrogen bonding, density (decided by the vacancy and void incorporation) and the microstructure factor (R*) was studied. T-S is found to be an active parameter in affecting the above said properties of the films. The films contain both vacancies and voids. At low hydrogen dilutions the films are vacancy dominated and at high hydrogen dilutions they are void dominated. It is found that T-S favors monohydride (Si-H) bonding at the cost of dihydride (Si-H-2) bonding. This dividing line is at C-H=14 at.% H for DC sputter deposited films. The microstructure structure factor R* is found to be zero for as deposited DC films at T-S=773K. The threshold C-H for void dominated region is found to be C-H=23 at.% H for RF, C-H=18 at.% H for PDC and C-H similar to 14 at.%H for DC sputter deposited films.

Energy funneling and macromolecular conformational dynamics: a 2D Raman correlation study of PEG melting

Cascading energy landscapes through funneling has been postulated as a mechanistic route for achieving the lowest energy configuration of a macromolecular system (such as proteins and polymers). In particular, understanding the molecular mechanism for the melting and crystallization of polymers is a challenging fundamental question. The structural modifications that lead to the melting of poly(ethylene glycol) (PEG) are investigated here. Specific Raman bands corresponding to different configurations of the PEG chain have been identified, and the molecular structural dynamics of PEG melting have been addressed using a combination of Raman spectroscopy, 2D Raman correlation and density functional theory (DFT) calculations. The melting dynamics of PEG have been unambiguously explained along the C-O bond rotation coordinate.

Multidimensional free energy surface of unfolding of HP-36: Microscopic origin of ruggedness

The protein folding funnel paradigm suggests that folding and unfolding proceed as directed diffusion in a multidimensional free energy surface where a multitude of pathways can be traversed during the protein's sojourn from initial to final state. However, finding even a single pathway, with the detail chronicling of intermediates, is an arduous task. In this work we explore the free energy surface of unfolding pathway through umbrella sampling, for a small globular a-helical protein chicken-villin headpiece (HP-36) when the melting of secondary structures is induced by adding DMSO in aqueous solution. We find that the unfolding proceeds through the initial separation or melting of aggregated hydrophobic core that comprises of three phenylalanine residues (Phe7, Phe11, and Phe18). This separation is accompanied by simultaneous melting of the second helix. Unfolding is found to be a multistage process involving crossing of three consecutive minima and two barriers at the initial stage. At a molecular level, Phe18 is observed to reorient itself towards other hydrophobic grooves to stabilize the intermediate states. We identify the configuration of the intermediates and correlate the intermediates with those obtained in our previous works. We also give an estimate of the barriers for different transition states and observe the softening of the barriers with increasing DMSO concentration. We show that higher concentration of DMSO tunes the unfolding pathway by destabilizing the third minimum and stabilizing the second one, indicating the development of a solvent modified, less rugged pathway. The prime outcome of this work is the demonstration that mixed solvents can profoundly transform the nature of the energy landscape and induce unfolding via a modified route. A successful application of Kramer's rate equation correlating the free energy simulation results shows faster rate of unfolding with increasing DMSO concentration. This work perhaps presents the first systematic theoretical study of the effect of a chemical denaturant on the microscopic free energy surface and rates of unfolding of HP-36. (C) 2014 AIP Publishing LLC.