Search

Doped cerium oxide is known for its reduction properties that are utilized in catalytic applications as well as in thermochemical cycling to produce solar fuels. Upon reduction of the lattice, oxygen vacancies and polarons are formed leading to a highly concentrated solution of defects in which the interactions of defects cannot be neglected anymore. In this study, the effect of defect interactions on the free energy of reduction for doped ceria with composition Ce1-x-yRExZryO2-x/2-δ was investigated by large scale Metropolis Monte Carlo multi-stage sampling simulations based on first-principles calculations. The simulations allowed the prediction of the relation between oxygen partial pressure and non-stoichiometry for the highly interacting, non-ideal system. The results show that the non-ideality observed in experiments can be traced back to the interactions of defects and allow prediction of the reduction behavior for various dopant types, dopant fractions, temperatures, and non-stoichiometries.

Hydrated acceptor-doped barium zirconate is a well-investigated proton conductor. In the analysis of most experimental studies, an ideal defect model is applied to fit the measured hydration data and obtain corresponding enthalpies and entropies. However, the data show a distinct deviation from ideal behaviour and thus defect interactions cannot be neglected. In the present contribution, the thermodynamics of water uptake into the yttrium-doped bulk material are investigated on the microscopic level with regards to ionic defect interactions. Metropolis Monte Carlo simulations using interaction models from first-principles energy calculations are applied to obtain an estimation of the free energy of interaction. The present results indicate that the ionic defect interactions are the primary reason for the non-ideality observed in experiments with varying yttrium fraction, proton fraction, and temperature. The activity coefficient quotients for the mass action law are obtained, which connect the ideal and real model and are of relevance to data evaluation and theoretical calculations.

The combination of density functional theory and Monte Carlo simulations is a powerful approach for the atomistic modeling of defect transport in solid electrolytes. The present contribution introduces the MOCASSIN software (Monte Carlo for Solid State Ionics) for kinetic and Metropolis Monte Carlo simulations of crystalline materials. MOCASSIN combines model building, visualization, and simulation, aiming to provide accessible MC for end users. Developed for the investigation of solid electrolytes, MOCASSIN is ideal for screening common variation parameters, such as temperature and doping fraction. The input effort is minimized using space groups for processing symmetry. The graphical interface for model building allows complex model input, including multiple mobile species, multiple migration paths, small polaron hopping, vehicle movements, multiple complex migration mechanisms, and custom interaction clusters. The software is provided free of charge for noncommercial usage., (© 2020 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

Acceptor-doped barium zirconate is a promising proton-conducting oxide for various applications, for example, electrolysers, fuel cells or methane-conversion cells. Despite many experimental and theoretical investigations there is, however, only a limited understanding as to how to connect the complex microscopic proton motion and the macroscopic proton conductivity for the full range of acceptor levels, from diluted acceptors to concentrated solid solutions. Here we show that a combination of density functional theory calculations and kinetic Monte Carlo simulations enables this connection. At low concentrations, acceptors trap protons, which results in a decrease of the average proton mobility. With increasing concentration, however, acceptors form nanoscale percolation pathways with low proton migration energies, which leads to a strong increase of the proton mobility and conductivity. Comparing our simulated proton conductivities with experimental values for yttrium-doped barium zirconate yields excellent agreement. We then predict that ordered dopant structures would not only strongly enhance the proton conductivities, but would also enable one- or two-dimensional proton conduction in barium zirconate. Finally, we show how the properties of other dopants influence the proton conductivity.

NASICON-structured materials of the composition Li 1+x Al x Ti 2-x (PO 4 ) 3 are regarded as solid electrolytes with high Li-ion conductivity applicable in all solid-state batteries. In this study we investigate the migration paths of constituting ions and monovalent charge carriers including K + , Na + , and, H + . The results proof that Li is the most mobile species in the investigated composition and that the formation of intrinsic defects is unlikely. In addition, we find surprisingly low migration energy for oxygen vacancies in the structure of the dedicated Li-ion conductor.

The lattice constant has a crucial effect on the defect chemistry and defect kinetics in solid state materials. However, within density functional theory, some functionals perform badly in reproducing the experimental lattice constant. In this study, energies of defect formation, interaction and migration in the model system ceria were calculated for different lattice constants to investigate the impact on the energies. The GGA+U functional in the PBE and PBEsol parametrization as well as the hybrid functional HSE06 were applied and results are compared among these three commonly applied functionals. The results suggest a strong influence of the lattice constant on the energies especially regarding oxygen ion migration. This influence has an impact on the accurate prediction of defect properties from first principles but can also be utilized for specific tailoring of material properties by chemo-mechanical design. In addition, the issue of the correct lattice constant, which should be used in the defect calculations, is discussed in this paper.

Materials with high oxygen ion conductivity and low electronic conductivity are required for electrolytes in solid oxide fuel cells (SOFC) and high-temperature electrolysis (SOEC). A potential candidate for the electrolytes, which separate oxidation and reduction processes, is rare-earth doped ceria. The prediction of the ionic conductivity of the electrolytes and a better understanding of the underlying atomistic mechanisms provide an important contribution to the future of sustainable and efficient energy conversion and storage. The central aim of this paper is the detailed investigation of the relationship between defect interactions at the microscopic level and the macroscopic oxygen ion conductivity in the bulk of doped ceria. By combining ab initio density functional theory (DFT) with Kinetic Monte Carlo (KMC) simulations, the oxygen ion conductivity is predicted as a function of the doping concentration. Migration barriers are analyzed for energy contributions, which are caused by the interactions of dopants and vacancies with the migrating oxygen vacancy. We clearly distinguish between energy contributions that are either uniform for forward and backward jumps or favor one migration direction over the reverse direction. If the presence of a dopant changes the migration energy identically for forward and backward jumps, the resulting energy contribution is referred to as blocking. If the change in migration energy due to doping is different for forward and backward jumps of a specific ionic configuration, the resulting energy contributions are referred to as trapping. The influence of both effects on the ionic conductivity is analyzed: blocking determines the dopant fraction where the ionic conductivity exhibits the maximum. Trapping limits the maximum ionic conductivity value. In this way, a deeper understanding of the underlying mechanisms determining the influence of dopants on the ionic conductivity is obtained and the ionic conductivity is predicted more accurately. The detailed results and insights obtained here for doped ceria can be generalized and applied to other ion conductors that are important for SOFCs and SOECs as well as solid state batteries.

Defects and their interactions play a crucial role for the properties of solid state materials. While energies of defect formation and defect interaction are commonly calculated from first principles, calculations of corresponding vibrational entropies in solids are less common due to the high computational demand. This is in particular true for defect pairs as the existence of multiple defects lowers the symmetry of the cell. Here, we determine the entropies of association of defects in doped cerium oxide by density functional theory and phonon calculations. We estimate the importance of the entropic contribution in the constant volume and constant pressure case and address the limitations of the calculations. The small values of the calculated entropies render the accurate determination difficult. However, even small entropy contributions could have a significant influence on defect interactions and ordering at high temperature.

The rate of oxygen ion jumps in a solid oxide depends not only on the activation energy but also on the pre-exponential factor of diffusion. In order to allow a fully ab initio prediction of the oxygen ion conductivity in pure and samarium doped ceria, we calculated the attempt frequency for an oxygen ion jump from first principles combining DFT+U, the NEB method, phonon calculations and the transition state theory. Different definitions of the jump attempt frequency are presented. The equivalence of the Eyring and the Vineyard method is shown without restriction to the Gamma point. Convergence checks of the phonon mesh reveal that the common reduction to the Gamma point is not sufficient to calculate the attempt frequency. Calculations of Sm doped ceria revealed an increase of the prefactor. The attempt frequency for the constant pressure case in quasi-harmonic approximation is larger than the attempt frequency at constant volume in harmonic approximation. The calculated electronic energies, enthalpies and entropies of migration are in agreement with the experimental diffusion coefficients and activation energies.

We investigate the interaction and distribution of defects in doped non-stoichometric ceria Ce1-xRExO2-x/2-δ (with RE = Lu, Y, Gd, Sm, Nd, and La) by combining DFT+U calculations and Monte Carlo simulations. The concentrated solution of defects in ceria is described by the pair interactions of dopant ions, oxygen vacancies, and small polarons. The calculated interaction energies for polarons and oxygen vacancies are in agreement with experimental results and previously reported calculations. Simulations reveal that in thermodynamic equilibrium the configurational energy decreases with increasing non-stoichiometry as well as increasing dopant fraction similar to the observed behavior of the enthalpy of reduction in experiments. This effect is attributed to the attractive interaction of oxygen vacancies with polarons and dopant ions.

We investigate the dopant distribution and its influence on the oxygen ion conductivity of ceria doped with rare earth oxides by combining density functional theory and Monte Carlo simulations. We calculate the association energies of dopant pairs, oxygen vacancy pairs and between dopant ions and oxygen vacancies by means of DFT + U including finite size corrections. The cation coordination numbers from ensuing Metropolis Monte Carlo simulations show remarkable agreement with experimental data. Combining Metropolis and Kinetic Monte Carlo simulations we find a distinct dependence of the ionic conductivity on the dopant distribution and predict long term degradation of electrolytes based on doped ceria.

We calculate entropies of formation for fully charged point defects, including the small polaron Ce'(Ce), in undoped fluorite-structured ceria by means of density functional theory in the GGA + U approximation. We discuss the behaviour of the entropy for the constant volume and the constant pressure case. Our results for constant pressure (p = 0) suggest that the change in volume, due to the formation of defects, dominates the entropy of formation. From the individual entropies of formation the entropies of Frenkel, anti-Frenkel and Schottky disorder as well as the entropy of reduction of ceria are obtained. At temperatures of about 1000 K the entropic contributions to the Gibbs energy are up to 0.9 eV per defect and thus are no longer negligible. For our calculated entropy of reduction of about 17 kB we find a remarkable agreement with experimental data from the literature.

The contralesional primary motor cortex (M1) has been suggested to be involved in the motor recovery after mirror therapy, but whether the ipsilesional M1 is influenced by the contralesional M1 via transcallosal interhemispheric inhibition (IHI) is still unclear. The present study investigated the change of IHI as well as the intracortical inhibition and intracortical facilitation of both M1 induced by training in a mirror with the use of transcranial magnetic stimulation (TMS). In this 2 × 2 factorial design (time × group), healthy subjects exercised standardized motor skills with their right hand on four consecutive days. Either a mirror (mirror group) or a board (control group) was positioned between their hands. Before and after training TMS was applied along with training tests of both hands. Tests were the same motor skills exercised daily by both groups. Tests of the untrained left hand improved significantly more in the mirror group than in the control group after training (P = 0.02) and showed a close correlation with an increase of intracortical inhibition of M1(left). IHI did not show any difference between investigation time points and groups. The present study confirms the previous suggestion of the involvement of the "contralesional" left-side (ipsilateral to the hand behind the mirror) M1 after mirror therapy, which is not mediated by IHI. Even with the same motor skill training (both groups performed same motor skills) but with different visual information, different networks are involved in training-induced plasticity.

Background: Previous studies demonstrated a posterior shift of activation toward the primary sensory cortex (S1) following stroke; however, any relationship between this posterior shift and clinical outcome measures for the affected hand function were unclear., Objective: The authors investigated the possible role of S1 in motor recovery., Methods: Assuming that previous studies examined inhomogeneous groups of patients, the authors selected participants with chronic stroke who had moderate hand paresis, normal sensory examination and somatosensory-evoked potentials, and no lesion within the S1, thalamus, or brain stem. Constraint-induced movement therapy (CIMT) was used to train the impaired hand. To relate fMRI (functional MRI) activation changes from baseline to post-CIMT, a correlation analysis was performed with changes of the Wolf Motor Function Test (WMFT) as a test for the hand function., Results: A close relationship was found between increases in hand function and peak changes in activation within the ipsilesional S1. With a better outcome, greater increases in activation within the S1 were evident (P < .03; r = 0.73)., Conclusion: In selected patients, the sensory network influences training-induced motor gains. This predictive knowledge of plasticity when applying CIMT may suggest strategies to enhance the effect of therapy, such as the addition of electrical stimulation to enhance S1 excitability.

Background: Stroke and mortality rates in patients with transient ischemic attack (TIA) differ widely between community-based studies and research cohorts. Our aim therefore was to provide a reliable estimate for TIA patients treated in German neurology departments with an acute stroke unit., Methods: A total of 1951 consecutively admitted TIA patients were prospectively documented in 13 centers and 1480 (75.9%) gave consent for long-term follow-up. During a mean follow-up of 23.4 months, we assessed recurrent cerebrovascular events and cause of death in 1448 patients via standardized telephone interview including confirmation of endpoint events by the treating physician., Results: Overall 94 patients (6.5%) suffered a stroke and 118 patients (8.1%) died, 21 due to stroke. The Kaplan-Meier estimate for stroke during the first year was 4.4% (95% CI 3.2-5.6%) which corresponds to a relative risk of 9.5 (95% CI 7.4-12.3) compared to the population-based stroke incidence in Germany. The annual rates after the first year were 2.2% (95% CI 1.7-2.7%) for stroke and 3.2% (95% CI 2.7-3.8%) for death. Independent predictors for stroke during follow-up were age and previous cerebrovascular events. The ABCD(2) score did not provide any meaningful prediction of stroke risk at 90 days., Conclusion: While the in-hospital risk of stroke was low, long-term stroke rates in our well-defined multicenter hospital-based cohort were comparable to a large randomized trial. In patients with a well-established diagnosis of TIA, only age and previous cerebrovascular events seem to constitute independent predictors for stroke during long-term follow-up.