Robust non-adiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of non-adiabatic electronic transitions, and test it against exact integration of the time-dependent Schroedinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the non-adiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.Comment: 36 pages, 5 figures. Accepted for publication in Journal of Chemical Physic

Permittivity measurement of thermoplastic composites at elevated temperature

[Abstract]: The material properties of greatest importance in microwave processing of a dielectric are the complex relative permittivity Epsilon = Epsilon' - jEpsilon'', and the loss tangent, tan Delta = Epsilon'/Epsilon''. This paper describes two convenient laboratory based methods to obtain Epsilon', Epsilon'' and hence tan Delta of fibre-reinforced thermoplastic (FRTP) composites. One method employs a microwave network analyser in conjunction with a waveguide transmission technique, chosen because it provides the widest possible frequency range with high accuracy. The values of the dielectric constant and dielectric loss of glass fibre reinforced (33%) low density polyethylene, LDPE/GF (33%), polystyrene, PS/GF (33%), and Nylon 66/GF (33%), were obtained. Results are compared with those obtained by another method using a high-temperature dielectric probe

Molecular dynamics study of structure and reactions at the hydroxylated Mg(0001)/bulk water interface

A molecular level understanding of the aqueous Mg corrosion mechanism will be essential in developing improved alloys for battery electrodes, automobile parts, and biomedical implants. The structure and reactivity of the hydroxylated surface is expected to be key to the overall mechanism because (i) it is predicted to be the metastable surface state (rather than the bare surface) under a range of conditions and (ii) it provides a reasonable model for the outer corrosion film/water interface. We investigate the structure, interactions, and reactivity at the hydroxylated Mg(0001)/water interface using a combination of static Density Functional Theory calculations and second-generation Car–Parrinello ab initio molecular dynamics. We carry out detailed structural analyses into, among other properties, near-surface water orientations, favored adsorption sites, and near-surface hydrogen bonding behavior. Despite the short timescale (tens of ps) of our molecular dynamics run, we observe a cathodic water splitting event; the rapid timescale for this reaction is explained in terms of near-surface water structuring lowering the reaction barrier. Furthermore, we observe oxidation of an Mg surface atom to effectively generate a univalent Mg species (Mg+). Results are discussed in the context of understanding the Mg corrosion mechanism: For example, our results provide an explanation for the catalytic nature of the Mg corrosion film toward water splitting and a feasible mechanism for the generation of the univalent Mg species often proposed as a key intermediate

Efficient simulations with electronic open boundaries

We present a reformulation of the Hairy Probe method for introducing electronic open boundaries that is appropriate for steady state calculations involving non-orthogonal atomic basis sets. As a check on the correctness of the method we investigate a perfect atomic wire of Cu atoms, and a perfect non-orthogonal chain of H atoms. For both atom chains we find that the conductance has a value of exactly one quantum unit, and that this is rather insensitive to the strength of coupling of the probes to the system, provided values of the coupling are of the same order as the mean inter-level spacing of the system without probes. For the Cu atom chain we find in addition that away from the regions with probes attached, the potential in the wire is uniform, while within them it follows a predicted exponential variation with position. We then apply the method to an initial investigation of the suitability of graphene as a contact material for molecular electronics. We perform calculations on a carbon nanoribbon to determine the correct coupling strength of the probes to the graphene, and obtain a conductance of about two quantum units corresponding to two bands crossing the Fermi surface. We then compute the current through a benzene molecule attached to two graphene contacts and find only a very weak current because of the disruption of the π-conjugation by the covalent bond between the benzene and the graphene. In all cases we find that very strong or weak probe couplings suppress the current

Surface-wave propagation on a grounded dielectric slab covered by a high-permittivity material

A grounded dielectric slab covered by a higher permittivity material would not normally be expected to support surface waves. This conclusion must be modified when the covering material is sufficiently lossy. Assuming a thin slab, approximate analysis shows that the fundamental TM0 surface wave is able to propagate if the cover loss is high enough. A numerical analysis has verified these conclusions. propagation of higher order TM and TE modes is found to be possible above cutoff frequencies, which reduce as cover loss is increased. These results are of significance when printed circuit transmission lines such as microstrip or slotline are used as contact sensors, e.g., for moist materials

Power dissipation in nanoscale conductors: classical, semi-classical and quantum dynamics

Modelling Joule heating is a difficult problem because of the need to introduce correct correlations between the motions of the ions and the electrons. In this paper we analyse three different models of current induced heating (a purely classical model, a fully quantum model and a hybrid model in which the electrons are treated quantum mechanically and the atoms are treated classically). We find that all three models allow for both heating and cooling processes in the presence of a current, and furthermore the purely classical and purely quantum models show remarkable agreement in the limit of high biases. However, the hybrid model in the Ehrenfest approximation tends to suppress heating. Analysis of the equations of motion reveals that this is a consequence of two things: the electrons are being treated as a continuous fluid and the atoms cannot undergo quantum fluctuations. A means for correcting this is suggested

Molecular recognition in olfaction

The mechanism by which the chemical identity of odourants is established by olfactory receptors is a matter of intense debate. Here we present an overview of recent ideas and data with a view to summarising what is known, and what has yet to be determined. We outline the competing theories, and summarise experimental results employing isotopes obtained for mammals, insects, and individual receptors that enable us to judge the relative correctness of the theories

Accelerating GW calculations through machine learned dielectric matrices

The GW approach produces highly accurate quasiparticle energies, but its application to large systems is computationally challenging due to the difficulty in computing the inverse dielectric matrix. To address this challenge, we develop a machine learning approach to efficiently predict density–density response functions (DDRF) in materials. An atomic decomposition of the DDRF is introduced, as well as the neighborhood density–matrix descriptor, both of which transform in the same way under rotations. The resulting DDRFs are then used to evaluate quasiparticle energies via the GW approach. To assess the accuracy of this method, we apply it to hydrogenated silicon clusters and find that it reliably reproduces HOMO–LUMO gaps and quasiparticle energy levels. The accuracy of the predictions deteriorates when the approach is applied to larger clusters than those in the training set. These advances pave the way for GW calculations of complex systems, such as disordered materials, liquids, interfaces, and nanoparticles

Structure and interactions at the Mg(0001)/water interface: An ab initio study

A molecular level understanding of metal/bulk water interface structure is key for a wide range of processes including aqueous corrosion, our focus, but their buried nature makes experimental investigation difficult and means we must mainly rely on simulations. We investigate the Mg(0001)/water interface using second generation Car-Parrinello molecular dynamics (MD) to gain structural information, combined with static density functional theory calculations to probe the atomic interactions and electronic structure (e.g calculating the potential of zero charge). By performing detailed structural analyses of both metal-surface atoms and the near-surface water we find, amongst other insights: i) water adsorption causes significant surface roughening, ii) strongly adsorbed water covers only one quarter of available surface sites and iii) adsorbed water avoids clustering on the surface. Static calculations are used to gain a deeper understanding of the structuring observed in MD. For example, we use an energy decomposition analysis combined with calculated atomic charges to show adsorbate clustering is unfavorable due to Coulombic repulsion between adsorption site surface atoms. Results are discussed in the context of previous simulations of metal/water interfaces. The largest differences for the Mg(0001)/water system appear to be the high degree of surface distortion and minimal difference between the metal work function and metal/water potential of zero charge. The structural information in this paper is important for understanding aqueous Mg corrosion, as the Mg(0001)/water interface is the starting point for key reactions. Furthermore, our focus on understanding the driving forces behind this structuring leads to important insights for general metal/water interfaces