Scripts for quantifying the effect of diamond nano-pillars on the fluorescence of NV centers

Summary

Scripts and data can be used to reproduce and build on the numerical results published under the title: "Optical and Spin Properties of NV Center Ensembles in Diamond Nano-Pillars" by Kseniia Volkova, Julia Heupel, Sergei Trofimov, Fridtjof Betz, Rémi Colom, Rowan W. MacQueen, Sapida Akhundzada, Meike Reginka, Arno Ehresmann, Johann P. Reithmaier, Sven Burger, Cyril Popov, and Boris Naydenov (Nanomaterials 12(9), 1516, 2022).

Method

The dipole emitters are assumed to be distributed uniformly 30 nm below the top surface of the nano-pillars. They are first integrated with a trapezoidal rule along the azimuth (because of the periodicity this results in a geometrical convergence) and with a 15 point Gauss-Kronrod quadrature rule in radial direction.

The main source of error results from the dipole positions being integrated only from 0 to R - min_dist, as it is challenging to model a dipole emitter located only few nanometers from the curved material interface. Further numerical parameters can be adjusted in the input files for JCMsuite. Both, a 3D setup and a 2D setup are provided. The letter exploits the rotational symmetry which results in a smaller memory footprint. Yet, as the distance of the dipole from the symmetry axis increases, many Fourier components are required which leads to long computation times.

For further quantitative studies we propose the 3D setup that is the default in the script 'integration.m', which allows to integrate closer to the side walls without increasing the costs. Furthermore in the second data set, shipped together with the data published in the paper, the height has been kept constant. In the paper the height has been chosen according to the fabricated samples. The 90° angle has been assigned to the [111] samples and the correspondingt height of 1400 nm and the 35.3° angle was assigned to the [100] samples and a height of 2200 nm.

Structure

The directories scattering2D, scattering3D and scatteringFlat contain input files for JCMsuite. The script 'integration.m' can be used to produce new data. With 'plotresults.m' you can either plot the results produced with 'integration.m' or those which were published in the related paper. Please note that the provided example produced with the script 'integration.m' differs from the published data, which has been computed with slightly different parameters.

Requirements

• JCMsuite 5.2.0
• Matlab R2019b

In order to produce new data, you must replace the corresponding place holders in the files by a path to your installation of JCMsuite. Free trial licenses are available, please refer to the homepage of JCMwave.