Logo: to the web site of Uppsala University

We present a method for large-scale DFT-based screening of ion diffusion in crystalline solids. This is accomplished by extending the Ionic TuTraSt method to sample the potential energy surface by using single-point DFT calculations. To drastically reduce the number of single-point DFT calculations, symmetry, interpolation, and exclusion of high-energy regions are employed. This approach is tested on a large data set of solid-state Li-ion conductors, for which the interpolation and high-energy exclusion are optimized to balance computational efficiency and accuracy of the obtained diffusion properties. Furthermore, the developed workflow is validated by comparison with ab initio molecular dynamics (AIMD) simulations on a set of known Li-ion superconducting materials.

We present an efficient method to compute diffusion coefficients of multiparticle systems with strong interactions directly from the geometry and topology of the potential energy field of the migrating particles. The approach is tested on Li-ion diffusion in crystalline inorganic solids, predicting Li-ion diffusion coefficients within 1 order of magnitude of molecular dynamics simulations at the same level of theory while being several orders of magnitude faster. The speed and transferability of our workflow make it well-suited for extensive and efficient screening studies of crystalline solid-state ion conductor candidates and promise to serve as a platform for diffusion prediction even up to the density functional level of theory.

Solid-state composite electrolytes have arisen as one of the most promising materials classes for next-generation Li-ion battery technology. These composites mix ceramic and solid-polymer ion conductors with the aim of combining the advantages of each material. The ion-transport mechanisms within such materials, however, remain elusive. This knowledge gap can to a large part be attributed to difficulties in studying processes at the ceramic–polymer interface, which are expected to play a major role in the overall ion transport through the electrolyte. Computational efforts have the potential of providing significant insight into these processes. One of the main challenges to overcome is then to understand how a sufficiently robust model can be constructed in order to provide reliable results. To this end, a series of molecular dynamics simulations are here carried out with a variation of certain structural (surface termination and polymer length) and pair potential (van der Waals parameters and partial charges) models of the Li₇La₃Zr₂O₁₂ (LLZO) poly(ethylene oxide) (PEO) system, in order to test how sensitive the outcome is to each variation. The study shows that the static and dynamic properties of Li-ion are significantly affected by van der Waals parameters as well as the surface terminations, while the thickness of the interfacial region – where the structure–dynamic properties are different as compared to the bulk-like regime – is the same irrespective of the simulation setup.

The structure and dynamics of a material are essentially determined by the complex combination of potential energy landscapes experienced by the individual atoms in the system. In turn, valuable information on the properties of the material is encoded in the shapes of the potential energy landscape. For example, configurations of particles within a solid are determined by the shapes and presence of energetic basins, and the self-diffusion of mobile particles is defined by the geometry of how these energetic basins are connected to form paths. Understanding diffusion processes in solids at the atomistic scale is crucial for many important applications such as predicting Li-ion conduction through a solid-state battery cell or membranes for separation processes including carbon capture and water purification. While modeling can facilitate such understanding, there are still many challenges to overcome in terms of reaching relevant length and time scales that capture the complexity of the material. In this Perspective, we will discuss state-of-the-art modeling methods for mass transport inside a solid-state material and how they relate to the geometry of the potential energy landscape. We believe that approaching diffusion from a geometrical standpoint offers great promise in advancing modeling methodologies while yielding a better understanding of the structure-dynamic properties relationship and rate-limiting processes.

The high pressure structural behavior of H-2 and Ne clathrate hydrates with approximate composition H-2/Ne center dot~4H(2)O and featuring cubic structure II (CS-II) was investigated by neutron powder diffraction using the deuterated analogues at ~95 K. CS-II hydrogen hydrate transforms gradually to isocompositional C-1 phase (filled ice II) at around 1.1 GPa but may be metastably retained up to 2.2 GPa. Above 3 GPa a gradual decomposition into C-2 phase (H-2 center dot H2O, filled ice I-c) and ice VIII' takes place. Upon heating to 200 K the CS-II to C-1 transition completes instantly whereas C-1 decomposition appears sluggish also at 200 K. C-1 was observed metastably up to 8 GPa. At 95 K C-1 and C-2 hydrogen hydrate can be retained below 1 GPa and yield ice II and ice I-c, respectively, upon complete release of pressure. In contrast, CS-II neon hydrate undergoes pressure-induced amorphization at 1.9 GPa, thus following the general trend for noble gas clathrate hydrates. Upon heating to 200 K amorphous Ne hydrate crystallizes as a mixture of previously unreported C-2 hydrate and ice VIII'.

The high-pressure structural behavior of the noble gas (Ng) clathrate hydrates Ar center dot 6.5 H2O and Xe center dot 7.2 H2O featuring cubic structures II and I, respectively, was investigated by neutron powder diffraction (using the deuterated analogues) at 95 K. Both hydrates undergo pressure-induced amorphization (PIA), indicated by the disappearance of Bragg diffraction peaks, but at rather different pressures, at 1.4 and above 4.0 GPa, respectively. Amorphous Ar hydrate can be recovered to ambient pressure when annealed at >1.5 GPa and 170 K and is thermally stable up to 120 K. In contrast, it was impossible to retain amorphous Xe hydrate at pressures below 3 GPa. Molecular dynamics (MD) simulations were used to obtain general insight into PIA of Ng hydrates, from Ne to Xe. Without a guest species, both cubic clathrate structures amorphize at 1.2 GPa, which is very similar to hexagonal ice. Filling of large-sized H cages does not provide stability toward amorphization for structure II, whereas filled small-sized dodecahedral D cages shift PIA successively to higher pressures with increasing size of the Ng guest. For structure I, filling of both kinds of cages, large-sized T and small-sized D, acts to stabilize in a cooperative fashion. Xe hydrate represents a special case. In MD, disordering of the guest hydration structure is already seen at around 2.5 GPa. However, the different coordination numbers of the two types of guests in the crystalline cage structure are preserved, and the state is shown to produce a Bragg diffraction pattern. The experimentally observed diffraction up to 4 GPa is attributed to this semicrystalline state.

Both polyesters and polycarbonates have been proposed as alternatives to polyethers as host materials for future polymer electrolytes for solid-state lithium-ion batteries. While being comparatively similar functional groups, the electron density on the coordinating carbonyl oxygen is different, thereby rendering different coordinating strength towards lithium ions. In this study, the transport properties of poly(epsilon-caprolactone) and poly(trimethylene carbonate) as well as random copolymers of systematically varied composition of the two have been investigated, in order to better elucidate the role of the coordination strength. The cationic transference number, a property well-connected with the complexing ability of the polymer, was shown to depend almost linearly on the ester content of the copolymer, increasing from 0.49 for the pure poly(epsilon-caprolactone) to 0.83 for pure poly(trimethylene carbonate). Contradictory to the transference number measurements that suggest a stronger lithium-to-ester coordination, DFT calculations showed that the carbonyl oxygen in the carbonate coordinates more strongly to the lithium ion than that of the ester. FT-IR measurements showed the coordination number to be higher in the polyester system, resulting in a higher total coordination strength and thereby resolving the paradox. This likely originates in properties that are specific of polymeric solvent systems, e.g. steric properties and chain dynamics, which influence the coordination chemistry. These results highlight the complexity in polymeric systems and their ion transport properties in comparison to low-molecular-weight analogues, and how polymer structure and steric effects together affect the coordination strength and transport properties.

Metal-organic frameworks (MOFs) are promising materials for the photocatalytic H-2 evolution reaction (HER) from water. To find the optimal MOF for a photocatalytic HER, one has to consider many different factors. For example, studies have emphasized the importance of light absorption capability, optical band gap, and band alignment. However, most of these studies have been carried out on very different materials. In this work, we present a combined experimental and computation study of the photocatalytic HER performance of a set of isostructural pyrene-based MOFs (M-TBAPy, where M = Sc, Al, Ti, and In). We systematically studied the effects of changing the metal in the node on the different factors that contribute to the HER rate (e.g., optical properties, the band structure, and water adsorption). In addition, for Sc-TBAPy, we also studied the effect of changes in the crystal morphology on the photocatalytic HER rate. We used this understanding to improve the photocatalytic HER efficiency of Sc-TBAPy, to exceed the one reported for Ti-TBAPy, in the presence of a co-catalyst.

Clathrate hydrates with the cubic structure II (CS-II) form typically with large guest molecules, such as tetrahydrofuran, trimethylamine oxide, or propane. However, CS-II is also realized for argon hydrate despite the comparatively small van der Waals diameter of the guest (around 3.8 angstrom). Here, the structure of deuterated argon hydrate was studied at ambient pressure in the temperature range 20-95 K using neutron diffraction and comparing natural Ar with Ar-36, which scatters neutrons more than 13 times more efficiently. The procedure allowed to unambiguously establish the positional disorder within the large cages of CS-H, while simultaneously refining host and guest structures. These cages are singly occupied and off-centered argon atoms distribute on two tetrahedron-shaped split positions with a ratio 3:1. Molecular dynamics (MD) simulations revealed that the crystallographic positional disorder structure is due to mobile argon atoms even at 20 K. The MD potential energy distribution confirmed the diffraction model. It is noted that the unit cell volumes of argon hydrate in the investigated temperature range are virtually identical to N-2 hydrate, which has a similar composition at ambient pressure, indicating a very similar (slightly attractive) host-guest interaction.

While solid polymer electrolytes (SPEs) have great potential for use in future lithium-based batteries, they do, however, not display conductivity at a sufficient level as compared to liquid electrolytes. To reach the needed requirements of lithium batteries it is therefore necessary to explore new materials classes to serve as novel polymer hosts. In this work, SPEs based on the polyketone poly(3,3-dimethylpentane-2,4-dione) were investigated. Polyketones are structurally similar to several polycarbonate and polyester SPE hosts investigated before but have, due to the lack of additional oxygen atoms in the coordinating motif, even more electronwithdrawing carbonyl groups and could therefore display better properties for coordination to the salt cation. In electrolyte compositions comprising 25-40 wt% LiTFSI salt, it was observed that this polyketone indeed conducts lithium ions with a high cation transference number, but that the ionic conductivity is limited by the semi-crystallinity of the polymer matrix. The crystallinity decreases with increasing salt content, and a fully amorphous SPE can be produced at 40 wt% salt, accompanied by an ionic conductivity of 3 x 10(-7) S cm(-1) at 32 degrees C. This opens up for further exploration of polyketone systems for SPE-based batteries.