Researcher profile

Claudia Draxl

Claudia Draxl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2026arXiv

From Knowledge to Action: Outcomes of the 2025 Large Language Model (LLM) Hackathon for Applications in Materials Science and Chemistry

Large language models (LLMs) are rapidly changing how researchers in materials science and chemistry discover, organize, and act on scientific knowledge. This paper analyzes a broad set of community-developed LLM applications in an effort to identify emerging patterns in how these systems can be used across the scientific research lifecycle. We organize the projects into two complementary categories: Knowledge Infrastructure, systems that structure, retrieve, synthesize, and validate scientific information; and Action Systems, systems that execute, coordinate, or automate scientific work across computational and experimental environments. The submissions reveal a shift from single-purpose LLM tools toward integrated, multi-agent workflows that combine retrieval, reasoning, tool use, and domain-specific validation. Prominent themes include retrieval-augmented generation as grounding infrastructure, persistent structured knowledge representations, multimodal and multilingual scientific inputs, and early progress toward laboratory-integrated closed-loop systems. Together, these results suggest that LLMs are evolving from general-purpose assistants into composable infrastructure for scientific reasoning and action. This work provides a community snapshot of that transition and a practical taxonomy for understanding emerging LLM-enabled workflows in materials science and chemistry.

preprint2022arXiv

A consistent picture of excitations in cubic BaSnO$_{3}$ revealed by combining theory and experiment

Among the transparent conducting oxides, the perovskite barium stannate is most promising for various electronic applications due to its outstanding carrier mobility achieved at room temperature. However, most of its important characteristics, such as band gaps, effective masses, and absorption edge, remain controversial. Here, we provide a fully consistent picture by combining state-of-the-art {\it ab initio} methodology with forefront electron energy-loss spectroscopy and optical absorption measurements. Valence electron energy-loss spectra, featuring signals originating from band gap transitions, are acquired on defect-free sample regions of a BaSnO$_{3}$ single crystal. These high-energy-resolution measurements are able to capture also very weak excitations below the optical gap, attributed to indirect transitions. By temperature-dependent optical absorption measurements, we assess band-gap renormalization effects induced by electron-phonon coupling. Overall, we find for the effective electronic mass, the direct and the indirect gap, the optical gap, as well as the absorption onsets and spectra, excellent agreement between both experimental techniques and the theoretical many-body results, supporting also the picture of a phonon-mediated mechanism where indirect transitions are activated by phonon-induced symmetry lowering. This work demonstrates a fruitful connection between different high-level theoretical and experimental methods for exploring the characteristics of advanced materials.

preprint2022arXiv

All-electron many-body approach to resonant inelastic x-ray scattering

We present a formalism for the resonant inelastic x-ray scattering (RIXS) cross section. The resulting compact expression in terms of polarizability matrix elements, particularly lends itself to the implementation in an all-electron many-body perturbation theory (MBPT) framework, which is realized in the full-potential package exciting. With the carbon K edge RIXS of diamond and the oxygen K edge RIXS of $β-\mathrm{Ga_2 O_3}$, respectively, we demonstrate the importance of electron-hole correlation and atomic coherence in the RIXS spectra.

preprint2022arXiv

Density-of-states similarity descriptor for unsupervised learning from materials data

We develop a materials descriptor based on the electronic density of states and investigate the similarity of materials based on it. As an application example, we study the Computational 2D Materials Database that hosts thousands of two-dimensional materials with their properties calculated by density-functional theory. Combining our descriptor with a clustering algorithm, we identify groups of materials with similar electronic structure. We characterize these clusters in terms of their crystal structure, their atomic composition, and the respective electronic configurations to rationalize the found (dis)similarities.

preprint2022arXiv

Ehrenfest dynamics implemented in the all-electron package exciting

Ehrenfest Dynamics combined with real-time time-dependent density functional theory has proven to be a reliable tool to study non-adiabatic molecular dynamics with a reasonable computational cost. Among other possibilities, it allows for assessing in real time electronic excitations generated by ultra-fast laser pulses, as e.g., in pump-probe spectroscopy, and their coupling to the nuclear vibrations even beyond the linear regime. In this work, we present its implementation in the all-electron full-potential package exciting. Three cases are presented as examples: diamond and cubic BN relaxed after an initial lattice distortion, and cubic BN exposed to a laser pulse. Comparison with the Octopus code exhibits good agreement.

preprint2022arXiv

FAIR data enabling new horizons for materials research

The prosperity and lifestyle of our society are very much governed by achievements in condensed matter physics, chemistry and materials science, because new products for sectors such as energy, the environment, health, mobility and information technology (IT) rely largely on improved or even new materials. Examples include solid-state lighting, touchscreens, batteries, implants, drug delivery and many more. The enormous amount of research data produced every day in these fields represents a gold mine of the twenty-first century. This gold mine is, however, of little value if these data are not comprehensively characterized and made available. How can we refine this feedstock; that is, turn data into knowledge and value? For this, a FAIR (findable, accessible, interoperable and reusable) data infrastructure is a must. Only then can data be readily shared and explored using data analytics and artificial intelligence (AI) methods. Making data 'findable and AI ready' (a forward-looking interpretation of the acronym) will change the way in which science is carried out today. In this Perspective, we discuss how we can prepare to make this happen for the field of materials science.

preprint2022arXiv

How ferroelectric BaTiO$_3$ can tune a two-dimensional electron gas at the interface of LaInO3 and BaSnO$_3$: a first-principles study

The emerging interest in two-dimensional electron gases (2DEGs), formed at interfaces between two insulating oxide perovskites poses crucial fundamental question in view of future electronic devices. In the framework of density-functional theory, we investigate the possibility to control the characteristics of the 2DEG formed at the LaInO$_3$/BaSnO$_3$ interface by including a ferroelectric BaTiO$_3$ layer. To do so, we examine how the orientation of the ferroelectric polarization impacts density and confinement of the 2DEG. We find that aligning the ferroelectric polarization toward (outward) the LaInO$_3$/BaSnO$_3$ interface leads to an accumulation (depletion) of the interfacial 2DEG. Varying its magnitude, we find a linear effect on the 2DEG charge density that is confined within the BaSnO$_3$ side. Analysis of the optimized geometries revels that inclusion of the BaTiO$_3$ block makes structural distortions at the LaInO$_3$/BaSnO$_3$ less pronounced, which, in turn, enhances the 2DEG density. Thicker ferroelectric layers allow for reaching higher polarization magnitudes. We discuss the mechanisms behind all these findings and rationalize how the characteristics of both 2DEG and 2D hole gases can be controlled in the considered heterostructures. Overall, our results can be generalized to other combinations of ferroelectric, polar, and nonpolar materials.

preprint2022arXiv

Influence of spin-orbit coupling on chemical bonding

The influence of spin-orbit interaction on chemical bonds in elemental solids and homonuclear dimers is analyzed by means of density-functional-theory calculations. Employing highly precise all-electron full-potential methodology, our results represent benchmark quality. Comparison of the scalar- and fully-relativistic approaches for elemental solids shows that the spin-orbit interaction may contract or expand the volume of the considered material. The largest variation of the volume is obtained for Au, Tl, I, Bi, Po and Hg, exhibiting changes between 1.0--7.6\%. Using the tight-binding model, we show for diatomic molecules that the nature of this effect lies in the angular rearrangement of bonding and antibonding orbitals introduced by spin-orbit coupling. Such an angular rearrangement appears in partially filled $p$- or $d$-orbitals in heavy elements. Finally, we discuss the impact of the relativistic effects on the chemical bonding in single-layer iodides and transition metal dichalcogenides

preprint2022arXiv

Numerical Quality Control for DFT-based Materials Databases

Electronic-structure theory is a strong pillar of materials science. Many different computer codes that employ different approaches are used by the community to solve various scientific problems. Still, the precision of different packages has only recently been scrutinized thoroughly, focusing on a specific task, namely selecting a popular density functional, and using unusually high, extremely precise numerical settings for investigating 71 monoatomic crystals. Little is known, however, about method- and code-specific uncertainties that arise under numerical settings that are commonly used in practice. We shed light on this issue by investigating the deviations in total and relative energies as a function of computational parameters. Using typical settings for basis sets and k-grids, we compare results for 71 elemental and 63 binary solids obtained by three different electronic-structure codes that employ fundamentally different strategies. On the basis of the observed trends, we propose a simple, analytical model for the estimation of the errors associated with the basis-set incompleteness. We cross-validate this model using ternary systems obtained from the NOMAD Repository and discuss how our approach enables the comparison of the heterogeneous data present in computational materials databases.

preprint2022arXiv

Robust excitons across the phase transition of two-dimensional hybrid perovskites

Two-dimensional halide perovskites are among intensely studied materials platforms profiting from solution based growth and chemical flexibility. They feature exceptionally strong interactions among electronic, optical as well as vibrational excitations and hold a great potential for future optoelectronic applications. A key feature for these materials is the occurrence of structural phase transitions that can impact their functional properties, including the electronic band gap and optical response dominated by excitons. However, to what extent the phase-transitions in two-dimensional perovskites alter the fundamental exciton properties remains barely explored so far. Here, we study the influence of the phase transition on both exciton binding energy and exciton diffusion, demonstrating their robust nature across the phase transition. These findings are unexpected in view of the associated substantial changes of the free carrier masses, strongly contrast broadly considered effective mass and drift-diffusion transport mechanisms, highlighting the unusual nature of excitons in two-dimensional perovskites.

preprint2022arXiv

Tuning two-dimensional electron (hole) gases at LaInO$_{3}$/BaSnO$_{3}$ interfaces: Impact of polar distortions, termination, and thickness

Two-dimensional election gases (2DEG), arising due to quantum confinement at interfaces between transparent conducting oxides, have received tremendous attention in view of electronic applications. The challenge is to find a material system that exhibits both a high charge-carrier density and mobility, at and even above room temperature. Here, we explore the potential of interfaces formed by two lattice-matched wide-gap oxides of emerging interest, $\textit{i.e.}$, the polar, orthorhombic perovskite LaInO$_{3}$ and the non-polar, cubic perovskite BaSnO$_{3}$, employing density-functional theory and many-body theory. We demonstrate that this material combination exhibits all the key features for reaching the goal. For periodic heterostructures, we find that the polar discontinuity at the interface is mainly compensated by electronic relaxation through charge transfer from the LaInO$_{3}$ to the BaSnO$_{3}$ side. This leads to the formation of a 2DEG hosted by the highly-dispersive Sn-$s$-derived conduction band and a 2D hole gas of O-$p$ character, strongly localized inside LaInO$_{3}$. Remarkably, structural distortions through octahedra tilts induce a depolarization field counteracting the polar discontinuity, and thus increasing the $critical$ (minimal) LaInO$_{3}$ thickness, $t_c$, required for the formation of a 2DEG. These polar distortions decrease with increasing LaInO$_{3}$ thickness, enhancing the polar discontinuity and leading to a 2DEG density of 0.5 electron per unit-cell surface. Interestingly, in non-periodic heterostructures, these distortions lead to a decrease of $t_c$, thereby enhancing and delocalizing the 2DEG. We rationalize how polar distortions, termination, and thickness can be exploited in view of tailoring the 2DEG characteristics, and why this material is superior to the most studied prototype LaAlO$_{3}$/SrTiO$_{3}$.

preprint2022arXiv

Two-dimensional electron gas at LaInO$_3$/BaSnO$_3$ interfaces controlled by a ferroelectric layer

With the example of LaInO$_{3}$/BaSnO$_3$, we demonstrate how both density and distribution of a two-dimensional electron gas (2DEG) formed at the interface between these perovskite oxides, can be efficiently controlled by a ferroelectric functional material. A polarization induced in a BaTiO$_3$ layer pointing toward the interface enhances the polar discontinuity which, in turn, significantly increases the 2DEG density and confinement, while, the opposite polarization depletes the 2DEG. Our predictions and analysis, based on first-principles calculations, can serve as a guide for designing such material combinations to be used in electronic devices.

preprint2021arXiv

Elastic stability of Ga$_2$O$_3$: Addressing the $β$ to $α$ phase transition from first principles

Elastic and structural properties of $β$-Ga$_2$O$_3$ and $α$-Ga$_2$O$_3$ are investigated from first principles. The full elastic tensors and elastic moduli of both phases at $0$ K are computed in the framework of semi-local density-functional theory. We determine mechanical instabilities of $β$-Ga$_2$O$_3$ by evaluating the full stiffness tensor under load for a range of hydrostatic pressure values. While a phase transition from the $β$ to $α$ phase is found to be energetically favored at $2.6$ GPa, we show that the $β$ phase is only mechanically unstable for much higher pressures ($>30$ GPa), which agrees well with experimental results. Our employed approach is based on the Born stability criterion, is independent of crystal symmetry, and thus can be readily applied to different materials.

preprint2021arXiv

Fingerprints of optical absorption in the perovskite LaInO$_{3}$: Insight from many-body theory and experiment

We provide a combined theoretical and experimental study of the electronic structure and the optical absorption edge of the orthorhombic perovskite LaInO$_{3}$. Employing density-functional theory and many-body perturbation theory, we predict a direct electronic quasiparticle band gap of about 5 eV and an effective electron (hole) mass of 0.31 (0.48) m$_{0}$. We find the lowest-energy excitation at 0.2 eV below the fundamental gap, reflecting a sizeable electron-hole attraction. Since the transition from the valence band maximum (VBM, $Γ$ point) is, however, dipole forbidden the onset is characterized by weak excitations from transitions around it. The first intense excitation appears about 0.32 eV above. Interestingly, this value coincides with an experimental value obtained by ellipsometry (4.80 eV) which is higher than the onset from optical absorption spectroscopy (4.35 eV). The latter discrepancy is attributed to the fact that the weak transitions that define the optical gap are not resolved by the ellipsometry measurement. The absorption edge shows a strong dependency on the light polarization, reflecting the character of the involved valence states. Temperature-dependent measurements show a redshift of the optical gap by about 120 meV by increasing the temperature from 5 to 300 K. Renormalization due to zero-point vibrations is extrapolated from the latter measurement to amount to 150 meV. By adding the excitonic binding energy of 0.2 eV obtained theoretically to the experimental optical absorption onset, we determine the fundamental band gap at room temperature to be 4.55 eV.

preprint2021arXiv

Rashba and Dresselhaus effects in 2D Pb-I-based perovskites

Bulk hybride halide perovskites are governed by significant Rashba and Dresselhaus splitting. This indicates that such effects will not only affect their optoelectronic properties but also those of their two dimensional layered relatives. This work aims at understanding how different ways of symmetry breaking influence these effects in those materials. For this purpose, model structures are adopted where the organic compounds are replaced by Cs atoms. Disregarding possible distortions in the inorganic layers, results in structures with composition Cs$_{n+1}$Pb$_n$I$_{3n+1}$. Using the all-electron full-potential density-functional-theory code \texttt{exciting}, the impact of atomic displacement on the band structure is systematically studied for $n=1$, 2, 3 and $\infty$. The displacement patterns that yield Rashba or Dresselhaus splitting are identified, and the amount of the splitting is determined as a function of displacement. Furthermore, the spin textures in the electronic states around the band gap are analyzed to differentiate between Rashba and Dresselhaus effects. This study reveals in-plane Pb displacements as the origin of the strongest effects.

preprint2020arXiv

Assessment of approaches for dispersive forces employing graphone as a case study

We have studied two interchange layer systems, (i) free standing partly hydrogenated graphene (graphone), and (ii) graphone on the Nickel (111) surface, to assess various density functional theory based computational schemes incorporating van der Waals forces. The various van der Waals methods have been employed ranging from the semiempirical force-field-like correction of Grimme, through non-local van der Waals density functionals, up to the functionals involving exact exchange and the random phase approximation for the correlation. Generally, all computational schemes lead to a similar qualitative picture of hydrogen layer physisorption and chemisorption to graphene. The largest discrepancies between the approaches emerge for the energetics of the investigated systems. Our studies shed light on the physical mechanisms of graphene hydrogenation both in vacuum and in the proximity of metallic surface. In particular, it is revealed that the adsorption of hydrogen atoms affects the nature of the bonding between graphene and the Ni(111) surface, from the weak to strong semi-covalent bonding. On the other hand, it turns out that the adsorption of hydrogen layer to graphene is stronger in the presence of the metallic surface.

preprint2020arXiv

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the $GW$-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperatures. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

preprint2020arXiv

High Performance Solution of Skew-symmetric Eigenvalue Problems with Applications in Solving the Bethe-Salpeter Eigenvalue Problem

We present a high-performance solver for dense skew-symmetric matrix eigenvalue problems. Our work is motivated by applications in computational quantum physics, where one solution approach to solve the so-called Bethe-Salpeter equation involves the solution of a large, dense, skew-symmetric eigenvalue problem. The computed eigenpairs can be used to compute the optical absorption spectrum of molecules and crystalline systems. One state-of-the art high-performance solver package for symmetric matrices is the ELPA (Eigenvalue SoLvers for Petascale Applications) library. We extend the methods available in ELPA to skew-symmetric matrices. This way, the presented solution method can benefit from the optimizations available in ELPA that make it a well-established, efficient and scalable library, such as GPU support. We compare performance and scalability of our method to the only available high-performance approach for skew-symmetric matrices, an indirect route involving complex arithmetic. In total, we achieve a performance that is up to 3.67 higher than the reference method using Intel's ScaLAPACK implementation. The runtime to solve the Bethe-Salpeter-Eigenvalue problem can be improved by a factor of 10. Our method is freely available in the current release of the ELPA library.

preprint2020arXiv

MoTe2 as a natural hyperbolic material across the visible and the ultraviolet region

Hyperbolic materials are of particular interest for the next generation of photonic and optoelectronic devices. Since artificial metamaterials are intrinsically limited by the size of their nanostructured components, there has been a hunt for natural hyperbolic materials in the last few years. In a first-principles work based on density-functional theory and many-body perturbation theory, we investigate the fundamental dielectric response of MoTe2 in monolayer, bilayer, and bulk form, and find that it is a natural type-II hyperbolic material with low losses between 3 and 6 eV. Going from the monolayer to the bulk, the energy window of hyperbolic dispersion is blue-shifted by a few tenths of an eV. We show that excitonic effects and optical anisotropy play a major role in the hyperbolic behavior of MoTe2. Our results confirm the potential of layered materials as hyperbolic media for opto-electronics, photonics, and nano-imaging applications.

preprint2020arXiv

Robust mixing in self-consistent linearized augmented planewave calculations

We devise a mixing algorithm for full-potential (FP) all-electron calculations in the linearized augmented planewave (LAPW) method. Pulay's direct inversion in the iterative subspace is complemented with the Kerker preconditioner and further improvements to achieve smooth convergence, avoiding charge sloshing and noise in the exchange-correlation potential. As the Kerker preconditioner was originally designed for the planewave basis, we have adapted it to the FP-LAPW method and implemented in the exciting code. Applications to the $2\times 2$ Au(111) surface with a vacancy and to the Pd(111) surface demonstrate that this approach and our implementation work reliably with both density and potential mixing.

preprint2020arXiv

Structural, electronic, and optical properties of periodic graphene/h-BN van der Waals heterostructures

The emerging interest in van der Waals heterostructures as new materials for opto-electronics and photonics poses questions about their stability and structure-property relations. In the framework of density-functional and many-body perturbation theory, we investigate the structural, electronic, and optical properties of periodic heterostructures formed by graphene and hexagonal boron nitride (h-BN). To understand how the constituents affect each other depending on the layer stacking, we examine 12 commensurate arrangements. We find that interaction with graphene improves the stability of bulk h-BN also in those configurations that are predicted to be energetically metastable. In return, the interaction with h-BN can open a band gap of a few hundred meV in graphene. Its actual size can be tuned by the arrangement of the layers. In the semiconducting configurations, the character and spatial distribution of optical excitations are affected by the specific stacking, that determines the electronic states involved in the transitions. Remarkably, six out of the 12 explored heterostructures remain semi-metallic.

preprint2020arXiv

The role of the interface in controlling the epitaxial relationship between orthorhombic $\text{LaInO}_\text{3}$ and cubic $\text{BaSnO}_\text{3}$

Epitaxial perovskite oxide interfaces with different symmetry of the epitaxial layers have attracted considerable attention due to the emergence of novel behaviors and phenomena. In this paper, we show by aberration corrected transmission electron microscopy that orthorhombic $\text{LaInO}_\text{3}$ films grow in form of three different types of domains on the cubic $\text{BaSnO}_\text{3}$ pseudosubstrate. Quantitative evaluation of our TEM data shows that $c_{pc}$-oriented and $a_{pc}/b_{pc}$-oriented domains are present with similar probability. While continuum elasticity theory suggests that $c_{pc}$-oriented domains should exhibit a significantly higher strain energy density than $a_{pc}/b_{pc}$-oriented domains, density functional calculations confirm that $c_{pc}$- and $a_{pc}$-oriented domains on $\text{BaSnO}_\text{3}$ have similar energies.

preprint2020arXiv

Work-function modification of PEG(thiol) adsorbed on the Au(111) surface: A first-principles study

The possibility of modifying the work function of electrodes is important for optimizing the energy barriers for charge-injection (extraction) at the interface to an organic material. In this study, we perform density-functional-theory calculations to investigate the impact of dithiol-terminated polyethylene glycol (PEG(thiol)) based self-assembled monolayers (SAMs) with different numbers of PEG repeat units on the work function of the Au(111) surface. We find that a monolayer of PEG(thiol) decreases the work function of the Au(111) surface, where the magnitude of this reduction strongly depends on the length of the PEG backbone. The main contribution arises from the dipole due to the adsorption-induced charge rearrangement at the interface. Our work reveals a pronounced odd-even effect, which can be traced back to the dipole moment of the PEG(thiol) layer.

preprint2019arXiv

Maximally localized Wannier functions within the (L)APW+LO method

We present a robust algorithm that computes (maximally localized) Wannier functions (WFs) without the need of providing an initial guess. Instead, a suitable starting point is constructed automatically from so-called local orbitals which are fundamental building blocks of the basis set within (linearized) augmented planewave methods. Our approach is applied to a vast variety of materials such as metals, bulk and low-dimensional semiconductors, and complex inorganic-organic hybrid interfaces. For the interpolation of electronic single-particle energies, an accuracy in the meV range can be easily achieved. We exemplify the capabilities of our method by the calculation of the joint density of states in aluminum, (generalized) Kohn-Sham and quasi-particle band structures in various semiconductors, and the electronic structure of $β$-Ga$_2$O$_3$, including electron and hole effective masses.

preprint2019arXiv

Photoemission Signatures of Non-Equilibrium Carrier Dynamics from First Principles

Time- and angle-resolved photoemission spectroscopy (tr-ARPES) constitutes a powerful tool to inspect the dynamics and thermalization of hot carriers. The identification of the processes that drive the dynamics, however, is challenging even for the simplest systems owing to the coexistence of several relaxation mechanisms. Here, we devise a Green's function formalism for predicting the tr-ARPES spectral function and establish the origin of carrier thermalization entirely from first principles. The predictive power of this approach is demonstrated by an excellent agreement with experiments for graphene over time scales ranging from a few tens of femtoseconds up to several picoseconds. Our work provides compelling evidence of a non-equilibrium dynamics dominated by the establishment of a hot-phonon regime.