Unexpected Self-Assembly of Nanographene Oxide Membranes upon Electron Beam Irradiation for Ultrafast Ion Sieving.

Nanographene oxide (nGO) flakes-graphene oxide with a lateral size of ≈100 nm or less-hold great promise for superior flux and energy-efficient nanofiltration membranes for desalination and precise ionic sieving owing to their unique high-density water channels with less tortuousness. However, their potential usage is currently limited by several challenges, including the tricky self-assembly of nano-sized flakes on substrates with micron-sized pores, severe swelling in aqueous solutions, and mechanical instability. Herein, the successful fabrication of a robust membrane stacked with nGO flakes on a substrate with a pore size of 0.22 µm by vacuum filtration is reported. This membrane achieved an unprecedented water permeance above 819.1 LMH bar-1, with a high rejection rate of 99.7% for multivalent metal ions. The nGO flakes prepared using an electron beam irradiation method, have uniquely pure hydroxyl groups and abundant aromatic regions. The calculations revealed the strong hydrogen bonds between two nGO flakes, which arise from hydroxyl groups, coupled with hydrophobic aromatic regions, greatly enhance the stability of stacked flakes in aqueous solutions and increase their effective lateral size. The research presents a simple yet effective approach toward the fabrication of advanced 2D nanographene membranes with superior performance for ion sieving applications.

p-block germanenes as a promising electrocatalysts for the oxygen reduction reaction.

The oxygen reduction reaction (ORR), a pivotal process in hydrogen fuel cells crucial for enhancing fuel cell performance through suitable catalysts, remains a challenging aspect of development. This study explores the catalytic potential of germanene on Al (111), taking advantage of the successful preparation of stable reconstructed germanene layers on Al (111) and the excellent catalytic performance exhibited by germanium-based nanomaterials. Through first-principles calculations, we demonstrate that the O2 molecule can be effectively activated on both freestanding and supported germanene nanosheets, featuring kinetic barriers of 0.40 and 0.04 eV, respectively. The presence of the Al substrate not only significantly enhances the stability of the reconstructed germanene but also preserves its exceptional ORR catalytic performance. These theoretical findings offer crucial insights into the substrate-mediated modulation of germanene stability and catalytic efficiency, paving the way for the design of stable and efficient ORR catalysts for future applications.

Confined and spontaneously transformed oxidation structures due to the intrinsic heterogeneous surface morphology of C3N monolayer.

Identifying the oxidation structure of two-dimensional interfaces is crucial to improve surface chemistry and electronic properties. Beyond graphene with only phenyl rings, a novel carbon-nitrogen material, C3N, presents an intrinsic heterogeneous surface morphology where each phenyl ring is encircled by six nitrogen atoms, yet its atomistic oxidation structure remains unclear. Here, combining a series of density functional theory calculations and ab initio molecular dynamics simulations, we demonstrate that thermodynamically favorable oxidation loci are confined to the phenyl ring, and kinetic transformations of oxidation structures are feasible along the phenyl ring, whereas those toward nitrogen atoms are proven to be extremely difficult. These results are attributed to the lower barrier of oxygen atom migration along the phenyl ring, while the significantly high barriers toward nitrogen atoms are due to the heterogeneous potential energy surface for oxygen-C3N interaction. This work highlights the significance of surface morphology on the characteristics of oxidation structure, offering insights into tunable electronic properties via confined interfacial oxidation.

Magnesium Ion Gated Ion Rejection through Carboxylated Graphene Oxide Nanopore: A Theoretical Study.

While nanoporous graphene oxide (GO) is recognized as one of the most promising reverse osmosis desalination membranes, limited attention has been paid to controlling desalination performance through the large GO pores, primarily due to significant ion leakage resulting in the suboptimal performance of these pores. In this study, we employed a molecular dynamics simulation approach to demonstrate that Mg2+ ions, adhered to carboxylated GO nanopores, can function as gates, regulating the transport of ions (Na+ and Cl-) through the porous GO membrane. Specifically, the presence of divalent cations near a nanopore reduces the concentration of salt ions in the vicinity of the pore and prolongs their permeation time across the pore. This subsequently leads to a notable enhancement in salt rejection rates. Additionally, the ion rejection rate increases with more adsorbed Mg2+ ions. However, the presence of the adsorbed Mg2+ ions compromises water transport. Here, we also elucidate the impact of graphene oxidation degree on desalination. Furthermore, we design an optimal combination of adsorbed Mg2+ ion quantity and oxidation degree to achieve high water flux and salt rejection rates. This work provides valuable insights for developing new nanoporous graphene oxide membranes for controlled water desalination.

Exceptional metal-semiconductor-metal transition of lead apatites via oxygen defect tuning.

Lead apatites, distinguished and compelling bulk materials with the stoichiometric arrangement as Pb10(POx)6Oy, are renowned for their structural complexity. Recently, the discovery of possible room-temperature superconductivity under ambient pressure in copper-substituted lead apatites has engendered considerable interest within both the physics community and beyond. Nevertheless, exploration of pristine Pb10(POx)6Oy parent structures has hitherto remained elusive. In this study, we employ density functional theory (DFT) calculations to investigate the effects of oxygen defects on the electronic structures of Pb10(POx)6Oy and Pb9Cu(POx)6Oy. We scrutinize two distinct categories of defects: oxygen atoms enmeshed within POx groups (Ox) and solitary oxygen atoms (Oy). Our investigation uncovers a profound influence of these defects on the band structure. Specifically, the introduction of Oy defects prompts a remarkable transition in Pb10(PO4)6Oy from a metal to semiconductor to metal state, accompanied by pivotal shifts in the principal electronic contributors from p orbitals of Oy to those of Pb atoms. Furthermore, the introduction of Ox defects in Pb10(POx)6O1 engenders metamorphosis in the band structure, transmuting it from a semiconductor to a metallic state. Significantly, our findings pinpoint the suitable range of x in the Pb10(POx)6O1 configuration as lying between 2 and 4. Additionally, our study also demonstrates that the oxygen defects (Ox/Oy) do not affect the metallic properties of copper-substituted lead apatites. This study elucidates the significant role of oxygen defects in modulating the electronic properties of apatite materials, offering insights into potential interdisciplinary applications. This establishes a crucial link between material composition and electronic behavior, revealing key mechanisms for engineering functionality in lead apatites and other advanced materials.

Rational Design of Full-Color Fluorescent C3N Quantum Dots.

Carbon-based quantum dots (QDs) exhibit unique photoluminescence due to size-dependent quantum confinement, giving rise to fascinating full-color emission properties. Accurate emission calculations using time-dependent density functional theory are a time-costing and expensive process. Herein, we employed an artificial neural network (ANN) combined with statistical learning to establish the relationship between geometrical/electronic structures of ground states and emission wavelength for C3N QDs. The emission energy of these QDs can be doubly modulated by size and edge effects, which are governed by the number of C4N2 rings and the CH group, respectively. Moreover, these two structural characteristics also determine the phonon vibration mode of C3N QDs to harmonize the emission intensity and lifetime of hot electrons in the electron-hole recombination process, as indicated by nonadiabatic molecular dynamics simulation. These computational results provide a general approach to atomically precise design the full-color fluorescent carbon-based QDs with targeted functions and high performance.

In situ and General Multidentate Ligand Passivation Achieves Efficient and Ultra-Stable CsPbX3 Perovskite Quantum Dots for White Light-Emitting Diodes.

Inorganic CsPbX3 perovskite quantum dots (PeQDs) show great potential in white light-emitting diodes (WLEDs) due to excellent optoelectronic properties, but their practical application is hampered by low photoluminescence quantum yield (PLQY) and especially poor stability. Herein,  we developed an in-situ and general multidentate ligand passivation strategy that allows for CsPbX3 PeQDs not only near-unit PLQY, but significantly improved stability against storage, heat, and polar solvent. The enhanced optical property arises from high effectiveness of the multidentate ligand, diethylenetriaminepentaacetic acid (DTPA) with five carboxyl groups, in passivating uncoordinated Pb2+ defects and suppressing nonradiative recombination. First-principles calculations reveal that the excellent stability is attributed to tridentate binding mode of DTPA that remarkably boosts the adsorption capacity to PeQD core. Finally, combining the green and red PeQDs with blue chip,  we demonstrated highly luminous WLEDs with distinctly enhanced operation stability, a wide color gamut of 121.3% of national television system committee, standard white light of (0.33,0.33) in CIE 1931, and tunable color temperatures from warm to cold white light readily by emitters' ratio. This study provides an operando yet general approach to achieve efficient and stable PeQDs for WLEDs and accelerates their progress to commercialization.

Graphene oxide-based large-area dynamic covalent interfaces.

Dynamic materials, being capable of reversible structural adaptation in response to the variation of external surroundings, have experienced significant advancements in the past several decades. In particular, dynamic covalent materials (DCMs), where the dynamic covalent bonds (DCBs) can reversibly break and reform under defined conditions, present superior dynamic characteristics, such as self-adaptivity, self-healing and shape memory. However, the dynamic characteristics of DCBs are mainly limited within the length scale of covalent bonds, due to the local position exchange or the inter-distance variation between the chemical compositions involved in the reversible covalent reactions. In this minireview, a discussion regarding the realization of long-range migration of chemical compositions along the interfaces of graphene oxide (GO)-based materials via the spatially connected and consecutive occurrence of DCB-based reversible covalent reactions is presented, and the interfaces are termed "large-area dynamic covalent interfaces (LDCIs)". The effective strategies, including water adsorption, interfacial curvature and metal-substrate support, as well as the potential applications of LDCIs in water dissociation and humidity sensing are summarized. Additionally, we also give an outlook on potential strategies to realize LDCIs on other 2D carbon-based materials, including the interfacial morphology and periodic element doping. This minireview provides insights into the realization of LDCIs on a wider range of 2D materials, and offers a theoretical perspective for advancing materials with long-range dynamic characteristics and improved performance, including controlled drug delivery/release and high-efficiency (bio)sensing.

Low lattice thermal conductivities and good thermoelectric performance of hexagonal antiperovskites X(Ba & Sr)3BiN with quartic anharmonicity.

Antiperovskites are a burgeoning class of semiconducting materials that showcase remarkable optoelectronic properties and catalytic properties. However, there has been limited research on their thermoelectric properties. Combining first-principles calculations, self-consistent phonon theory and the Boltzmann transport equation, we have discovered that the hexagonal antiperovskites X(Ba & Sr)3BiN exhibit strong quartic lattice anharmonicity, where the anharmonic vibrations of the light N atoms primarily affect the lattice thermal conductivity (κL) along the c-axis direction. As a result, the lattice thermal conductivities along the a(b)-axis direction are low. At 300 K, the κL values of Ba3BiN and Sr3BiN are only 1.27 W m-1 K-1 and 2.24 W m-1 K-1, respectively. Moreover, near the valence band maximum, the orbitals of the N atoms dominate. This dominance allows Sr3BiN to achieve high power factor under p-type doping, resulting in an impressive thermoelectric figure of merit (ZT) of 0.94 along the c-axis direction at 800 K. In the a(b)-axis direction, at 800 K, n-type doped Ba3BiN exhibits a ZT value of 1.47, surpassing that of traditional thermoelectric materials. Our research elucidates that the hexagonal antiperovskites X(Ba & Sr)3BiN represent a category of potential thermoelectric materials with pronounced anisotropy, low thermal conductivity, and high thermoelectric performance.

Solar driven CO2 hydrogenation to HCOOH on (TiO2)n (n = 1-6) atomic clusters.

Artificial photosynthesis is a crucial reaction that addresses energy and environmental challenges by converting CO2 into fuels and value-added chemicals. However, efficient catalytic activity using earth-abundant materials can be challenging due to intrinsic limitations. Herein, we explore neutral (TiO2)n (n = 1-6) atomic clusters for CO2 hydrogenation via comprehensive ab initio calculations combined with time-dependent functional theory. Our results show that these (TiO2)n clusters exhibit outstanding thermodynamic stabilities and decent surficial activities for CO2 activation and H2 dissociation, both of which possess kinetic barriers down to 0-0.74 eV. We establish a relationship between the binding strength of *CO2 species and electron characterization for these (TiO2)n clusters. These clusters, which have a wide energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccpied molecular orbital (LUMO) that allows them to harvest the solar light in the ultraviolet regime, enabling efficient catalysis for driving the catalysis of CO2 conversion. They provide exclusive reaction channels and high selectivity for yielding HCOOH products via the carboxyl mechanism, involving the kinetic barrier of the limiting step of 0.74-1.25 eV. We also investigated the substrate effect on supported (TiO2)n clusters, with non-metallic substrates featuring inert surfaces serving as suitable options for anchoring (TiO2)n clusters while preserving their intrinsic activity and selectivity. These computational results have significant implications not only for meeting energy demands but also for mitigating carbon emissions by utilizing CO2 as an alternative feedstock rather than considering it solely as a greenhouse gas.

Unexpectedly resisting protein adsorption on self-assembled monolayers terminated with two hydrophilic hydroxyl groups.

OH-terminated self-assembled monolayers, as protein-resistant surfaces, have significant potential in biocompatible implant devices, which can avoid or reduce adverse reactions caused by protein adhesion to biomaterial surfaces, such as thrombosis, immune response and inflammation. Here, molecular dynamics simulations were performed to evaluate the degree of protein adsorption on the self-assembled monolayer terminated with two hydrophilic OH groups ((OH)2-SAM) at packing densities (Σ) of 4.5 nm-2 and 6.5 nm-2, respectively. The results show that the structure of the (OH)2-SAM itself, i.e., a nearly perfect hexagonal-ice-like hydrogen bond structure in the OH matrix of the (OH)2-SAM at Σ = 4.5 nm-2 sharply reduces the number of hydrogen bonds (i.e., 0.7 ± 0.27) formed between the hydrophobic (OH)2-SAM surface and protein. While for Σ = 6.5 nm-2, the hydrophilic (OH)2-SAM surface can provide more hydrogen bonding sites to form hydrogen bonds (i.e., 6.2 ± 1.07) with protein. The number of hydrogen bonds formed between the (OH)2-SAM and protein at Σ = 6.5 nm-2 is ∼8 times higher than that at Σ = 4.5 nm-2, reflecting the excellent resistance to protein adsorption exhibited by the structure of the (OH)2-SAM itself at Σ = 4.5 nm-2. Compared with a traditional physical barrier effect formed by a large number of hydrogen bonds between the (OH)2-SAM and water at Σ = 6.5 nm-2, the structure of the (OH)2-SAM itself at Σ = 4.5 nm-2 proposed in this study significantly improves the performance of the (OH)2-SAM resistance to protein adsorption, which provides new insights into the mechanism of resistance to protein adsorption on the (OH)2-SAM.

High-Yield Synthesis of Sodium Chlorides of Unconventional Stoichiometries.

Abnormal salt crystals with unconventional stoichiometries, such as Na2 Cl, Na3 Cl, K2 Cl, and CaCl crystals that have been explored in reduced graphene oxide membranes (rGOMs) or diamond anvil cells, hold great promise in applications due to their unique electronic, magnetic, and optical properties predicted in theory. However, the low content of these crystals, only <1% in rGOM, limits their research interest and utility in applications. Here, a high-yield synthesis of 2D abnormal crystals with unconventional stoichiometries is reported, which is achieved by applying negative potential on rGOM. A more than tenfold increase in the abnormal Na2 Cl crystals is obtained using a potential of -0.6 V, resulting in an atomic content of 13.4 ± 4.7% for Na on rGOM. Direct observations by transmission electron microscopy and piezoresponse force microscopy demonstrates a unique piezoelectric behavior arising from 2D Na2 Cl crystals with square structure. The output voltage increases from 0 to ≈180 mV in the broad 0-150° bending angle regime, which meets the voltage requirement of most nanodevices in realistic applications. Density functional theory calculations reveal that the applied negative potential of the graphene surface can strengthen the effect of the Na+ -π interaction and reduce the electrostatic repulsion between cations, making more Na2 Cl crystals formed.

Ultralow lattice thermal conductivities and excellent thermoelectric properties of hypervalent triiodides XI3 (X = Rb, Cs) discovered by machine learning method.

Thermal conductivity and power factor are key factors in evaluating heat transfer performance and designing thermoelectric conversion devices. To search for materials with ultralow thermal conductivity and a high power factor, we proposed a set of universal statistical interaction descriptors (SIDs) and developed accurate machine learning models for the prediction of thermoelectric properties. For lattice thermal conductivity prediction, the SID-based model achieved the state-of-the-art results with an average absolute error of 1.76 W m-1 K-1. The well-performing models predicted that hypervalent triiodides XI3 (X = Rb, Cs) have ultralow thermal conductivities and high power factors. Combining first-principles calculations, the self-consistent phonon theory, and the Boltzmann transport equation, we obtained the anharmonic lattice thermal conductivities of 0.10 and 0.13 W m-1 K-1 for CsI3 and RbI3 in the c-axis direction at 300 K, respectively. Further studies show that the ultralow thermal conductivity of XI3 arises from the competition of vibrations between alkali metal atoms and halogen atoms. In addition, at 700 K, the thermoelectric figure of merit ZT values of CsI3 and RbI3 are 4.10 and 1.52, respectively, at the optimal hole doping level, which indicates hypervalent triiodides are potential high performance thermoelectric materials.

Nonselective Intercalation of G-Quadruplex-Targeting Ligands into Double-Stranded DNA Quantified by Single-Molecule Stretching.

Most G-quadruplex (G4)-targeting ligands reported so far contain planar heteroaromatic groups and can intercalate into adjacent base pairs of double-stranded DNA (dsDNA). However, quantitative data on the binding number γ (ligands/bp) of G4 ligands that intercalate into long dsDNA remain lacking, which are essential for understanding the selectivity of G4 ligands. Here, using a single-molecule stretching assay based on the lengthening of dsDNA, we analyzed the dissociation constants and the binding number of eight most commonly used G4 ligands that intercalate into dsDNA. We showed that five ligands (CX-5461, BRACO-19, RHPS4, TrisQ, and Phen-DC3) intercalate into dsDNA avidly (Kd = 0.5-2.1 µM, saturated γ > 0.2 ligands/bp), which was similar to the typical dsDNA intercalator EB. Two bisquinolines, PDS and 360A, showed moderate intercalation ability (Kd = 22.5 and 48.7 µM) and γ < 0.01 ligands/bp in the presence of 1 µM ligands. Porphyrin NMM showed no intercalative binding even at 200 µM. Molecular docking and molecular dynamics simulations were carried out to further evaluate the intercalative binding of these G4 ligands with dsDNA by calculating the binding energies and π-π stacking probability.

Locally spontaneous dynamic oxygen migration on biphenylene: a DFT study.

The dynamic oxygen migration at the interface of carbon allotropes dominated by the periodic hexagonal rings, including graphene and carbon nanotubes, has opened up a new avenue to realize dynamic covalent materials. However, for the carbon materials with hybrid carbon rings, such as biphenylene, whether the dynamic oxygen migration at its interface can still be found remains unknown. Using both density functional theory calculations and machine-learning-based molecular dynamics (MLMD) simulations, we found that the oxygen migration departing away from the four-membered carbon (C4) ring is hindered, and the oxygen atom prefers to spontaneously migrate toward/around the C4 ring. This locally spontaneous dynamic oxygen migration on the biphenylene is attributed to a high barrier of about 1.5 eV for the former process and a relatively low barrier of about 0.3 eV for the latter one, originating from the enhanced activity of the C-O bond near/around the C4 ring due to the hybrid carbon ring structure. Moreover, the locally spontaneous dynamic oxygen migration is further confirmed by MLMD simulations. This work sheds light on the potential of biphenylene as a catalyst for spatially controlled energy conversion and provides the guidance for realizing the dynamic covalent interface at other carbon-based or two-dimensional materials.

Ultralow lattice thermal conductivity of binary compounds A2B (A = Cs, Rb & B = Se, Te) with higher-order anharmonicity correction.

By employing first-principles calculations that integrate self-consistent phonon theory and the Boltzmann transport equation, we have delved into the thermal transport characteristics of hexagonal anisotropic materials A2B (A = Cs, Rb and B = Se, Te). Our computational results have disclosed that these A2B materials exhibit ultralow lattice thermal conductivity (κL) at room temperature. Specifically, in the case of Cs2Te, the κL values are a mere 0.15 W m-1 K-1 in the a(b) direction and 0.22 W m-1 K-1 in the c direction, both markedly less than the thermal conductivity of quartz glass, a conventional thermoelectric material (0.9 W m-1 K-1). Importantly, our calculations encompass higher-order anharmonic effects while computing the lattice thermal conductivities of these materials. This is essential since pronounced anharmonicity leads to a decrease in phonon group velocity, and consequently, lowers the κL values. Our results establish a theoretical foundation for exploring the thermal transport characteristics of anisotropic materials with substantial anharmonicity. Furthermore, the binary compounds A2B proffer a gamut of possibilities for a wide range of applications in thermoelectrics and thermal management, owing to their ultralow lattice thermal conductivity.

The single-atom catalytic activity of the hydrogen evolution reaction of the experimentally synthesized boridene 2D material: a density functional theory study.

CONTEXT: Previous theoretical studies have suggested that two-dimensional (2D) MBene materials might display adequate monatomic catalytic activity for the hydrogen evolution reaction (HER). Recently, a study reported the experimental synthesis of a 2D MBene (Mo4/3B2), re-defined as boridene, albeit no effort has been devoted to explore the single-atom catalytic activity for HER of this experimentally synthesized 2D material. Therefore, we herein investigate the single-atom HER performance of the boridene. Interestingly, with Mo defects mixed with single Au and Zn atoms shows excellent hydrogen evolution performance, and the change in the Gibbs free energy ([Formula: see text]) value is close to 0 eV, which can even match the performance of Pt-based materials. Through analysis of the charge density difference and density of states, the mechanism affecting the HER performance is explained at the electronic level. This work provides a new direction for the use of the Mo4/3B2 monolayer two-dimensional materials in the field of single-atom catalysis for HER. METHODS: This study used the DFT calculations in Vienna ab initio simulation package. The GGA-Perdew-Burke-Ernzerhof functional with DFT-D2 correction is used to describe the exchange-correlation interactions. The projection augmented wave is used with the plane wave cutoff of 500 eV. The convergences of energy and force are 10-5 eV and 0.01 eV/Å, respectively. A vacuum layer with a height of 20 Å is set in the Z direction. For geometry optimization, self-consistent, and DOS calculations, the k-point grids sampled in Brillouin zones are 3 × 3 × 1, 9 × 9 × 1, and 9 × 9 × 1, respectively. The AIMD simulation is performed in the canonical ensemble (NVT), and the temperature was maintained at 300 K by Nosé-Hoover thermostats with a time step of 2.0 fs.

Unexpectedly spontaneous water dissociation on graphene oxide supported by copper substrate.

Water dissociation is of fundamental importance in scientific fields and has drawn considerable interest in diverse technological applications. However, the high activation barrier of breaking the OH bond within the water molecule has been identified as the bottleneck, even for the water adsorbed on the graphene oxide (GO). Herein, using the density functional theory calculations, we demonstrate that the water molecule can be spontaneously dissociated on GO supported by the (111) surface of the copper substrate (Copper-GO). This process involves a proton transferring from water to the interfacial oxygen group, and a hydroxide covalently bonding to GO. Compared to that on GO, the water dissociation barrier on Copper-GO is significantly decreased to be less than or comparable to thermal fluctuations. This is ascribed to the orbital-hybridizing interaction between copper substrate and GO, which enhances the reaction activity of interfacial oxygen groups along the basal plane of GO for water dissociation. Our work provides a novel strategy to access water dissociation via the substrate-enhanced reaction activity of interfacial oxygen groups on GO and indicates that the substrate can serve as an essential key to tuning the catalytic performance of various two-dimensional material devices.

Unexpectedly efficient ion desorption of graphene-based materials.

Ion desorption is extremely challenging for adsorbents with superior performance, and widely used conventional desorption methods involve high acid or base concentrations and large consumption of reagents. Here, we experimentally demonstrate the rapid and efficient desorption of ions on magnetite-graphene oxide (M-GO) by adding low amounts of Al3+. The corresponding concentration of Al3+ used is reduced by at least a factor 250 compared to conventional desorption method. The desorption rate reaches ~97.0% for the typical radioactive and bivalent ions Co2+, Mn2+, and Sr2+ within ~1 min. We achieve effective enrichment of radioactive 60Co and reduce the volume of concentrated 60Co solution by approximately 10 times compared to the initial solution. The M-GO can be recycled and reused easily without compromising its adsorption efficiency and magnetic performance, based on the unique hydration anionic species of Al3+ under alkaline conditions. Density functional theory calculations show that the interaction of graphene with Al3+ is stronger than with divalent ions, and that the adsorption probability of Al3+ is superior than that of Co2+, Mn2+, and Sr2+ ions. This suggests that the proposed method could be used to enrich a wider range of ions in the fields of energy, biology, environmental technology, and materials science.

Remarkably enhanced dynamic oxygen migration on graphene oxide supported by copper substrate.

The dynamic covalent properties of graphene oxide (GO) are of fundamental interest to a broad range of scientific areas and technological applications. It remains a challenge to access feasible dynamic reactions for reversibly breaking/reforming the covalent bonds of oxygen functional groups on GO, although these reactions can be induced by photonic or mechanical routes, or mediated by adsorbed water. Here, using density functional theory calculations, we demonstrate the remarkably enhanced dynamic oxygen migration along the basal plane of GO supported by copper substrate (GO@copper), with C-O bond breaking reactions and proton transfer between neighboring epoxy and hydroxyl groups. Compared to reactions on GO, the energy barriers of oxygen migrations on GO@copper are sharply reduced to be less than or comparable to thermal fluctuations, and meanwhile the crystallographic match between GO and copper substrate induces new oxygen migration paths on GO@copper. This work sheds light on understanding of the metal substrate-enhanced dynamic properties of GO, and evidences the strategy to tune the activity of two-dimensional-interfacial oxygen groups for various potential applications.