Charge-transfer-driven inclusion of neutral TCNQ molecules in the galleries of a layered double hydroxide: an experimental and computational study

The layered double hydroxides (LDH) or anionic clays are an important class of ion-exchange materials. They consist of positively charged brucite-like inorganic sheets with charge-compensating exchangeable anions in the interlamellar space. Here we show how neutral TCNQ (7,7,8,8-tetracyanoquinodimethane) molecules can be included within the galleries of an LDH. To do so, we exploit the fact that TCNQ is a good electron acceptor that forms donor acceptor complexes with a variety of donors. The electron donor aniline was intercalated into a Mg-Al LDH as p-aminobenzoate (AB) ions by a conventional ion-exchange reaction. We show here that neutral TCNQ molecules may be driven into the galleries of the layered solid by charge-transfer complex formation with the intercalated p-aminobenzoate anions. We use diffraction and spectroscopic measurements in combination with molecular dynamics simulations and quantum chemical calculations to establish the nature of interactions and arrangement of the charge-transfer complex within the galleries of the layered double hydroxide. Electrostatic interactions between the TCNQ molecules and the anchored AB ions, subsequent to charge transfer, are the driving force for the inclusion of TCNQ molecules in the galleries of the LDH.

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Switching and release dynamics of an electrostatically actuated MEMS switch under the influence of squeeze-film damping

This paper reports analytical modeling, simulation and experimental validation for switching and release times of an electrostatically actuated micromachined switch. Presented work is an extension of our earlier work [1] that analytically argued, and numerically and experimentally demonstrated, why pull-in time is larger that pull-up time when the actuation voltage is less than twice of the pull-in voltage. In this paper, switching dynamics is investigated under the influence of squeeze-film damping. Tests were performed on SOI (silicon-on-insulator) based parallel beams structures.Typical voltage requirement for actuation is in the range of 10-30 V. All the experiments were performed in normal atmospheric pressure. Measurement results confirm that the quality factor Q has appreciable effect on the release time compared to the switching time. The quality factor Q is extracted from the response measurement and compared with the ANSYS simulation result. In addition, the dynamic pull-in effect has also been studied and reported in this paper. A contribution of this work includes the effect of various phenomena such as squeeze-film damping, dynamic pull-in, and frequency pull-in effects on the switching dynamics of a MEMS switch.

Lysis dynamics and membrane oligomerization pathways for Cytolysin A (ClyA) pore-forming toxin

Pore-forming toxins are known for their ability to efficiently form transmembrane pores which eventually leads to cell lysis. The dynamics of lysis and underlying self-assembly or oligomerization pathways leading to pore formation are incompletely understood. In this manuscript the pore-forming kinetics and lysis dynamics of Cytolysin-A (ClyA) toxins on red blood cells (RBCs) are quantified and compared with experimental lysis data. Lysis experiments are carried out on a fixed mass of RBCs, under isotonic conditions in phosphate-buffered saline, for different initial toxin concentrations ranging from 2.94-14.7 nM. Kinetic models which account for monomer binding, conformation and oligomerization to form the dodecameric ClyA pore complex are developed and lysis is assumed to occur when the number of pores per RBC (n(p)) exceeds a critical number, n(pc). By analysing the model in a sublytic regime (n(p) < n(pc)) the number of pores per RBC to initiate lysis is found to lie between 392 and 768 for the sequential oligomerization mechanism and between 5300 and 6300 for the non-sequential mechanism. Rupture rates which are first order in the number of RBCs are seen to provide the best agreement with the lysis experiments. The time constants for pore formation are estimated to lie between 1 and 20 s and monomer conformation time scales were found to be 2-4 times greater than the oligomerization times. Cell rupture takes places in 100s of seconds, and occurs predominantly with a steady number of pores ranging from 515 to 11 000 on the RBC surface for the sequential mechanism. Both the sequential irreversible and non-sequential kinetics provide similar predictions of the hemoglobin release dynamics, however the hemoglobin released as a function of the toxin concentration was accurately captured only with the sequential model. Each mechanism develops a distinct distribution of mers on the surface, providing a unique experimentally observable fingerprint to identify the underlying oligomerization pathways. Our study offers a method to quantify the extent and dynamics of lysis which is an important aspect of developing novel drug and gene delivery strategies based on pore-forming toxins.

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.

Experimental investigation of biofuel drop impact on stainless steel surface

Blends of conventional fuels such as Jet-A1 (aviation kerosene) and diesel with bio-derived components, referred to as biofttels, are gradually replacing the conventional fuels in aircraft and automobile engines. There is a lack of understanding on the interaction of spray drops of such biofuels with solid surfaces. The present study is an experimental investigation on the impact of biofuel drops onto a smooth stainless steel surface. The biofuel is a mixture of 90% commercially available camelina-derived biofuel and 10% aromatics. Biofuel drops were generated using a syringe-hypodermic needle arrangement. On demand, the needle delivers an almost spherical drop with drop diameter in the range 2.05-2.15 mm. Static wetting experiments show that the biofuel drop completely wets the stainless steel surface and exhibits an equilibrium contact angle of 5.6. High speed video camera was used to capture the impact dynamics of biofuel drops with Weber number ranging from 20 to 570. The spreading dynamics and maximum spreading diameter of impacting biofuel drops on the target surface were analyzed. For the impact of high Weber number biofuel drops, the spreading law suggests beta similar to tau(0.5) where beta is the spread factor and tau, the nondimensionalized time. The experimentally observed trend of maximum spread factor, beta(max) of camelina biofuel drop on the target surface with We compares well with the theoretically predicted trend from Ukiwe-Kwok model. After reaching beta(max), the impacting biofuel drop undergoes a prolonged sluggish spreading due to the high wetting nature of the camelina biofuel-stainless steel system. As a result, the final spread factor is found to be a little more than beta(max). (C) 2014 Elsevier Inc. All rights reserved.

Experimental signatures of a nonequilibrium phase transition governing the yielding of a soft glass

We present direct experimental signatures of a nonequilibrium phase transition associated with the yield point of a prototypical soft solid-a binary colloidal glass. By simultaneously quantifying single-particle dynamics and bulk mechanical response, we identified the threshold for the onset of irreversibility with the yield strain. We extracted the relaxation time from the transient behavior of the loss modulus and found that it diverges in the vicinity of the yield strain. This critical slowing down is accompanied by a growing correlation length associated with the size of regions of high Debye-Waller factor, which are precursors to yield events in glasses. Our results affirm that the paradigm of nonequilibrium critical phenomena is instrumental in achieving a holistic understanding of yielding in soft solids.

A molecular dynamics study of ambient and high pressure phases of silica: Structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, alpha-PbO2-type, and pyrite-type for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase. (C) 2014 AIP Publishing LLC.

Oscillation dynamics of sessile droplets subjected to substrate vibration

Sessile droplets on a vibrating substrate are investigated focusing on axisymmetric oscillations with pinned contact line. Proper orthogonal decomposition is employed to identify the different modes of droplet shape oscillation and quantitatively assess the droplet oscillation and spectral response. We offer the first experimental evidence for the analogy of an oscillating sessile droplet with a non-linear spring mass damper system. The qualitative and quantitative agreement of amplitude response and phase response curves and limit cycles of the model dynamical system with that observed experimentally suggest that the bulk oscillations in the fundamental mode of a sessile droplet can be very well modeled by a Duffing oscillator with a hard spring, especially near the resonance. The red shift of the resonance peak with an increase in the glycerol concentration is clearly evidenced by both the experimental and predicted amplitude response curves. The influence of various operational parameters such as excitation frequency and amplitude and fluid properties on the droplet oscillation characteristics is adequately captured by the model. (C) 2014 Elsevier Ltd. All rights reserved.

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.

De novo inference of protein function from coarse-grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high-throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 mu s coarse-grained (CG) molecular dynamics trajectories were used to compute normalized root-mean-square-fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three-dimensional autocorrelation vectors. Our in-house custom-built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics-signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof-of-principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom-made CG FF, useful to all. (C) 2014 Wiley Periodicals, Inc.

Investigation of selective sensing of a diamine for aldehyde by experimental and simulation studies

An organic molecule-o-phenylene diamine (OPD)-is selected as an aldehyde sensing material. It is studied for selectivity to aldehyde vapours both by experiment and simulation. A chemiresistor based sensor for detection of aldehyde vapours is fabricated. An o-phenylene diamine-carbon black composite is used as the sensing element. The amine groups in the OPD would interact with the carbonyl groups of the aldehydes. The selectivity and cross-sensitivity of the OPD-CB sensor to VOCs aldehyde, ketone and alcohol-are studied. The sensor shows good response to aldehydes compared to other VOCs. The higher response for aldehydes is attributed to the interaction of the carbonyl oxygen of aldehydes with-NH2 groups of OPD. The surface morphology of the sensing element is studied by scanning electron microscopy. The OPD-CB sensor is responsive to 10 ppm of formaldehyde. The interaction of the VOCs with the OPD-CB nanocomposite is investigated by molecular dynamics studies. The interaction energies of the analyte with the OPD-CB nanocomposite were calculated. It is observed that the interaction energies for aldehydes are higher than those for other analytes. Thus the OPD-CB sensor shows selectivity to aldehydes. The simulated radial distribution function is calculated for the O-H pair of analyte and OPD which further supports the finding that the amine groups are involved in the interaction. These results suggest that it is important and easy to identify appropriate sensing materials based on the understanding of analyte interaction properties.

SIMULATION STUDY OF THE ADSORPTION DYNAMICS OF CYLINDRICAL SILICA GEL PARTICLES

The paper presents a simulation study of loose cylindrically shaped particles packed within a copper plate and aluminum fins. The model presented solves coupled heat and mass transfer equations using the finite volume method based on ANSY S FLUENT medium. Three different arrangements of cylindrical particles are considered. The model is validated with experimental data. It is found that the arrangements which represented monolayer configurations are only marginally better in heat transfer and uptake efficiency than the tri-layer configuration in the presence of fins. However, there is an appreciable difference in the uptake curve between monoand tri-layer configurations in the absence of fins. Finally, it is found that the fin pitch also plays an important role in determining the time constant for the adsorber design.

Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.