Thermal conductivities of monolayer graphene oxide from machine learning molecular dynamics simulations

Graphene oxide (GO) exhibits rich chemical heterogeneity that strongly influences its structural, thermal, and mechanical properties, yet quantitatively linking reduction chemistry to heat transport remains challenging. In this study, we develop a machine-learned neuroevolution potential (NEP) trained on an existing density functional theory dataset [El-Machachi et al., Angew. Chem., Int. Ed. 63, e202410088 (2024)], achieving reasonable accuracy at a computational cost much lower than the existing machine-learned and empirical potentials. Leveraging this potential, we perform large-scale molecular dynamics (MD) simulations to model the thermal reduction of GO across realistic structural domains. Using the homogeneous nonequilibrium MD method with a proper quantum-statistical correction scheme, we find that reduced GO exhibits strongly suppressed thermal conductivities, ranging from a few to tens of Wm-1 K-1, substantially lower than pristine GO without defects and far below graphene. Moreover, the thermal conductivity of reduced GO increases moderately with increasing OH/O ratio, except at the highest oxidation level (O/C = 0.5), where this trend inverts, while decreasing significantly with increasing O/C ratio, a trend strongly correlated with the fraction of recovered graphene-like structures. Our study provides a computationally tractable and predictive atomistic machine learning framework for exploring how chemical structure governs heat transport in heterogeneous carbon materials.