Understanding the role of oxygen-vacancy defects in Cu2O(111) from first-principle calculations

The presence of defects, such as copper and oxygen vacancies, in cuprous oxide films determines their characteristic carrier conductivity and consequently their application as semiconducting systems. There are still open questions on the induced electronic re-distribution, including the formation of polarons. Indeed, to accurately reproduce the structural and electronic properties at the cuprous oxide surface, very large slab models and theoretical approaches that go beyond the standard generalized gradient corrected density functional theory are needed. In this work we investigate oxygen vacancies formed in proximity of a reconstructed Cu2O(111) surface, where the outermost unsaturated copper atoms are removed, thus forming non-stoichiometric surface layers with copper vacancies. We address simultaneously surface and bulk properties by modelling a thick and symmetric slab, to find that hybrid exchange-correlation functionals are needed to describe the oxygen vacancy in this system. Our simulations show that the formation of oxygen vacancies is favoured in the sub-surface layer. Moreover, the oxygen vacancy leads to a splitting and left-shift of the shallow hole states in the gap, which are associated with the deficiency of copper at the surface. These findings suggest that surface electronic structure and reactivity are sensitive to the presence of oxygen vacancies, also when the latter are formed deeper within the film

Understanding the role of oxygen-vacancy defects in Cu₂O(111) from first-principle calculations

The presence of defects, such as copper and oxygen vacancies, in cuprous oxide films determines their characteristic carrier conductivity and consequently their application as semiconducting systems. There are still open questions on the induced electronic re-distribution, including the formation of polarons. Indeed, to accurately reproduce the structural and electronic properties at the cuprous oxide surface, very large slab models and theoretical approaches that go beyond the standard generalized gradient corrected density functional theory are needed. In this work we investigate oxygen vacancies formed in proximity of a reconstructed Cu$_{2}$O(111) surface, where the outermost unsaturated copper atoms are removed, thus forming non-stoichiometric surface layers with copper vacancies. We address simultaneously surface and bulk properties by modelling a thick and symmetric slab, to find that hybrid exchange-correlation functionals are needed to describe the oxygen vacancy in this system. Our simulations show that the formation of oxygen vacancies is favoured in the sub-surface layer. Moreover, the oxygen vacancy leads to a splitting and left-shift of the shallow hole states in the gap, which are associated with the deficiency of copper at the surface. These findings suggest that surface electronic structure and reactivity are sensitive to the presence of oxygen vacancies, also when the latter are formed deeper within the film

Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative—particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods

Water on hexagonal boron nitride from diffusion Monte Carlo

Despite a recent flurry of experimental and simulation studies, an accurate estimate of the interaction strength of water molecules with hexagonal boron nitride is lacking. Here we report quantum Monte Carlo results for the adsorption of a water monomer on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction energy of -84 +/- 5 meV. We use the results to evaluate the performance of several widely used density functional theory (DFT) exchange correlation functionals, and find that they all deviate substantially. Differences in interaction energies between different adsorption sites are however better reproduced by DFT

Mechanisms of adsorbing hydrogen gas on metal decorated graphene

Hydrogen is a key player in global strategies to reduce greenhouse gas emissions. In order to make hydrogen a widely used fuel, we require more efficient methods of storing it than the current standard of pressurized cylinders. An alternative method is to adsorb H2 in a material and avoid the use of high pressures. Among many potential materials, layered materials such as graphene present a practical advantage as they are lightweight. However, graphene and other 2D materials typically bind H2 too weakly to store it at the typical operating conditions of a hydrogen fuel cell, meaning that high pressure would still be required. Modifying the material, for example by decorating graphene with adatoms, can strengthen the adsorption energy of H2 molecules, but the underlying mechanisms are still not well understood. In this work, we systematically screen alkali and alkaline-earth metal decorated graphene sheets for the static thermodynamic adsorption of hydrogen gas from first principles and focus on the mechanisms of binding. We show that there are three mechanisms of adsorption on metal decorated graphene and each leads to distinctly different hydrogen adsorption structures. The three mechanisms can be described as weak van der Waals physisorption, metal adatom facilitated polarization, and Kubas adsorption. Among these mechanisms, we find that Kubas adsorption is easily perturbed by an external electric field, providing a way to tune H2 adsorption. This work is foundational and builds our understanding of H2 adsorption under idealized conditions

Exploring water adsorption on isoelectronically doped graphene using alchemical derivatives

The design and production of novel 2-dimensional materials has seen great progress in the last decade, prompting further exploration of the chemistry of such materials. Doping and hydrogenating graphene is an experimentally realised method of changing its surface chemistry, but there is still a great deal to be understood on how doping impacts on the adsorption of molecules. Developing this understanding is key to unlocking the potential applications of these materials. High throughput screening methods can provide particularly effective ways to explore vast chemical compositions of materials. Here, alchemical derivatives are used as a method to screen the dissociative adsorption energy of water molecules on various BN doped topologies of hydrogenated graphene. The predictions from alchemical derivatives are assessed by comparison to density functional theory. This screening method is found to predict dissociative adsorption energies that span a range of more than 2 eV, with a mean absolute error $<0.1$ eV. In addition, we show that the quality of such predictions can be readily assessed by examination of the Kohn-Sham highest occupied molecular orbital in the initial states. In this way, the root mean square error in the dissociative adsorption energies of water is reduced by almost an order of magnitude (down to $\sim0.02$ eV) after filtering out poor predictions. The findings point the way towards a reliable use of first order alchemical derivatives for efficient screening procedures

Coulomb Interactions between Dipolar Quantum Fluctuations in van der Waals Bound Molecules and Materials

Mutual Coulomb interactions between electrons lead to a plethora of interesting physical and chemical effects, especially if those interactions involve many fluctuating electrons over large spatial scales. Here, we identify and study in detail the Coulomb interaction between dipolar quantum fluctuations in the context of van der Waals complexes and materials. Up to now, the interaction arising from the modification of the electron density due to quantum van der Waals interactions was considered to be vanishingly small. We demonstrate that in supramolecular systems and for molecules embedded in nanostructures, such contributions can amount to up to 6 kJ/mol and can even lead to qualitative changes in the long-range vdW interaction. Taking into account these broad implications, we advocate for the systematic assessment of so-called Coulomb singles in large molecular systems and discuss their relevance for explaining several recent puzzling experimental observations of collective behavior in nanostructured materials

Exploring dissociative water adsorption on isoelectronically BN doped graphene using alchemical derivatives

The design and production of novel 2-dimensional materials have seen great progress in the last decade, prompting further exploration of the chemistry of such materials. Doping and hydrogenating graphene are an experimentally realised method of changing its surface chemistry, but there is still a great deal to be understood on how doping impacts on the adsorption of molecules. Developing this understanding is key to unlocking the potential applications of these materials. High throughput screening methods can provide particularly effective ways to explore vast chemical compositions of materials. Here, alchemical derivatives are used as a method to screen the dissociative adsorption energy of water molecules on various BN doped topologies of hydrogenated graphene. The predictions from alchemical derivatives are assessed by comparison to density functional theory. This screening method is found to predict dissociative adsorption energies that span a range of more than 2 eV, with a mean absolute error <0.1 eV. In addition, we show that the quality of such predictions can be readily assessed by examination of the Kohn-Sham highest occupied molecular orbital in the initial states. In this way, the root mean square error in the dissociative adsorption energies of water is reduced by almost an order of magnitude (down to ∼0.02 eV) after filtering out poor predictions. The findings point the way towards a reliable use of first order alchemical derivatives for efficient screening procedures

Tuning dissociation using isoelectronically doped graphene and hexagonal boron nitride: water and other small molecules

Novel uses for 2-dimensional materials like graphene and hexagonal boron nitride (h-BN) are being frequently discovered especially for membrane and catalysis applications. Still however, a great deal remains to be understood about the interaction of environmentally and industrially elevant molecules such as water with these materials. Taking inspiration from advances in hybridising graphene and h-BN, we explore using density functional theory, the dissociation of water, hydrogen, methane, and methanol on graphene, h-BN, and their isoelectronic doped counterparts: BN doped graphene and C doped h-BN. We find that doped surfaces are considerably more reactive than their pristine counterparts and by comparing the reactivity of several small molecules we develop a general framework for dissociative adsorption. From this a particularly attractive consequence of isoelectronic doping emerges: substrates can be doped to enhance their reactivity specifically towards either polar or non-polar adsorbates. As such, these substrates are potentially viable candidates for selective catalysts and membranes, with the implication that a range of tuneable materials can be designed

Properties of the water to boron nitride interaction: from zero to two dimensions with benchmark accuracy

Molecular adsorption on surfaces plays an important part in catalysis, corrosion, desalination, and various other processes that are relevant to industry and in nature. As a complement to experiments, accurate adsorption energies can be obtained using various sophisticated electronic structure methods that can now be applied to periodic systems. The adsorption energy of water on boron nitride substrates, going from zero to 2-dimensional periodicity, is particularly interesting as it calls for an accurate treatment of polarizable electrostatics and dispersion interactions, as well as posing a practical challenge to experiments and electronic structure methods. Here, we present reference adsorption energies, static polarizabilities, and dynamic polarizabilities, for water on BN substrates of varying size and dimension. Adsorption energies are computed with coupled cluster theory, fixed-node quantum Monte Carlo (FNQMC), the random phase approximation (RPA), and second order M{\o}ller-Plesset (MP2) theory. These explicitly correlated methods are found to agree in molecular as well as periodic systems. The best estimate of the water/h-BN adsorption energy is $-107\pm7$ meV from FNQMC. In addition, the water adsorption energy on the BN substrates could be expected to grow monotonically with the size of the substrate due to increased dispersion interactions but interestingly, this is not the case here. This peculiar finding is explained using the static polarizabilities and molecular dispersion coefficients of the systems, as computed from time-dependent density functional theory (DFT). Dynamic as well as static polarizabilities are found to be highly anisotropic in these systems. In addition, the many-body dispersion method in DFT emerges as a particularly useful estimation of finite size effects for other expensive, many-body wavefunction based methods