Oscillating electric-field effects on adsorbed-water at rutile- and anatase-TiO2 surfaces

We have performed non-equilibrium molecular dynamics simulations of various TiO2/water interfaces at ambient temperature in presence of oscillating electric fields in frequency range 20–100 GHz and RMS intensities 0.05–0.25 V/Å. Although the externally applied fields are by one order of magnitude lower than the intrinsic electric field present on the interfaces (∼1.5–4.5 V/Å), significant non-thermal coupling of rotational and translational motion of water molecules was clearly observed. Enhancement of the motion, manifested by increase of diffusivity, was detected in the first hydration layer, which is known to be heavily confined by adsorption to the TiO2 surface. Interestingly, the diffusivity increases more rapidly on anatase than on rutile facets where the adsorbed water was found to be more organized and restrained. We observed that the applied oscillating field reduces number of hydrogen bonds on the interface. The remaining H-bonds are weaker than those detected under zero-field conditions; however, their lifetime increases on most of the surfaces when the low-frequency fields are applied. Reduction of adsorption interaction was observed also in IR spectra of interfacial water where the directional patterns are smeared as the intensities of applied fields increase

Communication: Influence of external static and alternating electric fields on water from long-time non-equilibrium ab initio molecular dynamics

The response of water to externally applied electric fields is of central relevance in the modern world, where many extraneous electric fields are ubiquitous. Historically, the application of external fields in non-equilibrium molecular dynamics has been restricted, by and large, to relatively inexpensive, more or less sophisticated, empirical models. Here, we report long-time non-equilibrium ab initio molecular dynamics in both static and oscillating (time-dependent) external electric fields, therefore opening up a new vista in rigorous studies of electric-field effects on dynamical systems with the full arsenal of electronic-structure methods. In so doing, we apply this to liquid water with state-ofthe-art non-local treatment of dispersion, and we compute a range of field effects on structural and dynamical properties, such as diffusivities and hydrogen-bond kinetics. Published by AIP Publishing

Perturbation of hydration layer in solvated proteins by external electric and electromagnetic fields: Insights from non-equilibrium molecular dynamics

Given the fundamental role of water in governing the biochemistry of enzymes, and in regulating their wider biological activity (e.g., by local water concentration surrounding biomolecules), the influence of extraneous electric and electromagnetic (e/m) fields thereon is of central relevance to biophysics and, more widely, biology. With the increase in levels of local and atmospheric microwave-frequency radiation present in modern life, as well as other electric-field exposure, the impact upon hydration-water layers surrounding proteins, and biomolecules generally, becomes a particularly pertinent issue. Here, we present a (non-equilibrium) molecular-dynamics-simulation study on a model protein (hen egg-white lysozyme) hydrated in water, in which we determine, inter alia, translational self-diffusivities for both hen egg-white lysozyme and its hydration layer together with relaxation dynamics of the hydrogen-bond network between the protein and its hydration-layer water molecules on a residue-per-residue basis. Crucially, we perform this analysis both above and below the dynamical-transition temperature (at ∼220 K), at 300 and 200 K, respectively, and we compare the effects of external static-electric and e/m fields with linear-response-régime (r.m.s.) intensities of 0.02 V/Å. It was found that the translational self-diffusivity of hen egg-white lysozyme and its hydration-water layer are increased substantially in static fields, primarily due to the induced electrophoretic motion, whilst the water-protein hydrogen-bond-network-rearrangement kinetics can also undergo rather striking accelerations, primarily due to the enhancement of a larger-amplitude local translational and rotational motion by charged and dipolar residues, which serves to promote hydrogen-bond breakage and re-formation kinetics. These external-field effects are particularly evident at 200 K, where they serve to induce the protein- and solvation-layer-response effects redolent of dynamical transition at a lower temperature (∼200 K) vis-à-vis the zero-field case (∼220 K)

Comparison of hydration reactions for "piano-stool" RAPTA-B and [Ru(eta(6) - arene)(en)Cl](+) complexes: Density functional theory computational study

The hydration process for two Ru(II) representative half-sandwich complexes: Ru(arene)(pta)Cl2 (from the RAPTA family) and [Ru(arene)(en)Cl]+ (further labeled as Ru_en) were compared with analogous reaction of cisplatin. In the study, quantum chemical methods were employed. All the complexes were optimized at the B3LYP/6-31G(d) level using Conductor Polarizable Continuum Model (CPCM) solvent continuum model and single-point (SP) energy calculations and determination of electronic properties were performed at the B3LYP/6-311++G(2df,2pd)/CPCM level. It was found that the hydration model works fairly well for the replacement of the first chloride by water where an acceptable agreement for both Gibbs free energies and rate constants was obtained. However, in the second hydration step worse agreement of the experimental and calculated values was achieved. In agreement with experimental values, the rate constants for the first step can be ordered as RAPTA-B > Ru_en > cisplatin. The rate constants correlate well with binding energies (BEs) of the Pt/Ru–Cl bond in the reactant complexes. Substitution reactions on Ru_en and cisplatin complexes proceed only via pseudoassociative (associative interchange) mechanism. On the other hand in the case of RAPTA there is also possible a competitive dissociation mechanism with metastable pentacoordinated intermediate. The first hydration step is slightly endothermic for all three complexes by 3–5 kcal/mol. Estimated BEs confirm that the benzene ligand is relatively weakly bonded assuming the fact that it occupies three coordination positions of the Ru(II) cation

Study of hydrogen-molecule guests in type II clathrate hydrates using a force-matched potential model parameterised from ab initio molecular dynamics

The force-matching method has been applied to parameterise an empirical potential model for water-water and water-hydrogen intermolecular interactions for use in clathrate-hydrate simulations containing hydrogen guest molecules. The underlying reference simulations constituted ab initio molecular dynamics (AIMD) of clathrate hydrates with various occupations of hydrogen-molecule guests. It is shown that the resultant model is able to reproduce AIMD-derived free-energy curves for the movement of a tagged hydrogen molecule between the water cages that make up the clathrate, thus giving us confidence in the model. Furthermore, with the aid of an umbrella-sampling algorithm, we calculate barrier heights for the force-matched model, yielding the free-energy barrier for a tagged molecule to move between cages. The barrier heights are reasonably large, being on the order of 30 kJ/mol, and are consistent with our previous studies with empirical models [C. J. Burnham and N. J. English, J. Phys. Chem. C 120, 16561 (2016) and C. J. Burnham et al., Phys. Chem. Chem. Phys. 19, 717 (2017)]. Our results are in opposition to the literature, which claims that this system may have very low barrier heights. We also compare results to that using the more ad hoc empirical model of Alavi et al. [J. Chem. Phys. 123, 024507 (2005)] and find that this model does very well when judged against the force-matched and ab initio simulation data

HAB79: A New Molecular Dataset for Benchmarking DFT and DFTB Electronic Couplings Against High-Level Ab-initio Calculations

A new molecular dataset called HAB79 is introduced to provide ab-initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is comprised of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUE) of 27.7% (26.3%) for the GGA functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost effective POD in combination with a GGA functional (PBE) and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude

Near-Shore Aggregation Mechanism of Electrolyte Decomposition Products to Explain Solid Electrolyte Interphase Formation

To get insight of the formation mechanism of solid electrolyte interphase (SEI) film in Lithium-ion battery (LIB), we examine a probable scenario, referred to as “surface growth mechanism,” for electrolyte involving ethylene carbonate (EC) solvent and vinylene carbonate (VC) additive by using density functional theory (DFT). We first extracted stable SEI film components (SFCs) for the EC/VC electrolyte and constructed probable SFC aggregates via DFT molecular dynamics. We then examined their solubility in the EC solution, their adhesion to a model graphite electrode, and the electronic properties. The results showed that the SFC aggregates are characterized by “unstable adhesion” to the graphite surface and “high electronic insulation” against the EC solution. These characteristics preclude explaining SEI growth up to a typical thickness of several tens of nanometers based on the surface growth mechanism. With the present results, we propose “near-shore aggregation” mechanism, where the SFCs formed at the electrode surface desorb into the near-shore region and form aggregates. The SFC aggregates coalesce and come into contact with the electrode to complete the SEI formation. The present model provides a novel perspective for the long-standing problem of SEI formation

Kinetics of trifurcated electron flow in the decaheme bacterial proteins MtrC and MtrF

The bacterium Shewanella oneidensis has evolved a sophisticated electron transfer (ET) machinery to export electrons from the cytosol to extracellular space during extracellular respiration. At the heart of this process are decaheme proteins of the Mtr pathway, MtrC and MtrF, located at the external face of the outer bacterial membrane. Crystal structures have revealed that these proteins bind 10 c-type hemes arranged in the peculiar shape of a staggered cross that trifurcates the electron flow, presumably to reduce extracellular substrates while directing electrons to neighboring multiheme cytochromes at either side along the membrane. Especially intriguing is the design of the heme junctions trifurcating the electron flow: they are made of coplanar and T-shaped heme pair motifs with relatively large and seemingly unfavorable tunneling distances. Here, we use electronic structure calculations and molecular simulations to show that the side chains of the heme rings, in particular the cysteine linkages inserting in the space between coplanar and T-shaped heme pairs, strongly enhance electronic coupling in these two motifs. This results in an [Formula: see text]-fold speedup of ET steps at heme junctions that would otherwise be rate limiting. The predicted maximum electron flux through the solvated proteins is remarkably similar for all possible flow directions, suggesting that MtrC and MtrF shuttle electrons with similar efficiency and reversibly in directions parallel and orthogonal to the outer membrane. No major differences in the ET properties of MtrC and MtrF are found, implying that the different expression levels of the two proteins during extracellular respiration are not related to redox function

Coherent Electron Transport across a 3 nm Bioelectronic Junction Made of Multi-Heme Proteins

Multi-heme cytochromes (MHCs) are fascinating proteins used by bacterial organisms to shuttle electrons within, between, and out of their cells. When placed in solid-state electronic junctions, MHCs support temperature-independent currents over several nanometers that are 3 orders of magnitude higher compared to other redox proteins of similar size. To gain molecular-level insight into their astonishingly high conductivities, we combine experimental photoemission spectroscopy with DFT+Σ current-voltage calculations on a representative Gold-MHC-Gold junction. We find that conduction across the dry, 3 nm long protein occurs via off-resonant coherent tunneling, mediated by a large number of protein valence-band orbitals that are strongly delocalized over heme and protein residues. This picture is profoundly different from the electron hopping mechanism induced electrochemically or photochemically under aqueous conditions. Our results imply that the current output in solid-state junctions can be even further increased in resonance, for example, by applying a gate voltage, thus allowing a quantum jump for next-generation bionanoelectronic devices

NMR spectroscopic detection of chirality and enantiopurity in referenced systems without formation of diastereomers

Enantiomeric excess of chiral compounds is a key parameter that determines their activity or therapeutic action. The current paradigm for rapid measurement of enantiomeric excess using NMR is based on the formation of diastereomeric complexes between the chiral analyte and a chiral resolving agent, leading to (at least) two species with no symmetry relationship. Here we report an effective method of enantiomeric excess determination using a symmetrical achiral molecule as the resolving agent, which is based on the complexation with analyte (in the fast exchange regime) without the formation of diastereomers. The use of N,N′-disubstituted oxoporphyrinogen as the resolving agent makes this novel method extremely versatile, and appropriate for various chiral analytes including carboxylic acids, esters, alcohols and protected amino acids using the same achiral molecule. The model of sensing mechanism exhibits a fundamental linear response between enantiomeric excess and the observed magnitude of induced chemical shift non-equivalence in the 1H NMR spectra