Arbitrary waveform synthesis by multiple harmonics generation and phasing in aperiodic optical superlattices

doi:10.1364/OE.17.016342

« Show journal navigation

Arbitrary waveform synthesis by multiple harmonics generation and phasing in aperiodic optical superlattices

Wei-Hsun Lin and A. H. Kung »View Author Affiliations

Optics Express, Vol. 17, Issue 18, pp. 16342-16351 (2009)
http://dx.doi.org/10.1364/OE.17.016342

View Full Text Article

Acrobat PDF (292 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Basic principle of the model
3. Design algorithm
4. Numerical examples
5. Discussion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (16)
Cited By
Figures (4)
Metrics

Abstract

The process of generating periodic optical waveforms includes the generation and phasing of several harmonics of a fundamental frequency. In this work, we show that simultaneous generation and phasing of the harmonics can be performed in a monolithic aperiodic optical superlattice (AOS). Stable periodic waveforms can thus be delivered to a predetermined location by simply sending a laser beam through a properly designed and fabricated AOS crystal. A detailed mathematical description for generating the domain pattern in such an AOS crystal is given and the process is numerically demonstrated. The waveform that is generated from a monolithic AOS is highly reproducible and phase stable. We also use propagation in air as an example to show how any predictable phase and amplitude modifications such as air dispersion that will alter the desired waveform can be pre-compensated in the design phase of the AOS crystal.

1. Introduction

In recent years there is a lot of interest in the synthesis of arbitrary optical waveforms because of their potential application in telecommunication [1

Z. Jiang, D. S. Seo, S.-D. Yang, D. E. Leaird, R. V. Roussev, C. Langrock, M. M. Fejer, and A. M. Weiner, “Four-User, 2.5-Gb/s, Spectrally Coded OCDMA System Demonstration Using Low-Power Nonlinear Processing,” J. Lightwave Technol. 23(1), 143–158 (2005). [CrossRef]

], quantum control [2

M. C. Stowe, F. C. Cruz, A. Marian, and J. Ye, “High resolution atomic coherent control via spectral phase manipulation of an optical frequency comb,” Phys. Rev. Lett. 96(15), 153001 (2006). [CrossRef] [PubMed]

], and ultrafast electronic processing [3

A. Baltuska, Th. Udem, M. Uiberacker, M. Hentschel, E. Goulielmakis, Ch. Gohle, R. Holzwarth, V. S. Yakovlev, A. Scrinzi, T. W. Hänsch, and F. Krausz, “Attosecond control of electronic processes by intense light fields,” Nature 421(6923), 611–615 (2003). [CrossRef] [PubMed]

]. It is now understood that there are two essential requirements toward achieving full control in the synthesis of periodic arbitrary optical waveforms. The first one is the generation of a broad spectrum consisting of a series of commensurate harmonic components. The generated frequency components serve as the building blocks of any desired waveform. The other requirement is to have the ability to manipulate the phase and amplitude of each harmonic to achieve the desired waveform. To obtain such a periodic waveform, the phase of each harmonic must be adjusted to satisfy the condition

φ_{n} = φ_{C E} + n φ_{m},

(1)

where

φ_{C E}

is the carrier-envelope phase and

φ_{m}

is a uniform phase difference between adjacent components. When both requirements are satisfied, an arbitrary periodic waveform can be produced by superposition of the components of the series [4

A. M. Weiner, J. P. Heritage, and E. M. Kirschner, “High-resolution femtosecond pulse shaping,” J. Opt. Soc. Am. B 5(8), 1563–1572 (1988). [CrossRef]

Generally these requirements are satisfied in consecutive steps. A phase coherent Fourier spectrum is first generated using methods such as self-phase modulation in nonlinear photonic crystal fiber [5

J. K. Ranka, R. S. Windeler, and A. J. Stentz, “Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm,” Opt. Lett. 25(1), 25–27 (2000). [CrossRef]

], high harmonic generation from ultrafast lasers [6

P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71(13), 1994–1997 (1993). [CrossRef] [PubMed]

], molecular modulation [7

S. E. Harris and A. V. Sokolov, “Subfemtosecond Pulse Generation by Molecular Modulation,” Phys. Rev. Lett. 81(14), 2894–2897 (1998). [CrossRef]

], or multiple harmonics generation [8

T. W. Hänsch, “A proposed sub-femtosecond pulse synthesizer using separate phase-locked laser oscillators,” Opt. Commun. 80(1), 71–75 (1990). [CrossRef]

I. V. Shutov and A. S. Chirkin, “Consecutive high-order harmonic generation and formation of subfemtosecond light pulses in aperiodical nonlinear photonic crystals,” Phys. Rev. A 78(1), 013827 (2008). [CrossRef]

]. In some of these methods periodic ultrashort pulse trains are formed at the location where the phase-coherent spectrum is produced. Yet material dispersion and propagation phase shifts distort the waveform after exiting the generator. Hence generally the generation is followed by some form of a phase compensation scheme that employs grating pairs, prisms, and/or spatial light modulators to adjust the phase and amplitude of the series components to arrive at the desired result [10

M. Y. Shverdin, D. R. Walker, D. D. Yavuz, G. Y. Yin, and S. E. Harris, “Generation of a single-cycle optical pulse,” Phys. Rev. Lett. 94(3), 033904 (2005). [CrossRef] [PubMed]

–12

Z. Jiang, C.-B. Huang, D. E. Leaird, and A. M. Weiner, “Optical arbitrary waveform processing of more than 100 spectral comb lines,” Nat. Photonics 1(8), 463–467 (2007). [CrossRef]

The optical Fourier synthesizer using multiple harmonics was first proposed by Hänsch [8

T. W. Hänsch, “A proposed sub-femtosecond pulse synthesizer using separate phase-locked laser oscillators,” Opt. Commun. 80(1), 71–75 (1990). [CrossRef]

], who phase-locked two lasers to construct a frequency comb, employing frequency mixing processes in nonlinear crystals. The frequency mixing process and the phase-locked lasers guarantee a stable phase difference between every pair of adjacent frequency components. Waveforms can then be synthesized subsequently by monitoring and shifting the phase of each frequency component. In this paper we describe an approach utilizing an aperiodic optical superlattice (AOS) to simultaneously generate and phase the spectrum such that periodic electromagnetic waveforms of any predetermined shape can be delivered to a specified location in a single setting.

It is well-known that quasi-phase-matched nonlinear interaction can be achieved by periodically modulating the second order nonlinearity

d^{(2)}

to create a nonlinear reciprocal lattice vector in a crystal that offsets the phase-mismatch of a nonlinear interaction in the medium [13

N. Bloembergen, Nonlinear Optics (World Scientific, 1996).

]. By introducing aperiodic or quasi-periodic modulation Zhu and associates first showed that two lattice vectors can be created in one crystal to provide simultaneous quasi-phase-matched second and third harmonic generation in a superlattice [14

S. N. Zhu, Y. Y. Zhu, and N. B. Ming, “Quasi–Phase-Matched Third-Harmonic Generation in a Quasi-Periodic Optical Superlattice,” Science 278(5339), 843–846 (1997). [CrossRef]

]. Recently, Shutov and Chirkin advanced the technique to show by simulation that several harmonics can be generated in an AOS to provide a broad comb of coherent frequencies [9

]. They proposed AOS designs that optimize the efficiency of each harmonic generation stage to produce up to 8 consecutive harmonics, the superposition of which gives pulses of ultimate duration at some fixed interaction length determined and limited by the dispersion properties of the AOS. However, additional steps are required to compensate any distortion to the pulse when delivering it to the location where it will be applied. Here we expand on the idea of using an AOS to propose and numerically demonstrate a device that allows the delivery of an ultrafast periodic waveform to a predetermined location. For example, a specific waveform can be delivered directly to a heterojunction of a nanoscale photoelectronic device where the ultrafast dynamics of the electron and hole carriers are to be characterized. This is accomplished by engineering the phase and amplitude of

d^{(2)}

of an AOS so that the AOS will provide the required comb of harmonics that Fourier synthesizes to the desired waveform at that location.

2. Basic principle of the model

We begin with the equation for second order nonlinear frequency conversion

\frac{\partial E_{r}}{\partial z} = - i \sum_{ω_{p} + ω_{q} = ω_{r}} \frac{ω_{r} d^{(2)} (z)}{c n_{r}} E_{p} E_{q} \exp (- i Δ k_{(r p q)} z),

(2)

where

E_{n = r, p, q}

is the complex amplitude

E_{n} (z) e^{i φ_{n}}

of the time-dependent electric field

E_{n} (z, t) = E_{n} (z) e^{i (ω_{n} t - k_{n} z + φ_{n})}

, z is the direction of propagation,

d^{(2)} (z)

is the second order nonlinear coefficient,

ω_{n}

and

k_{n} = ω_{n} n_{n} / c

are the frequency and wavevector of

E_{n}

n_{n}

is refractive index, and

Δ k_{(r p q)} = k_{p} + k_{q} - k_{r}

is the photon momentum mismatch.

Successful periodic waveform synthesis requires the formation of a commensurate series of harmonics that begins with the fundamental frequency

ω_{1}

and the nth component is

ω_{n} = n ω_{1}

. These harmonics can be generated in a series of processes described by Eq. (2). To calculate the phase of each harmonic at z, we express Eq. (2) as

\frac{\partial E_{r}}{\partial z} = - i \sum_{(r p q)} \frac{ω_{r}}{c n_{r}} {\tilde{d}}_{(r p q)} E_{p} E_{q},

(3)

where the RHS is the sum of all possible channels

ω_{p} + ω_{q} \to ω_{r}

which contribute to the generation of

E_{r}

, and

{\tilde{d}}_{(r p q)}

is defined as a channel specific effective second order nonlinear coefficient of the medium over a distance L,

{\tilde{d}}_{(r p q)} \equiv \tilde{d} (Δ k_{(r p q)}) = \frac{1}{L} \int d^{(2)} (z) \exp (- i Δ k_{(r p q)} z) d z .

(4)

In conventional quasi-phase-matched harmonic generation

d^{(2)} (z)

is designed to optimize the generation efficiency so that

{\tilde{d}}_{(r p q)}

is real. By constructing

{\tilde{d}}_{(r p q)}

to have a nonzero phase, the output

E_{r}

will carry an additional phase

φ_{{\tilde{d}}_{(r p q)}}

. Hence through the design of

d^{(2)} (z)

we can impose a phase that is required on each harmonic output to produce the designed waveform waveform.

The connection between

φ_{{\tilde{d}}_{(r p q)}}

and the phase of each harmonic is shown as follows. While there can be many possible combinations of p and q in Eq. (4) we shall simplify the situation by choosing

p = 1

and

q = r - 1

, that is

ω_{1} + ω_{r - 1} \to ω_{r}

, to be the dominating channels. With this choice, every harmonic is produced by sum mixing that includes

E_{1}

. The phase of the n th harmonic generated is given by:

\begin{array}{l} φ_{n}^{(f i n a l)} = φ_{n} + ϕ_{n} = - \frac{π}{2} + φ_{{\tilde{d}}_{(n (n - 1) 1)}} + φ_{n - 1} + φ_{1} + ϕ_{n}, \\ φ_{1}^{(f i n a l)} = φ_{1} + ϕ_{1}, \end{array}

(5)

where

φ_{n}

is the phase of

E_{n}

at the exit end of the crystal,

φ_{n}^{(f i n a l)}

is the required phase of

E_{n}

at the desired location and

ϕ_{n}

represents the additional phase shift modulo to 2π acquired by

E_{n}

during propagation from the crystal exit to the final destination. In a periodic waveform, the CEP is given by

φ_{C E} = 2 φ_{1} - φ_{2} = φ_{1} + φ_{2} - φ_{3} = \dots

, therefore if we design

d^{(2)} (z)

of a crystal such that

\frac{π}{2} - φ_{{\tilde{d}}_{(211)}} + 2 ϕ_{1} - ϕ_{2} = \frac{π}{2} - φ_{{\tilde{d}}_{(321)}} + ϕ_{1} + ϕ_{2} - ϕ_{3} = \dots,

(6)

then all the harmonics will be phase-coherent to give a CEP that equals

\frac{π}{2} - φ_{{\tilde{d}}_{(211)}} + 2 ϕ_{1} - ϕ_{2}

. Note that the CEP is not affected by the absolute phase of the fundamental frequency

φ_{1}

. As can be seen from Eq. (5), phase fluctuation of the input laser,

φ_{1} \to φ_{1} + Δ φ_{1}

, will only result in a phase shift,

n Δ φ_{1}

, for the nth component,

φ_{n}^{(f i n a l)} \to φ_{n}^{(f i n a l)} + n Δ φ_{1}

and then a time shift

Δ φ_{1} / ω_{1}

to the overall synthesized waveform. Since Δφ ₁ is less than 2π and ω ₁ is in the optical range, the shift is insignificant so long as the input consists of more than a few cycles.

3. Design algorithm

The design of an AOS that can deliver a desired waveform will amount to finding a solution for

d^{(2)} (z)

by solving numerically Eq. (2) with an initial input of

E_{1} = E (0)

and

E_{n} = 0

i \neq 1

and a given set of final E field amplitudes and phases as the boundary conditions. Shutov describes a generalized sign function

d^{(2)} (z) \propto sgn (\sum_{(r p q)}^{N} a_{(r p q)} \sin (Δ k_{(r p q)} z + ζ_{(r p q)})),

(7)

where

a_{(r p q)}

and

ζ_{(r p q)}

are amplitude and phase parameters of the harmonic modulation.

Δ k_{(r p q)}

in the sine function is to include the quasi-phase-matching of the channel

ω_{p} + ω_{q} \to ω_{r}

in the lattice periodicity to characterize the periodical dependence of

d^{(2)} (z)

for N simultaneous QPM processes. We take into consideration that there is a minimum domain length that can reliably be fabricated in the AOS and modify Eq. (7) to include this minimum length. We further add a criterion that the width of each domain is equal to an integral multiple of this minimum width in the superlattice to arrive at the following design function for the crystal:

d^{(2)} (z) \propto sgn (\int_{⌊ z / l - 1 ⌋ l}^{⌊ z / l ⌋ l} \sum_{(r p q)}^{N} a_{(r p q)} \cos (Δ k_{(r p q)} z' + ζ_{(r p q)}) d z'),

(8)

where l is the length of a block that equals to the minimum domain width,

⌊ z / l ⌋

is a floor function which returns the greatest integer less than or equal to

z / l

, and

a_{(r p q)}

and

ζ_{(r p q)}

are generation parameters which determine the amplitude and phase of the spectrum of the final poling sequence. The sum of all different frequency terms is integrated over the block length from

⌊ z / l - 1 ⌋ l

⌊ z / l ⌋ l

. The corresponding sign function returns the polarization direction of the block.

By iteratively adjusting the values of

a_{(r p q)}

and

ζ_{(r p q)}

, it is possible to produce a solution for

d^{(2)} (z)

that will result in predetermined output electric field amplitudes and phases. In this algorithm the calculated amplitude

| {\tilde{d}}_{(r p q)} |

is mainly affected by

a_{(r p q)}

and the phase

φ_{{\tilde{d}}_{(r p q)}}

ζ_{(r p q)}

. This can be seen by recognizing in Eq. (8) that inclusion of

ζ_{(r p q)}

is equivalent to introducing an offset to the initial phase

Δ k_{(r p q)}

of the frequency components in the domain pattern. This feature makes the algorithm very robust when searching for the final solution since the amplitude and phase of the output electric field can be considered separately. Note that the choice of the form of Eq. (8) is somewhat arbitrary. Any solution that has a basic sinusoidal form will work, as different forms of sinusoids will only result in a shift or a sign change between

ζ_{(r p q)}

and the calculated value of

φ_{{\tilde{d}}_{(r p q)}}

. For convenience we have chosen the cosine function so that the calculated

φ_{{\tilde{d}}_{(r p q)}}

will have the same sign as the generation parameters of the phase

ζ_{(r p q)}

A numerical procedure based on simulated annealing algorism is used to generate the poling pattern of the crystal by first selecting a length for the domain block. The value of each

φ_{{\tilde{d}}_{(r p q)}}

can be estimated from the desired value of the final CEP using Eq. (6). The initial generation parameters

a_{(r p q)}

are then arbitrarily chosen to calculate a

d^{(2)} (z)

with Eq. (8). The

d^{(2)} (z)

that results from this calculation is fed into Eq. (2) to obtain

E_{n}

for comparison with the values of the fields of the desired waveform. If the root mean square of the difference of the calculated fields and the desired fields is larger than a prescribed value the procedure is repeated using a randomly chosen set of

a_{(r p q)}

until the calculated and the desired E fields agree to within an acceptable value.

4. Numerical examples

We demonstrate the utility of the above procedure with a few numerical examples. We use lithium niobate at 100∘C as the nonlinear medium. In the first example we aim to produce a sawtooth waveform at the crystal exit by the superposition of the first five harmonics with channels:

ω_{1} + ω_{1} \to ω_{2}

ω_{1} + ω_{2} \to ω_{3}

ω_{1} + ω_{3} \to ω_{4}

ω_{1} + ω_{4} \to ω_{5}

. In correspondence to our previous work [11

W. J. Chen, Z. M. Hsieh, S. W. Huang, H. Y. Su, C. J. Lai, T. T. Tang, C. H. Lin, C. K. Lee, R. P. Pan, C. L. Pan, and A. H. Kung, “Sub-single-cycle optical pulse train with constant carrier envelope phase,” Phys. Rev. Lett. 100(16), 163906 (2008). [CrossRef] [PubMed]

], we choose 2.406 μm as the fundamental wavelength. The phase mismatch of each channel,

Δ k_{(r p q)}

, is listed in Table 1 .

Table 1 Phase mismatch of selected processes^a

Channel	$Δ k$ (1/μm)	$2 π / Δ k$ (μm)
$ω_{1} + ω_{1} \to ω_{2}$	0.22380283	28.074647
$ω_{1} + ω_{2} \to ω_{3}$	0.31866532	19.717192
$ω_{1} + ω_{3} \to ω_{4}$	0.53402295	11.765759
$ω_{1} + ω_{4} \to ω_{5}$	0.89204589	7.0435673
$ω_{2} + ω_{3} \to ω_{5}$	1.2022660	5.2261190
$ω_{2} + ω_{5} \to ω_{7}$	3.7292894	1.6848211

^a In lithium niobate at 100

° C

and

λ_{1}^{(v a c u u m)} = 2.406

μm, calculated using refractive index data from Ref [15

D. H. Jundt, M. M. Fejer, and R. L. Byer, “Optical Properties of Lithium-Rich Lithium Niobate Fabricated by Vapor Transport Equilibration,” IEEE J. Quantum Electron. 26(1), 135–138 (1990). [CrossRef]

We choose to set the length of each block at

7.04 / 3 = 2.35

μm where 7.04 μm is the smallest QPM period of all channels used in this process (as shown in the first four rows in Table 1), i.e.,

2 π / (k_{5} - k_{4} - k_{1})

. The total number of blocks is 2475. Hence the crystal length is approximately 6 mm. We use the plane wave approximation and assume a rectangular temporal input pulse with an intensity of 100 MW/cm² and assume the initial phase

φ_{1} = 0

. After several iterations we obtain a set of generation parameters

a_{(211)} = 1.000

a_{(321)} = 0.996

a_{(431)} = 1.174

a_{(541)} = 1.294

and

ζ_{(211)} = ζ_{(321)} = ζ_{(431)} = ζ_{(541)} = 0

that results in a calculated field amplitude ratio of 1.000:0.501:0.340:0.252:0.194 for the 1st to 5th harmonics, which is close to the desired 1:1/2:1/3:1/4:1/5 for a square waveform. The absolute values of the corresponding effective second order nonlinear coefficients

| \tilde{d} (Δ k) |

is calculated to be 42%, 41%, 45% and 39% of a perfectly quasi-phase-matched periodic poled nonlinear crystal. The evolution of the E fields and the calculated spectrum of

| \tilde{d} (Δ k) |

of the aperiodic pattern generated with this set of parameters are shown in Fig. 1 .

Fig. 1 (left frame) Absolute values of the normalized effective second order nonlinear coefficients

| \tilde{d} (Δ k) |

calculated from Eq. (4) of the designed aperiodic optical superlattice and (right frame) simulation results of the electric field evolution in the designed AOS.

From Eqs. (2) and (8) we expect that by setting all

ζ_{(r p q)} = 0

, then all

φ_{{\tilde{d}}_{(r p q)}}

will be 0. The calculated values of

φ_{{\tilde{d}}_{(211)}}

φ_{{\tilde{d}}_{(321)}}

φ_{{\tilde{d}}_{(431)}}

, and

φ_{{\tilde{d}}_{(541)}}

are 0.0002π, 0.0003π, 0.0061π, −0.0022π, respectively, which match well with the expected values of 0 to within 0.006π. From Eq. (6) we get

φ_{C E} = 0.5 π - φ_{\tilde{d}} = 0.5 π

and

φ_{m} = φ_{1} - φ_{C E} = - 0.5 π

. The predicted phases of the harmonics are

φ_{2} = - 0.5 π

φ_{3} = - 1.0 π

φ_{4} = - 1.5 π

and

φ_{5} = - 2.0 π

. The actual simulation gives

φ_{1} = - 0.0007 π

φ_{2} = - 0.5019 π

φ_{3} = - 1.0030 π

φ_{4} = - 1.4983 π

and

φ_{5} = - 2.0123 π

, which also match well with the prediction.

This algorithm works for any value of the CEP. We used several sets of values of

ζ_{(r p q)}

ζ_{(r p q)}

equaling to 0, or

π / 2

, or π, or

3 π / 2

respectively and repeated the calculation. Note that each set of

ζ_{(r p q)}

corresponds to a different AOS design. The phase evolution inside the lithium niobate crystal for each set of

ζ_{(p q r)}

and the corresponding envelope and electric field intensity of the waveforms are shown in Fig. 2 . We see that the phases converge rather quickly to a value that agrees with the desired values in all cases.

Fig. 2 (large frames) Simulation results of the phase evolution inside the designed AOS and (small frames) the corresponding waveforms (solid curve) and envelopes (dotted curve). (a) For AOS generation parameters of phase

ζ_{(211)} = ζ_{(321)} = ζ_{(431)} = ζ_{(541)} = 0

, simulation results agree with expected values

φ_{C E} = 0.5 π

and

φ_{m} = - 0.5 π

.For

ζ_{(211)} = ζ_{(321)} = ζ_{(431)} = ζ_{(541)} = 0.5 π

, 1.0π and 1.5 π from (b) to (d), the expected values are (b)

φ_{C E} = 0

φ_{m} = 0

, (c)

φ_{C E} = - 0.5 π

φ_{m} = 0.5 π

, (d)

φ_{C E} = - 1.0 π

φ_{m} = 1.0 π

, respectively. In a stable periodic waveform,

φ_{C E}

affects the shape of the waveform, and

φ_{m}

produces a timeshift. This is confirmed with the results shown in the small frame figures on the right where the envelopes shift sequentially by 1/4 period from (a) to (d), which is due to the 0.5π phase shift of

φ_{m}

from one figure to the next.

In the event the externally imposed phase shift

ϕ_{n}

of the n th harmonics is not zero, we can compensate that by changing the values of

ζ_{(r p q)}

in Eq. (8). This will be the case, for instance, if we want to have a sawtooth waveform 1 m away instead of right at the crystal exit. Assuming plane wave propagation in air, then using published data on the refractive index of air [16

E. Peck and K. Reeder, “Dispersion of Air,” J. Opt. Soc. Am. 62(8), 958–962 (1972). [CrossRef]

], the phase evolution of the above harmonics are calculated to be −831484.366π, −1662970.062π, −24944584.580π, −33259510.110π, −41574493.192π, which correspond to

ϕ_{n} = - 0.366 π

, −0.062π, −0.580π, −0.110π, −1.192π, for

n = 1

, 2, 3, 4, 5, respectively. From Eq. (6) we can compensate this dispersion by shifting

ζ_{(211)}

-0
                     .366 π \times 2 - (- 0.062 π)

ζ_{(321)}

-0
                     .366 π - 0.062 π - (- 0.580 π)

and so on. Solving Eq. (8) with this set of

ζ_{(r p q)}

then results in a new AOS design that gives a sawtooth waveform located at 1 m away from the crystal, as shown in Fig. 3 . The waveform at the crystal end immediately after the generation process has a very different form compared to Fig. 2.

Fig. 3 Generated electromagnetic waveform (solid curve) and the corresponding envelope (dotted curve) which are (a) at the crystal exit and (b) at 1 m away from the crystal exit. The generation amplitude parameters are set to

a_{(211)} = 1.000

a_{(321)} = 0.996

a_{(431)} = 1.174

a_{(541)} = 1.294

, which are the same as those for calculating Fig. 2. The generation phase parameters are set to

ζ_{(211)} = 1.332 π

ζ_{(321)} = 0.030 π

ζ_{(431)} = 0.186 π

ζ_{(541)} = 1.942 π

to compensate the phase distortion due to dispersion in air.

The algorithm we have described indeed is very versatile. For many common waveforms such as the square wave or the triangular waveform, only a subset of the harmonics is needed to synthesize the desired waveform. A suitable generation path can be chosen with just the necessary harmonics to optimize the field power usage and remove the need to filter out unnecessary harmonics. To demonstrate the point, we designed an AOS pattern for generating a square wave. The Fourier components of a square wave have only odd frequency terms. For our example we employ a total of five harmonics

ω_{n}

n = 1

, 2, 3, 5, 7 with the channels

ω_{1} + ω_{1} \to ω_{2}

ω_{1} + ω_{2} \to ω_{3}

ω_{2} + ω_{3} \to ω_{5}

ω_{2} + ω_{5} \to ω_{7}

ω_{2}

is required in order to generate the higher harmonics and is expected to exhaust its energy before propagating into free space. The phase mismatch

Δ k_{(r p q)}

of each channel is listed in Table 1. The simulation uses a lithium niobate crystal that has 104000 blocks. The block length is 0.6 μm. The input fundamental intensity is 100 MW/cm². After several iterations we obtain the generation amplitude parameters to be

a_{(211)} = 1.000

a_{(321)} = 2.188

a_{(532)} = 12.119

a_{(752)} = 21.416

. Simulation shows that the normalized amplitudes of the 1st, 2nd, 3rd, 5th and 7th harmonics at the output are 1.000, 0.002, 0.343, 0.204, 0.140, respectively. The generated 2nd harmonic amplitude is less than 0.2% of the fundamental, showing that it has been effectively converted to the higher odd harmonics for the square wave. The waveforms corresponding to a few different CEPs are plotted in Fig. 4 . Note that only

φ_{C E} = 0.5 π

gives a square waveform.

Fig. 4 Plots of simulated waveforms (solid curve) and envelopes (dotted curve). The generation amplitude parameters are set to

a_{(211)} = 1.000

a_{(321)} = 2.188

a_{(532)} = 12.119

a_{(752)} = 21.416

and the generation phase parameters are set to

ζ_{(211)} = ζ_{(321)} = ζ_{(431)} = ζ_{(541)} = 0

, 0.25π, 0.5π, 0.75π from (a) to (d) respectively. From Eq. (8) and with

φ_{1} = 0

, it can be predicted that

φ_{C E}

will change from 0.5π to −0.25π. The expected values of the CEP are (a)

φ_{C E} = 0.5 π

φ_{m} = - 0.5 π

. (b)

φ_{C E} = 0.25 π

φ_{m} = - 0.25 π

. (c)

φ_{C E} = 0

φ_{m} = 0

. (d)

φ_{C E} = - 0.25 π

φ_{m} = 0.25 π

, which agree with above waveforms.

5. Discussion

Many physical processes are sensitive to the form of the electric field. Meanwhile the CEP affects the time evolution of the electric field waveforms. These waveforms will be time-invariant if the CEP is constant. In general, electric fields experience phase shifts and amplitude modulations that distort the waveform. If the shifts and modulations are predictable, such as if they originate from propagation through a dispersive medium, Gouy phase shift, reflection from optics surfaces, etc. then they can be pre-compensated in the design phase of the AOS domain pattern by including the shifts and modulations in the iterative process. In this way we can deliver the desired waveform to any predetermined location.

An advantage of employing the AOS is that all the necessary steps in waveform synthesis, including broadband frequency generation, and amplitude and phase modulation are all completed in a monolithic crystal so that the resulting waveform will be very stable and reproducible. One can conceive many important applications for such highly stable and reproducible waveforms. For example, measurements that are highly repetitious, correlation measurements that require an ultrafast standard or probe, and heterodyne measurements that use an ultrafast oscillator would all benefit from having stable waveforms generated from such an AOS. On the other hand, the AOS approach has the severe restriction that only one particular waveform can be obtained from each AOS design and for one input intensity. As such, stable, CW or quasi-CW high power lasers are preferred to provide the fundamental frequency input. Precise knowledge of the refractive index of the crystal and a tight temperature control must be available for an accurate waveform. The maximum number of harmonics which can be generated simultaneously is limited by the transmission bandwidth of the crystal and the fact that the poling of crystal domain widths of less than 1 μm over any useful length is unrealistic with today’s poling technology. This last restriction can be somewhat mitigated by introducing high order QPM or using a longer fundamental wavelength.

In summary we have described the synthesis of periodic waveforms using a monolithic AOS. An algorithm to design the AOS is proposed and verified by numerical simulation. The AOS will enable us to create a desired electric waveform and deliver it to a predetermined location. The approach offers a simple and compact solution to consistently and reliably producing such waveforms for many interesting applications.

Acknowledgments

We thank Chao-Kuei Lee for helpful discussions. This work was supported by the Academia Sinica and the National Science Council of the R.O.C. under the Nanoscience and Technology Program.

References and links

1.	Z. Jiang, D. S. Seo, S.-D. Yang, D. E. Leaird, R. V. Roussev, C. Langrock, M. M. Fejer, and A. M. Weiner, “Four-User, 2.5-Gb/s, Spectrally Coded OCDMA System Demonstration Using Low-Power Nonlinear Processing,” J. Lightwave Technol. 23(1), 143–158 (2005). [CrossRef]
2.	M. C. Stowe, F. C. Cruz, A. Marian, and J. Ye, “High resolution atomic coherent control via spectral phase manipulation of an optical frequency comb,” Phys. Rev. Lett. 96(15), 153001 (2006). [CrossRef] [PubMed]
3.	A. Baltuska, Th. Udem, M. Uiberacker, M. Hentschel, E. Goulielmakis, Ch. Gohle, R. Holzwarth, V. S. Yakovlev, A. Scrinzi, T. W. Hänsch, and F. Krausz, “Attosecond control of electronic processes by intense light fields,” Nature 421(6923), 611–615 (2003). [CrossRef] [PubMed]
4.	A. M. Weiner, J. P. Heritage, and E. M. Kirschner, “High-resolution femtosecond pulse shaping,” J. Opt. Soc. Am. B 5(8), 1563–1572 (1988). [CrossRef]
5.	J. K. Ranka, R. S. Windeler, and A. J. Stentz, “Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm,” Opt. Lett. 25(1), 25–27 (2000). [CrossRef]
6.	P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71(13), 1994–1997 (1993). [CrossRef] [PubMed]
7.	S. E. Harris and A. V. Sokolov, “Subfemtosecond Pulse Generation by Molecular Modulation,” Phys. Rev. Lett. 81(14), 2894–2897 (1998). [CrossRef]
8.	T. W. Hänsch, “A proposed sub-femtosecond pulse synthesizer using separate phase-locked laser oscillators,” Opt. Commun. 80(1), 71–75 (1990). [CrossRef]
9.	I. V. Shutov and A. S. Chirkin, “Consecutive high-order harmonic generation and formation of subfemtosecond light pulses in aperiodical nonlinear photonic crystals,” Phys. Rev. A 78(1), 013827 (2008). [CrossRef]
10.	M. Y. Shverdin, D. R. Walker, D. D. Yavuz, G. Y. Yin, and S. E. Harris, “Generation of a single-cycle optical pulse,” Phys. Rev. Lett. 94(3), 033904 (2005). [CrossRef] [PubMed]
11.	W. J. Chen, Z. M. Hsieh, S. W. Huang, H. Y. Su, C. J. Lai, T. T. Tang, C. H. Lin, C. K. Lee, R. P. Pan, C. L. Pan, and A. H. Kung, “Sub-single-cycle optical pulse train with constant carrier envelope phase,” Phys. Rev. Lett. 100(16), 163906 (2008). [CrossRef] [PubMed]
12.	Z. Jiang, C.-B. Huang, D. E. Leaird, and A. M. Weiner, “Optical arbitrary waveform processing of more than 100 spectral comb lines,” Nat. Photonics 1(8), 463–467 (2007). [CrossRef]
13.	N. Bloembergen, Nonlinear Optics (World Scientific, 1996).
14.	S. N. Zhu, Y. Y. Zhu, and N. B. Ming, “Quasi–Phase-Matched Third-Harmonic Generation in a Quasi-Periodic Optical Superlattice,” Science 278(5339), 843–846 (1997). [CrossRef]
15.	D. H. Jundt, M. M. Fejer, and R. L. Byer, “Optical Properties of Lithium-Rich Lithium Niobate Fabricated by Vapor Transport Equilibration,” IEEE J. Quantum Electron. 26(1), 135–138 (1990). [CrossRef]
16.	E. Peck and K. Reeder, “Dispersion of Air,” J. Opt. Soc. Am. 62(8), 958–962 (1972). [CrossRef]

OCIS Codes
(190.4160) Nonlinear optics : Multiharmonic generation
(320.5540) Ultrafast optics : Pulse shaping
(230.5298) Optical devices : Photonic crystals

ToC Category:
Photonic Crystals

History
Original Manuscript: June 19, 2009
Revised Manuscript: July 28, 2009
Manuscript Accepted: July 28, 2009
Published: August 28, 2009

Citation
Wei-Hsun Lin and A. H. Kung, "Arbitrary waveform synthesis by multiple harmonics generation and phasing in aperiodic optical superlattices," Opt. Express 17, 16342-16351 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-18-16342

Sort: Author | Year | Journal | Reset

References

Z. Jiang, D. S. Seo, S.-D. Yang, D. E. Leaird, R. V. Roussev, C. Langrock, M. M. Fejer, and A. M. Weiner, “Four-User, 2.5-Gb/s, Spectrally Coded OCDMA System Demonstration Using Low-Power Nonlinear Processing,” J. Lightwave Technol. 23(1), 143–158 (2005). [CrossRef]
M. C. Stowe, F. C. Cruz, A. Marian, and J. Ye, “High resolution atomic coherent control via spectral phase manipulation of an optical frequency comb,” Phys. Rev. Lett. 96(15), 153001 (2006). [CrossRef] [PubMed]
A. Baltuska, Th. Udem, M. Uiberacker, M. Hentschel, E. Goulielmakis, Ch. Gohle, R. Holzwarth, V. S. Yakovlev, A. Scrinzi, T. W. Hänsch, and F. Krausz, “Attosecond control of electronic processes by intense light fields,” Nature 421(6923), 611–615 (2003). [CrossRef] [PubMed]
A. M. Weiner, J. P. Heritage, and E. M. Kirschner, “High-resolution femtosecond pulse shaping,” J. Opt. Soc. Am. B 5(8), 1563–1572 (1988). [CrossRef]
J. K. Ranka, R. S. Windeler, and A. J. Stentz, “Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm,” Opt. Lett. 25(1), 25–27 (2000). [CrossRef]
P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71(13), 1994–1997 (1993). [CrossRef] [PubMed]
S. E. Harris and A. V. Sokolov, “Subfemtosecond Pulse Generation by Molecular Modulation,” Phys. Rev. Lett. 81(14), 2894–2897 (1998). [CrossRef]
T. W. Hänsch, “A proposed sub-femtosecond pulse synthesizer using separate phase-locked laser oscillators,” Opt. Commun. 80(1), 71–75 (1990). [CrossRef]
I. V. Shutov and A. S. Chirkin, “Consecutive high-order harmonic generation and formation of subfemtosecond light pulses in aperiodical nonlinear photonic crystals,” Phys. Rev. A 78(1), 013827 (2008). [CrossRef]
M. Y. Shverdin, D. R. Walker, D. D. Yavuz, G. Y. Yin, and S. E. Harris, “Generation of a single-cycle optical pulse,” Phys. Rev. Lett. 94(3), 033904 (2005). [CrossRef] [PubMed]
W. J. Chen, Z. M. Hsieh, S. W. Huang, H. Y. Su, C. J. Lai, T. T. Tang, C. H. Lin, C. K. Lee, R. P. Pan, C. L. Pan, and A. H. Kung, “Sub-single-cycle optical pulse train with constant carrier envelope phase,” Phys. Rev. Lett. 100(16), 163906 (2008). [CrossRef] [PubMed]
Z. Jiang, C.-B. Huang, D. E. Leaird, and A. M. Weiner, “Optical arbitrary waveform processing of more than 100 spectral comb lines,” Nat. Photonics 1(8), 463–467 (2007). [CrossRef]
N. Bloembergen, Nonlinear Optics (World Scientific, 1996).
S. N. Zhu, Y. Y. Zhu, and N. B. Ming, “Quasi–Phase-Matched Third-Harmonic Generation in a Quasi-Periodic Optical Superlattice,” Science 278(5339), 843–846 (1997). [CrossRef]
D. H. Jundt, M. M. Fejer, and R. L. Byer, “Optical Properties of Lithium-Rich Lithium Niobate Fabricated by Vapor Transport Equilibration,” IEEE J. Quantum Electron. 26(1), 135–138 (1990). [CrossRef]
E. Peck and K. Reeder, “Dispersion of Air,” J. Opt. Soc. Am. 62(8), 958–962 (1972). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper