Publications

@article{Birschitzky_2024b,

title = {Machine Learning Small Polaron Dynamics},

author = {Viktor C. Birschitzky and Luca Leoni and Michele Reticcioli and Cesare Franchini},

url = {https://doi.org/10.1103/PhysRevLett.134.216301},

year  = {2025},

date = {2025-05-27},

urldate = {2024-09-24},

journal = {Physical Review Letters},

volume = {134},

issue = {21},

pages = {216301},

abstract = {Polarons are crucial for charge transport in semiconductors, significantly impacting material properties and device performance. The dynamics of small polarons can be investigated using first-principles molecular dynamics. However, the limited timescale of these simulations presents a challenge for adequately sampling infrequent polaron hopping events. Here, we introduce a message-passing neural network combined with first-principles molecular dynamics within the Born-Oppenheimer approximation that learns the polaronic potential energy surface by encoding the polaronic state, allowing for simulations of polaron hopping dynamics at the nanosecond scale. By leveraging the statistical significance of the long timescale, our framework can accurately estimate polaron (anisotropic) mobilities and activation barriers in prototypical polaronic oxides across different scenarios (hole polarons in rocksalt MgO and electron polarons in pristine and F-doped rutile TiO_{2}) within experimentally measured ranges.},

keywords = {P07},

pubstate = {published},

tppubtype = {article}

}

@article{Sokolovic_2025a,

title = {Duality and degeneracy lifting in two-dimensional electron liquids on SrTiO_{3}(001)},

author = {Igor Sokolović and Eduardo B. Guedes and Thomas P. van Waas and Fei Guo and Samuel Poncé and Craig Polley and Michael Schmid and Ulrike Diebold and Milan Radović and Martin Setvín and J. Hugo Dil},

url = {https://doi.org/10.1038/s41467-025-59258-4},

year  = {2025},

date = {2025-05-17},

journal = {Nature Communications},

volume = {16},

issue = {1},

pages = {4594},

abstract = {Two-dimensional electron liquids (2DELs) have increasing technological relevance for ultrafast electronics and spintronics, yet significant gaps in their fundamental understanding are exemplified on the prototypical SrTiO_{3}. We correlate the exact SrTiO_{3}(001) surface structure with distinct 2DELs through combined microscopic angle-resolved photoemission spectroscopy and non-contact atomic force microscopy on truly bulk-terminated surfaces that alleviate structural uncertainties inherent to this long-studied system. The SrO termination is shown to develop a 2DEL following the creation of oxygen vacancies, unlike the intrinsically metallic TiO_{2} termination. Differences in degeneracy of the 2DELs, with nearly the same band filling and identical band bending, are assigned to polar distortions of the Ti atoms in combination with spin order, supported with the extraction of fundamental electron-phonon coupling strength. These results not only resolve the ambiguities regarding 2DELs on SrTiO_{3} thus far, but also pave the way to manipulating band filling and spin order in oxide 2DELs in general.},

keywords = {P02},

pubstate = {published},

tppubtype = {article}

}

@article{Coretti2022,

title = {Learning Mappings between Equilibrium States of Liquid Systems Using Normalizing Flows},

author = {Alessandro Coretti and Sebastian Falkner and Phillip Geissler and Christoph Dellago},

doi = {10.1063/5.0253034},

year  = {2025},

date = {2025-05-08},

urldate = {2022-08-22},

journal = {The Journal of Chemical Physics},

volume = {162},

issue = {18},

pages = {184102},

abstract = {Generative models and, in particular, normalizing flows are a promising tool in statistical mechanics to address the sampling problem in condensed-matter systems. In this work, we investigate the potential of normalizing flows to learn a transformation to map different liquid systems into each other while allowing at the same time to obtain an unbiased equilibrium distribution. We apply this methodology to the mapping of a small system of fully repulsive disks modeled via the Weeks–Chandler–Andersen potential into a Lennard-Jones system in the liquid phase at different coordinates in the phase diagram. We obtain an improvement in the relative effective sample size of the generated distribution up to a factor of six compared to direct reweighting. We show that this factor can have a strong dependency on the thermodynamic parameters of the source and target system.},

keywords = {P12},

pubstate = {published},

tppubtype = {article}

}

@article{Cao_2025a,

title = {Quantum Delocalization Enables Water Dissociation on Ru(0001)},

author = {Yu Cao and Jiantao Wang and Mingfeng Liu and Yan Liu and Hui Ma and Cesare Franchini and Yan Sun and Georg Kresse and Xing-Qiu Chen and Peitao Liu},

doi = {10.1103/PhysRevLett.134.178001},

year  = {2025},

date = {2025-04-30},

journal = {Physical Review Letters},

volume = {134},

issue = {17},

pages = {178001},

abstract = {We revisit the long-standing question of whether water molecules dissociate on the Ru(0001) surface through nanosecond-scale path-integral molecular dynamics simulations on a sizable supercell. This is made possible through the development of an efficient and reliable machine-learning potential with near first-principles accuracy, overcoming the limitations of previous ab initio studies. We show that the quantum delocalization associated with nuclear quantum effects enables rapid and frequent proton transfers between water molecules, thereby facilitating the water dissociation on Ru(0001). This work provides the direct theoretical evidence of water dissociation on Ru(0001), resolving the enduring issue in surface sciences and offering crucial atomistic insights into water-metal interfaces.},

keywords = {P03, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Yigit_2024a,

title = {Preferential CO oxidation (PROX) on LaCoO_{3}–based catalysts: Effect of cobalt oxidation state on selectivity},

author = {Nevzat Yigit and Karin Föttinger and Johannes Bernardi and Günther Rupprechter},

url = {https://doi.org/10.1016/j.jcat.2025.115973},

year  = {2025},

date = {2025-01-20},

urldate = {2025-01-20},

journal = {Journal of Catalysis},

issue = {0021-9517},

pages = {115973},

abstract = {The perovskite LaCoO_{3} (LCO) was used as catalyst for preferential oxidation of CO (PROX). LCO was synthesized via the modified Pechini method and characterized by X-ray diffraction (XRD), high resolution transmission electron microscopy (HRTEM), and CO_{2–} and H_{2–} temperature programmed reduction (TPR). Different reductive and oxidative pretreatments were applied to systematically vary the Co oxidation state in order to examine its effect on catalytic performance and to single out active site requirements. Upon reduction at increasing temperature, LaCoO_{3} transformed to brownmillerite-type La_{2}Co_{2}O_{5}, exsolved Co^{0} nanoparticles supported on La_{2}O_{3} and, upon reoxidation, to Co_{3}O_{4}/La_{2}O_{3}. The Co oxidation state of the various catalysts correlated with their CO_{2} selectivity: LCO containing only Co^{3+} exhibited 100 % CO_{2} selectivity in a wide temperature window, whereas La_{2}Co_{2}O_{5}, Co/La_{2}O_{3} and Co_{3}O_{4}/La_{2}O_{3} had markedly lower selectivity. It is suggested that Co^{3+} is crucial and that the strong resistivity of LaCoO_{3} towards reduction is responsible for the high and stable CO_{2} selectivity over a temperature range of 100 °C–220 °C. Higher oxygen concentration further broadens the PROX window.},

keywords = {P08, P10},

pubstate = {published},

tppubtype = {article}

}

@article{Li_2024a,

title = {Sum frequency generation (SFG) spectroscopy at surfaces and interfaces: Adsorbate structure and molecular bond orientation},

author = {Xia Li and Günther Rupprechter},

url = {https://doi.org/10.1016/j.surfrep.2024.100645},

year  = {2024},

date = {2024-12-05},

urldate = {2024-12-05},

journal = {Surface Science Reports},

volume = {79},

issue = {4},

pages = {100645},

abstract = {Infrared (IR)-visible (Vis) sum frequency generation (SFG) is a second-order nonlinear optical process which is forbidden in centrosymmetric bulk media or isotropic phases, but allowed at (open) surfaces or (buried) interfaces where the inversion symmetry is broken. SFG spectroscopy is thus inherently surface- or interface-specific, providing information about the structure, orientation, surface number density, chirality, and dynamics of molecules, provided the system of interest is accessible by light. This review illustrates basic SFG concepts, theory, operation modes (e.g., frequency-domain, broadband, homodyne/heterodyne, time-resolved), and recent extensions and developments of SFG (e.g., doubly resonant, plasmon-enhanced, chiral, microscopy). To illustrate the wide range of SFG applications, selected case studies discuss the characterization of molecular structure and bond orientation at solid-gas, solid-liquid, liquid-air, liquid-liquid, and solid-solid interfaces.},

keywords = {P08},

pubstate = {published},

tppubtype = {article}

}

@article{Tangpakonsab_2024a,

title = {Mechanistic Insights into CO and H_{2} Oxidation on Cu/CeO_{2} Single Atom Catalysts: A Computational Investigation},

author = {Parinya Tangpakonsab and Alexander Genest and Gareth S. Parkinson and Günther Rupprechter},

url = {https://chemrxiv.org/engage/chemrxiv/article-details/673325d25a82cea2fadefd76},

year  = {2024},

date = {2024-11-14},

urldate = {2024-11-14},

journal = {ChemRxiv},

abstract = {Single atom catalysts (SACs) have attracted significant interest due to their unique properties and potential for enhancing catalytic performance in various chemical reactions. In this study, we atomistically explore adsorption properties and catalytic performance of single Cu atoms anchored at low-index CeO_{2} surfaces, focusing on the oxidation of CO and H_{2}. Utilizing density functional theory (DFT) calculations, we report that Cu adatoms bind favorably on different CeO_{2} surfaces, following a stability order of (100)>(110)>(111). The charge transfer from a single adsorbed Cu atom to Ce leads to the reduction of Ce^{4+} to Ce^{3+} and the oxidation of Cu^{0} to Cu^{+}. This strengthens molecular bonds at Cu sites, particularly for CO due to the less populated d-band, while H_{2} shows a by ~1 eV weaker adsorption. CO oxidation is energetically more favorable than H_{2} oxidation on the Cu/CeO_{2}(111) surface. The rate-controlling steps for the Mars-van Krevelen oxidation involve the formation of a bent CO_{2}^{-} intermediate for CO and H_{2}O for H_{2}. The lattice oxygen atom at the interface plays a key role for both oxidation processes. Our findings highlight the potential of single atom catalyst, Cu/CeO_{2}, for CO adsorption and oxidation in heterogeneous catalysis.},

keywords = {P04, P08},

pubstate = {published},

tppubtype = {article}

}

@article{Bichelmaier_2024a,

title = {Neural network enabled molecular dynamics study of HfO_{2} phase transitions},

author = {Sebastian Bichelmaier and Jesús Carrete and Georg K. H. Madsen},

url = {https://arxiv.org/abs/2408.02429},

doi = {10.1103/PhysRevB.110.174105},

year  = {2024},

date = {2024-11-07},

journal = {Phys. Rev. B},

volume = {110},

pages = {174105},

abstract = {The advances of machine-learned force fields have opened up molecular dynamics (MD) simulations for compounds for which ab initio MD is too resource intensive and phenomena for which classical force fields are insufficient. Here we describe a neural-network force field parametrized to reproduce the r2⁢SCAN potential energy landscape of HfO_{2}. Based on an automatic differentiable implementation of the isothermal-isobaric (𝑁⁢𝑃⁢𝑇) ensemble with flexible cell fluctuations, we study the phase space of HfO_{2}. We find excellent predictive capabilities regarding the lattice constants and experimental x-ray diffraction data. The phase transition away from monoclinic is clearly visible at a temperature around 2000 K, in agreement with available experimental data and previous calculations. Another abrupt change in lattice constants occurs around 3000 K. While the resulting lattice constants are closer to cubic, they exhibit a small tetragonal distortion, and there is no associated change in volume. We show that this high-temperature structure is in agreement with the available high-temperature diffraction data.},

keywords = {P09},

pubstate = {published},

tppubtype = {article}

}

@article{Redondo2024,

title = {Real-space investigation of polarons in hematite Fe_{2}O_{3}},

author = {Jesus Redondo and Michele Reticcioli and Vit Gabriel and Dominik Wrana and Florian Ellinger and Michele Riva and Giada Franceschi and Erik Rheinfrank and Igor Sokolovic and Zdenek Jakub and Florian Kraushofer and Aji Alexander and Laerte L. Patera and Jascha Repp and Michael Schmid and Ulrike Diebold and Gareth S. Parkinson and Cesare Franchini and Pavel Kocan and Martin Setvin},

url = {https://arxiv.org/abs/2303.17945

https://www.science.org/doi/10.1126/sciadv.adp7833},

year  = {2024},

date = {2024-11-01},

urldate = {2024-09-27},

journal = {Science Advances},

volume = {10},

issue = {44},

pages = {eadp7833},

abstract = {In polarizable materials, electronic charge carriers interact with the surrounding ions, leading to quasiparticle behaviour. The resulting polarons play a central role in many materials properties including electrical transport, optical properties, surface reactivity and magnetoresistance, and polaron properties are typically investigated indirectly through such macroscopic characteristics. Here, noncontact atomic force microscopy (nc-AFM) is used to directly image polarons in Fe_{2}O_{3} at the single quasiparticle limit. A combination of Kelvin probe force microscopy (KPFM) and kinetic Monte Carlo (KMC) simulations shows that Ti doping dramatically enhances the mobility of electron polarons, and density functional theory (DFT) calculations indicate that a metallic transition state is responsible for the enhancement. In contrast, hole polarons are significantly less mobile and their hopping is hampered further by the introduction of trapping centres.},

keywords = {P02, P04, P07},

pubstate = {published},

tppubtype = {article}

}

@article{Romano_2024b,

title = {Structure of the water/magnetite interface from sum frequency generation experiments and neural network based molecular dynamics simulations},

author = {Salvatore Romano and Harsharan Kaur and Moritz Zelenka and Pablo Montero De Hijes and Moritz Eder and Gareth S. Parkinson and Ellen H. G. Backus and Christoph Dellago},

url = {https://arxiv.org/abs/2410.12717},

year  = {2024},

date = {2024-10-16},

urldate = {2024-10-16},

journal = {arXiv},

abstract = {Magnetite, a naturally abundant mineral, frequently interacts with water in both natural settings and various technical applications, making the study of its surface chemistry highly relevant. In this work, we investigate the hydrogen bonding dynamics and the presence of hydroxyl species at the magnetite-water interface using a combination of neural network potential-based molecular dynamics simulations and sum frequency generation vibrational spectroscopy. Our simulations, which involved large water systems, allowed us to identify distinct interfacial species, such as dissociated hydrogen and hydroxide ions formed by water dissociation. Notably, water molecules near the interface exhibited a preference for dipole orientation towards the surface, with bulk-like water behavior only re-emerging beyond 60 Å from the surface. The vibrational spectroscopy results aligned well with the simulations, confirming the presence of a hydrogen bond network in the surface ad-layers. The analysis revealed that surface-adsorbed hydroxyl groups orient their hydrogen atoms towards the water bulk. In contrast, hydrogen-bonded water molecules align with their hydrogen atoms pointing towards the magnetite surface.},

keywords = {P04, P11, P12},

pubstate = {published},

tppubtype = {article}

}