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.

Investigations on doping induced changes in structural, electronic structure and magnetic behavior of spintronic Cr-ZnS nanoparticles

Dilute magnetic semiconducting Zn1-xCrxS (x = 0.00, 0.01, 0.03, 0.05, 0.07) nanoparticles were synthesized by the co-precipitation technique using thioglycerol as the capping agent. Powder X-ray diffraction studies showed that Zn1-xCrxS nanoparticles exhibit zinc blende structure with no secondary phase, indicating that Cr ions are substituted at the Zn sites. Photoluminescence and Raman studies show the incorporation of Cr in ZnS nanoparticles. X-ray absorption studies depict that the valence of Zn remains unchanged and maintained in the divalent state, upon doping with Cr. The M-H curves at room temperature indicate the presence of weak ferromagnetism at room temperature due to structural defects. The increase in ferromagnetism with increasing Cr content up to 3%, demonstrates the possibility of tailoring the weak ferromagnetism in ZnS by appropriate Cr doping. (C) 2015 Elsevier Ltd. All rights reserved.

Electronic structure origin of conductivity and oxygen reduction activity changes in low-level Cr-substituted (La, Sr)MnO3

The electronic structure of the (La0.8Sr0.2)(0.98)Mn1-xCrxO3 model series (x = 0, 0.05, or 0.1) was measured using soft X-ray synchrotron radiation at room and elevated temperature. O K-edge near-edge X-ray absorption fine structure (NEXAFS) spectra showed that low-level chromium substitution of (La, Sr)MnO3 resulted in lowered hybridisation between O 2p orbitals and M 3d and M 4sp valance orbitals. Mn L-3-edge resonant photoemission spectroscopy measurements indicated lowered Mn 3d-O 2p hybridisation with chromium substitution. Deconvolution of O K-edge NEXAFS spectra took into account the effects of exchange and crystal field splitting and included a novel approach whereby the pre-peak region was described using the nominally filled t(2g) up arrow state. 10% chromium substitution resulted in a 0.17 eV lowering in the energy of the t(2g) up arrow state, which appears to provide an explanation for the 0.15 eV rise in activation energy for the oxygen reduction reaction, while decreased overlap between hybrid O 2p-Mn 3d states was in qualitative agreement with lowered electronic conductivity. An orbital-level understanding of the thermodynamically predicted solid oxide fuel cell cathode poisoning mechanism involving low-level chromium substitution on the B-site of (La, Sr)MnO3 is presented. (C) 2015 AIP Publishing LLC.

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.

Ab initio pair potentials at metal-ceramic interfaces

A systematic approach is proposed to obtain the interfacial interatomic potentials. By inverting ab initio adhesive energy curves for the metal-MgO ceramic interfaces, We derive interfacial potentials between Ag and O2-, Ag and Mg2+, Al and O2-, Al and Mg2+. The interfacial potentials, obtained from this method, demonstrate general features of bondings between metal atoms and ceramic ions.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Electronic structure and surface properties of SrBi2Ta2O9 and related oxides

The electronic structure of SrBi2Ta2O9 and related oxides such as SrBi2Nb2O9, Bi2WO6 and Bi3Ti4O12 have been calculated by the tight-binding method. In each case, the band gap is about 4.1 eV and the band edge states occur on the Bi-O layers and consist of mixed O p/Bi s states at the top of the valence band and Bi p states at the bottom of the conduction band. The main difference between the compounds is that Nb 5d and Ti 4d states in the Nb and Ti compounds lie lower than the Ta 6d states in the conduction band. The surface pinning levels are found to pin Schottky barriers 0.8 eV below the conduction band edge.