JH

J.R. Heringa

info

Please Note

11 records found

A Python toolkit for site-resolved diffusion analysis in solid-state molecular dynamics

Journal article (2026) - Anastasia K. Lavrinenko, Theodosios Famprikis, Victor Landgraf, Jouke R. Heringa, Stef Smeets, Victor Azizi, Simone Ciarella, Marnix Wagemaker, Alexandros Vasileiadis
Molecular dynamics (MD) simulations have become essential for understanding diffusion mechanisms in solid-state materials such as ionic conductors, fuel cells, and gas sensors, yet most existing studies and software tools extract only standard metrics, leaving much of the information contained in the trajectories unused. Here we introduce GEMDAT, a user-friendly Python toolkit for site-resolved diffusion analysis of MD simulations of solid-state materials (https://github.com/GEMDAT-repos/GEMDAT). Beyond mean-squared displacements, radial distribution functions, and Arrhenius-based activation energies, GEMDAT provides jump rates, attempt frequencies, site-specific activation energies, rotational diffusion, and percolation. Our tool provides access to vibrational amplitudes, site geometries, and site occupancies—quantities that are also directly comparable to experimental diffraction data. Migration sites can be defined manually or identified automatically from the trajectory. A built-in caching approach, together with rapid visualization capabilities, makes the workflow fast and interactive. We demonstrate GEMDAT on a series of case studies involving crystalline Li- and Na-ion conductors, plastic crystals, amorphous structures, and surface configurations, showing how the code extracts atomic-level structural features and connects them to macroscopic transport properties, thereby guiding the optimization and development of solid-state materials. ...
By varying the bromine content and cooling method, we are able to induce site disorder in the Li6-xPS5-xBr1+x (x = 0, 0.3, 0.5) system via two routes, allowing us to disentangle the impact of site disorder and chemical composition on conductivity. Through solid-state nuclear magnetic resonance (NMR), we can explore the chemical environment as well as short-range lithium-ion dynamics and compare these to results obtained from neutron diffraction and electrochemical impedance spectroscopy (EIS). We find that the cooling method has a profound effect on the 7Li and 31P environment that cannot be explained through 4d site disorder alone. The configurational entropy (Sconf) is used as a more complete descriptor of structural disorder and linked to distortions in both the phosphorus and lithium environment. These distortions are correlated to increased intercage movement through 7Li T1 spin-lattice relaxation (SLR) NMR. Further analysis of the prefactors obtained from SLR NMR and EIS allows us to obtain the migrational entropy (ΔSm). For short-range SLR movement, the ΔSm correlates well with Sconf, implying that increased intercage movement is related to distortion of the lithium cages as well as a decrease of the intercage distance. Comparison to EIS shows that an increase in short-range movement translates into increased long-range movement in a straightforward manner for slow-cooled samples. However, for quench-cooled samples, this correlation is lost. Lattice softness and phonon-ion interactions are suggested to play an important role in long-range conduction which only becomes apparent when chemical composition and disorder are disentangled. This work shows that by altering one synthesis step, the relationship between site-occupancy-based descriptors (site disorder or Sconf) and lithium dynamics is changed profoundly. Furthermore, it shows that chemical composition and descriptors of site disorder cannot be seen as one and the same, as both play a role that changes with the length scale probed. Finally, it challenges the implicit assumption that increased short-range diffusivity automatically results in increased long-range diffusivity. ...
Journal article (2025) - Pascalle Jacobs, Tony Van Der Hiele, Jouke Heringa, Gabrielle P. Verbeeke, Linda Tonk, Remment Ter Hofstede
European oysters (Ostrea edulis) once covered large areas of the North Sea, but have disappeared due to a combination of overexploitation and the destruction of benthic habitats including hard settlement substrate. Offshore wind parks offer an opportunity for oyster restoration as fishing is banned inside these parks and scour protection provides hard settlement substrate. However, ecological restoration of marine systems is capital-intensive. The success of restoration projects is mainly determined by the choice of methods and techniques and consequently costs. Costs and cost-effectiveness information are therefore key in decision making processes concerning the selection of restoration efforts and techniques. So far, economic viability of marine ecosystem restoration have mainly focused on near-shore shallow habitats. The aim of this study was to provide insight into the most cost-effective deployment options to create a European flat oyster reef in an offshore wind farm in the North Sea. Within the current policy and legislation framework, several deployment scenarios were identified based on best practices, expert knowledge, and preliminary results of several pilots. The 9 scenarios included 'adults placed loose on the seafloor, 'adults glued on granite, 'spat settled on shells, 'spat settled on granite and a combined adult and spat scenario. Cost-effectiveness of the different scenarios was determined by modelling the expected reef biomass post-deployment both with and without the option to add additional settlement substrate post-deployment. The main conclusions from this exercise were that: 1. based on investment value, the scenarios adult loose on the seafloor, 'adults in cages and 'spat on shells had the highest revenues per Euro invested; 2. adding substrate in the years post-deployment increased cost-effectiveness in the model for all scenarios, and 3. the time post-deployment to reach a self-sustaining adult oyster population was, with 8-10 yr, shortest for the scenarios 'spat settled on shells' and the combined scenario of 'adults placed loose on the seafloor' and 'spat settled on shells'. ...

A unified view on the role of sulfur/halide distribution and local environments

Understanding diffusion mechanisms in solid electrolytes is crucial for advancing solid-state battery technologies. This study investigates the role of structural disorder in Li7−xPS6−xBrx argyrodites using ab initio molecular dynamics, focusing on the correlation between key structural descriptors and Li-ion conductivity. Commonly suggested parameters, such as configurational entropy, bromide site occupancy, and bromine content, correlate with Li-ion diffusivity but do not consistently explain conductivity trends. We find that a uniform distribution of bromine and sulfur ions across the 4a and 4d sublattices is critical for achieving high conductivity by facilitating optimal lithium jump activation energies, anion-lithium distances, and charge distribution. Additionally, we introduce the ionic potential as a simple descriptor that predicts argyrodite conductivity by assessing the interaction strength between cations and anions. By analyzing the correlation between ionic potential and conductivity for a range of argyrodite compositions published over the past decade, we demonstrate its broad applicability. Minimizing and equalizing ionic potentials across both sublattices enhances conductivity by reducing the strength of anion-lithium interactions. Our analysis of local environments coordinating Li jumps reveals that balancing high and low-energy pathways is crucial for enabling macroscopic diffusion, supported by investigating percolating pathways. This study highlights the significance of the anionic framework in lithium mobility and informs the design of solid electrolytes for improved energy storage systems. ...
Journal article (2024) - Qidi Wang, Jianlin Wang, Jouke R. Heringa, Xuedong Bai, Marnix Wagemaker
One of the primary challenges to improving lithium-ion batteries lies in comprehending and controlling the intricate interphases. However, the complexity of interface reactions and the buried nature make it difficult to establish the relationship between the interphase characteristics and electrolyte chemistry. Herein, we employ diverse characterization techniques to investigate the progression of electrode-electrolyte interphases, bringing forward opportunities to improve the interphase properties by what we refer to as high-entropy solvation disordered electrolytes. Through formulating an electrolyte with a regular 1.0 M concentration that includes multiple commercial lithium salts, the solvation interaction with lithium ions alters fundamentally. The participation of several salts can result in a weaker solvation interaction, giving rise to an anion-rich and disordered solvation sheath despite the low salt concentration. This induces a conformal, inorganic-rich interphase that effectively passivates electrodes, preventing solvent co-intercalation. Remarkably, this electrolyte significantly enhances the performance of graphite-containing anodes paired with high-capacity cathodes, offering a promising avenue for tailoring interphase chemistries. ...
Journal article (2021) - Bei Tian, Jouke R. Heringa, Wim G. Bouwman
Small angle scattering is frequently applied to study the anisotropy in complex soft matter systems. One emerging application is to probe the multi-scale structure in food matrices; while few models are available to describe the anisotropic scattering pattern in a quantitative, yet simple manner. For this purpose, anisotropy is introduced to the Guinier–Porod model to study the scattering from non-spherical objects with a preferred orientation. This generalised anisotropic Guinier–Porod model can be adapted to approximate the sector scattering from both cylinders and ellipsoids (both prolate and oblate). In practice, it is applied to describe the anisotropic scattering from fibres in a meat analogue made of calcium caseinate. A good agreement is found between fitted dimensions of the fibres and those observed from the microscopy image. The effect of orientation distribution on the shape and intensity of the scattering pattern is further discussed and three means to obtain the orientation distribution of the symmetry axis are proposed. Given the model is straightforward and the fitting remains phenomenological, it provides a novel approach to extract information from complex food systems. ...
Phase transitions play a crucial role in Li-ion battery electrodes being decisive for both the power density and cycle life. The kinetic properties of phase transitions are relatively unexplored and the nature of the phase transition in defective spinel Li4+ xTi5O12 introduces a controversy as the very constant (dis)charge potential, associated with a first-order phase transition, appears to contradict the exceptionally high rate performance associated with a solid-solution reaction. With the present density functional theory study, a microscopic mechanism is put forward that provides deeper insight in this intriguing and technologically relevant material. The local substitution of Ti with Li in the spinel Li4+ xTi5O12 lattice stabilizes the phase boundaries that are introduced upon Li-ion insertion. This facilitates a subnanometer phase coexistence in equilibrium, which although very similar to a solid solution should be considered a true first-order phase transition. The resulting interfaces are predicted to be very mobile due to the high mobility of the Li ions located at the interfaces. This highly mobile, almost liquid-like, subnanometer phase morphology is able to respond very fast to nonequilibrium conditions during battery operation, explaining the excellent rate performance in combination with a first-order phase transition. ...
Journal article (2017) - Z. Li, Swapna Ganapathy, Y. Xu, Jouke Heringa, Quanyao Zhu, Wen Chen, Marnix Wagemaker
The lithium air, or Li–O2, battery system is a promising electrochemical energy storage system because of its very high theoretical specific energy, as required by automotive applications. Fundamental research has resulted in much progress in mitigating detrimental (electro)chemical processes; however, the detailed structural evolution of the crystalline Li2O2 and LiOH discharge products, held at least partially responsible for the limited reversibility and poor rate performance, is hard to measure operando under realistic electrochemical conditions. This study uses Rietveld refinement of operando X-ray diffraction data during a complete discharge–charge cycle to reveal the detailed structural evolution of Li2O2 and LiOH crystallites in 1,2-dimethoxyethane (DME) and DME/LiI electrolytes, respectively. The anisotropic broadened reflections confirm and quantify the platelet crystallite shape of Li2O2 and LiOH and show how the average crystallite shape evolves during discharge and charge. Li2O2 is shown to form via a nucleation and growth mechanism, whereas the decomposition appears to start at the smallest Li2O2 crystallite sizes because of their larger exposed surface. In the presence of LiI, platelet LiOH crystallites are formed by a particle-by-particle nucleation and growth process, and at the end of discharge, H2O depletion is suggested to result in substoichiometric Li(OH)1–x, which appears to be preferentially decomposed during charging. Operando X-ray diffraction proves the cyclic formation and decomposition of the LiOH crystallites in the presence of LiI over multiple cycles, and the structural evolution provides key information for understanding and improving these highly relevant electrochemical systems. ...
Journal article (2016) - L.J. Bannenberg, H. Schreuders, L. van Eijck, J.R. Heringa, N.J. Steinke, R.M. Dalgliesh, B. Dam, F.M. Mulder, A.A. van Well
Nanostructuring is widely applied in both battery and hydrogen materials to improve the performance of these materials as energy carriers. Nanostructuring changes the diffusion length as well as the thermodynamics of materials. We studied the impact of nanostructuring on the hydrogenation in a model system consisting of a thin film of magnesium sandwiched between two titanium layers and capped with palladium. While we verified optically the coexistence of the metallic α-MgDx and the insulating β-MgD2–y phase, neutron reflectometry shows significant deviations from the thermodynamic solubility limits in bulk magnesium during the phase transformation. This suggests that the kinetics of the phase transformation in nanostructured battery and hydrogen storage systems is enhanced not only as a result of the reduced length scale but also due to the increased solubility in the parent phases. ...