Frenkel-Poole Mechanism Unveils Black Diamond as Quasi-Epsilon-Near-Zero Surface

Abstract

A recent innovation in diamond technology has been the development of the “black diamond” (BD), a material with very high optical absorption generated by processing the diamond surface with a femtosecond laser. In this work, we investigate the optical behavior of the BD samples to prove a near to zero dielectric permittivity in the high electric field condition, where the Frenkel-Poole (FP) effect takes place. Zero-epsilon materials (ENZ), which represent a singularity in optical materials, are expected to lead to remarkable developments in the fields of integrated photonic devices and optical interconnections. Such a result opens the route to the development of BD-based, novel, functional photonic devices.

1. Introduction

Since the first report in 2014 [1] about the enhanced optical absorption property in the visible range of the fs-laser treated diamond surface, considerable studies have been performed to analyze and optimize the treatments inducing surface nanostructuring (see Figure 1) [2,3]. This material, called “Black Diamond” (BD), is able to combine a very high optical absorptance, due to the formation of regular nanostructures over its surface, with a high mobility of the charge carriers [4], thus conceiving BD as a possible choice to photo-thermionic cathodes in high-temperature solar cells [5]. Furthermore, the possibility of deploying nanostructured diamonds for new applications as near-IR antireflection coatings [6,7] was demonstrated. On the other hand, this treatment introduces irregularities in the atomic arrangement of carbon atoms, and therefore the strain in its lattice, as it happens also in partially graphitized diamond substrates [8,9]. In optimized BD samples, where surface nanostructuring has a pseudo-bidimensional periodicity, an increased number of defects are present in the surface, as demonstrated by the increase of optical absorbance [10], and, correspondingly, the surface becomes orders of magnitude higher conductive than the bulk. Through thermal annealing up to 800 K and stress relaxation, it is possible to observe sheet resistance down to about 25 kΩ at high temperature and spintronic effects through the sequential current measurements in the cryogenic temperature range [11,12]. By enhancing the laser energy further, the spatial arrangement of the created defects percolation path breaks up and the black diamond surface becomes more insulating. On such highly defective BD samples it is possible to perform high electric field current measurements and observe FP conduction mechanism.

The application of large electric fields induces a considerable bending of the conduction band before reaching the breakdown condition. According to the Frenkel-Poole effect (FP), under a high electric field E, a trapped electron can jump to the conduction band without large thermal energy contributions. This induces a consequent measurable nonlinear enhancement of the current density J, following the equations:

(1)$J=J_0\left(T\right)·e^{\beta \sqrt{E}} ,\beta =\raisebox{1ex}{$q^{\raisebox{1ex}{$3$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$}\!\left/ \!\raisebox{-1ex}{$k_{B}T\sqrt{\pi {ϵ}_0{ϵ}_r}$}\right., E=\raisebox{1ex}{$V$}\!\left/ \!\raisebox{-1ex}{$d$}\right.$

where, apart from well-known physical and mathematical constants (k_B the Boltzmann constant, q the electron charge and pi), and easily measurable experimental parameters (T the absolute temperature and d the electrodes distance), the only unknown parameter influencing the current/voltage relation is ε_r, the material dielectric permittivity [13].

In this paper, we demonstrate through the FP effect that by enhancing the surface nanostructuring of a diamond (bulk relative permittivity equal to 5.5), it is possible to obtain its value diminution up to about 10⁻². This is equivalent to considering the black diamond surface almost as an Epsilon Near Zero material (ENZ), representing a singularity in optical materials [14,15,16,17,18]. As observed for other ENZ surfaces [19,20,21], the coupling of such surfaces with nanoantennas allows us to obtain optical non-linearities and strongly improve the directionality of the emitted electromagnetic field; it is expected to lead to remarkable developments in the fields of integrated photonic devices and optical interconnections through the modulation of the Purcell effect [22] (i.e., the increase in spontaneous emission rate generally used to classify nanoantennas). Thanks to the presence of NV-centers, simple diamond nanoparticles also act as active dielectric nanoantennas with a high Purcell-factor [23]. Importantly, we note that NV-centers are also allocable by controlled surface laser writing [24], and it is, therefore, possible to combine in diamonds new electromagnetic effects to quantum properties of single photon emitters, making diamonds a leading contender as election material for the future quantum computer industry.

Furthermore, the demonstrated strong reduction of the dielectric polarizability through surface nanostructuring and the defects ordering on a wide band-gap material, like a diamond, in which n-type doping is physically very difficult, opens the way to a new procedure to realize the ENZ surface for photonic devices. In fact, ENZ materials have different fabrication techniques [25], but the most widely adopted and technologically convenient technique is to strongly enhance the free carrier concentration by heavy n-type doping of silicon or III-V semiconductors and to deploy the charges’ collective motion at plasma frequency [26,27,28]. We infer that highly defected surfaces with shallow electron traps may behave towards electromagnetic waves in a similar manner to the previously mentioned degenerated semiconductors. Even if this physical phenomenon should be confirmed by specific electromagnetic experiments, this is the first scientific report on quasi-ENZ effect on a diamond substrate up to our knowledge, and it can be useful to the large community of scientists working on diamond devices, from biosensors [29,30,31,32] to dosimeters [33,34] and to quantum computing [35,36]. In fact, as easy examples, the localized implementation of ENZ surfaces on diamonds gives the opportunity to improve optical interconnections between individual quantum building blocks (see Figure 2) or to improve nanoantennas radiation emission efficiency [19].

2. Materials and Methods

Polycrystalline diamond substrates with a lateral dimension of 10 × 10 mm² and a thickness of 250 µm were purchased by Element Six Ltd. The sample surface was treated with a laser Spectra-Physics Spitfire Pro 100 F 1 K XP (pulse duration about 100 fs, wavelength 800 nm, repetition rate up to 1 kHz). The horizontally polarized laser beam was focused perpendicularly to the diamond plate surface in a high-vacuum chamber (pressure < 10⁻⁷ mbar) with a spot size on the sample of 150 µm diameter. The pulse energy was kept constant at 3.6 mJ/pulse, corresponding to a power density one hundred times larger than the diamond damage threshold (200 MW/µm²). The black diamond was created with double-nanotexturing with a total accumulated laser fluence of 6.25 kJ/cm² on the diamond surface along the X-Y directions, and was obtained following the two-step bustrophedic procedure described in [37], by means of a computer-controlled micrometric X-Y translational stage. In the first step, a one-dimensional LIPSS structure was defined along the X direction. In the second step, the LIPSS structure was redefined along the Y direction thanks to the incubation effect guaranteeing a higher laser light absorption than the sample before the first treatment step. After the laser treatment, both samples were immersed in a strongly oxidizing solution (H₂SO₄:HClO₄:HNO₃ in the 1:1:1 ratio) for 15 min at the boiling point in order to remove any debris, as well as non-diamond carbon phases. Finally, the samples were cleaned by ultrasound sonication and rinsed in deionized water.

Surface morphology was studied by a Zeiss LEO Supra 35 field emission gun scanning electron microscope (FEG SEM). Secondary electron micrographs were acquired with an in-lens detector. The electron beam energy was set to 10 keV.

Scanning Kelvin probe measurements of the central part of the black diamond sample were performed by an OmegaScope platform integrated into a LabRAM HR Evolution Raman microscope (HORIBA Ltd., Kyoto, Japan). Morphology and KPFM imaging were performed in tapping mode, setting the operational amplitude at 20 nm, and using a conductive Cr-Au coating silicon pyramidal tip (MikroMasch HQ: NSC14/Al BS, Wetzlar, Germany) with a characteristic radius of ~35 nm and about 110 kHz as the resonance frequency. The scan rate was fixed at 0.5 Hz. All the microscopies were acquired in kelvin probe mode and analyzed using the AIST-NT SPM control software.

Raman spectroscopy was performed in back-scattering geometry by using a Dilor XY triple spectrometer (1 cm⁻¹ resolution) equipped with an Ar+ laser (514.5 nm), a cooled CCD detector, and an adapted Olympus microscope arranged in confocal mode. The spot size was 2 mm, and 20 intensity points were collected to provide a suitable statistic. Line widths and intensities were determined using the Thermo-Grams Suite v9.2 software.

After Raman spectroscopy experiments, a windowed stainless-steel sheet that was 50 µm thick was used as shadow-mask for thermal evaporation of Aluminum contacts on each sample. The distance between contacts was fixed at 1 mm for two parallel rectangles (see Figure 3).

Conductivity measurements up to 500 V in both bias voltage polarities were performed in a vacuum chamber with substrate temperature regulation performed by hot wires heating a copper hot-plate, on which the diamond plate is glued by silver paste to maximize the thermal flux. The current was acquired by a Keithley 487 Picoammeter/Voltage Source through two screwed molybdenum metal tips adhered in the central region of both Aluminum pads, both glued always by silver paste, by the application of a first ramp from 0 V → 500 V → 0 V (instrument maximum source voltage) and a second ramp with reversed electric field direction from 0 V → −500 V → 0 V.

3. Results

3.1. Morphological and Chemical Characterization

The sample analyzed in this paper is produced through a double-texturing technique with the target to maximize the largest total accumulated absorbed laser fluence, since the accumulated laser fluence of the second step was kept as large as the first: i.e., 50–50%. Indeed, the increase of absorptance after the first treatment is coupled with a decrease of the diamond ablation threshold. Therefore, the total accumulated absorbed laser fluence, defined as the product of the diamond absorptance at the laser wavelength and the delivered fluence, is characterized by a greater weight in the second treatment. As can be seen in Figure 4a, the diamond sample shows almost regular LIPSS, but is oriented according to the direction expected from the polarization direction of the second treatment, even if it is possible to observe perpendicular shallow lines deriving from the first treatment light polarization direction. Clearly, even if we did not perform an SEM session in the middle of the two steps of the BD fabrication, it is reasonable to assume that LIPSS formed horizontally after the first treatment, but thanks to the largest absorptance of the diamond surface during the second step, they were swept away during the formation of the new LIPSS in the second treatment. Traces of the LIPSS formed during the first step are evident in the inset zoom. The LIPSS periodicity related to the ratio of laser wavelength and the material refractive index does not appreciably change between the first texturing (First-step LIPSS) and the second one (Second-step LIPSS). A strong reduction of the dielectric permittivity may change the LIPSS periodicity, however:

• The ε_r change is supposed to be almost 2D. Therefore, it is unable to significantly contribute to the overall laser-diamond substrate interaction of the second texturing.

• The first treatment is performed at half of the standard power to produce a high-quality black-diamond surface. Therefore, the superficial ε_r change is also limited in value.

• The ε_r reduction is further enhanced by the two-stage texturing due to the incubation effect (see Appendix B).

In Figure 4b we show the Raman spectrum performed after the debris surface cleaning (green curve) showing considerable tensile stress in the diamond lattice. On the other hand, the annealing, performed at a temperature slightly higher than 500 °C for more than 24 h, revealed that high-temperature conditions influence reassembling of the diamond lattice and stress release in the surface layer. Therefore, it was decided to first bring the sample to the maximum temperature and second analyze its electrical behavior during the cooling phase (see Section 3.2).

A main point regarding black diamond nanostructured surface is the relation between the trap’s disposition and the periodicity of the laser-induced surface nanowaves. As stated in the introduction, black diamond samples with high ordering of the surface waves exhibit a superior conductivity on behalf of samples that absorbed even more laser energy but suffered a complete reorganization of the surface waves like the one analyzed in this paper. Therefore, it is reasonable to suppose that the surface waves are characterized by distinct traps depth levels. We performed surface Kelvin-Probe measurements, widely employed to study defects’ energy levels in organic electronic materials. Indeed, as shown in Figure 5, the higher black diamond zones of the laser-induced nanowaves, estimated to be at about 45°, have a distinctly more positive surface potential than the lower ones, with edge potentials reaching values of a hundred millivolts, which is more than twice the RT thermal voltage. Positive potentials mean that electrons are trapped in those regions [38,39] through localized defects with energy close to the conduction band.

3.2. Electrical Measurements

We performed high electric-field J(E) measurements (up to 5 kV/cm) with 23 steps of temperature in both the polarities of the applied electric field, with typical results of the first voltage ramp (0 V → 500 V) shown in Figure 6a–d (four I/V curves at different temperatures with relative Matlab fitting lines, see next paragraph), while opposite polarity electric field results are shown in Appendix A. Usually, it is possible to reveal the coefficient β of the FP effect by plotting the logarithmic of the current density ($\log J$) vs. the square root of the voltage (√V) and extracting its slope S = β/√d, where d is the distance between the measurement electrodes. However, in order to be more accurate, we decided to operate the coefficient extraction directly by using the Matlab “Curve Fitting” tool with the imposition of Equation (2) directly to the original data (see Figure 6). We noticed that the switching point between the two exponential regimes always falls in the 230 V–380 V range. Therefore, we applied three different fittings for each measurement:

• A.. A Linear Fitting (I = A_Offset + G_Lin*V), excluding point acquired at voltages > 100 V (red lines), where A_Offset is the little initial bias error of the instrument and G_Lin is the overall conductance of the sample at low voltages.

• B.. A First Exponential Fitting (I = G₁^FP*V*exp(β₁√V)), excluding the point acquired at voltages > 230 V (blue lines), where G₁^FP represents the sample linear conductance to be combined with the exponential contribution and β₁ represents the FP coefficient divided by the contact distance square root.

• C.. A Second Exponential Fitting (I = G₂^FP*V*exp(β₂√V)), considering all the experimental points (green lines), where G₂^FP and β₂ have the same meaning of G₁^FP and β₁.

It is possible to notice that the difference between the two exponential regimes changes with temperature. Such a difference is evident at room temperature (RT), it increases in the range 50 °C → 200 °C and it diminishes down to vanish at the highest available temperature, close to 500 °C. The extracted parameters (fits B and C of the previous list) are reported in Figure 7a (conductance) and Figure 7b (FP exponential coefficient) as a function of the 23 temperatures of measurement. Such results show that conductance increases exponentially with temperature as it is expected for thermal activation of charges trapped in shallow defects [40] and that the first exponential regime is characterized by a lower conductance value. In Figure 7b, it is instead evident that the exponential coefficient is slightly higher in first low-voltage range. On the other hand, the curves of Figure A2 regarding this factor show a reversed situation with cases of temperatures showing exactly the same exponential coefficient extracted but with values that always range from about 0.6 at the lowest temperatures to less than 0.3 at the highest ones.

4. Discussion

As expressed in Equation (1), the behavior of the current-electric field curve at a given temperature is determined only by the optical dielectric constant of the insulating film. In the vast majority of articles, the identification of the FP effect as the dominant current mechanism is made by plotting the sample conductance as a function of the square root of the applied electric field, leading to a straight line from which the value of the high-frequency dielectric constant is calculated [13]. In diamond, the frequency of the extracted k should correspond to the timing related to the fast electron transport in the conduction band (ε_r~5.5) and not to the static dielectric response (ε_r~7.8). In fact, in literature relative to the FP effect for diamond samples, it is typically found that a value for the exponential coefficient of about 3.2 × 10⁻⁴ eVcm/V, properly corresponding to an ε_r = 5.5 [41,42], even though higher values have been found under irradiation [43].

By substituting in Equation (1) the inherent physical constants, we can evaluate @ RT the following relation between the absolute value of dielectric permittivity and the parameter (β_1,2) extracted by the I/V curves fitting:

(2)${\beta }_{1,2}\left({\epsilon }_r=1,\mathrm{RT}\right)=92.7*{10}^{-3} \left[\frac{1}{\sqrt{V}}\right]; \to {\epsilon }_r\left(T\right)={\left(\frac{92.7*{10}^{-3}}{{\beta }_{1,2}}*\frac{300K}{T\left(K\right)}\right)}^2;$

where a larger measured slope corresponds to lower values of dielectric permittivity, according to the correction due to temperature variation.

If we consider the results shown in the previous section, we immediately understand that the exponential increase of the current with the electric field is much larger than the one calculated for Vacuum/Air. In fact, the largest value of β_1,2, measured just close to RT, corresponds to a value of 0.6, almost one decade larger than the one corresponding at ${\epsilon }_r=1$, and therefore the measured slope corresponds to a relative dielectric permittivity close to 10⁻².

In Figure 8a we report the behavior of the slope extracted from the various I/V curves analyses as a function of the temperature at which the measurement was conducted. While the slopes shown in Figure 7b decrease with temperature, this does not correspond to a constant value of dielectric permittivity since the variation due to temperature according to Equation (2) is much larger. Therefore, as reported in the adjacent Figure 8b, the extracted value of the relative dielectric permittivity in reality halves to values close to the hundredth.

FP conduction was observed in the non-linear region at high voltages for many diamond samples regardless of deposition processes [44,45]. Amplification of the FP effect was obtained through X-ray irradiation in [43], where it was found for the exponential coefficient that a parabolic dependence on the density of occupied traps such to reach more than quadruple values in a specific electric field range (2 kV/cm → 7 kV/cm). Such behavior was attributed to the intense trap occupation by the promotion of electrons from valence bands, causing coulomb interaction between the trap centers with subsequent barrier lowering with an estimated mean distance (λ_traps) of the electron-filled trap states of about 70 nm. The energy reduction due to the coulombic interaction of close traps additive to the electric field contribution may be related to their distance (a) with the following equation based on the Vollmann model [46]:

(3)$F{P}_{energy}\left(T\right)=q*{\Phi }_{th}{}^2*{\beta }_{1,2}{}^2*\frac{d}{a} ;$

where ${\Phi }_{th}$ is the thermal potential (see Figure 9). Since we are not observing a transition from Poole to Frenekl-Poole conduction regimes, it is not possible to extrapolate the real traps’ distance from the I/V curves, but we can consider that, given the high extrapolated beta coefficient, FP electron emission happens at distances of hundreds of nm, thus corresponding to the spatial gap between laser induced nanowaves (see Figure 4a and Figure 5b). Obviously, such value is just a tentative estimation and new, more complex transport experiments should be performed in order to correctly evaluate it. We can also note that the FP effect, as shown in Figure 8, increases with temperature, allowing us to estimate the lower values of the dielectric permittivity. Such behavior perfectly agrees with the increase of the conductance due to the thermal activation of electrons in the conduction band, which interacts at an increased rate with the empty trap center positively charged, and enhances the traps filling percentage at the origin of the FP effect reduction in BD.

5. Conclusions

Defected surfaces have plenty of new physics phenomena to be analyzed, even more if the defects’ geometry follows a specific regular pattern, as in the case of femtosecond laser surface-induced LIPSS. In this paper, we have shown that the electrical conduction of double nanotextured BD films is strongly influenced by the applied electric field. By applying the FP model to the I/V measurements, it was possible to reveal a dielectric permittivity near zero for the BD sample. The effect of the electric field over the trapped charges (FP effect) is appreciable just after it reaches 1 kV/cm, which is, to our knowledge, far lower than any other reported value over the presence of the FP effect in diamonds. The dielectric permittivity furthermore reduces its absolute value up to reach values lower than 0.01 for temperatures between 400 °C and 450 °C with electric fields lower than 2.3 kV/cm.

Such a result demonstrates that it is possible to reach quasi-ENZ polarizability conditions on diamond surfaces thanks to laser surface nanostructuring, an approach that is very likely to be extendable to other materials’ surfaces, widening the class of ENZ materials. This new property found on diamonds demonstrates again that it should be considered even more than other materials for the continued developments of novel functional photonic devices, since its properties are unique in nature and it may combine very low surface polarizability to other quantum properties related to its optically-active defect centers.

Author Contributions

Conceptualization, A.O., S.P. and M.M.; methodology, A.O. and M.M.; software, A.O., M.M. and A.B.; validation, A.O. and M.M.; formal analysis, A.O. and E.B.; investigation, A.O., E.B., A.B., M.G., M.C.R., M.M., R.P., V.S., V.V. and S.O.; resources, D.M.T. and R.P.; data curation, A.O. and E.B.; writing—original draft preparation, A.O.; writing—review and editing, All; visualization, All.; supervision, S.S. and D.B.; project administration, D.M.T.; funding acquisition, D.B. and D.M.T. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The data presented in the current work are available on request from the corresponding authors.

Acknowledgments

The authors thank Alessio Monti and Giuseppe Prestopino for their useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. Schematic of the “Black Diamond” fabrication process.

Figure 2. Sketch of the possible combination of quantum blocks on diamond substrates through spatially controlled laser induced ENZ waveguides.

Figure 3. Sketch of the current measurements experimental setup.

Figure 4. Surface characterization of black diamond samples (a) SEM image of the surface; (b) Raman analysis of the diamond peak at 1331.2 cm−1 before (green curve) and after (red curve) of the annealing at temperature slightly higher than 500 °C for more than 24 h. Black line corresponds to natural diamond Raman peak.

Figure 5. Surface potential maps of black diamond laser induced periodic surface nanostructures and their role as electron traps. (a) AFM image of the surface; (b) SKPM surface potential images.

Figure 6. I/V curves relative to the positive voltage ramp with corresponding Matlab fitting for the linear part (red line), the first exponential part up to 230V (blue curve) and the total exponential behavior up to maximum voltage of 500 V (green curve). (a) measurement at about 50 °C; (b) measurement at about 200 °C; (c) measurement at about 350 °C; (d) measurement at about 500 °C.

Figure 7. Frenkel-Poole- parameters as extracted from the I/V curves: (a) conductance for V < 230 V (blue curve, number 1 subscript in the legend) and for V < 500 V (green curve, number 2 subscript in the legend); (b) Extraction of the exponential coefficient for V < 230 V (blue curve, number 1 subscript in the legend) and for V < 500 V (green curve, number 2 subscript in the legend).

Figure 8. Relative Dielectric Permittivity as extracted by the I/V curves reported in Figure 6. (a) Extraction of parameter with V < 230 V; (b) Extraction of parameter with V > 300 V.

Figure 9. (a) FP energy in standard materials with high density filled traps with reciprocal coulombic interaction with highlighted barrier lowering. (b) FP energy in quasi-ENZ BD surface.

Appendix A. Opposite Direction Voltage Ramp Measurements

We proceeded to perform accurate I/V curves measurements with 23 steps of temperature both in the positive and negative voltage direction. In Figure A1 we show four measurements with relative Matlab fitting corresponding to the same temperatures of Figure 4 but for the opposite sign voltage ramp. The sample is symmetric, so it does not meaningfully stress which contact is positive or negative. However, due to possible charge accumulation at the interfaces, it is important to stress that these values were collected with a voltage ramp successive to the first one shown in Figure 4.

Figure A1. I/V curves relative to the negative voltage ramp with corresponding Matlab fitting for the linear part (red line), the first exponential part up to 230 V (blue curve) and the total exponential behavior up to maximum voltage of 500 V (green curve). (a) n. 23 measurement at about 50 °C; (b) n. 15 measurement at about 200 °C; (c) n8 measurement at about 350 °C; (d) n. 1 measurement at about 500 °C.

We can then plot the obtained parameter of first FP conductance (red curve), first FP exponential coefficient β₁ (blue curve), and second FP exponential coefficient β₂ (green curve) as a function of the 23 measured temperature steps. Such results are reported in Figure 5a (linear conductivity) and Figure 5b (FPs).

Figure A2. Frenkel-Poole parameters as extracted by the I/V curves reported in Figure A1. (a) Extraction of Conductance for V < 230 V (blue curve, number 1 subscript in the legend) and for V < 500 V (green curve, number 2 subscript in the legend); (b) Extraction of Exponential Coefficient for V < 230 V (blue curve, number 1 subscript in the legend) and for V < 500 V (green curve, number 2 subscript in the legend).

As a last element, we show in Figure A3 the curves of the linear conductance as of the fitting case (A), where it is possible to observe the shallow exponential behavior with BD activation energies much lower that the untreated diamond substrate as reported previously [11].

Figure A3. Linear conductance of BD sample as extracted by the I/V curves reported in Figure 4 and Figure A1. (a) Extraction of sample conductance with V < 100 V for the first voltage ramp (b) Extraction of sample conductance with V < 100 V for the second voltage ramp with opposite biasing.

Appendix B. Analysis of Black-Diamond Sample with Reduced Total Accumulated Laser Fluence

In order to look for a trend in the dielectric permittivity reduction in the produced black diamonds, we analyzed a sample with the nanostructured surface in a single treatment (1T-BD). therefore, without the benefit of the increased optical energy absorption of the first laser treatment (incubation effect). The analysis of the FP effect according to temperature produced a similar result to what is shown in the main text, with a low dielectric permittivity decreasing with temperature. The absolute value of ε_r (Figure A4b), as extrapolated from the slope of the I vs. √V graphs (Figure A4a), is again more than fifty times lower than in bulk diamond:$${\epsilon }_r\left(1T-BD\right)<0.1 Vs {\epsilon }_r\left(Diamond\right)=5.5$$

On the other hand, we may note that the ε_r(1T − BD) minimum value is around 6 × 10⁻², which is about six times the minimum value shown in the main text for ε_r(2T − BD). This ratio perfectly agrees with a lower density of electron traps.

Figure A4. (a) Current of 1T-BD sample at the higher experimental temperatures. (b) Extraction of of 1T-BD Relative Permittivity according to Temperature by using Equation (3).

References

1. Calvani, P.; Bellucci, A.; Girolami, M.; Orlando, S.; Valentini, V.; Lettino, A.; Trucchi, D.M. Optical properties of femtosecond laser-treated diamond. Appl. Phys. A Mater. Sci. Proc.; 2014; 117, pp. 25-29. [DOI: https://dx.doi.org/10.1007/s00339-014-8311-9]

2. Girolami, M.; Bellucci, A.; Calvani, P.; Trucchi, D.M. Large single-crystal diamond substrates for ionizing radiation detection. Phys. Status Solidi; 2016; 213, pp. 2634-2640. [DOI: https://dx.doi.org/10.1002/pssa.201600262]

3. Trucchi, D.M.; Melosh, N.A. Electron-emission materials: Advances, applications, and models. MRS Bull.; 2017; 42, pp. 488-492. [DOI: https://dx.doi.org/10.1557/mrs.2017.142]

4. Girolami, M.; Bellucci, A.; Mastellone, M.; Orlando, S.; Orsini, A.; Salvatori, S.; Trucchi, D.M. Transport properties of photogenerated charge carriers in black diamond films. Ceram. Int.; 2018; 45, pp. 9544-9547. [DOI: https://dx.doi.org/10.1016/j.ceramint.2018.09.313]

5. Trucchi, D.M.; Bellucci, A.; Girolami, M.; Calvani, P.; Cappelli, E.; Orlando, S.; Polini, R.; Silvestroni, L.; Sciti, D.; Kribus, A. Solar Thermionic-Thermoelectric Generator (ST2G): Concept, Materials Engineering, and Prototype Demonstration. Adv. Energy Mater.; 2018; 8, pp. 1-14. [DOI: https://dx.doi.org/10.1002/aenm.201802310]

6. Granados, E.; Martinez-Calderon, M.; Gomez, M.; Rodriguez, A.; Olaizola, S.M. Photonic structures in diamond based on femtosecond UV laser induced periodic surface structuring (LIPSS). Opt. Express; 2017; 25, 15330. [DOI: https://dx.doi.org/10.1364/OE.25.015330]

7. Martínez-Calderon, M.; Azkona, J.J.; Casquero, N.; Rodríguez, A.; Domke, M.; Gómez-Aranzadi, M.; Olaizola, S.M.; Granados, E. Tailoring diamond’s optical properties via direct femtosecond laser nanostructuring. Sci. Rep.; 2018; 8, 14262. [DOI: https://dx.doi.org/10.1038/s41598-018-32520-0]

8. Salvatori, S.; Rossi, M.C.; Conte, G.; Kononenko, T.; Komlenok, M.; Khomich, A.; Ralchenko, V.; Konov, V.; Provatas, G.; Jaksic, M. Diamond Detector with Laser-Formed Buried Graphitic Electrodes: Micron-Scale Mapping of Stress and Charge Collection Efficiency. IEEE Sens. J.; 2019; 19, pp. 11908-11917. [DOI: https://dx.doi.org/10.1109/JSEN.2019.2939618]

9. Khomich, A.A.; Ashikkalieva, K.K.; Bolshakov, A.P.; Kononenko, T.V.; Ralchenko, V.G.; Konov, V.I.; Oliva, P.; Conte, G.; Salvatori, S. Very long laser-induced graphitic pillars buried in single-crystal CVD-diamond for 3D detectors realization. Diam. Relat. Mater.; 2018; 90, pp. 84-92. [DOI: https://dx.doi.org/10.1016/j.diamond.2018.10.006]

10. Girolami, M.; Bellucci, A.; Mastellone, M.; Orlando, S.; Valentini, V.; Montereali, R.M.; Vincenti, M.A.; Polini, R.; Trucchi, D.M. Optical characterization of double-nanotextured black diamond films. Carbon N. Y.; 2018; 138, pp. 384-389. [DOI: https://dx.doi.org/10.1016/j.carbon.2018.07.055]

11. Orsini, A.; Bellucci, A.; Girolami, M.; Mastellone, M.; Orlando, S.; Prestopino, G.; Valentini, V.; Salvatori, S.; Trucchi, D.M. Electrical conductivity of double textured black diamond films from RT to 800 K. Diam. Relat. Mater.; 2019; 93, pp. 1-7. [DOI: https://dx.doi.org/10.1016/j.diamond.2019.01.011]

12. Orsini, A.; Barettin, D.; Ercoli, F.; Rossi, M.C.; Pettinato, S.; Salvatori, S.; Mezzi, A.; Polini, R.; Bellucci, A.; Mastellone, M. et al. Charge Transport Mechanisms of Black Diamond at Cryogenic Temperatures. Nanomaterials; 2022; 12, 2253. [DOI: https://dx.doi.org/10.3390/nano12132253]

13. Jeong, D.S.; Park, H.B.; Hwang, C.S. Reasons for obtaining an optical dielectric constant from the Poole-Frenkel conduction behavior of atomic-layer-deposited HfO2 films. Appl. Phys. Lett.; 2005; 86, pp. 1-3. [DOI: https://dx.doi.org/10.1063/1.1865326]

14. Jahani, S.; Zhao, H.; Jacob, Z. Switching Purcell effect with nonlinear epsilon-near-zero media. Appl. Phys. Lett.; 2018; 113, 21103. [DOI: https://dx.doi.org/10.1063/1.5030023]

15. Alù, A.; Silveirinha, M.G.; Salandrino, A.; Engheta, N. Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern. Phys. Rev. B Condens. Matter Mater. Phys.; 2007; 75, pp. 4-15. [DOI: https://dx.doi.org/10.1103/PhysRevB.75.155410]

16. Monti, A.; Alù, A.; Toscano, A.; Bilotti, F. Metasurface-based anti-reflection coatings at optical frequencies. J. Opt.; 2018; 20, 055001. [DOI: https://dx.doi.org/10.1088/2040-8986/aab886]

17. Fruhnert, M.; Monti, A.; Fernandez-Corbaton, I.; Alù, A.; Toscano, A.; Bilotti, F.; Rockstuhl, C. Tunable scattering cancellation cloak with plasmonic ellipsoids in the visible. Phys. Rev. B; 2016; 93, 245127. [DOI: https://dx.doi.org/10.1103/PhysRevB.93.245127]

18. Monti, A.; Alù, A.; Toscano, A.; Bilotti, F. Optical invisibility through metasurfaces made of plasmonic nanoparticles. J. Appl. Phys.; 2015; 117, 123103. [DOI: https://dx.doi.org/10.1063/1.4916257]

19. Kim, J.; Dutta, A.; Naik, G.V.; Giles, A.J.; Bezares, F.J.; Ellis, C.T.; Tischler, J.G.; Mahmoud, A.M.; Caglayan, H.; Glembocki, O.J. et al. Role of epsilon-near-zero substrates in the optical response of plasmonic antennas. Optica; 2016; 3, 339. [DOI: https://dx.doi.org/10.1364/OPTICA.3.000339]

20. Liu, C.; Pang, K.; Manukyan, K.; Reshef, O.; Zhou, Y.; Patrow, J.; Pennathur, A.; Song, H.; Zhao, Z.; Zhang, R. et al. Resonance Splitting and Enhanced Optical Nonlinearities in ITO-Based Epsilon-Near-Zero Metasurface with Cross-Shaped Nanoantennas. Proceedings of the 2019 Conference on Lasers and Electro-Optics (CLEO); Munich, Germany, 23–27 June 2019; pp. 1-2.

21. Rocco, D.; De Angelis, C.; de Ceglia, D.; Carletti, L.; Scalora, M.; Vincenti, M.A. Dielectric nanoantennas on epsilon-near-zero substrates: Impact of losses on second order nonlinear processes. Opt. Commun.; 2020; 456, 124570. [DOI: https://dx.doi.org/10.1016/j.optcom.2019.124570]

22. Hasan, M.R.; Hellesø, O.G. Dielectric optical nanoantennas. Nanotechnology; 2021; 32, 202001. [DOI: https://dx.doi.org/10.1088/1361-6528/abdceb] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33461187]

23. Zalogina, A.S.; Savelev, R.S.; Ushakova, E.V.; Zograf, G.P.; Komissarenko, F.E.; Milichko, V.A.; Makarov, S.V.; Zuev, D.A.; Shadrivov, I.V. Purcell effect in active diamond nanoantennas. Nanoscale; 2018; 10, pp. 8721-8727. [DOI: https://dx.doi.org/10.1039/C7NR07953B] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29701731]

24. Chen, Y.-C.; Salter, P.S.; Knauer, S.; Weng, L.; Frangeskou, A.C.; Stephen, C.J.; Ishmael, S.N.; Dolan, P.R.; Johnson, S.; Green, B.L. et al. Laser writing of coherent colour centres in diamond. Nat. Photonics; 2017; 11, pp. 77-80. [DOI: https://dx.doi.org/10.1038/nphoton.2016.234]

25. Wu, J.; Xie, Z.T.; Sha, Y.; Fu, H.Y.; Li, Q. Epsilon-near-zero photonics: Infinite potentials. Photonics Res.; 2021; 9, 1616. [DOI: https://dx.doi.org/10.1364/PRJ.427246]

26. Kinsey, N.; DeVault, C.; Boltasseva, A.; Shalaev, V.M. Near-zero-index materials for photonics. Nat. Rev. Mater.; 2019; 4, pp. 742-760. [DOI: https://dx.doi.org/10.1038/s41578-019-0133-0]

27. Reshef, O.; De Leon, I.; Alam, M.Z.; Boyd, R.W. Nonlinear optical effects in epsilon-near-zero media. Nat. Rev. Mater.; 2019; 4, pp. 535-551. [DOI: https://dx.doi.org/10.1038/s41578-019-0120-5]

28. Niu, X.; Hu, X.; Chu, S.; Gong, Q. Epsilon-Near-Zero Photonics: A New Platform for Integrated Devices. Adv. Opt. Mater.; 2018; 6, 1701292. [DOI: https://dx.doi.org/10.1002/adom.201701292]

29. Pasquarelli, A.; Andrilli, L.H.S.; Bolean, M.; Reis Ferreira, C.; Cruz, M.A.E.; de Oliveira, F.A.; Ramos, A.P.; Millán, J.L.; Bottini, M.; Ciancaglini, P. Ultrasensitive Diamond Microelectrode Application in the Detection of Ca2+ Transport by AnnexinA5-Containing Nanostructured Liposomes. Biosensors; 2022; 12, 525. [DOI: https://dx.doi.org/10.3390/bios12070525]

30. Xie, H.; Li, X.; Luo, G.; Niu, Y.; Zou, R.; Yin, C.; Huang, S.; Sun, W.; Li, G. Nano-diamond modified electrode for the investigation on direct electrochemistry and electrocatalytic behavior of myoglobin. Diam. Relat. Mater.; 2019; 97, 107453. [DOI: https://dx.doi.org/10.1016/j.diamond.2019.107453]

31. Feng, Z.; Gao, N.; Liu, J.; Li, H. Boron-doped diamond electrochemical aptasensors for trace aflatoxin B1 detection. Anal. Chim. Acta; 2020; 1122, pp. 70-75. [DOI: https://dx.doi.org/10.1016/j.aca.2020.04.062]

32. Yang, K.-H.; Narayan, R.J. Biocompatibility and functionalization of diamond for neural applications. Curr. Opin. Biomed. Eng.; 2019; 10, pp. 60-68. [DOI: https://dx.doi.org/10.1016/j.cobme.2019.03.002]

33. Pettinato, S.; Girolami, M.; Olivieri, R.; Stravato, A.; Caruso, C.; Salvatori, S. Time-Resolved Dosimetry of Pulsed Photon Beams for Radiotherapy Based on Diamond Detector. IEEE Sens. J.; 2022; 22, pp. 12348-12356. [DOI: https://dx.doi.org/10.1109/JSEN.2022.3173892]

34. Pettinato, S.; Orsini, A.; Girolami, M.; Trucchi, M.D.; Rossi, C.M.; Salvatori, S. A High-Precision Gated Integrator for Repetitive Pulsed Signals Acquisition. Electronics; 2019; 8, 1231. [DOI: https://dx.doi.org/10.3390/electronics8111231]

35. Evans, R.E.; Bhaskar, M.K.; Sukachev, D.D.; Nguyen, C.T.; Sipahigil, A.; Burek, M.J.; Machielse, B.; Zhang, G.H.; Zibrov, A.S.; Bielejec, E. et al. Photon-mediated interactions between quantum emitters in a diamond nanocavity. Science; 2018; 362, pp. 662-665. [DOI: https://dx.doi.org/10.1126/science.aau4691]

36. Liu, K.; Zhang, S.; Ralchenko, V.; Qiao, P.; Zhao, J.; Shu, G.; Yang, L.; Han, J.; Dai, B.; Zhu, J. Tailoring of Typical Color Centers in Diamond for Photonics. Adv. Mater.; 2021; 33, 2000891. [DOI: https://dx.doi.org/10.1002/adma.202000891]

37. Bellucci, A.; Calvani, P.; Girolami, M.; Orlando, S.; Polini, R.; Trucchi, D.M. Optimization of black diamond films for solar energy conversion. Appl. Surf. Sci.; 2016; 380, pp. 8-11. [DOI: https://dx.doi.org/10.1016/j.apsusc.2016.02.107]

38. Tello, M.; Chiesa, M.; Duffy, C.M.; Sirringhaus, H. Charge Trapping in Intergrain Regions of Pentacene Thin Film Transistors. Adv. Funct. Mater.; 2008; 18, pp. 3907-3913. [DOI: https://dx.doi.org/10.1002/adfm.200800009]

39. He, T.; Wu, Y.; D’Avino, G.; Schmidt, E.; Stolte, M.; Cornil, J.; Beljonne, D.; Ruden, P.P.; Würthner, F.; Frisbie, C.D. Crystal step edges can trap electrons on the surfaces of n-type organic semiconductors. Nat. Commun.; 2018; 9, 2141. [DOI: https://dx.doi.org/10.1038/s41467-018-04479-z]

40. Chakraborty, A.; Orsini, A.; Kar, J.P.; Gatta, F.; Khan, U.; Falconi, C. Ultra-efficient thermo-convective solution-growth of vertically aligned ZnO nanowires. Nano Energy; 2022; 97, 107167. [DOI: https://dx.doi.org/10.1016/j.nanoen.2022.107167]

41. Gonon, P.; Prawer, S.; Boiko, Y.; Jamieson, D.N. Electrical conduction in polycrystalline diamond and the effects of UV irradiation. Diam. Relat. Mater.; 1997; 6, pp. 860-864. [DOI: https://dx.doi.org/10.1016/S0925-9635(96)00675-9]

42. Brozdnichenko, A.N.; Dolgintsev, D.M.; Castro, R.A. The dielectric properties of the diamond-like films grown by ion-plasma method. J. Phys. Conf. Ser.; 2014; 572, 012025. [DOI: https://dx.doi.org/10.1088/1742-6596/572/1/012025]

43. Girolami, M.; Bellucci, A.; Calvani, P.; Flammini, R.; Trucchi, D.M. Radiation-assisted Frenkel-Poole transport in single-crystal diamond. Appl. Phys. Lett.; 2013; 103, 083502. [DOI: https://dx.doi.org/10.1063/1.4818904]

44. Gonon, P.; Boiko, Y.; Prawer, S.; Jamieson, D. Poole–Frenkel conduction in polycrystalline diamond. J. Appl. Phys.; 1996; 79, pp. 3778-3780. [DOI: https://dx.doi.org/10.1063/1.361212]

45. Chou, H.-T.; Lee, C.-C.; Sun, C.-H. Conduction mechanisms in undoped polycrystalline diamond films. Optoelectronic Materials and Devices II; Su, Y.-K.; Bhattacharya, P. SPIE: Bellingham, WA, USA, 2000; Volume 4078, pp. 680-688.

46. Pillonnet, A.; Ongaro, R. Synthetic representations of the Poole Frenkel (PF) and Poole regimes. Rev. Phys. Appl.; 1990; 25, pp. 229-242. [DOI: https://dx.doi.org/10.1051/rphysap:01990002502022900]