Optical Expanders with
Applications in Optical Computing
John H. Reif Akitoshi Yoshida
July 20, 1999
Abstract
We describ e and investigate an optical system whichwe call an optical expander. An optical expander elec-
tro optically expands an optical b o olean pattern enco ded in d bits into an optical pattern of size N bits. Each
expanded pattern is one of the N mutually orthogonal b o olean patterns. We wish the expansion to b e exp onen-
tial, so wehaved=clog N for some constant c. An optical expander can b e viewed as either an electro optical
line deco der which converts d bits of optically enco ded binary information to up to N unique optical outputs;
or a digital b eam de ector which de ects an input laser b eam into one of N distinct directions with a control
signal of d bits. We show that an optical expander can not b e constructed by using linear optical systems, so a
non-linear optical lter must b e used. We describ e two di erent architectures to implement an optical expander.
One uses an optical matrix-vector multiplier and an arrayofN threshold devices. The other uses log N novel
re ection/transmission switching cells. We then analyze these architectures in terms of size, energy requirement,
and sp eed.
Our optical expander can help develop various applications in electro optical computing systems. In general,
b ecause of I/O constraints and the limited fan-in/fan-out of electrical circuits, the conventional VLSI technology
is not suitable for building a large line deco der. On the other hand, conventional acousto optic b eam de ectors are
bulky and limited by capacity-sp eed pro duct. We show that as an electro optical functional unit, the line deco der
nds many applications in optical computing such as optical interconnects and optical memory.Thus, the design
and development of optical expanders is vital to optical computing. Our optical expander utilizes high sp eed and
high space-bandwidth pro duct connections provided by optical b eams in a volume, so it o ers fast and accurate
op erations.
Holographic memory system and message routing systems are p otential applications of the optical expanders,
and are discussed to further motivate the design and development of optical expanders. Key words: Optical
computing, electro optical interconnections, holographic memory.
1 Intro duction
1.1 Potential of Optical Computing
Optical Computing has recently b ecome a very active research eld. Optics has b een used for image pro cessing
and long distance communications as well as for lo cal area networks. Recently, much attention was given to the
incorp oration of optics into VLSI electrical circuits. In uenced by the success of optical lo cal area network and the
progress of optical emitters and detectors, there has b een a large amount of attention fo cused on digital optical
computing utilizing optics to o er global interconnections with large fan-outs. The p ossibilities of this approach are
also contrasted to the limitation of current VLSI technology. The VLSI technology is not suitable for interconnection
intensive circuits due to its two dimensionality, I/O constraints on the chip b order, and electrical prop erties such
as resistance, capacitance, and inductance. In contrast, optics can utilize free-space interconnects as well as guided
wave technology, neither of which has the problems mentioned for the VLSI technology. Many researchers have b een
investigating suitable optical logic devices, interconnection schemes, and architectures for such optical computing
systems. Surveys and intro ductions are found in [1] [2] [3] [4] [5].
Department of Computer Science, Duke University, Durham, NC 27706. 919-660-6568 919-660-6561 Email: [email protected] and
[email protected]. Research supp orted in part by Air Force Contract No. AFOSR-87-0386, Oce of Naval Research, Contract No. N00014-
87-K0310, DARPA/ARO Contract No. DAAL03-88-K-0195, DARPA/ISTO Contract No. N00014-88-K-0458 1
In general, replacing electrical devices and interconnects with their functional equivalent in optics do es not lead
to a go o d result. Wemust nd a particular area where optics can b e advantageous over electronics. One such area
is in high sp eed interconnection requiring large fan-outs among pro cessing units, mo dules, or b oards. Because optics
can provide free interconnections in a volume without intro ducing various drawbacks found in the two dimensional
electrical wiring, it can provide an alternativeway to implement an interconnection network from a large number
of source units to a large number of destination units. Another area in which optics may play an advantageous
role is in holographic memory. Holographic memory may be advantageous over the conventional random access
memory RAM in terms of storage capacity, nonvolatility, and parallel readout capability. However, in spite of their
p otential, practical implementation has b een limited by the lack of electro optical interface which creates a mutually
orthogonal optical pattern of size N from an input of size log N . We call this interface the optical expander. We
discuss the optical expander and its applications in detail.
1.2 Disadvantage of VLSI implementation
Optical Expander can b e viewed as an electro optical implementation of a line deco der. The line deco der is a basic
comp onentofcombinatory logic, and it is used to select one of N devices with an input control signal of size log N
bits. The typical VLSI implementation of a large line deco der su ers from two problems. One is due to the top ological
prop erties of a VLSI chip. In VLSI, all wires must run on a two dimensional plane with a constantnumberoflayers,
and all I/O pads must be lo cated on the b order of the chip. As N b ecomes large, even though the logic gates
o ccupies a small area, the interconnections and the I/O o ccupy a large area, resulting in a serious area consumption
problem. The other problem is due to the small fan-in and fan-out of electrical devices. Atypical implementation
requires a tree structure of logic gates in log N stages. This will lead to log N time steps to generate an output.
A sp eedup maybe obtained by taking advantage of the capacitance of electrical wires. One can use a two-phase
clo ck with a bus con guration. In phase one, N buses are charged with Vdd. Each bus has its corresp onding log N
pulldown gates to discharge the voltage. In phase two, the pulldown gates are enabled to discharge the buses, and the
outputs are inverted to generate N outputs. However, this will still not solve the I/O constraints problem mentioned
earlier. These I/O constraints may force the chip to transmit N bits in bit sequential manner. Optics may provide
an alternativeway to implement a large line deco der by using its free space three dimensional interconnection with
large fan-in and fan-out capabilities.
1.3 Disadvantage of Beam De ector
Analog b eam de ectors based on acousto optic e ect have several drawbacks. First of all, they are bulky and
acousto optic mo dulators require high drive power. Secondly, they are limited by capacity-sp eed pro duct.[6] The
capacity-sp eed pro duct of an acousto optic mo dulator can b e expressed as
N 1=t=4f 1
where N is the numb er of resolvable spatial p oints, 4f is a frequency bandwidth of the acoustic wave, and t is a
switching time, resp ectively. A frequency band width 4f as high as 300 MHz can b e obtained by the acousto optic
material such as alpha-io dic acid.[7]Ifwewant to switch the de ector every 1sec, then with a safety factor of 2, the
numb er of resolvable p oints will b e at most 150. In order to overcome the disadvantage of the acousto optic b eam
de ector, a multistage digital b eam de ector was designed.[8] They demonstrated a 20-stage de ector consisting of a
series of nitrob enzene Kerr cells and birefringent calcite prisms. The laser b eam was de ected into a two dimensional
1024 1024 plane in every 2sec. This approach provided a great exibility and accuracy in controlling the de ection
angle. However, it required very high bias and switching voltage of several kilovolts, and also resulted in a large
power consumption.
1.4 Description of Optical Expanders
An optical expander takes as an input a b o olean pattern of size d = c log N bits, and expands it to a b o olean
pattern of size N bits, where c is a constant satisfying 1 c 2. Each expanded b o olean pattern must b e mutually
orthogonal to the others. Thus, the optical expander can b e viewed as either of the following: an electro optical line
deco der which converts d bits of optically enco ded binary information to up to N unique optical outputs, or a digital
b eam de ector which uses a control signal enco ded in d bits to de ect an input laser b eam into one of N directions. 2
More precisely, an optical expander takes as input one of N distinct b o olean vectors p ;p ;:::;p of length d.
0 1 N 1
We call these vectors the input patterns. Each input pattern is optically enco ded by using d pixels, each pixel b eing
either ON denoted by 1 or OFF denoted by 0. We will require that each input pattern has exactly d=2 pixels ON.
The optical expander pro duces a spatial output pattern r from given input pattern p . Each output pattern r is
i i i
one of N mutually orthogonal b o olean vectors of length N . In addition to our standard optical expander, we de ne
a generalized optical expander. The generalized optical expander is similar to the standard optical expander except
for one detail. Unlike the standard optical expander, each expanded b o olean pattern mayhave more than one ON in
its elements. In other words, the generalized optical expander creates a b o olean pattern of size N which is a bitwise
OR pro duct of some subset of the N mutually orthogonal b o olean patterns. The advantage of the generalized optical
expander b ecomes clear in certain applications. It can be used for broadcasting messages in a message switching
network. It can also be applied to a holographic memory system with a multiple readout capability, where the
bitwise OR, AND, or XOR pro ducts of several images data can b e directly obtained as a sup erimp osed output on
the detector array.
Our optical expander accepts an input pattern enco ded in d bits, and expands it into a pattern enco ded in N
bits. We wish to have an exp onential expansion, so d has to be represented by d = c log N for some constant c.
First, we describ e an optical expander with the constant c =2. Later, we will lo ok at an enco ding scheme with the
constant c 1 for a large d. Table 1 This allows us to pro duce a greater numb er of orthogonal patterns with the
same numb er of input bits. However, setting c = 2 o ers several advantages. First of all, it makes the co ding scheme
simple, since d = 2 log N o ers a co ding scheme where each p can be a concatenation of two binary strings: one
i
representing i in binary format, and the other representing i in one's complement binary format. Thus, p can b e
i
easily pro duced from the binary-co ded output from the electrical interface without any additional electrical mapping
interfaces. Second, it also makes the optical interconnection patterns from d optical inputs to the threshold array
regular, thus resulting in a simple implementation. Finally, it can provide an addressing scheme for the generalized
optical expander. We will discuss this enco ding issue later in detail.
1.5 Optical Expanders require Non-linear optical systems
A linear optical system can not b e used as an optical expander, since any linear mapping from an input of size d
creates no more than d linear indep endent output patterns. Thus, it is imp ossible to create a set of N>dmutually
orthogonal patterns byany linear optical system on d linear indep endent patterns.
1.5.1 Non Linear Optical Filters
Non-linearity can b e intro duced to an optical system bytwo metho ds. One can use a non-linear device. Thresholding
the input intensity at a certain level to pro duce output is a non-linear op eration. It can b e implemented by optical
non-linear devices such as optical logic etalon OLE[9] [10], or by electro optical non-linear devices such as self
electro optic e ect device SEED[11]. The other metho d is to translate an input into a spatial pattern, and then to
apply a linear lter at the fourier plane. An example is Theta mo dulation, where data are enco ded as a grating of
di erent orientations. In our optical expanders, we use non-linear devices.
1.5.2 Outline of Pap er
In section 2, we describ e applications of our optical expanders. In particular, section 2.1 describ es optical
interconnects in general. Sections 2.2 and 2.3 discuss holographic memory and message routing, resp ectively. We
start our discussion on the optical expanders in section 3. In section 3.1, we describ e our rst optical expander.
It consists of two parts: a linear part and a non-linear part. The linear part is a matrix-vector multiplier, and the
non-linear part is an arrayof thresholding devices. In section 3.2, we describ e our second optical expander. It
consists of log N identical switching cells. Finally, section 4 concludes the pap er.
2 Applications of Optical Expanders
Before we discuss sp eci c applications, we lo ok at the general issues. There are various applications which can b e
eciently implemented by use of optical expanders. Systems which require N distinct entry b eams either to an
N -sup erimp osed hologram or to an array of N devices may use an optical expander to generate N b eams from 3
optical input of c log N bits. Without our optical expander, such systems require either a b eam de ector to de ect
a laser b eam into one of N unique directions [12], or an electrically implemented line deco der which accepts log N
bits of binary information and creates one of the N mutually orthogonal b eams.
The conventional acouso- or electro-optical b eam de ectors have several drawbacks describ ed earlier. Electrically
implemented large line deco ders are not practical in terms of sp eed and wiring areas for a large N . As we mentioned
earlier, the I/O constraints limit the size of system which can b e practically implemented.
Our optical expander will provide an advantage to these devices by utilizing free space with exibility and accuracy
provided by digital op erations. The following sections describ e typical applications.
2.1 Optical Interconnects
Because of the availability of non-linear electrical devices as gates which are extensively used in the interconnection
network, electrically implemented interconnections are widely seen among many computer organizations.[13] [14]
However, the future of electric interconnections is not necessarily bright. The problem comes from its restricted
two dimensionality and RC delayoninterconnections.[15] This implies that with the current technology of electrical
interconnects, the interconnection area will so on o ccupy a large p ortion of a chip and the interconnection delay will
b ecome a b ottleneck in pro cessing.
These drawbacks do not exist in optical interconnections. Light b eams need not be con ned in a wave guide
such as an optical b er, but can travel freely through space. Light b eams can provide a great bandwidth, and the
propagation of light traveling through space or in a b er is not a ected by resistance, capacitance, or inductance.
Thus, optical interconnections may o er a high data transfer rate in a simple architecture by a set of light b eams
freely traveling through space. A large numb er of pap ers discuss the p otential of optical interconnections. Surveys
are found in [16] [2] [17].
Several theoretical studies have b een made to investigate the advantage of free space optical interconnects.[18] [19]
[20] The studies indicate that optical interconnects haveanadvantage over their electrical counterparts in terms of
area volume, sp eed, and p ower consumption for high sp eed communications except for the shortest distance within
chips. Furthermore, they suggest that optical interconnects may b e advantageous in large area VLSI circuits which
require high data rates and/or large fan-out. These results as well as work on various implementations of optical
interconnects, lead to our interest in designing our optical expander which can b e eciently used for a holographic
memory storage and a message routing network.
Finally, from the purely theoretical computational p oint of view, for a given problem, there is a lower b ound
2
on the circuit area and its computational time. One suchlower b ound on VLSI mo del called \AT b ounds" states
2 2 1
that AT = I , where A is the circuit area, T is the time used by the circuit, and I is information content of
the problem.See [21] In a three dimensional electro optical mo del describ ed by Barakat and Reif called VLSIO, the
3=2 3=2
similar lower b ound can b e expressed as VT = I .[22] This implies that as the information content b ecomes
larger, the VLSI circuit requires a larger and larger area to solve the problem in a xed amount of time. Using
optical interconnection as in VLSIO mo del overcomes this interconnection problem by utilizing space in a volume.
2.2 Holographic Memory Storage
Holograms can b e used to implement random access memory storage systems. [23] [24] [25] [26] [27] The basic idea of
holographic memory storage is that the data are arranged in blo cks which are stored in holograms. A large blo ckof
memory can b e retrieved at a any given time by using its corresp onding reconstruction b eam. This typ e of memory
is particularly suited for read-only applications, since the holograms can b e xed. However, dynamically mo di able
holograms such as photorefractive materials may be useful for active holographic memory storage systems. The
work in the 70s promised the advantage of holographic memory over other typ es of memory in terms of bit/volume
ratio, size, and throughput. However, the lack of appropriate recording materials and fast addressing metho ds kept
holographic memory b ehind the progress of conventional bip olar or MOS based memory.
Recently, the advance in recording materials suchasvarious photo crystals and the success in fabricating an array
of large numb er of micro lasers have provided a chance for holographic memory to b e eciently implemented. Several
prototyp es of such a memory storage system have b een develop ed at Micro electronics and Computer Technology [28]
and Bellcore [29].
1
Information content is the numb er of bits that must cross a b oundary in order to solve the problem. The b oundary separates the
circuit into two sides, each of which holds approximately half the input bits. 4
In a typical holographic memory storage system, the data are organized in blo cks.Fig. 1 Each blo ckisatwo
2
dimensional bit image consisting of L L pixels total of L pixels p er blo ck The detection of a bit can b e made
by using a PIN photo detector. N blo cks of images are stored in either a single multiple-exp osure volume hologram
or an arrayof N N holograms. Ifamultiple-exp osure volume hologram is used as a storage medium, we need N
mutually orthogonal patterns to retrieve each blo ck. In other words, in order to read out each blo ck, we need a total
of N b eams such that each b eam has a distinct incident angle to the holographic medium. When one of these b eams
illuminates the holographic medium, the corresp onding stored image is formed on the detector array. The detector
array can b e an arrayof LL PIN photo detectors which convert optical signals to electrical signals. If instead an
array of thin holograms is used, each hologram can b e separately illuminated by a distinct b eam. Again, N mutually
orthogonal patterns are necessary to retrieve N blo cks of image.
Using a laser b eam with wavelength , the numb er of bits which can b e stored in a unit volume is prop ortional to
3 2
1= , whereas in a unit area it is prop ortional to 1= . Therefore, the rst approach can achieve a large bit/volume
ratio. However, the di raction eciency of each image signi cantly decreases as the number of stored images
increases. The second approach is limited by the maximum di raction angle at which the image is formed on the
detector array.[25] As the array b ecomes large, it b ecomes dicult for the hologram at the corner to form its image
on the detector array. This may b e improved byhaving several detector windows on the back of the hologram array,
and merging these windows together into the detector array with some integrated optics. In any case, wemust create
N mutually orthogonal addressing patterns to address N images. These patterns can b e electro optically created by
our optical expander with input size d = c log N , where 1 c 2. Thus, the electrical interface to the holographic
memory system needs only a pattern enco ded in d pixels: this pattern is electro optically expanded into a pattern of
N pixels to retrieve one of N blo cks of images. Without our optical expanders, an alternative approach requires a
set of N orthogonal addressing patterns whichmust b e either electrically created by using a VLSI line deco der, or
acousto- or electro optically created by using a b eam de ector. As we mentioned earlier, these previous approaches
have disadvantages. Thus, our optical expanders o er an advantageous approach.
2.3 Message Routing
Interconnection networks in parallel pro cessing computers are very imp ortant sub jects. There are manyintercon-
nection networks for di erent applications, since di erent algorithms require di erent degree of globality of the
interconnects.
Message routing is a task where messages are to b e moved among the various pro cessing units. We assume there
are N pro cessors and N messages, where each pro cessor has a distinct message with a distinct destination address.
Then, simultaneously, each message is routed from its originating pro cessor to its destination pro cessor.
In parallel computing there are several classes of problems. Some problems require no communications or only
xed lo cal communications among pro cessing units; examples are image pro cessing and various matrix op erations
using systolic algorithms. Some problems require very intensive and sometimes dynamic communications among
distant pro cessing units; a fourier transform requires communications among distant pro cessing units, but can be
implemented by a xed connection; problems often found in forecasting and AI, on the other hand, require global
and dynamic interconnections.
In a general purp ose parallel computer, the full advantage of parallelism will only b e realized if each pro cessing
unit has a direct communication path to every other pro cessing unit. In such a condition, each pro cessing unit can
pro cess its data without having serious communication delays, thus resulting in high overall throughput rates. If this
condition is not satis ed, the communication cycle may far exceed the pro cessing cycle, and this will cause a serious
b ottleneck in the overall system sp eed. For example, the Connection Machine is a 65,536 pro cessor single instruction
multiple data stream SIMD computer develop ed at Thinking Machines [30]. Its data transfer time is at least 1000
times slower than the pro cessor instruction step, intro ducing a serious b ottleneck. Thus, we recognize that message
routing is a very crucial problem in designing ecient parallel computers.
Among various message routing networks the highest level of interconnection is a crossbar network which uses
2
N interconnects to connect N source units and N destination units. The number of electrical interconnection
wires required by each pro cessing unit to communicate with the other pro cessing unit on- and o -b oard will limit
the feasible size of the network. The prop ertyof light b eams which we brie y mentioned ab ovemay give a great
p otential for an alternative high-sp eed optical crossbar typ e of networks.
There are several optical interconnection networks which have already b een prop osed. One is optical crossbar
network. [31] [32] [33] [34] The optical crossbar network typically uses an N N spatial light mo dulator SLM to 5
connect N source pro cessors to N destination pro cessors. Each source pro cessor uses a column of the N N SLM to
address one of N distinct destination pro cessors. The advantage of this optical crossbar is that once all the entries
of the N N SLM are set, the message can be transmitted at avery high data rate, namely at an optical pulse
mo dulation rate. Fig. 2 This matrix-vector multiplier based crossbar network has two drawbacks. One is that
at most 1=N of the power incident on the SLM will reach the detector. The other is that it takes a long time to
electrically set an N N SLM.
Another network uses multiple stage optical switching networks. [31] [35] [36] [37] This is an optical implementa-
tion of various xed multiple stage electrical interconnection networks. The multiple stage networks typically require
the setting of only O N log N switches to connect N sources to N destinations. However, they require O log N
steps to route a message and some top ologies do not guarantee a non-blo cking routing.
Toovercome the drawbacks of previously prop osed systems, a network which used xed multiple-exp osure holo-
grams to connect N source units to N destination units was presented.[38] [39] Fig. 3 Unlike the matrix-vector
multiplier based optical crossbar networks, this holographic interconnection network uses xed holograms to steer
spatially enco ded light b eams transmitted from the source units toward the destination units.
The basic idea of steering each light b eam to its destination is to use holographic asso ciative matching. Each
pro cessor has a xed hologram which implements connections to other pro cessors. Each pro cessor can establish its
connection to other pro cessors by illuminating its hologram by a reference b eam. They demonstrated a 4 pro cessor
system with two di erent spatial enco ding schemes. The rst scheme used a set of N mutually orthogonal patterns
to enco de destination addresses. In this scheme, the system crosstalk caused by a false matching was minimized.
However, as with the traditional crossbar networks, an N N SLM was required. To gain an advantage over the
traditional matrix-vector multiplier based crossbar network, the second scheme employed a shorter address to reduce
the size of SLM, which further led to a shorter recon guration sp eed. In this scheme, each address was enco ded
p
as a non-orthogonal pattern of size 2N bits. The system crosstalk was distributed among di erent detectors by
increasing the dimensionality of the detector plane. As we see, if we can electro optically create N mutually orthogonal
patterns, an ecient crossbar network can b e implemented by using them as addressing patterns to the holograms.
Our optical expanders can create such N mutually orthogonal optical patterns from input optical patterns of size d.
3 Optical Expanders
The optical expander is a non-linear electro optical lter which creates a large number N of mutually orthogonal
patterns. It uses electoro optical devices with at most d b o olean input bits. Here d is no greater than 2 log N . There
maybevarious ways to implement an optical expander. Our rst approach uses matrix-vector multiplication followed
by threshold op eration. We call this the Matrix-Vector Multiplier MVM Optical Expander. Our second approach
is based on an novel idea to use only log N identical switches to digitally de ect the laser b eam. We call this the
Digital Beam De ector DBD Optical Expander. The following sections rst describ e the MVM Optical Expander,
and then the DBD Optical Expander.
3.1 MVM Optical Expander
First, we describ e two di erent enco ding schemes which can b e used in this mo del.
3.1.1 Enco ding of Input
Our optical expander uses an optical matrix-vector multiplier. The vector represents one of N distinct input patterns
p ;p ;:::;p of length d. Each pattern has exactly d=2 ONs denoted by 1 and OFFs denoted by 0, whichkeep
0 1 N 1
the total p ower of each pattern equal. A set of such patterns can b e easily obtained.
We presenttwo enco ding schemes, eachof which has an advantage and a disadvantage. The rst scheme uses
d = 2 log N bits. Each pattern p consists of two bit strings: one enco ding i in binary and the other enco ding i in
i
one's complement. This enco ding, which is called the dual rail co ding [40], is often used in optical computing [41].
If this is considered to o large, we can set d smaller than this value. To do this, we can recursively enumerate a set of
bit strings which has d=2 1s and 0s. We show that this gives d log N for large d. This can signi cantly reduce the
numb er of bits required to pro duce a large N . The analysis of this enco ding can b e done by using Stirling's formula
p
x x
x! 2xx e . To nd the asymptotic value of c as d b ecomes large, wehave; 6
r
p
d d
d! 2dd e 2
d
d=c d
2 = = = 2 2
2 d d
d=2
d=2! 2d=2d=2 e d
Thus, c can b e written as
d
c 3
1 2
d + log
2 d
>From this equation, we see that as d approaches in nity, c approaches 1. This enco ding pro duces a signi cantly
larger N as d increases. For example, d = 14 yields N = 3432 in this enco ding, whereas in the dual rail enco ding, we
7
have2 = 128. The advantage of the dual rail enco ding comes into play when it is used for the generalized optical
expander. We will explain this in the next section.
3.1.2 Matrix-Vector Multiplication
We describ e the general idea of how to create an orthogonal b o olean pattern of size N from a given input of size d.
We use matrix-vector multiplication as follows. The matrix which is of size N d consists of N rows, where the i-th
rowisp . We can consider this multiplication as a b o olean matrix-vector multiplication by using optical devices. Let
i
p b e an input to the matrix-vector multiplier. Then, the result of the multiplication is an output vector of length
k
N , where the i-th element of the vector is the inner pro duct of p and p .
k i
0 1 0 1
p :
0
B C B C
p :
1
B C B C
B C B C
::: :
B C B C
B C B C
0
p 0
k
T
B C B C
p = 4
k
B C B C
::: :
B C B C
B C B C
p d=2
k
B C B C
@ A @ A
::: :
p :
N 1
T
Here, p represents the transp osed vector of pattern p .
k k
0
0
The output vector has value d=2 at the k -th p osition and has a 0 at the k -th p osition, where p is the bit
k
wise complementofp . We apply a threshold op eration at a certain intensity level to pro duce one of N mutually
k
orthogonal patterns of length N .
In the generalized optical expander, input vector p may have more than d=2 1s. Thus, the output after the
k
threshold op erations mayhave more than one 1. A control signal of size d N 1. This is obvious, since the numberofsuch distinct patterns is 2 ,thus requiring N bits to cho ose one of them. If the enco ding with c = 2 is used, it is easy to set several p ositions of the output to 1s. We can consider each elementof the output to b e a leaf of a complete binary tree of hight d=2. In the dual rail enco ding, each bit in the input has its corresp onding complement. Setting a particular bit p osition to 1 and its complement to 0 corresp onds to traversing a no de to its right child. Similarly, setting a paricular bit p osition to 0 and its complement to 1 corresp onds to traversing a no de to its left child. Thus, only one leaf can b e nally reached. In order to reach several leaves, b oth children must b e traversed. Setting b oth bits to 1s corresp onds to traversing b oth children. In this way, the output can have several 1s. The following section describ es optical matrix-vector multipliers and threshold op erations in more detail. 3.1.3 Optical Matrix-Vector Multiplier In this section rst we review the traditional matrix-vector multiplier. After that, we describ e the matrix-vector multiplier whichwe use in our optical expander. Let X be a vector of length N , and let A b e a matrix of size N N . The optical implementation of a parallel, real-time matrix-vector multiplier is shown in Fig. 4. [17] The input vector is represented by an arrayof N light emitting dio des LEDs, where the lightintensity of each LED represents the value of each element in the vector. The matrix is represented by a transparency matrix mask of size N N . The transmittance of the i; j entry of the mask is prop ortional to the value of a . The light i;j from each LED is spread onto the corresp onding column in the transparency mask. Thus, the intensity of the light 7 passing the i; j entry is prop ortional to x a . The light from a whole row is collected, and fo cused onto the j i;j corresp onding p osition of the output vector Y . The intensity of the light at the i-th p osition of the output vector is then prop ortional to the i-th element of the pro duct. For the optical expanders we need to multiply a matrix A of size N d with a vector X of length d. We can arrange the input vector, the matrix mask, and the output vector p p into nearly square shap es. This means the vector of length d is formatted in d d, and the matrix of size d N p p p p is formatted in dJ N dJ N . The output vector of length N is formatted in N N . The examples are shown in gure 5. To clarify the idea of this multiplier, we describ e an example with size d =4, N =6. The input vector of length p p d is represented bya d d either LED or laser dio de LD arrayinrow ma jor order. If d is a square of some integer, the array can b e square. Otherwise, the arraymay not b e square. The matrix mask is divided into d blo cks p p of equal size. These blo cks are arranged in d d matrix form. The i-th column of the matrix is represented by the i-th blo ckinrow ma jor order. Each blo ck has N elements where the transparency of each element represents an elementofthe column. Again, if N is a square of some integer, the shap e of each blo ck b ecomes square. The elements of each column are arranged in row ma jor order in its corresp onding blo ck. In op eration, plane waves are pro duced from the dio de array by lenses. The plane wave from the i-th dio de illuminates the i-th blo ck which enco des the i-th column. Thus, the plane wave passing the i-th blo ck represents the vector pro duct of the i-th column of the matrix and the i-th element of the vector. All these plane waves are p p sup erimp osed at the detector array, and thus the pro duct is obtained. The output vector is obtained in N N matrix form, where the elements are organized in row ma jor order. We use this idea to construct a multiplier for our optical expander. We represent the input vector X byanLD array of size d formatted in a square shap e. The output vector Y of size N can also b e formatted in a square shap e. An arrayof lenses is placed in front of the LD arrayso that a set of plane waves are obtained by the lenses. Each plain wave then illuminates its corresp onding matrix mask. Again, an array of lenses with a fo cal length f is placed distance f away from the matrix masks. This makes a set 1 1 of fourier transform of the matrix mask patterns to b e pro duced at the back fo cal plane of the lens array. A lens with a fo cal length f is placed distance f away from the fourier plane, so that the matrix mask patterns will b e 2 2 sup erimp osed on the back fo cal plane of this lens. We note that there will b e a phase shift for each mask pattern at the sup erimp osed image. The phase shift is prop ortional to the lateral shift of the mask pattern from the optical axis of the second lens. It can b e ignored when the intensity is pro cessed at the image. This approach seems to require complex matrix mask construction. However, if the detectors are formatted in a two dimensional plane with enco ding scheme using c = 2, it will b ecome very regular. This is shown in gure 6. The interconnection patterns from an LD to the threshold array can b e group ed into d=4 di erent patterns. Then, each interconnection pattern is one of the d=4 patterns with one of four di erent orientations. As we see from the gure, these basic d=4 patterns are regular. Figure 6 compares an example of the two enco ding schemes. 3.1.4 Threshold Op eration The purp ose of p erforming a threshold op eration is to intro duce a non-linearity. A non-linearity is necessary to generate a set of N mutually orthogonal b o olean patterns from a set of N distinct d-linear indep endent b o olean patterns. Here, wehaved2 log N . As we recall, the output vector from the matrix-vector multiplier has a value 0 0 d=2atthe k-th p osition and has value 0 at the k -th p osition, where p is the bitwise complement of input pattern k p . i One p ossibilitytomake the output vector orthogonal is to set the threshold value at intensity I = I d=2, where 1 I is a factor which corresp onds the intensity of a one. The thresholded output b ecomes High, if the intensity is 1 I I d=2, and b ecomes Low otherwise. This will pro duce one of the mutually orthogonal patterns, since there is 1 exactly one p osition in the output vector which has intensity I = I d=2. This metho d will work as long as d is 1 small. However, when d b ecomes relatively large, there is a practical problem. The problem is caused by the physical limit of how nely we can threshold the intensity. When a threshold op eration at intensity I = I d=2 is p erformed, 1 intensity I = I d=2, and intensity I I d=2 1 must b e distinguished. This will b ecome dicult when d b ecomes 1 1 large. Our solution to this problem is the following. Since the output vector also has exactly one p osition which has intensity I =0,we threshold intensityatI =0. This means wemust distinguish intensity I = 0 from intensity I I . Then, the complement of the thresholded output b ecomes orthogonal. In this approachwe need not concern 1 the physical limit of the thresholding device mentioned ab ove, but rather the sensitivity and dynamic range of the device. 8 However, in this approach, twotyp es of intensity noises must still b e considered. One is caused by the voltage uctuation of the driving circuit and the other is a background intensity caused by di raction from di erentchannels. These noises must b e kept low enough so that I = 0 and I I can b e distinguished by the detectors. This can b e 1 improved by using di erential devices at the threshold devices. The generic mo del of a di erential device consists of two PIN photo dio des connected electrically in series.Fig. 7 The logical state of the device is de ned by the ratio of the intensity of two light b eams falling on the two photo dio des. In our threshold device, one photo dio de is used to detect the light b eam from the LD array, and the other is used to cancel the noise by biasing the voltage drop caused by the rst photo dio de. Several optical or electro optical switches have b een designed and demonstrated. The optical logic etalon OLE [9] [10] and the self electro optic e ect device SEED [11] are the most common devices. As for a di erential device, the Symmetric SEED S-SEED has b een demonstrated [42] and widely used as a basic comp onent.[43] [44] [45] Both the OLE and the SEED have a contrast ratio of approximately 5:1. Recently, low-threshold electrically pump ed vertical-cavity surface-emitting microlaser dio de arrays SELDA's have b een develop ed [29]. Each laser can b e as small as a few micrometers. The output light from each laser has high contrast, since an OFF state pro duces no light. The optical output p orts such as laser dio des or mo dulators are based on GaAs technology. The input p orts such as photo dio de detectors are based on Si technology. This makes it dicult to fabricate laser dio des and photo detectors on the same single substrate. Several metho ds of fabricating hybrid systems have b een investigated [46] [47] [48], and may lead to a large array of optical threshold devices with a very high contrast. 3.1.5 Overall Architecture and Analysis of MVM Optical Expander The overall architecture is shown in gure 8. The output b eam can b e used to either address a hologram, or send a message to a destination. We now analyze three issues concerning the MVM Optical Expander. The rst issue is the di raction limit imp osed by a nite size of the lenses. The second concerns the fan-out factor which dep ends on the minimum de- tectable p ower of the detectors at a given op erating condition. The third concerns the dynamic range of the detectors. Di raction Limit Here, we consider the di raction limit of the system. Figure 9 shows a side view of the system. We examine the image formation of a arbitrary matrix mask on the detector plane. In this gure, each mask has a dimension of M M . The rst lens L has an ap erture function a x; y and has 1 1 a fo cal length of f . The second lens L has an ap erture function a x; y and has a fo cal length f . Note there is 1 2 2 2 a lateral shift of the mask from the optical axis of the second lens L . Let s and s represent this shift along the 2 x y x-axis and the y -axis, resp ectively. For our analysis, we use rectangular ap ertures for b oth lenses. Thus, wehave x y a x; y = 5 1 D D 1 1 x y a x; y = 6 2 D D 2 2 Here, denotes a rectangular function, D and D are the size of each ap erture. 1 2 We assume that a plain wave is incident on the mask. We denote the mask pattern by g x; y . Then, the light distribution E x; y just in front of the lens L is given as [49] z 1 1 E x; y =[gx; y h x; y ]a x; y 7 z f 1 1 1 Here, represents a convolution op erator and h x; y isapoint-spread function of a distance f . This p oint- f 1 1 spread function is the light distribution of a p oint source at distance z = f away. Thus, wehave 1 2 2 1 x + y h x; y = exp jk z + 8 z jz 2z Using the Fresnel approximation, we can write the light distribution at the back fo cal plane of the lens L as 1 2 2 1 x + y x; y = x; y g 9 E exp jk f + FfE z 1 z 2 1 jf 2f 1 1 9 Here, wehave the spatial frequencies = x=f and = y =f . x 1 y 1 We use the following relations to expand the fourier transform term. " ! 2 2 2 z + x y Ffh x; y g = H ; = exp jk z 10 z z x y 2z Thus, wehave 2 2 1 x + y x y x y x y E G H A x; y = exp jk f + ; ; ; z f 1 1 2 1 jf 2f f f f f f f 1 1 1 1 1 1 1 1 1 y x y x A 11 = exp[j 2kf ]G ; ; 1 1 jf f f f f 1 1 1 1 1 Again, the light distribution just in front of the second lens L can b e written by using a p oint-spread function 2 as E x; y =[E x s ;y s h x; y ]a x; y 12 z z x y f 2 3 2 2 Here, note that we shift the optical axis by s and s . x y At the detector plane, 2 2 x + y 1 exp jk f + FfE x; y g 13 E x; y = 2 z z 3 4 jf 2f 2 2 Here, the spatial frequencies = x=f and = y =f . x 2 y 2 Expanding the fourier transformation term yields the following. 2 2 1 x + y x y x y E x; y = exp jk f + FfE x s ;y s gH ; A ; z 2 z x y f 2 4 2 2 jf 2f f f f f 2 2 2 2 2 2 f s x s y 1 x y = exp[j 2k f + f ] exp j 2 + 1 2 f f f 2 2 2 f f f f x y 1 1 1 1 g x; y a x; y A ; 14 1 2 f f f f f f 2 2 2 2 2 2 Since wehave a rectangular ap erture function. wehave 2 A ; =D sincD sincD 15 2 x y 2 x 2 y 2 s x f s y x 1 y j 2 E x; y = exp[j 2kf + f ] exp + z 1 2 4 f f f 2 2 2 D x f f f x f y D y 2 1 1 1 1 2 2 sinc g x; y D sinc 16 2 f f f D f D f f 2 2 2 1 2 1 2 2 Thus, the resolution of the mask image on the detector is limited by the convolution with the sinc functions. The ap erture at the rst lens has to b e at least as large as the mask size. Thus, wemust have: D M 17 1 The ap erture at the second lens determines the resolution of the image. When D approaches in nity, the 2 convolution term b ecomes a delta function, thus forming a p erfect image on the detector plane. A nite sized D 2 reduces the resolution of the image. Wenow examine the minimum detector spacing de ned as the minimum spacing between the centers of two detectors, the detector size and the total numb er of detectors N . First, the detector size and the minimum detector spacing are limited by the width of the principal lob e of the sinc function which is expressed as 10 2f 2 4x = 18 D 2 Assuming that the f -numb er of lens L is f= = 1 and the wave length is 0.8m, we get 4x as 1:6m. The detector 2 size and the spacing between the detectors are determined by numerically examining the convolution of the sinc function with an input image. We calculated the convolution term with several input images and the results are shown here.Fig. 10 We notice that a very small detector size with avery small spacing may cause a problem for a large N , since the in uence of the nonprincipal lob es of the sinc function is not negligible. As we increase the detector size, it will b ecome negligible. Using this analysis, we can determine the overall size. In order to minimize the signal skew caused by the di er- ence in the propagation delay,we minimize the lateral shift of the masks from the optical axis of the second lens. To minimize the lateral shift, wemust use the smallest p ossible mask, the size of which is determined by the resolution of the mask. The resolution of the mask patterns using electron-b eam lithography is ab out 0:5m.[50]We assume that we use an arrayof55m PIN photo dio de typ e of detectors such as SEEDs. As we lo ok at the plots, we notice the convolution of x=3:6m with the sinc function pro duces a sp ot of diameter ab out 5m. Further examining the graph, we nd that if we take the spacing larger than 5m, the e ect of the second lob e is negligible for N up to 100 2 . If the mask pattern has a sp ot expressed as x=0:5m, we need the ratio of the fo cal length as 3:6=0:5 8. This b ecomes also the magni cation factor of the mask image to the detector image. If the sp ots on the mask pattern are separated by 1:5m, this separation corresp onds to 12m on the detector plane. Each mask can be made in p p p p p 1:5 N 1:5 Nm. The rst lenses L shave a fo cal length 1:5 Nm with an ap erture 1:5 N 1:5 Nm. The 1 p p p second lens has a fo cal length 12 Nm with an ap erture 12 N 12 Nm. The detector array has N 5 5 m p p detectors fabricated in 12 N 12 Nm. Then, the light path from the input lasers to the detectors is no more p p than 30 Nm. The free space propagation delay of this distance is 0:1 N psec. We are using an o axis imaging. The di erence in the propagation delay caused by the o axis fourier transform is prop ortional to the lateral shift 0 0 0 of the mask from the optical axis. Thus, the maximum di erence is s + s =v . Here, s is the largest shift of x y c x p 0 0 0 the x-axis, s is of the y -axis, and v is the sp eed of light. In this case, wehaves =s =0:75 NJ log Nm, y c x y p thus the maximum di erence is 5 NJ log N fsec. For example, when N = 1024, this b ecomes 0.5 psec, and the total propagation delayis 3.2 psec. As we further increase N to 65,536, the di erence is 5.12 psec, and the total propagation delay is 25.6 psec. Thus, a bit rate around a giga hertz do es not cause a signal skew problem. Detectable Power Limit The largest N is limited by the minimum detectable power of the photo detectors. The optical signal emitted from one source is distributed to N p ositions of the mask. Let P denote the total radiation p ower from the source. 0 Then, ignoring other losses, each sp ot on the mask receives at most P =N optical power. In fact, if each sp ot is 0 0:5 0:5m with a spacing 1.5 m, it receives only 0:25P =N power. Fig. 11 This p ower must corresp ond to the 0 intensity I at the detector plane. 1 In this section, we analyze the relation b etween P , N , and a bit rate B . The bit error rate BER of a given 0 optical link is given by the integral [51] Z 1 2 1 x p P E = exp dx 19 2 2 Q where Q is a parameter relating to the desired error rate P E . For the case of a PIN photo dio de with a high- imp edance FET front-end ampli er, the detector sensitivity can b e approximated to p hv c 2 P = Q hi i 20 n e 2 where h is Planck's constant, v is the sp eed of light, is the wave length, e is the fundamental charge, and hi i is c n the noise currentpower at a given bit rate. At high bit rates, the noise currentpower is dominated by the thermal channel noise term. This can b e written as