Fabrication and microfluidic analysis of graphene-based molecular communication receiver for Internet of Nano Things (IoNT).

Bio-inspired molecular communications (MC), where molecules are used to transfer information, is the most promising technique to realise the Internet of Nano Things (IoNT), thanks to its inherent biocompatibility, energy-efficiency, and reliability in physiologically-relevant environments. Despite a substantial body of theoretical work concerning MC, the lack of practical micro/nanoscale MC devices and MC testbeds has led researchers to make overly simplifying assumptions about the implications of the channel conditions and the physical architectures of the practical transceivers in developing theoretical models and devising communication methods for MC. On the other hand, MC imposes unique challenges resulting from the highly complex, nonlinear, time-varying channel properties that cannot be always tackled by conventional information and communication tools and technologies (ICT). As a result, the reliability of the existing MC methods, which are mostly adopted from electromagnetic communications and not validated with practical testbeds, is highly questionable. As the first step to remove this discrepancy, in this study, we report on the fabrication of a nanoscale MC receiver based on graphene field-effect transistor biosensors. We perform its ICT characterisation in a custom-designed microfluidic MC system with the information encoded into the concentration of single-stranded DNA molecules. This experimental platform is the first practical implementation of a micro/nanoscale MC system with nanoscale MC receivers, and can serve as a testbed for developing realistic MC methods and IoNT applications.

Introduction

Nanotechnology is enabling us to devise ever-smaller devices to interact with the universe at molecular resolution. Though small enough to penetrate into cells, these devices individually are of limited capability. To unleash the full potential of nanotechnology, communication among nanomachines is a must. This would enable more complex applications, e.g., continuous health monitoring with intrabody nanosensor networks, theranostic applications with distributed sensor and actuator networks, and industrial nanosensor applications. Internet of Nano Things (IoNT), which defines these networks of nanomachines, such as nanobiosensors and engineered bacteria, integrated with the Internet infrastructure, has seen tremendous interest recently[1-3]. IoNT is set to transform the way we connect with and understand the world at the bottom. However, conventional electromagnetic (EM) communication techniques, proved impractical at nanoscale due to the antenna size and power limitations, and large propagation losses. On the other hand, Nature itself already provides a robust way of nanocommunication, i.e., Molecular Communications (MC), which is the common communication modality among living cells, ranging from our own neurons to bacteria[4]. Using the same language with living cells by encoding, transmitting and receiving information with molecules provides an energy efficient and reliable nanocommunication, even at harsh biological environments, where, we expect, the most impactful medical applications of IoNT would be implemented.

Although there is tremendous interest in this field accompanied by a large body of theoretical work to develop models and devise communication methods for MC[5-11], researchers rely on overly simplifying assumptions about the physical constraints of the transceivers and the channel conditions due to the lack of a testbed for MC at micro/nanoscale, where these models and methods can be validated. On the other hand, MC brings about unique challenges resulting from its highly complex, nonlinear, time-varying channel properties due to the discrete nature of information carriers (molecules), substantial channel memory and peculiarities of molecular interactions at nanoscale, that cannot be always tackled by conventional ICT tools[4,12,13]. This leaves a huge question mark over the reliability of the existing MC methods, which are mostly adopted from conventional EM communications and not validated with practical testbeds.

Few studies in MC literature have focused on macroscale implementation of MC systems considering the physical limitations of a receiver, although the utilized receivers are made of off-the-shelf macroscale components. In one of these studies, the isopropyl alcohol (IPA) is used as information carrier, and commercially available metal oxide semiconductor alcohol sensors are used as MC receiver[14]. This study provides a testbed for MC with macroscale dimensions, which is utilised in a later study to estimate its combined channel and receiver model[15]. This testbed is then extended to a molecular multiple-input multiple-output (MIMO) system to improve the achievable data rate[16]. In a recent study, information is encoded in pH level of the transmitted fluid, and a pH probe sensor is used as the MC receiver[17]. In a later study, on the grounds that the use of acids and bases for information transmission can adversely affect the other processes in the application environment, such as in the human body, magnetic nanoparticles (MNs) are employed as information-carrying molecules in microfluidic channels[18]. In the same work, a bulky susceptometer is used to detect the concentration of MNs and decode the transmitted messages. The performance of MN-based MC, where an external magnetic field is employed to attract the MNs to a passive receiver, is analysed in another work[19]. However, the focus of these works is on macroscale MC using off-the-shelf sensors as receiver. Therefore, these studies do not contribute to the development of a design and optimisation framework for practical nanoscale MC receivers that can actually be integrated into micro/nanoscale devices.

As the first step to overcome this challenge, in this work, we report on the first implementation of a nanoscale MC receiver based on graphene field-effect transistor-based DNA biosensors (graphene bioFETs), and its ICT performance tests in a custom-designed microfluidic MC system. The main objective of this work is to provide an experimental testbed at physically relevant dimensions for nanonetworks, which can be used to reveal and study the effects of intricate biochemical and physical processes on the MC performance, and develop practical and realistic communication methods, including new MC detection techniques.

Graphene, with its exceptional electrical[20], chemical and mechanical properties[21], such as high carrier mobility (up to V s at room-temperature[22,23]), one atomic layer thickness and two-dimensional geometry exposing all its atoms to the sensing environment, provides very high sensitivity towards biochemical molecules especially in a bioFET configuration[24-27]. Owing to these properties, graphene has been extensively studied for selective sensing of a wide range of biomolecules ranging from carbohydrates[28] to proteins[29] and oligonucleotides[30,31]. Meeting the fundamental requirements of an MC receiver, such as the capability of label-free and reversible detection and high sensitivity[32], graphene bioFET stands as an ideal candidate for the implementation of the MC receiver. Flexibility and nanoscale 2d geometry of graphene are particularly favourable for the integration of graphene-based MC receiver into functional nanoscale devices.

Functionalisation of graphene with biomolecular probes can provide the selectivity against target analytes, required for avoiding biochemical interference for MC applications in physiologically relevant environments. In this work, graphene is functionalised with single-stranded DNA (ssDNA) probes (pDNAs) which undergo reversible hybridisation reaction with the complementary target DNAs (tDNAs). DNA hybridisation has found applications in a wide range of areas as a biotechnological tool to obtain selective and addressable artificial processes exploiting its high specificity, such as autonomous cargo transport[33], smart drug delivery[34], and fabrication of DNA origami-based stimuli-responsive nanomachines[35]. However, the main reason for selecting DNA as the recognition element is that DNA can be easily customized with different base sequences of different lengths, and can be designed to bind not only complementary DNAs but also peptides, proteins, carbohydrates and small molecules[36]. As a result, the implemented testbed can serve as a model system to provide insight into a broad range of MC systems relying on detection of molecular messages through affinity-based ligand-receptor interactions[6,12]. Moreover, integration of the fabricated MC receiver into a pressure-regulated microfluidic testbed provides control over the fluid flow rate, and enables flexibility and practicality in testing different channel geometries, which can mimic the most promising application environments of the MC inside human body, e.g., circulatory system[37,38].

In the remainder of this paper, we elaborate on the fabrication process of single layer graphene (SLG) bioFET-based MC receiver and its integration into a microfluidic testbed. The electrical characterisation of the device is performed at each step of functionalisation. Sensing response characteristics are revealed to determine the affinity between the complementary tDNA-pDNA pair. Selectivity of the device against complementary tDNAs is examined through real-time sensing response to non-complementary target DNAs (ntDNAs). Following the fabrication and sensitivity/selectivity analysis, we provide an MC detection performance analysis based on the transmission of pseudo-random binary data encoded into the concentration of tDNAs. The time-varying response of the MC receiver is fitted by a previously developed microfluidic MC model. This analysis provides important insights particularly into the infamous intersymbol interference (ISI) problem of MC. In the theoretical MC literature, ISI is typically attributed to the spatial dispersion of MC signals due to the low-pass characteristics of the diffusion or convection-diffusion channel and to the nondegradable information molecules remaining in the MC channel after their corresponding signalling interval, causing interference over the subsequent transmissions. In this paper, we experimentally show that the slow ligand-receptor binding reactions (i.e., DNA hybridisation) on the receiver surface result in a significant level of additional ISI, such that the previously transmitted information molecules remaining bound to receptors for longer than a bit interval cause a deviation from the expected received signal and complicate the detection of the subsequent bits. Similar to the existing approaches in the MC literature[39,40], a concentration difference-based detection method is utilised to cope with the ISI effects by obviating the need for channel state information (CSI).

Methods

MC receiver is fabricated in three consecutive steps. First, a graphene field-effect transistor (GFET) with chemical vapour deposition (CVD)-grown SLG is fabricated on Si/ substrate through optical lithography techniques. Then, a polydimethylsiloxane (PDMS)-based microfluidic channel is produced to encapsulate the GFET for bio-functionalisation, and real-time microfluidic sensing and communication experiments. Finally, for selectivity of the MC receiver against information-carrying target DNAs, bio-functionalisation of the GFET channels with probe DNA molecules is performed inside the microfluidic channel connected to a pressure-regulated microfluidic setup.

Fabrication of graphene field-effect transistor

MC receiver is fabricated with a CVD-grown SLG polycrystalline domain on an n-type Si/ substrate (525 μm with 90 nm thermal oxide layer). The CVD-grown SLG on Cu with PMMA coating (60 nm, 495K, A2) is obtained from Graphenea Inc.

Wet transfer of CVD graphene on Si/ substrate

The Cu layer of the PMMA/SLG/Cu stack should be removed before transfer onto Si/ substrate. The backside is partially covered with a graphene film, which may result in poor Cu etching performance. This backside graphene on Cu is etched away with plasma at 3 W for 30 s in a low-power Reactive Ion Etching (RIE) process. Following this step, Cu is etched by placing the PMMA/SLG/Cu stack on the surface of a solution of ammonium persulphate (APS) [1.8 g of APS in 150 ml DI water (18.2 M cm)].

The Si/ substrate is cleaned before the transfer by means of sonication for 10 mins in acetone followed by immersion in IPA for 5 mins and drying with nitrogen (). Once the Cu is entirely dissolved, the floating PMMA/SLG stack is transferred onto the surface of DI water in a beaker by a glass slide to clean the APS residuals, and then, the stack is fished onto the Si/ substrate. The resulting sample (PMMA/SLG//Si) is left vertically to dry overnight, and then annealed over a hot plate at 150 °C for 2 h. The sample is then transferred into a beaker with acetone for PMMA removal for 2 h, and then immersed in IPA for 5 mins and dried with , leaving only the SLG film on the Si/ substrate (Fig. 1a).

Figure 1

Process flow for the patterning of SLG channels. (a) SLG transferred on Si/ substrate. (b) Coating of sample with photoresist layer. (c) SLG channel pattern defined by optical lithography. (d) SLG channel pattern after RIE. (e) Removal of residual photoresist layer from SLG surface, and the resulting SLG channel on the substrate.

Patterning of SLG channels

MC receiver is designed to contain 7 GFETs having isolated source and drain contacts but being exposed to a common electrolyte gate. The individual SLG channels are patterned via optical lithography with a laser writer according to the design shown in Supplementary Fig. 4. Prior to all laser writing processes in this work, the sample is spin-coated with a photoresist (AZ-5214E from Microchemicals GmbH) at 4000 rpm for 60 s (Fig. 1b), and baked at 110 °C for 50 s on a hot plate. The photoresist layer is exposed by direct laser writing (wavelength-405 nm 169 mJ/) via a laser writer system, and the pattern is successively developed in diluted developer solution (1:4, AZ-351B/DI Water) for 35-45 s followed by brief immersion of the sample in DI-water for 2 s and drying with (Fig. 1c). The patterning of the individual SLG channels is completed with RIE removing the undesired areas of SLG film, which are not covered with the photoresist layer, via plasma at 3 W for 60 s (Fig. 1d). Finally, to remove the photoresist layer, the sample is successively immersed in acetone and IPA for 20 mins and 5 mins, respectively, and dried with , leaving the patterned SLG film on the substrate (Fig. 1e).

Deposition of contacts

In the following step, metal contact areas (source and drain) are defined on the sample through another optical lithography process (Fig. 2a,b). Depositions of 5 nm Cr and 50 nm Au are performed successively over the sample covered with the patterned photoresist layer by thermal evaporation at 10−6 mbar (Fig. 2c). Once the metals are deposited uniformly, the sample is dipped in acetone for 2 h for lift-off process, during which the metals over the photoresist layer are removed leaving only the patterned contacts (Fig. 2d).

Figure 2

Process flow for the deposition of contacts. (a) Coating of sample with photoresist layer. (b) Contact pattern defined by optical lithography. (c) Deposition of Cr and Au metal films through thermal evaporation. (d) Patterned contacts after lift-off process.

Deposition of Insulator

Finally, the drain and source contact areas, which might be exposed to the electrolyte during microfluidic experiments, are insulated to prevent any parasitic current between metal contacts through the electrolyte. For this, a thin layer of (20 nm) is uniformly deposited over the sample through atomic layer deposition (ALD) via Cambridge Nanotech Savannah at 120 °C with a base pressure of 0.5 mbar using trimethylaluminum as a precursor and O as oxidant. Oxidant and precursors are alternately introduced into the chamber in a 20 sccm flow of carrier gas (Fig. 3a). This step is followed by another optical lithography process, which defines the windows over to expose only the SLG channels to the electrolyte, and to expose the source and drain contact pads, which remain outside of the microfluidic channel for electrical measurements (Fig. 3b,c). The exposed areas of are then wet-etched in Phosphoric Acid (85% wt. in O obtained from Sigma-Aldrich) at 60 °C for 1 min (Fig. 3d). For removing any excess photoresist, the sample is dipped in acetone for 20 mins and IPA for 5 mins followed by drying with , leaving the patterned film on top (Fig. 3e). The optical images of the fabricated GFET are shown in Fig. 4.

Figure 3

Process flow for the deposition of insulator. (a) Uniform film over the sample after ALD. (b) Coating of sample with photoresist layer. (c) Insulator pattern defined by optical lithography. (d) Exposed SLG channel and contacts after wet etching of . (e) Patterned film after removal of excess photoresist.

Figure 4

(a) Optical micrograph of the fabricated GFET channels after etching process (only six of the seven channels are visible). (b) A closer look into one of the GFET channels. (c) Overall view of the fabricated 7-channel GFET before the bonding of microfluidic PDMS layer.

Fabrication of microfluidic channels and device integration

Fabricated GFET is encapsulated with a PDMS microfluidic channel, as demonstrated in Fig. 5. To this end, a 3d-printed mould is designed to define the geometry of the rectangular fluidic channel within the PDMS layer (see Supplementary Fig. 5). The microfluidic channel has a width of μm and a height of μm. At one end, the channel bifurcates for connection with the two channel inlets, which are designated for connection to the fluid reservoirs containing the buffer and information-carrying tDNA solutions during the communication experiments.

Figure 5

(a) Microfluidic PDMS layer bonded to the GFET surface after the inlets and outlet are defined, and the Pt gate electrode is placed on top. (b) Cross-sectional view of the MC receiver after PDMS layer bonding.

PDMS prepolymer is prepared using a 10:1 mixture of PDMS base monomer (Sylgard 184 Silicone Elastomer) and PDMS curing agent (obtained from Dow Corning Corporation). Air bubbles inside the PDMS are removed by degassing in a desiccator for 1 hour. The degassed mixture is poured onto the 3d-printed mould and left for curing overnight at room temperature. The cured PDMS is then carefully peeled off from its mould (see Supplementary Fig. 5). The inlet and outlet holes are punched through the PDMS layer for microfluidic connections by a biopsy punch (1.25 mm radius). A platinum (Pt) wire having a diameter of 0.5 mm acting as the common solution gate is then mounted to the top of the PDMS channel right above the SLG channels (Fig. 5, and Supplementary Fig. 5). The length of the Pt wire inside the channel is set to 1 cm.

In the next step, the patterned PDMS with the Pt solution gate is bonded to the surface of the MC receiver, ensuring that the graphene channels are well-aligned with the microfluidic channel and not placed under the PDMS walls. The most common method for bonding PDMS on and glass substrates is based on the plasma activation of the PDMS surface and the target substrate. This requires the surfaces of both the PDMS and the target substrate to be smooth. In our case, however, the exposed graphene channels on the target substrate prevents the application of the plasma, as this would cause the removal of graphene channels through plasma etching. Moreover, the plasma activation of only the PDMS surface is not sufficient because curing in the 3d-printed moulds made of Polylactic acid (PLA) results in PDMS layers with a rough surface (see Supplementary Fig. 5) rendering plasma activation ineffective in bonding. Therefore, we apply an alternative method, which was introduced in a previous work[41], for bonding porous membranes into PDMS devices. In this method, a thin layer of PDMS prepolymer, which is in liquid form, is coated on the bonding surface of the cured PDMS layer. Then, the PDMS is carefully placed on the sample, which is cleaned off any dust with prior to bonding. After placement, it takes approximately 1 min for the PDMS prepolymer to spread uniformly and cover the entire area between the PDMS and substrate except for the empty area defining the microfluidic channel. Once a uniform PDMS prepolymer layer is observed, the temperature of the hot plate is increased to 150 °C, and the prepolymer, serving as mortar, is quickly cured, resulting in a strong bonding. This method has consistently yielded leakage-free PDMS-substrate bonding during the fabrication process.

After bonding process, the inlet and outlet tubes are placed on the predefined inlet/outlet holes, as shown in Fig. 6b. Here, Teflon PTFE tubing (1/16” OD x 1/32” ID, obtained from Darwin Microfluidics) is preferred because of its higher chemical stability compared to Tygon tubing, which reacts with DMF used in the functionalisation process. The placement of the inlet and outlet tubes is followed by the application of PDMS prepolymer around the connection points of inlet, outlet and Pt gate over the cured PDMS layer for the complete sealing of the device.

Figure 6

(a) Microfluidic measurement setup consisting of a 4-channel pressure regulator, a high-precision SMU, electrical probes and microfluidic accessories. (b) A closer look into the fabricated graphene-based MC receiver connected to the microfluidic setup. (c) Probe connections for electrical tests of the device.

The dimensions of the microfluidic channel together with the fluid flow rate and fluid properties determine the Reynolds number, which is a dimensionless variable indicating the fluid flow regime in the channel[42]. Reynolds number is the ratio of the inertial forces to the viscous forces, and can be given bywhere is the fluid density, u is the linear flow velocity of the fluid, is the viscosity of the fluid, and is the hydraulic diameter, which can be obtained for rectangular channels as followsHere is the cross-sectional area of the channel, and is the cross-sectional channel perimeter. In the sensing and communication experiments of this work, water-based solutions are flowed at a constant volumetric flow rate μl/min. The linear flow velocity can then be obtained as μm s−1. By using kg and Pa s for water, we can obtain the Reynolds number for the microfluidic MC system as , indicating a strong laminar flow regime, where viscous forces overcome the inertial forces resulting in non-crossing, parallel streamlines[42].

Microfluidic setup

For functionalisation and electrical characterisation in the next steps, the device was connected to a microfluidic test setup, as shown in Figs. 6a and 7. The setup consists of a pressure regulator (OB1 MK3 - Microfluidic flow control system, obtained from Elveflow) with four pressure outlets, two of which are connected to fluid reservoirs through the pressure inlets. The fluid outlets of the fluid reservoirs are connected to the device through PTFE tubing, as shown in Fig. 6b and c .

Figure 7

Conceptual drawing of the microfluidic measurement setup with the practical implementation shown in Fig. 6.

Throughout the bio-functionalisation, sensing and communication experiments, microfluidic flow sensors are partly utilised for feedback-controlled modulation of the inlet pressure, and mechanical flow switches are used in cases where immediate stop/start of the microfluidic flow is required.

Functionalisation of GFET

Due to its one atomic thickness and 2d structure, the electronic properties of the pristine SLG is highly sensitive to the biochemical environment in the vicinity of its surface. Therefore, it suffers from low-level selectivity. On the other hand, in order to suppress the interference from other biochemical processes in physiologically relevant applications of the MC receiver, selectivity against information-carrying molecules is a must. Selectivity of the graphene can be realised through bio-functionalisation with recognition elements such as DNA and antibodies. As ssDNAs are preferred as target information-carrying molecules, i.e., tDNAs, in this work, the fabricated GFET is functionalised with probe DNAs (pDNAs), which are complementary to tDNAs.

For increasing the strength of the probe DNA immobilisation, and reducing the effect of nonspecific binding, the pristine SLG channels are first functionalised with 1-Pyrenebutyric acid N-hydroxysuccinimide ester (PBASE, obtained from Cambridge Bioscience Ltd), which has been widely utilised in the literature as linker molecules between graphene surface and DNA molecules[30,31,43,44]. PBASE is an aromatic molecule having an aromatic pyrenyl group and an amine-reactive succinimide group (Fig. 8a). PBASE exhibits a strong affinity towards SLG as its aromatic pyrenyl group interacts with the basal plane of graphene through π–π interactions resulting in a strong noncovalent binding (Fig. 8b). The noncovalent attachment of the PBASE does not alter the inherent electronic structure and physical properties of the graphene[45]. For the functionalisation of the SLG with PBASE molecules, 10 mM solution of PBASE in N,N-Dimethylformamide (DMF, anhydrous, 99.8%, obtained from Sigma-Aldrich) is prepared in a glass bottle, and sonicated for 30 s for mixing. The prepared PBASE/DMF solution is flowed through the microfluidic channel until the entire channel is filled with the solution. Then the flow is stopped, and the SLG channels are exposed to steady PBASE/DMF solution for 2 h. After functionalisation with PBASE, unbound PBASE molecules are removed from the channel with pure DMF, followed by rinsing with phosphate buffered saline (PBS, pH 7.4).

Figure 8

(a) Molecular structure of PBASE, and conceptual drawing of probe DNA (pDNA) and complementary target DNA (tDNA). (b) Noncovalent binding of PBASE to graphene via π–π interaction. (c) Immobilisation of pDNA via conjugation reaction with the succinimide group of PBASE. (d) pDNA-tDNA hybridisation.

The next step is the immobilisation of 18-mer 5’-amine-modified probe DNAs, which have the base sequence N-()-5’-AGG ACT TCA CCG TAT TGC-3’. The DNAs are custom designed and obtained from Sigma Aldrich. 2 μM of probe DNAs, prepared in PBS, is flowed through the microfluidic channel over the SLG channels. The device is left for immobilisation with probe DNAs overnight at 4 °C inside a wet chamber. The amine group of the pDNA reacts with the succinimide group of PBASE through conjugation reaction (Fig. 8c). The excess pDNA is then removed from the channel with PBS rinsing. Note that although the ssDNAs and the double-stranded DNAs (dsDNAs) shown in Fig. 8c and Fig. 8d are depicted as vertically aligned over the linker molecules, the orientation of DNAs tethered to surfaces through their single end can be influenced by the electrical potential of the surface, electrolyte flow conditions, the length of the DNAs, temperature, pH, ionic strength of the electrolyte, and the existence of the linker molecules. It is shown through molecular dynamics simulations that under zero potential of the surface and in the absence of the solution gate potential and the linker molecules, the flexible nature of the ssDNAs results in tilted and near-parallel orientation on the surface[46]. On the other hand, dsDNAs attain more vertical alignment under zero potential due to their higher rigidity[46]. However, to the best of author’s knowledge, the effect of the solution gate potential and the linker molecules has not been studied in the literature.

After pDNA immobilisation, the passivation of the unbound PBASE molecules is necessary to prevent nonspecific binding of target DNAs (tDNAs). This is performed by flowing 100 mM ethanolamine (OH) solution prepared in DI water through the microfluidic channel. Ethanolamine reacts with amine-reactive succinimide group of unbound PBASE molecules. With the passivation of PBASE, the device becomes ready for sensing and communication experiments with tDNAs.

Results and discussion

Electrical characterisation of MC receiver

For the electrical characterisation of the fabricated devices, direct-current (DC) measurements are taken using a high-precision source measure unit (SMU, Keysight B2902A), which is connected to the device electrodes via high-impedance passive probes, as shown in Fig. 6. On the other hand, the mobility of the GFET channels is measured before functionalisation in a back-gate configuration using a probe station, as detailed in Supplementary Information. Based on the linear approximation of the transfer curve (Supplementary Fig. 1), the mobility is calculated as V s−1.

Transfer characteristics

After each step of functionalisation, transfer characteristics of the devices are obtained with a constant drain-to-source bias mV, and a solution gate potential varying between – 0.2 and 1.2 V. The sweep rate of is set to 140 mV/s. All data are obtained after removal of excessive functional molecules from the microfluidic channel and the SLG surfaces by rinsing with PBS. This ensures that no change occurs in transfer characteristics due to ongoing chemical reactions. The PBS (pH 7.4) is used as the electrolyte in all measurements of transfer characteristics.

The measurements are taken from four of the SLG channels in the MC receiver, which have the same geometry, and provided in Fig. 9a–d. Hysteresis was negligible for all channels (see Supplementary Fig. 2), thus, only the forward sweep of is shown. First measurement is taken with only the PBS electrolyte inside the microfluidic channel prior to the functionalisation process. The p-type behaviour and ambipolar characteristics of the SLG-based devices are revealed with the charge neutrality point (CNP), i.e., the gate voltage of the minimum conductance, observed at 0.57 V on average over four channels with a standard deviation of 0.02 V. Upon functionalisation with PBASE linker, a negative shift of the CNP by mV is observed, indicating n-type doping. The negative shift of the CNP after the PBASE functionalisation was previously attributed to the dominance of the n-type doping effect of DMF in competition with the p-type doping effect of the PBASE molecules over long incubation times[45,47]. Note that the large standard deviation of the CNP shift is mostly resulting from the atypical behaviour observed in the third GFET channel, the transfer characteristics of which are demonstrated in Fig. 9c. In this channel, only 67 mV-shift in CNP is observed with the functionalisation of PBASE. The significantly smaller shift compared to other channels, which manifest consistent transfer characteristics, can be indicative of the poor functionalisation of the PBASE linkers on this particular GFET channel. The poor functionalisation can be due to the residual polymers or insulator material on the graphene surface remaining from the fabrication process and preventing the non-covalent attachment of the PBASE molecules.

Figure 9

Transfer characteristics of four channels at different steps of functionalisation in terms of drain-source current as a function of varying gate voltage with sweep rate 140 mV/s. In all measurements drain-source voltage is held constant at V.

On the other hand, the immobilisation of pDNAs resulted in a positive shift of the CNP in consistence with the previous literature[30]. DNA molecules are negatively charged at pH 7.4, attracting hole carriers to the graphene surface, thus, contributing to the p-type doping[30]. The shift of the CNP is observed as mV over the four channels. Again, the large standard deviation can be attributed to the third channel given in Fig. 9c, where the CNP shift is only 11 mV. This is again indicative of the poor functionalisation of the PBASE molecules, which in turn results in a very low concentration of immobilised pDNAs. Other channels, on the other hand, show similar CNP shifts, indicating more consistent immobilisation of DNAs.

The last set of measurements are taken after the passivation of excess PBASE linkers with the ethanolamine, and the introduction of the 0.01xPBS to the microfluidic channel to be used for the following sensing and communication experiments. 0.01xPBS is the 100-fold diluted version of PBS with DI water (18.2 M cm). While the ethanolamine does not possess any charge, the observed positive shift of the CNP is consistent with the previous literature reporting increased p-type doping with decreasing ionic concentration of the buffer solution[48]. Also note that in all of the measured GFET channels, a mobility reduction is observed upon passivation with ethanolamine in 0.01xPBS. The reason for switching to the diluted version of PBS for sensing and communication experiments is to decrease the effect of the Debye screening for enhancing the sensitivity of the device for the hybridisation events on the SLG surface. Note that in all measurements, leakage current has been detected to be under 15 nA (see Supplementary Fig. 3, and therefore, its effect on the transfer curve characteristics is negligible.

With the help of the transfer characteristics, we can deepen our analysis by determining the surface density of the immobilised pDNAs. To this end, we need to first determine the electrolyte gate capacitance, which can be approximated by the overall capacitance of three parallel plate capacitors connected in series:where is the EDLC between the graphene and the electrolyte, is the EDLC between Pt gate electrode and electrolyte, and is the quantum capacitance of graphene[30]. The EDLC of graphene to electrolyte can be calculated as , with m m being the area of graphene surface exposed to electrolyte, is the vacuum permittivity, and is the relative permittivity of PBS electrolyte, which is only slightly lower than the one of water, thus, taken as [30,49]. Lastly, is the Debye length which gives the thickness of the EDLC, and it can be approximated in aqueous solutions as in nm, with being the ionic density in M[50]. For 1xPBS buffer, the ionic density is 150 mM, thus, nm. The resulting EDLC for graphene is nF. The EDLC between the Pt electrode and the electrolyte can be obtained similarly as . However, since the surface area of the Pt wire inside the electrolyte ( cm mm μm2, calculated as the area of an half sphere of 1cm-length and 0.5mm-diameter) is significantly larger than the graphene surface area ( μm2), can be neglected. Lastly, the quantum capacitance of graphene per unit area has been reported as μF[30,51], which gives nF. The overall gate capacitance given by (3) then becomes nF.

The effective electric charge of a single immobilised 18-mer pDNA screened by the EDL can be written aswhere is the elementary charge, r is the effective length of the pDNA taken as the half of its length, i.e., (basepairs) (nm/basepair) nm, by assuming a vertical orientation for the single-stranded pDNAs following the analysis in line with the literature[30], where similar solution gate potentials and the same type of linker molecules are used.

Finally, the surface density of the immobilised pDNAs can be written as a function of the average shift in upon pDNA functionalisation:which gives . A similar surface density ( ) for pDNA was reported in a previous study[30].

Sensing response

For determining the sensing characteristics of the MC receiver, complementary 18-mer target DNA (tDNA: 5’-GCA ATA CGG TGA AGT CCT-3’, obtained from Sigma Aldrich) is prepared in 0.01xPBS solution. tDNAs of varying concentrations (50 nM, 100 nM, … 10 μM) are successively flowed through the microfluidic channel in the order of increasing concentration, and is recorded in real time with mV and V. During the experiment, the volumetric flow rate is held constant at μl/min. The results of the measurements taken from forth channel (whose transfer characteristics are given in Fig. 9d) are provided in Fig. 10a, where a decrease in the drain-source current is observed with increasing tDNA concentration, implying n-type doping effect in contrast to the p-type doping of pDNAs. The n-type doping effect upon target DNA hybridisation or probe DNA immobilisation was previously reported in recent studies[31,43,47,52-54], where the effect is mainly attributed to the partial interaction of the DNAs with the graphene surface through π–π stacking of the nucleobase aromatic rings resulting in direct electron transfer to graphene instead of electrostatic gating. The electron transfer from DNA upon immobilisation was also reported for CNT transistors[55]. Moreover, it is known that DNA molecules immobilised on a surface with their single end can be stretched in parallel to the surface under lateral flow[56], and the extent of this conformational change can be increased by a positive surface potential attracting the negatively charged DNA molecules to the surface[46,54]. Note that in our case, the graphene surface is continuously exposed to a positive drain-source bias during the sensing measurements. Therefore, we speculate that the conformational change resulting from microfluidic flow, hybridisation, and the positive surface charge of graphene brings the hybridised DNA molecules closer to the graphene surface, causing their partial interaction and electron transfer. Also, compared to the measurements taken with the single-stranded pDNAs under no-flow conditions, the ionic strength of the electrolyte (0.01 × PBS), in which the sensing experiments are performed, is significantly lower (compared to 1 × PBS), such that the attractive electrostatic force caused by the positive surface potential extends more into the electrolyte without being significantly screened[57]. Given that the hybridised dsDNAs carry twice the amount of negative charge of the ssDNAs, it can be considered that the hybridised pDNA-tDNA pairs, in our case, are more strongly attracted to the graphene surface compared to pDNAs[58]. The stronger electrostatic attraction combined with the stretching effect of lateral microfluidic flow supports our argument. However, this requires further confirmation, potentially through molecular dynamics simulations of both ssDNAs and hybridised dsDNAs under similar conditions to understand the effect of lateral flow, surface potential, solution-gate potential, and the linker molecules.

Figure 10

(a) Real-time sensing response of the MC receiver in terms of drain-source current with varying concentration of complementary target DNAs (tDNAs). (b) Equilibrium sensing response fitted by the Langmuir adsorption isotherm. Resulting dissociation constant for pDNA-tDNA hybridisation is nM.

Each working concentration of tDNAs were propagated in the channel until reaches a plateau. The value of at these plateaus are used to construct the sensing response graph of the MC receiver, which is provided in Fig. 10b. The response curve is fitted by the Langmuir adsorption isotherm, i.e.,where is the receiver response in saturation, which occurs when all the probe DNAs are hybridised. is the dissociation constant of pDNA-tDNA hybridisation, and is the applied concentration of tDNAs. The curve fitting gives the dissociation constant as nM for the DNA hybridisation on the fabricated MC receiver, and the receiver response at saturation as μA.

Table 1

Kinetic constants of DNA hybridisation measured by MC receiver.

DNA	Sequence	Binding rate k^+	Unbinding rate k^-	Dissociation constant K_D
		(M^-1s^-1)	(×10^-4 s^-1)	(×10^-6 M)
tDNA	5′-GCA ATA CGG TGA AGT CCT-3′	1814.9	13.538	0.746
ntDNA_1	5′-GCA ATA CGG CGA AGT CCT-3′	355.3	12.454	3.506
ntDNA_2	5′-GCA CGT CGG CGT CGT CAT-3′	48.9	13.110	26.829

Specificity

The specificity of the MC receiver against the complementary tDNAs is evaluated by comparing the receiver’s response to different ssDNAs, which are not complementary to the pDNAs. The ultimate specificity can be determined with the application of an ssDNA having only one single base-pair mismatch. For this purpose, we use 18-mer with the base sequence 5′-GCA ATA CGG CGA AGT CCT-3′, which has the mismatch in its 10th base pair, where T to C mutation occurs. Another test is performed with the application of 18-mer (5′-GCA CGT CGG CGT CGT CAT-3′), which has 7 base-pair mismatches. Complementary tDNA is also applied for comparison. All DNAs are dissolved in 0.01xPBS with 1 M working concentration. The measurement results before and after a moving mean filter of 21-s window length is applied are provided in Fig. 11a.

Figure 11

Specificity analysis of the MC receiver. (a) Real-time sensing response for complementary tDNA, non-complementary with single base-pair mismatch, and non-complementary with 7 base-pair mismatches (see Table 1). At s, the DNA solutions are replaced with 0.01xPBS solution to allow dissociation of the hybridised DNAs. (b) Real-time sensing response fitted by the Langmuir adsorption/desorption model, Eqs. (7) and (8).

The response curve of the MC receiver for different DNAs is fitted by the Langmuir model of adsorption to determine the kinetic rates of the DNA hybridisation for the three DNA sequences. The solution of the Langmuir model gives the time-varying response of the MC receiver during the association and dissociation phases of DNA hybridisation[30] as followswhere, is the input concentration, which is set to μM for all DNAs. denotes the time of dissociation, and is the asymptotic value of the sensing response, which occurs when the hybridisation reaches equilibrium. The variables to be fitted are the binding rate , unbinding rate , and . The fitted response is plotted in Fig. 11b, and the resulting kinetic rates are provided in Table 1, which shows that the binding rate of the target DNAs substantially decreases with increasing number of base-pair mismatches.

The nonlinear curve fitting gives the dissociation constant of the tDNA as nM, which is very close to the value obtained by the fitting of the sensor response, i.e., 730 nM. On the other hand, we obtain higher dissociation constants for non-complementary DNAs, i.e., μM and μM, indicating the specificity of the MC receiver against the complementary tDNAs.

Molecular communication performance

Time-varying communication experiments are performed in the microfluidic testbed to reveal the detection performance of the MC receiver. For these experiments, both inlets are utilised as shown in Fig. 6. One of the inlets is connected to the reservoir containing the 0.01xPBS buffer solution, and the other one is connected to the reservoir containing the tDNA solution in 0.01xPBS buffer. Manually controlled mechanical switches are utilised as they proved more effective in stopping the fluid flow into the channel than the digital control due to the fact that the pressure controller can drift out of calibration as the experiments progress. Accordingly, for the transmission of tDNAs, the buffer flow is rapidly stopped through the mechanical switch along the buffer line, and the tDNA line is opened at the same time. When the transmission ends, flow in the tDNA line is stopped, and the buffer flow is simultaneously started again by means of mechanical switching.

Time-varying response

To determine the time-varying response of the receiver, varying-length pulses of 1 μM tDNAs are flowed through the microfluidic channel. The results of three independent measurements taken from the same channel are provided in Fig. 12a and b for 30-second and 60-second pulses, respectively. As with the previous cases, a 21-second-length moving mean filter is applied for each measurement.

Figure 12

Normalised pulse response of the MC receiver in terms of drain-source current with constant operating voltages set to V and V. (a, b) Three independent measurements taken from the same channel for 30-second-long and 60-second-long 1 μM tDNA pulses, respectively. (c, d) Pulse responses fitted by the normalised output of the transformed model given in (13).

The time-varying response of the MC receiver to finite-length concentration pulses can be approximated by the analytical microfluidic MC pulse transmission model developed in our previous study[59]. This approximate model gives the receiver response in terms of number of bound receptors, i.e.,where and denote the start times of association (i.e., pDNA-tDNA hybridisation in this case) and dissociation phases, respectively, is the number of bound receptors at time t, is the number of bound receptors at equilibrium, is the number of bound receptors when the dissociation starts. Here denotes the principal branch of the Lambert function, and denotes Heaviside step function with . The parameters , and are provided in[[51], Eqs. (37–39)] as functions of , which is the upper limit of the receiver response in terms of bound number of receptors. In theory, can be taken as equal to the total number of receptors. Lastly, the transport parameter incorporating the effect of the microfluidic channel geometry and flow velocity on the molecular transport dynamics is also provided in [51, Eq. (13)] .

To use this model for fitting the electrical measurements obtained by the graphene-based MC receiver, we need to transform the variables in (9) to obtain a function of . The number of bound receptors is proportional to the change in the drain-source current through the following relation:with . Here, is the transconductance of the device, is the effective charge of a single tDNA molecule, and is the total gate capacitance in 0.01xPBS buffer. Gate capacitance is obtained as nF using (3) with the new Debye length in diluted PBS electrolyte, i.e., nm. The effective charge of a single tDNA is obtained similarly using (4) with . The transconductance of the device before the communication experiments can be calculated by a linear approximation of the curve around V bias. Averaging over the transfer curves of the four graphene channels obtained at the last step of functionalisation in 0.01xPBS, shown in Fig. 9, we obtain μA/V. Similarly, the following transformation is made:

Here, is the response of the receiver to 1 μM tDNA at equilibrium, and obtained from the sensing response given in Fig. 10 as μA. Note that can also be calculated as a function of as followswith being the average tDNA concentration passing over the receiver surface[59]. The dissociation constant for pDNA-tDNA pair given in Table 1. Accordingly, (9) can be rewritten by substituting (10), (11), and (12) into (9) as followswith the transformed parameters given as

The microfluidic MC pulse transmission model developed in our previous study[59] includes many other input parameters concerning the microfluidic channel geometry and molecular transport dynamics. The linear flow velocity is already obtained as μm/s by transforming the volumetric flow rate μl/min. Diffusion coefficient of tDNAs is taken as μm2/s, which is in the range of previously reported values for ssDNAs of similar lengths[60,61].

For comparison to the empirical MC signals, the model response is normalised by the mean empirical baseline current ( μA, as plotted in Fig. 10a), i.e.,

The normalised output of the transformed model for pulse lengths s and s is shown in Fig. 12c and d, respectively, presented in comparison to the mean of three independent measurements taken for each pulse length. The error bars represent the standard deviation of the measurements. We can observe that the approximate model is highly accurate in calculating the propagation delay and during the association phase. Although its accuracy is lower during the dissociation phase, the model well captures the pulse amplitude and the pulse width, which are the two important parameters in evaluating the performance of a communication system. The same model will be applied in the next section for binary data transmission.

Data transmission

To reveal the detection performance of the MC receiver, 20-bit-long pseudorandom binary information encoded into the concentration of tDNAs is transmitted through the microfluidic channel and detected by the fabricated MC receiver located at the bottom of the channel. Here bit-1 is represented by a 30-second pulse of 1 μM tDNA at the beginning of a bit interval, and bit-0 is represented by no pulse transmission during the entire bit interval, i.e., the buffer line stays connected into the microfluidic channel. The results are provided for bit intervals of 60-, 120-, and 360-s length in Fig. 13. As expected, the slow rate of dissociation of the bound tDNAs from the immobilised pDNAs results in a significant amount of ISI, which can potentially complicate the decoding at the receiver. The ISI is more pronounced for shorter bit intervals, i.e., higher transmission rates, such that the baseline could not find enough time to recover. On the other hand, for 360-second-long bit interval (Fig. 13c), ISI is significantly lower and is able to return to the baseline.

Figure 13

Normalised receiver response for binary data transmission with fixed pulse length s and varying bit intervals: (a) s, (b) s, (c) s. The grey lines denote the received MC signals normalised by the baseline current, and the solid red lines are their low-pass filtered version by a moving average filter of 21-s length in MATLAB. Solid blue lines represent the normalised output of the model given in (17). Dashed orange lines indicate the time instants when bit-1 is transmitted, i.e., when the mechanical switch in the tDNA line is opened for 30 s.

Another important observation is that the propagation delay between the transmission time of the bit-1 and the receiver response is 55 s in each test, and the delay remains almost constant during the entire data transmission period. This implies that the application of a constant delay shift in the receiver response can be sufficient for the synchronisation of the receiver with the transmitter during the decoding process.

As the receiver response is far from saturation in the performed tests, we can assume that the response is linear and time-invariant (LTI). Based on this assumption, we can reconstruct the received signal by applying the superposition principle of LTI systems in the approximate pulse transmission model developed in our previous study[59]. Accordingly, using (13), the received signal can be written aswhere is the number of transmitted bits, is the start time of transmission, is the bit interval, and is the ith transmitted bit. Similar to (17), the model response is normalised by the mean empirical baseline current ( μA), i.e.,The normalised model response is plotted in Fig. 13 over the empirical MC signal for varying bit intervals: s in Fig. 13a, s in Fig. 13b, s in Fig. 13c. We can see that the approximate model is very accurate in capturing the transmission delay and the overall trend of , although it is not able to exactly reconstruct the original signal. The deviations can largely be attributed to the LTI approximation in writing (18), which neglects the time-varying and nonlinear effects of the previously transmitted bits on the subsequent bits. In other words, the LTI approximation is based on the assumption that the kinetics of the ligand-receptor binding reactions on the receiver surface are not altered by the number of receptors bound to the ligands transmitted in the previous bit intervals. This assumption is made based on the observation that the number of bound receptors is much lower than the total number of receptors and the receiver is away from saturation, where the nonlinearity is maximized. On the other hand, we know that the rate of the forward reaction in ligand-receptor binding process always depends nonlinearly on the fraction of bound receptors, translating into a nonlinear response curve for the receiver as shown in Fig. 10b. Therefore, the deviation is mainly caused by the superposition of nonlinear responses within an LTI model. Other practical factors that cause approximation error are the nonspecific adsorption of tDNAs to the graphene surface and the steric hindrance effect on DNA hybridisation, the extent of which depends on the surface receptor concentration[62]. Nonuniform receptor distribution on the receiver surface and the effects of mixing at the microfluidic channel entrance, as well as the defects introduced randomly in the graphene channel during the receiver fabrication can also be speculated as important factors that might have resulted in the approximation error.

Based on the observation of constant delay during data transmission, the delay-shifted version of the bit intervals is indicated with dashed lines in Fig. 14 with the transmitted bits written inside each bit interval. As the envisioned MC applications demand low-complexity communication techniques due to the resource and size limitations of the communicating nanodevices, constant threshold detection has been favoured in the literature. In this scheme, the received signal is sampled at a predefined sampling time instant, and its amplitude is compared to a threshold value for deciding between bit-0 and bit-1. However, the high ISI of the MC channel, and the resulting drift of baseline, observed in these experiments, render the constant threshold detection methods ineffective, especially at high communication rates. Therefore, in line with some of the previous work in the MC literature concerning ligand-receptor binding systems, we utilise a simple difference-based detection method, which decodes the information based on the difference in measurements taken at the start and end points of a bit interval. Recall that the receiver is synchronised with the transmitter through the application of a constant delay shift at the receiver side. The sampling time points are demonstrated in Fig. 14 with green dots. Accordingly, the decoded bits can be written as followswhere is the decoded bit, and is the discrete time sample of at decision points, and is the indicator function, which outputs 1 if the inside expression is true.

Figure 14

Delay-shifted version of the receiver response for binary data transmission for varying bit intervals: (a) s, (b) s, (c) s. Dashed orange line denotes the start time of the data transmission, and grey dashed lines demarcate the individual bit intervals. Green dots on the received MC signal indicates the sampled current values for difference-based detection method with the decoding rule formulated in (20). Transmitted/decoded bits are noted above each bit interval, with the red-coloured ones denoting the erroneously decoded bits.

Applying the difference-based detection method on the unfiltered MC signal yielded 5% bit error rate (BER). We observed 1 bit error for 60 s bit interval, and 2 bit errors for 360 s bit intervals. The decoding of the unfiltered MC signal of 120 s bit interval yielded no error. The erroneous transmissions are indicated in Fig. 14. As can be observed in the figure, the one bit error for the case of 60 s bit interval was due to the ISI, such that the receiver output was not able to return to its baseline during bit-0 transmission, and information molecules nearby the receiver surface remaining from the previous bit-1 transmission might have caused a further decrease in the output current. For the case of 360 s bit interval, both bit errors occurred during bit-1 transmission. These were due to the relatively fast recovery of the receiver output during a long bit interval after a short (i.e., 30 s-length) MC concentration pulse was flowed over the receiver surface. Although these bit errors observed during the transmission of short bit sequences are not statistically significant, they demonstrate the effect of ISI and the importance of the optimisation of the bit intervals given the pulse length as well as the need for more reliable detection techniques. Nonetheless, the simple difference-based detection applied on the low-pass filtered MC signal correctly decoded all the transmitted bits.

Discussion

Data transmission rates in our experimental testbed are indeed quite low (in the range of 0.003–0.17 bit/s) revealing the effect of ISI as highlighted in our study. Higher data rates have been reported for MC in earlier studies which propose testbeds of different scales and technologies. For example, 0.25 bit/s data rate was reported using magnetic nanoparticles as information carriers in a fluidic channel[18]. In a macro scale testbed using the pH level for encoding information, the authors obtained 4 bits/s using a neural network-based detection method[17,63]. In a more recent experimental study, the authors reported data rates in the range of 5-10 bits/s with 0.2–5% BER, by encoding information into the concentration of NaCl molecules in tubular channels and employing a Viterbi decoder in the receiver[64]. However, as repeatedly stated in our manuscript, these testbeds are not scalable to micro/nanoscales and typically utilize computationally-complex detection techniques, e.g., neural network-based and Viterbi-based, which may not be suitable for the envisioned resource-limited micro/nanomachines.

Although specifying a range for required data rates for certain applications is quite speculative at this stage, given that coding, networking, and multiple access schemes are yet to be developed and tested for practical applications, we can discuss the two most promising applications of MC and IoNT, which are the smart drug delivery and continuous health monitoring with intrabody distributed molecular sensor and actuator networks[65]. The communicated information in these autonomous applications could be the presence of a target biomolecule informative of cancer progression or the feedback received from a smart drug carrier upon the successful unloading of a drug cargo to the target tissue. Communication of such information is unlikely to necessitate a continuous stream of large amount of data and likely to tolerate delay as in the case of low-throughput networks in IoT[66]. Additionally, considering that these applications are dealing with biological systems, which typically rely on protein-protein interactions, enzymatic reactions and passive means of mass transfer, e.g., diffusion, which are inherently slow processes, we do not anticipate high data-rate requirements in these applications.

Nonetheless, there are many rooms for improvement of data transmission rates. One direction is to develop high data-rate modulation schemes, such as ratio-based modulation[67] and depleted molecule-shift keying (D-MoSK)[68], and reliable detection techniques that can cope with the ISI, such as Viterbi decoders, matched filters, and data-driven machine learning-based detection methods, as already demonstrated to a certain extent in previous experimental MC testbeds. This unfortunately substantially increases the overall computational complexity, which may render such methods impractical for resource-limited micro/nano-machines.

The second direction, which is more promising to the authors’ opinion, is to optimise the physical design of the MC devices. There might be a low level of control over the propagation channels in the envisioned MC and IoNT applications; however, the MC receivers can be optimised through the rational design of ligand-receptor pairs. The nature has evolved to make the molecular interactions, e.g., DNA hybridisation, antibody-antigen binding, typically highly specific and strong, such that once these high-affinity ligand-receptor pairs bind to each other, they stay bound for very long periods of time[69]. Although such high-affinity interactions reduce the interference of biochemically similar molecules in the environment, they lead to significant levels of ISI and complicate the detection of subsequent symbols once utilised as part of the sensing mechanism in the MC receiver. On the other hand, having low-affinity ligand-receptor pairs for signalling would decrease the sensitivity of the receiver and may result in biomolecular interference in physiologically relevant environments which are typically crowded with different types of molecules. This tradeoff between specificity, sensitivity, and ISI can be addressed through optimising the binding affinity and kinetic reaction rates of ligand-receptor pairs to maximise the reliable data transmission rate. Aptamers, which are short nucleic acids that can be rationally designed to bind a broad range of biomolecules and ions, are ideally suited for this purpose, as their binding affinity and kinetic rates of reaction with target molecules can be tuned by modifying only their base sequences[70]. Optimisation and rational design of aptamers as receptors on a graphene-based MC receiver for varying types of information molecules along with the validation of the optimal communication performance in our microfluidic testbed are left for future work.

Conclusion

This proof-of-concept study reports on the very first fabrication and characterisation of a nanoscale graphene bioFET-based MC receiver. The ICT tests of the MC receiver is performed in a custom-designed micro/nanoscale MC system using a PDMS-based microfluidic channel as the controllable propagation medium, and a pressure-regulated flow control system for transmitting binary data encoded into the concentration of target DNA molecules. The time-varying hybridisation of the information-carrying target DNAs with the probe DNAs immobilised on the SLG surface of the MC receiver is selectively transduced into a change in drain-source current over the SLG channel. In light of the experimental sensing results, this transduction mechanism is speculated to be related to both electrostatic gating and direct electron transfer. The ICT performance of the MC receiver is evaluated by binary data transmissions at different bit intervals. The response of the MC receiver is well-fitted with the previously developed approximate analytical model of microfluidic MC. The slow hybridisation kinetics of DNA molecules is revealed to cause significant ISI. A simple difference-based detection method is shown to overcome the ISI to a significant extent, and provide reliable detection performance. The fabricated graphene-based nanoscale MC receiver and the overall microfluidic MC system can be used as an experimental testbed for probing intricate dynamics of MC, and developing novel communication techniques, transceiver architectures, and applications for MC and IoNT.

Supplementary Information.