Inkjet-Printed Molybdenum Disulfide and Nitrogen-Doped Graphene Active Layer High On/Off Ratio Transistors.

Fully inkjet-printed device fabrication is a crucial goal to enable large-area printed electronics. The limited number of two-dimensional (2D) material inks, the bottom-gated structures, and the low current on/off ratio of thin-film transistors (TFTs) has impeded the practical applications of the printed 2D material TFTs. In the search for TFTs with high current ratios, we introduce a stable and efficient method of nitrogen-doped graphene (NDG) ink preparation for inkjet printing by liquid-phase exfoliation. The NDG thin film is print-stacked with molybdenum disulfide (MoS2) by multiple printing passes to construct a MoS2–NDG stack. We demonstrate top-gated fully inkjet-printed MoS2–NDG transistors with silver drain, source, and gate electrodes, and a barium titanate (BaTiO3) dielectric. A 100% inkjet-printed MoS2–NDG vertical 2D active heterostructure layer transistor with a current on/off ratio of 1200 is exhibited. The results may lead towards the development of all-printed 2D material-based transistor switches.


Introduction
Graphene and related two-dimensional (2D) materials remain at the center of research for more than a decade because of their exciting and unusual properties [1]. Printing of solution-processed 2D materials facilitates the fabrication of the next generation of flexible devices. Inkjet printing as a process for flexible and printed electronics, has become more popular for the fabrication of electronic devices because of its ability for large-area fabrication, the limited number of process steps involved, low-temperature processing, and low cost [2]. Graphene is a zero-band gap semimetal in which the valence and conduction bands are cone shaped, and they meet at the K point of the Brillouin zone. Due to the absence of band gap, the typical current on/off ratio of pristine graphene thin-film transistors (TFTs) is approximately 5 [3]. For switching applications of TFTs, a current on/off ratio between 10 4 and 10 7 and a band gap of thin film larger than 0.36 eV are desirable [4]. As with regular thin-film graphene TFTs, the inkjet-printed graphene transistors had current ratios of ≈10 [2,5]. Recently, all inkjet-printed graphene transistors were reported for wearable and textile electronics [6]. However, the current on/off ratio was still small (about 2.5). Molybdenum disulfide (MoS 2 ), a semiconducting layered transition metal dichalcogenide (TMD), is one of the most studied 2D materials after graphene. Mono and few-layered MoS 2 have a direct band gap of 1.8 eV, while bulk MoS 2 has an indirect band gap of 1.29 eV [7]. Despite having a band gap in MoS 2 , liquid phase-exfoliated (LPE) and inkjet-printed MoS 2

NDG and MoS 2 Nanosheets
The flake dimensions must match the requirements of drop-on-demand (DOD) printing. The lateral dimensions of the dispersed nanosheets should be at least 1/50 of the nozzle diameter (a ≈ 21.5 µm) to prevent the nozzle clogging. We purchased the NDG powder from Sigma-Aldrich (product no. 900527) and MoS 2 ultrafine powder from Graphene supermarket. The scanning electron microscopy (SEM) images of bulk nanosheets are shown in Figure 1a-b. Prior to dispersing the nanosheets into solvents, the nanosheets were sonicated for 10 h to break the bulk nanosheets. The lateral size and thickness of the individual NDG and MoS 2 nanosheets were extracted from atomic force microscopy (AFM) measurements for over 50 flakes as shown in Figure 1c-f. The flake size follows the normal distribution (peak ≈ 150 nm) for NDG and log-normal distribution (peak ≈ 100 nm) for MoS 2 , respectively, which fulfills the DOD requirements. The statistical thickness of normal distribution for NDG nanosheets peaked at 20 nm, and the Lorentz distribution of thickness for MoS 2 nanosheets peaked at 25 nm, respectively. Assuming a thickness of 0.37 nm for monolayer graphene and 0.615 nm for monolayer MoS 2 [17,18], the nanosheets consist of an average number of layers; N ≈ 54 for the NDG and N ≈ 40 for the MoS 2 . Raman is one of the most common vibrational spectroscopies for fingerprinting carbon species. After the sonication, the Raman spectrum of NDG nanosheets on Si substrate was acquired under ambient conditions, which is shown in Figure 2a. The NDG shows two intense peaks, which are the D (≈1341 cm −1 ) and G (≈1578 cm −1 ) peaks, respectively. There are two additional peaks, 2D (≈2670 cm −1 ) and D + D' (≈2904 cm −1 ), present in the NDG nanosheets. The G and 2D peaks are the signature of graphene-like structure in nanosheets. The G band corresponds to E 2g phonons at the Brillouin zone center. The G band arises from the stretching of C-C bonds in a sp 2 hybridized carbon system. The D band is related to defects in sp 2 hybridized carbon materials, which is activated by an inter-valley double resonance (DR) Raman process. The intensity ratio I D (intensity of D peak)/I G (intensity of G peak) of 1.005 indicates a high degree of disorder due to nitrogen doping. The strong D peak can also originate from the submicrometer flake boundaries and vacancies. The 2D peak is the secondary D peak, which appears at approximately 2700 cm −1 . For single-layer graphene, the 2D band shows a sharp peak and high intensity compared to the D and G peaks [19]. The 2D band broadens and the intensity reduces with the number of graphene layers. In the Raman spectrum of Figure 2a, the intensity of the 2D band is much smaller compared to the D and G bands. We attribute the low intensity of the 2D band to the presence of many layers (N ≈ 54) in NDG nanosheets. This Raman spectrum strongly agrees with the Raman spectra previously reported for NDG synthesized by chemical vapor deposition at 800 and 900 • C [20]. Figure 2b shows the Raman spectrum of MoS 2 flakes with characteristic E 1g and A 1g peaks at 378 and 405 cm −1 respectively, which are obtained by Lorentzian fits in the spectrum. Raman is one of the most common vibrational spectroscopies for fingerprinting carbon species. After the sonication, the Raman spectrum of NDG nanosheets on Si substrate was acquired under ambient conditions, which is shown in Figure 2a. The NDG shows two intense peaks, which are the D (≈ 1341 cm −1 ) and G (≈ 1578 cm −1 ) peaks, respectively. There are two additional peaks, 2D (≈ 2670 cm −1 ) and D + D' (≈ 2904 cm −1 ), present in the NDG nanosheets. The G and 2D peaks are the signature of graphene-like structure in nanosheets. The G band corresponds to E2g phonons at the Brillouin zone center. The G band arises from the stretching of C-C bonds in a sp 2 hybridized carbon system. The D band is related to defects in sp 2 hybridized carbon materials, which is activated by an intervalley double resonance (DR) Raman process. The intensity ratio ID (intensity of D peak)/IG (intensity of G peak) of 1.005 indicates a high degree of disorder due to nitrogen doping. The strong D peak can also originate from the submicrometer flake boundaries and vacancies. The 2D peak is the secondary D peak, which appears at approximately 2700 cm −1 . For single-layer graphene, the 2D band shows a sharp peak and high intensity compared to the D and G peaks [19]. The 2D band broadens and the intensity reduces with the number of graphene layers. In the Raman spectrum of Figure 2a, the intensity of the 2D band is much smaller compared to the D and G bands. We attribute the low intensity of the 2D band to the presence of many layers (N ≈ 54) in NDG nanosheets. This Raman spectrum strongly agrees with the Raman spectra previously reported for NDG synthesized by chemical vapor deposition at 800 and 900 °C [20]. Figure 2b shows the Raman spectrum of MoS2 flakes with characteristic E1g and A1g peaks at 378 and 405 cm −1 respectively, which are obtained by Lorentzian fits in the spectrum.

Printable Inks and Thin Films
In this study, we used a commercial drop-on-demand (DOD) Fujifilm Dimatix materials printer (DMP-2800). For DOD inkjet printing, the material inks must have the ability to generate droplets. The droplets ejecting out from a nozzle are influenced by ink viscosity η (mPa s), surface tension γ

Printable Inks and Thin Films
In this study, we used a commercial drop-on-demand (DOD) Fujifilm Dimatix materials printer (DMP-2800). For DOD inkjet printing, the material inks must have the ability to generate droplets. The droplets ejecting out from a nozzle are influenced by ink viscosity η (mPa s), surface tension γ (mN/m), particle density ρ (g/cm 3 ), and nozzle diameter a (µm). We prepared and characterized the NDG and MoS 2 inks (see 'Materials and Methods'). For inkjet printing of an ink, the inverse Ohnesorge number (Z = Oh −1 ) is used as the figure of merit (FOM) to characterize the drop formation and jettability of an ink. Z is calculated by using Z = (γρa) 1/2 η and an optimal range of Z between 1 and 14 is required for stable DOD printing [5]. For NDG ink, using η = 2.94 mPa s, γ = 20.98 mN/m, ρ = 1.054 g/cm −3 , a = 21.5 µm, we obtain Z = 7.5, which is within the conventional range for inkjet printing. There were no satellite drops following the primary drops, which was confirmed by capturing the dynamics of drop formation. The Z for MoS 2 and barium titanate (BaTiO 3 ) inks were 2 and 3.4, which are suitable for inkjet printing.
The concentration of nanosheets in the ink can be calculated from the UV-vis absorbance spectrum via the Beer-Lambert law. The Beer-Lambert law can be expressed as A = αcl, where A is the absorbance, α is the absorption coefficient (Lg −1 m −1 ), c is the concentration of dispersed graphene (g/L), and l is the light path length (m). The optical absorption of the NDG ink is shown in Figure 3a. The spectrum is mostly featureless due to the linear dispersion of the Dirac electrons, whereas the peak in the UV region is a signature of Hove singularity in the graphene density of states [6]. From A = 0.501, α = 2460 Lg −1 m −1 , l = 0.01 m, our estimated NDG concentration is 0.4 mg/mL [21]. Figure 3b shows the optical absorbance spectra for as prepared MoS 2 ink. This absorbance spectrum complies with the previously reported MoS 2 ink absorbance spectra [8]. For MoS 2 presence in dispersion, two characteristic peaks appear at approximately 600 and 672 nm wavelength respectively [10]. Two excitonic absorption peaks at 605 and 666 nm are observed, which arise from the K points of the Brillouin zone in MoS 2 nanosheets. The peak at 666 nm corresponds to the lowest optical bandgap of 1.86 eV for MoS 2 nanosheets, which is higher than the bulk MoS 2 bandgap of 1.3 eV [22]. The extinction coefficient of MoS 2 at 672 nm wavelength is~α 672 = 3400 mL/(mg.m) [8]. Using the Beer-Lambert law, the calculated final MoS 2 concentration is~0.1 mg/mL. (mN/m), particle density ρ (g/cm 3 ), and nozzle diameter a (µ m). We prepared and characterized the NDG and MoS2 inks (see 'Materials and Methods'). For inkjet printing of an ink, the inverse Ohnesorge number (Z = Oh −1 ) is used as the figure of merit (FOM) to characterize the drop formation and jettability of an ink. Z is calculated by using Z = (γρa) 1/2 η ⁄ and an optimal range of Z between 1 and 14 is required for stable DOD printing [5]. For NDG ink, using η = 2.94 mPa s, γ = 20.98 mN/m, ρ = 1.054 g/cm −3 , a = 21.5 µ m, we obtain Z = 7.5, which is within the conventional range for inkjet printing. There were no satellite drops following the primary drops, which was confirmed by capturing the dynamics of drop formation. The Z for MoS2 and barium titanate (BaTiO3) inks were 2 and 3.4, which are suitable for inkjet printing.
The concentration of nanosheets in the ink can be calculated from the UV-vis absorbance spectrum via the Beer-Lambert law. The Beer-Lambert law can be expressed as A = αcl, where A is the absorbance, α is the absorption coefficient (Lg −1 m −1 ), c is the concentration of dispersed graphene (g/L), and l is the light path length (m). The optical absorption of the NDG ink is shown in Figure 3a. The spectrum is mostly featureless due to the linear dispersion of the Dirac electrons, whereas the peak in the UV region is a signature of Hove singularity in the graphene density of states [6]. From A = 0.501, α = 2460 Lg −1 m −1 , l = 0.01 m, our estimated NDG concentration is 0.4 mg/mL [21] 1 . Figure  3b shows the optical absorbance spectra for as prepared MoS2 ink. This absorbance spectrum complies with the previously reported MoS2 ink absorbance spectra [8]. For MoS2 presence in dispersion, two characteristic peaks appear at approximately 600 and 672 nm wavelength respectively [10]. Two excitonic absorption peaks at 605 and 666 nm are observed, which arise from the K points of the Brillouin zone in MoS2 nanosheets. The peak at 666 nm corresponds to the lowest optical bandgap of 1.86 eV for MoS2 nanosheets, which is higher than the bulk MoS2 bandgap of 1.3 eV [22]. The extinction coefficient of MoS2 at 672 nm wavelength is ~α672 = 3400 mL/(mg.m) [8]. Using the Beer-Lambert law, the calculated final MoS2 concentration is ~ 0.1 mg/mL.    We also calculated the full width at half maximum (FWHM) of the 2D peaks for the Raman spectra shown in Figure 4. The FWHM for the NDG powder, 2, 5 and 10 printed layers was 77.4, 63.6, 62.6, and 92.9, respectively. This, along with 2D peak position shift, implies that the nitrogen doping has a strong effect on printed NDG samples.
Molecules 2020, 25, x FOR PEER REVIEW 6 of 13 peak indicates the presence of many layers of graphene in the printed thin films. Doping has strong effects on the position and FWHM of the 2D peak. The 2D peak position varied between 2667 and 2711 cm −1 . We also calculated the full width at half maximum (FWHM) of the 2D peaks for the Raman spectra shown in Figure 4. The FWHM for the NDG powder, 2, 5 and 10 printed layers was 77.4, 63.6, 62.6, and 92.9, respectively. This, along with 2D peak position shift, implies that the nitrogen doping has a strong effect on printed NDG samples.

MoS2-NDG Stack
The 2D materials can be arranged into a stack to create a materials assembly with novel properties [23]. In the stack, atomically thin 2D materials are bonded by weak van der Walls forces [24]. The MoS2 dispersions produced by LPE have the optical bandgap ranges from 1.8 eV to 1.3 eV [12]. We created and optimized the stack structures of printed MoS2 and NDG materials. To inspect the inside of a stack, a cross-sectional slice was produced by focused ion beam (FIB) milling for a thicker sample. Figure 5a shows the FIB-SEM cross-sectional image of the MoS2-NDG stack. The cross-sectional thickness of the NDG and MoS2 prints are 527 and 105 nm, respectively. The interfacial separation between the MoS2 and NDG is not abrupt and well-maintained along the cross-section. There is no identifiable large flake in the cross-section. This ensures that the nanosheets are dissolved well in the solvents, and our sonication and centrifugation processes are reliable in removing thicker flakes. Raman of the MoS2-NDG stack in Figure 5b shows the characteristic MoS2 peaks (E2g ∼ 378 cm −1 , A1g ∼ 403 cm −1 ) and NDG peaks (D ∼ 1357 cm −1 , G ∼ 1590 cm −1 , and 2D ∼ 2681 cm −1 ).

MoS 2 -NDG Stack
The 2D materials can be arranged into a stack to create a materials assembly with novel properties [23]. In the stack, atomically thin 2D materials are bonded by weak van der Walls forces [24]. The MoS 2 dispersions produced by LPE have the optical bandgap ranges from 1.8 eV to 1.3 eV [12]. We created and optimized the stack structures of printed MoS 2 and NDG materials. To inspect the inside of a stack, a cross-sectional slice was produced by focused ion beam (FIB) milling for a thicker sample. Figure 5a shows the FIB-SEM cross-sectional image of the MoS 2 -NDG stack. The cross-sectional thickness of the NDG and MoS 2 prints are 527 and 105 nm, respectively. The interfacial separation between the MoS 2 and NDG is not abrupt and well-maintained along the cross-section. There is no identifiable large flake in the cross-section. This ensures that the nanosheets are dissolved well in the solvents, and our sonication and centrifugation processes are reliable in removing thicker flakes. Raman of the MoS 2 -NDG stack in Figure 5b shows the characteristic MoS 2 peaks (E 2g~3 78 cm −1 , A 1g~4 03 cm −1 ) and NDG peaks (D~1357 cm −1 , G~1590 cm −1 , and 2D~2681 cm −1 ).

Transistor Fabrication and Characterization
The transistors themselves are all inkjet printed on the glass substrate without involving any photolithography patterning or surface pretreatment steps. The TFTs are top gated, consisting of the source, drain, gate electrodes, a channel, and a dielectric layer. The inset of Figure 6a shows the schematic structure of an NDG transistor. Highly conductive silver nanoparticle ink (Novacentrix Metalon ® JS-B40G) was used to print source, drain and gate terminals. As the very first layer, source and drain electrodes were printed on the glass substrate in a printing pass. The printed silver patterns were cured at 150 °C for 30 min. The separation between source and drain electrodes, i.e., the channel length (L) was ∼ 80 µ m. The ed silver contacts deposited uniformly with an average thickness of 0.50 µ m. For channel region, a percolation network of 100 nm NDG was printed with 40 printing passes followed by a MoS2 printing repeated 4 times on top of NDG. The channel layer had some overlapping with source and drain contacts to avoid discontinuity in the structure. Barium titanate (BaTiO3) ink (k ∼ 20.5) of a thickness of ∼ 2 µ m was used as the dielectric (see 'Materials and Methods'). The printed dielectric was first baked at 100 °C for 20 min, and subsequently baked at 230 °C for 30 min. Finally, the silver gate electrode was printed on top of the dielectric layer.

Transistor Fabrication and Characterization
The transistors themselves are all inkjet printed on the glass substrate without involving any photolithography patterning or surface pretreatment steps. The TFTs are top gated, consisting of the source, drain, gate electrodes, a channel, and a dielectric layer. The inset of Figure 6a shows the schematic structure of an NDG transistor. Highly conductive silver nanoparticle ink (Novacentrix Metalon ® JS-B40G) was used to print source, drain and gate terminals. As the very first layer, source and drain electrodes were printed on the glass substrate in a printing pass. The printed silver patterns were cured at 150 • C for 30 min. The separation between source and drain electrodes, i.e., the channel length (L) was~80 µm. The ed silver contacts deposited uniformly with an average thickness of 0.50 µm. For channel region, a percolation network of 100 nm NDG was printed with 40 printing passes followed by a MoS 2 printing repeated 4 times on top of NDG. The channel layer had some overlapping with source and drain contacts to avoid discontinuity in the structure. Barium titanate (BaTiO 3 ) ink (k~20.5) of a thickness of~2 µm was used as the dielectric (see 'Materials and Methods'). The printed dielectric was first baked at 100 • C for 20 min, and subsequently baked at 230 • C for 30 min.   Figure 6a shows the non-linear transfer curve of the transistor at V D = −2 V. The applied gate voltage changed from −20 to 10 V. The I D increases as the gate voltage changes from positive to negative, which suggests that the stack creates a p-type transport layer. A current on/off ratio of 1200 is calculated at V D = −2 V from the transfer characteristics. The electric field mobility (µ = L×g m W×C ox ×V D ) is calculated under the assumption of constant electric field along the channel and no velocity saturation. The device mobility was~0.01 cm 2 /V.s, which is comparable with liquid-exfoliated MoS 2 -based TFTs [8]. Moreover, the device has a good off state. We measured the gate leakage (I G ) from −10 to +10 V. The I G shows random signal in this range and does not indicate any gate leakage and the drain current (I D ) was independent of the gate current (Supplementary Information). Figure 6b shows the output characteristics of the device at a gate voltage sweep from 10 to −20 V. There is no Schottky barrier and the current started flowing upon applying V D implying the contact is ohmic and is the result of direct tunneling process. The transfer and output characteristics of the device exhibit p-type behaviors. To benchmark the on/off ratio results, we also investigated the inkjet-printed NDG transistors. The similar fabrication procedures were carried out for NDG transistors (L~80 µm) on the glass substrate with same Ag electrodes and BaTiO 3 gate dielectric. For NDG transistors, a channel of 100 nm thickness was created with 40 printing repetitions to ensure the availability of NDG nanosheets on the top surface of the NDG thin films. The channel layer had some overlapping with source and drain contacts to avoid discontinuity in the structure. Transfer characteristics of the NDG transistor are shown in Figure S4a. The gate voltages were changed from −20 to +20 V. The transfer characteristics were measured at V D = 5, 10, 15, and 20 V, respectively. The drain current increased as the V GS changed from negative to positive, which confirmed the n-type doping in NDG nanosheets due to nitrogen doping. However, the drain current did not change significantly over a wide range of gate bias. For example, at V D = 20 V, the drain current changed from 1105 to 1247 mA in the given gate voltage range. The current on/off ratio can be calculated from the transfer characteristics and it is the ratio of maximum to minimum drain current for a specific V D . At V D = 20, the on/off ratio of the NDG is~1.13. This on/off ratio is comparable to other inkjet-printed graphene TFTs, but not suitable for transistor switches. The drain currents were measured in V GS = −20, −10, 0, 10, and 20 V, respectively. Still the linear I D -V D curves were not distinguishable for various gate voltages, i.e., the field effect in the output curves were not apparent. A linear I D -V D curve for V GS = 10 V is displayed in Figure S4b as the output characteristics. We also fabricated several MoS 2 transistors following same fabrication processes for a channel length of~80 µm. However, we were unsuccessful in achieving a working transistor. Our interpretation in this matter is that the printed MoS 2 layers failed to form a semiconducting percolation network on a bare glass substrate. However, MoS 2 layers when printed on top of NDG layers, formed a continuous MoS 2 percolation network and maintained a clear separation from the NDG layers underneath due to hydrophobic nature of NDG surface. Some MoS 2 flakes may penetrate to the NDG layers, which reduced the overall conductivity of the NDG channel. Therefore, the high current on/off ratio originated from printed MoS 2 and the NDG layers underneath supported MoS 2 to start the nucleation when printed.

Discussion
We demonstrated a high current on/off ratio transistor based on the MoS 2 -NDG stack as the channel. The inks, thin films, and devices were characterized by various experimental characterization techniques. The NDG transistor was used as a control structure which showed low on/off ratios. Upon printing MoS 2 on top of NDG to create a MoS 2 -NDG stack, semiconducting properties were achieved, and on/off ratios increased significantly for a channel length of 80 µm.

Materials and Methods
Nitrogen-doped graphene (NDG) ink formulation. The NDG ink was prepared by dispersing NDG flakes into terpineol/cyclohexanone mixture. The NDG powder (≈5 mg, Sigma-Aldrich product no. 900527, Saint Louis, MO 63103, USA) was sonicated for 10 h to break the bigger particles. For a 10 mL solution, 1.5 mL of terpineol was used as the solvent. To curtail the viscosity of terpineol, 8.5 mL of cyclohexanone was added. Then, ethyl cellulose (as a surfactant) of 48 mg was mixed into the solution and sonicated for 1 h. Finally, the sonicated NDG powder (≈5 mg) was poured into the solution. The resulting suspension was further sonicated for 20 h, and then centrifuged at 5000 rpm for 30 min to remove the thicker flakes. After centrifugation, the suspension was left to settle overnight. The final ink was filtered out from the residue and ready for characterization.
Nitrogen-doped graphene (NDG) ink characterization. The surface tension of NDG ink was measured using the pendant drop method (First Ten Angstroms FTA10000B, Portsmouth, VA 23704, USA). A small amount of ink (3.8 µL) was suspended from a needle and surface tension was extracted from the software image. The ink density is ρ = m/V, m is the mass of ink in grams and V is volume of ink in mL. The ink density ρ for NDG ink was 1.054 g/cm 3 . The viscosity of the ink was measured with the RheoSense m-VROC viscometer. The viscosity of NDG ink was ≈2.94 mPa s. Similar characterization was carried out for the MoS 2 ink.
Atomic force microscopy (AFM). A Bruker atomic force microscope (Billerica, MA 01821, USA) was used to extract the flake thickness and lateral distribution. The flakes were scraped on a silicon substrate after sonication before mixing with solvents. The scan areas were 500 × 500 µm 2 and 200 × 200 µm 2 for the nanosheets.
Scanning electron microscopy. Scanning electron microscopy was carried out using a Helios Nanolab 400 FEI SEM (Hillsboro, OR 97124, USA). The field emission gun was operated at an accelerating voltage of 5 kV and gun current of 0.34 nA. The images were captured in the field-free mode of SEM using TLD detector at the indicated magnification scale in SEM images.
Raman spectroscopy. Raman measurements were performed on the sonicated bulk NDG powder and inkjet-printed layers using a Horiba LabRAM HR evolution Raman spectrometer (Austin, TX 78754, USA). The laser excitation of 532 nm was used at 100 × objective, with an incident power less than 1 mW to avoid sample damage.
Molybdenum disulfide (MoS 2 ) ink. The MoS 2 ultrafine powder was purchased from Graphene Supermarket (Ronkonkoma, NY 11779, USA). The NMP (1-Methyl-2-pyrrolidinone) solvent was diluted to DI water by a 45:55 mL ratio. The weight of NMP/DI water composition was measured and mixed with 2% weight equivalent of MoS 2 powder. The mixture was sonicated for 20 h and followed by a centrifugation at 5000 rpm for 30 min. The supernatant was filtered out from the top of the mixture after it rested overnight. The resultant MoS 2 ink was yellowish. The ink viscosity, surface tension, and density were 3.8 mPa s, 25 mN/m, and 0.1 g/cm 3 , respectively.
Optical absorption spectroscopy. UV-vis absorption spectroscopy was used to calculate the flakes concentration in NDG and MoS 2 inks utilizing the Beer-Lambert law. For the absorbance spectra measurements, the NDG ink was diluted in isopropyl alcohol (IPA) at a ratio of 1:20 and the MoS 2 ink was as prepared. The measurements were performed with a Shimadzu UV-2501 spectrophotometer (Kyoto, Japan) using a glass cuvette with an optical path length of 10 mm. For concentration calculation, absorption coefficients of α GR~2 460 for the graphene ink at 660 nm wavelength [21], and α MoS2~3 400 for the MoS 2 ink [8] at 672 nm wavelength were used.
Barium titanate (BaTiO 3 ) ink preparation. The dielectric ink includes 20% by weight of BaTiO 3 nanoparticles, 1% by weight of the Triton X-100, 14.5% by weight of Cyclopentanone, and 64.5% by weight of SU8 2005 photoresist. For ink preparation, SU8 2005 (MicroChem, Newton, MA 02464, USA), Cyclopen-tanone (Sigma-Aldrich, product-W391018), and Triton X-100 (Sigma-Aldrich, product-X100) were mixed and sonicated for 4 h. The BaTiO 3 nanoparticles (US Research Nanomaterials, product-US3827, Houston, TX 77084, USA) are then added to the mixture, and the composition was sonicated for an additional 6 h. The resultant mixture was left untouched for one day to form a stable dispersion. The ink was collected from the suspension, leaving the residue of thick powders at the bottom of the bottle. The density of BaTiO 3 was measured to be 1.4 g/cm 3 . The surface tension of 55 mN/m and viscosity of 12 mPa s were considered for inverse Ohnesorge number calculation [25]. The effective dielectric constant (k) of the printed dielectric was measured from the transmission coefficient S 12 with a T-shape microstrip line resonator. The effective dielectric constant (k) was~20.5. The k-value of printed dielectric is consistent with the previous report [26].
Viscosity measurement. The viscosity was measured at room temperature with a Rheosense m-VROC viscometer (San Ramon, CA 94583, USA). In m-VROC Rheosense software (version: 3.1.4), we entered the sample size in mL and an estimated viscosity for an ink and it automatically calculated minimum and maximum flow rate. We set five different flow rates within the displayed range and measurement times and dispersed the sample through the measurement sensor. The regression coefficient (R 2 ) was greater than 0.99 for all our viscosity measurements Inkjet Printing. A DoD Fujifilm Dimatix (DMP-2800, Santa Clara, CA 95050, USA) inkjet printer was used to print all the thin films and devices. It was equipped with a Dimatix materials cartridge (DMC-11610). The cartridge head has 16 piezoelectric jetting nozzles with a diameter of 21.5 µm each. The nozzles can dispense a droplet of the nominal volume of 10 pL from the cartridge head. The optimum drop spacing for the NDG and MoS2 inks was 30 µm, and six nozzles were used to print these inks. The BaTiO 3 ink was printed with a 35 µm drop spacing with one nozzle. There was no delay between subsequent printing passes except a 0.3 s purge in every ten printing passes. Additionally, to avoid the coffee ring effect, the platen temperature was set at 60 • C. The 2D device patterns were created in AutoCAD, where the design layout represents different parts of the transistor. The design patterns were converted into a printable version of the Fujifilm Dimatix printer with ACE 300 software.
Electrical characterization. The current-voltage measurements of the devices were performed with the Keysight B1500A semiconductor device analyzer.

Summary
In summary, we developed the NDG and MoS 2 inks for inkjet printing and studied the inks and thin-film properties through experimental characterizations. We devised a way to prepare a MoS 2 -NDG stack to use as the transistor channel. The devices were fully inkjet-printed without any involvement of Cleanroom facility. Based on the test results of the DC analysis of the device, the current on/off ratio of 1200 is reported in this work.
Supplementary Materials: The following are available online. Table S1: Data table for Raman measurements; Figure S1: Optical image of final transistor; Figure S2: Absorbance spectrum for solvents of (a) NDG ink, and (b) MoS 2 ink; Figure S3: Gate leakage measurement of the MoS 2 -NDG transistor; Figure S4: (a) Transfer characteristics, and (b) output curve of the NDG transistor; Table S2: Comparison of on/off ratios of 2D material TFTs; Figure S5: Comparison of current on/off ratios.