Nicolas Onofrio

Stability and migration of small copper clusters in amorphous dielectrics

We use density functional theory (DFT) to study the thermodynamic stability and migration of copper ions and small clusters embedded in amorphous silicon dioxide. We perform the calculations over an ensemble of statistically independent structures to quantify the role of the intrinsic atomic-level variability in the amorphous matrix affect the properties. The predicted formation energy of a Cu ion in the silica matrix is 2.7 ± 2.4 eV, significantly lower the value for crystalline SiO2. Interestingly, we find that Cu clusters of any size are energetically favorable as compared to isolated ions; showing that the formation of metallic clusters does not require overcoming a nucleation barrier as is often assumed. We also find a broad distribution of activation energies for Cu migration, from 0.4 to 1.1 eV. This study provides insights into the stability of nanoscale metallic clusters in silica of interest in electrochemical metallization cell memories and optoelectronics.

Voltage equilibration for reactive atomistic simulations of electrochemical processes

We introduce electrochemical dynamics with implicit degrees of freedom (EChemDID), a model to describe electrochemical driving force in reactive molecular dynamics simulations. The method describes the equilibration of external electrochemical potentials (voltage) within metallic structures and their effect on the self-consistent partial atomic charges used in reactive molecular dynamics. An additional variable assigned to each atom denotes the local potential in its vicinity and we use fictitious, but computationally convenient, dynamics to describe its equilibration within connected metallic structures on-the-fly during the molecular dynamics simulation. This local electrostatic potential is used to dynamically modify the atomic electronegativities used to compute partial atomic changes via charge equilibration. Validation tests show that the method provides an accurate description of the electric fields generated by the applied voltage and the driving force for electrochemical reactions. We demonstrate EChemDID via simulations of the operation of electrochemical metallization cells. The simulations predict the switching of the device between a high-resistance to a low-resistance state as a conductive metallic bridge is formed and resistive currents that can be compared with experimental measurements. In addition to applications in nanoelectronics, EChemDID could be useful to model electrochemical energy conversion devices.

Molecular dynamic simulation of tip-polymer interaction in tapping-mode atomic force microscopy

We present a molecular dynamic study of the interaction between an amorphous silica tip (SiO2) and an amorphous poly-(methyl-methacrylate) substrate under conditions relevant for tapping-mode atomic force microscopy. To capture the actual dynamics of the tip, we use the dynamic contact simulation method [Kim et al., J. Appl. Phys. 112, 094325 (2012)]. We obtain force-displacement relationships both for neat polymer substrates and a sample with a sub-surface nanotube and extract the local stiffness and energy dissipation per cycle. The simulations capture non-trivial aspects of the interaction that originate from the viscoelastic nature of the polymer including an increase in repulsive interaction force during approach with tip velocity and an increase in adhesion during retraction with decreasing tip velocity. Scans of local stiffness and dissipation over the samples reveal intrinsic variability in the amorphous polymer but also the effect of local surface topography on the extracted properties as well as the ability of the method to detect a sub-surface nanotube. This insight and quantitative data should be valuable to interpret the results of atomic force microscopy studies.

First principles investigation of copper and silver intercalated molybdenum disulfide

We characterize the energetics and atomic structures involved in the intercalation of copper and silver into the van der Waals gap of molybdenum disulfide as well as the resulting ionic and electronic transport properties using first-principles density functional theory. The intercalation energy of systems with formula (Cu,Ag)xMoS2 decreases with ion concentration and ranges from 1.2 to 0.8 eV for Cu; Ag exhibits a stronger concentration dependence from 2.2 eV for x = 0.014 to 0.75 eV for x = 1 (using the fcc metal as a reference). Partial atomic charge analysis indicates that approximately half an electron is transferred per metallic ion in the case of Cu at low concentrations and the ionicity decreases only slightly with concentration. In contrast, while Ag is only slightly less ionic than Cu for low concentrations, charge transfer reduces significantly to approximately 0.1 e for x = 1. This difference in ionicity between Cu and Ag correlates with their intercalation energies. Importantly, the predicted...

The dynamics of copper intercalated molybdenum ditelluride

Layered transition metal dichalcogenides are emerging as key materials in nanoelectronics and energy applications. Predictive models to understand their growth, thermomechanical properties, and interaction with metals are needed in order to accelerate their incorporation into commercial products. Interatomic potentials enable large-scale atomistic simulations connecting first principle methods and devices. We present a ReaxFF reactive force field to describe molybdenum ditelluride and its interactions with copper. We optimized the force field parameters to describe the energetics, atomic charges, and mechanical properties of (i) layered MoTe2, Mo, and Cu in various phases, (ii) the intercalation of Cu atoms and small clusters within the van der Waals gap of MoTe2, and (iii) bond dissociation curves. The training set consists of an extensive set of first principles calculations computed using density functional theory (DFT). We validate the force field via the prediction of the adhesion of a single layer MoTe2 on a Cu(111) surface and find good agreement with DFT results not used in the training set. We characterized the mobility of the Cu ions intercalated into MoTe2 under the presence of an external electric field via finite temperature molecular dynamics simulations. The results show a significant increase in drift velocity for electric fields of approximately 0.4 V/Å  and that mobility increases with Cu ion concentration.

Modeling resistive switching materials and devices across scales

Resistance switching devices based on electrochemical processes have attractive significant attention in the field of nanoelectronics due to the possibility of switching in nanosecond timescales, miniaturization to tens of nanometer and multi-bit storage. Their deceptively simple structures (metal-insulator-metal stack) hide a set of complex, coupled, processes that govern their operation, from electrochemical reactions at interfaces, diffusion and aggregation of ionic species, to electron and hole trapping and Joule heating. A combination of experiments and modeling efforts are contributing to a fundamental understanding of these devices, and progress towards a predictive understanding of their operation is opening the possibility for the rational optimization. In this paper we review recent progress in modeling resistive switching devices at multiple scales; we briefly describe simulation tools appropriate at each scale and the key insight that has been derived from them. Starting with ab initio electronic structure simulations that provide an understanding of the mechanisms of operation of valence change devices pointing to the importance of the aggregation of oxygen vacancies in resistance switching and how dopants affect performance. At slightly larger scales we describe reactive molecular dynamics simulations of the operation of electrochemical metallization cells. Here the dynamical simulations provide an atomic picture of the mechanisms behind the electrochemical formation and stabilization of conductive metallic filaments that provide a low-resistance path for electronic conduction. Kinetic Monte Carlo simulations are one step higher in the multiscale ladder and enable larger scale simulations and longer times, enabling, for example, the study of variability in switching speed and resistance. Finally, we discuss physics-based simulations that accurately capture subtleties of device behavior and that can be incorporated in circuit simulations.

Novel doping alternatives for single-layer transition metal dichalcogenides

Successful doping of single-layer transition metal dichalcogenides (TMDs) remains a formidable barrier to their incorporation into a range of technologies. We use density functional theory to study doping of molybdenum and tungsten dichalcogenides with a large fraction of the periodic table. An automated analysis of the energetics, atomic and electronic structure of thousands of calculations results in insightful trends across the periodic table and points out promising dopants to be pursued experimentally. Beyond previously studied cases, our predictions suggest promising substitutional dopants that result in p-type transport and reveal interesting physics behind the substitution of the metal site. Doping with early transition metals (TMs) leads to tensile strain and a significant reduction in the bandgap. The bandgap increases and strain is reduced as the d-states are filled into the mid TMs; these trends reverse as we move into the late TMs. Additionally, the Fermi energy increases monotonously as the ...

Exceptional catalytic effects of black phosphorus quantum dots in shuttling-free lithium sulfur batteries

Lithium sulfur batteries with high energy densities are promising next-generation energy storage systems. However, shuttling and sluggish conversion of polysulfides to solid lithium sulfides limit the full utilization of active materials. Physical/chemical confinement is useful for anchoring polysulfides, but not effective for utilizing the blocked intermediates. Here, we employ black phosphorus quantum dots as electrocatalysts to overcome these issues. Both the experimental and theoretical results reveal that black phosphorus quantum dots effectively adsorb and catalyze polysulfide conversion. The activity is attributed to the numerous catalytically active sites on the edges of the quantum dots. In the presence of a small amount of black phosphorus quantum dots, the porous carbon/sulfur cathodes exhibit rapid reaction kinetics and no shuttling of polysulfides, enabling a low capacity fading rate (0.027% per cycle over 1000 cycles) and high areal capacities. Our findings demonstrate application of a metal-free quantum dot catalyst for high energy rechargeable batteries.Lithium sulfur batteries are promising for next-generation energy storage, but are hindered by polysulfide shuttle effects. Here the authors use black phosphorus quantum dots to adsorb and catalyze the conversion of lithium polysulfides to lithium sulfide, achieving low capacity fade and high sulfur loading.

Building Predictive Chemistry Models

Density Functional Theory (DFT) simulations allow for sophisticated modeling of chemical interactions, but the extreme computational cost makes it inviable for large scale applications. Molecular dynamics models, specifically ReaxFF, can model much larger simulations with greater speed, but with lesser accuracy. The accuracy of ReaxFF can be improved by comparing predictions of both methods and tuning ReaxFF’s parameters. Molecular capabilities of ReaxFF were gauged by simulating copper complexes in water over a 200 ps range, and comparing energy predictions against ReaxFF. To gauge solid state capabilities, volumetric strain was applied to simulated copper bulk and the strain response functions used to predict elastic constants, which were then compared against experimental data and ReaxFF predictions. Results suggest ReaxFF’s predictions are fairly robust, making it useful for molecular simulations. Training ReaxFF with this data can improve the accuracy of molecular dynamics simulations, providing wider application of molecular modeling software. KEYWORDS Molecular Dynamics (MD), Density Functional Theory (DFT), ReaxFF, copper complexes, bulk modulus REFERENCES Van Duin, Adri, Siddharth Dasgupta, Francois Lorant, and William Goddard III. “ReaxFF: A Reactive Force Field for Hydrocarbons”. J. Phys. Chem. A. 2001. web. Van Duin, Adri, Vayacheslav Bryantsec, Mamadou Diallo, William Goddard, Obaidur Rahaman, Douglas Doren, David Raymand, Kersti Hermansson. “Development and Validation of a ReaxFF Reactive Force Field for Cu Cation/Water Interactions and Copper Metal/Metal Oxide/Metal Hydroxide Condensed Phases”. J. Phys. Chem. A. 2010. web.

Molecular Exploration Tool

Density Functional Theory (DFT) which is based on quantum mechanics theory has been broadly used to compute the energy and the structure of molecules and solids. However, the DFT method is limited when running calculations for a large system and only thousands of atoms can be solved. Alternatively, Molecular Dynamics (MD) simulation can be used to investigate the properties of the atomic system for large systems in the classical mechanics approximation. When running the MD simulation, the electronic structure is approximated by Force Fields (FF) which can be parameterized against DFT calculations. Nevertheless, the accuracy of the MD results and the FF is suspicious for the scientists because of the variety and complexity of the FF. Hence, a free web-browser based tool has been developed to allow the user upload a force field, run MD simulations and compare the results with the DFT calculations. Users can select desired molecules and solids in the database, run MD simulation, plot the corresponding energies and visualize the atomic structures. So that users can find out if they can trust the FF results according to the comparison with DFT calculations.