banner

Blog

Jul 07, 2023

Ultrahigh resistance of hexagonal boron nitride to mineral scale formation

Nature Communications volume 13, Article number: 4523 (2022) Cite this article

5948 Accesses

6 Citations

10 Altmetric

Metrics details

Formation of mineral scale on a material surface has profound impact on a wide range of natural processes as well as industrial applications. However, how specific material surface characteristics affect the mineral-surface interactions and subsequent mineral scale formation is not well understood. Here we report the superior resistance of hexagonal boron nitride (hBN) to mineral scale formation compared to not only common metal and polymer surfaces but also the highly scaling-resistant graphene, making hBN possibly the most scaling resistant material reported to date. Experimental and simulation results reveal that this ultrahigh scaling-resistance is attributed to the combination of hBN’s atomically-smooth surface, in-plane atomic energy corrugation due to the polar boron-nitrogen bond, and the close match between its interatomic spacing and the size of water molecules. The latter two properties lead to strong polar interactions with water and hence the formation of a dense hydration layer, which strongly hinders the approach of mineral ions and crystals, decreasing both surface heterogeneous nucleation and crystal attachment.

Interfacial interactions play a fundamental role in many aqueous processes including adsorption, catalytic reaction, corrosion, filtration, and scale formation. In particular, scale formation, i.e., the development of mineral deposits on a material surface due to precipitation from the bulk solution and/or crystal formation initiated by surface nucleation, has great impacts on interfacial transfer of mass, heat, electrons, and light. It causes profound performance decline in numerous industrial processes, such as impaired heat transfer in heat exchangers and boilers, increased pressure drop in pipes, flow blockage in filtration membranes, corrosion damage of steam turbines, decreased conductivity and activity of electrodes, premature failure of heating and electrochemical components etc.1,2,3,4, all leading to higher operation cost and safety risk. It is reported that economical loss due to mineral scaling in boilers, turbines, and heat exchangers accounts for 0.17–0.25% of the gross domestic product (GDP) in industrialized nations5. Understanding mineral scaling behavior is important for the development of next generation materials and technologies that address these critical challenges.

Scale formation can occur through the deposition of mineral crystals formed in the bulk solution, as well as through surface-induced heterogeneous nucleation with crystals growing from nucleation sites on a surface6. Both processes are strongly influenced by material surface properties. Similar to deposition of other particles, material properties affect the attachment of mineral crystals via hydrophobic and electrostatic interactions. Surface-induced heterogeneous nucleation is a more thermodynamically favorable process, but it is poorly understood because it happens on very small time and length scales7. Few previous studies have investigated the different surface properties that influence surface-induced heterogeneous nucleation: roughness, charge, and hydrophobicity6. Surface roughness is directly related to the number of nucleation sites; it is generally recognized that mineral crystallization increases with surface roughness. However, findings on the impact of charge and surface hydrophobicity have been inconsistent. For example, some studies found that surface charge influenced the heterogeneous nucleation via electrostatic interactions or complexation reactions with the mineral ions2,8, while others reported similar nucleation rates on surfaces of different charges9. Contradictory results have also been reported on the role of surface hydrophobicity in mineral scaling. Hydrophilic coatings such as graphene oxide (GO), graft polymers, and polyethylene glycol, were shown to delay the onset of CaCO3 scaling in some studies10,11,12, while other studies showed that hydrophilic surfaces promoted CaCO3 nucleation8, and GO had no anti-scaling effect13. One possible reason for such apparent contradiction is that modification of a surface property (e.g., hydrophobicity or charge) often leads to inevitable changes in other surface properties, which makes it difficult to discern the role of an individual surface property for development of anti-scaling materials.

Two dimensional (2D) materials possess atomically smooth surfaces, and have attracted tremendous interest for their potential applications in processes where interfacial interactions play a critical role. For example, graphene was shown to effectively inhibit nucleation of metals14 and metal oxides15 in nonaqueous phases. Theoretical and experimental research on water-surface interactions16,17,18,19,20,21 and related phenomena, such as ultrafast water transport22,23,24,25,26, de-icing27, and anti-fouling28, also suggest that the atomically smooth morphology plays an important role in preventing adhesion. However, no studies have investigated the scaling behavior, especially heterogeneous nucleation, on the surface of 2D nanomaterials such as graphene and hexagonal boron nitride (hBN) in aqueous solutions. Furthermore, 2D materials vary widely in surface chemistry despite their common feature of atomic-level smoothness: Graphene consists of a single layer of carbon atoms arranged in a honeycomb lattice nanostructure, featuring a small lattice constant, low in-plane polarity, and high hydrophobicity; hBN, another important 2D material with its lattice structure and lattice constant similar to those of graphene29,30, has high in-plane polarity due to the boron-nitrogen bond and higher hydrophilicity than graphene18. It is unknown how such differences in surface chemistry influence mineral scale formation on these atomically smooth surfaces.

Here we investigate mineral scale formation on graphene and hBN surface, and compare it to that on metal (i.e., titanium (Ti)) and polymer (i.e., polyvinylidene fluoride (PVDF)) surfaces commonly used for high scaling potential applications (Supplementary Note 1, Supplementary Figs. 1 and 2). Surface induced heterogeneous nucleation on theses surfaces are investigated using both real time, in situ measurements as well as ex situ characterization methods. We also directly quantify the binding force of mineral crystals grown from surface-induced heterogenous nucleation. Experimental measurements combined with density functional theory (DFT) illustrate the effect of surface chemistry on the formation of the hydration layer and its key role in mineral ion-surface interactions. Very importantly, we discover hBN’s outstanding anti-scaling properties and demonstrate its potential application as an anti-scaling coating in practical engineering systems.

To evaluate the scaling behaviors of graphene and hBN, graphene and hBN nanocoatings grown on flat Cu substrate were dipped in supersaturated CaSO4 solution (50 mmol L−1, saturation index (SI) of 3.28) and characterized using a video camera for 220 min (Fig. 1A). As the induction time of homogeneous nucleation for CaSO4 at this concentration (<40 min6) is much shorter than the operation time, both bulk precipitation due to homogeneous nucleation and surface induced heterogeneous nucleation contribute to scale formation. As shown in Fig. 1B and Supplementary Movie 1, CaSO4 crystals were found on the pristine Cu surface after 20 min, and severe scale formation (0.123 ± 0.014 mg mm−2) was observed after 220 min of testing. On the graphene surface, CaSO4 crystals occurred after 40 min, and increased to 0.015 ± 0.008 mg mm−2 at the end of the experiment (Fig. 1C). However, the hBN-coated Cu surface exhibited almost no scale formation (0.001 ± 0.001 mg mm−2) during the 220 min experiment, except at its edges where the hBN nanocoating was damaged due to sample cutting (Fig. 1B).

A Experimental setup and B video snapshots of the pristine Cu (pristine), graphene (Gra) coated Cu (Gra@Cu), and hBN coated Cu (hBN@Cu) exposed to a supersaturated CaSO4 solution. The scale bars are 3 mm in length; C Amount of CaSO4 precipitate formed on the three samples after 220 min of exposure. Error bars represent the standard deviation of precipitation amount.

After illustrating the macroscopic anti-scaling behavior of graphene and hBN using video camera, we investigated the scale formation caused by surface-induced heterogeneous nucleation on graphene, hBN, Ti, and PVDF using microscopic ex situ and in situ measurements with supersaturated CaCO3. CaCO3 and CaSO4 are both common scalants in water, but CaCO3 has a solubility (0.12 mmol L−1) two orders of magnitude lower than CaSO4 (19.10 mmol L−1). It allows using much lower concentration to prepare supersaturated CaCO3 solution with saturation index similar to supersaturated CaSO4, which would generate smaller amount of scalants on the surfaces and is conducive to performing microscopic ex situ or in situ characterizations. As shown in Fig. 2A, supersaturated CaCO3 (SI of 1.18) was introduced to the test surfaces immediately after preparation before homogeneous nucleation occurred (Supplementary Note 2, Supplementary Fig. 3). Vaterite was the main CaCO3 crystal formed on the test surfaces (Supplementary Fig. 4). Its growth was the fastest on Ti (Fig. 2B), followed by PVDF (Fig. 2C), graphene (Fig. 2D), and hBN (Fig. 2E). The same trend was observed at different reaction time (Fig. 2F). In over 15 h, almost no vaterite was observed on hBN (0.4 ± 0.2 mm−2), while 55.6 ± 3.3, 13.2 ± 4.7, and 3.6 ± 0.8 vaterite crystals mm−2 were found on Ti, PVDF, and graphene surfaces, respectively (Fig. 2G). The very few vaterite crystals formed on hBN were much smaller in size than those formed on the other surfaces. These results suggest prolonged vaterite heterogeneous nucleation induction time and slower vaterite formation kinetics on the hBN surface.

A A schematic of QCMD experimental setup for characterizing CaCO3 heterogeneous nucleation on four surfaces. SI: saturation index, PC: personal computer, QCMD: quartz crystal microbalance with dissipation; SEM images of CaCO3 crystals formed on B Ti, C PVDF, D graphene, and E hBN coatings on Cu substrate after 9 h of operation, with insets showing the size distribution of vaterite crystals formed. The particle size distribution shown in the B–E inserts are obtained using the ImageJ software and fitted with a Gaussian distribution function; Average number and size of vaterite crystals grown on the four surfaces (F) as a function of operation time and G after 15 h of operation; H–J Real time in-situ monitoring of CaCO3 scale formation using QCMD. H Ion adsorption. I Induction and growth of CaCO3 crystals. Three arrows marked on the Ti, PVDF, and graphene curves indicate the start and finish of two stages: the induction stage i and the growth stage ii. Data obtained using the hBN sample only exhibit one stage (stage i) during the operation. J Average frequency change rate (ΔF/Δt) during the two-stage scale-formation process and surface roughness of the four nanocoatings. Error bars represent the standard deviation.

Real time, in situ measurement of CaCO3 nucleation and growth rates using the quartz crystal microbalance with dissipation (QCMD) technique (Fig. 2A) further confirms that hBN greatly hinders surface-induced heterogeneous nucleation of CaCO3 (Fig. 2H, I). Upon introduction of an influent solution supersaturated with CaCO3, adsorption of ions (e.g., Ca2+, CO32−) led to an immediate albeit small change in frequency and dissipation, which reached equilibrium in ~1 h. The larger frequency and dissipation changes observed for the Ti and PVDF surfaces suggested more adsorption of Ca2+ and CO32− ions than that on the graphene and hBN surfaces. The similar frequency and dissipation changes on graphene and hBN surface suggested that they had similar ion adsorption, a result of their similar surface charge31,32. After the initial ion adsorption, greater frequency and dissipation changes occurred, and exhibited two distinct stages. The first stage (stage i) featured relatively slow changes in frequency (dF/dt) and dissipation (dD/dt), signaling the formation of vaterite nuclei, i.e., the induction process. This was followed by a much faster frequency shift (stage ii), which signaled crystal growth from the nuclei previously formed on the surface (Fig. 2I). CaCO3 mass accumulation rate calculated using the Sauerbrey equation (Supplementary Fig. 5) showed a very large increase in stage ii compared to that in stage i on Ti, PVDF and graphene surfaces. As shown in Fig. 2I, J, the Ti surface showed the shortest induction time and the fastest crystal growth kinetics, followed by PVDF. The heterogeneous nucleation induction time on graphene was significantly longer compared to Ti and PVDF (Fig. 2I), and a much slower crystal growth kinetics was observed (Fig. 2J and Supplementary Fig. 5). On the hBN surface, the frequency and dissipation shifts throughout the experiment were very small. It did not exhibit a distinct stage ii with fast frequency and dissipation changes. This suggests that the induction time on the hBN surface is longer than the approximately 8 h contact time with the test solution in the experiment; no detectable crystal growth occurred during this period (Fig. 2I, J and Supplementary Fig. 5). This may be attributed to the low ion adsorption on hBN (Fig. 2I), which prolongs the induction of heterogeneous nucleation and hinders vaterite crystal growth. After 10 h of operation, the system was flushed with ultrapure water. Frequency and dissipation signals approached the original DI water baseline at different rate for all samples (Fig. 2I), confirming that the frequency and dissipation changes observed earlier resulted from the nucleation and growth of vaterite, which was removed or dissolved when flushed with ultrapure water.

Although the four material surfaces are all nominally flat and smooth at the macroscopic length scale, their nanoscale morphologies differ significantly (Fig. 3A). The root-mean-square roughness values of Ti (2.0 ± 0.2 nm) and PVDF (1.2 ± 0.2 nm) are almost an order of magnitude higher than those of graphene (0.3 ± 0.1 nm) and hBN (0.2 ± 0.1 nm), which are atomically smooth with no observable sags, crests, or other surface irregularities. During the nucleation and crystal growth stages, the nanoscale rough features on Ti and PVDF serve as heterogeneous nucleation sites, provide cratered surface for nuclei attachment, and increases friction that may hinder local flow velocity33, all conducive to vaterite nucleation and crystal growth, resulting in a shorter induction stage and faster vaterite growth (Fig. 2I, J).

AFM image of A Ti, B PVDF, C graphene, D hBN surfaces; E–H Measurement of vaterite detaching force from various coatings. E Sectional view schematic showing the measurement method by adding a lateral force on a vaterite crystal grown on a sample surface. Top view SEM image of a vaterite crystal (F) before and G after being pushed off from the sample surface by the nanoindenter tip. The scale bars in F and G represent 10 µm. The dashed circle in G shows the contact area between the vaterite crystal and the underlying sample surface after detachment; H Applied force load increases linearly with the moving depth of the nanotip upon contact to the vaterite crystal; I The average vaterite detaching force and surface roughness of the four nanocoatings. Error bars represent the standard deviation of corresponding results.

During the crystal growth stage, the higher roughness of Ti and PVDF surfaces increases the contact area with the vaterite crystals, resulting in stronger binding forces. We quantify the lateral force required to detach vaterite crystals from the surface using a nanoindenter inside a SEM (Supplementary Note 3, Supplementary Fig. 6). The applied load increases linearly with the lateral travelling distance of the nanoindenter tip relative to the vaterite crystal, until the vaterite crystal is pushed off the surface (Fig. 3E–H, and Supplementary Movies 2 and 3), i.e., the “detaching point”. The applied load at the detaching point is defined as the detaching force, which reflects the binding force between the vaterite crystal and the sample surface. The contact area (Fig. 3G) normalized detaching force followed the order of Ti (44.9 ± 12.7 μN μm−2) > PVDF (27.6 ± 3.4 μN μm−2) > hBN (12.1 ± 4.3 μN μm−2) ≈ graphene (11.1 ± 6.1 μN μm−2), correlating well with the surface roughness (Figs. 2J and 3I). The atomically smooth graphene and hBN exhibit notably lower binding forces. Note that these forces are much higher than typical colloidal adhesion forces34, a notable distinction between scaling due to crystal formation from surface-induced heterogeneous nucleation and that due to deposition of mineral crystals formed in the bulk solution.

The lower binding force on the atomically-smooth graphene and hBN surfaces measured in our study is consistent with findings from previous studies that show atomically smooth carbon and boron-nitride materials exhibit super-lubricity at solid/solid rigid junctions35,36,37,38 and peculiar interactions with water in 1D nanotubes22,39,40,41 or 2D structures21,23,24,42,43,44 at their solid/liquid interfaces. The reduced binding force allows easier detachment of the mineral crystals and decreases the formation of scales under hydraulic shear.

The similarity in their atomically smooth surface morphology is apparently responsible for the similarly in low vaterite detachment forces on graphene and hBN. The lower nucleation rate on hBN compared to graphene, however, suggests the difference in surface chemistry may play an important role. This is supported by DFT calculations. As a result of the delocalized π system in graphene, the charge distribution on graphene surface is much more uniform than on hBN, with a small dipole moment (0.036 D) and moderate carbon atom electronegativity (2.55) (Fig. 4A, B, Supplementary Fig. 7, Supplementary Note 4). In hBN, N atoms share their lone pair electrons with B, and the electrons of the π system are more localized on the N atoms. The electronegativity of N (3.04) is therefore notably higher than B (2.04), leading to corrugated charge distribution on hBN surface with a dipole moment of 0.047 D (Fig. 4A). The highest occupied molecular orbital (HOMO) is more localized on N atoms, and the lowest unoccupied molecular orbital (LUMO) is more concentrated on B atoms.

DFT simulated charge distribution on A hBN and B graphene surfaces. Blue and red spheres in A represent boron and nitrogen atoms, respectively. Silver, light blue, and light red spheres in B represent carbon atoms (Supplementary Fig. 7). The edge of both hBN and graphene flakes are passivated with H atoms; DFT simulated interaction of H2O, Ca2+, CO32−, and CaCO3 with C hBN and D graphene surface (Supplementary Fig. 11); E DFT calculated distance between water and surface, between water molecules, and between solutes (Ca2+, CO32−, or CaCO3) and hBN or graphene surfaces. Gra: graphene; F Minimum energy needed for a Ca2+ ion to penetrate the hydration layer on hBN and graphene surface calculated by DFT. Detailed calculation and simulation processes and results can be seen in Supplementary Note 4. G–I The formation of hydration layer on hBN and graphene surface characterized by QCMD. G QCMD experimental setup. PC: personal computer; H Frequency and dissipation change on hBN and graphene surfaces with changes in gas flow humidity; I The effect of gas flow humidity on average frequency decrease and the absolute ratio of dissipation change over frequency decrease (|ΔD/ΔF|).

The non-uniform charge distribution of the hBN surface has an important influence on its interaction with polar water molecules. The alternating negative and positive charge regions on hBN attract H and O atoms in water, respectively, leading to stronger interaction with water molecules than on graphene surface, where the charge is more uniform. As a result, water molecules are more closely packed and located closer to the hBN surface: the average distance between water molecules is 2.45 Å, and the closest distance between water molecule and the hBN surface is around 2.92 Å (Fig. 4C–E and Supplementary Fig. 8), consistent with previously reported values from quantum Monte Carlo simulations45. In comparison, on graphene surface, the closest water molecules are located 2.95 Å away from the surface, and the distance between water molecules measures 2.55 Å on average. Similar results have been reported in previous studies. Utilizing ab initio molecular dynamics simulation, it was reported that the friction coefficient on hBN was about 3 times larger than that on graphene, and the adsorbed water molecule had faster slippage on graphene than on hBN18.

The water molecule packing on the surface has direct impact on ion-surface interactions. Simulation of interactions between Ca2+, CO32−, or CaCO3 with hBN or graphene in the presence and absence of water (Supplementary Figs. 9–11) show that hBN binds more strongly with Ca2+, CO32−, and CaCO3 than graphene in the absence of water molecules (Supplementary Fig. 9), due to the locally charged B and N having higher attraction to the charged Ca2+, CO32−, and polar CaCO3 species. However, when water molecules are introduced, the distance between the surface and Ca2+ or CaCO3 is larger on hBN than on graphene (Fig. 4E and Supplementary Fig. 11). This is attributed to the denser hydration layer on hBN, which hinders the approach of Ca2+ and CaCO3 to the surface. More energy (1.0 eV) is required for a Ca2+ ion to penetrate the hydration layer to reach the hBN surface than the graphene surface (Fig. 4F).

Experimental measurements of water vapor adsorption confirm the simulation results. Adsorption of water vapor on graphene and hBN surfaces was measured using the QCMD technique at varying relative humidity of the influent N2/H2O mixed gas (Fig. 4G). Both graphene and hBN coated quartz sensors exhibit a decrease in the resonate frequency with the increase in the relative humidity of the influent gas (Fig. 4H), indicating increasing adsorption of water molecules on graphene and hBN surfaces. The hBN coated sensor experiences larger frequency decreases than the graphene coated sensor at all humidity levels (Fig. 4H, I), suggesting more adsorption of water on hBN than on graphene. Furthermore, the ratio of dissipation increase over frequency change (ΔD/ΔF) is smaller on hBN than on graphene (Fig. 4I). As water slippage on surfaces has little impact on ΔD and ΔF46,47, the measured ΔD/ΔF directly reflects the viscoelastic properties of the adsorbed water layer48. The smaller ΔD/ΔF on hBN indicates that the water layer adsorbed on the polar hBN surface is more compact and viscous compared to that on the uniformly charged graphene, consistent with the simulation results that a denser and more compact water layer forms on hBN. Although several studies have hinted the higher affinity of hBN for water18,21,41, our study experimentally confirm the denser and more rigid interfacial water layer structure on hBN than graphene.

Condensation of water molecules due to van der Waals attraction to form a hydration layer with higher density than bulk water has been discussed for materials of the carbon and boron-nitride family17,18,21, nanotubes39,40,41, and layered lamellar structures23,42,49,50. Our DFT simulation and QCMD water vapor adsorption experiments demonstrate the alignment of water molecules and formation of a condensed hydration layer on hBN. Interestingly, the calculated distance between water molecules on the 2D graphene and hBN surfaces (2.45–2.65 Å) are smaller than those reported for conventional material surfaces, such as NaCl(001) (~3.0 Å51), Cu(110) (2.8 ± 0.2 Å52), and Au(111) (2.71–2.88 Å53), as well as the calculated intermolecular distance between water molecules in bulk water (2.81 Å). The denser interfacial water layer formed on the 2D ultra-smooth graphene and hBN surfaces may be related to their small atom size and lattice constant that closely matches the size of water molecule. Most conventional materials have lattice constants larger than 3.0 Å (Supplementary Table 1), significantly larger than the water molecule size (2.7 Å), causing loose distribution of water on their surfaces. The lattice constants of graphene (2.46 Å) and hBN (2.50 Å) match the water molecule size better, with hBN’s lattice constant being the closest to a water molecule size. The close match between the distribution of hBN surface energy corrugation (charge nonuniformity) and the di-pole of water molecules leads to strong polar interactions between the hBN surface and polar water molecules, and hence compact packing of water molecules on hBN surface. Combining an atomically smooth surface, large in-plane polarity, and a lattice constant matching best with the size of water molecules (Supplementary Table 1), hBN is possibly one of the most anti-scaling materials. The understanding of how hBN interacts with water can also provide insights to other processes that occur in aqueous solutions, such as absorption, lubrication, catalysis, and corrosion etc.

Besides the flat Cu substrate, we grew hBN coating on both inside and outside surfaces of stainless steel (SS) pipes using the chemical vapor deposition (CVD) method (Fig. 5A), and assessed scale formation behaviors of the pristine and hBN-coated SS pipes by pumping real produced water collected from an oil & gas production site in Texas through the pipes (Fig. 5B, Supplementary Note 5). These experiments represent performance of the hBN coating on more realistic substrates in complex solutions. As shown in Fig. 5C, significant scale formation occurred inside the pristine SS tube after 7 days of operation, and its effective inner diameter decreased by ~17.2% after 14 days of operation, resulting in a 31.3% reduction of cross-sectional area for water passage. In contrast, the hBN-coated SS tube only experienced slight fouling by organic matter in the produced water at the outlet (Fig. 5D), with no measurable scale formation throughout the tube. These results demonstrate the excellent anti-scaling potential of the hBN nanocoating for real water and wastewater applications.

A Photograph of stainless steel (SS) tube before (top) and after (bottom) the growth of hBN nanocoating (hBN@SS); B A photograph of produced water taken from an oil & gas production site in Texas, USA; Scaling behavior of the C pristine and D hBN coated SS tube during 14 days of exposure to real produced water. The scale bars are 1 mm in length. Red dash circles indicate the original inner diameter of the pristine SS tube.

Overall, the unprecedented scaling resistance of hBN combined with its other unique properties, including superior mechanical strength54,55, high dielectric constant56, high thermal conductivity57, and chemical and thermal stability58,59,60,61, make hBN an excellent candidate for multi-functional coatings in various industrial processes (e.g. fluid transmission, heat exchange, membrane separation). These coatings can protect substrate materials from scaling or corrosion62, as well as ensure energetically efficient fluid flow, high heat transfer rate, and long material lifetime.

In summary, we discover that the unique combination of atomically smooth surface, high in-plane polarity, and proper interatomic spacing (lattice constant) makes hBN possibly the most anti-scaling material known, with notable advantages over graphene. The atomic smoothness of both hBN and graphene surfaces reduces the number of surface heterogeneous nucleation sites as well as the binding force between the mineral crystal and the surface, resulting in much greater resistance to mineral scaling compared to conventional scale-resistant materials. Interestingly, hBN exhibits even greater resistance to surface-induced heterogeneous nucleation than graphene, which stems from the in-plane dipole of the boron-nitrogen bond at a scale closely matching that of water molecules, leading to in-plane atomic energy corrugation that favors interaction with water. Direct experimental measurements and DFT calculations show that such chemical structure results in the formation of a dense hydration layer on the hBN surface, which strongly hinders mineral ions from approaching the surface and hence suppresses the heterogenous nucleation process. The hBN nanocoating grown on a stainless-steel tube exhibits outstanding anti-scaling properties in real oil & gas produced water, demonstrating its potential application in practical engineering systems. The results of the study provide important insights for future development of novel functional materials by manipulating their interactions with surrounding media. On the other hand, the scalability, durability, and long-term scaling resistance of the hBN coating as well as the specific roles of substrate materials and defects in large scale coatings need to be evaluated before practical applications can be possible.

Ti, PVDF, graphene and hBN coatings were formed on Cu foil (25 µm in thickness, McMaster-CARR, USA). Before coating, the Cu substrate was electrochemically polished. Ti films of 40 nm in thickness were formed on the Cu substrate using a sputter coater (Denton Desk V, Denton Vacuum, USA) at a current of 30 mA. PVDF was coated on the Cu substrate via spin coating. A 1wt % PVDF in DMF solution was applied at a rotation speed of 2000 rpm for 1 min, followed by overnight drying at 60 °C before use.

Graphene and hBN films were formed on Cu substrate by chemical vapor deposition (CVD). During growth, the electrochemically polished Cu foil was first loaded into a tube furnace and heated to 1000 °C, followed by annealing at 1 Torr 15% H2/Ar for 20 min. The graphene growth was then carried out by feeding 15% H2/Ar (flow rate of 100 sccm) and methane (flow rate of 10 sccm) to the furnace for 20 min. For hBN growth, ammonia borane was used as the precursor. After Cu foil annealing, ammonia borane was heated to ~85 °C for evaporation upstream of the tube furnace, and was carried by H2 gas to the Cu substrate to grow hBN at 1000 °C for 30 min. After growth with graphene or hBN, the Cu foil was rapidly cooled down to room temperature for further use.

After fabrication, the morphology of the four coatings were characterized by scanning electron microscopy (SEM, Quanta FEG 250, Thermo Fisher Scientific, USA) and atomic force microscopy (AFM) (NX20, Park, Suwon, Korea). The chemical property of the graphene and hBN were characterized by Raman spectroscopy (Renishaw inVia, UK) and X-ray photoelectron spectroscopy (XPS) (PHI Quantera II, Physical Electronics, USA) (Supplementary Note 1, Supplementary Figs. 1 and 2).

The graphene and hBN nanocoatings were formed on flat Cu substrates, and their scaling resistance was evaluated in a solution supersaturated with CaSO4, a common scalant in ground water that was very difficult to remove. The supersaturated CaSO4 solution (50 mmol L−1) was prepared by mixing 100 mmol L−1 CaCl2 with 100 mmol L−1 Na2SO4 at a 1:1 volume ratio. The prepared solution had a pH of 7.0 and a saturation index (SI) of 3.28. SI is defined as the ratio between the chemical activity product of the mineral ions and their solubility product63. During the experiments, the pristine Cu foil, hBN- and graphene- grown Cu foils were immersed in the CaSO4 solution. Scale formation was monitored using a video camera for 220 min (Supplementary Movie 1). The mass of mineral precipitation was determined by measuring the sample mass before and after the experiment.

Coated samples were cut into round wafers and installed in a flow cell. Supersaturated CaCO3 solution was prepared by mixing CaCl2 and Na2CO3 solutions, which were prepared using ultrapure water and aerated with air for 24 h to reach equilibrium before use. As shown in Fig. 2A, the equilibrated CaCl2 and Na2CO3 solutions were continuously fed into an inline mixer at 0.1 mL min−1 and then the flow cell at a combined flow rate of 0.2 mL min−1.The chemistry of the mixed feed solution is simulated using Visual MINTEQ (version 3.1, KTH): Ca2+ and CO32− activity are 1.23 × 10−4 and 1.24 × 10−4 mol L−1, respectively; pH is 8.84, ionic strength is 5.56 × 10−3, and calcite, aragonite, and vaterite saturation indexes ({Ca2+}{CO32−}/Ksp) are 4.45, 3.18, and 1.18 respectively. The induction time for homogeneous nucleation was determined by dynamic light scattering (DLS) measurements using a NanoBrook Omni (Brookhaven Instrument, Holtsville, NY, USA) (Supplementary Note 2, Supplementary Fig. 3). Since the hydraulic retention time from the mixer to the samples is much shorter than the induction time for CaCO3 homogeneous nucleation, formation of CaCO3 crystals on the sample surface is attributed to surface-induced heterogeneous nucleation. After the experiments, samples were retrieved and characterized using scanning electron microscopy (SEM, Quanta FEG 250, Thermo Fisher Scientific, USA) (Supplementary Note 2, Fig. 2B–E and Supplementary Fig. 4). The number and size of the crystals formed were analyzed using ImageJ software. More than 10 SEM images were analyzed for each surface to collect the crystal number and size data. The particle size distribution data were fitted with a Gaussian distribution function to obtain the average and standard deviation of particle size.

The heterogeneous nucleation and crystal formation of CaCO3 was also characterized in situ and in real time using the quartz crystal microbalance with dissipation (QCMD) technique (Qsense E4 analyzer, Biolin Scientific, Sweden), which has a mass sensitivity and a dissipation sensitivity of 1.8 ng cm−2 and 0.1 × 10−6 in liquid, respectively64. To perform the QCMD experiments, the four test materials were first coated on QCMD sensors. Ti and PVDF were coated on QCMD sensors using the same methods as those used for coating the Cu substrate. Graphene and hBN were first grown on Cu substrate, and then transferred onto QCMD sensors using a PMMA-assisted method. In this method, a PMMA film was spin-coated on the graphene- or hBN- grown on the Cu foil. The sample was then immersed in a FeCl3 solution to dissolve the Cu substrate. The graphene or hBN layer immobilized on the PMMA film was then transferred to DI water, and subsequently collected onto the QCMD sensor surface with the graphene or hBN side attaching to the sensor surface. Finally, acetone and IPA were used to dissolve the PMMA layer, exposing the graphene or hBN surface. The coated sensors were then mounted in 4 parallel QCMD cells, and characterized for CaCO3 nucleation and growth at a temperature of 22 °C using the same supersaturated solution and flow rates as the offline scaling experiments described above. The frequency and dissipation data were continuously monitored for each coated sensor.

The increase of adsorbed or precipitated mass (Δm, ng cm−2) on the sensors was calculated from the frequency change using the Sauerbrey equation:

Here C is the mass sensitivity constant, which is related to the properties of the quartz sensors and equals 17.7 ng cm−2 Hz−1 in this study. Δf is the frequency change (Hz). n is the harmonic number.

To measure the binding force between the vaterite crystals and the various surfaces that they grow from, samples were prepared by flowing the mixed CaCl2 and Na2CO3 solution over the four surfaces for 15 h to allow vaterite growth. After crystal growth, the samples are transferred into a SEM assembled with nanoindenter equipment. To perform the measurement, the tip of the nanoindenter approaches and pushes a vaterite crystal grown on sample surface at a constant speed of 0.03 μm s−1, until the vaterite crystal is detached from the underlying surface (Supplementary Note 3, Supplementary Fig. 6, Supplementary Movies 2 and 3). During this process, the applied load is recorded, and the applied load at the detaching point for the vaterite crystal is referred to as the detaching force. After normalizing the detaching force by contact area between the vaterite crystal and the surface, the obtained force reflects the binding force between the vaterite crystal and the underlying coating.

As shown in Supplementary Fig. 6 and Supplementary Movies 2 and 3, the nanotip can push the vaterite crystal at a direction perpendicular (Supplementary Fig. 6a, Supplementary Movie 2) or parallel (Supplementary Fig. 6b, Supplementary Movie 3) to the plane of the semi-hexagonal vaterite crystal. The detaching force differs greatly depending on the direction in which the load is applied (Supplementary Fig. 6a4, b4). This is due to the anisotropy structure of the vaterite crystal. As shown in Supplementary Fig. 6a, when pushing the crystal at the perpendicular direction, the detaching force is also greatly affected by the location of the contact point (e.g., distance from the coating surface). Therefore, data reported in the manuscript are all obtained by pushing the vaterite crystal in the direction parallel to the vaterite plane as that shown in Supplementary Fig. 6b.

In this study, we used density function theory (DFT) to evaluate (i) the charge distribution on graphene and hBN, (ii) the interaction between water molecules and hBN or graphene surface, (iii) the interaction of Ca2+, CO32−, and CaCO3 with bare hBN and graphene, and (iv) the interaction between Ca2+, CO32−, and CaCO3 and hydrated hBN or graphene. Detailed calculation and simulation methods can be found in the Supplementary Note 4 and Supplementary Figs. 7–11.

In addition to calculation, we also used the QCMD technique to quantitively verify the simulation results by measuring the adsorption of water vapor on graphene and hBN surfaces utilizing N2/H2O mixed gas at different humidity levels. In these experiments, a customized humidity control device was installed between a N2 gas tank and the QCMD analyzer (Fig. 4G). N2 gas with different humidity was continuously flown into the QCMD analyzer at a flow rate of 7.5 mL min−1. Corresponding frequency and dissipation change were recorded at 22 °C.

We also grew a hBN nanocoating inside a stainless-steel tube, and compared its scale forming behavior with an uncoated stainless-steel tube when exposed to a real oil & gas produced water taken from Texas. The produced water has a total organic carbon (TOC) concentration of 108 mg L−1 and conductivity of 148 mS cm−1. Detailed water quality data are shown in Supplementary Table 2 (Supplementary Note 5). As minerals have precipitated out from the solution during transportation of the produced water sample, the received produced water is no longer supersaturated (Supplementary Table 2) and cannot be directly utilized for scaling experiment. Instead, the oil & gas produced water is spiked with 50 mM CaCl2 or Na2SO4 to prepare the respective Ca2+- or SO42−- rich solutions. During the experiment, the two solutions were continuously fed into an inline mixer at 0.1 mL min−1 and then flowed through the pristine and hBN coated SS tubes at a combined flow rate of 0.2 mL min−1. Scale formation was monitored using a camera during 14 days of operation.

All data generated or analyzed during this study are reported in this published article and its supplementary information files, and are available from the author upon request.

All codes supporting the findings of this study are available from the corresponding author upon request.

Kazi, S. N., Duffy, G. G. & Chen, X. D. Mineral scale formation and mitigation on metals and a polymeric heat exchanger surface. Appl. Therm. Eng. 30, 2236–2242 (2010).

Article CAS Google Scholar

Mi, B. X. & Elimelech, M. Gypsum scaling and cleaning in forward osmosis: measurements and mechanisms. Environ. Sci. Technol. 44, 2022–2028 (2010).

Article ADS CAS PubMed Google Scholar

Tong, T. Z., Zhao, S., Boo, C., Hashmi, S. M. & Elimelech, M. Relating silica scaling in reverse osmosis to membrane surface properties. Environ. Sci. Technol. 51, 4396–4406 (2017).

Article ADS CAS PubMed Google Scholar

Zuo, K. C., Yuan, L. L., Wei, J. C., Liang, P. & Huang, X. Competitive migration behaviors of multiple ions and their impacts on ion-exchange resin packed microbial desalination cell. Bioresour. Technol. 146, 637–642 (2013).

Article CAS PubMed Google Scholar

Muller-Steinhagen, H., Malayeri, M. R. & Watkinson, A. P. Fouling of heat exchangers-new approaches to solve an old problem. Heat. Transf. Eng. 26, 1–4 (2005).

Article ADS CAS Google Scholar

Huang, X., Li, C., Zuo, K. & Li, Q. Predominant effect of material surface hydrophobicity on gypsum scale formation. Environ. Sci. Technol. 54, 15395–15404 (2020).

Article ADS CAS PubMed Google Scholar

Nielsen, M. H., Aloni, S. & De Yoreo, J. J. In situ TEM imaging of CaCO3 nucleation reveals coexistence of direct and indirect pathways. Science 345, 1158–1162 (2014).

Article ADS CAS PubMed Google Scholar

Aizenberg, J., Black, A. J. & Whitesides, G. M. Control of crystal nucleation by patterned self-assembled monolayers. Nature 398, 495–498 (1999).

Article ADS CAS Google Scholar

Shaffer, D. L., Tousley, M. E. & Elimelech, M. Influence of polyamide membrane surface chemistry on gypsum scaling behavior. J. Membr. Sci. 525, 249–256 (2017).

Article CAS Google Scholar

Cao, B., Ansari, A., Yi, X., Rodrigues, D. F. & Hu, Y. Gypsum scale formation on graphene oxide modified reverse osmosis membrane. J. Membr. Sci. 552, 132–143 (2018).

Article CAS Google Scholar

Kim, M.-M., Lin, N. H., Lewis, G. T. & Cohen, Y. Surface nano-structuring of reverse osmosis membranes via atmospheric pressure plasma-induced graft polymerization for reduction of mineral scaling propensity. J. Membr. Sci. 354, 142–149 (2010).

Article CAS Google Scholar

Ray, J. R., Wong, W. & Jun, Y.-S. Antiscaling efficacy of CaCO3 and CaSO4 on polyethylene glycol (PEG)-modified reverse osmosis membranes in the presence of humic acid: interplay of membrane surface properties and water chemistry. Phys. Chem. Chem. Phys. 19, 5647–5657 (2017).

Article CAS PubMed Google Scholar

Ray, J. R. et al. Hydrophilic, bactericidal nanoheater-enabled reverse osmosis membranes to improve fouling resistance. ACS Appl. Mater. Interfaces 7, 11117–11126 (2015).

Article CAS PubMed Google Scholar

Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614–617 (2013).

Article ADS CAS PubMed Google Scholar

Wang, X. R., Tabakman, S. M. & Dai, H. J. Atomic layer deposition of metal oxides on pristine and functionalized graphene. J. Am. Chem. Soc. 130, 8152–8153 (2008).

Article CAS PubMed Google Scholar

Briscoe, W. H. et al. Boundary lubrication under water. Nature 444, 191–194 (2006).

Article ADS CAS PubMed Google Scholar

Rafiee, J. et al. Wetting transparency of graphene. Nat. Mater. 11, 217–222 (2012).

Article ADS CAS PubMed Google Scholar

Tocci, G., Joly, L. & Michaelides, A. Friction of water on graphene and hexagonal boron nitride from ab initio methods: very different slippage despite very similar interface structures. Nano Lett. 14, 6872–6877 (2014).

Article ADS CAS PubMed Google Scholar

Li, H. & Zeng, X. C. Wetting and interfacial properties of water nanodroplets in contact with graphene and monolayer boron-nitride sheets. ACS Nano 6, 2401–2409 (2012).

Article CAS PubMed Google Scholar

Gopinadhan, K. et al. Complete steric exclusion of ions and proton transport through confined monolayer water. Science 363, 145–147 (2019).

Article ADS CAS PubMed Google Scholar

Grosjean, B., Bocquet, M. L. & Vuilleumier, R. Versatile electrification of two-dimensional nanomaterials in water. Nat. Commun. 10, 1656 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Holt, J. K. et al. Fast mass transport through sub-2-nanometer carbon nanotubes. Science 312, 1034–1037 (2006).

Article ADS CAS PubMed Google Scholar

Joshi, R. K. et al. Precise and ultrafast molecular sieving through graphene oxide membranes. Science 343, 752–754 (2014).

Article ADS CAS PubMed Google Scholar

Nair, R. R., Wu, H. A., Jayaram, P. N., Grigorieva, I. V. & Geim, A. K. Unimpeded permeation of water through helium-leak-tight graphene-based membranes. Science 335, 442–444 (2012).

Article ADS CAS PubMed Google Scholar

Esfandiar, A. et al. Size effect in ion transport through angstrom-scale silts. Science 358, 511–513 (2017).

Article ADS CAS PubMed Google Scholar

Radha, B. et al. Molecular transport through capillaries made with atomic-scale precision. Nature 538, 222–225 (2016).

Article ADS CAS PubMed Google Scholar

Li, Y. L. et al. Laser-induced graphene in controlled atmospheres: from superhydrophilic to superhydrophobic surfaces. Adv. Mater. 29, 1700496 (2017).

Article CAS Google Scholar

Seo, D. H. et al. Anti-fouling graphene-based membranes for effective water desalination. Nat. Commun. 9, 683 (2018).

Article ADS PubMed PubMed Central CAS Google Scholar

Dean, C. R. et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 5, 722–725 (2010).

Article ADS CAS PubMed Google Scholar

Song, L. et al. Large scale growth and characterization of atomic hexagonal boron nitride layers. Nano Lett. 10, 3209–3215 (2010).

Article ADS CAS PubMed Google Scholar

Qu, J., Li, Q., Luo, C., Cheng, J. & Hou, X. Characterization of flake boron nitride prepared from the low temperature combustion synthesized precursor and its application for dye adsorption. Coatings 8, 214 (2018).

Article CAS Google Scholar

Li, D., Müller, M. B., Gilje, S., Kaner, R. B. & Wallace, G. G. Processable aqueous dispersions of graphene nanosheets. Nat. Nanotechnol. 3, 101–105 (2008).

Article ADS CAS PubMed Google Scholar

Jia, J., Song, Q., Liu, Z. & Wang, B. Effect of wall roughness on performance of microchannel applied in microfluidic device. Microsyst. Technol. 25, 2385–2397 (2019).

Article CAS Google Scholar

Hong, Y.-K. et al. The effect of additives in post-Cu CMP cleaning on particle adhesion and removal. J. Electrochem. Soc. 151, G756 (2004).

Article CAS Google Scholar

Liu, S. W. et al. Robust microscale superlubricity under high contact pressure enabled by graphene-coated microsphere. Nat. Commun. 8, 14029 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Choi, J. S. et al. Friction anisotropy–driven domain imaging on exfoliated monolayer graphene. Science 333, 607–610 (2011).

Article ADS CAS PubMed Google Scholar

Lee, C. et al. Frictional characteristics of atomically thin sheets. Science 328, 76–80 (2010).

Article ADS CAS PubMed Google Scholar

Song, Y. M. et al. Robust microscale superlubricity in graphite/hexagonal boron nitride layered heterojunctions. Nat. Mater. 17, 894–899 (2018).

Article ADS CAS PubMed Google Scholar

Secchi, E. et al. Massive radius-dependent flow slippage in carbon nanotubes. Nature 537, 210–213 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Siria, A. et al. Giant osmotic energy conversion measured in a single transmembrane boron nitride nanotube. Nature 494, 455–458 (2013).

Article ADS CAS PubMed Google Scholar

Won, C. Y. & Aluru, N. R. Water permeation through a subnanometer boron nitride nanotube. J. Am. Chem. Soc. 129, 2748–2749 (2007).

Article CAS PubMed Google Scholar

Chen, C. et al. Functionalized boron nitride membranes with ultrafast solvent transport performance for molecular separation. Nat. Commun. 9, 1902 (2018).

Article ADS PubMed PubMed Central CAS Google Scholar

Mouterde, T. et al. Molecular streaming and its voltage control in angstrom-scale channels. Nature 567, 87–90 (2019).

Article ADS CAS PubMed Google Scholar

Comtet, J. et al. Direct observation of water-mediated single-proton transport between hBN surface defects. Nat. Nanotechnol. 15, 598–604 (2020).

Article ADS CAS PubMed Google Scholar

Wu, Y., Wagner, L. K. & Aluru, N. R. Hexagonal boron nitride and water interaction parameters. J. Chem. Phys. 144, 164118 (2016).

Article ADS PubMed CAS Google Scholar

Huang, K. & Szlufarska, I. Friction and slip at the solid/liquid interface in vibrational systems. Langmuir 28, 17302–17312 (2012).

Article CAS PubMed Google Scholar

Johannsmann, D., Reviakine, I. & Richter, R. P. Dissipation in films of adsorbed nanospheres studied by quartz crystal microbalance (QCM). Anal. Chem. 81, 8167–8176 (2009).

Article CAS PubMed Google Scholar

Cho, N.-J., Frank, C. W., Kasemo, B. & Höök, F. Quartz crystal microbalance with dissipation monitoring of supported lipid bilayers on various substrates. Nat. Protoc. 5, 1096–1106 (2010).

Article CAS PubMed Google Scholar

Fumagalli, L. et al. Anomalously low dielectric constant of confined water. Science 360, 1339–1342 (2018).

Article ADS CAS PubMed Google Scholar

Algara-Siller, G. et al. Square ice in graphene nanocapillaries. Nature 519, 443–445 (2015).

Article ADS CAS PubMed Google Scholar

Guo, J. et al. Real-space imaging of interfacial water with submolecular resolution. Nat. Mater. 13, 184–189 (2014).

Article ADS CAS PubMed Google Scholar

Shiotari, A. & Sugimoto, Y. Ultrahigh-resolution imaging of water networks by atomic force microscopy. Nat. Commun. 8, 14313 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Ma, R. Z. et al. Atomic imaging of the edge structure and growth of a two-dimensional hexagonal ice. Nature 577, 60–63 (2020).

Article ADS CAS PubMed Google Scholar

Falin, A. et al. Mechanical properties of atomically thin boron nitride and the role of interlayer interactions. Nat. Commun. 8, 1–9 (2017).

Article ADS Google Scholar

Yang, Y. et al. Intrinsic toughening and stable crack propagation in hexagonal boron nitride. Nature 594, 57–61 (2021).

Article ADS CAS PubMed Google Scholar

Watanabe, K., Taniguchi, T. & Kanda, H. Direct-bandgap properties and evidence for ultraviolet lasing of hexagonal boron nitride single crystal. Nat. Mater. 3, 404–409 (2004).

Article ADS CAS PubMed Google Scholar

Jiang, P., Qian, X., Yang, R. & Lindsay, L. Anisotropic thermal transport in bulk hexagonal boron nitride. Phys. Rev. Mater. 2, 064005 (2018).

Article CAS Google Scholar

Chilkoor, G. et al. Hexagonal boron nitride: the thinnest insulating barrier to microbial corrosion. ACS Nano 12, 2242–2252 (2018).

Article CAS PubMed Google Scholar

Hu, S. et al. Proton transport through one-atom-thick crystals. Nature 516, 227–230 (2014).

Article ADS CAS PubMed Google Scholar

Liu, Z. et al. Ultrathin high-temperature oxidation-resistant coatings of hexagonal boron nitride. Nat. Commun. 4, 2541 (2013).

Article ADS PubMed CAS Google Scholar

Zuo, K. et al. Multifunctional nanocoated membranes for high-rate electrothermal desalination of hypersaline waters. Nat. Nanotechnol. 15, 1025–1032 (2020).

Article ADS CAS PubMed Google Scholar

Husain, E. et al. Marine corrosion protective coatings of hexagonal boron nitride thin films on stainless steel. ACS Appl. Mater. Interfaces 5, 4129–4135 (2013).

Article CAS PubMed Google Scholar

Sand, K. K. et al. Calcite growth kinetics: dependence on saturation Index, Ca2+:CO32– activity ratio, and surface atomic structure. Cryst. Growth Des. 16, 3602–3612 (2016).

Article CAS Google Scholar

Uddenberg, H. Q-Sense E4 operator manual. Sweden: Q-Sense AB; 2009. 57.

Download references

This work was supported by the NSF Nanosystems Engineering Research Center for Nanotechnology-Enabled Water Treatment (EEC-1449500) (K.Z., X.H. and Q.L.), and the NSF I/UCRC Center for Atomically Thin Multifunctional Coatings (ATOMIC) under award # IIP-1539999 (T.Z., H.G. and J.L.). W.W. thank the financial support from NSFC under grant number 51788104 and 52001183. E.F.O. thank the Brazilian agencies CNPq and FAPESP (Grants 2013/08293-7, 2016/18499-0, and 2019/07157-9) for financial support, computational support from the Center for Computational Engineering and Sciences at Unicamp through the FAPESP/CEPID Grant No. 2013/08293-7 and the Center for Scientific Computing (NCC/GridUNESP) of São Paulo State University (UNESP).

These authors contributed equally: Kuichang Zuo, Xiang Zhang.

The Key Laboratory of Water and Sediment Sciences, Ministry of Education; College of Environment Sciences and Engineering, Peking University, Beijing, 100871, China

Kuichang Zuo

Department of Civil and Environmental Engineering, Rice University, MS 519, 6100 Main Street, Houston, TX, 77005, USA

Kuichang Zuo, Xiaochuan Huang, Pedro J. J. Alvarez & Qilin Li

NSF Nanosystems Engineering Research Center Nanotechnology-Enabled Water Treatment, Rice University, MS 6398, 6100 Main Street, Houston, TX, 77005, USA

Kuichang Zuo, Xiang Zhang, Xiaochuan Huang, Pedro J. J. Alvarez, Menachem Elimelech, Pulickel M. Ajayan, Jun Lou & Qilin Li

Department of Materials Science and NanoEngineering, Rice University, 6100 Main Street, Houston, TX, 77005, USA

Xiang Zhang, Eliezer F. Oliveira, Hua Guo, Tianshu Zhai, Pulickel M. Ajayan, Jun Lou & Qilin Li

São Paulo State Department of Education, São Paulo, Brazil

Eliezer F. Oliveira

Key Laboratory of Advanced Materials (MOE), School of Materials Science and Engineering, Tsinghua University, Beijing, 100084, PR China

Weipeng Wang

Department of Chemical and Environmental Engineering, Yale University, New Haven, CT, 06520-8286, USA

Menachem Elimelech

Department of Chemistry, Rice University, 6100 Main Street, Houston, TX, 77005, USA

Jun Lou

Department of Chemical and Biomolecular Engineering, Rice University, 6100 Main Street, Houston, TX, 77005, USA

Qilin Li

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

All authors contributed intellectually to this paper. K.Z., X.H. designed and performed the experiments, analyzed data, and drafted the manuscript. Q.L., J.L., P.M.A., W.W., M.E. and P.J.J.A. conceived the idea, designed the experiments, revised the manuscript, and led the project. X.Z., T.Z., H.G. assisted sample growth and characterization. E.F.O. performed the simulation.

Correspondence to Weipeng Wang, Pulickel M. Ajayan, Jun Lou or Qilin Li.

The authors declare no competing interests.

Nature Communications thanks Jean Comtet, Elisabeth Mansfield and Wenrong Yang for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Zuo, K., Zhang, X., Huang, X. et al. Ultrahigh resistance of hexagonal boron nitride to mineral scale formation. Nat Commun 13, 4523 (2022). https://doi.org/10.1038/s41467-022-32193-4

Download citation

Received: 27 January 2022

Accepted: 20 July 2022

Published: 04 August 2022

DOI: https://doi.org/10.1038/s41467-022-32193-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

SHARE