Advances in Optical Technologies

Volume 2014 (2014), Article ID 569684, 8 pages

http://dx.doi.org/10.1155/2014/569684

Research Article

Demonstration of Controlling the Spatiotemporal Dynamics of Optical Near-Field Excitation Transfer in Y-Junction Structure Consisting of Randomly Distributed Quantum Dots

Wataru Nomura,¹ Makoto Naruse,² Masashi Aono,3,4 Song-Ju Kim,⁵ Tadashi Kawazoe,¹ Takashi Yatsui,¹ and Motoichi Ohtsu¹

¹School of Engineering, The Nanophotonics Research Center, The University of Tokyo, 2-11-16 Yayoi, Bunkyo-ku, Tokyo 113-8656, Japan

²Photonic Network Research Institute, National Institute of Information and Communications Technology, 4-2-1 Nukui-kita, Koganei, Tokyo 184-8795, Japan

³Earth-Life Science Institute, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan

⁴PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi-shi, Saitama 332-0012, Japan

⁵WPI Center for Materials Nanoarchitectonics (MANA), National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan

Received 10 March 2014; Accepted 24 March 2014; Published 27 April 2014

Academic Editor: Chennupati Jagadish

Copyright © 2014 Wataru Nomura et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Solution searching devices that operate on the basis of controlling the spatiotemporal dynamics of excitation transfer via dressed photon interactions between quantum dots have been proposed. Long-range excitation transfer based on dressed photon interactions between randomly distributed quantum dots is considered to be effective in realizing such devices. Here, we successfully controlled the spatiotemporal dynamics of excitation transfer using a Y-junction structure consisting of randomly dispersed CdSe/ZnS core-shell quantum dots. This Y-junction structure has two “output ends” and one “tap end.” By exciting one output end with control light, we observed increased excitation transfer to the other output end via a state-filling effect. Conversely, we observed reduced excitation transfer to the output ends by irradiating the tap end with control light, due to excitation of defect levels in the tap end. These results show the possibility of controlling the optical excitation transfer dynamics between multiple quantum dots.

1. Introduction

Light excitation in quantum dots (QDs) generates dressed photons, which are light fields localized in the vicinity of the QDs, giving rise to dressed photon interactions with other nearby matter, as well as excitation energy transfer via these interactions [1]. In particular, various optical functional devices, such as logic gates called nanophotonic devices [2–4], light-harvesting devices [5], and optical signal transmitting systems [6, 7], have been realized using QDs formed of CuCl, ZnO, InAs, CdSe, and so forth, based on optical near-field excitation transfer between QDs. Nanophotonic devices have been shown to function as logic gates, such as AND, NOT, and XOR logic gates [2–4]. These devices consist of two or three closely spaced QDs having different energy levels, and by inputting a light beam serving as a power supply and another light beam serving as a control signal, excitation energy transfer between the QDs is controlled so that the light emitted from one of the QDs serves as the output.

On the other hand, novel solution searching and decision making devices using a QD array provided with multiple output QDs have recently been proposed [8–10]. In these devices, by inputting control signals to the output QDs based on certain rules, the probability of the optical excitation being transferred uniformly to each QD is controlled to obtain a solution. In these operations, it is necessary to control the spatiotemporal dynamics of the optical excitation transfer between spatially distributed QDs. The features of these devices are that they operate at high speed with low energy consumption since they are driven by excitation transfer based on dressed photon interactions [11]. The proposed solution searching device finds a solution to a nondeterministic polynomial (NP) time complete problem, the satisfiability problem (SAT), much faster than the WalkSAT algorithm, which is one of the fastest stochastic local search algorithms [9, 12]. The proposed decision making device exhibits high efficiency and adaptability in solving decision making problems [10].

However, although it would be ideal to utilize the models discussed in [8–10], it would be technically challenging to directly implement their architecture. One reason for this is that, in order to implement and operate this device, it is essential to fabricate a QD array in which QDs of strictly selected dimensions and number must be precisely located with nanometer-order positioning precision, as reported in [8], which is technologically demanding from the viewpoint of fabrication as the first step. In addition, another reason is the technical challenge of having to guide the control light to multiple output QDs and separately observe the output signals coming from them.

In contrast to this approach, we distributed an extremely large number of identical-size QDs randomly, and by using a structure in which an output QD is provided at the terminal end, we experimentally demonstrated the possibility of optical excitation energy transfer over distances of micrometer order via repeated dressed photon interactions between the QDs [6]. Using such a randomly distributed QD structure overcomes the fabrication difficulties mentioned above, and because of the possibility of scaling up the technique to micrometer order, it is expected that it will also be possible to overcome the technical challenges related to control and observation.

In the work reported in this paper, using a structure in which multiple QDs are randomly distributed, we experimentally investigated control of the spatiotemporal dynamics of optical excitation transfer based on dressed photon interactions between the QDs, which is a requirement for the realization of solution searching and decision making devices. Section 2 describes fundamental principles. Section 3 describes the fabricated device, experimental setup, energy transfer measurement, and characterizations. Section 4 concludes the paper.

2. Principle

First, we will explain optical excitation transfer due to dressed photon interactions and suppression of this process. As shown in Figure 1, a small QD (QD1) and a large QD (QD2) are placed in close proximity. A first excited state E11 in QD1 is resonant with a second excited state E22 in QD2 due to their size ratio of 1 : 1.43 [13]. When QD1 is excited by excitation light , a dressed photon is generated in the vicinity of QD1. The excitation energy in QD1 is transferred to E22 by dressed photon interactions and then relaxes to E21, which is the first excited state in QD2, via intersubband relaxation. Thus, the excitation energy is transferred unidirectionally (Figure 1(a)). This is observed as emitted light from QD2. However, if E21, which is the level to which the excitation relaxes, has already been excited by control light , the excitation energy transferred to E22 cannot relax to E21, but it is transferred back to E11, and nutation occurs between E11 and E22 (Figure 1(b)). This is observed as emitted light from QD1, or when there is a QD in the vicinity, this leads to an increase in the probability of excitation transfer to that QD. Based on this principle, assuming that QD1, having multiple QD2’s in the vicinity, is excited, by irradiating one QD2 with , the probability of excitation transfer to the other QD2’s is changed; in other words, it is possible to control the spatiotemporal dynamics.

(a)

(b)

(c)

(d)

(e)

Figure 1: (a), (b) Schematic illustration of optical excitation transfer from QD1 to QD2 via dressed photon interactions, where (a) shows the case where energy transfer occurs from QD1 to QD2 and (b) shows the case where control light is incident on QD2, thus suppressing energy transfer. (c) Schematic diagram of long-range optical excitation transfer system. Green and yellow circles represent QD1’s and QD2’s, respectively. (d) Schematic diagram of experimental model. Green and yellow regions represent groups of QD1 and QD2 structures consisting of randomly distributed QD1’s and QD2’s, respectively. (e) Energy diagram schematically showing the experimental model from point B to point A.

Next, we will explain optical excitation transfer based on dressed photon interactions in a randomly distributed QD structure. In this structure, as shown in Figure 1(c), identical-size QDs (QD1) are randomly distributed in close proximity to each other, and QD2’s are placed at one location to serve as an output. The input light incident on the QD1’s at one end undergoes repeated dressed photon interactions between the QD1’s and eventually arrives at the QD2’s, where it is output as emitted light from E21. Because the optical excitation is transferred with high probability by QDs having strong dressed photon interactions, optical excitation transfer with high energy efficiency is possible regardless of the QD1 distribution, and it has been experimentally verified that it is possible to transfer the excitation over a distance of more than 10 μm [6].

The experimental model used in this paper is shown in Figure 1(d). The green color in the figure represents a Y-junction QD1 structure consisting of a group of QD1’s randomly distributed in close proximity to each other, in the form of strips with a width of 1 μm and a length of 7.5 μm extending in a radiating pattern in three directions from the center O. The ends of the Y-junction are labeled A, B, and C. QD2 structures serving as output ports, where QD2’s are grouped together in squares with a side length of 3 μm, are connected at A and B so as to overlap with the QD1 structure. From A and B, called “output ends”, we obtain output signals that are observed as radiation. In addition, A and B also serve as terminals for inputting control light. On the other hand, end C, which we call a “tap end,” is not connected to any QD2’s. The detailed experimental conditions used in this model are omitted. Figure 1(e) shows an energy band diagram in which the section including points A-O-B is illustrated one-dimensionally. When excitation light , serving as an input signal, is incident on the center point O in the QD1 structure, optical excitation is generated in QD1 and is transferred through the QD1 structure by undergoing repeated energy transfer due to dressed photon interactions between the QD1’s. Since the optical excitation transferred to QD2 reaches E21 via intersubband relaxation, no reverse flow of the transfer occurs, and therefore the optical excitation is transferred to the ends, A and B, where the QD2 structures are provided. At this time, if end B is irradiated with control light that excites QD2, QD2 is excited, and the excitation light transferred from QD1 does not relax to E21, and thus excitation transfer is not observed. Therefore, a reverse flow of the optical excitation occurs at end B, and the level of excitation transfer to A, which is the other output end, increases. At locations where the excitation transfer is increased, the emitted light and from QD1 and QD2 both increase, and, therefore, we can confirm that the excitation transfer is controlled by evaluating the increase in at end A when end B is irradiated with control light . Note that even if end C is irradiated with control light , a change like that occurring at points A and B is not expected to occur since the excitation state of the QDs does not change.

3. Experiments

3.1. Device

In our experiments, we used commercially available CdSe/ZnS core-shell QDs (manufactured by Quantum Design, Inc.). As QD1 and QD2 described in the previous section, we selected QDs with core-diameters of  nm and  nm, respectively. Their first excited states, estimated from absorption spectra, were  eV and  eV, respectively. As described above, E11 is resonant with E22, which is the second excited state of QD2 [14]. The outside of the shell was modified with a carboxyl group. The solution was dropped onto a substrate and allowed to dry, affording a group of randomly distributed QDs located in close proximity with a center-to-center distance of about 10 to 20 nm.

We fabricated samples using the procedure shown in Figure 2(a). First, by using e-beam lithography on a silica substrate, in e-beam resist (ZEP-520A, Zeon Corp.), we drew locations at which QD2’s were to be deposited (Figure 2(a)(i)). Next, the QD2 solution was dropped onto the substrate and allowed to dry (Figure 2(a)(ii)). Then, after removing the resist, we obtained the QD2 structure on the substrate (Figure 2(a)(iii)). Next, e-beam resist was applied on this structure (Figure 2(a)(iv)), and using e-beam lithography, we drew Y-junction locations where QD1’s were to be deposited by aligning the positions so as to contact the QD2 structure (Figure 2(a)(v)). After developing the resist, we dropped the QD1 solution onto the substrate (Figure 2(a)(vi)). Finally, the resist was removed to obtain the Y-junction QD1 structure (Figure 2(a)(vii)). Figure 2(b) shows a fluorescence microscope image obtained when the fabricated sample was excited with a mercury-vapor lamp. The QD1 structure and the QD2 structure were observed as green and yellow colors, respectively. Because points A and B were set at positions where the QD1 and QD2 structures were initially in contact, in practice the distances OA and OB of the fabricated samples were about 6.8 μm and 3.3 μm, respectively. Point C, which did not possess the QD2 structure, was set at a position 7 μm away from O.

(a)

(b)

Figure 2: (a) Cross-sectional illustrations of sample fabrication process. (i) E-beam resist is applied to SiO₂ substrate, and regions where QD2 structures are to be formed are drawn by e-beam lithography. (ii) QD2 solution is dropped onto substrate and allowed to dry. (iii) E-beam resist is removed to obtain QD2 structures. (iv) E-beam resist is applied again. (v) Region where QD1 structure is to be formed is drawn by e-beam lithography. (vi) QD1 solution is dropped onto resist and allowed to dry. (vii) E-beam resist is removed to obtain the sample. (b) Fluorescence microscope image of sample excited by mercury vapor lamp.

3.2. Experimental Setup

The experimental setup used for performing measurements is shown in Figure 3. For the excitation light , we used the 2nd harmonic wave of a mode-locked Ti : sapphire laser (Mira 900, Coherent, Inc.) with a photon energy  eV, a pulse width of 2 ps, and a repetition frequency of 80 MHz, and for the control light , we used a CW diode-pumped solid state (DPSS) laser with a photon energy  eV. The incident powers of and were 5 μW and 10 μW, respectively, and the beam spot diameters were both 5 μm. The beams from the two lasers were shaped by passing them through pinholes with 5 μm diameter apertures; then, using a half-mirror, the beams irradiated the sample from the back surface of the substrate via an objective lens with a numerical aperture NA = 0.4. As described in the previous section, because it is possible to observe the change in excitation transfer at the output end as a variation in the emitted light from QD1, we evaluated this in order to make it easier to separate it from the control light. The light emitted from QD1 was observed using an electron multiplying CCD camera (Hamamatsu ImagEM C9100-13H, Hamamatsu Photonics K.K.) after passing through an objective lens with NA = 0.55 and two bandpass filters with transmission wavelengths of  nm, which is close to the wavelength of . The acquired images were 512 × 512 pixel, 16-bit grayscale images, and the resolution was 0.37 μm/pixel.

Figure 3: Schematic diagram of experimental setup for measurement.

As an example of the measurement results, Figure 4(a) shows a CCD image acquired when O and A were irradiated with excitation light and control light , respectively. The green and yellow broken lines show the QD1 and QD2 structures, respectively. Point O was irradiated with , and we observed that the light emitted from the QD1 structure spread out in three directions due to transfer of the excitation. A cross-sectional brightness profile taken along the light-blue dotted line in this figure is shown by the red solid line in Figure 4(b). The black solid line shows the case where point O was irradiated with only . The red solid line is higher in a portion of the QD1 structure serving as the transfer path OB and the point B having the QD2 structure, confirming the tendency for the optical excitation directed towards point B to be increased due to the influence of incident on point A.

(a)

(b)

(c)

Figure 4: (a) Image of light emission from E11 in sample, acquired with CCD camera. Points O and A are irradiated with light and , respectively. Green and yellow broken lines represent the QD1 and QD2 structures, respectively. (b) Cross-sectional profile of emission intensity from E11, acquired with CCD. The red solid line is the cross-sectional profile taken along the light-blue dotted line in (a). The black solid line is the cross-sectional profile taken along the same location, when point O was irradiated with only. The 0 μm base point and the blue broken line indicate point O and point B, respectively. (c) Graph of relative change in emission intensity from E11. Blue, red, and green bars represent Cases  1,  2, and  3, respectively.

3.3. Energy Transfer in Y-Junction Structure

To evaluate the intensity of the emitted light, we used the brightness values acquired with the CCD camera, where the total brightness of 3 × 3 pixels was defined as emission intensity, . While irradiating point O with , we measured , where indicates the output end A or B, and indicates where the control-light is irradiated, that is, to the input terminal A, B, or C. Also, when no control light is irradiated, we represent this by . Concerning the characteristics of the device under study, described in Section 3.1, there are four representative categories:

(a)Reference: without any ;(b)Case 1: emission end and control-light input terminal are the same point ;(c)Case 2: emission end and control-light input terminal are opposite points (() = (A, B) or (B, A));(d)Case 3: control-light input terminal is point C ().

The evaluation was performed using the relative intensity change, , defined as . At this time, if is positive in Case 2, excitation transfer from one end at which the control light is input to the other end will be increased, and control of the excitation transfer will be observed.

Figure 4(c) shows a plot of obtained in Cases 1, 2, and 3. First, we evaluated Case 1 (i.e., ), as shown by the blue bars in Figure 4(c); they took positive values and , confirming that the emitted light was stronger. This is the effect of controlling the excitation transfer, and regarding energy transfer in the regions where QD1 and QD2 coexist, we obtained results showing that it is possible to control the excitation transfer also among multiple QDs. This indicates that an AND-gate operation was manifested with multiple randomly distributed QDs, as reported in [3].

Next, we evaluated Case 2 (i.e., = A, B), as shown by the red bars in Figure 4(c); they also took positive values and , confirming that the emitted light was stronger. was higher than because the wider overlapped area of QD1 and QD2 structures at point B provided more variation of excitation transfer than that at point A. This difference can be reduced by preparing finer structures and will not disturb the proper operation of solution exploring algorithms. From these results, we confirmed that the excitation transfer was suppressed by the incident control light. These results show that it is possible to control the optical excitation transfer dynamics in a randomly distributed QD structure. Since the relative change was reduced by losses involved in the optical excitation transfer, we expect that larger values will be possible by designing a structure having shorter distances OA and OB.

3.4. Discussion on Negative Intensity Change

Moreover, we evaluated Case 3 (i.e. ), as shown by the green bars in Figure 4(c); they took negative values and . Ideally, QD1 should be transparent to ; however, in our experiments, defect levels with energies lower than E11 in QD1 were excited by [15, 16]. We consider that this excitation prevented the excitation in E11 from slow and nonradiative relaxation via the defect levels.

To verify this, we performed photoluminescence measurements on the QD1 structure. In this experiment, we used samples having only the QD1 structure, fabricated with the same process as that used to fabricate the measurement samples, and we performed measurements with an experimental setup in which the bandpass filters and the CCD camera in the setup shown in Figure 3 were replaced with a spectrometer. The incident powers of and were 5 μW and 770 μW, respectively, and the beams irradiated the sample in spot diameters of 10 μm. The emission intensities were integrated in the range 2.25 < < 2.29 eV, corresponding to the peak wavelength of , and were normalized to the emission intensity when only was radiated and were compared. As a result, the emission intensity for the case where and were simultaneously radiated was 1.09, showing that the emitted light was increased by radiating (Figure 5(a)). This was the result of defect levels in QD1 being excited by , suppressing the nonradiative relaxation process from E11 via defect levels, which increased the emitted light (Figure 5(b)).

(a)

(b)

(c)

(d)

Figure 5: (a) Normalized intensity of emitted light from E11, showing cases where the QD1 structure was irradiated with only and with both and . (b) Schematic image of emission from QD1 with only or with both and . (c), (d) Energy diagram schematically showing the experimental sample from point C to point A: (c) without , (d) point C irradiated with .

From this result, the negative values in Case 3 can be explained as follows. Without , the excitation transferred to end C reflects and enhances at ends A and B (Figure 5(c)). By irradiating end C with , at C increases due to excitation of defect levels in QD1, and the reflection decreases (Figure 5(d)). As a result, at ends A and B decreases by irradiation of C with . This differs from the principle of excitation transfer suppression described in Section 2 and Figure 1(e) but indicates the possibility of another method of controlling the spatiotemporal dynamics of excitation transfer. Also, this increased light emission due to the fact that excitation of defect levels can also occur in QD1 at points A and B; however, the energy transfer suppression effect due to excitation of E21 should be larger than this because absorption of ( eV) in QD2 is 24 times higher than that in QD1. Thus, the effect of excitation of defect levels at ends A and B is negligible.

4. Summary

In summary, to show the possibility of controlling the spatiotemporal dynamics of optical excitation transfer based on dressed photon interactions between multiple randomly distributed quantum dots (QDs), we conducted experiments to control excitation transfer using randomly distributed QD structures. Using CdSe/ZnS core-shell QDs, we fabricated a Y-junction structure composed of randomly distributed QDs (QD1 structure). Optical excitation incident at the center of this structure was transferred to a QD2 structure serving as an output end located 3–7 μm away, and by irradiating a QD2 structure at another end with control light , we observed a maximum increase of . Thus, we have shown, for the first time, that it is possible to control the spatiotemporal dynamics of optical excitation transfer in a randomly distributed QD structure. Our findings will lead to simplified implementation and driving of solution searching and decision making devices based on the optical excitation transfer dynamics between QDs and will contribute to their practical realization. On the other hand, when control light was incident on the output end that did not have a QD2 structure, we found that the amount of excitation energy transferred to the output end was reduced by a maximum of . This was thought to be because the light emission from QD1 was increased due to excitation of defect levels in QD1, and the amount of excitation energy transferred to the end provided with the QD2 structure serving as an output was reduced. This effect acted in a direction causing a reduction in excitation transfer, opposite to the increased transfer level due to excitation of the QD2 structure, indicating the possibility of a different kind of control of the spatiotemporal dynamics. This result is expected to lead to the development of novel computing device architectures in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

A part of this work was supported by the “Development of Next-generation High-performance Technology for Photovoltaic Power Generation System” Program, NEDO, Japan.

References

M. Ohtsu, T. Kawazoe, T. Yatsui, and M. Naruse, “Nanophotonics: application of dressed photons to novel photonic devices and systems,” IEEE Journal on Selected Topics in Quantum Electronics, vol. 14, no. 6, pp. 1404–1417, 2008.

T. Kawazoe, K. Kobayashi, S. Sangu, and M. Ohtsu, “Demonstration of a nanophotonic switching operation by optical near-field energy transfer,” Applied Physics Letters, vol. 82, no. 18, pp. 2957–2959, 2003.

T. Kawazoe, M. Ohtsu, S. Aso et al., “Two-dimensional array of room-temperature nanophotonic logic gates using InAs quantum dots in mesa structures,” Applied Physics B: Lasers and Optics, vol. 103, no. 3, pp. 537–546, 2011.

A. Karimkhani and M. K. Moravvej-Farshi, “Design of three-input nanophotonic and gates,” Journal of the Optical Society of America B: Optical Physics, vol. 26, no. 5, pp. 1084–1090, 2009.

T. Kawazoe, K. Kobayashi, and M. Ohtsu, “Optical nanofountain: a biomimetic device that concentrates optical energy in a nanometric region,” Applied Physics Letters, vol. 86, no. 10, Article ID 103102, pp. 1–3, 2005.

W. Nomura, T. Yatsui, T. Kawazoe, M. Naruse, and M. Ohtsu, “Structural dependency of optical excitation transfer via optical near-field interactions between semiconductor quantum dots,” Applied Physics B: Lasers and Optics, vol. 100, pp. 181–187, 2010.

P. Holmström, L. Thylén, and A. Bratkovsky, “Composite metal/quantum-dot nanoparticle-array waveguides with compensated loss,” Applied Physics Letters, vol. 97, Article ID 073110, 2010.

M. Naruse, M. Aono, S.-J. Kim et al., “Spatiotemporal dynamics in optical energy transfer on the nanoscale and its application to constraint satisfaction problems,” Physical Review B, vol. 86, Article ID 125407, 2012.

M. Aono, M. Naruse, S.-J. Kim et al., “Amoeba-inspired nanoarchitectonic computing: solving intractable computational problems using nanoscale photoexcitation transfer dynamics,” Langmuir, vol. 29, no. 254, pp. 7557–7564, 2013.

S.-J. Kim, M. Naruse, M. Aono, M. Ohtsu, and M. Hara, “Decision maker based on nanoscale photo-excitation transfer,” Scientific Reports, vol. 3, Article ID 2370, 2013.

M. Naruse, H. Hori, K. Kobayashi, P. Holmström, L. Thylén, and M. Ohtsu, “Lower bound of energy dissipation in optical excitation transfer via optical near-field interactions,” Optics Express, vol. 18, no. 23, pp. A544–A553, 2010.

U. Schöning, “A probabilistic algorithm for k-SAT and constraint satisfaction problems,” in Proceedings of the 40th IEEE Symposium on Foundations of Computer Science, pp. 410–414, 1999.

N. Sakakura and Y. Masumoto, “Persistent spectral-hole-burning spectroscopy of CuCl quantum cubes,” Physical Review B—Condensed Matter and Materials Physics, vol. 56, no. 7, pp. 4051–4055, 1997.

C. Trallero-Giner, A. Debernardi, M. Cardona, E. Menéndez-Proupín, and A. I. Ekimov, “Optical vibrons in CdSe dots and dispersion relation of the bulk material,” Physical Review B—Condensed Matter and Materials Physics, vol. 57, no. 8, pp. 4664–4669, 1998.

M. Y. Valakh, Y. G. Sadfyev, N. O. Korsunska et al., “Deep-level defects in CdSe/ZnSe QDs and giant anti-stokes photoluminescence,” Semiconductor Physics, Quantum Electronics and Optoelectronics, vol. 5, no. 3, pp. 254–257, 2002.

M. Danek, K. F. Jensen, C. B. Murray, and M. G. Bawendi, “Synthesis of luminescent thin-film CdSe/ZnSe quantum dot composites using CdSe quantum dots passivated with an overlayer of ZnSe,” Chemistry of Materials, vol. 8, no. 1, pp. 173–180, 1996.