Efficient methods for the quantum chemical treatment of protein structures: The effects of London-dispersion and basis-set incompleteness on peptide and water-cluster geometries

We demonstrate how quantum chemical Hartree-Fock (HF) or density functional theory (DFT) optimizations with small basis sets of peptide and water cluster structures are decisively improved if London-dispersion effects, the basis-set-superposition error (BSSE), and other basis-set incompleteness errors are addressed. We concentrate on three empirical corrections to these problems advanced by Grimme and co-workers that lead to computational strategies that are both accurate and efficient. Our analysis encompasses a reoptimized version of Hobza's P26 set of tripeptide structures, a new test set of conformers of cysteine dimers, and isomers of the water hexamer. These systems reflect features commonly found in protein crystal structures. In all cases, we recommend Grimme's DFT-D3 correction for London-dispersion. We recommend usage of large basis sets such as cc-pVTZ whenever possible to reduce any BSSE effects and, if this is not possible, to use Grimme's gCP correction to account for BSSE when small basis sets are used. We demonstrate that S-S and C-S bond lengths are very prone to basis-set incompleteness and that polarization functions should always be used on S atoms. At the double-ζ level, the PW6B95-D3-gCP DFT method combined with the SVP and 6-31G* basis sets yields accurate results. Alternatively, the HF-D3-gCP/SV method is recommended, with inclusion of polarization functions for S atoms only. Minimal basis sets offer an intriguing route to highly efficient calculations, but due to significant basis-set incompleteness effects, calculated bond lengths are seriously overestimated, making applications to large proteins very difficult, but we show that Grimme's newest HF-3c correction overcomes this problem and so makes this computational strategy very attractive. Our results provide a useful guideline for future applications to the optimization, quantum refinement, and dynamics of large proteins. © 2013 American Chemical Society

Problems, successes and challenges for the application of dispersion-corrected density-functional theory combined with dispersion-based implicit solvent models to large-scale hydrophobic self-assembly and polymorphism

© 2015 Taylor & Francis. The recent advent of dispersion-corrected density-functional theory (DFT) methods allows for quantitative modelling of molecular self-assembly processes, and we consider what is required to develop applications to the formation of large self-assembled monolayers (SAMs) on hydrophobic surfaces from organic solution. Focus is on application of the D3 dispersion correction of Grimme combined with the solvent dispersion model of Floris, Tomasi and Pascual-Ahuir to simulate observed scanning-tunnelling microscopy (STM) images of various polymorphs of tetraalkylporphyrin SAMs on highly oriented pyrolytic graphite surfaces. The most significant problem is identified as the need to treat SAM structures that are incommensurate with those of the substrate, providing a challenge to the use of traditional periodic-imaging boundary techniques. Using nearby commensurate lattices introduces non-systematic errors into calculated lattice constants and free energies of SAM formation that are larger than experimental uncertainties and polymorph differences. Developing non-periodic methods for polymorph interface simulation also remains a challenge. Despite these problems, existing methods can be used to interpret STM images and SAM atomic structures, distinguishing between multiple feasible polymorph types. They also provide critical insight into the factors controlling polymorphism. All this stems from a delicate balance that the intermolecular D3 and solvent Floris, Tomasi and Pascual-Ahuir corrections provide. Combined optimised treatments should yield fully quantitative approaches in the future

Free energies for the coordination of ligands to the magnesium of chlorophyll-a in solvents

© 2015 Taylor & Francis. The coordination of bases to chlorophyll magnesium modifies spectroscopic properties in solution as well as in situ in reaction centres. We evaluate the free energies of complexation of one or two pyridine, 1-propanol, diethyl ether or water solvent molecules at 298 and 150 K to rationalise observed phenomena. Various a priori dispersion-corrected density functional theory calculations are performed as well as second-order Møller-Plesset calculations, focusing on the effects of dispersion modifying the intermolecular interactions, of dispersion modifying solvation energies, of entropy, and of basis-set superposition error. A process of particular interest is magnesium complexation in ether at low temperature that is often exploited to assign the Q-band visible absorption spectrum of chlorophyll. Recently, we demonstrated that trace water interferes with this process, but the nature of the resulting complex could not be uniquely determined; here, it is identified as ether.Chlorophyll-a.H2O, consistent with interpretations based on our authoritative 2013 assignment

Getting Past the Language Gap: Innovations in Machine Translation

In this chapter, we will be reviewing state of the art machine translation systems, and will discuss innovative methods for machine translation, highlighting the most promising techniques and applications. Machine translation (MT) has benefited from a revitalization in the last 10 years or so, after a period of relatively slow activity. In 2005 the field received a jumpstart when a powerful complete experimental package for building MT systems from scratch became freely available as a result of the unified efforts of the MOSES international consortium. Around the same time, hierarchical methods had been introduced by Chinese researchers, which allowed the introduction and use of syntactic information in translation modeling. Furthermore, the advances in the related field of computational linguistics, making off-the-shelf taggers and parsers readily available, helped give MT an additional boost. Yet there is still more progress to be made. For example, MT will be enhanced greatly when both syntax and semantics are on board: this still presents a major challenge though many advanced research groups are currently pursuing ways to meet this challenge head-on. The next generation of MT will consist of a collection of hybrid systems. It also augurs well for the mobile environment, as we look forward to more advanced and improved technologies that enable the working of Speech-To-Speech machine translation on hand-held devices, i.e. speech recognition and speech synthesis. We review all of these developments and point out in the final section some of the most promising research avenues for the future of MT

Ultraviolet photodissociation action spectroscopy of the N-pyridinium cation

© 2015 AIP Publishing LLC. The S1←S0 electronic transition of the N-pyridinium ion (C5H5NH+) is investigated using ultraviolet photodissociation (PD) spectroscopy of the bare ion and also the N2-tagged complex. Gas-phase N-pyridinium ions photodissociate by the loss of molecular hydrogen (H2) in the photon energy range 37 000-45 000 cm-1 with structurally diagnostic ion-molecule reactions identifying the 2-pyridinylium ion as the exclusive co-product. The photodissociation action spectra reveal vibronic details that, with the aid of electronic structure calculations, support the proposal that dissociation occurs through an intramolecular rearrangement on the ground electronic state following internal conversion. Quantum chemical calculations are used to analyze the measured spectra. Most of the vibronic features are attributed to progressions of totally symmetric ring deformation modes and out-of-plane modes active in the isomerization of the planar excited state towards the non-planar excited state global minimum

Electron-vibration entanglement in the Born-Oppenheimer description of chemical reactions and spectroscopy

This journal is © the Owner Societies. Entanglement is sometimes regarded as the quintessential measure of the quantum nature of a system and its significance for the understanding of coupled electronic and vibrational motions in molecules has been conjectured. Previously, we considered the entanglement developed in a spatially localized diabatic basis representation of the electronic states, considering design rules for qubits in a low-temperature chemical quantum computer. We extend this to consider the entanglement developed during high-energy processes. We also consider the entanglement developed using adiabatic electronic basis, providing a novel way for interpreting effects of the breakdown of the Born-Oppenheimer (BO) approximation. We consider: (i) BO entanglement in the ground-state wavefunction relevant to equilibrium thermodynamics, (ii) BO entanglement associated with low-energy wavefunctions relevant to infrared and tunneling spectroscopies, (iii) BO entanglement in high-energy eigenfunctions relevant to chemical reaction processes, and (iv) BO entanglement developed during reactive wavepacket dynamics. A two-state single-mode diabatic model descriptive of a wide range of chemical phenomena is used for this purpose. The entanglement developed by BO breakdown correlates simply with the diameter of the cusp introduced by the BO approximation, and a hierarchy appears between the various BO-breakdown correction terms, with the first-derivative correction being more important than the second-derivative correction which is more important than the diagonal correction. This simplicity is in contrast to the complexity of BO-breakdown effects on thermodynamic, spectroscopic, and kinetic properties. Further, processes poorly treated at the BO level that appear adequately treated using the Born-Huang adiabatic approximation are found to have properties that can only be described using a non-adiabatic description. For the entanglement developed between diabatic electronic states and the nuclear motion, qualitatively differently behavior is found compared to traditional properties of the density matrix and hence entanglement provides new information about system properties. For chemical reactions, this type of entanglement simply builds up as the transition-state region is crossed. It is robust to small changes in parameter values and is therefore more attractive for making quantum qubits than is the related fragile ground-state entanglement, provided that coherent motion at the transition state can be sustained

Auditory presentation and synchronization in Adobe Flash and HTML5/JavaScript Web experiments

Substantial recent research has examined the accuracy of presentation durations and response time measurements for visually presented stimuli in Web-based experiments, with a general conclusion that accuracy is acceptable for most kinds of experiments. However, many areas of behavioral research use auditory stimuli instead of, or in addition to, visual stimuli. Much less is known about auditory accuracy using standard Web-based testing procedures. We used a millisecond-accurate Black Box Toolkit to measure the actual durations of auditory stimuli and the synchronization of auditory and visual presentation onsets. We examined the distribution of timings for 100 presentations of auditory and visual stimuli across two computers with difference specs, three commonly used browsers, and code written in either Adobe Flash or JavaScript. We also examined different coding options for attempting to synchronize the auditory and visual onsets. Overall, we found that auditory durations were very consistent, but that the lags between visual and auditory onsets varied substantially across browsers and computer systems

Synthetically tuneable biomimetic artificial photosynthetic reaction centres that closely resemble the natural system in purple bacteria

© 2016 The Royal Society of Chemistry. Porphyrin-based photosynthetic reaction centre (PRC) mimics, ZnPQ-Q2HP-C60 and MP2Q-Q2HP-C60 (M = Zn or 2H), designed to have a similar special-pair electron donor and similar charge-separation distances, redox processes and photochemical reaction rates to those in the natural PRC from purple bacteria, have been synthesised and extensive photochemical studies performed. Mechanisms of electron-transfer reactions are fully investigated using femtosecond and nanosecond transient absorption spectroscopy. In benzonitrile, all models show picosecond-timescale charge-separations and the final singlet charge-separations with the microsecond-timescale. The established lifetimes are long compared to other processes in organic solar cells or other organic light harvesting systems. These rigid, synthetically flexible molecules provide the closest mimics to the natural PRC so far synthesised and present a future direction for the design of light harvesters with controllable absorption, redox, and kinetics properties

HE 0557-4840 - Ultra-Metal-Poor and Carbon-Rich

We report the discovery and high-resolution, high S/N, spectroscopic analysis of the ultra-metal-poor red giant HE 0557-4840, which is the third most heavy-element deficient star currently known. Its atmospheric parameters are T_eff = 4900 K, log g = 2.2, and [Fe/H]= -4.75. This brings the number of stars with [Fe/H] < -4.0 to three, and the discovery of HE 0557-4840 suggests that the metallicity distribution function of the Galactic halo does not have a "gap" between [Fe/H] = -4.0, where several stars are known, and the two most metal-poor stars, at [Fe/H] ~ -5.3. HE 0557-4840 is carbon rich - [C/Fe] = +1.6 - a property shared by all three objects with [Fe/H] < -4.0, suggesting that the well-known increase of carbon relative to iron with decreasing [Fe/H] reaches its logical conclusion - ubiquitous carbon richness - at lowest abundance. We also present abundances (nine) and limits (nine) for a further 18 elements. For species having well-measured abundances or strong upper limits, HE 0557-4840 is "normal" in comparison with the bulk of the stellar population at [Fe/H] ~ -4.0 - with the possible exception of Co. We discuss the implications of these results for chemical enrichment at the earliest times, in the context of single ("mixing and fallback") and two-component enrichment models. While neither offers a clear solution, the latter appears closer to the mark. Further data are required to determine the oxygen abundance and improve that of Co, and hence more strongly constrain the origin of this object.Comment: Submitted to Astrophysical Journal. 52 pages (41 text, 11 figures

Relating transition-state spectroscopy to standard chemical spectroscopic processes

© 2017 Elsevier B.V. Transition-state spectra are mapped out using generalized adiabatic electron-transfer theory. This simple model depicts diverse chemical properties, from aromaticity, through bound reactions such as isomerizations and atom-transfer processes with classic transition states, to processes often described as being “non-adiabatic”, to those in the “inverted” region that become slower as they are made more exothermic. Predictably, the Born-Oppenheimer approximation is found inadequate for modelling transition-state spectra in the weak-coupling limit. In this limit, the adiabatic Born-Huang approximation is found to perform much better than non-adiabatic surface-hopping approaches. Transition-state spectroscopy is shown to involve significant quantum entanglement between electronic and nuclear motion