The basis for the formation of stable metal clusters on an electrode surface

Metal nanoclusters can be produced cheaply and precisely in an electrochemical environment. Experimentally this method works in some systems, but not in others, and the unusual stability of the clusters has remained a mystery. We have simulated the deposition of the clusters using classical molecular dynamics and studied their stability by grand-canonical Monte Carlo simulations. We find that electrochemically stable clusters occur only in those cases where the two metals involved form stable alloys.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Dynamical Equations For The Period Vectors In A Periodic System Under Constant External Stress

The purpose of this paper is to derive the dynamical equations for the period vectors of a periodic system under constant external stress. The explicit starting point is Newton’s second law applied to halves of the system. Later statistics over indistinguishable translated states and forces associated with transport of momentum are applied to the resulting dynamical equations. In the final expressions, the period vectors are driven by the imbalance between internal and external stresses. The internal stress is shown to have both full interaction and kinetic-energy terms.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Dynamic resonances and tunnelling in the multiphoton ionization of argon

We present results of wavepacket simulations for multiphoton ionization in argon. A single active electron model is applied to estimate the single-electron ionization rates and photoelectron energy distributions for lambda = 390 nm light with intensities up to I = 2 x 10(14) W cm(-2). The multiphoton ionization rates are compared with R-matrix Floquet calculations and found to be in very good agreement. The photoelectron energy distribution is used to study the nature of ionization at the higher intensities. Our results are consistent with recent calculations and experiments which show the imprint of the tunnelling process in the multiphoton regime. For few-cycle intense pulses, we find that the strong modulation of intensity and increased bandwidth leads to dynamic mixing of the 3d and 5s resonances.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (Pi) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.

Correlated electron-ion dynamics in metallic systems

In this paper we briefly discuss the problem of simulating non-adiabatic processes in systems that are usefully modelled using molecular dynamics. In particular we address the problems associated with metals, and describe two methods that can be applied: the Ehrenfest approximation and correlated electron-ion dynamics (CEID). The Ehrenfest approximation is used to successfully describe the friction force experienced by an energetic particle passing through a crystal, but is unable to describe the heating of a wire by an electric current. CEID restores the proper heating.

A molecular model for drug binding to tandem repeats of telomeric G-quadruplexes.

The extreme 3'-ends of human telomeres consist of 150–250 nucleotides of single-stranded DNA sequence together with associated proteins. Small-molecule ligands can compete with these proteins and induce a conformational change in the DNA to a four-stranded quadruplex arrangement, which is also no longer a substrate for the telomerase enzyme. The modified telomere ends provide signals to the DNA-damage-response system and trigger senescence and apoptosis. Experimental structural data are available on such quadruplex complexes comprising up to four telomeric DNA repeats, but not on longer systems that are more directly relevant to the single-stranded overhang in human cells. The present paper reports on a molecular modelling study that uses Molecular Dynamics simulation methods to build dimer and tetramer quadruplex repeats. These incorporate ligand-binding sites and are models for overhang–ligand complexes.

Comparison of Dynamics of Extracellular Accesses to the β(1) and β(2) Adrenoceptors Binding Sites Uncovers the Potential of Kinetic Basis of Antagonist Selectivity.

From the molecular mechanism of antagonist unbinding in the ß(1) and ß(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of ß(1) -selective Esmolol and ß(2) -selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the ß(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the ß(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the ß(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the ß(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface. © 2012 John Wiley & Sons A/S.

High pressure CO2 absorption studies on imidazolium-based ionic liquids:Experimental and simulation approaches

A combined experimental-computational study on the CO absorption on 1-butyl-3-methylimidazolium hexafluophosphate, 1-ethyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide, and 1-butyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide ionic liquids is reported. The reported results allowed to infer a detailed nanoscopic vision of the absorption phenomena as a function of pressure and temperature. Absorption isotherms were measured at 318 and 338K for pressures up to 20MPa for ultrapure samples using a state-of-the-art magnetic suspension densimeter, for which measurement procedures are developed. A remarkable swelling effect upon CO absorption was observed for pressures higher than 10MPa, which was corrected using a method based on experimental volumetric data. The experimental data reported in this work are in good agreement with available literature isotherms. Soave-Redlich-Kwong and Peng-Robinson equations of state coupled with bi-parametric van der Waals mixing rule were used for successful correlations of experimental high pressure absorption data. Molecular dynamics results allowed to infer structural, energetic and dynamic properties of the studied CO+ionic liquids mixed fluids, showing the relevant role of the strength of anion-cation interactions on fluid volumetric properties and CO absorption. © 2012 Elsevier B.V.

Comparison of dynamics of wildtype and V94M human UDP-galactose 4-epimerase-A computational perspective on severe epimerase-deficiency galactosemia

UDP-galactose 4'-epimerase (GALE) catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. Type III galactosemia, an inherited metabolic disease, is associated with mutations in human GALE. The V94M mutation has been associated with a very severe form of type III galactosemia. While a variety of structural and biochemical studies have been reported that elucidate differences between the wildtype and this mutant form of human GALE, little is known about the dynamics of the protein and how mutations influence structure and function. We performed molecular dynamics simulations on the wildtype and V94M enzyme in different states of substrate and cofactor binding. In the mutant, the average distance between the substrate and both a key catalytic residue (Tyr157) and the enzyme-bound NAD(+) cofactor and the active site dynamics are altered making substrate binding slightly less stable. However, overall stability or dynamics of the protein is not altered. This is consistent with experimental findings that the impact is largely on the turnover number (kcat), with less substantial effects on Km. Active site fluctuations were found to be correlated in enzyme with substrate bound to just one of the subunits in the homodimer suggesting inter-subunit communication. Greater active site loop mobility in human GALE compared to the equivalent loop in Escherichia coli GALE explains why the former can catalyze the interconversion of UDP-N-acetylgalactosamine and UDP-N-acetylglucosamine while the bacterial enzyme cannot. This work illuminates molecular mechanisms of disease and may inform the design of small molecule therapies for type III galactosemia.

Viscous behavior of imidazolium-based ionic liquids

The viscosity of four imidazolium-based ionic liquids is analyzed as a function of pressure and temperature. Experimental measurements were carried out using an electromagnetic moving piston viscometer in the 303-353 K and 0.1-70 MPa ranges on synthesized ultrapure samples, and compared with available literature data. Molecular dynamics simulations were used to analyze the fluids' dynamic properties from a nanoscopic viewpoint, with special attention paid to self-diffusion coefficients and dynamic viscosity. Simulated properties are in excellent agreement with experimental results in spite of the glasslike dynamics of some of the studied fluids. © 2013 American Chemical Society.

Topical Review: The current understanding on the diamond machining of silicon carbide

The Glenn Research Centre of NASA, USA (www.grc.nasa.gov/WWW/SiC/, silicon carbide electronics) is in pursuit of realizing bulk manufacturing of silicon carbide (SiC), specifically by mechanical means. Single point diamond turning (SPDT) technology which employs diamond (the hardest naturally-occurring material realized to date) as a cutting tool to cut a workpiece is a highly productive manufacturing process. However, machining SiC using SPDT is a complex process and, while several experimental and analytical studies presented to date aid in the understanding of several critical processes of machining SiC, the current knowledge on the ductile behaviour of SiC is still sparse. This is due to a number of simultaneously occurring physical phenomena that may take place on multiple length and time scales. For example, nucleation of dislocation can take place at small inclusions that are of a few atoms in size and once nucleated, the interaction of these nucleations can manifest stresses on the micrometre length scales. The understanding of how stresses manifest during fracture in the brittle range, or dislocations/phase transformations in the ductile range, is crucial in understanding the brittle–ductile transition in SiC. Furthermore, there is a need to incorporate an appropriate simulation-based approach in the manufacturing research on SiC, owing primarily to the number of uncertainties in the experimental research that includes wear of the cutting tool, poor controllability of the nano-regime machining scale (effective thickness of cut), and coolant effects (interfacial phenomena between the tool, workpiece/chip and coolant), etc. In this review, these two problems are combined together to posit an improved understanding on the current theoretical knowledge on the SPDT of SiC obtained from molecular dynamics simulation.

On transition rates in surface hopping

Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work, we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely (i) non-zero electronic transition rate ("hopping probability") at avoided crossings; (ii) rescaling of the nuclear velocities to conserve total energy; (iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4770280]

Ultracold, radiative charge transfer in hybrid Yb ion – Rb atom traps

Ultracold hybrid ion–atom traps offer the possibility of microscopic manipulation of quantum coherences in the gas using the ion as a probe. However, inelastic processes, particularly charge transfer can be a significant process of ion loss and has been measured experimentally for the ${\rm Y}{{{\rm b}}^{+}}$ ion immersed in a Rb vapour. We use first-principles quantum chemistry codes to obtain the potential energy curves and dipole moments for the lowest-lying energy states of this complex. Calculations for the radiative decay processes cross sections and rate coefficients are presented for the total decay processes; ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm Yb}(6{{{\rm s}}^{2}}{{\;}^{1}}{\rm S})+{\rm R}{{{\rm b}}^{+}}(4{{{\rm p}}^{6}}{{\;}^{1}}{\rm S})+h\nu $ and ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm YbR}{{{\rm b}}^{+}}({{X}^{1}}{{\Sigma }^{+}})+h\nu $. Comparing the semi-classical Langevin approximation with the quantum approach, we find it provides a very good estimate of the background at higher energies. The results demonstrate that radiative decay mechanisms are important over the energy and temperature region considered. In fact, the Langevin process of ion–atom collisions dominates cold ion–atom collisions. For spin-dependent processes [1] the anisotropic magnetic dipole–dipole interaction and the second-order spin–orbit coupling can play important roles, inducing coupling between the spin and the orbital motion. They measured the spin-relaxing collision rate to be approximately five orders of magnitude higher than the charge-exchange collision rate [1]. Regarding the measured radiative charge transfer collision rate, we find that our calculation is in very good agreement with experiment and with previous calculations. Nonetheless, we find no broad resonances features that might underly a strong isotope effect. In conclusion, we find, in agreement with previous theory that the isotope anomaly observed in experiment remains an open question.