Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

Electronic Structure Evolution across the Peierls Metal-Insulator Transition in a Correlated Ferromagnet

Transition metal compounds often undergo spin-charge-orbital ordering due to strong electron-electron correlations. In contrast, low-dimensional materials can exhibit a Peierls transition arising from low-energy electron-phonon-coupling-induced structural instabilities. We study the electronic structure of the tunnel framework compound K2Cr8O16, which exhibits a temperature-dependent (T-dependent) paramagnetic-to-ferromagnetic- metal transition at T-C = 180 K and transforms into a ferromagnetic insulator below T-MI = 95 K. We observe clear T-dependent dynamic valence (charge) fluctuations from above T-C to T-MI, which effectively get pinned to an average nominal valence of Cr+3.75 (Cr4+:Cr3+ states in a 3:1 ratio) in the ferromagnetic-insulating phase. High-resolution laser photoemission shows a T-dependent BCS-type energy gap, with 2G(0) similar to 3.5(k(B)T(MI)) similar to 35 meV. First-principles band-structure calculations, using the experimentally estimated on-site Coulomb energy of U similar to 4 eV, establish the necessity of strong correlations and finite structural distortions for driving the metal-insulator transition. In spite of the strong correlations, the nonintegral occupancy (2.25 d-electrons/Cr) and the half-metallic ferromagnetism in the t(2g) up-spin band favor a low-energy Peierls metal-insulator transition.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Ab Initio MD Simulations of the Bronsted Acidity of Glutathione in Aqueous Solutions: Predicting pK(a) Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pK(a) value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides shift in the dissociation constant values as compared with the isolated accurate estimates of the pK(a) and correctly predicts the cysteine amino acid.

First principles DFT investigation Of Yttrium-decorated Boron-Nitride Nanotube: Electronic Structure and Hydrogen Storage

The electronic structure and hydrogen storage capability of Yttrium-doped BNNTs has been theoretically investigated using first principles density functional theory (DFT). Yttrium atom prefers the hollow site in the center of the hexagonal ring with a binding energy of 0.8048eV. Decorating by Y makes the system half-metallic and magnetic with a magnetic moment of 1.0 mu(B). Y decorated Boron-Nitride (8,0) nanotube can adsorb up to five hydrogen molecules whose average binding energy is computed as 0.5044eV. All the hydrogen molecules are adsorbed with an average desorption temperature of 644.708 K. Taking that the Y atoms can be placed only in alternate hexagons, the implied wt% comes out to be 5.31%, a relatively acceptable value for hydrogen storage materials. Thus, this system can serve as potential hydrogen storage medium.

Weak-coupling superconductivity in a strongly correlated iron pnictide

Iron-based superconductors have been found to exhibit an intimate interplay of orbital, spin, and lattice degrees of freedom, dramatically affecting their low-energy electronic properties, including superconductivity. Albeit the precise pairing mechanism remains unidentified, several candidate interactions have been suggested to mediate the superconducting pairing, both in the orbital and in the spin channel. Here, we employ optical spectroscopy (OS), angle-resolved photoemission spectroscopy (ARPES), ab initio band-structure, and Eliashberg calculations to show that nearly optimally doped NaFe0.978Co0.022As exhibits some of the strongest orbitally selective electronic correlations in the family of iron pnictides. Unexpectedly, we find that the mass enhancement of itinerant charge carriers in the strongly correlated band is dramatically reduced near the Gamma point and attribute this effect to orbital mixing induced by pronounced spin-orbit coupling. Embracing the true band structure allows us to describe all low-energy electronic properties obtained in our experiments with remarkable consistency and demonstrate that superconductivity in this material is rather weak and mediated by spin fluctuations.

Tailoring of Electronic Structure and Thermoelectric Properties of a Topological Crystalline Insulator by Chemical Doping

Topological crystalline insulators (TCIs) are a new quantum state of matter in which linearly dispersed metallic surface states are protected by crystal mirror symmetry. Owing to its vanishingly small bulk band gap, a TCI like Pb0.6Sn0.4Te has poor thermoelectric properties. Breaking of crystal symmetry can widen the band gap of TCI. While breaking of mirror symmetry in a TCI has been mostly explored by various physical perturbation techniques, chemical doping, which may also alter the electronic structure of TCI by perturbing the local mirror symmetry, has not yet been explored. Herein, we demonstrate that Na doping in Pb0.6Sn0.4Te locally breaks the crystal symmetry and opens up a bulk electronic band gap, which is confirmed by direct electronic absorption spectroscopy and electronic structure calculations. Na doping in Pb0.6Sn0.4Te increases p-type carrier concentration and suppresses the bipolar conduction (by widening the band gap), which collectively gives rise to a promising zT of 1 at 856 K for Pb0.58Sn0.40Na0.02Te. Breaking of crystal symmetry by chemical doping widens the bulk band gap in TCI, which uncovers a route to improve TCI for thermoelectric applications.

Investigations on structural, magnetic and electronic structure of Gd-doped ZnO nanostructures synthesized using sol-gel technique

GdxZn1-xO (x = 0, 0.02, 0.04 and 0.06) nanostructures have been synthesized using sol-gel technique and characterized to understand their structural and magnetic properties. X-ray diffraction (XRD) results show that Gd (0, 2, 4 and 6 %)-doped ZnO nanostructures crystallized in the wurtzite structure having space group C3(v) (P6(3)mc). Photoluminescence and Raman studies of Gd-doped ZnO powder show the formation of singly ionized oxygen vacancies. X-ray absorption spectroscopy reveals that Gd replaces the Zn atoms in the host lattice and maintains the crystal symmetry with slight lattice distortion. Gd L-3-edge spectra reveal charge transfer between Zn and Gd dopant ions. O K-edge spectra also depict the charge transfer through the oxygen bridge (Gd-O-Zn). Weak magnetic ordering is observed in all Gd-doped ZnO samples.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Effect of Structural Modification on the Quantum-Size Effect in II-VI Semiconducting Nanocrystals

The electronic structure of group II-VI semiconductors in the stable wurtzite form is analyzed using state-of-the-art ab initio approaches to extract a simple and chemically transparent tight-binding model. This model can be used to understand the variation in the bandgap with size, for nanoclusters of these compounds. Results complement similar information already available for same systems in the zinc blende structure. A comparison with all available experimental data on quantum size effects in group II-VI semiconductor nanoclusters establishes a remarkable agreement between theory and experiment in both structure types, thereby verifying the predictive ability of our approach. The significant dependence of the quantum size effect on the structure type suggests that the experimental bandgap change at a given size compared to the bulk bandgap, may be used to indicate the structural form of the nanoclusters, particularly in the small size limit, where broadening of diffraction features often make it difficult to unambiguously determine the structure.

Molecular structure-activity relationship study of some non-steroidal antiinflammatory agents using electrostatic potential mapping

A series of 6,11-dihydro-11-oxodibenz[b,e]oxepin-2-acetic acids (DOAA) which are known to be anti-inflammatory agents were studied. The geometries of some of the molecules obtained from X-ray crystallography were used in the calculations as such while the geometries of their derivatives were obtained by local, partial geometry optimization around the Sites of substitution employing the AMI method, keeping the remaining parts of the geometries the same as those in the parent molecules. Molecular electrostatic potential (MEP) mapping was performed for the molecules using optimized hybridization displacement charges (HDC) combined with Lowdin charges, as this charge distribution has been shown earlier to yield near ab initio quality results. A good correlation has been found between the MEP values near the oxygen atoms of the hydroxyl groups of the carboxy groups of the molecules and their anti-inflammatory activities. The result is broadly in agreement with the model proposed earlier by other authors regarding the structure-activity relationship for other similar molecules.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

Structure and stability of Li2F

The equilibrium geometries and fundamental vibration frequencies of the Li2F system were calculated by ab initio methods at the MP2 = full/6-311(+ +)G** and CCSD(T) levels. Two isomers were observed and are best described as salts of the Li-2(+) cation with F-. A linear isomer with an arrangement of atoms such as Li-Li-F and a bent C-2v structure are predicted. The stability of these structures are discussed in terms of charge resonance between Li and Li+. (C) 1999 Elsevier Science B.V. All rights reserved.

Studies on n=2 Aurivillius phases: structure of the series Bi3-xLaxTiNbO9 (0 <= x <= 1)

The crystal structures of the solid solutions of Bi3-xLaxTiNbO9 (0 less than or equal to x less than or equal to 1) have been analyzed by powder X-ray diffraction with supporting evidence from selected area electron diffraction (SAD). The structure of the starting member (x = 0) is verified to be in the orthorhombic space group A2(1) am while the end member (x = 1) is determined to crystallize in the centrosymmetric orthorhombic space group Pmcb. The structure of x = 1 phase is solved by ab initio powder diffraction. The intermediate compositions belong to the space group A2(1) am as confirmed by Rietveld refinements. Rietveld refinements on all the compositions reveal that the La3+ ion is disordered only in the A site and not in the [Bi2O2](2+) layer. The tilt in the Ti/NbO6 octahedra decreases with increasing x. (C) 2003 Elsevier B.V. All rights reserved.

Microwave spectrum and structure of C(6)H(5)CCH center dot center dot center dot H(2)S complex

Rotational spectra of C(6)H(5)CCH center dot center dot center dot H(2)S, C(6)H(5)CCH center dot center dot center dot H(2)(34)S, C(6)H(5)CCH center dot center dot center dot HDS, C(6)H(5)CCH center dot center dot center dot D(2)S and C(6) H(5)CCD center dot center dot center dot H(2)S complexes have been observed using a pulsed nozzle Fourier transform microwave spectrometer. The observed spectrum is consistent with a structure in which hydrogen sulfide is located over the phenyl ring pi cloud and the distance between the centers of masses of the two monomers is 3.74 +/- 0.01 angstrom. In the complex, the H(2)S unit is shifted from the phenyl ring center towards the acetylene group. The vibrationally averaged structure has an effective Cs symmetry. Ab initio calculations were performed at MP2/aug-cc-pVDZ level of theory to locate the possible geometries of the complex. The calculations reveal the experimentally observed structure to be more stable than a coplanar arrangement of the monomers, which was observed for the C(6)H(5)CCH center dot center dot center dot H(2)O complex. Atoms in molecule theoretical analysis shows the presence of S-H center dot center dot center dot pi hydrogen bond. For the parent isotopologue, each transition frequency was found to split into two resulting from an interchange of the equivalent hydrogens of H(2)S unit in the complex. (C) 2011 Elsevier Inc. All rights reserved.