# Fileset

[s41563-023-01658-2.pdf](https://mdr.nims.go.jp/filesets/15238101-68c0-4337-aa94-ad651cacd0d7/download)

## Creator

Nathan Ronceray, Yi You, Evgenii Glushkov, Martina Lihter, Benjamin Rehl, Tzu-Heng Chen, Gwang-Hyeon Nam, Fanny Borza, [Kenji Watanabe](https://orcid.org/0000-0003-3701-8119), [Takashi Taniguchi](https://orcid.org/0000-0002-1467-3105), Sylvie Roke, Ashok Keerthi, Jean Comtet, Boya Radha, Aleksandra Radenovic

## Rights

[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[Liquid-activated quantum emission from pristine hexagonal boron nitride for nanofluidic sensing](https://mdr.nims.go.jp/datasets/fc9c18f8-425d-4234-9d4a-059360fce920)

## Fulltext

Liquid-activated quantum emission from pristine hexagonal boron nitride for nanofluidic sensingNature Materials | Volume 22 | October 2023 | 1236–1242 1236nature materialsArticle https://doi.org/10.1038/s41563-023-01658-2Liquid-activated quantum emission from pristine hexagonal boron nitride for nanofluidic sensingNathan Ronceray    1,2  , Yi You3,4, Evgenii Glushkov    1, Martina Lihter    1,9, Benjamin Rehl2, Tzu-Heng Chen1, Gwang-Hyeon Nam    3,4, Fanny Borza1, Kenji Watanabe    5, Takashi Taniguchi    6, Sylvie Roke    2, Ashok Keerthi    4,7, Jean Comtet    8, Boya Radha    3,4   & Aleksandra Radenovic    1 Liquids confined down to the atomic scale can show radically new properties. However, only indirect and ensemble measurements operate in such extreme confinement, calling for novel optical approaches that enable direct imaging at the molecular level. Here we harness fluorescence originating from single-photon emitters at the surface of hexagonal boron nitride for molecular imaging and sensing in nanometrically confined liquids. The emission originates from the chemisorption of organic solvent molecules onto native surface defects, revealing single-molecule dynamics at the interface through the spatially correlated activation of neighbouring defects. Emitter spectra further offer a direct readout of the local dielectric properties, unveiling increasing dielectric order under nanometre-scale confinement. Liquid-activated native hexagonal boron nitride defects bridge the gap between solid-state nanophotonics and nanofluidics, opening new avenues for nanoscale sensing and optofluidics.Nanostructures made of two-dimensional (2D) materials have become prominent in nanofluidic research1,2. Liquid confinement to a few molecular  layers between atomically smooth walls has led to anomalies in molecular transport3–6 and structure7,8. Nevertheless, a direct observation of these emerging phenomena remains challenging due to the limitations of current techniques in extreme confinements, where even molecular fluorophores cannot penetrate9. This calls for the development of imag-ing methods that can access molecular properties in confinement1,2.To tackle this goal, solid-state optically active defects show promise.  Fluorescent defects in diamond have enabled the optical probing of nanoscale matter10, including liquids11, but they cannot be easily embed-ded in 2D nanostructures. Coincidentally, hexagonal boron nitride (hBN) has been used in both nanofluidics6,8,12 and nanophotonics, where various point defects within its 6 eV bandgap have been identified as room-temperature quantum emitters13–15. Although emitters in hBN have been artificially induced using techniques such as irradiation16,17 or carbon doping18,19, the potential of liquid treatments remains largely unexplored. Recent studies have combined liquid and irradiation treat-ments to activate plasma-induced surface defects in hBN using water20 and binary mixtures of water with organic solvents21. Post-treatment Received: 9 May 2022Accepted: 31 July 2023Published online: 31 August 2023 Check for updates1Laboratory of Nanoscale Biology, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland. 2Laboratory for Fundamental BioPhotonics, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland. 3Department of Physics and Astronomy, School of Natural Sciences, The University of Manchester, Manchester, UK. 4National Graphene Institute, The University of Manchester, Manchester, UK. 5Research Center for Electronic and Optical Materials, National Institute for Materials Science, Tsukuba, Japan. 6Research Center for Materials Nanoarchitectonics, National Institute for Materials Science, Tsukuba, Japan. 7Department of Chemistry, School of Natural Sciences, The University of Manchester, Manchester, UK. 8Soft Matter Sciences and Engineering, ESPCI Paris, PSL University, CNRS, Sorbonne Université, Paris, France. 9Present address: Institute of Physics, Zagreb, Croatia.   e-mail: nathan.ronceray@epfl.ch; radha.boya@manchester.ac.uk; aleksandra.radenovic@epfl.chhttp://www.nature.com/naturematerialshttps://doi.org/10.1038/s41563-023-01658-2http://orcid.org/0000-0001-8069-2411http://orcid.org/0000-0001-7715-4524http://orcid.org/0000-0003-1859-8453http://orcid.org/0009-0000-5041-2917http://orcid.org/0000-0003-3701-8119http://orcid.org/0000-0002-1467-3105http://orcid.org/0000-0002-6062-7871http://orcid.org/0000-0002-8479-4762http://orcid.org/0000-0003-2389-3879http://orcid.org/0000-0003-1345-7029http://orcid.org/0000-0001-8194-2785http://crossmark.crossref.org/dialog/?doi=10.1038/s41563-023-01658-2&domain=pdfmailto:nathan.ronceray@epfl.chmailto:radha.boya@manchester.ac.ukmailto:aleksandra.radenovic@epfl.chNature Materials | Volume 22 | October 2023 | 1236–1242 1237Article https://doi.org/10.1038/s41563-023-01658-2(Fig. 1b). The fluorescence intensity gradually decreased over continu-ous illumination and stabilized after several seconds, revealing emis-sion from sub-diffraction spots (Fig. 1c and Supplementary Video 1), which can be localized with ~10 nm precision using single-molecule localization microscopy29 (Methods). We attribute this emission to the activation of defects present in the as-exfoliated crystal through con-tact with the liquid, resulting in randomly distributed transient emit-ters on the surface (Supplementary Fig. 1). Under constant 3.5 kW cm–2 illumination, the number of emitters decreased to a stable value of approximately 0.4 per square micron (Fig. 1d). This process does not deteriorate the crystal (Supplementary Figs. 2 and 3), and the decrease observed in Fig. 1c is reversible: when left in the dark, the crystal fluores-cence recovered within tens of minutes (Supplementary Fig. 4) without inducing additional emitters in the steady state (Supplementary Fig. 5).Remarkably, the fluorescent activation of the surface occurred with the most common organic solvents, including n-alkanes (pentane  to hexadecane) and primary alcohols (methanol to 1-pentanol). However, we did not observe any emission in pure water, heavy water and hydrogen peroxide. To quantitatively compare the steady- state fluorescence in different liquid media, we imaged freshly cleaved  hBN crystals in several liquids under 3.5 kW cm–2 illumination. The  steady-state fluorescence can be quantified using crystal brightness Icrystal, defined as the sum of the localized emitter intensities per surface unit and time unit (Fig. 1e). From these observations, we can classify  solvents into three types based on the extent of hBN fluorescent acti-vation. Most organic solvents, such as primary alcohols, n-alkanes and chloroalkanes, exhibited intense fluorescence (Icrystal > 500 pho-tons μm–2 s–1), representing type-I activation. Glycerol and other high-boiling-point liquids (≥200 °C) exhibited a limited but measurable level of fluorescence, classified as type II. On the other hand, pure water showed no activation, falling into type III (Icrystal < 10 photons μm–2 s–1). Furthermore, the addition of 10 v/v% of water in ethanol reduced the of ion-beam-exposed hBN with liquids has also been shown to modify the defect emission properties22. Furthermore, plasma-induced surface defects have been utilized to study single interfacial charge dynamics, such as proton hopping20,21. However, even mild plasma treatment induces mechanical and chemical changes that result in hBN crystals no longer having atomically smooth surfaces23, thus preventing the integration of defects in ultraflat van der Waals heterostructures, which are crucial for advancements in ångström-scale fluidics4,8.In this study, we demonstrate that organic solvents can activate visible-range quantum emission from pristine high-quality hBN crystals.  We attribute this phenomenon to the interaction between organic molecules and native surface defects24–27. By employing spectral super-resolution microscopy28, we observe defect-mediated molecular  random walks and couplings between defect dipoles and the liquid medium, leading to tunable emission wavelengths through the dielec-tric properties of the liquid. Leveraging the intrinsic properties of widely used 2D materials and common solvents, the fluorescence activation mechanism reported here is utilized to image nanofluidic structures, with emitters serving as the nanoscale probes of the order and dynamics of liquid media confined to the nanoscale.Liquid-activated fluorescence from pristine hBNThe surface of untreated hBN crystals exhibits visible-range fluores-cence when in contact with common organic solvents like ethanol. To demonstrate this effect, we exfoliated high-quality hBN crystals24 onto a glass coverslip, which was placed in a liquid-filled chamber on an inverted microscope (Fig. 1a). Photoluminescence (PL) of the crystals under 561 nm wide-field laser illumination (0.35–3.50 kW cm–2) was  collected using a high-numerical-aperture objective and projected  onto a camera chip. Pristine hBN in air or water did not exhibit fluo-rescence under these illumination conditions. However, we observed intense fluorescence from as-exfoliated crystals in contact with ethanol High-numerical-aperture objectiveOrganic solventhBN crystalsCoverslipadLaser ON +10 sSteady stateType IType IIWater Ethanolbec101102103I crystal (photons µm–2 s–1)EthanolChloroform1-butanolDimethyl sulfoxideCyclohexaneDichloromethaneGlycerolWaterEthylene  glycol0 5 10Time (s)0.51.01.5Emitters per µm2Type III5 µmFig. 1 | Liquid-induced fluorescence from pristine hBN crystals. a, Sketch of the experimental setup. b, Wide-field fluorescence images of an hBN crystal under 3.5 kW cm–2 and 561 nm laser light illumination with 1 s exposure time. No fluorescence was observed in water, but in ethanol, the entire crystal surface became fluorescent. The images underwent linear contrast enhancement.  c, A zoomed-in view of the dashed yellow box in b reveals dense clusters of emission when the laser is turned ON, with 6 ms exposure time. After 10 s of  wide-field illumination, the crystal surface reached a stable number of diffraction-limited isolated emitters. d, Localization microscopy-based counting of the emitters as a function of illumination time: after 5 s, a steady state was reached. The dashed line is a fit to an offset exponential relaxation. e, Liquid dependency of the crystal fluorescence, showing strongly activating liquids (type I), mildly activating liquids (type II) and no activation in water (type III).http://www.nature.com/naturematerialsNature Materials | Volume 22 | October 2023 | 1236–1242 1238Article https://doi.org/10.1038/s41563-023-01658-2number of emitters, resulting in a fivefold reduction in the fluorescence signal (Supplementary Fig. 6). The distinction between activation types could not be solely explained by the physical parameters of the solvents (Supplementary Fig. 7), indicating chemical specificity. The steady-state emitter density was found to depend on the liquid medium (Supplementary Figs. 7 and 8), illumination power (Supplementary Fig. 9) and macroscopic flow over the crystal (Supplementary Fig. 10).Considering the dynamics of these fluorescent emitters, a strik-ing observation is the presence of fluorescent trajectories on the crystal surface (Supplementary Video 2). These trajectories indicate the correlated activation of neighbouring defects, corresponding to molecular random walks. By linking the super-resolved localiza-tions of emitters (Methods), we can extract the associated trajectories  (Fig. 2b). Previously observed trajectories on plasma-exposed hBN in water and binary mixtures of water and alcohols were attributed to proton hopping20,21. A similar phenomenology is thus expected here, with (1) defect activation due to reversible charge transfer from and to the solvent and (2) correlated activation of neighbouring defects, medi-ated by the lateral motion of charge-bearing solvent molecules that remain physisorbed on the crystal surface. However, the aprotic nature of some solvents used here, as well as the spectral differences, points to a distinct emitter type and reactivity. Since neither pristine hBN nor the liquids used here possess visible-range electronic transitions, the activation of emitters at the interface must arise from a rearrangement of the electronic structure that generates this new optically addressable electronic transition. This could be explained by the chemisorption of organic molecules onto hBN defects30. The chemical selectivity (Fig. 1e)  suggests that a necessary and sufficient condition for a pure liquid to activate native hBN defects is the presence of a carbon atom in its molecular structure. This finding aligns with recent research, indicating the crucial role of carbon in activating visible-range emitters in hBN18. Although a direct observation of the exact chemical structure of the emitters is challenging, the physicochemical interactions between the defects and the liquid along with the photophysical properties of the emitters will guide structural assignment.Analysis of emitter dynamicsTo quantitatively understand the trajectories, we acquired 50,000 6-ms-long frames of a 13 × 13 μm2 crystal area in the steady state, yield-ing 700,000 localizations leading to 100,000 trajectories. Figure 2a shows a subset of these trajectories overlaid with the super-resolved image obtained from 5,000 frames. This visualization reveals that although some trajectories exhibit a free-hopping behaviour, oth-ers remain trapped for extended periods, resulting in bright spots in  the super-resolved image. Most emitters are active for tens to  hundreds of milliseconds, but a fraction displays stability for several seconds (Fig. 2c), with some activations lasting for over a minute (Sup-plementary Fig. 11). As measurements were conducted at a low emit-ter density where trajectories do not split or merge with statistical  significance, we assign trajectories to single molecules binding to single defects.We analysed the observed molecular random walks through their trajectory residence times on the crystal surface TTres, which comprise complex information on both chemisorption energy at defect sites and physisorption energy on pristine hBN in between defects, as well as the residence times of molecules at single-defect sites TDres corre-sponding to chemisorption only. Both residence times were found to follow a double exponential decay (Fig. 2d), with slow exponential decay components τDres = 35 ± 1 ms and τTres = 82 ± 3 ms. Assuming  that single-defect residence times follow the Arrhenius equation τDres = ν−1eΔG/kT  where ν ≈ 1012−1013 s−1 is a molecular attempt rate21,  we obtain a desorption energy barrier ΔG ≈ 24–27kT ≈ 0.6–0.7 eV, which is larger than typical physisorption energies (tens of millielectronvolts) and smaller than covalent bonding energies (several electronvolts)31. This can be rationalized in terms of a lowered energy barrier under illumination20, consistent with the observed light-induced reduction −300 0 300One-dimensional displacement (nm)10–410–3PDF (nm–1)HoppingTrappedSingledefectcb0 300 600Residence time (ms)Probability (ms–1)10–610–410–20 50 100Lag time (ms)0150300Mode width (nm)Trajectoryτ = 6 msτ = 142 msTrappedHoppingdef500 nmTrappedHoppingaPristine hBNin liquid Native defects0 3 s00.50 0.55 0.60Time (s)0Intensity (photons per 6 ms)Residence timeTrapped10310310–2Fig. 2 | Surface of pristine hBN reveals interfacial molecular dynamics.  a, Overlay of a super-resolved image from 5,000 frames showing hopping emitters in isopropanol as the linked trajectories, as well as trapped spots.  b, Artist’s view of the correlated activation of neighbouring defects leading to the trajectories. c, Representative intensity traces from the same images, taken from 7 × 7 pixel bins around emitters (dashed yellow box in Fig. 1c) with 6 ms exposure time. The top right trace corresponds to a long defect activation. The top left trace corresponds to a short activation of the same defect, magnified in the bottom panel. d, Distribution of residence times on single defects and for the entire trajectories. The dotted lines are fits to a two-component exponential decay. e, Displacement probability density functions (PDF) of the trajectories after different lag times τ = 6 ms, 24 ms, 66 ms,142 ms. The dashed lines are fits to two-component Gaussians. f, Visualizing the evolution of the two modes of the Gaussian fit in e with increasing lag time. The central region, corresponding to the trapped state, remains of a constant width, whereas the tails, corresponding to hopping, enlarge with time. The solid line is a fit to a standard diffusion curve.http://www.nature.com/naturematerialsNature Materials | Volume 22 | October 2023 | 1236–1242 1239Article https://doi.org/10.1038/s41563-023-01658-2in the number of emitters (Fig. 1d) as well as the illumination power dependency of the density of emitters (Supplementary Fig. 9).Turning our attention to emitter motion, we find that their one-dimensional displacement probability density function32 PDF(x,τ) follows a two-component Gaussian distribution (Fig. 2e). The central part of the distribution remains of a constant size (~20 nm) correspond-ing to the localization uncertainty when a trajectory is trapped. The tails of the distribution, however, enlarge with increasing lag time, characterizing the hopping events. As shown in Fig. 2f, the hopping tail size scales as √2Dτ  with lag time τ, which corresponds to Brownian diffusion with a diffusion coefficient D = 9.1 × 10−14 m2 s–1, which is over four orders of magnitude slower than bulk liquid molecular diffusion coefficients. Although this observation makes it impossible to resolve the molecular travel time between defects, this slowdown enables  the detection of bright emission from localized spots, which we now propose as a spectral sensing tool.Spectral properties and solvatochromic sensingWe examined the spectral response of the emitters to their liquid environment using spectral single-molecule localization microscopy (sSMLM)28, which enables simultaneous localization and spectral characterization (Fig. 3a and Methods). Single-emitter spectra were found to be homogeneously distributed, with the appearance of a single  population of emitters when exposed to the same liquid environment. Ensemble-averaged spectra were consistently characterized by two peaks, classically attributed to the zero-phonon line (ZPL) and the phonon side band (PSB) for emitters embedded in a matrix (Fig. 3b).Interestingly, these ensemble spectra appeared to strongly depend on the activating liquid, more precisely on its static dielectric constant ϵliq. Figure 3b shows the spectra of emitters obtained in the following liquids of increasing polarity: pentane, tert-butanol and dime-thyl sulfoxide. A notable polarity-induced solvatochromic redshift in the emission was gradually observed from non-polar pentane (615 nm) to the more polar tert-butanol (626 nm) to the highly polar dimethyl sulfoxide (641 nm). Beyond the ZPL shift, we observed changes in the PSB, which is less clearly defined for polar solvents. Figure 3c shows the centre wavelengths of both peaks as obtained from fitting to a two-Lorentzian model for several solvents ordered by increasing ϵliq. We tested 1-pentanol, isopropanol and methanol on top of previously introduced liquids to interpolate dielectric constant values and found that both ZPL and PSB are redshifted by over 25 nm (~80 meV) in highly polar liquids compared with non-polar alkanes. In the range of ϵliq < 25, a linear dependence was observed between the ZPL wavelength and dielectric constant (Fig. 3c, dashed line), with a slope of approximately 1 nm per unit.The Jablonski diagram (Fig. 3d) illustrates the process giving rise to the observed spectra; by absorbing a photon (excitation; green arrow), the emitter is excited to a dipolar state, which interacts with the solvent and relaxes before radiating (dipole relaxation; curved arrow). Direct evidence of the dipolar nature of our liquid-activated emitters is presented in Fig. 3e where either the linearly polarized light used for excitation (green) or PL emission (orange) was rotated as the signal was 620 640 660 680Peak wavelength (nm)010203040Solvent dielectric constantsSMLMZPL PSBExcitationPSBDipolerelaxation+PhononZPLEmissionExcitationTert­butanol1­pentanolDimethyl sulfoxideIsopropanolEthanolMethanolAcetonitrilePentanesSMLM600 625 650 675 700Wavelength (nm)00.51.0Normalized PLTert­butanolDimethyl sulfoxidePentaneCrystal phononsMolecular dipolesabcdefPolarizationmeasurementsα0°60°Fig. 3 | Spectral properties of surface dipole emitters coupled to both solid and liquid environments. a, sSMLM splits the fluorescence signal from an emitter into a localization component (left; single-molecule localization microscopy (sSMLM)) and a spectral component (right) on the same camera chip. b, Ensemble spectra of liquid-activated emitters in different type-I solvents, exhibiting a clear ZPL and PSB. c, Visualizing the wavelength shifts of both ZPL (circles) and PSB (squares), which correlate with the dielectric constant of the liquid. The peak positions were obtained by fitting to a sum of Lorentzians. The dashed line indicates the linear solvatochromic range, with a slope of 1 nm per unit. The error bars correspond to standard deviations of fitting parameters  from groups of 100 single-molecule spectra. d, Jablonski diagram of processes at play: 561 nm laser excitation induces a dipolar excited state, which can directly emit (orange arrow; ZPL) or with the emission of a phonon (red arrow; PSB).  e, Normalized intensity as a function of input-light or output-light polarization angle α relative to the emitter axis. The solid lines correspond to fits to ideal electric dipole emission cos2α + constant. Supplementary Fig. 11 provides more details. f, Sketch of an excited emitter that can interact with the crystal through phonons as well as with the surrounding molecules (yellow ellipses).0 5 10Time (s)00.51.0Intensity (×106 photons s–1)−20 0 20Time delay (ns)00.51.0Normalized coincidencesa bc−25 0 25Time delay (ns)00.51.0Norm. coincidences∆λ ≈ 20 nmg(2)(0) = 0.45g(2)(0) = 0.25HexadecaneAcetonitrileSPD1SPD250%50%Emitter g(2)Fig. 4 | Quantum emission from liquid-activated emitters. a, Representative PL trace from an isolated emitter in acetonitrile under 0.7 mW confocal excitation. The shaded region corresponds to the 10-s-long trace used for photon statistics. b, Normalized coincidences g(2) measured from time-correlated single-photon counting in a Hanbury Brown and Twiss geometry, in hexadecane and acetonitrile. In both liquids, a pronounced antibunching is observed with g(2)(0) < 0.5 without background correction, proving the single-photon emission. The fluorescent lifetimes, corresponding to the width of the antibunching dip, were found to be 2.73 ± 0.09 ns and 2.20 ± 0.20 ns for acetonitrile and hexadecane, respectively. The spectra are shown in the inset, demonstrating liquid-tunable single-photon emission from 615 to 636 nm. c, Single-photon statistics under pulsed laser excitation, showing a suppression of the central peak due to antibunching.http://www.nature.com/naturematerialsNature Materials | Volume 22 | October 2023 | 1236–1242 1240Article https://doi.org/10.1038/s41563-023-01658-2monitored through an analyser (Supplementary Fig. 11). As sketched in Fig. 3f, the solvatochromic redshift can be described by the presence of liquid molecular dipoles stabilizing the excited-state dipole, thus lower-ing its energy and redshifting the ZPL. After this step (Fig. 3d), the tran-sition back to the ground state can occur in two ways: direct emission of a photon (ZPL; orange arrow) and phonon-assisted emission (PSB; red arrow), which is redshifted compared with the ZPL as a fraction of the energy leads to lattice vibrations (phonon; zigzag arrow). Our analysis of the PSBs revealed increasing phonon broadening and decreasing  PSB content with increasing solvent polarity (Supplementary  Discussion and Supplementary Fig. 12). In the case of a mixture of  polar ethanol and apolar heptane, the spectrum was very similar to that  of the polar liquid, demonstrating a strong affinity between the  emitters and polar molecules (Supplementary Figs. 13 and 14).Time-resolved measurements reveal quantum emissionTo prove that the measured fluorescence originates from single-photon emitters, we performed time-correlated photon counting. For this, a 0.7 mW continuous-wave 561 nm laser beam was focused to an ~1 μm2 spot onto hBN crystals in liquid, and the fluorescence signal was col-lected by two single-photon detectors in a Hanbury Brown and Twiss interferometer configuration (Fig. 4a, inset). The typical time trace of a stable single emitter in acetonitrile is shown in Fig. 4a. The analysis of the photon arrival times from a 10 s window shows a clear photon antibunching dip g(2)(0) = 0.25 ± 0.02 at zero delay time, demonstrat-ing single-photon emission (Fig. 4b). This result implies that the bright spots are single emitters and not clusters, and therefore, their optical readout truly reports on the nanoscale properties of the liquid. This activation of quantum emission through the strong chemisorption interaction between a single activating molecule and a single defect was observed in carbon nanotubes33, but the mechanism at play here exhibits the particularities of being transient and observed in liquid. Photon statistics under pulsed excitation show the suppression of the correlated pulse peak at zero delay time, confirming single-photon emission (Fig. 4c). This feature was also found in hexadecane with a measured g(2)(0) = 0.45 ± 0.04. We, thus, demonstrated liquid-tunable single-photon emission with a ZPL shift of 21 nm (Fig. 4b, inset). This shift is comparable with those achieved by hBN defects in response to strain34 or electric fields35, which shows the potential of liquid-activated emitters as dielectric environment sensors.Integration in single-digit nanofluidic systemsBuilding on the characterization of the emission in bulk liquids, we probed hBN–liquid interfaces in molecular confinement in two-dimensional nanoslits. As sketched in Fig. 5a, the nanoslits are obtained by the van der Waals assembly of heterostructures compris-ing three crystals: bottom, spacer and top. The top crystal was cho-sen to be muscovite mica for its transparency and lack of fluorescent properties, and the bottom crystal was pristine hBN to be activated by the liquid. The middle crystal, composed of few-layer graphene patterned by electron-beam lithography, acted as a spacer, defining a slit-shaped channel between the hBN and mica crystals, whose height 1 µm2.4 nmAtomic force microscopyNanoslitsSingle-molecule localization microscopy0 200Photons per localization0Counts103ab dMica tophBN bottomGraphene spacer DefectW = 150 nmcMasked hBNNanoslits1 µmNanoslitsMasked hBNMasked  hBNe hEthanolBare2.4 nm1.4 nm600 650 700Wavelength (nm)Normalized PLMaskedgAcetonitrile600 650 700Wavelength (nm)Normalized PLiBare2.4 nm1.4 nmMaskedh200 nmf1 µmNanoslitsh εconfεliq = 24εliq = 38εconf ≈ 11εconf ≈ 7ℓdipFig. 5 | Nanoslit-embedded liquid-activated emitters. a, Sketch of the heterostructure nanoslit device. The red glow indicates an emitter inside the nanoslit. b, Overlay of a super-resolved image of masked ethanol-activated hBN and the atomic force microscopy mapping of the graphene spacers. c, Optical micrograph of the heterostructure. On the purple-coloured part of the image, only graphene spacers on the hBN bottom crystal are present, which leads to masked hBN. The bottom blue region corresponds to the full heterostructure with slits. d, Super-resolved image of acetonitrile-activated emitters embedded in 2.4-nm-high nanoslits, from 30,000 frames with 20 ms exposure time and 1.4 kW cm–2 illumination. e, Comparison of localization intensity distributions for masked hBN and 2.4 nm nanoslits in acetonitrile, showing no loss of photons but an overall reduction in the number of localizations (Supplementary Fig. 18). f, Illustration of the effect of confinement: the liquid dielectric constant can be changed by confinement (ϵconf) and the defect dipole can interact with solvent molecules (yellow ellipses) within a range ℓdip, comparable with the confinement size h. g, Representative trajectories in the 2.4-nm-high nanoslits filled with ethanol, overlaid with the super-resolved image. h,i, sSMLM spectra of liquid-activated defects in nanoslits filled with ethanol (h) and acetonitrile (i). The confinement size is tuned from the open geometry to 2.4 nm down to 1.4 nm. The solid lines correspond to two-component Lorentzian fits, and the black dashes indicate the extracted ZPL position.http://www.nature.com/naturematerialsNature Materials | Volume 22 | October 2023 | 1236–1242 1241Article https://doi.org/10.1038/s41563-023-01658-2was set by the number of 3.4-Å-thick graphene layers4 (Methods and Supplementary Fig. 15).An optical micrograph of a device with h = 2.4 nm is provided in Fig. 5c. Here the bottom blue region corresponds to the full hetero-structure with nanoslits, and the top purple region corresponds to the open hBN crystal masked by graphene spacers but without encap-sulation by mica. We verified that covering the pristine hBN crystal with a patterned few-layer graphene crystal masks liquid-activated emitters, as was observed for other types of hBN emitters36,37. On bare hBN, emitters are randomly distributed (Supplementary Fig. 16), but the graphene mask allows for their precise positioning on the basal plane of hBN in liquid. An overlay of the graphene-spacer atomic force microscopy topography and the super-resolved image is shown in  Fig. 5b, demonstrating the correspondence between the lithographi-cally defined graphene pattern and the optically measured fluores-cence from masked hBN. We further verified that capping the masked hBN with mica does not quench its fluorescence, allowing the direct imaging of emitters in a confined liquid (Fig. 5d). Figure 5e shows that the localization intensity distributions with and without the confin-ing mica top are similar, but the number of emitters in confinement  are reduced by two-thirds. Single-defect residence times inside nanoslits were found to be slightly longer in nanoslits than in masked hBN (Supplementary Fig. 17). Therefore, the observed threefold decrease in the number of emitters comes not from faster photobleach-ing of the emitters but from the confinement-induced slowdown of their activation kinetics.Integrating emitters into nanofluidic structures allows probing the effect of confinement on the liquid structure and dynamics. We first confirm the observation of trajectories in 1.4 nm confinement (Fig. 5g), where a set of emitters is shown in the 150-nm-wide slit. We then focus on the spectral properties of the confined emitters, which can be robustly extracted through sSMLM with relatively low numbers  of localizations (<1,000). We present sSMLM spectra obtained in the bare, masked and confined geometries for ethanol (Fig. 5h) and ace-tonitrile (Fig. 5i). For both solvents, bringing the confinement size from 2.4 nm down to 1.4 nm leads to a clear blueshift (pinpointed by the dashes) indicating the ZPL position. For ethanol, the ZPL blueshifts from 637 to 624 ± 2 nm and from 636 to 621 ± 1 nm for acetonitrile under 1.4 nm confinement, bringing the spectral signature of a strongly polar solvent close to that of non-polar alkanes. This substantial confinement-induced blueshift suggests that emitters experience a reduced dielectric constant of ϵconf ≈ 11 in ethanol and 7 in acetonitrile when considering Fig. 3c as a calibration curve.Bringing the top mica wall close to the emitters can impact their emission in two ways. First, confinement by the opposite wall of low-ered dielectric constant ϵwall ≈ 8 could reduce the effective number of solvent molecules interacting with the emitter and confine the electric field lines within the slit38, potentially destabilizing the excited state. Second, the emission can be affected by reducing the out-of-plane component of the liquid dielectric tensor8. As depicted in Fig. 5f, the interaction range between a solvent molecule with dipole μS and the defect with dipole μD is given by ℓdip = (µSµD/2πϵokBT)1/3. Assum-ing μD ≈ a × e, where a ≈ 0.25 nm is the in-plane lattice parameter of hBN and e is the elementary charge, we estimate μD ≈ 12D. Using the dipole moments of ethanol (1.7D) and acetonitrile (3.4D), we find ℓdip ≈ 1.0 and 1.3 nm, respectively, which are smaller than the height of the nanoslit. Hence, the observed effect under 1.4 nm confinement might not solely arise from geometrical effects due to the proximity of the top wall, and could be explained by the confinement-induced reduction in the solvent dielectric constant, as observed for water8 and predicted for other liquids39. These results consolidate the picture of dipolar-environment-tuned emitters, whose properties are affected by changes in the sensing hemisphere with volume ~2/3πℓ3dip, which encloses fewer than 100 molecules in the case of acetonitrile. Beyond the passive diffusion and dielectric sensing of confined liquids, tracking emitter dynamics in confinement may be used to directly image the nanoscale flow and study its interplay with defects40.OutlookCrystals of hBN, already known for their exceptional optical properties, exhibit a peculiar interaction with liquids. When in contact with organic solvents, native point defects on the atomically smooth surface of the crystal become emissive. This unique system, where the encounter of a single defect with a single organic molecule yields a single-photon emit-ter, combines solid-state emitters and organic fluorophores, providing a new tool for studying solid–liquid interfaces. We demonstrated two sensing approaches using liquid-activated hBN: the activation dynam-ics provide insights into interfacial charge transfer between defects and single molecules, whereas the emission spectra of the emitters offer infor-mation about the nanoscale dielectric environment. These phenomena  were found to hold in confinement as small as a few nanometres, where only ensemble-averaged measurement techniques have been successful so far. As it relies on common samples and widely available single-molecule microscopy techniques, this approach could be readily applied for  optical imaging and sensing in nanofluidic systems in operando.Online contentAny methods, additional references, Nature Portfolio reporting sum-maries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contri-butions and competing interests; and statements of data and code avail-ability are available at https://doi.org/10.1038/s41563-023-01658-2.References1.  Faucher, S. et al. Critical knowledge gaps in mass transport through single-digit nanopores: a review and perspective.  J. Phys. Chem. C 123, 21309–21326 (2019).2.  Bocquet, L. Nanofluidics coming of age. Nat. Mater. 19, 254–256 (2020).3.  Majumder, M., Chopra, N., Andrews, R. & Hinds, B. J. Enhanced flow in carbon nanotubes. Nature 438, 44–44 (2005).4.  Radha, B. et al. Molecular transport through capillaries made with atomic-scale precision. Nature 538, 222–225 (2016).5.  Secchi, E. et al. Massive radius-dependent flow slippage in carbon nanotubes. Nature 537, 210–213 (2016).6.  Mouterde, T. et al. Molecular streaming and its voltage control in ångström-scale channels. Nature 567, 87–90 (2019).7.  Agrawal, K. V., Shimizu, S., Drahushuk, L. W., Kilcoyne, D. &  Strano, M. S. Observation of extreme phase transition temperatures of water confined inside isolated carbon nanotubes. Nat. Nanotechnol. 12, 267–273 (2017).8.  Fumagalli, L. et al. Anomalously low dielectric constant of confined water. Science 360, 1339–1342 (2018).9.  Grimm, J. B. & Lavis, L. D. Caveat fluorophore: an insiders’ guide to small-molecule fluorescent labels. Nat. Methods 19, 149–158 (2022).10.  Schirhagl, R., Chang, K., Loretz, M. & Degen, C. L. Nitrogen- vacancy centers in diamond: nanoscale sensors for physics and biology. Annu. Rev. Phys. Chem. 65, 83–105 (2014).11.  Aslam, N. et al. Nanoscale nuclear magnetic resonance with chemical resolution. Science 357, 67–71 (2017).12.  Siria, A. et al. Giant osmotic energy conversion measured in a single transmembrane boron nitride nanotube. Nature 494, 455–458 (2013).13.  Tran, T. T., Bray, K., Ford, M. J., Toth, M. & Aharonovich, I.  Quantum emission from hexagonal boron nitride monolayers. Nat. Nanotechnol. 11, 37–41 (2016).14.  Hayee, F. et al. Revealing multiple classes of stable quantum emitters in hexagonal boron nitride with correlated optical and electron microscopy. Nat. Mater. 19, 534–539 (2020).http://www.nature.com/naturematerialshttps://doi.org/10.1038/s41563-023-01658-2Nature Materials | Volume 22 | October 2023 | 1236–1242 1242Article https://doi.org/10.1038/s41563-023-01658-215.  Kianinia, M., Xu, Z.-Q., Toth, M. & Aharonovich, I. Quantum emitters in 2D materials: emitter engineering, photophysics, and integration in photonic nanostructures. Appl. Phys. Rev. 9, 011306 (2022).16.  Choi, S. et al. Engineering and localization of quantum emitters in large hexagonal boron nitride layers. ACS Appl. Mater. Interfaces 8, 29642–29648 (2016).17.  Fournier, C. et al. Position-controlled quantum emitters with reproducible emission wavelength in hexagonal boron nitride. Nat. Commun. 12, 3779 (2021).18.  Mendelson, N. et al. Identifying carbon as the source of visible single-photon emission from hexagonal boron nitride. Nat. Mater. 20, 321–328 (2021).19.  Koperski, M. et al. Midgap radiative centers in carbon-enriched hexagonal boron nitride. Proc. Natl Acad. Sci. 117, 13214–13219 (2020).20.  Comtet, J. et al. Direct observation of water-mediated single- proton transport between hBN surface defects. Nat. Nanotechnol. 15, 598–604 (2020).21.  Comtet, J. et al. Anomalous interfacial dynamics of single proton charges in binary aqueous solutions. Sci. Adv. 7, eabg8568 (2021).22.  Glushkov, E. et al. Engineering optically active defects in hexagonal boron nitride using focused ion beam and water.  ACS Nano 16, 3695–3703 (2022).23.  Na, Y. S. et al. Modulation of optical and electrical properties in hexagonal boron nitride by defects induced via oxygen plasma treatment. 2D Mater. 8, 045041 (2021).24.  Taniguchi, T. & Watanabe, K. Synthesis of high-purity boron  nitride single crystals under high pressure by using Ba–BN solvent. J. Cryst. Growth 303, 525–529 (2007).25.  Wong, D. et al. Characterization and manipulation of individual defects in insulating hexagonal boron nitride using scanning tunnelling microscopy. Nat. Nanotechnol. 10, 949–953 (2015).26.  Schué, L., Stenger, I., Fossard, F., Loiseau, A. & Barjon, J. Characterization methods dedicated to nanometer-thick hBN layers. 2D Mater. 4, 015028 (2016).27.  Henck, H. et al. Direct observation of the band structure in bulk hexagonal boron nitride. Phys. Rev. B 95, 085410 (2017).28.  Comtet, J. et al. Wide-field spectral super-resolution mapping of optically active defects in hexagonal boron nitride. Nano Lett. 19, 2516–2523 (2019).29.  Feng, J. et al. Imaging of optically active defects with nanometer resolution. Nano Lett. 18, 1739–1744 (2018).30.  Lvova, N. & Ananina, O. Y. Theoretical study of the adsorption properties of porous boron nitride nanosheets. Comput. Mater. Sci. 115, 11–17 (2016).31.  Huber, F. et al. Chemical bond formation showing a transition from physisorption to chemisorption. Science 366, 235–238 (2019).32.  Anthony, S., Zhang, L. & Granick, S. Methods to track single- molecule trajectories. Langmuir 22, 5266–5272 (2006).33.  He, X. et al. Tunable room-temperature single-photon emission  at telecom wavelengths from sp3 defects in carbon nanotubes. Nat. Photon. 11, 577–582 (2017).34.  Mendelson, N., Doherty, M., Toth, M., Aharonovich, I. & Tran, T. T.  Strain-induced modification of the optical characteristics of quantum emitters in hexagonal boron nitride. Adv. Mater. 32, 1908316 (2020).35.  Nikolay, N. et al. Very large and reversible Stark-shift tuning of single emitters in layered hexagonal boron nitride. Phys. Rev. Appl. 11, 041001 (2019).36.  Xu, Z.-Q. et al. Charge and energy transfer of quantum emitters in 2D heterostructures. 2D Mater. 7, 031001 (2020).37.  Stewart, J. C. et al. Quantum emitter localization in layer- engineered hexagonal boron nitride. ACS Nano 15, 13591–13603 (2021).38.  Kavokine, N., Marbach, S., Siria, A. & Bocquet, L. Ionic Coulomb blockade as a fractional Wien effect. Nat. Nanotechnol. 14, 573–578 (2019).39.  Motevaselian, M. H. & Aluru, N. R. Universal reduction in dielectric response of confined fluids. ACS Nano 14, 12761–12770 (2020).40.  Seal, A. & Govind Rajan, A. Modulating water slip using atomic-scale defects: friction on realistic hexagonal boron nitride surfaces. Nano Lett. 21, 8008–8016 (2021).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/.© The Author(s) 2023http://www.nature.com/naturematerialshttp://creativecommons.org/licenses/by/4.0/http://creativecommons.org/licenses/by/4.0/Nature MaterialsArticle https://doi.org/10.1038/s41563-023-01658-2MethodsSample preparationPristine hBN flakes from high-quality crystals24 were exfoliated onto borosilicate glass coverslips (no. 1.5 Micro Coverglass, Electron Micros-copy Sciences, 25 mm in diameter, 170 μm thick), using low-adhesion blue tape. The glass coverslips were cleaned either by (1) sonication in acetone followed by rinsing in isopropanol or (2) sonication in 2% Hellmanex III glassware cleaning agent solution, followed by sonication in deionized water. No differences were observed between (1) and (2) above. The coverslips were rinsed three times in the last solution, after which they were dried with a nitrogen gun covered with a 20-nm-pore filter. Adhesion to the substrate on exfoliation was promoted by oxygen-plasma cleaning of the coverslip (2 min, 100 W). High-quality crystals purchased from HQ Graphene and hBN crystals obtained by the crystallization of hBN from molten iron in a nitrogen–hydrogen atmosphere were also tested, and exhibited no notable difference with the samples grown by us. The crystals were immediately immersed fol-lowing exfoliation, and the PEEK chamber was thoroughly rinsed three times at each solvent exchange with a fresh solvent. The chamber was covered by an additional glass coverslip to prevent solvent evaporation and contamination. Supplementary Methods gives the fabrication details for the microfluidic flow cell used for flow measurements.  Supplementary Table 1 lists the solvents used.Nanoslits are made of a van der Waals heterostructure of a spacer layer sandwiched between a top layer and a bottom layer following the same protocol as previously reported4. Here the van der Waals stack is composed of top mica/graphene spacer/bottom hBN. In brief, thin (few atomic layers) graphene was first patterned via electron-beam lithography into parallel strips with a width of 1 μm and a separation of 150 nm. A mica crystal (~200 nm thick) was then transferred on top of the graphene spacer, via a polymethyl-methacrylate-based transfer method. Then, this mica–graphene spacer stack was lifted and trans-ferred onto a freshly exfoliated hBN layer. This whole stack was then transferred onto a glass coverslip for imaging. The slit dimensions of the final device are shown in Fig. 5a, with a width of approximately 150 nm, length of 20 μm and height equivalent to the thickness of the graphene spacer. Supplementary Fig. 15 provide details of the fabrica-tion flowchart and materials.Chemicals usedAll the chemicals were purchased with the maximum purity grade available, and some were purchased in their anhydrous forms. Supplementary Table 1 provides the full list. No effect of solvent purity or residual water traces on hBN fluorescence activation was observed.Optical microscopyWide-field imaging was performed on a custom wide-field fluores-cence microscope, as described elsewhere28. Briefly, the emitters are excited using a 561 nm laser (monolithic laser combiner 400B, Agilent Technologies), which is collimated and focused on the back focal plane  of a high-numerical-aperture oil-immersion microscope objective (Olympus TIRFM 100× with a numerical aperture of 1.45). This con-figuration leads to wide-field illumination of the sample in a circle with ~25 μm diameter. Fluorescence emission from the sample is col-lected by the same objective and spectrally separated from the exci-tation light using dichroic and emission filters (ZT488/561rpc-UF1 and ZET488/561m, Chroma) before being projected on an electron-multiplying charge-coupled device camera (Andor iXon Ultra 897) with an electron multiplication gain of 150. An additional spectral path, mounted in parallel to the localization path, allows for the simulta-neous measurements of the emission spectra from individual emitters (as described below). The sample itself is mounted in a sealed fluidic chamber, which is placed on a piezoelectric scanner (Nano-Drive, MadCityLabs) for a fine focus. The typical exposure time is 20–50 ms and the typical laser power is 10–100 mW for the wide-field excitation area of 2 × 103 μm2, resulting in a power density of 0.35–3.50 kW cm–2. Unless mentioned otherwise, the illumination power density was set to 3.5 kW cm–2. A typical acquired image stack contained 2,000–10,000 frames.Confocal microscopy measurements were performed in an inverted microscope configuration allowing to image samples in a liquid in custom-made chambers made of a glass ring glued onto the glass coverslip using an epoxy resin (Araldite). The excitation laser was a 561 nm picosecond diode laser (Sepia II, PicoQuant) used either in the pulsed or continuous-wave mode. The emission, after collection with a water-immersion microscope objective (Nikon SR Plan Apo 60× with numerical aperture of 1.27) was split between two fibre-coupled avalanche photodiodes (SPCM-AQRH, Excelitas) in a Hanbury Brown and Twiss configuration. Photon correlation measurements were performed using the PicoHarp time-correlated single-photon count-ing module (PicoQuant).sSMLM procedureThe image stacks acquired from the wide-field microscope were pro-cessed using ThunderSTORM41. The localization uncertainty of emitters is given by σloc ≈ σPSF/√Nloc , where σPSF represents the radial extent of the microscope point spread function (approximately 110 nm), and Nloc is the number of photons counted in the diffraction-limited spot. Supplementary Table 2 lists all the symbols used here. In our imaging conditions, localization uncertainties ranged from 10 to 20 nm. Locali-zations used for emitter counting, super-resolution image rendering and trajectory analysis were filtered to keep only the localization events with 30 < σPSF < 200 nm and filter out artifacts. To generate super-resolved images, individual single-molecule localization micros-copy localization events were rendered as 2D Gaussians with a standard deviation of 20 nm unless specified otherwise, reflecting the typical uncertainty on the emitter’s position.sSMLM was performed following a procedure described in our previous work28, which we summarize here. Emitters in the spatial channel are localized using the ImageJ plugin ThunderSTORM. Briefly, a wavelet filter is applied to each frame, and peaks are then fitted by 2D-integrated Gaussians. The obtained localizations (xloc, yloc) are matched to a reference spectral point (wavelength of 650 nm) through a matrix transformation of the form (xspec, yspec) = A × (xloc, yloc) + B. Here A is a 2 × 2 matrix and B is a vector. The spectrum is then extracted as the vertical profile around (xspec, yspec) in a box of 6 × 40 pixels centred at the emitter. The pixel values are translated as spectral intensities through the linear relationship Δyspec ≈ a × λ, with a ≈ 0.25 pixels nm–1. This value and transformation matrix coefficients were obtained using fiducial markers emitting at known wavelengths. As described in our previous work, this was achieved by using broadband emitters (fluo-rescent beads) and narrow bandpass filters in the emission path. The signal-to-noise ratio of the single-emitter spectrum is evaluated and all the spectra passing a threshold are averaged to yield an ensemble spectrum. We typically used 3,000 frames for each solvent and set the threshold to Iloc > 300 photons in all the measurements.Ensemble sSMLM spectra were maximum normalized and fit  using the Python 3.7 package LMFIT42, to a model composed of two Lorentzians corresponding to the ZPL and PSB and a linear background. As shown in Fig. 5g,h, only the part of the normalized signal emerging more than 20% over the background was fit, as the background-induced spectra tails could pose fitting issues. For spectra in nanoslits, a spatial filter was applied to avoid counting the signal from contamination between nanoslits (Supplementary Fig. 18). The error bars on the ZPL and PSB wavelengths (Fig. 3c) as well as ZPL–PSB detuning (Supple-mentary Fig. 12) were obtained by evaluating the standard deviation of the parameters obtained from fitting the average of randomly chosen groups of 100 single-emitter sSMLM spectra. This ensured a sufficient signal-to-noise ratio for fitting and reflecting variations within the ensemble spectra.http://www.nature.com/naturematerialsNature MaterialsArticle https://doi.org/10.1038/s41563-023-01658-2Trajectory analysisWide-field fluorescence frames acquired with the electron-multiplying charge-coupled device camera were first localized using Thunder-STORM41. Trajectories were obtained by applying the Crocker–Grier linking algorithm43 to the localization microscopy tables. We used the implementation of the algorithm provided by the Python pack-age trackpy44. Briefly, for each localization event at a given frame, the algorithm links another event if it is found at the next frame within a specified search range. This search range was set to 120 nm in Fig. 2, as the probability of having a one-dimensional displacement exceeding this value is less than 1% according to the analysis shown in Fig. 2e.In the statistical analysis of trajectories, the emitter position is described as a random variable of time X(t). The trajectory residence time on the crystal surface TTres corresponds to the duration of this single trajectory. The one-dimensional displacement probability dis-tribution function is given by the PDF (x,τ) = P(X(t + τ) – X(t) = x). The residence times of the molecules at single-defect sites TDres  were obtained using the same linking algorithm with a linking range of 35 nm (about twice the typical localization uncertainty) instead of 120 nm.Data availabilityThe linked localization table corresponding to Supplementary Video 2 and Fig. 2 and the corresponding raw frames are available via Zenodo at https://doi.org/10.5281/zenodo.8087398. Source data are provided with this paper.References41.  Ovesný, M., Křížek, P., Borkovec, J., Švindrych, Z. & Hagen, G. M.  ThunderSTORM: a comprehensive ImageJ plug-in for PALM and  STORM data analysis and super-resolution imaging. Bioinformatics 30, 2389–2390 (2014).42.  Newville, M., Stensitzki, T., Allen, D. B. & Ingargiola, A. LMFIT: non-linear least-square minimization and curve-fitting for Python. Zenodo https://doi.org/10.5281/zenodo.11813 (2014).43.  Crocker, J. C. & Grier, D. G. Methods of digital video microscopy for colloidal studies. J. Colloid Interface Sci. 179, 298–310 (1996).44.  Allan, D. B., Caswell, T., Keim, N. C., van der Wel, C. M. &  Verweij, R. W. soft-matter/trackpy: v0.6.1. Zenodo https://doi.org/ 10.5281/zenodo.7670439 (2023).AcknowledgementsA.R., N.R., E.G., M.L. and T.-H.C. acknowledge funding from the European Research Council (grant 101020445—2D-LIQUID), the Swiss National Science Foundation (grant 200021 192037) and the CCMX Materials Challenge grant ‘Large area growth of 2D materials for device integration’. B. Radha acknowledges funding from the Royal Society University Research Fellowship URF/R1/180127. B. Radha, Y.Y. and G.-H.N. acknowledge funding from the European Research Council (grant 852674—AngstroCAP) and RS enhancement award RF/ERE/210016. A.K. acknowledges funding from the Royal Society International Exchanges Award IES/R1/201028 and EPSRC New Horizons grant EP/V048112/1. K.W. and T.T. acknowledge support from JSPS KAKENHI (grant nos. 20H00354, 21H05233 and 23H02052) and World Premier International Research Center Initiative (WPI), MEXT, Japan. J.C. acknowledges funding from the ANR (grant ‘GUACAmole’ ANR-22-CE06-0003-01). S.R. thanks the Julia Jacobi fund. We thank  E. Lopriore and A. Kis for the atomic force microscopy characterization and A. Söll and Z. Sofer for providing the hBN crystals grown with a different approach. N.R. acknowledges useful conversations with M.-L. Bocquet and V. Artemov. We thank M. Weber and M. Zhang for reporting independent checks of the activation phenomenon.Author contributionsN.R. observed the effect. A.R., B. Radha, J.C., A.K. and N.R. conceived the defect integration in nanoslits. A.R., N.R., M.L., T.-H.C., B. Radha, A.K. and J.C. designed the experiments. E.G. fabricated the marked coverslips. Y.Y. and G.-H.N. fabricated the nanoslit devices under the supervision of B. Radha and A.K. F.B. fabricated the microfluidic flow cell. N.R. performed the sSMLM measurements. N.R. performed the confocal microscopy measurements with the assistance of E.G. N.R. analysed the data and interpreted it with the help of J.C., E.G., M.L.,  B. Rehl, S.R. and A.R. K.W. and T.T. contributed the materials.  N.R. wrote the paper with the assistance of J.C. and inputs from all authors. A.R. supervised the project. All authors discussed the results and commented on the manuscript.Competing interestsThe authors declare no competing interests.Additional informationSupplementary information The online version  contains supplementary material available at  https://doi.org/10.1038/s41563-023-01658-2.Correspondence and requests for materials should be addressed to Nathan Ronceray, Boya Radha or Aleksandra Radenovic.Peer review information Nature Materials thanks the anonymous reviewers for their contribution to the peer review of this work.Reprints and permissions information is available at  www.nature.com/reprints.http://www.nature.com/naturematerialshttps://doi.org/10.5281/zenodo.8087398https://doi.org/10.5281/zenodo.11813https://doi.org/10.5281/zenodo.7670439https://doi.org/10.5281/zenodo.7670439https://doi.org/10.1038/s41563-023-01658-2http://www.nature.com/reprints Liquid-activated quantum emission from pristine hexagonal boron nitride for nanofluidic sensing Liquid-activated fluorescence from pristine hBN Analysis of emitter dynamics Spectral properties and solvatochromic sensing Time-resolved measurements reveal quantum emission Integration in single-digit nanofluidic systems Outlook Online content Fig. 1 Liquid-induced fluorescence from pristine hBN crystals. Fig. 2 Surface of pristine hBN reveals interfacial molecular dynamics. Fig. 3 Spectral properties of surface dipole emitters coupled to both solid and liquid environments. Fig. 4 Quantum emission from liquid-activated emitters. Fig. 5 Nanoslit-embedded liquid-activated emitters.