Electronic Transport through DNA Nucleotides in Atomically Thin

ACS2GO © 2018. ← → → ←. loading. To add this web app to the home screen open the browser option menu and tap on Add to homescreen...
1 downloads 0 Views 980KB Size
Subscriber access provided by University of Otago Library

Biological and Medical Applications of Materials and Interfaces

Electronic Transport through DNA Nucleotides in Atomically Thin Phosphorene Electrodes for Rapid DNA Sequencing Rameshwar L Kumawat, Priyanka Garg, Sourabh Kumar, and Biswarup Pathak ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.8b17239 • Publication Date (Web): 12 Dec 2018 Downloaded from http://pubs.acs.org on December 17, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Electronic Transport through DNA Nucleotides in Atomically Thin Phosphorene Electrodes for Rapid DNA Sequencing Rameshwar L. Kumawat,† Priyanka Garg,# Sourabh Kumar,# Biswarup Pathak †,#,* †Discipline

of Metallurgy Engineering and Materials Science, #Discipline of Chemistry, School of Basic Sciences and, Indian Institute of Technology (IIT) Indore, Indore, Madhya Pradesh, 453552, India Email: [email protected] Abstract Rapid progresses in developing the fast, low-cost, and reliable methods for DNA sequencing are envisaged for development of personalized medicine. In this respect, nanotechnology has paved the role for the development of advanced DNA sequencing techniques including sequencing with solid-state nanopores or nanogaps. Herein, we have explored the application of a black phosphorene based nanogap-device for DNA sequencing. Using densityfunctional-theory (DFT) based non-equilibrium Green’s function (NEGF) approach, we have computed transverse transmission and current-voltage (I-V) characteristics of all the four DNA nucleotides (dAMP, dGMP, dTMP, and dCMP) as functions of applied bias voltages. We deduce that it is in principle; possible to differentiate between all the four nucleotides by three sequencing runs at distinct applied bias voltages, i.e., at 0.2 V, 1.4 V, and 1.6 V, where individual identification of all the four nucleotides may be possible. Hence, we believe our study might be helpful for experimentalist towards the development of a phosphorene based nanodevice for DNA sequencing to diagnose critical diseases. Keywords: DNA sequencing, black phosphorene, non-equilibrium Green’s function, density functional theory, I−V characteristics

1

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. Introduction Increasing demand for personalized medicine calls for a transformation in our health care system, which in turn demands for the development of new generation DNA sequencing based devices.1,2 In this respect, “$1000 Genome” project3 introduced by the National Institutes of Health (NIH) aims to develop personalized medicines4 based treatments by reducing the cost of the human genome sequence.5-8 Nanopore/Nanogap based DNA sequencing is one of the most growing third-generation DNA sequencing-based technologies that are presently foremost elegant to encounter such challenges as they have been growing rapidly to provide cheaper, faster, and accurate DNA sequencing.9-11 In 2012, Cherf and co-workers have reported that the information encoded in a DNA can be achieved through biological nanopores.12 However, in 2010, Schneider and co-workers and recently Garaj and co-workers have reported that the solid-state nanopore-based DNA sequencing is an emerging technology over biological nanopores.13,14 As such solid nanopore-based techniques help to achieve a single-nucleotide or base resolution.15 Furthermore, solid-state nanopores/nanogaps have better stability and high potential sequencing speeds over biological nanopores/nanogaps.16-18 In this regard, sequencing of the human genome (DNA) by solid-state nanopore/nanogap has progressed significantly. Several experimental and theoretical groups have used gold and various low dimensional materials based electrodes for rapid DNA sequencing.20-21 Graphene nanopores/nanogaps are excellent candidates to achieve single-nucleotide resolution due to its one atom-thickness.22 In a very recent report by Traversi and co-workers reported that the solid-state nanopore combined with the graphene nanoribbon can be used as a nanodevice for the nucleobase sensing.23 However, the surface properties of such graphene nanopores/nanogaps are very sensitive towards the target nucleotides. Moreover, such device

2

ACS Paragon Plus Environment

Page 2 of 21

Page 3 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

generates high noise ratio and geometry deformation after the insertion of a DNA molecule into the pore/gap.24,25 This has prompted the exploration of other potential low dimensional materials such as silicene (Si), molybdenum disulfide (MoS2), boron-nitride (BN), and graphene-hBN, to name a few.26-31 Two-dimensional (2D) monolayer black phosphorus, known as “Phosphorene”, has been considered as a new alternative to graphene for the next-generation semiconductor-based devices.32,48 It exhibits a puckered honeycomb-like structural geometry along the armchair direction. Phosphorene based materials have received a lot of attention due to their biocompatibility, low cytotoxicity, unique optical, electronic, and mechanical properties. For example, phosphorene has a high anisotropic electrical conductivity, carrier mobility (~1000 cm2V-1s-1), and a wide tunable band gap.32-34 Recent experimental studies have confirmed the stability of black phosphorus nanostructures in the air and water.35-37 There are several experimental reports on the stability of black phosphorus in air and water. Such experiential studies have reported that an effective encapsulation of phosphorene by other 2D materials39, such as h-BN,40 graphene,41 Poly(methyl methacrylate) (PMMA),42 and Al2O343, can prevent from degradation under ambient conditions. Such encapsulation technologies (packaging and passivation) lead to improved device reliability and performance, which is very necessary for its applications in current nano-electronics.38,44-48 For example, Pei and co-workers have demonstrated a highly controllable method for fabrication of high quality, air-stable phosphorene films with a selected number of layers ranging from a few down to the monolayer. After that in 2017, Cupo and co-workers have theoretically and experimentally studied the black phosphorene nano-hole/nanopore with different radius.58-59 Several studies show the potential biomedical applications of phosphorus nanostructure as it does not affect the protein structures.49 Therefore, we believe that these advances must facilitate its

3

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

development for DNA sequencing devices techniques.

All these reports indicate that

phosphorene can be a promising nanomaterial for bio-sensing30 applications. Furthermore, it would be interesting to find out whether the buckling nature in phosphorene is helpful for DNA sequencing or not. Inspired by all these reports, we have proposed phosphorene-based nanodevice for DNA sequencing. Using density-functional-theory (DFT) based non-equilibrium Green’s function (NEGF) approach,50-52 we have studied the transmission spectra and current as a function of bias voltage for the four nucleotides while passing through the phosphorene nanogap.

2. Model and Computational Methodology In this study, we have modeled a phosphorene nanodevice (Figure 1) for DNA sequencing where phosphorene used as electrodes as well as scattering region. The proposed device contains two phosphorene electrodes fixed oppositely in a nanogap and a target nucleotide as shown in Figure 1. A nanogap of ~1.45 nm is created along the armchair direction,53 as the electrical conductance in the armchair direction is one order higher than that in the zigzag direction.54 Such sized nanogap is good enough to allow the nucleotides (size ~ 1.0-1.2 nm) to pass through. The nucleotide is a part of single-stranded DNA (ssDNA) containing a base, a sugar, and a phosphate group. We have considered four target nucleotides for our study: (i) deoxy adenosine monophosphate (dAMP), (ii) deoxy guanidine monophosphate (dGMP), (iii) deoxy thymidine monophosphate (dTMP), and (iv) deoxy cytosine monophosphate (dCMP). The electronic calculations have been performed using the SIESTA code.50 Considering the optimized structures of armchair phosphorene nanogap and nucleotides, we placed each nucleotide inside the nanogap with H-terminated zigzag edges and the device is optimized.

4

ACS Paragon Plus Environment

Page 4 of 21

Page 5 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Left (L), and right (R) parts of the electrodes are periodic along the xz-plane. Furthermore, a periodic boundary condition has been adopted to avoid the electrostatic interactions. We have considered a large vacuum distance of ~16 Å, which is enough to avoid any interactions between the repeated images. We have used the GGA-PBE functional for the exchange-

Figure 1. Schematic of proposed phosphorene nanodevice with left (L), right (R) electrodes and scattering region. Top and side views of the fully optimized structures of the nanodevices with target nucleotides (a) dAMP, (b) dGMP, (c) dTMP, and (d) dCMP. Here, the transport direction is along the z-direction.

correlation-potentials. To model the atomic-core-electrons, the Troullier-Martins-normconserving51 pseudo-potentials are used, whereas double-zeta polarized (DZP) basis sets are used for valance electrons.52 A mesh cut off value of 200 Ry is utilized for the real-space integration. For the treatment of brillouin-zone sampling, a k-space grid of 4 × 1 × 2 has been employed. All the structures are optimized using the conjugate gradient algorithm (CG).

5

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 21

The convergence tolerance is set to 0.0001 eV, and the atomic forces are smaller than 0.01 eV/Å.31 The electronic transport calculations are done within the framework of the Landauer55,57 methodology (DFT+NEGF methods) as employed in the TranSiesta code.50 The basis-set and the real-space integration engaged in the transport calculations are same as used for geometry optimization. The current is calculated using the integration of the transmission spectrum as

2e μ

I(V𝑏) = h ∫μL T(E,V𝑏)[f(E ― μL) ― f(E ― μR)]dE R

(1)

where T(E,V𝑏) is the transmission probability of the electrons. E is the energy of the electrons travelling from L (left) to R (right) electrode in the presence of an applied bias voltage Vb. Here, f(E ― μL, R) and μL, R are the Fermi-Dirac distribution of electrons and chemicalpotential of the electrodes, respectively.

3. Results and Discussion The proposed device has phosphorene zigzag edges, which are terminated with H atoms. This has been done to improve the interactions between electrodes and nucleotides, which in turn can improve the transmission. We have calculated zero bias transmission for all the four nucleotides to understand the coupling between electrodes and nucleotides and Figure 2a shows the zero bias transmission for all the four nucleotides while placed between the electrodes. In the -0.5 to -2.0 eV energy region of the transmission curve, we find distinct peaks for all the four nucleotides. The maximum intensity comes from dAMP, dGMP and dTMP and lowest intensity peaks come from dCMP. This could be due to the highest occupied molecular orbital states (HOMO) of the device, which appears close to the Fermi energy. Similarly, in the 0.1 to 2.0 eV energy regions, peaks with high intensities appear for 6

ACS Paragon Plus Environment

Page 7 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

dGMP, dAMP and dTMP and low intensity for dCMP. This could be due to the contribution of the lowest unoccupied molecular orbitals (LUMO) of the device. Furthermore, we have done density of states analysis (Figure 2b) and plotted the molecular orbitals (Figure 2c) corresponds to these peaks appear around the Fermi energy. Our detailed molecular orbital analysis shows (Figure 2c) that HOMOs are mainly localized on the nucleotides, whereas LUMOs are localized on the electrodes. Figure 2c shows that there is a strong overlap

Figure 2. (a) Zero-bias transmission spectrum as a function of energy plotted for all the four DNA nucleotides (dAMP, dGMP, dTMP, and dCMP). (b) DOS at zero-bias plotted for the four DNA nucleotides. The Fermi level (Ef) has been set to zero. (c) Molecular orbitals responsible for these transmission peaks labels with letters A, T, G, and C have been presented for all the four devices. Color code: (P:orange), (O:red), (N:blue), (C:grey), and (H:white).

7

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

between electrode and nucleotide in case of dGMP, dTMP, and dAMP, where the overlap is less significant for the dCMP nucleotide. The zero-bias transmission as a function of energy is different for different nucleobases. This is because, there are two types of nucleobases: (i) purine (such as adenine and guanine), and (ii) pyrimidine based (cytosine and thymine). The pyrimidine based nucleobases are smaller in size than the purine based nucleobases and the number of hydrogen available for interactions are also different in these two cases. This is one of the reasons for weak coupling between dCMP and the phosphorene electrodes. As a result, pyrimidine base dCMP show less overlapping with phosphorene electrode. Our DOS plot in figure 2b also reveals that the orbital of dCMP does not significantly overlap with the phosphorene electrode. It can also be seen from figure 2c:C1 that the electron density is completely localized on the nucleobase which further confirms that there is a weak coupling between the phosphorene electrodes and dCMP. So, this is the reason that dCMP shows the lowest transmission intensity peaks. All these results indicate that all the three nucleotides (dGMP, dAMP, dTMP) are strongly coupling while dCMP shows somewhat weak coupling with the phosphorene electrodes. This is very much consistent with our transmission spectra. Furthermore, the existence of such peaks reveals that the DNA nucleotides can be identified under applied bias voltages, as these are the important factors for their tunneling transmittance. In addition, we have also investigated the effect of Hubbard U term (Ueffective = 1, 2, and 3) on the dGMP/dCMP + electrodes systems as they strongly/weakly couple with the phosphorene electrodes. However, we have not observed any change in their electronic structures (Figure S5, supporting information), which suggests that the single particle picture may be valid for such cases.

In order to electrically distinguish the four nucleotides, we have calculated the current signal for all the four DNA nucleotides (dAMP, dGMP, dTMP, and dCMP) while placed between

8

ACS Paragon Plus Environment

Page 8 of 21

Page 9 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

phosphorene electrodes. We have used up to 2V voltages to understand the current vs. voltage pattern for all the four nucleotides. Figure 3 shows the I-V pattern for a specific orientation of all the four nucleotides. We have used a flowchart (Figure S1) to show that we could distinguish them sequentially. The I-V pattern curve shows that the current increases as we increase the bias voltage. Besides, the four nucleotides can be distinguished into three different runs as shown in (Figure 3b-d, Figure S1). At 0.2 V bias voltage, the current (I) for dGMP is nearly 2 orders of amount higher than that for dCMP (Figure 3b). This further indicates the good sensitivity of the proposed nanodevice. Thereafter, at 1.4 V, the current for dGMP is nearly 1.5 orders of amount higher than that for dAMP or dTMP ( Figure 3c). Further, at 1.6 V, we have noted significant current differences between all the four DNA nucleotides for their individual identification. The magnitude of the currents at 1.6 V voltage (Figure 3d) is in the following order: IdGMP > IdTMP > IdAMP > IdCMP. Therefore at 1.6 V (Figure S2), we could distinguish dAMP and dTMP (dTMP > dTMP). Moreover, dGMP and dTMP nucleotides can be differentiated from the other two DNA nucleotides for their

9

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. (a) I-V pattern (current axis is presented on a logarithmic scale) curve for the four nucleotides (dAMP, dGMP, dTMP, and dCMP) while located in a nanogap. The bottom panel shows a view of the I-V pattern curve for the identification of nucleotides at three applied bias voltages, (b) 0.2 V, (c) 1.4 V, and (d) 1.6 V, used for identification of each DNA nucleotide. higher current values, which results in the closure of the Fermi-energy (Ef) of the proposed device to the broad HOMO peaks of dGMP and dTMP nucleotides (Figure 2). The other two nucleotides (dAMP, and dCMP), which have HOMO states further apart from the Fermienergy (Ef), unveil different current amounts, displaying tiny overlap with each other. Such differences in current between these two nucleobases are robust. In our investigation, the current increases the maximum for dGMP compared to the other three nucleotides. However, in all the four cases, the current values are very high (I > 10 ―10A), which is likely to avoid the electrical background noise in experimental investigations. Thus, the voltage window of 0.2 to 2 V will be good enough for a perfect differentiation of all the four DNA nucleotides. The 0.2-2 V bias window can be good for the device stability as graphene starts rupturing at 3.1 V.57 In this regard; we require a low bias voltage to distinguish all the four nucleotides. It is well known that the devices are driven out of the equilibrium condition under applied bias voltages. Therefore to understand their I-V pattern, we have analyzed the transmission spectra under applied biases (i.e. at 0.2, 1.4, and 1.6 V, respectively) for the four DNA nucleotides (Figure 4). These Figures show that the transmission peak shifts as we change the bias voltage. For example, all the transmission peaks move from high energy to low energy region as we increase the bias voltage. However, we find that the peak intensity is maximum for dGMP and dTMP, whereas low for dCMP. This is very much in agreement with our I-V curve too. Therefore, the observed molecular states of the target nucleotide and

10

ACS Paragon Plus Environment

Page 10 of 21

Page 11 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

the electrodes are contributing in the transmission spectrum, which shows significant changes under applied bias voltages. Hence, the I-V pattern analysis shows unique variation in the low applied bias voltages region. Thus, all the four DNA nucleotides convey unique current signs under applied biases.

Figure 4. Bias-dependent transmission spectra of the four devices (a) dAMP, (b) dGMP, (c) dTMP, and (d) dCMP) under applied bias voltages of 0.2, 1.4, and 1.6 V. The Fermi level has been set to zero. The broken lines depict the voltage window of ± 𝑉/2. Furthermore, in the context of our phosphorene nanogap electrode device performance, we have considered the effect of rotational and lateral translation for all the nucleobases on the transmission spectrums to account the dynamical effects. We understand that the transport will most likely be dominated by the tail of the distribution of transmission values and not 11

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

necessarily that the lowest energy structure will provide the highest conductance. For this, we have considered the rotations of all four DNA nucleobases around the y-axis in the xz plane at 60°, and 90°and the corresponding transmission spectrums are also shown in Figure S3. The relative energy of these rotated configurations with respect to the original (zero degrees) orientation is also presented in Table S1. By rotating the nucleotides inside the gap, the transmission spectrum of the device varies along with their orientations. Further, comparing these results against the case of zero-degree rotation, we find a variation in the transmission spectrums. Moreover, for the dGMP nucleobase the largest deviation in the transmission spectra corresponds to the orientation, which has the highest energy (60°). For all other case, there is a negligible difference due to rotation. In addition by laterally shifting the nucleobases (left and right side) by 0.5 Å, we note there is a change to the transmission peaks (Figure S4). Therefore, we concluded that there is some unavoidable effects of lateral translations on the transmission though it won’t effect much the I-V pattern observed for the ground state structure.15,31 Similarly, solution can play an important role. Recently, Feliciano et al., have shown the effect of solvent (water) medium for the zero-bias transmission of graphene-nucleobase system. The small changes observed in the transmission plots undoubtedly describe the effect of screening, assisted by the water molecules. However, this will only affect the magnitude of the transverse electrical current. Solvent medium can provide better or less transverse electrical current. But the trend observed within systems doesn’t change.60 In a similar report by Lagerqvist and co-workers, they have shown the structural fluctuations in the nucleotides causing the change in transmission function. However, they have concluded that the direct effect of water on transmission can be neglected.61,62 In conclusion, the trend in I-V pattern will not be changed due to the solvent effects though there can be change in the magnitude of the transverse electrical current.

12

ACS Paragon Plus Environment

Page 12 of 21

Page 13 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

4. Conclusion We have demonstrated the sensing application of phosphorene based nanogap-device for DNA sequencing. We find that phosphorene works as an ideal nano-electrode, and the zigzag H terminated phosphorene nanogap can interact with an incoming target nucleotide in a cooperative manner. By computing the I-V pattern and the transverse transmission spectrum of all the four target nucleotides in the nanodevice, we perceive that phosphorene nanogap based sequencing is indeed feasible to individually identify all the four DNA nucleotides. We find a noticeable difference in I-V signals into three different runs at 0.2, 1.4, and 1.6 V voltages, and observed the adequate sensitivity at 1.6 V. It is fascinating to observe that dAMP and dTMP show a comparatively smaller current value than the dGMP, which is significant enough to distinguish them. Moreover, dGMP and dCMP have nearly 2 orders of difference in current values. However, in the case of dGMP, the current trace is continuously increasing, whereas, in the case of dCMP, the current value is always lowest. This is due to the coupling of nucleobases with these electrodes; the induced density in the noticed nucleobases delivers a molecular level identification as noted over current-voltage. Hence, based on our results, we believe that it is possible to realize a phosphorene based nanodevice to distinguish between the four DNA nucleotides. 5. Supporting Information: Flowchart of I-V pattern to distinguish all the four nucleobases sequentially, The I-V pattern for individual identification of dAMP and dTMP at 1.6 V applied bias voltage, Effect of rotation and lateral translation of nucleobases on transmission spectra, Effect of Hubbard U term on the electronic structures of dGMP/dCMP. 6. Conflicts of interest There are no conflicts of interest to declare.

13

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

7. Acknowledgments We thank IIT Indore for the lab and computing facilities. This work is supported by DSTSERB, (Project Number: EMR/2015/002057) New Delhi and CSIR [Grant number: 01(2886)/17/EMR (II)]. R.L.K. and P.G. thanks MHRD for research fellowships. S.K. thanks UGC fellowship.

8. References (1) Saha, K. K.; Drndić, M.; Nikolić, B. K. DNA Base-Specific Modulation of Microampere Transverse Edge Currents through a Metallic Graphene Nanoribbon with a Nanopore. Nano Lett. 2012, 12, 50–55. (2) Heerema, S. J.; Dekker, C. Graphene Nanodevices for DNA Sequencing. Nat. Nanotechnol. 2016, 11, 127–136. (3) Check Hayden, E. The $1,000 genome. Nature 2014, 507 (7413), 294−295. (4) Steinbock, L. J.; Radenovic, A. The Emergence of Nanopores in Next-generation Sequencing. Nanotechnology 2015, 26 (7), 074003. (5) Mardis, E. R. A Decade’s Perspective on DNA Sequencing Technology. Nature 2011, 470 (7333), 198−203. (6) Rabbani, B.; Tekin, M.; Mahdieh, N. The Promise of Whole-exome Sequencing in Medical Genetics. J. Hum. Genet. 2014, 59 (1), 5−15. (7) Shendure, J.; Aiden, E. L. The Expanding Scope of DNA Sequencing. Nat. Biotechnol. 2012, 30 (11), 1084−1094. (8) Venkatesan, B. M.; Bashir, R. Nanopore Sensors for Nucleic Acid Analysis. Nat. Nanotechnol. 2011, 6, 615–624. (9) Branton, D.; Deamer, D. W.; Marziali, A.; Bayley, H.; Benner, S. A.; Butler, T.; Ventra, M. D.; Garaj, S.; Hibbs, A.; Huang, X.; Jovanovich, S. B.; Krstic, P. S.; Lindsay, S.; Ling, X.

14

ACS Paragon Plus Environment

Page 14 of 21

Page 15 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

S.; Mastrangelo, C. H.; Meller, A.; Oliver, J. S.; Pershin, Y. V.; Ramsey, J. M.; Riehn, R.; Soni, G. V.; Tabard-Cossa, V.; Wanunu, M.; Wiggin, M.; Schloss, J. A. The Potential and Challenges of Nanopore Sequencing. Nat. Biotechnol. 2008, 26, 1146–1153. (10) Nelson, T.; Zhang, B.; Prezhdo, O. V. Detection of Nucleic Acids with Graphene Nanopores: Ab initio Characterization of a Novel Sequencing Device. Nano Lett. 2010, 10, 3237–3242. (11) He, H.; Scheicher, R. H.; Pandey, R. Functionalized Nanopore-Embedded Electrodes for Rapid DNA Sequencing. J. Phys. Chem. C 2008, 112, 3456–3459. (12) Cherf, G. M.; Lieberman, K. R.; Rashid, H.; Lam, C. E.; Karplus, K.; Akeson, M. Automated Forward and Reverse Ratcheting of DNA in a Nanopore at 5-Å Precision. Nat. Biotechnol. 2012, 30, 344–348. (13) Schneider, G. F.; Kowalczyk, S. W.; Calado, V. E.; Pandraud, G.; Zandbergen, H. W.; Vandersypen, L. M. K.; Dekker, C. DNA Translocation through Grapheme Nanopores. Nano Lett. 2010, 10, 3163–3167. (14) Garaj, S.; Liu, S.; Golovchenko, J. A.; Branton, D. Molecule-hugging Graphene Nanopores. Proc. Natl. Acad. Sci. 2013, 110, 12192–12196. (15) Prasongkit, J.; Grigoriev, A.; Pathak, B.; Ahuja, R.; Scheicher, R. H. Transverse Conductance of DNA Nucleotides in a Graphene Nanogap from First Principles. Nano Lett. 2011, 11, 1941–1945. (16) Kim, H. S.; Kim, Y. H. Recent Progress in Atomistic Simulation of Electrical Current DNA Sequencing. Biosens. Bioelectron. 2015, 69, 186–198. (17) Wanunu, M. Nanopores: A Journey towards DNA Sequencing. Phys. Life Rev. 2012, 9, 125–158.

15

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(18) Kudr, J.; Skalickova, S.; Nejdl, L.; Moulick, A.; Ruttkay–Nedecky, B.; Adam, V.; Kizek, R. Fabrication of Solid-State Nanopores and its Perspectives. Electrophoresis 2015, 36, 2367–2379. (19) Dekker, C. Solid-state Nanopores. Nat. Nanotechnol. 2007, 2, 209–215. (20) Yang N.; Jiang X. Nanocarbons for DNA Sequencing: A Review. Carbon 2017, 115, 293-311. (21) Shendure J.; Ji H. Next-generation DNA Sequencing. Nat. Biotechnol. 2008, 26, 11351145. (22) Postma, H. W. Ch. Rapid Sequencing of Individual DNA Molecules in Graphene Nanogaps. Nano Lett. 2010, 10, 420–425. (23) Traversi, F.; Raillon, C.; Benameur, S. M.; Liu, K.; Khlybov, S.; Tosun, M.; Krasnozhon, D.; Kis, A.; Radenovic, A. Detecting the Translocation of DNA Through a Nanopore Using Graphene Nanoribbons. Nat. Nanotechnol. 2013, 8 (12), 939−945. (24) Merchant, C. A.; Healy, K.; Wanunu, M.; Ray, V.; Peterman, N.; Bartel, J.; Fischbein, M. D.; Venta, K.; Luo, Z.; Johnson, A. T. C.; Drndić, M. DNA Translocation through Graphene Nanopores. Nano Lett. 2010, 10 (8), 2915−2921. (25) Garaj, S.; Hubbard, W.; Reina, A.; Kong, J.; Branton, D.; Golovchenko, J. A. Graphene as a Subnanometre Trans-Electrode Membrane. Nature 2010, 467 (7312), 190−193. (26) Amorim, R. G.; Scheicher, R. H. Silicene as a New Potential DNA Sequencing Device. Nanotechnology 2015, 26, 154002. (27) Liu, K.; Feng, J.; Kis, A.; Radenovic, A. Atomically Thin Molybdenum Disulfide Nanopores with High Sensitivity for DNA Translocation. ACS Nano 2014, 8 (3), 2504−2511. (28) Farimani, A. B.; Min, K.; Aluru, N. R. DNA Base Detection Using a Single-Layer MoS2. ACS Nano 2014, 8, 7914–7922.

16

ACS Paragon Plus Environment

Page 16 of 21

Page 17 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(29) Liu, S.; Lu, B.; Zhao, Q.; Li, L.; Gao, T.; Chen, Y.; Zhang, Y.; Liu, Z.; Fan, Z.; Yang, F.; You, L.; Yu, D. Boron Nitride Nanopores: Highly Sensitive DNA Single-Molecule Detectors. Adv. Mater. 2013, 25, 4549-4554. (30) Gu, Z.; Zhang, Y.; Luan, B.; Zhou, R. DNA Translocation through Single-Layer Boron Nitride Nanopores. Soft Matter 2016, 12, 817-823. (31) Shukla, V.; Jena, N. K.; Grigoriev, A.; Ahuja, R. Prospects of Graphene−hBN Heterostructure Nanogap for DNA Sequencing. ACS Appl. Mater. Interfaces 2017, 9, 39945– 39952. (32) Lee G.; J. Pearton, S.; Ren, F.; Kim, J. Two-Dimensionally Layered p-Black Phosphorus/n-MoS2/p-Black

Phosphorus

Heterojunctions.

ACS

Appl.

Mater.

Interfaces 2018, 10 (12), 10347-10352.

(33) Castellanos-Gomez, A. Black Phosphorus: Narrow Gap, Wide Applications. J. Phys. Chem. Lett. 2015, 6, 4280–4291. (34) Kou, L.; Chen, C.; Smith, S. C. Phosphorene: Fabrication, Properties, and Applications. J. Phys. Chem. Lett. 2015, 6, 2794–2805. (35) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black Phosphorus Field-Effect Transistors. Nat. Nanotechnol. 2014, 9, 372–377. (36) Hanlon, D.; Backes, C.; Doherty, E.; Cucinotta, C. S.; Berner, N. C.; Boland, C.; Lee, K.; Harvey, A.; Lynch, P.; Gholamvand, Z.; Zhang, S.; Wang, K.; Moynihan, G.; Pokle, A.; Ramasse, Q. M.; McEvoy N.; Blau, W. J.; Wang, J.; Abellan, G.; Hauke, F.; Hirsch, A.; Sanvito, S.; O’Regan, D. D.; Duesberg, G. S.; Nicolosi, V.; Coleman, J. N. Liquid Exfoliation of Solvent-Stabilized Few-Layer Black Phosphorus for Applications Beyond Electronics. Nat. Commun. 2015, 6, 8563.

17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(37) Yasaei, P.; Kumar, B.; Foroozan, T.; Wang, C.; Asadi, M.; Tuschel, D.; Indacochea, J. E.; Klie, R. F.; Salehi-Khojin, A. High Quality Black Phosphorus Atomic Layers by LiquidPhase Exfoliation. Adv. Mater. 2015, 27, 1887−1892. (38) Kang, J.; Wood, J. D.; Wells, S. A.; Lee, J.-H.; Liu, X.; Chen, K. S.; Hersam, M. C. Solvent Exfoliation of Electronic-Grade, Two Dimensional Black Phosphorus. ACS Nano 2015, 9, 3596−3604. (39) Ling, X.; Wang, H.; Huang, S.; Xia, F.; Dresselhaus, M. S. The Renaissance of Black Phosphorus. Proc. Natl. Acad. Sci. 2015, 112 (15), 4523-4530. (40) Gillgren, N.; Wickramaratne, D.; Shi, Y.; Espiritu, T.; Yang, J.; Hu, J.; Wei, J.; Liu, X.; Mao, Z.; Watanabe, K.; Taniguchi, T.; Bockrath, M.; Barlas, Y.; Lake, R. K.; Lau, N. C. Gate Tunable Quantum Oscillations in Air-Stable and High Mobility Few-Layer Phosphorene Heterostructures. 2D Mater. 2015, 2, 011001. (41) Avsar A.; Vera-Marun I. J.; Tan J. Y.; Watanabe K.; Taniguchi T.; Castro Neto A. H.; O zyilmaz B. Air-Stable Transport in Graphene-Contacted, Fully Encapsulated Ultrathin Black Phosphorus-Based Field-Effect Transistors. ACS Nano, 2015, 9 (4), 4138-4145. (42) Koenig, S. P.; Doganov, R. A.; Schmidt, H.; Castro Neto, A. H.; Özyilmaz, B. Electric Field Effect in Ultrathin Black Phosphorus. Appl. Phys. Lett. 2014, 104(10), 103106. (43) Wood, J.D.; Wells, S.A.; Jariwala, D.; Chen, K.S.; Cho, E.; Sangwan, V.K.; Liu, X.; Lauhon, L.J.; Marks, T.J.; Hersam, M.C. Effective Passivation of Exfoliated Black Phosphorus Transistors Against Ambient Degradation. Nano lett. 2014, 14(12), 6964-6970. (44) Cai, Y.; Ke, Q.; Zhang, G.; Zhang, Y. W. Energetics, Charge Transfer, and Magnetism of Small Molecules Physisorbed on Phosphorene. J. Phys. Chem. C 2015, 119(6), 3102-3110. (45) Rojas, J. P.; Torres Sevilla, G. A.; Ghoneim, M. T.; Inayat, S. B.; Ahmed, S. M.; Hussain, A. M.; Hussain, M. M. Transformational Silicon Electronics. ACS nano 2014, 8(2), 1468-1474.

18

ACS Paragon Plus Environment

Page 18 of 21

Page 19 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(46) Okamoto, H.; Sugiyama, Y.; Nakano, H.; Synthesis and Modification of Silicon Nanosheets and Other Silicon Nanomaterials. Chem. Eur. J. 2011, 17, 9864–9887. (47) Zwanenburg, F. A.; Dzurak, A. S.; Morello, A.; Sim-mons, M. Y.; Hollenberg, L.C.L.; Klimeck, G.; Rogge, S.; Coppersmith, S. N.; Eriksson, M. A. Silicon Quantum Electronics. Rev. Mod. Phys. 2013, 85, 961-1019. (48) Sorkin, V.; Cai, Y.; Ong, Z.; Zhang, G.; Zhang, Y. W. Recent Advances in the Study of Phosphorene and Its Nanostructures. Critical Reviews in Solid State and Materials Sciences, 2017, 42(1), 1-82. (49) Cortés-arriagada, D. Phosphorene as a Template Material for Physisorption of DNA / RNA Nucleobases and Resembling of Base Pairs: A Cluster DFT Study and Comparisons with Graphene. J. Phys. Chem. C 2018, 122, 4870-4880. (50) Brandbyge, M.; Mozos, J. L.; Ordejón, P.; Taylor, J.; Stokbro, K. Density-Functional Method for Nonequilibrium Electron Transport. Phys. Rev. B - Condens. Matter Mater. Phys. 2002, 65 (16), 165401. (51) Pathak, B.; Lofa˚s, H.; Prasongkit, J.; Grigoriev, A.; Ahuja, R.; Scheicher, R. H. DoubleFunctionalized Nanopore-Embedded gold Electrodes for Rapid DNA Sequencing. Appl. Phys. Lett. 2012, 100, 023701. (52) Prasongkit, J.; Grigoriev, A.; Pathak, B.; Ahuja, R.; Scheicher, R. H. Theoretical Study of Electronic Transport through DNA Nucleotides in a Double-Functionalized Graphene Nanogap. J. Phys. Chem. C 2013, 117, 15421–15428. (53) Nelson, T.; Zhang, B.; Prezhdo, O. V. Detection of Nucleic Acids with Graphene Nanopores : Ab Initio Characterization of a Novel Sequencing Device. Nano Lett. 2010, 10, 3237–3242.

19

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(54) Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J. A.; Lo, C.; Yang, L. Enhanced Thermoelectric Efficiency via Orthogonal Electrical and Thermal Conductances in Phosphorene. Nano Lett. 2014, 14, 6393−6399. (55) Buttiker, M.; Imry, Y.; Landauer, R.; Pinhas, S. Generalized Many-Channel Conductance Formula with Application to Small Rings. Phys. Rev. conductance formula. 1985, 31, 6207–6215. (56) Martin T.; Landauer, R. Wave-Packet Approach to Noise in Multichannel Mesoscopic System. Phys. Rev. B 1992, 45, 1742. (57) Barreiro, A.; Börrnert, F.; Rümmeli, M. H.; Büchner, B.; Vandersypen, L. M. K. Graphene at High Bias: Cracking, Layer by Layer Sublimation, and Fusing. Nano Lett. 2012, 12 (4), 1873−1878. (58) Pei, J.; Gai, X.; Yang, J.; Wang, X.; Yu, X.; Choi, D.K.; Luther-Davies, B.; Lu, Y. Producing Air-Stable Monolayers of Phosphorene and Their Defect Engineering. Nat. Commun. 2016, 7, 10450. (59) Cupo, A.; Das P.M.; Chine, C.-C.; Danda, G.; Kharche, N.; Tristant, D.; Drndic, M.; Meunier, V. Periodic Arrays of Phosphorene Nanopores as Antidot Lattices with Tunable Properties. ACS Nano 2017, 11, 7494-7507. (60) Feliciano, G. T.; Sanz-Navarro, C.; Coutinho-Neto, M. D.; Ordejon, P.; Scheicher, R.H.; Rocha, A. R. Addressing the Environment Electrostatic Effect on Ballistic Electron Transport in Large Systems: A QM/MM-NEGF Approach. J. Phys. Chem. B 2018, 122, 2, 485-492. (61) Lagerqvist, J.; Zwolak, M.; Ventra, M. D. Fast DNA Sequencing via Transverse Electronic Transport. Nano Lett. 2006, 6, 4, 779-782.

20

ACS Paragon Plus Environment

Page 20 of 21

Page 21 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

(62) Lagerqvist, J.; Zwolak, M.; Ventra, M. D. Influence of the Environment and Probes on Rapid DNA Sequencing via Transverse Electronic Transport. Biophys J. 2007, 93(7), 2384– 2390.

Table of Content (TOC):

21

ACS Paragon Plus Environment