Search

Ferroelectrics are of practical interest for non-volatile data storage due to their reorientable, crystallographically defined polarization. Yet efforts to integrate conventional ferroelectrics into ultrathin memories have been frustrated by film-thickness limitations, which impede polarization reversal under low applied voltage. Wurtzite materials, including magnesium-substituted zinc oxide (Zn,Mg)O, have been shown to exhibit scalable ferroelectricity as thin films. In this work, we explain the origins of ferroelectricity in (Zn,Mg)O, showing that large strain fluctuations emerge locally in (Zn,Mg)O and can reduce local barriers to ferroelectric switching by more than 40%. We provide concurrent experimental and computational evidence of these effects by demonstrating polarization switching in ZnO/(Zn,Mg)O/ZnO heterostructures featuring built-in interfacial strain gradients. These results open up an avenue to develop scalable ferroelectrics by controlling strain fluctuations atomistically., Comment: 12 pages, 5 figures

The interface between NiCr alloys and FLiNaK molten salt exhibits complex corrosion behavior, mainly driven by intricate chemical interactions involving Cr and F$\mathrm{^-}$ ions. Understanding these dynamic reactions is crucial for developing effective corrosion mitigation strategies to ensure the long-term durability of Ni-based alloy components in molten salt technologies. However, obtaining molecular-level understanding through experiments is challenging. To address this, we utilize reactive molecular dynamics simulations enabled by a reactive force field, ReaxFF, to investigate detailed reaction dynamics at the atomic level. Since there is currently no available force field involving fluoride salt and Ni-based alloys, we first present the development of the ReaxFF parameter set for Ni/Cr/F/Li/Na/K based on extensive first-principles calculations. With this force field, we achieve a strong agreement for the structure of FLiNaK molten salt by comparing the pair distribution functions with experimental and simulation results. Furthermore, it successfully reproduces the experimental phenomenon of Cr dissolution in fluoride salt, with the corrosion rate depending on the alloy and salt compositions. Particularly, it reveals that increasing the concentration of Li can enhance the formation of a compact double layer, mitigating Cr dissolution. This work enables a fundamental understanding of the interfacial behavior between fluoride salt and NiCr alloys.

All-inorganic halide perovskites have received a lot of attention as attractive alternatives to overcome the stability issues of hybrid halide perovskites that are commonly associated with organic cations. To find a compromise between the optoelectronic properties of CsPbI$_{3}$ and CsPbBr$_{3}$, perovskites with CsPb(Br$_{\rm{x}}$I$_{\rm{1-x}}$)$_{3}$ mixed compositions are commonly used. An additional benefit is that, without sacrificing the optoelectronic properties for applications such as solar cells or LEDs, small amounts of Br in CsPbI$_{3}$ can prevent the inorganic perovskite from degrading to a photoinactive nonperovskite yellow phase. Despite indications that strain in the perovskite lattice plays a role in the stabilization of the material, a full understanding of such strain is lacking. Here we develop a reactive force field (ReaxFF) for perovskites starting from our previous work for CsPbI$_{3}$, we extend this force field to CsPbBr$_{3}$ and mixed CsPb(Br$_{\rm{x}}$I$_{\rm{1-x}}$)$_{3}$ compounds. This force field is used in large-scale molecular dynamics simulations to study perovskite phase transitions and the internal ion dynamics associated with the phase transitions. We find that an increase of the Br content lowers the temperature at which the perovskite reaches a cubic structure. Specifically, by substituting Br for I, the smaller ionic radius of Br induces a strain in the lattice that changes the internal dynamics of the octahedra. Importantly, this effect propagates through the perovskite lattice ranging up to distances of 2 nm, explaining why small concentrations of Br in CsPb(Br$_{\rm{x}}$I$_{\rm{1-x}}$)$_{3}$ (x $\leq$ 1/4) have a significant impact on the phase stability of mixed halide perovskites., Comment: 23 pages, 6 figures

Electrocatalysis is a catalytic process where the rate of an electrochemical reaction occurring at the electrode-electrolyte interface can be controlled by varying the electrical potential. Electrocatalysis can be applied to generate hydrogen which can be stored for future use in fuel cells for clean electricity. The use of solid oxide in electrocatalysis specially in hydrogen evolution reaction is promising. However, further improvements are essential in order to meet the ever-increasing global energy demand. Improvement of the performance of these high energy chemical systems is directly linked to the understanding and improving the complex physical and chemical phenomena and exchanges that take place at their different interfaces. To enable large length and time scale atomistic simulations of solid oxide electrocatalysis for hydrogen generation, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium (BZY20). All parameters for the eReaxFF were optimized to reproduce quantum mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, water adsorption, water splitting and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed zero-voltage molecular dynamics simulations to observe water adsorption and the eventual hydrogen production. Based on our simulation results, we conclude that this force field sets a stage for the introduction of explicit electron concept in order to simulate electron conductivity, electron leakage and non-zero-voltage effects on hydrogen generation. Overall, we demonstrate how atomistic-scale simulations can enhance our understanding of processes at interfaces in solid oxide materials., Comment: 1 file, 20 pages, 15 figures

The properties of $\mathrm{MoS_2}$ can be tuned or optimized through doping. In particular, Ni doping has been shown to improve the performance of $\mathrm{MoS_2}$ for various applications, including catalysis and tribology. To enable investigation of Ni-doped $\mathrm{MoS_2}$ with reactive molecular dynamics simulations, we developed a new ReaxFF force field to describe this material. The force field parameters were optimized to match a large set of density-functional theory (DFT) calculations of 2H-$\mathrm{MoS_2}$ doped with Ni, at four different sites (Mo-substituted, S-substituted, octahedral intercalation, and tetrahedral intercalation), under uniaxial, biaxial, triaxial, and shear strain. The force field was evaluated by comparing ReaxFF- and DFT-relaxed structural parameters, the tetrahedral/octahedral energy difference in doped 2H, energies of doped 1H and 1T monolayers, and doped 2H structures with vacancies. We demonstrated the application of the force field with reactive simulations of sputtering deposition and annealing of Ni-doped $\mathrm{MoS_2}$ films. Results show that the developed force field can successfully model the phase transition of Ni-doped $\mathrm{MoS_2}$ from amorphous to crystalline. The newly developed force field can be used in subsequent investigations to study the properties and behavior of Ni-doped $\mathrm{MoS_2}$ using reactive molecular dynamics simulations., Comment: 41 pages, 11 figures, 3 tables

Layered chalcogenides including Bi2Te3, Sb2Te3, Bi-Sb-Te ternary alloys and heterostructures are known as great thermoelectric, topological insulators and recently highlighted as plasmonic building blocks beyond noble metals. Here, we conduct a joint in situ transmission electron microscopy (in situ TEM) and density functional theory (DFT) calculations to investigate the temperature dependent nanoscale dynamics, interfacial properties and further identifying the role of native defects and edge configurations in anisotropic sublimations of Bi2Te3-Sb2Te3 in-plane heterostructure and Sb2-xBixTe3 alloy. Structural dynamics including edge evolution, formation, expansion, and coalescence of thermally induced polygonal nanopores are reported. The nanopores appear to be initiated by preferential dissociation of chalcogenide species (Te) from the center, heterointerface and edges in the heterostructure and only from the outer edges in the alloy counterpart. This results in a reduced thermal stability and significantly different sublimation pathways of the heterostructure. Furthermore, triangular and quasi hexagonal configurations are observed to be the dominant nanopores configurations in the heterostructure. Additionally, our DFT calculations provide a mechanistic understanding on the role of native defects and edge formation energies, revealing the antisite defects TeBi to be the dominant native defect in a Te-rich condition and playing a key role on the defects assisted sublimation. These findings significantly impact our understanding of controlling the nanoscale sublimation dynamics and can ultimately assist us in designing tunable low-dimensional chalcogenides.

The commercialization of perovskite solar cells is hindered by the poor long-term stability of the metal halide perovskite (MHP) light absorbing layer. Solution processing, the common fabrication method for MHPs, produces polycrystalline films with a wide variety of defects, such as point defects, surfaces, and grain boundaries. Although the optoelectronic effects of such defects have been widely studied, the evaluation of their impact on the long-term stability remains challenging. In particular, an understanding of the dynamics of degradation reactions at the atomistic scale is lacking. In this work, using reactive force field (ReaxFF) molecular dynamics simulations, we investigate the effects of defects, in the forms of surfaces, surface defects and grain boundaries, on the stability of the inorganic halide perovskite CsPbI$_{3}$. Our simulations establish a stability trend for a variety of surfaces, which correlates well with the occurrence of these surfaces in experiments. We find that a perovskite surface degrades by progressively changing the local geometry of PbI$_{\mathrm{x}}$ octahedra from corner- to edge- to face-sharing. Importantly, we find that Pb dangling bonds and the lack of steric hindrance of I species are two crucial factors that induce degradation reactions. Finally, we show that the stability of these surfaces can be modulated by adjusting their atomistic details, either by creating additional point defects or merging them to form grain boundaries. While in general additional defects, particularly when clustered, have a negative impact on the material stability, some grain boundaries have a stabilizing effect, primarily because of the additional steric hindrance., Comment: 29 pages, 7 figures

Ruthenium thin films can serve as protective caps for multi-layer extreme ultraviolet mirrors exposed to atomic hydrogen. Hydrogen permeation through ruthenium is problematic as it leads to blisters on the mirrors. H has been shown to exhibit low solubility in bulk Ru, and rapidly diffuses in and out of Ru. Therefore, the underlying mechanisms of the blistering effect remains unknown. This work makes use of reactive molecular dynamics simulations to study the influence of imperfections in a Ru film on the behaviour of H. For the Ru/H system, a ReaxFF force field was parametrised which reproduces structures and energies obtained from quantum-mechanical calculations. Molecular dynamics simulations have been performed with the newly-developed force field, to study the effect of tilt and twist grain boundaries on the overall diffusion behaviour of H in Ru. Our simulations show the tilt and twist grain boundaries provide energetically favourable sites for hydrogen atoms and act as sinks and highways for H. They therefore block H transport across their planes, and favour diffusion along their planes. This results in the accumulation of hydrogen at the grain boundaries. The strong effect of the grain boundaries on the hydrogen diffusion suggests tailoring the morphology of ruthenium thin films as a means to curb the rate of hydrogen permeation.

Chemically stable quantum-confined 2D metals are of interest in next-generation nanoscale quantum devices. Bottom-up design and synthesis of such metals could enable the creation of materials with tailored, on-demand, electronic and optical properties for applications that utilize tunable plasmonic coupling, optical non-linearity, epsilon-near-zero behavior, or wavelength-specific light trapping. In this work, we demonstrate that the electronic, superconducting and optical properties of air-stable two-dimensional metals can be controllably tuned by the formation of alloys. Environmentally robust large-area two-dimensional InxGa1-x alloys are synthesized by Confinement Heteroepitaxy (CHet). Near-complete solid solubility is achieved with no evidence of phase segregation, and the composition is tunable over the full range of x by changing the relative elemental composition of the precursor. The optical and electronic properties directly correlate with alloy composition, wherein the dielectric function, band structure, superconductivity, and charge transfer from the metal to graphene are all controlled by the indium/gallium ratio in the 2D metal layer.

Halide perovskites make efficient solar cells due to their exceptional optoelectronic properties, but suffer from several stability issues. The characterization of the degradation processes is challenging because of the limitations in the spatio-temporal resolution in experiments and the absence of efficient computational methods to study the reactive processes. Here, we present the first effort in developing reactive force fields for large scale molecular dynamics simulations of the phase instability and the defect-induced degradation reactions in inorganic CsPbI$_{3}$. We find that the phase transitions are driven by a combination of the anharmonicity of the perovskite lattice with the thermal entropy. At relatively low temperatures, the Cs cations tend to move away from the preferential positions with good contacts with the surrounding metal halide framework, potentially causing its conversion to a non-perovskite phase. Our simulations of defective structures reveal that, although both iodine vacancies and interstitials are very mobile in the perovskite lattice, the vacancies have a detrimental effect on the stability, initiating the decomposition reactions of perovskites to PbI$_{2}$. Our work puts ReaxFF forward as an effective computational framework to study reactive processes in halide perovskites., Comment: 11 pages, 6 figures

Transition metal impurities such as nickel, copper, and iron, in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon lattice. In this work, we performed ReaxFF reactive force field based molecular dynamics simulations, studying different configurations of individual and crystalline copper atoms inside bulk silicon by looking at the diffusional behavior of copper in silicon. The ReaxFF Cu/Si parameter set was developed by training against DFT data, including the energy barrier for an individual Cu-atom inside a silicon lattice. We found that the diffusion of copper atoms has a direct relationship with the temperature. Moreover, it is also shown that individual copper atoms start to clusterize inside bulk silicon at elevated temperatures. Our simulation results provide a comprehensive picture of the effects of temperature and copper concentration on the crystallization of individual copper inside silicon lattice. Finally, the stress-strain relationship of Cu/Si compounds under uniaxial tensile loading have been obtained. Our results indicate a decrease in the elastic modulus with increasing level of Cu-impurity concentration. We observe spontaneous microcracking of the Si during the stress-strain tests as a consequence of the formation of a small Cu clusters adjacent to the Si surface., Comment: 29 pages, 13 figures

Carbon Fibers (CFs) are the key solution for the future lightweight vehicle with enhanced fuel efficiency and reduced emissions owing to their ultrahigh strength to weight ratio. However, the high cost of the current dominant PAN-based CFs hinders their application. The use of low-cost alternative precursors may overcome this issue. Unfortunately, low-cost CFs derived from cheaper single component precursors suffer from poor mechanical properties. Developing composite CFs by adding nanoadditives is very promising for low-cost CFs. Therefore, a fundamental understanding of carbonization condition impacts and polymer/additives conversion mechanisms during whole CF production are essential to develop low-cost CFs. In this work, we have demonstrated how the carbonization temperature affects the PAN/graphene CFs properties by performing a series of ReaxFF based molecular dynamics simulations. We found that graphene edges along with the nitrogen and oxygen functional groups have a catalytic role and act as seeds for the graphitic structure growth. Our MD simulations unveil that the addition of the graphene to PAN precursor modifies all-carbon membered rings in CFs and enhances the alignments of 6-member carbon rings in carbonization which leads to superior mechanical properties compare to PAN-based CFs. These ReaxFF simulation results are validates by experimental structural and mechanical characterizations. Interestingly, mechanical characterizations indicate that PAN/graphene CFs carbonized at 1250 C demonstrate 90.9% increase in strength and 101.9% enhancement in Young's modulus compare to the PAN-based CFs carbonized at 1500 C. The superior mechanical properties of PAN/graphene CFs at lower carbonization temperatures offers a path to both energy savings and cost reduction by decreasing the carbonization temperature and could provide key insights for the development of low-cost CFs.

Engineering atomic-scale defects is crucial for realizing wafer-scale, single-crystalline transition metal dichalcogenide monolayers for electronic devices. However, connecting atomic-scale defects to larger morphologies poses a significant challenge. Using electron microscopy and atomistic simulations, we provide insights into WS2 crystal growth mechanisms, providing a direct link between synthetic conditions and the microstructure. Dark-field TEM imaging of coalesced monolayer WS2 films illuminates defect arrays that atomic-resolution STEM imaging identifies as translational grain boundaries. Imaging reveals the films to have nearly a single orientation with imperfectly stitched domains. Through atomic-resolution imaging and ReaxFF reactive force field-based molecular dynamics simulations, we observe two types of translational mismatch and discuss their atomic structures and origin. Our results indicate that the mismatch results from relatively fast growth rates. Through statistical analysis of >1300 facets, we demonstrate that the macrostructural features are constructed from nanometer-scale building blocks, describing the system across sub-{\AA}ngstrom to multi-micrometer length scales., Comment: 20 pages, 4 figures

Motivated by the problem of optimization of force-field systems in physics using large-scale computer simulations, we consider exploration of a deterministic complex multivariate response surface. The objective is to find input combinations that generate output close to some desired or "target" vector. In spite of reducing the problem to exploration of the input space with respect to a one-dimensional loss function, the search is nontrivial and challenging due to infeasible input combinations, high dimensionalities of the input and output space and multiple "desirable" regions in the input space and the difficulty of emulating the objective function well with a surrogate model. We propose an approach that is based on combining machine learning techniques with smart experimental design ideas to locate multiple good regions in the input space.

Human transport to Mars and deep space explorations demand the development of new materials with extraordinary high performance-to-mass ratios. Promising candidates to fulfill these requirements are ultrahigh strength lightweight (UHSL) materials, which consist of polymer matrices fortified by pristine carbon nanotubes (CNTs). Previous investigations have showed that with an increase in CNT diameter, its preferred configuration changes from a circular form to a flattened shape that can be obtained in high pressure or tension conditions. The ReaxFF reactive force field can reveal detailed chemical interactions at the atomistic scale. To enable ReaxFF simulations on CNT/polymer interfaces, we trained force field parameters to capture the proper structure of flattened carbon nanotubes (flCNTs), i.e. dumbbell-like shape CNTs, against available polymer consistent force field -- interface force field (PCFF-IFF) data which had good proximity to density functional theory (DFT) data. In this study we used accelerated ReaxFF molecular dynamics simulation using the optimized force field to study the polymerization of diglycidyl ether of bisphenol F (Bis F) and diethyltoluenediamine (DEDTA) molecules in vicinity of circular and flattened CNTs. Our results indicate that the flat regions of flCNT are more favorable spots for the polymers to settle compared to curved regions due to higher binding energies. Moreover, higher dimer generation around flCNT results in more effective coating of the CNT which leads to higher load transfer in compared to circular CNT. According to our results there is a high alignment between polymers and CNT surface which is due to strong pi-pi interactions of aromatic carbon rings in the polymers and CNTs. These ReaxFF simulations indicate the capability of this method to simultaneously observe the polymerization of monomers along with their interactions with CNTs., Comment: 18 Pages, 14 Figures

Ever since the discovery of graphene and subsequent explosion of interest in single atom thick materials, studying their mechanical properties has been an active area of research. New length scales often necessitate a rethinking of physical laws, making such studies crucial for understanding and ultimately utilizing novel material properties. Here we report on the investigation of nanoscale periodic ripples in suspended, single layer graphene sheets by scanning tunneling microscopy and atomistic scale simulations. Unlike the sinusoidal ripples found in classical fabrics, we find that graphene forms triangular ripples, where bending is limited to a narrow region on the order of a few unit cell dimensions at the apex of each ripple. This non-classical bending profile results in graphene behaving like a bizarre fabric, which regardless of how it is draped, always buckles at the same angle. Investigating the origin of such non-classical mechanical properties, we find that unlike a thin classical fabric, both in-plane and out-of-plane deformations occur in a graphene sheet. These two modes of deformation compete with each other, resulting in a strain-locked optimal buckling configuration when draped. Electronically, we see that this in-plane deformation generates pseudo electric fields creating a ~3 nm wide pnp heterojunction purely by strain modulation., Comment: 34 pages, 5 figures and supplementary section

Three-dimensional epitaxial heterostructures are based on covalently-bonded interfaces, whereas those from 2-dimensional (2D) materials exhibit van der Waals interactions. Under the right conditions, however, material structures with mixed interfacial van der Waals and covalent bonding may be realized. Atomically thin layers formed at the epitaxial graphene (EG)/silicon carbide (SiC) interface indicate that EG/SiC interfaces provide this unique environment and enable synthesis of a rich palette of 2D materials not accessible with traditional techniques. Here, we demonstrate a method termed confinement heteroepitaxy (CHet), to realize air-stable, structurally unique, crystalline 2D-Ga, In, and Sn at the EG/SiC interface. The first intercalant layer is covalently-bonded to the SiC, and is accompanied by a vertical bonding gradient that ends with van der Waals interactions. Such structures break out of plane centrosymmetry, thereby introducing atomically thin, non-centrosymmetric 2D allotropes of 3D materials as a foundation for tunable superconductivity, topological states, and plasmonic properties., Comment: 10 pages of main text with 3 figures, 21 pages of SI with 15 figures and 1 table

During the carbonization process of raw polymer precursors, graphitic structures can evolve. The presence of these graphitic structures affects mechanical properties of the carbonized carbon fibers. To gain a better understanding of the chemistry behind the evolution of these structures, we performed atomistic scale simulations using the ReaxFF reactive force field. Three different polymers were considered as a precursor: idealized ladder PAN (polyacrylonitrile), a proposed oxidized PAN and PBO (poly(p-phenylene-2,6-benzobisoxazole)). We determined the underlying molecular details of polymers conversion into a carbon fiber structure. Since these are C/H/O/N-based polymers, we first developed an improved force field for C/H/O/N chemistry based on the Density Functional Theory (DFT) data with a particular focus on N2 formation kinetics and its interactions with polymer-associated radicals formed during the carbonization process. Then, using this improved force field, we performed atomistic scale simulations of the initial stage of the carbonization process for the considered polymers. Based on our simulation data we determined the molecular pathways for the formation of low-molecular weight gas-species, all-carbon rings crucial for further graphitic structures evolution and possible alignment of the evolved all-carbon 6-membered rings clusters., Comment: 46 pages, 13 figures

Two-dimensional (2D) molybdenum disulfide (MoS2) has attracted significant attention because of its outstanding properties, suitable for application in several critical technologies like, solar cells, photocatalysis, lithium-ion batteries, nanoelectronics, and electrocatalysis. Similar to graphene and other 2D materials, the physical and chemical properties of MoS2 can be tuned by the chemical functionalization and defects. In this investigation, our objective is to explore the mechanical properties of single-layer MoS2 functionalized by the hydrogen atoms. We moreover analyze the effects of different types of defects on the mechanical response of MoS2 at the room temperature. To investigate these systems, we conducted reactive molecular dynamics simulations using the ReaxFF forcefield. We demonstrate that an increase in the hydrogen adatoms or defects contents significantly affects the critical mechanical characteristics of MoS2, elastic modulus, tensile strength, stretchability and failure behavior. Our reactive molecular dynamics results provide useful information concerning the mechanical response of hydrogenated and defective MoS2 and the design of nanodevices., Comment: Extreme Mechanics Letters 2018

Transition metal dichalcogenides (TMDs) represent a large family of high-quality 2D materials with attractive electronic, thermal, chemical, and mechanical properties. Chemical vapour deposition (CVD) technique is currently the most reliable route to synthesis few-atomic layer thick and large-scale TMDs films. However, the effects of grain boundaries formed during the CVD method on the properties of TMDs nanomembranes have remained less explored. In this study, we therefore aim to investigate the thermal conduction along polycrystalline molybdenum disulfide (MoS2) as the representative member of TMDs nanomembranes family. This goal was achieved by developing a combined atomistic-continuum multiscale method. In the proposed approach, reactive molecular dynamics simulations were carried out to assess thermal contact conductance of diverse grain boundaries with various defects configurations. The effective thermal conductivity along the CVD grown polycrystalline and single-layer MoS2 was finally acquired by conducting finite element modelling. Insight provided by this investigation can be useful to evaluate the effective thermal transport along a wide variety of 2D materials and structures.

Lithium-ion battery performance is strongly influenced by the ionic conductivity of the electrolyte, which depends on the speed at which Li ions migrate across the cell and relates to their solvation structure. The choice of solvent can greatly impact both solvation and diffusivity of Li ions. We use first principles molecular dynamics to examine the solvation and diffusion of Li ions in the bulk organic solvents ethylene carbonate (EC), ethyl methyl carbonate (EMC), and a mixture of EC/EMC. We find that Li ions are solvated by either carbonyl or ether oxygen atoms of the solvents and sometimes by the PF$_6^-$ anion. Li$^+$ prefers a tetrahedrally-coordinated first solvation shell regardless of which species are involved, with the specific preferred solvation structure dependent on the organic solvent. In addition, we calculate Li diffusion coefficients in each electrolyte, finding slightly larger diffusivities in the linear carbonate EMC compared to the cyclic carbonate EC. The magnitude of the diffusion coefficient correlates with the strength of Li$^+$ solvation. Corresponding analysis for the PF$_6^-$ anion shows greater diffusivity associated with a weakly-bound, poorly defined first solvation shell. These results may be used to aid in the design of new electrolytes to improve Li-ion battery performance., Comment: 31 pages, 12 figures, Published in the Journal of Physical Chemistry B

Proton transfer across single layer graphene is associated with large computed energy barriers and is therefore thought to be unfavorable at room temperature unless nanoscale holes or dopants are introduced, or a potential bias is applied. Here, we subject single layer graphene supported on fused silica to cycles of high and low pH and show that protons transfer reversibly from the aqueous phase through the graphene to the other side where they undergo acid-base chemistry with the silica hydroxyl groups. After ruling out diffusion through macroscopic pinholes, the protons are found to transfer through rare, naturally occurring atomic defects. Computer simulations reveal low energy barriers of 0.68 to 0.75 eV for aqueous proton transfer across hydroxyl-terminated atomic defects that participate in a Grotthuss-type relay, while pyrylium-like ether terminations shut down proton exchange. Unfavorable energy barriers to helium and hydrogen transfer indicate the transfer process is selective for aqueous protons., Comment: 80 pages, including Supporting Information, 3 Figures, 1 Table, final pre-edited version

Free standing silicene is predicted to display comparable electronic properties as graphene. However, the yet synthesized silicene-like structures have been only realized on different substrates which turned out to exhibit versatile crystallographic structures that are very different from the theoretically predicted buckled phase of freestanding silicene. This calls for a different approach where silicene is stabilized using very weakly interacting surfaces. We propose here a novel route by using graphene bilayer as a scaffold. The confinement between the flat graphene layers results in a planar clustering of Si atoms with small buckling, which is energetically unfavorable in vacuum. Buckled hexagonal arrangement of Si atoms similar to free-standing silicene is observed for large clusters, which, in contrast to Si atoms on metallic surfaces, is only very weakly van der Waals coupled to the graphene layers. These clusters are found to be stable well above room temperature. Our findings, which are supported by density functional tight-binding calculations, show that intercalating bilayer graphene with Si is a favorable route to realize silicene., Comment: to appear in Phys. Rev. B