Enhanced Second Harmonic Generation Efficiency via Wavefront

Jun 6, 2017 - state) of second harmonic light can be controlled by changing the phase profile of .... with intensities much brighter than those withou...
0 downloads 0 Views 6MB Size
Subscriber access provided by TRENT UNIV

Article

Enhanced second harmonic generation efficiency via wavefront shaping Jonathan V. Thompson, Brett H. Hokr, Graham A. Throckmorton, Dawei Wang, Marlan O Scully, and Vladislav V. Yakovlev ACS Photonics, Just Accepted Manuscript • Publication Date (Web): 06 Jun 2017 Downloaded from http://pubs.acs.org on June 9, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Photonics is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Enhanced second harmonic generation eciency via wavefront shaping Jonathan V. Thompson, Wang,





Brett H. Hokr,

Marlan O. Scully,

†, ‡



Graham A. Throckmorton,

and Vladislav V. Yakovlev



Dawei

∗, †

Texas A&M University, College Station, TX 77843, and Baylor University, Waco, TX 76798 E-mail: [email protected]

Abstract

Nonlinear Crystal

Optical second harmonic generation is a fundamental nonlinear eect with a large impact on laser technology and optical imaging/sensing. For most practical applications of second harmonic generation, high conversion eciency is required.

Nonlinear Crystal

However, many techniques used to

achieve high eciency are limited to fabrication methods and optical energy requirements. Here, we investigate and demonstrate substantial enhancement in the conversion eciency of second harmonic generation in a nonlinear crystal via the application of wavefront shaping. In

Figure 1: Control of the second harmonic gen-

a one-dimensional understanding of second har-

eration conversion eciency is achieved by ap-

monic generation, a phase oset applied to the

plication of a phase mask to the incident fun-

fundamental wave has no eect on the intensity

damental wave.

of the generated light. We show that when the

Keywords

pump eld is not a plane wave, enhanced conversion eciency can be controlled by the application of a phase mask to the fundamental

wavefront shaping,

beam. This investigation of the dependence of

second harmonic genera-

tion, optimization algorithm, nonlinear crystal,

conversion eciency upon the transverse phase

conversion eciency, nonlinear generation

prole of the incident pump laser yields the promise of new ways to enhance nonlinear gen-

Optical second harmonic generation (SHG)

eration with limited optical energy and without

was the rst discovered nonlinear optical ef-

the need for specic fabrication techniques.

fect,

1

and is now widely used for wavelength



conversion in green



sensing and chemical structure assessment,

To whom correspondence should be addressed 1 ‡ 2

laser pointers,

and biomedical imaging.

79

optical

26

For example, sec-

ond harmonic generation imaging microscopy has been used to image molecular structures such as collagen

ACS Paragon Plus Environment 1

10

7

in tissues. This is useful for

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the detection and research of cancer, and connective tissue disorders.

11

10,12

brosis,

Page 2 of 12

ing of the wavefront does have an eect on the nonlinear process. For example, the transverse

The conversion eciency of a nonlinear pro-

mode (beam state) of second harmonic light can

cess is a critical parameter for practical ap-

be controlled by changing the phase prole of

plications.

the incident fundamental beam.

Over the years, multiple methods

40

have been used to improve the eciency of

In this article, we present theoretical and ex-

second harmonic generation processes including

perimental results which indicate that shaping

the selection of the nonlinear material,

in-

of the phase prole (i.e. wavefront shaping) of

the use of

the fundamental beam not only controls the sec-

teraction length of the crystal, higher intensities, riodic poling,

18

16

13,15,16

photonic crystals,

13,14

17

or pe-

ond harmonic beam shape, but also substan-

to name a few. When it comes

tially impacts the intensity or conversion e-

to nonlinear biomedical imaging and sensing,

ciency of second harmonic generation, as con-

the above methods are not always applicable;

ceptually illustrated in Fig. 1.

however, the strength of the signal generated

demonstrated using second harmonic conver-

through a nonlinear optical interaction will de-

sion eciency from a nonlinear crystal (with

termine the ultimate sensitivity and imaging

no scattering present) as a feedback param-

speed and depth ultimate sensitivity, imaging

eter for wavefront shaping optimization algo-

speed, and depth .

rithms. Recent advances in wavefront shaping

19

It is therefore imperative

4143

This eect is

44,45

to explore new ways to increase the nonlinear

techniques

conversion eciency. These advances will likely

vided ecient and accurate tools to modify the

come from our through better understanding

phase prole of a laser beam.

of the nonlinear optical process in all three di-

clude search algorithms for feedback-based op-

mensions, and, more particularly, from better

timization

understanding the role of the transverse phase

utilization of these search algorithms with sec-

prole on of the incident fundamental beam.

ond harmonic generation as feedback produces

A one dimensional theory

20,21

41,46

and technologies

have pro-

These tools in-

of the applied phase mask. The

is often su-

a phase mask for which the generation eciency

cient to describe the generation of second har-

is enhanced. These techniques have been suc-

monic (SH) light.

Three-dimensional treat-

cessfully demonstrated in highly scattering en-

ment, however, is necessary when the incident

vironments with both second harmonic genera-

beam is either naturally disturbed by dirac-

tion

tion

back.

22,23

and scattering,

24

or through active

47

and two-photon uorescence

48

as feed-

Here, we show enhancements in sec-

preparation of particular beam states rang-

ond harmonic conversion eciency as high as

ing

and

an order of magnitude when no scattering is

Other

present by implementing similar techniques and

from

Gaussian

other sophisticated

22,25,26

3033

to

Bessel

2729

beams states.

instances that require three-dimensional treat-

iterative search algorithms.

ment include the treatment of articial gauge

show enhancements in second harmonic conver-

elds in photon uids,

SHG from photonic

sion eciency as high as an order of magnitude

intensity-dependent phase shifts due

by implementing an iterative search algorithm.

crystals,

35

34

to pump elds that are not plane-waves,

22,36

In particular, we

as

These results may impact laser technologies in

well as shaping of the second harmonic beam

the form of lower input powers required for fre-

through facet engraving,

quency doubling, as well as increasing the sen-

phy,

38,39

37

nonlinear hologra-

and plasmonic metasurfaces.

33

sitivity of second harmonic biomedical imaging

Of particular interest here is the treatment

and other detection/sensing techniques. These

of nonlinear generation when the wavefront, or

results should also translate to many other non-

phase prole, of the pump beam is non-uniform.

linear optical conversion methods.

In the one dimensional case (plane wave), a phase shift in the fundamental beam has no impact upon nonlinear generation.

However,

when the pump beam is not a plane wave, shap-

ACS Paragon Plus Environment 2

Page 3 of 12

Photodiodes

6ps Laser

SHGF

Lens Beamsplitter

CCD

Half-wave plate

Lens

Polarizing beam cube

Lens

I0

Glass Slide

Mirror

LCOS SLM

SHG Crytal

SHGB

Mirror

a)

b)

Figure 2: Experimental setup. (a) The spatial phase of a 6 picosecond laser pulse was modied via spatial light modulator (SLM). The phase mask was imaged onto the sample with a 4-f imaging setup. For enhanced generation from a nonlinear crystal (b), the forward and backward generated signals (SHG F , SHGB ) as well as a reference signal (I 0 ) were collected with photodiodes. These signals were used to create metrics that drove an optimization algorithm to enhance the second harmonic generation by wavefront shaping.

Numerical Simulations Normalized Conversion Efficiency (SHG/I0)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

where



10

represent the fundamental and

= 2ω1 ), respectively, nj is the refractive index for the j th eld, o is the permittivity of free space, and c is the

3π/2

8

j = 1, 2

second harmonic elds ( ω2

speed of light. The electric eld and nonlinear polarization are given by

z = 5 mm

6

π

Phase Mask

E˜j (r, t) = A˜j (r)ei(kj z−ωj t) + c.c. P˜j (r, t) = P˜j (r)e−iωj t + c.c.

4 π/2 2

z = 5 mm z =10 mm 0

200 400 600 800

The tilde indicates complex variables,

z =10 mm

0

Phase Mask

monic generation,

P˜j

is written as

is the

21

P˜1 (r) = 4o deff A˜2 A˜∗1 ei(k2 −k1 )z P˜2 (r) = 2o deff A˜2 e2ik1 z

Simulated enhancement of second

harmonic generation via wavefront shaping for dierent interaction lengths (left).



complex eld amplitude, and for second har-

Algorithm Iteration Figure 4:

(2)

(3)

1

The opti-

14

mized phase masks are shown on the right.

Here,

Results

of the medium. Solution of the wave equation

is the eective nonlinear coecient

(Eq. 1) yields the intensity of the two frequencies of light via the relation

Theoretical Model

2 Ij = 2nj 0 c A˜j (z) .

Second harmonic generation is modeled via the three-dimensional wave equation

∇2 E˜j −

deff

(4)

21,22,49

Substitution of Eq. 2 into Eq. 1, and applica-

tion of the slowly n2j ∂ 2 E˜j 1 ∂ 2 P˜j = (1) 2 Paragon Plus Environment ACS c2 ∂t2 o c2 ∂t 3

varying approximation yields

ACS Photonics

2.2

Enhancement in Backward Direction 2π

a)

Normalized Conversion Efficiency (SHGB/I0)

Enhancement in Forward Direction Normalized Conversion Efficiency (SHGF/I0)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 12

11



b)

10

2 3π/2 1.8 Phase Mask: Large

1.6

π

1.4 π/2 1.2 1

Large Pixels Small Pixels 0

0

200 400 600 800 Phase Mask: Small

9

3π/2

8 7 Phase Mask: Large

6

π

5 4 π/2

3 Large Pixels Small Pixels

2 1

Algorithm Iteration

0

200 400 600 800 Phase Mask: Small

0

Algorithm Iteration

Figure 3: Enhanced conversion eciency of second harmonic generation in in the forward (a) and backward (b) direction from a nonlinear crystal due to wavefront shaping of the fundamental. For both cases, the nal phase masks as determined by the search algorithm are shown for large (top) and small (bottom) pixel sizes.

a dierential equation for the complex amplitude,

where

22

F1 =

∂ A˜j (r) i 2 ˜ i(ωj )2 ˜ = Pj (r)e−ikj z ∇⊥ Aj (r) + ∂z 2kj 2kj o c2

F2 =

ω22 deff A21 . k 2 c2

(9)

the amplitude has a dependence upon the trans-

The dependence of the eld amplitude on the

verse gradient of the phase. As opposed to the

transverse phase is more apparent by explicitly

one-dimensional case (plane wave), these terms

writing the amplitude and phase as

A˜j (r) = Aj (r)e

2ω12 deff A2 A1 , k 1 c2

and

From the equations in this form, we see that

(5)

iφj (r)

∆k = k2 − 2k1 , ∆φ = φ2 − 2φ1 ,

will survive the absolute value in Eq. 4. Thus, a phase mask applied to the fundamental beam

(6)

will inuence the intensity or conversion eciency of the generated light.

and separating the real and imaginary parts of Eq. 5. This yields the following system of equa-

Experimental Results

tions.

  −1 2 (∂x Aj )(∂x φj ) ∂z Aj (r) = 2kj   2 + (∂y Aj )(∂y φj ) + Aj ∇⊥ φj  + (−1)j Fj sin ∆kz + ∆φ    −1 2 2 Aj ∂z φj (r) = Aj (∂x φj ) + (∂y φj ) 2kj  2 − ∇⊥ Aj  + Fj cos ∆kz + ∆φ

The

enhancement

of

second

harmonic

con-

version eciency was experimentally demonstrated (see Methods for details) by applying a

(7)

phase mask to a 6 picosecond laser pulse with a wavelength of 1064 nm (Attodyne, APLX-1064532). As depicted in Fig. 2, the phase mask was applied to the fundamental pulse with a spatial light modulator (LCOS-SLM, Hamamatsu X10468-08).

(8)

The applied phase mask was

then imaged with a 4-f telescope conguration (magnication of 1/8) to a Potassium Titanyl Phosphate (KTiOPO 4 , KTP) nonlinear crystal

ACS Paragon Plus Environment 4

Page 5 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

(type II, rotated 45



with respect to the beam

agation length, and from 2.3% to 23% for a

polarization). The second harmonic light gener-

10 mm propagation length.

ated in both the forward and backward (reec-

eciency from the back side of the crystal is

tion from back surface) directions was measured

taken into account, these values correspond to

with respect to the incident fundamental using

an enhancement in conversion eciency in the

large area photodiodes (Thorlabs, DET100A).

backward direction from 0.01% to 0.1%.

The conversion eciency was used as feedback

optimum phase masks calculated by the simu-

for a search algorithm, which determined an op-

lation are also shown in Fig. 4. Here, the phase

timum phase mask for enhancement.

pattern more clearly shows a checkerboard-like

In this

case, we used both a genetic algorithm a continuous sequential algorithm

41

46

and

behavior.

If the reection

The

While this pattern is not as clearly

in an out-

visible in the experimental results, possibly due

Both algorithms

to experimental factors such as noise and slight

produced similar results; therefore, only results

aberration eects, large phase jumps are clearly

from the genetic algorithm are shown here. Be-

observed.

cause enhancement may also depend upon the

the phase, gives sharp jumps in the transverse

resolution of the phase mask applied, we also

gradient of the phase which in turn eects the

show results for phase masks with low resolu-

conversion eciency (see Eqs. 7, 8).

ward spiral conguration.

42

This behavior, with large jumps in

tion (large macropixels) and higher resolution

The transverse intensity prole before and af-

(small macropixels). At the plane of the crystal,

ter wavefront shaping of the fundamental and

these macropixels had sizes of (150 ×180)

second harmonic pulses are shown in Fig. 5.

and (60×60)

µm,

µm

respectively.

Here, a grid-like pattern in both the fundamen-

The second harmonic conversion eciency in

tal and second harmonic proles is caused by

both the forward and backward directions as

the wavefront shaping.

a function of algorithm iteration are shown in

positions in the prole with intensities much

Fig. 3.

brighter than those without shaping.

The second harmonic conversion e-

This causes localized Numer-

ciency, calculated from the second harmonic di-

ical simulation results showing similar behavior

vided by the incident fundamental ( ISHG /I0 ),

are shown in Fig. 6a-b. Here, a grid-like pattern

is shown in Fig. 3 as a function of algorithm

caused by the shaped wavefront is also observed

iteration for both the forward and backward

in the fundamental and generated light at the

directions. For comparison, the conversion ef-

output of the interaction region.

ciencies were normalized by the initial (at

like patterns are likely due to the phase jumps

wavefront) value, which was 9.7% for the for-

in the transverse direction of the applied phase

ward generated light, and 0.02% for the back-

mask.

ward generated light. After 800 iterations, the

the enhancement of second harmonic conversion

conversion eciency was enhanced by a factor

is sensitive to the transverse phase gradients.

of two (19.5%) in the forward direction and

Therefore, we expect to see enhancement local-

as much as an order of magnitude (0.2%) in

ized to the regions where the phase gradient is

the backward direction. The nal phase masks

non-zero.

(both large and small macropixels) as deter-

Examination of Eqs. 7 & 8 shows that

Discussion

mined by the algorithm are also shown in the gures. Similar results from numerical calculations are shown in Fig. 4. As depicted in the g-

Further analysis of the theoretical model pro-

ure, numerical simulations were performed for

vides more insight into the eects of wave-

dierent crystal (propagation) thicknesses. Dif-

front shaping upon second harmonic genera-

ferent enhancements were achieved based upon

tion. Taking a cross-section of the second har-

the crystal lengths, but these enhancements

monic intensity as is propagates through the

are not necessarily linear with propagation distance.

These grid-

crystal (Fig. 6c-e) shows that after long propa-

Here, the conversion eciency was en-

gation distances, the energy in the second har-

hanced from 4.4% to 15.2% for a 5 mm prop-

ACS Paragon Plus Environment 5

ACS Photonics Fundamental Profile

Numerical Simulation Profile

1

1

a)

0.35

b)

0.3

0.8

0.8 0.25

0.6 Flat wavefront, Large pixels

0.6

Shaped wavefront, Large pixels

Fundamental

0.2

Second Harmonic

0.4

0.15

0.4

0.2

0.1 0.2

Flat wavefront, Small pixels

Shaped wavefront, Small pixels

0.05

0 Shaped Fundamental

Second Harmonic Profile

Shaped Second Harmonic

0

1

0

0.8 0.6 Shaped wavefront, Large pixels

0.4

c)

1 0.5 0 -0.5 -1 -1.5

Flat wavefront

0

0.2 Flat wavefront, Small pixels

Shaped wavefront, Small pixels

2.5

5

7.5

1.5

d)

1

1x106

0 -0.5

5x105

-1 -1.5

Shaped wavefront

0

Propagation (mm)

0

Figure 5: Transverse prole of the fundamental (top) and second harmonic (bottom) pulse af-

2.4 2 1.6 1.2 0.8 0.4 0

Flat Shaped

2.5

5

7.5

1.5x106

5

0

7.5

f)

1.0x106 5.0x105

10

Propagation (mm)

ter the nonlinear crystal with (right side) and

2.5

Propagation (mm)

e)

0

2x106

0.5

Intensity (W/cm2)

Flat wavefront, Large pixels

1.5

Transverse Position (mm)

Transverse Position (mm)

Propagation of Second Harmonic

SHG Power (kW)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 12

0

-1

0

1

Transverse Position (mm)

without (left side) wavefront shaping applied. Both large and small pixel sizes where used for

Figure 6:

wavefront shaping.

fundamental (a) and second harmonic (b) pulse

Transverse prole of the simulated

after the nonlinear crystal with and without wavefront shaping applied.

monic eld pumps back into the fundamental eld.

the simulated second harmonic intensity inside

At this point, it is important to recall

the interaction region is also shown for at (c)

that the simulation assumes the wave-vector mismatch is zero ( k2

= 2k1 ).

and shaped (d) incident wavefronts.

This behavior

was also observed by Sheng et al

22

tion gives the power as a function of propagation distance (e). A slice of the intensity at the

tween second harmonic generation and prop-

output of the crystal is also shown (f ).

agation eects when the pump beam is not a

36

Integra-

tion of the intensity over the transverse direc-

and is due

to diraction and the complex interplay be-

plane-wave.

A cross section of

When wavefront shaping is ap-

plied, this parametric pumping of energy back

the turnover in harmonic power occurs with re-

into the fundamental eld is counteracted, and

spect to the back surface of the crystal.

the harmonic eld continues to grow. Integra-

length dependence explains the dierence in en-

tion over the transverse direction shows this ef-

hancement between the forward and backward

fect in the average power as a function of prop-

direction experimental results. Because the sec-

agation distance (Fig. 6e). This indicates that

ond harmonic light detected in the backward

the enhancement attainable via wavefront shap-

direction was likely due to a reection from the

ing is dependent upon the length of the crys-

back surface of the crystal, this light propa-

tal.

gated through an eectively thicker nonlinear

This dependence is not necessarily linear

with crystal length, but it depends upon where

ACS Paragon Plus Environment 6

This

Page 7 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

medium. This indicates that while the enhance-

monic generation in nonlinear crystals. These

ment attainable via wavefront shaping is depen-

results should apply to other nonlinear eects,

dent upon the length of the crystal it is not nec-

and may also be related to articial gauge elds

essarily directly or linearly proportional to the

in photon uids.

34

crystal length. More specically, Fig. 6e shows

Because there is little to no scattering in-

that the enhancement depends upon the loca-

volved in the problem, certain identiable pat-

tion of the turnover in harmonic power with re-

terns in the phase mask may be responsible for

spect to the exit aperture of the crystal. Thus,

the majority of the enhancement. This has been

the enhancement of generated signal is depen-

observed, with the genetic algorithm converg-

dent upon the crystal length; however, a factor

ing to a phase mask with patterns consisting

of two longer crystal will not necessarily yield

of large (close to

a factor of two larger enhancement. For exam-

ing of this phase pattern for specic applica-

ple, according to Fig. 6e, signal from a crystal

tions will remove the need for a feedback fed

that is 5 mm long will be enhanced by a factor

search algorithm. Furthermore, more sophisti-

of 3.4 times, whereas signal from a 10 mm long

cated feedback metrics may also provide for the

crystal will be enhanced by a factor of 10. This

extension of this technique to biological applica-

dependence explains the dierence observed in

tions including the enhancement of second har-

the enhancement between the experimental re-

monic images of collagen structures in a mostly

sults in the forward and backward direction.

transparent sample.

Because the second harmonic light detected in

phase jumps. Understand-

Methods

the backward direction was likely due to a reection from the back surface of the crystal, this light propagated through an eectively thicker

Experimental Setup

nonlinear medium. A slice through the brightest points of the in-

The experimental setup is illustrated in Fig. 2.

tensity prole at the output of the crystal is

Here, a phase mask was applied to a 1064 nm, 6

shown in Fig. 6f for both shaped and at input wavefronts.

π)

picosecond laser pulse (Attodyne APLX-1064-

This gives better perspective

532) by a liquid crystal on silicon spatial light

to see that not only are the bright intensity

modulator (LCOS-SLM, Hamamatsu X10468-

spikes enhanced, but the darker regions are at

08). Prior to phase shaping, the beam was ex-

least as bright or brighter than generation from

panded by a telescope to ll more of the active

the unshaped fundamental. It also shows how

region of the spatial light modulator.

drastic the intensity spikes may be. While this

A half-

wave plate and polarizing beam cube were used

spatial prole may not be ideal for all appli-

to control the power and set the polarization

cations, there are instances, including optical

to the phase-only modulation orientation of the

detection, when the spatial prole is less impor-

LCOS-SLM. The phase mask applied to the

tant than the attainment of greater conversion

beam was imaged onto the face of a nonlinear

eciencies.

crystal using lenses in a 4-f imaging congura-

Conclusion

tion. The magnication of the lenses (Thorlabs,

In conclusion, we have theoretically and experi-

in diameter (0.7 mm, FWHM) than that of the

mentally demonstrated a dependence of second

nonlinear crystal (5 mm).

AC254-400-C-ML, AC254-050-C-ML) was chosen to ensure that the shaped pulse was smaller

The second harmonic light was generated

harmonic conversion eciency upon the phase prole of the incident fundamental wave.

by a Potassium Titanyl Phosphate (KTiOPO 4 , ◦

We

have shown that the utilization of wavefront

KTP) nonlinear crystal (type II, rotated 45

shaping techniques can yield a substantial in-

with respect to the beam polarization).

crease to the conversion eciency of second har-

pulse energy incident upon the crystal was ap-

ACS Paragon Plus Environment 7

The

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

proximately 2.3

µJ (uence of 100 µJ/cm2 ) with

a repetition rate of 10 kHz.

Page 8 of 12

els were (150 ×180)

The crystal was

and (60×60)

µm

plane of the crystal.

µm

at the

For the large or small

×

phase matched by hand and the orientation was

macropixels, there were a total of (10

rotated about the optical axis to optimize the

(25

conversion eciency of the unshaped pulses.

respectively. The genetic algorithm was allowed

The second harmonic light generated (532 nm)

to select a value between 0 and 2 π as the phase

was detected in both the forward and backward

delay for each of the macropixels.

×

11) or

33) macropixels in the modulator array,

directions using large area photodiodes (Thor-

After wavefront optimization via the search

labs DET100A) with glass color lters (Thor-

algorithm, an image of the output from the non-

labs FGB39) to reject the fundamental.

The

linear crystal in the forward direction was ac-

light detected in the backward direction was

quired using imaging lenses (Thorlabs, AC254-

likely due to reections from the back surface

035-A-ML, AC254-045-A-ML) and a charge-

of the crystal, and was separated from the inci-

coupled device (CCD) (Mightex CGE-B013-U).

dent fundamental for detection by reection o

The lenses were placed such that the back sur-

a glass slide (1 mm thick). This same glass slide

face of the crystal was imaged to the CCD. A

was also used to direct a small proportion of

glass color lter (Thorlabs FGB39) was used to

the incident beam into an additional large area

reject the fundamental light when imaging the

photodiode (Thorlabs DET100A) to be used as

generated second harmonic (532 nm).

a reference (I0 ). Both the reference beam and

mirror (Thorlabs M254C00) was used to reject

the forward generated beam were further at-

the second harmonic light when imaging the

tenuated by neutral density lters with 2 OD

fundamental.

(Thorlabs ND20B, ND10B) placed directly be-

without disruption to the alignment. Further-

fore the photodiodes.

more, neutral density lters were used to avoid

A cold

These lters were interchanged

The signal from the photodiodes was collected

saturation of the CCD in both cases. Care was

and averaged 128 times by an oscilloscope (Ag-

also taken at each stage to prevent reections

ilent DSO6034A). This signal was used as feed-

from optics in the forward direction from prop-

back (I2 /I1 ) (ISHG /I0 ) for an optimization al-

agating back into the nonlinear crystal.

gorithm which searched for a favorable phase

Numerical Simulations

mask to enhance the feedback metric; conversion eciency in this case.

46

A genetic algo-

with a population of 20 individuals cy-

Comparison of the experimental results with

cled over 800 generations was used as the op-

second harmonic theory were performed via nu-

timization algorithm.

merical simulation of Eq. 5.

rithm

The optimization algo-

rithm used was a genetic algorithm

46

Here, a spectral

method in the transverse direction,

with a

×

22

with a

population of 20 individuals per generation, cy-

grid of 128

cled over 800 generations.

The results from

tion with a nite dierence step in the propa-

this algorithm were compared with those from

gation (z ) direction. Enhancements for a total

a continuous sequential algorithm

in an out-

propagation distance of both 5 mm and 10 mm

and found to be

were computed with 128 and 256 total steps

ward spiral conguration similar.

42

41

in the

Therefore, only results from the ge-

z

128 points, was used in conjunc-

direction, respectively. This grid size

netic algorithm have been presented here. Be-

was sucient to produce results comparable to

cause the spatial resolution of the spatial light

established theory for at wavefronts.

modulator does have inuence on the magniused by grouping the pixels of the spatial light

we used an initial fundamental beam waist of 6 2 0.7 mm and intensity of 2 × 10 W/cm (9.8 × 103 W total power). We assumed the eective

modulator into macropixels composed of either

nonlinear coecient to be

tude of the results, two dierent pixel sizes were

(60

× 72) or (24 × 24) pixels. With an individµm and a magnication of

the index of refraction,

22

Here

deff = 1.4 cm/V, and n = 1.8, for both the

ual pixel size of 20

fundamental (1064 nm) and second harmonic

1/8 from the 4-f imaging setup, these macropix-

(532 nm) waves.

ACS Paragon Plus Environment 8

14

It is also important to note

Page 9 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

References

that the simulation assumes that phase matching conditions are satised (i.e. mismatch is zero).

wave-vector

(1) Franken, P. A.; Hill, A. E.; Peters, C. W.;

The simulation receives a

Weinreich, G. Generation of optical har-

phase mask as input, and calculates the output intensity distribution. A genetic algorithm

46

monics.

Phys. Rev. Lett. 1961, 7, 118119.

comparable to that

(2) Shen, Y. R. Surface properties probed by

used in the experiment was used to drive the

second-harmonic and sum-frequency gen-

simulation to nd an optimum phase mask for

eration.

enhanced second harmonic generation.

Here,

Nature 1989, 337, 519525.

the generated intensity of second harmonic from

(3) Corn, R. Optical second harmonic gen-

a population of 20 fundamental beams with ran-

eration as a probe of surface chemistry.

Chem. Rev. 1994, 25, 4113.

domly generated phase masks were evaluated and ranked. Each pixel in the phase mask was given a value between

0

and

2π .

(4) Eisenthal, K. B. Liquid interfaces probed

An individual

by second-harmonic and sum-frequency

phase mask in the next generation of masks was

spectroscopy.

created by mixing two masks of the current generation with using a random binary template, (n) (n) T . Thus, φ(n+1) = T φ1 + (1 − T )φ2 . The two

(5) Stupp,

parent phase masks were randomly chosen with

(6) Williams, Webb,

nth generation were al-

2π .

Here, the mutation rate decay,

Supramolecular

materi-

W.

R.

M.;

W.

Zipfel,

Interpreting

Science

W.

R.;

second-

harmonic generation images of collagen I

lowed to mutate to a random number between and

I.

als: self-organized nanostructures.

with highest rank. Furthermore, the a mutation (n) rate of R = [(0.088)e(−n/λ) +0.012] percent of

0

S.

1997, 276, 384389.

a triangular weight giving priority to the masks

the pixels of a mask per

Chem. Rev. 1996, 96, 1343

1360.

brils.

λ,

ensures that fewer mutations occur after many

Biophys. J. 2005, 88, 137786.

(7) Campagnola, P. J.; Loew, L. M. Second-

generations, and was set to 100 for both the

harmonic imaging microscopy for visual-

simulations and the experiment.

izing biomolecular arrays in cells, tissues and organisms.

Acknowledgements

Nat. Biotechnol. 2003, 21,

13561360. (8) Zoumi, A.;

We acknowledge the support of the National

vivo by using second-harmonic genera-

#1532188 and ECCS award #1509268), the

tion and two-photon excited uorescence.

Proc. Natl. Acad. Sci. U. S. A. 2002, 99,

A-

1261), the Oce of Naval Research (Grant No.

110149.

N00014-16-1-3054), and the US Department of Defense (grant #FA9550-15-1-0517).

Tromberg, B. J.

Imaging cells and extracellular matrix in

Science Foundation (DBI awards #1455671, Robert A. Welch Foundation (Grant No.

Yeh, A.;

(9) Brown, E.;

JVT is

McKee, T.;

DiTomaso, E.;

supported by the Herman F. Heep and Min-

Pluen, A.; Seed, B.; Boucher, Y.; Jain, R.

nie Belle Heep Texas A&M University Endowed

Dynamic imaging of collagen and its mod-

Fund held/administered by the Texas A&M

ulation in tumors in vivo using second-

Foundation. BHH acknowledges a graduate fel-

harmonic generation.

lowship from the Department of Defense Sci-

796800.

ence, Mathematics and Research for Transfor-

Nat. Med. 2003, 9,

(10) Chen, X.; Nadiarynkh, O.; Plotnikov, S.;

mation (SMART) fellowship program.

Campagnola, P. J. Second harmonic generation microscopy for quantitative analysis of collagen brillar structure.

toc. 2012, 7, 654669.

ACS Paragon Plus Environment 9

Nat. Pro-

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(11) Tilbury, K.; Campagnola, P. J. Applica-

2015,

7,

(22) Sheng, S. C.; Siegman, A. E. Nonlinearoptical calculations using fast-transform

2132.

methods:

(12) Mao, J. R.; Bristow, J. The Ehlers-Danlos syndrome: on beyond collagens.

Invest. 2001, 107, 10639. (13) Driscoll,

T.

A.;

Perkins,

with depletion and diraction.

J. Clin.

E.;

generation

A 1980, 21, 599606.

Ho-

harmonic

Phys. Rev.

McIver, J. K. Second-

generation

in

a

nonlinear,

man, H. J.; Stone, R. E. Ecient second-

anisotropic medium with diraction and

harmonic generation in KTP crystals.

depletion.

Opt. Soc. Am. B (14) Dmitriev, gosyan, D.

(15)

Second-harmonic

(23) Dreger, M. A.; P.

3rd ed.;

Academic Press: Burlington, MA, 2008.

ing microscopy in ovarian and breast can-

Perspect. Medicin. Chem.

Nonlinear Optics,

(21) Boyd, R. W.

tions of second-harmonic generation imagcer.

Page 10 of 12

V.;

1986, 3, 683.

Gurzadyan,

J.

G.;

Niko-

(24) Frostig, H.;

Handbook of nonlinear optical

Silberberg,

crystals ; Springer-Verlag: Berlin, 1991. Shen, Y. The principles of nonlinear optics ; Wiley: New York, 1984.

(17) Chen,

M.-L.;

Liu,

beams.

de

R.-J.;

crystals.

Light Sci.

A.

Sterke,

gett,

1968,

C.

A.;

Sheppard,

39,

M.

Vector

C.

3597

J.

treatment

R.; of

2004, 21, 2206.

Dholakia, K.;

M.

Eciency

of

Allen, L.;

A 1999, 60, 24382441. (28) Shinozaki,

Appl. Phys. Lett. 1991, 59, 2657

Kamijoh,

2659.

K.; T.

Xu, A

Pad-

second-harmonic

generation with Bessel beams.

for high-eciency second-harmonic gener-

C.-Q.;

comparison

Phys. Rev.

Sasaki, of

H.;

optical

second-harmonic generation eciency using Bessel and Gaussian beams in bulk

(19) Sacconi, L.; Dombeck, D. a.; Webb, W. W. Overcoming

photodamage

generation

in

crystals.

second-

microscopy:

Proc. Natl. Acad. Sci.

103, 31243129.

(20) Armstrong, J. A.;

Opt. Commun. 1997, 133, 300

304.

real-

time optical recording of neuronal action potentials.

J. Appl. Phys.

(27) Arlt, J.;

Byer, R. L. Periodically poled LiNbO3

harmonic

2016,

J. Opt. Soc. Am. B

(18) Jundt, D. H.; Magel, G. A.; Fejer, M. M.;

ation.

light

tightly focused vignetted Gaussian beams.

second and third harmonics from chirped photonic

coherent

second-harmonic generation produced by

Y. Simultaneous broadband generation of

Appl. 2014, 3, e189.

Focusing

3639.

Zhang, C.; Sheng, Y.; Ma, B.-Q.; Li, Z.-

nonlinear

Y.

Derevyanko, S.;

ric interaction of focused Gaussian light

Appl. Phys. Lett. 1977,

Ren,

7,

(25) Boyd, G. D.; Kleinman, D. A. Paramet-

(26) Asatryan,

B.-Q.;

Small, E.;

arXiv:1607.08105

monic generation from 532 to 266 nm in

30, 9193.

1990,

through a nonlinear scattering medium.

(16) Reintjes, J.; Eckardt, R. C. Ecient harADP and KD*P.

J. Opt. Soc. Am. B

776.

(29) Planchon, T. A.; Gao, L.; Milkie, D. E.;

2006,

Davidson, M. W.; Galbraith, J. A.; Galbraith, C. G.;

N., B.;

J., D.;

dimensional isotropic imaging of living

Per-

cells using Bessel beam plane illumination.

Nat. Methods 2011, 8, 417423.

shan, P. S. Interactions between Light Waves in a Nonlinear Dielectric.

Rev. 1962, 127, 19181939.

Betzig, E. Rapid three-

Phys.

ACS Paragon Plus Environment 10

Page 11 of 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

(30) Wang, H.; Shi, L.; Lukyanchuk, B.; Shep-

in second harmonic generation.

2015, 40, 19441947.

pard, C.; Chong, C. T. Creation of a needle of longitudinally polarized light in vacuum using binary optics.

2008, 2, 501505.

Nat. Photonics

(41) Vellekoop, I. M.; Mosk, A. P. Phase control algorithms for focusing light through turbid media.

(31) Abraham, N. B.; Firth, W. J. Overview of

J. Opt. Soc. Am. B

1990, 7, 951.

(42) Thompson, J.; Hokr, B.; Yakovlev, V. Optimization of focusing through scattering

(32) Magni, V. Optimum beam for second harmonic generation.

176, 245251.

Opt. Commun.

Opt. Commun. 2008, 281,

30713080.

transverse eects in nonlinear-optical systems.

Opt. Lett.

media using the continuous sequential al-

2000,

gorithm.

J. Mod. Opt. 2016, 63, 8084.

(43) Thompson, J. V.; Throckmorton, G. A.;

(33) Keren-Zur, S.; Avayu, O.; Michaeli, L.;

Hokr, B. H.; Yakovlev, V. V. Wavefront

Ellenbogen, T. Nonlinear beam shaping

shaping enhanced Raman scattering in a

with plasmonic metasurfaces.

turbid medium.

tonics 2016, 3, 117123.

ACS Pho-

Opt. Lett. 2016, 41, 1769

1772.

(34) Westerberg, N.; Maitland, C.; Faccio, D.;

(44) Vellekoop, I. M. Feedback-based wavefront

Wilson, K.; Öhberg, P.; Wright, E. M.

shaping.

Synthetic magnetism for photon uids.

12206.

Phys. Rev. A 2016, 94, 023805.

(45) Yu,

Opt. Express

H.;

Park,

J.;

23,

2015,

Lee,

K.;

12189

Yoon,

J.;

(35) Iliew, R.; Etrich, C.; Pertsch, T.; Led-

Kim, K.; Lee, S.; Park, Y. Recent ad-

erer, F.; Kivshar, Y. S. Huge enhancement

vances in wavefront shaping techniques

of backward second-harmonic generation

for biomedical applications.

with slow light in photonic crystals.

Rev. A 2010, 81, 110.

Phys.

Phys. 2015, 15, 632641.

Curr. Appl.

(46) Conkey, D. B.; Brown, A. N.; Caravaca-

(36) Smith, A. V.; Bowers, M. S. Phase dis-

Aguirre, A. M.; Piestun, R. Genetic al-

tortions in sum- and dierence-frequency

gorithm optimization for focusing through

mixing in crystals.

turbid media in noisy environments.

1995, 12, 49.

(37) Lightman, S.;

J. Opt. Soc. Am. B

Gvishi, R.;

Express 2012, 20, 48404849.

Hurvitz, G.;

(47) Katz, O.;

Small, E.;

Guan, Y.;

Opt.

Silber-

Arie, A. Shaping of light beams by 3D

berg, Y. Noninvasive nonlinear focusing

direct laser writing on facets of nonlinear

and imaging through strongly scattering

crystals.

turbid layers.

(38) Hong,

Opt. Lett. 2015, 40, 4460.

X.

H.;

Qin, Y. Q.;

Yang,

B.;

Zhang,

Optica 2014, 1, 170174.

C.;

(48) Aulbach, J.; Gjonaj, B.; Johnson, P.; La-

Zhu, Y. Y. Nonlinear vol-

gendijk, A. Spatiotemporal focusing in

ume holography for wave-front engineer-

opaque scattering media by wave front

ing.

shaping with nonlinear feedback.

Phys. Rev. Lett. 2014, 113, 15.

press 2012, 20, 2923729251.

(39) Trajtenberg-Mills, S.; Juwiler, I.; Arie, A. On-axis beams.

shaping

of

(49) Kleinman, D. A. Theory of second har-

second-harmonic

Laser Photonics Rev.

2015,

Opt. Ex-

44,

monic

generation

of

1962, 128, 17611775.

4044. (40) Libster-hershko, A.; Trajtenberg-mills, S.; Arie, A. Dynamic control of light beams

ACS Paragon Plus Environment 11

light.

Phys. Rev.

Graphical TOC Entry

Page 12 of 12

2.2

Nonlinear Crystal

Normalized Conversion Efficiency

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Nonlinear Crystal

2 1.8 1.6 1.4 1.2 1

0

200 400 600 800 Algorithm Iteration

ACS Paragon Plus Environment 12