Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

First hyperpolarizability of Ru-half-sandwich complexes: The effect of halogen atom substitution on the ancillary ligand

The first hyperpolarizability (beta) of a series of half-sandwich Ru complexes with a mercaptobenzothiazole ligand bearing a halogen atom substitution in the para-position has been investigated by hyper-Rayleigh scattering and quantum chemical calculations. The heterocyclic ligand with a bromine atom in the para position makes it a very good donor and charge flows to the Ru center enhancing the beta value of the complex by a factor of 2 compared to the complex with the ligand without the halogen substitution. The resonance (+R) and the inductive (-I) effects exerted by the halogen atom in the para position push electrons in opposing directions in the complex. For the Br and Cl atoms the resonance effect dominates which enables the ligand to donate electrons to the metal center thereby increasing the hyperpolarizability whereas for the fluorine atom, the inductive effect is dominant which reduces the charge flow to the metal and the hyperpolarizability drops even below that of the unsubstituted ligand. This unprecedented halogen atom effect on beta of metal complexes is reported. (C) 2015 Elsevier By. All rights reserved.

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.

Identification of local variations within secondary structures of proteins

Secondary-structure elements (SSEs) play an important role in the folding of proteins. Identification of SSEs in proteins is a common problem in structural biology. A new method, ASSP (Assignment of Secondary Structure in Proteins), using only the path traversed by the C atoms has been developed. The algorithm is based on the premise that the protein structure can be divided into continuous or uniform stretches, which can be defined in terms of helical parameters, and depending on their values the stretches can be classified into different SSEs, namely -helices, 3(10)-helices, -helices, extended -strands and polyproline II (PPII) and other left-handed helices. The methodology was validated using an unbiased clustering of these parameters for a protein data set consisting of 1008 protein chains, which suggested that there are seven well defined clusters associated with different SSEs. Apart from -helices and extended -strands, 3(10)-helices and -helices were also found to occur in substantial numbers. ASSP was able to discriminate non--helical segments from flanking -helices, which were often identified as part of -helices by other algorithms. ASSP can also lead to the identification of novel SSEs. It is believed that ASSP could provide a better understanding of the finer nuances of protein secondary structure and could make an important contribution to the better understanding of comparatively less frequently occurring structural motifs. At the same time, it can contribute to the identification of novel SSEs. A standalone version of the program for the Linux as well as the Windows operating systems is freely downloadable and a web-server version is also available at .

Online_DPI: a web server to calculate the diffraction precision index for a protein structure

An online computing server, Online_DPI (where DPI denotes the diffraction precision index), has been created to calculate the `Cruickshank DPI' value for a given three-dimensional protein or macromolecular structure. It also estimates the atomic coordinate error for all the atoms available in the structure. It is an easy-to-use web server that enables users to visualize the computed values dynamically on the client machine. Users can provide the Protein Data Bank (PDB) identification code or upload the three-dimensional atomic coordinates from the client machine. The computed DPI value for the structure and the atomic coordinate errors for all the atoms are included in the revised PDB file. Further, users can graphically view the atomic coordinate error along with `temperature factors' (i.e. atomic displacement parameters). In addition, the computing engine is interfaced with an up-to-date local copy of the Protein Data Bank. New entries are updated every week, and thus users can access all the structures available in the Protein Data Bank. The computing engine is freely accessible online at http://cluster.physics.iisc.ernet.in/dpi/.

Theoretical Insights to Niobium-Doped Monolayer MoS2-Gold Contact

We report a first principles study of the electronic properties for a contact formed between Nb-doped monolayer MoS2 and gold for different doping concentrations. We first focus on the shift of energy levels in band structure and the density of states with respect to the Fermi level for a geometrically optimized 5 x 5 MoS2 supercell for both pristine and Nb-doped structures. The doping is achieved by substituting Mo atoms with Nb atoms at random positions. It is observed that for an experimentally reported sheet hole doping concentration of (rho(2D)) 1.8 x 10(14) cm(-2), the pristine MoS2 converts to degenerate p-type semiconductor. Next, we interface this supercell with six layers of < 111 > cleaved surface of gold to investigate the contact nature of MoS2-Au system. By careful examination of projected band structure, projected density of states, effective potential and charge density difference, we demonstrate that the Schottky barrier nature observed for pure MoS2-Au contact can be converted from n-type to p-type by efficient Nb doping.

Remarkable Effect of Chalcogen Substitution on an Enzyme Mimetic for Deiodination of Thyroid Hormones

Iodothyronine deiodinases are selenoenzymes which regulate the thyroid hormone homeostasis by catalyzing the regioselective deiodination of thyroxine (T4). Synthetic deiodinase mimetics are important not only to understand the mechanism of enzyme catalysis, but also to develop therapeutic agents as abnormal thyroid hormone levels have implications in different diseases, such as hypoxia, myocardial infarction, critical illness, neuronal ischemia, tissue injury, and cancer. Described herein is that the replacement of sulfur/selenium atoms in a series of deiodinase mimetics by tellurium remarkably alters the reactivity as well as regioselectivity toward T4. The tellurium compounds reported in this paper represent the first examples of deiodinase mimetics which mediate sequential deiodination of T4 to produce all the hormone derivatives including T0 under physiologically relevant conditions.

Do we see what we should see? Describing non-covalent interactions in protein structures including precision

The power of X-ray crystal structure analysis as a technique is to `see where the atoms are'. The results are extensively used by a wide variety of research communities. However, this `seeing where the atoms are' can give a false sense of security unless the precision of the placement of the atoms has been taken into account. Indeed, the presentation of bond distances and angles to a false precision (i.e. to too many decimal places) is commonplace. This article has three themes. Firstly, a basis for a proper representation of protein crystal structure results is detailed and demonstrated with respect to analyses of Protein Data Bank entries. The basis for establishing the precision of placement of each atom in a protein crystal structure is non-trivial. Secondly, a knowledge base harnessing such a descriptor of precision is presented. It is applied here to the case of salt bridges, i.e. ion pairs, in protein structures; this is the most fundamental place to start with such structure-precision representations since salt bridges are one of the tenets of protein structure stability. Ion pairs also play a central role in protein oligomerization, molecular recognition of ligands and substrates, allosteric regulation, domain motion and alpha-helix capping. A new knowledge base, SBPS (Salt Bridges in Protein Structures), takes these structural precisions into account and is the first of its kind. The third theme of the article is to indicate natural extensions of the need for such a description of precision, such as those involving metalloproteins and the determination of the protonation states of ionizable amino acids. Overall, it is also noted that this work and these examples are also relevant to protein three-dimensional structure molecular graphics software.

The X-C---Y (X = O/F, Y = O, S, F, Cl, Br, N, P,) Carbon bond and hydrophobic interactions

While the tetrahedral face of methane has an electron rich centre and can act as a hydrogen bond acceptor, substitution of one of its hydrogens with some electron withdrawing group (such as -F/OH) can make the opposite face electron deficient. Electrostatic potential calculations confirm this and high level quantum calculations show interactions between the positive face of methanol/methyl fluoride and electron rich centers of other molecules such as H2O. Analysis of the wave functions of atoms in molecules shows the presence of an unusual C···Y interaction, which could be called 'carbon bonding'. NBO analysis and vibrational frequency shifts confirm the presence of this interaction. Given the properties of alkyl groups bonded to electronegative elements in biological molecules, such interactions could play a significant role, which is yet to be recognized. This and similar interactions could give an enthalpic contribution to what is called the 'hydrophobic interactions'.

Magnetometry using Ramsey interferometry in a Yb atomic beam

We use the Ramsey separated oscillatory fields technique in a 400 degrees C thermal beam of ytterbium (Yb) atoms to measure the Larmor precession frequency (and hence the magnetic field) with high precision. For the experiment, we use the strongly allowed S-1(0) P-1(1) transition at 399 nm, and choose the odd isotope Yb-171 with nuclear spin I = 1/2, so that the ground state has only two magnetic sublevels m(F) = +/- 1/2. With a magnetic field of 22.2 G and a separation of about 400 mm between the oscillatory fields, the central Ramsey fringe is at 16.64 kHz and has a width of 350 Hz. The technique can be readily adapted to a cold atomic beam, which is expected to give more than an order-of-magnitude improvement in precision. The signal-to-noise ratio is comparable to other techniques of magnetometry; therefore it should be useful for all kinds of precision measurements such as searching for a permanent electric dipole moment in atoms.

A new class of high strength high temperature Cobalt based gamma-gamma ` Co-Mo-Al alloys stabilized with Ta addition

The present paper reports a new class of Co based superalloys that has gamma-gamma' microstructure and exhibits much lower density compared to other commercially available Co superalloys including Co-Al-W based alloys. The basic composition is Co-10Al-5Mo (at%) with addition of 2 at% Ta for stabilization of gamma' phase. The gamma-gamma' microstructure evolves through solutionising and aging treatment. Using first principles calculations, we observe that Ta plays a crucial role in stabilizing gamma' phase. By addition of Ta in the basic stoichiometric composition Co-3(Al, Mo), the enthalpy of formation (Delta H-f) of L1(2) structure (gamma' phase) becomes more negative in comparison to DO19 structure. The All of the L12 structure becomes further more negative by the occupancy of Ni and Ti atoms in the lattice suggesting an increase in the stability of the gamma' precipitates. Among large number of alloys studied experimentally, the paper presents results of detailed investigations on Co-10Al-5Mo-2Ta, Co-30Ni-10Al-5Mo-2Ta and Co-30Ni-10Al-5Mo-2Ta-2Ti. To evaluate the role alloying elements, atom probe tomography investigations were carried out to obtain partition coefficients for the constituent elements. The results show strong partitioning of Ni, Al, Ta and Ti in ordered gamma' precipitates. 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

An enolato-bridged dinuclear Cu(II) complex with a coumarin-assisted precursor: a spectral, magnetic and biological study

A new, phenoxo-bridged Cu-II dinuclear complex Cu-2(L)(2)(DMF)(2)] (1) has been obtained by employing the coumarin-assisted tridentate precursor, H2L, benzoic acid(7-hydroxy-4-methyl-2-oxo-2H-chromen-8-ylmethylene)-hydrazide]. Complex 1 has been systematically characterized by FTIR, UV-Vis, fluorescence and PR spectrometry. The single crystal X-ray diffraction analysis of 1 shows that the geometry around each copper ion is square pyramidal, comprising two enolato oxygen atoms belonging to different ligands (which assemble the dimer bridging the two metal centers), one imine-N and one phenolic-O atoms of the Schiff base and one oxygen atom from the DMF molecule. The temperature dependent magnetic interpretation agrees with the existence of weak ferromagnetic interactions between the bridging dinuclear Cu(II) ions. Both the ligand and complex 1 exhibit anti-mycobacterial activity and considerable efficacy towards M. tuberculosis H37Rv ATCC 27294 and M. tuberculosis H37Ra ATCC 25177 strains. The cytotoxicity study on human adenocarcinoma cell lines (MCF7) suggests that the ligand and complex 1 have potential anticancer properties. Molecular docking of H2L with the enoyl acyl carrier protein reductase of M. tuberculosis H37R(v) (PDB ID: 4U0K) is examined and the best docked pose of H2L shows one hydrogen bond with Thr196 (1.99 angstrom).

A Split-and-Merge Dictionary Learning Algorithm for Sparse Representation: Application to Image Denoising

In big data image/video analytics, we encounter the problem of learning an over-complete dictionary for sparse representation from a large training dataset, which cannot be processed at once because of storage and computational constraints. To tackle the problem of dictionary learning in such scenarios, we propose an algorithm that exploits the inherent clustered structure of the training data and make use of a divide-and-conquer approach. The fundamental idea behind the algorithm is to partition the training dataset into smaller clusters, and learn local dictionaries for each cluster. Subsequently, the local dictionaries are merged to form a global dictionary. Merging is done by solving another dictionary learning problem on the atoms of the locally trained dictionaries. This algorithm is referred to as the split-and-merge algorithm. We show that the proposed algorithm is efficient in its usage of memory and computational complexity, and performs on par with the standard learning strategy, which operates on the entire data at a time. As an application, we consider the problem of image denoising. We present a comparative analysis of our algorithm with the standard learning techniques that use the entire database at a time, in terms of training and denoising performance. We observe that the split-and-merge algorithm results in a remarkable reduction of training time, without significantly affecting the denoising performance.

Microwave spectroscopic and theoretical investigations of the strongly hydrogen bonded hexafluoroisopropanol...water complex

This paper reports microwave spectroscopic and theoretical investigations on the interaction of water with hexafluoroisopropanol (HFIP). The HFIP monomer can exist in two conformations, antiperiplanar (AP) and synclinical (SC). The former is about 5 kJ mol(-1) more stable than the latter. Theoretical calculations predicted three potential minima for the complex, two having AP and one having SC conformations. Though, the binding energy for the HFIP(SC)...H2O turned out to be larger than that for the other two conformers having HFIP in the AP form, the global minimum for the complex in the potential energy hypersurface had HFIP in the AP form. Experimental rotational constants for four isotopologues measured using a pulsed nozzle Fourier transform microwave spectrometer, correspond to the global minimum in the potential energy hypersurface. The structural parameters and the internal dynamics of the complex could be determined from the rotational spectra of the four isotopologues. The global minimum has the HFIP(AP) as a hydrogen bond donor forming a strong hydrogen bond with H2O. To characterize the strength of the bonding and to probe the other interactions within the complex, atoms in molecules, non-covalent interaction index and natural bond orbital theoretical analyses have been performed.

Study of a four-level system in vee plus ladder configuration

We present the results of a theoretical study of a four-level atomic system in vee + ladder configuration using a density matrix analysis. The absorption and dispersion profiles are derived for a weak probe field and for varying strengths of the two strong control fields. For specificity, we choose energy levels of Rb-87, and present results for both stationary atoms and moving atoms in room temperature vapor. An electromagnetically induced absorption (EIA) peak with negative dispersion is observed at zero probe de-tuning when the control fields have equal strengths, which switches to electromagnetically induced transparency (ET) with positive dispersion (due to splitting of the EIA peak) when the control fields are unequal. There is significant linewidth narrowing in thermal vapor. (C) 2015 Elsevier B.V. All rights reserved.