1

Novel Modulation Techniques using Isomers as
Messenger Molecules for Molecular
Communication via Diffusion

arXiv:1201.0913v1 [q-bio.QM] 3 Jan 2012

Na-Rae Kim and Chan-Byoung Chae
School of Integrated Technology
Yonsei University, Korea
Email: {nrkim, cbchae}@yonsei.ac.kr

Abstract—In this paper, we propose novel modulation techniques using isomers as messenger molecules for nano communication via diffusion. To evaluate achievable rate performance, we
compare the proposed techniques with concentration-based and
molecular-type-based methods. Analytical and numerical results
confirm that the proposed modulation techniques achieve higher
data transmission rate performance than conventional insulin
based concepts.
Index Terms—Nano communication network, molecular communication, modulation technique, isomer, diffusion, messenger
molecule.

I. I NTRODUCTION
Since R. Feynman’s talk in 1959, people have eagerly
pursued practical work on smaller and smaller scales, notably
through nanotechnology [1]. Nanotechnology, recently, produced a new branch of research called nano communication
networks (NCNs). NCNs interconnect several nanoscale machines (nanomachines in short) to carry out more complex
tasks in a cooperative manner, or to perform simple tasks [2].
These networks are not just smaller versions of traditional
communication networks; they have their own features and
are applicable in many fields, including biomedical, industrial,
military, and environmental [3].
NCNs can be realized by several methods. We can rely on
traditional communication systems that use electromagnetic
fields or ultrasonic wave. We have to overcome, however,
some radio frequency (RF) device barriers [2]. On the practical
front, researchers are considering new materials, such as
carbon-nano tubes (CNTs) and graphene [4]. Little research,
however, has focused on channel models and human body
absorption with terahertz bands. In addition to theses activities,
a new concept utilizing diffusion especially for short range
communication has been introduced [5].
This new paradigm of communication, molecular communication, sends/receives information-encoded molecules between
nanomachines. One of its advantages is its biocompatibility.
Biocompatibility is the most important and challenging issue
confronting intra-body applications. By using biomolecules
we can create inherently biocompatible systems that require
This work was in part supported by the Ministry of Knowledge Economy
under the “IT Consilience Creative Program” (NIPA-2010-C1515-1001-0001)
and the Yonsei University Research Fund of 2011.

no inorganic harmful materials. Moreover, driven by chemical
reactions, molecular communication is energy efficient. For
these reasons, molecular communication is the main focus of
this paper.
The authors in [6], [7] have extensively studied the fundamentals of molecular communication via diffusion. A new
energy model has been investigated to understand how much
energy is required to generate messenger molecules in [6]
and concentration-based and molecular-type-based modulation
techniques have been introduced in [7]. By using a simple
symmetric binary channel model, the achievable transmission
rate (achievable rate hereafter) has been extensively compared.
In analyzing their modulation techniques, however, [6], [7]
did not clearly suggest concrete structures for the messenger
molecules. They considered insulin-based nano networks but
it is unclear how these could be utilized for molecular-typebased modulations. In this paper, to maximize the achievable
rate with less transmit power/energy, we propose using isomers
as messenger molecules. To consider the properties of isomers,
we slightly modify the energy model described in [6]. To
the best of our knowledge, the proposed method is the first
attempt to design appropriate messenger molecules for NCN
via diffusion.
This paper is organized as follows. Section II describes the
channel and energy models under consideration. Section III
explains isomers for messenger molecules and also proposes
two modulation techniques, i) concentration-based and ii)
molecular-type-based. We present the numerical results in
Section IV, and Section V offers our conclusions.

II. S YSTEM M ODEL
For simplicity, we consider a nano communication system
consisting of a single transmitter and a single receiver, as
illustrated in Fig. 1. For this system, the messenger molecules
diffuse through the medium1 (i.e. water) at body temperature.
This paper assumes that no collisions occur among the propagating messenger molecules and/or the molecules and the
medium.
1 In our future work, blood will be considered as the medium to get more
reasonable results.

2

total number of molecules during a symbol duration becomes
different depending on the current and previous symbol.
Binomial(n, Phit (d, Ts )),
Nc ∼ N (np1 , np1 (1 − p1 )),
Np ∼ N np2 , np2 (1 − p2 )) − N (np1 , np1 (1 − p1 ))

Fig. 1. System model of nano communication via diffusion with a single
transmitter and a single receiver.

where, Nc denotes the number of molecules transmitted and
received during the current symbol duration, and Np is the
number of molecules transmitted from the previous symbol
duration but received during the current symbol duration. For
simplicity, we use p1 for Phit during Ts and p2 for Phit during
2Ts .

A. Channel Model
The particles or messenger molecules released from a transmitter nanomachine spread out through the medium by Brownian motion [6]. Such motion is basically driven by diffusion,
meaning the particles move from areas of higher concentration
region to areas of lower concentration. A random process,
the displacement of messenger molecules follows a normal
distribution with zero mean. Thus, the standard deviation σ of
the displacement can be obtained as follows:
4 X ∼ N (0, σ 2 ),
∂2c
∂c
= D 2,
∂t
∂x
2
1
e−x /4Dt ,
c(x, t) =
1/2
(4πDt)
x2 = 2Dt,
√
σ = 2Dt,
Kb T
D=
6πηrmm

(1)
(2)

(3)

where, c indicates the concentration of Brownian particles at
time t at point x, and D represents the diffusion coefficient
of the particles calculated from the Boltzmann constant (Kb ),
temperature (T ), the viscosity of the medium (η) and the
radius of a messenger molecule (rmm ). Eq. (1) is Fick’s
second law of diffusion. By solving this partial differential
equation, we obtain the general equation for the concentration
of the particles represented in (2). The first moment of (2)
is zero, which indicates the displacement has a zero mean,
and the second moment becomes the variance.
Hence, the
√
displacement has a standard deviation of 2Dt, and finally
we have
4X ∼ N (0, 2Dt).

(4)

When n messenger molecules are transmitted by the transmitter nanomachine, the molecules have a probability of hitting the receiver nanomachine. We represent this as a binomial
distribution of n times of trials with a probability of Phit for
each. Phit is determined by the symbol duration (Ts ) and the
distance between the transmitter and the receiver (d), which are
both affected by the diffusion coefficient. If n is large enough,
and nPhit is not zero, binomial can be approximated as a
normal distribution. In addition, we have to take into account
the overflow molecules from the previous symbol [6]. The

B. Energy Model
We use an energy model to that described in [6]. To prevent
them from interacting with others during propagation, the
messenger molecules are encapsulated in vesicles. The vesicles
are then carried to the boundary of the machine (e.g. the
eukaryotic cell), and released into the propagation medium.
Here we use runit to denote the machine radius size.
The energy needed in these steps are calculated. Moreover,
we can obtain the cost of synthesizing messenger molecules
(hexoses as an example) by the enthalpy of formation ∆H.
Then the final energy model can then be presented as follows [6]:
n
ET = nES + (EV + EC + EE ),
cv
∆Hhexose
J per messenger molecule,
ES =
6.02 × 1023
EV = 83 × 5(4πrv 2 )zJ,


runit /2
zJ,
EC = 83 ×
8
EE = 83 × 10zJ,

3
rv
√
cv =
rmm 3
where, ET is the total energy cost required to transmit n number of molecules, ES is the synthesizing cost of one hexose
molecule calculated from the sum of bond energy (e.g., the
enthalpy change), EV is the vesicle-synthesizing cost having
a radius of rv . EC is the cost of intra-cellular transportation,
EE is for membrane fusion, and cv is the capacity of one
vesicle that is related to the radius of messenger molecules
rmm .
III. M ODULATION T ECHNIQUES
In molecular communication via diffusion, the unique properties of the messenger molecules can determine modulation
techniques. Two such techniques that have been proposed
include the use of concentration and type of messenger
molecules [7]. A usable messenger molecule itself, however,
has yet to be proposed. Thus this paper proposes practical
messenger molecules, as well as analyzes and compares their
achievable rates by applying different modulation techniques.

3

H

C

O

H

O

C

O

H

C

O

C

O

H

C

O

H

C

O

H

C

O

C

OH

HO

C

H

H

C

OH

HO

C

H

H

C

OH

H

C

OH

C

H

HO

C

H

H

C

OH

H

C

OH

HO

C

H

HO

C

H

H

H

C

OH

H

C OH

H

C

OH

H

C

OH

HO

C

H

HO

C

H

OH

H

C

OH

H

C

H

C

OH

H

C

OH

H

C

OH

H

C

OH

OH

HO

C

H

C

OH

H

C

OH

HO

C

H

HO

C

H

C

OH

H

C

CH2OH

CH2OH

D-Altrose

D-Glucose

OH

CH2OH
D-Mannose

CH2OH
D-Gulose

CH2OH
D-Idose

CH2OH
D-Galactose

CH2OH
D-Talose

D-form isomers of aldohexoses.

Mirror

H

C

O

H

C

O

HO

C

H

H

C

OH

OH

HO

C

H

OH

HO

C

H

H

C

OH

HO

C

H

H

C

H

C

CH2OH
D-Glucose

Fig. 3.

H

H

C

CH2OH

H

HO

H
H

D-Allose

Fig. 2.

C

CH2OH
L-Glucose

Mirror image of D- and L- glucose. They are enantiomers.

A. Isomers for Messenger Molecules
The most important thing we have to be careful when
designing messenger molecules is that has to be non-toxic
to human body. The messenger molecule suggested in [6],
however, is highly flammable (i.e., hydrofluorocarbon), which
means it may be inappropriate as a messenger molecule.
For several reasons, potential candidates for messenger
molecules are isomers, molecules with the same number and
types of atoms [8]. First of all, they consist of the same type
of atoms, lightening the burden on the transmitter nanomachine synthesizing the messengers. For numerical analysis,
this paper uses the isomers known as hexoses, especially
aldohexoses. Note that we can select from the aldoses family,
hexoses, pentoses, tetroses, or trioses based on the required
modulation order.
Aldohexoses are monosaccharides with the chemical formula C6 H12 O6 . They have four chiral carbon atoms that give
them 16 (=24 ) stereoisomers. Fig. 2 illustrates eight kinds of
D- form diastereomers. The enantiomers of each molecule
are another set of eight L- form diastereomers. Therefore,
aldohexoses have 16 different shapes. Here, diastereomers are
stereoisomers that are not enantiomers (mirror-image isomers).
For example, D-glucose and L-glucose as shown in Fig. 3
are enantiomers. D-glucose and D-galactose, isomers but not
mirror images, are diastereomers.
When each isomer is dissolved in aqueous solution, it
mostly exists as a cyclic form. D-glucose, for instance,
undergoes nucleophilic addition reaction generating four
cyclic anomers, α-, β- forms of D- glucopyranose and Dglucofuranose [9]. We consider, however, only two pyranose
forms since they predominate with 36 and 64 percentages. If
the functional groups attached to carbon number 1 (C1) and
the C5 shown in Fig. 4 have a trans-structure, it is called

α- form, and if a cis-structure, β- form. They interconvert
each other in solution through a process called mutarotation.
Therefore, by deploying hexoses group in the system, there
are a total of 32 different isomers. This means that we can
increase a modulation order up to 32, i.e., 5 bits per symbol.
B. Molecular Concentration Based
When the concentration of messenger molecules is used,
the technique is known as concentration-based modulation,
originally introduced in [7]. It is also referred to as CSK
(concentration shift keying). In binary CSK (BCSK), one
threshold exists, and a receiver nanomachine decodes the
symbol as ‘1’ if the number of received messenger molecules
exceeds the threshold, ‘0’ otherwise [6]. Generally, 2n -CSK
systems transmit n bits per symbol, and requires 2n -1 number
of thresholds. In this system, there is, theoretically, no limit
in the modulation order. As the modulation order increases,
however, so does the error probability since the minimum
distance between two neighboring thresholds decreases. It only
uses one kind of molecule, and D-glucopyranose is used for
analysis here. Below is a probabilistic analysis for BCSK
systems:
Pa (0, 0) = Pb (0, 0, 0) + Pb (1, 0, 0)
1
= [P (Nn < τ ) + P (Np + Nn < τ )]
4"
!#
τ 
1
τ − n(p2 − p1 )
=
−Q p
2−Q
,
4
σ
n[p2 (1 − p2 ) + p1 (1 − p1 )] + σ 2
Pa (0, 1) = Pb (0, 0, 1) + Pb (1, 0, 1)
1
= [P (Nn ≥ τ ) + P (Np + Nn ≥ τ )]
4"
!#
τ 
1
τ − n(p2 − p1 )
=
Q
+Q p
4
σ
n[p2 (1 − p2 ) + p1 (1 − p1 )] + σ 2
where, Pa (X, Y ) indicates the probability of X sent and Y
received, and Pb (Z, X, Y ) is the probability of X sent, Y
received, and Z previously sent. Q(·) is the tail probability of
the normal distribution.
C. Molecular Type Based
When different types of molecules represent different symbols, the technique is known as molecular-type-based modulation, referred to as molecule shift keying (MoSK). Unlike

4

H

1

C
2

H

C

HO

C

H

C

H

C

3

4
5

TABLE I
S IMULATION PARAMETERS .

O
OH

Parameters
Phit for Ts
Phit for 2Ts
Radius of the hexoses [12]
D of hexoses (3)
4Hhexose
Ts [6]
Viscosity of the water
Temperature (body temperature)

H
OH
OH

6

CH2OH
OH
4

α

6

O

5

HO
HO

2
3

OH

1

OH

OH
4

HO
HO

β

6

O

5
2
3

Values
0.6097
0.7208
0.38 nm
597.25 µm2 /sec
1271 kJ/mol
5.9sec
0.001kg/sec · m
36.5 ◦ C = 310K

1

OH

OH

Fig. 4.
D-glucose and its α- and β- anomers. They undergo rapid
interconversion in water, which is called mutarotation.

the work in [7], we here use a set of isomers. We thus
name it IMoSK (isomer-based MoSK). IMoSK requires only
one threshold, making it simpler than the CSK system. For
systematic analysis, we can choose one among the several
aforementioned isomer sets; a modulation order can be determined by the sets used. For example, hexoses have a
modulation order of up to 32, trioses of 4. For simplicity,
we apply, as did in [7], the addictive white Gaussian noise
(AWGN) model. The mutarotation effect is also considered.

error term. Thus, Pa (α, β) can be calculated as
Rtα − Req
0.99
=−
ts + 1,
Rα − Req
3600
αRα + βRβ
αRα + (n − α)Rβ
Rtα =
=
,
α+β
n
(Rtα − Rβ )n
α=
, β = n − α.
Rα − Rβ
β
If β ≥ τ, Pa (α, β) = Pa (α, β)AW GN +
n
where, Req = 52.7◦ , Rα =112.2◦ , and Rβ =18.7◦ [10]. The
parameters α and β indicate the number of α- form and βform molecules, respectively, and n is the number of total
molecules transmitted. Optical specific rotation of the chemical
compound is defined as the observed angle of optical rotation
when plane-polarized light passes through D-glucopyranoses.2

1) B-IMoSK-AWGN: Only AWGN considered,
Pa (α, α)AW GN = Pb (α, α, α) + Pb (β, α, α)
1
= [P (Np + Nc + Nn ≥ τ )P (Nn < τ )
4
+ P (Nc + Nn ≥ τ )P (Np + Nn < τ )],
Pa (α, β)AW GN = Pb (α, α, β) + Pb (β, α, β)
1
= [P (Nn ≥ τ )P (Np + Nc + Nn < τ )
4
+ P (Np + Nn ≥ τ )P (Nc + Nn < τ )].
Note that Pa (β, β) and Pa (β, α) can also be calculated
similarly. Due to the page limit, we omit it here.
2) B-IMoSK-Mutarotation Considered (α-D-glucopyranose
⇔ β-D-glucopyranose): When α- and β- D-glucopyranose
are chosen, there is a possibility of incorrect decoding, such
as α- sent, β- received, or β- sent, α- received due to
mutarotation. Thus we derive an error probability considering
the mutarotation process. The α- or β- form sent varies its
number with time [10], and the number can be calculated by
observing the specific optical rotation. Rtα or Rtβ (observed
specific optical rotation at time t) minus Req (at equilibrium)
divided by Rα or Rβ (at time 0) minus Req has a linear
relationship with time (Ts ) at about 36.5◦ C. Rtα and Rtβ
values are calculated by Rα and Rβ , and Req cab be found
in [10]. The number of α- and β form existed after time Ts
is obtained as below. If the number of β- form exceeds the
threshold value after Ts when α- sent, nβ value is added to the

3) 32-IMoSK: When using hexoses, the system has the
maximum modulation order of 32. Probabilities of X sent and
Y received for all X and Y values are obtained similarly and
thus omitted here (see [13]).
IV. N UMERICAL R ESULTS
Assume that all hexoses have the same physical properties;
size, diffusion coefficient, and enthalpy of formation. Phit
value is calculated by the same numerical calculation used
in [6]. We approximate the value for the hexoses by assuming
it to be proportional to the diffusion coefficient. Therefore,
we have the hitting probabilities shown in Table I and apply
these into the proposed scheme explained in Section III. Here,
we define the achievable rate R that maximizes the mutual
information I(X; Y ) as follows
X
P (X, Y )
I(X; Y ) =
P (X, Y ) log2
,
P (X)P (Y )
(5)
R = max I(X; Y )
τ

where, P (X, Y ) denotes the joint probability of X and Y
while P (X) and P (Y ) are the probabilities of events X and
Y . In (5), the threshold values vary from 1 to 1000.
Fig. 5 compares the achievable rate of the proposed method
using hexoses with the conventional insulin based method in
BCSK and B-IMoSK systems. It is remarkable that we obtain
about 8 dB of SNR gain in both BCSK and B-IMoSK systems.
2 It could be measured by a polarimeter, and there is a linear relationship
between the observed rotation and the concentration of the compound [11]

5

Fig. 5. Achievable rate comparisons of the conventional insulin based and
the proposed methods using two kinds of aldohexose isomers.

Fig. 7.

Achievable rate comparisons. Triose and hexose used.

messenger molecules. We first introduced energy and channel
models for our system. Next we proposed several modulation
methods able to support up to 5 bits per symbol. We also compared the achievable rate performance with existing modulation methods (concentration-based and molecular-type-based).
This work differs from prior work in that it proposes practical
messenger molecules and provides guidelines for selecting
from among several possible candidates. Future work, we will
consider, to achieve more modulation degrees of freedom,
the ratio of enatiomers. To increase the transmission data
rate further, we will also consider multiple sets of messenger
molecules.
R EFERENCES
Fig. 6. Achievable rates of 32-IMoSK and B-IMoSK. Thirty two kinds of
isomers of aldohexose used for 32-IMoSK.

In addition, the B-IMoSK system shows a better data rate
performance than BCSK. This is due mainly to the size of the
proposed messenger molecule being much smaller than that
of insulin. Hence, the transmit energy of the proposed method
is much less than that of the conventional method. This result
will be discussed in more detail in [13].
Fig. 6 shows the achievable rates of 32-IMoSK and BIMoSK using hexoses as messenger molecules. As can be
seen from Fig. 6, the 32-IMoSK system has the maximum
rate of 5 (bits per symbol), and the B-IMoSK of 1 (bit per
symbol). In Fig. 7, we compare the achievable rates using
trioses and hexoses. Obviously, trioses achieve higher SNR
gain than hexoses due to their smaller sizes. Trioses, however,
have the data transmission limit 1 bit per symbol. From this
result, we can conclude that trioses can be selected for a low
data rate system with a higher transmission reliability and
hexoses for a high data rate system.
V. C ONCLUSION
To make nano communication feasible in practice, this
work proposed novel modulation techniques using isomers as

[1] R. P. Feynman, “There’s plenty of room at the bottom (data storage),”
Journal of Microelectromechanical Systems, vol. 1, pp. 60–66, 3 1992.
[2] I. F. Akyildiz, F. Brunetti, and C. Blzquez, “Nanonetworks: A new
communication paradigm,” Computer Networks, vol. 52, pp. 2260 –
2279, Apr. 2008.
[3] L. P. Gin and I. F. Akyildiz, “Molecular communication options for long
range nanonetworks,” Computer Networks, 2009.
[4] B. Atakan and O. B. Akan, “Carbon nanotube-based nanoscale ad hoc
networks,” IEEE Comm. Mag., vol. 48, pp. 129–135, 6 2010.
[5] M. Moore, A. Enomoto, T. Nakano, R. Egashira, T. Suda, A. Kayasuga,
H. Kojima, H. Sakakibara, and K. Oiwa, “A design of a molecular
communication system for nanomachines using molecular motors,”
Proc. IEEE Int. Conf. Pervas. Comput. Commun. Workshops, 2006.
[6] M. S. Kuran, H. B. Yilmaz, T. Tugcu, and B. zerman, “Energy model for
communication via diffusion in nanonetworks,” Nano Communication
Networks, vol. 1, pp. 86 – 95, July. 2010.
[7] M. S. Kuran, H. B. Yilmaz, T. Tugcu, and I. F. Akyildiz, “Modulation
techniques for communication via diffusion in nanonetworks,” Proc.
IEEE Int. Conf. on Comm., 2011.
[8] T. McKee and J. R. Mcfee, Biochemisty. Oxford University Press,
2009.
[9] P. M. Dewick, Essentials of Organic Chemistry: For Students of Pharmacy, Medicinal Chemistry, and Biological Chemistry. Wiley, 2006.
[10] A. Carruthers and D. L. Melchior, “Transport of α- and β-D-glucose
by the intact human red cell,” Biochemistry, vol. 24, pp. 4244–4250,
1985.
[11] J. R. Mohrig, C. N. Hammond, and P. F. Schatz, Techniques in Organic
Chemistry. Freeman and Company, 2010.
[12] S. G. Schultz and A. K. Solomon, “Determination of the effective
hydrodynamic radii of small molecules by viscometry,” J. Gen. Physiol.,
1961.
[13] N.-R. Kim and C.-B. Chae, “Novel modulation techniques using isomers
as messenger molecules for molecular communication via diffusion,” in
preparation for journal publication.

