Abstract
In the context of ionic transport in solids, the variation of a migration barrier height under electric fields is traditionally assumed to be equal to the classical electric work of a point charge that carries the transport charge. However, how reliable is this phenomenological model and how does it fare with respect to Modern Theory of Polarization? In this work, we show that such a classical picture does not hold in general as collective dipole effects may be critical. Such effects are unraveled by an appropriate polarization decomposition and by an expression that we derive, which defines the equivalent polarization-work charge. The equivalent polarization-work charge is not equal neither to the transported charge, nor to the Born effective charge of the migrating atom alone, but it is defined by the total polarization change at the transition state. Our findings are illustrated by oxygen charged defects in MgO and in SiO2.
Similar content being viewed by others
Introduction
Ion motion in solids proceeds via jumps between different sites. In the absence of an external electric field, the probability that an ion will jump is proportional to the Boltzmann weight of the migration barrier. The field will, however, lower the barrier for motion along its direction, thus increasing the probability of movement and, therefore, the drift velocity of the ion.
Very recently, it has been formally demonstrated that the transported chargeâthe charge that is carried from one site to a symmetrically equivalent site (displacement of an atomic sublattice by a lattice vector)âis an integer number, thus providing a sound definition of the oxidation state of a migrating ion (see ref. 1). It might therefore appear natural to assume that the charge involved in the barrier height variation is the same charge as the transported one. This assumption was already implicit in Cabrera and Mottâs Theory of oxidation of metals2 (see also refs. 3,4,5), in which the following simple formula was accounting for the activation enthalpy change under electric fields:
WS is the electric work performed by a single ion of charge Q hopping a potential barrier that lies halfway (d/2) between its initial and final position (\(\frac{Qd}{2}\) is the dipole moment originated from the position change). In this early expression for WS, the electric field, E, was, in addition, assumed to be colinear with the jump coordinate. It is worth noting that the displacement, \(\frac{d}{2}\), is not a displacement of a whole lattice vector, as in the case of the transported charge.
Equation (1) has been widely adopted since its early formulation, and it is still nowadays applied in many different contexts, or for example, in the very active field of multiscale modeling6,7,8,9,10,11,12,13, or in the interpretations of experimental data related to oxidation refs. 14,15,16. Recent studies have started to question CabreraâMottâs theory of oxidation of metals in nanostructures17,18, but not Eq. (1). Even at the level of full first-principle calculations, discussions about the limit of the classical picture are only marginal. For instance, ref. 19 does not report any significant deviation of their results with respect to classical models in the case of MgO.
Materials are, however, composed of many interacting particles. Ionic motion occurs, in principle, in a multidimensional energy landscape in which the migration path up to the barrier might comprise configuration changes. It is, therefore, hard to believe that the electric work is, in general, only dependent on the dipole momentâthe position change times the chargeâof a single ion. A question naturally arises: is it mere coincidence or the signature of a deeper truth as in the aforementioned case of the oxidation state1?
While the classical picture is almost universally adopted in the literature, a few anomalous results have been reported, together with inconclusive attempts to explain them. For instance, ref. 20 revealed anomalous activation energies for the ionic conductivity in tantalum oxide. To fit the experimental data, the charge Q in Eq. (1) was multiplied by an effective parameter. More recently, because of discrepancies between the predicted ionic drift mobility in oxide-based memristors and the measured ones (discrepancies of several orders of magnitude in the ratio between the volatility and the switching time), it was proposed to correct Eq. (1), by including local field effects (the Lorentz field)21. Such corrections were, however, strongly criticized as lacking physical ground as commented and documented by Meuffels and Schroeder22.
Adopting the framework of Modern Theory of Polarization23,24,25 and decomposing the macroscopic polarization in terms of contributions from the migrating atom (MA)âthe hopping ionâand from the environment atoms, we derive an expression that allows a meaningful representation of the underlying many-particle screening mechanisms. This includes nontrivial contributions from the environment, and provides a definition of the equivalent polarization-work charge, that is, the charge to be multiplied by the MA displacement to obtain the polarization change. We then show that, in general, the equivalent polarization-work charge is not equal to the transported charge, it is not equal to the Born effective charge (BEC) of the MA, and it is not restricted to integer numbers. Our proposed new point of view demonstrates that the success of Eq. (1) can be traced back to subtle system-dependent and symmetry-related compensations between different effects. Environment contributions are crucial for understanding and modeling the polarization-work evolution up to the saddle point, with a corresponding global effect that is critically dependent on the material and the system. Our new point of view is illustrated through two examples: charged oxygen defects in an ionic metal oxide, MgO, and in an iono-covalent semiconductor oxide, SiO2.
Results
Equivalent polarization-work charge
In order to find a rational for the success of Eq. (1), we will exploit BECs and an appropriate polarization decomposition. Details on the derivation are provided in the Supplementary Note 2.
The thermodynamic potential, namely the electric enthalpy, H, related to the presence of an external electric field, at fixed atomic positions and up to the first order in the perturbation, is: \(H(\overrightarrow{E})={H}^{(0)}-\Omega \overrightarrow{P}\cdot \overrightarrow{E}\), where H(0) is the unperturbed potential, \(\overrightarrow{P}\) is the macroscopic polarization, and Ω is the volume of the system.
Along the migration path, the electric enthalpy varies proportionally to the variation of the polarization. This variation is called the polarization work (examples are provided in Supplementary note 1). The work from the initial, 0, to the saddle point, S, is therefore:
The BEC of an atom i, \({Z}_{i,\alpha \beta }^{* }\) is defined as the variation of the systemâs polarization, âPα, with respect to the variation of its position, âri,β: \({Z}_{i,\alpha \beta }^{* }=\frac{\partial {P}_{\alpha }}{\partial {r}_{i,\beta }}\). It is worth noting that the BEC is, in general, a nonsymmetric matrix. As a function of the BECs, the polarization variation up to the saddle point reads:
where λ is the reaction coordinate, N the total number of atoms, and \(d{\overrightarrow{r}}_{i}\) is the infinitesimal displacement vector of atom i.
At this point, the MA and environment contributions can be disentangled by using the BEC of the MA (an alternative decomposition that disentangles ionic and electronic contributions can be found in the Supplementary Note 3):
Equation (4) is still exact. The environment contribution can be efficiently computed as \({\mathrm{d}}{\overrightarrow{P}}_{{\rm{env}}}={\mathrm{d}}\overrightarrow{P}-{\mathrm{d}}{\overrightarrow{P}}_{\mathrm{MA}}\).
We will now perform a series of approximations to bridge between the Modern Theory of Polarization and the classical picture.
If the migration does not significantly change the chemical nature of the bonding between the species, including the MA, the BECs will remain almost constant along the whole reaction path. The BEC could, therefore, be taken out of the integrals in Eq. (4):
where \({\overrightarrow{d}}_{{\rm{MA}}}^{S}\) and \({\overrightarrow{d}}_{j}^{S}\) are, respectively, the integrated displacement vector of the MA and of environment atoms along the path from initial to the saddle point. The assumption of constant BECs along the reaction path is justified for the environment atoms for which chemical environment does not change significantly, while it might be more critical for the MA. If a more accurate polarization work is required, we expect that only a detailed evaluation of the MA BEC should be necessary.
Multiplying and dividing the environment contribution by \({\overrightarrow{d}}_{{\rm{MA}}}^{S}\), an expression for an equivalent polarization-work charge matrix can be derived. Component by component it reads
Because of the perturbation originated by the MA jump, it is more likely that uncompensated dipoles appear along the MA displacement direction. The major contribution to the work is, therefore, expected to come from the polarization component parallel to the MA displacement. Trivially, on the sole MA side, the larger dipole contribution will be along the direction of its displacement (jump), as it is along this direction that the MA moves the most. Similarly, MA neighbors are expected to be dragged or pushed along the same direction as a consequence of the Coulomb interaction. Assuming colinearity between the displacement of the MA and the infinitesimal polarization \({\mathrm{d}}\overrightarrow{P}\) and assuming that \({\mathrm{d}}\overrightarrow{P}\) keeps almost the same orientation along the path, up to the saddle point, we obtain an approximate form for the polarization work at the saddle point, that reminds Eq. (1)
with
Because of the colinearity assumption, θ is now the angle between the MA net displacement \({\overrightarrow{d}}_{{\rm{MA}}}\) and the field \(\overrightarrow{E}\). Equation (8) defines the equivalent polarization-work charge, Qeqv, to be plugged into Eq. (1). The sign is fixed by the sign of the macroscopic polarization.
Equation (6) clearly demonstrates that the variation of the barrier height, WS, and consequently of the drift velocity, under electric fields, is neither directly related to the transported charge nor to the BEC of the MA alone. Furthermore, it is clearly not restricted to integer numbers.
Charged defect migration in MgO and SiO2
Let us now see how our sample systems behave with respect to the classical theory, in light of our more complete picture.
For all cases, the projection of the macroscopic polarization variation along the displacement vector of the MA depends quasi-linearly on the MA displacement along the path, as shown in Fig. 1aâc (lower panels). The magnitude of the polarization projection is very close to the modulus of the total polarization, which also exhibits a quasi-linear dependency with dMA. The perpendicular components are, indeed, very small compared with the parallel one. Therefore, the latter should be very close to be colinear with the MA displacement vector. Details on how the projections are defined can be found in the Supplementary Note 2.
The top panels describe the atomic mechanism of the migration of a charged +2 oxygen vacancy, V\({}_{\mathrm{O}}^{2+}\), in MgO, and of a charged â2 oxygen interstitial, I\({}_{\mathrm{O}}^{2-}\), and a charged oxygen +2 vacancy, V\({}_{\mathrm{O}}^{2+}\), in SiO2. The remaining panels show the parallel (ÎPâ¥) and perpendicular projections (ÎPâ¥, ÎPÃ) of ÎP as a function of the net displacement of the migrating atom, dMA, along the reaction path, for, respectively, V\({}_{\mathrm{O}}^{2+}\) in MgO (a), I\({}_{\mathrm{O}}^{2-}\) (b), and V\({}_{\mathrm{O}}^{2+}\) (c) in SiO2. The parallel component corresponds to the projection along the MA displacement vector (\({\mathrm{d}}{\overrightarrow{r}}_{{\rm{MA}}}^{J}\)), which depends on the reaction path coordinates at a given point of the path. Further details on how the projections are defined can be found in Supplementary Note 2.
On highly symmetric, compact, and ionic materials, like MgO19,26 and similar rocksalts, because of the symmetry of the lattice, the movement of environment atoms results in a negligible net displacement from the minimum to the saddle point, \({\overrightarrow{d}}_{{\rm{env}}}^{S}\), see Fig. 1a (upper panel) and Fig. 2a. Furthermore, because of the strong ionic nature of the material, Mg and O environment atoms exhibit BECs, see upper panel of Fig. 3a, that are close to their corresponding oxidation level, +2 and â2, and that results in compensating dipoles along the whole path (the charges are opposite and the net displacements are equal and pointing in the same direction). Therefore, the environment contributions to the total dipole cancel out, and only the displacement and charge of the MA contribute to the polarization work and the variation of the migration barrier, see the bottom panel of Fig. 3a. In addition, because of the cubic symmetry, BEC matrices tend to be symmetric and isotropic with identical eigenvalues, see the upper panel of Fig. 3a. It can therefore be simplified to a scalar. The saddle point lies at half the distance between the initial (0) and final (F) configuration, and Eq. (7) simplifies into Eq. (1), with Qeqv equal to the MA-transported charge, Q. Therefore, the classical picture seems to hold because of symmetry-related/material-related compensations. As the vacancy moves in the opposite direction with respect to the MA, the polarization work associated with the defect has opposite sign, leading to Q â +2.
Projection of the total net displacement of environment atoms (for each environment species) along the MA displacement vector as a function of the MA net displacement for, respectively, V\({}_{\mathrm{O}}^{2+}\) in MgO (a), I\({}_{\mathrm{O}}^{2-}\) (b) and V\({}_{\mathrm{O}}^{2+}\), (c) in SiO2. The projection is computed as \({d}_{{\mathrm{env}},\parallel }^{(I)}=\mathop{\int}\nolimits_{0}^{I}| \mathop{\sum }\nolimits_{i}^{N-1}\frac{{\mathrm{d}}{\overrightarrow{r}}_{i}^{(I)}}{{\mathrm{d}}\lambda }{\mathrm{d}}\lambda \frac{{\mathrm{d}}{\overrightarrow{r}}_{{\rm{MA}}}^{(I)}}{| {\mathrm{d}}{\overrightarrow{r}}_{{\rm{MA}}}^{(I)}| }|\), where I tags a given point along the reaction path, i is the atom index, and N is the total number of atoms.
Top panels show the evolution of the eigenvalues and the trace (Tr(Z*)) of the BEC matrix \({Z}_{\alpha \beta }^{* }\) of the moving atom (MA) as function of its net displacement (dMA), for, respectively, V\({}_{\mathrm{O}}^{2+}\) in MgO (a, the MA is an oxygen atom), I\({}_{\mathrm{O}}^{2-}\) (b, the MA is an oxygen atom), and V\({}_{\mathrm{O}}^{2+}\) (c, the MA is a silicon atom) in SiO2. The inset in c shows the imaginary part of eigenvalue W1 (as the BEC is a real matrix but not symmetric, its diagonalization might lead to a real eigenvalue and two complex conjugate eigenvalues; see also ref. 38). Bottom panels display the evolution of the parallel projection of the total macroscopic polarization, ÎPâ¥, of the parallel projection of the MA contribution, ÎPMA,â¥, and of the environment contribution ÎPenv,â¥, as a function of the MA net displacement, dMA.
Net displacements related to environment atoms and BECs are rather more complex in SiO2 than in the simple cubic and ionic crystal, MgO, see Fig. 2b, c. SiO2 polymorphs (excluding high-pressure phases) are three-dimensional assemblies of corner-sharing SiO4 tetrahedra. Oxygen interstitials with charge â2 have the structure of two edge-sharing pentahedra (pentacoordination for the two silicon atoms). The migration of the oxygen interstitial involves the hopping of one of the oxygen atoms (MA) that defines the shared edge, see Fig. 1b (upper panel). For the case of the charge +2 oxygen vacancy, two separated units can be identified, both made from a back-projected positively charged silicon atom (each unit carries a +1 charge) weakly bonded to a three-coordinated back oxygen. The movement of one of these units involves the displacement of the back-projected silicon (MA) that gets attached to a different neighboring oxygen atom, see Fig. 1c (upper panel). Because of the tetrahedral site symmetry and the iono-covalent nature of the bonds, the BEC of oxygen and silicon atoms is strongly anisotropic and its trace does not correspond to the oxidation level of the species.
In the particular case of I\({}_{\mathrm{O}}^{2-}\), the BEC matrix is significantly nonisotropic, Fig. 3b (upper panel), with eigenvalue differences of up to 1.5e, as in the case of networking oxygen atoms. This anisotropy makes it crucial to explicitly account for the tensor nature of the BEC, and its corresponding matrixâvector product with the displacement vector, \({Z}^{* }\cdot {\overrightarrow{d}}^{S}\). Along the reaction path, the trace of the MA BEC remains almost constant with a signature of charge localization near the saddle point. This stronger localization contributes, however, small variations on the trace of the order of a few tens of e. Charge-transfer effects are, therefore, negligible. As a consequence, the MA BEC can be approximated by a constant matrix taking, for instance, the BEC at the initial configuration, Z*,0. Figure 2b clearly shows that in the saddle-point neighborhood, the MA displacement is combined with the recoil of mainly silicon environment atoms. The contribution to the polarization coming from the parallel component of the environment, ÎPenv,⥠is, indeed, very important and comparable to the contribution of the MA, ÎPMA,â¥, see Fig. 3b (bottom panel) and Table 1. If only the MA macroscopic polarization component is taken into account, the charge found from the ratio at the saddle point is of about â2.2e, in line with the classical ideas, as there is no charge transfer on the MA. The total macroscopic polarization work is, however, significantly higher, leading to a high negative equivalent polarization-work charge of â5.1e or of 4.9e, depending on whether the whole or the parallel component is considered. Such an anomalously high polarization work is originated from displacements of neighboring atoms, mainly of opposite charge, that is, Si atoms. This results in a significant increase of the polarization that apparently dresses the MA with a higher equivalent polarization-work charge, similar to a âQuasi-Particleâ/âQuasi-ionâ picture. In other words, the apparent âdressingâ is a many-particle effect that arises from the dielectric response of environment atoms that contributes to the overall macroscopic dipole moment.
Finally, for V\({}_{\mathrm{O}}^{2+}\), the displacements of oxygen and silicon environment atoms are not trivially compensating each other, see Fig. 2c as in the MgO case. From the minimum to the saddle point, O environment atoms recoil about twice the recoil distance of silicon atoms. However, as the trace of oxygen BECs is close to be half the opposite value of silicon BECs, the contributions to the total dipole of environment atoms cancel out, see also Fig. 3c (bottom panel) and Table 1. As such, only the BEC of the MA and its displacement contribute to the polarization work, therefore explaining the agreement between the calculated equivalent polarization-work charge, 3.5e, and the BEC of the MA.
Discussion
From the above analysis, we can draw a simple recipe to estimate the polarization work from a zero-field reaction path, short-cutting the burden of a full first-principles Berry Phase calculation over the whole path, at least as a first approximation. The macroscopic polarization change can be estimated assuming isotropic charges and a small path-curvature. Equation (8) then simplifies into a sum of species net displacements (initial and saddle point, only) times their respective oxidation state evaluated from the material isotropic BEC, 1/3Tr(Z*). If environment atoms net displacement at the saddle point times the oxidation state is close to zero, only the MA contribution is relevant, and Eq. (7) simplifies to the classical Eq. (1). Table 1 summarizes the change in the macroscopic polarization obtained at different levels of approximation. Our three examples show the possible different behaviors: (i) in V\({}_{\mathrm{O}}^{+2}\)@MgO, the classical picture would be adequate since the system is highly symmetric and environmental contribution cancels out, (ii) in I\({}_{\mathrm{O}}^{-2}\)@SiO2, the environment contribution is essential to obtain an accurate description and the classical picture fails completely, (iii) in V\({}_{\mathrm{O}}^{+2}\)@SiO2, the system has low symmetry and, although environment contribution cancels out, the tensor nature of the BEC is required for an accurate description, and a simplified version of Eq. (7) based on isotropic charges gives only semi-quantitative agreement.
In conclusion, we establish a clear distinction between transport charge and equivalent polarization-work charge beyond classical ideas, such as providing a solid ground to improve our understanding and modeling accuracy of electrochemical processes under electric fields. We show that, to understand the variation of the migration barrier under electric fields and the underlying microscopic dipole effects, it is necessary to go beyond the classical idea of a single-moving ion. The knowledge of nominal charges and zero-field reaction barriers alone is in general not sufficient to address/characterize ionic motion under electric fields. From the framework of the Modern Theory of Polarization, we demonstrate that the environment response plays, in general, a crucial role. In cases in which colinearity holds and in which BECs do not vary significantly along the reaction path, it is possible to generalize Eq. (1), still keeping a linear relation between a charge and a displacement, but with a revisited âchargeâ concept, that is, the equivalent polarization-work charge. The equivalent polarization-work charge is not equal to the transported charge or the BEC of the MA alone, and it is not restricted to integer numbers. The standard BEC charge definition becomes âdressedâ by the many-particle environment dipole response, leading to a semiclassical model of screened moving ions. It is worth noting that the âdressingâ is not a charge transfer but a pure dipole effect. Owing to the connection that we established between our revisited charge model and the Modern Theory of Polarization, the success of Eq. (1) and similar simplified classical pictures can be traced back to symmetry-related compensations between environment atoms that can only apply to specific systems. We expect, therefore, strong environment contributions to the ionic conductivity under electric fields in low-symmetry or low-density materials and heterostructures, where environment dipoles compensate each otherâfor instance, at surfaces and interfaces (cathodes, memristors, and OxRRAM as SiOx), in ceramics and glasses, in bulk phosphates and derivatives (Li-ion battery cathodes as LiFePO4), and in soft-matter solids. In particular, the challenging field of polymer-based electrolytes to be used in future flexible high-performing batteries27,28 could benefit from our framework in order to better understand the detailed atomic-scale mechanisms underlying the ionic conductivity in these materials. As the calculation of BECs and the polarization work is computationally demanding fully from first principles, in practice, it is possible to proceed step by step by first checking the environment net displacement, the local stoichiometry, and the nature of the bonding, and later by including (or not) each of the approximations that brings from the exact polarization work to our semiclassical approximation.
We hope that the present analysis will spur the inclusion of calculated BEC among the properties reported in popular open-source material database such as NOMAD29 or Material Clouds30 and will open the way to a more rational design of materials for electrochemistry.
Our findings will not only impact on theoretical/computational modeling, but also the way experimental measures are analyzed as high drift velocities might be tagged âanomalousâ only because of reminding old classical ideas.
Methods
The calculations presented in this work are all based on Density Functional Theory and Modern Theory of Polarization and the Berry Phase23,24 as implemented in Quantum ESPRESSO31,32. A PerdewâZunger33 exchange-correlation functional together with a Norm-conserving pseudopotential from the original QE library has been used. The valence wave functions are expanded into plane waves based on an energy cutoff of 80âRy (1088 eV)). The Brillouin zone is sampled only at Î k-point (https://www.quantum-espresso.org).
The SiO2, Quartz model is a supercell of 216 atoms with a group symmetry P3221. The model has been fully relaxed with a force and stress convergence threshold of 7âmeV/Ã . The initial and final defect configurations have been relaxed with the same force convergence threshold.
Reaction paths have been calculated by means of Climbing Image Nudged Elastic Band 34,35 (NEB). The path has been decomposed into nine frames. The convergence threshold for the NEB is 0.1âeV/à , and the distance between the frames is 0.3âà , while the mean displacement of the MA between consecutive frames is 0.4âà .
BECs are computed from the Berry Phase by means of finite differences with a real space grid of 0.1âà . In the calculation of the polarization, special care has been taken to refer all reported values to the same branch, see refs. 36,37.
Data availability
The software used for the present study (Quantum ESPRESSO 6.0) is released under GPL and is freely available at https://github.com/QEF/q-e. The pseudopotentials (O.pz-mt.UPF, Si.pz-vbc.UPF, Mg.pz-n-vbc.UPF) are available in the pseudopotential table at www.quantum-espresso.org under the sub-category âoriginal QE pp table.â Source data that support the findings of this study, including Quantum ESPRESSO inputs and outputs, and the post-processing script (python notebook) that computes BECs, displacements, and polarization changes, are available as a zip file as Supplementary information. Source data are provided with this paper.
References
Grasselli, F. & Baroni, S. Topological quantization and gauge invariance of charge transport in liquid insulators. Nat. Phys. https://doi.org/10.1038/s41567-019-0562-0 (2019).
Cabrera, N. & Mott, N. F. Theory of the oxidation of metals. Rep. Prog. Phys. 12, 163 (1949).
Mott, N. F. & Gurney, R. W. Electronic Processes in Ionic Crystals 2nd edn (Clarendon Press, 1953).
Lawless, K. R. The oxidation of metals. Rep. Prog. Phys. 37, 231â316 (1974).
Dignam, M. J. Ion transport in solids under conditions which include large electric fields. J. Phys. Chem. Solids 29, 249â260 (1968).
Sasikumar, K. et al. Evolutionary optimization of a charge transfer ionic potential model for Ta/Ta-Oxide heterointerfaces. Chem. Mater. 29, 3603â3614 (2017).
Martirosyan, K. S. & Zyskin, M. Reactive self-heating model of aluminum spherical nanoparticles. Appl. Phys. Lett. 102, 053112 (2013).
Raffone, F. & Cicero, G. Unveiling the fundamental role of temperature in RRAM Switching Mechanism by multiscale simulations. ACS Appl. Mater. Interfaces 10, 7512â7519 (2018).
Aldana, S. et al. An in-depth description of bipolar resistive switching in Cu/HfOx/Pt devices, a 3D kinetic monte carlo simulation approach. J. Appl. Phys. 123, 154501 (2018).
Menzel, S. Simulation and modeling of the switching dynamics in resistive switching devices. ECS Trans. 69, 19â32 (2015).
Dirkmann, S. & Mussenbrock, T. Resistive switching in memristive electrochemical metallization devices. AIP Adv. 7, 065006 (2017).
Sadi, T., Wang, L., Gerrer, L. & Asenov, A. Physical simulation of Si-based resistive random-access memory devices. In 2015 International Conference on Simulation of Semiconductor Processes and Devices (SISPAD) 385â388 (2015).
Sadi, T. et al. Investigation of resistance switching in SiOx RRAM cells using a 3D multi-scale kinetic monte carlo simulator. J. Phys. 30, 084005 (2018).
Han, L. et al. Interrogation of bimetallic particle oxidation in three dimensions at the nanoscale. Nat. Commun. 7, 13335 (2016).
Apolinário, A. et al. Modeling the growth kinetics of anodic TiO2 nanotubes. J. Phys. Chem. Lett. 6, 845â851 (2015).
Choksi, A. J., Lal, R. & Chandorkar, A. N. Growth kinetics of silicon dioxide on silicon in an inductively coupled rf plasma at constant anodization currents. J. Appl. Phys. 72, 1550â1557 (1992).
Pratt, A. et al. Enhanced oxidation of nanoparticles through strain-mediated ionic transport. Nat. Mater. 13, 26 (2013).
Liu, Y. et al. Enhanced oxidation resistance of active nanostructures via dynamic size effect. Nat. Commun. 8, 14459 (2017).
El-Sayed, A.-M., Watkins, M. B., Grasser, T. & Shluger, A. L. Effect of electric field on migration of defects in oxides: Vacancies and interstitials in bulk MgO. Phys. Rev. B 98, 064102 (2018).
Bean, C. P., Fisher, J. C. & Vermilyea, D. A. Ionic conductivity of tantalum oxide at very high fields. Phys. Rev. 101, 551â554 (1956).
Strukov, D. B. & Williams, R. S. Exponential ionic drift: fast switching and low volatility of thin-film memristors. Appl. Phys. A 94, 515â519 (2009).
Meuffels, P. & Schroeder, H. Comment on âexponential ionic drift: fast switching and low volatility of thin-film memristorsâ by D.B. Strukov and R.S. Williams in Appl. Phys. A (2009) 94, 515â519. Appl. Phys. A 105, 65 (2011).
Resta, R. Macroscopic polarization in crystalline dielectrics: the geometric phase approach. Rev. Mod. Phys. 66, 899â915 (1994).
Umari, P. & Pasquarello, A. Ab initio molecular dynamics in a finite homogeneous electric field. Phys. Rev. Lett. 89, 157602 (2002).
Nunes, R. W. & Gonze, X. Berry-phase treatment of the homogeneous electric field perturbation in insulators. Phys. Rev. B 63, 155107 (2001).
Youssef, M., Van Vliet, K. J. & Yildiz, B. Polarizing oxygen vacancies in insulating metal oxides under a high electric field. Phys. Rev. Lett. 119, 126002 (2017).
Tominaga, Y. Ion-conductive polymer electrolytes based on poly(ethylene carbonate) and its derivatives. Polym. J. 49, 291â299 (2017).
Liu, W. et al. Enhancing ionic conductivity in composite polymer electrolytes with well-aligned ceramic nanowires. Nat. Energy 2, 17035 (2017).
Nomad. https://nomad-coe.eu/ (2017).
Materialsclouds. https://www.materialscloud.org/home (2015).
Giannozzi, P. et al. Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. 21, 395502 (2009).
Giannozzi, P. et al. Advanced capabilities for materials modelling with Quantum ESPRESSO. J. Phys. 29, 465901 (2017).
Perdew, J. P. & Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B 23, 5048â5079 (1981).
Henkelman, G. & Jónsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J. Chem. Phys. 113, 9978â9985 (2000).
Henkelman, G., Uberuaga, B. P. & Jónsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 113, 9901â9904 (2000).
King-Smith, R. D. & Vanderbilt, D. Theory of polarization of crystalline solids. Phys. Rev. B 47, 1651â1654 (1993).
Vanderbilt, D. & King-Smith, R. D. Electric polarization as a bulk quantity and its relation to surface charge. Phys. Rev. B 48, 4442â4455 (1993).
Gonze, X., Allan, D. C. & Teter, M. P. Dielectric tensor, effective charges, and phonons in α-quartz by variational density-functional perturbation theory. Phys. Rev. Lett. 68, 3603â3606 (1992).
Acknowledgements
We are grateful to Stefano Baroni for useful discussions. High-performance computing resources have been provided by the Arctur company (https://hpc.arctur.si) and by CINECA super computing center (www.cineca.it, grant number ISCRA C HP10CQGOFM). N.S., L.G., and M.V. acknowledge the financial support from the Slovenian Research Agency (research core funding no. P2-0412). N.S., L.M.-S., S.d.G., L.G., and A.H. are active members of the Multiscale and Multi-Model Approach for Materials in Applied Science consortium (MAMMASMIAS consortium), and acknowledge the efforts of the consortium in fostering scientific collaboration.
Author information
Authors and Affiliations
Contributions
P.B. and N.R. proposed the topic. N.S. performed the DFT calculations and related analysis. N.S., L.M-S., and S.d.G. analyzed the results and wrote the paper. A.H., M.V., L.G., and B.S. helped to revise the paper. All authors discussed and commented on the results and on the paper.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Salles, N., Martin-Samos, L., de Gironcoli, S. et al. Collective dipole effects in ionic transport under electric fields. Nat Commun 11, 3330 (2020). https://doi.org/10.1038/s41467-020-17173-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-17173-w
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.