Next Article in Journal
Editorial for the Special Issue on Graphene-Nanocomposite-Based Flexible Supercapacitors
Previous Article in Journal
An Adaptive Noise Reduction Method for High Temperature and Low Voltage Electromagnetic Detection Signals Based on SVMD Combined with ICEEMDAN
Previous Article in Special Issue
Density Functional Theory Provides Insights into β-SnSe Monolayers as a Highly Sensitive and Recoverable Ozone Sensing Material
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First-Principles Insights into Highly Sensitive and Reusable MoS2 Monolayers for Heavy Metal Detection

1
Department of Engineering Technology, Guangdong Open University, Guangzhou 510091, China
2
Centre for Optical and Electromagnetic Research, South China Academy of Advanced Optoelectronics, South China Normal University, Guangzhou 510006, China
3
Ministry of Education Key Laboratory of Textile Fiber Products, School of Materials Science and Engineering, Wuhan Textile University, Wuhan 430220, China
4
School of Materials and Chemistry Engineering, Tongren University, Tongren 554300, China
5
School of Artificial Intelligence, Guangdong Vocational College of Post and Telecom, Guangzhou 510630, China
6
Big Data Engineering College, Kaili University, Kaili 556011, China
*
Authors to whom correspondence should be addressed.
Micromachines 2024, 15(8), 978; https://doi.org/10.3390/mi15080978
Submission received: 24 June 2024 / Revised: 23 July 2024 / Accepted: 28 July 2024 / Published: 30 July 2024
(This article belongs to the Special Issue Gas Sensors: From Fundamental Research to Applications)

Abstract

:
This study explores the potential of MoS2 monolayers as heavy metal sensors for As, Cd, Hg, and Pb using density functional theory (DFT) and Non-Equilibrium Green’s Function (NEGF) simulations. Our findings reveal that As and Pb adsorption significantly alters the surface structure and electronic properties of MoS2, introducing impurity levels and reducing the band gap. Conversely, Cd and Hg exhibit weaker interactions with the MoS2 surface. The MoS2 monolayer sensors demonstrate exceptional sensitivity for all four target heavy metals, with values reaching 126,452.28% for As, 1862.67% for Cd, 427.71% for Hg, and 83,438.90% for Pb. Additionally, the sensors demonstrate selectivity for As and Pb through distinct response peaks at specific bias voltages. As and Pb adsorption also induces magnetism in the MoS2 system, potentially enabling magnetic sensing applications. The MoS2 monolayer’s moderate adsorption energy facilitates rapid sensor recovery at room temperature for As, Hg, and Cd. Notably, Pb recovery time can be significantly reduced at elevated temperatures, highlighting the reusability of the sensor. These results underscore the potential of MoS2 monolayers as highly sensitive, selective, and regenerable sensors for real-time heavy metal detection.

1. Introduction

Heavy metals such as arsenic (As), cadmium (Cd), mercury (Hg), and lead (Pb) are known for their acute toxicity and carcinogenic properties, even at minute concentrations. Furthermore, they resist biodegradation in the environment. Found in water bodies, air, and soil, these pervasive environmental contaminants enter the human body through various pathways, including inhalation, the food chain, and absorption through the skin and gastrointestinal tract. Upon entry, heavy metals react with organic molecules in the body to form metal complexes, altering the original biochemical functions of substances and thus posing significant health risks.
In the environment, heavy metals exist in various forms, including ions, particle-bound (HMp), and elemental (HM0). Traditional detection methods have primarily focused on the ionic and HMp states. However, the elemental form, HM0, presents a challenge due to its low solubility in water or electrolyte solutions, rendering conventional detection methods ineffective [1]. For instance, elemental mercury (Hg0), emitted from coal-fired boiler flue gases, and gaseous cadmium (Cd0), from solid waste incineration, are difficult to measure and capture due to their low solubility and surface affinity, respectively [2,3]. Similarly, elemental lead (Pb0), released into the fuel flue gases with low solubility in water, poses a challenge for purify using wet flue gas desulfurization systems [4]. Elemental arsenic (As0) is also one of the predominant forms in which arsenic exists [5]. Therefore, there is an urgent need to develop sensors capable of real-time detection of HM0.
Recent studies have highlighted the potential of two-dimensional materials in heavy metal adsorption and sensing. M. Ghashghaee et al. discovered that black phosphorene monolayers could adsorb Cd, Pb, Hg, and As, effectively removing these heavy metals from the environment [6]. M. A. Hamed et al. reported that graphene-structured ZnO exhibited adsorption capabilities for Cd [7]. Carbon-based two-dimensional materials have also been employed for heavy metal adsorption. C. Zhao et al. demonstrated that pyridinic nitrogen-doped graphene modified with Pd clusters could effectively adsorb Hg and As [8]. I. Shtepliuk et al. utilized the selectivity of graphene/SiC Schottky diodes for highly sensitive detection of lead atoms [9]. J. Liao et al. observed that single vacancy defects significantly enhanced the adsorption of Cd on graphene [10]. Research by M. Srivastava et al. revealed that boron-modified graphene and copper-modified nitrogen-doped graphene could detect As in water [11]. Furthermore, J. Zou et al. modified glassy carbon electrodes with sulfur-doped graphitic carbon nitride nanosheets, which proved to be sensitive and selective for Pb detection [12].
Among two-dimensional chalcogenides, molybdenum disulfide (MoS2) monolayers have attracted attention for their high carrier mobility, tunable bandgap, and extensive specific surface area, positioning them as promising candidates for gas sensing applications [13,14,15,16,17,18,19,20,21]. X. Mu and Zhao et al. investigated the adsorption performance of 1T-MoS2 and 2H-MoS2 for elemental mercury [22,23]. However, reports on the sensing capabilities of MoS2 monolayers for heavy metals are still forthcoming.
To address this knowledge gap, our study employs density functional theory (DFT) to explore the heavy metal sensing capabilities of MoS2-based monolayer sensors. The adsorption energies, charge transfer, density of states (DOS), band structure, and Bader analysis of As, Cd, Hg, and Pb on MoS2 monolayers have been investigated. The desorption process is predicted to be governed by the recovery time. This paper also measures the current voltage characteristics and sensitivity of the MoS2 monolayer when detecting heavy metals. Our findings could pave the way for the development of novel nanosensors for environmental monitoring and facilitate the rapid, selective detection of heavy metal contaminants.

2. Materials and Methods

Calculations were conducted using the Vienna Ab initio Simulation Package (VASP) 5.4, which implements spin-polarized density functional theory (DFT) [24,25]. The Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA) was employed for the exchange-correlation functional [26]. To account for van der Waals interactions within the MoS2 monolayer, the DFT-D3 method was utilized [27]. Geometry optimization continued until the forces on all atoms converged below a threshold of 0.02 eV/Å. The total energy convergence criterion was set to 10−5 eV. A plane wave energy cutoff of 400 eV was used, and a vacuum layer of 20 Å was applied perpendicular to the 2H-MoS2 monolayer to minimize interlayer interactions. Brillouin zone integration was performed using a 4 × 4 × 1 k-point Monkhorst–Pack mesh.
The adsorption energy Ead was computed using the relation:
E a d = E t o t a l E M o S 2 E m e t a l
where Etotal is the total energy of the fully relaxed adsorbed system consisting of the MoS2 monolayer and the adsorbed metal atom, E M o S 2 is the energy of the pristine MoS2 monolayer, and Emetal is the energy of the isolated heavy metal atom.
Additionally, quantum transport calculations were conducted using the TRANSIESTA code within the SIESTA 4.1.5 framework, utilizing the non-equilibrium Green’s function method [28]. The Landauer–Büttiker formula was applied to calculate the current-voltage (I-V) characteristics of the MoS2 monolayer both before and following heavy metal adsorption under [29]:
I V b = 2 e h μ R μ L T E , V b f E μ L f E μ R
where e and h are the electron charge and Planck’s constant, respectively. I is the current through the central scattering region. μ and f E μ are the electrochemical potentials and the Fermi–Dirac distribution function, where the subscript L/R is employed to differentiate between the left and right electrodes, respectively. By employing Green’s function formalism, the transmission coefficient T E , V b at energy E and bias voltage Vb was calculated.
For the visualization of crystal structures and differential charge density, the software VESTA 3 was employed [30]. Data post-processing was performed using the tool vaspkit 1.3.3 [31].

3. Results and Discussion

3.1. Structures

As depicted in Figure 1, the MoS2 monolayer supercell exhibits a distinct structural arrangement, where each Mo atom is coordinated with six S atoms, and each S atom is coordinated with three Mo atoms. The Mo–S bond length is measured at 2.417 Å, aligning remarkably with previous studies, which reported bond lengths of 2.424 Å [19]. The bond angles within the MoS2 lattice, specifically the S-Mo-S and Mo-S-Mo angles, are 82.60°, forming a hexagonal configuration.
In this study, the energies of heavy metal atoms at four distinct sites on the MoS2 surface have been computed: the hexagonal center (H site), atop the Mo atom (TMo site), atop the S atom (TS site), and the bridge site of the Mo-S bond (B site). These calculations aim to ascertain the most energetically favorable adsorption sites for heavy metals on the MoS2 monolayer surface.
Figure 2 presents the most stable adsorption configurations of heavy metals—As, Pb, Cd, and Hg—on the MoS2 monolayer. Geometric optimization results indicate a predilection for As and Pb to adsorb at the TS site, where they effectively engage with the MoS2 surface, achieving the lowest adsorption energies. In contrast, Cd and Hg show a preference for adsorbing on the TMo site.
During the adsorption process, the heavy metals cause alterations in the underlying bond angles of the MoS2 structure. Specifically, following the adsorption of Hg and Cd, the S-Mo-S bond angles increase to 83.15° and 83.36°, respectively. Adsorption of As and Pb results in an increase in the Mo-S-Mo bond angles to 83.52° and 83.74°, respectively. These changes in bond angles indicate the extent of structural perturbation imposed by different heavy metals on the MoS2 surface, with Pb exerting the most significant impact, followed by As, Cd, and Hg, in that order. The adsorption distance data show that As and Pb are closer to the MoS2 surface, at 2.175 Å and 2.625 Å, respectively, compared to the greater distances for Cd (3.358 Å) and Hg (3.526 Å). This suggests a tighter interaction and stronger interfacial forces between As and Pb with the MoS2 surface, leading to more pronounced structural distortions. Conversely, the longer adsorption distances for Cd and Hg suggest relatively weaker interactions and a diminished interfacial coupling effect.
A detailed analysis of the adsorption characteristics of different heavy metal ions on the MoS2 monolayer surface has been performed, and the key parameters, including adsorption energy, adsorption distance, Bader charge transfer, and the corresponding changes in the material’s bandgap, are summarized in Table 1.
Among the listed data, the heavy metal elements Cd and Hg exhibit longer adsorption distances, lower corresponding adsorption energies, and smaller charge transfers. This indicates relatively weaker interactions between Cd and Hg with the MoS2 monolayer, which are not conducive to forming stable chemisorption states. Conversely, As and Pb demonstrate significantly enhanced adsorption characteristics on the MoS2 monolayer. Not only is their adsorption energy significantly reduced, indicating more stable adsorption, but their adsorption distances are also substantially shorter. This implies that the interactions between As and Pb with the MoS2 monolayer are more intense, leading to more efficient surface adsorption reactions.

3.2. Eletronic Properties

In the exploration of monolayer MoS2 as a sensing material, the analysis of its band structure is crucial, particularly the characteristics of the bands near the Fermi level (from −2 eV to 2 eV) which directly influence the material’s conductivity changes and sensing response sensitivity. The theoretical band gap of pristine monolayer MoS2 has been calculated to be 1.585 eV.
Upon adsorption of As atoms on the MoS2 monolayer surface, a significant alteration in the band structure is observed. As depicted in Figure 3, four primary impurity levels originating from the As atoms emerge between the valence band maximum and the conduction band minimum. The introduction of these impurity levels substantially narrows the original band gap, reducing the new band gap of MoS2 post-As adsorption to 0.242 eV. This represents an approximate 84.7% decrease compared to the baseline state, suggesting a marked enhancement in the electrical conductivity of the MoS2-As adsorption system, thereby improving its efficiency as a heavy metal sensor.
On the other hand, the adsorption of Cd or Pb on the MoS2 monolayer introduces a new impurity level below the Fermi energy, resulting in band gaps of the MoS2-Cd and MoS2-Pb adsorption systems being reduced to 1.14 eV and 0.326 eV, respectively. Compared to the baseline state, the band gaps have been reduced by 28.1% and 79.4%, respectively. While adsorption in both cases promotes an increase in conductivity, the extent is not as significant as with As adsorption. In contrast, Hg atom adsorption on the MoS2 monolayer does not induce noticeable impurity level formation, implying that post-Hg adsorption, the band gap of the MoS2 monolayer decreases by only 3.2%. This suggests that Hg adsorption has a minimal impact on the change in conductivity of the MoS2 monolayer, and compared to the adsorption of other elements, it does not significantly alter its conductive response characteristics as a sensor.
To further elucidate the adsorption characteristics of MoS2 towards heavy metals, the impact of heavy metal adsorption on the electronic properties of MoS2 was investigated by analyzing the total density of states (TDOS) and partial density of states (PDOS) as shown in Figure 4.
After As adsorption on MoS2, four new hybridized energy levels were generated near the Fermi level. These levels correspond to spin-up peaks at −0.809 eV and −0.292 eV, a spin-down peak at −0.070 eV, and a sharp, narrow peak just above the Fermi level at 0.152 eV. The p-orbitals of As play a dominant role in the formation of these hybridized levels, simultaneously inducing peak states in the density of states for Mo’s d-orbitals and S’s p-orbitals at corresponding energy positions. These results demonstrate strong orbital hybridization between As and the surface atoms of MoS2, and due to the asymmetry of the DOS curves for spin-up and spin-down, it is evident that As adsorption imparts magnetism to the MoS2 system, with a calculated magnetic moment of 1 µB.
For Cd adsorption on MoS2, a DOS peak is observed at −0.121 eV, resulting from the hybridization of Cd’s s-orbital electrons with Mo’s d-orbital electrons and S’s p-orbital electrons. However, Cd adsorption does not lead to significant spin polarization, as the DOS curves remain symmetrical for spin-up and spin-down, indicating that the Cd adsorption system is non-magnetic.
After Pb adsorption on MoS2, a DOS peak is observed at −0.121 eV, primarily attributed to the cooperative effect of Pb’s p-orbital electrons with Mo’s d-orbital electrons and S’s p-orbital electrons. Similar to As adsorption, the system post-Pb adsorption exhibits asymmetry in the DOS curves for spin-up and spin-down, resulting in a magnetic moment of 2 µB, confirming the magnetic nature of the Pb adsorption system.
As a result, As, Cd, and Pb all exhibit a certain degree of orbital hybridization effects during their interaction with MoS2. The adsorption of As and Pb leads to a transition of the MoS2 system to a magnetic state, suggesting that MoS2 holds promise as a magnetic sensing material for As and Pb detection, while the detection of Cd and Hg does not depend on a magnetic response mechanism.
To gain deeper insight into the bonding characteristics between the adsorbed metal atoms and the MoS2 monolayer surface, we analyzed the differential charge density. This quantity is calculated using Equation (3):
ρ = ρ t o t a l ρ M o S 2 ρ m e t a l
where ρ represents the charge density difference, ρ t o t a l is the charge density of the combined adsorption system, ρ M o S 2 is the charge density of the pristine MoS2 monolayer, and ρ m e t a l is the charge density of the isolated adsorbate atom. Figure 5 illustrates the differential charge density, with regions of electron accumulation depicted in yellow and areas of electron depletion shown in blue.
Through the study of As and Pb adsorption on the MoS2 surface, significant regions of charge exchange were identified within these adsorption systems, as further confirmed by Bader charge analysis depicted in Table 1. Notably, As and Pb exhibited substantial charge transfers to the MoS2, with values of −0.243 e and −0.410 e, respectively. In contrast, Cd exhibits a relatively weaker charge exchange with the MoS2 surface, with a value of −0.115 e. The charge transfer between Hg and MoS2 is even more negligible, amounting to only −0.023 e.

3.3. Recovery Time

The recovery time is a crucial property for sensor performance as it determines the speed and repeatability of the sensor’s response. The recovery time (τ) represents the time needed for the sensor surface to restore its initial adsorption capacity after becoming saturated with adsorbates. The recovery time, denoted by (τ), is deduced from Equation (4), which is an extrapolation from classical transition state theory references:
τ = exp E a d k B T / ω
where ω signifies the attempt frequency, postulated to be 1013 s−1 [32]. Ead represents the adsorption energy, kB is the Boltzmann constant, and T denotes the temperature. After adsorption of various heavy metals on MoS2, the calculated recovery times within the 298–498 K range are presented in Table S1. As delineated in Figure 6, the recovery times for heavy metals on the MoS2 monolayer are inversely proportional to the temperature increment. At 298 K, the desorption times for As, cadmium (Cd), and Hg from the MoS2 monolayer surface are 9.67 s, 3.54 × 10−6 s, and 1.56 × 10−7 s, respectively, indicating a swift recovery performance at ambient conditions. For Pb, the recovery time on a MoS2 monolayer reduces to 2.21 × 1010 s at 298 K, and further to 9.32 s at an elevated temperature of 498 K. This observation suggests that the MoS2 monolayer, when employed as a Pb sensor, exhibits reusability at the heightened temperature of 498 K. Consequently, the MoS2 monolayer emerges as a promising candidate for sensing As, Cd, Hg, and Pb, attributable to its expedited recovery dynamics.

3.4. Sensitivity

The electrical transport properties of the MoS2 monolayer upon heavy metal adsorption were investigated by constructing a double-probe model, comprising two electrodes with an intervening scattering region, as depicted in Figure 7. A bias voltage was applied across the electrodes. This quantum transport problem was solved using the NEGF-DFT method under open boundary conditions and in the presence of external electric fields.
Figure 7 also delineates the current-voltage (I-V) characteristics of MoS2 monolayers, both in pristine form and following the adsorption of heavy metals. At a bias voltage of 0.2 V, the MoS2 monolayer with adsorbed As demonstrated a significantly enhanced current, reaching 1.75 × 10−8 µA, which is markedly higher than the basal current of intrinsic MoS2 at 1.39 × 10−11 µA. Comparatively, the currents observed for MoS2 monolayers with adsorbed Cd, Hg, and Pb were 8.92 × 10−12 µA, 9.01 × 10−12 µA, and 1.35 × 10−10 µA, respectively, indicating a pronounced selectivity for arsenic.
For Pb adsorption, the current at low bias voltages (≤0.6 V) was markedly higher than that of intrinsic MoS2, while at higher bias voltages (≥0.8 V), the disparity was considerably reduced. Notably, at a bias voltage of 0.4 V, the MoS2 monolayer with adsorbed arsenic (As) manifested a current of 2.54 × 10−8 µA, which is considerably higher than the currents measured for MoS2 monolayers with adsorbed As (3.07 × 10−10 µA), Cd (1.95 × 10−11 µA), Hg (1.72 × 10−11 µA), and the intrinsic MoS2 (3.04 × 10−11 µA).
The MoS2 monolayer with Hg adsorption showed a relatively minor difference in current from intrinsic MoS2, except in the 0.5–0.6 V range. The current for Cd-adsorbed MoS2 was generally lower than that for As, but higher than for Hg.
The sensitivity (S) of MoS2 monolayers was also evaluated using Equation (5):
S = I I 0 / I 0
where I represents the current measured across the MoS2 monolayer upon heavy metal adsorption, and I0 denotes the intrinsic current in the absence of heavy metal contaminants.
As illustrated in Figure 8, the sensitivity of the MoS2-based monolayer sensors can be quantitatively estimated. Upon adsorption of As, the MoS2 sensor exhibits a remarkable sensitivity, reaching 126,452.28% at a bias voltage of 0.2 V. This sensitivity is significantly higher compared to other heavy metals, indicating a pronounced selectivity for As. Moreover, the MoS2 sensor maintains a sensitivity above 200% across the entire voltage range below 1.2 V after As adsorption, demonstrating exceptional sensitivity and superior sensing performance.
In the case of Pb, the MoS2 sensor achieves a sensitivity of 83,438.90% at 0.4 V, which is considerably higher than that for other heavy metals, suggesting good selectivity at this voltage. Within a voltage range of ≤0.7 V, the sensitivity for Pb remains above 200%.
For Cd, the peak sensitivity of 1862.67% is reached at 0.8 V, with the sensitivity exceeding 200% at bias voltages of 0.1 V, 0.3 V, and 0.6–0.8 V. Conversely, the sensitivity for Hg is relatively lower, peaking at 427.71% at 0.5 V and surpassing 200% within the range of 0.5–0.6 V.
The limit of detection (LOD) signifies the lowest analyte concentration a sensor can reliably detect. Determined via signal-to-noise ratio (SNR), LOD compares the analyte signal to a blank signal to establish the minimum detectable concentration. For MoS2 sensors, current serves as the output signal. Baseline current (I0) represents the MoS2 current without heavy metal adsorption, while current (I) reflects the post-adsorption state.
S = I I 0 / I 0 = I / / I 0 1 = SNR 1
Conventionally, LOD corresponds to an SNR of 3:1, which also translates to a sensitivity of 200%.
To quantify heavy metal adsorption on MoS2 monolayers, concentration is expressed as the atomic count per area. The calculations above employed a 25.54 Å × 22.12 Å scattering region (564.94 Å2). A single heavy metal atom within this area (concentration: 1/564.94 Å2) yields MoS2 sensitivities surpassing the LOD for As (126,452.28%), Cd (1862.67%), Hg (427.71%), and Pb (83,438.90%).
Doubling the scattering region while maintaining a single metal atom reduces concentration to 1/1129.88 Å2. Under these conditions, the sensitivities of MoS2 towards As, Cd, Hg, and Pb are 132.33%, 215.42%, 122.38%, and 828.19%, respectively. Here, MoS2 sensitivity significantly exceeds the LOD for Pb, marginally surpasses it for Cd, but falls below for As and Hg. Consequently, MoS2 detects Pb at concentrations below 1/1129.88 Å2, Cd near 1/1129.88 Å2, and As and Hg between 1/1129.88 Å2 and 1/564.94 Å2.
These findings underscore the high sensitivity and selectivity of MoS2-based monolayer sensors towards heavy metals, positioning them as promising candidates for future environmental monitoring applications.
A comparative investigation of different two-dimensional nanomaterials’ sensitivity to heavy metal detection [7,8,9,10,11,22,33,34,35,36,37,38,39,40] has been conducted in this study. The data are clearly listed in Table S2 and are presented in graphical form in Figure 9. As depicted in Figure 9a, the MoS2 monolayer exhibits a sensitivity of 126,452.28% to As. Although this is slightly less than that of zigzag ZnO nanoribbons (3.52 × 107%) [33] and SnSe (384,467%) [34], it is significantly superior to SnS (43,286%) [35], Cu-N-doped graphene-H2O (55%) [36], graphene (36.06%) [37], Cu-B-doped graphene-H2O (25%) [36], and graphene-H2O (19%) [36]. Additionally, MoS2 demonstrates a shorter recovery time under ambient conditions, facilitating rapid and continuous detection. In contrast, zigzag ZnO nanoribbons are lacking in adsorption energy data, which hinders the assessment of their reusability at room temperature.
As shown in Figure 9b, upon further investigation of Cd sensing, MoS2 demonstrates the highest sensitivity reported to date, reaching 1862.67%, exceeding SnSe (1462%) [34] and SnS (442%) [35], and far outperforming the graphene/SiC composite material (8.5702 × 10−8%) [9]. For Hg detection, although the sensitivity of MoS2 is 427.71%, lower than SnSe (1791%) [34], it is still significantly higher than SnS (10%) [35] and graphene/SiC (3.868 × 10−8%) [9], as Figure 9c illustrates.
In the context of Pb sensing depicted by Figure 9d, the sensitivity of MoS2 reaches 83,438.90%, which is second only to zigzag ZnO nanoribbons (2.86 × 107%) [33] and SnS (368,860%) [35], but far exceeds SnSe (26,160%) [34], graphene/SiC (13,600%) [9], N-doped graphene-H2O (329.04%), and B-doped graphene-H2O (16.12%) [38]. Notably, the recovery time of MoS2 at 498 K is only 9.32 s, compared to SnS, which requires 76.44 s at 598 K, and SnSe, which takes 38.3 s at 498 K to recover fully. This underscores the rapid response capability of MoS2.
A critical advantage of MoS2 monolayer sensors over other 2D material-based heavy metal sensors lies in their ability to selectively detect specific heavy metals in complex samples. Unlike other 2D materials, such as tin disulfide (SnS) and tin selenide (SnSe), where multiple heavy metals exhibit similar response peaks, MoS2 displays unique response peaks for each heavy metal.
As illustrated in Table S2, SnS exhibits response peaks for Cd, Hg, and Pb at the same voltage (0.5 V), leading to potential interference issues. In samples containing multiple heavy metals, the signal of more sensitive metals like Pb, Cd, and Hg can be masked by the stronger response of Pb. Similarly, SnSe exhibits response peaks for As, Cd, Hg, and Pb at distinct bias voltages (0.2 V, 0.1 V, 0.1 V, and 0.2 V, respectively). However, in mixed heavy metal samples, Cd and Hg, as well as As and Pb, can still interfere with each other.
In contrast, MoS2 demonstrates remarkable selectivity, with distinct response peaks for different heavy metals at varying bias voltages. For instance, As, Cd, Hg, and Pb exhibit response peaks at 0.2 V, 0.8 V, 0.5 V, and 0.4 V, respectively. This unique characteristic allows for the selective detection of specific heavy metals by adjusting the working voltage applied across the MoS2 monolayer.

4. Conclusions

The sensing performance of the MoS2 monolayer on heavy metals (As, Cd, Hg, and Pb) is investigated using DFT and NEGF simulations. As and Pb adsorption onto MoS2 monolayers causes significant surface structural modifications and intensive electrical interactions. These interactions introduce impurity energy levels within the adsorption system, leading to significant band gap alterations. Conversely, the MoS2 surface is less affected by Cd and Hg, as evidenced by a weaker electronic exchange. The sensitivities of MoS2 monolayers to As, Cd, Hg, and Pb are 126,452.28%, 1862.67%, 427.71%, and 83,438.90%, respectively. For As or Pb, the MoS2 monolayer sensors exhibit response peaks at different bias voltages that are significantly higher than those for other heavy metals, demonstrating a high degree of selectivity. When As or Pb are adsorbed on MoS2 monolayers, they acquire a non-zero magnetic moment of 1 µB or 2 µB, respectively, indicating that they may find use as magnetic sensors. The MoS2 monolayer’s moderate adsorption energy facilitates rapid sensor recovery for As, Hg, and Cd at room temperature. Furthermore, the recovery time for Pb can be remarkably reduced to just 9.32 s at 498 K, highlighting the MoS2 monolayer’s efficacy as a reusable sensor. The results of this study show that MoS2 monolayer has the potential to be a high sensitivity, selectivity, and reusability heavy metal sensing material that can be used to make heavy metal detectors that operate in real-time.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/mi15080978/s1, Table S1: Temperature-dependent recovery time of MoS2 for heavy metal adsorption; Table S2: Comparison between this work and previous reports.

Author Contributions

Conceptualization, J.W. (Jiayin Wu); methodology, J.W. (Jiayin Wu); software, Z.L.; investigation, J.W. (Jiayin Wu); resources, T.L.; writing—original draft preparation, J.W. (Jiayin Wu); writing—review and editing, Q.M., T.L., J.W. (Jingting Wei) and B.L.; supervision, Z.L. and X.X.; funding acquisition, Q.M., T.L., Z.L. and B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Foundation Research Project of Kaili University, grant number 2022ZD05; National Natural Science Foundation of China, grant number 22269019; the Key Research Platform and Research Project of General Colleges and Universities in Guangdong Province, grant number 2022KTSCX289, 2023ZDZX1069; 2019 Annual School-level Research Project of Guangdong Open University (Guangdong Polytechnic Institute), grant number ZD1905.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to thank Peiyuan Lai, Xianxu Li, Hongfei Qiu, Zhiqing Huang, and Pengchen Jian for their guidance on computational facilities.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kamalinahad, S.; Solimannejad, M.; Shakerzadeh, E. Sensing of Ozone (O3) Molecule via Pristine Singe-Walled Aluminum Nitride Nanotube: A DFT Study. Superlattices Microstruct. 2016, 89, 390–397. [Google Scholar] [CrossRef]
  2. Munsif, S.; Ayub, K.; Nur-e-Alam, M.; Ahmed, S.; Ahmad, A.; Ul-Haq, Z. Sensing of H2S, NO2, SO2, and O3 through Pristine and Ni-Doped Zn12O12 Nanocage. Comput. Theor. Chem. 2023, 1229, 114305. [Google Scholar] [CrossRef]
  3. Rad, A.S.; Ayub, K. O3 and SO2 Sensing Concept on Extended Surface of B12N12 Nanocages Modified by Nickel Decoration: A Comprehensive DFT Study. Solid State Sci. 2017, 69, 22–30. [Google Scholar] [CrossRef]
  4. Chen, Z.; Liu, H.; Wang, H.; Liu, Y.; Wei, Z. Flue Gas Pb0 Removal from Sludge Incineration through Biological Lead Oxidation Coupled Denitrification. Fuel 2024, 355, 129500. [Google Scholar] [CrossRef]
  5. Dhuldhaj, U.P.; Yadav, I.C.; Singh, S.; Sharma, N.K. Microbial Interactions in the Arsenic Cycle: Adoptive Strategies and Applications in Environmental Management. In Reviews of Environmental Contamination and Toxicology Volume 224; Whitacre, D.M., Ed.; Springer: New York, NY, USA, 2013; pp. 1–38. ISBN 978-1-4614-5882-1. [Google Scholar]
  6. Ghashghaee, M.; Ghambarian, M. Adsorption of Toxic Mercury, Lead, Cadmium, and Arsenic Ions on Black Phosphorous Nanosheet: First-Principles Calculations. Struct. Chem. 2019, 30, 85–96. [Google Scholar] [CrossRef]
  7. Hamed Mashhadzadeh, A.; Fathalian, M.; Ghorbanzadeh Ahangari, M.; Shahavi, M.H. DFT Study of Ni, Cu, Cd and Ag Heavy Metal Atom Adsorption onto the Surface of the Zinc-Oxide Nanotube and Zinc-Oxide Graphene-like Structure. Mater. Chem. Phys. 2018, 220, 366–373. [Google Scholar] [CrossRef]
  8. Zhao, C.; Wu, H. Density Functional Investigation of Mercury and Arsenic Adsorption on Nitrogen Doped Graphene Decorated with Palladium Clusters: A Promising Heavy Metal Sensing Material in Farmland. Appl. Surf. Sci. 2017, 399, 55–66. [Google Scholar] [CrossRef]
  9. Shtepliuk, I.; Eriksson, J.; Khranovskyy, V.; Iakimov, T.; Lloyd Spetz, A.; Yakimova, R. Monolayer Graphene/SiC Schottky Barrier Diodes with Improved Barrier Height Uniformity as a Sensing Platform for the Detection of Heavy Metals. Beilstein J. Nanotechnol. 2016, 7, 1800–1814. [Google Scholar] [CrossRef]
  10. Liao, J.; Yang, L.; Wang, C.Z.; Lin, S. Tuning Cd Adsorption Behaviours on Graphene by Introducing Defects: A First-Principles Study. Mater. Technol. 2017, 32, 840–844. [Google Scholar] [CrossRef]
  11. Srivastava, M.; Srivastava, A. First Principle Investigation of Boron Functionalized Graphene to Detect the Presence of “Arsenic” in Water. Mater. Today Proc. 2022, 48, 661–665. [Google Scholar] [CrossRef]
  12. Zou, J.; Mao, D.; Li, N.; Jiang, J. Reliable and Selective Lead-Ion Sensor of Sulfur-Doped Graphitic Carbon Nitride Nanoflakes. Appl. Surf. Sci. 2020, 506, 144672. [Google Scholar] [CrossRef]
  13. Cui, H.; Zhang, X.; Zhang, G.; Tang, J. Pd-Doped MoS2 Monolayer: A Promising Candidate for DGA in Transformer Oil Based on DFT Method. Appl. Surf. Sci. 2019, 470, 1035–1042. [Google Scholar] [CrossRef]
  14. Jasmine, J.M.; Aadhityan, A.; Preferencial kala, C.; Thiruvadigal, D.J. A First-Principles Study of Cl2, PH3, AsH3, BBr3 and SF4 Gas Adsorption on MoS2 Monolayer with S and Mo Vacancy. Appl. Surf. Sci. 2019, 489, 841–848. [Google Scholar] [CrossRef]
  15. Li, B.; Zhou, Q.; Peng, R.; Liao, Y.; Zeng, W. Adsorption of SF6 Decomposition Gases (H2S, SO2, SOF2 and SO2F2) on Sc-Doped MoS2 Surface: A DFT Study. Appl. Surf. Sci. 2021, 549, 149271. [Google Scholar] [CrossRef]
  16. Wang, J.; Zhou, Q.; Lu, Z.; Wei, Z.; Zeng, W. Gas Sensing Performances and Mechanism at Atomic Level of Au-MoS2 Microspheres. Appl. Surf. Sci. 2019, 490, 124–136. [Google Scholar] [CrossRef]
  17. Chakraborty, B.; Maity, I.; Chung, P.; Ho, M.; Bhattacharyya, P. Understanding the Highly Selective Methanol Sensing Mechanism of Electrodeposited Pristine MoS2 Using First Principle Analysis. IEEE Sens. J. 2021, 21, 16484–16491. [Google Scholar] [CrossRef]
  18. Gui, Y.; Liu, D.; Li, X.; Tang, C.; Zhou, Q. DFT-Based Study on H2S and SOF2 Adsorption on Si-MoS2 Monolayer. Results Phys. 2019, 13, 102225. [Google Scholar] [CrossRef]
  19. Qian, H.; Lu, W.; Wei, X.; Chen, W.; Deng, J. H2S and SO2 Adsorption on Pt-MoS2 Adsorbent for Partial Discharge Elimination: A DFT Study. Results Phys. 2019, 12, 107–112. [Google Scholar] [CrossRef]
  20. Jiang, T.; He, Q.; Bi, M.; Chen, X.; Sun, H.; Tao, L. First-Principles Calculations of Adsorption Sensitivity of Au-Doped MoS2 Gas Sensor to Main Characteristic Gases in Oil. J Mater Sci 2021, 56, 13673–13683. [Google Scholar] [CrossRef]
  21. Liu, H.; Wang, F.; Hu, K.; Li, T.; Yan, Y.; Li, J. The Adsorption and Sensing Performances of Ir-Modified MoS2 Monolayer toward SF6 Decomposition Products: A DFT Study. Nanomaterials 2021, 11, 100. [Google Scholar] [CrossRef]
  22. Mu, X.; Gao, X.; Zhao, H.; George, M.; Wu, T. Density Functional Theory Study of the Adsorption of Elemental Mercury on a 1T-MoS2 Monolayer. J. Zhejiang Univ. Sci. A 2018, 19, 60–67. [Google Scholar] [CrossRef]
  23. Zhao, H.; Mu, X.; Yang, G.; George, M.; Cao, P.; Fanady, B.; Rong, S.; Gao, X.; Wu, T. Graphene-like MoS2 Containing Adsorbents for Hg0 Capture at Coal-Fired Power Plants. Appl. Energy 2017, 207, 254–264. [Google Scholar] [CrossRef]
  24. Kresse, G.G.; Furthmüller, J. Furthmüller Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B Condens. Matter 1996, 54, 11169. [Google Scholar] [CrossRef]
  25. Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  26. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865, Erratum Phys. Rev. Lett. 1997, 78, 1396. [Google Scholar] [CrossRef]
  27. Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2010, 27, 1787–1799. [Google Scholar]
  28. Soler, J.M.; Artacho, E.; Gale, J.D.; García, A.; Junquera, J.; Ordejón, P.; Sánchez-Portal, D. The SIESTA Method for Ab Initio Order-N Materials Simulation. J. Phys. Condens. Matter 2002, 14, 2745. [Google Scholar] [CrossRef]
  29. Meir, Y.; Wingreen, N.S. Landauer Formula for the Current through an Interacting Electron Region. Phys. Rev. Lett. 1992, 68, 2512–2515. [Google Scholar] [CrossRef]
  30. Momma, K.; Izumi, F. VESTA3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  31. Wang, V.; Xu, N.; Liu, J.-C.; Tang, G.; Geng, W.-T. VASPKIT: A User-Friendly Interface Facilitating High-Throughput Computing and Analysis Using VASP Code. Comput. Phys. Commun. 2021, 267, 108033. [Google Scholar] [CrossRef]
  32. Qiu, P.; Qin, Y.; Bai, Y.; Xia, Q.; Zheng, A. Gas Selectivity Regulation of Monolayer SnS by Introducing Nonmetallic Dopants: A Combined Theoretical and Experimental Investigation. Appl. Surf. Sci. 2021, 570, 151155. [Google Scholar] [CrossRef]
  33. Sankush Krishna, M.; Singh, S.; Mohammed, M.K.A. Carcinogenic Heavy Metals Detection Based on ZnO Nanoribbons. IEEE Sens. J. 2022, 22, 16929–16937. [Google Scholar] [CrossRef]
  34. Wu, J.; Li, X.; Huang, L.; Liang, T.; Xing, X.; Luo, A. SnSe Monolayer-Based Heavy Metal Sensors with High Sensitivity, Selectivity, and Reusability: Insights from First Principle Calculation. Results Phys. 2023, 53, 106973. [Google Scholar] [CrossRef]
  35. Wu, J.; Li, X.; Liao, H.; Xue, S.; Huang, L.; Xing, X.; Luo, A. SnS Monolayers Based Heavy Metal Sensors: DFT and NEGF Analysis. Results Phys. 2023, 51, 106704. [Google Scholar] [CrossRef]
  36. Srivastava, M.; Srivastava, A. Cu Decorated Functionalized Graphene for Arsenic Sensing in Water: A First Principles Analysis. Appl. Surf. Sci. 2021, 560, 149700. [Google Scholar] [CrossRef]
  37. Srivastava, M.; Srivastava, A.; Pandey, S.K. Suitability of Graphene Monolayer as Sensor for Carcinogenic Heavy Metals in Water: A DFT Investigation. Appl. Surf. Sci. 2020, 517, 146021. [Google Scholar] [CrossRef]
  38. Srivastava, M.; Srivastava, A. DFT Analysis of Nitrogen and Boron Doped Graphene Sheet as Lead Detector. Mater. Sci. Eng. B 2021, 269, 115165. [Google Scholar] [CrossRef]
  39. Shtepliuk, I.; Yakimova, R. Interaction of Epitaxial Graphene with Heavy Metals: Towards Novel Sensing Platform. Nanotechnology 2019, 30, 294002. [Google Scholar] [CrossRef]
  40. Yang, T.; Jiang, X.; Yi, W.; Cheng, X.; Cheng, T. Enhanced Fast Response to Hg0 by Adsorption-Induced Electronic Structure Evolution of Ti2C Nanosheet. Appl. Surf. Sci. 2021, 544, 148925. [Google Scholar] [CrossRef]
Figure 1. (a) The top side view of MoS2 monolayer with the position sites (H, B, TS, and TMo); (b) the side view of MoS2 monolayer.
Figure 1. (a) The top side view of MoS2 monolayer with the position sites (H, B, TS, and TMo); (b) the side view of MoS2 monolayer.
Micromachines 15 00978 g001
Figure 2. The most stable configuration and Bader charge analysis for (a) As, (b) Cd, (c) Hg, and (d) Pb on the MoS2 monolayer. Negative numbers indicate loss of electrons, positive numbers indicate gain of electrons.
Figure 2. The most stable configuration and Bader charge analysis for (a) As, (b) Cd, (c) Hg, and (d) Pb on the MoS2 monolayer. Negative numbers indicate loss of electrons, positive numbers indicate gain of electrons.
Micromachines 15 00978 g002
Figure 3. Band structures of (a) As, (b) Cd, (c) Hg, and (d) Pb adsorbed onto the MoS2 monolayer.
Figure 3. Band structures of (a) As, (b) Cd, (c) Hg, and (d) Pb adsorbed onto the MoS2 monolayer.
Micromachines 15 00978 g003
Figure 4. Total and partial DOS of MoS2 monolayer: (a) As, (b) Cd, (c) Hg, and (d) Pb adsorption (Fermi level is set as 0 eV).
Figure 4. Total and partial DOS of MoS2 monolayer: (a) As, (b) Cd, (c) Hg, and (d) Pb adsorption (Fermi level is set as 0 eV).
Micromachines 15 00978 g004
Figure 5. Charge density differential results for (a) As, (b) Cd, (c) Hg, and (d) Pb adsorbed on monolayer MoS2. The isosurface value is taken as 0.0003 e/Å3.
Figure 5. Charge density differential results for (a) As, (b) Cd, (c) Hg, and (d) Pb adsorbed on monolayer MoS2. The isosurface value is taken as 0.0003 e/Å3.
Micromachines 15 00978 g005
Figure 6. The recovery time (τ) for As, Cd, Hg, and Pb adsorbed on MoS2 monolayers at different temperatures (T).
Figure 6. The recovery time (τ) for As, Cd, Hg, and Pb adsorbed on MoS2 monolayers at different temperatures (T).
Micromachines 15 00978 g006
Figure 7. The current-voltage characteristics of MoS2 monolayers with As, Cd, Hg, and Pb adsorption.
Figure 7. The current-voltage characteristics of MoS2 monolayers with As, Cd, Hg, and Pb adsorption.
Micromachines 15 00978 g007
Figure 8. The sensitivity for As, Cd, Hg, and Pb adsorbed on MoS2 monolayers.
Figure 8. The sensitivity for As, Cd, Hg, and Pb adsorbed on MoS2 monolayers.
Micromachines 15 00978 g008
Figure 9. The comparison of sensitivity for (a) As, (b) Cd, (c) Hg, and (d) Pb sensing materials.
Figure 9. The comparison of sensitivity for (a) As, (b) Cd, (c) Hg, and (d) Pb sensing materials.
Micromachines 15 00978 g009
Table 1. Adsorption energy, adsorption distance, charge transfer, and energy gap.
Table 1. Adsorption energy, adsorption distance, charge transfer, and energy gap.
Heavy
Metal
Adsorption Energy
Ead (eV)
Adsorption Distance
D (Å) 1
Charge Transfer
△Q (e)
StyleEnergy Gap
Eg (eV)
As−0.8272.175−0.243Donor0.242
Cd−0.4463.358−0.115Donor1.140
Hg−0.3663.526−0.023Donor1.534
Pb−1.3802.625−0.410Donor0.326
1 The adsorption distance is defined as the closest contact distance between the adsorbed atom and the nearest atom on the MoS2 surface.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, J.; Li, Z.; Liang, T.; Mo, Q.; Wei, J.; Li, B.; Xing, X. First-Principles Insights into Highly Sensitive and Reusable MoS2 Monolayers for Heavy Metal Detection. Micromachines 2024, 15, 978. https://doi.org/10.3390/mi15080978

AMA Style

Wu J, Li Z, Liang T, Mo Q, Wei J, Li B, Xing X. First-Principles Insights into Highly Sensitive and Reusable MoS2 Monolayers for Heavy Metal Detection. Micromachines. 2024; 15(8):978. https://doi.org/10.3390/mi15080978

Chicago/Turabian Style

Wu, Jiayin, Zongbao Li, Tongle Liang, Qiuyan Mo, Jingting Wei, Bin Li, and Xiaobo Xing. 2024. "First-Principles Insights into Highly Sensitive and Reusable MoS2 Monolayers for Heavy Metal Detection" Micromachines 15, no. 8: 978. https://doi.org/10.3390/mi15080978

APA Style

Wu, J., Li, Z., Liang, T., Mo, Q., Wei, J., Li, B., & Xing, X. (2024). First-Principles Insights into Highly Sensitive and Reusable MoS2 Monolayers for Heavy Metal Detection. Micromachines, 15(8), 978. https://doi.org/10.3390/mi15080978

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop