Abstract

By using first principle calculations, a simple model of salinity sensor based on graphene electrode is constructed and its electron transport property is systematically investigated. It is found that all saltwater clusters at different salinity exhibit an obvious increase of the current while the saltwater to be detected is passing through the device. Moreover, only changing one Na atom acted as the conductive medium, and the electron transport behaviors could be clearly distinguished among the saltwater by negative differential resistance phenomenon, which demonstrates that the graphene-based salinity sensor could be capable of distinguishing saltwater at different salinity efficiently and accurately. This study provides a new path for the creation of the novel salinity sensor by graphene and other 2D material electrode.

1. Introduction

It has been attracted more attention since the stable, one-atom layer graphene was fabricated in 2004 by exfoliating from its van der Waals solid graphite [1]. From then on, the atomic-scale thickness two-dimensional van der Waals layered materials started to intrigue a research boom because of their extremely large surface area [24] and through chemical vapor deposition (CVD) and other process methods, many graphene-like 2D nanomaterials have been fabricated with the similar hexagonal structure, such as boron nitride (h-BN) and transition-metal dichalcogenides (TMDCs) [3, 5, 6]. Due to quantum confinement effect, they exhibit marvelous electrical, optical, mechanical, and chemical properties. Based on these unique electronic properties, they have broad application prospects in electronics, optoelectronics, energy, environment, and other fields [26]. As a result of its wide bandgap insulating property, h-BN can be widely applied in nonlinear optics and ultraviolet lasers/detectors [7]. TMDCs have very rich electrical properties and a direct bandgap, making them show great potential for applications in areas such as optoeletronics and nanoelectronics [8].

Although various novel graphene-like two-dimensional materials have inspired a new wave of research and been studied experimentally and theoretically [912], graphene is the most popular candidate for the next generation of nanodevices because of its extremely high electron mobility which is beyond 200,000 cm2 V−1s−1 at electron densities of about [13, 14] while many other 2D materials do not possess. Furthermore, compared with other two-dimensional materials, it possesses superior mechanical strength and flexibility [15]. Therefore, graphene is the most popular choice to develop the electrodes of the sensor owing to the large surface area and semimetallic property. Another advantage for sensing purpose is that graphene has less impact on the environment than other electrode materials [16]. Many of these properties lead the graphene to sensor applications, such as electrochemical and electrical sensors. In recent years, because of the higher electron mobility and surface activity, graphene has been largely considered as transistors for the detection of biomolecules and gas reducing response time and manufacture cost [17, 18]. The high sensitivity of graphene to the local environment shows great advantage in sensing applications, where ultralow concentrations of adsorbed molecules react significantly to the electronic properties of graphene [1922].

To date, detecting polluted gas and gaseous water molecules with graphene has been investigated extensively [2325]. Yokoyama et al. [23] found that at a higher voltage of 0.9 V, suspended graphene sensor displayed the good selectivity to H2, and Leenaerts et al. [24] explored the graphene sensor adsorbed with various kinds of gas molecules theoretically which is consistent with the previous experimental work performed by Schedin et al. [25]. Salinity is the most important physical parameter in oceanography for understanding global climate and the dynamic feature of marine ecosystems, where differences in salinity drive ocean circulation [26]. Therefore, accurate salinity measurement is of great importance. Although the salinity sensor has been widely applied, there are still some problems to be solved. Slow reaction times of conductivity sensors which are installed in the glider produce erroneous salinities, especially on strong thermoclines, so the high sensitivity and miniaturization of the sensor are essential for accurate salinity measurements [27]. As 2D material exhibiting high conductivity and sensitivity, graphene is the most proper candidate to achieve accuracy and higher speed in detection where salinity sensing is required [20]. But detecting the salinity of the saltwater with the electrode of graphene has not been studied systematically yet. In this paper, we design a sensor with graphene electrode and perform a systematic study of the transport properties of a small drop of water with different salinity by using the Extended Huckel method and the nonequilibrium Green function.

2. Materials and Methods

Graphene nanoribbons perform different electronic properties; for the armchair edges, it shows semiconducting, and for the zigzag edges, it shows the metallic. Therefore, the graphene nanoribbons with the zigzag edges () are chosen as the electrode of the salinity sensor, which is cut from the graphene. Many other researchers chose the field effect transistor- (FET-) based structure similar to ours as the sensor platform both in experiment and theory, and they found that using graphene as electrode could achieve accuracy and higher speed in detection [20, 21, 28]. All the dangling bonds at the edges are very unstable and tend to bond with other atoms in the vicinity, so carbon atoms at the edges are passivated by the hydrogen atoms to avoid the effects of dangling bonds. The structures are fully optimized to make sure that the fore on each atom is less than 0.5 eVÅ-1.

Figure 1 shows the detected saltwater passing through the simulated devices in which the h-BN nanowalls are fabricated to allow the saltwater pass through the channel effectively. The water cluster is with 33 water molecules to ensure the simulated density of 1 g/cm3. Recently, the researchers found that the minimal structure for proper solvation is made of 21 water molecules [29]. It is well known that almost 99% of the elements in seawater solutes are sodium (Na) and chloride (Cl). To simulate the saltwater simply, the Na atom is chosen for the solutes of simulated saltwater and the salinity of saltwater could be further deduced by calculating the Na concentration. The Na concentration is calculated by relative atomic mass in which Na atom is 23 and H2O is 18. Therefore, the concentration of one, two, and three Na atoms in seawater solutes is 6%, 11%, and 15%, respectively. The water cluster is simulated enough time to an equilibrium state by dynamic simulation, and the trajectory recorded is 0.1 ps, which is repeated three times. The simulated water molecules in the middle of axis are substituted by Na atoms to ensure the solutes be detected effectively by graphene electrode in the channel. Figure 1(a) shows the top view of the device. The device consists of three parts: the left electrode, the central scattering region, and the right electrode, in which the width of electrodes is and central region is The length of the device is 44.8 nm, and the gap width of channel is 5.9 nm which allows the saltwater to pass freely through. Calculations of the electron transport properties were carried out in the bipolar system using the Extended Huckel method and the nonequilibrium Green function (NEGF) in the AtomistixToolKit (ATK) calculation package [30, 31], which only considers the electron-electron and not the electron-phonon interaction. However, for graphene electrode, this method is enough to study the electron transport properties since it is consistent with the experimental work [32]. The set of Cerda parameters applies to carbon and hydrogen, where the vacuum is set to -7.36577 eV and -6.2568 eV, respectively. Density mesh cut-off is 100 Ry, and electron temperature is 300 K. The maximum interaction range is set to be 10 Å. The criterion of electron energy convergence is 10-5 eV. The current is calculated by integrating a transmission function using the Landauer formula.

3. Results and Discussion

Figure 2 shows the simulated current-voltage (I-V) characteristics of the pure water drop and saltwater with different Na concentration. It is noticed that when the bias is applied, there is no current in the transport direction for the pure water drop. However, all the saltwater clusters result in the onset of the current where the negative differential resistance (NDR) phenomenon is observed. As the pure water cluster is added with Na atom, the current initially increases with bias and then fluctuates slightly around the state of equilibrium between the current of 0.3 and 0.5 nA before 0.4 V. When the bias voltage is further increased and exceeds 0.4 V, the current breaks the stable state and begins to come up. It is clear that in all three cases, the current climbs significantly as the bias voltage reaches 0.6 V and the peak currents are observed. Before that the growth trend becomes divided among the three cases, for the saltwater with 15% Na, the current ascends directly to the peak, while for the other two cases, the current increases fast first and slowly afterwards, even more approaching a peak value. The maximum obtained value of current for the saltwater with 6%, 11%, and 15% Na is 1.65 nA, 1.96 nA, and 2.35 nA, respectively. Further increase in the bias voltage gives rise to a current drop until the bias of 0.7 V. Similar to the growth tendency, except the case with 15% Na, the others drop quickly to a low value where the saltwater with 6% Na (only one Na atom) results in the current maintaining the value of 0.3 nA. For the saltwater with 15% Na, the current declines less than half which still remains the high value among all the cases. In the rest time, the cases with 11% and 15% Na appear a tight range and then the current continues to grow. The minimum current obtained for the saltwater with 6%, 11%, and 15% Na is 0.34 nA, 0.59 nA, and 1.43 nA individually. Therefore, an obvious NDR behavior appears in all the cases, and the peak-to-valley ratio (PVR) is calculated to compare the difference. The obtained PVR is 4.9, 3.3, and 1.6 for the saltwater with 6%, 11%, and 15% Na, respectively. According to the above analysis, it is clear that the saltwater with different Na concentration shows obvious differences in I-V characteristics; especially from the pure water to seawater, the current undergoes a dramatic change, which demonstrates that the graphene electrode is extremely sensitive to the change of the salinity. As mentioned above, through the simple saltwater models, the salinity could be deduced and the salinity of the detected saltwater could match the I-V characteristics, that is, the measurement principle of conductivity sensor in the salinometer. The peak, valley current, and PVR are all could be the parameters of the conductivity sensor because there are significant differences in these obtained values for the saltwater with different Na concentration. Moreover, the conductivity parameters differ obviously when altering one solute atom, indicating that the sensor with graphene electrode is a salinity sensor with high accuracy. In the follow-up works, we will consider the electron-phonon interaction and study the thermal conduction of the 2D materials by machine-learning interatomic potentials [33, 34].

In order to gain further insight into the parameters which shows the difference in the I-V characteristics among the cases with Na atom, transmission spectrum at some bias voltages is plotted in Figure 3. The transmission spectrum is influenced by both energy and bias voltage which indicates the probability of an electron transporting through the channel in certain energy. When a bias voltage is applied, the current can be calculated by the Landauer formula [35]: where is the bias voltage, is the transmission probability, and and are Fermi-Dirac distribution function of the right and left electrodes.

According to the Landauer formula, the current is calculated by integrating the transmission probability in the bias window, indicating that the transmission spectrum going to the bias window mainly contributes to the current. To understand the NDR effect, the case with 15% Na is taken for an example, as shown in Figure 3(a). At the bias of 0.4 V, the transmission spectrum exhibits two obvious peaks at 0.26 eV and -0.24 eV. As the peaks are not in the bias window, they could not contribute to the current resulting in a low current. However, when the bias is increased to 0.6 eV, there is also a dominant peak at -0.24 eV and another smaller peak at 0.26 eV, but the two peaks locate inside the bias window leading to a peak current. At 0.7 V, the transmission probability still shows two peaks with the same energy where the peak value at -0.24 eV drops significantly and almost increases double at 0.26 eV compared to the bias of 0.6 eV, resulting a valley current which is higher than that of the bias of 0.4 eV. Next, peak current is taken for exhibiting the difference in I-V characteristics among the models with different Na concentration by transmission spectrum. As can be seen from Figure 3(b), they all exhibit a dominant peak at energy -0.24 eV and small peaks at other energies inside the bias window. Although the dominant peak grows a little with the increase of Na concentration, it contributes significantly to the peak current, especially the saltwater with the concentration of Na from 6% to 11%. At the energy 0.26 eV, the contribution of small peak to current is increased with the higher growth in transmission spectrum. It can be seen from the above analysis that the current peak is driven by the dominant transmission peak at -0.24 eV.

Figure 4 shows the transmission pathways of saltwater with 15% Na at 0.6 V at energy -0.24 eV. The transmission pathways describe the paths of electrons transporting through the device where the arrow represents the direction and the color indicates the relative magnitude of the transmission probability. The transmission pathways are stronger at the left hand of the device than those at the right side. The pathways also occur in the saltwater cluster between the graphene electrodes where the Na atoms provide the transport paths in the saltwater cluster. And this is consistent with the above analysis that only the pure water cluster could not generate current, and the current is enhanced obviously by adding Na atom.

In the end, we construct a sensor device with longer graphene electrode in which the size of channel and Na concentration is the same with the abovementioned device, as shown in Figure 5(a). I-V characteristics of the pure water and saltwater with different Na concentration are calculated (Figure 5(b)). It is also found that pure water could not generate the current, and there is NDR effect in all the saltwater clusters. Among all the cases of different Na concentration, it also results in significant differences in I-V characteristics. Therefore, no matter the graphene electrode is long or short, the graphene sensor is sensitive to the Na concentration.

4. Conclusions

Using the nonequilibrium Green function together with the Extended Huckel method, a simple model of salinity sensor based on graphene electrode is constructed and the electron transport property of the simulated saltwater is investigated. It is found that except the pure water, all the saltwater clusters exhibit an obvious increase of the current and the NDR effect while the detected saltwater is passing through the device, because the Na atom acts as a conductive medium in the simulated saltwater cluster. Moreover, only changing one Na atom acted as the conductive medium, and the I-V characteristics could be clearly distinguished among the saltwater by peak current, valley current, and the obtained PVR, no matter the graphene electrode is long or short. This study may provide a new path for the creation of the novel salinity sensor by graphene and other 2D materials electrode.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the support from the National Natural Science Foundation of China (11804179) and Natural Science Foundation of Shandong Province (ZR2019PA018).