Quantum Simulation of the Silicene and Germanene for Sensing and Sequencing of DNA/RNA Nucleobases

Over the last decade, we have been witnessing the rise of two-dimensional (2D) materials. Several 2D materials with outstanding properties have been theoretically predicted and experimentally synthesized. 2D materials are good candidates for sensing and detecting various biomolecules because of their extraordinary properties, such as a high surface-to-volume ratio. Silicene and germanene are the monolayer honeycomb structures of silicon and germanium, respectively. Quantum simulations have been very effective in understanding the interaction mechanism of 2D materials and biomolecules and may play an important role in the development of effective and reliable biosensors. This article focuses on understanding the interaction of DNA/RNA nucleobases with silicene and germanane monolayers and obtaining the possibility of using silicene and germanane monolayers as a biosensor for DNA/RNA nucleobases’ sequencing using the first principle of Density Functional Theory (DFT) calculations with van der Waals (vdW) correction and nonequilibrium Green’s function method. Guanine (G), Cytosine (C), Adenine (A), Thymine (T), and Uracil (U) were examined as the analytes. The strength of adsorption between the DNA/RNA nucleobases and silicene and germanane is G > C > A > T > U. Moreover, our recent work on the investigation of Au- and Li-decorated silicene and germanane for detection of DNA/RNA nucleobases is presented. Our results show that it is possible to get remarkable changes in transmittance due to the adsorption of nucleobases, especially for G, A, and C. These results indicate that silicene and germanene are both good candidates for the applications in fast sequencing devices for DNA/RNA nucleobases. Additionally, our present results have the potential to give insight into experimental studies and can be valuable for advancements in biosensing and nanobiotechnology.

Over the recent years, graphene counterparts such as silicene and germanene hav gained valuable attention in the scientific community. Silicene was successfully grown experimentally on Ag [2,[30][31][32], Ir [33], and ZrB2 [34] substrate and germanene were also grown on the Pt (111) surface [35]. Although silicene and germanene stand out with thei properties similar to graphene, they are different from graphene in terms of their struc tural properties. They have a buckled honeycomb crystal geometry (see Figure 1b,c) which provides for sp 3 bonding. Because of their crystal structure, interaction of molecule and atoms with silicene and germanane monolayers is stronger than with graphene. Thi makes them potential candidates for the applications of sensing devices. Despite the importance of the interactions of 2D materials with biomolecules, ther are still limited studies with silicene and germanene. A covalent bonding mechanism be tween DNA/RNA nucleobases with silicene and germanene is absent in similar previou studies. Quantum simulations have been very effective in understanding the interaction mechanism of 2D materials and biomolecules and may play an important role in the de velopment of effective and reliable biosensors. With the increase and widespread use o high-performance computing resources, it has become clear how powerful the first prin ciple of the Density Functional Theory (DFT) is in modeling such systems. Metal-doped monolayers have better binding capabilities compared to pristine monolayers [17]. Gold atoms have some advantages over other metal atoms. The color of gold nanoparticles de pends on their shape and size; they are highly preferred in nanomedicine.
In this work, we have calculated the binding of DNA/RNA nucleobases (NB) with pristine silicene/germanene and Au-and Li-decorated silicene/germanene by using th quantum simulations technique based on the Density Functional Theory (DFT), with van der Waals (vdW) correction and electronic transport properties, such as changes in th transmittance and conductance produced by adsorption of the nucleobases for the firs time. Despite the importance of the interactions of 2D materials with biomolecules, there are still limited studies with silicene and germanene. A covalent bonding mechanism between DNA/RNA nucleobases with silicene and germanene is absent in similar previous studies. Quantum simulations have been very effective in understanding the interaction mechanism of 2D materials and biomolecules and may play an important role in the development of effective and reliable biosensors. With the increase and widespread use of high-performance computing resources, it has become clear how powerful the first principle of the Density Functional Theory (DFT) is in modeling such systems. Metaldoped monolayers have better binding capabilities compared to pristine monolayers [17]. Gold atoms have some advantages over other metal atoms. The color of gold nanoparticles depends on their shape and size; they are highly preferred in nanomedicine.
In this work, we have calculated the binding of DNA/RNA nucleobases (NB) with pristine silicene/germanene and Au-and Li-decorated silicene/germanene by using the quantum simulations technique based on the Density Functional Theory (DFT), with van der Waals (vdW) correction and electronic transport properties, such as changes in the transmittance and conductance produced by adsorption of the nucleobases for the first time.

Materials and Methods
The interactions' mechanism of nucleobases with the pristine and Au/Li-decorated silicene/germanene has been investigated using the Density Functional Theory (DFT)based SIESTA package, where the basis set was used as a linear combination of numerical atomic orbitals (LCNAO) [36]. Electronic transport calculations were performed with nonequilibrium Green's function (NEGF) formalism as implemented in the TranSiesta code [37]. The Perdew, Burke, and Ernzerhof functional within the generalized gradient approximation was used for the exchange-correlation functional [38]. The pseudopotentials that we used in this study are of the norm-conserving kind [39]. Conjugate gradient methods were used to optimize the geometry of the systems to get minimum energy configurations. Monkhorst-Pack meshes for Brillouin zone sampling were taken as 3 × 3 × 1 and 9 × 9 × 1 k-points for the relaxations and band structure calculations, respectively [40]. Furthermore, the Grimme-D2 method was used to describe long-range interactions [41]. For electronic transport calculations in the scattering region, the supercell had a size of 15.6 A 0 × 20 A 0 × 47.4 A 0 , while for the electrodes, we used a supercell with the size of 15.6 A 0 × 20 A 0 × 13.5 A 0 . A grid of 10 × 1 × 30 k-points was used for electronic transport calculations in the scatter region, while for the electrodes, the grid was composed of 60 × 1 × 30 k-points. The transport direction is aligned with the z-axis. The value of the mesh cut-off is taken as 300 Ry. We have employed periodic boundary conditions (PBC), so we have to give a space between two graphene layers in the adjacent supercells within 20 A 0 . Force tolerance was set to 0.04 eV/A 0 for convergence. Due to the heavy metal Au atom, we have included a spin-orbit (SO) contribution to the calculations. Visualization of the systems has been performed by using VESTA software [42].
In the electronic transport calculation, the system was divided into three regions: two semi-infinite electrodes and one central scattering region. In our calculations, the electrodes were taken as pristine silicene for the right-and left-hand sides (see Figure 2). The scattering region consisted of the silicene layer and the DNA/RNA nucleobases.
The interactions' mechanism of nucleobases with the pristine and Au/Li-dec silicene/germanene has been investigated using the Density Functional Theory based SIESTA package, where the basis set was used as a linear combination of nu atomic orbitals (LCNAO) [36]. Electronic transport calculations were performe nonequilibrium Green's function (NEGF) formalism as implemented in the Tra code [37]. The Perdew, Burke, and Ernzerhof functional within the generalized g approximation was used for the exchange-correlation functional [38]. The pseudo tials that we used in this study are of the norm-conserving kind [39]. Conjugate g methods were used to optimize the geometry of the systems to get minimum ener figurations. Monkhorst-Pack meshes for Brillouin zone sampling were taken as 3 and 9 × 9 × 1 k-points for the relaxations and band structure calculations, respective Furthermore, the Grimme-D2 method was used to describe long-range interactio For electronic transport calculations in the scattering region, the supercell had a 15.6 A 0 × 20 A 0 × 47.4 A 0 , while for the electrodes, we used a supercell with the size A 0 × 20 A 0 × 13.5 A 0 . A grid of 10 × 1 × 30 k-points was used for electronic transpor lations in the scatter region, while for the electrodes, the grid was composed of 60 k-points. The transport direction is aligned with the z-axis. The value of the mesh is taken as 300 Ry. We have employed periodic boundary conditions (PBC), so we give a space between two graphene layers in the adjacent supercells within 20 A tolerance was set to 0.04 eV/A 0 for convergence. Due to the heavy metal Au atom, w included a spin-orbit (SO) contribution to the calculations. Visualization of the s has been performed by using VESTA software [42].
In the electronic transport calculation, the system was divided into three regio semi-infinite electrodes and one central scattering region. In our calculations, th trodes were taken as pristine silicene for the right-and left-hand sides (see Figure  scattering region consisted of the silicene layer and the DNA/RNA nucleobases. where SS and HS are the overlap matrix and Hamiltonian, respectively, for the sca gion. The self-energies are given as Σ , where are the surface G The transmission can be obtained via the Green's functions for this system where the coupling matrices are given by Γ α = i Σ α − Σ † α with α ≡ {L, R}. Σ L/R are the self-energies of the semi-infinite electrodes and where S S and H S are the overlap matrix and Hamiltonian, respectively, for the scatter region. The self-energies are given as Σ α = V Sα g α V αS , where g α are the surface Green's functions for the semi-infinite leads and V Sα = V † αS are the coupling matrix elements between the scattering region and the electrodes. The Hamiltonian H S was obtained by the previous DFT calculations. Further information regarding NEGF can be found in [37,43].

Results and Discussion
Our recent results are detailed in the following sections. First of all, we investigated the binding mechanism of DNA/RNA nucleobases with bare silicene/germanene monolayers. By functionalizing Silicene/Germanene with Li and Au atoms, we have demonstrated that we can strongly binded Nuclebases to the functionalized surface. The obtained results indicate that Li and Au adatoms on silicene/germanene can significantly increase the binding energy of DNA/RNA nucleobases.

Binding Mechanism of Bare Silicene and Germanene Monolayers with Nucleobases
The geometric relaxation covers both different absorption heights and orientations of DNA/RNA NB on the silicene and germanene. We examined the several possible orientations of NB on silicene and germanene layers; the most favorable structures of NB on silicene and germanene sheets are shown in Figures 3 and 4, respectively.

Results and Discussion
Our recent results are detailed in the following sections. First of all, we investigated the binding mechanism of DNA/RNA nucleobases with bare silicene/germanene monolayers. By functionalizing Silicene/Germanene with Li and Au atoms, we have demonstrated that we can strongly binded Nuclebases to the functionalized surface. The obtained results indicate that Li and Au adatoms on silicene/germanene can significantly increase the binding energy of DNA/RNA nucleobases.

Binding Mechanism of Bare Silicene and Germanene Monolayers with Nucleobases
The geometric relaxation covers both different absorption heights and orientations of DNA/RNA NB on the silicene and germanene. We examined the several possible orientations of NB on silicene and germanene layers; the most favorable structures of NB on silicene and germanene sheets are shown in Figures 3 and 4, respectively.  It has been reported in earlier studies that being parallel is energetically more favorable because of the larger effective contact area for graphene [20]. According to our latest DFT calculations, all nucleobases have a tilted orientation on silicene. On the other hand, G, C, and T have a tilted orientation on germanene, but U and A show a parallel orientation. The tilting orientation provides sp 3 bonding to strengthen chemical bonding between biomolecules and monolayers. We show that there is a chemisorption between G, C, A, and silicene. The oxygen atom of G and C and the nitrogen atom of A interact with the silicone atom of silicene, and they form a chemical bond. Our results are in agree-Biosensors 2021, 11, 59 5 of 13 ment with the previous reported data except for A on silicene [14,20]. Amaorim et al. [14] and Hussain et al. [20] reported that A has a tilted orientation on silicene with no bonding. Our calculations predict bonding between the nitrogen atom of A and the silicone atom of silicene. However, there is no chemical bonding between nucleobases and germanene monolayers. It has been reported in earlier studies that being parallel is energetically more favorable because of the larger effective contact area for graphene [20]. According to our latest DFT calculations, all nucleobases have a tilted orientation on silicene. On the other hand, G, C, and T have a tilted orientation on germanene, but U and A show a parallel orientation. The tilting orientation provides sp 3 bonding to strengthen chemical bonding between biomolecules and monolayers. We show that there is a chemisorption between G, C, A, and silicene. The oxygen atom of G and C and the nitrogen atom of A interact with the silicone atom of silicene, and they form a chemical bond. Our results are in agreement with the previous reported data except for A on silicene [14,20]. Amaorim et al. [14] and Hussain et al. [20] reported that A has a tilted orientation on silicene with no bonding. Our calculations predict bonding between the nitrogen atom of A and the silicone atom of silicene. However, there is no chemical bonding between nucleobases and germanene monolayers.
The binding energy, as shown in Table 1, can be defined, EBinding, as follows: In terms of the total energies of the pristine silicene/germanene = monolayers (Li/Audecorated, in some cases) and of the NB (NB = Adenine (A), Guanine (G), Cytosine (C), Thymine (T) and Uracil (U)), respectively. The binding energy can be taken as an indicator of the strength of the interaction between biomolecules and monolayers.
Our recent results show that adsorption energies are in the range from −0.58 eV to −1.35 eV for silicene and −0.65 eV to −1.30 eV for germanene. The order of adsorption energies of the NB with silicene/germanene is G > C > A > T > U. We can conclude from the results that the interaction of Guanine is stronger than the other nucleobases for both bare silicene and germanene. Uracil also has the lowest binding with bare silicene and germanene. The order of binding energies of nucleobases on silicene and germanene is  Table 1, can be defined, E Binding , as follows: In terms of the total energies of the pristine silicene/germanene = monolayers (Li/Audecorated, in some cases) and of the NB (NB = Adenine (A), Guanine (G), Cytosine (C), Thymine (T) and Uracil (U)), respectively. The binding energy can be taken as an indicator of the strength of the interaction between biomolecules and monolayers.
Our recent results show that adsorption energies are in the range from −0.58 eV to −1.35 eV for silicene and −0.65 eV to −1.30 eV for germanene. The order of adsorption energies of the NB with silicene/germanene is G > C > A > T > U. We can conclude from the results that the interaction of Guanine is stronger than the other nucleobases for both bare silicene and germanene. Uracil also has the lowest binding with bare silicene and germanene. The order of binding energies of nucleobases on silicene and germanene is consistent with the previously reported data [20]. In the case of silicene, Amaorim et al. [14] reported that C has the strongest binding among the nucleobases, and that the binding energy of A is equal to T. The binding energy order is different from that of other 2D materials, such as graphene, MoS 2, and WS 2 (G > A > T > C > U), calculated by the DFT-D2 method [18]. The difference possibly comes from the buckled honeycomb structure of silicene and germanene.
The closest adsorption distances between the monolayer atom/adatom and the closest pair of atoms between the silicene/germanene monolayers and nucleobases are shown in Tables 2 and 3. Considering the most favorable adsorption geometries, in some cases, the reported distances correspond to bond length. There may be a correlation between binding energies and distances. More specifically, distances between Li/Au-supported germanene and the nucleobases are smaller than the distances between pristine germanene and the nucleobases.

Binding Mechanism of Li/Au-Doped Silicene and Germanene Monolayers with Nucleobases
According to our present results, in terms of binding energy, we confirm that the Li atom is chemisorbed onto the surface of silicene and germanene monolayers. The preferable adsorption site of the Au and Li atoms on silicene and germanene monolayers is a hollow site, see Figure 1a. Our calculated binding energy of the Li(Au) atom is −2.51 (−2.13) eV and −2.55 (−2.11) eV for silicene and germanene, respectively. Li and Au binding values are in good agreement with previous studies [44][45][46][47][48][49].
M. Sadeghi et al. [25] showed that the Li atom can improve the binding energy between the nucleobases and, MoS 2 . Song et al. [50] showed that the Au atom can increase the binding energy between the aromatic amino acids such as phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp) on MoS 2 monolayer. Additionally, Gorkan et al. [17] and Kadioglu et al. [51] reported that interactions of amino acids and nucleobases with Au clusters supported blue and black phosphorene, respectively.
In this part, we have calculated the binding of five nucleobases with the Li/Au adatom adsorbed to silicene and germanene monolayers. Li is an alkali metal and alkali metals are known as highly reactive metals. The equilibrium adsorption geometries with minimum total energies are shown in Figures 3 and 4. The oxygen atom of G and C, the nitrogen atom of A, and one of the oxygen atoms of T and U interact with the Li atom on silicene/germanene, and they form a covalent bonding.
According to our recent calculations for Li (Au), the binding energies range from −1.45 (−1.01) eV to −2.30 (−1.84) eV and −1.53 (−1.12) eV to −2.28 (−1.75) eV for silicene and germanene, respectively. Although U has the lowest binding energy for both bare silicene and germanene, the binding energy of U increases 150% and 135% on Li-supported silicene and germanene, respectively. G has the strongest binding energy among the nucleobases; both bare and Li/Au-supported silicene and germanene.
This energy range is much more diverse than that achieved in pristine silicene and germanene monolayers. We say that the adsorbed Li and Au atom to pristine silicene/germanene monolayers can greatly increase the binding between the nucleobases and silicene/germanene monolayers. Li/Au + silicene/germanene monolayers preserve the order of binding energies with the NB, and it is G > C > A > T > U. This is easily seen in the optimized structure in Figure 3c,d, the Li + Silicene/Cytosine and Li + Silicene/Thymine, where the Li atom moves toward the bridge and valley site, respectively. Diffusion of the Li atom can be explained by Coulomb interaction between the Li atom and C/T. This is shown in the optimized structure in Figure 4b, the Au + Germanene/Guanine, where the Au atom diffuses to the next hollow site. Diffusion of the Au atom can be explained by the Coulomb interaction between the Au atom and G. Unlike Li on silicene/germanane, the nitrogen atom of G and C interacts with the Au atom on silicene/germanene, and they form covalent bonding.

Electronic Properties of Li/Au-Doped Silicene and Germanene Monolayers with Nucleobases
We have investigated the effects of adsorption of Li/Au atoms in the electronic properties of silicene and germanene monolayers. In Figure 1e,f it is shown that pristine silicene and germanene are naturally semimetal and have a Dirac cone in K-point (see Figure 1d). Adsorption of the DNA/RNA nucleobases changes the band structure of bare silicene and germanene. In Figures 5a and 6a,f, as examples of the DNA/RNA nucleobases, it is shown that there is a bandgap opening in Silicene + G, Silicene + C, Germanene + G, and Germanene + C. Adsorption of the DNA/RNA nucleobases on bare silicene and germanene causes a bandgap opening that ranges from 22 meV to 107 meV. We showed the maximum bandgap opening of silicene and germanene in Figures 5 and 6, respectively. As a result of the adsorption of a Li and an Au atom, complexes of the semiconductor silicene and germanene/nucleobase complexes become metallic. This is due to the charge transfer from Li/Au atoms to silicene/germanene. Moreover, there is a small bandgap opening in the band below the Fermi level, and it occurs at the K symmetry point (see Figure 1d). For Au + Silicene/C see Figure 4e and for Au + Germanene/C see Figure 5e; there is a spin-orbit coupling that produces a band splitting. manene causes a bandgap opening that ranges from 22 meV to 107 meV. We showed the maximum bandgap opening of silicene and germanene in Figures 5 and 6, respectively. As a result of the adsorption of a Li and an Au atom, complexes of the semiconductor silicene and germanene/nucleobase complexes become metallic. This is due to the charge transfer from Li/Au atoms to silicene/germanene. Moreover, there is a small bandgap opening in the band below the Fermi level, and it occurs at the K symmetry point (see Figure 1d). For Au + Silicene/C see Figure 4e and for Au + Germanene/C see Figure 5e; there is a spin-orbit coupling that produces a band splitting.

Electronic Transport Properties
Here, we show the functionality of silicene and germanene for the DNA/RNA nucleobases' sequencing. When the DNA/RNA nucleobases are adsorbed on the surface of the silicene, electronic transport properties of silicene/germanene will change. Transmittance T(E) shows the probability of an electron to be passed from one electrode to the other through the scattering region in the middle. In order to obtain electron transport proper-

Electronic Transport Properties
Here, we show the functionality of silicene and germanene for the DNA/RNA nucleobases' sequencing. When the DNA/RNA nucleobases are adsorbed on the surface of the silicene, electronic transport properties of silicene/germanene will change. Transmittance T(E) shows the probability of an electron to be passed from one electrode to the other through the scattering region in the middle. In order to obtain electron transport properties, the transmission at zero bias voltage of the silicene/germanene attached to the semi-infinite electrodes is calculated, as shown in Figure 7a,d. Interaction of the DNA/RNA nucleobases with silicene/germanene monolayers causes reduction of energy-resolved transmission. The sensitivity histogram of bare silicene/germanene (left and middle panel) and Lisupported silicene/germanene (right panel) is shown in Figures 8 and 9, respectively. Our recent analysis shows that different conductance data from measurements at corresponding gate voltages allow distinguishing the nucleobases. At Vg = 0.62 V, it is possible to distinguish A from T and U. At Vg = −2.14 V, the theoretical silicene biosensor shows larger sensitivity for G and C and smaller sensitivity for U. At Vg = 2.70 V, and the bare silicene biosensor exhibits larger sensitivity for G, A, and C, while having smaller sensitivity for T and U. However, it is not possible to distinguish G, A, and C. A Li-supported silicene biosensor allows us to distinguish A from G and C for the same gate voltage. According to the sensitivity analysis of bare silicene distinguishing of T and U, it would not be as easy as G, C, and A. This problem can be solved by using a Li-supported silicene biosensor device. We note that the Li-supported silicene biosensor device can have be much more effective than the bare silicene biosensor device in distinguishing the DNA/RNA nucleobases. The main aim of this study is sensing and sequencing the DNA/RNA nucleobases by using silicene and germanene monolayers. We can design a theoretical biosensor for sensing the DNA/RNA nucleobases based on the variation in conductance. The sensitivity can be expressed as S(%) = |g−g 0 | g 0 , where g and g 0 are zero bias conductance of the biosensor with and without a nucleobase, respectively. For a small bias voltage, the conductance can be expressed as a function of gate voltage V g , g V g = g 0 T(E = µ), where g 0 = 2e 2 h is the quantum conductance, µ = E F − eV g is the chemical potential, h is the Planck's constant, and e is the charge of the electron. The applied gate voltage (V g ) can tune the chemical potential to variation in the transmission. Measuring conductance rather than the current allows a better distinction of the change in the signal [52] because conductance through the biosensor changes by a unit of quantum conductance g 0 in its interaction with the DNA/RNA nucleobases.
The sensitivity histogram of bare silicene/germanene (left and middle panel) and Li-supported silicene/germanene (right panel) is shown in Figures 8 and 9, respectively. Our recent analysis shows that different conductance data from measurements at corresponding gate voltages allow distinguishing the nucleobases. At V g = 0.62 V, it is possible to distinguish A from T and U. At V g = −2.14 V, the theoretical silicene biosensor shows larger sensitivity for G and C and smaller sensitivity for U. At V g = 2.70 V, and the bare silicene biosensor exhibits larger sensitivity for G, A, and C, while having smaller sensitivity for T and U. However, it is not possible to distinguish G, A, and C. A Li-supported silicene biosensor allows us to distinguish A from G and C for the same gate voltage. According to the sensitivity analysis of bare silicene distinguishing of T and U, it would not be as easy as G, C, and A. This problem can be solved by using a Li-supported silicene biosensor device. We note that the Li-supported silicene biosensor device can have be much more effective than the bare silicene biosensor device in distinguishing the DNA/RNA nucleobases. According to the sensitivity histogram of the theoretical biosensor based on germanene, there is a larger sensitivity for G and C, as expected. Distinguishing T and U would not be as easy as G and C for the pristine germanane biosensor, but it is possible to distinguish T with Li-supported germanene.

Conclusions
In this manuscript, the adsorption behaviors of the DNA/RNA nucleobases on pristine silicene/germanane monolayer surfaces and Li/Au-supported silicene/germanene monolayers were investigated using DFT with vdW correction and the nonequilibrium Green's function method. Our research covers all five DNA/RNA nucleobases in a nonaqueous environment. Our recent results show the sensing capability of pristine silicene/germanene and Li/Au-supported silicene/germanene and express their binding  According to the sensitivity histogram of the theoretical biosensor based on germanene, there is a larger sensitivity for G and C, as expected. Distinguishing T and U would not be as easy as G and C for the pristine germanane biosensor, but it is possible to distinguish T with Li-supported germanene.

Conclusions
In this manuscript, the adsorption behaviors of the DNA/RNA nucleobases on pristine silicene/germanane monolayer surfaces and Li/Au-supported silicene/germanene monolayers were investigated using DFT with vdW correction and the nonequilibrium Green's function method. Our research covers all five DNA/RNA nucleobases in a nonaqueous environment. Our recent results show the sensing capability of pristine silicene/germanene and Li/Au-supported silicene/germanene and express their binding mechanisms. Au/Li atoms strongly interact with silicene and germanane monolayers. Silicene and germanene accept electrons from Au/Li adatoms, according to their greater electronegativity. The order of adsorption between the DNA/RNA nucleobases and silicene/germanane monolayers is G > C > A > T > U. The interaction of T and U with pristine According to the sensitivity histogram of the theoretical biosensor based on germanene, there is a larger sensitivity for G and C, as expected. Distinguishing T and U would not be as easy as G and C for the pristine germanane biosensor, but it is possible to distinguish T with Li-supported germanene.

Conclusions
In this manuscript, the adsorption behaviors of the DNA/RNA nucleobases on pristine silicene/germanane monolayer surfaces and Li/Au-supported silicene/germanene mono-layers were investigated using DFT with vdW correction and the nonequilibrium Green's function method. Our research covers all five DNA/RNA nucleobases in a nonaqueous environment. Our recent results show the sensing capability of pristine silicene/germanene and Li/Au-supported silicene/germanene and express their binding mechanisms. Au/Li atoms strongly interact with silicene and germanane monolayers. Silicene and germanene accept electrons from Au/Li adatoms, according to their greater electronegativity. The order of adsorption between the DNA/RNA nucleobases and silicene/germanane monolayers is G > C > A > T > U. The interaction of T and U with pristine silicene and germanene can be described as weak and noncovalent. T and U are physisorbed, while G, C, and A are chemisorbed on silicene. Adsorbents (Li/Au atoms) did not change the order of binding. It is clear that the adsorbed Li/Au atoms can significantly increase the binding between the nucleobases and silicene/germanene monolayers. The calculated orientations and bonding mechanism of the DNA/RNA nucleobases on silicene and germanene may play an important role for experimental studies to identify them by using scanning tunneling microscopy (STM). Adsorption of the DNA/RNA nucleobases on pristine silicene and germanene causes a bandgap opening at the K symmetry point. Li-and Au-supported silicene/germanene + nucleobases complexes become metallic. Additionally, in the case of C on Au-supported silicene/germanene, band splitting occurs below the Fermi level because of spin-orbit coupling. We also focused on the changes in electronic transport properties of bare and Li-supported silicene/germanane due to adsorption of the DNA/RNA nucleobases. We simulated that it is possible to distinguish G, C, and A from T and U due to changes in electronic signals. For the specific gate voltages, changes in transmission caused by the nucleobases differ from those of pristine and Li-supported silicene/germanene. Thus, it is possible to distinguish all the DNA/RNA nucleobases. Our findings show that both bare and metal-doped silicene and germanene are possible candidates for DNA/RNA sequencing biosensor applications.
Strong binding capacities of both bare and metal-doped silicene and germanene and tunable electronic properties of such systems indicate the possible applications of nextgeneration nanobiosensor technology for biomolecules. The obtained results will provide reference to both theoretical and experimental studies.