Division of Advanced Electronics and Optical Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
Abstract
We estimate the phase velocity of a modulation microwave in a
quasi-velocity-matched (QVM) electro-optic (EO) phase modulator
(QVM-EOM) using EO sampling which is accurate and the most reliable
technique for measuring voltage waveforms at an electrode. The
substrate of the measured QVM-EOM is a stoichiometric periodically
domain-inverted LiTaO3 crystal. The electric field of a standing wave
in a resonant microstrip line (width: 0.5 mm, height: 0.5 mm) is measured by
employing a CdTe crystal as an EO sensor. The wavelength of the traveling
microwave at 16.0801 GHz is determined as 3.33 mm by fitting the theoretical
curve to the measured electric field distribution. The phase velocity is
estimated as vm=5.35×107 m/s, though there exists
about 5% systematic error due to the perturbation by the EO sensor.
Relative dielectric constant of εr=41.5 is led as the maximum
likelihood value that derives the estimated phase velocity.
1. Introduction
Wide expansion of an optical spectrum by a deep phase
modulation at a high-modulation frequency is the essential technique for many
electro-optic (EO) light controlling fields, such as ultrashort pulse
generation [1], comb
generation [2],
time-to-space mapping of an optical signal [3], and so on.
Quasi-velocity-matching (QVM) with periodic domain
inversion of an LiTaO3 (LT) or LiNbO3 (LN) crystal in a traveling-wave EOM is one of the promising
techniques to widely spread the optical spectrum. Especially, periodically
domain-inverted stoichiometric LT (SLT) crystals are attractive QVM-EOMs
because of their stronger EO effect and lower coercive fields than
congruent-melt LT crystals. On the basis of the QVM technique, Morimoto et al.
spread a 3-THz-wide modulation sidebands (comb) at a 16.25-GHz modulation
frequency [4]. The QVM
technique compensates for velocity mismatching between an optical group
velocity
and a microwave phase velocity
in the EO crystal by periodically inverting
the sign of the EO coefficient of the crystal. In the ideal QVM-EOM, the
modulation index is almost proportional to the interaction length and as a
result, large modulation index at the high modulation frequency can be achieved
with long device.
Accurate values of
and
are needed in designing the QVM-EOM, because
uncertainty of these values results in degradation of the modulation
efficiency, which is achievable modulation index per unit interaction length
per unit modulation electric field strength. The group velocity of the optical
wave can be calculated with group refractive index
of the EO crystal. Likewise, the phase
velocity of the modulation microwave is calculated as
,
where
is the effective relative dielectric constant
which is determined by the geometry of the transmission line and the relative
dielectric constant
of the substrate. The group refractive index
of SLT crystals is estimated accurately using Sellmeier equation with constants
listed in [5]. In
contrast with the accurate value of
,
very little data of
for LT crystal are available, which is in the
range
at microwave frequencies [6, 7]. To our knowledge, no
accurate data of
for SLT crystals at microwave frequencies are
available.
In this paper, we estimate the phase velocity of the
modulation microwave in SLT-based QVM-EOM on the basis of a direct observation
of an amplitude of a standing-wave in a resonant microstrip line using EO sampling
technique [8]. Because
the EO sampling is accurate and most reliable technique for measuring voltage
waveforms at an electrode, the phase velocity can be estimated accurately from
the relation of
,
where
is the wavelength of the traveling microwave
and
is the resonant frequency. In Section 2,
the principle of the QVM is summarized. Section 3 discusses influence that the
estimation error of the phase velocity gives to the modulation efficiency. In
Section 4, experimental results are presented.
2. Quasi-Velocity-Matching
Figures 1(a) and 1(b) show a schematic of a normal traveling-wave
EOM and a QVM-EOM, respectively. To avoid unnecessary complexities, we employ a
one-dimensional analysis for our devices. For a traveling-wave EOM, if the
axis is the direction in which the optical and
modulation wave propagate in a virgin (single-domain) EO crystal, the variation
of the refractive index induced by an electric field
is obtained at
as
(1)where
is the amplitude of index changes,
is the extraordinary refractive index of the
crystal,
is the EO coefficient of the crystal, and
.
A half-period
of the domain inversion is given
by
(2) where
is the phase velocity of the modulation wave,
and
is the group velocity of the optical wave.
Here, we assume that the optical wave arrives at point
at time
.
Figure 1: (a) Traveling-wave EOM, (b) QVM-EOM with periodic domain inversion.
For
,
the phase retardation of optical wave at the position
is expressed as
(3)where
(4)If there is the
so-called velocity mismatching between the modulation wave and the optical
wave, the modulation index
becomes a periodical function of
with period
.
When a traveling-wave EOM has a suitable
domain-inverted half-period of
,
QVM occurs and accordingly a large modulation index is achieved. In such a
situation, phase retardation
given to the optical wave passing through the
length of
,
,
in a periodically domain-inverted crystal as shown in Figure 1(b) is expressed
as
(5)
The modulation index
can be expressed by
(6)
Figure 2 shows modulation indices of (a) the QVM-EOM,
and (b) the typical traveling-wave EOM (velocity mismatching). The solid line
of Figure 2(a) is calculated using (6). The modulation index achieved by
QVM-EOM is almost proportional to the interaction length though it is lower
than perfect velocity-matched condition by a factor of 2/
.
The upper limit of the modulation index for non-domain-inverted EOM is
.
Figure 2: Modulation indices of EOMs. (a) QVM-EOM with periodic domain inversion, and (b) typical
traveling-wave EOM (velocity mismatching).
3. Efficiency Analysis for the QVM-EOM
The analysis presented in the former section was ideal
case in which there is no error in the period of the domain inversion.
Uncertainties in the value of the phase velocity and group velocity degrade the
modulation efficiency through the length error of the
domain-inversion period for the QVM. In this
section, we will discuss degradation of the modulation efficiency in the
presence of length error in the period of the domain inversion.
Figure 3 show the refractive index changes seen from
the optical wave in periodically domain-inverted EO crystal. The half period of
the domain inversion in Figures 3(a) and 3(b) is
and
,
where
,
respectively. The QVM occurs in the case of Figure 3(a), and as a result
the traveling optical wave of
sees the positive refractive index change
throughout the interaction. The case of 10% error in the domain-inversion length is shown
in Figure 3(b). If there is an error in length of domain inversion, the
optical wave of
sees not only positive but also a negative
refractive index change throughout the interaction.
Figure 3: Refractive index change seen from the optical wave in periodically domain-inverted EO
crystals. (a) Perfectly quasi-velocity-matched condition. Half period of the
domain inversion is

.
(b) Half period of the domain inversion is

,
where

.
Figure 4(a) shows the normalized modulation index
calculated for the interaction length of
.
The modulation efficiency almost proportionally decreases with the increase of
.
The modulation efficiency decreases to about 50% in the case of
.
Figure 4(b) shows the modulation index as a function of the interaction
length. The calculation is carried out for the length error in the half period
of the domain inversion of (i)
,
(ii)
,
and (iii)
.
An upper limit of the modulation index exists if
.
Because longer interaction length is essential for larger modulation index, it
is important to estimate the length of the half period of the domain inversion
in high accuracy.
Figure 4: (a) Modulation index as a function of

.
The interaction length is

.
(b) Modulation index as a function of the interaction length. (i)

,
(ii)

,
and (iii)

.
Propagation of uncertainties of
and
can be evaluated by
(7)
The influence on the modulation efficiency of
uncertainty of
is lower than half of the influence on the
modulation efficiency of uncertainty of
.
Moreover,
can be estimated more accurately than
by using Sellmeier equation with constants for
SLT crystal listed in [5]. Using Sellmeier equation,
at
nm is derived. From this reliable
value of
,
the group velocity is derived as
m/s.
The phase velocity of the modulation microwave can be
calculated as
.
The effective relative dielectric constant for microstrip line can be
calculated using Kobayashi's formula [9, 10]. Kobayashi's formula is claimed to predict dispersion
better than 0.6% in the range
and
,
where
and
are width and height of the microstrip line,
respectively. However, uncertainty of the estimated
is comparatively large, because the reported
value of
for LT crystal is inaccurate and ranging
[7]. If we employ values of
of the modulation microwave for 16.25 GHz can
be calculated as
m/s
m/s using Kobayashi's formula. In that case,
the half period of the domain inversion for QVM at 16.25 GHz is calculated as
2.7 mm
3.1 mm with the group velocity of
m/s. There is 14% of length uncertainty in the half period of
the domain inversion. From Figure 4(b),
results in the degradation of the modulation
efficiency to about 10%.
4. Estimation of the Phase Velocity of the Modulation Wave Using EO Sampling Technique
Figure 5 shows a schematic of our QVM-EOM. An SLT
crystal is used for an EO substrate. The modulation electrode structure is a
microstrip line. The width of the microstrip line is
mm and the height of the substrate is
mm. The strip line is open-terminated for
resonance.
Figure 5: A schematic of
our QVM-EOM. The substrate of the QVM-EOM is the stoichiometric LiTaO
3 crystal. The domain inversion is performed with the half period of

.
The modulation electrode structure is the microstrip line. The width of the
microstrip line is

mm. The height of the substrate is

mm.
Figure 6(a) shows an experimental setup. A pulsed
fiber laser (repetition frequency:
MHz) was used for probe pulses. The
repetition frequency of the pulsed laser and the frequency of the microwave
signal source are synchronized with each other. The frequency of the microwave
should be set to be
,
where
is an integer and
is an offset frequency. In the experiment, the
frequency of the microwave is set to be
GHz, where
and
kHz.
Figure 6: (a) Experimental setup. (b) Cross-section of the microstrip line under
the measurement.
Figure 6(b) shows a cross-section of the microstrip
line under the measurement. A CdTe crystal is used for the EO sensor. The
aperture size of the EO sensor is 3 mm
3 mm. The probe beam is focused on the CdTe
crystal by the object lens (
). The spot size is about 20
m. The probe beam is reflected by a dielectric
mirror attached to the EO sensor. Polarization of the reflected beam is
modulated by the electric field of the microwave. The polarization-modulated
beam propagates along the same path of the incident beam path. Two orthogonal
polarization components of the modulated beam are differentially detected by
two photodiodes. The detected signal component of
passes through a bandpass filter whose center
frequency is 100 kHz. The amplitude of the signal, which is proportional to the
amplitude of the electric field of the standing wave, is measured by a spectrum
analyzer.
In the experiment, we scanned the probe beam in the
-direction with 50
m step and measured electric field profile of
the resonant standing wave. Figure 7 shows an experimental result. The
theoretical curve of
(solid line) is fitted to the experimental
data. Fitting parameters are offset power
and
.
From the least square fitting,
,
and
are estimated as
dBm,
,
and
,
respectively. As a result, the wavelength of the traveling microwave is
estimated as
mm. The phase velocity at the modulation
frequency of
GHz is estimated as
m/s. From this phase velocity, relative
dielectric constant of
is derived using Kobayashi's formula. The half
period of the domain inversion at the modulation frequency of near 16 GHz is
derived as
mm with the phase velocity of
m/s and the group velocity of 1.24
m/s.
Figure 7: Magnitude of the voltage
standing wave on an open-terminated microstrip line
at

GHz.
The estimated phase velocity is slower than the
real-phase velocity because of the perturbation of the EO sensor. In our case,
this systematic error is estimated by numerical analysis (moment method) to be
5% or less. The 5% error in the phase velocity corresponds to the
domain-inversion length error of
which results in the degradation of the
modulation efficiency to about 50%.
The systematic error depends mainly on (1) the
difference of the dielectric constant between the device under test (DUT) and
the EO sensor, (2) the distance between the DUT and the EO sensor, and (3) the
thickness of the EO sensor. Using low-dielectric
constant materials such as a polymer [11] as a sensor, the perturbation can be reduced. By
loading a low-dielectric material between the DUT and the sensor to keep a
distance, the perturbation can also be reduced, though the sensitivity is
reduced at the same time. The use of a thinner sensor is an adequate plan to
reduce the perturbation, however it also reduces the sensitivity. For a
specific case, the tradeoff relationship between the perturbation strength and
the sensitivity has been analyzed in the literature [12].
By using SLT substrate of the EOM as the EO sensor, it
is also possible to measure the standing wave without using an external
superstrate sensor [13]. In this case, special geometry should be constructed
to transmit the probe beam between the microstrip line and the ground plane.
This type of EOS system is now under construction and the results including the
comparison with current technique will appear elsewhere.
5. Conclusion
We have estimated the phase velocity of the modulation
microwave of 16 GHz in the periodically domain-inverted SLT phase modulator.
The theoretical curve of the amplitude distribution of the standing wave in the
microstrip line was fitted to the experimental data measured by EO sampling
technique. From the curve fitting, the phase velocity of
m/s was estimated though there is about 5% systematic error due to the perturbation by
the EO sensor. Relative dielectric constant of
was led as the maximum likelihood value that
derives the estimated phase velocity.
Acknowledgments
The authors would like to thank Drs. H. Togo and N.
Kukutsu from NTT Microsystem Integration Laboratories for their cooperation and
support. They also thank Dr. T. Kobayashi for helpful discussions. This
research was partially supported by a grant from the Global COE Program,
“Center for Electronic Devices Innovation,” from the Ministry of
Education, Culture, Sports, Science, and Technology of Japan, and Grant-in-Aid
for Scientific Research on Priority Areas, 19023006, 2008.
References
- K. Shibuya, S. Hisatake, and T. Kobayashi, “Sub-picosecond optical pulse generation with repetition rate of several tens of gigahertz using electrooptic standing-wave phase grating with periodic domain inversion,” Japanese Journal of Applied Physics, vol. 44, no. 9A, pp. 6558–6563, 2005.
- S. Hisatake, Y. Nakase, K. Shibuya, and T. Kobayashi, “Generation of flat power-envelope terahertz-wide modulation sidebands from a continuous-wave laser based on an external electro-optic phase modulator,” Optics Letters, vol. 30, no. 7, pp. 777–779, 2005.
- S. Hisatake and T. Kobayashi, “Time-to-space mapping of a continuous light wave with picosecond time resolution based on an electrooptic beam deflection,” Optics Express, vol. 14, no. 26, pp. 12704–12711, 2006.
- A. Motimoto, M. Tamaru, Y. Matsuda, M. Arisawa, and T. Kobayashi, in 3 THz sideband generation by electro-optic modulation and the application to ultrashort pulse generation, p. 234, Institute of Electrical and Electronics Engineers, Piscataway, NJ, USA, July 1995.
- M. Nakamura, S. Higuchi, S. Takekawa, K. Terabe, Y. Furukawa, and K. Kitamura, “Refractive indices in undoped and MgO-doped near-stoichiometric crystals,” Japanese Journal of Applied Physics, vol. 41, no. 4, pp. L465–L467, 2002.
- M. V. Jacob, J. G. Hartnett, J. Mazierska, J. Krupka, and M. E. Tobar, “Lithium tantalate—a high permittivity dielectric material for microwave communication systems,” in Proceedings of IEEE Region 10 Annual International Conference on Convergent Technologies for Asia-Pacific Region (TENCON '03), vol. 4, pp. 1362–1366, Bangalore, India, October 2003.
- M. V. Jacob, J. G. Hartnett, J. Mazierska, V. Giordano, J. Krupka, and M. E. Tobar, “Temperature dependence of permittivity and loss tangent of lithium tantalate at microwave frequencies,” IEEE Transactions on Microwave Theory and Techniques, vol. 52, no. 2, pp. 536–541, 2004.
- J. A. Valdmanis, G. Mourou, and C. W. Gabel, “Picosecond electro-optic sampling system,” Applied Physics Letters, vol. 41, no. 3, pp. 211–212, 1982.
- M. Kobayashi and R. Terakado, “Accurately approximate formula of effective filling fraction for
microstrip line with isotropic substrate and its application to
the case with anisotropic substrate,” IEEE Transactions on Microwave Theory and Techniques, vol. 27, no. 9, pp. 776–778, 1979.
- M. Kobayashi, “A dispersion formula satisfying recent requirements in microstrip
CAD,” IEEE Transactions on Microwave Theory and Techniques, vol. 36, no. 8, pp. 1246–1250, 1988.
- M. Yaita and T. Nagatsuma, “Optical sampling of electrical signals in poled polymeric media,” IEICE Transactions on Electronics, vol. E76-C, no. 2, pp. 222–228, 1993.
- T. Nagatsuma, T. Shibata, E. Sano, and A. Iwata, “Subpicosecond sampling using a noncontact electro-optic probe,” Journal of Applied Physics, vol. 66, no. 9, pp. 4001–4009, 1989.
- B. H. Kolner and D. M. Bloom, “Direct electro-optic sampling of transmission-line signals propagating on a gaas substrate,” Electronics Letters, vol. 20, no. 20, pp. 818–819, 1984.