Nanodiamond arrays on glass for quantification and fluorescence characterisation

TL;DR

This paper presents a precise method for creating nanodiamond arrays on glass to statistically analyze their fluorescence properties, advancing the development of quantum and sensing devices.

Contribution

The study introduces a novel electron beam lithography-based technique for depositing nanodiamonds in predefined arrays for fluorescence characterization.

Findings

Achieved 90% deposition yield across thousands of array sites

Observed broad distribution of fluorescent lifetimes consistent with literature

Demonstrated potential for hybrid photonic device fabrication

Abstract

Quantifying the variation in emission properties of fluorescent nanodiamonds is important for developing their wide-ranging applicability. Directed self-assembly techniques show promise for positioning nanodiamonds precisely enabling such quantification. Here we show an approach for depositing nanodiamonds in pre-determined arrays which are used to gather statistical information about fluorescent lifetimes. The arrays were created via a layer of photoresist patterned with grids of apertures using electron beam lithography and then drop-cast with nanodiamonds. Electron microscopy revealed a 90% average deposition yield across 3,376 populated array sites, with an average of 20 nanodiamonds per site. Confocal microscopy, optimised for nitrogen vacancy fluorescence collection, revealed a broad distribution of fluorescent lifetimes in agreement with literature. This method for statistically quantifying fluorescent nanoparticles provides a step towards fabrication of hybrid photonic devices for applications from quantum cryptography to sensing.

Figures

Tables

Table 1. Table 1: Quantification of nanodiamond material coverage of six of the eight arrays indicated in Figure 1 . Each array was fabricated with 623 apertures developed in PMMA, and after deposition of nanodiamond material and removal of the photoresist, the yield was quantified by eye from SEM micrographs.

Array ID	Empty Sites	Populated Sites	Percent Yield
b	10	613	98
c	30	593	95
d	65	558	89
e	109	514	82
f	42	581	93
g	106	517	82
total	362	3376	90

Full text

Introduction

Diamond has long been studied for its remarkable properties including chemical inertness, biocompatibility, transparency from the ultraviolet to the infra-red range, high thermal conductivity and mechanical strength [1, 2, 3, 4]. Recently, a large amount of research has focused on investigating the optical properties of fluorescent defects, often referred to as colour centres [1]. Detection of individual colour centres by confocal microscopy was demonstrated 20 years ago [5], and over 500 active optical centres have so far been characterised [6], many of which show single photon emission characteristics [7].

One colour centre in particular has received the most attention: the nitrogen-vacancy (NV) centre, consisting of a substitutional nitrogen with a vacancy in its nearest neighbour lattice site [8]. The valence electrons of the nitrogen atom give rise to two possible charge configurations of this centre: neutral NV0, or negatively charged NV-. The emission spectrum of the NV- centre is centred around 700 nm and its standard fluorescent lifetime in bulk diamond is $\sim$12 ns; both of these aspects make it attractive for biological imaging because autofluorescence from surrounding media can be avoided with spectral filters and/or time-gating [9, 10]. Other aspects of the NV centre, including its long-lived and controllable spin state, make it attractive for quantum information purposes [11]. The NV centre has also been presented as a sensor, with changes in the environment (e.g. temperature [12], magnetic field [13], or nearby nuclear spins [14]) giving rise to variations in fluorescence, which is the basis of optically detected magnetic resonance (ODMR) measurements.

The excited state lifetime of the NV centre in diamond varies depending on several factors. The number of impurities in the diamond lattice, presence of non-diamond carbon material, size of the host crystal lattice (e.g. bulk single crystal diamond or nanodiamond), surface modification, and irradiation have all been reported to contribute to lifetime values ranging from 11.2 ns in the bulk to 25 ns [15, 16, 17, 18, 19, 20, 21, 22].

Integrating nanodiamonds with mature photonic devices is essential to using their properties in a range of applications from quantum cryptography to sensing. Here we present a method for positioning as-received ball-milled nanodiamonds in pre-determined locations on glass, and use it as a tool for measuring the fluorescent lifetime statistics of the deposited material.

A number of approaches to locate NV centres in diamond have been developed over the last several years which fall into two broad categories: top-down direct-write techniques, and bottom-up assembly techniques.

Top-down approaches involve the creation of colour centres directly in bulk diamond using high-energy ion, neutron, or electron irradiation followed by annealing [23, 24, 25]. Implanting colour centres with nanoscale precision is a challenge, which can be overcome by using a lithographic mask [26, 27] or a direct-write focussed beam [28, 21]. Vertical distribution in the substrate is another challenge for direct writing techniques [29, 30, 31], especially since the depth of the colour centre below the surface has an impact on the fluorescence and to sense things near the surface the the colour centres have to be proximal to the surface.

Some bottom-up techniques involve nanodiamond particles, using them either as-received from the manufacturer or optically activating them by irradiating the nanodiamond powder before positioning [32]. The nanodiamonds can be deposited by spin-coating a solution on a substrate (leading to single particles dispersed randomly), pre-characterised, then positioned by microscopic probe (colloquially known as ‘pick and place’) [33, 34]; however, this method is time-consuming and not suited to large-scale production.

An alternative to deterministic nanoparticle positioning is to use techniques of directed self-assembly. They can be more efficient, have high spatial control of potential locations, and are often compatible with parallel fabrication methods [35, 36, 37, 38, 39]. To optimise the particle surface properties for self-assembly, nanodiamonds often need significant chemical and physical processing. When suspended in solution, nanodiamonds acquire a zeta potential ($\sim$-30 mV as-received [39]), which leads to readily formed agglomerates: to counteract this, the material can be treated with liquid phase purification or high temperature annealing [3, 40]. Additional processes can be performed to enhance the fabrication, for example reactive ion etching has been used to attach a linker molecule before self-assembly [41]. The method presented here keeps nanodiamond processing to a minimum, does not require any additional steps to process the substrate, and is a low-temperature fabrication method which is compatible with complementary metal-oxide-semiconductor (CMOS) processes.

Directed self-assembly suffers from the fact that pre-selected nanodiamonds cannot be placed at a given location, and the attachment specificity isn’t always perfect. The practical harnessing of directed self-assembly therefore requires using defect-tolerant design considerations. If a device requires a certain number of elements to work, it may be more efficient to overcompensate during fabrication and have many defective elements than it is to laboriously fabricate every element perfectly and achieve 100 percent yield. In this context, self-assembly of nanodiamonds could lead to the creation of a large-scale many-qubit device [42]. Arrays such as the ones presented here are the first steps towards robust hybrid device fabrication on optically transparent substrates for quantum information processing, and towards novel large-area ODMR sensing and biological imaging.

Methods

Array Fabrication

A 3 inch [100] Si wafer was cleaned in an ultrasonic bath of acetone, then ethanol, then isopropyl alcohol (IPA), for 1 minute each. A 100 nm layer of polymethyl methacrylate (PMMA) A2 950 K was spin coated at 4000 rpm for 30 seconds, followed by a soft bake on a hot plate at 200*∘*C for 2 minutes. The wafer was then diced into $\sim$20 mm square samples.

We used Si substrates to optimise the exposure process, and ultimately translated the technique to a transparent substrate: a glass slide (ProSciTech G300 white glass) was diced into $\sim$20 mm squares, and cleaned and spin coated in a similar method to the Si.

For the electron beam lithography (EBL) exposure process, an FEI Nova NanoSEM equipped with the Nabity Pattern Generation System software was used. The current of the 30kV electron beam was measured with a Faraday Cup; in high-vacuum mode ($1\times 10^{-6}$ Torr for the silicon substrates) the current was measured to be 300 pA. The dwell time of the electron beam at each point was calculated by the NPGS software (using the measured value of the current). To mitigate the issue of electrical charge gathering on the surface of the glass slide samples, we used a low-vacuum mode in the SEM which introduces water vapour into the vacuum chamber (known as Variable Pressure Electron Beam Lithography, VP-EBL [43]), at a pressure of 5 Torr and a measured beam current of 101 pA.

The pattern was designed to take into account the diffraction-limited resolution of confocal microscopy which was to be utilised post fabrication: to be confocally resolvable, the array sites needed to be at least 300 nm apart. We designed apertures ranging from 50 nm to 1 $\mu$m in diameter, in arrays with pitches ranging from 1.5 $\mu$m to 10 $\mu$m, and repeated the arrays over areas up to 15 mm2.

The exposed patterns were developed in a 1:3 MIBK:IPA bath for 1 minute, then rinsed with IPA followed by water, and dried with pressurised air. To characterise the developed apertures without exposing the remaining photoresist, an atomic force microscope was used.

As-received nanodiamonds with a nominal diameter of 45 nm (Nabond Technologies, China) were suspended in Milli-Q water at a concentration of 1 $\mu$g/mL and sonicated for 10 minutes to break up the larger agglomerations. The zeta potential and average size of the remaining aggregates was measured with a Malvern Zetasizer, found to be -40 mV and and 124 nm respectively. The solution was drop cast on the developed photoresist and allowed to air dry overnight.

An acetone lift-off step was performed to strip the PMMA and top layer of nanodiamond material away, leaving behind nanodiamonds that had come into contact with the apertures. To remove remaining organic material (solvent or photoresist), the samples were placed in a plasma cleaner (Gatan Model 950, 65 W 13.56 MHz power supply) and subject to 2 minutes of H2/O2 plasma. Characterisation of the topography of the arrays was performed with an FEI Verios SEM, under low accelerating voltage and high stage bias without conductive coating.

Confocal Microscopy

A custom confocal microscope [44] was used to characterise the NV- fluorescence of the nanodiamond arrays at room temperature. A 532 nm continuous wave laser (LaserQuantum GEM 532) was used to excite the sample with 0.5 mW of excitation power. The beam was passed through a 532 nm line filter to filter out fluorescence from a single-mode fibre. A 561 nm dichroic mirror reflected the beam onto the back of a 100x objective with NA = 0.9. The same objective collected the fluorescence, which was passed through a 532 nm notch filter to filter out the reflected laser light. A 697 $\pm$ 75 nm band pass filter was used to block background light and transmit a region that would show the zero phonon lines and phonon sideband of the NV centre emission spectrum. The collected signal was split 90:10, with 10% going to a spectrometer (Princeton Instruments SpectraPro 2500i with PIXIS 100BR camera) to confirm that the emission spectrum matched that of the NV centre, and the remaining 90% being split 50:50 between two avalanche photodiode detectors.

For fluorescent lifetime measurements, a Fianium supercontinuum source was used, set to deliver 0.1 mW of 520 $\pm$ 10 nm light at a pulse rate of 20 MHz (corresponding to a repetition rate of 50 ns). The excitation and emission light was passed through the same filters as for the continuous wave excitation. The time-resolved direct fluorescence decay traces were obtained by a correlator card (Picoquant, TimeHarp 260).

We also excited the same sites with a 633 nm HeNe CW laser, and collected emission through a 633 nm notch filter and a 633 nm long pass filter, to look for silicon vacancy (SiV) centre emission in the as-received nanodiamond material.

Results and Discussion

The fabrication of nanodiamond arrays using EBL has been published elsewhere, with a subsequent diamond film growth step [39]. In this study, we used the same array fabrication method, but we present quantification of the deposition of nanodiamonds on a glass substrate, and we present fluorescent characterisation of the nanodiamond arrays without any other material processing.

Deformation of EBL patterns can be caused by electrical charge gathering on insulating substrates and deflecting the electron beam, however we found that the overall array pattern was not deformed, even on the glass substrate. We attribute the lack of deformation of our arrays to the use of VP-EBL and the simplicity of the overall pattern, each of which minimises the impact of charging. Charging can also be reduced by minimising the area of each aperture, but this requires a precisely focused electron beam. 200 nm circular apertures were written with the VP-EBL using a raster pattern for the exposure. The focus and alignment of the electron beam are critical requirements for accurate EBL, but attaining high precision of these in VP-EBL is challenging. The water vapour molecules that were introduced into the vacuum chamber to minimise charging can cause scattering in the electron beam and decrease the signal-to-noise ratio of the SEM image. Since the clarity of the image is essential to fine tuning the focus of the beam, a noisy image means a precisely focused beam is difficult to achieve.

The deposition process left nanodiamond residue on the PMMA layer in typical coffee-ring patterns as the water evaporated and the droplet shrank. After stripping the photoresist, arrays of nanodiamond material were visible by optical microscope, as were some areas of non-specific deposition, as shown in Figure 1.

Two theories are proposed based on previous literature [39] as to why the nanodiamonds remain at the aperture locations during the photoresist removal step. The negative zeta potential of the particles in solution may be responsible for a van der Waal attraction to the positively charged native silicon oxide layer on the wafer. Alternatively, the surface chemistry of the nanodiamonds (e.g. carboxylic groups) may result in covalent bonding with straggling polymer chains at the edges of the developed apertures in the photoresist. We often saw nanodiamonds in a circular pattern at each site, which may be evidence of the latter hypothesis; however this does not explain non-specific deposition, where nanodiamonds are attached to the substrate in places other than those pre-determined by EBL.

To quantify the yield of deposited nanodiamond material, six individual arrays were imaged by SEM (Figure 1) and the number of populated sites, displayed in Table 1 for each selected array, was deduced by visually counting the number of unpopulated sites. One of the arrays quantified showed 98% yield (613 out of 623 sites populated with nanodiamond material), and the average total yield was 90% over the 6 arrays considered. However, Figure 1 shows some variability and the extent of non-specific deposition.

Nanodiamonds at individual sites were quantified from SEM micrographs for a selection of sites in arrays (h) and (i) shown in Figure 1. Close up SEM micrographs of several sites from array (h) are shown in Figure 2, displaying distributions of number and size of nanodiamonds per site.

Confocal microscopy was performed after developing and before drop-casting the nanodiamond material, and again after the photoresist was removed. The developed PMMA showed no outstanding fluorescence above background levels ($\sim$ 5000 counts per second).

Figure 2 shows an SEM micrograph and a fluorescent map of a region of array (h) shown in Figure 1. Fluorescence spectra was collected from over 200 sites from arrays (h) and (i), for 30 seconds each. 153 of them displayed the characteristic NV centre broad phonon sideband (centred around 700 nm), although the zero phonon lines (NV 0 at 575 nm and NV - at 637 nm) were almost always too weak to distinguish. We also noted no obvious shifting of the phonon sideband.

The quantification of nanodiamonds that do host NV centres in contrast to those that do not could be achieved with an array that has individual particles at each site. As the arrays hold multiple particles per site, we restrict our statistical analysis not to the number of nanodiamonds that host NV centres, but the number of sites that host at least one NV centre, and the number of nanodiamonds over 30 nm present (estimated visually from SEM micrographs). We note that the sites with NV fluorescence and those without show similar numbers of particles per site. The average for the former is 21 with a standard deviation of 10, and the average for the latter is 17 with a standard deviation of 9.

To collect time-dependent fluorescence data, we used pulse laser excitation on those sites that displayed typical NV fluorescence (153 sites in total from 220 sites from arrays (h) and (i) shown in Figure 1). Measuring the fluorescent intensity $I$ as a function of time $t$ allowed us to calculate the fluorescent lifetime $\tau$ of a site by fitting the data (with the Matlab routine lsqcurvefit) to the basic exponential decay equation $I=\alpha e^{-t/\tau}+\beta$ (with $\alpha$ and $\beta$ as scaling parameters to account for total fluorescence and background). Three representative lifetime traces are shown in Figure 3 (a), along with a measurement taken at an off-array site to gauge the instrument response and background fluorescence of the substrate.

Single exponential fitting gives a limited interpretation of the data, so more rigorous analysis was performed with a fitting function of a simple sum of weighted exponentials of the form $a+\sum_{i=1}^{n}b_{i}\exp\left(-t/c_{i}\right)$ for $n$ exponentials with weighting $b_{i}$ and lifetime $c_{i}$, and background $a$. Each time series was analysed for $n=1,2,3,4,5$, and the best fit determined by choosing the number of exponential terms that maximised the adjusted R2 parameter for the series. As constraints, the first exponential was constrained $0<c_{1}100$ ns, and the other exponentials constrained $5<c_{i}<100$ ns $i>1$. We discarded lifetimes close to 5 ns and 100 ns as being artefacts of the fitting routine.

Figure 3 (b) shows a histogram of the lifetimes obtained from the fluorescence time series fits. The data set contained 153 fluorescent sites from a total of 220 analysed sites. The modal fluorescence lifetime was 24.5 ns, with a strongly asymmetric distribution towards longer lifetimes. As expected, all of the measured lifetimes were greater than the spontaneous emission lifetime of NV centres in bulk diamond. Interestingly, most sites showed a single lifetime, and only two sites showed two distinct lifetimes in the range expected for NV centres (between 12 and 50 nm). We do not interpret this as implying that most of the sites hosted single, isolated NV centres. Instead it is likely that NV centres within a single site see similar refractive index environments, and hence radiative density of states, and therefore exhibit similar lifetimes. Further analysis is required to fully understand the significance of this result. Fluorescent lifetime imaging experiments have shown similar lifetime distributions [18, 45], but the work presented here is the first time that this extent of statistical analysis has been conducted by direct measurement of NV centre lifetimes in as-received nanodiamonds.

Finally, we also performed a study of whether any SiV centres were present in the as received nanodiamond material. SiV centres are usually studied only after nanodiamond material has been irradiated [24], or after diamond films are grown via CVD under specific conditions to incorporate silicon atoms [46]. We observed no distinctive peaks at 737 nm corresponding to the zero phonon line of the SiV centre, which is unsurprising. Since the nanodiamonds used in this study were unprocessed, it is unlikely that they hosted many, if any, SiV centres.

Conclusion

We have demonstrated a directed self-assembly method to position nanodiamonds in pre-determined locations at room temperature with minimal material processing. Despite an average of 10% random assembly, the process is faster than manipulating individual particles, and the agility of EBL array pre-determination provides excellent spatial control.

To advance the technique to a point where device fabrication is viable, particle attachment must be studied further. Recently published work [41] suggests an optimal nanodiamond solution concentration to minimise agglomeration. Other researchers [47] have demonstrated ways to prevent agglomeration of nanodiamonds in solution and control their zeta potential, but they add degrees of complexity such as chemical processing and additive surface modification.

The ordered nature of arrays allows for identification of the same particles over different stages of an experiment. Potential future directions include characterising a set of nanodiamonds, deterministically changing their environment or modifying their surface, and re-characterising, with the aim of quantifying the modification.

Arrays of NV centres in bulk diamond have previously been used for magnetic imaging [48], and single nanodiamonds have been used for single-molecule nuclear magnetic resonance [14]. Our nanodiamond arrays provide initial steps for combining these two techniques for chemical sensing and imaging of biological systems, for example the surface of cells.

The work shown here also opens avenues for fabricating defect-tolerant hybrid photonic devices, consisting of, for example, nanodiamonds on dielectric substrates or waveguides.

Acknowledgements

The authors acknowledge the facilities, and the scientific and technical assistance of the RMIT University’s Microscopy & Microanalysis Facility, a linked laboratory of the Australian Microscopy & Microanalysis Research Facility. We also thank Dr. Desmond Lau for assistance with the confocal microscope studies. This work has been supported by ARC grants (FT110100225, FT160100357, LE140100131, CE140100003).

Author contributions statement

A.H. conducted experiments and wrote the manuscript, B.G. and A.G. conceived the experiments and reviewed the manuscript. All authors analysed and interpreted the data.

Additional information

Competing financial interests: The authors declare no competing financial interests.