Dependence of the surface parameters of the NCD samples on the boron doping level.

*Values adjusted to the surface potential of gold.For the NDP measurements, the accuracy can achieve 5% (which is the precision of the boron atoms in the etalon). However, the precision of the NDP technique also depends on other parameters, e.g. the stability of the neutron beam intensity, identical geometry of the etalon and the measured sample, etc. Realistically, the NDP data can be routinely measured with accuracy of 10% in our case.For roughness, potential, phase and contact angle, the data is presented as Mean ± S.D. (Standard Deviation). In the case of roughness and AFM phase, each rms value was determined from 65 536 data points on each sample type. The mean and S.D. of rms values were calculated from 5 such measurements across the sample. In the case of surface potential, the mean and S.D. values were calculated from 65 536 measurements across each sample. The contact angle was calculated from fitting the curve of the water droplet, and the mean and S.D. values were calculated from 16 measurements for each sample type.Statistical Analysis: ANOVA, Student-Newman-Keuls Method. Statistical significance: I, II, III, IV: p≤0.05 compared to the group labelled with the same Roman number.For the room temperature electrical resistivity measurements, the accuracy is better than 1%.