## Abstract

Size-dependent quantization of energy spectrum of conducting electrons in solids leads to oscillating dependence of electronic properties on corresponding dimension(s). In conventional metals with typical energy Fermi *E*
_{F} ~ 1 eV and the charge carrier's effective masses *m** of the order of free electron mass *m*
_{0}, the quantum size phenomena provide noticeable impact only at nanometer scales. Here we experimentally demonstrate that in single-crystalline semimetal bismuth nanostructures the electronic conductivity non-monotonously decreases with reduction of the effective diameter. In samples grown along the particular crystallographic orientation the electronic conductivity abruptly increases at scales of about 50 nm due to metal-to-insulator transition mediated by the quantum confinement effect. The experimental findings are in reasonable agreement with theory predictions. The quantum-size phenomena should be taken into consideration to optimize operation of the next generation of ultra-small quantum nanoelectronic circuits.

## Introduction

With reduction of dimension(s) classic physics gradually degenerates into quantum. For example, energy separation between the *n*-th and the (*n+1*)-th energy levels of a particle with mass *m* confined to a space of characteristic dimension *a* is Δ*E*(*n*) *~*
*n*/*ma*
^{2} given *n* >> 1. This key quantum feature is in the origin of discrete energy spectra of atoms. In metals size quantization affects basically all electronic properties: conductivity, Hall and thermoelectric coefficients, reflectivity, chemical surface activity, etc.^{1} For a typical metal with effective mass *m** ≈ *m*
_{0} at realistic temperatures *k*
_{
B
}
*T* < Δ*E* one should expect quantum size phenomena to manifest at difficult-to-achieve scales *a*
*~* 1 nm*.* A solution is to study materials with *m**<<*m*
_{
0
} such as bismuth, antimony or *BiSb* alloys.^{2} Peculiarity of the energy spectrum of bismuth is the co-existence of two types of charges: 'light' L-electrons with *m*
_{e}
***<<*m*
_{0} and 'heavy' T-holes with *m*
_{h}
***
*~*
*m*
_{0} (Fig. 1a, b). The heavy holes 'stabilize' Fermi energy *E*
_{F} leading to oscillating dependence of transport properties vs. characteristic dimension, and pronounced decrease of electron resistivity due to metal-to-insulator transition at the threshold scale *a*
_{0} so that *E*
_{F}
*~*
*1/m*
_{e}
**a*
_{0}
^{2}. As bismuth energy spectrum is highly anisotropic (Fig. 1a), the size-dependent behavior of characteristic parameter, e.g. electric conductivity, depends on particular crystallographic orientation.

At early stages of experimental studies^{3, 4} 2D films were of primary attention, while later other low-dimensional bismuth structures have been investigated: nanowires,^{5, 6} microcylinders^{7} and point contacts.^{8, 9} Experiments with wide polycrystalline films^{3, 4} or bundles of nanowires^{6} can be interpreted only qualitatively as the contributions of individual grains/nanowires are averaged. STM experiments^{8, 9} have demonstrated clear size-dependent variation of electric conductivity. However, comparison of STM data with theory is ambiguous as neither the crystallographic orientation, nor the shape of an STM-pulled point contact are clearly defined.

Here we present experimental study of electric conductivity of single-crystalline bismuth nanostructures (or alternatively they can be called *nanorods*, given their relatively short length compared to width and thickness). The oscillatory dependence of resistance on sample cross-section is in good agreement with model predictions. The pronounced increase of resistivity below a certain scale can be interpreted as manifestation of size-mediated metal-to-insulator transition.

## Results

There have been numerous approaches of various complexity to describe energy spectrum of bulk bismuth.^{10} For our purposes it is quite sufficient to consider rather simple model.^{11, 12} For T-point holes (Fig.1a, b) the energy *E* vs. momentum *h*
** k** dependence is:

where the diagonal components of hole effective mass tensor (in units of free electron mass *m*
_{
0
}) are *m*
_{
x
}
^{h} = *m*
_{
y
}
^{h} = 0.059, *m*
_{
z
}
^{h} = 0.634. For L- point electrons (“+”) and L-point holes (“−”):

where the corresponding masses along the bisectrix axis are *m*
_{
x
}
^{e} = 0.00139, *m*
_{
y
}
^{e} = 0.291, *m*
_{
z
}
^{e} = 0.0071. There is the consensus between various literature sources that the energy gap in L-point is *E*
_{g} = 13.7 ± 0.1 meV, and the overlapping between T and L-bands is Δ = 38 meV (Fig. 1b).

Quantum size phenomena in bismuth have been considered already at early ears of 'bismuthology'.^{2, 13, 14} For interpretation of our data we utilize rather straightforward model based on Boltzmann equation.^{15} Following^{15} and considering bismuth nanowire of rectangular cross section *w·t*, where *w* is width and *t* is thickness, one gets the oscillatory behavior of electronic conductivity on corresponding dimensions:

where *S* is the volume concentration of scatterers each of strength *V*
_{
0
}, *e* is electron charge; outer summation stands for contribution of different sub-bands (i*, j*); \({\Lambda }_{m'n'}^{mn}\equiv (2+{\delta }_{mm'})(2+{\delta }_{nn'})\) with *δ*
_{
ab
} being Kronecker delta function; \({U}_{ij}^{e}\equiv 1-{(i/{r}_{w})}^{2}+({\mu }_{y}^{e}/{\mu }_{x}^{e})[1-{(j/{r}_{t})}^{2}]\); [*z*] stands for integer part of argument *z*. Parameters *r*
_{
w
} ≡ *w*/*w*
_{
0
} and *r*
_{
t
} ≡ *t*/*t*
_{
0
} are width and thickness of nanowire, normalized by the relevant dimensions *w*
_{0} = (*h*/2)(*M*
_{
x
}Δ_{
x
})^{1/2} and *t*
_{
0
} = (*h*/2)(*M*
_{
y
}Δ_{
y
})^{1/2} corresponding to metal-to-insulator transition: when gap opens between the lowest L-electron and the top-most T-hole bands (Fig. 1b). Here *M*
_{
x,y
}≡(*m*
^{e}
_{
x,y
} m^{h}
_{
x,y
}) / (*m*
^{e}
_{
x,y
} + *m*
^{h}
_{
x,y
}) and *μ*
^{e}
_{
x,y
}≡(Δ_{
x,y
}
*m*
^{h}
_{
x,y
}) / (*m*
^{e}
_{
x,y
} + *m*
^{h}
_{
x,y
}), Δ_{
x,y
} is the energy bands overlap due to confinement of charge carriers in corresponding direction *x* or *y*. Parameters *S* and *V*
_{0} are not known with necessary accuracy and are difficult to be determined from transport experiments. Contributing to absolute value of electron conductivity, they can be considered as adjustable parameters to fit *σ*
^{e}(*w*,*t*) data for thick wires *w*→∞ and *t*→∞ when size phenomena can be neglected. Similar to Eq. 3 expression for hole conductivity can be obtained and should be taken into consideration to calculate the total electric conductivity. Given much larger values of T-hole masses compared to L-electrons along bisectrix axis (Fig. 1a), we consider the hole conductivity *σ*
^{h}(*w*,*t*) as constant and neglect the size-dependent component.

The *σ*
^{e}(*w*,*t*) characteristic depends on particular orientation of the sample against crystallographic axes. Simple *σ*
^{e}(*w*,*t*) dependence appears for high-symmetric directions, while along an arbitrary direction the corresponding characteristics might be rather complicated containing secondary maxima originating from contributions of all non-equivalent pockets of bismuth Fermi surface (Fig. 1a). Note that even for the most favorable orientation along bisectrix axis corresponding to the lowest electron masses *m*
^{e}
_{
x
}, the actual *σ*
^{e}(*w,t*) dependence is affected not only by the cross section (*w·t*), but also by the actual values of width *w* and thickness *t* (Fig. 2).

The experiment consisted of bismuth thin films and nanostructures resistance *R* measurements as function of transverse dimensions *w* and *t* between sessions of ion milling. All samples, including the co-deposited and co-sputtered 2D films, demonstrated non-monotonous *R*(*t,w*) dependencies (Fig. 3). Note that electron transport data obtained on wide granular 2D films (Fig. 3a), similar to early studies,^{2, 3, 4} can be interpreted only qualitatively as orientation of electric current vs. crystallographic axes in each grain is different. For sufficiently long nanowires with length *l* noticeably exceeding the characteristic grain size, as studied in ref. 16, an agreement between experiment and theory is speculative. The problem originates from random impact of each of the *i*-th grain connected in series. Superposition of all those contributions results in complicated overall \(R=\sum _{i}{R}_{i}({d}_{x}^{i},{d}_{y}^{i},{d}_{z}^{i})\) dependence, difficult to analyze as the orientation of each grain and its dimensions {*d*
^{i}
_{
x
}, *d*
^{i}
_{
y
}, *d*
^{i}
_{
z
}} are not known.

Here we focus on quasi-1D structures *w,t*<<*l* ~ 1 μm with the 'body' formed from a single grain (Fig. 1c), where the *R*(*t,w*) dependencies can be modeled with reasonable accuracy. For the most favorable orientation of the sample along bisectrix axis (Fig. 1a), the dominating impact comes from L-electrons with minimal electron effective mass *m*
_{
x
}
^{e} ≈ 0.0014 *m*
_{0}. In this geometry the size-dependent contributions from two other electronic L-pockets and T-holes with much larger effective masses are noticeably weaker. Certainly, the absolute value of total electric conductivity is the sum of contributions from all charge carriers.

*R*(*t,w*) data with corresponding theory dependence is presented in Fig. 3b. The experiment shows both the non-monotonous variation of sample resistance vs. transverse dimensions *R*(*t,w*), and the pronounced increase of resistance below the threshold (*t*
_{
0
}
*·w*
_{
0
})^{1/2} ~ 50 nm. The non-monotonous variation of resistance *R*(*t,w)* originates from the size-dependent quantization of the energy spectrum (Fig. 1b), with minima of resistance (maxima of conductivity) corresponding to matching of the Fermi energy *E*
_{F} and the *n*-th energy level. The increase of the resistance below 50 nm is related to opening of the gap in the energy spectrum: the size-dependent transition from metallic to dielectric state. Yet another manifestation of the same effect has been recently observed as non-monotonous dependencies of power factor and thermoelectric figure of merit in Bi nanowires at sub-300 nm scales.^{17}

## Discussion

The agreement between experiment and theory is acceptable, but not perfect: the main features do match, while the 'fine structure' predicted by theory is not observed. Presumably several factors contribute to the discrepancy. First, finite temperature *T* = 4.2 K corresponds to broadening of each size-quantized energy level by ~0.5 meV, leading to smearing of features predicted by the essentially zero-temperature model.^{15} The observation is supported by data in Fig. 3, which indeed show much weaker *R*(*t,w*) dependence at room temperature compared to *T* = 4.2 K. Second, due to requirement to fabricate single-crystalline samples of minimal size suitable to be ion milled to reduce transverse dimensions, the measuring probes are of comparable dimensions as the 'body' of the sample (Fig. 1c). Hence, the electrodes cannot be considered as non-invasive and size quantization phenomena inside the probes might provide non-negligible contribution to measured signal. Third is the uncertainty in orientation of crystallographic axes with respect to direction of bias current (*x*). For data presented in Fig. 3b the misorientation between the bisectrix axis and the sample axis *x*, deduced from theory fits with Eq. 3, is just about few degrees. The effect of such a small misorientation was formally taken into account as small adjustment (increase) of the effective mass m_{x}
^{e} in Eq. 2. Note that the model from ref. 15 assumes all current lines being parallel to each other. That would be the case if the length of the sample *l* is much larger than the transverse dimensions *w* and *t*, which is not the case particularly for 'thick' samples before ion milling (Fig. 1c). Better agreement with experiment formally can be achieved assuming size-dependent contribution of more than one crystallite (e.g. from nearby regions of the probes) to the measured integral signal. However, as those contributions are not clearly defined, such 'improvement' of fits would be speculative.

To conclude, we fabricated and measured resistance of single-crystalline bismuth nanostructures as function of transverse dimensions: thickness *t* and width *w*. The *R*(*t,w*) data showed clear non-monotonous dependence with pronounced increase of resistance approaching the threshold diameter ~ 50 nm. Both effects are more pronounced at low temperatures. Our findings are consistent with model calculations taking into consideration size-dependent quantization of the energy spectrum of charge carriers in bismuth. The observed quantum size effect is of universal origin and should be present in any metallic conductor of relevant dimensions, and, hence, should be taken into consideration to optimize operation of the next generation of ultra-small quantum nanoelectronic circuits.

## Methods

Experimental confirmation of a size phenomenon typically involves studies of multiple structures with various values of the relevant dimension: e.g. thickness of a film or diameter of a wire. The approach assumes that all other parameters of the system, except the dimension(s), are the same. Unfortunately in practice the requirement it is not easy to satisfy, particularly if one deals with various artificially fabricated state-of-the-art small nanostructures. An alternative approach is to repeatedly re-measure the same sample after successive steps of size reduction.^{18, 19} Low energy directional ion milling^{20, 21} was used to progressively and non-destructively reduce the diameter of bismuth nanostructures. For acceleration energy ~1 keV the penetration depth Ar^{+} ions into bismuth matrix is about 1.5 nm, which is comparable to thickness of just few top-most atomic layers. Hence, the method can be considered as virtually non-destructive.^{18} Extra information on ion milling is provided in Supplementary information file*.* AFM and SEM analyses were made *ex situ* at room temperatures and atmosphere after each sputtering and resistance measurement sequence. As bismuth is a noble material, contamination of the surface due oxidation is not an issue.

Quantitative analysis of quantum size effects in a material with highly anisotropic Fermi surface as bismuth (Fig. 1a) requires fabrication of single-crystalline nanostructure with well-defined orientation of crystallographic axes. Several approaches have been taken to fabricate single crystalline micro- or nanowires of various materials.^{6, 7, 22, 16} Here we used conventional e-beam lift-off lithography and high speed ~1 nm/s 99.9995% pure bismuth e-gun deposition in vacuum ~10^{−6} mBar on mica substrate heated to 140 °C. The recipe provides high-quality polycrystalline films with large grain size up to 1.5 μm. X-ray analyses of co-deposited 2D films revealed that trigonal axis C3 is normal to the film plane (Fig. 1c), which is in agreement with literature data.^{3, 4, 14} Combined with ion beam milling^{20, 21} the method enables sequential electron transport measurements of bismuth nanostructures with progressively reduced transverse dimension ~(*t·w*)^{1/2} down to sub-50 nm scales. All structures were designed to enable electric resistance measurements using 4-probe configuration to reduce the impact of the probes (Fig. 1c). The absence of grain boundaries containing foreign material(s) (e.g. impurities or oxide) is confirmed by the linearity of V–I dependencies, typical for Ohmic conductors.

## References

- 1.
Tringides, M. C., Jatochowski, M. & Bauer, E. Quantum size effects in metallic nanostructures.

*Phys. Today***60**, 50–54 (2007). - 2.
Lutski, V. N. Quantum size effects—present state and perspectives of experimental investigations.

*Phys. Stat. Sol.***1**, 199–220 (1970). - 3.
Ogrin, Y. V., Lutskii, V. N. & Elinson, M. I. Observation of quantum size effects in thin bismuth films.

*Zh. Exper. Teor. Fiz. Pisma.***3**, 114–118 (1966). - 4.
Ogrin, Y. F.

*et al.*The temperature dependence of the specific resistance and the Hall constant of dimension-quantized bismuth films.*JETP***53**, 1218–1224 (1967). - 5.
Sun, X., Zhang, Z. & Dresselhaus, M. S. Theoretical modeling of thermoelectricity in Bi nanowires.

*Appl. Phys. Lett.***74**, 4005–4007 (1999). - 6.
Heremans, J.

*et al.*Bismuth nanowire arrays: synthesis and galvanomagnetic properties.*Phys. Rev. B***61**, 2921–2930 (2000). - 7.
Nikolaeva, A., Huber, T., Konopko, L. & Tsurkan, A. Observation of the semiconductor-semimetal and semimetal-semiconductor transitions in Bi quantum wires induced by anisotropic deformation and magnetic field.

*J. Low Temp. Phys.***158**, 530–535 (2010). - 8.
Costa-Krämer, J. L., Garcia, N. & Olin, H. Conductance quantization in bismuth nanowires at 4 K.

*Phys. Rev. Lett.***78**, 4990–4993 (1997). - 9.
Rodrigo, J. G., García-Martín, A., Sáenz, J. J. & Vieira, S. Quantum conductance in semimetallic bismuth nanocontacts.

*Phys. Rev. Lett.***88**, 246801 (2002). - 10.
McClure, J. W. & Choi, K. H. Energy band model and properties of electrons in bismuth.

*Solid State Commun.***21**, 1015–1018 (1977). - 11.
Lax, B. and Mavroides, J. G.,

*Solid State Physics*(Academic, 1960). - 12.
Zhang, Z.

*et al.*Electronic transport properties of single-crystal bismuth nanowire arrays.*Phys. Rev. B***61**, 4850–4861 (2000). - 13.
Sandomirskii, V. B. Quantum size effect in a semimetal film.

*Sov. Phys. JETP***25**, 101–106 (1967). - 14.
Cottey, A. A. Quantum size effect in bismuth films normal to trigonal axis.

*Phys. Stat. Sol.***72**, K63–K66 (1975). - 15.
Farhangfar, S. Quantum size effects in a one-dimensional semimetal.

*Phys. Rev. B***74**, 205318 (2006). - 16.
Farhangfar, S. Quantum size effects in solitary wires of bismuth.

*Phys. Rev. B***76**, 205437 (2007). - 17.
Kim, J., Lee, S., Brovman, Y. M., Kim, P. & Lee, W. Diameter-dependent thermoelectric figure of merit in single-crystalline Bi nanowires.

*Nanoscale***7**, 5053–5059 (2015). - 18.
Lehtinen, J. S.

*et al*Evidence of quantum phase slip effect in titanium nanowires.*Phys. Rev. B***85**, 094508 (2012). - 19.
Arutyunov, K. Yu

*et al.*Quantum phase slip phenomenon in ultra-narrow superconducting nanorings.*Sci. Rep.***2**, 293 (2012). - 20.
Savolainen, M.

*et al.*Ion beam sputtering method for progressive reduction of nanostructures dimensions.*Appl. Phys. A***79**, 1769–1773 (2004). - 21.
Zgirski, M.

*et al.*Ion beam shaping and downsizing of nanostructures.*Nanotechnology***19**, 055301 (2008). - 22.
Arutyunov, K. Yu, Ryynanen, T. V., Pekola, J. P. & Pavolotski, A. B. Superconducting transition of single-crystal tin microstructures.

*Phys. Rev. B***63**, 092506 (2001).

## Acknowledgements

The article was prepared within the framework of the Academic Fund Program at the National Research University Higher School of Economics (HSE) in 2016 (grant No. 16-05-0029) and supported within the framework of a subsidy granted to the HSE by the Government of the Russian Federation for the implementation of the Global Competitiveness Program

## Author information

### Affiliations

### Contributions

E.A.S. (MSc student) made calculations and theory fits. K.-P.R. (MSc student) fabricated nanostructures, made microscopic analyses and performed transport measurements. K.Y.A. (Prof.) initiated the project, designed variable temperature sample holder in the vacuum chamber, built measuring set-up, participated in transport measurements, suggested interpretation of the results and wrote the text.

### Corresponding author

## Ethics declarations

### Competing interests

These authors do not have any competing financial interests.

## Electronic supplementary material

## Rights and permissions

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

## About this article

### Cite this article

Sedov, E.A., Riikonen, KP. & Arutyunov, K.Y. Quantum size phenomena in single-crystalline bismuth nanostructures.
*npj Quant Mater* **2, **18 (2017). https://doi.org/10.1038/s41535-017-0017-8

Received:

Revised:

Accepted:

Published: