Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Interferometric tunability of absorption

Open Access Open Access

Abstract

We propose an interferometric setup that permits to tune the quantity of radiation absorbed by an object illuminated by a fixed light source. The method can be used to selectively irradiate portions of an object based on their transmissivities or to accurately estimate the transmissivities from rough absorption measurements.

©2006 Optical Society of America

1. Introduction

When an object is illuminated, it will absorb radiation proportionally to its absorption coefficient: Darker portions of the object will absorb more light than more transparent ones. Is there a way around this? In this paper we analyze a setup which uses classical light sources (i.e. coherent beams) and permits to easily tune the quantity of a light absorbed by an object independently on its transparency, by appropriately tuning an interferometer phase. With the same setup, high efficiency measurements of the absorption coefficient can be performed via a feedback mechanism. The only underlying assumption is that the object introduces a negligible dephasing into a probe beam. Since we can employ quasi-monochromatic light, this assumption is met in a variety of systems. Moreover, in the case of objects that have a homogeneous phase image, the dephasing can be easily compensated with the interferometer phase.

Our proposal draws inspiration from the so called “interaction-free-measurement” setups, where a partially transparent object can be discriminated from a totally transparent one with asymptotically negligible radiation absorption [1, 2, 3, 4, 5]. Even though such proposals were originally based on single-photon light pulses, analogous results have been obtained also with classical light [6, 7].

The layout of the paper follows. We start by describing the proposed interferometric setup. We show how the absorption peak can be tuned and we analyze the irradiation selectivity. We then give the protocol for high efficiency estimation of η. We conclude by analyzing inhomogeneous objects, which incorporate different transmissivities. Since the process does not involve any quantum effects (such as entanglement or squeezing) one could also analyze it in terms of a classical theory of radiation, instead of the quantum formalism we use here for rigor.

2. The apparatus

The proposed apparatus is a modification of the experimental setup of Ref. [7]. It is obtained by concatenating a collection of N Mach-Zehnder (MZ) interferometers and is depicted in Fig. 1. Initially a coherent state |α〉 enters through one interferometer port (associated with the annihilation operator a 0), and no photons enter from the other port (associated with the annihilation operator b 0). As shown in Fig. 1, after each MZ, one of the two emerging beams (the R beam) is focused on the object. Then the two beams are recombined at the input port of the next MZ. After N of such steps, the radiation leaves the apparatus at the Nth interferometer outputs aN and bN . As we will show, appropriately tuning the interferometers phase ϕ and the number N of MZs it is possible to choose the value of the transmissivity η that will absorb the most radiation in the object. The transmissivity η of an object is the probability that a single photon will pass through it or, equivalently, the percentage of the transmitted intensity of an impinging coherent beam. Note that an apparatus employing a single MZ which is crossed N times by the light can also be employed, where N can be controlled by appropriately tilting one of the interferometer mirrors [7] or by using an acousto-optics switch.

The input-output relations of the interferometer can be obtained observing that when two coherent states |αn 〉 and |βn 〉 impinge, respectively, into the input ports an and bn of the n-th MZ (see Fig. 2), the corresponding outputs at ports a n+1 and b n+1 are still coherent states of amplitudes α n+1 and β n+1, given by

(αn+1βn+1)=S(αnβn),

with

S=eiϕ2(cos(ϕ2)isin(ϕ2)isin(ϕ2)cos(ϕ2)),

where ϕ is the interferometer phase. The output an is directly fed into a port of the MZ number n+1, while the output bn first passes through the object and then enters the other port of the same MZ. Since the object absorbs each photon with a probability η without introducing any phase factor, its action on the input coherent state |β 0+1〉 can be modeled as a beam splitter with transmissivity η that couples the input radiation to a vacuum state and then discards one of the two outputs. As a result the state |β n+1〉 is transformed into a coherent state of reduced amplitude √ηβ n+1 [8]. Thus, in the presence of the absorber, the amplitudes α n+1 and β n+1 of the coherent states at the input of the MZ interferometer number n+1 is given by

(αn+1βn+1)=S(η)(αnβn),

where

S(η)=eiϕ2(cos(ϕ2)isin(ϕ2)iηsin(ϕ2)ηcos(ϕ2)).
 figure: Fig. 1.

Fig. 1. Proposed apparatus. It consists of N Mach-Zehnder (MZ) interferometers concatenated so that the output ports of the nth MZ is fed into the input ports of the successive one (for the sake of clarity the first interferometer is graphically enhanced). All interferometers act on the radiation with the same phase shift ϕ. The object to be irradiated is placed outside the MZs and it interacts only with the R beams. Initially the radiation enters from the input a 0. After N round trips, it exits through the outputs aN and bN .

Download Full Size | PDF

Iterating Eq. (3) N times we can express the amplitude of the coherent states emerging from the whole apparatus as

(αNβN)=SN(η)(α00).

Some examples of such evolution are given in Fig. 3, and an analytic solution can be obtained by diagonalizing S(η) [7]. The light absorbed by the object is given by

Iab=α02(αN2+βN2)rα02,

i.e. the input intensity |α 0|2 minus the total output intensity |αN |2+|βN |2. The quantity r is a complicated function of N, ϕ and η which can be explicitly computed from Eq. (5). It measures the “effective” absorption constant of the object.

3. Discussion

The possibility of changing the absorption of the illuminated object from its natural value 1-η to an effective value r≃0 allows one to determine the presence of a completely opaque object (i.e. η=0) with only an asymptotically small fraction of the input radiation being absorbed [7, 3, 4, 5, 6]. Consider Eq. (5) for η=1 (e.g. completely transparent object) and η=0 (e.g. completely opaque object). In these cases simple analytical solutions can be obtained yielding

αN=α0eiNϕ2cos(Nϕ2)
βN=iα0eiNϕ2sin(Nϕ2),

for η=1, and

αN=α0eiNϕ2cosN(ϕ2)
βN=0,

for η=0. By choosing ϕ=π/N, from Eqs. (8) and (10), it is immediate to see that all radiation exits from the bN -port if η=1 and that most of the radiation (asymptotically all of it for N→∞) exits from the aN -port if η=0 [9]. In both cases the light absorption is minimal (i.e. exactly null in the first case and asymptotically null in the second one). Nonetheless they can be discriminated by simply looking from which interferometer ports (e.g. aN or bN ) the light emerges.

 figure: Fig. 2.

Fig. 2. Mach-Zehnder interferometer constituting the nth element in the Mach-Zehnder sequence in the apparatus of Fig 1. The first of the two 50–50 beam splitters transforms the input annihilation operators an, bn into c =(an+bn )/√2 and d=(bn-an )/√2 respectively. The second beam splitter transforms the annihilation operators c and d into a n+1=(c+d)/√2 and b n+1=(d-c)/√2.

Download Full Size | PDF

The possibility of controlling the effective absorption r of the object by changing the interferometer parameters is evident from Fig. 4 where we plot r as a function of the transmissivity η for different values of ϕ(choosing again N=π/ϕ): The function r exhibits a peak for η=ηmax which increases from ηmax ~0 to ηmax ~1 as ϕ decreases. This effect can be explained intuitively as follows. For small values of ϕ(i.e. high values of N) little radiation is leaked into the R modes at every round trip with the exception of the case when η is high. On the contrary, for small values of N (i.e. large values of ϕ) a larger amount of radiation is leaked into the R modes at every round trip, so that the absorption peak moves to lower values of η. The dependence of the absorption peak maximum as a function of ϕ and N is depicted in Fig. 5, left. This graph also shows the values of ηmax that can be attained in practice: it can be accurately fine-tuned only for ηmax ≲0.5 since only few discrete values of ηmax are achievable for low N, whereas high ηmax ~1 requires large N, which can be difficult to achieve practically. Finally, the selectivity of the absorption, i.e. the width of the effective absorption curve r as a function of η (see Fig. 4), is not constant when ϕ is varied: The value of the width-at-half-maximum is smaller for absorption curves peaked at ηmax ≃0,1 and larger for ηmax ≃0.5. In the limit ϕ→ 0 the r-curve becomes a very narrow spike peaked just below η=1.

 figure: Fig. 3.

Fig. 3. Plot of the (rescaled) output amplitudes of the nth MZ interferometer |αn /α 0|2 in the a-modes at L (continuous line) and |βn /α 0|2 in the b-modes at R (dashed line). Initially all the radiation is in mode a 0, but, as the evolutions progresses, more and more radiation is transferred to the b-modes, until (for n=π/ϕ) the radiation is entirely transferred. Here ϕ=π/10 so that the total transfer occurs for n=10 (vertical line). Left: The object is completely transparent (η=1), so that the total energy (dotted line) is constant; Right the object is semi-transparent (η=.9), so that the total energy decreases as the evolution progresses.

Download Full Size | PDF

4. High precision η-measurements

Our scheme can be easily adapted to high precision estimation of the absorption coefficient, starting from low-quality measurements of the absorption r. [The main idea is revealed by the lower right graph of Fig. 4.] The required iterative procedure is composed by the following steps: i) start by roughly estimating η through an absorption measurement and set the interferometer phase so that the r curve has a steep slope corresponding to such value of η; ii) perform another absorption measurement and estimate a better value of η; iii) again tune the interferometer phase, and so on. Since the absorption curve r for the values of η~1 can be very steep, a very good estimate of these ηs can be achieved even when the measurement of r contains a large error Δr (see Fig. 4). Notice that the high values of η~1 are the hardest to estimate reliably without strongly irradiating the object, since they are associated to the region of least transparency.

5. Inhomogeneous samples

In deriving Eq. (3) we implicitly assumed that η is constant, i.e. that we employ a spatially homogeneous object with uniform absorption within the waist of the light beams crossing it. Instead, if it has spatially inhomogeneous absorption, we can still use Eq. (3) to describe the ab- sorption of “portions” of the incoming beam. In fact, in the limit in which the scale of the spatial inhomogeneities is much larger than the wavelength λ of the source light, the diffraction of the propagating beam induced by these inhomogeneities can be neglected. In this regime the illuminating beam can be effectively decomposed in independent “sub-beams” of the wavelength λ and interacting independently with the different portions of the object. A typical example is when one performs imaging of a macroscopic object. Controlling the interferometer parameters we can then selectively choose which portion will absorb the most of the input radiation by appropriately shifting the position of the absorption peak of r.

Explicitly, to selectively irradiate an inhomogeneous sample one thus needs to: i) estimate the transmissivity of the various portions of the sample by conventional imaging techniques, identifying the transmissivity η̄ of the region that needs irradiation; ii) tune the phase ϕ=π/N so that the absorption is maximized for η̄.

S. L. acknowledges financial support by ARDA, DARPA, ARO, AFOSR, NSF, and CMI; V. G. was in part supported by the EC under contract IST-SQUBIT2 and by the Quantum Information research program of Centro di Ricerca Matematica Ennio De Giorgi of Scuola Normale Superiore; L. M. acknowledges financial support by the Ministero Italiano dell’Università e della Ricerca (MIUR) through FIRB (bando 2001) and PRIN 2005.

 figure: Fig. 4.

Fig. 4. Above left: Plot of the (rescaled) apparatus output amplitudes |αN0 |2 at the output aN (continuous line) and |βN0 |2 at bN (dashed line) as a function of the transmissivity η of the object with ϕ=π/60. For η=1 (total transparency) all the output radiation is at bN , whereas for η=0 (total absorption) all the output radiation is at aN and a small amount of radiation has been absorbed. The absorbed radiation is proportional to the effective absorption constant r, defined in Eq. (6), which is depicted by the dotted lines. Above right: Plot of the absorption r only, in the case ϕ=π/200. This plot illustrates how the setup can be used to measure the absorption coefficient η by measuring r, attaining a higher precision than with a direct measure (in which case the absorption is given by the dashed line): Starting from the same uncertainty Δr in the measurement of r with our procedure we can obtain a lower uncertainty Δη in the measurement of η than the one, Δη, obtainable from a direct measurement. Below: Plot of r as a function of η and N=π/ϕ. Notice how the absorption peak shifts as a function of ϕ. The absorption peak moves to higher η for decreasing ϕ.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Left: Maximum ηmaxN (circles) and average value ηav (squares) of the absorption peak r of Fig. 4 as a function of N=π/ϕ. (The maximum and the average follow different evolutions because of the asymmetry in the absorption curves). Increasing N (i.e. decreasing ϕ), the maximum in the absorption peak moves to higher values of η. The graph also details which are the actual values of η that can be achieved through the proposed setup as a function of N. The dotted line is the function [(N-1)/N]4 that gives a good interpolation of the peak evolution. Right: Selectivity in the irradiation as a function of the transmissivity peak. The width of the peaks of the dotted line in Fig. 4 is not uniform. Here we plot the Root Mean Square (stars) and the width at half maximum (circles) of the absorption curve as a function of the absorption curve maximum. Notice that the RMS is almost constant over the whole range.

Download Full Size | PDF

References and links

1. A. Elizur, L. Vaidman, and L. Vaidman, “Quantum mechanical interaction-free measurements,” Found. Phys.23, 987–997 (1993); “On the realization of interaction-free measurements,” Quantum Opt.6, 119 (1994). [CrossRef]  

2. P. Kwiat, H. Weinfurter, A. Zeilinger, P. Kwiat, H. Weinfurter, T. Herzog, A. Zeilinger, and M. A. Kasevich, “Quantum Seeing in the Dark,” Sci. Am.275, 72–77 (1996); “Interaction-Free Measurement,” Phys. Rev. Lett.74, 4763–4767 (1995). [CrossRef]  

3. A. G. White, J. R. Mitchell, O. Nairz, and P. G. Kwiat, “≪Interaction-free≫imaging,” Phys. Rev. A 58, 605–613 (1998). [CrossRef]  

4. P. Kwiat, “Experimental and theoretical progress in interaction-free measurements,” Physica Scripta T76, 115–121 (1998). [CrossRef]  

5. P. Facchi, Z. Hradil, G. Krenn, S. Pascazio, and J. Řeháček, “Quantum Zeno tomography,” Phys. Rev. A 66, 12110 (2002). [CrossRef]  

6. S. Inoue and G. Björk, “Experimental demonstration of exposure-free imaging and contrast amplification,” J. Opt. B: Quantum Semiclass. Opt. 2, 338–345 (2000). [CrossRef]  

7. J.-S. Jang, “Optical interaction-free measurement of semitransparent objects,” Phys. Rev. A 59, 2322–2329 (1999). [CrossRef]  

8. Rigorously speaking, the transformation βn+1→√ηβn+1 does not correspond to the linear mapping |βn+1〉→|√ηβn+1〉 but to its density matrix counterpart: |βn+1〉〈βn+1|→|√ηβn+1〉〈√ηβn+1|. The latter accounts for decoherence effect, whereas the former does not. Since in our analysis we are always dealing with factorized states, the two transformations coincide for us.

9. This effect can be explained intuitively as follows. If the object is transparent, at the first round trip a small amount of radiation leaks into the R modes, at the second round trip a higher amount leaks there and so on constantly increasing through constructive interference until all the radiation moves into such modes after N=π/ϕ round trips. If, instead, the object is opaque, the little radiation that has leaked into the R modes at the first round trip is absorbed and does not contribute to the constructive interference that would draw more radiation into these modes at the second round trip. As a result very little radiation transfers and most of it remains in the L modes.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Proposed apparatus. It consists of N Mach-Zehnder (MZ) interferometers concatenated so that the output ports of the nth MZ is fed into the input ports of the successive one (for the sake of clarity the first interferometer is graphically enhanced). All interferometers act on the radiation with the same phase shift ϕ. The object to be irradiated is placed outside the MZs and it interacts only with the R beams. Initially the radiation enters from the input a 0. After N round trips, it exits through the outputs aN and bN .
Fig. 2.
Fig. 2. Mach-Zehnder interferometer constituting the nth element in the Mach-Zehnder sequence in the apparatus of Fig 1. The first of the two 50–50 beam splitters transforms the input annihilation operators an, bn into c =(an+bn )/√2 and d=(bn-an )/√2 respectively. The second beam splitter transforms the annihilation operators c and d into a n+1=(c+d)/√2 and b n+1=(d-c)/√2.
Fig. 3.
Fig. 3. Plot of the (rescaled) output amplitudes of the nth MZ interferometer |αn /α 0|2 in the a-modes at L (continuous line) and |βn /α 0|2 in the b-modes at R (dashed line). Initially all the radiation is in mode a 0, but, as the evolutions progresses, more and more radiation is transferred to the b-modes, until (for n=π/ϕ) the radiation is entirely transferred. Here ϕ=π/10 so that the total transfer occurs for n=10 (vertical line). Left: The object is completely transparent (η=1), so that the total energy (dotted line) is constant; Right the object is semi-transparent (η=.9), so that the total energy decreases as the evolution progresses.
Fig. 4.
Fig. 4. Above left: Plot of the (rescaled) apparatus output amplitudes |αN0 |2 at the output aN (continuous line) and |βN0 |2 at bN (dashed line) as a function of the transmissivity η of the object with ϕ=π/60. For η=1 (total transparency) all the output radiation is at bN , whereas for η=0 (total absorption) all the output radiation is at aN and a small amount of radiation has been absorbed. The absorbed radiation is proportional to the effective absorption constant r, defined in Eq. (6), which is depicted by the dotted lines. Above right: Plot of the absorption r only, in the case ϕ=π/200. This plot illustrates how the setup can be used to measure the absorption coefficient η by measuring r, attaining a higher precision than with a direct measure (in which case the absorption is given by the dashed line): Starting from the same uncertainty Δr in the measurement of r with our procedure we can obtain a lower uncertainty Δη in the measurement of η than the one, Δη, obtainable from a direct measurement. Below: Plot of r as a function of η and N=π/ϕ. Notice how the absorption peak shifts as a function of ϕ. The absorption peak moves to higher η for decreasing ϕ.
Fig. 5.
Fig. 5. Left: Maximum ηmaxN (circles) and average value ηav (squares) of the absorption peak r of Fig. 4 as a function of N=π/ϕ. (The maximum and the average follow different evolutions because of the asymmetry in the absorption curves). Increasing N (i.e. decreasing ϕ), the maximum in the absorption peak moves to higher values of η. The graph also details which are the actual values of η that can be achieved through the proposed setup as a function of N. The dotted line is the function [(N-1)/N]4 that gives a good interpolation of the peak evolution. Right: Selectivity in the irradiation as a function of the transmissivity peak. The width of the peaks of the dotted line in Fig. 4 is not uniform. Here we plot the Root Mean Square (stars) and the width at half maximum (circles) of the absorption curve as a function of the absorption curve maximum. Notice that the RMS is almost constant over the whole range.

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

( α n + 1 β n + 1 ) = S ( α n β n ) ,
S = e i ϕ 2 ( cos ( ϕ 2 ) i sin ( ϕ 2 ) i sin ( ϕ 2 ) cos ( ϕ 2 ) ) ,
( α n + 1 β n + 1 ) = S ( η ) ( α n β n ) ,
S ( η ) = e i ϕ 2 ( cos ( ϕ 2 ) i sin ( ϕ 2 ) i η sin ( ϕ 2 ) η cos ( ϕ 2 ) ) .
( α N β N ) = S N ( η ) ( α 0 0 ) .
I ab = α 0 2 ( α N 2 + β N 2 ) r α 0 2 ,
α N = α 0 e iN ϕ 2 cos ( N ϕ 2 )
β N = i α 0 e iN ϕ 2 sin ( N ϕ 2 ) ,
α N = α 0 e iN ϕ 2 cos N ( ϕ 2 )
β N = 0 ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.