Structure, dynamics and thermodynamics of single-file water under confinement: effects of polarizability of water molecules

Various structural, dynamic and thermodynamic properties of water molecules confined in single-wall carbon nanotubes (CNTs) are investigated using both polarizable and non-polarizable water models. The inclusion of polarizability quantitatively affects the nature of hydrogen bonding, which governs many properties of confined water molecules. Polarizable water leads to tighter hydrogen bonding and makes the distance between neighboring water molecules shorter than that for non-polarizable water. Stronger hydrogen bonding also decreases the rotational entropy and makes the diffusion constant smaller than in TIP3P and TIP3PM water models. The reorientational dynamics of the water molecules is governed by a jump mechanism, the barrier for the jump being highest for the polarizable water model. Our results highlight the role of polarizability in governing the dynamics of confined water and demonstrate that the inclusion of polarizability is necessary to obtain agreement with the results of ab initio simulations for the distributions of waiting and jump times. The SPC/E water model is found to predict various water properties in close agreement with the results of polarizable water models with much lower computational costs.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

Amorphous W-S-N thin films: The atomic structure behind ultra-low friction

Amorphous W-S-N in the form of thin films has been identified experimentally as an ultra-low friction material, enabling easy sliding by the formation of a WS2 tribofilm. However, the atomic-level structure and bonding arrangements in amorphous W-S-N, which give such optimum conditions for WS2 formation and ultra-low friction, are not known. In this study, amorphous thin films with up to 37 at.% N are deposited, and experimental as well as state-of-the-art ab initio techniques are employed to reveal the complex structure of W-S-N at the atomic level. Excellent agreement between experimental and calculated coordination numbers and bond distances is demonstrated. Furthermore, the simulated structures are found to contain N bonded in molecular form, i.e. N-2, which is experimentally confirmed by near edge X-ray absorption fine structure and X-ray photoelectron spectroscopy analysis. Such N-2 units are located in cages in the material, where they are coordinated mainly by S atoms. Thus this ultra-low friction material is shown to be a complex amorphous network of W, S and N atoms, with easy access to W and S for continuous formation of WS2 in the contact region, and with the possibility of swift removal of excess nitrogen present as N-2 molecules. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Synergistic Effect of Mo plus Cu Codoping on the Photocatalytic Behavior of Metastable TiO2 Solid Solutions

Codoping with Cu and Mo is shown to have a synergistic effect on the photocatalytic activity of TiO2. The enhancement in activity is observed only if the synthesis route results in TiO2 in which (Cu, Mo) codopants are forced into the TiO2 lattice. Using X-ray photoelectron spectroscopy, Cu and Mo are shown to be present in the +2 and +6 oxidation states, respectively. A systematic study of the ternary system shows that TiO2 containing 6 mol % CuO and 1.5 mol % MoO3 is the most active ternary composition. Ab initio calculations show that codoping of TiO2 using (Mo, Cu) introduces levels above the valence band, and below the conduction band, resulting in a significant reduction in the band gap (similar to 0.8 eV). However, codoping also introduces deep defect states, which can have a deleterious impact on photoactivity. This helps rationalize the narrow compositional window over which the enhancement in photocatalytic activity is observed.

Dynamics of the chemical bond: inter- and intra-molecular hydrogen bond

In this discussion, we show that a static definition of a `bond' is not viable by looking at a few examples for both inter-and intra-molecular hydrogen bonding. This follows from our earlier work (Goswami and Arunan, Phys. Chem. Chem. Phys. 2009, 11, 8974) which showed a practical way to differentiate `hydrogen bonding' from `van der Waals interaction'. We report results from ab initio and atoms in molecules theoretical calculations for a series of Rg center dot center dot center dot HX complexes (Rg = He/Ne/Ar and X = F/Cl/Br) and ethane-1,2-diol. Results for the Rg center dot center dot center dot HX/DX complexes show that Rg center dot center dot center dot DX could have a `deuterium bond' even when Rg center dot center dot center dot HX is not `hydrogen bonded', according to the practical criterion given by Goswami and Arunan. Results for ethane-1,2-diol show that an `intra-molecular hydrogen bond' can appear during a normal mode vibration which is dominated by the O center dot center dot center dot O stretching, though a `bond' is not found in the equilibrium structure. This dynamical `bond' formation may nevertheless be important in ensuring the continuity of electron density across a molecule. In the former case, a vibration `breaks' an existing bond and in the later case, a vibration leads to `bond' formation. In both cases, the molecule/complex stays bound irrespective of what happens to this `hydrogen bond'. Both these cases push the borders on the recent IUPAC recommendation on hydrogen bonding (Arunan et al. Pure. Appl. Chem. 2011, 83 1637) and justify the inclusive nature of the definition.

The Self-Association of Graphane Is Driven by London Dispersion and Enhanced Orbital Interactions

We investigated the nature of the cohesive energy between graphane sheets via multiple CH center dot center dot center dot HC interactions, using density functional theory (DFT) including dispersion correction (Grimmes D3 approach) computations of n]graphane sigma dimers (n = 6-73). For comparison, we also evaluated the binding between graphene sheets that display prototypical pi/pi interactions. The results were analyzed using the block-localized wave function (BLW) method, which is a variant of ab initio valence bond (VB) theory. BLW interprets the intermolecular interactions in terms of frozen interaction energy (Delta E-F) composed of electrostatic and Pauli repulsion interactions, polarization (Delta E-pol), charge-transfer interaction (Delta E-CT), and dispersion effects (Delta E-disp). The BLW analysis reveals that the cohesive energy between graphane sheets is dominated by two stabilizing effects, namely intermolecular London dispersion and two-way charge transfer energy due to the sigma CH -> sigma*(HC) interactions. The shift of the electron density around the nonpolar covalent C-H bonds involved in the intermolecular interaction decreases the C-H bond lengths uniformly by 0.001 angstrom. The Delta E-CT term, which accounts for similar to 15% of the total binding energy, results in the accumulation of electron density in the interface area between two layers. This accumulated electron density thus acts as an electronic glue for the graphane layers and constitutes an important driving force in the self-association and stability of graphane under ambient conditions. Similarly, the double faced adhesive tape style of charge transfer interactions was also observed among graphene sheets in which it accounts for similar to 18% of the total binding energy. The binding energy between graphane sheets is additive and can be expressed as a sum of CH center dot center dot center dot HC interactions, or as a function of the number of C-H bonds.

Microwave Spectrum of Hexafluoroisopropanol and Torsional Behavior of Molecules with a CF3-C-CF3 Group

This paper presents the first microwave spectroscopic investigation on hexafluoroisopropanol (HFIP). A pulsed nozzle Fourier transform microwave spectrometer has been used to determine the rotational constants for HFIP as A = 2105.12166(18) MHz, B = 1053.99503(12) MHz, and C = 932.33959(13) MHz. In addition, five isotopologues of HFIP have been observed experimentally to determine the accurate structure of HFIP. The observed spectrum could be assigned to the most stable conformer of HFIP, called antiperiplanar. Available spectroscopic information and ab initio calculations on five prototype molecules helped in exploring the torsional behavior of molecules having a CF3-C-CF3 group. Two-dimensional potential energy surfaces have been analyzed for all molecules, which explained the presence/absence of doubling in the rotational transitions. With the help of natural bond orbital (NBO) analysis, reasons for the conformational preference of HFIP have been explained.

Design of symmetric TIM barrel proteins from first principles

Background: Computational protein design is a rapidly maturing field within structural biology, with the goal of designing proteins with custom structures and functions. Such proteins could find widespread medical and industrial applications. Here, we have adapted algorithms from the Rosetta software suite to design much larger proteins, based on ideal geometric and topological criteria. Furthermore, we have developed techniques to incorporate symmetry into designed structures. For our first design attempt, we targeted the (alpha/beta)(8) TIM barrel scaffold. We gained novel insights into TIM barrel folding mechanisms from studying natural TIM barrel structures, and from analyzing previous TIM barrel design attempts. Methods: Computational protein design and analysis was performed using the Rosetta software suite and custom scripts. Genes encoding all designed proteins were synthesized and cloned on the pET20-b vector. Standard circular dichroism and gel chromatographic experiments were performed to determine protein biophysical characteristics. 1D NMR and 2D HSQC experiments were performed to determine protein structural characteristics. Results: Extensive protein design simulations coupled with ab initio modeling yielded several all-atom models of ideal, 4-fold symmetric TIM barrels. Four such models were experimentally characterized. The best designed structure (Symmetrin-1) contained a polar, histidine-rich pore, forming an extensive hydrogen bonding network. Symmetrin-1 was easily expressed and readily soluble. It showed circular dichroism spectra characteristic of well-folded alpha/beta proteins. Temperature melting experiments revealed cooperative and reversible unfolding, with a T-m of 44 degrees C and a Gibbs free energy of unfolding (Delta G degrees) of 8.0 kJ/mol. Urea denaturing experiments confirmed these observations, revealing a C-m of 1.6 M and a Delta G degrees of 8.3 kJ/mol. Symmetrin-1 adopted a monomeric conformation, with an apparent molecular weight of 32.12 kDa, and displayed well resolved 1D-NMR spectra. However, the HSQC spectrum revealed somewhat molten characteristics. Conclusions: Despite the detection of molten characteristics, the creation of a soluble, cooperatively folding protein represents an advancement over previous attempts at TIM barrel design. Strategies to further improve Symmetrin-1 are elaborated. Our techniques may be used to create other large, internally symmetric proteins.

Anisotropic magnetic couplings and structure-driven canted to collinear transitions in Sr2IrO4 by magnetically constrained noncollinear DFT

We study the canted magnetic state in Sr2IrO4 using fully relativistic density functional theory (DFT) including an on-site Hubbard U correction. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation. Moreover, we put forward a scheme to study the anisotropic magnetic couplings by mapping magnetically constrained noncollinear DFT onto a general spin Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin, and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the structural distortions and the competition between isotropic exchange and DM interactions.

Na2.44Mn1.79(SO4)(3): a new member of the alluaudite family of insertion compounds for sodium ion batteries

Sodium-ion batteries have been extensively pursued as economic alternatives to lithium-ion batteries. Investigating the polyanion chemistry, alluaudite structured Na2Fe2II(SO4)(3) has been recently discovered as a 3.8 V positive electrode material (Barpanda et al., Nature Commun., 5: 4358, 2014). Registering the highest ever Fe-III/Fe-II redox potential (vs. Na/Na+) and formidable energy density, it has opened up a new polyanion family for sodium batteries. Exploring the alluaudite family, here we report isotypical Na2+2xMn2-xII(SO4)(3) (x = 0.22) as a novel high-voltage cathode material for the first time. Following low-temperature (ca. 350 degrees C) solid-state synthesis, the structure of this new alluaudite compound has been solved adopting a monoclinic framework (s.g. C2/c) showing antiferromagnetic ordering at 3.4 K. Synergising experimental and ab initio DFT investigation, Na2+2xMn2-xII(SO4)(3) has been found to be a potential high-voltage (ca. 4.4 V) cathode material for sodium batteries.

X-Ha center dot center dot center dot C hydrogen bonds in n-alkane-HX (X = F, OH) complexes are stronger than C-Ha <-X hydrogen bonds

Computational study of X-Ha <-C and C-Ha <-X hydrogen bonds in n-alkane-HX complexes (X =F,OH, alkane =propane, butane, pentane) has been carried out in this work. Ab initio and density functional theories were used for this study. For n-alkane-H2O complexes both Oa <-H-C and O-Ha <-C hydrogen bonded complex have been found, while for n-alkane-HF complexes, our attempt to optimize Fa <-H-C H-bond was not successful. Like most of the hydrogen bonded systems, strong correlation between binding energy and stretching frequency of H-F and O-H stretching mode was observed. The values of electron density and Laplacian of electron density are within the accepted range for hydrogen bonds. In all these cases, X-Ha <-C hydrogen bonds are found to be stronger than C-Ha <-X hydrogen bonds.

Reconstruction of the Disassembly Pathway of an Icosahedral Viral Capsid and Shape Determination of Two Successive Intermediates

Viral capsids derived from an icosahedral plant virus widely used in physical and nanotechnological investigations were fully dissociated into dimers by a rapid change of pH. The process was probed in vitro at high spatiotemporal resolution by time-resolved small-angle X-ray scattering using a high brilliance synchrotron source. A powerful custom-made global fitting algorithm allowed us to reconstruct the most likely pathway parametrized by a set of stoichiometric coefficients and to determine the shape of two successive intermediates by ab initio calculations. None of these two unexpected intermediates was previously identified in self-assembly experiments, which suggests that the disassembly pathway is not a mirror image of the assembly pathway. These findings shed new light on the mechanisms and the reversibility of the assembly/disassembly of natural and synthetic virus-based systems. They also demonstrate that both the structure and dynamics of an increasing number of intermediate species become accessible to experiments.

Common effect of chemical and external pressures on the magnetic properties of RCoPO (R = La, Pr, Nd, Sm). II.

We investigate the direct correspondence between Co band ferromagnetism and structural parameters in the pnictide oxides RCoPO for different rare-earth ions (R = La, Pr, Nd, Sm) by means of muon-spin spectroscopy and ab initio calculations, complementing our results published previously G. Prando et al., Common effect of chemical and external pressures on the magnetic properties of RCoPO (R = La, Pr), Phys. Rev. B 87, 064401 (2013)]. We find that both the transition temperature to the ferromagnetic phase T-C and the volume of the crystallographic unit cell V are conveniently tuned by the R ionic radius and/or external pressure. We report a linear correlation between T-C and V and our ab initio calculations unambiguously demonstrate a full equivalence of chemical and external pressures. As such, we show that R ions influence the ferromagnetic phase only via the induced structural shrinkage without involving any active role from the electronic f degrees of freedom, which are only giving a sizable magnetic contribution at much lower temperatures.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).

Synthesis, and crystal and electronic structure of sodium metal phosphate for use as a hybrid capacitor in non-aqueous electrolyte

Energy storage devices based on sodium have been considered as an alternative to traditional lithium based systems because of the natural abundance, cost effectiveness and low environmental impact of sodium. Their synthesis, and crystal and electronic properties have been discussed, because of the importance of electronic conductivity in supercapacitors for high rate applications. The density of states of a mixed sodium transition metal phosphate (maricite, NaMn1/3Co1/3Ni1/3PO4) has been determined with the ab initio generalized gradient approximation (GGA)+Hubbard term (U) method. The computed results for the mixed maricite are compared with the band gap of the parent NaFePO4 and the electrochemical experimental results are in good agreement. A mixed sodium transition metal phosphate served as an active electrode material for a hybrid supercapacitor. The hybrid device (maricite versus carbon) in a nonaqueous electrolyte shows redox peaks in the cyclic voltammograms and asymmetric profiles in the charge-discharge curves while exhibiting a specific capacitance of 40 F g(-1) and these processes are found to be quasi-reversible. After long term cycling, the device exhibits excellent capacity retention (95%) and coulombic efficiency (92%). The presence of carbon and the nanocomposite morphology, identified through X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) studies, ensures the high rate capability while offering possibilities to develop new cathode materials for sodium hybrid devices.