Source author record

Paul A. Madden

Paul A. Madden appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2015arXiv

Polarization effects in ionic solids and melts

Ionic solids and melts are compounds in which the interactions are dominated by electrostatic effects. However, the polarization of the ions also plays an important role in many respects as has been clarified in recent years thanks to the development of realistic polarizable interaction potentials. After detailing these models, we illustrate the importance of polarization effects on a series of examples concerning the structural properties, such as the stabilization of particular crystal structures or the formation of highly-coordinated multivalent ions in the melts, as well as the dynamic properties such as the diffusion of ionic species. The effects on the structure of molten salts interfaces (with vacuum and electrified metal) is also described. Although most of the results described here concern inorganic compounds (molten fluorides and chlorides, ionic oxides...), the particular case of the room-temperature ionic liquids, a special class of molten salts in which at least one species is organic, will also be briefly discussed to indicate how the ideas gained from the study of "simple" molten salts are being transferred to these more complex systems.

preprint2015arXiv

Water exchange at a hydrated platinum electrode is rare and collective

We use molecular dynamics simulations to study the exchange kinetics of water molecules at a model metal electrode surface -- exchange between water molecules in the bulk liquid and water molecules bound to the metal. This process is a rare event, with a mean residence time of a bound water of about 40 ns for the model we consider. With analysis borrowed from the techniques of rare-event sampling, we show how this exchange or desorption is controlled by (1) reorganization of the hydrogen bond network within the adlayer of bound water molecules, and by (2) interfacial density fluctuations of the bulk liquid adjacent to the adlayer. We define collective coordinates that describe the desorption mechanism. Spatial and temporal correlations associated with a single event extend over nanometers and tens of picoseconds.

preprint2014arXiv

Effect of dispersion interactions on the properties of LiF in condensed phases

Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap-repulsion and polarization terms of the interaction potential are derived on a purely non empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parameterizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained by direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the density of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.

preprint2014arXiv

Molecular Dynamics Simulation of Coherent Interfaces in Fluorite Heterostructures

The standard model of enhanced ionic conductivities in solid electrolyte heterostructures follows from a continuum mean-field description of defect distributions that makes no reference to crystalline structure. To examine ionic transport and defect distributions while explicitly accounting for ion-ion correlations and lattice effects, we have performed molecular dynamics simulations of a model coherent fluorite heterostructure without any extrinsic defects, with a difference in standard chemical potentials of mobile fluoride ions between phases induced by an external potential. Increasing the offset in fluoride ion standard chemical potentials across the internal interfaces decreases the activation energies for ionic conductivity and diffusion and strongly enhances fluoride ion mobilities and defect concentrations near the heterostructure interfaces. Non-charge-neutral "space-charge" regions, however, extend only a few atomic spacings from the interface, suggesting a continuum model may be inappropriate. Defect distributions are qualitatively inconsistent with the predictions of the continuum mean-field model, and indicate strong lattice-mediated defect-defect interactions. We identify an atomic-scale "Frenkel polarisation" mechanism for the interfacial enhancement in ionic mobility, where preferentially oriented associated Frenkel pairs form at the interface and promote local ion mobility via concerted diffusion processes.

preprint2014arXiv

Relationships Between Atomic Diffusion Mechanisms and Ensemble Transport Coefficients in Crystalline Polymorphs

Ionic transport in conventional ionic solids is generally considered to proceed via independent diffusion events or "hops''. This assumption leads to well-known Arrhenius expressions for transport coefficients, and is equivalent to assuming diffusion is a Poisson process. Using molecular dynamics simulations of the low-temperature B1, B3, and B4 AgI polymorphs, we have compared rates of ion-hopping with corresponding Poisson distributions to test the assumption of independent hopping in these common structure-types. In all cases diffusion is a non-Poisson process, and hopping is strongly correlated in time. In B1 the diffusion coefficient can be approximated by an Arrhenius expression, though the physical significance of the parameters differs from that commonly assumed. In low temperature B3 and B4 diffusion is characterised by concerted motion of multiple ions in short closed loops. Diffusion coefficients can not be expressed in a simple Arrhenius form dependent on single-ion free-energies, and intrinsic diffusion must be considered a many-body process.

preprint2014arXiv

Thermal Conductivity of Simple Liquids:Origin of Temperature and Packing Fraction Dependences

The origin of both weak temperature dependence and packing fraction dependence of $T^{1/4}η^{3/2}$ in the thermal conductivity of the simple Lennard-Jones (LJ) liquid is explored. In order to discuss the relative contributions from attractive or repulsive part of the interaction potential separately, the thermal conductivity of a series of Weeks-Chandler-Anderson (WCA) fluids is calculated by molecular dynamics simulations. The results show that the repulsive part plays the main role in the heat conduction, while the attractive part has no direct effect on the thermal conductivity for a given packing fraction. By investigating WCA fluids with potentials of varying softness, we explain the difference observed between the LJ liquids such as argon and Coulombic liquids such as NaCl.

preprint2013arXiv

Characterizing heterogeneous dynamics at hydrated electrode surfaces

In models of Pt 111 and Pt 100 surfaces in water, motions of molecules in the first hydration layer are spatially and temporally correlated. To interpret these collective motions, we apply quantitative measures of dynamic heterogeneity that are standard tools for considering glassy systems. Specifically, we carry out an analysis in terms of mobility fields and distributions of persistence times and exchange times. In so doing, we show that dynamics in these systems is facilitated by transient disorder in frustrated two-dimensional hydrogen bonding networks. The frustration is the result of unfavorable geometry imposed by strong metal-water bonding. The geometry depends upon the structure of the underlying metal surface. Dynamic heterogeneity of water on the Pt 111 surface is therefore qualitatively different than that for water on the Pt 100 surface. In both cases, statistics of this adlayer dynamic heterogeneity responds asymmetrically to applied voltage.

preprint2013arXiv

Charge fluctuations in nano-scale capacitors

The fluctuations of the charge on an electrode contain information on the microscopic correlations within the adjacent fluid and their effect on the electronic properties of the interface. We investigate these fluctuations using molecular dynamics simulations in a constant-potential ensemble with histogram reweighting techniques. This approach offers in particular an efficient, accurate and physically insightful route to the differential capacitance that is broadly applicable. We demonstrate these methods with three different capacitors: pure water between platinum electrodes, and a pure as well as a solvent-based organic electrolyte each between graphite electrodes. The total charge distributions with the pure solvent and solvent-based electrolytes are remarkably Gaussian, while in the pure ionic liquid the total charge distribution displays distinct non-Gaussian features, suggesting significant potential-driven changes in the organization of the interfacial fluid.

preprint2013arXiv

From molten salts to room temperature ionic liquids: Simulation studies on chloroaluminate systems

An interaction potential including chloride anion polarization effects, constructed from first-principles calculations, is used to examine the structure and transport properties of a series of chloroaluminate melts. A particular emphasis was given to the study of the equimolar mixture of aluminium chloride with 1-ethyl-3-methylimidazolium chloride, which forms a room temperature ionic liquid EMI-AlCl 4. The structure yielded by the classical simulations performed within the framework of the polarizable ion model is compared to the results obtained from entirely electronic structure-based simulations: An excellent agreement between the two flavors of molecular dynamics is observed. When changing the organic cation EMI+ by an inorganic cation with a smaller ionic radius (Li+, Na+, K+), the chloroaluminate speciation becomes more complex, with the formation of Al2Cl 7- in small amounts. The calculated transport properties (diffusion coefficients, electrical conductivity and viscosity) of EMI-AlCl4 are in good agreement with experimental data.

preprint2013arXiv

Structure and dynamics in yttrium-based molten rare earth alkali fluorides

The transport properties of molten LiF-YF3 mixtures have been studied by pulsed field gradient nuclear magnetic resonance spectroscopy, potentiometric experiments, and molecular dynamics simulations. The calculated diffusion coefficients and electric conductivities compare very well with the measurements accross a wide composition range. We then extract static (radial distribution functions, coordination numbers distributions) and dynamic (cage correlation functions) quantities from the simulations. Then, we discuss the interplay between the microscopic structure of the molten salts and their dynamic properties. It is often considered that variations in the diffusion coefficient of the anions are mainly driven by the evolution of its coordination with the metallic ion (Y3+ here). We compare this system with fluorozirconate melts and demonstrate that the coordination number is a poor indicator of the evolution of the diffusion coefficient. Instead, we propose to use the ionic bonds lifetime. We show that the weak Y-F ionic bonds in LiF-YF3 do not induce the expected tendency of the fluoride diffusion coefficient to converge toward the one of yttrium cation when the content in YF3 increases. Implications on the validity of the Nernst-Einstein relation for estimating the electrical conductivity are discussed.

preprint2012arXiv

Including many-body effects in models for ionic liquids

Realistic modeling of ionic systems necessitates taking explicitly account of many-body effects. In molecular dynamics simulations, it is possible to introduce explicitly these effects through the use of additional degrees of freedom. Here we present two models: The first one only includes dipole polarization effect, while the second also accounts for quadrupole polarization as well as the effects of compression and deformation of an ion by its immediate coordination environment. All the parameters involved in these models are extracted from first-principles density functional theory calculations. This step is routinely done through an extended force-matching procedure, which has proven to be very succesfull for molten oxides and molten fluorides. Recent developments based on the use of localized orbitals can be used to complement the force-matching procedure by allowing for the direct calculations of several parameters such as the individual polarizabilities.

preprint2012arXiv

Solvation at Aqueous Metal Electrodes

We present a study of the solvation properties of model aqueous electrode interfaces. The exposed electrodes we study strongly bind water and have closed packed crystalline surfaces, which template an ordered water adlayer adjacent to the interface. We find that these ordered water structures facilitate collective responses in the presence of solutes that are correlated over large lengthscales and across long timescales. Specifically, we show that the liquid water adjacent to the ordered adlayers forms a soft, liquid-vapor-like interface with concomitant manifestations of hydrophobicity. Temporal defects in the adlayer configurations create a dynamic heterogeneity in the degree to which different regions of the interface attract hydrophobic species. The structure and heterogeneous dynamics of the adlayer defects depend upon the geometry of the underlying ordered metal surface. For both 100 and 111 surfaces, the dynamical heterogeneity relaxes on times longer than nanoseconds. Along with analyzing time scales associated with these effects, we highlight implications for electrolysis and the particular catalytic efficiency of platinum.

preprint2008arXiv

Water at an electrochemical interface - a simulation study

The results of molecular dynamics simulations of the properties of water in an aqueous ionic solution close to an interface with a model metallic electrode are described. In the simulations the electrode behaves as an ideally polarizable hydrophilic metal, supporting image charge interactions with charged species, and it is maintained at a constant electrical potential with respect to the solution so that the model is a textbook representation of an electrochemical interface through which no current is passing. We show how water is strongly attracted to and ordered at the electrode surface. This ordering is different to the structure that might be imagined from continuum models of electrode interfaces. Further, this ordering significantly affects the probability of ions reaching the surface. We describe the concomitant motion and configurations of the water and ions as functions of the electrode potential, and we analyze the length scales over which ionic atmospheres fluctuate. The statistics of these fluctuations depend upon surface structure and ionic strength. The fluctuations are large, sufficiently so that the mean ionic atmosphere is a poor descriptor of the aqueous environment near a metal surface. The importance of this finding for a description of electrochemical reactions is examined by calculating, directly from the simulation, Marcus free energy profiles for transfer of charge between the electrode and a redox species in the solution and comparing the results with the predictions of continuum theories. Significant departures from the electrochemical textbook descriptions of the phenomenon are found and their physical origins are characterized from the atomistic perspective of the simulations.