Voltage equilibration for reactive atomistic simulations of electrochemical processes

Significant Power Consumption Reduction and Speed Boosting in Phase Change Memory with Nanocurrent Channels

The excessive power consumption is challenging for phase change memory (PCM) on its way to becoming universal memory in complex hierarchies of memory systems. Here, from the perspective of device structure, by adding a nanocurrent-channel (NCC) layer between the electrode layer and phase change layer, a RESET power consumption reduction by more than 95% and 10 times faster SET speed were realized simultaneously. Through the first principle calculations, Au and SiO2 were screened as the metal and insulating matrix material of NCC layer, respectively. Our PCM device with a Au-SiO2 NCC layer shows an ultralow RESET power consumption, down to 381 fJ, and an ultrafast SET speed (8 ns). Much higher current density near NCC in the phase change layer and thermal barrier effect of insulating matrix material were confirmed by finite element analysis (FEA), and the role of Au nanochannels was revealed by transmission electron microscopy (TEM). Our NCC layer structure provides a simple and practicable method to significantly decrease PCM power consumption.

Nanoionics enabled atomic point contact construction and quantum conductance effects

The miniaturization of electronic devices is important for the development of high-density and function-integrated information devices. Atomic-point-contact (APC) structures refer to narrow contact areas formed by one or more atoms between two conductive electrodes that produce quantum conductance effects when the electrons pass through the APC channel, providing a new development path for the miniaturization of information devices. Recently, nanoionics has enabled the electric field reconfiguration of APC structures in solid-state electrolytes, offering new approaches to controlling the quantum conductance states, which may lead to the development of emerging information technologies with low power consumption, high speed, and high density. This review provides an overview of APC structures with a focus on the fabrication methods enabled by nanoionics technology. In particular, the advantages of electric field-driven nanoionics in the construction of APC structures are summarized, and the influence of external fields on quantum conductance effects is discussed. Recent studies on electric field regulation of APC structures to achieve precise control of quantum conductance states are also reviewed. The potential applications of quantum conductance effects in memory, computing, and encryption-related information technologies are further explored. Finally, the challenges and future prospects of quantum conductance effects in APC structures are discussed.

Effect of Electric Fields on the Decomposition of Phosphate Esters

Phosphate esters decompose on metal surfaces and form protective polyphosphate films. For many applications, such as in lubricants for electric vehicles and wind turbines, an understanding of the effect of electric fields on molecular decomposition is urgently required. Experimental investigations have yielded contradictory results, with some suggesting that electric fields improve tribological performance, while others have reported the opposite effect. Here, we use nonequilibrium molecular dynamics (NEMD) simulations to study the decomposition of tri-n-butyl phosphate (TNBP) molecules nanoconfined between ferrous surfaces (iron and iron oxide) under electrostatic fields. The reactive force field (ReaxFF) method is used to model the effects of chemical bonding and molecular dissociation. We show that the charge transfer with the polarization current equalization (QTPIE) method gives more realistic behavior compared to the standard charge equilibration (QEq) method under applied electrostatic fields. The rate of TNBP decomposition via carbon-oxygen bond dissociation is faster in the nanoconfined systems than that in the bulk due to the catalytic action of the surfaces. In all cases, the application of an electric field accelerates TNBP decomposition. When electric fields are applied to the confined systems, the phosphate anions are pulled toward the surface with high electric potential, while the alkyl cations are pulled to the surface with lower potential, leading to asymmetric film growth. Analysis of the temperature- and electric field strength-dependent dissociation rate constants using the Arrhenius equation suggests that, on reactive iron surfaces, the increased reactivity under an applied electric field is driven mostly by an increase in the pre-exponential factor, which is linked to the number of molecule-surface collisions. Conversely, the accelerated decomposition of TNBP on iron oxide surfaces can be attributed to a reduction in the activation energy with increasing electric field strength. Single-molecule nudged-elastic band (NEB) calculations also show a linear reduction in the energy barrier for carbon-oxygen bond breaking with electric field strength, due to stabilization of the charged transition state. The simulation results are consistent with experimental observations of enhanced and asymmetric tribofilm growth under electrostatic fields.

Metal Nanogap Memory: Performances and Switching Mechanism

The nanogap memory (NGM) device, emerging as a promising nonvolatile memory candidate, has attracted increasing attention for its simple structure, nano/atomic scale size, elevated operating speed, and robustness to high temperatures. In this study, nanogap memories based on Pd, Au, and Pt were fabricated by combining nanofabrication with electromigration technology. Subsequent evaluations of the electrical characteristics were conducted under ambient air or vacuum conditions at room temperature. The investigation unveiled persistent challenges associated with metal NGM devices, including (1) prolonged SET operation time in comparison to RESET, (2) the potential generation of error bits when enhancing switching speeds, and (3) susceptibility to degradation during program/erase cycles. While these issues have been encountered by predecessors in NGM device development, the underlying causes have remained elusive. Employing molecular dynamics (MD) simulation, we have, for the first time, unveiled the dynamic processes of NGM devices during both SET and RESET operations. The MD simulation highlights that the adjustment of the tunneling gap spacing in nanogap memory primarily occurs through atomic migration or field evaporation. This dynamic process enables the device to transition between the high-resistance state (HRS) and the low-resistance state (LRS). The identified mechanism provides insight into the origins of the aforementioned challenges. Furthermore, the study proposes an effective method to enhance the endurance of NGM devices based on the elucidated mechanism.

Understanding electrochemical interfaces through comparing experimental and computational charge density-potential curves

Electrode-electrolyte interfaces play a decisive role in electrochemical charge accumulation and transfer processes. Theoretical modelling of these interfaces is critical to decipher the microscopic details of such phenomena. Different force field-based molecular dynamics protocols are compared here in a view to connect calculated and experimental charge density-potential relationships. Platinum-aqueous electrolyte interfaces are taken as a model. The potential of using experimental charge density-potential curves to transform cell voltage into electrode potential in force-field molecular dynamics simulations, and the need for that purpose of developing simulation protocols that can accurately calculate the double-layer capacitance, are discussed.

Electrochemical rewiring through quantum conductance effects in single metallic memristive nanowires

Memristive devices have been demonstrated to exhibit quantum conductance effects at room temperature. In these devices, a detailed understanding of the relationship between electrochemical processes and ionic dynamic underlying the formation of atomic-sized conductive filaments and corresponding electronic transport properties in the quantum regime still represents a challenge. In this work, we report on quantum conductance effects in single memristive Ag nanowires (NWs) through a combined experimental and simulation approach that combines advanced classical molecular dynamics (MD) algorithms and quantum transport simulations (DFT). This approach provides new insights on quantum conductance effects in memristive devices by unravelling the intrinsic relationship between electronic transport and atomic dynamic reconfiguration of the nanofilment, by shedding light on deviations from integer multiples of the fundamental quantum of conductance depending on peculiar dynamic trajectories of nanofilament reconfiguration and on conductance fluctuations relying on atomic rearrangement due to thermal fluctuations.

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

Co-Ion Desorption as the Main Charging Mechanism in Metallic 1T-MoS2 Supercapacitors

Metallic 1T-MoS₂ is a promising electrode material for supercapacitor applications. Its layered structure allows the efficient intercalation of ions, leading to experimental volumetric capacitance as high as 140 F/cm³. Molecular dynamics could in principle be used to characterize its charging mechanism; however, unlike conventional nanoporous carbon, 1T-MoS₂ is a multicomponent electrode. The Mo and S atoms have very different electronegativities so that 1T-MoS₂ cannot be simulated accurately using the conventional constant potential method. In this work, we show that controlling the electrochemical potential of the atoms allows one to recover average partial charges for the elements in agreement with electronic structure calculations for the material at rest, without compromising the ability to simulate systems under an applied voltage. The simulations yield volumetric capacitances in agreement with experiments. We show that due to the large electronegativity of S, the co-ion desorption is the main charging mechanism at play during the charging process. This contrasts drastically with carbon materials for which ion exchange and counterion adsorption usually dominate. In the future, our method can be extended to the study of a wide range of families of 2D layered materials such as MXenes.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

A molecular perspective on induced charges on a metallic surface

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behavior of a gold electrode interacting with sodium or chloride ions fixed at various positions, in a vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially under vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely, neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents, and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

Atomistic Insights on the Full Operation Cycle of a HfO2-Based Resistive Random Access Memory Cell from Molecular Dynamics

We characterize the atomic processes that underlie forming, reset, and set in HfO₂-based resistive random access memory (RRAM) cells through molecular dynamics (MD) simulations, using an extended charge equilibration method to describe external electric fields. By tracking the migration of oxygen ions and the change in coordination of Hf atoms in the dielectric, we characterize the formation and dissolution of conductive filaments (CFs) during the operation of the device with atomic detail. Simulations of the forming process show that the CFs form through an oxygen exchange mechanism, induced by a cascade of oxygen displacements from the oxide to the active electrode, as opposed to aggregation of pre-existing oxygen vacancies. However, the filament breakup is dominated by lateral, rather than vertical (along the filament), motion of vacancies. In addition, depending on the temperature of the system, the reset can be achieved through a redox effect (bipolar switch), where oxygen diffusion is governed by the applied bias, or by a thermochemical process (unipolar switch), where the diffusion is driven by temperature. Unlike forming and similar to reset, the set process involves lateral oxygen atoms as well. This is driven by field localization associated with conductive paths.

Effect of the metallicity on the capacitance of gold-aqueous sodium chloride interfaces

Electrochemistry experiments have established that the capacitance of electrode-electrolyte interfaces is much larger for good metals, such as gold and platinum, than for carbon-based materials. Despite the development of elaborate electrode interaction potentials, to date molecular dynamics simulations are not able to capture this effect. Here, we show that changing the width of the Gaussian charge distribution used to represent the atomic charges in gold is an effective way to tune its metallicity. Larger Gaussian widths lead to a capacitance of aqueous solutions (pure water and 1 M NaCl) in good agreement with recent ab initio molecular dynamics results. For pure water, the increase in the capacitance is not accompanied by structural changes, while in the presence of salt, the Na⁺ cations tend to adsorb significantly on the surface. For a strongly metallic gold electrode, these ions can even form inner sphere complexes on hollow sites of the surface.

Hybrid Polymer-Garnet Materials for All-Solid-State Energy Storage Devices

Hybrid electrolyte materials comprising polymer-ionic salt matrixes embedded with garnet particles constitute a promising class of materials for the realization of all-solid-state batteries. In addition to providing solutions to the safety issues inherent to current liquid electrolytes, hybrid polymer electrolytes offer advantages over other solid-state electrolytes. This is because their functional properties such as ionic conductivity, electrochemical stability, and mechanical and thermal properties can be tailored to a particular application by independently optimizing the properties of the constituent materials. This independent optimization permits the rational design of solid-state electrolytes, thereby solving the current bottlenecks that prevent their practical implementation into battery devices. This Mini-Review starts with a survey of solid-state electrolytes, focusing on their materials and ion transport limitations. Next, we summarize the current understanding of transport mechanisms in composite polymer electrolytes (CPEs) with the purpose of identifying materials' solutions for further improving their properties. The overall goal of the Mini-Review is to foster heightened research interest in these hybrid structures to rapidly advance development of future all-solid-state battery devices.

Molecular Simulation of Electrode-Solution Interfaces

Many key industrial processes, from electricity production, conversion, and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of an electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role in the structure, thermodynamics, dynamics, and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling because they are a place where quantum chemistry meets statistical physics. In the present review, we explore the recent advances in the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt electrolytes.

Controlling potential difference between electrodes based on self-consistent-charge density functional tight binding

A proper understanding and description of the electronic response of the electrode surfaces in electrochemical systems are quite important because the interactions between the electrode surface and electrolyte give rise to unique and useful interfacial properties. Atomistic modeling of the electrodes requires not only an accurate description of the electronic response under a constant-potential condition but also computational efficiency in order to deal with systems large enough to investigate the interfacial electrolyte structures. We thus develop a self-consistent-charge density functional tight binding based method to model a pair of electrodes in electrochemical cells under the constant-potential condition. The method is more efficient than the (ab initio) density functional theory calculations so that it can treat systems as large as those studied in classical atomistic simulations. It can also describe the electronic response of electrodes quantum mechanically and more accurately than the classical counterparts. The constant-potential condition is introduced through a Legendre transformation of the electronic energy with respect to the difference in the number of electrons in the two electrodes and their electrochemical potential difference, through which the Kohn-Sham equations for each electrode are variationally derived. The method is applied to platinum electrodes faced parallel to each other under an applied voltage. The electronic response to the voltage and a charged particle is compared with the result of a classical constant-potential method based on the chemical potential equalization principle.

Molecular Simulations of the Microstructure Evolution of Solid Electrolyte Interphase during Cyclic Charging/Discharging

Lithium (Li) metal is regarded as one of the most promising anode materials for use in next-generation high-energy-density rechargeable batteries because of its high volumetric and gravimetric specific capacity, as well as low reduction potential. Unfortunately, uncontrolled dendritic Li growth during cyclic charging/discharging leads to low columbic efficiency and critical safety issues. Hence, comprehensive understanding of the formation mechanism for Li-dendrite growth, particularly at the onset of dendrite formation, is essential for developing Li-metal anode batteries. In this study, reactive molecular dynamics (MD) simulations in combination with the electrochemical dynamics with implicit degrees of freedom (EChemDID) method were performed to investigate the formation and evolution of solid electrolyte interphase (SEI) films for a Li-metal anode under cyclic charging/discharging processes in two distinct dimensions, namely, electrolyte compositions and initial surface morphologies. Our simulations indicated that regardless of the electrolyte compositions and initial anode morphologies, inhomogeneous Li reduction, namely, the formation of Li-reduction "hotspots" during cyclic charging cycles, took place and could serve as the seed for subsequent dendrite growth. The fluorine-containing electrolyte additives could notably mitigate the Li-anode roughening processes by forming dense-SEI-layer products or suppressing electrolyte decomposition. A series of Li-ion-drifting simulations suggest that Li ions navigate through the SEI layer via pathways composed of low-density atoms and become reduced at these reduction hotspots, promoting inhomogeneous deposition and subsequent dendrite growth. The present study reveals atomistic details of the early stage of dendrite growth during cyclic loadings under different electrolyte compositions and anode morphologies, thereby providing insights for designing artificial SEI layers or electrolytes for long-life, high-capacity Li-ion batteries.

Glassy Li metal anode for high-performance rechargeable Li batteries

Lithium metal has been considered an ideal anode for high-energy rechargeable Li batteries, although its nucleation and growth process remains mysterious, especially at the nanoscale. Here, cryogenic transmission electron microscopy was used to reveal the evolving nanostructure of Li metal deposits at various transient states in the nucleation and growth process, in which a disorder-order phase transition was observed as a function of current density and deposition time. The atomic interaction over wide spatial and temporal scales was depicted by reactive molecular dynamics simulations to assist in understanding the kinetics. Compared to crystalline Li, glassy Li outperforms in electrochemical reversibility, and it has a desired structure for high-energy rechargeable Li batteries. Our findings correlate the crystallinity of the nuclei with the subsequent growth of the nanostructure and morphology, and provide strategies to control and shape the mesostructure of Li metal to achieve high performance in rechargeable Li batteries.

Insights into dynamic sliding contacts from conductive atomic force microscopy

Friction in nanoscale contacts is determined by the size and structure of the interface that is hidden between the contacting bodies. One approach to investigating the origins of friction is to measure electrical conductivity as a proxy for contact size and structure. However, the relationships between contact, friction and conductivity are not fully understood, limiting the usefulness of such measurements for interpreting dynamic sliding properties. Here, atomic force microscopy (AFM) was used to simultaneously acquire lattice resolution images of the lateral force and current flow through the tip-sample contact formed between a highly oriented pyrolytic graphite (HOPG) sample and a conductive diamond AFM probe to explore the underlying mechanisms and correlations between friction and conductivity. Both current and lateral force exhibited fluctuations corresponding to the periodicity of the HOPG lattice. Unexpectedly, while lateral force increased during stick events of atomic stick-slip, the current decreased exponentially. Molecular dynamics (MD) simulations of a simple model system reproduced these trends and showed that the origin of the inverse correlation between current and lateral force during atomic stick-slip was atom-atom distance across the contact. The simulations further demonstrated transitions between crystallographic orientation during slip events were reflected in both lateral force and current. These results confirm that the correlation between conduction and atom-atom distance previously proposed for stationary contacts can be extended to sliding contacts in the stick-slip regime.

Charge fluctuations from molecular simulations in the constant-potential ensemble

Statistical mechanics of constant-potential molecular simulations yields a new fluctuation–dissipation relation for the differential capacitance.

Simulations of the effect of an oxide on contact area measurements from conductive atomic force microscopy

Nanoscale contact area in conductive atomic force microscopy can be determined by analyzing current flow using electron transport theories. However, it is recognized that native oxides on the conductive tip will reduce current flow, thus degrading the accuracy of the measured contact area. To quantify the adverse effect of an oxide on contact area measurements, we use molecular dynamics simulations of an oxide-coated platinum tip and a crystalline platinum substrate, where both the contact size and conductance can be inferred from the positions of atoms in the interface. We develop a method to approximate conductance based on the distance between atoms in platinum channels across the contact. Then, the contact area calculated from conductance using ballistic transport and tunneling theories is compared to that obtained using the known positions of atoms in the contact. The difference is small for very thin (<0.1 nm) or very thick (>1.0 nm) oxides, where ballistic transport and tunneling theories work well; however, the difference is significant for oxides between these limits, which is expected to be the case for platinum in many practical applications.

Pulse Dynamics of Electric Double Layer Formation on All-Solid-State Graphene Field-Effect Transistors

Electric double layer (EDL) dynamics in graphene field-effect transistors (FETs) gated with polyethylene oxide (PEO)-based electrolytes are studied by molecular dynamics (MD) simulations from picoseconds to nanoseconds and experimentally from microseconds to milliseconds. Under an applied field of approximately mV/nm, EDL formation on graphene FETs gated with PEO:CsClO₄ occurs on the timescale of microseconds at room temperature and strengthens within 1 ms to a sheet carrier density of n_S ≈ 10¹³ cm^-2. Stronger EDLs (i.e., larger n_S) are induced experimentally by pulsing with applied voltages exceeding the electrochemical window of the electrolyte; electrochemistry is avoided using short pulses of a few milliseconds. Dynamics on picosecond to nanosecond timescales are accessed using MD simulations of PEO:LiClO₄ between graphene electrodes with field strengths of hundreds of mV/nm which is 100× larger than experiment. At 100 mV/nm, EDL formation initiates in sub-nanoseconds achieving charge densities up to 6 × 10¹³ cm^-2 within 3 nanoseconds. The modeling shows that under sufficiently high electric fields, EDLs with densities ∼10¹³ cm^-2 can form within a nanosecond, which is a timescale relevant for high-performance electronics such as EDL transistors (EDLTs). Moreover, the combination of experiment and modeling shows that the timescale for EDL formation ( n_S = 10¹³ to 10¹⁴ cm^-2) can be tuned by 9 orders of magnitude by adjusting the field strength by only 3 orders of magnitude.

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard-Jones potential

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard–Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.

Molecular dynamics models for predicting the behavior of metallic nanostructures typically do not take into account polarization effects in metals. Here, the authors introduce a polarizable Lennard–Jones potential that provides quantitative insight into the role of induced charges at metal surfaces and related complex material interfaces.

Atomistic simulations of contact area and conductance at nanoscale interfaces

Atomistic simulations were used to study conductance across the interface between a nanoscale gold probe and a graphite surface with a step edge. Conductance on the graphite terrace was observed to increase with load and be approximately proportional to contact area calculated from the positions of atoms in the interface. The relationship between area and conductance was further explored by varying the position of the contact relative to the location of the graphite step edge. These simulations reproduced a previously-reported current dip at step edges measured experimentally and the trend was explained by changes in both contact area and the distribution of distances between atoms in the interface. The novel approach reported here provides a foundation for future studies of the fundamental relationships between conductance, load and surface topography at the atomic scale.

Atomistic simulations of electrochemical metallization cells: mechanisms of ultra-fast resistance switching in nanoscale devices

We describe a new method that enables reactive molecular dynamics (MD) simulations of electrochemical processes and apply it to study electrochemical metallization cells (ECMs). The model, called EChemDID, extends the charge equilibration method to capture the effect of external electrochemical potential on partial atomic charges and describes its equilibration over connected metallic structures, on-the-fly, during the MD simulation. We use EChemDID to simulate resistance switching in nanoscale ECMs; these devices consist of an electroactive metal separated from an inactive electrode by an insulator and can be reversibly switched to a low-resistance state by the electrochemical formation of a conducting filament between electrodes. Our structures use Cu as the active electrode and SiO2 as the dielectric and have dimensions at the foreseen limit of scalability of the technology, with a dielectric thickness of approximately 1 nm. We explore the effect of device geometry on switching timescales and find that nanowires with an electroactive shell, where ions migrate towards a smaller inactive electrode core, result in faster switching than planar devices. We observe significant device-to-device variability in switching timescales and intermittent switching for these nanoscale devices. To characterize the evolution in the electronic structure of the dielectric as dissolved metallic ions switch the device, we perform density functional theory calculations on structures obtained from an EChemDID MD simulation. These results confirm the appearance of states around the Fermi energy as the metallic filament bridges the electrodes and show that the metallic ions and not defects in the dielectric contribute to the majority of those states.

Charge-transfer potentials for ionic crystals: Cauchy violation, LO-TO splitting, and the necessity of an ionic reference state

In this work, we study how including charge transfer into force fields affects the predicted elastic and vibrational Γ-point properties of ionic crystals, in particular those of rock salt. In both analytical and numerical calculations, we find that charge transfer generally leads to a negative contribution to the Cauchy pressure, P(C) ≡ C12 - C66, where C12 and C66 are elements of the elastic tensor. This contribution increases in magnitude with pressure for different charge-transfer approaches in agreement with results obtained with density functional theory (DFT). However, details of the charge-transfer models determine the pressure dependence of the longitudinal optical-transverse optical splitting and that for partial charges. These last two quantities increase with density as long as the chemical hardness depends at most weakly on the environment while experiments and DFT find a decrease. In order to reflect the correct trends, the charge-transfer expansion has to be made around ions and the chemical (bond) hardness has to increase roughly exponentially with inverse density or bond lengths. Finally, the adjustable force-field parameters only turn out meaningful, when the expansion is made around ions.