**Open Journal of Biophysics**

Vol.09 No.04(2019), Article ID:94625,15 pages

10.4236/ojbiphy.2019.94017

ISFET Based DNA Sensor: Current-Voltage Characteristic and Sensitivity to DNA Molecules

Lusine Gasparyan^{1,2}, Ilya Mazo^{1}, Vahan Simonyan^{3}, Ferdinand Gasparyan^{1,2 }

^{1}Argentis LLC, Gaithersburg, MD, USA

^{2}Yerevan State University, Yerevan, Armenia

^{3}The Department of Biochemistry & Molecular Medicine, The George Washington University, Washington DC, USA

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: August 7, 2019; Accepted: August 24, 2019; Published: August 27, 2019

ABSTRACT

Dependency of both source-drain current and current sensitivity of nanosize ISFET biosensor vs. concentration of DNA molecules in aqueous solution theoretically is investigated. In calculations it is carried out effects concerning charge carriers distribution in current channel and concerning carriers’ mobility behavior in high electrical fields in the channel. The influence of DNA molecules on the work of ISFET biosensors is manifested by a change in the magnitude of the gate surface charge. Starting with fairly low concentrations of DNA, ISFET sensors respond to the presence of DNA molecules in an aqueous solution which is manifested by modulation of channel conductance and therefore the source-drain current changes of the field-effect transistor. It is shown that the current sensitivity with respect to concentration of DNA molecules linearly depends on the source-drain voltage and reaches high values.

**Keywords:**

DNA, ISFET, Sensitivity, Nanosize, Current Channel

1. Introduction

The ion-sensitive field-effect transistor (ISFET) is one of the most popular semiconductor biosensors, and has been introduced as the first nanosized bio-chemical sensor. Currently, the use of ISFET technology encompasses a wide range of applications in a variety of areas, and those in the bioelectronic monitoring areas are particularly noteworthy. The ISFET sensor has been used to measure H^{+} or OH^{−} ions concentrations in aqueous solution, causing an interface potential on the gate insulator (oxide). Much attention has been paid to silicon based biosensors in the field of bio-analytical applications due to their favorable characteristics (Si-based technology, sensitivity, speed, signal-to-noise ratio, miniaturization, etc.). The introduction of the ISFET biosensor was in 1970 [1]; the first report regarding the use of an enzymatically modified ISFET for the direct detection of penicillin was in 1980 [2]. The operating principle of the ISFET device and recent advances and developments in the bio-analytical use of ISFET-based biosensors and detection strategies were covered in reviews [3] [4] [5]. Label-free detection of DNA using field-effect transistors (FET) with a real-time electrical readout system for rapid, cost-effective, and simple analysis of DNA samples has been proposed in [6]. When DNA strands bind to the free sites of gate surface of ISFETs, changes in surface electrical potential occur due to the negative charge of DNA molecule, shift in the source-drain current-voltage or capacitance-voltage characteristics thereby allowing for excellent performance of DNA sensing. Then conductivity of the semiconductor depletion layer will be modulated and source-drain signal current will be changed. It is clear that the capacitance of the system electrolyte-insulator-semiconductor depletion layer will be changed. In [7] DNA binding behavior was monitored using an ISFET biosensor, which was observed as changes in the threshold voltage. The change in DNA architecture and content, either due to hybridization or enzymatic reaction, yields a local pH charge variation and a rearrangement of ionic species near the sensor surface that modulate the sensor’s response. One of the main disadvantages in measuring the biomolecular recognition using FET-biosensors is the Debye screening length. For this reason, it is necessary for FET measurements that biological sensing should take place within the Debye length. Field-effect transistors are described as three-electrode devices in which the current flowing between the source and drain electrodes can be modulated by changing the potential applied to the gate and source electrodes [8]. The current-control mechanism is based on an electric field generated by the voltage applied to the gate layer [8]. The current is conducted by one type of carriers (electrons or holes) depending on the semiconductor type. In the case of a p-type semiconductor, when a positive gate voltage is applied, holes are repelled from the semiconductor-insulator interface creating a depletion layer. Applying a positive gate voltage attracts electrons to the semiconductor surface, when a sufficiently high concentration of electrons is accumulated in this region a conductive channel is created at the semiconductor-insulator interface allowing a current flow between source and drain. The gate voltage modulates the channel conductance. Some groups reported on the use of FETs for real-time amplification and detection of nucleic acid using pH sensing [9] [10] [11] [12]. Detail review of advances and developments in the bio-analytical use of ISFET-based biosensors presented in [3] [5]. An electrochemical sensor for detection of unlabeled ssDNA using peptide nucleic acid probes coupled to the FET gate is demonstrated in [13]. An application of ISFET technology for the detection of single nucleotide polymorphisms suggested in [14]. In this study authors developed a useful procedure for sequencing one base via the detection of single-base mismatch in DNA. Note that the sequencing mechanisms of sensors based on ISFETs can actually be used only for detecting nucleic acid using pH sensitivity and amplifying the useful signal in real time. Static, dynamic characteristics and pH sensitivity of bio FET sensors made on nanosize silicon (nanowire, nanoribbon) are detailed study by us in [15] [16] [17]. In [17] pH sensitivity of the biochemical sensors was introduced as
$\Delta {I}_{ds}/\Delta \text{pH}$ , where,
$\Delta {I}_{ds}$ and
$\Delta \text{pH}$ are the elementary changes in source-drain current and pH. Current-voltage characteristics, low-frequency noises, pH-sensitivity and SNR (signal-to-noise ratio) for liquid and back double gated Si NW (nanowire) sensors are investigated. It is shown that source-drain current substantially depends on pH-value and grows with increasing of the liquid gate voltage. The pH sensitivity increases as a result of the liquid gate voltage increase, giving opportunity to measure very low proton concentrations in the electrolyte medium at the high values of the front gate voltage. It is shown that SNR for Si NW based biochemical sensor has higher value, reaching up to 10^{5}. In [15] [16] it is shown that in Si nanosize FET biosensors pH sensitivity increases with the increase of current channel length approaching the Nernst limit value of 59.5 mV/pH, indicating that larger area devices are more suitable for the pH sensing. The pH sensitivity increases also with the increasing of the back-gate voltage and approaches to 59.5 mV/pH.

The analysis conducted above shows a high level of experimental research in this area. Very good results have been achieved in the field of sensors fabrication and pH-sensitivity; effective methods have been proposed for lowering the level of low-frequency noise in order to increase the signal-to-noise ratio. However, it seems to us that there are still many open questions when studying the physical mechanisms occurring in a semiconductor and interface semiconductor-insulator that determine and influence on the basic parameters of the sensors.

The aim of this research is the theoretical investigation, simulation and demonstration of the dependency of source-drain current of the silicon nanowire based ISFET biosensor vs. concentration of DNA molecules in aqueous solution. DNA detection mechanism and source-drain current sensitivity depending on the influence of the DNA molecules which occur in an aqueous solution over the Debye screening length will be investigated and discussed.

2. Physical Processes in the ISFET Biosensor

To study the source-drain current of the ISFET biosensor and its sensitivity to the presence of negatively charged DNA molecules in an aqueous solution, it is necessary to consider the physical processes occurring, in particular, at the interface between the gate insulator and the electrolyte. The main physical processes taking place in the ISFET biosensor for DNA molecule detection are sketched in Figure 1. It is presented a schematic representation of the device structure (a), distribution of the gate potential ${V}_{g}$ over the layered structure of the biosensor (b), and energy diagram of the electrolyte-insulator-semiconductor structure. In the left side of Figure 1 also shows the coordinate system used.

Figure 1. Schematic illustration of the ISFET biosensor (a), distribution of potentials (b), and energy diagram of ISFET in equilibrium ( ${V}_{g}=0$ ) (c). Here RE means reference electrode; ${V}_{g}$ is the gate voltage, ${V}_{sd}$ is the source-drain potential, ${E}_{F}$ is the semiconductor Fermi level energy, ${E}_{C}$ is the conduction and ${E}_{V}$ is the valence band energy, ${E}_{\text{redox}}$ is the redox potential. The colored circles in figure (a) show the DNA nucleotides adenine, thymine, guanine and cytosine.

Particularly silicon based structures and silicon oxide as an insulator will be discussed. In the Figure 1 RE is the reference electrode, ${V}_{g}$ is the applied gate voltage, ${\varphi}_{Si}$ , ${\varphi}_{ch}$ , ${\varphi}_{ox}$ and ${\varphi}_{dl}$ are potentials of the silicon substrate, current channel (semiconductor depletion layer), oxide layer and double layer, correspondingly.

The balance equation for the potentials according to Figure 1, can be represented as follows:

${V}_{g}={\varphi}_{s}+{\varphi}_{Si}+{\varphi}_{ch}+{\varphi}_{ox}+{\varphi}_{dl}$ . (1)

To estimate these potentials as well as the threshold voltage, ${V}_{th}$ , and flat-band voltage, ${V}_{FB}$ , we can use the following relations [8] [18] - [23] :

$\begin{array}{l}{V}_{th}={V}_{FB}+2{\phi}_{F}+{\varphi}_{ch};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{V}_{FB}={\varphi}_{bulk,sol}-{\varphi}_{ch}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}};\\ {\varphi}_{bulk,sol}\approx 0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\varphi}_{Si}={\varphi}_{bulk,Si}\approx 0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\phi}_{F}=2{\phi}_{T}\mathrm{ln}\frac{{N}_{A}}{{n}_{i}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\varphi}_{ch}=\sqrt{\frac{4q{\epsilon}_{0}{\epsilon}_{Si}{N}_{A}{\phi}_{T}}{{C}_{ox}^{2}}};\\ {\phi}_{T}=\frac{{k}_{B}T}{q};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\varphi}_{dl}=2{\phi}_{T}\left(\frac{{\epsilon}_{w}}{{\epsilon}_{r}}\frac{{N}_{sol}}{{K}_{AK}^{+}+{H}_{s}^{+}}\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\varphi}_{ox}=\frac{q{N}_{t}}{{C}_{ox}}.\end{array}$ (2)

Here q is the elementary charge;
${k}_{B}$ is the Boltzmann’s constant; T is the absolute temperature;
${\phi}_{T}$ is the thermal voltage;
${\phi}_{F}$ is the Fermi potential;
${\varphi}_{bulk,sol}$ and
${\varphi}_{bulk,Si}$ are the electric potentials of the bulk solution and the bulk silicon substrate;
${\varphi}_{dl}$ is the potential of double layer;
${\Phi}_{Si}$ and
${\Phi}_{ox}$ are the work functions of silicon and silicon oxide (SiO_{2}), correspondingly;
${Q}_{ox}$ is the oxide layer charge per unit area,
${C}_{ox}$ is the capacitance of the oxide layer per unit area;
${\epsilon}_{0}$ ,
${\epsilon}_{Si}$ ,
${\epsilon}_{ox}$ ,
${\epsilon}_{w}$ and
${\epsilon}_{r}$ are the dielectric permittivities of free space, silicon, silicon dioxide, water and electrolyte, respectively;
${N}_{A}$ is the doping acceptor concentration in p-Si substrate;
${n}_{i}$ is the intrinsic carrier concentration in bulk silicon;
${K}_{AK}^{+}$ is the molar concentration of the cations in the solution,
${H}_{s}^{+}$ is the molar concentration of the hydrogen ions at the oxide surface;
${N}_{sol}$ is the molar concentration of the solution;
${N}_{t}$ is the concentration of surface open electronic binding sites (traps) per oxide unit area. Note that the redox potential
${E}_{\text{redox}}$ is a measure of the ease with which a molecule will accept electrons and double layer in solution consist of IHL (Inner Helmholtz layer), OHL (Outer Helmholtz layer) and GCL (Gouy-Chapman layer) [24].

The main physical processes occur in the conductive channel. Therefore, for further calculations, it is necessary to determine the surface potential of the interface between the semiconductor layer (channel) and insulator ${\varphi}_{ch}$ . This can be calculated using Equation (1) and Equation (2) and expressions for the density of minority carriers in semiconductor. For ${\varphi}_{ch}$ we receive:

${\varphi}_{ch}={\phi}_{T}\mathrm{ln}\left(\frac{\eta {C}_{ox}{\phi}_{T}{N}_{A}}{qt{n}_{i}^{2}}\right)+{\phi}_{T}\mathrm{ln}\left\{\mathrm{ln}\left[1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}-{V}_{th}}{\eta {\phi}_{T}}\right)\right]\right\}$ . (3)

where

$\eta =1+\frac{{C}_{d}}{{C}_{ox}}\approx 1+\sqrt{\frac{q{\epsilon}_{0}{\epsilon}_{Si}{N}_{A}}{2{\phi}_{T}{C}_{ox}^{2}}}$

is the factor of the transistor non-ideality ( ${C}_{d}$ is the capacitance of the silicon depletion layer per unit area).

We consider the case of an inversion n-channel liquid-gated FET (Figure 1(c)). It is clear that the majority of processes in the structure are therefore determined by the electrons. The channel source-drain current consists of drift and diffusion components. It is well known that the diffusion component is dominant in the sub-threshold mode and the drift component is dominant in the over-threshold region. The channel source-drain current in Y direction ${I}_{sd}\left(y\right)$ can be calculated using the following equation [8] :

${I}_{sd}\left(y\right)={\mu}_{ef}w\left[{Q}_{ch}\frac{\text{d}{V}_{sd}}{\text{d}y}+{\phi}_{T}\frac{\text{d}{Q}_{ch}\left(y\right)}{\text{d}y}\right]$ . (4)

Here w is the channel width in Z direction, ${\mu}_{ef}$ is the effective mobility and ${Q}_{ch}$ is the charge density of the channel mobile carriers ( $\left[{Q}_{ch}\right]=\left[\text{C}/{\text{cm}}^{2}\right]$ ):

${Q}_{ch}={\displaystyle {\int}_{0}^{t}qn\left(x,{V}_{g}\right)\text{d}x}$ . (5)

Here t is the thickness of current channel in X direction, $n\left(x,{V}_{g}\right)$ is the electron’s concentration in the channel. The behavior of the source-drain current is defined by the distribution of the concentration of the mobile charge carriers over the conducting channel. Obviously, the concentration of mobile carriers in the channel depends on both the coordinate $x$ (see Figure 1(a)) and the applied gate voltage. At the same time, the charge surface concentration only depends on the gate voltage. Hence the overall concentration can be presented as follows:

$n\left(x,{V}_{g}\right)={n}_{s}\left({V}_{g}\right)\times f\left(x,{V}_{g}\right)$ . (6)

Here
${n}_{s}\left({V}_{g}\right)$ is the electron surface concentration per unit area at the oxide interface and
$f\left(x,{V}_{g}\right)$ in unit of [cm^{−1}] is the function which describes the charge carrier distribution in the X-Z plane of the channel (Figure 1) [25].

The surface concentration can be described using the unified charge control model from expression [26] :

${V}_{g}-{V}_{th}=\frac{q}{{C}_{ox}}\left({n}_{s}-{n}_{s,t}\right)+\eta {V}_{th}\mathrm{ln}\left(\frac{{n}_{s}}{{n}_{s,t}}\right)$ , (7)

where ${n}_{s,t}$ is the surface density of electrons per unit area at the threshold voltage: ${n}_{s,t}={n}_{s}$ at the ${V}_{g}={V}_{th}$ . It should be noted that the influence of the charge states of the electrolyte is determined by the value of ${V}_{th}$ (see Equation (2)). The concentration ${n}_{s,t}$ can be expressed as:

${n}_{s,t}=\frac{\eta {C}_{ox}{\phi}_{T}}{2q}$ . (8)

Equation (7) has no analytical solution for ${n}_{s}$ in terms of ${V}_{g}$ . The following approximate solution is suitable for strong inversion and sub-threshold regimes [22] :

${n}_{s}=2{n}_{s,t}\mathrm{ln}\left[1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}-{V}_{th}}{\eta {\phi}_{T}}\right)\right]$ . (9)

After determining ${n}_{s}\left({V}_{g}\right)$ , we must also calculate the function $f\left(x,{V}_{g}\right)$ in order to evaluate the influence of peculiarities of the carrier distribution on the physical processes taking place in the channel.

In order to find function $f\left(x,{V}_{g}\right)$ for the case of the quasi classical approach, we use the following dependence of $n\left(x\right)$ [8] :

$n\left(x\right)={N}_{c}\mathrm{exp}\left[-\frac{\left({E}_{c}-q\varphi \left(x\right)\right)-{E}_{F}}{{k}_{B}T}\right]={n}_{0}\mathrm{exp}\left[\frac{\varphi \left(x\right)}{{\phi}_{T}}\right]$ . (10)

Here ${N}_{c}$ is the density of states in the conduction band of a semiconductor, ${E}_{c}$ is the semiconductor conduction band energy, $\varphi \left(x\right)$ is the contact potential at the oxide-channel interface. To determine $\varphi \left(x\right)$ we have to solve the Poisson equation:

$\frac{{\text{d}}^{2}\varphi \left(x\right)}{\text{d}{x}^{2}}=-\frac{\rho \left(x\right)}{{\epsilon}_{0}{\epsilon}_{Si}}$ . (11)

Here $\rho \left(x\right)$ is the space charge density for the fully ionized acceptor centers in semiconductor (usually it is boron in silicon):

$\rho \left(x\right)=-q\left({N}_{A}^{-}+n-p\right)=-q{p}_{0}\left[1-\mathrm{exp}\left(-\frac{\varphi}{{\phi}_{T}}\right)+\frac{{n}_{0}}{{p}_{0}}\mathrm{exp}\left(\frac{\varphi}{{\phi}_{T}}\right)\right]$ . (12)

Here $n,p$ and ${n}_{0},{p}_{0}$ are the concentrations of the non-equilibrium and equilibrium electrons and holes, respectively, ${N}_{A}^{-}$ is the concentration of negatively charged acceptors. We can use following boundary conditions to solve Equation (11) (see Figure 1):

$x\to \infty \Rightarrow \varphi \to 0$ , $x\to 0\Rightarrow \varphi \to {\varphi}_{ch}$ . (13)

Using Equation (12) and boundary conditions (13), we obtain the following solution of Equation (11):

$\varphi \left(x\right)=\left\{{\varphi}_{ch}+\frac{q{n}_{0}}{{\epsilon}_{0}{\epsilon}_{Si}}\left[1-\mathrm{exp}\left(-\frac{x}{{l}_{s}}\right)\right]\right\}\mathrm{exp}\left(-\frac{x}{{l}_{s}}\right)$ , (14)

where

${l}_{s}=\frac{{L}_{D}}{1+{n}_{0}/{p}_{0}}$ , ${L}_{D}=\sqrt{\frac{{\epsilon}_{0}{\epsilon}_{Si}{\phi}_{T}}{q{p}_{0}}}$ , (15)

${L}_{D}$ is the Debye screening length.

Then using expression for ${\varphi}_{ch}$ from Equation (3) finally we have for the function $f\left(x,{V}_{g}\right)$ :

$\begin{array}{l}f\left(x,{V}_{g}\right)=\frac{{n}_{0}}{{n}_{s}}\times \mathrm{exp}\{[\mathrm{ln}\left(\frac{\eta {\phi}_{T}{C}_{ox}{N}_{A}}{qt{n}_{i}^{2}}\right)+\mathrm{ln}\left[\mathrm{ln}\left(1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}-{V}_{th}}{\eta {\phi}_{T}}\right)\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{q{t}^{2}{n}_{i}^{2}}{{\epsilon}_{0}{\epsilon}_{Si}{\phi}_{T}{N}_{A}}]\mathrm{exp}\left(-\frac{x}{{l}_{s}}\right)\}\end{array}$ . (16)

The field caused by the applied gate voltage in the inversion layer of liquid-gated FETs changes the transport behavior of the charge carriers and results in more frequent scattering events than in the absence of the gate voltage. The carrier’s mobility degrades as the result of scattering processes [27] [28] [29]. The mobility dependence on the transversal electric field (Y direction) at the applied gate voltage was taken into account using the following empiric equation [30] :

${\left({\mu}_{ef}\right)}_{x}={\mu}_{0}-\theta \left({V}_{g}+{V}_{th}\right)$ , (17)

where
${\mu}_{0}$ is the low-field magnitude of the mobility,
$\theta $ is the coefficient taken as 28 cm^{2}/(V^{2}s) [30] [31]. Since the modeling and the measurements are performed for low drain biases in linear mode, the effect of the electron velocity saturation on the drain current can be neglected. Other authors in [32] assume that when Coulomb scattering dominates the mobility dependence on the gate voltage
$\mu \left({V}_{g}\right)$ in the vicinity
${V}_{g}={V}_{th}$ can be linearized in the form [33] [34]

$\mu \left({V}_{g}\right)={\mu}_{0}+\theta \left({V}_{g}-{V}_{th}\right)$ , (18)

where ${\mu}_{0}={\mu |}_{{V}_{g}={V}_{th}}$ , $\theta ={\frac{\text{d}\mu \left({V}_{g}\right)}{\text{d}{V}_{g}}|}_{{V}_{g}={V}_{th}}$ , and in general case $\theta $ can be positive or negative.

In further calculations we will use expression (17) for the mobility of major carriers in the channel.

3. Source-Drain Current

Using Equation (4) and Equation (17) we can present the drift component of the source-drain current by following expression:

${I}_{sd}\left(y\right)\approx \frac{qw{n}_{s}{V}_{ds}}{l}\left[{\mu}_{0}-\theta \left({V}_{G}+{V}_{th}\right)\right]{\displaystyle {\int}_{0}^{t}f\left(x,{V}_{g}\right)\text{d}x}$ ,

where ${n}_{s}$ is determined from Equation (9).

For the integral

${\int}_{0}^{t}f\left(x,{V}_{g}\right)\text{d}x$

we have

${\int}_{0}^{t}f\left(x,{V}_{g}\right)\text{d}x}=\frac{{n}_{0}}{{n}_{s}}{\displaystyle {\int}_{0}^{t}\mathrm{exp}\left[G\mathrm{exp}\left(-\frac{x}{{l}_{s}}\right)\right]\text{d}x}\approx \frac{{n}_{0}}{{n}_{s}}t\left[1+\frac{{l}_{s}}{t}G\left(1-{\text{e}}^{-t/{l}_{s}}\right)\right]$ . (19)

Here

$G\equiv \frac{1}{B}+\mathrm{ln}B+\mathrm{ln}\left[\mathrm{ln}\left(1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}-{V}_{th}}{{\phi}_{T}}\right)\right)\right]$ , $B\equiv \frac{{\phi}_{T}{\epsilon}_{0}{\epsilon}_{ox}{N}_{A}}{q{t}^{2}{n}_{i}^{2}}$ . (20)

Thus

${I}_{sd}\left(y\right)\approx qw{n}_{0}{V}_{ds}\frac{t}{l}\left[{\mu}_{0}-\theta \left({V}_{G}+{V}_{th}\right)\right]\left[1+\frac{{l}_{s}}{t}G\left(1-{\text{e}}^{-t/{l}_{s}}\right)\right]$ . (21)

Behavior of dependency ${I}_{sd}$ from number of DNA molecules in solution ${N}_{\text{DNA}}$ can be determined from dependency ${V}_{th}\left({Q}_{ox}\right)$ .

For simplicity of further calculations assume that $\eta \approx 1$ and taking account that oxide layer capacitance for unit area

${C}_{ox}=\frac{{\epsilon}_{0}{\epsilon}_{ox}}{t}$ . (22)

For ${V}_{th}$ from (2) we have:

${V}_{th}=2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}$ . (23)

Substituting (23) into (21) gives:

$\begin{array}{l}{I}_{sd}\left({Q}_{ox}\right)\approx qw{n}_{0}{V}_{ds}\frac{t}{l}\left[{\mu}_{0}-\theta \left({V}_{G}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}\right)\right]\\ \times \left\{1+\frac{{l}_{s}}{t}\left\{\frac{1}{B}+\mathrm{ln}B+\mathrm{ln}\left[\mathrm{ln}\left(1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}}{{\phi}_{T}}\right)\right)\right]\right\}\left(1-{\text{e}}^{-t/{l}_{s}}\right)\right\}.\end{array}$ (24)

The influence of the oxide interface on the source-drain current of the transistor can be taken into account by the change of the charge of the oxide interface traps
${Q}_{ox}=q{N}_{t}$ (
${N}_{t}$ is the surface concentration of oxide interface all proton donors and proton acceptors traps in units cm^{−2}, Figure 2). Depending on presence of negatively charged DNA molecules the charge on the oxide surface will decrease (compensated) and become

${{Q}^{\prime}}_{ox}=q{N}_{t}^{+}\left(1-\delta \right)$ , $\delta \equiv \frac{{N}_{\text{DNA}}}{{N}_{t}^{+}}$ , (25)

where ${N}_{t}^{+}$ and ${N}_{\text{DNA}}$ are the surface concentrations of positively charged proton acceptor ${\text{OH}}_{2}^{+}$ traps (Figure 2) and DNA surface concentration in solution

Figure 2. Site binding on the interface SiO_{2}-solution and the process of binding of the negatively charged DNA on the
${\text{OH}}_{2}^{+}$ free bond (a), SiO_{2} tetrahedron structure (b).

near the oxide at a distance of the Debye length. It is clear that (25) is correct for the ${N}_{\text{DNA}}\le {N}_{t}^{+}$ , or $\delta \le 1$ . In the case of super compensation when ${N}_{\text{DNA}}>{N}_{t}^{+}$ the DNA additional molecules do not bind on the surface positively charged sites (proton acceptors, Figure 2) and ISFET biosenor do not sense this DNA additional molecules.

Figure 2 show site binding on the interface SiO_{2}-solution, process of negatively charged DNA molecule binding and silicon dioxide tetrahedron structure.

Consider the source-drain current sensitivity of the ISFET biosensor to DNA molecules S as a change in source-drain current $\Delta {I}_{sd}$ for a corresponding change in the proportion of DNA molecules in the solution $\Delta \delta $ :

$S=\left|\frac{\Delta {I}_{sd}}{\Delta \delta}\right|,\text{A}$ (26)

As $\Delta \delta $ is dimensionless parameter sensitivity will be measured by the Ampere.

4. Numerical Simulation and Discussion

For numerical computation, we use the following values, which correspond to the sample geometry and the parameters of the materials for the investigated nanosize structure at the room temperature: ${\mu}_{0}=260\text{\hspace{0.17em}}{\text{cm}}^{\text{2}}/\left(\text{V}\cdot \text{s}\right)$ [35] , $\theta =28\text{\hspace{0.17em}}{\text{cm}}^{\text{2}}/\left({\text{V}}^{2}\cdot \text{s}\right)$ [30] [31] , ${\phi}_{T}=0.026\text{\hspace{0.17em}}\text{V}$ , ${N}_{sol}=0.015\text{\hspace{0.17em}}\text{mol}/\text{l}$ [25] , ${K}_{AK}^{+}=0.001\text{\hspace{0.17em}}\text{mol}/\text{l}$ [25] , $w=150\text{\hspace{0.17em}}\text{nm}$ , $l=200\text{\hspace{0.17em}}\text{nm}$ , $t=20\text{\hspace{0.17em}}\text{nm}$ , ${\epsilon}_{Si}=11.6$ , ${\epsilon}_{ox}=3.9$ , ${\epsilon}_{w}\approx 80$ , ${\epsilon}_{r}\approx 78$ , ${\epsilon}_{0}=8.85\times {10}^{-14}\text{F}/\text{cm}$ , $q=1.6\times {10}^{-19}\text{C}$ , ${\Phi}_{Si}=4.85\text{\hspace{0.17em}}\text{eV}$ [36] , ${\Phi}_{ox}=5\text{\hspace{0.17em}}\text{eV}$ [36] , ${N}_{A}={10}^{15}{\text{cm}}^{-3}$ , ${N}_{V}=2.5\times {10}^{19}{\text{cm}}^{-3}$ ( ${m}_{n}^{*}={m}_{p}^{*}\equiv {m}_{0}=9.1\times {10}^{-31}\text{kg}$ , ${m}_{n}^{*}$ and ${m}_{p}^{*}$ are effective mass of electrons and holes, ${m}_{0}$ is the free electron mass), ${p}_{0}\approx {N}_{A}={10}^{15}{\text{cm}}^{-3}$ . Let’s consider the case when “power of hydrogen” of the solution is equal to 7 ( $\text{pH}=7$ ). As

$\text{pH}=-\mathrm{log}\left[{\text{H}}^{+}\right]$ ,

we get

$\left[{\text{H}}^{+}\right]={10}^{-\text{pH}}={10}^{-7}\text{mol}/\text{l}$ .

For the electron concentration in the inversion layer we can assume that it is equal to majority carrier’s concentration in p-Si, e.g.
${n}_{0}\left(\text{inv}\right)\approx {p}_{0}={10}^{15}{\text{cm}}^{-3}$ . For the
${N}_{t}^{+}$ we can do following estimation. It is assumed that traps concentration on the interface silicon oxide-electrolyte is the same as silicon oxide-Si. According to data [36] traps concentration in Si-SiO_{2} interface is about
$\left({10}^{10}\text{-}{10}^{11}\right){\text{cm}}^{-2}$ . In further calculation we will use
${N}_{t}^{+}\approx {10}^{11}{\text{cm}}^{-2}$ .

For numerical simulation let’s simplify expression for source-drain current assuming that:

1) At the noted above parameters $B\propto {10}^{10}$ and we can ignore term $\frac{1}{B}$ compared $\mathrm{ln}B$ as $\frac{1}{B}\ll \mathrm{ln}B$ ;

2) $\begin{array}{l}\left[\mathrm{ln}\left(1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}}{{\phi}_{T}}\right)\right)\right]\\ \approx \frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}}{{\phi}_{T}}\right);\end{array}$

3) $\begin{array}{l}\mathrm{ln}\left[\mathrm{ln}\left(1+\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}}{{\phi}_{T}}\right)\right)\right]\\ \approx \mathrm{ln}\frac{1}{2}\mathrm{exp}\left(\frac{{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}}{{\phi}_{T}}\right);\end{array}$

4) $1-{\text{e}}^{-t/{l}_{s}}\approx 1-\frac{t}{{l}_{s}}$ .

Thus source-drain current can be presented as follows:

$\begin{array}{l}{I}_{sd}\approx \frac{qtw{n}_{0}{V}_{ds}}{l{\phi}_{T}}\left[{\mu}_{0}-\theta \left({V}_{G}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\times \left[{\phi}_{T}\left(1+\mathrm{ln}\frac{B}{2}\right)+{V}_{g}+2{\phi}_{F}+{\varphi}_{dl}-\frac{{\Phi}_{Si}-{\Phi}_{ox}}{q}+\frac{{Q}_{ox}}{{C}_{ox}}\right].\end{array}$ (27)

Results of numerical calculations of source-drain current ${I}_{sd}$ vs source-drain voltage ${V}_{ds}$ and DNA concentration are presented in Figures 3-5. The error in plotting the dependencies in Figures 3-5 does not exceed (5 - 7)%. In order not to complicate the graphs, these errors are not shown in the figures. Note that fluctuations of values by (5 - 7)% particularly do not affect the course of dependencies and do not change the mechanisms for explaining their behavior. As expected, the dependence of ${I}_{sd}$ on the source-drain voltage is linear. Dependency ${I}_{sd}$ vs DNA concentration ${N}_{\text{DNA}}$ (or $\delta $ ) is very weak (see Figure 4). At the beginning ( $0<\delta <0.1$ ) source-drain current increases very weakly. This

Figure 3. Dependency of the source-drain current vs source-drain voltage.

Figure 4. Dependency of the source-drain current vs DNA molecules concentration in solution.

Figure 5. Dependency of the current sensitivity vs. source-drain voltage at the several values of DNA molecules concentration ${N}_{\text{DNA}}$ or $\delta $ .

behavior can be explained as follows. Assume that all proton donors (OH^{−}) sites in the interface oxide-electrolyte (Figure 2) fully accept protons (H^{+}) from solution and changes in their charges are not significant and can be neglected during the sensor operation. At low concentration of DNA molecules
$0<\delta <0.1$ , they are bind very weakly or not bind with free proton acceptor sites
${\text{OH}}_{2}^{+}$ due to its high diffusion activity and the presence of a Coulomb barrier near proton acceptor sites. With increasing DNA concentration and therefore decreasing their diffusion activity (when
$0.1<\delta <0.8$ ) negatively charged DNA molecules overcome Coulomb barrier near the proton acceptor sites and bond with them on the oxide surface. As a result the positive surface charge decreases and correspondingly decreases deepness of current channel and its conductivity (source-drain current). At the super compensation of the proton acceptor sites
${\text{OH}}_{2}^{+}$ (
$\delta \to \text{1}$ , high concentration of DNA molecules) in oxide-electrolyte interface
${I}_{sd}\left(\delta \right)$ dependency has increasing behavior. Probably it is conditioned by the ionic Coulomb blockade effect [37]. As it is known this effect appears in the electro-diffusive transport of charged ions (in our case DNA molecules) through sub-nanometer artificial nanopores [38] or biological ion channels [39]. DNA molecules cannot overcome Coulomb barrier of the proton acceptor sites and channel conductance again increase. Consequently, the current increases.

The dependence of current sensitivity on source-drain voltage for several values of DNA molecules concentration is presented in Figure 5. These dependencies are constructed according to formula (26) using the data in Figure 4. They are also linear dependency. Depending on values of DNA concentration in solution and source-drain applied voltage current sensitivity can reach comparatively high values (up to 5 … 9 μA/δ in the case of parameters chosen above).

Based on the above reasoning, we can draw the following conclusions.

• ISFET nanosized structures can be used for detecting charged DNA molecules.

• The influence of DNA molecules on the work of ISFET biosensors is manifested by a change in the magnitude of the surface charge of the gate electrode.

• Starting with fairly low concentrations of DNA, ISFET biosensors respond to the presence of DNA molecules in an aqueous solution which is manifested by modulation of channel conductance or the source-drain current.

• Current sensitivity linearly depends on the source-drain voltage and reaches high values.

Thus, ISFET nanosized silicon biosensors can be successfully used to detect very low concentrations of DNA molecules in an aqueous solution with high sensitivity. As an advantage note also the compatibility of silicon based devices with modern CMOS technology.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Gasparyan, L., Mazo, I., Simonyan, V. and Gasparyan, F. (2019) ISFET Based DNA Sensor: Current-Voltage Characteristic and Sensitivity to DNA Molecules. Open Journal of Biophysics, 9, 239-253. https://doi.org/10.4236/ojbiphy.2019.94017

References

- 1. Bergveld, P. (1970) Development of an Ion Sensitive Solid State Device for Neurophysiological Measurements. IEEE Transactions on Biomedical Engineering, BME-17, 70-71. https://doi.org/10.1109/TBME.1970.4502688
- 2. Caras, S. and Janata, J. (1980) Field Effect Transistor Sensitive to Penicillin. Analytical Chemistry, 52, 1935-1937. https://doi.org/10.1021/ac50062a035
- 3. Lee, C.-S., Kim, S.K. and Kim, M. (2009) Ion-Sensitive Field-Effect Transistor for Biological Sensing. Sensors, 9, 7111-7131. https://doi.org/10.3390/s90907111
- 4. Veigas, B., Fortunato, E. and Baptista, P.V. (2015) Field Effect Sensors for Nucleic Acid Detection: Recent Advances and Future Perspectives. Sensors, 15, 10380-10398. https://doi.org/10.3390/s150510380
- 5. Hu, Y. (2015) Advanced Sensing and Processing Methodologies for ISFET Based DNA Sequencing. Imperial College London, London.
- 6. Ingebrandt, S., Han, Y., Nakamura, A., Poghossian, A., Schöning, M.J. and Offenhausser, A. (2007) Label-Free Detection of Single Nucleotide Polymorphisms Utilizing the Differential Transfer Function of Field-Effect Transistors. Biosensors and Bioelectronics, 22, 2834-2840. https://doi.org/10.1016/j.bios.2006.11.019
- 7. Goncalves, D., Prazeres, D., Chu, V. and Conde, J. (2008) Detection of DNA and Proteins Using Amorphous Silicon Ion-Sensitive Thin-Film Field Effect Transistors. Biosensors and Bioelectronics, 24, 545-551.
- 8. Sze, S.M. (1986) Physics of Semiconductor Devices. Wiley, New York.
- 9. Veigas, B., Branquinho, R., Pinto, J.V., Wojcik, P.J., Martins, R., Fortunato, E. and Baptista, P.V. (2014) Ion Sensing (EIS) Real-Time Quantitative Monitorization of Isothermal DNA Amplification. Biosensors and Bioelectronics, 52, 50-55. https://doi.org/10.1016/j.bios.2013.08.029
- 10. Wong, M.L. and Medrano, J.F. (2005) Real-Time PCR for mRNA Quantitation. BioTechniques, 39, 75-85. https://doi.org/10.2144/05391RV01
- 11. Branquinho, R., Veigas, B., Pinto, J.V., Martins, R., Fortunato, E. and Baptista, P.V. (2011) Real-Time Monitoring of PCR Amplification of Proto-Oncogene c-MYC Using a Ta2O5 Electrolyte-Insulator-Semiconductor Sensor. Biosensors and Bioelectronics, 28, 44-49. https://doi.org/10.1016/j.bios.2011.06.039
- 12. Toumazou, C., Shepherd, L., Reed, S., Chen, G., et al. (2013) Simultaneous DNA Amplification and Detection Using a pH-Sensing Semiconductor System. Nature Methods, 10, 641-646.
- 13. Kaisti, M., Kerko, A., Aarikka, E., Saviranta, P., Boeva, Z., Soukka, T. and Lehmusvuori, A. (2017) Real-Time Wash-Free Detection of Unlabeled PNA-DNA Hybridization. Scientific Reports, 7, Article No. 15734. https://doi.org/10.1038/s41598-017-16028-7
- 14. Purushothaman, S., Toumazou, C. and Ou, C. (2006) Protons and Single Nucleotide Polymorphism Detection: A Simple Use for the Ion Sensitive Field Effect Transistor. Sensors and Actuators B: Chemical, 114, 964-968. https://doi.org/10.1016/j.snb.2005.06.069
- 15. Gasparyan, F., Zadorozhnyi, I., Khondkaryan, H., Arakelyan, A. and Vitusevich, S. (2018) Photoconductivity, pH Sensitivity, Noise, and Channel Length Effects in Si Nanowire FET Sensors. Nanoscale Research Letters, 13, 87-96. https://doi.org/10.1186/s11671-018-2494-5
- 16. Gasparyan, F., Zadorozhnyi, I., Khondkaryan, H., Arakelyan, A. and Vitusevich, S. (2017) Biochemical Sensors Based on Silicon Nanoribbon FETs. Part 1: Samples Fabrication, CVCs, pH-Sensitivity. Proceedings of 11th International Conference on Semiconductor Micro- & Nanoelectronics, Yerevan, Armenia, 23-25 June 2017.
- 17. Gasparyan, F.V. and Khondkaryan, H.D. (2015) Low-Frequency Noises and SNR of Double Gated Si NW ISFET Based Biochemical Sensor. Proceedings of 10th International Conference on Semiconductor Micro- & Nanoelectronics, Yerevan, Armenia, 11-13 September 2015, 68-71.
- 18. Gasparyan, F.V., Poghossian, A., Vitusevich, S.A., Petrychuk, M.V., Sydoruk, V.A., Siqueira, J.R., Oliveira, O.N., Offenhäusser, A. and Schöning, M.J. (2011) Low-Frequency Noise in Field-Effect Devices Functionalized with Dendrimer/Carbon-Nanotube Multilayers. IEEE Sensors Journal, 11, 142-149.https://doi.org/10.1109/JSEN.2010.2052355
- 19. Deen, M.J., Shinwari, M.W., Ranuárez, J.C. and Landheer, D. (2006) Noise Considerations in Field-Effect Biosensors. Journal of Applied Physics, 100, Article ID: 074703. https://doi.org/10.1063/1.2355542
- 20. Hassibi, A., Navid, R., Dutton, R.W. and Lee, T.H. (2004) Comprehensive Study of Noise Processes in Electrode Electrolyte Interfaces. Journal of Applied Physics, 96, 1074-1082. https://doi.org/10.1063/1.1755429
- 21. Janicki, M., Daniel, M., Szermer, M. and Napieralski, A. (2004) Ion Sensitive Field Effect Transistor Modelling for Multidomain Simulation Purposes. Microelectronics Journal, 35, 831-840. https://doi.org/10.1016/j.mejo.2004.06.015
- 22. Ytterdal, T., Cheng, Y. and Fjeldly, T.A. (2003) Device Modeling for Analog and RF CMOS Circuit Design. John Wiley & Sons, New York. https://doi.org/10.1002/0470863803
- 23. Tarasov, A., Wipf, M., Stoop, R.L., Bedner, K., Fu, W., Guzenko, V.A., Knopfmacher, O., Calame, M. and Scho, C. (2012) Understanding the Electrolyte Background for Biochemical Sensing with Ion-Sensitive Field-Effect Transistors. ACS Nano, 6, 9291-9298. https://doi.org/10.1021/nn303795r
- 24. Nakamura, M., Sato, N., Hoshi, N. and Sakata, O. (2011) Outer Helmholtz Plane of the Electrical Double Layer Formed at the Solid Electrode-Liquid Interface. ChemPhysChem, 12, 1430-1434. https://doi.org/10.1002/cphc.201100011
- 25. Pud, S., Gasparyan, F., Petrychuk, M., Li, J., Offenhausser, A. and Vitusevich, S.A. (2014) Single Trap Dynamics in Electrolyte-Gated Si-Nanowire Field Effect Transistors. Journal of Applied Physics, 115, Article ID: 233705. https://doi.org/10.1063/1.4883757
- 26. Gao, X.P.A., Zheng, G. and Lieber, C.M. (2010) Subthreshold Regime has the Optimal Sensitivity for Nanowire FET Biosensors. Nano Letters, 10, 547-552. https://doi.org/10.1021/nl9034219
- 27. Sabnis, A.G. and Clemens, J.T. (1979) Characterization of the Electron Mobility in the Inverted (100) Si Surface. International Electron Devices Meeting, 25, 18-21. https://doi.org/10.1109/IEDM.1979.189528
- 28. Woo, J. and Cheng, B. (1996) Measurement and Modeling of the N-Channel and P-Channel MOSFET’s Inversion Layer Mobility at Room and Low Temperature Operation. Journal de Physique IV Colloque, 6, C3-43-C3-47.
- 29. Zhou, X. and Lim, K.Y. (2000) A Novel Approach to Compact I-V Modeling for Deep-Submicron MOSFET’s Technology Development with Process Correlation. Technical Proceedings of the International Conference on Modeling and Simulation of Microsystems (MSM 2000), San Diego, CA, 27-29 March 2000, 333-336.
- 30. Park, C., Lee, C., Lee, K., Moon, B.-J., Byun, Y.H. and Shur, M. (1991) A Unified Current-Voltage Model for Long-Channel nMOSFETs. IEEE Transactions on Electron Devices, 38, 399-406. https://doi.org/10.1109/16.69923https://www.academia.edu/23535536/A_unified_current-voltage_model_for_long-channel_nMOSFETs
- 31. Gildenblat, G.S. and Huang, C.-L. (1989) Engineering Model of Inversion Channel Mobility for 60-300k Temperature Range. Electronics Letters, 25, 634-636. https://doi.org/10.1049/el:19890430
- 32. Rudenko, T., Kylchitska, V., Khairuddin, M. Arshad, M., Raskin, J.-P., Nazarov, A. and Flandre, D. (2011) On the MOSFET Threshold Voltage Extraction by Transconductance and Transconductance-to-Current Ratio Change Methods: Part I—Effect of Gate-Voltage-Dependent Mobility. IEEE Transactions on Electron Devices, 58, 4172-4178. https://doi.org/10.1109/TED.2011.2168226
- 33. Takagi, S., Iwase, M. and Toriumi, A. (1988) On the Universality of Inversion-Layer Mobility in N- and P-Channel MOSFETs. Technical Digest, International Electron Devices Meeting, San Francisco, CA, 11-14 December 1988, 398-401.https://doi.org/10.1109/IEDM.1988.32840
- 34. Tsuno, M., Suga, M., Tanaka, M., Shibahara, K., Miura-Mattausch, M. and Hirose, M. (1999) Physically-Based Threshold Voltage Determination for MOSFETs of All Gate Lengths. IEEE Transactions on Electron Devices, 46, 1429-1434. https://doi.org/10.1109/16.772487
- 35. Jeong, C., Antoniadis, D. and Lundstrom, M.S. (2009) On Backscattering and Mobility in Nanoscale Silicon MOSFETs. IEEE Transactions on Electron Devices, 56, 2762-2769. https://doi.org/10.1109/TED.2009.2030844
- 36. Pantelides, S.T., Wang, S., Franceschetti, A., et al. (2006) Si/SiO
_{2}and SiC/SiO_{2 }Interfaces for MOSFETs-Challenges and Advances. Materials Science Forum, 527-529, 935-948. https://doi.org/10.4028/http://www.scientific.net/MSF.527-529.935 - 37. Krems, M. and Di Ventra, M. (2013) Ionic Coulomb Blockade in Nanopores. Journal of Physics: Condensed Matter, 25, Article ID: 065101. https://doi.org/10.1088/0953-8984/25/6/065101
- 38. Feng, J., Liu, K., Graf, M., Dumcenco, D., Kis, A., Di Ventra, M. and Radenovich, A. (2016) Observation of Ionic Coulomb Blockade in Nanopores. Nature Materials, 15, 850-855. https://doi.org/10.1038/nmat4607
- 39. Kaufman, I.K., McClintock, P.V.E. and Eisenberg, R.S. (2015) Coulomb Blockade Model for Permeation and Selectivity in Biological Ion Channels. New Journal of Physics, 17, Article ID: 083021. https://doi.org/10.1088/1367-2630/17/8/083021